A method for fast energy estimation and visualization of protein-ligand interaction

SummaryA new computational and graphical method for facilitating ligand-protein docking studies is developed on a three-dimensional computer graphics display. Various physical and chemical properties inside the ligand binding pocket of a receptor protein, whose structure is elucidated by X-ray crystal analysis, are calculated on three-dimensional grid points and are stored in advance. By utilizing those tabulated data, it is possible to estimate the non-bonded and electrostatic interaction energy and the number of possible hydrogen bonds between protein and ligand molecules in real time during an interactive docking operation. The method also provides a comprehensive visualization of the local environment inside the binding pocket.With this method, it becomes easier to find a roughly stable geometry of ligand molecules, and one can therefore make a rapid survey of the binding capability of many drug candidates. The method will be useful for drug design as well as for the examination of protein-ligand interactions.

[1]  M Karplus,et al.  Picosecond dynamics of tyrosine side chains in proteins. , 1979, Biochemistry.

[2]  J M Burridge,et al.  Computer graphics in drug design: molecular modeling of thyroid hormone-prealbumin interactions. , 1982, Journal of medicinal chemistry.

[3]  J M Blaney,et al.  Comparison of the inhibition of Escherichia coli and Lactobacillus casei dihydrofolate reductase by 2,4-diamino-5-(substituted-benzyl)pyrimidines: quantitative structure-activity relationships, X-ray crystallography, and computer graphics in structure-activity analysis. , 1982, Journal of medicinal chemistry.

[4]  J Deisenhofer,et al.  Structure of the complex formed by bovine trypsin and bovine pancreatic trypsin inhibitor. II. Crystallographic refinement at 1.9 A resolution. , 1974, Journal of molecular biology.

[5]  J M Blaney,et al.  The use of crystallography, graphics, and quantitative structure--activity relationships in the analysis of the papain hydrolysis of X-phenyl hippurates. , 1982, Archives of biochemistry and biophysics.

[6]  Thomas E. Ferrin,et al.  Computer graphics in real‐time docking with energy calculation and minimization , 1985 .

[7]  Arieh Warshel,et al.  Calculations of chemical processes in solutions , 1979 .

[8]  Andrew Dearing,et al.  Studies of nucleotide conformations and interactions. The relative stabilities of double‐helical B‐DNA sequence isomers , 1981, Biopolymers.

[9]  J M Burridge,et al.  Refined crystal structures of Escherichia coli and chicken liver dihydrofolate reductase containing bound trimethoprim. , 1985, The Journal of biological chemistry.

[10]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[11]  J M Burridge,et al.  Dihydrofolate reductase. The stereochemistry of inhibitor selectivity. , 1985, The Journal of biological chemistry.

[12]  Peter A. Kollman,et al.  Molecular mechanics simulation of protein-ligand interactions: binding of thyroid hormone analogs to prealbumin , 1982 .

[13]  J M Blaney,et al.  Electrostatic potential molecular surfaces. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[14]  T F Koetzle,et al.  The crystal and molecular structure of the antifolate drug trimethoprim (2,4-diamino-5-(3,4,5-trimethoxybenzyl)pyrimidine). A neutron diffraction study. , 1976, Journal of the American Chemical Society.

[15]  T. L. Blundell,et al.  DOCKER, an interactive program for simulating protein receptor and substrate interactions , 1983 .

[16]  J. Bolin,et al.  Crystal structures of Escherichia coli and Lactobacillus casei dihydrofolate reductase refined at 1.7 A resolution. I. General features and binding of methotrexate. , 1982, The Journal of biological chemistry.

[17]  R Balasubramanian,et al.  Potential functions for hydrogen bond interactions. I. A modified Lippincott-Schroeder potential function for NH, O interaction between peptide groups. , 1970, Biochimica et biophysica acta.

[18]  Jack D. Dunitz,et al.  Lone-pair directionality in hydrogen-bond potential functions for molecular mechanics calculations: the inhibition of human carbonic anhydrase II by sulfonamides , 1985 .

[19]  E. Lippincott,et al.  One‐Dimensional Model of the Hydrogen Bond , 1955 .

[20]  Nobuhiro Kosugi,et al.  Efficient methods for solving the open-shell scf problem and for obtaining an initial guess. The , 1980 .

[21]  Peter Murray-Rust,et al.  Directional hydrogen bonding to sp2- and sp3-hybridized oxygen atoms and its relevance to ligand-macromolecule interactions , 1984 .

[22]  Graham M. Smith,et al.  Intermolecular Modeling Methods in Drug Design/Modeling the Mechanism of Peptide Cleavage by Thermolysin , 1985, Annals of the New York Academy of Sciences.

[23]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[24]  L. Kuyper,et al.  Receptor-based design of dihydrofolate reductase inhibitors: comparison of crystallographically determined enzyme binding with enzyme affinity in a series of carboxy-substituted trimethoprim analogues. , 1982, Journal of medicinal chemistry.

[25]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[26]  J. Samama,et al.  Design of inhibitors from the three-dimensional structure of alcohol dehydrogenase. Chemical synthesis and enzymic properties , 1984 .

[27]  G M Crippen,et al.  General distance geometry three-dimensional receptor model for diverse dihydrofolate reductase inhibitors. , 1984, Journal of medicinal chemistry.

[28]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .