Probing oligomerization of amyloid beta peptide in silico

Essential collective dynamics analysis reveals nonlocal interactions of alpha-synuclein38-95 monomers with fibrillar seeds

Alpha-synuclein (α-syn) is a cytoplasmic protein commonly found in the nervous system. In solution, α-syn adopts disordered unfolded conformations, although it can also form α-helices upon binding to lipid membranes. Under conditions that are not yet fully understood, α-syn can misfold and aggregate, giving rise to β-sheet rich amyloid fibrils, which then tend to accumulate in degenerating neurons. This leads to Parkinson's disease (PD) and several other conditions collectively termed synucleinopathies. Development of disease-modifying treatments requires detailed understanding of structure and dynamics of α-syn's misfolded aggregates. We have employed 1000 ns long, all-atom molecular dynamics simulations to investigate the interaction of monomeric α-syn_38-95 fragments, which contain the most important amyloidogenic regions, with preformed fibrillar seeds composed of staggered, β-sheet rich α-syn chains of matching length. The simulations indicate that α-syn_38-95 monomers tend to form aggregates with the fibrillar seeds, although we have not observed alignment of the monomeric chains with β-strands of the fibril. To analyze the stability of these aggregates, we have employed the essential collective dynamics method, which allows making accurate assessment of dynamical coupling across individual atoms in macromolecules and supramolecular complexes. The analysis revealed extensive dynamical coupling across initially monomeric α-syn chains and the fibrillar seeds including distal regions thereof that did not contact the monomer directly. We have discussed structural origins of these long-range interactions, their impacts for the stability of α-syn aggregates, and potential implications for the development of anti-PD treatments.

Significance of native PLGA nanoparticles in the treatment of Alzheimer's disease pathology

Alzheimer's disease (AD) is believed to be triggered by increased levels/aggregation of β-amyloid (Aβ) peptides. At present, there is no effective disease-modifying treatment for AD. Here, we evaluated the therapeutic potential of FDA-approved native poly(d,l-lactide-co-glycolide) (PLGA) nanoparticles on Aβ aggregation and in cellular/animal models of AD. Our results showed that native PLGA can not only suppress the spontaneous aggregation but can also trigger disassembly of preformed Aβ aggregates. Spectroscopic studies, molecular dynamics simulations and biochemical analyses revealed that PLGA, by interacting with the hydrophobic domain of Aβ_1-42, prevents a conformational shift towards the β-sheet structure, thus precluding the formation and/or triggering disassembly of Aβ aggregates. PLGA-treated Aβ samples can enhance neuronal viability by reducing phosphorylation of tau protein and its associated signaling mechanisms. Administration of PLGA can interact with Aβ aggregates and attenuate memory deficits as well as Aβ levels/deposits in the 5xFAD mouse model of AD. PLGA can also protect iPSC-derived neurons from AD patients against Aβ toxicity by decreasing tau phosphorylation. These findings provide unambiguous evidence that native PLGA, by targeting different facets of the Aβ axis, can have beneficial effects in mouse neurons/animal models as well as on iPSC-derived AD neurons - thus signifying its unique therapeutic potential in the treatment of AD pathology.

•

PLGA nanoparticles by interacting with hydrophobic domain inhibit Aβ aggregation.

•

PLGA-mediated inhibition of Aβ aggregation can increase viability of mouse neurons.

•

PLGA administration can attenuate cognitive deficits/pathology in 5xFAD AD mouse model.

•

PLGA can protect iPSC-derived neurons from AD patients against Aβ toxicity.

Impact of ultrafine particles and secondary inorganic ions on early onset and progression of amyloid aggregation: Insights from molecular simulations

