Enzymatic transition-state analysis and transition-state analogs

Enzymatic Synthesis of 2-Chloropurine Arabinonucleosides with Chiral Amino Acid Amides at the C6 Position and an Evaluation of Antiproliferative Activity In Vitro

A number of purine arabinosides containing chiral amino acid amides at the C6 position of the purine were synthesized using a transglycosylation reaction with recombinant E. coli nucleoside phosphorylases. Arsenolysis of 2-chloropurine ribosides with chiral amino acid amides at C6 was used for the enzymatic synthesis, and the reaction equilibrium shifted towards the synthesis of arabinonucleosides. The synthesized nucleosides were shown to be resistant to the action of E. coli adenosine deaminase. The antiproliferative activity of the synthesized nucleosides was studied on human acute myeloid leukemia cell line U937. Among all the compounds, the serine derivative exhibited an activity level (IC50 = 16 μM) close to that of Nelarabine (IC50 = 3 μM) and was evaluated as active.

Clostridioides difficile TcdB Toxin Glucosylates Rho GTPase by an SNi Mechanism and Ion Pair Transition State

Toxins TcdA and TcdB from Clostridioides difficile glucosylate human colon Rho GTPases. TcdA and TcdB glucosylation of RhoGTPases results in cytoskeletal changes, causing cell rounding and loss of intestinal integrity. Clostridial toxins TcdA and TcdB are proposed to catalyze glucosylation of Rho GTPases with retention of stereochemistry from UDP-glucose. We used kinetic isotope effects to analyze the mechanisms and transition-state structures of the glucohydrolase and glucosyltransferase activities of TcdB. TcdB catalyzes Rho GTPase glucosylation with retention of stereochemistry, while hydrolysis of UDP-glucose by TcdB causes inversion of stereochemistry. Kinetic analysis revealed TcdB glucosylation via the formation of a ternary complex with no intermediate, supporting an S_Ni mechanism with nucleophilic attack and leaving group departure occurring on the same face of the glucose ring. Kinetic isotope effects combined with quantum mechanical calculations revealed that the transition states of both glucohydrolase and glucosyltransferase activities of TcdB are highly dissociative. Specifically, the TcdB glucosyltransferase reaction proceeds via an S_Ni mechanism with the formation of a distinct oxocarbenium phosphate ion pair transition state where the glycosidic bond to the UDP leaving group breaks prior to attack of the threonine nucleophile from Rho GTPase.

Transition Path Sampling Based Calculations of Free Energies for Enzymatic Reactions: The Case of Human Methionine Adenosyl Transferase and Plasmodium vivax Adenosine Deaminase

Transition path sampling (TPS) is widely used for the calculations of reaction rates, transition state structures, and reaction coordinates of condensed phase systems. Here we discuss a scheme for the calculation of free energies using the ensemble of TPS reactive trajectories in combination with a window-based sampling technique for enzyme-catalyzed reactions. We calculate the free energy profiles of the reactions catalyzed by the human methionine S-adenosyltransferase (MAT2A) enzyme and the Plasmodium vivax adenosine deaminase (pvADA) enzyme to assess the accuracy of this method. MAT2A catalyzes the formation of S-adenosine-l-methionine following a S_N2 mechanism, and using our method, we estimate the free energy barrier for this reaction to be 16 kcal mol^-1, which is in excellent agreement with the experimentally measured activation energy of 17.27 kcal mol^-1. Furthermore, for the pvADA enzyme-catalyzed reaction we estimate a free energy barrier of 21 kcal mol^-1, and the calculated free energy profile is similar to that predicted from experimental observations. Calculating free energies by employing our simple method within TPS provides significant advantages over methods such as umbrella sampling because it is free from any applied external bias, is accurate compared to experimental measurements, and has a reasonable computational cost.

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions

We present a fast, accurate, and robust approach for determination of free energy profiles and kinetic isotope effects for RNA 2'-O-transphosphorylation reactions with inclusion of nuclear quantum effects. We apply a deep potential range correction (DPRc) for combined quantum mechanical/molecular mechanical (QM/MM) simulations of reactions in the condensed phase. The method uses the second-order density-functional tight-binding method (DFTB2) as a fast, approximate base QM model. The DPRc model modifies the DFTB2 QM interactions and applies short-range corrections to the QM/MM interactions to reproduce ab initio DFT (PBE0/6-31G*) QM/MM energies and forces. The DPRc thus enables both QM and QM/MM interactions to be tuned to high accuracy, and the QM/MM corrections are designed to smoothly vanish at a specified cutoff boundary (6 Å in the present work). The computational speed-up afforded by the QM/MM+DPRc model enables free energy profiles to be calculated that include rigorous long-range QM/MM interactions under periodic boundary conditions and nuclear quantum effects through a path integral approach using a new interface between the AMBER and i-PI software. The approach is demonstrated through the calculation of free energy profiles of a native RNA cleavage model reaction and reactions involving thio-substitutions, which are important experimental probes of the mechanism. The DFTB2+DPRc QM/MM free energy surfaces agree very closely with the PBE0/6-31G* QM/MM results, and it is vastly superior to the DFTB2 QM/MM surfaces with and without weighted thermodynamic perturbation corrections. ¹⁸O and ³⁴S primary kinetic isotope effects are compared, and the influence of nuclear quantum effects on the free energy profiles is examined.

Inhibition of Clostridium difficile TcdA and TcdB toxins with transition state analogues

Clostridium difficile causes life-threatening diarrhea and is the leading cause of healthcare-associated bacterial infections in the United States. TcdA and TcdB bacterial toxins are primary determinants of disease pathogenesis and are attractive therapeutic targets. TcdA and TcdB contain domains that use UDP-glucose to glucosylate and inactivate host Rho GTPases, resulting in cytoskeletal changes causing cell rounding and loss of intestinal integrity. Transition state analysis revealed glucocationic character for the TcdA and TcdB transition states. We identified transition state analogue inhibitors and characterized them by kinetic, thermodynamic and structural analysis. Iminosugars, isofagomine and noeuromycin mimic the transition state and inhibit both TcdA and TcdB by forming ternary complexes with Tcd and UDP, a product of the TcdA- and TcdB-catalyzed reactions. Both iminosugars prevent TcdA- and TcdB-induced cytotoxicity in cultured mammalian cells by preventing glucosylation of Rho GTPases. Iminosugar transition state analogues of the Tcd toxins show potential as therapeutics for C. difficile pathology.

The Clostridium difficile virulence factors TcdA and TcdB contain a glucosyltransferase domain (GTD), which has both glucohydrolase (GH) and glucosyltransferase (GT) activities. Here, the authors characterize the transition state features of the TcdA and TcdB GH reactions by measuring kinetic isotope effects and they identify two transition state analogues, isofagomine and noeuromycin that inhibit TcdA and TcdB. They also present the crystal structures of TcdB-GTD bound to these inhibitors and the reaction product UDP.

Enzymatic Transition States and Drug Design

