Ricin A-Chain: Kinetic Isotope Effects and Transition State Structure with Stem-Loop RNA†

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.

Oligonucleotide transition state analogues of saporin L3

Ribosome inactivating proteins (RIPs) are among the most toxic agents known. More than a dozen clinical trials against refractory cancers have been initiated using modified RIPs with impressive results. However, dose-limiting toxicity due to vascular leak syndrome limits success of the therapy. We have previously reported some tight-binding transition state analogues of Saporin L3 that mimic small oligonucleotide substrates in which the susceptible adenosine has been replaced by a 9-deazaadenyl hydroxypyrrolidinol derivative. They provide the first step in the development of rescue agents to prevent Saporin L3 toxicity on non-targeted cells. Here we report the synthesis, using solution phase chemistry, of these and a larger group of transition state analogues. They were tested for inhibition against Saporin L3 giving K_i values as low as 3.3 nM and indicating the structural requirements for inhibition.

Transition State Structure of RNA Depurination by Saporin L3

Saporin L3 from the leaves of the common soapwort is a catalyst for hydrolytic depurination of adenine from RNA. Saporin L3 is a type 1 ribosome inactivating protein (RIP) composed only of a catalytic domain. Other RIPs have been used in immunotoxin cancer therapy, but off-target effects have limited their development. In the current study, we use transition state theory to understand the chemical mechanism and transition state structure of saporin L3. In favorable cases, transition state structures guide the design of transition state analogues as inhibitors. Kinetic isotope effects (KIEs) were determined for an A14C mutant of saporin L3. To permit KIE measurements, small stem-loop RNAs that contain an AGGG tetraloop structure were enzymatically synthesized with the single adenylate bearing specific isotopic substitutions. KIEs were measured and corrected for forward commitment to obtain intrinsic values. A model of the transition state structure for depurination of stem-loop RNA (5'-GGGAGGGCCC-3') by saporin L3 was determined by matching KIE values predicted via quantum chemical calculations to a family of intrinsic KIEs. This model indicates saporin L3 displays a late transition state with the N-ribosidic bond to the adenine nearly cleaved, and the attacking water nucleophile weakly bonded to the ribosyl anomeric carbon. The transition state retains partial ribocation character, a feature common to most N-ribosyl transferases. However, the transition state geometry for saporin L3 is distinct from ricin A-chain, the only other RIP whose transition state is known.

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.

Probing the catalytic mechanism of bovine CD38/NAD+ glycohydrolase by site directed mutagenesis of key active site residues

Bovine CD38/NAD(+) glycohydrolase catalyzes the hydrolysis of NAD(+) to nicotinamide and ADP-ribose and the formation of cyclic ADP-ribose via a stepwise reaction mechanism. Our recent crystallographic study of its Michaelis complex and covalently-trapped intermediates provided insights into the modalities of substrate binding and the molecular mechanism of bCD38. The aim of the present work was to determine the precise role of key conserved active site residues (Trp118, Glu138, Asp147, Trp181 and Glu218) by focusing mainly on the cleavage of the nicotinamide-ribosyl bond. We analyzed the kinetic parameters of mutants of these residues which reside within the bCD38 subdomain in the vicinity of the scissile bond of bound NAD(+). To address the reaction mechanism we also performed chemical rescue experiments with neutral (methanol) and ionic (azide, formate) nucleophiles. The crucial role of Glu218, which orients the substrate for cleavage by interacting with the N-ribosyl 2'-OH group of NAD(+), was highlighted. This contribution to catalysis accounts for almost half of the reaction energy barrier. Other contributions can be ascribed notably to Glu138 and Asp147 via ground-state destabilization and desolvation in the vicinity of the scissile bond. Key interactions with Trp118 and Trp181 were also proven to stabilize the ribooxocarbenium ion-like transition state. Altogether we propose that, as an alternative to a covalent acylal reaction intermediate with Glu218, catalysis by bCD38 proceeds through the formation of a discrete and transient ribooxocarbenium intermediate which is stabilized within the active site mostly by electrostatic interactions.

Divergent mechanisms for enzymatic excision of 5-formylcytosine and 5-carboxylcytosine from DNA

