Range of messenger action of calcium ion and inositol 1,4,5-trisphosphate

Signalling switches maintain intercellular communication in the vascular endothelium

BACKGROUND AND PURPOSE

The single layer of cells lining all blood vessels, the endothelium, is a sophisticated signal co-ordination centre that controls a wide range of vascular functions including the regulation of blood pressure and blood flow. To co-ordinate activities, communication among cells is required for tissue level responses to emerge. While a significant form of communication occurs by the propagation of signals between cells, the mechanism of propagation in the intact endothelium is unresolved.

EXPERIMENTAL APPROACH

Precision signal generation and targeted cellular manipulation was used in conjunction with high spatiotemporal mesoscale Ca2+ imaging in the endothelium of intact blood vessels.

KEY RESULTS

Multiple mechanisms maintain communication so that Ca2+ wave propagation occurs irrespective of the status of connectivity among cells. Between adjoining cells, regenerative IP3-induced IP3 production transmits Ca2+ signals and explains the propagated vasodilation that underlies the increased blood flow accompanying tissue activity. The inositide is itself sufficient to evoke regenerative phospholipase C-dependent Ca2+ waves across coupled cells. None of gap junctions, Ca2+ diffusion or the release of extracellular messengers is required to support this type of intercellular Ca2+ signalling. In contrast, when discontinuities exist between cells, ATP released as a diffusible extracellular messenger transmits Ca2+ signals across the discontinuity and drives propagated vasodilation.

CONCLUSION AND IMPLICATIONS

These results show that signalling switches underlie endothelial cell-to-cell signal transmission and reveal how communication is maintained in the face of endothelial damage. The findings provide a new framework for understanding wave propagation and cell signalling in the endothelium.

Toward understanding the cellular control of vertebrate mineralization: The potential role of mitochondria

This review examines the possible role of mitochondria in maintaining calcium and phosphate ion homeostasis and participating in the mineralization of bone, cartilage and other vertebrate hard tissues. The paper builds on the known structural features of mitochondria and the documented observations in these tissues that the organelles contain calcium phosphate granules. Such deposits in mitochondria putatively form to buffer excessively high cytosolic calcium ion concentrations and prevent metabolic deficits and even cell death. While mitochondria protect cytosolic enzyme systems through this buffering capacity, the accumulation of calcium ions by mitochondria promotes the activity of enzymes of the tricarboxylic acid (TCA/Krebs) cycle, increases oxidative phosphorylation and ATP synthesis, and leads to changes in intramitochondrial pH. These pH alterations influence ion solubility and possibly the transitions and composition in the mineral phase structure of the granules. Based on these considerations, mitochondria are proposed to support the mineralization process by providing a mobile store of calcium and phosphate ions, in smaller cluster or larger granule form, while maintaining critical cellular activities. The rise in the mitochondrial calcium level also increases the generation of citrate and other TCA cycle intermediates that contribute to cell function and the development of extracellular mineral. This paper suggests that another key role of the mitochondrion, along with the effects just noted, is to supply phosphate ions, derived from the breakdown of ATP, to endolysosomes and autophagic vesicles originating in the endoplasmic reticulum and Golgi and at the plasma membrane. These many separate but interdependent mitochondrial functions emphasize the critical importance of this organelle in the cellular control of vertebrate mineralization.

Emergence of large-scale cell death through ferroptotic trigger waves

Large-scale cell death is commonly observed during organismal development and in human pathologies1-5. These cell death events extend over great distances to eliminate large populations of cells, raising the question of how cell death can be coordinated in space and time. One mechanism that enables long-range signal transmission is trigger waves6, but how this mechanism might be used for death events in cell populations remains unclear. Here we demonstrate that ferroptosis, an iron- and lipid-peroxidation-dependent form of cell death, can propagate across human cells over long distances (≥5 mm) at constant speeds (around 5.5 μm min-1) through trigger waves of reactive oxygen species (ROS). Chemical and genetic perturbations indicate a primary role of ROS feedback loops (Fenton reaction, NADPH oxidase signalling and glutathione synthesis) in controlling the progression of ferroptotic trigger waves. We show that introducing ferroptotic stress through suppression of cystine uptake activates these ROS feedback loops, converting cellular redox systems from being monostable to being bistable and thereby priming cell populations to become bistable media over which ROS propagate. Furthermore, we demonstrate that ferroptosis and its propagation accompany the massive, yet spatially restricted, cell death events during muscle remodelling of the embryonic avian limb, substantiating its use as a tissue-sculpting strategy during embryogenesis. Our findings highlight the role of ferroptosis in coordinating global cell death events, providing a paradigm for investigating large-scale cell death in embryonic development and human pathologies.

Differential behaviors of calcium-induced calcium release in one dimensional dendrite by Nernst-Planck equation, cable model and pure diffusion model

The source and dynamics of calcium is the key factor that regulates dendritic integration. Apart from the voltage-gated and ligand-gated calcium influx, an important source of calcium is from inner store of endoplasmic reticulum with a regenerative process of calcium-induced calcium release (CICR). To trigger this process, inositol 1,4,5-trisphosphate (IP3) and calcium are needed to satisfy certain requirements. The aim of our paper is to investigate how the CICR depends on the dynamics of membrane potential. We utilize one dimensional dendritic model to calculate membrane potential by Nernst-Planck Equation (NPE) and cable model and Pure Diffusion (PD) model, computational simulations are carried out to inject the calcium influx by synaptic stimulation and to predict subsequent CICR and calcium wave propagation. Our results demonstrate that CICR initiation and calcium wave propagation have much difference between electro-diffusion process of NPE and cable model. We find that cable model has lower threshold of IP3 stimulation to trigger CICR but is more difficult for calcium propagation than NPE, PD model requires even higher threshold of IP3 to initiate CICR process and calcium duration is shorter than NPE; the regenerative calcium wave propagates with faster speed in NPE than that in cable model and in PD model. Our work addresses the important role of electro-diffusion dynamics of charged ions in regulating CICR process in dendritic structure; and provides theoretical predictions for neurological process which requires sustaining calcium for downstream signaling processes.

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.

Neuroprotective gap-junction-mediated bystander transformations in the adult zebrafish spinal cord after injury

The adult zebrafish spinal cord displays an impressive innate ability to regenerate after traumatic insults, yet the underlying adaptive cellular mechanisms remain elusive. Here, we show that while the cellular and tissue responses after injury are largely conserved among vertebrates, the large-size fast spinal zebrafish motoneurons are remarkably resilient by remaining viable and functional. We also reveal the dynamic changes in motoneuron glutamatergic input, excitability, and calcium signaling, and we underscore the critical role of calretinin (CR) in binding and buffering the intracellular calcium after injury. Importantly, we demonstrate the presence and the dynamics of a neuron-to-neuron bystander neuroprotective biochemical cooperation mediated through gap junction channels. Our findings support a model in which the intimate and dynamic interplay between glutamate signaling, calcium buffering, gap junction channels, and intercellular cooperation upholds cell survival and promotes the initiation of regeneration.

Cellular communication among smooth muscle cells: The role of membrane potential via connexins

Communication via action potentials among neurons has been extensively studied. However, effective communication without action potentials is ubiquitous in biological systems, yet it has received much less attention in comparison. Multi-cellular communication among smooth muscles is crucial for regulating blood flow, for example. Understanding the mechanism of this non-action potential communication is critical in many cases, like synchronization of cellular activity, under normal and pathological conditions. In this paper, we employ a multi-scale asymptotic method to derive a macroscopic homogenized bidomain model from the microscopic electro-neutral (EN) model. This is achieved by considering different diffusion coefficients and incorporating nonlinear interface conditions. Subsequently, the homogenized macroscopic model is used to investigate communication in multi-cellular tissues. Our computational simulations reveal that the membrane potential of syncytia, formed by interconnected cells via connexins, plays a crucial role in propagating oscillations from one region to another, providing an effective means for fast cellular communication. Statement of Significance: In this study, we investigated cellular communication and ion transport in vascular smooth muscle cells, shedding light on their mechanisms under normal and abnormal conditions. Our research highlights the potential of mathematical models in understanding complex biological systems. We developed effective macroscale electro-neutral bi-domain ion transport models and examined their behavior in response to different stimuli. Our findings revealed the crucial role of connexinmediated membrane potential changes and demonstrated the effectiveness of cellular communication through syncytium membranes. Despite some limitations, our study provides valuable insights into these processes and emphasizes the importance of mathematical modeling in unraveling the complexities of cellular communication and ion transport.

