A Statistical Thermodynamic Model for Investigating the Stability of DNA Sequences from Oligonucleotides to Genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of −0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.     

A Statistical Thermodynamic Model for Investigating the Stability of DNA Sequences from Oligonucleotides to Genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of −0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.     

Parallel-Stranded DNA with Natural Base Sequences

Noncanonical parallel-stranded DNA double helices (ps-DNA) of natural nucleotide sequences are usually less stable than the canonical antiparallel-stranded DNA structures, which ensures reliable cell functioning. However, recent data indicate a possible role of ps-DNA in DNA loops or in regions of trinucleotide repeats connected with neurodegenerative diseases. The review surveys recent studies on the effect of nucleotide sequence on preference of one or other type of DNA duplex. (1) Ps-DNA of mixed AT/GC composition was found to have conformational and thermodynamic properties drastically different from those of a Watson–Crick double helix. Its stability depends strongly on the specific sequence in a manner peculiar to the ps double helix, because of the energy disadvantage of the AT/GC contacts. The AT/GC boundary facilitated flipping of A and T out of the ps double helix. Proton acceptor groups of bases are exposed into both grooves of the ps-DNA and are accessible to solvent and ligands, including proteins. (2) DNA regions containing natural minor bases isoguanine and isomethylisocytosine were shown to form ps-DNA with transAT-, trans isoGC, and transiso^5meCG pairs exceeding in stability a related canonical duplex. (3) Nucleotide sequence dG(GT)₄G from yeast telomeres and microsatellites was demonstrated to form novel ps-DNA with GG and TT base pairing. Unlike d(GT) _n - and d(G _n T _m ) sequences able to form quadruplexes, the dG(GT)₄G sequence formed no alternative double- or multistranded structures in a wide range of experimental conditions, thus suggesting that the nucleotide context governs the observed structural polymorphism of the d(GT) _n sequence. The possible biological role of ps-DNA and the prospects of its study are discussed.     

The nucleotide sequence of the mitochondrial DNA genome of an abundant petite mutant of Saccharomyces cerevisiae carrying the ori1 replication origin

We have determined the 903 bp nucleotide sequence of the mitochondrial DNA genome of a Saccharomyces cerevisiae petite mutant BB5. This petite, containing the 265 nucleotide ori1 region, is representative of a class of petites arising at exceptionally high frequency within the population of spontaneous petites derived from a particular mit- strain Mb12. The DNA sequences of both the ori1 region and the flanking intergenic regions have been compared to those of the corresponding regions of mtDNA in a previously reported petite strain, a1/1R/1 of Bernardi's laboratory, that has a similar (880 bp) repeat unit. The BB5 petite genome carries a canonical ori1 sequence that is identical in both petite mtDNAs, but the flanking intergenic sequences show significant differences between the two petite strains. The divergence is considered to arise from differences in the sequences flanking ori1 in the respective parent strains.     

Parallel-stranded DNA with natural base sequences

Noncanonical parallel-stranded DNA double helices (ps-DNA) comprising natural nucleotide sequences are usually second in stability to antiparallel-stranded (aps) canonical DNA structures, which ensures reliable cell functioning. However, recent data indicate a possible role of ps-DNA in DNA loops or in trinucleotide repeats connected with neurodegenerative diseases. The review surveys recent studies on the effect of nucleotide sequence on preference of one or other type of DNA duplex. (1) Ps-DNA with mixed AT/GC composition was found to have conformational and thermodynamic properties drastically different from those of Watson-Crick double helix. Its stability depends strongly on the specific sequence in a manner peculiar to the ps double helix, because of the energy disadvantage of the AT/GC contacts. The AT/GC boundary facilitated flipping of A and T out of the ps double helix. Proton acceptor groups of bases are exposed into the both grooves of the ps-DNA and are accessible to solvent and ligands, including proteins. (2) DNA regions containing natural minor bases isoguanine and isomethylcytosine were shown to form ps-DNA with transAT-, trans isoGC, and trans iso5meCG pairs exceeding in stability a related aps duplex. (3) Nucleotide sequence dG(GT)4G from yeast telomeres and microsatellites was demonstrated to form novel ps-DNA with GG and TT base pairing. Unlike d(GT)n and d(GnTm) sequences able to form quadruplexes, the dG(GT)4G sequence formed no alternative double- or multistranded structures in a wide range of experimental conditions, thus suggesting that the nucleotide context governs the observed structural polymorphism of the d(GT)n sequence. The possible biological role of ps-DNA and the prospects of its study are discussed.     

