Dimer formation enhances structural differences between amyloid β-protein (1-40) and (1-42): an explicit-solvent molecular dynamics study

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.

Aβ42 and Aβ40: similarities and differences

The abnormal accumulation of amyloid-β (Aβ) peptide in the brain is one of the most important hallmarks of Alzheimer's disease. Aβ is an aggregation-prone and toxic polypeptide with 39-43 residues, derived from the amyloid precursor protein proteolysis process. According to the amyloid hypothesis, abnormal accumulation of Aβ in the brain is the primary influence driving Alzheimer's disease pathologies. Among all kinds of Aβ isoforms, Aβ40 and Aβ42 are believed to be the most important ones. Although these two kinds of Aβ differ only in two amino acid residues, recent studies show that they differ significantly in their metabolism, physiological functions, toxicities, and aggregation mechanism. In this review, we mainly summarize the similarities and differences between Aβ42 and Aβ40, recent studies on selective inhibitors as well as probes will also be mentioned.

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.

Disordered amyloidogenic peptides may insert into the membrane and assemble into common cyclic structural motifs

Aggregation of disordered amyloidogenic peptides into oligomers is the causative agent of amyloid-related diseases. In solution, disordered protein states are characterized by heterogeneous ensembles. Among these, β-rich conformers self-assemble via a conformational selection mechanism to form energetically-favored cross-β structures, regardless of their precise sequences. These disordered peptides can also penetrate the membrane, and electrophysiological data indicate that they form ion-conducting channels. Based on these and additional data, including imaging and molecular dynamics simulations of a range of amyloid peptides, Alzheimer's amyloid-β (Aβ) peptide, its disease-related variants with point mutations and N-terminal truncated species, other amyloidogenic peptides, as well as a cytolytic peptide and a synthetic gel-forming peptide, we suggest that disordered amyloidogenic peptides can also present a common motif in the membrane. The motif consists of curved, moon-like β-rich oligomers associated into annular organizations. The motif is favored in the lipid bilayer since it permits hydrophobic side chains to face and interact with the membrane and the charged/polar residues to face the solvated channel pores. Such channels are toxic since their pores allow uncontrolled leakage of ions into/out of the cell, destabilizing cellular ionic homeostasis. Here we detail Aβ, whose aggregation is associated with Alzheimer's disease (AD) and for which there are the most abundant data. AD is a protein misfolding disease characterized by a build-up of Aβ peptide as senile plaques, neurodegeneration, and memory loss. Excessively produced Aβ peptides may directly induce cellular toxicity, even without the involvement of membrane receptors through Aβ peptide-plasma membrane interactions.

Structural diversity of Alzheimer's disease amyloid-β dimers and their role in oligomerization and fibril formation

Alzheimer's disease (AD) is associated with the formation of toxic amyloid-β (Aβ)42 oligomers, and recent evidence supports a role for Aβ dimers as building blocks for oligomers. Molecular dynamics simulation studies have identified clans for the dominant conformations of Aβ42 forming dimers; however, it is unclear if a larger spectrum of dimers is involved and which set(s) of dimers might evolve to oligomers verse fibrils. Therefore, for this study we generated multiple structural conformations of Aβ42, using explicit all-atom molecular dynamics, and then clustering the different structures based on key conformational similarities. Those matching a selection threshold were then used to model a process of oligomerization. Remarkably, we showed a greater diversity in Aβ dimers than previously described. Depending on the clan family, different types of Aβ dimers were obtained. While some had the tendency to evolve into oligomeric rings, others formed fibrils of diverse characteristics. Then we selected the dimers that would evolve to membranephilic annular oligomers. Nearly one third of the 28 evaluated annular oligomers had the dimer interfaces between the neighboring Aβ42 monomers with possible salt bridges between the residue K28 from one side and either residue E22 or D23 on the other. Based on these results, key amino acids were identified for point mutations that either enhanced or suppressed the formation and toxicity of oligomer rings. Our studies suggest a greater diversity of Aβ dimers. Understanding the structure of Aβ dimers might be important for the rationale design of small molecules that block formation of toxic oligomers.

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.

Discrete molecular dynamics can predict helical prestructured motifs in disordered proteins

