pith. sign in

arxiv: 2207.04568 · v1 · submitted 2022-07-10 · 🧬 q-bio.SC · cond-mat.stat-mech· q-bio.BM

Stochastic dynamics and ribosome-RNAP interactions in Transcription-Translation Coupling

Xiangting Li , Tom Chou This is my paper

Pith reviewed 2026-05-24 12:00 UTC · model grok-4.3

classification 🧬 q-bio.SC cond-mat.stat-mechq-bio.BM
keywords transcription-translation couplingstochastic modelribosome-RNAP interactionsdelay time distributionsprotection from terminationRNAP pausingprokaryotic gene expressionelongation rates
0
0 comments X

The pith

A stochastic model of ribosome-RNAP interactions predicts delay time distributions and new quantitative measures for transcription-translation coupling.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a continuous-time stochastic model that incorporates elongation rates, initiation and termination, RNAP pausing, and direct ribosome-RNAP exclusion plus binding interactions. The model derives how these molecular features shape the statistical distributions of delay times between transcription and translation. It further defines two new TTC metrics, the direct binding probability and the fraction of time protected from termination proteins, which vary differently with the underlying rates. These allow the identification of parameter regimes that produce acceleration or deceleration of transcription along with different levels of protection.

Core claim

By treating ribosome and RNAP as interacting stochastic particles with constant-rate binding and exclusion, the model yields explicit delay-time distributions and two additional coupling measures; the binding probability quantifies direct contact while the protected-time fraction quantifies resistance to termination attack, and both metrics together reveal when the coupled system accelerates, decelerates, or shields the processes.

What carries the argument

Continuous-time stochastic model of ribosome and RNAP motion that includes constant-rate exclusion and binding interactions, from which delay distributions and the two TTC metrics are computed.

If this is right

  • Delay time distributions are controlled by the specific values of elongation rates, pausing frequency, and binding parameters.
  • The direct ribosome-RNAP binding probability provides a measure of coupling strength distinct from timing statistics.
  • The protected-time fraction quantifies how much of the transcription process avoids attack by termination proteins.
  • Depending on parameter values the model produces either faster or slower overall transcription.
  • The two new metrics respond differently to changes in the rates of known molecular processes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If binding rates turn out to be sequence-dependent, extending the model to position-specific rates would be a direct next step.
  • Live-cell measurements of delay distributions under controlled crowding could test the constant-rate prediction.
  • The protection metric offers a way to interpret how termination-factor mutations affect coupled versus uncoupled transcription.
  • The same stochastic framework could be adapted to explore analogous coupling in systems where multiple polymerases interact on the same template.

Load-bearing premise

Binding and exclusion rates between ribosome and RNAP remain fixed constants that do not depend on sequence context or cellular crowding.

What would settle it

Single-molecule or population measurements of delay distributions and protected-time fractions that remain unchanged when sequence context or molecular crowding is varied would show the constant-rate assumption fails.

Figures

Figures reproduced from arXiv: 2207.04568 by Tom Chou, Xiangting Li.

