DNA-DNA interactions in tight supercoils are described by a small effective charge density

Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations

Inferring the 3D structure and conformation of disordered biomolecules, e.g., single stranded nucleic acids (ssNAs), remains challenging due to their conformational heterogeneity in solution. Here, we use escape-time electrometry (ETe) to measure with sub elementary-charge precision the effective electrical charge in solution of short to medium chain length ssNAs in the range of 5-60 bases. We compare measurements of molecular effective charge with theoretically calculated values for simulated molecular conformations obtained from Molecular Dynamics simulations using a variety of forcefield descriptions. We demonstrate that the measured effective charge captures subtle differences in molecular structure in various nucleic acid homopolymers of identical length, and also that the experimental measurements can find agreement with computed values derived from coarse-grained molecular structure descriptions such as oxDNA, as well next generation ssNA force fields. We further show that comparing the measured effective charge with calculations for a rigid, charged rod-the simplest model of a nucleic acid-yields estimates of molecular structural dimensions such as linear charge spacings that capture molecular structural trends observed using high resolution structural analysis methods such as X-ray scattering. By sensitively probing the effective charge of a molecule, electrometry provides a powerful dimension supporting inferences of molecular structural and conformational properties, as well as the validation of biomolecular structural models. The overall approach holds promise for a high throughput, microscopy-based biomolecular analytical approach offering rapid screening and inference of molecular 3D conformation, and operating at the single molecule level in solution.

Nucleosome spacing controls chromatin spatial structure and accessibility

Recent research highlights the significance of the three-dimensional structure of chromatin in regulating various cellular processes, particularly transcription. This is achieved through dynamic chromatin structures that facilitate long-range contacts and control spatial accessibility. Chromatin consists of DNA and a variety of proteins, of which histones play an essential structural role by forming nucleosomes. Extensive experimental and theoretical research in recent decades has yielded conflicting results about key factors that regulate the spatial structure of chromatin, which remains enigmatic. By using a computer model that allows us to simulate chromatin volumes containing physiological nucleosome concentrations, we investigated whether nucleosome spacing or nucleosome density is fundamental for three-dimensional chromatin accessibility. Unexpectedly, the regularity of the nucleosome spacing is crucial for determining the accessibility of the chromatin network to diffusive processes, whereas variation in nucleosome concentrations has only minor effects. Using only the basic physical properties of DNA and nucleosomes was sufficient to generate chromatin structures consistent with published electron microscopy data. Contrary to other work, we found that nucleosome density did not substantially alter the properties of chromatin fibers or contact probabilities of genomic loci. No breakup of fiber-like structures was observed at high molar density. These findings challenge previous assumptions and highlight the importance of nucleosome spacing as a key driver of chromatin organization. These results identified changes in nucleosome spacing as a tentative mechanism for altering the spatial chromatin structure and thus genomic functions.

Dynamic interplay between target search and recognition for a Type I CRISPR-Cas system

CRISPR-Cas effector complexes enable the defense against foreign nucleic acids and have recently been exploited as molecular tools for precise genome editing at a target locus. To bind and cleave their target, the CRISPR-Cas effectors have to interrogate the entire genome for the presence of a matching sequence. Here we dissect the target search and recognition process of the Type I CRISPR-Cas complex Cascade by simultaneously monitoring DNA binding and R-loop formation by the complex. We directly quantify the effect of DNA supercoiling on the target recognition probability and demonstrate that Cascade uses facilitated diffusion for its target search. We show that target search and target recognition are tightly linked and that DNA supercoiling and limited 1D diffusion need to be considered when understanding target recognition and target search by CRISPR-Cas enzymes and engineering more efficient and precise variants.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

A quantitative model for the dynamics of target recognition and off-target rejection by the CRISPR-Cas Cascade complex

CRISPR-Cas effector complexes recognise nucleic acid targets by base pairing with their crRNA which enables easy re-programming of the target specificity in rapidly emerging genome engineering applications. However, undesired recognition of off-targets, that are only partially complementary to the crRNA, occurs frequently and represents a severe limitation of the technique. Off-targeting lacks comprehensive quantitative understanding and prediction. Here, we present a detailed analysis of the target recognition dynamics by the Cascade surveillance complex on a set of mismatched DNA targets using single-molecule supercoiling experiments. We demonstrate that the observed dynamics can be quantitatively modelled as a random walk over the length of the crRNA-DNA hybrid using a minimal set of parameters. The model accurately describes the recognition of targets with single and double mutations providing an important basis for quantitative off-target predictions. Importantly the model intrinsically accounts for observed bias regarding the position and the proximity between mutations and reveals that the seed length for the initiation of target recognition is controlled by DNA supercoiling rather than the Cascade structure.

DNA fluctuations reveal the size and dynamics of topological domains

