High yield spin labeling of long RNAs for EPR spectroscopy

ABSTRACT: Site-directed spin labeling is a powerful tool for investigating the conformation and dynamics of biomacromolecules such as RNA. Here we int...
1 downloads 4 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Article

High yield spin labeling of long RNAs for EPR spectroscopy Mark Kerzhner, Hideto Matsuoka, Christine Wuebben, Michael Famulok, and Olav Schiemann Biochemistry, Just Accepted Manuscript • DOI: 10.1021/acs.biochem.8b00040 • Publication Date (Web): 01 May 2018 Downloaded from http://pubs.acs.org on May 2, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

High yield spin labeling of long RNAs for EPR spectroscopy Mark Kerzhner1, Hideto Matsuoka†, 2, Christine Wuebben2, Michael Famulok*,1,3 and Olav Schiemann*,2 1

Life & Medical Sciences Institute Chemical Biology & Medicinal Chemistry Unit c/o Kekulé-Institut für Organische Chemie und Biochemie University of Bonn, Gerhard-Domagk-Str. 1, 53121 Bonn, Germany 2 Institute of Physical and Theoretical Chemistry, University of Bonn, Wegelerstr. 12, 53115 Bonn, Germany 3

Max Planck Fellowship Chemical Biology Group, Stiftung caesar , Ludwig-Erhard-Allee 2, 53175 Bonn, Germany

ABSTRACT: Site-directed spin labeling is a powerful tool for investigating the conformation and dynamics of biomacromolecules such as RNA. Here we introduce a spin labeling strategy based on click-chemistry in solution that, in combination with enzymatic ligation, allows highly efficient labeling of complex and long RNAs with short reaction times and suppressed RNA degradation. With this approach a 34 nucleotides long aptamer domain of the preQ1 riboswitch and a 81 nucleotides long TPP riboswitch aptamer could be labeled with two labels in several positions. We then show that conformations of the preQ1 aptamer and its dynamics can be monitored in the absence and presence of Mg2+ and preQ1 ligand by continuous wave (cw) Electron Paramagnetic Resonance (EPR) spectroscopy at room temperature and pulsed electron-electron double resonance spectroscopy (PELDOR or DEER) in the frozen state.

INTRODUCTION RNA molecules are essential for many biological processes. Metabolite-sensing regulatory RNAs, or riboswitches, found in the non-coding regions of bacterial mRNAs control a broad range of genes in bacteria. In response to specific binding by cognate metabolites, riboswitches regulate protein-free gene expression. Due to their conformational variability, structural studies towards the dynamics of riboswitches are of great interest. Nuclear magnetic resonance (NMR)1-4 and Fluorescence Resonance Energy Transfer (FRET)5 are most commonly used to elucidate structural dynamics of riboswitches. Complementary to these methods, Electron Paramagnetic Resonance (EPR) spectroscopy in combination with site-directed spin labeling has also been employed to investigate local dynamics and conformational changes in riboswitches.6,7 Especially, pulsed EPR methods such as Pulsed Electron-Electron Double Resonance (PELDOR or DEER)8,9 enables one to measure distances in the range of 1.5 to 16 nm between sitespecifically incorporated spin labels in biomacromolecules10 and to follow their conformational dynamics.11,12 Two main strategies for site-directed spin labeling of RNA exist. One is the phosphoramidite method, in which spinlabeled phosphoramidites are incorporated at specific positions during the RNA solid phase synthesis.13,14,18 The other is postsynthetic spin labeling of reactive pre-functionalized sites in the RNA.15-18 The latter avoids exposing spin labels to the reducing environments during the oligonucleotide synthesis. This important advantage led to the development of various post-synthetic spin labeling approaches including the incorporation of paramagnetic nitroxides to all three sub-units of RNA: phosphate backbone,19,20 nucleobase21,22 and ribose sugar.23,24 For example, Grant et. al. used a deoxyribophosphorothiolate linkage in RNA as a reactive site for spin labeling with 1-oxyl-2,2,5,5-tetramethylpyrroline (R5).19 Another strategy, introduced by Ramos and Varani, attaches 3-(2iodoacetamido)-PROXYL (IAP) to 4-thiouridine modified

RNA.25 However, both spin labeling approaches are rather time consuming with reaction times of 12 – 20 h and have only limited yields. To improve the efficiency, a new approach was proposed on the basis of click chemistry on solid support and successfully applied to a 16 nucleotides long RNA in our previous work.26 The advantage of using click chemistry is that the reaction can be carried out in a straightforward way, provides high selectivity, specificity, high yields and convenient product purification. However, there are several disadvantages in spin labeling based on the solid support approach. Requiring 24 h of reaction time and two deprotection steps this method still remains rather tedious. In addition, all of the mentioned labeling methods are restricted to RNAs of lengths around 50 nucleotides (nt) that are amenable to chemical synthesis. This limit was overcome in some cases by spin labeling of two or more smaller RNA strands during solid phase synthesis followed by subsequent enzymatic ligation of the spin labeled RNA fragments. For example Duss et. al assembled a 70 kDa proteinRNA complex by spin labeling two 4-thiouridine modified RNA fragments and ligating them to each other with 20% – 40% efficiency using T4 DNA ligase.27 Thereby, a DTT-free ligation protocol was applied to prevent the reduction of the spin labels. Although this approach provides access to spin labeled RNA sequences of 70 nt or larger, the IAP spin label attached to the 4’-thio position drastically changes the hydrogen bonding pattern of the 4'-thiouridine base and changes the bonding situation within its ring system. Alternatively, Hoebartner and co-workers synthesized a 118 nt long SAM-I riboswitch containing a single spin-label, by site specific labeling of an 18 nt long RNA fragment at the exocyclic amino group of cytidine with 4-amino-TEMPO and subsequent ligation to an in vitro transcribed unmodified RNA using DNAzyme-catalysed ligation.7 The ligation yield of isolated spin labeled product amounted to 30% only and this approach

