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

Effect of caffeine on the aggregation of amyloid-β-A 3D RISM study

Alzheimer's disease is a detrimental neurological disorder caused by the formation of amyloid fibrils due to the aggregation of amyloid-β peptide. The primary therapeutic approaches for treating Alzheimer's disease are targeted to prevent this amyloid fibril formation using potential inhibitor molecules. The discovery of such inhibitor molecules poses a formidable challenge to the design of anti-amyloid drugs. This study investigates the effect of caffeine on dimer formation of the full-length amyloid-β using a combined approach of all-atom, explicit water molecular dynamics simulations and the three-dimensional reference interaction site model theory. The change in the hydration free energy of amyloid-β dimer, with and without the inhibitor molecules, is calculated with respect to the monomeric amyloid-β, where the hydration free energy is decomposed into energetic and entropic components, respectively. Dimerization is accompanied by a positive change in the partial molar volume. Dimer formation is spontaneous, which implies a decrease in the hydration free energy. However, a reverse trend is observed for the dimer with inhibitor molecules. It is observed that the negatively charged residues primarily contribute for the formation of the amyloid-β dimer. A residue-wise decomposition reveals that hydration/dehydration of the side-chain atoms of the charged amino acid residues primarily contribute to dimerization.

A comprehensive computer simulation insight into inhibitory mechanisms of EGCG and NQTrp ligands on amyloid-beta assemblies as the Alzheimer's disease insignia

Amyloid-β peptide with predominant presence in the senile plaques is the most common agent for Alzheimer's disease (AD) incidence. Assembly of the amyloid-β(1-42) (Aβ1-42) isoform is known as the main reason for the AD appearance. Epigallocatechin gallate (EGCG) and 1,4-naphthoquinone-2-yl-L-tryptophan (NQTrp) are two small molecules that inhibit the formation of the Aβ1-42 fibrils. The present study provides molecular insight to clarify the inhibitory mechanisms of the EGCG and NQTrp ligands on the Aβ1-42 assemblies by using molecular dynamics (MD) simulation. Hence, nine different Aβ1-42-containing systems including the monomer, dimer, and hexamer of Aβ1-42 considering each of them in a media with no ligands, in the presence of one EGCG ligand, and in the presence of one EGCG ligand were studied with a simulation time of 1 µs for each system. The precise investigation of the peptide-ligand distance, conformational factor (Pi), solvent accessible surface area (SASA), dictionary of secondary structure (DSSP), and Lys28-Ala42 salt bridge analyses confirmed that the hydroxyl-rich structure of the EGCG ligand applied its inhibitory effect on the aggregation of the peptides indirectly by involving water molecules. While the hydroxyl-free structure of the NQTrp ligand exposed its inhibitory effect through a direct interaction with the Aβ1-42 peptides. Besides, reduced density gradient (RDG) analysis clarified the hydrogen bonding interactions as the dominant ones for the peptide-EGCG systems, and also, steric and van der Waals interactions for the peptide-NQTrp systems.Communicated by Ramaswamy H. Sarma.

Hydroxytyrosol-Donepezil Hybrids Play a Protective Role in an In Vitro Induced Alzheimer's Disease Model and in Neuronal Differentiated Human SH-SY5Y Neuroblastoma Cells

Alzheimer's disease (AD) is the most common neurodegenerative pathology among progressive dementias, and it is characterized by the accumulation in the brain of extracellular aggregates of beta-amyloid proteins and neurofibrillary intracellular tangles consisting of τ-hyperphosphorylated proteins. Under normal conditions, beta-amyloid peptides exert important trophic and antioxidant roles, while their massive presence leads to a cascade of events culminating in the onset of AD. The fibrils of beta-amyloid proteins are formed by the process of fibrillogenesis that, starting from individual monomers of beta-amyloid, can generate polymers of this protein, constituting the hypothesis of the "amyloid cascade". To date, due to the lack of pharmacological treatment for AD without toxic side effects, chemical research is directed towards the realization of hybrid compounds that can act as an adjuvant in the treatment of this neurodegenerative pathology. The hybrid compounds used in this work include moieties of a hydroxytyrosol, a nitrohydroxytyrosol, a tyrosol, and a homovanillyl alcohol bound to the N-benzylpiperidine moiety of donepezil, the main drug used in AD. Previous experiments have shown different properties of these hybrids, including low toxicity and antioxidant and chelating activities. The purpose of this work was to test the effects of hybrid compounds mixed with Aβ1-40 to induce fibrillogenesis and mimic AD pathogenesis. This condition has been studied both in test tubes and by an in vitro model of neuronal differentiated human SH-SY5Y neuroblastoma cells. The results obtained from test tube experiments showed that some hybrids inhibit the activity of the enzymes AChE, BuChE, and BACE-1. Cell experiments suggested that hybrids could inhibit fibrillogenesis, negatively modulating caspase-3. They were also shown to exert antioxidant effects, and the acetylated hybrids were found to be more functional and efficient than nonacetylated forms.

