Conformational changes in orotidine 5'-monophosphate decarboxylase: a structure-based explanation for how the 5'-phosphate group activates the enzyme

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.

Estimating Electron Density Support for Individual Atoms and Molecular Fragments in X-ray Structures

Macromolecular structures resolved by X-ray crystallography are essential for life science research. While some methods exist to automatically quantify the quality of the electron density fit, none of them is without flaws. Especially the question of how well individual parts like atoms, small fragments, or molecules are supported by electron density is difficult to quantify. While taking experimental uncertainties correctly into account, they do not offer an answer on how reliable an individual atom position is. A rapid quantification of this atomic position reliability would be highly valuable in structure-based molecular design. To overcome this limitation, we introduce the electron density score EDIA for individual atoms and molecular fragments. EDIA assesses rapidly, automatically, and intuitively the fit of individual as well as multiple atoms (EDIA_m) into electron density accompanied by an integrated error analysis. The computation is based on the standard 2fo - fc electron density map in combination with the model of the molecular structure. For evaluating partial structures, EDIA_m shows significant advantages compared to the real-space R correlation coefficient (RSCC) and the real-space difference density Z score (RSZD) from the molecular modeler's point of view. Thus, EDIA abolishes the time-consuming step of visually inspecting the electron density during structure selection and curation. It supports daily modeling tasks of medicinal and computational chemists and enables a fully automated assembly of large-scale, high-quality structure data sets. Furthermore, EDIA scores can be applied for model validation and method development in computer-aided molecular design. In contrast to measuring the deviation from the structure model by root-mean-squared deviation, EDIA scores allow comparison to the underlying experimental data taking its uncertainty into account.

Utilization of Substrate Intrinsic Binding Energy for Conformational Change and Catalytic Function in Phosphoenolpyruvate Carboxykinase

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Previous work has demonstrated that the enzyme cycles between a catalytically inactive open state and a catalytically active closed state. The transition of the enzyme between these states requires the transition of several active site loops to shift from mobile, disordered structural elements to stable ordered states. The mechanism by which these disorder-order transitions are coupled to the ligation state of the active site however is not fully understood. To further investigate the mechanisms by which the mobility of the active site loops is coupled to enzymatic function and the transitioning of the enzyme between the two conformational states, we have conducted structural and functional studies of point mutants of E89. E89 is a proposed key member of the interaction network of mobile elements as it resides in the R-loop region of the enzyme active site. These new data demonstrate the importance of the R-loop in coordinating interactions between substrates at the OAA/PEP binding site and the mobile R- and Ω-loop domains. In turn, the studies more generally demonstrate the mechanisms by which the intrinsic ligand binding energy can be utilized in catalysis to drive unfavorable conformational changes, changes that are subsequently required for both optimal catalytic activity and fidelity.

Rate and Equilibrium Constants for an Enzyme Conformational Change during Catalysis by Orotidine 5'-Monophosphate Decarboxylase

The caged complex between orotidine 5'-monophosphate decarboxylase (ScOMPDC) and 5-fluoroorotidine 5'-monophosphate (FOMP) undergoes decarboxylation ∼300 times faster than the caged complex between ScOMPDC and the physiological substrate, orotidine 5'-monophosphate (OMP). Consequently, the enzyme conformational changes required to lock FOMP at a protein cage and release product 5-fluorouridine 5'-monophosphate (FUMP) are kinetically significant steps. The caged form of ScOMPDC is stabilized by interactions between the side chains from Gln215, Tyr217, and Arg235 and the substrate phosphodianion. The control of these interactions over the barrier to the binding of FOMP and the release of FUMP was probed by determining the effect of all combinations of single, double, and triple Q215A, Y217F, and R235A mutations on kcat/Km and kcat for turnover of FOMP by wild-type ScOMPDC; its values are limited by the rates of substrate binding and product release, respectively. The Q215A and Y217F mutations each result in an increase in kcat and a decrease in kcat/Km, due to a weakening of the protein-phosphodianion interactions that favor fast product release and slow substrate binding. The Q215A/R235A mutation causes a large decrease in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of OMP, which are limited by the rate of the decarboxylation step, but much smaller decreases in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of FOMP, which are limited by the rate of enzyme conformational changes. By contrast, the Y217A mutation results in large decreases in kcat/Km for ScOMPDC-catalyzed decarboxylation of both OMP and FOMP, because of the comparable effects of this mutation on rate-determining decarboxylation of enzyme-bound OMP and on the rate-determining enzyme conformational change for decarboxylation of FOMP. We propose that kcat = 8.2 s(-1) for decarboxylation of FOMP by the Y217A mutant is equal to the rate constant for cage formation from the complex between FOMP and the open enzyme, that the tyrosyl phenol group stabilizes the closed form of ScOMPDC by hydrogen bonding to the substrate phosphodianion, and that the phenyl group of Y217 and F217 facilitates formation of the transition state for the rate-limiting conformational change. An analysis of kinetic data for mutant enzyme-catalyzed decarboxylation of OMP and FOMP provides estimates for the rate and equilibrium constants for the conformational change that traps FOMP at the enzyme active site.

