Biophysics of protein-DNA interactions and chromosome organization

Nucleosomes play a dual role in regulating transcription dynamics

Transcription has a mechanical component, as the translocation of the transcription machinery or RNA polymerase (RNAP) on DNA or chromatin is dynamically coupled to the chromatin torsion. This posits chromatin mechanics as a possible regulator of eukaryotic transcription, however, the modes and mechanisms of this regulation are elusive. Here, we first take a statistical mechanics approach to model the torsional response of topology-constrained chromatin. Our model recapitulates the experimentally observed weaker torsional stiffness of chromatin compared to bare DNA and proposes structural transitions of nucleosomes into chirally distinct states as the driver of the contrasting torsional mechanics. Coupling chromatin mechanics with RNAP translocation in stochastic simulations, we reveal a complex interplay of DNA supercoiling and nucleosome dynamics in governing RNAP velocity. Nucleosomes play a dual role in controlling the transcription dynamics. The steric barrier aspect of nucleosomes in the gene body counteracts transcription via hindering RNAP motion, whereas the chiral transitions facilitate RNAP motion via driving a low restoring torque upon twisting the DNA. While nucleosomes with low dissociation rates are typically transcriptionally repressive, highly dynamic nucleosomes offer less of a steric barrier and enhance the transcription elongation dynamics of weakly transcribed genes via buffering DNA twist. We use the model to predict transcription-dependent levels of DNA supercoiling in segments of the budding yeast genome that are in accord with available experimental data. The model unveils a paradigm of DNA supercoiling-mediated interaction between genes and makes testable predictions that will guide experimental design.

Twist-stretch relations in nucleic acids

Nucleic acids are highly deformable helical molecules constantly stretched, twisted and bent in their biological functioning. Single molecule experiments have shown that double stranded (ds)-RNA and standard ds-DNA have opposite twist-stretch patterns and stretching properties when overwound under a constant applied load. The key structural features of the A-form RNA and B-form DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian model which accounts for the radial, bending and twisting fluctuations of the base pairs. Using path integral techniques which sum over the ensemble of the base pair fluctuations, I compute the average helical repeat of the molecules as a function of the load. The obtained twist-stretch relations and stretching properties, for short A- and B-helical fragments, are consistent with the opposite behaviors observed in kilo-base long molecules.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Nonaffine Mechanics of Entangled Networks Inspired by Intermediate Filaments

Inspired by massive intermediate filament (IF) reorganization in superstretched epithelia, we examine computationally the principles controlling the mechanics of a set of entangled filaments whose ends slide on the cell boundary. We identify an entanglement metric and threshold beyond which random loose networks respond nonaffinely and nonlinearly to stretch by self-organizing into structurally optimal star-shaped configurations. A simple model connecting cellular and filament strains links emergent mechanics to cell geometry, network topology, and filament mechanics. We identify a safety net mechanism in IF networks and provide a framework to harness entanglement in soft fibrous materials.

Direct observation of surface charge and stiffness of human metaphase chromosomes

Metaphase chromosomes in which both polynucleotides and proteins are condensed with hierarchies are closely related to life phenomena such as cell division, cancer development, and cellular senescence. Nevertheless, their nature is rarely revealed, owing to their structural complexity and technical limitations in analytical methods. In this study, we used surface potential and nanomechanics mapping technology based on atomic force microscopy to measure the surface charge and intrinsic stiffness of metaphase chromosomes. We found that extra materials covering the chromosomes after the extraction process were positively charged. With the covering materials, the chromosomes were positively charged (ca. 44.9 ± 16.48 mV) and showed uniform stiffness (ca. 6.23 ± 1.98 MPa). In contrast, after getting rid of the extra materials through treatment with RNase and protease, the chromosomes were strongly negatively charged (ca. -197.4 ± 77.87 mV) and showed relatively non-uniform and augmented stiffness (ca. 36.87 ± 17.56 MPa). The results suggested undulating but compact coordination of condensed chromosomes. Additionally, excessive treatment with RNase and protease could destroy the chromosomal structure, providing an exceptional opportunity for multiscale stiffness mapping of polynucleotides, nucleosomes, chromatin fibers, and chromosomes in a single image. Our approach offers a new horizon in terms of an analytical technique for studying chromosome-related diseases.

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.

