Koo2013 - Integrated shear stress induced NO production model

  public model
Model Identifier
Short description
Related Publication
  • In silico modeling of shear-stress-induced nitric oxide production in endothelial cells through systems biology.
  • Koo A, Nordsletten D, Umeton R, Yankama B, Ayyadurai S, García-Cardeña G, Dewey CF Jr
  • Biophysical journal , 5/ 2013 , Volume 104 , pages: 2295-2306 , PubMed ID: 23708369
  • Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA.
  • Nitric oxide (NO) produced by vascular endothelial cells is a potent vasodilator and an antiinflammatory mediator. Regulating production of endothelial-derived NO is a complex undertaking, involving multiple signaling and genetic pathways that are activated by diverse humoral and biomechanical stimuli. To gain a thorough understanding of the rich diversity of responses observed experimentally, it is necessary to account for an ensemble of these pathways acting simultaneously. In this article, we have assembled four quantitative molecular pathways previously proposed for shear-stress-induced NO production. In these pathways, endothelial NO synthase is activated 1), via calcium release, 2), via phosphorylation reactions, and 3), via enhanced protein expression. To these activation pathways, we have added a fourth, a pathway describing actual NO production from endothelial NO synthase and its various protein partners. These pathways were combined and simulated using CytoSolve, a computational environment for combining independent pathway calculations. The integrated model is able to describe the experimentally observed change in NO production with time after the application of fluid shear stress. This model can also be used to predict the specific effects on the system after interventional pharmacological or genetic changes. Importantly, this model reflects the up-to-date understanding of the NO system, providing a platform upon which information can be aggregated in an additive way.
Vijayalakshmi Chelliah, administrator

Metadata information

BioModels Database MODEL1308190000
BioModels Database BIOMD0000000468
BioModels Database BIOMD0000000466
BioModels Database BIOMD0000000464
BioModels Database BIOMD0000000465
BioModels Database BIOMD0000000467
PubMed 23708369
Taxonomy Homo sapiens

Curation status

Name Description Size Actions

Model files

BIOMD0000000468_url.xml SBML L2V4 representation of Koo2013 - Integrated shear stress induced NO production model 418.71 KB Preview | Download

Additional files

BIOMD0000000468-biopax2.owl Auto-generated BioPAX (Level 2) 106.99 KB Preview | Download
BIOMD0000000468.svg Auto-generated Reaction graph (SVG) 252.53 KB Preview | Download
BIOMD0000000468.sci Auto-generated Scilab file 205.00 Bytes Preview | Download
Integratedmodel_shearstress.cps Copasi file for the plot that corresponds to NO production under shear stress condition (red). 755.50 KB Preview | Download
BIOMD0000000468_urn.xml Auto-generated SBML file with URNs 416.65 KB Preview | Download
BIOMD0000000468.vcml Auto-generated VCML file 882.00 Bytes Preview | Download
BIOMD0000000468-biopax3.owl Auto-generated BioPAX (Level 3) 176.30 KB Preview | Download
BIOMD0000000468.pdf Auto-generated PDF file 563.28 KB Preview | Download
BIOMD0000000468.m Auto-generated Octave file 30.31 KB Preview | Download
BIOMD0000000468.png Auto-generated Reaction graph (PNG) 1.65 MB Preview | Download
BIOMD0000000468.xpp Auto-generated XPP file 22.70 KB Preview | Download
Integratedmodel_noshearstress.cps Copasi file for the plot that corresponds to NO production without shear stress condition (green). Without shear stress implies that the shear activating reactions (re5, re57, re102 and re103) are removed from the model. 725.58 KB Preview | Download

  • Model originally submitted by : Vijayalakshmi Chelliah
  • Submitted: 19-Aug-2013 12:44:22
  • Last Modified: 21-Dec-2018 17:19:59
  • Version: 3 public model Download this version
    • Submitted on: 21-Dec-2018 17:19:59
    • Submitted by: administrator
    • With comment: Include the additional files provided by the submitter in the original submission: Integratedmodel_noshearstress.cps, Integratedmodel_shearstress.cps
  • Version: 2 public model Download this version
    • Submitted on: 22-May-2014 18:39:16
    • Submitted by: Vijayalakshmi Chelliah
    • With comment: Current version of Koo2013 - Integrated shear stress induced NO production model
  • Version: 1 public model Download this version
    • Submitted on: 19-Aug-2013 12:44:22
    • Submitted by: Vijayalakshmi Chelliah
    • With comment: Original import of Koo2013 - Integrated shear stress induced NO production model

