Nucleosome structural features and intrinsic properties of the TATAAACGCC repeat sequence

HNSPPI: a hybrid computational model combing network and sequence information for predicting protein-protein interaction

Most life activities in organisms are regulated through protein complexes, which are mainly controlled via Protein-Protein Interactions (PPIs). Discovering new interactions between proteins and revealing their biological functions are of great significance for understanding the molecular mechanisms of biological processes and identifying the potential targets in drug discovery. Current experimental methods only capture stable protein interactions, which lead to limited coverage. In addition, expensive cost and time consuming are also the obvious shortcomings. In recent years, various computational methods have been successfully developed for predicting PPIs based only on protein homology, primary sequences of protein or gene ontology information. Computational efficiency and data complexity are still the main bottlenecks for the algorithm generalization. In this study, we proposed a novel computational framework, HNSPPI, to predict PPIs. As a hybrid supervised learning model, HNSPPI comprehensively characterizes the intrinsic relationship between two proteins by integrating amino acid sequence information and connection properties of PPI network. The experimental results show that HNSPPI works very well on six benchmark datasets. Moreover, the comparison analysis proved that our model significantly outperforms other five existing algorithms. Finally, we used the HNSPPI model to explore the SARS-CoV-2-Human interaction system and found several potential regulations. In summary, HNSPPI is a promising model for predicting new protein interactions from known PPI data.

Deciphering the mechanical code of the genome and epigenome

Diverse DNA-deforming processes are impacted by the local mechanical and structural properties of DNA, which in turn depend on local sequence and epigenetic modifications. Deciphering this mechanical code (that is, this dependence) has been challenging due to the lack of high-throughput experimental methods. Here we present a comprehensive characterization of the mechanical code. Utilizing high-throughput measurements of DNA bendability via loop-seq, we quantitatively established how the occurrence and spatial distribution of dinucleotides, tetranucleotides and methylated CpG impact DNA bendability. We used our measurements to develop a physical model for the sequence and methylation dependence of DNA bendability. We validated the model by performing loop-seq on mouse genomic sequences around transcription start sites and CTCF-binding sites. We applied our model to test the predictions of all-atom molecular dynamics simulations and to demonstrate that sequence and epigenetic modifications can mechanically encode regulatory information in diverse contexts.

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

DNA mechanics and its biological impact

Almost all nucleoprotein interactions and DNA manipulation events involve mechanical deformations of DNA. Extraordinary progresses in single-molecule, structural, and computational methods have characterized the average mechanical properties of DNA, such as bendability and torsional rigidity, in high resolution. Further, the advent of sequencing technology has permitted measuring, in high-throughput, how such mechanical properties vary with sequence and epigenetic modifications along genomes. We review these recent technological advancements, and discuss how they have contributed to the emerging idea that variations in the mechanical properties of DNA play a fundamental role in regulating, genome-wide, diverse processes involved in chromatin organization.

Measuring DNA mechanics on the genome scale

Mechanical deformations of DNA such as bending are ubiquitous and implicated in diverse cellular functions¹. However, the lack of high-throughput tools to directly measure the mechanical properties of DNA limits our understanding of whether and how DNA sequences modulate DNA mechanics and associated chromatin transactions genome-wide. We developed an assay called loop-seq to measure the intrinsic cyclizability of DNA – a proxy for DNA bendability – in high throughput. We measured the intrinsic cyclizabilities of 270,806 50 bp DNA fragments that span the entire length of S. cerevisiae chromosome V and other genomic regions, and also include random sequences. We discovered sequence-encoded regions of unusually low bendability upstream of Transcription Start Sites (TSSs). These regions disfavor the sharp DNA bending required for nucleosome formation and are co-centric with known Nucleosome Depleted Regions (NDRs). We show biochemically that low bendability of linker DNA located about 40 bp away from a nucleosome edge inhibits nucleosome sliding into the linker by the chromatin remodeler INO80. The observation explains how INO80 can create promoter-proximal nucleosomal arrays in the absence of any other factors² by reading the DNA mechanical landscape. We show that chromosome wide, nucleosomes are characterized by high DNA bendability near dyads and low bendability near the linkers. This contrast increases for nucleosomes deeper into gene bodies, suggesting that DNA mechanics plays a previously unappreciated role in organizing nucleosomes far from the TSS, where nucleosome remodelers predominate. Importantly, random substitution of synonymous codons does not preserve this contrast, suggesting that the evolution of codon choice has been impacted by selective pressure to preserve sequence-encoded mechanical modulations along genes. We also provide evidence that transcription through the TSS-proximal nucleosomes is impacted by local DNA mechanics. Overall, this first genome-scale map of DNA mechanics hints at a ‘mechanical code’ with broad functional implications.

Effects of Psoralen on Histone-DNA Interactions Studied by Using Atomic Force Microscopy

The investigation of the DNA-histone interactions and factors that affect such interactions in the nucleosome is essential for understanding the role of chromatin organization in all cellular processes involved in the repair, transcription, and replication of the eukaryotic genome. As a kind of photosensitive molecule, psoralen (PSO) is used in the treatment of skin disease with ultraviolet light (PSO and ultra violet light, type A). The effect of treatment is remarkable, but the side effect is also obvious. PSO can be embedded in a 5' TA sequence in double-stranded DNA (dsDNA), and dsDNA is mainly wrapped around a histone octamer to form a nucleosome structure in human cells. Therefore, it is very necessary to explore the influence of PSO on DNA-histone interactions. To this end, the binding specificity and mode of DNA and histone in the presence or absence of PSO are investigated systematically. The results show that the presence of PSO (no matter if there is ultra violet light treatment) can increase the overall probability of histone binding to dsDNA while lowering the selectivity of histone binding to the specific DNA sequence in vitro. In addition, the increase of solution ionic strength can lower the ratio of histone binding to nonspecific DNA.

LeNup: learning nucleosome positioning from DNA sequences with improved convolutional neural networks

