Nucleophilic Participation in the Transition State for Human Thymidine Phosphorylase

Engineering of the Recombinant Expression and PEGylation Efficiency of the Therapeutic Enzyme Human Thymidine Phosphorylase

Human thymidine phosphorylase (HsTP) is an enzyme with important implications in the field of rare metabolic diseases. Defective mutations of HsTP lead to mitochondrial neurogastrointestinal encephalomyopathy (MNGIE), a disease with a high unmet medical need that is associated with severe neurological and gastrointestinal complications. Current efforts focus on the development of an enzyme replacement therapy (ERT) using the Escherichia coli ortholog (EcTP). However, bacterial enzymes are counter-indicated for human therapeutic applications because they are recognized as foreign by the human immune system, thereby eliciting adverse immune responses and raising significant safety and efficacy risks. Thus, it is critical to utilize the HsTP enzyme as starting scaffold for pre-clinical drug development, thus de-risking the safety concerns associated with the use of bacterial enzymes. However, HsTP expresses very poorly in E. coli, whereas its PEGylation, a crucial chemical modification for achieving long serum persistence of therapeutic enzymes, is highly inefficient and negatively affects its catalytic activity. Here we focused on the engineering of the recombinant expression profile of HsTP in E. coli cells, as well as on the optimization of its PEGylation efficiency aiming at the development of an alternative therapeutic approach for MNGIE. We show that phylogenetic and structural analysis of proteins can provide important insights for the rational design of N’-terminus-truncation constructs which exhibit significantly improved recombinant expression levels. In addition, we developed and implemented a criteria-driven rational surface engineering strategy for the substitution of arginine-to-lysine and lysine-to-arginine residues to achieve more efficient, homogeneous and reproducible PEGylation without negatively affecting the enzymatic catalytic activity upon PEGylation. Collectively, our proposed strategies provide an effective way to optimize enzyme PEGylation and E. coli recombinant expression and are likely applicable for other proteins and enzymes.

Kinetic isotope effects reveal early transition state of protein lysine methyltransferase SET8

Protein lysine methyltransferases (PKMTs) catalyze the methylation of protein substrates, and their dysregulation has been linked to many diseases, including cancer. Accumulated evidence suggests that the reaction path of PKMT-catalyzed methylation consists of the formation of a cofactor(cosubstrate)-PKMT-substrate complex, lysine deprotonation through dynamic water channels, and a nucleophilic substitution (SN2) transition state for transmethylation. However, the molecular characters of the proposed process remain to be elucidated experimentally. Here we developed a matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS) method and corresponding mathematic matrix to determine precisely the ratios of isotopically methylated peptides. This approach may be generally applicable for examining the kinetic isotope effects (KIEs) of posttranslational modifying enzymes. Protein lysine methyltransferase SET8 is the sole PKMT to monomethylate histone 4 lysine 20 (H4K20) and its function has been implicated in normal cell cycle progression and cancer metastasis. We therefore implemented the MS-based method to measure KIEs and binding isotope effects (BIEs) of the cofactor S-adenosyl-l-methionine (SAM) for SET8-catalyzed H4K20 monomethylation. A primary intrinsic 13C KIE of 1.04, an inverse intrinsic α-secondary CD3 KIE of 0.90, and a small but statistically significant inverse CD3 BIE of 0.96, in combination with computational modeling, revealed that SET8-catalyzed methylation proceeds through an early, asymmetrical SN2 transition state with the C-N and C-S distances of 2.35-2.40 Å and 2.00-2.05 Å, respectively. This transition state is further supported by the KIEs, BIEs, and steady-state kinetics with the SAM analog Se-adenosyl-l-selenomethionine (SeAM) as a cofactor surrogate. The distinct transition states between protein methyltransferases present the opportunity to design selective transition-state analog inhibitors.

Mechanisms for enzymatic cleavage of the N-glycosidic bond in DNA

