Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach

Delving into the crucial role of the initial structure in the dynamic and self-assembly of amyloid beta

Alzheimer's disease involves the accumulation of amyloid beta (Aβ) monomers that form oligomers and fibrils in the brain. Studying the Aβ monomer is critical for understanding Aβ assembly and peptide behavior and has implications for drug design. Choosing a starting structure with a higher aggregation tendency for cost-effective MD studies and drug design is crucial. Previous studies have utilized distinct initial conformations, leading to varying results. Hence, this study was conducted to compare different initial conformations using the same MD simulation protocol to investigate the behavior and oligomerization propensity of different starting structures of Aβ during 1μs. The behavior of the monomers and their self-assembly systems were studied thoroughly, and the results revealed that highly helical Aβ monomers which used as starting structures retain high helix content during the simulation, and their tautomerization states did not cause significant changes in the structure. On the other hand, the Aβ extended and S-shaped monomers displayed the fingerprints of the fibril structure, which is believed to be more favorable for self-assembly. Self-assembly behaviors were seen for three S-shaped and three Aβ extended peptides. However, both conformations did not show stable β-sheet intermolecular interaction. For the Aβ16-22 monomer as a fragment of the Aβ that can assemble into fibrils, the impacts of capping and uncapping on the initial structure were also investigated. The results displayed that capped and uncapped structures can form oligomers with β-sheet at termini. However, in the capped state, β-sheet interactions were more stable and remained relatively longer than uncapped.

Survey of the Aβ-peptide structural diversity: molecular dynamics approaches

The review deals with the application of Molecular Dynamics (MD) to the structure modeling of beta-amyloids (Aβ), currently classified as intrinsically disordered proteins (IDPs). In this review, we strive to relate the main advances in this area but specifically focus on the approaches and methodology. All relevant papers on the Aβ modeling are cited in the Tables in Supplementary Data, including a concise description of the applied approaches, sorted according to the types of the studied systems: modeling of the monomeric Aβ and Aβ aggregates. Similar sections focused according to the type of modeled object are present in the review. In the final part of the review, novel methods of general IDP modeling not confined to Aβ are described.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01253-y.

Sequence-based identification of amyloidogenic β-hairpins reveals a prostatic acid phosphatase fragment promoting semen amyloid formation

β-Structure-rich amyloid fibrils are hallmarks of several diseases, including Alzheimer's (AD), Parkinson's (PD), and type 2 diabetes (T2D). While amyloid fibrils typically consist of parallel β-sheets, the anti-parallel β-hairpin is a structural motif accessible to amyloidogenic proteins in their monomeric and oligomeric states. Here, to investigate implications of β-hairpins in amyloid formation, potential β-hairpin-forming amyloidogenic segments in the human proteome were predicted based on sequence similarity with β-hairpins previously observed in Aβ, α-synuclein, and islet amyloid polypeptide, amyloidogenic proteins associated with AD, PD, and T2D, respectively. These three β-hairpins, established upon binding to the engineered binding protein β-wrapin AS10, are characterized by proximity of two sequence segments rich in hydrophobic and aromatic amino acids, with high β-aggregation scores according to the TANGO algorithm. Using these criteria, 2505 potential β-hairpin-forming amyloidogenic segments in 2098 human proteins were identified. Characterization of a test set of eight protein segments showed that seven assembled into Thioflavin T-positive aggregates and four formed β-hairpins in complex with AS10 according to NMR. One of those is a segment of prostatic acid phosphatase (PAP) comprising amino acids 185-208. PAP is naturally cleaved into fragments, including PAP(248-286) which forms functional amyloid in semen. We find that PAP(185-208) strongly decreases the protein concentrations required for fibril formation of PAP(248-286) and of another semen amyloid peptide, SEM1(86-107), indicating that it promotes nucleation of semen amyloids. In conclusion, β-hairpin-forming amyloidogenic protein segments could be identified in the human proteome with potential roles in functional or disease-related amyloid formation.

