Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment

Energetic contributions of residues to the formation of early amyloid-β oligomers

Low-weight amyloid-β (Aβ) oligomers formed at early stages of oligomerization rather than fibril assemblies seem to be the toxic components that drive neurodegeneration in Alzheimer's disease. Unfortunately, detailed knowledge of the structure of these early oligomers at the residue level is not yet available. In this study, we performed all-atom explicit solvent molecular dynamics simulations to examine the oligomerization process of Aβ10-35 monomers when forming dimers, trimers, tetramers and octamers, with four independent simulations of a total simulated time of 3 μs for each oligomer system. The decomposition of the stability free energy by MM-GBSA methodology allowed us to unravel the network of energetic interactions that stabilize such oligomers. The contribution of the intermonomeric van der Waals term is the most significant energy feature of the oligomerization process, consistent with the so-called hydrophobic effect. Furthermore, the decomposition of the stability free energy into residues and residue-pairwise terms revealed that it is mainly apolar interactions between the three specific hydrophobic fragments 31-35 (C-terminal region), 17-20 (central hydrophobic core) and 12-14 (N-terminal region) that are responsible for such a favourable effect. The conformation in which the hydrophobic cthr-chc interaction is oriented perpendicularly is particularly important. We propose three other model substructures that favour the oligomerization process and can thus be considered as molecular targets for future inhibitors. Understanding Aβ oligomerization at the residue level could lead to more efficient design of inhibitors of this process.

Interaction of amyloid inhibitor proteins with amyloid beta peptides: insight from molecular dynamics simulations

Knowledge of the detailed mechanism by which proteins such as human αB- crystallin and human lysozyme inhibit amyloid beta (Aβ) peptide aggregation is crucial for designing treatment for Alzheimer's disease. Thus, unconstrained, atomistic molecular dynamics simulations in explicit solvent have been performed to characterize the Aβ_17–42 assembly in presence of the αB-crystallin core domain and of lysozyme. Simulations reveal that both inhibitor proteins compete with inter-peptide interaction by binding to the peptides during the early stage of aggregation, which is consistent with their inhibitory action reported in experiments. However, the Aβ binding dynamics appear different for each inhibitor. The binding between crystallin and the peptide monomer, dominated by electrostatics, is relatively weak and transient due to the heterogeneous amino acid distribution of the inhibitor surface. The crystallin-bound Aβ oligomers are relatively long-lived, as they form more extensive contact surface with the inhibitor protein. In contrast, a high local density of arginines from lysozyme allows strong binding with Aβ peptide monomers, resulting in stable complexes. Our findings not only illustrate, in atomic detail, how the amyloid inhibitory mechanism of human αB-crystallin, a natural chaperone, is different from that of human lysozyme, but also may aid de novo design of amyloid inhibitors.

Tracking the mechanism of fibril assembly by simulated two-dimensional ultraviolet spectroscopy

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by the accumulation of plaque deposits in the human brain. The main component of these plaques consists of highly ordered structures called amyloid fibrils, formed by the amyloid β-peptide (Aβ). The mechanism connecting Aβ and AD is yet undetermined. In a previous study, a coarse-grained united-residue model and molecular dynamics simulations were used to model the growth mechanism of Aβ amyloid fibrils. On the basis of these simulations, a dock/lock mechanism was proposed, in which Aβ fibrils grow by adding monomers at either end of an amyloid fibril template. To examine the structures in the early time-scale formation and growth of amyloid fibrils, simulated two-dimensional ultraviolet spectroscopy is used. These early structures are monitored in the far ultraviolet regime (λ = 190-250 nm) in which the computed signals originate from the backbone nπ* and ππ* transitions. These signals show distinct cross-peak patterns that can be used, in combination with molecular dynamics, to monitor local dynamics and conformational changes in the secondary structure of Aβ-peptides. The protein geometry-correlated chiral xxxy signal and the non-chiral combined signal xyxy-xyyx were found to be sensitive to, and in agreement with, a dock/lock pathway.

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.

Synergistic applications of MD and NMR for the study of biological systems

Modern biological sciences are becoming more and more multidisciplinary. At the same time, theoretical and computational approaches gain in reliability and their field of application widens. In this short paper, we discuss recent advances in the areas of solution nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations that were made possible by the combination of both methods, that is, through their synergistic use. We present the main NMR observables and parameters that can be computed from simulations, and how they are used in a variety of complementary applications, including dynamics studies, model-free analysis, force field validation, and structural studies.

Probing the efficacy of peptide-based inhibitors against acid- and zinc-promoted oligomerization of amyloid-β peptide via single-oligomer spectroscopy

