From sarcomere to cell: An efficient algorithm for linking mathematical models of muscle contraction

Two classes of mathematical framework have previously been developed to model active tension generation in contracting muscle. Cross-bridge models of muscle are biophysically based but computationally expensive to solve, and thus unsuitable for embedding in spatially distributed continuum representations. Fading memory models are computationally efficient but provide limited biophysical insight. In this study a novel computational method is proposed for coupling these two frameworks such that biophysical events can be determined and computational tractability maintained. Within the cross-bridge model, the functional forms of the distribution of cross-bridges, as a function of strain in each state, are approximated using the distribution moment approach. Using the variables of area, mean and standard deviation of each distribution, analytic expressions are developed to calculate the temporal dynamics of stiffness, tension and energy. A root finding method is employed to adjust the variables such that the temporal dynamics of the cross-bridge model match those of an equivalent fading memory model. The method is demonstrated for sinusoidal perturbations in length at two frequencies, with an approximate 30-fold increase in computational efficiency over a conventional technique for finding a solution to the cross-bridge model.

[1]  Philip W. Kuchel,et al.  Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: computer simulation and Metabolic Control Analysis , 1999 .

[2]  H R Halvorson,et al.  Elementary steps of contraction probed by sinusoidal analysis technique in rabbit psoas fibers. , 1993, Advances in experimental medicine and biology.

[3]  Theo Arts,et al.  Cardiac Mechanoenergetics Replicated by Cross-Bridge Model , 2000, Annals of Biomedical Engineering.

[4]  P. Hunter,et al.  Integration from proteins to organs: the Physiome Project , 2003, Nature Reviews Molecular Cell Biology.

[5]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[6]  R. Winslow,et al.  Cardiac Ca2+ dynamics: the roles of ryanodine receptor adaptation and sarcoplasmic reticulum load. , 1998, Biophysical journal.

[7]  D. Noble,et al.  Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. , 1998, The Canadian journal of cardiology.

[8]  Nicolas P. Smith,et al.  Mathematical modelling of the heart: cell to organ , 2002 .

[9]  D. Smith,et al.  The theory of sliding filament models for muscle contraction. III. Dynamics of the five-state model. , 1990, Journal of theoretical biology.

[10]  G. Piazzesi,et al.  A cross-bridge model that is able to explain mechanical and energetic properties of shortening muscle. , 1995, Biophysical journal.

[11]  Masataka Kawai,et al.  Sinusoidal analysis: a high resolution method for correlating biochemical reactions with physiological processes in activated skeletal muscles of rabbit, frog and crayfish , 1980, Journal of Muscle Research & Cell Motility.

[12]  P W Kuchel,et al.  Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: in vivo kinetic characterization of 2,3-bisphosphoglycerate synthase/phosphatase using 13C and 31P NMR. , 1999, The Biochemical journal.

[13]  G. Zahalak A distribution-moment approximation for kinetic theories of muscular contraction , 1981 .

[14]  Martin Berzins,et al.  Developing software for time-dependent problems using the method of lines and differential-algebraic integrators , 1989 .

[15]  P. Hunter,et al.  Modelling the mechanical properties of cardiac muscle. , 1998, Progress in biophysics and molecular biology.

[16]  D. Smith,et al.  A strain-dependent ratchet model for [phosphate]- and [ATP]-dependent muscle contraction , 1998, Journal of Muscle Research & Cell Motility.

[17]  G. Zahalak The two-state cross-bridge model of muscle is an asymptotic limit of multi-state models. , 2000, Journal of theoretical biology.

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

[19]  T. L. Hill,et al.  Cross-bridge model of muscle contraction. Quantitative analysis. , 1980, Biophysical journal.

[20]  M Kawai,et al.  Crossbridge scheme and the kinetic constants of elementary steps deduced from chemically skinned papillary and trabecular muscles of the ferret. , 1993, Circulation research.

[21]  G. Zahalak,et al.  Muscle activation and contraction: constitutive relations based directly on cross-bridge kinetics. , 1990, Journal of biomechanical engineering.

[22]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. II. Afterdepolarizations, triggered activity, and potentiation. , 1994, Circulation research.

[23]  G I Zahalak,et al.  A distribution-moment model of deactivation in cardiac muscle. , 1998, Journal of biomechanics.

[24]  T. L. Hill,et al.  Theoretical formalism for the sliding filament model of contraction of striated muscle. Part I. , 1974, Progress in biophysics and molecular biology.