Nanojunctions: Specificity of Ca2+ signaling requires nano-scale architecture of intracellular membrane contact sites

Spatio-temporal definition of Ca2+ signals involves the assembly of signaling complexes within the nano-architecture of contact sites between the sarco/endoplasmic reticulum (SR/ER) and the plasma membrane (PM). While the requirement of precise spatial assembly and positioning of the junctional signaling elements is well documented, the role of the nano-scale membrane architecture itself, as an ion-reflecting confinement of the signalling unit, remains as yet elusive. Utilizing the Na+/Ca2+ Exchanger-1 / SR/ER Ca2+ ATPase-2-mediated ER Ca2+ refilling process as a junctional signalling paradigm, we provide here the first evidence for an indispensable cellular function of the junctional membrane architecture. Our stochastic modeling approach demonstrates that junctional ER Ca2+ refilling operates exclusively at nano-scale membrane spacing, with a strong inverse relationship between junctional width and signaling efficiency. Our model predicts a breakdown of junctional Ca2+ signaling with loss of reflecting membrane confinement. In addition we consider interactions between Ca2+ and the phospholipid membrane surface, which may support interfacial Ca2+ transport and promote receptor targeting. Alterations in the molecular and nano-scale membrane organization at organelle-PM contacts are suggested as a new concept in pathophysiology.

Amyloid β-Oligomers Inhibit the Nuclear Ca2+ Signals and the Neuroprotective Gene Expression Induced by Gabazine in Hippocampal Neurons

Hippocampal neuronal activity generates dendritic and somatic Ca2+ signals, which, depending on stimulus intensity, rapidly propagate to the nucleus and induce the expression of transcription factors and genes with crucial roles in cognitive functions. Soluble amyloid-beta oligomers (AβOs), the main synaptotoxins engaged in the pathogenesis of Alzheimer's disease, generate aberrant Ca2+ signals in primary hippocampal neurons, increase their oxidative tone and disrupt structural plasticity. Here, we explored the effects of sub-lethal AβOs concentrations on activity-generated nuclear Ca2+ signals and on the Ca2+-dependent expression of neuroprotective genes. To induce neuronal activity, neuron-enriched primary hippocampal cultures were treated with the GABAA receptor blocker gabazine (GBZ), and nuclear Ca2+ signals were measured in AβOs-treated or control neurons transfected with a genetically encoded nuclear Ca2+ sensor. Incubation (6 h) with AβOs significantly reduced the nuclear Ca2+ signals and the enhanced phosphorylation of cyclic AMP response element-binding protein (CREB) induced by GBZ. Likewise, incubation (6 h) with AβOs significantly reduced the GBZ-induced increases in the mRNA levels of neuronal Per-Arnt-Sim domain protein 4 (Npas4), brain-derived neurotrophic factor (BDNF), ryanodine receptor type-2 (RyR2), and the antioxidant enzyme NADPH-quinone oxidoreductase (Nqo1). Based on these findings we propose that AβOs, by inhibiting the generation of activity-induced nuclear Ca2+ signals, disrupt key neuroprotective gene expression pathways required for hippocampal-dependent learning and memory processes.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.

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.

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.

Indicator-dependent differences in detection of local intracellular Ca2+ release events evoked by voltage-gated Ca2+ entry in pancreatic β-cells

Genetically encoded Ca²⁺ indicators have become widely used in cell signalling studies as they offer advantages over cell-loaded dye indicators in enabling specific cellular or subcellular targeting. Comparing responses from dye and protein-based indicators may provide information about indicator properties and cell physiology, but side-by-side recordings in cells are scarce. In this study, we compared cytoplasmic Ca²⁺ concentration ([Ca²⁺]_i) changes in insulin-secreting β-cells recorded with commonly used dyes and indicators based on circularly permuted fluorescent proteins. Total internal reflection fluorescence (TIRF) imaging of K⁺ depolarization-triggered submembrane [Ca²⁺]_i increases showed that the dyes Fluo-4 and Fluo-5F mainly reported stable [Ca²⁺]_i elevations, whereas the proteins R-GECO1 and GCaMP5G more often reported distinct [Ca²⁺]_i spikes from an elevated level. [Ca²⁺]_i spiking occurred also in glucose-stimulated cells. The spikes reflected Ca²⁺ release from the endoplasmic reticulum, triggered by autocrine activation of purinergic receptors after exocytotic release of ATP and/or ADP, and the spikes were consequently prevented by SERCA inhibition or P2Y₁-receptor antagonism. Widefield imaging, which monitors the entire cytoplasm, increased the spike detection by the Ca²⁺ dyes. The indicator-dependent response patterns were unrelated to Ca²⁺ binding affinity, buffering and mobility, and probably reflects the much slower dissociation kinetics of protein compared to dye indicators. Ca²⁺ dyes thus report signalling within the submembrane space excited by TIRF illumination, whereas the protein indicators also catch Ca²⁺ events originating outside this volume. The study highlights that voltage-dependent Ca²⁺ entry in β-cells is tightly linked to local intracellular Ca²⁺ release mediated via an autocrine route that may be more important than previously reported direct Ca²⁺ effects on phospholipase C or on intracellular channels mediating calcium-induced calcium release.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca2+ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca2+ -handling proteins such as the L-type Ca2+ channel and Na+ -Ca2+ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca2+ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca2+ -handling and Ca2+ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca2+ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca2+ release. We found that varying TATS density and thus the localization of key Ca2+ -handling proteins has profound effects on Ca2+ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca2+ concentration through decreased Ca2+ extrusion. This elevated Ca2+ increases RyR open probability causing spontaneous Ca2+ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca2+ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca2+ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca2+ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca2+ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca2+ handling governed by the TATS. Simulated TATS loss causes diastolic Ca2+ and voltage instabilities through reduced Na+ -Ca2+ exchanger-mediated Ca2+ removal, cleft Ca2+ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca2+ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca2+ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca2+ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

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.

Cav2.2 Channels Sustain Vesicle Recruitment at a Mature Glutamatergic Synapse

The composition of voltage-gated Ca2+ channel (Cav) subtypes that gate action potential (AP)-evoked release changes during the development of mammalian CNS synapses. Cav2.2 and Cav2.3 lose their function in gating-evoked release during postnatal synapse maturation. In mature boutons, Cav2.1 currents provide the almost exclusive trigger for evoked release, and Cav2.3 currents are required for the induction of presynaptic long-term potentiation. However, the functional significance of Cav2.2 remained elusive in mature boutons, although they remain present at active zones and continue contributing significantly to presynaptic Ca2+ influx. Here, we addressed the functional significance of Cav2.2 and Cav2.3 at mature parallel-fiber (PF) to Purkinje neuron synapses of mice of either sex. These synapses are known to exhibit the corresponding developmental Cav subtype changes in gating release. We addressed two hypotheses, namely that Cav2.2 and Cav2.3 are involved in triggering spontaneous glutamate release and that they are engaged in vesicle recruitment during repetitive evoked release. We found that spontaneous miniature release is Ca2+ dependent. However, experiments with Cav subtype-specific blockers excluded the spontaneous opening of Cavs as the Ca2+ source for spontaneous glutamate release. Thus, neither Cav2.2 nor Cav2.3 controls spontaneous release from PF boutons. Furthermore, vesicle recruitment during brief bursts of APs was also independent of Ca2+ influx through Cav2.2 and Cav2.3. However, Cav2.2, but not Cav2.3, currents significantly boosted vesicle recruitment during sustained high-frequency synaptic transmission. Thus, in mature PF boutons Cav2.2 channels are specifically required to sustain synaptic transmission during prolonged neuronal activity.SIGNIFICANCE STATEMENT At young CNS synapses, action potential-evoked release is gated via three subtypes of voltage-gated Ca2+ channels: Cav2.1, Cav2.2, and Cav2.3. During postnatal maturation, Cav2.2 and Cav2.3 lose their function in gating evoked release, such that at mature synapses Cav2.1 provides the almost exclusive source for triggering evoked release. Cav2.3 currents are required for the induction of presynaptic long-term potentiation. However, the function of the still abundant Cav2.2 in mature boutons remained largely elusive. Here, we studied mature cerebellar parallel-fiber synapses and found that Cav2.2 does not control spontaneous release. However, Ca2+ influx through Cav2.2 significantly boosted vesicle recruitment during trains of action potentials. Thus, Cav2.2 in mature parallel-fiber boutons participate in sustaining synaptic transmission during prolonged activity.

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.