Equilibrium fluctuations of DNA plectonemes

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e., over- or underwound. They are ubiquitous in cellular DNA, and their physical properties have attracted significant interest both from the experimental side and from the modeling side. In this paper, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched phase and a plectonemic phase. A variety of mechanisms are found to contribute to extension fluctuations, characterized by the variance 〈Δz^{2}〉. We find the dominant contribution to 〈Δz^{2}〉 to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the twistable wormlike chain and analyze the fluctuation of various quantities, the results of which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations.

Shaping the genome via lengthwise compaction, phase separation, and lamina adhesion

The link between genomic structure and biological function is yet to be consolidated, it is, however, clear that physical manipulation of the genome, driven by the activity of a variety of proteins, is a crucial step. To understand the consequences of the physical forces underlying genome organization, we build a coarse-grained polymer model of the genome, featuring three fundamentally distinct classes of interactions: lengthwise compaction, i.e., compaction of chromosomes along its contour, self-adhesion among epigenetically similar genomic segments, and adhesion of chromosome segments to the nuclear envelope or lamina. We postulate that these three types of interactions sufficiently represent the concerted action of the different proteins organizing the genome architecture and show that an interplay among these interactions can recapitulate the architectural variants observed across the tree of life. The model elucidates how an interplay of forces arising from the three classes of genomic interactions can drive drastic, yet predictable, changes in the global genome architecture, and makes testable predictions. We posit that precise control over these interactions in vivo is key to the regulation of genome architecture.

DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases

Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or ‘push’ between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate by modifying the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the DNA torsional response. Using stochastic simulations, we demonstrate DNA supercoiling-mediated long-range cooperation between co-transcribing RNAPs. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern transcriptional dynamics.

Impact of linker histone in the formation of ambochlorin-calf thymus DNA complex: Multi-spectroscopic, stopped-flow, and molecular modeling approaches

Objective(s):

This study aimed to evaluate the role of the linker histone (H₁) in the binding interaction between ambochlorin (Amb), and calf thymus DNA (ctDNA) as binary and ternary systems.

Materials and Methods:

The project was accomplished through the means of absorbance, fluorescence, stopped-flow circular dichroism spectroscopy, viscosity, thermal melting, and molecular modeling techniques.

Results:

Spectroscopic analysis revealed that although Amb was strongly bound to both ctDNA and ctDNA-H₁, it showed a greater tendency to ctDNA in the presence of the linker histone. The obtained thermodynamic parameters revealed that both Amb-ctDNA and Amb-ctDNA-H₁ interactions were spontaneous, endothermic, and entropy-favored, and hydrophobic interactions played the main role in the formation and stabilization of complexes. Analysis of the stopped-flow circular dichroism results revealed that the binding process of Amb-ctDNA and Amb-ctDNA-H₁ required a time of more than 150 milliseconds to complete. Moreover, Amb-ctDNA complex formation was marginally decelerated in the presence of the linker histone. The docking results suggested that the presence of the linker histone may alter the binding sites of Amb from ctDNA minor grooves to major grooves.

Conclusion:

All quenching processes were governed by a dynamic mechanism. Additionally, Amb did not stabilize or induce considerable conformational changes in ctDNA and ctDNA-H₁ complex upon binding. In silico molecular docking results confirmed that Amb was bound to the double-helical ctDNA and ctDNA-H₁ via ctDNA grooves. In summary, some binding properties of the interactions between Amb and ctDNA change in the presence of the linker histone.

Base pair fluctuations in helical models for nucleic acids

