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

Computational Model of the Effect of Mitochondrial Dysfunction on Excitation–Contraction Coupling in Skeletal Muscle

It has become well established that mitochondria not only regulate myoplasmic calcium in skeletal muscle, but also use that calcium to stimulate oxidative phosphorylation (OXPHOS). While experimental approaches have allowed for imaging of mitochondrial calcium and membrane potentials in isolated fibers, capturing the role of mitochondria and the impact of mitochondrial impairments on excitation–contraction coupling (ECC) remains difficult to explore in intact muscle. Computational models have been widely used to examine the structure and function of skeletal muscle contraction; however, models of ECC to date lack communication between the myoplasm and mitochondria for regulating calcium and ATP during sustained contractions. To address this, a mathematical model of mitochondrial calcium handling and OXPHOS was integrated into a physiological model of ECC incorporating action potential propagation, calcium handling between the sarcoplasmic reticulum (SR) and the myoplasm, and crossbridge cycling. The model was used to examine the protective role of mitochondria during repeated stimulation and the impact of mitochondrial dysfunction on ECC resulting from progressive OXPHOS inhibition. Pathological myoplasmic calcium accumulation occurred through distinct mechanisms in the model in the case of either electron transport chain, F1F0 ATP synthase, or adenine nucleotide transporter impairments. To investigate the effect of each impairment on force, a model of calcium-stimulated apoptosis was utilized to capture dysfunction-induced reductions in muscle mass, driving whole muscle force loss. The model presented in this study can be used to examine the role of mitochondria in the regulation of calcium, ATP, and force generation during voluntary contraction.

Calcium Dynamics of Ex Vivo Long-Term Cultured CD8+ T Cells Are Regulated by Changes in Redox Metabolism

T cells reach a state of replicative senescence characterized by a decreased ability to proliferate and respond to foreign antigens. Calcium release associated with TCR engagement is widely used as a surrogate measure of T cell response. Using an ex vivo culture model that partially replicates features of organismal aging, we observe that while the amplitude of Ca2+ signaling does not change with time in culture, older T cells exhibit faster Ca2+ rise and a faster decay. Gene expression analysis of Ca2+ channels and pumps expressed in T cells by RT-qPCR identified overexpression of the plasma membrane CRAC channel subunit ORAI1 and PMCA in older T cells. To test whether overexpression of the plasma membrane Ca2+ channel is sufficient to explain the kinetic information, we adapted a previously published computational model by Maurya and Subramaniam to include additional details on the store-operated calcium entry (SOCE) process to recapitulate Ca2+ dynamics after T cell receptor stimulation. Simulations demonstrated that upregulation of ORAI1 and PMCA channels is not sufficient to explain the observed alterations in Ca2+ signaling. Instead, modeling analysis identified kinetic parameters associated with the IP3R and STIM1 channels as potential causes for alterations in Ca2+ dynamics associated with the long term ex vivo culturing protocol. Due to these proteins having known cysteine residues susceptible to oxidation, we subsequently investigated and observed transcriptional remodeling of metabolic enzymes, a shift to more oxidized redox couples, and post-translational thiol oxidation of STIM1. The model-directed findings from this study highlight changes in the cellular redox environment that may ultimately lead to altered T cell calcium dynamics during immunosenescence or organismal aging.

The sensitivity of fast muscle contractile function to the major components of the sarcomere Ca(2+)-cycling system

A reaction-diffusion model of a muscle sarcomere was developed to evaluate the sensitivity of force characteristics to diffusion and Ca(2+)-cycling components. The model compared well to experimental force measurements. Diffusion led to Ca(2+) gradients that enhanced maximal force and accelerated relaxation compared to when diffusion was infinitely fast. However, a modest increase in sarcomere length or radius led to a decrease in maximal force. Lowering the Ca(2+) release rate caused a lower maximal force, but increasing the rate led to only modest gains in maximal force while incurring much greater ATP costs associated with reuptake. Greater parvalbumin binding rates decreased maximal force but enhanced relaxation, and this effect was magnified when Ca(2+) uptake rates were lowered as may occur during fatigue. These results show a physiological set of parameters that lead to a functional sarcomere of known dimensions and contractile function, and the effects of parameter variation on muscle function.

Muscle-type specific autophosphorylation of CaMKII isoforms after paced contractions

We explored to what extent isoforms of the regulator of excitation-contraction and excitation-transcription coupling, calcium/calmodulin protein kinase II (CaMKII) contribute to the specificity of myocellular calcium sensing between muscle types and whether concentration transients in its autophosphorylation can be simulated. CaMKII autophosphorylation at Thr287 was assessed in three muscle compartments of the rat after slow or fast motor unit-type stimulation and was compared against a computational model (CaMuZclE) coupling myocellular calcium dynamics with CaMKII Thr287 phosphorylation. Qualitative differences existed between fast- (gastrocnemius medialis) and slow-type muscle (soleus) for the expression pattern of CaMKII isoforms. Phospho-Thr287 content of δA CaMKII, associated with nuclear functions, demonstrated a transient and compartment-specific increase after excitation, which contrasted to the delayed autophosphorylation of the sarcoplasmic reticulum-associated βM CaMKII. In soleus muscle, excitation-induced δA CaMKII autophosphorylation demonstrated frequency dependence (P = 0.02). In the glycolytic compartment of gastrocnemius medialis, CaMKII autophosphorylation after excitation was blunted. In silico assessment emphasized the importance of mitochondrial calcium buffer capacity for excitation-induced CaMKII autophosphorylation but did not predict its isoform specificity. The findings expose that CaMKII autophosphorylation with paced contractions is regulated in an isoform and muscle type-specific fashion and highlight properties emerging for phenotype-specific regulation of CaMKII.

