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

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

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.

Understanding protein diffusion on force-induced stretched DNA conformation

DNA morphology is subjected to environmental conditions and is closely coupled with its function. For example, DNA experiences stretching forces during several biological processes, including transcription and genome transactions, that significantly alter its conformation from that of B-DNA. Indeed, a well-defined 1.5 times extended conformation of dsDNA, known as Σ-DNA, has been reported in DNA complexes with proteins such as Rad51 and RecA. A striking feature in Σ-DNA is that the nucleobases are partitioned into triplets of three locally stacked bases separated by an empty rise gap of ∼ 5 Å. The functional role of such a DNA base triplet was hypothesized to be coupled with the ease of recognition of DNA bases by DNA-binding proteins (DBPs) and the physical origin of three letters (codon/anti-codon) in the genetic code. However, the underlying mechanism of base-triplet formation and the ease of DNA base-pair recognition by DBPs remain elusive. To investigate, here, we study the diffusion of a protein on a force-induced stretched DNA using coarse-grained molecular dynamics simulations. Upon pulling at the 3' end of DNA by constant forces, DNA exhibits a conformational transition from B-DNA to a ladder-like S-DNA conformation via Σ-DNA intermediate. The resulting stretched DNA conformations exhibit non-uniform base-pair clusters such as doublets, triplets, and quadruplets, of which triplets are energetically more stable than others. We find that protein favors the triplet formation compared to its unbound form while interacting non-specifically along DNA, and the relative population of it governs the ruggedness of the protein-DNA binding energy landscape and enhances the efficiency of DNA base recognition. Furthermore, we analyze the translocation mechanism of a DBP under different force regimes and underscore the significance of triplet formation in regulating the facilitated diffusion of protein on DNA. Our study, thus, provides a plausible framework for understanding the structure-function relationship between triplet formation and base recognition by a DBP and helps to understand gene regulation in complex regulatory processes.

Elucidating the Role of Topological Constraint on the Structure of Overstretched DNA Using Fluorescence Polarization Microscopy

The combination of DNA force spectroscopy and polarization microscopy of fluorescent DNA intercalator dyes can provide valuable insights into the structure of DNA under tension. These techniques have previously been used to characterize S-DNA-an elongated DNA conformation that forms when DNA overstretches at forces ≥ 65 pN. In this way, it was deduced that the base pairs of S-DNA are highly inclined, relative to those in relaxed (B-form) DNA. However, it is unclear whether and how topological constraints on the DNA may influence the base-pair inclinations under tension. Here, we apply polarization microscopy to investigate the impact of DNA pulling geometry, torsional constraint, and negative supercoiling on the orientations of intercalated dyes during overstretching. In contrast to earlier predictions, the pulling geometry (namely, whether the DNA molecule is stretched via opposite strands or the same strand) is found to have little influence. However, torsional constraint leads to a substantial reduction in intercalator tilting in overstretched DNA, particularly in AT-rich sequences. Surprisingly, the extent of intercalator tilting is similarly reduced when the DNA molecule is negatively supercoiled up to a critical supercoiling density (corresponding to ∼70% reduction in the linking number). We attribute these observations to the presence of P-DNA (an overwound DNA conformation). Our results suggest that intercalated DNA preferentially flanks regions of P-DNA rather than those of S-DNA and also substantiate previous suggestions that P-DNA forms predominantly in AT-rich sequences.

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.

Understanding enhanced mechanical stability of DNA in the presence of intercalated anticancer drug: Implications for DNA associated processes

Most of the anticancer drugs bind to double-stranded DNA (dsDNA) by intercalative-binding mode. Although experimental studies have become available recently, a molecular-level understanding of the interactions between the drug and dsDNA that lead to the stability of the intercalated drug is lacking. Of particular interest are the modifications of the mechanical properties of dsDNA observed in experiments. The latter could affect many biological functions, such as DNA transcription and replication. Here, we probe, via all-atom molecular dynamics (MD) simulations, the change in the mechanical properties of intercalated drug-DNA complexes for two intercalators, daunomycin and ethidium. We find that, upon drug intercalation, the stretch modulus of DNA increases significantly, whereas its persistence length and bending modulus decrease. Steered MD simulations reveal that it requires higher forces to stretch the intercalated dsDNA complexes than the normal dsDNA. Adopting various pulling protocols to study force-induced DNA melting, we find that the dissociation of dsDNA becomes difficult in the presence of intercalators. The results obtained here provide a plausible mechanism of function of the anticancer drugs, i.e., via altering the mechanical properties of DNA. We also discuss long-time consequences of using these drugs, which require further in vivo investigations.

