This Article OPEN ACCESS ARTICLE Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Greenbaum, J. A. Articles by Tullius, T. D. Search for Related Content PubMed PubMed Citation Articles by Greenbaum, J. A. Articles by Tullius, T. D. Social Bookmarking                   What's this?

Resource

Construction of a genome-scale structural map at single-nucleotide resolution

¹ Program in Bioinformatics, Boston University, Boston, Massachusetts 02215, USA; ² Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA

	Abstract

Top Abstract Results Discussion Methods Acknowledgments References

 
Few methods are available for mapping the local structure of DNA throughout a genome. The hydroxyl radical cleavage pattern is a measure of the local variation in solvent-accessible surface area of duplex DNA, and thus provides information on the local shape and structure of DNA. We report the construction of a relational database, ORChID (OH Radical Cleavage Intensity Database), that contains extensive hydroxyl radical cleavage data produced from two DNA libraries. We have used the ORChID database to develop a set of algorithms that are capable of predicting the hydroxyl radical cleavage pattern of a DNA sequence of essentially any length, to high accuracy. We have used the prediction algorithm to produce a structural map of the 30 Mb of the ENCODE regions of the human genome.

We report here the construction of a library of hydroxyl radical cleavage patterns of DNA, as a means of compiling structural information for a wide variety of DNA sequences. Although the hydroxyl radical cleavage pattern (Price and Tullius 1992) does not yield a high-resolution three-dimensional structure of a DNA molecule, it is a reflection of an important structural parameter, the solvent-accessible surface area of the DNA backbone (Balasubramanian et al. 1998). The cleavage pattern thus provides an image of the shape of the DNA backbone and how it varies with respect to nucleotide sequence. We describe the use of a fluorescence-based sequencer to obtain cleavage patterns, and introduce methods for normalization and quantitation of cleavage data. We present the design considerations of a relational database, ORChID (OH Radical Cleavage Intensity Database), to hold hydroxyl radical cleavage data, as well as its population and user interface.

We take advantage of the ORChID database to investigate how the nucleotide sequence of a DNA molecule affects its pattern of hydroxyl radical cleavage. We show here that, as expected, similar sequences yield similar cleavage patterns. However, we also find instances of long segments of DNA with low nucleotide sequence identity that produce nearly identical cleavage patterns. That is, it is possible for different DNA sequences to yield similar cleavage profiles. As the cleavage profile is a reflection of the underlying DNA structure, this indicates that considerably different DNA sequences can share a common structure. This leads to the intriguing possibility that DNA structure might be evolutionarily conserved, irrespective of the sequence of nucleotides.

Finally, we use the ORChID database to construct an algorithm that allows the prediction of the hydroxyl radical cleavage pattern of any DNA sequence to high accuracy. The speed of this algorithm makes it feasible to produce a structural map of a large genome. Here we report the use of this algorithm to construct a structural map of the 30 Mb of the ENCODE regions of the human genome.

	Results

Top Abstract Results Discussion Methods Acknowledgments References

 
We began the construction of a database of hydroxyl radical cleavage patterns by obtaining two different libraries of single-stranded DNA molecules (Supplemental Table S1). One library, R40, was a collection of 158-nt-long DNA molecules synthesized with a segment of 40 random nucleotides (nt) in the center. The other library, pentamer, consisted of 14 members (123 or 118 nt in length), each of which contained a subset of all the 1024 possible pentanucleotides (Supplemental Fig. S1). In both libraries, the test sequence at the center was flanked by a common palindromic sequence on either side (Fig. 1A), to aid in data normalization.

View larger version (34K): [in this window] [in a new window]   Figure 1. (A) Design of the DNA molecules used to construct the R40 library. The test sequence (an insert of 40 random nucleotides) is located near the center of the DNA strand, flanked by common sequences on both sides. The pentamer library was constructed in a similar manner, with individual pentamer sequences (see Supplemental Fig. S1) serving as the test sequence. (B) Typical electropherogram of a sample from the R40 library. The R40 test sequence is boxed, and the common flanking sequences are shaded. Each peak in the pattern represents cleavage by the hydroxyl radical at one nucleotide of the DNA molecule. The area of a peak is proportional to the extent of cleavage at that nucleotide (Shadle et al. 1997). Note that the cleavage patterns of the common palindromic flanking sequences at the 5' (left) and 3' (right) ends are similar for a particular member of the library. This also holds true for different library members (data not shown).

