Modelling of calcium handling in airway myocytes

The relative contributions of store-operated and voltage-gated Ca2+ channels to the control of Ca2+ oscillations in airway smooth muscle

KEY POINTS

Agonist-dependent oscillations in the concentration of free cytosolic calcium are a vital mechanism for the control of airway smooth muscle contraction and thus are a critical factor in airway hyper-responsiveness. Using a mathematical model, closely tied to experimental work, we show that the oscillations in membrane potential accompanying the calcium oscillations have no significant effect on the properties of the calcium oscillations. In addition, the model shows that calcium entry through store-operated calcium channels is critical for calcium oscillations, but calcium entry through voltage-gated channels has much less effect. The model predicts that voltage-gated channels are less important than store-operated channels in the control of airway smooth muscle tone.

ABSTRACT

Airway smooth muscle contraction is typically the key mechanism underlying airway hyper-responsiveness, and the strength of muscle contraction is determined by the frequency of oscillations of intracellular calcium (Ca²⁺ ) concentration. In airway smooth muscle cells, these Ca²⁺ oscillations are caused by cyclic Ca²⁺ release from the sarcoplasmic reticulum, although Ca²⁺ influx via plasma membrane channels is also necessary to sustain the oscillations over longer times. To assess the relative contributions of store-operated and voltage-gated Ca²⁺ channels to this Ca²⁺ influx, we generated a comprehensive mathematical model, based on experimental Ca²⁺ measurements in mouse precision-cut lung slices, to simulate Ca²⁺ oscillations and changes in membrane potential. Agonist-induced Ca²⁺ oscillations are accompanied by oscillations in membrane potential, although the membrane potential oscillations are too small to generate large Ca²⁺ currents through voltage-gated Ca²⁺ channels, and thus have little effect on the Ca²⁺ oscillations. Ca²⁺ entry through voltage-gated channels only becomes important when the cell is depolarized (e.g. by a high external K⁺ concentration). As a result, agonist-induced Ca²⁺ oscillations are critically dependent on Ca²⁺ entry through store-operated channels but do not depend strongly on Ca²⁺ entry though voltage-gated channels.

Activation of store-operated calcium entry in airway smooth muscle cells: insight from a mathematical model

Intracellular dynamics of airway smooth muscle cells (ASMC) mediate ASMC contraction and proliferation, and thus play a key role in airway hyper-responsiveness (AHR) and remodelling in asthma. We evaluate the importance of store-operated entry (SOCE) in these dynamics by constructing a mathematical model of ASMC signaling based on experimental data from lung slices. The model confirms that SOCE is elicited upon sufficient depletion of the sarcoplasmic reticulum (SR), while receptor-operated entry (ROCE) is inhibited in such conditions. It also shows that SOCE can sustain agonist-induced oscillations in the absence of other influx. SOCE up-regulation may thus contribute to AHR by increasing the oscillation frequency that in turn regulates ASMC contraction. The model also provides an explanation for the failure of the SERCA pump blocker CPA to clamp the cytosolic of ASMC in lung slices, by showing that CPA is unable to maintain the SR empty of . This prediction is confirmed by experimental data from mouse lung slices, and strongly suggests that CPA only partially inhibits SERCA in ASMC.

Inflammation alters regional mitochondrial Ca²+ in human airway smooth muscle cells

Regulation of cytosolic Ca(2+) concentration ([Ca(2+)](cyt)) in airway smooth muscle (ASM) is a key aspect of airway contractility and can be modulated by inflammation. Mitochondria have tremendous potential for buffering [Ca(2+)](cyt), helping prevent Ca(2+) overload, and modulating other intracellular events. Here, compartmentalization of mitochondria to different cellular regions may subserve different roles. In the present study, we examined the role of Ca(2+) buffering by mitochondria and mitochondrial Ca(2+) transport mechanisms in the regulation of [Ca(2+)](cyt) in enzymatically dissociated human ASM cells upon exposure to the proinflammatory cytokines TNF-α and IL-13. Cells were loaded simultaneously with fluo-3 AM and rhod-2 AM, and [Ca(2+)](cyt) and mitochondrial Ca(2+) concentration ([Ca(2+)](mito)) were measured, respectively, using real-time two-color fluorescence microscopy in both the perinuclear and distal, perimembranous regions of cells. Histamine induced a rapid increase in both [Ca(2+)](cyt) and [Ca(2+)](mito), with a significant delay in the mitochondrial response. Inhibition of the mitochondrial Na(+)/Ca(2+) exchanger (1 μM CGP-37157) increased [Ca(2+)](mito) responses in perinuclear mitochondria but not distal mitochondria. Inhibition of the mitochondrial uniporter (1 μM Ru360) decreased [Ca(2+)](mito) responses in perinuclear and distal mitochondria. CGP-37157 and Ru360 significantly enhanced histamine-induced [Ca(2+)](cyt). TNF-α and IL-13 both increased [Ca(2+)](cyt), which was associated with decreased [Ca(2+)](mito) in the case of TNF-α but not IL-13. The effects of TNF-α on both [Ca(2+)](cyt) and [Ca(2+)](mito) were affected by CGP-37157 but not by Ru360. Overall, these data demonstrate that in human ASM cells, mitochondria buffer [Ca(2+)](cyt) after agonist stimulation and its enhancement by inflammation. The differential regulation of [Ca(2+)](mito) in different parts of ASM cells may serve to locally regulate Ca(2+) fluxes from intracellular sources versus the plasma membrane as well as respond to differential energy demands at these sites. We propose that such differential mitochondrial regulation, and its disruption, may play a role in airway hyperreactivity in diseases such as asthma, where [Ca(2+)](cyt) is increased.

BK channel β1 subunits regulate airway contraction secondary to M2 muscarinic acetylcholine receptor mediated depolarization

The large conductance calcium- and voltage-activated potassium channel (BK channel) and its smooth muscle-specific β1 subunit regulate excitation–contraction coupling in many types of smooth muscle cells. However, the relative contribution of BK channels to control of M2- or M3-muscarinic acetylcholine receptor mediated airway smooth muscle contraction is poorly understood. Previously, we showed that knockout of the BK channel β1 subunit enhances cholinergic-evoked trachea contractions. Here, we demonstrate that the enhanced contraction of the BK β1 knockout can be ascribed to a defect in BK channel opposition of M2 receptor-mediated contractions. Indeed, the enhanced contraction of β1 knockout is eliminated by specific M2 receptor antagonism. The role of BK β1 to oppose M2 signalling is evidenced by a greater than fourfold increase in the contribution of L-type voltage-dependent calcium channels to contraction that otherwise does not occur with M2 antagonist or with β1 containing BK channels. The mechanism through which BK channels oppose M2-mediated recruitment of calcium channels is through a negative shift in resting voltage that offsets, rather than directly opposes, M2-mediated depolarization. The negative shift in resting voltage is reduced to similar extents by BK β1 knockout or by paxilline block of BK channels. Normalization of β1 knockout baseline voltage with low external potassium eliminated the enhanced M2-receptor mediated contraction. In summary, these findings indicate that an important function of BK/β1 channels is to oppose cholinergic M2 receptor-mediated depolarization and activation of calcium channels by restricting excitation–contraction coupling to more negative voltage ranges.

A mathematical analysis of agonist- and KCl-induced Ca(2+) oscillations in mouse airway smooth muscle cells