A stretched conformation of DNA with a biological role?

We have discovered a well-defined extended conformation of double-stranded DNA, which we call Σ-DNA, using laser-tweezers force-spectroscopy experiments. At a transition force corresponding to free energy change ΔG = 1·57 ± 0·12 kcal (mol base pair)-1 60 or 122 base-pair long synthetic GC-rich sequences, when pulled by the 3'-3' strands, undergo a sharp transition to the 1·52 ± 0·04 times longer Σ-DNA. Intriguingly, the same degree of extension is also found in DNA complexes with recombinase proteins, such as bacterial RecA and eukaryotic Rad51. Despite vital importance to all biological organisms for survival, genome maintenance and evolution, the recombination reaction is not yet understood at atomic level. We here propose that the structural distortion represented by Σ-DNA, which is thus physically inherent to the nucleic acid, is related to how recombination proteins mediate recognition of sequence homology and execute strand exchange. Our hypothesis is that a homogeneously stretched DNA undergoes a 'disproportionation' into an inhomogeneous Σ-form consisting of triplets of locally B-like perpendicularly stacked bases. This structure may ensure improved fidelity of base-pair recognition and promote rejection in case of mismatch during homologous recombination reaction. Because a triplet is the length of a gene codon, we speculate that the structural physics of nucleic acids may have biased the evolution of recombinase proteins to exploit triplet base stacks and also the genetic code.

DNA partitions into triplets under tension in the presence of organic cations, with sequence evolutionary age predicting the stability of the triplet phase

Using atomistic simulations, we show the formation of stable triplet structure when particular GC-rich DNA duplexes are extended in solution over a timescale of hundreds of nanoseconds, in the presence of organic salt. We present planar-stacked triplet disproportionated DNA (Σ DNA) as a possible solution phase of the double helix under tension, subject to sequence and the presence of stabilising co-factors. Considering the partitioning of the duplexes into triplets of base pairs as the first step of operation of recombinase enzymes like RecA, we emphasise the structure-function relationship in Σ DNA. We supplement atomistic calculations with thermodynamic arguments to show that codons for 'phase 1' amino acids (those appearing early in evolution) are more likely than a lower entropy GC-rich sequence to form triplets under tension. We further observe that the four amino acids supposed (in the 'GADV world' hypothesis) to constitute the minimal set to produce functional globular proteins have the strongest triplet-forming propensity within the phase 1 set, showing a series of decreasing triplet propensity with evolutionary newness. The weak form of our observation provides a physical mechanism to minimise read frame and recombination alignment errors in the early evolution of the genetic code.

Remarkable similarity of force induced dsRNA conformational changes to stretched dsDNA and their detection using electrical measurements

We show the emergence of S-RNA under stretching in analogy to S-DNA and propose a method for its detection using electrical measurement.

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.

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.

Rupture of DNA aptamer: New insights from simulations

Base-pockets (non-complementary base-pairs) in a double-stranded DNA play a crucial role in biological processes. Because of thermal fluctuations, it can lower the stability of DNA, whereas, in case of DNA aptamer, small molecules, e.g., adenosinemonophosphate and adenosinetriphosphate, form additional hydrogen bonds with base-pockets termed as "binding-pockets," which enhance the stability. Using the Langevin dynamics simulations of coarse grained model of DNA followed by atomistic simulations, we investigated the influence of base-pocket and binding-pocket on the stability of DNA aptamer. Striking differences have been reported here for the separation induced by temperature and force, which require further investigation by single molecule experiments.

Dramatic changes in DNA conductance with stretching: structural polymorphism at a critical extension

In order to interpret recent experimental studies of the dependence of conductance of ds-DNA as the DNA is pulled from the 3'end1-3'end2 ends, which find a sharp conductance jump for a very short (4.5%) stretching length, we carried out multiscale modeling to predict the conductance of dsDNA as it is mechanically stretched to promote various structural polymorphisms. We calculate the current along the stretched DNA using a combination of molecular dynamics simulations, non-equilibrium pulling simulations, quantum mechanics calculations, and kinetic Monte Carlo simulations. For 5'end1-5'end2 attachments we find an abrupt jump in the current within a very short stretching length (6 Å or 17%) leading to a melted DNA state. In contrast, for 3'end1-3'end2 pulling it takes almost 32 Å (84%) of stretching to cause a similar jump in the current. Thus, we demonstrate that charge transport in DNA can occur over stretching lengths of several nanometers. We find that this unexpected behaviour in the B to S conformational DNA transition arises from highly inclined base pair geometries that result from this pulling protocol. We found that the dramatically different conductance behaviors for two different pulling protocols arise from how the hydrogen bonds of DNA base pairs break.

