Approximation for Cooperative Interactions of a Spatially-Detailed Cardiac Sarcomere Model

We developed a novel ordinary differential equation (ODE) model, which produced results that correlated well with the Monte Carlo (MC) simulation when applied to a spatially-detailed model of the cardiac sarcomere. Configuration of the novel ODE model was based on the Ising model of myofilaments, with the "co-operative activation" effect introduced to incorporate nearest-neighbor interactions. First, a set of parameters was estimated using arbitrary Ca transient data to reproduce the combinational probability for the states of three consecutive regulatory units, using single unit probabilities for central and neighboring units in the MC simulation. The parameter set thus obtained enabled the calculation of the state transition of each unit using the ODE model with reference to the neighboring states. The present ODE model not only provided good agreement with the MC simulation results but was also capable of reproducing a wide range of experimental results under both steady-state and dynamic conditions including shortening twitch. The simulation results suggested that the nearest-neighbor interaction is a reasonable approximation of the cooperativity based on end-to-end interactions. Utilizing the modified ODE model resulted in a reduction in computational costs but maintained spatial integrity and co-operative effects, making it a powerful tool in cardiac modeling.