Dimerization Mechanism of Alzheimer Aβ40 Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation

The AβA2V paradigm: From molecular insights to therapeutic strategies in Alzheimer's disease and primary tauopathies

Alzheimer's disease, the leading cause of dementia globally, represents an unresolved clinical challenge due to its complex pathogenesis and the absence of effective treatments. Considering the multifactorial etiology of the disease, mainly characterized by the accumulation of amyloid β plaques and neurofibrillary tangles of tau protein, we discuss the A673V mutation in the gene coding for the amyloid precursor protein, which is associated with the familial form of Alzheimer's disease in a homozygous state. The mutation offers new insights into the molecular mechanisms of the disease, particularly regarding the contrasting roles of the A2V and A2T mutations in amyloid β peptide aggregation and toxicity. This review aims to describe relevant studies on A2V-mutated variants of the amyloid β peptide, revealing a protective effect against amyloid-β and tau pathology. Notably, special attention is given to the development of the peptide Aβ1-6A2V(D), which shows significant neuroprotective activity through inhibition of the assembly of amyloid β into amyloid fibrils. The therapeutic potential of this peptide emerges from its ability to reduce amyloid β-induced toxicity, with promising results from studies in human neuroblastoma cells and transgenic animal models.

Structures of Oligomeric States of Tau Protein, Amyloid-β, α-Synuclein and Prion Protein Implicated in Alzheimer's Disease, Parkinson's Disease and Prionopathies

In this review, we focus on the biophysical and structural aspects of the oligomeric states of physiologically intrinsically disordered proteins and peptides tau, amyloid-β and α-synuclein and partly disordered prion protein and their isolations from animal models and human brains. These protein states may be the most toxic agents in the pathogenesis of Alzheimer's and Parkinson's disease. It was shown that oligomers are important players in the aggregation cascade of these proteins. The structural information about these structural states has been provided by methods such as solution and solid-state NMR, cryo-EM, crosslinking mass spectrometry, AFM, TEM, etc., as well as from hybrid structural biology approaches combining experiments with computational modelling and simulations. The reliable structural models of these protein states may provide valuable information for future drug design and therapies.

Molecular Insight into the Effect of HIV-TAT Protein on Amyloid-β Peptides

Increased deposition of amyloid-β (Aβ) plaques in the brain is a frequent pathological feature observed in human immunodeficiency virus (HIV)-positive patients. Emerging evidence indicates that HIV regulatory proteins, particularly the transactivator of transcription (TAT) protein, could interact with Aβ peptide, accelerating the formation of Aβ plaques in the brain and potentially contributing to the onset of Alzheimer's disease in individuals with HIV infection. Nevertheless, the molecular mechanisms underlying these processes remain unclear. In the present study, we have used long all-atom molecular dynamics simulations to probe the direct interactions between the TAT protein and Aβ peptide at the molecular level. Sampling over 28.0 μs, our simulations show that TAT protein induces a shift in the Aβ monomer ensemble toward elongated conformations, exposing aggregation-prone regions on the surface and thereby inducing subsequent aggregation. TAT protein also appears to enhance the stability of preformed Aβ fibrils, while increasing the β-sheet content within these fibrils. Our atomistically detailed simulations qualitatively agree with previous in vitro and in vivo studies. Importantly, our simulations identify key interactions between Aβ and the TAT protein that drive the Aβ aggregation process and stabilize the preformed Aβ aggregates, which are particularly challenging to obtain through current experimental techniques.

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

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

A novel bio-inspired strategy to prevent amyloidogenesis and synaptic damage in Alzheimer's disease

