Comparing atomistic molecular mechanics force fields for a difficult target: a case study on the Alzheimer’s amyloid β-peptide

Structure Comparison of Beta Amyloid Peptide Aβ 1-42 Isoforms. Molecular Dynamics Modeling

Beta amyloid peptide Aβ 1-42 (Aβ42) has a unique dual role in the human organism, as both the peptide with an important physiological function and one of the most toxic biological compounds provoking Alzheimer's disease (AD). There are several known Aβ42 isoforms that we discuss here that are highly neurotoxic and lead to the early onset of AD. Aβ42 is an intrinsically disordered protein with no experimentally solved structure under physiological conditions. The objective of this research was to establish the appropriate molecular dynamics (MD) methodology and model a uniform set of structures for the Aβ42 isoforms that form the core of this study. For that purpose, force field selection and verification including convergence testing for MD simulations was made. Replica exchange MD and conventional MD modeling of several Aβ42 and Aβ16 isoforms that have neurotoxic and amyloidogenic effects impacting the severity of Alzheimer's disease were carried out with the optimal force field and solvent parameters. A standardized ensemble of structures for the Aβ42 and Aβ16 isoforms covering 30-50% of the conformational ensembles extracted from the free energy minima was calculated from MD trajectories. The resulting data set of modeled structures includes Aβ42 wild type, isoD7, pS8, D7H, and H6R-Aβ42 and Aβ16 wild type, isoD7, pS8, D7H, and H6R-Aβ16. The representative structures are given in the Supporting Information; they are open for public access. In the study, we also evaluated the differences between the structures of Aβ42 isoforms and speculate on their possible relevance to the known functions. Utilizing several representative structures for a single disordered protein for docking, with their subsequent averaging by conformations, would markedly increase the reliability of docking results.

Sterol-activated amyloid beta fibril formation

The metabolic processes that link Alzheimer's disease (AD) to elevated cholesterol levels in the brain are not fully defined. Amyloid beta (Aβ) plaque accumulation is believed to begin decades prior to symptoms and to contribute significantly to the disease. Cholesterol and its metabolites accelerate plaque formation through as-yet-undefined mechanisms. Here, the mechanism of cholesterol (CH) and cholesterol 3-sulfate (CS) induced acceleration of Aβ42 fibril formation is examined in quantitative ligand binding, Aβ42 fibril polymerization, and molecular dynamics studies. Equilibrium and pre-steady-state binding studies reveal that monomeric Aβ42•ligand complexes form and dissociate rapidly relative to oligomerization, that the ligand/peptide stoichiometry is 1-to-1, and that the peptide is likely saturated in vivo. Analysis of Aβ42 polymerization progress curves demonstrates that ligands accelerate polymer synthesis by catalyzing the conversion of peptide monomers into dimers that nucleate the polymerization reaction. Nucleation is accelerated ∼49-fold by CH, and ∼13,000-fold by CS - a minor CH metabolite. Polymerization kinetic models predict that at presumed disease-relevant CS and CH concentrations, approximately half of the polymerization nuclei will contain CS, small oligomers of neurotoxic dimensions (∼12-mers) will contain substantial CS, and fibril-formation lag times will decrease 13-fold relative to unliganded Aβ42. Molecular dynamics models, which quantitatively predict all experimental findings, indicate that the acceleration mechanism is rooted in ligand-induced stabilization of the peptide in non-helical conformations that readily form polymerization nuclei.

Preventing the amyloid-beta peptides accumulation on the cell membrane by applying GHz electric fields: A molecular dynamic simulation

Alzheimer's disease is associated with accumulating different amyloid peptides on the nerve cell membranes. The non-thermal effects of the GHz electric fields in this topic have yet to be well recognized. Hence, in this study, the impacts of 1 and 5 GHz electric fields on the amyloid peptide proteins accumulation on the cell membrane have been investigated, utilizing molecular dynamics (MD) simulation. The obtained results indicated that this range of electric fields did not significantly affect the peptide structure. Moreover, it was found that the peptide penetration into the membrane was increased as the field frequency was increased when the system was exposed to a 20 mv/nm oscillating electric field. In addition, it was observed that the protein-membrane interaction is reduced significantly in the presence of the 70 mv/nm electric field. The molecular level results reported in this study could be helpful in better understanding Alzheimer's disease.

Study on molecular mechanisms of destabilizing Aβ(1-42) protofibrils by licochalcone A and licochalcone B using molecular dynamics simulations

Amyloid-beta (Aβ) protofibrils are closely related to Alzheimer's disease. Their behaviors with or without the presence of Aβ fibrillization inhibitors have been intensively studied by molecular dynamics simulations. In this work, the molecular mechanisms of licochalcone A and licochalcone B on destabilizing Aβ(1-42) protofibrils are explored. It is found that both two licochalcones can disorder the configuration of the Aβ(1-42) protofibril. The stable interactions between the Aβ(1-42) protofibril and licochalcone A or licochalcone B are able to be formed. A reduction of the β-sheet structure contents and an increment of the random coil structures of Aβ(1-42) protofibril are observed in the presence of either licochalcone A or licochalcone B. The hydrogen bonds inside the Aβ(1-42) protofibril could be partially collapsed to varying degrees by two licochalcones. Furthermore, the van der Waals interactions between Aβ(1-42) protofibril and licochalcone A make an important contribution to the binding free energy, while the contribution of the electrostatic interactions between Aβ(1-42) protofibril and licochalcone B is more prominent in the binding affinity. Our work may help in the development of new drug candidates for disrupting the Aβ protofibril.