Simian virus 40 and polyomavirus large tumor antigens have different requirements for high-affinity sequence-specific DNA binding.

By using a DNA fragment immunoassay, the binding of simian virus 40 (SV40) and polyomavirus (Py) large tumor (T) antigens to regulatory regions at both viral origins of replication was examined. Although both Py T antigen and SV40 T antigen bind to multiple discrete regions on their proper origins and the reciprocal origin, several striking differences were observed. Py T antigen bound efficiently to three regions on Py DNA centered around an MboII site at nucleotide 45 (region A), a BglI site at nucleotide 92 (region B), and another MboII site at nucleotide 132 (region C). Region A is adjacent to the viral replication origin, and region C coincides with the major early mRNA cap site. Weak binding by Py T antigen to the origin palindrome centered at nucleotide 3 also was observed. SV40 T antigen binds strongly to Py regions A and B but only weakly to region C. This weak binding on region C was surprising because this region contains four tandem repeats of GPuGGC, the canonical pentanucleotide sequence thought to be involved in specific binding by T antigens. On SV40 DNA, SV40 T antigen displayed its characteristic hierarchy of affinities, binding most efficiently to site 1 and less efficiently to site 2. Binding to site 3 was undetectable under these conditions. In contrast, Py T antigen, despite an overall relative reduction of affinity for SV40 DNA, binds equally to fragments containing each of the three SV40 binding sites. Py T antigen, but not SV40 T antigen, also bound specifically to a region of human Alu DNA which bears a remarkable homology to SV40 site 1. However, both tumor antigens fail to precipitate DNA from the same region which has two direct repeats of GAGGC. These results indicate that despite similarities in protein structure and DNA sequence, requirements of the two T antigens for pentanucleotide configuration and neighboring sequence environment are different.     

A computer aided thermodynamic approach for predicting the formation of Z-DNA in naturally occurring sequences.

The ease with which a particular DNA segment adopts the left-handed Z-conformation depends largely on the sequence and on the degree of negative supercoiling to which it is subjected. We describe a computer program (Z-hunt) that is designed to search long sequences of naturally occurring DNA and retrieve those nucleotide combinations of up to 24 bp in length which show a strong propensity for Z-DNA formation. Incorporated into Z-hunt is a statistical mechanical model based on empirically determined energetic parameters for the B to Z transition accumulated to date. The Z-forming potential of a sequence is assessed by ranking its behavior as a function of negative superhelicity relative to the behavior of similar sized randomly generated nucleotide sequences assembled from over 80,000 combinations. The program makes it possible to compare directly the Z-forming potential of sequences with different base compositions and different sequence lengths. Using Z-hunt, we have analyzed the DNA sequences of the bacteriophage phi X174, plasmid pBR322, the animal virus SV40 and the replicative form of the eukaryotic adenovirus-2. The results are compared with those previously obtained by others from experiments designed to locate Z-DNA forming regions in these sequences using probes which show specificity for the left-handed DNA conformation.     

Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

DNA rearrangement associated with the integration of T-DNA in tobacco: an example for multiple duplications of DNA around the integration target

Transferred DNA (T-DNA) of the tumor-inducing (Ti) plasmid is transferred from Agrobacterium tumefaciens to plant cells and is stably integrated into the plant nuclear genome. By the inverse polymerase chain reaction DNA fragments were amplified that contained the T-DNA/plant DNA junctions from the total DNA of a transgenic tobacco plant that had a single copy of the T-DNA in a repetitive region of its genome. A DNA fragment containing the target site was amplified from the total DNA of non-transformed tobacco by the polymerase chain reaction using high-stringency conditions. Comparison of the nucleotide sequence of the target site with those of the T-DNA/plant DNA junctions revealed that various duplications of short stretches of nucleotide sequences around the target and in the incoming T-DNA had accompanied the integration of the T-DNA. A deletion of 16 bp at the target site was also found and the target site was similar, in terms of nucleotide sequence, to regions around the breakpoints of the T-DNA. This finding provides a clear example of the occurrence of complex rearrangements during the integration of T-DNA.     