(*) You might be seeing discontinuous revisions as only public revisions are displayed here. Any private revisions unpublished model revision of this model will only be shown to the submitter and their collaborators.

: Variable used inside SBML models

Reactions Rate Parameters
s48 => s47; s3, s3, s48, s47 k13*s3*s48-k14*s47 k14 = 1152.0 inv_sec; k13 = 0.08 inv_sec_sub
s5 => s2; s1, s2, s1 k_CCE*(fracK*Cao/(K3+Cao)-s2)*(s1-s2) fracK = 7071067.81 substance; Cao = 100.0 substance; K3 = 150.0 substance; k_CCE = 0.0 inv_sec_sub
s4 => s3; s3, s4 (-k6)*s3*(B_T-s4)+k7*s4 B_T = 120000.0 substance; k6 = 0.1 inv_sec_sub; k7 = 300.0 inv_sec
s3 => s9; s3 0.5*dot_Vhi*s3^4/(K_hi^4+s3^4) dot_Vhi = 4760.0 sub_sec; K_hi = 380.0 substance
s11 => s3 0.5*dot_q_inpass dot_q_inpass = 6000.0 sub_sec
s108 => s115; s107, s49, s107, s49 ktr1*s107+ktr1k2*s49 ktr1=1.2E-4 nM_inv_s; tr2=3.0E-6 nM_inv_s; ktr1k2=9.0E-6 nM_inv_s
s62 => s57 + s48; s62 kDD*s62 kDD = 9.45E-5 inv_sec
s105 => s49; s39, s39 exp(tau*(delay-s39/tc))/(1+2*exp(tau*(delay-s39/tc))+exp(2*tau*(delay-s39/tc)))/uc*29.256 uc=3600.0 dimensionless; tau=0.55 dimensionless; unity=20.0 nM_inv_s; tc=3600.0 substance; delay=5.0 dimensionless
s61 => s62; s27, s61, s27, s62 kp*s61*s27/(s61+Kmp)-Vdp*s62/(s62+Kmdp) Vdp = 4.0 sub_sec; kp = 0.1 inv_sec; Kmp = 5.0 substance; Kmdp = 20.0 substance
s62 => s65 + s48; s62, s65, s48 k18*s62-k17*s65*s48 k18 = 1.5 inv_sec; k17 = 1.5E-4 inv_sec_sub
s63 => s64; s45, s50, s61, s62, s58, s60, s65, s45, s58, s62, s65, s60 K_cam_no*(s45+s58)+K_p_no*(s62+s65)+K_pcam_no*s60 K_cam_no=17.0 inv_sec; K_p_no=5.0 inv_sec; K_pcam_no=17.0 inv_sec
s65 => s66; s65 Vdp*s65/(s65+Kmdp) Vdp = 4.0 sub_sec; Kmdp = 20.0 substance
s65 => s118 + s57; s65 kDD*s65 kDD = 9.45E-5 inv_sec
s95 => s93; s39, s119, s39 exp(unity-(s39/tf)^0.35)*normal*((s39+tiny_num)/unimol)^-0.65*(unity-(s39/tf)^0.35) unimol=1.0 substance; normal=4.0 nM_inv_s; tiny_num=1.0E-6 substance; tf=60.0 substance; unity=1.0 dimensionless
Curator's comment:
(added: 19 Aug 2013, 14:33:04, updated: 19 Aug 2013, 14:33:04)
The model is the integration of the all the sub-models (MODEL1302180003-6) described in the paper. Figure 6 of the paper is reproduced here. The model as such reproduces the plot that corresponds to NO production under shear-stress condition (red). To obtain the plot for NO production under no-shear-stress condition (green), shear activating reactions (re5, re57, re102 and re103) should be removed. The corresponding Copasi files can be downloaded below.