Chemically modified tetracyclines as inhibitors of MMP-2 matrix metalloproteinase: a molecular and structural study.

The present study focuses on the direct interaction of chemically modified tetracyclines (CMTs) with the active site of the matrix metalloproteinase 2 (MMP-2). Molecular docking, molecular dynamics (MD) simulations, and free energy calculations were accomplished for seven CMT derivatives. New sets of parameters are proposed for structural and catalytic zinc atoms in order to study MMPs and their complexes by means of the AMBER force field. Our computational results show that six CMTs studied bind to the catalytic zinc of the MMP-2 enzyme at the O11-O12 site as proposed experimentally. The exception was the CMT-3 analogue that is found embedding within the active site, enhancing the van der Waals and hydrophobic contacts with the hydrophobic S1' pocket in the MMP-2 enzyme. The binding energy calculated in solution predicts the CMT-3 complexes as the most favorable, followed by the CMT-7 and CMT-8 analogues, respectively, which is in line with experimental findings. This work is the first step toward understanding the mechanism of CMTs as MMP inhibitors at a molecular level.

[1]  P. Kollman,et al.  New-generation amber united-atom force field. , 2006, The journal of physical chemistry. B.

[2]  S. Balaz,et al.  A practical approach to docking of zinc metalloproteinase inhibitors. , 2004, Journal of molecular graphics & modelling.

[3]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[4]  Arieh Warshel,et al.  Free energy relationships in metalloenzyme-catalyzed reactions. Calculations of the effects of metal ion substitutions in staphylococcal nuclease , 1990 .

[5]  David L Mobley,et al.  Nonlinear scaling schemes for Lennard-Jones interactions in free energy calculations. , 2007, The Journal of chemical physics.

[6]  Wenfang Xu,et al.  Progress in the development of matrix metalloproteinase inhibitors. , 2008, Current medicinal chemistry.

[7]  Sophie Perrier,et al.  Crystal structures of novel non-peptidic, non-zinc chelating inhibitors bound to MMP-12. , 2004, Journal of molecular biology.

[8]  M. Mita,et al.  A phase II and pharmacological study of the matrix metalloproteinase inhibitor (MMPI) COL-3 in patients with advanced soft tissue sarcomas , 2007, Investigational New Drugs.

[9]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[10]  Thomas Simonson,et al.  Free energy of particle insertion , 1993 .

[11]  M Karplus,et al.  Zinc binding in proteins and solution: A simple but accurate nonbonded representation , 1995, Proteins.

[12]  L. Curtiss,et al.  Gaussian-3 (G3) theory for molecules containing first and second-row atoms , 1998 .

[13]  A. Simas,et al.  Importance of tautomers in the chemical behavior of tetracyclinesdagger. , 1999, Journal of pharmaceutical sciences.

[14]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[15]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[16]  L. Pantanowitz,et al.  Antimicrobial and non-antimicrobial tetracyclines in human cancer trials. , 2011, Pharmacological research.

[17]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[18]  Sérgio F. Sousa,et al.  Effective tailor-made force field parameterization of the several Zn coordination environments in the puzzling FTase enzyme: opening the door to the full understanding of its elusive catalytic mechanism , 2006 .

[19]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[20]  A. Sparreboom,et al.  Chemically modified tetracyclines as inhibitors of matrix metalloproteinases. , 2004, Drug resistance updates : reviews and commentaries in antimicrobial and anticancer chemotherapy.

[21]  Maria J. Ramos,et al.  vsLab—An implementation for virtual high‐throughput screening using AutoDock and VMD , 2011 .

[22]  Arnaud Blondel,et al.  Ensemble variance in free energy calculations by thermodynamic integration: Theory, optimal “Alchemical” path, and practical solutions , 2004, J. Comput. Chem..

[23]  A. Henney,et al.  Crystal structure of human MMP9 in complex with a reverse hydroxamate inhibitor. , 2002, Journal of molecular biology.

[24]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[25]  M. Zerner,et al.  Conformational analysis of the anhydrotetracycline molecule: a toxic decomposition product of tetracycline. , 1998, Journal of pharmaceutical sciences.

[26]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[27]  W. B. Almeida,et al.  The conformational and tautomeric equilibrium of 5a,6-anhydrotetracycline in aqueous solution at pH 7 , 2003 .

[28]  Armando Rossello,et al.  Amber force field implementation, molecular modelling study, synthesis and MMP-1/MMP-2 inhibition profile of (R)- and (S)-N-hydroxy-2-(N-isopropoxybiphenyl-4-ylsulfonamido)-3-methylbutanamides. , 2006, Bioorganic & medicinal chemistry.

