Amyloid β Fibril Elongation by Monomers Involves Disorder at the Tip

Molecular Insights into the Dynamics of Amyloid Fibril Growth: Elongation and Lateral Assembly of GNNQQNY Protofibrils

The self-assembly of peptides and proteins into β-sheet rich amyloid fibrils is linked to both functional and pathological states. In this study, the growth of fibrillar structures of the short peptide GNNQQNY, a fragment from the yeast prion Sup35 protein, was examined. Molecular dynamics simulations were used to study alternative mechanisms of fibril growth, including elongation through binding of monomers as well as fibril self-assembly into larger, more mature structures. It was found that after binding, monomers diffused along preformed fibrils toward the ends, supporting the mechanism of fibril growth via elongation. Lateral assembly of protofibrils was found to occur readily, suggesting that this could be the key to transitioning from isolated fibrils to mature multilayer structures. Overall, the work provides mechanistic insights into the competitive pathways that govern amyloid fibril growth.

Interactions between Lipid Vesicle Membranes and Single Amino Acid Fibrils: Probable Origin of Specific Neurological Disorders

Amyloid fibrils are known to be responsible for several neurological disorders, like Alzheimer's disease (AD), Parkinson's disease (PD), etc. For decades, mostly proteins and peptide-based amyloid fibrils have been focused on, and the topic has acknowledged the rise, development, understanding of, and controversy, as well. However, the single amino acid based amyloid fibrils, responsible for several disorders, such as phenylketonuria, tyrosenimia type II, hypermethioninemia, etc., have gotten scientific attention lately. To understand the molecular level pathogenesis of such disorders originated due to the accumulation of single amino acid-based amyloid fibrils, interaction of these fibrils with phospholipid vesicle membranes is found to be an excellent cell-free in vitro setup. Based on such an in vitro setup, these fibrils show a generic mechanism of membrane insertion driven by electrostatic and hydrophobic effects inside the membrane that reduces the integral rigidity of the membrane. Alteration of such fundamental properties of the membrane, therefore, might be referred to as one of the prime pathological factors for the development of these neurological disorders. Hence, such interactions must be investigated in cellular and intracellular compartments to design suitable therapeutic modulators against fibrils.

Protein Corona Formation and Aggregation of Amyloid β 1-40-Coated Gold Nanocolloids

Amyloid fibrillogenesis is a pathogenic protein aggregation process that occurs through a highly ordered process of protein-protein interactions. To better understand the protein-protein interactions involved in amyloid fibril formation, we formed nanogold colloid aggregates by stepwise additions of ∼2 nmol of amyloid β 1-40 peptide (Aβ1-40) at pH ∼3.7 and ∼25 °C. The processes of protein corona formation and building of gold colloid [diameters (d) of 20 and 80 nm] aggregates were confirmed by a red-shift of the surface plasmon resonance (SPR) band, λpeak, as the number of Aβ1-40 peptides [N(Aβ1-40)] increased. The normalized red-shift of λpeak, Δλ, was correlated with the degree of protein aggregation, and this process was approximated as the adsorption isotherm explained by the Langmuir-Freundlich model. As the coverage fraction (θ) was analyzed as a function of ϕ, which is the N(Aβ1-40) per total surface area of nanogold colloids available for adsorption, the parameters for explaining the Langmuir-Freundlich model were in good agreement for both 20 and 80 nm gold, indicating that ϕ could define the stage of the aggregation process. Surface-enhanced Raman scattering (SERS) imaging was conducted at designated values of ϕ and suggested that a protein-gold surface interaction during the initial adsorption stage may be dependent on the nanosize. The 20 nm gold case seems to prefer a relatively smaller contacting section, such as a -C-N or C═C bond, but a plane of the benzene ring may play a significant role for 80 nm gold. Regardless of the size of the particles, the β-sheet and random coil conformations were considered to be used to form gold colloid aggregates. The methodology developed in this study allows for new insights into protein-protein interactions at distinct stages of aggregation.

Key Residue for Aggregation of Amyloid-β Peptides