Orotidine Monophosphate Decarboxylase--A Fascinating Workhorse Enzyme with Therapeutic Potential

Orotidine 5'-monophosphate decarboxylase (ODCase) is known as one of the most proficient enzymes. The enzyme catalyzes the last reaction step of the de novo pyrimidine biosynthesis, the conversion from orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate. The enzyme is found in all three domains of life, Bacteria, Eukarya and Archaea. Multiple sequence alignment of 750 putative ODCase sequences resulted in five distinct groups. While the universally conserved DxKxxDx motif is present in all the groups, depending on the groups, several characteristic motifs and residues can be identified. Over 200 crystal structures of ODCases have been determined so far. The structures, together with biochemical assays and computational studies, elucidated that ODCase utilized both transition state stabilization and substrate distortion to accelerate the decarboxylation of its natural substrate. Stabilization of the vinyl anion intermediate by a conserved lysine residue at the catalytic site is considered the largest contributing factor to catalysis, while bending of the carboxyl group from the plane of the aromatic pyrimidine ring of OMP accounts for substrate distortion. A number of crystal structures of ODCases complexed with potential drug candidate molecules have also been determined, including with 6-iodo-uridine, a potential antimalarial agent.

Investigating the role of a backbone to substrate hydrogen bond in OMP decarboxylase using a site-specific amide to ester substitution

Hydrogen bonds between backbone amide groups of enzymes and their substrates are often observed, but their importance in substrate binding and/or catalysis is not easy to investigate experimentally. We describe the generation and kinetic characterization of a backbone amide to ester substitution in the orotidine 5'-monophosphate (OMP) decarboxylase from Methanobacter thermoautotrophicum (MtOMPDC) to determine the importance of a backbone amide-substrate hydrogen bond. The MtOMPDC-catalyzed reaction is characterized by a rate enhancement (∼10(17)) that is among the largest for enzyme-catalyzed reactions. The reaction proceeds through a vinyl anion intermediate that may be stabilized by hydrogen bonding interaction between the backbone amide of a conserved active site serine residue (Ser-127) and oxygen (O4) of the pyrimidine moiety and/or electrostatic interactions with the conserved general acidic lysine (Lys-72). In vitro translation in conjunction with amber suppression using an orthogonal amber tRNA charged with L-glycerate ((HO)S) was used to generate the ester backbone substitution (S127(HO)S). With 5-fluoro OMP (FOMP) as substrate, the amide to ester substitution increased the value of Km by ∼1.5-fold and decreased the value of kcat by ∼50-fold. We conclude that (i) the hydrogen bond between the backbone amide of Ser-127 and O4 of the pyrimidine moiety contributes a modest factor (∼10(2)) to the 10(17) rate enhancement and (ii) the stabilization of the anionic intermediate is accomplished by electrostatic interactions, including its proximity of Lys-72. These conclusions are in good agreement with predictions obtained from hybrid quantum mechanical/molecular mechanical calculations.

Substrate distortion contributes to the catalysis of orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate (UMP) by 17 orders of magnitude. Eight new crystal structures with ligand analogues combined with computational analyses of the enzyme's short-lived intermediates and the intrinsic electronic energies to distort the substrate and other ligands improve our understanding of the still controversially discussed reaction mechanism. In their respective complexes, 6-methyl-UMP displays significant distortion of its methyl substituent bond, 6-amino-UMP shows the competition between the K72 and C6 substituents for a position close to D70, and the methyl and ethyl esters of OMP both induce rotation of the carboxylate group substituent out of the plane of the pyrimidine ring. Molecular dynamics and quantum mechanics/molecular mechanics computations of the enzyme-substrate complex also show the bond between the carboxylate group and the pyrimidine ring to be distorted, with the distortion contributing a 10-15% decrease of the ΔΔG(⧧) value. These results are consistent with ODCase using both substrate distortion and transition-state stabilization, primarily exerted by K72, in its catalysis of the OMP decarboxylation reaction.

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

