1
|
Bím D, Navrátil M, Gutten O, Konvalinka J, Kutil Z, Culka M, Navrátil V, Alexandrova AN, Bařinka C, Rulíšek L. Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II. J Phys Chem B 2022; 126:132-143. [PMID: 34978450 DOI: 10.1021/acs.jpcb.1c09240] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol-1 can be inferred for the known kcat ≈ 1 s-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (kcat = 8.6 × 10-5 s-1 to 2.7 s-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.
Collapse
Affiliation(s)
- Daniel Bím
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Michal Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Jan Konvalinka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Biochemistry, Faculty of Science, Charles University, Hlavova 2030, 2120 00 Prague, Czech Republic
| | - Zsófia Kutil
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Martin Culka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Václav Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Cyril Bařinka
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
2
|
Barbosa ED, Lima Neto JX, Bezerra KS, Oliveira JIN, Machado LD, Fulco UL. Quantum Biochemical Investigation of Lys49-PLA 2 from Bothrops moojeni. J Phys Chem B 2021; 125:12972-12980. [PMID: 34793159 DOI: 10.1021/acs.jpcb.1c07298] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Envenomation via snakebites occurs largely in areas where it is harder to access the hospital. Its mortality rate and sequelae acquired by the survivors symbolize a big challenge for antivenom therapy. In particular, the homologous phospholipase A2 (Lys49-PLA2) proteins can induce myonecrosis and are not effectively neutralized by current treatments. Thus, by taking advantage of crystallographic structures of Bothrops moojeni Lys49-PLA2 complexed with VRD (varespladib) and AIN (aspirin), a quantum biochemistry study based on the molecular fractionation with conjugate cap scheme within the density functional theory formalism is performed to unveil these complexes' detailed interaction energies. The calculations revealed that important interactions between ligands and the Lys49-PLA2 pocket could occur up to a pocket radius of r = 6.5 (5.0 Å) for VRD (AIN), with the total interaction energy of the VRD ligand being higher than that of the AIN ligand, which is well-correlated with the experimental binding affinity. Furthermore, we have identified the role played by the amino acids LYS0069, LYS0049, LEU0005, ILE0009, CYS0029, GLY0030, HIS0048, PRO0018, ALA0019, CYS0045, TYR0052, TYR0022, PRO0125*, and PHE0126* (LYS0069, LYS0049, GLY0032, LEU0002, and LEU0005) in the VRD↔Lys49-PLA2 (AIN↔Lys49-PLA2) complex. Our simulations are a valuable tool to support the big challenge for neutralizing the damages in victims of snakebites.
Collapse
Affiliation(s)
- E D Barbosa
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| | - J X Lima Neto
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| | - K S Bezerra
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| | - J I N Oliveira
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| | - L D Machado
- Departamento de Física Teórica e Experimental, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| | - U L Fulco
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal 59072-970, Rio Grande do Norte, Brazil
| |
Collapse
|
3
|
Barbosa ED, Neto JXL, Teixeira DG, Bezerra KS, do Amaral VS, Oliveira JIN, Lima JPMS, Machado LD, Fulco UL. Exploring human porphobilinogen synthase metalloprotein by quantum biochemistry and evolutionary methods. Metallomics 2021; 13:6206860. [PMID: 33791795 DOI: 10.1093/mtomcs/mfab017] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 03/16/2021] [Accepted: 03/21/2021] [Indexed: 12/16/2022]
Abstract
Previous studies have shown the porphobilinogen synthase (PBGS) zinc-binding mechanism and its conservation among the living cells. However, the precise molecular interaction of zinc with the active center of the enzyme is unknown. In particular, quantum chemistry techniques within the density functional theory (DFT) framework have been the key methodology to describe metalloproteins, when one is looking for a compromise between accuracy and computational feasibility. Considering this, we used DFT-based models within the molecular fractionation with conjugate caps scheme to evaluate the binding energy features of zinc interacting with the human PBGS. Besides, phylogenetic and clustering analyses were successfully employed in extracting useful information from protein sequences to identify groups of conserved residues that build the ions-binding site. Our results also report a conservative assessment of the relevant amino acids, as well as the benchmark analysis of the calculation models used. The most relevant intermolecular interactions in Zn2+-PBGS are due to the amino acids CYS0122, CYS0124, CYS0132, ASP0169, SER0168, ARG0221, HIS0131, ASP0120, GLY0133, VAL0121, ARG0209, and ARG0174. Among these residues, we highlighted ASP0120, GLY0133, HIS0131, SER0168, and ARG0209 by co-occurring in all clusters generated by unsupervised clustering analysis. On the other hand, the triple cysteines at 2.5 Å from zinc (CYS0122, CYS0124, and CYS0132) have the highest energy attraction and are absent in the taxa Viridiplantae, Sar, Rhodophyta, and some Bacteria. Additionally, the performance of the DFT-based models shows that the processing time-dependence is more associated with the choice of the basis set than the exchange-correlation functional.
Collapse
Affiliation(s)
- E D Barbosa
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - J X Lima Neto
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - D G Teixeira
- Institute of Tropical Medicine, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - K S Bezerra
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - V S do Amaral
- Departamento de Bioquímica, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - J I N Oliveira
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - J P M Santos Lima
- Departamento de Bioquímica, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - L D Machado
- Departamento de Física Teórica e Experimental, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| | - U L Fulco
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, 59072-970 Natal-RN, Brazil
| |
Collapse
|
4
|
Long Z. Computational analysis of the metal selectivity of matrix metalloproteinase 8. PLoS One 2020; 15:e0243321. [PMID: 33275641 PMCID: PMC7717551 DOI: 10.1371/journal.pone.0243321] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Accepted: 11/18/2020] [Indexed: 11/18/2022] Open
Abstract
Matrix metalloproteinase (MMP) is a class of metalloenzyme that cleaves peptide bonds in extracellular matrices. Their functions are important in both health and disease of animals. Here using quantum mechanics simulations of the MMP8 protein, the coordination chemistry of different metal cofactors is examined. Structural comparisons reveal that Jhan-Teller effects induced by Cu(II) coordination distorts the wild-type MMP8 active site corresponding to a significant reduction in activity observed in previous experiments. In addition, further analysis suggests that a histidine to glutamine mutation at residue number 197 can potentially allow the MMP8 protein to utilize Cu(II) in reactions. Simulations also demonstrates the requirement of a conformational change in the ligand before enzymatic cleavage. The insights provided here will assist future protein engineering efforts utilizing the MMP8 protein.
Collapse
Affiliation(s)
- Zheng Long
- Department of Chemistry and Biochemistry, University of California San Diego, San Diego, California, United States of America
- * E-mail:
| |
Collapse
|