DNA supercoiling is a key regulatory mechanism that orchestrates DNA readout, recombination, and genome maintenance. DNA-binding proteins often mediate these processes by bringing two distant DNA sites together, thereby inducing (transient) topological domains. In order to understand the dynamics and molecular architecture of protein-induced topological domains in DNA, quantitative and time-resolved approaches are required. Here, we present a methodology to determine the size and dynamics of topological domains in supercoiled DNA in real time and at the single-molecule level. Our approach is based on quantifying the extension fluctuations-in addition to the mean extension-of supercoiled DNA in magnetic tweezers (MT). Using a combination of high-speed MT experiments, Monte Carlo simulations, and analytical theory, we map out the dependence of DNA extension fluctuations as a function of supercoiling density and external force. We find that in the plectonemic regime, the extension variance increases linearly with increasing supercoiling density and show how this enables us to determine the formation and size of topological domains. In addition, we demonstrate how the transient (partial) dissociation of DNA-bridging proteins results in the dynamic sampling of different topological states, which allows us to deduce the torsional stiffness of the plectonemic state and the kinetics of protein-plectoneme interactions. We expect our results to further the understanding and optimization of magnetic tweezer measurements and to enable quantification of the dynamics and reaction pathways of DNA processing enzymes in the context of physiologically relevant forces and supercoiling densities.

Opto-Electrostatic Determination of Nucleic Acid Double-Helix Dimensions and the Structure of the Molecule–Solvent Interface

A DNA molecule is highly electrically charged in solution. The electrical potential at the molecular surface is known to vary strongly with the local geometry of the double helix and plays a pivotal role in DNA–protein interactions. Further out from the molecular surface, the electrical field propagating into the surrounding electrolyte bears fingerprints of the three-dimensional arrangement of the charged atoms in the molecule. However, precise extraction of the structural information encoded in the electrostatic “far field” has remained experimentally challenging. Here, we report an optical microscopy-based approach that detects the field distribution surrounding a charged molecule in solution, revealing geometric features such as the radius and the average rise per basepair of the double helix with up to sub-Angstrom precision, comparable with traditional molecular structure determination techniques like X-ray crystallography and nuclear magnetic resonance. Moreover, measurement of the helical radius furnishes an unprecedented view of both hydration and the arrangement of cations at the molecule–solvent interface. We demonstrate that a probe in the electrostatic far field delivers structural and chemical information on macromolecules, opening up a new dimension in the study of charged molecules and interfaces in solution.

Differences in nanoscale organization of regulatory active and inactive human chromatin

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.

Characterizing the free-energy landscapes of DNA origamis

We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.

Multi-resolution simulation of DNA transport through large synthetic nanostructures

Modeling and simulation has become an invaluable partner in development of nanopore sensing systems. The key advantage of the nanopore sensing method - the ability to rapidly detect individual biomolecules as a transient reduction of the ionic current flowing through the nanopore - is also its key deficiency, as the current signal itself rarely provides direct information about the chemical structure of the biomolecule. Complementing experimental calibration of the nanopore sensor readout, coarse-grained and all-atom molecular dynamics simulations have been used extensively to characterize the nanopore translocation process and to connect the microscopic events taking place inside the nanopore to the experimentally measured ionic current blockades. Traditional coarse-grained simulations, however, lack the precision needed to predict ionic current blockades with atomic resolution whereas traditional all-atom simulations are limited by the length and time scales amenable to the method. Here, we describe a multi-resolution framework for modeling electric field-driven passage of DNA molecules and nanostructures through to-scale models of synthetic nanopore systems. We illustrate the method by simulating translocation of double-stranded DNA through a solid-state nanopore and a micron-scale slit, capture and translocation of single-stranded DNA in a double nanopore system, and modeling ionic current readout from a DNA origami nanostructure passage through a nanocapillary. We expect our multi-resolution simulation framework to aid development of the nanopore field by providing accurate, to-scale modeling capability to research laboratories that do not have access to leadership supercomputer facilities.

Writhing and hockling instabilities in twisted elastic fibers

The buckling and twisting of slender, elastic fibers is a deep and well-studied field. A slender elastic rod that is twisted with respect to a fixed end will spontaneously form a loop, or hockle, to relieve the torsional stress that builds. Further twisting results in the formation of plectonemes-a helical excursion in the fiber that extends with additional twisting. Here we use an idealized, micron-scale experiment to investigate the energy stored, and subsequently released, by hockles and plectonemes as they are pulled apart, in analogy with force spectroscopy studies of DNA and protein folding. Hysteresis loops in the snapping and unsnapping inform the stored energy in the twisted fiber structures.

Probing the stability of the SpCas9-DNA complex after cleavage

CRISPR-Cas9 is a ribonucleoprotein complex that sequence-specifically binds and cleaves double-stranded DNA. Wildtype Cas9 and its nickase and cleavage-incompetent mutants have been used in various biological techniques due to their versatility and programmable specificity. Cas9 has been shown to bind very stably to DNA even after cleavage of the individual DNA strands, inhibiting further turnovers and considerably slowing down in-vivo repair processes. This poses an obstacle in genome editing applications. Here, we employed single-molecule magnetic tweezers to investigate the binding stability of different Streptococcus pyogenes Cas9 variants after cleavage by challenging them with supercoiling. We find that different release mechanisms occur depending on which DNA strand is cleaved. After initial target strand cleavage, supercoils are only removed after the collapse of the R-loop. We identified several states with different stabilities of the R-loop. Most importantly, we find that the post-cleavage state of Cas9 exhibits a higher stability than the pre-cleavage state. After non-target strand cleavage, supercoils are immediately but slowly released by swiveling of the non-target strand around Cas9 bound to the target strand. Consequently, Cas9 and its non-target strand nicking mutant stay stably bound to the DNA for many hours even at elevated torsional stress.

Nanomechanics of negatively supercoiled diaminopurine-substituted DNA