It is known that oligomers of amyloid-β (Aβ) peptide are associated with Alzheimer's disease. Aβ has two isoforms: Aβ40 and Aβ42. Although the difference between Aβ40 and Aβ42 is only two additional C-terminal residues, Aβ42 aggregates much faster than Aβ40. It is unknown what role the C-terminal two residues play in accelerating aggregation. Since Aβ42 is more toxic than Aβ40, its oligomerization process needs to be clarified. Moreover, clarifying the differences between the oligomerization processes of Aβ40 and Aβ42 is essential to elucidate the key factors of oligomerization. Therefore, to investigate the dimerization process, which is the early oligomerization process, Hamiltonian replica-permutation molecular dynamics simulations were performed for Aβ40 and Aβ42. We identified a key residue, Arg5, for the Aβ42 dimerization. The two additional residues in Aβ42 allow the C-terminus to form contact with Arg5 because of the electrostatic attraction between them, and this contact stabilizes the β-hairpin. This β-hairpin promotes dimer formation through the intermolecular β-bridges. Thus, we examined the effects of amino acid substitutions of Arg5, thereby confirming that the mutations remarkably suppressed the aggregation of Aβ42. Moreover, the mutations of Arg5 suppressed the Aβ40 aggregation. It was found by analyzing the simulations that Arg5 is important for Aβ40 to form intermolecular contacts. Thus, it was clarified that the role of Arg5 in the oligomerization process varies due to the two additional C-terminal residues.

Antibody binding modulates the dynamics of the membrane-bound prion protein

Misfolding of the cellular prion protein (PrP^C) is associated with lethal neurodegeneration. PrP^C consists of a flexible tail (residues 23-123) and a globular domain (residues 124-231) whose C-terminal end is anchored to the cell membrane. The neurotoxic antibody POM1 and the innocuous antibody POM6 recognize the globular domain. Experimental evidence indicates that POM1 binding to PrP^C emulates the influence on PrP^C of the misfolded prion protein (PrP^Sc) while the binding of POM6 has the opposite biological response. Little is known about the potential interactions between flexible tail, globular domain, and the membrane. Here, we used atomistic simulations to investigate how these interactions are modulated by the binding of the Fab fragments of POM1 and POM6 to PrP^C and by interstitial sequence truncations to the flexible tail. The simulations show that the binding of the antibodies restricts the range of orientations of the globular domain with respect to the membrane and decreases the distance between tail and membrane. Five of the six sequence truncations influence only marginally this distance and the contact patterns between tail and globular domain. The only exception is a truncation coupled to a charge inversion mutation of four N-terminal residues, which increases the distance of the flexible tail from the membrane. The interactions of the flexible tail and globular domain are modulated differently by the two antibodies.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Single-molecule Observation of Self-Propagating Amyloid Fibrils

The assembly of misfolded proteins into amyloid fibrils is associated with amyloidosis, including neurodegenerative diseases, such as Alzheimer’s, Parkinson’s, and prion diseases. The self-propagation of amyloid fibrils is widely observed in the aggregation pathways of numerous amyloidogenic proteins. This propensity with plasticity in primary nucleation allows amyloid fibril polymorphism, which is correlated with the pathology/phenotypes of patients. Because the interference with the nucleation and replication processes of amyloid fibrils can alter the amyloid structure and the outcome of the disease, these processes can be a target for developing clinical drugs. Single-molecule observation of amyloid fibril replication can be an experimental system to provide the kinetic parameters for simulation studies and confirm the effect of clinical drugs. Here, we review single-molecule observation of the amyloid fibril replication process using fluorescence microscopy and time-lapse atomic force microscopy, including high-speed atomic force microscopy. We discussed the amyloid fibril replication process and combined single-molecule observation results with molecular dynamics simulations.

Mini Abstract Structural dynamics in amyloid aggregation is related with various Alzheimer’s and Parkinson’s disease symptoms. Single-molecule observation using high-speed atomic force microscopy can directly visualize the structural dynamics of individual amyloid aggregate assemblies. Here, we review historical and recent studies of single-molecule observation of amyloid aggregation with supportive molecular dynamics simulation.

