Probing the Dynamics of the P1 Helix within the Tetrahymena Group I

Jun 18, 2009 - Comparisons of the dynamics of the P1 helix in wild type and mutant ribozymes and in model constructs reveal a highly immobilized docke...
0 downloads 8 Views 2MB Size
Published on Web 06/18/2009

Probing the Dynamics of the P1 Helix within the Tetrahymena Group I Intron Xuesong Shi,† Emilia T. Mollova,† Goran Pljevaljcˇic´,‡ David P. Millar,‡ and Daniel Herschlag*,† Department of Biochemistry, Stanford UniVersity, Stanford, California 94305, and Department of Molecular Biology, The Scripps Research Institute, La Jolla, California 92037 Received April 7, 2009; E-mail: [email protected]

Abstract: RNA conformational transformations are integral to RNA’s biological functions. Further, structured RNA molecules exist as a series of dynamic intermediates in the course of folding or complexation with proteins. Thus, an understanding of RNA folding and function will require deep and incisive understanding of its dynamic behavior. However, existing tools to investigate RNA dynamics are limited. Here, we introduce a powerful fluorescence polarization anisotropy approach that utilizes a rare base analogue that retains substantial fluorescence when incorporated into helices. We show that 6-methylisoxanthopterin (6-MI) can be used to follow the nanosecond dynamics of individual helices. We then use 6-MI to probe the dynamics of an individual helix, referred to as P1, within the 400nt Tetrahymena group I ribozyme. Comparisons of the dynamics of the P1 helix in wild type and mutant ribozymes and in model constructs reveal a highly immobilized docked state of the P1 helix, as expected, and a relatively mobile “open complex” or undocked state. This latter result rules out a model in which slow docking of the P1 helix into its cognate tertiary interactions arises from a stable alternatively docked conformer. The results are consistent with a model in which stacking and tertiary interactions of the A3 tether connecting the P1 helix to the body of the ribozyme limit P1 mobility and slow its docking, and this model is supported by cross-linking results. The ability to isolate the nanosecond motions of individual helices within complex RNAs and RNA/protein complexes will be valuable in distinguishing between functional models and in discerning the fundamental behavior of important biological species.

Introduction

It has become increasingly apparent that RNAs are involved in numerous complex biological processes. Interestingly, many of these processes involve dynamic rearrangements of RNA structure. Examples include the helix rearrangements in spliceosome assembly, the processive catalytic cycle of protein synthesis by the ribosome, and the multistep self-splicing reaction carried out by more primitive RNA catalysts.1-5 RNA rearrangements are also integral to regulation of gene expression and accompany protein binding to regulatory elements and metabolite binding to riboswitches.6,7 While powerful approaches have been developed to follow conformational transitions in RNA, including rapid chemical †

Stanford University. The Scripps Research Institute. (1) Narlikar, G. J.; Herschlag, D. Nat. Struct. Biol. 1996, 3, 701–710. (2) Wilson, K. S.; Noller, H. F. Cell 1998, 92, 337–349. (3) Zhuang, X. W.; Bartley, L. E.; Babcock, H. P.; Russell, R.; Ha, T. J.; Herschlag, D.; Chu, S. Science 2000, 288, 2048–2051. (4) Smith, D. J.; Query, C. C.; Konarska, M. M. Mol. Cell 2008, 30, 657– 666. (5) The RNA World, 3rd ed.; Gesteland, R. F., Cech, T. R., Atkins, J. F., Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, 2006. (6) Legiewicz, M.; Badorrek, C. S.; Turner, K. B.; Fabris, D.; Hamm, T. E.; Rekosh, D.; Hammarskjold, M. L.; Le Grice, S. F. J. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 14365–14370. (7) Montange, R. K.; Batey, R. T. Annu. ReV. Biophys. 2008, 37, 117– 133. (8) Sclavi, B.; Sullivan, M.; Chance, M. R.; Brenowitz, M.; Woodson, S. A. Science 1998, 279, 1940–1943. ‡

10.1021/ja902797j CCC: $40.75  2009 American Chemical Society

probing and bulk and single molecule fluorescence,3,8-10 these approaches typically follow transitions between states separated by significant barriers for interconversion and not the underlying rapid conformational excursions of the individual ensembles. Furthermore, the tools available to follow such dynamic processes are limited. Recent elegant applications of NMR to RNA molecules of moderate size have led to detailed information about local motions on the nanosecond time scale.11,12 Nevertheless, these approaches are currently limited to moderately sized molecules at high concentration and cannot be applied to transient kinetic intermediates. Recently, applied sitespecific paramagnetic adducts can probe nanosecond dynamics at a particular site in a complex RNA13,14 and should also be highly useful in future studies. However, this approach also cannot interrogate transient intermediates and requires specialized experimental apparatus. Thus, additional dynamic probes of RNA would be valuable. (9) Joo, C.; Balci, H.; Ishitsuka, Y.; Buranachai, C.; Ha, T. Annu. ReV. Biochem. 2008, 77, 51–76. (10) Klostermeier, D.; Millar, D. P. Methods 2001, 23, 240–254. (11) Zhang, Q.; Stelzer, A. C.; Fisher, C. K.; Al-Hashimi, H. M. Nature 2007, 450, 1263–U14. (12) Zhang, Q.; Sun, X. Y.; Watt, E. D.; Al-Hashimi, H. M. Science 2006, 311, 653–656. (13) Edwards, T. E.; Sigurdsson, S. T. Nat. Protocol 2007, 2, 1954–1962. (14) Grant, G. P. G.; Boyd, N.; Herschlag, D.; Qin, P. Z. J. Am. Chem. Soc. 2009, 131, 3136–3137. J. AM. CHEM. SOC. 2009, 131, 9571–9578

