Helicoidal transfer matrix model for inhomogeneous DNA melting

Synchronized Oscillations in Double-Helix B-DNA Molecules with Mirror-Symmetric Codons

A fully analytical treatment of the base-pair and codon dynamics in double-stranded DNA molecules is introduced, by means of a realistic treatment that considers different mass values for G, A, T, and C nucleotides and takes into account the intrinsic three-dimensional, helicoidal geometry of DNA in terms of a Hamitonian in cylindrical coordinates. Within the framework of the Peyrard–Dauxois–Bishop model, we consider the coupling between stretching and stacking radial oscillations as well as the twisting motion of each base pair around the helix axis. By comparing the linearized dynamical equations for the angular and radial variables corresponding to the bp local scale with those of the longer triplet codon scale, we report an underlying hierarchical symmetry. The existence of synchronized collective oscillations of the base-pairs and their related codon triplet units are disclosed from the study of their coupled dynamical equations. The possible biological role of these correlated, long-range oscillation effects in double standed DNA molecules containing mirror-symmetric codons of the form XXX, XX’X, X’XX’, YXY, and XYX is discussed in terms of the dynamical equations solutions and their related dispersion relations.

Base-Pairs' Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules

By introducing a suitable renormalization process, the charge carrier and phonon dynamics of a double-stranded helical DNA molecule are expressed in terms of an effective Hamiltonian describing a linear chain, where the renormalized transfer integrals explicitly depend on the relative orientations of the Watson-Crick base pairs, and the renormalized on-site energies are related to the electronic parameters of consecutive base pairs along the helix axis, as well as to the low-frequency phonons' dispersion relation. The existence of synchronized collective oscillations enhancing the π-π orbital overlapping among different base pairs is disclosed from the study of the obtained analytical dynamical equations. The role of these phonon-correlated, long-range oscillation effects on the charge transfer properties of double-stranded DNA homopolymers is discussed in terms of the resulting band structure.

Thermal and mechanical denaturation properties of a DNA model with three sites per nucleotide

In this paper, we show that the coarse grain model for DNA, which has been proposed recently by Knotts et al. [J. Chem. Phys. 126, 084901 (2007)], can be adapted to describe the thermal and mechanical denaturation of long DNA sequences by adjusting slightly the base pairing contribution. The adjusted model leads to (i) critical temperatures for long homogeneous sequences that are in good agreement with both experimental ones and those obtained from statistical models, (ii) a realistic step-like denaturation behaviour for long inhomogeneous sequences, and (iii) critical forces at ambient temperature of the order of 10 pN, close to measured values. The adjusted model furthermore supports the conclusion that the thermal denaturation of long homogeneous sequences corresponds to a first-order phase transition and yields a critical exponent for the critical force equal to σ = 0.70. This model is both geometrically and energetically realistic, in the sense that the helical structure and the grooves, where most proteins bind, are satisfactorily reproduced, while the energy and the force required to break a base pair lie in the expected range. It therefore represents a promising tool for studying the dynamics of DNA-protein specific interactions at an unprecedented detail level.

Compatible solute influence on nucleic acids: many questions but few answers

Compatible solutes are small organic osmolytes including but not limited to sugars, polyols, amino acids, and their derivatives. They are compatible with cell metabolism even at molar concentrations. A variety of organisms synthesize or take up compatible solutes for adaptation to extreme environments. In addition to their protective action on whole cells, compatible solutes display significant effects on biomolecules in vitro. These include stabilization of native protein and nucleic acid structures. They are used as additives in polymerase chain reactions to increase product yield and specificity, but also in other nucleic acid and protein applications.

Interactions of compatible solutes with nucleic acids and protein-nucleic acid complexes are much less understood than the corresponding interactions of compatible solutes with proteins. Although we may begin to understand solute/nucleic acid interactions there are only few answers to the many questions we have. I summarize here the current state of knowledge and discuss possible molecular mechanisms and thermodynamics.

Stability domains of actin genes and genomic evolution

In eukaryotic genes, the protein coding sequence is split into several fragments, the exons, separated by noncoding DNA stretches, the introns. Prokaryotes do not have introns in their genomes. We report calculations of the stability domains of actin genes for various organisms in the animal, plant, and fungi kingdoms. Actin genes have been chosen because they have been highly conserved during evolution. In these genes, all introns were removed so as to mimic ancient genes at the time of the early eukaryotic development, i.e., before intron insertion. Common stability boundaries are found in evolutionarily distant organisms, which implies that these boundaries date from the early origin of eukaryotes. In general, the boundaries correspond with intron positions in the actins of vertebrates and other animals, but not much for plants and fungi. The sharpest boundary is found in a locus where fungi, algae, and animals have introns in positions separated by one nucleotide only, which identifies a hot spot for insertion. These results suggest that some introns may have been incorporated into the genomes through a thermodynamically driven mechanism, in agreement with previous observations on human genes. They also suggest a different mechanism for intron insertion in plants and animals.

The human genomic melting map

In a living cell, the antiparallel double-stranded helix of DNA is a dynamically changing structure. The structure relates to interactions between and within the DNA strands, and the array of other macromolecules that constitutes functional chromatin. It is only through its changing conformations that DNA can organize and structure a large number of cellular functions. In particular, DNA must locally uncoil, or melt, and become single-stranded for DNA replication, repair, recombination, and transcription to occur. It has previously been shown that this melting occurs cooperatively, whereby several base pairs act in concert to generate melting bubbles, and in this way constitute a domain that behaves as a unit with respect to local DNA single-strandedness. We have applied a melting map calculation to the complete human genome, which provides information about the propensities of forming local bubbles determined from the whole sequence, and present a first report on its basic features, the extent of cooperativity, and correlations to various physical and biological features of the human genome. Globally, the melting map covaries very strongly with GC content. Most importantly, however, cooperativity of DNA denaturation causes this correlation to be weaker at resolutions fewer than 500 bps. This is also the resolution level at which most structural and biological processes occur, signifying the importance of the informational content inherent in the genomic melting map. The human DNA melting map may be further explored at http://meltmap.uio.no.

Use of DNA melting simulation software for in silico diagnostic assay design: targeting regions with complex melting curves and confirmation by real-time PCR using intercalating dyes

BACKGROUND

DNA melting curve analysis using double-stranded DNA-specific dyes such as SYTO9 produce complex and reproducible melting profiles, resulting in the detection of multiple melting peaks from a single amplicon and allowing the discrimination of different species. We compare the melting curves of several Naegleria and Cryptosporidium amplicons generated in vitro with in silico DNA melting simulations using the programs POLAND and MELTSIM., then test the utility of these programs for assay design using a genetic marker for toxin production in cyanobacteria.

RESULTS

The SYTO9 melting curve profiles of three species of Naegleria and two species of Cryptosporidium were similar to POLAND and MELTSIM melting simulations, excepting some differences in the relative peak heights and the absolute melting temperatures of these peaks. MELTSIM and POLAND were used to screen sequences from a putative toxin gene in two different species of cyanobacteria and identify regions exhibiting diagnostic melting profiles. For one of these diagnostic regions the POLAND and MELTSIM melting simulations were observed to be different, with POLAND more accurately predicting the melting curve generated in vitro. Upon further investigation of this region with MELTSIM, inconsistencies between the melting simulation for forward and reverse complement sequences were observed. The assay was used to accurately type twenty seven cyanobacterial DNA extracts in vitro.

CONCLUSION

Whilst neither POLAND nor MELTSIM simulation programs were capable of exactly predicting DNA dissociation in the presence of an intercalating dye, the programs were successfully used as tools to identify regions where melting curve differences could be exploited for diagnostic melting curve assay design. Refinements in the simulation parameters would be required to account for the effect of the intercalating dye and salt concentrations used in real-time PCR. The agreement between the melting curve simulations for different species of Naegleria and Cryptosporidium and the complex melting profiles generated in vitro using SYTO9 verified that the complex melting profile of PCR amplicons was solely the result of DNA dissociation. Other data outputs from these simulations were also used to identify the melting domains that contributed to the observed melting peaks for each of the different PCR amplicons.