Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Area of Science:

• Computational Biology
• Biophysics
• Molecular Dynamics

Background:

• Protein folding is crucial for biological function but computationally challenging.
• Accurate calculation of protein folding energetics is vital for drug design and genetic interpretation.
• First-principles calculation of protein folding state functions remains an open problem.

Purpose of the Study:

• To present a simple yet accurate method for calculating protein folding energetics.
• To enable precise determination of conformational stability (ΔG).
• To validate the approach on various protein types and conditions.

Main Methods:

• Employing molecular dynamics simulations to compute energies of folded and unfolded protein states at varying temperatures.
• Deriving key thermodynamic quantities, enthalpy change (ΔH) and heat capacity change (ΔCp).
• Utilizing these quantities to calculate the Gibbs free energy of folding (ΔG).

Main Results:

• Successfully calculated the energetics for two- and three-state proteins across major structural classes.
• Accurately determined small differences in protein stability (ΔΔG) under altered solution conditions.
• Demonstrated the method's efficacy in assessing effects of single amino acid residue variations.

Conclusions:

• The presented method provides an accurate and accessible approach to protein folding energetics.
• This technique can significantly advance protein and drug design.
• It offers a valuable tool for interpreting genetic variations affecting protein stability.