Computational modelling identifies the impact of subtle anatomical variations between amphibian and mammalian skeletal muscle on spatiotemporal calcium dynamics.

The physical sites of calcium entry and exit in the skeletal muscle cell are distinct and highly organised in space. It was investigated whether the highly structured spatial organisation of sites of Ca(2+) release, uptake and action in skeletal muscle cells substantially impacts the dynamics of cytosolic Ca(2+) handling and thereby the physiology of the cell. Hereto, the spatiotemporal dynamics of the free calcium distribution in a fast-twitch (FT) muscle sarcomere was studied using a reaction-diffusion computational model for two genotypes with known anatomical differences. A computational model of a murine FT muscle sarcomere is developed, de novo including a closed calcium mass balance to simulate spatiotemporal high stimulation frequency calcium dynamics at 35 degrees C. Literature data on high-frequency calcium dye measurements were used as a first step towards model validation. The murine and amphibian sarcomere models were phenotypically distinct to capture known differences in positions of troponin C, actin-myosin overlap and calcium release within the sarcomere between frog and mouse. The models predicted large calcium gradients throughout the myoplasm as well as differences in calcium concentrations near the mitochondria of frog and mouse. Furthermore, the predicted Ca(2+) concentration was high at positions where Ca(2+) has a regulatory function, close to the mitochondria and troponin C.

Computational model of excitable cell indicates ATP free energy dynamics in response to calcium oscillations are undampened by cytosolic ATP buffers.

Mitochondria in excitable cells are recurrently exposed to pulsatile calcium gradients that activate cell function. Rapid calcium uptake by the mitochondria has previously been shown to cause uncoupling of oxidative phosphorylation. To test (i) if periodic nerve firing may cause oscillation of the cytosolic thermodynamic potential of ATP hydrolysis and (ii) if cytosolic adenylate (AK) and creatine kinase (CK) ATP buffering reactions dampen such oscillations, a lumped kinetic model of an excitable cell capturing major aspects of the physiology has been developed. Activation of ATP metabolism by low-frequency calcium pulses caused large oscillation of the cytosolic, but not mitochondrial ATP/ADP, ratio. This outcome was independent of net ATP synthesis or hydrolysis during mitochondrial calcium uptake. The AK/CK ATP buffering reactions dampened the amplitude and rate of cytosolic ATP/ADP changes on a timescale of seconds, but not milliseconds. These model predictions suggest that alternative sources of capacitance in neurons and striated muscles should be considered to protect ATP-free energy-driven cell functions.