Free energy for protein folding from nonequilibrium simulations using the Jarzynski equality

Abstract

The equilibrium free energy difference between two long-lived molecular species or "conformational states" of a protein (or any other molecule) can in principle be estimated by measuring the work needed to shuttle the system between them, independent of the irreversibility of the process. This is the meaning of the Jarzynski equality (JE), which we test in this paper by performing simulations that unfold a protein by pulling two atoms apart. Pulling is performed fast relative to the relaxation time of the molecule and is thus far from equilibrium. Choosing a simple protein model for which we can independently compute its equilibrium properties, we show that the free energy can be exactly and effectively estimated from nonequilibrium simulations. To do so, one must carefully and correctly determine the ensemble of states that are pulled, which is more important the farther from equilibrium one performs simulations; this highlights a potential problem in using the JE to extract the free energy from forced unfolding experiments. The results presented here also demonstrate that the free energy difference between the native and denatured states of a protein measured in solution is not always equal to the free energy profile that can be estimated from forced unfolding simulations (or experiments) using the JE.