Small Oligomers of Aβ42 Protein in the Bulk Solution with AlphaFold2

A Targetable Self-association Surface of the Huntingtin exon1 Helical Tetramer Required for Assembly of Amyloid Pre-nucleation Oligomers

Polyglutamine (polyQ) sequences undergo repeat-length dependent formation of disease-associated, amyloid-like cross-β core structures with kinetics and aggregate morphologies often influenced by the flanking sequences. In Huntington's disease (HD), the httNT segment on the polyQ's N-terminal flank enhances aggregation rates by changing amyloid nucleation from a classical homogeneous mechanism to a two-step process requiring an ɑ-helix-rich oligomeric intermediate. A folded, helix-rich httNT tetrameric structure suggested to be this critical intermediate was recently reported. Here we employ single alanine replacements along the httNT sequence to assess this proposed structure and refine the mechanistic model. We find that Ala replacement of hydrophobic residues within simple httNT peptides greatly suppresses helicity, supporting the tetramer model. These same helix-disruptive replacements in the httNT segment of an exon-1 analog greatly reduce aggregation kinetics, suggesting that an ɑ-helix rich multimer - either the tetramer or a larger multimer - plays an on-pathway role in nucleation. Surprisingly, several other Ala replacements actually enhance helicity and/or amyloid aggregation. The spatial localization of these residues on the tetramer surface suggests a self-association interface responsible for formation of the octomers and higher-order multimers most likely required for polyQ amyloid nucleation. Multimer docking of the tetramer, using the protein-protein docking algorithm ClusPro, predicts this symmetric surface to be a viable tetramer dimerization interface. Intriguingly, octomer formation brings the emerging polyQ chains into closer proximity at this tetramer-tetramer interface. Further supporting the potential importance of tetramer super-assembly, computational docking with a known exon-1 aggregation inhibitor predicts ligand contacts with residues at this interface.

Peptide Self-Assembly into Amyloid Fibrils: Unbiased All-Atom Simulations

Protein self-assembly plays an important role in biological systems, accounting for the formation of mesoscopic structures that can be highly symmetric as in the capsid of viruses or disordered as in molecular condensates or exhibit a one-dimensional fibrillar morphology as in amyloid fibrils. Deposits of the latter in tissues of individuals with degenerative diseases like Alzheimer's and Parkinson's has motivated extensive efforts to understand the sequence of molecular events accounting for their formation. These studies aim to identify on-pathway intermediates that may be the targets for therapeutic intervention. This detailed knowledge of fibril formation remains obscure, in part due to challenges with experimental analyses of these processes. However, important progress is being achieved for short amyloid peptides due to advances in our ability to perform completely unbiased all-atom simulations of the self-assembly process. This perspective discusses recent developments, their implications, and the hurdles that still need to be overcome to further advance the field.

Recent Computational Advances Regarding Amyloid-β and Tau Membrane Interactions in Alzheimer's Disease

The interactions of amyloid proteins with membranes have been subject to many experimental and computational studies, as these interactions contribute in part to neurodegenerative diseases. In this review, we report on recent simulations that have focused on the adsorption and insertion modes of amyloid-β and tau proteins in membranes. The atomistic-resolution characterization of the conformational changes of these amyloid proteins upon lipid cell membrane and free lipid interactions is of interest to rationally design drugs targeting transient oligomers in Alzheimer's disease.

Multistep molecular mechanisms of Aβ16-22 fibril formation revealed by lattice Monte Carlo simulations

As a model of self-assembly from disordered monomers to fibrils, the amyloid-β fragment Aβ16-22 was subject to past numerous experimental and computational studies. Because dynamics information between milliseconds and seconds cannot be assessed by both studies, we lack a full understanding of its oligomerization. Lattice simulations are particularly well suited to capture pathways to fibrils. In this study, we explored the aggregation of 10 Aβ16-22 peptides using 65 lattice Monte Carlo simulations, each simulation consisting of 3 × 109 steps. Based on a total of 24 and 41 simulations that converge and do not converge to the fibril state, respectively, we are able to reveal the diversity of the pathways leading to fibril structure and the conformational traps slowing down the fibril formation.

Molecular Insights into the Self-Assembly of a Full-Length hIAPP Trimer: β-Protofibril Formed by β-Hairpin Lateral or Longitudinal Association

