Parameter Estimation in a Crossbridge Muscle Model

[+] Author and Article Information
David C. Lin

Departments of Biological Systems Engineering and Veterinary and Comparative Anatomy, Pharmacology, and Physiology, Washington State University, Pullman, WA 99163

T. Richard Nichols

Department of Physiology, Emory School of Medicine, Atlanta, GA 30322

J Biomech Eng 125(1), 132-140 (Feb 14, 2003) (9 pages) doi:10.1115/1.1537262 History: Received December 01, 2001; Revised August 01, 2002; Online February 14, 2003
Copyright © 2003 by ASME
Your Session has timed out. Please sign back in to continue.


Sandercock, T. G., Lin, D. C., and Rymer, W. Z., “Muscle Models,” in Handbook of Brain Theory and Neural Networks, M. A. Arbib, Editor. In press.
Riener,  R., 1999, “Model-Based Development of Neuroprosthesis for Paraplegic Patients,” Philos. Trans. R. Soc. London, 354(1385), pp. 877–894.
Pette,  D., and Staron,  R. S., 2000, “Myosin Isoforms, Muscle Fiber Types, and Transitions,” Microsc. Res. Tech., 50(6), pp. 500–509.
Morgan,  D. L., and Allen,  D. G., 1999, “Early Events in Stretch-Induced Muscle Damage,” J. Appl. Physiol., 87(6), pp. 2007–2015.
Brown,  I. E., Scott,  S. H., and Loeb,  G. E., 1996, “Mechanics of Feline Soleus: II. Design and Validation of a Mathematical Model,” J. Muscle Res. Cell Motil., 17(2), pp. 221–233.
Zajac,  F. E., 1989, “Muscle and Tendon: Properties, Models, Scaling, and Application to Biomechanics and Motor Control,” Crit. Rev. Biomed. Eng., 17(4), pp. 359–411.
Hill,  A. V., 1938, “The Heat of Shortening and the Dynamic Constants of Muscle,” Proc. R. Soc. London, 126, pp. 136–195.
Iwamoto,  H., Sugaya,  R., and Sugi,  H., 1990, “Force-Velocity Relation of Frog Skeletal Muscle Fibres Shortening under Continuously Changing Load,” J. Physiol. (London), 422, pp. 185–202.
Talmadge,  R. J., Grossman,  E. J., and Roy,  R. R., 1996, “Myosin Heavy Chain Composition of Adult Feline (Felis Catus) Limb and Diaphragm Muscles,” J. Exp. Zool., 275(6), pp. 413–420.
Malamud,  J. G., Godt,  R. E., and Nichols,  T. R., 1996, “Relationship between Short-Range Stiffness and Yielding in Type-Identified, Chemically Skinned Muscle Fibers from the Cat Triceps Surae Muscles,” J. Neurophysiol., 76(4), pp. 2280–2289.
Godt,  R. E., and Nosek,  T. M., 1989, “Changes of Intracellular Milieu with Fatigue or Hypoxia Depress Contraction of Skinned Rabbit Skeletal and Cardiac Muscle,” J. Physiol. (London), 412, pp. 155–180.
Scott,  S. H., Brown,  I. E., and Loeb,  G. E., 1996, “Mechanics of Feline Soleus: I. Effect of Fascicle Length and Velocity on Force Output,” J. Muscle Res. Cell Motil., 17(2), pp. 207–219.
Mashima,  H., Akazawa,  K., Kushima,  H., and Fujii,  K., 1972, “The Force-Load-Velocity Relation and the Viscous-Like Force in the Frog Skeletal Muscle,” Jpn. J. Physiol., 22(1), pp. 103–120.
Huxley,  A. F., 1957, “Muscle Structure and Theories of Contraction,” Prog. Biophys. Biophys. Chem., 7, pp. 257–318.
Zahalak,  G. I., and Ma,  S. P., 1990, “Muscle Activation and Contraction: Constitutive Relations Based Directly on Cross-Bridge Kinetics,” ASME J. Biomech. Eng., 112(1), pp. 52–62.
Zahalak,  G. I., 1981, “A Distribution-Moment Approximation for Kinetic Theories of Muscular Contraction,” Math. Biosci., 55, pp. 89–114.
Ma,  S. P., and Zahalak,  G. I., 1991, “A Distribution-Moment Model of Energetics in Skeletal Muscle,” J. Biomech., 24(1), pp. 21–35.
Brown,  I. E., and Loeb,  G. E., 1999, “Measured and Modeled Properties of Mammalian Skeletal Muscle. I. The Effects of Post-Activation Potentiation on the Time Course and Velocity Dependencies of Force Production,” J. Muscle Res. Cell Motil., 20(5–6), pp. 443–456.
Joyce,  G. C., Rack,  P. M., and Westbury,  D. R., 1969, “The Mechanical Properties of Cat Soleus Muscle During Controlled Lengthening and Shortening Movements,” J. Physiol. (London), 204(2), pp. 461–474.
Manly, B. F. J., 1997, Randomization, Bootstrap, and Monte Carlo Methods in Biology, Chapman and Hall, London.
Morgan,  D. L., 1990, “New Insights into the Behavior of Muscle During Active Lengthening,” Biophys. J., 57(2), pp. 209–221.
Herzog,  W., and Leonard,  T. R., 2000, “The History Dependence of Force Production in Mammalian Skeletal Muscle Following Stretch-Shortening and Shortening-Stretch Cycles,” J. Biomech., 33(5), pp. 531–542.
Rassier,  D. E., MacIntosh,  B. R., and Herzog,  W., 1999, “Length Dependence of Active Force Production in Skeletal Muscle,” J. Appl. Physiol., 86(5), pp. 1445–1457.
Edman,  K. A., Caputo,  C., and Lou,  F., 1993, “Depression of Tetanic Force Induced by Loaded Shortening of Frog Muscle Fibres,” J. Physiol. (London), 466, pp. 535–552.
Edman,  K. A., and Tsuchiya,  T., 1996, “Strain of Passive Elements During Force Enhancement by Stretch in Frog Muscle Fibres,” J. Physiol. (London), 490(Pt 1), pp. 191–205.
Woledge,  R. C., Curtin,  N. A., and Homsher,  E., 1985, “Energetic Aspects of Muscle Contraction,” Monogr. Physiol. Soc., 41, pp. 1–357.
Spector,  S. A., Gardiner,  P. F., Zernicke,  R. F., Roy,  R. R., and Edgerton,  V. R., 1980, “Muscle Architecture and Force-Velocity Characteristics of Cat Soleus and Medial Gastrocnemius: Implications for Motor Control,” J. Neurophysiol., 44(5), pp. 951–960.
Sandercock,  T. G., and Heckman,  C. J., 1997, “Force from Cat Soleus Muscle During Imposed Locomotor-Like Movements: Experimental Data Versus Hill-Type Model Predictions,” J. Neurophysiol., 77(3), pp. 1538–1552.
Widrick,  J. J., Romatowski,  J. G., Karhanek,  M., and Fitts,  R. H., 1997, “Contractile Properties of Rat, Rhesus Monkey, and Human Type I Muscle Fibers,” Am. J. Physiol., 272(1 Pt 2), pp. R34–42.
Galler,  S., Schmitt,  T. L., and Pette,  D., 1994, “Stretch Activation, Unloaded Shortening Velocity and Myosin Heavy Chain Isoforms of Rat Skeletal Muscle Fibres,” J. Physiol. (London), 478(Pt 3), pp. 513–521.
Fitts,  R. H., and Widrick,  J. J., 1996, “Muscle Mechanics: Adaptations with Exercise-Training,” Exerc Sport Sci. Rev., 24, pp. 427–473.
Stienen,  G. J., Versteeg,  P. G., Papp,  Z., and Elzinga,  G., 1992, “Mechanical Properties of Skinned Rabbit Psoas and Soleus Muscle Fibres During Lengthening: Effects of Phosphate and Ca2+,” J. Physiol. (London), 451, pp. 503–523.
Cole,  G. K., van den Bogert,  A. J., Herzog,  W., and Gerritsen,  K. G., 1996, “Modelling of Force Production in Skeletal Muscle Undergoing Stretch,” J. Biomech., 29(8), pp. 1091–1104.
Potma,  E. J., van Graas,  I. A., and Stienen,  G. J., 1994, “Effects of pH on Myofibrillar Atpase Activity in Fast and Slow Skeletal Muscle Fibers of the Rabbit,” Biophys. J., 67(6), pp. 2404–2410.
Reggiani,  C., Potma,  E. J., Bottinelli,  R., Canepari,  M., Pellegrino,  M. A., and Stienen,  G. J., 1997, “Chemo-Mechanical Energy Transduction in Relation to Myosin Isoform Composition in Skeletal Muscle Fibres of the Rat,” J. Physiol. (London), 502(Pt 2), pp. 449–460.
Zahalak,  G. I., 1986, “A Comparison of the Mechanical Behavior of the Cat Soleus Muscle with a Distribution-Moment Model,” J. Biomech. Eng., 108(2), pp. 131–140.
Lin,  D. C., and Rymer,  W. Z., 1993, “Mechanical Properties of Cat Soleus Muscle Elicited by Sequential Ramp Stretches: Implications for Control of Muscle,” J. Clin. Neurophysiol., 70(3), pp. 997–1008.
Harry,  J. D., Ward,  A. W., Heglund,  N. C., Morgan,  D. L., and McMahon,  T. A., 1990, “Cross-Bridge Cycling Theories Cannot Explain High-Speed Lengthening Behavior in Frog Muscle,” Biophys. J., 57(2), pp. 201–208.
Brenner,  B., 1986, “The Cross-Bridge Cycle in Muscle. Mechanical, Biochemical, and Structural Studies on Single Skinned Rabbit Psoas Fibers to Characterize Cross-Bridge Kinetics in Muscle for Correlation with the Actomyosin-Atpase in Solution,” Basic Res. Cardiol., 81(Supplement 1), pp. 1–15.
Getz,  E. B., Cooke,  R., and Lehman,  S. L., 1998, “Phase Transition in Force During Ramp Stretches of Skeletal Muscle,” Biophys. J., 75(6), pp. 2971–2983.
Huxley,  A. F., and Simmons,  R. M., 1971, “Proposed Mechanism of Force Generation in Striated Muscle,” Nature (London), 233(5321), pp. 533–538.
Campbell,  K. S., and Moss,  R. L., 2000, “A Thixotropic Effect in Contracting Rabbit Psoas Muscle: Prior Movement Reduces the Initial Tension Response to Stretch,” J. Physiol. (London), 525(Pt 2), pp. 531–548.


