Multiscale modeling of calcium signaling in the cardiac dyad

Raimond Winslow, Antti Tanskanen, Mindao Chen, Joseph L. Greenstein

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Calcium (Ca2+)-induced Ca2+-release (CICR) takes place in spatially restricted microdomains known as dyads. The length scale over which CICR occurs is on the order of nanometers and relevant time scales range from micro- to milliseconds. Quantitative understanding of CICR therefore requires development of models that are applicable over a range of spatio-temporal scales. We will present several new approaches for multiscale modeling of CICR. First, we present a model of dyad Ca2+ dynamics in which the Fokker-Planck equation (FPE) is solved for the probability P(x, t) that a Ca2+ ion is located at dyad position x at time t. Using this model, we demonstrate that (a) Ca2+ signaling in the dyad is mediated by approximately tens of Ca2+ ions; (b) these signaling events are noisy due to the small number of ions involved; and (c) the geometry of the RyR (ryanodine receptors) protein may function to restrict the diffusion of and to "funnel" Ca2+ ions to activation-binding sites on the RyR, thus increasing RyR open probability and excitation-contraction (EC) coupling gain. Simplification of this model to one in which the dyadic space is represented using a single compartment yields the stochastic local-control model of CICR developed previously. We have shown that this model captures fundamental properties of CICR, such as graded release and voltage-dependent gain, may be integrated within a model of the myocyte and may be simulated in reasonable times using a combination of efficient numerical methods and parallel computing, but remains too complex for general use in cell simulations. To address this problem, we show how separation of time scales may be used to formulate a model in which nearby L-type Ca2+ channels (LCCs) and RyRs gate as a coupled system that may be described using low-dimensional systems of ordinary differential equations, thus reducing computational complexity while capturing fundamentally important properties of CICR. The simplified model may be solved many orders of magnitude faster than can either of the more detailed models, thus enabling incorporation into tissue-level simulations.

Original languageEnglish (US)
Title of host publicationAnnals of the New York Academy of Sciences
Pages362-375
Number of pages14
Volume1080
DOIs
StatePublished - Oct 2006

Publication series

NameAnnals of the New York Academy of Sciences
Volume1080
ISSN (Print)00778923
ISSN (Electronic)17496632

Fingerprint

Calcium Signaling
Ryanodine Receptor Calcium Release Channel
Calcium
Ions
Excitation Contraction Coupling
Muscle Cells
Binding Sites
Fokker Planck equation
Parallel processing systems
Ordinary differential equations
Computational complexity
Numerical methods
Proteins
Chemical activation
Tissue

Keywords

  • Calcium-induced
  • Calcium-release
  • Dyad
  • Excitation-contraction coupling
  • Modeling

ASJC Scopus subject areas

  • Biochemistry, Genetics and Molecular Biology(all)

Cite this

Winslow, R., Tanskanen, A., Chen, M., & Greenstein, J. L. (2006). Multiscale modeling of calcium signaling in the cardiac dyad. In Annals of the New York Academy of Sciences (Vol. 1080, pp. 362-375). (Annals of the New York Academy of Sciences; Vol. 1080). https://doi.org/10.1196/annals.1380.027

Multiscale modeling of calcium signaling in the cardiac dyad. / Winslow, Raimond; Tanskanen, Antti; Chen, Mindao; Greenstein, Joseph L.

Annals of the New York Academy of Sciences. Vol. 1080 2006. p. 362-375 (Annals of the New York Academy of Sciences; Vol. 1080).

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Winslow, R, Tanskanen, A, Chen, M & Greenstein, JL 2006, Multiscale modeling of calcium signaling in the cardiac dyad. in Annals of the New York Academy of Sciences. vol. 1080, Annals of the New York Academy of Sciences, vol. 1080, pp. 362-375. https://doi.org/10.1196/annals.1380.027
Winslow R, Tanskanen A, Chen M, Greenstein JL. Multiscale modeling of calcium signaling in the cardiac dyad. In Annals of the New York Academy of Sciences. Vol. 1080. 2006. p. 362-375. (Annals of the New York Academy of Sciences). https://doi.org/10.1196/annals.1380.027
Winslow, Raimond ; Tanskanen, Antti ; Chen, Mindao ; Greenstein, Joseph L. / Multiscale modeling of calcium signaling in the cardiac dyad. Annals of the New York Academy of Sciences. Vol. 1080 2006. pp. 362-375 (Annals of the New York Academy of Sciences).
@inproceedings{dfd09f0a132d4fb6847449dcbb2495cc,
title = "Multiscale modeling of calcium signaling in the cardiac dyad",
abstract = "Calcium (Ca2+)-induced Ca2+-release (CICR) takes place in spatially restricted microdomains known as dyads. The length scale over which CICR occurs is on the order of nanometers and relevant time scales range from micro- to milliseconds. Quantitative understanding of CICR therefore requires development of models that are applicable over a range of spatio-temporal scales. We will present several new approaches for multiscale modeling of CICR. First, we present a model of dyad Ca2+ dynamics in which the Fokker-Planck equation (FPE) is solved for the probability P(x, t) that a Ca2+ ion is located at dyad position x at time t. Using this model, we demonstrate that (a) Ca2+ signaling in the dyad is mediated by approximately tens of Ca2+ ions; (b) these signaling events are noisy due to the small number of ions involved; and (c) the geometry of the RyR (ryanodine receptors) protein may function to restrict the diffusion of and to {"}funnel{"} Ca2+ ions to activation-binding sites on the RyR, thus increasing RyR open probability and excitation-contraction (EC) coupling gain. Simplification of this model to one in which the dyadic space is represented using a single compartment yields the stochastic local-control model of CICR developed previously. We have shown that this model captures fundamental properties of CICR, such as graded release and voltage-dependent gain, may be integrated within a model of the myocyte and may be simulated in reasonable times using a combination of efficient numerical methods and parallel computing, but remains too complex for general use in cell simulations. To address this problem, we show how separation of time scales may be used to formulate a model in which nearby L-type Ca2+ channels (LCCs) and RyRs gate as a coupled system that may be described using low-dimensional systems of ordinary differential equations, thus reducing computational complexity while capturing fundamentally important properties of CICR. The simplified model may be solved many orders of magnitude faster than can either of the more detailed models, thus enabling incorporation into tissue-level simulations.",
keywords = "Calcium-induced, Calcium-release, Dyad, Excitation-contraction coupling, Modeling",
author = "Raimond Winslow and Antti Tanskanen and Mindao Chen and Greenstein, {Joseph L.}",
year = "2006",
month = "10",
doi = "10.1196/annals.1380.027",
language = "English (US)",
isbn = "1573316512",
volume = "1080",
series = "Annals of the New York Academy of Sciences",
pages = "362--375",
booktitle = "Annals of the New York Academy of Sciences",

}

