Monomeric Aβ(1-40) and Aβ(1-42) Peptides in Solution Adopt Very Similar Ramachandran Map Distributions That Closely Resemble Random Coil

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.

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.

NMR spectroscopy captures the essential role of dynamics in regulating biomolecular function

Biomolecules are in constant motion. To understand how they function, and why malfunctions can cause disease, it is necessary to describe their three-dimensional structures in terms of dynamic conformational ensembles. Here, we demonstrate how nuclear magnetic resonance (NMR) spectroscopy provides an essential, dynamic view of structural biology that captures biomolecular motions at atomic resolution. We focus on examples that emphasize the diversity of biomolecules and biochemical applications that are amenable to NMR, such as elucidating functional dynamics in large molecular machines, characterizing transient conformations implicated in the onset of disease, and obtaining atomic-level descriptions of intrinsically disordered regions that make weak interactions involved in liquid-liquid phase separation. Finally, we discuss the pivotal role that NMR has played in driving forward our understanding of the biomolecular dynamics-function paradigm.

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.

Systematic Differences between Current Molecular Dynamics Force Fields To Represent Local Properties of Intrinsically Disordered Proteins

The prevalence of intrinsically disordered proteins (IDPs) and protein regions in structural biology has prompted the recent development of molecular dynamics (MD) force fields for the more realistic representations of such systems. Using experimental nuclear magnetic resonance backbone scalar ³J-coupling constants of the intrinsically disordered proteins α-synuclein and amyloid-β in their native aqueous environment as a metric, we compare the performance of four recent MD force fields, namely, AMBER ff14SB, CHARMM C36m, AMBER ff99SB-disp, and AMBER ff99SBnmr2, by partitioning the polypeptides into an overlapping series of heptapeptides for which a cumulative total of 276 μs MD simulations were performed. The results show substantial differences between the different force fields at the individual residue level. Except for ff99SBnmr2, the force fields systematically underestimate the scalar ³J(H^N,Hα)-couplings due to an underrepresentation of β-conformations and an overrepresentation of either α- or PP_II conformations. The study demonstrates that the incorporation of coil library information in modern MD force fields, as shown here for ff99SBnmr2, provides substantially improved performance and more realistic sampling of the local backbone dihedral angles of IDPs as reflected by the good accuracy of the computed scalar ³J(H^N,Hα)-couplings with less than 0.5 Hz error. Such force fields will enable a better understanding of how structural dynamics and thermodynamics influence the IDP function. Although the methodology based on heptapeptides used here does not allow the assessment of potential intramolecular long-range interactions, its computational affordability permits well-converged simulations that can be easily parallelized. This should make the quantitative validation of intrinsic disorder observed in MD simulations of polypeptides with experimental scalar J-couplings widely applicable.

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.

An α-helix mimetic oligopyridylamide, ADH-31, modulates Aβ42 monomer aggregation and destabilizes protofibril structures: insights from molecular dynamics simulations

Alzheimer's disease (AD), an epidemic growing worldwide due to no effective medical aid available in the market, is a neurological disorder. AD is known to be directly associated with the toxicity of amyloid-β (Aβ) aggregates. In search of potent inhibitors of Aβ aggregation, Hamilton and co-workers reported an α-helix mimetic, ADH-31, which acts as a powerful antagonist of Aβ42 aggregation. To identify the key interactions between protein-ligand complexes and to gain insights into the inhibitory mechanism of ADH-31 against Aβ42 aggregation, molecular dynamics (MD) simulations were performed in the present study. The MD simulations highlighted that ADH-31 showed distinct binding capabilities with residues spanning from the N-terminal to the central hydrophobic core (CHC) region of Aβ42 and restricted the conformational transition of the helix-rich structure of Aβ42 into another form of secondary structures (coil/turn/β-sheet). Hydrophobic contacts, hydrogen bonding and π-π interaction contribute to the strong binding between ADH-31 and Aβ42 monomer. The Dictionary of Secondary Structure of Proteins (DSSP) analysis highlighted that the probability of helical content increases from 38.5% to 50.2% and the turn content reduces from 14.7% to 6.2% with almost complete loss of the β-sheet structure (4.5% to 0%) in the Aβ42 monomer + ADH-31 complex. The per-residue binding free energy analysis demonstrated that Arg5, Tyr10, His14, Gln15, Lys16, Val18, Phe19 and Lys28 residues of Aβ42 are responsible for the favourable binding free energy in Aβ42 monomer + ADH-31 complex, which is consistent with the 2D HSQC NMR of the Aβ42 monomer that depicted a change in the chemical shift of residues spanning from Glu11 to Phe20 in the presence of ADH-31. The MD simulations highlighted the prevention of sampling of amyloidogenic β-strand conformations in Aβ42 trimer in the presence of ADH-31 as well as the ability of ADH-31 to destabilize Aβ42 trimer and protofibril structures. The lower binding affinity between Aβ42 trimer chains in the presence of ADH-31 highlights the destabilization of the Aβ42 trimer structure. Overall, MD results highlighted that ADH-31 inhibited Aβ42 aggregation by constraining Aβ peptides into helical conformation and destabilized Aβ42 trimer as well as protofibril structures. The present study provides a theoretical insight into the atomic level details of the inhibitory mechanism of ADH-31 against Aβ42 aggregation as well as protofibril destabilization and could be implemented in the structure-based drug design of potent therapeutic agents for AD.

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.

α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.

Structural Characteristics of Monomeric Aβ42 on Fibril in the Early Stage of Secondary Nucleation Process

