Effects of isoproterenol#

using ModelingToolkit
using OrdinaryDiffEq, SteadyStateDiffEq, DiffEqCallbacks
using Plots
using CSV
using DataFrames
using Dates
using CaMKIIModel
using CaMKIIModel: second, μM
Plots.default(lw=1.5)

Setup model#

sys = build_neonatal_ecc_sys(simplify=true, reduce_iso=true, reduce_camk=true)
tend = 205second
prob = ODEProblem(sys, [], tend)
stimstart = 30second
stimend = 120second
alg = KenCarp47()
KenCarp47(; linsolve = nothing, nlsolve = OrdinaryDiffEqNonlinearSolve.NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}, Rational{Int64}}(1//100, 10, 1//5, 1//5, false, true, 0//1), precs = DEFAULT_PRECS, smooth_est = true, extrapolant = linear, controller = PI, autodiff = ADTypes.AutoForwardDiff(),)

Without isoproterenol#

@unpack Istim = sys
callback = build_stim_callbacks(Istim, stimend; period=1second, starttime=stimstart)
@time sol = solve(prob, alg; callback)
  1.041096 seconds (65.98 k allocations: 212.864 MiB, 3.59% gc time)
retcode: Success
Interpolation: 3rd order Hermite
t: 4071-element Vector{Float64}:
      0.0
      0.019333123284849308
      0.09098129546757994
      0.20735582943983855
      0.45611910276862533
      1.0472562928406122
      2.0342667657223696
      3.451776641148454
      6.733792423912806
     12.450117081023194
      ⋮
 153525.25012366346
 158165.8598019807
 163776.6411118618
 169611.6642461836
 176342.19489586214
 183684.1140466068
 191258.90946281631
 199427.37350661293
 205000.0
u: 4071-element Vector{Vector{Float64}}:
 [0.0026, 830.0, 830.0, 0.00702, 0.966, 0.22156, 0.09243, 0.00188, 0.00977, 0.26081  …  0.12113, 0.12113, 0.12113, 0.12113, 0.12113, 0.12113, 0.12113, -68.79268, 13838.37602, 150952.75035000002]
 [0.0025985659590295374, 829.9994975880693, 829.9999654017583, 0.007020208561323646, 0.9660018276402699, 0.2215639938807744, 0.09242766669618695, 0.001879826176081169, 0.009769882842069916, 0.2608055799112056  …  0.12113000000000405, 0.12113000000038378, 0.12113000001389383, 0.12113000037367234, 0.12113000811850358, 0.1211301421535992, 0.1211319185255175, -68.79732946186226, 13838.375711519551, 150952.7504077194]
 [0.002593310524681686, 829.9976422525165, 829.9998348635061, 0.00702008986027517, 0.9660085993318799, 0.2215787945110049, 0.09241901973901424, 0.0018789042446514422, 0.009769459049732313, 0.26078941009382306  …  0.12113000089551408, 0.12113000693481112, 0.12113004704056851, 0.12113027926811751, 0.12113143375871602, 0.12113624253259858, 0.12115248970210647, -68.81607005265552, 13838.374513299192, 150952.75062244057]
 [0.0025849478335103777, 829.9946506157202, 829.9996133960318, 0.007017234889558555, 0.9660195931792828, 0.22160283236040595, 0.09240497563401104, 0.0018766436937193197, 0.009768805322556195, 0.2607638707001999  …  0.12113010483012127, 0.12113038966346883, 0.12113131833063846, 0.12113403405679991, 0.12114109135582982, 0.12115723910142134, 0.12118951551963696, -68.84681913285226, 13838.372555328333, 150952.7509740238]
 [0.00256770780104998, 829.9883445519497, 829.9990964774663, 0.0070030143052343926, 0.9660430731013478, 0.2216542075240971, 0.0923749578279845, 0.0018699196262940568, 0.00976754294689657, 0.2607122383588733  …  0.12113292110288301, 0.12113633873395535, 0.12114295934422187, 0.12115494411090462, 0.12117519662647316, 0.12120713813483534, 0.12125418995384041, -68.91207316777846, 13838.368384399068, 150952.7517373371]
 [0.0025298083713907005, 829.9738229154025, 829.9976227528975, 0.006946470397388757, 0.9660987611411453, 0.22177624785138977, 0.09230364339951613, 0.001850533769446942, 0.009765228548607825, 0.26060507831375035  …  0.12115849269665879, 0.1211733598588045, 0.12119414515426227, 0.12122229800810434, 0.12125925929399752, 0.1213063255235844, 0.12136450208851983, -69.06429960994137, 13838.358559292059, 150952.75361533565]
 [0.0024748079782158975, 829.9509077855928, 829.9944193668179, 0.0068337541728754475, 0.9661914345949096, 0.2219798989175973, 0.09218462776986637, 0.0018178174419014683, 0.009763403076884485, 0.26047114911179464  …  0.12122454121844703, 0.12125242557826561, 0.12128628091507473, 0.12132672249013034, 0.12137426379919634, 0.12142927821043337, 0.12149196513254239, -69.31024285744343, 13838.34240329337, 150952.75694915713]
 [0.0024106893964811337, 829.9205247545176, 829.9883533251192, 0.006675512062656749, 0.9663239343003682, 0.22227216398019892, 0.09201383579531924, 0.001773815690753224, 0.009765046552800477, 0.2603648729749554  …  0.12132150441596078, 0.12135855488207913, 0.12140052118634269, 0.12144759471003595, 0.121499889830345, 0.12155743299711905, 0.12162015426876957, -69.64637448767432, 13838.319723120678, 150952.76215620188]
 [0.0023106557277343394, 829.8589858006806, 829.9689941817725, 0.006355039713677859, 0.9666286957812559, 0.22294828828976576, 0.09161908329520764, 0.0016855899094922325, 0.009786689200123853, 0.26041507667541997  …  0.12150439833798571, 0.12154722283384205, 0.12159304920072048, 0.1216418099553358, 0.12169340184085986, 0.12174768443532376, 0.12180447948880598, -70.35063788657428, 13838.2694938995, 150952.77598869594]
 [0.002228995040224517, 829.7704141396295, 829.923524981244, 0.005926150085676071, 0.9671554512230358, 0.22412560303273427, 0.09093422465365385, 0.0015688568862426045, 0.009873892913308327, 0.26104417804603025  …  0.12171274673876392, 0.1217515266716739, 0.12179171693004434, 0.12183318573340766, 0.12187579038290373, 0.1219193775332836, 0.121963783644082, -71.35304509818071, 13838.18882346682, 150952.80529264064]
 ⋮
 [0.0020427822499612955, 784.5705928361841, 784.9343641366844, 0.006285742961033888, 0.9997476285219717, 0.9996108819600831, 0.0015108660008011888, 0.0024264390093091466, 0.0015675234985916397, 0.0011157701814057516  …  0.13694598688419785, 0.13694596415399624, 0.13694593696562474, 0.13694590539945886, 0.13694586953267518, 0.13694582943940806, 0.13694578519089703, -70.43082440918127, 13607.749833006837, 151185.30233112013]
 [0.0020259750078751217, 784.6163365838529, 784.9757120846311, 0.006308692511330401, 0.9997446880935218, 0.9997062624414006, 0.0015226265701695494, 0.002436123054204679, 0.001573785829356736, 0.0011251534126126535  …  0.1366359698803179, 0.13663578906356108, 0.1366356064038893, 0.13663542189932926, 0.13663523554797946, 0.13663504734800627, 0.1366348572976408, -70.3772315171145, 13579.42650863126, 151213.6576494461]
 [0.001996646756046244, 784.362779202432, 784.7160349787207, 0.006335401714923246, 0.9997412411483091, 0.999734172546418, 0.0015363302060718079, 0.0024473974939953746, 0.0015810714524005405, 0.0011360309269089165  …  0.13623415561768498, 0.13623388970007247, 0.13623362344661866, 0.13623335680902646, 0.13623308974091072, 0.13623282219770474, 0.1362325541365715, -70.31510224724333, 13545.998515545509, 151247.1533714861]
 [0.00196232491987867, 783.9403317537941, 784.2869668994612, 0.006362431269318514, 0.9997377231755972, 0.9997378359259566, 0.001550237926945716, 0.0024588114887596946, 0.0015884495451699231, 0.0011471355458345597  …  0.1358106196440906, 0.13581032009374672, 0.13581002090542768, 0.1358097220105924, 0.13580942334341553, 0.13580912484065422, 0.1358088264415219, -70.2524921184015, 13512.139339623423, 151281.09701808463]
 [0.0019219897448604375, 783.3744083454476, 783.7134909839983, 0.006392942940491108, 0.9997337138514424, 0.999735422516465, 0.0015660024919756782, 0.002471700867704182, 0.0015967773659851202, 0.0011596790808541344  …  0.13533011366146125, 0.1353298029287261, 0.13532949292788785, 0.1353291835809063, 0.13532887481283415, 0.1353285665516654, 0.13532825872819246, -70.18213292915085, 13474.18574761773, 151319.15383306387]
 [0.0018789903164332137, 782.7245669481246, 783.0557065813091, 0.006425558662512325, 0.9997293823003752, 0.9997313503264378, 0.0015829377999889318, 0.002485484826637793, 0.0016056845198664101, 0.0011732157523120393  …  0.13482293289552036, 0.1348226228043163, 0.1348223136336112, 0.13482200530180374, 0.13482169773052602, 0.13482139084448463, 0.1348210845713118, -70.10728981816608, 13434.084114116084, 151359.36895294566]
 [0.0018365180658432869, 782.047149487092, 782.3704791996265, 0.006458544300802321, 0.9997249530520705, 0.999726943082142, 0.001600157158632046, 0.002499431602630893, 0.0016146973480228165, 0.0011869691561009399  …  0.13432065713674843, 0.13432035241861146, 0.1343200487300905, 0.13431974598855057, 0.1343194441146313, 0.1343191430320865, 0.13431884266763305, -70.03198079056111, 13394.082755688252, 151399.48570823882]
 [0.0017930492874439364, 781.3207468475036, 781.6361029305967, 0.0064933841889831265, 0.9997202209587296, 0.9997221973027168, 0.0016184463726731549, 0.0025141688480070073, 0.0016242200432758193, 0.0012015959628772666  …  0.13380304024010845, 0.13380274309271578, 0.13380244705424066, 0.13380215204203513, 0.13380185797672653, 0.13380156478205651, 0.13380127238473014, -69.95285215574584, 13352.452891572524, 151441.23701638903]
 [0.0017647619023658855, 780.830437838852, 781.1406117188316, 0.006516724021088335, 0.999717019859629, 0.9997189832869458, 0.0016307574425662278, 0.0025240455006058254, 0.0016306023325550296, 0.0012114432486847863  …  0.13346383404217976, 0.133463542406556, 0.13346325192122885, 0.13346296250384374, 0.13346267407530898, 0.1334623865596356, 0.13346209988378677, -69.90007814149098, 13324.920610871459, 151468.85017708497]
i = (sys.t / 1000, sys.vm)
plot(sol, idxs=i, tspan=(100second, 101second), title="Action potential", xlabel="Time (s)")
_images/b190ca3a982f55a403368609cdf6fe61a4533dca62ac65b6a50d883eda173a26.png
plot(sol, idxs=(sys.t / 1000, [sys.Cai_sub_SR, sys.Cai_sub_SL, sys.Cai_mean]), tspan=(100second, 101second), title="Calcium transient", xlabel="Time (s)", ylabel="Conc. (μM)", label=["Ca (SR)" "Ca (SL)" "Ca (avg)"])
_images/f2ae0291809546749ce7b6b2d2d27fb2c92bb674ae2856816276820116c7f0de.png
plot(sol, idxs=(sys.t / 1000, sys.CaMKAct * 100), title="Active CaMKII", label=false, ylabel="Active fraction (%)", xlabel="Time (s)")
_images/87094db655d7faa2b8cc50192c9697906db4e325dd22053b852ae79b24a361ea.png

0.1uM isoproterenol#

prob2 = remake(prob, p=[sys.ISO => 0.1μM])
sol2 = solve(prob2, alg; callback)
retcode: Success
Interpolation: 3rd order Hermite
t: 4138-element Vector{Float64}:
      0.0
      0.019332211071584192
      0.09097534627605569
      0.20734339686524408
      0.4561200538541287
      1.0549952720719427
      2.0687580593126613
      3.5324261824253345
      6.954473747659736
     12.881140201769806
      ⋮
 152270.482708526
 157174.07587050938
 162194.01343084572
 168138.29627423352
 174916.26987407374
 181898.84288967986
 189304.81198883272
 197588.1225791665
 205000.0
u: 4138-element Vector{Vector{Float64}}:
 [0.0026, 830.0, 830.0, 0.00702, 0.966, 0.22156, 0.09243, 0.00188, 0.00977, 0.26081  …  0.12113, 0.12113, 0.12113, 0.12113, 0.12113, 0.12113, 0.12113, -68.79268, 13838.37602, 150952.75035000002]
 [0.002598561776347308, 829.9994983969723, 830.0014907812235, 0.007020212947405008, 0.9660018275573545, 0.221563993694571, 0.0924276668078433, 0.0018798283012559313, 0.009769882789165586, 0.2608055790255411  …  0.12113000000000056, 0.12113000000008267, 0.12113000000364073, 0.1211300001166239, 0.12113000299821446, 0.12113006189494423, 0.1211309809540452, -68.79688468304944, 13838.37575462767, 150952.75040560117]
 [0.0025932333906583892, 829.9976597572193, 830.0069923279007, 0.007020182507782679, 0.9660085988422271, 0.22157879333123556, 0.0924190204915213, 0.0018789474208576141, 0.009769457815727975, 0.26078938725640693  …  0.12113000020174947, 0.12113000185234848, 0.12113001474369577, 0.12113010200803158, 0.12113060640007198, 0.12113303174274752, 0.12114238287990785, -68.81397720473576, 13838.374716348011, 150952.75061240065]
 [0.002584621446261353, 829.9947407323419, 830.0158577065781, 0.007017680221388922, 0.9660195923730465, 0.22160283004278358, 0.09240497731310347, 0.0018768406238241547, 0.009768798977675412, 0.26076374828883014  …  0.12113003007199945, 0.12113012844945044, 0.12113049616083325, 0.12113171992467786, 0.1211353044921302, 0.12114443083154028, 0.12116438512395096, -68.84205163448622, 13838.373018773336, 150952.75095089723]
 [0.0025665030871097365, 829.9887756080547, 830.0345419613458, 0.007004853088813848, 0.9660430748851897, 0.22165420888269807, 0.09237495857236798, 0.0018706615190431774, 0.009767513260722714, 0.26071163727178975  …  0.12113109689132587, 0.12113262872712134, 0.12113588865795849, 0.12114231020892302, 0.12115399268330548, 0.1211735958345763, 0.12120393253792845, -68.90160259182096, 13838.369406865842, 150952.7516853767]
 [0.002524878616001687, 829.9758982053276, 830.07814410011, 0.006952675605532075, 0.9660994972023292, 0.22177785078777174, 0.09230271445722965, 0.0018527422088444536, 0.009765055027054551, 0.26060065813180305  …  0.12114357362996761, 0.12115176909538257, 0.12116376521507281, 0.12118067468437907, 0.12120363670655872, 0.12123369231137104, 0.12127164038451914, -69.0421184107563, 13838.360815589964, 150952.7535137264]
 [0.00246109996477572, 829.9585434198331, 830.1479502978362, 0.0068472036120339344, 0.9661946906443825, 0.22198703092009986, 0.09218048699222062, 0.0018222096890642774, 0.009762833375314839, 0.26045565430474177  …  0.12118006281992985, 0.12119663497761278, 0.12121718646328616, 0.12124217405786958, 0.12127196546348169, 0.12130680044052823, 0.12134675555621606, -69.27157622080777, 13838.346581408403, 150952.75679928518]
 [0.002382758683301369, 829.9422438990155, 830.2409658671871, 0.0066990037339929405, 0.9663315006590419, 0.22228882452599166, 0.09200417115231563, 0.0017811002406865244, 0.009763716676849107, 0.26032881924894974  …  0.12123448815456439, 0.12125657408666768, 0.1212817846915382, 0.12131020885887851, 0.12134186045923348, 0.12137666702212292, 0.12141446087873138, -69.58534388478576, 13838.326681500257, 150952.76194966168]
 [0.002249831541002491, 829.9364424528658, 830.4298664467993, 0.0063983355059790605, 0.9666490840575152, 0.22299376432344578, 0.09159274697262061, 0.0016985405358012685, 0.009783000857592001, 0.2603166532077402  …  0.12132652897086686, 0.12134992079665141, 0.12137464276885106, 0.12140056903109953, 0.12142753330789122, 0.12145532836849124, 0.12148370617698084, -70.2425024040955, 13838.282994292276, 150952.77573103318]
 [0.0021239182944550387, 829.9982399691283, 830.6928613527689, 0.0060038545810325245, 0.9671939916127574, 0.22421393887512645, 0.09088324181611296, 0.0015916044910949695, 0.009862701804147046, 0.2607890677215021  …  0.12137592664517363, 0.12139282572293232, 0.1214091612631606, 0.12142480441004644, 0.12143960641561434, 0.1214533998189489, 0.12146599978828718, -71.1550700755306, 13838.215505646505, 150952.80458617522]
 ⋮
 [0.0073146875004273795, 921.3111350424876, 922.8372442673206, 0.00621632953768987, 0.9997563697872119, 0.9995674489205083, 0.0014756873624894675, 0.0023971671168631684, 0.0015486035494228147, 0.001087827177827683  …  0.12690731956341642, 0.1269071553885708, 0.12690698535411543, 0.12690680954319294, 0.1269066280356576, 0.12690644090823658, 0.1269062482346817, -70.59411328190205, 13855.447607237578, 150927.44243962405]
 [0.00722290818564516, 921.4901697334809, 922.9973924318351, 0.0062433363972850256, 0.9997529999609079, 0.9997032071128521, 0.0014892812418353196, 0.0024085528978078563, 0.0015559585586829552, 0.001098556090082494  …  0.12644089018468604, 0.12644072850764165, 0.12644056118440006, 0.1264403882943622, 0.1264402099137888, 0.1264400261159544, 0.12643983697129246, -70.53036749029351, 13820.442136600705, 150962.49883417104]
 [0.007132108627847399, 921.6620339611582, 923.150557185851, 0.006270727254031783, 0.9997495493507649, 0.9997381113366643, 0.0015031331537414113, 0.002420104889598678, 0.00156342418330898, 0.0011095708614125453  …  0.12597799318211722, 0.12597783414482552, 0.12597766966765508, 0.1259774998264211, 0.12597732469394182, 0.12597714434018578, 0.12597695883240992, -70.46599490741474, 13785.406967581663, 150997.58514476087]
 [0.007028583082949953, 921.8514954087509, 923.3186796523768, 0.006302819980298003, 0.9997454633663905, 0.9997449540510915, 0.0015194458293461555, 0.002433645223031485, 0.00157217239732809, 0.0011225355991132883  …  0.12544833103095102, 0.1254481749666311, 0.12544801370186628, 0.1254478473083215, 0.125447675854829, 0.12544749940752703, 0.12544731802999107, -70.39092616937515, 13744.94321311874, 151038.10815449725]
 [0.006915552382464886, 922.0500067523041, 923.4938704460396, 0.006338953528049358, 0.9997408078707524, 0.9997426766644008, 0.0015379203360093734, 0.00244889711467549, 0.0015820251283820395, 0.0011372511426259075  …  0.1248678046342907, 0.12486765205464238, 0.12486749452901216, 0.12486733212471628, 0.12486716490641103, 0.12486699293622192, 0.12486681627386738, -70.30685771446443, 13700.121315694741, 151082.99594042663]
 [0.006804459843946287, 922.2361636601723, 923.6570847387662, 0.006375667093771092, 0.9997360201470834, 0.9997382835800523, 0.0015568059542802523, 0.0024644027648597347, 0.001592046998596164, 0.0011523065215193354  …  0.12429479326782757, 0.12429464424155347, 0.12429449052612665, 0.12429433218447485, 0.12429416927703935, 0.12429400186189714, 0.1242938299948755, -70.22192860931602, 13655.37030331296, 151127.81298682958]
 [0.0066922417332850185, 922.4148307507354, 923.8125540268298, 0.00641403017111655, 0.9997309526599988, 0.9997332832912299, 0.0015766664531288673, 0.002480613061087455, 0.0016025195303765618, 0.0011680748522190849  …  0.12371356472390317, 0.12371341943311376, 0.1237132697054624, 0.12371311559962239, 0.12371295717194818, 0.12371279447658921, 0.1237126275655968, -70.13370023408227, 13609.434641207708, 151173.81666220704]
 [0.006573203256253084, 922.5934983659067, 923.9665888758834, 0.0064562278992018616, 0.9997253010148252, 0.9997276295977502, 0.001598658841927146, 0.0024984529975030883, 0.001614047794745646, 0.0011856577970248887  …  0.12309422660091429, 0.12309408550324269, 0.12309394023740704, 0.12309379085759026, 0.12309363741583354, 0.12309347996214129, 0.12309331854458028, -70.03725612195532, 13559.858385460602, 151223.46657862776]
 [0.0064721693428980796, 922.7358937568729, 924.088057654229, 0.006493344733516329, 0.9997202628390558, 0.9997225836604192, 0.001618131395065923, 0.0025141532138503696, 0.0016241926698836691, 0.0012012301821853206  …  0.1225663038391262, 0.12256616646164024, 0.12256602513991986, 0.12256587992446184, 0.12256573086376693, 0.1225655780044376, 0.12256542139127033, -69.95294083205617, 13517.05790484642, 151266.3308482879]
plot(sol2, idxs=(sys.t / 1000, sys.vm), tspan=(100second, 101second), title="Action potential", xlabel="Time (s)")
_images/82abf60cd4110d279b8c0bea881452ff9428320206a2edecd1ed933f435d6b66.png
plot(sol2, idxs=(sys.t / 1000, [sys.Cai_sub_SR, sys.Cai_sub_SL, sys.Cai_mean]), tspan=(100second, 101second), title="Calcium transcient", xlabel="Time (s)", ylabel="Conc. (μM)", label=["Ca (SR)" "Ca (SL)" "Ca (avg)"])
_images/8396ed7435d858b7c06d74d688d3d40555e92d49bc1951be46750f3b184ffc34.png
plot(sol2, idxs=(sys.t / 1000, sys.CaMKAct * 100), title="Active CaMKII", label=false, ylabel="Active fraction (%)", xlabel="Time (s)")
_images/c0592c7164c01edbb39075d389ad1b9dcc7cd0912071d1c3ecf4d3a1c344c52f.png

Comparison#

i = (sys.t / 1000, sys.Cai_mean)
tspan = (100second, 101second)
plot(sol, idxs=i, title="Calcium transcient", lab="ISO (-)"; tspan)
plot!(sol2, idxs=i, lab="ISO (0.1uM)", xlabel="Time (s)", ylabel="Concentration (μM)"; tspan)
_images/f8881b20b73c069d36d29de5badc77a116d6f333379ba08b50ec45bf33dd92fc.png

Compare the numbers

ca_ctl = sol(tspan[1]:1:tspan[2], idxs=sys.Cai_mean)
println(extrema(ca_ctl))
(0.23505056446792008, 0.5781214777886791)
ca_iso = sol2(tspan[1]:1:tspan[2], idxs=sys.Cai_mean)
println(extrema(ca_iso))
(0.19033626526983588, 0.711093945093875)
savefig("iso-caT.pdf")
savefig("iso-caT.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/iso-caT.png"
i = (sys.t / 1000, sys.CaMKAct * 100)
plot(sol, idxs=i, title="Active CaMKII", lab="ISO (-)")
plot!(sol2, idxs=i, lab="ISO (0.1uM)", ylabel="Active fraction (%)", xlabel="Time (s)")
_images/4ab3bc3617e3a5bce825e65dabbea725d7cd155fffb3e30f9e506bab71e46295.png
savefig("iso-camkact.pdf")
savefig("iso-camkact.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/iso-camkact.png"
i = (sys.t / 1000, sys.vm)
tspan = (100second, 101second)
plot(sol, idxs=i, title="Action potential", lab="ISO (-)"; tspan)
plot!(sol2, idxs=i, lab="ISO (0.1uM)", xlabel="Time (ms)", ylabel="Voltage (mV)"; tspan)
_images/7e23b9b339879ceacc079ccb7b1e589f057c05dac0c9a39fde163fa9fade18ea.png

Experimental data#

chemicaldf = CSV.read(joinpath(@__DIR__, "data/CaMKAR-chemical.csv"), DataFrame)
ts = Dates.value.(chemicaldf[!, "Time"]) ./ 10^9
ctl = chemicaldf[!, "Ctrl Mean"]
ctl_error = chemicaldf[!, "Ctrl SD"] ./ sqrt.(chemicaldf[!, "Ctrl N"])

iso = chemicaldf[!, "isoproterenol 100nM Mean"]
iso_error = chemicaldf[!, "isoproterenol 100nM SD"] ./ sqrt.(chemicaldf[!, "isoproterenol 100nM N"])

plot(ts, ctl, yerr=ctl_error, lab="Control", color=:blue, markerstrokecolor=:blue)
plot!(ts, iso, yerr=iso_error, lab="ISO 100nM", color=:red, markerstrokecolor=:red)
plot!(xlabel="Time (sec.)", ylabel="CaMKII activity (A.U.)")
_images/e151e9cc176acf7e0374a623a6bdf51ee537fba88b436053999e3ebc4d9d0996.png
savefig("iso-exp.pdf")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/iso-exp.pdf"

This notebook was generated using Literate.jl.