Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study

Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1

Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.

Catalytic Mechanism of Collagen Hydrolysis by Zinc(II)-Dependent Matrix Metalloproteinase-1

Human matrix metalloproteinase-1 (MMP-1) is a zinc(II)-dependent enzyme that catalyzes collagenolysis. Despite the availability of extensive experimental data, the mechanism of MMP-1-catalyzed collagenolysis remains poorly understood due to the lack of experimental structure of a catalytically productive enzyme-substrate complex of MMP-1. In this study, we apply molecular dynamics and combined quantum mechanics/molecular mechanics to reveal the reaction mechanism of MMP-1 based on a computationally modeled structure of the catalytically competent complex of MMP-1 that contains a large triple-helical peptide substrate. Our proposed mechanism involves the participation of an auxiliary (second) water molecule (wat2) in addition to the zinc(II)-coordinated water (wat1). The reaction initiates through a proton transfer to Glu219, followed by a nucleophilic attack by a zinc(II)-coordinated hydroxide anion nucleophile at the carbonyl carbon of the scissile bond, leading to the formation of a tetrahedral intermediate (IM2). The process continues with a hydrogen-bond rearrangement to facilitate proton transfer from wat2 to the amide nitrogen of the scissile bond and, finally, C-N bond cleavage. The calculations indicate that the rate-determining step is the water-mediated nucleophilic attack with an activation energy barrier of 22.3 kcal/mol. Furthermore, the calculations show that the hydrogen-bond rearrangement/proton-transfer step can proceed in a consecutive or concerted manner, depending on the conformation of the tetrahedral intermediate, with the consecutive mechanism being energetically preferable. Overall, the study reveals the crucial role of a second water molecule and the dynamics for effective MMP-1-catalyzed collagenolysis.

Effects of the Nature of the Metal Ion, Protein and Substrate on the Catalytic Center in Matrix Metalloproteinase-1: Insights from Multilevel MD, QM/MM and QM Studies

Matrix metalloproteinase-1 (MMP-1) is a Zn(II) dependent endopeptidase involved in the degradation of collagen, the most abundant structural protein in the extracellular matrix of connective tissues and the human body. Herein we performed a multilevel computational analysis including molecular dynamics (MD), combined quantum mechanics/molecular mechanics (QM/MM), and quantum mechanics (QM) calculations to characterize the structure and geometry of the catalytic Zn(II) within the MMP-1 protein environment in comparison to crystallographic and spectroscopic data. The substrate's removal fine-tuned impact on the conformational dynamics and geometry of the catalytic Zn(II) center was also explored. Finally, the study examined the effect of substituting catalytic Zn(II) by Co(II) on the overall structure and dynamics of the MMP-1 THP complex and specifically on the geometry of the catalytic metal center. Overall our QM/MM and QM studies were in good agreement with the MM description of the Zn(II) centers in the MD simulations.

Hydrolysis Mechanism of the Linkers by Matrix Metalloproteinase-9 Using QM/MM Calculations

Activatable cell-penetrating peptides (ACPPs) are known to be able to decrease the cytotoxicity of cell-penetrating peptide (CPP)-based drug delivery systems. Furthermore, they can improve the targeting of CPPs when specifically recognized and hydrolyzed by characteristic proteases. A comprehensive and profound understanding of the recognition and hydrolysis process will provide a better design of the ACPP-based drug delivery system. Previous studies have clearly described how ACPPs are recognized and bound by MMPs. However, the hydrolysis mechanism of ACPPs is still unsolved. This work focuses on a proteinase-sensitive cleavable linker of ACPPs (PLGLAG), the key structure for recognition and hydrolysis, trying to determine the mechanism by which MMP-9 hydrolyzes its substrate PLGLAG. The quantum mechanics/molecular mechanics (QM/MM) calculations herein show that MMP-9 proteolysis is a water-mediated four-step reaction. More specifically, it consists of (i) nucleophilic attack, (ii) hydrogen-bond rearrangement, (iii) proton transfer, and finally (iv) amide bond rupture. Considering the reversibility of multistep reaction, the second step (i.e., hydrogen-bond rearrangement) has the highest barrier and is the rate-limiting step in the hydrolysis of PLGLAG. The possible design and improvement of the key P1 and P1' sites are also explored through mutations. The present results indicate that, while the mutations affect the reaction energy barriers and the rate-limiting steps, all mutants considered could be hydrolyzed by MMP-9. To provide further insights, the hydrolysis mechanism of MMP-2, which has a similar hydrolysis process to that of MMP-9 but with different reaction barriers, is also studied and compared. As a result, this work provides detailed insights into the hydrolysis mechanism of ACPPs by MMP-9 and, thus, also possible insights for the development of new strategies for ACPP-based delivery systems.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Peptide Hydrolysis by Metal (Oxa)cyclen Complexes: Revisiting the Mechanism and Assessing Ligand Effects