DNA glycosylases remove damaged or enzymatically modified nucleobases from DNA, thereby initiating the base excision repair (BER) pathway, which is found in all forms of life. These ubiquitous enzymes promote genomic integrity by initiating repair of mutagenic and/or cytotoxic lesions that arise continuously due to alkylation, deamination, or oxidation of the normal bases in DNA. Glycosylases also perform essential roles in epigenetic regulation of gene expression, by targeting enzymatically-modified forms of the canonical DNA bases. Monofunctional DNA glycosylases hydrolyze the N-glycosidic bond to liberate the target base, while bifunctional glycosylases mediate glycosyl transfer using an amine group of the enzyme, generating a Schiff base intermediate that facilitates their second activity, cleavage of the DNA backbone. Here we review recent advances in understanding the chemical mechanism of monofunctional DNA glycosylases, with an emphasis on how the reactions are influenced by the properties of the nucleobase leaving-group, the moiety that varies across the vast range of substrates targeted by these enzymes.

Alternative substrates reveal catalytic cycle and key binding events in the reaction catalysed by anthranilate phosphoribosyltransferase from Mycobacterium tuberculosis

AnPRT (anthranilate phosphoribosyltransferase), required for the biosynthesis of tryptophan, is essential for the virulence of Mycobacterium tuberculosis (Mtb). AnPRT catalyses the Mg2+-dependent transfer of a phosphoribosyl group from PRPP (5'-phosphoribosyl-1'-pyrophosphate) to anthranilate to form PRA (5'-phosphoribosyl anthranilate). Mtb-AnPRT was shown to catalyse a sequential reaction and significant substrate inhibition by anthranilate was observed. Antimycobacterial fluoroanthranilates and methyl-substituted analogues were shown to act as alternative substrates for Mtb-AnPRT, producing the corresponding substituted PRA products. Structures of the enzyme complexed with anthranilate analogues reveal two distinct binding sites for anthranilate. One site is located over 8 Å (1 Å=0.1 nm) from PRPP at the entrance to a tunnel leading to the active site, whereas in the second, inner, site anthranilate is adjacent to PRPP, in a catalytically relevant position. Soaking the analogues for variable periods of time provides evidence for anthranilate located at transient positions during transfer from the outer site to the inner catalytic site. PRPP and Mg2+ binding have been shown to be associated with the rearrangement of two flexible loops, which is required to complete the inner anthranilate-binding site. It is proposed that anthranilate first binds to the outer site, providing an unusual mechanism for substrate capture and efficient transfer to the catalytic site following the binding of PRPP.

In vitro selection of DNA-cleaving deoxyribozyme with site-specific thymidine excision activity

Single-nucleotide polymorphisms, either inherited or due to spontaneous DNA damage, are associated with numerous diseases. Developing tools for site-specific nucleotide modification may one day provide a way to alter disease polymorphisms. Here, we describe the in vitro selection and characterization of a new deoxyribozyme called F-8, which catalyzes nucleotide excision specifically at thymidine. Cleavage by F-8 generates 3'- and 5'-phosphate ends recognized by DNA modifying enzymes, which repair the targeted deoxyribonucleotide while maintaining the integrity of the rest of the sequence. These results illustrate the potential of DNAzymes as tools for DNA manipulation.

The kinetic mechanism of Human Thymidine Phosphorylase - a molecular target for cancer drug development

Human Thymidine Phosphorylase (HTP), also known as the platelet-derived endothelial cell growth factor (PD-ECGF) or gliostatin, catalyzes the reversible phosphorolysis of thymidine (dThd) to thymine and 2-deoxy-α-d-ribose-1-phosphate (2dR1P). HTP is a key enzyme in the pyrimidine salvage pathway involved in dThd homeostasis in cells. HTP is a target for anticancer drug development as its enzymatic activity promotes angiogenesis. Here, we describe cloning, expression, and purification to homogeneity of recombinant TYMP-encoded HTP. Peptide fingerprinting and the molecular mass value of the homogenous protein confirmed its identity as HTP assessed by mass spectrometry. Size exclusion chromatography showed that HTP is a dimer in solution. Kinetic studies revealed that HTP displayed substrate inhibition for dThd. Initial velocity and isothermal titration calorimetry (ITC) studies suggest that HTP catalysis follows a rapid-equilibrium random bi-bi kinetic mechanism. ITC measurements also showed that dThd and Pi binding are favorable processes. The pH-rate profiles indicated that maximal enzyme activity was achieved at low pH values. Functional groups with apparent pK values of 5.2 and 9.0 are involved in dThd binding and groups with pK values of 6.1 and 7.8 are involved in phosphate binding.

Glycosidic bond cleavage in deoxynucleotides: effects of solvent and the DNA phosphate backbone in the computational model

Density functional theory (B3LYP) was employed to examine the hydrolysis of the canonical 2'-deoxynucleotides in varied environments (gas phase or water) using different computational models for the sugar residue (methyl or phosphate group at C5') and nucleophile (water activated through full or partial proton abstraction). Regardless of the degree of nucleophile activation, our results show that key geometrical parameters along the reaction pathway are notably altered upon direct inclusion of solvent effects in the optimization routine, which leads to significant changes in the reaction energetics and better agreement with experiment. Therefore, despite the wide use of gas-phase calculations in the literature, small model computational work, as well as large-scale enzyme models, that strive to understand nucleotide deglycosylation must adequately describe the environment. Alternatively, although inclusion of the phosphate group at C5' also affects the geometries of important stationary points, the effects cancel to yield unchanged deglycosylation barriers, and therefore smaller computational models can be used to estimate the energy associated with nucleotide deglycosylation, with the 5' phosphate group included if full (geometric) details of the reaction are desired. Hydrogen-bonding interactions with the nucleobase can significantly reduce the barrier to deglycosylation, which supports suggestions that discrete hydrogen-bonding interactions with active-site amino acid residues can play a significant role in enzyme-catalyzed nucleobase excision. Taken together with previous studies, the present work provides vital clues about the components that must be included in future studies of the deglycosylation of isolated noncanonical nucleotides, as well as the corresponding enzyme-catalyzed reactions.

Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines

Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(‡) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.

