arxiv: 2604.25701 · v1 · submitted 2026-04-28 · ⚛️ physics.bio-ph · physics.data-an· q-bio.BM· q-bio.MN· q-bio.PE

Recognition: unknown

Bayesian Rate Inference for Sequence Motif Dynamics in Systems of Reactive Nucleic Acids

Johannes Harth-Kitzerow , Ulrich Gerland , Torsten A. En{\ss}lin

Authors on Pith no claims yet

Pith reviewed 2026-05-07 14:03 UTC · model grok-4.3

classification ⚛️ physics.bio-ph physics.data-anq-bio.BMq-bio.MNq-bio.PE

keywords Bayesian inferencemotif rate equationsnucleic acid reactionsligation countsRNA worldstrand reactor simulationssequence motifsreaction kinetics

0 comments

The pith

A Bayesian framework infers parameters of motif rate equations from ligation count data in nucleic acid simulations.

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

The paper develops a Bayesian inference method to determine the parameters of motif rate equations that simplify nucleic acid reaction dynamics. These equations project full strand interactions such as hybridization, dehybridization, templated ligation, and cleavage onto sequence motif space. Using ligation count data from detailed strand reactor simulations, the framework calibrates the reduced model and supplies uncertainty estimates. This matters because it supports efficient exploration of the high-dimensional parameter spaces relevant to the RNA world hypothesis, where short RNA strands might build toward catalytic systems. The method also serves as a bridge for eventually inferring rates directly from experimental observations.

Core claim

The authors present a Bayesian inference framework that infers the parameters of motif rate equations from ligation count data generated by strand reactor simulations. The motif rate equations reduce complex nucleic acid strand dynamics to interactions among sequence motifs. This matching procedure between the simpler description and the full simulations includes rigorous uncertainty quantification and acts as a step toward inferring reaction rate constants from laboratory data.

What carries the argument

Bayesian inference framework that extracts motif rate parameters from simulated ligation count data

If this is right

Motif rate equations can be calibrated to reproduce key statistics from detailed strand reactor simulations.
Parameter uncertainties are quantified rigorously through the Bayesian posterior.
High-dimensional parameter spaces for RNA world models become feasible to scan efficiently.
The approach supplies a route to infer reaction rates from experimental ligation data with uncertainty estimates.

Where Pith is reading between the lines

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

The framework could highlight which sequence motifs most strongly influence replication dynamics in prebiotic conditions.
Applying the method across varied environmental parameters in simulations might identify conditions that favor longer strand formation.
Extensions could incorporate additional reaction channels such as cleavage or template-specific effects while retaining the inference procedure.

Load-bearing premise

The motif rate equations must capture enough of the essential strand reactor dynamics for ligation count data to meaningfully constrain the parameters.

What would settle it

Running the inferred motif rates through the rate equations and finding that the predicted ligation counts deviate systematically from the original simulation data outside the estimated uncertainties would falsify the framework.

Figures

Figures reproduced from arXiv: 2604.25701 by Johannes Harth-Kitzerow, Torsten A. En{\ss}lin, Ulrich Gerland.

**Figure 1.** Figure 1: Simulation model of the strand reactor and its projection onto motif rate equation. view at source ↗

**Figure 2.** Figure 2: Inference model: Prior parameters (dashed circle), inferred parameters (solid circles), view at source ↗

**Figure 3.** Figure 3: (a) Ligation counts mean and (b) standard deviation of parameter set 1 that serve as view at source ↗

**Figure 4.** Figure 4: (a) Backprojected ligation counts, (b) prior mean rate constants, (c) theoretical mean view at source ↗

**Figure 5.** Figure 5: (a) Prior and (b) posterior rate constants standard deviations in scenario 1. view at source ↗

**Figure 6.** Figure 6: Sequence motif concentrations of the original simulation (gold) and the integrated view at source ↗

**Figure 7.** Figure 7: (a) Backprojected ligation counts, (b) theoretical mean rate constants, (c) posterior view at source ↗

**Figure 8.** Figure 8: Zebraness of the strand reactor simulation (gold) and the integrated motif rate equations view at source ↗

