Effect of the Tottori familial disease mutation (D7N) on the monomers and dimers of Aβ40 and Aβ42

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.

Fe(2+) binding on amyloid β-peptide promotes aggregation

The metal ions Zn(2+) , Cu(2+) , and Fe(2+) play a significant role in the aggregation mechanism of Aβ peptides. However, the nature of binding between metal and peptide has remained elusive; the detailed information on this from the experimental study is very difficult. Density functional theory (dft) (M06-2X/6-311++G (2df,2pd) +LANL2DZ) has employed to determine the force field resulting due to metal and histidine interaction. We performed 200 ns molecular dynamics (MD) simulation on Aβ1-42 -Zn(2+) , Aβ1-42 -Cu(2+) , and Aβ1-42 -Fe(2+) systems in explicit water with different combination of coordinating residues including the three Histidine residues in the N-terminal. The present investigation, the Aβ1-42 -Zn(2+) system possess three turn conformations separated by coil structure. Zn(2+) binding caused the loss of the helical structure of N-terminal residues which transformed into the S-shaped conformation. Zn(2+) has reduced the coil and increases the turn content of the peptide compared with experimental study. On the other hand, the Cu(2+) binds with peptide, β sheet formation is observed at the N-terminal residues of the peptide. Fe(2+) binding is to promote the formation of Glu22-Lys28 salt-bridge which stabilized the turn conformation in the Phe19-Gly25 residues, subsequently β sheets were observed at His13-Lys18 and Gly29-Gly37 residues. The turn conformation facilitates the β sheets are arranged in parallel by enhancing the hydrophobic contact between Gly25 and Met35, Lys16 and Met35, Leu17 and Leu34, Val18 and Leu34 residues. The Fe(2+) binding reduced the helix structure and increases the β sheet content in the peptide, which suggested, Fe(2+) promotes the oligomerization by enhancing the peptide-peptide interaction. Proteins 2016; 84:1257-1274. © 2016 Wiley Periodicals, Inc.

Structure of ring-shaped Aβ₄₂ oligomers determined by conformational selection

The oligomerization of amyloid beta (Aβ) peptides into soluble non-fibrillar species plays a critical role in the pathogenesis of Alzheimer’s disease. However, it has been challenging to characterize the tertiary and quaternary structures of Aβ peptides due to their disordered nature and high aggregation propensity. In this work, replica exchange molecular dynamics simulations were used to explore the conformational space of Aβ₄₂ monomer. Among the most populated transient states, we identified a particular conformation which was able to generate ring-shaped pentamers and hexamers, when docked onto itself. The structures of these aggregates were stable during microsecond all-atom MD simulations in explicit solvent. In addition to high resolution models of these oligomers, this study provides support for the conformational selection mechanism of Aβ peptide self-assembly.

Allosteric stabilization of the amyloid-β peptide hairpin by the fluctuating N-terminal

Immobilized ions modulate nearby hydrophobic interactions and influence molecular recognition and self-assembly. We simulated disulfide bond-locked double mutants (L17C/L34C) and observed allosteric modulation of the peptide's intra-molecular interactions by the N-terminal tail. We revealed that the non-contacting charged N-terminal residues help the transfer of entropy to the surrounding solvation shell and stabilizing β-hairpin.

Alzheimer's protective A2T mutation changes the conformational landscape of the Aβ₁₋₄₂ monomer differently than does the A2V mutation

