SmallErkSubnetwork.cpp Source File

00001 // ErkSubnetwork.cpp: implementation of the ErkSubnetwork class.
00002 //
00004 
00005 #ifdef _WINDOWS
00006 #include "stdafx.h"
00007 #endif
00008 
00009 #include "SmallErkSubnetwork.h"
00010 
00012 // Construction/Destruction
00014 
00015 // Changed stoichiometry on 3/8/01
00016 // - ERK and MEK are at CFactor
00017 // - NGFR = 10000, EGFR = 80000
00018 // - all else at 0.2*CFactor
00019 
00020 SmallErkSubnetwork::SmallErkSubnetwork()
00021 {
00022         // CFactor = #molecules/1uM
00023         double CFactor = 600000.0, SecToMin = 60.0;
00024 
00025         // 1
00026         Chemical *EGF = new Chemical(chemicals.size(), 0.0, "EGF");
00027         chemicals.push_back(EGF);
00028 
00029         // 2
00030         Chemical *NGF = new Chemical(chemicals.size(), 0.0, "NGF");
00031         chemicals.push_back(NGF);
00032 
00033         // 3 (ref - constraint on [EGFR]_T)
00034         Chemical *freeEGFReceptor = new Chemical(chemicals.size(), 80000, "freeEGFReceptor");
00035         chemicals.push_back(freeEGFReceptor);
00036 
00037         // 4 (ref)
00038         Chemical *boundEGFReceptor = new Chemical(chemicals.size(), 0.0, "boundEGFReceptor");
00039         chemicals.push_back(boundEGFReceptor);
00040 
00041         // 5 (ref)
00042         Chemical *freeNGFReceptor = new Chemical(chemicals.size(), 10000, "freeNGFReceptor");
00043         chemicals.push_back(freeNGFReceptor);
00044 
00045         // 6
00046         Chemical *boundNGFReceptor = new Chemical(chemicals.size(), 0.0, "boundNGFReceptor");
00047         chemicals.push_back(boundNGFReceptor);
00048 
00049         // 7
00050         Chemical *SosInactive = new Chemical(chemicals.size(), 0.2*CFactor, "SosInactive");
00051         chemicals.push_back(SosInactive);
00052                 
00053         // 8
00054         Chemical *SosActive = new Chemical(chemicals.size(), 0.0, "SosActive");
00055         chemicals.push_back(SosActive);
00056 
00057         // 9
00058         Chemical *P90RskInactive = new Chemical(chemicals.size(), 0.2*CFactor, "P90RskInactive");
00059         chemicals.push_back(P90RskInactive);
00060 
00061         // 10
00062         Chemical *P90RskActive = new Chemical(chemicals.size(), 0.0, "P90RskActive");
00063         chemicals.push_back(P90RskActive);
00064 
00065         // 11
00066         Chemical *RasInactive = new Chemical(chemicals.size(), 0.2*CFactor,"RasInactive");
00067         chemicals.push_back(RasInactive);
00068 
00069         // 12
00070         Chemical *RasActive = new Chemical(chemicals.size(), 0.0, "RasActive");
00071         chemicals.push_back(RasActive);
00072 
00073         // 13
00074         Chemical *RasGapActive = new Chemical(chemicals.size(), CFactor*0.2, "RasGapActive");
00075         chemicals.push_back(RasGapActive);
00076 
00077         
00078         // 14 (ref - constraint on [Raf1]_T = CFactor*0.003)
00079         Chemical *Raf1Inactive = new Chemical(chemicals.size(), CFactor*0.2, "Raf1Inactive");
00080         chemicals.push_back(Raf1Inactive);
00081                 
00082         // 15 (ref)
00083         Chemical *Raf1Active = new Chemical(chemicals.size(), 0.0, "Raf1Active");
00084         chemicals.push_back(Raf1Active);
00085 
00086         // 16 (ref- constraint on [BRaf]_T = CFactor*0.003)
00087         Chemical *BRafInactive = new Chemical(chemicals.size(), CFactor*0.2, "BRafInactive");
00088         chemicals.push_back(BRafInactive);
00089                 
00090         // 17 (ref)
00091         Chemical *BRafActive = new Chemical(chemicals.size(), 0.0, "BRafActive");
00092         chemicals.push_back(BRafActive);
00093 
00094         // 18 (ref - constraint on [Mek]_T)
00095         Chemical *MekInactive = new Chemical(chemicals.size(), CFactor, "MekInactive");
00096         chemicals.push_back(MekInactive);
00097 
00098         // 19 (ref)
00099         Chemical *MekActive = new Chemical(chemicals.size(), 0.0, "MekActive");
00100         chemicals.push_back(MekActive);
00101 
00102         // 20 (ref - constraint on [Erk]_T)
00103         Chemical *ErkInactive = new Chemical(chemicals.size(), CFactor,"ErkInactive");
00104         chemicals.push_back(ErkInactive);
00105 
00106         // 21 (ref)
00107         Chemical *ErkActive = new Chemical(chemicals.size(), 0.0, "ErkActive");
00108         chemicals.push_back(ErkActive);
00109         
00110         
00111         // 22
00112         Chemical *PI3KInactive = new Chemical(chemicals.size(), 0.2*CFactor, "PI3KInactive");
00113         chemicals.push_back(PI3KInactive);
00114 
00115         // 23
00116         Chemical *PI3KActive = new Chemical(chemicals.size(), 0.0, "PI3KActive");
00117         chemicals.push_back(PI3KActive);
00118 
00119         // 24
00120         Chemical *AktInactive = new Chemical(chemicals.size(), 0.2*CFactor, "AktInactive");
00121         chemicals.push_back(AktInactive);
00122 
00123         // 25
00124         Chemical *AktActive = new Chemical(chemicals.size(), 0.0, "AktActive");
00125         chemicals.push_back(AktActive);
00126         
00127 
00128         // 26
00129         Chemical *C3GInactive = new Chemical(chemicals.size(), 0.2*CFactor, "C3GInactive");
00130         chemicals.push_back(C3GInactive);
00131                 
00132         // 27
00133         Chemical *C3GActive = new Chemical(chemicals.size(), 0.0, "C3GActive");
00134         chemicals.push_back(C3GActive);
00135 
00136         // 28
00137         Chemical *Rap1Inactive = new Chemical(chemicals.size(), 0.2*CFactor, "Rap1Inactive");
00138         chemicals.push_back(Rap1Inactive);
00139 
00140         // 29
00141         Chemical *Rap1Active = new Chemical(chemicals.size(), 0.0, "Rap1Active");
00142         chemicals.push_back(Rap1Active);
00143 
00144         // 30
00145         Chemical *RapGapActive = new Chemical(chemicals.size(), 0.2*CFactor, "RapGapActive");
00146         chemicals.push_back(RapGapActive);
00147         
00148         // 31
00149         Chemical *PP2AActive = new Chemical(chemicals.size(), 0.2*CFactor, "PP2AActive");
00150         chemicals.push_back(PP2AActive);
00151 
00152         // 32 
00153         Chemical *Raf1PPtase = new Chemical(chemicals.size(), 0.2*CFactor, "Raf1PPtase");
00154         chemicals.push_back(Raf1PPtase);
00155 
00156         /*
00157         // 33 
00158         Chemical *PIP2 = new Chemical(chemicals.size(), CFactor, "PIP2");
00159         chemicals.push_back(PIP2);
00160 
00161         // 34 
00162         Chemical *PIP3 = new Chemical(chemicals.size(), 0.0, "PIP3");
00163         chemicals.push_back(PIP3);
00164 
00165         // 35
00166         Chemical *PTEN = new Chemical(chemicals.size(), 0.2*CFactor, "PTEN");
00167         chemicals.push_back(PTEN);
00168         */
00169 
00170         RateConstant *k_rb_EGF = new RateConstant(SecToMin*0.00001, "krbEGF");
00171         rateConstants.push_back(k_rb_EGF);
00172         RateConstant *k_ru_EGF = new RateConstant(SecToMin*0.01, "kruEGF");
00173         rateConstants.push_back(k_ru_EGF);
00174         RateConstant *k_rb_NGF = new RateConstant(SecToMin*0.00005, "krbNGF");
00175         rateConstants.push_back(k_rb_NGF);
00176         RateConstant *k_ru_NGF = new RateConstant(SecToMin*0.02, "kruNGF");
00177         rateConstants.push_back(k_ru_NGF);
00178 
00179         RateConstant *k_EGF = new RateConstant(SecToMin*0.2, "kEGF");
00180         rateConstants.push_back(k_EGF);
00181         RateConstant *Km_EGF = new RateConstant(CFactor*0.8333, "KmEGF");
00182         rateConstants.push_back(Km_EGF);
00183         RateConstant *k_NGF = new RateConstant(SecToMin*0.4, "kNGF");
00184         rateConstants.push_back(k_NGF);
00185         RateConstant *Km_NGF = new RateConstant(CFactor*0.8333, "KmNGF");
00186         rateConstants.push_back(Km_NGF);
00187 
00188         RateConstant *k_dSos = new RateConstant(0.01*SecToMin, "kdSos");
00189         rateConstants.push_back(k_dSos);
00190         RateConstant *Km_dSos = new RateConstant(100000.0, "KmdSos");
00191         rateConstants.push_back(Km_dSos);       
00192         RateConstant *k_Sos = new RateConstant(0.075*SecToMin, "kSos");
00193         rateConstants.push_back(k_Sos);
00194         RateConstant *Km_Sos = new RateConstant(0.25*CFactor, "KmSos");
00195         rateConstants.push_back(Km_Sos);
00196 
00197         RateConstant *k_RasGap = new RateConstant(0.2*SecToMin, "kRasGap");
00198         rateConstants.push_back(k_RasGap);
00199         RateConstant *Km_RasGap = new RateConstant(CFactor*1.0104, "KmRasGap");
00200         rateConstants.push_back(Km_RasGap);
00201         RateConstant *k_RasToRaf1 = new RateConstant(0.08*SecToMin, "kRasToRaf1");
00202         rateConstants.push_back(k_RasToRaf1);
00203         RateConstant *Km_RasToRaf1 = new RateConstant(150000.0, "KmRasToRaf1");
00204         rateConstants.push_back(Km_RasToRaf1);
00205         
00206         RateConstant *k_pRaf1 = new RateConstant(0.05*SecToMin, "kpRaf1");
00207         rateConstants.push_back(k_pRaf1);
00208         // (ref)
00209         RateConstant *Km_pRaf1 = new RateConstant(0.6*CFactor, "KmpRaf1");
00210         rateConstants.push_back(Km_pRaf1);
00211         RateConstant *k_pBRaf = new RateConstant(0.025*SecToMin, "kpBRaf");
00212         rateConstants.push_back(k_pBRaf);
00213         RateConstant *Km_pBRaf = new RateConstant(100000.0, "KmpBRaf");
00214         rateConstants.push_back(Km_pBRaf);
00215 
00216         RateConstant *k_dMek = new RateConstant(0.03*SecToMin, "kdMek");
00217         rateConstants.push_back(k_dMek);
00218         RateConstant *Km_dMek = new RateConstant(125000.0, "KmdMek");
00219         rateConstants.push_back(Km_dMek);
00220         
00221         RateConstant *k_pMek_Cytoplasmic = new RateConstant(0.01*SecToMin, "kpMekCytoplasmic");
00222         rateConstants.push_back(k_pMek_Cytoplasmic);
00223         // (ref)
00224         RateConstant *Km_pMek_Cytoplasmic = new RateConstant(125000.0, "KmpMekCytoplasmic");
00225         rateConstants.push_back(Km_pMek_Cytoplasmic);
00226         RateConstant *k_dErk = new RateConstant(0.01*SecToMin, "kdErk");
00227         rateConstants.push_back(k_dErk);
00228         RateConstant *Km_dErk = new RateConstant(300000.0, "KmdErk");
00229         rateConstants.push_back(Km_dErk);
00230         RateConstant *k_pP90Rsk = new RateConstant(0.01*SecToMin, "kpP90Rsk");
00231         rateConstants.push_back(k_pP90Rsk);
00232         RateConstant *Km_pP90Rsk = new RateConstant(300000.0, "KmpP90Rsk");
00233         rateConstants.push_back(Km_pP90Rsk);
00234                 
00235         
00236         RateConstant *k_PI3K = new RateConstant(0.05*SecToMin, "kPI3K");
00237         rateConstants.push_back(k_PI3K);
00238         RateConstant *Km_PI3K = new RateConstant(100000.0, "KmPI3K");
00239         rateConstants.push_back(Km_PI3K);
00240         RateConstant *k_PI3K_Ras = new RateConstant(0.025*SecToMin, "kPI3KRas");
00241         rateConstants.push_back(k_PI3K_Ras);
00242         RateConstant *Km_PI3K_Ras = new RateConstant(300000.0, "KmPI3KRas");
00243         rateConstants.push_back(Km_PI3K_Ras);
00244         RateConstant *k_Akt = new RateConstant(0.01*SecToMin, "kAkt");
00245         rateConstants.push_back(k_Akt);
00246         RateConstant *Km_Akt = new RateConstant(300000.0, "KmAkt");
00247         rateConstants.push_back(Km_Akt);
00248         RateConstant *k_dRaf1ByAkt = new RateConstant(0.1*SecToMin, "kdRaf1ByAkt");
00249         rateConstants.push_back(k_dRaf1ByAkt);
00250         RateConstant *Km_Raf1ByAkt = new RateConstant(300000.0, "KmRaf1ByAkt"); 
00251         rateConstants.push_back(Km_Raf1ByAkt);
00252         
00253         
00254         RateConstant *k_C3G_NGF = new RateConstant(0.05*SecToMin, "kC3GNGF");
00255         rateConstants.push_back(k_C3G_NGF);
00256         RateConstant *Km_C3G_NGF = new RateConstant(200000.0, "KmC3GNGF");
00257         rateConstants.push_back(Km_C3G_NGF);
00258         RateConstant *k_C3G = new RateConstant(0.01*SecToMin, "kC3G"); 
00259         rateConstants.push_back(k_C3G);
00260         RateConstant *Km_C3G = new RateConstant(300000.0, "KmC3G");
00261         rateConstants.push_back(Km_C3G);
00262         RateConstant *k_RapGap = new RateConstant(0.002*SecToMin, "kRapGap");
00263         rateConstants.push_back(k_RapGap);
00264         RateConstant *Km_RapGap = new RateConstant(300000.0, "KmRapGap");
00265         rateConstants.push_back(Km_RapGap);
00266         RateConstant *k_Rap1ToBRaf = new RateConstant(0.05*SecToMin, "kRap1ToBRaf");
00267         rateConstants.push_back(k_Rap1ToBRaf);
00268         RateConstant *Km_Rap1ToBRaf = new RateConstant(300000.0, "KmRap1ToBRaf");
00269         rateConstants.push_back(Km_Rap1ToBRaf);
00270 
00271         RateConstant *k_dRaf1 = new RateConstant(0.1*SecToMin, "kdRaf1");
00272         rateConstants.push_back(k_dRaf1);
00273         RateConstant *Km_dRaf1 = new RateConstant(300000.0, "KmdRaf1");
00274         rateConstants.push_back(Km_dRaf1);
00275 
00276         RateConstant *k_dBRaf = new RateConstant(0.05*SecToMin, "kdBRaf");
00277         rateConstants.push_back(k_dBRaf);
00278         RateConstant *Km_dBRaf = new RateConstant(300000.0, "KmdBRaf");
00279         rateConstants.push_back(Km_dBRaf);
00280 
00281         /*
00282         RateConstant *k_PIP2 = new RateConstant(SecToMin*2, "k_PIP2");
00283         rateConstants.push_back(k_PIP2);
00284         RateConstant *Km_PIP2 = new RateConstant(0.25*CFactor,"Km_PIP2");
00285         rateConstants.push_back(Km_PIP2);
00286         RateConstant *k_dPIP3 = new RateConstant(SecToMin, "k_PIP3");
00287         rateConstants.push_back(k_dPIP3);
00288         RateConstant *Km_dPIP3 = new RateConstant(0.25*CFactor,"Km_dPIP3");
00289         rateConstants.push_back(Km_dPIP3);
00290         */
00291 
00292         ReceptorBindingReaction *EGFBindingReaction = new ReceptorBindingReaction(EGF, freeEGFReceptor, boundEGFReceptor, k_rb_EGF, "EGFBindingReaction");
00293         reactions.push_back(EGFBindingReaction);
00294         ReceptorUnbindingReaction *EGFUnbindingReaction = new ReceptorUnbindingReaction(boundEGFReceptor, freeEGFReceptor, EGF, k_ru_EGF, "EGFUnbindingReaction");
00295         reactions.push_back(EGFUnbindingReaction);
00296         ReceptorBindingReaction *NGFBindingReaction = new ReceptorBindingReaction(NGF, freeNGFReceptor, boundNGFReceptor, k_rb_NGF, "NGFBindingReaction");
00297         reactions.push_back(NGFBindingReaction);
00298         ReceptorUnbindingReaction *NGFUnbindingReaction = new ReceptorUnbindingReaction(boundNGFReceptor, freeNGFReceptor, NGF, k_ru_NGF, "NGFUnbindingReaction");
00299         reactions.push_back(NGFUnbindingReaction);
00300 
00301         MichaelisMentenPhosphorylationReaction *SosActivationByEGFReaction = new MichaelisMentenPhosphorylationReaction(boundEGFReceptor, SosInactive, SosActive, k_EGF, Km_EGF, "SosActivationbyEGFReaction");
00302         reactions.push_back(SosActivationByEGFReaction);
00303         MichaelisMentenPhosphorylationReaction *SosActivationByNGFReaction = new MichaelisMentenPhosphorylationReaction(boundNGFReceptor, SosInactive, SosActive, k_NGF, Km_NGF, "SosActivationbyNGFReaction");
00304         reactions.push_back(SosActivationByNGFReaction);
00305         MichaelisMentenDephosphorylationReaction *SosDeactivationReaction = new MichaelisMentenDephosphorylationReaction(P90RskActive, SosActive, SosInactive, k_dSos, Km_dSos, "SosDeactivationReaction");
00306         reactions.push_back(SosDeactivationReaction);
00307         MichaelisMentenGProteinActivationReaction *RasActivationReaction = new MichaelisMentenGProteinActivationReaction(SosActive, RasInactive, RasActive, k_Sos, Km_Sos, "RasActivationReaction");
00308         reactions.push_back(RasActivationReaction);
00309         
00310         MichaelisMentenGProteinDeactivationReaction *RasDeactivationReaction = new MichaelisMentenGProteinDeactivationReaction(RasGapActive, RasActive, RasInactive, k_RasGap, Km_RasGap, "RasDeactivationReaction");
00311         reactions.push_back(RasDeactivationReaction);
00312         MichaelisMentenProteinProteinActivationReaction *Raf1ByRasActivationReaction = new MichaelisMentenProteinProteinActivationReaction(RasActive, Raf1Inactive, Raf1Active, k_RasToRaf1, Km_RasToRaf1, "Raf1ByRasActivationReaction");
00313         reactions.push_back(Raf1ByRasActivationReaction);
00314 
00315         MichaelisMentenPhosphorylationReaction *MekbyRaf1ActivationReaction = new MichaelisMentenPhosphorylationReaction(Raf1Active, MekInactive, MekActive, k_pRaf1, Km_pRaf1, "MekbyRaf1ActivationReaction");
00316         reactions.push_back(MekbyRaf1ActivationReaction);
00317         MichaelisMentenPhosphorylationReaction *MekbyBRafActivationReaction = new MichaelisMentenPhosphorylationReaction(BRafActive, MekInactive, MekActive, k_pBRaf, Km_pBRaf, "MekbyBRafActivationReaction");
00318         reactions.push_back(MekbyBRafActivationReaction);
00319         MichaelisMentenPhosphorylationReaction *ErkActivationReaction = new MichaelisMentenPhosphorylationReaction(MekActive, ErkInactive, ErkActive, k_pMek_Cytoplasmic, Km_pMek_Cytoplasmic, "ErkActivationReaction");
00320         reactions.push_back(ErkActivationReaction);
00321 
00322         
00323         MichaelisMentenDephosphorylationReaction *MekDeactivationReaction = new MichaelisMentenDephosphorylationReaction(PP2AActive, MekActive, MekInactive, k_dMek, Km_dMek, "MekDeactivationReaction");
00324         reactions.push_back(MekDeactivationReaction);
00325         MichaelisMentenDephosphorylationReaction *ErkDeactivationReaction = new MichaelisMentenDephosphorylationReaction(PP2AActive, ErkActive, ErkInactive, k_dErk, Km_dErk, "ErkDeactivationReaction");
00326         reactions.push_back(ErkDeactivationReaction);
00327 
00328         MichaelisMentenDephosphorylationReaction *Raf1byPPtaseDeactivationReaction = new MichaelisMentenDephosphorylationReaction(Raf1PPtase, Raf1Active, Raf1Inactive, k_dRaf1, Km_dRaf1, "Raf1byPPtaseDeactivationReaction");
00329         reactions.push_back(Raf1byPPtaseDeactivationReaction);
00330 
00331         MichaelisMentenPhosphorylationReaction *P90RskActivationReaction = new MichaelisMentenPhosphorylationReaction(ErkActive, P90RskInactive, P90RskActive, k_pP90Rsk, Km_pP90Rsk, "P90RskActivationReaction");
00332         reactions.push_back(P90RskActivationReaction);
00333         
00334         
00335         MichaelisMentenPhosphorylationReaction *PI3KbyEGFRActivationReaction = new MichaelisMentenPhosphorylationReaction(boundEGFReceptor, PI3KInactive, PI3KActive, k_PI3K, Km_PI3K, "PI3KbyEGFRActivationReaction");
00336         reactions.push_back(PI3KbyEGFRActivationReaction);
00337         MichaelisMentenProteinProteinActivationReaction *PI3KbyRasActivationReaction = new MichaelisMentenProteinProteinActivationReaction(RasActive, PI3KInactive, PI3KActive, k_PI3K_Ras, Km_PI3K_Ras, "PI3KbyRasActivationReaction");
00338         reactions.push_back(PI3KbyRasActivationReaction);
00339         
00340         /*
00341         MichaelisMentenPhosphorylationReaction *PIP2PIP3Conversion = new MichaelisMentenPhosphorylationReaction(PI3KActive,PIP2,PIP3,k_PIP2,Km_PIP2,"PIP2PIP3Conversion");
00342         reactions.push_back(PIP2PIP3Conversion);
00343         MichaelisMentenDephosphorylationReaction *PIP3PIP2Conversion = new MichaelisMentenDephosphorylationReaction(PTEN,PIP3,PIP2,k_dPIP3,Km_dPIP3,"PIP3PIP2Conversion");
00344         reactions.push_back(PIP3PIP2Conversion);
00345         */
00346         
00347         MichaelisMentenPhosphorylationReaction *AktActivationReaction = new MichaelisMentenPhosphorylationReaction(PI3KActive, AktInactive, AktActive, k_Akt, Km_Akt, "AktActivationReaction");
00348         reactions.push_back(AktActivationReaction);
00349         MichaelisMentenProteinProteinDeactivationReaction *Raf1ByAktDeactivationReaction = new MichaelisMentenProteinProteinDeactivationReaction(AktActive, Raf1Active, Raf1Inactive, k_dRaf1ByAkt, Km_Raf1ByAkt, "RafByAktDeactivationReaction");
00350         reactions.push_back(Raf1ByAktDeactivationReaction);     
00351         
00352         MichaelisMentenPhosphorylationReaction *C3GActivationReaction = new MichaelisMentenPhosphorylationReaction(boundNGFReceptor, C3GInactive, C3GActive, k_C3G_NGF, Km_C3G_NGF, "C3GActivationReaction");
00353         reactions.push_back(C3GActivationReaction);
00354         MichaelisMentenGProteinActivationReaction *Rap1ActivationReaction = new MichaelisMentenGProteinActivationReaction(C3GActive, Rap1Inactive, Rap1Active, k_C3G, Km_C3G, "Rap1ActivationReaction");
00355         reactions.push_back(Rap1ActivationReaction);
00356         MichaelisMentenGProteinDeactivationReaction *Rap1DeactivationReaction = new MichaelisMentenGProteinDeactivationReaction(RapGapActive, Rap1Active, Rap1Inactive, k_RapGap, Km_RapGap, "Rap1DeactivationReaction");
00357         reactions.push_back(Rap1DeactivationReaction);
00358         MichaelisMentenProteinProteinActivationReaction *BRafByRap1ActivationReaction = new MichaelisMentenProteinProteinActivationReaction(Rap1Active, BRafInactive, BRafActive, k_Rap1ToBRaf, Km_Rap1ToBRaf, "BRafByRap1ActivationReaction");
00359         reactions.push_back(BRafByRap1ActivationReaction);
00360 
00361         MichaelisMentenDephosphorylationReaction *BRafbyPPtaseDeactivationReaction = new MichaelisMentenDephosphorylationReaction(Raf1PPtase, BRafActive, BRafInactive, k_dBRaf, Km_dBRaf, "BRafbyPPtaseDeactivationReaction");
00362         reactions.push_back(BRafbyPPtaseDeactivationReaction);
00363 }
00364 
00365 SmallErkSubnetwork::~SmallErkSubnetwork()
00366 {
00367 
00368 }