Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations

Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations

Model Status

This cellML model is known to run in both PCEnv and COR. The units have been checked and they are consistent. This second version of the CellML model has been modified by Steve Niederer. The CellML model was compared with Jeremy Rice's original Matlab code taken from his web page. Changes to the CellML model include: body temperature in the Rice Matlab code is acurate to 3sf ie body temp is 310 not 310.15. Also in the Rice Matlab code derivative 12 the force integral is fixed to be zero. To fix the CellML model Steve introduced:

SEon = 1.0; To cellml model and change afterload definition to: if (SEon == 1) afterload = KSE*(SLset-Y(5)); else afterload = 0.0; end;

and changed xbmodsp = 4/3 (which is the mouse specific modification for XB cycling rates) to xbmodsp = 1 (to represent the rat model and not the mouse)

And finally Steve also reset all of the initial conditions after solving the model a coupl of times for 1000ms.

Model Structure

Abstract: We develop a point model of the cardiac myofilament (MF) to simulate wide variety of experimental muscle characterizations including Force-Ca relations and twitches under isometric, isosarcometric, isotonic and auxotonic conditions. Complex MF behaviors are difficult to model because spatial interactions cannot be directly implemented as ordinary differential equations (ODEs). We therefore allow phenomenological approximations with careful consideration to the relationships with the underlying biophysical mechanisms. We describe new formulations that avoid mean-field approximations found in most existing MF models. To increase the scope and applicability of the model, we include length- and temperature-dependent effects that play important roles in MF response. We have also included a representation of passive restoring forces to simulate isolated cell shortening protocols. Possessing both computational efficiency and the ability to simulate a wide variety of muscle responses, the MF representation is well-suited for coupling to existing cardiac cell models of electrophysiology and Ca-handling mechanisms. To illustrate this suitability, the MF model is coupled to the Chicago rabbit cardiomyocyte model. The combined model generates realistic appearing action potentials, intracellular Ca transients, and cell shortening signals. The combined model also demonstrates that the feedback effects of force on Ca binding to troponin can modify the cystolic Ca transient.

Note that the current CellMl model describes the rat model of myofilamentcontraction. One of the goals of the published work by Rice et al. was to develop a model of the myofilaments which could be coupled to existing models of cardiac electrophysiology and Ca-handling mechanisms. Similarly, it is our intention to be able to couple the CellML model of myofilament contration with cardiac models of electrophysiology and calcium change. To achieve this we will first have to modify the rat myofilament model to convert it to a rabbit myofilament model. This can be achieved by decreasing the transition rates in the crossbridge cycle by a factor of 5 to simulate the changes in myosin isoforms (rat is predominately V1 while rabbit is V3) and species-based changes in rate. We can then couple the model to the Chicago rabbit ventricular myocyte model by using the cytosolic Ca concentration ([Ca]c) from the Chicago model as the input to the myofilament model. Other minor modifications are also required, as described in the paper below:

Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations, John Jeremy Rice, Fei Wang, Donald M. Bers and Pieter P. de Tombe, 2008, Biophysical Journal . (A PDF version of the article is available to journal subscribers on the Biophysical Journal website.) PubMed ID: 18234826

Schematic diagram of the model.