Bubbles and denaturation in DNA

Force-induced unzipping of DNA in the presence of solvent molecules

The melting of double-stranded DNA (dsDNA) in the presence of solvent molecules is a fundamental process with significant implications for understanding the thermal and mechanical behavior of DNA and its interactions with the surrounding environment. The solvents play an essential role in the structural transformation of DNA subjected to a pulling force. In this study, we simulate the thermal and force induced denaturation of dsDNA and elucidate the solvent dependent melting behavior, identifying key factors that influence the stability of DNA melting in presence of solvent molecules. Using a statistical model, we first find the melting profile of short heterogeneous DNA molecules in the presence of solvent molecules in Force ensemble. We also investigate the effect of solvent's strengths on the melting profile of DNA. In the force ensemble, we consider two homogeneous DNA chains and apply the force on different locations along the chain in the presence of solvent molecules. Different pathways manifest the melting of the molecule in both ensembles, and we found several interesting features of melting DNA in a constant force ensemble, such as lower critical force when the chain is pulled from the base pair close to a solvent molecule. The results provide new insights into the force-induced unzipping of DNA and could be used to develop new methods for controlling the unzipping process. By providing a better understanding of melting and unzipping of dsDNA in the presence of solvent molecules, this study provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA nanostructures.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Non-canonical DNA structures in the human ribosomal DNA

Non-canonical structures (NCS) refer to the various forms of DNA that differ from the B-conformation described by Watson and Crick. It has been found that these structures are usual components of the genome, actively participating in its essential functions. The present review is focused on the nine kinds of NCS appearing or likely to appear in human ribosomal DNA (rDNA): supercoiling structures, R-loops, G-quadruplexes, i-motifs, DNA triplexes, cruciform structures, DNA bubbles, and A and Z DNA conformations. We discuss the conditions of their generation, including their sequence specificity, distribution within the locus, dynamics, and beneficial and detrimental role in the cell.

Equilibrium melting probabilities of a DNA molecule with a defect: An exact solution of the Poland-Scheraga model

In this study we derive analytically the equilibrium melting probabilities for basepairs of a DNA molecule with a defect site. We assume that the defect is characterized by a change in the Watson-Crick basepair energy of the defect basepair, and in the associated two stacking energies for the defect, as compared to the remaining parts of the DNA. The defect site could, for instance, occur due to DNA basepair mismatching, cross-linking, or by the chemical modifications when attaching fluorescent labels, such as fluorescent-quencher pairs, to DNA. Our exact solution of the Poland-Scheraga model for DNA melting provides the probability that the labeled basepair, and its neighbors, are open at different temperatures. Our work is of direct importance, for instance, for studies where fluorophore-quencher pairs are used for studying single basepair fluctuations of designed DNA molecules.

Melting of dsDNA attached with AuNPs

DNA-linked gold nanoparticles (DNA-AuNPs) are combined nanomaterials that contain the optical and electronic properties of AuNPs with the unique functions of DNA. These hybrid systems are used in various nanobiotechnology, medical, and pharmaceutical sciences (Löwe et al. in FEBS J 287(23):5039, 2020; Speer et al. in Annu Rev Biophys 51:267, 2022). In recent years, there has been an increasing interest in studying the behavior of DNA-AuNPs in the presence of molecular solvents. In the present work, we study the thermal melting of DNA-linked gold nanoparticles (DNA-AuNP). In the first part of the study, we find the melting profile of short heterogeneous DNA-linked AuNP in the presence of solvent in the solution. We also study the effect of the location of the gold nanoparticle attached to the DNA molecule. In this case, we move the location of the AuNP from one end to the other. We found that while the melting temperature is susceptible to the location of the AuNP when it is near the ends, there is a region in the middle section of the chain where the melting temperature remains constant.

Pulling short DNA with mismatch base pairs