9

9571

ARTICLES

Shi et al.

Traditionally, fluorescence has provided powerful probes of protein dynamics. Changes in fluorescence intensity or lifetimes provide evidence for conformational transitions, and the observed fluorescence changes are often related to changes in local dynamics, although those changes can be difficult to sort out. Fluorescence polarization anisotropy (FPA) provides a more direct readout of dynamics, as the degree of anisotropy is a function of the fluorophore’s motion over its fluorescence lifetime, typically in the low nanosecond regime. FPA has been applied extensively in studies of protein dynamics and the formation of protein complexes (e.g., refs 15-17). However, the use of FPA for nucleic acid studies has been severely limited, as described below.18-21 Complex RNAs are made up of short helical regions connected by junctions, and these regions are ultimately brought together or positioned by conformational preferences of the junctions, long-range tertiary interactions, and/or bound proteins. As RNA helices typically form rapidly and stably, the folding of complex RNAs often represents rearrangements of preformed helical elements. Thus, we would like to probe the dynamics of individual RNA helices within the context of larger RNA molecules and complexes. Fluorescent dyes that are base analogues and incorporated into RNA are typically severely quenched upon duplex formation22-25 and are therefore unusable for FPA studies of helices. The alternative of tethering a dye away from the helix is suboptimal in several respects. The flexibility in the tether decouples the motion of the dye from the helix and reduces the sensitivity of the anisotropy of the dye to the motion of the helix. The bulkiness of a tethered dye can also produce unwanted interactions or steric effects. Further, the fluorescence lifetime and anisotropy of such a dye can be affected by changes in the motions of the helix to which it is attached or by changes in the local environment of the dye. Attachment of base analogue dyes in single-stranded regions of RNA can be powerful probes of local conformational changes but cannot cleanly report on regional dynamics. Also, the previous approach of using a nucleic acid intercalator such as ethidium bromide18,19 is not generally tenable, as the dye is not localized to a single helix and may also distort the helical geometry and thus the overall RNA structure. In this study, we utilized the unusual property of a fluorescent guanosine analogue, 6-methylisoxanthopterin (6-MI; Figure 1A). This base maintains significant fluorescence when incorporated into duplexes,26 thereby allowing FPA measurements of individual helix dynamics. After we initiated our studies with this dye, another class of dyes with similar properties has been (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) 9572

Bucci, E.; Steiner, R. F. Biophys. Chem. 1988, 30, 199–224. Weber, G. Biochem. J. 1952, 51, 155–168. LiCata, V. J.; Wowor, A. J. Methods Cell Biol. 2008, 84, 243–262. Millar, D. P.; Robbins, R. J.; Zewail, A. H. J. Chem. Phys. 1981, 74, 4200–4201. Thomas, J. C.; Allison, S. A.; Appellof, C. J.; Schurr, J. M. Biophys. Chem. 1980, 12, 177–188. Chirico, G.; Collini, M.; Toth, K.; Brun, N.; Langowski, J. Eur. Biophys. J. 2001, 29, 597–606. Duhamel, J.; Kanyo, J.; DinterGottlieb, G.; Lu, P. Biochemistry 1996, 35, 16687–16697. Asseline, U. Curr. Org. Chem. 2006, 10, 491–518. Rist, M. J.; Marino, J. P. Curr. Org. Chem. 2002, 6, 775–793. Ward, D. C.; Reich, E.; Stryer, L. J. Biol. Chem. 1969, 244, 1228– 1237. Wilson, J. N.; Kool, E. T. Org. Biomol. Chem. 2006, 4, 4265–4274. Hawkins, M. E.; Pfleiderer, W.; Balis, F. M.; Porter, D.; Knutson, J. R. Anal. Biochem. 1997, 244, 86–95. J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

Figure 1. The Tetrahymena group I ribozyme (A) and docking equilibrium of its P1 duplex (B). Part B shows the tertiary structure model,30 colored as in part A. (C) Position of cross-linking sites between the P1 duplex and the ribozyme (Table 21,31).