A statistical method is developed to estimate the maximum amplitude of the base pair fluctuations in a three dimensional mesoscopic model for nucleic acids. The base pair thermal vibrations around the helix diameter are viewed as a Brownian motion for a particle embedded in a stable helical structure. The probability to return to the initial position is computed, as a function of time, by integrating over the particle paths consistent with the physical properties of the model potential. The zero time condition for the first-passage probability defines the constraint to select the integral cutoff for various macroscopic helical conformations, obtained by tuning the twist, bending, and slide motion between adjacent base pairs along the molecule stack. Applying the method to a short homogeneous chain at room temperature, we obtain meaningful estimates for the maximum fluctuations in the twist conformation with ∼10.5 base pairs per helix turn, typical of double stranded DNA helices. Untwisting the double helix, the base pair fluctuations broaden and the integral cutoff increases. The cutoff is found to increase also in the presence of a sliding motion, which shortens the helix contour length, a situation peculiar of dsRNA molecules.

Study on the interactions between B-norcholesteryl benzimidazole compounds with ct-DNA

The study of molecule-DNA interaction is very important for designing an improved therapeutic agent. In previous studies, we synthesized some B-norcholesteryl benzimidazole compounds, and the tests on cancer cells showed that these compounds had good in vitro anti-cancer activities. In order to further investigate mechanism of their actions, three different B-norcholesteryl benzimidazole compounds were selected and interaction of these compounds with the calf thymus DNA (ct-DNA) was monitored by using various methods including UV-Vis and fluorescence spectroscopic techniques, viscosity measurement, and circular dichroism (CD). The results proved a hypochromic effect accompanied with a slight red-shift due to the interaction of the molecules with ct-DNA. According to the UV-Vis and fluorescence spectra, the mentioned compounds were bound to DNA, preferentially through partial intercalation into the DNA helix. Moreover, the ethidium bromide (EB) and Hoechst 33258 competitive binding experiments were also used to confirm the interaction mode of the compounds with ct-DNA. In the Hoechst 33258 displacement experiment, no significant change in the fluorescence intensity was observed. Additional assays such as iodide quenching, viscosity, and CD spectroscopy further confirmed that intercalation should be the major binding mode of the selected compounds with DNA. The cytotoxicity of these three compounds was also evaluated by MTT method, and the results confirmed that binding ability of these compounds to DNA was consistent with their cytotoxicity behavior. The experimental results indicated a higher binding affinity for compound 3 compared to the other compounds. This research provided a better understanding on the molecular mechanism of the interaction between B-norcholesteryl benzimidazole compounds and tumor cells, and offered a beneficial perspective to the designation of novel B-norsteroidal anticancer compounds.

Dynamical control of denaturation bubble nucleation in supercoiled DNA minicircles

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

Properties of gene expression and chromatin structure with mechanically regulated elongation

In recent years, physical elements of transcription have emerged as central in our understanding of gene expression. Recent work has been done introducing a simple description of the basic physical elements of transcription where RNA elongation, RNA polymerase (RNAP) rotation and DNA supercoiling are coupled (1). Here we generalize this framework to accommodate the behavior of many RNAPs operating on multiple genes on a shared piece of DNA. The resulting framework is combined with well-established stochastic processes of transcription resulting in a model which characterizes the impact of the mechanical properties of transcription on gene expression and DNA structure. Transcriptional bursting readily emerges as a common phenomenon with origins in the geometric nature of the genetic system and results in the bounding of gene expression statistics. Properties of a multiple gene system are examined with special attention paid to the role that genome composition (gene orientation, size and intergenic distance) plays in the ability of genes to transcribe. The role of transcription in shaping DNA structure is examined and the possibility of transcription driven domain formation is discussed.

Twist-bend coupling and the statistical mechanics of the twistable wormlike-chain model of DNA: Perturbation theory and beyond

The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling G between twisting and bending, known to arise from the groove asymmetry of the DNA double helix. Here the effect of G on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness C_{eff} for small twist-bend coupling. We find that the "bare" G is "screened" by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with G=0 but with long-wavelength bending and twisting rigidities that are shifted by G-dependent amounts. Using results for torsional and bending rigidities for freely fluctuating DNA, we show how our perturbative results can be extended to a nonperturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo "triad" and molecular dynamics "oxDNA" models, characterized by different degrees of coarse graining, validating the perturbative and nonperturbative analyses. While our theory is in generally good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely used coarse-grained models.

