ROS effects#

using ModelingToolkit
using OrdinaryDiffEq, 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
@unpack Istim = sys
callback = build_stim_callbacks(Istim, stimend; period=1second, starttime=stimstart)
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(),)

Comparisons#

@time sol = solve(prob, alg; callback)
  3.634825 seconds (7.42 M allocations: 566.133 MiB, 3.72% gc time, 71.27% compilation 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]

ROS 0.1uM

prob2 = remake(prob, p=[sys.ROS => 0.1μM])
@time sol2 = solve(prob2, alg; callback)
  1.024565 seconds (65.99 k allocations: 212.910 MiB, 5.64% gc time)
retcode: Success
Interpolation: 3rd order Hermite
t: 4072-element Vector{Float64}:
      0.0
      0.019333123284849367
      0.09098129546758021
      0.20735582943983916
      0.45611910276862666
      1.0472562928406153
      2.0342667657223754
      3.451776641148464
      6.733792423912826
     12.45011708102323
      ⋮
 152613.47237072335
 157039.85798296845
 162443.11402596254
 168307.23424107043
 175796.8868001054
 183286.53935914036
 192001.04954351514
 201547.44731144857
 205000.0
u: 4072-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.0025849478335103772, 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.0025677078010499796, 829.9883445519497, 829.9990964774663, 0.007003014305234393, 0.9660430731013478, 0.2216542075240971, 0.0923749578279845, 0.0018699196262940572, 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.0025298083713907, 829.9738229154025, 829.9976227528975, 0.006946470397388756, 0.9660987611411453, 0.22177624785138977, 0.09230364339951613, 0.0018505337694469418, 0.009765228548607825, 0.26060507831375035  …  0.12115849269665879, 0.1211733598588045, 0.12119414515426227, 0.12122229800810433, 0.12125925929399752, 0.12130632552358436, 0.1213645020885198, -69.06429960994139, 13838.358559292059, 150952.75361533565]
 [0.0024748079782158953, 829.9509077855928, 829.9944193668179, 0.006833754172875444, 0.9661914345949096, 0.2219798989175973, 0.09218462776986637, 0.001817817441901466, 0.009763403076884485, 0.26047114911179464  …  0.12122454121844703, 0.12125242557826563, 0.12128628091507475, 0.12132672249013038, 0.12137426379919639, 0.12142927821043341, 0.1214919651325424, -69.31024285744344, 13838.34240329337, 150952.75694915713]
 [0.002410689396481131, 829.9205247545176, 829.9883533251192, 0.006675512062656741, 0.9663239343003682, 0.22227216398019892, 0.09201383579531924, 0.0017738156907532213, 0.009765046552800477, 0.2603648729749554  …  0.12132150441596079, 0.12135855488207915, 0.12140052118634272, 0.121447594710036, 0.12149988983034504, 0.12155743299711909, 0.12162015426876961, -69.64637448767434, 13838.319723120678, 150952.76215620188]
 [0.0023106557277343376, 829.8589858006806, 829.9689941817725, 0.006355039713677853, 0.9666286957812559, 0.22294828828976576, 0.09161908329520764, 0.0016855899094922301, 0.009786689200123853, 0.26041507667541997  …  0.1215043983379857, 0.12154722283384202, 0.12159304920072048, 0.12164180995533583, 0.12169340184085987, 0.12174768443532379, 0.12180447948880604, -70.35063788657429, 13838.2694938995, 150952.77598869594]
 [0.002228995040224516, 829.7704141396295, 829.923524981244, 0.005926150085676065, 0.9671554512230358, 0.22412560303273427, 0.09093422465365385, 0.0015688568862426025, 0.009873892913308327, 0.26104417804603025  …  0.1217127467387638, 0.12175152667167378, 0.12179171693004416, 0.12183318573340747, 0.12187579038290351, 0.12191937753328334, 0.12196378364408175, -71.35304509818073, 13838.18882346682, 150952.80529264064]
 ⋮
 [0.002044641089120985, 784.5152825814785, 784.8797814954723, 0.006281094863346002, 0.9997482220905666, 0.999573487470715, 0.0015084831575798863, 0.0024244782078248765, 0.0015662571468244466, 0.001113880853443376  …  0.13700214574472938, 0.1370021699420155, 0.13700218891537494, 0.13700220276943417, 0.13700221160466064, 0.1370022155175667, 0.13700221460090123, -70.44170248252901, 13613.396671443046, 151179.65359574076]
 [0.002030956786338835, 784.6341298264952, 784.9946574701265, 0.006303207012166781, 0.9997453925585061, 0.9996928143068631, 0.001519815302117068, 0.00243380795004504, 0.001572288603781846, 0.0011229124516400149  …  0.1367141063317371, 0.13671395397049182, 0.13671379929019276, 0.13671364230376487, 0.1367134830236133, 0.13671332146164883, 0.13671315762931188, -70.39002390703453, 13586.248015434106, 151206.8255624215]
 [0.002004126329169713, 784.4417909857252, 784.7965426204061, 0.006329124334760273, 0.9997420538594752, 0.9997310900052673, 0.0015331064056942013, 0.0024447473220491046, 0.0015793603823942518, 0.0011334961830436147  …  0.13633087422574827, 0.13633062208346844, 0.1363303693489091, 0.1363301159815109, 0.13632986194231905, 0.1363296071939048, 0.13632935170029142, -70.32968086166287, 13553.871270220992, 151239.2627412616]
 [0.0019701358288284065, 784.0428050346817, 784.390922701769, 0.006356435020101513, 0.9997385060753253, 0.9997378193199946, 0.001547148719505271, 0.0024562787196191265, 0.001586810276542772, 0.0011446602976260564  …  0.1359051208826579, 0.13590482607599602, 0.13590453151996326, 0.13590423714913138, 0.135903942900663, 0.13590364871418528, 0.13590335453167085, -70.26635877471438, 13519.636746410692, 151273.58001041217]
 [0.001925240872306776, 783.4219208380392, 783.7616079881939, 0.006390488794131708, 0.999734037783438, 0.999735698268376, 0.0015647314933619795, 0.0024706637930244153, 0.0015961088001366808, 0.001158713551158885  …  0.1353686250188678, 0.1353683146198888, 0.13536800493027448, 0.1353676958725128, 0.1353673873721632, 0.13536707935770667, 0.13536677176040357, -70.18777978062509, 13477.223912420175, 151316.1071603]
 [0.0018812836857426158, 782.7601830847561, 783.0917442171016, 0.0064238047181704476, 0.999729616561336, 0.9997315806453955, 0.0015820246944354779, 0.0024847436263076353, 0.0016052060822519388, 0.0011724751879512829  …  0.13484996593266108, 0.1348496557197288, 0.1348493464179446, 0.1348490379459314, 0.13484873022553506, 0.13484842318166848, 0.13484811674216168, -70.11130498246344, 13436.227663661131, 151357.21921926073]
 [0.001832478499265976, 781.9809820820941, 782.3035694084764, 0.006461735819121486, 0.9997245218279068, 0.9997265118964155, 0.0016018284513508651, 0.0025007813368116997, 0.0016155682137272425, 0.0011882783431881449  …  0.13427267940605297, 0.13427237536217165, 0.13427207235604702, 0.13427177030504303, 0.13427146912979765, 0.13427116875406286, 0.13427086910455377, -70.02471455854614, 13390.242488731152, 151403.33707317826]
 [0.001782168115163597, 781.1337975600096, 781.4471592386413, 0.006502300080339709, 0.9997190010675766, 0.9997209724198294, 0.0016231436816816081, 0.0025179413961018103, 0.0016266580320545324, 0.001205357338220176  …  0.13367276227687255, 0.13367246719126905, 0.13367217323132796, 0.13367188031448698, 0.13367158836145496, 0.13367129729605173, 0.13367100704505722, -69.93267010157932, 13341.90220076345, 151451.8186038869]
 [0.00176476913605272, 780.8305763558983, 781.1407513244353, 0.0065167190202319895, 0.999717020458502, 0.9997189839270085, 0.0016307550440619973, 0.0025240432778758115, 0.0016305998246819008, 0.0012114327900439692  …  0.13346390373438016, 0.13346361208661045, 0.1334633215899223, 0.13346303216196684, 0.13346274372365854, 0.13346245619901542, 0.13346216951500847, -69.90008938888167, 13324.926015454637, 151468.84470879182]

