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

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.

Cysteine Radical and Glutamate Collaboratively Enable C-H Bond Activation and C-N Bond Cleavage in a Glycyl Radical Enzyme HplG

Hydroxyprolines are abundant in nature and widely utilized by many living organisms. Isomerization of trans-4-hydroxy-d-proline (t4D-HP) to generate 2-amino-4-ketopentanoate has been found to need a glycyl radical enzyme HplG, which catalyzes the cleavage of the C-N bond, while dehydration of trans-4-hydroxy-l-proline involves a homologous enzyme of HplG. Herein, molecular dynamics simulations and quantum mechanics/molecular mechanics (QM/MM) calculations are employed to understand the reaction mechanism of HplG. Two possible reaction pathways of HplG have been explored to decipher the origin of its chemoselectivity. The QM/MM calculations reveal that the isomerization proceeds via an initial hydrogen shift from the Cγ site of t4D-HP to a catalytic cysteine radical, followed by cleavage of the Cδ-N bond in t4D-HP to form a radical intermediate that captures a hydrogen atom from the cysteine. Activation of the Cδ-H bond in t4D-HP to bring about dehydration of t4D-HP possesses an extremely high energy barrier, thus rendering the dehydration pathway implausible in HplG. On the basis of the current calculations, conserved residue Glu429 plays a pivotal role in the isomerization pathway: the hydrogen bonding between it and t4D-HP weakens the hydroxyalkyl Cγ-Hγ bond, and it acts as a proton acceptor to trigger the cleavage of the C-N bond in t4D-HP. Our current QM/MM calculations rationalize the origin of the experimentally observed chemoselectivity of HplG and propose an H-bond-assisted bond activation strategy in radical-containing enzymes. These findings have general implications on radical-mediated enzymatic catalysis and expand our understanding of how nature wisely and selectively activates the C-H bond to modulate catalytic selectivity.

Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins

The photochemically active sites of the proteins sfGFP66azF and Venus66azF, members of the green fluorescent protein (GFP) family, contain a non-canonical amino acid residue p-azidophenylalanine (azF) instead of Tyr66. The light-induced decomposition of azF at these sites leads to the formation of reactive arylnitrene (nF) intermediates followed by the formation of phenylamine-containing chromophores. We report the first study of the reaction mechanism of the reduction of the arylnitrene intermediates in sfGFP66nF and Venus66nF using molecular modeling methods. The Gibbs energy profiles for the elementary steps of the chemical reaction in sfGFP66nF are computed using molecular dynamics simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. Structures and energies along the reaction pathway in Venus66nF are evaluated using a QM/MM approach. According to the results of the simulations, arylnitrene reduction is coupled with oxidation of the histidine side chain on the His148 residue located near the chromophore.

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.

Molecular Dynamics Simulations of Matrix Metalloproteinase 13 and the Analysis of the Specificity Loop and the S1'-Site

The specificity loop of Matrix Metalloproteinases (MMPs) is known to regulate recognition of their substrates, and the S1'-site surrounded by the loop is a unique place to address the selectivity of ligands toward each MMP. Molecular dynamics (MD) simulations of apo-MMP-13 and its complex forms with various ligands were conducted to identify the role of the specificity loop for the ligand binding to MMP-13. The MD simulations showed the dual role of T247 as a hydrogen bond donor to the ligand, as well as a contributor to the formation of the van der Waal surface area, with T245 and K249 on the S1'-site. The hydrophobic surface area mediated by T247 blocks the access of water molecules to the S1'-site of MMP-13 and stabilizes the ligand in the site. The F252 residue is flexible in order to search for the optimum location in the S1'-site of the apo-MMP-13, but once a ligand binds to the S1'-site, it can form offset π-π or edge-to-π stacking interactions with the ligand. Lastly, H222 and Y244 provide the offset π-π and π-CH(Cβ) interactions on each side of the phenyl ring of the ligand, and this sandwiched interaction could be critical for the ligand binding to MMP-13.

Unraveling the Reaction Mechanism of Russell's Viper Venom Factor X Activator: A Paradigm for the Reactivity of Zinc Metalloproteinases?

Snake venom metalloproteinases (SVMPs) are important drug targets against snakebite envenoming, the neglected tropical disease with the highest mortality worldwide. Here, we focus on Russell's viper (Daboia russelii), one of the "big four" snakes of the Indian subcontinent that, together, are responsible for ca. 50,000 fatalities annually. The "Russell's viper venom factor X activator" (RVV-X), a highly toxic metalloproteinase, activates the blood coagulation factor X (FX), leading to the prey's abnormal blood clotting and death. Given its tremendous public health impact, the WHO recognized an urgent need to develop efficient, heat-stable, and affordable-for-all small-molecule inhibitors, for which a deep understanding of the mechanisms of action of snake's principal toxins is fundamental. In this study, we determine the catalytic mechanism of RVV-X by using a density functional theory/molecular mechanics (DFT:MM) methodology to calculate its free energy profile. The results showed that the catalytic process takes place via two steps. The first step involves a nucleophilic attack by an in situ generated hydroxide ion on the substrate carbonyl, yielding an activation barrier of 17.7 kcal·mol^-1, while the second step corresponds to protonation of the peptide nitrogen and peptide bond cleavage with an energy barrier of 23.1 kcal·mol^-1. Our study shows a unique role played by Zn²⁺ in catalysis by lowering the pK_a of the Zn²⁺-bound water molecule, enough to permit the swift formation of the hydroxide nucleophile through barrierless deprotonation by the formally much less basic Glu140. Without the Zn²⁺ cofactor, this step would be rate-limiting.

Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1

