大肠杆菌 (E. coli) 因其遗传信息简单且容易培养而被广泛用作各种生物化学、生物技术研究等的模式微生物。 本文首先通过 COBRA 工具包构建了大肠杆菌碳中心代谢系统的小规模网络,然后分别模拟敲除了其中的 25 种酶,并利用代谢通量平衡分析和代谢通量变异分析研究了基因敲除后的代谢网络。通过这些结果的综合分析,我们可以推断大肠杆菌碳中心代谢系统的不同的酶对整个体系的作用。 E. coli has clear genetic background and is easy to be cultured, so it has been widely used in a variety of biochemistry, biotechnology research as the mode microorganisms. In this work, a small E. coli metabolic system has been built in the cobra workspace. 25 enzymes in the metabolic system were separately knocked out and the metabolic system was investigated by FBA (flux balance analysis) and FVA (flux variability analysis). By combing the outcome of the analysis, we could infer the role played by the various enzymes in the E. coli metabolic system.
大肠杆菌(Escherichia coli, E. coli)革兰氏阴性短杆菌,大小0.5 × 1~3微米。周身鞭毛,能运动,无芽孢。能发酵多种糖类产酸、产气,是人和动物肠道中的正常栖居菌,婴儿出生后即随哺乳进入肠道,与人终身相伴,其代谢活动能抑制肠道内分解蛋白质的微生物生长,减少蛋白质分解产物对人体的危害,还能合成维生素B和K,以及有杀菌作用的大肠杆菌素[
大肠杆菌是研究微生物遗传的重要材料,如局限性转导就是1954年在大肠杆菌K12菌株中发现的。莱德伯格(Lederberg)采用两株大肠杆菌的营养缺陷型进行实验,奠定了研究细菌接合方法学上的基础,以及基因工程的研究。大肠杆菌作为外源基因表达的宿主,遗传背景清楚,技术操作简单,培养条件简单,大规模发酵经济,倍受遗传工程专家的重视。目前大肠杆菌是应用最广泛,最成功的表达体系,常作高效表达的首选体系。
系统生物学着重以系统的思路和方法定量研究研究生物过程[
同时,进一步,在另一个层面上,对代谢网络的研究也逐渐演化出来单独的静态分析和动态分析。静态分析主要包括Schilling等学者提出的极端途径(extreme pathway analysis),基元模式(elementary flux mode analysis),FBA流平衡分析(flux balance analysis)和jeong等学者提出的图论法(graph theory)[3-6]。而动态分析方法,依靠基于动力学信息的,它包括代谢控制分析(metabolic control analysis),控制论模拟方法(cybernetic modeling)。当前在基因组规模的代谢网络来说,许多反应的动力学参数和反应物浓度参数都是不清楚的,故静态的研究方法占据了目前代谢网络研究的主流位置[4,7-10]。
本篇论文在matlab cobra工具箱的基础上对大肠杆菌代谢系统进行FBA和FVA分析,将分析结果进行汇总后,为显示在敲除和不敲除各种酶的结果变化[
Matlab版本7.0以上cobra toolbox版本2.0以上libsbml版本4.0.1以上。
Sbml toolbox版本3.1.1以上gurobi 5.0以上。
大肠杆菌小型中心代谢模型“ecoli_core_model. Xml”,此模型中包含大肠杆菌中心代谢的全部信息[
在matlab命令窗口中执行指令: model = readCbModel(‘文件名’)其中文件名为模型保存目录[14,15]。
执行指令:
model = changeRxnBounds(model,'EX_glc(e)',-16.7,'l');其中'EX_glc(e)'代表了底物摄取速率,“l”是“lower”的缩写,表示模型的下限[
执行指令:
model = removeRxns(model,'rxnRemoveList'); 'rxnRemoveList'代表的是含有所敲除反应的反应信息的单元矩阵,酶的敲除中,rxnRemoveList可直接用所需敲除酶的英文缩写代替。如敲除烯醇化酶即ENO酶,直接在命令窗口中输入指令model = removeRxns(model,'ENO')即可,执行命令后表示该酶所直接参与的反应的电子信息从模型中去除,该反应在模型中不再起作用。
执行指令: changeCobraSolver('gurobi5','LP'); Fluxdistribution = optimizeCbModel(model); flux = geometricFBA(model)
其中第一步指令是修改cobra solver;第二步是执行标准FBA命令。
执行命令:
[minFlux,maxFlux] = fluxVariability(model,90);其中“90”代表一个百分比,为模型所能达到的目标函数相对于最优值得比例。FVA命令的结果有minFlux和maxFlux,分别代表模型中95个反应的通量分布情况。
作图从两个方向进行,分别为横向和纵向,横向作图在于显示模型中每个单独反应在不同酶被敲出后受到的影响,作图时横坐标代表被敲除的酶,为更好的显示结果,分别用数字代表各种酶,对应情况为:1) 没有敲除任何酶;2) EDD;3) GLUUPT;4) GdH;5) GND;6) TK2;7) PFK;8) HXI;9) TAL;10) TK1;11) SER;12) SHMT;13) GEV;14) ENO;15) PYK;16) PDH;17) AC-EX;18) GLTA;19) ICD;20) AKD;21) FUM;22) MDH;23) GS1;24) GS2;25) PPC;26) MEZ[
论文中只将反应ATP synthase,Biomass,glutamate transport的流量变化图附上(如图1~6),其他反应可以同理进行分析。
横向分析后,会发现在酶14(ENO),15(PYK),16(PDH),21(FUM),25(PPC)等被敲出后,反应受到的影响最大,本次论文只将酶ENO即烯醇化酶,酶
图1. 反应ATP synthase在敲除不同的酶后最优流量值的变化
图2. 反应ATP synthase在敲除不同的酶后流量范围的变化
图3. 反应Biomass在敲除不同的酶后最优流量值的变化
图4. 反应Biomass在敲除不同的酶后流量范围的变化
图5. 反应glutamate transport在敲除不同的酶后最优流量 值的变化
图6. 反应glutamate transport在敲除不同的酶后流量范围的变化
PDH即丙酮酸脱氢酶,酶FUM即延胡索酸酶,酶PPC即磷酸烯醇式丙酮酸羧化酶的纵向分析图附上(如图7~12),其它酶可以用同种方法进行分析。
综合横向分析和纵向分析,我们可以确定酶ENO即烯醇化酶,酶PDH即丙酮酸脱氢酶,酶FUM即延胡索酸酶,酶PPC即磷酸烯醇式丙酮酸羧化酶等对整个代谢系统的影响最大,以酶ENO为例,通过纵向分析来看,会发现,在ENO敲出后,在大肠杆菌代谢系统中有相当多的反应都受到了影响,而酶的催化作用具有专一性,即只催化一种或一类底物,因此我们可以推测出酶ENO所参与的反应应该是位于大肠杆菌代谢系统比较靠前的的位置,可以为代谢系统中之后的反应提供底物等。再返回汇总数据,可以发现,酶ENO,被敲出后反应代谢系统中反应就会少一个,说明在这些酶所直接参与的反应中,他们的作用是不可替代的,即大肠杆菌代谢系统中ENO酶无同工酶,由此我们也可以基本判定定酶ENO所催化的反应。再观察横向分析图与汇总的数据,会发现酶ENO对反应的影响可以分为两类:
1) ENO酶敲除后,反应被抑制,这又可以被分为
图7. 敲除酶ENO代谢系统中全部反应的FBA结果变化
图8. 敲除酶ENO代谢系统中全部反应FVA结果变化
图9. 敲除酶PDH代谢系统中全部反应的FBA结果变化
图10. 敲除酶PDH代谢系统中全部反应FVA结果变化
两类:a) ENO酶敲出后该反应被抑制,但是并未停止,如反应ATP synthease;b) ENO酶被敲出后该反应完
图11. 敲除酶FUM代谢系统中全部反应的FBA结果变化
图12. 敲除酶FUM代谢系统中全部反应FVA结果变化
全停止,如反应Biomass。
这说明代谢系统中反应ATP synthease和反应Biomass在酶ENO直接参与的反应后面进行。而且反应Biomass离ENO酶直接参与的反应可能更近一点,这样有利于我们对肠杆菌代谢网络结构的理解。
2) ENO酶被敲出后,反应被促进,如反应glutamate transport。考虑到大肠杆菌是兼性厌氧性微生物,因此猜测:酶ENO可抑制反应glutamate transport的进行或者是抑制反应glutamate transport底物的合成,也可能是酶ENO的产物对反应有抑制作用,即反应glutamate transport可能是大肠杆菌在另一种不同的环境才能进行的。
同理可以应用于大肠杆菌代谢系统中的其他酶和反应的分析,这将有利于我们对大肠杆菌代谢系统的认识。在此所用的大肠杆菌小型中心代谢系统模型,是比较简单的的代谢系统,这样的分析方案不仅适用于这样的小型代谢系统,也可用于对比较复杂的代谢网络分析。
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[