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

A turn for the worse: Aβ β-hairpins in Alzheimer's disease

Amyloid-β (Aβ) oligomers are a cause of neurodegeneration in Alzheimer's disease (AD). These soluble aggregates of the Aβ peptide have proven difficult to study due to their inherent metastability and heterogeneity. Strategies to isolate and stabilize homogenous Aβ oligomer populations have emerged such as mutations, covalent cross-linking, and protein fusions. These strategies along with molecular dynamics simulations have provided a variety of proposed structures of Aβ oligomers, many of which consist of molecules of Aβ in β-hairpin conformations. β-Hairpins are intramolecular antiparallel β-sheets composed of two β-strands connected by a loop or turn. Three decades of research suggests that Aβ peptides form several different β-hairpin conformations, some of which are building blocks of toxic Aβ oligomers. The insights from these studies are currently being used to design anti-Aβ antibodies and vaccines to treat AD. Research suggests that antibody therapies designed to target oligomeric Aβ may be more successful at treating AD than antibodies designed to target linear epitopes of Aβ or fibrillar Aβ. Aβ β-hairpins are good epitopes to use in antibody development to selectively target oligomeric Aβ. This review summarizes the research on β-hairpins in Aβ peptides and discusses the relevance of this conformation in AD pathogenesis and drug development.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Resolution exchange with tunneling for enhanced sampling of protein landscapes

Simulations of protein folding and protein association happen on timescales that are orders of magnitude larger than what can typically be covered in all-atom molecular dynamics simulations. Use of low-resolution models alleviates this problem but may reduce the accuracy of the simulations. We introduce a replica-exchange-based multiscale sampling technique that combines the faster sampling in coarse-grained simulations with the potentially higher accuracy of all-atom simulations. After testing the efficiency of our Resolution Exchange with Tunneling (ResET) in simulations of the Trp-cage protein, an often used model to evaluate sampling techniques in protein simulations, we use our approach to compare the landscape of wild-type and A2T mutant Aβ_{1-42} peptides. Our results suggest a mechanism by that the mutation of a small hydrophobic alanine (A) into a bulky polar threonine (T) may interfere with the self-assembly of Aβ fibrils.

Anatomy and Formation Mechanisms of Early Amyloid-β Oligomers with Lateral Branching: Graph Network Analysis on Large-Scale Simulations

Oligomeric amyloid-β aggregates (AβOs) effectively trigger Alzheimer's disease-related toxicity, generating great interest in understanding their structures and formation mechanisms. However, AβOs are heterogeneous and transient, making their structure and formation difficult to study. Here, we performed graph network analysis of tens of microsecond massive simulations of early amyloid-β (Aβ) aggregations at near-atomic resolution to characterize AβO structures with sizes up to 20-mers. We found that AβOs exhibit highly curvilinear, irregular shapes with occasional lateral branches, consistent with recent cryo-electron tomography experiments. We also found that Aβ40 oligomers were more likely to develop branches than Aβ42 oligomers, explaining an experimental observation that only Aβ40 was trapped in network-like aggregates and exhibited slower fibrillization kinetics. Moreover, AβO architecture dissection revealed that their curvilinear appearance is related to the local packing geometries of neighboring peptides and that Aβ40's greater branching ability originates from specific C-terminal interactions at branching interfaces. Finally, we demonstrate that whether Aβ oligomerization causes oligomers to elongate or to branch depends on the sizes and shapes of colliding aggregates. Collectively, this study provides bottom-up structural information for understanding early Aβ aggregation and AβO toxicity.

Graph network analysis on large-scale simulations uncovers the differential branching behaviours of large Aβ40 and Aβ42 oligomers.

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.

Alzheimer's protection effect of A673T mutation may be driven by lower Aβ oligomer binding affinity