Due to misincorporation during gene replication, the accuracy of the gene expression is often compromised. This results in a mismatch or defective pair in the DNA molecule (James et al. 2016). Here, we present our study of the stability of DNA with defects in the thermal and force ensembles. We consider DNA with a different number of defects from 2to16 and study how the denaturation process differs in both ensembles. Using a statistical model, we calculate the melting point of the DNA chain in both the ensemble. Our findings display different manifestations of DNA denaturation in thermal and force ensembles. While the DNA with defects denatures at a lower temperature than the intact DNA, the point from which the DNA is pulled is important in force ensemble.

A predictive model for the thermomechanical melting transition of double stranded DNA

By extending the classical Peyrard-Bishop model, we are able to obtain a fully analytical description for the mechanical response of DNA under stretching at variable values of temperature, number of base pairs and intrachains and interchains bonds stiffness. In order to compare elasticity and temperature effects, we first analyze the system in the zero temperature mechanical limit, important to describe several experimental effects including possible hysteresis. We then analyze temperature effects in the framework of equilibrium Statistical Mechanics. In particular, we obtain an analytical expression for the temperature-dependent melting force and unzipping assigned displacement in the thermodynamical limit, also depending on the relative stability of intra vs. inter molecular bonds. Such results coincide with the purely mechanical model in the limit of zero temperature and with the denaturation temperature that we obtain with the classical transfer integral method. Based on our analytical results, we obtain explicitly phase diagrams and cooperativity parameters, where also discreteness effect can be accounted for. The obtained results are successfully applied in reproducing the thermomechanical experimental melting of DNA and the response of DNA hairpins. Due to the generality of the model, exemplified in the proposed analysis of both overstretching and unzipping experiments, we argue that the proposed approach can be extended to other thermomechanically induced molecular melting phenomena. STATEMENT OF SIGNIFICANCE: We obtain a fully analytical description of the complex wiggly energy landscape of two stranded macromolecules under unzipping loading. Based on Equilibrium Statistical Mechanics, we describe the combined thermomechanical effects and the melting transition of double stranded molecules such as nucleic acids. This is proved by quantitatively predicting the experimental behavior of both melting of DNA and DNA hairpins opening. While analytical results have been previously attained under special conditions on the relative stiffness of the covalent vs. non-covalent bonds of the base pairs, our model is completely general in this respect, thus representing a tool in the perspective of the design at the molecular scale. We show that the obtained model can be fully inscribed in the theory of phase transitions giving a new interpretation of the thermomechanical behavior of double stranded molecules.

Nonlinear physics opens a new paradigm for accurate transcription start site prediction

There is evidence that DNA breathing (spontaneous opening of the DNA strands) plays a relevant role in the interactions of DNA with other molecules, and in particular in the transcription process. Therefore, having physical models that can predict these openings is of interest. However, this source of information has not been used before either in transcription start sites (TSSs) or promoter prediction. In this article, one such model is used as an additional information source that, when used by a machine learning (ML) model, improves the results of current methods for the prediction of TSSs. In addition, we provide evidence on the validity of the physical model, as it is able by itself to predict TSSs with high accuracy. This opens an exciting avenue of research at the intersection of statistical mechanics and ML, where ML models in bioinformatics can be improved using physical models of DNA as feature extractors.

A review on nonlinear DNA physics

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo, possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field.

Melting of DNA in confined geometries

The stability of DNA molecules during viral or biotechnological encapsulation is a topic of active current research. We studied the thermal stability of double-stranded DNA molecules of different lengths in a confined space. Using a statistical model, we evaluate the melting profile of DNA molecules in two geometries: conical and cylindrical. Our results show that not only the confinement, but also the geometry of the confined space plays a prominent role in the stability and opening of the DNA duplex. We find that for more confined spaces, cylindrical confinement stabilizes the DNA, but for less confined spaces conical geometry stabilizes the DNA overall. We also analyse the interaction between DNA sequence and stability, and the evenness with which strand separation occurs. Cylindrical and conical geometries enable a better controlled tuning of the stability of DNA encapsulation and the efficiency of its eventual release, compared to spherical or quasi-spherical geometries.

