Interactions between Connected Half-Sarcomeres Produce Emergent Mechanical Behavior in a Mathematical Model of Muscle

Most reductionist theories of muscle attribute a fiber's mechanical properties to the scaled behavior of a single half-sarcomere. Mathematical models of this type can explain many of the known mechanical properties of muscle but have to incorporate a passive mechanical component that becomes ∼300% stiffer in activating conditions to reproduce the force response elicited by stretching a fast mammalian muscle fiber. The available experimental data suggests that titin filaments, which are the mostly likely source of the passive component, become at most ∼30% stiffer in saturating Ca2+ solutions. The work described in this manuscript used computer modeling to test an alternative systems theory that attributes the stretch response of a mammalian fiber to the composite behavior of a collection of half-sarcomeres. The principal finding was that the stretch response of a chemically permeabilized rabbit psoas fiber could be reproduced with a framework consisting of 300 half-sarcomeres arranged in 6 parallel myofibrils without requiring titin filaments to stiffen in activating solutions. Ablation of inter-myofibrillar links in the computer simulations lowered isometric force values and lowered energy absorption during a stretch. This computed behavior mimics effects previously observed in experiments using muscles from desmin-deficient mice in which the connections between Z-disks in adjacent myofibrils are presumably compromised. The current simulations suggest that muscle fibers exhibit emergent properties that reflect interactions between half-sarcomeres and are not properties of a single half-sarcomere in isolation. It is therefore likely that full quantitative understanding of a fiber's mechanical properties requires detailed analysis of a complete fiber system and cannot be achieved by focusing solely on the properties of a single half-sarcomere.

[1]  W. Steffen,et al.  Single-molecule measurement of the stiffness of the rigor myosin head. , 2008, Biophysical journal.

[2]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[3]  D L Morgan,et al.  The effect on tension of non‐uniform distribution of length changes applied to frog muscle fibres. , 1979, The Journal of physiology.

[4]  K. Edman,et al.  Redistribution of sarcomere length during isometric contraction of frog muscle fibres and its relation to tension creep. , 1984, The Journal of physiology.

[5]  K. Campbell,et al.  History-dependent mechanical properties of permeabilized rat soleus muscle fibers. , 2002, Biophysical journal.

[6]  Walter Herzog,et al.  Dynamics of individual sarcomeres during and after stretch in activated single myofibrils , 2003, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[7]  Kenneth S Campbell,et al.  Filament compliance effects can explain tension overshoots during force development. , 2006, Biophysical journal.

[8]  M. Lakie,et al.  A cross‐bridge mechanism can explain the thixotropic short‐range elastic component of relaxed frog skeletal muscle , 1998, The Journal of physiology.

[9]  K. Campbell,et al.  Myocardial short-range force responses increase with age in F344 rats. , 2009, Journal of molecular and cellular cardiology.

[10]  K. Edman,et al.  Absence of plateau of the sarcomere length-tension relation in frog muscle fibres. , 1984, Acta physiologica Scandinavica.

[11]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[12]  A. Huxley Muscular contraction. Review lecture , 1974 .

[13]  Marco Linari,et al.  Stiffness and fraction of Myosin motors responsible for active force in permeabilized muscle fibers from rabbit psoas. , 2007, Biophysical journal.

[14]  Thomas L. Daniel,et al.  Sarcomere Lattice Geometry Influences Cooperative Myosin Binding in Muscle , 2007, PLoS Comput. Biol..

[15]  Michael P. Sheetz,et al.  Stretching Single Talin Rod Molecules Activates Vinculin Binding , 2009, Science.

[16]  Andrew P. Weir,et al.  Function and genetics of dystrophin and dystrophin-related proteins in muscle. , 2002, Physiological reviews.

[17]  D. Rassier,et al.  The mechanical behavior of individual sarcomeres of myofibrils isolated from rabbit psoas muscle. , 2009, American journal of physiology. Cell physiology.

[18]  T L Daniel,et al.  Compliant realignment of binding sites in muscle: transient behavior and mechanical tuning. , 1998, Biophysical journal.

[19]  D. Morgan New insights into the behavior of muscle during active lengthening. , 1990, Biophysical journal.

[20]  S. Lehman,et al.  Phase transition in force during ramp stretches of skeletal muscle. , 1998, Biophysical journal.

[21]  J J Fredberg,et al.  On the theory of muscle contraction: filament extensibility and the development of isometric force and stiffness. , 1996, Biophysical journal.