Effect of carbonaceous ultrafine particles on the structure and oligomerization of Aβ42 peptide

The impact of pervasive air pollutants on human health is a growing concern in scientific communities. Among different air pollutants, ultrafine particles (UFPs; with aerodynamic diameter <100 nm) might pass through biological barriers and have a severe impact on human health, including early progression of neurodegenerative diseases such as Alzheimer's disease (AD). A significant fraction of UFPs consists of carbonaceous compounds, composed of elemental and organic carbon (EC and OC). While in-vivo experimental studies showed the neurotoxicity of typical OC and polycyclic aromatic hydrocarbons (PAHs), the molecular interactions involved in the progression of AD remain unclear. In this study, molecular dynamics simulations were performed to investigate the impact of carbonaceous UFPs on the structure of the Aβ₄₂ monomer and the oligomerization of four Aβ₄₂ peptides, associated with the development of AD. For the simulations, a fullerene (C₆₀) was used for the modeling of EC, while benzo [a]pyrene (B[a]P) was used for the modeling of OC. The results revealed that the presence of C₆₀ accelerated the tetramerization of Aβ₄₂ peptides by 2.5 times, while C₆₀/B[a]P promoted the unfolding of the peptide monomer showing the strongest interactions with the Aβ₄₂ monomer. Similarly, C₆₀/4B[a]P decreased the number of helices in the secondary structure of the peptide monomer by 60%. The simplified UFP models in this study, promoted the early aggregation of peptides to dimers, suggesting the progression of AD.

Modeling structural interconversion in Alzheimers' amyloid beta peptide with classical and intrinsically disordered protein force fields

A comprehensive understanding of the aggregation mechanism in amyloid beta 42 (Aβ42) peptide is imperative for developing therapeutic drugs to prevent or treat Alzheimer's disease. Because of the high flexibility and lack of native tertiary structures of Aβ42, molecular dynamics (MD) simulations may help elucidate the peptide's dynamics with atomic details and collectively improve ensembles not seen in experiments. We applied microsecond-timescale MD simulations to investigate the dynamics and conformational changes of Aβ42 by using a newly developed Amber force field (ff14IDPSFF). We compared the ff14IDPSFF and the regular ff14SB force field by examining the conformational changes of two distinct Aβ42 monomers in explicit solvent. Conformational ensembles obtained by simulations depend on the force field and initial structure, Aβ42^α-helix or Aβ42^β-strand. The ff14IDPSFF sampled a high ratio of disordered structures and diverse β-strand secondary structures; in contrast, ff14SB favored helicity during the Aβ42^α-helix simulations. The conformations obtained from Aβ42^β-strand simulations maintained a balanced content in the disordered and helical structures when simulated by ff14SB, but the conformers clearly favored disordered and β-sheet structures simulated by ff14IDPSFF. The results obtained with ff14IDPSFF qualitatively reproduced the NMR chemical shifts well. In-depth peptide and cluster analysis revealed some characteristic features that may be linked to early onset of the fibril-like structure. The C-terminal region (mainly M35-V40) featured in-registered anti-parallel β-strand (β-hairpin) conformations with tested systems. Our work should expand the knowledge of force field and structure dependency in MD simulations and reveals the underlying structural mechanism-function relationship in Aβ42 peptides. Communicated by Ramaswamy H. Sarma.

Ambient Benzo[a]pyrene's Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer's Disease

Long-time exposure to ambient ultrafine particles is associated with an increased risk of neurodegenerative diseases such as Alzheimer's disease (AD), which is triggered by the aggregation of Aβ peptide monomers into toxic oligomers. Among different ultrafine air pollutants, polycyclic aromatic hydrocarbons (PAHs) are known to have a negative neural impact; however, the impact mechanism remains obscure. We herein examined the effect of Benzo[a]Pyrene (B[a]P), one of the typical PAHs on Aβ₄₂ oligomerization using all-atom molecular dynamics simulations. In particular, the simulations were performed using four molecules of Aβ₄₂ in the presence of 5.00 mM, 12.5 mM, and 50.0 mM of B[a]P. The results revealed strong hydrophobic interactions between Aβ₄₂ peptides and B[a]P, which in turn resulted in increased interpeptide electrostatic interactions. Furthermore, 5.00 mM of B[a]P accelerated the kinetics of the formation of peptide tetramer by 30%, and stabilized C-terminus in Aβ₄₂ peptides, suggesting consequent progression of AD in the presence of 5.00 mM B[a]P. In contrast, 12.5 mM and 50.0 mM of B[a]P decreased interpeptide interactions and H-bonding due to the aggregation of numerous B[a]P clusters with the peptides, suppressing oligomerization kinetics of Aβ₄₂ peptides by 13% and 167%, respectively. While the study elucidates the effect of small environmental hydrophobic molecules on the formation of Aβ oligomers, the impact of ambient ultrafine particles on AD in the complex composition of the environmental realm requires further systematic delving into the field.

Insights into Molecular Mechanisms of EGCG and Apigenin on Disrupting Amyloid-Beta Protofibrils Based on Molecular Dynamics Simulations

