Measurement of the phase diagram of DNA unzipping in the temperature-force plane

Sheetlike structure in the proximity of compact DNA

We determine the phase diagram of DNA with inter- and intrastrand native-pair interactions that mimic the compaction of DNA. We show that DNA takes an overall sheetlike structure in the region where an incipient transition to a compact phase would have occurred. The stability of this phase is due to the extra entropy from the folding of the sheet, which is absent in the remaining polymerlike states of the phase diagram.

Appearance of de Gennes length in force-induced transitions

Using Langevin dynamic simulations, a simple coarse-grained model of a DNA protein construct is used to study the DNA rupture and the protein unfolding. We identify three distinct states: (i) zipped DNA and collapsed protein, (ii) unzipped DNA and stretched protein, and (iii) unzipped DNA and collapsed protein. Here, we find a phase diagram that shows these states depending on the size of the DNA handle and the protein. For a less stable protein, unfolding is solely governed by the size of the linker DNA, whereas if the protein's stability increases, complete unfolding becomes impossible because the rupture force for DNA has reached a saturation regime influenced by the de Gennes length. We show that unfolding occurs via a few intermediate states by monitoring the force-extension curve of the entire protein. We extend our study to a heterogeneous protein system, where similar intermediate states in two systems can lead to different protein unfolding paths.

Stochastic resonance in a model of a periodically driven DNA: Multiple transitions, scaling, and sequence dependence

We numerically study stochastic resonance in the unzipping of a model double-stranded DNA by a periodic force. We observe multiple peaks in stochastic resonance in the output signal as the driving force frequency is varied for different force amplitudes, temperature, chain length, and chain heterogeneity. Multiple peaks point to the existence of multiple stable and metastable states, which correspond to dynamical states of partially zipped and unzipped conformations and transitions between them. We quantify such transitions by looking at the time evolution of the fraction of bound base pairs. We obtain phase diagrams in the force amplitude-temperature plane both in the resonance frequency of the primary peak and the output signal at the peak value. We further obtain an excellent scaling behavior of the output signal for changing lengths of the DNA. Resonance behavior is also affected by chain heterogeneity as it depends strongly on which base pair the periodic forcing is applied.

Force-induced melting of DNA hairpin: Unfolding pathways and phase diagrams

Using the exact enumeration technique, we have studied the force-induced melting of a DNA hairpin on the face centered cubic lattice for two different sequences which differ in terms of loop closing base pairs. The melting profiles obtained from the exact enumeration technique is consistent with the Gaussian network model and Langevin dynamics simulations. Probability distribution analysis based on the exact density of states revealed the microscopic details of the opening of the hairpin. We showed the existence of intermediate states near the melting temperature. We further showed that different ensembles used to model single-molecule force spectroscopy setups may give different force-temperature diagrams. We delineate the possible reasons for the observed discrepancies.

Hybrid nucleobases as new and efficient unnatural genetic letters

To expand the existing genetic letters beyond the natural four nucleotides, such as G, C, A, and T, it is necessary to design robust nucleotides that can not only produce stable and unperturbed DNA but also function naturally in living cells. Although hydrophobic bases, such as d5SICS (2,6-dimethyl-2H-isoquiniline-1-thione) and dNaM (2-methoxy-3-methylnaphthalene) were shown to be replicated in bacterial cells, the d5SICS:dNaM base-pair was found to perturb the structure of the duplex DNA. Therefore, it is necessary to design nucleobases that can form base pairs like the natural G:C and A:T pairs. Here, a reliable dispersion-corrected density functional theory has been used to design several nucleobases that can produce three-hydrogen-bonded base pairs like the G:C pair. In doing so, the Watson-Crick faces of d5SICS and dNaM were modified by replacing the hydrophobic groups with hydrogen bond donors and acceptors. As dNaM contains an unnatural C-glycosidic bond (C-dNaM), it was also modified to contain the natural N-glycosidic bond (N-dNaM). This technique produced 91 new bases (N-d5SICS-X (X = 1-33), C-dNaM-X (X = 1-35), and N-dNaM-X (X = 1-23), where X is the different types of modifications applied to d5SICS and dNaM) and 259 base-pairs. Among these base pairs, 76 base pairs are found to be more stable than the G:C pair. Interestingly, the N-d5SICS-32:C-dNaM-32 and N-d5SICS-32:N-dNaM-20 pairs are found to be the most stable with binding energies of about -28.0 kcal/mol. The base-pair patterns of these pairs are also analogous to that of the G:C pair. Hence, it is proposed that N-d5SICS-32, C-dNaM-32, and N-dNaM-20 would act as efficient new genetic letters to produce stable and unperturbed artificial DNA.Communicated by Ramaswamy H. Sarma.

