1
|
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
|
2
|
Abstract
In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.
Collapse
|
3
|
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
|
4
|
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
|
5
|
The Role of the Substrate Phosphodianion in Catalysis by Orotidine 5'-Monophosphate Decarboxylase. Biochemistry 2023; 62:969-970. [PMID: 36791154 PMCID: PMC10052792 DOI: 10.1021/acs.biochem.3c00031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
|
6
|
Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations. J Mol Graph Model 2023; 119:108402. [PMID: 36610324 DOI: 10.1016/j.jmgm.2022.108402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 12/17/2022] [Accepted: 12/26/2022] [Indexed: 12/31/2022]
Abstract
The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.
Collapse
|
7
|
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
|
8
|
Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
|
9
|
An Enzyme with High Catalytic Proficiency Utilizes Distal Site Substrate Binding Energy to Stabilize the Closed State but at the Expense of Substrate Inhibition. ACS Catal 2022; 12:3149-3164. [PMID: 35692864 PMCID: PMC9171722 DOI: 10.1021/acscatal.1c05524] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 02/10/2022] [Indexed: 02/05/2023]
Abstract
Understanding the factors that underpin the enormous catalytic proficiencies of enzymes is fundamental to catalysis and enzyme design. Enzymes are, in part, able to achieve high catalytic proficiencies by utilizing the binding energy derived from nonreacting portions of the substrate. In particular, enzymes with substrates containing a nonreacting phosphodianion group coordinated in a distal site have been suggested to exploit this binding energy primarily to facilitate a conformational change from an open inactive form to a closed active form, rather than to either induce ground state destabilization or stabilize the transition state. However, detailed structural evidence for the model is limited. Here, we use β-phosphoglucomutase (βPGM) to investigate the relationship between binding a phosphodianion group in a distal site, the adoption of a closed enzyme form, and catalytic proficiency. βPGM catalyzes the isomerization of β-glucose 1-phosphate to glucose 6-phosphate via phosphoryl transfer reactions in the proximal site, while coordinating a phosphodianion group of the substrate(s) in a distal site. βPGM has one of the largest catalytic proficiencies measured and undergoes significant domain closure during its catalytic cycle. We find that side chain substitution at the distal site results in decreased substrate binding that destabilizes the closed active form but is not sufficient to preclude the adoption of a fully closed, near-transition state conformation. Furthermore, we reveal that binding of a phosphodianion group in the distal site stimulates domain closure even in the absence of a transferring phosphoryl group in the proximal site, explaining the previously reported β-glucose 1-phosphate inhibition. Finally, our results support a trend whereby enzymes with high catalytic proficiencies involving phosphorylated substrates exhibit a greater requirement to stabilize the closed active form.
Collapse
|
10
|
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
|
11
|
How enzymes harness highly unfavorable proton transfer reactions. Protein Sci 2021; 30:735-744. [PMID: 33554401 DOI: 10.1002/pro.4037] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Revised: 02/01/2021] [Accepted: 02/05/2021] [Indexed: 11/12/2022]
Abstract
Acid-base reactions that are exceedingly unfavorable under standard conditions can be catalytically important at enzyme active sites. For example, in triose phosphate isomerase, a glutamate side chain (nominal pKa ≈ 4 in solution) can in fact deprotonate a CH group that is vicinal to a carbonyl (pKa ≈ 18 in solution). This is true because of three distinct interactions: (a) ground state pKa shifts due to environment polarity and electrostatics; (b) dramatic increases in effective molarity due to optimization of proximity and orientation; and (c) transition state pKa shifts due to binding interactions and the formation of strong low barrier hydrogen bonds. In this report, we review the literature showing that the sum of these three effects supplies more than enough free energy to push forward proton transfer reactions that under standard conditions are exceedingly nonspontaneous and slow.
Collapse
|
12
|
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
|
13
|
Modeling the Role of a Flexible Loop and Active Site Side Chains in Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase. ACS Catal 2020; 10:11253-11267. [PMID: 33042609 PMCID: PMC7536716 DOI: 10.1021/acscatal.0c02757] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Revised: 08/31/2020] [Indexed: 11/30/2022]
Abstract
![]()
Glycerol-3-phosphate
dehydrogenase is a biomedically important
enzyme that plays a crucial role in lipid biosynthesis. It is activated
by a ligand-gated conformational change that is necessary for the
enzyme to reach a catalytically competent conformation capable of
efficient transition-state stabilization. While the human form (hlGPDH) has been the subject of extensive structural and
biochemical studies, corresponding computational studies to support
and extend experimental observations have been lacking. We perform
here detailed empirical valence bond and Hamiltonian replica exchange
molecular dynamics simulations of wild-type hlGPDH
and its variants, as well as providing a crystal structure of the
binary hlGPDH·NAD R269A variant where the enzyme
is present in the open conformation. We estimated the activation free
energies for the hydride transfer reaction in wild-type and substituted hlGPDH and investigated the effect of mutations on catalysis
from a detailed structural study. In particular, the K120A and R269A
variants increase both the volume and solvent exposure of the active
site, with concomitant loss of catalytic activity. In addition, the
R269 side chain interacts with both the Q295 side chain on the catalytic
loop, and the substrate phosphodianion. Our structural data and simulations
illustrate the critical role of this side chain in facilitating the
closure of hlGPDH into a catalytically competent
conformation, through modulating the flexibility of a key catalytic
loop (292-LNGQKL-297). This, in turn, rationalizes a tremendous 41,000
fold decrease experimentally in the turnover number, kcat, upon truncating this residue, as loop closure is
essential for both correct positioning of key catalytic residues in
the active site, as well as sequestering the active site from the
solvent. Taken together, our data highlight the importance of this
ligand-gated conformational change in catalysis, a feature that can
be exploited both for protein engineering and for the design of allosteric
inhibitors targeting this biomedically important enzyme.
Collapse
|
14
|
Abstract
Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.
Collapse
|
15
|
Enhancing a de novo enzyme activity by computationally-focused ultra-low-throughput screening. Chem Sci 2020; 11:6134-6148. [PMID: 32832059 PMCID: PMC7407621 DOI: 10.1039/d0sc01935f] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2020] [Accepted: 05/18/2020] [Indexed: 01/02/2023] Open
Abstract
Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers (∼2 × 104 M-1 s-1 and ∼102 s-1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within ∼2 kcal mol-1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stability-guidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design.
Collapse
|
16
|
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
|
17
|
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: 11] [Impact Index Per Article: 2.2] [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
|
18
|
Computational Studies of Catalytic Loop Dynamics in Yersinia Protein Tyrosine Phosphatase Using Pathway Optimization Methods. J Phys Chem B 2019; 123:7840-7851. [PMID: 31437399 PMCID: PMC6752976 DOI: 10.1021/acs.jpcb.9b06759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Yersinia Protein Tyrosine Phosphatase (YopH) is the most efficient enzyme among all known PTPases and relies on its catalytic loop movements for substrate binding and catalysis. Fluorescence, NMR, and UV resonance Raman (UVRR) techniques have been used to study the thermodynamic and dynamic properties of the loop motions. In this study, a computational approach based on the pathway refinement methods nudged elastic band (NEB) and harmonic Fourier beads (HFB) has been developed to provide structural interpretations for the experimentally observed kinetic processes. In this approach, the minimum potential energy pathways for the loop open/closure conformational changes were determined by NEB using a one-dimensional global coordinate. Two dimensional data analyses of the NEB results were performed as an efficient method to qualitatively evaluate the energetics of transitions along several specific physical coordinates. The free energy barriers for these transitions were then determined more precisely using the HFB method. Kinetic parameters were estimated from the energy barriers using transition state theory and compared against experimentally determined kinetic parameters. When the calculated energy barriers are calibrated by a simple "scaling factor", as have been done in our previous vibrational frequency calculations to explain the ligand frequency shift upon its binding to protein, it is possible to make structural interpretations of several observed enzyme dynamic rates. For example, the nanosecond kinetics observed by fluorescence anisotropy may be assigned to the translational motion of the catalytic loop and microsecond kinetics observed in fluorescence T-jump can be assigned to the loop backbone dihedral angle flipping. Furthermore, we can predict that a Trp354 conformational conversion associated with the loop movements would occur on the tens of nanoseconds time scale, to be verified by future UVRR T-jump studies.
Collapse
|
19
|
Improving biocatalysis of cefaclor with penicillin acylase immobilized on magnetic nanocrystalline cellulose in deep eutectic solvent based co-solvent. BIORESOURCE TECHNOLOGY 2019; 288:121548. [PMID: 31152956 DOI: 10.1016/j.biortech.2019.121548] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Revised: 05/21/2019] [Accepted: 05/22/2019] [Indexed: 06/09/2023]
Abstract
Deep eutectic solvent (DES), has been considered as a new type of green solvent applied in enzymatic systems. Here, we reported DES-buffer co-solvent as a novel reaction medium for high efficient synthesis of cefaclor by penicilin acylase immobilized on magnetic nanocrystalline cellulose. Effect of DES composition, DES-buffer ratio, temperature, pH, substrate ratio and substrate concentration was systematically investigated. In co-solvent consisting of choline chloride (ChCl):glycol-buffer (7:3, v/v), conversion of 7-ACCA was 94%, synthesis to hydrolysis ratio was 1.8, and yield of cefaclor reached 91%, higher than that in aqueous buffer with optimized yield of 84%, showing the great potential of DES as organic solvent alternative. To the best of our knowledge, this is the first example of biosynthesis of cefaclor in the DES-buffer co-solvent.
Collapse
|
20
|
|
21
|
Thermodynamic analysis of remote substrate binding energy in 3α-hydroxysteroid dehydrogenase/carbonyl reductase catalysis. Chem Biol Interact 2019; 302:183-189. [PMID: 30794798 DOI: 10.1016/j.cbi.2019.02.011] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2018] [Revised: 01/30/2019] [Accepted: 02/14/2019] [Indexed: 11/16/2022]
Abstract
The binding energy of enzyme and substrate is used to lower the activation energy for the catalytic reaction. 3α-HSD/CR uses remote binding interactions to accelerate the reaction of androsterone with NAD+. Here, we examine the enthalpic and entropic components of the remote binding energy in the 3α-HSD/CR-catalyzed reaction of NAD+ with androsterone versus the substrate analogs, 2-decalol and cyclohexanol, by analyzing the temperature-dependent kinetic parameters through steady-state kinetics. The effects of temperature on kcat/Km for 3α-HSD/CR acting on androsterone, 2-decalol, and cyclohexanol show the reactions are entropically favorable but enthalpically unfavorable. Thermodynamic analysis from the temperature-dependent values of Km and kcat shows the binding of the E-NAD+ complex with either 2-decalol or cyclohexanol to form the ternary complex is endothermic and entropy-driven, and the subsequent conversion to the transition state is both enthalpically and entropically unfavorable. Hence, solvation entropy may play an important role in the binding process through both the desolvation of the solute molecules and the release of bound water molecules from the active site into bulk solvent. As compared to the thermodynamic parameters of 3α-HSD/CR acting on cyclohexanol, the hydrophobic interaction of the B-ring of steroids with the active site of 3α-HSD/CR contributes to catalysis by increasing exclusively the entropy of activation (ΔTΔS‡ = 1.8 kcal/mol), while the BCD-ring of androsterone significantly lowers ΔΔH‡ by 10.4 kcal/mol with a slight entropic penalty of -1.9 kcal/mol. Therefore, the remote non-reacting sites of androsterone may induce a conformational change of the substrate binding loop with an entropic cost for better interaction with the transition state to decrease the enthalpy of activation, significantly increasing catalytic efficiency.
Collapse
|
22
|
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
|
23
|
Computational physical organic chemistry using the empirical valence bond approach. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.07.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
24
|
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
|
25
|
Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case. J Am Chem Soc 2018; 140:15889-15903. [PMID: 30362343 DOI: 10.1021/jacs.8b09378] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.
Collapse
|
26
|
Evolutionary Effects on Bound Substrate pKa in Dihydrofolate Reductase. J Am Chem Soc 2018; 140:16650-16660. [DOI: 10.1021/jacs.8b09089] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
27
|
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
|