pith. sign in

arxiv: 2601.18073 · v2 · submitted 2026-01-26 · ⚛️ physics.bio-ph

Effects of stimulation frequencies on energy efficiency of a muscle fiber during contraction

Jiaxiang Xu , Bin Chen This is my paper

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

classification ⚛️ physics.bio-ph
keywords muscle contractionenergy efficiencystimulation frequencycross-bridge cycleshortening velocitybiophysical modelingmyosin kineticscalcium signaling
0
0 comments X

The pith

Emergent shortening velocity primarily determines cross-bridge efficiency in muscle fibers, peaking at an optimal speed across stimulation frequencies.

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

This paper creates a biophysical model that links calcium-driven excitation to the detailed steps of the cross-bridge cycle in a single muscle fiber. The simulations reveal that the speed at which the fiber shortens during contraction is the main factor controlling how efficiently the muscle converts chemical energy into mechanical work. Efficiency reaches its highest point at one particular shortening velocity and falls away when the fiber moves slower or faster than that speed. Stimulation frequency plays a lesser role, mainly by influencing what shortening velocity emerges. This approach helps explain why experiments have given conflicting results on frequency effects and points to velocity as the key control point for muscle energy use.

Core claim

Our model predictions indicate that the emergent shortening velocity is the primary determinant of cross-bridge efficiency: efficiency peaks at an optimal velocity and declines at higher or lower velocities, while frequency appears to exert secondary influence. Critically, the velocity yielding peak efficiency remains almost consistent across frequencies, with a slight upward shift at higher frequencies in most of our parametric studies. Interestingly, elevated inorganic phosphate ([Pi]) appears to amplify the efficiency disparity between high- and low-frequency regimes in our analysis. Our work suggests that stimulation frequency modulates efficiency predominantly through its regulation of

What carries the argument

Biophysical model integrating calcium-mediated excitation with a detailed cross-bridge cycle for single-fiber simulations, showing that shortening velocity governs efficiency through myosin power stroke kinetics.

If this is right

  • Efficiency peaks at an optimal shortening velocity that stays nearly the same across different stimulation frequencies.
  • Higher frequencies cause only a minor upward shift in the velocity that maximizes efficiency.
  • Raised levels of inorganic phosphate widen the efficiency difference seen between high-frequency and low-frequency stimulation.
  • Frequency controls efficiency indirectly by setting the shortening velocity rather than by changing cross-bridge steps directly.

Where Pith is reading between the lines

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

  • The single-fiber model could be tested in whole-muscle preparations where shortening velocity varies with external load.
  • In conditions with chronically elevated phosphate, the efficiency gap between stimulation rates may become more pronounced.
  • Training or rehabilitation strategies might target specific shortening velocities to improve energy use independent of neural firing rate.

Load-bearing premise

The detailed cross-bridge cycle and calcium-mediated excitation components accurately reproduce the observed velocity dependence of efficiency without additional regulatory mechanisms.

What would settle it

Directly measure efficiency at controlled shortening velocities in single fibers under different stimulation frequencies; if the peak efficiency velocity varies substantially with frequency, the primary role of velocity would be falsified.

read the original abstract

Contradictory experimental reports on the relationship between efficiency and stimulation frequency have hindered mechanistic understanding in converting neural activity into mechanical work during muscle contraction. To resolve this issue, we develop a biophysical model integrating calcium-mediated excitation with a detailed cross-bridge cycle to enable single-fiber simulations. Our model predictions indicate that the emergent shortening velocity is the primary determinant of cross-bridge efficiency: efficiency peaks at an optimal velocity and declines at higher or lower velocities, while frequency appears to exert secondary influence. Critically, the velocity yielding peak efficiency remains almost consistent across frequencies, with a slight upward shift at higher frequencies in most of our parametric studies. Interestingly, elevated inorganic phosphate ([Pi]) appears to amplify the efficiency disparity between high- and low-frequency regimes in our analysis. Our work suggests that stimulation frequency modulates efficiency predominantly through its regulation of shortening velocity, which primarily governs the kinetics of the myosin power stroke. This work may help clarify neural control of muscle energetics, and provide a quantitative foundation for studying muscle function in physiological and pathological contexts.

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 develops a biophysical model integrating calcium-mediated excitation with a detailed cross-bridge cycle to perform single-fiber simulations of muscle contraction. Model predictions indicate that the emergent shortening velocity is the primary determinant of cross-bridge efficiency, which peaks at an optimal velocity and declines at higher or lower velocities, while stimulation frequency exerts only secondary influence by modulating the velocity distribution. The velocity for peak efficiency remains largely consistent across frequencies (with slight upward shifts at higher frequencies in most parametric studies), and elevated inorganic phosphate amplifies efficiency differences between high- and low-frequency regimes.

