Amyloid oligomer structure characterization from simulations: A general method

Molecular Insights into Cu/Zn Metal Response to the Amyloid β-Peptide (1-42)

Aβ1-40 peptide and Aβ1-42 peptide are the building units of beta-amyloid plaques present in Alzheimer's disease (AD)-affected brain. The binding affinity of various divalent metal ions such as Cu and Zn present in AD-affected brain with different amino acids available in Aβ-peptide became the focus to explore their role in soluble neurotoxic oligomer formation. Cu2+ metal ions are known to enhance the neurotoxicity of the Aβ1-42 peptide by catalyzing the formation of soluble neurotoxic oligomers. The competitive preference of both Cu2+ and Zn2+ simultaneously to interact with the Aβ-peptide is unknown. The divalent Cu and Zn ions were inserted in explicit aqueous Aβ1-42 peptide configurations to get insights into the binding competence of these metal ions with peptides using classical molecular dynamics (MD) simulations. The metal-ion interactions reveal that competitive binding preferences of various peptide sites become metal-ion-specific and differ significantly. For Cu2+, interactions are found to be more significant with respect to those of Asp-7, His-6, Glu-11, and His-14. Asp-1, Glu-3, Asp-7, His-6, Glu-11, and His-13 amino acid residues show higher affinity toward Zn2+ ions. MD simulations show notable variation in the solvent-accessible surface area in the hydrophobic region of the peptide. Infinitesimal mobility was obtained for Zn2+ compared to Cu2+ in an aqueous solution and Cu2+ diffusivity deviated significantly at different time scales, proving its labile features in aqueous Aβ1-42 peptides.

Insights into the Mixture of Aβ24 and Aβ42 Peptides from Atomistic Simulations

Amyloid-β (Aβ) oligomers play a central role in Alzheimer's disease (AD). Plaques of AD patients consist of Aβ40 and Aβ42 peptides and truncated Aβ peptides. The Aβ24 peptide, identified in human AD brains, was found to impair Aβ42 clearance through the brain-blood barrier. The Aβ24 peptide was also shown to reduce Aβ42 aggregation kinetics in pure buffer, but the underlying mechanism is unknown at atomistic level. In this study, we explored the conformational ensemble of the equimolar mixture of Aβ24 and Aβ42 by replica exchange molecular dynamics simulations and compared it to our previous results on the pure Aβ42 dimer. Our simulations demonstrate that the truncation at residue 24 changes the secondary, tertiary, and quaternary structures of the dimer, offering an explanation of the slower aggregation kinetics of the mixture.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Molecular Mechanism and Kinetics of Amyloid-β42 Aggregate Formation: A Simulation Study

As an important neuropathological hallmark of Alzheimer's disease (AD), the oligomerization of amyloid-β (Aβ) peptides has been intensively investigated in both theoretical and experimental studies. However, the oligomerization space in terms of the kinetics, molecular mechanism, and oligomer structures remains mysterious to us. An equation that can quantitatively describe the time it takes for Aβ oligomers to appear in the human brain at a given Aβ monomer concentration is extremely vital for us to understand the development and disease progression of AD. In this study, we utilized molecular dynamics (MD) simulations to investigate the oligomerization of Aβ₄₂ peptides at five different monomer concentrations. We have elucidated the formation pathways of Aβ tetramers, characterized the oligomer structures, estimated the oligomerization time for Aβ dimers, trimers, and tetramers, and for the first-time derived equations that could quantitatively describe the relationship between the oligomerization time and the monomer concentration. Applying these equations, our prediction of oligomerization time agrees well with the experimental and clinical findings, in spite of the limitations of our oligomerization simulations. We have found that the Aβ oligomerization time depends on the monomer concentration by a power of -2.4. The newly established equations will enable us to quantitatively estimate the risk score of AD, which is a function of age. Moreover, we have identified the most dominant pathway of forming Aβ tetramers, probably the most important and toxic Aβ oligomer. Our results have shown that the structures of Aβ₄₂ dimer, trimer, and tetramer, which are distinguishable from each other, depend on the monomer concentration at which the oligomers form. Representative oligomer structures, which can serve as potential drug targets, have been identified by clustering analysis. The MD sampling adequacy has been validated by the excellent agreement between the calculated and measured collisional cross section (CCS) parameters (the prediction errors are within 2%). In a conclusion, this study provides the kinetics and structure basis for developing inhibitors to decelerate the Aβ oligomerization process.