Quantum theoretical study of cleavage of the glycosidic bond of 2'-deoxyadenosine: base excision-repair mechanism of DNA by MutY

The enzyme adenine DNA glycosylase, also called MutY, is known to catalyze base excision repair by removal of adenine from the abnormal 2'-deoxyadenosine:8-oxo-2'-deoxyguanosine pair in DNA. The active site of the enzyme was considered to consist of a glutamic acid residue along with two water molecules. The relevant reaction mechanism involving different barrier energies was studied theoretically. Molecular geometries of the various molecules and complexes involved in the reaction, e.g., the reactant, intermediate, and product complexes as well as transition states, were optimized employing density functional theory at the B3LYP/6-31G(d,p) level in the gas phase. It was followed by single-point energy calculations at the B3LYP/AUG-cc-pVDZ, BHandHLYP/AUG-cc-pVDZ, and MP2/AUG-cc-pVDZ levels in the gas phase. Single-point energy calculations were also carried out at the B3LYP/AUG-cc-pVDZ and BHandHLYP/AUG-cc-pVDZ levels in aqueous media as well as in the solvents chlorobenzene and dichloroethane. For the solvation calculations, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. It is found that glutamic acid along with two water molecules would effectively cleave the glycosidic bond of adenosine by a new two-step reaction mechanism proposed here which is different from the three-step mechanism proposed by other authors earlier regarding the working mechanism of MutY.

Transition state analysis of the arsenolytic depyrimidination of thymidine by human thymidine phosphorylase