TY - GEN

T1 - Multiscale modeling of calcium signaling in the cardiac dyad

AU - Winslow, Raimond

AU - Tanskanen, Antti

AU - Chen, Mindao

AU - Greenstein, Joseph L.

PY - 2006/10

Y1 - 2006/10

N2 - Calcium (Ca2+)-induced Ca2+-release (CICR) takes place in spatially restricted microdomains known as dyads. The length scale over which CICR occurs is on the order of nanometers and relevant time scales range from micro- to milliseconds. Quantitative understanding of CICR therefore requires development of models that are applicable over a range of spatio-temporal scales. We will present several new approaches for multiscale modeling of CICR. First, we present a model of dyad Ca2+ dynamics in which the Fokker-Planck equation (FPE) is solved for the probability P(x, t) that a Ca2+ ion is located at dyad position x at time t. Using this model, we demonstrate that (a) Ca2+ signaling in the dyad is mediated by approximately tens of Ca2+ ions; (b) these signaling events are noisy due to the small number of ions involved; and (c) the geometry of the RyR (ryanodine receptors) protein may function to restrict the diffusion of and to "funnel" Ca2+ ions to activation-binding sites on the RyR, thus increasing RyR open probability and excitation-contraction (EC) coupling gain. Simplification of this model to one in which the dyadic space is represented using a single compartment yields the stochastic local-control model of CICR developed previously. We have shown that this model captures fundamental properties of CICR, such as graded release and voltage-dependent gain, may be integrated within a model of the myocyte and may be simulated in reasonable times using a combination of efficient numerical methods and parallel computing, but remains too complex for general use in cell simulations. To address this problem, we show how separation of time scales may be used to formulate a model in which nearby L-type Ca2+ channels (LCCs) and RyRs gate as a coupled system that may be described using low-dimensional systems of ordinary differential equations, thus reducing computational complexity while capturing fundamentally important properties of CICR. The simplified model may be solved many orders of magnitude faster than can either of the more detailed models, thus enabling incorporation into tissue-level simulations.

AB - Calcium (Ca2+)-induced Ca2+-release (CICR) takes place in spatially restricted microdomains known as dyads. The length scale over which CICR occurs is on the order of nanometers and relevant time scales range from micro- to milliseconds. Quantitative understanding of CICR therefore requires development of models that are applicable over a range of spatio-temporal scales. We will present several new approaches for multiscale modeling of CICR. First, we present a model of dyad Ca2+ dynamics in which the Fokker-Planck equation (FPE) is solved for the probability P(x, t) that a Ca2+ ion is located at dyad position x at time t. Using this model, we demonstrate that (a) Ca2+ signaling in the dyad is mediated by approximately tens of Ca2+ ions; (b) these signaling events are noisy due to the small number of ions involved; and (c) the geometry of the RyR (ryanodine receptors) protein may function to restrict the diffusion of and to "funnel" Ca2+ ions to activation-binding sites on the RyR, thus increasing RyR open probability and excitation-contraction (EC) coupling gain. Simplification of this model to one in which the dyadic space is represented using a single compartment yields the stochastic local-control model of CICR developed previously. We have shown that this model captures fundamental properties of CICR, such as graded release and voltage-dependent gain, may be integrated within a model of the myocyte and may be simulated in reasonable times using a combination of efficient numerical methods and parallel computing, but remains too complex for general use in cell simulations. To address this problem, we show how separation of time scales may be used to formulate a model in which nearby L-type Ca2+ channels (LCCs) and RyRs gate as a coupled system that may be described using low-dimensional systems of ordinary differential equations, thus reducing computational complexity while capturing fundamentally important properties of CICR. The simplified model may be solved many orders of magnitude faster than can either of the more detailed models, thus enabling incorporation into tissue-level simulations.

KW - Calcium-induced

KW - Calcium-release

KW - Dyad

KW - Excitation-contraction coupling

KW - Modeling

UR - http://www.scopus.com/inward/record.url?scp=33845667860&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33845667860&partnerID=8YFLogxK

U2 - 10.1196/annals.1380.027

DO - 10.1196/annals.1380.027

M3 - Conference contribution

C2 - 17132795

AN - SCOPUS:33845667860

SN - 1573316512

SN - 9781573316514

VL - 1080

T3 - Annals of the New York Academy of Sciences

SP - 362

EP - 375

BT - Annals of the New York Academy of Sciences

ER -