GPU-Enabled Spatiotemporal Model of Stochastic Cardiac Calcium Dynamics and Arrhythmias
|
|
- Primrose Adams
- 5 years ago
- Views:
Transcription
1 GPU-Enabled Spatiotemporal Model of Stochastic Cardiac Calcium Dynamics and Arrhythmias M. Saleet Jafri Hoang Trong Minh Tuan Department of Bioinformatics and Computational Biology George Mason University Institute of Computational Medicine, The Johns Hopkins University Department of Biomedical Engineering and Technology The University of Maryland Baltimore
2 Understanding Ca 2+ -Dependent Cardiac Arrhythmias Heart disease is the leading cause of death in developed nations and an increasing problem in the developing world. The contraction of heart muscle pumps blood throughout the body to supply oxygen and nutrients to and remove carbon dioxide and waste products from the body tissues. Death often occurs by cardiac arrhythmias that prevent this normal function of the heart.
3 Basic Components of Cardiac E-C Coupling Membrane Currents Calcium Handling Force Generation Actio n Potential Calcium Transient Force Transient
4 Understanding Ca 2+ -Dependent Cardiac Arrhythmias As cardiac myocyte calcium dynamics are essential for the contraction of the heart, the dysfunction of normal calcium dynamics is often a major factor in cardiac arrhythmias. Computational models are an essential tool to understand the complex dynamcs of cardiac myocyte calcium signaling as the explain mechanisms by integrating the known information about this system.
5 Ca 2+ Transient and Ca 2+ Sparks in Ventricular Myocytes [Ca 2+ ] "sparks" are the elementary release events. They are synchronized by the electrical signal of the cell to produce the elevation of [Ca 2+ ].
6 What is a Ca 2+ spark? seconds cell images at 0.5 sec per image sparks line-scan image at 2 ms per line location time (from Cheng, Lederer & Cannell (1993), Science 262:740) spark
7 (From L. Fernando Santana, unpublished) Heart Cell
8 Sarcomere Geometry Many diads Diads separated NSR connections All sarcomeres shorten uniformly
9 T-Tubules and SR Apposition Z-line SR M RyRs T-tubule TT-SR junction Modified from J. Frank (1990)
10 Ryanodine Receptor Organization from Baddeley, et al., micron
11 RYR'S RYR'S RyR Channel Properties Three properties of RyR gating need to be included in the model. 1. Large number of RyRs (Franzini-Armstrong et al., 1998) 2. SR lumenal [Ca 2+ ] 3. Coupled gating of RyRs [Ca 2+ ] lumen and RyR gating Coupled gating of RyRs control trans [Ca 2+ ]=20 mm trans [Ca 2+ ]=5 mm FK506 from Gyorke & Gyorke (1998) Biophys J. 75:280 Skeletal Muscle RyRs: Marx et al., (1998) Science 281:818. Heart RyRs: Gaburjakova et al. (2001) Biophys. J. 80:380A.
12 Approach Contraction of the heart is caused by the summation of calcium sparks. Build a model of contraction in the heart starting with these stochastic events. Use this model to answer fundamental questions about the mechanisms of arrhythmia that could not be answered with previous modeling efforts.
13 Questions What is the molecular basis of the calcium leak from the SR? How does the leak play a role in arrhythmogenisis?
14 Computational Challenges Stochastic Simulation is very computationally expensive. Various approaches have been suggested to address these. Reduction methods Reduction methods make assumptions about the system to make the reduction possible. These might not be valid under all relevant conditions. Sometimes the reduction require simplification of the system by reducing model complexity. This limits the details that can be included. Monte Carlo Simulation of simpler system Other attempts have simplified the dynamics of ryanodine receptor gating or the system. Using slower kinetics, fewer channels, release units, omitting physiological/biophysical detail reduced veracity of the model. We have developed the Ultrafast Monte Carlo Method that makes the computation possible
15 New Compartment Model Benchmarks Ultra-fast Ultra-fast Method Monte Carlo Monte Carlo Speedup On CPU On GPU No I/O 69 min 3:40 min 19 x I/O 70:52 min 4:26 min 16 x 20,000 release units 49 RyRs 6 DHPR 1 second physiological simulation time NOTE: Colleagues running MPI spatial cardiac code use at
16 Ultrafast Monte Carlo Method CONDITION: Row sums at each matrix are zero
17 Ultrafast Monte Carlo Method (cont.) Single channel state: [x] with x=1..n, N = 5 Cluster state: [y1,y2,,yn] with y = number of i channels in state i-th Heterogeneous cluster state: [y1z1,y1z2,,y1zn,,ynz1,ynz2,,ynzn] with y i = number of channels in state i-th, z j = number of channels in state j-th
18 Ultrafast Monte Carlo Method (cont.) The method works Non-stationary state transition, i.e. functional transition rate Heterogeneous clusters where channels are modeled as Markov-chain Memory efficient using compact form representation of state transition matrix Lazy approach: only calculate probability of next state transition if needed
19 SR Ca 2+ Leak The leak of calcium out of the SR helps maintain calcium homeostasis. It balances the SR Ca 2+ -ATPase flux. It increases when Ca 2+ in the SR increases limiting SR loading. Increases in SR Ca 2+ can lead to larger Ca 2+ release events. In some conditions, such as heart failure, disease, and Ca 2+ overload, the leak has been suggested to increase the generation of cardiac arrhythmias.
20 Ca 2+ Leak Mechanisms Calcium release through the RyRs in the form of calcium sparks has been suggested to account for part of the SR Ca 2+ leak. Calcium spark rate increases with increasing SR load. However, there remains a certain amount of leak, called invisible leak that is not yet measured. Various sources have been suggested: Backflux through the SR Ca 2+ ATP ase Other ion channels IP 3 receptors Non-junctional RyRs
21 Ca 2+ Leak Mechanisms Backflux through the SR Ca2+ ATPase Backflux through the SR Ca 2+ ATPase has not been observed under physiological conditions. Only extreme experimental manipulation can do so. Other ion channels IP3 receptors Other ion channels have not been found. Calcium flux through IP 3 receptors has not been observed in adult myocytes where they are low in number <5% the number of RyRs. Non-junctional RyRs Non-junctional channels are very small in number less than 5% of the total number of RyRs. They see bulk myoplasmic calcium which does not reach the high levels needed to trigger calcium release in the diad. Flux through these is likely small. We propose and alternative mechanism to account for invisible leak
22 Model: Revised Sticky Cluster
23 Model Equations
24 SERCA Models
25 Model Solution RyR open state calculated using our Ultrafast Monte Carlo Method Fluxes calculated to determine derivatives Differential equations solved using a Euler Method Programmed in Fortran 90/CUDA Fortran (Portland Group Compiler) on a HP z800 Linux Workstation with NVIDIA Fermi 2050 GPUs
26 Ca 2+ Dependence of Open Probability
27 Calcium Spark Mechanism
28 Calcium Transient
29 Resting Ca 2+ spark behavior 20,000 CRU - 1% plotted
30 Individual Ca 2+ Spark
31 Spark and Quark Visualization
32 RyR Opening to Spark Transition
33 SR Ca 2+ Leak - Experiment (From Zima et al., 2010)
34 SR Ca 2+ Leak - Simulation
35 Leak Analysis
36 Effects of Phosphorylation
37 New Spatial Model Benchmarks Ultra-fast Method Monte Carlo On GPU No I/O I/O 3:09 hr 3:50 hr 20,000 release units 49 RyRs 6 DHPR > 4,000,000 grid elements 1 second physiological simulation time
38 Whole-cell Modeling The cell of size 100x20x18 mm 3 is modeled with a rectangular grid with a mesh of 0.2 mm At each grid point, it contains calcium in the myoplasm, calcium in the network SR.The T-tubule is assumed to be everywhere.
39 Whole-cell modelling (cont.) Euler method with forward difference in time and central difference in space Neumann boundary condition
40 Resting Myocyte Activity Experimental Spontaneous Calcium Sparks
41 Resting Myocyte Activity Simulated Spontaneous Calcium Sparks
42 Calcium Entrained Arrhythmias Experimental Calcium overload
43 Calcium Entrained Arrhythmias Simulated Calcium Overload
44 Conclusions Our Ultrafast Markov chain Monte Carlo method make stochastic simulation of calcium dynamics possible. Calcium release initiation occurs stochastically with the opening of one RyR that can trigger additional RyRs to open. One a critical number (~6 RyRs) opens, the remaining channels open causing a spark. Calcium release termination occurs through a combination of reduced SR Ca 2+ that results in reduced RyR opening, stochastic closure, and coupled gating. Calcium leak is comprised of spontaneous calcium sparks, and the opening of one or a few RyR channels in the release sites (invisible leak). Propagation between release sites depends upon calcium load, and release site placement.
45 Co-workers: Modeling Ca 2+ sparks and leak W. Jonathan Lederer BioMET / Johns Hopkins University Aristide Chikando Univ. Maryland Baltimore George (Blair) Williams George Mason U. / Univ. Maryland Baltimore Greg Smith William & Mary College Hoang Trong Minh Tuan George Mason University Eric A. Sobie Mt. Sinai School of Medicine This work was supported by the National Science Foundation, the National Institutes of Health, and the European Union 7 th Framework Program. Thanks to Steve Worley for use of his Pseudo-Random Number Generator GPU code.
46 Traditional Monte Carlo Algorithm (A) State Diagram with states X, Y, Z and transition probabilities p and q (B) A Uniform random number [0, 1] determines state transition.
47 Ultra-Fast Markov Chain Monte Carlo Algorithm Submitted for Publication and Patent Pending Properties: Fast Exact stochastic method Low memory usage Can be used for any Monte Carlo simulation.
48 Q-Matrix of Transition Probabilities Ryanodine Receptor M=2 states minimal model x = f(ca) : Ca-dependent: C O y = g(*) : Ca-independent: O C Single-channel rate-transition matrix Chapman-Kolmogorow equation
49 State Matrix A cluster of RyR State: (c 1, c 2 ) with c 1 = number of Closed RyRs Sˆ ( M N 1)! c 2 = number of Open RyRs N!( M 1)! E.g: N=5 RyRs Cluster rate-transition matrix AR(:,:), BR(:,:) of size 6x6 This reduces the state space from N M states increasing computational efficiency.
50 Transition Matrix for Cluster Exact simulation: In a small time-step, only a SINGLE channel can change state NOTE: rate out + rate in = 0 This allows use of vector-matrix algebra to perform Monte Carlo leveraging CPU/GPU design
51 Adaptive Time Step t min 0.1 P min The time step is chosen so that transitions only occur 10% of the time. P min is the most negative row sum in the transition matrix. Using the adaptive time step decreases simulation time by about 100x
52 Heterogenous Cluster Transition Matrix A heterogeneous cluster: e.g. release site with DHPR + RyR m = # of k-state RyR cluster states n = # of j-state DHPR cluster states
53 Motivation Complexity 50 2-state RyR: cluster of 51 states 7 6-state DHPR: cluster of 924 states Release site: 47,124 states Memory demand (double-precision): 16GB K matrix: Highly sparse Question - How to handle the computation with such sparse matrices in the GPU? Answer - We have developed a novel compact form representation and compact form Kronecker product.
54 Compact Form Matrix Representation Compact form for K: Use two separate matrices: Kcomp(:,:) = keep ONLY non-zero rate transition Kidx(:,:) = keep true column index of Kcomp(i,j) A Acomp Aidx x x
55 Compact Form Matrix Representation A homogeneous cluster: Full matrix Compact form Less number of conditional comparisons Less memory demand
56 Compact Form Matrix Representation A heterogeneous cluster, e.g. RyR + DHPR: Full matrix Compact form
57 Ultra-fast Monte Carlo Method Pros: Based upon matrix representation of the release units and their state space derived in part from the theory of solution of stochastic automata networks. Uses an adaptive time step. Only considers possible transitions from current state of each release units (scales as the number of possible transitions rather than the size of the state space) Reduces number of conditional comparisons. Amenable to parallelization.
58 Ultra-fast Monte Carlo Method Caveats: Reduces number of conditional comparisons by treating channels as identical so the number of channels in a particular configuration is counted. The transition from one state to another should be single agent dependent, e.g. either V or [Ca], but not both. m This does not allow for more than one event in a time step. In our simulations this only occurs ~6% of the time and this fraction can be reduced if needed.
59 Presentation Outline 1. Introduction to Excitation-Contraction Coupling 2. Calcium Sparks Experiment and Model 3. Ultrafast Monte Carlo Algorithm 4. GPU Implementation 5. Calcium-Entrained Arrhythmias 6. Conclusions
60 Computational Model of Resting Ca 2+ Dynamics RyR Minimal model with 2 states Release site: Each site has 50 RyR Whole-cell model 10,000 release sites
61 GPU Issues How GPU fit to our problem & algorithm Highly independent of release site computation Large amount of computation can be done in parallel State space is reused at every computational step Low memory demand makes it fit to the limited device memory (4GB in Tesla 1060, 3GB in Fermi) Need to minimize transfers between CPU and GPU due to Memory access latency Large amounts of I/O
Modeling 3-D Calcium Waves from Stochastic Calcium sparks in a Sarcomere Using COMSOL
Modeling 3-D Calcium Waves from Stochastic Calcium sparks in a Sarcomere Using COMSOL Zana Coulibaly 1, Leighton T. Izu 2 and Bradford E. Peercy 1 1 University of Maryland, Baltimore County 2 University
More informationA spatial Monte Carlo model of IP 3 -dependent calcium. oscillations. Evan Molinelli. Advisor: Gregry D. Smith. Department of Applied Science
A spatial Monte Carlo model of IP 3 -dependent calcium oscillations Evan Molinelli Advisor: Gregry D. Smith Department of Applied Science The College of William and Mary Williamsburg, VA 23187 4/14/06
More informationA model for diffusion waves of calcium ions in a heart cell is given by a system of coupled, time-dependent reaction-diffusion equations
Long-Time Simulation of Calcium Waves in a Heart Cell to Study the Effects of Calcium Release Flux Density and of Coefficients in the Pump and Leak Mechanisms on Self-Organizing Wave Behavior Zana Coulibaly,
More informationModeling. EC-Coupling and Contraction
Bioeng 6460 Electrophysiology and Bioelectricity Modeling of EC-Coupling and Contraction Frank B. Sachse fs@cvrti.utah.edu Overview Quiz Excitation-Contraction Coupling Anatomy Cross Bridge Binding Coupling
More informationParameter Sensitivity Analysis of Stochastic Models Provides Insights into Cardiac Calcium Sparks
1142 Biophysical Journal Volume 104 March 2013 1142 1150 Parameter Sensitivity Analysis of Stochastic Models Provides Insights into Cardiac Calcium Sparks Young-Seon Lee, Ona Z. Liu, Hyun Seok Hwang, Bjorn
More informationA calcium-induced calcium release mechanism mediated by calsequestrin
A calcium-induced calcium release mechanism mediated by calsequestrin Young-Seon Lee and James P. Keener Department of Biomedical Sciences, Cornell University Ithaca, New York 14853, USA Department of
More informationBasic mechanisms of arrhythmogenesis and antiarrhythmia
EHRA EDUCATIONAL REVIEW AND PREPARATORY COURSE ON INVASIVE CARDIAC ELECTROPHYSIOLOGY EUROPEAN HEART HOUSE, February 2011 Basic mechanisms of arrhythmogenesis and antiarrhythmia Antonio Zaza Università
More informationRole of coupled gating between cardiac ryanodine receptors in the genesis of triggered arrhythmias
Am J Physiol Heart Circ Physiol 297: H171 H180, 2009. First published May 8, 2009; doi:10.1152/ajpheart.00098.2009. Role of coupled gating between cardiac ryanodine receptors in the genesis of triggered
More informationLecture 13, 05 October 2004 Chapter 10, Muscle. Vertebrate Physiology ECOL 437 University of Arizona Fall instr: Kevin Bonine t.a.
Lecture 13, 05 October 2004 Chapter 10, Muscle Vertebrate Physiology ECOL 437 University of Arizona Fall 2004 instr: Kevin Bonine t.a.: Nate Swenson Vertebrate Physiology 437 18 1. Muscle A. Sarcomere
More informationExcitation-Contraction Coupling Gain and Cooperativity of the Cardiac Ryanodine Receptor: A Modeling Approach
Biophysical Journal Volume 89 November 2005 3017 3025 3017 Excitation-Contraction Coupling Gain and Cooperativity of the Cardiac Ryanodine Receptor: A Modeling Approach Kai Wang,* Yuhai Tu, y Wouter-Jan
More informationComputer Science Department Technical Report. University of California. Los Angeles, CA
Computer Science Department Technical Report University of California Los Angeles, CA 995-1596 COMPUTER SIMULATION OF THE JAFRI-WINSLOW ACTION POTENTIAL MODEL Sergey Y. Chernyavskiy November 1998 Boris
More informationThe Department of Electrical Engineering. nkrol Mentors: Dr. Mohammad Imtiaz, Dr. Jing Wang & Dr. In Soo Ahn
Bradley University The Department of Electrical Engineering nkrol Mentors: Dr. Mohammad Imtiaz, Dr. Jing Wang & Dr. In Soo Ahn AN ELECTRICAL ENGINEERING PERSPECTIVE OF THE HUMAN HEART This research project
More informationCellular Electrophysiology and Biophysics
BIOEN 6003 Cellular Electrophysiology and Biophysics Modeling of Force Development in Myocytes II Frank B. Sachse, University of Utah Overview Experimental Studies Sliding Filament Theory Group work Excitation-Contraction
More informationMuscle tissue. Types. Functions. Cardiac, Smooth, and Skeletal
Types Cardiac, Smooth, and Skeletal Functions movements posture and body position Support soft tissues Guard openings body temperature nutrient reserves Muscle tissue Special Characteristics of Muscle
More informationSpontaneous Calcium Release in Cardiac Myocytes: Store Overload and Electrical Dynamics
Spora: A Journal of Biomathematics Volume 1 Issue 1 Article 6 215 Spontaneous Calcium Release in Cardiac Myocytes: Store Overload and Electrical Dynamics Amanda M. Alexander Western Washington University,
More informationAnalysing cardiac excitation contraction coupling with mathematical models of local control
Progress in Biophysics & Molecular Biology 85 (2004) 141 162 Analysing cardiac excitation contraction coupling with mathematical models of local control Christian Soeller*, Mark B. Cannell Department of
More informationPHYSIOLOGY CHAPTER 9 MUSCLE TISSUE Fall 2016
PHYSIOLOGY CHAPTER 9 MUSCLE TISSUE Fall 2016 2 Chapter 9 Muscles and Muscle Tissue Overview of Muscle Tissue types of muscle: are all prefixes for muscle Contractility all muscles cells can Smooth & skeletal
More informationOur patient for the day...
Muscles Ch.12 Our patient for the day... Name: Eddy Age: Newborn Whole-body muscle contractions No relaxation Severe difficulty breathing due to inadequate relaxation of breathing muscles Diagnosed with
More informationCALIFORNIA STATE UNIVERSITY, NORTHRIDGE. Cardiac Cells. For the degree of Master of Science in Physics. Gonzalo Hernandez Hernandez
CALIFORNIA STATE UNIVERSITY, NORTHRIDGE Role of Network Connectivity and Fluctuations in the Nucleation of Calcium Waves in Cardiac Cells A thesis submitted in partial fulfillment of the requirements For
More informationUNIT 6 THE MUSCULAR SYSTEM
UNIT 6 THE MUSCULAR SYSTEM I. Functions of Muscular System A. Produces Movement Internal vs. External «locomotion & manipulation «circulate blood & maintain blood pressure «move fluids, food, baby B. Maintaining
More informationIntroduction to Physiology IV - Calcium Dynamics
Introduction to Physiology IV - Calcium Dynamics J. P. Keener Mathematics Department Introduction to Physiology IV - Calcium Dynamics p.1/26 Introduction Previous lectures emphasized the role of sodium
More informationCoombes, Stephen and Hinch, Robert and Timofeeva, Yulia (2004) Receptors, sparks and waves in a firediffuse-fire framework for calcium release.
Coombes, Stephen and Hinch, Robert and Timofeeva, Yulia (24) Receptors, sparks and waves in a firediffuse-fire framework for calcium release. Access from the University of Nottingham repository: http://eprints.nottingham.ac.uk/17/1/review.pdf
More informationAnatomy & Physiology Standards and Benchmarks
Anatomy & Standards and Standard 1: Understands and applies principles of scientific inquiry Power : Identifies questions and concepts that guide science investigations Uses technology and mathematics
More informationRole of the Na -Ca 2 Exchanger as an Alternative Trigger of CICR in Mammalian Cardiac Myocytes
Biophysical Journal Volume 82 March 2002 1483 1496 1483 Role of the Na -Ca 2 Exchanger as an Alternative Trigger of CICR in Mammalian Cardiac Myocytes Chunlei Han, Pasi Tavi, and Matti Weckström Department
More informationThe Journal of Physiology
J Physiol 590.18 (2012) pp 4403 4422 4403 Modelling cardiac calcium sparks in a three-dimensional reconstruction of a calcium release unit Johan Hake 1,2, Andrew G. Edwards 1,ZeyunYu 3, Peter M. Kekenes-Huskey
More informationEvolution of Cardiac Calcium Waves from Stochastic Calcium Sparks
Biophysical Journal Volume 80 January 2001 103 120 103 Evolution of Cardiac Calcium Waves from Stochastic Calcium Sparks Leighton T. Izu,* W. Gil Wier, and C. William Balke* Departments of *Medicine and
More informationلجنة الطب البشري رؤية تنير دروب تميزكم
1) Hyperpolarization phase of the action potential: a. is due to the opening of voltage-gated Cl channels. b. is due to prolonged opening of voltage-gated K + channels. c. is due to closure of the Na +
More informationBiology September 2015 Exam One FORM G KEY
Biology 251 17 September 2015 Exam One FORM G KEY PRINT YOUR NAME AND ID NUMBER in the space that is provided on the answer sheet, and then blacken the letter boxes below the corresponding letters of your
More informationBiology September 2015 Exam One FORM W KEY
Biology 251 17 September 2015 Exam One FORM W KEY PRINT YOUR NAME AND ID NUMBER in the space that is provided on the answer sheet, and then blacken the letter boxes below the corresponding letters of your
More informationCoupled dynamics of voltage and calcium in paced cardiac cells
Coupled dynamics of voltage and calcium in paced cardiac cells Yohannes Shiferaw, 1,2 Daisuke Sato, 1 and Alain Karma 1 1 Department of Physics and Center for Interdisciplinary Research on Complex Systems,
More informationApplication of Maxwell Equations to Human Body Modelling
Application of Maxwell Equations to Human Body Modelling Fumie Costen Room E3, E0c at Sackville Street Building, fc@cs.man.ac.uk November 6, 00 Fumie Costen Room E3, E0c at Sackville Street Building, Application
More information2013 NSF-CMACS Workshop on Atrial Fibrillation
2013 NSF-CMACS Workshop on A Atrial Fibrillation Flavio H. Fenton School of Physics Georgia Institute of Technology, Atlanta, GA and Max Planck Institute for Dynamics and Self-organization, Goettingen,
More informationIonic Charge Conservation and Long-Term Steady State in the Luo Rudy Dynamic Cell Model
3324 Biophysical Journal Volume 81 December 2001 3324 3331 Ionic Charge Conservation and Long-Term Steady State in the Luo Rudy Dynamic Cell Model Thomas J. Hund,* Jan P. Kucera,* Niels F. Otani,* and
More informationIn-Silico Modeling of Glycosylation Modulation Dynamics in K+ Ion Channels and Cardiac Signaling
In-Silico Modeling of Glycosylation Modulation Dynamics in K+ Ion Channels and Cardiac Signaling Dongping Du, Hui Yang Industrial and Management Systems Engineering Andrew R. Ednie, Eric S. Bennett Molecular
More informationCalcium dynamics, 7. Calcium pumps. Well mixed concentrations
Calcium dynamics, 7 Calcium is an important ion in the biochemistry of cells Is used a signal carrier, i.e. causes contraction of muscle cells Is toxic at high levels The concentration is regulated through
More informationIntroduction to numerical computations on the GPU
Introduction to numerical computations on the GPU Lucian Covaci http://lucian.covaci.org/cuda.pdf Tuesday 1 November 11 1 2 Outline: NVIDIA Tesla and Geforce video cards: architecture CUDA - C: programming
More informationPractical Combustion Kinetics with CUDA
Funded by: U.S. Department of Energy Vehicle Technologies Program Program Manager: Gurpreet Singh & Leo Breton Practical Combustion Kinetics with CUDA GPU Technology Conference March 20, 2015 Russell Whitesides
More informationInvariant Manifold Reductions for Markovian Ion Channel Dynamics
Invariant Manifold Reductions for Markovian Ion Channel Dynamics James P. Keener Department of Mathematics, University of Utah, Salt Lake City, UT 84112 June 4, 2008 Abstract We show that many Markov models
More informationA note on discontinuous rate functions for the gate variables in mathematical models of cardiac cells
Procedia Computer Science (2) (22) 6 945 95 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science ICCS 2 A note on discontinuous rate functions for
More informationInvariant Manifolds in Markovian Ion Channel Dynamics
Invariant Manifolds in Markovian Ion Channel Dynamics J. P. Keener Department of Mathematics University of Utah October 10, 2006 Ion channels are multiconformational proteins that open or close in a stochastic
More informationFinite Volume Model to Study Calcium Diffusion in Neuron Involving J RYR, J SERCA and J LEAK
JOURNAL OF COMPUTING, VOLUME 3, ISSUE 11, NOVEMB 011, ISSN 151-9617 WWW.JOURNALOFCOMPUTING.ORG 41 Finite Volume Model to Study lcium Diffusion in Neuron Involving J RYR, J SCA and J LEAK Amrita Tripathi
More informationFundamentals of Neurosciences. Smooth Muscle. Dr. Kumar Sambamurti 613-SEI; ;
Fundamentals of Neurosciences Smooth Muscle Dr. Kumar Sambamurti 613-SEI; 792-4315; sambak@musc.edu 1 Smooth Muscle Structure Cells much smaller than skeletal muscle (2-5µM diam, 100-400µM long) Single
More informationFrontiers in CardioVascular Biology
Frontiers in CardioVascular Biology Young Investigator Award Mitofusin 2 controls calcium transmission between the SR and mitochondria and regulates the bioenergetic feedback response in cardiac myocytes
More informationPhysiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES
Physiology Unit 2 MEMBRANE POTENTIALS and SYNAPSES In Physiology Today Ohm s Law I = V/R Ohm s law: the current through a conductor between two points is directly proportional to the voltage across the
More informationModel of Intracellular Calcium Cycling in Ventricular Myocytes
3666 Biophysical Journal Volume 85 December 2003 3666 3686 Model of Intracellular Calcium Cycling in Ventricular Myocytes Y. Shiferaw,* M. A. Watanabe,* A. Garfinkel, y J. N. Weiss, y and A. Karma* *Department
More informationRapid Report. Fundamental calcium release events revealed by two-photon excitation photolysis of caged calcium in guinea-pig cardiac myocytes
Keywords: Calcium release; Cardiac myocyte; Sarcoplasmic reticulum 7700 Journal of Physiology (1998), 508.3, pp. 801 809 801 Rapid Report Fundamental calcium release events revealed by two-photon excitation
More informationDynamical Systems for Biology - Excitability
Dynamical Systems for Biology - Excitability J. P. Keener Mathematics Department Dynamical Systems for Biology p.1/25 Examples of Excitable Media B-Z reagent CICR (Calcium Induced Calcium Release) Nerve
More informationEffectiveness of numerical techniques for calculating the quantity of calcium ion species during calcium sparksin heart muscle
Effectiveness of numerical techniques for calculating the quantity of calcium ion species during calcium sparksin heart muscle D. A. Lott Departments of Mathematics and Biotechnology, Delaware State University,
More informationEfficient implementation of the overlap operator on multi-gpus
Efficient implementation of the overlap operator on multi-gpus Andrei Alexandru Mike Lujan, Craig Pelissier, Ben Gamari, Frank Lee SAAHPC 2011 - University of Tennessee Outline Motivation Overlap operator
More informationAlgorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method
Algorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method Ilya B. Labutin A.A. Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 3, acad. Koptyug Ave., Novosibirsk
More informationReview. Watch the following video How Stuff Works - Cells video
Types of Cells Review Watch the following video How Stuff Works - Cells video Multi-cellular Organisms Living things can be as simple as single-celled organisms such as a Paramecium or an Amoeba that lives
More informationStochastic Modelling of Electron Transport on different HPC architectures
Stochastic Modelling of Electron Transport on different HPC architectures www.hp-see.eu E. Atanassov, T. Gurov, A. Karaivan ova Institute of Information and Communication Technologies Bulgarian Academy
More informationChem Lecture 9 Pumps and Channels Part 1
Chem 45 - Lecture 9 Pumps and Channels Part 1 Question of the Day: What two factors about a molecule influence the change in its free energy as it moves across a membrane? Membrane proteins function as
More informationIntroduction to electrophysiology. Dr. Tóth András
Introduction to electrophysiology Dr. Tóth András Topics Transmembran transport Donnan equilibrium Resting potential Ion channels Local and action potentials Intra- and extracellular propagation of the
More informationTwo Dimensional Coaxial Circular Elements in FEM to Study Calcium Diffusion in Neuron Cells
Applied Mathematical Sciences, Vol. 6, 2012, no. 10, 455-466 Two Dimensional Coaxial Circular Elements in FEM to Study Calcium Diffusion in Neuron Cells Amrita Tripathi and Neeru Adlakha Department of
More informationCELL BIOLOGY - CLUTCH CH. 9 - TRANSPORT ACROSS MEMBRANES.
!! www.clutchprep.com K + K + K + K + CELL BIOLOGY - CLUTCH CONCEPT: PRINCIPLES OF TRANSMEMBRANE TRANSPORT Membranes and Gradients Cells must be able to communicate across their membrane barriers to materials
More informationMembrane Potentials, Action Potentials, and Synaptic Transmission. Membrane Potential
Cl Cl - - + K + K+ K + K Cl - 2/2/15 Membrane Potentials, Action Potentials, and Synaptic Transmission Core Curriculum II Spring 2015 Membrane Potential Example 1: K +, Cl - equally permeant no charge
More informationIntroduction to electrophysiology 1. Dr. Tóth András
Introduction to electrophysiology 1. Dr. Tóth András Topics Transmembran transport Donnan equilibrium Resting potential Ion channels Local and action potentials Intra- and extracellular propagation of
More informationComputational Cell Biology
Computational Cell Biology Course book: Fall, Marland, Wagner, Tyson: Computational Cell Biology, 2002 Springer, ISBN 0-387-95369-8 (can be found in the main library, prize at amazon.com $59.95). Synopsis:
More informationUsing a CUDA-Accelerated PGAS Model on a GPU Cluster for Bioinformatics
Using a CUDA-Accelerated PGAS Model on a GPU Cluster for Bioinformatics Jorge González-Domínguez Parallel and Distributed Architectures Group Johannes Gutenberg University of Mainz, Germany j.gonzalez@uni-mainz.de
More informationMEMBRANE STRUCTURE. Lecture 9. Biology Department Concordia University. Dr. S. Azam BIOL 266/
MEMBRANE STRUCTURE Lecture 9 BIOL 266/4 2014-15 Dr. S. Azam Biology Department Concordia University RED BLOOD CELL MEMBRANE PROTEINS The Dynamic Nature of the Plasma Membrane SEM of human erythrocytes
More informationChapter 8 Calcium-Induced Calcium Release
Chapter 8 Calcium-Induced Calcium Release We started in Chap. 2 by assuming that the concentrations of the junctional sarcoplasmic reticulum (JSR) and the network sarcoplasmic reticulum (NSR) are identical
More informationPhysiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES
Physiology Unit 2 MEMBRANE POTENTIALS and SYNAPSES Neuron Communication Neurons are stimulated by receptors on dendrites and cell bodies (soma) Ligand gated ion channels GPCR s Neurons stimulate cells
More informationA Dynamic Model of Excitation-Contraction Coupling during Acidosis in Cardiac Ventricular Myocytes
3074 Biophysical Journal Volume 90 May 2006 3074 3090 A Dynamic Model of Excitation-Contraction Coupling during Acidosis in Cardiac Ventricular Myocytes Edmund J. Crampin and Nicolas P. Smith Bioengineering
More informationA Model Based Analysis of Steady-State versus Dynamic Elements in the Relationship between Calcium and Force
A Model Based Analysis of Steady-State versus Dynamic Elements in the Relationship between Calcium and Force Casey L. Overby, Sanjeev G. Shroff BACKGROUND Cardiac contraction and calcium. Intracellular
More informationLecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers
Lecture Notes 8C2 Inleiding Meten en Modelleren Cellular electrophysiology: modeling and simulation Nico Kuijpers nico.kuijpers@bf.unimaas.nl February 9, 2 2 8C2 Inleiding Meten en Modelleren Extracellular
More informationSlide 1. Slide 2. Slide 3. Muscles general information. Muscles - introduction. Microtubule Function
Slide 1 Muscles general information Vertebrates and many invertebrates have three main classes of muscle Skeletal muscle connect bones are are used for complex coordianted activities. Smooth muscles surround
More informationComputational Biology Course Descriptions 12-14
Computational Biology Course Descriptions 12-14 Course Number and Title INTRODUCTORY COURSES BIO 311C: Introductory Biology I BIO 311D: Introductory Biology II BIO 325: Genetics CH 301: Principles of Chemistry
More information1. Control a. In a negative feedback system, a is a device that measures the current value of the.
031804 Quiz Questions Week #1 1. Control a. In a negative feedback system, a is a device that measures the current value of the. b. Define an effector? c. What does it mean when a negative feedback system
More informationEffects of Ryanoids on Spontaneous and Depolarization-Evoked Calcium Release Events in Frog Muscle
Biophysical Journal Volume 87 July 2004 243 255 243 Effects of Ryanoids on Spontaneous and Depolarization-Evoked Calcium Release Events in Frog Muscle Chiu Shuen Hui,* Henry R. Besch Jr., y and Keshore
More informationA framework for detailed multiphase cloud modeling on HPC systems
Center for Information Services and High Performance Computing (ZIH) A framework for detailed multiphase cloud modeling on HPC systems ParCo 2009, 3. September 2009, ENS Lyon, France Matthias Lieber a,
More informationAccording to the diagram, which of the following is NOT true?
Instructions: Review Chapter 44 on muscular-skeletal systems and locomotion, and then complete the following Blackboard activity. This activity will introduce topics that will be covered in the next few
More informationPhysiology Coloring Book: Panels 29, 32, 33,
ELEC4623/ELEC9734: Semester 2, 2009 Dr Stephen Redmond School of Electrical Engineering & Telecommunications Email: s.redmond@unsw.edu.au Ph: 9385 6101 Rm: 458, ELECENG (G17) Physiology Coloring Book:
More informationAPPLICATION OF CUDA TECHNOLOGY FOR CALCULATION OF GROUND STATES OF FEW-BODY NUCLEI BY FEYNMAN'S CONTINUAL INTEGRALS METHOD
APPLICATION OF CUDA TECHNOLOGY FOR CALCULATION OF GROUND STATES OF FEW-BODY NUCLEI BY FEYNMAN'S CONTINUAL INTEGRALS METHOD M.A. Naumenko, V.V. Samarin Joint Institute for Nuclear Research, Dubna, Russia
More informationAccelerating Model Reduction of Large Linear Systems with Graphics Processors
Accelerating Model Reduction of Large Linear Systems with Graphics Processors P. Benner 1, P. Ezzatti 2, D. Kressner 3, E.S. Quintana-Ortí 4, Alfredo Remón 4 1 Max-Plank-Institute for Dynamics of Complex
More informationCardiac cell-cell Communication Part 1 Alonso P. Moreno D.Sc. CVRTI, Cardiology
Bioengineering 6003 Cellular Electrophysiology and Biophysics Cardiac cell-cell Communication Part 1 Alonso P. Moreno D.Sc. CVRTI, Cardiology moreno@cvrti.utah.edu November 2010 poster Physiological Relevance
More informationPerformance Evaluation of MPI on Weather and Hydrological Models
NCAR/RAL Performance Evaluation of MPI on Weather and Hydrological Models Alessandro Fanfarillo elfanfa@ucar.edu August 8th 2018 Cheyenne - NCAR Supercomputer Cheyenne is a 5.34-petaflops, high-performance
More informationTR A Comparison of the Performance of SaP::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems
TR-0-07 A Comparison of the Performance of ::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems Ang Li, Omkar Deshmukh, Radu Serban, Dan Negrut May, 0 Abstract ::GPU is a
More informationMembrane Potential. 1. Resting membrane potential (RMP): 2. Action Potential (AP):
Membrane Potential 1. Resting membrane potential (RMP): 2. Action Potential (AP): Resting Membrane Potential (RMP) It is the potential difference across the cell membrane. If an electrode of a voltmeter
More informationMuscle regulation and Actin Topics: Tropomyosin and Troponin, Actin Assembly, Actin-dependent Movement
1 Muscle regulation and Actin Topics: Tropomyosin and Troponin, Actin Assembly, Actin-dependent Movement In the last lecture, we saw that a repeating alternation between chemical (ATP hydrolysis) and vectorial
More informationLecture 3 13/11/2018
Lecture 3 13/11/2018 1 Plasma membrane ALL cells have a cell membrane made of proteins and lipids. protein channel Cell Membrane Layer 1 Layer 2 lipid bilayer protein pump Lipid bilayer allows water, carbon
More informationNonuniformity of Sarcomere Shortenings in the Isolated Rat Ventricular Myocyte
Japanese Journal of Physiology, 52, 371 381, 2002 Nonuniformity of Sarcomere Shortenings in the Isolated Rat Ventricular Myocyte Nobuaki SARAI, Yasuki KIHARA*, Toshiaki IZUMI*, Tamotsu MITSUIYE, Satoshi
More informationFaster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs
Faster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs Christopher P. Stone, Ph.D. Computational Science and Engineering, LLC Kyle Niemeyer, Ph.D. Oregon State University 2 Outline
More informationScalable and Power-Efficient Data Mining Kernels
Scalable and Power-Efficient Data Mining Kernels Alok Choudhary, John G. Searle Professor Dept. of Electrical Engineering and Computer Science and Professor, Kellogg School of Management Director of the
More informationIntroduction. The study of animal form and function is integrated by the common set of problems that all animals must solve.
Introduction The study of animal form and function is integrated by the common set of problems that all animals must solve. These include how to extract oxygen from the environment, how to nourish themselves,
More informationthe axons of the nerve meet with the muscle cell.
Steps to Contraction 1. A nerve impulse travels to the neuromuscular junction on a muscle cell. The neuromuscular junction is the point where the axons of the nerve meet with the muscle cell. 2. Ach is
More informationTransport of ions across plasma membranes
Transport of ions across plasma membranes Plasma Membranes of Excitable tissues Ref: Guyton, 13 th ed: pp: 61-71. 12 th ed: pp: 57-69. 11th ed: p57-71, Electrical properties of plasma membranes Part A:
More informationTISSUE. A) Types. (i)
MUSCLES & MUSCLE TISSUE I. OVERVIEW - Muscle ( little mouse ) - tissue designed to cause movementt thru contraction ( shortening ). A) Types - There are some SIMILARITIES between muscle types: (i) All
More informationThe Journal of Physiology
J Physiol 590.24 (2012) pp 6353 6361 6353 Do calcium waves propagate between cells and synchronize alternating calcium release in rat ventricular myocytes? Y. Li, D. A. Eisner and S. C. O Neill Unit of
More informationMembrane Protein Pumps
Membrane Protein Pumps Learning objectives You should be able to understand & discuss: Active transport-na + /K + ATPase ABC transporters Metabolite transport by lactose permease 1. Ion pumps: ATP-driven
More informationAPPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics
APPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics Due: Thursday, December 6, 2018 by 4:59 p.m. Submit as a PDF to Assignments on Canvas 1 Introduction Cardiac Arrhythmia, or irregular
More informationA CUDA Solver for Helmholtz Equation
Journal of Computational Information Systems 11: 24 (2015) 7805 7812 Available at http://www.jofcis.com A CUDA Solver for Helmholtz Equation Mingming REN 1,2,, Xiaoguang LIU 1,2, Gang WANG 1,2 1 College
More informationNeurons and Nervous Systems
34 Neurons and Nervous Systems Concept 34.1 Nervous Systems Consist of Neurons and Glia Nervous systems have two categories of cells: Neurons, or nerve cells, are excitable they generate and transmit electrical
More informationAn Efficient Calculation Algorithm in Continuous-Time Markov Analysis Using Large-Scale Numerical Calculation
International Mathematical Forum, Vol. 7, 2012, no. 10, 455-467 An Efficient Calculation Algorithm in Continuous-Time Markov Analysis Using Large-Scale Numerical Calculation Nobuko Kosugi and Koichi Suyama
More informationChapter 2 Basic Cardiac Electrophysiology: Excitable Membranes
Chapter Basic Cardiac Electrophysiology: Excitable Membranes Deborah A. Jaye, Yong-Fu Xiao, and Daniel C. Sigg Abstract Cardiomyocytes are excitable cells that have the ability to contract after excitation;
More informationFinite-Horizon Statistics for Markov chains
Analyzing FSDT Markov chains Friday, September 30, 2011 2:03 PM Simulating FSDT Markov chains, as we have said is very straightforward, either by using probability transition matrix or stochastic update
More informationAn Application of Reverse Quasi Steady State Approximation
IOSR Journal of Mathematics (IOSR-JM) ISSN: 78-578. Volume 4, Issue (Nov. - Dec. ), PP 37-43 An Application of Reverse Quasi Steady State Approximation over SCA Prashant Dwivedi Department of Mathematics,
More informationComputational Modeling of the Cardiovascular and Neuronal System
BIOEN 6900 Computational Modeling of the Cardiovascular and Neuronal System Modeling of Force Development in Myocytes Overview Recapitulation Modeling of Conduction Modeling of Force in Skeletal Muscle
More informationThe dynamics of Ca 2+ puffs: a study of instantaneously-coupled intracellular Ca 2+ channels. Author: Vien Nguyen
The dynamics of Ca 2+ puffs: a study of instantaneously-coupled intracellular Ca 2+ channels A thesis submitted in partial fulfillment for the degree of Bachelor of Science in Physics at the College of
More informationSolving PDEs with CUDA Jonathan Cohen
Solving PDEs with CUDA Jonathan Cohen jocohen@nvidia.com NVIDIA Research PDEs (Partial Differential Equations) Big topic Some common strategies Focus on one type of PDE in this talk Poisson Equation Linear
More information