Effect of the English familial disease mutation (H6R) on the monomers and dimers of Aβ40 and Aβ42

Prediction of the 3D conformation of a small peptide vaccine targeting Aβ42 oligomers

The original etiology of Alzheimer's disease (AD) is the deposition of amyloid-beta (Aβ) proteins, which starts from the aggregation of the Aβ oligomers. The optimal therapeutic strategy targeting Aβ oligomer aggregation is the development of AD vaccines. Despite the fact that positive progress has been made for experimental attempts at AD vaccines, the physicochemical and even structural properties of these AD vaccines remain unclear. In this study, through immunoinformatic and molecular dynamics (MD) simulations, we first designed and simulated an alternative of vaccine TAPAS and found that the structure of the alternative can reproduce the 3D conformation of TAPAS determined experimentally. Meanwhile, immunoinformatic methods were used to analyze the physicochemical properties of TAPAS, including immunogenicity, antigenicity, thermal stability, and solubility, which confirm well the efficacy and safety of the vaccine, and validate the scheme reliability of immunoinformatic and MD simulations in designing and simulating the TAPAS vaccine. Using the same scheme, we predicted the 3D conformation of the optimized ACI-24 peptide vaccine, an Aβ peptide with the first 15 residues of Aβ42 (Aβ1-15). The vaccine was verified once to be effective against both full-length Aβ1-42 and truncated Aβ4-42 aggregates, but an experimental 3D structure was absent. We have also explored the immune mechanism of the vaccine at the molecular level and found that the optimized ACI-24 and its analogues can block the growth of either full-length Aβ1-42 or truncated Aβ4-42 pentamer by contacting the hydrophobic residues within the N-terminus and β1 region on the contact surface of either pentamer. Additionally, residues (D1, D7, S8, H13, and Q15) were identified as the key residues of the vaccine to contact either of the two Aβ oligomers. This work provides a feasible implementation scheme of immunoinformatic and MD simulations for the development of AD small peptide vaccines, validating the power of the scheme as a parallel tool to the experimental approaches and injecting molecular-level information into the understanding and design of anti-AD vaccines.

Binding mechanism of full-length Aβ40 peptide to a mixed lipid bilayer

The destructive effect of Aβ peptides on membranes is an important source of its cytotoxicity in the pathogenesis of Alzheimer's disease. We have investigated the binding mechanism between the Aβ42 peptide and bilayer in our former work. However, as another abundant form of Aβ peptides in the physiological environment, the binding mechanism between Aβ40 peptide and the lipid bilayer still remains ambiguous. Hence, we performed all-atom simulations on the Aβ40 peptides with the lipid bilayer herein using replica exchange with the solute tempering 2 method. We obtained four major binding models with the hydrophobic C-terminus as the most preferable binding region. Hydrophobic residues and positively charged residues are the principal residues involved in the peptide-bilayer interactions. Aβ40 peptides in our simulation mainly adopt a β-rich conformation in both bound and unbound states. Besides, we determined peptide-water interactions and found that bound peptides prefer forming hydrogen bonds with water molecules than unbound peptides. Our findings herein may provide new insights for the in-depth understanding of the membrane-destructive mechanism of Aβ peptides.

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.

Binding Models of Aβ42 Peptide with Membranes Explored by Molecular Simulations

One of the factors contributing to the toxicity of amyloid-β (Aβ) peptides is the destruction of membrane integrity through Aβ peptide-membrane interactions. The binding of Aβ peptides to membranes has been studied by experiments and theoretical simulations extensively. The exact binding mechanism, however, still remains elusive. In the present study, the molecular basis of the peptide-bilayer binding mechanism of the full-length Aβ42 monomer with POPC/POPS/CHOL bilayers is investigated by all-atom (AA) simulations. Three main binding models in coil, bend, and turn structures are obtained. Model 1 of the three models with the central hydrophobic core (CHC) buried inside the membrane is the dominant binding model. The structural features of the peptide, the peptide-bilayer interacting regions, the intrapeptide interactions, and peptide-water interactions are studied. The binding of the Aβ42 monomer to the POPC/POPS/CHOL bilayer is also explored by coarse-grained (CG) simulations as a complement. Both the AA and CG simulations show that residues in CHC prefer forming interactions with the bilayer, indicating the crucial role of CHC in peptide-bilayer binding. Our results can provide new insights for the investigation of the peptide-bilayer binding mechanism of the Aβ peptide.

Recognition of Aβ oligomer by LilrB2 acceptor: a tetracoordinated zipper mechanism