Motivation

Nucleosome positioning plays significant roles in proper genome packing and its accessibility to execute transcription regulation. Despite a multitude of nucleosome positioning resources available on line including experimental datasets of genome-wide nucleosome occupancy profiles and computational tools to the analysis on these data, the complex language of eukaryotic Nucleosome positioning remains incompletely understood.

Results

Here, we address this challenge using an approach based on a state-of-the-art machine learning method. We present a novel convolutional neural network (CNN) to understand nucleosome positioning. We combined Inception-like networks with a gating mechanism for the response of multiple patterns and long term association in DNA sequences. We developed the open-source package LeNup based on the CNN to predict nucleosome positioning in Homo sapiens, Caenorhabditis elegans, Drosophila melanogaster as well as Saccharomyces cerevisiae genomes. We trained LeNup on four benchmark datasets. LeNup achieved greater predictive accuracy than previously published methods.

Availability and implementation

LeNup is freely available as Python and Lua script source code under a BSD style license from https://github.com/biomedBit/LeNup.

Supplementary information

Supplementary data are available at Bioinformatics online.

Hydroxyl-radical footprinting combined with molecular modeling identifies unique features of DNA conformation and nucleosome positioning

Nucleosomes are the most abundant protein–DNA complexes in eukaryotes that provide compaction of genomic DNA and are implicated in regulation of transcription, DNA replication and repair. The details of DNA positioning on the nucleosome and the DNA conformation can provide key regulatory signals. Hydroxyl-radical footprinting (HRF) of protein–DNA complexes is a chemical technique that probes nucleosome organization in solution with a high precision unattainable by other methods. In this work we propose an integrative modeling method for constructing high-resolution atomistic models of nucleosomes based on HRF experiments. Our method precisely identifies DNA positioning on nucleosome by combining HRF data for both DNA strands with the pseudo-symmetry constraints. We performed high-resolution HRF for Saccharomyces cerevisiae centromeric nucleosome of unknown structure and characterized it using our integrative modeling approach. Our model provides the basis for further understanding the cooperative engagement and interplay between Cse4p protein and the A-tracts important for centromere function.

The 1-Particle-per-k-Nucleotides (1PkN) Elastic Network Model of DNA Dynamics with Sequence-Dependent Geometry

Coarse-grained models of DNA have made important contributions to the determination of the physical properties of genomic DNA, working as a molecular machine for gene regulation. In this study, to analyze the global dynamics of long DNA sequences with consideration of sequence-dependent geometry, we propose elastic network models of DNA where each particle represents k nucleotides (1-particle-per-k-nucleotides, 1PkN). The models were adjusted according to profiles of the anisotropic fluctuations obtained from our previous 1-particle-per-1-nucleotide (1P1N) model, which was proven to reproduce such profiles of all-atom models. We confirmed that the 1P3N and 1P4N models are suitable for the analysis of detailed dynamics such as local twisting motion. The models are intended for the analysis of large structures, e.g., 10-nm fibers in the nucleus, and nucleoids of mitochondrial or phage DNA at low computational costs. As an example, we surveyed the physical characteristics of the whole mitochondrial human and Plasmodium falciparum genomes.

DNA Structural Correlation in Short and Long Ranges

Recent single-molecule measurements have revealed the DNA allostery in protein/DNA binding. MD simulations showed that this allosteric effect is associated with the deformation properties of DNA. In this study, we used MD simulations to further investigate the mechanism of DNA structural correlation, its dependence on DNA sequence, and the chemical modification of the bases. Besides a random sequence, poly d(AT) and poly d(GC) are also used as simpler model systems, which show the different bending and twisting flexibilities. The base-stacking interactions and the methyl group on the 5-carbon site of thymine causes local structures and flexibility to be very different for the two model systems, which further lead to obviously different tendencies of the conformational deformations, including the long-range allosteric effects.

iNuc-PseKNC: a sequence-based predictor for predicting nucleosome positioning in genomes with pseudo k-tuple nucleotide composition

MOTIVATION

Nucleosome positioning participates in many cellular activities and plays significant roles in regulating cellular processes. With the avalanche of genome sequences generated in the post-genomic age, it is highly desired to develop automated methods for rapidly and effectively identifying nucleosome positioning. Although some computational methods were proposed, most of them were species specific and neglected the intrinsic local structural properties that might play important roles in determining the nucleosome positioning on a DNA sequence.

RESULTS

Here a predictor called 'iNuc-PseKNC' was developed for predicting nucleosome positioning in Homo sapiens, Caenorhabditis elegans and Drosophila melanogaster genomes, respectively. In the new predictor, the samples of DNA sequences were formulated by a novel feature-vector called 'pseudo k-tuple nucleotide composition', into which six DNA local structural properties were incorporated. It was observed by the rigorous cross-validation tests on the three stringent benchmark datasets that the overall success rates achieved by iNuc-PseKNC in predicting the nucleosome positioning of the aforementioned three genomes were 86.27%, 86.90% and 79.97%, respectively. Meanwhile, the results obtained by iNuc-PseKNC on various benchmark datasets used by the previous investigators for different genomes also indicated that the current predictor remarkably outperformed its counterparts.

AVAILABILITY

A user-friendly web-server, iNuc-PseKNC is freely accessible at http://lin.uestc.edu.cn/server/iNuc-PseKNC.

Poly(dA:dT)-rich DNAs are highly flexible in the context of DNA looping