Inequality in the Frequency of the Open States Occurrence Depends on Single 2H/1H Replacement in DNA

In the present study, the effect of 2H/1H isotopic exchange in hydrogen bonds between nitrogenous base pairs on occurrence and open states zones dynamics is investigated. These processes are studied using mathematical modeling, taking into account the number of open states between base pairs. The calculations of the probability of occurrence of open states in different parts of the gene were done depending on the localization of the deuterium atom. The mathematical modeling study demonstrated significant inequality (dependent on single 2H/1H replacement in DNA) among three parts of the gene similar in length of the frequency of occurrence of the open states. In this paper, the new convenient approach of the analysis of the abnormal frequency of open states in different parts of the gene encoding interferon alpha 17 was presented, which took into account both rising and decreasing of them that allowed to make a prediction of the functional instability of the specific DNA regions. One advantage of the new algorithm is diminishing the number of both false positive and false negative results in data filtered by this approach compared to the pure fractile methods, such as deciles or quartiles.

Improving the mesoscopic modeling of DNA denaturation dynamics

Although previously developed mesoscopic DNA models have successfully reproduced thermodynamic denaturation data, recent studies show that these overestimate the rate of denaturation by orders of magnitude. Using adapted Peyrard-Bishop-Dauxois (PBD) models, we have calculated the denaturation rates of several DNA hairpins and made comparison with experimental data. We show that the addition of a barrier at the onsite potential of the PBD model gives a more accurate description of the unzipping dynamics of short DNA sequences. The new models provide a refined theoretical insight on the dynamical mechanisms of unzipping which can have implications for the understanding of transcription and replication processes. Still, this class of adapted PBD models seems to have a fundamental limitation which implies that it is not possible to get agreement with available experimental results on the dynamics of long DNA sequences and at the same time maintain the good agreement regarding its thermodynamics. The reason for this is that the denaturation rate of long DNA chains is not dramatically lowered by the additional barrier-as the base-pairs that open are more likely to remain open, facilitating the opening of the full DNA molecule. Some care has to be taken, since experimental techniques suitable to the study of denaturation rates of long sequences seem not to agree with other experimental data on short DNA sequences. Further research, both theoretical and experimental, is therefore needed to resolve these inconsistencies-which will be a starting point for new minimalistic models that are able to describe both thermodynamics and dynamics at a predictive level.

Kinetic theory for DNA melting with vibrational entropy

By treating DNA as a vibrating nonlinear lattice, an activated kinetic theory for DNA melting is developed to capture the breakage of the hydrogen bonds and subsequent softening of torsional and bending vibration modes. With a coarse-grained lattice model, we identify a key bending mode with GHz frequency that replaces the hydrogen vibration modes as the dominant out-of-phase phonon vibration at the transition state. By associating its bending modulus to a universal in-phase bending vibration modulus at equilibrium, we can hence estimate the entropic change in the out-of-phase vibration from near-equilibrium all-atom simulations. This and estimates of torsional and bending entropy changes lead to the first predictive and sequence-dependent theory with good quantitative agreement with experimental data for the activation energy of melting of short DNA molecules without intermediate hairpin structures.

Direct observation of breathing dynamics at the mismatch induced DNA bubble with nanometre accuracy: a smFRET study

The detailed conformational dynamics of the melted region in double-stranded DNA has been studied using a combination of ensemble and single-molecule FRET techniques. We monitored the millisecond time scale fluctuation kinetics of the two strands at the bubble region that varies with the size of the bubble. As the individual strands at the melting bubble behave as single-stranded DNA, and hence fluctuate dynamically to attain energetically favored configurations, the rates of these fluctuations increase with increase in the bubble size. In different short DNAs under investigation, the two strands never cross each other to form a knot, irrespective of the number of base pair mismatches present. Rather, they prefer to stay apart from each other, as the size of the bubble increases and follow exactly an opposite trend for bubbles of smaller size. The range within which the bubble strands fluctuate are monitored with great accuracy in the nanometre resolution from the single-molecule FRET measurements. The shape of the bubble that plays a crucial role in determining the activity of the DNA was speculated. These results shall be useful in quantifying the chemical processes within DNA as well as to develop a deeper understanding of the activity of the DNA due to induced mismatches.

