Side chain interactions can impede amyloid fibril growth: replica exchange simulations of Abeta peptide mutant

Revealing the Nanoarchitectonics of Amyloid β-Aggregation on Two-Dimensional Biomimetic Membranes by Surface-Enhanced Infrared Absorption Spectroscopy

The in vivo folding of amyloid β (Aβ) is influenced by many factors among which biomembrane interfaces play an important role. Here, using surface-enhanced infrared absorption (SEIRA) spectroscopy and atomic force microscopy (AFM), the adsorption, structure, and morphology of Aβ42 aggregating on different two-dimensional interfaces were investigated. Results show that interfaces facilitate the aggregation of Aβ42 and are conducive to the formation of homogeneous aggregates, while the aggregates vary on different interfaces. On hydrophobic interfaces, strong hydrophobic interactions with the C-terminus of Aβ42 result in the formation of small oligomers with a small proportion of the β-sheet structure. On hydrophilic interfaces, hydrogen-bonding interactions and electrostatic interactions promote the formation of large aggregate particles with β-sheet structure. The hydration repulsion plays an important role in the interaction of Aβ42 with interfaces. These findings help to understand the nature of Aβ42 adsorption and aggregation on the biomembrane interface and the origin of heterogeneity and polymorphism of Aβ42 aggregates.

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.

Oligomerization of FVFLM peptides and their ability to inhibit beta amyloid peptides aggregation: consideration as a possible model

Preeclampsia, a pregnancy-specific disorder, shares typical pathophysiological features with protein misfolding disorders including Alzheimer's disease. Characteristic for preeclampsia is the involvement of multiple proteins of which fragments of SERPINA1 and β-amyloid co-aggregate in urine and placenta of preeclamptic women. To explore the biophysical basis of this interaction, we investigated the multidimensional efficacy of the FVFLM sequence in SERPINA1, as a model inhibitory agent of β-amyloid aggregation. After studying the oligomerization of FVFLM peptides using all-atom molecular dynamics simulations with the GROMOS43a1 force field and explicit water, we report that FVFLM can aggregate and its aggregation is spontaneous with a remarkably faster rate than that recorded for KLVFF (aggregation "hot-spot" from β-amyloid). The fast kinetics of FVFLM aggregation was found to be driven primarily by core-like aromatic interactions originating from the anti-parallel orientation of complementarily uncharged strands. The conspicuously stable aggregation mechanism observed for FVFLM peptides is found not to conform to the popular 'dock-lock' scheme. We also found high propensity of FVFLM for KLVFF binding. When present, FVFLM disrupts the β-amyloid aggregation pathway and we propose that FVFLM-like peptides might be used to prevent the assembly of full-length Aβ or other pro-amyloidogenic peptides into amyloid fibrils.

Is the Conformational Ensemble of Alzheimer's Aβ10-40 Peptide Force Field Dependent?

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 ³J_HNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of ³J_HNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble.

Dependence of protein conformational ensembles on force field parameterizations limits the predictive power of molecular dynamics simulations. To address this problem, we evaluated five all-atom force fields for their consistency in reproducing the conformational ensemble of Alzheimer’s Aβ10-40 peptide. To generate conformational ensembles, we have used replica exchange molecular dynamics and computed Aβ10-40 secondary and tertiary structures. We found that, although all force fields predict Aβ10-40 unfolded structure, they strongly disagree on helix and β propensities and tertiary structure distributions. We have also calculated Aβ10-40 J-coupling and residual dipolar coupling constants and compared them with the experimental data for the full-length Aβ1-40 peptide. Unexpectedly, we determined that in silico Aβ10-40 and experimental Aβ1-40 constants are in better agreement than these quantities computed and measured previously for identical peptides, such as Aβ1-40 or Aβ1-42. We then concluded that the conformational ensembles of Aβ10-40 and Aβ1-40 are similar and on this basis argue that CHARMM36 force field with standard TIP3P water model provides the most accurate description of Aβ10-40. Although our objective was not to evaluate the biomolecular force fields in general, our study is expected to facilitate their proper selection for the simulations of Alzheimer’s peptides.

Chirality-assisted ring-like aggregation of aβ(1-40) at liquid-solid interfaces: a stereoselective two-step assembly process

Molecular chirality is introduced at liquid-solid interfaces. A ring-like aggregation of amyloid Aβ(1-40) on N-isobutyryl-L-cysteine (L-NIBC)-modified gold substrate occurs at low Aβ(1-40) concentration, while D-NIBC modification only results in rod-like aggregation. Utilizing atomic force microscope controlled tip-enhanced Raman scattering, we directly observe the secondary structure information for Aβ(1-40) assembly in situ at the nanoscale. D- or L-NIBC on the surface can guide parallel or nonparallel alignment of β-hairpins through a two-step process based on electrostatic-interaction-enhanced adsorption and subsequent stereoselective recognition. Possible electrostatic interaction sites (R5 and K16) and a chiral recognition site (H14) of Aβ(1-40) are proposed, which may provide insight into the understanding of this effect.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