Transition state theory teaches that chemically stable mimics of enzymatic transition states will bind tightly to their cognate enzymes. Kinetic isotope effects combined with computational quantum chemistry provides enzymatic transition state information with sufficient fidelity to design transition state analogues. Examples are selected from various stages of drug development to demonstrate the application of transition state theory, inhibitor design, physicochemical characterization of transition state analogues, and their progress in drug development.

Transition State Analogue Inhibitors of 5'-Deoxyadenosine/5'-Methylthioadenosine Nucleosidase from Mycobacterium tuberculosis

Mycobacterium tuberculosis 5'-deoxyadenosine/5'-methylthioadenosine nucleosidase (Rv0091) catalyzes the N-riboside hydrolysis of its substrates 5'-methylthioadenosine (MTA) and 5'-deoxyadenosine (5'-dAdo). 5'-dAdo is the preferred substrate, a product of radical S-adenosylmethionine-dependent enzyme reactions. Rv0091 is characterized by a ribocation-like transition state, with low N-ribosidic bond order, an N7-protonated adenine leaving group, and an activated but weakly bonded water nucleophile. DADMe-Immucillins incorporating 5'-substituents of the substrates 5'-dAdo and MTA were synthesized and characterized as inhibitors of Rv0091. 5'-Deoxy-DADMe-Immucillin-A was the most potent among the 5'-dAdo transition state analogues with a dissociation constant of 640 pM. Among the 5'-thio substituents, hexylthio-DADMe-Immucillin-A was the best inhibitor at 87 pM. The specificity of Rv0091 for the Immucillin transition state analogues differs from those of other bacterial homologues because of an altered hydrophobic tunnel accepting the 5'-substituents. Inhibitors of Rv0091 had weak cell growth effects on M. tuberculosis or Mycobacterium smegmatis but were lethal toward Helicobacter pylori, where the 5'-methylthioadenosine nucleosidase is essential in menaquinone biosynthesis. We propose that Rv0091 plays a role in 5'-deoxyadenosine recycling but is not essential for growth in these Mycobacteria.

Catalytic-site design for inverse heavy-enzyme isotope effects in human purine nucleoside phosphorylase

Heavy-enzyme isotope effects (¹⁵N-, ¹³C-, and ²H-labeled protein) explore mass-dependent vibrational modes linked to catalysis. Transition path-sampling (TPS) calculations have predicted femtosecond dynamic coupling at the catalytic site of human purine nucleoside phosphorylase (PNP). Coupling is observed in heavy PNPs, where slowed barrier crossing caused a normal heavy-enzyme isotope effect (k_chemlight/k_chemheavy > 1.0). We used TPS to design mutant F159Y PNP, predicted to improve barrier crossing for heavy F159Y PNP, an attempt to generate a rare inverse heavy-enzyme isotope effect (k_chemlight/k_chemheavy < 1.0). Steady-state kinetic comparison of light and heavy native PNPs to light and heavy F159Y PNPs revealed similar kinetic properties. Pre-steady-state chemistry was slowed 32-fold in F159Y PNP. Pre-steady-state chemistry compared heavy and light native and F159Y PNPs and found a normal heavy-enzyme isotope effect of 1.31 for native PNP and an inverse effect of 0.75 for F159Y PNP. Increased isotopic mass in F159Y PNP causes more efficient transition state formation. Independent validation of the inverse isotope effect for heavy F159Y PNP came from commitment to catalysis experiments. Most heavy enzymes demonstrate normal heavy-enzyme isotope effects, and F159Y PNP is a rare example of an inverse effect. Crystal structures and TPS dynamics of native and F159Y PNPs explore the catalytic-site geometry associated with these catalytic changes. Experimental validation of TPS predictions for barrier crossing establishes the connection of rapid protein dynamics and vibrational coupling to enzymatic transition state passage.

Kinetic Isotope Effects and Transition State Structure for Human Phenylethanolamine N-Methyltransferase

Phenylethanolamine N-methyltransferase (PNMT) catalyzes the S-adenosyl-l-methionine (SAM)-dependent conversion of norepinephrine to epinephrine. Epinephrine has been associated with critical processes in humans including the control of respiration and blood pressure. Additionally, PNMT activity has been suggested to play a role in hypertension and Alzheimer's disease. In the current study, labeled SAM substrates were used to measure primary methyl-14C and 36S and secondary methyl-3H, 5'-3H, and 5'-14C intrinsic kinetic isotope effects for human PNMT. The transition state of human PNMT was modeled by matching kinetic isotope effects predicted via quantum chemical calculations to intrinsic values. The model provides information on the geometry and electrostatics of the human PNMT transition state structure and indicates that human PNMT catalyzes the formation of epinephrine through an early SN2 transition state in which methyl transfer is rate-limiting.

Effect of short term external perturbations on bacterial ecology and activities in a partial nitritation and anammox reactor

This research investigated the short term effects of temperature changes (lasting 2-4weeks each) from 35±2°C to 21±2°C and 13±2°C and sulfide toxicity on partial nitrification-anammox (PN/A) system. Temperatures below 20°C and sulfide content as low as 5mgSL(-1) affected both aerobic and anaerobic catabolic activities of ammonia oxidation and the expression of related functional gene markers. The activity of AOB was inversely correlated with ammonium monooxygenase (amoA) gene expression. In contrast, the activity of AMX bacteria was positively correlated with the expression of their hydrazine synthase (hzsA) gene. Although the overall activities of AMX bacteria decreased at lower temperatures, the AMX bacteria were still active at the low temperatures. The inverse correlation between amoA gene expressions and the corresponding AOB activities was surprising. 16S rDNA based high throughput amplicon sequencing revealed the dominance of Chloroflexi, Planctomycetes and Proteobacteria phyla the distribution of which changed with temperature changes.

Transition State Structure and Inhibition of Rv0091, a 5'-Deoxyadenosine/5'-methylthioadenosine Nucleosidase from Mycobacterium tuberculosis

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) is a bacterial enzyme that catalyzes the hydrolysis of the N-ribosidic bond in 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (SAH). MTAN activity has been linked to quorum sensing pathways, polyamine biosynthesis, and adenine salvage. Previously, the coding sequence of Rv0091 was annotated as a putative MTAN in Mycobacterium tuberculosis. Rv0091 was expressed in Escherichia coli, purified to homogeneity, and shown to be a homodimer, consistent with MTANs from other microorganisms. Substrate specificity for Rv0091 gave a preference for 5'-deoxyadenosine relative to MTA or SAH. Intrinsic kinetic isotope effects (KIEs) for the hydrolysis of [1'-(3)H], [1'-(14)C], [5'-(3)H2], [9-(15)N], and [7-(15)N]MTA were determined to be 1.207, 1.038, 0.998, 1.021, and 0.998, respectively. A model for the transition state structure of Rv0091 was determined by matching KIE values predicted via quantum chemical calculations to the intrinsic KIEs. The transition state shows a substantial loss of C1'-N9 bond order, well-developed oxocarbenium character of the ribosyl ring, and weak participation of the water nucleophile. Electrostatic potential surface maps for the Rv0091 transition state structure show similarity to DADMe-immucillin transition state analogues. DADMe-immucillin transition state analogues showed strong inhibition of Rv0091, with the most potent inhibitor (5'-hexylthio-DADMe-immucillinA) displaying a Ki value of 87 pM.

