Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations

Regulatory disturbances in the dynamical signaling systems of C a 2 + and NO in fibroblasts cause fibrotic disorders

Studying the calcium dynamics within a fibroblast cell individually has provided only a restricted understanding of its functions. However, research efforts focusing on systems biology approaches for such investigations have been largely neglected by researchers until now. Fibroblast cells rely on signaling from calcium ( C a 2 + ) and nitric oxide (NO) to maintain their physiological functions and structural stability. Various studies have demonstrated the correlation between NO and the control of C a 2 + dynamics in cells. However, there is currently no existing model to assess the disruptions caused by various factors in regulatory dynamics, potentially resulting in diverse fibrotic disorders. A mathematical model has been developed to investigate the effects of changes in parameters such as buffer, receptor, sarcoplasmic endoplasmic reticulum C a 2 + -ATPase (SERCA) pump, and source influx on the regulation and dysregulation of spatiotemporal calcium and NO dynamics in fibroblast cells. This model is based on a system of reaction-diffusion equations, and numerical simulations are conducted using the finite element method. Disturbances in key processes related to calcium and nitric oxide, including source influx, buffer mechanism, SERCA pump, and inositol trisphosphate ( I P 3 ) receptor, may contribute to deregulation in the calcium and NO dynamics within fibroblasts. The findings also provide new insights into the extent and severity of disorders resulting from alterations in various parameters, potentially leading to deregulation and the development of fibrotic disease.

Modelling Cross Talk in the Spatiotemporal System Dynamics of Calcium, IP3 and Nitric Oxide in Neuron Cells

The bioenergetic system of calcium ([Ca2+]), inositol 1, 4, 5-trisphophate (IP3) and nitric oxide (NO) regulate the diverse mechanisms in neurons. The dysregulation in any or all of the calcium, IP3 and nitric oxide dynamics may cause neurotoxicity and cell death. Few studies are noted in the literature on the interactions of two systems like [Ca2+] with IP3 and [Ca2+] with nitric oxide in neuron cells, which gives limited insights into regulatory and dysregulatory processes in neuron cells. But, no study is available on the cross talk in dynamics of three systems [Ca2+], IP3 and NO in neurons. Thus, the cross talk in the system dynamics of [Ca2+], IP3 and NO regulation processes in neurons have been studied using mathematical model. The two-way feedback process between [Ca2+] and IP3 and two-way feedback process between [Ca2+] and NO through cyclic guanosine monophosphate (cGMP) with plasmalemmal [Ca2+]-ATPase (PMCA) have been incorporated in the proposed model. This coupling handles the indirect two-way feedback process between IP3 and nitric oxide in neuronal cells automatically. The numerical outcomes were acquired by employing the finite element method (FEM) with the Crank-Nicholson scheme (CNS). The present model incorporating the sodium-calcium exchanger (NCX) and voltage-gated calcium channel (VGCC) provides novel insights into the various regulatory and dysregulatory processes due to buffer, IP3-receptor, ryanodine receptor, cGMP kinetics through PMCA channel, etc. and their impacts on the interactive spatiotemporal system dynamics of [Ca2+], IP3 and NO in neurons. It is concluded that the behavior of different crucial mechanisms is quite different for interactions of two systems of [Ca2+] and NO and the interactions of three systems of [Ca2+], IP3 and nitric oxide in neuronal cell due to mutual regulatory adjustments. The association of several neurological disorders with the alterations in calcium, IP3 and NO has been explored in neurons.

Model of Calcium Dynamics Regulating [Formula: see text], ATP and Insulin Production in a Pancreatic [Formula: see text]-Cell

The calcium signals regulate the production and secretion of many signaling molecules like inositol trisphosphate ([Formula: see text]) and adenosine triphosphate (ATP) in various cells including pancreatic [Formula: see text]-cells. The calcium signaling mechanisms regulating [Formula: see text], ATP and insulin responsible for various functions of [Formula: see text]-cells are still not well understood. Any disturbance in these mechanisms can alter the functions of [Formula: see text]-cells leading to diabetes and metabolic disorders. Therefore, a mathematical model is proposed by incorporating the reaction-diffusion equation for calcium dynamics and a system of first-order differential equations for [Formula: see text], ATP-production and insulin secretion with initial and boundary conditions. The model incorporates the temporal dependence of [Formula: see text]-production and degradation, ATP production and insulin secretion on calcium dynamics in a [Formula: see text]-cell. The piecewise linear finite element method has been used for the spatial dimension and the Crank-Nicolson scheme for the temporal dimension to obtain numerical results. The effect of changes in source influxes and buffers on calcium dynamics and production of [Formula: see text], ATP and insulin levels in a [Formula: see text]-cell has been analyzed. It is concluded that the dysfunction of source influx and buffers can cause significant variations in calcium levels and dysregulation of [Formula: see text], ATP and insulin production, which can lead to various metabolic disorders, diabetes, obesity, etc. The proposed model provides crucial information about the changes in mechanisms of calcium dynamics causing proportionate disturbances in [Formula: see text], ATP and insulin levels in pancreatic cells, which can be helpful for devising protocols for diagnosis and treatment of various metabolic diseases.

Bound Ca2+ moves faster and farther from single open channels than free Ca2

A concept of Ca2+ nanodomains established in the cytoplasm after opening single-calcium channels helps mechanistically understand the physiological mechanisms of Ca2+ signaling. It predicts standing gradients of cytoplasmic free Ca2+ around single channels in the plasma membrane. The fate of bound Ca2+ attracted much less attention. This study aimed to examine the profiles of Ca2+ bound to low-mobility buffers such as bulky Ca2+-binding proteins. The solution of non-linear PDEs for an immobile buffer predicts fast decay of free [Ca2+] from the channel lumen and the traveling wave for bound Ca2+. For low-mobility buffers like calmodulin, the calculated profiles of free and bound Ca2+ are similar. Theoretical predictions are tested by imaging 1D profiles of Ca2+ bound to low-mobility fluo-4-dextran. The traveling waves of bound Ca2+ are observed that develop during the opening of single channels. The findings tempt to propose that Ca2+ signaling may not be solely related by the absolute free [Ca2+] at the sensor location, which is extremely localized, but determined by the time when a wave of bound Ca2+ reaches a threshold needed for sensor activation.

An integrate-and-fire approach to Ca2+ signaling. Part II: Cumulative refractoriness

Inositol 1,4,5-trisphosphate-induced Ca2+ signaling is a second messenger system used by almost all eukaryotic cells. The agonist concentration stimulating Ca2+ signals is encoded in the frequency of a Ca2+ concentration spike sequence. When a cell is stimulated, the interspike intervals (ISIs) often show a distinct transient during which they gradually increase, a system property we refer to as cumulative refractoriness. We extend a previously published stochastic model to include the Ca2+ concentration in the intracellular Ca2+ store as a slow adaptation variable. This model can reproduce both stationary and transient statistics of experimentally observed ISI sequences. We derive approximate expressions for the mean and coefficient of variation of the stationary ISIs. We also consider the response to the onset of a constant stimulus and estimate the length of the transient and the strength of the adaptation of the ISI. We show that the adaptation sets the coefficient of variation in agreement with current ideas derived from experiments. Moreover, we explain why, despite a pronounced transient behavior, ISI correlations can be weak, as often observed in experiments. Finally, we fit our model to reproduce the transient statistics of experimentally observed ISI sequences in stimulated HEK cells. The fitted model is able to qualitatively reproduce the relationship between the stationary interval correlations and the number of transient intervals, as well as the strength of the ISI adaptation. We also find positive correlations in the experimental sequence that cannot be explained by our model.

Impact of Interdependent Ca2+ and IP3 Dynamics On ATP Regulation in A Fibroblast Model

The vital participation of Ca2+ in human organ functions such as muscular contractions, heartbeat, brain functionality, skeletal activity, etc, motivated the scientists to thoroughly research the mechanisms of calcium (Ca2+) signalling in distinct human cells. Ca2+, inositol triphosphate (IP3), and adenosine triphosphate (ATP) play important roles in cell signaling and physiological processes. ATP and its derivatives are hypothesized to be important in the pathogenic process that leads to fibrotic illnesses like fibrosis. Fluctuations in Ca2+ and IP3 in a fibroblast cell influence ATP production. To date, no evidence of coupled Ca2+ and IP3 mechanics regulating ATP generation in a fibroblast cell during fibrotic disease has been found. The current work suggests an integrated mechanism for Ca2+ and IP3 dynamics in a fibroblast cell that regulates ATP generation. Simulation has been carried out using the finite element approach. The mechanics of interdependent systems findings vary dramatically from the results of basic independent system mechanics and give fresh information about the two systems' activities. The numerical results provide new insights into the impacts of disturbances in source influx, the serca pump, and buffers on interdependent Ca2+ and IP3 dynamics and ATP synthesis in a fibroblast cell. According to the findings of this study, fibrotic disorders cannot be attributed solely to disruptions in the processes of calcium signaling mechanics but also to disruptions in IP3 regulation mechanisms affecting the regulation of calcium in the fibroblast cell and ATP release.