Maturation of siRNA by strand separation: Steered molecular dynamics study

RNA interference, particularly siRNA induced gene silencing is becoming an important avenue of modern therapeutics. The siRNA is delivered to the cells as short double helical RNA which becomes single stranded for forming the RISC complex. Significant experimental evidence is available for most of the steps except the process of the separation of the two strands. We have attempted to understand the pathway for double stranded siRNA (dsRNA) to single stranded (ssRNA) molecules using steered molecular dynamics simulations. As the process is completely unexplored we have applied force from all possible directions restraining all possible residues to convert dsRNA to ssRNA. We found pulling one strand along the helical axis direction restraining the far end of the other strand demands excessive force for ssRNA formation. Pulling a central residue of one strand, in a direction perpendicular to the helix axis, while keeping the base paired residue fixed requires intermediate force for strand separation. Moreover, we found that in this process the force requirement is quite high for the first bubble formation (nucleation energy) and the bubble propagation energies are quite small. We believe the success rate of the design of siRNA sequences for gene silencing may increase if this mechanistic knowledge is utilized for such a design process.Communicated by Ramaswamy H. Sarma.

Antitumor drugs effect on the stability of double-stranded DNA: steered molecular dynamics analysis

Denaturation of the DNA double helix inside the cell is essential for cellular processes such as replication and transcription for the growth of the cells. However, the growth of unwanted cells, which are responsible for cancerous kind of disease, is one of the biggest challenges of modern therapeutics. DNA cross-linking agents may kill cancer cells by damaging their DNA and stopping them from dividing. In the present study, we have carried out steered molecular dynamics simulations to study the effects of rupture and unzipping forces on the stability of dsDNA in the absence and presence of covalently bonded drugs. We have found that the stability of dsDNA increases strongly in the presence of covalently bonded drugs. The microscopic study of disruption of hydrogen-bonds associated with base-pairs of the dsDNA and the study of the variation of stacking overlap parameters gives evidence of symmetry during the rupture and asymmetry in the unzip event. The significance of the mechanism of force-induced melting study of the dsDNA in the absence and presence of antitumor drugs might have a biological relevance as it provides a pathway to open the double helix in a specific position and may help for the pharmaceutical design of drugs.Communicated by Ramaswamy H. Sarma.

5-Methyl-cytosine stabilizes DNA but hinders DNA hybridization revealed by magnetic tweezers and simulations

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.

Hysteresis loop area scaling exponents in DNA unzipping by a periodic force: A Langevin dynamics simulation study

Using Langevin dynamics simulations, we study the hysteresis in unzipping of longer double-stranded DNA chains whose ends are subjected to a time-dependent periodic force with frequency ω and amplitude G keeping the other end fixed. We find that the area of the hysteresis loop, A_{loop}, scales as 1/ω at higher frequencies, whereas it scales as (G-G_{c})^{α}ω^{β} with exponents α=1 and β=1.25 in the low-frequency regime. These values are same as the exponents obtained in Monte Carlo simulation studies of a directed self-avoiding walk model of a homopolymer DNA [R. Kapri, Phys. Rev. E 90, 062719 (2014)10.1103/PhysRevE.90.062719], and the block copolymer DNA [R. K. Yadav and R. Kapri, Phys. Rev. E 103, 012413 (2021)2470-004510.1103/PhysRevE.103.012413] on a square lattice, and differs from the values reported earlier using Langevin dynamics simulation studies on a much shorter DNA hairpins.