The aggregation of amyloid-β (Aβ) peptides plays a crucial role in the etiology of Alzheimer's disease (AD). Recently, it has been reported that an A2T mutation in Aβ can protect against AD. Interestingly, a nonpolar A2V mutation also has been found to offer protection against AD in the heterozygous state, although it causes early-onset AD in homozygous carriers. Since the conformational landscape of the Aβ monomer is known to directly contribute to the early-stage aggregation mechanism, it is important to characterize the effects of the A2T and A2V mutations on Aβ₁₋₄₂ monomer structure. Here, we have performed extensive atomistic replica-exchange molecular dynamics simulations of the solvated wild-type (WT), A2V, and A2T Aβ₁₋₄₂ monomers. Our simulations reveal that although all three variants remain as collapsed coils in solution, there exist significant structural differences among them at shorter timescales. A2V exhibits an enhanced double-hairpin population in comparison to the WT, similar to those reported in toxic WT Aβ₁₋₄₂ oligomers. Such double-hairpin formation is caused by hydrophobic clustering between the N-terminus and the central and C-terminal hydrophobic patches. In contrast, the A2T mutation causes the N-terminus to engage in unusual electrostatic interactions with distant residues, such as K16 and E22, resulting in a unique population comprising only the C-terminal hairpin. These findings imply that a single A2X (where X = V or T) mutation in the primarily disordered N-terminus of the Aβ₁₋₄₂ monomer can dramatically alter the β-hairpin population and switch the equilibrium toward alternative structures. The atomistically detailed, comparative view of the structural landscapes of A2V and A2T variant monomers obtained in this study can enhance our understanding of the mechanistic differences in their early-stage aggregation.

The Dependence of Amyloid-β Dynamics on Protein Force Fields and Water Models

We studied the dynamics of Aβ40 , involved in Alzheimer's disease, by using 21 methods combined from Amber03, Amber99sb-ILDN, Charmm27, Charmm22*, OPLS-2001, OPLS-2006, OPLS-2008, Gromos96-43a1, Gromos96-53a6, Gromos96-54a7, and the water models SPC, TIP3P, TIP4P. Major differences in the structural ensembles were systematized: Amber03, Charmm27, and Gromos96-54a7 stabilize the helices; Gromos96-43a1 and Gromos53a6 favor the β-strands (with Charmm22* and Amber99sb-ILDN in between), and OPLS produces unstructured ensembles. The accuracy of the NMR chemical shifts was in the order: Charmm22*>Amber99sb-ILDN>OPLS-2008≈Gromos96-43a1>Gromos96-54a7≈OPLS-2001>OPLS-2006>Gromos96-53a6>Charmm27>Amber03. The computed (3) JHNHα -coupling constants were sensitive to experiment type and Karplus parameterization. Overall, the ensembles of Charmm22* and Amber99sb-ILDN provided the best agreement with experimental NMR and circular dichroism data, providing a model for the real Aβ monomer ensemble. Also, the polar water model TIP3P significantly favored helix and compact conformations.

Polymorphic cross-seeding amyloid assemblies of amyloid-β and human islet amyloid polypeptide

Epidemiological studies have shown that the development of Alzheimer's disease (AD) is associated with type 2 diabetes (T2D), but it still remains unclear how AD and T2D are connected. Heterologous cross-seeding between the causative peptides of Aβ and hIAPP may represent a molecular link between AD and T2D. Here, we computationally modeled and simulated a series of cross-seeding double-layer assemblies formed by Aβ and hIAPP peptides using all-atom and coarse-gained molecular dynamics (MD) simulations. The cross-seeding Aβ-hIAPP assemblies showed a wide range of polymorphic structures via a combination of four β-sheet-to-β-sheet interfaces and two packing orientations, focusing on a comparison of different matches of β-sheet layers. Two cross-seeding Aβ-hIAPP assemblies with different interfacial β-sheet packings exhibited high structural stability and favorable interfacial interactions in both oligomeric and fibrillar states. Both Aβ-hIAPP assemblies displayed interfacial dehydration to different extents, which in turn promoted Aβ-hIAPP association depending on interfacial polarity and geometry. Furthermore, computational mutagenesis studies revealed that disruption of interfacial salt bridges largely disfavor the β-sheet-to-β-sheet association, highlighting the importance of salt bridges in the formation of cross-seeding assemblies. This work provides atomic-level information on the cross-seeding interactions between Aβ and hIAPP, which may be involved in the interplay between these two disorders.

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.

Exploring the Alzheimer amyloid-β peptide conformational ensemble: A review of molecular dynamics approaches