Leukocyte immunoglobulin-like receptor B2 (LilrB2) is one of discovered cell surface β-amyloid (Aβ) receptors and taken as a promising therapeutic target for the treatment of Alzheimer's disease (AD). Aβ42 oligomer rather than monomer is toxic to neuronal cells and can directly bind to LilrB2, resulting in synaptic loss and cognitive impairment in the development of AD. Therefore, uncovering the mechanism of interaction between Aβ42 oligomer and LilrB2 becomes the first step to obtain a clear drug target and specific binding sites. Herein, a tetracoordinated mechanism for the Aβ oligomer-LilrB2 binding was first put forward by employing Aβ42 dimer mimic-antiparallel copies of Aβ42 core fragment ¹⁶KLVFFA²¹, to bind LilrB2 as models, in which four key residues (F5/F6/L12/F14) in the Aβ42 mimic are bound strongly with LilrB2 residue(s) or accommodated by four hydrophobic cavities in LilrB2 to generate a stable complex. Bi-dentate binding, however, cannot keep the complex Aβ42 mimic-LilrB2 stable. The inhibitor fluspirilene can disturb the binding of four key residues of Aβ42 to LilrB2, justifying the tetracoordinated zipper mechanism on the other hand.

Impact of A2T and D23N mutations on C99 homodimer conformations

The proteolytic cleavage of C99 by γ-secretase is the last step in the production of amyloid-β (Aβ) peptides. Previous studies have shown that membrane lipid composition, cholesterol concentration, and mutation in the transmembrane helix modified the structures and fluctuations of C99. In this study, we performed atomistic molecular dynamics simulations of the homodimer of the 55-residue congener of the C-terminal domain of the amyloid protein precursor, C99(1-55), in a POPC-cholesterol lipid bilayer, and we compared the conformational ensemble of WT sequence to those of the A2T and D23N variants. These mutations are particularly interesting as the protective Alzheimer's disease (AD) A2T mutation is known to decrease Aβ production, whereas the early onset AD D23N mutation does not affect Aβ production. We found noticeable differences in the structural ensembles of the three sequences. In particular, A2T varies from both WT and D23N by having long-range effects on the population of the extracellular justamembrane helix, the interface between the G29xxx-G33xxx-G37 motifs and the fluctuations of the transmembrane helical topologies.

Effect of Cholesterol Molecules on Aβ1-42 Wild-Type and Mutants Trimers

Alzheimer’s disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain, and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with a solute tempering simulation of an Aβ42 trimer with cholesterol and compare it with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex, involving hydrophobic residues L17–F20 and L30–M35 with a non-negligible contribution of loop residues D22–K28 and N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations moderately impact cholesterol-binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into aggregates but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates; this provides a mechanistic link between cholesterol and Alzheimer’s disease.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure

Amyloid beta (Aβ) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aβ peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aβ experiments, as it allows to explore the binding mechanism between metal ions and Aβ oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aβ peptides. Aβ oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aβ binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aβ. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aβ dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.

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.

The Pathogenesis Mechanism, Structure Properties, Potential Drugs and Therapeutic Nanoparticles against the Small Oligomers of Amyloid-β

Alzheimer's Disease (AD) is a devastating neurodegenerative disease that affects millions of people in the world. The abnormal aggregation of amyloid β protein (Aβ) is regarded as the key event in AD onset. Meanwhile, the Aβ oligomers are believed to be the most toxic species of Aβ. Recent studies show that the Aβ dimers, which are the smallest form of Aβ oligomers, also have the neurotoxicity in the absence of other oligomers in physiological conditions. In this review, we focus on the pathogenesis, structure and potential therapeutic molecules against small Aβ oligomers, as well as the nanoparticles (NPs) in the treatment of AD. In this review, we firstly focus on the pathogenic mechanism of Aβ oligomers, especially the Aβ dimers. The toxicity of Aβ dimer or oligomers, which attributes to the interactions with various receptors and the disruption of membrane or intracellular environments, were introduced. Then the structure properties of Aβ dimers and oligomers are summarized. Although some structural information such as the secondary structure content is characterized by experimental technologies, detailed structures are still absent. Following that, the small molecules targeting Aβ dimers or oligomers are collected; nevertheless, all of these ligands have failed to come into the market due to the rising controversy of the Aβ-related "amyloid cascade hypothesis". At last, the recent progress about the nanoparticles as the potential drugs or the drug delivery for the Aβ oligomers are present.

Pulsed Hydrogen-Deuterium Exchange Reveals Altered Structures and Mechanisms in the Aggregation of Familial Alzheimer's Disease Mutants