Conformational stability of fibrillar amyloid-beta oligomers via protofilament pair formation - a systematic computational study

Amyloid-[Formula: see text] (A[Formula: see text]) oligomers play a crucial role in Alzheimer's disease due to their neurotoxic aggregation properties. Fibrillar A[Formula: see text] oligomerization can lead to protofilaments and protofilament pairs via oligomer elongation and oligomer association, respectively. Small fibrillar oligomers adopt the protofilament topology, whereas fibrils contain at least protofilament pairs. To date, the underlying growth mechanism from oligomers to the mature fibril still remains to be elucidated. Here, we performed all-atom molecular dynamics simulations in explicit solvent on single layer-like protofilaments and fibril-like protofilament pairs of different size ranging from the tetramer to the 48-mer. We found that the initial U-shaped topology per monomer is maintained over time in all oligomers. The observed deviations of protofilaments from the starting structure increase significantly with size due to the twisting of the in-register parallel [Formula: see text]-sheets. This twist causes long protofilaments to be unstable and leads to a breakage. Protofilament pairs, which are stabilized by a hydrophobic interface, exhibit more fibril-like properties such as the overall structure and the twist angle. Thus, they can act as stable conformational templates for further fibril growth. Key properties like the twist angle, shape complementarity, and energetics show a size-dependent behavior so that small oligomers favor the protofilament topology, whereas large oligomers favor the protofilament pair topology. The region for this conformational transition is at the size of approximately twelve A[Formula: see text] monomers. From that, we propose the following growth mechanism from A[Formula: see text] oligomers to fibrils: (1) elongation of short protofilaments; (2) breakage of large protofilaments; (3) formation of short protofilament pairs; and (4) elongation of protofilament pairs.

Template induced conformational change of amyloid-β monomer

Population of aggregation-prone conformers for the monomeric amyloid-β (Aβ) can dramatically speed up its fibrillar aggregation. In this work, we study the effect of preformed template on the conformational distributions of the monomeric Aβ by replica exchange molecular dynamics. Our results show that the template consisting of Aβ peptides with cross-β structure can induce the formation of β-rich conformations for the monomeric Aβ, which is the key feature of the aggregation-prone conformers. Similar effect is observed when the hIAPP peptides and poly alanine peptides were used as templates, suggesting that the template effect is insensitive to the sequence details of the template peptides. In comparison, the template with helical structure has no significant effects on the β-propensity of the monomeric Aβ. Analysis to the interaction details revealed that the template tends to disrupt the intrapeptide interactions of the monomeric Aβ, which are absent in the fibrillar state, suggesting that the preformed template can reorganize the intrapeptide interactions of the monomeric Aβ during the capturing stage and reduce the energy frustrations for the fibrillar aggregations.

Biochemistry of amyloid β-protein and amyloid deposits in Alzheimer disease

Progressive cerebral deposition of the amyloid β-protein (Aβ) in brain regions serving memory and cognition is an invariant and defining feature of Alzheimer disease. A highly similar but less robust process accompanies brain aging in many nondemented humans, lower primates, and some other mammals. The discovery of Aβ as the subunit of the amyloid fibrils in meningocerebral blood vessels and parenchymal plaques has led to innumerable studies of its biochemistry and potential cytotoxic properties. Here we will review the discovery of Aβ, numerous aspects of its complex biochemistry, and current attempts to understand how a range of Aβ assemblies, including soluble oligomers and insoluble fibrils, may precipitate and promote neuronal and glial alterations that underlie the development of dementia. Although the role of Aβ as a key molecular factor in the etiology of Alzheimer disease remains controversial, clinical trials of amyloid-lowering agents, reviewed elsewhere in this book, are poised to resolve the question of its pathogenic primacy.

Coarse-grained simulations of protein aggregation

Protein aggregation is believed to be responsible for a number of human diseases and limits the yields of pharmaceutical proteins during production. Computer simulations can be used to develop novel experimentally testable hypotheses pertaining to aggregation. While all-atom simulations with explicit solvent are too computationally intensive to address the multitude of relevant time scales, coarse-grained models make it possible to observe the transition of monomers to an equilibrium containing aggregates. Here, we provide the reader with background information and a list of steps for setting up, performing, and analyzing computer simulations of aggregating coarse-grained (CG) proteins.

Mapping conformational ensembles of aβ oligomers in molecular dynamics simulations

Although the oligomers formed by Aβ peptides appear to be the primary cytotoxic species in Alzheimer's disease, detailed information about their structures appears to be lacking. In this article, we use exhaustive replica exchange molecular dynamics and an implicit solvent united-atom model to study the structural properties of Aβ monomers, dimers, and tetramers. Our analysis suggests that the conformational ensembles of Aβ dimers and tetramers are very similar, but sharply distinct from those sampled by the monomers. The key conformational difference between monomers and oligomers is the formation of β-structure in the oligomers occurring together with the loss of intrapeptide interactions and helix structure. Our simulations indicate that, independent of oligomer order, the Aβ aggregation interface is largely confined to the sequence region 10-23, which forms the bulk of interpeptide interactions. We show that the fractions of β structure computed in our simulations and measured experimentally are in good agreement.