Alzheimer's disease is one of the most common dementia among elderly worldwide. There is no therapeutic drugs until now to treat effectively this disease. One main reason is due to the poorly understood mechanism of Aβ peptide aggregation, which plays a crucial role in the development of Alzheimer's disease. It remains challenging to experimentally or theoretically characterize the secondary and tertiary structures of the Aβ monomer because of its high flexibility and aggregation propensity, and its conformations that lead to the aggregation are not fully identified. In this review, we highlight various structural ensembles of Aβ peptide revealed and characterized by computational approaches in order to find converging structures of Aβ monomer. Understanding how Aβ peptide forms transiently stable structures prior to aggregation will contribute to the design of new therapeutic molecules against the Alzheimer's disease.

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.

Coupling of Zinc-Binding and Secondary Structure in Nonfibrillar Aβ40 Peptide Oligomerization

Nonfibrillar neurotoxic amyloid β (Aβ) oligomer structures are typically rich in β-sheets, which could be promoted by metal ions like Zn(2+). Here, using molecular dynamics (MD) simulations, we systematically examined combinations of Aβ40 peptide conformations and Zn(2+) binding modes to probe the effects of secondary structure on Aβ dimerization energies and kinetics. We found that random conformations do not contribute to dimerization either thermodynamically or kinetically. Zn(2+) couples with preformed secondary structures (α-helix and β-hairpin) to speed dimerization and stabilize the resulting dimer. Partial α-helices increase the dimerization speed, and dimers with α-helix rich conformations have the lowest energy. When Zn(2+) coordinates with residues D1, H6, H13, and H14, Aβ40 β-hairpin monomers have the fastest dimerization speed. Dimers with experimentally observed zinc coordination (E11, H6, H13, and H14) form with slower rate but have lower energy. Zn(2+) cannot stabilize fibril-like β-arch dimers. However, Zn(2+)-bound β-arch tetramers have the lowest energy. Collectively, zinc-stabilized β-hairpin oligomers could be important in the nucleation-polymerization of cross-β structures. Our results are consistent with experimental findings that α-helix to β-structural transition should accompany Aβ aggregation in the presence of zinc ions and that Zn(2+) stabilizes nonfibrillar Aβ oligomers and, thus, inhibits formation of less toxic Aβ fibrils.

Folding Atomistic Proteins in Explicit Solvent Using Simulated Tempering

Following a previous report on a coarse-grained protein model in implicit solvent, we applied simulated tempering (ST) with on-the-fly Helmholtz free energy (weight factors) determination to the folding or aggregation of seven proteins with the CHARMM, OPLS, and AMBER protein, and the SPC and TIP3P water force fields. For efficiency and reliability, we also performed replica exchange molecular dynamics (REMD) simulations on the alanine di- and deca-peptide, and the dimer of the Aβ16-22 Alzheimer's fragment, and used experimental data and previous simulation results on the chignolin, beta3s, Trp-cage, and WW domain peptides of 10-37 amino acids. The sampling with ST is found to be more efficient than with REMD for a much lower CPU cost. Starting from unfolded or extended conformations, the WW domain and the Trp-cage peptide fold to their NMR structures with a backbone RMSD of 2.0 and 1 Å. Remarkably, the ST simulation explores transient non-native topologies for Trp-cage that have been rarely discussed by other simulations. Our ST simulations also show that the CHARMM22* force field has limitations in describing accurately the beta3s peptide. Taken together, these results open the door to the study of the configurations of single proteins, protein aggregates, and any molecular systems at atomic details in explicit solvent using a single normal CPU. They also demonstrate that our ST scheme can be used with any force field ranging from quantum mechanics to coarse-grain and atomistic.

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.

Stability of amyloid oligomers