Single molecule identification of homology-dependent interactions between long ssRNA and dsDNA

Long non-coding RNAs (lncRNAs) are prominently associated with chromosomes in an ever-increasing diversity of roles. To provide further insight into the potential nature of these associations, we have explored, for the first time, the interaction of long single-stranded (ss) RNAs with cognate homologous double-stranded (ds) DNA in vitro. Using magnetic tweezers, we measured the effects of ssRNA on force extension curves for dsDNA. We observe that the presence of ssRNA impedes the extension of dsDNA, specifically at low forces, dependent on homology between the RNA and DNA species, and dependent on ssRNA lengths (≥1 kb). The observed effect also depends on the concentration of ssRNA and is abolished by overstretching of the dsDNA. These findings show that significant homologous contacts can occur between long ssRNA and dsDNA in the absence of protein and that these contacts alter the mechanical properties of the dsDNA. We propose that long ssRNA interacts paranemically with long dsDNA via periodic short homologous interactions, e.g. mediated by RNA/DNA triplex-formation, and that dsDNA extension is impeded by formation of RNA secondary structure in the intervening unbound regions. Analogous interactions in vivo would permit lncRNAs to mediate the juxtaposition of two or more DNA regions on the same or different chromosomes.

Single-molecule mechanics of protein-labelled DNA handles

DNA handles are often used as spacers and linkers in single-molecule experiments to isolate and tether RNAs, proteins, enzymes and ribozymes, amongst other biomolecules, between surface-modified beads for nanomechanical investigations. Custom DNA handles with varying lengths and chemical end-modifications are readily and reliably synthesized en masse, enabling force spectroscopic measurements with well-defined and long-lasting mechanical characteristics under physiological conditions over a large range of applied forces. Although these chemically tagged DNA handles are widely used, their further individual modification with protein receptors is less common and would allow for additional flexibility in grabbing biomolecules for mechanical measurements. In-depth information on reliable protocols for the synthesis of these DNA-protein hybrids and on their mechanical characteristics under varying physiological conditions are lacking in literature. Here, optical tweezers are used to investigate different protein-labelled DNA handles in a microfluidic environment under different physiological conditions. Digoxigenin (DIG)-dsDNA-biotin handles of varying sizes (1000, 3034 and 4056 bp) were conjugated with streptavidin or neutravidin proteins. The DIG-modified ends of these hybrids were bound to surface-modified polystyrene (anti-DIG) beads. Using different physiological buffers, optical force measurements showed consistent mechanical characteristics with long dissociation times. These protein-modified DNA hybrids were also interconnected in situ with other tethered biotinylated DNA molecules. Electron-multiplying CCD (EMCCD) imaging control experiments revealed that quantum dot-streptavidin conjugates at the end of DNA handles remain freely accessible. The experiments presented here demonstrate that handles produced with our protein-DNA labelling procedure are excellent candidates for grasping single molecules exposing tags suitable for molecular recognition in time-critical molecular motor studies.

The poor homology stringency in the heteroduplex allows strand exchange to incorporate desirable mismatches without sacrificing recognition in vivo

RecA family proteins are responsible for homology search and strand exchange. In bacteria, homology search begins after RecA binds an initiating single-stranded DNA (ssDNA) in the primary DNA-binding site, forming the presynaptic filament. Once the filament is formed, it interrogates double-stranded DNA (dsDNA). During the interrogation, bases in the dsDNA attempt to form Watson–Crick bonds with the corresponding bases in the initiating strand. Mismatch dependent instability in the base pairing in the heteroduplex strand exchange product could provide stringent recognition; however, we present experimental and theoretical results suggesting that the heteroduplex stability is insensitive to mismatches. We also present data suggesting that an initial homology test of 8 contiguous bases rejects most interactions containing more than 1/8 mismatches without forming a detectable 20 bp product. We propose that, in vivo, the sparsity of accidental sequence matches allows an initial 8 bp test to rapidly reject almost all non-homologous sequences. We speculate that once the initial test is passed, the mismatch insensitive binding in the heteroduplex allows short mismatched regions to be incorporated in otherwise homologous strand exchange products even though sequences with less homology are eventually rejected.

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.

Force-induced melting of DNA--evidence for peeling and internal melting from force spectra on short synthetic duplex sequences