5-Methylcytosine (mC) is an epigenetic mark that impacts transcription, development, and genome stability, and aberrant DNA methylation contributes to aging and cancer. Active DNA demethylation involves stepwise oxidation of mC to 5-hydroxymethylcytosine, 5-formylcytosine (fC), and potentially 5-carboxylcytosine (caC), excision of fC or caC by thymine DNA glycosylase (TDG), and restoration of cytosine via follow-on base excision repair. Here, we investigate the mechanism for TDG excision of fC and caC. We find that 5-carboxyl-2'-deoxycytidine ionizes with pK(a) values of 4.28 (N3) and 2.45 (carboxyl), confirming that caC exists as a monoanion at physiological pH. Calculations do not support the proposal that G·fC and G·caC base pairs adopt a wobble structure that is recognized by TDG. Previous studies show that N-glycosidic bond hydrolysis follows a stepwise (S(N)1) mechanism, and that TDG activity increases with pyrimidine N1 acidity, that is, leaving group quality of the target base. Calculations here show that fC and the neutral tautomers of caC are acidic relative to other TDG substrates, but the caC monoanion exhibits poor acidity and likely resists TDG excision. While fC activity is independent of pH, caC excision is acid-catalyzed, and the pH profile indicates that caC ionizes in the enzyme-substrate complex with an apparent pKa of 5.8, likely at N3. Mutational analysis reveals that Asn191 is essential for excision of caC but dispensable for fC activity, indicating that N191 may stabilize N3-protonated forms of caC to facilitate acid catalysis and suggesting that N191A-TDG could potentially be useful for studying DNA demethylation in cells.

Depurination of N7-methylguanine by DNA glycosylase AlkD is dependent on the DNA backbone

DNA glycosylase AlkD excises N7-methylguanine (7mG) by a unique but unknown mechanism, in which the damaged nucleotide is positioned away from the protein and the phosphate backbone is distorted. Here, we show by methylphosphonate substitution that a phosphate proximal to the lesion has a significant effect on the rate enhancement of 7mG depurination by the enzyme. Thus, instead of a conventional mechanism whereby protein side chains participate in N-glycosidic bond cleavage, AlkD remodels the DNA into an active site composed exclusively of DNA functional groups that provide the necessary chemistry to catalyze depurination.

Nucleic acid bases in anionic 2'-deoxyribonucleotides: a DFT/B3LYP study of structures, relative stability, and proton affinities

Protonation of nucleobases in anions of canonical 2'-deoxyribonucleotides has been investigated by the DFT computational study at the B3LYP/aug-cc-pvdz level of theory. It is demonstrated that the protonation leads to a significant decrease of conformational space of purine nucleotides while almost all conformers found for non-protonated molecules correspond to minima of the potential energy surface for protonated mdTMP and mdCMP. However, in all nucleotides, only one conformer is populated. This applies to all tautomers of protonated molecules except the mdTMP and mdCMP with the proton attached to the carbonyl group where a minor population of second conformer is observed. Protonation of nucleobase leads to significant elongation of the N-glycosidic bond. These findings agree well with suggestions that protonation of nucleobase is a first step in cleavage of the glycosidic bond. The oxygen atoms of both carbonyl groups of thymine and the N3 atom of the pyrimidine ring of cytosine, guanine, and adenine represent the most preferable sites for protonation of anions of 2'-deoxyrobonucleotides. The highest proton affinity is observed for the base in mdGMP and the lowest for the thymine moiety in mdTMP. It should be noted that calculated values of the proton affinities in anionic nucleotides are significantly higher (by 2-3 eV) than for nucleosides and neutral nucleotides. This allows assuming that the proton affinity of the base in DNA macromolecule may be tuned by changing the extent of shielding or neutralization of negative charge of the phosphate group.

Transition state analysis of acid-catalyzed hydrolysis of an enol ether, enolpyruvylshikimate 3-phosphate (EPSP)