Database layout and design
Careful consideration was taken to store the data in a meaningful, intuitive manner, while minimizing redundancy. To achieve this, one central table was established to which most other tables have a relationship (Supplemental Fig. S2). Within the central Sample Table, information, including the sequence, individual peak areas of the hydroxyl radical cleavage pattern, gel id, and gel lane, is stored. Both raw and normalized cleavage data are incorporated into the database.

To enhance the usability of the ORChID database, some of the more commonly queried data were combined into views, exemplified by trimers (Supplemental Table S2) and trimer summary (Supplemental Table S3). Corresponding views exist for N-mers ranging from monomers through septamers. An important use of these views involves the hydroxyl radical cleavage prediction algorithm that is discussed below. The algorithm is based on a sliding N-mer window model and thus requires the mean area for each peak of each N-mer for its calculations. Collection of these mean peak areas in the Summary views greatly reduces the complexity of the SQL statement that is used to access the relevant data.

Database usage
The ORChID database can be accessed directly through PostgreSQL’s psql interface or through a graphical Web interface at http://dna.bu.edu/orchid. From the Web interface, the user has the option to query the database by sequence ID, by nucleotide sequence, or by entering any properly formed SELECT statement. A separate page allows the user to make plots of cleavage intensity data for any sequence or pair of sequences that exist in the database. This functionality permits one to explore the database and investigate the properties of hydroxyl radical cleavage data.

The most widely useful feature of the Web interface is the ability to calculate a predicted hydroxyl radical cleavage pattern for any given sequence. From the Prediction Page, the user inputs a DNA sequence of nearly any length, and receives tabular and graphical output of the predicted hydroxyl radical cleavage pattern in a few seconds. Several options for prediction are provided, including the use of different prediction algorithms and output settings. Details of the prediction algorithm are discussed below.

Reproducibility of hydroxyl radical cleavage data
To quantitatively assess the quality of the hydroxyl radical cleavage data that we obtain by our experimental scheme, we took advantage of the pair of identical flanking sequences that are present in each member of the two DNA libraries (Fig. 1). We extracted from the ORChID database 112 examples of independently measured cleavage patterns for the common dodecanucleotide sequence d(CGCGAATTCGCG) and evaluated the reproducibility of the cleavage pattern. Figure 2 depicts the mean cleavage intensity for each position in the sequence, with the standard deviations indicated as error bars. The per-position standard deviations range from 0.14 to 0.35, which is between 5% and 13% of the complete range of values in the cleavage pattern, 2.67. The majority of the positions have a standard deviation of 0.26 or less, which translates to <10% of the range. The range of standard deviations we observe is consistent with previous work on quantitative analysis of hydroxyl radical cleavage patterns (Shadle et al. 1997).

View larger version (15K): [in this window] [in a new window]   Figure 2. Reproducibility of the cleavage pattern of the common flanking sequence (see Fig. 1). The mean cleavage intensities, taken from 112 instances of the common flanking sequence cleavage pattern, are plotted as bars, with standard deviations shown as error bars.

To investigate this question, we divided the sequences in the ORChID database into overlapping N-mers ranging from 8 to 34 nucleotides (nt) in length, and calculated Pearson correlations for all pairwise cleavage pattern comparisons. Similarly, for each pair of N-mers, we calculated the degree of nucleotide sequence identity. We then determined the relationship between sequence identity and cleavage pattern similarity (Supplemental Table S4). Given the notion that similar DNA sequences share a common structure, one would expect that sequences with a high degree of identity would also exhibit similar cleavage patterns. However, we found that overall, the Pearson correlation of sequence identity and cleavage similarity is rather low, 0.36 for the N-mer lengths we studied.

By dividing these data into subsets of similar sequence identity and then binning into discrete levels of cleavage similarity (Supplemental Tables S5–S7), we obtain a clearer sense of the relationship between cleavage pattern similarity and sequence identity. Despite the low Pearson correlation between these two parameters (Supplemental Table S4), the heatmaps shown in Figure 3 clearly illustrate that they are tightly linked.

