{"id":"8844970f-df49-41ae-b8a0-ac2d07b4323e","arxiv_id":"2501.09584","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For long semiflexible DNA chains under tension, the looping time scales as the inverse of the equilibrium looping probability, recovering the quadratic-exponential force dependence at low forces.","lead":"This paper proposes a simple formula for how fast a piece of DNA can form a loop when pulled by a force, and tests it against computer simulations. It argues that the looping time is simply the inverse of the probability that the DNA happens to be looped, which works well for long DNA chains.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim rests on an unproven equivalence between mean first-passage time from the unlooped state and inverse equilibrium looping probability; the cited Kac theorem concerns recurrence times, and the simulation validation only matches shifted shapes, not absolute values.","rationale":"The reader's weakest assumption identifies exactly the load-bearing step of the paper. I independently arrive at the same point: the derivation of Eq. 8 is not a theorem for the measured first-passage time, and the validation in Fig. 3 only checks the force-scaling after an arbitrary vertical shift, so it cannot confirm the absolute inverse proportionality τ ∼ 1/P_loop. My two-state Markov-chain example shows the assumption is not generally valid, which raises the correctness risk from 'medium' to a level requiring an explicit test. However, the paper may still be correct for the specific semiflexible chains studied; the simulation appears to agree in shape, and the paper includes a zero-force length-scaling comparison in the Supplementary Material. The right disposition is therefore to keep the reader's CONDITIONAL verdict rather than reject outright: the condition is a direct test of τ(f)P_loop(f) without rescaling. I agree with the reader's recommendation that the derivation of the inverse scaling and broader validation would strengthen the claim.","tokens_in":9283,"tokens_out":13175,"duration_ms":141401,"concrete_test":"Use the same LAMMPS model: from long equilibrium trajectories, measure P_loop(f) directly as the fraction of time with r ≤ rc; in separate runs measure τ(f) as the mean first-passage time to r ≤ rc from an equilibrated unlooped start. Plot τ(f) P_loop(f) versus f for N = 23, 30, 40 with no vertical rescaling. If this product is not flat in f, the missing prefactor in Eq. 8 depends on f and the model's predicted exp(f^2) form is not supported. As a minimal analytical check, solve the two-state Markov chain with rates k_on, k_off matching the observed looping rate and P_loop; verify that 1/k_on ≠ (k_on+k_off)/k_on, showing Kac's theorem does not apply to the first-passage quantity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive step sits in the first paragraph of Section III: the paper asserts τ(f) ∼ 1/P_loop(f) by invoking the Markov-chain recurrence theorem, then combines this with P_loop(f) ∼ 1/Z(f) to obtain Eq. 8, τ(f) ∼ exp(−βF(⟨x(f)⟩)). The cited theorem is Kac's recurrence theorem: the mean recurrence time to a state, starting from that state, is the inverse of its stationary probability. The measured quantity is not a recurrence time; it is the mean first-passage time to the looped set starting from an equilibrated unlooped configuration. For a two-state chain with U→L rate k_on and L→U rate k_off, the first-passage time from U is 1/k_on, whereas 1/P_loop = (k_on+k_off)/k_on; these agree only in the absorbing-loop limit k_off ≪ k_on, which is not the situation in the simulations, where a 'loop' is merely r ≤ rc and can break immediately. The equivalence therefore carries an extra dynamical assumption that is not derived. The numerical support in Fig. 3 also does not test the proportionality: each theory is vertically shifted to one simulation point, so only the f-dependence of τ(f)/τ(0) is compared, not τ(f)P_loop(f). The close relation of Eq. 8 to the two-state model (recovered for xc→0) makes this assumption the sole load-bearing distinction of the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reviews two existing theories for tension-dependent DNA looping — the two-state model of Blumberg et al. and the barrier-escape approach of Shin et al. — and compares them with optical-tweezer data and with new Langevin simulations of semiflexible chains. The authors propose a third model in which the looping time is inversely proportional to the equilibrium looping probability, leading for a wormlike chain to Eq. (8), τ(f) ∼ exp[-βF(⟨x(f)⟩)], and hence to exponential-in-f^2 growth at low force. They find that this model agrees with simulations for long chains (L ≳ 4l_p) over the whole force range tested, while the prior theories fail outside a narrower parameter regime.","tokens_in":9652,"tokens_out":6721,"duration_ms":67788,"significance":"The paper has several strengths: it provides a compact comparison of existing theories, identifies a simple and analytically evaluable prediction, tests the prediction over a wider force range than the original experiment, and includes a finite-size check for the shortest chain. If Eq. (8) is established, the low-force f^2 prediction would resolve a concrete controversy and be directly useful for designing DNA-looping experiments. However, the central inverse-scaling assumption is imported from previous work rather than derived, and the simulation validation is shape-based rather than absolute; these points need to be addressed before the main claim can be considered established.","major_comments":[{"comment":"The central scaling τ(f) ∼ 1/P_loop(f) is introduced by invoking the Markov-chain recurrence theorem, but the theorem gives the mean recurrence time to a state starting from that state, whereas the simulations measure the mean first-passage time from an equilibrated unlooped configuration to the looped set r ≤ r_c. In a two-state description these two quantities are 1/k_on and (k_on + k_off)/k_on respectively; they coincide only in the absorbing-loop limit k_off ≪ k_on, which is not established for the capture-radius definition used here. The equivalence should either be derived for semiflexible chains or explicitly identified as an additional dynamical assumption, with evidence that any missing prefactor is force-independent over the tested range.","section":"Section III, first paragraph and Eq. (8)"},{"comment":"The simulation comparison does not actually test the proportionality τ ∼ 1/P_loop, because each theory curve is vertically shifted to coincide with one simulation point ('we made them coincide with the simulation for the most upper data point at every length'; the supplementary figure states that the lines are 'shifted arbitrarily'). This compares only the shape of τ(f)/τ(0), not the absolute relation τ(f)P_loop(f). A plot of the unscaled product τ(f)P_loop(f) versus f for the long-chain data, or an absolute comparison with a single independently justified normalization across all lengths and forces, is needed to support Eq. (8).","section":"Fig. 3 and Supplementary Fig. 4"},{"comment":"The statement G(f) = F(⟨x(f)⟩) is used to pass from the partition function to the force-extension free energy, but this identity is not generally valid for a wormlike chain; it would follow from a saddle-point evaluation of Z(f), whose accuracy is not examined. Without a check against the simulated ⟨x(f)⟩ or a numerical evaluation of Z(f), Eq. (8) is not a parameter-free consequence of the inverse-scaling assumption, and the claimed reduction to the Blumberg result for x_c → 0 remains a limiting statement rather than a derivation.","section":"Section III, after Eq. (7)"},{"comment":"The approximation P_loop(f) ∼ 1/Z(f) neglects exp(β f x) in the looped region, which requires β f r_c ≪ 1. With the simulation parameters β = 1 and r_c = 1, this condition is violated for f ≳ 0.3, i.e., for a substantial part of the force range in which the model is claimed to work. The authors should quantify the error introduced by this approximation or restrict the low-force claim accordingly.","section":"Section III, Eqs. (6)–(7)"}],"minor_comments":[{"comment":"'Inversionally proportional' should read 'inversely proportional'; the same typo appears in the Conclusion.","section":"Section III, after Eq. (8)"},{"comment":"The citation 'the two-state model of Blumberg et al. [10]' appears to be a numbering error: Ref. [9] is the two-state paper, while Ref. [10] is Shin et al.","section":"Introduction, paragraph 2"},{"comment":"'Most upper data point' should be 'uppermost data point', and the vertical rescaling constants C for each chain length should be stated in the main text, not only in the caption, because they are part of the comparison procedure.","section":"Fig. 3 caption"},{"comment":"'A loop is defined as the mean first passage time for r to become smaller than or equal to the capture radius' should be 'The looping time is defined as the mean first passage time for r to become smaller than or equal to the capture radius'.","section":"Section IV, Simulation"},{"comment":"'Jhon Willey & Sons, Inc.' should be 'John Wiley & Sons, Inc.'.","section":"Supplementary Material, Ref. [14]"},{"comment":"The phrase 'quadratic-exponential increase' would be clearer if it were stated once as 'an increase exponential in f^2', since this is the specific functional form being claimed.","section":"Section V, Results"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a statistical-mechanics/soft-matter journal. The main risk is that the central result is presented as a derivation but in fact rests on an unstated dynamical assumption, and the simulation validation in Fig. 3 is only shape-based. I would welcome a revision that adds an absolute test of τ(f)P_loop(f) and that quantifies the approximations leading to Eq. (8). The self-citation to Ref. 13 is relevant, but the simulation results for semiflexible chains are new and useful."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the Laeremans–Ellenbroek Letter. The useful contribution is a clean demonstration that for long semiflexible chains (L ≳ 4lp), the looping time under tension follows the inverse equilibrium looping probability, giving τ(f) ~ exp(βF(⟨x⟩)) which grows as exp(const × f²) at low force. They validate this across the full force range with Langevin simulations for several chain lengths, and they correctly point out that the two-state and barrier-escape theories fail in this regime for different reasons. The agreement with the Chen et al. experiment is honestly qualified—they note the experiment probes a narrow regime where the barrier-escape theory accidentally matches.\n\nThe main soft spot is the theoretical justification of the inverse scaling. The paper invokes Kac's recurrence theorem, but τ as measured is a first-passage time from an unlooped configuration, not a recurrence time from the looped state. The two-state calculation in the stress-test note shows these differ unless the residence time in the loop is negligible, which is plausible for a small capture radius but not shown. So the scaling is an assumption, supported by the zero-force results in the supplement and by the shape agreement in Fig. 3, but not derived. The paper also treats G(f)=F(⟨x(f)⟩) as a thermodynamic identity; it's a saddle-point approximation, which is fine for long chains but should be stated. And because the theories are vertically shifted, the simulations validate the force dependence of τ(f)/τ(0), not the absolute proportionality constant. For a Letter, this is acceptable if the claim is about scaling.\n\nThe paper is honest: it explicitly says the theories only provide a scaling, it tests finite-size effects, and it spells out when each model works. The reduction to the two-state model for xc→0 is a nice clarification.\n\nI'd send this to peer review. A referee should push for a more careful derivation of the inverse scaling (or at least a statement of when it holds) and for a direct test of τ(f) P_loop(f) rather than just shifted shapes. But the core result is likely correct and useful for people modeling DNA looping. I'd bring it to a reading group as an example of how simulation can discriminate between theories.","headline":"Useful, honest Letter; inverse-probability scaling for looping time under tension is well supported for long chains, but the scaling is assumed rather than derived.","tokens_in":10128,"tokens_out":7398,"would_cite":true,"duration_ms":76201,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82D60","82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"For long semiflexible chains, the tension-dependent DNA looping time is the inverse of the equilibrium looping probability, giving exponential-in-$f^2$ growth at low force and agreement with simulation over the full force range.","keywords":["DNA looping","semiflexible polymer","looping time","equilibrium looping probability","worm-like chain","force-dependent kinetics","mean first passage time","two-state model"],"falsifier":"A direct test would be an optical-tweezer measurement of $\\tau(f)$ on a DNA construct with $L \\gtrsim 4l_p$ at forces below roughly 80 fN: the model predicts that $\\log \\tau$ grows quadratically with force, while the barrier-escape theory predicts almost linear growth in force, so the curvature of the measured curve decides.","tokens_in":9100,"feed_emoji":"🧬","tokens_out":8727,"duration_ms":88027,"temperature":0.7,"pith_summary":"The paper sets out to settle which theoretical description correctly predicts how tension slows DNA looping. It argues that for long semiflexible chains, meaning about four persistence lengths or more, the looping time $\\tau(f)$ is simply the inverse of the equilibrium looping probability, $\\tau(f) \\sim 1/P_{\\rm loop}(f)$, over the whole force range. That relation yields exponential-in-$f^2$ growth at low force, matching the older two-state model's scaling but for a different reason, and it avoids the single-reaction-coordinate free energy on which the rival barrier-escape theory depends. The paper supports the claim with bead-spring simulations and uses it to explain why the barrier-escape theory only seemed to fit existing optical-tweezer data in a narrow force window. If correct, it gives a practical route to predict looping kinetics from equilibrium loop-closure probabilities and makes low-force predictions that future experiments can test.","feed_headline":"For long DNA, looping time tracks inverse loop probability","feed_subtitle":"New model matches simulations over all forces and settles whether looping slows exponentially in force or in force squared.","key_machinery":"The load-bearing object is the inverse scaling relation $\\tau(f) \\sim 1/P_{\\rm loop}(f)$, where $P_{\\rm loop}(f)$ is the equilibrium probability that the chain's end-to-end distance is within a capture radius. It is obtained by writing the equilibrium distribution as $e^{-\\beta F(0,f)}/Z(f) \\sim 1/Z(f)$, using the Markov-chain recurrence-time relation, and then translating $Z(f)$ into $F(\\langle x(f)\\rangle)$ with the thermodynamic identity $G(f)=F(\\langle x(f)\\rangle)$. The force-extension curve of the worm-like chain, evaluated analytically, closes the calculation. This machinery matters because it replaces the barrier-escape reaction coordinate with a property that can be computed from equilibrium loop-closure probability alone.","core_discovery":"The paper's central claim is that, for a semiflexible chain under tension, the mean looping time is inversely proportional to the equilibrium looping probability for long chains, and that this proportionality is the correct organizing principle for the force dependence. Concretely, because a loop corresponds to an end-to-end distance near zero, $P_{\\rm loop}(f) \\sim 1/Z(f)$, and then $\\tau(f) \\sim e^{-\\beta F(\\langle x(f)\\rangle)}$. This reproduces the quadratic-exponential growth of the two-state model in the low-force limit but interprets it as an equilibrium-probability statement rather than a two-state energy gap. The paper validates the claim with simulations for $L/l_p$ from 2.2 to 8, finding that the inverse-probability model tracks the whole force range for $L/l_p \\gtrsim 4.6$, while the two prior theories fail except in narrow regimes.","pith_inferences":["Editorial extension: any method that computes equilibrium loop-closure probabilities, such as J-factor calculations for sequence-dependent or supercoiled DNA, becomes a kinetic prediction for long chains; this link is implicit in the paper but not developed there.","Editorial extension: the paper's explanation of the short-chain breakdown, that the capture radius is not small relative to chain length, suggests that a finite-capture-radius correction to $1/P_{\\rm loop}$ could extend the theory to shorter DNA; the authors do not propose one.","Editorial extension: because the reduction to the two-state model is described as a mathematical coincidence rather than a physical mechanism, fitting experimental looping-time data to a two-state Arrhenius form could be misleading; the paper hints at this risk when discussing protein size, and the warning could be sharpened by fitting simulated data to both functional forms."],"forward_implications":["For chains with $L \\gtrsim 4l_p$, looping kinetics can be predicted directly from equilibrium looping probabilities, avoiding the fragile reaction-coordinate free energy of the barrier-escape route.","The low-force scaling debate is settled in favor of exponential-in-$f^2$ growth, matching the two-state model's functional form but for a different reason.","The apparent experimental agreement with the barrier-escape model in the 60--180 fN range is not strong evidence for that model, since the simulations show deviations just outside that range.","The model gives force-dependent predictions that are within reach of optical-tweezer experiments at sub-piconewton forces.","Because the derivation is not tied to a specific polymer model, the same inverse-probability relation should hold for other semiflexible chains once their equilibrium looping probability is known."],"supporting_citations":[{"why":"Supplies the inverse scaling between looping time and looping probability for a freely jointed chain under tension, and the evidence that the barrier-escape approach can fail because local equilibrium is not satisfied.","marker":"[13]"},{"why":"Establishes the zero-force result that semiflexible-chain looping time is inversely proportional to the equilibrium looping probability, the idea extended here to finite tension.","marker":"[15]"},{"why":"Provides the two-state model whose exponential-in-$f^2$ prediction is recovered in the lowest-force limit and whose meaning is reinterpreted.","marker":"[9]"},{"why":"Provides the barrier-escape model, the main rival prediction, and the comparison baseline whose force-range validity is shown to be limited.","marker":"[10]"},{"why":"Provides the optical-tweezer experimental data against which the prior theories and the new model are compared.","marker":"[8]"},{"why":"Supplies the force-extension relation used in constructing the free energy for both the two-state model and the inverse-probability model.","marker":"[12]"},{"why":"Supplies the analytical force-extension approximation that lets the new model be evaluated fully analytically.","marker":"[30]"},{"why":"Supplies the numerical scheme for computing equilibrium loop-closure probability used to validate the inverse scaling at zero force.","marker":"[27]"}],"fun_headline_variants":["New model matches DNA looping time across all forces","Long DNA looping time: inverse loop probability wins","DNA tension: looping time equals inverse loop chance","Semiflexible DNA looping time follows probability","Predicting DNA looping time from loop probability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the assumption, imported from earlier work rather than proved here, that the average time for an unlooped chain to first form a loop equals the inverse of the equilibrium looping probability; if that equality fails, the model's predictions collapse.","fun_headline_variants_meta":{"raw":{"variants":["New model matches DNA looping time across all forces","Long DNA looping time: inverse loop probability wins","DNA tension: looping time equals inverse loop chance","Semiflexible DNA looping time follows probability","Predicting DNA looping time from loop probability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000176,"raw_usage":{"total_tokens":1220,"prompt_tokens":806,"completion_tokens":414,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":422,"completion_tokens_details":{"reasoning_tokens":344}},"tokens_in":422,"tokens_out":414,"duration_ms":4411,"temperature":1.0,"reasoning_tokens":344,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:52:21.710657+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be an optical-tweezer measurement of $\\tau(f)$ on a DNA construct with $L \\gtrsim 4l_p$ at forces below roughly 80 fN: the model predicts that $\\log \\tau$ grows quadratically with force, while the barrier-escape theory predicts almost linear growth in force, so the curvature of the measured curve decides.","supporting_citations":[{"cited_title":"Polymer dynamics under tension: mean first passage time for looping","cited_arxiv_id":"2410.01347","evidence_quote":"Supplies the inverse scaling between looping time and looping probability for a freely jointed chain under tension, and the evidence that the barrier-escape approach can fail because local equilibrium is not satisfied."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the zero-force result that semiflexible-chain looping time is inversely proportional to the equilibrium looping probability, the idea extended here to finite tension."},{"cited_title":"Blumberg, A","cited_arxiv_id":null,"evidence_quote":"Provides the two-state model whose exponential-in-$f^2$ prediction is recovered in the lowest-force limit and whose meaning is reinterpreted."},{"cited_title":"Shin and W","cited_arxiv_id":null,"evidence_quote":"Provides the barrier-escape model, the main rival prediction, and the comparison baseline whose force-range validity is shown to be limited."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the optical-tweezer experimental data against which the prior theories and the new model are compared."},{"cited_title":"Petrosyan, Improved approximations for some polymer ex- tension models, Rheologica Acta 56, 21 (2017)","cited_arxiv_id":null,"evidence_quote":"Supplies the analytical force-extension approximation that lets the new model be evaluated fully analytically."},{"cited_title":"Sinha and S","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical scheme for computing equilibrium loop-closure probability used to validate the inverse scaling at zero force."}],"review_version":1}