Combined in vivo and in silico investigations of activation of glycolysis in contracting skeletal muscle

The hypothesis was tested that the variation of in vivo glycolytic flux with contraction frequency in skeletal muscle can be qualitatively and quantitatively explained by calcium-calmodulin activation of phosphofructokinase (PFK-1). Ischemic rat tibialis anterior muscle was electrically stimulated at frequencies between 0 and 80 Hz to covary the ATP turnover rate and calcium concentration in the tissue. Estimates of in vivo glycolytic rates and cellular free energetic states were derived from dynamic changes in intramuscular pH and phosphocreatine content, respectively, determined by phosphorus magnetic resonance spectroscopy ((31)P-MRS). Computational modeling was applied to relate these empirical observations to understanding of the biochemistry of muscle glycolysis. Hereto, the kinetic model of PFK activity in a previously reported mathematical model of the glycolytic pathway (Vinnakota KC, Rusk J, Palmer L, Shankland E, Kushmerick MJ. J Physiol 588: 1961-1983, 2010) was adapted to contain a calcium-calmodulin binding sensitivity. The two main results were introduction of regulation of PFK-1 activity by binding of a calcium-calmodulin complex in combination with activation by increased concentrations of AMP and ADP was essential to qualitatively and quantitatively explain the experimental observations. Secondly, the model predicted that shutdown of glycolytic ATP production flux in muscle postexercise may lag behind deactivation of PFK-1 (timescales: 5-10 s vs. 100-200 ms, respectively) as a result of accumulation of glycolytic intermediates downstream of PFK during contractions.

Molecules in motion: influences of diffusion on metabolic structure and function in skeletal muscle

Metabolic processes are often represented as a group of metabolites that interact through enzymatic reactions, thus forming a network of linked biochemical pathways. Implicit in this view is that diffusion of metabolites to and from enzymes is very fast compared with reaction rates, and metabolic fluxes are therefore almost exclusively dictated by catalytic properties. However, diffusion may exert greater control over the rates of reactions through: (1) an increase in reaction rates; (2) an increase in diffusion distances; or (3) a decrease in the relevant diffusion coefficients. It is therefore not surprising that skeletal muscle fibers have long been the focus of reaction-diffusion analyses because they have high and variable rates of ATP turnover, long diffusion distances, and hindered metabolite diffusion due to an abundance of intracellular barriers. Examination of the diversity of skeletal muscle fiber designs found in animals provides insights into the role that diffusion plays in governing both rates of metabolic fluxes and cellular organization. Experimental measurements of metabolic fluxes, diffusion distances and diffusion coefficients, coupled with reaction-diffusion mathematical models in a range of muscle types has started to reveal some general principles guiding muscle structure and metabolic function. Foremost among these is that metabolic processes in muscles do, in fact, appear to be largely reaction controlled and are not greatly limited by diffusion. However, the influence of diffusion is apparent in patterns of fiber growth and metabolic organization that appear to result from selective pressure to maintain reaction control of metabolism in muscle.

Magnitude and control of mitochondrial sensitivity to ADP

The transduction function for ADP stimulation of mitochondrial ATP synthesis in skeletal muscle was reconstructed in vivo and in silico to investigate the magnitude and origin of mitochondrial sensitivity to cytoplasmic ADP concentration changes. Dynamic in vivo measurements of human leg muscle phosphocreatine (PCr) content during metabolic recovery from contractions were performed by (31)P-NMR spectroscopy. The cytoplasmic ADP concentration ([ADP]) and rate of oxidative ATP synthesis (Jp) at each time point were calculated from creatine kinase equilibrium and the derivative of a monoexponential fit to the PCr recovery data, respectively. Reconstructed [ADP]-Jp relations for individual muscles containing more than 100 data points were kinetically characterized by nonlinear curve fitting yielding an apparent kinetic order and ADP affinity of 1.9 +/- 0.2 and 0.022 +/- 0.003 mM, respectively (means +/- SD; n = 6). Next, in silico [ADP]-Jp relations for skeletal muscle were generated using a computational model of muscle oxidative ATP metabolism whereby model parameters corresponding to mitochondrial enzymes were randomly changed by 50-150% to determine control of mitochondrial ADP sensitivity. The multiparametric sensitivity analysis showed that mitochondrial ADP ultrasensitivity is an emergent property of the integrated mitochondrial enzyme network controlled primarily by kinetic properties of the adenine nucleotide translocator.