Intrinsically disordered proteins (IDPs) lack a stable tertiary structure, but their short binding regions termed Pre-Structured Motifs (PreSMo) can form transient secondary structure elements in solution. Although disordered proteins are crucial in many biological processes and designing strategies to modulate their function is highly important, both experimental and computational tools to describe their conformational ensembles and the initial steps of folding are sparse. Here we report that discrete molecular dynamics (DMD) simulations combined with replica exchange (RX) method efficiently samples the conformational space and detects regions populating α-helical conformational states in disordered protein regions. While the available computational methods predict secondary structural propensities in IDPs based on the observation of protein-protein interactions, our ab initio method rests on physical principles of protein folding and dynamics. We show that RX-DMD predicts α-PreSMos with high confidence confirmed by comparison to experimental NMR data. Moreover, the method also can dissect α-PreSMos in close vicinity to each other and indicate helix stability. Importantly, simulations with disordered regions forming helices in X-ray structures of complexes indicate that a preformed helix is frequently the binding element itself, while in other cases it may have a role in initiating the binding process. Our results indicate that RX-DMD provides a breakthrough in the structural and dynamical characterization of disordered proteins by generating the structural ensembles of IDPs even when experimental data are not available.

Comparison of three competing dynamic force spectroscopy models to study binding forces of amyloid-β (1-42)

We performed single molecule dynamic force spectroscopy experiments to study the dimerization of two amyloid-β (1-42) peptides and compared three different theoretical models used to fit experimental data: Bell-Evans, Dudko-Hummer-Szabo, and Friddle-De Yoreo. Using these models we extracted values of the dissociation rate at zero force, k0, and height and the width of the energy barrier, ΔG and xβ. We show the importance of including the effect of the linker molecule. All three models corrected for the linker effect give comparable results for xβ and show more discrepancy for k0 and ΔG values, ΔG parameter correlates well between Dudko-Hummer-Szabo and Friddle-De Yoreo models but differs for the Bell-Evans model.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

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.

A kinetic approach to the sequence-aggregation relationship in disease-related protein assembly

It is generally accepted that oligomers of aggregating proteins play an important role in the onset of neurodegenerative diseases. While in silico aggregation studies of full length amyloidogenic proteins are computationally expensive, the assembly of short protein fragments derived from these proteins with similar aggregating properties has been extensively studied. In the present work, molecular dynamics simulations are performed to follow peptide aggregation on the microsecond time scale. By defining aggregation states, we identify transition networks, disconnectivity graphs, and first passage time distributions to describe the kinetics of the assembly process. This approach unravels differences in the aggregation into hexamers of two peptides with different primary structures. The first is GNNQQNY, a hydrophilic fragment from the prion protein Sup35, and the second is KLVFFAE, a fragment from amyloid-β protein, with a hydrophobic core delimited by two charged amino acids. The assembly of GNNQQNY suggests a mechanism of monomer addition, with a bias toward parallel peptide pairs and a gradual increase in the amount of β-strand content. For KLVFFAE, a mechanism involving dimers rather than monomers is revealed, involving a generally higher β-strand content and a transition toward a larger number of antiparallel peptide pairs during the rearrangement of the hexamer. The differences observed for the aggregation of the two peptides suggests the existence of a sequence-aggregation relationship.

Mechanical resistance in unstructured proteins

Single-molecule pulling experiments on unstructured proteins linked to neurodegenerative diseases have measured rupture forces comparable to those for stable folded proteins. To investigate the structural mechanisms of this unexpected force resistance, we perform pulling simulations of the amyloid β-peptide (Aβ) and α-synuclein (αS), starting from simulated conformational ensembles for the free monomers. For both proteins, the simulations yield a set of rupture events that agree well with the experimental data. By analyzing the conformations occurring shortly before rupture in each event, we find that the mechanically resistant structures share a common architecture, with similarities to the folds adopted by Aβ and αS in amyloid fibrils. The disease-linked Arctic mutation of Aβ is found to increase the occurrence of highly force-resistant structures. Our study suggests that the high rupture forces observed in Aβ and αS pulling experiments are caused by structures that might have a key role in amyloid formation.

Effect of metals on kinetic pathways of amyloid-β aggregation

Metal ions, including copper and zinc, have been implicated in the pathogenesis of Alzheimer’s disease through a variety of mechanisms including increased amyloid-β affinity and redox effects. Recent reports have demonstrated that the amyloid-β monomer does not necessarily travel through a definitive intermediary en-route to a stable amyloid fibril structure. Rather, amyloid-β misfolding may follow a variety of pathways resulting in a fibrillar end-product or a variety of oligomeric end-products with a diversity of structures and sizes. The presence of metal ions has been demonstrated to alter the kinetic pathway of the amyloid-β peptide which may lead to more toxic oligomeric end-products. In this work, we review the contemporary literature supporting the hypothesis that metal ions alter the reaction pathway of amyloid-β misfolding leading to more neurotoxic species.

