The mechanism of catalysis of the chorismate to prephenate reaction by the Escherichia coli mutase enzyme

Computational redesign of cytochrome P450 CYP102A1 for highly stereoselective omeprazole hydroxylation by UniDesign

Cytochrome P450 CYP102A1 is a prototypic biocatalyst that has great potential in chemical synthesis, drug discovery, and biotechnology. CYP102A1 variants engineered by directed evolution and/or rational design are capable of catalyzing the oxidation of a wide range of organic compounds. However, it is difficult to foresee the outcome of engineering CYP102A1 for a compound of interest. Here, we introduce UniDesign as a computational framework for enzyme design and engineering. We tested UniDesign by redesigning CYP102A1 for stereoselective metabolism of omeprazole (OMP), a proton pump inhibitor, starting from an active but nonstereoselective triple mutant (TM: A82F/F87V/L188Q). To shift stereoselectivity toward (R)-OMP, we computationally scanned three active site positions (75, 264, and 328) for mutations that would stabilize the binding of the transition state of (R)-OMP while destabilizing that of (S)-OMP and picked three variants, namely UD1 (TM/L75I), UD2 (TM/A264G), and UD3 (TM/A328V), for experimentation, based on computed energy scores and models. UD1, UD2, and UD3 exhibit high turnover rates of 55 ± 4.7, 84 ± 4.8, and 79 ± 5.7 min-1, respectively, for (R)-OMP hydroxylation, whereas the corresponding rates for (S)-OMP are only 2.2 ± 0.19, 6.0 ± 0.68, and 14 ± 2.8 min-1, yielding an enantiomeric excess value of 92, 87, and 70%, respectively. These results suggest the critical roles of L75I, A264G, and A328V in steering OMP in the optimal orientation for stereoselective oxidation and demonstrate the utility of UniDesign for engineering CYP102A1 to produce drug metabolites of interest. The results are discussed in the context of protein structures.

Mechanism of the Stereoselective Catalysis of Diels-Alderase PyrE3 Involved in Pyrroindomycin Biosynthesis

The biosynthesis of pyrroindomycins A and B features a complexity-building [4 + 2] cycloaddition cascade, which generates the spirotetramate core under the catalytic effects of monofunctional Diels-Alderases PyrE3 and PyrI4. We recently showed that the main functions of PyrI4 include acid catalysis and induced-fit/conformational selection. We now present quantum mechanical and molecular dynamics studies implicating a different mode of action by PyrE3, which prearranges an anionic polyene substrate into a high-energy reactive conformation at which an inverse-electron-demand Diels-Alder reaction can occur with a low barrier. Stereoselection is realized by strong binding interactions at the endo stereochemical relationship and a local steric constraint on the endo-1,3-diene unit. These findings, illustrating distinct mechanisms for PyrE3 and PyrI4, highlight how nature has evolved multiple ways to catalyze Diels-Alder reactions.

Unraveling the mechanism of enantio-controlling switches of an alcohol dehydrogenase toward sterically small ketone

Efficient synthesis of chiral compounds under mild conditions is highly desirable in the chemical and pharmaceutical communities, but it often faces difficulties. Although various enzymes have been harnessed as biocatalysts...

Protein Dynamics and Substrate Protonation States Mediate the Catalytic Action of trans-4-Hydroxy-l-Proline Dehydratase

The enzyme trans-4-hydroxy-l-proline (Hyp) dehydratase (HypD) is among the most abundant glycyl radical enzymes (GREs) in the healthy human gut microbiome and is considered a promising antibiotic target for the prominent antibiotic-resistant pathogen Clostridium difficile. Although an enzymatic mechanism has been proposed, the role of the greater HypD protein environment in mediating radical reactivity is not well understood. To fill this gap in understanding, we investigate HypD across multiple time- and length-scales using electronic structure modeling and classical molecular dynamics. We observe that the Hyp substrate protonation state significantly alters both its enzyme-free reactivity and its dynamics within the enzyme active site. Accurate coupled-cluster modeling suggests the deprotonated form of Hyp to be the most reactive protonation state for C5-H_pro-S activation. In the protein environment, hydrophobic interactions modulate the positioning of the Cys434 radical to enhance the reactivity of C5-H_pro-S abstraction. Long-time dynamics reveal that changing Hyp protonation states triggers the switching of a Leu643-gated water tunnel, a functional feature that has not yet been observed for members of the GRE superfamily.

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.