View larger version (40K): [in this window] [in a new window]   Figure 3. Cleavage/sequence correlation at various levels of sequence identity. Heat maps were created with Matrix2png (Pavlidis and Noble 2003), using the data from Supplemental Tables S5, S6, and S7 as input. (A) 10-mers; (B) 20-mers; (C) 30-mers. The intensity of each rectangle indicates the number of pairs of sequences in that cell, ranging from black (lowest) to white (highest).

We next examined particular examples of cleavage patterns of some of the outlier sequences. Figure 4 depicts two 10-mers having completely different sequences (i.e., 0% identity), but with a Pearson coefficient of 0.94 between their cleavage patterns. Supplemental Figure S3 shows another example, two 20-mers with 10% sequence identity, yet a cleavage pattern correlation coefficient of 0.81. (The sequences of these four DNA molecules are listed in Supplemental Table S8.) These two examples illustrate the idea that two or more divergent sequences can share similar cleavage patterns over a relatively long stretch of DNA.

View larger version (9K): [in this window] [in a new window]   Figure 4. High similarity in cleavage patterns for two sequences with low sequence identity. Plotted are the hydroxyl radical cleavage patterns of two 10-mer sequences that share no common nucleotides (sequence identity = 0%). Note the significant correlation (R = 0.94) of the two patterns. (See Supplemental Table S8 for sequences.)

The prediction algorithms all involve treating a DNA sequence as being made up of overlapping N-mers. As an example, we discuss the Sliding Trimer Window algorithm. We also have implemented several related higher-order prediction algorithms. An overview of the Sliding Trimer Window algorithm is presented in Figure 5. This algorithm works by dividing the target sequence into overlapping trimers, and then retrieving the corresponding cleavage data from the ORChID database. We obtain the predicted hydroxyl radical cleavage intensity for each nucleotide in the sequence by taking the average of the three cleavage intensities at each position that are contributed by the three overlapping trinucleotides that are associated with that nucleotide. Figure 6 depicts the predicted and observed patterns for one sample from the ORChID database, for which the Pearson correlation between experimental and predicted pattern is 0.91. The correlation between the predicted and experimental datasets is striking, particularly given the simplicity of the model.

View larger version (24K): [in this window] [in a new window]   Figure 5. Sliding trimer window algorithm. The sequence to be predicted is shown below the bar graph. It is divided into overlapping trinucleotides. The hydroxyl radical cleavage intensity data for each trimer are retrieved from the ORChID database and are listed below each nucleotide. The values in each column are averaged to produce a predicted hydroxyl radical cleavage intensity, which is represented as a bar at the top. Note that the two terminal nucleotides at each end rely on data from only one or two trinucleotides, rather than three.

View larger version (20K): [in this window] [in a new window]   Figure 6. Predicted hydroxyl radical cleavage pattern. The hydroxyl radical cleavage pattern of sample ID 25201 from the ORChID database was predicted using the Sliding Trimer Window algorithm (see Fig. 5). The predicted pattern (broken line) is compared to the experimental pattern (solid line). The experimental and predicted patterns have a correlation coefficient of 0.91.

View this table: [in this window] [in a new window]   Table 1. Correlation between predicted and experimental cleavage patterns

	Discussion

Top Abstract Results Discussion Methods Acknowledgments References

 
While methods like X-ray crystallography and NMR produce detailed structural information for DNA, they are severely limited both by the time required and the length of the DNA molecule that can be studied. Although many crystal and NMR structures of DNA are currently available (Berman et al. 1992, 2000), most of these are no longer than 12 base pairs (bp) in length. Furthermore, it has been suggested that crystallization conditions enforce a static structure that may not be truly representative of DNA in solution (Dickerson et al. 1987; DiGabriele et al. 1989; Dlakic et al. 1996; Ganunis et al. 1996; Berman 1997; Ghosh and Bansal 2001). Despite recent advances in NMR (Zhou et al. 1999), the length of a DNA molecule that can be studied by this method is still rather limited.

