Dual effects of familial Alzheimer's disease mutations (D7H, D7N, and H6R) on amyloid β peptide: Correlation dynamics and zinc binding

An N-terminal acidic β-sheet domain is responsible for the metal-accumulation properties of amyloid-β protofibrils: a molecular dynamics study

The influence of metal ions on the structure of amyloid- β (Aβ) protofibril models was studied through molecular dynamics to explore the molecular mechanisms underlying metal-induced Aβ aggregation relevant in Alzheimer's disease (AD). The models included 36-, 48-, and 188-mers of the Aβ42 sequence and two disease-modifying variants. Primary structural effects were observed at the N-terminal domain, as it became susceptible to the presence of cations. Specially when β-sheets predominate, this motif orients N-terminal acidic residues toward one single face of the β-sheet, resulting in the formation of an acidic region that attracts cations from the media and promotes the folding of the N-terminal region, with implications in amyloid aggregation. The molecular phenotype of the protofibril models based on Aβ variants shows that the AD-causative D7N mutation promotes the formation of N-terminal β-sheets and accumulates more Zn2+, in contrast to the non-amyloidogenic rodent sequence that hinders the β-sheets and is more selective for Na+ over Zn2+ cations. It is proposed that forming an acidic β-sheet domain and accumulating cations is a plausible molecular mechanism connecting the elevated affinity and concentration of metals in Aβ fibrils to their high content of β-sheet structure at the N-terminal sequence.

Structure Comparison of Beta Amyloid Peptide Aβ 1-42 Isoforms. Molecular Dynamics Modeling

Beta amyloid peptide Aβ 1-42 (Aβ42) has a unique dual role in the human organism, as both the peptide with an important physiological function and one of the most toxic biological compounds provoking Alzheimer's disease (AD). There are several known Aβ42 isoforms that we discuss here that are highly neurotoxic and lead to the early onset of AD. Aβ42 is an intrinsically disordered protein with no experimentally solved structure under physiological conditions. The objective of this research was to establish the appropriate molecular dynamics (MD) methodology and model a uniform set of structures for the Aβ42 isoforms that form the core of this study. For that purpose, force field selection and verification including convergence testing for MD simulations was made. Replica exchange MD and conventional MD modeling of several Aβ42 and Aβ16 isoforms that have neurotoxic and amyloidogenic effects impacting the severity of Alzheimer's disease were carried out with the optimal force field and solvent parameters. A standardized ensemble of structures for the Aβ42 and Aβ16 isoforms covering 30-50% of the conformational ensembles extracted from the free energy minima was calculated from MD trajectories. The resulting data set of modeled structures includes Aβ42 wild type, isoD7, pS8, D7H, and H6R-Aβ42 and Aβ16 wild type, isoD7, pS8, D7H, and H6R-Aβ16. The representative structures are given in the Supporting Information; they are open for public access. In the study, we also evaluated the differences between the structures of Aβ42 isoforms and speculate on their possible relevance to the known functions. Utilizing several representative structures for a single disordered protein for docking, with their subsequent averaging by conformations, would markedly increase the reliability of docking results.

Molecular Mechanism of Zinc-Dependent Oligomerization of Alzheimer's Amyloid-β with Taiwan (D7H) Mutation

Amyloid-β (Aβ) is a peptide formed by 39-43 amino acids, heterogenous by the length of its C-terminus. Aβ constitutes a subnanomolar monomeric component of human biological fluids; however, in sporadic variants of Alzheimer's disease (AD), it forms soluble neurotoxic oligomers and accumulates as insoluble extracellular polymeric aggregates (amyloid plaques) in the brain tissues. The plaque formation is controlled by zinc ions; therefore, abnormal interactions between the ions and Aβ seem to take part in the triggering of sporadic AD. The amyloid plaques contain various Aβ isoforms, among which the most common is Aβ with an isoaspartate in position 7 (isoD7). The spontaneous conversion of D7 to isoD7 is associated with Aβ aging. Aβ molecules with isoD7 (isoD7-Aβ) easily undergo zinc-dependent oligomerization, and upon administration to transgenic animals (mice, nematodes) used for AD modeling, act as zinc-dependent seeds of the pathological aggregation of Aβ. The formation of zinc-bound homo- and hetero-oligomers with the participation of isoD7-Aβ is based on the rigidly structured segment 11-EVHH-14, located in the Aβ metal binding domain (Aβ₁₆). Some hereditary variants of AD are associated with familial mutations within the domain. Among these, the most susceptible to zinc-dependent oligomerization is Aβ with Taiwan (D7H) mutation (D7H-Aβ). In this study, the D7H-Aβ metal binding domain (D7H-Aβ₁₆) has been used as a model to establish the molecular mechanism of zinc-induced D7H-Aβ oligomerization through turbidimetry, dynamic light scattering, isothermal titration calorimetry, mass spectrometry, and computer modelling. Additionally, the modeling data showed that a molecule of D7H-Aβ, as well as isoD7-Aβ in combination with two Aβ molecules, renders a stable zinc-induced heterotrimer. The trimers are held together by intermolecular interfaces via zinc ions, with the primary interfaces formed by 11-EVHH-14 sites of the interacting trimer subunits. In summary, the obtained results confirm the role of the 11-EVHH-14 region as a structure and function determinant for the zinc-dependent oligomerization of all known Aβ species (including various chemically modified isoforms and AD-associated mutants) and point at this region as a potent target for drugs aimed to stop amyloid plaque formation in both sporadic and hereditary variants of AD.

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.