reported27-29 and should be useful for studies analogous to those presented herein. We show that 6-MI accurately reports on overall nucleic dynamics using a series of model constructs, and we use this dye to interrogate the dynamics of a single RNA helix within the ∼400 nucleotide Tetrahymena group I ribozyme (Figure 1). Experimental Section DNA and RNA oligonucleotides made via solid-phase synthesis were obtained from Integrated DNA Technologies (Coralville, IA). 6-MI containing DNA and RNA oligonucleotides were obtained from TriLink (San Diego, CA) or Fidelity Systems (Gaithersburg, MD). All oligonucleotides were purified by gel electrophoresis or HPLC. Ribozymes were prepared by in vitro transcription with T7 RNA polymerase and purified by gel electrophoresis as previously described for the wild type and mutant L-16 ScaI ribozyme32 and for the L-21 ScaI ribozyme used in cross-linking experiments.31 Lac repressor (LacI) was expressed from plasmid pMDB1 (a gift from Michael Brenowitz, Albert Einstein College of Medicine) in E. coli and was purified as described.33 Assembly of Model Constructs and Ribozyme Complexes. Model constructs (Figure 3A) were assembled in two steps. The partially self-complementary long oligonucleotide, 5′-r(GGACAGGA(27) Sandin, P.; Borjesson, K.; Li, H.; Martensson, J.; Brown, T.; Wilhelmsson, L. M.; Albinsson, B. Nucleic Acids Res. 2008, 36, 157– 167. (28) Sandin, P.; Wilhelmsson, L. M.; Lincoln, P.; Powers, V. E. C.; Brown, T.; Albinsson, B. Nucleic Acids Res. 2005, 33, 5019–5025. (29) Wilhelmsson, L. M.; Holmen, A.; Lincoln, P.; Nielson, P. E.; Norden, B. J. Am. Chem. Soc. 2001, 123, 2434–2435. (30) Lehnert, V.; Jaeger, L.; Michel, F.; Westhof, E. Chem. Biol. 1996, 3, 993–1009. (31) Wang, J. F.; Downs, W. D.; Cech, T. R. Science 1993, 260, 504–508. (32) Karbstein, K.; Lee, J.; Herschlag, D. Biochemistry 2007, 46, 4861– 4875.

P1 Helix within the Tetrahymena Group I Intron

Figure 2. Predicted (black) and observed (gray) fluorescence polarization anisotropy of DNA duplexes of various lengths (Table 1A; see Supporting Information for description of duplexes and predictions).

GGG-X-AGUUA) d(GCGAATTGTGAGCGCTCACAATTCGC)r(UAACU)-3′, where X ) none, r(AAA), r(UUU), or -(CH2CH2-O)6-, was annealed to form a duplex with short hangover region on both sides by annealing at 95 °C for 5 min, cooling to 50 °C over 40 min, and quickly cooling to 4 °C. As the partially self-complementary long strand can also form stable hairpin structures, the purity of the annealed products was assessed by nondenaturing gel electrophoresis, which revealed negligible amounts (2-4%) of hairpin in all cases (data not shown). The annealed long strand duplexes were then hybridized with a 6-MI-labeled short 11mer RNA strand to form the complete model constructs (Figure 3A) to be used in FPA measurements. A mixture of long and short strands at a ratio of 4:1 was annealed at 38 °C for 30 min, and complete hybridization (>98%) was confirmed by nondenaturing gel electrophoresis (data not shown). Each model construct can bind up to two 6-MI labeled short strands. Further addition of duplexed long strands beyond the 4:1 ratio did not change the observed anisotropy or total fluorescence intensity (data not shown), providing additional evidence for complete hybridization and evidence against potential complexities from the fraction of complexes with two versus one fluorescent strand attached. Three oligonucleotide RNA strands were used: SF, 5′-(CCC UCC UFU CC)-3′, SFc, 5′-(CCC UCU UFU CC)-3′, and SFo, 5′-(CCCm UCU UFU CC)-3′, where “F” refers to the 6-MI deoxyribose residue and “m” refers to a 2′-OCH3 substitution. SFc and SFo are used in ribozyme experiments, as SFc and SFo favor the closed and open complexes, respectively (refs 34, 35; see Results and Discussion). SF was used in all model construct studies except for model constructs directly compared to ribozymes; in those cases, SFo was used. SF differs from SFc in only one residue: U6 of SFc forms a conserved G · U wobble base pair with the ribozyme’s internal guide sequence (IGS), and this G · U pair is important for docking.36-38 SF replaces U6 with C6 to form a G · C Watson-Crick base pair with the opposite strand. Complete binding of LacI to the model constructs was ensured by titrating concentrated LacI solution into hybridized model constructs until saturation (Figure S1). Experiments with LacI were carried out with a constant Na+ concentration of 300 mM. Extensive control experiments showed that this concentration of Na+ prevented any significant nonspecific binding of LacI, which could otherwise affect FPA values (Supporting Information, Figure S1; data not shown). (33) Hsieh, M.; Hensley, P.; Brenowitz, M.; Fetrow, J. S. J. Biol. Chem. 1994, 269, 13825–13835. (34) Pyle, A. M.; Cech, T. R. Nature 1991, 350, 628–631. (35) Herschlag, D.; Eckstein, F.; Cech, T. R. Biochemistry 1993, 32, 8299– 8311. (36) Knitt, D. S.; Narlikar, G. J.; Herschlag, D. Biochemistry 1994, 33, 13864–13879. (37) Strobel, S. A.; Cech, T. R. Science 1995, 267, 675–679. (38) Bartley, L. E.; Zhuang, X. W.; Das, R.; Chu, S.; Herschlag, D. J. Mol. Biol. 2003, 328, 1011–1026.

ARTICLES