Human thymidine phosphorylase (hTP) is responsible for thymidine (dT) homeostasis, promotes angiogenesis, and is involved in metabolic inactivation of antiproliferative agents that inhibit thymidylate synthase. Understanding its transition state structure is on the path to design transition state analogues. Arsenolysis of dT by hTP permits kinetic isotope effect (KIE) analysis of the reaction by forming thymine and the chemically unstable 2-deoxyribose 1-arsenate. The transition state for the arsenolytic reaction was characterized using multiple KIEs and computational analysis. Transition state analysis revealed a concerted bimolecular (A(N)D(N)) mechanism. A transition state constrained to match the intrinsic KIE values was found using density functional theory (B3LYP/6-31G*). An active site histidine is implicated as the catalytic base responsible for activation of the arsenate nucleophile and stabilization of the thymine leaving group during the isotopically sensitive step. At the transition state, the deoxyribose ring exhibits significant oxocarbenium ion character with bond breaking (r(C-N) = 2.45 Å) nearly complete and minimal bond making to the attacking nucleophile (r(C-O) = 2.95 Å). The transition state model predicts a deoxyribose conformation with a 2'-endo ring geometry. Transition state structure for the slow hydrolytic reaction of hTP involves a stepwise mechanism [Schwartz, P. A., Vetticatt, M. J., and Schramm, V. L. (2010) J. Am. Chem. Soc. 132, 13425-13433], in contrast to the concerted mechanism described here for arsenolysis.

Transition state analysis of thymidine hydrolysis by human thymidine phosphorylase

Human thymidine phosphorylase (hTP) is responsible for thymidine (dT) homeostasis, and its action promotes angiogenesis. In the absence of phosphate, hTP catalyzes a slow hydrolytic depyrimidination of dT yielding thymine and 2-deoxyribose (dRib). Its transition state was characterized using multiple kinetic isotope effect (KIE) measurements. Isotopically enriched thymidines were synthesized enzymatically from glucose or (deoxy)ribose, and intrinsic KIEs were used to interpret the transition state structure. KIEs from [1'-(14)C]-, [1-(15)N]-, [1'-(3)H]-, [2'R-(3)H]-, [2'S-(3)H]-, [4'-(3)H]-, and [5'-(3)H]dTs provided values of 1.033 ± 0.002, 1.004 ± 0.002, 1.325 ± 0.003, 1.101 ± 0.004, 1.087 ± 0.005, 1.040 ± 0.003, and 1.033 ± 0.003, respectively. Transition state analysis revealed a stepwise mechanism with a 2-deoxyribocation formed early and a higher energetic barrier for nucleophilic attack of a water molecule on the high energy intermediate. An equilibrium exists between the deoxyribocation and reactants prior to the irreversible nucleophilic attack by water. The results establish activation of the thymine leaving group without requirement for phosphate. A transition state constrained to match the intrinsic KIEs was found using density functional theory. An active site histidine (His116) is implicated as the catalytic base for activation of the water nucleophile at the rate-limiting transition state. The distance between the water nucleophile and the anomeric carbon (r(C-O)) is predicted to be 2.3 A at the transition state. The transition state model predicts that deoxyribose adopts a mild 3'-endo conformation during nucleophilic capture. These results differ from the concerted bimolecular mechanism reported for the arsenolytic reaction (Birck, M. R.; Schramm, V. L. J. Am. Chem. Soc. 2004, 126, 2447-2453).

Transition state of ADP-ribosylation of acetyllysine catalyzed by Archaeoglobus fulgidus Sir2 determined by kinetic isotope effects and computational approaches

Sirtuins are protein-modifying enzymes distributed throughout all forms of life. These enzymes bind NAD(+), a universal metabolite, and react it with acetyllysine residues to effect deacetylation of protein side chains. This NAD(+)-dependent deacetylation reaction has been observed for sirtuin enzymes derived from archaeal, eubacterial, yeast, metazoan, and mammalian species, suggesting conserved chemical mechanisms for these enzymes. The first chemical step of deacetylation is the reaction of NAD(+) with an acetyllysine residue which forms an enzyme-bound ADPR-peptidylimidate intermediate and nicotinamide. In this manuscript, the transition state for the ADP-ribosylation of acetyllysine is solved for an Archaeoglobus fulgidus sirtuin (Af2Sir2). Kinetic isotope effects (KIEs) were obtained by the competitive substrate method and were [1(N)-(15)N] = 1.024(2), [1'(N)-(14)C] = 1.014(4), [1'(N)-(3)H] = 1.300(3), [2'(N)-(3)H] = 1.099(5), [4'(N)-(3)H] = 0.997(2), [5'(N)-(3)H] = 1.020(5), [4'(N)-(18)O] = 0.984(5). KIEs were calculated for candidate transition state structures using computational methods (Gaussian 03 and ISOEFF 98) in order to match computed and experimentally determined KIEs to solve the transition state. The results indicate that the enzyme stabilizes a highly dissociated oxocarbenium ionlike transition state with very low bond orders to the leaving group nicotinamide and the nucleophile acetyllysine. A concerted yet highly asynchronous substitution mechanism forms the ADPR-peptidylimidate intermediate of the sirtuin deacetylation reaction.