Overtwisting induces polygonal shapes in bent DNA

By combining analytical results and simulations of various coarse-grained models, we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternation of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, rather than the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants.

DNA size in confined environments

For short DNA molecules in crowded environments, we evaluate macroscopic parameters such as the average end-to-end distance and the twist conformation by tuning the strength of the site specific confinement driven by the crowders.

Energy required to pinch a DNA plectoneme

DNA supercoiling plays an important role from a biological point of view. One of its consequences at the supramolecular level is the formation of DNA superhelices named plectonemes. Normally separated by a distance on the order of 10 nm, the two opposite double strands of a DNA plectoneme must be brought closer if a protein or protein complex implicated in genetic regulation is to be bound simultaneously to both strands, as if the plectoneme was locally pinched. We propose an analytic calculation of the energetic barrier, of elastic nature, required to bring closer the two loci situated on the opposed double strands. We examine how this energy barrier scales with the DNA supercoiling. For physically relevant values of elastic parameters and of supercoiling density, we show that the energy barrier is in the k_{B}T range under physiological conditions, thus demonstrating that the limiting step to loci encounter is more likely the preceding plectoneme slithering bringing the two loci side by side.

Bend-Induced Twist Waves and the Structure of Nucleosomal DNA

Recent work indicates that twist-bend coupling plays an important role in DNA micromechanics. Here we investigate its effect on bent DNA. We provide an analytical solution of the minimum-energy shape of circular DNA, showing that twist-bend coupling induces sinusoidal twist waves. This solution is in excellent agreement with both coarse-grained simulations of minicircles and nucleosomal DNA data, which is bent and wrapped around histone proteins in a superhelical conformation. Our analysis shows that the observed twist oscillation in nucleosomal DNA, so far attributed to the interaction with the histone proteins, is an intrinsic feature of free bent DNA, and should be observable in other protein-DNA complexes.

Ionic liquids make DNA rigid

Persistence length of double-stranded DNA (dsDNA) is known to decrease with an increase in ionic concentration of the solution. In contrast to this, here we show that the persistence length of dsDNA increases dramatically as a function of ionic liquid (IL) concentration. Using all atom explicit solvent molecular dynamics simulations and theoretical models, we present, for the first time, a systematic study to determine the mechanical properties of dsDNA in various hydrated ILs at different concentrations. We find that dsDNA in 50 wt % ILs have lower persistence length and stretch modulus in comparison to 80 wt % ILs. We further observe that both the persistence length and stretch modulus of dsDNA increase as we increase the concentration of ILs. The present trend of the stretch modulus and persistence length of dsDNA with IL concentration supports the predictions of the macroscopic elastic theory, in contrast to the behavior exhibited by dsDNA in monovalent salt. Our study further suggests the preferable ILs that can be used for maintaining DNA stability during long-term storage.

DNA elasticity from coarse-grained simulations: The effect of groove asymmetry

It is well established that many physical properties of DNA at sufficiently long length scales can be understood by means of simple polymer models. One of the most widely used elasticity models for DNA is the twistable worm-like chain (TWLC), which describes the double helix as a continuous elastic rod with bending and torsional stiffness. An extension of the TWLC, which has recently received some attention, is the model by Marko and Siggia, who introduced an additional twist-bend coupling, expected to arise from the groove asymmetry. By performing computer simulations of two available versions of oxDNA, a coarse-grained model of nucleic acids, we investigate the microscopic origin of twist-bend coupling. We show that this interaction is negligible in the oxDNA version with symmetric grooves, while it appears in the oxDNA version with asymmetric grooves. Our analysis is based on the calculation of the covariance matrix of equilibrium deformations, from which the stiffness parameters are obtained. The estimated twist-bend coupling coefficient from oxDNA simulations is G=30±1 nm. The groove asymmetry induces a novel twist length scale and an associated renormalized twist stiffness κ_t≈80 nm, which is different from the intrinsic torsional stiffness C≈110 nm. This naturally explains the large variations on experimental estimates of the intrinsic stiffness performed in the past.