The fibrillization and deposition of amyloid-beta (Aβ) protofibrils are one of the important factors leading to Alzheimer's disease. Molecular dynamics simulations can offer information on intermolecular interaction mechanisms between Aβ protofibrils and Aβ fibrillization inhibitors. Here, in this work, we explore the early molecular mechanisms of (-)-epigallocatechin-3-gallate (EGCG) and apigenin on disrupting Aβ42 protofibrils based on molecular simulations. The binding modes of EGCG and apigenin with the Aβ42 protofibril are obtained. Furthermore, we compare the behavioral mechanisms of EGCG and apigenin on disturbing the Aβ42 protofibril. Both EGCG and apigenin are able to decrease the proportion of the β-sheet and bend structures of the Aβ42 protofibril while inducing random coil structures. The results of hydrogen bonds and D23-K28 salt bridges illustrate that EGCG and apigenin have the ability of destabilizing the Aβ42 protofibril. Meanwhile, the van der Waals interactions between the EGCG and Aβ42 protofibril are shown to be larger than that of apigenin with the Aβ42 protofibril, but the electrostatic interactions between apigenin and the Aβ42 protofibril are dominant in the binding affinity. Our findings may help in designing effective drug candidates for disordering the Aβ protofibril and impeding Aβ fibrillization.

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.

Molecular Dynamics Simulations Indicate Aromaticity as a Key Factor in the Inhibition of IAPP(20-29) Aggregation

Islet amyloid polypeptide (IAPP) is a 37-residue amyloidogenic hormone implicated in the progression of Type II Diabetes (T2D). T2D affects an estimated 422 million people yearly and is a comorbidity with numerous diseases. IAPP forms toxic oligomers and amyloid fibrils that reduce pancreatic β-cell mass and exacerbate the T2D disease state. Toxic oligomer formation is attributed, in part, to the formation of interpeptide β-strands comprised of residues 20-29 (IAPP_(20-29)). Flavonoids, a class of polyphenolic natural products, have been found experimentally to inhibit IAPP aggregate formation. Many of these small flavonoids differ structurally only slightly; the influence of functional group placement on inhibiting the aggregation of the IAPP_(20-29) has yet to be explored. To probe the role of small-molecule structural features that impede IAPP aggregation, molecular dynamics simulations were performed to observe trimer formation on a model fragment of IAPP_(20-29) in the presence of morin, quercetin, dihydroquercetin, epicatechin, and myricetin. Contacts between Phe23 residues were critical to oligomer formation, and small-molecule contacts with Phe23 were a key predictor of β-strand reduction. Structural properties influencing the ability of compounds to disrupt Phe23-Phe23 contacts included aromaticity and carbonyl and hydroxyl group placement. This work provides key information on design considerations for T2D therapeutics that target IAPP aggregation.

Simulations of Cross-Amyloid Aggregation of Amyloid-β and Islet Amyloid Polypeptide Fragments

Amyloid-beta (Aβ) and islet amyloid polypeptide (IAPP) are small peptides, classified as amyloids, that have the potential to self-assemble and form cytotoxic species, such as small soluble oligomers and large insoluble fibrils. The formation of Aβ aggregates facilitates the progression of Alzheimer's disease (AD), while IAPP aggregates induce pancreatic β-cell apoptosis, leading to exacerbation of Type 2 diabetes (T2D). Cross-amyloid interactions between Aβ and IAPP have been described both in vivo and in vitro, implying the role of Aβ or IAPP as modulators of cytotoxic self-aggregation of each species, and suggesting that Aβ-IAPP interactions are a potential molecular link between AD and T2D. Using molecular dynamics simulations, "hot spot" regions of the two peptides were studied to understand the formation of hexamers in a heterogenous and homogenous peptide-containing environment. Systems of only Aβ_(16-22) peptides formed antiparallel, β-barrel-like structures, while systems of only IAPP_(20-29) peptides formed stacked, parallel beta sheets and had relatively unstable aggregation structures after 2 μs of simulation time. Systems containing both Aβ and IAPP (1:1 ratio) hexamers showed antiparallel, β-barrel-like structures, with an interdigitated arrangement of Aβ_(16-22) and IAPP_(20-29). These β-barrel structures have features of cytotoxic amyloid species identified in previous literature. Ultimately, this work seeks to provide atomistic insight into both the mechanism behind cross-amyloid interactions and structural morphologies of these toxic amyloid species.

Impact of Electronic Polarization on Preformed, β-Strand Rich Homogenous and Heterogeneous Amyloid Oligomers