Amyloid-β (Aβ) aggregates are believed to be one of the main causes of Alzheimer's disease. Aβ peptides form fibrils having cross β-sheet structures mainly through primary nucleation, secondary nucleation, and elongation. In particular, self-catalyzed secondary nucleation is of great interest. Here, we investigate the adsorption of Aβ42 peptides to the Aβ42 fibril to reveal a role of adsorption as a part of secondary nucleation. We performed extensive molecular dynamics simulations based on replica exchange with solute tempering 2 (REST2) to two systems: a monomeric Aβ42 in solution and a complex of an Aβ42 peptide and Aβ42 fibril. Results of our simulations show that the Aβ42 monomer is extended on the fibril. Furthermore, we find that the hairpin structure of the Aβ42 monomer decreases but the helix structure increases by adsorption to the fibril surface. These structural changes are preferable for forming fibril-like aggregates, suggesting that the fibril surface serves as a catalyst in the secondary nucleation process. In addition, the stabilization of the helix structure of the Aβ42 monomer on the fibril indicates that the strategy of a secondary nucleation inhibitor design for Aβ40 can also be used for Aβ42.

Hyaluronan-carnosine conjugates inhibit Aβ aggregation and toxicity

Alzheimer’s disease is the most common neurodegenerative disorder. Finding a pharmacological approach that cures and/or prevents the onset of this devastating disease represents an important challenge for researchers. According to the amyloid cascade hypothesis, increases in extracellular amyloid-β (Aβ) levels give rise to different aggregated species, such as protofibrils, fibrils and oligomers, with oligomers being the more toxic species for cells. Many efforts have recently been focused on multi-target ligands to address the multiple events that occur concurrently with toxic aggregation at the onset of the disease. Moreover, investigating the effect of endogenous compounds or a combination thereof is a promising approach to prevent the side effects of entirely synthetic drugs. In this work, we report the synthesis, structural characterization and Aβ antiaggregant ability of new derivatives of hyaluronic acid (Hy, 200 and 700 kDa) functionalized with carnosine (Car), a multi-functional natural dipeptide. The bioactive substances (HyCar) inhibit the formation of amyloid-type aggregates of Aβ₄₂ more than the parent compounds; this effect is proportional to Car loading. Furthermore, the HyCar derivatives are able to dissolve the amyloid fibrils and to reduce Aβ-induced toxicity in vitro. The enzymatic degradation of Aβ is also affected by the interaction with HyCar.

Experimentally Consistent Simulation of Aβ21-30 Peptides with a Minimal NMR Bias

Misfolded amyloid peptides are neurotoxic molecules associated with Alzheimer's disease. The Aβ_21-30 peptide fragment is a decapeptide fragment of the complete Aβ42 peptide which is a hypothesized cause of Alzheimer's disease via amyloid fibrillogenesis. Aβ_21-30 is investigated here with a combination of NMR (nuclear magnetic resonance) spectroscopy experiments and molecular dynamics simulations with experiment directed simulation (EDS). EDS is a maximum entropy biasing method that augments a molecular dynamics simulation with experimental data (NMR chemical shifts) to improve agreement with experiments and thus accuracy. EDS molecular dynamics shows that the Aβ_21-30 monomer has a β turn stabilized by the following interactions: S26-K28, D23-S26, and D23-K28. NMR, total correlation spectroscopy, and rotating frame Overhauser effect spectroscopy experiments provide independent agreement. Subsequent two- and four-monomer EDS simulations show aggregation. Diffusion coefficients calculated from molecular simulation also agreed with experimentally measured values only after using EDS, providing independent assessment of accuracy. This work demonstrates how accuracy can be improved by directly using experimental data in molecular dynamics of complex processes like self-assembly.

A rationally designed bicyclic peptide remodels Aβ42 aggregation in vitro and reduces its toxicity in a worm model of Alzheimer's disease

Bicyclic peptides have great therapeutic potential since they can bridge the gap between small molecules and antibodies by combining a low molecular weight of about 2 kDa with an antibody-like binding specificity. Here we apply a recently developed in silico rational design strategy to produce a bicyclic peptide to target the C-terminal region (residues 31–42) of the 42-residue form of the amyloid β peptide (Aβ42), a protein fragment whose aggregation into amyloid plaques is linked with Alzheimer’s disease. We show that this bicyclic peptide is able to remodel the aggregation process of Aβ42 in vitro and to reduce its associated toxicity in vivo in a C. elegans worm model expressing Aβ42. These results provide an initial example of a computational approach to design bicyclic peptides to target specific epitopes on disordered proteins.

Curcumin Inhibits the Primary Nucleation of Amyloid-Beta Peptide: A Molecular Dynamics Study

The amyloid plaques are a key hallmark of neurodegenerative diseases such as Alzheimer’s disease and Parkinson’s disease. Amyloidogenesis is a complex long-lasting multiphase process starting with the formation of nuclei of amyloid peptides: a process assigned as a primary nucleation. Curcumin (CU) is a well-known inhibitor of the aggregation of amyloid-beta (Aβ) peptides. Even more, CU is able to disintegrate preformed Aβ firbils and amyloid plaques. Here, we simulate by molecular dynamics the primary nucleation process of 12 Aβ peptides and investigate the effects of CU on the process. We found that CU molecules intercalate among the Aβ chains and bind tightly to them by hydrogen bonds, hydrophobic, π–π, and cation–π interactions. In the presence of CU, the Aβ peptides form a primary nucleus of a bigger size. The peptide chains in the nucleus become less flexible and more disordered, and the number of non-native contacts and hydrogen bonds between them decreases. For comparison, the effects of the weaker Aβ inhibitor ferulic acid (FA) on the primary nucleation are also examined. Our study is in good agreement with the observation that taken regularly, CU is able to prevent or at least delay the onset of neurodegenerative disorders.

