Single DNA denaturation and bubble dynamics

Stretching multistate flexible chains and loops

Polymer loop structure commonly appears in biological phenomena, such as DNA looping and DNA denaturation. When a chain forms a loop, its elastic behavior differs from that of an open chain due to the loss of entropy. In the case of reversible loop formation, interesting behavior emerges related to the multistate nature of the conformations. In this study, we model a multistate reversible loop as a looping Gaussian chain, which can bind (close) reversibly at one or several points to form a loop, or a zipping Gaussian loop, which can zip reversibly to form a double-stranded chain. For each model, we calculate the force-extension relations in the fixed-extension (Helmholtz) and the fixed-force (Gibbs) statistical ensembles. Unlike the single Gaussian chain or loop, the multilevel systems demonstrate qualitatively distinct tensile elasticity and ensemble inequivalence. In addition, we investigate a Gaussian necklace consisting of reversible alternating blocks of the zipped chain and loop and obtain the force-temperature phase diagram. The phase diagram implies a force-induced phase transition from a completely looped (denatured) state to a mixed (bound) state.

Influence of a Single Deuterium Substitution for Protium on the Frequency Generation of Different-Size Bubbles in IFNA17

The influence of a single 2H/1H replacement on the frequency generation of different-size bubbles in the human interferon alpha-17 gene (IFNA17) under various energies was studied by a developed algorithm and mathematical modeling without simplifications or averaging. This new approach showed the efficacy of researching DNA bubbles and open states both when all hydrogen bonds in nitrogenous base pairs are protium and after an 2H-substitution. After a single deuterium substitution under specific energies, it was demonstrated that the non-coding region of IFNA17 had a more significant regulatory role in bubble generation in the whole gene than the promoter had. It was revealed that a single deuterium substitution for protium has an influence on the frequency generation of DNA bubbles, which also depends on their size and is always higher for the smaller bubbles under the largest number of the studied energies. Wherein, compared to the natural condition under the same critical value of energy, the bigger raises of the bubble frequency occurrence (maximums) were found for 11-30 base pair (bp) bubbles (higher by 319%), 2-4 bp bubbles (higher by 300%), and 31 bp and over ones (higher by 220%); whereas the most significant reductions of the indicators (minimums) were observed for 11-30 bp bubbles (lower by 43%) and bubbles size over 30 bp (lower by 82%). In this study, we also analyzed the impact of several circumstances on the AT/GC ratio in the formation of DNA bubbles, both under natural conditions and after a single hydrogen isotope exchange. Moreover, based on the obtained data, substantial positive and inverse correlations were revealed between the AT/GC ratio and some factors (energy values, size of DNA bubbles). So, this modeling and variant of the modified algorithm, adapted for researching DNA bubbles, can be useful to study the regulation of replication and transcription in the genes under different isotopic substitutions in the nucleobases.

Melting of confined DNA: static and dynamic properties

We study dsDNA (double strand DNA) melting in detail within varying strip-like confinement in a two-dimensional lattice model. The interplay between reduced configurational entropy and attractive base-pairing energy results in a non-monotonic melting profile of DNA. Structural transitions associated with confined DNA melting reveal a stretched or extended state for very strong confinement. By using the exact enumeration method, we investigate the emergence of a local denatured zone e.g. bubbles during DNA melting. The survival time of a single bubble within varying strip width is studied from the Fokker-Planck formalism by considering the bubble size as a reaction co-ordinate. We show that a simple lattice model can capture the sequence heterogeneity effect on DNA melting and bubble dynamics within the strip. Different time scales of bubble zipping for different DNA sequences are found, which may have potential applications in denaturation mapping.

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.

Dynamical control of denaturation bubble nucleation in supercoiled DNA minicircles

We examine the behavior of supercoiled DNA minicircles containing between 200 and 400 base-pairs, also named microDNA, in which supercoiling favors thermally assisted DNA denaturation bubbles of nanometer size and controls their lifetime. Mesoscopic modeling and accelerated dynamics simulations allow us to overcome the limitations of atomistic simulations encountered in such systems, and offer detailed insight into the thermodynamic and dynamical properties associated with the nucleation and closure mechanisms of long-lived thermally assisted denaturation bubbles which do not stem from bending- or torque-driven stress. Suitable tuning of the degree of supercoiling and size of specifically designed microDNA is observed to lead to the control of opening characteristic times in the millisecond range, and closure characteristic times ranging over well distinct timescales, from microseconds to several minutes. We discuss how our results can be seen as a dynamical bandwidth which might enhance selectivity for specific DNA binding proteins.