Transition state for the NSD2-catalyzed methylation of histone H3 lysine 36

Nuclear receptor SET domain containing protein 2 (NSD2) catalyzes the methylation of histone H3 lysine 36 (H3K36). It is a determinant in Wolf-Hirschhorn syndrome and is overexpressed in human multiple myeloma. Despite the relevance of NSD2 to cancer, there are no potent, selective inhibitors of this enzyme reported. Here, a combination of kinetic isotope effect measurements and quantum chemical modeling was used to provide subangstrom details of the transition state structure for NSD2 enzymatic activity. Kinetic isotope effects were measured for the methylation of isolated HeLa cell nucleosomes by NSD2. NSD2 preferentially catalyzes the dimethylation of H3K36 along with a reduced preference for H3K36 monomethylation. Primary Me-(14)C and (36)S and secondary Me-(3)H3, Me-(2)H3, 5'-(14)C, and 5'-(3)H2 kinetic isotope effects were measured for the methylation of H3K36 using specifically labeled S-adenosyl-l-methionine. The intrinsic kinetic isotope effects were used as boundary constraints for quantum mechanical calculations for the NSD2 transition state. The experimental and calculated kinetic isotope effects are consistent with an SN2 chemical mechanism with methyl transfer as the first irreversible chemical step in the reaction mechanism. The transition state is a late, asymmetric nucleophilic displacement with bond separation from the leaving group at (2.53 Å) and bond making to the attacking nucleophile (2.10 Å) advanced at the transition state. The transition state structure can be represented in a molecular electrostatic potential map to guide the design of inhibitors that mimic the transition state geometry and charge.

Transition States, analogues, and drug development

Enzymes achieve their transition states by dynamic conformational searches on the femtosecond to picosecond time scale. Mimics of reactants at enzymatic transition states bind tightly to enzymes by stabilizing the conformation optimized through evolution for transition state formation. Instead of forming the transient transition state geometry, transition state analogues convert the short-lived transition state to a stable thermodynamic state. Enzymatic transition states are understood by combining kinetic isotope effects and computational chemistry. Analogues of the transition state can bind millions of times more tightly than substrates and show promise for drug development for several targets.

Double-headed sulfur-linked amino acids as first inhibitors for betaine-homocysteine S-methyltransferase 2

Betaine-homocysteine S-methyltransferase 2 (BHMT-2) catalyzes the transfer of a methyl group from S-methylmethionine to l-homocysteine, yielding two molecules of l-methionine. It is one of three homocysteine methyltransferases in mammals, but its overall contribution to homocysteine remethylation and sulfur amino acid homeostasis is not known. Moreover, recombinant BHMT-2 is highly unstable, which has slowed research on its structural and catalytic properties. In this study, we have prepared the first series of BHMT-2 inhibitors to be described, and we have tested them with human recombinant BHMT-2 that has been stabilized by copurification with human recombinant BHMT. Among the compounds synthesized, (2S,8RS,11RS)-5-thia-2,11-diamino-8-methyldodecanedioic acid (11) was the most potent (K(i)(app) ∼77 nM) and selective inhibitor of BHMT-2. Compound 11 only weakly inhibited human BHMT (IC(50) about 77 μM). This compound (11) may be useful in future in vivo studies to probe the physiological significance of BHMT-2 in sulfur amino acid metabolism.

Transition states of native and drug-resistant HIV-1 protease are the same

HIV-1 protease is an important target for the treatment of HIV/AIDS. However, drug resistance is a persistent problem and new inhibitors are needed. An approach toward understanding enzyme chemistry, the basis of drug resistance, and the design of powerful inhibitors is to establish the structure of enzymatic transition states. Enzymatic transition structures can be established by matching experimental kinetic isotope effects (KIEs) with theoretical predictions. However, the HIV-1 protease transition state has not been previously resolved using these methods. We have measured primary (14)C and (15)N KIEs and secondary (3)H and (18)O KIEs for native and multidrug-resistant HIV-1 protease (I84V). We observed (14)C KIEs ((14)V/K) of 1.029 ± 0.003 and 1.025 ± 0.005, (15)N KIEs ((15)V/K) of 0.987 ± 0.004 and 0.989 ± 0.003, (18)O KIEs ((18)V/K) of 0.999 ± 0.003 and 0.993 ± 0.003, and (3)H KIEs ((3)V/K) KIEs of 0.968 ± 0.001 and 0.976 ± 0.001 for the native and I84V enzyme, respectively. The chemical reaction involves nucleophilic water attack at the carbonyl carbon, proton transfer to the amide nitrogen leaving group, and C-N bond cleavage. A transition structure consistent with the KIE values involves proton transfer from the active site Asp-125 (1.32 Å) with partial hydrogen bond formation to the accepting nitrogen (1.20 Å) and partial bond loss from the carbonyl carbon to the amide leaving group (1.52 Å). The KIEs measured for the native and I84V enzyme indicate nearly identical transition states, implying that a true transition-state analogue should be effective against both enzymes.

Structural basis of the substrate specificity of Bacillus cereus adenosine phosphorylase

Purine nucleoside phosphorylases catalyze the phosphorolytic cleavage of the glycosidic bond of purine (2'-deoxy)nucleosides, generating the corresponding free base and (2'-deoxy)-ribose 1-phosphate. Two classes of PNPs have been identified: homotrimers specific for 6-oxopurines and homohexamers that accept both 6-oxopurines and 6-aminopurines. Bacillus cereus adenosine phosphorylase (AdoP) is a hexameric PNP; however, it is highly specific for 6-aminopurines. To investigate the structural basis for the unique substrate specificity of AdoP, the active-site mutant D204N was prepared and kinetically characterized and the structures of the wild-type protein and the D204N mutant complexed with adenosine and sulfate or with inosine and sulfate were determined at high resolution (1.2-1.4 Å). AdoP interacts directly with the preferred substrate through a hydrogen-bond donation from the catalytically important residue Asp204 to N7 of the purine base. Comparison with Escherichia coli PNP revealed a more optimal orientation of Asp204 towards N7 of adenosine and a more closed active site. When inosine is bound, two water molecules are interposed between Asp204 and the N7 and O6 atoms of the nucleoside, thus allowing the enzyme to find alternative but less efficient ways to stabilize the transition state. The mutation of Asp204 to asparagine led to a significant decrease in catalytic efficiency for adenosine without affecting the efficiency of inosine cleavage.

Kinetic Isotope Effects from QM/MM Subset Hessians: "Cut-Off" Analysis for SN2 Methyl Transfer in Solution