Reducing the assemblies of amyloid-beta multimers by sodium dodecyl sulfate surfactant at concentrations lower than critical micelle concentration: molecular dynamics simulation exploration

Amyloid-β peptide, the predominant proteinaceous component of senile plaques, is responsible for the incidence of Alzheimer's disease (AD), an age-associated neurodegenerative disorder. Specifically, the amyloid-β(1-42) (Aβ1-42) isoform, known for its high toxicity, is the predominant biomarker for the preliminary diagnosis of AD. The aggregation of the Aβ1-42 peptides can be affected by the components of the cellular medium through changing their structures and molecular interactions. In this study, we investigated the effect of sodium dodecyl sulfate (SDS) at much lower concentrations than the critical micelle concentration (CMC) on Aβ1-42 aggregation. For this purpose, we studied mono-, di-, tri- and tetramers of Aβ1-42 peptide in two different concentrations of SDS molecules (10 and 40 molecules) using a 300 ns molecular dynamics simulation for each system. The distance between the center of mass (COM) of Aβ1-42 peptides confirms that an increase in the number of SDS molecules decreases their aggregation probability due to greater interaction with SDS molecules. Besides, the less compactness parameter reveals the reduced aggregation probability of Aβ1-42 peptides. Based on the energetic FEL landscapes, SDS molecules with the concentration closer to the CMC are an effective inhibitory agent to prevent the formation of Aβ1-42 fibrils. Also, the aggregation direction of the peptide pairs can be predicted by determining the direction of the accumulation-deterrent forces.Communicated by Ramaswamy H. Sarma.

Exploring the aggregation of amyloid-β 42 through Monte Carlo simulations

Coarse-grained Monte Carlo simulations are performed for a disordered protein, amyloid-β 42 to identify the interactions and understand the mechanism of its aggregation. A statistical potential is developed from a selected dataset of intrinsically disordered proteins, which accounts for the respective contributions of the bonded and non-bonded potentials. While, the bonded potential comprises the bond, bend, and dihedral constraints, the nonbonded interactions include van der Waals interactions, hydrogen bonds, and the two-body potential. The two-body potential captures the features of both hydrophobic and electrostatic interactions that brings the chains at a contact distance, while the repulsive van der Waals interactions prevent them from a collapse. Increased two-body hydrophobic interactions facilitate the formation of amorphous aggregates rather than the fibrillar ones. The formation of aggregates is validated from the interchain distances, and the total energy of the system. The aggregate is structurally characterized by the root-mean-square deviation, root-mean-square fluctuation and the radius of gyration. The aggregates are characterized by a decrease in SASA, an increase in the non-local interactions and a distinct free energy minimum relative to that of the monomeric state of amyloid-β 42. The hydrophobic residues help in nucleation, while the charged residues help in oligomerization and aggregation.

How do salt and lipids affect conformational dynamics of Aβ42 monomers in water?

It is well established that amyloid β-protein (Aβ) self-assembly is involved in triggering of Alzheimer's disease. On the other hand, evidence of physiological function of Aβ interacting with lipids has only begun to emerge. Details of Aβ-lipid interactions, which may underlie physiological and pathological activities of Aβ, are not well understood. Here, the effects of salt and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids on conformational dynamics of Aβ42 monomer in water are examined by all-atom molecular dynamics (MD). We acquired six sets of 250 ns long MD trajectories for each of the three lipid concentrations (0, 27, and 109 mM) in the absence and presence of 150 mM salt. Ten replica trajectories per set are used to enhance sampling of Aβ42 conformational space. We show that salt facilitates long-range tertiary contacts in Aβ42, resulting in more compact Aβ42 conformations. By contrast, addition of lipids results in lipid-concentration dependent Aβ42 unfolding concomitant with enhanced stability of the turn in the A21-A30 region. At the high lipid concentration, salt enables the N-terminal region of Aβ42 to form long-range tertiary contacts and interact with lipids, which results in formation of a parallel β-strand. Aβ42 forms stable lipid-protein complexes whereby the protein is adhered to the lipid cluster rather than embedded into it. We propose that the inability of Aβ42 monomer to get embedded into the lipid cluster may be important for facilitating repair of leaks in the blood-brain barrier without penetrating and damaging cellular membranes.