Triggered isothermal PCR by denaturation bubble-mediated strand exchange amplification

Here, we introduced the concept of strand exchange amplification (SEA) mediated by denaturation bubbles.

Pulling short DNA molecules having defects on different locations

We present a study on the role of defects on the stability of short DNA molecules. We consider short DNA molecules (16 base pairs) and investigate the thermal as well as mechanical denaturation of these molecules in the presence of defects that occur anywhere in the molecule. For the investigation, we consider four different kinds of chains. Not only are the ratios of AT to GC different in these molecules but also the distributions of AT and GC along the molecule are different. With suitable modifications in the statistical model to show the defect in a pair, we investigate the denaturation of short DNA molecules in thermal as well as constant force ensembles. In the force ensemble, we pulled the DNA molecule from each end (keeping other end free) and observed some interesting features of opening of the molecule in the presence of defects in the molecule. We calculate the probability of opening of the DNA molecule in the constant force ensemble to explain the opening of base pairs and hence the denaturation of molecules in the presence of defects.

Twisting and bending stress in DNA minicircles

The interplay between bending of the molecule axis and appearance of disruptions in circular DNA molecules, with ∼100 base pairs, is addressed. Three minicircles with different radii and almost equal contents of AT and GC pairs are investigated. The DNA sequences are modeled by a mesoscopic Hamiltonian which describes the essential interactions in the helix at the level of the base pair and incorporates twisting and bending degrees of freedom. Helix unwinding and bubble formation patterns are consistently computed by a path integral method that sums over a large number of molecule configurations compatible with the model potential. The path ensembles are determined, as a function of temperature, by minimizing the free energy of the system. Fluctuational openings appear along the helix to release the stress due to the bending of the molecule backbone. In agreement with the experimental findings, base pair disruptions are found with larger probability in the smallest minicircle of 66 bps whose bending angle is ∼6°. For this minicircle, a sizeable untwisting is obtained with the helical repeat showing a step-like increase at T = 315 K. The method can be generalized to determine the bubble probability profiles of open ends linear sequences.

Modeling the interaction of DNA with alternating fields

We study the influence of a terahertz field on thermal properties of DNA molecules. A Peyrard-Bishop-Dauxois model with the inclusion of a solvent interaction term is considered. The terahertz field is included as a sinusoidal driven force in the equation of motion. We show how under certain field and system parameters, the melting transition and bubble formation are modified.

Effect of temperature on DNA double helix: An insight from molecular dynamics simulation

The three-dimensional structure of DNA contains various sequence-dependent structural information, which control many cellular processes in life, such as replication, transcription, DNA repair, etc. For the above functions, DNA double helices need to unwind or melt locally, which is different from terminal melting, as often seen in molecular dynamics (MD) simulations or even in many DNA crystal structures. We have carried out detailed MD simulations of DNA double helices of regular oligonucleotide fragments as well as in polymeric constructs with water and charge-neutralizing counter-ions at several different temperatures. We wanted to eliminate the end-effect or terminal melting propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The polymeric construct is expected to allow us to see local melting at elevated temperatures. Comparative structural analysis of oligonucleotides and its corresponding virtual polymer at various temperatures ranging from 300 K to 400 K is discussed. The general behaviour, such as volume expansion coefficients of both the simulations show high similarity, indicating polymeric construct, does not give many artificial constraints. Local melting of a polymer, even at elevated temperature, may need a high nucleation energy that was not available in the short (7 ns) simulations. We expected to observe such nucleation followed by cooperative melting of the polymers in longer MD runs. Such simulations of different polymeric sequences would facilitate us to predict probable melting origins in a polymeric DNA.

