Theoretical examination of Mg(2+)-mediated hydrolysis of a phosphodiester linkage as proposed for the hammerhead ribozyme

The hammerhead ribozyme is an RNA molecule capable of self-cleavage at a unique site within its sequence. Hydrolysis of this phosphodiester linkage has been proposed to occur via an in-line attack geometry for nucleophilic displacement by the 2'-hydroxyl on the adjoining phosphorus to generate a 2',3'-cyclic phosphate ester with elimination of the 5'-hydroxyl group, requiring a divalent metal ion under physiological conditions. The proposed S(N)2(P) reaction mechanism was investigated using density functional theory calculations incorporating the hybrid functional B3LYP to study this metal ion-dependent reaction with a tetraaquo magnesium (II)-bound hydroxide ion. For the Mg(2+)-catalyzed reaction, the gas-phase geometry optimized calculations predict two transition states with a kinetically insignificant, yet clearly defined, pentacoordinate intermediate. The first transition state located for the reaction is characterized by internal nucleophilic attack coupled to proton transfer. The second transition state, the rate-determining step, involves breaking of the exocyclic P-O bond where a metal-ligated water molecule assists in the departure of the leaving group. These calculations demonstrate that the reaction mechanism incorporating a single metal ion, serving as a Lewis acid, functions as a general base and can afford the necessary stabilization to the leaving group by orienting a water molecule for catalysis.     

Effect of metal ion exchange on the photocurrent response from bacteriorhodopsin on tin oxide electrodes

The transient photocurrent response from bacteriorhodopsin (bR) on tin oxide electrodes was strongly influenced by metal ions bound to bR molecules. The photocurrent polarity reversal pH, which corresponded to the pH value for the reversal of the proton release/uptake sequence in the bR photocycle, of cation-substituted purple membrane (PM) was shifted to lower pH with the increase in the cation affinities to carboxyl groups and a close correlation was noted between the two values. This suggests that the metal ion present in the extracellular region of a bR molecule modulates the pK(a) of proton release groups of bR by stabilizing the ionized state of the proton-releasing glutamic acids. The behavior of photocurrents at light-off in alkaline media, reflecting the proton uptake by bR, was unchanged by binding monovalent (Na(+) and K(+)) or divalent cations (Mg(2+) and Ca(2+)), but was drastically changed by binding La(3+) ions. This can be explained by invoking a substantial slowing of the proton uptake process in the presence of La(3+).     

Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Br?nsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Br?nsted base.     

Direct pK(a) measurement of the active-site cytosine in a genomic hepatitis delta virus ribozyme

Hepatitis delta virus ribozymes have been proposed to perform self-cleavage via a general acid/base mechanism involving an active-site cytosine, based on evidence from both a crystal structure of the cleavage product and kinetic measurements. To determine whether this cytosine (C75) in the genomic ribozyme has an altered pK(a) consistent with its role as a general acid or base, we used (13)C NMR to determine its microscopic pK(a) in the product form of the ribozyme. The measured pK(a) is moderately shifted from that of a free nucleoside or a base-paired cytosine and has the same divalent metal ion dependence as the apparent reaction pK(a)'s measured kinetically. However, under all conditions tested, the microscopic pK(a) is lower than the apparent reaction pK(a), supporting a model in which C75 is deprotonated in the product form of the ribozyme at physiological pH. While additional results suggest that the pK(a) is not shifted in the reactant state of the ribozyme, these data cannot rule out elevation of the C75 pK(a) in an intermediate state of the transesterification reaction.     

How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

We have carried out density functional theory QM/MM calculations on the catalytic subunit of cAMP-dependent protein kinase (PKA). The QM/MM calculations indicate that the phosphorylation reaction catalyzed by PKA is mainly dissociative, and Asp166 serves as the catalytic base to accept the proton delivered by the substrate peptide. Among the key interactions in the active site, the Mg(2+) ions, glycine rich loop, and Lys72 are found to stabilize the transition state through electrostatic interactions. On the other hand, Lys168, Asn171, Asp184, and the conserved waters bound to Mg(2+) ions do not directly contribute to lower the energy barrier of the phosphorylation reaction, and possible roles for these residues are proposed. The QM/MM calculations with different QM/MM partition schemes or different initial structures yield consistent results. In addition, we have carried out 12 ns molecular dynamics simulations on both wild type and K168A mutated PKA, respectively, to demonstrate that the catalytic role of Lys168 is to keep ATP and substrate peptide in the near-attack reactive conformation.     