The mechanism responsible for peptide bond hydrolysis by Co(III) and Cu(II) complexes with (oxa)cyclen ligands has been revisited by means of computational tools. We propose that the mechanism starts by substrate coordination and an outer-sphere attack on the amide C atom of a solvent water molecule assisted by the metal hydroxo moiety as a general base, which occurs through six-membered ring transition states. This new mechanism represents a more likely scenario than the previously proposed mechanisms that involved an inner-sphere nucleophilic attack through more strained four-membered rings transition states. The corresponding computed overall free-energy barrier of 25.2 kcal mol^-1 for hydrolysis of the peptide bond in Phe-Ala by a cobalt(III) oxacyclen catalyst (1) is consistent with the experimental values obtained from rate constants. Also, we assessed the influence of the nature of the ligand throughout a systematic replacement of N by O atoms in the (oxa)cyclen ligand. Increasing the number of coordinating O atoms accelerates the reaction by increasing the Lewis acidity of the metal ion. On the other hand, the higher reactivity observed for the copper(II) oxacyclen catalyst with respect to the analogous Co(III) complex can be attributed to the larger Brönsted basicity of the copper(II) hydroxo ligand. Ultimately, the detailed understanding of the ligand and metal nature effects allowed us to identify the double role of the metal hydroxo complexes as Lewis acids and Brönsted bases and to rationalize the observed reactivity trends.

Conformation and Domain Movement Analysis of Human Matrix Metalloproteinase-2: Role of Associated Zn2+ and Ca2+ Ions

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn²⁺ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn²⁺ ion 1, Ca²⁺ ion 1, and Ca²⁺ ion 3 with significant conformational stability at the binding sites of Zn²⁺ ion 1 and Ca²⁺ ion 1. Ca²⁺ ion 2 diffuses freely away from its crystallographically defined binding site. Zn²⁺ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca²⁺ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca²⁺ ion 3′s coordination sphere.

Solvent water interactions within the active site of the membrane type I matrix metalloproteinase

Matrix metalloproteinases (MMP) are an important family of proteases which catalyze the degradation of extracellular matrix components. While the mechanism of peptide cleavage is well established, the process of enzyme regeneration, which represents the rate limiting step of the catalytic cycle, remains unresolved. This step involves the loss of the newly formed N-terminus (amine) and C-terminus (carboxylate) protein fragments from the site of catalysis coupled with the inclusion of one or more solvent waters. Here we report a novel crystal structure of membrane type I MMP (MT1-MMP or MMP-14), which includes a small peptide bound at the catalytic Zn site via its C-terminus. This structure models the initial product state formed immediately after peptide cleavage but before the final proton transfer to the bound amine; the amine is not present in our system and as such proton transfer cannot occur. Modeling of the protein, including earlier structural data of Bertini and coworkers [I. Bertini, et al., Angew. Chem., Int. Ed., 2006, 45, 7952-7955], suggests that the C-terminus of the peptide is positioned to form an H-bond network to the amine site, which is mediated by a single oxygen of the functionally important Glu240 residue, facilitating efficient proton transfer. Additional quantum chemical calculations complemented with magneto-optical and magnetic resonance spectroscopies clarify the role of two additional, non-catalytic first coordination sphere waters identified in the crystal structure. One of these auxiliary waters acts to stabilize key intermediates of the reaction, while the second is proposed to facilitate C-fragment release, triggered by protonation of the amine. Together these results complete the enzymatic cycle of MMPs and provide new design criteria for inhibitors with improved efficacy.

[(H2O)Zn(Imidazole)n]2+: the vital roles of coordination number and geometry in Zn–OH2 acidity and catalytic hydrolysis

