Replica Exchange and Multicanonical Algorithms with the coarse-grained UNRES force field

Area of Science:

• Computational chemistry
• Biophysics
• Protein dynamics

Background:

• Accurate simulation of protein folding requires robust algorithms capable of exploring vast conformational space.
• Coarse-grained models, like UNRES, offer a computationally efficient approach to studying large biomolecules.

Purpose of the Study:

• To implement and compare three Replica Exchange (REM) based algorithms (REM, REMUCA, REMUCAREM) using the UNRES force field.
• To assess the performance of these algorithms in simulating the folding thermodynamics of a peptide and two small proteins.

Main Methods:

• Implementation of REM, REMUCA, and REMUCAREM algorithms in both Monte Carlo and Molecular Dynamics (MD) versions.
• Application of algorithms to poly-L-alanine (ala(20)), staphylococcal protein A B-domain (1BDD), and E. coli Mltd Lysm Domain (1E0G).
• Calculation of thermodynamic averages, heat capacity, and free energy maps (RMSD, radius of gyration) at various temperatures.

Main Results:

• All algorithms showed good agreement for poly-L-alanine, indicating correct implementation and equal effectiveness for small systems.
• REM demonstrated superior performance for the complex α+β protein (1E0G), yielding reliable thermodynamic statistics where others failed.
• Free energy calculations confirmed correct folding for ala(20) and 1BDD, while 1E0G showed native structure stability only at low temperatures due to higher entropy.

Conclusions:

• The choice of algorithm is crucial for accurate protein folding simulations, especially for complex structures.
• Replica Exchange (REM) is a reliable method for obtaining thermodynamic properties of proteins, particularly those with significant entropic contributions.
• Entropy plays a critical role in protein folding, influencing the difficulty of conformational space coverage and the stability of the native state.