ACS Paragon Plus Environment

Biochemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

is restricted by the sequence requirements at the ends of RNA sites. Here, we employ click chemistry (Scheme 1a) in solution combined with a ligation strategy to provide access to long spin labeled RNAs in high yields. The effectiveness of the insolution click-labeling procedure is shown on the aptamer domain of a preQ1 riboswitch from Fusobacterium Nucleatum consisting of 34 nt (Scheme 1b) and its applicability for labeling long RNAs on the 81 nt long TPP riboswitch aptamer domain (Scheme 1c). The remarkably high yields that this approach achieves enabled not only cw EPR but also PELDOR studies on these functional RNAs. RESULTS AND DISCUSSION Spin labeling by click chemistry in solution Click-labeling of DNA has been reported previously28,29 and we showed recently that it can be extended to RNA with the oligonucleotide still on solid support and fully protected.26 Especially the protection of the chemically unstable 2’hydroxy group was crucial. However, click labeling would be easier and more straightforward to use if it were possible to perform on the fully deprotected RNA in solution. In a first step, this was tested on the aptamer domain of a class I preQ1 riboswitch, which specifically binds 7-aminomethyl-7deazaguanine (preQ1). The preQ1 riboswitch contains only 34 nt and thus can be fully synthesized by conventional solid phase RNA synthesis. As shown in scheme 1, the click reaction was carried out with the RNA functionalized with the alkyne group and the nitroxide label functionalized with the azide group. The reason for using the functionalities in this order is that the 5-ethynyl2’-deoxyuridine is commercially available and that the alkyne group is chemically more stable under the reducing conditions of automated RNA synthesis than the azide. In addition, the N3-functionalized label 2 can be easily synthesized.30 Furthermore, due to the previously reported partial oxidation of ethynyl moieties during alkaline nucleobase deprotection31, commercially available triisopropylsilyl (TIPS)-protected 5ethynyl-2’-dU phosphoramidites were used for RNA synthesis. The reaction was performed in the following way: First the catalytic Cu(I)/THPTA complex was prepared by mixing Cu(I) with an excess of tris(3-hydroxypropyltriazolylmethyl)amine (THPTA) in a solution of DMSO. The resulting mixture was incubated for five minutes at room temperature. Using an excess of THPTA ensured that no free Cu(I) was present, which could damage the RNA. This solution was then added to a DMSO/water solution of alkynylated RNA and spin label 2, and the resulting mixture was incubated for 30 min at 37 °C. Using the solvent mixture DMSO/water ensured that all reactants were soluble. In contrast to previous reports on clicklabeling of DNAs,32 degassing of the solvent was not necessary (Fig. S1). Using these reaction conditions, the singly spin labeled preQ1 constructs dU4, dU23 and dU32 as well as the doubly labeled ones dU4-dU32, dU4-dU23 and dU4-dU32(G11A) were prepared for cw EPR and PELDOR measurements. Mutant dU4-dU32(G11A) was chosen as a control because the point mutation G11A leads to a drastically reduced affinity of the riboswitch to the ligand preQ1.33 LCMS monitoring of the reactions revealed that all RNAs were fully spin-labeled after 30 minutes. This reaction time is considerably shorter then reported previously for click-labeling of DNA.29,32 Quantita-

Page 2 of 8

tive conversion of the starting material was achieved even for the doubly labeled constructs and RNA degradation was almost completely suppressed, as shown exemplarily for preQ1 dŲ4-dŲ32 in figure 1 (see Fig. S2 for the others).

Scheme 1 a) Spin labeling of 5-ethynyl-2’-dU modified RNA 1 with the azide-functionalized nitroxide 2 in solution using click chemistry. b) Secondary structure of the preQ1 aptamer in the folded state. Spin labeled and mutated positions are indicated as X1, X2, X3 and X4. c) Secondary structure of the TPP aptamer domain in the metabolite-bound state. Spin labeled positions are indicated as dŲ. The black arrow indicates the ligation site.