DNA breathing dynamics under periodic forcing: Study of several distribution functions of relevant Brownian functionals

In this paper, we study DNA breathing dynamics in the presence of an external periodic force by proposing and inspecting several probability distribution functions (PDFs) of relevant Brownian functionals which specify the bubble lifetime, reactivity, and average size. We model the bubble dynamics process by an overdamped Langevin equation of broken base pairs on the Poland-Scheraga free energy landscape. Introducing an effective time-independent description for timescales larger than T[over ̃]=2π/ω (where ω is the frequency of external periodic force) and using an elegant backward Fokker-Planck method we derive closed form expressions of several PDFs associated with such stochastic processes. For instance, with an initial bubble size of x_{0}, we derive the following analytical expressions: (i) the PDF P(t_{f}|x_{0}) of the first passage time t_{f} which specifies the lifetime of the DNA breathing process, (ii) the PDF P(A|x_{0}) of the area A until the first passage time, and it provides much valuable information about the average bubble size and reactivity of the process, and (iii) the PDF P(M) associated with the maximum bubble size M of the breathing process before complete denaturation. Our analysis is limited to two limits: (a) large bubble size and (b) small bubble size. We further confirm our analytical predictions by computing the same PDFs with direct numerical simulations of the corresponding Langevin equations. We obtain very good agreement of our theoretical predictions with the numerically simulated results. Finally, several nontrivial scaling behaviors in the asymptotic limits for the above-mentioned PDFs are predicted, which can be verified further from experimental observation. Our main conclusion is that the large bubble dynamics is unaffected by the rapidly oscillating force, but the small bubble dynamics is significantly affected by the same periodic force.

Freezing shortens the lifetime of DNA molecules under tension

DNA samples are commonly frozen for storage. However, freezing can compromise the integrity of DNA molecules. Considering the wide applications of DNA molecules in nanotechnology, changes to DNA integrity at the molecular level may cause undesirable outcomes. However, the effects of freezing on DNA integrity have not been fully explored. To investigate the impact of freezing on DNA integrity, samples of frozen and non-frozen bacteriophage lambda DNA were studied using optical tweezers. Tension (5–35 pN) was applied to DNA molecules to mimic mechanical interactions between DNA and other biomolecules. The integrity of the DNA molecules was evaluated by measuring the time taken for single DNA molecules to break under tension. Mean lifetimes were determined by maximum likelihood estimates and variances were obtained through bootstrapping simulations. Under 5 pN of force, the mean lifetime of frozen samples is 44.3 min with 95% confidence interval (CI) between 36.7 min and 53.6 min while the mean lifetime of non-frozen samples is 133.2 min (95% CI: 97.8–190.1 min). Under 15 pN of force, the mean lifetimes are 10.8 min (95% CI: 7.6–12.6 min) and 78.5 min (95% CI: 58.1–108.9 min). The lifetimes of frozen DNA molecules are significantly reduced, implying that freezing compromises DNA integrity. Moreover, we found that the reduced DNA structural integrity cannot be restored using regular ligation process. These results indicate that freezing can alter the structural integrity of the DNA molecules.

Twist-stretch profiles of DNA chains

Helical molecules change their twist number under the effect of a mechanical load. We study the twist-stretch relation for a set of short DNA molecules modeled by a mesoscopic Hamiltonian. Finite temperature path integral techniques are applied to generate a large ensemble of possible configurations for the base pairs of the sequence. The model also accounts for the bending and twisting fluctuations between adjacent base pairs along the molecules stack. Simulating a broad range of twisting conformation, we compute the helix structural parameters by averaging over the ensemble of base pairs configurations. The method selects, for any applied force, the average twist angle which minimizes the molecule's free energy. It is found that the chains generally over-twist under an applied stretching and the over-twisting is physically associated to the contraction of the average helix diameter, i.e. to the damping of the base pair fluctuations. Instead, assuming that the maximum amplitude of the bending fluctuations may decrease against the external load, the DNA molecule first over-twists for weak applied forces and then untwists above a characteristic force value. Our results are discussed in relation to available experimental information albeit for kilo-base long molecules.