Here we have described an alternative method for the acquisition of DNA structural information based on the collection of a library of hydroxyl radical cleavage patterns of DNA. Our approach can be made into a high-throughput method for determining structural features of DNA molecules. By organizing this information into a database, we have been able to study the effect of the sequence of a DNA molecule on its structure and make several key observations. We discovered that it is possible for highly divergent DNA sequences to produce closely related hydroxyl radical cleavage patterns; this is an indication that these stretches of DNA have a similar backbone shape. For example, there are 36 pairs of decamer sequences in the ORChID database that have only one nucleotide in common, yet have a Pearson coefficient of >0.9 when their cleavage patterns are compared (Supplemental Table S5). Lowering the Pearson coefficient threshold to 0.8 increases the number of structurally similar sequence pairs nearly 10-fold, to 355. These results indicate that the structural similarity of highly divergent sequences is common.

An intriguing implication of this finding is the role that DNA structural similarity may play in the binding of proteins to DNA. Most known transcription factor-binding sites (Matys et al. 2006; Vlieghe et al. 2006) are between 8 and 12 bp in length, and the vast majority of these sites are highly variable. Our work has shown the potential for divergent DNA sequences to adopt a similar backbone conformation over a stretch of 20 nt or so. It therefore is reasonable to propose that some DNA-binding proteins may be backbone conformation-specific (i.e., using indirect readout), rather than DNA sequence-specific. The development of algorithms to identify sequences with degenerate cleavage patterns may therefore prove useful in understanding how transcription factors locate their binding sites in genomic DNA. As well, since it has previously been shown that certain DNA sequences can act by virtue of their structure (Bracco et al. 1989; Kim et al. 1995; Angermayr et al. 2002), methods to identify such structures would benefit from incorporating hydroxyl radical cleavage data.

By studying the general features of hydroxyl radical cleavage data and organizing the data in meaningful ways, we have successfully developed algorithms for the prediction of hydroxyl radical cleavage patterns of DNA. Given the speed of our prediction algorithms, along with the addition of more sequences to the ORChID database and the concurrent development of higher-order predictive models, the experimental determination of hydroxyl radical cleavage patterns of DNA will become unnecessary. We look forward to the further use of ORChID data to help to understand functional regions of genomes in terms of their local DNA structural features.

	Methods

Top Abstract Results Discussion Methods Acknowledgments References

 
Hydroxyl radical cleavage of DNA
Twenty microliters of fluorescently labeled DNA (see Supplemental Material) and 50 µL of buffer (20 mM Tris, 20 mM NaCl at pH 8) were pipetted into a 1.5-mL Eppendorf tube. Next, 10-µL drops of [Fe(EDTA)]^2– (50 µM) and ascorbate (1 mM) were pipetted onto the wall of the tube, but not mixed. To initiate the reaction, 10 µL of H₂O₂ (0.03%) was combined with the other two reagents and mixed into the DNA sample in buffer by vigorous pipetting. The reaction was quenched after 2 min by the addition of 400 µL of ethanol (100%) and vortexing. The DNA was ethanol-precipitated. The dried DNA pellet was dissolved in 6 µL of formamide loading dye, and the sample was electrophoresed on a denaturing polyacrylamide gel using a Visible Genetics Long-Read Tower automated sequencer.

Data quantitation
The flanking common palindromic sequences (Fig. 1) were located by visual inspection of the data set, and their positions were recorded. Next, the text file containing the fluorescence intensity data was parsed and converted into a format readable by Origin (OriginLab Software). The fluorescence trace was plotted in Origin, and a value for the baseline was determined by noting the fluorescence intensity just before the first DNA fragment was observed. The raw data were then read into PeakFit (Systat Software), and the baseline was subtracted. The peaks in the data set were simultaneously fit using the EMG + GMG peak function. The area of a peak was taken to represent the amount of cleavage at that nucleotide. We then used the peak areas of the flanking common palindromic sequences to normalize the data for the test sequence. Further details are provided in the Supplemental Material.

Sliding N-mer window algorithms
The method illustrated in Figure 5 is applicable to all of our sliding N-mer algorithms. As an example, we describe here the Sliding Trimer Window algorithm. First, the sequence for which the cleavage pattern is to be predicted is divided into overlapping trinucleotide sequences. The corresponding cleavage intensities for each trinucleotide are retrieved from the ORChID database. The average cleavage intensity is calculated for each position in the sequence of interest, using Equation 1:

(1)

