To facilitate the measurement of intramolecular distances in solvated RNA systems, a combination of spin-labeling, electron paramagnetic resonance (EPR), and molecular dynamics (MD) simulation is presented. The fairly rigid spin label 2,2,5,5-tetramethyl-pyrrolin-1-yloxyl-3-acetylene (TPA) was base and site specifically introduced into RNA through a Sonogashira palladium catalyzed cross-coupling on column. For this purpose 5-iodo-uridine, 5-iodo-cytidine and 2-iodo-adenosine phosphoramidites were synthesized and incorporated into RNA-sequences. Application of the recently developed ACE® chemistry presented the main advantage to limit the reduction of the nitroxide to an amine during the oligonucleotide automated synthesis and thus to increase substantially the reliability of the synthesis and the yield of labeled oligonucleotides. 4-Pulse Electron Double Resonance (PELDOR) was then successfully used to measure the intramolecular spin–spin distances in six doubly labeled RNA-duplexes. Comparison of these results with our previous work on DNA showed that A- and B-Form can be differentiated. Using an all-atom force field with explicit solvent, MD simulations gave results in good agreement with the measured distances and indicated that the RNA A-Form was conserved despite a local destabilization effect of the nitroxide label. The applicability of the method to more complex biological systems is discussed.
RNA and its structural diversity has gained major attention in the last 10 years, in particular through the finding of RNAi as a natural antiviral mechanism of cells (
In order to rationally approach RNA as a 3D target, simple, fast and accurate methods to gain structural and dynamical information are necessary. Electron Paramagnetic Resonance (EPR) has already proved its efficiency in characterizing the structural environment of paramagnetic centers (
Despite these advances, each of the RNA spin labeling strategies mentioned above has its disadvantages and limitations. The major disadvantage of spin labeling RNA bases is the restriction to uridine. Either a 5-iodouridine (
Here, we report an extension of RNA base specific labeling to cytosine and the purine adenine, in addition to a considerable increase of the yield of TPA labeled RNAs using ACetoxyEthyl orthoester (ACE®) chemistry. PELDOR measurements for each duplex-RNA yielded a dipolar modulation, from which the distance between the spin labels distances could be extracted. Furthermore, molecular dynamics (MD) simulations gave results in good agreement with the measured distances and indicated that the TPA label induces only a small and local structural distortion.
As we already had learned from our DNA work ( Phosphoramidites prepared for incorporation of iodinated bases into protected RNA oligomers.
In addition to the pyrimidine bases U and C iodinated at the 5 position, we also modified a purine base A at the 2 position, which brings more flexibility in the choice of the labeled nucleotide and therefore of the spin-label position in the strand. This flexibility is particularly important for future EPR studies or NMR based ‘Paramagnetic Relaxation Enhancement’ measurements on biological systems. We took also into consideration the orientation of the nitroxide spin label in the RNA, which can be directed for duplexes either into the major (for U, C and potentially for A) or into the minor groove [for A and G (
Initially, we decided to use the current standard phosphoramidite chemistry with the acid-labile 4,4′-dimethoxytrityl group (DMT) and the fluoride-labile
The phosphoramidites Sonogashira cross-coupling during the solid-phase synthesis (example of 2-iodoadenosine).
Direct transfer of the procedure reported for DNA failed for RNA, in part due to its lower reactivity: three successive cross-couplings were necessary to achieve nearly quantitative yields in the case of 5-iodo-uridine in RNA (
As we recently have had excellent results with the newly developed ACE chemistry ( Synthetic pathway to 2-iodoadenosine phosphoramidite (ACE® chemistry); (
The 5′ and 3′ hydroxyl groups were simultaneously protected with the Markiewicz silyl group (
Phosphoramidites Spin-labeled RNA and their corresponding masses; the chemistry used for the preparation is indicated by a cross, a: RNA6 is non-self-complementaryRNA Sequence Method Calc. mass (g mol−1) Meas. mass (g mol−1) TBDMS ACE 1 3′ CGA CUA UAG UCG × 3958.6 3963.4 5′ GCU GAU AUC AGC 2 3′ CUG ACU AGU CAG × 3958.6 3962.1 5′ GAC UGA UCA GUC 3 3′ GC × 5244.3 5248.4 5′ CGA CUG AUA UCA G 4 3′ CGA C × × 3958.6 3960.0 5′ GCU GUA U 5 3′ GCA × 5244.3 5244.5 5′ CGU GUA UGC AUA 6a 3′ CGA G × 4938.1 4938.5 5′ GCU CAC UAU GU 4915.1 4920.6
To stress the advantage of using the ACE chemistry, we compared the results obtained by both methods for RNA (
To determine if the spin label disturbs the RNA structure, UV-melting and Circular Dichroism (CD) studies were performed. The UV-melting curves showed a destabilization of the duplexes between 1.5 and 5.1°C, slightly higher than for spin labeled DNA-duplexes (
We then subjected RNA-duplexes (
Also all other RNAs exhibit visible oscillations (supporting information) and the extracted mean distances are listed in RNA sequences and the corresponding N–O/N–O distances from PELDOR and Molecular Dynamics (MD) simulations The PELDOR distances are given along with the experimental error determined from the full width of the peak at half height. The MD distances are given including the width of the distance distribution in brackets. For RNA RNA Sequence 1 3′CGA C 19.3 ± 1.2 [DNA 23.3 ± 0.6] 18.0 (2.4) [DNA 21.4 (1.6)] 5′GCU GAU6 A 2 3′C 33.7 ± 3.9 [DNA 34.7 ± 1.4] 30.5 (2.4) [DNA 33.0 (2.7)] 5′GAC UGA UCA G 3 3′GC 38.7 ± 1.3 [DNA 44.8 ± 5.0] 36.2 (3.1) [DNA 43.3 (2.5)] 5′CGACUGAUAUCAG 4 3′CGAC 21.9 ± 0.8 24.7 (0.8) 5′GCUGUAU 5 3′GCA 33.6 ± 2.6 34.3 (1.8) 5′CGUGUAUGCAUA 6 3′GCAG 26.9 ± 1.3 24.6 (2.4) 5′GCUCACUAUGU
To be able to translate the measured N–O/N–O distances into RNA structures, we performed 50 ns all-atom MD simulations of all doubly labeled RNAs in explicit water solvent. The simulations yielded single-peaked distance distributions (measured between the oxygen atoms), whose mean and width are reported in Correlation of the PELDOR and MD distances for RNAs
Apart from the distances, MD simulations may provide detailed information on the structure and conformational dynamics of these RNA systems ( Molecular dynamics simulation results obtained for RNA
To study to what extent the spin labels affect the overall structure of RNA
In conclusion, the spin-label TPA was introduced into RNA by Sonogashira cross-coupling on column during oligonucleotide solid phase synthesis utilizing the bases A, U and C.
Application of the recently developed ACE chemistry presented the main advantage to limit the reduction of the nitroxide TPA into the corresponding amine during the oligonucleotide synthesis and thereby to increase substantially the reliability of the synthesis and the yield of labeled oligonucleotides. Thus, we are able now to introduce the spin label either into the major or minor groove of duplex RNA and to adenine, uridine and cytosine. The combination of site specific labeling with the advantage of the rigid label renders this method advantageous for distance measurements.
PELDOR experiments enabled us to measure the intramolecular distances in the six doubly labeled RNA-duplexes with significant modulations depth of up to 40% and allowed us to distinguish A-form RNA from B-form DNA duplexes. MD simulations on the same oligonucleotides gave results in good agreement with the measured distances and showed that the destabilization effect of the label is only local. Thus the combination of this spin-label strategy with PELDOR and MD opens a way to study structures in complex RNA folds and how they change upon binding of metals, small organic ligands or proteins.
The reactions were monitored by thin-layer chromatography (TLC) analysis on silica gel aluminum plates (silica gel 60
Commercially available from Glen Research, Sterling, VA, USA and Dharmacon, respectively.
All reactions were carried out under a protective argon atmosphere.
A mixture of 35 g (0.12 mol) guanosine, 70 ml (0.74 mol, 6.2 eq) acetic anhydride in 140 ml of dried dimethylformamide/pyridine 5/2 was heated at 75°C for 4 h. The resulting solution was cooled to 4°C, at which the product crystallized overnight. After filtration, washing with isopropanol and drying overnight, 46 g (91%) 2′,3′,5′-tri-
To a solution of pre-dried 2′,3′,5′-tri-
Iodine (11.3 g, 44 mmol), diiodomethane (36 ml, 10 eq), copper iodide (9.3 g, 49 mmol) and isopentyl nitrite (17.8 ml, 0.13 mol) were added to a solution of 19.0 g of 2′,3′,5′-tri-
To 750 ml of abs ethanol saturated with ammonia at 0°C was added 2′,3′,5′-tri-
A mixture of 2-iodo-adenosine (4.44 g, 11.3 mmol), 7.5 ml
A solution of
TEMED/HF was freshly and separately prepared: 0.89 ml TEMED (5.9 mmol) was dissolved in 2 ml acetonitrile at 0°C and 0.15 ml hydrofluoric acid (48% in water, 4.1 mmol) was added slowly in 2 min. This mixture was stirred for 5 min at 0°C and added dropwise in 5 min at RT to a solution of 1.08 g 2′-
Compound
About 216 μl Methyl-
About 1.50 g
To a solution of 1.82 g of 5′-
A solution of 705 mg
Acetylchloride (17.5 ml, 0.24 mol) was added to a solution of 10 g cytidine (40 mmol) in 65 ml acetic acid. After stirring overnight at RT, the solution was concentrated and the obtained white powder recrystallized in ethanol. Yield: 13 g (78%).
Iodine (4.11 g, 14.4 mmol) was added to a suspension of 2′,3′,5′-tri-
About 7.61 g 2′,3′,5′-tri-
5-Iodo-cytidine
A solution of
TEMED, HF was freshly and separately prepared: 0.85 ml tetramethylethylendiamine (5.7 mmol) was dissolved in 2 ml acetonitrile at 0°C and 0.14 ml hydrofluoric acid (48% in water, 4.1 mmol, 3.5 eq) was added slowly in 2 min. This mixture was stirred for 5 min at 0°C and added dropwise in 5 min at RT to a solution of 1.0 g 2′-
About 0.8 ml of a 0.45 M tetrazole solution in acetonitrile was added to a solution of methyl-
The oligonucleotides were synthesized on a 1 μmol scale on a EXPEDITE synthesizer from Perseptive Biosystems, Foster City, CA, USA, with phosphoramidites purchased from Biospring, Frankfurt am Main, Germany (TBDMS chemistry) or on a 0.2 μmol scale on a rebuilt ABI 392 synthesizer (Applied Biosystems, Foster City, CA, USA) with phosphoramidites purchased from Dharmacon (ACE chemistry). Every RNA synthesis was stopped after incorporation of the iodinated phosphoramidite without deprotecting the 5′-hydroxyl group (DMTon). The column was removed from the synthesizer and maintained under argon atmosphere. In the mean time 9.5 mg copper(I) iodide were dissolved in dried and deoxygenated CH2Cl2/Et3N (1.75/0.75 ml). About 150 μl of this solution were added under argon to a mixture of PdII(PPh3)2Cl2 (2.1 mg) and TPA (2 mg). The orange solution was given into the column and moved in it back and forth using two syringes. After a reaction time of 2.5 h the column was washed with 10 ml abs CH2Cl2, dried for 10 min under vacuum and flushed with argon. The Sonogashira cross-coupling was performed twice or three times depending on the labeled base and on the chosen chemistry for the oligonucleotide synthesis. Note that the amounts of reagents were not reduced for the 0.2 μmol synthesis. Then the column was reinstalled on the synthesizer to end the synthesis of the oligonucleotide.
When using the TBDMS chemistry, the oligonucleotides were cleaved from the controlled pore glass (CPG) and the amino groups deprotected with a mixture of ammonia (32%)/MeOH (3/1) over 24 h. The TBDMS groups were cleaved with triethylamine, HF over 24 h. After precipitation in abs. ethanol (−20°C over night), the RNA strands were purified via anion-exchange chromatography (Dionex NucleoPacTM PA 100 column, 250 × 9 mm, flow 5 ml min−1) on a JASCO-HPLC.
When using the ACE chemistry, the methyl group on the phosphate was first cleaved on-column with a 0.4 M solution of disodium-2-carbamoyl-2-cyanoethylene-1,1-dithiolate-trihydrate (S2Na2) in DMF/H2O: 98/2 in 30 min. Then, the oligonucleotides were cleaved from the solid support and deprotected with methylamine (40% in water): 10 min at 55°C for the unmodified RNA-strands and 12 h at RT for the spin labeled ones. Followed a purification through anion-exchange HPLC.
The oligonucleotides were desalted with PD-10 Sephadex columns from Amersham Biosciences, Piscataway, NJ, USA and finally characterized with a MALDI-Tof VOYAGER DE-PRO mass spectrometer from Applied Biosystems. In case of the ACE chemistry the final deprotection of the 2′-ACE groups was performed under sterile conditions with a TEMED-acetic acid buffer pH 3.8, 30 min at 60°C for all the RNA strands.
Calf intestine alkaline phosphatase (Sigma-Aldrich, St. Louis, MO, USA) and
All EPR samples had a volume of 100 μl and contained 0.1 mM duplex in phosphate buffer (140 mM NaCl, 10 mM Na2HPO4, 10 mM NaH2PO4, 20% ethylene glycol, pH 7). The solutions were transferred into sterile standard quartz EPR tubes and shock frozen in liquid nitrogen. The 4-Pulse ELDOR experiments were performed on an ELEXSYS E580 pulsed X-band EPR spectrometer using a flex line probehead housing a dielectric ring resonator all from Bruker, Fällanden, Switzerland. The temperature was adjusted with a temperature control system (ITC) in combination with a helium cryostat both from Oxford. For the PELDOR measurements a second microwave source was coupled into the microwave bridge using a commercially available set-up from Bruker. The pumping pulse was applied at the resonance frequency ν0 = νB of the resonator and the detection pulses at a frequency νA, 80 MHz higher than νB. The resonator used exhibits in overcoupled conditions a resonance frequency ν0 of 9.7 GHz, a quality factor Q of about 100, a conversion factor κ of 4 μT/W and a bandwidth of 97 MHz. Accordingly, even with a frequency offset of 80 MHz between νA = ν0 and νB both pulses are still within the bandwidth of the resonator. The pulse lengths used for the detection pulses were 32 ns and for the pumping pulse 12 ns. The amplitude of the detection pulses was chosen to optimize the refocused echo (global power 1 kW, attenuation in the +x-channel 8.0 (π/2-pulse), attenuation in the −x-channel 8.1 (π-pulses)). The amplitude of the inversion pulse with frequency νB was adjusted to a π-pulse using the pulse sequence π(νB-second source)-T-π/2(νB-main source)-τ-π(νB-main source) as far as possible, which is usually at about 8 dB of the second source. B0 was set to a field value so that the detection pulses excited the low field site of the nitroxide field sweep spectrum, whereas the inversion pulse selected the central
Tikhonov regularizations of the acquired time traces were performed with the program ‘DEER 2006’ from G. Jeschke available at
All MD simulations were performed with the GROMACS (
Details about the oligonucleotide analyses (enzymatic digestion,
Funding to pay the Open Access publication charges for this article was provided by Deutsche Forschungsgemeinschaft (SFB 579)
This work was supported by the SFB 579 RNA-ligand-interactions. Y.M. acknowledges support from a Singapore Ministry of Education (MoE), University Research Committee (URC, RG65/06) grant, and G. S. gratefully acknowledge support by the Frankfurt Center for Scientific Computing and the Fonds der Chemischen Industrie.