Summary
Derivatives of free energy differences have been calculated by molecular dynamics techniques. The systems under study were ternary complexes of Trimethoprim (TMP) with dihydrofolate reductases of E. coli and chicken liver, containing the cofactor NADPH. Derivatives are taken with respect to modification of TMP, with emphasis on altering the 3-, 4- and 5-substituents of the phenyl ring. A linear approximation allows the encompassing of a whole set of modifications in a single simulation, as opposed to a full perturbation calculation, which requires a separate simulation for each modification. In the case considered here, the proposed technique requires a factor of 1000 less computing effort than a full free energy perturbation calculation. For the linear approximation to yield a significant result, one has to find ways of choosing the perturbation evolution, such that the initial trend mirrors the full calculation. The generation of new atoms requires a careful treatment of the singular terms in the non-bonded interaction. The result can be represented by maps of the changed molecule, which indicate whether complex formation is favoured under movement of partial charges and change in atom polarizabilities. Comparison with experimental measurements of inhibition constants reveals fair agreement in the range of values covered. However, detailed comparison fails to show a significant correlation. Possible reasons for the most pronounced deviations are given.
Similar content being viewed by others
References
Beveridge, D.L. and DiCapua, F.M., Annu. Rev. Biophys. Biophys. Chem., 18 (1989) 431.
Wong, C.F., J. Am. Chem. Soc., 113 (1991) 3208.
Gao, J., Kuczera, A., Tidor, B. and Karplus, M., Science, 244 (1989) 1069.
Cross, A.J., Chem. Phys. Lett., 128 (1986) 198.
McDonald, J.J. and Brooks III, C.L., J. Am. Chem. Soc., 114 (1992) 2062.
Matthews, D.A., Bolin, J.T., Burridge, J.M., Filman, D.J., Volz, K.W., Kaufman, B.T. Bedell, C.R., Champness, J.N., Stammers, D.K. and Kraut, J., J. Biol. Chem., 260 (1985) 381.
Selassie, C.D., Fang, Z-X., Li, R-L., Hansch, C., Debnath, G., Klein, T.E., Langridge, R. and Kaufman, B.T., J. Med. Chem., 32 (1989) 1895.
Selassie, C.D., Li, R-L., Poe, M. and Hansch, C., J. Med. Chem., 34 (1991) 46.
Bystroff, C., Oatley, S.J. and Kraut, J., Biochemistry, 29 (1990) 3263.
van Gunsteren, W.F. and Berendsen, H.J.C., Groningen Molecular Similation (GROMOS) Library Manual, Biomos, Nijenborgh 16, 9747 AG Groningen, The Netherlands, 1987.
Ryckaert, J.-P., Ciccotti, G. and Berendsen, H.J.C., J. Comput. Phys., 2 (1977) 327.
vanGunsteren, W.F. and Berendsen, H.J.C., Angew. Chem. Int. Ed. Engl., 29 (1990) 992.
Berendsen, H.J.C., Postma, J.P.M., vanGunsteren, W.F., DiNola, A. and Haak, J.R., J. Che Phys., 92 (1984) 3684.
Berendsen, H.J.C., Postma, J.P.M., vanGunsteren, W.F. and Hermans, J., In Pullman, B., (Ed.) Intermolecular Forces. Reidel, Dordrecht, 1981, pp. 331–342.
Gasteiger, J. and Marsili, M., Tetrahedron, 36 (1980) 3219.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Gerber, P.R., Mark, A.E. & van Gunsteren, W.F. An approximate but efficient method to calculate free energy trends by computer simulation: Application to dihydrofolate reductase-inhibitor complexes. J Computer-Aided Mol Des 7, 305–323 (1993). https://doi.org/10.1007/BF00125505
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF00125505