Alzheimer's disease (AD) is a neurodegenerative disorder, associated with the aggregation of amyloid beta (Aβ) peptides and formation of plaques. The impact of airborne particulate matter (PM) and ultrafine particles (UFPs), on early onset and progression of AD has been recently hypothesized. Considering their small size, carbon black nanoparticles and UFPs can penetrate into human organism and affect Alzheimer's progression. While experiments show that the exposure of PM and UFPs can lead to enhanced concentrations of Aβ peptides, the interactions between the peptides and UFPs remain obscured. Particularly, the impact of UFPs on the initial rate of aggregation of the peptides is ambiguous. Herein, we perform molecular dynamics simulations to investigate the aggregation of Aβ_16-21 peptides, an aggregation-prone segment of Aβ, in the presence of UFPs, mimicked by C₆₀, under different salt solutions suggesting the presence of the inorganic constituents of PM in the blood. In particular, the simulations were performed in the presence of Na⁺, Cl^- and CO₃^-2 ions to characterize typical buffer environments and electrolytes present in human blood. Furthermore, NH₄⁺_, NO₃^- and SO₄^-2 ions, found in PM, were used in the simulations. The results revealed high propensity for the aggregation of Aβ_16-21 peptides. Moreover, the peptides made clusters with C₆₀ molecules, that would be expected to act as a nucleation site for the formation of amyloid plaques. Taken together, the results showed that UFPs affected the peptide aggregation differently, depending on the type of ions present in the simulation environment. In the presence of C₆₀, SO₄^-2 and NO₃^- ions accelerated the aggregation of Aβ_16-21 peptides, however, NH₄⁺ ions decelerated their aggregation. In addition, UFP lowered β-sheets amounts at all environments, except NaCl solution.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

Investigation of Cyc1 protein structure stability after H53I mutation using computational approaches to improve redox potential

Acidithiobacillus ferrooxidans (Af) is an acidophilic bacterium that grows in rigid surroundings and gets its own energy from the oxidation of Fe²⁺ to Fe³⁺. These bacteria are involved in the bioleaching process. Cyc₁ is a periplasmic protein with a crucial role in electron transportation in the respiratory chain. His53 of the Cyc₁ protein, involved in electron transfer to CoxB, was selected for mutation and bioinformatics studies. His53 was substituted by Ile using PyMol software. Molecular dynamics simulations were performed for wild and mutant types of Cyc₁ protein. The conformational changes of mutated protein were studied by analyzing RMSD, RMSF, SASA, Rg, H Bond, and DSSP. The results of the RMSF analysis indicated an increase in the flexibility of the ligand in the mutant. Finally, active site instability leads to an increase in the value of E₀ at the mutation point and improving electron transfer. On the other, His53 in Cyc₁ is interconnected to Glu126 in CoxB through the water molecule (W76) and hydrogen bonding. In the H53I mutation, there was a decrease in the distance between H2O 2030, 2033, and isoleucine 53, and subsequently, the distance to the water molecule 76 between the two proteins was reduced and strengthens the hydrogen bond between Cyc₁ and CoxB, finally improves electron transfer and the bioleaching process.

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.

Structural stability predictions and molecular dynamics simulations of RBD and HR1 mutations associated with SARS-CoV-2 spike glycoprotein

The COVID-19 pandemic is caused by human transmission and infection of Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2). There is no trusted drug against the virus; hence, efforts are on discovering novel inhibitors for the virus. The entry of a SARS-CoV-2 virus particle into a host cell is initiated by its spike glycoprotein and host Angiotensin-Converting Enzyme 2 (ACE2) receptor interaction. Spike glycoprotein domains, namely, the Receptor Binding Domain (RBD) and Heptad Repeat (HR) domains, are essential for this activity. We have studied the impact of mutations such as A348T, N354D, D364Y, G476S, V483A, S494D in the RBD (319-591), and S939F, S940T, T941A, S943P (912-984) in the HR1 domains of spike glycoprotein. Summarily, we utilized the computational screening algorithms to rank the deleterious, damaging and disease-associated spike glycoprotein mutations. Subsequently, to understand the changes in conformation, flexibility and function of the spike glycoprotein mutants, Molecular Dynamics (MD) simulations were performed. The computational predictions and analysis of the MD trajectories suggest that the RBD and HR1 mutations induce significant phenotypic effects on the pre-binding spike glycoprotein structure, which are presumably consequential to its binding to the receptor and provides lead to design inhibitors against the binding.Communicated by Ramaswamy H. Sarma.

A molecular journey to check the conformational dynamics of tau tubulin kinase 2 mutations associated with Alzheimer's disease

