The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.

The calculation of binding affinities for flexible ligands has hitherto required the availability of reliable molecular mechanics parameters for the ligands, a restriction that can in principle be lifted by using a mixed quantum mechanics/molecular mechanics (QM/MM) representation in which the ligand is treated quantum mechanically. The feasibility of this approach is evaluated here, combining QM/MM with the Poisson-Boltzmann/surface area model of continuum solvation and testing the method on a set of 47 benzamidine derivatives binding to trypsin. The experimental range of the absolute binding energy (DeltaG = -3.9 to -7.6 kcal/mol) is reproduced well, with a root-mean-square (RMS) error of 1.2 kcal/mol. When QM/MM is applied without reoptimization to the very different ligands of FK506 binding protein the RMS error is only 0.7 kcal/mol. The results show that QM/MM is a promising new avenue for automated docking and scoring of flexible ligands. Suggestions are made for further improvements in accuracy.

The molecular motor myosin converts chemical energy from ATP hydrolysis into mechanical work, thus driving a variety of essential motility processes. Although myosin function has been studied extensively, the catalytic mechanism of ATP hydrolysis and its chemomechanical coupling to the motor cycle are not completely understood. Here, the catalysis mechanism in myosin II is examined using quantum mechanical/molecular mechanical reaction path calculations. The resulting reaction pathways, found in the catalytically competent closed/closed conformation of the Switch-1/Switch-2 loops of myosin, are all associative with a pentavalent bipyramidal oxyphosphorane transition state but can vary in the activation mechanism of the attacking water molecule and in the way the hydrogens are transferred between the heavy atoms. The coordination bond between the Mg2+ metal cofactor and Ser237 in the Switch-1 loop is broken in the product state, thereby facilitating the opening of the Switch-1 loop after hydrolysis is completed, which is required for subsequent strong rebinding to actin. This reveals a key element of the chemomechanical coupling that underlies the motor cycle, namely, the modulation of actin unbinding or binding in response to the ATP or ADP x P(i) state of nucleotide-bound myosin.

No abstract