Monitoring early-stage β-amyloid dimer aggregation by histidine site-specific two-dimensional infrared spectroscopy in a simulation study

Monitoring early-stage β-amyloid (Aβ) dimerization is a formidable challenge for understanding neurological diseases. We compared β-sheet formation and histidine site-specific two-dimensional infrared (2D IR) spectroscopic signatures of Aβ dimers with different histidine states (δ; N^δ1-H, ε; N^ε2-H, or π; both protonated). Molecular dynamics (MD) simulations revealed that β-sheet formation is favored for the δδδ:δδδ and πππ:πππ tautomeric isomers showing strong couplings and frequent contacts between the central hydrophobic core and C-terminus compared with the εεε:εεε isomer. Characteristic blue-shifts in the 2D IR central bands were observed upon monomer-dimer transformation. The εεε:εεε dimer exhibited larger frequency shifts than δδδ:δδδ and πππ:πππ implying that the red-shift may have a correlation with N^δ1-H(δ) protonation. Our results support the tautomerization/protonation hypothesis that attributes Aβ misfolding to histidine tautomers as a possible primary initiator for Aβ aggregation and facilitates the application of histidine site-specific 2D IR spectroscopy for studying early-stage Aβ self-assembly.

Binary Matrix Method to Enumerate, Hierarchically Order, and Structurally Classify Peptide Aggregation

Protein aggregation is a common and complex phenomenon in biological processes, yet a robust analysis of this aggregation process remains elusive. The commonly used methods such as center-of-mass to center-of-mass (COM-COM) distance, the radius of gyration (R_g), hydrogen bonding (HB), and solvent accessible surface area do not quantify the aggregation accurately. Herein, a new and robust method that uses an aggregation matrix (AM) approach to investigate peptide aggregation in a MD simulation trajectory is presented. An nxn two-dimensional AM is created by using the interpeptide C_α-C_α cutoff distances, which are binarily encoded (0 or 1). These aggregation matrices are analyzed to enumerate, hierarchically order, and structurally classify the aggregates. Comparison of the present AM method suggests that it is superior to the HB method since it can incorporate nonspecific interactions and the R_g and COM-COM methods since the cutoff distance is independent of the length of the peptide. More importantly, the present method can structurally classify the peptide aggregates, which the conventional R_g, COM-COM, and HB methods fail to do. The unique selling point of this method is its ability to structurally classify peptide aggregates using two-dimensional matrices.

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.

Mechanistic Insights into the Co-Aggregation of Aβ and hIAPP: An All-Atom Molecular Dynamic Study

Patients with Alzheimer's disease (AD) have a high risk of developing Type II diabetes (T2D). The co-aggregation of the two disease-related proteins, Aβ and hIAPP, has been proposed as a potential molecular mechanism. However, the detailed Aβ-hIAPP interactions and structural characteristics of co-aggregates are mostly unknown at atomic level. Here, we explore the conformational ensembles of the Aβ-hIAPP heterodimer and Aβ or hIAPP homodimer by performing all-atom explicit-solvent replica exchange molecular dynamic simulations. Our simulations show that the interaction propensity of Aβ-hIAPP in the heterodimer is comparable with that of Aβ-Aβ/hIAPP-hIAPP in the homodimer. Similar hot spot residues of Aβ/hIAPP in the homodimer and heterodimer are identified, indicating that both Aβ and hIAPP have similar molecular recognition sites for self-aggregation and co-aggregation. Aβ in the heterodimer possesses three high β-sheet probability regions: the N-terminal region E3-H6, the central hydrophobic core region K16-E22, and the C-terminal hydrophobic region I31-A41, which is highly similar to Aβ in the homodimer. More importantly, in the heterodimer, the regions E3-H6, F19-E22, and I31-M35 of Aβ and the amyloid core region N20-T30 of hIAPP display higher β-sheet probability than they do in homodimer, implying their crucial roles in the formation of β-sheet-rich co-aggregates. Our study sheds light on the co-aggregation of Aβ and hIAPP at an atomic level, which will be helpful for an in-depth understanding of the molecular mechanism for epidemiological correlation of AD and T2D.