Mutations of the Amyloid Precursor Protein, from which the amyloid β peptide Aβ42 is cleaved, are associated with familial Alzheimer's disease. The disease-relevant familial mutations include the Arctic (E22G), Iowa (D23N), Italian (E22K), Dutch (E22Q), Japanese (D7N), English (D6R), and Flemish (A21G) variants. A detailed mechanistic understanding of the aggregation behavior of the mutant peptides at the residue level is, however, still lacking. We report here a study of the aggregation kinetics of these mutants in vitro by pulsed hydrogen-deuterium exchange mass spectrometry (HDX-MS) to obtain a temporally and sequence resolved picture of their self-assembly. For all variants, HDX occurs to give a bimodal distribution representing two soluble classes of aggregates, one protected and one solvent-exposed. There is no evidence of other classes of structural intermediates within the detection limits of the HDX approach. The fractional changes in the bimodal exchange profiles for several regions of Aβ42 reveal that the central and C-terminal peptides gain protection upon fibril formation, whereas the N-terminal regions remain largely solvent-accessible. For these mutants, all peptide fragments follow the same kinetics, acquiring solvent protection at the same time, further supporting that there are no significant populations of intermediate species under our experimental conditions. The results demonstrate the potential of pulsed HDX-MS for resolving the region-specific aggregation behavior of Aβ42 isoforms in solution where X-ray crystallography and solid-state NMR (ssNMR) are challenged.

Small-Molecule Targeted Aβ42 Aggregate Degradation: Negatively Charged Small Molecules Are More Promising than the Neutral Ones

Heavy evidence has confirmed that Aβ42 oligomers are the most neurotoxic aggregates and play a critical role in the occurrence and development of Alzheimer's disease by causing functional neuron death, cognitive damage, and dementia. Disordered Aβ42 oligomers are challenging therapeutic targets, and no drug is currently in clinical use that modifies the properties of their monomeric states. Here, a negatively charged molecule (ER), rather than the neutral TS1 one, is identified by a molecular dynamics simulation method to be more capable of binding and sequestering the intrinsically disordered amyloid-β peptide Aβ42 in its soluble pentameric state as well as its monomeric components. Results reveal that the ERs interact with Aβ and inhibit the primary nucleation pathways in its aggregation process in entropic expansion mechanism for both Aβ42 and Aβ40 oligomers but with opposite characteristics of hydrophobic surface area (HSA). The interaction between Aβ42 oligomer and either charged ER or neutral TS1/TS0 characterizes decreased HSA, and the decrease in ER-involved case is highly visible, consistent with the observations from in silico and in vitro studies. By contrast, the presence of these inhibitors causes the HSA of Aβ40 oligomer to change undetectably and there is even a bit of increase in the histidine isomerized Aβ40 oligomer. The HSA distinction between Aβ42 and Aβ40 oligomer is possibly derived from the different effects of M35-inhibitor interaction, which is analogous to the effect of M35 oxidation. In comparison with the neutral TS1/TS0 inhibitors, ER is more prone to bind the residues located in the central (β1) and C-terminal (β2) regions of Aβ42 peptide, two key nucleation regions for Aβ intramolecular folding, intermolecular aggregation, and assembly. Notably, ER can strongly bind the charged residues, such as K16, K28, D23, to greatly disturb the potential stabilizer (e.g., salt-bridge, etc.) in metastable Aβ42 oligomers and protofibrils. These results illustrate the strategy of overcoming Alzheimer's disease from inhibiting its early stage Aβ aggregation with two kinds of small molecules to alter their behavior for therapeutic purposes and strongly recommend paying more attention to the engineering and development of negatively charged inhibitors, the long-term underappreciated ones, targeting the early stage Aβ aggregates.

Impact of Mutations on the Conformational Transition from α-Helix to β-Sheet Structures in Arctic-Type Aβ40: Insights from Molecular Dynamics Simulations

The amyloid-β (Aβ) protein aggregation into toxic oligomers and fibrils has been recognized as a key player in the pathogenesis of Alzheimer's disease. Recent experiments reported that a double alanine mutation (L17A/F19A) in the central hydrophobic core (CHC) region of [G22]Aβ₄₀ (familial Arctic mutation) diminished the self-assembly propensity of [G22]Aβ₄₀. However, the molecular mechanism behind the decreased aggregation tendency of [A17/A19/G22]Aβ₄₀ is not well understood. Herein, we carried out molecular dynamics simulations to elucidate the structure and dynamics of [G22]Aβ₄₀ and [A17/A19/G22]Aβ₄₀. The results for the secondary structure analysis reveal a significantly increased amount of the helical content in the CHC and C-terminal region of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀. The bending free-energy analysis of D23-K28 salt bridge suggests that the double alanine mutation in the CHC region of [G22]Aβ₄₀ has the potential to reduce the fibril formation rate by 0.57 times of [G22]Aβ₄₀. Unlike [G22]Aβ₄₀, [A17/A19/G22]Aβ₄₀ largely sampled helical conformation, as determined by the minimum energy conformations extracted from the free-energy landscape. The present study provided atomic level details into the experimentally observed diminished aggregation tendency of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀.

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.

