From 2009.igem.org
model = changeRxnBounds(model,{'EX_glc(e)','EX_gly(e)','EX_ala_D(e)','EX_ala_L(e)','EX_arg_L(e)','EX_asn_L(e)',
'EX_asp_L(e)','EX_cys_D(e)','EX_cys_L(e)','EX_gln_L(e)','EX_glu_L(e)','EX_his_L(e)','EX_ile_L(e)',
'EX_leu_L(e)','EX_lys_L(e)','EX_met_D(e)','EX_met_L(e)','EX_phe_L(e)','EX_pro_L(e)','EX_ser_D(e)',
'EX_ser_L(e)','EX_thr_L(e)','EX_trp_L(e)','EX_tyr_L(e)','EX_val_L(e)','EX_adn(e)','EX_gsn(e)',
'EX_cytd(e)','EX_thymd(e)','EX_uri(e)','EX_thm(e)','EX_pnto_R(e)','EX_xyl_D(e)','EX_skm(e)',
'EX_26dap_M(e)'},[-99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99
-99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99],'l');
solutionM=optimizeCbModel(model);
InitialGrowthRate = solution.f
"