Fine-Tuning of Sequence Specificity by Near Attack Conformations in Enzyme-Catalyzed Peptide Hydrolysis

The catalytic role of near attack conformations (NACs), molecular states that lie on the pathway between the ground state (GS) and transition state (TS) of a chemical reaction, is not understood completely. Using a computational approach that combines Bürgi–Dunitz theory with all-atom molecular dynamics simulations, the role of NACs in catalyzing the first stages of HIV-1 protease peptide hydrolysis was previously investigated using a substrate that represents the recognized SP1-NC cleavage site of the HIV-1 Gag polyprotein. NACs were found to confer no catalytic effect over the uncatalyzed reaction there ( Δ Δ G N ‡ ∼ 0 kcal/mol). Here, using the same approach, the role of NACs across multiple substrates that each represent a further recognized cleavage site is investigated. Overall rate enhancement varies by | Δ Δ G ‡ | ∼ 12–15 kcal/mol across this set, and although NACs contribute a small and approximately constant barrier to the uncatalyzed reaction (< Δ G N ‡ u > = 4.3 ± 0.3 kcal/mol), they are found to contribute little significant catalytic effect ( | Δ Δ G N ‡ | ∼ 0–2 kcal/mol). Furthermore, no correlation is exhibited between NAC contributions and the overall energy barrier ( R 2 = 0.01). However, these small differences in catalyzed NAC contributions enable rates to match those required for the kinetic order of processing. Therefore, NACs may offer an alternative and subtle mode compared to non-NAC contributions for fine-tuning reaction rates during complex evolutionary sequence selection processes—in this case across cleavable polyproteins whose constituents exhibit multiple functions during the virus life-cycle.

Crystal structure of the multifunctional SAM-dependent enzyme LepI provides insights into its catalytic mechanism

Pericyclic reactions are among the most powerful synthetic transformations widely applied in the synthesis of multiple regioselective and stereoselective carbon-carbon bonds. LepI is a recently identified S-adenosyl-l-methionine (SAM)-dependent enzyme, which could catalyze dehydration, Diels-Alder reaction, and the retro-Claisen rearrangement reactions. However, the mechanism underlying these reactions by LepI remains elusive. Here we report the structure of LepI in complex with SAM as its co-factor, which adopts a typical class I methyltransferase fold. Docking studies are performed to investigate the binding modes of various substrates/products and provide insights into the catalytic mechanism of the multiple reactions catalyzed by LepI. Our study suggests that the dehydration reaction may start from the deprotonation of the hydroxyl group on the pyridone ring of the substrate by LepI^H133, during which R295 and D296 play important roles in substrate binding and stabilizing the reaction intermediate. The stereoselective dehydration is accomplished through the trans-conformer of the leaving alcohol group which is trapped by nearby residues. The pericyclic reactions following dehydration are facilitated by the hydrophobic and hydrophilic interactions in the binding pocket. H133 and R295, two residues not conserved in other methyltransferases, might account for the unique activity of LepI among the SAM-dependent methyltransferase family. Together, this study provides important structural insights into the unique reactions catalyzed by LepI and will shed light on the knowledge of mechanisms of pericyclic reactions.

Crystal structure of LepI, a multifunctional SAM-dependent enzyme which catalyzes pericyclic reactions in leporin biosynthesis

LepI is a novel multifunctional enzyme that catalyzes stereoselective dehydration, Diels-Alder reaction, and retro-Claisen rearrangement. Here we report the crystal structure of LepI in complex with its co-factor S-adenosyl methionine (SAM). LepI forms a tetramer via the N-terminal helical domain and binds to a SAM molecule in the C-terminal catalytic domain. The binding modes of various LepI substrates are investigated by docking simulations, which suggest that the substrates are bound via both hydrophobic and hydrophilic forces, as well as cation-π interactions with the positively charged SAM. The reaction starts with a dehydration step in which H133 possibly deprotonates the pyridone hydroxyl group of the substrate, while D296 might protonate an alkyl-chain hydroxyl group. Subsequent pericyclization may be facilitated by the correct fold of the substrate's alkyl chain and a thermodynamic driving force towards σ-bonds at the expense of π-bonds. These results provide structural insights into LepI catalysis and are important in understanding the mechanism of enzymatic pericyclization.

