Abstract
A simulated annealing method for finding important ligand fragments is described. At a given temperature, ligand fragments are randomly selected and randomly placed within the given receptor cavity, often replacing or forming bonds with existing ligand fragments. For each new ligand fragment combination, the bonded, nonbonded, polarization and solvation energies of the new ligand–receptor system are compared to the previous configuration. Acceptance or rejection of the new system is decided using the Boltzmann distribution\({\text{e}}^{{\text{ - E/kT}}}\), where E is the energy difference between the old and new systems, k is the Boltzmann constant and T is the temperature. Thus, energetically unfavorable fragment switches are sometimes accepted, sacrificing immediate energy gains in the interest of finding a system with minimum energy. By lowering the temperature, the rate of unfavorable switches decreases and energetically favorable combinations become more difficult to change. The process is terminated when the frequency of switches becomes too small. As a test, the method predicted positions and types of important ligand fragments for neuraminidase that were in accord with the known ligand, sialic acid.
Similar content being viewed by others
References
Andrews, P., In Foye, W. (Ed.) Principles of Medicinal Chemistry, 3rd ed., Lea & Febiger, Philadelphia, PA, U.S.A., 1990, pp. 855– 860.
Bugg, C., Carson, W. and Montgomery, J., Sci. Am., 269 (1993) 92.
Von Itzstein, M., Wu, W., Kok, G., Pegg, M., Dyason, J., Jin, B., Pahn, T., Smythe, M., White, H., Oliver, S., Colman, P., Varghese, J., Ryan, M., Woods, J., Bethell, R., Hotham, V., Cameron, J. and Penn, C., Nature, 363 (1993) 418.
Kuntz, I., Blaney, J., Oatley, S., Langridge, R. and Ferrin, T., J. Mol. Biol., 161 (1982) 269.
DesJarlais, R., Sheridan, R., Seibel, J., Dixon, S., Kuntz, I. and Venkataraghavan, R., J. Med. Chem., 31 (1988) 722.
Martin, Y., J. Med. Chem., 35 (1992) 2145.
Rotstein, S. and Murcko, M., J. Med. Chem., 36 (1993) 1700.
Namasivayam, S. and Dean, P., J. Mol. Graph., 4 (1986) 46.
Lewis, R., Roe, D., Huang, C., Ferrin, T., Langridge, R. and Kuntz, I., J. Mol. Graph., 10 (1992) 66.
Kato, Y., Akiko, I. and Iitaka, Y., Tetrahedron, 43 (1987) 5229.
Meng, E., Shoichet, B. and Kuntz, I., J. Comput. Chem., 13 (1992) 505.
Karfunkel, H., J. Comput. Chem., 7 (1986) 113.
DesJarlais, R., Sheridan, R., Dixon, S., Kuntz, I. and Venkataraghavan, R., J. Med. Chem., 29 (1986) 2149.
Billeter, M., Havel, T. and Kuntz, I., Biopolymers, 26 (1987) 777.
Goodsell, D. and Olson, A., Proteins, 8 (1990) 195.
Hart, T. and Read, R., Proteins, 13 (1992) 206.
Goodford, P., J. Med. Chem., 28 (1985) 849.
Boobbyer, D., Goodford, P., McWhinni, P. and Wade, R., J. Med. Chem., 32 (1989) 1083.
Wade, R., Clark, K. and Goodford, P., J. Med. Chem., 36 (1993) 140.
Wade, R. and Goodford, P., J. Med. Chem., 36 (1993) 148.
Kellogg, G. and Abraham, D., J. Mol. Graph., 10 (1992) 212.
Miranker, A. and Karplus, M., Proteins, 1 (1991) 29.
Caflisch, A., Miranker, A. and Karplus, M., J. Med. Chem., 35 (1993) 2142.
Gillet, V., Johnson, P., Mata, P., Sike, S. and Williams, P., J. Comput.-Aided Mol. Design, 7 (1993) 127.
Nishibata, Y. and Akiko, I., Tetrahedron, 47 (1991) 8985.
Nishibata, Y. and Akiko, I., J. Med. Chem., 36 (1993) 2921.
Lawrence, M. and Davis, P., Proteins Struct. Funct. Genet., 12 (1992) 31.
Böhm, H.-J., J. Comput.-Aided Mol. Design, 6 (1992) 61.
Eisen, M., Wiley, D., Karplus, M. and Hubbard, R., Proteins Struct. Funct. Genet., 19 (1994) 199.
Lewis, R., J. Comput.-Aided Mol. Design, 4 (1990) 205.
Lewis, R.A. and Dean, P.M., Proc. R. Soc. London, B236 (1989) 125.
Singh, J., Saldanha, J. and Thornton, J., Protein Eng., 4 (1991) 251.
Moon, J. and Howe, W., Proteins Struct. Funct. Genet., 11 (1991) 314.
Rotstein, S. and Murcko, M., J. Comput.-Aided Mol. Design, 7 (1993) 23.
Pickett, S.D. and Sternberg, M.J.E., J. Mol. Biol., 231 (1993) 825.
Zielinski, P.J., An Annealing Algorithm for Designing Ligands from Receptor Structures, UMI, Ann Arbor, MI, U.S.A., 1994, p. 1.
Mayo, S., Olafson, B. and Goddard, W., J. Phys. Chem., 94 (1990) 8897.
Cramer, C. and Truhlar, D., J. Am. Chem. Soc., 113 (1991) 8305.
Wang, H. and Levinthal, C., J. Comput. Chem., 12 (1991) 868.
Cramer, C. and Truhlar, D., J. Am. Chem. Soc., 113 (1991) 8552.
Cramer, C. and Truhlar, D., Chem. Phys. Lett., 198 (1992) 74.
Cramer, C. and Truhlar, D., J. Comput. Chem., 13 (1992) 1089.
Cramer, C. and Truhlar, D., J. Comput.-Aided Mol. Design, 6 (1992) 629.
Cramer, C. and Truhlar, D., Science, 256 (1992) 213.
Cramer, C. and Truhlar, D., J. Am. Chem. Soc., 114 (1992) 8226.
Allen, L., J. Am. Chem. Soc., 111 (1989) 9115.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Burt, S., Hutchins, C. & Zielinski, P.J. A Monte Carlo method for finding important ligand fragments from receptor data. J Comput Aided Mol Des 11, 243–255 (1997). https://doi.org/10.1023/A:1007952511172
Issue Date:
DOI: https://doi.org/10.1023/A:1007952511172