Several mutations conferring protection against Alzheimer's disease (AD) have been described, none as profound as the A673T mutation, where carriers are four times less likely to get AD compared to noncarriers. This mutation results in reduced amyloid beta (Aβ) protein production in vitro and lower lifetime Aβ concentration in carriers. Better understanding of the protective mechanisms of the mutation may provide important insights into AD pathophysiology and identify productive therapeutic intervention strategies for disease modification. Aβ(1-42) protein forms oligomers that bind saturably to a single receptor site on neuronal synapses, initiating the downstream toxicities observed in AD. Decreased formation, toxicity, or stability of soluble Aβ oligomers, or reduction of synaptic binding of these oligomers, may combine with overall lower Aβ concentration to underlie A673T's disease protecting mechanism. To investigate these possibilities, we compared the formation rate of soluble oligomers made from Icelandic A673T mutant and wild type (wt) Aβ(1-42) synthetic protein, the amount and intensity of oligomer bound to mature primary rat hippocampal/cortical neuronal synapses, and the potency of bound oligomers to impact trafficking rate in neurons in vitro using a physiologically relevant oligomer preparation method. At equal protein concentrations, mutant protein forms approximately 50% or fewer oligomers of high molecular weight (>50 kDa) compared to wt protein. Mutant oligomers are twice as potent at altering the cellular vesicle trafficking rate as wt at equivalent concentrations, however, mutant oligomers have a >4-fold lower binding affinity to synaptic receptors (Kd  = 1,950 vs. 442 nM). The net effect of these differences is a lower overall toxicity at a given concentration. This study demonstrates for the first time that mutant A673T Aβ oligomers prepared with this method have fundamentally different assembly characteristics and biological impact from wt protein and indicates that its disease protecting mechanism may result primarily from the mutant protein's much lower binding affinity to synaptic receptors. This suggests that therapeutics that effectively reduce oligomer binding to synapses in the brain may be beneficial in AD.

The role of amyloid oligomers in neurodegenerative pathologies

Many neurodegenerative diseases are rooted in the activities of amyloid-like proteins which possess conformations that spread to healthy proteins. These include Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD) and amyotrophic lateral sclerosis (ALS). While their clinical manifestations vary, their protein-level mechanisms are remarkably similar. Aberrant monomeric proteins undergo conformational shifts, facilitating aggregation and formation of solid fibrils. However, there is growing evidence that intermediate oligomeric stages are key drivers of neuronal toxicity. Analysis of protein dynamics is complicated by the fact that nucleation and growth of amyloid-like proteins is not a linear pathway. Feedback within this pathway results in exponential acceleration of aggregation, but activities exerted by oligomers and fibrils can alter cellular interactions and the cellular environment as a whole. The resulting cascade of effects likely contributes to the late onset and accelerating progression of amyloid-like protein disorders and the widespread effects they have on the body. In this review we explore the amyloid-like proteins associated with AD, PD, HD and ALS, as well as the common mechanisms of amyloid-like protein nucleation and aggregation. From this, we identify core elements of pathological progression which have been targeted for therapies, and which may become future therapeutic targets.

Effects of Ions and Small Compounds on the Structure of Aβ42 Monomers

Aggregation of amyloid-β (Aβ) proteins in the brain is a hallmark of Alzheimer's disease. This phenomenon can be promoted or inhibited by adding small molecules to the solution where Aβ is embedded. These molecules affect the ensemble of conformations sampled by Aβ monomers even before aggregation starts. Here, we perform extensive all-atom replica exchange molecular dynamics (REMD) simulations to provide a comparative study of the ensemble of conformations sampled by Aβ₄₂ monomers in solutions that promote (i.e., aqueous solution containing NaCl) and inhibit (i.e., aqueous solutions containing scyllo-inositol or 4-aminophenol) aggregation. Simulations performed in pure water are used as our reference. We find that secondary-structure content is only affected in an antagonistic manner by promoters and inhibitors at the C-terminus and the central hydrophilic core. Moreover, the end of the C-terminus binds more favorably to the central hydrophobic core region of Aβ₄₂ in NaCl adopting a type of strand-loop-strand structure that is disfavored by inhibitors. Nonpolar residues that form the dry core of larger aggregates of Aβ₄₂ (e.g., PDB ID 2BEG) are found at close proximity in these strand-loop-strand structures, suggesting that their formation could play an important role in initiating nucleation. In the presence of inhibitors, the C-terminus binds the central hydrophilic core with a higher probability than in our reference simulation. This sensitivity of the C-terminus, which is affected in an antagonistic manner by inhibitors and promoters, provides evidence for its critical role in accounting for aggregation.

Acetylation of Aβ40 Alters Aggregation in the Presence and Absence of Lipid Membranes

