Mechanism of proteolysis in matrix metalloproteinase‐2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase‐2 (MMP‐2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all‐atom three‐dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace‐Gln‐Gly∼Ile‐Ala‐Gly‐Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four‐step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and CN bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate‐limiting step in MMP‐2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations. © 2015 Wiley Periodicals, Inc.

[1]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[2]  S. C. Rogers,et al.  QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis , 2003 .

[3]  Frank Neese,et al.  Natural triple excitations in local coupled cluster calculations with pair natural orbitals. , 2013, The Journal of chemical physics.

[4]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[5]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[6]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[7]  M. Coletta,et al.  pH- and temperature-dependence of functional modulation in metalloproteinases. A comparison between neutrophil collagenase and gelatinases A and B. , 2000, Biophysical journal.

[8]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[9]  M. Milla,et al.  Key feature of the catalytic cycle of TNF-α converting enzyme involves communication between distal protein sites and the enzyme catalytic core , 2007, Proceedings of the National Academy of Sciences.

[10]  Dimas Suárez,et al.  Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study. , 2008, The journal of physical chemistry. B.

[11]  Frank Neese,et al.  The ORCA program system , 2012 .

[12]  T L Blundell,et al.  X-ray structure of a hydroxamate inhibitor complex of stromelysin catalytic domain and its comparison with members of the zinc metalloproteinase superfamily. , 1996, Structure.

[13]  J. Cha,et al.  Site-directed mutagenesis of the active site glutamate in human matrilysin: investigation of its role in catalysis. , 1997, Biochemistry.

[14]  A. Eisen,et al.  Human skin fibroblast collagenase. Assessment of activation energy and deuterium isotope effect with collagenous substrates. , 1981, The Journal of biological chemistry.

[15]  Mark A. Ratner,et al.  6-31G * basis set for atoms K through Zn , 1998 .

[16]  Per E M Siegbahn,et al.  Catalytic mechanism of matrix metalloproteinases: two-layered ONIOM study. , 2002, Inorganic chemistry.

[17]  Walter Thiel,et al.  Linear scaling geometry optimisation and transition state search in hybrid delocalised internal coordinates , 2000 .

[18]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[19]  Michele Parrinello,et al.  A hybrid Gaussian and plane wave density functional scheme , 1997 .

[20]  J. O'Connell,et al.  Mutation of the active site glutamic acid of human gelatinase A: effects on latency, catalysis, and the binding of tissue inhibitor of metalloproteinases-1. , 1994, Biochemistry.

[21]  A. Nemukhin,et al.  Exploration of the Zinc Finger Motif in Controlling Activity of Matrix Metalloproteinases , 2014, The journal of physical chemistry. B.

[22]  Michele Parrinello,et al.  Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach , 2005, Comput. Phys. Commun..

[23]  Teter,et al.  Separable dual-space Gaussian pseudopotentials. , 1996, Physical review. B, Condensed matter.

[24]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[25]  Christopher M. Overall,et al.  Validating matrix metalloproteinases as drug targets and anti-targets for cancer therapy , 2006, Nature Reviews Cancer.

[26]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[27]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[28]  H. Birkedal‐Hansen,et al.  Comparative sequence specificities of human 72- and 92-kDa gelatinases (type IV collagenases) and PUMP (matrilysin). , 1993, Biochemistry.

[29]  Alessandro Laio,et al.  An Efficient Real Space Multigrid QM/MM Electrostatic Coupling. , 2005, Journal of chemical theory and computation.

[30]  D. Suárez,et al.  Kinetic and binding effects in peptide substrate selectivity of matrix metalloproteinase‐2: Molecular dynamics and QM/MM calculations , 2010, Proteins.

[31]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[32]  A. Nemukhin,et al.  Computational characterization of ketone-ketal transformations at the active site of matrix metalloproteinases. , 2014, The journal of physical chemistry. B.

[33]  I. Bertini,et al.  Snapshots of the reaction mechanism of matrix metalloproteinases. , 2006, Angewandte Chemie.

[34]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[35]  J. Rivail,et al.  Insights in the Peptide Hydrolysis Mechanism by Thermolysin: A Theoretical QM/MM study , 2000 .

[36]  Alexey A. Sokol,et al.  ChemShell—a modular software package for QM/MM simulations , 2014 .

[37]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[38]  M. Stack,et al.  The effect of pH, temperature, and D2O on the activity of porcine synovial collagenase and gelatinase. , 1990, Archives of biochemistry and biophysics.

[39]  W Smith,et al.  DL_POLY_2.0: a general-purpose parallel molecular dynamics simulation package. , 1996, Journal of molecular graphics.