Free Energy Computation by Monte Carlo Integration

TL;DR

This paper details a Monte Carlo integration method for computing protein-ligand binding free energies, emphasizing its simplicity, objectivity, and direct calculation advantages over perturbation methods.

Contribution

It introduces a Monte Carlo approach that reduces parameter dependence and assumptions, providing a more straightforward and deterministic way to calculate binding free energies.

Findings

Method yields gas-phase binding free energies.

Results are consistent with experimental data.

Minimal processing needed for free energy determination.

Abstract

The principles behind the computation of protein-ligand binding free energies by Monte Carlo integration are described in detail. The simulation provides gas-phase binding free energies that can be converted to aqueous energies by solvation corrections. The direct integration simulation has several characteristics beneficial to free-energy calculations. One is that the number of parameters that must be set for the simulation is small and can be determined objectively, making the outcome more deterministic, with respect to choice of input conditions, as compared to perturbation methods. Second, the simulation is free from assumptions about the starting pose or nature of the binding site. A final benefit is that binding free energies are a direct outcome of the simulation, and little processing is required to determine them. The well-studied T4 lysozyme experimental free energy data and crystal structures were used to evaluate the method.