Alzheimer's disease (AD) is an irreversible neurodegenerative disorder that affects millions of people worldwide. AD pathogenesis is intricate. It primarily involves two main molecular players-amyloid-β (Aβ) and tau-which actually have an intrinsic trend to generate molecular assemblies that are toxic to neurons. Incomplete knowledge of the molecular mechanisms inducing the onset and sustaining the progression of the disease, as well as the lack of valid models to fully recapitulate the pathogenesis of human disease, have until now hampered the development of a successful therapy for AD. The overall experience with clinical trials with a number of potential drugs-including the recent outcomes of studies with monoclonal antibodies against Aβ-seems to indicate that Aβ-targeting is not effective if it is not accompanied by an efficient challenge of Aβ neurotoxic properties. We took advantage from the discovery of a naturally-occurring variant of Aβ (AβA2V) that has anti-amyloidogenic properties, and designed a novel bio-inspired strategy for AD based on the intranasal delivery of a six-mer peptide (Aβ1-6A2V) retaining the anti-amyloidogenic abilities of the full-length AβA2V variant. This approach turned out to be effective in preventing the aggregation of wild type Aβ and averting the synaptic damage associated with amyloidogenesis in a mouse model of AD. The results of our preclinical studies inspired by a protective model already existing in nature, that is the human heterozygous AβA2V carriers which seem to be protected from AD, open the way to an unprecedented and promising approach for the prevention of the disease in humans.

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.

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

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

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.

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.

Manganese promotes α-synuclein amyloid aggregation through the induction of protein phase transition

α-Synuclein (α-Syn) is the major protein component of Lewy bodies, a key pathological feature of Parkinson’s disease (PD). The manganese ion Mn²⁺ has been identified as an environmental risk factor of PD. However, it remains unclear how Mn²⁺ regulates α-Syn aggregation. Here, we discovered that Mn²⁺accelerates α-Syn amyloid aggregation through the regulation of protein phase separation. We found that Mn²⁺ not only promotes α-Syn liquid-to-solid phase transition but also directly induces soluble α-Syn monomers to form solid-like condensates. Interestingly, the lipid membrane is integrated into condensates during Mn²⁺-induced α-Syn phase transition; however, the preformed Mn²⁺/α-syn condensates can only recruit lipids to the surface of condensates. In addition, this phase transition can largely facilitate α-Syn amyloid aggregation. Although the Mn²⁺-induced condensates do not fuse, our results demonstrated that they could recruit soluble α-Syn monomers into the existing condensates. Furthermore, we observed that a manganese chelator reverses Mn²⁺-induced α-Syn aggregation during the phase transition stage. However, after maturation, α-Syn aggregation becomes irreversible. These findings demonstrate that Mn²⁺ facilitates α-Syn phase transition to accelerate the formation of α-Syn aggregates and provide new insights for targeting α-Syn phase separation in PD treatment.

Formaldehyde-Crosslinked Nontoxic Aβ Monomers to Form Toxic Aβ Dimers and Aggregates: Pathogenicity and Therapeutic Perspectives

Alzheimer's disease (AD) is characterized by the presence of senile plaques in the brain. However, medicines targeting amyloid-beta (Aβ) have not achieved the expected clinical effects. This review focuses on the formation mechanism of the Aβ dimer (the basic unit of oligomers and fibrils) and its tremendous potential as a drug target. Recently, age-associated formaldehyde and Aβ-derived formaldehyde have been found to crosslink the nontoxic Aβ monomer to form the toxic dimers, oligomers and fibrils. Particularly, Aβ-induced formaldehyde accumulation and formaldehyde-promoted Aβ aggregation form a vicious cycle. Subsequently, formaldehyde initiates Aβ toxicity in both the early-and late-onset AD. These facts also explain why AD drugs targeting only Aβ do not have the desired therapeutic effects. Development of the nanoparticle-based medicines targeting both formaldehyde and Aβ dimer is a promising strategy for improving the drug efficacy by penetrating blood-brain barrier and extracellular space into the cortical neurons in AD patients.

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

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

Dimerization of α-Synuclein Fragments Studied by Isothermal-Isobaric Replica-Permutation Molecular Dynamics Simulation