Single molecule experiments have demonstrated a progressive transition from a B- to an L-form helix as DNA is gently stretched and progressively unwound. The particular sequence of a DNA segment defines both base stacking and hydrogen bonding that affect the partitioning and conformations of the two phases. Naturally or artificially modified bases alter H-bonds and base stacking and DNA with diaminopurine (DAP) replacing adenine was synthesized to produce linear fragments with triply hydrogen-bonded DAP:T base pairs. Both unmodified and DAP-substituted DNA transitioned from a B- to an L-helix under physiological conditions of mild tension and unwinding. This transition avoids writhing and the ease of this transition may prevent cumbersome topological rearrangements in genomic DNA that would require topoisomerase activity to resolve. L-DNA displayed about tenfold lower persistence length than B-DNA. However, left-handed DAP-substituted DNA was twice as stiff as unmodified L-DNA. Unmodified DNA and DAP-substituted DNA have very distinct mechanical characteristics at physiological levels of negative supercoiling and tension.

Force spectroscopy with electromagnetic tweezers

Force spectroscopy using magnetic tweezers (MTs) is a powerful method to probe the physical characteristics of single polymers. Typically, molecules are functionalized for specific attachment to a glass surface at one end and a micrometer-scale paramagnetic bead at the other end. By applying an external magnetic field, multiple molecules can be stretched and twisted simultaneously without exposure to potentially damaging radiation. The majority of MTs utilize mobile, permanent magnets to produce forces on the beads (and the molecule under test). However, translating and rotating the permanent magnets may require expensive precision actuators, limit the rate at which force can be changed, and may induce vibrations that disturb tether dynamics and bead tracking. Alternatively, the magnetic field can be produced with an electromagnet, which allows fast force modulation and eliminates motor-associated vibration. Here, we describe a low-cost quadrapolar electromagnetic tweezer design capable of manipulating DNA-tethered MyOne paramagnetic beads with forces as high as 15 pN. The solid-state nature of the generated B-field modulated along two axes is convenient for accessing the range of forces and torques relevant for studying the activity of DNA motor enzymes like polymerases and helicases. Our design specifically leverages technology available at an increasing number of university maker spaces and student-run machine shops. Thus, it is an accessible tool for undergraduate education that is applicable to a wide range of biophysical research questions.

Single-molecule insights into torsion and roadblocks in bacterial transcript elongation

During transcription, RNA polymerase (RNAP) translocates along the helical template DNA while maintaining high transcriptional fidelity. However, all genomes are dynamically twisted, writhed, and decorated by bound proteins and motor enzymes. In prokaryotes, proteins bound to DNA, specifically or not, frequently compact DNA into conformations that may silence genes by obstructing RNAP. Collision of RNAPs with these architectural proteins, may result in RNAP stalling and/or displacement of the protein roadblock. It is important to understand how rapidly transcribing RNAPs operate under different levels of supercoiling or in the presence of roadblocks. Given the broad range of asynchronous dynamics exhibited by transcriptional complexes, single-molecule assays, such as atomic force microscopy, fluorescence detection, optical and magnetic tweezers, etc. are well suited for detecting and quantifying activity with adequate spatial and temporal resolution. Here, we summarize current understanding of the effects of torsion and roadblocks on prokaryotic transcription, with a focus on single-molecule assays that provide real-time detection and readout.

Radiation-Induced Effect on Spin-Selective Electron Transfer through Self-Assembled Monolayers of ds-DNA

Stability of the DNA molecule is essential for the proper functioning and sustainability of all living organisms. In this study, we investigate the effect of gamma radiation (γ-radiation) on spin-selective electron transfer through double strand (ds)DNA molecules. Self-assembled monolayers (SAMs) of 21-base long DNA are prepared on Au-coated Ni thin film. We measure the spin polarization (%) of the SAMs of ds-DNA using the spin-dependent electrochemical technique. We use a Cs-based γ-radiation source to expose the SAMs of ds-DNA immobilized on thin films for various time intervals ranging from 0–30 min. The susceptibility of DNA to γ-radiation is measured by spin-dependent electrochemistry. We observe that the efficiency of spin filtering by ds-DNA gradually decreases when exposure (to γ-radiation) time increases, and drops below 1% after 30 min of exposure. The change in spin polarization value is related either to the conformational perturbation in DNA or to structural damage in DNA molecules caused by ionizing radiation.

WASP: a software package for correctly characterizing the topological development of ribbon structures

We introduce the Writhe Application Software Package (WASP) which can be used to characterisze the topology of ribbon structures, the underlying mathematical model of DNA, Biopolymers, superfluid vorticies, elastic ropes and magnetic flux ropes. This characterization is achieved by the general twist–writhe decomposition of both open and closed ribbons, in particular through a quantity termed the polar writhe. We demonstrate how this decomposition is far more natural and straightforward than artificial closure methods commonly utilized in DNA modelling. In particular, we demonstrate how the decomposition of the polar writhe into local and non-local components distinctly characterizes the local helical structure and knotting/linking of the ribbon. This decomposition provides additional information not given by alternative approaches. As example applications, the WASP routines are used to characterise the evolving topology (writhe) of DNA minicircle and open ended plectoneme formation magnetic/optical tweezer simulations, and it is shown that the decomponsition into local and non-local components is particularly important for the detection of plectonemes. Finally it is demonstrated that a number of well known alternative writhe expressions are actually simplifications of the polar writhe measure.

Decision-Making in Cascade Complexes Harboring crRNAs of Altered Length