Novel Disassembly Mechanisms of Sigmoid Aβ42 Protofibrils by Introduced Neutral and Charged Drug Molecules

Alzheimer's disease (AD) is characterized by fibrillar deposits of amyloid-β (Aβ) peptides and neurofibrillary tangles of Tau proteins. Aβ peptides are composed of 37-49 residues, among which the Aβ₄₂ isoform is particularly toxic and aggregation-prone and is enriched in the plaques of AD brains and thus considered central to the development of AD. Therefore, disaggregation and disruption provide potential therapeutic approaches to reduce, inhibit, and even reverse Aβ aggregation. Here we capture the atomic-level details of the interactions between sigmoid Aβ₄₂ fibril 2MXU or 5KK3 and either natural tanshinone compounds TS1 or TS0 or negatively charged ER, proposing two unprecedented disassembly mechanisms. Natural TS1 or TS0 prefers to insert into the cavity together with part at the surface of the 2MXU to open up the mouth and twist the conformation, destroying the ordered growth of subsequent monomers along the fibril axis. For the more compact two-fold 5KK3 , attachment of TS1 or TS0 at the surface including some inserted in cavity results in the separation of the two folds. In the two sigmoid fibril systems, it is no longer applicable for the routine criteria to assess Aβ₄₂ fibril disassembly by introduction of these drugs, such as either reduced H-bond number, decreased β-sheet contents, or both. ER, like-charged to Aβ₄₂ fibril, is especially exceptional, and departs utterly from the neutral ones to disassemble Aβ₄₂ fibril. Besides the inapplicable routine criteria, positive binding energy between ER and Aβ₄₂ fibril also deviates from the hypotheses of "ligands exhibiting greater affinity for the β-amyloid peptide are effective at altering its aggregation and inhibiting cell toxicity" ( Cairo et al. , Biochemistry 2002 , 41 , 8620 - 8629 ) but results in stronger disassembly effect on the two kinds of sigmoid Aβ₄₂ fibrils than neutral TS0 or TS1. The disassembly power of charged ER molecules derives from its stronger deformation ability to the conformation of Aβ₄₂ fibril than the neutral ones, twisting the one-fold 2MXU into tapered-shape and separating two-fold 5KK3 in two parts further, which is in great agreement with experimental observations ( Irwin et al. Biomacromolecules 2013 , 14 ( 1 ), 264 - 274 ). The unusual disassembly mechanisms fill the gaps and offer an alternative direction in engineering new inhibitors to treat AD.

Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.

Intrinsically disordered proteins are proteins that have no well-defined structure in at least one functional form. Mutations in one amino acid may still affect their function significantly, especially in subtle ways with cumulative adverse effects on health. Here we report on molecular dynamics simulations of a protein that is critical for neuronal health throughout adulthood (brain-derived neurotrophic factor). We investigate the effects of a mutation carried by 30% of human population, which has been widely studied for its association with aging-related and stress-related disorders, reduced volume of the hippocampus, and variations in episodic memory. We identify a molecular mechanism in which the mutation may change the global conformations of the protein and its ability to bind to receptors.

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.

How ARVC-Related Mutations Destabilize Desmoplakin: An MD Study

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a familial heart disease linked to mutations in several desmosomal proteins, but the specific effects of these mutations on the molecular level are poorly understood. Among the many documented ARVC-related genetic variants, a striking hotspot of nine mutations has been identified in the plakin domain of desmoplakin. This hotspot can be found at the meeting point of three different subdomains of desmoplakin: two spectrin repeats and a Src homology 3 domain. We set out to understand the effect of these mutations. We determine, using molecular dynamics simulations, how these mutations affect the mechanics of this interface, performing two different classes of simulations. First, we sample the dynamics of the plakin domain, in particular the tendency of the interdomain hinge to buckle, and then we apply an external force onto the constructs and determine the force necessary to break them. We find that surface-exposed mutations are not affecting the dynamics to a very large degree but that most buried mutations make the junction more flexible and decrease the rupture forces observed. Our data suggest that buried ARVC mutations destabilize desmoplakin and thereby impair desmosome integrity under tension.

