Free energy calculations simulation methods

We assume that the unbinding reaction takes place on a time scale long ( ompared to the relaxation times of all other degrees of freedom of the system, so that the friction coefficient can be considered independent of time. This condition is difficult to satisfy on the time scales achievable in MD simulations. It is, however, the most favorable case for the reconstruction of energy landscapes without the assumption of thermodynamic reversibility, which is central in the majority of established methods for calculating free energies from simulations (McCammon and Harvey, 1987 Elber, 1996) (for applications and discussion of free energy calculation methods see also the chapters by Helms and McCammon, Hermans et al., and Mark et al. in this volume). 

For reasons of space and because of their prime importance, we focus here on free energy calculations based on detailed molecular dynamics (MD) or Monte Carlo (MC) simulations. However, several other computational approaches exist to calculate free energies, including continuum dielectric models and integral equation methods [4,14]. 

The titration coordinates evolve along with the dynamics of the conformational degrees of freedom, r, in simulations with GB implicit solvent models [37, 57], An extended Hamiltonian formalism, in analogy to the A dynamics technique developed for free energy calculations [50], is used to propagate the titration coordinates. The deprotonated and protonated states are those, for which the A value is approximately 1 or 0 (end-point states), respectively. Thus, in contrast to the acidostat method, where A represents the extent of deprotonation, is estimated from the relative occupancy of the states with A 1 (see later discussions). The extended Hamiltonian in the CPHMD method is a sum of the following terms [42],... 

Oberhofer, H. Dellago, C. Geissler, P. L., Biased sampling of nonequilibrium trajectories. Can fast switching simulations outperform conventional free energy calculation methods, J. Phys. Chem. B 2005,109, 6902-6915... 

Understanding and Improving Free Energy Calculations in Molecular Simulations Error Analysis and Reduction Methods... 

We can expect to see future research directed at QM/MM and ab initio simulation methods to handle these electronic structure effects coupled with path integral or approximate quantum free energy methods to treat nuclear quantum effects. These topics are broadly reviewed in [32], Nuclear quantum effects for the proton in water have already received some attention [30, 76, 77]. Utilizing the various methods briefly described above (and other related approaches), free energy calculations have been performed for a wide range of problems involving proton motion [30, 67-69, 71, 72, 78-80]. 

---