The multi-subunit type I CRISPR-Cas surveillance complex Cascade uses its crRNA to recognize dsDNA targets. Recognition involves DNA unwinding and base-pairing between the crRNA spacer region and a complementary DNA strand, resulting in formation of an R-loop structure. The modular Cascade architecture allows assembly of complexes containing crRNAs with altered spacer lengths that promise increased target specificity in emerging biotechnological applications. Here we produce type I-E Cascade complexes containing crRNAs with up to 57-nt-long spacers. We show that these complexes form R-loops corresponding to the designed target length, even for the longest spacers tested. Furthermore, the complexes can bind their targets with much higher affinity compared with the wild-type form. However, target recognition and the subsequent Cas3-mediated DNA cleavage do not require extended R-loops but already occur for wild-type-sized R-loops. These findings set important limits for specificity improvements of type I CRISPR-Cas systems.

Primed CRISPR adaptation in Escherichia coli cells does not depend on conformational changes in the Cascade effector complex detected in Vitro

In type I CRISPR-Cas systems, primed adaptation of new spacers into CRISPR arrays occurs when the effector Cascade-crRNA complex recognizes imperfectly matched targets that are not subject to efficient CRISPR interference. Thus, primed adaptation allows cells to acquire additional protection against mobile genetic elements that managed to escape interference. Biochemical and biophysical studies suggested that Cascade-crRNA complexes formed on fully matching targets (subject to efficient interference) and on partially mismatched targets that promote primed adaption are structurally different. Here, we probed Escherichia coli Cascade-crRNA complexes bound to matched and mismatched DNA targets using a magnetic tweezers assay. Significant differences in complex stabilities were observed consistent with the presence of at least two distinct conformations. Surprisingly, in vivo analysis demonstrated that all mismatched targets stimulated robust primed adaptation irrespective of conformational states observed in vitro. Our results suggest that primed adaptation is a direct consequence of a reduced interference efficiency and/or rate and is not a consequence of distinct effector complex conformations on target DNA.

Anticooperative Binding Governs the Mechanics of Ethidium-Complexed DNA

DNA intercalators bind nucleic acids by stacking between adjacent basepairs. This causes a considerable elongation of the DNA backbone as well as untwisting of the double helix. In the past few years, single-molecule mechanical experiments have become a common tool to characterize these deformations and to quantify important parameters of the intercalation process. Parameter extraction typically relies on the neighbor-exclusion model, in which a bound intercalator prevents intercalation into adjacent sites. Here, we challenge the neighbor-exclusion model by carefully quantifying and modeling the force-extension and twisting behavior of single ethidium-complexed DNA molecules. We show that only an anticooperative ethidium binding that allows for a disfavored but nonetheless possible intercalation into nearest-neighbor sites can consistently describe the mechanical behavior of intercalator-bound DNA. At high ethidium concentrations and elevated mechanical stress, this causes an almost complete occupation of nearest-neighbor sites and almost a doubling of the DNA contour length. We furthermore show that intercalation into nearest-neighbor sites needs to be considered when estimating intercalator parameters from zero-stress elongation and twisting data. We think that the proposed anticooperative binding mechanism may also be applicable to other intercalating molecules.

"End-to-end" stacking of small dsRNA

