Nonproductive Binding Modes as a Prominent Feature of Aβ40 Fiber Elongation: Insights from Molecular Dynamics Simulation

Delineating the Role of GxxxG Motif in Amyloidogenesis: A New Perspective in Targeting Amyloid-Beta Mediated AD Pathogenesis

The pursuit of a novel structural motif that can shed light on the key functional attributes is a primary focus in the study of protein folding disorders. Decades of research on Alzheimer's disease (AD) have centered on the Amyloid β (Aβ) pathway, highlighting its significance in understanding the disorder. The diversity in the Aβ pathway and the possible silent tracks which are yet to discover, makes it exceedingly intimidating to the interdisciplinary scientific community. Over the course of AD research, Aβ has consistently been at the forefront of scientific inquiry and discussion. In this review, we epitomize the role of a potential structural motif (GxxxG motif) that may provide a new horizon to the Aβ conflict. We emphasize on how comprehensive understanding of this motif from a structure-function perspective may pave the way for designing novel therapeutics intervention in AD and related diseases.

Solvent Relaxation NMR: A Tool for Real-Time Monitoring Water Dynamics in Protein Aggregation Landscape

Solvent dynamics strongly induce the fibrillation of an amyloidogenic system. Probing the solvation mechanism is crucial as it enables us to predict different proteins' functionalities, such as the aggregation propensity, structural flexibility, and toxicity. This work shows that a straightforward NMR method in conjunction with phenomenological models gives a global and qualitative picture of water dynamics at different concentrations and temperatures. Here, we study amyloid system Aβ40 and its fragment AV20 (A21-V40) and G37L (mutation at Gly37 → Leu of AV20), having different aggregation and toxic properties. The independent validation of this method is elucidated using all-atom classical MD simulation. These two state-of-the-art techniques are pivotal in linking the effect of solvent environment in the near hydration-shell to their aggregation nature. The time-dependent modulation in solvent dynamics probed with the NMR solvent relaxation method can be further adopted to gain insight into amyloidogenesis and link with their toxicity profiles.

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

The metastable states of proteins

The intriguing process of protein folding comprises discrete steps that stabilize the protein molecules in different conformations. The metastable state of protein is represented by specific conformational characteristics, which place the protein in a local free energy minimum state of the energy landscape. The native-to-metastable structural transitions are governed by transient or long-lived thermodynamic and kinetic fluctuations of the intrinsic interactions of the protein molecules. Depiction of the structural and functional properties of metastable proteins is not only required to understand the complexity of folding patterns but also to comprehend the mechanisms of anomalous aggregation of different proteins. In this article, we review the properties of metastable proteins in context of their stability and capability of undergoing atypical aggregation in physiological conditions.

The hazardous effects of the environmental toxic gases on amyloid beta-peptide aggregation: A theoretical perspective

Alzheimer's disease (AD) is one of the leading causes of dementia in elderly people. It has been well documented that the exposure to environmental toxins such as CO, CO₂, SO₂ and NO₂ that are present in the air is considered as a hallmark for the progression of Alzheimer's disease. However, their actual mechanism by which environmental toxin triggers the aggregation of Aβ42 peptide at the molecular and atomic levels remain unknown. In this study, molecular dynamics simulation was carried out to study the aggregation mechanism of the Aβ₄₂ peptide due to its interaction of toxic gas (CO, CO₂, SO₂ and NO₂). During the 400 ns simulation, all the Aβ₄₂ interacted toxic gas (CO, CO₂, SO_2, and NO₂) complexes have smaller Root Mean Square Deviation values when compared to the Aβ₄₂ peptide, which shows that the interaction of toxic gases (CO, CO₂, SO_2, and NO₂) would increase the Aβ₄₂ peptide structural stability. The radius of gyration analysis also supports that Aβ₄₂ interacted CO₂ and SO₂ complexes have the minimum value in the range of 0.95 nm and 1.5 nm. It is accounted that the Aβ₄₂ interacted CO₂ and SO₂ complexes have a greater compact structure in comparison to Aβ₄₂ interacted CO and NO₂ complexes. Furthermore, all the Aβ₄₂ interacted toxic gas (CO, CO₂, SO_2, and NO₂) complexes exhibited an enhanced secondary structural probability for coil and turn regions with a reduced α-helix probability, which indicates that the interaction of toxic gases may enhance the toxicity and aggregation of Aβ₄₂.

Distinct Binding Dynamics, Sites and Interactions of Fullerene and Fullerenols with Amyloid-β Peptides Revealed by Molecular Dynamics Simulations

The pathology Alzheimer's disease (AD) is associated with the self-assembly of amyloid-β (Aβ) peptides into β-sheet enriched fibrillar aggregates. A promising treatment strategy is focused on the inhibition of amyloid fibrillization of Aβ peptide. Fullerene C60 is proved to effectively inhibit Aβ fibrillation while the poor water-solubility restricts its use as a biomedicine agent. In this work, we examined the interaction of fullerene C60 and water-soluble fullerenol C60(OH)6/C60(OH)12 (C60 carrying 6/12 hydroxyl groups) with preformed Aβ40/42 protofibrils by multiple molecular dynamics simulations. We found that when binding to the Aβ42 protofibril, C60, C60(OH)6 and C60(OH)12 exhibit distinct binding dynamics, binding sites and peptide interaction. The increased number of hydroxyl groups C60 carries leads to slower binding dynamics and weaker binding strength. Binding free energy analysis demonstrates that the C60/C60(OH)6 molecule primarily binds to the C-terminal residues 31-41, whereas C60(OH)12 favors to bind to N-terminal residues 4-14. The hydrophobic interaction plays a critical role in the interplay between Aβ and all the three nanoparticles, and the π-stacking interaction gets weakened as C60 carries more hydroxyls. In addition, the C60(OH)6 molecule has high affinity to form hydrogen bonds with protein backbones. The binding behaviors of C60/C60(OH)6/C60(OH)12 to the Aβ40 protofibril resemble with those to Aβ42. Our work provides a detailed picture of fullerene/fullerenols binding to Aβ protofibril, and is helpful to understand the underlying inhibitory mechanism.

Probing transient non-native states in amyloid beta fiber elongation by NMR

Using NMR to probe transient binding of Aβ1-40 monomers to fibers, we find partially bound conformations with the highest degree of interaction near F19-K28 and a lesser degree of interaction near the C-terminus (L34-G37). This represents a shift away from the KLVFFA recognition sequence (residues 16-21) currently used for inhibitor design.