Significance. If the integrated model faithfully captures the relevant kinetics, the work supplies a mechanistic account of how neural stimulation frequency affects muscle energy efficiency through velocity-dependent cross-bridge attachment/detachment rates. This offers a quantitative framework that could reconcile contradictory experimental reports on frequency-efficiency relationships and supports further study of muscle function in physiological and pathological settings. The parametric robustness checks across frequency and [Pi] are a clear strength.

major comments (2)
  1. [Abstract and Results] The central claim that velocity is the primary determinant rests entirely on simulation outputs; the manuscript provides no direct comparison of predicted efficiency values or velocity-efficiency curves to experimental measurements at matched stimulation frequencies, which limits the ability to substantiate resolution of contradictory reports (Abstract and Results sections).
  2. [Methods] Methods: the precise definition and calculation of efficiency (e.g., work output relative to ATP hydrolysis rate) must be shown explicitly to confirm it derives from independent kinetic rates rather than being partly constrained by the velocity parameters used in the Huxley-style cycle.
minor comments (2)
  1. [Methods] Add a table or supplementary section listing all parameter values, initial conditions, and numerical integration settings to ensure full reproducibility of the single-fiber simulations.
  2. [Figures] Figure legends should explicitly state the frequency range, [Pi] values, and number of parameter sweeps performed so readers can assess the scope of the reported robustness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and recommendation for minor revision. We address the major comments point by point below, with revisions where appropriate to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Abstract and Results] The central claim that velocity is the primary determinant rests entirely on simulation outputs; the manuscript provides no direct comparison of predicted efficiency values or velocity-efficiency curves to experimental measurements at matched stimulation frequencies, which limits the ability to substantiate resolution of contradictory reports (Abstract and Results sections).

    Authors: We agree that the absence of direct experimental comparisons at matched frequencies limits the strength of claims regarding resolution of contradictory reports. As a modeling study, our work derives predictions from an integrated biophysical model built on literature kinetic parameters rather than new experimental data. We will revise the Abstract and Results to explicitly frame the velocity-efficiency relationship as an emergent model prediction, clarify that frequency influences efficiency indirectly via velocity, and add discussion of how these predictions may inform future experiments to test reconciliation of experimental discrepancies. revision: partial

  2. Referee: [Methods] Methods: the precise definition and calculation of efficiency (e.g., work output relative to ATP hydrolysis rate) must be shown explicitly to confirm it derives from independent kinetic rates rather than being partly constrained by the velocity parameters used in the Huxley-style cycle.

    Authors: We appreciate this suggestion for improved clarity. We will expand the Methods section to include an explicit subsection defining efficiency as the ratio of mechanical work (force times shortening distance) to the cumulative ATP hydrolysis events counted from cross-bridge detachment transitions. This formulation will be shown to follow directly from the state-dependent kinetic rates in the Huxley-style cycle without additional velocity constraints imposed beyond the model's attachment/detachment probabilities. revision: yes

Circularity Check

0 steps flagged

No significant circularity: efficiency emerges from independent cross-bridge kinetics

full rationale

The derivation integrates calcium excitation with standard Huxley-style cross-bridge attachment/detachment rates to produce an emergent shortening velocity; efficiency is then computed directly from those rates and power-stroke energetics. No parameter is fitted to the target efficiency-velocity curve, no self-citation supplies a uniqueness theorem, and no ansatz is smuggled in. The reported velocity dependence is a direct, falsifiable output of the kinetic scheme rather than a renaming or redefinition of the inputs. Parametric sweeps over frequency and [Pi] further confirm the result is not forced by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract alone does not enumerate specific free parameters, axioms, or invented entities; typical muscle models contain many rate constants and scaling factors whose values are not stated here.