PELDOR (pulsed electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the "end-to-end" stacking of small double-stranded (ds) RNAs. For this study, the dsRNA molecules were only singly labeled with the spin label TPA to avoid multispin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From these data, the dissociation constant K_d was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we demonstrate by PELDOR experiments quantitatively.

Discontinuous growth of DNA plectonemes due to atomic scale friction

We develop a model to explain discontinuities in the increase of the length of a DNA plectoneme when the DNA filament is continuously twisted under tension. We account for DNA elasticity, electrostatic interactions and entropic effects due to thermal fluctuation. We postulate that a corrugated energy landscape that contains energy barriers is the cause of jumps in the length of the plectoneme as the number of turns is increased. Thus, our model is similar to the Prandtl-Tomlinson model of atomic scale friction. The existence of a corrugated energy landscape can be justified due to the close proximity of the neighboring pieces of DNA in a plectoneme. We assume the corrugated energy landscape to be sinusoidal since the plectoneme has a periodic helical structure and rotation of the bead is a form of periodic motion. We perform calculations with different tensile forces and ionic concentrations, and show that rotation-extension curves manifest stair-step shapes under relatively high ionic concentrations and high forces. We show that the jump in the plectonemic growth is caused by the flattening of the energy barrier in the corrugated landscape.

The pH-dependent elastic properties of nanoscale DNA films and the resultant bending signals for microcantilever biosensors

The diverse mechanical properties of nanoscale DNA films on solid substrates have a close correlation with complex detection signals of micro-/nano-devices. This paper is devoted to formulating several multiscale models to study the effect of pH-dependent ionic inhomogeneity on the graded elastic properties of nanoscale DNA films and the resultant bending deflections of microcantilever biosensors. First, a modified inverse Debye length is introduced to improve the classical Poisson-Boltzmann equation for the electrical potential of DNA films to consider the inhomogeneous effect of hydrogen ions. Second, the graded characteristics of the particle distribution are taken into consideration for an improvement in Parsegian's mesoscopic potential for both attraction-dominated and repulsion-dominated films. Third, by the improved interchain interaction potential and the thought experiment about the compression of a macroscopic continuum DNA bar, we investigate the diversity of the elastic properties of single-stranded DNA (ssDNA) films due to pH variations. The relevant theoretical predictions quantitatively or qualitatively agree well with the relevant DNA experiments on the electrical potential, film thickness, condensation force, elastic modulus, and microcantilever deflections. The competition between attraction and repulsion among the fixed charges and the free ions endows the DNA film with mechanical properties such as a remarkable size effect and a non-monotonic behavior, and a negative elastic modulus is first revealed in the attraction-dominated ssDNA film. There exists a transition between the pH-sensitive parameter interval and the pH-insensitive one for the bending signals of microcantilevers, which is predominated by the initial stress effect in the DNA film.

How proteins bind to DNA: target discrimination and dynamic sequence search by the telomeric protein TRF1

Target search as performed by DNA-binding proteins is a complex process, in which multiple factors contribute to both thermodynamic discrimination of the target sequence from overwhelmingly abundant off-target sites and kinetic acceleration of dynamic sequence interrogation. TRF1, the protein that binds to telomeric tandem repeats, faces an intriguing variant of the search problem where target sites are clustered within short fragments of chromosomal DNA. In this study, we use extensive (>0.5 ms in total) MD simulations to study the dynamical aspects of sequence-specific binding of TRF1 at both telomeric and non-cognate DNA. For the first time, we describe the spontaneous formation of a sequence-specific native protein-DNA complex in atomistic detail, and study the mechanism by which proteins avoid off-target binding while retaining high affinity for target sites. Our calculated free energy landscapes reproduce the thermodynamics of sequence-specific binding, while statistical approaches allow for a comprehensive description of intermediate stages of complex formation.

Probing the salt dependence of the torsional stiffness of DNA by multiplexed magnetic torque tweezers

The mechanical properties of DNA fundamentally constrain and enable the storage and transmission of genetic information and its use in DNA nanotechnology. Many properties of DNA depend on the ionic environment due to its highly charged backbone. In particular, both theoretical analyses and direct single-molecule experiments have shown its bending stiffness to depend on salt concentration. In contrast, the salt-dependence of the twist stiffness of DNA is much less explored. Here, we employ optimized multiplexed magnetic torque tweezers to study the torsional stiffness of DNA under varying salt conditions as a function of stretching force. At low forces (<3 pN), the effective torsional stiffness is ∼10% smaller for high salt conditions (500 mM NaCl or 10 mM MgCl₂) compared to lower salt concentrations (20 mM NaCl and 100 mM NaCl). These differences, however, can be accounted for by taking into account the known salt dependence of the bending stiffness. In addition, the measured high-force (6.5 pN) torsional stiffness values of C = 103 ± 4 nm are identical, within experimental errors, for all tested salt concentration, suggesting that the intrinsic torsional stiffness of DNA does not depend on salt.

Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.

Understanding the origin of liquid crystal ordering of ultrashort double-stranded DNA

Recent experiments have shown that short double-stranded DNA (dsDNA) fragments having six- to 20-base pairs exhibit various liquid crystalline phases. This violates the condition of minimum molecular shape anisotropy that analytical theories demand for liquid crystalline ordering. It has been hypothesized that the liquid crystalline ordering is the result of end-to-end stacking of dsDNA to form long supramolecular columns which satisfy the shape anisotropy criterion necessary for ordering. To probe the thermodynamic feasibility of this process, we perform molecular dynamics simulations on ultrashort (four base pair long) dsDNA fragments, quantify the strong end-to-end attraction between them, and demonstrate that the nematic ordering of the self-assembled stacked columns is retained for a large range of temperature and salt concentration.

Single-Molecule Insight Into Target Recognition by CRISPR-Cas Complexes

Ribonucleoprotein (RNP) complexes from CRISPR-Cas systems have attracted enormous interest since they can be easily and flexibly reprogrammed to target any desired locus for genome engineering and gene regulation applications. Basis for the programmability is a short RNA (crRNA) inside these complexes that recognizes the target nucleic acid by base pairing. For CRISPR-Cas systems that target double-stranded DNA this results in local DNA unwinding and formation of a so-called R-loop structure. Here we provide an overview how this target recognition mechanism can be dissected in great detail at the level of a single molecule. Specifically, we demonstrate how magnetic tweezers are applied to measure the local DNA unwinding at the target in real time. To this end we introduce the technique and the measurement principle. By studying modifications of the consensus target sequence, we show how different sequence elements contribute to the target recognition mechanism. From these data, a unified target recognition mechanism can be concluded for the RNPs Cascade and Cas9 from types I and II CRISPR-Cas systems. R-loop formation is hereby initiated on the target at an upstream element, called protospacer adjacent motif (PAM), from which the R-loop structure zips directionally toward the PAM-distal end of the target. At mismatch positions, the R-loop propagation stalls and further propagation competes with collapse of the structure. Upon full R-loop zipping conformational changes within the RNPs trigger degradation of the DNA target. This represents a shared labor mechanism in which zipping between nucleic acid strands is the actual target recognition mechanism while sensing of the R-loop arrival at the PAM-distal end just verifies the success of the full zipping.

Evidence of protein-free homology recognition in magnetic bead force-extension experiments

Earlier theoretical studies have proposed that the homology-dependent pairing of large tracts of dsDNA may be due to physical interactions between homologous regions. Such interactions could contribute to the sequence-dependent pairing of chromosome regions that may occur in the presence or the absence of double-strand breaks. Several experiments have indicated the recognition of homologous sequences in pure electrolytic solutions without proteins. Here, we report single-molecule force experiments with a designed 60 kb long dsDNA construct; one end attached to a solid surface and the other end to a magnetic bead. The 60 kb constructs contain two 10 kb long homologous tracts oriented head to head, so that their sequences match if the two tracts fold on each other. The distance between the bead and the surface is measured as a function of the force applied to the bead. At low forces, the construct molecules extend substantially less than normal, control dsDNA, indicating the existence of preferential interaction between the homologous regions. The force increase causes no abrupt but continuous unfolding of the paired homologous regions. Simple semi-phenomenological models of the unfolding mechanics are proposed, and their predictions are compared with the data.

Probing the mechanical properties, conformational changes, and interactions of nucleic acids with magnetic tweezers

Nucleic acids are central to the storage and transmission of genetic information. Mechanical properties, along with their sequence, both enable and fundamentally constrain the biological functions of DNA and RNA. For small deformations from the equilibrium conformations, nucleic acids are well described by an isotropic elastic rod model. However, external forces and torsional strains can induce conformational changes, giving rise to a complex force-torque phase diagram. This review focuses on magnetic tweezers as a powerful tool to precisely determine both the elastic parameters and conformational transitions of nucleic acids under external forces and torques at the single-molecule level. We review several variations of magnetic tweezers, in particular conventional magnetic tweezers, freely orbiting magnetic tweezers and magnetic torque tweezers, and discuss their characteristic capabilities. We then describe the elastic rod model for DNA and RNA and discuss conformational changes induced by mechanical stress. The focus lies on the responses to torque and twist, which are crucial in the mechanics and interactions of nucleic acids and can directly be measured using magnetic tweezers. We conclude by highlighting several recent studies of nucleic acid-protein and nucleic acid-small-molecule interactions as further applications of magnetic tweezers and give an outlook of some exciting developments to come.

TrmBL2 from Pyrococcus furiosus Interacts Both with Double-Stranded and Single-Stranded DNA

In many hyperthermophilic archaea the DNA binding protein TrmBL2 or one of its homologues is abundantly expressed. TrmBL2 is thought to play a significant role in modulating the chromatin architecture in combination with the archaeal histone proteins and Alba. However, its precise physiological role is poorly understood. It has been previously shown that upon binding TrmBL2 covers double-stranded DNA, which leads to the formation of a thick and fibrous filament. Here we investigated the filament formation process as well as the stabilization of DNA by TrmBL2 from Pyroccocus furiosus in detail. We used magnetic tweezers that allow to monitor changes of the DNA mechanical properties upon TrmBL2 binding on the single-molecule level. Extended filaments formed in a cooperative manner and were considerably stiffer than bare double-stranded DNA. Unlike Alba, TrmBL2 did not form DNA cross-bridges. The protein was found to bind double- and single-stranded DNA with similar affinities. In mechanical disruption experiments of DNA hairpins this led to stabilization of both, the double- (before disruption) and the single-stranded (after disruption) DNA forms. Combined, these findings suggest that the biological function of TrmBL2 is not limited to modulating genome architecture and acting as a global repressor but that the protein acts additionally as a stabilizer of DNA secondary structure.

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

Introducing improved structural properties and salt dependence into a coarse-grained model of DNA

We introduce an extended version of oxDNA, a coarse-grained model of deoxyribonucleic acid (DNA) designed to capture the thermodynamic, structural, and mechanical properties of single- and double-stranded DNA. By including explicit major and minor grooves and by slightly modifying the coaxial stacking and backbone-backbone interactions, we improve the ability of the model to treat large (kilobase-pair) structures, such as DNA origami, which are sensitive to these geometric features. Further, we extend the model, which was previously parameterised to just one salt concentration ([Na(+)] = 0.5M), so that it can be used for a range of salt concentrations including those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by presenting results from simulations of the structure of large DNA objects and by using the model to investigate some salt-dependent properties of DNA.

Extending the range for force calibration in magnetic tweezers

Magnetic tweezers are a wide-spread tool used to study the mechanics and the function of a large variety of biomolecules and biomolecular machines. This tool uses a magnetic particle and a strong magnetic field gradient to apply defined forces to the molecule of interest. Forces are typically quantified by analyzing the lateral fluctuations of the biomolecule-tethered particle in the direction perpendicular to the applied force. Since the magnetic field pins the anisotropy axis of the particle, the lateral fluctuations follow the geometry of a pendulum with a short pendulum length along and a long pendulum length perpendicular to the field lines. Typically, the short pendulum geometry is used for force calibration by power-spectral-density (PSD) analysis, because the movement of the bead in this direction can be approximated by a simple translational motion. Here, we provide a detailed analysis of the fluctuations according to the long pendulum geometry and show that for this direction, both the translational and the rotational motions of the particle have to be considered. We provide analytical formulas for the PSD of this coupled system that agree well with PSDs obtained in experiments and simulations and that finally allow a faithful quantification of the magnetic force for the long pendulum geometry. We furthermore demonstrate that this methodology allows the calibration of much larger forces than the short pendulum geometry in a tether-length-dependent manner. In addition, the accuracy of determination of the absolute force is improved. Our force calibration based on the long pendulum geometry will facilitate high-resolution magnetic-tweezers experiments that rely on short molecules and large forces, as well as highly parallelized measurements that use low frame rates.

Pore Translocation of Knotted Polymer Chains: How Friction Depends on Knot Complexity

Knots can affect the capability of polymers to translocate through narrow pores in complex and counterintuitive ways that are still relatively unexplored. We report here on a systematic theoretical and computational investigation of the driven translocation of flexible chains accommodating a large repertoire of knots trapped at the pore entrance. These include composite knots, which are the most common form of spontaneous entanglement in long polymers. Two unexpected results emerge from this study. First, the high force translocation compliance does not decrease systematically with knot complexity. Second, the response of composite knots is so dependent on the order of their factor knots, that their hindrance can even be lower than some of their prime components. We show that the resulting rich and seemingly disparate phenomenology can be captured in a seamless framework based on the mechanism by which the tractive force is propagated along and past the knots. The quantitative scheme can be viably used for predictive purposes and, hence, ought to be useful in applicative contexts, too.

Topological patterns in two-dimensional gel electrophoresis of DNA knots

Gel electrophoresis is a powerful experimental method to probe the topology of DNA and other biopolymers. Although there is a large body of experimental work that allows us to accurately separate different topoisomers of a molecule, a full theoretical understanding of these experiments has not yet been achieved. Here we show that the mobility of DNA knots depends crucially and subtly on the physical properties of the gel and, in particular, on the presence of dangling ends. The topological interactions between these and DNA molecules can be described in terms of an "entanglement number" and yield a nonmonotonic mobility at moderate fields. Consequently, in 2D electrophoresis, gel bands display a characteristic arc pattern; this turns into a straight line when the density of dangling ends vanishes. We also provide a novel framework to accurately predict the shape of such arcs as a function of molecule length and topological complexity, which may be used to inform future experiments.

Changing chromatin fiber conformation by nucleosome repositioning

Chromatin conformation is dynamic and heterogeneous with respect to nucleosome positions, which can be changed by chromatin remodeling complexes in the cell. These molecular machines hydrolyze ATP to translocate or evict nucleosomes, and establish loci with regularly and more irregularly spaced nucleosomes as well as nucleosome-depleted regions. The impact of nucleosome repositioning on the three-dimensional chromatin structure is only poorly understood. Here, we address this issue by using a coarse-grained computer model of arrays of 101 nucleosomes considering several chromatin fiber models with and without linker histones, respectively. We investigated the folding of the chain in dependence of the position of the central nucleosome by changing the length of the adjacent linker DNA in basepair steps. We found in our simulations that these translocations had a strong effect on the shape and properties of chromatin fibers: i), Fiber curvature and flexibility at the center were largely increased and long-range contacts between distant nucleosomes on the chain were promoted. ii), The highest destabilization of the fiber conformation occurred for a nucleosome shifted by two basepairs from regular spacing, whereas effects of linker DNA changes of ?10 bp in phase with the helical twist of DNA were minimal. iii), A fiber conformation can stabilize a regular spacing of nucleosomes inasmuch as favorable stacking interactions between nucleosomes are facilitated. This can oppose nucleosome translocations and increase the energetic costs for chromatin remodeling. Our computational modeling framework makes it possible to describe the conformational heterogeneity of chromatin in terms of nucleosome positions, and thus advances theoretical models toward a better understanding of how genome compaction and access are regulated within the cell.

