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

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.

Polysaccharide-based tumor microenvironment-responsive drug delivery systems for cancer therapy

The biochemical indicators of tumor microenvironment (TME) that are different from normal tissues provide the possibility for constructing intelligent drug delivery systems (DDSs). Polysaccharides with good biocompatibility, biodegradability, and unique biological properties are ideal materials for constructing DDSs. Nanogels, micelles, organic-inorganic nanocomposites, hydrogels, and microneedles (MNs) are common polysaccharide-based DDSs. Polysaccharide-based DDSs enable precise control of drug delivery and release processes by incorporating TME-specific biochemical indicators. The classification and design strategies of polysaccharide-based TME-responsive DDSs are comprehensively reviewed. The advantages and challenges of current polysaccharide-based DDSs are summarized and the future directions of development are foreseen. The polysaccharide-based TME-responsive DDSs are expected to provide new strategies and solutions for cancer therapy and make important contributions to the realization of precision medicine.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

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.

Optimized Recombinant Expression and Characterization of Collagenase in Bacillus subtilis WB600

Background: The collagenase encoding gene col was cloned into a pP43NMK vector and amplified in Escherichia coli JM109 cells. The shuttle vector pP43NMK was used to sub-clone the col gene to obtain the vector pP43NMK-col for the expression of collagenase in Bacillus subtilis WB600. The enzyme was characterized and the composition of the expression medium and culture conditions were optimized. Methods: The expressed recombinant enzyme was purified by ammonium sulfate, ultrafiltration, and through a nickel column. The purified collagenase had an activity of 9405.54 U/mg. Results: The recombinant enzyme exhibited optimal activity at pH 9.0 and 50 °C. Catalytic efficiency of the recombinant collagenase was inhibited by Fe3+ and Cu2+, but stimulated by Co2+, Ca2+, Zn2+, and Mg2+. The optimal conditions for its growth were at pH 7.0 and 35 °C, using 15 g/L of fructose and 36 g/L of yeast powder and peptone mixture (2:1) at 260 rpm with 11% inoculation. The maximal extracellular activity of the recombinant collagenase reached 2746.7 U/mL after optimization of culture conditions, which was 2.4-fold higher than that before optimization. Conclusions: This study is a first attempt to recombinantly express collagenase in B. subtilis WB600 and optimize its expression conditions, its production conditions, and possible scale-up.

Promiscuous Catalytic Activity of a Binuclear Metallohydrolase: Peptide and Phosphoester Hydrolyses

In this study, chemical promiscuity of a binuclear metallohydrolase Streptomyces griseus aminopeptidase (SgAP) has been investigated using DFT calculations. SgAP catalyzes two diverse reactions, peptide and phosphoester hydrolyses, using its binuclear (Zn-Zn) core. On the basis of the experimental information, mechanisms of these reactions have been investigated utilizing leucine p-nitro aniline (Leu-pNA) and bis(4-nitrophenyl) phosphate (BNPP) as the substrates. The computed barriers of 16.5 and 16.8 kcal/mol for the most plausible mechanisms proposed by the DFT calculations are in good agreement with the measured values of 13.9 and 18.3 kcal/mol for the Leu-pNA and BNPP hydrolyses, respectively. The former was found to occur through the transfer of two protons, while the latter with only one proton transfer. They are in line with the experimental observations. The cleavage of the peptide bond was the rate-determining process for the Leu-pNA hydrolysis. However, the creation of the nucleophile and its attack on the electrophile phosphorus atom was the rate-determining step for the BNPP hydrolysis. These calculations showed that the chemical nature of the substrate and its binding mode influence the nucleophilicity of the metal bound hydroxyl nucleophile. Additionally, the nucleophilicity was found to be critical for the Leu-pNA hydrolysis, whereas double Lewis acid activation was needed for the BNPP hydrolysis. That could be one of the reasons why peptide hydrolysis can be catalyzed by both mononuclear and binuclear metal cofactors containing hydrolases, while phosphoester hydrolysis is almost exclusively by binuclear metallohydrolases. These results will be helpful in the development of versatile catalysts for chemically distinct hydrolytic reactions.