Effect of loop sequence on unzipping of short DNA hairpins

The dependence of stability on the sequence of a DNA hairpin has been investigated through atomistic simulations. For this, a sequence of 16 bases of a hairpin, which consists of a loop of four bases and a stem of six base pairs, has been considered. We have taken eight different sequences, where the first five base pairs were kept fixed in all sequences, whereas the loop sequence and the identity of the duplex base pair closing the loop have been varied. For these hairpin structures, force-induced melting (unzipping) studies were carried out to investigate the effect of the variables on the stability of hairpin. The temperature at which half of the base pairs are open is termed the melting temperature. We defined the unzipping force F_{h} (half of the base pairs are open) and showed that it may not provide the effect of closing the base pair or loop sequence on the stability of the DNA hairpin. In order to have a better understanding of the stability of a DNA hairpin, the closing base pair or hairpin loop must be open. This requires complete opening of the stem. We defined a force F_{c} at which all base pairs of the stem are open, and we showed that the F_{c} gives better understanding of DNA hairpin stability.

Elasticity of a DNA chain dotted with bubbles under force

The flexibility and the extension along the direction of the force are shown to be related to the bubble number fluctuation and the average number of bubbles, respectively, when the strands of the DNA are subjected to a force along the same direction, here called a stretching force. The force-temperature phase diagram shows the existence of a tricritical point, where the first-order force-induced zipping transition becomes continuous. On the other hand, when the forces are being applied in opposite directions, here called an unzipping force, the transition remains first order, with the possibility of vanishing of the low-temperature reentrant phase for a semiflexible DNA. Moreover, we found that the bulk elasticity changes only if an external force penetrates the bound phase and affects the bubble states.

Anomalous Laterally Stressed Kinetically Trapped DNA Surface Conformations

DNA kinking is inevitable for the highly anisotropic 1D-1D electrostatic interaction with the one-dimensionally periodically charged surface. The double helical structure of the DNA kinetically trapped on positively charged monomolecular films comprising the lamellar templates is strongly laterally stressed and extremely perturbed at the nanometer scale. The DNA kinetic trapping is not a smooth 3D-> 2D conformational flattening but is a complex nonlinear in-plane mechanical response (bending, tensile and unzipping) driven by the physics beyond the scope of the applicability of the linear worm-like chain approximation. Up to now, the DNA molecule adsorbed on a surface was believed to always preserve its native structure. This belief implies a negligible contribution of lateral surface forces during and after DNA adsorption although their impact has never been elucidated. High-resolution atomic force microscopy was used to observe that stiff DNA molecules kinetically trapped on monomolecular films comprising one-dimensional periodically charged lamellar templates as a single layer or as a sublayer are oversaturated by sharp discontinuous kinks and can also be locally melted and supercoiled. We argue that kink/anti-kink pairs are induced by an overcritical lateral bending stress (> 30 pNnm) inevitable for the highly anisotropic 1D-1D electrostatic interaction of DNA and underlying rows of positive surface charges. In addition, the unexpected kink-inducing mechanical instability in the shape of the template-directed DNA confined between the positively charged lamellar sides is observed indicating the strong impact of helicity. The previously reported anomalously low values of the persistence length of the surface-adsorbed DNA are explained by the impact of the surface-induced low-scale bending. The sites of the local melting and supercoiling are convincingly introduced as other lateral stress-induced structural DNA anomalies by establishing a link with DNA high-force mechanics. The results open up the study in the completely unexplored area of the principally anomalous kinetically trapped DNA surface conformations in which the DNA local mechanical response to the surface-induced spatially modulated lateral electrostatic stress is essentially nonlinear. The underlying rich and complex in-plane nonlinear physics acts at the nanoscale beyond the scope of applicability of the worm-like chain approximation.

Unzipping of a double-stranded block copolymer DNA by a periodic force

