An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.     

The Ambiguous Base-Pairing and High Substrate Efficiency of T-705 (Favipiravir) Ribofuranosyl 5′-Triphosphate towards Influenza A Virus Polymerase

T-705 (Favipiravir) is a broad-spectrum antiviral molecule currently in late stage clinical development for the treatment of influenza virus infection. Although it is believed that T-705 potency is mediated by its ribofuranosyl triphosphate (T-705 RTP) metabolite that could be mutagenic, the exact molecular interaction with the polymerase of influenza A virus (IAVpol) has not been elucidated. Here, we developed a biochemical assay to measure the kinetics of nucleotide incorporation by IAVpol in the elongation mode. In this assay, T-705 RTP was recognized by IAVpol as an efficient substrate for incorporation to the RNA both as a guanosine and an adenosine analog. Compared to natural GTP and ATP, the discrimination of T-705 RTP was about 19- and 30-fold, respectively. Although the single incorporation of the ribonucleotide monophosphate form of T-705 did not efficiently block RNA synthesis, two consecutive incorporation events prevented further primer extension. In comparison, 3′-deoxy GTP caused immediate chain termination but was incorporated less efficiently by the enzyme, with a discrimination of 4,900-fold relative to natural GTP. Collectively, these results provide the first detailed biochemical characterization to evaluate the substrate efficiency and the inhibition potency of nucleotide analogs against influenza virus polymerase. The combination of ambiguous base-pairing with low discrimination of T-705 RTP provides a mechanistic basis for the in vitro mutagenic effect of T-705 towards influenza virus.     

4-thio-U cross-linking identifies the active site of the VS ribozyme

To identify nucleotides in or near the active site, we have used a circularly permuted version of the VS ribozyme capable of cleavage and ligation to incorporate a single photoactive nucleotide analog, 4-thio- uridine, immediately downstream of the scissile bond. Exposure to UV light produced two cross-linked RNAs, in which the 4-thio-uridine was cross-linked to A756 in the 730 loop of helix VI. The cross-links formed only under conditions that support catalytic activity, suggesting that they reflect functionally relevant conformations of the RNA. One of the cross-linked RNAs contains a lariat, indicative of intramolecular cross-linking in the ligated RNA; the other is a branched molecule in which the scissile phosphodiester bond is cleaved, but occupies the same site in the ribozyme-substrate complex. These are the two forms of the RNA expected to be the ground state structures on either side of the transition state. This localization of the active site is consistent with previous mutational, biochemical and biophysical data, and provides direct evidence that the cleavage site in helix I interacts with the 730 loop in helix VI.     

Investigation of adenosine base ionization in the hairpin ribozyme by nucleotide analog interference mapping.

Tertiary structure in globular RNA folds can create local environments that lead to pKa perturbation of specific nucleotide functional groups. To assess the prevalence of functionally relevant adenosine-specific pKa perturbation in RNA structure, we have altered the nucleotide analog interference mapping (NAIM) approach to include a series of a phosphorothioate-tagged adenosine analogs with shifted N1 pKa values. We have used these analogs to analyze the hairpin ribozyme, a small self-cleaving/ligating RNA catalyst that is proposed to employ a general acid-base reaction mechanism. A single adenosine (A10) within the ribozyme active site displayed an interference pattern consistent with a functionally significant base ionization. The exocyclic amino group of a second adenosine (A38) contributes substantially to hairpin catalysis, but ionization of the nucleotide does not appear to be important for activity. Within the hairpin ribozyme crystal structure, A10 and A38 line opposite edges of a solvent-excluded cavity adjacent to the 5'-OH nucleophile. The results are inconsistent with the model of ribozyme chemistry in which A38 acts as a general acid-base catalyst, and suggest that the hairpin ribozyme uses an alternative mechanism to achieve catalytic rate enhancement that utilizes functional groups within a solvent-excluded cleft in the ribozyme active site.     

Nucleotide analog interference mapping of the hairpin ribozyme: implications for secondary and tertiary structure formation.