Fluorescence-Labeled Amyloid Beta Monomer: A Molecular Dynamical Study

The aggregation process of the Amyloidβ (Aβ) peptide is one of the central questions in Alzheimers’s research. Fluorescence-labeled single-molecule detection is a novel technique concerning the early stage investigation of Aβ aggregation, where the labeling dyes are covalently bound to the Aβ monomer. As the influence of the dye on the conformational space of the Aβ monomer can be significant, its effect on the seeding process is an open question. The applied fluorescent molecule continuously switches between an active (ON) and an inactive (OFF) state, where the latter supports an extra rotational restriction at many commercially available dyes. However, only a few theoretical studies simulated the Aβ monomer in the presence of a dye and none of them considered the difference between the ON and the OFF states. Therefore, we examined the impact of a selected fluorescence dye (Alexa 568) on the conformational space of the monomeric Aβ(1–42) peptide in its ON and OFF state by replica exchange molecular dynamic simulations. Investigations on secondary structure elements as well as dye-peptide contact analysis for the monomers are presented. Experimental and theoretical NMR shifts were contrasted to qualify the calculation protocol and theoretical values of the labeled and the non-labeled peptide were also compared. We found that the first five residues have higher helical propensity in the presence of the dye, and electrostatic properties could strongly affect the connection between the dye and the peptide parts.

Differences in the free energies between the excited states of Aβ40 and Aβ42 monomers encode their aggregation propensities

The early events in the aggregation of the intrinsically disordered peptide, amyloid-β (Aβ), involve transitions from the disordered free energy ground state to assembly-competent states. Are the fingerprints of order found in the amyloid fibrils encoded in the conformations that the monomers access at equilibrium? If so, could the enhanced aggregation rate of Aβ42 compared to Aβ40 be rationalized from the sparsely populated high free energy states of the monomers? Here, we answer these questions in the affirmative using coarse-grained simulations of the self-organized polymer-intrinsically disordered protein (SOP-IDP) model of Aβ40 and Aβ42. Although both the peptides have practically identical ensemble-averaged properties, characteristic of random coils (RCs), the conformational ensembles of the two monomers exhibit sequence-specific heterogeneity. Hierarchical clustering of conformations reveals that both the peptides populate high free energy aggregation-prone ([Formula: see text]) states, which resemble the monomers in the fibril structure. The free energy gap between the ground (RC) and the [Formula: see text] states of Aβ42 peptide is smaller than that for Aβ40. By relating the populations of excited states of the two peptides to the fibril formation time scales using an empirical formula, we explain nearly quantitatively the faster aggregation rate of Aβ42 relative to Aβ40. The [Formula: see text] concept accounts for fibril polymorphs, leading to the prediction that the less stable [Formula: see text] state of Aβ42, encoding for the U-bend fibril, should form earlier than the structure with the S-bend topology, which is in accord with Ostwald's rule rationalizing crystal polymorph formation.

Comparison of Synthetic Neuronal Model Membrane Mimics in Amyloid Aggregation at Atomic Resolution

Alzheimer's disease (AD) is a severe neurodegenerative disorder caused by abnormal accumulation of toxic amyloid plaques of the amyloid-beta (Aβ) or the tau proteins in the brain. The plaque deposition leading to the collapse of the cellular integrity is responsible for a myriad of surface phenomena acting at the neuronal lipid interface. Recent years have witnessed dysfunction of the blood-brain barriers (BBB) associated with AD. Several studies support the idea that BBB acts as a platform for the formation of misfolded Aβ peptide, promoting oligomerization and fibrillation, compromising the overall integrity of the central nervous system. While the amyloid plaque deposition has been known to be responsible for the collapse of the BBB membrane integrity, the causal effect relationship between BBB and Aβ amyloidogenesis remains unclear. In this study, we have used physiologically relevant synthetic model membrane systems to gain atomic insight into the functional aspects of the lipid interface. Here, we have used a minimalist BBB mimic, POPC/POPG/cholesterol/GM1, to compare with the native BBB (total lipid brain extract (TLBE)), to understand the molecular events occurring in the membrane-induced Aβ₄₀ amyloid aggregation. Our study showed that the two membrane models accelerated the Aβ₄₀ aggregation kinetics with differential secondary structural transitions of the peptide. The observed structural transitions are defined by the lipid compositions, which in turn undermines the differences in lipid surface phenomena, leading to peptide induced cellular toxicity in the neuronal membrane.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

How Does the Mono-Triazole Derivative Modulate Aβ42 Aggregation and Disrupt a Protofibril Structure: Insights from Molecular Dynamics Simulations