Insights into Mechanism of Aβ42 Fibril Growth on Surface of Graphene Oxides: Oxidative Degree Matters

The filamentous β-amyloid deposition has been regarded as the hallmark pathology of Alzheimer's disease (AD). Nanomaterials such as graphene oxides (GOs) have achieved significant progress in the therapy of AD, but the molecular pathway of the growth propagation remains challenging to investigate, especially on the surfaces of materials. The thermodynamics and kinetics of fibril elongation on GO surfaces with different oxidative degrees have been investigated by a combination of in vitro experiments and simulations. ThT kinetics, calorimetric measurements, and TEM observations suggest that low oxidative GO-10 promotes the fibril elongation, while both high oxidative GO-20 and GO-40 inhibit the fibril elongation. Computational results reveal that the apparent regulation behaviors of GOs on filament growth depend on the balance between the promoting effect by templating the incoming of monomers and the retarding effect by capturing the monomer during docking and locking phases through hydrogen bonding. This work will promote the understanding of the interplay between biomolecules and materials, thus providing new thoughts for the rational design of novel materials for amyloidosis therapy.

Molecular dynamics study of water channels in natural and synthetic amyloid-β fibrils

We compared all-atom explicit solvent molecular dynamics simulations of three types of Aβ(1-40) fibrils: brain-seeded fibrils (2M4J, with a threefold axial symmetry) and the other two, all-synthetic fibril polymorphs (2LMN and 2LMP, made under different fibrillization conditions). Fibril models were constructed using either a finite or an infinite number of layers made using periodic images. These studies yielded four conclusions. First, finite fibrils tend to unravel in a manner reminiscent of fibril dissolution, while infinite fibrils were more stable during simulations. Second, salt bridges in these fibrils remained stable in those fibrils that contained them initially, and those without salt bridges did not develop them over the time course of the simulations. Third, all fibrils tended to develop a "stagger" or register shift of β-strands along the fibril axis. Fourth and most importantly, the brain-seeded, 2M4J, infinite fibrils allowed bidirectional transport of water in and out of the central longitudinal core of the fibril by rapidly developing gaps at the fibril vertices. 2LMP fibrils also showed this behavior, although to a lesser extent. The diffusion of water molecules in the fibril core region involved two dynamical states: a localized state and directed diffusion in the presence of obstacles. These observations provided support for the hypothesis that Aβ fibrils could act as nanotubes. At least some Aβ oligomers resembled fibrils structurally in having parallel, in-register β-sheets and a sheet-turn-sheet motif. Thus, our findings could have implications for Aβ cytotoxicity, which may occur through the ability of oligomers to form abnormal water and ion channels in cell membranes.

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

Redox-Dependent Copper Ion Modulation of Amyloid-β (1-42) Aggregation In Vitro

Plaque deposits composed of amyloid-β (Aβ) fibrils are pathological hallmarks of Alzheimer's disease (AD). Although copper ion dyshomeostasis is apparent in AD brains and copper ions are found co-deposited with Aβ peptides in patients' plaques, the molecular effects of copper ion interactions and redox-state dependence on Aβ aggregation remain elusive. By combining biophysical and theoretical approaches, we here show that Cu2+ (oxidized) and Cu+ (reduced) ions have opposite effects on the assembly kinetics of recombinant Aβ(1-42) into amyloid fibrils in vitro. Cu2+ inhibits both the unseeded and seeded aggregation of Aβ(1-42) at pH 8.0. Using mathematical models to fit the kinetic data, we find that Cu2+ prevents fibril elongation. The Cu2+-mediated inhibition of Aβ aggregation shows the largest effect around pH 6.0 but is lost at pH 5.0, which corresponds to the pH in lysosomes. In contrast to Cu2+, Cu+ ion binding mildly catalyzes the Aβ(1-42) aggregation via a mechanism that accelerates primary nucleation, possibly via the formation of Cu+-bridged Aβ(1-42) dimers. Taken together, our study emphasizes redox-dependent copper ion effects on Aβ(1-42) aggregation and thereby provides further knowledge of putative copper-dependent mechanisms resulting in AD.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Steady, Symmetric, and Reversible Growth and Dissolution of Individual Amyloid-β Fibrils

