pith. sign in

arxiv: 2208.11805 · v2 · submitted 2022-08-25 · 🧮 math.DS · q-bio.QM

On the Diffusion Time Evolution of Folding Chains in the Heteropolymer Model

Okezue Bell This is my paper

Pith reviewed 2026-05-24 11:26 UTC · model grok-4.3

classification 🧮 math.DS q-bio.QM
keywords heteropolymer modelprotein foldingdiffusionpower lawLennard-Jones potentialrandomnesstime evolution
0
0 comments X

The pith

Folding amino acid chains in the heteropolymer model follow a power-law diffusion D ∼ t^ν, with the exponent dropping from two-thirds to one-half as Lennard-Jones coupling randomness increases.

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

The paper shows that Iori et al.'s heteropolymer model produces a scaling law for the time evolution of folding chains. Diffusion distance grows as a power of time, and the exponent itself depends on how random the interaction strengths are. A reader would care because this turns a complex folding process into a simple, tunable scaling relation that can be checked in simulations. The result links the degree of disorder in the potential directly to the speed of structural relaxation.

Core claim

In the heteropolymer model, the folding amino acid chain evolves according to a power law D ∼ t^ν. The power ν decreases from 0.66̅ to 0.5 when the randomness of the coupling constants in the Lennard-Jones potential increases.

What carries the argument

The heteropolymer model with a Lennard-Jones potential whose coupling constants have tunable randomness, used to track the diffusion distance D of the folding chain over time.

If this is right

  • The power-law form persists across the full range of coupling randomness examined.
  • Higher randomness systematically lowers the effective diffusion exponent.
  • The model supplies a direct mathematical relation between potential disorder and folding time scale.

Where Pith is reading between the lines

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

  • If the scaling survives in more detailed force fields, folding times for real sequences could be estimated from a single randomness parameter.
  • The 0.5–0.66̅ window might be compared with anomalous diffusion exponents measured in other disordered polymer systems.
  • Varying the temperature or chain length in the same model would test whether the exponent remains pinned to the randomness parameter.

Load-bearing premise

The heteropolymer model from Iori et al. captures the essential dynamics of real protein folding and the reported power law is not produced by the choice of simulation parameters or numerical method.

What would settle it

A simulation run with the same model but different integration step size, potential cutoff, or random seed sequence that yields an exponent outside the 0.5–0.66̅ interval or no power law at all.

Figures

Figures reproduced from arXiv: 2208.11805 by Okezue Bell.