Clinical studies have identified that abnormal self-assembly of amyloid-β (Aβ) peptide into toxic fibrillar aggregates is associated with the pathology of Alzheimer's disease (AD). The most acceptable therapeutic approach to stop the progression of AD is to inhibit the formation of β-sheet-rich structures. Recently, we designed and evaluated a series of novel mono-triazole derivatives 4(a-x), where compound 4v was identified as the most potent inhibitor of Aβ₄₂ aggregation and disaggregates preformed Aβ₄₂ fibrils significantly. Moreover, 4v strongly averts the Cu²⁺-induced Aβ₄₂ aggregation and disaggregates the preformed Cu²⁺-induced Aβ₄₂ fibrils, halts the generation of reactive oxygen species, and shows neuroprotective effects in SH-SY5Y cells. However, the underlying molecular mechanism of inhibition of Aβ₄₂ aggregation by 4v and disaggregation of preformed Aβ₄₂ fibrils remains obscure. In this work, molecular dynamics (MD) simulations have been performed to explore the conformational ensemble of the Aβ₄₂ monomer and a pentameric protofibril structure of Aβ₄₂ in the presence of 4v. The MD simulations highlighted that 4v binds preferentially at the central hydrophobic core region of the Aβ₄₂ monomer and chains D and E of the Aβ₄₂ protofibril. The dictionary of secondary structure of proteins analysis indicated that 4v retards the conformational conversion of the helix-rich structure of the Aβ₄₂ monomer into the aggregation-prone β-sheet conformation. The binding free energy calculated by the molecular mechanics Poisson-Boltzmann surface area method revealed an energetically favorable process with ΔG _binding = -44.9 ± 3.3 kcal/mol for the Aβ₄₂ monomer-4v complex. The free energy landscape analysis highlighted that the Aβ₄₂ monomer-4v complex sampled conformations with significantly higher helical contents (35 and 49%) as compared to the Aβ₄₂ monomer alone (17%). Compound 4v displayed hydrogen bonding with Gly37 (chain E) and π-π interactions with Phe19 (chain D) of the Aβ₄₂ protofibril. Further, the per-residue binding free energy analysis also highlighted that Phe19 (chain D) and Gly37 (chain E) of the Aβ₄₂ protofibril showed the maximum contribution in the binding free energy. The decreased binding affinity and residue-residue contacts between chains D and E of the Aβ₄₂ protofibril in the presence of 4v indicate destabilization of the Aβ₄₂ protofibril structure. Overall, the structural information obtained through MD simulations indicated that 4v stabilizes the native helical conformation of the Aβ₄₂ monomer and persuades a destabilization in the protofibril structure of Aβ₄₂. The results of the study will be useful in the rational design of potent inhibitors against amyloid aggregation.

High-Speed Atomic Force Microscopy Reveals the Structural Dynamics of the Amyloid-β and Amylin Aggregation Pathways

Individual Alzheimer's disease (AD) patients have been shown to have structurally distinct amyloid-β (Aβ) aggregates, including fibrils, in their brain. These findings suggest the possibility of a relationship between AD progression and Aβ fibril structures. Thus, the characterization of the structural dynamics of Aβ could aid the development of novel therapeutic strategies and diagnosis. Protein structure and dynamics have typically been studied separately. Most of the commonly used biophysical approaches are limited in providing substantial details regarding the combination of both structure and dynamics. On the other hand, high-speed atomic force microscopy (HS-AFM), which simultaneously visualizes an individual protein structure and its dynamics in liquid in real time, can uniquely link the structure and the kinetic details, and it can also unveil novel insights. Although amyloidogenic proteins generate heterogeneously aggregated species, including transient unstable states during the aggregation process, HS-AFM elucidated the structural dynamics of individual aggregates in real time in liquid without purification and isolation. Here, we review and discuss the HS-AFM imaging of amyloid aggregation and strategies to optimize the experiments showing findings from Aβ and amylin, which is associated with type II diabetes, shares some common biological features with Aβ, and is reported to be involved in AD.

Metal ion coordination delays amyloid-β peptide self-assembly by forming an aggregation-inert complex

A detailed understanding of the molecular pathways for amyloid-β (Aβ) peptide aggregation from monomers into amyloid fibrils, a hallmark of Alzheimer's disease, is crucial for the development of diagnostic and therapeutic strategies. We investigate the molecular details of peptide fibrillization in vitro by perturbing this process through addition of differently charged metal ions. Here, we used a monovalent probe, the silver ion, that, similarly to divalent metal ions, binds to monomeric Aβ peptide and efficiently modulates Aβ fibrillization. On the basis of our findings, combined with our previous results on divalent zinc ions, we propose a model that links the microscopic metal-ion binding to Aβ monomers to its macroscopic impact on the peptide self-assembly observed in bulk experiments. We found that substoichiometric concentrations of the investigated metal ions bind specifically to the N-terminal region of Aβ, forming a dynamic, partially compact complex. The metal-ion bound state appears to be incapable of aggregation, effectively reducing the available monomeric Aβ pool for incorporation into fibrils. This is especially reflected in a decreased fibril-end elongation rate. However, because the bound state is significantly less stable than the amyloid state, Aβ peptides are only transiently redirected from fibril formation, and eventually almost all Aβ monomers are integrated into fibrils. Taken together, these findings unravel the mechanistic consequences of delaying Aβ aggregation via weak metal-ion binding, quantitatively linking the contributions of specific interactions of metal ions with monomeric Aβ to their effects on bulk aggregation.

Oligomerization and Conformational Change Turn Monomeric β-Amyloid and Tau Proteins Toxic: Their Role in Alzheimer's Pathogenesis

The structural polymorphism and the physiological and pathophysiological roles of two important proteins, β-amyloid (Aβ) and tau, that play a key role in Alzheimer's disease (AD) are reviewed. Recent results demonstrate that monomeric Aβ has important physiological functions. Toxic oligomeric Aβ assemblies (AβOs) may play a decisive role in AD pathogenesis. The polymorph fibrillar Aβ (fAβ) form has a very ordered cross-β structure and is assumed to be non-toxic. Tau monomers also have several important physiological actions; however, their oligomerization leads to toxic oligomers (TauOs). Further polymerization results in probably non-toxic fibrillar structures, among others neurofibrillary tangles (NFTs). Their structure was determined by cryo-electron microscopy at atomic level. Both AβOs and TauOs may initiate neurodegenerative processes, and their interactions and crosstalk determine the pathophysiological changes in AD. TauOs (perhaps also AβO) have prionoid character, and they may be responsible for cell-to-cell spreading of the disease. Both extra- and intracellular AβOs and TauOs (and not the previously hypothesized amyloid plaques and NFTs) may represent the novel targets of AD drug research.