The Precursor to Glutathione (GSH), γ-Glutamylcysteine (GGC), Can Ameliorate Oxidative Damage and Neuroinflammation Induced by Aβ40 Oligomers in Human Astrocytes

Glutathione (GSH) is one of the most abundant thiol antioxidants in cells. Many chronic and age-related diseases are associated with a decline in cellular GSH levels or impairment in the catalytic activity of the GSH biosynthetic enzyme glutamate cysteine ligase (GCL). γ-glutamylcysteine (GGC), a precursor to glutathione (GSH), can replenish depleted GSH levels under oxidative stress conditions, by circumventing the regulation of GSH biosynthesis and providing the limiting substrate. Soluble amyloid-β (Aβ) oligomers have been shown to induce oxidative stress, synaptic dysfunction and memory deficits which have been reported in Alzheimer’s disease (AD). Calcium ions, which are increased with age and in AD, have been previously reported to enhance the formation of Aβ₄₀ oligomers, which have been casually associated with the pathogenesis of the underlying neurodegenerative condition. In this study, we examined the potential beneficial effects of GGC against exogenous Aβ₄₀ oligomers on biomarkers of apoptosis and cell death, oxidative stress, and neuroinflammation, in human astrocytes. Treatment with Aβ₄₀ oligomers significantly reduced the cell viability and apoptosis of astrocyte brain cultures and increased oxidative modifications of DNA, lipids, and protein, enhanced pro-inflammatory cytokine release and increased the activity of the proteolytic matrix metalloproteinase enzyme, matric metalloproteinase (MMP)-2 and reduced the activity of MMP-9 after 24 h. Co-treatment of Aβ₄₀ oligomers with GGC at 200 μM increased the activity of the antioxidant enzymes superoxide dismutase (SOD) and glutathione peroxidase (GPx) and led to significant increases in the levels of the total antioxidant capacity (TAC) and GSH and reduced the GSSG/GSH ratio. GGC also upregulated the level of the anti-inflammatory cytokine IL-10 and reduced the levels of the pro-inflammatory cytokines (TNF-α, IL-6, and IL-1β) and attenuated the changes in metalloproteinase activity in oligomeric Aβ₄₀-treated astrocytes. Our data provides renewed insight on the beneficial effects of increased GSH levels by GGC in human astrocytes, and identifies yet another potential therapeutic strategy to attenuate the cytotoxic effects of Aβ oligomers in AD.

Understanding the Twisted Structure of Amyloid Fibrils via Molecular Simulations

Accumulation and aggregation of amyloid are associated with the pathogenesis of many human diseases, such as Alzheimer's disease and Type 2 diabetes mellitus. Therefore, a quantitative understanding of the molecular mechanisms causing different aggregated structures and biomechanical properties of amyloid fibrils could shed some light into the progression of these diseases. In this work, we develop coarse-grained molecular dynamics (CGMD) models to simulate the dynamic self-assembly of two types of amyloids (amylin and amyloid β (Aβ)). We investigate the structural and mechanical properties of different types of aggregated amyloid fibrils. Our simulations demonstrate that amyloid fibrils could result from longitudinal growth of protofilament bundles, confirming one of the hypotheses on the fibril formation. In addition, we find that the persistence length of amylin fibrils increases concurrently with their pitch length, suggesting that the bending stiffness of amylin fibrils becomes larger when the amylin fibrils are less twisted. Similar results are observed for Aβ fibrils. These findings quantify the connection between the structural and the biomechanical properties of the fibrils. The CGMD models developed in this work can be potentially used to examine efficacy of anti-aggregation drugs, which could help in developing new treatments.

Glycation induces conformational changes in the amyloid-β peptide and enhances its aggregation propensity: molecular insights