Using Monte Carlo simulations, we study the hysteresis in unzipping of a double-stranded block copolymer DNA with -A_{n}B_{n}- repeat units. Here A and B represent two different types of base pairs having two and three bonds, respectively, and 2n represents the number of such base pairs in a unit. The end of the DNA are subjected to a time-dependent periodic force with frequency (ω) and amplitude (g_{0}) keeping the other end fixed. We find that the equilibrium force-temperature phase diagram for the static force is independent of the DNA sequence. For a periodic force case, the results are found to be dependent on the block copolymer DNA sequence and on the base pair type on which the periodic force is acting. We observe hysteresis loops of various shapes and sizes and obtain the scaling of loop area both at low- and high-frequency regimes.

Polymer in wedge-shaped confinement: Effect on the θ temperature

The equilibrium properties of a finite-length linear polymer chain confined in an infinite wedge composed of two perfectly reflecting hard walls meeting at a variable apex angle (α) are presented. One end of the polymer is anchored a distance y from the apex on the conical axis of symmetry, while the other end is free. We report here, the nonmonotonic behavior of θ temperature as a function of y for a finite-length chain. Data collapse for different chain lengths indicates that such behavior will exist for all finite lengths. We delineate the origin of such nonmonotonic behavior, which may have potential applications in understanding the cellular process occurring in nanoconfined geometries.

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.

A Universal Assay for Making DNA, RNA, and RNA-DNA Hybrid Configurations for Single-Molecule Manipulation in Two or Three Steps without Ligation

Despite having a great variety of topologies, most DNA, RNA, and RNA-DNA hybrid (RDH) configurations for single-molecule manipulation are composed of several single-stranded (ss) DNA and ssRNA strands, with functional labels at the two ends for surface tethering. On this basis, we developed a simple, robust, and universal amplification-annealing (AA) assay for making all these configurations in two or three steps without inefficient digestion and ligation reactions. As examples, we made ssDNA, short ssDNA with double-stranded (ds) DNA handles, dsDNA with ssDNA handles, replication-fork shaped DNA/RDH/RNA, DNA holiday junction, three-site multiple-labeled and nicked DNA, torsion-constrained RNA/RDH, and short ssRNA with RDH handles. In addition to single-molecule manipulation techniques including optical tweezers, magnetic tweezers, and atomic force microscopy, these configurations can be applied in other surface-tethering techniques as well.

Determination of the Thermodynamic Parameters of DNA Double Helix Unwinding with the Help of Mechanical Methods

For the first time, rupture event scanning (REVS) procedure based on quartz crystal microbalance (QCM) and involving only mechanical action was used to determine the height of the energy barrier for dsDNA unwinding. Melting point was determined with the help of this procedure. To determine the thermodynamic parameters including enthalpy, DNA denaturation was represented as a unimolecular process. This allowed us to recover the energy profiles from the experimental data obtained by force measurements at different scanning times (reaction times) for different temperatures. The thus obtained results were compared with the data obtained with the help of another mechanical method, namely, atomic force microscopy. The mechanism of DNA unwinding in QCM-based experiments through the unzipping mode, as proposed by us in previous works, was confirmed. Thus, we demonstrated that REVS procedure may be used to assess the thermodynamic parameters of dsDNA unwinding.

Manipulation of double-stranded DNA melting by force

By integrating elasticity-as described by the Gaussian network model-with bond binding energies that distinguish between different base-pair identities and stacking configurations, we study the force induced melting of a double-stranded DNA (dsDNA). Our approach is a generalization of our previous study of thermal dsDNA denaturation [J. Chem. Phys. 145, 144101 (2016)JCPSA60021-960610.1063/1.4964285] to that induced by force at finite temperatures. It allows us to obtain semimicroscopic information about the opening of the chain, such as whether the dsDNA opens from one of the ends or from the interior, forming an internal bubble. We study different types of force manipulation: (i) "end unzipping," with force acting at a single end base pair perpendicular to the helix, (ii) "midunzipping," with force acting at a middle base pair perpendicular to the helix, and (iii) "end shearing," where the force acts at opposite ends along the helix. By monitoring the free-energy landscape and probability distribution of intermediate denaturation states, we show that different dominant intermediate states are stabilized depending on the type of force manipulation used. In particular, the bubble state of the sequence L60B36, which we have previously found to be a stable state during thermal denaturation, is absent for end unzipping and end shearing, whereas very similar bubbles are stabilized by midunzipping, or when the force location is near the middle of the chain. Ours results offer a simple tool for stabilizing bubbles and loops using force manipulations at different temperatures, and may implicate on the mechanism in which DNA enzymes or motors open regions of the chain.