Lattice model of oligonucleotide hybridization in solution. I. Model and thermodynamics

A coarse-grained lattice model of DNA oligonucleotides is proposed to investigate the general mechanisms by which single-stranded oligonucleotides hybridize to their complementary strands in solution. The model, based on a high-coordination cubic lattice, is simple enough to allow the direct simulation of DNA solutions, yet capturing how the fundamental thermodynamic processes are microscopically encoded in the nucleobase sequences. Physically relevant interactions are considered explicitly, such as interchain excluded volume, anisotropic base-pairing and base-stacking, and single-stranded bending rigidity. The model is studied in detail by a specially adapted Monte Carlo simulation method, based on parallel tempering and biased trials, which is designed to overcome the entropic and enthalpic barriers associated with the sampling of hybridization events of multiple single-stranded chains in solution. This methodology addresses both the configurational complexity of bringing together two complementary strands in a favorable orientation (entropic barrier) and the energetic penalty of breaking apart multiple associated bases in a double-stranded state (enthalpic barrier). For strands with sequences restricted to nonstaggering association and homogeneous pairing and stacking energies, base-pairing is found to dominate the hybridization over the translational and conformational entropy. For strands with sequence-dependent pairing corresponding to that of DNA, the complex dependence of the model's thermal stability on concentration, sequence, and degree of complementarity is shown to be qualitatively and quantitatively consistent both with experiment and with the predictions of statistical mechanical models.

Stacking interactions in denaturation of DNA fragments

A mesoscopic model for heterogeneous DNA denaturation is developed in the framework of the path integral formalism. The base pair stretchings are treated as one-dimensional, time-dependent paths contributing to the partition function. The size of the paths ensemble, which measures the degree of cooperativity of the system, is computed versus temperature consistently with the model potential physical requirements. It is shown that the ensemble size strongly varies with the molecule backbone stiffness providing a quantitative relation between stacking and features of the melting transition. The latter is an overall smooth crossover which begins from the adenine-thymine-rich portions of the fragment. The harmonic stacking coupling shifts, along the T -axis, the occurrence of the multistep denaturation but it does not change the character of the crossover. The methods to compute the fractions of open base pairs versus temperature are discussed: by averaging the base pair displacements over the path ensemble, we find that such fractions signal the multisteps of the transition in good agreement with the indications provided by the specific heat plots.

Structural correlations and melting of B-DNA fibers

Despite numerous attempts, understanding the thermal denaturation of DNA is still a challenge due to the lack of structural data on the transition since standard experimental approaches to DNA melting are made in solution and do not provide spatial information. We report a measurement using neutron scattering from oriented DNA fibers to determine the size of the regions that stay in the double-helix conformation as the melting temperature is approached from below. A Bragg peak from the B form of DNA is observed as a function of temperature and its width and integrated intensity are measured. These results, complemented by a differential calorimetry study of the melting of B-DNA fibers as well as electrophoresis and optical observation data, are analyzed in terms of a one-dimensional mesoscopic model of DNA.

Modeling DNA response to terahertz radiation

Collective response of DNA to terahertz electric fields is studied in a simple pair bond model. We confirm, with some caveats, a previous observation of destabilizing DNA breather modes and explore the parameter dependence of these modes. It is shown that breather modes are eliminated under reasonable physical conditions and that thermal effects are significant.

Computing DNA duplex instability profiles efficiently with a two-state model: trends of promoters and binding sites

Background

DNA instability profiles have been used recently for predicting the transcriptional start site and the location of core promoters, and to gain insight into promoter action. It was also shown that the use of these profiles can significantly improve the performance of motif finding programs.

Results