Metalloproteinase-1 (MMP-1) catalyzed collagen degradation is essential for a wide variety of normal physiological processes, while at the same time contributing to several diseases in humans. Therefore, a comprehensive understanding of this process is of great importance. Although crystallographic and spectroscopic studies provided fundamental information about the structure and function of MMP-1, the precise mechanism of collagen degradation especially considering the complex and flexible structure of the substrate, remains poorly understood. In addition, how the protein environment dynamically reorganizes at the atomic scale into a catalytically active state capable of collagen hydrolysis remains unknown. In this study, we applied experimentally-guided multiscale molecular modeling methods including classical molecular dynamics (MD), well-tempered (WT) classical metadynamics (MetD), combined quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM MetD simulations to explore and characterize the early catalytic events of MMP-1 collagenolysis. Importantly the study provided a complete atomic and dynamic description of the transition from the open to the closed form of the MMP-1•THP complex. Notably, the formation of catalytically active Michaelis complex competent for collagen cleavage was characterized. The study identified the changes in the coordination state of the catalytic zinc(II) associated with the conformational transformation and the formation of catalytically productive ES complex. Our results confirm the essential role of the MMP-1 catalytic domain's α-helices (hA, hB and hC) and the linker region in the transition to the catalytically competent ES complex. Overall, the results provide unique mechanistic insight into the conformational transformations and associated changes in the coordination state of the catalytic zinc(II) that would be important for the design of effective MMP-1 inhibitors.

Demystifying DPP III Catalyzed Peptide Hydrolysis—Computational Study of the Complete Catalytic Cycle of Human DPP III Catalyzed Tynorphin Hydrolysis

Dipeptidyl peptides III (DPP III) is a dual-domain zinc exopeptidase that hydrolyzes peptides of varying sequence and size. Despite attempts to elucidate and narrow down the broad substrate-specificity of DPP III, there is no explanation as to why some of them, such as tynorphin (VVYPW), the truncated form of the endogenous heptapeptide spinorphin, are the slow-reacting substrates of DPP III compared to others, such as Leu-enkephalin. Using quantum molecular mechanics calculations followed by various molecular dynamics techniques, we describe for the first time the entire catalytic cycle of human DPP III, providing theoretical insight into the inhibitory mechanism of tynorphin. The chemical step of peptide bond hydrolysis and the substrate binding to the active site of the enzyme and release of the product were described for DPP III in complex with tynorphin and Leu-enkephalin and their products. We found that tynorphin is cleaved by the same reaction mechanism determined for Leu-enkephalin. More importantly, we showed that the product stabilization and regeneration of the enzyme, but not the nucleophilic attack of the catalytic water molecule and inversion at the nitrogen atom of the cleavable peptide bond, correspond to the rate-determining steps of the overall catalytic cycle of the enzyme.

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.

Boronic Acids as Prospective Inhibitors of Metallo-β-Lactamases: Efficient Chemical Reaction in the Enzymatic Active Site Revealed by Molecular Modeling

Boronic acids are prospective compounds in inhibition of metallo-β-lactamases as they form covalent adducts with the catalytic hydroxide anion in the enzymatic active site upon binding. We compare this chemical reaction in the active site of the New Delhi metallo-β-lactamase (NDM-1) with the hydrolysis of the antibacterial drug imipenem. The nucleophilic attack occurs with the energy barrier of 14 kcal/mol for imipenem and simultaneously upon binding a boronic acid inhibitor. A boron atom of an inhibitor exhibits stronger electrophilic properties than the carbonyl carbon atom of imipenem in a solution that is quantified by atomic Fukui indices. Upon forming the prereaction complex between NDM-1 and inhibitor, the lone electron pair of the nucleophile interacts with the vacant p-orbital of boron that facilitates the chemical reaction. We analyze a set of boronic acid compounds with the benzo[b]thiophene core complexed with the NDM-1 and propose quantitative structure-sroperty relationship (QSPR) equations that can predict IC50 values from the calculated descriptors of electron density. These relations are applied to classify other boronic acids with the same core found in the database of chemical compounds, PubChem, and proposed ourselves. We demonstrate that the IC50 values for all considered benzo[b]thiophene-containing boronic acid inhibitors are 30-70 μM.

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.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Proof of concept for poor inhibitor binding and efficient formation of covalent adducts of KRASG12C and ARS compounds