Ionic effects on the temperature-force phase diagram of DNA

Double-stranded DNA (dsDNA) undergoes a structural transition to single-stranded DNA (ssDNA) in many biologically important processes such as replication and transcription. This strand separation arises in response either to thermal fluctuations or to external forces. The roles of ions are twofold, shortening the range of the interstrand potential and renormalizing the DNA elastic modulus. The dsDNA-to-ssDNA transition is studied on the basis that dsDNA is regarded as a bound state while ssDNA is regarded as an unbound state. The ground state energy of DNA is obtained by mapping the statistical mechanics problem to the imaginary time quantum mechanics problem. In the temperature-force phase diagram the critical force F _c (T) increases logarithmically with the Na⁺ concentration in the range from 32 to 110 mM. Discussing this logarithmic dependence of F _c (T) within the framework of polyelectrolyte theory, it inevitably suggests a constraint on the difference between the interstrand separation and the length per unit charge during the dsDNA-to-ssDNA transition.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Multiplexed single-molecule force spectroscopy using a centrifuge

We present a miniature centrifuge force microscope (CFM) that repurposes a benchtop centrifuge for high-throughput single-molecule experiments with high-resolution particle tracking, a large force range, temperature control and simple push-button operation. Incorporating DNA nanoswitches to enable repeated interrogation by force of single molecular pairs, we demonstrate increased throughput, reliability and the ability to characterize population heterogeneity. We perform spatiotemporally multiplexed experiments to collect 1,863 bond rupture statistics from 538 traceable molecular pairs in a single experiment, and show that 2 populations of DNA zippers can be distinguished using per-molecule statistics to reduce noise.

Single molecule force spectroscopy (SMFS) has limitations in throughput and the ability to repeatedly interrogate single bonds. Here the authors repurpose a benchtop centrifuge and use DNA nanoswitches to enable high throughput SMFS capable of repeatedly measuring forces of single molecular pairs.

Free-energy landscape and characteristic forces for the initiation of DNA unzipping

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.

Force-Induced Rupture of a DNA Duplex: From Fundamentals to Force Sensors

The rupture of double-stranded DNA under stress is a key process in biophysics and nanotechnology. In this article, we consider the shear-induced rupture of short DNA duplexes, a system that has been given new importance by recently designed force sensors and nanotechnological devices. We argue that rupture must be understood as an activated process, where the duplex state is metastable and the strands will separate in a finite time that depends on the duplex length and the force applied. Thus, the critical shearing force required to rupture a duplex depends strongly on the time scale of observation. We use simple models of DNA to show that this approach naturally captures the observed dependence of the force required to rupture a duplex within a given time on duplex length. In particular, this critical force is zero for the shortest duplexes, before rising sharply and then plateauing in the long length limit. The prevailing approach, based on identifying when the presence of each additional base pair within the duplex is thermodynamically unfavorable rather than allowing for metastability, does not predict a time-scale-dependent critical force and does not naturally incorporate a critical force of zero for the shortest duplexes. We demonstrate that our findings have important consequences for the behavior of a new force-sensing nanodevice, which operates in a mixed mode that interpolates between shearing and unzipping. At a fixed time scale and duplex length, the critical force exhibits a sigmoidal dependence on the fraction of the duplex that is subject to shearing.

Unzipping DNA by a periodic force: hysteresis loop area and its scaling