Airway hyperresponsiveness is a major characteristic of asthma and is generally ascribed to excessive airway narrowing associated with the contraction of airway smooth muscle cells (ASMCs). ASMC contraction is initiated by a rise in intracellular calcium concentration ([Ca(2+)](i)), observed as oscillatory Ca(2+) waves that can be induced by either agonist or high extracellular K(+) (KCl). In this work, we present a model of oscillatory Ca(2+) waves based on experimental data that incorporate both the inositol trisphosphate receptor and the ryanodine receptor. We then combined this Ca(2+) model and our modified actin-myosin cross-bridge model to investigate the role and contribution of oscillatory Ca(2+) waves to contractile force generation in mouse ASMCs. The model predicts that: 1), the difference in behavior of agonist- and KCl-induced Ca(2+) waves results principally from the fact that the sarcoplasmic reticulum is depleted during agonist-induced oscillations, but is overfilled during KCl-induced oscillations; 2), regardless of the order in which agonist and KCl are added into the cell, the resulting [Ca(2+)](i) oscillations will always be the short-period, agonist-induced-like oscillations; and 3), both the inositol trisphosphate receptor and the ryanodine receptor densities are higher toward one end of the cell. In addition, our results indicate that oscillatory Ca(2+) waves generate less contraction than whole-cell Ca(2+) oscillations induced by the same agonist concentration. This is due to the spatial inhomogeneity of the receptor distributions.

Towards a virtual lung: multi-scale, multi-physics modelling of the pulmonary system

The essential function of the lung, gas exchange, is dependent on adequate matching of ventilation and perfusion, where air and blood are delivered through complex branching systems exposed to regionally varying transpulmonary and transmural pressures. Structure and function in the lung are intimately related, yet computational models in pulmonary physiology usually simplify or neglect structure. The geometries of the airway and vascular systems and their interaction with parenchymal tissue have an important bearing on regional distributions of air and blood, and therefore on whole lung gas exchange, but this has not yet been addressed by modelling studies. Models for gas exchange have typically incorporated considerable detail at the level of chemical reactions, with little thought for the influence of structure. To date, relatively little attention has been paid to modelling at the cellular or subcellular level in the lung, or to linking information from the protein structure/interaction and cellular levels to the operation of the whole lung. We review previous work in developing anatomically based models of the lung, airways, parenchyma and pulmonary vasculature, and some functional studies in which these models have been used. Models for gas exchange at several spatial scales are briefly reviewed, and the challenges and benefits from modelling cellular function in the lung are discussed.

Quantitative cellular description of gastric slow wave activity

Interstitial cells of Cajal (ICC) are responsible for the spontaneous and omnipresent electrical activity in the stomach. A quantitative description of the intracellular processes whose coordinated activity is believed to generate electrical slow waves has been developed and is presented here. In line with recent experimental evidence, the model describes how the interplay between the mitochondria and the endoplasmic reticulum in cycling intracellular Ca(2+) provides the primary regulatory signal for the initiation of the slow wave. The major ion channels that have been identified as influencing slow wave activity have been modeled according to data obtained from isolated ICC. The model has been validated by comparing the simulated profile of the slow waves with experimental recordings and shows good correspondence in terms of frequency, amplitude, and shape in both control and pharmacologically altered conditions.

Spatio-temporal modelling explains the effect of reduced plasma membrane Ca2+ efflux on intracellular Ca2+ oscillations in hepatocytes

In many non-excitable eukaryotic cells, including hepatocytes, Ca(2+) oscillations play a key role in intra- and intercellular signalling, thus regulating many cellular processes from fertilisation to death. Therefore, understanding the mechanisms underlying these oscillations, and consequently understanding how they may be regulated, is of great interest. In this paper, we study the influence of reduced Ca(2+) plasma membrane efflux on Ca(2+) oscillations in hepatocytes. Our previous experiments with carboxyeosin show that a reduced plasma membrane Ca(2+) efflux increases the frequency of Ca(2+) oscillations, but does not affect the duration of individual transients. This phenomenon can be best explained by taking into account not only the temporal, but also the spatial dynamics underlying the generation of Ca(2+) oscillations in the cell. Here we divide the cell into a grid of elements and treat the Ca(2+) dynamics as a spatio-temporal phenomenon. By converting an existing temporal model into a spatio-temporal one, we obtain theoretical predictions that are in much better agreement with the experimental observations.