Study of disorders in regulatory spatiotemporal neurodynamics of calcium and nitric oxide

Experimental studies have reported the dependence of nitric oxide (NO) on the regulation of neuronal calcium ([Ca2+]) dynamics in neurons. But, there is no model available to estimate the disorders caused by various parameters in their regulatory dynamics leading to various neuronal disorders. A mathematical model to analyze the impacts due to alterations in various parameters like buffer, ryanodine receptor, serca pump, source influx, etc. leading to regulation and dysregulation of the spatiotemporal calcium and NO dynamics in neuron cells is constructed using a system of reaction-diffusion equations. The numerical simulation is performed with the finite element approach. The disturbances in the different constitutive processes of [Ca2+] and nitric oxide including source influx, buffer mechanism, ryanodine receptor, serca pump, IP3 receptor, etc. can be responsible for the dysregulation in the [Ca2+] and NO dynamics in neurons. Also, the results reveal novel information about the magnitude and intensity of disorders in response to a range of alterations in various parameters of this neuronal dynamics, which can cause dysregulation leading to neuronal diseases like Parkinson's, cerebral ischemia, trauma, etc.

Analysis of dim-light responses in rod and cone photoreceptors with altered calcium kinetics

Rod and cone photoreceptors in the retina of vertebrates are the primary sensory neurons underlying vision. They convert light into an electrical current using a signal transduction pathway that depends on Ca[Formula: see text] feedback. It is known that manipulating the Ca[Formula: see text] kinetics affects the response shape and the photoreceptor sensitivity, but a precise quantification of these effects remains unclear. We have approached this task in mouse retina by combining numerical simulations with mathematical analysis. We consider a parsimonious phototransduction model that incorporates negative Ca[Formula: see text] feedback onto the synthesis of cyclic GMP, and fast buffering reactions to alter the Ca[Formula: see text] kinetics. We derive analytic results for the photoreceptor functioning in sufficiently dim light conditions depending on the photoreceptor type. We exploit these results to obtain conceptual and quantitative insight into how response waveform and amplitude depend on the underlying biophysical processes and the Ca[Formula: see text] feedback. With a low amount of buffering, the Ca[Formula: see text] concentration changes in proportion to the current, and responses to flashes of light are monophasic. With more buffering, the change in the Ca[Formula: see text] concentration becomes delayed with respect to the current, which gives rise to a damped oscillation and a biphasic waveform. This shows that biphasic responses are not necessarily a manifestation of slow buffering reactions. We obtain analytic approximations for the peak flash amplitude as a function of the light intensity, which shows how the photoreceptor sensitivity depends on the biophysical parameters. Finally, we study how changing the extracellular Ca[Formula: see text] concentration affects the response.

Calcium messages in flagella are faster than messenger particles

Calcium is one of the most versatile messengers for intracellular signaling. In the case of cilia and flagella calcium has the central role in transfer of communications between extracellular stimuli and intracellular formation of frequency modulated signal and their deciphering by target proteins. In this paper, the diffusion of fluorescently or otherwise tagged and un-tagged Ca2＋ particles is analyzed by solving the system of pertaining reaction-diffusion equations. We used Fourier transform tools to get asymptotic eigenfunctions for tagged (un-tagged) free and buffered Ca2＋ ions. We made some numerical estimations for diffusion coefficients corroborating the fact that messages diffuse faster than Ca2＋ messengers. From the best of our knowledge, this is the first time that Ca2＋ signaling in living cells is biophysically elaborated within the framework of model presented here. We suggest the experimental assay on the basis of radioactive Ca2＋ as tagged probe.

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

On the Ca2＋ elevation in vascular endothelial cells due to inositol trisphosphate-sensitive store receptors activation: A data-driven modeling approach

Agonist-induced Ca2+ signaling is essential for the regulation of many vital functions in endothelial cells (ECs). A broad range of stimuli elevate the cytosolic Ca2＋ concentration by promoting a pathway mediated by inositol 1,4,5 trisphosphate (IP3) which causes Ca2＋ release from intracellular stores. Despite its importance, there are very few studies focusing on the quantification of such dynamics in the vascular endothelium. Here, by using data from isolated ECs, we established a minimalistic modeling framework able to quantitatively capture the main features (averaged over a cell population) of the cytosolic Ca2＋ response to different IP3 stimulation levels. A suitable description of Ca2＋-regulatory function of inositol 1,4,5 trisphosphate receptors (IP3Rs) and corresponding parameter space are identified by comparing the different model variants against experimental mean population data. The same approach is used to numerically assess the relevance of cytosolic Ca2＋ buffering, as well as Ca2＋ store IP3-sensitivity in the overall cell dynamics. The variability in the dynamics' features observed across the population can be explained (at least in part) through variation of certain model parameters (such as buffering capacity or Ca2＋ store sensitivity to IP3). The results, in terms of experimental fitting and validation, support the proposed minimalistic model as a reference framework for the quantification of the EC Ca2＋ dynamics induced by IP3Rs activation.

Simulation of biochemical dynamics of [Formula: see text] and [Formula: see text] in fibroblast cell

Calcium dynamics is not only responsible for maintaining the framework and functions of the cell but also plays a role in the dynamics of other biochemical systems in the cell. Phospholipase C-[Formula: see text] l ([Formula: see text]) has a crucial role in the function of fibroblast cells. Experiments have shown that [Formula: see text] and [Formula: see text] have interdependent dynamics in fibroblast cells. However, no reaction-diffusion model exists for the two-way feedback system dynamics of [Formula: see text] and [Formula: see text] in fibroblasts till date. The computational model is designed to investigate the impact of variations in several processes, such as the [Formula: see text] pump, buffer process, source inflow, etc., on the system dynamics of [Formula: see text] and [Formula: see text] in fibroblast cells. The computational findings are obtained using finite element techniques, and the consequences of dysregulation in various processes on the spatiotemporal calcium and [Formula: see text] dynamics in fibroblasts are investigated. The results lead to the conclusion that the effects of buffer, source influx, diffusion, and [Formula: see text] pump can cause fluctuations in the dynamics of [Formula: see text] and [Formula: see text] in fibroblasts. Disruptions in these constitutive processes can result in changes in the dynamics of calcium and [Formula: see text]. Thus, the current model provides new/novel information regarding the precise dysregulatory constitutive systems that regulate calcium and [Formula: see text] kinetics, such as source inflow, diffusion, [Formula: see text], and buffer, can be responsible for excessive calcium and [Formula: see text] concentrations leading to fibrotic illnesses such as cancer and fibrosis.

IP3 receptors: An "elementary" journey from structure to signals

Inositol 1,4,5-trisphosphate receptors (IP3Rs) are large tetrameric channels which sit mostly in the membrane of the endoplasmic reticulum (ER) and mediate Ca2+ release from intracellular stores in response to extracellular stimuli in almost all cells. Dual regulation of IP3Rs by IP3 and Ca2+ itself, upstream "licensing", and the arrangement of IP3Rs into small clusters in the ER membrane, allow IP3Rs to generate spatially and temporally diverse Ca2+ signals. The characteristic biphasic regulation of IP3Rs by cytosolic Ca2+ concentration underpins regenerative Ca2+ signals by Ca2+-induced Ca2+-release, while also preventing uncontrolled explosive Ca2+ release. In this way, cells can harness a simple ion such as Ca2+ as a near-universal intracellular messenger to regulate diverse cellular functions, including those with conflicting outcomes such as cell survival and cell death. High-resolution structures of the IP3R bound to IP3 and Ca2+ in different combinations have together started to unravel the workings of this giant channel. Here we discuss, in the context of recently published structures, how the tight regulation of IP3Rs and their cellular geography lead to generation of "elementary" local Ca2+ signals known as Ca2+ "puffs", which form the fundamental bottleneck through which all IP3-mediated cytosolic Ca2+ signals must first pass.

Cellular nitric oxide synthesis is affected by disorders in the interdependent [Formula: see text] and [Formula: see text] dynamics during cystic fibrosis disease

