ExampleNetwork.cpp Source File

00001 // ExampleNetwork.cpp: implementation of the CExampleNetwork class.
00002 //
00004 
00005 #include "ExampleNetwork.h"
00006 
00008 // Construction/Destruction
00010 
00011 CExampleNetwork::CExampleNetwork()
00012 {
00013         double CFactor = 600000;
00014         double SecPerMin = 60.0;
00015 
00016         Chemical *Ligand = new Chemical(chemicals.size(),0.0,"Ligand");
00017         chemicals.push_back(Ligand);
00018 
00019         Chemical *FreeReceptor = new Chemical(chemicals.size(),1000,"FreeReceptor");
00020         chemicals.push_back(FreeReceptor);
00021 
00022         Chemical *BoundReceptor = new Chemical(chemicals.size(),0.0,"BoundReceptor");
00023         chemicals.push_back(BoundReceptor);
00024 
00025         Chemical *InactiveGEF = new Chemical(chemicals.size(),10000,"InactiveGEF");
00026         chemicals.push_back(InactiveGEF);
00027 
00028         Chemical *ActiveGEF = new Chemical(chemicals.size(),0.0,"ActiveGEF");
00029         chemicals.push_back(ActiveGEF);
00030 
00031         Chemical *SmallGGDP = new Chemical(chemicals.size(),10000,"SmallGGDP");
00032         chemicals.push_back(SmallGGDP);
00033 
00034         Chemical *SmallGGTP = new Chemical(chemicals.size(),0.0,"SmallGGTP");
00035         chemicals.push_back(SmallGGTP);
00036 
00037         RateConstant *k_b = new RateConstant(0.1*SecPerMin,"k_b");
00038         rateConstants.push_back(k_b);
00039 
00040         RateConstant *k_GEF = new RateConstant(0.05*SecPerMin,"k_GEF");
00041         rateConstants.push_back(k_GEF);
00042 
00043         RateConstant *Km_GEF = new RateConstant(CFactor*0.5,"Km_GEF");
00044         rateConstants.push_back(Km_GEF);
00045 
00046         RateConstant *k_G = new RateConstant(0.075*SecPerMin,"k_G");
00047         rateConstants.push_back(k_G);
00048 
00049         RateConstant *Km_G = new RateConstant(CFactor*0.25,"Km_G");
00050         rateConstants.push_back(Km_G);
00051 
00052         Reaction *LigandBinding = new HeterodimerizationReaction(Ligand,FreeReceptor,BoundReceptor,k_b,"LigandBinding");
00053         reactions.push_back(LigandBinding);
00054 
00055         Reaction *GEFActivation = new MichaelisMentenPhosphorylationReaction(BoundReceptor,InactiveGEF,ActiveGEF,k_GEF,Km_GEF,"GEFActivation");
00056         reactions.push_back(GEFActivation);
00057 
00058         Reaction *SmallGActivation = new MichaelisMentenGProteinActivationReaction(ActiveGEF,SmallGGDP,SmallGGTP,k_G,Km_G,"GProteinActivation");
00059         reactions.push_back(SmallGActivation);
00060 }
00061 
00062 CExampleNetwork::~CExampleNetwork()
00063 {
00064 
00065 }