Proteins are one of the most vital components of biological functions. Proteins have evolutionarily conserved structures as the shape and folding pattern predominantly determine their function. Considerable research efforts have been made to study the protein folding mechanism. The misfolding of protein intermediates of large groups form polymers with unwanted aggregates that may initiate various diseases. Amongst the diseases caused by misfolding of proteins, Alzheimer's disease (AD) is one of the most prevalent neuro-disorders which has a worldwide impact on human health. The disease is associated with several vital proteins and single amino acid mutations. Tau tubulin kinase 2 (TTBK2) is one of the kinases which is known to phosphorylate tau and tubulin. The literature strongly supports that the mutations-K50E, D163A, R181E, A184E and K143E are associated with multiple important cellular processes of TTBK2. In this study, to understand the molecular basis of the functional effects of the mutations, we have performed structural modeling for TTBK2 and its mutations, using computational prediction algorithms and Molecular Dynamics (MD) simulations. The MD simulations highlighted the impact of the mutations on the Wild Type (WT) by the conformational dynamics, Free Energy Landscape (FEL) and internal molecular motions, indicating the structural de-stabilization which may lead to the disruption of its biological functions. The destabilizing effect of TTBK2 upon mutations provided valuable information about individuals carrying this mutant which could be used as a diagnostic marker in AD.

Screening Malaria-box compounds to identify potential inhibitors against SARS-CoV-2 Mpro, using molecular docking and dynamics simulation studies

Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) Main protease (M^pro) is one of the vital drug targets amongst all the coronaviruses, as the protein is indispensable for virus replication. The study aimed to identify promising lead molecules against M^pro enzyme through virtual screening of Malaria Venture (MMV) Malaria Box (MB) comprising of 400 experimentally proven compounds. The binding affinities were studied using virtual screening based molecular docking, which revealed five molecules having the highest affinity scores compared to the reference molecules. Utilizing the established 3D structure of M^pro the binding affinity conformations of the docked complexes were studied by Molecular Dynamics (MD) simulations. The MD simulation trajectories were analysed to monitor protein deviation, relative fluctuation, atomic gyration, compactness covariance, residue-residue map and free energy landscapes. Based on the present study outcome, we propose three Malaria_box (MB) compounds, namely, MB_241, MB_250 and MB_266 to be the best lead compounds against M^pro activity. The compounds may be evaluated for their inhibitory activities using experimental techniques.

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.

Folding dynamics of Aβ42 monomer at pH 4.0–7.5 with and without physiological salt conditions – does the β1 or β2 region fold first?

The seeding region of Aβ42 monomer is jointly affected by the solution acidity, ionic distribution of the salt, and charged residues.

Tautomerization Effect of Histidines on Oligomer Aggregation of β-Amyloid(1-40/42) during the Early Stage: Tautomerism Hypothesis for Misfolding Protein Aggregation

As the intrinsic origin of the hypothesis for β-amyloid (Aβ) from Alzheimer's disease, histidine behaviors were found to play a crucial role in Aβ aggregation. To investigate the histidine behaviors during the early stage of aggregation, Aβ40/42 pentamers with different histidine isomer states were simulated at the atomic level. Results show that five Aβ40 (δδδ) and Aβ42 (εδδ) monomers can rapidly decrease the aggregation threshold, promote stable pentamer formation, and increase pentamer contents by 51.8% and 56.7%, respectively, as compared with the values of their wild-type (εεε) counterparts. Additionally, pentamers of Aβ40 (δδδ) and Aβ42 (εδδ) have different aggregation pathways and disassembly species, Tr+D and Te+M, during the growth of the pentamer. This work discloses the significance of histidine tautomerization in Aβ aggregation, implying a potential way to control Aβ aggregation and develop the assembly inhibitors.

Elucidating Self-Assembling Peptide Aggregation via Morphoscanner: A New Tool for Protein-Peptide Structural Characterization