The detection of T-Nos, a genetic element present in GMOs, by cross-priming isothermal amplification with real-time fluorescence

An isothermal cross-priming amplification (CPA) assay for Agrobacterium tumefaciens nopaline synthase terminator (T-Nos) was established and investigated in this work. A set of six specific primers, recognizing eight distinct regions on the T-Nos sequence, was designed. The CPA assay was performed at a constant temperature, 63 °C, and detected by real-time fluorescence. The results indicated that real-time fluorescent CPA had high specificity, and the limit of detection was 1.06 × 10(3) copies of rice genomic DNA, which could be detected in 40 min. Comparison of real-time fluorescent CPA and conventional polymerase chain reaction (PCR) was also performed. Results revealed that real-time fluorescent CPA had a comparable sensitivity to conventional real-time PCR and had taken a shorter time. In addition, different contents of genetically modified (GM)-contaminated rice seed powder samples were detected for practical application. The result showed real-time fluorescent CPA could detect 0.5 % GM-contaminated samples at least, and the whole reaction could be finished in 35 min. Real-time fluorescent CPA is sensitive enough to monitor labeling systems and provides an attractive method for the detection of GMO.

Magnetic tweezers measurements of the nanomechanical stability of DNA against denaturation at various conditions of pH and ionic strength

The opening of DNA double strands is extremely relevant to several biological functions, such as replication and transcription or binding of specific proteins. Such opening phenomenon is particularly sensitive to the aqueous solvent conditions in which the DNA molecule is dispersed, as it can be observed by considering the classical dependence of DNA melting temperature on pH and salt concentration. In the present work, we report a single-molecule study of the stability of DNA against denaturation when subjected to changes in solvent. We investigated the appearance of DNA instability under specific external applied force and imposed twist values, which was revealed by an increase in the temporal fluctuations in the DNA extension. These fluctuations occur in the presence of a continuous interval of equilibrium states, ranging from a plectonemic state to a state characterized by denaturation bubbles. In particular, we observe the fluctuations only around a characteristic force value. Moreover, this characteristic force is demonstrated to be notably sensitive to variations in the pH and ionic strength. Finally, an extension of a theoretical model of plectoneme formation is used to estimate the average denaturation energy, which is found to be linearly correlated to the melting temperature of the DNA double strands.

Single-molecule study of the DNA denaturation phase transition in the force-torsion space

We use the "magnetic tweezers" technique to show the structural transitions that the DNA undergoes in the force-torsion space. In particular, we focus on the regions corresponding to negative supercoiling. These regions are characterized by the formation of the so-called denaturation bubbles, which play an essential role in the replication and transcription of DNA. We experimentally map the region of the force-torsion space where the denaturation takes place. We observe that large fluctuations in DNA extension occur at one of the boundaries of this region, i.e., when the formation of denaturation bubbles and of plectonemes compete. To describe the experiments, we introduce a suitable extension of the classical model. The model correctly describes the position of the denaturation regions, the transition boundaries, and the measured values of the DNA extension fluctuations.

Anharmonic stacking in supercoiled DNA

Multistep denaturation in a short circular DNA molecule is analyzed by a mesoscopic Hamiltonian model which accounts for the helicoidal geometry. Computation of melting profiles by the path integral method suggests that stacking anharmonicity stabilizes the double helix against thermal disruption of the hydrogen bonds. Twisting is essential in the model to capture the importance of nonlinear effects on the thermodynamical properties. In a ladder model with zero twist, anharmonic stacking scarcely affects the thermodynamics. Moderately untwisted helices, with respect to the equilibrium conformation, show an energetic advantage against the overtwisted ones. Accordingly moderately untwisted helices better sustain local fluctuational openings and make more unlikely the thermally driven complete strand separation.

Cross priming amplification: mechanism and optimization for isothermal DNA amplification