The active site of a zinc-dependent metalloproteinase influences the computed pK(a) of ligands coordinated to the catalytic zinc ion

TNF-alpha converting enzyme (TACE) is a multidomain, membrane-anchored protein that includes a Zn-dependent protease domain. It releases the soluble form of cytokine tumor necrosis factor-alpha (TNF-alpha) from its membrane-bound precursor. TACE is a metalloprotease containing a catalytic glutamic acid, Glu-406, and a Zn(2+) ion ligated to three imidazoles. The protonation states of the active site glutamic acid and inhibitors are important factors in understanding the potency of inhibitors with acidic zinc-ligating groups such as hydroxamic and carboxylic acids. Density functional methods were utilized to compute pK(a) values using a model of the catalytic site of TACE and to predict a concomitant mechanism of binding, consistent with lowering the pK(a) of the bound ligand and raising the pK(a) of the active site Glu-406. Weak acids, such as hydroxamic acids, bind in their neutral form and then transfer an acidic proton to Glu-406. Stronger acids, such as carboxylic acids, bind in their anionic form and require preprotonation of Glu-406. Similar binding events would be expected for other zinc-dependent proteases.     

Competition between Li+ and Mg2+ in metalloproteins. Implications for lithium therapy

Lithium is used (in the form of soluble salts) to treat bipolar disorder and has been considered as a possible drug in treating chronic neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. One of the proposed mechanisms of Li(+) action involves a competition between the alien Li(+) and native Mg(2+) for metal-binding sites and subsequent inhibition of key enzymes involved in specific neurotransmission pathways, but not vital Mg(2+) proteins in the cell. This raises the following intriguing questions: Why does Li(+) replace Mg(2+) only in enzymes involved in bipolar disorder, but not in Mg(2+) proteins essential to cells? In general, what factors allow monovalent Li(+) to displace divalent Mg(2+) in proteins? Specifically, how do the composition, overall charge, and solvent exposure of the metal-binding site as well as a metal-bound phosphate affect the selectivity of Li(+) over Mg(2+)? Among the many possible factors, we show that the competition between Mg(2+) and Li(+) depends on the net charge of the metal complex, which is determined by the numbers of metal cations and negatively charged ligands, as well as the relative solvent exposure of the metal cavity. The protein itself is found to select Mg(2+) over the monovalent Li(+) by providing a solvent-inaccessible Mg(2+)-binding site lined by negatively charged Asp/Glu, whereas the cell machinery was found to select Mg(2+) among other competing divalent cations in the cellular fluids such as Ca(2+) and Zn(2+) by maintaining a high concentration ratio of Mg(2+) to its biogenic competitor in various biological compartments. The calculations reveal why Li(+) replaces Mg(2+) only in enzymes that are known targets of Li(+) therapy, but not in Mg(2+) enzymes essential to cells, and also reveal features common to the former that differ from those in the latter proteins.     

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.     

Density functional theory investigation on the mechanism of the hepatitis delta virus ribozyme

Density functional theory methods have been used to investigate the hepatitis delta virus (HDV) ribozyme and its catalyzed phosphodiester cleavage. In particular, the effects of the environment's polarity and/or specific hydrogen-bond interactions on the proton affinity of the active site cytosine's N3 ring center have been considered. In addition, the basicities of possible hydrated Mg2+ ion species were also examined. The mechanism previously proposed for the HDV ribozyme in which the active site cytosine (C75) is protonated and thus acts as an acid while the Mg2+ species acts as the complementary base was then investigated. The possible role of tautomerization of C75 is also discussed.     

Theoretical study toward understanding the catalytic mechanism of pyruvate decarboxylase