Overstretching of DNA occurs at about 60-70 pN when a torsionally unconstrained double-stranded DNA molecule is stretched by its ends. During the transition, the contour length increases by up to 70% without complete strand dissociation. Three mechanisms are thought to be involved: force-induced melting into single-stranded DNA where either one or both strands carry the tension, or a B-to-S transition into a longer, still base-paired conformation. We stretch sequence-designed oligonucleotides in an effort to isolate the three processes, focusing on force-induced melting. By introducing site-specific inter-strand cross-links in one or both ends of a 64 bp AT-rich duplex we could repeatedly follow the two melting processes at 5 mM and 1 M monovalent salt. We find that when one end is sealed the AT-rich sequence undergoes peeling exhibiting hysteresis at low and high salt. When both ends are sealed the AT sequence instead undergoes internal melting. Thirdly, the peeling melting is studied in a composite oligonucleotide where the same AT-rich sequence is concatenated to a GC-rich sequence known to undergo a B-to-S transition rather than melting. The construct then first melts in the AT-rich part followed at higher forces by a B-to-S transition in the GC-part, indicating that DNA overstretching modes are additive.

Magnetic force micropiston: an integrated force/microfluidic device for the application of compressive forces in a confined environment

Cellular biology takes place inside confining spaces. For example, bacteria grow in crevices, red blood cells squeeze through capillaries, and chromosomes replicate inside the nucleus. Frequently, the extent of this confinement varies. Bacteria grow longer and divide, red blood cells move through smaller and smaller passages as they travel to capillary beds, and replication doubles the amount of DNA inside the nucleus. This increase in confinement, either due to a decrease in the available space or an increase in the amount of material contained in a constant volume, has the potential to squeeze and stress objects in ways that may lead to changes in morphology, dynamics, and ultimately biological function. Here, we describe a device developed to probe the interplay between confinement and the mechanical properties of cells and cellular structures, and forces that arise due to changes in a structure's state. In this system, the manipulation of a magnetic bead exerts a compressive force upon a target contained in the confining space of a microfluidic channel. This magnetic force microfluidic piston is constructed in such a way that we can measure (a) target compliance and changes in compliance as induced by changes in buffer, extract, or biochemical composition, (b) target expansion force generated by changes in the same parameters, and (c) the effects of compression stress on a target's structure and function. Beyond these issues, our system has general applicability to a variety of questions requiring the combination of mechanical forces, confinement, and optical imaging.

An effective mesoscopic model of double-stranded DNA

Watson and Crick's epochal presentation of the double helix structure in 1953 has paved the way to intense exploration of DNA's vital functions in cells. Also, recent advances of single molecule techniques have made it possible to probe structures and mechanics of constrained DNA at length scales ranging from nanometers to microns. There have been a number of atomistic scale quantum chemical calculations or molecular level simulations, but they are too computationally demanding or analytically unfeasible to describe the DNA conformation and mechanics at mesoscopic levels. At micron scales, on the other hand, the wormlike chain model has been very instrumental in describing analytically the DNA mechanics but lacks certain molecular details that are essential in describing the hybridization, nano-scale confinement, and local denaturation. To fill this fundamental gap, we present a workable and predictive mesoscopic model of double-stranded DNA where the nucleotides beads constitute the basic degrees of freedom. With the inter-strand stacking given by an interaction between diagonally opposed monomers, the model explains with analytical simplicity the helix formation and produces a generalized wormlike chain model with the concomitant large bending modulus given in terms of the helical structure and stiffness. It also explains how the helical conformation undergoes overstretch transition to the ladder-like conformation at a force plateau, in agreement with the experiment.

Coarse-grained simulations of DNA overstretching

We use a recently developed coarse-grained model to simulate the overstretching of duplex DNA. Overstretching at 23 °C occurs at 74 pN in the model, about 6-7 pN higher than the experimental value at equivalent salt conditions. Furthermore, the model reproduces the temperature dependence of the overstretching force well. The mechanism of overstretching is always force-induced melting by unpeeling from the free ends. That we never see S-DNA (overstretched duplex DNA), even though there is clear experimental evidence for this mode of overstretching under certain conditions, suggests that S-DNA is not simply an unstacked but hydrogen-bonded duplex, but instead probably has a more exotic structure.

Effect of undulations on spontaneous braid formation

