Studying the unfolding kinetics of proteins under pressure using long molecular dynamic simulation runs

The usefulness of computational methods such as molecular dynamics simulation has been extensively established for studying systems in equilibrium. Nevertheless, its application to complex non-equilibrium biological processes such as protein unfolding has been generally regarded as producing results which cannot be interpreted straightforwardly. In the present study, we present results for the kinetics of unfolding of apomyoglobin, based on the analysis of long simulation runs of this protein in solution at 3 kbar (1 atm = 1.01325, bar = 101,325 Pa). We hereby demonstrate that the analysis of the data collected within a simulated time span of 0.18 mus suffices for producing results, which coincide remarkably with the available unfolding kinetics experimental data. This not only validates molecular dynamics simulation as a valuable alternative for studying non-equilibrium processes, but also enables a detailed analysis of the actual structural mechanism which underlies the unfolding process of proteins under elusive denaturing conditions such as high pressure.