ROS 0.5uM

prob3 = remake(prob, p=[sys.ROS => 0.5μM])
@time sol3 = solve(prob3, alg; callback)
  1.008109 seconds (65.96 k allocations: 212.775 MiB, 3.88% gc time)
retcode: Success
Interpolation: 3rd order Hermite
t: 4070-element Vector{Float64}:
      0.0
      0.01933312328485172
      0.09098129546759129
      0.20735582943986441
      0.4561191027686822
      1.047256292840743
      2.0342667657226237
      3.451776641148885
      6.7337924239136475
     12.45011708102475
      ⋮
 148407.71053766683
 152963.08286658817
 157518.4551955095
 163268.1543188535
 169694.74092358912
 177941.12476481107
 186187.50860603302
 195845.67458090663
 205000.0
u: 4070-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.002598565959029538, 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.001878904244651442, 0.009769459049732313, 0.26078941009382306  …  0.12113000089551408, 0.12113000693481112, 0.12113004704056851, 0.12113027926811751, 0.12113143375871602, 0.12113624253259858, 0.12115248970210647, -68.81607005265553, 13838.374513299192, 150952.75062244057]
 [0.0025849478335103764, 829.9946506157202, 829.9996133960318, 0.007017234889558553, 0.9660195931792828, 0.22160283236040595, 0.09240497563401104, 0.0018766436937193184, 0.009768805322556195, 0.2607638707001999  …  0.12113010483012127, 0.12113038966346883, 0.12113131833063846, 0.12113403405679991, 0.12114109135582982, 0.12115723910142134, 0.12118951551963694, -68.84681913285227, 13838.372555328333, 150952.7509740238]
 [0.002567707801049977, 829.9883445519497, 829.9990964774663, 0.007003014305234387, 0.9660430731013478, 0.2216542075240971, 0.0923749578279845, 0.0018699196262940549, 0.00976754294689657, 0.2607122383588733  …  0.12113292110288301, 0.12113633873395534, 0.12114295934422187, 0.12115494411090462, 0.12117519662647318, 0.12120713813483536, 0.12125418995384044, -68.91207316777847, 13838.368384399068, 150952.7517373371]
 [0.0025298083713906926, 829.9738229154025, 829.9976227528975, 0.006946470397388735, 0.9660987611411453, 0.22177624785138977, 0.09230364339951613, 0.0018505337694469353, 0.009765228548607825, 0.26060507831375035  …  0.12115849269665878, 0.12117335985880448, 0.12119414515426227, 0.12122229800810433, 0.12125925929399756, 0.12130632552358439, 0.1213645020885198, -69.06429960994141, 13838.358559292059, 150952.75361533565]
 [0.002474807978215884, 829.9509077855928, 829.9944193668179, 0.006833754172875412, 0.9661914345949096, 0.2219798989175974, 0.09218462776986637, 0.0018178174419014564, 0.009763403076884485, 0.26047114911179464  …  0.12122454121844704, 0.12125242557826563, 0.12128628091507479, 0.12132672249013038, 0.12137426379919637, 0.12142927821043342, 0.12149196513254244, -69.31024285744353, 13838.34240329337, 150952.75694915713]
 [0.0024106893964811167, 829.9205247545176, 829.9883533251192, 0.006675512062656691, 0.9663239343003682, 0.22227216398019903, 0.09201383579531924, 0.0017738156907532087, 0.00976504655280048, 0.2603648729749554  …  0.1213215044159608, 0.12135855488207917, 0.12140052118634272, 0.12144759471003601, 0.12149988983034508, 0.12155743299711913, 0.12162015426876965, -69.64637448767446, 13838.319723120678, 150952.76215620188]
 [0.0023106557277343207, 829.8589858006806, 829.9689941817725, 0.006355039713677761, 0.966628695781256, 0.22294828828976596, 0.09161908329520757, 0.001685589909492206, 0.009786689200123865, 0.2604150766754201  …  0.12150439833798583, 0.12154722283384219, 0.12159304920072066, 0.12164180995533604, 0.12169340184086011, 0.12174768443532405, 0.1218044794888063, -70.35063788657448, 13838.2694938995, 150952.77598869594]
 [0.0022289950402245053, 829.7704141396295, 829.923524981244, 0.0059261500856759575, 0.967155451223036, 0.2241256030327346, 0.09093422465365371, 0.0015688568862425737, 0.009873892913308362, 0.2610441780460306  …  0.12171274673876394, 0.12175152667167392, 0.12179171693004433, 0.12183318573340762, 0.12187579038290368, 0.12191937753328352, 0.12196378364408192, -71.35304509818098, 13838.18882346682, 150952.80529264064]
 ⋮
 [0.0020427963561344326, 783.9341478970643, 784.3009800191974, 0.006258823893386904, 0.9997510563304198, 0.9992129260155501, 0.0014970539272434094, 0.002415084086567328, 0.0015601763447900638, 0.0011047375624272996  …  0.13722398335053676, 0.137224334219969, 0.1372246745950339, 0.1372250047480601, 0.13722532494055414, 0.1372256354237306, 0.1372259364390119, -70.49393530669363, 13639.748604866658, 151153.32454080155]
 [0.0020439940696006086, 784.5385372411694, 784.9027633010024, 0.006282881542686192, 0.9997479938239393, 0.9995887201790881, 0.0015093993782631398, 0.0024252317386472527, 0.001566742152498455, 0.0011146010144969364  …  0.13698085088117393, 0.13698085638263988, 0.13698085696606255, 0.13698085272638427, 0.13698084375476877, 0.13698083013878587, 0.13698081196258624, -70.4375201658923, 13611.229537457635, 151181.82102815836]
 [0.0020288997178260094, 784.6282590750978, 784.988302964019, 0.006305543651933786, 0.9997450927562128, 0.9996989883262268, 0.0015210127084674483, 0.002434794277042567, 0.0015729270318733041, 0.0011238554031907787  …  0.13668105895078114, 0.13668089389747587, 0.13668072673715345, 0.13668055747608374, 0.13668038612028185, 0.13668021267552047, 0.13668003714734162, -70.38457338896309, 13583.345473641277, 151209.73220898723]
 [0.001999525491567839, 784.3938616809376, 784.7476903413053, 0.006333010877620636, 0.9997415509465567, 0.999733135029437, 0.0015351020533405334, 0.0024463881677782386, 0.0015804204476755229, 0.0011350687412090702  …  0.13627109460155937, 0.1362708335942634, 0.13627057215943822, 0.1362703102515394, 0.13627004782682536, 0.13626978484326882, 0.13626952126047384, -70.32065303692754, 13548.9972914711, 151244.1474439959]
 [0.0019618334693176625, 783.9337455110241, 784.280286599746, 0.006362807821045566, 0.9997376738580532, 0.9997378278532045, 0.001550432003982062, 0.0024589703956224966, 0.0015885523918487254, 0.0011473069652382362  …  0.13580468448520872, 0.1358043846205855, 0.1358040851265144, 0.13580378593420947, 0.1358034869776085, 0.13580318819323955, 0.13580288952009778, -70.25162177440158, 13511.668616954526, 151281.56871601718]
 [0.0019124976932602198, 783.2345281910868, 783.5718515209029, 0.006400099177288416, 0.999732767644386, 0.9997345875789513, 0.0015697104260781222, 0.0024747248500681523, 0.001598732970677498, 0.0011626546444861892  …  0.13521808341373998, 0.13521777212898586, 0.1352174616270332, 0.13521715182876154, 0.13521684265818593, 0.13521653404230286, 0.13521622591094634, -70.16567912133726, 13465.343184619527, 151328.02088514858]
 [0.001864731944598966, 782.5008032181006, 782.8293177187606, 0.006436530497663171, 0.9997279147225948, 0.9997298997368163, 0.001588654816177273, 0.002490123587176036, 0.0016086831900824965, 0.0011777565888509798  …  0.1346545949888579, 0.13465428645344027, 0.13465397887882755, 0.13465367218299695, 0.13465336628717492, 0.134653061115678, 0.13465275659576353, -70.0821976642932, 13420.716792987181, 151372.77439442652]
 [0.001811818584950727, 781.6384050809147, 781.9572020208129, 0.006478196504189125, 0.9997222906497952, 0.9997242745687103, 0.0016104609350097014, 0.0025077436976412215, 0.0016200674880725236, 0.001195180636679787  …  0.13402702785261933, 0.13402672726335887, 0.13402642775179274, 0.13402612923518106, 0.1340258316340631, 0.1340255348720963, 0.13402523887590523, -69.98729456394896, 13370.522166863897, 151423.1146276122]
 [0.0017647781668277318, 780.8307127865006, 781.140889733984, 0.006516711396823504, 0.9997170216064906, 0.9997189847944175, 0.0016307509259814936, 0.002524040119234729, 0.0016305987411095327, 0.0012114429564518557  …  0.13346401863898513, 0.13346372700143774, 0.1334634365146879, 0.13346314709640952, 0.13346285866753982, 0.13346257115211874, 0.13346228447713912, -69.90010637539936, 13324.935163461343, 151468.83532441987]