Ribozyme complexes for FPA experiments were assembled as follows. The ribozyme was folded at 50 °C for 30 min in 50 mM NaMOPS, pH 7.0, and 10 mM MgCl2.39 A concentrated solution of the folded ribozyme was then titrated into 200 nM of short fluorescent substrate strands until saturation of the FPA signal. Typically, an excess of ∼2-fold was used. Saturation was confirmed before each experiment and was assumed to be maintained upon further addition of salt. Fluorescence Polarization Anisotropy (FPA). FPA measurements were performed on a Fluorolog-3 spectrometer using a T-format. Excitation and emission wavelengths were set to 350 and 425 nm, respectively. Each data point presented herein is the average of measurements from 2-4 independent samples. For each independent sample, samples were generally allowed at least 4 min to equilibrate before the first measurement, and a total of 3-4 consecutive measurements were made over ∼10 min, with a standard error that was generally smaller than 0.0015. The standard errors among different independent samples was generally less than 0.002. We estimate the total error in FPA value as (0.00152 + 0.0022)0.5 ) 0.0025, and the total error has been used as error limits in figures and calculations. All FPA experiments were carried out with 50 mM NaMOPS, pH 7.0, with specified amounts of added NaCl and MgCl2. The FPA of the model constructs were measured at 15 °C unless noted otherwise. The ribozyme measurements were carried out at 25 °C. Azidophenacyl Cross-Linking. The radiolabeled 5′-azidophenacyl L-21 ScaI ribozyme used in cross-linking studies was the same as that used in previous cross-linking studies and was prepared following the same protocol.1,31 Photo-cross-linking was carried out by irradiating for 10 min with a 312 nm Foto/Prep transilluminator (Fotodyne Inc., Hartland, WI) at room temperature. The cross-linking products were separated on a 6% denaturing acrylamide gel with 7 M urea and imaged using a Molecular Dynamics phosphorimager (Amersham Biosciences Inc., Sunnyvale, CA). Results and Discussion

As described in the Introduction, an ability to use FPA to assess motion of individual helices within a structured RNA or RNA/protein complex would greatly aid investigation of the properties and dynamics of these species. Required for such studies is a fluorescent base analogue that, when incorporated into a helix, maintains substantial fluorescence so that the dye will have the same motions as the helix and be able to report on these motions. Initial reports from Hawkins and co-workers suggested that 6-MI (Figure 1A), a fluorescent guanosine analogue that can base pair effectively with C residues, might be such an analogue.26 We first determined the fluorescent properties of 6-MI in simple duplexes to determine whether 6-MI FPA can accurately report on helix motions. We then extended these tests to model constructs (MC) of helices connected by single-stranded tethers. Having established the ability of 6-MI to report on helix motions, we then sought to isolate motions emanating from particular tether sequences to compare the behavior of model tethers with tethers of the same sequence embedded in a complex folded RNA, the group I ribozyme from Tetrahymena thermophila. Testing 6-MI as a Probe of Nucleic Acid Dynamics Using Simple Duplexes. Nucleic acid duplexes of moderate length are

rigid elements with defined, regular shapes so that dynamic properties can be predicted and compared to experimental measures.21,40 We therefore determined the FPA of 6-MI-labeled DNA duplexes of varying length and compared the experimental (39) Herschlag, D.; Cech, T. R. Biochemistry 1990, 29, 10159–10171. (40) Smith, A. L.; Cekan, P.; Rangel, D. P.; Sigurdsson, S. T.; Mailer, C.; Robinson, B. H. J. Phys. Chem. B 2008, 112, 9219–9236. J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

9573

ARTICLES

Shi et al.

Table 1. Sequences of Oligonucleotides Used in This Studya name

sequence

SF SFc SFo SC SO

(A) RNA Substrates or the Short Strand in Model Constructs 5′-r(CCC UCC UFU CC)-3′ 5′-r(CCC UCU UFU CC)-3′ 5′-r(CCCm UCU UFU CC)-3′ 5′-r(CCC UCU)-3′ 5′-r(CCCm UCU)-3′

MC-nicked MC-3A MC-3U MC-PEG

5′-r(GGA 5′-r(GGA 5′-r(GGA 5′-r(GGA

11mer 23mer 35mer 47mer

CAG CAG CAG CAG

GAG GAG GAG GAG

(B) The Long Strand of Model Constructs GGA GUU A)d(GC GAA TTG TGA GCG CTC ACA ATT CGC)r(UAACU)-3′ GGA AAA GUU A)d(GC GAA TTG TGA GCG CTC ACA ATT CGC)r(UAA CU)-3′ GGUUUA GUU A)d(GC GAA TTG TGA GCG CTC ACA ATT CGC)r(UAA CU)-3′ GG-(CH2-CH2-O)6-AGU UA)d(GC GAA TTG TGA GCG CTC ACA ATT CGC)r(UAA CU)-3′

(C) DNA Duplexes of Different Length 5′-d(CCC TCT TFT CC)-3′ 3′-d(GGG AGG ACA GG)-5′ 5′-d(AGG TTG ATT TTG CCC TCT TFT CC)-3′ 3′-d(TCC AAC TAA AAC GGG AGG ACA GG)-5′ 5′-d(TGT GTA AGT TTT AGG TTG ATT TTG/CCC TCT TFT CC)-3′ 3′-d(ACA CAT TCA AAA TCC AAC TAA AAC GGG AGG ACA GG)-5′ 5′-d(GCG GCT CCA ATG TGT GTA AGT TTT AGG TTG ATT TTG CGC TGT TFT CC)-3′ 3′-d(CGC CGA GGT TAC ACA CAT TCA AAA TCC AAC TAA AAC GGG AGG ACA GG)-5′

(A) 6-MI containing oligonucleotides used as substrates for the ribozyme or used as the short strand in model constructs (Figure 3A). F ) 6-MI, m substitution of 2′-OH with 2′-OCH3. (B) The long strand of the model constructs (Figure 3A). (C) A series of model DNA duplexes of different length (Figure 2). The “/” represents the position of the nick in the 35mer (see Supporting Information). a

