REVIEW 5 major objections 5 minor 37 references
Modal Geometry Governs Proteoform Dynamics
T0 review · 5 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read This paper claims that conserved modal curvature, not raw energy alone, dictates which proteoform transitions happen and in what order, with least-action paths selected on a Boolean lattice.
desk verdict Ambitious but unsupported: the universal geometric law is a tautology and the headline regression is circular, though the GAPDH bitwise PDE weights are a concrete, testable seed. 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 modal manifold, the Boolean hypercube whose vertices are all possible proteoform states, equipped with a conserved scalar curvature field defined as a zero-sum Laplacian of occupancy and shape. This curvature is the mechanism that differentiates iso-potential transitions: in the MGF field equation, curvature acts as a geometric penalty, and least-action selection on the lattice chooses shortest Hamming paths such as 000 to 100 to 110 to 111. Supporting machinery includes the four axioms, PDE-derived bitwise weights from the atomic structure, and a bounded Ricci-flow degeneracy that quantifies oscillations between ordered and chaotic occupancy distributions.
What would settle it
Watch single GAPDH molecules as they oxidize: the model predicts Cys152 is modified first and the fully oxidized 111 state is reached mainly through 000 to 100 to 110 to 111; if a different first cysteine or a different dominant route is observed, the theory's core claim fails.
Extended reading notes
Core claim
MGF Theory asserts that proteoform dynamics is governed by conserved modal curvature. On the modal manifold of all possible modification states, each mode is one configuration, e.g., the eight cysteine redox states of GAPDH from 000 to 111. Axioms one through three fix the lattice structure, volume conservation, and first-order Hamming-1 transitions. Axiom four defines curvature as a zero-sum Laplacian of occupancy and shape, so total curvature always vanishes. The MGF field equation adds a curvature penalty to the energy cost and entropy gain of a transition: concentrated occupancy creates wells that suppress transitions, while distributed occupancy flattens the landscape and facilitates mo
Load-bearing premise
The argument stands on the premise that the curvature computed from a molecule's occupancy and shape is a physical agent that changes transition rates, rather than merely a mathematical summary of the occupancy itself.
Editorial extensions
If this is right
- If the theory is right, activation-energy barriers in proteoform networks are not fixed: they shift as occupancy redistributes curvature, producing hysteresis and path-dependent histories.
- For GAPDH specifically, the theory predicts Cys152 oxidizes first and full oxidation traverses 000 to 100 to 110 to 111, giving the hyperoxidized 111 mode a definite preferred route.
- If curvature governs transitions, moves with equal energy cost are distinguished geometrically, so no extra energy scale is needed to explain why one modification site is favored over another.
- The same axioms apply to any number of modification sites and any post-translational modification basis, so preferred geodesics should exist for tyrosine phosphorylation, lysine acetylation, and other proteoform systems.
Reading between the lines
- My extension: the total-curvature identity sum zero holds for any occupancy because the graph Laplacian has zero-sum rows; the theory's physical content therefore depends on whether the curvature field itself, rather than the occupancy it is derived from, changes transition rates.
- My extension: the Poisson and Dirichlet-energy recipe can be treated as a general predictor of oxidation order; applying it to other multi-cysteine proteins with known experimental oxidation hierarchies would either generalize or localize the GAPDH result.
- My extension: the predicted Ricci-flow oscillations imply proteoform populations should show wave-like spreading in time, so time-resolved measurements of modification-state distributions could look for periodic ordering and disordering of occupancy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes 'Modal Geometric Field (MGF) Theory' for proteoform dynamics. Four axioms define the proteoform state space as a Boolean lattice (Axiom 1), impose volume invariance (Axiom 2), restrict transitions to Hamming-1 moves (Axiom 3), and introduce a graph-Laplacian 'curvature' field that is claimed to be a conserved geometric quantity controlling transitions (Axiom 4). The central field equation couples an energy term, this curvature, and an entropy term. The manuscript reports simulations of the resulting action integral, claims that curvature and anisotropy significantly 'shape' the action, analyzes geodesics and commutator algebras, and presents a GAPDH case study in which PDE-derived bitwise weights predict Cys152 as the primary oxidation site and give an activation energy near an empirical value. The paper concludes that modal geometry governs proteoform dynamics universally and scale-invariantly.
Significance. If the central claims were correct, the paper would offer a new, universal geometric framework for proteoform dynamics, with concrete predictions such as the preferred GAPDH oxidation pathway 000 → 100 → 110 → 111. Strengths include the clear presentation of the Boolean/combinatorial state space, the availability of the Julia source code, and the attempt to connect a geometric model to a specific experimental system. However, the load-bearing evidence is not independent: the 'conservation of curvature' is a mathematical identity of the graph Laplacian, the regression showing that curvature predicts the action is circular because the action includes the same curvature and anisotropy terms as additive components, and the GAPDH activation-energy comparison relies on an unconstrained coupling constant. At present the manuscript does not establish the claimed geometric governance; the central law is asserted and then fitted rather than tested.
major comments (5)
- [The MGF Theory Field Equation (p. 5–6)] The field equation T_{i→j} = ΔE_{i→j}·exp[ρ(x_i)·ΔH1] − R(x_i) + ΔS is dimensionally inconsistent. ΔE is an energy, ΔS is an entropy (or an energy only if multiplied by T), and R(x_i) is defined in Axiom 4 as a graph Laplacian of a dimensionless occupancy/shape field, hence dimensionless. The three terms cannot be added in a physical equation. Moreover, the stationary condition ∑_{k=0}^R P(i→j)·T_{i→j}=0 is ill-defined: the summation index k does not appear in the summand, and T is not a probability. Because this equation is the central law of the paper, the manuscript does not provide a testable physical equation.
- [Conservation of Curvature (p. 10; Supplemental Axiom 4)] The 'conservation of curvature' is a tautology. With R = Lρ for the graph Laplacian L, the identity 1^T L = 0 implies ∑_x R(x) = 0 for every vector ρ at every instant. The text's own derivation, ∑ R = 1^T L ρ = 0, makes this explicit. This is the zero-sum property of the Laplacian, not a dynamical conservation law, and it does not depend on Axiom 2 or on any symmetry. The statement that 'curvature cannot be created or destroyed, only transported' therefore carries no physical content beyond the definition of R. The Noether analogy and the abstract's claim that conserved curvature governs dynamics are not established.
- [Modelling MGF Theory / Methods: Action integral (p. 6, p. 14)] The regression in Figure 2A is circular. The response is defined as 𝒜 = Σ_s[ΔS(s) + λ_R R(x_s)² + λ_A ||A(x_s)||²] in the main text and as 𝒜 = Σ_t[α_mass Σ|Δρ| + α_geom(⟨R²⟩ + ⟨||A||²⟩) + α_entropy S_deg] in Methods. The predictors in the regression are R and A, which appear additively inside the response. For any trajectory ensemble with nonzero variance in R or A, the general linear model must return nonzero coefficients; the sign of each coefficient is determined by covariance with the other additive terms. Figure 2A therefore cannot support the conclusion that 'geometry shaped proteoform dynamics'; it is a property of the estimator, not an empirical finding.
- [Curvature, Energy, and Entropy / Numerical example (p. 12–13)] The GAPDH activation-energy comparison is not a parameter-free prediction. The text states a symmetric barrier of ≈33 kJ/mol per cysteine bit, then defines ΔG_eff^‡(b) = κ_modal |K|(1 − ln w_E(b)) with |K| ≈ 33 kJ/mol. For Cys152 (w_E = 0.388) this gives roughly 64 kJ/mol, not the claimed 76 kJ/mol. To obtain 76 kJ/mol one must choose κ_modal ≈ 1.18, an unstated free parameter. Since κ_modal is not determined independently, the reported agreement with the empirical 67 kJ/mol is a fit, not a validation of the theory.
- [The Geometric Primitive (p. 11–12)] The geometric-primitive calculation depends on several arbitrary choices: Gaussian kernel width σ = 4.0 Å, cutoff r_c = 6.5 Å, Cα-only representation, the specific construction of the substrate Laplacian, the choice of source terms, and rescaling of Forman-Ricci curvature. No sensitivity analysis is provided. This calculation is the only molecular-level input that breaks the iso-potential symmetry of the three cysteine bits, so the prediction that Cys152 is the geometric primitive (000 → 100) is load-bearing. Without a demonstration that the result is robust to these modeling choices, the prediction may be an artifact of the chosen parameters.
minor comments (5)
- [Axiom 3 vs. simulations] Axiom 3 defines single-molecule occupancy as ρ(x)=1 at the occupied mode and 0 elsewhere, but the simulations evolve 100 molecules per run and use fractional occupancies. The relationship between the single-molecule axiom and the ensemble-level field equation needs clarification.
- [Field equation notation] The underbrace 'Energy' covers only ΔE·exp[...], while −R(x_i) and +ΔS are also treated as energy-like terms. The notation suggests that R and ΔS are not energy contributions, which is inconsistent with the dimensional issue noted above.
- [Ricci Flow definitions] Two definitions of Φ(ρ) are given: Φ(ρ)=max_{x,y}|R(x)−R(y)| and Φ(ρ)=∑(ρ_x−1/|ℳ|)². These are not equivalent. The text calls them 'two equivalent forms' (p. 15), which is incorrect.
- [Scale-invariance and universality] The scale-invariance test consists of running the same code for different R values. This does not test universality across different PTM types or protein structures, nor does it establish that the same κ_modal applies across systems.
- [References and typos] The GAPDH empirical oxidation claim cites references [22–24], but the mapping from text citations to the reference list is inconsistent (e.g., references 24–27 are about unrelated topics). There are also typos: 'Erying' for Eyring, 'In' for ln, and 'communicative' for commutative.
Circularity Check
Conservation is definitional, the regression evidence is tautological, and the GAPDH activation energy is calibrated by a free coupling constant.
-
self definitional
[Axiom 4, Results: 'Occupancy Curves the Modal Manifold'; Methods: 'Conservation of Curvature']
"Formally, if ℱ: (𝜌, 𝜑) ↦ 𝑅 denotes this projection, then curvature is 𝑅(𝑥) = (𝐿𝜑𝜎)(𝑥), ∑ 𝑅(𝑥) = 0 ... where 𝜎 𝐺(𝜌, 𝜑) encodes occupancy and shape, and 𝐿𝜑 is a symmetric zero-sum Laplacian. Hence, curvature cannot be created or destroyed: it is the conserved geometric expression of molecular identity"
A zero-sum Laplacian satisfies 1^T L = 0 algebraically, so ∑ R(x) = 0 for every σ by construction. The claimed Noether-like conservation law is an identity of the definition, not a derived physical conservation. The Methods section even states 'Since 𝐿⊺1=0, total curvature vanishes for any occupancy' and then calls this a 'Noether-like discrete analogue.' Thus the central assertion that curvature is 'always conserved' is imposed by the chosen representation rather than demonstrated.
-
self definitional
[Results: 'Modelling MGF Theory'; Methods: 'Action integral']
"To model the MGF equation, the stepwise action integral (𝒜) was defined as: 𝒜 = ∑[Δ𝑆(𝑠) + 𝜆𝑅𝑅(𝑥𝑠)2 + 𝜆𝐴‖𝐴(𝑥𝑠)‖2] ... Regression analysis revealed that geometry shaped proteoform dynamics (Figure 2A). High curvature suppressed 𝒜 (β = –1059.7, p < 1e–94), whereas directional anisotropy strongly amplified it (β = +2112.8, p < 1e–88). Hence, modal geometry governs proteoform dynamics."
The action integral is defined to include squared curvature and squared anisotropy as additive components. Regressing 𝒜 on R and A then simply recovers the same terms that were put into 𝒜. The significant β coefficients are forced by construction; they cannot provide independent evidence that geometry 'shaped' dynamics. The empirical content reduces to the arbitrary λ/α weights chosen by the authors, so the conclusion is a tautology of the estimator rather than a discovery.
1 more flagged steps
-
fitted input called prediction
[Results: 'Curvature, Energy, and Entropy'; Methods: 'Curvature & Energy']
"∆𝐺𝑐𝑢𝑟𝑣‡ = 𝜅𝑚𝑜𝑑𝑎𝑙|𝐾𝑖𝑗|, Where 𝜅𝑚𝑜𝑑𝑎𝑙 is the geometric coupling constant (𝐽 ⋅ 𝑚𝑜𝑙−1). ... Numerical example: ... occupancy induced curvature produced a symmetric barrier height of ≈ 33 𝑘𝐽⋅𝑚𝑜𝑙−1 per cysteine bit. By construction, this value is the same for all three bits, since the curvature spectrum is invariant under axiom 4. ... Applied to Cys152, this equation yielded an activation energy of 76 𝑘𝐽⋅𝑚𝑜𝑙−1, which is congruent with 67 𝑘𝐽⋅𝑚𝑜𝑙−1 as computed from empirical kinetic data31."
κ_modal is a free coupling constant; the paper provides no independent determination of it. The baseline barrier is explicitly stated to be fixed 'by construction', and the final 76 kJ/mol value is then called 'congruent' with the empirical 67 kJ/mol. Whether or not the numbers happen to agree, this is a calibration check, not a prediction: the energy scale and baseline are set by an unspecified constant, and the text even contains an internal inconsistency (Cys152's w_E gives ~64 kJ/mol, not 76, under the stated 33 kJ/mol baseline).
full rationale
The paper's central law—'curvature, which is always conserved, governs proteoform dynamics'—rests on two constructional identities and one calibration. First, curvature is defined as the output of a zero-sum Laplacian, so its global conservation is a mathematical identity (1^T L = 0), not an empirical law. Second, the simulation evidence for 'geometry shaped dynamics' comes from regressing the action integral on curvature and anisotropy, but the action integral was defined to include squared curvature and squared anisotropy terms; the regression merely rediscovers the construction. Third, the GAPDH activation-energy agreement relies on an unspecified coupling constant κ_modal, with the baseline barrier described as fixed 'by construction'; the 'prediction' is thus calibrated, not independent. The PDE-based bit weights and the empirically known oxidation-proneness of Cys152 retain some independent content, but they do not salvage the universal claim, which is not supported by any non-circular test. The paper's self-citations are not the main problem; the definitions and free parameters are.
Assumptions & free parameters
free parameters (4)
- kappa_modal =
not stated; implied by symmetric barrier of about 33 kJ/mol
- alpha_mass, alpha_geom, alpha_entropy =
0.01, 0.1, 0.1
- Gaussian kernel parameters sigma and r_c =
sigma = 4.0 Angstrom, r_c = 6.5 Angstrom
- symmetric barrier height |K| per cysteine =
about 33 kJ/mol
assumptions (6)
- standard math Boolean lattice {0,1}^R enumerates all proteoform states
- standard math Probability mass over modes is conserved, sum rho = 1
- domain assumption Only Hamming distance 1 transitions are allowed
- ad hoc to paper Zero-sum graph Laplacian curvature is a conserved physical field that governs transitions
- ad hoc to paper The modal manifold is non-symplectic because tangent spaces are {0}
- ad hoc to paper Scale invariance and universality across all R and all PTM types
invented entities (3)
-
Modal curvature R(x)
-
Anisotropic fibre bundle A(x)
-
Occupancy-induced metric deformation g = e^(2 sigma) g0
Cite this review
Pith. "Pith review of Modal Geometry Governs Proteoform Dynamics." pith.science (2026). https://pith.science/paper/BQ4TMC3B
@misc{pith2026250820004,
author = {Pith},
title = {Pith review of: Modal Geometry Governs Proteoform Dynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/BQ4TMC3B}},
note = {Machine review of arXiv:2508.20004}
}
read the original abstract
The fundamental laws governing proteoform dynamics have yet to be formulated. As a result, it is unclear how a specific proteoform, a distinct molecular variant of a protein, dynamically shapes its own future by evolving into new modes that exist only in potential until realised. Here, Modal Geometric Field (MGF) Theory couples real and abstract proteoform transitions through four axioms. Axioms 1 to 3 (invariant) dictate that only first-order transitions occur on the discrete, volume-invariant, non symplectic modal manifold. Axiom 4 (mutable) projects the occupancy and shape of a real, instantiated molecule into the modal manifold, generating occupancy-induced curvature. By coupling what is real to what is abstract, curvature, which is always conserved, governs proteoform dynamics by dictating the least-action modal transition. Because curvature distribution renders activation energy relative, barriers are mutable, and entropy emerges inevitably from curvature transport. This unification of energy, entropy, and curvature yields hysteresis, path dependence, fractal self similarity, and trajectories that oscillate between order and chaos. As a scale invariant and universal framework, MGF Theory reveals how modal geometry governs proteoform dynamics
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589 (2021)
2021
-
[2]
Smith, L. M. et al. Proteoform: a single term describing protein complexity. Nat Methods 10, 186–187 (2013)
work page 2013
-
[3]
Smith, L. M. & Kelleher, N. L. Proteoforms as the next proteomics currency. Science 359, 1106–1107 (2018)
work page 2018
-
[4]
Carbonara, K., Andonovski, M. & Coorssen, J. R. Proteomes Are of Proteoforms: Embracing the Complexity. Proteomes 9, 38 (2021)
work page 2021
-
[5]
Burnum-Johnson, K. E. et al. New Views of Old Proteins: Clarifying the Enigmatic Proteome. Mol Cell Proteomics 21, 100254 (2022)
work page 2022
-
[6]
Aebersold, R. et al. How many human proteoforms are there? Nat Chem Biol 14, 206–214 (2018)
work page 2018
-
[7]
Cobley, J. N., Chatzinikolaou, P. N. & Schmidt, C. A. The Nonlinear Cysteine Redox Dynamics in the i-Space: A Proteoform-centric Theory of Redox Regulation. Redox Biol. 103523 (2025) doi:10.1016/j.redox.2025.103523
-
[8]
Huang, K. Statistical Mechanics. vol. 2nd Edition (Wiley-VCH, 1987)
work page 1987
Show all 37 references
-
[9]
Die Grundlage der allgemeinen Relativitätstheorie
Einstein, A. Die Grundlage der allgemeinen Relativitätstheorie. Ann Phys 49, 769–822 (1916)
1916
-
[10]
Cobley, J. N. Oxiforms: Unique cysteine residue‐ and chemotype‐specified chemical combinations can produce functionally‐distinct proteoforms. Bioessays 45, (2023). 20
2023
-
[11]
Cobley, J. N. Exploring the unmapped cysteine redox proteoform landscape. Am. J. Physiol.-Cell Physiol. (2024) doi:10.1152/ajpcell.00152.2024
2024
-
[12]
Stanley, R. P. Enumerative Combinatorics. (2011). doi:10.1017/cbo9781139058520
2011 doi
-
[13]
Parvez, S., Long, M. J. C., Poganik, J. R. & Aye, Y. Redox Signaling by Reactive Electrophiles and Oxidants. Chem Rev 118, 8798–8888 (2018)
2018
-
[14]
C., Keske, J
Moser, C. C., Keske, J. M., Warncke, K., Farid, R. S. & Dutton, P. L. Nature of biological electron transfer. Nature 355, 796–802 (1992)
1992
-
[15]
Cobley, J. N. Information, Geometry, and Chaos: Revealing Latent Cysteine Butterflies on Fractal Redox Shapes in the Proteomic Spectra. Preprints doi:10.20944/preprints202507.2120.v1
- [16]
-
[17]
Choudhary, D., Foster, K. R. & Uphoff, S. Chaos in a bacterial stress response. Curr. Biol. 33, 5404-5414.e9 (2023)
2023
-
[18]
Chaos: Making the New Science
Gleick, J. Chaos: Making the New Science. (Penguin Books, 2008)
2008
-
[19]
LYAPUNOV, A. M. The general problem of the stability of motion. Int. J. Control 55, 531–534 (1992)
1992
-
[20]
B., Swinney, H
Wolf, A., Swift, J. B., Swinney, H. L. & Vastano, J. A. Determining Lyapunov exponents from a time series. Phys. D: Nonlinear Phenom. 16, 285–317 (1985)
1985
- [21]
-
[22]
Mandelbrot, B. B. Fractals in physics: Squig clusters, diffusions, fractal measures, and the unicity of fractal dimensionality. J. Stat. Phys. 34, 895–930 (1984)
1984
- [23]
-
[24]
R., Davies, M
Glover, M. R., Davies, M. J. & Fuentes-Lemus, E. Oxidation of the active site cysteine residue of glyceraldehyde-3-phosphate dehydrogenase to the hyper-oxidized sulfonic acid form is favored under crowded conditions. Free Radic. Biol. Med. (2023) doi:10.1016/j.freeradbiomed.20...
2023 doi
-
[25]
Talwar, D. et al. The GAPDH redox switch safeguards reductive capacity and enables survival of stressed tumour cells. Nat Metabolism 1–17 (2023) doi:10.1038/s42255-023- 00781-3
2023 doi
-
[26]
Ralser, M. et al. Dynamic rerouting of the carbohydrate flux is key to counteracting oxidative stress. J Biology 6, 10 (2007)
2007
-
[27]
Peralta, D. et al. A proton relay enhances H2O2 sensitivity of GAPDH to facilitate metabolic adaptation. Nat Chem Biol 11, 156–163 (2015)
2015
-
[28]
von et al
Montfort, C. von et al. The role of GAPDH in the selective toxicity of CNP in melanoma cells. PLOS ONE 19, e0300718 (2024)
2024
-
[29]
Brown, G. C. Bioenergetic myths of energy transduction in eukaryotic cells. Front. Mol. Biosci. 11, 1402910 (2024)
2024
-
[30]
Cobley, J. N. 50 shades of oxidative stress: A state-specific cysteine redox pattern hypothesis. Redox Biol. 67, 102936 (2023)
2023
-
[31]
The Activated Complex in Chemical Reactions
Eyring, H. The Activated Complex in Chemical Reactions. J. Chem. Phys. 3, 107–115 (1935)
1935
-
[32]
C., Peskin, A
Winterbourn, C. C., Peskin, A. V., Kleffmann, T., Radi, R. & Pace, P. E. Carbon dioxide/bicarbonate is required for sensitive inactivation of mammalian glyceraldehyde-3- phosphate dehydrogenase by hydrogen peroxide. Proc National Acad Sci 120, e2221047120 (2023)
2023
-
[33]
Roesch, E. et al. Julia for biologists. Nat Methods 1–10 (2023) doi:10.1038/s41592-023- 01832-z. 21
2023 doi
-
[34]
& Ulam, S
Metropolis, N. & Ulam, S. The Monte Carlo Method. J. Am. Stat. Assoc. 44, 335–341 (1949)
1949
-
[35]
Kalos, M. H. & Whitlock, P. A. Monte Carlo Methods. (2022) doi:10.1002/9783527626212. Supplemental Information Axiom 1 Extended to Multinomial Theorem. The structure of the entire proteoform state space ( ΩP)—all theoretically possible molecular variants from a single copy of ...
2022 doi
-
[36]
Horizontal Plane – Encodes PTM state-space structure using multinomial coefficients, stratifying proteoforms by their modification-level degeneracy
-
[37]
two-site jump
Vertical Plane – Encodes sequence-level divergence, placing near -identical PS (e.g., isoforms or mutations) closer together. This multinomial Pascal simplex that structures the entire proteoform landscape across sequence variability and PTM states. All proteoform modes—whethe...
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.