Molecular simulations are now commonly used to complement experimental techniques in investigating amyloids and their role in human diseases. In this chapter, we will summarize techniques and approaches often used in amyloid simulations and will present recent success stories. Our examples will be focused on lessons learned from molecular dynamics simulations in aqueous environments that start from preformed aggregates. These studies explore the limitations that arise from the choice of force field, the role of mutations in the growth of amyloid aggregates, segmental polymorphism, and the importance of cross-seeding. Furthermore, they give evidence for potential toxicity mechanisms. We finally discuss the role of molecular simulations in the search for aggregation inhibitors.

Effect of the English familial disease mutation (H6R) on the monomers and dimers of Aβ40 and Aβ42

The self-assembly of the amyloid beta (Aβ) peptides into senile plaques is the hallmark of Alzheimer's disease. Recent experiments have shown that the English familial disease mutation (H6R) speeds up the fibril formation process of alloforms Aβ40 and Aβ42 peptides altering their toxicity to cells. We used all-atom molecular dynamics simulations at microsecond time scales with the OPLS-AA force field and TIP4P explicit water model to study the structural dynamics of the monomer and dimer of H6R sequences of both peptides. The reason behind the self-assembly acceleration is common that upon mutation the net charge is reduced leading to the weaker repulsive interaction between chains that facilitates the peptide association. In addition, our estimation of the solvation free energy shows that the mutation enhances the hydrophobicity of both peptides speeding up their aggregation. However, we can show that the acceleration mechanisms are different for different peptides: the rate of fibril formation of Aβ42 increases due to increased β-structure at the C-terminal in both monomer and dimer and enhanced stability of salt bridge Asp23-Lys28 in monomer, while the enhancement of turn at residues 25-29 and reduction of coil in regions 10-13, 26-19, and 30-34 would play the key role for Aβ40. Overall, our study provides a detailed atomistic picture of the H6R-mediated conformational changes that are consistent with the experimental findings and highlights the important role of the N-terminal in Aβ peptide aggregation.

Alzheimer's disease--a panorama glimpse

The single-mutation of genes associated with Alzheimer's disease (AD) increases the production of Aβ peptides. An elevated concentration of Aβ peptides is prone to aggregation into oligomers and further deposition as plaque. Aβ plaques and neurofibrillary tangles are two hallmarks of AD. In this review, we provide a broad overview of the diverses sources that could lead to AD, which include genetic origins, Aβ peptides and tau protein. We shall discuss on tau protein and tau accumulation, which result in neurofibrillary tangles. We detail the mechanisms of Aβ aggregation, fibril formation and its polymorphism. We then show the possible links between Aβ and tau pathology. Furthermore, we summarize the structural data of Aβ and its precursor protein obtained via Nuclear Magnetic Resonance (NMR) or X-ray crystallography. At the end, we go through the C-terminal and N-terminal truncated Aβ variants. We wish to draw reader's attention to two predominant and toxic Aβ species, namely Aβ4-42 and pyroglutamate amyloid-beta peptides, which have been neglected for more than a decade and may be crucial in Aβ pathogenesis due to their dominant presence in the AD brain.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Molecular structure of the NQTrp inhibitor with the Alzheimer Aβ1-28 monomer

The self-assembly of the amyloid-β (Aβ) peptide of various amino acid lengths into senile plaques is one hallmark of Alzheimer's disease pathology. In the past decade, many small molecules, including NQTrp, have been identified to reduce aggregation and toxicity. However, due to the heterogeneity of the conformational ensemble of Aβ with drugs, we lack detailed structures of the transient complexes. Following our previous simulation of the monomer of Aβ1-28, here we characterize the equilibrium ensemble of the Aβ1-28 monomer with NQTrp by means of extensive atomistic replica exchange molecular dynamics simulations using a force field known to fold diverse proteins correctly. While the secondary structure content and the intrinsic disorder of the whole peptides are very similar and the lifetimes of the salt-bridges remain constant, the population of β-hairpin is reduced by a factor of 1.5 and the population of α-helix in the region 17-24 is increased by a factor of two upon NQTrp binding. These two factors, which impact the free energy barrier for nucleation, provide a first explanation for the reported reduced Aβ1-40/1-42 aggregation kinetics in the presence of NQTrp. Backbone and side-chain interactions of Aβ with NQTrp may also inhibit Aβ-Aβ contacts. The fraction of free Aβ1-28 monomer is, however, on the order of 20-25% at 17.5 mM, and this shows that the affinity of NQTrp is low and hence its inhibitory activity is not very strong. This inhibitor can be improved to reduce the formation of dimer, a critical step in aggregation and toxicity.