Effects of nucleophile, oxidative damage, and nucleobase orientation on the glycosidic bond cleavage in deoxyguanosine

Deglycosylation of nucleotides occurs during many essential biological processes, including DNA repair, and is initiated by a variety of nucleophiles. In the present work, density functional theory (B3LYP) was used to investigate the thermodynamics and kinetics of the glycosidic bond cleavage reaction in the model nucleoside forms of guanine and its major oxidation product, 8-oxoguanine. Base excision facilitated by four different nucleophiles (hydroxyl anion (fully activated water), formate-water complex (partially activated water), lysine, and proline) was considered, which spans nucleophiles involved in a collection of spontaneous and enzyme-catalyzed processes. Because some enzymes that catalyze deglycosylation can accommodate more than one orientation of the base with respect to the sugar moiety, the effects of the (anti/syn) base orientation on the barrier height were also considered. We find that the nucleophile has a very large effect on the overall (gas-phase) reaction energetics. Although this effect decreases in different (polar) environments, the nucleophile has the greatest influence on the overall reaction as compared to whether the base is damaged or to the base orientation. Furthermore, the effects are significant in environments that most closely resemble (nonpolar) enzymatic active sites. Our results provide a greater understanding of the relative effects of the nucleophile, damage to the nucleobase, and the nucleobase orientation with respect to the sugar moiety on the deglycosylation pathway, which provide qualitative explanations for relative base excision rates observed in some biological systems.

Modeling the dissociative hydrolysis of the natural DNA nucleosides

Two-dimensional PCM-B3LYP/6-31+G(d) potential energy surfaces for the hydrolysis of the four natural 2'-deoxyribonucleosides (2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and thymidine) are characterized using a model that includes both implicit (bulk) solvent effects and (three or four) explicit water molecules in the optimization routine. For the first time, the experimentally predicted dissociative (S(N)1) mechanism is found to be favored over the synchronous (S(N)2) pathway for all nucleosides studied. Due to the success of our model in stabilizing the charge-separated intermediates along the S(N)1 pathway, it is proposed that the new model presented here is the smallest system capable of generating the experimentally predicted oxacarbenium cation intermediate. We therefore stress that dissociative mechanisms should be studied with methodologies that account for the (bulk) environment in the optimization routine, where these effects are often only included as a correction to the energy in the current literature. In addition to accounting for charge stabilization through implicit solvation, nucleophile activation and leaving group stabilization should also be explicitly introduced into the model to further stabilize the system. Our work also emphasizes the importance of studying the Gibbs surface, which in some cases provides a better description of chemically important regions of the reaction surface or changes the calculated trend in the magnitude of dissociative barriers. In addition, it is proposed that the methodology presented in this study can be used to calculate uncatalyzed deglycosylation barriers for a range of DNA nucleosides, which when compared to the corresponding enzyme-catalyzed reactions, will allow the prediction of the rate enhancement (barrier reduction) due to the enzyme.

The role of phosphate in the action of thymidine phosphorylase inhibitors: Implications for the catalytic mechanism

The design and synthesis of 5-fluoro-6-[(2-aminoimidazol-1-yl)methyl]uracil (AIFU), a potent inhibitor of thymidine phosphorylase (TP) with K(i)-values of 11nM (ecTP) and 17nM (hTP), are described. Kinetic studies established that the type of inhibition of TP by AIFU is uncompetitive with respect to inorganic phosphate (or arsenate). The results obtained suggest that AIFU and other zwitterionic thymine analog inhibitors of TP act as transition state analogs, mimicking the anionic thymine leaving group, consistent with an S(N)2-type catalytic mechanism, and anchored by their protonated side chains to the enzyme-bound phosphate by electrostatic and H-bonding interactions.

