1
|
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
|
2
|
Cousido-Siah A, Ruiz FX, Fanfrlík J, Giménez-Dejoz J, Mitschler A, Kamlar M, Veselý J, Ajani H, Parés X, Farrés J, Hobza P, Podjarny AD. IDD388 Polyhalogenated Derivatives as Probes for an Improved Structure-Based Selectivity of AKR1B10 Inhibitors. ACS Chem Biol 2016; 11:2693-2705. [PMID: 27359042 DOI: 10.1021/acschembio.6b00382] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Human enzyme aldo-keto reductase family member 1B10 (AKR1B10) has evolved as a tumor marker and promising antineoplastic target. It shares high structural similarity with the diabetes target enzyme aldose reductase (AR). Starting from the potent AR inhibitor IDD388, we have synthesized a series of derivatives bearing the same halophenoxyacetic acid moiety with an increasing number of bromine (Br) atoms on its aryl moiety. Next, by means of IC50 measurements, X-ray crystallography, WaterMap analysis, and advanced binding free energy calculations with a quantum-mechanical (QM) approach, we have studied their structure-activity relationship (SAR) against both enzymes. The introduction of Br substituents decreases AR inhibition potency but improves it in the case of AKR1B10. Indeed, the Br atoms in ortho position may impede these drugs to fit into the AR prototypical specificity pocket. For AKR1B10, the smaller aryl moieties of MK181 and IDD388 can bind into the external loop A subpocket. Instead, the bulkier MK184, MK319, and MK204 open an inner specificity pocket in AKR1B10 characterized by a π-π stacking interaction of their aryl moieties and Trp112 side chain in the native conformation (not possible in AR). Among the three compounds, only MK204 can make a strong halogen bond with the protein (-4.4 kcal/mol, using QM calculations), while presenting the lowest desolvation cost among all the series, translated into the most selective and inhibitory potency AKR1B10 (IC50 = 80 nM). Overall, SAR of these IDD388 polyhalogenated derivatives have unveiled several distinctive AKR1B10 features (shape, flexibility, hydration) that can be exploited to design novel types of AKR1B10 selective drugs.
Collapse
Affiliation(s)
- Alexandra Cousido-Siah
- Department
of Integrated Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire, CNRS, INSERM, UdS, 1
rue Laurent Fries 67404 CEDEX Illkirch, France
| | - Francesc X. Ruiz
- Department
of Integrated Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire, CNRS, INSERM, UdS, 1
rue Laurent Fries 67404 CEDEX Illkirch, France
- Department
of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| | - Joan Giménez-Dejoz
- Department
of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain
| | - André Mitschler
- Department
of Integrated Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire, CNRS, INSERM, UdS, 1
rue Laurent Fries 67404 CEDEX Illkirch, France
| | - Martin Kamlar
- Department
of Organic Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - Jan Veselý
- Department
of Organic Chemistry, Charles University in Prague, Hlavova 2030, 128 43 Prague 2, Czech Republic
| | - Haresh Ajani
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, Olomouc, 771 46 Olomouc, Czech Republic
| | - Xavier Parés
- Department
of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain
| | - Jaume Farrés
- Department
of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, Olomouc, 771 46 Olomouc, Czech Republic
| | - Alberto D. Podjarny
- Department
of Integrated Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire, CNRS, INSERM, UdS, 1
rue Laurent Fries 67404 CEDEX Illkirch, France
| |
Collapse
|
3
|
Sato W, Hitaoka S, Inoue K, Imai M, Saio T, Uchida T, Shinzawa-Itoh K, Yoshikawa S, Yoshizawa K, Ishimori K. Energetic Mechanism of Cytochrome c-Cytochrome c Oxidase Electron Transfer Complex Formation under Turnover Conditions Revealed by Mutational Effects and Docking Simulation. J Biol Chem 2016; 291:15320-31. [PMID: 27226541 DOI: 10.1074/jbc.m115.708065] [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: 12/03/2015] [Indexed: 11/06/2022] Open
Abstract
Based on the mutational effects on the steady-state kinetics of the electron transfer reaction and our NMR analysis of the interaction site (Sakamoto, K., Kamiya, M., Imai, M., Shinzawa-Itoh, K., Uchida, T., Kawano, K., Yoshikawa, S., and Ishimori, K. (2011) Proc. Natl. Acad. Sci. U.S.A. 108, 12271-12276), we determined the structure of the electron transfer complex between cytochrome c (Cyt c) and cytochrome c oxidase (CcO) under turnover conditions and energetically characterized the interactions essential for complex formation. The complex structures predicted by the protein docking simulation were computationally selected and validated by the experimental kinetic data for mutant Cyt c in the electron transfer reaction to CcO. The interaction analysis using the selected Cyt c-CcO complex structure revealed the electrostatic and hydrophobic contributions of each amino acid residue to the free energy required for complex formation. Several charged residues showed large unfavorable (desolvation) electrostatic interactions that were almost cancelled out by large favorable (Columbic) electrostatic interactions but resulted in the destabilization of the complex. The residual destabilizing free energy is compensated by the van der Waals interactions mediated by hydrophobic amino acid residues to give the stabilized complex. Thus, hydrophobic interactions are the primary factors that promote complex formation between Cyt c and CcO under turnover conditions, whereas the change in the electrostatic destabilization free energy provides the variance of the binding free energy in the mutants. The distribution of favorable and unfavorable electrostatic interactions in the interaction site determines the orientation of the binding of Cyt c on CcO.
Collapse
Affiliation(s)
- Wataru Sato
- From the Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628
| | - Seiji Hitaoka
- the Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka 819-0315
| | - Kaoru Inoue
- the Division of Chemistry, Graduate School of Science, and
| | - Mizue Imai
- From the Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628
| | - Tomohide Saio
- From the Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628, Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, and
| | - Takeshi Uchida
- From the Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628, Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, and
| | - Kyoko Shinzawa-Itoh
- the Department of Life Science, Graduate School of Life Science, University of Hyogo, Hyogo 678-1297, Japan
| | - Shinya Yoshikawa
- the Department of Life Science, Graduate School of Life Science, University of Hyogo, Hyogo 678-1297, Japan
| | - Kazunari Yoshizawa
- the Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka 819-0315
| | - Koichiro Ishimori
- From the Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628, Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, and
| |
Collapse
|
4
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
5
|
Yoshida T, Hayashi T, Mashima A, Sasahara K, Chuman H. A simple and efficient dispersion correction to the Hartree–Fock theory (3): A comprehensive performance comparison of HF-Dtq with MP2 and DFT-Ds. Bioorg Med Chem Lett 2016; 26:589-593. [DOI: 10.1016/j.bmcl.2015.11.061] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2015] [Revised: 11/14/2015] [Accepted: 11/19/2015] [Indexed: 11/25/2022]
|
6
|
Hitaoka S, Chuman H, Yoshizawa K. A QSAR study on the inhibition mechanism of matrix metalloproteinase-12 by arylsulfone analogs based on molecular orbital calculations. Org Biomol Chem 2015; 13:793-806. [DOI: 10.1039/c4ob01843e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The inhibition mechanism of matrix metalloproteinase-12 by arylsulfone analogs is revealed using a comprehensive computational approach including docking simulations, molecular orbital calculations, and QSAR.
Collapse
Affiliation(s)
- Seiji Hitaoka
- Institute for Materials Chemistry and Engineering and International Research Center for Molecular System
- Kyushu University
- Fukuoka 819-0395
- Japan
| | - Hiroshi Chuman
- Institute of Health Biosciences
- The University of Tokushima Graduate School
- Tokushima 770-8505
- Japan
| | - Kazunari Yoshizawa
- Institute for Materials Chemistry and Engineering and International Research Center for Molecular System
- Kyushu University
- Fukuoka 819-0395
- Japan
| |
Collapse
|
7
|
Multiscale quantum chemical approaches to QSAR modeling and drug design. Drug Discov Today 2014; 19:1921-7. [DOI: 10.1016/j.drudis.2014.09.024] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2014] [Revised: 08/01/2014] [Accepted: 09/26/2014] [Indexed: 12/24/2022]
|
8
|
Mashima A, Kurahashi M, Sasahara K, Yoshida T, Chuman H. Connecting Classical QSAR and LERE Analyses Using Modern Molecular Calculations, LERE-QSAR (VI): Hydrolysis of Substituted Hippuric Acid Phenyl Esters by Trypsin. Mol Inform 2014; 33:802-14. [PMID: 27485426 DOI: 10.1002/minf.201400099] [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: 07/14/2014] [Accepted: 09/14/2014] [Indexed: 11/08/2022]
Abstract
The reaction mechanism of trypsin was studied by applying DFT and ab initio molecular orbital (MO) calculations to complexes of trypsin with a congeneric series of eight para-substituted hippuric acid phenyl esters, for which a previous quantitative structureactivity relationship (QSAR) study revealed nice linearity of Hammett substitution constant σ(-) with logarithmic values of the MichaelisMenten and catalytic rate constants. Based on the LERE procedure, we performed QSAR analyses on each elementary reaction step during the acylation process. The present calculations showed that the rate-determining step during the acylation process is the transition state (TS) between the enzymesubstrate complex (ES) and tetrahedral intermediate (TET), and that the proton transfer occurs from Ser195 to His57, not between His57 and Asp102. The LERE-QSAR analysis statistically suggested that the variation of overall free-energy changes leading to formation of TS is governed mostly by that of activation energies required to form TS from ES. In spite of a very limited number of congeneric ligands in the current work, it is critically essential to clarify and verify physicochemical meanings of a typical QSAR/Chemoinformatics parameter, Hammett σ(-) based on quantum chemical calculations on the proteinligand kinetics; how Hammett σ(-) behaves in terms of proteinligand interaction energies.
Collapse
Affiliation(s)
- Akira Mashima
- Institute of Health Biosciences, The University of Tokushima Graduate School, 1-78 Shomachi, Tokushima 770-8505, Japan phone/fax: +81-88-633-7257/+81-88-633-9508
| | - Masahiro Kurahashi
- Institute of Health Biosciences, The University of Tokushima Graduate School, 1-78 Shomachi, Tokushima 770-8505, Japan phone/fax: +81-88-633-7257/+81-88-633-9508
| | - Katsunori Sasahara
- Institute of Health Biosciences, The University of Tokushima Graduate School, 1-78 Shomachi, Tokushima 770-8505, Japan phone/fax: +81-88-633-7257/+81-88-633-9508
| | - Tatsusada Yoshida
- Institute of Health Biosciences, The University of Tokushima Graduate School, 1-78 Shomachi, Tokushima 770-8505, Japan phone/fax: +81-88-633-7257/+81-88-633-9508
| | - Hiroshi Chuman
- Institute of Health Biosciences, The University of Tokushima Graduate School, 1-78 Shomachi, Tokushima 770-8505, Japan phone/fax: +81-88-633-7257/+81-88-633-9508.
| |
Collapse
|
9
|
Khrenova M, Savitsky AP, Topol IA, Nemukhin AV. Exploration of the zinc finger motif in controlling activity of matrix metalloproteinases. J Phys Chem B 2014; 118:13505-12. [PMID: 25375834 PMCID: PMC4254000 DOI: 10.1021/jp5088702] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Revised: 11/05/2014] [Indexed: 12/19/2022]
Abstract
Discovering ways to control the activity of matrix metalloproteinases (MMPs), zinc-dependent enzymes capable of degrading extracellular matrix proteins, is an important field of cancer research. We report here a novel strategy for assembling MMP inhibitors on the basis of oligopeptide ligands by exploring the pattern known as the zinc finger motif. Advanced molecular modeling tools were used to characterize the structural binding motifs of experimentally tested MMP inhibitors, as well as those of newly proposed peptidomimetics, in their zinc-containing active sites. The results of simulations based on the quantum mechanics/molecular mechanics (QM/MM) approach and Car-Parrinello molecular dynamics with QM/MM potentials demonstrate that, upon binding of Regasepin1, a known MMP-9 inhibitor, the Zn(2+)(His3) structural element is rearranged to the Zn(2+)(Cys2His2) zinc finger motif, in which two Cys residues are borrowed from the ligand. Following consideration of the crystal structure of MMP-2 with its inhibitor, the oligopeptide APP-IP, we proposed a new peptidomimetic with two replacements in the substrate, Tyr3Cys and Asp6Cys. Simulations show that this peptide variant blocks an enzyme active site by the Zn(2+)(Cys2His2) zinc finger construct. Similarly, a natural substrate of MMP-2, Ace-Gln-Gly ∼ Ile-Ala-Gly-Nme, can be converted to an inhibiting compound by two replacements, Ile by Cys and Gly by the d isomer of Cys, favoring formation of the zinc finger motif.
Collapse
Affiliation(s)
- Maria
G. Khrenova
- A.N. Bach Institute
of Biochemistry of the Russian Academy of Science, Leninsky Prospect, 33, Moscow 119071, Russian Federation
- Chemistry
Department, M.V. Lomonosov Moscow State
University, Leninskie
Gory 1/3, Moscow, 119991, Russian Federation
| | - Alexander P. Savitsky
- A.N. Bach Institute
of Biochemistry of the Russian Academy of Science, Leninsky Prospect, 33, Moscow 119071, Russian Federation
| | - Igor A. Topol
- Advanced
Biomedical Computing Center, Information Systems Program, Leidos Biomedical
Research Inc., Frederick National Laboratory
for Cancer Research, Frederick, Maryland 21702, United States
| | - 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
|
10
|
A simple and efficient dispersion correction to the Hartree–Fock theory. Bioorg Med Chem Lett 2014; 24:1037-42. [DOI: 10.1016/j.bmcl.2014.01.020] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2013] [Revised: 01/08/2014] [Accepted: 01/08/2014] [Indexed: 11/18/2022]
|
11
|
Zhu T, Xiao X, Ji C, Zhang JZH. A New Quantum Calibrated Force Field for Zinc-Protein Complex. J Chem Theory Comput 2013; 9:1788-98. [PMID: 26587635 DOI: 10.1021/ct301091z] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A quantum calibrated polarizable-charge transfer force field (QPCT) has been proposed to accurately describe the interaction dynamics of zinc-protein complexes. The parameters of the QPCT force field were calibrated by quantum chemistry calculation and capture the polarization and charge transfer effect. QPCTs are validated by molecular dynamic simulation of the hydration shell of the zinc ion, five proteins containing the most common zinc-binding sites (ZnCys2His2, ZnCys3His1, ZnCys4, Zn2Cys6), as well as protein-ligand binding energy in zinc protein MMP3. The calculated results show excellent agreement with the experimental measurement and with results from QM/MM simulation, demonstrating that QPCT is accurate enough to maintain the correct structural integrity of the zinc binding pocket and provide accurate interaction dynamics of the zinc-residue complex. The current approach can also be extended to the study of interaction dynamics of other metal-containing proteins by recalibrating the corresponding parameters to the specific complexes.
Collapse
Affiliation(s)
- Tong Zhu
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | - Xudong Xiao
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
| | - Changge Ji
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
| | - John Z H Zhang
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|