Deoxyribozymes: useful DNA catalysts in vitro and in vivo.

Deoxyribozymes (DNA enzymes; DNAzymes) are catalytic DNA sequences. Using the technique of in vitro selection, individual deoxyribozymes have been identified that catalyze RNA cleavage, RNA ligation, and a growing range of other chemical reactions. DNA enzymes have been used in vitro for applications such as biochemical RNA manipulation and analytical assays for metal ions, small organic compounds, oligonucleotides, and proteins. Deoxyribozymes have also been utilized as in vivo therapeutic agents to destroy specific mRNA targets. Although many conceptual and practical challenges remain to be addressed, deoxyribozymes have substantial promise to contribute meaningfully for applications both in vitro and in vivo.

An early transition state for folding of the P4-P6 RNA domain.

Tertiary folding of the 160-nt P4-P6 domain of the Tetrahymena group I intron RNA involves burying of substantial surface area, providing a model for the folding of other large RNA domains involved in catalysis. Stopped-flow fluorescence was used to monitor the Mg2+-induced tertiary folding of pyrene-labeled P4-P6. At 35 degrees C with [Mg2+] approximately 10 mM, P4-P6 folds on the tens of milliseconds timescale with k(obs) = 15-31 s(-1). From these values, an activation free energy deltaG(double dagger) of approximately 8-16 kcal/mol is calculated, where the large range for deltaG(double dagger) arises from uncertainty in the pre-exponential factor relating k(obs) and delta G(double dagger). The folding rates of six mutant P4-P6 RNAs were measured and found to be similar to that of the wild-type RNA, in spite of significant thermodynamic destabilization or stabilization. The ratios of the kinetic and thermodynamic free energy changes phi = delta deltaG(double dagger)/delta deltaG(o') are approximately 0, implying a folding transition state in which most of the native-state tertiary contacts are not yet formed (an early folding transition state). The k(obs) depends on the Mg2+ concentration, and the initial slope of k(obs) versus [Mg2+] suggests that only approximately 1 Mg2+ ion is bound in the rate-limiting folding step. This is consistent with an early folding transition state, because folded P4-P6 binds many Mg2+ ions. The observation of a substantial deltaG(double dagger) despite an early folding transition state suggests that a simple two-state folding diagram for Mg2+-induced P4-P6 folding is incomplete. Our kinetic data are some of the first to provide quantitative values for an activation barrier and location of a transition state for tertiary folding of an RNA domain.

Multiple folding pathways for the P4-P6 RNA domain.

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.

Quantifying the energetic interplay of RNA tertiary and secondary structure interactions.

To understand the RNA-folding problem, we must know the extent to which RNA structure formation is hierarchical (tertiary folding of preformed secondary structure). Recently, nuclear magnetic resonance (NMR) spectroscopy was used to show that Mg2+-dependent tertiary interactions force secondary structure rearrangement in the 56-nt tP5abc RNA, a truncated subdomain of the Tetrahymena group I intron. Here we combine mutagenesis with folding computations, nondenaturing gel electrophoresis, high-resolution NMR spectroscopy, and chemical-modification experiments to probe further the energetic interplay of tertiary and secondary interactions in tP5abc. Point mutations predicted to destabilize the secondary structure of folded tP5abc greatly disrupt its Mg2+-dependent folding, as monitored by nondenaturing gels. Imino proton assignments and sequential NOE walks of the two-dimensional NMR spectrum of one of the tP5abc mutants confirm the predicted secondary structure, which does not change in the presence of Mg2+. In contrast to these data on tP5abc, the same point mutations in the context of the P4-P6 domain (of which P5abc is a subdomain) shift the Mg2+ dependence of P4-P6 folding only moderately, and dimethyl sulfate (DMS) modification experiments demonstrate that Mg2+ does cause secondary structure rearrangement of the P4-P6 mutants' P5abc subdomains. Our data provide experimental support for two simple conclusions: (1) Even single point mutations at bases involved only in secondary structure can be enough to tip the balance between RNA tertiary and secondary interactions. (2) Domain context must be considered in evaluating the relative importance of tertiary and secondary contributions. This tertiary/secondary interplay is likely relevant to the folding of many large RNA and to bimolecular snRNA-snRNA and snRNA-intron RNA interactions.

RNA tertiary folding monitored by fluorescence of covalently attached pyrene.

The pathways by which large RNAs adopt tertiary structure are just beginning to be explored, and new methods that reveal RNA folding are highly desirable. Here we report an assay for RNA tertiary folding in which the fluorescence of a covalently incorporated chromophore is monitored. Folding of the 160-nucleotide Tetrahymena group I intron P4-P6 domain was used as a test system. Guided by the P4-P6 X-ray crystal structure, we chose a nucleotide (U107) for which derivatization at the 2'-position should not perturb the folded conformation. A 15-mer RNA oligonucleotide with a 2'-amino substitution at U107 was derivatized with a pyrene chromophore on a variable-length tether, and then ligated to the remainder of P4-P6, providing a site-specifically pyrene-labeled P4-P6 derivative. Upon titration of the pyrene-derivatized P4-P6 with Mg(2+), the equilibrium fluorescence intensity reversibly increased several-fold, as expected if the probe's chemical microenvironment changes as the RNA to which it is attached folds. The concentration and specificity of divalent ions required to induce the fluorescence change (Mg(2+) approximately Ca(2+) > Sr(2+)) correlated well with biochemical folding assays that involve nondenaturing gel electrophoresis. Furthermore, mutations in P4-P6 remote from the chromophore that shifted the Mg(2+) folding requirement on nondenaturing gels also affected in a predictable way the Mg(2+) requirement for the fluorescence increase. Initial stopped-flow studies with millisecond time resolution suggest that this fluorescence method will be useful for following the kinetics of P4-P6 tertiary folding. We conclude that a single site-specifically tethered chromophore can report the formation of global structure of a large RNA molecule, allowing one to monitor both the equilibrium progress and the real-time kinetics of RNA tertiary folding.

Energetics and cooperativity of tertiary hydrogen bonds in RNA structure.

Tertiary interactions that allow RNA to fold into intricate three-dimensional structures are being identified, but little is known about the thermodynamics of individual interactions. Here we quantify the tertiary structure contributions of individual hydrogen bonds in a "ribose zipper" motif of the recently crystallized Tetrahymena group I intron P4-P6 domain. The 2'-hydroxyls of P4-P6 nucleotides C109/A184 and A183/G110 participate in forming the "teeth" of the zipper. These four nucleotides were substituted in all combinations with their 2'-deoxy and (separately) 2'-methoxy analogues, and thermodynamic effects on the tertiary folding DeltaG degrees ' were assayed by the Mg2+ dependence of electrophoretic mobility in nondenaturing gels. The 2'-deoxy series showed a consistent trend with an average contribution to the tertiary folding DeltaG degrees' of -0.4 to -0.5 kcal/mol per hydrogen bond. Contributions were approximately additive, reflecting no cooperativity among the hydrogen bonds. Each "tooth" of the ribose zipper (comprising two hydrogen bonds) thus contributes about -1.0 kcal/mol to the tertiary folding DeltaG degrees'. Single 2'-methoxy substitutions destabilized folding by approximately 1 kcal/mol, but the trend reversed with multiple 2'-methoxy substitutions; the folding DeltaG degrees' for the quadruple 2'-methoxy derivative was approximately unchanged relative to wild-type. On the basis of these data and on temperature-gradient gel results, we conclude that entropically favorable hydrophobic interactions balance enthalpically unfavorable hydrogen bond deletions and steric clashes for multiple 2'-methoxy substitutions. Because many of the 2'-deoxy derivatives no longer have the characteristic hydrogen-bond patterns of the ribose zipper motif but simply have individual long-range ribose-base or ribose-ribose hydrogen bonds, we speculate that the energetic value of -0.4 to -0.5 kcal/mol per tertiary hydrogen bond may be more generally applicable to RNA folding.

Asymmetrical contributions of subunit pore regions to ion selectivity in an inward rectifier K+ channel.

We have investigated aspects of ion selectivity in K+ channels by functional expression of wild-type and mutant heteromultimeric G protein-coupled inward-rectifier K+ (GIRK) channels in Xenopus oocytes. Within the K+ channel pore (P) region signature sequence, a large number of point mutations in GIRK1 and GIRK4 subunits have been made at a key tyrosine residue--the "signature" tyrosine of the GYG. Studies of mutant GIRK1/GIRK4 heteromultimers reveal that the GIRK1 and GIRK4 subunits contribute asymmetrically to K+ selectivity. The signature tyrosine of GIRK1 can be mutated to many different residues while retaining selectivity; in contrast, the analogous position in GIRK4 must be tyrosine for maximum selectivity. Other residues of the P region also contribute to selectivity, and studies with GIRK1/GIRK4 chimeras reveal that an intact, heteromultimeric P region is necessary and sufficient for optimal K+ selectivity. We propose that the GIRK1 and GIRK4 P regions play roles similar to the two P regions of an emerging family of K+ channels whose subunits each have two P regions connected in tandem. We find different consequences between similar mutations in inward-rectifier and voltage-gated K+ channels, which suggests that the pore structures and selectivity mechanisms in the two classes of channel may not be identical. We confirm that GIRK4 subunits alone can form functional channels in oocytes, but we find that these channels are measurably permeable to Na2+ and Ca2+.

In vivo incorporation of unnatural amino acids into ion channels in Xenopus oocyte expression system.

A general method for the incorporation of unnatural amino acids into ion channels and membrane receptors using a Xenopus oocyte expression system has been described. A large number of unnatural amino acids have been incorporated into the nAChR, GIRK, and Shaker K+ channels. Continuing efforts focus on incorporating unnatural amino acids that differ substantially from the natural amino acids, for example, residues that include fluorophores. In addition, we are addressing the feasibility of incorporating unnatural amino acids into ion channels and membrane receptors in mammalian cells.

Flash decaging of tyrosine sidechains in an ion channel.

A nonsense codon suppression technique was employed to incorporate ortho-nitrobenzyl tyrosine, "caged tyrosine," in place of tyrosine at any of three positions (93, 127, or 198) in the alpha subunit of the muscle nicotinic ACh receptor (nAChR) expressed in Xenopus oocytes. The ortho-nitrobenzyl group was then removed by 1 ms flashes at 300-350 nm to yield tyrosine itself while macroscopic currents were recorded during steady ACh exposure. Responses to multiple flashes showed (1) that each flash decages up to 17% of the tyrosines and (2) that two tyrosines must be decaged per receptor for a response. The conductance relaxations showed multiple kinetic components; rate constants (<0.1 s(-1) to 10(3) s(-1)) depended on pH and the site of incorporation, and relative amplitudes depended on the number of prior flashes. This method, which is potentially quite general, (1) provides a time-resolved assay for the behavior of a protein when a mutant sidechain is abruptly changed to the wild-type residue and (2) will also allow for selective decaging of sidechains that are candidates for covalent modification (such as phosphorylation) in specific proteins in intact cells.

A regenerative link in the ionic fluxes through the weaver potassium channel underlies the pathophysiology of the mutation.

The homozygous weaver mouse displays neuronal degeneration in several brain regions. Previous experiments in heterologous expression systems showed that the G protein-gated inward rectifier K+ channel (GIRK2) bearing the weaver pore-region GYG-to-SYG mutation (i) is not activated by G beta gamma subunits, but instead shows constitutive activation, and (ii) is no longer a K(+)-selective channel but conducts Na+ as well. The present experiments on weaverGIRK2 (wvGIRK2) expressed in Xenopus oocytes show that the level of constitutive activation depends on intracellular Na+ concentration. In particular, manipulations that decrease intracellular Na+ produce a component of Na(+)-permeable current activated via a G protein pathway. Therefore, constitutive activation may not arise because the weaver mutation directly alters the gating transitions of the channel protein. Instead, there may be a regenerative cycle of Na+ influx through the wvGIRK2 channel, leading to additional Na+ activation. We also show that the wvGIRK2 channel is permeable to Ca2+, providing an additional mechanism for the degeneration that characterizes the weaver phenotype. We further demonstrate that the GIRK4 channel bearing the analogous weaver mutation has properties similar to those of the wvGIRK2 channel, providing a glimpse of the selective pressures that have maintained the GYG sequence in nearly all known K+ channels.

Dose-response relations for unnatural amino acids at the agonist binding site of the nicotinic acetylcholine receptor: tests with novel side chains and with several agonists.

Structure-function relations in the nicotinic acetylcholine receptor are probed using a recently developed method based on chemical synthesis of nonsense suppressor tRNAs with unnatural amino acid residues, site-directed incorporation at nonsense codons in Xenopus laevis oocytes, and electrophysiological measurements. A broad range of unnatural amino acids, as many as 14 at a given site, are incorporated at three sites, alpha 93, alpha 190, and alpha 198, all of which are tyrosine in the wild-type receptor and are thought to contribute to the agonist binding site. Confirming and expanding upon earlier studies using conventional mutagenesis, the three tyrosines are shown to be in substantially different structural microenvironments. In particular, a crucial role is established for the hydroxyl group of alpha Tyr93, whereas a variety of substituents are functional at the analogous position of alpha Tyr198. Interestingly, consideration of three different agonists (acetylcholine, nicotine, and tetramethylammonium) does not discriminate between these two best-characterized binding site residues. In addition, double-mutation studies establish the independent effects of mutations at the pore region (second transmembrane region) and at the agonist binding site, and this observation leads to a novel strategy for adjusting EC50 values. These results establish the broad generality and great potential of the unnatural amino acid methodology for illuminating subtle structural distinctions in neuroreceptors and related integral membrane proteins.

Subunit stoichiometry of a heteromultimeric G protein-coupled inward-rectifier K+ channel.

We investigated the stoichiometry of the heteromultimeric G protein-coupled inward-recitfier K+ channel (GIRK) formed from GIRK1 and GIRK4 subunits. Multimeric GIRK constructs with several concatenated channel subunits were expressed in Xenopus oocytes. Coexpression of various trimeric constructs with different monomers clearly showed that the functional channel has stoichiometry (GIRK1)2(GIRK4)2. Efforts to establish a preferred arrangement of subunits around the channel pore suggest that more than one arrangement may be viable.

Nicotinic receptor binding site probed with unnatural amino acid incorporation in intact cells.

The nonsense codon suppression method for unnatural amino acid incorporation has been applied to intact cells and combined with electrophysiological analysis to probe structure-function relations in the nicotinic acetylcholine receptor. Functional receptors were expressed in Xenopus oocytes when tyrosine and phenylalanine derivatives were incorporated at positions 93, 190, and 198 in the binding site of the alpha subunit. Subtle changes in the structure of an individual side chain produced readily detectable changes in the function of this large channel protein. At each position, distinct features of side chain structure dominated the dose-response relation, probably by governing the agonist-receptor binding.