Compact fibril-like structure of amyloid β-peptide (1-42) monomers

Amyloid β (Aβ) monomers are the smallest assembly units, and play an important role in most of the individual processes involved in amyloid fibril formation. An important question is whether the monomer can adopt transient fibril-like conformations in solution. Here we use enhanced sampling simulations to study the Aβ42 monomer structural flexibility. We show that the monomer frequently adopts conformations with the N-terminus region structured very similarly to the conformation it adopts inside the fibril. This intrinsic propensity of monomeric Aβ to adopt fibril-like conformations could explain the low free energy barrier for Aβ42 fibril elongation.

Cross-Linked Amyloid β-Protein Oligomers: A Missing Link in Alzheimer's Disease Pathology?

Amyloid β-protein (Aβ) oligomers are broadly viewed as the proximate mediators of toxicity in Alzheimer's disease (AD). Recent studies, however, provide substantial evidence that Aβ is involved in protection and repair of the central nervous system whereby Aβ oligomer and subsequent fibril formation are integral to its normal antimicrobial and antiviral function. These developments raise a question of what exactly makes Aβ oligomers toxic in the context of AD. This Perspective describes a paradigm shift in the search for toxic Aβ oligomer species that involves oxidative-stress-induced stabilization of Aβ oligomers via cross-linking and reviews most recent research elucidating structural aspects of cross-linked Aβ oligomers and potential inhibition of their toxicity.

Distinct Primary Nucleation of Polymorphic Aβ Dimers Yields to Distinguished Fibrillation Pathways

Polymorphic Aβ dimers are the smallest toxic species that play a role in the pathology of Alzheimer's disease. There is great interest in understanding the malfunctions that yield to these toxic species and in providing insights into the molecular mechanisms of the primary nucleation. Herein, we present a first work that demonstrates two distant edges states of Aβ dimers. The first is the so-called "random coil" state dimer that mimics the primary seeding/nucleation that is far from a fibrillation state. The second is the "fibril-like" state dimer that is structurally in close proximity to the fibril, a well-organized state into a fibril-like structure. We show for the first time that a conformational change of one monomer within the dimer impedes primary nucleation, while less fluctuations and relatively large number of interactions in nucleation domains induce the primary nucleation to produce toxic stable species. Overall, the current study exhibits a diversity of primary nucleation in each dimer state, suggesting distinct molecular mechanisms of fibril formation. The conformations of the early stage Aβ dimers that were achieved may provide crucial data for designing inhibitors to impede the primary nucleation.

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.

Exploring Amyloid-β Dimer Structure Using Molecular Dynamics Simulations

A major hallmark of Alzheimer's disease (AD) is the aggregation of amyloid-β peptides in the brains of people afflicted by the disease. The exact pathway to this catastrophic event is unknown. In this work, a total of 9.5 μs molecular dynamics simulations have been performed to investigate the structure and dynamics of the smallest form of toxic Aβ oligomers, i.e., the Aβ dimers. This study suggests that specific hydrophobic regions are vital in the aggregation process. Different possible structures for Aβ dimers are reported along with their relative binding affinity. These data may be used to design better Aβ-aggregation inhibitors. The diversity of the dimer structures suggests several aggregation pathways.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Top-down Multiscale Approach To Simulate Peptide Self-Assembly from Monomers

Modeling peptide assembly from monomers on large time and length scales is often intractable at the atomistic resolution. To address this challenge, we present a new approach which integrates coarse-grained (CG), mixed-resolution, and all-atom (AA) modeling in a single simulation. We simulate the initial encounter stage with the CG model, while the further assembly and reorganization stages are simulated with the mixed-resolution and AA models. We have implemented this top-down approach with new tools to automate model transformations and to monitor oligomer formations. Further, a theory was developed to estimate the optimal simulation length for each stage using a model peptide, melittin. The assembly level, the oligomer distribution, and the secondary structures of melittin simulated by the optimal protocol show good agreement with prior experiments and AA simulations. Finally, our approach and theory have been successfully validated with three amyloid peptides (β-amyloid 16-22, GNNQQNY fragment from the yeast prion protein SUP35, and α-synuclein fibril 35-55), which highlight the synergy from modeling at multiple resolutions. This work not only serves as proof of concept for multiresolution simulation studies but also presents practical guidelines for further self-assembly simulations at more physically and chemically relevant scales.