Large-scale DNA deformation is ubiquitous in transcriptional regulation in prokaryotes and eukaryotes alike. Though much is known about how transcription factors and constellations of binding sites dictate where and how gene regulation will occur, less is known about the role played by the intervening DNA. In this work we explore the effect of sequence flexibility on transcription factor-mediated DNA looping, by drawing on sequences identified in nucleosome formation and ligase-mediated cyclization assays as being especially favorable for or resistant to large deformations. We examine a poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong to a class of highly inflexible DNAs; two strong nucleosome positioning sequences that share a set of particular sequence features common to nucleosome-preferring DNAs; and a CG-rich sequence representative of high G+C-content genomic regions that correlate with high nucleosome occupancy in vivo. To measure the flexibility of these sequences in the context of DNA looping, we combine the in vitro single-molecule tethered particle motion assay, a canonical looping protein, and a statistical mechanical model that allows us to quantitatively relate the looping probability to the looping free energy. We show that, in contrast to the case of nucleosome occupancy, G+C content does not positively correlate with looping probability, and that despite sharing sequence features that are thought to determine nucleosome affinity, the two strong nucleosome positioning sequences behave markedly dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich DNA that is often characterized as highly inflexible in fact exhibits one of the highest propensities for looping that we have measured. These results argue for a need to revisit our understanding of the mechanical properties of DNA in a way that will provide a basis for understanding DNA deformation over the entire range of biologically relevant scenarios that are impacted by DNA deformability.

Nucleosome structural changes induced by binding of non-histone chromosomal proteins HMGN1 and HMGN2

Interactions between the nucleosome and the non-histone chromosomal proteins (HMGN1 and HMGN2) were studied by circular dichroism (CD) spectroscopy to elucidate structural changes in the nucleosome induced by HMGN binding. Unlike previous studies that used a nucleosome extracted from living cells, in this study we utilized a nucleosome reconstituted from unmodified recombinant histones synthesized in Escherichia coli and a 189-bp synthetic DNA fragment harboring a nucleosome positioning sequence. This DNA fragment consists of 5′-TATAAACGCC-3′ repeats that has a high affinity to the histone octamer. A nucleosome containing a unique octamer-binding sequence at a specific location on the DNA was produced at sufficiently high yield for spectroscopic analysis. CD data have indicated that both HMGN1 and HMGN2 can increase the winding angle of the nucleosome DNA, but the extent of the structural changes induced by these proteins differs significantly. This suggests HMGN1 and HMGN2 would have different abilities to facilitate nucleosome remodeling.

•

A nucleosome was reconstituted from recombinant histones and a synthetic DNA.

•

Nucleosomes were produced at sufficiently high yield for spectroscopic analysis.

•

A nucleosome with and without HMGN proteins was analyzed using CD spectroscopy.

•

CD data indicate that HMGN proteins increase the winding angle of the nucleosome DNA.

•

HMGN1 and HMGN2 may have different abilities to facilitate nucleosome remodeling.

Protein sliding and DNA denaturation are essential for DNA organization by human mitochondrial transcription factor A

Mitochondria organize their genome in protein-DNA complexes called nucleoids. The mitochondrial transcription factor A (TFAM), a protein that regulates mitochondrial transcription, is abundant in these nucleoids. TFAM is believed to be essential for mitochondrial DNA compaction, yet the exact mechanism has not been resolved. Here we use a combination of single-molecule manipulation and fluorescence microscopy to show the nonspecific DNA-binding dynamics and compaction by TFAM. We observe that single TFAM proteins diffuse extensively over DNA (sliding) and, by collisions, form patches on DNA in a cooperative manner. Moreover, we demonstrate that TFAM induces compaction by changing the flexibility of the DNA, which can be explained by local denaturation of the DNA (melting). Both sliding of TFAM and DNA melting are also necessary characteristics for effective, specific transcription regulation by TFAM. This apparent connection between transcription and DNA organization clarifies how TFAM can accomplish two complementary roles in the mitochondrial nucleoid at the same time.

Predicting nucleosome positions in yeast: using the absolute frequency

Nucleosome is the basic structure of chromatin in eukaryotic cells, and they form the chromatin fiber interconnected by sections of linker DNA. Nucleosome positioning is of great significance for gene transcription regulation. In this paper, we consider the difference of absolute frequency of nucleotides between the nucleosome forming and nucleosome inhibiting sequences. Based on the 2-mer absolute frequency of nucleotides in genome, a new model is constructed to distinguish nucleosome DNA and linker DNA. When used to predict DNA potential for forming nucleosomes in S. cerevisiae, the model achieved a high accuracy of 96.05%. Thus, the model is very useful for predicting nucleosome positioning.

The divalent cations Ca2+ and Mg2+ play specific roles in stabilizing histone-DNA interactions within nucleosomes that are partially redundant with the core histone tail domains

We previously reported that reconstituted nucleosomes undergo sequence-dependent translational repositioning upon removal of the core histone tail domains under physiological conditions, indicating that the tails influence the choice of position. We report here that removal of the core histone tail domains increases the exposure of the DNA backbone in nucleosomes to hydroxyl radicals, a nonbiased chemical cleavage reagent, indicative of an increase in the motility of the DNA on the histone surface. Moreover, we demonstrate that the divalent cations Mg(2+) and Ca(2+) can replace the role of the tail domains with regard to stabilization of histone-DNA interactions within the nucleosome core and restrict repositioning of nucleosomes upon tail removal. However, when nucleosomes were incubated with Mg(2+) after tail removal, the original distribution of translational positions was not re-established, indicating that divalent cations increase the energy barrier between translational positions rather than altering the free energy differences between positions. Interestingly, other divalent cations such as Zn(2+), Fe(2+), Co(2+), and Mn(2+) had little or no effect on the stability of histone-DNA interactions within tailless nucleosomes. These results support the idea that specific binding sites for Mg(2+) and Ca(2+) ions exist within the nucleosome and play a critical role in nucleosome stability that is partially redundant with the core histone tail domains.

Human epigenome data reveal increased CpG methylation in alternatively spliced sites and putative exonic splicing enhancers