CPA is a class of isothermal amplification reactions that is carried out by a strand displacement DNA polymerase and does not require an initial denaturation step or the addition of a nicking enzyme. At the assay temperature of 63°C, the formation of a primer-template hybrid at transient, spontaneous denaturation bubbles in the DNA template is favored over re-annealing of the template strands by the high concentration of primer relative to template DNA. Strand displacement is encouraged by the annealing of cross primers with 5' ends that are not complementary to the template strand and the binding of a displacement primer upstream of the crossing primer. The resulting exponential amplification of target DNA is highly specific and highly sensitive, producing amplicons from as few as four bacterial cells. Here we report on the basic CPA mechanism - single crossing CPA - and provide details on alternative mechanisms.

Bubble solitons in an inhomogeneous, helical DNA molecular chain with flexible strands

Base pair opening in an inhomogeneous, DNA double helical molecular chain with flexible strands is investigated by studying its internal dynamics. For the study, a generalized model which takes into account the energies involved in stacking and hydrogen bonds along with inhomogeneity, helicity, and phonons coupled to the stacking and hydrogen bonds is proposed. The internal dynamics of the proposed DNA model is governed by a perturbed nonlinear Schrödinger equation. The unperturbed, completely integrable nonlinear Schrödinger equation which admits soliton solutions and forming a bubble corresponds to DNA dynamics with homogeneous and rigid strands. The results of the soliton perturbation analysis show that the inhomogeneity in stacking and hydrogen bonds in localized and periodic forms and the helicity do not alter the amplitude under perturbation. However, the flexibility of the strands diminishes the perturbed amplitude. On the other hand, the velocity of the soliton and bubble are unaltered due to all the above effects. However, the position and phase of the soliton and the bubble vary linearly in time. While the position of the soliton depends on the initial velocity, the phase depends on both the initial velocity and the initial amplitude of the soliton. The above effects introduce small fluctuation in the tail of the soliton, without affecting the robust nature of the soliton and the bubble during propagation. The soliton and the bubble obtained as solutions of the internal dynamics of the DNA molecule represent an opening of the base pairs which is essential for the transcription process.

DNA breathing dynamics: analytic results for distribution functions of relevant Brownian functionals

We investigate DNA breathing dynamics by suggesting and examining several Brownian functionals associated with bubble lifetime and reactivity. Bubble dynamics is described as an overdamped random walk in the number of broken base pairs. The walk takes place on the Poland-Scheraga free-energy landscape. We suggest several probability distribution functions that characterize the breathing process, and adopt the recently studied backward Fokker-Planck method and the path decomposition method as elegant and flexible tools for deriving these distributions. In particular, for a bubble of an initial size x₀, we derive analytical expressions for (i) the distribution P(t{f}|x₀) of the first-passage time t{f}, characterizing the bubble lifetime, (ii) the distribution P(A|x₀) of the area A until the first-passage time, providing information about the effective reactivity of the bubble to processes within the DNA, (iii) the distribution P(M) of the maximum bubble size M attained before the first-passage time, and (iv) the joint probability distribution P(M,t{m}) of the maximum bubble size M and the time t{m} of its occurrence before the first-passage time. These distributions are analyzed in the limit of small and large bubble sizes. We supplement our analytical predictions with direct numericalsimulations of the related Langevin equation, and obtain a very good agreement in the appropriate limits. The nontrivial scaling behavior of the various quantities analyzed here can, in principle, be explored experimentally.

Label-free single-molecule detection of DNA-hybridization kinetics with a carbon nanotube field-effect transistor

Single-molecule measurements of biomolecules can provide information about the molecular interactions and kinetics that are hidden in ensemble measurements. However, there is a requirement for techniques with improved sensitivity and time resolution for use in exploring biomolecular systems with fast dynamics. Here, we report the detection of DNA hybridization at the single-molecule level using a carbon nanotube field-effect transistor. By covalently attaching a single-stranded probe DNA sequence to a point defect in a carbon nanotube, we are able to measure two-level fluctuations in the conductance of the nanotube in the presence of a complementary DNA target. The kinetics of the system are studied as a function of temperature, allowing the measurement of rate constants, melting curves and activation energies for different sequences and target concentrations. The kinetics demonstrate non-Arrhenius behaviour, in agreement with DNA hybridization experiments using fluorescence correlation spectroscopy. This technique is label-free and could be used to probe single-molecule dynamics at microsecond timescales.