Amyloids are a subset of intrinsically disordered proteins (IDPs) that self-assemble into cross-[Formula: see text] oligomers and fibrils. The structural plasticity of amyloids leads to sampling of metastable, low-molecular-weight oligomers that contribute to cytotoxicity. Of interest are amyloid-[Formula: see text] (A[Formula: see text] and islet amyloid polypeptide (IAPP), which are involved in the pathology of Alzheimer’s disease and Type 2 diabetes mellitus, respectively. In addition to forming homogenous oligomers and fibrils, these species have been found to cross-aggregate in heterogeneous structures. Biophysical properties, including electronic effects, that are unique or conserved between homogenous and heterogeneous amyloids oligomers are thus far unexplored. Here, we simulated homogenous and heterogeneous amyloid oligomers of A[Formula: see text] and IAPP[Formula: see text] fragments using the Drude oscillator model to investigate the impact of electronic polarization on the structural morphology and stability of preformed hexamers. Upon simulation of preformed, [Formula: see text]-strand rich oligomers with Drude, structural rearrangement occurred causing some loss of [Formula: see text]-strand structure in favor of random coil content for all oligomers. Homogenous A[Formula: see text] was the most stable system, deriving stability from low polarization in hydrophobic residues and through salt bridge formation. Changes in polarization were observed primarily for A[Formula: see text] residues in heterogeneous cross-amyloid systems, displaying a decrease in charged residue dipole moments and an increase in hydrophobic sidechain dipole moments. This work is the first study utilizing the Drude-2019 force field with amyloid oligomers, providing insight into the impact of electronic effects on oligomer structure and highlighting the importance of different microenvironments on amyloid oligomer stability.

Molecular dynamics study of conformation transition from helix to sheet of Aβ42 peptide

Aβ42 peptides can form helix and sheet structure under different conditions. The conformational conversion is closely associated with Aβ peptides aggregation and their neurotoxicity. But the transition from helix to sheet is not be clearly understood. In this study we performed microsecond timescale MD simulations of Aβ42 peptide to investigate the conformation transition from α-helix to β-sheet. Markov state model (MSM) was built to facilitate identification of crucial intermediate states and possible transition pathway. Based on the analysis, we found that the region Y10-A21 in the middle of Aβ42 peptide plays an initial role in this transition. MSM model revealed that the collapse of helical structure in this region might trigger the formation of sheet structure. Moreover, we further simulated the aggregation of Aβ42 peptides with different conformations. We found that the Aβ42 peptides forming sheet structure have higher aggregation potential compared with peptides with helix structure. These results demonstrate that we can prevent the aggregation of Aβ42 peptides by stabilizing the helix structure in the region of Y10-A21. In addition, this study provides new insight into better understanding the conformational transition and aggregation of Aβ42 peptides.

Adsorption of the amyloid β40 monomer on charged gold nanoparticles and slabs: a molecular dynamics study

Negatively charged nanoparticles are known to inhibit the fibrillation of amyloidogenic protein amyloid β (Aβ40), though the overall charge on the protein is negative. In this work a molecular dynamics study is reported to investigate the interaction of Aβ40 on negatively charged gold nanoparticles (3-5 nm) and charged (positive and negative) and neutral gold slabs. The equilibrium structures of Aβ40 on gold surfaces are characterized using residue-specific contacts on the gold surface, secondary structure analysis and binding free energy calculations. The simulation results reveal that the Aβ40 protein in water interconverts into β-sheets, which are building blocks of the mature fibrils, whereas on gold nanoparticles Aβ40 unfolds and adsorbs. Both the negatively charged gold nanoparticles and gold slabs arrest the formation of β-sheets in Aβ40, whereas the positively charged gold slab does not inhibit the formation of β-sheets. The residue-specific interactions between Aβ40 and the gold surfaces are important in governing the adsorption of Aβ40 on charged surfaces.

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.

Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment

The amlyoid-β peptide (Aβ) is closely linked to the development of Alzheimer's disease. Molecular dynamics (MD) simulations have become an indispensable tool for studying the behavior of this peptide at the atomistic level. General key aspects of MD simulations are the force field used for modeling the peptide and its environment, which is important for accurate modeling of the system of interest, and the length of the simulations, which determines whether or not equilibrium is reached. In this study we address these points by analyzing 30-μs MD simulations acquired for Aβ40 using seven different force fields. We assess the convergence of these simulations based on the convergence of various structural properties and of NMR and fluorescence spectroscopic observables. Moreover, we calculate Markov state models for the different MD simulations, which provide an unprecedented view of the thermodynamics and kinetics of the amyloid-β peptide. This further allows us to provide answers for pertinent questions, like: which force fields are suitable for modeling Aβ? (a99SB-UCB and a99SB-ILDN/TIP4P-D); what does Aβ peptide really look like? (mostly extended and disordered) and; how long does it take MD simulations of Aβ to attain equilibrium? (at least 20-30 μs). We believe the analyses presented in this study will provide a useful reference guide for important questions relating to the structure and dynamics of Aβ in particular, and by extension other similar disordered proteins.

Folding and self-assembly of short intrinsically disordered peptides and protein regions

Proteins and peptide fragments are highly relevant building blocks in self-assembly for nanostructures with plenty of applications. Intrinsically disordered proteins (IDPs) and protein regions (IDRs) are defined by the absence of a well-defined secondary structure, yet IDPs/IDRs show a significant biological activity. Experimental techniques and computational modelling procedures for the characterization of IDPs/IDRs are discussed. Directed self-assembly of IDPs/IDRs allows reaching a large variety of nanostructures. Hybrid materials based on the derivatives of IDPs/IDRs show a promising performance as alternative biocides and nanodrugs. Cell mimicking, in vivo compartmentalization, and bone regeneration are demonstrated for IDPs/IDRs in biotechnological applications. The exciting possibilities of IDPs/IDRs in nanotechnology with relevant biological applications are shown.

Inhibition of tau amyloid formation and disruption of its preformed fibrils by Naphthoquinone-Dopamine hybrid

Misfolding and aggregation of tau protein, into pathological amyloids, are hallmarks of a group of neurodegenerative diseases collectively termed tauopathies and their modulation may be therapeutically valuable. Herein, we describe the synthesis and characterization of a dopamine-based hybrid molecule, naphthoquinone-dopamine (NQDA). Using thioflavin S assay, CD, transmission electron microscopy, dynamic light scattering, Congo Red birefringence, and large unilamellar vesicle leakage assays, we demonstrated its efficacy in inhibiting the in vitro aggregation of key tau-derived amyloidogenic fragments, PHF6 (VQIVYK) and PHF6* (VQIINK), prime drivers of aggregation of full-length tau in disease pathology. Isothermal titration calorimetry analysis revealed that the interaction between NQDA and PHF6 is spontaneous and has significant binding efficiency driven by both entropic and enthalpic processes. Furthermore, NQDA efficiently disassembled preformed fibrils of PHF6 and PHF6* into nontoxic species. Molecular dynamic simulations supported the in vitro results and provided a plausible mode of binding of NQDA with PHF6 fibril. NQDA was also capable of inhibiting the aggregation of full-length tau protein and disrupting its preformed fibrils in vitro in a dose-dependent manner. In a comparative study, the IC₅₀ value (50% inhibition of fibril formation) of NQDA in inhibiting the aggregation of PHF6 (25 µm) was ~ 17 µm, which is lower than for other bona fide amyloid inhibitors, naphthoquinone-tryptophan, rosmarinic acid, epigallocatechin gallate, ~ 21, ~ 77, or ~ 19 µm, respectively. Comparable superiority of NQDA was observed for inhibition of PHF6*. These findings suggest that NQDA can be a useful scaffold for designing new therapeutics for Alzheimer's disease and other tauopathies.

An α-helix mimetic oligopyridylamide, ADH-31, modulates Aβ42 monomer aggregation and destabilizes protofibril structures: insights from molecular dynamics simulations

Alzheimer's disease (AD), an epidemic growing worldwide due to no effective medical aid available in the market, is a neurological disorder. AD is known to be directly associated with the toxicity of amyloid-β (Aβ) aggregates. In search of potent inhibitors of Aβ aggregation, Hamilton and co-workers reported an α-helix mimetic, ADH-31, which acts as a powerful antagonist of Aβ42 aggregation. To identify the key interactions between protein-ligand complexes and to gain insights into the inhibitory mechanism of ADH-31 against Aβ42 aggregation, molecular dynamics (MD) simulations were performed in the present study. The MD simulations highlighted that ADH-31 showed distinct binding capabilities with residues spanning from the N-terminal to the central hydrophobic core (CHC) region of Aβ42 and restricted the conformational transition of the helix-rich structure of Aβ42 into another form of secondary structures (coil/turn/β-sheet). Hydrophobic contacts, hydrogen bonding and π-π interaction contribute to the strong binding between ADH-31 and Aβ42 monomer. The Dictionary of Secondary Structure of Proteins (DSSP) analysis highlighted that the probability of helical content increases from 38.5% to 50.2% and the turn content reduces from 14.7% to 6.2% with almost complete loss of the β-sheet structure (4.5% to 0%) in the Aβ42 monomer + ADH-31 complex. The per-residue binding free energy analysis demonstrated that Arg5, Tyr10, His14, Gln15, Lys16, Val18, Phe19 and Lys28 residues of Aβ42 are responsible for the favourable binding free energy in Aβ42 monomer + ADH-31 complex, which is consistent with the 2D HSQC NMR of the Aβ42 monomer that depicted a change in the chemical shift of residues spanning from Glu11 to Phe20 in the presence of ADH-31. The MD simulations highlighted the prevention of sampling of amyloidogenic β-strand conformations in Aβ42 trimer in the presence of ADH-31 as well as the ability of ADH-31 to destabilize Aβ42 trimer and protofibril structures. The lower binding affinity between Aβ42 trimer chains in the presence of ADH-31 highlights the destabilization of the Aβ42 trimer structure. Overall, MD results highlighted that ADH-31 inhibited Aβ42 aggregation by constraining Aβ peptides into helical conformation and destabilized Aβ42 trimer as well as protofibril structures. The present study provides a theoretical insight into the atomic level details of the inhibitory mechanism of ADH-31 against Aβ42 aggregation as well as protofibril destabilization and could be implemented in the structure-based drug design of potent therapeutic agents for AD.

DSS1 allosterically regulates the conformation of the tower domain of BRCA2 that has dsDNA binding specificity for homologous recombination

DSS1 is an evolutionary conserved, small intrinsically disordered protein that regulates various cellular functions. Although several studies have elucidated the role of DSS1 in stabilizing BRCA2 and its importance in homologous recombination repair (HRR), yet the structural mechanism behind the stability and HRR remains elusive. In this study, using molecular dynamics simulation we show that DSS1 stabilizes linearly arranged DNA/DSS1 binding domains of BRCA2 with many native contacts. These contacts are absent in the complexes with two missense DSS1 mutants associated with germline breast cancer and somatic mouth carcinoma. Most importantly, our protein energy-based network models show DSS1 allosterically regulates the conformation of the distant tower domain of BRCA2 that has dsDNA binding specificity for HRR. We further postulate that the unique conformation of the tower domain with kinked-helices might be responsible for DNA strand invasion and initiation of HRR. Induced conformation of the tower domain by the kinked-helices is absent in the unbound BRCA2, as well as in the two mutant DSS1-BRCA2 complexes. This suggests that DSS1 allosterically regulates the tower domain conformations of BRCA2 that affects dsDNA binding, essential for HRR. Our results add a new dimension to the function of DSS1 and its role in regulating HRR.

Purpurin modulates Tau-derived VQIVYK fibrillization and ameliorates Alzheimer's disease-like symptoms in animal model

Neurofibrillary tangles of the Tau protein and plaques of the amyloid β peptide are hallmarks of Alzheimer's disease (AD), which is characterized by the conversion of monomeric proteins/peptides into misfolded β-sheet rich fibrils. Halting the fibrillation process and disrupting the existing aggregates are key challenges for AD drug development. Previously, we performed in vitro high-throughput screening for the identification of potent inhibitors of Tau aggregation using a proxy model, a highly aggregation-prone hexapeptide fragment 306VQIVYK311 (termed PHF6) derived from Tau. Here we have characterized a hit molecule from that screen as a modulator of Tau aggregation using in vitro, in silico, and in vivo techniques. This molecule, an anthraquinone derivative named Purpurin, inhibited ~ 50% of PHF6 fibrillization in vitro at equimolar concentration and disassembled pre-formed PHF6 fibrils. In silico studies showed that Purpurin interacted with key residues of PHF6, which are responsible for maintaining its β-sheets conformation. Isothermal titration calorimetry and surface plasmon resonance experiments with PHF6 and full-length Tau (FL-Tau), respectively, indicated that Purpurin interacted with PHF6 predominantly via hydrophobic contacts and displayed a dose-dependent complexation with FL-Tau. Purpurin was non-toxic when fed to Drosophila and it significantly ameliorated the AD-related neurotoxic symptoms of transgenic flies expressing WT-FL human Tau (hTau) plausibly by inhibiting Tau accumulation and reducing Tau phosphorylation. Purpurin also reduced hTau accumulation in cell culture overexpressing hTau. Importantly, Purpurin efficiently crossed an in vitro human blood-brain barrier model. Our findings suggest that Purpurin could be a potential lead molecule for AD therapeutics.

Molecular insights into the inhibitory mechanisms of gallate moiety on the Aβ1-40 amyloid aggregation: A molecular dynamics simulation study

Alzheimer's disease is the most common form of neurodegenerative disease and the formation of Aβ amyloid aggregates has been widely demonstrated to be the principal cause of Alzheimer's disease. Our previous study and other studies suggested that the gallate moiety played an obligatory role in the inhibition process of naturally occurring polyphenols on Aβ amyloid fibrils formation. However, the detailed mechanisms were still unknown. Thus, in the present study, the gallic acid (GA) was specially selected and the molecular recognition mechanisms between GA molecules and Aβ_1-40 monomer were examined and analyzed by molecular dynamics simulation. The in silico experiments revealed that GA significantly prevented the conformational changes of Aβ_1-40 monomer with no β-sheet structure during the whole 100 ns. By analyzing the binding sites of GA molecules to Aβ_1-40 monomer, we found that both hydrophilic and hydrophobic amino acid residues were participated in the binding of GA molecules to Aβ_1-40 monomer. Moreover, results from the binding free energy analysis further demonstrated that the strength of polar interactions was significantly stronger than that of nonpolar interactions. We believed that our results could help to elucidate the underlying mechanisms of gallate moiety on the anti-amyloidogenic effects of polyphenols at the atomic level.

Modelling of interactions between Aβ(25-35) peptide and phospholipid bilayers: effects of cholesterol and lipid saturation

Aggregation of amyloid beta (Aβ) peptides in neuronal membranes is a known promoter of Alzheimer’s disease. To gain insight into the molecular details of Aβ peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aβ(25–35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most “active” amino-acid for peptide–lipid and peptide–cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer’s disease.

Molecular dynamics simulations of Aβ(25–35) peptides in phospholipid bilayers are carried out to investigate the effect of polyunsaturated lipids and cholesterol on aggregation of the peptides.

Pairwise-additive and polarizable atomistic force fields for molecular dynamics simulations of proteins

Protein force fields have been undergoing continual development since the first complete parameter sets were introduced nearly four decades ago. The functional forms that underlie these models have many common elements for the treatment of bonded and nonbonded forces, which are reviewed here. The most widely used force fields to date use a fixed-charge convention in which electronic polarization effects are treated via a mean-field approximation during partial charge assignment. Despite success in modeling folded proteins over many years, the fixed-charge assumption has limitations that cannot necessarily be overcome within their potential energy equations. To overcome these limitations, several force fields have recently been derived that explicitly treat electronic polarization effects with straightforward extensions of the potential energy functions used by nonpolarizable force fields. Here, we review the history of the most popular nonpolarizable force fields (AMBER, CHARMM, OPLS, and GROMOS) as well as studies that have validated them and applied them to studies of protein folding and misfolding. Building upon these force fields are more recent polarizable interaction potentials, including fluctuating charge models, POSSIM, AMOEBA, and the classical Drude oscillator. These force fields differ in their implementations but all attempt to model electronic polarization in a computationally tractable manner. Despite their recent emergence in the field of protein folding, several studies have already applied these polarizable models to challenging problems in this domain, including the role of polarization in folding free energies and sequence-specific effects on the stability of α-helical structures.

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.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.

Effect of lipid saturation on amyloid-beta peptide partitioning and aggregation in neuronal membranes: molecular dynamics simulations

Aggregation of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β (Aβ) peptides, cleaved from the amyloid precursor protein, is known as a precursor of the Alzheimer’s disease (AD). It is also known that Alzheimer’s disease is characterized by a substantial decrease of the amount of polyunsaturated lipids in the neuronal membranes of the frontal gray matter. To get insight into possible interconnection of these phenomena, we have carried out molecular dynamics simulations of two fragments of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptide, A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{1-28}$$\end{document}1-28 and A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{26-40}$$\end{document}26-40, in four different lipid bilayers: two monocomponent ones (14:0-14:0 PC, 18:0-22:6 PC), and two bilayers containing mixtures of 18:0-18:0 PE, 22:6-22:6 PE, 16:0-16:0 PC and 18:1-18:1 PC lipids of composition mimicking neuronal membranes in a “healthy” and “AD” brain. The simulations showed that the presence of lipids with highly unsaturated 22:6cis fatty acids chains strongly affects the interaction of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptides with lipid membranes. The polyunsaturated lipids cause stronger adsorption of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-peptides by the membrane and lead to weaker binding between peptides when the latter form aggregates. This difference in the behaviour observed in monocomponent bilayers is propagated in a similar fashion to the mixed membranes mimicking composition of neuronal membranes in “healthy” and “AD” brains, with “healthy” membrane having higher fraction of polyunsaturated lipids. Our simulations give strong indication that it can be physical–chemical background of the interconnection between amyloid fibrillization causing Alzheimer’s disease, and content of polyunsaturated lipids in the neuronal membranes.

Polyproline chains destabilize the Alzheimer's amyloid-β protofibrils: A molecular dynamics simulation study

Alzheimer's is a fatal neurodegenerative disease for which there is no cure at present. The disease is characterized by the presence of plaques, principally comprising the amyloid-β peptide (viz., β-sheet) in the brains of a patient. In our present work, we study the interaction of these β-sheets with a different number of repeating units of proline (β-sheet breaker) by docking and all atom molecular dynamics simulations. Our results indicate that proline can break the amyloid protofibrils apart, cause them to break their β-sheet structure, and in some cases even induce the formation of 3₁₀ helices, which may be intermediates in the unfolding of these β-sheets. We have also observed that some of the important hydrogen bonds and salt bridges between chains were disrupted by proline and the tight interatomic packing of atoms in the fibrils was made relatively loose. Proline chains had a tendency to make several contacts with charged residues. Proline chains binded well to the fibrils by strong electrostatic interactions while hydrophobic interactions played a less important role. This leads to the conclusion that proline can break the amyloid fibrils apart and can be considered in the design of novel peptide-based drugs to treat Alzheimer's disease and potentially other diseases caused by the misfolding of proteins into β-sheets.

Exploring Amyloid-β Dimer Structure Using Molecular Dynamics Simulations

A major hallmark of Alzheimer's disease (AD) is the aggregation of amyloid-β peptides in the brains of people afflicted by the disease. The exact pathway to this catastrophic event is unknown. In this work, a total of 9.5 μs molecular dynamics simulations have been performed to investigate the structure and dynamics of the smallest form of toxic Aβ oligomers, i.e., the Aβ dimers. This study suggests that specific hydrophobic regions are vital in the aggregation process. Different possible structures for Aβ dimers are reported along with their relative binding affinity. These data may be used to design better Aβ-aggregation inhibitors. The diversity of the dimer structures suggests several aggregation pathways.

Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles

The concept of intrinsically disordered proteins (IDPs) has emerged relatively slowly, but over the past 20 years, it has become an intense research area in structural biology. Indeed, because of their considerable flexibility and structural heterogeneity, the determination of IDP conformational ensemble is particularly challenging and often requires a combination of experimental measurements and computational approaches. With the improved accuracy of all-atom force fields and the increasing computing performances, molecular dynamics (MD) simulations have become more and more reliable to generate realistic conformational ensembles. And the combination of MD simulations with experimental approaches, such as nuclear magnetic resonance (NMR) and/or small-angle X-ray scattering (SAXS) allows one to converge toward a more accurate and exhaustive description of IDP structures. In this Review, we discuss the state of the art of MD simulations of IDP conformational ensembles, with a special focus on studies that back-calculated and directly compared theoretical and experimental NMR or SAXS observables, such as chemical shifts (CS), ³J-couplings (³Jc), residual dipolar couplings (RDC), or SAXS intensities. We organize the review in three parts. In the first section, we discuss the studies which used NMR and/or SAXS data to test and validate the development of force fields or enhanced sampling techniques. In the second part, we explore different methods for the refinement of MD-derived structural ensembles, such as NMR or SAXS data-restrained MD simulations or ensemble reweighting to better fit experiments. Finally, we survey some recent studies combining MD simulations with NMR and/or SAXS measurements to investigate the relationship between IDP conformational ensemble and biological activity, as well as their implication in human diseases. From this review, we noticed that quite a few studies compared MD-generated conformational ensembles with both NMR and SAXS measurements to validate IDP structures at both local and global levels. Yet, beside the IDP propensity to form local secondary structures, their dynamic extension or compactness also appears important for their activity. Thus, we believe that a close synergy between MD simulations, NMR, and SAXS experiments would be greatly appropriate to address the challenges of characterizing the disordered structures of proteins and their complexes, relative to their biological functions.

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ16-22 Dimer

We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.

Integrating in vitro and in silico approaches to evaluate the "dual functionality" of palmatine chloride in inhibiting and disassembling Tau-derived VQIVYK peptide fibrils

BACKGROUND

Alzheimer's disease (AD) is the most common neurodegenerative disorder which is characterized by the deposits of intra-cellular tau protein and extra-cellular amyloid-β (Aβ) peptides in the human brain. Understanding the mechanism of protein aggregation and finding compounds that are capable of inhibiting its aggregation is considered to be highly important for disease therapy.

METHODS

We used an in vitro High-Throughput Screening for the identification of potent inhibitors of tau aggregation using a proxy model; a highly aggregation-prone hexapeptide fragment ³⁰⁶VQIVYK³¹¹ derived from tau. Using ThS fluorescence assay we screened a library of 2401 FDA approved, bio-active and natural compounds in attempt to find molecules which can efficiently modulate tau aggregation.

RESULTS

Among the screened compounds, palmatine chloride (PC) alkaloid was able to dramatically reduce the aggregation propensity of PHF6 at sub-molar concentrations. PC was also able to disassemble preformed aggregates of PHF6 and reduce the amyloid content in a dose-dependent manner. Insights obtained from MD simulation showed that PC interacted with the key residues of PHF6 responsible for β-sheet formation, which could likely be the mechanism of inhibition and disassembly. Furthermore, PC could effectively inhibit the aggregation of full-length tau and disassemble preformed aggregates.

CONCLUSIONS

We found that PC possesses "dual functionality" towards PHF6 and full-length tau, i.e. inhibit their aggregation and disassemble pre-formed fibrils.

GENERAL SIGNIFICANCE

The "dual functionality" of PC is valuable as a disease modifying strategy for AD, and other tauopathies, by inhibiting their progress and reducing the effect of fibrils already present in the brain.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

Characterization of the structural ensembles of p53 TAD2 by molecular dynamics simulations with different force fields

The conformations of p53 TAD2 in complexes and sampled in simulations with five force fields.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

Envisaging the Structural Elevation in the Early Event of Oligomerization of Disordered Amyloid β Peptide

In Alzheimer's disease (AD), amyloid β (Aβ) protein plays a detrimental role in neuronal injury and death. Recent in vitro and in vivo studies suggest that soluble oligomers of the Aβ peptide are neurotoxic. Structural properties of the oligomeric assembly, however, are largely unknown. Our present investigation established that the 40-residue-long Aβ peptide (Aβ40) became more helical, ordered, and compact in the oligomeric state, and both the helical and β-sheet components were found to increase significantly in the early event of oligomerization. The band-selective two-dimensional NMR analysis suggested that majority of the residues from sequence 12 to 22 gained a higher-ordered secondary structure in the oligomeric condition. The presence of a significant amount of helical conformation was confirmed by Raman bands at 1650 and 1336 cm^-1. Other residues remained mostly in the extended polyproline II (PPII) and less compact β-conformation space. In the event of maturation of the oligomers into an amyloid fiber, both the helical content and the PPII-like structural components declined and ∼72% residues attained a compact β-sheet structure. Interestingly, however, some residues remained in the collagen triple helix/extended 2.5₁-helix conformation as evidenced by the amide III Raman signature band at 1272 cm^-1. Molecular dynamics analysis using an optimized potential for liquid simulation force field with the peptide monomer indicated that some of the residues may have preferences for helical conformation and this possibly contributed in the event of oligomer formation, which eventually became a β-sheet-rich amyloid fiber.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimer's disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimer's disease.

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.

The molecular behavior of a single β-amyloid inside a dipalmitoylphosphatidylcholine bilayer at three different temperatures: An atomistic simulation study: Aβ interaction with DPPC: Atomistic simulation

The behavior of a single Aβ40 molecule within a dipalmitoylphosphatidylcholine (DPPC) bilayer was studied by all-atom molecular dynamics simulations. The effect of membrane structure was investigated on Aβ40 behavior, secondary structure, and insertion depth. Simulations were performed at three temperatures (323, 310, and 300 K) to probe three different bilayer fluidities. Results show that at all above temperatures, the peptide contains two short helices, coil, bend, and turn structures. At 300 K, the peptide contains a region with β structure in C-terminal region. Our results also show that Aβ decreases the bilayer thickness and the order of lipids in its vicinity which leads to water insertion into the bilayer and concomitant increase in the local fluidity. The peptide remains embedded in the bilayer at all temperatures, and become inserted into the bilayer up to several residues at 323 and 310 K. At 310 and 300 K, the dominant interaction energy between Aβ and bilayer changes from electrostatic to van der Waals. It can be proposed that at higher temperatures (e.g., 323 K), Lys28 and the C-terminal region of the peptide play the role of two anchors that keep Aβ inside the top leaflet. This study demonstrates that Aβ molecule can perturb the integrity of cellular membranes. Proteins 2017; 85:1298-1310. © 2017 Wiley Periodicals, Inc.

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13-26) with hydrogen bonds and π-π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16-20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM-PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (-43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.