A variant-dependent molecular clock with anomalous diffusion models SARS-CoV-2 evolution in humans.

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans has been monitored at an unprecedented level due to the public health crisis, yet the stochastic dynamics underlying such a process is dubious. Here, considering the number of acquired mutations as the displacement of the viral particle from the origin, we performed biostatistical analyses from numerous whole genome sequences on the basis of a time-dependent probabilistic mathematical model. We showed that a model with a constant variant-dependent evolution rate and nonlinear mutational variance with time (i.e., anomalous diffusion) explained the SARS-CoV-2 evolutionary motion in humans during the first 120 wk of the pandemic in the United Kingdom. In particular, we found subdiffusion patterns for the Primal, Alpha, and Omicron variants but a weak superdiffusion pattern for the Delta variant. Our findings indicate that non-Brownian evolutionary motions occur in nature, thereby providing insight for viral phylodynamics.

Gene regulation by a protein translation factor at the single-cell level.

Gene expression is inherently stochastic and pervasively regulated. While substantial work combining theory and experiments has been carried out to study how noise propagates through transcriptional regulations, the stochastic behavior of genes regulated at the level of translation is poorly understood. Here, we engineered a synthetic genetic system in which a target gene is down-regulated by a protein translation factor, which in turn is regulated transcriptionally. By monitoring both the expression of the regulator and the regulated gene at the single-cell level, we quantified the stochasticity of the system. We found that with a protein translation factor a tight repression can be achieved in single cells, noise propagation from gene to gene is buffered, and the regulated gene is sensitive in a nonlinear way to global perturbations in translation. A suitable mathematical model was instrumental to predict the transfer functions of the system. We also showed that a Gamma distribution parameterized with mesoscopic parameters, such as the mean expression and coefficient of variation, provides a deep analytical explanation about the system, displaying enough versatility to capture the cell-to-cell variability in genes regulated both transcriptionally and translationally. Overall, these results contribute to enlarge our understanding on stochastic gene expression, at the same time they provide design principles for synthetic biology.

Repurposing the mammalian RNA-binding protein Musashi-1 as an allosteric translation repressor in bacteria.

The RNA recognition motif (RRM) is the most common RNA-binding protein domain identified in nature. However, RRM-containing proteins are only prevalent in eukaryotic phyla, in which they play central regulatory roles. Here, we engineered an orthogonal post-transcriptional control system of gene expression in the bacterium Escherichia coli with the mammalian RNA-binding protein Musashi-1, which is a stem cell marker with neurodevelopmental role that contains two canonical RRMs. In the circuit, Musashi-1 is regulated transcriptionally and works as an allosteric translation repressor thanks to a specific interaction with the N-terminal coding region of a messenger RNA and its structural plasticity to respond to fatty acids. We fully characterized the genetic system at the population and single-cell levels showing a significant fold change in reporter expression, and the underlying molecular mechanism by assessing the in vitro binding kinetics and in vivo functionality of a series of RNA mutants. The dynamic response of the system was well recapitulated by a bottom-up mathematical model. Moreover, we applied the post-transcriptional mechanism engineered with Musashi-1 to specifically regulate a gene within an operon, implement combinatorial regulation, and reduce protein expression noise. This work illustrates how RRM-based regulation can be adapted to simple organisms, thereby adding a new regulatory layer in prokaryotes for translation control.

Nonequilibrium thermodynamics of the RNA-RNA interaction underlying a genetic transposition program.

Thermodynamic descriptions are powerful tools to formally study complex gene expression programs evolved in living cells on the basis of macromolecular interactions. While transcriptional regulations are often modeled in the equilibrium, other interactions that occur in the cell follow a more complex pattern. Here, we adopt a nonequilibrium thermodynamic scheme to explain the RNA-RNA interaction underlying IS10 transposition. We determine the energy landscape associated with such an interaction at the base-pair resolution, and we present an original scaling law for expression prediction that depends on different free energies characterizing that landscape. Then, we show that massive experimental data of the IS10 RNA-controlled expression are better explained by this thermodynamic description in nonequilibrium. Overall, these results contribute to better comprehend the kinetics of post-transcriptional regulations and, more broadly, the functional consequences of processes out of the equilibrium in biology.

CRISPR-Mediated Strand Displacement Logic Circuits with Toehold-Free DNA.

DNA nanotechnology, and DNA computing in particular, has grown extensively over the past decade to end with a variety of functional stable structures and dynamic circuits. However, the use as designer elements of regular DNA pieces, perfectly complementary double strands, has remained elusive. Here, we report the exploitation of CRISPR-Cas systems to engineer logic circuits based on isothermal strand displacement that perform with toehold-free double-stranded DNA. We designed and implemented molecular converters for signal detection and amplification, showing good interoperability between enzymatic and nonenzymatic processes. Overall, these results contribute to enlarge the repertoire of substrates and reactions (hardware) for DNA computing.