Difference in dimer conformation between amyloid-β(1–42) and (1–43) proteins: Replica exchange molecular dynamics simulations in water

Identification of key stabilizing interactions of amyloid-β oligomers based on fragment molecular orbital calculations on macrocyclic β-hairpin peptides

Analyzing the electronic states and inter-/intra-molecular interactions of amyloid oligomers expand our understanding of the molecular basis of Alzheimer's disease and other amyloid diseases. In the current study, several high-resolution crystal structures of oligomeric assemblies of Aβ-derived peptides have been studied by the ab initio fragment molecular orbital (FMO) method. The FMO method provides comprehensive details of the molecular interactions between the residues of the amyloid oligomers at the quantum mechanical level. Based on the calculations, two sequential aromatic residues (F19 and F20) and negatively charged E22 on the central region of Aβ have been identified as key residues in oligomer stabilization and potential interesting pharmacophores for preventing oligomer formation.

Energy landscapes of protein aggregation and conformation switching in intrinsically disordered proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of 'simple' protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.

Molecular insight into the early stage of amyloid-β(1-42) Homodimers aggregation influenced by histidine tautomerism

Aggregated amyloid β-peptide (Aβ) in small oligomeric forms inside the brain causes synaptic function disruption and the development of Alzheimer's disease (AD). Histidine is an important amino acid that may lead to structural changes. Aβ42 monomer chain includes 3 histidine residues that considering two ε and δ tautomers 8 isomers, including (εεε) and (εδδ) could be formed. Molecular dynamics simulation on homodimerization of (εεε) (the most common type of tautomers) and (εδδ) tautomers with different initial configurations using monomer chains from our previous work were performed to uncover the tautomeric behavior of histidine on Aβ42 aggregation in a physiological pH which is still largely unknown and impossible to observe experimentally. We found a higher propensity of forming β-sheet in (εδδ) homodimers and specifically in a greater amount from Aβ42 than from Aβ40. A smaller amount of β-sheet formation was observed for (εεε) homodimers compared with (εδδ). Additionally, interactions in (εδδ) homodimers may indicate the importance of the hydrophobic core and C-/N-terminals during oligomerization. Our findings indicate the important role of the tautomeric effect of histidine and (εδδ) homodimers at the early stage of Aβ aggregation.

Binary Structure of Amyloid Beta Oligomers Revealed by Dual Recognition Mapping

Amyloid beta (Aβ) oligomers are widely considered to be the causative agent of Alzheimer's disease (AD), a progressive neurodegenerative disorder. Determining the structure of oligomers is, therefore, important for understanding the disease and developing therapeutic agents; however, elucidating the structure has been proven difficult due to heterogeneity, noncrystallinity, and variability. Herein, we investigated homo- and hetero-oligomers of Aβ40 and Aβ42 using atomic force microscopy (AFM) and revealed characteristics of the molecular structure. By examining the surface of individual oligomers with sequential N- and C-terminus specific antibody-tethered tips, we simultaneously mapped the N- and C-terminus distributions and the elastic modulus. Interestingly, both the N- and C-termini of Aβ peptides were recognized on the oligomer surface, and the termini detected pixel regions exhibited a lower elastic modulus than silent pixel regions. These two types of regions were randomly distributed on the oligomer surface.

Understanding the Twisted Structure of Amyloid Fibrils via Molecular Simulations

Accumulation and aggregation of amyloid are associated with the pathogenesis of many human diseases, such as Alzheimer's disease and Type 2 diabetes mellitus. Therefore, a quantitative understanding of the molecular mechanisms causing different aggregated structures and biomechanical properties of amyloid fibrils could shed some light into the progression of these diseases. In this work, we develop coarse-grained molecular dynamics (CGMD) models to simulate the dynamic self-assembly of two types of amyloids (amylin and amyloid β (Aβ)). We investigate the structural and mechanical properties of different types of aggregated amyloid fibrils. Our simulations demonstrate that amyloid fibrils could result from longitudinal growth of protofilament bundles, confirming one of the hypotheses on the fibril formation. In addition, we find that the persistence length of amylin fibrils increases concurrently with their pitch length, suggesting that the bending stiffness of amylin fibrils becomes larger when the amylin fibrils are less twisted. Similar results are observed for Aβ fibrils. These findings quantify the connection between the structural and the biomechanical properties of the fibrils. The CGMD models developed in this work can be potentially used to examine efficacy of anti-aggregation drugs, which could help in developing new treatments.

Initial Structural Models of the Aβ42 Dimer from Replica Exchange Molecular Dynamics Simulations

Experimental characterization of the molecular structure of small amyloid (A)β oligomers that are currently considered as toxic agents in Alzheimer's disease is a formidably difficult task due to their transient nature and tendency to aggregate. Such structural information is of importance because it can help in developing diagnostics and an effective therapy for the disease. In this study, molecular simulations and protein-protein docking are employed to explore a possible connection between the structure of Aβ monomers and the properties of the intermonomer interface in the Aβ42 dimer. A structurally diverse ensemble of conformations of the monomer was sampled in microsecond timescale implicit solvent replica exchange molecular dynamics simulations. Representative structures with different solvent exposure of hydrophobic residues and secondary structure content were selected to build structural models of the dimer. Analysis of these models reveals that formation of an intramonomer salt bridge (SB) between Asp23 and Lys28 residues can prevent the building of a hydrophobic interface between the central hydrophobic clusters (CHCs) of monomers upon dimerization. This structural feature of the Aβ42 dimer is related to the difference in packing of hydrophobic residues in monomers with the Asp23-Lys28 SB in on and off states, in particular, to a lower propensity to form hydrophobic contacts between the CHC domain and C-terminal residues in monomers with a formed SB. These findings could have important implications for understanding the difference between aggregation pathways of Aβ monomers leading to neurotoxic oligomers or inert fibrillar structures.