Calcium ([Formula: see text]), inositol trisphosphate ([Formula: see text]), and nitric oxide (NO) signaling are essential to maintain the structural integrity and physiological activity of fibroblast cells. The accumulation of excess quantity of NO for longer periods can lead to a variety of fibrotic disorders, including heart disease, penile fibrosis in Peyronie's disease (PD), and cystic fibrosis. The dynamics of these three signaling processes and their interdependence in fibroblast cells are not clearly known to date. A systems biology model is proposed using reaction-diffusion equations for calcium, [Formula: see text], and calcium-dependent NO synthesis in fibroblast cells. The finite element method (FEM) is used to examine [Formula: see text], [Formula: see text], and NO regulation and dysregulation in cells. The results throw light on the conditions that disturb the coupled [Formula: see text] and [Formula: see text] dynamics and the influence of these factors on the levels of NO concentration in the fibroblast cell. The findings suggest that changes in source inflow, buffers, and diffusion coefficient might induce an increase or reduction in nitric oxide and [Formula: see text] synthesis, resulting in fibroblast cell diseases. Furthermore, the findings provide new information regarding the size and intensity of diseases in response to changes in several factors of their dynamics, which has been linked to the development of cystic fibrosis and cancer. This knowledge could be valuable for developing novel approaches to the diagnosis of diseases and therapies for various disorders of fibroblast cells.

Mechanistic insights of neuronal calcium and IP3 signaling system regulating ATP release during ischemia in progression of Alzheimer's disease

The mechanisms of calcium ([Ca²⁺]) signaling in various human cells have been widely analyzed by scientists due to its crucial role in human organs like the heartbeat, muscle contractions, bone activity, brain functionality, etc. No study is reported for interdependent [Ca²⁺] and IP₃ mechanics regulating the release of ATP in neuron cells during Ischemia in Alzheimer's disease advancement. In the present investigation, a finite element method (FEM) is framed to explore the interdependence of spatiotemporal [Ca²⁺] and IP₃ signaling mechanics and its role in ATP release during Ischemia as well as in the advancement of Alzheimer's disorder in neuron cells. The results provide us insights of the mutual spatiotemporal impacts of [Ca²⁺] and IP₃ mechanics as well as their contributions to ATP release during Ischemia in neuron cells. The results obtained for the mechanics of interdependent systems differ significantly from the results of simple independent system mechanics and provide new information about the processes of the two systems. From this study, it is concluded that neuronal disorders cannot only be simply attributed to the disturbance caused directly in the processes of calcium signaling mechanics, but also to the disturbances caused in IP₃ regulation mechanisms impacting the calcium regulation in the neuron cell and ATP release.

An integrate-and-fire approach to Ca2+ signaling. Part I: Renewal model

In computational neuroscience integrate-and-fire models capture the spike generation by a subthreshold dynamics supplemented by a simple fire-and-reset rule; they allow for a numerically efficient and analytically tractable description of stochastic single cell as well as network dynamics. Stochastic spiking is also a prominent feature of Ca2+ signaling which suggests to adopt the integrate-and-fire approach for this fundamental biophysical process. The model introduced here consists of two components describing 1) activity of clusters of inositol-trisphosphate receptor channels and 2) dynamics of the global Ca2+ concentrations in the cytosol. The cluster dynamics is given in terms of a cyclic Markov chain, capturing the puff, i.e., the punctuated release of Ca2+ from intracellular stores. The cytosolic Ca2+ concentration is described by an integrate-and-fire dynamics driven by the puff current. For the cyclic Markov chain we derive expressions for the statistics of the interpuff interval, the single-puff strength and the puff current assuming constant cytosolic Ca2+. The latter condition is often well approximated because cytosolic Ca2+ varies much slower than the cluster activity does. Furthermore, because the detailed two-component model is numerically expensive to simulate and difficult to treat analytically, we develop an analytical framework to approximate the driving puff current of the stochastic cytosolic Ca2+ dynamics by a temporally uncorrelated Gaussian noise. This approximation reduces our two-component system to an integrate-and-fire model with a nonlinear drift function and a multiplicative Gaussian white noise, a model that is known to generate a renewal spike train, i.e., a point process with statistically independent interspike intervals. The model allows for fast numerical simulations, permits to derive analytical expressions for the rate of Ca2+ spiking and the coefficient of variation of the interspike interval, and to approximate the interspike interval density and the spike train power spectrum. Comparison of these statistics to experimental data is discussed.

Effect of disturbances in neuronal calcium and IP3 dynamics on β-amyloid production and degradation

Overproduction and accumulation of β-amyloid and its improper clearance can cause neurotoxicity leading to Alzheimer's disease. The production and degradation of β-amyloid depend on the calcium ([Ca2+]) and IP3 dynamics in the nerve cells. Thus, there is a need to understand the impacts of disturbances in the processes of [Ca2+] and IP3 dynamics on β-amyloid production and its degradation. Here, a model is proposed to investigate the role of [Ca2+] and IP3 dynamics on β-amyloid production and degradation. The problem is formulated in terms of the initial boundary value problem involving the system of two reaction-diffusion equations respectively for [Ca2+] and IP3 in the nerve cell. The solution is obtained by employing the Finite element approach. The numerical results are used to analyze the impact of various mechanisms of calcium and IP3 dynamics on β-amyloid production and degradation in a neuron cell. The results indicate that disturbances in any of the constitutive processes of interdependent calcium and IP3 dynamics like source influx, buffering, serca pump, and IP3 dynamics, etc. can cause dynamic changes in β-amyloid production and degradation, which in turn can be the cause of neurotoxicity and neuronal disorders like Alzheimer's disease. Thus, the relationships obtained by the proposed model among various mechanisms can be useful in addressing the challenges of identifying specific constitutive processes causing neuronal disorders like Alzheimer's disease, etc., and developing the framework for their diagnosis and treatment.

Mathematical modeling of intracellular calcium in presence of receptor: a homeostatic model for endothelial cell

Calcium is a ubiquitous molecule and second messenger that regulates many cellular functions ranging from exocytosis to cell proliferation at different time scales. In the vasculature, a constant adenosine triphosphate (ATP) concentration is maintained because of ATP released by red blood cells (RBCs). These ATP molecules continuously react with purinergic receptors on the surface of endothelial cells (ECs). Consequently, a cascade of chemical reactions are triggered that result in a transient cytoplasmic calcium (Ca[Formula: see text]), followed by return to its basal concentration. The mathematical models proposed in the literature are able to reproduce the transient peak. However, the trailing concentration is always higher than the basal cytoplasmic Ca[Formula: see text] concentrations, and the Ca[Formula: see text] concentration in endoplasmic reticulum (ER) remains lower than its initial concentration. This means that the intracellular homeostasis is not recovered. We propose, herein, a minimal model of calcium kinetics. We find that the desensitization of EC surface receptors due to phosphorylation and recycling plays a vital role in maintaining calcium homeostasis in the presence of a constant stimulus (ATP). The model is able to capture several experimental observations such as refilling of Ca[Formula: see text] in the ER, variation of cytoplasmic Ca[Formula: see text] transient peak in ECs, the resting cytoplasmic Ca[Formula: see text] concentration, the effect of removing ATP from the plasma on Ca[Formula: see text] homeostasis, and the saturation of cytoplasmic Ca[Formula: see text] transient peak with increase in ATP concentration. Direct confrontation with several experimental results is conducted. This work paves the way for systematic studies on coupling between blood flow and chemical signaling, and should contribute to a better understanding of the relation between (patho)physiological conditions and Ca[Formula: see text] kinetics.

Effect of Buffers with Multiple Binding Sites on Calcium Waves

The existence and properties of intracellular waves of increased free cytoplasmic calcium concentration (calcium waves) are strongly affected by the binding and unbinding of calcium ions to a multitude of different buffers in the cell. These buffers can be mobile or immobile and, in general, have multiple binding sites that are not independent. Previous theoretical studies have focused on the case when each buffer molecule binds a single calcium ion. In this study, we analyze how calcium waves are affected by calcium buffers with two non-independent binding sites, and show that the interactions between the calcium binding sites can result in the emergence of new behaviors. In particular, for certain combinations of kinetic parameters, the profiles of buffer molecules with one calcium ion bound can be non-monotone.

Encoding, transmission, decoding, and specificity of calcium signals in plants

Calcium acts as a signal and transmits information in all eukaryotes. Encoding machinery consisting of calcium channels, stores, buffers, and pumps can generate a variety of calcium transients in response to external stimuli, thus shaping the calcium signature. Mechanisms for the transmission of calcium signals have been described, and a large repertoire of calcium binding proteins exist that can decode calcium signatures into specific responses. Whilst straightforward as a concept, mysteries remain as to exactly how such information processing is biochemically implemented. Novel developments in imaging technology and genetically encoded sensors (such as calcium indicators), in particular for multi-signal detection, are delivering exciting new insights into intra- and intercellular calcium signaling. Here, we review recent advances in characterizing the encoding, transmission, and decoding mechanisms, with a focus on long-distance calcium signaling. We present technological advances and computational frameworks for studying the specificity of calcium signaling, highlight current gaps in our understanding and propose techniques and approaches for unravelling the underlying mechanisms.

