High-throughput single-molecule experiments reveal heterogeneity, state switching, and three interconnected pause states in transcription

Pausing by bacterial RNA polymerase (RNAp) is vital in the recruitment of regulatory factors, RNA folding, and coupled translation. While backtracking and intra-structural isomerization have been proposed to trigger pausing, our mechanistic understanding of backtrack-associated pauses and catalytic recovery remains incomplete. Using high-throughput magnetic tweezers, we examine the Escherichia coli RNAp transcription dynamics over a wide range of forces and NTP concentrations. Dwell-time analysis and stochastic modeling identify, in addition to a short-lived elemental pause, two distinct long-lived backtrack pause states differing in recovery rates. We identify two stochastic sources of transcription heterogeneity: alterations in short-pause frequency that underlies elongation-rate switching, and variations in RNA cleavage rates in long-lived backtrack states. Together with effects of force and Gre factors, we demonstrate that recovery from deep backtracks is governed by intrinsic RNA cleavage rather than diffusional Brownian dynamics. We introduce a consensus mechanistic model that unifies our findings with prior models.

A kinetic model predicts SpCas9 activity, improves off-target classification, and reveals the physical basis of targeting fidelity

The S. pyogenes (Sp) Cas9 endonuclease is an important gene-editing tool. SpCas9 is directed to target sites based on complementarity to a complexed single-guide RNA (sgRNA). However, SpCas9-sgRNA also binds and cleaves genomic off-targets with only partial complementarity. To date, we lack the ability to predict cleavage and binding activity quantitatively, and rely on binary classification schemes to identify strong off-targets. We report a quantitative kinetic model that captures the SpCas9-mediated strand-replacement reaction in free-energy terms. The model predicts binding and cleavage activity as a function of time, target, and experimental conditions. Trained and validated on high-throughput bulk-biochemical data, our model predicts the intermediate R-loop state recently observed in single-molecule experiments, as well as the associated conversion rates. Finally, we show that our quantitative activity predictor can be reduced to a binary off-target classifier that outperforms the established state-of-the-art. Our approach is extensible, and can characterize any CRISPR-Cas nuclease - benchmarking natural and future high-fidelity variants against SpCas9; elucidating determinants of CRISPR fidelity; and revealing pathways to increased specificity and efficiency in engineered systems.

Hybridization Kinetics Explains CRISPR-Cas Off-Targeting Rules

Due to their specificity, efficiency, and ease of programming, CRISPR-associated nucleases are popular tools for genome editing. On the genomic scale, these nucleases still show considerable off-target activity though, posing a serious obstacle to the development of therapies. Off targeting is often minimized by choosing especially high-specificity guide sequences, based on algorithms that codify empirically determined off-targeting rules. A lack of mechanistic understanding of these rules has so far necessitated their ad hoc implementation, likely contributing to the limited precision of present algorithms. To understand the targeting rules, we kinetically model the physics of guide-target hybrid formation. Using only four parameters, our model elucidates the kinetic origin of the experimentally observed off-targeting rules, thereby rationalizing the results from both binding and cleavage assays. We favorably compare our model to published data from CRISPR-Cas9, CRISPR-Cpf1, CRISPR-Cascade, as well as the human Argonaute 2 system.

Global DNA Compaction in Stationary-Phase Bacteria Does Not Affect Transcription

In stationary-phase Escherichia coli, Dps (DNA-binding protein from starved cells) is the most abundant protein component of the nucleoid. Dps compacts DNA into a dense complex and protects it from damage. Dps has also been proposed to act as a global regulator of transcription. Here, we directly examine the impact of Dps-induced compaction of DNA on the activity of RNA polymerase (RNAP). Strikingly, deleting the dps gene decompacted the nucleoid but did not significantly alter the transcriptome and only mildly altered the proteome during stationary phase. Complementary in vitro assays demonstrated that Dps blocks restriction endonucleases but not RNAP from binding DNA. Single-molecule assays demonstrated that Dps dynamically condenses DNA around elongating RNAP without impeding its progress. We conclude that Dps forms a dynamic structure that excludes some DNA-binding proteins yet allows RNAP free access to the buried genes, a behavior characteristic of phase-separated organelles.