Directional R-Loop Formation by the CRISPR-Cas Surveillance Complex Cascade Provides Efficient Off-Target Site Rejection

CRISPR-Cas systems provide bacteria and archaea with adaptive immunity against foreign nucleic acids. In type I CRISPR-Cas systems, invading DNA is detected by a large ribonucleoprotein surveillance complex called Cascade. The crRNA component of Cascade is used to recognize target sites in foreign DNA (protospacers) by formation of an R-loop driven by base-pairing complementarity. Using single-molecule supercoiling experiments with near base-pair resolution, we probe here the mechanism of R-loop formation and detect short-lived R-loop intermediates on off-target sites bearing single mismatches. We show that R-loops propagate directionally starting from the protospacer-adjacent motif (PAM). Upon reaching a mismatch, R-loop propagation stalls and collapses in a length-dependent manner. This unambiguously demonstrates that directional zipping of the R-loop accomplishes efficient target recognition by rapidly rejecting binding to off-target sites with PAM-proximal mutations. R-loops that reach the protospacer end become locked to license DNA degradation by the auxiliary Cas3 nuclease/helicase without further target verification.

Close encounters with DNA

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

Driving knots on DNA with AC/DC electric fields: topological friction and memory effects

The dynamical properties of entangled polyelectrolytes are investigated theoretically and computationally for a proposed novel micromanipulation setup. Specifically, we investigate the effects of DC and AC electric fields acting longitudinally on knotted DNA chains, modelled as semiflexible chains of charged beads, under mechanical tension. We consider various experimentally accessible values of the field amplitude and frequency as well as several of the simplest knot types. In particular, we consider both torus and twist knots because they are respectively known to be able or unable to slide along macroscopic threads and ropes. Strikingly, this qualitative distinction disappears in this microscopic context because all the considered knot types acquire a systematic drift in the direction of the electric force. Notably, the knot drift velocity and diffusion coefficient in zero field (both measurable also experimentally) can be used to define a characteristic "frictional" lengthscale for the various knot types. This previously unexplored length provides valuable information on the extent of self-interactions in the nominal knotted region. It is finally observed that the motion of a knot can effectively follow the AC field only if the driving period is larger than the knot relaxation time (for which the self-diffusion time provides an upper bound). These results suggest that salient aspects of the intrinsic dynamics of knots in DNA chains could be probed experimentally by means of external, time-dependent electric fields.

