Quantum and classical simulations of orotidine monophosphate decarboxylase: support for a direct decarboxylation mechanism

Adenylate Kinase-Catalyzed Reactions of AMP in Pieces: Specificity for Catalysis at the Nucleoside Activator and Dianion Catalytic Sites

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.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

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.

Enzyme Architecture: Breaking Down the Catalytic Cage that Activates Orotidine 5'-Monophosphate Decarboxylase for Catalysis

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 -CH₂OH 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 k_cat/ K_m 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.

Re-refinement of Plasmodium falciparum orotidine 5'-monophosphate decarboxylase provides a clearer picture of an important malarial drug target

The development of antimalarial drugs remains a public health priority, and the orotidine 5'-monophosphate decarboxylase from Plasmodium falciparum (PfOMPDC) has great potential as a drug target. The crystallization of PfOMPDC with substrate bound represents an important advance for structure-based drug-design efforts [Tokuoka et al. (2008), J. Biochem. 143, 69-78]. The complex of the enzyme bound to the substrate OMP (PDB entry 2za1) would be of particular utility in this regard. However, re-refinement of this structure of the Michaelis complex shows that the bound ligand is the product rather than the substrate. Here, the re-refinement of a set of three structures, the apo enzyme and two versions of the product-bound form (PDB entries 2za1, 2za2 and 2za3), is reported. The improved geometry and fit of these structures to the observed electron density will enhance their utility in antimalarial drug design.

Orotidine 5'-Monophosphate Decarboxylase: Probing the Limits of the Possible for Enzyme Catalysis

The mystery associated with catalysis by what were once regarded as protein black boxes, diminished with the X-ray crystallographic determination of the three-dimensional structures of enzyme–substrate complexes. The report that several high-resolution X-ray crystal structures of orotidine 5′-monophosphate decarboxylase (OMPDC) failed to provide a consensus mechanism for enzyme-catalyzed decarboxylation of OMP to form uridine 5′-monophosphate, therefore, provoked a flurry of controversy. This controversy was fueled by the enormous 10²³-fold rate acceleration for this enzyme, which had “jolted many biochemists’ assumptions about the catalytic potential of enzymes.” Our studies on the mechanism of action of OMPDC provide strong evidence that catalysis by this enzyme is not fundamentally different from less proficient catalysts, while highlighting important architectural elements that enable a peak level of performance. Many enzymes undergo substrate-induced protein conformational changes that trap their substrates in solvent occluded protein cages, but the conformational change induced by ligand binding to OMPDC is incredibly complex, as required to enable the development of 22 kcal/mol of stabilizing binding interactions with the phosphodianion and ribosyl substrate fragments of OMP. The binding energy from these fragments is utilized to activate OMPDC for catalysis of decarboxylation at the orotate fragment of OMP, through the creation of a tight, catalytically active, protein cage from the floppy, open, unliganded form of OMPDC. Such utilization of binding energy for ligand-driven conformational changes provides a general mechanism to obtain specificity in transition state binding. The rate enhancement that results from the binding of carbon acid substrates to enzymes is partly due to a reduction in the carbon acid pK_a that is associated with ligand binding. The binding of UMP to OMPDC results in an unusually large >12 unit decrease in the pK_a = 29 for abstraction of the C-6 substrate hydrogen, due to stabilization of an enzyme-bound vinyl carbanion, which is also an intermediate of OMPDC-catalyzed decarboxylation. The protein–ligand interactions operate to stabilize the vinyl carbanion at the enzyme active site compared to aqueous solution, rather than to stabilize the transition state for the concerted electrophilic displacement of CO₂ by H⁺ that avoids formation of this reaction intermediate. There is evidence that OMPDC induces strain into the bound substrate. The interaction between the amide side chain of Gln-215 from the phosphodianion gripper loop and the hydroxymethylene side chain of Ser-154 from the pyrimidine umbrella of ScOMPDC position the amide side chain to interact with the phosphodianion of OMP. There are no direct stabilizing interactions between dianion gripper protein side chains Gln-215, Tyr-217, and Arg-235 and the pyrimidine ring at the decarboxylation transition state. Rather these side chains function solely to hold OMPDC in the catalytically active closed conformation. The hydrophobic side chains that line the active site of OMPDC in the region of the departing CO₂ product may function to stabilize the decarboxylation transition state by providing hydrophobic solvation of this product.

A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer

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.

Enzyme Architecture: Erection of Active Orotidine 5'-Monophosphate Decarboxylase by Substrate-Induced Conformational Changes

Orotidine 5′-monophosphate decarboxylase (OMPDC) catalyzes the decarboxylation of 5-fluoroorotate (FO) with k_cat/K_m = 1.4 × 10^–7 M^–1 s^–1. Combining this and related kinetic parameters shows that the 31 kcal/mol stabilization of the transition state for decarboxylation of OMP provided by OMPDC represents the sum of 11.8 and 10.6 kcal/mol stabilization by the substrate phosphodianion and the ribosyl ring, respectively, and an 8.6 kcal/mol stabilization from the orotate ring. The transition state for OMPDC-catalyzed decarboxylation of FO is stabilized by 5.2, 7.2, and 9.0 kcal/mol, respectively, by 1.0 M phosphite dianion, d-glycerol 3-phosphate and d-erythritol 4-phosphate. The stabilization is due to the utilization of binding interactions of the substrate fragments to drive an enzyme conformational change, which locks the orotate ring of the whole substrate, or the substrate pieces in a caged complex. We propose that enzyme-activation is a possible, and perhaps probable, consequence of any substrate-induced enzyme conformational change.

Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects

Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.

Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