Non-inositol 1,4,5-trisphosphate (IP3) receptor IP3-binding proteins

Conventionally, myo-D-inositol 1, 4,5-trisphosphate (IP₃) is thought to exert its second messenger effects through the gating of IP₃R Ca²⁺ release channels, located in Ca²⁺-storage organelles like the endoplasmic reticulum. However, there is considerable indirect evidence to support the concept that IP₃ might interact with other, non-IP₃R proteins within cells. To explore this possibility further, the Protein Data Bank was searched using the term "IP3". This resulted in the retrieval of 203 protein structures, the majority of which were members of the IP₃R/ryanodine receptor superfamily of channels. Only 49 of these structures were complexed with IP₃. These were inspected for their ability to interact with the carbon-1 phosphate of IP₃, since this is the least accessible phosphate group of its precursor, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂). This reduced the number of structures retrieved to 35, of which 9 were IP₃Rs. The remaining 26 structures represent a diverse range of proteins, including inositol-lipid metabolizing enzymes, signal transducers, PH domain containing proteins, cytoskeletal anchor proteins, the TRPV4 ion channel, a retroviral Gag protein and fibroblast growth factor 2. Such proteins may impact on IP₃ signalling and its effects on cell-biology. This represents an area open for exploration in the field of IP₃ signalling.

Regulation of cardiac function by cAMP nanodomains

Cyclic adenosine monophosphate (cAMP) is a diffusible intracellular second messenger that plays a key role in the regulation of cardiac function. In response to the release of catecholamines from sympathetic terminals, cAMP modulates heart rate and the strength of contraction and ease of relaxation of each heartbeat. At the same time, cAMP is involved in the response to a multitude of other hormones and neurotransmitters. A sophisticated network of regulatory mechanisms controls the temporal and spatial propagation of cAMP, resulting in the generation of signaling nanodomains that enable the second messenger to match each extracellular stimulus with the appropriate cellular response. Multiple proteins contribute to this spatiotemporal regulation, including the cAMP-hydrolyzing phosphodiesterases (PDEs). By breaking down cAMP to a different extent at different locations, these enzymes generate subcellular cAMP gradients. As a result, only a subset of the downstream effectors is activated and a specific response is executed. Dysregulation of cAMP compartmentalization has been observed in cardiovascular diseases, highlighting the importance of appropriate control of local cAMP signaling. Current research is unveiling the molecular organization underpinning cAMP compartmentalization, providing original insight into the physiology of cardiac myocytes and the alteration associated with disease, with the potential to uncover novel therapeutic targets. Here, we present an overview of the mechanisms that are currently understood to be involved in generating cAMP nanodomains and we highlight the questions that remain to be answered.

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.

Computational investigation of IP3 diffusion

Inositol 1,4,5-trisphosphate (IP₃) plays a key role in calcium signaling. After stimulation, it diffuses from the plasma membrane where it is produced to the endoplasmic reticulum where its receptors are localized. Based on in vitro measurements, IP₃ was long thought to be a global messenger characterized by a diffusion coefficient of ~ 280 μm²s^-1. However, in vivo observations revealed that this value does not match with the timing of localized Ca²⁺ increases induced by the confined release of a non-metabolizable IP₃ analog. A theoretical analysis of these data concluded that in intact cells diffusion of IP₃ is strongly hindered, leading to a 30-fold reduction of the diffusion coefficient. Here, we performed a new computational analysis of the same observations using a stochastic model of Ca²⁺ puffs. Our simulations concluded that the value of the effective IP₃ diffusion coefficient is close to 100 μm²s^-1. Such moderate reduction with respect to in vitro estimations quantitatively agrees with a buffering effect by non-fully bound inactive IP₃ receptors. The model also reveals that IP₃ spreading is not much affected by the endoplasmic reticulum, which represents an obstacle to the free displacement of molecules, but can be significantly increased in cells displaying elongated, 1-dimensional like geometries.

Localisation of Intracellular Signals and Responses during Phagocytosis

Phagocytosis is one of the most polarised of all cellular activities. Both the stimulus (the target for phagocytosis) and the response (its internalisation) are focussed at just one part of the cell. At the locus, and this locus alone, pseudopodia form a phagocytic cup around the particle, the cytoskeleton is rearranged, the plasma membrane is reorganised, and a new internal organelle, the phagosome, is formed. The effect of signals from the stimulus must, thus, both be complex and yet be restricted in space and time to enable an effective focussed response. While many aspects of phagocytosis are being uncovered, the mechanism for the restriction of signalling or the effects of signalling remains obscure. In this review, the details of the problem of restricting chemical intracellular signalling are presented, with a focus on diffusion into the cytosol and of signalling lipids along the plasma membrane. The possible ways in which simple diffusion is overcome so that the restriction of signalling and effective phagocytosis can be achieved are discussed in the light of recent advances in imaging, biophysics, and cell biochemistry which together are providing new insights into this area.

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.

Diffusive exit rates through pores in membrane-enclosed structures

The function of many membrane-enclosed intracellular structures relies on release of diffusing particles that exit through narrow pores or channels in the membrane. The rate of release varies with pore size, density, and length of the channel. We propose a simple approximate model, validated with stochastic simulations, for estimating the effective release rate from cylinders, and other simple-shaped domains, as a function of channel parameters. The results demonstrate that, for very small pores, a low density of channels scattered over the boundary is sufficient to achieve substantial rates of particle release. Furthermore, we show that increasing the length of passive channels will both reduce release rates and lead to a less steep dependence on channel density. Our results are compared to previously-measured local calcium release rates from tubules of the endoplasmic reticulum, providing an estimate of the relevant channel density responsible for the observed calcium efflux.

Involvement of Ca2+ in Signaling Mechanisms Mediating Muscarinic Inhibition of M Currents in Sympathetic Neurons

Acetylcholine can excite neurons by suppressing M-type (KCNQ) potassium channels. This effect is mediated by M1 muscarinic receptors coupled to the Gq protein. Although PIP2 depletion and PKC activation have been strongly suggested to contribute to muscarinic inhibition of M currents (IM), direct evidence is lacking. We investigated the mechanism involved in muscarinic inhibition of IM with Ca2+ measurement and electrophysiological studies in both neuronal (rat sympathetic neurons) and heterologous (HEK cells expressing KCNQ2/KCNQ3) preparations. We found that muscarinic inhibition of IM was not blocked either by PIP2 or by calphostin C, a PKC inhibitor. We then examined whether muscarinic inhibition of IM uses multiple signaling pathways by blocking both PIP2 depletion and PKC activation. This maneuver, however, did not block muscarinic inhibition of IM. Additionally, muscarinic inhibition of IM was not prevented either by sequestering of G-protein βγ subunits from Gα-transducin or anti-Gβγ antibody or by preventing intracellular trafficking of channel proteins with blebbistatin, a class-II myosin inhibitor. Finally, we re-examined the role of Ca2+ signals in muscarinic inhibition of IM. Ca2+ measurements showed that muscarinic stimulation increased intracellular Ca2+ and was comparable to the Ca2+ mobilizing effect of bradykinin. Accordingly, 20-mM of BAPTA significantly suppressed muscarinic inhibition of IM. In contrast, muscarinic inhibition of IM was completely insensitive to 20-mM EGTA. Taken together, these data suggest a role of Ca2+ signaling in muscarinic modulation of IM. The differential effects of EGTA and BAPTA imply that Ca2+ microdomains or spatially local Ca2+ signals contribute to inhibition of IM.