Early amyloid β-protein aggregation precedes conformational change

The aggregation of amyloid-β protein (1–42) is studied at experimental concentrations using all-atom molecular dynamics simulations.

Spontaneous dimer states of the Aβ21–30decapeptide

Computational examination of the spontaneous dimerization of Aβ_21–30and stability measures of the resulting parallel and anti-parallel aligned dimers.

Amyloid β-protein oligomers and Alzheimer's disease

The oligomer cascade hypothesis, which states that oligomers are the initiating pathologic agents in Alzheimer’s disease, has all but supplanted the amyloid cascade hypothesis, which suggested that fibers were the key etiologic agents in Alzheimer’s disease. We review here the results of in vivo, in vitro and in silico studies of amyloid β-protein oligomers, and discuss important caveats that should be considered in the evaluation of these results. This article is divided into four sections that mirror the main approaches used in the field to better understand oligomers: (1) attempts to locate and examine oligomers in vivo in situ; that is, without removing these species from their environment; (2) studies involving oligomers extracted from human or animal tissues and the subsequent characterization of their properties ex vivo; (3) studies of oligomers that have been produced synthetically and studied using a reductionist approach in relatively simple in vitro biophysical systems; and (4) computational studies of oligomers in silico. These multiple orthogonal approaches have revealed much about the molecular and cell biology of amyloid β-protein. However, as informative as these approaches have been, the amyloid β-protein oligomer system remains enigmatic.

Effects of Zn2+ binding on the structural and dynamic properties of amyloid β peptide associated with Alzheimer's disease: Asp1 or Glu11?

Extensive experimental and computational studies have suggested that multiple Zn(2+) binding modes in amyloid β (Aβ) peptides could exist simultaneously. However, consistent results have not been obtained for the effects of Zn(2+) binding on Aβ structure, dynamics, and kinetics in particular. Some key questions such as why it is so difficult to distinguish the polymorphic states of metal ions binding to Aβ and what the underlying rationale is, necessitate elucidation. In this work, two 3N1O Zn(2+) binding modes were constructed with three histidines (His(6), His(13), and His(14)), and Asp(1)/Glu(11) of Aβ40 coordinated to Zn(2+). Results from molecular dynamics simulations reveal that the conformational ensembles of different Zn(2+)-Aβ40 complexes are nonoverlapping. The formation of turn structure and, especially, the salt bridge between Glu(22)/Asp(23) and Lys(28) is dependent on specific Zn(2+) binding mode. Agreement with available NMR observations of secondary and tertiary structures could be better achieved if the two simulation results are considered together. The free energy landscape constructed by combining both conformations of Aβ40 indicates that transitions between distinct Aβ40 conformations thar are ready for Zn(2+) binding could be possible in aqueous solution. Markov state model analyses reveal the complex network of conformational space of Aβ40 modeulated by Zn(2+) binding, suggesting various misfolding pathways. The binding free energies evaluated using a combination of quantum mechanics calculations and the MM/3D-RISM method suggest that Glu(11) is the preferred oxygen ligand of Zn(2+). However, such preference is dependent on the relative populations of different conformations with specific Zn(2+) binding modes, and therefore could be shifted when experimental or simulation conditions are altered.

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

Recent experiments have shown that the mutation Tottori (D7N) alters the toxicity, assembly and rate of fibril formation of the wild type (WT) amyloid beta (Aβ) Aβ40 and Aβ42 peptides. We used all-atom molecular dynamics simulations in explicit solvent of the monomer and dimer of both alloforms with their WT and D7N sequences. The monomer simulations starting from a random coil and totaling 3 μs show that the D7N mutation changes the fold and the network of salt bridges in both alloforms. The dimer simulations starting from the amyloid fibrillar states and totaling 4.4 μs also reveal noticeable changes in terms of secondary structure, salt bridge, and topology. Overall, this study provides physical insights into the enhanced rate of fibril formation upon D7N mutation and an atomic picture of the D7N-mediated conformational change on Aβ40 and Aβ42 peptides.

Mechanism of amyloid β-protein dimerization determined using single-molecule AFM force spectroscopy