Exploring Structural Insights of Aβ42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-β 42 (Aβ42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aβ42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aβ42 dimer. The binding of α-Syn to Aβ42 reduces the propensity of Aβ42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aβ42 significantly increases the β-content of α-Syn, consistent with the experiments showing that Aβ42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aβ42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aβ42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

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

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

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

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

Molecular Insights into Cu/Zn Metal Response to the Amyloid β-Peptide (1-42)

Aβ1-40 peptide and Aβ1-42 peptide are the building units of beta-amyloid plaques present in Alzheimer's disease (AD)-affected brain. The binding affinity of various divalent metal ions such as Cu and Zn present in AD-affected brain with different amino acids available in Aβ-peptide became the focus to explore their role in soluble neurotoxic oligomer formation. Cu2+ metal ions are known to enhance the neurotoxicity of the Aβ1-42 peptide by catalyzing the formation of soluble neurotoxic oligomers. The competitive preference of both Cu2+ and Zn2+ simultaneously to interact with the Aβ-peptide is unknown. The divalent Cu and Zn ions were inserted in explicit aqueous Aβ1-42 peptide configurations to get insights into the binding competence of these metal ions with peptides using classical molecular dynamics (MD) simulations. The metal-ion interactions reveal that competitive binding preferences of various peptide sites become metal-ion-specific and differ significantly. For Cu2+, interactions are found to be more significant with respect to those of Asp-7, His-6, Glu-11, and His-14. Asp-1, Glu-3, Asp-7, His-6, Glu-11, and His-13 amino acid residues show higher affinity toward Zn2+ ions. MD simulations show notable variation in the solvent-accessible surface area in the hydrophobic region of the peptide. Infinitesimal mobility was obtained for Zn2+ compared to Cu2+ in an aqueous solution and Cu2+ diffusivity deviated significantly at different time scales, proving its labile features in aqueous Aβ1-42 peptides.

Protein Corona Formation and Aggregation of Amyloid β 1-40-Coated Gold Nanocolloids

Amyloid fibrillogenesis is a pathogenic protein aggregation process that occurs through a highly ordered process of protein-protein interactions. To better understand the protein-protein interactions involved in amyloid fibril formation, we formed nanogold colloid aggregates by stepwise additions of ∼2 nmol of amyloid β 1-40 peptide (Aβ1-40) at pH ∼3.7 and ∼25 °C. The processes of protein corona formation and building of gold colloid [diameters (d) of 20 and 80 nm] aggregates were confirmed by a red-shift of the surface plasmon resonance (SPR) band, λpeak, as the number of Aβ1-40 peptides [N(Aβ1-40)] increased. The normalized red-shift of λpeak, Δλ, was correlated with the degree of protein aggregation, and this process was approximated as the adsorption isotherm explained by the Langmuir-Freundlich model. As the coverage fraction (θ) was analyzed as a function of ϕ, which is the N(Aβ1-40) per total surface area of nanogold colloids available for adsorption, the parameters for explaining the Langmuir-Freundlich model were in good agreement for both 20 and 80 nm gold, indicating that ϕ could define the stage of the aggregation process. Surface-enhanced Raman scattering (SERS) imaging was conducted at designated values of ϕ and suggested that a protein-gold surface interaction during the initial adsorption stage may be dependent on the nanosize. The 20 nm gold case seems to prefer a relatively smaller contacting section, such as a -C-N or C═C bond, but a plane of the benzene ring may play a significant role for 80 nm gold. Regardless of the size of the particles, the β-sheet and random coil conformations were considered to be used to form gold colloid aggregates. The methodology developed in this study allows for new insights into protein-protein interactions at distinct stages of aggregation.

Computational insights into the cross-talk between medin and Aβ: implications for age-related vascular risk factors in Alzheimer's disease

The aggregation of medin forming aortic medial amyloid is linked to arterial wall degeneration and cerebrovascular dysfunction. Elevated levels of arteriolar medin are correlated with an increased presence of vascular amyloid-β (Aβ) aggregates, a hallmark of Alzheimer's disease (AD) and vascular dementia. The cross-interaction between medin and Aβ results in the formation of heterologous fibrils through co-aggregation and cross-seeding processes both in vitro and in vivo. However, a comprehensive molecular understanding of the cross-interaction between medin and Aβ-two intrinsically disordered proteins-is critically lacking. Here, we employed atomistic discrete molecular dynamics simulations to systematically investigate the self-association, co-aggregation and also the phenomenon of cross-seeding between these two proteins. Our results demonstrated that both Aβ and medin were aggregation prone and their mixture tended to form β-sheet-rich hetero-aggregates. The formation of Aβ-medin hetero-aggregates did not hinder Aβ and medin from recruiting additional Aβ and medin peptides to grow into larger β-sheet-rich aggregates. The β-barrel oligomer intermediates observed in the self-aggregations of Aβ and medin were also present during their co-aggregation. In cross-seeding simulations, preformed Aβ fibrils could recruit isolated medin monomers to form elongated β-sheets. Overall, our comprehensive simulations suggested that the cross-interaction between Aβ and medin may contribute to their pathological aggregation, given the inherent amyloidogenic tendencies of both medin and Aβ. Targeting medin, therefore, could offer a novel therapeutic approach to preserving brain function during aging and AD by improving vascular health.

Sterol-activated amyloid beta fibril formation

The metabolic processes that link Alzheimer's disease (AD) to elevated cholesterol levels in the brain are not fully defined. Amyloid beta (Aβ) plaque accumulation is believed to begin decades prior to symptoms and to contribute significantly to the disease. Cholesterol and its metabolites accelerate plaque formation through as-yet-undefined mechanisms. Here, the mechanism of cholesterol (CH) and cholesterol 3-sulfate (CS) induced acceleration of Aβ42 fibril formation is examined in quantitative ligand binding, Aβ42 fibril polymerization, and molecular dynamics studies. Equilibrium and pre-steady-state binding studies reveal that monomeric Aβ42•ligand complexes form and dissociate rapidly relative to oligomerization, that the ligand/peptide stoichiometry is 1-to-1, and that the peptide is likely saturated in vivo. Analysis of Aβ42 polymerization progress curves demonstrates that ligands accelerate polymer synthesis by catalyzing the conversion of peptide monomers into dimers that nucleate the polymerization reaction. Nucleation is accelerated ∼49-fold by CH, and ∼13,000-fold by CS - a minor CH metabolite. Polymerization kinetic models predict that at presumed disease-relevant CS and CH concentrations, approximately half of the polymerization nuclei will contain CS, small oligomers of neurotoxic dimensions (∼12-mers) will contain substantial CS, and fibril-formation lag times will decrease 13-fold relative to unliganded Aβ42. Molecular dynamics models, which quantitatively predict all experimental findings, indicate that the acceleration mechanism is rooted in ligand-induced stabilization of the peptide in non-helical conformations that readily form polymerization nuclei.

Co-aggregation of α-synuclein with amyloid-β stabilizes β-sheet-rich oligomers and enhances the formation of β-barrels

Alzheimer's disease (AD) and Parkinson's disease (PD) are the two most common neurodegenerative diseases with markedly different pathological features of β-amyloid (Aβ) plaques and α-synuclein (αS) Lewy bodies (LBs), respectively. However, clinical overlaps in symptoms and pathologies between AD and PD are commonly observed caused by the cross-interaction between Aβ and αS. To uncover the molecular mechanisms behind their overlapping symptoms and pathologies, we computationally investigated the impact of αS on an Aβ monomer and dimerization using atomistic discrete molecular dynamics simulations (DMD). Our results revealed that αS could directly interact with Aβ monomers and dimers, thus forming β-sheet-rich oligomers, including potentially toxic β-barrel intermediates. The binding hotspot involved the second half of the N-terminal domain and NAC region in αS, along with residues 10-21 and 31-42 in Aβ. In their hetero-complex, the binding hotspot primarily assumed a β-sheet core buried inside, which was dynamically shielded by the highly charged, amyloid-resistant C-terminus of αS. Because the amyloid prion region was the same as the binding hotspot being buried, their fibrillization may be delayed, causing the toxic oligomers to increase. This study sheds light on the intricate relationship between Aβ and αS and provides insights into the overlapping pathology of AD and PD.

Nucleation of α-Synuclein Amyloid Fibrils Induced by Cross-Interaction with β-Hairpin Peptides Derived from Immunoglobulin Light Chains

Heterologous interactions between different amyloid-forming proteins, also called cross-interactions, may have a critical impact on disease-related amyloid formation. β-hairpin conformers of amyloid-forming proteins have been shown to affect homologous interactions in the amyloid self-assembly process. Here, we applied two β-hairpin-forming peptides derived from immunoglobulin light chains as models to test how heterologous β-hairpins modulate the fibril formation of Parkinson's disease-associated protein α-synuclein (αSyn). The peptides SMAhp and LENhp comprise β-strands C and C' of the κ4 antibodies SMA and LEN, which are associated with light chain amyloidosis and multiple myeloma, respectively. SMAhp and LENhp bind with high affinity to the β-hairpin-binding protein β-wrapin AS10 according to isothermal titration calorimetry and NMR spectroscopy. The addition of SMAhp and LENhp affects the kinetics of αSyn aggregation monitored by Thioflavin T (ThT) fluorescence, with the effect depending on assay conditions, salt concentration, and the applied β-hairpin peptide. In the absence of agitation, substoichiometric concentrations of the hairpin peptides strongly reduce the lag time of αSyn aggregation, suggesting that they support the nucleation of αSyn amyloid fibrils. The effect is also observed for the aggregation of αSyn fragments lacking the N-terminus or the C-terminus, indicating that the promotion of nucleation involves the interaction of hairpin peptides with the hydrophobic non-amyloid-β component (NAC) region.

SEVI Inhibits Aβ Amyloid Aggregation by Capping the β-Sheet Elongation Edges

Inhibiting the aggregation of amyloid peptides with endogenous peptides has broad interest due to their intrinsically high biocompatibility and low immunogenicity. Here, we investigated the inhibition mechanism of the prostatic acidic phosphatase fragment SEVI (semen-derived enhancer of viral infection) against Aβ42 fibrillization using atomistic discrete molecular dynamic simulations. Our result revealed that SEVI was intrinsically disordered with dynamic formation of residual helices. With a high positive net charge, the self-aggregation tendency of SEVI was weak. Aβ42 had a strong aggregation propensity by readily self-assembling into β-sheet-rich aggregates. SEVI preferred to interact with Aβ42, rather than SEVI themselves. In the heteroaggregates, Aβ42 mainly adopted β-sheets buried inside and capped by SEVI in the outer layer. SEVI could bind to various Aβ aggregation species─including monomers, dimers, and proto-fibrils─by capping the exposed β-sheet elongation edges. The aggregation processes Aβ42 from the formation of oligomers to conformational nucleation into fibrils and fibril growth should be inhibited as their β-sheet elongation edges are being occupied by the highly charged SEVI. Overall, our computational study uncovered the molecular mechanism of experimentally observed inhibition of SEVI against Aβ42 aggregation, providing novel insights into the development of therapeutic strategies against Alzheimer's disease.

Governing dynamics and preferential binding of the AXH domain influence the aggregation pathway of Ataxin-1

The present state of understanding the mechanism of Spinocerebellar Ataxia-1, a fatal neurodegenerative disease linked to the protein Ataxin-1 (ATXN1), is baffled by a set of self-contradictory, and hence, inconclusive observations. This fallacy poses a bottleneck to the effective designing of curable drugs as the field is currently missing the specific druggable site. To understand the fundamentals of pathogenesis, we tried to decipher the intricacies of the extremely complicated landscape by targeting the relevant species that supposedly dictate the structure-function paradigm. The atomic-level description and characterization of the dynamism of the systems reveal the existence of structural polymorphism in all the leading stakeholders of the overall system. The very existence of conformational heterogeneity in every species creates numerous possible combinations of favorable interactions because of the variability in segmental cross-talks and hence claims its role in the choice of routes between functional activity and dysfunctional disease-causing aggregation. Despite this emergent configurational diversity, there is a common mode of operative intermolecular forces that dictates the extent of stability of all the multimeric complexes due to the localized population of a specific type of residue. The present research proposes a dynamic switch mechanism between aggregability and functional activity, based on the logical interpretation of the estimated variables, which is practically dictated by the effective concentration of the interacting species involved in the cell.

Structural properties of Aβ (1-40) peptide in protonation stage of one, two, and three: New insights from the histidine protonation behaviors

Structural properties and aggregation tendency can be significantly influenced by histidine behaviors (histidine on N^δ-H state is defined as δ, likewise, N^ε-H: ε and both N^δ-H and N^ε-H: p). In current study, we investigated structural properties of Aβ(1-40) peptide during protonation evolution stage of one, two, and three by total 19 independent replica exchange molecular dynamics simulations using implicit solvent. Our results show that any kind of protonated state will promote β-sheet structure formation in comparison with deprotonated (εεε). With increase in number of protonation, the lowest β-sheet content increased. The highest averaged β-sheet structure content was detected in (δpδ) (46.0 %), (εpp) (36.8 %), and (ppp) (16.0 %) in each protonation stage. With three β-strand structures, (δpδ) shows more stable features and high hydrophobic properties. Further analysis confirmed that H13 and H14 are more important than H6. Specifically, H13 and H14 have a synergistic effect for structural formations by controlling H-bond networks in H13(p) with V39/V40 and H14(p/δ) with G37/G38. Finally, the Pearson correlation coefficient results confirmed that experimental result (ref. 44) is corresponding to our (εpp) system. Our current study will be conducive to understanding the effects of the histidine behaviors, it provides new insights for exploration protein folding and misfolding processes.

Insights into the Mixture of Aβ24 and Aβ42 Peptides from Atomistic Simulations

Amyloid-β (Aβ) oligomers play a central role in Alzheimer's disease (AD). Plaques of AD patients consist of Aβ40 and Aβ42 peptides and truncated Aβ peptides. The Aβ24 peptide, identified in human AD brains, was found to impair Aβ42 clearance through the brain-blood barrier. The Aβ24 peptide was also shown to reduce Aβ42 aggregation kinetics in pure buffer, but the underlying mechanism is unknown at atomistic level. In this study, we explored the conformational ensemble of the equimolar mixture of Aβ24 and Aβ42 by replica exchange molecular dynamics simulations and compared it to our previous results on the pure Aβ42 dimer. Our simulations demonstrate that the truncation at residue 24 changes the secondary, tertiary, and quaternary structures of the dimer, offering an explanation of the slower aggregation kinetics of the mixture.

Key Residue for Aggregation of Amyloid-β Peptides

It is known that oligomers of amyloid-β (Aβ) peptide are associated with Alzheimer's disease. Aβ has two isoforms: Aβ40 and Aβ42. Although the difference between Aβ40 and Aβ42 is only two additional C-terminal residues, Aβ42 aggregates much faster than Aβ40. It is unknown what role the C-terminal two residues play in accelerating aggregation. Since Aβ42 is more toxic than Aβ40, its oligomerization process needs to be clarified. Moreover, clarifying the differences between the oligomerization processes of Aβ40 and Aβ42 is essential to elucidate the key factors of oligomerization. Therefore, to investigate the dimerization process, which is the early oligomerization process, Hamiltonian replica-permutation molecular dynamics simulations were performed for Aβ40 and Aβ42. We identified a key residue, Arg5, for the Aβ42 dimerization. The two additional residues in Aβ42 allow the C-terminus to form contact with Arg5 because of the electrostatic attraction between them, and this contact stabilizes the β-hairpin. This β-hairpin promotes dimer formation through the intermolecular β-bridges. Thus, we examined the effects of amino acid substitutions of Arg5, thereby confirming that the mutations remarkably suppressed the aggregation of Aβ42. Moreover, the mutations of Arg5 suppressed the Aβ40 aggregation. It was found by analyzing the simulations that Arg5 is important for Aβ40 to form intermolecular contacts. Thus, it was clarified that the role of Arg5 in the oligomerization process varies due to the two additional C-terminal residues.

Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins

For intrinsically disordered proteins (IDPs), a pressing question is how sequence codes for function. Dynamics serves as a crucial link, reminiscent of the role of structure in sequence-function relations of structured proteins. To define general rules governing sequence-dependent backbone dynamics, we carried out long molecular dynamics simulations of eight IDPs. Blocks of residues exhibiting large amplitudes in slow dynamics are rigidified by local inter-residue interactions or secondary structures. A long region or an entire IDP can be slowed down by long-range contacts or secondary-structure packing. On the other hand, glycines promote fast dynamics and either demarcate rigid blocks or facilitate multiple modes of local and long-range inter-residue interactions. The sequence-dependent backbone dynamics endows IDPs with versatile response to binding partners, with some blocks recalcitrant while others readily adapting to intermolecular interactions.

Characterization of the Conformations of Amyloid Beta 42 in Solution That May Mediate Its Initial Hydrophobic Aggregation

Intrinsically disordered peptides, such as amyloid β42 (Aβ42), lack a well-defined structure in solution. Aβ42 can undergo abnormal aggregation and amyloidogenesis in the brain, forming fibrillar plaques, a hallmark of Alzheimer's disease. The insoluble fibrillar forms of Aβ42 exhibit well-defined, cross β-sheet structures at the molecular level and are less toxic than the soluble, intermediate disordered oligomeric forms. However, the mechanism of initial interaction of monomers and subsequent oligomerization is not well understood. The structural disorder of Aβ42 adds to the challenges of determining the structural properties of its monomers, making it difficult to understand the underlying molecular mechanism of pathogenic aggregation. Certain regions of Aβ42 are known to exhibit helical propensity in different physiological conditions. NMR spectroscopy has shown that the Aβ42 monomer at lower pH can adopt an α-helical conformation and as the pH is increased, the peptide switches to β-sheet conformation and aggregation occurs. CD spectroscopy studies of aggregation have shown the presence of an initial spike in the amount of α-helical content at the start of aggregation. Such an increase in α-helical content suggests a mechanism wherein the peptide can expose critical non-polar residues for interaction, leading to hydrophobic aggregation with other interacting peptides. We have used molecular dynamics simulations to characterize in detail the conformational landscape of monomeric Aβ42 in solution to identify molecular properties that may mediate the early stages of oligomerization. We hypothesized that conformations with α-helical structure have a higher probability of initiating aggregation because they increase the hydrophobicity of the peptide. Although random coil conformations were found to be the most dominant, as expected, α-helical conformations are thermodynamically accessible, more so than β-sheet conformations. Importantly, for the first time α-helical conformations are observed to increase the exposure of aromatic and hydrophobic residues to the aqueous solvent, favoring their hydrophobically driven interaction with other monomers to initiate aggregation. These findings constitute a first step toward characterizing the mechanism of formation of disordered, low-order oligomers of Aβ42.

Scan-Find-Scan-Model: Discrete Site-Targeted Suppressor Design Strategy for Amyloid-β

Alzheimer's disease is undoubtedly the most well-studied neurodegenerative disease. Consequently, the amyloid-β (Aβ) protein ranks at the top in terms of getting attention from the scientific community for structural property-based characterization. Even after decades of extensive research, there is existing volatility in terms of understanding and hence the effective tackling procedures against the disease that arises due to the lack of knowledge of both specific target- and site-specific drugs. Here, we develop a multidimensional approach based on the characterization of the common static-dynamic-thermodynamic trait of the monomeric protein, which efficiently identifies a small target sequence that contains an inherent tendency to misfold and consequently aggregate. The robustness of the identification of the target sequence comes with an abundance of a priori knowledge about the length and sequence of the target and hence guides toward effective designing of the target-specific drug with a very low probability of bottleneck and failure. Based on the target sequence information, we further identified a specific mutant that showed the maximum potential to act as a destabilizer of the monomeric protein as well as enormous success as an aggregation suppressor. We eventually tested the drug efficacy by estimating the extent of modulation of binding affinity existing within the fibrillar form of the Aβ protein due to a single-point mutation and hence provided a proof of concept of the entire protocol.

Multi-eGO: An in silico lens to look into protein aggregation kinetics at atomic resolution

Protein aggregation into amyloid fibrils is the archetype of aberrant biomolecular self-assembly processes, with more than 50 associated diseases that are mostly uncurable. Understanding aggregation mechanisms is thus of fundamental importance and goes in parallel with the structural characterization of the transient oligomers formed during the process. Oligomers have been proven elusive to high-resolution structural techniques, while the large sizes and long time scales, typical of aggregation processes, have limited the use of computational methods to date. To surmount these limitations, we here present multi-eGO, an atomistic, hybrid structure-based model which, leveraging the knowledge of monomers conformational dynamics and of fibril structures, efficiently captures the essential structural and kinetics aspects of protein aggregation. Multi-eGO molecular dynamics simulations can describe the aggregation kinetics of thousands of monomers. The concentration dependence of the simulated kinetics, as well as the structural features of the resulting fibrils, are in qualitative agreement with in vitro experiments carried out on an amyloidogenic peptide from Transthyretin, a protein responsible for one of the most common cardiac amyloidoses. Multi-eGO simulations allow the formation of primary nuclei in a sea of transient lower-order oligomers to be observed over time and at atomic resolution, following their growth and the subsequent secondary nucleation events, until the maturation of multiple fibrils is achieved. Multi-eGO, combined with the many experimental techniques deployed to study protein aggregation, can provide the structural basis needed to advance the design of molecules targeting amyloidogenic diseases.

A Small Molecule Stabilizes the Disordered Native State of the Alzheimer's Aβ Peptide

The stabilization of native states of proteins is a powerful drug discovery strategy. It is still unclear, however, whether this approach can be applied to intrinsically disordered proteins. Here, we report a small molecule that stabilizes the native state of the Aβ42 peptide, an intrinsically disordered protein fragment associated with Alzheimer's disease. We show that this stabilization takes place by a disordered binding mechanism, in which both the small molecule and the Aβ42 peptide remain disordered. This disordered binding mechanism involves enthalpically favorable local π-stacking interactions coupled with entropically advantageous global effects. These results indicate that small molecules can stabilize disordered proteins in their native states through transient non-specific interactions that provide enthalpic gain while simultaneously increasing the conformational entropy of the proteins.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Single-molecule fluorescence imaging and deep learning reveal highly heterogeneous aggregation of amyloid-β 42

There are various diseases caused by protein aggregation such as Alzheimer’s, Parkinson’s, and Huntington’s diseases. From the diversity in the fibril structure, aggregation is expected to occur via heterogeneous pathways. However, characterization of this heterogeneity is extremely difficult because it requires following individual fibril formation in a mixture from early oligomerization stages. In this work, we investigated aggregation of the 42-residue isoform of amyloid β (Aβ42) using single-molecule fluorescence imaging and deep learning. We could track the growth of individual fibrils, which allows for a quantitative description of heterogeneous fibril formation and discovery of a new fibril nucleation mechanism. Further characterization of heterogeneity involving Aβ42 will be important for better understanding the disease mechanism.

Polymorphism in the structure of amyloid fibrils suggests the existence of many different assembly pathways. Characterization of this heterogeneity is the key to understanding the aggregation mechanism and toxicity, but in practice it is extremely difficult to probe individual aggregation pathways in a mixture. Here, we present development of a method combining single-molecule fluorescence lifetime imaging and deep learning for monitoring individual fibril formation in real time and their high-throughput analysis. A deep neural network (FNet) separates an image of highly overlapping fibrils into single fibril images, which allows for tracking the growth and changes in characteristics of individual fibrils. Using this method, we investigated aggregation of the 42-residue amyloid-β peptide (Aβ42). We demonstrate that highly heterogeneous fibril formation can be quantitatively characterized in terms of the number of cross-β subunits, elongation speed, growth polarity, and conformation of fibrils. Tracking individual fibril formation and growth also leads to the discovery of a general nucleation mechanism (termed heterogeneous secondary nucleation), where a fibril is formed on the surface of an oligomer with a different structure. Our development will be broadly applicable to characterization of heterogeneous aggregation processes of other proteins.

A novel non‑selective atypical PKC agonist could protect neuronal cell line from Aβ‑oligomer induced toxicity by suppressing Aβ generation

Atypical protein kinase C (aPKCs) serve key functions in embryonic development by regulating apical-basal polarity. Previous studies have shed light on their roles during adulthood, especially in the development of Alzheimer's disease (AD). Although the crystal structure of PKCι has been resolved, an agonist of aPKCs remains to be discovered. In the present study, by using the Discovery Studio program and LibDock methodology, a small molecule library (K66-X4436 KINA Set) of compounds were screened for potential binding to PKCι. Subsequently, the computational docking results were validated using affinity selection-mass spectrometry, before in vitro kinase activity was used to determine the function of the hit compounds. A cell-based model assay that can mimic the pathology of AD was then established and used to assess the function of these hit compounds. As a result, the aPKC agonist Z640 was identified, which could bind to PKCι in silico, in vitro and in this cell-based model. Z640 was further confirmed as a non-selective aPKC agonist that can activate the kinase activity of both PKCι and PKCζ. In the cell-based assay, Z640 was found to protect neuronal cell lines from amyloid-β (Aβ) oligomer-induced cell death by reducing reactive oxygen species production and restore mitochondrial function. In addition, Z640 could reduce Aβ40 generation in a dose-dependent manner and shift amyloid precursor protein processing towards the non-amyloid pathway. To conclude, the present study is the first, to the best of the authors' knowledge to identify an aPKC agonist by combining computer-assisted drug discovery and cell-based assays. The present study also revealed that aPKC agonists have therapeutic potential for the treatment of AD.

Epitope alteration by small molecules and applications in drug discovery

Small molecules and antibodies are normally considered separately in drug discovery, except in the case of covalent conjugates. We unexpectedly discovered several small molecules that could inhibit or enhance antibody–epitope interactions which opens new possibilities in drug discovery and therapeutic modulation of auto-antibodies. We first discovered a small molecule, CRANAD-17, that enhanced the binding of an antibody to amyloid beta (Aβ), one of the major hallmarks of Alzheimer's disease, by stable triplex formation. Next, we found several small molecules that altered antibody–epitope interactions of tau and PD-L1 proteins, demonstrating the generality of this phenomenon. We report a new screening technology for ligand discovery, screening platform based on epitope alteration for drug discovery (SPEED), which is label-free for both the antibody and small molecule. SPEED, applied to an Aβ antibody, led to the discovery of a small molecule, GNF5837, that inhibits Aβ aggregation and another, obatoclax, that binds Aβ plaques and can serve as a fluorescent reporter in brain slices of AD mice. We also found a small molecule that altered the binding between Aβ and auto-antibodies from AD patient serum. SPEED reveals the sensitivity of antibody–epitope interactions to perturbation by small molecules and will have multiple applications in biotechnology and drug discovery.

A screening platform based on epitope alteration for drug discovery (SPEED).

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

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

Functional Gallic Acid-Based Dendrimers as Synthetic Nanotools to Remodel Amyloid-β-42 into Noncytotoxic Forms

The self-assembly of amyloid-β (Aβ) generates cytotoxic oligomers linked to the onset and progression of Alzheimer's disease (AD). As many fundamental molecular pathways that control Aβ aggregation are yet to be unraveled, an important strategy to control Aβ cytotoxicity is the development of bioactive synthetic nanotools capable of interacting with the heterogeneous ensemble of Aβ species and remodel them into noncytotoxic forms. Herein, the synthesis of nanosized, functional gallic acid (Ga)-based dendrimers with a precise number of Ga at their surface is described. It is shown that these Ga-terminated dendrimers interact by H-bonding with monomeric/oligomeric Aβ species at their Glu, Ala, and Asp residues, promoting their remodeling into noncytotoxic aggregates in a process controlled by the Ga units. The multivalent presentation of Ga on the dendrimer surface enhances their ability to interact with Aβ, inhibiting the primary and secondary nucleation of Aβ fibrillization and disrupting the Aβ preformed fibrils.

Data Mining of Molecular Simulations Suggest Key Amino Acid Residues for Aggregation, Signaling and Drug Action

Alzheimer's disease, the most common form of dementia, currently has no cure. There are only temporary treatments that reduce symptoms and the progression of the disease. Alzheimer's disease is characterized by the prevalence of plaques of aggregated amyloid β (Aβ) peptide. Recent treatments to prevent plaque formation have provided little to relieve disease symptoms. Although there have been numerous molecular simulation studies on the mechanisms of Aβ aggregation, the signaling role has been less studied. In this study, a total of over 38,000 simulated structures, generated from molecular dynamics (MD) simulations, exploring different conformations of the Aβ42 mutants and wild-type peptides were used to examine the relationship between Aβ torsion angles and disease measures. Unique methods characterized the data set and pinpointed residues that were associated in aggregation and others associated with signaling. Machine learning techniques were applied to characterize the molecular simulation data and classify how much each residue influenced the predicted variant of Alzheimer's Disease. Orange3 data mining software provided the ability to use these techniques to generate tables and rank the data. The test and score module coupled with the confusion matrix module analyzed data with calculations of specificity and sensitivity. These methods evaluating frequency and rank allowed us to analyze and predict important residues associated with different phenotypic measures. This research has the potential to help understand which specific residues of Aβ should be targeted for drug development.

Met35 Oxidation Hinders Aβ25-35 Peptide Aggregation within the Dimyristoylphosphatidylcholine Bilayer

Using all-atom explicit solvent replica exchange molecular dynamics simulations, we studied the aggregation of oxidized (ox) Aβ25-35 peptides into dimers mediated by the zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer. By comparing oxAβ25-35 aggregation with that observed for reduced and phosphorylated Aβ25-35 peptides, we elucidated plausible impact of post-translational modifications on cytotoxicity of Aβ peptides involved in Alzheimer's disease. We found that Met35 oxidation reduces helical propensity in oxAβ25-35 peptides bound to the lipid bilayer and enhances backbone fluctuations. These factors destabilize the wild-type head-to-tail dimer interface and lower the aggregation propensity. Met35 oxidation diversifies aggregation pathways by adding monomeric species to the bound conformational ensemble. The oxAβ25-35 dimer becomes partially expelled from the DMPC bilayer and as a result inflicts limited disruption to the bilayer structure compared to wild-type Aβ25-35. Interestingly, the effect of Ser26 phosphorylation is largely opposite, as it preserves the wild-type head-to-tail aggregation interface and strengthens, not weakens, aggregation propensity. The differing effects can be attributed to the sequence locations of these post-translational modifications, since in contrast to Ser26 phosphorylation, Met35 oxidation directly affects the wild-type C-terminal aggregation interface. A comparison with experimental data is provided.

Role of the English (H6R) Mutation on the Structural Properties of Aβ40 and Aβ42 Owing to the Histidine Tautomeric Effect

As an intrinsic origin cause, histidine behaviors play a critical role in protein misfolding processes. Generally, the English (H6R) mutation will disrupt H6 interactions. However, the structural properties of Aβ40 H6R and Aβ42 H6R under the complex influence of a histidine tautomeric effect and an H6R mutation remain unclear. Therefore, we performed a replica exchange molecular dynamics simulation to unveil such structural properties. Our result showed that the H6R substitute could promote the generation of β-sheet structures in comparison to the wild type. Three β-strand structure properties were observed in Aβ40 (rδδ), Aβ42 (rεε), Aβ42 (rεδ), and Aβ42 (rδδ) with β-sheet contents of 47.5%, 37.2%, 46.9%, and 38.6%, respectively, and the dominant conformational properties of Aβ40 (rδδ), Aβ42 (rεε), Aβ42 (rεδ), and Aβ42 (rδδ) had top conformational states of 86.0%, 73.2%, 67.0%, and 56.5%, respectively. Further analysis confirmed that R6 had different mechanisms for controlling the conformational features in Aβ40 H6R and Aβ42 H6R. In the Aβ40 systems, H14 H-bond networks played a critical role in controlling the structural properties. However, in the Aβ42 systems, R6 was more important because it was directly involved in the β-strand formation and maintained the β-sheet between the N-terminus and the central hydrophobic core region. Our current study helps to elucidate the histidine tautomeric behaviors in H6R mutations, which will present opportunities to understand the correlation between with/without H6 and the Aβ40/Aβ42 H6R misfolding mechanisms.

Spontaneous Formation of β-sheet Nano-barrels during the Early Aggregation of Alzheimer's Amyloid Beta

Soluble low-molecular-weight oligomers formed during the early aggregation of amyloid peptides have been hypothesized as a major toxic species of amyloidogenesis. Herein, we performed the first synergic in silico, in vitro and in vivo validations of the structure, dynamics and toxicity of Aβ42 oligomers. Aβ peptides readily assembled into β-rich oligomers comprised of extended β-hairpins and β-strands. Nanosized β-barrels were observed with certainty with simulations, transmission electron microscopy and Fourier transform infrared spectroscopy, corroborated by immunohistochemistry, cell viability, apoptosis, inflammation, autophagy and animal behavior assays. Secondary and tertiary structural proprieties of these oligomers, such as the sequence regions with high β-sheet propensities and inter-residue contact frequency patterns, were similar to the properties known for Aβ fibrils. The unambiguous spontaneous formation of β-barrels in the early aggregation of Aβ42 supports their roles as the common toxic intermediates in Alzheimer's pathobiology and a target for Alzheimer's therapeutics.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

Cholesterol Molecules Alter the Energy Landscape of Small Aβ1-42 Oligomers

Small amyloid-β (Aβ) oligomers are believed to be key pathogenic species in Alzheimer's disease (AD). One suggested toxicity mechanism is the detergent model where oligomers remove lipid molecules from the bilayer. Senile plaques of AD patients also accumulate a 1:1 ratio of cholesterol/Aβ. What are the dominant structures of small Aβ42 oligomers with cholesterol molecules in aqueous solution? Here, we answer this question by performing atomistic replica exchange molecular dynamics simulations of Aβ42 dimers and trimers. Our simulations demonstrate that the interactions with cholesterol molecules change completely the energy landscape of small Aβ42 oligomers. This result shows that simulations in the bulk solution cannot recapitulate aggregation in the brain extracellular space.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Secondary structure dependence of amyloid-β(1-40) on simulation techniques and force field parameters

Our recent studies revealed that none of the selected widely used force field parameters and molecular dynamics simulation techniques yield structural properties for the intrinsically disordered α-synuclein that are in agreement with various experiments via testing different force field parameters. Here, we extend our studies on the secondary structure properties of the disordered amyloid-β(1-40) peptide in aqueous solution. For these purposes, we conducted extensive replica exchange molecular dynamics simulations and obtained extensive molecular dynamics simulation trajectories from David E. Shaw group. Specifically, these molecular dynamics simulations were conducted using various force field parameters and obtained results are compared to our replica exchange molecular dynamics simulations and experiments. In this study, we calculated the secondary structure abundances and radius of gyration values for amyloid-β(1-40) that were simulated using varying force field parameter sets and different simulation techniques. In addition, the intrinsic disorder propensity, as well as sequence-based secondary structure predisposition of amyloid-β(1-40) and compared the findings with the results obtained from molecular simulations using various force field parameters and different simulation techniques. Our studies clearly show that the epitope region identification of amyloid-β(1-40) depends on the chosen simulation technique and chosen force field parameters. Based on comparison with experiments, we find that best computational results in agreement with experiments are obtained using the a99sb*-ildn, charmm36m, and a99sb-disp parameters for the amyloid-β(1-40) peptide in molecular dynamics simulations without parallel tempering or via replica exchange molecular dynamics simulations.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

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

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

β-Turn exchanges in the α-synuclein segment 44-TKEG-47 reveal high sequence fidelity requirements of amyloid fibril elongation

The folding of turns and β-hairpins has been implicated in amyloid formation, with diverse potential consequences such as promotion or inhibition of fibril nucleation, fibril elongation, or off-pathway oligomer formation. In the Parkinson's disease-associated protein α-synuclein (αS), a β-hairpin comprised of residues 36-56 was detected in complex with an engineered binding protein, with a turn formed by the αS sequence segment 44-TKEG-47. Molecular dynamics simulations revealed extensive populations of transient β-hairpin conformations in this region in free, monomeric αS. Here, we investigated potential effects of turn formation on αS fibril formation by studying the aggregation kinetics of an extensive set of αS variants with between two and four amino acid exchanges in the 44-TKEG-47 segment. The exchanges were chosen to specifically promote formation of β1-, β1'-, or β2'-turns. All variants assembled into amyloid fibrils, with increased β1'- or β2'-turn propensity associated with faster aggregation and increased β1-turn propensity with slower aggregation compared to wild-type (WT) αS. Atomic force microscopy demonstrated that β-turn exchanges altered fibril morphology. In cross-elongation experiments, the turn variants showed a low ability to elongate WT fibril seeds, and, vice versa, WT monomer did not efficiently elongate turn variant fibril seeds. This demonstrates that sequence identity in the turn region is crucial for efficient αS fibril elongation. Elongation experiments of WT fibril seeds in the presence of both WT and turn variant monomers suggest that the turn variants can bind and block WT fibril ends to different degrees, but cannot efficiently convert into the WT fibril structure. Our results indicate that modifications in the 44-TKEG-47 segment strongly affect amyloid assembly by driving αS into alternative fibril morphologies, whose elongation requires high sequence fidelity.

Hydration Thermodynamics of the N-Terminal FAD Mutants of Amyloid-β

The hydration thermodynamics of amyloid-β (Aβ) and its pathogenic familial Alzheimer's disease (FAD) mutants such as A2V, Taiwan (D7H), Tottori (D7N), and English (H6R) and the protective A2T mutant is investigated by a combination of all-atom, explicit water molecular dynamics (MD) simulations and the three-dimensional reference interaction site model (3D-RISM) theory. The change in the hydration free energy on mutation is decomposed into the energetic and entropic components, which comprise electrostatic and nonelectrostatic contributions. An increase in the hydration free energy is observed for A2V, D7H, D7N, and H6R mutations that increase the aggregation propensity of Aβ and lead to an early onset of Alzheimer's disease, while a reverse trend is noted for the protective A2T mutation. An antiphase correlation is found between the change in the hydration energy and the internal energy of Aβ upon mutation. A residue-wise decomposition analysis shows that the change in the hydration free energy of Aβ on mutation is primarily due to the hydration/dehydration of the side-chain atoms of the negatively charged residues. The decrease in the hydration of the negatively charged residues on mutation may decrease the solubility of the mutant, which increases the observed aggregation propensity of the FAD mutants. Results obtained from the theory show an excellent match with the experimentally reported data.

A kinetic ensemble of the Alzheimer's Aβ peptide

The conformational and thermodynamic properties of disordered proteins are commonly described in terms of structural ensembles and free energy landscapes. To provide information on the transition rates between the different states populated by these proteins, it would be desirable to generalize this description to kinetic ensembles. Approaches based on the theory of stochastic processes can be particularly suitable for this purpose. Here, we develop a Markov state model and apply it to determine a kinetic ensemble of Aβ42, a disordered peptide associated with Alzheimer's disease. Through the Google Compute Engine, we generated 315-µs all-atom molecular dynamics trajectories. Using a probabilistic-based definition of conformational states in a neural network approach, we found that Aβ42 is characterized by inter-state transitions on the microsecond timescale, exhibiting only fully unfolded or short-lived, partially folded states. Our results illustrate how kinetic ensembles provide effective information about the structure, thermodynamics and kinetics of disordered proteins.

A lowly populated, transient β-sheet structure in monomeric Aβ1-42 identified by multinuclear NMR of chemical denaturation

Chemical denaturation is a well-established approach for probing the equilibrium between folded and unfolded states of proteins. We demonstrate applicability of this method to the detection of a small population of a transiently folded structural element in a system that is often considered to be intrinsically fully disordered. The ¹H^N, ¹⁵N, ¹³C^α, and ¹³C' chemical shifts of Aβ^1-40 and Aβ^1-42 peptides and their M35-oxidized variants were monitored as a function of urea concentration and compared to analogous urea titrations of synthetic pentapeptides of homologous sequence. Fitting of the chemical shift titrations yields a 10 ± 1% population for a structured element at the C-terminus of Aβ^1-42 that folds with a cooperativity of m = 0.06 kcal/mol·M. The fit also yields the chemical shifts of the folded state and, using a database search, for Aβ^1-42 these shifts identified an antiparallel intramolecular β-sheet for residues I32-A42, linked by a type I' β-turn at G37 and G38. The structure is destabilized by oxidation of M35. Paramagnetic relaxation rates and two previously reported weak, medium-range NOE interactions are consistent with this transient β-sheet. Introduction of the requisite A42C mutation and tagging with MTSL resulted in a small stabilization of this β-sheet. Chemical shift analysis suggests a C-terminal β-sheet may be present in Aβ^1-40 too, but the turn type at G37 is not type I'. The approach to derive Transient Structure from chemical Denaturation by NMR (TSD-NMR), demonstrated here for Aβ peptides, provides a sensitive tool for identifying the presence of lowly populated, transiently ordered elements in proteins that are considered to be intrinsically disordered, and permits extraction of structural data for such elements.

Disclosing the Interaction of Gold Nanoparticles with Aβ(1-40) Monomers through Replica Exchange Molecular Dynamics Simulations

Amyloid-β aggregation is one of the principal causes of amyloidogenic diseases that lead to the loss of neuronal cells and to cognitive impairments. The use of gold nanoparticles treating amyloidogenic diseases is a promising approach, because the chemistry of the gold surface can be tuned in order to have a specific binding, obtaining effective tools to control the aggregation. In this paper, we show, by means of Replica Exchange Solute Tempering Molecular Simulations, how electrostatic interactions drive the absorption of Amyloid-β monomers onto citrates-capped gold nanoparticles. Importantly, upon binding, amyloid monomers show a reduced propensity in forming β-sheets secondary structures that are characteristics of mature amyloid fibrils.

Predicting Genetic Variation Severity Using Machine Learning to Interpret Molecular Simulations

Distinct missense mutations in a specific gene have been associated with different diseases as well as differing severity of a disease. Current computational methods predict the potential pathogenicity of a missense variant but fail to differentiate between separate disease or severity phenotypes. We have developed a method to overcome this limitation by applying machine learning to features extracted from molecular dynamics simulations, creating a way to predict the effect of novel genetic variants in causing a disease, drug resistance, or another specific trait. As an example, we have applied this novel approach to variants in calmodulin associated with two distinct arrhythmias as well as two different neurodegenerative diseases caused by variants in amyloid-β peptide. The new method successfully predicts the specific disease caused by a gene variant and ranks its severity with more accuracy than existing methods. We call this method molecular dynamics phenotype prediction model.

Emergence of Barrel Motif in Amyloid-β Trimer: A Computational Study

Amyloid-β (Aβ) peptides form assemblies that are pathological hallmarks of Alzheimer's disease. Aβ oligomers are soluble, mobile, and toxic forms of the peptide that act in the extracellular space before assembling into protofibrils and fibrils. Therefore, oligomers play an important role in the mechanism of Alzheimer's disease. Since it is difficult to determine by experiment the atomic structures of oligomers, which accumulate fast and are polymorphic, computer simulation is a useful tool to investigate elusive oligomers' structures. In this work, we report extended all-atom molecular dynamics simulations, both canonical and replica exchange, of Aβ(1-42) trimer starting from two different initial conformations: (i) the pose produced by the best docking of a monomer aside of a dimer (simulation 1), representing oligomers freshly formed by assembling monomers, and (ii) a configuration extracted from an experimental mature fibril structure (simulation 2), representing settled oligomers in equilibrium with extended fibrils. We showed that in simulation 1, regions with small β-barrels are populated, indicating the chance of spontaneous formation of domains resembling channel-like structures. These structural domains are alternative to those more representative of mature fibrils (simulation 2), the latter showing a stable bundle of C-termini that is not sampled in simulation 1. Moreover, trimer of Aβ(1-42) can form internal pores that are large enough to be accessed by water molecules and Ca2+ ions.

Allosteric inhibition explained through conformational ensembles sampling distinct "mixed" states

Allosteric modulation provides an effective avenue for selective and potent enzyme inhibition. Here, we summarize and critically discuss recent advances on the mechanisms of allosteric partial agonists for three representative signalling enzymes activated by cyclic nucleotides: the cAMP-dependent protein kinase (PKA), the cGMP-dependent protein kinase (PKG), and the exchange protein activated by cAMP (EPAC). The comparative analysis of partial agonism in PKA, PKG and EPAC reveals a common emerging theme, i.e. the sampling of distinct “mixed” conformational states, either within a single domain or between distinct domains. Here, we show how such “mixed” states play a crucial role in explaining the observed functional response, i.e. partial agonism and allosteric pluripotency, as well as in maximizing inhibition while minimizing potency losses. In addition, by combining Nuclear Magnetic Resonance (NMR), Molecular Dynamics (MD) simulations and Ensemble Allosteric Modeling (EAM), we also show how to map the free-energy landscape of conformational ensembles containing “mixed” states. By discussing selected case studies, we illustrate how MD simulations and EAM complement NMR to quantitatively relate protein dynamics to function. The resulting NMR- and MD-based EAMs are anticipated to inform not only the design of new generations of highly selective allosteric inhibitors, but also the choice of multidrug combinations.

αB-Crystallin Chaperone Inhibits Aβ Aggregation by Capping the β-Sheet-Rich Oligomers and Fibrils

Inhibiting the cytotoxicity of amyloid aggregation by endogenous proteins is a promising strategy against degenerative amyloid diseases due to their intrinsically high biocompatibility and low immunogenicity. In this study, we investigated the inhibition mechanism of the structured core region of αB-crystallin (αBC) against Aβ fibrillization using discrete molecular dynamics simulations. Our computational results recapitulated the experimentally observed Aβ binding sites in αBC and suggested that αBC could bind to various Aβ aggregate species during the aggregation process-including monomers, dimers, and likely other high molecular weight oligomers, protofibrils, and fibrils-by capping the exposed β-sheet elongation surfaces. Thus, the nucleation of Aβ oligomers into fibrils and the fibril growth could be inhibited. Mechanistic insights obtained from our systematic computational studies may aid in the development of novel therapeutic strategies to modulate the aggregation of pathological, amyloidogenic protein in degenerative diseases.