Suspended Collagen Hydrogels to Replicate Human Colonic Epithelial Cell Interactions with Immune Cells

The human colon plays a critical role in fluid and salt absorption and harbors the largest immune compartment. There is a widespread need for in vitro models of human colon physiology with its innate immune system. A method is described to produce a cassette with a network of struts supporting a suspended, non-chemically cross-linked collagen hydrogel scaffold compatible with the co-culture of primary gastrointestinal epithelium and migratory inflammatory cells. The epithelial monolayer cultured on the suspended collagen possesses a population of polarized and differentiated cells similar to that present in vivo. This epithelial layer displays proper barrier function with a transepithelial electrical resistance (TEER) ≥ 1,500 Ω cm2 and an apparent permeability ≤10-5 cm2 s-1 . Immune cells plated on the basal face of the scaffold transmigrated over a period of 24 h to the epithelial layer in response to epithelial production of IL-8 induced by luminal stimulation of Clostridium difficile Toxin A. These studies demonstrate that this in vitro platform possesses a functional primary colonic epithelial layer with an immune cell compartment capable of recruitment in response to pro-inflammatory cues coming from the epithelium.

Computational simulations of the effects of gravity on lymphatic transport

Physical forces, including mechanical stretch, fluid pressure, and shear forces alter lymphatic vessel contractions and lymph flow. Gravitational forces can affect these forces, resulting in altered lymphatic transport, but the mechanisms involved have not been studied in detail. Here, we combine a lattice Boltzmann-based fluid dynamics computational model with known lymphatic mechanobiological mechanisms to investigate the movement of fluid through a lymphatic vessel under the effects of gravity that may either oppose or assist flow. Regularly spaced, mechanical bi-leaflet valves in the vessel enforce net positive flow as the vessel walls contract autonomously in response to calcium and nitric oxide (NO) levels regulated by vessel stretch and shear stress levels. We find that large gravitational forces opposing flow can stall the contractions, leading to no net flow, but transient mechanical perturbations can re-establish pumping. In the case of gravity strongly assisting flow, the contractions also cease due to high shear stress and NO production, which dilates the vessel to allow gravity-driven flow. In the intermediate range of oppositional gravity forces, the vessel actively contracts to offset nominal gravity levels or to modestly assist the favorable hydrostatic pressure gradients.

Simulation of Calcium Dynamics in Realistic Three-Dimensional Domains

The cytosolic concentration of free calcium ions ([Ca2+]) is an important intracellular messenger in most cell types, and the spatial distribution of [Ca2+] is often critical. In a salivary gland acinar cell, a polarised epithelial cell, whose principal function is to transport water and thus secrete saliva, [Ca2+] controls the secretion of primary saliva, but increases in [Ca2+] are localised to the apical regions of the cell. Hence, any quantitative explanation of how [Ca2+] controls saliva secretion must take into careful account the spatial distribution of the various Ca2+ sources, Ca2+ sinks, and Ca2+-sensitive ion channels. Based on optical slices, we have previously constructed anatomically accurate three-dimensional models of seven salivary gland acinar cells, and thus shown that a model in which Ca2+ responses are confined to the apical regions of the cell is sufficient to provide a quantitative and predictive explanation of primary saliva secretion. However, reconstruction of such anatomically accurate cells is extremely time consuming and inefficient. Here, we present an alternative, mostly automated method of constructing three-dimensional cells that are approximately anatomically accurate and show that the new construction preserves the quantitative accuracy of the model.

Ca2+ signalling system initiated by endoplasmic reticulum stress stimulates PERK activation

The accumulation of unfolded proteins within the Endoplasmic Reticulum (ER) activates a signal transduction pathway termed the unfolded protein response (UPR), which attempts to restore ER homoeostasis. If this cannot be done, UPR signalling ultimately induces apoptosis. Ca²⁺ depletion in the ER is a potent inducer of ER stress. Despite the ubiquity of Ca²⁺ as an intracellular messenger, the precise mechanism(s) by which Ca²⁺ release affects the UPR remains unknown. Tethering a genetically encoded Ca²⁺ indicator (GCamP6) to the ER membrane revealed novel Ca²⁺ signalling events initiated by Ca²⁺ microdomains in human astrocytes under ER stress, induced by tunicamycin (Tm), an N-glycosylation inhibitor, as well as in a cell model deficient in all three inositol triphosphate receptor isoforms. Pharmacological and molecular studies indicate that these local events are mediated by translocons and that the Ca²⁺ microdomains impact (PKR)-like-ER kinase (PERK), an UPR sensor, activation. These findings reveal the existence of a Ca²⁺ signal mechanism by which stressor-mediated Ca²⁺ release regulates ER stress.

Cell-To-Cell Communication in the Resistance Vasculature

The arterial vasculature can be divided into large conduit arteries, intermediate contractile arteries, resistance arteries, arterioles, and capillaries. Resistance arteries and arterioles primarily function to control systemic blood pressure. The resistance arteries are composed of a layer of endothelial cells oriented parallel to the direction of blood flow, which are separated by a matrix layer termed the internal elastic lamina from several layers of smooth muscle cells oriented perpendicular to the direction of blood flow. Cells within the vessel walls communicate in a homocellular and heterocellular fashion to govern luminal diameter, arterial resistance, and blood pressure. At rest, potassium currents govern the basal state of endothelial and smooth muscle cells. Multiple stimuli can elicit rises in intracellular calcium levels in either endothelial cells or smooth muscle cells, sourced from intracellular stores such as the endoplasmic reticulum or the extracellular space. In general, activation of endothelial cells results in the production of a vasodilatory signal, usually in the form of nitric oxide or endothelial-derived hyperpolarization. Conversely, activation of smooth muscle cells results in a vasoconstriction response through smooth muscle cell contraction. © 2022 American Physiological Society. Compr Physiol 12: 1-35, 2022.

Temporal and spatial summation of laser heat stimuli in cultured nociceptive neurons of the rat

We studied the efficacy of a near-infrared laser (1475 nm) to activate rat dorsal root ganglion (DRG) neurons with short punctate radiant heat pulses (55 µm diameter) and investigated temporal and spatial summation properties for the transduction process for noxious heat at a subcellular level. Strength-duration curves (10–80 ms range) indicated a minimum power of 30.2mW for the induction of laser-induced calcium transients and a chronaxia of 13.9 ms. However, threshold energy increased with increasing stimulus duration suggesting substantial radial cooling of the laser spot. Increasing stimulus duration demonstrated suprathreshold intensity coding of calcium transients with less than linear gains (Stevens exponents 0.29/35mW, 0.38/60mW, 0.46/70mW). The competitive TRPV1 antagonist capsazepine blocked responses to short near-threshold stimuli and significantly reduced responses to longer duration suprathreshold heat. Heating 1/3 of the soma of a neuron was sufficient to induce calcium transients significantly above baseline (p < 0.05), but maximum amplitude was only achieved by centering the laser over the entire neuron. Heat-induced calcium increase was highest in heated cell parts but rapidly reached unstimulated areas reminiscent of spreading depolarization and opening of voltage-gated calcium channels. Full intracellular equilibrium took about 3 s, consistent with a diffusion process. In summary, we investigated transduction mechanisms for noxious laser heat pulses in native sensory neurons at milliseconds temporal and subcellular spatial resolution and characterized strength duration properties, intensity coding, and spatial summation within single neurons. Thermal excitation of parts of a nociceptor spread via both membrane depolarization and intracellular calcium diffusion.