Quantum chemical modeling of the reaction path of chorismate mutase based on the experimental substrate/product complex

Chorismate mutase is a well-known model enzyme, catalyzing the Claisen rearrangement of chorismate to prephenate. Recent high-resolution crystal structures along the reaction coordinate of this enzyme enabled computational analyses at unprecedented detail. Using quantum chemical simulations, we investigated how the catalytic reaction mechanism is affected by electrostatic and hydrogen-bond interactions. Our calculations showed that the transition state (TS) was mainly stabilized electrostatically, with Arg90 playing the leading role. The effect was augmented by selective hydrogen-bond formation to the TS in the wild-type enzyme, facilitated by a small-scale local induced fit. We further identified a previously underappreciated water molecule, which separates the negative charges during the reaction. The analysis includes the wild-type enzyme and a non-natural enzyme variant, where the catalytic arginine was replaced with an isosteric citrulline residue.

Mycobacterium tuberculosis chorismate mutase: A potential target for TB

Mycobacterium tuberculosis chorismate mutase (MtbCM) catalyzes the rearrangement of chorismate to prephenate in the shikimate biosynthetic pathway to form the essential amino acids, phenylalanine and tyrosine. Two genes encoding chorismate mutase have been identified in Mtb. The secretory form,∗MtbCM (encoded by Rv1885c) is assumed to play a key role in pathogenesis of tuberculosis. Also, the inhibition of MtbCM may hinder the supply of nutrients to the organism. Indeed, the existence of chorismate mutase (CM) in bacteria, fungi and higher plants but not in human and low sequence homology among known CM makes it an interesting target for the discovery of anti-tubercular agents. The present article mainly focuses on the recent developments in the structure, function and inhibition of MtbCM. The understanding of various aspects of MtbCM as presented in the current article may facilitate the design and subsequent chemical synthesis of new inhibitors against ∗MtbCM, that could lead to the discovery and development of novel and potent anti-tubercular agents in future.

Interdomain Conformational Changes Provide Allosteric Regulation en Route to Chorismate

Multifunctional proteins play a variety of roles in metabolism. Here, we examine the catalytic function of the combined 3-deoxy-d-arabino heptulosonate-7-phosphate synthase (DAH7PS) and chorismate mutase (CM) from Geobacillus sp. DAH7PS operates at the start of the biosynthetic pathway for aromatic metabolites, whereas CM operates in a dedicated branch of the pathway for the biosynthesis of amino acids tyrosine and phenylalanine. In line with sequence predictions, the two catalytic functions are located in distinct domains, and these two activities can be separated and retain functionality. For the full-length protein, prephenate, the product of the CM reaction, acts as an allosteric inhibitor for the DAH7PS. The crystal structure of the full-length protein with prephenate bound and the accompanying small angle x-ray scattering data reveal the molecular mechanism of the allostery. Prephenate binding results in the tighter association between the dimeric CM domains and the tetrameric DAH7PS, occluding the active site and therefore disrupting DAH7PS function. Acquisition of a physical gating mechanism to control catalytic function through gene fusion appears to be a general mechanism for providing allostery for this enzyme.

Mechanistic studies of new oximes reactivators of human butyryl cholinesterase inhibited by cyclosarin and sarin

Butyryl cholinesterase (BChE) has been seen as a key enzyme in the search for new strategies in the treatment of poisoning by organophosphates (OPs), since human BChE (HssBChE), complexed with the appropriate oxime, can be a suitable scavenger and deactivator for OPs in the blood stream. However, the efficacy of HssBChE is limited by its strict stoichiometric scavenging, slow reactivation, and propensity for aging. The improvement of the reactivation rate by new and more efficient oximes could contribute to mitigate this problem and increase the HssBChE efficiency as scavenger. Several oximes have been synthesized and tested with this goal, some with promising results, but the mechanistic aspects of the reactivation reaction are not fully understood yet. In order to better investigate this mechanism, docking and mixed quantum and molecular mechanics combined with principal components analysis were performed here to evaluate the capacity of reactivation and determine the preferred route for the reactivation reaction of two new oximes on HssBChE inhibited by the neurotoxic agents cyclosarin and sarin. Plots of potential energies were calculated and all the transition states of the reactional mechanism were determined. Our results showed a good correlation with experimental data and pointed to the most efficient oxime with both OPs. The protocol used could be a suitable tool for a preliminary evaluation of the HssBChE reactivation rates by new oximes.