The use of selective covalent inhibitors with low binding affinity and high reactivity with the target enzyme is a promising way to solve a long-standing problem of the "undruggable" RAS-like proteins. Specifically, compounds of the ARS family that prevent the activation of the GDP-bound G12C mutant of Kirsten RAS (KRAS) are in the focus of recent experimental research. We report the first computational characterization of the entire reaction mechanism of the covalent binding of ARS-853 to the KRASG12C·GDP complex. The application of molecular dynamics, molecular docking and quantum mechanics/molecular mechanics approaches allowed us to model the inhibitor binding to the protein and the chemical reaction of ARS-853 with Cys12 in the enzyme binding site. We estimated a full set of kinetic constants and carried out numerical kinetic analysis of the process. Thus, we were able to compare directly the physicochemical parameters of the reaction obtained in silico and the macroscopic parameters observed in experimental studies. From our computational results, we explain the observed unusual dependence of the rate constant of covalent complex formation, kobs, on the ARS concentration. The latter depends both on the non-covalent binding step with the equilibrium constant, Ki, and on the rate constant of covalent adduct formation, kinact. The calculated ratio kinact/Ki = 213 M-1 s-1 reproduces the corresponding experimental value of 250 ± 30 M-1 s-1 for the interaction of ARS-853 with KRASG12C. Electron density analysis in the reactive region demonstrates that covalent bond formation occurs efficiently according to the Michael addition mechanism, which assumes the activation of the C[double bond, length as m-dash]C bond of ARS-853 by a water molecule and Lys16 in the binding site of KRASG12C. We also refine the kinact and Ki constants of the ARS-107 compound, which shares common features with ARS-853, and show that the decrease in the kinact/Ki ratio in the case of ARS-107 is explained by changes in both Ki and kinact constants.

Dynamical properties of enzyme-substrate complexes disclose substrate specificity of the SARS-CoV-2 main protease as characterized by the electron density descriptors

A dynamical approach is proposed to discriminate between reactive (rES) and nonreactive (nES) enzyme-substrate complexes taking the SARS-CoV-2 main protease (Mpro) as an important example. Molecular dynamics simulations with the quantum mechanics/molecular mechanics potentials (QM(DFT)/MM-MD) followed by the electron density analysis are employed to evaluate geometry and electronic properties of the enzyme with different substrates along MD trajectories. We demonstrate that mapping the Laplacian of the electron density and the electron localization function provides easily visible images of the substrate activation that allow one to distinguish rES and nES. The computed fractions of reactive enzyme-substrate complexes along MD trajectories well correlate with the findings of recent experimental studies on the substrate specificity of Mpro. The results of our simulations demonstrate the role of the theory level used in QM subsystems for a proper description of the nucleophilic attack of the catalytic cysteine residue in Mpro. The activation of the carbonyl group of a substrate is correctly characterized with the hybrid DFT functional PBE0, whereas the use of a GGA-type PBE functional, that lacks the admixture of the Hartree-Fock exchange fails to describe activation.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Novel flavin-based fluorescent proteins with red-shifted emission bands: a computational study

The iLOV protein is a promising member of the class of flavin mononucleotide (FMN) based fluorescent proteins (FbFPs). It is becoming a popular tool for bioanalytical applications and bioimaging as a competitor of the well-known green fluorescent protein and its analogues. The main limitation of FbFPs is that all the members have close values of their absorption and emission band maxima. Therefore the upcoming challenge is to introduce novel variants of FbFPs to extend their color palette. We report the results of computational studies of iLOV variants, introducing point mutations and chromophore analogues. We found that point mutations of the apoprotein and substitution of FMN with either 8-amino-FMN or 8-methylamino-FMN lead to the red shift of emission bands up to 100 nm. Substitution with 1-deaza-FMN and the point mutations of the apoprotein result in a set of novel fluorescent proteins with emission bands in the "transparent" window where light readily penetrates through mammalian tissues. Newly suggested FbFPs can be used for multicolor imaging and also as components of FRET pairs.

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.

Competition between two cysteines in covalent binding of biliverdin to phytochrome domains

In this work, we disclose a mechanism of competing chemical reactions of protein assembly for a bacterial phytochrome using modern methods of molecular modeling.

[(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).

Distinct Protein Hydration Water Species Defined by Spatially Resolved Spectra of Intermolecular Vibrations

In this molecular dynamics simulation study, we analyze intermolecular vibrations in the hydration shell of a solvated enyzme, the membrane type 1–matrix metalloproteinase, with high spatial resolution. Our approach allows us to characterize vibrational signatures of the local hydrogen bond network, the translational mobility of water molecules, as well as the molecular entropy, in specific local environments. Our study demonstrates the heterogeneity of water properties within the hydration shell of a complex biomolecule. We define a classification scheme based on the vibrational density of states that allows us to distinguish separate classes of hydration water species and facilitates the description of hydration water properties at distinct hydration sites. The results demonstrate that no single characteristic of the protein surface is sufficient to determine the properties of nearby water. The protein surface geometry, quantified here by the number of protein atoms in the vicinity of a hydration water molecule, as well as the chemical nature of a solvated protein functional group, influences dynamic and thermodynamic properties of solvating water molecules.

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.