/Bioelectrics/Bidomain/jrw-1998.xml

http://github.com/xyan075/examples · XML · 5333 lines · 4738 code · 290 blank · 305 comment · 0 complexity · 37f7f29a040fabee36fb7c2b8debc8e2 MD5 · raw file

Large files are truncated click here to view the full file

  1. <?xml version="1.0" encoding="iso-8859-1"?>
  2. <!-- FILE : jafri_rice_winslow_model_1998.xml
  3. CREATED : September 2001
  4. LAST MODIFIED : 30th July 2003
  5. AUTHOR : Catherine Lloyd
  6. Department of Engineering Science
  7. The University of Auckland
  8. MODEL STATUS : This model conforms to the CellML 1.0 Specification released on
  9. 10th August 2001, and the CellML Metadata 1.0 Specification released on 16th
  10. January, 2002.
  11. DESCRIPTION : This file contains a CellML description of the mammalian
  12. ventricular action potential based on the Jafri-Rice-Winslow model, 1998. This
  13. model is a development of the LR-II model. In particular, it makes an accurate
  14. model of the membrane currents and adds a more sophisticated model of calcium
  15. ion handling.
  16. CHANGES:
  17. 19/10/2001 - CML - Removed document type definition as this is declared as
  18. optional according to the W3C recommendation.
  19. 24/10/2001 - CML - Made changes to some of the metadata, bringing them up to
  20. date with the most recent working draft (26th September) of
  21. the Metadata specification.
  22. 07/12/2001 - CML - Changed tau_y_calculation after checking mathml using the
  23. validator.
  24. 04/01/2002 - CML - Altered some of the connections.
  25. 21/01/2002 - AAC - Updated metadata to conform to the 16/1/02 CellML Metadata
  26. 1.0 Specification.
  27. 25/02/2002 - CML - Corrected several equations.
  28. 28/02/2002 - CML - Corrected several equations, variable units and their
  29. initial values.
  30. 06/05/2002 - CML - Added some initial values.
  31. 22/07/2002 - CML - Added more metadata.
  32. 09/04/2003 - AAC - Added publication date information.
  33. 04/06/2003 - CML - Fixed MathML in a few components.
  34. 30/07/2003 - CML - Altered a few equations.
  35. -->
  36. <model
  37. name="jafri_rice_winslow_model_1998"
  38. cmeta:id="jafri_rice_winslow_model_1998"
  39. xmlns="http://www.cellml.org/cellml/1.0#"
  40. xmlns:cellml="http://www.cellml.org/cellml/1.0#"
  41. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  42. <rdf:RDF
  43. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  44. xmlns:cmeta="http://www.cellml.org/metadata/1.0#"
  45. xmlns:bqs="http://www.cellml.org/bqs/1.0#"
  46. xmlns:dc="http://purl.org/dc/elements/1.1/"
  47. xmlns:dcterms="http://purl.org/dc/terms/"
  48. xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
  49. <!--
  50. The following RDF block contains metadata that applies to this document
  51. as a whole, as indicated by the empty about attribute on the
  52. <rdf:Description> element.
  53. -->
  54. <rdf:Description rdf:about="">
  55. <!--
  56. The Model Builder Metadata. The Dublin Core "creator" element is used
  57. to indicate the person who translated the model into CellML.
  58. -->
  59. <dc:creator rdf:parseType="Resource">
  60. <vCard:N rdf:parseType="Resource">
  61. <vCard:Family>Lloyd</vCard:Family>
  62. <vCard:Given>Catherine</vCard:Given>
  63. <vCard:Other>May</vCard:Other>
  64. </vCard:N>
  65. <vCard:EMAIL rdf:parseType="Resource">
  66. <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  67. <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet" />
  68. </vCard:EMAIL>
  69. <vCard:ORG rdf:parseType="Resource">
  70. <vCard:Orgname>The University of Auckland</vCard:Orgname>
  71. <vCard:Orgunit>The Bioengineering Research Group</vCard:Orgunit>
  72. </vCard:ORG>
  73. </dc:creator>
  74. <!--
  75. The Creation Date metadata. This is the date on which the model
  76. was translated into CellML.
  77. -->
  78. <dcterms:created rdf:parseType="Resource">
  79. <dcterms:W3CDTF>2001-09-24</dcterms:W3CDTF>
  80. </dcterms:created>
  81. <!--
  82. The Modification History metadata. This lists the changes that have been
  83. made to the document, who made the changes, and when they were made.
  84. -->
  85. <cmeta:modification rdf:parseType="Resource">
  86. <rdf:value>
  87. Corrected equations.
  88. </rdf:value>
  89. <cmeta:modifier rdf:parseType="Resource">
  90. <vCard:N rdf:parseType="Resource">
  91. <vCard:Family>Lloyd</vCard:Family>
  92. <vCard:Given>Catherine</vCard:Given>
  93. <vCard:Other>May</vCard:Other>
  94. </vCard:N>
  95. </cmeta:modifier>
  96. <dcterms:modified rdf:parseType="Resource">
  97. <dcterms:W3CDTF>2003-07-30</dcterms:W3CDTF>
  98. </dcterms:modified>
  99. </cmeta:modification>
  100. <cmeta:modification rdf:parseType="Resource">
  101. <rdf:value>
  102. Corrected equations: alpha_j_calculation and beta_j_calculation in
  103. fast_sodium_current_j_gate, alpha_X_calculation and beta_X_calculation in time_dependent_potassium_current_X_gate, and f_NaK_calculation and
  104. i_NaK_calculation in Ca_release_current_from_JSR.
  105. </rdf:value>
  106. <cmeta:modifier rdf:parseType="Resource">
  107. <vCard:N rdf:parseType="Resource">
  108. <vCard:Family>Lloyd</vCard:Family>
  109. <vCard:Given>Catherine</vCard:Given>
  110. <vCard:Other>May</vCard:Other>
  111. </vCard:N>
  112. </cmeta:modifier>
  113. <dcterms:modified rdf:parseType="Resource">
  114. <dcterms:W3CDTF>2003-06-04</dcterms:W3CDTF>
  115. </dcterms:modified>
  116. </cmeta:modification>
  117. <cmeta:modification rdf:parseType="Resource">
  118. <rdf:value>
  119. Added some initial values from Penny Noble's documentation.
  120. </rdf:value>
  121. <cmeta:modifier rdf:parseType="Resource">
  122. <vCard:N rdf:parseType="Resource">
  123. <vCard:Family>Lloyd</vCard:Family>
  124. <vCard:Given>Catherine</vCard:Given>
  125. <vCard:Other>May</vCard:Other>
  126. </vCard:N>
  127. </cmeta:modifier>
  128. <dcterms:modified rdf:parseType="Resource">
  129. <dcterms:W3CDTF>2002-05-06</dcterms:W3CDTF>
  130. </dcterms:modified>
  131. </cmeta:modification>
  132. <cmeta:modification rdf:parseType="Resource">
  133. <rdf:value>
  134. Corrected several equations, variable units and their initial values.
  135. </rdf:value>
  136. <cmeta:modifier rdf:parseType="Resource">
  137. <vCard:N rdf:parseType="Resource">
  138. <vCard:Family>Lloyd</vCard:Family>
  139. <vCard:Given>Catherine</vCard:Given>
  140. <vCard:Other>May</vCard:Other>
  141. </vCard:N>
  142. </cmeta:modifier>
  143. <dcterms:modified rdf:parseType="Resource">
  144. <dcterms:W3CDTF>2002-02-28</dcterms:W3CDTF>
  145. </dcterms:modified>
  146. </cmeta:modification>
  147. <cmeta:modification rdf:parseType="Resource">
  148. <rdf:value>
  149. Corrected several equations.
  150. </rdf:value>
  151. <cmeta:modifier rdf:parseType="Resource">
  152. <vCard:N rdf:parseType="Resource">
  153. <vCard:Family>Lloyd</vCard:Family>
  154. <vCard:Given>Catherine</vCard:Given>
  155. <vCard:Other>May</vCard:Other>
  156. </vCard:N>
  157. </cmeta:modifier>
  158. <dcterms:modified rdf:parseType="Resource">
  159. <dcterms:W3CDTF>2002-02-25</dcterms:W3CDTF>
  160. </dcterms:modified>
  161. </cmeta:modification>
  162. <cmeta:modification rdf:parseType="Resource">
  163. <rdf:value>
  164. Altered some of the connections.
  165. </rdf:value>
  166. <cmeta:modifier rdf:parseType="Resource">
  167. <vCard:N rdf:parseType="Resource">
  168. <vCard:Family>Lloyd</vCard:Family>
  169. <vCard:Given>Catherine</vCard:Given>
  170. <vCard:Other>May</vCard:Other>
  171. </vCard:N>
  172. </cmeta:modifier>
  173. <dcterms:modified rdf:parseType="Resource">
  174. <dcterms:W3CDTF>2002-01-04</dcterms:W3CDTF>
  175. </dcterms:modified>
  176. </cmeta:modification>
  177. <cmeta:modification rdf:parseType="Resource">
  178. <rdf:value>
  179. Changed tau_y_calculation after checking mathml using the validator.
  180. </rdf:value>
  181. <cmeta:modifier rdf:parseType="Resource">
  182. <vCard:N rdf:parseType="Resource">
  183. <vCard:Family>Lloyd</vCard:Family>
  184. <vCard:Given>Catherine</vCard:Given>
  185. <vCard:Other>May</vCard:Other>
  186. </vCard:N>
  187. </cmeta:modifier>
  188. <dcterms:modified rdf:parseType="Resource">
  189. <dcterms:W3CDTF>2001-12-07</dcterms:W3CDTF>
  190. </dcterms:modified>
  191. </cmeta:modification>
  192. <cmeta:modification rdf:parseType="Resource">
  193. <rdf:value>
  194. Removed document type definition as this is declared as optional
  195. according to the W3C recommendation.
  196. </rdf:value>
  197. <cmeta:modifier rdf:parseType="Resource">
  198. <vCard:N rdf:parseType="Resource">
  199. <vCard:Family>Lloyd</vCard:Family>
  200. <vCard:Given>Catherine</vCard:Given>
  201. <vCard:Other>May</vCard:Other>
  202. </vCard:N>
  203. </cmeta:modifier>
  204. <dcterms:modified rdf:parseType="Resource">
  205. <dcterms:W3CDTF>2001-10-19</dcterms:W3CDTF>
  206. </dcterms:modified>
  207. </cmeta:modification>
  208. <!-- The Publisher metadata. -->
  209. <dc:publisher>
  210. The University of Auckland, Bioengineering Research Group
  211. </dc:publisher>
  212. </rdf:Description>
  213. <!--
  214. The following metadata refers to the model itself, as indicated by the
  215. reference to the ID "jafri_rice_winslow_model_1998", which is declared on
  216. the <model> element.
  217. -->
  218. <rdf:Description rdf:about="#jafri_rice_winslow_model_1998">
  219. <!-- A human readable name for the model. -->
  220. <dc:title>
  221. The Jafri-Rice-Winslow Model for Calcium Regulation in the Ventricular
  222. Myocyte, 1997
  223. </dc:title>
  224. <!-- A comment regarding the model. -->
  225. <cmeta:comment rdf:parseType="Resource">
  226. <rdf:value>
  227. This is the CellML description of Jafri, Rice and Winslow's
  228. mathematical model for calcium regulation in the ventricular myocyte.
  229. It is based on an accurate model of the membrane currents and adds a
  230. more sophisticated model of calcium handling. The JRW model is based
  231. on the LR-II model for ventricular action potentials, with several
  232. modifications.
  233. </rdf:value>
  234. <!-- The creator of the comment. -->
  235. <dc:creator rdf:parseType="Resource">
  236. <vCard:FN>Catherine Lloyd</vCard:FN>
  237. </dc:creator>
  238. </cmeta:comment>
  239. <!-- The model is relevant to all mammals. -->
  240. <cmeta:species>Mammalia</cmeta:species>
  241. <cmeta:bio_entity>Ventricular Myocyte</cmeta:bio_entity>
  242. <!--
  243. The CellML Metadata Specification recommends that bibliographic metadata
  244. is used to provide information about the original model reference. The
  245. "identifier" attribute on the "BibliographicReference" class provides an
  246. elegant way to identify a cited reference using a database identifier such as Pubmed. All associated data such as author, journal title,
  247. date, etc can be looked up on the database.
  248. -->
  249. <bqs:reference rdf:parseType="Resource">
  250. <bqs:Pubmed_id>9512016</bqs:Pubmed_id>
  251. <bqs:JournalArticle rdf:parseType="Resource">
  252. <dc:creator>
  253. <rdf:Seq>
  254. <rdf:li rdf:parseType="Resource">
  255. <bqs:Person rdf:parseType="Resource">
  256. <vCard:N rdf:parseType="Resource">
  257. <vCard:Family>Jafri</vCard:Family>
  258. <vCard:Given>M</vCard:Given>
  259. <vCard:Other>Saleet</vCard:Other>
  260. </vCard:N>
  261. </bqs:Person>
  262. </rdf:li>
  263. <rdf:li rdf:parseType="Resource">
  264. <bqs:Person rdf:parseType="Resource">
  265. <vCard:N rdf:parseType="Resource">
  266. <vCard:Family>Rice</vCard:Family>
  267. <vCard:Given>John</vCard:Given>
  268. <vCard:Other>Jeremy</vCard:Other>
  269. </vCard:N>
  270. </bqs:Person>
  271. </rdf:li>
  272. <rdf:li rdf:parseType="Resource">
  273. <bqs:Person rdf:parseType="Resource">
  274. <vCard:N rdf:parseType="Resource">
  275. <vCard:Family>Winslow</vCard:Family>
  276. <vCard:Given>Raimond</vCard:Given>
  277. <vCard:Other>L</vCard:Other>
  278. </vCard:N>
  279. </bqs:Person>
  280. </rdf:li>
  281. </rdf:Seq>
  282. </dc:creator>
  283. <dc:title>
  284. Cardiac Ca2+ Dynamics: The Roles of Ryanodine Receptor Adaptation
  285. and Sarcoplasmic Reticulum Load
  286. </dc:title>
  287. <dcterms:issued rdf:parseType="Resource">
  288. <dcterms:W3CDTF>1998-03</dcterms:W3CDTF>
  289. </dcterms:issued>
  290. <bqs:Journal rdf:parseType="Resource">
  291. <dc:title>Biophysical Journal</dc:title>
  292. </bqs:Journal>
  293. <bqs:volume>74</bqs:volume>
  294. <bqs:first_page>1149</bqs:first_page>
  295. <bqs:last_page>1168</bqs:last_page>
  296. </bqs:JournalArticle>
  297. </bqs:reference>
  298. </rdf:Description>
  299. </rdf:RDF>
  300. <!--
  301. Below, we define some additional units for association with variables and
  302. constants within the model. The identifiers are fairly self-explanatory.
  303. -->
  304. <units name="ms">
  305. <unit prefix="milli" units="second" />
  306. </units>
  307. <units name="per_ms">
  308. <unit prefix="milli" units="second" exponent="-1" />
  309. </units>
  310. <units name="per_mm">
  311. <unit prefix="milli" units="metre" exponent="-1" />
  312. </units>
  313. <units name="mV">
  314. <unit prefix="milli" units="volt" />
  315. </units>
  316. <units name="per_mV">
  317. <unit prefix="milli" units="volt" exponent="-1" />
  318. </units>
  319. <units name="per_mV_ms">
  320. <unit units="mV" exponent="-1" />
  321. <unit units="ms" exponent="-1" />
  322. </units>
  323. <units name="mS_per_mm2">
  324. <unit prefix="milli" units="siemens" />
  325. <unit prefix="milli" units="metre" exponent="-2" />
  326. </units>
  327. <units name="mm2">
  328. <unit prefix="milli" units="metre" exponent="2" />
  329. </units>
  330. <units name="mM">
  331. <unit prefix="nano" units="mole" />
  332. <unit prefix="milli" units="metre" exponent="-3" />
  333. </units>
  334. <units name="mM_per_ms">
  335. <unit units="mM" />
  336. <unit units="ms" exponent="-1" />
  337. </units>
  338. <units name="per_mM_per_ms">
  339. <unit units="mM" exponent="-1" />
  340. <unit units="ms" exponent="-1" />
  341. </units>
  342. <units name="per_mM3_per_ms">
  343. <unit units="mM" exponent="-3" />
  344. <unit units="ms" exponent="-1" />
  345. </units>
  346. <units name="per_mM4_per_ms">
  347. <unit units="mM" exponent="-4" />
  348. <unit units="ms" exponent="-1" />
  349. </units>
  350. <units name="uF_per_mm2">
  351. <unit prefix="micro" units="farad" />
  352. <unit prefix="milli" units="metre" exponent="-2" />
  353. </units>
  354. <units name="uA_per_mm2">
  355. <unit prefix="micro" units="ampere" />
  356. <unit prefix="milli" units="metre" exponent="-2" />
  357. </units>
  358. <units name="mm_per_ms">
  359. <unit prefix="milli" units="metre" />
  360. <unit units="ms" exponent="-1" />
  361. </units>
  362. <units name="gas_constant_units">
  363. <unit prefix="pico" units="joule" />
  364. <unit prefix="nano" units="mole" exponent="-1" />
  365. <unit units="kelvin" exponent="-1" />
  366. </units>
  367. <units name="faradays_constant_units">
  368. <unit prefix="nano" units="coulomb" />
  369. <unit prefix="nano" units="mole" exponent="-1" />
  370. </units>
  371. <!--
  372. The "environment" component is used to declare variables that are used by
  373. all or most of the other components, in this case just "time".
  374. -->
  375. <component name="environment">
  376. <variable name="time" public_interface="out" units="ms" />
  377. </component>
  378. <!--
  379. The "membrane" component is really the `root' node of our model.
  380. It defines the action potential variable "V" among other things.
  381. -->
  382. <component name="membrane" cmeta:id="membrane">
  383. <rdf:RDF
  384. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  385. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  386. <rdf:Description rdf:about="#membrane">
  387. <cmeta:comment rdf:parseType="Resource">
  388. <rdf:value>
  389. The main component of the model which defines the action potential.
  390. </rdf:value>
  391. </cmeta:comment>
  392. </rdf:Description>
  393. </rdf:RDF>
  394. <!-- These variables are defined here and used in other components. -->
  395. <variable
  396. name="V" public_interface="out"
  397. initial_value="-84.1638" units="mV" />
  398. <variable
  399. name="R" public_interface="out"
  400. initial_value="8.3145e3" units="gas_constant_units" />
  401. <variable
  402. name="T" public_interface="out"
  403. initial_value="310.0" units="kelvin" />
  404. <variable
  405. name="F" public_interface="out"
  406. initial_value="9.6845e4" units="faradays_constant_units" />
  407. <!-- These variables are imported from other components. -->
  408. <variable name="time" public_interface="in" units="ms" />
  409. <variable name="i_Na" public_interface="in" units="uA_per_mm2" />
  410. <variable name="i_Ca_L_Ca" public_interface="in" units="uA_per_mm2" />
  411. <variable name="i_Ca_L_K" public_interface="in" units="uA_per_mm2" />
  412. <variable name="i_K" public_interface="in" units="uA_per_mm2" />
  413. <variable name="i_K1" public_interface="in" units="uA_per_mm2" />
  414. <variable name="i_NaCa" public_interface="in" units="uA_per_mm2" />
  415. <variable name="i_Kp" public_interface="in" units="uA_per_mm2" />
  416. <variable name="i_p_Ca" public_interface="in" units="uA_per_mm2" />
  417. <variable name="i_Na_b" public_interface="in" units="uA_per_mm2" />
  418. <variable name="i_Ca_b" public_interface="in" units="uA_per_mm2" />
  419. <variable name="i_NaK" public_interface="in" units="uA_per_mm2" />
  420. <variable name="i_ns_Ca" public_interface="in" units="uA_per_mm2" />
  421. <!-- These variables is defined here and only used internally. -->
  422. <variable name="Cm" initial_value="0.01" units="uF_per_mm2" />
  423. <variable name="I_stim" units="uA_per_mm2" />
  424. <rdf:RDF
  425. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  426. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  427. <rdf:Description rdf:about="#membrane_voltage_diff_eq">
  428. <cmeta:comment rdf:parseType="Resource">
  429. <rdf:value>
  430. The kinetics of the transmembrane potential, defined as the sum of
  431. all the sarcolemmal currents and an applied stimulus current.
  432. </rdf:value>
  433. </cmeta:comment>
  434. </rdf:Description>
  435. </rdf:RDF>
  436. <math xmlns="http://www.w3.org/1998/Math/MathML"
  437. cmeta:id="membrane_voltage_diff_eq">
  438. <apply id="membrane_voltage_diff_eq"><eq />
  439. <apply><diff />
  440. <bvar><ci> time </ci></bvar>
  441. <ci> V </ci>
  442. </apply>
  443. <apply><divide />
  444. <apply><minus />
  445. <ci>I_stim</ci>
  446. <apply><plus />
  447. <ci> i_Na </ci>
  448. <ci> i_Ca_L_Ca </ci>
  449. <ci> i_Ca_L_K </ci>
  450. <ci> i_K </ci>
  451. <ci> i_NaCa </ci>
  452. <ci> i_K1 </ci>
  453. <ci> i_Kp </ci>
  454. <ci> i_p_Ca </ci>
  455. <ci> i_Na_b </ci>
  456. <ci> i_Ca_b </ci>
  457. <ci> i_NaK </ci>
  458. <ci> i_ns_Ca </ci>
  459. </apply>
  460. </apply>
  461. <ci> Cm </ci>
  462. </apply>
  463. </apply>
  464. </math>
  465. <rdf:RDF
  466. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  467. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  468. <rdf:Description rdf:about="#IStim_for_cmiss_eq">
  469. <cmeta:comment rdf:parseType="Resource">
  470. <rdf:value>
  471. This is a dummy equation that we simply use to make grabbing the
  472. value in CMISS much easier.
  473. </rdf:value>
  474. </cmeta:comment>
  475. </rdf:Description>
  476. </rdf:RDF>
  477. <variable name="IStimC" units="uA_per_mm2" public_interface="out"/>
  478. <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="IStim_for_cmiss_eq">
  479. <apply id="IStim_for_cmiss"><eq/>
  480. <ci>IStimC</ci>
  481. <ci>I_stim</ci>
  482. </apply>
  483. </math>
  484. </component>
  485. <component name="fast_sodium_current" cmeta:id="fast_sodium_current">
  486. <rdf:RDF
  487. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  488. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  489. <rdf:Description rdf:about="#fast_sodium_current">
  490. <cmeta:comment rdf:parseType="Resource">
  491. <rdf:value>
  492. The fast sodium current component contains the differential
  493. equations governing the influx of sodium ions through the cell
  494. surface membrane into the cell.
  495. </rdf:value>
  496. </cmeta:comment>
  497. </rdf:Description>
  498. </rdf:RDF>
  499. <!-- These variables are defined here and used in other components. -->
  500. <variable name="i_Na" public_interface="out" units="uA_per_mm2" />
  501. <variable name="E_Na" public_interface="out" units="mV" />
  502. <!-- This variable is defined here and only used internally. -->
  503. <variable name="g_Na" initial_value="0.128" units="mS_per_mm2" />
  504. <!--
  505. Time is imported from the environment, and membrane potential is
  506. imported from the membrane component. These variables are used in the
  507. "sodium_current" parent component, which also acts as an interface,
  508. passing the variables to its encapsulated gate components.
  509. -->
  510. <variable
  511. name="time" public_interface="in"
  512. private_interface="out" units="ms" />
  513. <variable
  514. name="V" public_interface="in"
  515. private_interface="out" units="mV" />
  516. <variable name="R" public_interface="in" units="gas_constant_units" />
  517. <variable name="T" public_interface="in" units="kelvin" />
  518. <variable name="F" public_interface="in" units="faradays_constant_units" />
  519. <variable name="Nai" public_interface="in" units="mM" />
  520. <variable name="Nao" public_interface="in" units="mM" />
  521. <!-- These variables are imported from encapsulated components. -->
  522. <variable name="m" private_interface="in" units="dimensionless" />
  523. <variable name="h" private_interface="in" units="dimensionless" />
  524. <variable name="j" private_interface="in" units="dimensionless" />
  525. <rdf:RDF
  526. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  527. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  528. <rdf:Description rdf:about="#i_Na_calculation_eq">
  529. <cmeta:comment rdf:parseType="Resource">
  530. <rdf:value>
  531. Calculation of the fast sodium current using the three
  532. Hodkin-Huxley type voltage-dependent gating variables m, h, and j.
  533. </rdf:value>
  534. </cmeta:comment>
  535. </rdf:Description>
  536. </rdf:RDF>
  537. <math xmlns="http://www.w3.org/1998/Math/MathML"
  538. cmeta:id="i_Na_calculation_eq">
  539. <apply id="i_Na_calculation"><eq />
  540. <ci> i_Na </ci>
  541. <apply><times />
  542. <ci> g_Na </ci>
  543. <apply><power />
  544. <ci> m </ci>
  545. <cn cellml:units="dimensionless"> 3.0 </cn>
  546. </apply>
  547. <ci> h </ci>
  548. <ci> j </ci>
  549. <apply><minus />
  550. <ci> V </ci>
  551. <ci> E_Na </ci>
  552. </apply>
  553. </apply>
  554. </apply>
  555. </math>
  556. <rdf:RDF
  557. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  558. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  559. <rdf:Description rdf:about="#E_Na_calculation_eq">
  560. <cmeta:comment rdf:parseType="Resource">
  561. <rdf:value>
  562. The sodium reversal potential.
  563. </rdf:value>
  564. </cmeta:comment>
  565. </rdf:Description>
  566. </rdf:RDF>
  567. <math xmlns="http://www.w3.org/1998/Math/MathML"
  568. cmeta:id="E_Na_calculation_eq">
  569. <apply id="E_Na_calculation"><eq />
  570. <ci> E_Na </ci>
  571. <apply><times />
  572. <apply><divide />
  573. <apply><times />
  574. <ci> R </ci>
  575. <ci> T </ci>
  576. </apply>
  577. <ci> F </ci>
  578. </apply>
  579. <apply><ln />
  580. <apply><divide />
  581. <ci> Nao </ci>
  582. <ci> Nai </ci>
  583. </apply>
  584. </apply>
  585. </apply>
  586. </apply>
  587. </math>
  588. </component>
  589. <!--
  590. The "fast_sodium_current_m_gate" is the m gate encapsulated inside the fast
  591. sodium current.
  592. -->
  593. <component name="fast_sodium_current_m_gate"
  594. cmeta:id="fast_sodium_current_m_gate">
  595. <rdf:RDF
  596. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  597. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  598. <rdf:Description rdf:about="#fast_sodium_current_m_gate">
  599. <cmeta:comment rdf:parseType="Resource">
  600. <rdf:value>
  601. The voltage-dependent activation gate for the fast sodium current -
  602. the m gate.
  603. </rdf:value>
  604. </cmeta:comment>
  605. </rdf:Description>
  606. </rdf:RDF>
  607. <!-- This variable is defined here and used in other components. -->
  608. <variable
  609. name="m" public_interface="out"
  610. initial_value="0.0328302" units="dimensionless" />
  611. <!-- These variables are defined here and only used internally. -->
  612. <variable name="alpha_m" units="per_ms" />
  613. <variable name="beta_m" units="per_ms" />
  614. <!--
  615. These variables are imported from the environment and the membrane via
  616. the "fast_sodium_current" component.
  617. -->
  618. <variable name="V" public_interface="in" units="mV" />
  619. <variable name="time" public_interface="in" units="ms" />
  620. <rdf:RDF
  621. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  622. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  623. <rdf:Description rdf:about="#alpha_m_calculation_eq">
  624. <cmeta:comment rdf:parseType="Resource">
  625. <rdf:value>
  626. The opening rate of the m gate.
  627. </rdf:value>
  628. </cmeta:comment>
  629. </rdf:Description>
  630. </rdf:RDF>
  631. <math xmlns="http://www.w3.org/1998/Math/MathML"
  632. cmeta:id="alpha_m_calculation_eq">
  633. <apply id="alpha_m_calculation"><eq />
  634. <ci> alpha_m </ci>
  635. <apply><divide />
  636. <apply><times />
  637. <cn cellml:units="per_mV_ms"> 0.32 </cn>
  638. <apply><plus />
  639. <ci> V </ci>
  640. <cn cellml:units="mV"> 47.13 </cn>
  641. </apply>
  642. </apply>
  643. <apply><minus />
  644. <cn cellml:units="dimensionless"> 1.0 </cn>
  645. <apply><exp />
  646. <apply><times />
  647. <cn cellml:units="per_mV"> -0.1 </cn>
  648. <apply><plus />
  649. <ci> V </ci>
  650. <cn cellml:units="mV"> 47.13 </cn>
  651. </apply>
  652. </apply>
  653. </apply>
  654. </apply>
  655. </apply>
  656. </apply>
  657. </math>
  658. <rdf:RDF
  659. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  660. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  661. <rdf:Description rdf:about="#beta_m_calculation_eq">
  662. <cmeta:comment rdf:parseType="Resource">
  663. <rdf:value>
  664. The closing rate of the m gate.
  665. </rdf:value>
  666. </cmeta:comment>
  667. </rdf:Description>
  668. </rdf:RDF>
  669. <math xmlns="http://www.w3.org/1998/Math/MathML"
  670. cmeta:id="beta_m_calculation_eq">
  671. <apply id="beta_m_calculation"><eq />
  672. <ci> beta_m </ci>
  673. <apply><times />
  674. <cn cellml:units="per_ms"> 0.08 </cn>
  675. <apply><exp />
  676. <apply><divide />
  677. <apply><minus />
  678. <ci> V </ci>
  679. </apply>
  680. <cn cellml:units="mV"> 11.0 </cn>
  681. </apply>
  682. </apply>
  683. </apply>
  684. </apply>
  685. </math>
  686. <rdf:RDF
  687. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  688. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  689. <rdf:Description rdf:about="#dm_dt_eq">
  690. <cmeta:comment rdf:parseType="Resource">
  691. <rdf:value>
  692. The kinetics of the m gate.
  693. </rdf:value>
  694. </cmeta:comment>
  695. </rdf:Description>
  696. </rdf:RDF>
  697. <math xmlns="http://www.w3.org/1998/Math/MathML"
  698. cmeta:id="dm_dt_eq">
  699. <apply id="dm_dt"><eq />
  700. <apply><diff />
  701. <bvar><ci> time </ci></bvar>
  702. <ci> m </ci>
  703. </apply>
  704. <apply><minus />
  705. <apply><times />
  706. <ci> alpha_m </ci>
  707. <apply><minus />
  708. <cn cellml:units="per_mV"> 1.0 </cn>
  709. <ci> m </ci>
  710. </apply>
  711. </apply>
  712. <apply><times />
  713. <ci> beta_m </ci>
  714. <ci> m </ci>
  715. </apply>
  716. </apply>
  717. </apply>
  718. </math>
  719. </component>
  720. <!--
  721. The "fast_sodium_current_h_gate" component is the h gate encapsulated in
  722. the fast sodium current.
  723. -->
  724. <component name="fast_sodium_current_h_gate"
  725. cmeta:id="fast_sodium_current_h_gate">
  726. <rdf:RDF
  727. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  728. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  729. <rdf:Description rdf:about="#fast_sodium_current_h_gate">
  730. <cmeta:comment rdf:parseType="Resource">
  731. <rdf:value>
  732. The voltage-dependent inactivation gate for the fast sodium current
  733. - the h gate.
  734. </rdf:value>
  735. </cmeta:comment>
  736. </rdf:Description>
  737. </rdf:RDF>
  738. <!-- This variable is defined here and used in other components. -->
  739. <variable
  740. name="h" public_interface="out"
  741. initial_value="0.988354" units="dimensionless" />
  742. <!-- These variables are defined here and only used internally. -->
  743. <variable name="alpha_h" units="per_ms" />
  744. <variable name="beta_h" units="per_ms" />
  745. <!--
  746. These variables are imported from the environment and the membrane via
  747. the "fast_sodium_current" component.
  748. -->
  749. <variable name="V" public_interface="in" units="mV" />
  750. <variable name="time" public_interface="in" units="ms" />
  751. <rdf:RDF
  752. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  753. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  754. <rdf:Description rdf:about="#alpha_h_calculation_eq">
  755. <cmeta:comment rdf:parseType="Resource">
  756. <rdf:value>
  757. The opening rate of the h gate.
  758. </rdf:value>
  759. </cmeta:comment>
  760. </rdf:Description>
  761. </rdf:RDF>
  762. <math xmlns="http://www.w3.org/1998/Math/MathML"
  763. cmeta:id="alpha_h_calculation_eq">
  764. <apply id="alpha_h_calculation"><eq />
  765. <ci> alpha_h </ci>
  766. <piecewise>
  767. <piece>
  768. <apply><times />
  769. <cn cellml:units="per_ms"> 0.135 </cn>
  770. <apply><exp />
  771. <apply><divide />
  772. <apply><plus />
  773. <cn cellml:units="mV"> 80.0 </cn>
  774. <ci> V </ci>
  775. </apply>
  776. <cn cellml:units="mV"> -6.8 </cn>
  777. </apply>
  778. </apply>
  779. </apply>
  780. <apply><lt />
  781. <ci> V </ci>
  782. <cn cellml:units="mV"> -40.0 </cn>
  783. </apply>
  784. </piece>
  785. <otherwise>
  786. <cn cellml:units="per_ms"> 0.0 </cn>
  787. </otherwise>
  788. </piecewise>
  789. </apply>
  790. </math>
  791. <rdf:RDF
  792. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  793. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  794. <rdf:Description rdf:about="#beta_h_calculation_eq">
  795. <cmeta:comment rdf:parseType="Resource">
  796. <rdf:value>
  797. The closing rate of the h gate.
  798. </rdf:value>
  799. </cmeta:comment>
  800. </rdf:Description>
  801. </rdf:RDF>
  802. <math xmlns="http://www.w3.org/1998/Math/MathML"
  803. cmeta:id="beta_h_calculation_eq">
  804. <apply id="beta_h_calculation"><eq />
  805. <ci> beta_h </ci>
  806. <piecewise>
  807. <piece>
  808. <apply><plus />
  809. <apply><times />
  810. <cn cellml:units="per_ms"> 3.56 </cn>
  811. <apply><exp />
  812. <apply><times />
  813. <cn cellml:units="mV"> 0.079 </cn>
  814. <ci> V </ci>
  815. </apply>
  816. </apply>
  817. </apply>
  818. <apply><times />
  819. <cn cellml:units="per_ms"> 310000.0 </cn>
  820. <apply><exp />
  821. <apply><times />
  822. <cn cellml:units="per_mV"> 0.35 </cn>
  823. <ci> V </ci>
  824. </apply>
  825. </apply>
  826. </apply>
  827. </apply>
  828. <apply><lt />
  829. <ci> V </ci>
  830. <cn cellml:units="mV"> -40.0 </cn>
  831. </apply>
  832. </piece>
  833. <otherwise>
  834. <apply><divide />
  835. <cn cellml:units="dimensionless"> 1.0 </cn>
  836. <apply><times />
  837. <cn cellml:units="ms"> 0.13 </cn>
  838. <apply><plus />
  839. <cn cellml:units="dimensionless"> 1.0 </cn>
  840. <apply><exp />
  841. <apply><divide />
  842. <apply><plus />
  843. <ci> V </ci>
  844. <cn cellml:units="mV"> 10.66 </cn>
  845. </apply>
  846. <cn cellml:units="mV"> -11.1 </cn>
  847. </apply>
  848. </apply>
  849. </apply>
  850. </apply>
  851. </apply>
  852. </otherwise>
  853. </piecewise>
  854. </apply>
  855. </math>
  856. <rdf:RDF
  857. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  858. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  859. <rdf:Description rdf:about="#dh_dt_eq">
  860. <cmeta:comment rdf:parseType="Resource">
  861. <rdf:value>
  862. The kinetics of the h gate.
  863. </rdf:value>
  864. </cmeta:comment>
  865. </rdf:Description>
  866. </rdf:RDF>
  867. <math xmlns="http://www.w3.org/1998/Math/MathML"
  868. cmeta:id="dh_dt_eq">
  869. <apply id="dh_dt"><eq />
  870. <apply><diff />
  871. <bvar><ci> time </ci></bvar>
  872. <ci> h </ci>
  873. </apply>
  874. <apply><minus />
  875. <apply><times />
  876. <ci> alpha_h </ci>
  877. <apply><minus />
  878. <cn cellml:units="dimensionless"> 1.0 </cn>
  879. <ci> h </ci>
  880. </apply>
  881. </apply>
  882. <apply><times />
  883. <ci> beta_h </ci>
  884. <ci> h </ci>
  885. </apply>
  886. </apply>
  887. </apply>
  888. </math>
  889. </component>
  890. <!--
  891. The "fast_sodium_current_j_gate" component is the j gate encapsulated in
  892. the fast sodium current.
  893. -->
  894. <component name="fast_sodium_current_j_gate"
  895. cmeta:id="fast_sodium_current_j_gate">
  896. <rdf:RDF
  897. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  898. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  899. <rdf:Description rdf:about="#fast_sodium_current_j_gate">
  900. <cmeta:comment rdf:parseType="Resource">
  901. <rdf:value>
  902. The voltage-dependent slow inactivation gate for the fast sodium
  903. current - the j gate.
  904. </rdf:value>
  905. </cmeta:comment>
  906. </rdf:Description>
  907. </rdf:RDF>
  908. <!-- This variable is defined here and used in other components. -->
  909. <variable
  910. name="j" public_interface="out"
  911. initial_value="0.992540" units="dimensionless" />
  912. <!-- These variables are defined here and only used internally. -->
  913. <variable name="alpha_j" units="per_ms" />
  914. <variable name="beta_j" units="per_ms" />
  915. <!--
  916. These variables are imported from the environment and the membrane via
  917. the "fast_sodium_current" component.
  918. -->
  919. <variable name="V" public_interface="in" units="mV" />
  920. <variable name="time" public_interface="in" units="ms" />
  921. <rdf:RDF
  922. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  923. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  924. <rdf:Description rdf:about="#alpha_j_calculation_eq">
  925. <cmeta:comment rdf:parseType="Resource">
  926. <rdf:value>
  927. The opening rate of the j gate.
  928. </rdf:value>
  929. </cmeta:comment>
  930. </rdf:Description>
  931. </rdf:RDF>
  932. <math xmlns="http://www.w3.org/1998/Math/MathML"
  933. cmeta:id="alpha_j_calculation_eq">
  934. <apply id="alpha_j_calculation"><eq />
  935. <ci> alpha_j </ci>
  936. <piecewise>
  937. <piece>
  938. <apply><times />
  939. <apply><minus />
  940. <apply><times />
  941. <cn cellml:units="per_mV_ms"> -127140.0 </cn>
  942. <apply><exp />
  943. <apply><times />
  944. <cn cellml:units="per_mV"> 0.2444 </cn>
  945. <ci> V </ci>
  946. </apply>
  947. </apply>
  948. </apply>
  949. <apply><times />
  950. <cn cellml:units="per_mV_ms"> 0.00003474 </cn>
  951. <apply><exp />
  952. <apply><times />
  953. <cn cellml:units="per_mV"> -0.04391 </cn>
  954. <ci> V </ci>
  955. </apply>
  956. </apply>
  957. </apply>
  958. </apply>
  959. <apply><divide />
  960. <apply><plus />
  961. <ci> V </ci>
  962. <cn cellml:units="mV"> 37.78 </cn>
  963. </apply>
  964. <apply><plus />
  965. <cn cellml:units="dimensionless"> 1.0 </cn>
  966. <apply><exp />
  967. <apply><times />
  968. <cn cellml:units="per_mV"> 0.311 </cn>
  969. <apply><plus />
  970. <ci> V </ci>
  971. <cn cellml:units="mV"> 79.23 </cn>
  972. </apply>
  973. </apply>
  974. </apply>
  975. </apply>
  976. </apply>
  977. </apply>
  978. <apply><lt />
  979. <ci> V </ci>
  980. <cn cellml:units="mV"> -40.0 </cn>
  981. </apply>
  982. </piece>
  983. <otherwise>
  984. <cn cellml:units="per_ms"> 0.0 </cn>
  985. </otherwise>
  986. </piecewise>
  987. </apply>
  988. </math>
  989. <rdf:RDF
  990. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  991. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  992. <rdf:Description rdf:about="#beta_j_calculation_eq">
  993. <cmeta:comment rdf:parseType="Resource">
  994. <rdf:value>
  995. The closing rate of the j gate.
  996. </rdf:value>
  997. </cmeta:comment>
  998. </rdf:Description>
  999. </rdf:RDF>
  1000. <math xmlns="http://www.w3.org/1998/Math/MathML"
  1001. cmeta:id="beta_j_calculation_eq">
  1002. <apply id="beta_j_calculation"><eq />
  1003. <ci> beta_j </ci>
  1004. <piecewise>
  1005. <piece>
  1006. <apply><times />
  1007. <cn cellml:units="per_ms"> 0.1212 </cn>
  1008. <apply><divide />
  1009. <apply><exp />
  1010. <apply><times />
  1011. <cn cellml:units="per_mV"> -0.01052 </cn>
  1012. <ci> V </ci>
  1013. </apply>
  1014. </apply>
  1015. <apply><plus />
  1016. <cn cellml:units="dimensionless"> 1.0 </cn>
  1017. <apply><exp />
  1018. <apply><times />
  1019. <cn cellml:units="per_mV"> -0.1378 </cn>
  1020. <apply><plus />
  1021. <ci> V </ci>
  1022. <cn cellml:units="mV"> 40.14 </cn>
  1023. </apply>
  1024. </apply>
  1025. </apply>
  1026. </apply>
  1027. </apply>
  1028. </apply>
  1029. <apply><lt />
  1030. <ci> V </ci>
  1031. <cn cellml:units="mV"> -40.0 </cn>
  1032. </apply>
  1033. </piece>
  1034. <otherwise>
  1035. <apply><times />
  1036. <cn cellml:units="per_ms"> 0.3 </cn>
  1037. <apply><divide />
  1038. <apply><exp />
  1039. <apply><times />
  1040. <cn cellml:units="per_mV"> -0.0000002535 </cn>
  1041. <ci> V </ci>
  1042. </apply>
  1043. </apply>
  1044. <apply><plus />
  1045. <cn cellml:units="dimensionless"> 1.0 </cn>
  1046. <apply><exp />
  1047. <apply><times />
  1048. <cn cellml:units="per_mV"> -0.1 </cn>
  1049. <apply><plus />
  1050. <ci> V </ci>
  1051. <cn cellml:units="mV"> 32.0 </cn>
  1052. </apply>
  1053. </apply>
  1054. </apply>
  1055. </apply>
  1056. </apply>
  1057. </apply>
  1058. </otherwise>
  1059. </piecewise>
  1060. </apply>
  1061. </math>
  1062. <rdf:RDF
  1063. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  1064. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  1065. <rdf:Description rdf:about="#dj_dt_eq">
  1066. <cmeta:comment rdf:parseType="Resource">
  1067. <rdf:value>
  1068. The kinetics of the j gate.
  1069. </rdf:value>
  1070. </cmeta:comment>
  1071. </rdf:Description>
  1072. </rdf:RDF>
  1073. <math xmlns="http://www.w3.org/1998/Math/MathML"
  1074. cmeta:id="dj_dt_eq">
  1075. <apply id="dj_dt"><eq />
  1076. <apply><diff />
  1077. <bvar><ci> time </ci></bvar>
  1078. <ci> j </ci>
  1079. </apply>
  1080. <apply><minus />
  1081. <apply><times />
  1082. <ci> alpha_j </ci>
  1083. <apply><minus />
  1084. <cn cellml:units="dimensionless"> 1.0 </cn>
  1085. <ci> j </ci>
  1086. </apply>
  1087. </apply>
  1088. <apply><times />
  1089. <ci> beta_j </ci>
  1090. <ci> j </ci>
  1091. </apply>
  1092. </apply>
  1093. </apply>
  1094. </math>
  1095. </component>
  1096. <component name="L_type_Ca_channel" cmeta:id="L_type_Ca_channel">
  1097. <rdf:RDF
  1098. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  1099. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  1100. <rdf:Description rdf:about="#L_type_Ca_channel">
  1101. <cmeta:comment rdf:parseType="Resource">
  1102. <rdf:value>
  1103. The JWR model creates a new mathematical model to describe the
  1104. L-type calcium channel that is based on the experimentally observed
  1105. mode-switching behaviour of the channel. Inactivation occurs as
  1106. calcium ion binding induces the channel to switch (from mode normal)
  1107. to a mode in which transitions to open states are extremely slow
  1108. (mode Ca). The channel has one voltage inactivation gate, y. As well
  1109. as Ca, the channel is assumed permeable to K ions also.
  1110. </rdf:value>
  1111. </cmeta:comment>
  1112. </rdf:Description>
  1113. </rdf:RDF>
  1114. <!-- These variables are defined here and used in other components. -->
  1115. <variable name="i_Ca_L_Ca" public_interface="out" units="uA_per_mm2" />
  1116. <variable name="i_Ca_L_K" public_interface="out" units="uA_per_mm2" />
  1117. <!-- These variables are defined here and only used internally. -->
  1118. <variable name="P_Ca" initial_value="33.75e-6" units="mm_per_ms" />
  1119. <variable name="P_K" initial_value="1.0e-9" units="mm_per_ms" />
  1120. <variable name="p_k" units="mm_per_ms" />
  1121. <variable name="i_Ca_L_Ca_half" initial_value="-4.58e-3" units="uA_per_mm2" />
  1122. <variable name="i_Ca_L_Ca_max" units="uA_per_mm2" />
  1123. <variable name="O" initial_value="9.84546e-21" units="dimensionless" />
  1124. <variable name="O_Ca" initial_value="0.0" units="dimensionless" />
  1125. <variable name="alpha" units="per_ms" />
  1126. <variable name="beta" units="per_ms" />
  1127. <variable name="gamma" units="per_ms" />
  1128. <variable name="alpha_a" units="per_ms" />
  1129. <variable name="beta_b" units="per_ms" />
  1130. <variable name="a" initial_value="2.0" units="dimensionless" />
  1131. <variable name="b" initial_value="2.0" units="dimensionless" />
  1132. <variable name="g" initial_value="2.0" units="per_ms" />
  1133. <variable name="f" initial_value="0.3" units="per_ms" />
  1134. <variable name="g_" initial_value="0.0" units="per_ms" />
  1135. <variable name="f_" initial_value="0.0" units="per_ms" />
  1136. <variable name="omega" initial_value="0.01" units="per_ms" />
  1137. <variable name="C0" initial_value="0.997208" units="dimensionless" />
  1138. <variable name="C1" initial_value="6.38897e-5" units="dimensionless" />
  1139. <variable name="C2" initial_value="1.535e-9" units="dimensionless" />
  1140. <variable name="C3" initial_value="1.63909e-14" units="dimensionless" />
  1141. <variable name="C4" initial_value="6.56337e-20" units="dimensionless" />
  1142. <variable name="C_Ca0" initial_value="2.72826e-3" units="dimensionless" />
  1143. <variable name="C_Ca1" initial_value="6.99215e-7" units="dimensionless" />
  1144. <variable name="C_Ca2" initial_value="6.71989e-11" units="dimensionless" />
  1145. <variable name="C_Ca3" initial_value="2.87031e-15" units="dimensionless" />
  1146. <variable name="C_Ca4" initial_value="4.59752e-20" units="dimensionless" />
  1147. <!-- These variables are imported from other components. -->
  1148. <variable name="time" public_interface="in" private_interface="out" units="ms" />
  1149. <variable name="V" public_interface="in" private_interface="out" units="mV" />
  1150. <variable name="Ca_SS" public_interface="in" units="mM" />
  1151. <variable name="Cao" public_interface="in" units="mM" />
  1152. <variable name="Ko" public_interface="in" units="mM" />
  1153. <variable name="Ki" public_interface="in" units="mM" />
  1154. <variable name="R" public_interface="in" units="gas_constant_units" />
  1155. <variable name="T" public_interface="in" units="kelvin" />
  1156. <variable name="F" public_interface="in" units="faradays_constant_units" />
  1157. <!-- These variables are imported from encapsulated components. -->
  1158. <variable name="y" private_interface="in" units="dimensionless" />
  1159. <rdf:RDF
  1160. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  1161. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  1162. <rdf:Description rdf:about="#i_Ca_L_Ca_calculation_eq">
  1163. <cmeta:comment rdf:parseType="Resource">
  1164. <rdf:value>
  1165. Calculation of the calcium current component of the total channel
  1166. current, given as the maximal current multiplied by the
  1167. voltage-dependent inactivation gate and the open probability of the
  1168. channel based on the mode-switching model.
  1169. </rdf:value>
  1170. </cmeta:comment>
  1171. </rdf:Description>
  1172. </rdf:RDF>
  1173. <math xmlns="http://www.w3.org/1998/Math/MathML"
  1174. cmeta:id="i_Ca_L_Ca_calculation_eq">
  1175. <apply id="i_Ca_L_Ca_calculation"><eq />
  1176. <ci> i_Ca_L_Ca </ci>
  1177. <apply><times />
  1178. <ci> i_Ca_L_Ca_max </ci>
  1179. <ci> y </ci>
  1180. <apply><plus />
  1181. <ci> O </ci>
  1182. <ci> O_Ca </ci>
  1183. </apply>
  1184. </apply>
  1185. </apply>
  1186. </math>
  1187. <rdf:RDF
  1188. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  1189. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  1190. <rdf:Description rdf:about="#i_Ca_L_K_calculation_eq">
  1191. <cmeta:comment rdf:parseType="Resource">
  1192. <rdf:value>
  1193. Calculation of the potassium current component of the total channel
  1194. current.
  1195. </rdf:value>
  1196. </cmeta:comment>
  1197. </rdf:Description>
  1198. </rdf:RDF>
  1199. <math xmlns="http://www.w3.org/1998/Math/MathML"
  1200. cmeta:id="i_Ca_L_K_calculation_eq">
  1201. <apply id="i_Ca_L_K_calculation"><eq />
  1202. <ci> i_Ca_L_K </ci>
  1203. <apply><times />
  1204. <ci> p_k </ci>
  1205. <ci> y </ci>
  1206. <apply><plus />
  1207. <ci> O </ci>
  1208. <ci> O_Ca </ci>
  1209. </apply>
  1210. <apply><divide />
  1211. <apply><times />
  1212. <ci> V </ci>
  1213. <apply><power/>
  1214. <ci> F </ci>
  1215. <cn cellml:units="dimensionless"> 2.0 </cn>
  1216. </apply>
  1217. </apply>
  1218. <apply><times />
  1219. <ci> R </ci>
  1220. <ci> T </ci>
  1221. </apply>
  1222. </apply>
  1223. <apply><divide />
  1224. <apply><minus />
  1225. <apply><times />
  1226. <ci> Ki </ci>
  1227. <apply><exp />
  1228. <apply><divide />
  1229. <apply><times />
  1230. <ci> V </ci>
  1231. <ci> F </ci>
  1232. </apply>
  1233. <apply><times />
  1234. <ci> R </ci>
  1235. <ci> T </ci>
  1236. </apply>
  1237. </apply>
  1238. </apply>
  1239. </apply>
  1240. <ci> Ko </ci>
  1241. </apply>
  1242. <apply><minus />
  1243. <apply><exp />
  1244. <apply><divide />
  1245. <apply><times />
  1246. <ci> V </ci>
  1247. <ci> F </ci>
  1248. </apply>
  1249. <apply><times />
  1250. <ci> R </ci>
  1251. <ci> T </ci>
  1252. </apply>
  1253. </apply>
  1254. </apply>
  1255. <cn cellml:units="dimensionless"> 1.0 </cn>
  1256. </apply>
  1257. </apply>
  1258. </apply>
  1259. </apply>
  1260. </math>
  1261. <rdf:RDF
  1262. xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  1263. xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
  1264. <rdf:Description rdf:about="#p_K_calculation_eq">
  1265. <cmeta:comment rdf:parseType="Resource">
  1266. <rdf:value>
  1267. The potassium permeability of the channel, which depends on the
  1268. calcium current component.
  1269. </rdf:value>
  1270. </cmeta:comment>
  1271. </rdf:Description>
  1272. </rdf:RDF>
  1273. <math xmlns="http://www.w3.org/1998/Math/MathML"
  1274. cmeta:id="p_K_calculation_eq">
  1275. <apply id="p_k_calculation"><eq />
  1276. <ci> p_k </ci>
  1277. <apply><divide />
  1278. <ci> P_K </ci>
  1279. <apply><plus />
  1280. <cn cellml:units="dimensionless"> 1.0 </cn>
  1281. <apply><divide />
  1282. <ci> i_Ca_L_Ca…