The cytotoxicity of the amyloid beta (Aβ) peptide, implicated in the pathogenesis of Alzheimer's disease (AD), can be enhanced by its post-translational glycation, a series of non-enzymatic reactions with reducing sugars and reactive dicarbonyls. However, little is known about the underlying mechanisms that potentially enhance the cytotoxicity of the advanced glycation modified Aβ. In this work, fully atomistic molecular dynamics (MD) simulations are exploited to obtain direct molecular insights into the process of early Aβ self-assembly in the presence and absence of glycated lysine residues. Analyses of data exceeding cumulative timescales of 1 microsecond for each system reveal that glycation results in a stronger enthalpy of association between Aβ monomers and lower conformational entropy, in addition to a sharp overall increase in the beta-sheet content. Further analyses reveal that the enhanced interactions originate, in large part, due to markedly stronger, as well as new, inter-monomer salt bridging propensities in the glycated variety. Interestingly, these conformational and energetic effects are broadly reflected in preformed protofibrillar forms of Aβ small oligomers modified with glycation. Our combined results imply that glycation consolidates Aβ self-assembly regardless of its point of occurrence in the pathway. They provide a basis for further mechanistic studies and therapeutic endeavors that could potentially result in novel ways of combating AGE related AD progression.

High-Resolution Structures of the Amyloid-β 1-42 Dimers from the Comparison of Four Atomistic Force Fields

The dimer of the amyloid-β peptide Aβ of 42 residues is the smallest toxic species in Alzheimer's disease, but its equilibrium structures are unknown. Here we determined the equilibrium ensembles generated by the four atomistic OPLS-AA, CHARMM22*, AMBER99sb-ildn, and AMBERsb14 force fields with the TIP3P water model. On the basis of 144 μs replica exchange molecular dynamics simulations (with 750 ns per replica), we find that the four force fields lead to random coil ensembles with calculated cross-collision sections, hydrodynamics properties, and small-angle X-ray scattering profiles independent of the force field. There are, however, marked differences in secondary structure, with the AMBERsb14 and CHARMM22* ensembles overestimating the CD-derived helix content, and the OPLS-AA and AMBER99sb-ildn secondary structure contents in agreement with CD data. Also the intramolecular beta-hairpin content spanning residues 17-21 and 30-36 varies between 1.5% and 13%. Overall, there are significant differences in tertiary and quaternary conformations among all force fields, and the key finding, irrespective of the force field, is that the dimer is stabilized by nonspecific interactions, explaining therefore its possible transient binding to multiple cellular partners and, in part, its toxicity.

Conformational Ensembles of the Wild-Type and S8C Aβ1-42 Dimers

We characterized the dimer of the amyloid-β wild-type (WT) peptide, Aβ, of 42 residues and its disulfide-bond-locked double mutant (S8C) by replica exchange molecular dynamics simulations. Aβ dimers are known to be the smallest toxic species in Alzheimer's disease, and the S8C mutant has been shown experimentally to form an exclusive homogeneous and neurotoxic dimer. Our 50 μs all-atom simulations reveal similar secondary structures and collision cross-sections but very different intramolecular and intermolecular conformations upon double S8C mutation. Both dimers are very dynamic with hundreds of free-energy minima that differ from the U-shape and S-shape conformations of the peptides in the fibrils. The only common structural feature, shared by both species with a probability of 4% in WT and 12% in S8C-S8C, is a three-stranded β-sheet spanning the 17-23, 29-36, and 39-41 residues, which does not exist in the Aβ40 WT dimers.

Replica exchange molecular dynamics study of the truncated amyloid beta (11–40) trimer in solution

The structure of the 3Aβ_11–40 oligomer is determined for the first time using T-REMD simulations.

Dimerization Mechanism of Alzheimer Aβ40 Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation

Amyloid beta (Aβ) oligomerization is associated with the origin and progression of Alzheimer's disease (AD). While the A2V mutation enhances aggregation kinetics and toxicity, mixtures of wild-type (WT) and A2V, and also WT and A2T, peptides retard fibril formation and protect against AD. In this study, we simulate the equilibrium ensemble of WT:A2T Aβ₄₀ dimer by means of extensive atomistic replica exchange molecular dynamics and compare our results with previous equivalent simulations of A2V:A2V, WT:WT, and WT:A2V Aβ₄₀ dimers for a total time scale of nearly 0.1 ms. Qualitative comparison of the resulting thermodynamic properties, such as the relative binding free energies, with the reported experimental kinetic and thermodynamic data affords us important insight into the conversion from slow-pathway to fast-pathway dimer conformations. The crucial reaction coordinate or driving force of such transformation turns out to be related to hydrophobic interpeptide interactions. Analysis of the equilibrium ensembles shows that the fast-pathway conformations contain interpeptide out-of-register antiparallel β-sheet structures at short interpeptide distances. In contrast, the slow-pathway conformations are formed by the association of peptides at large interpeptide distances and high intrapeptide compactness, such as conformations containing intramolecular three-stranded β-sheets which sharply distinguish fast (A2V:A2V and WT:WT) and slow (WT:A2T and WT:A2V) amyloid-forming sequences. Also, this analysis leads us to predict that a molecule stabilizing the intramolecular three-stranded β-sheet or inhibiting the formation of an interpeptide β-sheet spanning residues 17-20 and 31-37 would further reduce fibril formation and probably the cytotoxicity of Aβ species.

Impact of the A2V Mutation on the Heterozygous and Homozygous Aβ1-40 Dimer Structures from Atomistic Simulations

The A2V mutation was reported to protect from Alzheimer's disease in its heterozygous form and cause an early Alzheimer's disease type dementia in its homozygous form. Experiments showed that the aggregation rate follows the order A2V > WT (wild-type) > A2V-WT. To understand the impact of this mutation, we carried out replica exchange molecular dynamics simulations of Aβ1-40 WT-A2V and A2V-A2V dimers and compared to the WT dimer. Our atomistic simulations reveal that the mean secondary structure remains constant, but there are substantial differences in the intramolecular and intermolecular conformations upon single and double A2V mutation. Upon single mutation, the intrinsic disorder is reduced, the intermolecular potential energies are reduced, the population of intramolecular three-stranded β-sheets is increased, and the number of all α dimer topologies is decreased. Taken together, these results offer an explanation for the reduced aggregation rate of the Aβ1-40 A2V-WT peptides and the protective effect of A2V in heterozygotes.

Structures of the Alzheimer's Wild-Type Aβ1-40 Dimer from Atomistic Simulations

We have studied the dimer of amyloid beta peptide Aβ of 40 residues by means of all-atom replica exchange molecular dynamics. The Aβ-dimers have been found to be the smallest toxic species in Alzheimer's disease, but their inherent flexibilities have precluded structural characterization by experimental methods. Though the 24-μs-scale simulation reveals a mean secondary structure of 18% β-strand and 10% α helix, we find transient configurations with an unstructured N-terminus and multiple β-hairpins spanning residues 17-21 and 30-36, but the antiparallel and perpendicular peptide orientations are preferred over the parallel organization. Short-lived conformational states also consist of all α topologies, and one compact peptide with β-sheet structure stabilized by a rather extended peptide with α-helical content. Overall, this first all-atom study provides insights into the equilibrium structure of the Aβ1-40 dimer in aqueous solution, opening a new avenue for a comprehensive understanding of the impact of pathogenic and protective mutations in early-stage Alzheimer's disease on a molecular level.

Combined experimental and simulation studies suggest a revised mode of action of the anti-Alzheimer disease drug NQ-Trp

Inhibition of the aggregation of the monomeric peptide β-amyloid (Aβ) into oligomers is a widely studied therapeutic approach in Alzheimer's disease (AD). Many small molecules have been reported to work in this way, including 1,4-naphthoquinon-2-yl-L-tryptophan (NQ-Trp). NQ-Trp has been reported to inhibit aggregation, to rescue cells from Aβ toxicity, and showed complete phenotypic recovery in an in vivo AD model. In this work we investigated its molecular mechanism by using a combined approach of experimental and theoretical studies, and obtained converging results. NQ-Trp is a relatively weak inhibitor and the fluorescence data obtained by employing the fluorophore widely used to monitor aggregation into fibrils can be misinterpreted due to the inner filter effect. Simulations and NMR experiments showed that NQ-Trp has no specific "binding site"-type interaction with mono- and dimeric Aβ, which could explain its low inhibitory efficiency. This suggests that the reported anti-AD activity of NQ-Trp-type molecules in in vivo models has to involve another mechanism. This study has revealed the potential pitfalls in the development of aggregation inhibitors for amyloidogenic peptides, which are of general interest for all the molecules studied in the context of inhibiting the formation of toxic aggregates.