The hairpin ribozyme is a small, naturally occurring RNA capable of folding into a distinct three-dimensional structure and catalyzing a specific phosphodiester transfer reaction. We have adapted a high throughput screening procedure entitled nucleotide analog interference mapping (NAIM) to identify functional groups important for proper folding and catalysis of this ribozyme. A total of 18 phosphorothioate-tagged nucleotide analogs were used to determine the contribution made by individual ribose 2'-OH and purine functional groups to the hairpin ribozyme ligation reaction. Substitution with 2'-deoxy-nucleotide analogs disrupted activity at six sites within the ribozyme, and a unique interference pattern was observed at each of the 11 conserved purine nucleotides. In most cases where such information is available, the NAIM data agree with the previously reported single-site substitution results. The interference patterns are interpreted in comparison to the isolated loop A and loop B NMR structures and a model of the intact ribozyme. These data provide biochemical evidence in support of many, but not all, of the non-canonical base-pairs observed by NMR in each loop, and identify the functional groups most likely to participate in the tertiary interface between loop A and loop B. These groups include the 2'-OH groups of A10, G11, U12, C25, and A38, the exocyclic amine of G11, and the minor groove edge of A9 and A24. The data also predict non-A form sugar pucker geometry at U39 and U41. Based upon these results, a revised model for the loop A tertiary interaction with loop B is proposed. This work defines the chemical basis of purine nucleotide conservation in the hairpin ribozyme, and provides a basis for the design and interpretation of interference suppression experiments.     

Autolytic processing of a phosphorothioate diester bond.

A small satellite RNA of tobacco ringspot virus replicates in tissues infected with tobacco ringspot virus and accumulates in virus capsids, forming virus-like particles. Previous research showed that multimeric forms of this satellite RNA have tandem repeats of the "monomeric" satellite RNA sequence of 359 or 360 nucleotide residues. The multimeric RNAs undergo autolytic processing at a specific CpA phosphodiester bond, the junction, to generate the monomeric RNA. We substituted phosphorothioate diester bonds for various sets of phosphodiester bonds, in dimeric and truncated forms of the satellite RNA. The degree of reduction in autolytic cleavage varied both with the sites of substitution and the size of the RNA molecules. Analyses of a product of the autolysis reaction suggest that one phosphorothioate diester bond most strongly interferes with processing, the one introduced at the CpA junction during its synthesis from adenosine-5'-0-(1-thiotriphosphate). However, extensive introduction of phosphorothioate diester bonds elsewhere in the molecule also decreased processing, possibly by altering conformation.     

Biochemical detection of monovalent metal ion binding sites within RNA

Many RNAs, including the ribosome, RNase P, and the group II intron, explicitly require monovalent cations for activity in vitro. Although the necessity of monovalent cations for RNA function has been known for more than a quarter of a century, the characterization of specific monovalent metal sites within large RNAs has been elusive. Here we describe a biochemical approach to identify functionally important monovalent cations in nucleic acids. This method uses thallium (Tl+), a soft Lewis acid heavy metal cation with chemical properties similar to those of the physiological alkaline earth metal potassium (K+). Nucleotide analog interference mapping (NAIM) with the sulfur-substituted nucleotide 6-thioguanosine in combination with selective metal rescue of the interference with Tl+ provides a distinct biochemical signature for monovalent metal ion binding. This approach has identified a K+ binding site within the P4-P6 domain of the Tetrahymena group I intron that is also present within the X-ray crystal structure. The technique also predicted a similar binding site within the Azoarcus group I intron where the structure is not known. The approach is applicable to any RNA molecule that can be transcribed in vitro and whose function can be assayed.     

Crystal structure of an RNA duplex containing phenyl-ribonucleotides, hydrophobic isosteres of the natural pyrimidines

Chemically modified nucleotide analogs have gained widespread popularity for probing structure-function relationships. Among the modifications that were incorporated into RNAs for assessing the role of individual functional groups, the phenyl nucleotide has displayed surprising effects both in the contexts of the hammerhead ribozyme and pre-mRNA splicing. To examine the conformational properties of this hydrophobic base analog, we determined the crystal structure of an RNA double helix with incorporated phenyl ribonucleotides at 1.97 A resolution. In the structure, phenyl residues are engaged in self-pairing and their arrangements suggest energetically favorable stacking interactions with 3'-adjacent guanines. The presence of the phenyl rings in the center of the duplex results in only moderate changes of the helical geometry. This finding is in line with those of earlier experiments that showed the phenyl analog to be a remarkably good mimetic of natural base function. Because the stacking interactions displayed by phenyl residues appear to be similar to those for natural bases, reduced conformational restriction due to the lack of hydrogen bonds with phenyl as well as alterations in its solvent structure may be the main causes of the activity changes with phenyl-modified RNAs.     

In vitro selection of an archaeal RNase P RNA mimics natural variation

Archaeal and bacterial RNase P RNAs are similar in sequence and secondary structure, but in the absence of protein, the archaeal RNAs are much less active and require extreme ionic conditions for activity. To assess how readily the activity of the archaeal RNA alone could be improved by small changes in sequence, in vitro selection was used to generate variants of a Methanobacterium formicicum RNase P RNA: Bacillus subtilus pre-tRNA(Asp) self-cleaving conjugate RNA. Functional variants were generated with a spectrum of mutations that were predominately consistent with natural variation in this RNA. Variants generated from the selection had cleavage rates comparable to that of wild type; variants with improved cleavage rates or lower ionic requirements were not obtained. This suggests that the RNase P RNAs of Bacteria and Archaea are globally optimized and the basis for the large biochemical differences between these two types of RNase P RNA is distributed in the molecule.     