Molecular dynamics simulations reveal the mechanism of graphene oxide nanosheet inhibition of Aβ1–42 peptide aggregation

Graphene oxide nanosheets inhibit Aβ1–42 aggregation by weakening inter-peptide interactions and reducing β-sheet contents mostly via salt bridge, hydrogen bonding and cation–π interactions with charged residues.

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.

Aβ under stress: the effects of acidosis, Cu2+-binding, and oxidation on amyloid β-peptide dimers

In light of the high affinity of Cu2+ for Alzheimer's Aβ1-42 and its ability to subsequently catalyze the formation of radicals, we examine the effects of Cu2+ binding, Aβ oxidation, and an acidic environment on the conformational dynamics of the smallest Aβ1-42 oligomer, the Aβ1-42 dimer. Transition networks calculated from Hamiltonian replica exchange molecular dynamics (H-REMD) simulations reveal that the decreased pH considerably increased the β-sheet content, whereas Cu2+ binding increased the exposed hydrophobic surface area, both of which can contribute to an increased oligomerization propensity and toxicity.

Recent Progress in Alzheimer's Disease Research, Part 1: Pathology

The field of Alzheimer's disease (AD) research has grown exponentially over the past few decades, especially since the isolation and identification of amyloid-β from postmortem examination of the brains of AD patients. Recently, the Journal of Alzheimer's Disease (JAD) put forth approximately 300 research reports which were deemed to be the most influential research reports in the field of AD since 2010. JAD readers were asked to vote on these most influential reports. In this 3-part review, we review the results of the 300 most influential AD research reports to provide JAD readers with a readily accessible, yet comprehensive review of the state of contemporary research. Notably, this multi-part review identifies the "hottest" fields of AD research providing guidance for both senior investigators as well as investigators new to the field on what is the most pressing fields within AD research. Part 1 of this review covers pathogenesis, both on a molecular and macro scale. Part 2 review genetics and epidemiology, and part 3 covers diagnosis and treatment. This part of the review, pathology, reviews amyloid-β, tau, prions, brain structure, and functional changes with AD and the neuroimmune response of AD.

Pathways of Amyloid-β Aggregation Depend on Oligomer Shape

One of the main research topics related to Alzheimer's disease is the aggregation of the amyloid-β peptide, which was shown to follow different pathways for the two major alloforms of the peptide, Aβ40 and the more toxic Aβ42. Experimental studies emphasized that oligomers of specific sizes appear in the early aggregation process in different quantities and might be the key toxic agents for each of the two alloforms. We use transition networks derived from all-atom molecular dynamics simulations to show that the oligomers leading to the type of oligomer distributions observed in experiments originate from compact conformations. Extended oligomers, on the other hand, contribute more to the production of larger aggregates thus driving the aggregation process. We further demonstrate that differences in the aggregation pathways of the two Aβ alloforms occur as early as during the dimer stage. The higher solvent-exposure of hydrophobic residues in Aβ42 oligomers contributes to the different aggregation pathways of both alloforms and also to the increased cytotoxicity of Aβ42.

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.

Self-Assembly Pathways of β-Sheet-Rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations

Early oligomerization during amyloid-β (Aβ) aggregation is essential for Aβ neurotoxicity. Understanding how unstructured Aβs assemble into oligomers, especially those rich in β-sheets, is essential but remains challenging as the assembly process is too transient for experimental characterization and too slow for molecular dynamics simulations. So far, atomic simulations are limited only to studies of either oligomer structures or assembly pathways for short Aβ segments. To overcome the computational challenge, we combine in this study a hybrid-resolution model and adaptive sampling techniques to perform over 2.7 ms of simulations of formation of full-length Aβ40 dimers that are the earliest toxic oligomeric species. The Markov state model is further employed to characterize the transition pathways and associated kinetics. Our results show that for two major forms of β-sheet-rich structures reported experimentally, the corresponding assembly mechanisms are markedly different. Hairpin-containing structures are formed by direct binding of soluble Aβ in β-hairpin-like conformations. Formation of parallel, in-register structures resembling fibrils occurs ∼100-fold more slowly and involves a rapid encounter of Aβ in arbitrary conformations followed by a slow structural conversion. The structural conversion proceeds via diverse pathways but always requires transient unfolding of encounter complexes. We find that the transition kinetics could be affected differently by intra-/intermolecular interactions involving individual residues in a conformation-dependent manner. In particular, the interactions involving Aβ's N-terminal part promote the assembly into hairpin-containing structures but delay the formation of fibril-like structures, thus explaining puzzling observations reported previously regarding the roles of this region in the early assembly process.

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.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimer's disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimer's disease.

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.

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.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Self-assembly of the full-length amyloid Aβ42 protein in dimers