Proton transfer to carbon represents a significant catalytic challenge because of the large intrinsic energetic barrier and the frequently unfavorable thermodynamics. Multiple kinetic isotope effects (KIEs) were measured for acid-catalyzed hydrolysis of the enol ether functionality of enolpyruvylshikimate 3-phosphate (EPSP) as a nonenzymatic analog of the EPSP synthase (AroA) reaction. The large solvent deuterium KIE demonstrated that protonating C3 was the rate-limiting step, and the lack of solvent hydron exchange into EPSP demonstrated that protonation was irreversible. The reaction mechanism was stepwise, with C3, the methylene carbon, being protonated to form a discrete oxacarbenium ion intermediate before water attack at the cationic center, that is, an AH(‡)*AN (or AH(‡) + AN) mechanism. The calculated 3-(14)C and 3,3-(2)H2 KIEs varied as a function of the extent of proton transfer at the transition state, as reflected in the C3-H(+) bond order, nC3-H+. The calculated 3-(14)C KIE was a function primarily of C3 coupling with the movement of the transferring proton, as reflected in the reaction coordinate contribution ((light)ν(‡)/(heavy)ν(‡)), rather than of changes in bonding. Coupling was strongest in early and late transition states, where the reaction coordinate frequency was lower. The other calculated (14)C and (18)O KIEs were more sensitive to interactions with counterions and solvation in the model structures than nC3-H+. The KIEs revealed a moderately late transition state with significant oxacarbenium ion character and with a C3-H(+) bond order ≈0.6.

Role of environment for catalysis of the DNA repair enzyme MutY

Control of the N-glycosylase reaction by the DNA repair enzyme, MutY, entails the organization of solvent molecules. Classical molecular dynamics and QM/MM simulations were used to investigate the solvent and environment effects contributing to catalysis. Our findings suggest that the entire reaction is an energetically neutral process, in which the first step is rate determining, requiring protonation of adenine (N(7)) to initiate cleavage, and the second step is strongly exothermic, involving hydrolysis of an oxacarbenium ion intermediate. Although water molecules are catalytically active in both steps, the first step requires an entirely different level of solvent organization compared to the second. Needed to secure protonation at N(7), a long-term solvation pattern is established which facilitates the involvement of three out of the five structured water molecules in the active site. This persistent arrangement coordinates the catalytically active water molecules into prime positions to assist the proton transfer: (i) a water molecule frequently bridges the catalytic residues and (ii) the bridging water molecule is assisted by 1-2 other 'supporting' water molecules. To maintain this configuration, MutY, surprisingly, uses hydrophobic residues in combination with hydrophilic residues to tune the microenvironment into a 'water trap'. Hydrophilic residues prolong solvent residence times by maintaining hydrogen-bonding networks, whereas the hydrophobic residues constrain the positioning of the catalytic water molecules that assist the proton-transfer event. In this way, the enzyme uses both entropic and enthalpic considerations to guide the water-assisted reaction.

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.

7-Substituted pterins provide a new direction for ricin A chain inhibitors

Ricin is a potent toxin found in castor seeds. The A chain, RTA, enzymaticlly depurinates a specific adenosine in ribosomal RNA, inhibiting protein synthesis. Ricin is a known chemical weapons threat having no effective antidote. This makes the discovery of new inhibitors of great importance. We have previously used 6-substituted pterins, such as pteroic acid, as an inhibitor platform with moderate success. We now report the success of 7-carboxy pterin (7CP) as an RTA inhibitor; its binding has been monitored using both kinetic and temperature shift assays and by X-ray crystallography. We also discuss the synthesis of various derivatives of 7CP, and their binding affinity and inhibitory effects, as part of a program to make effective RTA inhibitors.

Small-molecule inhibitor leads of ribosome-inactivating proteins developed using the doorstop approach

