Theoretical Study of the Function of the IP3 Receptor / BK Channel Complex in a Single Neuron


  • María Eugenia Pérez Bonilla Benemérita Universidad Autónoma de Puebla
  • Marleni Reyes Monreal Benemérita Universidad Autónoma de Puebla
  • Jessica Quintero Pérez Benemérita Universidad Autónoma de Puebla
  • Miguel Pérez Escalera Benemérita Universidad Autónoma de Puebla
  • Arturo Reyes Lazalde Benemérita Universidad Autónoma de Puebla


BK channels, IP3 receptor, Ca2 microdomains, Single neuron simulations, SERCA pump


Large conductance calcium-activated potassium (BK) channels carry out many functions in the central nervous system. The opening of BK channels requires a rise in the cytosolic calcium ([Ca2+]cyt) concentration, which can occur in two ways: calcium influx from voltage-gated calcium channels (VGCCs) located on the plasma membrane and calcium efflux through the endoplasmic reticulum (ER) membrane to the cytosol triggered by inositol 1,4,5-trisphosphate (IP3) receptors (IP3-Rs) and ryanodine receptors (RyRs). The BK channel/IP3-R/RyR interaction has been widely reported in smooth muscle but scarce information exist on neurons, where its presence is uncertain. The aim of this study was to develop a computational model of a neuron to replicate the interaction between the release of Ca2+ from the ER (through IP3-Rs and RyRs) and the opening of BK channels on the plasma membrane to regulate the level of [Ca2+]cyt, based on the Hodgkin-Huxley formalism and the Goldbeter model. The mathematical models were implemented on Visual Basic® and differential equations were solved numerically. Various conditions of BK conductance and the efflux of endoplasmic Ca2+ were explored. The results show that an abrupt increase in [Ca2+]cyt (≥ 5 mM) activates the BK channels and either pauses or stops the action potential train.


Download data is not yet available.

Author Biographies

María Eugenia Pérez Bonilla, Benemérita Universidad Autónoma de Puebla

María Eugenia Pérez Bonilla got her Bachelor of Medicine and Master of Physiology from the Benemérita Universidad Autónoma de Puebla. She has a doctorate in experimental pathology from the Center for Research and Advanced Studies of the Instituto Politecnico Nacional. She is part of a multidisciplinary group with an interest in development of interactive software for education and research in biomedical sciences.

Marleni Reyes Monreal, Benemérita Universidad Autónoma de Puebla

Marleni Reyes Monreal is a computer technician, she received her Bachelor Degree in Graphic Design and her Master Degree in Aesthetics and Art from Benemérita Universidad Autónoma de Puebla. She has a Master Degree in Educational Technology (multimedia design and simulators) and she is currently working to obtain her doctorate in Ecoeducation (simulators) and a second one in History (Art). She is part of a multidisciplinary group with an interest in the design and aesthetic development of interactive software for education and research.

Jessica Quintero Pérez, Benemérita Universidad Autónoma de Puebla

Jessica Quintero Pérez has a degree in Physiotherapy from the Benemérita Universidad Autónoma de Puebla. He has a master's degree in Manual Physiotherapy of the Locomotor System from the Universidad de Alcalá in Madrid, Spain. She is a PhD candidate in the Biomechanics and Bioengineering applied to health program at the University of Alcalá in Madrid, Spain. Active teacher of the Degree in Physiotherapy at the Benemérita Universidad Autónoma de Puebla. Professor of the Diploma in Physiotherapy in Cervical-Temporomandibular Disorders at the Universidad Cuauhtémoc, Querétaro. Teacher and monitor of the Master's degree in Manual Physiotherapy of the Locomotor System, Mexico 2017-2020 headquarters. Speaker and workshop leader in national and international Physiotherapy congresses. She is part of a multidisciplinary group with an interest in the development of interactive software for education and research in Biomedical Sciences.

Miguel Pérez Escalera, Benemérita Universidad Autónoma de Puebla

Miguel Pérez Escalera got his Bachelor Degree in Computer Science from Benemérita Universidad Autónoma de Puebla, he has a Master Degree in Computer Science from the same institution and he is working for his doctorate, also in Computer Science in Universidad de las Americas Puebla. He is part of a multidisciplinary group with an interest in the development of interactive software and 3D systems for education and research.

Arturo Reyes Lazalde, Benemérita Universidad Autónoma de Puebla

Arturo Reyes Lazalde received his Bachelor of Medicine and Master of Physiology from Benemérita Universidad Autónoma de Puebla. He has a doctorate in basic biomedical research (neuroscience) from Universidad Nacional Autónoma de México. He is part of a multidisciplinary group with an interest in the development of interactive software for education and research in biomedical sciences.


Sausbier U, Sausbier M, Sailer CA, Arntz C, Knaus HG, Neuhuber W, et al. Ca2+-activated K+ channels of the BK-type in the mouse brain. Histochem Cell Biol [Internet]. 2005;125(6):725–41. Available from:

Faber ESL, Sah P. Calcium-Activated Potassium Channels : Multiple Contributions to Neuronal Function. Neurosci [Internet]. 2003;9(3):181–94. Available from:

Meech RW. Intracellular calcium injection causes increased potassium conductance in Aplysia nerve cells. Comp Biochem Physiol Part A Physiol [Internet]. 1972;42(2):493–9. Available from:

Latorre R, Castillo K, Carrasquel-Ursulaez W, Sepulveda RV, Gonzalez-Nilo F, Gonzalez C, et al. Molecular determinants of BK channel functional diversity and functioning. Physiol Rev [Internet]. 2017;97(1):39–87. Available from:

Kshatri AS, Gonzalez-Hernandez A, Giraldez T. Physiological Roles and Therapeutic Potential of Ca2+ Activated Potassium Channels in the Nervous System. Front Mol Neurosci [Internet]. 2018;11:1–18. Available from:

Lee US, Cui J. BK channel activation : structural and functional insights. Trends Neurosci [Internet]. 2010;33(9):415–23. Available from:

Horrigan FT, Aldrich RW. Coupling between voltage sensor activation , Ca2+ binding and channel opening in large conductance (BK) potassium channels. J Gen Physiol [Internet]. 2002;120:267–305. Available from:

Ha GE, Cheong E. Spike frequency adaptation in neurons of the central nervous system. Exp Neurobiol [Internet]. 2017;26(4):179–85. Available from:

Bock T, Stuart GJ. The Impact of BK Channels on Cellular Excitability Depends on their Subcellular Location. Front Cell Neurosci [Internet]. 2016;10:1–8. Available from:

Gu N, Vervaeke K, Storm JF. BK potassium channels facilitate high-frequency firing and cause early spike frequency adaptation in rat CA1 hippocampal pyramidal cells. J Physiol [Internet]. 2007;3:859–82. Available from:

Berkefeld H, Fakler B, Schulte U. Ca2+-Activated K+ Channels: From Protein Complexes to Function. Physiol Rev [Internet]. 2010;90(4):1437–59. Available from:

Misonou H, Menegola M, Buchwalder L, Park E, Meredith A, Rhodes K, et al. Immunolocalization of the Ca2+ -activated K+ channel Slo1 in axons and nerve terminals of mammalian brain and cultured neurons. J Comp Neurol [Internet]. 2006;496(3):289–302. Available from:

Wang Z-W. Regulation of Synaptic Transmission by Presynaptic CaMKII and BK Channels. Mol Neurobiol [Internet]. 2008;38(2):153–66. Available from:

Griguoli M, Sgritta M, Cherubini E. Presynaptic BK channels control transmitter release: physiological relevance and potential therapeutic implications. J Physiol [Internet]. 2016;13:3489–500. Available from:

Bygrave FL, Benedetti A. What is the concentration of calcium ions in the endoplasmic reticulum? Cell Calcium [Internet]. 1996;19(6):547–51. Available from:

Berridge MJ. Calcium microdomains: Organization and function. Cell Calcium [Internet]. 2006;40(5–6):405–12. Available from:

Dopico AM, Bukiya AN, Jaggar JH. Calcium- and voltage-gated BK channels in vascular smooth muscle. Pflügers Archiv [Internet]. 2018;470(9):1271–89. Available from:

Fakler B, Adelman JP. Control of KCa Channels by Calcium Nano/Microdomains. Neuron [Internet]. 2008;59(6):873–81. Available from:

Samtleben S, Jaepel J, Fecher C, Andreska T, Rehberg M, Blum R. Direct imaging of ER calcium with targeted-esterase induced dye loading (TED). J Vis Exp [Internet]. 2013;(75):1–17. Available from:

Verkhratsky A. Physiology and Pathophysiology of the Calcium Store in the Endoplasmic Reticulum of Neurons. Physiol Rev [Internet]. 2005;85(1):201–79. Available from:

Karagas NE, Venkatachalam K. Roles for the Endoplasmic Reticulum in Regulation of Neuronal Calcium Homeostasis. Cells [Internet]. 2019;8(10):1232. Available from:

Foskett JK, White C, Cheung K-H, Mak D-OD. Inositol Trisphosphate Receptor Ca2+ Release Channels. Physiol Rev [Internet]. 2007;87(2):593–658. Available from:

Irie T, Trussell LO. Double-Nanodomain Coupling of Calcium Channels, Ryanodine Receptors, and BK Channels Controls the Generation of Burst Firing. Neuron [Internet]. 2017;96(4):856-870.e4. Available from:

Wu Y, Whiteus C, Xu CS, Hayworth KJ, Weinberg RJ, Hess HF, et al. Contacts between the endoplasmic reticulum and other membranes in neurons. Proc Natl Acad Sci U S A [Internet]. 2017;114(24):E4859–67. Available from:

Terasaki M, Shemesh T, Kasthuri N, Klemm RW, Hayworth KJ, Hand AR, et al. Stacked endoplasmic reticulum sheets are connected by helicoidal membrane motifs. Cell [Internet]. 2013;154(2):285–96. Available from:

Takahashi K, Wood RL. Subsurface cisterns in the Purkinje cells of cerebellum of Syrian hamster. Zeitschrift für Zellforsch und Mikroskopische Anat [Internet]. 1970;110(3):311–20. Available from:

Pan NC, Bai YF, Yang Y, Hökfelt T, Xu ZQD. Activation of galanin receptor 2 stimulates large conductance Ca2+ -dependent K+ (BK) channels through the IP3 pathway in human embryonic kidney (HEK293) cells. Biochem Biophys Res Commun [Internet]. 2014;446(1):316–21. Available from:

Mujica PE, González FG. Interaction between IP3 receptors and BK channels in arterial smooth muscle: Non-canonical IP3 signaling at work. J Gen Physiol [Internet]. 2011;137(5):473–7. Available from:

Bootman MD, Collins TJ, Peppiatt CM, Prothero LS, MacKenzie L, De Smet P, et al. Calcium signalling - An overview. Semin Cell Dev Biol [Internet]. 2001;12(1):3–10. Available from:

Larkum ME, Watanabe S, Nakamura T, Lasser-Ross N, Ross WN. Synaptically Activated Ca2+ Waves in Layer 2/3 and Layer 5 Rat Neocortically Pyramidal Neurons. J Physiol [Internet]. 2003;549(2):471–88. Available from:

Ross WN. Understanding calcium waves and sparks in central neurons. Nat Rev Neuoscience [Internet]. 2015;13(3):157–68. Available from:

Goldbeter A, Dupont G, Berridge MJ. Minimal model for signal-induced Ca2+ oscillations and for their frequency encoding through protein phosphorylation. Proc Natl Acad Sci U S A [Internet]. 1990;87(4):1461–5. Available from:

De Young GW, Keizer J. A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration. Proc Natl Acad Sci U S A [Internet]. 1992;89(20):9895–9. Available from:

Li Y-X, Rinzel J. Equations for InsP3 Receptor-mediated [Ca2+]i Oscillations Derived from a Detailed Kinetic Model: A Hodgkin-Huxley Like Formalism. J Theor Biol [Internet]. 1994;166(4):461–73. Available from:

Bezprozvanny I, Ehrlich BE. Inositol (1,4,5)-trisphosphate (Insp3)-gated Ca channels from cerebellum: Conduction properties for divalent cations and regulation by intraluminal calcium. J Gen Physiol [Internet]. 1994;104(5):821–56. Available from:

Keizer J, Li YX, Stojilkovic S, Rinzel J. Essay InsP3-induced Ca2+ excitability of the endoplasmic reticulum. Mol Biol Cell [Internet]. 1995;6:945–51. Available from:

Tang Y, Stephenson JL, Othmer HG. Simplification and Analysis of Models of Calcium Dynamics Based on IP3-Sensitive Calcium Channel Kinetics. Biophys J [Internet]. 1996;70(1):246–63. Available from:

Sneyd J, Falcke M. Models of the inositol trisphosphate receptor. Prog Biophys Mol [Internet]. 2005;89(3):207–45. Available from:

Kötter R. Neuroscience Databases: A Practical Guide. New York: Springer; 2003. 310 p.

Hituri K, Linne M-L. Comparison of Models for IP3 Receptor Kinetics Using Stochastic Simulations. PLoS One [Internet]. 2013;8(4): e59618. Available from:

Edgerton JR, Reinhart PH. Distinct contributions of small and large conductance Ca2+- activated K+ channels to rat Purkinje neuron function. J Physiol [Internet]. 2003;548(1):53–69. Available from:

Blaustein MP, Golovina VA. Structural complexity and functional diversity of endoplasmic reticulum Ca2+ stores. Trends Neurosci [Internet]. 2001;24(10):602–8. Available from:

Sterratt D, Graham B, Gillies A, Willshaw D. Principles of Computational Modelling in Neuroscience. Cambridge: Cambridge University Press; 2011. 300p.

Traub RD, Miles R. Neuronal Networks Hippocampus. Cambridge: Cambridge University Press; 1991. 281p.

Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in Nerve. J Physiol [Internet]. 1952;117(4):500–44. Available from:

McCarron JG, Chalmers S, Bradley KN, MacMillan D, Muir TC. Ca2+ microdomains in smooth muscle. Cell Calcium [Internet]. 2006;40(5–6):461–93. Available from:

Blackwell KT. Approaches and tools for modeling signaling pathways and calcium dynamics in neurons. J Neurosci Methods [Internet]. 2013;220(2):131-140. Available from:

Chen-Engerer HJ, Hartmann J, Karl RM, Yang J, Feske S, Konnerth A. Two types of functionally distinct Ca2+ stores in hippocampal neurons. Nat Commun [Internet]. 2019;10(1):1–8. Available from:

Dupont G, Goldbeter A. One-pool model for Ca2+ oscillations involving Ca2+ and inositol 1,4,5-triphosphate as co-agonists for Ca2+ release. Cell Calcium [Internet]. 1993;14(4):311–22. Available from:

Blackwell KT. Modeling calcium concentration and biochemical reactions. Brains, Minds, and Media [Internet]. 2005;1:1–27. Available from:

Zill DG. Ecuaciones diferenciales con aplicaciones. México: Grupo Editorial Iberoamérica; 1988. 516p. Spanish.

Constantinides A. Applied Numerical Methods with Personal Computers. New York, USA: McGraw-Hill International Editions; 1987. 626p.

Foskett JK, White C, Cheung K-H, Mak D-OD. Inositol Trisphosphate Receptor Ca2+ release channels. Physiol Rev [Internet]. 2007;87(2):593–658. Available from:

Khodakhah K, Ogden D. Fast activation and inactivation of inositol trisphosphate‐evoked Ca2+ release in rat cerebellar Purkinje neurones. J Physiol [Internet]. 1995;487(2):343–58. Available from:

Cheng H, Lederer WJ. Calcium Sparks. Physiol Rev [Internet]. 2008;88(4):1491–545. Available from:

Naraghi M, Neher E. Linearized Buffered Ca2+ Diffusion in Microdomains and Its Implications for Calculation of [Ca2+] at the Mouth of a Calcium Channel. J Neurosci [Internet]. Available from: 1997;17(18):6961–73. Available from:

Rizzuto R, Pozzan T. Microdomains of Intracellular Ca2+: Molecular Determinants and Functional Consequences. Physiol Rev [Internet]. 2006;86(1):369–408. Available from:

Weaver AK, Olsen ML, McFerrin MB, Sontheimer H. BK Channels Are Linked to Inositol 1,4,5-Triphosphate Receptors via Lipid Rafts: A novel mechanism for coupling [Ca2+]i to ion channel activation. J Biol Chem [Internet]. 2007;282(43):31558–68. Available from:

Koch C, Segev I. Methods in Neuronal Modeling. Cambridge: The MIT Press; 1999. 671p.

Berkefeld H, Fakler B. Repolarizing Responses of BKCa-Cav Complexes Are Distinctly Shaped by Their Cav Subunits. J Neurosci [Internet]. 2008;28(33):8238–45. Available from:

Hou P, Xiao F, Liu H, Yuchi M, Zhang G, Wu Y, et al. Extrapolating microdomain Ca2+ dynamics using BK channels as a Ca2+ sensor. Sci Rep [Internet]. 2016;6:17343. Available from:

Meyer T, Stryer L. Molecular model for receptor-stimulated calcium spiking. Proc Natl Acad Sci U S A [Internet]. 1988;85(14):5051–5. Available from:

Clements MA, Swapna I, Morikawa H. Inositol 1,4,5-Triphosphate Drives Glutamatergic and Cholinergic Inhibition Selectively in Spiny Projection Neurons in the Striatum. J Neurosci [Internet]. 2013;33(6):2697–708. Available from:




How to Cite

Pérez Bonilla, M. E., Reyes Monreal, M., Quintero Pérez, J., Pérez Escalera, M., & Reyes Lazalde, A. (2021). Theoretical Study of the Function of the IP3 Receptor / BK Channel Complex in a Single Neuron. Mexican Journal of Biomedical Engineering, 42(2), 15–31. Retrieved from



Research Articles

Share on:

Dimensions Citation