Designing nucleosomal force sensors

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

Extreme bendability of DNA double helix due to bending asymmetry

Experimental data of the DNA cyclization (J-factor) at short length scales exceed the theoretical expectation based on the wormlike chain (WLC) model by several orders of magnitude. Here, we propose that asymmetric bending rigidity of the double helix in the groove direction can be responsible for extreme bendability of DNA at short length scales and it also facilitates DNA loop formation at these lengths. To account for the bending asymmetry, we consider the asymmetric elastic rod (AER) model which has been introduced and parametrized in an earlier study [B. Eslami-Mossallam and M. R. Ejtehadi, Phys. Rev. E 80, 011919 (2009)]. Exploiting a coarse grained representation of the DNA molecule at base pair (bp) level and using the Monte Carlo simulation method in combination with the umbrella sampling technique, we calculate the loop formation probability of DNA in the AER model. We show that the DNA molecule has a larger J-factor compared to the WLC model which is in excellent agreement with recent experimental data.

Nucleosome dynamics: Sequence matters

About three quarter of all eukaryotic DNA is wrapped around protein cylinders, forming nucleosomes. Even though the histone proteins that make up the core of nucleosomes are highly conserved in evolution, nucleosomes can be very different from each other due to posttranslational modifications of the histones. Another crucial factor in making nucleosomes unique has so far been underappreciated: the sequence of their DNA. This review provides an overview of the experimental and theoretical progress that increasingly points to the importance of the nucleosomal base pair sequence. Specifically, we discuss the role of the underlying base pair sequence in nucleosome positioning, sliding, breathing, force-induced unwrapping, dissociation and partial assembly and also how the sequence can influence higher-order structures. A new view emerges: the physical properties of nucleosomes, especially their dynamical properties, are determined to a large extent by the mechanical properties of their DNA, which in turn depends on DNA sequence.

Why Do Nucleosomes Unwrap Asymmetrically?

Nucleosomes, DNA spools with a protein core, engage about three-quarters of eukaryotic DNA and play a critical role in chromosomal processes, ranging from gene regulation, recombination, and replication to chromosome condensation. For more than a decade, micromanipulation experiments where nucleosomes are put under tension, as well as the theoretical interpretations of these experiments, have deepened our understanding of the stability and dynamics of nucleosomes. Here we give a theoretical explanation for a surprising new experimental finding: nucleosomes wrapped onto the 601 positioning sequence (the sequence used in most laboratories) respond highly asymmetrically to external forces by always unwrapping from the same end. Using a computational nucleosome model, we show that this asymmetry can be explained by differences in the DNA mechanics of two very short stretches on the wrapped DNA portion. Our finding suggests that the physical properties of nucleosomes, here the response to forces, can be tuned locally by the choice of the underlying base-pair sequence. This leads to a new view of nucleosomes: a physically highly varied set of DNA-protein complexes whose properties can be tuned on evolutionary time scales to their specific function in the genomic context.

Confinement dynamics of a semiflexible chain inside nano-spheres

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.

Definition of the persistence length in the coarse-grained models of DNA elasticity

By considering the detailed structure of DNA in the base pair level, two possible definitions of the persistence length are compared. One definition is related to the orientation of the terminal base pairs, and the other is based on the vectors which connect two adjacent base pairs at each end of the molecule. It is shown that although these definitions approach each other for long DNA molecules, they are dramatically different on short length scales. We show analytically that the difference mostly comes from the shear flexibility of the molecule and can be used to measure the shear modulus of DNA.

Contribution of nonlocal interactions to DNA elasticity

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.

Asymmetric elastic rod model for DNA

In this paper we consider the anharmonic corrections to the anisotropic elastic rod model for DNA. Our model accounts for the difference between the bending energies of positive and negative rolls, which comes from the asymmetric structure of the DNA molecule. We will show that the model can explain the high flexibility of DNA at small length scales, as well as kink formation at high deformation limit.