scholarly journals Effects of Force Fields on the Conformational and Dynamic Properties of Amyloid β(1-40) Dimer Explored by Replica Exchange Molecular Dynamics Simulations

2017 ◽  
Author(s):  
Charles R. Watts ◽  
Andrew Gregory ◽  
Cole Frisbie ◽  
Sándor Lovas
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Dynamic Properties ◽  
Amyloid Β ◽  
Md Simulations ◽  
Conformational Space ◽  
Random Coil ◽  
Replica Exchange ◽  
Β Sheet

AbstractAlzheimer’s disease is histologically marked by fibrils of Amyloid beta (Aβ) peptide within the extracellular matrix. Fibrils themselves are benign compared to the cytotoxicity of the oligomers and pre-fibrillary aggregates. The conformational space and structural ensembles of Aβ peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1–40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER-ff99sb-ILDN, AMBER-ff99sb*-ILDN, AMBER-ff99sb-NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER-ff99sb-ILDN (8.8±6.4%) and CHARMM36 (2.7±4.2%), tend to overestimate the α-helical content compared to experimental CD (5.3±5.2%). Using the AMBER-ff99sb-NMR force field resulted in the greatest degree of variance (41.3±12.9%). Except for the AMBER-ff99sb-NMR force field, the others tended to under estimate the expected amount of β-sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER-ff99sb-NMR, reproduce a theoretically expected β-sheet-turn-β-sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C-terminal hydrophobic cores from residues 17-21 and 30-36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.

scholarly journals High-resolution mining of the SARS-CoV-2 main protease conformational space: supercomputer-driven unsupervised adaptive sampling

2021 ◽  
Author(s):  
Théo Jaffrelot Inizan ◽  
Frédéric Célerse ◽  
Olivier Adjoua ◽  
Dina El Ahdab ◽  
Luc-Henri Jolly ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
High Resolution ◽  
Adaptive Sampling ◽  
Force Fields ◽  
Sampling Strategy ◽  
Md Simulations ◽  
Conformational Space ◽  
Many Body ◽  
Polarizable Force Fields ◽  
Main Protease

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs).

scholarly journals Different force fields give rise to different amyloid aggregation pathways in molecular dynamics simulations

2020 ◽  
Author(s):  
Suman Samantray ◽  
Feng Yin ◽  
Batuhan Kav ◽  
Birgit Strodel
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Intrinsically Disordered Proteins ◽  
Force Fields ◽  
Md Simulations ◽  
Peptide Sequence ◽  
Disordered Proteins ◽  
Test Case ◽  
Peptide Aggregation ◽  
Intrinsically Disordered

AbstractThe progress towards understanding the molecular basis of Alzheimers’s disease is strongly connected to elucidating the early aggregation events of the amyloid-β (Aβ) peptide. Molecular dynamics (MD) simulations provide a viable technique to study the aggregation of Aβ into oligomers with high spatial and temporal resolution. However, the results of an MD simulation can only be as good as the underlying force field. A recent study by our group showed that none of the force fields tested can distinguish between aggregation-prone and non-aggregating peptide sequences, producing the same and in most cases too fast aggregation kinetics for all peptides. Since then, new force fields specially designed for intrinsically disordered proteins such as Aβ were developed. Here, we assess the applicability of these new force fields to studying peptide aggregation using the Aβ16−22 peptide and mutations of it as test case. We investigate their performance in modeling the monomeric state, the aggregation into oligomers, and the stability of the aggregation end product, i.e., the fibrillar state. A main finding is that changing the force field has a stronger effect on the simulated aggregation pathway than changing the peptide sequence. Also the new force fields are not able to reproduce the experimental aggregation propensity order of the peptides. Dissecting the various energy contributions shows that AMBER99SB-disp overestimates the interactions between the peptides and water, thereby inhibiting peptide aggregation. More promising results are obtained with CHARMM36m and especially its version with increased protein–water interactions. It is thus recommended to use this force field for peptide aggregation simulations and base future reparameterizations on it.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane

2021 ◽  
Author(s):  
Hebah Fatafta ◽  
Mohammed Khaled ◽  
Abdallah Sayyed-Ahmad ◽  
Birgit Strodel
Keyword(s):  
Aqueous Phase ◽  
Amyloid Β ◽  
Neuronal Cell ◽  
Md Simulations ◽  
Random Coil ◽  
Neuronal Membrane ◽  
Amyloid Β Peptide ◽  
Neuronal Membranes ◽  
Amyloid Oligomers ◽  
Β Sheet

The aggregation of amyloid β-peptides into neurotoxic oligomers is a key feature in the development of Alzheimer's disease. Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity. To gain a detailed understanding of the mutual effects of amyloid-β oligomers and the neuronal membrane, we carried out a total of 12 μs all-atom molecular dynamics (MD) simulations of the dimerization of the full-length Aβ42 peptide in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The conformational changes of Aβ42 resulting from its dimerization and interactions with the neuronal membrane are compared to those occurring upon its dimerization in the aqueous phase, which is also tested by 12 μs of MD simulations. We find that the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. Aβ42 dimerization in solution, on the other hand, is characterized by a random coil to β-sheet transition that seems to be on-pathway to amyloid aggregation. As the neurotoxic activity of amyloid oligomers increases with oligomer order, the results suggest that GM1 is neuroprotective against Aβ-mediated toxicity by inhibiting the formation of ordered amyloid oligomers.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Developing a molecular dynamics force field for both folded and disordered protein states

2018 ◽  
Vol 115 (21) ◽  
pp. E4758-E4766 ◽  
Author(s):  
Paul Robustelli ◽  
Stefano Piana ◽  
David E. Shaw
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Md Simulation ◽  
State Of The Art ◽  
Force Fields ◽  
Md Simulations ◽  
Disordered Proteins ◽  
Disordered Protein ◽  
Data Points ◽  
Folded Proteins

Molecular dynamics (MD) simulation is a valuable tool for characterizing the structural dynamics of folded proteins and should be similarly applicable to disordered proteins and proteins with both folded and disordered regions. It has been unclear, however, whether any physical model (force field) used in MD simulations accurately describes both folded and disordered proteins. Here, we select a benchmark set of 21 systems, including folded and disordered proteins, simulate these systems with six state-of-the-art force fields, and compare the results to over 9,000 available experimental data points. We find that none of the tested force fields simultaneously provided accurate descriptions of folded proteins, of the dimensions of disordered proteins, and of the secondary structure propensities of disordered proteins. Guided by simulation results on a subset of our benchmark, however, we modified parameters of one force field, achieving excellent agreement with experiment for disordered proteins, while maintaining state-of-the-art accuracy for folded proteins. The resulting force field, a99SB-disp, should thus greatly expand the range of biological systems amenable to MD simulation. A similar approach could be taken to improve other force fields.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

Sign in / Sign up

Export Citation Format

Share Document