1
|
Zhou S, Wang L. Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects. Phys Chem Chem Phys 2020; 22:4884-4895. [DOI: 10.1039/c9cp06840f] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Electronic and nuclear quantum effects determine the symmetry and highly downfield 1H NMR chemical shifts of short hydrogen bonds.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology
- Institute for Quantitative Biomedicine
- Rutgers University
- Piscataway
- USA
| | - Lu Wang
- Department of Chemistry and Chemical Biology
- Institute for Quantitative Biomedicine
- Rutgers University
- Piscataway
- USA
| |
Collapse
|
2
|
Welborn VV, Head-Gordon T. Fluctuations of Electric Fields in the Active Site of the Enzyme Ketosteroid Isomerase. J Am Chem Soc 2019; 141:12487-12492. [DOI: 10.1021/jacs.9b05323] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Valerie Vaissier Welborn
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
3
|
Zhou S, Wang L. Unraveling the structural and chemical features of biological short hydrogen bonds. Chem Sci 2019; 10:7734-7745. [PMID: 31588321 PMCID: PMC6764281 DOI: 10.1039/c9sc01496a] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Accepted: 06/30/2019] [Indexed: 02/06/2023] Open
Abstract
Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.
The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| | - Lu Wang
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| |
Collapse
|
4
|
Afonin AV, Pavlov DV, Vashchenko AV. Case study of 2-vinyloxypyridine: Quantitative assessment of the intramolecular C H⋯N hydrogen bond energy and its contribution to the one-bond 13C1H coupling constant. J Mol Struct 2019. [DOI: 10.1016/j.molstruc.2018.08.050] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
5
|
Zhou F, Liu Y, Wang Z, Lu T, Yang Q, Liu Y, Zheng B. A new type of halogen bond involving multivalent astatine: an ab initio study. Phys Chem Chem Phys 2019; 21:15310-15318. [PMID: 31241070 DOI: 10.1039/c9cp02406a] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Theoretical studies on the dimers formed by CO with the halides of multivalent astatine as a Lewis-acid center are carried out to examine the typical characteristics of supervalent halogen bonds. Calculations at the MP2/aug-cc-pVTZ level reveal that the multiple nucleophilic sites of multivalent halide monomers can promote the formation of various types of halogen bonds, among which the most stable ones are At-halogen bond complexes with multivalent astatine as a Lewis acid center, followed by the π-halogen bond dimers, and the weakest ones are the X-halogen bonds. Compared with multivalent Cl-, Br-, and I-centers, At, as the heaviest halogen, exhibits the highest halogen-bond donating ability. We found that the electrostatic term and the dispersion term play an important role in the overall attractive interaction energy, and the smallest attraction term for all complexes is the polarization term (ΔEpol). Moreover, the tri and pentavalent halides analyzed here possess very "flexible" tautomerism in which the transformation occurs during the formation of the dimers. AIM theory and NBO analysis are also employed here.
Collapse
Affiliation(s)
- Fengxiang Zhou
- Key Laboratory of Theoretical Organic Chemistry and Functional Molecular, Ministry of Education, Hunan University of Science and Technology, Xiangtan, 411201, China. and Laboratory of Controllable Preparation and Functional Application of Fine Polymers, School of Chemistry and Chemical Engineering, Hunan University of Science and Technology, Xiangtan, 411201, China and Department of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, China.
| | - Yuan Liu
- Key Laboratory of Theoretical Organic Chemistry and Functional Molecular, Ministry of Education, Hunan University of Science and Technology, Xiangtan, 411201, China. and Laboratory of Controllable Preparation and Functional Application of Fine Polymers, School of Chemistry and Chemical Engineering, Hunan University of Science and Technology, Xiangtan, 411201, China
| | - Zhaoxu Wang
- Key Laboratory of Theoretical Organic Chemistry and Functional Molecular, Ministry of Education, Hunan University of Science and Technology, Xiangtan, 411201, China. and Laboratory of Controllable Preparation and Functional Application of Fine Polymers, School of Chemistry and Chemical Engineering, Hunan University of Science and Technology, Xiangtan, 411201, China
| | - Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing 100022, China
| | - Qingyuan Yang
- Department of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, China.
| | - Yi Liu
- Key Laboratory of Theoretical Organic Chemistry and Functional Molecular, Ministry of Education, Hunan University of Science and Technology, Xiangtan, 411201, China.
| | - Baishu Zheng
- Key Laboratory of Theoretical Organic Chemistry and Functional Molecular, Ministry of Education, Hunan University of Science and Technology, Xiangtan, 411201, China.
| |
Collapse
|
6
|
Li C, Guo W, Zhou P, Tang Z. The effects of the heteroatom and position on excited-state intramolecular proton transfer of new hydroxyphenyl benzoxazole derivatives: a time-dependent density functional theory study. Org Chem Front 2019. [DOI: 10.1039/c9qo00295b] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The effects of the heteroatom and position on excited-state intramolecular proton transfer (ESIPT) of 2-[4′-(N-4,6-dichloro-1,3,5-triazi-n-2-yl)2′hydroxyphenyl]benzoxazole (4THBO) have been investigated via time-dependent density functional theory studies.
Collapse
Affiliation(s)
- Changming Li
- School of Electrical Engineering
- University of South China
- Hengyang 421001
- China
- State Key Lab of Molecular Reaction Dynamics
| | - Wei Guo
- School of Electrical Engineering
- University of South China
- Hengyang 421001
- China
| | - Panwang Zhou
- Institute of Molecular Sciences and Engineering
- Shandong University
- Qingdao 266235
- P. R. China
| | - Zhe Tang
- Institute of Molecular Sciences and Engineering
- Shandong University
- Qingdao 266235
- P. R. China
| |
Collapse
|
7
|
Kim JJ, Lim CR, Reddy BM, Park SE. Hierarchical porous organic polymer as an efficient metal-free catalyst for acetalization of carbonyl compounds with alcohols. MOLECULAR CATALYSIS 2018. [DOI: 10.1016/j.mcat.2017.10.022] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
8
|
Hao J, Yang Yang YY. The effects of different heterocycles and solvents on the ESIPT mechanisms of three novel photoactive mono-formylated benzoxazole derivatives. Org Chem Front 2018. [DOI: 10.1039/c8qo00349a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The detailed effects of different heterocycles and solvents on the dynamical ESIPT mechanisms of three novel mono-formylated benzoxazole derivatives A–C in two different surroundings have been expounded by the TDDFT method at the B3LYP/TZVP level.
Collapse
Affiliation(s)
- Jiaojiao Hao
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Yang Yang Yang Yang
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| |
Collapse
|
9
|
Hao J, Yang Y. Exploring the ESIPT dynamical processes of two novel chromophores: symmetrical structure CHC and asymmetric structure CHN. Org Chem Front 2018. [DOI: 10.1039/c7qo01127j] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The detailed ESIPT dynamical processes of CHC (symmetrical structure) and CHN (asymmetric structure) chromophores were revealed and compared using the TDDFT method at the B3LYP/6-31+G(d,p) level.
Collapse
Affiliation(s)
- Jiaojiao Hao
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Yang Yang
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| |
Collapse
|
10
|
Wang L, Fried SD, Markland TE. Proton Network Flexibility Enables Robustness and Large Electric Fields in the Ketosteroid Isomerase Active Site. J Phys Chem B 2017; 121:9807-9815. [DOI: 10.1021/acs.jpcb.7b06985] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Lu Wang
- Department
of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Stephen D. Fried
- Medical Research
Council Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge CB2 0QH, U.K
| | - Thomas E. Markland
- Department
of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
11
|
Hanaya K, Yoshioka S, Ariyasu S, Aoki S, Shoji M, Sugai T. Development of a novel sulfonate ester-based prodrug strategy. Bioorg Med Chem Lett 2016; 26:545-550. [DOI: 10.1016/j.bmcl.2015.11.074] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Revised: 11/14/2015] [Accepted: 11/20/2015] [Indexed: 11/16/2022]
|
12
|
Samet M, Kass SR. Preorganized Hydrogen Bond Donor Catalysts: Acidities and Reactivities. J Org Chem 2015; 80:7727-31. [DOI: 10.1021/acs.joc.5b01475] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Masoud Samet
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Steven R. Kass
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
13
|
Brodkin HR, DeLateur NA, Somarowthu S, Mills CL, Novak WR, Beuning PJ, Ringe D, Ondrechen MJ. Prediction of distal residue participation in enzyme catalysis. Protein Sci 2015; 24:762-78. [PMID: 25627867 PMCID: PMC4420525 DOI: 10.1002/pro.2648] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2014] [Revised: 01/10/2015] [Accepted: 01/26/2015] [Indexed: 11/09/2022]
Abstract
A scoring method for the prediction of catalytically important residues in enzyme structures is presented and used to examine the participation of distal residues in enzyme catalysis. Scores are based on the Partial Order Optimum Likelihood (POOL) machine learning method, using computed electrostatic properties, surface geometric features, and information obtained from the phylogenetic tree as input features. Predictions of distal residue participation in catalysis are compared with experimental kinetics data from the literature on variants of the featured enzymes; some additional kinetics measurements are reported for variants of Pseudomonas putida nitrile hydratase (ppNH) and for Escherichia coli alkaline phosphatase (AP). The multilayer active sites of P. putida nitrile hydratase and of human phosphoglucose isomerase are predicted by the POOL log ZP scores, as is the single-layer active site of P. putida ketosteroid isomerase. The log ZP score cutoff utilized here results in over-prediction of distal residue involvement in E. coli alkaline phosphatase. While fewer experimental data points are available for P. putida mandelate racemase and for human carbonic anhydrase II, the POOL log ZP scores properly predict the previously reported participation of distal residues.
Collapse
Affiliation(s)
- Heather R Brodkin
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
- Department of Biochemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
- Department of Chemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
| | - Nicholas A DeLateur
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
| | - Srinivas Somarowthu
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
| | - Caitlyn L Mills
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
| | - Walter R Novak
- Department of Biochemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
- Department of Chemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
| | - Penny J Beuning
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
| | - Dagmar Ringe
- Department of Biochemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
- Department of Chemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis UniversityWaltham, Massachusetts, 02454–9110
| | - Mary Jo Ondrechen
- Department of Chemistry and Chemical Biology, Northeastern UniversityBoston, Massachusetts, 02115
| |
Collapse
|
14
|
Sigala PA, Ruben EA, Liu CW, Piccoli PMB, Hohenstein EG, Martínez TJ, Schultz AJ, Herschlag D. Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability. J Am Chem Soc 2015; 137:5730-40. [PMID: 25871450 DOI: 10.1021/ja512980h] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.
Collapse
Affiliation(s)
| | | | | | - Paula M B Piccoli
- §Intense Pulsed Neutron Source, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | | | | | - Arthur J Schultz
- §Intense Pulsed Neutron Source, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | | |
Collapse
|
15
|
Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site. Proc Natl Acad Sci U S A 2014; 111:18454-9. [PMID: 25503367 DOI: 10.1073/pnas.1417923111] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.
Collapse
|
16
|
Platzer G, Okon M, McIntosh LP. pH-dependent random coil (1)H, (13)C, and (15)N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements. JOURNAL OF BIOMOLECULAR NMR 2014; 60:109-129. [PMID: 25239571 DOI: 10.1007/s10858-014-9862-y] [Citation(s) in RCA: 149] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Accepted: 09/09/2014] [Indexed: 06/03/2023]
Abstract
The pK a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the (1)H, (13)C, and (15)N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain (15)N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.
Collapse
Affiliation(s)
- Gerald Platzer
- Department of Biochemistry and Molecular Biology, Life Sciences Centre, 2350 Health Sciences Mall, University of British Columbia, Vancouver, BC, V6T 1Z3, Canada
| | | | | |
Collapse
|
17
|
Samet M, Wang XB, Kass SR. A Preorganized Hydrogen Bond Network and Its Effect on Anion Stability. J Phys Chem A 2014; 118:5989-93. [DOI: 10.1021/jp505308v] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Masoud Samet
- Department
of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Xue-Bin Wang
- Physical
Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Steven R. Kass
- Department
of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
18
|
Natarajan A, Schwans JP, Herschlag D. Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site. J Am Chem Soc 2014; 136:7643-54. [PMID: 24787954 PMCID: PMC4046884 DOI: 10.1021/ja413174b] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2013] [Indexed: 02/04/2023]
Abstract
Hydrogen bonds are ubiquitous in enzyme active sites, providing binding interactions and stabilizing charge rearrangements on substrate groups over the course of a reaction. But understanding the origin and magnitude of their catalytic contributions relative to hydrogen bonds made in aqueous solution remains difficult, in part because of complexities encountered in energetic interpretation of traditional site-directed mutagenesis experiments. It has been proposed for ketosteroid isomerase and other enzymes that active site hydrogen bonding groups provide energetic stabilization via "short, strong" or "low-barrier" hydrogen bonds that are formed due to matching of their pKa or proton affinity to that of the transition state. It has also been proposed that the ketosteroid isomerase and other enzyme active sites provide electrostatic environments that result in larger energetic responses (i.e., greater "sensitivity") to ground-state to transition-state charge rearrangement, relative to aqueous solution, thereby providing catalysis relative to the corresponding reaction in water. To test these models, we substituted tyrosine with fluorotyrosines (F-Tyr's) in the ketosteroid isomerase (KSI) oxyanion hole to systematically vary the proton affinity of an active site hydrogen bond donor while minimizing steric or structural effects. We found that a 40-fold increase in intrinsic F-Tyr acidity caused no significant change in activity for reactions with three different substrates. F-Tyr substitution did not change the solvent or primary kinetic isotope effect for proton abstraction, consistent with no change in mechanism arising from these substitutions. The observed shallow dependence of activity on the pKa of the substituted Tyr residues suggests that the KSI oxyanion hole does not provide catalysis by forming an energetically exceptional pKa-matched hydrogen bond. In addition, the shallow dependence provides no indication of an active site electrostatic environment that greatly enhances the energetic response to charge accumulation, consistent with prior experimental results.
Collapse
Affiliation(s)
- Aditya Natarajan
- Department of Biochemistry, Stanford University School of Medicine, Stanford, California 94305, United States
| | | | - Daniel Herschlag
- Department of Biochemistry, Stanford University School of Medicine, Stanford, California 94305, United States
| |
Collapse
|
19
|
Schwans JP, Hanoian P, Lengerich BJ, Sunden F, Gonzalez A, Tsai Y, Hammes-Schiffer S, Herschlag D. Experimental and computational mutagenesis to investigate the positioning of a general base within an enzyme active site. Biochemistry 2014; 53:2541-55. [PMID: 24597914 PMCID: PMC4004248 DOI: 10.1021/bi401671t] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The
positioning of catalytic groups within proteins plays an important
role in enzyme catalysis, and here we investigate the positioning
of the general base in the enzyme ketosteroid isomerase (KSI). The
oxygen atoms of Asp38, the general base in KSI, were previously shown
to be involved in anion–aromatic interactions with two neighboring
Phe residues. Here we ask whether those interactions are sufficient,
within the overall protein architecture, to position Asp38 for catalysis
or whether the side chains that pack against Asp38 and/or the residues
of the structured loop that is capped by Asp38 are necessary to achieve
optimal positioning for catalysis. To test positioning, we mutated
each of the aforementioned residues, alone and in combinations, in
a background with the native Asp general base and in a D38E mutant
background, as Glu at position 38 was previously shown to be mispositioned
for general base catalysis. These double-mutant cycles reveal positioning
effects as large as 103-fold, indicating that structural
features in addition to the overall protein architecture and the Phe
residues neighboring the carboxylate oxygen atoms play roles in positioning.
X-ray crystallography and molecular dynamics simulations suggest that
the functional effects arise from both restricting dynamic fluctuations
and disfavoring potential mispositioned states. Whereas it may have
been anticipated that multiple interactions would be necessary for
optimal general base positioning, the energetic contributions from
positioning and the nonadditive nature of these interactions are not
revealed by structural inspection and require functional dissection.
Recognizing the extent, type, and energetic interconnectivity of interactions
that contribute to positioning catalytic groups has implications for
enzyme evolution and may help reveal the nature and extent of interactions
required to design enzymes that rival those found in biology.
Collapse
Affiliation(s)
- Jason P Schwans
- Department of Biochemistry, Stanford University , B400 Beckman Center, 279 Campus Drive, Stanford, California 94305, United States
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Gu J, Liu M, Guo F, Xie W, Lu W, Ye L, Chen Z, Yuan S, Yu H. Virtual screening of mandelate racemase mutants with enhanced activity based on binding energy in the transition state. Enzyme Microb Technol 2014; 55:121-7. [DOI: 10.1016/j.enzmictec.2013.10.008] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2013] [Revised: 10/22/2013] [Accepted: 10/23/2013] [Indexed: 11/30/2022]
|
21
|
Wanninayake U, Walker KD. A Bacterial Tyrosine Aminomutase Proceeds through Retention or Inversion of Stereochemistry To Catalyze Its Isomerization Reaction. J Am Chem Soc 2013; 135:11193-204. [DOI: 10.1021/ja403918w] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Udayanga Wanninayake
- Department
of Chemistry, and ‡Department of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| | - Kevin D. Walker
- Department
of Chemistry, and ‡Department of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
22
|
Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site. Proc Natl Acad Sci U S A 2013; 110:E2552-61. [PMID: 23798390 DOI: 10.1073/pnas.1302191110] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.
Collapse
|
23
|
van der Kamp MW, Mulholland AJ. Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. Biochemistry 2013; 52:2708-28. [PMID: 23557014 DOI: 10.1021/bi400215w] [Citation(s) in RCA: 399] [Impact Index Per Article: 36.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.
Collapse
Affiliation(s)
- Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| | | |
Collapse
|
24
|
Layfield JP, Hammes-Schiffer S. Calculation of vibrational shifts of nitrile probes in the active site of ketosteroid isomerase upon ligand binding. J Am Chem Soc 2012; 135:717-25. [PMID: 23210919 DOI: 10.1021/ja3084384] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The vibrational Stark effect provides insight into the roles of hydrogen bonding, electrostatics, and conformational motions in enzyme catalysis. In a recent application of this approach to the enzyme ketosteroid isomerase (KSI), thiocyanate probes were introduced in site-specific positions throughout the active site. This paper implements a quantum mechanical/molecular mechanical (QM/MM) approach for calculating the vibrational shifts of nitrile (CN) probes in proteins. This methodology is shown to reproduce the experimentally measured vibrational shifts upon binding of the intermediate analogue equilinen to KSI for two different nitrile probe positions. Analysis of the molecular dynamics simulations provides atomistic insight into the roles that key residues play in determining the electrostatic environment and hydrogen-bonding interactions experienced by the nitrile probe. For the M116C-CN probe, equilinen binding reorients an active-site water molecule that is directly hydrogen-bonded to the nitrile probe, resulting in a more linear C≡N--H angle and increasing the CN frequency upon binding. For the F86C-CN probe, equilinen binding orients the Asp103 residue, decreasing the hydrogen-bonding distance between the Asp103 backbone and the nitrile probe and slightly increasing the CN frequency. This QM/MM methodology is applicable to a wide range of biological systems and has the potential to assist in the elucidation of the fundamental principles underlying enzyme catalysis.
Collapse
Affiliation(s)
- Joshua P Layfield
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | | |
Collapse
|
25
|
Abstract
This brief review analyzes the underlying physical principles of enzyme catalysis, with an emphasis on the role of equilibrium enzyme motions and conformational sampling. The concepts are developed in the context of three representative systems, namely, dihydrofolate reductase, ketosteroid isomerase, and soybean lipoxygenase. All of these reactions involve hydrogen transfer, but many of the concepts discussed are more generally applicable. The factors that are analyzed in this review include hydrogen tunneling, proton donor-acceptor motion, hydrogen bonding, pKa shifting, electrostatics, preorganization, reorganization, and conformational motions. The rate constant for the chemical step is determined primarily by the free energy barrier, which is related to the probability of sampling configurations conducive to the chemical reaction. According to this perspective, stochastic thermal motions lead to equilibrium conformational changes in the enzyme and ligands that result in configurations favorable for the breaking and forming of chemical bonds. For proton, hydride, and proton-coupled electron transfer reactions, typically the donor and acceptor become closer to facilitate the transfer. The impact of mutations on the catalytic rate constants can be explained in terms of the factors enumerated above. In particular, distal mutations can alter the conformational motions of the enzyme and therefore the probability of sampling configurations conducive to the chemical reaction. Methods such as vibrational Stark spectroscopy, in which environmentally sensitive probes are introduced site-specifically into the enzyme, provide further insight into these aspects of enzyme catalysis through a combination of experiments and theoretical calculations.
Collapse
Affiliation(s)
- Sharon Hammes-Schiffer
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana-Champaign , Urbana, Illinois 61801, United States
| |
Collapse
|
26
|
Tian W, Huang X, Li Q, Li W, Cheng J, Gong B. Effect of superalkali substituents on the strengths and properties of hydrogen and halogen bonds. J Mol Model 2012. [DOI: 10.1007/s00894-012-1685-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
27
|
Hanoian P, Hammes-Schiffer S. Water in the active site of ketosteroid isomerase. Biochemistry 2011; 50:6689-700. [PMID: 21710970 DOI: 10.1021/bi200703y] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Classical molecular dynamics simulations were utilized to investigate the structural and dynamical properties of water in the active site of ketosteroid isomerase (KSI) to provide insight into the role of these water molecules in the enzyme-catalyzed reaction. This reaction is thought to proceed via a dienolate intermediate that is stabilized by hydrogen bonding with residues Tyr16 and Asp103. A comparative study was performed for the wild-type (WT) KSI and the Y16F, Y16S, and Y16F/Y32F/Y57F (FFF) mutants. These systems were studied with three different bound ligands: equilenin, which is an intermediate analog, and the intermediate states of two steroid substrates. Several distinct water occupation sites were identified in the active site of KSI for the WT and mutant systems. Three additional sites were identified in the Y16S mutant that were not occupied in WT KSI or the other mutants studied. The number of water molecules directly hydrogen bonded to the ligand oxygen was approximately two in the Y16S mutant and one in the Y16F and FFF mutants, with intermittent hydrogen bonding of one water molecule in WT KSI. The molecular dynamics trajectories of the Y16F and FFF mutants reproduced the small conformational changes of residue 16 observed in the crystal structures of these two mutants. Quantum mechanical/molecular mechanical calculations of (1)H NMR chemical shifts of the protons in the active site hydrogen-bonding network suggest that the presence of water in the active site does not prevent the formation of short hydrogen bonds with far-downfield chemical shifts. The molecular dynamics simulations indicate that the active site water molecules exchange much more frequently for WT KSI and the FFF mutant than for the Y16F and Y16S mutants. This difference is most likely due to the hydrogen-bonding interaction between Tyr57 and an active site water molecule that is persistent in the Y16F and Y16S mutants but absent in the FFF mutant and significantly less probable in WT KSI.
Collapse
Affiliation(s)
- Philip Hanoian
- Department of Chemistry, 104 Chemistry Building, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | |
Collapse
|
28
|
Chen M, Jiang M, Sun Y, Guo ZF, Guo Z. Stabilization of the second oxyanion intermediate by 1,4-dihydroxy-2-naphthoyl-coenzyme A synthase of the menaquinone pathway: spectroscopic evidence of the involvement of a conserved aspartic acid. Biochemistry 2011; 50:5893-904. [PMID: 21627110 DOI: 10.1021/bi200376x] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
1,4-Dihydroxy-2-naphthoyl-coenzyme A (DHNA-CoA) synthase, or MenB, catalyzes an intramolecular Claisen condensation involving two oxyanion intermediates in the biosynthetic pathway of menaquinone, an essential respiration electron transporter in many microorganisms. Here we report the finding that the DHNA-CoA product and its analogues bind and inhibit the synthase from Escherichia coli with significant ultraviolet--visible spectral changes, which are similar to the changes induced by deprotonation of the free inhibitors in a basic solution. Dissection of the structure--affinity relationships of the inhibitors identifies the hydroxyl groups at positions 1 (C1-OH) and 4 (C4-OH) of DHNA-CoA or their equivalents as the dominant and minor sites, respectively, for the enzyme--ligand interaction that polarizes or deprotonates the bound ligands to cause the observed spectral changes. In the meantime, spectroscopic studies with active site mutants indicate that C4-OH of the enzyme-bound DHNA-CoA interacts with conserved polar residues Arg-91, Tyr-97, and Tyr-258 likely through a hydrogen bonding network that also includes Ser-161. In addition, site-directed mutation of the conserved Asp-163 to alanine causes a complete loss of the ligand binding ability of the protein, suggesting that the Asp-163 side chain is most likely hydrogen-bonded to C1-OH of DHNA-CoA to provide the dominant polarizing effect. Moreover, this mutation also completely eliminates the enzyme activity, strongly supporting the possibility that the Asp-163 side chain provides a strong stabilizing hydrogen bond to the tetrahedral oxyanion, which takes a position similar to that of C1-OH of the enzyme-bound DHNA-CoA and is the second high-energy intermediate in the intracellular Claisen condensation reaction. Interestingly, both Arg-91 and Tyr-97 are located in a disordered loop forming part of the active site of all available DHNA-CoA synthase structures. Their involvement in the interaction with the small molecule ligands suggests that the disordered loop is folded in interaction with the substrates or reaction intermediates, supporting an induced-fit catalytic mechanism for the enzyme.
Collapse
Affiliation(s)
- Minjiao Chen
- Department of Chemistry and State Key Laboratory for Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong SAR, China
| | | | | | | | | |
Collapse
|