The ForceLAB simulator: Application to the comparison of current models of cardiomyocyte contraction

Mathematical models are useful tools in the study of physiological phenomena. However, due to differences in assumptions and formulations, discrepancy in simulations may occur. Among the models for cardiomyocyte contraction based on Huxley's cross-bridge cycling, those proposed by Negroni and Lascano (NL) and Rice et al. (RWH) are the most frequently used. This study was aimed at developing a computational tool, ForceLAB, which allows implementing different contraction models and modifying several functional parameters. As an application, electrically-stimulated twitches triggered by an equal Ca²⁺ input and steady-state force x pCa relationship (pCa = -log of the molar free Ca²⁺ concentration) simulated with the NL and RWH models were compared. The equilibrium Ca²⁺-troponin C (TnC) dissociation constant (K_d) was modified by changing either the association (k_on) or the dissociation (k_off) rate constant. With the NL model, raising K_d by either maneuver decreased monotonically twitch amplitude and duration, as expected. With the RWH model, in contrast, the same K_d variation caused increase or decrease of peak force depending on which rate constant was modified. Additionally, force x pCa curves simulated using Ca²⁺ binding constants estimated in cardiomyocytes bearing wild-type and mutated TnC were compared to curves previously determined in permeabilized fibers. Mutations increased k_on and k_off, and decreased K_d. Both models produced curves fairly comparable to the experimental ones, although sensitivity to Ca²⁺ was greater, especially with RWH model. The NL model reproduced slightly better the qualitative changes associated with the mutations. It is expected that this tool can be useful for teaching and investigation.