A hallmark of Alzheimer's disease (AD) is the formation of senile plaques comprised of the β-amyloid (Aβ) peptide. Aβ fibrillization is a complex nucleation-dependent process involving a variety of metastable intermediate aggregates and features the formation of inter- and intramolecular salt bridges involving lysine residues, K16 and K28. Cationic lysine residues also mediate protein-lipid interactions via association with anionic lipid headgroups. As several toxic mechanisms attributed to Aβ involve membrane interactions, the impact of acetylation on Aβ₄₀ aggregation in the presence and absence of membranes was determined. Using chemical acetylation, varying mixtures of acetylated and nonacetylated Aβ₄₀ were produced. With increasing acetylation, fibril and oligomer formation decreased, eventually completely arresting fibrillization. In the presence of total brain lipid extract (TBLE) vesicles, acetylation reduced the interaction of Aβ₄₀ with membranes; however, fibrils still formed at near complete levels of acetylation. Additionally, the combination of TBLE and acetylated Aβ promoted annular aggregates. Finally, toxicity associated with Aβ₄₀ was reduced with increasing acetylation in a cell culture assay. These results suggest that in the absence of membranes that the cationic character of lysine plays a major role in fibril formation. However, acetylation promotes unique aggregation pathways in the presence of lipid membranes.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

Oleuropein aglycone and hydroxytyrosol interfere differently with toxic Aβ1-42 aggregation

Oleuropein aglycone (OleA), the most abundant polyphenol in extra virgin olive oil (EVOO), and Hydroxythyrosol (HT), the OleA main metabolite, have attracted our interest due to their multitarget effects, including the interference with amyloid aggregation path. However, the mechanistic details of their anti-amyloid effect are not known yet. We report here a broad biophysical approach and cell biology techniques that enabled us to characterize the different molecular mechanisms by which OleA and HT modulate the Aβ_1-42 fibrillation, a main histopathological feature of Alzheimer's disease (AD). In particular, OleA prevents the growth of toxic Aβ_1-42 oligomers and blocks their successive growth into mature fibrils following its interaction with the peptide N-terminus, while HT speeds up harmless fibril formation. Our data demonstrate that, by stabilizing oligomers and fibrils, both polyphenols reduce their seeding activity and aggregate/membrane interaction on human neuroblastoma SH-SY5Y cells. These findings highlight the great potential of EVOO polyphenols and offer the possibility to validate and to optimize their use for possible AD prevention and therapy.

Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.

Significantly different contact patterns between Aβ40 and Aβ42 monomers involving the N-terminal region

Aβ42 peptide, with two additional residues at C-terminus, aggregates much faster than Aβ40. We performed equilibrium replica-exchange molecular dynamics simulations of their monomers using our residue-specific force field. Simulated ³ J_{HNH α} -coupling constants agree excellently with experimental data. Aβ40 and Aβ42 have very similar local conformational features, with considerable β-strand structures in the segments: A2-H6 (A), L17-A21 (B), A30-V36 (C) of both peptides and V39-I41 (D) of Aβ42. Both peptides have abundant A-B and B-C contacts, but Aβ40 has much more contacts between A and C than Aβ42, which may retard its aggregation. Only Aβ42 has considerable A-B-C-D topology. Decreased probability of A-C contact in Aβ42 relates to the competition from C-D contact. Increased A-C contact probability may also explain the slower aggregation of A2T and A2V mutants of Aβ42.

Cu and Zn coordination to amyloid peptides: From fascinating chemistry to debated pathological relevance

Several diseases share misfolding of different peptides and proteins as a key feature for their development. This is the case of important neurodegenerative diseases such as Alzheimer's and Parkinson's diseases and type II diabetes mellitus. Even more, metal ions such as copper and zinc might play an important role upon interaction with amyloidogenic peptides and proteins, which could impact their aggregation and toxicity abilities. In this review, the different coordination modes proposed for copper and zinc with amyloid-β, α-synuclein and IAPP will be reviewed as well as their impact on the aggregation, and ROS production in the case of copper. In addition, a special focus will be given to the mutations that affect metal binding and lead to familial cases of the diseases. Different modifications of the peptides that have been observed in vivo and could be relevant for the coordination of metal ions are also described.

Weaker N-Terminal Interactions for the Protective over the Causative Aβ Peptide Dimer Mutants