Comparisons#

i = (sys.t / 1000, sys.CaMKAct * 100)
plot(sol, idxs=i, lab="ROS (-)")
plot!(sol2, idxs=i, lab="ROS 0.1uM")
plot!(sol3, idxs=i, lab="ROS 0.5uM")
plot!(xlabel="Time (s)", ylabel="Active fraction (%)", title="Simulation")
_images/5aa000ea4cf1aefdba3dee4796cc6d21c1bcfe77e213f38a6b45614d43403e4c.png
savefig("ros-camkii.pdf")
savefig("ros-camkii.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/ros-camkii.png"

Oxidized fraction

i = (sys.t / 1000, 100 * (sys.CaMKBOX + sys.CaMKPOX + sys.CaMKAOX + sys.CaMKOX ))
plot(sol, idxs=i, lab="ROS (-)")
plot!(sol2, idxs=i, lab="ROS 0.1uM")
plot!(sol3, idxs=i, lab="ROS 0.5uM")
plot!(xlabel="Time (s)", ylabel="Oxidized fraction (%)", title="Simulation")
_images/48b34e3dced39bbf75e72308a3d0f1a984ad94379e1d349e913841708d3dd8d8.png
savefig("ros-camkiiox.pdf")
savefig("ros-camkiiox.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/ros-camkiiox.png"

Autophosphorylated fraction

i = (sys.t / 1000, 100 * (sys.CaMKP + sys.CaMKA + sys.CaMKA2))
plot(sol, idxs=i, lab="ROS (-)")
plot!(sol2, idxs=i, lab="ROS 0.1uM")
plot!(sol3, idxs=i, lab="ROS 0.5uM")
plot!(xlabel="Time (s)", ylabel="Phosphorylated fraction (%)", title="Simulation")
_images/2011ab9d2a559ca41105e4af55ee05e9d00aa92475b8ad7942821f58c654a3dd.png
savefig("ros-camkiip.pdf")
savefig("ros-camkiip.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/ros-camkiip.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"])

ros50 = chemicaldf[!, "H2O2 50uM Mean"]
ros50_error = chemicaldf[!, "H2O2 50uM SD"] ./ sqrt.(chemicaldf[!, "H2O2 50uM N"])
ros200 = chemicaldf[!, "H2O2 200uM Mean"]
ros200_error = chemicaldf[!, "H2O2 200uM SD"] ./ sqrt.(chemicaldf[!, "H2O2 200uM N"])

plot(ts, ctl, yerr=ctl_error, lab="Control", color=:blue, markerstrokecolor=:blue)
plot!(ts, ros50, yerr=ros50_error, lab="H2O2 50uM", color=:red, markerstrokecolor=:red)
plot!(ts, ros200, yerr=ros200_error, lab="H2O2 200uM", color=:green, markerstrokecolor=:green)
plot!(xlabel="Time (s)", ylabel="CaMKII activity (A.U.)", title="Experiment")
_images/1e3a59a78a9924a02e0a0ba0a0f8eed32c4a86d56e918566566f6ddf4dcc17d4.png
savefig("ros-exp.pdf")
savefig("ros-exp.png")
"/home/runner/work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/.cache/docs/ros-exp.png"

This notebook was generated using Literate.jl.