Understanding the Dynamics of the Transient and Permanent Opening Events of the Mitochondrial Permeability Transition Pore with a Novel Stochastic Model

The mitochondrial permeability transition pore (mPTP) is a non-selective pore in the inner mitochondrial membrane (IMM) which causes depolarization when it opens under conditions of oxidative stress and high concentrations of Ca2+. In this study, a stochastic computational model was developed to better understand the dynamics of mPTP opening and closing associated with elevated reactive oxygen species (ROS) in cardiomyocytes. The data modeled are from "photon stress" experiments in which the fluorescent dye TMRM (tetramethylrhodamine methyl ester) is both the source of ROS (induced by laser light) and sensor of the electrical potential difference across the IMM. Monte Carlo methods were applied to describe opening and closing of the pore along with the Hill Equation to account for the effect of ROS levels on the transition probabilities. The amplitude distribution of transient mPTP opening events, the number of transient mPTP opening events per minute in a cell, the time it takes for recovery after transient depolarizations in the mitochondria, and the change in TMRM fluorescence during the transition from transient to permanent mPTP opening events were analyzed. The model suggests that mPTP transient open times have an exponential distribution that are reflected in TMRM fluorescence. A second multiple pore model in which individual channels have no permanent open state suggests that 5-10 mPTP per mitochondria would be needed for sustained mitochondrial depolarization at elevated ROS with at least 1 mPTP in the transient open state.

Cardiac Alternans Occurs through the Synergy of Voltage- and Calcium-Dependent Mechanisms

Cardiac alternans is characterized by alternating weak and strong beats of the heart. This signaling at the cellular level may appear as alternating long and short action potentials (APs) that occur in synchrony with alternating large and small calcium transients, respectively. Previous studies have suggested that alternans manifests itself through either a voltage dependent mechanism based upon action potential restitution or as a calcium dependent mechanism based on refractoriness of calcium release. We use a novel model of cardiac excitation-contraction (EC) coupling in the rat ventricular myocyte that includes 20,000 calcium release units (CRU) each with 49 ryanodine receptors (RyR2s) and 7 L-type calcium channels that are all stochastically gated. The model suggests that at the cellular level in the case of alternans produced by rapid pacing, the mechanism requires a synergy of voltage- and calcium-dependent mechanisms. The rapid pacing reduces AP duration and magnitude reducing the number of L-type calcium channels activating individual CRUs during each AP and thus increases the population of CRUs that can be recruited stochastically. Elevated myoplasmic and sarcoplasmic reticulum (SR) calcium, [Ca2+]myo and [Ca2+]SR respectively, increases ryanodine receptor open probability (Po) according to our model used in this simulation and this increased the probability of activating additional CRUs. A CRU that opens in one beat is less likely to open the subsequent beat due to refractoriness caused by incomplete refilling of the junctional sarcoplasmic reticulum (jSR). Furthermore, the model includes estimates of changes in Na+ fluxes and [Na+]i and thus provides insight into how changes in electrical activity, [Na+]i and sodium-calcium exchanger activity can modulate alternans. The model thus tracks critical elements that can account for rate-dependent changes in [Na+]i and [Ca2+]myo and how they contribute to the generation of Ca2+ signaling alternans in the heart.

Finite element analysis predicts Ca2+ microdomains within tubular-sarcoplasmic reticular junctions of amphibian skeletal muscle