The fibrillar protein deposits of the human islet amyloid polypeptide (hIAPP) in the pancreatic islet of Langerhans are pathological hallmark of type II diabetes. Extensive experimental studies have revealed that the oligomeric formations of the hIAPP are more toxic than the mature fibrils. Exploring the oligomeric conformations in the early aggregation state is valuable for effective therapeutics. In this work, using the all-atom explicit-solvent replica exchange molecular dynamic (REMD) simulations, we investigated the structural features and the assembly mechanisms of the full-length hIAPP trimer in solution. The hIAPP trimer adopted more β-sheets than a-helix conformations, and three types of ordered conformations including open β-barrel, single-layer, and double-layer U-shaped β-sheet structures with five β-strands were captured in our simulations. A representative single-layer β-sheet conformation with a CCS value of 1400 Ų in our simulations matches exactly the experimentally ESI-IMS-MS-derived hIAPP trimer sample. These five β-strand conformations formed via the β-hairpin lateral and longitudinal association, respectively, showing two β-protofibril formation models. To the best of our knowledge, it is the first time to reveal two routes to β-sheet formation in the hIAPP trimers on the atomic level. The contact probabilities between pairs of the β-stranded residue show that the hydrophobic interactions between the residues F₁₅ ∼ V₁₇ and A₂₅ ∼ L₂₇ are responsible for the inter- and intra-peptide β-hairpin formations. All of these results indicate that the β-sheet formation is the first step in the conformational changes toward pathological aggregation and provides evidence of the β-sheet assembly mechanism into hIAPP aggregation.

Molecular Integrative Analysis of the Inhibitory Effects of Dipeptides on Amyloid β Peptide 1-42 Polymerization

The major pathological feature of Alzheimer's disease (AD) is the aggregation of amyloid β peptide (Aβ) in the brain. Inhibition of Aβ₄₂ aggregation may prevent the advancement of AD. This study employed molecular dynamics, molecular docking, electron microscopy, circular dichroism, staining of aggregated Aβ with ThT, cell viability, and flow cytometry for the detection of reactive oxygen species (ROS) and apoptosis. Aβ₄₂ polymerizes into fibrils due to hydrophobic interactions to minimize free energy, adopting a β-strand structure and forming three hydrophobic areas. Eight dipeptides were screened by molecular docking from a structural database of 20 L-α-amino acids, and the docking was validated by molecular dynamics (MD) analysis of binding stability and interaction potential energy. Among the dipeptides, arginine dipeptide (RR) inhibited Aβ₄₂ aggregation the most. The ThT assay and EM revealed that RR reduced Aβ₄₂ aggregation, whereas the circular dichroism spectroscopy analysis showed a 62.8% decrease in β-sheet conformation and a 39.3% increase in random coiling of Aβ₄₂ in the presence of RR. RR also significantly reduced the toxicity of Aβ₄₂ secreted by SH-SY5Y cells, including cell death, ROS production, and apoptosis. The formation of three hydrophobic regions and polymerization of Aβ₄₂ reduced the Gibbs free energy, and RR was the most effective dipeptide at interfering with polymerization.

Metastable alpha-rich and beta-rich conformations of small Aβ42 peptide oligomers

Probing the structures of amyloid-β (Aβ) peptides in the early steps of aggregation is extremely difficult experimentally and computationally. Yet, this knowledge is extremely important as small oligomers are the most toxic species. Experiments and simulations on Aβ42 monomer point to random coil conformations with either transient helical or β-strand content. Our current conformational description of small Aβ42 oligomers is funneled toward amorphous aggregates with some β-sheet content and rare high energy states with well-ordered assemblies of β-sheets. In this study, we emphasize another view based on metastable α-helix bundle oligomers spanning the C-terminal residues, which are predicted by the machine-learning AlphaFold2 method and supported indirectly by low-resolution experimental data on many amyloid polypeptides. This finding has consequences in developing novel chemical tools and to design potential therapies to reduce aggregation and toxicity.

Novel machine learning approaches revolutionize protein knowledge

Breakthrough methods in machine learning (ML), protein structure prediction, and novel ultrafast structural aligners are revolutionizing structural biology. Obtaining accurate models of proteins and annotating their functions on a large scale is no longer limited by time and resources. The most recent method to be top ranked by the Critical Assessment of Structure Prediction (CASP) assessment, AlphaFold 2 (AF2), is capable of building structural models with an accuracy comparable to that of experimental structures. Annotations of 3D models are keeping pace with the deposition of the structures due to advancements in protein language models (pLMs) and structural aligners that help validate these transferred annotations. In this review we describe how recent developments in ML for protein science are making large-scale structural bioinformatics available to the general scientific community.