Glycosidic bond cleavage in deoxynucleotides — A density functional study

Density functional theory was used to study the glycosidic bond cleavage in deoxynucleotides with the main goal to determine the effects of the nucleobase, hydrogen bonding with the nucleobase, and the (bulk) environment on the reaction energetics. Since direct glycosidic bond cleavage is a high-energy process, two nucleophile models were considered (HCOO–···H2O and HO–), which represent different stages of activation of a water nucleophile. The glycosidic bond cleavage barriers were found to decrease, while the reaction exothermicity increases, with an increase in the nucleobase acidity. The gas-phase barriers and reaction energies for bond cleavage in all deoxynucleotides were found to be significantly affected by hydrogen-bonding interactions with the nucleobase (by up to 30 kJ mol–1 depending on the nucleophile). Although the barriers increase and reaction energies become less exothermic in enzymatic and aqueous environments, the effects of the bulk environment are similar in the presence and absence of small molecules bound to the nucleobase. Therefore, the effects of hydrogen bonding with the bases are approximately the same in all environments. Our results suggest that hydrogen bonding with the nucleobase may play an important role in the glycosidic bond cleavage in both pyrimidine and purine nucleotides in a variety of environments.

A kinetic, modeling and mechanistic re-analysis of thymidine phosphorylase and some related enzymes

Thymidine phosphorylase (TP) is an important target enzyme for cancer chemotherapy but currently available inhibitors lack in vivo potency. Related enzymes also are therapeutic targets. A greater understanding of enzyme structure and mechanism may help in the design of improved drugs and this work assists in that regard. Also important is the correct identification of the ionization states and tautomeric forms of substrates and products when bound to the enzyme and during the course of the reaction. Approximate methods for estimating some deltapK(a)s between aqueous and protein-bound substrates are exemplified for nucleobases and nucleosides. The estimates demonstrate that carbonyl-protonated thymidine and hydroxy tautomers of thymine are not involved in TP's actions. Other estimates indicate that purine nucleoside phosphorylase binds inosine and guanosine as zwitterionic tautomers and that phosphorolysis proceeds through these forms. Extensive molecular modeling based on an X-ray structure of human TP indicates that TP is likely to be mechanistically similar to all other natural members of the class in proceeding through a alpha-oxacarbenium-like transition state or states.

Thymidine phosphorylase fromEscherichia coli: Tight-binding inhibitors as enzyme active-site titrants

Thymidine phosphorylase (EC 2.4.2.4) catalyses the reversible phosphorolysis of pyrimidine 2'-deoxynucleosides, forming 2-deoxyribose-1-phosphate and pyrimidine. 5-Chloro-6-(2-imino-pyrrolidin-1-yl)methyl-uracil hydrochloride (TPI, 1) and its 5-bromo analogue (2), 6-(2-amino-imidazol-1-yl)methyl-5-bromo-uracil (3) and its 5-chloro analogue (4) act as tight-binding stoichiometric inhibitors of recombinant E. coli thymidine phosphorylase, and thus can be used as the first active-site titrants for it using either thymidine or 5-nitro-2'-deoxyuridine as substrate.

Transition-state analysis of the DNA repair enzyme MutY

The transition state (TS) structure of MutY-catalyzed DNA hydrolysis was solved using multiple kinetic isotope effect (KIE) measurements. MutY is a base excision repair enzyme which cleaves adenine from 8-oxo-G:A mismatches in vivo, and also from G:A mismatches in vitro. TS analysis of G:A-DNA hydrolysis revealed a stepwise S(N)1 (D(N)*A(N)(double dagger)) mechanism proceeding through a highly reactive oxacarbenium ion intermediate which would have a lifetime in solution of <10(-10) s. C-N bond cleavage is reversible; the N-glycoside bond breaks and reforms repeatedly before irreversible water attack on the oxacarbenium ion. KIEs demonstrated that MutY uses general acid catalysis by protonating N7. It enforces a 3'-exo sugar ring conformation and other sugar ring distortions to stabilize the oxacarbenium ion. Combining the experimental TS structure with the previously reported crystal structure of an abortive Michaelis complex elucidates the step-by-step catalytic sequence.