values to predicted values (see Supporting Information). Duplexes ranging in length from 11 to 47 base pairs gave experimental FPA values that increase with helix length in a fashion that closely mirrors the predicted values, suggesting that 6-MI anisotropy can be used to monitor dynamic motion of helices (Figure 2). The observation of slightly smaller values relative to the predicted values may arise from residual motion of 6-MI within the helix, from a slight offset of the transmission dipole relative to a perpendicular from the helical axis,41 or from a combination of these factors. The duplexes also all have a single G · T mismatch (Table 1) that could lead to greater motion of 6-MI and lessen the observed anisotropy. Regardless of the origin of this small difference, the trend observed in Figure 2 and results described below provide strong support for the ability of 6-MI to report on helix dynamics. Further Testing 6-MI as a Probe of Nucleic Acid Dynamics Using Model Nucleic Acid Constructs with Intrinsic Flexibility.

A further prediction concerning the FPA of 6-MI-containing helices is that the introduction of single-stranded junction sequences should increase the dynamic motion, relative to constructs with the junctional residues engaged in base pairing or omitted. In addition, structured RNAs consist of short helical regions connected by formally single-stranded regions, so that constructs such as those studied herein can provide informative models for the constraints present for analogous sequences and structures within complex structured RNAs (see below). Figure 3A depicts a series of model constructs (MC) with two helices connected by different junctions. Constructs with single-stranded junctions of three A residues (MC-3A), three U residues (MC-3U), and a PEG tether of the length of three residues (MC-PEG) were compared to a construct lacking a tether (MC-nicked). Adding a single-stranded region in the model constructs is predicted to allow segmental motion and thus may lead to decreased anisotropy relative to a continuous helix. As predicted, the junction-containing constructs showed reduced anisotropies (Figure 4). (41) Lakowicz, J. R. Principles of Fluorescence Spectroscopy, 3rd ed.; Springer: New York, 2006. 9574

J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

Figure 3. Model constructs used to study junction dynamics. (A) Model

constructs (MC) with tethered RNA helices used in this study. The RNA portion of the chimera oligonucleotides is in blue; the DNA portion is in magenta; and the position of the 6-MI fluorescence probe is shown as an orange line (Table 1, see the Experimental Section for the assembly of the model constructs). (B) Binding of LacI reduces the overall tumbling of the constructs, isolating motions of the 6-MI-containing helix caused by the connecting junction.

The results of Figure 4 also agree with the expectation of greater stacking interactions, and thus less motion, for the construct containing a 3A junction than for the construct with a 3U junction.42,43 It is less intuitive that the 3U and PEG

P1 Helix within the Tetrahymena Group I Intron

Figure 4. Effect of junction identity on anisotropy. Anisotropy of model

oligonucleotide constructs (Figure 3A, Table 1) without LacI bound. The oligonucleotide used is SF (see text and Table 1A). The anisotropy for the 11mer duplex with the same sequence as the 11mer region of the model constructs is shown for comparison. Conditions: 50 mM NaMOPS, pH 7.0, 280 mM NaCl, 15 °C.

junctions should give similar dynamics. Small-angle X-ray scattering and electrostatic computations indicate that the helices for each of these constructs are largely coaxial under the ionic conditions used in this experiment, and the preferred coaxial orientation has been attributed to electrostatic repulsion between the helices.44,45 Perhaps the restrictions from electrostatic repulsion result in similar motions on the nanosecond time scale for the 3U and PEG constructs despite the greater intrinsic flexibility of PEG. This observation underscores the importance of probing dynamics on multiple time scales in future studies. The results of this and the previous section provide strong support for the ability of site-specifically incorporated 6-MI to probe the dynamics of individual helices. However, additional experiments gave paradoxical results that revealed a complicating factor for direct comparison of model construct mobility with mobilities observed for helices within large structured RNAs and RNA/protein complexes. We describe this complexity in the next section and the solution to it before probing the dynamics of a critical helical element of the Tetrahymena group I ribozyme. Isolation of Junctional Motion in Model Constructs Using a Protein “Anchor”. We were at first surprised to see that the

anisotropy of model constructs with the 3U and PEG junctions approaches that for the 3A construct with added ionic screening from Mg2+ (Figure 5B). However, results of small-angle X-ray scattering (SAXS) experiments with related constructs helped resolve this paradox and exposed the need for a helical “anchor”. SAXS studies revealed that constructs with 3U and PEG junctions more readily adopt a bent conformation than those with 3A junctions as ionic screening is increased.44,45 Apparently, the electrostatic relaxation from higher screening that allows bent conformers is opposed by the linear preference that arises from stacking of A residues.43,46 Thus, the observed FPA of the 3U and PEG constructs is a function of both the junction flexibility and the molecule’s overall tumbling and rotation, and (42) Tso, P. O. P.; Melvin, I. S.; Olson, A. C. J. Am. Chem. Soc. 1963, 85, 1289–1296. (43) Holcomb, D. N.; Tinoco, I. Biopolymers 1965, 3, 121–133. (44) Bai, Y. Ph.D. Dissertation, 2007. (45) Bai, Y.; Chu, V. B.; Lipfert, J.; Pande, V. S.; Herschlag, D.; Doniach, S. J. Am. Chem. Soc. 2008, 130, 12334–12341. (46) Brahms, J.; Michelso, A. M.; Van Holde, K. E. J. Mol. Biol. 1966, 15, 467–488.

ARTICLES