Here, P_i is the predicted cleavage intensity at position i, and T_i is the trinucleotide beginning at position i. The superscript on T indicates the nucleotide in the trimer (e.g., T¹ is the first nucleotide of the trimer). The ends of the pattern are calculated similarly, except that cleavage data are retrieved from only one or two trimers, rather than three (see Fig. 5).

Correlation of predicted patterns with experimental patterns
Hydroxyl radical cleavage patterns were predicted for each of 78 sequences in the ORChID database: 56 pentamer library sequences (14 different DNA duplexes, both strands, two independent experimental examples of each), and 22 sequences from the R40 library. Before a particular cleavage pattern was predicted, its corresponding cleavage data were removed from the database; the data were replaced in the database upon completion of the prediction. The ends of the predicted and observed patterns were trimmed by several nucleotides in order to include only the most accurate predictions. Pearson correlation coefficients were calculated for the trimmed patterns versus the experimental patterns (see Table 1).

	Acknowledgments

Top Abstract Results Discussion Methods Acknowledgments References

 
This work was funded by an ENCODE Technology Development grant from the National Human Genome Research Institute of the National Institutes of Health (R01 HG003541). J.A.G. was supported by an IGERT training grant from the National Science Foundation (DGE-9870710). We thank Steve Parker and Eric Bishop for helpful discussions.

	Footnotes

 
³ Corresponding author.

E-mail tullius{at}bu.edu; fax (617) 353-6466.

[Supplemental material is available online at www.genome.org.]

Article is online at http://www.genome.org/cgi/doi/10.1101/gr.6073107

	References

Top Abstract Results Discussion Methods Acknowledgments References

 

Angermayr, M., Oechsner, U., Gregor, K., Schroth, G.P., and Bandlow, W. 2002. Transcription initiation in vivo without classical transactivators: DNA kinks flanking the core promoter of the housekeeping yeast adenylate kinase gene, AKY2, position nucleosomes and constitutively activate transcription. Nucleic Acids Res. 30: 4199–4207.[Abstract/Free Full Text]

Balasubramanian, B., Pogozelski, W.K., and Tullius, T.D. 1998. DNA strand breaking by the hydroxyl radical is governed by the accessible surface areas of the hydrogen atoms of the DNA backbone. Proc. Natl. Acad. Sci. 95: 9738–9743.[Abstract/Free Full Text]

Barbic, A., Zimmer, D.P., and Crothers, D.M. 2003. Structural origins of adenine-tract bending. Proc. Natl. Acad. Sci. 100: 2369–2373.[Abstract/Free Full Text]

Berman, H.M. 1997. Crystal studies of B-DNA: The answers and the questions. Biopolymers 44: 23–44.[CrossRef][Medline]

Berman, H.M., Olson, W.K., Beveridge, D.L., Westbrook, J., Gelbin, A., Demeny, T., Hsieh, S.H., Srinivasan, A.R., and Schneider, B. 1992. The nucleic acid database. A comprehensive relational database of three-dimensional structures of nucleic acids. Biophys. J. 63: 751–759.[Medline]

Berman, H.M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T.N., Weissig, H., Shindyalov, I.N., and Bourne, P.E. 2000. The Protein Data Bank. Nucleic Acids Res. 28: 235–242.[Abstract/Free Full Text]

Bhattacharyya, D. and Bansal, M. 1990. Local variability and base sequence effects in DNA crystal structures. J. Biomol. Struct. Dyn. 8: 539–572.[Medline]

Bracco, L., Kotlarz, D., Kolb, A., Diekmann, S., and Buc, H. 1989. Synthetic curved DNA sequences can act as transcriptional activators in Escherichia coli. EMBO J. 8: 4289–4296.[Medline]

Calladine, C.R. 1982. Mechanics of sequence-dependent stacking of bases in B-DNA. J. Mol. Biol. 161: 343–352.[CrossRef][Medline]

Calladine, C.R. and Drew, H.R. 1986. Principles of sequence-dependent flexure of DNA. J. Mol. Biol. 192: 907–918.[CrossRef][Medline]

Dickerson, R.E. 1983. Base sequence and helix structure variation in B and A DNA. J. Mol. Biol. 166: 419–441.[Medline]

Dickerson, R.E. 1992. DNA structure from A to Z. Methods Enzymol. 211: 67–111.[Medline]