According to ESI mass spectrometric analysis only small traces of attached but reduced nitroxide could be observed. Pure spin labeled oligos where obtained after one HPLC run. Starting with 2.5 nmol RNA, 1.8 nmol of labeled RNA were obtained on average, which corresponds to a yield of about 72%. The identity of the modified oligos was confirmed by ESI mass analysis (Fig. 1 and Table S1). For the sake of comparison, the single mutant preQ1 dŲ32 and the double mutant preQ1 dŲ4-dŲ32 were also labeled by click reaction on solid support.26 For this approach, the RNA on solid support was incubated with the label and the catalyst for a total duration of 24 hours, during which the solid support was washed, dried and treated with freshly prepared reaction mixture three times. Even after these three successive spin labeling procedures, considerably lower conversion of the starting material (36%) was observed for both RNAs (Figure S3). This observation clearly demonstrates that not only the reaction time is shorter for the in-solution approach but also the handling is easier and the yields are twice as large. The in-solution labeling strategy was then applied to the 81 nt long TPP riboswitch (Scheme 1c). Due to the large RNA size two singly ethynyl pre-functionalized RNA fragments of the TPP aptamer, consisting of 45 nt (TPP-dŲ20) and 36 nt (TPP-dŲ68), respectively, were spin labeled quantitatively using click chemistry in solution (Fig. 2). After labeling the two fragments, they were assembled by splinted enzymatic T4

2 ACS Paragon Plus Environment

Page 3 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry nor effect of the dŲ label on the secondary structures of the preQ1 and TPP riboswitches (Fig. S4). In addition, only small decreases in Tm were observed, where preQ1 dŲ4 exhibits the largest decrease at 8 °C (Tab. 1), supporting the finding that the dŲ label has only a minor and likely just local effect on the preQ1 and TPP structures. For dŲ32, the Tm decrease by 5 °C correlates with the observation, that this label site is fairly occluded according to the mtsslWizard results (SI page S4).

Figure 1 HPLC and ESI mass analysis of preQ1 dŲ4-dŲ32 after 30 min of reaction. a) HPLC profile with the UV peaks recorded at 260 nm identified as follows: I = two reduced nitroxides, II = one reduced nitroxide, III = double spin labeled RNA. The correlation of UV signals is based on LCMS analysis. b) ESI mass spectrometric analysis of the UV peak III (collected from 9.710.1min). Mcalcd(preQ1 dŲ4-dŲ32) = 11376.78 m/z and Mfound(preQ1 dŲ4-dŲ32) = 11376.71 m/z.

DNA ligation analogously to previously published protocols.34,35 Afterwards, the doubly spin labeled full-length TPP aptamer was separated from the unreacted RNA fragments at 80 °C on an anion-exchange HPLC. Performing the ligation on a 5 nmol scale, 2.5 nmol of doubly labeled full-length TPP aptamer were obtained after purification (yield 50%). The purity of the doubly labeled product was verified by HPLC and ESI spectroscopic analysis (Fig. 2). The presence of dithiothreitol (DTT) in the ligation buffer did not affect the spin labels demonstrating the high stability of the dŲ label towards reduction.

Figure 2 HPLC profiles of a) TPP dŲ68 and of b) TPP dŲ20 after 30 min of reaction. c) HPLC profile of TPP dŲ20-dŲ68 after purification. The major UV signal in all HPLC chromatograms corresponds to the spin labeled RNA. d) ESI mass spectrometric analysis of TPP dŲ20-dŲ68. Mcalcd(TPP dŲ20-dŲ68) = 26596.17 m/z and Mfound(TPP dŲ20-dŲ68) = 26595.14 m/z.

Tm and CD control In a previous study reporting the influence of the dŲ label on the structure of a RNA duplex26 only minor differences were found between the unlabeled and labeled construct with respect to their Tm values and CD spectra. Also here, the CD spectra of the labeled riboswitches revealed only a slight decrease of the extrema at 210 nm and 260 nm indicating a mi-

Table 1. Differences in Tm-values for the spin labeled riboswitches as compared to the unlabeled ones.a

Construct preQ1 dŲ4 preQ1 dŲ23 preQ1 dŲ32 preQ1 dŲ4-dŲ23 preQ1 dŲ4-dŲ32 preQ1 dŲ4-dŲ32 (G11A) TPP dŲ20-dŲ68

∆Tm [°C]b w/o Mg2+

∆Tm [°C]b,c with Mg2+

-8 -3 -5 -2 -5 -5 -2

-4 +3 +1 -2 -6 -2 ±0

a

The absolute values are collected in the SI. Values are given per label. Melting curves of the constructs with the preQ1 ligand bound could not be recorded due to the UV absorbance of the ligand at 260 nm. b c

pre Q1 riboswitch: CW-EPR All labeled preQ1 riboswitch constructs display cw X-band EPR spectra at room temperature that are typical for partially immobilized nitroxide spectra and show no evidence for the presence of free label. Spin counting reveals quantitative spin labeling (Table S4). As an example, the cw EPR spectra of three samples of preQ1 dŲ32 (preQ1 dŲ32 alone, with Mg2+, and with Mg2+ plus the preQ1 ligand) are shown in Figure 3a (for the others see SI). Each of the three spectra is a superposition of two spectra, one with broad spectral width yielding the broad shoulder on the low field site and the other with narrower spectral width contributing mainly to the three intense lines. Only changing the rotational correlation time τc and the ratio of both spectra but using the same parameters otherwise can simulate both spectra. Figure 3 shows the best-fit simulated EPR spectra overlaid with the experimental ones. The narrower spectrum corresponds to a more mobile nitroxide label with a τc of 0.7 ns (Fig. 3, blue line) and the broader one to a more immobile one (Fig. 3, green line) with a τc of 6.3 ns. In the absence of Mg2+ and ligand, the fraction of the mobile species is 71% and that of the more immobile one 29% (Tab. 2). Upon addition of Mg2+ the fraction of the mobile species is reduced to 37% and further to 19% after adding the ligand. Simultaneously, the fraction of the immobile species increases to 63% and 81% after adding Mg2+ and then the ligand, respectively. This observation can be explained in accordance with the model of Breaker39 that was later confirmed by Micura using NMR36 in which the preQ1 riboswitch adopts a stem-loop structure with a single stranded 3’ overhang (Scheme 2, conformation I) in the absence of Mg2+.36 Upon addition of Mg2+, the riboswitch pre-organizes into a pseudoknot-like structure (conformation II). Subsequent addition of the preQ1 ligand leads only to small structural changes (conformation III).37