Using Monte Carlo simulations, we study the hysteresis in the unzipping of double-stranded DNA whose ends are subjected to a time-dependent periodic force with frequency (ω) and amplitude (G). For the static force, i.e., ω→0, the DNA is in equilibrium with no hysteresis. On increasing ω, the area of the hysteresis loop initially increases and becomes maximum at frequency ω*(G), which depends on the force amplitude G. If the frequency is increased further, we find that for lower amplitudes the loop area decreases monotonically to zero, but for higher amplitudes it has an oscillatory component. The height of subsequent peaks decreases, and finally the loop area becomes zero at very high frequencies. The number of peaks depends on the length of the DNA. We give a simple analysis to estimate the frequencies at which maxima and minima occur in the loop area. We find that the area of the hysteresis loop scales as 1/ω in the high-frequency regime, whereas it scales as G(α)ω(β) with exponents α=1 and β=5/4 at low frequencies. The values of the exponents α and β are different from the exponents reported earlier based on the hysteresis of small hairpins.

Statistical mechanics of DNA rupture: theory and simulations

We study the effects of the shear force on the rupture mechanism on a double stranded DNA. Motivated by recent experiments, we perform the atomistic simulations with explicit solvent to obtain the distributions of extension in hydrogen and covalent bonds below the rupture force. We obtain a significant difference between the atomistic simulations and the existing results in the literature based on the coarse-grained models (theory and simulations). We discuss the possible reasons and improve the coarse-grained model by incorporating the consequences of semi-microscopic details of the nucleotides in its description. The distributions obtained by the modified model (simulations and theoretical) are qualitatively similar to the one obtained using atomistic simulations.

Statistical mechanics of DNA unzipping under periodic force: scaling behavior of hysteresis loops

A simple model of DNA based on two interacting polymers has been used to study the unzipping of a double stranded DNA subjected to a periodic force. We propose a dynamical transition where, without changing the physiological condition, it is possible to bring DNA from the zipped or unzipped state to a new dynamic (hysteretic) state by varying the frequency of the applied force. Our studies reveal that the area of the hysteresis loop grows with the same exponents as of the isotropic spin systems. These exponents are amenable to verification in the force spectroscopic experiments.

Dynamical phase transition of a periodically driven DNA

Replication and transcription are two important processes in living systems. To execute such processes, various proteins work far away from equilibrium in a staggered way. Motivated by this, aspects of hysteresis during unzipping of DNA under a periodic drive are studied. A steady-state phase diagram of a driven DNA is proposed which is experimentally verifiable. As a two-state system, we also compare the results of DNA with that of an Ising magnet under an asymmetrical variation of the magnetic field.

Hysteresis and nonequilibrium work theorem for DNA unzipping

We study by using Monte Carlo simulations the hysteresis in unzipping and rezipping of a double stranded DNA (dsDNA) by pulling its strands in opposite directions in the fixed force ensemble. The force is increased at a constant rate from an initial value g(0) to some maximum value g(m) that lies above the phase boundary and then decreased back again to g(0). We observed hysteresis during a complete cycle of unzipping and rezipping. We obtained probability distributions of work performed over a cycle of unzipping and rezipping for various pulling rates. The mean of the distribution is found to be close (the difference being within 10%, except for very fast pulling) to the area of the hysteresis loop. We extract the equilibrium force versus separation isotherm by using the work theorem on repeated nonequilibrium force measurements. Our method is capable of reproducing the equilibrium and the nonequilibrium force-separation isotherms for the spontaneous rezipping of dsDNA.

Micromechanics of base pair unzipping in the DNA duplex

