Bromodomains,
protein domains involved in epigenetic regulation,
are able to bind small molecules with high affinity. In the present
study, we report free energy calculations for the binding of seven
ligands to the first BRD4 bromodomain, using the attach-pull-release
(APR) method to compute the reversible work of removing the ligands
from the binding site and then allowing the protein to relax conformationally.
We test three different water models, TIP3P, TIP4PEw, and SPC/E, as
well as the GAFF and GAFF2 parameter sets for the ligands. Our simulations
show that the apo crystal structure of BRD4 is only metastable, with
a structural transition happening in the absence of the ligand typically
after 20 ns of simulation. We compute the free energy change for this
transition with a separate APR calculation on the free protein and
include its contribution to the ligand binding free energies, which
generally causes an underestimation of the affinities. By testing
different water models and ligand parameters, we are also able to
assess their influence in our results and determine which one produces
the best agreement with the experimental data. Both free energies
associated with the conformational change and ligand binding are affected
by the choice of water model, with the two sets of ligand parameters
affecting their binding free energies to a lesser degree. Across all
six combinations of water model and ligand potential function, the
Pearson correlation coefficients between calculated and experimental
binding free energies range from 0.55 to 0.83, and the root-mean-square
errors range from 1.4–3.2 kcal/mol. The current protocol also
yields encouraging preliminary results when used to assess the relative
stability of ligand poses generated by docking or other methods, as
illustrated for two different ligands. Our method takes advantage
of the high performance provided by graphics processing units and
can readily be applied to other ligands as well as other protein systems.