A reaction's transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.

Challenges in computational studies of enzyme structure, function and dynamics

In this review we give an overview of the field of Computational enzymology. We start by describing the birth of the field, with emphasis on the work of the 2013 chemistry Nobel Laureates. We then present key features of the state-of-the-art in the field, showing what theory, accompanied by experiments, has taught us so far about enzymes. We also briefly describe computational methods, such as quantum mechanics-molecular mechanics approaches, reaction coordinate treatment, and free energy simulation approaches. We finalize by discussing open questions and challenges.

Enzyme architecture: deconstruction of the enzyme-activating phosphodianion interactions of orotidine 5'-monophosphate decarboxylase

The mechanism for activation of orotidine 5'-monophosphate decarboxylase (OMPDC) by interactions of side chains from Gln215 and Try217 at a gripper loop and R235, adjacent to this loop, with the phosphodianion of OMP was probed by determining the kinetic parameters k(cat) and K(m) for all combinations of single, double, and triple Q215A, Y217F, and R235A mutations. The 12 kcal/mol intrinsic binding energy of the phosphodianion is shown to be equal to the sum of the binding energies of the side chains of R235 (6 kcal/mol), Q215 (2 kcal/mol), Y217 (2 kcal/mol), and hydrogen bonds to the G234 and R235 backbone amides (2 kcal/mol). Analysis of a triple mutant cube shows small (ca. 1 kcal/mol) interactions between phosphodianion gripper side chains, which are consistent with steric crowding of the side chains around the phosphodianion at wild-type OMPDC. These mutations result in the same change in the activation barrier to the OMPDC-catalyzed reactions of the whole substrate OMP and the substrate pieces (1-β-D-erythrofuranosyl)orotic acid (EO) and phosphite dianion. This shows that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The 12 kcal/mol intrinsic phosphodianion binding energy of OMP is divided between the 8 kcal/mol of binding energy, which is utilized to drive a thermodynamically unfavorable conformational change of the free enzyme, resulting in an increase in (k(cat))(obs) for OMPDC-catalyzed decarboxylation of OMP, and the 4 kcal/mol of binding energy, which is utilized to stabilize the Michaelis complex, resulting in a decrease in (K(m))(obs).

How Accurate Are Transition States from Simulations of Enzymatic Reactions?

The rate expression of traditional transition state theory (TST) assumes no recrossing of the transition state (TS) and thermal quasi-equilibrium between the ground state and the TS. Currently, it is not well understood to what extent these assumptions influence the nature of the activated complex obtained in traditional TST-based simulations of processes in the condensed phase in general and in enzymes in particular. Here we scrutinize these assumptions by characterizing the TSs for hydride transfer catalyzed by the enzyme Escherichia coli dihydrofolate reductase obtained using various simulation approaches. Specifically, we compare the TSs obtained with common TST-based methods and a dynamics-based method. Using a recently developed accurate hybrid quantum mechanics/molecular mechanics potential, we find that the TST-based and dynamics-based methods give considerably different TS ensembles. This discrepancy, which could be due equilibrium solvation effects and the nature of the reaction coordinate employed and its motion, raises major questions about how to interpret the TSs determined by common simulation methods. We conclude that further investigation is needed to characterize the impact of various TST assumptions on the TS phase-space ensemble and on the reaction kinetics.

Design of inhibitors of ODCase

ODCase is a highly proficient enzyme responsible for the decarboxylation of orotidine monophosphate to generate uridine monophosphate. ODCase has attracted early attention due to its interesting mechanism of catalysis. In order to exploit therapeutic advantages due to the inhibition of ODCase, one must have selective inhibitors of this enzyme from the pathogen, or a dysregulated molecular mechanism involving ODCase. ODCase inhibitors have potential applications as anticancer agents, antiviral agents, antimalarial agents and potentially act against other parasitic diseases. A variety of C6-substituted uridine monophosphate derivatives have shown excellent inhibition of ODCase. 6-iodouridine is a potent inhibitor of the malaria parasite, and its monophosphate form covalently inhibits ODCase. A variety of inhibitors of ODCase with potential applications as therapeutic agents are discussed in this review.

Enzyme architecture: the activating oxydianion binding domain for orotidine 5'-monophophate decarboxylase

Orotidine 5'-monophosphate decarboxylase catalyzes the decarboxylation of truncated substrate (1-β-D-erythrofuranosyl)orotic acid to form (1-β-D-erythrofuranosyl)uracil. This enzyme-catalyzed reaction is activated by tetrahedral oxydianions, which bind weakly to unliganded OMPDC and tightly to the enzyme-transition state complex, with the following intrinsic oxydianion binding energies (kcal/mol): SO3(2-), -8.3; HPO3(2-), -7.7; S2O3(2-), -4.6; SO4(2-), -4.5; HOPO3(2-), -3.0; HOAsO3(2-), no activation detected. We propose that the oxydianion and orotate binding domains of OMPDC perform complementary functions in catalysis of decarboxylation reactions: (1) The orotate binding domain carries out decarboxylation of the orotate ring. (2) The activating oxydianion binding domain has the cryptic function of utilizing binding interactions with tetrahedral inorganic oxydianions to drive an enzyme conformational change that results in the stabilization of transition states at the distant orotate domain.