Conversion of aminodeoxychorismate synthase into anthranilate synthase with Janus mutations: mechanism of pyruvate elimination catalyzed by chorismate enzymes

The central importance of chorismate enzymes in bacteria, fungi, parasites, and plants combined with their absence in mammals makes them attractive targets for antimicrobials and herbicides. Two of these enzymes, anthranilate synthase (AS) and aminodeoxychorismate synthase (ADCS), are structurally and mechanistically similar. The first catalytic step, amination at C2, is common between them, but AS additionally catalyzes pyruvate elimination, aromatizing the aminated intermediate to anthranilate. Despite prior attempts, the conversion of a pyruvate elimination-deficient enzyme into an elimination-proficient one has not been reported. Janus, a bioinformatics method for predicting mutations required to functionally interconvert homologous enzymes, was employed to predict mutations to convert ADCS into AS. A genetic selection on a library of Janus-predicted mutations was performed. Complementation of an AS-deficient strain of Escherichia coli grown on minimal medium led to several ADCS mutants that allow growth in 6 days compared to 2 days for wild-type AS. The purified mutant enzymes catalyze the conversion of chorismate to anthranilate at rates that are ∼50% of the rate of wild-type ADCS-catalyzed conversion of chorismate to aminodeoxychorismate. The residues mutated do not contact the substrate. Molecular dynamics studies suggest that pyruvate elimination is controlled by the conformation of the C2-aminated intermediate. Enzymes that catalyze elimination favor the equatorial conformation, which presents the C2-H to a conserved active site lysine (Lys424) for deprotonation and maximizes stereoelectronic activation. Acid/base catalysis of pyruvate elimination was confirmed in AS and salicylate synthase by showing incorporation of a solvent-derived proton into the pyruvate methyl group and by solvent kinetic isotope effects on pyruvate elimination catalyzed by AS.

Methyltransferases do not work by compression, cratic, or desolvation effects, but by electrostatic preorganization

The enzyme catechol O-methyltransferase (COMT) catalyzes the transfer of a methyl group from S-adenosylmethionine to dopamine and related catechols. The search for the origin of COMT catalysis has led to different proposals and hypothesis, including the entropic, the NAC, and the compression proposals as well as the more reasonable electrostatic idea. Thus, it is important to understand the catalytic power of this enzyme and to examine the validity of different proposals and in particular the repeated recent implication of the compression idea. The corresponding analysis should be done by well-defined physically-based considerations that involve computations rather than circular interpretations of experimental results. Thus, we explore here the origin of the catalytic efficiency of COMT by using the empirical valence bond and the linear response approximation approaches. The results demonstrate that the catalytic effect of COMT is mainly due to electrostatic preorganization effects. It is also shown that the compression, NAC and entropic proposals do not account for the catalytic effect.

Electrostatic transition state stabilization rather than reactant destabilization provides the chemical basis for efficient chorismate mutase catalysis

For more than half a century, transition state theory has provided a useful framework for understanding the origins of enzyme catalysis. As proposed by Pauling, enzymes accelerate chemical reactions by binding transition states tighter than substrates, thereby lowering the activation energy compared with that of the corresponding uncatalyzed process. This paradigm has been challenged for chorismate mutase (CM), a well-characterized metabolic enzyme that catalyzes the rearrangement of chorismate to prephenate. Calculations have predicted the decisive factor in CM catalysis to be ground state destabilization rather than transition state stabilization. Using X-ray crystallography, we show, in contrast, that a sluggish variant of Bacillus subtilis CM, in which a cationic active-site arginine was replaced by a neutral citrulline, is a poor catalyst even though it effectively preorganizes chorismate for the reaction. A series of high-resolution molecular snapshots of the reaction coordinate, including the apo enzyme, and complexes with substrate, transition state analog and product, demonstrate that an active site, which is only complementary in shape to a reactive substrate conformer, is insufficient for effective catalysis. Instead, as with other enzymes, electrostatic stabilization of the CM transition state appears to be crucial for achieving high reaction rates.