The Zn(ii)–(Imidazole(ate))_n coordination motif occurs in numerous biochemical systems, including carbonic anhydrase and the matrix metalloproteinases (MMPs).

Molecular Dynamics Studies of Matrix Metalloproteases

Matrix metalloproteases are multidomain enzymes with a remarkable proteolytic activity located in the extracellular environment. Their catalytic activity and structural properties have been intensively studied during the last few decades using both experimental and theoretical approaches, but many open questions still remain. Extensive molecular dynamics simulations enable the sampling of the configurational space of a molecular system, thus contributing to the characterization of the structure, dynamics, and ligand binding properties of a particular MMP. Based on previous computational experience, we provide in this chapter technical and methodological guidelines that may be useful to and stimulate other researchers to perform molecular dynamics simulations to help address unresolved questions concerning the molecular mode of action of MMPs.

Effects of Mutations on Structure-Function Relationships of Matrix Metalloproteinase-1

Matrix metalloproteinase-1 (MMP-1) is one of the most widely studied enzymes involved in collagen degradation. Mutations of specific residues in the MMP-1 hemopexin-like (HPX) domain have been shown to modulate activity of the MMP-1 catalytic (CAT) domain. In order to reveal the structural and conformational effects of such mutations, a molecular dynamics (MD) study was performed of in silico mutated residues in the X-ray crystallographic structure of MMP-1 complexed with a collagen-model triple-helical peptide (THP). The results indicate an important role of the mutated residues in MMP-1 interactions with the THP and communication between the CAT and the HPX domains. Each mutation has a distinct impact on the correlated motions in the MMP-1•THP. An increased collagenase activity corresponded to the appearance of a unique anti-correlated motion and decreased correlated motions, while decreased collagenase activity corresponded both to increased and decreased anti-correlated motions.

Pressure and Temperature Effects on the Activity and Structure of the Catalytic Domain of Human MT1-MMP

Membrane type 1-matrix metalloproteinase (MT1-MMP or MMP-14) is a zinc-transmembrane metalloprotease involved in the degradation of extracellular matrix and tumor invasion. While changes in solvation of MT1-MMP have been recently studied, little is known about the structural and energetic changes associated with MT1-MMP while interacting with substrates. Steady-state kinetic and thermodynamic data (including activation energies and activation volumes) were measured over a wide range of temperatures and pressures by means of a stopped-flow fluorescence technique. Complementary temperature- and pressure-dependent Fourier-transform infrared measurements provided corresponding structural information of the protein. MT1-MMP is stable and active over a wide range of temperatures (10-55 °C). A small conformational change was detected at 37 °C, which is responsible for the change in activity observed at the same temperature. Pressure decreases the enzymatic activity until complete inactivation occurs at 2 kbar. The inactivation is associated with changes in the rate-limiting step of the reaction caused by additional hydration of the active site upon compression and/or minor conformational changes in the active site region. Based on these data, an energy and volume diagram could be established for the various steps of the enzymatic reaction.

Methodological aspects of QM/MM calculations: A case study on matrix metalloproteinase-2

We address methodological issues in quantum mechanics/molecular mechanics (QM/MM) calculations on a zinc-dependent enzyme. We focus on the first stage of peptide bond cleavage by matrix metalloproteinase-2 (MMP-2), that is, the nucleophilic attack of the zinc-coordinating water molecule on the carbonyl carbon atom of the scissile fragment of the substrate. This step is accompanied by significant charge redistribution around the zinc cation, bond cleavage, and bond formation. We vary the size and initial geometry of the model system as well as the computational protocol to demonstrate the influence of these choices on the results obtained. We present QM/MM potential energy profiles for a set of snapshots randomly selected from QM/MM-based molecular dynamics simulations and analyze the differences in the computed profiles in structural terms. Since the substrate in MMP-2 is located on the protein surface, we investigate the influence of the thickness of the water layer around the enzyme on the QM/MM energy profile. Thin water layers (0-2 Å) give unrealistic results because of structural reorganizations in the active-site region at the protein surface. A 12 Å water layer appears to be sufficient to capture the effect of the solvent; the corresponding QM/MM energy profile is very close to that obtained from QM/MM/SMBP calculations using the solvent macromolecular boundary potential (SMBP). We apply the optimized computational protocol to explain the origin of the different catalytic activity of the Glu116Asp mutant: the energy barrier for the first step is higher, which is rationalized on structural grounds. © 2016 Wiley Periodicals, Inc.