Self-Assembly of Amyloid-Beta (Aβ) Peptides from Solution to Near In Vivo Conditions

Understanding the atomistic resolution changes during the self-assembly of amyloid peptides or proteins is important to develop compounds or conditions to alter the aggregation pathways and suppress the toxicity and potentially aid in the development of drugs. However, the complexity of protein aggregation and the transient order/disorder of oligomers along the pathways to fibril are very challenging. In this Perspective, we discuss computational studies of amyloid polypeptides carried out under various conditions, including conditions closely mimicking in vivo and point out the challenges in obtaining physiologically relevant results, focusing mainly on the amyloid-beta Aβ peptides.

Exploring the binding kinetics and behaviors of self-aggregated beta-amyloid oligomers to phase-separated lipid rafts with or without ganglioside-clusters

Lipid binding kinetics and energetics of self-aggregated and disordered beta-amyloid oligomers of various sizes, from solution to lipid raft surfaces, were investigated using MD simulations. Our systems include small (monomers to tetramers) and larger (octamers and dodecamers) oligomers. Our lipid rafts contain saturated and unsaturated phosphatidylcholine (PC), cholesterol, and with or without asymmetrically distributed monosialotetrahexosylganglioside (GM1). All rafts exhibited dynamic and structurally diversified domains including liquid-ordered (Lo), liquid-disordered (Ld), and interfacial Lod domains. For rafts without GM1, all oligomers bound to the Lod domain. For GM1-containing rafts, all small oligomers and most larger oligomers bound specifically to the GM1-clusters embedded in the Lo domain. Lipid-protein binding energies followed an order of GM1 >> unsaturated PC > saturated PC > cholesterol for all rafts. In addition, protein-induced membrane structural disruption increased progressively with the size of the oligomer for the annular lipids surrounding the membrane-bound protein in non-GM1-containing rafts. We propose that the tight binding of beta-amyloid oligomers to the GM1-clusters and the structural perturbation of lipids surrounding the membrane-bound proteins at the Lod domain are early molecular events of the beta-amyloid aggregation process on neuronal membrane surfaces that trigger the onset of Alzheimer's.

Big dynorphin is a neuroprotector scaffold against amyloid β-peptide aggregation and cell toxicity

Amyloid β-peptide (Aβ) misfolding into β-sheet structures triggers neurotoxicity inducing Alzheimer’s disease (AD). Molecules able to reduce or to impair Aβ aggregation are highly relevant as possible AD treatments since they should protect against Aβ neurotoxicity. We have studied the effects of the interaction of dynorphins, a family of opioid neuropeptides, with Aβ₄₀ the most abundant species of Aβ. Biophysical measurements indicate that Aβ₄₀ interacts with Big Dynorphin (BigDyn), lowering the amount of hydrophobic aggregates, and slowing down the aggregation kinetics. As expected, we found that BigDyn protects against Aβ₄₀ aggregates when studied in human neuroblastoma cells by cell survival assays. The cross-interaction between BigDyn and Aβ₄₀ provides insight into the mechanism of amyloid pathophysiology and may open up new therapy possibilities.

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

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

Inhibition Mechanisms of (-)-Epigallocatechin-3-gallate and Genistein on Amyloid-beta 42 Peptide of Alzheimer's Disease via Molecular Simulations

The misfolding and self-assembly of amyloid-beta (Aβ) peptides are one of the most important factors contributing to Alzheimer's disease (AD). This study aims to reveal the inhibition mechanisms of (-)-epigallocatechin-3-gallate (EGCG) and genistein on the conformational changes of Aβ42 peptides by using molecular docking and molecular dynamics (MD) simulation. The results indicate that both EGCG and genistein have inhibitory effects on the conformational transition of Aβ42 peptide. EGCG and genistein reduce the ratio of β-sheet secondary structures of Aβ42 peptide while inducing random coil structures. In terms of hydrophobic interactions in the central hydrophobic core of Aβ42 peptide, the binding affinities of EGCG are significantly larger in comparison with that of genistein. Our findings illustrate the inhibition mechanisms of EGCG and genistein on the Aβ42 peptides and prove that EGCG is a very promising inhibitor in impeding the conformational change of Aβ42 peptide.