Ribosome-inactivating proteins (RIPs) are toxic because they bind to 28S rRNA and depurinate a specific adenine residue from the α-sarcin/ricin loop (SRL), thereby inhibiting protein synthesis. Shiga-like toxins (Stx1 and Stx2), produced by Escherichia coli, are RIPs that cause outbreaks of foodborne diseases with significant morbidity and mortality. Ricin, produced by the castor bean plant, is another RIP lethal to mammals. Currently, no US Food and Drug Administration-approved vaccines nor therapeutics exist to protect against ricin, Shiga-like toxins, or other RIPs. Development of effective small-molecule RIP inhibitors as therapeutics is challenging because strong electrostatic interactions at the RIP•SRL interface make drug-like molecules ineffective in competing with the rRNA for binding to RIPs. Herein, we report small molecules that show up to 20% cell protection against ricin or Stx2 at a drug concentration of 300 nM. These molecules were discovered using the doorstop approach, a new approach to protein•polynucleotide inhibitors that identifies small molecules as doorstops to prevent an active-site residue of an RIP (e.g., Tyr80 of ricin or Tyr77 of Stx2) from adopting an active conformation thereby blocking the function of the protein rather than contenders in the competition for binding to the RIP. This work offers promising leads for developing RIP therapeutics. The results suggest that the doorstop approach might also be applicable in the development of other protein•polynucleotide inhibitors as antiviral agents such as inhibitors of the Z-DNA binding proteins in poxviruses. This work also calls for careful chemical and biological characterization of drug leads obtained from chemical screens to avoid the identification of irrelevant chemical structures and to avoid the interference caused by direct interactions between the chemicals being screened and the luciferase reporter used in screening assays.

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.

The rate of spontaneous cleavage of the glycosidic bond of adenosine

Previous estimates of the rate of spontaneous cleavage of the glycosidic bond of adenosine were determined by extrapolating the rates of the acid- and base-catalyzed reactions to neutral pH. Here we show that cleavage also proceeds through a pH-independent mechanism. Rate constants were determined as a function of temperature at pH 7 and a linear Arrhenius plot was constructed. Uncatalyzed cleavage occurs with a rate constant of 3.7x10(-12)s(-1) at 25 degrees C, and the rate enhancement generated by the corresponding glycoside hydrolase is approximately 5x10(12)-fold.

Identification of new classes of ricin toxin inhibitors by virtual screening

We used two virtual screening programs, ICM and GOLD, to dock nearly 50,000 compounds into each of two conformations of the target protein ricin A chain (RTA). A limited control set suggests that candidates scored highly by two programs may have a higher probability of being ligands than those in a list from a single program. Based on the virtual screens, we purchased 306 compounds that were subjected to a kinetic assay. Six compounds were found to give modest, but significant, inhibition of RTA. They also tended to inhibit Shiga toxin A chain, with roughly the same IC(50). The compounds generally represent novel chemical platforms that do not resemble RTA substrates, as currently known inhibitors do. These six were also tested in a cell-based assay for their ability to protect cells from intact ricin. Two compounds were effective in this regard, showing modest to strong ricin inhibition, but also showing some cytotoxicity. RTA, with its large, polar active site is a difficult drug design target which is expected to bind small molecules only weakly. The ability of the method to find these novel platforms is encouraging and suggests virtual screening can contribute to the search for ricin and Shiga toxin inhibitors.

Transition state analogues in structures of ricin and saporin ribosome-inactivating proteins

Ricin A-chain (RTA) and saporin-L1 (SAP) catalyze adenosine depurination of 28S rRNA to inhibit protein synthesis and cause cell death. We present the crystal structures of RTA and SAP in complex with transition state analogue inhibitors. These tight-binding inhibitors mimic the sarcin-ricin recognition loop of 28S rRNA and the dissociative ribocation transition state established for RTA catalysis. RTA and SAP share unique purine-binding geometry with quadruple pi-stacking interactions between adjacent adenine and guanine bases and 2 conserved tyrosines. An arginine at one end of the pi-stack provides cationic polarization and enhanced leaving group ability to the susceptible adenine. Common features of these ribosome-inactivating proteins include adenine leaving group activation, a remarkable lack of ribocation stabilization, and conserved glutamates as general bases for activation of the H(2)O nucleophile. Catalytic forces originate primarily from leaving group activation evident in both RTA and SAP in complex with transition state analogues.

Transition state analogues rescue ribosomes from saporin-L1 ribosome inactivating protein

