caffeine-effects

ntumitolab · Nov 21, 2024 · 20a0841 · 20a0841
1 parent 830ebaf
commit 20a0841
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 8 deletions.
diff --git a/docs/_toc.yml b/docs/_toc.yml
@@ -11,3 +11,4 @@ chapters:
   - file: effects-bcl
   - file: iso-effects
   - file: ros-effects
+  - file: caffeine-effects
diff --git a/docs/caffeine-effects.jl b/docs/caffeine-effects.jl
@@ -1,6 +1,6 @@
 # # Caffeine Effects
 # Caffeine increase RyR opening sensitivity to luminal and subspace calcium
-# In this model, we decrease the mid saturation sub-SR calcium concentration for the opeing rate
+# In this model, we decrease the mid saturation sub-SR calcium concentration for the opening rate
 using ModelingToolkit
 using OrdinaryDiffEq
 using DiffEqCallbacks
@@ -20,25 +20,34 @@ function add_coffee_affect!(integrator)
 end
 
 @unpack Istim = sys
-callback = CallbackSet(build_stim_callbacks(Istim, stimend; period=1, starttime=stimstart), PresetTimeCallback(200.0, add_coffee_affect!))
+callback = build_stim_callbacks(Istim, stimend; period=1, starttime=stimstart)
+# Add caffeine at t = 200
+callback_caf = CallbackSet(build_stim_callbacks(Istim, stimend; period=1, starttime=stimstart), PresetTimeCallback(200.0, add_coffee_affect!))
 
 #---
 sol = solve(prob, alg; callback, abstol=1e-6, reltol=1e-6)
+sol_caf = solve(prob, alg; callback = callback_caf, abstol=1e-6, reltol=1e-6)
 
 #---
-plot(sol, idxs=sys.vm*1000, tspan=(199, 210), title="Action potential", lab=false)
+plot(sol, idxs=sys.vm*1000, title="Action potential", lab="Ctl")
+plot!(sol_caf, idxs=sys.vm*1000, lab="Caf", tspan=(199, 210), ylabel="Voltage (mV)")
 
 #---
-plot(sol, idxs=sys.Cai_sub_SR*1E6, tspan=(195, 210), title="Calcium transient", lab=false, ylabel="Subspace calcium (nM)")
+plot(sol, idxs=sys.Cai_sub_SR*1E6, title="Calcium transient (During caffeine addition)", lab="Ctl")
+plot!(sol_caf, idxs=sys.Cai_sub_SR*1E6, tspan=(195, 210), lab="Caf", ylabel="Subspace calcium (nM)")
 
 #---
-plot(sol, idxs=sys.Cai_sub_SR*1E6, tspan=(295, 300), title="Calcium transient", lab=false, ylabel="Subspace calcium (nM)")
+plot(sol, idxs=sys.Cai_sub_SR*1E6, title="Calcium transient (After caffeine addition)", lab="Ctl", ylabel="Subspace calcium (nM)")
+plot!(sol_caf, idxs=sys.Cai_sub_SR*1E6, tspan=(295, 300), lab="Caf")
 
 #---
-plot(sol, idxs=[sys.CaNSR*1E3, sys.CaJSR*1E3], tspan=(195, 210), title="SR Calcium", lab=["NSR" "JSR"], ylabel="SR calcium (μM)")
+plot(sol, idxs=sys.CaJSR*1E3, title="SR Calcium (During caffeine addition)", lab="Ctl", ylabel="SR calcium (μM)")
+plot!(sol_caf, idxs=sys.CaJSR*1E3, tspan=(195, 210), lab="Caf", ylims=(200, 850))
 
 #---
-plot(sol, idxs=sys.CaMKAct, title="Active CaMKII")
+plot(sol, idxs=sys.Jrel, title="RyR Ca flux", lab="Ctl")
+plot!(sol_caf, idxs=sys.Jrel, tspan=(195, 210), lab="Caf", ylabel="μM/ms")
 
 #---
-plot(sol, idxs=sys.CaMKAct, tspan=(195, 210), title="Active CaMKII")
+plot(sol, idxs=sys.CaMKAct, title="Active CaMKII", lab="Ctl")
+plot!(sol_caf, idxs=sys.CaMKAct, lab="Caf")