Monte Carlo simulations of polyalanine using a reduced model and statistics-based interaction potentials

Area of Science:

• Computational Biology
• Biophysics
• Protein Science

Background:

• Understanding protein folding, specifically the coil-to-helix transition, is crucial for predicting protein structure and function.
• Existing models often require significant computational resources for accurate simulations.

Purpose of the Study:

• To develop and validate a coarse-grained model for simulating the coil-to-helix transition in polyalanine.
• To analyze the thermodynamics of this transition using a novel residue-residue interaction potential.

Main Methods:

• Utilized a coarse-grained residue-residue interaction potential derived from Protein Data Bank statistical analysis.
• Incorporated interaction potentials dependent on radial distance and relative orientation of interaction sites.
• Modeled two types of interaction sites: side chain and virtual backbone (for hydrogen bonding).
• Studied two polyalanine chain lengths and analyzed thermodynamic results using the Zimm-Bragg model.

Main Results:

• The coarse-grained model successfully investigated the coil-to-helix transition for polyalanine.
• Thermodynamic results of the transition were analyzed and showed good agreement with the Zimm-Bragg model.
• Model predictions were quantitatively and qualitatively consistent with all-atom and other reduced-model Monte Carlo simulations.

Conclusions:

• The developed coarse-grained potential provides a computationally efficient yet accurate method for studying protein secondary structure transitions.
• This reduced model effectively captures the essential physics governing the coil-to-helix transition in polyalanine.
• The approach holds promise for broader applications in simulating protein folding dynamics.