Corrigendum: Fibril Surface-Dependent Amyloid Precursors Revealed by Coarse-Grained Molecular Dynamics Simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Ma YW
PROVIDER: S-EPMC9252589 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril-secondary nucleation-as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.
| S-EPMC8378332 | biostudies-literature
Project description:Fibrils formed by the β-amyloid (Aβ) peptide play a central role in the development of Alzheimer's disease. In this study, the principles governing their growth and stability are investigated by analyzing canonical and replica-exchange molecular dynamics trajectories of Aβ(9-40) fibrils. In particular, an unstructured monomer was allowed to interact freely with an Aβ fibril template. Trajectories were generated with the coarse-grained united-residue force field, and one- and two-dimensional free-energy landscapes (FELs) along the backbone virtual-bond angle θ and backbone virtual-bond-dihedral angle γ of each residue and principal components, respectively, were analyzed. Also, thermal unbinding (unfolding) of an Aβ peptide from the fibril template was investigated. These analyses enable us to illustrate the entire process of Aβ fibril elongation and to elucidate the key residues involved in it. Several different pathways were identified during the search for the fibril conformation by the monomer, which finally follows a dock-lock mechanism with two distinct locking stages. However, it was found that the correct binding, with native hydrogen bonds, of the free monomer to the fibril template at both stages is crucial for fibril elongation. In other words, if the monomer is incorrectly bound (with nonnative hydrogen bonds) to the fibril template during the first "docking" stage, it can remain attached to it for a long time before it dissociates and either attempts a different binding or allows another monomer to bind. This finding is consistent with an experimentally observed "stop-and-go" mechanism of fibril growth.
| S-EPMC5245180 | biostudies-literature
Project description:The mechanism of curvature generation in membranes has been studied for decades due to its important role in many cellular functions. However, it is not clear if, or how, aggregates of lipid-anchored proteins might affect the geometry and elastic property of membranes. As an initial step toward addressing this issue, we performed structural, geometrical, and stress field analyses of coarse-grained molecular dynamics trajectories of a domain-forming bilayer in which an aggregate of lipidated proteins was asymmetrically bound. The results suggest a general mechanism whereby asymmetric incorporation of lipid-modified protein aggregates curve multidomain membranes primarily by expanding the surface area of the monolayer in which the lipid anchor is inserted.
| S-EPMC3999789 | biostudies-literature
Project description:Molecular dynamics (MD) simulation has remained the most indispensable tool in studying equilibrium/non-equilibrium conformational dynamics since its advent 30 years ago. With advances in spectroscopy accompanying solved biocomplexes in growing sizes, sampling their dynamics that occur at biologically interesting spatial/temporal scales becomes computationally intractable; this motivated the use of coarse-grained (CG) approaches. CG-MD models are used to study folding and conformational transitions in reduced resolution and can employ enlarged time steps due to the absence of some of the fastest motions in the system. The Boltzmann-Inversion technique, heavily used in parameterizing these models, provides a smoothed-out effective potential on which molecular conformation evolves at a faster pace thus stretching simulations into tens of microseconds. As a result, a complete catalytic cycle of HIV-1 protease or the assembly of lipid-protein mixtures could be investigated by CG-MD to gain biological insights. In this review, we survey the theories developed in recent years, which are categorized into Folding-based and Molecular-Mechanics-based. In addition, physical bases in the selection of CG beads/time-step, the choice of effective potentials, representation of solvent, and restoration of molecular representations back to their atomic details are systematically discussed.
| S-EPMC2735960 | biostudies-literature
Project description:Many free-energy sampling and quantum mechanics/molecular mechanics (QM/MM) computations on protein complexes have been performed where, by necessity, a single component is studied isolated in solution while its overall configuration is kept in the complex-like state by either rigid restraints or harmonic constraints. A drawback in these studies is that the system's native fluctuations are lost, both due to the change of environment and the imposition of the extra potential. Yet, we know that having both accurate structure and fluctuations is likely crucial to achieving correct simulation estimates for the subsystem within its native larger protein complex context. In this work, we provide a new approach to this problem by drawing on ideas developed to incorporate experimental information into a molecular simulation by relative entropy minimization to a target system. We show that by using linear biases on coarse-grained (CG) observables (such as distances or angles between large subdomains within a protein), we can maintain the protein in a particular target conformation while also preserving the correct equilibrium fluctuations of the subsystem within its larger biomolecular complex. As an application, we demonstrate this algorithm by training a bias that causes an actin monomer (and trimer) in solution to sample the same average structure and fluctuations as if it were embedded within a much larger actin filament. Additionally, we have developed a number of algorithmic improvements that accelerate convergence of the on-the-fly relative entropy minimization algorithms for this type of application. Finally, we have contributed these methods to the PLUMED open source free energy sampling software library.
| S-EPMC5649387 | biostudies-literature
Project description:The pathological growth of amyloid fibrils in neurons underlies the progression of neurodegenerative diseases including Alzheimer's and Parkinson's disease. Fibrils form when soluble monomers oligomerise in the cytoplasm. Their subsequent growth occurs via nucleated polymerization mechanisms involving the free ends of the fibrils augmented by secondary nucleation of new oligomers at their surface. Amyloid fibrils possess a complex interactome with diffusing cytoplasmic proteins that regulates many aspects of their growth, seeding capacity, biochemical activity and transition to pathological inclusions in diseased brains. Changes to their surface are also expected to modify their interactome, pathogenicity and spreading in the brain. Many assays visualise fibril formation, growth and inclusion formation by decorating monomeric proteins with fluorescent tags such as GFP. Recent studies from our group suggest that tags with sizes comparable to the fibril radius may modify the fibril surface accessibility and thus their PTM pattern, interactome and ability to form inclusions. Using coarse-grained molecular simulations of a single alpha synuclein fibril tagged with GFP we find that thermal fluctuations of the tags create a non-monotonic, size-dependent sieve around the fibril that perturbs its interactome with diffusing species. Our results indicate that experiments using tagged and untagged monomers to study the growth and interactome of fibrils should be compared with caution, and the confounding effects of the tags are more complex than a reduction in surface accessibility. The prevalence of fluorescent tags in amyloid fibril growth experiments suggests this has implications beyond the specific alpha synuclein fibrils we model here.
| S-EPMC8753129 | biostudies-literature
Project description:Topological constraints, such as those associated with DNA supercoiling, play an integral role in genomic regulation and organization in living systems. However, physical understanding of the principles that underlie DNA organization at biologically relevant length scales remains a formidable challenge. We develop a coarse-grained simulation approach for predicting equilibrium conformations of supercoiled DNA. Our methodology enables the study of supercoiled DNA molecules at greater length scales and supercoiling densities than previously explored by simulation. With this approach, we study the conformational transitions that arise due to supercoiling across the full range of supercoiling densities that are commonly explored by living systems. Simulations of ring DNA molecules with lengths at the scale of topological domains in the Escherichia coli chromosome (∼10 kilobases) reveal large-scale conformational transitions elicited by supercoiling. The conformational transitions result in three supercoiling conformational regimes that are governed by a competition among chiral coils, extended plectonemes, and branched hyper-supercoils. These results capture the nonmonotonic relationship of size versus degree of supercoiling observed in experimental sedimentation studies of supercoiled DNA, and our results provide a physical explanation of the conformational transitions underlying this behavior. The length scales and supercoiling regimes investigated here coincide with those relevant to transcription-coupled remodeling of supercoiled topological domains, and we discuss possible implications of these findings in terms of the interplay between transcription and topology in bacterial chromosome organization.
| S-EPMC5052444 | biostudies-literature
Project description:Coarse-grained (CG) models in molecular dynamics (MD) are powerful tools to simulate the dynamics of large biomolecular systems on micro- to millisecond timescales. However, the CG model, potential energy terms, and parameters are typically not transferable between different molecules and problems. So parameterizing CG force fields, which is both tedious and time-consuming, is often necessary. We present RedMDStream, a software for developing, testing, and simulating biomolecules with CG MD models. Development includes an automatic procedure for the optimization of potential energy parameters based on metaheuristic methods. As an example we describe the parameterization of a simple CG MD model of an RNA hairpin.
| S-EPMC4407251 | biostudies-literature
Project description:Coarse-grained (CG) force fields have become promising tools for studies of protein behavior, but the balance of speed and accuracy is still a challenge in the research of protein coarse graining methodology. In this work, 20 CG beads have been designed based on the structures of amino acid residues, with which an amino acid can be represented by one or two beads, and a CG solvent model with five water molecules was adopted to ensure the consistence with the protein CG beads. The internal interactions in protein were classified according to the types of the interacting CG beads, and adequate potential functions were chosen and systematically parameterized to fit the energy distributions. The proposed CG force field has been tested on eight proteins, and each protein was simulated for 1000 ns. Even without any extra structure knowledge of the simulated proteins, the Cα root mean square deviations (RMSDs) with respect to their experimental structures are close to those of relatively short time all atom molecular dynamics simulations. However, our coarse grained force field will require further refinement to improve agreement with and persistence of native-like structures. In addition, the root mean square fluctuations (RMSFs) relative to the average structures derived from the simulations show that the conformational fluctuations of the proteins can be sampled.
| S-EPMC3509591 | biostudies-literature
Project description:A molecular dynamics investigation and coarse-grained analysis of inactivated actin-related protein (Arp) 2/3 complex is presented. It was found that the nucleotide binding site within Arp3 remained in a closed position with bound ATP or ADP, but opened when simulation with no nucleotide was performed. In contrast, simulation of the isolated Arp3 subunit with bound ATP, showed a fast opening of the nucleotide binding cleft. A homology model for the missing subdomains 1 and 2 of Arp2 was constructed, and it was also found that the Arp2 binding cleft remained closed with bound nucleotide. Within the nucleotide binding cleft a distinct opening and closing period of 10 ns was observed in many of the simulations of Arp2/3 as well as isolated Arp3. Substitution studies were employed, and several alanine substitutions were found to induce a partial opening of the ATP binding cleft in Arp3 and Arp2, whereas only a single substitution was found to induce opening of the ADP binding cleft. It was also found that the nucleotide type did not cause a substantial change on interfacial contacts between Arp3 and the ArpC2, ArpC3 and ArpC4 subunits. Nucleotide-free Arp3 had generally less stable contacts, but the overall contact architecture was constant. Finally, nucleotide-dependent coarse-grained models for Arp3 are developed that serve to further highlight the structural differences induced in Arp3 by nucleotide hydrolysis.
| S-EPMC2586551 | biostudies-literature