This paper is an extension of a recent study where it was shown that forces dependent on molecular helical structure may cause two DNA molecules to spontaneously braid [R. Cortini et al., Biophys. J. 101, 875 (2011)]. Here, bending fluctuations of molecular center lines about the braid axis are incorporated into the braiding theory, which may be generalized to other helix-dependent interactions and other helical molecules. The free energy of the pair of molecules is recalculated and compared to its value without incorporating undulations. We find that the loss of configurational entropy due to confinement of the molecules in the braid is quite high. This contribution to the free energy increases the amount of attraction needed for spontaneous braiding due to helix-dependent forces. The theory will be further developed for plectonemes and braids under mechanical forces in later work.

Tension induces a base-paired overstretched DNA conformation

Mixed-sequence DNA molecules undergo mechanical overstretching by approximately 70% at 60-70 pN. Since its initial discovery 15 y ago, a debate has arisen as to whether the molecule adopts a new form [Cluzel P, et al. (1996) Science 271:792-794; Smith SB, Cui Y, Bustamante C (1996) Science 271:795-799], or simply denatures under tension [van Mameren J, et al. (2009) Proc Natl Acad Sci USA 106:18231-18236]. Here, we resolve this controversy by using optical tweezers to extend small 60-64 bp single DNA duplex molecules whose base content can be designed at will. We show that when AT content is high (70%), a force-induced denaturation of the DNA helix ensues at 62 pN that is accompanied by an extension of the molecule of approximately 70%. By contrast, GC-rich sequences (60% GC) are found to undergo a reversible overstretching transition into a distinct form that is characterized by a 51% extension and that remains base-paired. For the first time, results proving the existence of a stretched basepaired form of DNA can be presented. The extension observed in the reversible transition coincides with that produced on DNA by binding of bacterial RecA and human Rad51, pointing to its possible relevance in homologous recombination.

Dynamical model for the full stretching curve of DNA

We present a phenomenological dynamical model able to describe the stretching features of the curve of DNA length vs applied force. As concerns the chain, the model is based on the discrete wormlike chain model with elastic modifications, which properly describes the elongation features at low and intermediate forces. The dynamics is developed under a double-well potential with a linear term, which, at high forces, accounts for the narrow transition present in the DNA elongation (overstretching). A quite good agreement between simulation and experiment is obtained.

Complementary strand relocation may play vital roles in RecA-based homology recognition

RecA-family proteins mediate homologous recombination and recombinational DNA repair through homology search and strand exchange. Initially, the protein forms a filament with the incoming single-stranded DNA (ssDNA) bound in site I. The RecA–ssDNA filament then binds double-stranded DNA (dsDNA) in site II. Non-homologous dsDNA rapidly unbinds, whereas homologous dsDNA undergoes strand exchange yielding heteroduplex dsDNA in site I and the leftover outgoing strand in site II. We show that applying force to the ends of the complementary strand significantly retards strand exchange, whereas applying the same force to the outgoing strand does not. We also show that crystallographically determined binding site locations require an intermediate structure in addition to the initial and final structures. Furthermore, we demonstrate that the characteristic dsDNA extension rates due to strand exchange and free RecA binding are the same, suggesting that relocation of the complementary strand from its position in the intermediate structure to its position in the final structure limits both rates. Finally, we propose that homology recognition is governed by transitions to and from the intermediate structure, where the transitions depend on differential extension in the dsDNA. This differential extension drives strand exchange forward for homologs and increases the free energy penalty for strand exchange of non-homologs.

Two distinct overstretched DNA structures revealed by single-molecule thermodynamics measurements

Double-stranded DNA is a dynamic molecule whose structure can change depending on conditions. While there is consensus in the literature about many structures DNA can have, the state of highly-stretched DNA is still not clear. Several groups have shown that DNA in the torsion-unconstrained B-form undergoes an "overstretching" transition at a stretching force of around 65 pN, which leads to approximately 1.7-fold elongation of the DNA contour length. Recent experiments have revealed that two distinct structural transitions are involved in the overstretching process: (i) a hysteretic "peeling" off one strand from its complementary strand, and (ii) a nonhysteretic transition that leads to an undetermined DNA structure. We report the first simultaneous determination of the entropy (ΔS) and enthalpy changes (ΔH) pertaining to these respective transitions. For the hysteretic peeling transition, we determined ΔS ∼ 20 cal/(K.mol) and ΔH ∼ 7 kcal/mol. In the case of the nonhysteretic transition, ΔS ∼ -3 cal/(K.mol) and ΔH ∼ 1 kcal/mol. Furthermore, the response of the transition force to salt concentration implies that the two DNA strands are spatially separated after the hysteretic peeling transition. In contrast, the corresponding response after the nonhysteretic transition indicated that the strands remained in close proximity. The selection between the two transitions depends on DNA base-pair stability, and it can be illustrated by a multidimensional phase diagram. Our results provide important insights into the thermodynamics of DNA overstretching and conformational structures of overstretched DNA that may play an important role in vivo.