A finite element analysis modelled diffusional generation of steady-state Ca²⁺ microdomains within skeletal muscle transverse (T)-tubular-sarcoplasmic reticular (SR) junctions, sites of ryanodine receptor (RyR)-mediated SR Ca²⁺ release. It used established quantifications of sarcomere and T-SR anatomy (radial diameter [Formula: see text]; axial distance [Formula: see text]). Its boundary SR Ca²⁺ influx densities,[Formula: see text], reflected step impositions of influxes, [Formula: see text] deduced from previously measured Ca²⁺ signals following muscle fibre depolarization. Predicted steady-state T-SR junctional edge [Ca²⁺], [Ca²⁺]_edge, matched reported corresponding experimental cytosolic [Ca²⁺] elevations given diffusional boundary efflux [Formula: see text] established cytosolic Ca²⁺ diffusion coefficients [Formula: see text] and exit length [Formula: see text]. Dependences of predicted [Ca²⁺]_edge upon [Formula: see text] then matched those of experimental [Ca²⁺] upon Ca²⁺ release through their entire test voltage range. The resulting model consistently predicted elevated steady-state T-SR junctional ~ µM-[Ca²⁺] elevations radially declining from maxima at the T-SR junction centre along the entire axial T-SR distance. These [Ca²⁺] heterogeneities persisted through 10⁴- and fivefold, variations in D and w around, and fivefold reductions in d below, control values, and through reported resting muscle cytosolic [Ca²⁺] values, whilst preserving the flux conservation ([Formula: see text] condition, [Formula: see text]. Skeletal muscle thus potentially forms physiologically significant ~ µM-[Ca²⁺] T-SR microdomains that could regulate cytosolic and membrane signalling molecules including calmodulin and RyR, These findings directly fulfil recent experimental predictions invoking such Ca²⁺ microdomains in observed regulatory effects upon Na⁺ channel function, in a mechanism potentially occurring in similar restricted intracellular spaces in other cell types.

Stationary Ca2+ nanodomains in the presence of buffers with two binding sites

We examine closed-form approximations for the equilibrium Ca²⁺ and buffer concentrations near a point Ca²⁺ source representing a Ca²⁺ channel, in the presence of a mobile buffer with two Ca²⁺ binding sites activated sequentially and possessing distinct binding affinities and kinetics. This allows us to model the impact on Ca²⁺ nanodomains of realistic endogenous Ca²⁺ buffers characterized by cooperative Ca²⁺ binding, such as calretinin. The approximations we present involve a combination or rational and exponential functions, whose parameters are constrained using the series interpolation method that we recently introduced for the case of simpler Ca²⁺ buffers with a single Ca²⁺ binding site. We conduct extensive parameter sensitivity analysis and show that the obtained closed-form approximations achieve reasonable qualitative accuracy for a wide range of buffer's Ca²⁺ binding properties and other relevant model parameters. In particular, the accuracy of the derived approximants exceeds that of the rapid buffering approximation in large portions of the relevant parameter space.

Buffering and total calcium levels determine the presence of oscillatory regimes in cardiac cells

Calcium oscillations and waves induce depolarization in cardiac cells which are believed to cause life-threathening arrhythimas. In this work, we study the conditions for the appearance of calcium oscillations in both a detailed subcellular model of calcium dynamics and a minimal model that takes into account just the minimal ingredients of the calcium toolkit. To avoid the effects of homeostatic changes and the interaction with the action potential we consider the somewhat artificial condition of a cell without pacing and with no calcium exchange with the extracellular medium. Both the full subcellular model and the minimal model present the same scenarios depending on the calcium load: two stationary states, one with closed ryanodine receptors (RyR) and most calcium in the cell stored in the sarcoplasmic reticulum (SR), and another, with open RyRs and a depleted SR. In between, calcium oscillations may appear. The robustness of these oscillations is determined by the amount of calsequestrin (CSQ). The lack of this buffer in the SR enhances the appearance of oscillations. The minimal model allows us to relate the stability of the oscillating state to the nullcline structure of the system, and find that its range of existence is bounded by a homoclinic and a Hopf bifurcation, resulting in a sudden transition to the oscillatory regime as the cell calcium load is increased. Adding a small amount of noise to the RyR behavior increases the parameter region where oscillations appear and provides a gradual transition from the resting state to the oscillatory regime, as observed in the subcellular model and experimentally.

In cardiac cells, calcium plays a very important role. An increase in calcium levels is the trigger used by the cell to initiate contraction. Besides, calcium modulates several transmembrane currents, affecting the cell transmembrane potential. Thus, dysregulations in calcium handling have been associated with the appearance of arrhythmias. Often, this dysregulation results in the appearance of periodic calcium waves or global oscillations, providing a pro-arrhythmic substrate. In this paper, we study the onset of calcium oscillations in cardiac cells using both a detailed subcellular model of calcium dynamics and a minimal model that takes into account the essential ingredients of the calcium toolkit. Both reproduce the main experimental results and link this behavior with the presence of different steady-state solutions and bifurcations that depend on the total amount of calcium in the cell and in the level of buffering present. We expect that this work will help to clarify the conditions under which calcium oscillations appear in cardiac myocytes and, therefore, will represent a step further in the understanding of the origin of cardiac arrhythmias.

Optical Mapping of cAMP Signaling at the Nanometer Scale

Cells relay a plethora of extracellular signals to specific cellular responses by using only a few second messengers, such as cAMP. To explain signaling specificity, cAMP-degrading phosphodiesterases (PDEs) have been suggested to confine cAMP to distinct cellular compartments. However, measured rates of fast cAMP diffusion and slow PDE activity render cAMP compartmentalization essentially impossible. Using fluorescence spectroscopy, we show that, contrary to earlier data, cAMP at physiological concentrations is predominantly bound to cAMP binding sites and, thus, immobile. Binding and unbinding results in largely reduced cAMP dynamics, which we term "buffered diffusion." With a large fraction of cAMP being buffered, PDEs can create nanometer-size domains of low cAMP concentrations. Using FRET-cAMP nanorulers, we directly map cAMP gradients at the nanoscale around PDE molecules and the areas of resulting downstream activation of cAMP-dependent protein kinase (PKA). Our study reveals that spatiotemporal cAMP signaling is under precise control of nanometer-size domains shaped by PDEs that gate activation of downstream effectors.

Efficient Approximations for Stationary Single-Channel Ca2+ Nanodomains across Length Scales

We consider the stationary solution for the Ca²⁺ concentration near a point Ca²⁺ source describing a single-channel Ca²⁺ nanodomain in the presence of a single mobile Ca²⁺ buffer with 1:1 Ca²⁺ binding. We present computationally efficient approximants that estimate stationary single-channel Ca²⁺ nanodomains with great accuracy in broad regions of parameter space. The presented approximants have a functional form that combines rational and exponential functions, which is similar to that of the well-known excess buffer approximation and the linear approximation but with parameters estimated using two novel, to our knowledge, methods. One of the methods involves interpolation between the short-range Taylor series of the free buffer concentration and its long-range asymptotic series in inverse powers of distance from the channel. Although this method has already been used to find Padé (rational-function) approximants to single-channel Ca²⁺ and buffer concentrations, extending this method to interpolants combining exponential and rational functions improves accuracy in a significant fraction of the relevant parameter space. A second method is based on the variational approach and involves a global minimization of an appropriate functional with respect to parameters of the chosen approximations. An extensive parameter-sensitivity analysis is presented, comparing these two methods with previously developed approximants. Apart from increased accuracy, the strength of these approximants is that they can be extended to more realistic buffers with multiple binding sites characterized by cooperative Ca²⁺ binding, such as calmodulin and calretinin.

Ca2+ Release via IP3 Receptors Shapes the Cardiac Ca2+ Transient for Hypertrophic Signaling

Calcium (Ca²⁺) plays a central role in mediating both contractile function and hypertrophic signaling in ventricular cardiomyocytes. L-type Ca²⁺ channels trigger release of Ca²⁺ from ryanodine receptors for cellular contraction, whereas signaling downstream of G-protein-coupled receptors stimulates Ca²⁺ release via inositol 1,4,5-trisphosphate receptors (IP₃Rs), engaging hypertrophic signaling pathways. Modulation of the amplitude, duration, and duty cycle of the cytosolic Ca²⁺ contraction signal and spatial localization have all been proposed to encode this hypertrophic signal. Given current knowledge of IP₃Rs, we develop a model describing the effect of functional interaction (cross talk) between ryanodine receptor and IP₃R channels on the Ca²⁺ transient and examine the sensitivity of the Ca²⁺ transient shape to properties of IP₃R activation. A key result of our study is that IP₃R activation increases Ca²⁺ transient duration for a broad range of IP₃R properties, but the effect of IP₃R activation on Ca²⁺ transient amplitude is dependent on IP₃ concentration. Furthermore we demonstrate that IP₃-mediated Ca²⁺ release in the cytosol increases the duty cycle of the Ca²⁺ transient, the fraction of the cycle for which [Ca²⁺] is elevated, across a broad range of parameter values and IP₃ concentrations. When coupled to a model of downstream transcription factor (NFAT) activation, we demonstrate that there is a high correspondence between the Ca²⁺ transient duty cycle and the proportion of activated NFAT in the nucleus. These findings suggest increased cytosolic Ca²⁺ duty cycle as a plausible mechanism for IP₃-dependent hypertrophic signaling via Ca²⁺-sensitive transcription factors such as NFAT in ventricular cardiomyocytes.

Modelling dendritic spines with the finite element method, investigating the impact of geometry on electric and calcic responses

Understanding the relationship between shape and function of dendritic spines is an elusive topic. Several modelling approaches have been used to investigate the interplay between spine geometry, calcium diffusion and electric signalling. We here use a second order finite element method to solve the Poisson-Nernst-Planck equations and describe electrodiffusion in dendritic spines. With this, we obtain relationships between dendritic geometry and calcic as well as electric responses to synaptic events. Our findings support the hypothesis that spine geometry plays a role shaping the electrical responses to synaptic events. Our method was also able to reveal the fine scale distribution of calcium in spines with irregular shapes.

A mathematical model for persistent post-CSD vasoconstriction

Cortical spreading depression (CSD) is the propagation of a relatively slow wave in cortical brain tissue that is linked to a number of pathological conditions such as stroke and migraine. Most of the existing literature investigates the dynamics of short term phenomena such as the depolarization and repolarization of membrane potentials or large ion shifts. Here, we focus on the clinically-relevant hour-long state of neurovascular malfunction in the wake of CSDs. This dysfunctional state involves widespread vasoconstriction and a general disruption of neurovascular coupling. We demonstrate, using a mathematical model, that dissolution of calcium that has aggregated within the mitochondria of vascular smooth muscle cells can drive an hour-long disruption. We model the rate of calcium clearance as well as the dynamical implications on overall blood flow. Based on reaction stoichiometry, we quantify a possible impact of calcium phosphate dissolution on the maintenance of F0F1-ATP synthase activity.

Two-variable nullcline analysis of ionic general equilibrium predicts calcium homeostasis in ventricular myocytes

Ventricular contraction is roughly proportional to the amount of calcium released from the Sarcoplasmic Reticulum (SR) during systole. While it is rather straightforward to measure calcium levels and contractibility under different physiological conditions, the complexity of calcium handling during systole and diastole has made the prediction of its release at steady state impossible. Here we approach the problem analyzing the evolution of intracellular and extracellular calcium fluxes during a single beat which is away from homeostatic balance. Using an in-silico subcellular model of rabbit ventricular myocyte, we show that the high dimensional nonlinear problem of finding the steady state can be reduced to a two-variable general equilibrium condition where pre-systolic calcium level in the cytosol and in the SR must fulfill simultaneously two different equalities. This renders calcium homeostasis as a problem that can be studied in terms of its equilibrium structure, leading to precise predictions of steady state from single-beat measurements. We show how changes in ion channels modify the general equilibrium, as shocks would do in general equilibrium macroeconomic models. This allows us to predict when an enhanced entrance of calcium in the cell reduces its contractibility and explain why SERCA gene therapy, a change in calcium handling to treat heart failure, might fail to improve contraction even when it successfully increases SERCA expression.

Cardiomyocytes, upon voltage excitation, release calcium, which leads to cell contraction. However, under some pathological conditions, calcium handling is impaired. Recently, SERCA gene therapy, whose aim is to improve Ca²⁺ sequestration by the Sarcoplasmic Reticulum (SR), has failed to improve the prognosis of patients with Heart Failure. This, together with recent counterintuitive results in calcium handling, has highlighted the need for a framework to understand calcium homeostasis across species and pathologies. We show here that the proper framework is a general equilibrium approach of two independent variables. The development of this framework allows us to find a possible mechanism for the failure of SERCA gene therapy even when it manages to increase Ca SERCA expression.

Modeling the mechanics of calcium regulation in T lymphocyte: A finite element method approach

Changes in cellular Ca[Formula: see text] concentration control a variety of physiological activities including hormone and neurotransmitter release, muscular contraction, synaptic plasticity, ionic channel permeability, apoptosis, enzyme activity, gene transcription and reproduction process. Spatial–temporal Ca[Formula: see text] dynamics due to Ca[Formula: see text] release, buffering and re-uptaking plays a central role in studying the Ca[Formula: see text] regulation in T lymphocytes. In most cases, Ca[Formula: see text] has its major signaling function when it is elevated in the cytosolic compartment. In this paper, a two-dimensional mathematical model to study spatiotemporal variations of intracellular Ca[Formula: see text] concentration in T lymphocyte cell is proposed and investigated. The cell is assumed to be a circular shaped geometrical domain for the representation of properties of Ca[Formula: see text] dynamics within the cell including important parameters. Ca[Formula: see text] binding proteins for the dynamics of Ca[Formula: see text] are itself buffer and other physiological parameters located in Ca[Formula: see text] stores. The model incorporates the important biophysical processes like diffusion, reaction, voltage-gated Ca[Formula: see text] channel, leak from endoplasmic reticulum (ER), efflux from cytosol to ER via sarco–ER Ca[Formula: see text] adenosine triphosphate (SERCA) pumps, buffers and Na[Formula: see text] exchanger. The proposed mathematical model is solved using a finite difference method and the finite element method. Appropriate initial and boundary conditions are incorporated in the model based on biophysical conditions of the problem. Computer simulations in MATLAB R2019b are employed to investigate mathematical models of reaction–diffusion equation. The effect of source, buffer, Na[Formula: see text]/Ca[Formula: see text] exchanger, etc. on spatial and temporal patterns of Ca[Formula: see text] in T lymphocyte has been studied with the help of numerical results. From the obtained results, it is observed that, the coordinated combination of the incorporated parameters plays a significant role in Ca[Formula: see text] regulation in T lymphocytes. ER leak and voltage-gated Ca[Formula: see text] channel provides the necessary Ca[Formula: see text] to the cell when required for its proper functioning, while on the other side buffers, SERCA pump and Na[Formula: see text]/Ca[Formula: see text] exchanger makes balance in the Ca[Formula: see text] concentration, so as to prevent the cell from death as higher concentration for longer time is harmful for the cell and can cause cell death.

Simulating diffusion from a cluster of point sources using propagation integrals

A computational methodology to simulate the diffusion of ions from point sources (e.g., ion channels) is described. The outlined approach computes the ion concentration from a cluster of many ion channels at pre-specified locations as a function of time using the theory of propagation integrals. How the channels' open/closed states evolve in time does not need to be known at the start of the simulation, but can be updated on-the-fly as the simulation goes along. The technique uses analytic formulas for the solutions of the diffusion equation for three common geometries: (1) ions diffusing from a membrane (planar symmetry); (2) ions diffusing into a narrow cleft for effective two-dimensional diffusion (cylindrical symmetry); and (3) ions diffusing into open space like the cytosol (spherical symmetry). Because these formulas are exact solutions valid for arbitrarily long timesteps, no spatial or time discretizations are necessary. The only discrete locations are where the ion concentration is computed, and the only discrete timesteps are when the channels' open/closed states are updated. Beyond pure diffusion, the technique is generalized to the Excess Buffer Approximation of ion chelation to give an analytic solution of this approximation of the full reaction/diffusion system. Both the pure diffusion and the diffusion/buffering algorithms scale linearly with the number of channels and the number of ion concentration locations.

Modeling and simulation of spatial-temporal calcium distribution in T lymphocyte cell by using a reaction-diffusion equation

T lymphocytes are white blood cells that play a central role in cell-mediated immunity. Ca2+ has its major signaling function when it is elevated in the cytosolic compartment. The free cytosolic Ca2+ dynamics plays a very important role in the activation, and fate decision process in the T lymphocytes. Here, we develop a quantitative spatio-temporal Ca2+ dynamic model which includes, the Ca2+ releasing channels ER leak and voltage-gated Ca2+ channel, buffering and re-uptaking mechanism in the T lymphocytes. In this model, the cell is represented as a circular-shaped geometrical domain. This representation introduces modeling flexibility needed for detailed representation of the properties of Ca2+ dynamics in the cell including important parameters. The proposed mathematical model is solved using a finite difference method and the finite element method. Appropriate initial and boundary conditions are incorporated in the model based on biophysical conditions of the problem. Computer simulations in MATLAB R2010a are employed to investigate mathematical models of reaction-diffusion equation. The estimation is based on reaction-diffusion equation associated with biophysical and biochemical reactions taking place in the cell. From our results, it is observed that, the coordinated combination of the incorporated parameters plays a significant role in Ca2+ regulation in T lymphocytes. ER leak and voltage-gated Ca2+ channel provides the necessary Ca2+ to the cell when required for its proper functioning, while on the other side buffers and Na+/Ca2+ exchanger makes balance in the Ca2+ concentration, so as to prevent the cell from death as higher concentration for longer time is harmful for the cell and can cause cell death. These results have been used to study the relationship of Ca2+ concentration with parameters like VGCC, Na+/Ca2+ exchanger, ER leak and buffers. The significance of the study reveals that there is a significant variation in Ca2+ profiles due to the effect of VGCC, Na+/Ca2+ exchanger, ER leak, and buffers. The results give us better insights of coordinated effect of VGCC, Na+/Ca2+ exchanger, ER leak, and buffers on Ca2+ distribution in T lymphocytes. T lymphocytes are the primary host cells to receive the viral infections which transmits the signal then to other cell types. The proper quantity of Ca2+ concentration makes T lymphocytes more active and healthier to fight the infection properly and can protect the immune system from various fatal viral infections. Thus, the application of the study lies in the field of immunology to protect a susceptible from various viral infectious diseases like HIV, HBV, HINI, etc. by strengthening the immune system. The outcomes of the study reveal that the applied finite element method is computationally very strong and effective to analyze differential equations that arise in Ca2+ dynamics.

Ca2+ signaling in T lymphocytes: the interplay of the endoplasmic reticulum, mitochondria, membrane potential, and CRAC channels on transcription factor activation

T cell receptor stimulation initiates a cascade of reactions that cause an increase in intracellular calcium (Ca2+) concentration mediated through inositol 1,4,5-trisphosphate (IP3). To understand the basic mechanisms by which the immune response in T cells is activated, it is useful to understand the signaling pathways that contain important targets for drugs in a quantitative fashion. A computational model helps us to understand how the selected elements in the pathways interact with each other, and which component plays the crucial role in systems. We have developed a mathematical model to explore the mechanism for controlling transcription factor activity, which regulates gene expression, by the modulation of calcium signaling triggered during T cell activation. The model simulates the activation and modulation of Ca2+ release-activated Ca2+ (CRAC) channels by mitochondrial dynamics and depletion of endoplasmic reticulum (ER) store, and also includes membrane potential in T-cells. The model simulates the experimental finding that increases in Ca2+ current enhances the activation of transcription factors and the Ca2+ influx through CRAC is also essential for the NFAT and NFκB activation. The model also suggests that plasma membrane Ca2+-ATPase (PMCA) controls a majority of the extrusion of Ca2+ and modulates the activation of CRAC channels. Furthermore, the model simulations explain how the complex interaction of the endoplasmic reticulum, membrane potential, mitochondria, and ion channels such as CRAC channels control T cell activation.

Intercellular Bridge Mediates Ca2+ Signals between Micropatterned Cells via IP3 and Ca2+ Diffusion

Intercellular bridges are plasma continuities formed at the end of the cytokinesis process that facilitate intercellular mass transport between the two daughter cells. However, it remains largely unknown how the intercellular bridge mediates Ca²⁺ communication between postmitotic cells. In this work, we utilize BV-2 microglial cells planted on dumbbell-shaped micropatterned assemblies to resolve spatiotemporal characteristics of Ca²⁺ signal transfer over the intercellular bridges. With the use of such micropatterns, considerably longer and more regular intercellular bridges can be obtained than in conventional cell cultures. The initial Ca²⁺ signal is evoked by mechanical stimulation of one of the daughter cells. A considerable time delay is observed between the arrivals of passive Ca²⁺ diffusion and endogenous Ca²⁺ response in the intercellular-bridge-connected cell, indicating two different pathways of the Ca²⁺ communication. Extracellular Ca²⁺ and the paracrine pathway have practically no effect on the endogenous Ca²⁺ response, demonstrated by application of Ca²⁺-free medium, exogenous ATP, and P2Y₁₃ receptor antagonist. In contrast, the endoplasmic reticulum Ca²⁺-ATPase inhibitor thapsigargin and inositol trisphosphate (IP₃) receptor blocker 2-aminoethyl diphenylborate significantly inhibit the endogenous Ca²⁺ increase, which signifies involvement of IP₃-sensitive calcium store release. Notably, passive Ca²⁺ diffusion into the connected cell can clearly be detected when IP₃-sensitive calcium store release is abolished by 2-aminoethyl diphenylborate. Those observations prove that both passive Ca²⁺ diffusion and IP₃-mediated endogenous Ca²⁺ response contribute to the Ca²⁺ increase in intercellular-bridge-connected cells. Moreover, a simulation model agreed well with the experimental observations.

Cytoplasmic Calcium Buffering: An Integrative Crosstalk

Calcium (Ca²⁺) buffering is part of an integrative crosstalk between different mechanisms and elements involved in the control of free Ca²⁺ ions persistence in the cytoplasm and hence, in the Ca²⁺-dependence of many intracellular processes. Alterations of Ca²⁺ homeostasis and signaling from systemic to subcellular levels also play a pivotal role in the pathogenesis of many diseases.Compared with Ca²⁺ sequestration towards intracellular Ca²⁺ stores, Ca²⁺ buffering is a rapid process occurring in a subsecond scale. Any molecule (or binding site) with the ability to bind Ca²⁺ ions could be considered, at least in principle, as a buffer. However, the term Ca²⁺ buffer is applied only to a small subset of Ca²⁺ binding proteins containing acidic side-chain residues.Ca²⁺ buffering in the cytoplasm mainly relies on mobile and immobile or fixed buffers controlling the diffusion of free Ca²⁺ ions inside the cytosol both temporally and spatially. Mobility of buffers depends on their molecular weight, but other parameters as their concentration, affinity for Ca²⁺ or Ca²⁺ binding and dissociation kinetics next to their diffusional mobility also contribute to make Ca²⁺ signaling one of the most complex signaling activities of the cell.The crosstalk between all the elements involved in the intracellular Ca²⁺ dynamics is a process of extreme complexity due to the diversity of structural and molecular elements involved but permit a highly regulated spatiotemporal control of the signal mediated by Ca²⁺ ions. The basis of modeling tools to study Ca²⁺ dynamics are also presented.

Analytic solutions for calcium ion fertilisation waves on the surface of eggs

The evolution of calcium fertilisation waves on the cortex of amphibian eggs can be described by a nonlinear reaction-diffusion process on the surface of a sphere. Here, we use the nonclassical symmetry technique to find an exact analytic solution that describes the evolution of the calcium concentration. The solutions presented compare well with published experimental results. The analytic solution can be used to give insight into the processes governing the fertilisation wave, such as the flow of calcium ions from the sperm entry point. By finding a spiral solution to an approximate equation linearised near saturation, we also demonstrate how solutions with other properties may be constructed using this technique.

Rethinking calcium profiles around single channels: the exponential and periodic calcium nanodomains

Many fundamental calcium-dependent physiological processes are triggered by high local calcium levels that are established around the sites of calcium entry into the cell (channels). They are dubbed as calcium nanodomains but their exact profiles are still elusive. The concept of calcium nanodomains stems from a linear model of calcium diffusion and is only valid when calcium increases are smaller than the concentration of cytoplasmic buffers. Recent data indicates that much higher calcium levels cause buffer saturation. Therefore, I sought explicit solutions of a nonlinear reaction-diffusion model and found a dichotomous solution. For small fluxes, the steady state calcium profile is quasi-exponential, and when calcium exceeds buffer concentration a spatial periodicity appears. Analytical results are supported by Monte-Carlo simulations. I also imaged 1D- and radial calcium distributions around single α-synuclein channels in cell-free conditions. Measured Ca profiles are consistent with theoretical predictions. I propose that the periodic calcium patterns may well arise under certain conditions and their specific functional role has to be established.

Extension of Rapid Buffering Approximation to Ca2+ Buffers with Two Binding Sites

Fundamental cell processes such as synaptic neurotransmitter release, endocrine hormone secretion, and myocyte contraction are controlled by highly localized calcium (Ca²⁺) signals resulting from brief openings of trans-membrane Ca²⁺ channels. On short temporal and spatial scales, the corresponding local Ca²⁺ nanodomains formed in the vicinity of a single or several open Ca²⁺ channels can be effectively approximated by quasi-stationary solutions. The rapid buffering approximation (RBA) is one of the most powerful of such approximations, and is based on the assumption of instantaneous equilibration of the bimolecular Ca²⁺ buffering reaction, combined with the conservation condition for the total Ca²⁺ and buffer molecule numbers. Previously, RBA has been generalized to an arbitrary arrangement of Ca²⁺ channels on a flat membrane, in the presence of any number of simple Ca²⁺ buffers with one-to-one Ca²⁺ binding stoichiometry. However, many biological buffers have multiple binding sites. For example, buffers and sensors phylogenetically related to calmodulin consist of two Ca²⁺-binding domains (lobes), with each domain binding two Ca²⁺ ions in a cooperative manner. Here we consider an extension of RBA to such buffers with two interdependent Ca²⁺ binding sites. We show that in the presence of such buffers, RBA solution is given by the solution to a cubic equation, analogous to the quadratic equation describing RBA in the case of a simple, one-to-one Ca²⁺ buffer. We examine in detail the dependence of RBA accuracy on buffering parameters, to reveal conditions under which RBA provides sufficient precision.

Research Progress on the Relationship Between Acute Pancreatitis and Calcium Overload in Acinar Cells

Acute pancreatitis is a human disease with multiple causes that leads to autodigestion of the pancreas. There is sufficient evidence to support the key role of sustained increase in cytosolic calcium concentrations in the early pathogenesis of the disease. To clarify the mechanism of maintaining calcium homeostasis in the cell and pathological processes caused by calcium overload would help to research directly targeted therapeutic agents. We will specifically review the following: intracellular calcium homeostasis and regulation, the occurrence of calcium overload in acinar cells, the role of calcium overload in the pathogenesis of AP, the treatment strategy proposed for calcium overload.

Simulation of P2X-mediated calcium signalling in microglia

KEY POINTS

A computational model of P2X channel activation in microglia was developed that includes downfield Ca²⁺ -dependent signalling pathways. This model provides quantitative insights into how diverse signalling pathways in microglia converge to control microglial function.

ABSTRACT

Microglia function is orchestrated through highly coupled signalling pathways that depend on calcium (Ca²⁺ ). In response to extracellular ATP, transient increases in intracellular Ca²⁺ driven through the activation of purinergic receptors, P2X and P2Y, are sufficient to promote cytokine synthesis. Although the steps comprising the pathways bridging purinergic receptor activation with transcriptional responses have been probed in great detail, a quantitative model for how these steps collectively control cytokine production has not been established. Here we developed a minimal computational model that quantitatively links extracellular stimulation of two prominent ionotropic purinergic receptors, P2X4 and P2X7, with the graded production of a gene product, namely the tumour necrosis factor α (TNFα) cytokine. In addition to Ca²⁺ handling mechanisms common to eukaryotic cells, our model includes microglia-specific processes including ATP-dependent P2X4 and P2X7 activation, activation of nuclear factor of activated T-cells (NFAT) transcription factors, and TNFα production. Parameters for this model were optimized to reproduce published data for these processes, where available. With this model, we determined the propensity for TNFα production in microglia, subject to a wide range of ATP exposure amplitudes, frequencies and durations that the cells could encounter in vivo. Furthermore, we have investigated the extent to which modulation of the signal transduction pathways influence TNFα production. Our results suggest that pulsatile stimulation of P2X4 via micromolar ATP may be sufficient to promote TNFα production, whereas high-amplitude ATP exposure is necessary for production via P2X7. Furthermore, under conditions that increase P2X4 expression, for instance, following activation by pathogen-associated molecular factors, P2X4-associated TNFα production is greatly enhanced. Given that Ca²⁺ homeostasis in microglia is profoundly important to its function, this computational model provides a quantitative framework to explore hypotheses pertaining to microglial physiology.

A 3D diffusional-compartmental model of the calcium dynamics in cytosol, sarcoplasmic reticulum and mitochondria of murine skeletal muscle fibers

Variations of free calcium concentration ([Ca2+]) are powerful intracellular signals, controlling contraction as well as metabolism in muscle cells. To fully understand the role of calcium redistribution upon excitation and contraction in skeletal muscle cells, the local [Ca2+] in different compartments needs to be taken into consideration. Fluorescent probes allow the determination of [Ca2+] in the cytosol where myofibrils are embedded, the lumen of the sarcoplasmic reticulum (SR) and the mitochondrial matrix. Previously, models have been developed describing intracellular calcium handling in skeletal and cardiac muscle cells. However, a comprehensive model describing the kinetics of the changes in free calcium concentration in these three compartments is lacking. We designed a new 3D compartmental model of the half sarcomere with radial symmetry, which accounts for diffusion of Ca2+ into the three compartments and simulates its dynamics at rest and at various rates of stimulation in mice skeletal muscle fibers. This model satisfactorily reproduces both the amplitude and time course of the variations of [Ca2+] in the three compartments in mouse fast fibers. As an illustration of the applicability of the model, we investigated the effects of Calsequestrin (CSQ) ablation. CSQ is the main Ca2+ buffer in the SR, localized in close proximity of its calcium release sites and near to the mitochondria. CSQ knock-out mice muscles still preserve a near-normal contractile behavior, but it is unclear whether this is caused by additional SR calcium buffering or a significant contribution of calcium entry from extracellular space, via stored-operated calcium entry (SOCE). The model enabled quantitative assessment of these two scenarios by comparison to measurements of local calcium in the cytosol, the SR and the mitochondria. In conclusion, the model represents a useful tool to investigate the impact of protein ablation and of pharmacological interventions on intracellular calcium dynamics in mice skeletal muscle.

Modeling the influence of co-localized intracellular calcium stores on the secretory response of bovine chromaffin cells

Catecholamines secretion from chromaffin cells is mediated by a Ca²⁺-dependent process in the submembrane space where the exocytotic machinery is located and high-Ca²⁺ microdomains (HCMDs) are formed by the coordinated activity of a functional triad composed of Ca²⁺ channels, endoplasmic reticulum (ER) and mitochondria. It has been observed experimentally that subpopulations of cortical mitochondria and ER associate to secretory sites in bovine chromaffin cells. Here, we study the effect of the geometrical distribution of the co-localized cortical organelles both in the formation of HCMDs in the vicinity of Ca²⁺ channels and on the secretory activity of bovine chromaffin cells in response to a single voltage pulse. Our simulations indicate that co-localized organelles have a dual role in the formation of HCMDs, having, on the one hand, an amplification effect due to the Ca²⁺-induced Ca²⁺-release mechanism from the ER and, on the other, acting as physical barriers to Ca²⁺ diffusion. In addition, our simulations suggest that the increased levels of Ca²⁺ in the microdomain enhances the secretion of the vesicles co-localized to the Ca²⁺ channels. As a whole, our results support the idea that the functional triads formed by Ca²⁺ channels, subplasmalemma ER and mitochondria have a positive effect on the secretion of catecholamines in bovine chromaffin cells.

A biophysically detailed computational model of urinary bladder small DRG neuron soma

Bladder small DRG neurons, which are putative nociceptors pivotal to urinary bladder function, express more than a dozen different ionic membrane mechanisms: ion channels, pumps and exchangers. Small-conductance Ca2+-activated K+ (SKCa) channels which were earlier thought to be gated solely by intracellular Ca2+ concentration ([Ca]i) have recently been shown to exhibit inward rectification with respect to membrane potential. The effect of SKCa inward rectification on the excitability of these neurons is unknown. Furthermore, studies on the role of KCa channels in repetitive firing and their contributions to different types of afterhyperpolarization (AHP) in these neurons are lacking. In order to study these phenomena, we first constructed and validated a biophysically detailed single compartment model of bladder small DRG neuron soma constrained by physiological data. The model includes twenty-two major known membrane mechanisms along with intracellular Ca2+ dynamics comprising Ca2+ diffusion, cytoplasmic buffering, and endoplasmic reticulum (ER) and mitochondrial mechanisms. Using modelling studies, we show that inward rectification of SKCa is an important parameter regulating neuronal repetitive firing and that its absence reduces action potential (AP) firing frequency. We also show that SKCa is more potent in reducing AP spiking than the large-conductance KCa channel (BKCa) in these neurons. Moreover, BKCa was found to contribute to the fast AHP (fAHP) and SKCa to the medium-duration (mAHP) and slow AHP (sAHP). We also report that the slow inactivating A-type K+ channel (slow KA) current in these neurons is composed of 2 components: an initial fast inactivating (time constant ∼ 25-100 ms) and a slow inactivating (time constant ∼ 200-800 ms) current. We discuss the implications of our findings, and how our detailed model can help further our understanding of the role of C-fibre afferents in the physiology of urinary bladder as well as in certain disorders.

A computational model of spatio-temporal cardiac intracellular calcium handling with realistic structure and spatial flux distribution from sarcoplasmic reticulum and t-tubule reconstructions

Intracellular calcium cycling is a vital component of cardiac excitation-contraction coupling. The key structures responsible for controlling calcium dynamics are the cell membrane (comprising the surface sarcolemma and transverse-tubules), the intracellular calcium store (the sarcoplasmic reticulum), and the co-localisation of these two structures to form dyads within which calcium-induced-calcium-release occurs. The organisation of these structures tightly controls intracellular calcium dynamics. In this study, we present a computational model of intracellular calcium cycling in three-dimensions (3-D), which incorporates high resolution reconstructions of these key regulatory structures, attained through imaging of tissue taken from the sheep left ventricle using serial block face scanning electron microscopy. An approach was developed to model the sarcoplasmic reticulum structure at the whole-cell scale, by reducing its full 3-D structure to a 3-D network of one-dimensional strands. The model reproduces intracellular calcium dynamics during control pacing and reveals the high-resolution 3-D spatial structure of calcium gradients and intracellular fluxes in both the cytoplasm and sarcoplasmic reticulum. We also demonstrated the capability of the model to reproduce potentially pro-arrhythmic dynamics under perturbed conditions, pertaining to calcium-transient alternans and spontaneous release events. Comparison with idealised cell models emphasised the importance of structure in determining calcium gradients and controlling the spatial dynamics associated with calcium-transient alternans, wherein the probabilistic nature of dyad activation and recruitment was constrained. The model was further used to highlight the criticality in calcium spark propagation in relation to inter-dyad distances. The model presented provides a powerful tool for future investigation of structure-function relationships underlying physiological and pathophysiological intracellular calcium handling phenomena at the whole-cell. The approach allows for the first time direct integration of high-resolution images of 3-D intracellular structures with models of calcium cycling, presenting the possibility to directly assess the functional impact of structural remodelling at the cellular scale.

The organisation of the membrane and sub-cellular structures of cells in the heart closely controls the coupling between its electrical and mechanical function. Computational models of the cellular calcium handling system, which is responsible for this electro-mechanical coupling, have been developed in recent years to study underlying structure-function relationships. Previous models have been largely idealised in structure; we present a new model which incorporates experimental data describing the high-resolution organisation of the primary structures involved in calcium dynamics. Significantly, the structure of the intracellular calcium store is modelled for the first time. The model is shown to reproduce calcium dynamics in control cells in both normal and abnormal conditions, demonstrating its suitability for future investigation of structure-function relationships. Thus, the model presented provides a powerful tool for the direct integration of experimentally acquired structural data in healthy and diseased cells and assessment of the role of structure in regulating normal and abnormal calcium dynamics.

Padé Approximation of a Stationary Single-Channel Ca2+ Nanodomain

We consider the stationary solution for the Ca²⁺ concentration near a point Ca²⁺ source describing a single-channel Ca²⁺ nanodomain, in the presence of a single mobile buffer with one-to-one Ca²⁺ binding stoichiometry. Previously, a number of Ca²⁺ nanodomains approximations have been developed, for instance the excess buffer approximation (EBA), the rapid buffering approximation (RBA), and the linear approximation (LIN), each valid for appropriate buffering conditions. Apart from providing a simple method of estimating Ca²⁺ and buffer concentrations without resorting to computationally expensive numerical solution of reaction-diffusion equations, such approximations proved useful in revealing the dependence of nanodomain Ca²⁺ distribution on crucial parameters such as buffer mobility and its Ca²⁺ binding properties. In this study, we present a different form of analytic approximation, which is based on matching the short-range Taylor series of the nanodomain concentration with the long-range asymptotic series expressed in inverse powers of distance from channel location. Namely, we use a "dual" Padé rational function approximation to simultaneously match terms in the short- and the long-range series, and we show that this provides an accurate approximation to the nanodomain Ca²⁺ and buffer concentrations. We compare this approximation with the previously obtained approximations and show that it yields a better estimate of the free buffer concentration for a wide range of buffering conditions. The drawback of our method is that it has a complex algebraic form for any order higher than the lowest bilinear order, and cannot be readily extended to multiple Ca²⁺ channels. However, it may be possible to extend the Padé method to estimate Ca²⁺ nanodomains in the presence of cooperative Ca²⁺ buffers with two Ca²⁺ binding sites, the case that existing methods do not address.

A model of auto immune response

Background

In this work, we develop a theoretical model of an auto immune response. This is based on modifications of standard second messenger trigger models using both signalling pathways and diffusion and a macro level dynamic systems approximation to the response of a triggering agent such as a virus, bacteria or environmental toxin.

Results

We show that there, in general, will be self damage effects whenever the triggering agent’s effect on the host can be separated into two distinct classes of cell populations.

In each population, the trigger acts differently and this behavior is mediated by the nonlinear interactions between two signalling agents.

Conclusion

If these interactions satisfy certain critical assumptions this will lead to collateral damage. If the initial triggering agent’s action involves any critical host cell population whose loss can lead to serious host health issues, then there is a much increased probability of host death.

Our model also shows that if the nonlinear interaction assumptions are satisfied, there is a reasonable expectation of oscillatory behavior in host health; i.e. periods of remission.