Applications of Magnetic Tweezers to Studies of NAPs

Nucleoid-associated proteins (NAPs) are important factors in shaping bacterial nucleoid and regulating global gene expression. A great deal of insights into NAPs can be obtained through studies using single DNA molecule, which has been made possible owing to recent rapid development of single-DNA manipulation techniques. These studies provide information on modes of binding to DNA, which shed light on the mechanism underlying the regulatory function of NAPs. In addition, how NAPs organize DNA and thus their contribution to chromosomal DNA packaging can be determined. In this chapter, we introduce transverse magnetic tweezers that allows for convenient manipulation of long DNA molecules, and its applications in studies of NAPs as exemplified by the E. coli H-NS protein. We describe how transverse magnetic tweezers is a powerful tool that can be used to characterize the DNA binding and organization modes of NAPs and how such information leads to better understanding of its roles in DNA packaging of bacterial nucleoid and transcription regulation.

Large-Scale Conformational Transitions in Supercoiled DNA Revealed by Coarse-Grained Simulation

Topological constraints, such as those associated with DNA supercoiling, play an integral role in genomic regulation and organization in living systems. However, physical understanding of the principles that underlie DNA organization at biologically relevant length scales remains a formidable challenge. We develop a coarse-grained simulation approach for predicting equilibrium conformations of supercoiled DNA. Our methodology enables the study of supercoiled DNA molecules at greater length scales and supercoiling densities than previously explored by simulation. With this approach, we study the conformational transitions that arise due to supercoiling across the full range of supercoiling densities that are commonly explored by living systems. Simulations of ring DNA molecules with lengths at the scale of topological domains in the Escherichia coli chromosome (∼10 kilobases) reveal large-scale conformational transitions elicited by supercoiling. The conformational transitions result in three supercoiling conformational regimes that are governed by a competition among chiral coils, extended plectonemes, and branched hyper-supercoils. These results capture the nonmonotonic relationship of size versus degree of supercoiling observed in experimental sedimentation studies of supercoiled DNA, and our results provide a physical explanation of the conformational transitions underlying this behavior. The length scales and supercoiling regimes investigated here coincide with those relevant to transcription-coupled remodeling of supercoiled topological domains, and we discuss possible implications of these findings in terms of the interplay between transcription and topology in bacterial chromosome organization.

Twist-stretch profiles of DNA chains

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

Twist-Bend Coupling and the Torsional Response of Double-Stranded DNA

Recent magnetic tweezers experiments have reported systematic deviations of the twist response of double-stranded DNA from the predictions of the twistable wormlike chain model. Here we show, by means of analytical results and computer simulations, that these discrepancies can be resolved if a coupling between twist and bend is introduced. We obtain an estimate of 40±10  nm for the twist-bend coupling constant. Our simulations are in good agreement with high-resolution, magnetic-tweezers torque data. Although the existence of twist-bend coupling was predicted long ago [J. Marko and E. Siggia, Macromolecules 27, 981 (1994)MAMOBX0024-929710.1021/ma00082a015], its effects on the mechanical properties of DNA have been so far largely unexplored. We expect that this coupling plays an important role in several aspects of DNA statics and dynamics.

Large-scale chromosome folding versus genomic DNA sequences: A discrete double Fourier transform technique

Using state-of-the-art techniques combining imaging methods and high-throughput genomic mapping tools leaded to the significant progress in detailing chromosome architecture of various organisms. However, a gap still remains between the rapidly growing structural data on the chromosome folding and the large-scale genome organization. Could a part of information on the chromosome folding be obtained directly from underlying genomic DNA sequences abundantly stored in the databanks? To answer this question, we developed an original discrete double Fourier transform (DDFT). DDFT serves for the detection of large-scale genome regularities associated with domains/units at the different levels of hierarchical chromosome folding. The method is versatile and can be applied to both genomic DNA sequences and corresponding physico-chemical parameters such as base-pairing free energy. The latter characteristic is closely related to the replication and transcription and can also be used for the assessment of temperature or supercoiling effects on the chromosome folding. We tested the method on the genome of E. coli K-12 and found good correspondence with the annotated domains/units established experimentally. As a brief illustration of further abilities of DDFT, the study of large-scale genome organization for bacteriophage PHIX174 and bacterium Caulobacter crescentus was also added. The combined experimental, modeling, and bioinformatic DDFT analysis should yield more complete knowledge on the chromosome architecture and genome organization.