How Fluorescent Tags Modify Oligomer Size Distributions of the Alzheimer Peptide

Within the complex aggregation process of amyloidogenic peptides into fibrils, early stages of aggregation play a central role and reveal fundamental properties of the underlying mechanism of aggregation. In particular, low-molecular-weight aggregates of the Alzheimer amyloid-β peptide (Aβ) have attracted increasing interest because of their role in cytotoxicity and neuronal apoptosis, typical of aggregation-related diseases. One of the main techniques used to characterize oligomeric stages is fluorescence spectroscopy. To this end, Aβ peptide chains are functionalized with fluorescent tags, often covalently bound to the disordered N-terminus region of the peptide, with the assumption that functionalization and presence of the fluorophore will not modify the process of self-assembly nor the final fibrillar structure. In this investigation, we systematically study the effects of four of the most commonly used fluorophores on the aggregation of Aβ (1–40). Time-resolved and single-molecule fluorescence spectroscopy have been chosen to monitor the oligomer populations at different fibrillation times, and transmission electron microscopy, atomic force microscopy and x-ray diffraction to investigate the structure of mature fibrils. Although the structures of the fibrils were only slightly affected by the fluorescent tags, the sizes of the detected oligomeric species varied significantly depending on the chosen fluorophore. In particular, we relate the presence of high-molecular-weight oligomers of Aβ (1–40) (as found for the fluorophores HiLyte 647 and Atto 655) to net-attractive, hydrophobic fluorophore-peptide interactions, which are weak in the case of HiLyte 488 and Atto 488. The latter leads for Aβ (1–40) to low-molecular-weight oligomers only, which is in contrast to Aβ (1–42). The disease-relevant peptide Aβ (1–42) displays high-molecular-weight oligomers even in the absence of significant attractive fluorophore-peptide interactions. Hence, our findings reveal the potentially high impact of the properties of fluorophores on transient aggregates, which needs to be included in the interpretation of experimental data of oligomers of fluorescently labeled peptides.

Compound CID 9998128 Is a Potential Multitarget Drug for Alzheimer's Disease

We have probed small molecule compound CID 9998128 as a potential multitarget drug for the Alzheimer's disease (AD) using in silico and in vitro experiments. By all-atom simulation and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method, we have demonstrated that this compound strongly binds to both amyloid β42 (Aβ₄₂) fibrils and β-secretase, and the van der Waals interaction dominates over the electrostatic interaction in binding affinity. A detailed analysis at the atomic level revealed that indazole in CID 99998128 structure made a major contribution to instability of all studied complexes. In vitro experiments have shown that CID 9998128 inhibits the Aβ₄₂ amyloid fibrillization and is capable to clear Aβ₄₂ fibrils. Moreover, the compound dose-dependently decreases β-site amyloid precursor protein cleaving enzyme (BACE-1) activity with EC₅₀ value in micromolar range. Thus, our study has revealed that CID 9998128 is a good candidate for AD treatment through preventing production of Aβ peptides and degrading their aggregates. For drug design, we predict that the chemical structure of potent AD multitarget inhibitors should not contain indazole.

Glycation induces conformational changes in the amyloid-β peptide and enhances its aggregation propensity: molecular insights

The cytotoxicity of the amyloid beta (Aβ) peptide, implicated in the pathogenesis of Alzheimer's disease (AD), can be enhanced by its post-translational glycation, a series of non-enzymatic reactions with reducing sugars and reactive dicarbonyls. However, little is known about the underlying mechanisms that potentially enhance the cytotoxicity of the advanced glycation modified Aβ. In this work, fully atomistic molecular dynamics (MD) simulations are exploited to obtain direct molecular insights into the process of early Aβ self-assembly in the presence and absence of glycated lysine residues. Analyses of data exceeding cumulative timescales of 1 microsecond for each system reveal that glycation results in a stronger enthalpy of association between Aβ monomers and lower conformational entropy, in addition to a sharp overall increase in the beta-sheet content. Further analyses reveal that the enhanced interactions originate, in large part, due to markedly stronger, as well as new, inter-monomer salt bridging propensities in the glycated variety. Interestingly, these conformational and energetic effects are broadly reflected in preformed protofibrillar forms of Aβ small oligomers modified with glycation. Our combined results imply that glycation consolidates Aβ self-assembly regardless of its point of occurrence in the pathway. They provide a basis for further mechanistic studies and therapeutic endeavors that could potentially result in novel ways of combating AGE related AD progression.

N-Terminus Binding Preference for Either Tanshinone or Analogue in Both Inhibition of Amyloid Aggregation and Disaggregation of Preformed Amyloid Fibrils-Toward Introducing a Kind of Novel Anti-Alzheimer Compounds