Effect of pH on the Aggregation of α-syn12 Dimer in Explicit Water by Replica-Exchange Molecular Dynamics Simulation

The dimeric structure of the N-terminal 12 residues drives the interaction of α-synuclein protein with membranes. Moreover, experimental studies indicated that the aggregation of α-synuclein is faster at low pH than neutral pH. Nevertheless, the effects of different pH on the structural characteristics of the α-syn12 dimer remain poorly understood. We performed 500 ns temperature replica exchange molecular dynamics (T-REMD) simulations of two α-syn12 peptides in explicit solvent. The free energy surfaces contain ten highly populated regions at physiological pH, while there are only three highly populated regions contained at acidic pH. The anti-parallel β-sheet conformations were found as the lowest free energy state. Additionally, these states are nearly flat with a very small barrier which indicates that these states can easily transit between themselves. The dimer undergoes a disorder to order transition from physiological pH to acidic pH and the α-syn12 dimer at acidic pH involves a faster dimerization process. Further, the Lys6–Asp2 contact may prevent the dimerization.

Effects of the Pathogenic Mutation A117V and the Protective Mutation H111S on the Folding and Aggregation of PrP106-126: Insights from Replica Exchange Molecular Dynamics Simulations

The fragment 106-126 of prion protein exhibits similar properties to full-length prion. Experiments have shown that the A117V mutation enhances the aggregation of PrP106-126, while the H111S mutation abolishes the assembly. However, the mechanism of the change in the aggregation behavior of PrP106-126 upon the two mutations is not fully understood. In this study, replica exchange molecular dynamics simulations were performed to investigate the conformational ensemble of the WT PrP106-126 and its two mutants A117V and H111S. The obtained results indicate that the three species are all intrinsically disordered but they have distinct morphological differences. The A117V mutant has a higher propensity to form β-hairpin structures than the WT, while the H111S mutant has a higher population of helical structures. Furthermore, the A117V mutation increases the hydrophobic solvent accessible surface areas of PrP106-126 and the H111S mutation reduces the exposure of hydrophobic residues. It can be concluded that the difference in populations of β-hairpin structures and the change of hydrophobic solvent accessible areas may induce the different aggregation behaviors of the A117V and the H111S mutated PrP106-126. Understanding why the two mutations have contrary effects on the aggregation of PrP106-126 is very meaningful for further elucidation of the mechanism underlying aggregation and design of inhibitor against aggregation process.

Cross dimerization of amyloid-β and αsynuclein proteins in aqueous environment: a molecular dynamics simulations study

Self-assembly of the intrinsically unstructured proteins, amyloid beta (Aβ) and alpha synclein (αSyn), are associated with Alzheimer's Disease, and Parkinson's and Lewy Body Diseases, respectively. Importantly, pathological overlaps between these neurodegenerative diseases, and the possibilities of interactions between Aβ and αSyn in biological milieu emerge from several recent clinical reports and in vitro studies. Nevertheless, there are very few molecular level studies that have probed the nature of spontaneous interactions between these two sequentially dissimilar proteins and key characteristics of the resulting cross complexes. In this study, we have used atomistic molecular dynamics simulations to probe the possibility of cross dimerization between αSyn1-95 and Aβ1-42, and thereby gain insights into their plausible early assembly pathways in aqueous environment. Our analyses indicate a strong probability of association between the two sequences, with inter-protein attractive electrostatic interactions playing dominant roles. Principal component analysis revealed significant heterogeneity in the strength and nature of the associations in the key interaction modes. In most, the interactions of repeating Lys residues, mainly in the imperfect repeats 'KTKEGV' present in αSyn1-95 were found to be essential for cross interactions and formation of inter-protein salt bridges. Additionally, a hydrophobicity driven interaction mode devoid of salt bridges, where the non-amyloid component (NAC) region of αSyn1-95 came in contact with the hydrophobic core of Aβ1-42 was observed. The existence of such hetero complexes, and therefore hetero assembly pathways may lead to polymorphic aggregates with variations in pathological attributes. Our results provide a perspective on development of therapeutic strategies for preventing pathogenic interactions between these proteins.