Facilitated dissociation of transcription factors from single DNA binding sites

The binding of transcription factors (TFs) to DNA controls most aspects of cellular function, making the understanding of their binding kinetics imperative. The standard description of bimolecular interactions posits that TF off rates are independent of TF concentration in solution. However, recent observations have revealed that proteins in solution can accelerate the dissociation of DNA-bound proteins. To study the molecular basis of facilitated dissociation (FD), we have used single-molecule imaging to measure dissociation kinetics of Fis, a key Escherichia coli TF and major bacterial nucleoid protein, from single dsDNA binding sites. We observe a strong FD effect characterized by an exchange rate [Formula: see text], establishing that FD of Fis occurs at the single-binding site level, and we find that the off rate saturates at large Fis concentrations in solution. Although spontaneous (i.e., competitor-free) dissociation shows a strong salt dependence, we find that FD depends only weakly on salt. These results are quantitatively explained by a model in which partially dissociated bound proteins are susceptible to invasion by competitor proteins in solution. We also report FD of NHP6A, a yeast TF with structure that differs significantly from Fis. We further perform molecular dynamics simulations, which indicate that FD can occur for molecules that interact far more weakly than those that we have studied. Taken together, our results indicate that FD is a general mechanism assisting in the local removal of TFs from their binding sites and does not necessarily require cooperativity, clustering, or binding site overlap.

Theoretical Methods for Studying DNA Structural Transitions under Applied Mechanical Constraints

Recent progress in single-molecule manipulation technologies has made it possible to exert force and torque on individual DNA biopolymers to probe their mechanical stability and interaction with various DNA-binding proteins. It was revealed in these experiments that the DNA structure and formation of nucleoprotein complexes by DNA-architectural proteins can be strongly modulated by an intricate interplay between the entropic elasticity of DNA and its global topology, which is closely related to the mechanical constraints applied to the DNA. Detailed understanding of the physical processes underlying the DNA behavior observed in single-molecule experiments requires the development of a general theoretical framework, which turned out to be a rather challenging task. Here, we review recent advances in theoretical methods that can be used to interpret single-molecule manipulation experiments on DNA.

Dynamic organization of mitotic chromosomes

The assembly of rod-shaped chromosomes during mitosis is an essential prerequisite for faithful segregation of genetic information into daughter cells. Despite the long history of chromosome research, it is only recently that we have acquired powerful approaches and crucial tools that help to unlock the secret of this seemingly complex process. In particular, in vitro assays, mammalian genetics, Hi-C analyses and computer simulations have provided valuable information during the past two years. These studies are now beginning to elucidate how the core components of mitotic chromosomes, namely, histones, topoisomerase IIα and condensins, cooperate with each other to convert very long stretches of DNA into rod-shaped chromosomes.

Molecular Mechanism of Facilitated Dissociation of Fis Protein from DNA

Fis protein is a nucleoid-associated protein that plays many roles in transcriptional regulation and DNA site-specific recombination. In contrast to the naïve expectation based on stoichiometry, recent single-molecule studies have shown that the dissociation of Fis protein from DNA is accelerated by increasing the concentration of the Fis protein. Because the detailed molecular mechanism of facilitated dissociation is still not clear, in this study, we employ computational methods to explore the binding landscapes of Fis:DNA complexes with various stoichiometries. When two Fis molecules are present, simulations uncover a ternary complex, where the originally bound Fis protein is partially dissociated from DNA. The simulations support a three-state sequential kinetic model (N ⇄ I → D) for facilitated dissociation, thus explaining the concentration-dependent dissociation.

