How sequence defines structure: a crystallographic map of DNA structure and conformation

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.

Directing macromolecular conformation through halogen bonds

The halogen bond, a noncovalent interaction involving polarizable chlorine, bromine, or iodine molecular substituents, is now being exploited to control the assembly of small molecules in the design of supramolecular complexes and new materials. We demonstrate that a halogen bond formed between a brominated uracil and phosphate oxygen can be engineered to direct the conformation of a biological molecule, in this case to define the conformational isomer of a four-stranded DNA junction when placed in direct competition against a classic hydrogen bond. As a result, this bromine interaction is estimated to be approximately 2-5 kcal/mol stronger than the analogous hydrogen bond in this environment, depending on the geometry of the halogen bond. This study helps to establish halogen bonding as a potential tool for the rational design and construction of molecular materials with DNA and other biological macromolecules.

Does topoisomerase II specifically recognize and cleave hairpins, cruciforms and crossovers of DNA?

DNA topoisomerase II is an enzyme that specializes in DNA disentanglement. It catalyzes the interconversion of DNA between different topological states. This event requires the passage of one duplex through another one via a transient double-strand break. Topoisomerase II is able to process any type of DNA, including structures such as DNA juxtapositions (crossovers), DNA hairpins or cruciforms, which are recognized with high specificity. In this review, we focused our attention on topoisomerase II recognizing DNA substrates that possess particular geometries. A strong cleavage site, as we identified in pBR322 DNA in the presence of ellipticine (site 22), appears to be characterized by a cruciform structure formed from two stable hairpins. The same sequence could also constitute a four-way junction structure stabilized by interactions involving ATC sequences. The latter have been shown to be able to promote Holliday junctions. We reviewed the recent literature that deals with the preferential recognition of crossovers by various topoisomerases. The single molecule relaxation experiments have demonstrated the differential abilities of the topoisomerases to recognize crossovers. It appears that enzymes, which distinguish the chirality of the crossovers, possess specialized domains dedicated to this function. We also stress that the formation of crossovers is dependent on the presence of adequate stabilizing sequences. Investigation of the impact of such structures on enzyme activity is important in order to both improve our knowledge of the mechanism of action of the topoisomerase II and to develop new inhibitors of this enzyme.

Programmable oligomers targeting 5'-GGGG-3' in the minor groove of DNA and NF-kappaB binding inhibition

A series of hairpin oligomers containing benzimidazole (Bi) and imidazopyridine (Ip) rings were synthesized and screened to target 5'-WGGGGW-3', a core sequence in the DNA-binding site of NF-kappaB, a prolific transcription factor important in biology and disease. Five Bi and Ip containing oligomers bound to the 5'-WGGGGW-3' site with high affinity. One of the oligomers (Im-Im-Im-Im-gamma-Py-Bi-Py-Bi-beta-Dp) was able to inhibit DNA binding by the transcription factor NF-kappaB.

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND

Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis.

RESULTS

HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.

High guanine and cytosine content increases mRNA levels in mammalian cells

Mammalian genes are highly heterogeneous with respect to their nucleotide composition, but the functional consequences of this heterogeneity are not clear. In the previous studies, weak positive or negative correlations have been found between the silent-site guanine and cytosine (GC) content and expression of mammalian genes. However, previous studies disregarded differences in the genomic context of genes, which could potentially obscure any correlation between GC content and expression. In the present work, we directly compared the expression of GC-rich and GC-poor genes placed in the context of identical promoters and UTR sequences. We performed transient and stable transfections of mammalian cells with GC-rich and GC-poor versions of Hsp70, green fluorescent protein, and IL2 genes. The GC-rich genes were expressed several-fold to over a 100-fold more efficiently than their GC-poor counterparts. This effect was not due to different translation rates of GC-rich and GC-poor mRNA. On the contrary, the efficient expression of GC-rich genes resulted from their increased steady-state mRNA levels. mRNA degradation rates were not correlated with GC content, suggesting that efficient transcription or mRNA processing is responsible for the high expression of GC-rich genes. We conclude that silent-site GC content correlates with gene expression efficiency in mammalian cells.