Direct observation of R-loop formation by single RNA-guided Cas9 and Cascade effector complexes

Clustered, regularly interspaced, short palindromic repeats (CRISPR)/CRISPR-associated (Cas) systems protect bacteria and archaea from infection by viruses and plasmids. Central to this defense is a ribonucleoprotein complex that produces RNA-guided cleavage of foreign nucleic acids. In DNA-targeting CRISPR-Cas systems, the RNA component of the complex encodes target recognition by forming a site-specific hybrid (R-loop) with its complement (protospacer) on an invading DNA while displacing the noncomplementary strand. Subsequently, the R-loop structure triggers DNA degradation. Although these reactions have been reconstituted, the exact mechanism of R-loop formation has not been fully resolved. Here, we use single-molecule DNA supercoiling to directly observe and quantify the dynamics of torque-dependent R-loop formation and dissociation for both Cascade- and Cas9-based CRISPR-Cas systems. We find that the protospacer adjacent motif (PAM) affects primarily the R-loop association rates, whereas protospacer elements distal to the PAM affect primarily R-loop stability. Furthermore, Cascade has higher torque stability than Cas9 by using a conformational locking step. Our data provide direct evidence for directional R-loop formation, starting from PAM recognition and expanding toward the distal protospacer end. Moreover, we introduce DNA supercoiling as a quantitative tool to explore the sequence requirements and promiscuities of orthogonal CRISPR-Cas systems in rapidly emerging gene-targeting applications.