the overall motions are slower for a stably bent than for a linear molecule.20,47 The Mg2+-induced increase in FPA of the 3U and PEG constructs relative to the 3A construct may thus not represent a relative change in junction flexibility but rather differences in the overall shapes of the molecules. (There is also an effect of Mg2+ on the anisotropy of a duplex control, presumably because of an effect on local mobility of 6-MI; comparisons of FPA values obtained under different ionic conditions have therefore been normalized for this effect (see, e.g., Supporting Information, Figure S4).) Our aim with the model constructs, in addition to establishing that 6-MI FPA can be used to monitor helical motions, was to provide benchmarks for motions of individual helices within complex RNAs. We therefore needed to isolate the helix motion due to the junction flexibility in the model constructs. To do that, we incorporated the Lac repressor binding sequence in a central DNA region of the model constructs (Figure 3A). This sequence allowed binding of the LacI protein, which, due to its size (MW ) 154 kDa), greatly restricts overall tumbling of the model construct on the nanosecond anisotropy time scale (Figure 3B, Supporting Information, and Figure 5). As shown in Figure 5, addition of LacI increases the observed anisotropy. The Mg2+induced decrease in the anisotropy difference between A3 and U3/PEG is less significant with LacI bound than without, which corresponds to a reduction in impact of the molecular shape and overall tumbling on the observed anisotropy upon LacI binding. The residual Mg2+ effect could arise from different effects of Mg2+ on junction flexibility for the different junction sequences but would require further study to resolve. Exploration of Local Helix Dynamics for the P1 Duplex of the Tetrahymena Ribozyme. The Tetrahymena group I ribozyme

is one of the best studied complex RNAs.48-50 In the reaction catalyzed by this ribozyme, an oligonucleotide substrate (S) mimicks the 5′-splice site in the natural self-splicing form of the intron base pairs with the ribozyme’s internal guide sequence (IGS) to form the P1 helix (Figure 1A). The P1 helix is connected to the rest of the ribozyme through a single-stranded junction region called J1/2. The P1 duplex docks into tertiary interactions with the ribozyme’s conserved catalytic core. The states without and with these tertiary interactions are referred to as the “open” and “closed” complexes, respectively, and formation of the closed complex precedes the reaction’s chemical step.3,51,52 P1 docking is slow, with a rate constant on the order of 1 s-1, and is rate-limiting under some reaction conditions.38,51,53 The slow P1 docking is one of the remaining mysteries of the Tetrahymena ribozyme, and slow conformational steps appear to be common in RNA systems.38,54,55 To learn more about the docked and undocked complexes and how the mobility of the P1 helix might affect the docking process, we incorporated 6-MI into this helix. (47) Levene, S. D.; Wu, H. M.; Crothers, D. M. Biochemistry 1986, 25, 3988–3995. (48) Woodson, S. A. Curr. Opin. Struct. Biol. 2005, 15, 324–330. (49) Hougland, J. L.; Piccirilli, J. A.; Forconi, M.; Lee, J.; Herschlag, D. In The RNA World, 3rd ed.; Gesteland, R. F., Cech, T. R., Atkins, J. F., Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, 2006; pp 133-205. (50) Vicens, Q.; Cech, T. R. Trends Biochem. Sci. 2006, 31, 41–51. (51) Bevilacqua, P. C.; Kierzek, R.; Johnson, K. A.; Turner, D. H. Science 1992, 258, 1355–1357. (52) Herschlag, D. Biochemistry 1992, 31, 1386–1399. (53) Narlikar, G. J.; Bartley, L. E.; Khosla, M.; Herschlag, D. Biochemistry 1999, 38, 14192–14204. (54) Zarrinkar, P. P.; Williamson, J. R. Science 1994, 265, 918–924. (55) Zhuang, X. W.; Kim, H.; Pereira, M. J. B.; Babcock, H. P.; Walter, N. G.; Chu, S. Science 2002, 296, 1473–1476. J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

9575

ARTICLES

Shi et al.

Figure 5. Effect of LacI binding on anisotropy. Anisotropy of model oligonucletide constructs (Figure 3) without (black) or with (gray) LacI bound obtained in the absence (A) or presence (B) of 10 mM MgCl2. Conditions: 50 mM NaMOPS, pH 7.0, 280 mM NaCl, 15 °C. The oligonucleotide used is SF (Table 1A). The LacI experiments require a NaCl concentration of 280 mM or higher. Nonspecific binding of LacI with the model constructs becomes noticeable at lower NaCl concentration (see Supporting Information).

Figure 6. Comparison of the anisotropy of the L-16 ScaI ribozyme with

the P1 duplex containing 6-MI and model constructs (Figure 1A). Conditions: 50 mM MOPS, pH 7.0, 10 mM Mg2+ at 25 °C, with (A) 280 mM NaCl, (B) 0 mM NaCl, and (C) 280 mM added NaCl. R3A and R3U are ribozymes with J1/2 sequence of AAA and UUU, respectively. Oligonucleotide SFc is used for the closed complex, and SFo is used for open complex and model constructs. (SFc and SFo are defined in the Experimental Section and Table 1A.) (A) Comparison of anisotropy of the ribozyme’s P1 duplex in the closed complex with model constructs containing the same junction sequence of AAA. Under the salt condition in (A) (280 mM added NaCl), even with SFc as the substrate, P1 is partially undocked (∼52%38), which lowers the observed anisotropy value for the closed complex (after correction for the fraction undocked, the anisotropy of the docked P1 is calculated to be 0.296). (B) Comparison of anisotropy of the ribozymes in the closed and open complex. Under the salt condition in (B) (no added NaCl), even with SFo as the oligonucleotide, R3A is not completely in the open complex (Kdock ) 0.4 ( 0.1) (L. E. Bartley and D. H., unpublished results). The value in the figure is corrected accordingly. (The uncorrected value is 0.279.) (C) Comparison of anisotropy of the ribozyme’s P1 duplex in the open complex with model constructs containing the same junction sequence of AAA and UUU. The data for the model constructs with LacI bound in Figure 6 differ from Figure 5 because the experiments were done at different temperature and use different oligonucleotides (see Experimental Section and Table 1).