Figure 1
Figure 1. Figure 1: Two comparative D4(t) calculation plots. Here, tMC is the number of full Monte Carlo sweeps of the chain, the scattered dots are the purely harmonic chain, Langevin dynamics normal mode expansion of Equation 7 is the continuous curve. We encounter an exclusion principle due to the discontinuous dynamics of the Metropolis algorithm: when it is known that the correct asymptotic equilibrium distribution is ob… view at source ↗
Figure 2
Figure 2. Figure 2: Distribution plots highlighting key features in quenched noise realizations in varying [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
read the original abstract

In this paper, we mathematically describe the time evolution of protein folding features via Iori et al.'s heteropolymer model. More specifically, we identify that the folding amino acid chain evolve according to a power law $D \sim t^{\nu}$. The power $\nu$ decreases from $0.\overline{66}$ to $0.5$ when the randomness of the coupling constants in the Lennard-Jones potential increases.

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

2 major / 2 minor

Summary. The manuscript applies Iori et al.'s heteropolymer model to the time evolution of folding amino-acid chains. It reports that the diffusion quantity D obeys the power-law scaling D ∼ t^ν, with the exponent ν decreasing from 2/3 to 1/2 as the randomness of the Lennard-Jones coupling constants is increased.

Significance. If the scaling regime is robust and the exponent extraction is insensitive to numerical parameters, the result would supply a concrete, disorder-dependent characterization of anomalous diffusion within a standard heteropolymer model. This could be of interest to mathematical biology and dynamical systems provided the derivation or simulation protocol is fully documented.

major comments (2)
  1. The manuscript provides no description of the fitting procedure used to extract ν, the time interval over which the power law is observed, or any error estimates on the reported values 2/3 and 1/2. Because the central claim is the specific dependence of ν on randomness, the absence of these details makes it impossible to assess whether the quoted exponents reflect a genuine asymptotic regime.
  2. No controls are reported for the sensitivity of the observed scaling to integration timestep, total trajectory length, chain length N, or the number of disorder realizations. The claim that ν drops from 2/3 to 1/2 therefore rests on unverified numerical stability.
minor comments (2)
  1. The abstract uses the non-standard notation 0.66̅; a conventional decimal or fraction (2/3) would improve readability.
  2. The symbol D is introduced without an explicit definition (e.g., mean-squared displacement of which coordinate).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed report and the opportunity to clarify our numerical procedures. We address each major comment below and will revise the manuscript to incorporate the requested documentation and controls.

read point-by-point responses

  1. Referee: The manuscript provides no description of the fitting procedure used to extract ν, the time interval over which the power law is observed, or any error estimates on the reported values 2/3 and 1/2. Because the central claim is the specific dependence of ν on randomness, the absence of these details makes it impossible to assess whether the quoted exponents reflect a genuine asymptotic regime.

    Authors: We agree that the original manuscript lacked a sufficient description of the fitting procedure, time windows, and error estimates. In the revised version we will add an explicit subsection in the Methods that details the linear regression performed on log-log plots of D(t), the precise time intervals over which the power-law regime is fitted for each disorder strength, and the error estimates derived from the variance across independent realizations. This will allow direct evaluation of whether the quoted values of ν correspond to an asymptotic regime. revision: yes

  2. Referee: No controls are reported for the sensitivity of the observed scaling to integration timestep, total trajectory length, chain length N, or the number of disorder realizations. The claim that ν drops from 2/3 to 1/2 therefore rests on unverified numerical stability.

    Authors: The original submission did not include explicit sensitivity tests. We will add a new appendix or subsection presenting additional simulation results that vary the integration timestep, extend total trajectory lengths, consider multiple chain lengths N, and increase the number of disorder realizations. These controls will be used to verify that the reported decrease in ν remains stable within the stated precision. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation is observational from model.

full rationale

The provided abstract and context present the power-law claim as an identification from the Iori et al. heteropolymer model simulations rather than a derivation that reduces to fitted inputs or self-citations by construction. No equations, self-citations, or ansatzes are quoted that would make the exponent ν equivalent to its inputs. The central claim remains an empirical observation from dynamics, with no load-bearing step shown to collapse into a fit or prior author result. This is the normal case of a self-contained numerical study.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the heteropolymer model from prior work and the numerical identification of the power-law scaling; no new entities are introduced.

free parameters (1)
  • randomness level of coupling constants
    The degree of randomness is varied as a control parameter whose effect on ν is reported.
axioms (1)
  • domain assumption The heteropolymer model equations from Iori et al. govern the system dynamics.
    The paper builds directly on this model without re-deriving it.

pith-pipeline@v0.9.0 · 5585 in / 1231 out tokens · 27030 ms · 2026-05-24T11:26:04.085418+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

  1. [1]

    Banerjee

    Ibraheem Alshareedah, Taranpreet Kaur, Jason Ngo, Hannah Seppala, Liz-Audrey Djomnang Kounatse, Wei Wang, Mahdi Muhammad Moosa, and Priya R. Banerjee. Interplay be- tween short-range attraction and long-range repulsion controls reentrant liquid condensation of ribonucleoprotein–rna complexes. Journal of the American Chemical Society , 141(37):14593–14602,...

    work page doi:10.1021/jacs.9b03689 2019

  2. [2]

    Chakraborty

    Arup K. Chakraborty. Disordered heteropolymers: models for biomimetic polymers and polymers with frustrating quenched disorder. Physics Reports , 342(1):1–61, 2001. ISSN 0370-1573. doi: https://doi.org/10.1016/S0370-1573(00)00006-5. URL https://www.sciencedirect.com/ science/article/pii/S0370157300000065

    work page doi:10.1016/s0370-1573(00)00006-5 2001

  3. [3]

    Diffusion of a heteropolymer in a Multi-Interface medium

    Frank den Hollander and Mario V W¨ uthrich. Diffusion of a heteropolymer in a Multi-Interface medium. Journal of Statistical Physics , 114(3):849–889, February 2004

    work page 2004

  4. [4]

    Accurate protein structure prediction by embeddings and deep learning representations

    Iddo Drori, Darshan Thaker, Arjun Srivatsa, Daniel Jeong, Yueqi Wang, Linyong Nan, Fan Wu, Dimitri Leggas, Jinhao Lei, Weiyi Lu, Weilong Fu, Yuan Gao, Sashank Karri, Anand Kannan, Antonio Moretti, Mohammed AlQuraishi, Chen Keasar, and Itsik Pe’er. Accurate protein structure prediction by embeddings and deep learning representations. 2019. doi:10.48550/ARX...

    work page doi:10.48550/arxiv.1911 2019

  5. [5]

    Wick’s theorem at finite temperature

    T.S Evans and D.A Steer. Wick’s theorem at finite temperature. Nuclear Physics B , 474(2): 481–496, 1996. ISSN 0550-3213. doi:https://doi.org/10.1016/0550-3213(96)00286-6. URL https://www.sciencedirect.com/science/article/pii/0550321396002866

    work page doi:10.1016/0550-3213(96)00286-6 1996

  6. [6]

    Kinematics and thermodynamics of a folding het- eropolymer

    M Fukugita, D Lancaster, and M G Mitchard. Kinematics and thermodynamics of a folding het- eropolymer. Proc. Natl. Acad. Sci. U. S. A. , 90(13):6365–6368, July 1993

    work page 1993

  7. [7]

    Vittoria Struglia

    Giulia Iori, Enzo Marinari, Giorgio Parisi, and M. Vittoria Struglia. Statistical mechanics of heteropolymer folding. Physica A: Statistical Mechanics and its Applications , 185(1):98–103, 6 This work is shared under an MIT License unless otherwise noted

    work page

  8. [8]

    doi:https://doi.org/10.1016/0378-4371(92)90442-S

    ISSN 0378-4371. doi:https://doi.org/10.1016/0378-4371(92)90442-S. URL https: //www.sciencedirect.com/science/article/pii/037843719290442S

    work page doi:10.1016/0378-4371(92)90442-s

  9. [9]

    John Jumper, Richard Evans, Alexander Pritzel, Tim Green, Michael Figurnov, Olaf Ronneberger, Kathryn Tunyasuvunakool, Russ Bates, Augustin ˇZ´ ıdek, Anna Potapenko, Alex Bridgland, Clemens Meyer, Simon A. A. Kohl, Andrew J. Ballard, Andrew Cowie, Bernardino Romera-Paredes, Stanislav Nikolov, Rishub Jain, Jonas Adler, Trevor Back, Stig Petersen, David Rei...

    work page doi:10.1038/s41586-021-03819-2 2021

  10. [10]

    Simulated tempering: A new monte carlo scheme

    E Marinari and G Parisi. Simulated tempering: A new monte carlo scheme. Europhysics Letters (EPL), 19(6):451–458, jul 1992. doi:10.1209/0295-5075/19/6/002. URL https://doi.org/10. 1209%2F0295-5075%2F19%2F6%2F002

    work page doi:10.1209/0295-5075/19/6/002 1992

  11. [11]

    Shakhnovich and A.M

    E.I. Shakhnovich and A.M. Gutin. Formation of microdomains in a quenched disordered heteropoly- mer. Journal de Physique , 50(14):1843–1850, 1989. doi:10.1051/jphys:0198900500140184300. URL https://hal.archives-ouvertes.fr/jpa-00211034

    work page doi:10.1051/jphys:0198900500140184300 1989

  12. [12]

    Monte carlo solution of structural dynamics

    Masanobu Shinozuka. Monte carlo solution of structural dynamics. Computers & Structures , 2(5): 855–874, 1972. ISSN 0045-7949. doi:https://doi.org/10.1016/0045-7949(72)90043-0. URL https://www.sciencedirect.com/science/article/pii/0045794972900430

    work page doi:10.1016/0045-7949(72)90043-0 1972

  13. [13]

    Portman, and Peter G

    Shoji Takada, John J. Portman, and Peter G. Wolynes. An elementary mode cou- pling theory of random heteropolymer dynamics. Proceedings of the Na- tional Academy of Sciences , 94(6):2318–2321, 1997. doi:10.1073/pnas.94.6.2318. https://www.pnas.org/doi/pdf/10.1073/pnas.94.6.2318. URL https://www.pnas.org/ doi/abs/10.1073/pnas.94.6.2318. 7

    work page doi:10.1073/pnas.94.6.2318 1997