The role of gene body methylation, which represents a major part of methylation in DNA, remains mostly unknown. Evidence based on the CpG distribution associates its presence with nucleosome positioning and alternative splicing. Recently, it was also shown that cytosine methylation influences splicing. However, to date, there is no methylation-based data on the association of methylation with alternative splicing and the distribution in exonic splicing enhancers (ESEs). We presently report that, based on the computational analysis of the Human Epigenome Project data, CpG hypermethylation (>80%) is frequent in alternatively spliced sites (particularly in noncanonical) but not in alternate promoters. The methylation frequency increases in sequences containing multiple putative ESEs. However, significant differences in the extent of methylation are observed among different ESEs. Specifically, moderate levels of methylation, ranging from 20% to 80%, are frequent in SRp55-binding elements, which are associated with response to extracellular conditions, but not in SF2/ASF, primarily responsible for alternative splicing, or in CpG islands. Finally, methylation is more frequent in the presence of AT repeats and CpGs separated by 10 nucleotides and lower in adjacent CpGs, probably indicating its dependence on helical formations and on the presence of nucleosome positioning-related sequences. In conclusion, our results show the regulation of methylation in ESEs and support its involvement in alternative splicing.

Prediction of nucleosome DNA formation potential and nucleosome positioning using increment of diversity combined with quadratic discriminant analysis

In this work, a novel method was developed to distinguish nucleosome DNA and linker DNA based on increment of diversity combined with quadratic discriminant analysis (IDQD), using k-mer frequency of nucleotides in genome. When used to predict DNA potential for forming nucleosomes, the model achieved a high accuracy of 94.94%, 77.60%, and 86.81%, respectively, for Saccharomyces cerevisiae, Homo sapiens, and Drosophila melanogaster. The area under the receiver operator characteristics curve of our classifier was 0.982 for S. cerevisiae. Our results indicate that DNA sequence preference is critical for nucleosome formation potential and is likely conserved across eukaryotes. The model successfully identified nucleosome-enriched or nucleosome-depleted regions in S. cerevisiae genome, suggesting nucleosome positioning depends on DNA sequence preference. Thus, IDQD classifier is useful for predicting nucleosome positioning.

Crystal structures of nucleosome core particles containing the '601' strong positioning sequence

Nucleosome positioning plays a key role in genomic regulation by defining histone-DNA context and by modulating access to specific sites. Moreover, the histone-DNA register influences the double-helix structure, which in turn can affect the association of small molecules and protein factors. Analysis of genomic and synthetic DNA has revealed sequence motifs that direct nucleosome positioning in vitro; thus, establishing the basis for the DNA sequence dependence of positioning would shed light on the mechanics of the double helix and its contribution to chromatin structure in vivo. However, acquisition of well-diffracting nucleosome core particle (NCP) crystals is extremely dependent on the DNA fragment used for assembly, and all previous NCP crystal structures have been based on human α-satellite sequences. Here, we describe the crystal structures of Xenopus NCPs containing one of the strongest known histone octamer binding and positioning sequences, the so-called '601' DNA. Two distinct 145-bp 601 crystal forms display the same histone-DNA register, which coincides with the occurrence of DNA stretching-overtwisting in both halves of the particle around five double-helical turns from the nucleosome center, giving the DNA an 'effective length' of 147 bp. As we have found previously with stretching around two turns from the nucleosome center for a centromere-based sequence, the terminal stretching observed in the 601 constructs is associated with extreme kinking into the minor groove at purine-purine (pyrimidine-pyrimidine) dinucleotide steps. In other contexts, these step types display an overall nonflexible behavior, which raises the possibility that DNA stretching in the nucleosome or extreme distortions in general have unique sequence dependency characteristics. Our findings indicate that DNA stretching is an intrinsically predisposed site-specific property of the nucleosome and suggest how NCP crystal structures with diverse DNA sequences can be obtained.

Structural insight into the sequence dependence of nucleosome positioning

Nucleosome positioning displays sequence dependency and contributes to genomic regulation in a site-specific manner. We solved the structures of nucleosome core particle composed of strong positioning TTTAA elements flanking the nucleosome center. The positioning strength of the super flexible TA dinucleotide is consistent with its observed central location within minor groove inward regions, where it can contribute maximally to energetically challenging minor groove bending, kinking and compression. The marked preference for TTTAA and positioning power of the site 1.5 double helix turns from the nucleosome center relates to a unique histone protein motif at this location, which enforces a sustained, extremely narrow minor groove via a hydrophobic "sugar clamp." Our analysis sheds light on the basis of nucleosome positioning and indicates that the histone octamer has evolved not to fully minimize sequence discrimination in DNA binding.

Computational prediction of nucleosome positioning by calculating the relative fragment frequency index of nucleosomal sequences

We developed an accurate method to predict nucleosome positioning from genome sequences by refining the previously developed method of Peckham et al. (2007). Here, we used the relative fragment frequency index we developed and a support vector machine to screen for nucleosomal and linker DNA sequences. Our twofold cross-validation revealed that the accuracy of our method based on the area under the receiver operating characteristic curve was 81%, whereas that of Peckham's method was 75% when both of two nucleosomal sequence data obtained from independent experiments were used for validation. We suggest that our method is more effective in predicting nucleosome positioning.

A comparative approach for the investigation of biological information processing: an examination of the structure and function of computer hard drives and DNA

BACKGROUND

The robust storage, updating and utilization of information are necessary for the maintenance and perpetuation of dynamic systems. These systems can exist as constructs of metal-oxide semiconductors and silicon, as in a digital computer, or in the "wetware" of organic compounds, proteins and nucleic acids that make up biological organisms. We propose that there are essential functional properties of centralized information-processing systems; for digital computers these properties reside in the computer's hard drive, and for eukaryotic cells they are manifest in the DNA and associated structures.

METHODS

Presented herein is a descriptive framework that compares DNA and its associated proteins and sub-nuclear structure with the structure and function of the computer hard drive. We identify four essential properties of information for a centralized storage and processing system: (1) orthogonal uniqueness, (2) low level formatting, (3) high level formatting and (4) translation of stored to usable form. The corresponding aspects of the DNA complex and a computer hard drive are categorized using this classification. This is intended to demonstrate a functional equivalence between the components of the two systems, and thus the systems themselves.