Ribosome inactivating proteins (RIPs) catalyze the hydrolytic depurination of one or more adenosine residues from eukaryotic ribosomes. Depurination of the ribosomal sarcin-ricin tetraloop (GAGA) causes inhibition of protein synthesis and cellular death. We characterized the catalytic properties of saporin-L1 from Saponaria officinalis (soapwort) leaves, and it demonstrated robust activity against defined nucleic acid substrates and mammalian ribosomes. Transition state analogue mimics of small oligonucleotide substrates of saporin-L1 are powerful, slow-onset inhibitors when adenosine is replaced with the transition state mimic 9-deazaadenine-9-methylene-N-hydroxypyrrolidine (DADMeA). Linear, cyclic, and stem-loop oligonucleotide inhibitors containing DADMeA and based on the GAGA sarcin-ricin tetraloop gave slow-onset tight-binding inhibition constants (K(i)*) of 2.3-8.7 nM under physiological conditions and bind up to 40000-fold tighter than RNA substrates. Saporin-L1 inhibition of rabbit reticulocyte translation was protected by these inhibitors. Transition state analogues of saporin-L1 have potential in cancer therapy that employs saporin-L1-linked immunotoxins.

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.

Detecting ricin: sensitive luminescent assay for ricin A-chain ribosome depurination kinetics

Ricin is a family member of the lethal ribosome-inactivating proteins (RIP) found in plants. Ricin toxin A-chain (RTA) from castor beans catalyzes the hydrolytic depurination of a single base from a GAGA tetraloop of eukaryotic rRNA to release a single adenine from the sarcin-ricin loop (SRL). Protein synthesis is inhibited by loss of the elongation factor binding site resulting in cell death. We report a sensitive coupled assay for the measurement of adenine released from ribosomes or small stem-loop RNAs by RTA catalysis. Adenine phosphoribosyl transferase (APRTase) and pyruvate orthophosphate dikinase (PPDK) convert adenine to ATP for quantitation by firefly luciferase. The resulting AMP is cycled to ATP to give sustained luminescence proportional to adenine concentration. Subpicomole adenine quantitation permits the action of RTA on eukaryotic ribosomes to be followed in continuous, high-throughput assays. Facile analysis of RIP catalytic activity will have applications in plant toxin detection, inhibitor screens, mechanistic analysis of depurinating agents on oligonucleotides and intact ribosomes, and in cancer immunochemotherapy. Kinetic analysis of the catalytic action of RTA on rabbit reticulocyte 80S ribosomes establishes a catalytic efficiency of 2.6 x 10(8) M(-1) s(-1), a diffusion limited reaction indicating catalytic perfection even with large reactants.

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.

Concentration-dependent processivity of multiple glutamate ligations catalyzed by folylpoly-gamma-glutamate synthetase

Folylpoly-gamma-glutamate synthetase (FPGS, EC 6.3.2.17) is an ATP-dependent ligase that catalyzes formation of poly-gamma-glutamate derivatives of reduced folates and antifolates such as methotrexate and 5,10-dideaza-5,6,7,8-tetrahydrofolate (DDAH 4PteGlu 1). While the chemical mechanism of the reaction catalyzed by FPGS is known, it is unknown whether single or multiple glutamate residues are added following each folate binding event. A very sensitive high-performance liquid chromatography method has been used to analyze the multiple ligation reactions onto radiolabeled DDAH 4PteGlu 1 catalyzed by FPGS to distinguish between distributive or processive mechanisms of catalysis. Reaction time courses, substrate trapping, and pulse-chase experiments were used to assess folate release during multiple glutamate additions. Together, the results of these experiments indicate that hFPGS can catalyze the processive addition of approximately four glutamate residues to DDAH 4PteGlu 1. The degree of processivity was determined to be dependent on the concentration of the folate substrate, thus suggesting a mechanism for the regulation of folate polyglutamate synthesis in cells.

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.

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.

A new interpretation of chlorine leaving group kinetic isotope effects; a theoretical approach

The chlorine leaving group kinetic isotope effects (KIEs) for the S(N)2 reactions between methyl chloride and a wide range of anionic, neutral, and radical anion nucleophiles were calculated in the gas phase and, in several cases, using a continuum solvent model. In contrast to the expected linear dependence of the chlorine KIEs on the C(alpha)-Cl bond order in the transition state, the KIEs fell in a very small range (1.0056-1.0091), even though the C(alpha)-Cl transition state bond orders varied widely from approximately 0.32 to 0.78, a range from reactant-like to very product-like. This renders chlorine KIEs, and possibly other leaving-group KIEs, less useful for studies of reaction mechanisms than commonly assumed. A partial explanation for this unexpected relationship between the C(alpha)-Cl transition state bond order and the magnitude of the chlorine KIE is presented.