Sclerotiorin Stabilizes the Assembly of Nonfibrillar Abeta42 Oligomers with Low Toxicity, Seeding Activity, and Beta-sheet Content

The self-assembly of the 42-residue amyloid-β peptide, Aβ42, into fibrillar aggregates is associated with neuronal dysfunction and toxicity in Alzheimer's disease (AD) patient brains, suggesting that small molecules acting on this process might interfere with pathogenesis. Here, we present experimental evidence that the small molecule sclerotiorin (SCL), a natural product belonging to the group of azaphilones, potently delays both seeded and nonseeded Aβ42 polymerization in cell-free assays. Mechanistic biochemical studies revealed that the inhibitory effect of SCL on fibrillogenesis is caused by its ability to kinetically stabilize small Aβ42 oligomers. These structures exhibit low β-sheet content and do not possess seeding activity, indicating that SCL acts very early in the amyloid formation cascade before the assembly of seeding-competent, β-sheet-rich fibrillar aggregates. Investigations with NMR WaterLOGSY experiments confirmed the association of Aβ42 assemblies with SCL in solution. Furthermore, using ion mobility-mass spectrometry, we observed that SCL directly interacts with a small fraction of Aβ42 monomers in the gas phase. In comparison to typical amyloid fibrils, small SCL-stabilized Aβ42 assemblies are inefficiently taken up into mammalian cells and have low toxicity in cell-based assays. Overall, these mechanistic studies support a pathological role of stable, β-sheet-rich Aβ42 fibrils in AD, while structures with low β-sheet content may be less relevant.

High-yield Production of Amyloid-β Peptide Enabled by a Customized Spider Silk Domain

During storage in the silk gland, the N-terminal domain (NT) of spider silk proteins (spidroins) keeps the aggregation-prone repetitive region in solution at extreme concentrations. We observe that NTs from different spidroins have co-evolved with their respective repeat region, and now use an NT that is distantly related to previously used NTs, for efficient recombinant production of the amyloid-β peptide (Aβ) implicated in Alzheimer’s disease. A designed variant of NT from Nephila clavipes flagelliform spidroin, which in nature allows production and storage of β-hairpin repeat segments, gives exceptionally high yields of different human Aβ variants as a solubility tag. This tool enables efficient production of target peptides also in minimal medium and gives up to 10 times more isotope-labeled monomeric Aβ peptides per liter bacterial culture than previously reported.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.

Mercury and Alzheimer's Disease: Hg(II) Ions Display Specific Binding to the Amyloid-β Peptide and Hinder Its Fibrillization

Brains and blood of Alzheimer's disease (AD) patients have shown elevated mercury concentrations, but potential involvement of mercury exposure in AD pathogenesis has not been studied at the molecular level. The pathological hallmark of AD brains is deposition of amyloid plaques, consisting mainly of amyloid-β (Aβ) peptides aggregated into amyloid fibrils. Aβ peptide fibrillization is known to be modulated by metal ions such as Cu(II) and Zn(II). Here, we study in vitro the interactions between Aβ peptides and Hg(II) ions by multiple biophysical techniques. Fluorescence spectroscopy and atomic force microscopy (AFM) show that Hg(II) ions have a concentration-dependent inhibiting effect on Aβ fibrillization: at a 1:1 Aβ·Hg(II) ratio only non-fibrillar Aβ aggregates are formed. NMR spectroscopy shows that Hg(II) ions interact with the N-terminal region of Aβ(1-40) with a micromolar affinity, likely via a binding mode similar to that for Cu(II) and Zn(II) ions, i.e., mainly via the histidine residues His6, His13, and His14. Thus, together with Cu(II), Fe(II), Mn(II), Pb(IV), and Zn(II) ions, Hg(II) belongs to a family of metal ions that display residue-specific binding interactions with Aβ peptides and modulate their aggregation processes.

Structural Interaction of Apolipoprotein A-I Mimetic Peptide with Amyloid-β Generates Toxic Hetero-oligomers

Apolipoproteins are involved in pathological conditions of Alzheimer's disease (AD), and it has been reported that truncated apolipoprotein fragments and β-amyloid (Aβ) peptides coexist as neurotoxic heteromers within the plaques. Therefore, it is important to investigate these complexes at the molecular level to better understand their properties and roles in the pathology of AD. Here, we present a mechanistic insight into such heteromerization using a structurally homologue apolipoprotein fragment of apoA-I (4F) complexed with Aβ(M1-42) and characterize their toxicity. The 4F peptide slows down the aggregation kinetics of Aβ(M1-42) by constraining its structural plasticity. NMR and CD experiments identified 4F-Aβ(M1-42) heteromers comprised of unstructured Aβ(M1-42) and helical 4F. A uniform two-fold reduction in ¹⁵N/¹H NMR signal intensities of Aβ(M1-42) with no observable chemical shift perturbation indicated the formation of a large complex, which was further confirmed by diffusion NMR experiments. Microsecond-scale atomistic molecular dynamics simulations showed that 4F interaction with Aβ(M1-42) is electrostatically driven and induces unfolding of Aβ(M1-42). Neurotoxicity profiling of Aβ(M1-42) complexed with 4F confirms a significant reduction in cell viability and neurite growth. Thus, the molecular architecture of heteromerization between 4F and Aβ(M1-42) discovered in this study provides evidence toward our understanding of the role of apolipoproteins or their truncated fragments in exacerbating AD pathology.

Extensive tests and evaluation of the CHARMM36IDPSFF force field for intrinsically disordered proteins and folded proteins

Intrinsically disordered proteins (IDPs) have received increasing attention in recent studies due to their structural heterogeneity and critical biological functions. To fully understand the structural properties and determine accurate ensembles of IDPs, molecular dynamics (MD) simulation was widely used to sample diverse conformations and reveal the structural dynamics. However, the classical state-of-the-art force fields perform well for folded proteins while being unsatisfactory for the simulations of disordered proteins reported in many previous studies. Thus, improved force fields were developed to precisely describe both folded proteins and disordered proteins. Preliminary tests show that our newly developed CHARMM36IDPSFF (C36IDPSFF) force field can well reproduce the experimental observables of several disordered proteins, but more tests on different types of proteins are needed to further evaluate the performance of C36IDPSFF. Here, we extensively simulate short peptides, disordered proteins, and fast-folding proteins as well as folded proteins, and compare the simulated results with the experimental observables. The simulation results show that C36IDPSFF could substantially reproduce the experimental observables for most of the tested proteins but some limitations are also found in the radius of gyration of large disordered proteins and the stability of fast-folding proteins. This force field will facilitate large scale studies of protein structural dynamics and functions using MD simulations.

Spontaneous self-assembly of amyloid β (1-40) into dimers

The self-assembly and fibrillation of amyloid β (Aβ) proteins is the neuropathological hallmark of Alzheimer's disease. However, the molecular mechanism of how disordered monomers assemble into aggregates remains largely unknown. In this work, we characterize the assembly of Aβ (1-40) monomers into dimers using long-time molecular dynamics simulations. Upon interaction, the monomers undergo conformational transitions, accompanied by change of the structure, leading to the formation of a stable dimer. The dimers are stabilized by interactions in the N-terminal region (residues 5-12), in the central hydrophobic region (residues 16-23), and in the C-terminal region (residues 30-40); with inter-peptide interactions focused around the N- and C-termini. The dimers do not contain long β-strands that are usually found in fibrils.

Modulation of Amyloid-β42 Conformation by Small Molecules Through Nonspecific Binding

Aggregation of amyloid-β (Aβ) peptides is a crucial step in the progression of Alzheimer's disease (AD). Identifying aggregation inhibitors against AD has been a great challenge. We report an atomistic simulation study of the inhibition mechanism of two small molecules, homotaurine and scyllo-inositol, which are AD drug candidates currently under investigation. We show that both small molecules promote a conformational change of the Aβ42 monomer toward a more collapsed phase through a nonspecific binding mechanism. This finding provides atomistic-level insights into designing potential drug candidates for future AD treatments.

Importance of time-ordered non-uniform sampling of multi-dimensional NMR spectra of Aβ1-42 peptide under aggregating conditions

Although the order of the time steps in which the non-uniform sampling (NUS) schedule is implemented when acquiring multi-dimensional NMR spectra is of limited importance when sample conditions remain unchanged over the course of the experiment, it is shown to have major impact when samples are unstable. In the latter case, time-ordering of the NUS data points by the normalized radial length yields a reduction of sampling artifacts, regardless of the spectral reconstruction algorithm. The disadvantage of time-ordered NUS sampling is that halting the experiment prior to its completion will result in lower spectral resolution, rather than a sparser data matrix. Alternatively, digitally correcting for sample decay prior to reconstruction of randomly ordered NUS data points can mitigate reconstruction artifacts, at the cost of somewhat lower sensitivity. Application of these sampling schemes to the Alzheimer's amyloid beta (Aβ1-42) peptide at an elevated concentration, low temperature, and 3 kbar of pressure, where approximately 75% of the peptide reverts to an NMR-invisible state during the collection of a 3D 15N-separated NOESY spectrum, highlights the improvement in artifact suppression and reveals weak medium-range NOE contacts in several regions, including the C-terminal region of the peptide.

Diverse Structural Conversion and Aggregation Pathways of Alzheimer's Amyloid-β (1-40)

Complex amyloid aggregation of amyloid-β (1-40) (Aβ_1-40) in terms of monomer structures has not been fully understood. Herein, we report the microscopic mechanism and pathways of Aβ_1-40 aggregation with macroscopic viewpoints through tuning its initial structure and solubility. Partial helical structures of Aβ_1-40 induced by low solvent polarity accelerated cytotoxic Aβ_1-40 amyloid fibrillation, while predominantly helical folds did not aggregate. Changes in the solvent polarity caused a rapid formation of β-structure-rich protofibrils or oligomers via aggregation-prone helical structures. Modulation of the pH and salt concentration transformed oligomers to protofibrils, which proceeded to amyloid formation. We reveal diverse molecular mechanisms underlying Aβ_1-40 aggregation with conceptual energy diagrams and propose that aggregation-prone partial helical structures are key to inducing amyloidogenesis. We demonstrate that context-dependent protein aggregation is comprehensively understood using the macroscopic phase diagram, which provides general insights into differentiation of amyloid formation and phase separation from unfolded and folded structures.

Epitope Mapping by NMR of a Novel Anti-Aβ Antibody (STAB-MAb)