All-atom molecular dynamics (MD) simulations of DNA duplex unzipping in a water environment were performed. The investigated DNA double helix consists of a Drew-Dickerson dodecamer sequence and a hairpin (AAG) attached to the end of the double-helix chain. The considered system is used to examine the process of DNA strand separation under the action of an external force. This process occurs in vivo and now is being intensively investigated in experiments with single molecules. The DNA dodecamer duplex is consequently unzipped pair by pair by means of the steered MD. The unzipping trajectories turn out to be similar for the duplex parts with G·C content and rather distinct for the parts with A·T content. It is shown that during the unzipping each pair experiences two types of motion: relatively quick rotation together with all the duplex and slower motion in the frame of the unzipping fork. In the course of opening, the complementary pair passes through several distinct states: (i) the closed state in the double helix, (ii) the metastable preopened state in the unzipping fork and (iii) the unbound state. The performed simulations show that water molecules participate in the stabilization of the metastable states of the preopened base pairs in the DNA unzipping fork.

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.

Role of loop entropy in the force induced melting of DNA hairpin

Dynamics of a single stranded DNA, which can form a hairpin have been studied in the constant force ensemble. Using Langevin dynamics simulations, we obtained the force-temperature diagram, which differs from the theoretical prediction based on the lattice model. Probability analysis of the extreme bases of the stem revealed that at high temperature, the hairpin to coil transition is entropy dominated and the loop contributes significantly in its opening. However, at low temperature, the transition is force driven and the hairpin opens from the stem side. It is shown that the elastic energy plays a crucial role at high force. As a result, the force-temperature diagram differs significantly with the theoretical prediction.

Effect of shear force on the separation of double-stranded DNA

Using the Langevin dynamics simulation, we have studied the effects of shear force on the rupture of short double-stranded DNA at different temperatures. We show that the rupture force increases linearly with chain length and approaches the asymptotic value in accordance with the experiment. The qualitative nature of these curves remains almost the same for different temperatures but with a shift in the force. We observe three different regimes in the extension of covalent bonds (backbone) under shear force.

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.

Capture and manipulation of hybrid DNAs by carbon nanotube bundles

When approached from both sides, a piece of single-duplex-single DNA may be drawn into the inlets of two bundles of carbon nanotubes. This provides opportunities to manipulate the DNA by two bundles of nanotubes. The capture and manipulation processes envisaged above are simulated by molecular dynamics in this work. The radius of the carbon nanotube and the ambient temperature show the effects on the spontaneous insertion of DNA strands. This procedure, if successful, could be used for capturing expectant sdsDNAs, with subsequent manipulation to pull or to unzip the captured DNA.

Lambda exonuclease digestion of CGG trinucleotide repeats

Fragile X syndrome and other trinucleotide diseases are characterized by an elongation of a repeating DNA triplet. The ensemble-averaged lambda exonuclease digestion rate of different substrates, including one with an elongated FMR1 gene containing 120 CGG repeats, was measured using absorption and fluorescence spectroscopy. By use of magnetic tweezers sequence-dependent digestion rates and pausing was measured for individual lambda exonucleases. Within the triplet repeats a lower average and narrower distribution of rates and a higher frequency of pausing was observed.

Polymers in crowded environment under stretching force: Globule-coil transitions

We study flexible polymer macromolecules in a crowded (porous) environment, modeling them as self-attracting self-avoiding walks on site-diluted percolative lattices in space dimensions d=2,3 . The influence of stretching force on the polymer folding and the properties of globule-coil transitions are analyzed. Applying the pruned-enriched Rosenbluth chain-growth method, we estimate the transition temperature TTheta between collapsed and extended polymer configurations and construct the phase diagrams of the globule-coil coexistence when varying temperature and stretching force. The transition to a completely stretched state, caused by applying force, is discussed as well.

Nonlinearity in DNA and its relation to specific functions

In memory of Alwyn Scott, we discuss the connection between the nonlinear dynamics of double-stranded DNA, experimental findings, and specific DNA functions. We begin by discussing how thermally induced localized openings (bubbles) of the DNA double-strand are important for interpreting dynamic force spectroscopy data. Then, we demonstrate a correlation between the sequence-dependent propensity for bubble formation and transcription initiation and other regulatory effects in viral DNA. Finally, we discuss the possibility of a connection between DNA dynamics and the ability of repair proteins to recognize ultraviolet-radiation damage.

The structure of DNA overstretched from the 5'5' ends differs from the structure of DNA overstretched from the 3'3' ends