Orotidine 5'-monophosphate (OMP) decarboxylase (ODCase) catalyzes the decarboxylation of OMP to uridine 5'-monophosphate (UMP). Numerous studies of this reaction have suggested a plethora of mechanisms including covalent addition, ylide or carbene formation, and concerted or stepwise protonation. Recent experiments and simulations present strong evidence for a direct decarboxylation mechanism, although direct comparison between experiment and theory is still lacking. In the current work we present hybrid quantum mechanics-molecular mechanics simulations that address the detailed decarboxylation mechanisms for OMP and 5-fluoro-OMP by ODCase. Multidimensional potentials of mean force are computed as functions of structural progress coordinates for the Methanobacterium thermoautotrophicum ODCase reaction: the decarboxylation reaction coordinate, an orbital rehybridization coordinate, and the proton transfer coordinate between Lys72 and the substrate. The computed free energy profiles are in accord with the available experimental data. To facilitate further direct comparison with experiment, we compute the kinetic isotope effects (KIEs) for the enzyme-catalyzed reactions using a mass-perturbation-based path-integral method. The computed KIE provide further support for a direct decarboxylation mechanism. In agreement with experiment, the data suggest a role for Lys72 in stabilizing the transition state in the catalysis of OMP and, to a somewhat lesser extent, in 5-fluoro-OMP.

A mutational analysis of the active site loop residues in cis-3-Chloroacrylic acid dehalogenase

cis-3-Chloroacrylic acid dehalogenase (cis-CaaD) from Pseudomonas pavonaceae 170 and a homologue from Corynebacterium glutamicum designated Cg10062 are 34% identical in sequence (54% similar). The former catalyzes a key step in a bacterial catabolic pathway for the nematocide 1,3-dichloropropene, whereas the latter has no known biological activity. Although Cg10062 has the six active site residues (Pro-1, His-28, Arg-70, Arg-73, Tyr-103, and Glu-114) that are critical for cis-CaaD activity, it shows only a low level cis-CaaD activity and lacks the specificity of cis-CaaD: Cg10062 processes both isomers of 3-chloroacrylate with a preference for the cis isomer. The basis for these differences is unknown, but a comparison of the crystal structures of the enzymes covalently modified by an adduct resulting from their incubation with the same inhibitor offers a possible explanation. A six-residue active site loop in cis-CaaD shows a conformation strikingly different from that observed in Cg10062: the loop closes down on the active site of cis-CaaD, but not on that of Cg10062. To examine what this loop might contribute to cis-CaaD catalysis and specificity, the residues were changed individually to those found in Cg10062. Subsequent kinetic and mechanistic analysis suggests that the T34A mutant of cis-CaaD is more Cg10062-like. The mutant enzyme shows a 4-fold increase in Km (using cis-3-bromoacrylate), but not to the degree observed for Cg10062 (687-fold). The mutation also causes a 4-fold decrease in the burst rate (compared to that of wild-type cis-CaaD), whereas Cg10062 shows no burst rate. More telling is the reaction of the T34A mutant of cis-CaaD with the alternate substrate, 2,3-butadienoate. In the presence of NaBH4 and the allene, cis-CaaD is completely inactivated after one turnover because of the covalent modification of Pro-1. The same experiment with Cg10062 does not result in the covalent modification of Pro-1. The different outcomes are attributed to covalent catalysis (using Pro-1) followed by hydrolysis of the enamine or imine tautomer in cis-CaaD versus direct hydration of the allene to yield acetoacetate in the case of Cg10062. The T34A mutant shows partial inactivation, requiring five turnovers of the substrate per monomer, which suggests that the direct hydration route is favored 80% of the time. However, the mutation does not alter the stereochemistry at C-2 of [2-D]acetoacetate when the reaction is conducted in D2O. Both cis-CaaD and the T34 mutant generate (2R)-[2-D]acetoacetate, whereas Cg10062 generates mostly the 2S isomer. The combined observations are consistent with a role for the loop region in cis-CaaD specificity and catalysis, but the precise role remains to be determined.

Atomic resolution structure of the orotidine 5'-monophosphate decarboxylase product complex combined with surface plasmon resonance analysis: implications for the catalytic mechanism

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of its substrate by 17 orders of magnitude. One argument brought forward against steric/electrostatic repulsion causing substrate distortion at the carboxylate substituent as part of the catalysis has been the weak binding affinity of the decarboxylated product (UMP). The crystal structure of the UMP complex of ODCase at atomic resolution (1.03 Å) shows steric competition between the product UMP and the side chain of a catalytic lysine residue. Surface plasmon resonance analysis indicates that UMP binds 5 orders of magnitude more tightly to a mutant in which the interfering side chain has been removed than to wild-type ODCase. These results explain the low affinity of UMP and counter a seemingly very strong argument against a contribution of substrate distortion to the catalytic reaction mechanism of ODCase.

Long-range interactions in the α subunit of tryptophan synthase help to coordinate ligand binding, catalysis, and substrate channeling

