- Research article
- Open Access
Sequence dependent variations in RNA duplex are related to non-canonical hydrogen bond interactions in dinucleotide steps
BMC Research Notes volume 7, Article number: 83 (2014)
Sequence determines the three-dimensional structure of RNAs, and thereby plays an important role in carrying out various biological functions. RNA duplexes containing Watson-Crick (WC) basepairs, interspersed with non-Watson-Crick basepairs, are the dominant structural unit and form the scaffold for the 3-dimensional structure of RNA. It is therefore crucial to understand the geometric variation in the dinucleotide steps that form the helices. We have carried out a detailed analysis of the dinucleotide steps formed by AU and GC Watson-Crick basepairs in RNA structures (both free and protein bound) and compared the results to that seen in DNA. Further, the effect of protein binding on these steps was examined by comparing steps in free RNA structures with protein bound RNA structures.
Characteristic sequence dependent geometries are observed for the RR, RY and YR type of dinucleotide steps in RNA. Their geometric parameters show correlated variations that are different from those observed in B-DNA helices. Subtle, but statistically significant differences are seen in roll, slide and average propeller-twist values, between the dinucleotide steps of free RNA and protein bound RNA structures. Many non-canonical cross-strand and intra-strand hydrogen bonds were identified that can stabilise the RNA dinucleotide steps, among which YR steps show presence of many new unreported interactions.
Our work provides for the first time a detailed analysis of the conformational preferences exhibited by Watson-Crick basepair containing steps in RNA double helices. Overall, the WC dinucleotide steps show considerable conformational variability. Furthermore, we have identified hydrogen bond interactions in several of the dinucleotide steps that could play a role in determining the preferred geometry, in addition to the intra-basepair hydrogen bonds and stacking interactions. Protein binding affects the conformation of the steps that are in direct contact, as well as allosterically affect the steps that are not in direct physical contact.
The double helical structure of nucleic acids exists in various polymorphic sub-states, of which DNA prefers B-form conformation and RNA prefers A-form. The A-form in RNA is a right-handed helix formed by stacking of Watson-Crick (WC) basepairs along with a few non-Watson-Crick (NWC) basepairs. The overall conformation of the helices is dictated by the geometry of successive dinucleotide steps, which in-turn is dictated by the chemical nature of bases involved in forming these step. A-form helix is characterized by large roll angle and negative slide values, compared to B-form and is accompanied by a narrow but deep major groove and a wide, shallow minor groove. Significant progress has been made in the understanding of sequence dependent conformational preference of dinucleotide steps in DNA[2–6], while the geometric preference in RNA helices remains largely unexplored. Information on the geometric preference at the step level, can contribute to the overall understanding of the structural organization of RNA. A database of all possible dinucleotide steps with their step parameter values is available in the public domain, but there is very little explicit discussion about the conformational features of various dinucleotide steps observed in RNA.
The helical regions in DNA are comprised almost exclusively of canonical Watson-Crick (AT and GC) basepairs, which form 10 unique dinucleotide steps. On the contrary, an RNA duplex consists of canonical WC basepairs (AU and GC), that are paired along Watson-Crick edge in cis orientation (cisWW family), as well as basepairs involving the Hoogsteen and Sugar edges. In addition, other non-canonical basepairs from the cisWW family e.g. GU, GA, UU also occur frequently in the RNA duplex[8–11]. Thus, the type of dinucleotide steps that can occur in RNA is much larger and complicates any analysis of their sequence dependent conformations. In the present work, we focus only on the dinucleotide steps formed by the two canonical WC basepairs (AU and GC) of the cisWW family, which constitute a major proportion of the helical steps and compare them to the corresponding A-like and B-like DNA steps.
The diverse RNA structural motifs and equally diverse proteins that interact with RNA suggest that the conformational changes that occur during RNA-protein interaction can be characterized by changes in the protein, the RNA, or both. An analysis of crystal structures of free and RNA bound proteins, suggested that the proteins do not show any significant change on binding to RNA. A number of studies on RNA-protein interaction have focused on the RNA-protein interface and recognition mechanism[14–16], but none have analyzed the protein induced conformational change, if any, in RNA. On the other hand, the effect of protein binding on dinucleotide steps of DNA is well documented[4, 17–19]. Specific RNA-binding domains recognize the sequence and shape of the interacting region, while others interact in a non-specific manner. The wide and shallow minor groove present in RNA helix allows easy access for interaction with protein. However, sometimes, the narrow and deep major groove can also interact with proteins owing to presence of mismatch basepairs and bulges along the helix that lead to widening of the major groove. Understanding the conformational changes induced by the interaction of protein, on the dinucleotide steps of the RNA helix, can help deduce general mechanism involved in RNA-protein recognition.
The geometry of dinucleotide step is mainly influenced by basepair hydrogen bonds and stacking interactions[22–24]. Apart from the standard hydrogen bonds involved in basepairing, additional interactions involving base, ribose sugar (especially O2′ group) and phosphate atoms in the RNA backbone were reported in RNA structures[25–27]. These interactions mainly occur in hairpin, internal or junction loops or as part of tertiary interactions. In the light of the developments in our understanding of ‘weak hydrogen bonds’[28, 29], the importance of such interactions in the structure formation, folding and stability of various macromolecules are also being investigated[30–33]. Presence of potentially weak cross-strand and intra-strand hydrogen bond interactions in dinucleotide steps of B-DNA crystal structures have been reported and analyzed[30, 34, 35]. In RNA crystal structures, such interactions between bases in a dinucleotide step have not been reported. Moreover, the characteristic A-like geometry seen in RNA helices can possibly prevent their formation or favor some novel interactions. Recent molecular dynamic (MD) simulation studies on modeled RNA duplexes suggest that potential interactions in a dinucleotide steps are present between exo-cyclic atoms. Hence, identification of these additional stabilizing forces in various dinucleotide steps in RNA crystal structures can help understand their preferred geometry and in ab initio modeling of RNA structures.
In this work, a non-redundant RNA crystal structure dataset has been created and we have examined the intrinsic geometries of all WC basepairs and the dinucleotide steps formed by them in the helical regions. To understand the extent of conformational variability that can occur in the dinucleotide steps, the effect of protein binding on these steps was examined by comparing the helices in free RNA with protein bound RNA dataset. Further, we have carried out a systematic analysis of dinucleotide steps to identify potential hydrogen bond interactions between all four bases and correlated the occurrence of such bonds with the dinucleotide step geometry.
Preparation of dinucleotide dataset
The x-ray crystal structure dataset was created by extracting structures with resolution better than 3.0 Å from the Protein Data Bank. The dataset was made non-redundant using the web servers HD-RNAS and FR3D. Non-standard bases and other chemically modified bases were not included in this study. RNA structures that are not bound to proteins were grouped as ‘free-RNA’ dataset and those that were in complex with a protein were grouped as ‘bound-RNA’ dataset (Table 1). A non-redundant free DNA dataset containing structures with resolution better than 2.0 Å was created to compare it with RNA datasets. The 10 dinucleotide steps in DNA helices were grouped into A-like and B-like steps, based on their Zp value. They are referred to as ‘ADNA’ and ‘BDNA’ dataset respectively. The parameters obtained from crystal datasets are also compared with those of fibre diffraction models. A standard B-DNA fibre model (fibre-BDNA) was generated using NUCGEN. Unlike RNA crystal structures, the basepairs in RNA fibre models in the literature have small negative propeller-twist (-2.1°). 3DNA v2.1 has the option to generate a uniform A-RNA double helix with large negative propeller-twist (-10.5°)[41, 42]. We have used this RNA model (referred henceforth as ‘ModelRNA’) in our analysis for a more realistic comparison with crystal RNA datasets.
Intra-basepair and dinucleotide step parameters
The helical structures in each of the datasets were subjected to geometry based identification and classification of basepairs using BPFind program with default criteria. Helical stems containing 4 or more basepairs alone were included in the datasets. We analyzed steps formed by AU and GC basepair combinations of the cisWW family. The dinucleotide steps were grouped based on their sequence. The six intra-basepair and six dinucleotide step parameters were calculated using NUPARM program[40, 44]. All the intra-basepair and dinucleotide step parameters were calculated using the default option of line joining C6 - C8 atoms as y-axis. Additionally, Zp and cup parameters were also calculated for each step. ‘Zp’ relates the basepairs of the step to their backbone geometries. In a dinucleotide step, it gives the displacement of phosphate atoms in each strand from the midplane between the two stacked basepairs and is the best discriminator between A-form and B-form conformation[4, 19]. ‘Cup’ is the difference in buckle parameter between the two basepairs that form the dinucleotide step. Correlation between parameters was analysed and their statistically significant difference (P < 0.01) were checked using Pearson- Correlation coefficient (r) value. The distribution of the data was shown using a Mahalanobis ellipse fitted on to datapoints with the mean as centre and to cover 90% of the datapoints in each group. The stacking area overlap between the basepairs forming a dinucleotide step was calculated using 3DNA program. Stacking area overlap was calculated between bases by including all atoms, as well as the ring atoms alone (excluding exo-cyclic atoms). The total stacking area overlap is the sum of two intra-strand and two cross-strand overlap between the 4 bases involved in the step.
Hydrogen bond analysis
Hydrogen atoms coordinates was added to all the crystal structures using REDUCE program. The following criteria were used to identify hydrogen bonds (i) donor to acceptor distance (D..A) ≤ 3.8 Å, (ii) Angle D - H..A ≥ 90°. Only bonds that are observed in more than 50% of the cases in each of the 10 dinucleotide steps are discussed.
Dinucleotide steps interacting with protein
In the bound-RNA dataset, the interactions between the atoms in a dinucleotide steps and the protein atoms were identified using CONTACT program in CCP4 program suite. A contact distance of ≤ 4 Å between any pair of amino acid and RNA atom was considered to be interacting. Thus, the steps in the bound RNA dataset were sub-classified into two datasets, those that are in contact with protein (cont) and those that are not in contact (non-cont). In order to assess if the difference in step parameter values between the various datasets was significant, an unpaired student-t-test was carried out.
MATLAB was used for all statistical analysis and for plotting graphs.
The free-RNA dataset consists of 88 protein-free x-ray crystal structures, while bound-RNA dataset includes 127 structures (Table 1). Canonical WC basepairs (AU and GC) constitute more than 83% of the total basepairs in these structures and ~ 74% of dinucleotide steps are comprised of these basepairs (Table 2). This work focuses on the sequence dependent conformational preferences of the 10 dinucleotide steps formed by WC basepairs and are compared with those observed in A-DNA and B-DNA helices.
Dinucleotide step geometries
The intra-basepair parameters of WC basepairs present in the RNA datasets are comparable to that of the model structure but are characterized by large variations. The only noticeable features are that GC basepairs have higher negative buckle compared to AU basepairs, while AU basepairs show higher open angle value compared to GC (Table 3). Overall, AU and AT basepairs have lower buckle, slightly larger negative propeller-twist and open angle values compared to GC basepairs, in both RNA and ADNA crystal structures. RNA helices are GC rich but the datasets contain a good representation of all 10 unique dinucleotide steps. These dinucleotide steps are generally sub-grouped into three broad categories; Purine-Purine (RR)/ Pyrimidine-Pyrimidine (YY), Purine-Pyrimidine (RY) and Pyrimidine-Purine (YR). The mean and standard deviation values of the step parameters, along with the corresponding average propeller-twist, cup and Zp values, for the three types of dinucleotide steps in the free-RNA and bound-RNA dataset are tabulated in Tables 4,5 and6. The parameter values for the ModelRNA are also listed. A comparison of the crystal structure geometries with the model structure values indicates that these steps show some characteristic sequence dependent preferences in their geometries.
In general, average propeller-twist value is higher for dinucleotide steps containing only AU basepairs (AA/UU, AU/AU and UA/UA) and lower for steps with only GC basepairs (GG/CC, GC/GC and CG/ CG). Roll value differs between RR, RY and YR steps (YR > RR > RY), though the overall average roll value for WC basepair containing steps is lower than ModelRNA value (Table 4). It is interesting to note that among the RR sequences, the GG/CC step show slightly larger negative slide and positive Zp value, while AA/UU step has the smallest negative slide value, which is also reflected in their smaller Zp value. All other parameters have similar values.
Among RY steps, AC/GU and GC/GC steps have smallest positive roll angles and negative slide values, among all dinucleotide steps (Table 5). However, AU/AU steps have high roll and large negative average propeller-twist and low slide values, with correspondingly small Zp value, as compared to AC/GU and GC/GC steps. This finding was specific to RNA since the equivalent AT/AT steps in both ADNA and BDNA did not show any such difference when compared to other RY steps. All three YR steps have larger roll and negative slide values when compared to the RR and RY steps, as well as the ModelRNA (Table 6). The UA/UA steps show particularly high mean roll angle (14.1°) compared to CA/UG (11.0°) and CG/CG (11.5°). In addition, the slide and cup values for CG/CG steps have larger negative values, than those for CA/UG and UA/UA steps. Most of these sequence dependent features are unique to RNA helices, however the trends observed for them, seem to be similar to those reported earlier for DNA, particularly the trend observed for roll angle values (YR > RR > RY). Hence, an analysis of these trends and correlation between various parameters in freeRNA, ADNA and BDNA datasets has been carried out to identify any specific structure based features.
Correlation between dinucleotide step parameters
We have carried out a pair-wise correlation analysis between the various step parameters in the free-RNA dataset and compared the correlation coefficient values (r) with those of ADNA and BDNA datasets in order to identify correlations that are specific to A-form helices. The parameters that show statistically significant correlations at confidence level > 99.9% being discussed further (Figure 1, Additional file1). The well-characterized strong correlation between roll and twist that is observed for BDNA is not seen for either ADNA or free-RNA dataset. On the other hand, correlations between shift and tilt and slide and twist are present in both RNA and DNA datasets. Interestingly, the major differences between A and B-form structures are seen for correlations between basepair geometry dependent step parameters and other dinucleotide step parameters. For instance, average-propeller twist is positively correlated with roll, slide and rise in BDNA. However, in ADNA and RNA datasets it is negatively correlated with roll, twist and slide. Twist shows a significant negative correlation with cup value in BDNA, which is absent in A-form structures. Instead, roll and rise show negative correlations with cup in A-form structures, while slide shows a positive correlation with cup. Thus, overall, dinucleotide step parameters in RNA and ADNA datasets show similar correlations that are distinct from those in BDNA dataset.
We have also carried out a correlation analysis, for free-RNA dataset, considering the dinucleotide steps in each of the three sub-groups, RR, RY and YR separately (Figure 2). No significant correlation is seen between roll and twist, for any of the three sub-groups, confirming that this is absent in RNA. Several correlations between step parameter values showed same trend for RR, RY and YR type of steps, e.g. shift with tilt and slide with twist. Similarly, average propeller-twist shows significant negative correlation with slide and roll for all three dinucleotide step types. Cup also shows a negative correlation with roll in all three sub-groups. Thus, a comparison of correlations between parameters of the three step types suggests that the correlations seen in majority of the step are similar to those seen in the pooled dataset and are characteristic of A-form structure.
Effect of Protein binding on the dinucleotide step geometry
The mean and standard deviation values of intra-basepair parameters for canonical AU and GC basepairs in the bound-RNA dataset were compared with that of free-RNA dataset (Table 3). Almost all the parameters have larger standard deviation values in the bound dataset, indicating larger conformational sampling. However, mean values of intra-basepair parameters in bound-RNA dataset show no significant difference from free-RNA values, except for propeller-twist values. Though well within standard deviation, basepair propeller-twist in bound-RNA (AU: - 9.5°; GC: - 8.9°) have smaller negative values compared to free-RNA (AU: -12.7°; GC: - 10.9°).
Approximately 62% of the total dinucleotide steps in bound-RNA dataset are in direct contact with protein (Table 2). To examine the direct and indirect effect of protein binding, the dataset was divided into two sub-datasets: those steps that are in contact with protein (cont) and those that do not contact the protein (non-cont). The mean and standard deviation values of the step parameters and the corresponding average propeller-twist, cup and Zp value, for each of the 10 dinucleotide steps in the non-cont and cont dataset are tabulated separately in Tables 4,5 and6. The mean values of steps parameters of non-cont and cont dataset are quite similar to each other, but differ slightly from free-RNA values. Interestingly both cont and non-cont data show large standard deviation values. In addition, the correlation analysis for dinucleotide steps in the bound-RNA dataset does not show any significant difference from that of free-RNA and ADNA datasets (Additional file1).
To check for statistical significance of the differences between the step parameters of the three datasets (free-RNA, non-cont and cont) unpaired student t-test was carried out (Figure 3). No significant difference was found between cont and non-cont datasets. However, roll, slide and average propeller-twist values for several of the dinucleotide steps in the cont dataset show significant difference (P value < 0.05) from the free-RNA dataset. Some of the steps also showed significant difference between free-RNA and non-cont dataset.
Base overlap and formation of non-canonical hydrogen bonds in dinucleotide steps
Since the parameters for RR, RY and YR steps show some significant differences from Model-RNA, we calculated base stacking overlap for these steps and compared it for the various crystal datasets and the corresponding model structures (Table 7). Figure 4 illustrates the nomenclature used to refer the bases involved in basepair overlap area calculation and the different stacking patterns for RR, RY and YR steps, in A and B-form structures. In case of BDNA dataset, all three types of steps show only intra-strand base overlap, but negligible cross-strand overlap, with RR > RY > YR and the major contribution coming from exo-cyclic atoms, in all cases. In general, the overlap increases in the crystal structure steps as compared to the fibre-BDNA model. In A-form helices, the RR, RY and YR steps show distinctly different features, with RY> > RR ≈ YR. In RR steps, high intra-strand base overlap is seen in strand I (Pur-Pur stacking) while, unlike in BDNA, there is very little overlap in strand II (Pyr-Pyr stacking) and no overlap between the cross-strand bases. RY steps show high intra-strand base overlap along both strands I and II and no cross-strand overlap, with the exo-cyclic atoms making substantial contribution. YR steps in RNA and ADNA datasets are characterized by very small intra-strand contribution and stacking arises mainly due to cross-strand overlap of purine bases, with contributions from both ring and exo-cyclic atoms. Interestingly the overlap in free-RNA crystal structure steps is smaller than in the ModelRNA, for RR and RY steps. Our findings suggest that the combined effect of large negative slide and lower twist value contribute towards these overlap patterns, indicating that the interactions that determine the base stacking preferences of dinucleotide steps in RNA helices are different from those seen in DNA helices. We have therefore analyzed various dinucleotide steps to see whether the base overlap patterns are related to formation of some potential non-canonical, intra-strand or cross-strand hydrogen bonds.
Some of the geometric preferences seen in dinucleotide steps of DNA have been attributed to the presence of additional hydrogen bond interactions between the bases, particularly in oligo-A tracts[30, 34, 35]. Similarly, non-canonical hydrogen bonds between RNA bases involved in forming a dinucleotide steps can arise due to favourable intra-strand or cross-strand interactions on both major groove and minor groove side. Many such potential hydrogen bonds are possible in RNA model structure and are found to occur in crystal structures, but only those interactions that occur in more than 50% of each of the steps are discussed here. A list of such cross-strand and intra-strand interactions, along with the mean values of donor-acceptor (DA) distance, hydrogen-acceptor distance (HA) and hydrogen bond angle (DHA) in each dinucleotide step in free-RNA dataset is given in Table 8. Stick drawings of dinucleotide steps with hydrogen bonds marked for selected example (Additional file2) from crystal structures for RR, RY and YR steps are shown in Figures 5,6, and7 respectively. A complete list of hydrogen bonds present identified in RNA and DNA crystal datasets and fibre model structures is given in Additional file3. It is observed that the number of non-canonical hydrogen bonds is more in ModelRNA as compared to fibre-BDNA model. Many of these are retained, with improved hydrogen bond parameters, in the RNA crystal structures, while some potential interactions are found to occur in specific dinucleotide steps.
Among RR steps, cross-strand C-H..O hydrogen bonds are found in 85% and 59% of AA/UU and GA/UC steps respectively in freeRNA, on the minor groove side (Figure 5 and Additional file3). Similar interaction is also observed, though in smaller numbers, for bound RNA steps (non-cont and cont), as well as AA/TT and GA/TC steps in BDNA. The cross-strand N-H..O interaction between 6-amino group of Adenine and O4 atom of Uracil that is commonly seen in AA/TT steps of BDNA (89%) is favoured only in 40% of the AA/UU steps in freeRNA dataset. In GG/CC steps, intra-strand N-H..N interaction between the two Cytosine 4-amino groups shows significant presence in all A-form structures. It is not present in fibre-BDNA model but is seen in BDNA crystal structures with a slightly longer donor-acceptor (DA) distance, as compared to RNA structures. Though a similar pair of 6-amino groups of Adenine is present in AA/UU, they do not have favourable hydrogen bond geometry.
RY steps in A-form helices are characterized by large intra-strand overlap and the exo-cyclic atoms in the major groove are positioned almost above each other, leading to unfavorable N-H..O angles, in both AU/AU and GC/GC steps. However, AC/GU step show presence of a weak intra-strand N-H..N interaction between the exo-cyclic amino groups in Adenine (N6) and Cytosine (N4) (Figure 6 and Additional file3). This interaction is seen in RNA as well as in DNA helices, though the ModelRNA and fibre-BDNA model does not have a favourable geometry for this hydrogen bond. In addition, cross-strand interactions between the two purine bases are highly favoured in AC/GU and AU/AU steps in all RNA datasets and equivalent steps in DNA. A cross-strand N-H..O hydrogen bond is present in 83% of AC/GU in free-RNA, between 6-amino group of Adenine and O6 atom of Guanine. An even larger number (~90%) of AU/AU steps show cross-strand N-H..N interaction between the 6-amino group of Adenines, in both free-RNA and BDNA datasets. A combination of high negative average propeller-twist, smaller slide and positive roll values, seen in AU/AU steps in RNA favors this cross-strand interaction, while the intra-strand N-H..O hydrogen bond is relatively infrequent. GC/GC step does not show significant occurrence of any intra-strand or cross-strand interaction between their exo-cyclic groups.
In A-form helices, the large negative slide leads to high cross-strand overlap between the purine bases in YR steps (Table 7). Interestingly, the relative displacement of neighbouring bases within a strand, along with large positive roll, gives rise to favourable orientation of exo-cyclic groups in the A-form structure and hence all possible N-H..O and N-H..N hydrogen bonds are seen in large numbers. An intra-strand N-H..N interaction is present in CA/UG step between the 6-amino group of Adenine and 4-amino group of Cytosine; in 91% of the steps in freeRNA (Figure 7 and Additional file3). However, the relative orientation of the 6-amino group of Adenine and O6 oxygen atom of Guanine in the major groove does not favour a cross-strand hydrogen bond between them. Intra-strand N-H..O interaction between the 4-amino group of Cytosine and O6 oxygen atom of Guanine is seen in >60% of CG/CG steps, in both strands of RNA structures, while they are absent in BDNA dataset. Almost 100% of UA/UA steps in freeRNA and 80-95% in protein bound-RNA helices, form intra-strand N-H..O interaction between 6-amino group of Adenine and O4 oxygen atom of Uracil. A cross-strand N-H..N interaction between the two 6- amino groups of Adenine is also present in more than 60% of the A-like steps.
A rather unusual cross-strand N-H..N interaction is frequently observed between the 2-amino group of Guanine and N9 atom of the Purine base in CA/UG and CG/CG steps in A-like structures (Table 8 and Additional file3). Unlike other hydrogen bonds that are present in both model structure and RNA datasets, these N2..N9 interactions are much more favourable in the crystal dataset, with the mean Donor-Acceptor distance (DA) being ~3.6 Å (while it is ~4.1 Å in ModelRNA structure). This type of hydrogen bond is observed in ~65% of CA/UG steps, (Figure 7a). A similar type of hydrogen bond is seen in ~50% of CG/CG steps, with 31% showing a pair of reciprocal hydrogen bonds (Figure 7b). A combination of relatively higher values for negative cup, negative propeller twist and negative slide is characteristic of steps with reciprocal interactions between the 2-amino groups and N9 atoms of Guanines in CG/CG steps.
Thus, a number of non-canonical hydrogen bonds are present in the WC steps of both free as well as bound RNA crystal structures and their presence can be related to the sequence dependent geometries seen in the various dinucleotide steps. Overall, the percentage occurrence of these non-WC hydrogen bonds is smaller in bound dataset compared to free-RNA dataset.
Contrary to the generally accepted view that RNA helices are uniform and rigid, the various dinucleotide steps have characteristic features and can contribute to heterogeneity in the RNA helical regions. The intra-basepair parameters, propeller-twist and buckle, of the AU and GC basepairs show the usual preferences (large propeller-twist in AT and larger buckle in GC basepairs), which can influence the dinucleotide step geometry. Roll values differentiate the three types of steps in RNA. RY steps have small roll (except for AU/AU steps), YR have high roll and RR have intermediate roll values. Interestingly, while roll values in DNA vary from small negative to small positive, a similar trend is seen with roll for YR > RR > RY steps. In B-DNA, the difference in parameters between RR, RY and YR is attributed to the effect of exocyclic groups on slide, roll and twist. Unlike the dinucleotide steps in B-DNA, which show a large variation in twist value that is strongly correlated with roll and moderately with slide, the twist values of steps in RNA helices cluster within a small range and show a significant correlation only with slide. The larger positive roll and negative slide values lead to a large number of favourable intra as well as cross-strand interactions involving the exo-cyclic groups, particularly in YR steps. Interestingly the slide values of WC steps show a significant negative correlation with average propeller-twist in RNA and a positive correlation for B-DNA. Steps with large average propeller-twist (AA/UU, AU/AU and UA/UA) have smaller negative slide values in RNA. The proximity of atoms on major groove side, arising due to high roll and negative propeller-twist prevents large negative slide. Thus, the slide parameter is directly influenced by propeller-twist of the constituent basepairs, but as the basepairs become near planar (average propeller twist ≈ 0°), slide becomes large negative in RNA and small positive in B-DNA steps. Overall, the six dinucleotide step parameters in RNA helices show mean values as well as correlated variations that are different from those observed in B-form DNA[2, 4, 6, 19, 49].
The specificity in RNA-protein interaction is thought to be mainly brought about by the exposed bases which are present at helix termini, in bulges and loops. In our analysis, we find that the interactions occur mainly with the phosphate backbones and very few interactions are seen between proteins and the base atoms. In our study, majority of steps in the protein contacting (cont) and non-contacting (non-cont) dataset, showed statistically significant difference in comparison to free-RNA dataset for roll, slide and average propeller-twist. This suggests that, apart from protein induced conformation change on direct contact, it can allosterically affect the steps that are not in direct physical contact.
Recently, various additional hydrogen bond interactions have been reported from the RNA double helical regions, between the base and phosphate group oxygens (BPh)[26, 27, 43]. Also, the presence of weak hydrogen bonds between cross-strand amino groups in AA/TT and GA/TC steps of B-form DNA are well documented in crystal structure and supported by theoretical quantum chemical calculations. Similarly the presence of C-H..O interactions were reported in B-DNA crystal structures. The presence of cross-strand C-H..O interactions in AA/UU and GA/UC steps and N-H..N interaction in AU/AU steps have been reported from MD simulations of A-RNA duplex sequences. Two other cross-strand N-H..O interactions on the minor groove side in AG/CU and GG/CC steps, reported in the MD studies, occur in less than 20% of these steps in our RNA dataset. However, our analysis has confirmed the presence of other cross-strand interactions and identified novel cross-strand and intra-strand hydrogen bonds that can potentially provide added stability to the RNA dinucleotide step. The cross-strand C-H..O interaction in AT basepair containing steps, AA/TT and GA/TC, reported on the minor groove side of B-DNA crystal structures are surprisingly also found to occur in a majority of AA/UU and GA/UC steps in RNA. In B-DNA helices the AA/TT and GA/TC steps have large negative average propeller twist, but near zero roll and slide, while the AA/UU and GA/UC steps in RNA have large negative average propeller-twist, but moderately positive roll and large negative slide values. Thus, it appears that the same interaction is brought about by a combination of high negative propeller-twist and two different roll-slide geometries, both of which bring the pairing atoms close. Similarly, intra-strand N-H..N interactions in GG/CC and AC/GU steps, cross-strand N-H..O interactions in AC/GU step and the cross-strand N-H..N interactions in AU/AU step are present in both A-form and B-form helices. In RNA, when compared to RR and RY steps, the YR steps show a larger number of these potential hydrogen bonds, due to their unique cross-strand stacking. The weak interaction identified in this study can contribute to stacking and thus to overall stability of the steps. This is in agreement with the results of stacking energy calculated using QM, where the values for RY and YR steps are comparable though the base overlaps as shown in Table 7 are considerably lower for YR steps.
Among YR steps, CG/CG and CA/UG steps, in addition to intra-strand N-H..O and N-H..N hydrogen bonds, they also have an unusual N-H..N interaction between 2-amino groups of Guanines and N9 atoms of Purine bases in these dinucleotide steps. Hydrogen bonds are generally associated with the electronegative character of the donor and acceptor atoms. Electrostatic potential (ESP) derived charges, as well as partial charges in AMBER and CHARMM force fields assign near zero charges to the N9 atoms in Adenosine and Guanosine[52–54]. However partial charges calculated for Adenosine and Guanosine using Natural Bond Orbital (NBO) analysis indicate that the N9 atom is quite negative (Additional file4). Thus, the presence of this potential hydrogen bond needs to be further examined by quantum chemical methods.
It should also be mentioned that x-ray determined crystal structures do not have coordinates of hydrogen atoms. Programs that add hydrogen atoms to the nucleotide ring atoms as well as the N atoms in the pendent amino group, place the hydrogen atoms in the plane of the base, though many QM studies suggest that these amino groups can have pyramidal geometry with hydrogen atoms being out-of-plane[35, 56–59]. The introduction of non-planar pyramidal amino hydrogen atoms can facilitate further improvement in the geometry of N-H..N as well as N-H..O hydrogen bonds discussed here. The hydrogen bonds reported here are more prevalent in free RNA structures than protein-bound RNA helices, where they may be replaced by interactions with proteins or nearby water molecules. Flanking basepairs can also affect the formation of these weak hydrogen bonds. This sequence dependency can be studied by analysing all possible tetramer sequences in helices with the dinucleotide step of interest in the centre. However, the currently available RNA crystal structures do not have sufficient representation of all possible tetramer sequences, for a meaningful analysis.
Apart from using the well-known Watson-Crick edge, a base can pair with other bases using the Hoogsteen or Sugar edges. In our crystal dataset we focused only on the dinucleotide steps formed by WC basepair that belong to the cisWW family (~74%), to compare with equivalent dinucleotide steps in DNA helices. However, more than 32 types of cisWW steps, containing at least one non-Watson-Crick basepair, such as GU, AG and UU, are also present and constitute ~17% of the total number of steps, while ~9% of the steps contain bases that are paired along the Hoogsteen or Sugar edge. Hence, to get a complete picture of the RNA helical geometries, the non-canonical basepair containing steps were analysed, but the small number present for each of these step types in crystal structures poses a challenge in arriving at any statistically significant result.
Our analysis of the non-redundant RNA crystal structures shows that the RNA dinucleotide steps have characteristic sequence dependent variations. Overall, the steps show features that are attributes of their being of RR, RY or YR type. Several cross-strand and intra-strand potential hydrogen bonds are found to be highly prevalent in RNA helices and can be attributed to the observed geometrical preferences of various dinucleotide steps. Unusual cross-strand interactions are found to be present in the CA/UG and CG/CG steps, between 2-amino groups of Guanines and N9 atoms of Purine bases that are associated with the unique geometry of YR steps. The various dinucleotide steps in RNA bound to proteins show some significant differences in their dinucleotide parameters, from those in free RNA, while retaining most of the gross features, as well as the non-canonical cross-strand and intra-strand interactions.
Calladine C, Drew H, Luisi B, Travers A: Understanding DNA, Third Edition: The Molecule and How it Works. 2004, Academic Press
Bhattacharyya D, Bansal M: Local variability and base sequence effects in DNA crystal structures. J Biomol Struct Dyn. 1990, 8 (3): 539-572. 10.1080/07391102.1990.10507828.
Gorin AA, Zhurkin VB, Olson WK: B-DNA twisting correlates with base-pair morphology. J Mol Biol. 1995, 247 (1): 34-48. 10.1006/jmbi.1994.0120.
Marathe A, Karandur D, Bansal M: Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs. BMC Struct Biol. 2009, 9: 24-10.1186/1472-6807-9-24.
Packer MJ, Hunter CA: Sequence-dependent DNA structure: the role of the sugar-phosphate backbone. J Mol Biol. 1998, 280 (3): 407-420. 10.1006/jmbi.1998.1865.
Suzuki M, Amano N, Kakinuma J, Tateno M: Use of a 3D structure data base for understanding sequence-dependent conformational aspects of DNA. J Mol Biol. 1997, 274 (3): 421-435. 10.1006/jmbi.1997.1406.
Esguerra M, Olson WK: Albany 2011, Conversation 17: 2011. RNASTEPS, An Online Database of Sequence-dependent Deformability of RNA Helical Regions. 2011, Albany: Adenine Press
Das J, Mukherjee S, Mitra A, Bhattacharyya D: Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis. J Biomol Struct Dyn. 2006, 24 (2): 149-161. 10.1080/07391102.2006.10507108.
Nagaswamy U, Larios-Sanz M, Hury J, Collins S, Zhang Z, Zhao Q, Fox GE: NCIR: a database of non-canonical interactions in known RNA structures. Nucleic Acids Res. 2002, 30 (1): 395-397. 10.1093/nar/30.1.395.
Nagaswamy U, Voss N, Zhang Z, Fox GE: Database of non-canonical base pairs found in known RNA structures. Nucleic Acids Res. 2000, 28 (1): 375-376. 10.1093/nar/28.1.375.
Xin Y, Olson WK: BPS: a database of RNA base-pair structures. Nucleic Acids Res. 2009, 37 (Database issue): D83-88.
Williamson JR: Induced fit in RNA-protein recognition. Nat Struct Biol. 2000, 7 (10): 834-837. 10.1038/79575.
Ellis JJ, Jones S: Evaluating conformational changes in protein structures binding RNA. Proteins. 2008, 70 (4): 1518-1526.
Jones S, Daley DT, Luscombe NM, Berman HM, Thornton JM: Protein-RNA interactions: a structural analysis. Nucleic Acids Res. 2001, 29 (4): 943-954. 10.1093/nar/29.4.943.
Morozova N, Allers J, Myers J, Shamoo Y: Protein-RNA interactions: exploring binding patterns with a three-dimensional superposition analysis of high resolution structures. Bioinformatics. 2006, 22 (22): 2746-2752. 10.1093/bioinformatics/btl470.
Ellis JJ, Broom M, Jones S: Protein-RNA interactions: structural analysis and functional classes. Proteins. 2007, 66 (4): 903-911.
Lu XJ, Shakked Z, Olson WK: A-form conformational motifs in ligand-bound DNA structures. J Mol Biol. 2000, 300 (4): 819-840. 10.1006/jmbi.2000.3690.
Marathe A, Bansal M: An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression. BMC Struct Biol. 2011, 11: 1-10.1186/1472-6807-11-1.
Olson WK, Gorin AA, Lu XJ, Hock LM, Zhurkin VB: DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. Proc Natl Acad Sci U S A. 1998, 95 (19): 11163-11168. 10.1073/pnas.95.19.11163.
Stefl R, Skrisovska L, Allain FH: RNA sequence- and shape-dependent recognition by proteins in the ribonucleoprotein particle. EMBO Rep. 2005, 6 (1): 33-38. 10.1038/sj.embor.7400325.
Draper DE: Themes in RNA-protein recognition. J Mol Biol. 1999, 293 (2): 255-270. 10.1006/jmbi.1999.2991.
Xia T, SantaLucia J, Burkard ME, Kierzek R, Schroeder SJ, Jiao X, Cox C, Turner DH: Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with watson - crick base pairs. Biochemistry. 1998, 37 (42): 14719-14735. 10.1021/bi9809425.
Turner DH, Sugimoto N, Freier SM: RNA structure prediction. Annu Rev Biophys Biophys Chem. 1988, 17: 167-192. 10.1146/annurev.bb.17.060188.001123.
Turner DH, Bevilacqua PC: Thermodynamic Considerations for Evolution by RNA. The RNA World. Edited by: Atkins JF, Gesteland RF. 1993, New York: Cold Spring Harbor Press, 447-464.
Lu XJ, Olson WK, Bussemaker HJ: The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex. Nucleic Acids Res. 2010, 38 (14): 4868-4876. 10.1093/nar/gkq155.
Zgarbová M, Jurečka P, Banáš P, Otyepka M, Šponer JE, Leontis NB, Zirbel CL, Šponer J: Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of Key base–phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations. J Phys Chem. 2011, 115 (41): 11277-11292. 10.1021/jp204820b.
Zirbel CL, Šponer JE, Šponer J, Stombaugh J, Leontis NB: Classification and energetics of the base-phosphate interactions in RNA. Nucleic Acids Res. 2009, 37 (15): 4898-4918. 10.1093/nar/gkp468.
Desiraju GR, Steiner T: The Weak Hydrogen Bond: In Structural Chemistry and Biology (International Union of Crystallography Monographs on Crystallography). 2001, USA: Oxford University Press
Weiss MS, Brandl M, Suhnel J, Pal D, Hilgenfeld R: More hydrogen bonds for the (structural) biologist. Trends Biochem Sci. 2001, 26 (9): 521-523. 10.1016/S0968-0004(01)01935-1.
Ghosh A, Bansal M: C-H.O hydrogen bonds in minor groove of A-tracts in DNA double helices. J Mol Biol. 1999, 294 (5): 1149-1158. 10.1006/jmbi.1999.3323.
Mandel-Gutfreund Y, Margalit H, Jernigan RL, Zhurkin VB: A role for CH…O interactions in protein-DNA recognition. J Mol Biol. 1998, 277 (5): 1129-1140. 10.1006/jmbi.1998.1660.
Panigrahi SK, Desiraju GR: Strong and weak hydrogen bonds in the protein-ligand interface. Proteins. 2007, 67 (1): 128-141. 10.1002/prot.21253.
Singh SK, Babu MM, Balaram P: Registering alpha-helices and beta-strands using backbone C-H…O interactions. Proteins. 2003, 51 (2): 167-171. 10.1002/prot.10245.
Nelson HC, Finch JT, Luisi BF, Klug A: The structure of an oligo(dA).oligo(dT) tract and its biological implications. Nature. 1987, 330 (6145): 221-226. 10.1038/330221a0.
Shatzky-Schwartz M, Arbuckle ND, Eisenstein M, Rabinovich D, Bareket-Samish A, Haran TE, Luisi BF, Shakked Z: X-ray and solution studies of DNA oligomers and implications for the structural basis of A-tract-dependent curvature. J Mol Biol. 1997, 267 (3): 595-623. 10.1006/jmbi.1996.0878.
Besseova I, Banas P, Kuhrova P, Kosinova P, Otyepka M, Sponer J: Simulations of A-RNA duplexes. The effect of sequence, solute force field, water model, and salt concentration. J Phys Chem B. 2012, 116 (33): 9899-9916. 10.1021/jp3014817.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The protein data bank. Nucleic Acids Res. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.
Ray SS, Halder S, Kaypee S, Bhattacharyya D: HD-RNAS: An automated hierarchical database of RNA structures. 2012, Genetics: Frontiers in, 3-
Leontis N, Zirbel C: Nonredundant 3D Structure Datasets for RNA Knowledge Extraction and Benchmarking. RNA 3D Structure Analysis and Prediction, Volume 27. Edited by: Leontis N, Westhof E. 2012, Springer Berlin Heidelberg, 281-
Bansal M, Bhattacharyya D, Ravi B: NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures. Comput Appl Biosci. 1995, 11 (3): 281-287.
Arnott S: Polynucleotide secondary structures: an historical perspective. Oxford Handbook of Nucleic Acid Structure. Edited by: Neidle S. 1999, Oxford Press, 1-38.
Lu XJ, Olson WK: 3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures. Nucleic Acids Res. 2003, 31 (17): 5108-5121. 10.1093/nar/gkg680.
Leontis NB, Westhof E: Geometric nomenclature and classification of RNA base pairs. RNA. 2001, 7 (4): 499-512. 10.1017/S1355838201002515.
Bhattacharyya D, Bansal M: A self-consistent formulation for analysis and generation of non-uniform DNA structures. J Biomol Struct Dyn. 1989, 6 (4): 635-653. 10.1080/07391102.1989.10507727.
El Hassan M, Calladine CR: Conformational characteristics of DNA: empirical classifications and a hypothesis for the conformational behaviour of dinucleotide steps. Philos Trans Royal Soc A. 1997, 355 (1722): 43-100. 10.1098/rsta.1997.0002.
Word JM, Lovell SC, Richardson JS, Richardson DC: Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation. J Mol Biol. 1999, 285 (4): 1735-1747. 10.1006/jmbi.1998.2401.
Collaborative Computational Project N: The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr. 1994, 50: 760-763. 10.1107/S0907444994003112.
MATLAB: Version R2010b. 2010, Natick, Massachusetts: The MathWorks Inc
Svozil D, Kalina J, Omelka M, Schneider B: DNA conformations and their sequence preferences. Nucleic Acids Res. 2008, 36 (11): 3690-3706. 10.1093/nar/gkn260.
Šponer J, Kypr J: Close mutual contacts of the amino groups in DNA. Int J Biol Macromol. 1994, 16 (1): 3-6. 10.1016/0141-8130(94)90003-5.
Sponer J, Hobza P: Bifurcated hydrogen bonds in DNA crystal structures. An ab initio quantum chemical study. J Am Chem Soc. 1994, 116 (2): 709-714. 10.1021/ja00081a036.
Brooks BR CLB, Jr ADM, Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, et al: CHARMM: the biomolecular simulation program. J Comput Chem. 2009, 30 (10): 1545-1614. 10.1002/jcc.21287.
Case D, Cheatham T, Darden T, Gohlke H, Luo R, Merz K, Onufriev A, Simmerling C, Wang B, Woods R: The Amber biomolecular simulation programs. J Comput Chem. 2005, 26 (16): 1668-1688. 10.1002/jcc.20290.
Williams DE: Net Atomic Charge and Multipole Models for the ab Initio Molecular Electric Potential. Reviews in Computational Chemistry. 2007, John Wiley & Sons, Inc, 219-271.
Carpenter JE, Weinhold F: Analysis of the geometry of the hydroxymethyl radical by the "different hybrids for different spins" natural bond orbital procedure. J Mol Struct: THEOCHEM. 1988, 169: 41-62.
Sponer J, Leszczynski J, Hobza P: Electronic properties, hydrogen bonding, stacking, and cation binding of DNA and RNA bases. Biopolymers. 2001, 61 (1): 3-31. 10.1002/1097-0282(2001)61:1<3::AID-BIP10048>3.0.CO;2-4.
Sponer J, Hobza P: Nonplanar geometries of DNA bases. Ab initio second-order Moeller-Plesset study. J Phys Chem. 1994, 98 (12): 3161-3164. 10.1021/j100063a019.
Mukherjee S, Majumdar S, Bhattacharyya D: Role of hydrogen bonds in protein-DNA recognition: effect of nonplanar amino groups. J Phys Chem B. 2005, 109 (20): 10484-10492. 10.1021/jp0446231.
Bandyopadhyay D, Bhattacharyya D: Estimation of strength in different extra Watson-Crick hydrogen bonds in DNA double helices through quantum chemical studies. Biopolymers. 2006, 83 (3): 313-325. 10.1002/bip.20542.
The authors are also grateful to Dr. Arvind Marathe for useful discussions.
The authors declare that they have no competing interest.
DB and MB conceived the project. SK carried out the analysis. All authors participated in the writing of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Overall cross-correlation between dinucleotide step parameters for all steps comprising of canonical WC basepairs.
Additional file 1: For free-RNA dataset (N = 797): correlation coefficient (r) values ≥ 0.15 are significant at 99.9% confidence level. For bound-RNA dataset (N = 2531): correlation coefficients (r) values ≥ 0.10 are significant at 99.9% confidence level. For ADNA dataset (N = 195): correlation coefficients (r) values ≥ 0.32 are significant at 99.9%. For BDNA dataset (N = 212): correlation coefficients (r) values ≥ 0.23 are significant at 99.9% confidence level. ‘r’ values significant at 99.9% level, in each dataset, are shown in bold. (DOC 48 KB)
Additional file 3:Comparison of non Watson-Crick hydrogen bonds observed in all RNA and DNA helices. Mean and SD values of donor-acceptor distances and angles observed for cross-strand and intra-strand hydrogen bond in more than 50% of the 10 dinucleotide steps are listed. Hydrogen bond parameters in various model structures are also given. (XLS 51 KB)
Additional file 4: The partial charges derived from Natural Bond Orbital (NBO) and Electrostatic potential (ESP) calculations are listed, along with those used in AMBER and CHARMM force fields. The large negative charge assigned by NBO calculation is highlighted in bold for the N9 atoms of adenine and guanine bases, which are involved in an unusual cross-strand hydrogen bond with the 2-amino group of guanine in the CA/UG and CG/CG steps respectively (shown in Figure 7a and b). (DOC 84 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kailasam, S., Bhattacharyya, D. & Bansal, M. Sequence dependent variations in RNA duplex are related to non-canonical hydrogen bond interactions in dinucleotide steps. BMC Res Notes 7, 83 (2014). https://doi.org/10.1186/1756-0500-7-83
- Hydrogen bond