Alzheimer´s Disease (AD) is one of the most common neurodegenerative disorders worldwide. Excess of β-amyloid (Aβ), a peptide with a high propensity to misfold and self-aggregate, is believed to be the major contributor to the observed neuronal degeneration and cognitive decline in AD. Here, we characterize the epitope of a novel anti-Aβ monoclonal antibody, the STAB-MAb, which has previously demonstrated picomolar affinities for both monomers (K_D = 80 pM) and fibrils (K_D = 130 pM) of Aβ(1–42) and has shown therapeutic efficacy in preclinical mouse models of AD. Our findings reveal a widespread epitope that embraces several key Aβ residues that have been previously described as important in the Aβ fibrillation process. Of note, STAB-MAb exhibits a stronger affinity for the N-terminus of Aβ and stabilizes an α-helix conformation in the central to N-terminal region of the peptide, in addition to disrupting a characteristic salt-bridge of a hairpin structure present in fibrils. The NMR derived epitope supports the observed results from ThT-monitored fluorescence and electron microscopy experiments, in which STAB-MAb was shown to inhibit the formation of aggregates and promote disruption of pre-formed fibrils. In combination with the published in vitro and in vivo assays, our study highlights STAB-MAb as a rare and versatile antibody with analytical, diagnostic and therapeutic efficacy.

Preparation and characterization of a highly soluble Aβ1-42 peptide variant

Alzheimer's disease (AD) is a progressive neurological disease marked by the accumulation and deposition of misfolded amyloid beta or Abeta (Aβ) peptide. Two species of Aβ peptides are found in amyloid plaques, Aβ_1-40 and Aβ_1-42, with the latter being the more amyloidogenic of the two. Understanding how and why Aβ peptides misfold, oligomerize and form amyloid plaques requires a detailed understanding of their structure and dynamics. The poor solubility and strong aggregation tendencies of Aβ_1-42 has made the isolation and characterization of its different structural isoforms (monomer, dimer, oligomer, amyloid) exceedingly difficult. Furthermore, while synthetic Aβ_1-42 peptides (Aβ₄₂syn) are readily available, the cost of isotopically labeled peptide is substantial, making their characterization by NMR spectroscopy cost prohibitive. Here we describe the design, cloning, high-level production, isotopic labeling and biophysical characterization of a modified (solubility-tagged) Aβ_1-42 variant that exhibits excellent water solubility and shares similar aggregation properties as wildtype Aβ_1-42. Specifically, we attached six lysines (6K) to the C-terminus of native Aβ_1-42 to create a more soluble, monomeric form of Aβ_1-42 called Aβ₄₂C6K. A gene for the Aβ₄₂C6K was designed, synthesized and cloned into Escherichia coli (E. coli) and the peptide was expressed at milligram levels. The Aβ₄₂C6K peptide was characterized using circular dichroism (CD), NMR, electron microscopy and thioflavin T fluorescence. Its ability to form stable monomers, oligomers and fibrils under different conditions was assessed. Our results indicate that Aβ₄₂C6K stays monomeric at high concentrations (unlike Aβ_1-42) and can be induced to oligomerize and form fibrils like Aβ_1-42. Our novel construct could be used to explore the structure and dynamics of Aβ_1-42 as well as the interaction of ligands with Aβ_1-42 via NMR.

Computational investigation of retro-isomer equilibrium structures: Intrinsically disordered, foldable, and cyclic peptides

We use all-atom modeling and advanced-sampling molecular dynamics simulations to investigate quantitatively the effect of peptide bond directionality on the equilibrium structures of four linear (two foldable, two disordered) and two cyclic peptides. We find that the retro forms of cyclic and foldable linear peptides adopt distinctively different conformations compared to their parents. While the retro form of a linear intrinsically disordered peptide with transient secondary structure fails to reproduce a secondary structure content similar to that of its parent, the retro form of a shorter disordered linear peptide shows only minor differences compared to its parent.

Key Peptides and Proteins in Alzheimer's Disease

Alzheimer's Disease (AD) is a form of progressive dementia involving cognitive impairment, loss of learning and memory. Different proteins (such as amyloid precursor protein (APP), β- amyloid (Aβ) and tau protein) play a key role in the initiation and progression of AD. We review the role of the most important proteins and peptides in AD pathogenesis. The structure, biosynthesis and physiological role of APP are shortly summarized. The details of trafficking and processing of APP to Aβ, the cytosolic intracellular Aβ domain (AICD) and small soluble proteins are shown, together with other amyloid-forming proteins such as tau and α-synuclein (α-syn). Hypothetic physiological functions of Aβ are summarized. The mechanism of conformational change, the formation and the role of neurotoxic amyloid oligomeric (oAβ) are shown. The fibril formation process and the co-existence of different steric structures (U-shaped and S-shaped) of Aβ monomers in mature fibrils are demonstrated. We summarize the known pathogenic and non-pathogenic mutations and show the toxic interactions of Aβ species after binding to cellular receptors. Tau phosphorylation, fibrillation, the molecular structure of tau filaments and their toxic effect on microtubules are shown. Development of Aβ and tau imaging in AD brain and CSF as well as blood biomarkers is shortly summarized. The most probable pathomechanisms of AD including the toxic effects of oAβ and tau; the three (biochemical, cellular and clinical) phases of AD are shown. Finally, the last section summarizes the present state of Aβ- and tau-directed therapies and future directions of AD research and drug development.

Effects of in vivo conditions on amyloid aggregation