Significant structural differences between transient amyloid-β oligomers and less-toxic fibrils in regions known to harbor familial Alzheimer's mutations

Small oligomers of the amyloid β (Aβ) peptide, rather than the monomers or the fibrils, are suspected to initiate Alzheimer's disease (AD). However, their low concentration and transient nature under physiological conditions have made structural investigations difficult. A method for addressing such problems has been developed by combining rapid fluorescence techniques with slower two-dimensional solid-state NMR methods. The smallest Aβ40 oligomers that demonstrate a potential sign of toxicity, namely, an enhanced affinity for cell membranes, were thus probed. The two hydrophobic regions (residues 10-21 and 30-40) have already attained the conformation that is observed in the fibrils. However, the turn region (residues 22-29) and the N-terminal tail (residues 1-9) are strikingly different. Notably, ten of eleven known Aβ mutants that are linked to familial AD map to these two regions. Our results provide potential structural cues for AD therapeutics and also suggest a general method for determining transient protein structures.

The conformational stability of nonfibrillar amyloid-β peptide oligomers critically depends on the C-terminal peptide length

The amyloid-β (Aβ) peptide is one key molecule in the pathogenesis of Alzheimer's disease. We investigated the conformational stability of a nonfibrillar tetrameric Aβ structure by molecular dynamics (MD) simulations revealing that the stability of the Aβ tetramer depends critically on the C-terminal length. In contrast to the Aβ17-40 tetramer, which proved to be instable, the simulations demonstrate structural integrity of the Aβ17-42 and Aβ17-43 tetramers. These differences in stability can be attributed to an extension of the middle strand of a three-stranded antiparallel β sheet through residues 41-43, only present in the longer Aβ species that aggregate faster and are more neurotoxic. Additional MD simulations demonstrate that this higher stability is also present in the monomers forming the tetramer. In conclusion, our findings suggest the existence of a nonfibrillar oligomer topology that is significantly more stable for the longer Aβ species, thus offering a structural explanation for their higher neurotoxicity.

Atomic and dynamic insights into the beneficial effect of the 1,4-naphthoquinon-2-yl-L-tryptophan inhibitor on Alzheimer's Aβ1-42 dimer in terms of aggregation and toxicity

Aggregation of the amyloid β protein (Aβ) peptide with 40 or 42 residues is one key feature in Alzheimer's disease (AD). The 1,4-naphthoquinon-2-yl-L-tryptophan (NQTrp) molecule was reported to alter Aβ self-assembly and reduce toxicity. Though nuclear magnetic resonance experiments and various simulations provided atomic information about the interaction of NQTrp with Aβ peptides spanning the regions of residues 12-28 and 17-42, none of these studies were conducted on the full-length Aβ1-42 peptide. To this end, we performed extensive atomistic replica exchange molecular dynamics simulations of Aβ1-42 dimer with two NQTrp molecules in explicit solvent, by using a force field known to fold diverse proteins correctly. The interactions between NQTrp and Aβ1-42, which change the Aβ interface by reducing most of the intermolecular contacts, are found to be very dynamic and multiple, leading to many transient binding sites. The most favorable binding residues are Arg5, Asp7, Tyr10, His13, Lys16, Lys18, Phe19/Phe20, and Leu34/Met35, providing therefore a completely different picture from in vitro and in silico experiments with NQTrp with shorter Aβ fragments. Importantly, the 10 hot residues that we identified explain the beneficial effect of NQTrp in reducing both the level of Aβ1-42 aggregation and toxicity. Our results also indicate that there is room to design more efficient drugs targeting Aβ1-42 dimer against AD.