Density functional calculations are employed to explore the mechanisms of all elementary reaction steps involved in the catalytic cycle of pyruvate decarboxylase (PDC). Different models are constructed for mimicking the involvement of some key residues in a certain step. The effect of the protein framework on the potential energy profiles of active site models is approximately modeled by fixing some freedoms, based on the crystal structure of the PDC enzyme from Saccharomyces cerevisiae (ScPDC). Our calculations confirm that Glu51 is the most important residue in the formation of the ylide and the release of acetaldehyde via the proton relay between Glu51, N1', and the 4'-amino group of thiamine diphosphate. The presence of Glu477 and Asp28 residues makes the decarboxylation of lactylthiamin diphosphate (LThDP) an endothermic process with a significant free energy barrier. The protonation of the alpha-carbanion to form 2-(1-hydroxyethyl)-thiamin diphosphate is found to go through a concerted double proton transfer transition state involving both Asp28 and His115 residues. The final step, acetaldehyde release, is likely to proceed through a concerted transition state involving carbon-carbon bond-breaking and the deprotonation of the alpha-hydroxyl group. The decarboxylation of LThDP and the protonation of the alpha-carbanion are two rate-limiting steps, relative to the facile occurrence of the ylide formation and acetaldehyde release. The catalytic roles of residues Glu51, Glu477, Asp28, and Gly417 in the active site of ScPDC in individual steps elucidated from the present study are in good agreement with those derived from site-directed mutagenesis.     

Interactions of the hormone oxytocin with divalent metal ions

The interaction of the cyclic nonapeptide oxytocin (OT) with a number of alkaline earth and divalent transition metal ions (X(2+)) was examined employing mass spectrometry (MS) and ion mobility spectrometry (IMS) techniques in combination with molecular dynamics (MD) and density functional theory (DFT) calculations. Under acidic conditions it was found that OT exhibits an exceptionally strong affinity for all divalent metal ions resulting in strong [OT + X](2+) peaks in the mass spectrum. Under basic conditions only Cu(2+) and Ni(2+)-OT complexes were detected and these were singly, doubly, triply, or quadruply deprotonated. Collision-induced dissociation of the [OT - 3H + Cu](-) complex yielded exclusively C-terminal Cu(2+)-containing fragments (Cu(2+)fragment(3-)), suggesting that the Cu(2+) ligation site includes deprotonated C-terminal backbone amide nitrogen atoms and the N-terminal amino nitrogen atom in [OT - 3H + Cu](-). MD and DFT calculations indicate a square-planar complex is consistent with these observations and with experimental collision cross sections. MD and DFT calculations also indicate either an octahedral or trigonal-bipyramidal complex between Zn(2+) and OT is lowest in energy with carbonyl oxygens being the primary ligation sites. Both complexes yield cross sections in agreement with experiment. The biological impact of the structural changes induced in OT by divalent metal ion coodination is discussed.     

Elucidating the mechanism for the reduction of nitrite by copper nitrite reductase--a contribution from quantum chemical studies

Density functional methods have been applied to investigate the properties of the active site of copper-containing nitrite reductases and possible reaction mechanisms for the enzyme catalysis. The results for a model of the active site indicate that a hydroxyl intermediate is not formed during the catalytic cycle, but rather a state with a protonated nitrite bound to the reduced copper. Electron affinity calculations indicate that reduction of the T2 copper site does not occur immediately after nitrite binding. Proton affinity calculations are indicative of substantial pK(a) differences between different states of the T2 site. The calculations further suggest that the reaction does not proceed until uptake of a second proton from the bulk solution. They also indicate that Asp-92 may play both a key role as a proton donor to the substrate, and a structural role in promoting catalysis. In the D92N mutant another base, presumably a nearby histidine (His-249) may take the role as the proton donor. On the basis of these model calculations and available experimental evidence, an ordered reaction mechanism for the reduction of nitrite is suggested. An investigation of the binding modes of the nitric oxide product and the nitrite substrate to the model site has also been made, indicating that nitric oxide prefers to bind in an end-on fashion to the reduced T2 site.     

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.     

Intramolecular proton transfer induced by divalent alkali earth metal cation in the gas state

