TY - JOUR

T1 - Extraction of kinetics from equilibrium distributions of states using the Metropolis Monte Carlo method

AU - Chekmarev, Sergei F.

N1 - Funding Information: I thank Irina Gopich for valuable discussions of the results of the FRET experiments for folding of protein, and Hoi Sung Chung for sharing the corresponding FRET-efficiency histograms. The work was supported by a contract with the IT SB RAS. Publisher Copyright: © 2022 American Physical Society.

PY - 2022/3

Y1 - 2022/3

N2 - The Metropolis Monte Carlo (MMC) method is used to extract reaction kinetics from a given equilibrium distribution of states of a complex system. The approach is illustrated by the folding/unfolding reaction for two proteins: a model ß-hairpin and a-helical protein a3D. For the ß-hairpin, the free energy surfaces (FESs) and free energy profiles (FEPs) are employed as the equilibrium distributions of states, playing a role of the potentials of mean force to determine the acceptance probabilities of new states in the MMC simulations. Based on the FESs and PESs for a set of temperatures that were simulated with the molecular dynamics (MD) method, the MMC simulations are performed to extract folding/unfolding rates. It has been found that the rate constants and first-passage time (FPT) distributions obtained in the MMC simulations change with temperature in good agreement with those from the MD simulations. For a3D protein, whose equilibrium folding/unfolding was studied with the single-molecule FRET method [Chung et al., J. Phys. Chem. A 115, 3642 (2011)1089-563910.1021/jp1009669], the FRET-efficiency histograms at different denaturant concentrations were used as the equilibrium distributions of protein states. It has been found that the rate constants for folding and unfolding obtained in the MMC simulations change with denaturant concentration in reasonable agreement with the constants that were extracted from the photon trajectories on the basis of theoretical models. The simulated FPT distributions are single-exponential, which is consistent with the assumption of two-state kinetics that was made in the theoretical models. The promising feature of the present approach is that it is based solely on the equilibrium distributions of states, without introducing any additional parameters to perform simulations, which suggests its applicability to other complex systems.

AB - The Metropolis Monte Carlo (MMC) method is used to extract reaction kinetics from a given equilibrium distribution of states of a complex system. The approach is illustrated by the folding/unfolding reaction for two proteins: a model ß-hairpin and a-helical protein a3D. For the ß-hairpin, the free energy surfaces (FESs) and free energy profiles (FEPs) are employed as the equilibrium distributions of states, playing a role of the potentials of mean force to determine the acceptance probabilities of new states in the MMC simulations. Based on the FESs and PESs for a set of temperatures that were simulated with the molecular dynamics (MD) method, the MMC simulations are performed to extract folding/unfolding rates. It has been found that the rate constants and first-passage time (FPT) distributions obtained in the MMC simulations change with temperature in good agreement with those from the MD simulations. For a3D protein, whose equilibrium folding/unfolding was studied with the single-molecule FRET method [Chung et al., J. Phys. Chem. A 115, 3642 (2011)1089-563910.1021/jp1009669], the FRET-efficiency histograms at different denaturant concentrations were used as the equilibrium distributions of protein states. It has been found that the rate constants for folding and unfolding obtained in the MMC simulations change with denaturant concentration in reasonable agreement with the constants that were extracted from the photon trajectories on the basis of theoretical models. The simulated FPT distributions are single-exponential, which is consistent with the assumption of two-state kinetics that was made in the theoretical models. The promising feature of the present approach is that it is based solely on the equilibrium distributions of states, without introducing any additional parameters to perform simulations, which suggests its applicability to other complex systems.

UR - http://www.scopus.com/inward/record.url?scp=85126701323&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.105.034407

DO - 10.1103/PhysRevE.105.034407

M3 - Article

C2 - 35428044

AN - SCOPUS:85126701323

VL - 105

JO - Physical Review E

JF - Physical Review E

SN - 2470-0045

IS - 3

M1 - 034407

ER -