Transfer-matrix calculations of DNA polymer micromechanics under tension and torque constraints

Recent development of single-molecule manipulation technologies has made it possible to exert constant force and torque on individual DNA biopolymers to probe their elastic characteristics and structural stability. It has been previously shown that depending on the nature of applied mechanical constraints, DNA can exist in several forms including B-, L-, and P-DNA. However, there is still a lack of understanding of how structural heterogeneity of DNA, which may naturally arise due to sequence-dependent DNA properties, protein binding, or DNA damage, influences local stability of the above DNA states. To provide a more complete and detailed description of the DNA mechanics, we developed a theoretical framework based on transfer-matrix calculations and demonstrated how it can be used to predict the DNA behavior upon application of a wide range of force and torque constraints. The resulting phase diagram shows DNA structural transitions that are in good agreement with previous experimental and theoretical studies. We further discuss how the constructed formalism can be extended to include local inhomogeneities in the DNA physical properties, thus making it possible to investigate the effect of DNA sequence as well as protein binding on DNA structural stability.

Single-molecule fluorescence studies on DNA looping

Structure and dynamics of DNA impact how the genetic code is processed and maintained. In addition to its biological importance, DNA has been utilized as building blocks of various nanomachines and nanostructures. Thus, understanding the physical properties of DNA is of fundamental importance to basic sciences and engineering applications. DNA can undergo various physical changes. Among them, DNA looping is unique in that it can bring two distal sites together, and thus can be used to mediate interactions over long distances. In this paper, we introduce a FRET-based experimental tool to study DNA looping at the single molecule level. We explain the connection between experimental measurables and a theoretical concept known as the J factor with the intent of raising awareness of subtle theoretical details that should be considered when drawing conclusions. We also explore DNA looping-assisted protein diffusion mechanism called intersegmental transfer using protein induced fluorescence enhancement (PIFE). We present some preliminary results and future outlooks.

Condensin-Based Chromosome Organization from Bacteria to Vertebrates

Condensins are large protein complexes that play a central role in chromosome organization and segregation in the three domains of life. They display highly characteristic, rod-shaped structures with SMC (structural maintenance of chromosomes) ATPases as their core subunits and organize large-scale chromosome structure through active mechanisms. Most eukaryotic species have two distinct condensin complexes whose balanced usage is adapted flexibly to different organisms and cell types. Studies of bacterial condensins provide deep insights into the fundamental mechanisms of chromosome segregation. This Review surveys both conserved features and rich variations of condensin-based chromosome organization and discusses their evolutionary implications.

DNA-Segment-Facilitated Dissociation of Fis and NHP6A from DNA Detected via Single-Molecule Mechanical Response

The rate of dissociation of a DNA-protein complex is often considered to be a property of that complex, without dependence on other nearby molecules in solution. We study the kinetics of dissociation of the abundant Escherichia coli nucleoid protein Fis from DNA, using a single-molecule mechanics assay. The rate of Fis dissociation from DNA is strongly dependent on the solution concentration of DNA. The off-rate (k(off)) of Fis from DNA shows an initially linear dependence on solution DNA concentration, characterized by an exchange rate of k(ex)≈9×10(-4) (ng/μl)(-1) s(-1) for 100 mM univalent salt buffer, with a very small off-rate at zero DNA concentration. The off-rate saturates at approximately k(off,max)≈8×10(-3) s(-1) for DNA concentrations above ≈20 ng/μl. This exchange reaction depends mainly on DNA concentration with little dependence on the length of the DNA molecules in solution or on binding affinity, but this does increase with increasing salt concentration. We also show data for the yeast HMGB protein NHP6A showing a similar DNA-concentration-dependent dissociation effect, with faster rates suggesting generally weaker DNA binding by NHP6A relative to Fis. Our results are well described by a model with an intermediate partially dissociated state where the protein is susceptible to being captured by a second DNA segment, in the manner of "direct transfer" reactions studied for other DNA-binding proteins. This type of dissociation pathway may be important to protein-DNA binding kinetics in vivo where DNA concentrations are large.