Proteolytic cascade in human blood coagulation system. Blood coagulation is a defense mechanism which prevent bleeding due to vessel damage by plug formation at the injury site. Any perturbation of the coagulation cascade can result in uncontrolled clotting or hemorrhage. The delicate regulatory system is controlled by a number of plasma protein-protein interactions as well as interaction with cellular receptors. Here, we present a primary blood coagulation model which is based on absolute minimal transformations involved in the process. The current model comprises of an extended set(108 in number) of molecular species and complexes. The model can be freely downloaded and used for further user specific investigation.
- Initial Parameters:
- -------------------
- M 1
- N 0
- kp1 1e1
- km1 0.02
- Parameters ends
- -------------------
- Initial Molecules
- -------------------
- II(loc~lip~flu,bs,state~u~a)
- V(loc~lip~flu,bs2,bs10,state~u~a)
- VII(loc~lip~flu,tf,bsf,state~u~a)
- VIII(loc~lip~flu,pl,bs9,bs10,state~u~a)
- X(loc~lip~flu,bs,state~u~a)
- IX(loc~lip~flu,bs,state~u~a)
- XI(loc~lip~flu,bs)
- TF(bs)
- AT(bs)
- LPD()
- Molecules ends
- -------------------
- Start Species
- -------------------
- II(loc~lip,at,bs,state~u) 1400
- V(loc~lip,bs2,bs10,state~u) 20
- VII(loc~lip,tf,bsf,state~a) 0.1
- VIII(loc~lip,pl,bs9,bs10,state~u) 0.7
- IX(loc~lip,at,bs,state~u) 90
- X(loc~lip,at,bs,state~u). 170
- XI(loc~lip,bs) 30
- TF(bs) 0.005
- AT(bs) 3400
- LPD() M
- Species ends
- -------------------
- Expanded Species (Model predicted)
- -------------------
- TF(bs!1).VII(bsf,loc~lip,state~a,tf!1) 0
- VII(bsf!1,loc~lip,state~a,tf).X(at,bs!1,loc~lip,state~u) 0par
- IX(at,bs!1,loc~lip,state~u).VII(bsf!1,loc~lip,state~a,tf) 0
- VIII(bs10,bs9,loc~lip,pl,state~a) 0
- V(bs10,bs2,loc~lip,state~a) 0
- TF(bs!1).VII(bsf!2,loc~lip,state~a,tf!1).X(at,bs!2,loc~lip,state~u) 0
- IX(at,bs!1,loc~lip,state~u).TF(bs!2).VII(bsf!1,loc~lip,state~a,tf!2) 0
- VII(bsf!1,loc~lip,state~a,tf).X(at,bs!1,loc~lip,state~a) 0
- IX(at,bs!1,loc~lip,state~a).VII(bsf!1,loc~lip,state~a,tf) 0
- TF(bs!1).VII(bsf!2,loc~lip,state~a,tf!1).X(at,bs!2,loc~lip,state~a) 0
- IX(at,bs!1,loc~lip,state~a).TF(bs!2).VII(bsf!1,loc~lip,state~a,tf!2) 0
- AT(bs!1).IX(at!1,bs!2,loc~lip,state~a).VII(bsf!2,loc~lip,state~a,tf) 0
- AT(bs!1).VII(bsf!2,loc~lip,state~a,tf).X(at!1,bs!2,loc~lip,state~a) 0
- X(at,bs,loc~lip,state~a) 0
- IX(at,bs,loc~lip,state~a) 0
- AT(bs!1).IX(at!1,bs!2,loc~lip,state~a).TF(bs!3).VII(bsf!2,loc~lip,state~a,tf!3) 0
- AT(bs!1).TF(bs!2).VII(bsf!3,loc~lip,state~a,tf!2).X(at!1,bs!3,loc~lip,state~a) 0
- AT(bs!1).IX(at!1,bs,loc~lip,state~a) 0
- AT(bs!1).X(at!1,bs,loc~lip,state~a) 0
- IX(at,bs!1,loc~lip,state~a).VIII(bs10,bs9!1,loc~lip,pl,state~a) 0
- V(bs10!1,bs2,loc~lip,state~a).X(at,bs!1,loc~lip,state~a) 0
- AT(bs!1).IX(at!1,bs!2,loc~lip,state~a).VIII(bs10,bs9!2,loc~lip,pl,state~a) 0
- AT(bs!1).V(bs10!2,bs2,loc~lip,state~a).X(at!1,bs!2,loc~lip,state~a) 0
- IX(at,bs!1,loc~lip,state~a).VIII(bs10!2,bs9!1,loc~lip,pl,state~a).X(at,bs!2,loc~lip,state~u) 0
- IX(at,bs!1,loc~lip,state~a).VIII(bs10!2,bs9!1,loc~lip,pl,state~a).X(at,bs!2,loc~lip,state~a) 0
- AT(bs!1).IX(at,bs!2,loc~lip,state~a).VIII(bs10!3,bs9!2,loc~lip,pl,state~a).X(at!1,bs!3,loc~lip,state~a) 0
- II(at,bs!1,loc~lip,state~u).V(bs10!2,bs2!1,loc~lip,state~a).X(at,bs!2,loc~lip,state~a) 0
- AT(bs!1).IX(at!1,bs!2,loc~lip,state~a).VIII(bs10!3,bs9!2,loc~lip,pl,state~a).X(at,bs!3,loc~lip,state~u) 0
- AT(bs!1).IX(at!1,bs!2,loc~lip,state~a).VIII(bs10!3,bs9!2,loc~lip,pl,state~a).X(at,bs!3,loc~lip,state~a) 0
- AT(bs!1).AT(bs!2).IX(at!1,bs!3,loc~lip,state~a).VIII(bs10!4,bs9!3,loc~lip,pl,state~a).X(at!2,bs!4,loc~lip,state~a) 0
- AT(bs!1).II(at,bs!2,loc~lip,state~u).V(bs10!3,bs2!2,loc~lip,state~a).X(at!1,bs!3,loc~lip,state~a) 0
- VIII(bs10!1,bs9,loc~lip,pl,state~a).X(at,bs!1,loc~lip,state~u) 0
- VIII(bs10!1,bs9,loc~lip,pl,state~a).X(at,bs!1,loc~lip,state~a) 0
- AT(bs!1).VIII(bs10!2,bs9,loc~lip,pl,state~a).X(at!1,bs!2,loc~lip,state~a) 0
- II(at,bs!1,loc~lip,state~a).V(bs10!2,bs2!1,loc~lip,state~a).X(at,bs!2,loc~lip,state~a) 0
- AT(bs!1).II(at!1,bs!2,loc~lip,state~a).V(bs10!3,bs2!2,loc~lip,state~a).X(at,bs!3,loc~lip,state~a) 0
- AT(bs!1).II(at,bs!2,loc~lip,state~a).V(bs10!3,bs2!2,loc~lip,state~a).X(at!1,bs!3,loc~lip,state~a) 0
- AT(bs!1).AT(bs!2).II(at!1,bs!3,loc~lip,state~a).V(bs10!4,bs2!3,loc~lip,state~a).X(at!2,bs!4,loc~lip,state~a) 0
- Expanded Species ends
- -------------------
- Rules
- -------------------
- II(bs,loc~flu) + LPD <-> II(bs,loc~lip) kp1,km1
- IX(bs,loc~flu) + LPD <-> IX(bs,loc~lpd) kp1,km1
- X(bs,loc~flu) + LPD <-> X(bs,loc~lip) kp1,km1
- V(bs2,bs10,loc~flu) + LPD <-> V(bs2,bs10,loc~lip) kp1,km1
- VII(bsf,loc~flu) + LPD <-> VII(bsf,loc~lpd) kp1,km1
- VIII(bs9,bs10,,loc~flu) + LPD <-> VIII(bs9,bs10,loc~lip) kp1,km1
- TF(bs) + VII(tf) <-> TF(bs!1).VII(tf!1) kp1,km1
- AT(bs) + II(at,state~a) <-> AT(bs!1).II(at!1,state~a) kp1,km1
- AT(bs) + IX(at,state~a) <-> AT(bs!1).IX(at!1,state~a) kp1,km1
- AT(bs) + X(at,state~a) <-> AT(bs!1).X(at!1,state~a) kp1,km1
- VII(state~a,bsf,loc~flu) + X(bs,loc~flu) <-> VII(state~a,bsf!1,loc~flu).X(bs!1,loc~flu) kp1,km1
- VII(state~a,bsf,loc~lip) + X(bs,loc~lip) <-> VII(state~a,bsf!1,loc~lip).X(bs!1,loc~lip) kp1,km1
- VII(state~a,bsf!1).X(bs!1,state~u) -> VII(state~a,bsf!1).X(bs!1,state~a) kp1
- VII(state~a,bsf,loc~flu) + IX(bs,loc~flu) <-> VII(state~a,bsf!1,loc~flu).IX(bs!1,loc~flu) kp1,km1
- VII(state~a,bsf,loc~lip) + IX(bs,loc~lip) <-> VII(state~a,bsf!1,loc~lip).IX(bs!1,loc~lip) kp1,km1
- VII(state~a,bsf!1).IX(bs!1,state~u) -> VII(state~a,bsf!1).IX(bs!1,state~a) kp1
- VIII(state~u) -> VIII(state~a) kp1
- VIII(state~a,bs9,loc~lip) + IX(bs,state~a,loc~lip) <-> VIII(state~a,bs9!1,loc~lip).IX(bs!1,state~a,loc~lip) kp1,km1
- VIII(state~a,bs9,loc~flu) + IX(bs,state~a,loc~flu) <-> VIII(state~a,bs9!1,loc~flu).IX(bs!1,state~a,loc~flu) kp1,km1
- VIII(state~a,bs10,loc~lip) + X(bs,state~a,loc~lip) <-> VIII(state~a,bs10!1,loc~lip).X(bs!1,state~a,loc~lip) kp1,km1
- VIII(state~a,bs10,loc~flu) + X(bs,state~a,loc~flu) <-> VIII(state~a,bs10!1,loc~flu).X(bs!1,state~a,loc~flu) kp1,km1
- VIII(state~a,bs9!2,bs10!1).X(bs!1,state~u).IX(state~a,bs!2) -> VIII(state~a,bs9!2,bs10!1).X(bs!1,state~a).IX(state~a,bs!2) kp1
- V(state~u) -> V(state~a) kp1
- V(state~a,bs10,loc~lip) + X(state~a,bs,loc~lip) -> V(state~a,bs10!1,loc~lip).X(state~a,bs!1,loc~lip) kp1
- V(state~a,bs10,loc~flu) + X(state~a,bs,loc~flu) -> V(state~a,bs10!1,loc~flu).X(state~a,bs!1,loc~flu) kp1
- V(state~a,bs2,bs10!1,loc~lip).X(state~a,bs!1,loc~lip) + II(bs,loc~lip) -> V(state~a,bs2!2,bs10!1,loc~lip).X(state~a,bs!1,loc~lip).II(bs!2,loc~lip) kp1
- V(state~a,bs2!2,bs10!1).X(state~a,bs!1).II(bs!2,state~u) -> V(state~a,bs2!2,bs10!1).X(state~a,bs!1).II(bs!2,state~a) kp1
- Rules ends
- -------------------
The current model for coagulation cascade encompasses the basic events involved in thrombin generation and various coagulation factors transformations. The user has the ability to download our quantitative rule based model from the PathwayDB pages for further modifications and extensions.
System requirements
Model format: BNGL
Analysis system: BioNetGen standalone or VirtualCell enviornment http://vcell.org/bionetgen.
Rule-based model file has the following structure
Each block in the model file starts with "begin" and ends with "end" keywords. There are mainly four blocks in BNGL file namely, 1) Parameters 2) Species 3) Reaction rules and 4) Observables.
Parameters Block
Parameters refer to numerical values for the variables that are associated with species(eg. molecular concentrations, M) and rules (eg. reaction rates, k).
parameter_index/tab/parameter_name/tab/parameter_value
Index | Params | Values |
|---|---|---|
1 | M | 1 |
2 | kp1 | 1e+1 |
Species Block
Species block refers to Seed species i.e. the species which are defined at the start of network generation. Each species which exist in different states(active/inactive) must be defined separately in the Species block. Each species name is defined as a structured object with all the associated components for example, binding sites, state of activity and localization information.
species_name(component1, component2) Molecular_concentration
II(loc~lip,at,bs,state~u) M
Observables Block
Observable block refers to the molecules and species that the user wish to map using the model understudy. Each observable is associated with the corresponding components.
Molecules/space/Molecule_species_name/space/species_components
Molecules V V()
Molecules XXX VII(state~a,bsf!1).IX(bs!1,state~a)
GDAT file consists of timecourse generated for each observable.
Rules Block
Rules are the basic transformations that are needed for generating species and reactions in the rule-based model. Each line of a rule has reactant(s), an arrow(->,<->), product(s) and the reaction rates.
By definition, each rule generates single or multiple species(expanded species) and reactions.
index/tab/Reactant1 + Reactant2/space/ <->/space/Product1 + Product2/space/rates
Rule 1: II(bs,loc~flu) + LPD <-> II(bs,loc~lip) kp1,km1
Rule 2: IX(bs,loc~flu) + LPD <-> IX(bs,loc~lpd) kp1,km1
These four blocks are the basic framework for defining a rule-based model. Once, the user defines the model, there are following commands which can be used to run the model:
- CutDB | PathwayDB | ProteaseDB | SubstrateDB | ProfileDB | Tools
- Literature Track for CutDB | PMAP | CPP | User Forum | Burnham Institute for Medical Research
- E-mail: pmap_info at burnham dot org
- System last updated Jun.24.2009