Experimental maps of DNA structure at nucleotide resolution distinguish intrinsic from protein-induced DNA deformations.

Recognition of DNA by proteins depends on DNA sequence and structure. Often unanswered is whether the structure of naked DNA persists in a protein-DNA complex, or whether protein binding changes DNA shape. While X-ray structures of protein-DNA complexes are numerous, the structure of naked cognate DNA is seldom available experimentally. We present here an experimental and computational analysis pipeline that uses hydroxyl radical cleavage to map, at single-nucleotide resolution, DNA minor groove width, a recognition feature widely exploited by proteins. For 11 protein-DNA complexes, we compared experimental maps of naked DNA minor groove width with minor groove width measured from X-ray co-crystal structures. Seven sites had similar minor groove widths as naked DNA and when bound to protein. For four sites, part of the DNA in the complex had the same structure as naked DNA, and part changed structure upon protein binding. We compared the experimental map with minor groove patterns of DNA predicted by two computational approaches, DNAshape and ORChID2, and found good but not perfect concordance with both. This experimental approach will be useful in mapping structures of DNA sequences for which high-resolution structural data are unavailable. This approach allows probing of protein family-dependent readout mechanisms.

GBshape: a genome browser database for DNA shape annotations.

Many regulatory mechanisms require a high degree of specificity in protein-DNA binding. Nucleotide sequence does not provide an answer to the question of why a protein binds only to a small subset of the many putative binding sites in the genome that share the same core motif. Whereas higher-order effects, such as chromatin accessibility, cooperativity and cofactors, have been described, DNA shape recently gained attention as another feature that fine-tunes the DNA binding specificities of some transcription factor families. Our Genome Browser for DNA shape annotations (GBshape; freely available at http://rohslab.cmb.usc.edu/GBshape/) provides minor groove width, propeller twist, roll, helix twist and hydroxyl radical cleavage predictions for the entire genomes of 94 organisms. Additional genomes can easily be added using the GBshape framework. GBshape can be used to visualize DNA shape annotations qualitatively in a genome browser track format, and to download quantitative values of DNA shape features as a function of genomic position at nucleotide resolution. As biological applications, we illustrate the periodicity of DNA shape features that are present in nucleosome-occupied sequences from human, fly and worm, and we demonstrate structural similarities between transcription start sites in the genomes of four Drosophila species.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution.

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

The relationship between fine scale DNA structure, GC content, and functional elements in 1% of the human genome.

GC content has been shown to be an important aspect of human genomic function. Extending beyond the scope of GC content alone, there is a class of regions in the genome that have especially high GC content and are enriched for the CG dinucleotide--called CpG islands. CpG islands have been linked to biologically functional genomic elements. DNA structure also contributes to biological function. Recent studies found that some DNA structural properties are correlated with CpG island functionality. Here, we use hydroxyl radical cleavage patterns as a measure of DNA structure, to explore the relationship between GC content and fine-scale DNA structure. We show that there is a positive correlation between GC content and the solvent-accessible structural properties of a DNA sequence, and that the strength of this correlation decreases as genomic resolution increases. We demonstrate that regions of the genome that have highly solvent-accessible DNA structure tend to overlap functional genomic elements. Our results suggest that fine-scale DNA structural properties that are encoded in the genome are important for biological function, and that the highly solvent-accessible nature of high GC content regions and some CpG islands may account for some of their functional properties.

Mapping nucleic acid structure by hydroxyl radical cleavage.

Hydroxyl radical footprinting is a widely used method for following the folding of RNA molecules in solution. This method has the unique ability to provide experimental information on the solvent accessibility of each nucleotide in an RNA molecule, so that the folding of all domains of the RNA species can be followed simultaneously at single-nucleotide resolution. In recent work, hydroxyl radical footprinting has been used, often in combination with other global measures of structure, to work out detailed folding pathways and three-dimensional structures for increasingly large and complicated RNA molecules. These include synthetic ribozymes, and group I and group II ribozymes, from yeast, the Azoarcus cyanobacterium and Tetrahymena thermophila. Advances have been made in methods for analysis of hydroxyl radical data, so that the large datasets that result from kinetic folding experiments can be analyzed in a semi-automated and quantitative manner.

Interchange of DNA-binding modes in the deformed and ultrabithorax homeodomains: a structural role for the N-terminal arm.

The deformed (Dfd) and ultrabithorax (Ubx) homeoproteins regulate developmental gene expression in Drosophila melanogaster by binding to specific DNA sequences within its genome. DNA binding is largely accomplished via a highly conserved helix-turn-helix DNA-binding domain that is known as a homeodomain (HD). Despite nearly identical DNA recognition helices and similar target DNA sequence preferences, the in vivo functions of the two proteins are quite different. We have previously revealed differences between the two HDs in their interactions with DNA. In an effort to define the individual roles of the HD N-terminal arm and recognition helix in sequence-specific binding, we have characterized the structural details of two Dfd/Ubx chimeric HDs in complex with both the Dfd and Ubx-optimal-binding site sequences. We utilized hydroxyl radical cleavage of DNA to assess the positioning of the proteins on the binding sites. The effects of missing nucleosides and purine methylation on HD binding were also analyzed. Our results show that both the Dfd and Ubx HDs have similar DNA-binding modes when in complex with the Ubx-optimal site. There are subtle but reproducible differences in these modes that are completely interchanged when the Dfd N-terminal arm is replaced with the corresponding region of the Ubx HD. In contrast, we showed previously that the Dfd-optimal site sequence elicits a very different binding mode for the Ubx HD, while the Dfd HD maintains a mode similar to that elicited by the Ubx-optimal site. Our current methylation interference studies suggest that this alternate binding mode involves interaction of the Ubx N-terminal arm with the minor groove on the opposite face of DNA relative to the major groove that is occupied by the recognition helix. As judged by hydroxyl radical footprinting and the missing nucleoside experiment, it appears that interaction of the Ubx recognition helix with the DNA major groove is reduced. Replacing the Dfd N-terminal arm with that of Ubx does not elicit a complete interchange of the DNA-binding mode. Although the position of the chimera relative to DNA, as judged by hydroxyl radical footprinting, is similar to that of the Dfd HD, the missing nucleoside and methylation interference patterns resemble those of the Ubx HD. Repositioning of amino acid side-chains without wholesale structural alteration in the polypeptide appears to occur as a function of N-terminal arm identity and DNA-binding site sequence. Complete interchange of binding modes was achieved only by replacement of the Dfd N-terminal arm and the recognition helix plus 13 carboxyl-terminal residues with the corresponding residues of Ubx. The position of the N-terminal arm in the DNA minor groove appears to differ in a manner that depends on the two base-pair differences between the Dfd and Ubx-optimal-binding sites. Thus, N-terminal arm position dictates the binding mode and the interaction of the recognition helix with nucleosides in the major groove.

A computational method to search for DNA structural motifs in functional genomic elements.

The rapidly increasing availability of DNA sequence data from modern high-throughput experimental techniques has created the need for computational algorithms to aid in motif discovery in genomic DNA. Such algorithms are typically used to find a statistical representation of the nucleotide sequence of the target site of a DNA-binding protein within a collection of DNA sequences that are thought to contain segments to which the protein is bound. A major assumption of these algorithms is that the protein recognizes the primary order of nucleotides in the sequence. However, proteins can also recognize the three-dimensional shape and structure of DNA. To account for this, we developed a computational method to predict the local structural profiles of any set of DNA sequences and then to search within these profiles for common DNA structural motifs. Here we describe the details of this method and use it to find a DNA structural motif in the Saccharomyces cerevisiae yeast genome that is associated with binding of the transcription factor RLM1, a component of the protein kinase C-mediated MAP kinase pathway.

DNA shape, genetic codes, and evolution.

Although the three-letter genetic code that maps nucleotide sequence to protein sequence is well known, there must exist other codes that are embedded in the human genome. Recent work points to sequence-dependent variation in DNA shape as one mechanism by which regulatory and other information could be encoded in DNA. Recent advances include the discovery of shape-dependent recognition of DNA that depends on minor groove width and electrostatics, the existence of overlapping codes in protein-coding regions of the genome, and evolutionary selection for compensatory changes in nucleotide composition that facilitate nucleosome occupancy. It is becoming clear that DNA shape is important to biological function, and therefore will be subject to evolutionary constraint.

A map of minor groove shape and electrostatic potential from hydroxyl radical cleavage patterns of DNA.

DNA shape variation and the associated variation in minor groove electrostatic potential are widely exploited by proteins for DNA recognition. Here we show that the hydroxyl radical cleavage pattern is a quantitative measure of DNA backbone solvent accessibility, minor groove width, and minor groove electrostatic potential, at single nucleotide resolution. We introduce maps of DNA shape and electrostatic potential as tools for understanding how proteins recognize binding sites in a genome. These maps reveal periodic structural signals in yeast and Drosophila genomic DNA sequences that are associated with positioned nucleosomes.

Local DNA topography correlates with functional noncoding regions of the human genome.

The three-dimensional molecular structure of DNA, specifically the shape of the backbone and grooves of genomic DNA, can be dramatically affected by nucleotide changes, which can cause differences in protein-binding affinity and phenotype. We developed an algorithm to measure constraint on the basis of similarity of DNA topography among multiple species, using hydroxyl radical cleavage patterns to interrogate the solvent-accessible surface area of DNA. This algorithm found that 12% of bases in the human genome are evolutionarily constrained-double the number detected by nucleotide sequence-based algorithms. Topography-informed constrained regions correlated with functional noncoding elements, including enhancers, better than did regions identified solely on the basis of nucleotide sequence. These results support the idea that the molecular shape of DNA is under selection and can identify evolutionary history.

Footprinting protein-DNA complexes using the hydroxyl radical.

Hydroxyl radical footprinting has been widely used for studying the structure of DNA and DNA-protein complexes. The high reactivity and lack of base specificity of the hydroxyl radical makes it an excellent probe for high-resolution footprinting of DNA-protein complexes; this technique can provide structural detail that is not achievable using DNase I footprinting. Hydroxyl radical footprinting experiments can be carried out using readily available and inexpensive reagents and lab equipment. This method involves using the hydroxyl radical to cleave a nucleic acid molecule that is bound to a protein, followed by separating the cleavage products on a denaturing electrophoresis gel to identify the protein-binding sites on the nucleic acid molecule. We describe a protocol for hydroxyl radical footprinting of DNA-protein complexes, along with a troubleshooting guide, that allows researchers to obtain efficient cleavage of DNA in the presence and absence of proteins. This protocol can be completed in 2 d.

Detection of DNA structural motifs in functional genomic elements.

The completion of the human genome project has fueled the search for regulatory elements by a variety of different approaches. Many successful analyses have focused on examining primary DNA sequence and/or chromatin structure. However, it has been difficult to detect common sequence motifs within the feature of chromatin structure most closely associated with regulatory elements, DNase I hypersensitive sites (DHSs). Considering just the nucleotide sequence and/or the chromatin structure of regulatory elements may neglect a critical feature of what is recognized by the regulatory machinery--DNA structure. We introduce a new computational method to detect common DNA structural motifs in a large collection of DHSs that are found in the ENCODE regions of the human genome. We show that DHSs have common DNA structural motifs that show no apparent sequence consensus. One such structural motif is much more highly enriched in experimentally identified DHSs that are in CpG islands and near transcription start sites (TSSs), compared to DHSs not in CpG islands and farther from TSSs, suggesting that DNA structural motifs may participate in the formation of functional regulatory elements. We propose that studies of the conservation of DNA structure, independent of sequence conservation, will provide new information about the link between the nucleotide sequence of a DNA molecule and its experimentally demonstrated function.

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

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.

Effects of discontinuities in the DNA template on abortive initiation and promoter escape by Escherichia coli RNA polymerase.

Using singly gapped or nicked templates containing the T7A1 promoter, we have measured several kinetic parameters related to the process of transcription initiation by Escherichia coli RNA polymerase, confirming and extending previous results using a population of randomly gapped templates. A reduced probability of transcript abortion at RNA lengths of 6 and 7 nucleotides and a lower ratio of abortive to productive initiation events was observed for some discontinuous templates, consistent with models attributing abortive initiation to the accumulation of strain in the initiating complex. The effect of DNA discontinuity on abortion of shorter RNA transcripts (2-3 nucleotides) was less pronounced; abortion at these short chain lengths may primarily be attributed to the low stability of the RNA-DNA hybrid. Certain discontinuities had significant effects on the intrinsic catalytic capacity of the open complex and also on the partitioning between productive and unproductive complexes, suggesting that subtle changes in the conformation of the open complex can profoundly affect its function. The rate and efficiency of promoter escape were not correlated with the stability of the open promoter complex despite previous suggestions to the contrary. We conclude that the stability of the open promoter complex is only one of several factors that contribute to the overall rate of promoter escape.

Platinum anticancer drug damage enforces a particular rotational setting of DNA in nucleosomes.

We constructed two site-specifically modified nucleosomes containing an intrastrand cis-{Pt(NH3)2}2+ 1,3-d(GpTpG) cross-link, similar to one formed by the anticancer drugs carboplatin and cisplatin on DNA, and investigated their structures by hydroxyl radical footprinting and exonuclease III digestion. Hydroxyl radical footprinting demonstrated that the presence of the platinum cross-link selects out a specific rotational setting of DNA on the histone octamer core in each of two reconstituted nucleosomes in which the platinum positions differ by half a DNA helical turn. The {Pt(NH3)2}2+ cross-link is situated in a structurally similar location, with the undamaged strand projecting outward, forcing the DNA to adopt opposite rotational settings in its wrapping around the histone octamer in the two nucleosomes. Enzymatic digestion by exonuclease III of the nucleosome substrates revealed that the platinum cross-link affects the translational positioning of the DNA, forcing it into an asymmetric arrangement with respect to the core histone proteins. We suggest that these phasing phenomena may be central to the recognition and processing of platinum-DNA adducts in cancer cells treated with these drugs and possibly may be common to other DNA damaging events.

Probing DNA structure with hydroxyl radicals.

The hydroxyl radical is a useful probe for studying the shape of the surface of a DNA molecule. Using this technique, fine details of DNA structure can potentially be revealed. This unit describes how to use the hydroxyl radical to generate a random cleavage pattern at the surface of the molecule, separate the broken DNA strands by polyacrylamide gel electrophoresis, and analyze the cleavage pattern to give an image of the surface of the molecule.

Gapped DNA is anisotropically bent.

Ionizing radiation damages DNA in several ways, including through formation of a single-nucleoside gap in one DNA strand. We have developed a two-dimensional gel electrophoresis method to investigate the effect of a strand gap on DNA structure. We generate a library of gapped DNA molecules by treating a DNA restriction fragment with the hydroxyl radical, generated by the reaction of Fe(II) EDTA with hydrogen peroxide. The DNA molecule studied contains a fixed bend produced by a set of phased adenine tracts. The A-tract bend serves as a reference bend for investigating the conformational nature of a strand gap. In the first electrophoretic dimension, a bent DNA molecule that has been treated with the hydroxyl radical is electrophoresed on a native gel. Smearing of the band on the native gel indicates that the library of gapped DNA molecules contains a variety of DNA conformations. In the second electrophoretic dimension, gapped DNA molecules having different native gel mobilities are electrophoresed on separate lanes of a denaturing gel to reveal how each strand gap affects the native gel mobility (and thus shape) of the DNA. Our results demonstrate that a single-nucleoside gap in a DNA duplex leads to an anisotropic, directional bend in the DNA helix axis. The implications of our findings for recognition of this lesion by DNA repair proteins are discussed.

A general synthesis of specifically deuterated nucleotides for studies of DNA and RNA.

An efficient procedure is described for the preparation of ribonucleotides and deoxyribonucleotides with deuterium incorporated at the 1', 4', or 5' position. Three intermediates-[1-2H]-D-ribose, [4-2H]-D-ribose, and [5-2H(2)]-D-ribose-were prepared by chemical synthesis and subsequently converted to ribonucleotides and deoxyribonucleotides via enzymatic reactions. Milligram quantities of the desired products were obtained with an average deuterium content of 96+/-1%.