Importance of the Linker Region in Matrix Metalloproteinase-1 Domain Interactions

Collagenolysis is catalyzed by enzymes from the matrix metalloproteinase (MMP) family, where one of the most studied is MMP-1. The X-ray crystallographic structure of MMP-1 complexed with a collagen-model triple-helical peptide (THP) provided important atomistic information, but few details on the effects of the conformational flexibility on catalysis. In addition, the role of the linker region between the catalytic (CAT) and hemopexin-like (HPX) domains was not defined. In order to reveal the dynamics and correlations of MMP-1 comprehensive atomistic molecular dynamics simulations of an MMP-1•THP complex was performed. To examine the role of the linker region for MMP-1 function simulations with linker regions from MT1-MMP/MMP-14 and MMP-13 replacing the MMP-1 linker region were performed. The MD studies were in good agreement with the experimental observation that in the MMP-1•THP X-ray crystallographic structure MMP-1 is in a "closed" conformation. MD revealed that the interactions of the THP with the both the CAT and HPX domains of MMP-1 are dynamic in nature, and the linker region of MMP-1 influences the interactions and dynamics of both the CAT and HPX domains and collagen binding to MMP-1.

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.

Exploration of the zinc finger motif in controlling activity of matrix metalloproteinases

Discovering ways to control the activity of matrix metalloproteinases (MMPs), zinc-dependent enzymes capable of degrading extracellular matrix proteins, is an important field of cancer research. We report here a novel strategy for assembling MMP inhibitors on the basis of oligopeptide ligands by exploring the pattern known as the zinc finger motif. Advanced molecular modeling tools were used to characterize the structural binding motifs of experimentally tested MMP inhibitors, as well as those of newly proposed peptidomimetics, in their zinc-containing active sites. The results of simulations based on the quantum mechanics/molecular mechanics (QM/MM) approach and Car-Parrinello molecular dynamics with QM/MM potentials demonstrate that, upon binding of Regasepin1, a known MMP-9 inhibitor, the Zn(2+)(His3) structural element is rearranged to the Zn(2+)(Cys2His2) zinc finger motif, in which two Cys residues are borrowed from the ligand. Following consideration of the crystal structure of MMP-2 with its inhibitor, the oligopeptide APP-IP, we proposed a new peptidomimetic with two replacements in the substrate, Tyr3Cys and Asp6Cys. Simulations show that this peptide variant blocks an enzyme active site by the Zn(2+)(Cys2His2) zinc finger construct. Similarly, a natural substrate of MMP-2, Ace-Gln-Gly ∼ Ile-Ala-Gly-Nme, can be converted to an inhibiting compound by two replacements, Ile by Cys and Gly by the d isomer of Cys, favoring formation of the zinc finger motif.

Development and in vitro assessment of enzymatically-responsive poly(ethylene glycol) hydrogels for the delivery of therapeutic peptides

Despite the recent expansion of peptide drugs, delivery remains a challenge due to poor localization and rapid clearance. Therefore, a hydrogel-based platform technology was developed to control and sustain peptide drug release via matrix metalloproteinase (MMP) activity. Specifically, hydrogels were composed of poly(ethylene glycol) and peptide drugs flanked by MMP substrates and terminal cysteine residues as crosslinkers. First, peptide drug bioactivity was investigated in expected released forms (e.g., with MMP substrate residues) in vitro prior to incorporation into hydrogels. Three peptides (Qk (from Vascular Endothelial Growth Factor), SPARC113, and SPARC118 (from Secreted Protein Acidic and Rich in Cysteine)) retained bioactivity and were used as hydrogel crosslinkers in full MMP degradable forms. Upon treatment with MMP2, hydrogels containing Qk, SPARC113, and SPARC118 degraded in 6.7, 6, and 1 days, and released 5, 8, and, 19% of peptide, respectively. Further investigation revealed peptide drug size controlled hydrogel swelling and degradation rate, while hydrophobicity impacted peptide release. Additionally, Qk, SPARC113, and SPARC118 releasing hydrogels increased endothelial cell tube formation 3.1, 1.7, and 2.8-fold, respectively. While pro-angiogenic peptides were the focus of this study, the design parameters detailed allow for adaptation of hydrogels to control peptide release for a variety of therapeutic applications.