PicoNewton-millisecond force steps reveal the transition kinetics and mechanism of the double-stranded DNA elongation

We study the kinetics of the overstretching transition in λ-phage double-stranded (ds) DNA from the basic conformation (B state) to the 1.7-times longer and partially unwound conformation (S state), using the dual-laser optical tweezers under force-clamp conditions at 25°C. The unprecedented resolution of our piezo servo-system, which can impose millisecond force steps of 0.5-2 pN, reveals the exponential character of the elongation kinetics and allows us to test the two-state nature of the B-S transition mechanism. By analyzing the load-dependence of the rate constant of the elongation, we find that the elementary elongation step is 5.85 nm, indicating a cooperativity of ~25 basepairs. This mechanism increases the free energy for the elementary reaction to ~94 k(B)T, accounting for the stability of the basic conformation of DNA, and explains why ds-DNA can remain in equilibrium as it overstretches.

RecA homology search is promoted by mechanical stress along the scanned duplex DNA

A RecA–single-stranded DNA (RecA–ssDNA) filament searches a genome for sequence homology by rapidly binding and unbinding double-stranded DNA (dsDNA) until homology is found. We demonstrate that pulling on the opposite termini (3′ and 5′) of one of the two DNA strands in a dsDNA molecule stabilizes the normally unstable binding of that dsDNA to non-homologous RecA–ssDNA filaments, whereas pulling on the two 3′, the two 5′, or all four termini does not. We propose that the ‘outgoing’ strand in the dsDNA is extended by strong DNA–protein contacts, whereas the ‘complementary’ strand is extended by the tension on the base pairs that connect the ‘complementary’ strand to the ‘outgoing’ strand. The stress resulting from different levels of tension on its constitutive strands causes rapid dsDNA unbinding unless sufficient homology is present.

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.

B-DNA to zip-DNA: simulating a DNA transition to a novel structure with enhanced charge-transport characteristics

The forced extension of a DNA segment is studied in a series of steered molecular dynamics simulations, employing a broad range of pulling forces. Throughout the entire force range, the formation of a zipper-like (zip-) DNA structure is observed. In that structure, first predicted by Lohikoski et al., the bases of the DNA strands interdigitate with each other and form a single-base aromatic stack. Similar motifs, albeit only a few base pairs in extent, have been observed in experimental crystal structures. Analysis of the dynamics of structural changes in pulled DNA shows that S-form DNA, thought to be adopted by DNA under applied force, serves as an intermediate between B-DNA and zip-DNA. Therefore, the phase transition plateau observed in force-extension curves of DNA is suggested to reflect the B-DNA to zip-DNA structural transition. Electronic structure analysis of purine bases in zip-DNA indicates a several-fold to order of magnitude increase in the π-π electronic coupling among nearest-neighbor nucleobases, compared to B-DNA. We further observe that zip-DNA does not require base pair complementarity between DNA strands, and we predict that the increased electronic coupling in zip-DNA will result in a much higher rate of charge transfer through an all-purine zip-DNA compared to B-DNA of equal length.

Changes in the tension in dsDNA alter the conformation of RecA bound to dsDNA-RecA filaments

The RecA protein is an ATPase that mediates recombination via strand exchange. In strand exchange a single-stranded DNA (ssDNA) bound to RecA binding site I in a RecA/ssDNA filament pairs with one strand of a double-stranded DNA (dsDNA) and forms heteroduplex dsDNA in site I if homology is encountered. Long sequences are exchanged in a dynamic process in which initially unbound dsDNA binds to the leading end of a RecA/ssDNA filament, while heteroduplex dsDNA unbinds from the lagging end via ATP hydrolysis. ATP hydrolysis is required to convert the active RecA conformation, which cannot unbind, to the inactive conformation, which can unbind. If dsDNA extension due to RecA binding increases the dsDNA tension, then RecA unbinding must decrease tension. We show that in the presence of ATP hydrolysis decreases in tension induce decreases in length whereas in the absence of hydrolysis, changes in tension have no systematic effect. These results suggest that decreases in force enhance dissociation by promoting transitions from the active to the inactive RecA conformation. In contrast, increases in tension reduce dissociation. Thus, the changes in tension inherent to strand exchange may couple with ATP hydrolysis to increase the directionality and stringency of strand exchange.