Isotopic partition-function ratios and kinetic isotope effects for reaction of S-adenosylmethionine with catecholate in water are evaluated using a subset of 324 atoms within its surrounding aqueous environment at the AM1/TIP3P level. Two alternative methods for treating motion in the six librational degrees of freedom of the subset atoms relative to their environment are compared. A series of successively smaller subset Hessians are generated by cumulative deletion of rows and columns from the initial 972 × 972 Hessian. We find that it is better to treat these librations as vibrations than as translations and rotations and that there is no need to invoke the Teller-Redlich product rule. The validity of "cut-off" procedures for computation of isotope effects with truncated atomic subsets is assessed: to ensure errors in ln(KIE) < 1% (or 2% for the quantum-corrected KIE) for all isotopic substitutions considered, it is necessary to use a less-restrictive procedure than is suggested by the familiar two-bond cutoff rule.

Free energy study of the catalytic mechanism of Trypanosoma cruzi trans-sialidase. From the Michaelis complex to the covalent intermediate

Trypanosoma cruzi trans-sialidase (TcTS) is a crucial enzyme for the infection of Trypanosoma cruzi, the protozoa responsible for Chagas' disease in humans. It catalyzes the transfer of sialic acids from the host's glycoconjugates to the parasite's glycoconjugates. Based on kinetic isotope effect (KIE) studies, a strong nucleophilic participation at the transition state could be determined, and recently, elaborate experiments used 2-deoxy-2,3-difluorosialic acid as substrate and were able to trap a long-lived covalent intermediate (CI) during the catalytic mechanism. In this paper, we compute the KIE and address the entire mechanistic pathway of the CI formation step in TcTS using computational tools. Particularly, the free energy results indicate that in the transition state there is a strong nucleophilic participation of Tyr342, and after this, the system collapsed into a stable CI. We find that there is no carbocation intermediate for this reaction. By means of the energy decomposition method, we identify the residues that have the biggest influence on catalysis. This study facilitates the understanding of the catalytic mechanism of TcTS and can serve as a guide for future inhibitor design studies.

Enzymatic transition states, transition-state analogs, dynamics, thermodynamics, and lifetimes

Experimental analysis of enzymatic transition-state structures uses kinetic isotope effects (KIEs) to report on bonding and geometry differences between reactants and the transition state. Computational correlation of experimental values with chemical models permits three-dimensional geometric and electrostatic assignment of transition states formed at enzymatic catalytic sites. The combination of experimental and computational access to transition-state information permits (a) the design of transition-state analogs as powerful enzymatic inhibitors, (b) exploration of protein features linked to transition-state structure, (c) analysis of ensemble atomic motions involved in achieving the transition state, (d) transition-state lifetimes, and (e) separation of ground-state (Michaelis complexes) from transition-state effects. Transition-state analogs with picomolar dissociation constants have been achieved for several enzymatic targets. Transition states of closely related isozymes indicate that the protein's dynamic architecture is linked to transition-state structure. Fast dynamic motions in catalytic sites are linked to transition-state generation. Enzymatic transition states have lifetimes of femtoseconds, the lifetime of bond vibrations. Binding isotope effects (BIEs) reveal relative reactant and transition-state analog binding distortion for comparison with actual transition states.

Uridine phosphorylase from Trypanosoma cruzi: kinetic and chemical mechanisms