The effect of nucleotide composition on gene transcription is investigated for Hsp70, GFP, and IL-2, which all show increased expression to correlate with increased GC content at codon position 3.

The stacked-X DNA Holliday junction and protein recognition

The crystal structure of the four-stranded DNA Holliday junction has now been determined in the presence and absence of junction binding proteins, with the extended open-X form of the junction seen in all protein complexes, but the more compact stacked-X structure observed in free DNA. The structures of the stacked-X junction were crystallized because of an unexpected sequence dependence on the stability of this structure. Inverted repeat sequences that contain the general motif NCC or ANC favor formation of stacked-X junctions, with the junction cross-over occurring between the first two positions of the trinucleotides. This review focuses on the sequence dependent structure of the stacked-X junction and how it may play a role in structural recognition by a class of dimeric junction resolving enzymes that themselves show no direct sequence recognition.

Solution formation of Holliday junctions in inverted-repeat DNA sequences

The structure of Holliday junctions has now been well characterized at the atomic level through single-crystal X-ray diffraction in symmetric (inverted-repeat) DNA sequences. At issue, however, is whether the formation of these four-stranded complexes in solution is truly sequence dependent in the manner proposed or is an artifact of the crystallization process and, therefore, has no relevance to the behavior of this central intermediate in homologous recombination and recombination-dependent cellular processes. Here, we apply analytical ultracentrifugation to demonstrate that the sequence d(CCGGTACCGG), which crystallizes in the stacked-X form of the junction, assembles into four-stranded junctions in solution in a manner that is dependent on the DNA and cation concentrations, with an equilibrium established between the junction and duplex forms at 100-200 microM DNA duplex. In contrast, the sequence d(CCGCTAGCGG), which has been crystallized as B-DNA, is seen to adopt only the double-helical form at all DNA and salt concentrations that were tested. Thus, the ACC trinucleotide core is now shown to be important for the formation of Holliday junctions in both crystals and in solution and can be estimated to contribute approximately -4 kcal/mol to stabilizing this recombination intermediate in inverted-repeat sequences.

Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction

We report here the release of a web-based tool (MDDNA) to study and model the fine structural details of DNA on the basis of data extracted from a set of molecular dynamics (MD) trajectories of DNA sequences involving all the unique tetranucleotides. The dynamic web interface can be employed to analyze the first neighbor sequence context effects on the 10 unique dinucleotide steps of DNA. Functionality is included to build all atom models of any user-defined sequence based on the MD results. The backend of this interface is a relational database storing the conformational details of DNA obtained in 39 different MD simulation trajectories comprising all the 136 unique tetranucleotide steps. Examples of the use of this data to predict DNA structures are included.

AVAILABILITY

http://humphry.chem.wesleyan.edu:8080/MDDNA.

SUPPLEMENTARY INFORMATION

Supplementary data including color figures are available at Bioinformatics online.

Study of single-nucleotide polymorphisms by means of electrical conductance measurements

Understanding the complexities of DNA has been a hallmark of science for over a half century, and one of the important topics in DNA research is recognizing the occurrence of mutations in the base-stack. In this article, we present a study of SNPs by direct-contact electrical measurements to a single DNA duplex. We have used short, 11- and 12-bp dsDNA to investigate the change in conductance that occurs if a single base pair, a single base, or two separate bases in the stack are modified. All measurements are carried out in aqueous solution with the DNA chemically bound to the electrodes. These measurements demonstrate that the presence of a single base pair mismatch can be identified by the conductance of the molecule and can cause a change in the conductance of dsDNA by as much as an order of magnitude, depending on the specific details of the double helix and the single nucleotide polymorphism.

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

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.