Knowing that abeta amyloid peptide (Aβ₄₂) dimers are the smallest and most abundant neurotoxic oligomers for Alzheimer's disease (AD), we used molecular simulations with advanced sampling methods (replica-exchange) to characterize and compare interactions between the N-termini (residues 1-16) of wild type (WT-WT) and five mutant dimers under constrained and unconstrained conditions. The number of contacts and distances between the N-termini, and contact maps of their conformational landscape illustrate substantial differences for a single residue change. The N-terminal contacts are significantly diminished for the dimers containing the monomers that protect against (WT-A2T) as compared with those that predispose toward (A2V-A2V) AD and for the control WT-WT dimers. The reduced number of N-terminal contacts not only occurs at or near the second residue mutations but also is distributed through to the 10th residue. These findings provide added support to the accumulating evidence for the "N-terminal hypothesis of AD" and offer an alternate mechanism for the cause of protection from the A2T mutant.

Looking at the Disordered Proteins through the Computational Microscope

Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure-function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

Emergence of Alternative Structures in Amyloid Beta 1-42 Monomeric Landscape by N-terminal Hexapeptide Amyloid Inhibitors

Alzheimer’s disease (AD) is characterized by deposition of amyloid beta (Aβ) peptides into senile plaques in the brain. While most familial mutations are associated with early-onset AD, recent studies report the AD-protective nature of two genetic human Aβ variants, i.e. A2T and A2V, in the heterozygous state. The mixture of A2V Aβ1-6 (Aβ₆) hexapeptide and WT Aβ1–42 (Αβ₄₂) is also found neuroprotective. Motivated by these findings, in this study we investigate the effects of WT, A2V, and A2T Aβ₆ hexapeptide binding on the monomeric WT Aβ₄₂ landscape. For this purpose, we have performed extensive atomistic Replica Exchange Molecular Dynamics simulations, elucidating preferential binding of Aβ₄₂ with the A2V and A2T hexapeptides compared to WT Aβ₆. A notable reorganization of the Aβ₄₂ landscape is revealed due to hexapeptide association, as manifested by lowering of transient interactions between the central and C-terminal hydrophobic patches. Concurrently, Aβ₆-bound Aβ₄₂ monomer exhibits alternative structural features that are strongly dependent on the hexapeptide sequence. For example, a central helix is more frequently populated within the A2T-bound monomer, while A2V-bound Aβ₄₂ is often enhanced in overall disorder. Taken together, the present simulations offer novel molecular insights onto the effect of the N-terminal hexapeptide binding on the Aβ₄₂ monomer structure, which might help in explaining their reported amyloid inhibition properties.

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.

Pathogenic Aβ A2V versus protective Aβ A2T mutation: Early stage aggregation and membrane interaction

We investigated the effects of punctual A-to-V and A-to-T mutations in the amyloid precursor protein APP, corresponding to position 2 of Aβ1-42. Those mutations had opposite effects on the onset and progression of Alzheimer disease, the former inducing early AD pathology and the latter protecting against the onset of the disease. We applied Static and Dynamic Light Scattering and Circular Dichroism, to study the different mutants in the early stages of the aggregation process, essential for the disease. Comparative results showed that the aggregation pathways differ in the kinetics and extent of the process, in the size of the aggregates and in the evolution of the secondary structure, resulting in fibrils of different morphology, as seen by AFM. Mutated peptides had comparable toxic effects on N2a cells. Moreover, as assessed by X-ray scattering, all of them displayed disordering effects on the internal structure of mixed phospholipids-gangliosides model membranes.

Pyroglutamate-Modified Amyloid-β(3-42) Shows α-Helical Intermediates before Amyloid Formation

Pyroglutamate-modified amyloid-β (pEAβ) has been described as a relevant Aβ species in Alzheimer's-disease-affected brains, with pEAβ (3-42) as a dominant isoform. Aβ (1-40) and Aβ (1-42) have been well characterized under various solution conditions, including aqueous solutions containing trifluoroethanol (TFE). To characterize structural properties of pEAβ (3-42) possibly underlying its drastically increased aggregation propensity compared to Aβ (1-42), we started our studies in various TFE-water mixtures and found striking differences between the two Aβ species. Soluble pEAβ (3-42) has an increased tendency to form β-sheet-rich structures compared to Aβ (1-42), as indicated by circular dichroism spectroscopy data. Kinetic assays monitored by thioflavin-T show drastically accelerated aggregation leading to large fibrils visualized by electron microscopy of pEAβ (3-42) in contrast to Aβ (1-42). NMR spectroscopy was performed for backbone and side-chain chemical-shift assignments of monomeric pEAβ (3-42) in 40% TFE solution. Although the difference between pEAβ (3-42) and Aβ (1-42) is purely N-terminal, it has a significant impact on the chemical environment of >20% of the total amino acid residues, as revealed by their NMR chemical-shift differences. Freshly dissolved pEAβ (3-42) contains two α-helical regions connected by a flexible linker, whereas the N-terminus remains unstructured. We found that these α-helices act as a transient intermediate to β-sheet and fibril formation of pEAβ (3-42).