The first simple prediction is that the closed complex, with the P1 helix docked into its tertiary interactions, will exhibit low mobility. This prediction was confirmed, as the anisotropy of 6-MI incorporated in the P1 helix in the closed complex is considerably higher than the value for the same helix and same connecting sequence (J1/2 ) AAA) in model construct MC3A under identical conditions, 0.296 versus 0.243 (Figure 6A and legend). The model construct measurements for this and the other ribozyme comparisons were carried out with bound Lac repressor to eliminate contributions from the faster tumbling of the free model construct, and the model construct used for this comparison contained an internal G · U wobble pair to match the conserved G · U pair within the P1 helix (Figure 6 and Table 9576

J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

1). With Lac repressor bound, the overall size of the model construct is similar to the ribozyme (LacI + model construct ) 154 + 34 ) 188 kDa; ribozyme ) 130 kDa). One model for the slow docking of P1 is that the open complex is engaged in alternative nonproductive interactions that must be broken prior to correct docking. This model is supported by the formation of a specific cross-link between the P1 duplex and the ribozyme core from the open complex.1,31 However, we observed that the anisotropy of the open complex is considerably lower than that for the closed complex (Figure 6B) (the open complex anisotropy was determined by using an oligonucleotide with a 2′-methoxy substitution that disrupts the closed complex but has no effect on the docking rate constant or on the anisotropy of the remote 6-MI (ref 38; data not shown)). The low anisotropy value for the open complex provides evidence against models in which the open complex is trapped in a stable alternative docked conformer. The open complex is slightly less mobile than the corresponding model complex (Figure 6C), an effect that could arise from an interaction of one of the A residues of the J1/2 sequence. This model is consistent with the differential reactivity of the J1/2 A residues that is observed with dimethyl sulfate (DMS) in the ribozyme’s open complex.56 Also, local steric effects could alter the J1/2 conformation and mobility. However, an effect from local electrostatics is argued against by the observation that a ribozyme with a J1/2 junction consisting of three U residues does not exhibit higher anisotropy than its corresponding model control (Figure 6C and see below). Mirroring the greater mobility of U residues observed in the model constructs, a ribozyme with three U residues for its J1/2 sequence gave lower anisotropy for its open complex than that observed with the wild-type A3 sequence (Figure 6C). In this case, the ribozyme was slightly more mobile than the corresponding model construct. It is possible that there is residual motion within the ribozyme that increases the observed relaxation of the P1 helix relative to the model construct anchored by rigid binding to the Lac repressor. Consistent with this model is the observation that mutations that disrupt the ribozyme’s long-range tertiary interactions increase motion of the P1 duplex (i.e., decrease the anisotropy of the bound oligonucleotide; X. S., T. Benz-Moy, and D. H., unpublished results). To learn more about the open complexes, we compared crosslinking of ribozymes containing the wild-type A3 and the mutant (56) Engelhardt, M. A.; Doherty, E. A.; Knitt, D. S.; Doudna, J. A.; Herschlag, D. Biochemistry 2000, 39, 2639–2651.

P1 Helix within the Tetrahymena Group I Intron Table 2. Sites of Photo-Cross-Linking between the P1 Duplex with a 5′-Azido Group Attached and the Wild-type L-21 ScaI Ribozyme Reported in the Literature1,31 species

location

1 2b 2a

A114 U300 A88

* **

complex giving cross-link

closed complex no oligonucleotide substrate no oligonucleotide substrate, or open complex with oligonucleotide substrate a photoproduct that does not involve the 5′-azido group aggregation

U3 J1/2 sequences. Attachment of an azido-phenacyl cross-linker to the 5′-end of the L-21 ScaI version of the Tetrahymena group I ribozyme has been shown to give different cross-links in the open and closed complexes as well as with the free internal guide sequence (IGS) (Table 2; Figure 1C1,31). The cross-linking results support the conclusions from the anisotropy measurements of additional flexibility of the ribozyme containing the U3 J1/2 tether. The results with wildtype ribozyme reproduce earlier published results, indicating that different cross-linked species are preferentially formed in the unbound, open, and closed states (Figure 7, AAA1,31). With the U3-containing ribozyme, the closed complex does not form to an appreciable extent (L. E. Bartley and D. H., unpublished results), so that the same cross-linking pattern is observed with the oligonucleotides that favor and disfavor closed complex formation for the wild-type ribozyme (Figure 7, SC vs SO). Remarkably, the open complex of the U3-containing ribozyme gives all three major cross-links that are observed preferentially in the unbound, open, or closed states of the wild type ribozyme (Figure 7, SO). Further, the open complex cross-linking pattern for the mutant open complex resembles that of the unbound state (lanes 8 vs 10), consistent with extensive conformational mobility of this open complex. In contrast, for the wild-type