3 ACS Paragon Plus Environment

Biochemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 8

Figure 3 Room temperature cw X-band EPR spectra of a) preQ1 dŲ32 and b) preQ1 dŲ4-dŲ32(G11A) in the absence of preQ1 and Mg2+ (red), after addition of Mg2+ (black) and after addition of Mg2+ and preQ1 (blue). The signal intensity was normalized to the central peak height. The peak marked by an asterisk is due an impurity in the glass of the tube. Experimental (red) and calculated (black) CW-EPR spectra of c) preQ1 dŲ32 and d) preQ1 dŲ4-dŲ32(G11A) in the absence of preQ1 and Mg2+, after addition of Mg2+ and after addition of preQ1 at room temperature. The calculated spectra composed of two contributions with different mobility: more mobile (blue) and less mobile (green) spin labels. Conditions: cpreQ1 dŲ32 = 100 µM, cpreQ1 dŲ4-dŲ32(G11A) = 50 µM, 10 mM cacodylate buffer; cMg(II) = 2 mM, cpreQ1 = 1.2 mM.

In conformation I, the spin label at position U32 is located in a mobile, single stranded region. This conformation yields the narrow/mobile cw EPR spectrum, which is the dominating fraction in the absence of Mg2+. Adding Mg2+, this fraction reduces and the broader, immobile fraction increases in weight, which is assigned to conformation II. In conformation II, the single stranded overhang is bound to the loop and thus mobility at position U32 is greatly reduced giving rise to the EPR spectrum corresponding to the more immobile conformation. Obviously, this motional model assumes that pseudoknot formation leads to a spectrum that is dominated by slower global tumbling and that the local faster motion can be Table 2. Quantifications of mobile and immobile components of cw EPR spectra for the preQ1 constructs. w/o [%]a

Mg2+ [%]a

I

II+III

I

II+III

I

II+III

71 12 71 71 12 59

29 88 29 29 88 41

37 10 37 37 10 29

63 90 63 63 90 71

19 10 19 19 10 17

81 90 81 81 90 83

24

76

20

80

15

85

neglected.The small structural changes after addition of the ligand are not mirrored in the rotational correlation time. Instead, the ligand leads to a larger population of the two more immobile conformations II and III. Interestingly, the pseudoknot-like conformations are already populated in the absence of Mg2+. This has not been seen by room temperature NMR36 but maybe a consequence of the lower sensitivity of NMR with respect to EPR and the small extent to which these conformations are populated in the absence of Mg2+. The cw EPR spectra of preQ1 dŲ4 (Fig. S8) are very similar to that of preQ1 dŲ32 and the slow motion fraction also increases upon addition of Mg2+ and the ligand indicating that also position U4 can sense the structural rearrangement of the preQ1 riboswitch (Tab. 2).

preQ1 [%]a

preQ constructs dŲ4 dŲ23 dŲ32 dŲ4-dŲ23b

dŲ4-dŲ32c dŲ4-dŲ32 (G11A)d a

The number in front of the back slash denotes the fraction of the mobile, the number behind the slash, the fraction of the immobile component. b Simulations were done by assuming four components, fast and slow for each label site. The weight ratio was set to 50:50 for the two sites. c In these two cases only one mobile and one immobile component was sufficient to simulate the spectra, because preQ1 dŲ4 as well as preQ1 dŲ32 have identical τc values. d Simulations were done by assuming only two components, fast and slow.

Scheme 2 Schematic view of the conformational exchange proposed for the preQ1 riboswitch aptamer.36,38, 39 Spin labeled positions are depicted by dŲ.

This result is somewhat unexpected because position U4 locates within the stem region for which the movement of the label should be more restricted. However, the Tm studies revealed the largest drop of -8 °C in Tm for dŲ4, which may indicate that the stem is disrupted at this site leading to a higher mobility than expected. Although the U4 site is not involved in the loop/overhang interaction the sensing of the conformational rearrangement is traced back to the 3D structure (Fig. 4), where the folding in conformations II and III restricts the mobility of the label. Also for preQ1 dŲ23, the cw EPR spectra are a superposition of two spectra (Tab. 2 and

4 ACS Paragon Plus Environment

Page 5 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

Figure 4 Top: The derived distance distributions (red lines) for doubly labeled preQ1 dŲ4-dŲ32 a) without Mg2+ and preQ1, b) in presence of Mg2+ and c) after addition of preQ1. Each distance distribution was deconvoluted by two Gaussian components (gray and green dashed lines). Bottom: d) MtsslWizard distance distributions for preQ1-bound states of preQ1 dŲ4-dŲ23 (black) and preQ1 dŲ4dŲ32(blue). Next are shown the crystal structures of the preQ1 riboswitch for (e) preQ1 dŲ4-dŲ32 and (f) preQ1 dŲ4-dŲ23 in the preQ1bound state (PDB 3Q50), where the spin labels have been attached by means of the program mtsslWizard, are drawn as blue sticks with a red-coloured oxygen atom. Conditions: 20% deuterated ethylene glycol/buffer (10 mM cacodylate buffer), cRNA = 50 µM, cMg(II) = 2 mM, cpreQ1 = 1.2 mM.