Aggregates and fibrils of intrinsically disordered α-synuclein are associated with Parkinson's disease. Within a non-amyloid β component (NAC) spanning from the 61st to the 95th residue of α-synuclein, an 11-residue segment called NACore (₆₈GAVVTGVTAVA₇₈) is an essential region for both fibril formation and cytotoxicity. Although NACore peptides alone are known to form aggregates and amyloid fibrils, the mechanisms of aggregation and fibrillation remain unknown. This study investigated the dimerization process of NACore peptides as the initial stage of the aggregation and fibrillation processes. We performed an isothermal-isobaric replica-permutation molecular dynamics simulation, which is one of the efficient sampling methods, for the two NACore peptides in explicit water over 96 μs. The simulation succeeded in sampling a variety of dimer structures. An analysis of secondary structure revealed that most of the NACore dimers form intermolecular β-bridges. In particular, more antiparallel β-bridges were observed than parallel β-bridges. We also found that intramolecular secondary structures such as α-helix and antiparallel β-bridge are stabilized in the pre-dimer state. However, we identified that the intermolecular β-bridges tend to form directly between residues with no specific structure rather than via the intramolecular β-bridges. This is because the NACore peptides still have a low propensity to form the intramolecular secondary structures even though they are stabilized in the pre-dimer state.

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.

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

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

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.

Symmetry-breaking transitions in the early steps of protein self-assembly

Protein misfolding and subsequent self-association are complex, intertwined processes, resulting in development of a heterogeneous population of aggregates closely related to many chronic pathological conditions including Type 2 Diabetes Mellitus and Alzheimer's disease. To address this issue, here, we develop a theoretical model in the general framework of linear stability analysis. According to this model, self-assemblies of peptides with pronounced conformational flexibility may become, under particular conditions, unstable and spontaneously evolve toward an alternating array of partially ordered and disordered monomers. The predictions of the theory were verified by atomistic molecular dynamics (MD) simulations of islet amyloid polypeptide (IAPP) used as a paradigm of aggregation-prone polypeptides (proteins). Simulations of dimeric, tetrameric, and hexameric human-IAPP self-assemblies at physiological electrolyte concentration reveal an alternating distribution of the smallest domains (of the order of the peptide mean length) formed by partially ordered (mainly β-strands) and disordered (turns and coil) arrays. Periodicity disappears upon weakening of the inter-peptide binding, a result in line with the predictions of the theory. To further probe the general validity of our hypothesis, we extended the simulations to other peptides, the Aβ(1-40) amyloid peptide, and the ovine prion peptide as well as to other proteins (SOD1 dimer) that do not belong to the broad class of intrinsically disordered proteins. In all cases, the oligomeric aggregates show an alternate distribution of partially ordered and disordered monomers. We also carried out Surface Enhanced Raman Scattering (SERS) measurements of hIAPP as an experimental validation of both the theory and in silico simulations.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Stability of Aβ11-40 Trimers with Parallel and Antiparallel β-Sheet Organizations in a Membrane-Mimicking Environment by Replica Exchange Molecular Dynamics Simulation

The aggregation of the amyloid (Aβ) peptide of 39-43 amino acids into plaques is observed in the brain of Alzheimer's disease (AD) patients, but the mechanisms underlying the neurotoxicity of Aβ oligomers are still elusive. One suggested initial mechanism is related to the implications of amyloid membrane interactions, but characterization of these assemblies is challenging by experimental means. In this study, we have explored the stability of a trimer of Aβ11-40 in parallel and antiparallel β-sheet structures for the wild-type sequence and its F20W mutant in a dipalmitoylphosphatidylcholine membrane using atomistic replica exchange molecular dynamic simulations. We show that both the U-shape organization and the assembly of β-hairpins are maintained in the membrane and are resistant to the mutation F20W. In contrast the models are destabilized by the F19P mutation. Overall, our results indicate that these two assemblies represent minimal seeds or nuclei for the formation of either amyloid fibrils, a variety of β-barrel pores, or various aggregates for many Aβ sequences in a membrane-mimicking environment.