read the original abstract

The RNA world hypothesis suggests a pathway of how life emerged on early earth. It assumes that life started with RNA based systems, capable of storing, transmitting and replicating information, envisioning that monomers and short RNA oligomers interact to form longer strands, eventually becoming catalytically active ribozymes. Key reactions in RNA pools are hybridization, dehybridization, templated ligation, and cleavage. Those reactions depend on many environmental parameters and the wide range of possible configurations among interacting strands. In order to scan such high dimensional parameter spaces, efficient descriptions are needed. Motif rate equations project complex strand reactor dynamics onto sequence motif space. Here we present a Bayesian inference framework to infer their parameters from ligation count data produced by strand reactor simulations. This provides a framework to match the simpler motif rate equations to more complex simulations. Additionally, it is a step towards inferring reaction rate constants directly from experimental data, including rigorous uncertainty estimation. This could be an essential procedure to connect theory and experiment, and deepen our understanding of the essential features necessary for life to emerge.

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.

Desk Editor's Note private letter to a colleague

This paper gives a Bayesian method to infer motif rate parameters from ligation counts in RNA strand simulations.

read the letter

This paper gives a Bayesian method to infer motif rate parameters from ligation counts in RNA strand simulations. That's the core contribution. It does well by setting up a clear inference pipeline that takes data from the detailed strand reactor runs and fits the simpler motif equations. The use of ligation counts as the key observable is a reasonable choice for this system, and the Bayesian approach brings in uncertainty quantification, which is valuable for high-dimensional models. It also positions the work as a step toward experimental data fitting. The new element is applying this inference specifically to match motif dynamics to full simulations in the context of reactive nucleic acids and the RNA world hypothesis. The workflow seems practical for scanning parameter spaces efficiently. Soft spots include limited discussion of how accurately the motif model captures the full dynamics after fitting. The paper would benefit from more explicit checks, such as forward simulations with the inferred parameters to compare against held-out data or direct metrics of approximation error. The assumption that ligation data constrains the parameters sufficiently could be tested more thoroughly, though no circularity or inconsistency shows up. This is for modelers in origins-of-life chemistry and systems biology who need reduced descriptions of complex reaction networks. Readers working on Bayesian methods for biophysical simulations would also find the setup relevant. It deserves a serious referee because it tackles a practical issue with a statistically grounded method and has clear potential applications. I would recommend sending it for peer review.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a Bayesian inference framework to extract rate parameters for sequence motif dynamics from ligation count data generated by strand-reactor simulations of reactive nucleic acids. The approach is motivated by the RNA-world hypothesis and aims to project complex hybridization, dehybridization, templated ligation, and cleavage dynamics onto a lower-dimensional motif space, thereby enabling efficient parameter scanning and providing a route toward inference from experimental data with uncertainty quantification.

Significance. If the inference procedure is shown to recover parameters accurately and the motif projection remains faithful across the tested regimes, the work supplies a practical bridge between detailed stochastic simulations and reduced-order rate equations. The explicit use of Bayesian methods for uncertainty estimation is a clear methodological strength that could support falsifiable predictions and experimental calibration in prebiotic chemistry.

major comments (2)

[Results] The central claim that ligation counts from strand-reactor simulations contain sufficient information to constrain motif parameters rests on the projection step; without an explicit demonstration (e.g., in the results section) that the inferred rates reproduce the original simulation statistics within posterior predictive checks, the sufficiency of the motif equations cannot be verified.
[Methods] The likelihood function that maps simulated ligation counts to the motif-rate parameters is not specified in sufficient detail (e.g., Poisson vs. binomial assumption, handling of multiple ligation events per motif) to assess whether the posterior is properly calibrated or whether identifiability issues arise for certain motif classes.

minor comments (2)

[Introduction] Notation for motif states and rate constants should be introduced once in a dedicated table or equation block rather than scattered across the text.
[Discussion] The abstract states that the framework is 'a step towards inferring reaction rate constants directly from experimental data'; a short paragraph in the discussion clarifying the additional experimental observables required would strengthen this claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment below and have revised the manuscript to incorporate the suggested improvements where they strengthen the presentation and validation of our Bayesian inference framework.

read point-by-point responses

Referee: [Results] The central claim that ligation counts from strand-reactor simulations contain sufficient information to constrain motif parameters rests on the projection step; without an explicit demonstration (e.g., in the results section) that the inferred rates reproduce the original simulation statistics within posterior predictive checks, the sufficiency of the motif equations cannot be verified.

Authors: We agree that an explicit posterior predictive check would provide stronger evidence for the sufficiency of the motif projection. The original manuscript included parameter recovery tests and comparisons of inferred rates to simulation outputs, but did not contain a dedicated posterior predictive analysis. In the revised version we have added a new subsection to the Results section that performs these checks: posterior samples of the motif rates are used to forward-simulate ligation count distributions via the motif rate equations, which are then compared directly to the strand-reactor data. The added figures demonstrate that the observed ligation statistics lie within the posterior predictive intervals for the tested regimes, thereby confirming that the ligation counts contain sufficient information to constrain the motif parameters. revision: yes
Referee: [Methods] The likelihood function that maps simulated ligation counts to the motif-rate parameters is not specified in sufficient detail (e.g., Poisson vs. binomial assumption, handling of multiple ligation events per motif) to assess whether the posterior is properly calibrated or whether identifiability issues arise for certain motif classes.

Authors: We acknowledge that the likelihood specification was insufficiently detailed. In the revised Methods section we now explicitly state that ligation counts are modeled as independent Poisson random variables whose means are given by the time-integrated motif rate equations. This choice is motivated by the count-valued nature of the data and the low per-motif event probabilities in the simulated regimes; we also note that the Poisson assumption approximates a binomial process when the number of possible ligation sites is large. We have added a short paragraph discussing identifiability, showing that the motif parameters remain identifiable for the classes examined (supported by the shape of the posterior and the condition number of the observed information matrix). These clarifications allow readers to assess calibration and potential identifiability limitations. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper introduces a Bayesian inference procedure to extract motif rate parameters from ligation count data generated by separate strand-reactor simulations. The central claim is that the motif rate equations can be matched to these independent simulation outputs via standard Bayesian updating, providing a bridge between simplified and complex models. No load-bearing step reduces by construction to its own inputs: the motif projection is presented as an explicit modeling choice, the likelihood is constructed from simulation observables, and no parameter is fitted to a subset then re-labeled as a prediction of the same quantity. Self-citations, if present, are not invoked to justify uniqueness or to smuggle in an ansatz that would force the result. The derivation remains self-contained against external simulation benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the framework presumably relies on standard Bayesian priors and likelihoods for rate inference without additional postulates.

pith-pipeline@v0.9.0 · 5510 in / 943 out tokens · 60339 ms · 2026-05-07T14:03:47.209976+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

20 extracted references

[1]

Gradient estimators for parameter inference in discrete stochastic kinetic models.arXiv e-prints, 2604:02121, 2026

Ludwig Burger, Annalena Kofler, Lukas Heinrich, and Ulrich Gerland. Gradient estimators for parameter inference in discrete stochastic kinetic models.arXiv e-prints, 2604:02121, 2026

2026
[2]

Replication elongates short DNA, reduces sequence bias and develops trimer structure.Nucleic Acids Research, 52(3):1290–1297, 2023

Adriana Calacca Serr˜ ao, Felix T D¨ anekamp, Zs´ ofia Meggyesi, and Dieter Braun. Replication elongates short DNA, reduces sequence bias and develops trimer structure.Nucleic Acids Research, 52(3):1290–1297, 2023

2023
[3]

Francis H. C. Crick. The origin of the genetic code.Journal of Molecular Biology, 38(3):367– 379, 1968. 16

1968
[4]

Ferr´ e-D’Amar´ e and William G

Adrian R. Ferr´ e-D’Amar´ e and William G. Scott. Small self-cleaving ribozymes.Cold Spring Harbor Perspectives in Biology, 2(10), 2010

2010
[5]