Fig. S10). But the mobile and immobile species have slower rotational correlation times than found for dŲ32 and dŲ4. In addition, the 88% level of the fraction of the more immobile species is fairly high in the absence of Mg2+ and increases only marginally to 90% upon adding Mg2+ and preQ1. This steady level indicates that the label at position U23 is located within a rather rigid region that is not much influenced by the overall conformational change. Also these observations fit to the Micura model, where U23 locates right next to the stem and is not involved in the loop/single strand interaction. CW EPR measurements were also performed on the doubly labeled constructs preQ1 dŲ4-dŲ23 and preQ1 dŲ4-dŲ32 (Fig. S12 and S13). The spectra are a superposition of the spectra of the respective singly labeled constructs with slight differences in the fractions. To confirm that the observed spectral change after addition of Mg2+ and the ligand is due to the conformational change, cw EPR measurements were also performed on the mutant preQ1 dŲ4-dŲ32(G11A). Interestingly, its cw EPR spectrum contains a large fraction of immobile component already without Mg2+ and ligand, which only increases marginally upon addition of Mg2+ or Mg2+ plus preQ1 (Fig. 3b). This result is in nice agreement with a previous ligand binding study, which reported that the binding of the ligand was severely reduced because the loop overhang interaction is hindered.33 This supports the interpretation that the spectral change observed for the unmutated constructs is correlated with the conformational change. However, the reason why the immobile component dominates the spectra of preQ1 dŲ4-dŲ32(G11A) still remains unclear since the structure of this mutant is unknown.

preQ1 riboswitch: PELDOR If Mg2+ binding leads to a conformational change in the riboswitch, it should also be possible to observe this in PELDOR experiments. PELDOR measurements between spin labels yield label-to-label distance distributions, are wellestablished and have been used frequently on proteins11,40 but are rarely applied to riboswitches.6 In contrast to the cw EPR data, PELDOR data are collected in the frozen state. Previous PELDOR studies of a dŲ modified RNA26 duplex revealed orientation selection, however in the case of preQ1 orientation selectivity was found to be negligible most likely due to the higher conformational flexibility of preQ1 as compared to a duplex (Fig. S15). Figure 4 shows the PELDOR derived distance distributions for preQ1 dŲ4-dŲ32. The distance distribution of dŲ4-dŲ32 in the absence of Mg2+ can be deconvoluted with two Gaussians that peak at 1.9 nm and 2.8 nm with weights of 59% and 41%, respectively. In order to determine the appropriate conformation for the two Gaussian components, the X-ray structures of conformation II and III were labeled in silico with dŲ at positions U4 and U32 within the program mtsslWizard41-43 and the distance distributions were generated (Fig. 4). In the preQ1-bound state, conformation III, a mean distance of 2.7 nm was obtained, which fits nicely to the experimental distance at 2.8 nm, as shown in Figure 4e. For conformation II, the in silico labeling in mtsslWizard reveals that this site in conformation II is fairly occluded. This might indicate, that the introduction of the label at this site comes along with local structural changes. However, conformation II yields within mtsslWizard under these conditions a mean distance of 1.7 nm (Fig. S25), which fits nicely to the experimental distance of 1.9 nm. The secondary structure of

5 ACS Paragon Plus Environment

Biochemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

conformation I suggests a longer distance than for conformations II or III, but such longer distance was not observed in the PELDOR measurements. NMR showed that the pseudoknot formation occurs at lower temperature even in the absence of Mg2+,36 indicating that conformation I is mostly changed to conformation II at 50 K. After addition of preQ1, the weights of the two components at 1.9 nm and 2.8 nm change to 21% and 79%, respectively. The result indicates that only the fraction of conformation III increases after addition of preQ1, which is in accordance with the NMR data. As described above, the cw EPR measurements allowed us to discriminate between conformation I and II/III. On the other hand, the PELDOR measurements enabled to investigate the difference of conformation II and III. Table 3. Deconvolutions of the distance distributions for doubly labeled preQ1 mutants w/o Mg2+, with Mg2+ only and with Mg2+ and preQ1. w/o

preQ constructs dŲ4-dŲ23 dŲ4-dŲ32 dŲ4-dŲ32 (G11A)

Mg

preQ1

nm

%

nm

%

nm

%

II

1.9

100

1.9

100

1.9

100

II III II III

1.9 2.7 1.9 3

59 41 61 39

1.9 2.7 1.9 3

53 47 70 30

1.9 2.7 1.9 3

21 79 89 11