Alzheimer's amyloid-β A2T variant and its N-terminal peptides inhibit amyloid-β fibrillization and rescue the induced cytotoxicity

Alzheimer’s disease (AD) is the most common dementia affecting tens of million people worldwide. The primary neuropathological hallmark in AD is amyloid plaques composed of amyloid-β peptide (Aβ). Several familial mutations found in Aβ sequence result in early onset of AD. Previous studies showed that the mutations located at N-terminus of Aβ, such as the English (H6R) and Tottori (D7N) mutations, promote fibril formation and increase cytotoxicity. However, A2T mutant located at the very N-terminus of Aβ shows low-prevalence incidence of AD, whereas, another mutant A2V causes early onset of AD. To understand the molecular mechanism of the distinct effect and develop new potential therapeutic strategy, here, we examined the effect of full-length and N-terminal A2V/T variants to wild type (WT) Aβ40 by fibrillization assays and NMR studies. We found that full-length and N-terminal A2V accelerated WT fibrillization and induced large chemical shifts on the N-terminus of WT Aβ, whereas, full-length and N-terminal A2T retarded the fibrillization. We further examined the inhibition effect of various N-terminal fragments (NTFs) of A2T to WT Aβ. The A2T NTFs ranging from residue 1 to residue 7 to 10, but not 1 to 6 or shorter, are capable to retard WT Aβ fibrillization and rescue cytotoxicity. The results suggest that in the presence of full-length or specific N-terminal A2T can retard Aβ aggregation and the A2T NTFs can mitigate its toxicity. Our results provide a novel targeting site for future therapeutic development of AD.

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.

Alzheimer's Protective Cross-Interaction between Wild-Type and A2T Variants Alters Aβ42 Dimer Structure

Whole genome sequencing has recently revealed the protective effect of a single A2T mutation in heterozygous carriers against Alzheimer's disease (AD) and age-related cognitive decline. The impact of the protective cross-interaction between the wild-type (WT) and A2T variants on the dimer structure is therefore of high interest, as the Aβ dimers are the smallest known neurotoxic species. Toward this goal, extensive atomistic replica exchange molecular dynamics simulations of the solvated WT homo- and A2T hetero- Aβ_1-42 dimers have been performed, resulting into a total of 51 μs of sampling for each system. Weakening of a set of transient, intrachain contacts formed between the central and C-terminal hydrophobic residues is observed in the heterodimeric system. The majority of the heterodimers with reduced interaction between central and C-terminal regions lack any significant secondary structure and display a weak interchain interface. Interestingly, the A2T N-terminus, particularly residue F4, is frequently engaged in tertiary and quaternary interactions with central and C-terminal hydrophobic residues in those distinct structures, leading to hydrophobic burial. This atypical involvement of the N-terminus within A2T heterodimer revealed in our simulations implies possible interference on Aβ₄₂ aggregation and toxic oligomer formation, which is consistent with experiments. In conclusion, the present study provides detailed structural insights onto A2T Aβ₄₂ heterodimer, which might provide molecular insights onto the AD protective effect of the A2T mutation in the heterozygous state.

N-Terminal Hypothesis for Alzheimer's Disease

Although the amyloid (abeta peptide, Aβ) hypothesis is 25 years old, is the dominant model of Alzheimer's disease (AD) pathogenesis, and guides the development of potential treatments, it is still controversial. One possible reason is a lack of a mechanistic path from the cleavage products of the amyloid precursor protein (APP) such as soluble Aβ monomer and soluble molecular fragments to the deleterious effects on synaptic form and function. From a review of the recent literature and our own published work including aggregation kinetics and structural morphology, Aβ clearance, molecular simulations, long-term potentiation measurements with inhibition binding, and the binding of a commercial monoclonal antibody, aducanumab, we hypothesize that the N-terminal domains of neurotoxic Aβ oligomers are implicated in causing the disease.

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

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.