The N-Terminal Domain of Aβ40-Amyloid Fibril: The MOMD Perspective of its Dynamic Structure from NMR Lineshape Analysis

We have developed the stochastic microscopic-order-macroscopic-disorder (MOMD) approach for elucidating dynamic structures in the solid-state from ²H NMR lineshapes. In MOMD, the probe experiences an effective/collective motional mode. The latter is described by a potential, u, which represents the local spatial-restrictions, a local-motional diffusion tensor, R, and key features of local geometry. Previously we applied MOMD to the well-structured core domain of the 3-fold-symmetric twisted polymorph of the Aβ₄₀-amyloid fibril. Here, we apply it to the N-terminal domain of this fibril. We find that the dynamic structures of the two domains are largely similar but differ in the magnitude and complexity of the key physical parameters. This interpretation differs from previous multisimple-mode (MSM) interpretations of the same experimental data. MSM used for the two domains different combinations of simple motional modes taken to be independent. For the core domain, MOMD and MSM disagree on the character of the dynamic structure. For the N-terminal domain, they even disagree on whether this chain segment is structurally ordered (MOMD finds that it is), and whether it undergoes a phase transition at 260 K where bulklike water located in the fibril matrix freezes (MOMD finds that it does not). These are major differences associated with an important system. While the MOMD description is a physically sound one, there are drawbacks in the MSM descriptions. The results obtained in this study promote our understanding of the dynamic structure of protein aggregates. Thus, they contribute to the effort to pharmacologically control neurodegenerative disorders believed to be caused by such aggregates.

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.

The cellular modifier MOAG-4/SERF drives amyloid formation through charge complementation

While aggregation-prone proteins are known to accelerate aging and cause age-related diseases, the cellular mechanisms that drive their cytotoxicity remain unresolved. The orthologous proteins MOAG-4, SERF1A, and SERF2 have recently been identified as cellular modifiers of such proteotoxicity. Using a peptide array screening approach on human amyloidogenic proteins, we found that SERF2 interacted with protein segments enriched in negatively charged and hydrophobic, aromatic amino acids. The absence of such segments, or the neutralization of the positive charge in SERF2, prevented these interactions and abolished the amyloid-promoting activity of SERF2. In protein aggregation models in the nematode worm Caenorhabditis elegans, protein aggregation and toxicity were suppressed by mutating the endogenous locus of MOAG-4 to neutralize charge. Our data indicate that MOAG-4 and SERF2 drive protein aggregation and toxicity by interactions with negatively charged segments in aggregation-prone proteins. Such charge interactions might accelerate primary nucleation of amyloid by initiating structural changes and by decreasing colloidal stability. Our study points at charge interactions between cellular modifiers and amyloidogenic proteins as potential targets for interventions to reduce age-related protein toxicity.

The English (H6R) Mutation of the Alzheimer's Disease Amyloid-β Peptide Modulates Its Zinc-Induced Aggregation

