Why are MD simulated protein folding times wrong?

Dmitry Nerukhdn1

  • 1Unilever Centre for Molecular Sciences Informatics, Department of Chemistry, Cambridge University, Cambridge CB2 1EW, UK. 232@cam.ac.uk

Insights

Protein folding times from molecular dynamics simulations can significantly deviate from experimental values. This sensitivity to simulation parameters like forcefield and temperature means results must be interpreted cautiously.

Area of Science:

  • Computational biology
  • Biophysics
  • Molecular dynamics simulations

Background:

  • Accurate simulation of protein folding is crucial for understanding biological function.
  • Molecular dynamics (MD) simulations are a key tool for studying protein folding dynamics.

Purpose of the Study:

  • To investigate the causes of significant deviations between simulated and experimental protein folding times.
  • To assess the sensitivity of protein folding time predictions to simulation model parameters.

Main Methods:

  • Utilized Markov State Models (MSM) to describe conformational dynamics of peptides and proteins.
  • Analyzed the impact of varying simulation parameters, including forcefield and temperature, on folding times.
  • Examined two specific peptide systems to illustrate observed deviations.

Main Results:

  • Protein folding times simulated using MSMs are highly sensitive to model parameters.
  • Modest variations in forcefield and temperature can lead to order-of-magnitude differences in folding times.
  • Observed significant discrepancies between simulated and experimental folding rates for the studied peptides.

Conclusions:

  • Simulated protein folding rates derived from molecular dynamics require careful interpretation due to parameter sensitivity.
  • The choice of forcefield and temperature critically influences the accuracy of simulated folding times.
  • Further refinement of simulation methodologies is needed to improve the reliability of folding time predictions.

Related Concept Videos

Protein Folding01:25

Protein Folding

Proteins are chains of amino acids linked together by peptide bonds. Upon synthesis, a protein folds into a three-dimensional conformation, critical to its biological function. Interactions between its constituent amino acids guide protein folding, and hence the protein structure is primarily dependent on its amino acid sequence.
Protein Structure Is Critical to Its Biological Function
Proteins perform a wide range of biological functions such as catalyzing chemical reactions, providing...
Protein Folding01:22

Protein Folding

Overview
Protein Folding01:22

Protein Folding

Overview
Molecular Chaperones and Protein Folding03:00

Molecular Chaperones and Protein Folding

The native conformation of a protein is formed by interactions between the side chains of its constituent amino acids. When the amino acids cannot form these interactions, the protein cannot fold by itself and needs chaperones. Notably, chaperones do not relay any additional information required for the folding of polypeptides; the native conformation of a protein is determined solely by its amino acid sequence. Chaperones catalyze protein folding without being a part of the folded protein.
The...
Amyloid Fibrils03:03

Amyloid Fibrils

Amyloid fibrils are aggregates of misfolded proteins.  Under most circumstances, misfolded proteins are either refolded by chaperone proteins or degraded by the proteasome. However, in the case of a mutation or a disease, these proteins can accumulate to form large clusters and often further assemble to form elongated fibers, called fibrils. 
Amyloid deposits were observed as early as 1639 in the liver and the spleen.   In 1854, Rudolph Virchow performed iodine staining, normally used to...
Improving Translational Accuracy02:07

Improving Translational Accuracy

Base complementarity between the three base pairs of mRNA codon and the tRNA anticodon is not a failsafe mechanism. Inaccuracies can range from a single mismatch to no correct base pairing at all. The free energy difference between the correct and nearly correct base pairs can be as small as 3 kcal/ mol. With complementarity being the only proofreading step, the estimated error frequency would be one wrong amino acid in every 100 amino acids incorporated. However, error frequencies observed in...