Out-of-Register Aβ42 Assemblies as Models for Neurotoxic Oligomers and Fibrils

Structures and Dynamics of β-Rich Oligomers of ATTR (105-115) Assembly

Transthyretin (TTR) is a tetrameric homologous protein that can dissociate into monomers. Misfolding and aggregation of TTR can lead to amyloid transthyretin amyloidosis (ATTR), which can cause many diseases (e.g., senile systemic amyloidosis, familial amyloid cardiomyopathy, and familial amyloid polyneuropathy). Despite growing evidence indicating that small oligomers play a critical role in regulating cytotoxicity, the structures of these oligomeric intermediates and their conformational transformations are still unclear, impeding our understanding of neurodegenerative mechanisms and the development of therapeutics targeting early aggregation species. The TTR monomer protein consists of various fragments prone to self-aggregation, including the residue 105-115 sequence. Therefore, our study investigated the assembly progress of ATTR (105-115) peptides using all-atom molecular dynamics simulations. The findings indicate that the probability of β-sheet content increases with increasing numbers of peptides. Additionally, interactions between hydrophobic residues L110 and L111 are crucial for the formation of a β-rich oligomer formation. These β-rich oligomers may adopt β-barrel conformations, potentially toxic oligomer species. Free-energy analysis reveals that β-barrel conformations serve as intermediates for these β-rich oligomers. Our insights into the structural ensemble dynamics of ATTR (105-115) contribute to understanding the physical mechanisms underlying the β-barrel oligomers of ATTR. These findings may shed light on the pathological role of ATTR in neurodegenerative diseases and offer potential therapeutic targets.

Effect of Terahertz Waves on the Structure of the Aβ42 Monomer, Dimer, and Protofibril: Insights from Molecular Dynamics Simulations

Amyloid-β (Aβ) and its assemblies play important roles in the pathogenesis of Alzheimer's disease (AD). Recent studies conducted by experimental and computational researchers have extensively explored the structure, assembly, and influence of biomolecules and cell membranes on Aβ. However, the impact of terahertz waves on the structures of Aβ monomers and aggregates remains largely unexplored. In this study, we systematically investigate the molecular mechanisms by which terahertz waves affect the structure of the Aβ42 monomer, dimer, and tetramer through all-atom molecular dynamics (MD) simulations. Our findings indicate that terahertz waves at a specific frequency (42.55 THz) can enhance intramolecular and intermolecular interactions in the Aβ42 monomer and dimer, respectively, by resonating with the symmetric stretching mode of the -COO- groups and the symmetric bending/stretching mode of -CH3 groups. Consequently, the β-structure content of the Aβ42 monomer is greatly increased, and the binding energy between the monomers in the Aβ42 dimer is significantly enhanced. Additionally, our observations suggest that terahertz waves can mildly stabilize the structure of tetrameric protofibrils by enhancing the interactions among peripheral peptides. Furthermore, we also investigated the effect of the frequency of terahertz waves on the structure of Aβ42. The present study contributes to a better understanding of the impact of external fields on the biobehavior of Aβ42 peptides and may shed some light on the potential risks associated with electromagnetic field radiation.

Amyloid β Dodecamer Disrupts the Neuronal Membrane More Strongly than the Mature Fibril: Understanding the Role of Oligomers in Neurotoxicity

The amyloid cascade hypothesis states that senile plaques, composed of amyloid β (Aβ) fibrils, play a key role in Alzheimer’s disease (AD). However, recent experiments have shown that Aβ oligomers are more toxic to neurons than highly ordered fibrils. The molecular mechanism underlying this observation remains largely unknown. One of the possible scenarios for neurotoxicity is that Aβ peptides create pores in the lipid membrane that allow Ca²⁺ ions to enter cells, resulting in a signal of cell apoptosis. Hence, one might think that oligomers are more toxic due to their higher ability to create ion channels than fibrils. In this work, we study the effect of Aβ42 dodecamer and fibrils on a neuronal membrane, which is similar to that observed in AD patients, using all-atom molecular dynamics simulations. Due to short simulation times, we cannot observe the formation of pores, but useful insight on the early events of this process has been obtained. Namely, we showed that dodecamer distorts the lipid membrane to a greater extent than fibrils, which may indicate that ion channels can be more easily formed in the presence of oligomers. Based on this result, we anticipate that oligomers are more toxic than mature fibrils, as observed experimentally. Moreover, the Aβ–membrane interaction was found to be governed by the repulsive electrostatic interaction between Aβ and the ganglioside GM1 lipid. We calculated the bending and compressibility modulus of the membrane in the absence of Aβ and obtained good agreement with the experiment. We predict that the dodecamer will increase the compressibility modulus but has little effect on the bending modulus. Due to the weak interaction with the membrane, fibrils insignificantly change the membrane elastic properties.

Mechanistic insights into the mitigation of Aβ aggregation and protofibril destabilization by a D–enantiomeric decapeptide rk10

According to clinical studies, the development of Alzheimer’s disease (AD) is linked to the abnormal aggregation of amyloid-β (Aβ) peptides into toxic soluble oligomers, protofibrils as well as mature fibrils....

Effect of Lauric Acid on the Stability of Aβ42 Oligomers

While Alzheimer's disease is correlated with the presence of Aβ fibrils in patient brains, the more likely agents are their precursors, soluble oligomers that may form pores or otherwise distort cell membranes. Using all-atom molecular dynamics simulation, we study how the presence of fatty acids such as lauric acid changes the stability of pore-forming oligomers built from three-stranded Aβ₄₂ chains. Such a change would alter the distribution of amyloids in the fatty acid-rich brain environment and therefore could explain the lower polymorphism observed in Aβ fibrils derived from brains of patients with Alzheimer's disease. We find that lauric acid stabilizes both ring-like and barrel-shaped models, with the effect being stronger for barrel-like models than for ring-like oligomers.

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.

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

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

