- Research article
- Open Access
Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs
BMC Structural Biology volume 9, Article number: 24 (2009)
An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.
Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.
The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.
Watson and Crick proposed the double-helical structure for DNA in 1953, and almost simultaneously two forms were postulated on the basis of fibre diffraction analysis – the B-form DNA corresponded to the Watson and Crick structure and was found to occur in conditions of high humidity and low salt concentration, while the A-form occurred in conditions of lower humidity and higher salt concentration. Gross structural features for these as well as other polymorphic forms of DNA were refined during the next 25 years, using fibre diffraction data [1, 2]. The two forms of DNA were mainly characterised in terms of features such as sugar pucker [3–5], glycosyl torsion angle , base pair orientation and the groove widths , apart from the helical parameters rise and twist. However, it was only in the 1980s that the atomic details of the two forms were characterised. The first crystal structure of a B-form DNA was solved in 1981 , and was found to have significant sequence-dependant variability, with an average roll per dinucleotide step of 0.5 ± 5.2°, an average local helical twist of 35.6 ± 4.4° and an average slide of 0.2 ± 0.5 Å. Subsequent analyses of other crystal structures confirmed the sequence dependent effects observed here [9–12]. A-DNA, which was first crystallised by McCall , was found to have an average roll of 6.8 ± 2.6°, an average local helical twist of 30.8 ± 1.2° and an average slide of -1.5 ± 0.4 Å [10, 12, 14]. However, since the overall features of the two crystal structures were close to the fibre models of B and A forms, it was assumed that the two forms correspond to two stable minima that the DNA could assume and transition from one form to another would involve some energetic costs. As crystallographic methods improved and the number and variety of x-ray crystal structures of DNA increased, this idea began to lose ground. While most oligonucleotide structures solved during 1980–2000 had roll and twist values that ranged from exclusively A-like to exclusively B-like, a few appeared to show features intermediate between A-DNA and B-DNA, to a mixture of both types. [15–18]. Thus it appeared that A and B-form DNA were not well separated stable minima, and the dinucleotide steps in oligomeric DNA could assume conformations that ranged from B-like to intermediate to A-like . In addition, several other forms of synthetic DNA were also solved, which did not fit the canonical A-like or B-like conformation . Against this wide ranging polymorphism of the double-helical DNA molecule, particularly at the dinucleotide step level, the RNA duplex crystal structures, that were solved around the same time [19–24], stood out for their rigidity, and their conformational proximity to the A-RNA fibre model, independent of the sequence. In this study, we have analysed a large dataset of free RNA oligomers to verify the conformational rigidity of Watson-Crick basepaired RNA duplexes and then used it as a template against which to measure the A-like characteristics of each dinucleotide step as well as overall structure of both free and protein-bound DNA.
Several studies in the late 1990s also suggested that not only the classical B-form of DNA but also the A-form had biological relevance. Based on a comparison of free B-DNA oligomers and protein-bound DNA, it was suggested [25, 26] that protein binding causes DNA to assume A-like or an A-B intermediate conformation in terms of roll and twist. Subsequently it was shown that a new parameter, Z p , could be used to discriminate between A-like or B-like dinucleotide steps more reliably than roll or twist, and that entire structures could be classified as A-like or B-like in terms of their Z p values, irrespective of the local variations in their roll and twist values [14, 27]. Lu et al  highlighted the fact that in DNA structures bound to a few prominent protein families, the protein-bound region was induced to take up an A-like conformation as defined by Z p . However, the above mentioned studies [25, 26], that compared free and bound DNA, considered the overall B-form of the free oligomers as a reference, and not the inherent 'A-philicity' [28–32] of dinucleotide steps in the bound region. Given that atleast in a few cases, the putative binding region is known to assume an A-like conformation in its free form [33, 34], inclusion of A-DNA oligomer structures also in the analysis might provide better insights into the intrinsic preferences of a DNA sequence and help distinguish these from protein induced structural effects. Only one study compares the free and protein-bound forms of DNA, taking into consideration the A-form of DNA . Several other studies have implicated the variations in roll, especially at pyrimidine-purine steps, to be responsible for DNA bending and curvature [25, 36–41], in ways critical for the binding of the protein.
While the DNA dinucleotide steps were under scrutiny for their role in specifically binding to a protein, the DNA backbone was also shown to be involved in more than 50% of all the contacts between amino acids and the DNA in regulatory protein-DNA complexes . Hence several studies have also focused on how the variations in the DNA backbone might act as an indirect readout signal for protein recognition and binding [43–49]. In DNA oligomers, the sugar phosphate backbone was believed to be rigid, compared to the variation in local step geometry, defined by two neighbouring basepairs. The sugar ring assumed C3'-endo conformation in A-DNA  and C2'-endo in B-DNA . The related backbone torsion angle δ was found to assume values of about 84° for A-DNA  and about 128° for B-DNA [4, 5]. The torsion angles ε and ζ were observed to assume two conformations-BI and BII in B-DNA  but only the BI conformation in A-DNA . α and γ were found to show anticorrelated variation in A-form duplexes [3, 24, 50, 51], but were generally found to take up the g-, g+ conformation in B-DNA . However, recent studies have shown that unlike oligomers, the backbone in a significant proportion of nucleotides in bound DNA assumes non-classical conformations . There have also been attempts to analyse the backbone torsion angles, taking into account the correlation between more than two torsion angles and group them into seven distinct states [53, 54]. In this study, we have adapted this methodology  and analysed the variation in backbone parameters with respect to variation in dinucleotide step parameters across different datasets.
A crucial question of biological relevance is how the variations in DNA structure at the basepair, base-step and backbone level contribute to the overall structure of the molecule, and its implications for protein binding. A related question is how changes caused by protein binding at the local structural level affect the overall DNA structure. There have been efforts to go beyond dinucleotide steps and analyse the properties of all possible tetranucleotide, hexanucleotide and octanucleotide fragments using molecular dynamics simulation studies [53, 55, 56]. However, most of the high resolution DNA double helical crystal structures, especially those of free DNA, are too short in length, to allow a meaningful statistical analysis of all possible trinucleotide or higher order steps. The other approach is to try and quantify the overall DNA structure, in terms of parameters such as DNA curvature, bendability or stability. The importance of DNA curvature was first realized when it was observed that even unbound genomic DNA could have a well-defined, inherent curvature [57, 58]. Since most of the curved DNA observed in the early days were observed to have stretches of adenines, the initial models of DNA curvature, such as the 'wedge model'  and the 'junction model' , traced the origin of curvature to the presence of A-tracts, in phase with the DNA helical repeat. However, these models had to be abandoned when it was shown that sequences lacking in AA dinucleotides also adopted a curved structure . Thus new models which took into account variation in the geometries of all ten dinucleotide steps were proposed [62–64]. However, owing to difficulties in tracing a uniform path for the DNA axis in three dimensions, there is no standard methodology for calculating DNA curvature, despite its obvious importance in biological functions. Various measures of quantifying DNA curvature such as the radius of a circle fitted to the basepair centres projected onto a plane [63, 65], the ratio of the end-to-end distance of the DNA molecule to the actual path traced by the DNA axis [63, 65–69], ratio of the moments of inertia of an ellipsoid fitted to the molecule [63, 65–71] as well as the angles between two local helix axes vectors corresponding to two successive dinucleotide steps have been proposed and implemented [64, 69, 72]. However, each of these methods has its advantages and limitations, and no single method can unambiguously quantify all possible curved conformations adopted by DNA molecules. Hence a combination of all or several of these methods along with a close inspection of the local level distortions is required to fully understand the curvature of any given structure.
In this study, we have analysed an exhaustive dataset of protein-DNA complexes, and compared it with a complete, high resolution dataset of free DNA oligomers, without pre-classifying them as A-DNA or B-DNA. We have also separately analysed a dataset of DNA bound by proteins via a Helix-Turn-Helix (HTH) motif. The HTH motif is not only the most well-characterised, but also the most commonly occurring DNA-binding motif, and is present in a wide range of transcription factors. The HTH motif consists of two alpha helices linked by a turn region that protrudes out of the surface of the protein [73–75]. The second helix, usually referred to as the 'recognition helix', fits into the major groove of the DNA, and is involved in direct or indirect interactions with the DNA [74–76]. While the HTH motif has been studied extensively, the structural features of the DNA to which it binds have not been analysed in detail. The present analysis provides some interesting insights into the conformational flexibility of the DNA molecule, and reveals that many of the conformations observed in bound DNA, both at the local dinucleotide step level, and the gross structural level, are also accessible to unbound DNA, while a few conformations are solely induced by protein binding.
The structural parameters of three datasets of DNA – free oligomers, protein-bound DNA (excluding HTH motif-bound DNA) and HTH motif-bound DNA, and one dataset of RNA oligomers were analysed in order to gain a complete perspective of the features of DNA both within each set and also across the sets. As RNA is known to assume only A-like conformation, the RNA dataset was used as a reference point for A-like conformation and also to characterise the basepair effects from those due to the ribose sugar ring in RNA. The RNA dataset was observed to be rich in steps containing the G:C basepair, and had remarkably low percentage of steps containing only the A:U basepair (table 1). The free DNA dataset consists of a large proportion of the steps GG (23.9%) and CG (18.0%) (table 1). A significant number of these steps were found to occur in structures which were classified by the Nucleic Acid Database  as "A-DNA". A large number of these steps are indeed found to have high Z p values in the present analysis, matching our criteria for an A-like dinucleotide step, as defined in the next section. The free dataset also contains 5 structures with the Drew-Dickerson sequence d(CGCGAATTCGCG). These and other A-tract containing sequences primarily contribute to the high occurrence of AA steps (11.1%) in the free dataset. The HTH dataset consists of DNA bound by a wide variety of proteins ranging across 22 SCOP  classes, and includes 3 ternary TATA Binding Protein-Transcription Factor IIB-TATA-box (hereafter referred to as TBP-TFIIB-TA-DNA) complexes and 6 Catabolite Activator Protein-DNA (hereafter referred to as CAP-DNA) complexes (additional file 1). In the TBP-TFIIB-TA-box DNA ternary complexes, the HTH motif is present in the transcription factor TFIIB, which binds to the DNA immediately upstream of the TATA-box region. The complex dataset contains 8 TATA Binding Protein-DNA (hereafter referred to as TBP-DNA) complexes, which lack TFIIB, and hence have been excluded from the HTH dataset (additional file 1). Interestingly, the protein-bound datasets also have a significant proportion of CA/TG steps, which have been implicated in the kinks observed in several structures .
Variations of the dinucleotide step parameters
Among the six dinucleotide step parameters that measure the relative rotational and translational motions between adjacent basepairs about the x, y and z-axis (see 'Methods'), tilt, shift and rise were observed to have very little variation within and across the three DNA datasets, and so have not been reported here. On the other hand, in conformity with earlier studies [25–27], the parameters roll, twist and slide, as well as the parameter Z p (described in the 'Methods' section), were found to be excellent indicators for analysing the sequence dependent conformational flexibility of a DNA molecule. To highlight the characteristic features of each dinucleotide step in free as well as bound DNA, the dinucleotide step parameters Z p and slide are listed in tables 2 and 3, while figure 1 shows the variation of Z p versus slide. The corresponding values for roll and twist are listed in tables 4 and 5, while the variation of Z p versus roll is shown in figure 2.
The RNA oligomer dataset assumes mean values of high positive Z p (2.2 ± 0.3 Å), negative slide (-1.5 ± 0.4 Å), high roll (9.0 ± 4.0°) and low twist (31.1 ± 3.4°), all close to the values observed for the fibre models of A-form DNA helices . The low values of the standard deviations for all four parameters for individual steps, as well as for the entire dataset, confirms the conformational rigidity of the RNA structures. The sugar-phosphate backbone torsion angles χ and δ and the phase angle P were also observed to assume A-DNA fibre model-like conformation, for the entire dataset. Even steps that have previously been reported to be A-phobic in DNA are observed to be entirely A-like in RNA. This confirms the observation that the presence of even a single ribose sugar causes the entire structure to assume A-like conformation , while the presence of uracil in place of thymine also facilitates the A-form structure, particularly for AA/UU, GA/UC and AG/CU type of steps. Thus, the RNA dataset, with its well-defined and rigid boundaries stands in sharp contrast to the free and protein-bound DNA datasets and its parameters can be used as a criteria to define A-like conformation in DNA.
For the free DNA dataset, as seen in figure 1, a distinct bimodal distribution is observed for Z p and slide. The two distinct clusters for this dataset arise primarily because Z p assumes two distinctly different values with a clear separation between them. Using the RNA dataset as a template, we assigned as A-like, those DNA steps that lie within three standard deviations of the mean Z p value for the RNA dataset viz. Z p > 1.3 Å. The boundary for B-like conformation was assigned at 0.8 Å, based on visual inspection of figure 1. The points with Z p between 0.8 and 1.3 Å were considered to have intermediate conformation. With the exception of one datapoint, all free DNA oligomer steps can be classified as A-like or B-like in terms of Z p . The A-like cluster has a mean of 2.1 Å, close to the value for RNA oligomers and the other cluster has a mean value of -0.5 Å, which corresponds to the fibre model value for B-DNA (-0.6 Å) . For all the four datasets, slide was observed to correlate well with Z p for the overall data (figure 1), as well as for individual dinucleotide sequences. In contrast, roll does not show a significant correlation with Z p (figure 2), nor do the roll and twist parameters (additional file 2, figure 1) show any bimodal character.
A-philicity of dinucleotides in DNA structures
Efforts have been made several decades back to characterise individual dinucleotide steps as being A-phobic (or B-philic) (AA, CA and GA) or A-philic (GG, AG and AC) [28–31] on the basis of their ability to induce a B to A transition in solution. A more recent study on a larger dataset  reclassified the GA step as neutral and the AG step as B-philic. Our analysis confirms that the AA and GG steps are highly B-philic and A-philic respectively, in terms of their preference for Z p and slide values. The GA and AG steps in our dataset also show high preference for B-like conformation, except for very few steps which have A-like parameters. AC, which has earlier been reported to be A-philic, as well as AT, GC, TA and CG steps were observed to display both A-like or B-like values, though AT and TA show a preference for B-form.
CA/TG steps assume both A and B types of conformations in terms of Z p and slide, the steps with B-like values occurring in structures with A-tracts that are overall B-like, and steps with A-like values occurring in structures with a large number of C:G basepair containing steps that have high, or A-like, Z p values. B-like CA steps themselves assume two types of conformations in terms of roll, twist and the backbone torsion angles ε and ζ, namely BI and BII, thus confirming the highly flexible nature of this dinucleotide step, with no marked preference for A or B like geometries.
Thus it appears that in free DNA oligomers, the overall structure assumes A-like or B-like conformation depending on its sequence, particularly the proportion of AA/TT and GG/CC steps. Only AA, and to a lesser extent GA steps show strong preference for B-form, while GG is truly A-philic. All the other dinucleotide steps do not appear to have a strong intrinsic preference for A-like or B-like conformation, but assume a particular conformation depending on the conformation of neighbouring steps, as suggested by recent solution studies .
In the protein-bound DNA datasets, most of the structures were found to exclusively have B-like values for Z p , if the above-mentioned criteria for A-like and B-like DNA is used. Unlike the free dataset, no structure from the complex or HTH dataset was observed to have entirely A-like conformation. Even for an A-philic step such as GG, for which 92.8% of the steps in the free dataset take up an A-like geometry, about 90.0% of the datapoints in the complex dataset and about 71.6% of the datapoints in the HTH dataset were observed to have B-like values of Z p , with only 5.0% and 12.7% of datapoints respectively, showing an A-like value for Z p . Only a few steps in the DNA-binding region of some structures were observed to have A-like or near A-like characteristics. These complexes belong to a few specific families, such as the polymerases, endonucleases and transposases, and the structural features of these duplexes have been described in the 'Discussion' section.
Roll and twist are not good discriminators of A-form versus B-form
Roll and twist span a very wide range of values for the three DNA datasets, as evident from their values listed in tables 4 and 5. Unlike Z p and slide, there is no clear bimodal distribution for roll and twist for the free dataset, with the values varying in a continuous negatively correlated fashion, from high negative roll and very large twist to positive roll and low twist (additional file 2, figure 1). In the free DNA dataset, steps which have been classified as A-like or B-like based on their Z p values, have been listed separately in tables 4 and 5.
As mentioned above, CA steps show three types of conformations-one in which Z p , slide, roll and twist have typical A-like values and two different conformations, wherein Z p is B-like. CA steps with B-like Z p values are observed to assume either normal slide and twist with positive roll or high positive slide, large twist and negative roll. This bimodal distribution of the B-like CA steps has been observed in several previous studies [10, 65, 81, 82]. These steps also show a correlated variation in the backbone torsion angles, ε and ζ in both strands, with the low twist and positive roll steps having ε and ζ in the t, g- (or BI) conformation, while the large twist and negative roll steps have ε and ζ in the g-, t (or BII) conformation [9, 83]. When a CA step in the BII conformation occurs adjacent to an AG step such that it forms a CAG triplet, the AG step is often observed to have a high roll and a very low twist. This feature is observed in several DNA structures irrespective of whether the steps have bound ions [84, 85] or are present free . These CA steps and the adjacent AG steps do not show any correlated variation in Z p and slide, which have B-like values, with these CA steps being characterised by large positive slide values. The AG steps occurring adjacent to other steps do not assume this conformation. The high roll and low twist values of these AG steps, which are B-like in terms of Z p , skew the averages for roll (3.5 ± 4.7°) and twist (32.5 ± 7.0°) to A-like values.
In addition to the CA step, 8 of the 13 'B-like' GC steps are also observed to assume the BII conformation for one or both of the guanine backbone torsion angles, and have a corresponding negative value of roll and a large value of twist. As a result, 'B-like' GC steps have a negative average value for roll (-3.7 ± 5.7°) and a large average value for twist (38.2 ± 2.8°). AA steps, which are exclusively B-like in terms of Z p , have mean roll and twist values of 0.2 ± 4.0 and 36.1 ± 3.9 respectively, indicating that these steps are B-like in terms of roll and twist also. Among the other dinucleotide steps, the GG and CG steps have A-like mean values for both roll and twist, irrespective of whether their Z p value is A-like or B-like. For the remaining steps, mean values for roll and twist follow the trend set by Z p . However, the large values of standard deviations for all the steps, including B-philic steps such as AA and GA, and an A-philic step such as GG, indicate that a significant number of steps have intermediate conformation in terms of roll and twist. This is also illustrated by the Z p versus roll plot in figure 2, which does not show any clear demarcation between the A and B like steps.
The large, continuous variation in roll and twist has been observed earlier  and is also evident in the twist versus roll plot for the bound-DNA datasets (additional file 2, figure 1), where a large number of the mean roll and twist values are intermediate between those assumed by the A and B-DNA fibre models (tables 4, 5). The higher standard deviations for all the parameters in most of the steps in the bound datasets, when compared to the free DNA dataset, prompted us to individually examine the structures that are responsible for the high standard deviations. For the complex dataset, nearly all the datapoints with more than 3σ deviation from the mean roll or twist values of the free B-like DNA oligomer dataset were found to occur in structures belonging to three families-the TBP-bound DNA, the endonuclease-bound DNA and the hyperthermophile SAC7D protein-bound DNA. DNA bound to the integration host factor also undergoes significant distortions in roll and twist. For the HTH dataset, nearly all the datapoints with more than 3σ deviation from the mean roll or twist values of the free B-like DNA oligomer dataset are contributed by the TATA-box-TFIIB and CAP-bound DNA structures. On excluding these structures, the mean values are much closer to B-DNA fibre model values, with low standard deviations, and comparable to those obtained for B-like steps in the free dataset. Significantly, the exclusion of the above mentioned structural families made no significant difference in the mean values of Z p and slide for any of the steps (tables 2 and 3), indicating that the B-like DNA structure can accommodate large variations in roll and twist parameters, with no corresponding change in Z p and slide. This is further corroborated by the low correlation between either roll or twist with Z p or slide, for all the steps across the three DNA datasets. The low correlation between Z p and roll is clearly evident in figure 2.
Interestingly, several CAG triplets in the nucleosome structures  show the same unusual combination of parameters observed for the CAG triplets in some oligomers, with the CA step in BII conformation while the AG step has high roll and low twist values so that the overall roll and twist values for the two steps are similar to that in canonical B-DNA.
Free DNA oligomers can be classified as A-DNA or B-DNA in terms of Z p
At the overall structural level, most of the DNA duplexes in the free dataset can be entirely classified as A-like or B-like in terms of Z p , with the exception of 5 structures, 196D, 1P4Z, 1ZFA, 399D and 441D, wherein one or two of the steps show Z p values which differ significantly from that seen for the overall structure. Even the crystal structure of the G:C rich sequence d(CATGGGCCCATG) (1DC0), reported as an A ↔ B intermediate , assumes A-like Z p and slide values for all the steps, and hence can be described as A-like, though the roll and twist values show considerable variation.
The global x-displacement, helical rise, inclination and helical twist as well as the major and minor groove widths, described in the 'Methods' section, are also considered to be indicators of the overall A-like or B-like nature of a DNA structure, and we compared the average values of these parameters across the datasets (additional file 2, table 1). Since entire structures in the free dataset could be assigned as A-DNA or B-DNA on the basis of Z p , the averages of the global x-displacement, helical rise, inclination and helical twist for all the non-terminal basepairs within all the A-DNA structures were classified as 'A-DNA' values for the respective parameters. Similar procedure was adopted for the basepair orientation parameters within all the B-DNA structures to obtain 'B-DNA' values. As expected, the RNA dataset assumes A-like values for all the parameters, while the values for the A-like and B-like free DNA datasets being very close to their corresponding fibre model values reaffirms that the overall free DNA oligomer structures can be classified as A-like or B-like.
For both the bound DNA datasets, while the global helical rise is observed to be strongly B-DNA like, with very little variation, the global x-displacement, inclination and helical twist take up values between those for the 'A-DNA' and 'B-DNA' datasets, but closer to B-DNA. The groove width values for the bound DNA datasets for both the major and minor grooves are 'B-DNA' like. The rather large values of standard deviation for inclination and helical twist in case of the free 'B-DNA' dataset implies that B-DNA, in its free form, might be able to access the conformations observed in bound DNA.
Variations of the DNA backbone
The backbone torsion angle δ, defined by C5'-C4'-C3'-O3', the pseudorotation phase angle P , which characterises the sugar ring pucker, and the glycosidic torsion angle χ, defined by O4'-C1'-N1-C2 in pyrimidines and O4'-C1'-N9-C4 in purines, have the most characteristically distinct values in A and B-DNA [4, 5]. Figure 3 shows the variation of Z p with respect to the angle P. The two torsion angles χ (additional file 2, figure 2) and δ (additional file 2, figure 3) show similar behaviour. Note that each dinucleotide step described by a single Z p value encompasses four values of sugar pucker and glycosidic torsions, corresponding to the 4 bases constituting a dinucleotide step. As expected, the entire RNA dataset shows A-like conformation. Free DNA shows two clusters that correspond to A-like and B-like regions described by previous studies [4, 5]. An inspection of the four values of χ, δ and P that constitute each step in the free DNA dataset reveals that for a step with A-like Z p , all four values for all three angles were A-like and for a B-like step, all four values were B-like. A few exceptions were also observed in a few structures, where in a single step with A-like Z p value, one of the four P angles was observed to be B-like (1ZEY, 1ZF6, 1ZF87, 1ZFA) and vice versa (1EHV, 1DUO, 1ENN, 1IKK, 1SK5, 1ZFA, 307D, 423D, 463D, 477D, 7BNA) (see additional file 1 for detailed references corresponding to all the PDB id's). The B-like nature of backbone parameters also holds true with respect to B-like steps in the bound datasets, where all four values for all the three angles are usually B-like. Exceptions occur in the structures that displayed unusual behaviour in the local step parameters, and these have been described in detail in the relevant section.
We have also analysed the conformationally flexible torsion angles α, γ, ε-ζ, using a modified version of the algorithm of Dixit et al  such that it applied to torsion angles across a step. Table 6 and figure 4 show the distribution of the seven states described by this algorithm, across all dinucleotide steps. The RNA dataset displays classical behaviour, with an overwhelming majority of the steps assuming canonical values for α, γ, ε-ζ, viz. g-, g+, BI (state 1). For the three DNA datasets, there is much greater conformational flexibility, with α, γ, ε-ζ = g-, g+, BII (state 7), being the predominant non-canonical conformation. However, there is a significantly lower occurrence of the state 7 conformation in the bound datasets. Protein binding seems to induce a few B-like steps to assume the α, γ = t, t (state 3 or state 5) conformation, that is not preferred by free B-form DNA.
The most noteworthy difference between the free and the bound datasets was observed in the case of state 6, where the allowed ranges for α and γ occur between 0–220° and 240–270° respectively . However, all the datapoints belonging to state 6, in our four datasets, occupy a much narrower range, close to the α, γ = g+, g- conformation, that has been reported to be energetically unfavourable . While its presence is negligible in the free dataset, a substantial proportion of the steps in the HTH dataset assume this conformation (table 6). A significant proportion (22.4%) of the steps in the HTH dataset that assumed this conformation were observed to be AT, with one or both of the thymine bases in these steps taking up the unusual α and γ values. Although no restriction was placed on the ε-ζ value for state 6, almost all the datapoints for this state were observed to have the BI conformation. The steps with state 6 conformation occurred with equal frequency in the bound as well as the unbound regions of the DNA, and were not observed to assume unusual values for any other structural parameter.
Variations at the trinucleotide level
The absence of large protein-induced DNA distortion is also apparent when one examines the successive bending angles (figure 5). The successive bending angle is directly proportional to the difference in successive roll values, and can be considered to be a measure of the local bending at the trinucleotide level. The RNA dataset generally shows small successive bending angle values (with 96.9% of the values < 20°), as would be expected from a dataset with nearly uniform roll values. Of the protein-bound datasets, the HTH dataset shows surprising results. 55.7% of the triplets in the HTH dataset have bending angles between 0–10° when compared to 48.5% in the free dataset, indicating that a majority of the HTH-bound triplets tend to be less distorted than even the free triplets. 46.9% of the triplets in the complex dataset occur in this range. The trend is reversed for the range between 20–30°, which could be considered to indicate moderately 'distorted' triplets, with 14.8% of the free triplets occurring in this range when compared to only 6.8% of the HTH triplets and 10.6% of the triplets in the complex dataset. However, as noted before, binding by proteins belonging to a few specific families appears to cause large distortions in roll and twist values of a few dinucleotide steps in both the bound datasets. For example, bending angles for the steps that are distorted by TBP and CAP, in both the protein-bound datasets, range from 50° to 80°. An inspection of the stretches of DNA in the regions with high bending angles in free DNA oligomers revealed that dinucleotides with very high magnitude of roll and very large or very small twist are almost completely absent, yet a series of successive near normal roll and twist values frequently give rise to reasonably high bending angles at the triplet level. Protein binding, and especially HTH binding, does not appear to distort the DNA anymore than when it is in the free state, except in the case of a few special families.
Protein-induced distortions in DNA structure
Dinucleotide step level
Most of the significant distortions in the protein-bound DNA datasets were observed in terms of unusual roll and twist values, which occur in DNA bound to a small group of protein families. Most of these protein-DNA complexes have been exhaustively studied because of their biological relevance, leading to the perception that protein-bound DNA structure very often differs significantly from free DNA structure. The protein-bound DNA structures that are observed to be distorted can be classified into three classes: the first where the DNA structure is distorted but the distortions do not lead to strand break or strand separation, the second where the distortion leads to a nick in the DNA backbone, and third where the distortion leads to strand separation. The first class consists of DNA bound by proteins belonging to the hyperthermophile and integration host factor families in the complex dataset, and to the CAP and lac repressor families in the HTH dataset. The second class consists of DNA bound to proteins belonging to the endonuclease family in the complex dataset, and to the transposase and recombinase families in the HTH dataset. The third class consists of DNA bound to polymerases and TBPs, and occur in both the bound DNA datasets. In the following two paragraphs, the structural features of these distorted DNA are described briefly.
We have classified a step as distorted if its roll or twist value deviates by more than 3σ from the mean roll and twist values of the free B-like DNA oligomer dataset. Additional file 3 gives the base-step parameters and Z p values for the distorted steps in the DNA structures bound to different protein families. It is clear that there are a wide variety of distortions in DNA structure, depending on specific family, to which the bound protein belongs. Most of the kinks lead to significant bending of the overall structure. In the CAP-DNA complexes (additional file 4, figure 1a) and the integration host factor-DNA complex (additional file 4, figure 1b), there is a nick in the DNA backbone that seems essential for the crystallization to succeed. It is quite likely that the presence of the nick facilitates the curvature of the DNA duplex, especially since other complexes (an ARAC family transcriptional activator-DNA complex (1BL0) and a CENP-B protein-DNA complex (1HLV)) with similar modes of protein-binding as that of CAP-bound DNA, display lesser degree of bending, as estimated by end-to-end bending angle as well as d/l local values (table 7). It might also be speculated that the presence of the nick allows the protein to distort the DNA at the local level to a greater extent, thus causing a few steps to assume unusually distorted roll and twist values, though, there is no direct evidence that this occurs.
In most of these distorted DNA structures, the values of Z p and slide, as well as the backbone parameters retain classical B-DNA like values, indicating that B-DNA can accommodate considerable variations in roll and twist with little or no change in other parameters. Only a few steps in some of the endonuclease-bound DNA and the lac repressor-bound DNA have A-like Z p values. In these cases, one or two, and occasionally three of the four corresponding bases are observed to take up a C3'-endo sugar pucker. With the exception of a few endonuclease-bound DNA structures (1B94, 1B97, 1BGB, 1D02), the distorted steps themselves do not assume A-like Z p values, as evident from additional file 3. However, some of the distorted steps in the endonuclease-bound DNA have intermediate Z p values.
Gross structural level
Since the free DNA molecules are relatively short in length, it is difficult to ascertain whether the distortions observed at the local level add up to give a smooth global curvature. However both the protein-bound datasets contain several structures of greater length, hence we analysed the overall curvature of DNA structures from the complex and HTH datasets that consist of atleast 20 contiguous basepairs. We used the measures d/l local [63, 65–68], the RMSD from circle fit and the ratio of the RMSD from circle fit to the RMSD from line fit, to characterise DNA curvature, as described in the 'Methods' section. Tables 7 and 8 give the values of the ROC for long duplexes (length ≥ 18 basepairs, excluding the terminal ones) in the protein-bound datasets. For the DNA duplex that is curved, we have also estimated whether the curved helix axis is planar or has an out-of-plane component. A segment of the nucleosome structure (1KX5 ) consisting of 30 basepairs gives a d/l local value of 0.77 and a ROC value of 37.9 Å, with circle fit standard deviation of 0.6 Å, indicating that the measures used are quite adequate to define curvature of DNA molecules of this length.
Of the 28 HTH-bound DNA structures with length ≥ 20 basepairs (table 7), 7 meet the RMSD criteria that allows for reliable geometry assignment. Of these, 5 were found to have d/l local ≤ 0.98. A further 21 duplexes that did not meet the RMSD criteria were classified as unassigned. Of the structures where the DNA is curved, the CAP-DNA complexes (1CGP, 1J59, 1RUN), an ARAC family transcriptional activator-DNA complex (1BL0) and a CENP-B protein-DNA complex (1HLV), all consist of a dimeric protein that binds to two successive major grooves of the DNA, approximately one helix turn apart and the DNA is essentially curved due to two major in-phase kinks (additional file 4, figure 1a). All these duplexes have a negative out-of-plane bending angle.
Among the 21 DNA duplexes that have been classified as unassigned, more complex types of protein-binding is observed, indicating that there are several different modes of curvature for a DNA bound to the HTH motif. In some of these structures (CRE recombinase protein (4CRX), γδ resolvase-DNA complex (1GDT)), bending appears to be a result of large kinks at one or two steps in the duplex, as is evident from the higher maximum bending angle values obtained for these structures.
Among the HTH-bound DNA duplexes that are found to be curved, only the CAP binding duplexes have large values for the average successive bending angle, while for all the other duplexes, it is < 12°, indicating that curvature of the duplex can arise due to the cumulative effect of small amount of bending along the entire helix. In contrast, some of the unassigned duplexes that have d/l local > 0.98, and so could be considered linear, were observed to have large average successive bending angle values (lambda CI-NTD – sigma-region4 – DNA complex (1RIO, shown in additional file 4, figure 1c), PAX5-DNA complex (6PAX)), indicating that large distortions at the local level can cancel each other out and so need not cause the entire duplex to bend significantly.
In the complex dataset, there are only 9 structures (excluding the reference nucleosome structure) with DNA length ≥ 20 basepairs (table 7). Of these, 4 were found to be curved by the criteria of d/l local and ROC. The curved duplexes all comprise of DNA bound by I-Cre I endonucleases that bind to and have interactions along the entire length of the duplex. Among the unassigned duplexes, the DNA assumes a U-shaped structure in the integration host factor-DNA complex 1OWF (additional file 4, figure 1b), wherein a nick has been introduced into the DNA backbone to facilitate crystallisation. The duplex upto the nick has a d/l local value of 0.81 (given in table 8), but the entire 34 basepair duplex has a d/l local value of 0.33. In two endonuclease-bound complexes 1A73 and 1CYQ, the presence of junctions between two A-like regions separated by a B-like region leads to a non-linear geometry, as indicated by the d/l local values (~0.90), but the duplexes are not smoothly curved. A lambda-CII-DNA complex (1ZS4), and a human TBOX-protein 3-DNA complex (1H6F), are linear as per d/l local values. Some representative examples, showing the 3-dimensional path of the basepair centres in the DNA duplex, with different amounts of curvature are shown in figure 6, while the cartoon diagrams of few DNA-protein crystal structure complexes are shown in additional file 4, figure 1. They clearly illustrate the different extent of curvature (or lack of it), adopted by protein-bound DNA molecules.
No correlation was observed between the curvature and occurrence of the various backbone geometries (even the energetically unfavourable state 6 conformation) in these structures.
From our analysis, it is clear that for the individual dinucleotide steps in the free oligomer dataset, the dinucleotide step parameters Z p and to a lesser extent, slide, as well as the pseudorotation phase angle P for the sugar ring, the backbone torsion angle χ and the glycosidic torsion angle δ for the individual bases in a step are better indicators of A-like or B-like conformation than the traditionally used parameters of roll and twist, confirming the findings by earlier studies [14, 27]. A few of the dinucleotide steps seem to have a distinct preference for a particular conformation-AA and GA steps are strongly B-philic, while only the GG step is strongly A-philic.
In an earlier study to characterize how the DNA sequence defines conformation, Hays et al  have reported crystal structures of all the permutations of the inverted repeat d(CCnnnN6N78GG) under well-defined crystallographic conditions, which take up A-form, B-form and Holliday junction structures. Several of the structures reported in that work fit our selection criteria and are also part of the free dataset in the present study. The authors observed that the set of A-DNA crystal structures reported in their study are conformationally more uniform than the B-DNA structures. This also seems to be the case for the larger and sequentially more heterogeneous dataset analysed in this study vis-à-vis Z p , slide, roll and twist, not only for the entire free dataset, but also for the individual basepair dinucleotide steps. With the exception of the slide parameter for AT step, and the roll parameter for GG step, the standard deviation is always higher for the B-like free DNA steps as against the A-like steps, for all the four parameters. However, for a step like CA, one must consider the fact that it is in fact trimodal, with the B-like steps further subdivided into BI and BII conformations, making it difficult to compare the variation between A-like and B-like conformations. It must also be noted that with the exception of slide in case of A-like steps, the standard deviations obtained in our study are lower compared to those observed by Hays et al  for the slide, roll and twist parameters examined in both the studies. Hays et al  observed that the trinucleotide motifs GGN, NGG and CC(C/G) favour a transition to A-DNA conformation. Our analysis supports this conclusion, with most of the GC rich structures taking up an A-DNA conformation. The only exception to this rule was observed to be 1ZFB, for the sequence d(CCGCCGGCGG). However, an earlier structure (382D ) of the same sequence (not included in this study owing to resolution cutoff criteria), was observed to have an A-like conformation. Hence A-DNA definitely seems to be favoured by GC rich DNA, especially those with oligo-G tracts. Since the GG step is observed to be the most A-philic, this is to be expected. However a GC rich sequence which does not have an oligo-G tract does not necessarily favour an A-DNA conformation, since GC and CG steps seem to assume A-like and B-like conformation with nearly equal ease. In addition, with the exception of AC, which seems to equally favour the A and B-form, and CA, which is trimodal, all other steps where one or both basepairs are A:T seem B-philic in the free dataset.
For B-DNA structures, Gorin et al  have correlated the extent of B-DNA twisting with the basepair morphology and clash between the exo-cyclic groups in the four bases. The average values for slide, roll and twist, obtained by Gorin et al  for the dinucleotide steps in their dataset, comprising of B-DNA structures with a resolution cutoff of 3.0Å, are quite similar to the 'B-like' average values for different dinucleotide steps in the present study. (tables 3, 4, 5). The overall average values in the two studies are also observed to be similar. However, interestingly the low twist CA (BI) and high twist TA steps in Gorin et al dataset converge to nearly similar values in our high resolution dataset (36.5° and 34.7° respectively). This positions the CA step in a favoured conformation with minimal clash as predicted by the clash strength function designed by these authors, while the TA step is positioned in a less favourable conformation.
How different is bound DNA from free B-DNA?
In case of bound DNA, the DNA duplexes are almost entirely B-like in conformation in terms of Z p and slide, while roll and twist predominantly show variation that is similar to that of free DNA. The average values for slide, roll and twist, obtained in an earlier study by Olson et al  for different dinucleotide steps in a dataset of protein-DNA crystal structure complexes, are comparable to those obtained in this study for the complex dataset (excl. TEH, tables 3, 4, 5) and the HTH dataset (excl. TC, tables 3, 4, 5). This is expected, since these authors also considered only the step parameter values within 3σ deviation of the mean, for their 'B-like' protein-DNA dataset, essentially excluding the distorted steps. The overall average values for slide, roll and twist reported for the 'B-like' protein-DNA dataset are also observed to be similar to those obtained for the complex (excl. TEH) dataset and the HTH (excl. TC) dataset in this study.
Only a few DNA structures, bound to proteins belonging to a small group of families have highly unusual structural parameters, principally roll and twist. Apart from those structures, other structures in both the protein-bound datasets principally take up free B-DNA like values for all the dinucleotide step parameters. The nucleosome structure , considered a classic case of a highly curved structure, does not have highly unusual Z p , slide, roll or twist values, with only 5 out of 146 roll values lying just outside the 3σ deviation range of the B-DNA like oligomer dataset. The lack of sharp kinks gives the nucleosome structure a smooth curvature, with very small out-of-plane component for ~30 basepair fragments. For any randomly selected 30 basepair fragment of this structure, the RMSD from a plane fit was always < 0.15 Å, as against a value of 0.65 Å for the highly distorted CAP-bound DNA structure 1J59. Even for a randomly selected fragment of 76 basepairs (which nearly completes a full circle), the RMSD from a plane fit was observed to be only 0.26 Å. Similarly, the out-of-plane torsion angle values, for several random fragments of ~30 basepairs, were always observed to be < 10°, another indicator of the smooth curvature and gentle, regular pitch of the superhelix. It is also interesting to note that the ROC calculated for the 76 basepair fragments in 1KX3 and 1KX5 are 39.8 Å and 39.4 Å respectively, while that for a 2.8 Å resolution structure (1AOI ) is calculated as 41.5 Å, indicating that the DNA in different nucleosome structures has small variations in curvature.
A-like steps are limited to DNA bound to proteins from a few specific families
Protein-bound DNA structures, apart from being perceived as distorted, have also been characterised as being predominantly A-like . Our analysis clearly refutes this characterisation. The only protein-bound DNA structures having few steps with A-like values of Z p are those bound to some of the endonucleases, DNA polymerases, transposases and the homeodomains. Of the endonuclease bound-DNA structures, the IPpo I endonuclease-bound DNA structures 1A73 and 1CYQ have two separate A-like half turns , leading to a non-linear geometry, as explained above. The PvuII endonuclease-bound DNA structure 3PVI has an entire A-like stretch with only few bases at one end having a B-like geometry in terms of Z p .
Polymerase-bound DNA (1L3S, 1L3U, 1L3V, 1L5U, 1L3T, 1LV5, 1NJY, 1NJZ, 1NK5, 1NK6, 1NK7, 1NKC, 1U45, 1UA1, 2BDP, 4BDP) undergo gradual transition in Z p and slide, with B-like values at one end to intermediate to A-like values towards the nucleotide incorporation end, but do not show large variation in their B-like roll and twist values. At least two, and in several instances three, of the values for P, δ and χ are A-like. The DNA duplexes in these complexes have a variety of sequences indicating that the A-like nature of these duplexes is not sequence dependent, but is a result of protein-binding and the polymerisation process itself.
Transposase-bound DNA (1TC3, 1U78) have a 4–7 basepair long G:C rich region at one terminal of the duplex that assumes A-DNA like conformation in terms of all the parameters. The other end of the DNA duplex is A:T rich, with a narrow minor groove  and curvature characteristic of DNA containing oligo A-tracts [58, 59]. The transposase protein binds to these two regions-the G:C rich A-like region and the A:T rich region, via two HTH motifs that are connected by a long linker . However, the features observed for the G:C rich region and the oligo A-tract in the transposase-bound DNA are similar to those observed for free G:C rich oligomers and free oligo A-tracts, hence these features can be said to be intrinsic to the DNA sequence.
There are 18 homeodomain-bound DNA structures in the protein-bound datasets. Of these, only five steps (3 GG and 2 CG steps), occurring in four different structures, are observed to take up A-like values of Z p .
Another class of DNA, often cited as an example of protein-induced B ↔ A transition, are the Zn-finger-bound DNA structures. Nekludova et al  have shown that for a variety of protein-bound DNA molecules, including Zn-finger-bound DNA, a distinctive conformation with an enlarged major groove when compared to B-DNA, was observed. In our study, all the Zn-finger-bound DNA structures (1A1G, 1A1H, 1A1I, 1A1J, 1A1K, 1A1L, 1AAY, 1G2F, 1JK1, 1JK2, 1LLM) (see additional file 1 for the references for all these structures) assume B-like values of Z p for all the steps. Similarly the backbone sugar pucker and the torsion angle χ are observed to be close to B-DNA fibre model values. Several other parameters, indicating the A-like or B-like nature of a DNA structure, also assume B-like or intermediate values. The average global helical rise (3.3 ± 0.1Å) is entirely B-like. The slide values are B-like to intermediate, whereas the average values of global x-displacement (-1.4 ± 0.5 Å), inclination (7.1 ± 4.7°) and helical twist (33.5 ± 4.9°) are intermediate but closer to B-DNA fibre model values. The average interstrand P-P distance across the major groove is 18.0 ± 1.6 Å, very close to the B-DNA fibre model value while the average interstrand P-P distance across the minor groove is 12.6 ± 1.3 Å, again only marginally larger than the B-DNA fibre model value.
It is also interesting that while the 7 nucleotide long runt-domain binding site DNA sequences in the free form are reported to assume A-DNA like conformation (1XJX, 1XJY)  as well as near B-DNA like conformation [34, 93], depending on the flanking bases, the same sequences bound to the runt-domain protein are found to assume B-like values of Z p , slide, the backbone torsion angles and the groove widths. X-displacement, inclination and helical twist take up intermediate, but closer to B-DNA values – a behaviour similar to that observed for the bound DNA datasets in this study. Thus, in this case, while the free DNA sequences assume both the A and B-forms, the protein-bound DNA take up the B-form.
It has been suggested that the TBP-bound DNA conformation is closer to an A-DNA and the inherent A-philicity of the TATA sequence might facilitate the transition to the near A-like bound-conformation . Our analysis indicates that these assertions are not always valid. For example, the oligomer structure 1VJ4  for the sequence d(GGTATACC), takes up an A-DNA conformation, but the free DNA structures 1D56 and 1D57 , for the decamer d(CGATATATCG), both take up a B-DNA conformation, despite encompassing the TATA stretch. The TATA stretch in the TBP-bound DNA structures also take up entirely B-like Z p values. Though some of the other parameters such as roll, twist and rise do not have classical B-DNA values, this is more indicative of a distortion from the B-form, but not necessarily to an A-like conformation. The B-like nature of the TBP-bound DNA in terms of Z p and slide is also observed for a couple of hexamer sequences, which occur in both the free DNA dataset and some of the TBP-bound DNA structures. While the sequence TTTAAA takes up B-like Z p and slide values in the free (1IKK , 1SK5) as well as the TBP-bound DNA (1D3U , 1QNA ), the hexamer stretch GGCGCC takes up an A-DNA conformation in the free DNA structure 414D  as expected, but is observed to take up B-like Z p and slide values in the TBP-TFIIB-DNA complex 1C9B . It is also noteworthy that unlike the TBP-bound DNA from the complex dataset, the backbone parameters P and δ as well as χ take up entirely B-like values for the TBP-TFIIB-bound DNA from the HTH dataset, consistent with our observation throughout this study that HTH-bound DNA tends to be more B-like than other protein-bound DNA molecules.
There have been studies of protein-DNA complexes, using backbone conformational parameters such as sugar pucker  or the χ and δ torsion angles  to classify the DNA nucleotides as A-like or B-like. The Z p versus sugar pucker (figure 3) plot as well as the Z p versus δ (additional file 2, figure 2) and Z p versus χ (additional file 2, figure 3) plots clearly indicate that a C3'-endo conformation or A-like values of χ or δ do not necessarily imply an A-like conformation in the protein-bound datasets. Significantly, Tolstorukov et al  find only 12% of the protein-interacting nucleotides with a C3'-endo sugar pucker conformation. On the other hand, Lejeune et al  conclude that "A-DNA is more frequently implicated in protein-DNA interactions than the classical B-DNA conformation". We do not find this claim to be valid, using any of the backbone parameters for 'A versus B' discrimination.
HTH-bound DNA, while remaining B-like, occassionally takes up an unfavourable backbone conformation
The only effect that can be unambiguously ascribed to protein binding in the predominantly B-DNA like protein-bound duplexes occurs in the DNA backbone. The DNA backbone in the free dataset is quite uniform, with the angles α and γ almost completely in the canonical g-, g+ conformation in B-DNA, and (α, γ, ε-ζ) ranging from (g-, g+, BI) to (t, t, BI) conformation in A-DNA. On the other hand, backbone torsion angles in protein-bound DNA are observed to be considerably distorted. Steps that are B-like in terms of Z p and slide are observed to assume a wide variety of backbone conformations that are highly unusual, and in some cases, energetically unfavourable. In particular, HTH binding causes α and γ angles in DNA to assume the energetically unfavourable g+, g- conformation in much higher proportion (11.7%) than in unbound DNA. As described in the 'Results' section, the steps taking up this energetically unfavourable conformation occur with equal frequency in the bound as well as the unbound regions of the DNA, and are not observed to assume unusual values for any other structural parameter. Overall, 57 out of the 97 HTH-motif bound DNA structures are observed to adopt this unfavourable backbone conformation for some of the steps. Of these, 24 structures have 5 or more occurrences of the unfavourable backbone conformation. Thus it is seen that there are a large number of structures with atleast a few steps in this conformation. These structures have been solved in a variety of space groups. The proteins binding to the 24 HTH-bound DNA structures with 5 or more occurrences of state 6 get classified into 15 different SCOP classes. Thus it appears that binding by the HTH motif allows the DNA backbone to assume this energetically unfavourable conformation, even when there is no direct contact between them.
At the tri-nucleotide level, bound DNA, and especially HTH-bound DNA appears to have less distortion than free DNA. At the gross structural level, nearly half of the DNA structures of length ≥ 20 basepairs and bound to the HTH motif were observed to have moderate curvature. It was observed that in several of these cases, the DNA was bound by a dimer of 2 HTH motifs, with the two monomers binding to DNA at regions one helix turn apart and bending it in the same direction so that there was a net overall curvature. However there are other modes by which the DNA bound to the HTH motif was observed to be curved, such as in the case of the MAT alpha2-MCM1-DNA ternary complex 1MNM and the CRE-recombinase-DNA complex (4CRX). Yet other modes of curvature of protein-bound DNA are revealed in the complex dataset. Thus it is not possible to determine a uniform mode and mechanism for the DNA curvature observed in the bound datasets. With the exception of a few structures where it was difficult to determine three uniformly undistorted regions separated by a large kink, all the curved DNA structures have a negative out-of-plane component. With no long free DNA oligomers in the dataset, it is difficult to conclude whether free DNA by itself can attain such conformations and the protein merely 'locks' it in that conformation or the protein actually bends it to that state. Most of these curved structures, however, do not have highly unusual step parameter values and hence it is possible that longer free DNA oligomers with similar sequences might be able to achieve such curved conformations without the aid of proteins. Even a few steps with unusual parameters might occur in long free DNA oligomers, as indicated by the spontaneous development of one or two sharp kinks in the molecular dynamics simulations of 94 basepair free DNA minicircles . This has interesting implications especially for the HTH-binding DNA, since a majority of the proteins in this dataset are transcription activators or repressors, whose function on binding to the DNA is to cause structural changes in the DNA that allow or prevent other proteins of the transcription machinery to bind to the DNA and carry out transcription. It is tempting to speculate that these proteins merely increase the 'lifetime' of those conformations, as against inducing unfavourable conformations, which involves a much higher energetic cost. However, this needs to be verified using experimental and theoretical methods that trace the dynamic evolution of DNA structures under different conditions.
The free DNA oligomers, even in the crystalline state, sample a large conformational space, but each molecule is found to be entirely in the A or B form, depending primarily on its sequence. In case of protein-bound DNA, the claim that protein-binding generally favours the A-form of DNA , as well as the perception that it induces an energetically unfavourable conformation, are invalid. We find that the role of A-form is limited to the DNA structures bound to a few specific protein families such as transposases and DNA polymerases. Protein-induced distortion in DNA can occur via one of several different modes, such as a few steps taking up high positive roll and a smaller twist, a BII like transition of the backbone, leading to a negative roll and large twist, or in some cases, the two strands in the helix being pulled apart. However, these large, induced deviations from the free B-form are observed only in the DNA structures bound to the proteins such as CAP, TBP, integration host factor and Cre recombinase. It is to be noted that, even in these structures, the distortions are limited to a few steps and the remainder of the duplex shows B-DNA like features. In a large number of cases of the HTH motif-bound DNA, protein-binding does not induce any distortion in the dinucleotide step geometry, but the duplex takes up an energetically unfavourable backbone conformation, even when there are no contacts between the protein and the DNA backbone. Barring these exceptions, the average parameters at the level of dinucleotide step, trinucleotide and the backbone of protein-bound DNA structures, across a large and diverse set of protein families, are quite close to the free B-DNA oligomer values. Interestingly, this is observed even though very few hexamer or longer sequence motifs are common to the free and bound datasets, and the free DNA dataset is significantly smaller than the bound DNA datasets in terms of size. It is also striking to note that even a duplex structure as far away from a 'straight' DNA as seen in the 147 basepair long nucleosome, has very few (≤ 5) steps with highly distorted local parameters, indicating that 'normal' B-like parameters at the local level can cumulatively give rise to double helical structures with a wide range of geometries. These observations highlight the amazing adaptability of this structural form, and may explain why it has evolved to be biologically the most relevant design for double-helical DNA.
Crystallographic dataset generation
The four X-ray crystallographic datasets used in the analysis are (i) RNA oligomers dataset (hereafter referred to as the RNA dataset), (ii) DNA oligomer dataset (hereafter referred to as the free dataset), (iii) DNA-protein complexes dataset (excluding DNA bound by the HTH protein) (hereafter referred to as the complex dataset), and (iii) DNA-HTH protein complexes dataset (hereafter referred to as the HTH dataset). The RNA, free and complex datasets were extracted from the Protein Data Bank (PDB) . All three datasets contain structures with a resolution of 2.0 Å or better. Structures in the PDB that have the DNA-binding HTH motif were identified using the tool PredictDNAHTH, developed by McLaughlin et al . Since only 33 DNA-HTH protein complexes with a resolution of 2.0 Å or better were identified, the resolution cut-off for the HTH dataset was increased to 3.0 Å. There was no significant difference between the results obtained for the dataset with a cutoff of 2.0 Å and the dataset with a cutoff of 3.0 Å. Therefore the larger dataset, with a cutoff of 3.0 Å was used. In the three DNA datasets, only fragments of the DNA consisting of atleast 8 contiguous Watson-Crick basepairs were considered. The RNA dataset had much shorter structures, hence the length cutoff was reduced to five contiguous basepairs. Also steps with non-Watson-Crick basepairs, present in significant numbers in the RNA dataset, were not included in this analysis. In the free dataset, structures with any ligands other than ions or water were excluded. Identical basepairs from structures with a two-fold symmetry were considered only once. The RNA dataset consists of 52 structures (additional file 1) (75 individual duplexes which contain 276 dinucleotide steps comprising of Watson-Crick base pairs). The free dataset consists of 76 structures (additional file 1) (77 individual duplexes which contain 406 basepaired dinucleotide steps comprising of Watson-Crick basepairs).
The complex dataset consists of 85 structures (additional file 1) (112 duplexes which contain 1227 dinucleotide steps comprising of Watson-Crick basepairs). The HTH dataset (at ≤ 3.0 Å) consists of 97 structures (additional file 1) (126 duplexes which contain 1559 dinucleotide steps comprising of Watson-Crick base pairs).
Evaluation of dinucleotide step parameters and global helical parameters
The structural parameters of the duplexes i.e. the basepair parameters propeller twist, buckle, opening angle, shear, stretch and stagger as well as the dinucleotide step parameters tilt, roll, twist, shift, slide, and rise were determined by the NUPARM program [103–105], for all the four datasets. The parameter Z p , defined as the mean z-coordinate of the backbone phosphate atoms of the basepair with respect to the basepair dimer reference frame, was also calculated using the revised NUPARM program .
The dinucleotide step parameters tilt, roll and twist measure the relative rotational motion between adjacent basepairs about the x, y and z-axis respectively of a local basepair doublet coordinate system, whereas the dinucleotide step parameters shift, slide and rise measure relative translational motion between adjacent basepairs along the local doublet x, y and z-directions respectively.
The global helical parameters viz. the rotational parameters inclination, tip and the helical twist and the translational parameters x-displacement, y-displacement and z-displacement were also calculated using the NUPARM program. Inclination denotes the rotation of the basepair about the x-axis, tip denotes rotation about the y-axis and helical twist denotes rotation about the helical axis. Similarly the translational parameters denote displacement along the three axes. The mean of the global x-displacement, helical rise, inclination and helical twist for all the non-terminal basepairs within all the structures in a dataset were classified as the average values for the respective dataset. The protein-bound DNA sructures in which the roll or twist value for atleast one step deviated by more than 3σ from the mean roll and twist values of the free B-like DNA oligomer dataset, and also those structures which were curved or whose geometry of curvature could not be assigned (as given in tables 7, 8), were excluded from the calculation of mean values of global helical parameters, since fitting a single linear helical axis would be untenable in these cases. Overall, 49 structures from the complex dataset and 62 structures from the HTH dataset were included for these calculations.
Evaluation of groove widths
The minor groove width and the major groove width were calculated as the smallest interstrand phosphate separations along the two grooves, using the NUPARM program. Please note that the groove widths as defined here also include the phosphate diameter value.
Calculation and classification of backbone torsion angles
Backbone torsion angles α, β, γ, δ, ε, ζ; the glycosidic torsion angle χ and the pseudo rotation angle P  were calculated using the NUPARM program. Backbone torsion angles for a basepaired dinucleotide step, i.e., across the phosphodiester bond, were clustered and analysed. ε: C4' n -C3' n -O3' n -P n , ζ: C3' n -O3' n -Pn+1-O5'n+1, α: O3' n -Pn+1-O5'n+1-C5'n+1and γ: O5'n+1-C5'n+1-C4'n+1-C3'n+1are classified into seven states as per the algorithm proposed by Dixit et al . Since ε and ζ assume two related conformations, ε, ζ = t, g- being the canonical conformation, known as BI, and ε, ζ = g-, t being the non-canonical conformation, known as BII, a value of ε - ζ ≤ 0 has been classified as the BI conformation, and a value of ε - ζ > 0 has been classified as the BII conformation.
Calculation of bending/curvature
The calculation of the radius of curvature using a least square circle fit method and the ratio of end-to-end distance to the contour length (d/l local or d/l max ) were done as described previously in . The measure d/l max is reasonably independent of the length of the DNA sequence (data not shown), except for highly curved long DNA molecules, as in nucleosomal DNA, but does not distinguish between different types of bending for sequences with fewer than 30 basepairs (data not shown). The radius of curvature (ROC) is calculated by fitting a circle to the basepair centres of the DNA molecules. Smaller the radius of this circle, the more curved the DNA is. However, the quality of the fit to a circle is affected to a large extent by distortions at the local level in the duplex i. e. the successive bending angles. Thus the presence of several triplets that are distorted, even to a small degree, would lead to a poor circle fit and consequently an inaccurate value of radius of curvature (ROC). Thus ROC is only reported when the RMSD for a circle fit is ≤ 1.0 Å, and the ratio of RMSD for a circle fit to that for a line fit is ≤ 0.6.
When the d/l local value is ≤ 0.98, the RMSD for a circle fit is ≤ 1.0 Å, and the ratio of RMSD for a circle fit to that for a line fit is ≤ 0.6, we have assigned the DNA molecule geometry to be curved. When the d/l local value is > 0.98, the RMSD for a line fit is ≤ 1.0 Å, and the ratio of RMSD for a circle fit to that for a line fit is > 1.6, we have assigned the DNA molecule geometry to be linear. When neither the 'curved' nor 'linear' criteria are satisfied, the geometry of the DNA duplex is considered as 'unassigned'. For the DNA duplex that is curved, the out-of-plane component of DNA curvature was calculated as the torsion angle between the global helix axes vectors fitted to three relatively straight sections of the DNA molecule, separated by large kinks.
A local helix axis vector corresponding to each dinucleotide step is defined as the vector pointing in the direction of the cross-product of the differences of the x and y-vectors of the constituent basepair planes. The angle between two local helix axes vectors corresponding to overlapping dinucleotide steps, described as the successive bending angle, as well as the angle between the vectors corresponding to the dinucleotide steps at the two ends of the molecule, and described as the end-to-end bending angle, were also calculated using NUPARM and used as measures of curvature.
The entire analysis of the dinucleotide step parameters, backbone torsion angle parameters, the successive bending angles, the radius of curvature, d/l local and out-of-plane components of DNA curvature has been carried out excluding the terminal basepairs to eliminate end effects. The end-to-end bending angle has also been measured as the angle between the local helix axes vectors corresponding to the penultimate dinucleotide steps.
All the plots were generated using the MATLAB-7.4 package.
The values of the basepair parameters, base-step parameters as well as the backbone torsion angles obtained using the NUPARM package were compared to those obtained by the X3DNA package . The general trend of the parameters was observed to be similar. The parameters calculated by the two programs were different for the distorted regions of a few protein-bound DNA structures.
Saenger W: Principles of nucleic acid structure. New York: Springer; 1984.
Ghosh A, Bansal M: A glossary of DNA structures from A to Z. Acta Cryst 2003, 59: 620–626. 10.1107/S0108270103019577
Conner BN, Yun C, Dickerson JL, Dickerson RE: Helix geometry and hydration in an A-DNA tetramer: IC-C-G-G. J Mol Biol 1984, 174: 663–695. 10.1016/0022-2836(84)90089-5
Dickerson RE: Base sequence and helix structure variation in B and A DNA. J Mol Biol 1983, 166: 419–441. 10.1016/S0022-2836(83)80093-X
Dickerson RE, Kopka ML, Pjura P: J Mol Biol. Volume 2. Wiley, NewYork; 1984.
Schneider B, Neidle S, Berman HM: Conformations of the sugar-phosphate backbone in helical DNA crystal structures. Biopolymers 1997, 42: 113–124. 10.1002/(SICI)1097-0282(199707)42:1<113::AID-BIP10>3.0.CO;2-O
Bhattacharyya D, Bansal M: Groove width and depth of B-DNA structures depend on local variation in slide. J Biomol Struct Dyn 1992, 10: 213–226.
Drew HR, Dickerson RE: Structure of a B-DNA dodecamer. III. Geometry of hydration. J Mol Biol 1981, 151: 535–556. 10.1016/0022-2836(81)90009-7
Fratini AV, Kopka ML, Drew HR, Dickerson RE: Reversible bending and helix geometry in a B-DNA dodecamer: CGCGAATTBrCGCG. J Biol Chem 1982, 257: 14686–14707.
Bhattacharyya D, Bansal M: Local variability and base sequence effects in DNA crystal structures. J Biomol Struct Dyn 1990, 8: 539–572.
Gorin A, Zhurkin VB, Olson WK: B-DNA twisting correlates with base-pair morphology. J Mol Biol 1995, 247: 34–48. 10.1006/jmbi.1994.0120
Olson WK, Bansal M, Burley SK, Dickerson RE, Gerstein M, Harvey SC, Heinemann U, Lu XJ, Neidle S, Shakked Z, Sklenar H, Suzuki M, Tung CS, Westhof E, Wolberger C, Berman HM: A standard reference frame for the description of nucleic acid base-pair geometry. J Mol Biol 2001, 313: 229–237. 10.1006/jmbi.2001.4987
McCall M, Brown T, Kennard O: The crystal structure of d(G-G-G-G-C-C-C-C). a model for poly(dG).poly(dC). J Mol Biol 1985, 183: 385–396. 10.1016/0022-2836(85)90009-9
Lu X, Shakked Z, Olson WK: A-form conformational motifs in ligand-bound DNA structures. J Mol Biol 2000, 300: 819–840. 10.1006/jmbi.2000.3690
El Hassan MA, Calladine CR: Conformational characteristics of DNA: empirical classifcations and a hypothesis for the conformational behaviour of dinucleotide steps. Phil Trans R Soc Lond A 1997, 355: 43–100. 10.1098/rsta.1997.0002
Hays F, Teegarden A, Jones ZJ, Harms M, Raup D, Watson J, Cavaliere E, Ho P: How sequence defines structure: a crystallographic map of DNA structure and conformation. Proc Natl Acad Sci 2005, 102: 7157–7162. 10.1073/pnas.0409455102
Dickerson RE, Ng HL: DNA structure from A to B. Proc Natl Acad Sci 2001, 98: 6986–6988. 10.1073/pnas.141238898
Ng HL, Dickerson RE: Mildly eccentric 'E-DNA'. Nat Struct Biol 2001, 8: 107–108. 10.1038/84078
Seeman NC, Rosenberg JM, Suddath FL, Kim JJ, Rich A: RNA double-helical fragments at atomic resolution. I. the crystal and molecular structure of sodium adenylyl-3',5'-uridine hexahydrate. J Mol Biol 1976, 104: 109–144. 10.1016/0022-2836(76)90005-X
Rosenberg JM, Seeman NC, Day RO, Rich A: RNA double-helical fragments at atomic resolution. II. the crystal structure of sodium guanylyl-3',5'-cytidine nonahydrate. J Mol Biol 1976, 104: 145–167. 10.1016/0022-2836(76)90006-1
Hingerty B, Subramanian E, Stellman SD, Sato T, Broyde SB, Langridge R: The crystal and molecular structure of a calcium salt of guanylyl-3',5'-cytidine (GpC). Acta Cryst 1976, B32: 2998–3013.
Aggarwal A, Islam SA, Kuroda R, Sanderson MR, Neidle S, Berman HM: The structure of the ribodinucleoside monophosphate guanylyl-3',5'-cytidine as its ammonium octahydrate salt. Acta Cryst 1983, B39: 98–104.
Dock-Bregeon AC, Chevrier B, Podjarnyl A, Johnson J, de Bear JS, Gough GR, Gilham PT, Moras D: High resolution structure of the RNA duplex [U(U-A) 6 A] 2 . Nature 1988, 335: 375–378. 10.1038/335375a0
Dock-Bregeon AC, Chevrier B, Podjarnyl A, Johnson J, de Bear JS, Gough GR, Gilham PT, Moras D: Crystallographic structure of an RNA helix: [U(UA) 6 A] 2 . J Mol Biol 1989, 209: 459–474. 10.1016/0022-2836(89)90010-7
Olson WK, Gorin AA, Lu X, Hock LM, Zhurkin VB: DNA sequence-dependent deformabilty deduced from protein-DNA crystal complexes. Proc Natl Acad Sci 1998, 95: 11163–11168. 10.1073/pnas.95.19.11163
Jones S, van Heyningen P, Berman HM, Thornton JM: Protein-DNA interactions: a structural analysis. J Mol Biol 1999, 287: 877–896. 10.1006/jmbi.1999.2659
Lu X, Olson WK: 3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures. Nucleic Acids Res 2003, 31: 5108–5121. 10.1093/nar/gkg680
Ivanov VI, Minchenkova LE, Minyat EE, Frank-Kamenetskii MD, Schyolkina AK: The B to A transition of DNA in solution. J Mol Biol 1974, 87: 817–833. 10.1016/0022-2836(74)90086-2
Zavriev SK, Minchenkova LE, Frank-Kamenetskii MD, Ivanov VI: On the flexibility of the boundaries between the A-form and B-form sections in DNA molecule. Nucleic Acids Res 1978, 5: 2657–2663. 10.1093/nar/5.7.2657
Minchenkova LE, Schyolkina AK, Chernov BK, Ivanov VI: CC/CG contacts facilitate the B to A transition of DNA in solution. J Biomol Struct Dyn 1986, 4: 463–476.
Peticolas WL, Y Wang Y, Thomas GA: Some rules for predicting the base-sequence dependence of DNA conformation. Proc Natl Acad Sci 1988, 85: 2579–2583. 10.1073/pnas.85.8.2579
Tolstorukov MY, Ivanov VI, Malenkov GG, Jernigan RL, Zhurkin VB: Sequence-dependent B↔A transition in DNA evaluated with dimeric and trimeric scales. Biophys J 2001, 81: 3409–3421. 10.1016/S0006-3495(01)75973-5
Shakked Z, Rabinovich D, Cruse WB, Egert E, Kennard O, Sala G, Salisbury SA, Viswamitra MA: Crystalline A-DNA: the X-ray analysis of the fragment d(G-G-T-A-T-A-C-C). Proc R Soc Lond 1981, B213: 479–487. 10.1098/rspb.1981.0076
Kitayner M, H R, Rabinovich D, Shakked Z: Structures of the DNA-binding site of Runt-domain transcription regulators. Acta Crystallogr D Biol Crystallogr 2005, 61: 236–246. 10.1107/S0907444904032378
Berman H, Olson W: Molecular biology in the new era. Faircount, LLC, London; 2003.
Hagerman PJ: Pyrimidine 5-methyl groups influence the magnitude of DNA curvature. Biochemistry 1990, 29: 1980–1983. 10.1021/bi00460a003
Suzuki M, Yagi N: Stereochemical basis of DNA bending by transcription factors. Nucleic Acids Res 1995, 23: 2083–2091. 10.1093/nar/23.12.2083
Young MA, Ravishanker G, Beveridge DL, Berman HM: Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes. Biophys J 1995, 68: 2454–2468. 10.1016/S0006-3495(95)80427-3
Guzikevich-Guerstein G, Shakked Z: A novel form of the DNA double helix imposed on the TATA-box by the TATA-binding protein. Nat Struct Biol 1996, 3: 32–37. 10.1038/nsb0196-32
Dickerson RE, Chiu TK: Helix bending as a factor in protein/DNA recognition. Biopolymers 1997, 44: 361–403. 10.1002/(SICI)1097-0282(1997)44:4<361::AID-BIP4>3.0.CO;2-X
Dickerson RE: DNA bending: the prevalence of kinkiness and the virtues of normality. Nucleic Acids Res 1998, 26: 1906–1926. 10.1093/nar/26.8.1906
Mandel-Gutfreund Y, Schueler O, Margalit H: Comprehensive analysis of hydrogen bonds in regulatory protein DNA-complexes: in search of common principles. J Mol Biol 1995, 253: 370–382. 10.1006/jmbi.1995.0559
Otwinowski Z, Schevitz RW, Zhang RG, Lawson CL, Joachimiak A, Marmorstein RQ, Luisi BF, Sigler PB: Crystal structure of trp repressor/operator complex at atomic resolution. Nature 1988, 335: 321–329. 10.1038/335321a0
Beamer LJ, Pabo CO: Refined 1.8 Å crystal structure of the lambda repressor-operator complex. J Mol Biol 1992, 227: 177–196. 10.1016/0022-2836(92)90690-L
Thorogood H, Waters TR, Parker AW, Wharton CW, Connolly BA: Resonance Raman spectroscopy of 4-thiothymidine and oligodeoxynucleotides containing this base both free in solution and bound to the restriction endonuclease EcoRV. Biochemistry 1996, 35: 8723–8733. 10.1021/bi960230j
Wenz C, Jeltsch A, Pingoud A: Probing the indirect readout of the restriction enzyme EcoRV. mutational analysis of contacts to the DNA backbone. J Biol Chem 1996, 271: 5565–5573. 10.1074/jbc.271.10.5565
Hartmann B, Sullivan MR, Harris LF: Operator recognition by the phage 434 cI repressor: MD simulations of free and bound 50-bp DNA reveal important differences between the OR1 and OR2 sites. Biopolymers 2003, 68: 250–264. 10.1002/bip.10243
Flader W, Wellenzohn B, Winger RH, Hallbrucker A, Mayer E, Liedl KR: Stepwise induced fit in the pico- to nanosecond time scale governs the complexation of the even-skipped transcriptional repressor homeodomain to DNA. Biopolymers 2003, 68: 139–149. 10.1002/bip.10242
Zhang Y, Zhiqun X, Hegde RS, Shakked Z, Crothers DM: Predicting indirect readout effects in protein-DNA interactions. Proc Natl Acad Sci 2004, 101: 8337–8341. 10.1073/pnas.0402319101
Shakked Z, Rabinovich D: Sequence-dependent conformation of an A-DNA double helix: the crystal structure of the octamer d(G-G-T-A-T-A-C-C). J Mol Biol 1983.
Westhof E, Dumas P, Moras D: Crystallographic refinement of yeast aspartic acid transfer RNA. J Mol Biol 1985, 184: 119–145. 10.1016/0022-2836(85)90048-8
Varnai P, Djuranovic D, Lavery R, Hartmann B: α - γ transitions in the B-DNA backbone. Nucleic Acids Res 2002, 30: 5398–5406. 10.1093/nar/gkf680
Dixit SB, Beveridge DL, Case DA, Cheatham TE III, Giudice E, Lankas F, Lavery R, Maddocks JH, Osman R, Sklenar H, Thayer KM, Varnai P: Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps. Biophys J 2005, 89: 3721–3740. 10.1529/biophysj.105.067397
Madhumalar A, Bansal M: Sequence preference for BI/BII Conformations in DNA: MD and crystal structure data analysis. J Biomol Struct Dyn 2005, 23: 13–27.
Gardiner EJ, Hunter CA, Packer MJ, Palmer DS, Willett P: Sequence-dependent DNA structure: a database of octamer structural parameters. J Mol Biol 2003, 332: 1025–1035. 10.1016/j.jmb.2003.08.006
Beveridge DL, Barreiro G, Byun KS, Case DA, Cheatham TEI, Dixit SB, Giudice E, Lankas F, Lavery R, Maddocks JH, Osman R, Seibert E, Sklenar H, Stoll G, Thayer KM, Varnai P, Young MA: Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps. J Biomol Struct Dyn 2004, 10: 213–226.
Marini JC, Levene SD, Crothers DM, Englund PT: Bent helical structure in kinetoplast DNA. Proc Natl Acad Sci 1982, 79: 7664–7668. 10.1073/pnas.79.24.7664
Marini JC, Effron PN, Goodman TC, Singleton CK, Wells RD, Wartell MR, Englund PT: Physical characterization of a kinetoplast DNA fragment with unusual properties. J Biol Chem 1984, 259: 8974–8979.
Trifonov EN, Sussman JL: The pitch of chromatin DNA is reflected in its nucleotide sequence. Proc Natl Acad Sci 1980, 77: 3816–3820. 10.1073/pnas.77.7.3816
Levene SD, Crothers DM: A computer graphics study of sequence directed bending in DNA. J Biomol Struct Dyn 1983, 1: 429–435.
Brukner I, Jurukovski V, Konstantinovic M, Savic A: Curved DNA without AA/TT dinucleotide step. Nucleic Acids Res 1991, 19: 3549–3551. 10.1093/nar/19.13.3549
Bolshoy A, McNamara P, Harrington RE, Trifonov EN: Curved DNA without A-A: experimental estimation of all 16 DNA wedge angles. Proc Natl Acad Sci 1991, 88: 2312–2316. 10.1073/pnas.88.6.2312
Bansal M: Structural variations observed in DNA crystalstructures and their implications for protein-DNA interaction. In Proceedings of the Ninth Conversation: Biological Structure and Dynamics. New York: Adenine press; 1996:121–134.
Liu Y, Beveridge DL: A refined prediction method for gel retardation of DNA oligonucleotides from dinucleotide step parameters: reconciliation of DNA bending models with crystal structure data. J Biomol Struct Dyn 2001, 18: 505–526.
Bansal M, Bhattacharyya D, Mohanty D: Molecular Conformation and Biological Interactions. Bangalore: Indian Academy of Sciences 1991 chap. DNA bending: a natural consequence of base sequence dependent variability; 347–362.
Srinivasan AR, Torres R, Clark W, Olson WK: Base sequence effects in double helical DNA. I. Potential energy estimates of local base morphology. J Biomol Struct Dyn 1987, 5: 459–496.
Tung C, Garcia AE, Marr TG, Mishra SK, Burks C: Measures of DNA curvature. DNA Bending and Curvature 1988, 3: 265–277.
Dlakic M, Harrington RE: Unconventional helical phasing of repetitive DNA motifs reveals their relative bending contributions. Nucleic Acids Res 1998, 26: 4274–4279. 10.1093/nar/26.18.4274
Kanhere A, Bansal M: An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data. Nucleic Acids Res 2003, 31: 2647–2658. 10.1093/nar/gkg362
Jernigan RL, Sarai A, Shapiro B, Nussinov R: Relationship between curved DNA conformation and slow gel migration. J Biomol Struct Dyn 1987, 4: 561–567.
Olson WK, Marky NL, Jernigan RL, Zhurkin VB: Influence of flutuations on DNA curvature. A comparison of flexible and static wedge models of intrinsically bent DNA. J Mol Biol 1993, 232: 530–554. 10.1006/jmbi.1993.1409
Goodsell DS, Dickerson RE: Bending and curvature calculations in B-DNA. Nucleic Acids Res 1994, 22: 5497–5503. 10.1093/nar/22.24.5497
Harrison SC, Aggarwal AK: DNA recognition by proteins with the helix-turn-helix motif. Annu Rev Biochem 1990, 59: 933–969. 10.1146/annurev.bi.59.070190.004441
Pabo CO, Sauer RT: Protein-DNA recognition. Annu Rev Biochem 1984, 53: 293–321. 10.1146/annurev.bi.53.070184.001453
Pabo CO, Sauer RT: Transcription factors: structural families and principles of DNA recognition. Annu Rev Biochem 1992, 61: 1053–1095. 10.1146/annurev.bi.61.070192.005201
Wintjens R, Rooman M: Structural classification of HTH DNA-binding domains and protein-DNA interaction modes. J Mol Biol 1996, 262: 294–313. 10.1006/jmbi.1996.0514
Berman HM, Olson WK, Beveridge DL, Westbrook J, Gelbin A, Demeny T, Hsieh SH, Srinivasan AR, Schneider B: The Nucleic Acid Database: a comprehensive relational database of three-dimensional structures of nucleic acids. Biophys J 1992, 63: 751–759. 10.1016/S0006-3495(92)81649-1
Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol 1995, 247: 536–540.
Nagaich AK, Bhattacharyya D, Brahmachari SK, Bansal M: CA/TG sequence at the 5' end of oligo(A)-tracts strongly modulates DNA curvature. J Biol Chem 1994, 269: 7824–7833.
Arnott S, Hukins DW: Optimised parameters for A-DNA and B-DNA. Biochem Biophys Res Commun 1972, 47: 1504–1509. 10.1016/0006-291X(72)90243-4
Grzeskowiak K, Yanagi K, Prive GG, Dickerson RE: The structure of B-helical C-G-A-T-C-G-A-T-C-G and comparison with C-C-A-A-C-G-T-T-G-G. The effect of base pair reversals. J Biol Chem 1991, 266: 8861–8883.
Yanagi K, Prive GG, Dickerson R: Analysis of local helix geometry in three B-DNA decamers and eight dodecamers. J Mol Biol 1991, 217: 201–214. 10.1016/0022-2836(91)90620-L
Prive GG, Yanagi K, Dickerson RE: Structure of the B-DNA decamer C-C-A-A-C-G-T-T-G-G and comparison with isomorphous decamers C-C-A-A-G-A-T-T-G-G and C-C-A-G-G-C-C-T-G-G. J Mol Biol 1991, 217: 177–199. 10.1016/0022-2836(91)90619-H
Kielkopf CL, Ding S, Kuhn P, Rees DC: Conformational flexibility of B-DNA at 0.74 Å resolution: d(CCAGTACTGG) 2 . J Mol Biol 2000, 296: 787–801. 10.1006/jmbi.1999.3478
Chiu TK, Dickerson R: 1 Å crystal structure of B-DNA reveal sequence specific binding and groove specific bending of DNA by magnesium and calcium. J Mol Biol 2000, 301: 915–945. 10.1006/jmbi.2000.4012
Davey CA, Sargent DF, Luger K, Maeder AW, Richmond TJ: Solvent Mediated Interactions in the Structure of the Nucleosome Core Particle at 1.9 Å Resolution. J Mol Biol 2002, 319: 1097–1113. 10.1016/S0022-2836(02)00386-8
Ng HL, Kopka ML, Dickerson RE: The structure of a stable intermediate in the A ↔ B DNA helix transition. Proc Natl Acad Sci 2000, 97: 2035–2039. 10.1073/pnas.040571197
Mayer-Jung C, Moras D, Timsit Y: Hydration and recognition of methylated CpG steps in DNA. EMBO J 1998, 17: 2709–2718. 10.1093/emboj/17.9.2709
Luger K, Mader AW, Richmond RK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particle at 2.8 Å resolution. Nature 1997, 389: 251–260. 10.1038/38444
Lejeune D, Delsaux N, Charloteaux B, Thomas A, Brasseur R: Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure. Proteins: Str Func Bioinf 2005, 61: 258–271. 10.1002/prot.20607
Watkins S, van Pouderoyen G, Sixma TK: Structural analysis of the bipartite DNA-binding domain of Tc3 transposase bound to transposon DNA. Nucleic Acids Res 2004, 32: 4306–4312. 10.1093/nar/gkh770
Nekludova L, Pabo CO: Distinctive DNA conformation with enlarged major groove is found in Zn-finger-DNA and other protein-DNA complexes. Proc Natl Acad Sci 1994, 91: 6948–6952. 10.1073/pnas.91.15.6948
Bartfeld D, Shimon L, Couture GC, Rabinovich D, Frolow F, Levanon D, Groner Y, Shakked Z: DNA recognition by the RUNX1 transcription factor is mediated by an allosteric transition in the RUNT domain and by DNA bending. Structure 2002, 10: 1395–1407. 10.1016/S0969-2126(02)00853-5
Mack DR, Chiu TK, Dickerson RE: Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts. J Mol Biol 2001, 312: 1037–1049. 10.1006/jmbi.2001.4994
Littlefield O, Korkhin Y, Sigler PB: The structural basis for the oriented assembly of a TBP/TFB/promoter complex. Proc Natl Acad Sci 1999, 96: 13668–13673. 10.1073/pnas.96.24.13668
Patikoglou G, Kim JL, Sun L, Yang SH, Kodadek T, Burley SK: TATA element recognition by the TATA box-binding protein has been conserved throughout evolution. Genes Dev 1999, 13: 3217–3230. 10.1101/gad.13.24.3217
Savitha G, Viswamitra MA: An A-DNA structure with two independent duplexes in the asymmetric unit. Acta Crystallogr D Biol Crystallogr 1999, 55: 1136–1143. 10.1107/S0907444999003182
Tsai FTF, Sigler PB: Structural basis of preinitiation complex assembly on human pol II promoters. EMBO J 2000, 19: 25–36. 10.1093/emboj/19.1.25
Tolstorukov MY, Jernigan RL, Zhurkin VB: Protein-DNA recognition in the minor groove is facilitated by sugar switching. J Mol Biol 2004, 337: 65–76. 10.1016/j.jmb.2004.01.011
Lankas F, Lavery R, Maddocks JH: Kinking occurs during molecular dynamics simulations of small DNA minicircles. Structure 2006, 14: 1527–1534. 10.1016/j.str.2006.08.004
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: 235–242. 10.1093/nar/28.1.235
McLaughlin WA, Berman H: Statistical models for discerning protein structures containing the DNA-binding helix-turn-helix motif. J Mol Biol 2003, 330: 43–55. 10.1016/S0022-2836(03)00532-1
Bhattacharyya D, Bansal M: A self-consistent formulation for the analysis and generation of non-uniform DNA structures. J Biomol Struct Dyn 1989, 6: 635–653.
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: 281–287.
Mukherjee S, Bansal M, Bhattacharyya D: Conformational specificity of non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis. J Comput Aided Mol Des 2006, 20: 629–645. 10.1007/s10822-006-9083-x
We are grateful to Satish Gupta for providing the RNA dataset for this analysis. We would like to thank Senthil Kumar D. K. for helping with the preparation of additional file 1. This work was supported by Department of Biotechnology, India.
AM and MB conceived the study. AM developed the various tools for analysis and carried out the analysis of the free and complex datasets. DK carried out the analysis of the HTH-bound dataset. AM, MB and DK drafted the manuscript.
Electronic supplementary material
Additional file 1: PDB id's and corresponding references for the X-ray crystal structures used in this study. This file contains the PDB id's of the crystal structures and the citation for each of these structure. (XLS 78 KB)
Additional file 2: Variation of different local and global parameters for the four datasets. This file contains three figures and a table. The figures show the correlated variation of different parameters, namely; twist versus roll, Z p versus backbone torsion angle χ and Z p versus backbone torsion angle δ, for the four datasets used in this study. The table shows the values of the global helical basepair orientation parameters and the minor and major groove widths for the four datasets. (PDF 2 MB)
Additional file 3: Parameters for the distorted dinucleotide steps. This file gives the dinucleotide step parameters and Z p values for the distorted steps in the DNA structures bound to different protein families. (TXT 15 KB)
Additional file 4: Cartoon diagrams of few DNA-protein crystal structure complexes. This file contains the cartoon diagrams of four DNA-protein X-ray crystal structure complexes with different amounts of curvature. (PDF 247 KB)
Authors’ original submitted files for images
About this article
Cite this article
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 9, 24 (2009). https://doi.org/10.1186/1472-6807-9-24
- Backbone Torsion Angle
- Dinucleotide Step
- Helical Twist
- Sugar Pucker
- Glycosidic Torsion Angle