In this work we introduce a new method for computing DNA instability profiles. The model that we use is a modified Ising-type model and it is implemented via statistical mechanics. Our linear time algorithm computes the profile of a 10,000 base-pair long sequence in less than one second. The method we use also allows the computation of the probability that several consecutive bases are unpaired simultaneously. This is a feature that is not available in other linear-time algorithms. We use the model to compare the thermodynamic trends of promoter sequences of several genomes. In addition, we report results that associate the location of local extrema in the instability profiles with the presence of core promoter elements at these locations and with the location of the transcription start sites (TSS). We also analyzed the instability scores of binding sites of several human core promoter elements. We show that the instability scores of functional binding sites of a given core promoter element are significantly different than the scores of sites with the same motif occurring outside the functional range (relative to the TSS).

Conclusions

The time efficiency of the algorithm and its genome-wide applications makes this work of broad interest to scientists interested in transcriptional regulation, motif discovery, and comparative genomics.

Thermal and mechanical properties of a DNA model with solvation barrier

We study the thermal and mechanical behaviors of DNA denaturation in the frame of the mesoscopic Peyrard-Bishop-Dauxois model with the inclusion of solvent interaction. By analyzing the melting transition of a homogeneous A-T sequence, we are able to set suitable values of the parameters of the model and study the formation and stability of bubbles in the system. Then, we focus on the case of the P5 promoter sequence and use the principal component analysis of the trajectories to extract the main information on the dynamical behavior of the system. We find that this analysis method gives an excellent agreement with previous biological results.

Denaturation patterns in heterogeneous DNA

The thermodynamical properties of heterogeneous DNA sequences are computed by path-integral techniques applied to a nonlinear model Hamiltonian. The base pairs relative displacements are interpreted as time-dependent paths whose amplitudes are consistent with the model potential for the hydrogen bonds between complementary strands. The portion of configuration space contributing to the partition function is determined, at any temperature, by selecting the ensemble of paths which fulfill the second law of thermodynamics. For a short DNA fragment, the denaturation is signaled by a succession of peaks in the specific-heat plots while the entropy grows continuously versus T. Thus, the opening of the double strand with bubble formation appears as a smooth crossover due to base pair fluctuation effects which are accounted for by the path-integral method. The multistep transition is driven by the adenine-thymine- (AT) rich regions of the DNA fragment. The base pairs path ensemble shows an enhanced degree of cooperativity at about the same temperatures for which the specific-heat peaks occur. These findings establish a link between microscopic and macroscopic signatures of the transition. The fractions of mean base pair stretchings are computed by varying the AT base pairs content and taking some threshold values for the occurrence of the molecule denaturation.

Localization and delocalization of energy in a Peyrard-Bishop chain

This paper studies energy localization conditions in lattices of the type proposed by Peyrard and Bishop. Homogeneous and inhomogeneous lattices are analyzed and the role of interfaces in the latter is emphasized. Simulations allowed us to identify critical energy values for the existence of localization. After a certain energy value, it is possible to observe the loss of energy localization along the chain.

Genome wide application of DNA melting analysis

Correspondences between functional and thermodynamic melting properties in a genome are being increasingly employed for ab initio gene finding and for the interpretation of the evolution of genomes. Here we present the first systematic genome wide comparison between biologically coding domains and thermodynamically stable regions. In particular, we develop statistical methods to estimate the reliability of the resulting predictions. Not surprisingly, we find that the success of the approach depends on the difference in GC content between the coding and the non-coding parts of the genome and on the percentage of coding base-pairs in the sequence. These prerequisites vary strongly between species, where we observe no systematic differences between eukaryotes and prokaryotes. We find a number of organisms in which the strong correlation of coding domains and thermodynamically stable regions allows us to identify putative exons or genes to complement existing approaches. In contrast to previous investigations along these lines we have not employed the Poland-Scheraga (PS) model of DNA melting but use the earlier Zimm-Bragg (ZB) model. The Ising-like form of the ZB model can be viewed as an approximation to the PS model, with averaged loop entropies included into the cooperative factor [Formula: see text]. This results in a speed-up by a factor of 20-100 compared to the Fixman-Freire algorithm for the solution of the PS model. We show that for genomic sequences the resulting systematic errors are negligible compared to the parameterization uncertainty of the models. We argue that for limited computing resources, available CPU power is better invested in broadening the statistical base for genomic investigations than in marginal improvements of the description of the physical melting behavior.