Undulations in a weakly interacting mechanically generated molecular braid under tension

We consider mechanically generated molecular braids composed of two molecules where long range interactions between them can be considered to be very weak. We describe a model that takes account of the thermal fluctuations of the braid, steric interactions between the molecules, and external mechanical forces. In this model, both sets of ends, of the two molecules, are considered to be separated by a fixed distance much larger than the radius of the braid. One set of ends is rotated to generate a braid of a certain number of pitches (or turns), while the other set remains fixed. This model may describe the situation in which the ends of each molecule are attached to a substrate and a magnetic bead; to the latter a pulling force and rotational torque can be applied. We discuss various aspects of our model. Most importantly, an expression for the free energy is given, from which equations, determining the various geometric parameters of the braid, can be obtained. By numerically solving these equations, we give predictions from the model for the external torque needed to produce a braid with a certain number of turns per bending persistence length, as well as the end to end extension of the two molecules for a given pulling force. Other geometric parameters, as well as the lateral force required to keep the ends of the two molecules apart, are also calculated.

Self-consistent treatment of electrostatics in molecular DNA braiding through external forces

In this paper we consider a physical system in which two DNA molecules braid about each other. The distance between the two molecular ends, on either side of the braid, is held at a distance much larger than supercoiling radius of the braid. The system is subjected to an external pulling force, and a moment that induces the braiding. In a model, developed for understanding such a system, we assume that each molecule can be divided into a braided and unbraided section. We also suppose that the DNA is nicked so that there is no constraint of the individual linking numbers of the molecules. Included in the model are steric and electrostatic interactions, thermal fluctuations of the braided and unbraided sections of the molecule, as well as the constraint on the braid linking (catenation) number. We compare two approximations used in estimating the free energy of the braided section. One is where the amplitude of undulations of one molecule with respect to the other is determined only by steric interactions. The other is a self-consistent determination of the mean-squared amplitude of these undulations. In this second approximation electrostatics should play an important role in determining this quantity, as suggested by physical arguments. We see that if the electrostatic interaction is sufficiently large there are indeed notable differences between the two approximations. We go on to test the self-consistent approximation-included in the full model-against experimental data for such a system, and we find good agreement. However, there seems to be a slight left-right-handed braid asymmetry in some of the experimental results. We discuss what might be the origin of this small asymmetry.

Understanding nucleic acid-ion interactions

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

Origin of current blockades in nanopore translocation experiments

We present a detailed investigation of the ionic current in a cylindrical model nanopore in the absence and the presence of a double stranded DNA homopolymer. Our atomistic simulations are capable of reproducing almost exactly the experimental data obtained by Smeets et al., including notably the crossover salt concentration that yields equal current measurements in both situations. We can rule out that the observed current blockade is due to the steric exclusion of charge carriers from the DNA, since for all investigated salt concentrations the charge carrier density is higher when the DNA is present. Calculations using a mean-field electrokinetic model proposed by van Dorp et al. fail quantitatively in predicting this effect. We can relate the shortcomings of the mean-field model to a surface related molecular drag that the ions feel in the presence of the DNA. This drag is independent of the salt concentration and originates from electrostatic, hydrodynamic, and excluded volume interactions.

Multiplectoneme phase of double-stranded DNA under torsion

We use the wormlike chain model to study supercoiling of DNA under tension and torque. The model reproduces experimental data for a broad range of forces, salt concentrations, and contour lengths. We find a plane of first-order phase transitions ending in a smeared-out line of critical points, the multiplectoneme phase, which is characterized by a fast twist-mediated diffusion of plectonemes and a torque that rises after plectoneme formation with increasing linking number. The discovery of this phase at the same time resolves the discrepancies between existing models and experiment.