[22]  L. McLoon,et al.  Complex three-dimensional patterns of myosin isoform expression: differences between and within specific extraocular muscles , 1999, Journal of Muscle Research & Cell Motility.

[23]  K. Campbell,et al.  A thixotropic effect in contracting rabbit psoas muscle: prior movement reduces the initial tension response to stretch , 2000, The Journal of physiology.

[24]  Y Ueno,et al.  X-ray diffraction evidence for the extensibility of actin and myosin filaments during muscle contraction. , 1994, Biophysical journal.

[25]  John Trinick,et al.  Titin: properties and family relationships , 2003, Nature Reviews Molecular Cell Biology.

[26]  Mauro Birattari,et al.  Swarm Intelligence , 2012, Lecture Notes in Computer Science.

[27]  H E Huxley,et al.  X-ray diffraction measurements of the extensibility of actin and myosin filaments in contracting muscle. , 1994, Biophysical journal.

[28]  S. Ishiwata,et al.  Inter-sarcomere coordination in muscle revealed through individual sarcomere response to quick stretch , 2009, Proceedings of the National Academy of Sciences.

[29]  K. Edman,et al.  Strain of passive elements during force enhancement by stretch in frog muscle fibres. , 1996, The Journal of physiology.

[30]  H. Higuchi,et al.  Mechanism of action of 2, 3-butanedione 2-monoxime on contraction of frog skeletal muscle fibres , 1988, Journal of Muscle Research & Cell Motility.

[31]  R L Lieber,et al.  Desmin knockout muscles generate lower stress and are less vulnerable to injury compared with wild-type muscles. , 2000, American journal of physiology. Cell physiology.

[32]  M. Noble,et al.  Residual force enhancement after stretch of contracting frog single muscle fibers , 1982, The Journal of general physiology.

[33]  A. Huxley Muscle structure and theories of contraction. , 1957, Progress in biophysics and biophysical chemistry.

[34]  D L Morgan,et al.  Intersarcomere dynamics during fixed‐end tetanic contractions of frog muscle fibres. , 1979, The Journal of physiology.

[35]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[36]  A. Arner,et al.  Desmin filaments influence myofilament spacing and lateral compliance of slow skeletal muscle fibers. , 2005, Biophysical journal.

[37]  A. Huxley,et al.  The maximum length for contraction in vertebrate striated muscle , 1961, The Journal of physiology.

[38]  Jachen Denoth,et al.  A multisegmental cross-bridge kinetics model of the myofibril. , 2009, Journal of theoretical biology.

[39]  G. Piazzesi,et al.  The contractile response during steady lengthening of stimulated frog muscle fibres. , 1990, The Journal of physiology.

[40]  Jachen Denoth,et al.  Half-sarcomere dynamics in myofibrils during activation and relaxation studied by tracking fluorescent markers. , 2006, Biophysical journal.

[41]  A. Huxley,et al.  Rapid regeneration of power stroke in contracting muscle by attachment of second myosin head , 1997, Journal of Muscle Research & Cell Motility.

[42]  S Mochon,et al.  A quantitative model of intersarcomere dynamics during fixed-end contractions of single frog muscle fibers. , 1982, Biophysical journal.

[43]  S. Tideswell,et al.  Filament compliance and tension transients in muscle , 1996, Journal of Muscle Research & Cell Motility.

[44]  Kenneth S Campbell,et al.  SLControl: PC-based data acquisition and analysis for muscle mechanics. , 2003, American journal of physiology. Heart and circulatory physiology.

[45]  A. Straight,et al.  Specificity of blebbistatin, an inhibitor of myosin II , 2004, Journal of Muscle Research & Cell Motility.

[46]  Yiming Wu,et al.  Calcium-dependent molecular spring elements in the giant protein titin , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[47]  R. Moss,et al.  Characteristics of troponin C binding to the myofibrillar thin filament: extraction of troponin C is not random along the length of the thin filament. , 1997, Biophysical journal.

[48]  Roger Cooke,et al.  A model of crossbridge action: the effects of ATP, ADP and Pi , 1989, Journal of Muscle Research & Cell Motility.

[49]  D. Noble Music of life : biology beyond the genome , 2006 .

[50]  A. Huxley,et al.  A system for fast recording of longitudinal displacement of a striated muscle fibre , 1981 .