One of the grand challenges of biophysical chemistry is to understand the principles that govern protein misfolding and aggregation, which is a highly complex process that is sensitive to initial conditions, operates on a huge range of length- and timescales, and has products that range from protein dimers to macroscopic amyloid fibrils. Aberrant aggregation is associated with more than 25 diseases, which include Alzheimer's, Parkinson's, Huntington's, and type II diabetes. Amyloid aggregation has been extensively studied in the test tube, therefore under conditions that are far from physiological relevance. Hence, there is dire need to extend these investigations to in vivo conditions where amyloid formation is affected by a myriad of biochemical interactions. As a hallmark of neurodegenerative diseases, these interactions need to be understood in detail to develop novel therapeutic interventions, as millions of people globally suffer from neurodegenerative disorders and type II diabetes. The aim of this review is to document the progress in the research on amyloid formation from a physicochemical perspective with a special focus on the physiological factors influencing the aggregation of the amyloid-β peptide, the islet amyloid polypeptide, α-synuclein, and the hungingtin protein.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Memory Loss and the Onset of Alzheimer's Disease Could Be Under the Control of Extracellular Heat Shock Proteins

Alzheimer's disease (AD) is a major contemporary and escalating malady in which amyloid-β (Aβ) peptides are the most likely causative agent. Aβ peptides spontaneously tend to aggregate in extracellular fluids following a progression from a monomeric state, through intermediate forms, ending in amyloid fibers and plaques. It is generally accepted now that the neurotoxic agents leading to cellular death, memory loss, and other AD characteristics are the Aβ intermediate aggregated states. However, Aβ peptides are continuously produced, released into the extracellular space, and rapidly cleared from healthy brains. Coincidentally, members of the heat shock proteins (hsp) family are present in the extracellular medium of healthy cells and body fluids, opening the possibility that hsps and Aβ could meet and interact in the extracellular milieu of the brain. In this perspective and reflection article, we place our investigation showing that the presence of Hsp70s mitigate the formation of low molecular weight Aβ peptide oligomers resulting in a reduction of cellular toxicity, in context of the current understanding of the disease. We propose that it may be an inverse relationship between the presence of Hsp70, the stage of Aβ oligomers, neurotoxicity, and the incidence of AD, particularly since the expression and circulating levels of hsp decrease with aging. Combining these observations, we propose that changes in the dynamics of Hsp70s and Aβ concentrations in the circulating brain fluids during aging defines the control of the formation of Aβ toxic aggregates, thus determining the conditions for neuron degeneration and the incidence of AD.

Interference with Amyloid-β Nucleation by Transient Ligand Interaction

Amyloid-β peptide (Aβ) is an intrinsically disordered protein (IDP) associated with Alzheimer's disease. The structural flexibility and aggregation propensity of Aβ pose major challenges for elucidating the interaction between Aβ monomers and ligands. All-D-peptides consisting solely of D-enantiomeric amino acid residues are interesting drug candidates that combine high binding specificity with high metabolic stability. Here we characterized the interaction between the 12-residue all-D-peptide D3 and Aβ42 monomers, and how the interaction influences Aβ42 aggregation. We demonstrate for the first time that D3 binds to Aβ42 monomers with submicromolar affinities. These two highly unstructured molecules are able to form complexes with 1:1 and other stoichiometries. Further, D3 at substoichiometric concentrations effectively slows down the β-sheet formation and Aβ42 fibrillation by modulating the nucleation process. The study provides new insights into the molecular mechanism of how D3 affects Aβ assemblies and contributes to our knowledge on the interaction between two IDPs.

The combined force field-sampling problem in simulations of disordered amyloid-β peptides

Molecular dynamics simulations of intrinsically disordered proteins (IDPs) can provide high resolution structural ensembles if the force field is accurate enough and if the simulation sufficiently samples the conformational space of the IDP with the correct weighting of sub-populations. Here, we investigate the combined force field-sampling problem by testing a standard force field as well as newer fixed charge force fields, the latter specifically motivated for better description of unfolded states and IDPs, and comparing them with a standard temperature replica exchange (TREx) protocol and a non-equilibrium Temperature Cool Walking (TCW) sampling algorithm. The force field and sampling combinations are used to characterize the structural ensembles of the amyloid-beta peptides Aβ42 and Aβ43, which both should be random coils as shown recently by experimental nuclear magnetic resonance (NMR) and 2D Förster resonance energy transfer (FRET) experiments. The results illustrate the key importance of the sampling algorithm: while the standard force field using TREx is in poor agreement with the NMR J-coupling and nuclear Overhauser effect and 2D FRET data, when using the TCW method, the standard and optimized protein-water force field combinations are in very good agreement with the same experimental data since the TCW sampling method produces qualitatively different ensembles than TREx. We also discuss the relative merit of the 2D FRET data when validating structural ensembles using the different force fields and sampling protocols investigated in this work for small IDPs such as the Aβ42 and Aβ43 peptides.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Effect of Familial Mutations on the Interconversion of α-Helix to β-Sheet Structures in an Amyloid-Forming Peptide: Insight from Umbrella Sampling Simulations

Understanding the initial events of aggregation of amyloid β monomers to form β-sheet rich fibrils is useful for the development of therapeutics for Alzheimer's disease. In this context, the changes in energetics involved in the aggregation of helical amyloid β monomers into β-sheet rich dimers have been investigated using umbrella sampling simulations and density functional theory calculations. The results from umbrella sampling simulations for the free energy profile for the interconversion closely agree with the results of density functional theory calculations. The results reveal that helical peptides converted to β-sheet structures through coil-like conformations as intermediates that are mostly stabilized by intramolecular hydrogen bonds. The stabilization of intermediate structures could be a possible way to inhibit fibril formation. Mutations substantially decrease the height of the energy barrier for interconversion from α-helix to β-sheet structure when compared to that of the wild type, something that is attributed to an increase in the number of intramolecular hydrogen bonds between backbone atoms in the coil structures that correspond to a maximum value on the free energy surface. The reduction of the energy barrier leads to an enhancement of the rate of aggregation of amyloid β monomers upon introduction of various familial mutations, which is consistent with previous experimental reports.