1
|
Gorantla K, Krishnan A, Waheed SO, Varghese A, DiCastri I, LaRouche C, Paik M, Fields GB, Karabencheva-Christova TG. Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1. Biochemistry 2024; 63:1925-1940. [PMID: 38963231 PMCID: PMC11309001 DOI: 10.1021/acs.biochem.4c00076] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Revised: 06/14/2024] [Accepted: 06/24/2024] [Indexed: 07/05/2024]
Abstract
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.
Collapse
Affiliation(s)
- Koteswara
Rao Gorantla
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Anandhu Krishnan
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sodiq O. Waheed
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Ann Varghese
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Isabella DiCastri
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Ciara LaRouche
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Meredith Paik
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Gregg B. Fields
- Department
of Chemistry and Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida 33458, United States
| | | |
Collapse
|
2
|
Deng WH, Liao RZ. Cysteine Radical and Glutamate Collaboratively Enable C-H Bond Activation and C-N Bond Cleavage in a Glycyl Radical Enzyme HplG. J Chem Inf Model 2024; 64:4168-4179. [PMID: 38745447 DOI: 10.1021/acs.jcim.4c00122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
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.
Collapse
Affiliation(s)
- Wen-Hao Deng
- Key Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Rong-Zhen Liao
- Key Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| |
Collapse
|
3
|
Grigorenko BL, Khrenova MG, Jones DD, Nemukhin AV. Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins. Org Biomol Chem 2024; 22:337-347. [PMID: 38063860 DOI: 10.1039/d3ob01450a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
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.
Collapse
Affiliation(s)
- Bella L Grigorenko
- Chemistry Department, Lomonosov Moscow State University, Moscow, Russian Federation.
- Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russian Federation
| | - Maria G Khrenova
- Chemistry Department, Lomonosov Moscow State University, Moscow, Russian Federation.
- Bach Institute of Biochemistry, Moscow, Russian Federation
| | - D Dafydd Jones
- School of Biosciences, Molecular Biosciences Division, Cardiff University, Cardiff, UK
| | - Alexander V Nemukhin
- Chemistry Department, Lomonosov Moscow State University, Moscow, Russian Federation.
- Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russian Federation
| |
Collapse
|
4
|
Varghese A, Waheed SO, Gorantla K, DiCastri I, LaRouche C, Kaski B, Fields GB, Karabencheva-Christova TG. Catalytic Mechanism of Collagen Hydrolysis by Zinc(II)-Dependent Matrix Metalloproteinase-1. J Phys Chem B 2023; 127:9697-9709. [PMID: 37931179 PMCID: PMC10659029 DOI: 10.1021/acs.jpcb.3c04293] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 10/16/2023] [Accepted: 10/18/2023] [Indexed: 11/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ann Varghese
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sodiq O. Waheed
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Koteswararao Gorantla
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Isabella DiCastri
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Ciara LaRouche
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Brendan Kaski
- Department
of Kinesiology and Integrative Physiology, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Gregg B. Fields
- Department
of Chemistry and Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida 33458, United States
| | | |
Collapse
|
5
|
Choi JY, Chung E. Molecular Dynamics Simulations of Matrix Metalloproteinase 13 and the Analysis of the Specificity Loop and the S1'-Site. Int J Mol Sci 2023; 24:10577. [PMID: 37445757 PMCID: PMC10342107 DOI: 10.3390/ijms241310577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2023] [Revised: 05/05/2023] [Accepted: 06/20/2023] [Indexed: 07/15/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Jun Yong Choi
- Department of Chemistry and Biochemistry, Queens College, Flushing, NY 11367, USA
- Ph.D. Programs in Chemistry and Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, USA
| | - Eugene Chung
- Department of Chemistry and Biochemistry, Queens College, Flushing, NY 11367, USA
| |
Collapse
|
6
|
Castro-Amorim J, Oliveira A, Mukherjee AK, Ramos MJ, Fernandes PA. Unraveling the Reaction Mechanism of Russell's Viper Venom Factor X Activator: A Paradigm for the Reactivity of Zinc Metalloproteinases? J Chem Inf Model 2023. [PMID: 37092784 DOI: 10.1021/acs.jcim.2c01156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
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 Zn2+ in catalysis by lowering the pKa of the Zn2+-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 Zn2+ cofactor, this step would be rate-limiting.
Collapse
Affiliation(s)
- Juliana Castro-Amorim
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, Porto 4169-007, Portugal
| | - Ana Oliveira
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, Porto 4169-007, Portugal
| | - Ashis K Mukherjee
- Institute of Advanced Study in Science and Technology, Vigyan Path Garchuk, Paschim Boragaon, Guwahati 781035, Assam, India
| | - Maria J Ramos
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, Porto 4169-007, Portugal
| | - Pedro A Fernandes
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, Porto 4169-007, Portugal
| |
Collapse
|
7
|
Waheed SO, Varghese A, DiCastri I, Kaski B, LaRouche C, Fields GB, Karabencheva-Christova TG. Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1. Chemphyschem 2023; 24:e202200649. [PMID: 36161746 DOI: 10.1002/cphc.202200649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2022] [Revised: 09/23/2022] [Indexed: 02/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sodiq O Waheed
- Department of Chemistry, Michigan Technological University, Houghton, Michigan, 49931, USA
| | - Ann Varghese
- Department of Chemistry, Michigan Technological University, Houghton, Michigan, 49931, USA
| | - Isabella DiCastri
- Department of Chemistry, Michigan Technological University, Houghton, Michigan, 49931, USA
| | - Brenden Kaski
- Department of Kinesiology and Integrative Physiology, Michigan Technological University, Houghton, Michigan, 49931, USA
| | - Ciara LaRouche
- Department of Chemical Engineering, Michigan Technological University, Houghton, Michigan, 49931, USA
| | - Gregg B Fields
- Department of Chemistry & Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida, 33458, USA
| | | |
Collapse
|
8
|
Tomić A, Tomić S. Demystifying DPP III Catalyzed Peptide Hydrolysis—Computational Study of the Complete Catalytic Cycle of Human DPP III Catalyzed Tynorphin Hydrolysis. Int J Mol Sci 2022; 23:ijms23031858. [PMID: 35163780 PMCID: PMC8836397 DOI: 10.3390/ijms23031858] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 01/28/2022] [Accepted: 02/03/2022] [Indexed: 12/04/2022] Open
Abstract
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.
Collapse
|
9
|
Varghese A, Chaturvedi SS, DiCastri B, Mehler E, Fields GB, Karabencheva-Christova TG. 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. Chemphyschem 2021; 23:10.1002/cphc.202100680. [PMID: 35991515 PMCID: PMC9387770 DOI: 10.1002/cphc.202100680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Indexed: 11/06/2022]
Abstract
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.
Collapse
Affiliation(s)
- Ann Varghese
- Department of Chemistry, Michigan Technological University, Houghton, Michigan 49931
| | - Shobhit S Chaturvedi
- Department of Chemistry, Michigan Technological University, Houghton, Michigan 49931
| | - Bella DiCastri
- Department of Chemistry, Michigan Technological University, Houghton, Michigan 49931
| | - Emerald Mehler
- Department of Chemical Engineering, Michigan Technological University, Houghton, Michigan 49931
| | - Gregg B Fields
- Department of Chemistry and Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida 33458
| | | |
Collapse
|
10
|
Chen B, Kang Z, Zheng E, Liu Y, Gauld JW, Wang Q. Hydrolysis Mechanism of the Linkers by Matrix Metalloproteinase-9 Using QM/MM Calculations. J Chem Inf Model 2021; 61:5203-5211. [PMID: 34649435 DOI: 10.1021/acs.jcim.1c00825] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Binbin Chen
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Zhengzhong Kang
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - En Zheng
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Yingchun Liu
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - James W Gauld
- Department of Chemistry and Biochemistry, University of Windsor, Windsor N9B 3P4, Canada
| | - Qi Wang
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
11
|
Krivitskaya AV, Khrenova MG. Boronic Acids as Prospective Inhibitors of Metallo-β-Lactamases: Efficient Chemical Reaction in the Enzymatic Active Site Revealed by Molecular Modeling. Molecules 2021; 26:2026. [PMID: 33918209 PMCID: PMC8038151 DOI: 10.3390/molecules26072026] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Revised: 03/31/2021] [Accepted: 04/01/2021] [Indexed: 01/08/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Alexandra V. Krivitskaya
- Bach Institute of Biochemistry, Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences, 119071 Moscow, Russia;
| | - Maria G. Khrenova
- Bach Institute of Biochemistry, Federal Research Centre “Fundamentals of Biotechnology” of the Russian Academy of Sciences, 119071 Moscow, Russia;
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| |
Collapse
|
12
|
Elsässer B, Goettig P. Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies. Int J Mol Sci 2021; 22:3232. [PMID: 33810118 PMCID: PMC8004986 DOI: 10.3390/ijms22063232] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 02/24/2021] [Accepted: 02/26/2021] [Indexed: 12/11/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Peter Goettig
- Structural Biology Group, Department of Biosciences, University of Salzburg, Billrothstrasse 11, 5020 Salzburg, Austria;
| |
Collapse
|
13
|
De-shielding of activatable cell-penetrating peptides: recognizing and releasing in activation process. RESEARCH ON CHEMICAL INTERMEDIATES 2021. [DOI: 10.1007/s11164-020-04339-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
14
|
Palermo G, Spinello A, Saha A, Magistrato A. Frontiers of metal-coordinating drug design. Expert Opin Drug Discov 2020; 16:497-511. [PMID: 33874825 DOI: 10.1080/17460441.2021.1851188] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, Riverside, United States
| | - Angelo Spinello
- National Research Council (CNR) of Italy, Institute of Material (IOM) @ International School for Advanced Studies (SISSA), Trieste, Italy
| | - Aakash Saha
- Department of Bioengineering, University of California Riverside, Riverside, United States
| | - Alessandra Magistrato
- National Research Council (CNR) of Italy, Institute of Material (IOM) @ International School for Advanced Studies (SISSA), Trieste, Italy
| |
Collapse
|
15
|
Khrenova MG, Kulakova AM, Nemukhin AV. Proof of concept for poor inhibitor binding and efficient formation of covalent adducts of KRAS G12C and ARS compounds. Org Biomol Chem 2020; 18:3069-3081. [PMID: 32101243 DOI: 10.1039/d0ob00071j] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
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.
Collapse
Affiliation(s)
- Maria G Khrenova
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation. and Bach Institute of Biochemistry, Research Center of Biotechnology, Russian Academy of Sciences, Moscow, 119071, Russian Federation
| | - Anna M Kulakova
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation.
| | - Alexander V Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation. and Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, 119334, Russian Federation
| |
Collapse
|
16
|
Khrenova MG, Tsirelson VG, Nemukhin AV. Dynamical properties of enzyme-substrate complexes disclose substrate specificity of the SARS-CoV-2 main protease as characterized by the electron density descriptors. Phys Chem Chem Phys 2020; 22:19069-19079. [PMID: 32812956 DOI: 10.1039/d0cp03560b] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
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.
Collapse
Affiliation(s)
- Maria G Khrenova
- Bach Institute of Biochemistry, Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Leninsky Prospect, 33, bld. 2, Moscow, 119071, Russia and Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia.
| | - Vladimir G Tsirelson
- Mendeleev University of Chemical Technology, Miusskaya Square, 9, Moscow, 125047, Russia
| | - Alexander V Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia. and Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, 119334, Russia
| |
Collapse
|
17
|
Khrenova MG, Nemukhin AV, Tsirelson VG. Discrimination of enzyme–substrate complexes by reactivity using the electron density analysis: peptide bond hydrolysis by the matrix metalloproteinase-2. MENDELEEV COMMUNICATIONS 2020. [DOI: 10.1016/j.mencom.2020.09.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
18
|
Yao C, Chen B, Kang Z, Liu Y, Qi X, Wang Q. Binding, selectivity and sequence recognition of matrix metalloproteinase-2 to oligopeptides. J Biomol Struct Dyn 2019; 38:275-282. [DOI: 10.1080/07391102.2019.1578694] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Cai Yao
- Department of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou, China
| | - Binbin Chen
- Department of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou, China
| | - Zhengzhong Kang
- Department of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou, China
| | - Yingchun Liu
- Department of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou, China
| | - Xianrong Qi
- School of Pharmaceutical Sciences, Peking University, Beijing, China
| | - Qi Wang
- Department of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou, China
| |
Collapse
|
19
|
Lu Y, Farrow MR, Fayon P, Logsdail AJ, Sokol AA, Catlow CRA, Sherwood P, Keal TW. Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment. J Chem Theory Comput 2019; 15:1317-1328. [PMID: 30511845 DOI: 10.1021/acs.jctc.8b01036] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- You Lu
- Scientific Computing Department , STFC Daresbury Laboratory , Keckwick Lane, Daresbury , Warrington WA4 4AD , United Kingdom
| | - Matthew R Farrow
- Kathleen Lonsdale Materials Chemistry, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom
| | - Pierre Fayon
- Scientific Computing Department , STFC Daresbury Laboratory , Keckwick Lane, Daresbury , Warrington WA4 4AD , United Kingdom
| | - Andrew J Logsdail
- Kathleen Lonsdale Materials Chemistry, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom.,Cardiff Catalysis Institute, School of Chemistry , Cardiff University , Cardiff CF10 3AT , United Kingdom
| | - Alexey A Sokol
- Kathleen Lonsdale Materials Chemistry, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom
| | - C Richard A Catlow
- Kathleen Lonsdale Materials Chemistry, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom.,Cardiff Catalysis Institute, School of Chemistry , Cardiff University , Cardiff CF10 3AT , United Kingdom.,UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory , Harwell Science and Innovation Campus , Oxon OX11 0QX , United Kingdom
| | - Paul Sherwood
- Scientific Computing Department , STFC Daresbury Laboratory , Keckwick Lane, Daresbury , Warrington WA4 4AD , United Kingdom
| | - Thomas W Keal
- Scientific Computing Department , STFC Daresbury Laboratory , Keckwick Lane, Daresbury , Warrington WA4 4AD , United Kingdom
| |
Collapse
|
20
|
Meteleshko YI, Nemukhin AV, Khrenova MG. Novel flavin-based fluorescent proteins with red-shifted emission bands: a computational study. Photochem Photobiol Sci 2018; 18:177-189. [PMID: 30403258 DOI: 10.1039/c8pp00361k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Yulia I Meteleshko
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow, 119991, Russian Federation.
| | - Alexander V Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow, 119991, Russian Federation. and Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow, 119334, Russian Federation
| | - Maria G Khrenova
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow, 119991, Russian Federation. and Federal Research Center of Biotechnology, Bach Institute of Biochemistry, Russian Academy of Sciences, Leninskiy Prospect 33, 119071 Moscow, Russian Federation
| |
Collapse
|
21
|
Kapusta DP, Firsov DA, Khrenova MG, Grigorenko BL, Nemukhin AV. Effect of solvation water shells on enzyme active sites in zinc-dependent hydrolases. Struct Chem 2018. [DOI: 10.1007/s11224-018-1206-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
22
|
Decaneto E, Vasilevskaya T, Kutin Y, Ogata H, Grossman M, Sagi I, Havenith M, Lubitz W, Thiel W, Cox N. Solvent water interactions within the active site of the membrane type I matrix metalloproteinase. Phys Chem Chem Phys 2018; 19:30316-30331. [PMID: 28951896 DOI: 10.1039/c7cp05572b] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
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.
Collapse
Affiliation(s)
- Elena Decaneto
- Max Planck Institute for Chemical Energy Conversion, Stiftstraße. 34-36, D-45470, Mülheim an der Ruhr, Germany.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
23
|
Khrenova MG, Kulakova AM, Nemukhin AV. Competition between two cysteines in covalent binding of biliverdin to phytochrome domains. Org Biomol Chem 2018; 16:7518-7529. [DOI: 10.1039/c8ob02262c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this work, we disclose a mechanism of competing chemical reactions of protein assembly for a bacterial phytochrome using modern methods of molecular modeling.
Collapse
Affiliation(s)
- Maria G. Khrenova
- Department of Chemistry
- Lomonosov Moscow State University
- Moscow
- Russian Federation
- Federal Research Center of Biotechnology
| | - Anna M. Kulakova
- Department of Chemistry
- Lomonosov Moscow State University
- Moscow
- Russian Federation
| | - Alexander V. Nemukhin
- Department of Chemistry
- Lomonosov Moscow State University
- Moscow
- Russian Federation
- Emanuel Institute of Biochemical Physics
| |
Collapse
|
24
|
Linder DP, Baker BE, Rodgers KR. [(H2O)Zn(Imidazole)n]2+: the vital roles of coordination number and geometry in Zn–OH2 acidity and catalytic hydrolysis. Phys Chem Chem Phys 2018; 20:24979-24991. [DOI: 10.1039/c8cp03121e] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Zn(ii)–(Imidazole(ate))n coordination motif occurs in numerous biochemical systems, including carbonic anhydrase and the matrix metalloproteinases (MMPs).
Collapse
Affiliation(s)
- Douglas P. Linder
- Department of Chemistry and Physics
- Southwestern Oklahoma State University
- Weatherford
- USA
| | - Brett E. Baker
- Department of Chemistry and Physics
- Southwestern Oklahoma State University
- Weatherford
- USA
| | - Kenton R. Rodgers
- Department of Chemistry and Biochemistry
- North Dakota State University
- Fargo
- USA
| |
Collapse
|
25
|
Khrenova MG, Meteleshko YI, Nemukhin AV. Mutants of the Flavoprotein iLOV as Prospective Red-Shifted Fluorescent Markers. J Phys Chem B 2017; 121:10018-10025. [DOI: 10.1021/acs.jpcb.7b07533] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maria G. Khrenova
- Department
of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia
| | - Yulia I. Meteleshko
- Department
of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia
| | - Alexander V. Nemukhin
- Department
of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia
- Emanuel
Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, 119334, Russia
| |
Collapse
|
26
|
Grigorenko BL, Krylov AI, Nemukhin AV. Molecular Modeling Clarifies the Mechanism of Chromophore Maturation in the Green Fluorescent Protein. J Am Chem Soc 2017; 139:10239-10249. [DOI: 10.1021/jacs.7b00676] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Bella L. Grigorenko
- Department
of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
- Emanuel
Institute of Biochemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| | - Anna I. Krylov
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Alexander V. Nemukhin
- Department
of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
- Emanuel
Institute of Biochemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
| |
Collapse
|
27
|
Pattni V, Vasilevskaya T, Thiel W, Heyden M. Distinct Protein Hydration Water Species Defined by Spatially Resolved Spectra of Intermolecular Vibrations. J Phys Chem B 2017. [PMID: 28636363 PMCID: PMC5607456 DOI: 10.1021/acs.jpcb.7b03966] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
![]()
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.
Collapse
Affiliation(s)
- Viren Pattni
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, DE-45470 Mülheim an der Ruhr, Germany
| | - Tatiana Vasilevskaya
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, DE-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, DE-45470 Mülheim an der Ruhr, Germany
| | - Matthias Heyden
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, DE-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
28
|
Sparta M, Retegan M, Pinski P, Riplinger C, Becker U, Neese F. Multilevel Approaches within the Local Pair Natural Orbital Framework. J Chem Theory Comput 2017; 13:3198-3207. [DOI: 10.1021/acs.jctc.7b00260] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Manuel Sparta
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| | - Marius Retegan
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| | - Peter Pinski
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| | - Christoph Riplinger
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| | - Ute Becker
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical
Energy
Conversion, Stiftstrasse 34-36, D−45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
29
|
Ganguly A, Luong TQ, Brylski O, Dirkmann M, Möller D, Ebbinghaus S, Schulz F, Winter R, Sanchez-Garcia E, Thiel W. Elucidation of the Catalytic Mechanism of a Miniature Zinc Finger Hydrolase. J Phys Chem B 2017. [DOI: 10.1021/acs.jpcb.7b05027] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Abir Ganguly
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Trung Quan Luong
- Fakultät
für Chemie und Chemische Biologie, Physikalische Chemie I, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Oliver Brylski
- Fakultät
für Chemie und Biochemie, Physikalische Chemie II, Ruhr-Universität Bochum, 44780, Bochum, Germany
| | - Michael Dirkmann
- Fakultät
für Chemie und Biochemie, Organische Chemie I, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - David Möller
- Fakultät
für Chemie und Biochemie, Organische Chemie I, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Simon Ebbinghaus
- Fakultät
für Chemie und Biochemie, Physikalische Chemie II, Ruhr-Universität Bochum, 44780, Bochum, Germany
| | - Frank Schulz
- Fakultät
für Chemie und Biochemie, Organische Chemie I, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Roland Winter
- Fakultät
für Chemie und Chemische Biologie, Physikalische Chemie I, Technische Universität Dortmund, 44227 Dortmund, Germany
| | | | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
30
|
Khrenova MG, Solovyev ID, Lapshin GD, Savitsky AP. Molecular mechanism of interactions between MMP-2 and its oligopeptide-based inhibitors. MENDELEEV COMMUNICATIONS 2017. [DOI: 10.1016/j.mencom.2017.03.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
31
|
Karabencheva-Christova TG, Christov CZ, Fields GB. Collagenolytic Matrix Metalloproteinase Structure–Function Relationships: Insights From Molecular Dynamics Studies. STRUCTURAL AND MECHANISTIC ENZYMOLOGY 2017; 109:1-24. [DOI: 10.1016/bs.apcsb.2017.04.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
32
|
Vasilevskaya T, Khrenova MG, Nemukhin AV, Thiel W. Methodological aspects of QM/MM calculations: A case study on matrix metalloproteinase-2. J Comput Chem 2016; 37:1801-9. [PMID: 27140531 DOI: 10.1002/jcc.24395] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2016] [Revised: 04/06/2016] [Accepted: 04/07/2016] [Indexed: 01/15/2023]
Abstract
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.
Collapse
Affiliation(s)
| | - Maria G Khrenova
- Chemistry Department, Lomonosov Moscow State University, Moscow, Russia.,A.N. Bach Institute of Biochemistry, Russian Academy of Sciences, Moscow, Russia
| | - Alexander V Nemukhin
- Chemistry Department, Lomonosov Moscow State University, Moscow, Russia.,Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
33
|
Reaction mechanism of matrix metalloproteinases with a catalytically active zinc ion studied by the QM(DFTB)/MM simulations. MENDELEEV COMMUNICATIONS 2016. [DOI: 10.1016/j.mencom.2016.04.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
34
|
Kots ED, Khrenova MG, Lushchekina SV, Varfolomeev SD, Grigorenko BL, Nemukhin AV. Modeling the Complete Catalytic Cycle of Aspartoacylase. J Phys Chem B 2016; 120:4221-31. [DOI: 10.1021/acs.jpcb.6b02542] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ekaterina D. Kots
- Chemistry
Department, M.V. Lomonosov Moscow State University, Leninskie
Gory 1/3, Moscow 119991, Russian Federation
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow 119334, Russian Federation
| | - Maria G. Khrenova
- Chemistry
Department, M.V. Lomonosov Moscow State University, Leninskie
Gory 1/3, Moscow 119991, Russian Federation
| | - Sofya V. Lushchekina
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow 119334, Russian Federation
| | - Sergei D. Varfolomeev
- Chemistry
Department, M.V. Lomonosov Moscow State University, Leninskie
Gory 1/3, Moscow 119991, Russian Federation
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow 119334, Russian Federation
| | - Bella L. Grigorenko
- Chemistry
Department, M.V. Lomonosov Moscow State University, Leninskie
Gory 1/3, Moscow 119991, Russian Federation
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow 119334, Russian Federation
| | - Alexander V. Nemukhin
- Chemistry
Department, M.V. Lomonosov Moscow State University, Leninskie
Gory 1/3, Moscow 119991, Russian Federation
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow 119334, Russian Federation
| |
Collapse
|