Sequence of the Saccharomyces cerevisiae PHR1 gene and homology of the PHR1 photolyase to E. coli photolyase.

The nucleotide sequence of a 2301 base pair region of Saccharomyces cerevisiae DNA containing the PHR1 gene is reported. Within this region a single open reading frame of 1695 base pairs was found; using the insertional inactivation technique it was shown that part or all of this open reading frame specifies the PHR1-encoded photolyase. The amino acid sequence of the 565 amino acid long polypeptide predicted from the PHR1 nucleotide sequence was compared to the amino acid sequence of E. coli photolyase. Overall the sequence homology was 36.5%; however, two short regions near the amino terminus as well as the carboxy-terminal 150 amino acids display significantly greater sequence homology. The presence of these strongly conserved regions suggests that the yeast and E. coli photolyase possess common structural and functional domains involved in substrate and/or chromophore binding.     

DNA analysis servers: plot.it,bend.it,model.it and IS

The WWW servers at http://www.icgeb.trieste.it/dna/ are dedicated to the analysis of user-submitted DNA sequences; plot.it creates parametric plots of 45 physicochemical, as well as statistical, parameters; bend.it calculates DNA curvature according to various methods. Both programs provide 1D as well as 2D plots that allow localisation of peculiar segments within the query. The server model.it creates 3D models of canonical or bent DNA starting from sequence data and presents the results in the form of a standard PDB file, directly viewable on the user's PC using any molecule manipulation program. The recently established introns server allows statistical evaluation of introns in various taxonomic groups and the comparison of taxonomic groups in terms of length, base composition, intron type etc. The options include the analysis of splice sites and a probability test for exon-shuffling.     

Nucleotide sequence of the R.meliloti nitrogenase reductase (nifH) gene.

The nucleotide sequence of the structural gene (nifH) of nitrogenase reductase (Fe protein) from R.meliloti 41 with its flanking ends is reported. The amino acid sequence of nitrogenase reductase was deduced from the DNA sequence. The predicted R.meliloti nitrogenase reductase protein consists of 297 amino acid residues, has a molecular weight of 32,740 daltons and contains 5 cysteine residues. The codon usage in the nifH gene is presented. In the 5' flanking region, sequences resembling to consensus sequences of bacterial control regions were found. Comparison of the R.meliloti nifH nucleotide and amino acid sequences with those from different nitrogen-fixing organisms showed that the amino acid sequences are more conserved than the nucleotide sequences. This structural conservation of nitrogenase reductase may be related to its function and may explain the conservation of the nifH gene during evolution.     

A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.     

Secondary structure prediction and structure-specific sequence analysis of single-stranded DNA

DNA sequence analysis by oligonucleotide binding is often affected by interference with the secondary structure of the target DNA. Here we describe an approach that improves DNA secondary structure prediction by combining enzymatic probing of DNA by structure-specific 5′-nucleases with an energy minimization algorithm that utilizes the 5′-nuclease cleavage sites as constraints. The method can identify structural differences between two DNA molecules caused by minor sequence variations such as a single nucleotide mutation. It also demonstrates the existence of long-range interactions between DNA regions separated by >300 nt and the formation of multiple alternative structures by a 244 nt DNA molecule. The differences in the secondary structure of DNA molecules revealed by 5′-nuclease probing were used to design structure-specific probes for mutation discrimination that target the regions of structural, rather than sequence, differences. We also demonstrate the performance of structure-specific ‘bridge’ probes complementary to non-contiguous regions of the target molecule. The structure-specific probes do not require the high stringency binding conditions necessary for methods based on mismatch formation and permit mutation detection at temperatures from 4 to 37°C. Structure-specific sequence analysis is applied for mutation detection in the Mycobacterium tuberculosis katG gene and for genotyping of the hepatitis C virus.