Picomolar inhibitors as transition-state probes of 5'-methylthioadenosine nucleosidases

Transition states can be predicted from an enzyme's affinity to related transition-state analogues. 5'-Methylthioadenosine nucleosidases (MTANs) are involved in bacterial quorum sensing pathways and thus are targets for antibacterial drug design. The transition-state characteristics of six MTANs are compared by analyzing dissociation constants (K(d)) with a small array of representative transition-state analogues. These inhibitors mimic early or late dissociative transition states with K(d) values in the picomolar range. Our results indicate that the K(d) ratio for mimics of early and late transition states are useful in distinguishing between these states. By this criterion, the transition states of Neisseria meningitides and Helicobacter pylori MTANs are early dissociative, whereas Escherichia coli, Staphylococcus aureus, Streptococcus pneumoniae, and Klebsiella pneumoniae MTANs have late dissociative characters. This conclusion is confirmed independently by the characteristic [1'- (3)H] and [1'- (14)C] kinetic isotope effects (KIEs) of these enzymes. Large [1'- (3)H] and unity [1'- (14)C] KIEs are observed for late dissociative transition states, whereas early dissociative states showed close-to-unity [1'- (3)H] and significant [1'- (14)C] KIEs. K d values of various MTANs for individual transition-state analogues provide tentative information about transition-state structures due to varying catalytic efficiencies of enzymes. Comparing K d ratios for mimics of early and late transition states removes limitations inherent to the enzyme and provides a better predictive tool in discriminating between possible transition-state structures.

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.

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.

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.

Transition-state analysis of S. pneumoniae 5'-methylthioadenosine nucleosidase

Kinetic isotope effects (KIEs) and computer modeling are used to approximate the transition state of S. pneumoniae 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN). Experimental KIEs were measured and corrected for a small forward commitment factor. Intrinsic KIEs were obtained for [1'-3H], [1'-14C], [2'-3H], [4'-3H], [5'-3H(2)], [9-15N] and [Me-3H(3)] MTAs. The intrinsic KIEs suggest an SN1 transition state with no covalent participation of the adenine or the water nucleophile. The transition state was modeled as a stable ribooxacarbenium ion intermediate and was constrained to fit the intrinsic KIEs. The isotope effects predicted a 3-endo conformation for the ribosyl oxacarbenium-ion corresponding to H1'-C1'-C2'-H2' dihedral angle of 70 degrees. Ab initio Hartree-Fock and DFT calculations were performed to study the effect of polarization of ribosyl hydroxyls, torsional angles, and the effect of base orientation on isotope effects. Calculations suggest that the 4'-3H KIE arises from hyperconjugation between the lonepair (n(p)) of O4' and the sigma* (C4'-H4') antibonding orbital owing to polarization of the 3'-hydroxyl by Glu174. A [methyl-3H(3)] KIE is due to hyperconjugation between np of sulfur and sigma* of methyl C-H bonds. The van der Waal contacts increase the 1'-3H KIE because of induced dipole-dipole interactions. The 1'-3H KIE is also influenced by the torsion angles of adjacent atoms and by polarization of the 2'-hydroxyl. Changing the virtual solvent (dielectric constant) does not influence the isotope effects. Unlike most N-ribosyltransferases, N7 of the leaving group adenine is not protonated at the transition state of S. pneumoniae MTAN. This feature differentiates the S. pneumoniae and E. coli transition states and explains the 10(3)-fold decrease in the catalytic efficiency of S. pneumoniae MTAN relative to that from E. coli.

Transition-state structure of human 5'-methylthioadenosine phosphorylase

Kinetic isotope effects (KIEs) and computer modeling using density functional theory were used to approximate the transition state of human 5'-methylthioadenosine phosphorylase (MTAP). KIEs were measured on the arsenolysis of 5'-methylthioadenosine (MTA) catalyzed by MTAP and were corrected for the forward commitment to catalysis. Intrinsic KIEs were obtained for [1'-(3)H], [1'-(14)C], [2'-(3)H], [4'-(3)H], [5'-(3)H(2)], [9-(15)N], and [Me-(3)H(3)] MTAs. The primary intrinsic KIEs (1'-(14)C and 9-(15)N) suggest that MTAP has a dissociative S(N)1 transition state with its cationic center at the anomeric carbon and insignificant bond order to the leaving group. The 9-(15)N intrinsic KIE of 1.039 also establishes an anionic character for the adenine leaving group, whereas the alpha-primary 1'-(14)C KIE of 1.031 indicates significant nucleophilic participation at the transition state. Computational matching of the calculated EIEs to the intrinsic isotope effects places the oxygen nucleophile 2.0 Angstrom from the anomeric carbon. The 4'-(3)H KIE is sensitive to the polarization of the 3'-OH group. Calculations suggest that a 4'-(3)H KIE of 1.047 is consistent with ionization of the 3'-OH group, indicating formation of a zwitterion at the transition state. The transition state has cationic character at the anomeric carbon and is anionic at the 3'-OH oxygen, with an anionic leaving group. The isotope effects predicted a 3'-endo conformation for the ribosyl zwitterion, corresponding to a H1'-C1'-C2'-H2' torsional angle of 33 degrees. The [Me-(3)H(3)] and [5'-(3)H(2)] KIEs arise predominantly from the negative hyperconjugation of the lone pairs of sulfur with the sigma (C-H) antibonding orbitals. Human MTAP is characterized by a late S(N)1 transition state with significant participation of the phosphate nucleophile.