Improved high-force magnetic tweezers for stretching and refolding of proteins and short DNA

Although magnetic tweezers have many unique advantages in terms of specificity, throughput, and force stability, this tool has had limited application on short tethers because accurate measurement of force has been difficult for short tethers under large tension. Here, we report a method that allows us to apply magnetic tweezers to stretch short biomolecules with accurate force calibration over a wide range of up to 100 pN. We demonstrate the use of the method by overstretching of a short DNA and unfolding/refolding a protein of filamin A immunoglobulin domains 1-8. Other potential applications of this method are also discussed.

Single-molecule studies of the stringency factors and rates governing the polymerization of RecA on double-stranded DNA

RecA is a key protein in homologous recombination. During recombination, one single-stranded DNA (ssDNA) bound to site I in RecA exchanges Watson-Crick pairing with a sequence-matched ssDNA that was part of a double-stranded DNA molecule (dsDNA) bound to site II in RecA. After strand exchange, heteroduplex dsDNA is bound to site I. In vivo, direct polymerization of RecA on dsDNA through site I does not occur, though it does in vitro. The mechanisms underlying the difference have been unclear. We use single-molecule experiments to decouple the two steps involved in polymerization: nucleation and elongation. We find that elongation is governed by a fundamental clock that is insensitive to force and RecA concentration from 0.2 and 6 µM, though rates depend on ionic conditions. Thus, we can probe nucleation site stability by creating nucleation sites at high force and then measuring elongation as a function of applied force. We find that in the presence of ATP hydrolysis a minimum force is required for polymerization. The minimum force decreases with increasing RecA or ATP concentrations. We propose that force reduces the off-rate for nucleation site binding and that nucleation site stability is the stringency factor that prevents in vivo polymerization.

Transition dynamics and selection of the distinct S-DNA and strand unpeeling modes of double helix overstretching

Recent studies have revealed two distinct pathways for the DNA overstretching transition near 65 pN: ‘unpeeling’ of one strand from the other, and a transition from B-DNA to an elongated double-stranded ‘S-DNA’ form. However, basic questions concerning the dynamics of these transitions, relative stability of the two competing overstretched states, and effects of nicks and free DNA ends on overstretching, remain open. In this study we report that: (i) stepwise extension changes caused by sequence-defined barriers occur during the strand-unpeeling transition, whereas rapid, sequence-independent extension fluctuations occur during the B to S transition; (ii) the secondary transition that often occurs following the overstretching transition is strand-unpeeling, during which the extension increases by 0.01–0.02 nm per base pair of S-DNA converted to single-stranded DNA at forces between 75 and 110 pN; (iii) even in the presence of nicks or free ends, S-DNA can be stable under physiological solution conditions; (iv) distribution of small GC-rich islands in a large DNA plays a key role in determining the transition pathways; and (v) in the absence of nicks or free ends, torsion-unconstrained DNA undergoes the overstretching transition via creation of S-DNA. Our study provides a new, high-resolution understanding of the competition between unpeeling and formation of S-DNA.

Direct observation of multiple pathways of single-stranded DNA stretching

We observed multiple pathways of stretching single-stranded polydeoxynucleotides, poly(dA). Poly(dA) has been shown to undergo unique transitions under mechanical force, and such transitions were attributed to the stacking characteristics of poly(dA). Using single-molecule manipulation studies, we found that poly(dA) has two stretching pathways at high forces. The previously observed pathway has a free energy that is less than what is expected of single-stranded DNA with a random sequence, indicating the existence of a novel conformation of poly(dA) at large extensions. We also observed stepwise transitions between the two pathways by pulling the molecule with constant force, and found that the transitions are cooperative. These results suggest that the unique mechanical property of poly(dA) may play an important role in biological processes such as gene expression.

Quantitative modeling of forces in electromagnetic tweezers

This paper discusses numerical simulations of the magnetic field produced by an electromagnet for generation of forces on superparamagnetic microspheres used in manipulation of single molecules or cells. Single molecule force spectroscopy based on magnetic tweezers can be used in applications that require parallel readout of biopolymer stretching or biomolecular binding. The magnetic tweezers exert forces on the surface-immobilized macromolecule by pulling a magnetic bead attached to the free end of the molecule in the direction of the field gradient. In a typical force spectroscopy experiment, the pulling forces can range between subpiconewton to tens of piconewtons. In order to effectively provide such forces, an understanding of the source of the magnetic field is required as the first step in the design of force spectroscopy systems. In this study, we use a numerical technique, the method of auxiliary sources, to investigate the influence of electromagnet geometry and material parameters of the magnetic core on the magnetic forces pulling the target beads in the area of interest. The close proximity of the area of interest to the magnet body results in deviations from intuitive relations between magnet size and pulling force, as well as in the force decay with distance. We discuss the benefits and drawbacks of various geometric modifications affecting the magnitude and spatial distribution of forces achievable with an electromagnet.

