Determining equilibrium constants for dimerization reactions from molecular dynamics simulations

Abstract

With today's available computer power, free energy calculations from equilibrium molecular dynamics simulations "via counting" become feasible for an increasing number of reactions. An example is the dimerization reaction of transmembrane alpha-helices. If an extended simulation of the two helices covers sufficiently many dimerization and dissociation events, their binding free energy is readily derived from the fraction of time during which the two helices are observed in dimeric form. Exactly how the correct value for the free energy is to be calculated, however, is unclear, and indeed several different and contradictory approaches have been used. In particular, results obtained via Boltzmann statistics differ from those determined via the law of mass action. Here, we develop a theory that resolves this discrepancy. We show that for simulation systems containing two molecules, the dimerization free energy is given by a formula of the form G ∝ ln(P 1 /P 0 ). Our theory is also applicable to high concentrations that typically have to be used in molecular dynamics simulations to keep the simulation system small, where the textbook dilute approximations fail. It also covers simulations with an arbitrary number of monomers and dimers and provides rigorous error estimates. Comparison with test simulations of a simple Lennard Jones system with various particle numbers as well as with reference free energy values obtained from radial distribution functions show full agreement for both binding free energies and dimerization statistics.

Insights into Noncovalent Binding Obtained from Molecular Dynamics Simulations

“…The association constant K = exp(–Δ G / RT ) was obtained from integrating the PMF . The results are in qualitative agreement with experimental data while the magnitude of the aggregation strength is somewhat overpredicted by the simulations.…”

“…In previous studies , , umbrella sampling simulations were employed for the physical path and a combination of thermodynamic integration and enveloping distribution sampling for the alchemical path and their equivalency for binding of primary alcohols to α ‐cyclodextrin was shown. For the present work, this analysis was extended by direct counting to obtain the standard binding free enthalpy according to …”

Dynamic Cholesterol-Conditioned Dimerization of the G Protein Coupled Chemokine Receptor Type 4

“…The dissociation constants K D and the corresponding binding free energies Δ G were approximated according to [66]: c ⊘ is the standard concentration of 1 mol/l, N Av the Avogadro constant, and V the volume of the lipid-protein double layer. P 0,1 are the fractions of simulation time in monomeric and dimeric states, respectively.…”

Using Interpolation for Fast and Accurate Calculation of Ion–Ion Interactions

“…This, however, requires running long explicit-solvent MD simulations for a given ion pair in order to obtain reasonable statistics. 65−70 An alternative is to calculate the association constants from the radial distribution function. Both approaches imply that a definition of associated/dissociated state is given.…”