Cobalt-based paramagnetic probe to study RNA-protein interactions by NMR
Paramagnetic resonance enhancement (PRE) is an NMR technique that allows studying three-dimensional structures of RNA-protein complexes in solution. RNA strands are typically spin labeled using nitroxide reagents, which provide minimal perturbation to the native structure. The current work describes an alternative approach, which is based on a Co2+-based probe that can be covalently attached to RNA in the vicinity of the protein\xe2\x80\x99s binding site using \xe2\x80\x98click\xe2\x80\x99 chemistry. Similar to nitroxide spin labels, the transition metal based probe is capable of attenuating NMR signal intensities from protein residues localized <40 \xc3\x85 away. The extent of attenuation is related to the probe\xe2\x80\x99s distance, thus allowing for construction of the protein\xe2\x80\x99s contact surface map. This new paradigm has been applied to study binding of HIV-1 nucleocapsid protein 7, NCp7, to a model RNA pentanucleotide.
1. Introduction<!>2.1. Materials and instrumentation<!>2.2.1. Synthesis of 5<!>2.2.2. Synthesis of 6<!>2.2.3. Synthesis of 7<!>2.2.4. Synthesis of 1<!>2.3.1. Crystals of 1 (CCDC deposition number: 1438951)<!>2.4. Synthesis and preparation of modified RNA<!>2.5. Preparation and purification of protein<!>3.1. Synthesis and characterization<!>3.2. Incorporation of the probe into RNA pentanucleotide<!>3.3. Molecular dynamics simulations<!>3.4. PRE NMR experiments<!>4. Conclusions
<p>There is a strong interest in the biochemical community to develop a better fundamental understanding of how proteins interact with RNA [1]. RNA-protein interactions are ubiquitous in biology and critical at many regulatory steps of gene expression and stages of cell development [2,3]. X-ray crystallography has been the main tool in structural biology having paved way towards understanding a number of key ribonucleoprotein assemblies, such as ribosomes [4–6]. Nuclear magnetic resonance (NMR) spectroscopy is a complementary technique that, in addition to structure determination, allows studying dynamics interactions between biomolecules in solution and inside the live cells [7,8]. NMR is especially useful when the structural complexity precludes crystallization of a ribonucleoprotein complex.</p><p>The classical two-dimensional 1H and 15N correlation experiments (i.e. heteronuclear single quantum coherence (15N-HSQC)) correlate the 1H and 15N nuclei, giving rise to cross peaks reflective of labeled protein's environment [8]. The resonance frequencies of individual residues are determined by the amino acid specific chemical environments and reflect the overall 3-D conformation of an isotope labeled protein [9]. These two-dimensional NMR techniques have been used to study interactions between proteins and RNA aptamers [10–12]. Short distance atomic coupling, <6 Å, can be quantified to determine close range interactions in these compact structures. 15N-HSQC however does not provide long range structural information, which is important for larger or more flexible RNA-protein complexes.</p><p>Paramagnetic relaxation enhancement (PRE) is an NMR technique which can provide valuable structural information about RNA-protein complexes relying on distance constraints up to <40 Å away [13,14]. Extrinsic paramagnetic centers, incorporated onto RNA oligonucleotides via a linker, are capable of attenuating HSQC cross-peak signal intensities of the protein's backbone once the latter is attached to RNA. The extent of attenuation is related to the proximity of the protein's domains to the paramagnetic center, thus allowing construction of molecular maps describing protein's interactive surface. RNA has typically been modified with nitroxide spin labels due to their small size which minimally perturbs the native RNA structure [15]. The technique termed site-directed spin labeling (SDSL) has also allowed studying RNA's dynamic behavior in solution [16,17]. Nitroxide labels provide strong PRE effects due to slow electron relaxation, and NMR experiments which require these stronger PRE effects benefit from these types of paramagnetic tags. Clore and coworkers used a nitroxide spin-label to detect a rare conformational state in apo maltose-binding protein undergoing chemical exchange [18]. Qin and coworkers use SDSL to study dynamics of the substrate recognition duplex within the Tetrahymena ribozyme [19]. Additional reported nitroxide labels have provided important structural data related to RNA-protein binding in solution [20,21].</p><p>Paramagnetic metal ions can also be used to study long range molecular interactions with the benefits of additional versatility both through the choice of metal ion and ligand design. The latter can be engineered for site-specific incorporation into protein's metal binding domain [22]. Two-armed ligands can be employed for stronger and more rigid protein binding [23]. Analogous to nitroxide spin labels, paramagnetic centers with isotropic electron density, such as Mn2+, and Gd3+, can produce PRE effects which arise from slow electron relaxation that occurs during NMR experiments [24]. On the other hand, paramagnetic metal centers, such as high spin Co2+, Yb3+, and Pr3+, have anisotropic electron density which, in addition to PRE, can give structural information through pseudocontact shifts (PCS), and residual dipolar coupling (RDC).</p><p>The metal ion-based paramagnetic probes have been explored to study protein structure and dynamics. For example, a paramagnetic Cu2+ cyclen tag has been used to differentiate between two conformational states of human Ras protein [25]. Meanwhile, Ni2+ was utilized by Led and co-workers to obtain the structural information of thioredoxin [26,27]. Otting and co-workers genetically encoded Co2+-binding amino acid, bipyridylalanine, into West Nile virus NS2B-NS3 protease [22]. Work described herein is the first example where a paramagnetic transition metal-based probe is applied towards understanding binding interactions between RNA and proteins by NMR.</p><p>Inspired by the work of Graham [28], the Co2+-based paramagnetic NMR probe, shown in Scheme 2, has been designed for facile attachment to RNA strands using 'click' chemistry. The placement of the probe on the RNA has been strategically chosen to be in proximity to the protein's binding site and yet to be minimally disruptive of the RNA-protein interaction. Short linker reduces the probe's rotational flexibility thus simplifying interpretation of the NMR data [29]. Herein, we demonstrate that the Co2+-based probe is suitable for studying RNA-protein interactions by PRE NMR.</p><!><p>All chemicals were received from commercial sources and used without further purification. Chromatographic purifications were conducted using SiliaSphere™ spherical silica gel 5 μm, 60 Å silica gel (Silicycle). Thin layer chromatography (TLC) was performed on SiliaPlate™ silica gel TLC plates (250 μm thickness) purchased from Silicycle. Preparative TLC was performed on SiliaPlate™ silica gel TLC plates (1000 μm thickness). HPLC purification was performed using Phenomenex Luna 5μ C18(2) semi-preparative column (250 × 10 mm). 1H and 13C NMR spectroscopy was performed on a Bruker NMR at 400 (1H), and 100 (13C) MHz.</p><p>The RNA oligonucleotide synthesis was carried out on a 1.0 μmol scale using MerMade 8 DNA synthesizer. All the natural nucleoside phosphoramidites (TBDMS as the 2′-OH protecting group) and the accessory reagents were purchased from ChemGenes. After synthesis, the unmodified RNA oligo was cleaved from the beads and deprotected by the treatment with AMA solution (a 1:1 aqueous solution of methyl-amine and concentrated ammonium hydroxide) at 65 °C for 2 h. After evaporating the resulting solution to dryness, the 2′-TBDMS deprotection was performed using NEt3·HF solution at 65 °C for 2.5 h. The RNA was finally precipitated by using ammonium acetate and butanol, and re-dissolved in water.</p><p>The samples were analyzed on a Thermo Fisher Scientific (West Palm Beach, CA) LTQ Orbitrap Velos Mass spectrometer, using quartz capillary emitters. To facilitate spray optimization, 10% isopropyl alcohol was added to each sample prior to MS analysis. Cyclam and cyclen derivative were dissolved in MeOH and analyzed in the positive mode and the RNA was analyzed in a solution containing 150 mM ammonium acetate in the negative mode.</p><p>All 2-D NMR experiments were performed at RT on a 500 MHz Bruker Avance III NMR spectrometer equipped with the TCI cryoprobe. To prepare the NMR sample, 0.15 mM of [U-15N] HIV- nucleocapsid protein 7, NCp7, was dissolved in the NMR buffer 10 mM potassium phosphate, pH 6.5, and 10% D2O. RNA1 or RNA2 were added to the molar ratio of 1:1 between the NCp7 and RNA. 1H-15N HSQC with Watergate water suppression was used to monitor protein chemical shifts changes due to RNA binding [30]. 1024 × 256 points in proton and nitrogen dimensions, respectively, were acquired with 128 transients. The spectra were processed by using Topspin 2.1 (Bruker Inc) and analyzed by using Cara software [31]. The peaks were assigned based on published results [32]. The peak volumes were normalized by dividing the peak volume of each residue by the peak volume of E42 for NCp7-RNA1 complex.</p><p>The peaks intensity changes were quantified by using the ratio between the peak volumes: 
normalizedpeakvolumeofresiduefromRNA1normalizedpeakvolumeofresiduefromRNA2=RatioofPeakVolumeChange</p><!><p>A solution of Boc anhydride (2.95 g, 13.50 mmol) in anhydrous CH2Cl2 (60 mL) was added dropwise over 1 h to the solution of cyclam (1.50 g, 7.48 mmol) and DIPEA (6.50 mL, 37.3 mmol), in CH2Cl2 (250 mL) under nitrogen atmosphere. After addition, the mixture was cooled to −15 °C and a second solution of Boc anhydride (1.96 g, 8.98 mmol) in CH2Cl2 (60 mL) was added dropwise over 1 h. The solution was slowly warmed up to rt and stirred for 18 h. The solution was then washed with aqueous 0.5 M K2CO3 (2 × 150 mL), dried with Na2SO4 and concentrated under reduced pressure. The title product was obtained as a white foam by flash chromatography using a gradient of 0–100% EtOAc in hexanes, followed by 1–10% MeOH in EtOAc. RF = 0.32 in 1:9 MeOH:EtOAc with 0.1% Et3N) Yield = 1.77 g (47.32%); 1H NMR (CDCl3, 400 MHz) δ 3.39–3.28 (m, 12H), 2.79 (t, J = 5.5 Hz, 2H), 2.62 (t, J = 5.5 Hz, 2H), 1.92 (br s, 2H), 1.70 (br s, 2H), 1.45 (s, 27H); HRMS (DART) m/z: calcd. for C25H49N4O6 [M + H]+ 501.3652; found 501.3661. The observed spectra were consistent with the ones reported previously [33].</p><!><p>Sodium carbonate (1.50 g, 14.16 mmol) and propargyl bromide (~80% toluene, 268 μL, 3.54 mmol) were added to a solution of 5 (1.77 g, 3.54 mmol) in anhydrous acetonitrile (55 mL). The solution was stirred at reflux under nitrogen for 18 h. The insoluble salts were filtered and the solution was concentrated under reduced pressure. The title product was obtained as a white foam by flash chromatography using a gradient of 5% - 80% EtOAc in Hexanes. RF = 0.58 in 7:3 EtOAc:Hexanes. Yield = 0.98 g (51.31%); 1H NMR (CDCl3, 400 MHz) δ 3.40–3.29 (m, 14H), 2.67 (br s, 2H), 2.51 (t, J = 5.5 Hz, 2H), 2.16 (s, 1H), 1.89 (br s, 2H), 1.69 (br s, 2H), 1.46 (s, 27H); HRMS (ESI) m/z: calcd. for C28H51N4O6 [M + H]+ 539.3809; found 539.3778. The observed spectra were consistent with the ones reported previously [33].</p><!><p>Compound 6 (0.98 g, 1.82 mmol) was dissolved in a 1:1 mixture of TFA and CH2Cl2 (10 mL) and stirred for 1.5 h at room temperature. The solution was concentrated under reduced pressure and co-evaporated with MeOH (5 × 10 mL). The TFA salt was dissolved in a 1:9 mixture of MeOH:CH2Cl2 (50 mL) and washed with aqueous 1 M NaOH (150 mL), and brine (150 mL). The organic layer was dried with Na2SO4 and concentrated under reduced pressure. The title product was obtained as a light yellow solid.</p><p>Yield = 0.43 g (60.78%); 1H NMR (CDCl3, 400 MHz) δ 3.50–3.25 (m, 14H), 2.91 (t, J = 5.4 Hz, 2H), 2.84 (t, J = 5.5 Hz, 2H), 2.62 (s, 1H), 2.18 (t, J = 5.5 Hz, 2H), 1.93 (t, J = 5.5 Hz, 2H); 13C NMR (CD3OH, 100 MHz) δ 76.67, 76.46, 51.15, 50.00, 46.02, 45.73, 45.63, 45.23, 44.70, 42.49, 39.07, 23.59, 23.41; HRMS (ESI) m/z: calcd. for C11H23N4 [M + H]+ 211.1923; found 211.1903.</p><!><p>Co2+ chloride hexahydrate (59.7 mg, 0.252 mmol) was added to a solution of 7 (50.0 mg, 0.210 mmol) in EtOH (5 mL) and refluxed for1.5 h. The metal complex was precipitated with ether and filtered. The crude product was dissolved in MeOH (1 mL). X-ray quality crystals were obtained by slow diffusion of ether. HRMS (ESI) m/z: calcd. for C13H26ClCoN4 [Co(12)]3+ 166.0584; found 166.0575. IR (neat) 3159, 2982.05, 2937.76, 2515.59, 2030.18, 2105.99, 1977.50, 1610.60, 1454.87, 1422.51, 1383.28, 1253.78, 1240.95, 1136.13, 1107.25, 1092.25, 1065.06. 1H NMR (D2O, 400 MHz) δ 3.73 (t, J = 6.8 Hz, 2H), 3.58 (t, J = 6.9 Hz, 4H), 2.98–2.89 (m, 4H), 2.77 (t, J = 6.8 Hz, 3H), 2.53 (t, J = 6.9 Hz, 4H) 2.44–2.35 (m, 1H), 2.03–1.66 (m, 2H). Anal. Calcd. for C13H28Cl5CoN4 (M + 2HCl) C, 32.76; H, 5.92; N, 11.76. Found: C, 32.43; H, 6.06; N, 11.33.</p><!><p>Data collection was performed on a Bruker D8 VENTURE X-ray diffractometer with PHOTON 100 CMOS detector and mirror-monochromated Cu-Kα radiation (λ = 1.54178 Å) at T = 100(2) K. Data was corrected for absorption effects using the empirical method SADABS. The structures were solved with the SHELXTL-97 software package by direct methods and refined against F2 [34–36]. In general, all non-hydrogen atoms were located on the difference-Fourier map and refined anisotropically. Hydrogen atoms were placed in idealized locations and given isotropic thermal parameters equivalent to either 1.5 (terminal CH3 or NH3 hydrogen atoms) or 1.2 times the thermal parameter of the atom to which they were attached. At the final stages of refinement, the structures were checked for higher symmetry using PLATON [37].</p><p>All non-hydrogen atoms were refined anisotropically. Hydrogen atoms of two sumanenyl anions without disorder were located in the difference-Fourier map and refined with the C—H distances constrained to 1.01 ± 0.02 Å for CH and 0.95 ± 0.02 Å for CH2. All hydrogen atoms of the disordered sumanenyl dianion were included in idealized positions for structure factor calculations with Uiso(H) = 1.2 Ueq(C). This anion and K7 are disordered over two orientations that were resolved and treated with an occupancy ratio of 0.66:0.34. Both orientations were refined with anisotropic thermal parameters but all corresponding atoms' anisotropic thermal parameters were constrained to be the same.</p><!><p>Compound 2 was synthesized and immobilized on GE Healthcare Custom Primer Support™ Amino resin (GE Healthcare cat#17-5214-98) by following the previously reported procedure [38]. The functionalized resin was packed into empty Bioautomation MerMade columns (MM-1000-1) (4.00 mg resin per column). The RNA oligonucleotide synthesis was performed as described above. After synthesis, the RNA was cleaved from the resin as previously described [38]. After evaporating the resulting solution to dryness, the RNA was resuspended in DMSO (anhydrous, 100 μL). The 2′-TBDMS deprotection was performed using NEt3·HF (1:1, 125 μL) solution at 65 °C for 2.5 h. The RNA was cooled briefly and precipitated with 3 M NaOAc (25 μL), vortexed for 15 s, followed by butanol (1 mL) vortexed for 30 s. The solution was cooled at −80 °C for 30 min, then centrifuged for 15 min at 12500 rpm. The supernatant was removed and the precipitate was dried and resuspended in H2O. The synthetic RNA strands were purified by RP-HPLC using Phenomenex, Luna 5u C18(2) 100A column. Buffer A (100mM TEAA and 5% acetonitrile); Buffer B (100 mM TEAA 50% acetonitrile). The purification was achieved using an 18 min gradient of 0–85% Buffer B. After purification the solvents were lyophilized and the sample was desalted, re-suspended in water and quantified using nano-drop. 3 - HRMS (ESI) m/z: calcd. for C49H61N21O33P4 [M-2]2−797.6346; found 797.6353. RNA2 - HRMS (ESI) m/z: calcd. for C47H58N18O33P4 [M-2]2−763.1182; found 763.1222.</p><p>Click was done following the previously reported procedure [39]. Reaction progress was monitored by ESI-MS. Upon completion, the reaction mixture was lyophilized, resuspended in H2O, and purified by RP-HPLC using Phenomenex, Luna 5u C18(2) 100A column. Buffer A (100 mM TEAA and 5% acetonitrile); Buffer B (100 mM TEAA 50% acetonitrile). The purification was achieved using an 18 min gradient of 085% Buffer B. After purification the solvents were lyophilized and the sample was desalted, re-suspended in water and quantified using nano-drop. RNA1 - HRMS (ESI) m/z: calcd. for C62H85CoN25O33P4 [Co(RNA1)-4]2−944.6905; found 944.6933.</p><p>The clicked RNA strands were purified by RP-HPLC using Phenomenex, Luna 5u C18(2) 100A column. Buffer A (100 mM TEAA and 5% acetonitrile); Buffer B (100 mM TEAA 50% acetonitrile). The purification was achieved using an 18 min gradient of 0–85% Buffer B. After purification the solvents were lyophilized and the sample was desalted, re-suspended in water and quantified using nano-drop.</p><p>Prior to the 2-D NMR experiments, the RNA strands were ion exchanged by dissolving in 300 μL of aqueous 100 μM KCl. After 20 min equilibration, the solutions were loaded into SpinOUT™ columns (G biosciences GT-100, cat# 786-866) and the supernatants containing RNA were collected after centrifugation at 1000 G. RNA concentrations were determined by nano-drop.</p><!><p>NCp7 was isotope labeled by overexpression in BL21 pLysE cells grown in 15NH4Cl enriched media and purified by FPLC following the procedure from previous literature [32].</p><!><p>We designed an air-stable Co(III) complex, 1, that can be easily attached to an RNA strand of interest using click chemistry under aerobic conditions. Compound 1 is diamagnetic, thus allowing full characterization by NMR, in addition to mass spectrometry and X-ray crystallography. The metal can be reduced to paramagnetic Co(II) using sodium ascorbate after the probe is covalently attached to RNA and just prior to PRE NMR experiments [40]. The tetravalent macrocyclic ligand, cyclam, was chosen for its small size (relative to other known cobalt-binding ligands) and strong affinity to the metal (logK = 12.7 for Co(II)-cyclam) [41,42]. Macrocyclic ligands are known to provide higher kinetic stability to coordinating metals in comparison to their acyclic counterparts [43].</p><p>Compound 1 was synthesized in four steps from commercially available cyclam, as illustrated in Scheme 1. The synthesis commenced with selective protection of three secondary amine groups with Boc [33]. The remaining secondary amine group was alkylated with propargyl bromide. The Boc groups were cleaved using trifluoroacetic acid. Compound 1 was achieved upon refluxing CoCl2 in the ethanolic solution of 7 under aerobic conditions. Oxidation state of cobalt was evident based on the green color of the crystals which agrees with previous literature [44].</p><p>An X-ray crystallographic analysis was undertaken to confirm the identity of 1. Slow diffusion of ether vapor into 10 mM methanolic solution of 1 yielded green blocks suitable for X-ray analysis. The molecular structure, solved in monoclinic space group C 2/c, was formulated as [Co(7)Cl2]Cl. An ORTEP representation of the X-ray structure is shown in Fig. 1. The hexacoordinate Co3+ cation resides in a distorted octahedral geometry. The cobalt ion is coordinated in equatorial plane by the four nitrogen atoms of the macrocycle with interatomic distances of 1.965, 1.971, 1.975, and 2.066 Å. The apical positions are coordinated by the chloride ions with the observed interatomic distances of 2.2588 and 2.2678 Å. The observed crystallographic parameters are in agreement with previously reported Co3+ structure [44].</p><p>Diamagnetic Co3+ was reduced back to paramagnetic Co2+ using sodium ascorbate prior to NMR experiments. There is considerable amount of literature describing octahedral Co2+ cyclam complexes that were obtained by reduction of the corresponding Co3+ compounds. The reported magnetic susceptibility data suggests that octahedral Co2+ complexes containing modified cyclam ligands are typically in a high spin state with three unpaired electrons [45–47].</p><!><p>As a proof-of-principle, 1 was utilized to study binding of the RNA pentanucleotide.</p><p>5′-ACGCU*-3′ (U* represents the labeled nucleotide), and HIV-1 nucleocapsid protein, NCp7. The latter plays a number of key roles in the pathogenic lifecycle of the virus, including packaging of viral RNA into budding virions [48]. NCp7 has higher affinity for single stranded DNA or RNA than for double stranded DNA and has specificity for particular oligonucleotide sequences. The azide group was introduced at the 3′-end of the RNA strand using previously reported nucleoside, 2, immobilized on polystyrene resin [38]. The RNA strand, RNA1, having the sequence 5′-ACGCU*-3′ was synthesized by solid phase synthesis, as shown in Scheme 2. The paramagnetic probe was subsequently attached via 'click' chemistry using conditions that minimize Cu(l)-catalyzed hydrolysis of the RNA's phosphodiester backbone [39,49]. The later click reaction was done with excess sodium ascorbate to reduce cobalt back to its paramagnetic state, Co2+ [40,50,51].</p><p>There is a wealth of data describing oligonucleotide binding by NCp7, including an NMR structure of the protein's zinc finger domain bound to the single-stranded DNA pentanucleotide, 5′-ACGCC-3′ [52, 53]. The key binding interaction in this structure is pi-stacking between the indole ring of W25 inserted between the C2 and G3 bases, and the guanosine residue. Previously reported binding studies indicated that NCp7 binds equally tightly to the DNA and RNA pentanucleotide of the sequence described above [54,55]. This reported structure further predicts that attachment of the paramagnetic probe at the 3′-end of the pentanucleotide would place it in proximity of the protein, while minimally affecting its binding affinity. We replaced the cytidine residue at the 3′-end of the pentanucleotide by uridine to simplify the synthesis.</p><!><p>Molecular dynamics (MD) simulations, based on the published NMR structure of the NCp7-d(ACGCC) complex, were carried out to predict which protein residues would be most affected by the paramagnetic probe placed at the 3′-uridine. Ten parallel MD simulations initiated with different starting structures of the protein-nucleic acid complex were performed, generating a total of 1 microsecond simulation data. The simulation data was analyzed to assess the proximity of the paramagnetic probe to different amino acid residues of the protein. Fig. 2B shows the predicted interaction surface map of the RNA binding domain of NCp7 color-coded to reflect the proximity to the paramagnetic probe. A distance cut-off of 0.5 nm between the cobalt ion to the protein backbone was used for this purpose. Due to its attachment to the RNA's 3′-end nucleotide, which is not directly involved in protein binding, the paramagnetic probe will have some rotational flexibility bringing it within the contact distance with several sites on the protein surface, as shown in Fig. 2. From the figure it is evident that the amino acid residues R29, G35, and D48 would be in closest proximity to the probe and therefore will experience the greatest attenuation effects.</p><!><p>The proposed paradigm was tested by titrating 15N-labeled NCp7 with either paramagnetically tagged RNA1 or an unlabeled control RNA pentanucleotide, RNA2.15N-HSQC NMR spectra were obtained by titrating substoichiometric solutions of RNA to 500 μM [U-15N] NCp7 (450 μL) in 10 mM potassium phosphate buffer, 90% H2O/10% D2O. HSQC spectra of 15N-labeled NCp7 protein, as well as NCp7 bound to the unlabeled RNA were consistent with the ones reported previously [54,55].</p><p>Presence of the Co2+ probe resulted in global attenuation of the cross peak intensities. The extent of attenuation was found to be dependent on the predicted proximity of the protein's backbone to the paramagnetic center. Fig. 3A shows an overlay of the HSQC spectrum corresponding to NCp7 bound to RNA1 and NCp7 bound to the unlabeled RNA, RNA2. The peak volumes from the spectra were normalized against the residue most unaffected by the paramagnetic probe: E42. The extent of attenuation, illustrated in Fig. 3B, was quantitated by computing the ratios of each corresponding cross peak area between the two spectra.</p><p>The observed attenuation is consistent with the MD predictions. The strongest attenuation of intensity was observed for the G35 residue. Residues R29 and D48 were also attenuated, as predicted by the simulations. In addition, strong attenuation was observed for the following protein residues: K14, C15, N17, C18, A25, K26, N27, C28, H44, and R52. The graph bars were color coded to indicate the relative extent of attenuation indicative of the spatial proximity of the paramagnetic probe. Based on these results we were able to construct the protein's interactive surface map, shown in Fig. 3C, which is strikingly similar to that obtained by using molecular simulations (Fig. 2B). Residues shown in red are the most attenuated and thus the closest to the paramagnetic probe, while the ones shown in blue are the least attenuated.</p><p>The molecular map, shown in Fig. 3C, provides a visual representation of the residues most attenuated by proximity of the paramagnetic probe. Fig. 3D is a reduced representation of Fig. 3C, showing specifically the attenuation experienced by the N-H cross peaks from the protein backbone. Both of these describe a unique spectroscopic signature representing NCp7 bound to the paramagnetically labeled RNA1. Their distinguishing feature is the HSQC cross peaks whose intensities are attenuated as a function of proximity to the paramagnetic probe.</p><p>Paramagnetic metal-based probes can potentially offer three NMR experimental observables that yield long-range structural information: PREs, PCSs, and RDCs. While the PRE effect can be detected in any paramagnetic system, PCSs and RDCs can only be observed in systems with an anisotropic electron g-factor and there are a number of reports describing high spin Co(II)-based probes having PCS effects. The latter are dependent on the magnetic susceptibility tensor (usually referred to as the χ tensor). In general, if the electron g tensor is anisotropic, the χ tensor is also anisotropic [14]. The relationship between g-tensor and χ-tensor has been mathematically described [56]. The paramagnetic probe, 1, is attached to the RNA pentanucleotide through a flexible linker that allows the principal axes of the χ tensor to fluctuate within the frame of the macromolecule. In addition to the flexible linker, the model system also experiences a degree of conformational mobility from the oligonucleotide. The MD calculations, shown in Fig. 2, predicted proximity of the probe to several residues of NCp7. Because of the flexibility from both the linker and the oligonucleotide, the magnitude of the PCSs has been significantly reduced and only the PRE effects have been observed.</p><!><p>This work describes a new approach for studying RNA-protein interactions using two-dimensional correlation NMR experiments. We have shown that labeling a model RNA pentanucleotide with a Co(II)-based paramagnetic NMR probe allowed for creation of visual maps describing the RNA's binding to HIV-1 nucleocapsid protein, NCp7. The probe resulted in distance-dependent attenuation of HSQC crosspeaks corresponding to amide backbone of specific amino acid residues of NCp7 protein. The signal attenuation observed in the two-dimensional NMR experiments was consistent with the predictions obtained by MD simulation. The structural data obtained from this proof-of-concept study is in agreement with the previously described NOE experiments that elucidated binding of the unlabeled pentanucleotide, 5′-ACGCC-3′, and NCp7 [55]. This serves as a confirmation that the paramagnetic probe was attached just outside of the binding site, having minimal impact on the RNA-protein binding. In the future, the Co(II)-based paramagnetic probe could be used to study more complex systems involving RNA-protein interactions with contact distances further than the NOE distance constraint of 6 Å.</p><p>Future work will be focused on redesigning the probe to reduce its rotational flexibility. A Co(II)-based probe with a rigid attachment will be able to provide PCS data, in addition to PRE. Our ultimate goal is to apply the paramagnetic metal-based NMR probes to study how NCp7 binds viral RNA in live human cells [57,58]. To achieve this, the probe will have to be covalently attached at a strategically chosen internal positions of significantly larger RNA strands, such as 20-nt long stem loop SL3 [57].</p>