Exploring the aggregation free energy landscape of the amyloid-β protein (1-40)

A predictive coarse-grained protein force field [associative memory, water-mediated, structure, and energy model for molecular dynamics (AWSEM)-MD] is used to study the energy landscapes and relative stabilities of amyloid-β protein (1-40) in the monomer and all of its oligomeric forms up to an octamer. We find that an isolated monomer is mainly disordered with a short α-helix formed at the central hydrophobic core region (L17-D23). A less stable hairpin structure, however, becomes increasingly more stable in oligomers, where hydrogen bonds can form between neighboring monomers. We explore the structure and stability of both prefibrillar oligomers that consist of mainly antiparallel β-sheets and fibrillar oligomers with only parallel β-sheets. Prefibrillar oligomers are polymorphic but typically take on a cylindrin-like shape composed of mostly antiparallel β-strands. At the concentration of the simulation, the aggregation free energy landscape is nearly downhill. We use umbrella sampling along a structural progress coordinate for interconversion between prefibrillar and fibrillar forms to identify a conversion pathway between these forms. The fibrillar oligomer only becomes favored over its prefibrillar counterpart in the pentamer where an interconversion bottleneck appears. The structural characterization of the pathway along with statistical mechanical perturbation theory allow us to evaluate the effects of concentration on the free energy landscape of aggregation as well as the effects of the Dutch and Arctic mutations associated with early onset of Alzheimer's disease.

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.

A2T and A2V Aβ peptides exhibit different aggregation kinetics, primary nucleation, morphology, structure, and LTP inhibition

The histopathological hallmark of Alzheimer's disease (AD) is the aggregation and accumulation of the amyloid beta peptide (Aβ) into misfolded oligomers and fibrils. Here we examine the biophysical properties of a protective Aβ variant against AD, A2T, and a causative mutation, A2T, along with the wild type (WT) peptide. The main finding here is that the A2V native monomer is more stable than both A2T and WT, and this manifests itself in different biophysical behaviors: the kinetics of aggregation, the initial monomer conversion to an aggregation prone state (primary nucleation), the abundances of oligomers, and extended conformations. Aggregation reaction modeling of the conversion kinetics from native monomers to fibrils predicts the enhanced stability of the A2V monomer, while ion mobility spectrometry-mass spectrometry measures this directly confirming earlier predictions. Additionally, unique morphologies of the A2T aggregates are observed using atomic force microscopy, providing a basis for the reduction in long term potentiation inhibition of hippocampal cells for A2T compared with A2V and the wild type (WT) peptide. The stability difference of the A2V monomer and the difference in aggregate morphology for A2T (both compared with WT) are offered as alternate explanations for their pathological effects.

Amylin-Aβ oligomers at atomic resolution using molecular dynamics simulations: a link between Type 2 diabetes and Alzheimer's disease

Clinical studies have identified Type 2 diabetes (T2D) as a risk factor of Alzheimer's disease (AD). One of the potential mechanisms that link T2D and AD is the loss of cells associated with degenerative changes. Amylin1-37 aggregates (the pathological species in T2D) were found to be co-localized with those of Aβ1-42 (the pathological species in AD) to form the Amylin1-37-Aβ1-42 plaques, promoting aggregation and thus contributing to the etiology of AD. However, the mechanisms by which Amylin1-37 co-aggregates with Aβ1-42 are still elusive. This work presents the interactions between Amylin1-37 oligomers and Aβ1-42 oligomers at atomic resolution applying extensive molecular dynamics simulations for relatively large ensemble of cross-seeding Amylin1-37-Aβ1-42 oligomers. The main conclusions of this study are first, Aβ1-42 oligomers prefer to interact with Amylin1-37 oligomers to form single layer conformations (in-register interactions) rather than double layer conformations; and second, in some double layer conformations of the cross-seeding Amylin1-37-Aβ1-42 oligomers, the Amylin1-37 oligomers destabilize the Aβ1-42 oligomers and thus inhibit Aβ1-42 aggregation, while in other double layer conformations, the Amylin1-37 oligomers stabilize Aβ1-42 oligomers and thus promote Aβ1-42 aggregation.

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.

Characterization of Aβ Monomers through the Convergence of Ensemble Properties among Simulations with Multiple Force Fields