The self-assembly of amyloid (Aβ) proteins into nano-aggregates is a hallmark of Alzheimer's disease (AD) development, yet the mechanism of how disordered monomers assemble into aggregates remains elusive. Here, we applied long-time molecular dynamics simulations to fully characterize the assembly of Aβ42 monomers into dimers. Monomers undergo conformational changes during their interaction, but the resulting dimer structures do not resemble those found in fibril structures. To identify natural conformations of dimers among a set of simulated ones, validation approaches were developed and applied, and a subset of dimer conformations were characterized. These dimers do not contain long β-strands that are usually found in fibrils. The dimers are stabilized primarily by interactions within the central hydrophobic regions and the C-terminal regions, with a contribution from local hydrogen bonding. The dimers are dynamic, as evidenced by the existence of a set of conformations and by the quantitative analyses of the dimer dissociation process.

Impact of Cu(II) Binding on Structures and Dynamics of Aβ42 Monomer and Dimer: Molecular Dynamics Study

The classical force field, which is compatible with the Amber force field 99SB, has been obtained for the interaction of Cu(II) with monomer and dimers of amyloid-β peptides using the coordination where Cu(II) is bound to His6, His13 (or His14), and Asp1 with distorted planar geometry. The newly developed force field and molecular dynamics simulation were employed to study the impact of Cu(II) binding on structures and dynamics of Aβ₄₂ monomer and dimers. It was shown that in the presence of Cu(II) the β content of monomer is reduced substantially compared with the wild-type Aβ₄₂ suggesting that, in accord with experiments, metal ions facilitate formation of amorphous aggregates rather than amyloid fibrils with cross-β structures. In addition, one possible mechanism for amorphous assembly is that the Asp23-Lys28 salt bridge, which plays a crucial role in β sheet formation, becomes more flexible upon copper ion binding to the Aβ N-terminus. The simulation of dimers was conducted with the Cu(II)/Aβ stoichiometric ratios of 1:1 and 1:2. For the 1:1 ratio Cu(II) delays the Aβ dimerization process as observed in a number of experiments. The mechanism underlying this phenomenon is associated with slow formation of interchain salt bridges in dimer as well as with decreased hydrophobicity of monomer upon Cu-binding.

Assembly of Stefin B into Polymorphic Oligomers Probed by Discrete Molecular Dynamics

Assembly of an amyloidogenic protein stefin B into molten globule oligomers is studied by efficient discrete molecular dynamics. Consistent with in vitro findings, tetramers form primarily through dimer association, resulting in a decreased trimer abundance. Oligomers up to heptamers display elongated rod-like morphologies akin to protofibrils, whereas larger oligomers, decamers through dodecamers, form elongated, branched, as well as annular structures, providing structural insights into pore forming ability and toxicity of amyloidogenic proteins.

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.

Understanding Amyloid-β Oligomerization at the Molecular Level: The Role of the Fibril Surface

The aggregation of the amyloid β-peptide into fibrils is a complex process that involves mechanisms such as primary and secondary nucleation, fibril elongation and fibril fragmentation. Some of these processes generate neurotoxic Aβ oligomers, which are involved in the development of Alzheimer's disease. Recent experimental studies have emphasized the role of the fibril as a catalytic surface for the production of highly toxic oligomers during secondary nucleation. By using molecular dynamics simulations, we show that it is the hydrophobic fibril region that causes the structural changes required for catalyzing the formation of β-sheet-rich Aβ1-42 oligomers on the fibril surface. These results reveal, for the first time, the molecular basis of the secondary nucleation pathway.

Amino acid substitutions [K16A] and [K28A] distinctly affect amyloid β-protein oligomerization