Oligomers and fibrils of the amyloid-β (Aβ) peptide are implicated in the pathology of Alzheimer's disease. Here, we monitor the growth of individual Aβ40 fibrils by time-resolved in situ atomic force microscopy and thereby directly measure fibril growth rates. The measured growth rates in a population of fibrils that includes both single protofilaments and bundles of filaments are independent of the fibril thickness, indicating that cooperation between adjacent protofilaments does not affect incorporation of monomers. The opposite ends of individual fibrils grow at similar rates. In contrast to the "stop-and-go" kinetics that has previously been observed for amyloid-forming peptides, growth and dissolution of the Aβ40 fibrils are relatively steady for peptide concentration of 0-10 μM. The fibrils readily dissolve in quiescent peptide-free solutions at a rate that is consistent with the microscopic reversibility of growth and dissolution. Importantly, the bimolecular rate coefficient for the association of a monomer to the fibril end is significantly smaller than the diffusion limit, implying that the transition state for incorporation of a monomer into a fibril is associated with a relatively high free energy.

Network-Based Classification and Modeling of Amyloid Fibrils

Amyloid fibrils are locally ordered protein aggregates that self-assemble under a variety of physiological and in vitro conditions. Their formation is of fundamental interest as a physical chemistry problem and plays a central role in Alzheimer's disease, Type II diabetes, and other human diseases. As the number of known amyloid fibril structures has grown, the need has arisen for a nomenclature for describing and classifying fibril types, as well as a theoretical description of the physics that gives rise to the self-assembly of these structures. Here, we introduce a systematic nomenclature and coarse-graining methodology for describing the topology of fibrils and other protein aggregates, along with a computational methodology for simulating protein aggregation. Both have mathematical underpinnings in graph theory and statistical mechanics and are consistent with available experimental data on the fibril structure and aggregation kinetics. Our graph representation of the fibril topology enables us to define a network Hamiltonian based on connectivity patterns among monomers rather than detailed intermolecular interactions, greatly speeding up the simulation of large ensembles. Our simulation strategy is capable of recapitulating the formation of all currently known amyloid fibril topologies found in the Protein Data Bank, as well as the formation kinetics of fibrils and oligomers.

On the removal of initial state bias from simulation data

Classical atomistic simulations of biomolecules play an increasingly important role in molecular life science. The structure of current computing architectures favors methods that run multiple trajectories at once without requiring extensive communication between them. Many advanced sampling strategies in the field fit this mold. These approaches often rely on an adaptive logic and create ensembles of comparatively short trajectories whose starting points are not distributed according to the correct Boltzmann weights. This type of bias is notoriously difficult to remove, and Markov state models (MSMs) are one of the few strategies available for recovering the correct kinetics and thermodynamics from these ensembles of trajectories. In this contribution, we analyze the performance of MSMs in the thermodynamic reweighting task for a hierarchical set of systems. We show that MSMs can be rigorous tools to recover the correct equilibrium distribution for systems of sufficiently low dimensionality. This is conditional upon not tampering with local flux imbalances found in the data. For a real-world application, we find that a pure likelihood-based inference of the transition matrix produces the best results. The removal of the bias is incomplete, however, and for this system, all tested MSMs are outperformed by an alternative albeit less general approach rooted in the ideas of statistical resampling. We conclude by formulating some recommendations for how to address the reweighting issue in practice.

Disorder at the Tips of a Disease-Relevant Aβ42 Amyloid Fibril: A Molecular Dynamics Study

We present a simulation study of the early events of peptide dissociation from a fibril of the Alzheimer's Aβ42 peptide. The fibril consists of layers of two adjacent Aβ42 peptides each folded in an S-shaped structure which has been determined by solid state NMR spectroscopy of a monomorphic disease-relevant species. Multiple molecular dynamics runs (16 at 310 K and 15 at 370 K) were carried out starting from an 18-peptide protofibril for a cumulative sampling of about 15 μs. The simulations show structural stability of the fibrillar core and an overall increase in the twist to about 3 degrees. The N-terminal segment 1-14 is disordered in all peptides. At both ends of the fibril, the central segment 21-29, which includes part of the β₂ strand, dissociates in some of the simulations. The β₁ and β₃ strands, residues 15-20 and 35-41, respectively, are structurally stable. The transient binding of the N-terminal stretch to the β₃ strand of the adjacent peptide at the tip is likely to contribute to the arrest phase of the stop-and-go mechanism.