Multiscale modeling of presynaptic dynamics from molecular to mesoscale

Chemical synapses exhibit a diverse array of internal mechanisms that affect the dynamics of transmission efficacy. Many of these processes, such as release of neurotransmitter and vesicle recycling, depend strongly on activity-dependent influx and accumulation of Ca²⁺. To model how each of these processes may affect the processing of information in neural circuits, and how their dysfunction may lead to disease states, requires a computationally efficient modelling framework, capable of generating accurate phenomenology without incurring a heavy computational cost per synapse. Constructing a phenomenologically realistic model requires the precise characterization of the timing and probability of neurotransmitter release. Difficulties arise in that functional forms of instantaneous release rate can be difficult to extract from noisy data without running many thousands of trials, and in biophysical synapses, facilitation of per-vesicle release probability is confounded by depletion. To overcome this, we obtained traces of free Ca²⁺ concentration in response to various action potential stimulus trains from a molecular MCell model of a hippocampal Schaffer collateral axon. Ca²⁺ sensors were placed at varying distance from a voltage-dependent calcium channel (VDCC) cluster, and Ca²⁺ was buffered by calbindin. Then, using the calcium traces to drive deterministic state vector models of synaptotagmin 1 and 7 (Syt-1/7), which respectively mediate synchronous and asynchronous release in excitatory hippocampal synapses, we obtained high-resolution profiles of instantaneous release rate, to which we applied functional fits. Synchronous vesicle release occurred predominantly within half a micron of the source of spike-evoked Ca²⁺ influx, while asynchronous release occurred more consistently at all distances. Both fast and slow mechanisms exhibited multi-exponential release rate curves, whose magnitudes decayed exponentially with distance from the Ca²⁺ source. Profile parameters facilitate on different time scales according to a single, general facilitation function. These functional descriptions lay the groundwork for efficient mesoscale modelling of vesicular release dynamics.

Most information transmission between neurons in the brain occurs via release of neurotransmitter from synaptic vesicles. In response to a presynaptic spike, calcium influx at the active zone of a synapse can trigger the release of neurotransmitter with a certain probability. These stochastic release events may occur immediately after a spike or with some delay. As calcium accumulates from one spike to the next, the probability of release may increase (facilitate) for subsequent spikes. This process, known as short-term plasticity, transforms the spiking code to a release code, underlying much of the brain’s information processing. In this paper, we use an accurate, detailed model of presynaptic molecular physiology to characterize these processes at high precision in response to various spike trains. We then apply model reduction to the results to obtain a phenomenological model of release timing, probability, and facilitation, which can perform as accurately as the molecular model but with far less computational cost. This mesoscale model of spike-evoked release and facilitation helps to bridge the gap between microscale molecular dynamics and macroscale information processing in neural circuits. It can thus benefit large scale modelling of neural circuits, biologically inspired machine learning models, and the design of neuromorphic chips.

Influence of T-Bar on Calcium Concentration Impacting Release Probability

The relation of form and function, namely the impact of the synaptic anatomy on calcium dynamics in the presynaptic bouton, is a major challenge of present (computational) neuroscience at a cellular level. The Drosophila larval neuromuscular junction (NMJ) is a simple model system, which allows studying basic effects in a rather simple way. This synapse harbors several special structures. In particular, in opposite to standard vertebrate synapses, the presynaptic boutons are rather large, and they have several presynaptic zones. In these zones, different types of anatomical structures are present. Some of the zones bear a so-called T-bar, a particular anatomical structure. The geometric form of the T-bar resembles the shape of the letter “T” or a table with one leg. When an action potential arises, calcium influx is triggered. The probability of vesicle docking and neurotransmitter release is superlinearly proportional to the concentration of calcium close to the vesicular release site. It is tempting to assume that the T-bar causes some sort of calcium accumulation and hence triggers a higher release probability and thus enhances neurotransmitter exocytosis. In order to study this influence in a quantitative manner, we constructed a typical T-bar geometry and compared the calcium concentration close to the active zones (AZs). We compared the case of synapses with and without T-bars. Indeed, we found a substantial influence of the T-bar structure on the presynaptic calcium concentrations close to the AZs, indicating that this anatomical structure increases vesicle release probability. Therefore, our study reveals how the T-bar zone implies a strong relation between form and function. Our study answers the question of experimental studies (namely “Wichmann and Sigrist, Journal of neurogenetics 2010”) concerning the sense of the anatomical structure of the T-bar.

Calcium modeling of spine apparatus-containing human dendritic spines demonstrates an “all-or-nothing” communication switch between the spine head and dendrite

Dendritic spines are highly dynamic neuronal compartments that control the synaptic transmission between neurons. Spines form ultrastructural units, coupling synaptic contact sites to the dendritic shaft and often harbor a spine apparatus organelle, composed of smooth endoplasmic reticulum, which is responsible for calcium sequestration and release into the spine head and neck. The spine apparatus has recently been linked to synaptic plasticity in adult human cortical neurons. While the morphological heterogeneity of spines and their intracellular organization has been extensively demonstrated in animal models, the influence of spine apparatus organelles on critical signaling pathways, such as calcium-mediated dynamics, is less well known in human dendritic spines. In this study we used serial transmission electron microscopy to anatomically reconstruct nine human cortical spines in detail as a basis for modeling and simulation of the calcium dynamics between spine and dendrite. The anatomical study of reconstructed human dendritic spines revealed that the size of the postsynaptic density correlates with spine head volume and that the spine apparatus volume is proportional to the spine volume. Using a newly developed simulation pipeline, we have linked these findings to spine-to-dendrite calcium communication. While the absence of a spine apparatus, or the presence of a purely passive spine apparatus did not enable any of the reconstructed spines to relay a calcium signal to the dendritic shaft, the calcium-induced calcium release from this intracellular organelle allowed for finely tuned “all-or-nothing” spine-to-dendrite calcium coupling; controlled by spine morphology, neck plasticity, and ryanodine receptors. Our results suggest that spine apparatus organelles are strategically positioned in the neck of human dendritic spines and demonstrate their potential relevance to the maintenance and regulation of spine-to-dendrite calcium communication.

During the past decade it has become increasingly clear that abnormal synaptic plasticity is a major hallmark of neurological and cognitive disorders. Developing a better understanding of the synaptic plasticity process, which describes the ability of neurons to adapt their contacts in an activity-dependent manner, will lead to improved treatment of many neurological and cognitive disorders. It is known that calcium-dependent events such as synaptic transmission, intracellular calcium release, and calcium wave propagation, are required for many types of synaptic plasticity expression. However, the biological significance of these processes in neurons of the adult human cortex remains unknown. Due to technical limitations and ethical concerns, experimental data addressing this biologically and clinically relevant topic are not available. Therefore, we have implemented a computational model to study the intracellular calcium dynamics in realistic human dendritic spines based on detailed morphological reconstructions. With our model and simulations, we have established the morphological and biological requirements for the propagation of calcium from spines into the dendrites. Our results suggest a critical role for the calcium-storing spine apparatus organelle in regulating calcium homeostasis and propagation in human dendritic spines.

Dynamics and physiological meaning of complexes between ion channels and integrin receptors: the case of Kv11.1