Substrate orientation and specificity in xanthine oxidase: crystal structures of the enzyme in complex with indole-3-acetaldehyde and guanine

Xanthine oxidase is a molybdenum-containing hydroxylase that catalyzes the hydroxylation of sp(2)-hybridized carbon centers in a variety of aromatic heterocycles as well as aldehydes. Crystal structures of the oxidase form of the bovine enzyme in complex with a poor substrate indole-3-acetaldehyde and the nonsubstrate guanine have been determined, both at a resolution of 1.6 Å. In each structure, a specific and unambiguous orientation of the substrate in the active site is observed in which the hydroxylatable site is oriented away from the active site molybdenum center. The orientation seen with indole-3-acetaldehyde has the substrate positioned with the indole ring rather than the exocyclic aldehyde nearest the molybdenum center, indicating that the substrate must rotate some 30° in the enzyme active site to permit hydroxylation of the aldehyde group (as observed experimentally), accounting for the reduced reactivity of the enzyme toward this substrate. The principal product of hydroxylation of indole-3-acetaldehyde by the bovine enzyme is confirmed to be indole-3-carboxylic acid based on its characteristic UV-vis spectrum, and the kinetics of enzyme reduction are reported. With guanine, the dominant orientation seen crystallographically has the C-8 position that might be hydroxylated pointed away from the active site molybdenum center, in a configuration resembling that seen previously with hypoxanthine (a substrate that is effectively hydroxylated at position 2). The ∼180° reorientation required to permit reaction is sterically prohibited, indicating that substrate (mis)orientation in the active site is a major factor precluding formation of the highly mutagenic 8-hydroxyguanine.

Modification of residue 42 of the active site loop with a lysine-mimetic side chain rescues isochorismate-pyruvate lyase activity in Pseudomonas aeruginosa PchB

PchB is an isochorismate-pyruvate lyase from Pseudomonas aeruginosa. A positively charged lysine residue is located in a flexible loop that behaves as a lid to the active site, and the lysine residue is required for efficient production of salicylate. A variant of PchB that lacks the lysine at residue 42 has a reduced catalytic free energy of activation of up to 4.4 kcal/mol. Construction of a lysine isosteric residue bearing a positive charge at the appropriate position leads to the recovery of 2.5-2.7 kcal/mol (about 60%) of the 4.4 kcal/mol by chemical rescue. Exogenous addition of ethylamine to the K42A variant leads to a neglible recovery of activity (0.180 kcal/mol, roughly 7% rescue), whereas addition of propylamine caused an additional modest loss in catalytic power (0.056 kcal/mol, or 2% loss). This is consistent with the view that (a) the lysine-42 residue is required in a specific conformation to stabilize the transition state and (b) the correct conformation is achieved for a lysine-mimetic side chain at site 42 in the course of loop closure, as expected for transition-state stabilization by the side chain ammonio function. That the positive charge is the main effector of transition state stabilization is shown by the construction of a lysine-isosteric residue capable of exerting steric effects and hydrogen bonding but not electrostatic effects, leading to a modest increase of catalytic power (0.267-0.505 kcal/mol of catalytic free energy, or roughly 6-11% rescue).

Pericyclic reactions catalyzed by chorismate-utilizing enzymes

One of the fundamental questions of enzymology is how catalytic power is derived. This review focuses on recent developments in the structure--function relationships of chorismate-utilizing enzymes involved in siderophore biosynthesis to provide insight into the biocatalysis of pericyclic reactions. Specifically, salicylate synthesis by the two-enzyme pathway in Pseudomonas aeruginosa is examined. The isochorismate-pyruvate lyase is discussed in the context of its homologues, the chorismate mutases, and the isochorismate synthase is compared to its homologues in the MST family (menaquinone, siderophore, or tryptophan biosynthesis) of enzymes. The tentative conclusion is that the activities observed cannot be reconciled by inspection of the active site participants alone. Instead, individual activities must arise from unique dynamic properties of each enzyme that are tuned to promote specific chemistries.

pH Dependence of catalysis by Pseudomonas aeruginosa isochorismate-pyruvate lyase: implications for transition state stabilization and the role of lysine 42