To verify the assignment of the two distance components, PELDOR measurements were also performed on the mutant preQ1 dŲ4-dŲ32(G11A). As shown in Fig. 5a, the resulting distance distribution does not change after adding Mg2+ and preQ1 ligand. This finding strongly supports the assignment made here because the G11A mutation inhibits ligand binding and pseudoknot-formation.33 After addition of Mg2+ and preQ1 a small distance shift of 0.2 nm is observed, which may be connected to a small structural change upon ligand binding but the change is too small as to be regarded as reliable. For conformation I, the result indicates that the distance between the spin labels might be outside the range of the current PELDOR measurements, or that the interspin distance is widely distributed due to the very high conformational flexibility. This speculation would be supported by the shallow modulation depth shown in Fig. S18b. PELDOR measurements were also performed for the preQ1 dŲ4-dŲ23 construct, which should be insensitive to the pseudoknot formation according to the cw EPR data and the model. Indeed, independent of the presence or absence of Mg2+ and preQ1, the peaks of the distance distributions were only observed around 1.8 nm, as shown in the Fig.4d, e and 5b.

Figure 5 Distance distributions derived from the PELDOR time traces for a) preQ1 dŲ4-dŲ32(G11A) and b) preQ1 dŲ4-dŲ23

Page 6 of 8

in the absence of preQ1 and Mg2+ (red), after addition of Mg2+ (black) and after addition of preQ1 (blue). Conditions: 20% deuterated ethylene glycol/buffer (10 mM cacodylate buffer), cRNA = 50 µM, cMg(II) = 2 mM, cpreQ1 = 1.2 mM.

TPP riboswitch For the TPP riboswitch, the cw EPR and PELDOR data are only shown as proof-of-principle and only for the riboswitch without Mg2+ and TPP ligand. The cw EPR spectrum of TPP dŲ20-dŲ68 in the absence of Mg2+ and TPP ligand is shown in Fig. 6a. The spectrum is typical for a partially immobilized nitroxide demonstrating the successful incorporation of the spin labels into the TPP riboswitch. Computer simulations revealed again two a more mobile and a more immobile fraction, where the fraction of the more mobile species has a τc of 5.0 ns and amounts to 63.5% while the fraction of the less mobile species has a τc of 101.9 ns and a weight of 36.5% (see Table S13 for parameters). The PELDOR measurement on the same construct yields the time trace in Figure 6b. The time trace reveals one long oscillation period with a modulation depth of 20%. The time window of 12 µs is long enough to yield reliable distances of up to 8 nm. In the obtained distance distribution, two broad peaks are observed at 5.1 nm and 7.2 nm, indicating the presence of at least two conformations in the absence of both, Mg2+ and ligand (Fig. 6b, inset).

Figure 6 a) Experimental (red) and simulated (black) cw X-band EPR spectra of TPP dŲ20-dŲ68 in the absence of both TPP and Mg2+ at room temperature. The calculated spectra are composed of a more mobile (blue) and a less mobile (green) spin label. b) The PELDOR derived distance distribution for the same construct. Conditions: 20% deuterated ethylene glycol/buffer (50 mM HEPES buffer), cRNA = 50 µM.

CONCLUSIONS The spin labeling of a preQ1 and a TPP riboswitch was accomplished with unprecedented efficiency by click chemistry in solution. The new spin labeling strategy enabled a significantly reduced reaction time, a high labeling efficiency, and a significant suppression of RNA degradation during the spin labeling procedure. For the preQ1 riboswitch, the conformational equilibrium and pseudoknot formation could be monitored by cw EPR and PELDOR measurements. The correlation between the conformational change of the preQ1 riboswitch and EPR/PELDOR results were confirmed through the EPR and PELDOR measurements of the preQ1 dŲ4-dŲ32(G11A) mutant. The application of the spin labeling strategy was extended to the large TPP riboswitch that cannot be synthesized on solid support in sufficiently high quantity. The distance measurement showed that two stable conformations for the TPP riboswitch exist

6 ACS Paragon Plus Environment

Page 7 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry even in the free state. Thus, cw EPR and PELDOR in combination with site directed spin labeling of RNA via clickchemistry is well suited for studying the conformational dynamics of large functional RNAs.

ASSOCIATED CONTENT Supporting Information Original PELDOR time traces, CD-spectra, Tm-values etc. This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION Corresponding Author *Email: [email protected]; [email protected] Present Addresses † Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan

Author Contributions All authors have given approval to the final version of the manuscript.

ACKNOWLEDGMENT Funding from the DFG via the SFB813 “Chemistry at Spin Centers” project A6, and by the Max-Planck Society is gratefully acknowledged.