The cellular functions are regulated by a complex interplay of diffuse and local signals. Experimental work in cell physiology has led to recognize that understanding a cell's dynamics requires a deep comprehension of local fluctuations of cytosolic regulators. Macromolecular complexes are major determinants of local signaling. Multi-enzyme assemblies limit the diffusion restriction to reaction kinetics by direct exchange of metabolites. Likewise, close coupling of ion channels and transporters modulate the ion concentration around a channel mouth or transporter binding site. Extreme signal locality is brought about by conformational coupling between membrane proteins, as is typical of mechanotransduction. A paradigmatic case is integrin-mediated cell adhesion. Sensing the extracellular microenvironment and providing an appropriate response is essential in growth and development and has innumerable pathological implications. The process involves bidirectional signal transduction by complex supra-molecular structures that link integrin receptors to ion channels and transporters, growth factor receptors, cytoskeletal elements and other regulatory elements. The dynamics of such complexes is only beginning to be understood. A thoroughly studied example is the association between integrin receptors and the voltage-gated K+ channels Kv11.1. These channels are widely expressed in early embryos, where their physiological roles are poorly understood and apparently different from the shaping of action potential firing in the adult. Hints about these roles come from studies in cancer cells, where Kv11.1 is often overexpressed and appears to re-assume functions, such as controlling cell proliferation/differentiation, apoptosis and migration. Kv11.1 is implicated in these processes through its linking to integrin subunits.

Impact of β-Amyloids Induced Disruption of Ca2+ Homeostasis in a Simple Model of Neuronal Activity

Alzheimer’s disease is characterized by a marked dysregulation of intracellular Ca²⁺ homeostasis. In particular, toxic β-amyloids (Aβ) perturb the activities of numerous Ca²⁺ transporters or channels. Because of the tight coupling between Ca²⁺ dynamics and the membrane electrical activity, such perturbations are also expected to affect neuronal excitability. We used mathematical modeling to systematically investigate the effects of changing the activities of the various targets of Aβ peptides reported in the literature on calcium dynamics and neuronal excitability. We found that the evolution of Ca²⁺ concentration just below the plasma membrane is regulated by the exchanges with the extracellular medium, and is practically independent from the Ca²⁺ exchanges with the endoplasmic reticulum. Thus, disruptions of Ca²⁺ homeostasis interfering with signaling do not affect the electrical properties of the neurons at the single cell level. In contrast, the model predicts that by affecting the activities of L-type Ca²⁺ channels or Ca²⁺-activated K⁺ channels, Aβ peptides promote neuronal hyperexcitability. On the contrary, they induce hypo-excitability when acting on the plasma membrane Ca²⁺ ATPases. Finally, the presence of pores of amyloids in the plasma membrane can induce hypo- or hyperexcitability, depending on the conditions. These modeling conclusions should help with analyzing experimental observations in which Aβ peptides interfere at several levels with Ca²⁺ signaling and neuronal activity.

A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation

Calcium (Ca2+) plays a central role in the excitation and contraction of cardiac myocytes. Experiments have indicated that calcium release is stochastic and regulated locally suggesting the possibility of spatially heterogeneous calcium levels in the cells. This spatial heterogeneity might be important in mediating different signaling pathways. During more than 50 years of computational cell biology, the computational models have been advanced to incorporate more ionic currents, going from deterministic models to stochastic models. While periodic increases in cytoplasmic Ca2+ concentration drive cardiac contraction, aberrant Ca2+ release can underly cardiac arrhythmia. However, the study of the spatial role of calcium ions has been limited due to the computational expense of using a three-dimensional stochastic computational model. In this paper, we introduce a three-dimensional stochastic computational model for rat ventricular myocytes at the whole-cell level that incorporate detailed calcium dynamics, with (1) non-uniform release site placement, (2) non-uniform membrane ionic currents and membrane buffers, (3) stochastic calcium-leak dynamics and (4) non-junctional or rogue ryanodine receptors. The model simulates spark-induced spark activation and spark-induced Ca2+ wave initiation and propagation that occur under conditions of calcium overload at the closed-cell condition, but not when Ca2+ levels are normal. This is considered important since the presence of Ca2+ waves contribute to the activation of arrhythmogenic currents.

Engaging biological oscillators through second messenger pathways permits emergence of a robust gastric slow-wave during peristalsis

Peristalsis, the coordinated contraction—relaxation of the muscles of the stomach is important for normal gastric motility and is impaired in motility disorders. Coordinated electrical depolarizations that originate and propagate within a network of interconnected layers of interstitial cells of Cajal (ICC) and smooth muscle (SM) cells of the stomach wall as a slow-wave, underly peristalsis. Normally, the gastric slow-wave oscillates with a single period and uniform rostrocaudal lag, exhibiting network entrainment. Understanding of the integrative role of neurotransmission and intercellular coupling in the propagation of an entrained gastric slow-wave, important for understanding motility disorders, however, remains incomplete. Using a computational framework constituted of a novel gastric motility network (GMN) model we address the hypothesis that engaging biological oscillators (i.e., ICCs) by constitutive gap junction coupling mechanisms and enteric neural innervation activated signals can confer a robust entrained gastric slow-wave. We demonstrate that while a decreasing enteric neural innervation gradient that modulates the intracellular IP₃ concentration in the ICCs can guide the aboral slow-wave propagation essential for peristalsis, engaging ICCs by recruiting the exchange of second messengers (inositol trisphosphate (IP₃) and Ca²⁺) ensures a robust entrained longitudinal slow-wave, even in the presence of biological variability in electrical coupling strengths. Our GMN with the distinct intercellular coupling in conjunction with the intracellular feedback pathways and a rostrocaudal enteric neural innervation gradient allows gastric slow waves to oscillate with a moderate range of frequencies and to propagate with a broad range of velocities, thus preventing decoupling observed in motility disorders. Overall, the findings provide a mechanistic explanation for the emergence of decoupled slow waves associated with motility impairments of the stomach, offer directions for future experiments and theoretical work, and can potentially aid in the design of new interventional pharmacological and neuromodulation device treatments for addressing gastric motility disorders.

The coordinated contraction and relaxation of the muscles of the stomach, known as peristalsis is important for normal gastric motility and primarily governed by electrical depolarizations that originate and propagate within a network of interconnected layers of interstitial cells of Cajal (ICCs) and smooth muscle cells of the stomach wall as a slow-wave. Under normal conditions, a gastric slow-wave oscillates with a single period and uniform rostrocaudal lag, exhibiting network entrainment. However, the understanding of intrinsic and extrinsic mechanisms that ensure propagation of a robust entrained slow-wave remains incomplete. Here, using a computational framework, we show that in conjunction with an enteric neural innervation gradient along the rostrocaudal ICC chain, and intercellular electrical coupling, the intercellular exchange of inositol trisphosphate between ICCs prevents decoupling by extending the longitudinal entrainment range along the stomach wall, even when variability in intercellular coupling exists. The findings from our study indicate ways that ensure the rostrocaudal spread of a robust gastric slow-wave and provide a mechanistic explanation for the emergence of decoupled slow waves associated with motility impairments of the stomach.

Spatial synaptic modulation through IP3 diffusion triggered by ECB: a computational study with an astrocyte-neurons model

Recently, functional interactions between neurons and astrocytes have been steadily clarified. In particular, the effects of presynaptic depolarization-induced suppression of excitation (DSE) through endocannabinoid (ECB) and endocannabinoids-mediated synaptic potentiation (eSP) by an astrocyte have been used as an evidence of global heterosynaptic modulation. However, the mechanism of occurrence of spatial modulation in a neural network remains unknown. Although the Ca²⁺ density in astrocytes is strongly related to eSP through ECB, the mechanism of the rise in the ECB receptor in Ca²⁺ remains unclear. Since Ca²⁺ is closely related to inositol-1,4,5-trisphosphate (IP₃), it is believed that the released IP₃ affects Ca²⁺ in astrocytes that receive ECB. Therefore, this study approximately showed the spatial distribution of DSE or eSP with astrocyte-neuron computational models, assuming that the IP₃ caused by ECB is transmitted in an astrocyte. The results showed doughnut-shaped DSE, eSP, and DSE regions from the central ECB released points to the surroundings. They suggested that IP₃ diffusion plays an important role in mediating this synaptic modulation.

High rates of calcium-free diffusion in the cytosol of living cells