Figure 1
Figure 1. Figure 1: Schematic of translation of a nascent mRNA tran [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (A) State space of the stochastic model defined in terms of the leading ribosome and RNAP positions [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Internal state space associated with different values [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of different TTC indices. Common [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Slowdown induced by molecular coupling. The common parameters used are the same as those in Fig. [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Effects of molecular coupling. For those parameters not varied, we use the same values used to generate Figs. [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: TTC and performance under genomic variability. Again, we use the standard set of fixed parameter values as in Figs. [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
read the original abstract

Under certain cellular conditions, transcription and mRNA translation in prokaryotes appear to be "coupled," in which the formation of mRNA transcript and production of its associated protein are temporally correlated. Such transcription-translation coupling (TTC) has been evoked as a mechanism that speeds up the overall process, provides protection during the transcription, and/or regulates the timing of transcript and protein formation. What molecular mechanisms underlie ribosome-RNAP coupling and how they can perform these functions have not been explicitly modeled. We develop and analyze a continuous-time stochastic model that incorporates ribosome and RNAP elongation rates, initiation and termination rates, RNAP pausing, and direct ribosome and RNAP interactions (exclusion and binding). Our model predicts how distributions of delay times depend on these molecular features of transcription and translation. We also propose additional measures for TTC: a direct ribosome-RNAP binding probability and the fraction of time the translation-transcription process is "protected" from attack by transcription-terminating proteins. These metrics quantify different aspects of TTC and differentially depend on parameters of known molecular processes. We use our metrics to reveal how and when our model can exhibit either acceleration or deceleration of transcription, as well as protection from termination. Our detailed mechanistic model provides a basis for designing new experimental assays that can better elucidate the mechanisms of TTC.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a continuous-time stochastic model of transcription-translation coupling (TTC) that incorporates ribosome and RNAP elongation rates, initiation/termination rates, RNAP pausing, and direct ribosome-RNAP interactions via exclusion and binding. The model is used to predict how delay-time distributions depend on these molecular features and to define two new TTC metrics (direct binding probability and the fraction of time the process is protected from transcription-terminating proteins). These metrics are then applied to identify parameter regimes that produce acceleration or deceleration of transcription as well as protection from termination.

Significance. If the internal consistency of the stochastic construction holds, the work supplies a mechanistic, parameter-explicit framework for quantifying distinct aspects of TTC. The explicit treatment of stochastic pausing, binding, and exclusion, together with the two new metrics that differentially depend on known molecular rates, offers a basis for designing targeted experiments. The modeling approach is a clear strength relative to purely phenomenological descriptions.

minor comments (3)
  1. [model description] The abstract states that binding and exclusion rates are treated as constant parameters; the main text should explicitly state the numerical ranges or functional forms adopted for these rates and the justification for treating them as sequence- and crowding-independent (model-description section).
  2. [results figures] Figure legends and/or table captions should indicate the number of stochastic realizations used to generate each delay-time distribution and the convergence criterion applied.
  3. [TTC metrics] The protection metric is defined in terms of time spent in a 'protected' state; the precise state-space definition (which combinations of ribosome-RNAP configurations count as protected) should be stated as an equation or enumerated list.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, including the recognition of its mechanistic framework, explicit treatment of stochastic features, and potential for guiding experiments. The recommendation of minor revision is noted. No major comments were provided in the report, so we have no specific points to address point-by-point.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper constructs an explicit continuous-time stochastic model with stated parameters for elongation, initiation, termination, pausing, exclusion, and binding rates. Delay distributions and TTC metrics (binding probability, protection fraction) are computed directly from the model dynamics and parameter values. No equations reduce by construction to fitted inputs, no self-definitional re-use of outputs as inputs, and no load-bearing self-citations or imported uniqueness theorems appear in the abstract or model description. The derivation is therefore self-contained and independent of its own predictions.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. The model implicitly treats molecular rates as constant inputs and assumes direct binding is possible, but these cannot be audited without the full text.

pith-pipeline@v0.9.0 · 5766 in / 1171 out tokens · 18364 ms · 2026-05-24T12:00:07.992297+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

50 extracted references · 50 canonical work pages

  1. [1]

    Woodgate, D

    Stevenson-Jones, F., J. Woodgate, D. Castro-Roa, and N. Zenkin, 2020. Ribosome reactivates transcription by physically pushing RNA polymerase out of transcription arrest. Proceedings of the National Academy of Sciences 117:8462–8467

    work page 2020

  2. [2]

    Muteeb, N

    Chalissery, J., G. Muteeb, N. C. Kalarickal, S. Mohan, V. Jisha, and R. Sen, 2011. Interaction surface of the transcription terminator Rho required to form a complex with the C-terminal domain of the antiterminator NusG. Journal of Molecular Biology405:49–64

    work page 2011

  3. [3]

    Lawson, M. R., W. Ma, M. J. Bellecourt, I. Artsimovitch, A. Martin, R. Landick, K. Schulten, and J. M. Berger,

    work page

  4. [4]

    Molecular cell71:911–922

    Mechanism for the regulated control of bacterial transcription termination by a universal adaptor protein. Molecular cell71:911–922

    work page

  5. [5]

    Kohler, R., R. A. Mooney, D. J. Mills, R. Landick, and P.Cramer,2017.Architectureofatranscribing-translating expressome. Science 356:194–197

    work page 2017

  6. [6]

    Mobli, X

    Ma, C., M. Mobli, X. Yang, A. N. Keller, G. F. King, and P.J.Lewis,2015. RNApolymerase-inducedremodelling ofNusAproducesapauseenhancementcomplex. Nucleic Acids Research43:2829–2840

    work page 2015

  7. [7]

    Proshkin,S.,A.R.Rahmouni,A.Mironov,andE.Nudler,

    work page

  8. [8]

    Cooperation Between Translating Ribosomes and RNA Polymerase in Transcription Elongation.Science 328:504–508

    work page

  9. [9]

    Iyer, S., D. Le, B. R. Park, and M. Kim, 2018. Distinct mechanisms coordinate transcription and translation un- der carbon and nitrogen starvation in Escherichia coli. Nature microbiology3:741–748

    work page 2018

  10. [10]

    Vogel, U., and K. F. Jensen, 1994. The RNA chain elongationrateinEscherichiacolidependsonthegrowth rate. Journal of bacteriology176:2807–2813

    work page 1994

  11. [11]

    Fan, H., A. B. Conn, P. B. Williams, S. Diggs, J. Hahm, J. Gamper, Howard B., Y.-M. Hou, S. E. O’Leary, Y. Wang, and G. M. Blaha, 2017. Transcrip- tion–translation coupling: direct interactions of RNA polymerase with ribosomes and ribosomal subunits.Nu- cleic Acids Research45:11043–11055

    work page 2017

  12. [12]

    Mooney, R. A., S. E. Davis, J. M. Peters, J. L. Rowland, A.Z.Ansari,andR.Landick,2009. RegulatorTrafficking on Bacterial Transcription Units In Vivo.Molecular Cell 33:97–108

    work page 2009

  13. [13]

    Saxena, S., K. K. Myka, R. Washburn, N. Costantino, D. L. Court, and M. E. Gottesman, 2018. Escherichia coli transcription factor NusG binds to 70S ribosomes. Molecular Microbiology108:495–504

    work page 2018

  14. [14]

    Stitt,M.E.Gottesman,andP.Rösch,2010

    Burmann,B.M.,K.Schweimer,X.Luo,M.C.Wahl,B.L. Stitt,M.E.Gottesman,andP.Rösch,2010. ANusE:NusG Complex Links Transcription and Translation.Science 328:501–504. 12 Manuscript submitted to Biophysical Journal Stochastic ribosome-RNAP coupling model

    work page 2010

  15. [15]

    Chou, 2022

    Zuo, X., and T. Chou, 2022. Density- and elongation speed-dependenterrorcorrectioninRNApolymerization. Physical Biology19:026001

    work page 2022

  16. [16]

    Lakatos, 2004

    Chou, T., and G. Lakatos, 2004. Clustered Bottlenecks in mRNA Translation and Protein Synthesis.Phys. Rev. Lett.93:198101

    work page 2004

  17. [17]

    Fredrick, 2018

    Chen, M., and K. Fredrick, 2018. Measures of single- versus multiple-round translation argue against a mech- anism to ensure coupling of transcription and transla- tion. Proceedings of the National Academy of Sciences 115:10774–10779

    work page 2018

  18. [18]

    Zhu, M., M. Mori, T. Hwa, and X. Dai, 2019. Disruption of transcription–translation coordination in Escherichia colileadstoprematuretranscriptionaltermination. Nature microbiology4:2347–2356

    work page 2019

  19. [19]

    Mäkelä,J.,J.Lloyd-Price,O.Yli-Harja,andA.S.Ribeiro,

    work page

  20. [20]

    Stochastic sequence-level model of coupled tran- scription and translation in prokaryotes.BMC bioinfor- matics 12:1–13

    work page

  21. [21]

    Johnson,G.E.,J.-B.Lalanne,M.L.Peters,andG.-W.Li,

    work page

  22. [22]

    Functionally uncoupled transcription–translation in Bacillus subtilis.Nature585:124–128

    work page

  23. [23]

    Dai, X., M. Zhu, M. Warren, R. Balakrishnan, V. Patsalo, H.Okano,J.R.Williamson,K.Fredrick,Y.-P.Wang,and T.Hwa,2016. Reductionoftranslatingribosomesenables Escherichia coli to maintain elongation rates during slow growth. Nature Microbiology2

    work page 2016

  24. [24]

    Tuller, 2017

    Shaham, G., and T. Tuller, 2017. Genome scale analysis of Escherichia coli with a comprehensive prokaryotic sequence-basedbiophysicalmodeloftranslationinitiation and elongation.DNA Research25:195–205

    work page 2017

  25. [25]

    Riezman, 1977

    Kennell, D., and H. Riezman, 1977. Transcription and translation initiation frequencies of the Escherichia coli lac operon.Journal of Molecular Biology114:1–21

    work page 1977

  26. [26]

    Xu, L., H. Chen, X. Hu, R. Zhang, Z. Zhang, and Z. W. Luo,2006. AverageGeneLengthIsHighlyConservedin Prokaryotes and Eukaryotes and Diverges Only Between the Two Kingdoms.Molecular Biology and Evolution 23:1107–1108

    work page 2006

  27. [27]

    Bremer, 1976

    Young, R., and H. Bremer, 1976. Polypeptide-chain- elongation rate in Escherichia coli B/r as a function of growth rate.The Biochemical journal160:185–94

    work page 1976

  28. [28]

    Realtimedetermi- nationofbacterialinvivoribosometranslationelongation speedbasedonLacZ 𝛼 complementationsystem

    Zhu,M.,X.Dai,andY.-P.Wang,2016. Realtimedetermi- nationofbacterialinvivoribosometranslationelongation speedbasedonLacZ 𝛼 complementationsystem. Nucleic Acids Research44:e155–e155

    work page 2016

  29. [29]

    Nudler, 2003

    Epshtein, V., and E. Nudler, 2003. Cooperation Between RNA Polymerase Molecules in Transcription Elongation. Science 300:801–805

    work page 2003

  30. [30]

    Neuman, K. C., E. A. Abbondanzieri, R. Landick, J. Gelles, and S. M. Block, 2003. Ubiquitous Tran- scriptional Pausing Is Independent of RNA Polymerase Backtracking. Cell 115:437–447

    work page 2003

  31. [31]

    Su, and R

    Wang,C.,V.Molodtsov,E.Firlar,J.T.Kaelber,G.Blaha, M. Su, and R. H. Ebright, 2020. Structural basis of transcription-translation coupling. Science 369:1359– 1365

    work page 2020

  32. [32]

    Zielenkiewicz, 2013

    Siwiak, M., and P. Zielenkiewicz, 2013. Transimulation - Protein Biosynthesis Web Service.PLoS ONE8:e73943

    work page 2013

  33. [33]

    Roberts, 2008

    Hatoum, A., and J. Roberts, 2008. Prevalence of RNA polymerase stalling at Escherichia coli promoters after opencomplexformation. Molecularmicrobiology 68:17– 28

    work page 2008

  34. [34]

    Peeling and sliding in nucleosome repositioning

    Chou, T., 2007. Peeling and sliding in nucleosome repositioning. Physical Review Letters99:058105

    work page 2007

  35. [35]

    Spaulding, and A

    Teimouri, H., C. Spaulding, and A. B. Kolomeisky, 2022. Optimal pathways control fixation of multiple mutations during cancer initiation.Biophysical Journal

    work page 2022

  36. [36]

    Bortz, A. B., M. H. Kalos, and J. L. Lebowitz, 1975. A new algorithm for Monte Carlo simulation of Ising spin systems. Journal of Computational Physics17:10–18

    work page 1975

  37. [37]

    T., 1977

    Gillespie, D. T., 1977. Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry 81:2340–2361

    work page 1977

  38. [38]

    Koslover, D. J., F. M. Fazal, R. A. Mooney, R. Landick, and S. M. Block, 2012. Binding and Translocation of Termination Factor Rho Studied at the Single-Molecule Level. Journal of Molecular Biology423:664–676

    work page 2012

  39. [39]

    Journal of Biological Chemistry272:15329–15338

    Komissarova,N.,andM.Kashlev,1997.RNApolymerase switches between inactivated and activated states by translocatingbackandforthalongtheDNAandtheRNA. Journal of Biological Chemistry272:15329–15338

    work page 1997

  40. [40]

    RNA polymerase active center: the molecular engine of transcription

    Nudler, E., 2009. RNA polymerase active center: the molecular engine of transcription. Annual review of biochemistry78:335–361

    work page 2009

  41. [41]

    Davenport, R. J., G. J. Wuite, R. Landick, and C. Busta- mante, 2000. Single-Molecule Study of Transcriptional Pausing and Arrest byE. coliRNA Polymerase.Science 287:2497–2500

    work page 2000

  42. [42]

    Nayak, C

    Larson,M.H.,R.A.Mooney,J.M.Peters,T.Windgassen, D. Nayak, C. A. Gross, S. M. Block, W. J. Greenleaf, R. Landick, and J. S.Weissman, 2014. A pause sequence enriched at translation start sites drives transcription dynamics in vivo.Science 344:1042–1047. Manuscript submitted to Biophysical Journal13 Li and Chou

    work page 2014

  43. [43]

    Ribosome recycling, diffusion, and mRNA loop formation in translational regulation.Bio- physical Journal85:755–773

    Chou, T., 2003. Ribosome recycling, diffusion, and mRNA loop formation in translational regulation.Bio- physical Journal85:755–773

    work page 2003

  44. [44]

    Klumpp, S., J. J. Dong, and T. Hwa, 2012. On ribosome load, codon bias and protein abundance.PLoS ONE 7:e48542

    work page 2012

  45. [45]

    T., and J

    MacDonald, C. T., and J. H. Gibbs, 1969. Concerning the Kinetics of Polypeptide Synthesis on Polyribosomes. Biopolymers. Biopolymers7:707–725

    work page 1969

  46. [46]

    Duc,K.D.,andY.S.Song,2018.Theimpactofribosomal interference, codon usage, and exit tunnel interactions on translation elongation rate variation.PLOS Genetics 14:e1007166

    work page 2018

  47. [47]

    center-of-mass

    Morisaki, T., K. Lyon, K. F. DeLuca, J. G. DeLuca, B. P. English, Z. Zhang, L. D. Lavis, J. B. Grimm, S. Viswanathan, L. L. Looger, T. Lionnet, and T. J. Sta- sevich, 2016. Real-time quantification of single RNA translation dynamics in living cells.Science352:1425– 1429. 14 Manuscript submitted to Biophysical Journal Stochastic ribosome-RNAP coupling model...

    work page 2016

  48. [48]

    Then, the ribosome is always within close range of the RNAP and the system freely cycles among the four internal macrostates

    The instantaneous speeds satisfy𝑝 𝑞. Then, the ribosome is always within close range of the RNAP and the system freely cycles among the four internal macrostates. We may assume that the binding and unbinding rates𝑘a and 𝑘d are much larger than the pausing and unstalling rates𝑘 and 𝑘¸

    work page

  49. [49]

    This system maintains an appreciable probability of being coupled

    The instantaneous speeds satisfy𝑝 Ÿ 𝑞 and the rate of uncoupling𝑘d is slower than the rate of pausing𝑘. This system maintains an appreciable probability of being coupled. When the RNAP is bound and processive, the distance quickly increases until𝑑 ℓ. Because𝑘 ¡ 𝑘 d, the RNAP pauses often before it can break free from the ribosome. When the internal states...

    work page

  50. [50]

    𝑎 states

    The instantaneous speeds satisfy𝑝 Ÿ 𝑞 , but the dissociation rate𝑘d is larger than the pausing rate𝑘. This scenario is essentially the same as the previous one, with the only difference that the transition from the bound, processive state to an unbound processive state is fast and effectively irreversible. These scenarios can be coarse-grained into different...

    work page