Abstract : Gene expression consists in the synthesis of proteins from the information encoded on DNA. One of the two main steps of gene expression is the translation of messenger RNA (mRNA) into polypeptide sequences of amino acids. Here, by taking into account mRNA degradation, we model the motion of ribosomes along mRNA with a ballistic model where particles advance along a filament without excluded volume interactions. Unidirectional models of transport have previously been used to fit the average density of ribosomes obtained by the experimental ribosome sequencing (Ribo-Seq) technique. In this case an inverse fit gives access to the kinetic rates: the position-dependent speeds and the entry rate of ribosomes into mRNA. The degradation rate is not, however, accounted for and experimental data from different experiments are needed to have enough parameters for the fit. Here, we propose an entirely novel experimental setup and theoretical framework consisting in splitting the mRNAs into categories depending on the number of ribosomes from one to four. We analytically solve the ballistic model for a fixed number of ribosomes per mRNA, study the different regimes of degradation, and propose a criterion for the quality of the inverse fit. The proposed method provides a high sensitivity to the mRNA degradation rate. The additional equations coming from using the monosome (single ribosome) and polysome (arbitrary number) Ribo-Seq profiles enables us to determine all the kinetic rates in terms of the experimentally accessible mRNA degradation rate.

The bacterial flagellar motor (BFM) is the membrane-embedded rotary molecular motor which turns the flagellum that provides thrust to many bacterial species. This large multimeric complex, composed of a few dozen constituent proteins, has emerged as a hallmark of dynamic subunit exchange. The stator units are inner-membrane ion channels which dynamically bind and unbind to the peptidoglycan at the rotor periphery, consuming the ion motive force (IMF) and applying torque to the rotor when bound. The dynamic exchange is known to be a function of the viscous load on the flagellum, allowing the bacterium to dynamically adapt to its local viscous environment, but the molecular mechanisms of exchange and mechanosensitivity remain to be revealed. Here, by actively perturbing the steady-state stator stoichiometry of individual motors, we reveal a stoichiometry-dependent asymmetry in stator remodeling kinetics. We interrogate the potential effect of next-neighbor interactions and local stator unit depletion and find that neither can explain the observed asymmetry. We then simulate and fit two mechanistically diverse models which recapitulate the asymmetry, finding stator assembly dynamics to be particularly well described by a two-state catch-bond mechanism.

erratum

ParABS, the most widespread bacterial DNA segregation system, is composed of a centromeric sequence, parS, and two proteins, the ParA ATPase and the ParB DNA binding proteins. Hundreds of ParB proteins assemble dynamically to form nucleoprotein parS-anchored complexes that serve as substrates for ParA molecules to catalyze positioning and segregation events. The exact nature of this ParBS complex has remained elusive, what we address here by revisiting the Stochastic Binding model (SBM) introduced to explain the non-specific binding profile of ParB in the vicinity of parS. In the SBM, DNA loops stochastically bring loci inside a sharp cluster of ParB. However, previous SBM versions did not include the negative supercoiling of bacterial DNA, leading to use unphysically small DNA persistences to explain the ParB binding profiles. In addition, recent super-resolution microscopy experiments have revealed a ParB cluster that is significantly smaller than previous estimations and suggest that it results from a liquid-liquid like phase separation. Here, by simulating the folding of long (>30 kb) supercoiled DNA molecules calibrated with realistic DNA parameters and by considering different possibilities for the physics of the ParB cluster assembly, we show that the SBM can quantitatively explain the ChIP-seq ParB binding profiles without any fitting parameter, aside from the supercoiling density of DNA, which, remarkably, is in accord with independent measurements. We also predict that ParB assembly results from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., ParB clusters behave like liquid-like protein condensates with unconventional “leaky” boundaries.

Bacterial ParB partitioning proteins involved in chromosomes and low-copy-number plasmid segregation are CTP-dependent molecular switches. CTP-binding converts ParB dimers to DNA clamps, allowing unidimensional diffusion along the DNA. This sliding property has been proposed to explain the ParB spreading over large distances from parS centromere sites where ParB is specifically loaded. We modeled such a ‘Clamping & sliding’ mechanism as a typical reaction-diffusion system, compared it to the F-plasmid ParB DNA binding pattern, and found that it can account neither for the long range of ParB binding to DNA, nor for the rapid assembly kinetics observed in vivo after parS duplication. Also, it predicts a strong effect on the F-plasmid ParB binding pattern from the presence of a roadblock that is not observed in ChIP-seq. We conclude that although ‘Clamping & sliding’ can occur at short distances from parS, another mechanism must apply for ParB recruitment at larger genomic distances.

We investigate the measurement of DNA supercoiling density ($\sigma$) along chromosomes using interaction frequencies between DNA and DNA-anchored clusters of proteins. Specifically, we show how the physics of DNA supercoiling leads, in bacteria, to the quantitative modeling of binding properties of ParB proteins around their centromere-like site, {\it parS}. Using this framework, we provide an upper bound for $\sigma$ in the {\it Escherichia coli} chromosome, consistent with plasmid values, and offer a proof of concept for a high accuracy measurement. To reach these conclusions, we revisit the problem of the formation of ParB clusters. We predict, in particular, that they result from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., they behave like liquid-like protein condensates with unconventional ``leaky'' boundaries.

Translation is one of the main steps in the synthesis of proteins. It consists of ribosomes that translate sequences of nucleotides encoded on mRNA into polypeptide sequences of amino acids. Ribosomes bound to mRNA move unidirectionally, while unbound ribosomes diffuse in the cytoplasm. It has been hypothesized that finite diffusion of ribosomes plays an important role in ribosome recycling and that mRNA circularization enhances the efficiency of translation. In order to estimate the effect of cytoplasmic diffusion on the rate of translation, we consider a Totally Asymmetric Simple Exclusion Process (TASEP) coupled to a finite diffusive reservoir, which we call the Ribosome Transport model with Diffusion (RTD). In this model, we derive an analytical expression for the rate of protein synthesis as a function of the diffusion constant of ribosomes, which is corroborated with results from continuous-time Monte Carlo simulations. Using a wide range of biological relevant parameters, we conclude that diffusion in biological cells is fast enough so that it does not play a role in controlling the rate of translation initiation.