Aggregation of disease-related peptides

Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.

Role of hydrophobic residues for the gaseous formation of helical motifs

The secondary structure content of proteins and their complexes may change significantly on passing from aqueous solution to the gas phase (as in mass spectrometry experiments). In this work, we investigate the impact of hydrophobic residues on the formation of the secondary structure of a real protein complex in the gas phase. We focus on a well-studied protein complex, the amyloid-β (1-40) dimer (2Aβ). Molecular dynamics simulations reproduce the results of ion mobility-mass spectrometry experiments. In addition, a helix (not present in the solution) is identified involving 19FFAED23, consistent with infrared spectroscopy data on an Aβ segment. Our simulations further point to the role of hydrophobic residues in the formation of helical motifs - hydrophobic sidechains "shield" helices from being approached by residues that carry hydrogen bond sites. In particular, two hydrophobic phenylalanine residues, F19 and F20, play an important role for the helix, which is induced in the gas phase in spite of the presence of two carboxyl-containing residues.

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.

Looking at the Disordered Proteins through the Computational Microscope

Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure-function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.

Computational Study of the Aza-Michael Addition of the Flavonoid (+)-Taxifolin in the Inhibition of β-Amyloid Fibril Aggregation

Inhibition of abnormal protein self-aggregation is an attractive strategy against amyloidogenic diseases, but has found limited success due to the complexity of protein self-assembly, the absence of fully reproducible aggregation assays, and the scarce knowledge of the inhibition mechanisms by small molecules. In this context, catechol-containing compounds may lead to covalent adducts with amyloid fibrils that interfere with the aggregation process. In particular, the covalent adduct formed between the oxidized form of (+)-taxifolin and an β-amyloid (Aβ42) suggests the involvement of a specific recognition motif that enables the chemical reaction with Aβ42. In this study, we have examined the mechanisms implicated in the aza-Michael addition of the o-quinone species of (+)-taxifolin with Aβ42 fibrils. The results support the binding of (+)-taxifolin to the hydrophobic groove delimited by the edges defined by Lys16 and Glu22 residues in the fibril. The chemical reaction proceeds through the nucleophilic attack of the deprotonated amino group of a Lys16 residue in a process activated by the interaction between the o-quinone ring with a vicinal Lys16 residue, as well as by a water-assisted proton transfer, which is the rate-limiting step of the reaction. This specific inhibition mechanism, which may explain the enhanced anti-aggregating activity of oxidized flavonoids compared to fresh compounds, holds promise for developing disease-modifying therapies.

Dihydrochalcone molecules destabilize Alzheimer's amyloid-β protofibrils through binding to the protofibril cavity

Dihydrochalcone molecules destabilize Aβ_17–42protofibrils by disrupting the N-terminal β1 region and the turn region through binding to the protofibril cavity.

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

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

Major Reaction Coordinates Linking Transient Amyloid-β Oligomers to Fibrils Measured at Atomic Level

The structural underpinnings for the higher toxicity of the oligomeric intermediates of amyloidogenic peptides, compared to the mature fibrils, remain unknown at present. The transient nature and heterogeneity of the oligomers make it difficult to follow their structure. Here, using vibrational and solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we show that freely aggregating Aβ₄₀ oligomers in physiological solutions have an intramolecular antiparallel configuration that is distinct from the intermolecular parallel β-sheet structure observed in mature fibrils. The intramolecular hydrogen-bonding network flips nearly 90°, and the two β-strands of each monomeric unit move apart, to give rise to the well-known intermolecular in-register parallel β-sheet structure in the mature fibrils. Solid-state nuclear magnetic resonance distance measurements capture the interstrand separation within monomer units during the transition from the oligomer to the fibril form. We further find that the D23-K28 salt-bridge, a major feature of the Aβ₄₀ fibrils and a focal point of mutations linked to early onset Alzheimer's disease, is not detectable in the small oligomers. Molecular dynamics simulations capture the correlation between changes in the D23-K28 distance and the flipping of the monomer secondary structure between antiparallel and parallel β-sheet architectures. Overall, we propose interstrand separation and salt-bridge formation as key reaction coordinates describing the structural transition of the small Aβ₄₀ oligomers to fibrils.

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.