Calcium (Ca²⁺) is a universal second messenger that participates in the regulation of innumerous physiological processes. The way in which local elevations of the cytosolic Ca²⁺ concentration spread in space and time is key for the versatility of the signals. Ca²⁺ diffusion in the cytosol is hindered by its interaction with proteins that act as buffers. Depending on the concentrations and the kinetics of the interactions, there is a large range of values at which Ca²⁺ diffusion can proceed. Having reliable estimates of this range, particularly of its highest end, which corresponds to the ions free diffusion, is key to understand how the signals propagate. In this work, we present the first experimental results with which the Ca²⁺-free diffusion coefficient is directly quantified in the cytosol of living cells. By means of fluorescence correlation spectroscopy experiments performed in Xenopus laevis oocytes and in cells of Saccharomyces cerevisiae, we show that the ions can freely diffuse in the cytosol at a higher rate than previously thought.

RyR-mediated Ca2+ release elicited by neuronal activity induces nuclear Ca2+ signals, CREB phosphorylation, and Npas4/RyR2 expression

The expression of several hippocampal genes implicated in learning and memory processes requires that Ca²⁺ signals generated in dendritic spines, dendrites, or the soma in response to neuronal stimulation reach the nucleus. The diffusion of Ca²⁺ in the cytoplasm is highly restricted, so neurons must use other mechanisms to propagate Ca²⁺ signals to the nucleus. Here, we present evidence showing that Ca²⁺ release mediated by the ryanodine receptor (RyR) channel type-2 isoform (RyR2) contributes to the generation of nuclear Ca²⁺ signals induced by gabazine (GBZ) addition, glutamate uncaging in the dendrites, or high-frequency field stimulation of primary hippocampal neurons. Additionally, GBZ treatment significantly increased cyclic adenosine monophosphate response element binding protein (CREB) phosphorylation-a key event in synaptic plasticity and hippocampal memory-and enhanced the expression of Neuronal Per Arnt Sim domain protein 4 (Npas4) and RyR2, two central regulators of these processes. Suppression of RyR-mediated Ca²⁺ release with ryanodine significantly reduced the increase in CREB phosphorylation and the enhanced Npas4 and RyR2 expression induced by GBZ. We propose that RyR-mediated Ca²⁺ release induced by neuronal activity, through its contribution to the sequential generation of nuclear Ca²⁺ signals, CREB phosphorylation, Npas4, and RyR2 up-regulation, plays a central role in hippocampal synaptic plasticity and memory processes.

The roles of calcium and ATP in the physiology and pathology of the exocrine pancreas

This review deals with the roles of calcium ions and ATP in the control of the normal functions of the different cell types in the exocrine pancreas as well as the roles of these molecules in the pathophysiology of acute pancreatitis. Repetitive rises in the local cytosolic calcium ion concentration in the apical part of the acinar cells not only activate exocytosis but also, via an increase in the intramitochondrial calcium ion concentration, stimulate the ATP formation that is needed to fuel the energy-requiring secretion process. However, intracellular calcium overload, resulting in a global sustained elevation of the cytosolic calcium ion concentration, has the opposite effect of decreasing mitochondrial ATP production, and this initiates processes that lead to necrosis. In the last few years it has become possible to image calcium signaling events simultaneously in acinar, stellate, and immune cells in intact lobules of the exocrine pancreas. This has disclosed processes by which these cells interact with each other, particularly in relation to the initiation and development of acute pancreatitis. By unraveling the molecular mechanisms underlying this disease, several promising therapeutic intervention sites have been identified. This provides hope that we may soon be able to effectively treat this often fatal disease.

Stochastic reaction-diffusion modeling of calcium dynamics in 3D dendritic spines of Purkinje cells

Calcium (Ca²⁺) is a second messenger assumed to control changes in synaptic strength in the form of both long-term depression and long-term potentiation at Purkinje cell dendritic spine synapses via inositol trisphosphate (IP₃)-induced Ca²⁺ release. These Ca²⁺ transients happen in response to stimuli from parallel fibers (PFs) from granule cells and climbing fibers (CFs) from the inferior olivary nucleus. These events occur at low numbers of free Ca²⁺, requiring stochastic single-particle methods when modeling them. We use the stochastic particle simulation program MCell to simulate Ca²⁺ transients within a three-dimensional Purkinje cell dendritic spine. The model spine includes the endoplasmic reticulum, several Ca²⁺ transporters, and endogenous buffer molecules. Our simulations successfully reproduce properties of Ca²⁺ transients in different dynamical situations. We test two different models of the IP₃ receptor (IP₃R). The model with nonlinear concentration response of binding of activating Ca²⁺ reproduces experimental results better than the model with linear response because of the filtering of noise. Our results also suggest that Ca²⁺-dependent inhibition of the IP₃R needs to be slow to reproduce experimental results. Simulations suggest the experimentally observed optimal timing window of CF stimuli arises from the relative timing of CF influx of Ca²⁺ and IP₃ production sensitizing IP₃R for Ca²⁺-induced Ca²⁺ release. We also model ataxia, a loss of fine motor control assumed to be the result of malfunctioning information transmission at the granule to Purkinje cell synapse, resulting in a decrease or loss of Ca²⁺ transients. Finally, we propose possible ways of recovering Ca²⁺ transients under ataxia.

Effect of crista morphology on mitochondrial ATP output: A computational study

Folding of the mitochondrial inner membrane (IM) into cristae greatly increases the ATP-generating surface area, S _IM, per unit volume but also creates diffusional bottlenecks that could limit reaction rates inside mitochondria. This study explores possible effects of inner membrane folding on mitochondrial ATP output, using a mathematical model for energy metabolism developed by the Jafri group and two- and three-dimensional spatial models for mitochondria, implemented on the Virtual Cell platform. Simulations demonstrate that cristae are micro-compartments functionally distinct from the cytosol. At physiological steady states, standing gradients of ADP form inside cristae that depend on the size and shape of the compartments, and reduce local flux (rate per unit area) of the adenine nucleotide translocase. This causes matrix ADP levels to drop, which in turn reduces the flux of ATP synthase. The adverse effects of membrane folding on reaction fluxes increase with crista length and are greater for lamellar than tubular crista. However, total ATP output per mitochondrion is the product of flux of ATP synthase and S _IM which can be two-fold greater for mitochondria with lamellar than tubular cristae, resulting in greater ATP output for the former. The simulations also demonstrate the crucial role played by intracristal kinases (adenylate kinase, creatine kinase) in maintaining the energy advantage of IM folding.

A first-passage approach to diffusion-influenced reversible binding and its insights into nanoscale signaling at the presynapse

Synaptic transmission between neurons is governed by a cascade of stochastic calcium ion reaction–diffusion events within nerve terminals leading to vesicular release of neurotransmitter. Since experimental measurements of such systems are challenging due to their nanometer and sub-millisecond scale, numerical simulations remain the principal tool for studying calcium-dependent neurotransmitter release driven by electrical impulses, despite the limitations of time-consuming calculations. In this paper, we develop an analytical solution to rapidly explore dynamical stochastic reaction–diffusion problems based on first-passage times. This is the first analytical model that accounts simultaneously for relevant statistical features of calcium ion diffusion, buffering, and its binding/unbinding reaction with a calcium sensor for synaptic vesicle fusion. In particular, unbinding kinetics are shown to have a major impact on submillisecond sensor occupancy probability and therefore cannot be neglected. Using Monte Carlo simulations we validated our analytical solution for instantaneous calcium influx and that through voltage-gated calcium channels. We present a fast and rigorous analytical tool that permits a systematic exploration of the influence of various biophysical parameters on molecular interactions within cells, and which can serve as a building block for more general cell signaling simulators.

Modeling of Astrocyte Networks: Toward Realistic Topology and Dynamics

