1
|
Bowling PE, Dasgupta S, Herbert JM. Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites. J Chem Inf Model 2024; 64:3912-3922. [PMID: 38648614 DOI: 10.1021/acs.jcim.4c00221] [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: 04/25/2024]
Abstract
In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.
Collapse
Affiliation(s)
- Paige E Bowling
- Biophysics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Chemistry and Biochemistry, University of California-San Diego, La Jolla, California 92093, United States
| | - John M Herbert
- Biophysics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
2
|
Zhang M, Yang L, Ding W, Zhang H. The His23 and Lys79 pair determines the high catalytic efficiency of the inorganic pyrophosphatase of the haloacid dehalogenase superfamily. Biochim Biophys Acta Gen Subj 2022; 1866:130128. [PMID: 35278619 DOI: 10.1016/j.bbagen.2022.130128] [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: 09/25/2021] [Revised: 03/07/2022] [Accepted: 03/07/2022] [Indexed: 10/18/2022]
Abstract
Haloacid dehalogenase (HAD) superfamily members are mainly phosphomonoesterases, while BT2127 from Bacteroides thetaiotaomicron of the HAD superfamily is identified as an inorganic pyrophosphatase. In this study, to explore the roles of the Lys79 and His23 pair in the hydrolysis reaction of inorganic pyrophosphate (PPi) catalyzed by BT2127, a series of models were designed. Calculations were performed by using the density functional theory (DFT) method with the dispersion energy D3-B3LYP. The His23 and Lys79 pair plays a key role in the high catalytic efficiency of BT2127 with PPi. First, the His23 and Lys79 pair prompts Asp13 to easily provide a proton to the leaving group, which remarkably reduces the energy barrier of the phospho-transfer step; then, Lys79 provides a proton to the first leaving phosphate group via His23, produces a more electrically stabilized phosphate (H3PO4), makes this step exothermal, and further promotes the subsequent phospho-enzyme intermediate hydrolysis. The results suggest that the Lys79-His23 pair helps BT2127 reach high catalytic efficiency by strengthening the acid catalysis. Our study provides detailed chemical insights into the evolution of the inorganic pyrophosphatase function of BT2127 from the phosphomonoesterase of the HAD superfamily and the biomimetic enzyme design.
Collapse
Affiliation(s)
- Mingming Zhang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, PR China
| | - Ling Yang
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou 325001, PR China.
| | - Wanjian Ding
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, PR China.
| | - Hao Zhang
- Biomedical Research Center, College of Life Science and Engineering, Northwest Minzu University, Lanzhou, 730030, PR China.
| |
Collapse
|
3
|
Yang L, Lu Y, Tian W, Feng Y, Bai J, Zhang H. Insights into the functional divergence of the haloacid dehalogenase superfamily from phosphomonoesterase to inorganic pyrophosphatase. Arch Biochem Biophys 2021; 705:108896. [PMID: 33940035 DOI: 10.1016/j.abb.2021.108896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 04/25/2021] [Accepted: 04/26/2021] [Indexed: 11/20/2022]
Abstract
The evolution of enzyme catalytic structures and mechanisms has drawn increasing attention. In this study, we investigate the functional divergence from phosphomonoesterase to inorganic pyrophosphatase in the haloacid dehalogenase (HAD) superfamily. In this study, a series of models was constructed, and calculations were performed by using density functional theory with the B3LYP functional. The calculations suggest that in most HAD members, the active-site structure is unstable due to the binding of the substrate inorganic pyrophosphate (PPi), and reactions involving PPi cannot be catalyzed. In BT2127, which is a unique member of the HAD superfamily, the Mg2+-coordinating residues Asn172 and Glu47 play a role in stabilizing the active-site structure to adapt to the substrate PPi by providing much stronger coordination interactions with the Mg2+ ion. The calculation results suggest that Asn172 and Glu47 are crucial in the evolution of the inorganic pyrophosphatase activity in the HAD superfamily. Our study provides definitive chemical insight into the functional divergence of the HAD superfamily, and helps in understanding the evolution of enzyme catalytic structures and mechanisms.
Collapse
Affiliation(s)
- Ling Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, Institute of Theoretical and Simulation Chemistry, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin, 150001, PR China
| | - Yajie Lu
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, Institute of Theoretical and Simulation Chemistry, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin, 150001, PR China
| | - Weiquan Tian
- Chongqing Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Chongqing University, Huxi Campus, Chongqing 401331, PR China
| | - Yulan Feng
- Biomedical Research Center, College of Life Science and Engineering, Northwest Minzu University, Lanzhou, 730030, PR China
| | - Jialin Bai
- Biomedical Research Center, College of Life Science and Engineering, Northwest Minzu University, Lanzhou, 730030, PR China
| | - Hao Zhang
- Biomedical Research Center, College of Life Science and Engineering, Northwest Minzu University, Lanzhou, 730030, PR China.
| |
Collapse
|
4
|
Theoretical Studies on Catalysis Mechanisms of Serum Paraoxonase 1 and Phosphotriesterase Diisopropyl Fluorophosphatase Suggest the Alteration of Substrate Preference from Paraoxonase to DFP. Molecules 2018; 23:molecules23071660. [PMID: 29986514 PMCID: PMC6100192 DOI: 10.3390/molecules23071660] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2018] [Revised: 07/04/2018] [Accepted: 07/05/2018] [Indexed: 12/18/2022] Open
Abstract
The calcium-dependent β-propeller proteins mammalian serum paraoxonase 1 (PON1) and phosphotriesterase diisopropyl fluorophosphatase (DFPase) catalyze the hydrolysis of organophosphorus compounds and enhance hydrolysis of various nerve agents. In the present work, the phosphotriesterase activity development between PON1 and DFPase was investigated by using the hybrid density functional theory method B3LYP. Based on the active-site difference between PON1 and DFPase, both the wild type and the mutant (a water molecule replacing Asn270 in PON1) models were designed. The results indicated that the substitution of a water molecule for Asn270 in PON1 had little effect on the enzyme activity in kinetics, while being more efficient in thermodynamics, which is essential for DFP hydrolysis. Structure comparisons of evolutionarily related enzymes show that the mutation of Asn270 leads to the catalytic Ca2+ ion indirectly connecting the buried structural Ca2+ ion via hydrogen bonds in DFPase. It can reduce the plasticity of enzymatic structure, and possibly change the substrate preference from paraoxon to DFP, which implies an evolutionary transition from mono- to dinuclear catalytic centers. Our studies shed light on the investigation of enzyme catalysis mechanism from an evolutionary perspective.
Collapse
|
5
|
What roles do the residue Asp229 and the coordination variation of calcium play of the reaction mechanism of the diisopropyl-fluorophosphatase? A DFT investigation. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1896-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
6
|
Liao RZ, Siegbahn PEM. Phosphate Hydrolysis by the Fe2–Ca3-Dependent Alkaline Phosphatase PhoX: Mechanistic Insights from DFT calculations. Inorg Chem 2015; 54:11941-7. [DOI: 10.1021/acs.inorgchem.5b02268] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rong-Zhen Liao
- Key Laboratory of Material Chemistry for
Energy Conversion and Storage, Ministry of Education, School of Chemistry
and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Per E. M. Siegbahn
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
7
|
Valdez CE, Smith QA, Nechay MR, Alexandrova AN. Mysteries of metals in metalloenzymes. Acc Chem Res 2014; 47:3110-7. [PMID: 25207938 DOI: 10.1021/ar500227u] [Citation(s) in RCA: 101] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn(2+) in the active site but becomes much more active with Fe(2+). For β-lactamases, the replacement of the native Zn(2+) with Ni(2+) was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe(2+)- and Mn(2+)-dependent ribonucleotide reductases and W(4+)- and Mo(4+)-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe(2+)) by creating protein structures that utilize another metal (e.g., Mn(2+)) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.
Collapse
Affiliation(s)
- Crystal E. Valdez
- Department
of Chemistry and Biochemistry, and ‡California NanoSystems Institute, University of California, Los Angeles, California 90095, United States
| | - Quentin A. Smith
- Department
of Chemistry and Biochemistry, and ‡California NanoSystems Institute, University of California, Los Angeles, California 90095, United States
| | - Michael R. Nechay
- Department
of Chemistry and Biochemistry, and ‡California NanoSystems Institute, University of California, Los Angeles, California 90095, United States
| | - Anastassia N. Alexandrova
- Department
of Chemistry and Biochemistry, and ‡California NanoSystems Institute, University of California, Los Angeles, California 90095, United States
| |
Collapse
|