An isochorismate-pyruvate lyase with adventitious chorismate mutase activity from Pseudomonas aerugionsa (PchB) achieves catalysis of both pericyclic reactions in part by the stabilization of reactive conformations and in part by electrostatic transition-state stabilization. When the active site loop Lys42 is mutated to histidine, the enzyme develops a pH dependence corresponding to a loss of catalytic power upon deprotonation of the histidine. Structural data indicate that the change is not due to changes in active site architecture, but due to the difference in charge at this key site. With loss of the positive charge on the K42H side chain at high pH, the enzyme retains lyase activity at ∼100-fold lowered catalytic efficiency but loses detectable mutase activity. We propose that both substrate organization and electrostatic transition state stabilization contribute to catalysis. However, the dominant reaction path for catalysis is dependent on reaction conditions, which influence the electrostatic properties of the enzyme active site amino acid side chains.

Entropic and enthalpic components of catalysis in the mutase and lyase activities of Pseudomonas aeruginosa PchB

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions, demonstrating the eponymous activity and also chorismate mutase activity. The thermodynamic parameters for these enzyme-catalyzed activities, as well as the uncatalyzed isochorismate decomposition, are reported from temperature dependence of k(cat) and k(uncat) data. The entropic effects do not contribute to enzyme catalysis as expected from previously reported chorismate mutase data. Indeed, an entropic penalty for the enzyme-catalyzed mutase reaction (ΔS(++) = -12.1 ± 0.6 cal/(mol K)) is comparable to that of the previously reported uncatalyzed reaction, whereas that of the enzyme-catalyzed lyase reaction (ΔS(++) = -24.3 ± 0.2 cal/(mol K)) is larger than that of the uncatalyzed lyase reaction (-15.77 ± 0.02 cal/(mol K)) documented here. With the assumption that chemistry is rate-limiting, we propose that a reactive substrate conformation is formed upon loop closure of the active site and that ordering of the loop contributes to the entropic penalty for converting the enzyme substrate complex to the transition state.

Examining the case for the effect of barrier compression on tunneling, vibrationally enhanced catalysis, catalytic entropy and related issues

The idea that tunneling is enhanced by the compression of the donor-acceptor distance has attracted significant interest. In particular, recent studies argued that this proposal is consistent with pressure effects on enzymatic reactions, and that the observed pressure effects support the idea of vibrationally enhanced catalysis. However, a careful analysis of the current works reveals serious inconsistencies in the evidence presented to support these hypotheses. Apparently, tunneling decreases upon compression, and external pressure does not lead to the applicable compression of the free energy surface. Additionally, pressure experiments do not provide actual evidence for vibrationally enhanced catalysis. Finally, the temperature dependence of the entropy change in hydride transfer reactions is shown to reflect simple electrostatic effects.

Structure-function analyses of isochorismate-pyruvate lyase from Pseudomonas aeruginosa suggest differing catalytic mechanisms for the two pericyclic reactions of this bifunctional enzyme

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions in a single active site. PchB physiologically produces salicylate and pyruvate from isochorismate for ultimate incorporation of the salicylate into the siderophore pyochelin. PchB also produces prephenate from chorismate, most likely due to structural homology to the Escherchia coli chorismate mutase. The molecular basis of catalysis among enzymatic pericyclic reactions is a matter of debate, one view holding that catalysis may be derived from electrostatic transition state stabilization and the opposing view that catalysis is derived from the generation of a reactive substrate conformation. Mutant forms of PchB were generated by site-directed mutagenesis at the site (K42) hypothesized to be key for electrostatic transition state stabilization (K42A, K42Q, K42E, and K42H). The loop containing K42 is mobile, and a mutant to slow loop dynamics was also designed (A43P). Finally, a previously characterized mutation (I87T) was also produced. Circular dichroism was used to assess the overall effect on secondary structure as a result of the mutations, and X-ray crystallographic structures are reported for K42A with salicylate and pyruvate bound and for apo-I87T. The data illustrate that the active site architecture is maintained in K42A-PchB, which indicates that differences in activity are not caused by secondary structural changes or by differences in active site loop conformation but rather by the chemical nature of this key residue. In contrast, the I87T structure demonstrates considerable mobility, suggesting that loop dynamics and conformational plasticity may be important for efficient catalysis. Finally, the mutational effects on k(cat) provide evidence that the two activities of PchB are not covariant and that a single hypothesis may not provide a sufficient explanation for catalysis.

