1
|
Hegazy R, Cristobal JR, Richard JP. Glycerol 3-Phosphate Dehydrogenase Catalyzed Hydride Transfer: Enzyme Activation by Cofactor Pieces. Biochemistry 2024. [PMID: 39319842 DOI: 10.1021/acs.biochem.4c00324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/26/2024]
Abstract
Glycerol 3-phosphate dehydrogenase catalyzes reversible hydride transfer from glycerol 3-phosphate (G3P) to NAD+ to form dihydroxyacetone phosphate; from the truncated substrate ethylene glycol to NAD+ in a reaction activated by the phosphite dianion substrate fragment; and from G3P to the truncated nicotinamide riboside cofactor in a reaction activated by adenosine 5'-diphosphate, adenosine 5'-monophosphate, and ribose 5-phosphate cofactor fragments. The sum of the stabilization of the transition state for GPDH-catalyzed hydride transfer reactions of the whole substrates by the phosphodianion fragment of G3P and the ADP fragment of NAD+ is 25 kcal/mol. Fourteen kcal/mol of this transition state stabilization is recovered as phosphite dianion and AMP activation of the reactions of the substrate and cofactor fragments. X-ray crystal structures for unliganded GPDH, for a binary GPDH·NAD+ complex, and for a nonproductive ternary GPDH·NAD+·DHAP complex show that the ligand binding energy is utilized to drive an extensive protein conformational change that creates a caged complex for these ligands. The phosphite dianion and AMP fragments are proposed to activate GPDH for the catalysis of hydride transfer by stabilization of this active caged complex. The closure of a conserved loop [292-LNGQKL-297] during substrate binding stabilizes the G3P and NAD+ complexes by interactions, respectively, with the Q295 and K296 loop side chains. The appearance and apparent conservation of two side chains that interact with the hydride donor and acceptor to stabilize the active closed enzyme are proposed to represent a significant improvement in the catalytic performance of GPDH.
Collapse
Affiliation(s)
- Rania Hegazy
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Judith R Cristobal
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| |
Collapse
|
2
|
Hegazy R, Richard JP. Triosephosphate Isomerase: The Crippling Effect of the P168A/I172A Substitution at the Heart of an Enzyme Active Site. Biochemistry 2023; 62:2916-2927. [PMID: 37768194 PMCID: PMC10586322 DOI: 10.1021/acs.biochem.3c00414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2023] [Revised: 09/01/2023] [Indexed: 09/29/2023]
Abstract
The P168 and I172 side chains sit at the heart of the active site of triosephosphate isomerase (TIM) and play important roles in the catalysis of the isomerization reaction. The phosphodianion of substrate glyceraldehyde 3-phosphate (GAP) drives a conformational change at the TIM that creates a steric interaction with the P168 side chain that is relieved by the movement of P168 that carries the basic E167 side chain into a clamp that consists of the hydrophobic I172 and L232 side chains. The P168A/I172A substitution at TIM from Trypanosoma brucei brucei (TbbTIM) causes a large 120,000-fold decrease in kcat for isomerization of GAP that eliminates most of the difference in the reactivity of TIM compared to the small amine base quinuclidinone for deprotonation of catalyst-bound GAP. The I172A substitution causes a > 2-unit decrease in the pKa of the E167 carboxylic acid in a complex to the intermediate analog PGA, but the P168A substitution at the I172A variant has no further effect on this pKa. The P168A/I172A substitutions cause a 5-fold decrease in Km for the isomerization of GAP from a 0.9 kcal/mol stabilization of the substrate Michaelis complexes. The results show that the P168 and I172 side chains play a dual role in destabilizing the ground-state Michaelis complex to GAP and in promoting stabilization of the transition state for substrate isomerization. This is consistent with an important role for these side chains in an induced fit reaction mechanism [Richard, J. P. (2022) Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution. Biochemistry 61, 1533-1542].
Collapse
Affiliation(s)
- Rania Hegazy
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| | - John P. Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| |
Collapse
|
3
|
Cristobal J, Nagorski RW, Richard JP. Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces. Biochemistry 2023; 62:2314-2324. [PMID: 37463347 PMCID: PMC10399567 DOI: 10.1021/acs.biochem.3c00290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 06/29/2023] [Indexed: 07/20/2023]
Abstract
The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.
Collapse
Affiliation(s)
- Judith
R. Cristobal
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| | - Richard W. Nagorski
- Department
of Chemistry, Illinois State University, Normal, Illinois 61790-4160, United
States
| | - John P. Richard
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| |
Collapse
|
4
|
Hegazy R, Cordara G, Wierenga RK, Richard JP. The Role of Asn11 in Catalysis by Triosephosphate Isomerase. Biochemistry 2023; 62:1794-1806. [PMID: 37162263 PMCID: PMC10249627 DOI: 10.1021/acs.biochem.3c00133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2023] [Revised: 04/14/2023] [Indexed: 05/11/2023]
Abstract
Four catalytic amino acids at triosephosphate isomerase (TIM) are highly conserved: N11, K13, H95, and E167. Asparagine 11 is the last of these to be characterized in mutagenesis studies. The ND2 side chain atom of N11 is hydrogen bonded to the O-1 hydroxyl of enzyme-bound dihydroxyacetone phosphate (DHAP), and it sits in an extended chain of hydrogen-bonded side chains that includes T75' from the second subunit. The N11A variants of wild-type TIM from Trypanosoma brucei brucei (TbbTIM) and Leishmania mexicana (LmTIM) undergo dissociation from the dimer to monomer under our assay conditions. Values of Kas = 8 × 103 and 1 × 106 M-1, respectively, were determined for the conversion of monomeric N11A TbbTIM and LmTIM into their homodimers. The N11A substitution at the variant of LmTIM previously stabilized by the E65Q substitution gives the N11A/E65Q variant that is stable to dissociation under our assay conditions. The X-ray crystal structure of N11A/E65Q LmTIM shows an active site that is essentially superimposable on that for wild-type TbbTIM, which also has a glutamine at position 65. A comparison of the kinetic parameters for E65Q LmTIM and N11A/E65Q LmTIM-catalyzed reactions of (R)-glyceraldehyde 3-phosphate (GAP) and (DHAP) shows that the N11A substitution results in a (13-14)-fold decrease in kcat/Km for substrate isomerization and a similar decrease in kcat for DHAP but only a 2-fold decrease in kcat for GAP.
Collapse
Affiliation(s)
- Rania Hegazy
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| | - Gabriele Cordara
- Biocenter
Oulu, University of Oulu, P.O. Box 5000, FIN-90014 Oulu, Finland
| | - Rik K. Wierenga
- Faculty of
Biochemistry and Molecular Medicine, University
of Oulu, P.O. Box 5400, FIN-90014 Oulu, Finland
| | - John P. Richard
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| |
Collapse
|
5
|
Cristobal JR, Richard JP. Kinetics and mechanism for enzyme-catalyzed reactions of substrate pieces. Methods Enzymol 2023; 685:95-126. [PMID: 37245916 PMCID: PMC10251411 DOI: 10.1016/bs.mie.2023.03.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
The most important difference between enzyme and small molecule catalysts is that only enzymes utilize the large intrinsic binding energies of nonreacting portions of the substrate in stabilization of the transition state for the catalyzed reaction. A general protocol is described to determine the intrinsic phosphodianion binding energy for enzymatic catalysis of reactions of phosphate monoester substrates, and the intrinsic phosphite dianion binding energy in activation of enzymes for catalysis of phosphodianion truncated substrates, from the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates. The enzyme-catalyzed reactions so-far documented that utilize dianion binding interactions for enzyme activation; and, their phosphodianion truncated substrates are summarized. A model for the utilization of dianion binding interactions for enzyme activation is described. The methods for the determination of the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates, from initial velocity data, are described and illustrated by graphical plots of kinetic data. The results of studies on the effect of site-directed amino acid substitutions at orotidine 5'-monophosphate decarboxylase, triosephosphate isomerase, and glycerol-3-phosphate dehydrogenase provide strong support for the proposal that these enzymes utilize binding interactions with the substrate phosphodianion to hold the protein catalysts in reactive closed conformations.
Collapse
Affiliation(s)
- Judith R Cristobal
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, NY, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, NY, United States.
| |
Collapse
|
6
|
Fernandez P, Richard JP. Adenylate Kinase-Catalyzed Reactions of AMP in Pieces: Specificity for Catalysis at the Nucleoside Activator and Dianion Catalytic Sites. Biochemistry 2022; 61:2766-2775. [PMID: 36413937 PMCID: PMC9731266 DOI: 10.1021/acs.biochem.2c00531] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Revised: 10/31/2022] [Indexed: 11/23/2022]
Abstract
The pressure to optimize the enzymatic rate acceleration for adenylate kinase (AK)-catalyzed phosphoryl transfer has led to the evolution of an induced-fit mechanism, where the binding energy from interactions between the protein and substrate adenosyl group is utilized to drive a protein conformational change that activates the enzyme for catalysis. The adenine group of adenosine contributes 11.8 kcal mol-1 to the total ≥14.7 kcal mol-1 adenosine stabilization of the transition state for AK-catalyzed phosphoryl transfer to AMP. The relative third-order rate constants for activation of adenylate kinase, by the C-5 truncated adenosine 1-(β-d-erythrofuranosyl)adenine (EA), for catalysis of phosphoryl transfer from ATP to phosphite dianion (HP, kcat/KHPKAct = 260 M-2 s-1), fluorophosphate (47 M-2 s-1), and phosphate (9.6 M-2 s-1), show that substitution of -F for -H and of -OH for -H at HP results, respectively, in decreases in the reactivity of AK for catalysis of phosphoryl transfer due to polar and steric effects of the -F and -OH substituents. The addition of a 5'-CH2OH to the EA activator results in a 3.0 kcal mol-1 destabilization of the transition state for AK-activated phosphoryl transfer to HP due to a steric effect. This is smaller than the 8.3 kcal mol-1 steric effect of the 5'-CH2OH substituent at OMP on HP-activated OMPDC-catalyzed decarboxylation of 1-(β-d-erythrofuranosyl)orotate. The 2'-OH ribosyl substituent shows significant interactions with the transition states for AK-catalyzed phosphoryl transfer from ATP to AMP and for adenosine-activated AK-catalyzed phosphoryl transfer from ATP to HP.
Collapse
Affiliation(s)
- Patrick
L. Fernandez
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York14260−3000, United States
| | - John P. Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York14260−3000, United States
| |
Collapse
|
7
|
Identification and Mutation Analysis of Nonconserved Residues on the TIM-Barrel Surface of GH5_5 Cellulases for Catalytic Efficiency and Stability Improvement. Appl Environ Microbiol 2022; 88:e0104622. [PMID: 36000858 PMCID: PMC9469711 DOI: 10.1128/aem.01046-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Exploring the potential functions of nonconserved residues on the outer side of α-helices and systematically optimizing them are pivotal for their application in protein engineering. Based on the evolutionary structural conservation analysis of GH5_5 cellulases, a practical molecular improvement strategy was developed. Highly variable sites on the outer side of the α-helices of the GH5_5 cellulase from Aspergillus niger (AnCel5A) were screened, and 14 out of the 34 highly variable sites were confirmed to exert a positive effect on the activity. After the modular combination of the positive mutations, the catalytic efficiency of the mutants was further improved. By using CMC-Na as the substrate, the catalytic efficiency and specific activity of variant AnCel5A_N193A/T300P/D307P were approximately 2.0-fold that of AnCel5A (227 ± 21 versus 451 ± 43 ml/s/mg and 1,726 ± 19 versus 3,472 ± 42 U/mg, respectively). The half-life (t1/2) of variant AnCel5A_N193A/T300P/D307P at 75°C was 2.36 times that of AnCel5A. The role of these sites was successfully validated in other GH5_5 cellulases. Computational analyses revealed that the flexibility of the loop 6-loop 7-loop 8 region was responsible for the increased catalytic performance. This work not only illustrated the important role of rapidly evolving positions on the outer side of the α-helices of GH5_5 cellulases but also revealed new insights into engineering the proteins that nature left as clues for us to find. IMPORTANCE A comprehensive understanding of the residues on the α-helices of the GH5_5 cellulases is important for catalytic efficiency and stability improvement. The main objective of this study was to use the evolutionary conservation and plasticity of the TIM-barrel fold to probe the relationship between nonconserved residues on the outer side of the α-helices and the catalytic efficiency of GH5_5 cellulases by conducting structure-guided protein engineering. By using a four-step nonconserved residue screening strategy, the functional role of nonconserved residues on the outer side of the α-helices was effectively identified, and a variant with superior performance and capability was constructed. Hence, this study proved the effectiveness of this strategy in engineering GH5_5 cellulases and provided a potential competitor for industrial applications. Furthermore, this study sheds new light on engineering TIM-barrel proteins.
Collapse
|
8
|
Abstract
Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.
Collapse
Affiliation(s)
- John P Richard
- Department of Chemistry, University at Buffalo, the State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
9
|
Cristobal JR, Brandão TAS, Reyes AC, Richard JP. Protein-Ribofuranosyl Interactions Activate Orotidine 5'-Monophosphate Decarboxylase for Catalysis. Biochemistry 2021; 60:3362-3373. [PMID: 34726391 DOI: 10.1021/acs.biochem.1c00589] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The role of a global, substrate-driven, enzyme conformational change in enabling the extraordinarily large rate acceleration for orotidine 5'-monophosphate decarboxylase (OMPDC)-catalyzed decarboxylation of orotidine 5'-monophosphate (OMP) is examined in experiments that focus on the interactions between OMPDC and the ribosyl hydroxyl groups of OMP. The D37 and T100' side chains of OMPDC interact, respectively, with the C-3' and C-2' hydroxyl groups of enzyme-bound OMP. D37G and T100'A substitutions result in 1.4 kcal/mol increases in the activation barrier ΔG⧧ for catalysis of decarboxylation of the phosphodianion-truncated substrate 1-(β-d-erythrofuranosyl)orotic acid (EO) but result in larger 2.1-2.9 kcal/mol increases in ΔG⧧ for decarboxylation of OMP and for phosphite dianion-activated decarboxylation of EO. This shows that these substitutions reduce transition-state stabilization by the Q215, Y217, and R235 side chains at the dianion binding site. The D37G and T100'A substitutions result in <1.0 kcal/mol increases in ΔG⧧ for activation of OMPDC-catalyzed decarboxylation of the phosphoribofuranosyl-truncated substrate FO by phosphite dianions. Experiments to probe the effect of D37 and T100' substitutions on the kinetic parameters for d-glycerol 3-phosphate and d-erythritol 4-phosphate activators of OMPDC-catalyzed decarboxylation of FO show that ΔG⧧ for sugar phosphate-activated reactions is increased by ca. 2.5 kcal/mol for each -OH interaction eliminated by D37G or T100'A substitutions. We conclude that the interactions between the D37 and T100' side chains and ribosyl or ribosyl-like hydroxyl groups are utilized to activate OMPDC for catalysis of decarboxylation of OMP, EO, and FO.
Collapse
Affiliation(s)
- Judith R Cristobal
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Tiago A S Brandão
- Department of Chemistry, ICEx, Federal University of Minas Gerais, Belo Horizonte, MG 31270-901, Brazil
| | - Archie C Reyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| |
Collapse
|
10
|
Kelpšas V, Caldararu O, Blakeley MP, Coquelle N, Wierenga RK, Ryde U, von Wachenfeldt C, Oksanen E. Neutron structures of Leishmania mexicana triosephosphate isomerase in complex with reaction-intermediate mimics shed light on the proton-shuttling steps. IUCRJ 2021; 8:633-643. [PMID: 34258011 PMCID: PMC8256706 DOI: 10.1107/s2052252521004619] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Accepted: 04/30/2021] [Indexed: 06/13/2023]
Abstract
Triosephosphate isomerase (TIM) is a key enzyme in glycolysis that catalyses the interconversion of glyceraldehyde 3-phosphate and dihydroxy-acetone phosphate. This simple reaction involves the shuttling of protons mediated by protolysable side chains. The catalytic power of TIM is thought to stem from its ability to facilitate the deprotonation of a carbon next to a carbonyl group to generate an enediolate intermediate. The enediolate intermediate is believed to be mimicked by the inhibitor 2-phosphoglycolate (PGA) and the subsequent enediol intermediate by phosphoglycolohydroxamate (PGH). Here, neutron structures of Leishmania mexicana TIM have been determined with both inhibitors, and joint neutron/X-ray refinement followed by quantum refinement has been performed. The structures show that in the PGA complex the postulated general base Glu167 is protonated, while in the PGH complex it remains deprotonated. The deuteron is clearly localized on Glu167 in the PGA-TIM structure, suggesting an asymmetric hydrogen bond instead of a low-barrier hydrogen bond. The full picture of the active-site protonation states allowed an investigation of the reaction mechanism using density-functional theory calculations.
Collapse
Affiliation(s)
- Vinardas Kelpšas
- Department of Biology, Lund University, Sölvegatan 35, 223 62 Lund, Sweden
| | - Octav Caldararu
- Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Matthew P. Blakeley
- Large-Scale Structures Group, Institut Laue–Langevin, 71 Avenue des Martyrs, 38042 Grenoble, France
| | - Nicolas Coquelle
- Large-Scale Structures Group, Institut Laue–Langevin, 71 Avenue des Martyrs, 38042 Grenoble, France
| | - Rikkert K. Wierenga
- Faculty of Biochemistry and Molecular Medicine, University of Oulu, Pentti Kaiteran katu 1, 90570 Oulu, Finland
| | - Ulf Ryde
- Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | | | - Esko Oksanen
- European Spallation Source Consortium ESS ERIC, Odarslövsvägen 113, 224 84 Lund, Sweden
| |
Collapse
|
11
|
Richard JP, Cristobal JR, Amyes TL. Linear Free Energy Relationships for Enzymatic Reactions: Fresh Insight from a Venerable Probe. Acc Chem Res 2021; 54:2532-2542. [PMID: 33939414 PMCID: PMC8157535 DOI: 10.1021/acs.accounts.1c00147] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
![]()
Linear free energy relationships (LFERs) for substituent effects on reactions that
proceed through similar transition states provide insight into transition state
structures. A classical approach to the analysis of LFERs showed that differences in the
slopes of Brønsted correlations for addition of substituted alkyl alcohols to
ring-substituted 1-phenylethyl carbocations and to the β-galactopyranosyl
carbocation intermediate of reactions catalyzed by β-galactosidase provide
evidence that the enzyme catalyst modifies the curvature of the energy surface at the
saddle point for the transition state for nucleophile addition. We have worked to
generalize the use of LFERs in the determination of enzyme mechanisms. The defining
property of enzyme catalysts is their specificity for binding the transition state with
a much higher affinity than the substrate. Triosephosphate isomerase (TIM), orotidine
5′-monophosphate decarboxylase (OMPDC), and glycerol 3-phosphate dehydrogenase
(GPDH) show effective catalysis of reactions of phosphorylated substrates and strong
phosphite dianion activation of reactions of phosphodianion truncated substrates, with
rate constants kcat/Km
(M–1 s–1) and
kcat/KdKHPi
(M–2 s–1), respectively. Good linear logarithmic
correlations, with a slope of 1.1, between these kinetic parameters determined for
reactions catalyzed by five or more variant forms of each catalyst are observed, where
the protein substitutions are mainly at side chains which function to stabilize the cage
complex between the enzyme and substrate. This shows that the enzyme-catalyzed reactions
of a whole substrate and substrate pieces proceed through transition states of similar
structures. It provides support for the proposal that the dianion binding energy of
whole phosphodianion substrates and of phosphite dianion is used to drive the conversion
of these protein catalysts from flexible and entropically rich ground states to stiff
and catalytically active Michaelis complexes that show the same activity toward
catalysis of the reactions of whole and phosphodianion truncated substrates. There is a
good linear correlation, with a slope of 0.73, between values of the dissociation
constants log Ki for release of the transition state analog
phosphoglycolate (PGA) trianion and log
kcat/Km for isomerization of
GAP for wild-type and variants of TIM. This correlation shows that the substituted amino
acid side chains act to stabilize the complex between TIM and the PGA trianion and that
ca. 70% of this stabilization is observed at the transition state for
substrate deprotonation. The correlation provides evidence that these side chains
function to enhance the basicity of the E165 side chain of TIM, which deprotonates the
bound carbon acid substrate. There is a good linear correlation, with a slope of 0.74,
between the values of ΔG‡ and
ΔG° determined by electron valence bond (EVB) calculations
to model deprotonation of dihydroxyacetone phosphate (DHAP) in water and when bound to
wild-type and variant forms of TIM to form the enediolate reaction intermediate. This
correlation provides evidence that the stabilizing interactions of the transition state
for TIM-catalyzed deprotonation of DHAP are optimized by placement of amino acid side
chains in positions that provide for the maximum stabilization of the charged reaction
intermediate, relative to the neutral substrate.
Collapse
Affiliation(s)
- John P. Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Judith R. Cristobal
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Tina L. Amyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| |
Collapse
|
12
|
Fernandez PL, Nagorski RW, Cristobal JR, Amyes TL, Richard JP. Phosphodianion Activation of Enzymes for Catalysis of Central Metabolic Reactions. J Am Chem Soc 2021; 143:2694-2698. [PMID: 33560827 PMCID: PMC7919737 DOI: 10.1021/jacs.0c13423] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
The activation barriers ΔG⧧ for
kcat/Km for the reactions of
whole substrates catalyzed by 6-phosphogluconate dehydrogenase, glucose 6-phosphate
dehydrogenase, and glucose 6-phosphate isomerase are reduced by 11–13 kcal/mol by
interactions between the protein and the substrate phosphodianion. Between 4 and 6
kcal/mol of this dianion binding energy is expressed at the transition state for
phosphite dianion activation of the respective enzyme-catalyzed reactions of truncated
substrates d-xylonate or d-xylose. These and earlier results from
studies on β-phosphoglucomutase, triosephosphate isomerase, and glycerol
3-phosphate dehydrogenase define a cluster of six enzymes that catalyze reactions in
glycolysis or of glycolytic intermediates, and which utilize substrate dianion binding
energy for enzyme activation. Dianion-driven conformational changes, which convert
flexible open proteins to tight protein cages for the phosphorylated substrate, have
been thoroughly documented for five of these six enzymes. The clustering of metabolic
enzymes which couple phosphodianion-driven conformational changes to enzyme activation
suggests that this catalytic motif has been widely propagated in the proteome.
Collapse
Affiliation(s)
- Patrick L Fernandez
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Richard W Nagorski
- Department of Chemistry, Illinois State University, Normal, Illinois 61790-4160, United States
| | - Judith R Cristobal
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Tina L Amyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| |
Collapse
|
13
|
Singh S, Kumar K, Nath P, Goyal A. Role of glycine 256 residue in improving the catalytic efficiency of mutant endoglucanase of family 5 glycoside hydrolase from Bacillus amyloliquefaciens SS35. Biotechnol Bioeng 2020; 117:2668-2682. [PMID: 32484905 DOI: 10.1002/bit.27448] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2020] [Revised: 05/24/2020] [Accepted: 05/31/2020] [Indexed: 11/07/2022]
Abstract
Wild-type, BaGH5-WT and mutant, BaGH5-UV2 (aspartate residue mutated to glycine), endoglucanases belonging to glycoside hydrolase family 5 (GH5), from wild-type, and UV2 mutant strain of Bacillus amyloliquefaciens SS35, respectively, were earlier cloned in pHTP0 cloning vector. In this study, genes encoding BaGH5-WT or BaGH5-UV2 were cloned into pET28a(+) expression-vector and expressed in Escherichia coli BL-21(DE3)pLysS cells. BaGH5-UV2 showed 10-fold (43.6 U/mg) higher specific activity against carboxymethylcellulose sodium salt (CMC-Na), higher optimal temperature by 10°C at 65°C, and 22-fold higher catalytic efficiency against CMC-Na, than BaGH5-WT. BaGH5-UV2 showed stability in wider acidic pH range (5.0-7.0) unlike BaGH5-WT in narrow basic pH range (7.0-7.5). BaGH5-UV2 displayed a mutation, Asp256Gly in L11 loop, connecting β6 -sheet with α6 -helix, near active site toward the domain surface of (α/β)8 -TIM barrel fold. Molecular dynamics simulation studies showed more stable structure, accessibility of substrate for a catalytic site, and increased flexibility of loop L11 of BaGH5-UV2 than the wild type, suggesting enhanced catalysis by BaGH5-UV2. Molecular docking analysis displayed enhanced hydrogen bond interactions of cello-oligosaccharides with BaGH5-UV2, unlike BaGH5-WT. Thus, Gly256 residue of loop L11 plays an important role in enhancing catalytic efficiency, and pH stability of GH5 endoglucanase. Therefore, these results help in protein engineering of GH5 endoglucanase for improved biochemical properties.
Collapse
Affiliation(s)
- Shweta Singh
- Carbohydrate Enzyme Biotechnology Laboratory, Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India.,DBT PAN-IIT Centre for Bioenergy, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| | - Krishan Kumar
- Carbohydrate Enzyme Biotechnology Laboratory, Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| | - Priyanka Nath
- Carbohydrate Enzyme Biotechnology Laboratory, Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India.,DBT PAN-IIT Centre for Bioenergy, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| | - Arun Goyal
- Carbohydrate Enzyme Biotechnology Laboratory, Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India.,DBT PAN-IIT Centre for Bioenergy, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| |
Collapse
|
14
|
The role of ligand-gated conformational changes in enzyme catalysis. Biochem Soc Trans 2020; 47:1449-1460. [PMID: 31657438 PMCID: PMC6824834 DOI: 10.1042/bst20190298] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Revised: 10/03/2019] [Accepted: 10/07/2019] [Indexed: 11/17/2022]
Abstract
Structural and biochemical studies on diverse enzymes have highlighted the importance of ligand-gated conformational changes in enzyme catalysis, where the intrinsic binding energy of the common phosphoryl group of their substrates is used to drive energetically unfavorable conformational changes in catalytic loops, from inactive open to catalytically competent closed conformations. However, computational studies have historically been unable to capture the activating role of these conformational changes. Here, we discuss recent experimental and computational studies, which can remarkably pinpoint the role of ligand-gated conformational changes in enzyme catalysis, even when not modeling the loop dynamics explicitly. Finally, through our joint analyses of these data, we demonstrate how the synergy between theory and experiment is crucial for furthering our understanding of enzyme catalysis.
Collapse
|
15
|
Kulkarni YS, Amyes TL, Richard JP, Kamerlin SCL. Uncovering the Role of Key Active-Site Side Chains in Catalysis: An Extended Brønsted Relationship for Substrate Deprotonation Catalyzed by Wild-Type and Variants of Triosephosphate Isomerase. J Am Chem Soc 2019; 141:16139-16150. [PMID: 31508957 PMCID: PMC7032883 DOI: 10.1021/jacs.9b08713] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
We report results of detailed empirical valence bond simulations that model the effect of several amino acid substitutions on the thermodynamic (ΔG°) and kinetic activation (ΔG⧧) barriers to deprotonation of dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) bound to wild-type triosephosphate isomerase (TIM), as well as to the K12G, E97A, E97D, E97Q, K12G/E97A, I170A, L230A, I170A/L230A, and P166A variants of this enzyme. The EVB simulations model the observed effect of the P166A mutation on protein structure. The E97A, E97Q, and E97D mutations of the conserved E97 side chain result in ≤1.0 kcal mol-1 decreases in the activation barrier for substrate deprotonation. The agreement between experimental and computed activation barriers is within ±1 kcal mol-1, with a strong linear correlation between ΔG⧧ and ΔG° for all 11 variants, with slopes β = 0.73 (R2 = 0.994) and β = 0.74 (R2 = 0.995) for the deprotonation of DHAP and GAP, respectively. These Brønsted-type correlations show that the amino acid side chains examined in this study function to reduce the standard-state Gibbs free energy of reaction for deprotonation of the weak α-carbonyl carbon acid substrate to form the enediolate phosphate reaction intermediate. TIM utilizes the cationic side chain of K12 to provide direct electrostatic stabilization of the enolate oxyanion, and the nonpolar side chains of P166, I170, and L230 are utilized for the construction of an active-site cavity that provides optimal stabilization of the enediolate phosphate intermediate relative to the carbon acid substrate.
Collapse
Affiliation(s)
- Yashraj S Kulkarni
- Science for Life Laboratory, Department of Chemistry - BMC , Uppsala University, BMC , Box 576, S-751 23 Uppsala , Sweden
| | - Tina L Amyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Shina C L Kamerlin
- Science for Life Laboratory, Department of Chemistry - BMC , Uppsala University, BMC , Box 576, S-751 23 Uppsala , Sweden
| |
Collapse
|
16
|
Roland BP, Richards KR, Hrizo SL, Eicher S, Barile ZJ, Chang TC, Savon G, Bianchi P, Fermo E, Ricerca BM, Tortorolo L, Vockley J, VanDemark AP, Palladino MJ. Missense variant in TPI1 (Arg189Gln) causes neurologic deficits through structural changes in the triosephosphate isomerase catalytic site and reduced enzyme levels in vivo. Biochim Biophys Acta Mol Basis Dis 2019; 1865:2257-2266. [PMID: 31075491 DOI: 10.1016/j.bbadis.2019.05.002] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Revised: 12/20/2018] [Accepted: 01/06/2019] [Indexed: 12/27/2022]
Abstract
Mutations in the gene triosephosphate isomerase (TPI) lead to a severe multisystem condition that is characterized by hemolytic anemia, a weakened immune system, and significant neurologic symptoms such as seizures, distal neuropathy, and intellectual disability. No effective therapy is available. Here we report a compound heterozygous patient with a novel TPI pathogenic variant (NM_000365.5:c.569G>A:p.(Arg189Gln)) in combination with the common (NM_000365.5:c.315G>C:p.(Glu104Asp)) allele. We characterized the novel variant by mutating the homologous Arg in Drosophila using a genomic engineering system, demonstrating that missense mutations at this position cause a strong loss of function. Compound heterozygote animals were generated and exhibit motor behavioural deficits and markedly reduced protein levels. Furthermore, examinations of the TPIArg189Gln/TPIGlu104Asp patient fibroblasts confirmed the reduction of TPI levels, suggesting that Arg189Gln may also affect the stability of the protein. The Arg189 residue participates in two salt bridges on the backside of the TPI enzyme dimer, and we reveal that a mutation at this position alters the coordination of the substrate-binding site and important catalytic residues. Collectively, these data reveal a new human pathogenic variant associated with TPI deficiency, identify the Arg189 salt bridge as critical for organizing the catalytic site of the TPI enzyme, and demonstrates that reduced TPI levels are associated with human TPI deficiency. These findings advance our understanding of the molecular pathogenesis of the disease, and suggest new therapeutic avenues for pre-clinical trials.
Collapse
Affiliation(s)
- Bartholomew P Roland
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Department of Pharmacology, Vanderbilt University, Nashville, TN 37232, USA
| | - Kristen R Richards
- Biological Sciences and Structural Biology, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | - Stacy L Hrizo
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Department of Biology, Slippery Rock University, Slippery Rock, PA 16057, USA
| | - Samantha Eicher
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Zackery J Barile
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Tien-Chien Chang
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Grace Savon
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Paola Bianchi
- Fondazione IRCCS Ca' Granda Ospedale Maggiore Policlinico, UOC Ematologia, UOS Fisiopatologia delle Anemie, Via F Sforza, 35, 20122 Milan, Italy
| | - Elisa Fermo
- Fondazione IRCCS Ca' Granda Ospedale Maggiore Policlinico, UOC Ematologia, UOS Fisiopatologia delle Anemie, Via F Sforza, 35, 20122 Milan, Italy
| | - Bianca Maria Ricerca
- Hematology Institute, Universitary Hospital A. Gemelli, Largo A. Gemelli 8, 00168 Rome, Italy
| | - Luca Tortorolo
- Pediatric Intensive Care Unit, Universitary Hospital A. Gemelli, Largo A. Gemelli 8, 00168 Rome, Italy
| | - Jerry Vockley
- Department of Pediatrics and Human Genetics, University of Pittsburgh Schools of Medicine and Public health, Pittsburgh, PA 15261, USA
| | - Andrew P VanDemark
- Biological Sciences and Structural Biology, University of Pittsburgh, Pittsburgh, PA 15260, USA.
| | - Michael J Palladino
- Department of Pharmacology & Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA; Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA.
| |
Collapse
|
17
|
Abstract
![]()
The enormous rate accelerations observed
for many enzyme catalysts
are due to strong stabilizing interactions between the protein and
reaction transition state. The defining property of these catalysts
is their specificity for binding the transition state with a much
higher affinity than substrate. Experimental results are presented
which show that the phosphodianion-binding energy of phosphate monoester
substrates is used to drive conversion of their protein catalysts
from flexible and entropically rich ground states to stiff and catalytically
active Michaelis complexes. These results are generalized to other
enzyme-catalyzed reactions. The existence of many enzymes in flexible,
entropically rich, and inactive ground states provides a mechanism
for utilization of ligand-binding energy to mold these catalysts into
stiff and active forms. This reduces the substrate-binding energy
expressed at the Michaelis complex, while enabling the full and specific
expression of large transition-state binding energies. Evidence is
presented that the complexity of enzyme conformational changes increases
with increases in the enzymatic rate acceleration. The requirement
that a large fraction of the total substrate-binding energy be utilized
to drive conformational changes of floppy enzymes is proposed to favor
the selection and evolution of protein folds with multiple flexible
unstructured loops, such as the TIM-barrel fold. The effect of protein
motions on the kinetic parameters for enzymes that undergo ligand-driven
conformational changes is considered. The results of computational
studies to model the complex ligand-driven conformational change in
catalysis by triosephosphate isomerase are presented.
Collapse
Affiliation(s)
- John P Richard
- Department of Chemistry , SUNY, University at Buffalo , Buffalo , New York 14260-3000 , United States
| |
Collapse
|
18
|
Mydy LS, Cristobal JR, Katigbak RD, Bauer P, Reyes AC, Kamerlin SCL, Richard JP, Gulick AM. Human Glycerol 3-Phosphate Dehydrogenase: X-ray Crystal Structures That Guide the Interpretation of Mutagenesis Studies. Biochemistry 2019; 58:1061-1073. [PMID: 30640445 DOI: 10.1021/acs.biochem.8b01103] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Human liver glycerol 3-phosphate dehydrogenase ( hlGPDH) catalyzes the reduction of dihydroxyacetone phosphate (DHAP) to form glycerol 3-phosphate, using the binding energy associated with the nonreacting phosphodianion of the substrate to properly orient the enzyme-substrate complex within the active site. Herein, we report the crystal structures for unliganded, binary E·NAD, and ternary E·NAD·DHAP complexes of wild type hlGPDH, illustrating a new position of DHAP, and probe the kinetics of multiple mutant enzymes with natural and truncated substrates. Mutation of Lys120, which is positioned to donate a proton to the carbonyl of DHAP, results in similar increases in the activation barrier to hlGPDH-catlyzed reduction of DHAP and to phosphite dianion-activated reduction of glycolaldehyde, illustrating that these transition states show similar interactions with the cationic K120 side chain. The K120A mutation results in a 5.3 kcal/mol transition state destabilization, and 3.0 kcal/mol of the lost transition state stabilization is rescued by 1.0 M ethylammonium cation. The 6.5 kcal/mol increase in the activation barrier observed for the D260G mutant hlGPDH-catalyzed reaction represents a 3.5 kcal/mol weakening of transition state stabilization by the K120A side chain and a 3.0 kcal/mol weakening of the interactions with other residues. The interactions, at the enzyme active site, between the K120 side chain and the Q295 and R269 side chains were likewise examined by double-mutant analyses. These results provide strong evidence that the enzyme rate acceleration is due mainly or exclusively to transition state stabilization by electrostatic interactions with polar amino acid side chains.
Collapse
Affiliation(s)
- Lisa S Mydy
- Department of Structural Biology , University at Buffalo, SUNY , Buffalo , New York 14203 , United States
| | - Judith R Cristobal
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Roberto D Katigbak
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Paul Bauer
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden
| | - Archie C Reyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Shina Caroline Lynn Kamerlin
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Andrew M Gulick
- Department of Structural Biology , University at Buffalo, SUNY , Buffalo , New York 14203 , United States
| |
Collapse
|
19
|
Castro-Torres E, Jimenez-Sandoval P, Fernández-de Gortari E, López-Castillo M, Baruch-Torres N, López-Hidalgo M, Peralta-Castro A, Díaz-Quezada C, Sotelo-Mundo RR, Benitez-Cardoza CG, Espinoza-Fonseca LM, Ochoa-Leyva A, Brieba LG. Structural Basis for the Limited Response to Oxidative and Thiol-Conjugating Agents by Triosephosphate Isomerase From the Photosynthetic Bacteria Synechocystis. Front Mol Biosci 2018; 5:103. [PMID: 30538993 PMCID: PMC6277545 DOI: 10.3389/fmolb.2018.00103] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2018] [Accepted: 11/05/2018] [Indexed: 11/18/2022] Open
Abstract
In plants, the ancestral cyanobacterial triosephosphate isomerase (TPI) was replaced by a duplicated version of the cytosolic TPI. This isoform acquired a transit peptide for chloroplast localization and functions in the Calvin-Benson cycle. To gain insight into the reasons for this gene replacement in plants, we characterized the TPI from the photosynthetic bacteria Synechocystis (SyTPI). SyTPI presents typical TPI enzyme kinetics profiles and assembles as a homodimer composed of two subunits that arrange in a (β-α)8 fold. We found that oxidizing agents diamide (DA) and H2O2, as well as thiol-conjugating agents such as oxidized glutathione (GSSG) and methyl methanethiosulfonate (MMTS), do not inhibit the catalytic activity of SyTPI at concentrations required to inactivate plastidic and cytosolic TPIs from the plant model Arabidopsis thaliana (AtpdTPI and AtcTPI, respectively). The crystal structure of SyTPI revealed that each monomer contains three cysteines, C47, C127, and C176; however only the thiol group of C176 is solvent exposed. While AtcTPI and AtpdTPI are redox-regulated by chemical modifications of their accessible and reactive cysteines, we found that C176 of SyTPI is not sensitive to redox modification in vitro. Our data let us postulate that SyTPI was replaced by a eukaryotic TPI, because the latter contains redox-sensitive cysteines that may be subject to post-translational modifications required for modulating TPI's enzymatic activity.
Collapse
Affiliation(s)
- Eduardo Castro-Torres
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Pedro Jimenez-Sandoval
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Eli Fernández-de Gortari
- Division of Cardiovascular Medicine, Department of Internal Medicine, Center for Arrhythmia Research, University of Michigan, Ann Arbor, MI, United States
| | - Margarita López-Castillo
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Noe Baruch-Torres
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Marisol López-Hidalgo
- Laboratorio de Investigación Bioquímica, Programa Institucional en Biomedicina Molecular ENMyH-IPN, Ciudad de México, Mexico
| | - Antolín Peralta-Castro
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Corina Díaz-Quezada
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| | - Rogerio R Sotelo-Mundo
- Laboratorio de Estructura Biomolecular, Centro de Investigación en Alimentación y Desarrollo, A.C., Hermosillo, Mexico
| | - Claudia G Benitez-Cardoza
- Laboratorio de Investigación Bioquímica, Programa Institucional en Biomedicina Molecular ENMyH-IPN, Ciudad de México, Mexico
| | - L Michel Espinoza-Fonseca
- Division of Cardiovascular Medicine, Department of Internal Medicine, Center for Arrhythmia Research, University of Michigan, Ann Arbor, MI, United States
| | - Adrian Ochoa-Leyva
- Departamento de Microbiología Molecular, Instituto de Biotecnología, Universidad Nacional Autónoma de México, Cuernavaca, Mexico
| | - Luis G Brieba
- Laboratorio Nacional de Genómica para la Biodiversidad, Centro de Investigación y de Estudios Avanzados del IPN, Guanajuato, Mexico
| |
Collapse
|
20
|
Reyes AC, Plache DC, Koudelka AP, Amyes TL, Gerlt JA, Richard JP. Enzyme Architecture: Breaking Down the Catalytic Cage that Activates Orotidine 5'-Monophosphate Decarboxylase for Catalysis. J Am Chem Soc 2018; 140:17580-17590. [PMID: 30475611 PMCID: PMC6317530 DOI: 10.1021/jacs.8b09609] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report the results of a study of the catalytic role of a network of four interacting amino acid side chains at yeast orotidine 5'-monophosphate decarboxylase ( ScOMPDC), by the stepwise replacement of all four side chains. The H-bond, which links the -CH2OH side chain of S154 from the pyrimidine umbrella loop of ScOMPDC to the amide side chain of Q215 in the phosphodianion gripper loop, creates a protein cage for the substrate OMP. The role of this interaction in optimizing transition state stabilization from the dianion gripper side chains Q215, Y217, and R235 was probed by determining the kinetic parameter kcat/ Km for 16 enzyme variants, which include all combinations of single, double, triple, and quadruple S154A, Q215A, Y217F, and R235A mutations. The effects of consecutive Q215A, Y217F, and R235A mutations on Δ G⧧ for wild-type enzyme-catalyzed decarboxylation sum to 11.6 kcal/mol, but to only 7.6 kcal/mol when starting from S154A mutant. This shows that the S154A mutation results in a (11.6-7.6) = 4.0 kcal/mol decrease in transition state stabilization from interactions with Q215, Y217, and R235. Mutant cycles show that ca. 2 kcal/mol of this 4 kcal/mol effect is from the direct interaction between the S154 and Q215 side chains and that ca. 2 kcal/mol is from a tightening in the stabilizing interactions of the Y217 and R235 side chains. The sum of the effects of individual A154S, A215Q, F217Y and A235R substitutions at the quadruple mutant of ScOMPDC to give the corresponding triple mutants, 5.6 kcal/mol, is much smaller than 16.0 kcal/mol, the sum of the effects of the related four substitutions in wild-type ScOMPDC to give the respective single mutants. The small effect of substitutions at the quadruple mutant is consistent with a large entropic cost to holding the flexible loops of ScOMPDC in the active closed conformation.
Collapse
Affiliation(s)
- Archie C Reyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - David C Plache
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Astrid P Koudelka
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Tina L Amyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - John A Gerlt
- Department of Chemistry and Biochemistry , University of Illinois , Urbana , Illinois 61801 , United States
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| |
Collapse
|
21
|
Dik DA, Batuecas MT, Lee M, Mahasenan KV, Marous DR, Lastochkin E, Fisher JF, Hermoso JA, Mobashery S. A Structural Dissection of the Active Site of the Lytic Transglycosylase MltE from Escherichia coli. Biochemistry 2018; 57:6090-6098. [PMID: 30256085 DOI: 10.1021/acs.biochem.8b00800] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Lytic transglycosylases (LTs) are bacterial enzymes that catalyze the cleavage of the glycan strands of the bacterial cell wall. The mechanism of this cleavage is a remarkable intramolecular transacetalization reaction, accomplished by an ensemble of active-site residues. Because the LT reaction occurs in parallel with the cell wall bond-forming reactions catalyzed by the penicillin-binding proteins, simultaneous inhibition of both enzymes can be particularly bactericidal to Gram-negative bacteria. The MltE lytic transglycosylase is the smallest of the eight LTs encoded by the Escherichia coli genome. Prior crystallographic and computational studies identified four active-site residues-E64, S73, S75, and Y192-as playing roles in catalysis. Each of these four residues was individually altered by mutation to give four variant enzymes (E64Q, S73A, S75A, and Y192F). All four variants showed reduced catalytic activity [soluble wild type (100%) > soluble Y192F and S75A (both 40%) > S73A (4%) > E64Q (≤1%)]. The crystal structure of each variant protein was determined at the resolution of 2.12 Å for E64Q, 2.33 Å for Y192F, 1.38 Å for S73A, and 1.35 Å for S75A. These variants show alteration of the hydrogen-bond interactions of the active site. Within the framework of a prior computational study of the LT mechanism, we suggest the mechanistic role of these four active-site residues in MltE catalysis.
Collapse
Affiliation(s)
- David A Dik
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - María T Batuecas
- Department of Crystallography and Structural Biology , Inst. Química-Física "Rocasolano", CSIC , Serrano 119 , 28006 Madrid , Spain
| | - Mijoon Lee
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - Kiran V Mahasenan
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - Daniel R Marous
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - Elena Lastochkin
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - Jed F Fisher
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| | - Juan A Hermoso
- Department of Crystallography and Structural Biology , Inst. Química-Física "Rocasolano", CSIC , Serrano 119 , 28006 Madrid , Spain
| | - Shahriar Mobashery
- Department of Chemistry and Biochemistry , University of Notre Dame , 352 McCourtney Hall , Notre Dame , Indiana 46556 , United States
| |
Collapse
|
22
|
Zhai X, Reinhardt CJ, Malabanan MM, Amyes TL, Richard JP. Enzyme Architecture: Amino Acid Side-Chains That Function To Optimize the Basicity of the Active Site Glutamate of Triosephosphate Isomerase. J Am Chem Soc 2018; 140:8277-8286. [PMID: 29862813 PMCID: PMC6037162 DOI: 10.1021/jacs.8b04367] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
We report pH rate profiles for kcat and Km for the
isomerization reaction
of glyceraldehyde 3-phosphate catalyzed by wildtype triosephosphate
isomerase (TIM) from three organisms and by ten mutants of TIM; and,
for Ki for inhibition of this reaction
by phosphoglycolate trianion (I3–). The pH profiles for Ki show
that the binding of I3– to TIM (E) to form EH·I3– is accompanied by
uptake of a proton by the carboxylate side-chain of E165, whose function
is to abstract a proton from substrate. The complexes for several
mutants exist mainly as E–·I3– at high pH, in which cases the pH profiles define the pKa for deprotonation of EH·I3–. The linear
free energy correlation, with slope of 0.73 (r2 = 0.96), between kcat/Km for TIM-catalyzed isomerization and the disassociation
constant of PGA trianion for TIM shows that EH·I3– and the
transition state are stabilized by similar interactions with the protein
catalyst. Values of pKa = 10–10.5
were estimated for deprotonation of EH·I3– for wildtype TIM.
This pKa decreases to as low as 6.3 for
the severely crippled Y208F mutant. There is a correlation between
the effect of several mutations on kcat/Km and on pKa for EH·I3–. The results support a model where the strong basicity of
E165 at the complex to the enediolate reaction intermediate is promoted
by side-chains from Y208 and S211, which serve to clamp loop 6 over
the substrate; I170, which assists in the creation of a hydrophobic
environment for E165; and P166, which functions in driving the carboxylate
side-chain of E165 toward enzyme-bound substrate.
Collapse
Affiliation(s)
- Xiang Zhai
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 United States
| | - Christopher J Reinhardt
- Department of Chemistry , University of Illinois at Urbana-Champaign , 600 S Mathews Avenue , Urbana , Illinois 61801 , United States
| | - M Merced Malabanan
- Department of Biochemistry , Vanderbilt University , 842 Robinson Research Building , Nashville , Tennessee 37205 , United States
| | - Tina L Amyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 United States
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 United States
| |
Collapse
|
23
|
Kulkarni YS, Liao Q, Byléhn F, Amyes TL, Richard JP, Kamerlin SCL. Role of Ligand-Driven Conformational Changes in Enzyme Catalysis: Modeling the Reactivity of the Catalytic Cage of Triosephosphate Isomerase. J Am Chem Soc 2018. [PMID: 29516737 PMCID: PMC5867644 DOI: 10.1021/jacs.8b00251] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
![]()
We have previously performed empirical
valence bond calculations
of the kinetic activation barriers, ΔG‡calc, for the deprotonation of complexes
between TIM and the whole substrate glyceraldehyde-3-phosphate (GAP, Kulkarni et al.2017, 139, 10514–1052528683550). We now extend
this work to also study the deprotonation of the substrate pieces
glycolaldehyde (GA) and GA·HPi [HPi = phosphite
dianion]. Our combined calculations provide activation barriers, ΔG‡calc, for the TIM-catalyzed
deprotonation of GAP (12.9 ± 0.8 kcal·mol–1), of the substrate piece GA (15.0 ± 2.4 kcal·mol–1), and of the pieces GA·HPi (15.5 ± 3.5 kcal·mol–1). The effect of bound dianion on ΔG‡calc is small (≤2.6 kcal·mol–1), in comparison to the much larger 12.0 and 5.8 kcal·mol–1 intrinsic phosphodianion and phosphite dianion binding
energy utilized to stabilize the transition states for TIM-catalyzed
deprotonation of GAP and GA·HPi, respectively. This
shows that the dianion binding energy is essentially fully expressed
at our protein model for the Michaelis complex, where it is utilized
to drive an activating change in enzyme conformation. The results
represent an example of the synergistic use of results from experiments
and calculations to advance our understanding of enzymatic reaction
mechanisms.
Collapse
Affiliation(s)
- Yashraj S Kulkarni
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden
| | - Qinghua Liao
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden
| | - Fabian Byléhn
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden.,Department of Chemical Engineering , University College London , Torrington Place , London WC1E 7JE , United Kingdom
| | - Tina L Amyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Shina C L Kamerlin
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 24 Uppsala , Sweden
| |
Collapse
|
24
|
Reyes AC, Amyes TL, Richard JP. A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer. Org Biomol Chem 2018; 15:8856-8866. [PMID: 28956050 DOI: 10.1039/c7ob01652b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
There is no consensus of opinion on the origin of the large rate accelerations observed for enzyme-catalyzed hydride transfer. The interpretation of recent results from studies on hydride transfer reactions catalyzed by alcohol dehydrogenase (ADH) focus on the proposal that the effective barrier height is reduced by quantum-mechanical tunneling through the energy barrier. This interpretation contrasts sharply with the notion that enzymatic rate accelerations are obtained through direct stabilization of the transition state for the nonenzymatic reaction in water. The binding energy of the dianion of substrate DHAP provides 11 kcal mol-1 stabilization of the transition state for the hydride transfer reaction catalyzed by glycerol-3-phosphate dehydrogenase (GPDH). We summarize evidence that the binding interactions between (GPDH) and dianion activators are utilized directly for stabilization of the transition state for enzyme-catalyzed hydride transfer. The possibility is considered, and then discounted, that these dianion binding interactions are utilized for the stabilization of a tunnel ready state (TRS) that enables efficient tunneling of the transferred hydride through the energy barrier, and underneath the energy maximum for the transition state. It is noted that the evidence to support the existence of a tunnel-ready state for the hydride transfer reactions catalyzed by ADH is ambiguous. We propose that the rate acceleration for ADH is due to the utilization of the binding energy of the cofactor NAD+/NADH in the stabilization of the transition state for enzyme-catalyzed hydride transfer.
Collapse
Affiliation(s)
- Archie C Reyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, USA.
| | | | | |
Collapse
|
25
|
He R, Reyes AC, Amyes TL, Richard JP. Enzyme Architecture: The Role of a Flexible Loop in Activation of Glycerol-3-phosphate Dehydrogenase for Catalysis of Hydride Transfer. Biochemistry 2018; 57:3227-3236. [PMID: 29337541 PMCID: PMC6001809 DOI: 10.1021/acs.biochem.7b01282] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
![]()
The side chain of Q295 of glycerol-3-phosphate
dehydrogenase from
human liver (hlGPDH) lies in a flexible loop, that
folds over the phosphodianion of substrate dihydroxyacetone phosphate
(DHAP). Q295 interacts with the side-chain cation from R269, which
is ion-paired to the substrate phosphodianion. Kinetic parameters kcat/Km (M–1 s–1) and kcat/KGAKHPi (M–2 s–1) were determined, respectively, for catalysis
of the reduction of DHAP and for dianion activation of catalysis of
reduction of glycolaldehyde (GA) catalyzed by wild-type, Q295G, Q295S,
Q295A, and Q295N mutants of hlGPDH. These mutations
result in up to a 150-fold decrease in (kcat/Km)DHAP and up to a 2.7 kcal/mol
decrease in the intrinsic phosphodianion binding energy. The data
define a linear correlation with slope 1.1, between the intrinsic
phosphodianion binding energy and the intrinsic phosphite dianion
binding energy for activation of hlGPDH-catalyzed
reduction of GA, that demonstrates a role for Q295 in optimizing this
dianion binding energy. The R269A mutation of wild-type GPDH results
in a 9.1 kcal/mol destabilization of the transition state for reduction
of DHAP, but the same R269A mutation of N270A and Q295A mutants result
in smaller 5.9 and 4.9 kcal/mol transition-state destabilization.
Similarly, the N270A or Q295A mutations of R269A GPDH each result
in large falloffs in the efficiency of rescue of the R269A mutant
by guanidine cation. We conclude that N270, which interacts for the
substrate phosphodianion and Q295, which interacts with the guanidine
side chain of R269, function to optimize the apparent
transition-state stabilization provided by the cationic side chain
of R269.
Collapse
Affiliation(s)
- Rui He
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Archie C Reyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - Tina L Amyes
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| | - John P Richard
- Department of Chemistry , University at Buffalo, SUNY , Buffalo , New York 14260-3000 , United States
| |
Collapse
|
26
|
Zheng F, Tu T, Wang X, Wang Y, Ma R, Su X, Xie X, Yao B, Luo H. Enhancing the catalytic activity of a novel GH5 cellulase GtCel5 from Gloeophyllum trabeum CBS 900.73 by site-directed mutagenesis on loop 6. BIOTECHNOLOGY FOR BIOFUELS 2018; 11:76. [PMID: 29588661 PMCID: PMC5863444 DOI: 10.1186/s13068-018-1080-5] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2017] [Accepted: 03/14/2018] [Indexed: 05/17/2023]
Abstract
BACKGROUND Cellulases of glycosyl hydrolase (GH) family 5 share a (β/α)8 TIM-barrel fold structure with eight βα loops surrounding the catalytic pocket. These loops exposed on the surface play a vital role in protein functions, primarily due to the interactions of some key amino acids with solvent and ligand molecules. It has been reported that motions of these loops facilitate substrate access and product release, and loops 6 and 7 located at the substrate entrance of the binding pocket promote proton transfer reaction at the catalytic site motions. However, the role of these flexible loops in catalysis of GH5 cellulase remains to be explored. RESULTS In the present study, an acidic, mesophilic GH5 cellulase (with optimal activity at pH 4.0 and 70 °C), GtCel5, was identified in Gloeophyllum trabeum CBS 900.73. The specific activities of GtCel5 toward CMC-Na, barley β-glucan, and lichenan were 1117 ± 43, 6257 ± 26 and 5318 ± 54 U/mg, respectively. Multiple sequence alignment indicates that one amino acid residue at position 233 on the loop 6 shows semi-conservativeness and might contribute to the great catalytic performance. Saturation mutagenesis at position 233 was then conducted to reveal the vital roles of this position in enzyme properties. In comparison to the wild type, variants N233A and N233G showed decreased optimal temperature (- 10 °C) but increased activities (27 and 70%) and catalytic efficiencies (kcat/Km; 45 and 52%), respectively. The similar roles of position 233 in catalytic performance were also verified in the other two GH5 homologs, TeEgl5A and PoCel5, by reverse mutation. Further molecular dynamics simulations suggested that the substitution of asparagine with alanine or glycine may introduce more hydrogen bonds, increase the flexibility of loop 6, enhance the interactions between enzyme and substrate, and thus improve the substrate affinity and catalytic efficiency. CONCLUSION This study proposed a novel cellulase with potentials for industrial application. A specific position was identified to play key roles in cellulase-substrate interactions and enzyme catalysis. It is of great importance for understanding the binding mechanism of GH5 cellulases, and provides an effective strategy to improve the catalytic performance of cellulases.
Collapse
Affiliation(s)
- Fei Zheng
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
- College of Biological Sciences and Biotechnology, Beijing Forestry University, Beijing, 100083 People’s Republic of China
| | - Tao Tu
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| | - Xiaoyu Wang
- College of Biological Sciences and Biotechnology, Beijing Forestry University, Beijing, 100083 People’s Republic of China
| | - Yuan Wang
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| | - Rui Ma
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| | - Xiaoyun Su
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| | - Xiangming Xie
- College of Biological Sciences and Biotechnology, Beijing Forestry University, Beijing, 100083 People’s Republic of China
| | - Bin Yao
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| | - Huiying Luo
- Key Laboratory for Feed Biotechnology of the Ministry of Agriculture, Feed Research Institute, Chinese Academy of Agricultural Sciences, No. 12 South Zhongguancun Street, Haidian, Beijing, 100081 People’s Republic of China
| |
Collapse
|
27
|
Kulkarni Y, Liao Q, Petrović D, Krüger DM, Strodel B, Amyes TL, Richard JP, Kamerlin SCL. Enzyme Architecture: Modeling the Operation of a Hydrophobic Clamp in Catalysis by Triosephosphate Isomerase. J Am Chem Soc 2017; 139:10514-10525. [PMID: 28683550 PMCID: PMC5543394 DOI: 10.1021/jacs.7b05576] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Indexed: 12/22/2022]
Abstract
Triosephosphate isomerase (TIM) is a proficient catalyst of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde phosphate (GAP), via general base catalysis by E165. Historically, this enzyme has been an extremely important model system for understanding the fundamentals of biological catalysis. TIM is activated through an energetically demanding conformational change, which helps position the side chains of two key hydrophobic residues (I170 and L230), over the carboxylate side chain of E165. This is critical both for creating a hydrophobic pocket for the catalytic base and for maintaining correct active site architecture. Truncation of these residues to alanine causes significant falloffs in TIM's catalytic activity, but experiments have failed to provide a full description of the action of this clamp in promoting substrate deprotonation. We perform here detailed empirical valence bond calculations of the TIM-catalyzed deprotonation of DHAP and GAP by both wild-type TIM and its I170A, L230A, and I170A/L230A mutants, obtaining exceptional quantitative agreement with experiment. Our calculations provide a linear free energy relationship, with slope 0.8, between the activation barriers and Gibbs free energies for these TIM-catalyzed reactions. We conclude that these clamping side chains minimize the Gibbs free energy for substrate deprotonation, and that the effects on reaction driving force are largely expressed at the transition state for proton transfer. Our combined analysis of previous experimental and current computational results allows us to provide an overview of the breakdown of ground-state and transition state effects in enzyme catalysis in unprecedented detail, providing a molecular description of the operation of a hydrophobic clamp in triosephosphate isomerase.
Collapse
Affiliation(s)
- Yashraj
S. Kulkarni
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, Uppsala S-751 24, Sweden
| | - Qinghua Liao
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, Uppsala S-751 24, Sweden
| | - Dušan Petrović
- Institute
of Complex Systems: Structural Biochemistry, Forschungszentrum Jülich, Jülich 52425, Germany
| | - Dennis M. Krüger
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, Uppsala S-751 24, Sweden
| | - Birgit Strodel
- Institute
of Complex Systems: Structural Biochemistry, Forschungszentrum Jülich, Jülich 52425, Germany
- Institute
of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitätsstrasse 1, Düsseldorf 40225, Germany
| | - Tina L. Amyes
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - John P. Richard
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United States
| | - Shina C. L. Kamerlin
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, Uppsala S-751 24, Sweden
| |
Collapse
|
28
|
Amyes TL, Malabanan MM, Zhai X, Reyes AC, Richard JP. Enzyme activation through the utilization of intrinsic dianion binding energy. Protein Eng Des Sel 2017; 30:157-165. [PMID: 27903763 DOI: 10.1093/protein/gzw064] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2016] [Accepted: 11/14/2016] [Indexed: 11/12/2022] Open
Abstract
43 We consider 'the proposition that the intrinsic binding energy that results from the noncovalent interaction of a specific substrate with the active site of the enzyme is considerably larger than is generally believed. An important part of this binding energy may be utilized to provide the driving force for catalysis, so that the observed binding energy represents only what is left over after this utilization' [Jencks,W.P. (1975) Adv. Enzymol. Relat. Areas. Mol. Biol. , , 219-410]. The large ~12 kcal/mol intrinsic substrate phosphodianion binding energy for reactions catalyzed by triosephosphate isomerase (TIM), orotidine 5'-monophosphate decarboxylase and glycerol-3-phosphate dehydrogenase is divided into 4-6 kcal/mol binding energy that is expressed on the formation of the Michaelis complex in anchoring substrates to the respective enzyme, and 6-8 kcal/mol binding energy that is specifically expressed at the transition state in activating the respective enzymes for catalysis. A structure-based mechanism is described where the dianion binding energy drives a conformational change that activates these enzymes for catalysis. Phosphite dianion plays the active role of holding TIM in a high-energy closed active form, but acts as passive spectator in showing no effect on transition-state structure. The result of studies on mutant enzymes is presented, which support the proposal that the dianion-driven enzyme conformational change plays a role in enhancing the basicity of side chain of E167, the catalytic base, by clamping the base between a pair of hydrophobic side chains. The insight these results provide into the architecture of enzyme active sites and the development of strategies for the de novo design of protein catalysts is discussed.
Collapse
Affiliation(s)
- T L Amyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, NY 14260-3000, USA
| | - M M Malabanan
- Department of Biochemistry, Vanderbilt University, Nashville, TN37205-0146, USA
| | - X Zhai
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX77843-2128, USA
| | - A C Reyes
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, NY 14260-3000, USA
| | - J P Richard
- Department of Chemistry, University at Buffalo, SUNY, Buffalo, NY 14260-3000, USA
| |
Collapse
|
29
|
Amrein BA, Steffen-Munsberg F, Szeler I, Purg M, Kulkarni Y, Kamerlin SCL. CADEE: Computer-Aided Directed Evolution of Enzymes. IUCRJ 2017; 4:50-64. [PMID: 28250941 PMCID: PMC5331465 DOI: 10.1107/s2052252516018017] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/12/2016] [Accepted: 11/09/2016] [Indexed: 05/10/2023]
Abstract
The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.
Collapse
Affiliation(s)
- Beat Anton Amrein
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Fabian Steffen-Munsberg
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Ireneusz Szeler
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Miha Purg
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Yashraj Kulkarni
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Shina Caroline Lynn Kamerlin
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| |
Collapse
|
30
|
Reyes AC, Amyes TL, Richard JP. Enzyme Architecture: Self-Assembly of Enzyme and Substrate Pieces of Glycerol-3-Phosphate Dehydrogenase into a Robust Catalyst of Hydride Transfer. J Am Chem Soc 2016; 138:15251-15259. [PMID: 27792325 PMCID: PMC5291162 DOI: 10.1021/jacs.6b09936] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The stabilization of the transition
state for hlGPDH-catalyzed reduction of DHAP due
to the action of the phosphodianion
of DHAP and the cationic side chain of R269 is between 12.4 and 17
kcal/mol. The R269A mutation of glycerol-3-phosphate dehydrogenase
(hlGPDH) results in a 9.1 kcal/mol destabilization
of the transition state for enzyme-catalyzed reduction of dihydroxyacetone
phosphate (DHAP) by NADH, and there is a 6.7 kcal/mol stabilization of this transition state by 1.0 M guanidine cation (Gua+) [J. Am. Chem. Soc.2015, 137, 5312–5315]. The R269A mutant shows no detectable
activity toward reduction of glycolaldehyde (GA), or activation of
this reaction by 30 mM HPO32–. We report
the unprecedented self-assembly of R269A hlGPDH,
dianions (X2– = FPO32–, HPO32–, or SO42–), Gua+ and GA into a functioning catalyst of the reduction
of GA, and fourth-order reaction rate constants kcat/KGAKXKGua. The linear logarithmic correlation
(slope = 1.0) between values of kcat/KGAKX for dianion
activation of wildtype hlGPDH-catalyzed reduction
of GA and kcat/KGAKXKGua shows that the electrostatic interaction between exogenous dianions
and the side chain of R269 is not significantly perturbed by cutting hlGPDH into R269A and Gua+ pieces. The advantage
for connection of hlGPDH (R269A mutant + Gua+) and substrate pieces (GA + HPi) pieces, (ΔGS‡)HPi+E+Gua = 5.6 kcal/mol, is nearly equal to the sum
of the advantage to connection of the substrate pieces, (ΔGS‡)GA+HPi = 3.3 kcal/mol, for wildtype hlGPDH-catalyzed reaction of GA + HPi, and for connection
of the enzyme pieces, (ΔGS‡)E+Gua = 2.4
kcal/mol, for Gua+ activation of the R269A hlGPDH-catalyzed reaction of DHAP.
Collapse
Affiliation(s)
- Archie C Reyes
- Department of Chemistry, University at Buffalo, SUNY , Buffalo, New York 14260-3000, United States
| | - Tina L Amyes
- Department of Chemistry, University at Buffalo, SUNY , Buffalo, New York 14260-3000, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo, SUNY , Buffalo, New York 14260-3000, United States
| |
Collapse
|
31
|
Richard JP, Amyes TL, Malabanan MM, Zhai X, Kim KJ, Reinhardt CJ, Wierenga RK, Drake EJ, Gulick AM. Structure-Function Studies of Hydrophobic Residues That Clamp a Basic Glutamate Side Chain during Catalysis by Triosephosphate Isomerase. Biochemistry 2016; 55:3036-47. [PMID: 27149328 PMCID: PMC4934371 DOI: 10.1021/acs.biochem.6b00311] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Kinetic
parameters are reported for the reactions of whole substrates
(kcat/Km,
M–1 s–1) (R)-glyceraldehyde
3-phosphate (GAP) and
dihydroxyacetone phosphate (DHAP) and for the substrate pieces [(kcat/Km)E·HPi/Kd, M–2 s–1] glycolaldehyde (GA) and phosphite dianion
(HPi) catalyzed by the I172A/L232A mutant of triosephosphate
isomerase
from Trypanosoma brucei brucei (TbbTIM). A comparison with the corresponding parameters for wild-type,
I172A, and L232A TbbTIM-catalyzed reactions shows
that the effect of I172A and L232A mutations on ΔG⧧ for the wild-type TbbTIM-catalyzed
reactions of the substrate pieces is nearly the same
as the effect of the same mutations on TbbTIM previously
mutated at the second side chain. This provides strong evidence that
mutation of the first hydrophobic side chain does not affect the functioning
of the second side chain in catalysis of the reactions of the substrate
pieces. By contrast, the effects of I172A and L232A mutations on ΔG⧧ for wild-type TbbTIM-catalyzed
reactions of the whole substrate are different from
the effect of the same mutations on TbbTIM previously
mutated at the second side chain. This is due to the change in the
rate-determining step that determines the barrier to the isomerization
reaction. X-ray crystal structures are reported for I172A, L232A,
and I172A/L232A TIMs and for the complexes of these mutants to the
intermediate analogue phosphoglycolate (PGA). The structures of the
PGA complexes with wild-type and mutant enzymes are nearly superimposable,
except that the space opened by replacement of the hydrophobic side
chain is occupied by a water molecule that lies ∼3.5 Å
from the basic side chain of Glu167. The new water at I172A mutant TbbTIM provides a simple rationalization for the increase
in the activation barrier ΔG⧧ observed for mutant enzyme-catalyzed
reactions of the whole substrate and substrate pieces. By contrast,
the new water at the L232A mutant does not predict the decrease in
ΔG⧧ observed for the mutant
enzyme-catalyzed
reactions of the substrate piece GA.
Collapse
Affiliation(s)
- John P Richard
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - Tina L Amyes
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - M Merced Malabanan
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - Xiang Zhai
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - Kalvin J Kim
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - Christopher J Reinhardt
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260, United States
| | - Rik K Wierenga
- Department of Biochemistry and Biocenter, University of Oulu , P.O. Box 3000, FIN-90014 Oulu, Finland
| | - Eric J Drake
- Hauptman-Woodward Institute , 700 Ellicott Street, Buffalo, New York 14203, United States.,Department of Structural Biology, University at Buffalo, State University of New York , Buffalo, New York 14203, United States
| | - Andrew M Gulick
- Hauptman-Woodward Institute , 700 Ellicott Street, Buffalo, New York 14203, United States.,Department of Structural Biology, University at Buffalo, State University of New York , Buffalo, New York 14203, United States
| |
Collapse
|
32
|
Reyes AC, Amyes TL, Richard JP. Enzyme Architecture: A Startling Role for Asn270 in Glycerol 3-Phosphate Dehydrogenase-Catalyzed Hydride Transfer. Biochemistry 2016; 55:1429-32. [PMID: 26926520 PMCID: PMC4812618 DOI: 10.1021/acs.biochem.6b00116] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The
side chains of R269 and N270 interact with the phosphodianion
of dihydroxyacetone phosphate (DHAP) bound to glycerol 3-phosphate
dehydrogenase (GPDH). The R269A, N270A, and R269A/N270A mutations
of GPDH result in 9.1, 5.6, and 11.5 kcal/mol destabilization, respectively,
of the transition state for GPDH-catalyzed reduction of DHAP by the
reduced form of nicotinamide adenine dinucleotide. The N270A mutation
results in a 7.7 kcal/mol decrease in the intrinsic phosphodianion
binding energy, which is larger than the 5.6 kcal/mol effect of the
mutation on the stability of the transition state for reduction of
DHAP; a 2.2 kcal/mol stabilization of the transition state for unactivated
hydride transfer to the truncated substrate glycolaldehyde (GA); and
a change in the effect of phosphite dianion on GPDH-catalyzed reduction
of GA, from strongly activating to inhibiting. The N270A mutation
breaks the network of hydrogen bonding side chains, Asn270, Thr264,
Asn205, Lys204, Asp260, and Lys120, which connect the dianion activation
and catalytic sites of GPDH. We propose that this disruption dramatically
alters the performance of GPDH at these sites.
Collapse
Affiliation(s)
- Archie C Reyes
- Department of Chemistry, University at Buffalo , Buffalo, New York 14260-3000, United States
| | - Tina L Amyes
- Department of Chemistry, University at Buffalo , Buffalo, New York 14260-3000, United States
| | - John P Richard
- Department of Chemistry, University at Buffalo , Buffalo, New York 14260-3000, United States
| |
Collapse
|
33
|
Xiong YY, Li JQ, Yan CS, Gao HY, Zhou JP, Gong LL, Luo MB, Zhang L, Meng PP, Luo F. MOF catalysis of FeII-to-FeIII reaction for an ultrafast and one-step generation of the Fe2O3@MOF composite and uranium(vi) reduction by iron(ii) under ambient conditions. Chem Commun (Camb) 2016; 52:9538-41. [DOI: 10.1039/c6cc04597a] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Surface is for the first time attested to be one of outstanding merits of metal–organic frameworks (MOFs).
Collapse
|