The α-subunit of tryptophan synthase (αTS) catalyzes the conversion of indole-3-glycerol phosphate to d-glyceraldehyde-3-phosphate and indole. We propose that allosteric networks intrinsic to αTS are modulated by the binding of the β-subunit to regulate αTS function. Understanding these long-range amino acid networks in αTS thus gives insight into the coordination of the two active sites within TS. In this study, we have used Ala residues as probes for structural and dynamic changes of αTS throughout its catalytic cycle, in the absence of the β-subunit. Projection analysis of the chemical shift changes by site-specific amino acid substitutions and ligand titrations indicates that αTS has three important conformational states: ligand-free, glyceraldehyde-3-phosphate-bound(like), and the active states. The amino acid networks within these conformations are different, as suggested by chemical shift correlation analysis. In particular, there are long-range connections, only in the active state, between Ala47, which reports on structural and dynamic changes associated with the general acid/base Glu49, and residues within the β2α2 loop, which contains the catalytically important Asp60 residue. These long-range interactions are likely important for coordinating chemical catalysis. In the free state, but not in the active state, there are connections between the β2α2 and β6α6 loops that likely help to coordinate substrate binding. Changes in the allosteric networks are also accompanied by protein dynamic changes. During catalytic turnover, the protein becomes more rigid on the millisecond timescale and the active-site dynamics are driven to a faster nanosecond timescale.

Catalysis by orotidine 5'-monophosphate decarboxylase: effect of 5-fluoro and 4'-substituents on the decarboxylation of two-part substrates

The syntheses of two novel truncated analogs of the natural substrate orotidine 5'-monophosphate (OMP) for orotidine 5'-monophosphate decarboxylase (OMPDC) with enhanced reactivity toward decarboxylation are reported: 1-(β-d-erythrofuranosyl)-5-fluoroorotic acid (FEO) and 5'-deoxy-5-fluoroorotidine (5'-dFO). A comparison of the second-order rate constants for the OMPDC-catalyzed decarboxylations of FEO (10 M⁻¹ s⁻¹) and 1-(β-d-erythrofuranosyl)orotic acid (EO, 0.026 M⁻¹ s⁻¹) shows that the vinyl carbanion-like transition state is stabilized by 3.5 kcal/mol by interactions with the 5-F substituent of FEO. The OMPDC-catalyzed decarboxylations of FEO and EO are both activated by exogenous phosphite dianion (HPO₃²⁻), but the 5-F substituent results in only a 0.8 kcal stabilization of the transition state for the phosphite-activated reaction of FEO. This provides strong evidence that the phosphite-activated OMPDC-catalyzed reaction of FEO is not limited by the chemical step of decarboxylation of the enzyme-bound substrate. Evidence is presented that there is a change in the rate-limiting step from the chemical step of decarboxylation for the phosphite-activated reaction of EO, to closure of the phosphate gripper loop and an enzyme conformational change at the ternary E•FEO•HPO₃²⁻ complex for the reaction of FEO. The 4'-CH₃ and 4'-CH₂OH groups of 5'-dFO and orotidine, respectively, result in identical destabilizations of the transition state for the unactivated decarboxylation of 2.9 kcal/mol. By contrast, the 4'-CH₃ group of 5'-dFO and the 4'-CH₂OH group of orotidine result in very different 4.7 and 8.3 kcal/mol destabilizations of the transition state for the phosphite-activated decarboxylation. Here, the destabilizing effect of the 4'-CH₃ substituent at 5'-dFO is masked by the rate-limiting conformational change that depresses the third-order rate constant for the phosphite-activated reaction of the parent substrate FEO.

The Ω-loop lid domain of phosphoenolpyruvate carboxykinase is essential for catalytic function

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Recent studies have demonstrated that the enzyme contains a mobile active site lid domain that undergoes a transition between an open, disorded conformation and a closed, ordered conformation as the enzyme progresses through the catalytic cycle. The understanding of how this mobile domain functions in catalysis is incomplete. Previous studies showed that the closure of the lid domain stabilizes the reaction intermediate and protects the reactive intermediate from spurious protonation and thus contributes to the fidelity of the enzyme. To more fully investigate the roles of the lid domain in PEPCK function, we introduced three mutations that replaced the 11-residue lid domain with one, two, and three glycine residues. Kinetic analysis of the mutant enzymes demonstrates that none of the enzyme constructs exhibit any measurable kinetic activity, resulting in a decrease in the catalytic parameters of at least 10(6). Structural characterization of the mutants in complexes representing the catalytic cycle suggests that the inactivity is due to a role for the lid domain in the formation of the fully closed state of the enzyme that is required for catalytic function. In the absence of the lid domain, the enzyme is unable to achieve the fully closed state and is rendered inactive despite possessing all of the residues and substrates required for catalytic function. This work demonstrates how enzyme catalytic function can be abolished through the alteration of conformational equilibria despite all the elements required for chemical conversion of substrates to products remaining intact.