Self-assembling and molecular folding are ubiquitous in Nature: they drive the organization of systems ranging from living creatures to DNA molecules. Elucidating the complex dynamics underlying these phenomena is of crucial importance. However, a tool for the analysis of the various phenomena involved in protein/peptide aggregation is still missing. Here, an innovative software is developed and validated for the identification and visualization of b-structuring and b-sheet formation in both simulated systems and crystal structures of proteins and peptides. The novel software suite, dubbed Morphoscanner, is designed to identify and intuitively represent b-structuring and b-sheet formation during molecular dynamics trajectories, paying attention to temporary strand-strand alignment, suboligomer formation and evolution of local order. Self-assembling peptides (SAPs) constitute a promising class of biomaterials and an interesting model to study the spontaneous assembly of molecular systems in vitro. With the help of coarse-grained molecular dynamics the self-assembling of diverse SAPs is simulated into molten aggregates. When applied to these systems, Morphoscanner highlights different b-structuring schemes and kinetics related to SAP sequences. It is demonstrated that Morphoscanner is a novel versatile tool designed to probe the aggregation dynamics of self-assembling systems, adaptable to the analysis of differently coarsened simulations of a variety of biomolecules.

Effect of pH on Aβ42 Monomer and Fibril-like Oligomers-Decoding in Silico of the Roles of pK Values of Charged Residues

Amyloid beta-peptide (Aβ) is the key to developing Alzheimer's disease. Experiments have confirmed that different acidity influences directly not only the structural morphology and population of Aβ oligomers, but also the toxicity. The atomic-level association between the pH, charged residues, and Aβ properties remains obscure. Herein, conformational changes of Aβ₄₂ monomer, fibril-like trimer, and pentamer in the medium pH range of 4.0-7.5 are studied. The results reveal that, as the pH changes from 7.5 to the isoelectric pH, His6, His13, and His14 are protonated in turn, successively approach the center of mass of folded Aβ monomer, trigger ionic interactions and changes of neighboring turns (Asp7-Ser8, His14-Lys16) and even a distant one (Leu34-Met35), as well as concomitant changes of secondary structure, and promote the conformation transition from unfolded to folded. This observation discloses that protonation can convert these charged residues from originally hydrophilic to "hydrophobic-like". For fibril-like oligomers, the pH susceptibility essentially stems from the pK values of charged residues in the context of the Aβ fibril, and in turn one can predict the dynamic behavior of these residues in the processes of dissociation or stabilization of a fibril by comparing the pK values of residues involved in salt bridges in the normal state with those in the current context. This idea is justified by two fibril models and appears to be applicable to other peptides and their fibril systems.

Modulation of Alzheimer's amyloid β peptide oligomerization and toxicity by extracellular Hsp70

Alzheimer's disease (AD) is a progressive neurodegenerative disorder leading to dementia caused by advanced neuronal dysfunction and death. The most significant symptoms of AD are observed at late stages of the disease when interventions are most likely too late to ameliorate the condition. Currently, the predominant theory for AD is the "amyloid hypothesis," which states that abnormally increased levels of amyloid β (Aβ) peptides result in the production of a variety of aggregates that are neurotoxic. The specific mechanisms for Aβ peptide-induced cytotoxicity have not yet been completely elucidated. However, since the majority of Aβ is released into the extracellular milieu, it is reasonable to assume that toxicity begins outside the cells and makes its way inside where it disrupts the basic cellular process resulting in cell death. There is increasing evidence that hsp, particularly Hsp70, are exported into the extracellular milieu by an active export mechanism independent of cell death. Therefore, both Aβ peptides and Hsp70 may coexist in a common environment during pathological conditions. We observed that Hsp70 affected the Aβ assembling process in vitro preventing oligomer formation. Moreover, the presence of Hsp70 reduced the Aβ peptide-induced toxicity of cultured neurons (N2A cells). These results suggest a potential mechanism for the reduction of the detrimental effects of Aβ peptides in AD.

Membrane-Accelerated Amyloid-β Aggregation and Formation of Cross-β Sheets

Amyloid- β aggregates play a causative role in Alzheimer's disease. These aggregates are a product of the physical environment provided by the basic neuronal membrane, composed of a lipid bilayer. The intrinsic properties of the lipid bilayer allow amyloid- β peptides to nucleate and form well-ordered cross- β sheets within the membrane. Here, we correlate the aggregation of the hydrophobic fragment of the amyloid- β protein, A β 25 - 35 , with the hydrophobicity, fluidity, and charge density of a lipid bilayer. We summarize recent biophysical studies of model membranes and relate these to the process of aggregation in physiological systems.