RESULTS

Both the DNA complex and the computer hard drive contain components that fulfill the essential properties of a centralized information storage and processing system. The functional equivalence of these components provides insight into both the design process of engineered systems and the evolved solutions addressing similar system requirements. However, there are points where the comparison breaks down, particularly when there are externally imposed information-organizing structures on the computer hard drive. A specific example of this is the imposition of the File Allocation Table (FAT) during high level formatting of the computer hard drive and the subsequent loading of an operating system (OS). Biological systems do not have an external source for a map of their stored information or for an operational instruction set; rather, they must contain an organizational template conserved within their intra-nuclear architecture that "manipulates" the laws of chemistry and physics into a highly robust instruction set. We propose that the epigenetic structure of the intra-nuclear environment and the non-coding RNA may play the roles of a Biological File Allocation Table (BFAT) and biological operating system (Bio-OS) in eukaryotic cells.

CONCLUSIONS

The comparison of functional and structural characteristics of the DNA complex and the computer hard drive leads to a new descriptive paradigm that identifies the DNA as a dynamic storage system of biological information. This system is embodied in an autonomous operating system that inductively follows organizational structures, data hierarchy and executable operations that are well understood in the computer science industry. Characterizing the "DNA hard drive" in this fashion can lead to insights arising from discrepancies in the descriptive framework, particularly with respect to positing the role of epigenetic processes in an information-processing context. Further expansions arising from this comparison include the view of cells as parallel computing machines and a new approach towards characterizing cellular control systems.

The most frequent short sequences in non-coding DNA

The purpose of this work is to determine the most frequent short sequences in non-coding DNA. They may play a role in maintaining the structure and function of eukaryotic chromosomes. We present a simple method for the detection and analysis of such sequences in several genomes, including Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens. We also study two chromosomes of man and mouse with a length similar to the whole genomes of the other species. We provide a list of the most common sequences of 9–14 bases in each genome. As expected, they are present in human Alu sequences. Our programs may also give a graph and a list of their position in the genome. Detection of clusters is also possible. In most cases, these sequences contain few alternating regions. Their intrinsic structure and their influence on nucleosome formation are not known. In particular, we have found new features of short sequences in C. elegans, which are distributed in heterogeneous clusters. They appear as punctuation marks in the chromosomes. Such clusters are not found in either A. thaliana or D. melanogaster. We discuss the possibility that they play a role in centromere function and homolog recognition in meiosis.

Evidence that nucleosomes inhibit mismatch repair in eukaryotic cells

The influence of chromatin structure on DNA metabolic processes, including DNA replication and repair, has been a matter of intensive studies in recent years. Although the human mismatch repair (MMR) reaction has been reconstituted using purified proteins, the influence of chromatin structure on human MMR is unknown. This study examines the interaction between human MutSalpha and a mismatch located within a nucleosome or between two nucleosomes. The results show that, whereas MutSalpha specifically recognizes both types of nucleosomal heteroduplexes, the protein bound the mismatch within a nucleosome with much lower efficiency than a naked heteroduplex or a heterology free of histone proteins but between two nucleosomes. Additionally, MutSalpha displays reduced ATPase- and ADP-binding activity when interacting with nucleosomal heteroduplexes. Interestingly, nucleosomes block ATP-induced MutSalpha sliding along the DNA helix when the mismatch is in between two nucleosomes. These findings suggest that nucleosomes may inhibit MMR in eukaryotic cells. The implications of these findings for our understanding of eukaryotic MMR are discussed.

Sequence-specific physical properties of African green monkey alpha-satellite DNA contribute to centromeric heterochromatin formation