Toward a molecular theory of early and late events in monomer to amyloid fibril formation

Quantitative understanding of the kinetics of fibril formation and the molecular mechanism of transition from monomers to fibrils is needed to obtain insights into the growth of amyloid fibrils and more generally self-assembly multisubunit protein complexes. Significant advances using computations of protein aggregation in a number of systems have established generic and sequence-specific aspects of the early steps in oligomer formation. Theoretical considerations, which view oligomer and fibril growth as diffusion in a complex energy landscape, and computational studies, involving minimal lattice and coarse-grained models, have revealed general principles governing the transition from monomeric protein to ordered fibrillar aggregates. Detailed atomistic calculations have explored the early stages of the protein aggregation pathway for a number of amyloidogenic proteins, most notably amyloid β- (Aβ-) protein and fragments from proteins linked to various diseases. These computational studies have provided insights into the role of sequence, role of water, and specific interatomic interactions underlying the thermodynamics and dynamics of elementary kinetic steps in the aggregation pathway. Novel methods are beginning to illustrate the structural basis for the production of Aβ-peptides through interactions with secretases in the presence of membranes. We show that a variety of theoretical approaches, ranging from scaling arguments to minimal models to atomistic simulations, are needed as a complement to experimental studies probing the principles governing protein aggregation.

Computational backbone mutagenesis of Abeta peptides: probing the role of backbone hydrogen bonds in aggregation

Using replica exchange molecular dynamics (REMD) and united atom implicit solvent model we examine the role of backbone hydrogen bonds (HBs) in Abeta aggregation. The importance of HBs appears to depend on the aggregation stage. The backbone HBs have little effect on the stability of Abeta dimers or on their aggregation interface. The HBs also do not play a critical role in initial binding of Abeta peptides to the amyloid fibril. Their elimination does not change the continuous character of Abeta binding nor its temperature. However, cancellation of HBs forming between incoming Abeta peptides and the fibril disrupts the locked fibril-like states in the bound peptides. Without the support of HBs, bound Abeta peptides form few long beta-strands on the fibril edge. As a result, the deletion of peptide-fibril HBs is expected to impede fibril growth. As for the peptides bound to Abeta fibril the deletion of interpeptide HBs reduces the beta propensity in the dimers making them less competent for amyloid assembly. These simulation findings together with the backbone mutagenesis experiments suggest that a viable strategy for arresting fibril growth is the disruption of interpeptide HBs.

Globular state in the oligomers formed by Abeta peptides

Replica exchange molecular dynamics and implicit solvent model are used to study two oligomeric species of Abeta peptides, dimer and tetramer, which are typically observed in in vitro experiments. Based on the analysis of free energy landscapes, density distributions, and chain flexibility, we propose that the oligomer formation is a continuous transition occurring without metastable states. The density distribution computations suggest that Abeta oligomer consists of two volume regions-the core with fairly flat density profile and the surface layer with rapidly decreasing density. The core is mostly formed by the N-terminal residues, whereas the C-terminal tends to occur in the surface layer. Lowering the temperature results in the redistribution of peptide atoms from the surface layer into the core. Using these findings, we argue that Abeta oligomer resembles polymer globule in poor solvent. Abeta dimers and tetramers are found to be structurally similar suggesting that the conformations of Abeta peptides do not depend on the order of small oligomers.

Amyloid-beta42 oligomer structures from fibrils: a systematic molecular dynamics study

Recent experimental data demonstrate that small, soluble amyloid-beta42 oligomers play an important role in Alzheimer's disease because they exhibit neurotoxic properties and also act as seed for fibril growth. We performed all-atom molecular dynamics simulations in explicit solvent of 0.7 micros in total on five Abeta9-42 oligomers (monomer through pentamer) starting from the fibril conformation. The initial conformation proves to be stable in the trimer to pentamer, and the two parallel in-register beta-sheets as well as the connecting turn are preserved. The dimer undergoes larger conformational changes in its C-terminus, and the predominant conformation detected exhibits an additional antiparallel beta-sheet in one of the subunits. This conformational rearrangement allows efficient shielding of hydrophobic residues from the solvent, which is not possible for a dimer in the fibril conformation. In addition to the presence of the hydrogen bonds in the beta-sheets, the larger oligomers are stabilized by interchain D23-K28 salt bridges, whereas a D23-N27 interaction is found in the dimer. The degree of structural similarity to the fibril conformation detected for the oligomers in the simulation may also offer a structural explanation for the experimental finding that trimers and tetramers act as more potent seeds in fibril formation than dimers because only small conformational changes will be required for fibril growth. The fact that the dimer predominantly exists in conformations distinct from the larger oligomers and the fibril is also interesting for the design of anti-Alzheimer drugs, because it suggests that multiple drugs might be required to target the structurally different neurotoxic oligomers.