The coordination of zinc ions by histidine residues of amyloid-beta peptide (Aβ) plays a critical role in the zinc-induced Aβ aggregation implicated in Alzheimer's disease (AD) pathogenesis. The histidine to arginine substitution at position 6 of the Aβ sequence (H6R, English mutation) leads to an early onset of AD. Herein, we studied the effects of zinc ions on the aggregation of the Aβ42 peptide and its isoform carrying the H6R mutation (H6R-Aβ42) by circular dichroism spectroscopy, dynamic light scattering, turbidimetric and sedimentation methods, and bis-ANS and thioflavin T fluorescence assays. Zinc ions triggered the occurrence of amorphous aggregates for both Aβ42 and H6R-Aβ42 peptides but with distinct optical properties. The structural difference of the formed Aβ42 and H6R-Aβ42 zinc-induced amorphous aggregates was also supported by the results of the bis-ANS assay. Moreover, while the Aβ42 peptide demonstrated an increase in the random coil and β-sheet content upon complexing with zinc ions, the H6R-Aβ42 peptide showed no appreciable structural changes under the same conditions. These observations were ascribed to the impact of H6R mutation on a mode of zinc/peptide binding. The presented findings further advance the understanding of the pathological role of the H6R mutation and the role of H6 residue in the zinc-induced Aβ aggregation.

Molecular dynamics simulations of copper binding to N-terminus mutants of amyloid-β

We report results of molecular dynamic (MD) simulations on N-terminus mutants of the copper-bound, amyloid-β (Aβ) peptide. Eight structures of Aβ were modelled, including seven mutant peptides in addition to the unaltered wild-type (WT). Trajectories analysed for each individual system were all approximately 1.4 μs in length, yielding a total of over 11 μs in total. The impact of these mutations are marked and varied compared to the wild-type peptide, including effects on secondary structure, stability and conformational changes. Each system showed differing levels of stability with some showing consistent, compact conformations whereas others displayed more flexible structures. Contrasts between comparable mutations at similar sites, such as A2T/A2V and D7H/D7N, show the location as well as the type of mutation have effects on protein structure observed in Ramachandran plots. We also report notable changes in peptide structure at residues remote to the site of substitution showing these mutations influence the entirety of Aβ. Salt-bridge profiles show this most clearly: addition or removal of charged residues affecting all salt-bridge interactions present in WT, even those remote from the site of mutation. Effects on secondary structure differ between mutations, most notably a change in incidence of β-strand, which has been linked to enhanced aggregational properties for the peptide. GFN2-xTB semi-empirical calculations show clear differences in binding energies of the copper-centre for each system.Communicated by Ramaswamy H. Sarma.

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.

Solid-state NMR reveals a comprehensive view of the dynamics of the flexible, disordered N-terminal domain of amyloid-β fibrils

Amyloid fibril deposits observed in Alzheimer's disease comprise amyloid-β (Aβ) protein possessing a structured hydrophobic core and a disordered N-terminal domain (residues 1-16). The internal flexibility of the disordered domain is likely essential for Aβ aggregation. Here, we used ²H static solid-state NMR methods to probe the dynamics of selected side chains of the N-terminal domain of Aβ_1-40 fibrils. Line shape and relaxation data suggested a two-state model in which the domain's free state undergoes a diffusive motion that is quenched in the bound state, likely because of transient interactions with the structured C-terminal domain. At 37 °C, we observed freezing of the dynamics progressively along the Aβ sequence, with the fraction of the bound state increasing and the rate of diffusion decreasing. We also found that without solvation, the diffusive motion is quenched. The solvent acted as a plasticizer reminiscent of its role in the onset of global dynamics in globular proteins. As the temperature was lowered, the fraction of the bound state exhibited sigmoidal behavior. The midpoint of the freezing curve coincided with the bulk solvent freezing for the N-terminal residues and increased further along the sequence. Using ²H R _1ρ measurements, we determined the conformational exchange rate constant between the free and bound states under physiological conditions. Zinc-induced aggregation leads to the enhancement of the dynamics, manifested by the faster conformational exchange, faster diffusion, and lower freezing-curve midpoints.

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.

Kinetics of the Interactions between Copper and Amyloid-β with FAD Mutations and Phosphorylation at the N terminus

Mutations and post‐translational modifications of amyloid‐β (Aβ) peptide in its N terminus have been shown to increase fibril formation, yet the molecular mechanism is not clear. Here we investigated the kinetics of the interactions of copper with two Aβ peptides containing Familial Alzheimer's disease (FAD) mutations (English (H6R) and Tottori (D7N)), as well as with Aβ peptide phosphorylated at serine 8 (pS8). All three peptides bind to copper with a similar rate as the wild‐type (wt). The dissociation rates follow the order pS8>H6R>wt>D7N; the interconversion between the two coordinating species occurs 50 % faster for H6R and pS8, whereas D7N had only a negligible effect. Interestingly, the rate of ternary complex (copper‐bridged heterodimer) formation for the modified peptides was significantly faster than that for wt, thus leading us to propose that FAD and sporadic AD might share a kinetic origin for the enhanced oligomerisation of Aβ.

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.

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

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

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.