Amyloid β-protein (A β) assembles into oligomers that play a seminal role in Alzheimer's disease (AD), a leading cause of dementia among the elderly. Despite undisputed importance of A β oligomers, their structure and the basis of their toxicity remain elusive. Previous experimental studies revealed that the [K16A] substitution strongly inhibits toxicity of the two predominant A β alloforms in the brain, A β 40 and A β 42, whereas the [K28A] substitution exerts only a moderate effect. Here, folding and oligomerization of [A16]A β 40, [A28]A β 40, [A16]A β 42, and [A28]A β 42 are examined by discrete molecular dynamics (DMD) combined with a four-bead implicit solvent force field, DMD4B-HYDRA, and compared to A β 40 and A β 42 oligomer formation. Our results show that both substitutions promote A β 40 and A β 42 oligomerization and that structural changes to oligomers are substitution- and alloform-specific. The [K28A] substitution increases solvent-accessible surface area of hydrophobic residues and the intrapeptide N-to-C terminal distance within oligomers more than the [K16A] substitution. The [K16A] substitution decreases the overall β-strand content, whereas the [K28A] substitution exerts only a modest change. Substitution-specific tertiary and quaternary structure changes indicate that the [K16A] substitution induces formation of more compact oligomers than the [K28A] substitution. If the structure-function paradigm applies to A β oligomers, then the observed substitution-specific structural changes in A β 40 and A β 42 oligomers are critical for understanding the structural basis of A β oligomer toxicity and correct identification of therapeutic targets against AD.

The inhibitory mechanism of a fullerene derivative against amyloid-β peptide aggregation: an atomistic simulation study

Alzheimer's disease (AD) is associated with the pathological self-assembly of amyloid-β (Aβ) peptides into β-sheet enriched fibrillar aggregates. Aβ dimers formed in the initial step of Aβ aggregation were reported to be the smallest toxic species. Inhibiting the formation of β-sheet-rich oligomers and fibrils is considered as the primary therapeutic strategy for AD. Previous studies reported that fullerene derivatives strongly inhibit Aβ fibrillation. However, the underlying inhibitory mechanism remains elusive. As a first step to understand fullerene-modulated full-length Aβ aggregation, we investigated the conformational ensemble of the Aβ1-42 dimer with and without 1,2-(dimethoxymethano)fullerene (DMF) - a more water-soluble fullerene derivative - by performing a 340 ns explicit-solvent replica exchange molecular dynamics simulation. Our simulations show that although disordered states are the most abundant conformations of the Aβ1-42 dimer, conformations containing diverse extended β-hairpins are also populated. The first most-populated β-hairpins involving residues L17-D23 and A30-V36 strongly resemble the engineered β-hairpin which is a building block of toxic Aβ oligomers. We find that the interaction of DMFs with Aβ peptides greatly impedes the formation of such β-hairpins and inter-peptide β-sheets. Binding energy analyses demonstrate that DMF preferentially binds not only to the central hydrophobic motif LVFFA of the Aβ peptide as suggested experimentally, but also to the aromatic residues including F4 and Y10 and the C-terminal hydrophobic region I31-V40. This study reveals a complete picture of the inhibitory mechanism of full-length Aβ1-42 aggregation by fullerenes, providing theoretical insights into the development of drug candidates against AD.

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.

Soluble Prion Protein Binds Isolated Low Molecular Weight Amyloid-β Oligomers Causing Cytotoxicity Inhibition

A growing number of observations indicate that soluble amyloid-β (Aβ) oligomers play a major role in Alzheimer's disease. Recent studies strongly suggest that at least some of the neurotoxic effects of these oligomers are mediated by cellular, membrane-anchored prion protein and that Aβ neurotoxicity can be inhibited by soluble recombinant prion protein (rPrP) and its fragments. However, the mechanism by which rPrP interacts with Aβ oligomers and prevents their toxicity is largely unknown, and studies in this regard are hindered by the large structural heterogeneity of Aβ oligomers. To overcome this difficulty, here we used photoinduced cross-linking of unmodified proteins (PICUP) to isolate well-defined oligomers of Aβ42 and characterize these species with regard to their cytotoxicity and interaction with rPrP, as well the mechanism by which rPrP inhibits Aβ42 cytotoxicity. Our data shows that the addition of rPrP to the assembling Aβ42 results in a shift in oligomer size distribution, decreasing the population of toxic tetramers and higher order oligomers and increasing the population of nontoxic (and possibly neuroprotective) monomers. Isolated oligomeric species of Aβ42 are cytotoxic to primary neurons and cause permeation of model lipid bilayers. These toxic effects, which are oligomer size-dependent, can be inhibited by the addition of rPrP, and our data suggest potential mechanisms of this inhibitory action. This insight should help in current efforts to develop PrP-based therapeutics for Alzheimer's disease.