Aβ42 and Aβ40 are the two primary alloforms of human amyloid β−protein (Aβ). The two additional C−terminal residues of Aβ42 result in elevated neurotoxicity compared with Aβ40, but the molecular mechanism underlying this effect remains unclear. Here, we used single−molecule force microscopy to characterize interpeptide interactions for Aβ42 and Aβ40 and corresponding mutants. We discovered a dramatic difference in the interaction patterns of Aβ42 and Aβ40 monomers within dimers. Although the sequence difference between the two peptides is at the C−termini, the N−terminal segment plays a key role in the peptide interaction in the dimers. This is an unexpected finding as N−terminal was considered as disordered segment with no effect on the Aβ peptide aggregation. These novel properties of Aβ proteins suggests that the stabilization of N−terminal interactions is a switch in redirecting of amyloids form the neurotoxic aggregation pathway, opening a novel avenue for the disease preventions and treatments.

Binding, conformational transition and dimerization of amyloid-β peptide on GM1-containing ternary membrane: insights from molecular dynamics simulation

Interactions of amyloid-β (Aβ) with neuronal membrane are associated with the progression of Alzheimer's disease (AD). Ganglioside GM1 has been shown to promote the structural conversion of Aβ and increase the rate of peptide aggregation; but the exact nature of interaction driving theses processes remains to be explored. In this work, we have carried out atomistic-scale computer simulations (totaling 2.65 µs) to investigate the behavior of Aβ monomer and dimers in GM1-containing raft-like membrane. The oligosaccharide head-group of GM1 was observed to act as scaffold for Aβ-binding through sugar-specific interactions. Starting from the initial helical peptide conformation, a β-hairpin motif was formed at the C-terminus of the GM1-bound Aβ-monomer; that didn't appear in absence of GM1 (both in fluid POPC and liquid-ordered cholesterol/POPC bilayers and also in aqueous medium) within the simulation time span. For Aβ-dimers, the β-structure was further enhanced by peptide-peptide interactions, which might influence the propensity of Aβ to aggregate into higher-ordered structures. The salt-bridges and inter-peptide hydrogen bonds were found to account for dimer stability. We observed spontaneous formation of intra-peptide D(23)-K(28) salt-bridge and a turn at V(24)GSN(27) region - long been accepted as characteristic structural-motifs for amyloid self-assembly. Altogether, our results provide atomistic details of Aβ-GM1 and Aβ-Aβ interactions and demonstrate their importance in the early-stages of GM1-mediated Aβ-oligomerisation on membrane surface.

Discrete molecular dynamics study of oligomer formation by N-terminally truncated amyloid β-protein

In Alzheimer's disease (AD), amyloid β-protein (Aβ) self-assembles into toxic oligomers. Of the two predominant Aβ alloforms, Aβ1-40 and Aβ1-42, the latter is particularly strongly linked to AD. N-terminally truncated and pyroglutamated Aβ peptides were recently shown to seed Aβ aggregation and contribute significantly to Aβ-mediated toxicity, yet their folding and assembly were not explored computationally. Discrete molecular dynamics approach previously captured in vitro-derived distinct Aβ1-40 and Aβ1-42 oligomer size distributions and predicted that the more toxic Aβ1-42 oligomers had more flexible and solvent-exposed N-termini than Aβ1-40 oligomers. Here, we examined oligomer formation of Aβ3-40, Aβ3-42, Aβ11-40, and Aβ11-42 by the discrete molecular dynamics approach. The four N-terminally truncated peptides showed increased oligomerization propensity relative to the full-length peptides, consistent with in vitro findings. Conformations formed by Aβ3-40/42 had significantly more flexible and solvent-exposed N-termini than Aβ1-40/42 conformations. In contrast, in Aβ11-40/42 conformations, the N-termini formed more contacts and were less accessible to the solvent. The compactness of the Aβ11-40/42 conformations was in part facilitated by Val12. Two single amino acid substitutions that reduced and abolished hydrophobicity at position 12, respectively, resulted in a proportionally increased structural variability. Our results suggest that Aβ11-40 and Aβ11-42 oligomers might be less toxic than Aβ1-40 and Aβ1-42 oligomers and offer a plausible explanation for the experimentally observed increased toxicity of Aβ3-40 and Aβ3-42 and their pyroglutamated forms.

Molecular mechanism of the inhibition of EGCG on the Alzheimer Aβ(1-42) dimer

