PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of springeropenLink to Publisher's site
Journal of Biomolecular Nmr
 
J Biomol NMR. 2013 January; 55(1): 33–46.
Published online 2012 November 23. doi:  10.1007/s10858-012-9683-9
PMCID: PMC3555346

Database proton NMR chemical shifts for RNA signal assignment and validation

Abstract

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the 1H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson–Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 43 possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA 1H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.

Electronic supplementary material

The online version of this article (doi:10.1007/s10858-012-9683-9) contains supplementary material, which is available to authorized users.

Keywords: RNA, Chemical shift, A-form helices, NMR signal assignment and validation

Introduction

RNA molecules participate in a large and expanding array of known biological functions including gene regulation, maintenance of sub-cellular and viral structure, intracellular trafficking, antiviral restriction, catalysis, and, of course, propagation of genetic information (Korostelev and Noller 2007; Steitz 2008; Bessonov et al. 2008; Boisvert et al. 2007; Wakeman et al. 2007; Edwards et al. 2007; Bartel 2004; Kim 2005; Hassouna et al. 1984; Brodersen and Voinnet 2006; Doudna and Rath 2002; Ponting et al. 2009). Like proteins, the functional activities of most RNAs are intrinsically linked to their structures. Unfortunately, although a wealth of structural information is currently available for functionally active proteins and protein domains, structural information for functionally relevant RNAs remains relatively limited. Thus, the Protein Data Bank (PDB; http://www.rcsb.org/pdb/home/home.do) currently contains more than 55,000 protein structure depositions, whereas the Nucleic Acid Database (NDB; http://ndbserver.rutgers.edu/) contains atomic coordinate depositions for fewer than ~2,100 RNAs and protein/ligand-RNA complexes, of which ~1,600 were determined by X-ray crystallography and ~500 by NMR spectroscopy. Conformational heterogeneity and the presence of a relatively uniform, negative surface charge can hinder structural studies by X-ray crystallography, and as discussed below, difficulties associated primarily with limited chemical shift dispersion have generally limited NMR applications to relatively small RNAs. For these reasons, much of what is known about the structures of biologically functional RNAs (primarily secondary structure information) has been obtained by chemical and enzymatic accessibility mapping experiments, coupled with phylogenetic and free energy calculations. Although RNA probing methodologies are potentially very powerful and have been widely applied (Peattie and Gilbert 1980; Ehresmann et al. 1987; Stern et al. 1988; Forconi and Herschlag 2009; Weeks 2010), interpretation of the data can be problematic, particularly for RNAs that exist as equilibrium mixtures of multiple conformational species (see for example, Kladwang et al. 2011; Houck-Loomis et al. 2011; Lu et al. 2011a,b; Miyazaki et al. 2010).

NMR is a potentially powerful tool for probing RNA structure (Wüthrich 1986; Allain and Varani 1997; Lukavsky and Puglisi 2005), but its application to larger RNAs can be complicated by a number of factors. Inter-residue scalar couplings are generally weak, limiting the utility of “through bond” inter-residue connectivity experiments for signal assignment. The most commonly used assignment approach involves identification of sequential inter-residue NOE connectivities (Wüthrich 1986), but even this approach can be problematic for modest sized RNAs (ca. 25–60 nucleotides). Although resolution can be increased by 1H–13C heteronuclear spectral editing (Peterson et al. 2004; D’Souza et al. 2004; D’Souza and Summers 2004; Davis et al. 2005; Batey et al. 1995; Batey et al. 1992; Nikonowicz and Pardi 1992; Nikonowicz et al. 1992; Michnicka et al. 1993; Kim et al. 1995; Xu et al. 1996; Kim et al. 2002; Lukavsky et al. 2003; Lu et al. 2009), chemical shift dispersion is relatively limited (Allain and Varani 1997; Lukavsky and Puglisi 2005), and severe dipolar broadening of the aromatic 1H–13C signals that are critical for structural analysis can preclude detection of 1H–13C correlation NMR signals in larger RNAs (Lu et al. 2011a). In addition, interproton distances between elements of secondary structure in larger RNAs typically exceed those required for NOE detection (Lu et al. 2009; Tolbert et al. 2010). Thus, high-resolution NMR-based structural studies have been applied mainly to relatively small RNAs: Of the 496 RNA NMR structures that have been deposited in the NDB, only 19 contain 60 or more nucleotides; the largest is a symmetrical dimer of 132 nucleotides (two 66 nucleotide subunits), and the average size is ~27 nucleotides.

One approach for addressing issues of signal degeneracy involves the application of traditional 2D NOESY experiments to RNA samples that are site- and/or nucleotide-specifically labeled with deuterium (Miyazaki et al. 2010; D’Souza et al. 2004; Davis et al. 2005; Kim et al. 1995; Lu et al. 2009; Zhou et al. 2006; Nelissen et al. 2008; Heng et al. 2012; Duss et al. 2012). 2H-isotope edited 2D NMR has enabled nearly complete assignment of the aromatic, H1′, H2′, and H3′ ribose signals of RNAs containing up to 132 nucleotides (Miyazaki et al. 2010), and has also enabled assignment of selected residues within a 720 nucleotide RNA (Lu et al. 2011a; Heng et al. 2012). This approach, which involves comparison of high resolution 2D NOESY spectra obtained for multiple, differentially 2H-labeled samples, avoids relaxation problems associated with aromatic 1H–13C spectral editing and enables observation of signals in 2D 1H–1H NOESY spectra for protons with T2 values as short as 8 ms (Lu et al. 2011a). Although resolution and sensitivity can be improved dramatically by nucleotide-specific deuteration, signal overlap can still hinder the assignment process for RNAs comprising more than 150 nucleotides (Summers and coworkers, unpublished).

NMR chemical shifts have been widely utilized for NMR signal assignment and structural studies of proteins (for examples see: Grzesiek and Bax 1993; Wishart and Sykes 1994; Wishart et al. 1991, 1992; Cavalli et al. 2007; Shen et al. 2008; Wishart et al. 2008). Although relationships between 13C chemical shifts and RNA structure have been identified (Ebrahimi et al. 2001; Fares et al. 2007; Ohlenschlager et al. 2008), and 15N NMR chemical shifts have been incorporated into a probabilistic approach for automated assignment of RNA imino groups (Bahrami et al. 2012), heteronuclear NMR chemical shifts have not been widely exploited for RNA studies (Lam and Chi 2010; Aeschbacher et al. 2012). On the other hand, Wijmenga and co-workers showed that non-exchangeable 1H NMR chemical shifts for A-form helical residues could be back-calculated from a given 3D RNA structure (Cromsigt et al. 2001). For 28 examples tested, the back-calculated shifts were in good agreement with shifts reported in the Biological Magnetic Resonance Bank (BMRB; www.bmrb.wisc.edu), and some general 1H NMR chemical shift trends were identified (Cromsigt et al. 2001). Here we report a detailed analysis of the H8, H2, H6, H5, H1′, H2′, and H3′ proton NMR chemical shifts that have been deposited in the BMRB. After correcting for differences in chemical shift referencing and sample conditions, excellent correlations were observed, despite the fact that the data were obtained over a wide range of sample conditions. Our findings confirm and quantify previously identified trends and identify new sequence- and structure-dependent chemical shift correlations that can be used for assignment and/or validation of non-exchangeable 1H NMR chemical shifts and for the identification of non-canonical RNA structural features and intermolecular interaction sites.

Methods

NMR data were analyzed using “RNAShifts”, a program designed to download and analyze RNA 1H NMR chemical shifts that have been deposited in the BMRB. (Locally derived shifts that have yet to be deposited can also be analyzed). All 131 depositions available in the BMRB were used in the current analysis except BMRB ID 5170, 6814, 4816, 15697, 15915, 5023, 4253, 4894, and 15257, which could not be reliably used because either the BMRB assignments didn’t match the published PDB assignments, or because there was no associated publication or PDB file that could be used to identify RNA secondary structure. As additional input, files were manually generated for each deposition, based on published structural studies, that identify for each residue (1) whether or not the residue is base-paired, (2) the nature of the base-pairing partner, (3) any long-range intra- and/or inter-molecular interactions (e.g., sites of protein binding or participation in A-minor or other RNA–RNA contacts), (4) participation in structured (e.g., GNRA; G/g = guanosine, N/n = any nucleotide; R/r = purine; A/a = adenosine) or unstructured loops. A representative input file is shown in Supplementary Table S1.

The analysis focused on shifts reported for the non-exchangeable H8, H2, H6, H5, H1′, H2′ and H3′ protons of the central base pair of three consecutive canonical Watson–Crick base-pairs (WC-BPs) (here called WC-BP triplets: ([5′−n(i−1)−Ni−n(i+1)]:[5′−n(j−1)−nj−n(j+1)]; Ni = nucleotide for which the NMR shifts are being evaluated; n = neighboring nucleotides), Fig. 1a. As additional parameters, we denoted if the n(i−1):n(j+1) or n(i+1):n(j−1) base pairs were at terminating positions in the RNA, and we identified the secondary structural elements adjacent to the WC-BP triplets (canonical or non-canonical WC-BP, bulges, loops, long-range RNA–RNA interactions, and RNA–protein/ligand interactions), Table S1.

Fig. 1
a Definitions used for base pair triplets. The chemical shifts of the N(i) residue are analyzed in this work, and this strand may be preceded by a base-paired (WC or GU wobble) nucleotide (pre_n) or a non-base paired residue (5loop), or followed by a ...

We chose a relatively conservative approach in modeling the effect of the neighborhood of each central base pair. This was done because there are still, especially in comparison to proteins, relatively few chemical shift assignment sets for RNA deposited at the BMRB. Rather than using any non-linear or neural network approach we used an approach similar to the chemical shift increment method of Pretsch as used in predicting spectra of small organic molecules (Pretsch et al. 2009). Thus, for the central residue of each WC-BP triplet, we defined the attributes describing the neighborhood of the central nucleotide as described above, and calculated the contribution that each attribute makes to the predicted chemical shift. The predicted chemical shift is then a base chemical shift plus the linear contribution of the value corresponding to each attribute present in that nucleotide’s environment. The contribution of each attribute was calculated by linear regression of the chemical shifts in our database of RNA chemical shifts with the set of explanatory variables represented by the neighborhood attributes. The constant term of our regression model corresponds to a nucleotide embedded in a triplet of Watson–Crick base pairs with a U (uridine) flanking it on both the 5′ and 3′ sides and Watson–Crick base-paired nucleotides at the 5′ and 3′ ends of the triplet.

Our analysis included a total of 15 potential variables, Table 1, of which only some might potentially contribute significantly to the shift of a specific atom in a given central nucleotide. Because the approach includes a large number of independent variables relative to the chemical shift datasets, there was a significant danger of over-fitting using a conventional linear regression algorithm. Over-fitting can lead to excellent prediction of the training set, but poor predictive capability on novel datasets. To minimize the risk of over-fitting we chose an algorithm, Pace Regression (Projection Adjustment by Contribution Estimation), that is capable of assessing the importance of each of the parameters. Calculations were performed using the Weka Machine Learning and Data Mining Library system, which allowed us to perform a statistical analysis of the prediction model (Witten et al. 2011). Pace Regression is a linear regression system that uses various information criteria to assess the degree of importance of the regression variables (Wang and Witten 2002). Thus it provides one solution to the subset selection problem: which subset of a set of potential regressors is the appropriate set to explain the data, and thereby minimize the risk of overfitting and maximize the predictive capability on previously unseen data.

Table 1
Sequence variables and chemical shift corrections calculated by Pace regression

Use of Weka provided not only access to Pace Regression, but also various assessments of the quality of the predictions. In particular, we used 10-fold stratified cross-validation during our analysis. Rather than providing correlation coefficients and root mean squared (rms) deviations of the predictions using all the data in the prediction, this technique trains the model on 90 % of the data and then assesses the results of predicting the remaining 10 % of the data. The process is repeated 10 times, using a different subset of the data each time and derives the correlation coefficients and rms deviations based on the whole process. Pace regression was used independently on each atom type present in each of the four central nucleotides for a total of 19 regression calculations.

We were unable in our analysis to adequately identify and control for sample conditions (pH, temperature, ionic strength, etc.) and unusual molecular conformation, and there is a significant possibility of misassignment, especially of some atom types. Therefore, after dropping a single obvious major outlier, we minimized these effects by automatically trimming outliers and automatically adjusting the reference for the chemical shift sets. Automated outlier elimination was performed by running two passes of the Pace Regression for each atom/central nucleotide. In the first pass, the rms deviations between the experimental and predicted values were calculated using all of the data. Any data values that deviated from the predicted values by more than three times the rms deviation value were dropped, and a second pass of the Pace Regression was performed on the now trimmed dataset. Automatic re-referencing was achieved by performing the above analysis (including outlier detection) twice. In the first of these passes, the mean error of prediction was calculated for all the shifts from each BMRB file. Prior to the second pass, each shift was corrected by the mean deviation calculated for the corresponding BMRB file. The chemical shift corrections determined by this approach are listed in Table S3.

The RNAShifts program was written using JTcl (http://jtcl.kenai.com) and Swank (http://swank.kenai.com), which are the Java implementations of the Tcl programming language and Tk graphical user interface toolkit (Ousterhout and Jones 2010). The analysis mode is run in three stages. The first loads BMRB files (fetching them from http://bmrb.wisc.edu if necessary), extracts chemical shifts, and then uses the input template to assign attributes to each shift. The second stage reads the output of the first stage and generates input files in the format used by Weka. The third executes Weka multiple times for each proton type, manages the two passes used for outlier detection and generates various statistical output files. The graphical interface module allows plotting predicted and experimental data subject to various criteria for choosing subsets of the data and attributes for plotting. The RNAShifts program is available upon request from the author (BAJ).

Results and discussion

Outlier chemical shifts

The statistical analysis described above identified 65 chemical shift assignments from the full BMRB database that, after automated re-referencing, deviated from expected values by more than 3 standard deviations. Seven of these assignments were associated with earlier publications from the M.F.S. laboratory, and inspection of the original NMR spectra revealed that these signals had been erroneously assigned (corrections to BMRB files 15113 and 17083 have now been made). We also discovered relatively large systematic chemical shift variations for one of our earlier depositions (BMRB ID 6094) that were associated with improper chemical shift referencing (the residual water signal at 35 °C was erroneously assigned a chemical shift of 4.792 ppm). We therefore updated the BMRB with the modified values, which were used in the present analysis. Based on examination of published NMR spectra, we were able to correct 19 additional assignments in the BMRB—in many cases, the signals had been properly assigned in the published spectra but improperly recorded in the BMRB files. In all cases, the re-assigned (or typo-corrected) shifts were well within the 3-standard deviation cutoff. We were unable to determine the nature of the deviations observed for the remaining 38 outliers because relevant regions of the NMR spectra were not provided in the original publications, and these 38 assignments were not used in subsequent analyses. The majority of these outliers were associated with ribose protons, of which 17 were for highly overlapping H2′ and H3′ proton signals. Thus, of the 3,796 available chemical shifts, 3,758 were retained for analysis and 38 (1 %; mostly ribose assignments) were excluded.

Chemical shifts that were either re-assigned or excluded are summarized in Supplementary Table S2, and referencing corrections employed for all of the utilized depositions are summarized in Supplementary Table S3. The final dataset included values for the central base pairs of all of the 43 possible combinations of WC-BP triplets, with as few as one, and as many as 23, assignments for each of the possible combinations. A total of 137 additional triplets that contain G:U base pairs were also included in the analysis. As shown in Fig. 1b, the retained and re-referenced BMRB shifts (δ) were in good agreement with predicted shifts (δpred) (rms deviation for the entire dataset = 0.056). Good agreement was also obtained when training was performed using a two-fold cross-validation analysis, in which half of the data were used for training and half for validation (rms = 0.069 ppm), and when training was performed with 60 % of randomly-ordered BMRB entries and validation assessed with the remaining 40 % of the data (rms = 0.063, averaged over all atom types).

Chemical shift trends for canonical triplets

The re-referenced NMR chemical shifts (δ) were generally in good agreement with the mean shifts calculated for each unique sequence/atom type (left angle bracketδright angle bracket). For example, excellent correlations were observed in a plot of δ versus left angle bracketδright angle bracket for the central residues of “canonical triplets,” defined here as a triplet that contains only GC and/or AU base pairs and are both preceded and followed by canonical GC or AU base pairs (rms deviation = 0.043), Fig. 1c. The database utilized does not contain chemical shift values for aAa and uCa canonical triplets, nor for the H2′ and/or H3′ protons of the following canonical triplets: aAu (H2′, H3′), uGa (H2′, H3′), aUu (H2′, H3′), gGu (H3′), aCc (H3′). (Note that data were available for non-canonical forms of these triplets and were included in the analysis). There were no significant differences in correlation coefficients obtained upon fitting δ versus left angle bracketδright angle bracket for the A, G, C and U nucleotides, but as observed in plots of δ versus δpred, greater scatter was generally observed for the H2′ and H3′ protons, Fig. 2a.

Fig. 2
Plots of re-referenced 1H NMR chemical shifts (δ) reported for the central adenosine residues within canonical triplets (as defined in text and Fig. 1 caption) versus mean shifts calculated for canonical triplets (left angle bracketδright angle bracket ...

1H NMR chemical shift trends were readily observed in plots that compare δ with the mean shift calculated for canonical triplets (left angle bracketδright angle bracketcan), and with the coefficients obtained with the Pace Regression analysis. Plots of δ versus left angle bracketδright angle bracketcan for the n-A-n canonical triplets are shown in Fig. 2, and data for the n-G-n, n-C-n and n-U-n canonical triplets are plotted in Fig. 3. The contributions of the attributes calculated by Pace Regression are plotted in Fig. 4. The observed trends are consistent with several generalized correlations identified by Wijmenga and co-workers (Cromsigt et al. 2001). For example, δ values for purine-H8 protons in canonical triplets are highly sensitive to the nature of the 5′-residue within the triplet, with 5′-purines associated with more upfield chemical shifts. We further observe that 5′-uridines induce a larger downfield H8 shift than 5′-cytidines (Figs. 2c, c,3a),3a), and that the H8 chemical shift is also sensitive to the nature of the 3′-residue, Figs. 2c and and3a.3a. For example, the A-H8 left angle bracketδright angle bracketcan values observed for n-A-a canonical triplets are consistently downfield relative to those observed for n-A-g canonical triplets, Fig. 2c, and a similar trend is observed for n-G-a versus n-G-g triplets, Fig. 3a.

Fig. 3
Plots of re-referenced 1H NMR chemical shifts (δ) reported for the central guanosine (a), cytosine (b) and uracil (c) residues within canonical triplets (as defined in text and Fig. 1) versus mean shifts calculated for canonical triplets ...
Fig. 4
Plot of the chemical shift contributions (δcontrib) of each attribute relative to a canonical uNu triplet as obtained via Pace Regression for aromatic (a) and ribose (b) proton assignments (positive values denote downfield shifts). Data in these ...

The adenosine-H2 proton is sensitive to the nature of both the 5′- and 3′-nucleotides (Cromsigt et al. 2001) and exhibits a large chemical shift range of ~6.4–8.0 ppm. Importantly, the simultaneous presence of a 5′-pyrimidine and 3′-purine is associated with a significant upfield A-H2 NMR chemical shift, to a less crowded region of the RNA NMR spectrum (6.4–7.1 ppm, Fig. 2b) where they are potentially useful for structural characterization of large RNAs Lu et al. (2011). In contrast, significant downfield shifts are observed for the H2 protons of adenosines that are preceded by a purine and followed by a pyrimidine, Fig. 2b. The H5 protons of the C and U are sensitive to the nature of the preceding residue of the triplet but exhibit almost no detectable sensitivity to the nature of the following residue, Fig. 3c, d. The pyrimidine H6 protons are also more sensitive to the nature of the 5′ residue, but exhibit some sensitivity to the 3′ residue as well (Fig. 3c, e). The ribose protons appear to be sensitive to the nature of both the 5′ and 3′ residues, although the limited chemical shift dispersion and uncertainties regarding some of H2′ and H3′ assignments make it more difficult to identify clear chemical shift trends.

Influence of 5′- and 3′-terminal base pairs within the WC-BP triplet

The presence of 5′- and/or 3′-terminating base pairs within the WC-BP triplet has a significant influence on the chemical shifts of the central residue. As shown in Fig. 5a, the aromatic, H1′, H2′ and H3′ protons of the central residue exhibit small but significant downfield shifts relative to left angle bracketδright angle bracketcan values when adjacent to a 5′-terminating base-paired residue (the single H3′ outlier is most likely due to a misassignment or typo). The most significant perturbations are observed for the aromatic protons, which exhibit deviations in the range of 0.15–0.45 ppm. In contrast, most signals for residues that reside next to a 3′-terminal WC-BP exhibit smaller but nevertheless consistent upfield shifts relative to the left angle bracketδright angle bracketcan values, Fig. 5b. The most significant shifts are observed for H2′ protons which have a mean upfield shift of 0.2 ppm.

Fig. 5
Plots of δ versus left angle bracketδright angle bracketcan (defined in Figs. 1, ,22 captions) for the central residues of WC-BP triplets that contain a 5′-terminal base pair (5ter), a 3′-terminal base pair (3ter), or are ...

Influence of non-canonical elements adjacent to the WCBP triplets

Our analysis assessed the influence of non-canonical structural elements that reside immediately upstream (5loop) or downstream (3loop) of the WC-BP triplets. We defined these elements to include internally stacked residues that are not involved in Watson–Crick base pairing, looped or bulge residues believed to be flexible or structured (e.g., K-turns), and residues involved in base-triples or long-distance RNA–RNA interactions. As shown in Fig. 5d, the presence of non-canonical RNA structures at the 3′-end of the WC-BP triplet does not appear to significantly influence any of the proton shifts associated with the central residue of the triplet. On the other hand, the presence of non-canonical structure on the 5′-side of the WC-BP triplet results in small but significant upfield shifts relative to left angle bracketδright angle bracketcan values for the aromatic and H1′ protons, Fig. 5c.

Influence of G:U base pairing within the triplet

Because GU base pairs are both prevalent and functionally important (Varani and McClain 2000), we also assessed the influence of this class of base pairing on 1H NMR chemical shifts. Systematic variations are apparent for some protons of the central U of triplets when they are base paired with G. Considering only canonical triplets in which the central U:A base pair is substituted by U:G, the H5 protons exhibit a downfield shift and the H1′ and H2′ protons exhibit small upfield shifts, whereas the H6 and H3′ chemical shifts appear to be relatively unperturbed, Fig. 6a. If the central residue of the canonical triplet is a G, base pairing with U results in a small downfield shift of the H2′ NMR signal and upfield shift of H3′ (relative to base pairing with C) but does not significantly affect the shifts of the other G protons, Fig. 6b.

Fig. 6
Plots showing the sensitivity of the 1H NMR chemical shifts to GU and UG wobble pairing within the canonical WC-BP triplet. a The central U of an otherwise canonical triplet is paired with G. b The central G of an otherwise canonical triplet is paired ...

The presence of GU (or UG) base pairs at the n(i−1) or n(i+1) positions can significantly influence the signals of the central residue, and data for otherwise canonical triplets are shown in Fig. 6c–f. For triplets in which the central residues is a pyrimidine, the H1′ and H3′ are relatively unaffected by the presence of a preceding GU wobble, Fig. 6c. However, the H6, H5 and H2′ protons are systematically perturbed, with the u(wob)-U/C-n H6 signal shifted downfield, the g(wob)-U-n H6 signal shifted upfield, and the g(wob)-C-n C-H6 signal shifted downfield relative to the average canonical shifts, Fig. 6c. Interestingly, the u(wob)-C/U-n H5 shifts are relatively unperturbed relative to canonical shifts, whereas g(wob)-C/U-n H5 shifts are generally shifted downfield relative to the signals observed for the canonical triplets, Fig. 6c. Also, H2′ shifts of the central pyrimidine are shifted downfield when preceded by a UG wobble, but are shifted upfield when preceded by a GU wobble, Fig. 6c. When the central residue is a purine, the H1′ and H3′ proton shifts are relatively unaffected by a preceding wobble, but the H8, H2, and H2′ protons generally exhibit systematic downfield shifts, with the magnitude of the shift being somewhat greater for a preceding U(wob) compared to a preceding G(wob), Fig. 6d.

The presence of a subsequent GU wobble can also result in systematic chemical shift perturbations. For triplets in which the central residue is a pyrimidine followed by a U(wob) mismatch, the H6 and H3′ signals exhibits small upfield shifts but the remaining signals do not appear to be significantly perturbed, Fig. 6e. In contrast, the presence of a subsequent G(wob) mismatch does not appear to lead to any detectable perturbations, Fig. 6e. For triplets in which the central residue is a purine, a subsequent G(wob) leads to a small systematic downfield shift of the H1′ proton but does not significantly perturb the other NMR signals, whereas a subsequent U(wob) pair results in small upfield shifts of the H6 and H5 signals and a small downfield shift of the H2′ signal, Fig. 6f.

Chemical shift predictions

The Pace regression approach described above provided predicted chemical shift values for all possible combinations of WC-BP triplet parameters used in the present study, Table 1. 1H NMR chemical shifts observed for the canonical triplets are in good agreement with the shifts predicted using the Pace regression approach described above (δpred), Fig. 7a (rms deviation = 0.050). Excellent agreement was also observed for triplets that contained only a single modifying element (e.g., only a 5ter but no other non-canonical elements), with the greatest deviations observed for a few of the H2′ and H3′ assignments, Fig. 7b–h (rms deviation in the range 0.057–0.057). Good fits were also observed for triplets that contained more than one modifying element (rms deviation for all canonical and non-canonical triplets = 0.056), Fig. 7i. As observed in other fits, the largest deviations are observed for the H2′ and H3′ proton assignments.

Fig. 7
Plots of δ versus predicted chemical shift (δpred), calculated by Pace regression as described in the text. ah Data for triplets that are fully canonical (a) (rmsd = 0.050) or include a single non-canonical element, ...

The data in Table 1 can be used in computer programs such as NMRView (Johnson 2004; Johnson and Blevins 1994) or ad hoc calculations to predict chemical shifts. The constant term represents the value of the given atom in nucleotide i, when the i − 1 and i + 1 nucleotides are both U, and all nucleotides from i − 2 through i + 2 are present and in canonical Watson–Crick base pairs. For example, an A-H2 proton, in a canonical uAu triplet would be at 7.0299. Calculating the shift of the A-H2 proton in a different environment is done by adding to the constant term the contributions from any applicable columns in the A-H2 row of Table 1. For example, the chemical shift of an A-H2 proton in a gAc triplet, in which the i − 2 residue is in a loop, would be: 7.8469 ppm (7.0299 + 0.6899 + 0.0658 + 0.0622). If the i − 1 G is in a GU (rather than GC) base pair, the A-H2 proton chemical shift would be: 7.9217 ppm (7.0299 + 0.7637 + 0.0658 + 0.0622).

Conclusions

The present studies provide the first quantitative analysis of the RNA non-exchangeable 1H NMR chemical shifts in the BMRB. Our studies identify sequence-dependent chemical shift correlations and establish the influence of terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs). Excellent correlations were observed despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A relatively small number of outliers that were not utilized in the analysis, mainly ribose H2′ and H3′ assignments, are likely due to assignment or typographical errors and should be re-examined. Assignments for some triplet combinations were either limited or lacking; for example, the database does not include assignments for two of the 64 possible “canonical triplets.” Although shifts for these triplets could be predicted from assignments made for non-canonical triplets (e.g., WC-BP triplets adjacent to non-canonical structures or that contain terminal or GC base pairs), future studies of oligonucleotides with the missing sequences are clearly in order.

The statistics indicate that the protocol employed for chemical shift predictions, assigning attributes to different triplet environments and then conducting selection and linear model fitting with Pace Regression, performed very well for the data used in this study. However, as we move forward with this research and the number of attributes is expanded, alternative fitting methods such as Neural Networks and allowing attributes to contribute in non-linear modes may be required. The protocol used here can, of course, also be applied to nitrogen and carbon nuclei, and it will be interesting to determine if these nuclei exhibit similar environment- and structure-dependent sensitivities.

The 1H NMR shifts observed for residues that participate directly in long-range RNA–RNA interactions or interactions with ligands or proteins, as identified in the associated publications and/or the structure coordinate (PDB) files, generally deviated from the A-form helical triplet shifts. For example, the H6 and H5 NMR chemical shifts observed for residue U5 of the ScYLV P-1-P2 frameshifting pseudoknot (7.93 and 5.25 ppm, respectively) (Cornish et al. 2005), deviate by 0.24 and 0.29 ppm from the expected values (7.69 and 5.54 ppm, respectively) and are well outside the rms range calculated for canonical gUg triplets (rms = 0.06 and 0.03 ppm, respectively). Significant deviations were also observed for otherwise canonical A-form helical residues that interact with protein elements. In future studies of RNAs with unknown structures, the observation of outlier chemical shifts may serve as useful indicators of potential long-range RNA:RNA or RNA:protein interactions. In addition, the trends identified in the present studies should facilitate the refinement of algorithms used to calculate 1H NMR chemical shifts on the basis of RNA structural coordinates alone (Cromsigt et al. 2001; Case 1995, 2002; Dejaegere et al. 1999; Case et al. 2005), thereby making the 1H NMR chemical shift a more useful parameter for RNA structure refinement. Because the variations in chemical shifts observed for atoms of a given triplet are small, variations in the 3D structures of the triplets should also be small. This observation lends support for refinement approaches that utilize residual dipolar couplings and/or small angle X-ray scattering (SAXS) data to orient idealized A-form helices (Funari et al. 2000; Walsh et al. 2004; Zuo et al. 2008; Grishaev et al. 2008; Wang et al. 2009, 2010; Burke et al. 2012).

In the course of these studies, chemical shift trends were tentatively identified for a number of non-A-form helical structures that are well represented in the BMRB, particularly those of conserved tetraloops (e.g., GNRA). Future studies that include parameterizations for tetraloops, base triples, and other conserved and well-defined RNA sub-structures will likely lead to the identification of additional trends useful for 1H NMR assignment and verification. In addition, it should now be possible to incorporate the approach into software programs to enable semi-automated assignment of RNA, including large RNAs with different combinations of 2H-labeled or segmentally-labeled nucleotides (underway).

Electronic supplementary material

Acknowledgments

Support from the National Institute for General Medical Sciences (GM42561 to M.F.S. and GM103297 to B.J.) is gratefully acknowledged. S.B. was supported by grants from the National Institute of General Medical Sciences for enhancing minority access to research careers (MARC U*STAR 2T34 GM008663) and the Howard Hughes Medical Institute for enhancing undergraduate research training.

Abbreviations

BMRB
Biological Magnetic Resonance Data Bank
NDB
Nucleic Acid Database
PDB
Protein Data Bank
WC-BP
Watson–Crick base pair
G/g
Guanosine
N/n
Any nucleotide
R/r
Purine
Py/py
Pyrimidine
A/a
Adenosine
C/c
Cytosine
U/u
Uridine
left angle bracketδright angle bracket
Mean NMR chemical shift
left angle bracketδright angle bracketcan
Mean NMR chemical shift determined for a canonical triplet, defined here as a stretch of three sequential canonical base pairs that is both preceded and followed by at least one canonical base pair
δpred
Predicted NMR chemical shift

Contributor Information

Bruce A. Johnson, bruce/at/onemoonscientific.com.

Michael F. Summers, summers/at/hhmi.umbc.edu.

References

  • Aeschbacher T, Schubert M, Allain FH. A procedure to validate and correct the 13C chemical shift calibration of RNA datasets. J Biomol NMR. 2012;52:179–190. doi: 10.1007/s10858-011-9600-7. [PubMed] [Cross Ref]
  • Allain FH-T, Varani G. How accurately and precisely can RNA structure be determined by NMR? J Mol Biol. 1997;267:338–351. doi: 10.1006/jmbi.1996.0855. [PubMed] [Cross Ref]
  • Bahrami A, Clos LJ, 2nd, Markley JL, Butcher SE, Eghbalnia HR. RNA-PAIRS: RNA probabilistic assignment of imino resonance shifts. J Biomol NMR. 2012;52:289–302. doi: 10.1007/s10858-012-9603-z. [PMC free article] [PubMed] [Cross Ref]
  • Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–297. doi: 10.1016/S0092-8674(04)00045-5. [PubMed] [Cross Ref]
  • Batey RT, Inada M, Kujawinski E, Puglisi JD, Williamson JR. Preparation of isotopically labeled ribonucleotides for multidimensional NMR spectroscopy of RNA. Nucleic Acids Res. 1992;20:4515–4523. doi: 10.1093/nar/20.17.4515. [PMC free article] [PubMed] [Cross Ref]
  • Batey RT, Battiste JL, Williamson JR. Preparation of isotopically enriched RNAs for heteronuclear NMR. Methods Enzymol. 1995;261:300–322. doi: 10.1016/S0076-6879(95)61015-4. [PubMed] [Cross Ref]
  • Bessonov S, Anokhina M, Will CL, Urlaub H, Luhrmann R. Isolation of an active step I spliceosome and composition of its RNP core Nature 2008. 452846–850.8502008Natur.452..846B. doi: 10.1038/nature06842. [PubMed] [Cross Ref]
  • Boisvert FM, van Koningsbruggen S, Navascues J, Lamond AI. The multifunctional nucleolus. Nat Rev Mol Cell Biol. 2007;8:574–585. doi: 10.1038/nrm2184. [PubMed] [Cross Ref]
  • Brodersen P, Voinnet O. The diversity of RNA silencing pathways in plants. Trends Genet. 2006;22:268–280. doi: 10.1016/j.tig.2006.03.003. [PubMed] [Cross Ref]
  • Burke JE, Sashital DG, Zuo X, Wang YX, Butcher SE. Structure of the yeast U2/U6 snRNA complex. RNA. 2012;18:673–683. doi: 10.1261/rna.031138.111. [PubMed] [Cross Ref]
  • Case DA. Calibration of ring-current effects in proteins and nucleic acids. J Biomol NMR. 1995;6:341–346. doi: 10.1007/BF00197633. [PubMed] [Cross Ref]
  • Case DA (2002) AMBER 7 users manual. University of California San Francisco.
  • Case DA, Cheatham TEI, Darden T, Gohlke H, Luo R, Merz KMJ, Onufriev A, Simmerling C, Wang B, Woods R. The Amber biomolecular simulation programs. J Computat Chem. 2005;26:1668–1688. doi: 10.1002/jcc.20290. [PMC free article] [PubMed] [Cross Ref]
  • Cavalli A, Salvatella X, Dobson CM, Vendruscolo M. Protein structure determination from NMR chemical shifts Proc Natl Acad Sci USA 2007. 1049615–9620.96202007PNAS..104.9615C. doi: 10.1073/pnas.0610313104. [PubMed] [Cross Ref]
  • Cornish PV, Hennig M, Giedroc DP. A loop 2 cytidine-stem 1 minor groove interaction as a positive determinant for pseudoknot-stimulated -1 ribosomal frameshifting Proc Natl Acad Sci USA 2005. 10212694–12699.126992005PNAS..10212694C. doi: 10.1073/pnas.0506166102. [PubMed] [Cross Ref]
  • Cromsigt JAMTC, Hilbers CW, Wijmenga SS. Prediction of proton chemical shifts in RNA. Their use in structure refinement and validation. J Biomol NMR. 2001;21:11–29. doi: 10.1023/A:1011914132531. [PubMed] [Cross Ref]
  • Davis JH, Tonelli M, Scott LG, Jaeger L, Williamson JR, Butcher SE. RNA helical packing in solution: NMR structure of a 30 kDa GAAA tetraloop-receptor complex. J Mol Biol. 2005;351:371–382. doi: 10.1016/j.jmb.2005.05.069. [PubMed] [Cross Ref]
  • Dejaegere A, Bryce RA, Case DA (1999) An empirical analysis of proton chemical shifts in nucleic acids, In: Facelli J et al (eds) ACS symposium series: modeling NMR chemical shifts. American Chemical Society, Washington, DC, pp 194–206.
  • Doudna JA, Rath VL. Structure and function of the eukaryotic ribosome: the next frontier. Cell. 2002;109:153–156. doi: 10.1016/S0092-8674(02)00725-0. [PubMed] [Cross Ref]
  • D’Souza V, Summers MF. Structural basis for packaging the dimeric genome of Moloney Murine Leukaemia Virus Nature 2004. 431586–590.5902004Natur.431..586D. doi: 10.1038/nature02944. [PubMed] [Cross Ref]
  • D’Souza V, Dey A, Habib D, Summers MF. NMR structure of the 101 nucleotide core encapsidation signal of the Moloney Murine Leukemia Virus. J Mol Biol. 2004;337:427–442. doi: 10.1016/j.jmb.2004.01.037. [PubMed] [Cross Ref]
  • Duss O, Lukavsky PJ, Allain FH. Isotope labeling and segmental labeling of larger RNAs for NMR structural studies. Adv Exp Med Biol. 2012;992:121–144. doi: 10.1007/978-94-007-4954-2_7. [PubMed] [Cross Ref]
  • Ebrahimi M, Rossi P, Rogers C, Harbison GS. Dependence of 13C NMR chemical shifts on conformations of RNA nucleosides and nucleotides J Magn Reson 2001. 1501–9.92001JMagR.150....1E. doi: 10.1006/jmre.2001.2314. [PubMed] [Cross Ref]
  • Edwards TE, Klein DJ, Ferre-D’Amare AR. Riboswitches: small-molecule recognition by gene regulatory RNAs. Curr Opin Struct Biol. 2007;17:273–279. doi: 10.1016/j.sbi.2007.05.004. [PubMed] [Cross Ref]
  • Ehresmann C, Baudin F, Mougel M, Romby P, Ebel JP, Ehresmann B. Probing the structure of RNAs in solution. Nucleic Acids Res. 1987;15:9109–9128. doi: 10.1093/nar/15.22.9109. [PMC free article] [PubMed] [Cross Ref]
  • Fares C, Amata I, Carlomagno T. 13C-detection in RNA bases: revealing structure-chemical shift relationships. J Am Chem Soc. 2007;129:15814–15823. doi: 10.1021/ja0727417. [PubMed] [Cross Ref]
  • Forconi M, Herschlag D. Metal ion-based RNA cleavage as a structural probe. Methods Enzymol. 2009;468:91–106. doi: 10.1016/S0076-6879(09)68005-8. [PubMed] [Cross Ref]
  • Funari SS, Rapp G, Perbandt M, Dierks K, Vallazza M, Betzel C, Erdmann VA, Svergun DI. Structure of free Thermus flavus 5 S rRNA at 1.3 nm resolution from synchrotron X-ray solution scattering. J Biol Chem. 2000;275:31283–31288. doi: 10.1074/jbc.M004974200. [PubMed] [Cross Ref]
  • Grishaev A, Ying J, Canny MD, Pardi A, Bax A. Solution structure of tRNAVal from refinement of homology model against residual dipolar coupling and SAXS data. J Biomol NMR. 2008;42:99–109. doi: 10.1007/s10858-008-9267-x. [PMC free article] [PubMed] [Cross Ref]
  • Grzesiek S, Bax A. Amino acid type determination in the sequential assignment procedure of uniformly 13C/15N-enriched proteins. J Biomol NMR. 1993;3:185–204. [PubMed]
  • Hassouna N, Michot B, Bachellerie JP. The complete nucleotide sequence of mouse 28S rRNA gene. Implications for the process of size increase of the large subunit rRNA in higher eukaryotes. Nucleic Acids Res. 1984;12:3563–3583. doi: 10.1093/nar/12.8.3563. [PMC free article] [PubMed] [Cross Ref]
  • Heng X, Kharytonchyk S, Garcia EL, Lu K, Sachin Divakaruni S, LaCotti C, Edme K, Telesnitsky A, Summers MF. Identification of a minimal HIV-1 RNA packaging signal. J Mol Biol. 2012;417:224–239. doi: 10.1016/j.jmb.2012.01.033. [PMC free article] [PubMed] [Cross Ref]
  • Houck-Loomis B, Durney MA, Salguero C, Shankar N, Nagle JM, Goff SP, D’Souza VM. An equilibrium-dependent retroviral mRNA switch regulates translational recoding Nature 2011. 480561–564.5642011Natur.480..561H. [PMC free article] [PubMed]
  • Johnson BA. Using NMRView to visualize and analyze the NMR spectra of macromolecules. Methods Mol Biol. 2004;278:313–352. [PubMed]
  • Johnson BA, Blevins RA. NMRview: a computer program for the visualization and analysis of NMR data. J Biomol NMR. 1994;4:603–614. doi: 10.1007/BF00404272. [PubMed] [Cross Ref]
  • Kim VN. Small RNAs: classification, biogenesis, and function. Mol Cells. 2005;19:1–15. doi: 10.1016/j.molcel.2005.05.026. [PubMed] [Cross Ref]
  • Kim I, Watanabe S, Muto Y, Hosono K, Takai K, Takaku H, Kawai G, Watanabe K, Yokoyama S. Selective deuteration of RNA for NMR signal assignment. Nucleic Acids Symp Ser. 1995;34:123–124. [PubMed]
  • Kim I, Lukavsky PJ, Puglisi JD. NMR study of 100 kDa HCV IRES RNA using segmental isotope labeling. J Am Chem Soc. 2002;124:9338–9339. doi: 10.1021/ja026647w. [PubMed] [Cross Ref]
  • Kladwang W, Vanlang CC, Cordero P, Das R. Understanding the errors of SHAPE-directed RNA structure modeling. Biochemistry. 2011;50:8049–8056. doi: 10.1021/bi200524n. [PMC free article] [PubMed] [Cross Ref]
  • Korostelev A, Noller HF. The ribosome in focus: new structures bring new insights. Trends Biochem Sci. 2007;32:434–441. doi: 10.1016/j.tibs.2007.08.002. [PubMed] [Cross Ref]
  • Lam SL, Chi LM. Use of chemical shifts for structural studies of nucleic acids. Prog Nucl Magn Reson Spectrosc. 2010;56:289–310. doi: 10.1016/j.pnmrs.2010.01.002. [PubMed] [Cross Ref]
  • Lu K, Miyazaki Y, Summers MF. Isotope labeling strategies for NMR studies of RNA. J Biomol NMR. 2009;46:113–125. doi: 10.1007/s10858-009-9375-2. [PMC free article] [PubMed] [Cross Ref]
  • Lu K, Heng X, Garyu L, Monti S, Garcia E, Kharytonchyk S, Dorjsuren B, Kulandaivel G, Jones S, Hiremath A, Sachin Divakaruni S, LaCotti C, Barton S, Tummillo D, Hosic A, Edme K, Albrecht S, Telesnitsky A, Summers MF. NMR detection of structures in the HIV-1 5′-leader RNA that regulate genome packaging Science 2011. 344242–245.2452011Sci...334..242L. doi: 10.1126/science.1210460. [PMC free article] [PubMed] [Cross Ref]
  • Lu K, Heng X, Summers MF. Structural determinants and mechanism of HIV-1 genome packaging. J Mol Biol. 2011;410:609–633. doi: 10.1016/j.jmb.2011.04.029. [PMC free article] [PubMed] [Cross Ref]
  • Lukavsky PJ, Puglisi JD. Structure determination of large biological RNAs. Methods Enzymol. 2005;394:399–415. doi: 10.1016/S0076-6879(05)94016-0. [PubMed] [Cross Ref]
  • Lukavsky PJ, Kim I, Otto GA, Puglisi JD. Structure of HCV IRES domain II determined by NMR. Nat Struct Biol. 2003;10:1033–1038. doi: 10.1038/nsb1004. [PubMed] [Cross Ref]
  • Michnicka MJ, Harper JW, King GC. Selective isotopic enrichment of synthetic RNA: application to the HIV-1 TAR element. Biochemistry. 1993;32:395–400. doi: 10.1021/bi00053a002. [PubMed] [Cross Ref]
  • Miyazaki Y, Irobalieva RN, Tolbert BS, Smalls-Manty A, Iyalla K, Loeliger K, D’Souza V, Khant H, Schmid MF, Garcia E, Telesnitsky A, Chiu W, Summers MF. Structure of a conserved retroviral RNA packaging element by NMR spectroscopy and cryo-electron tomography. J Mol Biol. 2010;404:751–772. doi: 10.1016/j.jmb.2010.09.009. [PMC free article] [PubMed] [Cross Ref]
  • Nelissen FHT, van Gammeren AJ, Tessari M, Girard FC, Heus HA, Wijmenga SS. Multiple segmental and selective isotope labeling of large RNA for NMR structural studies. Nucl Acids Res. 2008;36:e89. doi: 10.1093/nar/gkn397. [PMC free article] [PubMed] [Cross Ref]
  • Nikonowicz EP, Pardi A. Three-dimensional heteronuclear NMR studies of RNA Nature 1992. 355184–186.1861992Natur.355..184N. doi: 10.1038/355184a0. [PubMed] [Cross Ref]
  • Nikonowicz EP, Sirr A, Legault P, Jucker FM, Baer LM, Pardi A. Preparation of 13C and 15N labelled RNAs for heteronuclear multi-dimensional NMR studies. Nucleic Acids Res. 1992;20:4507–4513. doi: 10.1093/nar/20.17.4507. [PMC free article] [PubMed] [Cross Ref]
  • Ohlenschlager O, Haumann S, Ramachandran R, Gorlach M. Conformational signatures of 13C chemical shifts in RNA ribose. J Biomol NMR. 2008;42:139–142. doi: 10.1007/s10858-008-9271-1. [PubMed] [Cross Ref]
  • Ousterhout JK, Jones K. Tcl and the Tk Toolkit. 2. Upper Saddle River, NJ: Addison-Wesley; 2010.
  • Peattie DA, Gilbert W. Chemical probes for higher-order structure in RNA Proc Natl Acad Sci USA 1980. 774679–4682.46821980PNAS...77.4679P. doi: 10.1073/pnas.77.8.4679. [PubMed] [Cross Ref]
  • Peterson RD, Theimer CA, Wu H, Feigon J. New applications of 2D filtered/edited NOESY for assignment and structure elucidation of RNA and RNA-protein complexes. J Biomol NMR. 2004;28:59–67. doi: 10.1023/B:JNMR.0000012861.95939.05. [PubMed] [Cross Ref]
  • Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136:629–641. doi: 10.1016/j.cell.2009.02.006. [PubMed] [Cross Ref]
  • Pretsch E, Bühlmann P, Badertscher M. Structure determination of organic compounds: tables of spectral data. 4. Berlin: Springer; 2009.
  • Shen Y, Lange O, Delaglio F, Rossi P, Aramini JM, Liu G, Eletsky A, Wu Y, Singarapu KK, Lemak A, Ignatchenko A, Arrowsmith CH, Szyperski T, Montelione GT, Baker D, Bax A. Consistent blind protein structure generation from NMR chemical shift data Proc Natl Acad Sci USA 2008. 1054685–4690.46902008PNAS..105.4685S. doi: 10.1073/pnas.0800256105. [PubMed] [Cross Ref]
  • Steitz TA. A structural understanding of the dynamic ribosome machine. Nat Rev Mol Cell Biol. 2008;9:242–253. doi: 10.1038/nrm2352. [PubMed] [Cross Ref]
  • Stern S, Moazed D, Noller HF. Structural analysis of RNA using chemical and enzymatic probing monitored by primer extension. Methods Enzymol. 1988;164:481–489. doi: 10.1016/S0076-6879(88)64064-X. [PubMed] [Cross Ref]
  • Tolbert BS, Miyazaki Y, Barton S, Kinde B, Starck P, Singh R, Bax A, Case DA, Summers MF. Major groove width variations in RNA structures determined by NMR and impact of 13C residual chemical shift anisotropy and 1H–13C residual dipolar coupling on refinement. J Biomol NMR. 2010;47:205–219. doi: 10.1007/s10858-010-9424-x. [PMC free article] [PubMed] [Cross Ref]
  • Varani G, McClain WH. The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems. EMBO Rep. 2000;1:18–23. doi: 10.1093/embo-reports/kvd001. [PubMed] [Cross Ref]
  • Wakeman CA, Winkler WC, Dann CE., III Structural features of metabolite-sensing riboswitches. Trends Biochem Sci. 2007;32:415–424. doi: 10.1016/j.tibs.2007.08.005. [PMC free article] [PubMed] [Cross Ref]
  • Walsh JD, Cabello-Villegas J, Wang YX. Periodicity in residual dipolar couplings and nucleic acid structures. J Am Chem Soc. 2004;126:1938–1939. doi: 10.1021/ja039446g. [PubMed] [Cross Ref]
  • Wang Y, Witten IH (2002) Modeling for optimal probability prediction. In: Proceedings of the 19th International Conference in Machine Learning, Sydney, Australia, pp 650–657.
  • Wang J, Zuo X, Yu P, Xu H, Starich MR, Tiede DM, Shapiro BA, Schwieters CD, Wang YX. A method for helical RNA global structure determination in solution using small-angle x-ray scattering and NMR measurements. J Mol Biol. 2009;393:717–734. doi: 10.1016/j.jmb.2009.08.001. [PMC free article] [PubMed] [Cross Ref]
  • Wang YX, Zuo X, Wang J, Yu P, Butcher SE (2010) Rapid global structure determination of large RNA and RNA complexes using NMR and small-angle X-ray scattering. Methods 52:180–191. [PubMed]
  • Weeks KM. Advances in RNA structure analysis by chemical probing. Curr Opin Struct Biol. 2010;20:295–304. doi: 10.1016/j.sbi.2010.04.001. [PMC free article] [PubMed] [Cross Ref]
  • Wishart DS, Sykes BD. The 13C chemical-shift index: a simple method for the identification of protein secondary structure using 13C chemical-shift data. J Biomol NMR. 1994;4:171–180. doi: 10.1007/BF00175245. [PubMed] [Cross Ref]
  • Wishart DS, Sykes BD, Richards FM. Relationship between nuclear magnetic resonance chemical shift and protein secondary structure. J Mol Biol. 1991;222:311–333. doi: 10.1016/0022-2836(91)90214-Q. [PubMed] [Cross Ref]
  • Wishart DS, Sykes BD, Richards FM. The chemical shift index: a fast and simple method for the assignment of protein secondary structure through NMR spectroscopy. Biochemistry. 1992;31:1647–1651. doi: 10.1021/bi00121a010. [PubMed] [Cross Ref]
  • Wishart DS, Arndt D, Berjanskii M, Tang P, Zhou J, Lin G. CS23D: a web server for rapid protein structure generation using NMR chemical shifts and sequence data. Nucleic Acids Res. 2008;36:W496–W502. doi: 10.1093/nar/gkn305. [PMC free article] [PubMed] [Cross Ref]
  • Witten IH, Frank E, Hall MA. Data mining: practical machine learning tools and techniques. 3. USA: Morgan Kaupmann; 2011.
  • Wüthrich K. NMR of proteins and nucleic acids. New York: Wiley; 1986.
  • Xu J, Lapham J, Crothers DM. Determining RNA solution structure by segmental isotopic labeling and NMR: application to Caenorhabditis elegans spliced leader RNA 1 Proc Natl Acad Sci USA 1996. 9344–48.481996PNAS...93...44X. doi: 10.1073/pnas.93.1.44. [PubMed] [Cross Ref]
  • Zhou J, Bean RL, Vogt VM, Summers MF. Solution structure of the Rous sarcoma virus nucleocapsid protein:uY RNA packaging signal complex. J Mol Biol. 2006;365:453–467. doi: 10.1016/j.jmb.2006.10.013. [PMC free article] [PubMed] [Cross Ref]
  • Zuo X, Wang J, Foster TR, Schwieters CD, Tiede DM, Butcher SE, Wang Y-X. Global molecular structure and interfaces: refining an RNA:RNA complex structure using solution X-ray scattering data. J Am Chem Soc. 2008;130:3292–3293. doi: 10.1021/ja7114508. [PubMed] [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer