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

An 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 160,000 atoms

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 160,000 atoms

MULTISCALE MOLECULAR MODELING STUDIES OF THE DYNAMICS AND CATALYTIC MECHANISMS OF IRON(II)- AND ZINC(II)-DEPENDENT METALLOENZYMES

Enzymes are biological systems that aid in specific biochemical reactions. They lower the reaction barrier, thus speeding up the reaction rate. A detailed knowledge of enzymes will not be achievable without computational modeling as it offers insight into atomistic details and catalytic species, which are crucial to designing enzyme-specific inhibitors and impossible to gain experimentally. This dissertation employs advanced multiscale computational approaches to study the dynamics and reaction mechanisms of non-heme Fe(II) and 2-oxoglutarate (2OG) dependent oxygenases, including AlkB, AlkBH2, TET2, and KDM4E, involved in DNA and histone demethylation. It also focuses on Zn(II) dependent matrix metalloproteinase-1 (MMP-1), which helps collagen degradation. Chapter 2 investigates the substrate selectivity and dynamics on the enzyme-substrate complexes of DNA repair enzymes, AlkB and FTO. Chapter 3 unravels the mechanisms and effects of dynamics on the demethylation of 3-methylcytosine substrate by AlkB and AlkBH2 enzymes. The results imply that the nature of DNA and conformational dynamics influence the electronic structure of the iron center during demethylation. Chapter 4 delineates how second-coordination and long-range residue mutations affect the oxidation of 5-methylcytosine substrate to 5-hydroxymethylcytosine by TET2 enzyme. The results reveal that mutations affect DNA binding/interactions and the energetic contributions of residues stabilizing key catalytic species. Chapter 5 describes the reparation of unnatural alkylated substrates by TET2, their effects on second-coordination interactions and long-range correlated motions in TET2. The study reveals that post-hydroxylation reactions occur in aqueous solution outside the enzyme environment. Chapter 6 establishes how applying external electric fields (EEFs) enhances specificity of KDM4E for C—H over N—H activation during dimethylated arginine substrate demethylation. The results reveal that applying positive EEFs parallel to Fe=O bond enhances C—H activation rate, while inhibiting the N—H one. Chapter 7 addresses the formation of catalytically competent MMP-1·THP complex of MMP-1. The studies reveal the role of MMP-1’s catalytic domain a-helices, the linker, and changes in coordination states of catalytic Zn(II) during the transition. Overall, the presented results contribute to the in-depth understanding of the fundamental mechanisms of the studied enzymes and provide a background for developing enzyme-specific inhibitors against the associated disorders and diseases