Experimental and theoretical studies of sequence effects on the fluctuation and melting of short DNA molecules

Understanding the melting of short DNA sequences probes DNA at the scale of the genetic code and raises questions which are very different from those posed by very long sequences, which have been extensively studied. We investigate this problem by combining experiments and theory. A new experimental method allows us to make a mapping of the opening of the guanines along the sequence as a function of temperature. The results indicate that non-local effects may be important in DNA because an AT-rich region is able to influence the opening of a base pair which is about 10 base pairs away. An earlier mesoscopic model of DNA is modified to correctly describe the timescales associated with the opening of individual base pairs well below melting, and to properly take into account the sequence. Using this model to analyze some characteristic sequences for which detailed experimental data on the melting is available (Montrichok et al 2003 Europhys. Lett. 62 452), we show that we have to introduce non-local effects of AT-rich regions to get acceptable results. This brings a second indication that the influence of these highly fluctuating regions of DNA on their neighborhood can extend to some distance.

Bubbles in DNA melting

This paper presents a review of largely our own work on the DNA melting transition, and some new measurements of the elastic energy of sharp bends in single-stranded DNA and RNA. The purpose is to present the point of view that studying the transition of intermediate size oligomers leads to valuable tests of the models, in particular the ingredients most important for a reduced-degrees-of-freedom description, such as the different role of base pairing and base stacking. We make the case that, with intermediate size oligomers, one can actually measure the bubble length, which exhibits a more interesting behavior than the fraction of dissociated bases alone. Here is where more work seems necessary, both on the experimental and the modeling side, to understand the differences between theory and experiments. We summarize our previous results on the cooperativity parameters, which suggest that the transition is never exactly two-state no matter how short the molecule, or in other words the nucleation size for bubbles opening at the ends of the molecule is essentially 1 base pair. We briefly discuss our own modification of the nearest-neighbor model which treats pairing and stacking separately, as a way to fit the experimental melting profiles in this intermediate length regime. Finally we go on to present some new measurements on the stability of DNA and RNA hairpins with very short loops.

Simplified Langevin approach to the Peyrard-Bishop-Dauxois model of DNA

A simple Langevin approach is used to study stationary properties of the Peyrard-Bishop-Dauxois model for DNA, allowing known properties to be recovered in an easy way. Results are shown for the denaturation transition in homogeneous samples, for which some implications, so far overlooked, of an analogy with equilibrium wetting transitions are highlighted. This analogy implies that the order parameter, asymptotically, exhibits a second-order transition even if it may be very abrupt for nonzero values of the stiffness parameter. Not surprisingly, we also find that, for heterogeneous DNA, within this model the largest bubbles in the premelting stage appear in adenine-thymine-rich regions, while we suggest the possibility of some sort of not strictly local effects owing to the merging of bubbles.

Theoretical investigation of finite size effects at DNA melting

We investigated how the finiteness of the length of a sequence affects the phase transition that takes place at the DNA melting temperature. For this purpose, we modified the transfer integral method to adapt it to the calculation of both extensive (partition function, entropy, specific heat, etc.) and nonextensive (order parameter and average separation between paired bases) thermodynamic quantities of finite sequences with open boundary conditions, and applied the modified procedure to two different dynamical models. We characterized in some detail the three effects that take place when the length of the sequence is decreased, namely, (i) the decrease of the critical temperature, (ii) the decrease of the peak values of all quantities that diverge at the thermodynamic limit but remain finite for finite sequences, like the specific heat and the correlation length, and (iii) the broadening of the temperature range over which the transition affects the dynamics of the system. We also performed a finite size scaling analysis of the two models and showed that the singular part of the free energy can indeed be expressed in terms of a homogeneous function. However, Josephson's identity is satisfied for none of the investigated models, so that the derivation of the characteristic exponents which appear, for example, in the expression of the specific heat requires some care.