Focused conformational sampling in proteins

A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.

Thermodynamics of Amyloid-β Fibril Elongation: Atomistic Details of the Transition State

Amyloid-β (Aβ) fibrils and plaques are one of the hallmarks of Alzheimer’s disease. While the kinetics of fibrillar growth of Aβ have been extensively studied, several vital questions remain. In particular, the atomistic origins of the Arrhenius barrier observed in experiments have not been elucidated. Employing the familiar thermodynamic integration method, we have directly simulated the dissociation of an Aβ_(15–40) (D23N mutant) peptide from the surface of a filament along its most probable path (MPP) using all-atom molecular dynamics. This allows for a direct calculation of the free energy profile along the MPP, revealing a multipeak energetic barrier between the free peptide state and the aggregated state. By definition of the MPP, this simulated unbinding process represents the reverse of the physical elongation pathway, allowing us to draw biophysically relevant conclusions from the simulation data. Analyzing the detailed atomistic interactions along the MPP, we identify the atomistic origins of these peaks as resulting from the dock-lock mechanism of filament elongation. Careful analysis of the dynamics of filament elongation could prove key to the development of novel therapeutic strategies for amyloid-related diseases.

Energy Landscapes for the Aggregation of Aβ17-42

The aggregation of the Aβ peptide (Aβ_1-42) to form fibrils is a key feature of Alzheimer's disease. The mechanism is thought to be a nucleation stage followed by an elongation process. The elongation stage involves the consecutive addition of monomers to one end of the growing fibril. The aggregation process proceeds in a stop-and-go fashion and may involve off-pathway aggregates, complicating experimental and computational studies. Here we present exploration of a well-defined region in the free and potential energy landscapes for the Aβ_17-42 pentamer. We find that the ideal aggregation process agrees with the previously reported dock-lock mechanism. We also analyze a large number of additional stable structures located on the multifunnel energy landscape, which constitute kinetic traps. The key contributors to the formation of such traps are misaligned strong interactions, for example the stacking of F19 and F20, as well as entropic contributions. Our results suggest that folding templates for aggregation are a necessity and that aggregation studies could employ such species to obtain a more detailed description of the process.

Self-Assembly Pathways of β-Sheet-Rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations

Early oligomerization during amyloid-β (Aβ) aggregation is essential for Aβ neurotoxicity. Understanding how unstructured Aβs assemble into oligomers, especially those rich in β-sheets, is essential but remains challenging as the assembly process is too transient for experimental characterization and too slow for molecular dynamics simulations. So far, atomic simulations are limited only to studies of either oligomer structures or assembly pathways for short Aβ segments. To overcome the computational challenge, we combine in this study a hybrid-resolution model and adaptive sampling techniques to perform over 2.7 ms of simulations of formation of full-length Aβ40 dimers that are the earliest toxic oligomeric species. The Markov state model is further employed to characterize the transition pathways and associated kinetics. Our results show that for two major forms of β-sheet-rich structures reported experimentally, the corresponding assembly mechanisms are markedly different. Hairpin-containing structures are formed by direct binding of soluble Aβ in β-hairpin-like conformations. Formation of parallel, in-register structures resembling fibrils occurs ∼100-fold more slowly and involves a rapid encounter of Aβ in arbitrary conformations followed by a slow structural conversion. The structural conversion proceeds via diverse pathways but always requires transient unfolding of encounter complexes. We find that the transition kinetics could be affected differently by intra-/intermolecular interactions involving individual residues in a conformation-dependent manner. In particular, the interactions involving Aβ's N-terminal part promote the assembly into hairpin-containing structures but delay the formation of fibril-like structures, thus explaining puzzling observations reported previously regarding the roles of this region in the early assembly process.