Amyloid-β (Aβ₄₀/Aβ₄₂) peptide with a length of 40 or 42 residues is naturally secreted as cleavage product of the amyloid precursor protein, and formation of Aβ aggregates in a patient's brain is a hallmark of Alzheimer's disease (AD). Therefore, disaggregation and disruption provide potential therapeutic approaches to reduce, inhibit, and even reverse Aβ aggregation. The disaggregation/inhibition effect of the inhibitors applies generally to both Aβ₄₀ and Aβ₄₂ aggregations. Here we capture the atomic-level details of the interaction between Aβ₄₀/Aβ₄₂ and either natural tanshinone compound TS1 or its derivative TS0, and observe novel results by using molecular dynamics simulations. We observe that the natural TS1 indeed inhibits the monomolecular Aβ₄₂ (mAβ₄₂) aggregation and disaggregates Aβ₄₂ amyloid fibrils, being in good agreement with the experimental results. TS1 is favorable to stabilize mAβ₄₀ and even Aβ₄₀ fibril, playing an opposite role to that in the Aβ₄₂ counterpart, however. TS0 can inhibit the misfolding of either mAβ₄₀ or mAβ₄₂ and disaggregate Aβ₄₂ fibril but stabilize the Aβ₄₀ fibril. Using a combination of secondary structural analysis, MM-PBSA binding energy calculations, and radial distribution functions computations, we find that both TS0 and TS1, especially the former, prefer to bind at the charged residues within disordered N-terminus with a scarce positive binding energy and disappear the characteristic C-terminal bend region of Aβ₄₂ fibril, as well as twist the Aβ₄₂ fibril seriously. It turns out to destabilize the Aβ₄₂ fibril and enable the conversion of U-shaped Aβ₄₂ fibril from the onefold to the twofold morphologies. The N-terminal binding preference helps us to identify N-terminal region as the specific epitope for specific inhibitors/drugs (such as TS0 and analogues), heralding unusual inhibition/disaggregation or stabilization mechanisms, and offering an alternative direction in engineering new inhibitors to treat AD.

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

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

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

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

Kinetics of the Interactions between Copper and Amyloid-β with FAD Mutations and Phosphorylation at the N terminus

Mutations and post‐translational modifications of amyloid‐β (Aβ) peptide in its N terminus have been shown to increase fibril formation, yet the molecular mechanism is not clear. Here we investigated the kinetics of the interactions of copper with two Aβ peptides containing Familial Alzheimer's disease (FAD) mutations (English (H6R) and Tottori (D7N)), as well as with Aβ peptide phosphorylated at serine 8 (pS8). All three peptides bind to copper with a similar rate as the wild‐type (wt). The dissociation rates follow the order pS8>H6R>wt>D7N; the interconversion between the two coordinating species occurs 50 % faster for H6R and pS8, whereas D7N had only a negligible effect. Interestingly, the rate of ternary complex (copper‐bridged heterodimer) formation for the modified peptides was significantly faster than that for wt, thus leading us to propose that FAD and sporadic AD might share a kinetic origin for the enhanced oligomerisation of Aβ.

Critical Nucleus Structure and Aggregation Mechanism of the C-terminal Fragment of Copper-Zinc Superoxide Dismutase Protein

The aggregation of the copper-zinc superoxide dismutase (SOD1) protein is linked to familial amyotrophic lateral sclerosis, a progressive neurodegenerative disease. A recent experimental study has shown that the (147)GVIGIAQ(153) SOD1 C-terminal segment not only forms amyloid fibrils in isolation but also accelerates the aggregation of full-length SOD1, while substitution of isoleucine at site 149 by proline blocks its fibril formation. Amyloid formation is a nucleation-polymerization process. In this study, we investigated the oligomerization and the nucleus structure of this heptapeptide. By performing extensive replica-exchange molecular dynamics (REMD) simulations and conventional MD simulations, we found that the GVIGIAQ hexamers can adopt highly ordered bilayer β-sheets and β-barrels. In contrast, substitution of I149 by proline significantly reduces the β-sheet probability and results in the disappearance of bilayer β-sheet structures and the increase of disordered hexamers. We identified mixed parallel-antiparallel bilayer β-sheets in both REMD and conventional MD simulations and provided the conformational transition from the experimentally observed parallel bilayer sheets to the mixed parallel-antiparallel bilayer β-sheets. Our simulations suggest that the critical nucleus consists of six peptide chains and two additional peptide chains strongly stabilize this critical nucleus. The stabilized octamer is able to recruit additional random peptides into the β-sheet. Therefore, our simulations provide insights into the critical nucleus formation and the smallest stable nucleus of the (147)GVIGIAQ(153) peptide.