The reversible phosphorolysis of uridine to generate uracil and ribose 1-phosphate is catalyzed by uridine phosphorylase and is involved in the pyrimidine salvage pathway. We define the reaction mechanism of uridine phosphorylase from Trypanosoma cruzi by steady-state and pre-steady-state kinetics, pH-rate profiles, kinetic isotope effects from uridine, and solvent deuterium isotope effects. Initial rate and product inhibition patterns suggest a steady-state random kinetic mechanism. Pre-steady-state kinetics indicated no rate-limiting step after formation of the enzyme-products ternary complex, as no burst in product formation is observed. The limiting single-turnover rate constant equals the steady-state turnover number; thus, chemistry is partially or fully rate limiting. Kinetic isotope effects with [1'-(3)H]-, [1'-(14)C]-, and [5'-(14)C,1,3-(15)N(2)]uridine gave experimental values of (α-T)(V/K)(uridine) = 1.063, (14)(V/K)(uridine) = 1.069, and (15,β-15)(V/K)(uridine) = 1.018, in agreement with an A(N)D(N) (S(N)2) mechanism where chemistry contributes significantly to the overall rate-limiting step of the reaction. Density functional theory modeling of the reaction in gas phase supports an A(N)D(N) mechanism. Solvent deuterium kinetic isotope effects were unity, indicating that no kinetically significant proton transfer step is involved at the transition state. In this N-ribosyl transferase, proton transfer to neutralize the leaving group is not part of transition state formation, consistent with an enzyme-stabilized anionic uracil as the leaving group. Kinetic analysis as a function of pH indicates one protonated group essential for catalysis and for substrate binding.

Transition-state analysis of Trypanosoma cruzi uridine phosphorylase-catalyzed arsenolysis of uridine

Uridine phosphorylase catalyzes the reversible phosphorolysis of uridine and 2'-deoxyuridine to generate uracil and (2-deoxy)ribose 1-phosphate, an important step in the pyrimidine salvage pathway. The coding sequence annotated as a putative nucleoside phosphorylase in the Trypanosoma cruzi genome was overexpressed in Escherichia coli , purified to homogeneity, and shown to be a homodimeric uridine phosphorylase, with similar specificity for uridine and 2'-deoxyuridine and undetectable activity toward thymidine and purine nucleosides. Competitive kinetic isotope effects (KIEs) were measured and corrected for a forward commitment factor using arsenate as the nucleophile. The intrinsic KIEs are: 1'-(14)C = 1.103, 1,3-(15)N(2) = 1.034, 3-(15)N = 1.004, 1-(15)N = 1.030, 1'-(3)H = 1.132, 2'-(2)H = 1.086, and 5'-(3)H(2) = 1.041 for this reaction. Density functional theory was employed to quantitatively interpret the KIEs in terms of transition-state structure and geometry. Matching of experimental KIEs to proposed transition-state structures suggests an almost synchronous, S(N)2-like transition-state model, in which the ribosyl moiety possesses significant bond order to both nucleophile and leaving groups. Natural bond orbital analysis allowed a comparison of the charge distribution pattern between the ground-state and the transition-state models.

Arsenate and phosphate as nucleophiles at the transition states of human purine nucleoside phosphorylase

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of 6-oxypurine (2'-deoxy)ribonucleosides, generating (2-deoxy)ribose 1-phosphate and the purine base. Transition-state models for inosine cleavage have been proposed with bovine, human, and malarial PNPs using arsenate as the nucleophile, since kinetic isotope effects (KIEs) are obscured on phosphorolysis due to high commitment factors. The Phe200Gly mutant of human PNP has low forward and reverse commitment factors in the phosphorolytic reaction, permitting the measurement of competitive intrinsic KIEs on both arsenolysis and phosphorolysis of inosine. The intrinsic 1'-(14)C, 1'-(3)H, 2'-(2)H, 9-(15)N, and 5'-(3)H(2) KIEs for inosine were measured for arsenolysis and phosphorolysis. Except for the remote 5'-(3)H(2), and some slight difference between the 2'-(2)H KIEs, all isotope effects originating in the reaction coordinate are the same within experimental error. Hence, arsenolysis and phosphorolysis proceed through closely related transition states. Although electrostatically similar, the volume of arsenate is greater than phosphate and supports a steric influence to explain the differences in the 5'-(3)H(2) KIEs. Density functional theory calculations provide quantitative models of the transition states for Phe200Gly human PNP-catalyzed arsenolysis and phosphorolysis, selected upon matching calculated and experimental KIEs. The models confirm the striking resemblance between the transition states for the two reactions.

Pyrophosphate interactions at the transition states of Plasmodium falciparum and human orotate phosphoribosyltransferases

Orotate phosphoribosyltransferases from Plasmodium falciparum and human sources (PfOPRT and HsOPRT) use orotidine as a slow substrate in the pyrophosphorolysis reaction. With orotidine, intrinsic kinetic isotope effects (KIEs) can be measured for pyrophosphorolysis, providing the first use of pyrophosphate (PPi) in solving an enzymatic transition state. Transition-state structures of PfOPRT and HsOPRT were solved through quantum chemical matching of computed and experimental intrinsic KIEs and can be compared to transition states solved with pyrophosphate analogues as slow substrates. PfOPRT and HsOPRT are characterized by late transition states with fully dissociated orotate, well-developed ribocations, and weakly bonded PPi nucleophiles. The leaving orotates are 2.8 A distant from the anomeric carbons at the transition states. Weak participation of the PPi nucleophiles gives C1'-O(PPi) bond distances of approximately 2.3 A. These transition states are characterized by C2'-endo ribosyl pucker, based on the beta-secondary [2'-(3)H] KIEs. The geometry at the 5'-region is similar for both enzymes, with C3'-C4'-C5'-O5' dihedral angles near -170 degrees . These novel phosphoribosyltransferase transition states are similar to but occur earlier in the reaction coordinate than those previously determined with orotidine 5'-monophosphate and phosphonoacetic acid as substrates. The similarity between the transition states with different substrate analogues supports similar transition state structures imposed by PfOPRT and HsOPRT even with distinct reactants. We propose that the transition state similarity with different nucleophiles is determined, in part, by the geometric constraints imposed by the catalytic sites.

Transition states of Plasmodium falciparum and human orotate phosphoribosyltransferases

Orotate phosphoribosyltransferases (OPRT) catalyze the formation of orotidine 5'-monophosphate (OMP) from alpha-D-phosphoribosylpyrophosphate (PRPP) and orotate, an essential step in the de novo biosynthesis of pyrimidines. Pyrimidine de novo biosynthesis is required in Plasmodium falciparum , and thus OPRT of the parasite (PfOPRT) is a target for antimalarial drugs. De novo biosynthesis of pyrimidines is also a feature of rapidly proliferating cancer cells. Human OPRT (HsOPRT) is therefore a target for neoplastic and autoimmune diseases. One approach to the inhibition of OPRTs is through analogues that mimic the transition states of PfOPRT and HsOPRT. The transition state structures of these OPRTs were analyzed by kinetic isotope effects (KIEs), substrate specificity, and computational chemistry. With phosphonoacetic acid (PA), an analogue of pyrophosphate, the intrinsic KIEs of [1'-(14)C], [1, 3-(15)N(2)], [3-(15)N], [1'-(3)H], [2'-(3)H], [4'-(3)H], and [5'-(3)H(2)] are 1.034, 1.028, 0.997, 1.261, 1.116, 0.974, and 1.013 for PfOPRT and 1.035, 1.025, 0.993, 1.199, 1.129, 0.962, and 1.019 for HsOPRT, respectively. Transition state structures of PfOPRT and HsOPRT were determined computationally by matching the calculated and intrinsic KIEs. The enzymes form late associative D(N)*A(N)(double dagger) transition states with complete orotate loss and partially associative nucleophile. The C1'-O(PA) distances are approximately 2.1 A at these transition states. The modest [1'-(14)C] KIEs and large [1'-(3)H] KIEs are characteristic of D(N)*A(N)(double dagger) transition states. The large [2'-(3)H] KIEs indicate a ribosyl 2'-C-endo conformation at the transition states. p-Nitrophenyl beta-D-ribose 5'-phosphate is a poor substrate of PfOPRT and HsOPRT but is a nanomolar inhibitor, supporting a reaction coordinate with strong leaving group activation.

Synthesis of isotopically labeled P-site substrates for the ribosomal peptidyl transferase reaction

Isotopomers of the ribosomal P-site substrate, the trinucleotide peptide conjugate CCA-pcb (Zhong, M.; Strobel, S. A. Org. Lett. 2006, 8, 55-58), have been designed and synthesized in 26-35 steps. These include individual isotopic substitution at the alpha-hydrogen, carbonyl carbon, and carbonyl oxygen of the amino acid, the O2' and O3' of the adenosine, and a remote label in the N3 and N4 of both cytidines. These isotopomers were synthesized by coupling cytidylyl-(3',5')-cytidine phosphoramidite isotopomers as the common synthetic intermediates, with isotopically substituted A-Phe-cap-biotin (A-pcb). The isotopic enrichment is higher than 99% for 1-13C (Phe), 2-2H (Phe), and 3,4-15N2 (cytidine), 93% for 2'/3'-18 O (adenosine), and 64% for 1-18 O (Phe). A new synthesis of highly enriched [1-18 O2]phenylalanine has been developed. The synthesis of [3'-18 O]adenosine was improved by Lewis acid aided regioselective ring opening of the epoxide and by an economical SN2-SN2 method with high isotopic enrichment (93%). Such substrates are valuable for studies of the ribosomal peptidyl transferase reaction by complete kinetic isotope effect analysis and of other biological processes catalyzed by nucleic acid related enzymes, including polymerases, reverse transcriptases, ligases, nucleases, and ribozymes.

Enzymes as a special class of therapeutic target: clinical drugs and modes of action

Enzymes catalyze multistep chemical reactions and achieve phenomenal rate accelerations by matching protein and substrate chemical groups in the transition state. Inhibitors that take advantage of these chemical interactions are among the most potent and effective drugs known. Recently, three new enzyme targets have been validated by FDA approval of new enzyme inhibitor drugs. These include mitogen-activated protein kinase, renin, and dipeptidyl peptidase IV. The drugs against these enzymes engage important enzyme functional groups, such as the active site serine in dipeptidyl peptidase IV. Clinical and pre-clinical discovery programs also demonstrate the same theme, as evidenced by pM and fM transition state inhibitors of purine nucleoside phosphorylase, methylthioadenosine phosphorylase, and 5-methylthioadenosine/S-adenosylhomocysteine nucleosidase, and covalent substrate trapping in leu-tRNA synthetase. The catalytic chemistry of enzymes is the key to designing potent inhibitors and makes them a special class of drug target.

Binding isotope effects: boon and bane

Kinetic isotope effects are increasingly applied to investigate enzyme reactions and have been used to understand transition state structure, reaction mechanisms, quantum mechanical hydride ion tunneling and to design transition state analogue inhibitors. Binding isotope effects are an inherent part of most isotope effect measurements but are usually assumed to be negligible. More detailed studies have established surprisingly large binding isotope effects with lactate dehydrogenase, hexokinase, thymidine phosphorylase, and purine nucleoside phosphorylase. Binding reactants into catalytic sites immobilizes conformationally flexible groups, polarizes bonds, and distorts bond angle geometry, all of which generate binding isotope effects. Binding isotope effects are easily measured and provide high-resolution and detailed information on the atomic changes resulting from ligand-macromolecular interactions. Although binding isotope effects complicate kinetic isotope effect analysis, they also provide a powerful tool for finding atomic distortion in molecular interactions.

Inaccuracies in selected ion monitoring determination of isotope ratios obviated by profile acquisition: nucleotide 18O/16O measurements

Precise and accurate measurements of isotopologue distributions (IDs) in biological molecules are needed for determination of isotope effects, quantitation by isotope dilution, and quantification of isotope tracers employed in both metabolic and biophysical studies. While single ion monitoring (SIM) yields significantly greater sensitivity and signal/noise than profile-mode acquisitions, we show that small changes in the SIM window width and/or center can alter experimentally determined isotope ratios by up to 5%, resulting in significant inaccuracies. This inaccuracy is attributed to mass granularity, the differential distribution of digital data points across the m/z ranges sampled by SIM. Acquiring data in the profile mode and fitting the data to an equation describing a series of equally spaced and identically shaped peaks eliminates the inaccuracies associated with mass granularity with minimal loss of precision. Additionally a method of using the complete ID profile data that inherently corrects for "spillover" and for the natural-abundance ID has been used to determine 18O/16O ratios for 5',3'-guanosine bis-[18O1]phosphate and TM[18O1]P with precisions of approximately 0.005. The analysis protocol is also applied to quadrupole time-of-flight tandem mass spectrometry using [2-(18)O] arabinouridine and 3'-UM[18O1]P which enhances signal/noise and minimizes concerns for background contamination.

Transition state analogue discrimination by related purine nucleoside phosphorylases

Transition state analogues of PNP, the Immucillins and DADMe-Immucillins, were designed to match transition state features of bovine and human PNPs, respectively. The inhibitors with or without the hydroxyl and hydroxymethyl groups of the substrate demonstrate that inhibitor geometry mimicking that of the transition state confers binding affinity discrimination. This finding is remarkable since crystallographic analysis indicates complete conservation of active site residues and contacts to ligands in human and bovine PNPs.

Transition-state variation in human, bovine, and Plasmodium falciparum adenosine deaminases

Adenosine deaminases (ADAs) from human, bovine, and Plasmodium falciparum sources were analyzed by kinetic isotope effects (KIEs) and shown to have distinct but related transition states. Human adenosine deaminase (HsADA) is present in most mammalian cells and is involved in B- and T-cell development. The ADA from Plasmodium falciparum (PfADA) is essential in this purine auxotroph, and its inhibition is expected to have therapeutic effects for malaria. Therefore, ADA is of continuing interest for inhibitor design. Stable structural mimics of ADA transition states are powerful inhibitors. Here we report the transition-state structures of PfADA, HsADA, and bovine ADA (BtADA) solved using competitive kinetic isotope effects (KIE) and density functional calculations. Adenines labeled at [6-13C], [6-15N], [6-13C, 6-15N], and [1-15N] were synthesized and enzymatically coupled with [1'-14C] ribose to give isotopically labeled adenosines as ADA substrates for KIE analysis. [6-13C], [6-15N], and [1-15N]adenosines reported intrinsic KIE values of (1.010, 1.011, 1.009), (1.005, 1.005, 1.002), and (1.004, 1.001, 0.995) for PfADA, HsADA, and BtADA, respectively. The differences in intrinsic KIEs reflect structural alterations in the transition states. The [1-15N] KIEs and computational modeling results indicate that PfADA, HsADA, and BtADA adopt early SNAr transition states, where N1 protonation is partial and the bond order to the attacking hydroxyl nucleophile is nearly complete. The key structural variation among PfADA, HsADA, and BtADA transition states lies in the degree of N1 protonation with the decreased bond lengths of 1.92, 1.55, and 1.28 A, respectively. Thus, PfADA has the earliest and BtADA has the most developed transition state. This conclusion is consistent with the 20-36-fold increase of kcat in comparing PfADA with HsADA and BtADA.

Influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine. A theoretical study

The influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine has been studied using density functional theory (DFT) methods. For the neutral system, two different pathways (with retention and inversion of configuration at the C1' anomeric carbon) have been found, both of them consisting of two steps and involving the formation of a dihydrofurane-like intermediate. The Gibbs free energy barrier for the first step is very high in both cases (53 and 46 kcal/mol for the process with inversion and with retention, respectively). However, the N7-protonated system shows a very different mechanism which consists of two steps. The first one leads to the formation of an oxacarbenium ion intermediate, with a Gibbs free energy barrier of 27 kcal/mol, and the second one corresponds to the nucleophilic attack of the water molecule to the oxacarbenium ion and takes place with a barrier of 1.3 kcal/mol. Thus, these results agree with a stepwise SN1 mechanism (DN*AN), with a discrete intermediate formed between the leaving group and the nucleophile approach, and show that N7 protonation strongly catalyzes the hydrolysis of the N-glycosidic bond, making the guanine a better leaving group. Finally, kinetic isotope effects have been calculated for the protonated system, and the results obtained are in very good agreement with experimental data for analogous systems.

Neighboring group participation in the transition state of human purine nucleoside phosphorylase

The X-ray crystal structures of human purine nucleoside phosphorylase (PNP) with bound inosine or transition-state analogues show His257 within hydrogen bonding distance of the 5'-hydroxyl. The mutants His257Phe, His257Gly, and His257Asp exhibited greatly decreased affinity for Immucillin-H (ImmH), binding this mimic of an early transition state as much as 370-fold (Km/Ki) less tightly than native PNP. In contrast, these mutants bound DADMe-ImmH, a mimic of a late transition state, nearly as well as the native enzyme. These results indicate that His257 serves an important role in the early stages of transition-state formation. Whereas mutation of His257 resulted in little variation in the PNP x DADMe-ImmH x SO4 structures, His257Phe x ImmH x PO4 showed distortion at the 5'-hydroxyl, indicating the importance of H-bonding in positioning this group during progression to the transition state. Binding isotope effect (BIE) and kinetic isotope effect (KIE) studies of the remote 5'-(3)H for the arsenolysis of inosine with native PNP revealed a BIE of 1.5% and an unexpectedly large intrinsic KIE of 4.6%. This result is interpreted as a moderate electronic distortion toward the transition state in the Michaelis complex with continued development of a similar distortion at the transition state. The mutants His257Phe, His257Gly, and His257Asp altered the 5'-(3)H intrinsic KIE to -3, -14, and 7%, respectively, while the BIEs contributed 2, 2, and -2%, respectively. These surprising results establish that forces in the Michaelis complex, reported by the BIEs, can be reversed or enhanced at the transition state.

Crystal structure of tryptophanyl-tRNA synthetase complexed with adenosine-5' tetraphosphate: evidence for distributed use of catalytic binding energy in amino acid activation by class I aminoacyl-tRNA synthetases

Tryptophanyl-tRNA synthetase (TrpRS) is a functionally dimeric ligase, which specifically couples hydrolysis of ATP to AMP and pyrophosphate to the formation of an ester bond between tryptophan and the cognate tRNA. TrpRS from Bacillus stearothermophilus binds the ATP analogue, adenosine-5' tetraphosphate (AQP) competitively with ATP during pyrophosphate exchange. Estimates of binding affinity from this competitive inhibition and from isothermal titration calorimetry show that AQP binds 200 times more tightly than ATP both under conditions of induced-fit, where binding is coupled to an unfavorable conformational change, and under exchange conditions, where there is no conformational change. These binding data provide an indirect experimental measurement of +3.0 kcal/mol for the conformational free energy change associated with induced-fit assembly of the active site. Thermodynamic parameters derived from the calorimetry reveal very modest enthalpic changes, consistent with binding driven largely by a favorable entropy change. The 2.5 A structure of the TrpRS:AQP complex, determined de novo by X-ray crystallography, resembles that of the previously described, pre-transition state TrpRS:ATP complexes. The anticodon-binding domain untwists relative to the Rossmann-fold domain by 20% of the way toward the orientation observed for the Products complex. An unexpected tetraphosphate conformation allows the gamma and deltad phosphate groups to occupy positions equivalent to those occupied by the beta and gamma phosphates of ATP. The beta-phosphate effects a 1.11 A extension that relocates the alpha-phosphate toward the tryptophan carboxylate while the PPi mimic moves deeper into the KMSKS loop. This configuration improves interactions between enzyme and nucleotide significantly and uniformly in the adenosine and PPi binding subsites. A new hydrogen bond forms between S194 from the class I KMSKS signature sequence and the PPi mimic. These complementary thermodynamic and structural data are all consistent with the conclusion that the tetraphosphate mimics a transition-state in which the KMSKS loop develops increasingly tight bonds to the PPi leaving group, weakening linkage to the Palpha as it is relocated by an energetically favorable domain movement. Consistent with extensive mutational data on Tyrosyl-tRNA synthetase, this aspect of the mechanism develops high transition-state affinity for the adenosine and pyrophosphate moieties, which move significantly, relative to one another, during the catalytic step.

Dependence of Transition State Structure on Substrate: The Intrinsic C-13 Kinetic Isotope Effect Is Different for Physiological and Slow Substrates of the Ornithine Decarboxylase Reaction Because of Different Hydrogen Bonding Structures

Ornithine decarboxylase is the first and the rate-controlling enzyme in polyamine biosynthesis; it decarboxylates l-ornithine to form the diamine putrescine. We present calculations performed using a combined quantum mechanical and molecular mechanical (QM/MM) method with the AM1 semiempirical Hamiltonian for the wild-type ornithine decarboxylase reaction with ornithine (the physiological substrate) and lysine (a "slow" substrate) and for mutant E274A with ornithine substrate. The dynamical method is variational transition state theory with quantized vibrations. We employ a single reaction coordinate equal to the carbon-carbon distance of the dissociating bond, and we find a large difference between the intrinsic kinetic isotope effect for the physiological substrate, which equals 1.04, and that for the slow substrate, which equals 1.06. This shows that, contrary to a commonly accepted assumption, kinetic isotope effects on slow substrates are not always good models of intrinsic kinetic isotope effects on physiological substrates. Furthermore, analysis of free-energy-based samples of transition state structures shows that the differences in kinetic isotope effects may be traced to different numbers of hydrogen bonds at the different transition states of the different reactions.

Femtomolar transition state analogue inhibitors of 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase from Escherichia coli

Escherichia coli 5'-methylthioadenosine/S-adenosyl-homocysteine nucleosidase (MTAN) hydrolyzes its substrates to form adenine and 5-methylthioribose (MTR) or S-ribosylhomocysteine (SRH). 5'-Methylthioadenosine (MTA) is a by-product of polyamine synthesis and SRH is a precursor to the biosynthesis of one or more quorum sensing autoinducer molecules. MTAN is therefore involved in quorum sensing, recycling MTA from the polyamine pathway via adenine phosphoribosyltransferase and recycling MTR to methionine. Hydrolysis of MTA by E. coli MTAN involves a highly dissociative transition state with ribooxacarbenium ion character. Iminoribitol mimics of MTA at the transition state of MTAN were synthesized and tested as inhibitors. 5'-Methylthio-Immucillin-A (MT-ImmA) is a slow-onset tight-binding inhibitor giving a dissociation constant (K(i)(*)) of 77 pm. Substitution of the methylthio group with a p-Cl-phenylthio group gives a more powerful inhibitor with a dissociation constant of 2 pm. DADMe-Immucillins are better inhibitors of E. coli MTAN, since they are more closely related to the highly dissociative nature of the transition state. MT-DADMe-Immucillin-A binds with a K(i)(*) value of 2 pm. Replacing the 5'-methyl group with other hydrophobic groups gave 17 transition state analogue inhibitors with dissociation constants from 10(-12) to 10(-14) m. The most powerful inhibitor was 5'-p-Cl-phenylthio-DADMe-Immucillin-A (pClPhT-DADMe-ImmA) with a K(i)(*) value of 47 fm (47 x 10(-15) m). These are among the most powerful non-covalent inhibitors reported for any enzyme, binding 9-91 million times tighter than the MTA and SAH substrates, respectively. The inhibitory potential of these transition state analogue inhibitors supports a transition state structure closely resembling a fully dissociated ribooxacarbenium ion. Powerful inhibitors of MTAN are candidates to disrupt key bacterial pathways including methylation, polyamine synthesis, methionine salvage, and quorum sensing. The accompanying article reports crystal structures of MTAN with these analogues.

Soybean epoxide hydrolase: identification of the catalytic residues and probing of the reaction mechanism with secondary kinetic isotope effects

Soybean epoxide hydrolase catalyzes the oxirane ring opening of 9,10-epoxystearate via a two-step mechanism involving the formation of an alkylenzyme intermediate, which, in contrast to most epoxide hydrolases studied so far, was found to be the rate-limiting step. We have probed residues potentially involved in catalysis by site-directed mutagenesis. Mutation of His(320), a residue predicted from sequence analysis to belong to the catalytic triad of the enzyme, considerably slowed down the second half-reaction. This kinetic manipulation provoked an accumulation of the reaction intermediate, which could be trapped and characterized by electrospray ionization mass spectrometry. As expected, mutation of Asp(126) totally abolished the activity of the enzyme from its crucial function as nucleophile involved in the formation of the alkylenzyme. In line with its role as the partner of His(320) in the "charge relay system," mutation of Asp(285) dramatically reduced the rate of catalysis. However, the mutant D285L still exhibited a very low residual activity, which, by structural analysis and mutagenesis, has been tentatively attributed to Glu(195), another acidic residue of the active site. Our studies have also confirmed the fundamental role of the conserved Tyr(175) and Tyr(255) residues, which are believed to activate the oxirane ring. Finally, we have determined the secondary tritium kinetic isotope effects on the epoxide opening step of 9,10-epoxystearate. The large observed values, i.e. (T)(V/K(m)) approximately 1.30, can be interpreted by the occurrence of a very late transition state in which the epoxide bond is broken before the nucleophilic attack by Asp(126) takes place.

New catalytic mechanism for human purine nucleoside phosphorylase

Human purine nucleoside phosphorylase has been submitted to intensive structure-based design of inhibitors, most of them using low-resolution structures of human PNP. Recently, several structures of human PNP have been reported, which allowed redefinition of the active site and understanding of the structural basis for inhibition of PNP by acyclovir and immucillin-H. Based on previously solved human PNP structures, we proposed here a new catalytic mechanism for human PNP, which is supported by crystallographic studies and explains previously determined kinetic data.

Enzymatic transition states: thermodynamics, dynamics and analogue design

Kinetic isotope effects and computational chemistry have defined the transition state structures for several members of the N-ribosyltransferase family. Transition state analogues designed to mimic their cognate transition state structures are among the most powerful enzyme inhibitors. In complexes of N-ribosyltransferases with their transition state analogues, the dynamic nature of the transition state is converted to an ordered, thermodynamic structure closely related to the transition state. This phenomenon is documented by peptide bond H/D exchange, crystallography and computational chemistry. Complexes with substrate, transition state and product analogues reveal reaction coordinate motion and catalytic interactions. Isotope-edited spectroscopic analysis and binding specificity of these complexes provides information about specific enzyme-transition state contacts. In combination with protein dynamic QM/MM models, it is proposed that the transition state is reached by stochastic dynamic excursions of the protein groups near the substrates in the closed conformation. Examples from fully dissociated (D(N) *A(N)), hybrid (D(N)A(N)) and symmetric nucleophilic displacement (A(N)D(N)) transition states are found in the N-ribosyltransferases. The success of transition state analogue inhibitor design based on kinetic isotope effects validates this approach to understanding enzymatic transition states.

Crystallographic structure of PNP from Mycobacterium tuberculosis at 1.9Å resolution

Even being a bacterial purine nucleoside phosphorylase (PNP), which normally shows hexameric folding, the Mycobacterium tuberculosis PNP (MtPNP) resembles the mammalian trimeric structure. The crystal structure of the MtPNP apoenzyme was solved at 1.9 A resolution. The present work describes the first structure of MtPNP in complex with phosphate. In order to develop new insights into the rational drug design, conformational changes were profoundly analyzed and discussed. Comparisons over the binding sites were specially studied to improve the discussion about the selectivity of potential new drugs.

An oxocarbenium-ion intermediate of a ribozyme reaction indicated by kinetic isotope effects

Many of the enzymes that catalyze reactions at nucleotide glycosidic linkages proceed through either a reactive oxocarbenium-ion intermediate or a transition state with considerable oxocarbenium character. To investigate how an RNA active site deals with the catalytic challenge of nucleotide synthesis, we probed the transition state of a ribozyme able to promote the formation of a pyrimidine nucleotide. Primary and secondary kinetic isotope effects indicate that this ribozyme stabilizes a highly dissociative reaction with considerable sp2 hybridization and negligible bond order between the departing pyrophosphate leaving group and the anomeric carbon. The small primary 13C isotope effect of 1.002 +/- 0.003 indicates that the reaction is likely to be less concerted than that observed for protein nucleotide synthesis enzymes, which typically have primary 13C isotope effects of 1.02-1.03. The dissociative nature of the ribozyme reaction most resembles the reaction of some hydrolytic enzymes, such as uracil DNA glycosylase, which uses the negative charges found in the phosphodiester backbone of its DNA substrate to transiently stabilize an oxocarbenium ion during hydrolysis. The detectable hydrolysis observed in the ribozyme reaction indicates that shielding of this reactive intermediate from water is a significant challenge for RNA, which protein enzymes that synthesize nucleotides have managed to overcome during evolution, apparently by the utilization of more concerted chemistry.

Adenine release is fast in MutY-catalyzed hydrolysis of G:A and 8-Oxo-G:A DNA mismatches

MutY, a DNA repair enzyme, is unusual in that it binds exceedingly tightly to its products after the chemical steps of catalysis. Until now it was not known whether the product being released in the rate-limiting step was DNA, adenine, or both. MutY hydrolyzes adenine from 8-oxo-G:A (OG:A) base pair mismatches as the first step in the base excision repair pathway, as well as from G:A mismatches. The products are adenine and DNA containing an apurinic (AP) site. Tight product binding may have a physiological role in preventing further damage at the OG:AP site. We developed a rate assay using [8-14C]adenine in OG:A or G:A mismatches that distinguishes between adenine hydrolysis and adenine release. [8-14C]Adenine was released quickly from the MutY.AP-DNA.[8-14C]adenine complex, with a rate constant greater than 5 min-1. This was much faster than the rate-limiting step, at 0.006-0.015 min-1. Gel retardation experiments showed that AP-DNA release was very slow, consistent with it being the rate-limiting step. Thus, the kinetic mechanism involves fast adenine release after hydrolysis followed by rate-limiting AP-DNA release. Adenine appears to be buried deep in the protein.DNA interface, but there is enough flexibility or open space for it to dissociate from the MutY.APDNA.adenine complex. These results have implications for the catalytic mechanism of MutY.

Biosynthesis of the 7-deazaguanosine hypermodified nucleosides of transfer RNA

Transfer RNA (tRNA) is structurally unique among nucleic acids in harboring an astonishing diversity of post-transcriptionally modified nucleoside. Two of the most radically modified nucleosides known to occur in tRNA are queuosine and archaeosine, both of which are characterized by a 7-deazaguanosine core structure. In spite of the phylogenetic segregation observed for these nucleosides (queuosine is present in Eukarya and Bacteria, while archaeosine is present only in Archaea), their structural similarity suggested a common biosynthetic origin, and recent biochemical and genetic studies have provided compelling evidence that a significant portion of their biosynthesis may in fact be identical. This review covers current understanding of the physiology and biosynthesis of these remarkable nucleosides, with particular emphasis on the only two enzymes that have been discovered in the pathways: tRNA-guanine transglycosylase (TGT), which catalyzes the insertion of a modified base into the polynucleotide with the concomitant elimination of the genetically encoded guanine in the biosynthesis of both nucleosides, and S-adenosylmethionine:tRNA ribosyltransferase-isomerase (QueA), which catalyzes the penultimate step in the biosynthesis of queuosine, the construction of the carbocyclic side chain.