EGCG inhibits the oligomerization of amyloid beta (16-22) hexamer: Theoretical studies

An extensive replica exchange molecular dynamics (REMD) simulation was performed to investigate the progress patterns of the inhibition of (-)-epigallocatechin-3-gallate (EGCG) on the Aβ_16-22 hexamer. Structural variations of the oligomers without and with EGCG were monitored and analyzed in detail. It has been found that EGCG prevents the formation of Aβ oligomer through two different ways by either accelerating the Aβ oligomerization or reducing the β-content of the hexamer. It also decreases the potential "highly toxic" conformations of Aβ oligomer, which is related to the conformations having high order β-sheet sizes. Both electrostatic and van der Waals interaction energies are found to be involved to the binding process. Computed results using quantum chemical methods show that the π-π stacking is a critical factor of the interaction between EGCG and the peptides. As a result, the binding free energy of the EGCG to the Aβ peptides is slightly larger than that of the curcumin.

Hydrodynamic effects on β-amyloid (16-22) peptide aggregation

Computer simulations based on simplified representations are routinely used to explore the early steps of amyloid aggregation. However, when protein models with implicit solvent are employed, these simulations miss the effect of solvent induced correlations on the aggregation kinetics and lifetimes of metastable states. In this work, we apply the multi-scale Lattice Boltzmann Molecular Dynamics technique (LBMD) to investigate the initial aggregation phases of the amyloid Aβ16-22 peptide. LBMD includes naturally hydrodynamic interactions (HIs) via a kinetic on-lattice representation of the fluid kinetics. The peptides are represented by the flexible OPEP coarse-grained force field. First, we have tuned the essential parameters that control the coupling between the molecular and fluid evolutions in order to reproduce the experimental diffusivity of elementary species. The method is then deployed to investigate the effect of HIs on the aggregation of 100 and 1000 Aβ16-22 peptides. We show that HIs clearly impact the aggregation process and the fluctuations of the oligomer sizes by favouring the fusion and exchange dynamics of oligomers between aggregates. HIs also guide the growth of the leading largest cluster. For the 100 Aβ16-22 peptide system, the simulation of ∼300 ns allowed us to observe the transition from ellipsoidal assemblies to an elongated and slightly twisted aggregate involving almost the totality of the peptides. For the 1000 Aβ16-22 peptides, a system of unprecedented size at quasi-atomistic resolution, we were able to explore a branched disordered fibril-like structure that has never been described by other computer simulations, but has been observed experimentally.

Replica exchange molecular dynamics study of the amyloid beta (11–40) trimer penetrating a membrane

The transmembrane Aβ_11–40 trimer is investigated for the first time using REMD and FEP.

Replica exchange molecular dynamics study of the truncated amyloid beta (11–40) trimer in solution

The structure of the 3Aβ_11–40 oligomer is determined for the first time using T-REMD simulations.

An Account of Amyloid Oligomers: Facts and Figures Obtained from Experiments and Simulations

The deposition of amyloid in brain tissue in the context of neurodegenerative diseases involves the formation of intermediate species-termed oligomers-of lower molecular mass and with structures that deviate from those of mature amyloid fibrils. Because these oligomers are thought to be primarily responsible for the subsequent disease pathogenesis, the elucidation of their structure is of enormous interest. Nevertheless, because of the high aggregation propensity and the polydispersity of oligomeric species formed by the proteins or peptides in question, the preparation of appropriate samples for high-resolution structural methods has proven to be rather difficult. This is why theoretical approaches have been of particular importance in gaining insights into possible oligomeric structures for some time. Only recently has it been possible to achieve some progress with regard to the experimentally based structural characterization of defined oligomeric species. Here we discuss how theory and experiment are used to determine oligomer structures and what can be done to improve the integration of the two disciplines.

Advances in the Simulation of Protein Aggregation at the Atomistic Scale

Protein aggregation into highly structured amyloid fibrils is associated with various diseases including Alzheimer's disease, Parkinson's disease, and type II diabetes. Amyloids can also have normal biological functions and, in the future, could be used as the basis for novel nanoscale materials. However, a full understanding of the physicochemical forces that drive protein aggregation is still lacking. Such understanding is crucial for the development of drugs that can effectively inhibit aberrant amyloid aggregation and for the directed design of functional amyloids. Atomistic simulations can help understand protein aggregation. In particular, atomistic simulations can be used to study the initial formation of toxic oligomers which are hard to characterize experimentally and to understand the difference in aggregation behavior between different amyloidogenic peptides. Here, we review the latest atomistic simulations of protein aggregation, concentrating on amyloidogenic protein fragments, and provide an outlook for the future in this field.