Leaving group activation by aromatic stacking: an alternative to general acid catalysis

General acid catalysis is a powerful and widely used strategy in enzymatic nucleophilic displacement reactions. For example, hydrolysis/phosphorolysis of the N-glycosidic bond in nucleosides and nucleotides commonly involves the protonation of the leaving nucleobase concomitant with nucleophilic attack. However, in the nucleoside hydrolase of the parasite Trypanosoma vivax, crystallographic and mutagenesis studies failed to identify a general acid. This enzyme binds the purine base of the substrate between the aromatic side-chains of Trp83 and Trp260. Here, we show via quantum chemical calculations that face-to-face stacking can raise the pKa of a heterocyclic aromatic compound by several units. Site-directed mutagenesis combined with substrate engineering demonstrates that Trp260 catalyzes the cleavage of the glycosidic bond by promoting the protonation of the purine base at N-7, hence functioning as an alternative to general acid catalysis.

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.

Catalysis by nucleoside hydrolases

Nucleoside hydrolases cleave the N-glycosidic bond of ribonucleosides. Because of their vital role in the protozoan purine salvage pathway, nucleoside hydrolases from parasitic protozoa in particular have been studied extensively by X-ray crystallography, kinetic methods and site-directed mutagenesis. An elaborate network of conserved interactions between the metalloenzyme and the ribose enables steric and electrostatic stabilisation of the oxocarbenium-ion-like transition state. Activation of the leaving group by protonation before the formation of the transition state is a recurring catalytic strategy of enzymes that cleave N-glycosidic bonds. However, the mechanisms underlying leaving group activation are still the subject of debate for the nucleoside hydrolases.

Combinatorial minimization and secondary structure determination of a nucleotide synthase ribozyme

We previously isolated from random sequences ribozymes able to form a glycosidic linkage between a ribose sugar and 4-thiouracil in a reaction that mimics protein-catalyzed nucleotide synthesis. Here we report on two serial in vitro selection experiments that defined the core motif of one of the nucleotide synthase ribozymes and provided improved versions of this ribozyme. The first selection experiment started from a degenerate sequence pool based on the previously isolated sequence and used a selection-amplification protocol that allowed the sequence requirements at the 3' terminus of the ribozyme to be interrogated. Comparing the active sequences identified in this experiment revealed the complicated secondary structure of the nucleotide synthase ribozyme. A second selection was then performed to remove nonessential sequence from the ribozyme. This selection started with a pool with variation introduced in both the sequence and the length of the nonconserved loops and joining regions. This pool was generated using a partial reblocking/deblocking strategy on a DNA synthesizer, allowing the combinatorial synthesis of both point deletions and point substitutions. The consensus ribozyme motif that emerged was an approximately 71 nt pseudoknot structure with five stems and two important joining segments. Comparative sequence analysis and a cross-linking experiment point to the probable location of nucleotide synthesis. The prototype isolate from the second selection was nearly 35 times more efficient than the initial isolate and at least 10(8) times more efficient than an upper limit of an as-yet undetectable uncatalyzed reaction, supporting the idea that RNA-catalyzed nucleotide synthesis might have been important in an RNA world.

Glycosidase mechanisms

The three-dimensional structure of glycosidases and of their complexes and the study of transition-state mimics reveal structural details that correlate with mechanism. Of particular interest are the transition-state conformations harnessed by individual enzymes and the substrate distortion observed in enzyme-ligand complexes. 3D-structure in synergy with transition-state mimicry opens the way for mechanistic interpretation of enzyme inhibition and for the development of therapeutic agents.

Transition state variation in enzymatic reactions

Experimental analysis of enzymatic transition states by kinetic isotope effect methods has established geometric variation in related transition state structures. Differences are apparent in development of the reaction coordinate, in solvolytic transition states relative to those in enzymatic catalytic sites, in the stereochemistry of related substrates at the transition state, and in reactions catalyzed by related enzymes.