pith-pipeline@v0.9.0 · 5470 in / 1028 out tokens · 27718 ms · 2026-05-16T11:26:48.191944+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages

  1. [1]

    Bruce M. and M. D. Carlson. 2019. Chapter 5 The Muscular System. The Human Body, Academic Press. 111-136

    work page 2019

  2. [2]

    De Luca, C. J., R. S. LeFever, M. P. McCue, and A. P. Xenakis. 1982. Control scheme governing concurrently active human motor units during voluntary contractions. J. Physiol. 329:129–142

    work page 1982

  3. [3]

    J., and Z

    De Luca, C. J., and Z. Erim. 1994. Common drive of motor units in regulation of muscle force. Trends. Neurosci. 17(7):299–305

    work page 1994

  4. [4]

    M., and J

    Enoka, R. M., and J. Duchateau. 2017. Rate coding and the control of muscle force. Cold. Spring. Harb. Perspect. Med. 7(10):a029702

    work page 2017

  5. [5]

    Farina, D., A. M. Castronovo, I. Vujaklija, A. Sturma, S. Salminger, C. Hofer, and O. Aszmann. 2017. Common Synaptic Input to Motor Neurons and Neural Drive to Targeted Reinnervated Muscles. J. Neurosci. 37(46):11285–11292

    work page 2017

  6. [6]

    J., and N

    Barclay, C. J., and N. A. Curtin. 2023. Advances in understanding the energetics of muscle contraction. J. Biomech. 156:111669

    work page 2023

  7. [7]

    Russell, A. J., J. J. Hartman, A. C. Hinken, A. R. Muci, R. Kawas, L. Driscoll, G. Godinez, K. H. Lee, D. Marquez, W. F. 4th Browne, M. M. Chen, D. Clarke, S. E. Collibee, M. Garard, R. Hansen, Z. Jia, P. P. Lu, H. Rodriguez, K. G. Saikali, J. Schaletzky, and F. I. Malik. 2012. Activation of fast skeletal muscle troponin as a potential therapeutic approac...

    work page 2012

  8. [8]

    Hill, A. V . 1939. The mechanical efficiency of frog muscle. Proc. R. Soc. Lond. B. 127 (1939):434-451

    work page 1939

  9. [9]

    Barclay, C. J. 1996. Mechanical efficiency and fatigue of fast and slow muscles of the mouse. J. Physiol. 497 ( Pt 3)(Pt 3):81–794

    work page 1996

  10. [10]

    Efficiency

    Steiger, G. J., and J. C. Rüegg. 1969. Energetics and “Efficiency” in the isolated contractile machinery of an insect fibrillar muscle at various frequencies of oscillation. Pflugers. Arch. 307:1–21

    work page 1969

  11. [11]

    Buschman, H. P. J., G. Elzinga, and R. C. Woledge. 1996. The effects of the level of activation and shortening velocity on energy output in type 3 muscle fibres from Xenopus laevis. Pflugers. Arch. 433:153–159

    work page 1996

  12. [12]

    B., and C

    Lewis, D. B., and C. J. Barclay. 2014. Efficiency and cross -bridge work output of skeletal muscle is decreased at low levels of activation. PFLUG. ARCH. EUR. J. PHY. 466(3):599–609

    work page 2014

  13. [13]

    Abbate, F., C. J. De Ruiter, C. Offringa, A. J. Sargeant, and A. De Haan. 2002. In situ rat fast skeletal muscle is more efficient at submaximal than at maximal activation levels. J. Appl. Physiol. 92:2089–2096

    work page 2002

  14. [14]

    Fenn, W. O. 1923. A quantitative comparison between the energy liberated and the work performed by the isolated sartorius muscle of the frog. J. Physiol. 58(2- 3):175–203

    work page 1923

  15. [15]

    Nataraj, C., and S. Hemmati. 2025. Harnessing Nonlinear Mechanics to Transform Medical Diagnostics. ASME. Appl. Mech. Rev. 77(6):061001

    work page 2025

  16. [16]

    Energetics of muscle contraction: further trials

    Yamada, K., 2017. Energetics of muscle contraction: further trials. J. Physiol. Sci. 67(1):19–43

    work page 2017

  17. [17]

    Barclay, C. J. 2017. The basis of differences in thermodynamic efficiency among skeletal muscles. Clin. Exp. Pharmacol. Physiol. 44(12):1279-1286

    work page 2017

  18. [18]

    Tao, and B

    Xu, J., J. Tao, and B. Chen. 2025. Delayed phosphate release can highly improve energy efficiency of muscle contraction. Biophys. J. 124(14):2303–2316

    work page 2025

  19. [19]

    Shorten, P. R., P. O'Callaghan, J. B. Davidson, and T. K. Soboleva. 2007. A mathematical model of fatigue in skeletal muscle force contraction. J. Muscle Res. Cell Motil. 28(6):293–313

    work page 2007

  20. [20]

    Röhrle, O., J. B. Davidson, and A. J. Pullan. 2008. Bridging scales: A three - dimensional electromechanical finite element model of skeletal muscle. SIAM. J. Sci. Comput. 30(6):2882–2904

    work page 2008

  21. [21]

    Heidlauf, T., and O. Röhrle. 2014. A multiscale chemo-electro-mechanical skeletal muscle model to analyze muscle contraction and force generation for different muscle fiber arrangements. Front. Physiol. 5(498):1–14

    work page 2014

  22. [22]

    Somjen, and D

    Henneman, E., G. Somjen, and D. O. Carpenter. 1965. Excitability and inhibitability of motoneurons of different sizes. J. NEUROPHYSIOL. 28(3):599–620

    work page 1965

  23. [23]

    Fill, M., and J. A. Copello. 2002. Ryanodine receptor calcium release channels. Physiol. Rev. 82(4):893–922

    work page 2002

  24. [24]

    Rios, E., and G. Brum. 1987. Involvement of dihydropyridine receptors in excitation-contraction coupling in skeletal muscle. Nature. 325(6106):717–720

    work page 1987

  25. [25]

    Yan, and Z

    Wu, J., N. Yan, and Z. Yan. 2017. Structure -function relationship of the voltage - gated calcium channel cav1.1 complex. ADV . EXP . MED. BIOL. 981:23–39

    work page 2017

  26. [26]

    Block, B. A., T. Imagawa, K. P. Campbell, and C. Franzini -Armstrong. 1988. Structural evidence for direct interaction between the molecular components of the transverse tubule/sarcoplasmic reticulum junction in skeletal muscle. J. Cell. Biol. 107(6 Pt 2):2587–2600

    work page 1988

  27. [27]

    Otazu, G. H., R. Futami, and N. Hoshimiya. 2001. A muscle activation model of variable stimulation frequency response and stimulation history, based on positive feedback in calcium dynamics. Biol. Cybern. 84(3):193–206

    work page 2001

  28. [28]

    B., and D

    Cannell, M. B., and D. G. Allen. 1984. Model of calcium movements during activation in the sarcomere of frog skeletal muscle. Biophys. J. 45(5):913–925

    work page 1984

  29. [29]

    M., and S

    Baylor, S. M., and S. Hollingworth. 1998. Model of sarcomeric Ca2+ movements, including ATP Ca2+ binding and diffusion, during activation of frog skeletal muscle. J. Gen. Physiol. 112(3):297–316

    work page 1998

  30. [30]

    Winegrad, S. 1970. The intracellular site of calcium activaton of contraction in frog skeletal muscle. J. Gen. Physiol. 55(1):77–88

    work page 1970

  31. [31]

    Primeau, J. O., G. P. Armanious, M. E. Fisher, and H. S. Young. 2018. The sarcoendoplasmic reticulum calcium ATPase. Subcell. Biochem. 87:229–258

    work page 2018

  32. [32]

    Kurebayashi, A

    Ogawa, Y., N. Kurebayashi, A. Irimajiri, and T. Hanai. 1981. Transient kinetics for Ca uptake by fragmented sarcoplasmic reticulum from bullfrog skeletal muscle with reference to the rate of relaxation of living muscle. Molecular and Cellular Aspects of Muscle Function. 417-435

    work page 1981

  33. [33]

    Razumova, M. V., A. E. Bukatina, and K. B. Campbell. 1999. Stiffness -distortion sarcomere model for muscle simulation. J. APPL. PHYSIOL. 87(5):1861–1876

    work page 1999

  34. [34]

    Campbell, K. B., M. V. Razumova, R. D. Kirkpatrick, and B. K. Slinker. 2001. Myofilament kinetics in isometric twitch dynamics. Ann. Biomed. Eng. 29(5):384– 405

    work page 2001

  35. [35]

    Moult, and M

    Herzberg, O., J. Moult, and M. N. James. 1986. A model for the Ca2+ -induced conformational transition of troponin C. A trigger for muscle contraction. J. Biol. Chem. 261(6):2638–2644

    work page 1986

  36. [36]

    L., and A

    Nastuk, W. L., and A. L. Hodgkin. 1950. The electrical activity of single muscle fibers. J. Cell. Comp. Physiol. 35(1):39-73

    work page 1950

  37. [37]

    Saeki, and Y

    Kawai, M., Y . Saeki, and Y . Zhao. 1993. Crossbridge scheme and the kinetic constants of elementary steps deduced from chemically skinned papillary and trabecular muscles of the ferret. CIRC. RES. 73(1):35–50

    work page 1993

  38. [38]

    W., and E

    Lymn, R. W., and E. W. Taylor. 1971. Mechanism of adenosine triphosphate hydrolysis by actomyosin. BioChemistry. 10:4617-4624

    work page 1971

  39. [39]

    Brunello, M

    Linari, M., E. Brunello, M. Reconditi, L. Fusi, M. Caremani, T. Narayanan, G. Piazzesi, V. Lombardi, and M. Irving. 2015. Force generation by skeletal muscle is controlled by mechanosensing in myosin filaments. Nature. 528:276- 279

    work page 2015

  40. [40]

    Chen, B. 2013. Self-regulation of motor force through chemomechanical coupling in skeletal muscle contraction. J. Appl. Mech. 80:051013

    work page 2013

  41. [41]

    Brenner, B. 1991. Rapid dissociation and reassociation of actomyosin cross-bridges during force generation: a newly observed facet of cross -bridge action in muscle. Proc. Natl. Acad. Sci. U. S. A. 88:10490-10494

    work page 1991

  42. [42]

    Bell, G. I. 1978. Models for the specific adhesion of cells to cells. Science. 200:618- 627

    work page 1978

  43. [43]

    Mackerle, J. 1998. A finite element bibliography for biomechanics. ASME. Appl. Mech. Rev. 51(10):587–634

    work page 1998

  44. [44]

    Chen, B., and H. Gao. 2011. Motor force homeostasis in skeletal muscle contraction. Biophys. J. 101:396-403

    work page 2011

  45. [45]

    L., and B

    Dong, C. L., and B. Chen. 2015. Catch -slip bonds can be dispensable for motor force regulation during skeletal muscle contraction. Phys. Rev. E. Stat. Nonlin. Soft. Matter. Phys. 92:012723

    work page 2015

  46. [46]

    Gillespie, D. T. 1976. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J. Comput. Phys. 22(4):403-434

    work page 1976

  47. [47]

    Gillespie, D. T. 1977. Exact stochastic simulation of coupled chemical reactions. J. Phys. Chem. 81:2340-2361

    work page 1977

  48. [48]

    Yoneda, J

    Washio, T., K. Yoneda, J. Okada, T. Kariya, S. Sugiura, and T. Hisada. 2016. Ventricular fiber optimization utilizing the branching structure. Int. J. Numer. Method. Biomed. Eng. 32(7):e02753

    work page 2016

  49. [49]

    E., and W

    Huxley, H. E., and W. Brown. 1967. The low -angle x-ray diagram of vertebrate striated muscle and its behaviour during contraction and rigor. J. Mol. Biol. 30(2):383–434

    work page 1967

  50. [50]

    Howard, J. 2001. Mechanics of Motor Proteins and the Cytoskeleton . Sinauer Associates, Sunderland, MA

    work page 2001

  51. [51]

    Reconditi, M

    Piazzesi, G., M. Reconditi, M. Linari, L. Lucii, P. Bianco, E. Brunello, V . Decostre, A. Stewart, D. B. Gore, T. C. Irving, M. Irving, and V . Lombardi. 2007. Skeletal muscle performance determined by modulation of number of myosin motors rather than motor force or stroke size. Cell. 131:784-95

    work page 2007

  52. [52]

    Trivedi, D. V ., J. M. Muretta, A. M. Swenson, J. P. Davis, D. D. Thomas, and C. M. Yengo. 2015. Direct measurements of the coordination of lever arm swing and the catalytic cycle in myosin V . Proc. Natl. Acad. Sci. U. S. A. 112:14593-14598

    work page 2015

  53. [53]

    Barclay, C. J., R. C. Woledge, and N. A. Curtin. 2010. Inferring crossbridge properties from skeletal muscle energetics. Prog. Biophys. Mol. Biol. 102(1):53-71

    work page 2010

  54. [54]

    de Haan, A. 1998. The influence of stimulation frequency on force -velocity characteristics of in situ rat medial gastrocnemius muscle. Exp. Physiol. 83(1):77– 84

    work page 1998

  55. [55]

    Buschman, H. P. J., M. Linari, G. Elzinga, R. C. Woledge. 1997. Mechanical and energy characteristics during shortening in isolated type -1 muscle fibres from Xenopus laevis studied at maximal and submaximal activation. Pflügers. Arch. 435:145–150

    work page 1997

  56. [56]

    Hill, A. V . 1913. The absolute mechanical efficiency of the contraction of an isolated muscle. J. Physiol. 46 (1913):435-469

    work page 1913

  57. [57]

    Hill, A. V . 1964. The efficiency of mechanical power development during muscular shortening and its relation to load. Proc. R. Soc. Lond. B. 159 (1964):319-324

    work page 1964

  58. [58]

    Xie, X. S. 2013. Biochemistry. Enzyme kinetics, past and present. Science. 342:1457-1459

    work page 2013

  59. [59]

    Fryer, M. W., V . J. Owen, G. D. Lamb, and D. G. Stephenson. 1995. Effects of creatine phosphate and Pi on Ca2+ movements and tension development in rat skinned skeletal muscle fibres. J. Physiol. 482 ( Pt 1)(Pt 1):123–140

    work page 1995

  60. [60]

    Esposito, E

    Danieli-Betto, D., A. Esposito, E. Germinario, D. Sandonà, T. Martinello, A. Jakubiec-Puka, D. Biral, and R. Betto. 2005. Deficiency of alpha -sarcoglycan differently affects fast - and slow -twitch skeletal muscles. Am. J. Physiol. Regul. Integr. Comp. Physiol. 289(5):R1328–R1337. Supplemental Information a) Model validation We have performed a series of...

    work page 2005

  61. [61]

    + 3.32/𝑠 exp( 𝑓 1.2) 𝑓 ≤ 5𝑝𝑁 178.45/s + 166/𝑠 exp(− 𝑓

    work page

  62. [62]

    Table S2 Default values of parameters for rate formulas used in the simulations Parameter Value Parameter Value 𝐾 4.5 mV (24) 𝛼 0.2 /ms (24) 𝑉̅ -20 mV (24) 𝑘B𝑇 4.14×10−21 J

    𝑓 > 5𝑝𝑁 and 𝑘𝑏𝑟𝑒𝑎𝑘2 = 199/𝑠 exp(− 𝑓 1.5) + 8.3/𝑠 exp( 𝑓 5), respectively (18). Table S2 Default values of parameters for rate formulas used in the simulations Parameter Value Parameter Value 𝐾 4.5 mV (24) 𝛼 0.2 /ms (24) 𝑉̅ -20 mV (24) 𝑘B𝑇 4.14×10−21 J

    work page