It has been suggested that the structure that results when double-stranded DNA (dsDNA) is pulled from the 3'3' ends differs from that which results when it is pulled from the 5'5' ends. In this work, we demonstrate, using lambda phage dsDNA, that the overstretched states do indeed show different properties, suggesting that they correspond to different structures. For 3'3' pulling versus 5'5' pulling, the following differences are observed: (i) the forces at which half of the molecules in the ensemble have made a complete force-induced transition to single stranded DNA are 141 +/- 3 pN and 122 +/- 4 pN, respectively; (ii) the extension vs. force curve for overstretched DNA has a marked change in slope at 127 +/- 3 pN for 3'3' and 110 +/- 3 pN for 5'5'; (iii) the hysteresis (H) in the extension vs. force curves at 150 mM NaCl is 0.3 +/- 0.8 pN microm for 3'3' versus 13 +/- 8 pN for 5'5'; and (iv) 3'3' and 5'5' molecules show different changes in hysteresis due to interactions with beta-cyclodextrin, a molecule that is known to form stable host-guest complexes with rotated base pairs, and glyoxal that is known to bind stably to unpaired bases. These differences and additional findings are well-accommodated by the corresponding structures predicted on theoretical grounds.

Looking into DNA breathing dynamics via quantum physics

We study generic aspects of bubble dynamics in DNA under time-dependent perturbations, for example, temperature change, by mapping the associated Fokker-Planck equation to a quantum time-dependent Schrödinger equation with imaginary time. In the static case we show that the eigenequation is exactly the same as that of the beta-deformed nuclear liquid drop model, without the issue of noninteger angular momentum. A universal breathing dynamics is demonstrated by using an approximate method in quantum mechanics. The calculated bubble autocorrelation function qualitatively agrees with experimental data. Under time-dependent modulations, utilizing the adiabatic approximation, bubble properties reveal memory effects.

DNA unzipping phase diagram calculated via replica theory

We show how single-molecule unzipping experiments can provide strong evidence that the zero-force melting transition of long molecules of natural dsDNA should be classified as a phase transition of the higher-order type (continuous). Toward this end, we study a statistical-mechanics model for the fluctuating structure of a long molecule of dsDNA, and compute the equilibrium phase diagram for the experiment in which the molecule is unzipped under applied force. We consider a perfect-matching dsDNA model, in which the loops are volume-excluding chains with arbitrary loop exponent c . We include stacking interactions, hydrogen bonds, and main-chain entropy. We include sequence heterogeneity at the level of random sequences; in particular, there is no correlation in the base-pairing (bp) energy from one sequence position to the next. We present heuristic arguments to demonstrate that the low-temperature macrostate does not exhibit degenerate ergodicity breaking. We use this claim to understand the results of our replica-theoretic calculation of the equilibrium properties of the system. As a function of temperature, we obtain the minimal force at which the molecule separates completely. This critical-force curve is a line in the temperature-force phase diagram that marks the regions where the molecule exists primarily as a double helix versus the region where the molecule exists as two separate strands. We compare our random-sequence model to magnetic tweezer experiments performed on the 48 502 bp genome of bacteriophage lambda . We find good agreement with the experimental data, which is restricted to temperatures between 24 and 50 degrees C . At higher temperatures, the critical-force curve of our random-sequence model is very different for that of the homogeneous-sequence version of our model. For both sequence models, the critical force falls to zero at the melting temperature T_{c} like |T-T_{c}|;{alpha} . For the homogeneous-sequence model, alpha=1/2 almost exactly, while for the random-sequence model, alpha approximately 0.9 . Importantly, the shape of the critical-force curve is connected, via our theory, to the manner in which the helix fraction falls to zero at T_{c} . The helix fraction is the property that is used to classify the melting transition as a type of phase transition. In our calculation, the shape of the critical-force curve holds strong evidence that the zero-force melting transition of long natural dsDNA should be classified as a higher-order (continuous) phase transition. Specifically, the order is 3rd or greater.