Figure 7. Photo-cross-linking products for 5′-azido-labeled L-21 ScaI ribozymes with AAA (R3A) and UUU (R3U) in the J1/2 junction. The species listed are defined in Table 2. Conditions: 50 mM NaMOPS, pH 7.0, 10 mM Mg2+, 23 °C, with varying concentrations of NaCl. Oligonucleotide SC, 5′-r(CCCUCU)-3′, favors the closed complex. Oligonucleotide SO, 5′-r(CCCmUCU)-3′, favors the open complex;34,35 “m” refers to a 2′OCH3 substitution. The varying NaCl concentration was used to destabilize residual closed complex.1,38 As expected, the intensity of band 1, the band associated with the closed complex in R3A, decreases with increasing NaCl concentration (lanes 3-5). However, the intensity of band 2a, the band associated with the open complex in R3A, does not increase detectably with increasing NaCl for R3A (lanes 3-5). For R3U, band 2a weakens with increasing NaCl concentration. A possible explanation is that increasing NaCl concentration renders the conformation associated with forming crosslinking at 2a less favorable as compared to conformations that give crosslinking at other positions or do not give cross-linking.

ARTICLES

ribozyme, cross-link 2b disappears and cross-link 1 appears in the open complex (lanes 2 vs 4). Conclusions and Implications

Docking of the P1 duplex occurs on the second time scale, much slower than expected for a free conformational search.38 The greater mobility of the P1 helix in the open complex as compared to the closed complex indicates that slow docking does not originate from trapping of P1 in an alternatively docked state. Nevertheless, the preferred cross-link formed in the open complex (Figure 71) is consistent with restricted motion in this state, as is the slightly lower anisotropy of this state relative to the corresponding model construct control. There is a modest increase in the rate of docking upon mutation of J1/2 to three U residues (L. E. Bartley and D. H., unpublished results). Correspondingly, when P1 is attached to U3 in the J1/2 mutant, there is an increase in the motion of the P1 duplex and the cross-linking positions that are accessible. Nevertheless, this J1/2 mutant destabilizes the docked state relative to the open complex (L. E. Bartley and D. H., unpublished results; see also Figure 7). On the basis of these results, we suggest the following model. The U3 J1/2 sequence allows additional conformer sampling of the P1 helix in the open complex, relative to the A3 sequence. This broader sampling increases the probability of collisions with regions of the ribozyme that make tertiary interactions with P1 in the closed complex, thereby increasing the rate of docking (kdock ) 3 s-1 for A3 and 6 s-1 for U3, L. E. Bartley and D. H., unpublished results). However, the docked complex is destabilized by the U3 mutation because tertiary and stacking interactions of the A3 sequence that are formed subsequent to the docking transition state cannot be made with the U residues (Kdock ) 30 for A3 and 1 for U3, L. E. Bartley and D. H., unpublished results). The change in protection of one of the A residues within J1/2 upon docking is consistent with an alternative preferred conformation of the wild-type open complex or with formation of tertiary interactions with J1/2 in the wild-type closed complex.56 Relating motions within complex RNAs to defined models will aid in analysis and interpretation, as will comparison to results from other experimental approaches. Herein, we have brought together fluorescence anisotropy results for wild type and mutant RNAs and model constructs with cross-linking, chemical protection, and single molecule fluorescence assays of P1 docking kinetics and equilibria. Although our finding in itself cannot fully explain why P1 docking happens on the second time scale, these results rule out a previous model for slow docking and lead to a more sophisticated docking model that integrates the results from these distinct experimental modalities. Our docking model is compatible with the previously proposed scenario in which the catalytic site is not prealigned for docking.38 Thus, even though the U3-tethered P1 duplex may attempt to dock more frequently than A3-tethered P1, the vast majority of the attempts for both species are likely to be unsuccessful because of the need of the catalytic core to rearrange, either prior to or following initial encounter with P1. Future challenges, in addition to testing and refining this model, are to better understand the molecular origins of RNA motions on the nanosecond time scale and the integration of techniques that provide dynamic and conformation information over a wide spectrum of time scales. We have explored the nanosecond time scale dynamics of an individual helix within different conformational states of a J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

9577

ARTICLES

Shi et al.

complex functional RNA. We suggest that 6-MI and related species26-29 will be powerful and incisive probes of RNA dynamics. Effective controls for overall tumbling effects, via the addition of protein anchors, and for local effects from different ionic environments can be carried out. Nonetheless, only an average motion over the nanosecond time scale is reported, and a full understanding of RNA conformational dynamics will require additional experimental modalities likely coupled to computation. Acknowledgment. We thank M. Brenowitz for the pMDB1 plasmid, L. E. Bartley and T. Benz-Moy for permission to cite unpublished results, and members of the Herschlag lab for

9578

J. AM. CHEM. SOC.

9

VOL. 131, NO. 27, 2009

comments on the manuscript. This work was supported by NIH Grants PO1 GM066275 and GM49243. Supporting Information Available: Saturation of LacI with model constructs, thermal and motional properties of 6-MI in helices, fluorescence lifetime measurement and base sequence requirements for the 6-MI experiments, FPA of DNA duplexes of different length, estimation of the anisotropy of simple duplexes with LacI bound, and normalization of FPA of model constructs. This material is available free of charge via the Internet at http://pubs.acs.org. JA902797J