Dickerson, R.E. 1997. Sequence-dependent helix deformability in the recognition of B-DNA. Biopolymers 44: 321.[CrossRef][Medline]

Dickerson, R.E. and Drew, H.R. 1981. Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure. J. Mol. Biol. 149: 761–786.[CrossRef][Medline]

Dickerson, R.E., Goodsell, D.S., Kopka, M.L., and Pjura, P.E. 1987. The effect of crystal packing on oligonucleotide double helix structure. J. Biomol. Struct. Dyn. 5: 557–579.[Medline]

Diekmann, S., von Kitzing, E., McLaughlin, L., Ott, J., and Eckstein, F. 1987. The influence of exocyclic substituents of purine bases on DNA curvature. Proc. Natl. Acad. Sci. 84: 8257–8261.[Abstract/Free Full Text]

DiGabriele, A.D. and Steitz, T.A. 1993. A DNA dodecamer containing an adenine tract crystallizes in a unique lattice and exhibits a new bend. J. Mol. Biol. 231: 1024–1039.[CrossRef][Medline]

DiGabriele, A.D., Sanderson, M.R., and Steitz, T.A. 1989. Crystal lattice packing is important in determining the bend of a DNA dodecamer containing an adenine tract. Proc. Natl. Acad. Sci. 86: 1816–1820.[Abstract/Free Full Text]

Dlakic, M., Park, K., Griffith, J.D., Harvey, S.C., and Harrington, R.E. 1996. The organic crystallizing agent 2-methyl-2,4-pentanediol reduces DNA curvature by means of structural changes in A-tracts. J. Biol. Chem. 271: 17911–17919.[Abstract/Free Full Text]

El Hassan, M.A. and Calladine, C.R. 1997. Conformational characteristics of DNA: Empirical classifications and a hypothesis for the conformational behavior of dinucleotide steps. Philos. Trans. R. Soc. Lond. A 355: 43–100.[CrossRef]

The ENCODE Project Consortium, 2004. The ENCODE (ENCyclopedia Of DNA Elements) Project. Science 306: 636–640.[Abstract/Free Full Text]

Ganunis, R.M., Guo, H., and Tullius, T.D. 1996. Effect of the crystallizing agent 2-methyl-2,4-pentanediol on the structure of adenine tract DNA in solution. Biochemistry 35: 13729–13732.[CrossRef][Medline]

Gardiner, E.J., Hunter, C.A., Packer, M.J., Palmer, D.S., and Willett, P. 2003. Sequence-dependent DNA structure: A database of octamer structural parameters. J. Mol. Biol. 332: 1025–1035.[CrossRef][Medline]

Ghosh, A. and Bansal, M. 2001. Structural features of B-DNA dodecamer crystal structures: Influence of crystal packing versus base sequence. Indian J. Biochem. Biophys. 38: 7–15.[Medline]

Greenbaum, J.A., Parker, S.C.J., and Tullius, T.D. 2007. Detection of DNA structural motifs in functional genomic elements. Genome Res. (this issue) doi: 10.1101/gr.5602807.[Abstract/Free Full Text]

Grzeskowiak, K. 1996. Sequence-dependent structural variation in B-DNA. Chem. Biol. 3: 785–790.[CrossRef][Medline]

Hays, F.A., Teegarden, A., Jones, Z.J., Harms, M., Raup, D., Watson, J., Cavaliere, E., and Ho, P.S. 2005. How sequence defines structure: A crystallographic map of DNA structure and conformation. Proc. Natl. Acad. Sci. 102: 7157–7162.[Abstract/Free Full Text]

Johansson, E., Parkinson, G., and Neidle, S. 2000. A new crystal form for the dodecamer C-G-C-G-A-A-T-T-C-G-C-G: Symmetry effects on sequence-dependent DNA structure. J. Mol. Biol. 300: 551–561.[CrossRef][Medline]

Karolchik, D., Baertsch, R., Diekhans, M., Furey, T.S., Hinrichs, A., Lu, Y.T., Roskin, K.M., Schwartz, M., Sugnet, C.W., Thomas, D.J., et al. 2003. The UCSC Genome Browser Database. Nucleic Acids Res. 31: 51–54.[Abstract/Free Full Text]