Allosteric stabilization of the amyloid-β peptide hairpin by the fluctuating N-terminal

Immobilized ions modulate nearby hydrophobic interactions and influence molecular recognition and self-assembly. We simulated disulfide bond-locked double mutants (L17C/L34C) and observed allosteric modulation of the peptide's intra-molecular interactions by the N-terminal tail. We revealed that the non-contacting charged N-terminal residues help the transfer of entropy to the surrounding solvation shell and stabilizing β-hairpin.

Alzheimer's protective A2T mutation changes the conformational landscape of the Aβ₁₋₄₂ monomer differently than does the A2V mutation

The aggregation of amyloid-β (Aβ) peptides plays a crucial role in the etiology of Alzheimer's disease (AD). Recently, it has been reported that an A2T mutation in Aβ can protect against AD. Interestingly, a nonpolar A2V mutation also has been found to offer protection against AD in the heterozygous state, although it causes early-onset AD in homozygous carriers. Since the conformational landscape of the Aβ monomer is known to directly contribute to the early-stage aggregation mechanism, it is important to characterize the effects of the A2T and A2V mutations on Aβ₁₋₄₂ monomer structure. Here, we have performed extensive atomistic replica-exchange molecular dynamics simulations of the solvated wild-type (WT), A2V, and A2T Aβ₁₋₄₂ monomers. Our simulations reveal that although all three variants remain as collapsed coils in solution, there exist significant structural differences among them at shorter timescales. A2V exhibits an enhanced double-hairpin population in comparison to the WT, similar to those reported in toxic WT Aβ₁₋₄₂ oligomers. Such double-hairpin formation is caused by hydrophobic clustering between the N-terminus and the central and C-terminal hydrophobic patches. In contrast, the A2T mutation causes the N-terminus to engage in unusual electrostatic interactions with distant residues, such as K16 and E22, resulting in a unique population comprising only the C-terminal hairpin. These findings imply that a single A2X (where X = V or T) mutation in the primarily disordered N-terminus of the Aβ₁₋₄₂ monomer can dramatically alter the β-hairpin population and switch the equilibrium toward alternative structures. The atomistically detailed, comparative view of the structural landscapes of A2V and A2T variant monomers obtained in this study can enhance our understanding of the mechanistic differences in their early-stage aggregation.

Polymorphic cross-seeding amyloid assemblies of amyloid-β and human islet amyloid polypeptide

Epidemiological studies have shown that the development of Alzheimer's disease (AD) is associated with type 2 diabetes (T2D), but it still remains unclear how AD and T2D are connected. Heterologous cross-seeding between the causative peptides of Aβ and hIAPP may represent a molecular link between AD and T2D. Here, we computationally modeled and simulated a series of cross-seeding double-layer assemblies formed by Aβ and hIAPP peptides using all-atom and coarse-gained molecular dynamics (MD) simulations. The cross-seeding Aβ-hIAPP assemblies showed a wide range of polymorphic structures via a combination of four β-sheet-to-β-sheet interfaces and two packing orientations, focusing on a comparison of different matches of β-sheet layers. Two cross-seeding Aβ-hIAPP assemblies with different interfacial β-sheet packings exhibited high structural stability and favorable interfacial interactions in both oligomeric and fibrillar states. Both Aβ-hIAPP assemblies displayed interfacial dehydration to different extents, which in turn promoted Aβ-hIAPP association depending on interfacial polarity and geometry. Furthermore, computational mutagenesis studies revealed that disruption of interfacial salt bridges largely disfavor the β-sheet-to-β-sheet association, highlighting the importance of salt bridges in the formation of cross-seeding assemblies. This work provides atomic-level information on the cross-seeding interactions between Aβ and hIAPP, which may be involved in the interplay between these two disorders.

Structures of the Alzheimer's Wild-Type Aβ1-40 Dimer from Atomistic Simulations

We have studied the dimer of amyloid beta peptide Aβ of 40 residues by means of all-atom replica exchange molecular dynamics. The Aβ-dimers have been found to be the smallest toxic species in Alzheimer's disease, but their inherent flexibilities have precluded structural characterization by experimental methods. Though the 24-μs-scale simulation reveals a mean secondary structure of 18% β-strand and 10% α helix, we find transient configurations with an unstructured N-terminus and multiple β-hairpins spanning residues 17-21 and 30-36, but the antiparallel and perpendicular peptide orientations are preferred over the parallel organization. Short-lived conformational states also consist of all α topologies, and one compact peptide with β-sheet structure stabilized by a rather extended peptide with α-helical content. Overall, this first all-atom study provides insights into the equilibrium structure of the Aβ1-40 dimer in aqueous solution, opening a new avenue for a comprehensive understanding of the impact of pathogenic and protective mutations in early-stage Alzheimer's disease on a molecular level.