Single-molecule stretching studies of RNA chaperones

RNA chaperone proteins play significant roles in diverse biological contexts. The most widely studied RNA chaperones are the retroviral nucleocapsid proteins (NC), also referred to as nucleic acid (NA) chaperones. Surprisingly, the biophysical properties of the NC proteins vary significantly for different viruses, and it appears that HIV-1 NC has optimal NA chaperone activity. In this review we discuss the physical nature of the NA chaperone activity of NC. We conclude that the optimal NA chaperone must saturate NA binding, leading to strong NA aggregation and slight destabilization of all NA duplexes. Finally, rapid kinetics of the chaperone protein interaction with NA is another primary component of its NA chaperone activity. We discuss these characteristics of HIV-1 NC and compare them with those of other NA binding proteins and ligands that exhibit only some characteristics of NA chaperone activity, as studied by single molecule DNA stretching.

Study of force induced melting of dsDNA as a function of length and conformation

We measure the constant force required to melt double-stranded (ds) DNA as a function of length for lengths from 12 to 100,000 base pairs, where the force is applied to the 3'3' or 5'5' ends of the dsDNA. Molecules with 32 base pairs or fewer melt before overstretching. For these short molecules, the melting force is independent of the ends to which the force is applied and the shear force as a function of length is well described by de Gennes theory with a de Gennes length of less than 10 bp. Molecules with lengths of 500 base pairs or more overstretch before melting. For these long molecules, the melting force depends on the ends to which the force is applied. The melting force as a function of length increases even when the length exceeds 1000 bp, where the length dependence is inconsistent with de Gennes theory. Finally, we expand de Gennes melting theory to 3'5' pulling and compare the predictions with experimental results.

DNA structural changes under different stretching methods studied by molecular dynamics simulations

We present a molecular dynamics simulation study of 22-mer DNA conformational variations obtained by stretching both 3'-termini and both 5'-termini. Stretching 3'-termini by 3.5 nm required 142 kJ mol(-1) and the force plateau was approximately 80 pN, whereas stretching 5'-termini by the same length required 190 kJ mol(-1) and the force plateau was approximately 100 pN. Stretching 3'-termini led to a larger untwisting of the double helix and the successive base pairs rolled to the side of the DNA minor groove, while stretching 5'-termini resulted in the base pairs rolling to the major groove side and reducing of the diameter of DNA molecule. The most distinctive difference between stretching 3'-termini and 5'-termini was that at the force plateau region stretching the 5'-termini resulted in breakage of the base pairs, which considerably disturbed the structure of the DNA double helix. All of the variations of base rotation and translation for both stretching methods took place when the relative length of DNA l was longer than 1.2, which was the point the force plateau appeared.

Undulations enhance the effect of helical structure on DNA interactions

During the past decade, theory and experiments have provided clear evidence that specific helical patterns of charged groups and adsorbed (condensed) counterions on the DNA surface are responsible for many important features of DNA-DNA interactions in hydrated aggregates. The effects of helical structure on DNA-DNA interactions result from a preferential juxtaposition of the negatively charged sugar phosphate backbone with counterions bound within the grooves of the opposing molecule. Analysis of X-ray diffraction experiments confirmed the mutual alignment of parallel molecules in hydrated aggregates required for such juxtaposition. However, it remained unclear how this alignment and molecular interactions might be affected by intrinsic and thermal fluctuations, which cause structural deviations away from an ideal double helical conformation. We previously argued that the torsional flexibility of DNA allows the molecules to adapt their structure to accommodate a more electrostatically favorable alignment between molecules, partially compensating disruptive fluctuation effects. In the present work, we develop a more comprehensive theory, incorporating also stretching and bending fluctuations of DNA. We found the effects of stretching to be qualitatively and quantitatively similar to those of twisting fluctuations. However, this theory predicts more dramatic and surprising effects of bending. Undulations of DNA in hydrated aggregates strongly amplify rather than weaken the helical structure effects. They enhance the structural adaptation, leading to better alignment of neighboring molecules and pushing the geometry of the DNA backbone closer to that of an ideal helix. These predictions are supported by a quantitative comparison of the calculated and measured osmotic pressures in DNA.