Satellite DNA, a major component of eukaryotic centromeric heterochromatin, is potentially associated with the processes ensuring the faithful segregation of the genetic material during cell division. Structural properties of alpha-satellite DNA (AS) from African green monkey (AGM) were studied. Atomic force microscopy imaging showed smaller end-to-end distances of AS fragments than would be expected for the persistence length of random sequence DNA. The apparent persistence length of the AS was determined as 35nm. Gel-electrophoresis indicated only a weak contribution of intrinsic curvature to the DNA conformations suggesting an additional contribution of an elevated bending flexibility to the reduced end-to-end distances. Next, the force-extension behavior of the naked AS and in complex with nucleosomes was studied using optical tweezers. The naked AS showed a reduced overstretching transition force (-18% the value determined for random DNA) and higher forces required to straighten the DNA. Finally, reconstituted AS nucleosomes disrupted at significantly higher forces as compared with random DNA nucleosomes which is probably due to structural properties of the AS which stabilize the nucleosomes. The data support that the AS plays a role in the formation of centromeric heterochromatin due to specific structural properties and suggest that a relatively higher mechanical stability of nucleosomes is important in AGM-AS chromatin.

Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (approximately 10.3 bp) at the two ends, with larger (approximately 11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames' start and end sites were occupied by nucleosomes.

Genome-wide analysis of T-DNA integration into the chromosomes of Magnaporthe oryzae

Agrobacterium tumefaciens-mediated transformation (ATMT) has become a prevalent tool for functional genomics of fungi, but our understanding of T-DNA integration into the fungal genome remains limited relative to that in plants. Using a model plant-pathogenic fungus, Magnaporthe oryzae, here we report the most comprehensive analysis of T-DNA integration events in fungi and the development of an informatics infrastructure, termed a T-DNA analysis platform (TAP). We identified a total of 1110 T-DNA-tagged locations (TTLs) and processed the resulting data via TAP. Analysis of the TTLs showed that T-DNA integration was biased among chromosomes and preferred the promoter region of genes. In addition, irregular patterns of T-DNA integration, such as chromosomal rearrangement and readthrough of plasmid vectors, were also observed, showing that T-DNA integration patterns into the fungal genome are as diverse as those of their plant counterparts. However, overall the observed junction structures between T-DNA borders and flanking genomic DNA sequences revealed that T-DNA integration into the fungal genome was more canonical than those observed in plants. Our results support the potential of ATMT as a tool for functional genomics of fungi and show that the TAP is an effective informatics platform for handling data from large-scale insertional mutagenesis.

Nucleosome Positioning Determinants

A previous report demonstrated that one site in a nucleosome assembled onto a synthetic positioning sequence known as Fragment 67 is hypersensitive to permanganate. The site is required for positioning activity and is located 1.5 turns from the dyad, which is a region of high DNA curvature in the nucleosome. Here, the permanganate sensitivity of the nucleosome positioning Fragment 601 was examined in order to expand the dataset of nucleosome sequences containing KMnO(4) hypersensitive sites. The hyperreactive T residue in the six sites detected as well as the one in Fragment 67 and three in the 5 S rDNA positioning sequence were contained within a TA step. Seven of the ten sequences were of the form CTAGPuG or the related sequence TTAAPu. These motifs were also found in the binding sites of several transcriptional regulatory proteins that kink DNA. In order to assess the significance of these sites, the 10 bp positioning determinant in Fragment 67 was removed and replaced by the nine sequences from the 5 S rDNA and Fragment 601. The results demonstrated that these derivative fragments promoted high nucleosome stability and positioning as compared to a control sequence that contained an AT step in place of the TA step. The importance of the TA step was further tested by making single base-pair substitutions in Fragment 67 and the results revealed that stability and positioning activity followed the order: TA>TG>TT>/=TC approximately GG approximately GA approximately AT. Sequences flanking the TA step were also shown to be critical for nucleosome stability and positioning. Nucleosome positioning was restored to near wild-type levels with (CTG)(3), which can form slipped stranded structures and with one base bulges that kink DNA. The results of this study suggest that local DNA structures are important for positioning and that single base-pair changes at these sites could have profound effects on those genomic functions that depend on ordered nucleosomes.

Nucleosome positioning signals in genomic DNA

Although histones can form nucleosomes on virtually any genomic sequence, DNA sequences show considerable variability in their binding affinity. We have used DNA sequences of Saccharomyces cerevisiae whose nucleosome binding affinities have been experimentally determined (Yuan et al. 2005) to train a support vector machine to identify the nucleosome formation potential of any given sequence of DNA. The DNA sequences whose nucleosome formation potential are most accurately predicted are those that contain strong nucleosome forming or inhibiting signals and are found within nucleosome length stretches of genomic DNA with continuous nucleosome formation or inhibition signals. We have accurately predicted the experimentally determined nucleosome positions across a well-characterized promoter region of S. cerevisiae and identified strong periodicity within 199 center-aligned mononucleosomes studied recently (Segal et al. 2006) despite there being no periodicity information used to train the support vector machine. Our analysis suggests that only a subset of nucleosomes are likely to be positioned by intrinsic sequence signals. This observation is consistent with the available experimental data and is inconsistent with the proposal of a nucleosome positioning code. Finally, we show that intrinsic nucleosome positioning signals are both more inhibitory and more variable in promoter regions than in open reading frames in S. cerevisiae.

The Core Histone Tail Domains Contribute to Sequence-dependent Nucleosome Positioning

The precise positioning of nucleosomes plays a critical role in the regulation of gene expression by modulating the DNA binding activity of trans-acting factors. However, molecular determinants responsible for positioning are not well understood. We examined whether the removal of the core histone tail domains from nucleosomes reconstituted with specific DNA fragments led to alteration of translational positions. Remarkably, we find that removal of tail domains from a nucleosome assembled on a DNA fragment containing a Xenopus borealis somatic-type 5S RNA gene results in repositioning of nucleosomes along the DNA, including two related major translational positions that move about 20 bp further upstream with respect to the 5S gene. In a nucleosome reconstituted with a DNA fragment containing the promoter of a Drosophila alcohol dehydrogenase gene, several translational positions shifted by about 10 bp along the DNA upon tail removal. However, the positions of nucleosomes assembled with a DNA fragment known to have one of the highest binding affinities for core histone proteins in the mouse genome were not altered by removal of core histone tail domains. Our data support the notion that the basic tail domains bind to nucleosomal DNA and influence the selection of the translational position of nucleosomes and that once tails are removed movement between translational positions occurs in a facile manner on some sequences. However, the effect of the N-terminal tails on the positioning and movement of a nucleosome appears to be dependent on the DNA sequence such that the contribution of the tails can be masked by very high affinity DNA sequences. Our results suggest a mechanism whereby sequence-dependent nucleosome positioning can be specifically altered by regulated changes in histone tail-DNA interactions in chromatin.

A whole genome analysis of 5' regulatory regions of human genes for putative cis-acting modulators of nucleosome positioning

We have carried out in silico analysis of upstream regions of 23,034 genes from the human genome for sequence motifs, which can potentially affect nucleosome positioning. Nucleosome exclusion elements (NEE) occur in 12% of the genes while less than 1% contain nucleosome positioning elements (NPE). NEE are significantly higher in 5' regions of certain categories of genes, namely, genes with active promoters, genes localised to gene-rich chromosomes 16, 17 and 19, genes having significantly higher expression levels and higher levels of occupancy of general transcription machinery proteins. NEE are also enriched in housekeeping and TATA-less genes, but are significantly under-represented in the upstream region of genes functionally classified under 'organ development' and 'morphogenesis' categories. Further, DNase I hypersensitive sites which co-localise with NEE, preferentially occur in 5' regulatory regions. Considering the positioning sequences identified so far, we speculate that low affinity nucleosome positioning in the upstream sequences of genes in the human genome is the default state requiring activation through chromatin remodelling, while, there appears to be a selection for nucleosome excluding sequences in the upstream sequences of genes that are ubiquitously expressed.

Why repetitive DNA is essential to genome function

There are clear theoretical reasons and many well-documented examples which show that repetitive, DNA is essential for genome function. Generic repeated signals in the DNA are necessary to format expression of unique coding sequence files and to organise additional functions essential for genome replication and accurate transmission to progeny cells. Repetitive DNA sequence elements are also fundamental to the cooperative molecular interactions forming nucleoprotein complexes. Here, we review the surprising abundance of repetitive DNA in many genomes, describe its structural diversity, and discuss dozens of cases where the functional importance of repetitive elements has been studied in molecular detail. In particular, the fact that repeat elements serve either as initiators or boundaries for heterochromatin domains and provide a significant fraction of scaffolding/matrix attachment regions (S/MARs) suggests that the repetitive component of the genome plays a major architectonic role in higher order physical structuring. Employing an information science model, the 'functionalist' perspective on repetitive DNA leads to new ways of thinking about the systemic organisation of cellular genomes and provides several novel possibilities involving repeat elements in evolutionarily significant genome reorganisation. These ideas may facilitate the interpretation of comparisons between sequenced genomes, where the repetitive DNA component is often greater than the coding sequence component.

Large scale preparation of nucleosomes containing site-specifically chemically modified histones lacking the core histone tail domains

The core histone tail domains are critical regulators of chromatin structure and function and modifications such as acetylation of lysine residues within the tails are central to this regulation. Studies have shown that the removal of core histone tail domains by trypsinization in which one-half to two-thirds of each core histone tail domain are removed in gross aspects mimics the acetylation of core histone tails. In addition, removal of the tails has been useful in understanding general tail function. Thus, removal of native core histone tails by trypsinization is a widely used method. In addition, many in vitro studies now employ core histones site-specifically modified with photo activatable cross-linking probes or fluorescent probes. However, in our experience, standard methods employing trypsinized donor chromatin for reconstitution of nucleosomes containing certain chemically modified histones lacking the core histone tail domains are not uniformly applicable. Here, we describe various methods for preparing nucleosomes containing a core histone modified with a cross-linking agent, APB, and lacking the core histone tail domains.

Structural properties of promoters: similarities and differences between prokaryotes and eukaryotes

During the process of transcription, RNA polymerase can exactly locate a promoter sequence in the complex maze of a genome. Several experimental studies and computational analyses have shown that the promoter sequences apparently possess some special properties, such as unusual DNA structures and low stability, which make them distinct from the rest of the genome. But most of these studies have been carried out on a particular set of promoter sequences or on promoter sequences from similar organisms. To examine whether the promoters from a wide variety of organisms share these special properties, we have carried out an analysis of sets of promoters from bacteria, vertebrates and plants. These promoters were analyzed with respect to the prediction of three different properties, such as DNA curvature, bendability and stability, which are relevant to transcription. All the promoter sequences are predicted to share certain features, such as stability and bendability profiles, but there are significant differences in DNA curvature profiles and nucleotide composition between the different organisms. These similarities and differences are correlated with some of the known facts about transcription process in the promoters from the three groups of organisms.

The influence of DNA stiffness upon nucleosome formation

The rotational and translational positioning of nucleosomes on DNA is dependent to a significant extent on the physicochemical properties of the double helix. We have investigated the influence of the axial flexibility of the molecule on the affinity for the histone octamer by substituting selected DNA sequences with either inosine for guanosine or diaminopurine for adenine. These substitutions, respectively, remove or add a purine 2-amino group exposed in the minor groove and, respectively, decrease and increase the apparent persistence length. We observe that for all sequences tested inosine substitution, with one exception, increases the affinity for histone binding. Conversely diaminopurine substitution decreases the affinity. In the sole example where replacement of guanosine with inosine decreases the persistence length as well as the affinity for histones, the substitution concomitantly removes an intrinsic curvature of the DNA molecule. We show that, to a first approximation, the binding energy of DNA to histones at 1M NaCl is directly proportional to the persistence length. The data also indicate that a high local flexibility of DNA can favour strong rotational positioning.

Site-specific Tn7 transposition into the human genome

The bacterial transposon, Tn7, inserts into a single site in the Escherichia coli chromosome termed attTn7 via the sequence-specific DNA binding of the target selector protein, TnsD. The target DNA sequence required for Tn7 transposition is located within the C-terminus of the glucosamine synthetase (glmS) gene, which is an essential, highly conserved gene found ubiquitously from bacteria to humans. Here, we show that Tn7 can transpose in vitro adjacent to two potential targets in the human genome: the gfpt-1 and gfpt-2 sequences, the human analogs of glmS. The frequency of transposition adjacent to the human gfpt-1 target is comparable with the E.coli glmS target; the human gfpt-2 target shows reduced transposition. The binding of TnsD to these sequences mirrors the transposition activity. In contrast to the human gfpt sequences, Tn7 does not transpose adjacent to the gfa-1 sequence, the glmS analog in Saccharomyces cerevisiae. We also report that a nucleosome core particle assembled on the human gfpt-1 sequence reduces Tn7 transposition by likely impairing the accessibility of target DNA to the Tns proteins. We discuss the implications of these findings for the potential use of Tn7 as a site-specific DNA delivery agent for gene therapy.

A determining influence for CpG dinucleotides on nucleosome positioning in vitro

DNA sequence information that directs the translational positioning of nucleosomes can be attenuated by cytosine methylation when a short run of CpG dinucleotides is located close to the dyad axis of the nucleosome. Here, we show that point mutations introduced to re-pattern methylation at the (CpG)3 element in the chicken betaA-globin promoter sequence themselves strongly influenced nucleosome formation in reconstituted chromatin. The disruptive effect of cytosine methylation on nucleosome formation was found to be determined by the sequence context of CpG dinucleotides, not just their location in the positioning sequence. Additional mutations indicated that methylation can also promote the occupation of certain nucleosome positions. DNase I analysis demonstrated that these genetic and epigenetic modifications altered the structural characteristics of the (CpG)3 element. Our findings support a proposal that the intrinsic structural properties of the DNA at the -1.5 site, as occupied by (CpG)3 in the nucleosome studied, can be decisive for nucleosome formation and stability, and that changes in anisotropic DNA bending or flexibility at this site explain why nucleosome positioning can be exquisitely sensitive to genetic and epigenetic modification of the DNA sequence.

Information content and complexity in the high-order organization of DNA

Nucleic acids are characterized by a vast structural variability. Secondary structural conformations include the main polymorphs A, B, and Z, cruciforms, intrinsic curvature, and multistranded motifs. DNA secondary motifs are stabilized and regulated by the primary base sequence, contextual effects, environmental factors, as well as by high-order DNA packaging modes. The high-order modes are, in turn, affected by secondary structures and by the environment. This review is concerned with the flow of structural information among the hierarchical structural levels of DNA molecules, the intricate interplay between the various factors that affect these levels, and the regulation and physiological significance of DNA high-order structures.

The structural basis of DNA flexibility

Although the average physico-chemical properties of a long DNA molecule may approximate to those of a thin isotropic homogeneous rod, DNA behaves more locally as an anisotropic heterogeneous rod. This bending anisotropy is sequence dependent and to a first approximation reflects both the geometry and stability of individual base steps. The biological manipulation and packaging of the molecule often depend crucially on local variations in both bending and torsional flexibility. However, whereas the probability of DNA untwisting can be strongly correlated with a high bending flexibility, DNA bending, especially when the molecule is tightly wrapped on a protein surface, may be energetically favoured by a less flexible sequence whose preferred configuration conforms more closely to that of the complementary protein surface. In the latter situation the lower bending flexibility may be more than compensated for on binding by a reduced required deformation energy relative to a fully isotropic DNA molecule.

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

Monovalent Cations Regulate DNA Sequence Recognition by 434 Repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.

Histone-DNA binding free energy cannot be measured in dilution-driven dissociation experiments

Despite decades of study on nucleosomes, there has been no experimental determination of the free energy of association between histones and DNA. Instead, only the relative free energy of association of the histone octamer for differing DNA sequences has been available. Recently, a method was developed based on quantitative analysis of nucleosome dissociation in dilution experiments that provides a simple practical measure of nucleosome stability. Solution conditions were found in which nucleosome dissociation driven by dilution fit well to a simple model involving a noncooperative nucleosome assembly/disassembly equilibrium, suggesting that this approach might allow absolute equilibrium affinity of the histone octamer for DNA to be measured. Here, we show that the nucleosome assembly/disassembly process is not strictly reversible in these solution conditions, implying that equilibrium affinities cannot be obtained from these measurements. Increases in [NaCl] or temperature, commonly employed to suppress kinetic bottlenecks in nucleosome assembly, lead to cooperative behavior that cannot be interpreted with the simple assembly/disassembly equilibrium model. We conclude that the dilution experiments provide useful measures of kinetic but not equilibrium stability. Kinetic stability is of practical importance: it may govern nucleosome function in vivo, and it may (but need not) parallel absolute thermodynamic stability.

DNA Deformability at the Base Pair Level

A complete set of harmonic force constants describing the DNA deformation energetics at the base pair level was obtained using unrestrained atomic-resolution molecular dynamics simulations of selected duplex oligonucleotides and subsequent analysis of structural fluctuations from the simulated trajectories. The deformation was described by the six base pair conformational parameters (buckle, propeller, opening, shear, stretch, stagger). The results for 13 AT pairs and 11 GC pairs in different sequence contexts suggest that buckle and propeller are very flexible (more than roll in TA dinucleotide steps), while stretch is exceptionally stiff. Only stretch and opening stiffness were found to depend unambiguously on the base pair identity (AT vs GC). The relationship of the results to a simple plates-and-springs model of base-base interactions is discussed.

Analysis of Chemical and Enzymatic Cleavage Frequencies in Supercoiled DNA

Chemical and enzymatic probing methods are powerful techniques for examining details of sequence-dependent structure in DNA and RNA. Reagents that cleave nucleic acid molecules in a structure-specific, but relatively sequence-non-specific manner, such as hydroxyl radical or DNase I, have been used widely to probe helical geometry in nucleic acid structures, nucleic acid-drug complexes, and in nucleoprotein assemblies. Application of cleavage-based techniques to structures present in superhelical DNA has been hindered by the fact that the cleavage pattern attributable to supercoiling-dependent structures is heavily mixed with non-specific cleavage signals that are inevitable products of multiple cleavage events. We present a rigorous mathematical procedure for extracting the cleavage pattern specific to supercoiled DNA and use this method to investigate the hydroxyl radical cleavage pattern in a cruciform DNA structure formed by a 60 bp inverted repeat sequence embedded in a negatively supercoiled plasmid. Our results support the presence of a stem-loop structure in the expected location and suggest that the helical geometry of the cruciform stem differs from that of the normal duplex form.

An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data

Curved DNA fragments are often found near functionally important sites such as promoters and origins of replication, and hence sequence-dependent DNA curvature prediction is of great utility in genomics and bioinformatics. In light of this, an assessment of three different dinucleotide step parameters (based on gel retardation as well as crystal structure data) is carried out. These parameters (BMHT, LB and CS) are evaluated quantitatively for their ability to predict correctly the experimental results of a large set of nucleic acid sequences containing A-tracts as well as GC-rich motifs. This set contained around 40 synthetic as well as natural sequences whose solution properties have been well characterized experimentally. All three models could account reasonably well for curvature in the various DNA sequences. The CS model, where dinucleotide parameters are calculated from crystal structure data, consistently shows slightly better correlation with experimental data. Our simple analysis also indicates that presently available trinucleotide parameters fail to predict curvature in some of the well-characterized sequences. The study shows that the dinucleotide parameters with some further refinement can be used to predict sequence-dependent curvature correctly in genomic sequences.