Amyloid β (Aβ) monomers represent a base state in the pathways of aggregation that result in the fibrils and oligomers implicated in the pathogenesis of Alzheimer's disease (AD). The structural properties of these intrinsically disordered peptides remain unclear despite extensive experimental and computational investigations. Further, there are mutations within Aβ that change the way the peptide aggregates and are known to cause familial AD (FAD). Here, we analyze the ensembles of different isoforms (Aβ42 and Aβ40) and mutants (E22Δ, D23N, E22K, E22G, and A2T in Aβ40) of Aβ generated with all-atom replica exchange molecular dynamics (REMD) simulations on the μs/replica time scale. These were run using three different force field/water model combinations: OPLS-AA/L and TIP3P ("OPLS"), AMBER99sb-ILDN and TIP4P-Ew ("ILDN"), as well as CHARMM22* and TIP3SP ("CHARMM"). Despite fundamental changes in simulation parameters, we find that the resulting ensembles demonstrate a strong convergence in structural properties. In particular, antiparallel contacts between L17-A21 and A30-L34 are prevalent in ensembles of Aβ40, directly forming β sheets in the OPLS and ILDN combinations. A21-A30 commonly forms an interceding region that rarely interacts with the rest of the peptide. Further, Aβ42 contributes new β hairpin motifs involving V40-I41 in both OPLS and ILDN. However, the structural flexibility of the central region and the electrostatic interactions that characterize it are notably different between the different conditions. Further, for OPLS, each of the FAD mutations disrupts central bend character and increases the polymorphism of antiparallel contacts across the central region. However, the studied mutations in the ILDN set primarily encourage more global contacts involving the N-terminus and the central region, and promote the formation of new β topologies that may seed different aggregates involved in disease phenotypes. These differences aside, the large degree of agreement between simulation sets across multiple force fields provides a generalizable characterization of Aβ that is also consistent with experimental data and models.

Amyloid β-Protein Assembly: Differential Effects of the Protective A2T Mutation and Recessive A2V Familial Alzheimer's Disease Mutation

Oligomeric states of the amyloid β-protein (Aβ) appear to be causally related to Alzheimer's disease (AD). Recently, two familial mutations in the amyloid precursor protein gene have been described, both resulting in amino acid substitutions at Ala2 (A2) within Aβ. An A2V mutation causes autosomal recessive early onset AD. Interestingly, heterozygotes enjoy some protection against development of the disease. An A2T substitution protects against AD and age-related cognitive decline in non-AD patients. Here, we use ion mobility-mass spectrometry (IM-MS) to examine the effects of these mutations on Aβ assembly. These studies reveal different assembly pathways for early oligomer formation for each peptide. A2T Aβ42 formed dimers, tetramers, and hexamers, but dodecamer formation was inhibited. In contrast, no significant effects on Aβ40 assembly were observed. A2V Aβ42 also formed dimers, tetramers, and hexamers, but it did not form dodecamers. However, A2V Aβ42 formed trimers, unlike A2T or wild-type (wt) Aβ42. In addition, the A2V substitution caused Aβ40 to oligomerize similar to that of wt Aβ42, as evidenced by the formation of dimers, tetramers, hexamers, and dodecamers. In contrast, wt Aβ40 formed only dimers and tetramers. These results provide a basis for understanding how these two mutations lead to, or protect against, AD. They also suggest that the Aβ N-terminus, in addition to the oft discussed central hydrophobic cluster and C-terminus, can play a key role in controlling disease susceptibility.

Oligomer Formation of Toxic and Functional Amyloid Peptides Studied with Atomistic Simulations

Amyloids are associated with diseases, including Alzheimer's, as well as functional roles such as storage of peptide hormones. It is still unclear what differences exist between aberrant and functional amyloids. However, it is known that soluble oligomers formed during amyloid aggregation are more toxic than the final fibrils. Here, we perform molecular dynamics simulations to study the aggregation of the amyloid-β peptide Aβ25-35, associated with Alzheimer's disease, and two functional amyloid-forming tachykinin peptides: kassinin and neuromedin K. Although the three peptides have similar primary sequences, tachykinin peptides, in contrast to Aβ25-35, form nontoxic amyloids. Our simulations reveal that the charge of the C-terminus is essential to controlling the aggregation process. In particular, when the kassinin C-terminus is not amidated, the aggregation kinetics decreases considerably. In addition, we observe that the monomeric peptides in extended conformations aggregate faster than those in collapsed hairpin-like conformations.

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.