High-Resolution Structures of the Amyloid-β 1-42 Dimers from the Comparison of Four Atomistic Force Fields

The dimer of the amyloid-β peptide Aβ of 42 residues is the smallest toxic species in Alzheimer's disease, but its equilibrium structures are unknown. Here we determined the equilibrium ensembles generated by the four atomistic OPLS-AA, CHARMM22*, AMBER99sb-ildn, and AMBERsb14 force fields with the TIP3P water model. On the basis of 144 μs replica exchange molecular dynamics simulations (with 750 ns per replica), we find that the four force fields lead to random coil ensembles with calculated cross-collision sections, hydrodynamics properties, and small-angle X-ray scattering profiles independent of the force field. There are, however, marked differences in secondary structure, with the AMBERsb14 and CHARMM22* ensembles overestimating the CD-derived helix content, and the OPLS-AA and AMBER99sb-ildn secondary structure contents in agreement with CD data. Also the intramolecular beta-hairpin content spanning residues 17-21 and 30-36 varies between 1.5% and 13%. Overall, there are significant differences in tertiary and quaternary conformations among all force fields, and the key finding, irrespective of the force field, is that the dimer is stabilized by nonspecific interactions, explaining therefore its possible transient binding to multiple cellular partners and, in part, its toxicity.

Conformational Ensembles of the Wild-Type and S8C Aβ1-42 Dimers

We characterized the dimer of the amyloid-β wild-type (WT) peptide, Aβ, of 42 residues and its disulfide-bond-locked double mutant (S8C) by replica exchange molecular dynamics simulations. Aβ dimers are known to be the smallest toxic species in Alzheimer's disease, and the S8C mutant has been shown experimentally to form an exclusive homogeneous and neurotoxic dimer. Our 50 μs all-atom simulations reveal similar secondary structures and collision cross-sections but very different intramolecular and intermolecular conformations upon double S8C mutation. Both dimers are very dynamic with hundreds of free-energy minima that differ from the U-shape and S-shape conformations of the peptides in the fibrils. The only common structural feature, shared by both species with a probability of 4% in WT and 12% in S8C-S8C, is a three-stranded β-sheet spanning the 17-23, 29-36, and 39-41 residues, which does not exist in the Aβ40 WT dimers.

Stereoisomers Probe Steric Zippers in Amyloid-β

Shape complementarity between close-packed residues plays a critical role in the amyloid aggregation process. Here, we probe such "steric zipper" interactions in amyloid-β (Aβ₄₀), whose aggregation is linked to Alzheimer's disease, by replacing natural residues by their stereoisomers. Such mutations are expected to specifically destabilize the shape sensitive "packing" interactions, which may potentially increase their solubility and change other properties. We study the stereomutants DF19 and DL34 and also the DA2/DF4/DH6/DS8 mutant of Aβ₄₀. F19-L34 is a critical contact in a tightly packed region of Aβ, while residues 1-9 are known to be disordered. While both DF19 and DL34 slow down the kinetics of aggregation and form amyloid fibrils efficiently, only DL34 increases the final solubility. DF19 gives rise to additional off-pathway aggregation which results in large, kinetically stable aggregates, and has lower net solubility. DA2/DF4/DH6/DS8 does not have an effect on the kinetics or the solubility. Notably, both DF19 and DL34 oligomers have a significantly lower level of interactions with lipid vesicles and live cells. We conclude that stereoisomers can cause complex site dependent changes in amyloid properties, and provide an effective tool to determine the role of individual residues in shaping the packed interiors of amyloid aggregates.