Grahic Jump Location
The distribution-moment crossbridge model parameters. The solid line describes the attachment rate (f(x) in Eq. (4)) as a function of crossbridge position, and the dashed lines describe the detachment rate (g(x)). The constant “h” is often called the crossbridge reach.
Grahic Jump Location
The distribution of parameter estimates calculated from the bootstrapping 100 different samples. Shown are the distributions which appear the most normal and the least normal.
Grahic Jump Location
The mean normalized standard deviation of each parameter estimate, calculated by dividing the standard deviation by the mean of the parameter estimates derived from the bootstrapping method.
Grahic Jump Location
Block diagram of the experimental servo system which implemented force controlled perturbations. The dSpace board was a digital signal processing (DSP) board with analog to digital and digital to analog capabilities which implemented the blocks within the dotted outline. Fd(t) was the desired time-course of force which was generated by software. G1 and G2 are the gains of the force controller (see Eq. (1)), Ld(t) was the input to the length-servo, L(t) was the actual fiber length, and F(t) was the force generated by the muscle fiber.
Grahic Jump Location
The staircase length perturbation which was used to record the F-L characteristics of a single fiber. Shown are the active and passive forces recorded from one fiber during the shortening (heavy line) and the lengthening (light line) perturbation. The active force was estimated by subtracting the force recorded in the passive trial from the force recorded in the active trial. Fiber lengths were normalized to the initial fiber length (Lo).
Grahic Jump Location
The force ramp perturbation which was used to record the F-V characteristics of a single fiber. Shown are two trials from one fiber, one shortening (heavy line) and one lengthening (light line). The rates of the force ramp were ±1.5 Po/s. The active force was estimated by subtracting the force recorded in the passive trial from the force recorded in the active trial.
Grahic Jump Location
Averaged F-L data across all fiber data, plotted with standard deviation. The force data of each fiber were first normalized to the initial active isometric force (Fo) and then averaged across fibers.
Grahic Jump Location
The F-V data of an individual fiber. Shown are the filtered experimental data (solid line) and the best fit of the hyperbolic Eqs. (2) and (3) (dashed-dot line).
Grahic Jump Location
Averaged F-V data across all fiber data, plotted with standard deviation (dashed line). The F-V data of each fiber were first fitted to hyperbolic equations, the velocity evaluated at specific forces, and the velocities of all trials averaged at those forces.
Grahic Jump Location
The averaged F-V data (solid line) of Fig. 7 overlaid with the DM model simulation (dashed line) with the best fitting parameters. The simplex algorithm minimized the sum of the squared velocity error difference between the simulation and experimental data.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In