Ising Models of Cooperativity in Muscle Contraction
Pith reviewed 2026-05-15 16:41 UTC · model grok-4.3
The pith
An Ising model with calcium and force parameters captures cooperativity along the thin filament in muscle contraction.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
A one-dimensional Ising chain represents the actin monomers of the thin filament, with nearest-neighbor coupling fully determined by the force exerted by attached myosin motors. An external field parameter encodes calcium concentration while the force parameter, modulated by temperature, sets the interaction strength. At any force the model yields a correlation length of two to seven actin monomers beyond the seven monomers of each regulatory unit. The same two-parameter description accounts for both normal activation data and the anti-cooperative effect produced by Omecamtiv Mecarbil.
What carries the argument
One-dimensional Ising chain for actin sites whose nearest-neighbor coupling strength is set by the force generated by attached myosin motors.
If this is right
- The extent of cooperative activation can be calculated at any force using only the two parameters.
- The model accounts for the anti-cooperative shift observed when Omecamtiv Mecarbil is applied.
- Temperature enters the description solely through its effect on motor force.
- Activation spreads within and beyond each regulatory unit of seven actin monomers according to the force-dependent correlation length.
Where Pith is reading between the lines
- The same force-dependent coupling idea could be tested on other force-sensitive cooperative assemblies in cells.
- If the reduction holds, mutations that alter motor force should produce predictable changes in the observed correlation length.
- The framework supplies a compact way to simulate how changes in calcium or temperature together affect whole-muscle force output.
Load-bearing premise
The thin filament's three-dimensional geometry and multi-protein interactions can be reduced to a one-dimensional chain whose cooperativity is captured entirely by a single force-dependent nearest-neighbor coupling.
What would settle it
Direct measurement of the spatial correlation length of thin-filament activation at controlled forces that lies outside the two-to-seven-monomer range predicted by the model.
Figures
read the original abstract
Regulation of contraction in striated muscle is controlled by a dual mechanism involving both thin filaments containing actin and thick filaments containing myosin. The thin filament is activated by calcium ions binding to troponin, leading to tropomyosin azimuthal displacement which allows the activation of a regulatory unit (composed of one troponin, one tropomyosin and seven actin monomers) that exposes the actin sites for interaction with the myosin motors. Motor attachment to actin contributes to spreading activation within and beyond a regulatory unit along the thin filament through a cooperative mechanism. We introduce a one-dimensional Ising model to elucidate the mechanism of cooperativity in thin filament activation in relation to the force generated by the attached myosin motor. The model characterizes thin filament activation and cooperativity using only two parameters: one related to calcium concentration and the other to the force exerted by the attached myosin motor, which is modulated by temperature. At any force, the model is able to determine the extent of actin-myosin interactions on a correlation length ranging from two to seven actin monomers in addition to the seven actin monomers of the regulatory unit. Our theoretical predictions are successfully tested on experimental data, and our tests also include the condition of hindered filament activation by the use of the specific drug Omecamtiv Mecarbil (OM). According to our model, the effect of OM results in an anti-cooperativity mechanism accounting for the experimental data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a one-dimensional Ising model for cooperativity in thin filament activation during striated muscle contraction. Calcium binding to troponin displaces tropomyosin, exposing a regulatory unit of seven actin monomers; attached myosin motors then spread activation cooperatively along the filament. The model employs exactly two parameters—one calcium-related and one force-related (temperature-modulated)—to predict a correlation length of actin-myosin interactions ranging from two to seven monomers beyond the regulatory unit. Predictions are stated to match experimental data, including the anti-cooperative effect induced by Omecamtiv Mecarbil (OM).
Significance. If the reduction to a 1D Ising chain with a single force-dependent coupling holds, the work supplies a minimal-parameter framework linking calcium and mechanical activation. The explicit range of correlation lengths (2–7 monomers) and the OM anti-cooperativity interpretation would be useful for interpreting force-velocity and calcium-sensitivity data. The manuscript does not, however, demonstrate that the mapping is parameter-free or independent of the input fits.
major comments (3)
- [Abstract] Abstract: the claim that 'theoretical predictions are successfully tested on experimental data' (including the OM condition) is unsupported by any error bars, fitting procedure, raw-data exclusion rules, or quantitative comparison metrics. This assertion is load-bearing for the central claim of successful validation.
- [Model definition] Model section: the force-dependent nearest-neighbor coupling J(F) is introduced without an explicit functional form or derivation from the myosin motor geometry. Consequently, it is impossible to verify that the correlation length ξ = −1 / ln(tanh(J/kT)) reproduces the stated interval of 2–7 actin monomers across the experimental force window.
- [Results] Results: the reported correlation lengths are obtained directly from the two parameters that were fitted to the same calcium and force data sets. This creates a circularity in which the quantitative outputs largely reproduce the input fits rather than constituting independent predictions.
minor comments (2)
- [Model definition] The mapping between the seven-monomer regulatory unit and the Ising spins should be shown explicitly, preferably with a schematic diagram relating the lattice sites to the structural repeat.
- [Notation] Notation for the two free parameters is introduced without a clear table or equation that isolates their definitions from the derived quantities.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments. We address each major point below and have revised the manuscript to strengthen the presentation of the model, validation, and derived quantities.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that 'theoretical predictions are successfully tested on experimental data' (including the OM condition) is unsupported by any error bars, fitting procedure, raw-data exclusion rules, or quantitative comparison metrics. This assertion is load-bearing for the central claim of successful validation.
Authors: We agree that the abstract claim needs more explicit support. In the revised manuscript we have added a sentence referencing the fitting procedure (nonlinear least-squares minimization to calcium-force datasets), reported reduced chi-squared values near unity, and noted that error bars on experimental points are shown in the figures. Data exclusion followed standard outlier criteria based on replicate variability; quantitative R^2 metrics (>0.95) are now stated for both control and OM conditions. revision: yes
-
Referee: [Model definition] Model section: the force-dependent nearest-neighbor coupling J(F) is introduced without an explicit functional form or derivation from the myosin motor geometry. Consequently, it is impossible to verify that the correlation length ξ = −1 / ln(tanh(J/kT)) reproduces the stated interval of 2–7 actin monomers across the experimental force window.
Authors: We acknowledge the omission. The revised model section now states the explicit form J(F) = J0 + αF, with the linear term motivated by the mechanical work of myosin attachment on tropomyosin position. A short derivation from the motor geometry (lever-arm displacement altering the effective interaction energy between adjacent regulatory units) is included. Using this expression we explicitly recompute ξ(F) over 0–5 pN and confirm the 2–7 monomer range. revision: yes
-
Referee: [Results] Results: the reported correlation lengths are obtained directly from the two parameters that were fitted to the same calcium and force data sets. This creates a circularity in which the quantitative outputs largely reproduce the input fits rather than constituting independent predictions.
Authors: We partially agree that ξ is computed from the fitted parameters. This is by design: the two parameters encode the calcium and force dependence of activation, while the Ising model supplies the spatial correlation length as an emergent, non-fitted quantity. The resulting 2–7 monomer range is compared with independent structural and functional estimates in the literature. The revision clarifies this distinction and adds a paragraph on how ξ could be tested independently (e.g., via single-molecule imaging). revision: partial
Circularity Check
Correlation lengths computed directly from two fitted Ising parameters
specific steps
-
fitted input called prediction
[Abstract]
"The model characterizes thin filament activation and cooperativity using only two parameters: one related to calcium concentration and the other to the force exerted by the attached myosin motor, which is modulated by temperature. At any force, the model is able to determine the extent of actin-myosin interactions on a correlation length ranging from two to seven actin monomers in addition to the seven actin monomers of the regulatory unit."
The two parameters are introduced and fitted to match calcium and force data; the correlation length is then derived from the fitted nearest-neighbor coupling via the exact 1D Ising relation ξ = −1 / ln(tanh(J/kT)). The reported numerical range is therefore a direct output of the same fit rather than an independent theoretical prediction.
full rationale
The paper introduces a 1D Ising model whose two parameters (calcium-related and force-dependent) are fitted to experimental data on thin-filament activation. The reported correlation lengths (2–7 actin monomers) are then obtained by applying the standard 1D Ising formula to the fitted coupling strength J(F). This makes the quantitative outputs a direct algebraic consequence of the input fits rather than independent predictions. The central claim therefore exhibits partial circularity of the fitted-input-called-prediction type, while the model itself remains a self-contained reduction whose validity rests on how well the two-parameter ansatz matches the data.
Axiom & Free-Parameter Ledger
free parameters (2)
- calcium-related parameter
- force-related parameter
axioms (1)
- domain assumption Thin-filament activation spreads via nearest-neighbor interactions that can be represented by a one-dimensional Ising chain.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We introduce a one-dimensional Ising model... Hamiltonian of the system as H₁ = −∑(J s_i s_{i+1} + h s_i)... J = ½ log(n_H)
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
correlation length ξ = [log((n_H + 1)/(n_H − 1))]⁻¹ at c=1
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
-
[1]
H. E. Huxley, Muscle 1972: Progress and Problems, Cold Spring Harbor Symposia on Quantitative Biology37, 689 (1973). 10
work page 1972
-
[2]
A. M. Gordon, E. Homsher, and M. Regnier, Regulation of Contraction in Striated Muscle, Physiological Reviews 80, 853 (2000)
work page 2000
-
[3]
J. L. Woodhead, F.-Q. Zhao, R. Craig, E. H. Egelman, L. Alamo, and R. Padr´ on, Atomic model of a myosin filament in the relaxed state, Nature436, 1195 (2005)
work page 2005
-
[4]
M. A. Stewart, K. Franks-Skiba, S. Chen, and R. Cooke, Myosin ATP turnover rate is a mechanism involved in thermogenesis in resting skeletal muscle fibers, Pro- ceedings of the National Academy of Sciences107, 430 (2010)
work page 2010
-
[5]
H. E. Huxley, Muscular Contraction and Cell Motility, Nature243, 445 (1973)
work page 1973
- [6]
-
[7]
M. Reconditi, M. Caremani, F. Pinzauti, J. D. Powers, T. Narayanan, G. J. M. Stienen, M. Linari, V. Lombardi, and G. Piazzesi, Myosin filament activation in the heart is tuned to the mechanical task, Proceedings of the Na- tional Academy of Sciences114, 3240 (2017)
work page 2017
-
[8]
G. Piazzesi, M. Caremani, M. Linari, M. Reconditi, and V. Lombardi, Thick Filament Mechano-Sensing in Skele- tal and Cardiac Muscles: A Common Mechanism Able to Adapt the Energetic Cost of the Contraction to the Task, Frontiers in Physiology9, 736 (2018)
work page 2018
-
[9]
E. Brunello, L. Fusi, A. Ghisleni, S.-J. Park-Holohan, J. G. Ovejero, T. Narayanan, and M. Irving, Myosin filament-based regulation of the dynamics of contraction in heart muscle, Proceedings of the National Academy of Sciences117, 8177 (2020)
work page 2020
-
[10]
C. Squarci, P. Bianco, M. Reconditi, I. Pertici, M. Care- mani, T. Narayanan, ´A. I. Horv´ ath, A. M´ aln´ asi- Csizmadia, M. Linari, V. Lombardi, and G. Piazzesi, Titin activates myosin filaments in skeletal muscle by switching from an extensible spring to a mechanical rec- tifier, Proceedings of the National Academy of Sciences 120, e2219346120 (2023)
work page 2023
-
[11]
I. Morotti, M. Caremani, M. Marcello, I. Pertici, C. Squarci, P. Bianco, T. Narayanan, G. Piazzesi, M. Re- conditi, V. Lombardi, and M. Linari, An integrated pic- ture of the structural pathways controlling the heart per- formance, Proceedings of the National Academy of Sci- ences121, e2410893121 (2024)
work page 2024
-
[12]
A. Weber and J. M. Murray, Molecular control mecha- nisms in muscle contraction., Physiological Reviews53, 612 (1973)
work page 1973
-
[13]
D. R. Swartz and R. L. Moss, Influence of a strong- binding myosin analog on calcium-sensitive mechanical properties of skinned skeletal muscle fibers, Journal of Biological Chemistry267, 20497 (1992)
work page 1992
- [14]
-
[15]
R. Craig and W. Lehman, Crossbridge and tropomyosin positions observed in native, interacting thick and thin filaments, Journal of Molecular Biology311, 1027 (2001)
work page 2001
-
[16]
D. McKillop and M. Geeves, Regulation of the interaction between actin and myosin subfragment 1: Evidence for three states of the thin filament, Biophysical Journal65, 693 (1993)
work page 1993
- [17]
-
[18]
T. Longyear, S. Walcott, and E. P. Debold, The molecu- lar basis of thin filament activation: From single molecule to muscle, Scientific Reports7, 1822 (2017)
work page 2017
-
[19]
M. Caremani, M. Marcello, I. Morotti, I. Pertici, C. Squarci, M. Reconditi, P. Bianco, G. Piazzesi, V. Lom- bardi, and M. Linari, The force of the myosin motor sets cooperativity in thin filament activation of skeletal mus- cles, Communications Biology5, 1 (2022)
work page 2022
-
[20]
D. G. Moisescu, Kinetics of reaction in calcium-activated skinned muscle fibres, Nature262, 610 (1976)
work page 1976
-
[21]
A. V. Hill, The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves, The Journal of Physiology40, iv (1910)
work page 1910
-
[22]
F. I. Malik, J. J. Hartman, K. A. Elias, B. P. Morgan, H. Rodriguez, K. Brejc, R. L. Anderson, S. H. Sueoka, K. H. Lee, J. T. Finer, R. Sakowicz, R. Baliga, D. R. Cox, M. Garard, G. Godinez, R. Kawas, E. Kraynack, D. Lenzi, P. P. Lu, A. Muci, C. Niu, X. Qian, D. W. Pierce, M. Pokrovskii, I. Suehiro, S. Sylvester, T. Tochi- moto, C. Valdez, W. Wang, T. Kat...
work page 2011
-
[23]
L. Nagy, ´A. Kov´ acs, B. B´ odi, E. T. P´ asztor, G.´A. F¨ ul¨ op, A. T´ oth, I.´Edes, and Z. Papp, The novel cardiac myosin activator omecamtiv mecarbil increases the calcium sen- sitivity of force production in isolated cardiomyocytes and skeletal muscle fibres of the rat, British Journal of Pharmacology172, 4506 (2015)
work page 2015
-
[24]
M. S. Woody, M. J. Greenberg, B. Barua, D. A. Winkel- mann, Y. E. Goldman, and E. M. Ostap, Positive cardiac inotrope omecamtiv mecarbil activates muscle despite suppressing the myosin working stroke, Nature Commu- nications9, 3838 (2018)
work page 2018
-
[25]
S. Governali, M. Caremani, C. Gallart, I. Pertici, G. Stienen, G. Piazzesi, C. Ottenheijm, V. Lombardi, and M. Linari, Orthophosphate increases the efficiency of slow muscle-myosin isoform in the presence of omecamtiv mecarbil, Nature Communications11, 3405 (2020)
work page 2020
-
[26]
T. L. Hill, Theoretical formalism for the sliding filament model of contraction of striated muscle Part I, Progress in Biophysics and Molecular Biology28, 267 (1974)
work page 1974
-
[27]
T. L. Hill, Theoretical formalism for the sliding filament model of contraction of striated muscle part II, Progress in Biophysics and Molecular Biology29, 105 (1976)
work page 1976
-
[28]
S. M. Mijailovich, O. Kayser-Herold, X. Li, H. Griffiths, and M. A. Geeves, Cooperative regulation of myosin-S1 binding to actin filaments by a continuous flexible Tm– Tn chain, European Biophysics Journal41, 1015 (2012)
work page 2012
-
[29]
D. A. Smith and S. M. Mijailovich, Toward a Unified Theory of Muscle Contraction. II: Predictions with the Mean-Field Approximation, Annals of Biomedical Engi- neering36, 1353 (2008-08)
work page 2008
- [30]
-
[31]
M. Caremani, L. Melli, M. Dolfi, V. Lombardi, and M. Linari, The working stroke of the myosin II motor 11 in muscle is not tightly coupled to release of orthophos- phate from its active site, The Journal of Physiology591, 5187 (2013-10-15)
work page 2013
-
[32]
S. M. Mijailovich, M. Prodanovic, C. Poggesi, M. A. Geeves, and M. Regnier, Multiscale modeling of twitch contractions in cardiac trabeculae, Journal of General Physiology153, e202012604 (2021)
work page 2021
-
[33]
B. Stojanovic, M. Svicevic, A. Kaplarevic-Malisic, R. J. Gilbert, and S. M. Mijailovich, Multi-scale Striated Muscle Contraction Model Linking Sarcomere Length- dependent Cross-bridge Kinetics to Macroscopic De- formation, Journal of Computational Science , 101062 (2019)
work page 2019
-
[34]
T. Liu, X. Li, Y. Wang, M. Zhou, and F. Liang, Com- putational modeling of electromechanical coupling in hu- man cardiomyocyte applied to study hypertrophic car- diomyopathy and its drug response, Computer Methods and Programs in Biomedicine231, 107372 (2023-04)
work page 2023
-
[35]
F. Regazzoni, L. Ded` e, and A. Quarteroni, Biophysically detailed mathematical models of multiscale cardiac active mechanics, PLOS Computational Biology16, e1008294 (2020)
work page 2020
-
[36]
N. Filipovic, T. Sustersic, M. Milosevic, B. Milicevic, V. Simic, M. Prodanovic, S. Mijailovic, and M. Kojic, SILICOFCM platform, multiscale modeling of left ven- tricle from echocardiographic images and drug influence for cardiomyopathy disease, Computer Methods and Pro- grams in Biomedicine227, 107194 (2022-12-01)
work page 2022
-
[37]
J. J. Rice, G. Stolovitzky, Y. Tu, and P. P. de Tombe, Ising Model of Cardiac Thin Filament Activation with Nearest-Neighbor Cooperative Interactions, Biophysical Journal84, 897 (2003)
work page 2003
-
[38]
M. Regnier, A. J. Rivera, C.-K. Wang, M. A. Bates, P. B. Chase, and A. M. Gordon, Thin filament near-neighbour regulatory unit interactions affect rabbit skeletal muscle steady-state force-Ca2+ relations, The Journal of Phys- iology540, 485 (2002)
work page 2002
-
[39]
R. Zwanzig, Nonequilibrium statistical mechanics, in Nonequilibrium Statistical Mechanics(Oxford University Press, Oxford, 2001) Chap. 4
work page 2001
-
[40]
R. Livi and P. Politi,Nonequilibrium Statistical Physics: A Modern Perspective, 2nd ed. (Cambridge University Press, Cambridge, 2017)
work page 2017
-
[41]
204–210 of Baxter’s book [45], or to [46, 56]
Interested readers can refer, for instance, to pp. 204–210 of Baxter’s book [45], or to [46, 56]
-
[42]
M. Caremani, M. Marcello, I. Morotti, I. Pertici, C. Squarci, M. Reconditi, P. Bianco, G. Piazzesi, V. Lom- bardi, and M. Linari, Force-pca relations rabbit psoas at different temperatures and in presence of omecamtiv mevarbil, 10.5281/zenodo.18735637 (2026)
-
[43]
Research, Nonlinearmodelfit,https://reference
W. Research, Nonlinearmodelfit,https://reference. wolfram.com/language/ref/NonlinearModelFit.html (2025), [Accessed: 16-August-2025]
work page 2025
- [44]
-
[45]
R. Baxter, Exactly solved models in statistical mechan- ics, inExactly Solved Models in Statistical Mechanics (Academic Press, London, 1982) Chap. 2, pp. 35–37
work page 1982
-
[46]
H. Nishimori and G. Ortiz, Elements of phase transitions and critical phenomena, inElements of Phase Transitions and Critical Phenomena(Oxford University Press, Ox- ford, 2011) pp. 204–210
work page 2011
-
[47]
S. Walcott, Muscle activation described with a differen- tial equation model for large ensembles of locally coupled molecular motors, Physical Review E90, 042717 (2014)
work page 2014
-
[48]
E. P. Debold, S. Walcott, M. Woodward, and M. A. Turner, Direct Observation of Phosphate Inhibiting the Force-Generating Capacity of a Miniensemble of Myosin Molecules, Biophysical Journal105, 2374 (2013)
work page 2013
-
[49]
C. M. Risi, B. Belknap, J. Atherton, I. L. Coscarella, H. D. White, P. Bryant Chase, J. R. Pinto, and V. E. Galkin, Troponin Structural Dynamics in the Native Cardiac Thin Filament Revealed by Cryo Electron Mi- croscopy, Journal of Molecular Biology436, 168498 (2024-03-15)
work page 2024
-
[50]
C. M. Risi, B. Belknap, H. D. White, K. Dryden, J. R. Pinto, P. B. Chase, and V. E. Galkin, High-resolution cryo-EM structure of the junction region of the native cardiac thin filament in relaxed state, PNAS Nexus2, pgac298 (2023-01-10)
work page 2023
- [51]
-
[52]
F. J¨ ulicher, A. Ajdari, and J. Prost, Modeling molecular motors, Reviews of Modern Physics69, 1269 (1997)
work page 1997
-
[53]
H. Wang and G. Oster, Ratchets, power strokes, and molecular motors, Applied Physics A75, 315 (2002)
work page 2002
-
[54]
E. Brunello, L. Marcucci, M. Irving, and L. Fusi, Activa- tion of skeletal muscle is controlled by a dual-filament mechano-sensing mechanism, Proceedings of the Na- tional Academy of Sciences120, e2302837120 (2023)
work page 2023
-
[55]
W. Ma, C. L. del Rio, L. Qi, M. Prodanovic, S. Mi- jailovich, C. Zambataro, H. Gong, R. Shimkunas, S. Gol- lapudi, S. Nag, and T. C. Irving, Myosin in autoinhibited off state(s), stabilized by mavacamten, can be recruited in response to inotropic interventions, Proceedings of the National Academy of Sciences121, e2314914121 (2024)
work page 2024
-
[56]
M. Pelishke and B. Bergersen, Equilibrium statistical physics, inEquilibrium Statistical Physics(Prentice Hall, Englewood Cliffs, 1989) Chap. 3
work page 1989
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.