Interactions between divalent alkali earth metal (DAEM) ions M (M?Be, Mg, Ca, Sr, Ba) and the second stable glycine conformer in the gas phase, which can transfer into the ground‐state glycine‐M²⁺ (except the glycine–Be²⁺) among each corresponding isomers when these divalent metal ions are bound, are studied at the hybrid three‐parameter B3LYP level with three different basis sets. Proton transfers from the hydroxyl to the amino nitrogen of the glycine without energy barriers have been first observed in the gas phase in these glycine–M²⁺ systems. The interaction between the glycine and these DAEM ions except beryllium and magnesium ion only create an amino hydrogen pointing to the original hydroxyl due to their weaker interaction relative to those divalent transition metal (DTM) ion‐bound glycine derivatives, being obviously different from that between the glycine and DTM ions, in which two amino hydrogens point to the original hydroxyl oxygen when these metal‐chelated glycine derivatives are produced. The interaction energy between the glycine and divalent magnesium would be the boundary of one or two amino hydrogens pointing to the hydrogyl oxygen, i.e., the ?170.3 kcal/mol of binding energy is a critical point. Similar intramolecular proton transfer has also been predicted for those DTM ion‐chelated glycine systems; however, that in the gas state has not been observed in the monovalent metal ion‐coordinated glycine systems. The binding energy between some monovalent TM ion and the glycine is similar to that of the glycine–Ba²⁺, which has the lowest binding strength among these DAEM–ion chelated glycine complexes. The difference among them only lies in the larger electrostatic and polarized effects in the latter, which favor the stability of the zwitterionic glycine form in the gas phase. According to these observations, we predict that the zwitterionic glycine would exist in the field of two positive charges in the gas phase. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 205–214, 2003     

Quantum chemical studies of a model for peptide bond formation. 3. Role of magnesium cation in formation of amide and water from ammonia and glycine

The SN2 reaction between glycine and ammonia molecules with magnesium cation Mg2+ as a catalyst has been studied as a model reaction for Mg(2+)-catalyzed peptide bond formation using the ab initio Hartree-Fock molecular orbital method. As in previous studies of the uncatalyzed and amine-catalyzed reactions between glycine and ammonia, two reaction mechanisms have been examined, i.e., a two-step and a concerted reaction. The stationary points of each reaction including intermediate and transition states have been identified and free energies calculated for all geometry-optimized reaction species to determine the thermodynamics and kinetics of each reaction. Substantial decreases in free energies of activation were found for both reaction mechanisms in the Mg(2+)-catalyzed amide bond formation compared with those in the uncatalyzed and amine-catalyzed amide bond formation. The catalytic effect of the Mg2+ cation is to stabilize both the transition states and intermediate, and it is attributed to the neutralization of the developing negative charge on the electrophile and formation of a conformationally flexible nonplanar five-membered chelate ring structure.     

Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.     

Effects of cations on the hydrogen bond network of liquid water: new results from X-ray absorption spectroscopy of liquid microjets

Oxygen K-edge X-ray absorption spectra (XAS) of aqueous chloride solutions have been measured for Li(+), Na(+), K(+), NH(4)(+), C(NH(2))(3)(+), Mg(2+), and Ca(2+) at 2 and 4 M cation concentrations. Marked changes in the liquid water XAS are observed upon addition of the various monovalent cation chlorides that are nearly independent of the identity of the cation. This indicates that interactions with the dissolved monovalent cations do not significantly perturb the unoccupied molecular orbitals of water molecules in the vicinity of the cations and that water-chloride interactions are primarily responsible for the observed spectral changes. In contrast, the addition of the divalent cations engenders changes unique from the case of the monovalent cations, as well as from each other. Density functional theory calculations suggest that the ion-specific spectral variations arise primarily from direct electronic perturbation of the unoccupied orbitals due to the presence of the ions, probably as a result of differences in charge transfer from the water molecules onto the divalent cations.     

Detection of a proton-transfer process by kinetic solvent isotope effects in NH(4)(+)-mediated reactions catalyzed by a hammerhead ribozyme

Hammerhead ribozymes have been considered to be metalloenzymes. However, this proposal was recently questioned by the finding that the reaction proceeds in the presence of high concentrations of monovalent ions such as NH(4)(+) ions and in the absence of any divalent metal ions. Our present analysis based on solvent isotope effects indicates that (1) a proton transfer(s) occurs only in the NH(4)(+)-mediated reaction but not in metal-ion-mediated reactions such as Mg(2+)- and Li(+)-mediated reactions, (2) the catalyst that stabilizes the 5' leaving group in the NH(4)(+)-mediated reaction is different from that in the metal-ion-mediated HH ribozyme reactions, (3) an NH(4)(+) ion seems to act as a general acid catalyst, and (4) a nucleobase alone should not be the catalyst.