Quantitation of free energy profiles in RNA-ligand interactions by nucleotide analog interference mapping

RNA interactions with protein and small molecule ligands serve a wide variety of biochemical functions in the cell. To best understand the specificity and affinity of these interactions, the free energy contribution made by individual function groups in the RNA must be determined. As an efficient method for obtaining such energetic profiles, we report quantitative nucleotide analog interference mapping (QNAIM). This extension of the NAIM methodology uses the magnitude of analog interference as a function of ligand concentration to calculate binding constants for RNA with individual analog substitutions. In this way, QNAIM not only defines which functional groups are important to an interaction but simultaneously determines the energetic contribution made by each occurrence of that functional group within the RNA polymer. To establish the utility of this approach, QNAIM was used to quantify functional group interactions within the signal recognition particle (SRP), specifically the 4.5S RNA with the M domain of Ffh. In each of the cases in which energetic data were available from previous site-specific substitution analyses, QNAIM provided nearly equivalent results. These experiments on a model system demonstrate that QNAIM is an efficient method to establish a chemically detailed free energy profile for a wide variety of RNA-ligand interactions.     

Co-optimization of ribozyme substrate stacking and L-arginine binding.

A model of the Tetrahymena catalytic site predicts that nucleotide 262 (nt262) caps an RNA pocket in which nucleoside substrates and arginine-like competitive inhibitors reside. Here we show that substituted RNAs behave as if nt262 stacks on nucleoside substrates, supporting the model. The more frequent an nt262 is in natural sequences, the more reactive the corresponding Tetrahymena RNA is for both cognate and non-cognate nucleoside substrates. These more reactive RNAs with the majority nt262 also bind arginine more strongly, stereoselect more strongly in favor of L-arginine, and make a greater distinction between the somewhat similar side-chains of L-arginine and L-lysine. These parallel changes in interaction with nucleosides and arginine analogs seem best explained by stacking of the arginine's guanidino group under the nt262 base. One consequence is that selection for improved Tetrahymena catalysis with nucleosides should also yield an improved arginine site.     

Synthesis and properties of fluorescent nucleotide substrates for DNA-dependent RNA polymerases.

A new class of fluorescent nucleotide analogs which contain the fluorophore 1-aminonaphthalene-5-sulfonate attached via a gamma-phosphoamidate bond has been synthesized. Both the purine and pyrimidine analogs have fluorescence emission maxima at 460 nm. Cleavage of the alpha-beta-phosphoryl bond produces change in both the absorption and fluorescence emission spectra. The fluorescence of the pyrimidine analogs is quenched; cleavage of the alpha-beta-phosphoryl bond of the UTP analog produces about a 14-fold increase in fluorescence intensity at 500 nm. Under the same conditions the fluorescence of the CTP analog increases about 8-fold, whereas the fluorescence of the purine analogs shows only a slight change. These derivatives are good substrates for Escherichia coli RNA polymerase with only slightly increased Km values and with Vmax values about 50 to 70% that of the normal nucleotides. They are used less efficiently by wheat germ RNA polymerase II. The ATP analog can be used by E. coli RNA polymerase to initiate RNA chains.     

The influence of arm length asymmetry and base substitution on the activity of the 10-23 DNA enzyme

A small oligodeoxyribonucleotide derived from in vitro selection has been shown to be capable of efficient sequence-specific cleavage of RNA at purine-pyrimidine junctions. As the reaction readily takes place under simulated physiologic conditions, this molecule described as the 10-23 general purpose RNA-cleaving DNA enzyme, has potential as a therapeutic agent. To further explore the character of this prototype, we examined the influence of base substitution and binding arm length asymmetry on its RNA cleaving activity. Surprisingly, substitution of the proximal nucleotide on the 3'-arm, to allow nonstandard Watson-Crick interactions, was found in some instances to improve the cleavage reaction rate. Although the identity of the unpaired purine in the RNA substrate cleavage site was found to have only a subtle influence on the rate of catalysis, with a slight decrease observed when a G at this position was changed to an A, nucleotide substitution (G to C) in the core motif at position 14 was found to completely abolish catalysis. The effect of arm length reduction varied with RNA substrate sequence and extent of helix asymmetry. Where the cleavage rate of one substrate was impaired by truncation of the deoxyribozymes 5'-arm (6 bp), the same modification in reactions with a different sequence produced a rate enhancement. Truncation of the 3'-arm, however, had no effect on the reaction rate of the one substrate tested yet nearly halved the cleavage rate in another substrate.