6-substituted 5-fluorouracil derivatives as transition state analogue inhibitors of thymidine phosphorylase

A combination of mechanism-based and structure-based design strategies led to the synthesis of a series of 5- and 6-substituted uracil derivatives as potential inhibitors of thymidine phosphorlase/platelet derived endothelial cell growth factor (TP/PD-ECGF). Among those tested, 6-imidazolylmethyl-5-fluorouracil was found to be the most potent inhibitor with a Ki-value of 51 nM, representing a new class of 5-fluoropyrimidines with a novel mechanism of action.

Kinetic isotope effects of nucleoside hydrolase from Escherichia coli

rihC is one of a group of three ribonucleoside hydrolases found in Escherichia coli (E. coli). The enzyme catalyzes the hydrolysis of selected nucleosides to ribose and the corresponding base. A family of Vmax/Km kinetic isotope effects using uridine labeled with stable isotopes, such as 2H, 13C, and 15N, were determined by liquid chromatography/mass spectrometry (LC/MS). The kinetic isotope effects were 1.012+/-0.006, 1.027+/-0.005, 1.134+/-0.007, 1.122+/-0.008, and 1.002+/-0.004 for [1'-13C], [1-15N], [1'-2H], [2'-2H], and [5'-2H2] uridine, respectively. A transition state based upon a bond-energy bond-order vibrational analysis (BEBOVIB) of the observed kinetic isotope effects is proposed. The main features of this transition state are activation of the heterocyclic base by protonation of/or hydrogen bonding to O2, an extensively broken C-N glycosidic bond, formation of an oxocarbenium ion in the ribose ring, C3'-exo ribose ring conformation, and almost no bond formation to the attacking nucleophile. The proposed transition state for the prokaryotic E. coli nucleoside hydrolase is compared to that of a similar enzyme isolated from Crithidia fasciculata (C. fasciculata).

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.

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.

The Synthesis of N-Ribosyl Transferase Inhibitors Based on a Transition State Blueprint

A quarter of a century ago transition state analysis and transition state analogue design promised the prospect of extraordinarily potent enzyme inhibitors. The present overview describes the transition state analysis of a variety of N-ribosyl transferases, the design and synthesis of extremely powerful transition state analogue inhibitors of these nucleoside processing enzymes, and their current therapeutic uses and potentials.