One avenue for prevention and treatment of Alzheimer's disease involves inhibiting the aggregation of amyloid-β peptide (Aβ). Given the deleterious effects reported for Aβ dimers and trimers, it is important to investigate inhibition of the earliest association steps. We have employed quantized photobleaching of dye-labeled Aβ peptides to characterize four peptide-based inhibitors of fibrillogenesis and/or cytotoxicity, assessing their ability to inhibit association in the smallest oligomers (n=2-5). Inhibitors were tested at acidic pH and in the presence of zinc, conditions that may promote oligomerization in vivo. Distributions of peptide species were constructed by examining dozens of surface-tethered monomers and oligomers, one at a time. Results show that all four inhibitors shift the distribution of Aβ species toward monomers; however, efficacies vary for each compound and sample environment. Collectively, these studies highlight promising design strategies for future oligomerization inhibitors, affording insight into oligomer structures and inhibition mechanisms in two physiologically significant environments.

Comparison of molecular dynamics simulation methods for amyloid β(1-42) monomers containing D-aspartic acid residues for predicting retention times in chromatography

Molecular dynamics simulations of amyloid β(1-42) containing D-aspartic acid residues were performed using several continuous solvent models to investigate the usefulness of simulation methods for D-amino acid-containing proteins and peptides. Normal molecular dynamics simulations and replica exchange molecular dynamics simulations, which are one of the generalized-ensemble algorithms, were performed. Because the β-structure contents of amyloid β(1-42) peptides obtained by replica exchange molecular dynamics simulations with Onufriev-Bashford-Case generalized Born implicit solvent were qualitatively consistent with experimental data, replica exchange molecular dynamics rather than other methods appeared to be more reasonable for calculations of amyloid β(1-42) containing D-aspartic acid residues. Computational results revealed that peptides with stereoinversion of Asp23 tend to form β-sheet structures by themselves, in contrast to the wild-type peptides that form β-sheet structures only after aggregation. These results are expected to be useful for computational investigations of proteins and peptides such as prediction of retention time of peptides and proteins containing D-aspartic acid residues.

The use of mass spectrometry to study amyloid-β peptides

Amyloid-β peptide (Aβ) varies in size from 39 to 43 amino acids and arises from sequential β- and γ-secretase processing of the amyloid precursor protein. Whereas the non-pathological role for Aβ is yet to be established, there is no disputing that Aβ is now widely regarded as central to the development of Alzheimer's disease (AD). The so named "amyloid cascade hypothesis" states that disease progression is the result of an increased Aβ burden in affected areas of the brain. To elucidate the Aβ role in AD, many analytical approaches have been proposed as suitable tools to investigate not only the total Aβ load but also many other issues that are considered crucial for AD, such as: (i) the aggregation state in which Aβ is present; (ii) its interaction with other species or metals; (iii) its ability to induce oxidative stress; and (iv) its degradative pathways. This review provides an insight into the use of mass spectrometry (MS) in the field of Aβ investigation aimed to assess its role in AD. In particular, the different MS-based approaches applied in vitro and in vivo that can provide detailed information on the above-mentioned issues are reviewed. Moreover, the advantages offered by the MS methods over all the other techniques are highlighted, together with the recent developments and uses of combined analytical approaches to detect and characterize Aβ.

Molecular dynamic simulation of wild type and mutants of the polymorphic amyloid NNQNTF segments of elk prion: structural stability and thermodynamic of association

A hexapeptide with amino acid sequence NNQNTF from the elk prion protein forms amyloid fibrils. Here we use molecular dynamic simulations of the oligomers and their single point glycine mutants to study their stability. In an effort to probe the structural stability and association thermodynamic in a realistic environment, all wildtype of NNQNTF polymorphic forms with different size and their corresponding double layer 5 strands single point glycine mutants were subjected to a total of 500 ns of explicit-solvent molecular dynamics (MD) simulation. Our results show that the structural stability of the NNQNTF oligomers increases with increasing the number of β-strands for double layers. Our results also demonstrated that hydrophobic interaction is the principle driving force to stabilize the adjacent β-strands while the steric zipper is responsible for holding the neighboring β-sheet layers together. We used MM-PBSA approach free energy calculations to determine the role of nonpolar effects, electrostatics and entropy in binding. Nonpolar effects remained consistently more favorable in wild type and mutants reinforcing the importance of hydrophobic effects in protein-protein binding. While entropy systematically opposed binding in all cases, there was no observed trend in the entropy difference between wildtype and glycine mutant. Free energy decomposition shows residues situated at the interface were found to make favorable contributions to the peptide-peptide association. The study of the wild type and mutants in an explicit solvent may provide valuable insight for amyloid aggregation inhibitor design efforts.

Membrane-targeted strategies for modulating APP and Abeta-mediated toxicity

Alzheimer's disease (AD) is a progressive neurodegenerative disease characterized by numerous pathological features including the accumulation of neurotoxic amyloid-β (Aβ) peptide. There is currently no effective therapy for AD, but the development of therapeutic strategies that target the cell membrane is gaining increased interest. The amyloid precursor protein (APP) from which Aβ is formed is a membrane-bound protein, and Aβ production and toxicity are both membrane mediated events. This review describes the critical role of cell membranes in AD with particular emphasis on how the composition and structure of the membrane and its specialized regions may influence toxic or benign Aβ/APP pathways in AD. The putative role of copper (Cu) in AD is also discussed, and we highlight how targeting the cell membrane with Cu complexes has therapeutic potential in AD.