Computational characterization of ketone-ketal transformations at the active site of matrix metalloproteinases

We modeled the first steps of hydrolysis reactions of a natural oligopeptide substrate of matrix metalloproteinase MMP-2 as well as of a substrate analogue. In the latter, the scissile amide group is substituted by a ketomethylene group which can be transformed to the ketal group upon binding of this compound to the enzyme active site. According to our quantum mechanical-molecular mechanical (QM/MM) calculations, the reaction of the ketone-ketal transformation proceeds with a low energy barrier (3.4 kcal/mol) and a high equilibrium constant (10(4)). The reaction product with the ketal group formed directly at the active site of the enzyme works as an inhibitor that chelates the zinc ion. On the other hand, the oligopeptide mimetic retains molecular groups responsible for binding of this compound to the enzyme active site. This example illustrates a strategy to design MMP inhibitors in situ by using data on binding specificity of substrates to a particular type of MMP and details of the reaction mechanism.

Unraveling the molecular structure of the catalytic domain of matrix metalloproteinase-2 in complex with a triple-helical peptide by means of molecular dynamics simulations

Herein, we present the results of a computational study that employed various simulation methodologies to build and validate a series of molecular models of a synthetic triple-helical peptide (fTHP-5) both in its native state and in a prereactive complex with the catalytic domain of the MMP-2 enzyme. First, the structure and dynamical properties of the fTHP-5 substrate are investigated by means of molecular dynamics (MD) simulations. Then, the propensity of each of the three peptide chains in fTHP-5 to be distorted around the scissile peptide bond is assessed by carrying out potential of mean force calculations. Subsequently, the distorted geometries of fTHP-5 are docked within the MMP-2 active site following a semirigid protocol, and the most stable docked structures are fully relaxed and characterized by extensive MD simulations in explicit solvent. Following a similar approach, we also investigate a hypothetical ternary complex formed between two MMP-2 catalytic units and a single fTHP-5 molecule. Overall, our models for the MMP-2/fTHP-5 complexes unveil the extent to which the triple helix is distorted to allow the accommodation of an individual peptide chain within the MMP active site.

Comparison of Three Chain-of-States Methods: Nudged Elastic Band and Replica Path with Restraints or Constraints

Chain-of-state methods are becoming important tools in studying the chemical reaction mechanisms, especially for biomacromolecules. In this article, three chain-of-state methods, nudged elastic band (NEB) method and the replica path method with restraints or constraints, were tested and compared using three model systems with various sizes and at different levels of theory: alanine dipeptide isomerization, β-alanine intramolecular condensation, and the matrix metalloproteinase 2 inhibition mechanism. The levels of theory used to describe the three model systems include molecular mechanics (MM), quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM). All three methods could correctly determine a reaction path with reasonable estimation of reaction barriers in most cases. The RMSD measurement with additional weighting schemes provides practically infinite choices of reaction coordinates to describe the reaction progress. These findings demonstrate that the chain-of-state methods are powerful tools when being used carefully to generate a plausible reaction mechanism with full pathway for complex systems at an affordable computational cost.

Kinetic and binding effects in peptide substrate selectivity of matrix metalloproteinase-2: Molecular dynamics and QM/MM calculations

Herein, we examine computationally the binding and hydrolysis reaction of the MMP-2 enzyme with two peptide substrates selected by the enzyme from a phage peptide library. Molecular dynamics simulations of the Michaelis complexes (25 ns) allow us to characterize the main enzyme/substrate contacts. Subsequently MM-PBSA calculations using independent trajectories for the complexes and the free substrates provide relative binding energies in good agreement with the experimental K(M) results. Computational alanine scanning analyses of the enzyme/substrate interaction energies confirm the relevance of the P(3), P(2), and P(1)' side chains for ligand binding. Finally, the hydrolysis of both peptides taking place at the MMP-2 active site is explored by means of hybrid quantum mechanical/molecular mechanics calculations. The computed reaction mechanisms result in rate-determining energy barriers being in consonance with the experimental k(cat) values. Overall, the computational protocol seems to capture the subtle differences in binding and catalysis experimentally observed for the two peptide substrates. Some implications of our results for the future design of novel and more specific MMP-2 inhibitors are also discussed.