Binding affinities for sulfonamide inhibitors with human thrombin using Monte Carlo simulations with a linear response method.

The binding of sulfonamide inhibitors to human thrombin is examined to evaluate the viability of calculating free energies of binding, deltaGb, utilizing Monte Carlo (MC) statistical mechanics with a linear response approach. Coulombic and van der Waals energy components determined from MC simulations of the bound and unbound inhibitors solvated in water plus a solvent-accessible surface area term, as an index for cavity formation, were correlated with the free energies of binding for the inhibitor MD-805 and six derivatives. The best correlations yield an average error of 0.8 kcal/mol for the seven binding affinities, which cover an observed range of 6.0 kcal/mol. The MC simulations also provided insights into the interactions occurring in the active site and the origins of variations in deltaGb. Equatorial placement of the carboxylate group at C2 in the piperidine ring of the inhibitors causes electrostatic destabilization with the side chain of Glu-H192, while axial disposition of the C4-methyl group reduces favorable hydrophobic interactions in the P-pocket of the enzyme.

[1]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[2]  T. Tilley,et al.  Reversible Cycloaddition of Isocyanates to Ruthenium Silylene Complexes , 1997 .

[3]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[4]  R. Ornstein,et al.  Binding free energy calculations for P450cam-substrate complexes. , 1996, Protein engineering.

[5]  T. Hansson,et al.  On the Validity of Electrostatic Linear Response in Polar Solvents , 1996 .

[6]  M Karplus,et al.  Functionality map analysis of the active site cleft of human thrombin , 1996, J. Comput. Aided Mol. Des..

[7]  T. Hansson,et al.  Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations. , 1995, Protein engineering.

[8]  E. Purisima,et al.  Calculation of relative binding free energies and configurational entropies: a structural and thermodynamic analysis of the nature of non-polar binding of thrombin inhibitors based on hirudin55-65. , 1995, Journal of molecular biology.

[9]  W. Leitner,et al.  σ Metathesis as a Critical Step for the Transition Metal Catalyzed Formation of Formic Acid from CO2 and H2? An Ab Initio Investigation , 1995 .

[10]  F. Diederich,et al.  Design of Novel, Nonpeptidic Thrombin Inhibitors and Structure of a Thrombin–Inhibitor Complex , 1995 .

[11]  E. Iwanowicz,et al.  Molecular modeling studies of novel retro-binding tripeptide active-site inhibitors of thrombin. , 1995, Bioorganic & medicinal chemistry.

[12]  E. Purisima,et al.  Interactions of hirudin-based inhibitor with thrombin: critical role of the IleH59 side chain of the inhibitor. , 1995, Biochemistry.

[13]  A. Tulinsky,et al.  Active-site mimetic inhibition of thrombin. , 1995, Acta crystallographica. Section D, Biological crystallography.

[14]  P. Grootenhuis,et al.  Correlation of binding affinities with non-bonded interaction energies of thrombin-inhibitor complexes. , 1995, Acta crystallographica. Section D, Biological crystallography.

[15]  K. Merz,et al.  Application of the free energy perturbation method to human carbonic anhydrase II inhibitors. , 1995, Journal of medicinal chemistry.

[16]  W. L. Jorgensen,et al.  AN EXTENDED LINEAR RESPONSE METHOD FOR DETERMINING FREE ENERGIES OF HYDRATION , 1995 .

[17]  Sherry L. Mowbray,et al.  Sugar Recognition by a Glucose/Galactose Receptor , 1995, The Journal of Biological Chemistry.

[18]  William L. Jorgensen,et al.  Free Energies of Hydration and Pure Liquid Properties of Hydrocarbons from the OPLS All-Atom Model , 1994 .

[19]  K Gubernator,et al.  Design and synthesis of potent and highly selective thrombin inhibitors. , 1994, Journal of medicinal chemistry.

[20]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[21]  Y Konishi,et al.  Crystal structure of the complex of human α‐thrombin and nonhydrolyzable bifunctional inhibitors, hirutonin—2 and hirutonin—6 , 1993, Proteins.

[22]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[23]  William L. Jorgensen,et al.  Monte Carlo simulations of the hydration of substituted benzenes with OPLS potential functions , 1993, J. Comput. Chem..

[24]  D. Turk,et al.  Crystallographic determination of thrombin complexes with small synthetic inhibitors as a starting point for the receptor-based design of antithrombotics. , 1993, Seminars in thrombosis and hemostasis.

[25]  H Brandstetter,et al.  Refined 2.3 A X-ray crystal structure of bovine thrombin complexes formed with the benzamidine and arginine-based thrombin inhibitors NAPAP, 4-TAPAP and MQPA. A starting point for improving antithrombotics. , 1992, Journal of molecular biology.

[26]  D. Banner,et al.  Crystallographic analysis at 3.0-A resolution of the binding to human thrombin of four active site-directed inhibitors. , 1994, The Journal of biological chemistry.

[27]  C. Esmon,et al.  The region of the thrombin receptor resembling hirudin binds to thrombin and alters enzyme specificity. , 1991, The Journal of biological chemistry.

[28]  V. Wheaton,et al.  Molecular cloning of a functional thrombin receptor reveals a novel proteolytic mechanism of receptor activation , 1991, Cell.

[29]  J. Maraganore,et al.  Thrombin structure and function: why thrombin is the primary target for antithrombotics. , 1991, Blood coagulation & fibrinolysis : an international journal in haemostasis and thrombosis.

[30]  William L. Jorgensen,et al.  Aromatic-aromatic interactions: free energy profiles for the benzene dimer in water, chloroform, and liquid benzene , 1990 .

[31]  C. Paradisi,et al.  Investigation of the electronic excited states of C2H3S+ and C2H3O+ by means of collision spectroscopy , 1990 .

[32]  R. Huber,et al.  The refined 1.9 A crystal structure of human alpha‐thrombin: interaction with D‐Phe‐Pro‐Arg chloromethylketone and significance of the Tyr‐Pro‐Pro‐Trp insertion segment. , 1989, The EMBO journal.

[33]  J W Fenton,et al.  Regulation of Thrombin Generation and Functions , 1988, Seminars in thrombosis and hemostasis.

[34]  Michael L. McKee,et al.  Theoretical investigation of the thermal ring opening of bicyclobutane to butadiene. Evidence for a nonsynchronous process , 1988 .

[35]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[36]  J. Singh,et al.  The geometries of interacting arginine‐carboxyls in proteins , 1987, FEBS letters.

[37]  B Honig,et al.  Internal cavities and buried waters in globular proteins. , 1986, Biochemistry.

[38]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[39]  A. Rashin,et al.  Buried surface area, conformational entropy, and protein stability , 1984, Biopolymers.

[40]  David E. Cane,et al.  Biosynthetic origin of the carbon skeleton and oxygen atoms of nargenicin A1 , 1984 .

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

[42]  R. Kikumoto,et al.  Selective inhibition of thrombin by (2R,4R)-4-methyl-1-[N2-[(3-methyl-1,2,3,4-tetrahydro-8-quinolinyl++ +) sulfonyl]-l-arginyl)]-2-piperidinecarboxylic acid. , 1984, Biochemistry.

[43]  R. Doolittle.  Fibrinogen and fibrin. , 1984, Annual review of biochemistry.

[44]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[45]  H. Ekelund,et al.  Fibrinogen to fibrin transformation in umbilical cord blood and purified neonatal fibrinogen. , 1974, Thrombosis research.

[46]  R. Doolittle,et al.  Isolation, characterization, and location of a donor-acceptor unit from cross-linked fibrin. , 1970, Proceedings of the National Academy of Sciences of the United States of America.