Philipp Frank, Reimar Leike, and Torsten A. Enßlin. Geometric variational inference.Entropy, 23(7), 2021

2021
[6]

Edoardo Gianni, Samantha L. Y. Kwok, Christopher J. K. Wan, Kevin Goeij, Bryce E. Clifton, Enrico S. Colizzi, James Attwater, and Philipp Holliger. A small polymerase ribozyme that can synthesize itself and its complementary strand.Science, 391(6789):1022–1028, 2026

2026
[7]

Rosenberger, Bernhard Altaner, and Ulrich Gerland

Tobias G¨ oppel, Joachim H. Rosenberger, Bernhard Altaner, and Ulrich Gerland. Thermo- dynamic and kinetic sequence selection in enzyme-free polymer self-assembly inside a non- equilibrium RNA reactor.Life, 12(4), 2022

2022
[8]

Enzyme-free replication with two or four bases.Angewandte Chemie International Edition, 57(29):8911–8915, 2018

Elena H¨ anle and Clemens Richert. Enzyme-free replication with two or four bases.Angewandte Chemie International Edition, 57(29):8911–8915, 2018

2018
[9]

Johannes Harth-Kitzerow, Ulrich Gerland, and Torsten A. Enßlin. MoRSAIK: Sequence Motif Reactor Simulation, Analysis and Inference Kit in python.arXiv e-prints, 2512:02204, 2025

2025
[10]

Enßlin, and Ulrich Gerland

Johannes Harth-Kitzerow, Tobias G¨ oppel, Ludwig Burger, Torsten A. Enßlin, and Ulrich Gerland. Sequence motif dynamics in RNA pools.Phys. Rev. E, 113:024407, 2026

2026
[11]

The RNA world: molecular cooperation at the origins of life.Nat Rev Genet, 16(1):7–17, 2014

Paul G Higgs and Niles Lehman. The RNA world: molecular cooperation at the origins of life.Nat Rev Genet, 16(1):7–17, 2014

2014
[12]

Jakob Knollm¨ uller and Torsten A. Enßlin. Encoding prior knowledge in the structure of the likelihood.arXiv e-prints, 1812:04403, 2018

2018
[13]

Jakob Knollm¨ uller and Torsten A. Enßlin. Metric Gaussian Variational Inference.arXiv e-prints, 1901:11033, 2019

1901
[14]

Matthew Levy and Stanley L. Miller. The stability of the RNA bases: Implications for the origin of life.Proceedings of the National Academy of Sciences, 95(14):7933–7938, 1998

1998
[15]

Leslie E. Orgel. Evolution of the genetic apparatus.Journal of Molecular Biology, 38(3):381– 393, 1968

1968
[16]

Rosenberger, Tobias G¨ oppel, Patrick W

Joachim H. Rosenberger, Tobias G¨ oppel, Patrick W. Kudella, Dieter Braun, Ulrich Gerland, and Bernhard Altaner. Self-assembly of informational polymers by templated ligation.Phys. Rev. X, 11:031055, 2021

2021
[17]

Taylor, Vitor B

Alexander I. Taylor, Vitor B. Pinheiro, Matthew J. Smola, Alexey S. Morgunov, Sew Peak- Chew, Christopher Cozens, Kevin M. Weeks, Piet Herdewijn, and Philipp Holliger. Catalysts from synthetic genetic polymers.Nature (London), 518(7539):427–430, 2015. 17

2015
[18]

Tkachenko and Sergei Maslov

Alexei V. Tkachenko and Sergei Maslov. Spontaneous emergence of autocatalytic information- coding polymers.The Journal of Chemical Physics, 143(4):045102, 2015

2015
[19]

Tkachenko and Sergei Maslov

Alexei V. Tkachenko and Sergei Maslov. Onset of natural selection in populations of auto- catalytic heteropolymers.The Journal of Chemical Physics, 149(13):134901, 2018

2018
[20]

Zaug and Thomas R

Arthur J. Zaug and Thomas R. Cech. The intervening sequence RNA ofTetrahymenais an enzyme.Science, 231(4737):470–475, 1986. 18

1986