Distinguishing the Effect on the Rate and Yield of Aβ42 Aggregation by Green Tea Polyphenol EGCG

Deposition of Aβ42 aggregates in the form of amyloid plaques is a pathological hallmark of Alzheimer's disease. A desired avenue of intervention is the inhibition of Aβ42 aggregation. Epigallocatechin gallate (EGCG), the main polyphenol in green tea, has been generally considered an inhibitor of Aβ aggregation. However, previous experiments focused on the reduction of the amount of Aβ42 aggregates, while the effect of EGCG on the rate of Aβ42 aggregation was not critically analyzed. Here we performed an experimental evaluation of Aβ42 aggregation kinetics in the absence and presence of EGCG at a wide range of concentrations. We found that EGCG reduced thioflavin T fluorescence in an EGCG concentration-dependent manner, suggesting that EGCG reduced the amount of Aβ42 fibrils. The effect of EGCG on the rate of Aβ42 aggregation appears to be bimodal. We found that higher EGCG-to-Aβ42 ratios promoted the rate of Aβ42 aggregation, while lower EGCG-to-Aβ42 ratios inhibited the aggregation rate. To confirm that the reduction of thioflavin T fluorescence is due to the lowered aggregate quantity, but not due to perturbation of thioflavin T binding to Aβ42 fibrils, we probed the effect of EGCG on Aβ42 aggregation using site-directed spin labeling. Electron paramagnetic resonance of spin-labeled Aβ42 aggregates suggests that high EGCG-to-Aβ42 ratios led to a greatly reduced amount of Aβ42 fibrils, and these aggregates adopt similar structures as the fibrils in the no-EGCG sample. Potential implications of this work in designing prevention or therapeutic strategies using EGCG are discussed.

The Role of Structural Polymorphism in Driving the Mechanical Performance of the Alzheimer's Beta Amyloid Fibrils

Alzheimer's Disease (AD) is related with the abnormal aggregation of amyloid β-peptides Aβ₁₋₄₀ and Aβ₁₋₄₂, the latter having a polymorphic character which gives rise to U- or S-shaped fibrils. Elucidating the role played by the nanoscale-material architecture on the amyloid fibril stability is a crucial breakthrough to better understand the pathological nature of amyloid structures and to support the rational design of bio-inspired materials. The computational study here presented highlights the superior mechanical behavior of the S-architecture, characterized by a Young's modulus markedly higher than the U-shaped architecture. The S-architecture showed a higher mechanical resistance to the enforced deformation along the fibril axis, consequence of a better interchain hydrogen bonds' distribution. In conclusion, this study, focusing the attention on the pivotal multiscale relationship between molecular phenomena and material properties, suggests the S-shaped Aβ₁₋₄₂ species as a target of election in computational screen/design/optimization of effective aggregation modulators.

Large fatty acid-derived Aβ42 oligomers form ring-like assemblies

As the primary toxic species in the etiology of Alzheimer disease (AD) are low molecular weight oligomers of Aβ, it is crucial to understand the structure of Aβ oligomers for gaining molecular insights into AD pathology. We have earlier demonstrated that in the presence of fatty acids, Aβ42 peptides assemble as 12-24mer oligomers. These Large Fatty Acid-derived Oligomers (LFAOs) exist predominantly as 12mers at low and as 24mers at high concentrations. The 12mers are more neurotoxic than the 24mers and undergo self-replication, while the latter propagate to morphologically distinct fibrils with succinct pathological consequences. In order to glean into their functional differences and similarities, we have determined their structures in greater detail by combining molecular dynamic simulations with biophysical measurements. We conjecture that the LFAO are made of Aβ units in an S-shaped conformation, with the 12mers forming a double-layered hexamer ring (6 × 2) while the structure of 24mers is a double-layered dodecamer ring (12 × 2). A closer inspection of the (6 × 2) and (12 × 2) structures reveals a concentration and pH dependent molecular reorganization in the assembly of 12 to 24mers, which seems to be the underlying mechanism for the observed biophysical and cellular properties of LFAOs.

Molecular dynamics simulations reveal the mechanism of graphene oxide nanosheet inhibition of Aβ1–42 peptide aggregation

Graphene oxide nanosheets inhibit Aβ1–42 aggregation by weakening inter-peptide interactions and reducing β-sheet contents mostly via salt bridge, hydrogen bonding and cation–π interactions with charged residues.

Exploring Configuration Space and Path Space of Biomolecules Using Enhanced Sampling Techniques-Searching for Mechanism and Kinetics of Biomolecular Functions

To understand functions of biomolecules such as proteins, not only structures but their conformational change and kinetics need to be characterized, but its atomistic details are hard to obtain both experimentally and computationally. Here, we review our recent computational studies using novel enhanced sampling techniques for conformational sampling of biomolecules and calculations of their kinetics. For efficiently characterizing the free energy landscape of a biomolecule, we introduce the multiscale enhanced sampling method, which uses a combined system of atomistic and coarse-grained models. Based on the idea of Hamiltonian replica exchange, we can recover the statistical properties of the atomistic model without any biases. We next introduce the string method as a path search method to calculate the minimum free energy pathways along a multidimensional curve in high dimensional space. Finally we introduce novel methods to calculate kinetics of biomolecules based on the ideas of path sampling: one is the Onsager⁻Machlup action method, and the other is the weighted ensemble method. Some applications of the above methods to biomolecular systems are also discussed and illustrated.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ₄₀ and Aβ₄₂ and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.