Neuronal firing and neuron-to-neuron synaptic wiring are currently widely described as orchestrated by astrocytes—elaborately ramified glial cells tiling the cortical and hippocampal space into non-overlapping domains, each covering hundreds of individual dendrites and hundreds thousands synapses. A key component to astrocytic signaling is the dynamics of cytosolic Ca²⁺ which displays multiscale spatiotemporal patterns from short confined elemental Ca²⁺ events (puffs) to Ca²⁺ waves expanding through many cells. Here, we synthesize the current understanding of astrocyte morphology, coupling local synaptic activity to astrocytic Ca²⁺ in perisynaptic astrocytic processes and morphology-defined mechanisms of Ca²⁺ regulation in a distributed model. To this end, we build simplified realistic data-driven spatial network templates and compile model equations as defined by local cell morphology. The input to the model is spatially uncorrelated stochastic synaptic activity. The proposed modeling approach is validated by statistics of simulated Ca²⁺ transients at a single cell level. In multicellular templates we observe regular sequences of cell entrainment in Ca²⁺ waves, as a result of interplay between stochastic input and morphology variability between individual astrocytes. Our approach adds spatial dimension to the existing astrocyte models by employment of realistic morphology while retaining enough flexibility and scalability to be embedded in multiscale heterocellular models of neural tissue. We conclude that the proposed approach provides a useful description of neuron-driven Ca²⁺-activity in the astrocyte syncytium.

A calcium optimum for cytotoxic T lymphocyte and natural killer cell cytotoxicity

Cytotoxic T lymphocytes (CTL) and natural killer (NK) cells are required for host defense. They destroy malignant target cells like cancer cells. Among metal cations, Ca²⁺ plays a prescinded role for CTL and NK cytotoxicity as it is the only cation used as ubiquitous second messenger. Measuring intracellular Ca²⁺ concentrations [Ca²⁺]_int in single cells has greatly changed our understanding of Ca²⁺ signaling. Yet, comparing the role of Ca²⁺ in the pre-[Ca²⁺]_int and [Ca²⁺]_int measurement era reveals that even in the pre-[Ca²⁺]_int measurement era (before 1980), the functions of Ca²⁺ and some other metal cations for the cytotoxic immune response were well established. It was even shown that Ca²⁺ influx across the plasma membrane but not Ca²⁺ release from intracellular sources is relevant for lymphocyte cytotoxicity and that very little Ca²⁺ is needed for efficient lymphocyte cytotoxicity against cancer cells. In the [Ca²⁺]_int measurement era after 1980, many of the important findings were better and more quantitatively refined and in addition the molecules important for Ca²⁺ transport were defined. The unexpected finding that there is a Ca²⁺ optimum of CTL and NK cell cytotoxicity deserves some attention and may be important for anti-cancer therapy.

Quantification of fluctuations from fluorescence correlation spectroscopy experiments in reaction-diffusion systems

Fluorescence correlation spectroscopy (FCS) is commonly used to estimate diffusion and reaction rates. In FCS the fluorescence coming from a small volume is recorded and the autocorrelation function (ACF) of the fluorescence fluctuations is computed. Scaling out the fluctuations due to the emission process, this ACF can be related to the ACF of the fluctuations in the number of observed fluorescent molecules. In this paper the ACF of the molecule number fluctuations is studied theoretically, with no approximations, for a reaction-diffusion system in which the fluorescence changes with binding and unbinding. Theoretical ACFs are usually derived assuming that fluctuations in the number of molecules of one species are instantaneously uncorrelated to those of the others and obey Poisson statistics. Under these assumptions, the ACF derived in this paper is characterized only by the diffusive timescale of the fluorescent species and its total weight is the inverse of the mean number of observed fluorescent molecules. The theory is then scrutinized in view of previous experimental results which, for a similar system, gave a different total weight and correct estimates of other diffusive timescales. The total weight mismatch is corrected by assuming that the variance of the number of fluorescent molecules depends on the variance of the particle numbers of the other species, as in the variance decomposition formula. Including the finite acquisition time in its computation, it is shown that the ACF depends on various timescales of the system and that its total weight coincides with the one obtained with the variance decomposition formula. This calculation implies that diffusion coefficients of nonobservable species can be estimated with FCS experiments performed in reaction-diffusion systems. Ways to proceed in future experiments are also discussed.

Calcium Oscillations in Pancreatic α-cells Rely on Noise and ATP-Driven Changes in Membrane Electrical Activity

In pancreatic α-cells, intracellular Ca²⁺ ([Ca²⁺]_i) acts as a trigger for secretion of glucagon, a hormone that plays a key role in blood glucose homeostasis. Intracellular Ca²⁺ dynamics in these cells are governed by the electrical activity of voltage-gated ion channels, among which ATP-sensitive K⁺ (K_ATP) channels play a crucial role. In the majority of α-cells, the global Ca²⁺ response to lowering external glucose occurs in the form of oscillations that are much slower than electrical activity. These Ca²⁺ oscillations are highly variable as far as inter-spike intervals, shapes and amplitudes are concerned. Such observations suggest that Ca²⁺ dynamics in α-cells are much influenced by noise. Actually, each Ca²⁺ increase corresponds to multiple cycles of opening/closing of voltage gated Ca²⁺ channels that abruptly become silent, before the occurrence of another burst of activity a few tens of seconds later. The mechanism responsible for this intermittent activity is currently unknown. In this work, we used computational modeling to investigate the mechanism of cytosolic Ca²⁺ oscillations in α-cells. Given the limited population of K_ATP channels in this cell type, we hypothesized that the stochastic activity of these channels could play a key role in the sporadic character of the action potentials. To test this assumption, we extended a previously proposed model of the α-cells electrical activity (Diderichsen and Göpel, 2006) to take Ca²⁺ dynamics into account. Including molecular noise on the basis of a Langevin type description as well as realistic dynamics of opening and closing of K_ATP channels, we found that stochasticity at the level of the activity of this channel is on its own not able to produce Ca²⁺ oscillations with a time scale of a few tens of seconds. However, when taking into account the intimate relation between Ca²⁺ and ATP changes together with the intrinsic noise at the level of the K_ATP channels, simulations displayed Ca²⁺ oscillations that are compatible with experimental observations. We analyzed the detailed mechanism and used computational simulations to identify the factors that can affect Ca²⁺ oscillations in α-cells.

Vascular Kv7 channels control intracellular Ca2+ dynamics in smooth muscle

Voltage-gated Kv7 (or KCNQ) channels control activity of excitable cells, including vascular smooth muscle cells (VSMCs), by setting their resting membrane potential and controlling other excitability parameters. Excitation-contraction coupling in muscle cells is mediated by Ca2+ but until now, the exact role of Kv7 channels in cytosolic Ca2+ dynamics in VSMCs has not been fully elucidated. We utilised microfluorimetry to investigate the impact of Kv7 channel activity on intracellular Ca2+ levels and electrical activity of rat A7r5 VSMCs and primary human internal mammary artery (IMA) SMCs. Both, direct (XE991) and G protein coupled receptor mediated (vasopressin, AVP) Kv7 channel inhibition induced robust Ca2+ oscillations, which were significantly reduced in the presence of Kv7 channel activator, retigabine, L-type Ca2+ channel inhibitor, nifedipine, or T-type Ca2+ channel inhibitor, NNC 55-0396, in A7r5 cells. Membrane potential measured using FluoVolt exhibited a slow depolarisation followed by a burst of sharp spikes in response to XE991; spikes were temporally correlated with Ca2+ oscillations. Phospholipase C inhibitor (edelfosine) reduced AVP-induced, but not XE991-induced Ca2+ oscillations. AVP and XE991 induced a large increase of [Ca2+]i in human IMA, which was also attenuated with retigabine, nifedipine and NNC 55-0396. RT-PCR, immunohistochemistry and electrophysiology suggested that Kv7.5 was the predominant Kv7 subunit in both rat and human arterial SMCs; CACNA1C (Cav1.2; L-type) and CACNA1 G (Cav3.1; T-type) were the most abundant voltage-gated Ca2+ channel gene transcripts in both types of VSMCs. This study establishes Kv7 channels as key regulators of Ca2+ signalling in VSMCs with Kv7.5 playing a dominant role.