Effect of pH on the Aggregation of α-syn12 Dimer in Explicit Water by Replica-Exchange Molecular Dynamics Simulation

The dimeric structure of the N-terminal 12 residues drives the interaction of α-synuclein protein with membranes. Moreover, experimental studies indicated that the aggregation of α-synuclein is faster at low pH than neutral pH. Nevertheless, the effects of different pH on the structural characteristics of the α-syn12 dimer remain poorly understood. We performed 500 ns temperature replica exchange molecular dynamics (T-REMD) simulations of two α-syn12 peptides in explicit solvent. The free energy surfaces contain ten highly populated regions at physiological pH, while there are only three highly populated regions contained at acidic pH. The anti-parallel β-sheet conformations were found as the lowest free energy state. Additionally, these states are nearly flat with a very small barrier which indicates that these states can easily transit between themselves. The dimer undergoes a disorder to order transition from physiological pH to acidic pH and the α-syn12 dimer at acidic pH involves a faster dimerization process. Further, the Lys6–Asp2 contact may prevent the dimerization.

Coupling of Zinc-Binding and Secondary Structure in Nonfibrillar Aβ40 Peptide Oligomerization

Nonfibrillar neurotoxic amyloid β (Aβ) oligomer structures are typically rich in β-sheets, which could be promoted by metal ions like Zn(2+). Here, using molecular dynamics (MD) simulations, we systematically examined combinations of Aβ40 peptide conformations and Zn(2+) binding modes to probe the effects of secondary structure on Aβ dimerization energies and kinetics. We found that random conformations do not contribute to dimerization either thermodynamically or kinetically. Zn(2+) couples with preformed secondary structures (α-helix and β-hairpin) to speed dimerization and stabilize the resulting dimer. Partial α-helices increase the dimerization speed, and dimers with α-helix rich conformations have the lowest energy. When Zn(2+) coordinates with residues D1, H6, H13, and H14, Aβ40 β-hairpin monomers have the fastest dimerization speed. Dimers with experimentally observed zinc coordination (E11, H6, H13, and H14) form with slower rate but have lower energy. Zn(2+) cannot stabilize fibril-like β-arch dimers. However, Zn(2+)-bound β-arch tetramers have the lowest energy. Collectively, zinc-stabilized β-hairpin oligomers could be important in the nucleation-polymerization of cross-β structures. Our results are consistent with experimental findings that α-helix to β-structural transition should accompany Aβ aggregation in the presence of zinc ions and that Zn(2+) stabilizes nonfibrillar Aβ oligomers and, thus, inhibits formation of less toxic Aβ fibrils.

Curcumin induces structural change and reduces the growth of amyloid-β fibrils: a QCM-D study

Curcumin inhibited Aβ fibril growth through leading to the structural conversion of the growing fibril to a more loosely constructed aggregate.

Crucial role of nonspecific interactions in amyloid nucleation

Protein oligomers have been implicated as toxic agents in a wide range of amyloid-related diseases. However, it has remained unsolved whether the oligomers are a necessary step in the formation of amyloid fibrils or just a dangerous byproduct. Analogously, it has not been resolved if the amyloid nucleation process is a classical one-step nucleation process or a two-step process involving prenucleation clusters. We use coarse-grained computer simulations to study the effect of nonspecific attractions between peptides on the primary nucleation process underlying amyloid fibrillization. We find that, for peptides that do not attract, the classical one-step nucleation mechanism is possible but only at nonphysiologically high peptide concentrations. At low peptide concentrations, which mimic the physiologically relevant regime, attractive interpeptide interactions are essential for fibril formation. Nucleation then inevitably takes place through a two-step mechanism involving prefibrillar oligomers. We show that oligomers not only help peptides meet each other but also, create an environment that facilitates the conversion of monomers into the β-sheet-rich form characteristic of fibrils. Nucleation typically does not proceed through the most prevalent oligomers but through an oligomer size that is only observed in rare fluctuations, which is why such aggregates might be hard to capture experimentally. Finally, we find that the nucleation of amyloid fibrils cannot be described by classical nucleation theory: in the two-step mechanism, the critical nucleus size increases with increases in both concentration and interpeptide interactions, which is in direct contrast with predictions from classical nucleation theory.