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.

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.

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.

Multiplexing Genetic and Nucleosome Positioning Codes: A Computational Approach.

Eukaryotic DNA is strongly bent inside fundamental packaging units: the nucleosomes. It is known that their positions are strongly influenced by the mechanical properties of the underlying DNA sequence. Here we discuss the possibility that these mechanical properties and the concomitant nucleosome positions are not just a side product of the given DNA sequence, e.g. that of the genes, but that a mechanical evolution of DNA molecules might have taken place. We first demonstrate the possibility of multiplexing classical and mechanical genetic information using a computational nucleosome model. In a second step we give evidence for genome-wide multiplexing in Saccharomyces cerevisiae and Schizosacharomyces pombe. This suggests that the exact positions of nucleosomes play crucial roles in chromatin function.

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.

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.