Kim, J., Klooster, S., and Shapiro, D.J. 1995. Intrinsically bent DNA in a eukaryotic transcription factor recognition sequence potentiates transcription activation. J. Biol. Chem. 270: 1282–1288.[Abstract/Free Full Text]

Koo, H.S. and Crothers, D.M. 1987. Chemical determinants of DNA bending at adenine-thymine tracts. Biochemistry 26: 3745–3748.[CrossRef][Medline]

Matys, V., Kel-Margoulis, O.V., Fricke, E., Liebich, I., Land, S., Barre-Dirrie, A., Reuter, I., Chekmenev, D., Krull, M., Hornischer, K., et al. 2006. TRANSFAC and its module TRANSCompel: Transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 34: D108–D110.[Abstract/Free Full Text]

Nelson, H.C., Finch, J.T., Luisi, B.F., and Klug, A. 1987. The structure of an oligo(dA) • oligo(dT) tract and its biological implications. Nature 330: 221–226.[CrossRef][Medline]

Ng, H. and Dickerson, R.E. 2001. Mildly eccentric ‘E-DNA.’ Nat. Struct. Biol. 8: 107–108.[Medline]

Olson, W.K., Gorin, A.A., Lu, X.J., Hock, L.M., and Zhurkin, V.B. 1998. DNA sequence-dependent deformability deduced from protein–DNA crystal complexes. Proc. Natl. Acad. Sci. 95: 11163–11168.[Abstract/Free Full Text]

Packer, M.J., Dauncey, M.P., and Hunter, C.A. 2000a. Sequence-dependent DNA structure: Dinucleotide conformational maps. J. Mol. Biol. 295: 71–83.[CrossRef][Medline]

Packer, M.J., Dauncey, M.P., and Hunter, C.A. 2000b. Sequence-dependent DNA structure: Tetranucleotide conformational maps. J. Mol. Biol. 295: 85–103.[CrossRef][Medline]

Pavlidis, P. and Noble, W.S. 2003. Matrix2png: A utility for visualizing matrix data. Bioinformatics 19: 295–296.[Abstract/Free Full Text]

Price, M.A. and Tullius, T.D. 1992. Using hydroxyl radical to probe DNA structure. Methods Enzymol. 212: 194–219.[Medline]

Shadle, S.E., Allen, D.F., Guo, H., Pogozelski, W.K., Bashkin, J.S., and Tullius, T.D. 1997. Quantitative analysis of electrophoresis data: Novel curve fitting methodology and its application to the determination of a protein–DNA binding constant. Nucleic Acids Res. 25: 850–860.[Abstract/Free Full Text]

Vlieghe, D., Sandelin, A., De Bleser, P.J., Vleminckx, K., Wasserman, W.W., van Roy, F., and Lenhard, B. 2006. A new generation of JASPAR, the open-access repository for transcription factor binding site profiles. Nucleic Acids Res. 34: D95–D97.[Abstract/Free Full Text]

Yanagi, K., Prive, G.G., and Dickerson, R.E. 1991. Analysis of local helix geometry in three B-DNA decamers and eight dodecamers. J. Mol. Biol. 217: 201–214.[CrossRef][Medline]

Zhou, H., Vermeulen, A., Jucker, F.M., and Pardi, A. 1999. Incorporating residual dipolar couplings into the NMR solution structure determination of nucleic acids. Biopolymers 52: 168–180.[CrossRef][Medline]

Received October 24, 2006; accepted in revised format January 29, 2007.

CiteULike

Connotea

Del.icio.us

Digg

Reddit

Technorati

This article has been cited by other articles:

				T. Abeel, Y. Saeys, E. Bonnet, P. Rouze, and Y. Van de Peer Generic eukaryotic core promoter prediction using structural features of DNA Genome Res., February 1, 2008; 18(2): 310 - 323. [Abstract] [Full Text] [PDF]

				J. A. Greenbaum, S. C.J. Parker, and T. D. Tullius Detection of DNA structural motifs in functional genomic elements Genome Res., June 1, 2007; 17(6): 940 - 946. [Abstract] [Full Text] [PDF]

This Article OPEN ACCESS ARTICLE Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Greenbaum, J. A. Articles by Tullius, T. D. Search for Related Content PubMed PubMed Citation Articles by Greenbaum, J. A. Articles by Tullius, T. D. Social Bookmarking                   What's this?