Tunneling and Classical Paths for Proton Transfer in an Enzyme Reaction Dominated by Tunneling: Oxidation of Tryptamine by Aromatic Amine Dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.

Design, synthesis, and evaluation of aza inhibitors of chorismate mutase

A series of aza inhibitors (4-9) of chorismate mutase (E.C. 5.4.99.5) was designed, prepared, and evaluated against the enzyme by monitoring the direct inhibition of the chorismate, 1, to prephenate, 2, conversion. None of these aza inhibitors displayed tighter binding to the enzyme than the native substrate chorismate or greater inhibitory action than the previously reported ether analogue, 3. Furthermore, no time-dependent loss of enzyme activity was observed in the presence of the two potentially reactive aza inhibitors (7 and 9). These results in conjunction with inhibition data from a broader series of chorismate mutase inhibitors allowed a novel proposal for the mechanistic role of chorismate mutase to be developed. This proposed mechanism was computationally verified and correlated with crystallographic studies of various chorismate mutases.

Computational Study of the Ground State of Thermophilic Indole Glycerol Phosphate Synthase: Structural Alterations at the Active Site with Temperature

Hyperthermophlic indole-3-glycerol phosphate synthase (IGPS) catalyzes the terminal ring-closure step in tryptophan biosynthesis. In this paper, we compare the results from the molecular dynamics (MD) simulation of enzyme-bound substrate at 298 K (E.S298) and 385 K (E.S385) solvated in the TIP3P water box using the CHARMM force field to address the question of the structural change of the Enzyme. Substrate complex with temperature. The population of the reactive Enzyme. Substrate conformers (near attack conformers or NACs) increases by approximately 1100-fold in going from room temperature (E.S298) to high temperature (E.S385). This increased population of NAC conformers in the Michaelis complex correlates well with the increase in rate in going from 298 to 385 K. The positioning of the two active site residues Lys53 and Lys110 controls binding of the substrate in the favorable orientation for general acid-catalyzed intramolecular ring formation reaction. It can be concluded that the NAC formation allowing general acid catalysis has much to do with the temperature dependence of the free energy of reaction.

The near attack conformation approach to the study of the chorismate to prephenate reaction

Standard free energies (DeltaGN degree) for formation of near attack conformers, those ground state conformers that can convert directly to the transition state, were calculated for the Claisen rearrangement of chorismate to prephenate in six different environments: water, wild-type enzymes from Bacillus subtilis and Escherichia coli, their Arg90Cit and Glu52Ala mutants, and the 1F7 catalytic antibody. Values of the calculated DeltaGN degrees and the experimentally determined activation energies (DeltaG++) are linearly related with the slope of approximately equal to 1. This demonstrates that the relative rate of the chorismate --> prephenate reaction is overwhelmingly dependent on the efficiency of formation of near attack conformers in the ground state.

A chorismate mutase from the soybean cyst nematode Heterodera glycines shows polymorphisms that correlate with virulence

Parasitism genes from phytoparasitic nematodes are thought to be essential for nematode invasion of the host plant, to help the nematode establish feeding sites, and to aid nematodes in the suppression of host plant defenses. One gene that may play several roles in nematode parasitism is chorismate mutase (CM). This secreted enzyme is produced in the nematode's esophageal glands and appears to function within the plant cell to manipulate the plant's shikimate pathway, which controls plant cell growth, development, structure, and pathogen defense. Using degenerate polymerase chain reaction primers, we amplified and cloned a chorismate mutase (Hg-cm-1) from Heterodera glycines, the soybean cyst nematode (SCN), and showed it had CM activity. RNA in situ hybridization of Hg-cm-1 cDNA to SCN sections confirms that it is specifically expressed in the nematodes' esophageal glands. DNA gel blots of genomic DNA isolated from SCN inbred lines that have differing virulence on SCN resistant soybean show Hg-cm-1 is a member of a polymorphic gene family. Some Hg-cm family members predominate in SCN inbred lines that are virulent on certain SCN resistant soybean cultivars. The same polymorphisms and correlation with virulence are seen in the Hg-cm-1 expressed in the SCN second-stage juveniles. Based on the enzymatic activity of Hg-cm-1 and the observation that different forms of the mutase are expressed in virulent nematodes, we hypothesize that the Hg-cm-1 is a virulence gene, some forms of which allow SCN to parasitize certain resistant soybean plants.