REFERENCES (1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

Buck, J., Furtig, B., Noeske, J., Wohnert, J., and Schwalbe, H. (2007) Time-resolved NMR methods resolving ligandinduced RNA folding at atomic resolution, Proc. Natl. Acad. Sci. U S A. 104, 15699-15704. Lee, M. K., Gal, M., Frydman, L., and Varani, G.(2010) Real-time multidimensional NMR follows RNA folding with second resolution, Proc. Natl. Acad. Sci. USA. 107, 91929197 Noeske, J., Buck, J., Furtig, B., Nasiri, H. R., Schwalbe, H., and Wohnert, J. (2007) Interplay of 'induced fit' and preorganization in the ligand induced folding of the aptamer domain of the guanine binding riboswitch, Nucleic Acids Res. 35 , 572-583. Noeske, J., Richter, C., Grundl, M. A., Nasiri, H. R., Schwalbe, H., and Wohnert, J. (2005) An intermolecular base triple as the basis of ligand specificity and affinity in the guanine- and adenine-sensing riboswitch RNAs, Proc. Natl. Acad. Sci. USA 102, 1372-1377. Lemay, J. F., Penedo, J. C., Tremblay, R., Lilley, D. M., and Lafontaine, D. A. (2006) Folding of the adenine riboswitch, Chem. Biol. 13, 857-868. Krstic, I., Frolow, O., Sezer, D., Endeward, B., Weigand, J. E., Suess, B., Engels, J. W., and Prisner, T. F. (2010) PELDOR Spectroscopy Reveals Preorganization of the Neomycin-Responsive Riboswitch Tertiary Structure, J. Am. Chem. Soc. 132, 1454-1455. Buttner, L., Seikowski, J., Wawrzyniak, K., Ochmann, A., and Hobartner, C. (2013) Synthesis of spin-labeled riboswitch RNAs using convertible nucleosides and DNA-catalyzed RNA ligation, Bioorg. Med. Chem.. 21, 6171-6180. Milov, A.D., Salikhov, K. M., and Shirov, M. D. (1981) Application of Eldor in Electron-Spin Echo for Paramagnetic Center Space Distribution in Solids Fiz. Tverd. Tela 23, 975982. Martin, R. E., Pannier, M., Diederich, F., Gramlich, V., Hubrich, M., and Spiess, H. W. (1998) Determination of endto-end distances in a series of TEMPO diradicals of up to 2.8 nm length with a new four-pulse double electron electron resonance experiment, Angew. Chem. Int. Ed.37, 2833-2837.

(10) El Mkami, H., and Norman, D. G. (2015) EPR Distance Measurements in Deuterated Proteins, Methods Enzymol. 564, 125-152. (11) Jeschke, G. (2012) DEER Distance Measurements on Proteins, Annu. Rev. Phys. Chem. 63, 419-446. (12) Schiemann, O., and Prisner, T. F. (2007) Long-range distance determinations in biomacromolecules by EPR spectroscopy, Quart. Rev. Biophys. 40, 1-53. (13) Spaltenstein, A.., Robinson, B. H., and Hopkins, P. B. (1988) A Rigid and Nonperturbing Probe for Duplex DNA Motion, J. Am. Chem. Soc. 110, 1299-1301. (14) Bannwarth, W., and Schmidt, D. (1994) Oligonucleotides Containing Spin-Labeled 2'-Deoxycytidine and 5-Methyl-2'Deoxycytidine as Probes for Structural Motifs of DNA, Bioorg. Med. Chem. Lett. 4, 977-980. (15) Schiemann, O., Piton, N., Plackmeyer, J., Bode, B. E., Prisner, T. F., and Engels, J. W. (2007) Spin labeling of oligonucleotides with the nitroxide TPA and use of PELDOR, a pulse EPR method, to measure intramolecular distances, Nat. Protoc. 2, 904-923. (16) Schiemann, O., Piton, N., Mu, Y. G., Stock, G., Engels, J. W., and Prisner, T. F. (2004) A PELDOR-based nanometer distance ruler for oligonucleotides, J. Am. Chem. Soc. 126, 5722-5729. (17) Piton, N., Mu, Y., Stock, G., Prisner, T. F., Schiemann, O., and Engels, J. W. (2007) Base-specific spin-labeling of RNA for structure determination, Nucleic Acids Res. 35, 31283143. (18) Shelke, S.A, and Sigurdsson, S. Th. (2012) Site-Directed Spin Labelling of Nucleic Acids, Eur. J. Org. Chem., 2291– 2301. (19) Qin, P. Z., Butcher, S. E., Feigon, J., and Hubbell, W. L. (2001) Quantitative analysis of the isolated GAAA tetraloop/receptor interaction in solution: A site-directed spin labeling study, Biochemistry 40, 6929-6936. (20) Grant, G. P., Boyd, N., Herschlag, D., and Qin, P. Z. (2009) Motions of the substrate recognition duplex in a group I intron assessed by site-directed spin labeling, J. Am. Chem. Soc. 131, 3136-3137. (21) Qin, P. Z., Hideg, K., Feigon, J., and Hubbell, W. L. (2003) Monitoring RNA base structure and dynamics using sitedirected spin labeling, Biochemistry 42, 6772-6783. (22) Sicoli, G., Wachowius, F., Bennati, M., and Hobartner, C. (2010) Probing secondary structures of spin-labeled RNA by pulsed EPR spectroscopy, Angew. Chem. Int. Ed. Engl. 49, 6443-6447. (23) Schiemann, O, Weber, A, Edwards, T. E., Prisner, T. F, and Sigurdsson, S. T. (2003) Nanometer distance measurements on RNA using PELDOR, J. Am. Chem. Soc. 125, 3434-3435 (24) Saha, S., Jagtap, A. P., Sigurdsson, S. T. (2015) Site-Directed Spin Labeling of RNA by Postsynthetic Modification of 2'Amino Groups, Chem. Commun.51, 13142-13145. (25) Ramos, A., and Varani, G. (1998) A new method to detect long-range protein-RNA contacts: NMR detection of electron-proton relaxation induced by nitroxide spin-labeled RNA, J. Am. Chem. Soc. 120, 10992-10993. (26) Kerzhner, M., Abdullin, D., Wiecek, J., Matsuoka, H., Hagelueken, G., Schiemann, O., and Famulok, M. (2016) Post-synthetic Spin-Labeling of RNA through Click Chemistry for PELDOR Measurements, Chem. Eur. J. 22, 12113-12121. (27) Duss, O., Yulikov, M., Jeschke, G., and Allain, F. H. (2014) EPR-aided approach for solution structure determination of large RNAs or protein-RNA complexes, Nat. Commun. 5, 19. (28) Wu, P., and Fokin, V. V. (2007) Catalytic azide-alkyne cycloaddition: Reactivity and applications, Aldrichim. Acta 40, 7-17. (29) Ding, P., Wunnicke, D., Steinhoff, H. J., and Seela, F. (2010) Site-directed spin-labeling of DNA by the azide-alkyne 'click' reaction: nanometer distance measurements on 7-deaza-2'-

7 ACS Paragon Plus Environment

Biochemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(30)

(31)

(32)

(33)

(34)

(35)

(36)

deoxyadenosine and 2'-deoxyuridine nitroxide conjugates spatially separated or linked to a 'dA-dT' base pair, Chem. Eur. J. 16 , 14385-14396. Jakobsen, U., Shelke, S. A., Vogel, S., and Sigurdsson, S. T. (2010) Site-directed spin-labeling of nucleic acids by click chemistry: detection of abasic sites in duplex DNA by EPR spectroscopy, J Am Chem Soc. 132 (30), 10424-14288. Ingale, S. A., Mei, H., Leonard, P., and Seela, F. (2013) Ethynyl side chain hydration during synthesis and workup of "clickable" oligonucleotides: bypassing acetyl group formation by triisopropylsilyl protection, J. Org. Chem.78, 11271-11282. Haugland, M. M., El-Sagheer, A. H., Porter, R. J., Pena, J., Brown, T., Anderson, E. A., and Lovett, J. E. (2016) 2'Alkynylnucleotides: A Sequence- and Spin Label-Flexible Strategy for EPR Spectroscopy in DNA, J. Am. Chem. Soc. 138, 9069-9072. Rieder, U., Lang, K., Kreutz, C., Polacek, N., and Micura, R. (2009) Evidence for Pseudoknot Formation of Class I preQ(1) Riboswitch Aptamers, ChemBioChem. 10, 11411144. Lang, K., and Micura, R. (2008) The preparation of sitespecifically modified riboswitch domains as an example for enzymatic ligation of chemically synthesized RNA fragments, Nat. Protoc. 3, 1457-1466. Haller, A., Altman, R. B., Souliere, M. F., Blanchard, S. C., and Micura, R. (2013) Folding and ligand recognition of the TPP riboswitch aptamer at single-molecule resolution, Proc. Natl. Acad. Sci. USA 110, 4188-4199. Santner, T., Rieder, U., Kreutz, C., and Micura, R. (2012) Pseudoknot Preorganization of the PreQ(1) Class I Riboswitch, J. Am. Chem. Soc. 134, 11928-11931.

Page 8 of 8

(37) Jenkins, J. L., Krucinska, J., McCarty, R. M., Bandarian, V., and Wedekind, J. E. (2011) Comparison of a preQ1 riboswitch aptamer in metabolite-bound and free states with implications for gene regulation, J. Biol. Chem. 286, 2462624637. (38) Rieder, U., Kreutz, C., and Micura, R. (2010) Folding of a transcriptionally acting PreQ(1) riboswitch, Proc. Natl. Acad. Sci. USA 107, 10804-10809. (39) Roth, A., Winkler, W. C., Regulski, E. E., Lee, B. W. K., Lim, J., Jona, I., Barrick, J. E., Ritwik, A., Kim, J. N., Welz, R., Iwata-Reuyl, D., and Breaker, R. R. (2007) A riboswitch selective for the queuosine precursor preQ(1) contains an unusually small aptamer domain, Nat. Struct. Mol. Biol. 14, 308-317. (40) Yin, D. M., Hannam, J. S., Schmitz, A., Schiemann, O., Hagelueken, G., and Famulok, M. (2017) Studying the Conformation of a Receptor Tyrosine Kinase in Solution by Inhibitor-Based Spin Labeling, Angew. Chem. 129, 8537 – 8541. (41) Hagelueken, G., Ward, R., Naismith, J. H., and Schiemann, O. (2012) MtsslWizard: In Silico Spin-Labeling and Generation of Distance Distributions in PyMOL, Appl. Magn. Reson. 42 , 377-391. (42) Hagelueken, G., Abdullin, D., Ward, R., and Schiemann, O. (2013) mtsslSuite: In silico spin labelling, trilateration and distance-constrained rigid body docking in PyMOL, Mol. Phys. 111, 2757-2766. (43) Hagelueken, G., Abdullin, D., Ward, R., and Schiemann, O. (2013) mtsslSuite: In silico spin labelling, trilateration and distance-constrained rigid body docking in PyMOL, Mol. Phys. 111, 2757-2766.

SYNOPSIS TOC Click chemistry in combination with ligation is applied for site-selective and high yield spin labeling of large riboswitches. The incorporated nitroxide dŲ allows monitoring the conformational change within the preQ1 riboswitch using cw EPR and PELDOR spectroscopy. TOC

8 ACS Paragon Plus Environment