REVIEW 3 major objections 5 minor 35 references
Acceleration of enzymatic reaction-diffusion kinetics by intermediate state
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Intermediate chemical states speed up molecular motors by lowering the effective energy barrier, and the speedup grows under an external load.
desk verdict A clean reaction-diffusion analysis of intermediate states in motors; the finite-w̃ numerics need opening up, but the analytic limits and design principle are worth refereeing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the effective barrier height $\Delta G^\ddagger=k_BT\ln[(e^{\Delta G_p^\ddagger/k_BT}+e^{\Delta G_i^\ddagger/k_BT})/2]$, where $\Delta G_p^\ddagger$ and $\Delta G_i^\ddagger$ are the barrier heights from the primary and intermediate minima to their intersection point. It converts the two-step reaction into a single effective Arrhenius time, and the paper shows that the intermediate position minimizing $\Delta G^\ddagger$ is close to the position maximizing the steady-state flux $J$. The supporting machinery is an overdamped Langevin description of diffusion on harmonic potentials $U_{p,n}$ and $U_{i,n}$; switching rates satisfying local detailed balance with a single asymmetry parameter $q$; the potential of mean force $V(x)=-k_BT\ln\sum_n(e^{-U_{p,n}/k_BT}+e^{-U_{i,n}/k_BT})$, the effective potential seen after chemical equilibration in the fast-switching limit; and closed-form effective rates of the Markov-jump limit, where spatial diffusion equilibrates before each reaction, that explain the slow-switching slowdown.
What would settle it
A direct test is to repeat the calculation with two intermediate states per step instead of one: if the optimized flux still grows under load, the design rule generalizes, but if the acceleration reverses or vanishes, the single-intermediate construction is the actual source of the effect. A complementary experiment would measure $J^*/J_0$ for a motor with a tunable intermediate position in the slow-switching, strongly asymmetric regime, where the paper predicts a slowdown; observing a speedup there would refute the claimed boundary.
Extended reading notes
Core claim
The paper's central claim is that, for a reaction-diffusion motor whose cycle is $p\to i\to p$ with one intermediate state per step, the intermediate state accelerates the motor whenever the bare switching rate $\tilde{w}$ is large enough: the optimized flux $J^*$ exceeds the direct flux $J_0$ for every asymmetry parameter $q\in[0,1]$. The reason is geometric and energetic: the intermediate potential well sits between the primary wells, splitting the spatial displacement $\Delta x$ and the chemical free-energy drop $\Delta G$ in two, so the harmonic barriers $\Delta G_p^\ddagger$ and $\Delta G_i^\ddagger$ are each smaller than the direct barrier. The effective barrier height $\Delta G^\ddagger=k_BT\ln[(e^{\Delta G_p^\ddagger/k_BT}+e^{\Delta G_i^\ddagger/k_BT})/2]$ is introduced, and the numerical results show that the position minimizing $\Delta G^\ddagger$ nearly coincides with the position maximizing the flux. Under an external load, modeled as a constant tilt of the potential, the load distribution factor $\theta$ (the fraction of the load acting on the forward barrier) is less than half for the two-step route, so the direct barrier grows faster with load than the intermediate-state barrier; consequently $J^*/J_0$ rises with the work $W_{\rm ext}$ per cycle, and the maximum Stokes efficiency stays high where the no-intermediate motor's efficiency collapses. The exception is the small-$\tilde{w}$ limit: for $q=0$ and $q=1$ the effective rates reduce to closed forms such as $J^*/J_0=1/(1+e^{\Delta G/2k_BT})$, which is always less than one, so in that limit the intermediate state slows the motor.
Load-bearing premise
The model assumes exactly one intermediate state per cycle, visited in strict order, with harmonic potentials for both primary and intermediate states and a single asymmetry parameter $q$; real motors with multiple intermediates, anharmonic potentials, or different rate splits may not obey the same design rule.
Editorial extensions
If this is right
- With equal spring constants and a large enough bare switching rate, the optimal intermediate position is the midpoint of the step, $\Delta G_i=\Delta G/2$ and $\Delta x_i=\Delta x/2$, and this position accelerates the motor for every $q$ between 0 and 1.
- Softer intermediate potentials (smaller $k_i$) give larger maximum flux and a wider region of intermediate positions that still accelerate the motor, so the acceleration is less sensitive to fine-tuning.
- Under external load the flux ratio $J^*/J_0$ and the maximum Stokes efficiency advantage over the no-intermediate motor both grow, so intermediate states are what let a motor keep working while doing mechanical work.
- In the slow-switching limit the intermediate state is beneficial only for intermediate $q$; for $q=0$ and $q=1$ the closed-form ratios show it always slows the motor, and for $q=1$ the slowdown is severe.
- The correlation between flux and $\Delta G^\ddagger$ supports a design rule: to maximize speed, place the intermediate state at the point that minimizes the sum of the two substep barriers.
Reading between the lines
- Beyond the paper, the load-distribution argument suggests a lever-like mechanism: an intermediate well near the midpoint cuts the fraction of the load acting on each individual barrier, so adding more substeps could further soften load sensitivity, though the paper does not compute this.
- Beyond the paper, fitting $q$ and $k_i$ to dwell-time data from motors such as $F_1$-ATPase could turn the design rule into a quantitative test across species with different substep sizes.
- Beyond the paper, the model gives synthetic-motor designers a concrete recipe—insert a weakly bound intermediate potential near the midpoint—and a measurable signature: the speedup over the direct motor should grow when an external force opposes the motion.
- Beyond the paper, the small-$\tilde{w}$ slowdown for extreme $q$ warns that the benefit is kinetic, not thermodynamic: when reactions are slow relative to diffusion, splitting $\Delta G$ reduces the driving force without the compensating barrier-lowering, so rate asymmetry must be engineered deliberately.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies a one-dimensional reaction-diffusion model of a molecular motor that can visit an intermediate chemical state between adjacent primary potential wells. The primary and intermediate potentials are harmonic, and the potential-switching rates obey local detailed balance with a single asymmetry parameter q. The authors compare the steady-state flux with and without the intermediate state, optimize the intermediate position (Δx_i, ΔG_i), and also examine the Stokes efficiency. They report that, for moderately large bare switching rates w̃, the optimized intermediate-state flux J* exceeds the direct flux J0 for every q between 0 and 1, and that the acceleration is larger under an external load. In the small-w̃ limit the intermediate state is shown analytically to reduce the flux for the extreme values q=0 and q=1. The central finite-w̃ results are obtained with a numerical steady-state Fokker-Planck solver cited to an unpublished same-group preprint, Ref. [24], and no code or data are supplied.
Significance. If the finite-w̃ results are correct, the paper offers a concrete and potentially useful design principle for placing intermediate states in molecular motors: position them so as to reduce the effective barrier height, which is particularly beneficial under load. The analytic small-w̃ formulas in Eqs. (23)-(30) are clean and correct, and the potential-of-mean-force construction in the large-w̃ limit is standard. The paper also clearly states all model assumptions and fits no parameters to the flux data, which is a strength. However, the headline claims for the finite-w̃ regime, including the universal 'any q' statement and the load-dependent speedup, rest on a numerical solver that is not described, whose implementation is not released, and for which no convergence checks are shown. As a result, the quantitative reach of the conclusions is currently not independently verifiable, which limits the paper's immediate impact despite the plausibility and partial analytic support of its central idea.
major comments (3)
- [II.B and III.A (Eqs. (8)-(11), Figs. 3-8)] The finite-w̃ steady-state flux is computed by 'a method previously reported in [24]', but Ref. [24] is an unpublished same-group preprint and the current manuscript gives no code, data, discretization parameters, boundary conditions, or convergence checks. Since the central claims—acceleration 'in most cases', acceleration 'at any given q', and the load-dependent effect—are all made for finite w̃ (Figs. 3b, 4c, 5, 6, and 8), this numerical core is a black box. I request a full description of the solver, a convergence study, and at least one nontrivial benchmark against the small-w̃ or large-w̃ analytic limits, or release of the code and data.
- [III.A, paragraph after Fig. 3] The statement that for relatively large w̃ 'J* is larger than J0 at any given q between 0 and 1' is a universal claim about q, but the manuscript does not display a finite-w̃ q-sweep. The finite-w̃ results shown in Figs. 4c, 5, and 6 use q=0, and Fig. 3b does not specify the q values or line styles used for the finite-w̃ curves. Please add a finite-w̃ J*/J0 versus q sweep, or otherwise delimit the claim to the particular q values actually simulated.
- [III.C, Eqs. around Fig. 7 and Fig. 8] The external-load model is not written down. The text says the load is introduced as 'the inclination of the potential landscape' and that the work per cycle is W_ext, but the explicit tilted potential (for instance U(x)-Fx) and the definition of the load distribution factor θ in terms of that potential are not given. Without these expressions, neither the analytic large-w̃ curves in Fig. 7 nor the finite-w̃ curves in Fig. 8 can be reproduced. Please state the loaded potential and the precise relation between W_ext and the tilt amplitude.
minor comments (5)
- [Eq. (13)] The integrand in the definition of w^{-,eff} is written as w^-_n(x), but the backward transition considered in Eq. (2) is from state (p,n+1) to (p,n). Please clarify the index convention or the shift used to rewrite the integral over P^{eq}_{p,n}(x).
- [Eq. (21)] The text states that the total turnover time is proportional to e^{ΔG‡p/kBT}+e^{ΔG‡i/kBT}, but Eq. (21) defines ΔG‡ using the average of the two exponentials rather than their sum. The factor of 2 does not affect the optimization, but the definition should be stated as a convention and justified.
- [Fig. 6 caption] There is a typo in the caption: 'sigfnificantly' should be 'significantly'.
- [Figs. 3b and 8] The captions should specify the q values and line styles used for each curve. In particular, Fig. 3b mixes small-w̃ and finite-w̃ results without stating which q is used, and Fig. 8 says 'various q' but does not enumerate them.
- [II, paragraph after Eq. (6)] The statement that 'the conclusions of this paper do not qualitatively change with the choice of these parameters' is broader than what is demonstrated. Please either provide supporting parameter sweeps or soften the claim to the parameter ranges actually tested.
Circularity Check
No circularity: flux is computed from stated Fokker-Planck dynamics; effective barrier is a post-hoc diagnostic, not a fitted input.
full rationale
The paper's central derivation is self-contained. The steady-state flux with and without intermediate states is defined by Eqs. (7)-(11) directly from the stated Fokker-Planck equations, with all parameters (ΔG, Δx, γ, kp, q, w̃, ΔGi, Δxi, ki) given as model inputs; no constants are fitted to any target flux, and J* is obtained by searching over (ΔGi, Δxi) rather than by imposing the acceleration claim. The effective barrier height ΔG‡ in Eq. (21) is defined from the potential intersections after the flux is computed, and is used only as a heuristic explanation; the optimization of flux is not performed by minimizing ΔG‡, so the 'lowering the effective barrier' statement is not definitionally equivalent to the computed flux increase. The small-w̃ analytic results (Eqs. (23)-(30)) are derived by quadrature from the same rates and are independent checks, not inputs. The self-citations (rate model [21], numerical method [24]) are methodological and contextual; [24] is a same-group preprint without released code, which is a reproducibility limitation for the finite-w̃ numerics, but it is not used as a logical premise that forces the central claim. There is no step where a fitted parameter is renamed a prediction, no uniqueness theorem imported from the authors, and no ansatz smuggled via citation. Therefore no significant circularity.
Assumptions & free parameters
free parameters (8)
- chemical free energy drop per step ΔG =
10 k_B T
- stepping distance Δx =
6 nm
- friction coefficient γ =
100 k_B T s/nm^2
- primary spring constant k_p =
2.4 k_B T/nm^2
- bare switching rate w̃ =
10^-2 s^-1 for the finite case; varied in Fig. 3a
- asymmetry parameter q =
varied 0 to 1
- intermediate spring constant k_i =
varied relative to k_p
- intermediate position (Δx_i, ΔG_i) =
optimized over range
assumptions (7)
- standard math Overdamped Langevin dynamics with Gaussian white noise and constant friction coefficient (Eq. 1).
- domain assumption Local detailed balance determines the ratio of forward and backward switching rates (Eqs. 2, 4, 5).
- ad hoc to paper Switching rate x-dependence is prescribed by a single asymmetry parameter q, with the same q for both substeps (Eqs. 3, 6).
- ad hoc to paper Comparison uses bare rate 2w̃ for the two-step intermediate cycle versus w̃ for the direct step.
- domain assumption Primary and intermediate potentials are harmonic wells with a single intermediate state per cycle.
- domain assumption In the large-w̃ limit, the dynamics is described by the potential of mean force (Eq. 20).
- ad hoc to paper Arrhenius form for barrier crossing, used to define the effective barrier height ΔG‡ (Eq. 21).
Cite this review
Pith. "Pith review of Acceleration of enzymatic reaction-diffusion kinetics by intermediate state." pith.science (2026). https://pith.science/paper/74SGOWUG
@misc{pith2026250517130,
author = {Pith},
title = {Pith review of: Acceleration of enzymatic reaction-diffusion kinetics by intermediate state},
year = {2026},
howpublished = {\url{https://pith.science/paper/74SGOWUG}},
note = {Machine review of arXiv:2505.17130}
}
read the original abstract
Biological molecular motors are high-performance nanomachines that convert chemical energy into mechanical motion via chemomechanical coupling. Their reaction cycles typically comprise a series of intermediate chemical states between the initial and final primary states. However, the influence of these intermediate states on motor performance has not yet been fully explored. In this study, we investigate the impact of intermediate states on the motor kinetics using a reaction-diffusion model. In most cases, the intermediate states accelerate the motor by lowering the effective barrier height. This acceleration is particularly pronounced when an external load is applied to the motor, implying the practical importance of the intermediate states. The intermediate states can also slow down the reaction in some cases, such as the slow reaction limit with asymmetric kinetics. Our findings provide practical insights into the design principles behind the high performance of biological molecular motors, as well as the development of efficient artificial molecular motors.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[24]
T. L. Hill,Free Energy Transduction and Biochemical Cycle Kinetics(Springer-Verlag, 1988)
work page 1988
-
[1]
7b] and the increase in ∆G ‡ 0 −∆G ‡ under an external load [Fig
We observed the positive correlation between the acceleration [Fig. 7b] and the increase in ∆G ‡ 0 −∆G ‡ under an external load [Fig. 7c]. In general, the bar- rier height approximately increases at a rate ofW extθby external load. Here,θis called the load distribution fac- tor [2, 29] and is defined here as the distance between the potential minimum and ...
-
[2]
B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, and P. Walter,Molecular Biology of the Cell. 6th ed. (W.W.Norton, 2014)
work page 2014
-
[3]
Howard,Mechanics of Motor Proteins and the Cy- toskeleton(Sinauer Associates Inc, 2001)
J. Howard,Mechanics of Motor Proteins and the Cy- toskeleton(Sinauer Associates Inc, 2001)
work page 2001
-
[4]
K. Svoboda, C. F. Schmidt, B. J. Schnapp, and S. M.470 Block, Direct observation of kinesin stepping by optical trapping interferometry, Nature365, 721 (1993)
work page 1993
-
[5]
T. Funatsu, Y. Harada, M. Tokunaga, K. Saito, and T. Yanagida, Imaging of single fluorescent molecules and individual ATP turnovers by single myosin molecules in aqueous solution, Nature374, 555 (1995)
work page 1995
- [6]
-
[7]
J. P. Abrahams, A. G. W. Leslie, R. Lutter, and J. E.480 Walker, Structure at 2.8 ˚A resolution of F1-ATPase from bovine heart mitochondria, Nature370, 621 (1994)
work page 1994
Show all 35 references
-
[8]
H. Noji, R. Yasuda, M. Yoshida, and K. Kinosita, Direct observation of the rotation of F 1-ATPase, Nature386, 299 (1997)
1997
-
[9]
Yasuda, H
R. Yasuda, H. Noji, K. Kinosita, and M. Yoshida, F 1- ATPase is a highly efficient molecular motor that rotates with discrete 120 ◦ steps, Cell93, 1117 (1998)
1998
-
[10]
Noji and H
H. Noji and H. Ueno, How does F 1-ATPase generate torque?: Analysis from cryo-electron microscopy and ro-490 tational catalysis of thermophilic F 1, Front. Microbiol. 13, 904084 (2022)
2022
-
[11]
Yasuda, H
R. Yasuda, H. Noji, M. Yoshida, K. Kinosita, and H. Itoh, Resolution of distinct rotational substeps by submillisec- ond kinetic analysis of F 1-ATPase, Nature410, 898 (2001)
2001
-
[12]
Zarco-Zavala, R
M. Zarco-Zavala, R. Watanabe, D. G. G. McMillan, T. Suzuki, H. Ueno, F. Mendoza-Hoffmann, J. J. Garc ´ ıa- Trejo, and H. Noji, The 3×120 ◦ rotary mechanism of Paracoccus denitrificans F1-ATPase is different from that500 of the bacterial and mitochondrial F 1-ATPases, Proc. Nat...
2020
-
[13]
R. R. Watanabe, B. T. Kiper, M. Zarco-Zavala, M. Hara, R. Kobayashi, H. Ueno, J. J. Garc ´ ıa-Trejo, C.-B. Li, and H. Noji, Rotary properties of hybrid F1-ATPases consist- ing of subunits from different species, iScience26, 106626 (2023)
2023
-
[14]
Suzuki, K
T. Suzuki, K. Tanaka, C. Wakabayashi, E. Saita, and M. Yoshida, Chemomechanical coupling of human mito- chondrial F 1-ATPase motor, Nat. Chem. Biol.10, 930510 (2014)
2014
-
[15]
A. I. Brown and D. A. Sivak, Allocating dissipation across a molecular machine cycle to maximize flux, Proc. Nat. Acad. Sci.114, 11057 (2017)
2017
-
[16]
J. A. Wagoner and K. A. Dill, Mechanisms for achieving high speed and efficiency in biomolecular machines, Proc. Nat. Acad. Sci.116, 5902 (2019)
2019
-
[17]
J¨ ulicher, A
F. J¨ ulicher, A. Ajdari, and J. Prost, Modeling molecular motors, Rev. Mod. Phys.69, 1269 (1997)
1997
-
[18]
Wang and G
H. Wang and G. Oster, Energy transduction in the F 1520 motor of ATP synthase, Nature396, 279 (1998)
1998
-
[19]
Elston, H
T. Elston, H. Wang, and G. Oster, Energy transduction in ATP synthase, Nature391, 510 (1998)
1998
-
[20]
Keller and C
D. Keller and C. Bustamante, The mechanochemistry of molecular motors, Biophys. J.78, 541 (2000)
2000
-
[21]
Toyabe, H
S. Toyabe, H. Ueno, and E. Muneyuki, Recovery of state- specific potential of molecular motor from single-molecule trajectory, EPL97, 40004 (2012)
2012
-
[22]
Kawaguchi, S.-i
K. Kawaguchi, S.-i. Sasa, and T. Sagawa, Nonequilibrium dissipation-free transport in F 1-ATPase and the thermo-530 dynamic role of asymmetric allosterism, Biophys. J.106, 2450 (2014)
2014
-
[23]
R. A. Marcus, Electron transfer reactions in chemistry. theory and experiment, Rev. Mod. Phys.65, 599 (1993)
1993
-
[25]
Nakayama and S
Y. Nakayama and S. Toyabe, Bypass mechanism of F1-ATPase for asymmetric enzyme kinetics (2024), arXiv:2402.17326 [physics.bio-ph]
2024 arXiv
-
[26]
Michaelis and M
L. Michaelis and M. Menten, Die Gerak der Invertin-540 wirkung, Biochem. Z49, 333 (1913)
1913
-
[27]
Wang and G
H. Wang and G. Oster, The Stokes efficiency for molec- ular motors and its applications, EPL57, 134 (2002)
2002
-
[28]
Li and S
C.-B. Li and S. Toyabe, Efficiencies of molecular motors: a comprehensible overview, Biophys. Rev.12, 419 (2020)
2020
-
[29]
Okuno, R
D. Okuno, R. Iino, and H. Noji, Rotation and structure of FoF1-ATP synthase, J. Biochem.149, 655 (2011)
2011
-
[30]
A. B. Kolomeisky and M. E. Fisher, Molecular motors: a theorist’s perspective, Annu. Rev. Phys. Chem.58, 675 (2007).550
2007
-
[31]
E. M. Geertsema, S. J. van der Molen, M. Martens, and B. L. Feringa, Optimizing rotary processes in synthetic molecular motors, Proc. Nat. Acad. Sci.106, 16919 (2009)
2009
-
[32]
J. A. Wagoner and K. A. Dill, Molecular motors: Power strokes outperform Brownian ratchets, J. Phys. Chem. B 120, 6327 (2016)
2016
-
[33]
Kassem, T
S. Kassem, T. van Leeuwen, A. S. Lubbe, M. R. Wil- son, B. L. Feringa, and D. A. Leigh, Artificial molecular motors, Chem. Soc. Rev.46, 2592 (2017).560
2017
-
[34]
Bath and A
J. Bath and A. J. Turberfield, DNA nanomachines, Nat. Nanotech.2, 275 (2007)
2007
-
[35]
Z. Wang, R. Hou, and I. Y. Loh, Track-walking molecular motors: a new generation beyond bridge-burning designs, Nanoscale11, 9240 (2019)
2019
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.