[29]  F. Villarreal,et al.  Tetracycline compounds with non-antimicrobial organ protective properties: Possible mechanisms of action , 2010, Pharmacological Research.

[30]  M. Hendrix,et al.  Chemically modified tetracyclines inhibit human melanoma cell invasion and metastasis , 1998, Clinical & Experimental Metastasis.

[31]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[32]  M. Zimmer,et al.  BIOINORGANIC MOLECULAR MECHANICS , 1995 .

[33]  Thomas Simonson,et al.  Molecular mechanics models for tetracycline analogs , 2009, J. Comput. Chem..

[34]  L. Golub,et al.  Excessive matrix metalloproteinase activity in diabetes: inhibition by tetracycline analogues with zinc reactivity. , 2001, Current medicinal chemistry.

[35]  R. Schmid,et al.  Modeling domino effects in enzymes: molecular basis of the substrate specificity of the bacterial metallo-beta-lactamases IMP-1 and IMP-6. , 2003, Biochemistry.

[36]  B. Lokeshwar,et al.  Cytotoxic activity and inhibition of tumor cell invasion by derivatives of a chemically modified tetracycline CMT-3 (COL-3). , 2001, Current medicinal chemistry.

[37]  D. Suárez,et al.  Molecular dynamics simulations of matrix metalloproteinase 2: role of the structural metal ions. , 2007, Biochemistry.

[38]  W. B. Almeida,et al.  Interaction of chemically modified tetracyclines with catalytic Zn(II) ion in matrix metalloproteinase: evidence for metal coordination sites , 2011 .

[39]  Tingjun Hou,et al.  Binding affinity of hydroxamate inhibitors of matrix metalloproteinase-2 , 2004, Journal of molecular modeling.

[40]  P A Fernandes,et al.  Substrate recognition in HIV-1 protease: a computational study. , 2010, The journal of physical chemistry. B.

[41]  Junmei Wang,et al.  Junmei Wang, Romain M. Wolf, James W. Caldwell, Peter A. Kollman, and David A. Case, "Development and testing of a general amber force field"Journal of Computational Chemistry(2004) 25(9) 1157–1174 , 2005, J. Comput. Chem..

[42]  A. W. Addison,et al.  Synthesis, structure, and spectroscopic properties of copper(II) compounds containing nitrogen–sulphur donor ligands; the crystal and molecular structure of aqua[1,7-bis(N-methylbenzimidazol-2′-yl)-2,6-dithiaheptane]copper(II) perchlorate , 1984 .

[43]  K. Appelt,et al.  Matrix Metalloproteinase Inhibitors for Cancer Therapy , 2000, Cancer Drug Discovery and Development.

[44]  P A Kollman,et al.  Calculation and prediction of binding free energies for the matrix metalloproteinases. , 2000, Journal of medicinal chemistry.

[45]  Robert Huber,et al.  Crystal structures of MMP-9 complexes with five inhibitors: contribution of the flexible Arg424 side-chain to selectivity. , 2007, Journal of molecular biology.

[46]  H. Stoll,et al.  Energy-adjustedab initio pseudopotentials for the second and third row transition elements , 1990 .

[47]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[48]  Tom L. Blundell,et al.  X-RAY STRUCTURE OF GELATINASE A CATALYTIC DOMAIN COMPLEXED WITH A HYDROXAMATE INHIBITOR , 1999 .

[49]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[50]  Xiaojie Xu,et al.  Predictions of Binding of a Diverse Set of Ligands to Gelatinase-A by a Combination of Molecular Dynamics and Continuum Solvent Models , 2002 .

[51]  Samuel H. Wilson,et al.  Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects. , 2007, Journal of molecular biology.

[52]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

[53]  A. Warshel,et al.  Computer simulation of the initial proton transfer step in human carbonic anhydrase I. , 1992, Journal of molecular biology.

[54]  T. Steinbrecher,et al.  Towards accurate free energy calculations in ligand protein-binding studies. , 2010, Current medicinal chemistry.

[55]  Mark A. Ratner,et al.  6‐31G* basis set for third‐row atoms , 2001, J. Comput. Chem..

[56]  Xiaojie Xu,et al.  Binding Affinities for a Series of Selective Inhibitors of Gelatinase-A Using Molecular Dynamics with a Linear Interaction Energy Approach , 2001 .

[57]  Arieh Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

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

[59]  A. Warshel,et al.  On the origin of the catalytic power of carboxypeptidase A and other metalloenzymes , 2009, Proteins.