Growing evidence supports that amyloid β (Aβ) oligomers are the major causative agents leading to neural cell death in Alzheimer's disease. The polyphenol (-)-epigallocatechin gallate (EGCG) was recently reported to inhibit Aβ fibrillization and redirect Aβ aggregation into unstructured, off-pathway oligomers. Given the experimental challenge to characterize the structures of Aβ/EGCG complexes, we performed extensive atomistic replica exchange molecular dynamics simulations of Aβ1-42 dimer in the present and absence of EGCG in explicit solvent. Our equilibrium Aβ dimeric structures free of EGCG are consistent with the collision cross section from ion-mobility mass spectrometry and the secondary structure composition from circular dichroism experiment. In the presence of EGCG, the Aβ structures are characterized by increased inter-center-of-mass distances, reduced interchain and intrachain contacts, reduced β-sheet content, and increased coil and α-helix contents. Analysis of the free energy surfaces reveals that the Aβ dimer with EGCG adopts new conformations, affecting therefore its propensity to adopt fibril-prone states. Overall, this study provides, for the first time, insights on the equilibrium structures of Aβ1-42 dimer in explicit aqueous solution and an atomic picture of the EGCG-mediated conformational change on Aβ dimer.

Cu(2+) affects amyloid-β (1-42) aggregation by increasing peptide-peptide binding forces

The link between metals, Alzheimer's disease (AD) and its implicated protein, amyloid-β (Aβ), is complex and highly studied. AD is believed to occur as a result of the misfolding and aggregation of Aβ. The dyshomeostasis of metal ions and their propensity to interact with Aβ has also been implicated in AD. In this work, we use single molecule atomic force spectroscopy to measure the rupture force required to dissociate two Aβ (1–42) peptides in the presence of copper ions, Cu²⁺. In addition, we use atomic force microscopy to resolve the aggregation of Aβ formed. Previous research has shown that metal ions decrease the lag time associated with Aβ aggregation. We show that with the addition of copper ions the unbinding force increases notably. This suggests that the reduction of lag time associated with Aβ aggregation occurs on a single molecule level as a result of an increase in binding forces during the very initial interactions between two Aβ peptides. We attribute these results to copper ions acting as a bridge between the two peptide molecules, increasing the stability of the peptide-peptide complex.

Comparative molecular dynamics study of human islet amyloid polypeptide (IAPP) and rat IAPP oligomers

Human islet amyloid polypeptide (hIAPP or amylin) is a causative agent in pancreatic amyloid deposits found in patients with type 2 diabetes. The aggregation of full-length hIAPP(1-37) into small oligomeric species is increasingly believed to be responsible for cell dysfunction and death. However, rat IAPP (rIAPP(1-37)), which differs from hIAPP in only six of 37 residues, loses its aggregation ability to form toxic amyloid species. Atomic details of the effect of sequence on the structure and toxicity between the amyloidogenic, toxic hIAPP peptide and the nonamyloidogenic, nontoxic rIAPP peptide remain unclear. Here, we probe sequence-induced differences in structural stability, conformational dynamics, and driving forces between different hIAPP and rIAPP polymorphic forms from monomer to pentamer using molecular dynamics simulations. Simulations show that hIAPP forms from trimer to pentamer exhibit high structural stability with well-preserved in-register parallel β-sheet and the U-bend conformation. The hIAPP trimer appears to be a smallest minimal seed in solution. The stabilities of parallel hIAPP oligomers increase with the number of peptides. Conversely, replacement of hIAPP sequence by rIAPP sequence causes a significant loss of favorable interpeptide interactions in all rIAPP oligomers, destabilizing the C-terminal β-sheet, turn conformation, and overall stability. A less β-sheet-rich structure and a disturbed U-shaped topology exert a large energy penalty on the self-assemble of the rIAPP peptides into highly ordered, in-register β-sheet-rich protofibrils and fibrils, which explains the nonamyloidogenic activity of rIAPP. Moreover, the absence of interior water within the U-turn region in the well-packed higher-order hIAPP oligomers, not in the poorly packed rIAPP oligomers, also stabilizes peptide association. This work provides atomic details of the sequence-structure relationship between the amyloidogenic hIAPP and its analogues such as the nonamyloidogenic rIAPP and some mutants, which could help in the development of novel therapeutic agents to block the formation of toxic hIAPP oligomeric species for type 2 diabetes.

Modeling amyloid-beta as homogeneous dodecamers and in complex with cellular prion protein

Soluble amyloid beta (Aβ) peptide has been linked to the pathology of Alzheimer's disease. A variety of soluble oligomers have been observed to be toxic, ranging from dimers to protofibrils. No tertiary structure has been identified as a single biologically relevant form, though many models are comprised of highly ordered β-sheets. Evidence exists for much less ordered toxic oligomers. The mechanism of toxicity remains highly debated and probably involves multiple pathways. Interaction of Aβ oligomers with the N-terminus of the cellular form of the prion protein (PrP(c)) has recently been proposed. The intrinsically disordered nature of this protein and the highly polymorphic nature of Aβ oligomers make structural resolution of the complex exceptionally challenging. In this study, molecular dynamics simulations are performed for dodecameric assemblies of Aβ comprised of monomers having a single, short antiparallel β-hairpin at the C-terminus. The resulting models, devoid of any intermolecular hydrogen bonds, are shown to correlate well with experimental data and are found to be quite stable within the hydrophobic core, whereas the α-helical N-termini transform to a random coil state. This indicates that highly ordered assemblies are not required for stability and less ordered oligomers are a viable component in the population of soluble oligomers. In addition, a tentative model is proposed for the association of Aβ dimers with a double deletion mutant of the intrinsically disordered N-terminus of PrP(c). This may be useful as a conceptual working model for the binding of higher order oligomers and in the design of further experiments.

Folding of pig gastric mucin non-glycosylated domains: a discrete molecular dynamics study

Mucin glycoproteins consist of tandem-repeating glycosylated regions flanked by non-repetitive protein domains with little glycosylation. These non-repetitive domains are involved in polymerization of mucin and play an important role in the pH-dependent gelation of gastric mucin, which is essential for protecting the stomach from autodigestion. We examine folding of the non-repetitive sequence of PGM-2X (242 amino acids) and the von Willebrand factor vWF-C1 domain (67 amino acids) at neutral and low pH using discrete molecular dynamics (DMD) in an implicit solvent combined with a four-bead peptide model. Using the same implicit solvent parameters, folding of both domains is simulated at neutral and low pH. In contrast to vWF-C1, PGM-2X folding is strongly affected by pH as indicated by changes in the contact order, radius of gyration, free-energy landscape, and the secondary structure. Whereas the free-energy landscape of vWF-C1 shows a single minimum at both neutral and low pH, the free-energy landscape of PGM-2X is characterized by multiple minima that are more numerous and shallower at low pH. Detailed structural analysis shows that PGM-2X partially unfolds at low pH. This partial unfolding is facilitated by the C-terminal region GLU236-PRO242, which loses contact with the rest of the domain due to effective "mean-field" repulsion among highly positively charged N- and C-terminal regions. Consequently, at low pH, hydrophobic amino acids are more exposed to the solvent. In vWF-C1, low pH induces some structural changes, including an increased exposure of CYS at position 67, but these changes are small compared to those found in PGM-2X. For PGM-2X, the DMD-derived average β-strand propensity increases from 0.26 ± 0.01 at neutral pH to 0.38 ± 0.01 at low pH. For vWF-C1, the DMD-derived average β-strand propensity is 0.32 ± 0.02 at neutral pH and 0.35 ± 0.02 at low pH. The DMD-derived structural information provides insight into pH-induced changes in the folding of two distinct mucin domains and suggests plausible mechanisms of the aggregation/gelation of mucin.

The toxicity of amyloid β oligomers

In this review, we elucidate the mechanisms of Aβ oligomer toxicity which may contribute to Alzheimer's disease (AD). In particular, we discuss on the interaction of Aβ oligomers with the membrane through the process of adsorption and insertion. Such interaction gives rises to phase transitions in the sub-structures of the Aβ peptide from α-helical to β-sheet structure. By means of a coarse-grained model, we exhibit the tendency of β-sheet structures to aggregate, thus providing further insights to the process of membrane induced aggregation. We show that the aggregated oligomer causes membrane invagination, which is a precursor to the formation of pore structures and ion channels. Other pathological progressions to AD due to Aβ oligomers are also covered, such as their interaction with the membrane receptors, and their direct versus indirect effects on oxidative stress and intraneuronal accumulation. We further illustrate that the molecule curcumin is a potential Aβ toxicity inhibitor as a β-sheet breaker by having a high propensity to interact with certain Aβ residues without binding to them. The comprehensive understanding gained from these current researches on the various toxicity mechanisms show promises in the provision of better therapeutics and treatment strategies in the near future.