{"id":"17334b3b-e265-41c8-b7cb-352ea47eb9a5","arxiv_id":"2505.05053","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A review of the author's mesoscopic path-integral model for nucleic acids, with applications to twist-stretch behavior, nanochannel stretching, and cyclization of short DNA.","lead":"This paper reviews a coarse-grained model of DNA and RNA that treats each base pair explicitly and uses path integrals to average over thermal fluctuations. It applies the model to show why DNA overwinds and RNA unwinds when stretched, and why short DNA loops more easily than elastic rod models predict.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Twist-stretch explanation is not isolated: no ablation separates tilt, slide, and diameter, so the claim that inclination drives RNA untwisting is untested.","rationale":"I agree with the reader's weakest-assumption identification: the load-bearing premise is that average geometric inputs (gamma, S, and R0) are a sufficient effective representation of the A-form versus B-form difference in a point-nucleotide Hamiltonian. The main advertised result, the opposite twist-stretch sign, is asserted to follow from these structural features, yet the paper never performs the control experiment of varying one geometric input at a time. The Section 4.A limitation paragraph explicitly concedes that the real origin is sugar pucker at the atomic scale, so the explanatory claim stands or falls on whether tilt, slide, and diameter are doing the causal work claimed. A simple ablation of the three A-form inputs would settle whether the sign is genuinely due to gamma or is an artifact of the combined input set. The cyclization calibration issue noted by the reader is real, but it is openly disclosed as a parameter-setting procedure and does not bear on the paper's first and most distinctive claim about twist-stretch. The preprint is honest and internally coherent, so I would keep the reader's CONDITIONAL verdict rather than escalate to rejection; the missing ablation is a concrete, testable gap rather than a demonstrated contradiction.","tokens_in":20124,"tokens_out":10408,"duration_ms":116333,"concrete_test":"Run a four-point ablation for the A-form parameters in the Section 4.A free-energy minimization, keeping all Hamiltonian parameters fixed: (i) R0 = 24 Å, gamma = 15 degrees, S = 0; (ii) R0 = 24 Å, gamma = 0 degrees, S = 0; (iii) R0 = 20 Å, gamma = 15 degrees, S = 0; (iv) R0 = 24 Å, gamma = 15 degrees, S set to a negative value as in Fig. 3(b). If untwisting (d<h>*/dF > 0) occurs in case (ii), the attribution to tilt is false and R0 is the effective control; if it disappears in case (iii), the claim that inclination is primary is supported; if it appears only with both gamma and S, the compound geometric input is what matters. Report the zero-force equilibrium <h>* alongside each case, since the absolute helical repeat is also part of the model's quantitative credibility.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline mechanistic claim is that the opposite twist-stretch response of B-DNA and A-RNA 'follows from the different structural features of the two helices' (Abstract; Section 4.A). In the model these features enter only as fixed inputs: B-form has R0 = 20 Å, gamma = 0, S = 0, while A-form has R0 = 24 Å, gamma = 15 degrees, and S < 0 in Fig. 3(b). The free-energy minimization then yields <h>*(F). The problem is attribution: no calculation is shown that varies gamma, S, and R0 independently. The text asserts that 'the base pair inclination is the primary cause' because all A-form curves with different |S|/d untwist, but a curve with gamma = 0 and A-form R0 is never shown. It is therefore possible that the sign change is produced by the larger radius R0 = 24 Å or by the altered force coupling d_S(gamma, S) in the term -F_ex d_S cos(phi_n), rather than by tilt per se. The limitation paragraph in Section 4.A concedes that a point-nucleotide model 'cannot capture the ultimate origin' of the effect, which is sugar pucker; this makes the sufficiency of (gamma, S, R0) as effective controls the load-bearing premise. If that premise fails, the model is doing no more than encoding the answer in its inputs.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a three-dimensional mesoscopic Hamiltonian model of double-stranded nucleic acids, treated by a finite-temperature path-integral method, and applies it to three problems: the twist-stretch response of B-DNA versus A-RNA, DNA stretching in a cylindrical nanochannel, and the cyclization J-factor of short DNA fragments. The paper claims that the model reproduces the experimentally observed opposite behavior (DNA overtwists, RNA untwists under tension) and attributes this to the different structural features of the two helical forms, in particular base-pair inclination. It also reports that the model yields cyclization J-factors at approximately 100 base pairs that are much larger than worm-like-chain predictions and consistent with FRET measurements. The manuscript is written as a review of the author's prior work, with the present version adding the A-RNA/B-DNA comparison and the cyclization discussion.","tokens_in":20478,"tokens_out":2210,"duration_ms":24863,"significance":"If the claims were fully supported, the model would be a useful coarse-grained alternative to worm-like-chain descriptions for short nucleic acid fragments, and the path-integral cutoff construction is a potentially valuable way to regularize base-pair fluctuation integrals. The cyclization calculation, however, is not an independent test because the parameters rho_n and alpha_n are calibrated to the experimental J-factor that the model then reproduces. The twist-stretch attribution is also not tested by an ablation of the structural parameters. The paper's strengths are the internally consistent computational framework and the explicit treatment of finite-size effects, but the central mechanistic claims need additional support before the results can be accepted as demonstrated.","major_comments":[{"comment":"The cyclization prediction is not independent. The text states that rho_n and alpha_n are estimated by comparing the model J-factor with the experimentally available looping probability, with the J-factor order of magnitude at N approximately 100 used as a reference point and the choice 'not unique.' Consequently, the agreement in Fig. 6(b) is built into the parameter selection, and the claim that short fragments 'maintain a sizeable bendability in line with the experimental data' is not a test of the model. Please either present a parameter-free prediction, show a systematic sensitivity analysis over a physically justified range of (rho_n, alpha_n), or demonstrate that the order-of-magnitude enhancement over WLC survives without any calibration to the FRET data.","section":"Section 4.C, Eq. (4.2)"},{"comment":"The assertion that base-pair inclination is the primary cause of RNA untwisting is not supported by the calculations shown. All A-form curves in Fig. 3(b) combine gamma = 15 degrees with R0 = 24 Angstrom, and no curve with gamma = 0 at the A-form radius is displayed; likewise, the effect of slide is only shown in the presence of tilt. The observed sign difference could therefore arise from the larger radius or from the modified force coupling d_S in the term -F_ex d_S cos(phi_n), rather than from tilt per se. Please provide an ablation varying gamma, S, and R0 independently, or a formal argument isolating the contribution of each parameter.","section":"Section 4.A, Fig. 3(b)"},{"comment":"The paper's headline mechanism is weakened by its own stated limitation. The paragraph near the end of Section 4.A concedes that the ultimate origin of the DNA/RNA difference is the sugar-pucker conformation at atomic scale, which the point-nucleotide Hamiltonian cannot represent. If the model's (gamma, S, R0) inputs are merely effective proxies, then the conclusion that the opposite twist-stretch pattern 'follows from the different structural features' needs to be tempered: the model shows that certain geometric inputs produce the observed sign, but it does not establish that those inputs are the cause rather than a fitted representation of it. Please clarify what, precisely, is being claimed about causality in the twist-stretch relation.","section":"Section 4.A, limitation paragraph"}],"minor_comments":[{"comment":"In the definition of H_b, the kinetic term is written with the index i (mu/2 dot-r_i^2) while the sum runs over n; this should be dot-r_n^2 for consistency.","section":"Section 2, Eq. (2.1)"},{"comment":"The bracket structure in Eq. (3.2) is difficult to parse, especially the product over n ≠ j and the single-trajectory integral. A clearer display separating the closed and open measures would help the reader verify the normalization and the role of the Heaviside constraint.","section":"Section 3.A, Eq. (3.2)"},{"comment":"The figure caption and text state that the calculation is performed at helical repeat h = 10, chosen as an input so N/h is an integer. This choice suppresses the twist-dependent oscillatory J-factor behavior that is standard in cyclization literature; the paper should note more prominently that the comparison to WLC in Fig. 6(b) is made under this special condition.","section":"Section 4.C, Fig. 6"},{"comment":"The abstract describes the work as a review, while Sections 4.A and 4.C present new calculations. Please align the wording (for instance, 'we review and extend') so that the contribution is stated accurately.","section":"Abstract and conclusions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is largely a synthesis of the author's previously published work, and the novelty of the present version is concentrated in the A-RNA/B-DNA twist-stretch comparison and the cyclization discussion. The cyclization result is calibrated rather than predictive, and the twist-stretch attribution needs an ablation study. These are fixable if the authors are willing to reframe the claims and add the missing tests; the underlying computational machinery is credible and worth publishing in revised form. I would not recommend rejection, but acceptance in the current form would overstate the evidentiary status of the central conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a self-review of Zoli's own mesoscopic Hamiltonian model for nucleic acids. The figures and equations are reprinted from his earlier papers; there is no new computation or independent validation. Take the headline claims with that frame.\n\nWhat the paper does well: it gives a clear, compact account of the path-integral method and the fluctuation cutoff, and it is unusually honest about what the model cannot do. The limitation in Section 4.A that a point-nucleotide model cannot capture the ultimate sugar-pucker origin is stated plainly. The channel-stretching part is illustrative but self-consistent.\n\nThe soft spots are real. The cyclization 'prediction' is not independent: Section 4.C says rho_n and alpha_n are set by matching the model J-factor to the experimental looping probability at N~100, and the same data are then used in Fig. 6(b) to claim agreement. That is calibration, not validation. The twist-stretch attribution also has a load-bearing gap. The abstract says the opposite pattern follows from the different structural features, and the text calls base-pair inclination the primary cause. But the A-form inputs differ in three ways at once—R0=24 Å, gamma=15°, and S<0—and no calculation separates them. A curve with A-form radius and gamma=0 is never shown. So the sign change could be driven by the larger diameter or by the force-coupling term, not by tilt. The stress-test note is right about this. The limitation paragraph concedes the model can't represent the atomic source, which makes the 'explains it' language in the abstract stronger than what is actually demonstrated.\n\nWho gets value: someone wanting a single-source summary of Zoli's model and its earlier results. But as a claim of explaining DNA/RNA mechanical differences, it is not there yet. A referee should ask for an ablation (gamma vs S vs R0) and for out-of-sample cyclization tests, ideally with code and data.\n\nMy call: it deserves a serious referee—it's a legitimate review with formal machinery, not a crackpot or a sloppy paper—but I would not cite it as independent evidence for the flexibility claims.","headline":"A self-review of the author's own mesoscopic model—no new results, the cyclization comparison is calibrated to the same data, and the twist-stretch attribution lacks an ablation.","tokens_in":20994,"tokens_out":3000,"would_cite":false,"duration_ms":27357,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A three-dimensional base-pair model reproduces the opposite twist-stretch response of B-DNA and A-RNA and predicts short-DNA cyclization J-factors far above worm-like-chain values.","keywords":["nucleic acids properties","mesoscopic Hamiltonian models","path integral methods","base pair fluctuations","twist-stretch coupling","A-RNA","B-DNA","cyclization J-factor"],"falsifier":"Run the path-integral calculation for the A-form RNA parameter set with $\\gamma = S = 0$ while keeping all other parameters fixed: the model must then produce B-DNA-like overtwisting. Any experiment or atomistic simulation showing that RNA-like untwisting persists in this geometric setting would refute the central claim that the opposite twist-stretch pattern follows from the helical-form tilt and slide.","tokens_in":19869,"feed_emoji":"🧬","tokens_out":7198,"duration_ms":68931,"temperature":0.7,"pith_summary":"This review article argues that a three-dimensional mesoscopic model, in which each nucleotide is a point and base pairs breathe through a finite ensemble of radial fluctuations, can reproduce two experimentally observed properties of short nucleic acids that elastic-rod models miss. First, under tension the model makes B-DNA overtwist and A-RNA untwist, tracing the opposite behavior to the helical-form geometry, base-pair tilt and slide in the A-form, rather than to sequence. Second, for roughly 100-base-pair DNA the model gives cyclization J-factors orders of magnitude larger than worm-like-chain predictions, in line with single-molecule FRET measurements. The author presents this as evidence that modeling flexibility at the level of individual base pairs is needed at short length scales.","feed_headline":"Why DNA overtwists and RNA untwists when pulled","feed_subtitle":"The same 3D base-pair model gives short DNA a looping probability far above worm-like-chain predictions.","key_machinery":"The load-bearing object is the three-dimensional mesoscopic Hamiltonian of Eq. (2.1): a Morse potential $V_1$ for each base pair's hydrogen bonds, a solvent barrier term, and a nonlinear stacking potential $V_2$ between adjacent base pairs that depends on the bending angle $\\varphi_n$ and twist angle $\\theta_n$ through the dimer distance $d_{n,n-1}$. The statistical method is the finite-temperature path integral: each radial coordinate is expanded in Fourier modes, and the integration cutoff on the Fourier coefficients is fixed by a first-passage probability benchmark, $P_j(R_0, 0) \\simeq 1/2$, which also prevents the partition-function divergence that plagues untwisted one-dimensional models. For twist-stretch curves, the paper adds a force term $-F_{\\rm ex}\\, d_S \\cos\\varphi_n$, computes the free energy for a set of helical-repeat values $h$, and selects the equilibrium $\\langle h \\rangle^*$ by free-energy minimization. The A-form geometry enters as average tilt $\\gamma = 15^\\circ$ and slide $S$, which shorten the rise distance $d_S$ and, the calculation shows, drive untwisting under load.","core_discovery":"The central claim is that a discrete Hamiltonian with one radial coordinate per base pair, plus bending and twist angles between adjacent pairs, contains enough physics to explain the twist-stretch dichotomy and short-scale bendability of double-stranded nucleic acids. For the B-form (zero tilt and slide), the equilibrium helical repeat decreases as the applied force increases, meaning DNA overtwists; for the A-form, the same calculation with a base-pair tilt $\\gamma = 15^\\circ$ and a negative slide $S$ produces the opposite trend, meaning RNA untwists under tension. The same partition-function machinery, with a loop-closure constraint, yields J-factors for chains of 60 to 140 base pairs that stay appreciable down to about $10^{-12}$ mol/litre at $N = 80$, whereas the worm-like-chain model drives the looping probability toward zero at short lengths. The paper concludes that the mesoscopic Hamiltonian captures the experimentally observed bendability of about 100-base-pair fragments because it allows large bending angles and independent base-pair fluctuations at every site.","pith_inferences":["One could test the geometric explanation directly by running the same path-integral calculation with $\\gamma = S = 0$ while keeping all other RNA parameters fixed; the model would then predict DNA-like overtwisting, whereas experiments on genuine A-form duplexes with normal ribose pucker would still show untwisting if the atomic-scale sugar pucker is the true cause.","If tilt and slide are indeed sufficient effective coordinates, the same machinery could predict twist-stretch coupling for other helical forms, such as Z-DNA or DNA-RNA hybrids, simply by changing those average geometric inputs.","The first-passage cutoff method is a self-contained way to set fluctuation amplitudes, so it could in principle be transferred to other coarse-grained helical polymer models that suffer from unbounded phase spaces."],"forward_implications":["If the model is right, the opposite twist-stretch responses of double-stranded DNA and RNA are a structural consequence of helical form, not of base sequence or of fine atomic interactions.","The J-factor results imply that fragments around 80 to 100 base pairs are considerably more bendable than worm-like-chain estimates, supporting the idea that short DNA can wrap and loop with ease.","Open-end and fraying effects raise the equilibrium helical repeat and enhance flexibility under moderate loads, so terminal base pairs should be included when predicting short-fragment mechanics.","For channel confinement, the model predicts that reducing the channel diameter straightens the molecule, increasing the end-to-end distance by up to a factor of three at strong confinement."],"supporting_citations":[{"why":"Establishes the twist-stretch relations in nucleic acids from which the A-form geometric inputs, tilt and slide, and the B-form benchmark are taken.","marker":"[47]"},{"why":"Magnetic-tweezers experiment showing that B-DNA overwinds when stretched; the model must reproduce this behavior.","marker":"[92]"},{"why":"Magnetic-tweezers experiment showing that double-stranded RNA untwists under force; sets the opposite target for the model.","marker":"[93]"},{"why":"Comparative molecular dynamics study used to justify the A-form tilt angle of 15 degrees.","marker":"[95]"},{"why":"Single-molecule FRET cyclization data giving a J-factor near $10^{-9}$ mol/litre for about 100 base-pair DNA; the model is calibrated against it.","marker":"[28]"},{"why":"Independent FRET measurement of short-DNA looping probability used as the second experimental reference.","marker":"[29]"},{"why":"Worm-like-chain ring-closure calculation that supplies the contrasting stiff predictions at short lengths.","marker":"[105]"},{"why":"Derives the first-passage probability criterion that sets the base-pair fluctuation cutoff in the path-integral method.","marker":"[81]"},{"why":"Earlier twisted-DNA thermodynamics with solvent interaction from which the J-factor computation and comparison figure are drawn.","marker":"[41]"}],"fun_headline_variants":["Discrete base-pair model reveals DNA overtwist, RNA untwist","Short DNA bendability beyond worm-like chain predictions","Statistical path integral for RNA and DNA mechanics","One model, two opposite twist-stretch patterns","Looping odds for 100-bp DNA from base-pair fluctuations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fragile premise is that the A-form versus B-form difference in twist-stretch response can be represented by two average geometric inputs, base-pair tilt $\\gamma = 15^\\circ$ and slide $S$, in a point-nucleotide Hamiltonian, even though the paper concedes that the true cause is the atomic-scale sugar-pucker conformation, which the model cannot represent.","fun_headline_variants_meta":{"raw":{"variants":["Discrete base-pair model reveals DNA overtwist, RNA untwist","Short DNA bendability beyond worm-like chain predictions","Statistical path integral for RNA and DNA mechanics","One model, two opposite twist-stretch patterns","Looping odds for 100-bp DNA from base-pair fluctuations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000212,"raw_usage":{"total_tokens":1455,"prompt_tokens":1018,"completion_tokens":437,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":355}},"tokens_in":634,"tokens_out":437,"duration_ms":5047,"temperature":1.0,"reasoning_tokens":355,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:15:01.093157+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the path-integral calculation for the A-form RNA parameter set with $\\gamma = S = 0$ while keeping all other parameters fixed: the model must then produce B-DNA-like overtwisting. Any experiment or atomistic simulation showing that RNA-like untwisting persists in this geometric setting would refute the central claim that the opposite twist-stretch pattern follows from the helical-form tilt and slide.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the twist-stretch relations in nucleic acids from which the A-form geometric inputs, tilt and slide, and the B-form benchmark are taken."},{"cited_title":"Duguet, The helical repeat of DNA at high temperature","cited_arxiv_id":null,"evidence_quote":"Magnetic-tweezers experiment showing that B-DNA overwinds when stretched; the model must reproduce this behavior."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Magnetic-tweezers experiment showing that double-stranded RNA untwists under force; sets the opposite target for the model."},{"cited_title":"Slocombe, M","cited_arxiv_id":null,"evidence_quote":"Comparative molecular dynamics study used to justify the A-form tilt angle of 15 degrees."},{"cited_title":"Dickerson, The DNA Helix and How It Is Read","cited_arxiv_id":null,"evidence_quote":"Worm-like-chain ring-closure calculation that supplies the contrasting stiff predictions at short lengths."},{"cited_title":"Zgarbov´ a, M","cited_arxiv_id":null,"evidence_quote":"Derives the first-passage probability criterion that sets the base-pair fluctuation cutoff in the path-integral method."},{"cited_title":"Zoli, J- factors of short DNA molecules J","cited_arxiv_id":null,"evidence_quote":"Earlier twisted-DNA thermodynamics with solvent interaction from which the J-factor computation and comparison figure are drawn."}],"review_version":1}