CN105574358B - 体外再现体内反应的方法 - Google Patents
体外再现体内反应的方法 Download PDFInfo
- Publication number
- CN105574358B CN105574358B CN201510926583.5A CN201510926583A CN105574358B CN 105574358 B CN105574358 B CN 105574358B CN 201510926583 A CN201510926583 A CN 201510926583A CN 105574358 B CN105574358 B CN 105574358B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- reaction
- metabolic
- cofactor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 231
- 238000000338 in vitro Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000002503 metabolic effect Effects 0.000 claims abstract description 228
- 239000000758 substrate Substances 0.000 claims abstract description 153
- 230000003834 intracellular effect Effects 0.000 claims abstract description 132
- 102000004190 Enzymes Human genes 0.000 claims abstract description 112
- 108090000790 Enzymes Proteins 0.000 claims abstract description 112
- 239000002207 metabolite Substances 0.000 claims abstract description 67
- 230000004060 metabolic process Effects 0.000 claims abstract description 52
- 230000004044 response Effects 0.000 claims abstract description 39
- 238000002474 experimental method Methods 0.000 claims abstract description 35
- 230000007246 mechanism Effects 0.000 claims abstract description 19
- 238000003198 gene knock in Methods 0.000 claims abstract description 8
- 239000000543 intermediate Substances 0.000 claims description 103
- 230000005764 inhibitory process Effects 0.000 claims description 83
- 238000001727 in vivo Methods 0.000 claims description 68
- 230000000694 effects Effects 0.000 claims description 28
- 238000006241 metabolic reaction Methods 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 22
- 238000006911 enzymatic reaction Methods 0.000 claims description 17
- 230000035484 reaction time Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000002401 inhibitory effect Effects 0.000 claims description 10
- 230000003197 catalytic effect Effects 0.000 claims description 9
- 238000003209 gene knockout Methods 0.000 claims description 9
- 230000000052 comparative effect Effects 0.000 claims description 8
- 230000035945 sensitivity Effects 0.000 claims description 8
- 230000001413 cellular effect Effects 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 6
- 230000002427 irreversible effect Effects 0.000 claims description 6
- 239000000376 reactant Substances 0.000 claims description 6
- 230000036632 reaction speed Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000013632 homeostatic process Effects 0.000 claims description 4
- 230000002829 reductive effect Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 claims description 3
- 210000002429 large intestine Anatomy 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000003204 osmotic effect Effects 0.000 claims description 3
- 230000036962 time dependent Effects 0.000 claims description 2
- 230000001629 suppression Effects 0.000 abstract 2
- 230000001419 dependent effect Effects 0.000 abstract 1
- 230000005055 memory storage Effects 0.000 abstract 1
- 108090000623 proteins and genes Proteins 0.000 description 10
- 241000282414 Homo sapiens Species 0.000 description 8
- 238000004088 simulation Methods 0.000 description 6
- 230000001351 cycling effect Effects 0.000 description 4
- 230000004907 flux Effects 0.000 description 4
- 238000010353 genetic engineering Methods 0.000 description 4
- 102000004169 proteins and genes Human genes 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000012269 metabolic engineering Methods 0.000 description 3
- 241000282412 Homo Species 0.000 description 2
- CXGXDTSWFKGIHD-CLCIBRLLSA-N [(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl] (2S,4aS,6aR,6aS,6bR,8aR,10S,12aS,14bR)-10-hydroxy-2,4a,6a,6b,9,9,12a-heptamethyl-13-oxo-3,4,5,6,6a,7,8,8a,10,11,12,14b-dodecahydro-1H-picene-2-carboxylate Chemical compound CC1(C)[C@@H](O)CC[C@@]2(C)[C@H]1CC[C@]1(C)[C@@H]2C(=O)C=C2[C@@H]3C[C@](C)(CC[C@]3(C)CC[C@@]12C)C(=O)O[C@@H]1O[C@H](CO)[C@@H](O)[C@H](O)[C@H]1O CXGXDTSWFKGIHD-CLCIBRLLSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000004883 computer application Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000000126 in silico method Methods 0.000 description 2
- 230000000977 initiatory effect Effects 0.000 description 2
- 230000037353 metabolic pathway Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 102000003960 Ligases Human genes 0.000 description 1
- 108090000364 Ligases Proteins 0.000 description 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001195 anabolic effect Effects 0.000 description 1
- 230000019522 cellular metabolic process Effects 0.000 description 1
- 230000036755 cellular response Effects 0.000 description 1
- 239000007806 chemical reaction intermediate Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002773 nucleotide Substances 0.000 description 1
- 125000003729 nucleotide group Chemical group 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000001243 protein synthesis Methods 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- -1 small molecule compounds Chemical class 0.000 description 1
- 150000003384 small molecules Chemical class 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000014616 translation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明属于生物代谢技术领域,提供了体外再现体内反应的方法,包括:打开计算机,调取数据库内存储的相关体内代谢的动力学参数;针对只含有辅因子代谢中间物CI、酶、底物和产物的最小代谢网络系统进行胞内和胞外的对比干实验;针对含有非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的无抑制代谢网络系统进行胞内和胞外的对比实验;在高、中、低三种抑制水平下,得出CI的浓度高低对NCI稳态水平的影响;通过调整该反应过程中酶的含量来模拟基因敲入、敲出实验得出酶量对NCI和反应速率的影响以及抑制机制、辅因子代谢中间物CI、酶的含量、非辅因子代谢中间物NCI对细胞内反应具有影响,进而得出胞外代谢反应系统。本发明弥补了行业空白。
Description
技术领域
本发明属于生物代谢技术领域,特别涉及一种体外再现体内反应的方法以及在生物等领域内的应用。
背景技术
生命体内的代谢过程是表现生命特征的唯一的方式,从基因表达,蛋白质的合成,到个体的生长,所有生命特征都是通过细胞代谢表现出来。因此了解代谢反应的内在机制对于揭示生命的全部过程是非常有必要的。
另一方面,体内(in vivo)反应和体外(in vitro)反应的区别是几百年来跨在人类面前的一条几乎难以逾越的鸿沟。科学家们可以清清楚楚地看见差不多所有的体内反应,组学数据已经大量积累,但是就是还没有人能够在体外重现一个体内反应系统。人类已经有能力(至少)部分组装基因组系统,但是还是无法复制在化学上了解近百年的代谢系统。目前,最前沿的代谢网络理论和方法顶多就是设法逼近细胞反应的某些性质,对“代谢启动/合成”和“代谢简化”这样的“合成生物学”的方向几乎没有实质性的促进。
目前,“合成生物学”着眼于基因组组装。即便如此,对于独立的细胞而言,基因复制所需的核苷酸和DNA合成酶只有通过合成代谢方能获得。目前,合成生物学研究者尽量绕开这一问题,代谢系统采用“借腹怀胎”的方法。
总的来说,胞内代谢反应的理解依然停留在静态理解代谢反应途径方面,而对于代谢反应网络动态特征的描述很少。
代谢网络包含几千个酶和大量的小分子化合物,这些小分子化合物的数量是在一千至一万之间。平均来说,每个小分子在胞内的含量都很低,为32uM。胞内必然有某种机制在保证反应速度的条件下使得胞内代谢物摩尔浓度很低。了解这种机制,将有助于人工复制代谢网络反应(可用于代谢工程改造),甚至构建非天然的代谢网络系统(可用于合成生物学)。此外代谢网络大都处于稳态,稳态对于细胞的存活具有重要的意义,在生物工程领域,稳态是实现生物产品高产的首要前提。
因此,为了解决上述问题,生物代谢技术领域急需一种体外再现体内反应的方法。
发明内容
本发明提供了一种体外再现体内反应的方法,技术方案如下:
体外再现体内反应的方法,包括如下实验步骤:
步骤一,打开计算机,调取数据库内存储的相关酶反应动力学参数,然后根据体内反应机理,在计算机中实现体内代谢网络反应;
步骤二,针对只含有辅因子代谢中间物CI、酶、底物和产物的最小代谢网络系统进行胞内和胞外的对比干实验,并且绘制出辅因子代谢中间物CI、底物以及产物其对应的浓度随时间的变化曲线图,得出辅因子代谢中间物CI在体内代谢反应中发挥着全局性的作用,而且辅因子代谢中间物CI在胞内的摩尔含量位于不太高的稳态,胞内代谢与胞外代谢相比的产物得率至少100-1000倍;
步骤三,针对含有非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的代谢网络系统在零抑制、低抑制、中抑制和高抑制四种抑制水平下进行胞内和胞外的对比干实验,并且绘制出非辅因子代谢中间物NCI、辅因子代谢中间物CI、底物以及产物其对应浓度随时间的变化曲线图;
步骤四,通过将步骤三条件下得出的多幅浓度随时间的变化曲线图进行比较,得出如下结论:(1)引入抑制后的胞内反应的非辅因子代谢中间物NCI都能维持稳态,稳态水平的高低由抑制强度的高低决定;(2)辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素;综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可;
步骤五,在高、中、低三种抑制水平下,通过调整CI在系统中的摩尔含量,得出CI的浓度高低对NCI稳态水平的影响,得出CI浓度越低,NCI稳态时的水平越低;
步骤六,针对步骤三的实验,通过调整该反应过程中酶的含量来模拟基因敲入、敲出实验,基因敲出实验能够通过设置酶量为零来模拟,得出酶量对非辅因子代谢中间物NCI的稳态摩尔含量和反应速率具有影响,而且能够定量的计算出代谢流对于不同酶的敏感程度;
步骤七,根据步骤二、步骤四至步骤六得出的结论,得出与胞内反应相一致的胞外代谢反应系统,具体如下:
其中,S1、S2代表底物,P1、P2代谢产物,B和Q代表辅因子代谢物CI,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7。
优选的,在上述体外再现体内反应的方法中,步骤一中根据体内反应机理,在计算机中实现体内代谢网络反应的具体步骤为:
已知在体内代谢网络中,双底物反应的比例为70-80%,第一、第二底物都过量,那么在体内代谢网络的反应中能够得到2个底物对酶的亲和性常数Km、1个酶的催化常数kcat以及最大反应速度Vmax;
Vmax=e1,0×kcat;
双底物的不可逆反应式为:S1+B→P1+Q,该反应中包含的基元反应为:
其中,k1,1,k1,-1,k1,2,k1,-2,k1,3,k1,-3k1,4,k1,-4均为酶反应的速率常数,S1和B表示反应的底物,P1和Q表示反应的产物,E1表示酶,E1B、E1BS1、E1Q表示酶复合物;
根据酶定量守恒定量,总酶量e1,0满足如下公式:
e1,0=e1+e1b+e1bs1+e1q
e1代表E1的摩尔量,e1b代表E1B的摩尔量,e1bs1代表E1BS1的摩尔量,e1q代表E1Q的摩尔量;
由于上述双底物反应式中和为不可逆反应,因此,k1,-3=k1,-4=0;
P1的生成速率的具体计算公式为:
其中,rp1代表P1的生成速率,s1代表底物S1的摩尔量,b代表底物B的摩尔量,由公式(1)可知,rp1除了由底物B和S1确定外,还由k1,1,k1,-1,k1,2,k1,-2,k1,3,k1,4和e1,0确定;
其中,Kmb表示底物B对酶的亲和性常数,并且Kmb越大,代表底物B对酶的亲和性越低;
其中,kcat表示酶的催化常数,并且kcat越大,代表没得催化效率越高;
其中,Kms1表示底物S1对酶的亲和性常数,并且Kms1越大,代表底物S1对酶的亲和性越低;
首先,假设s1>>b,同时将公式(2)和(3)代入公式公式(1)中,得出:
然后,将公式(3)除以公式(2),得出k1,1,具体如下:
由公式(6)可知,k1,1为固定值;
进一步地,假设b>>s1,同时将公式(2)和(3)代入公式公式(1)中,得出:
进一步地,将公式(2)、(3)和(4)代入公式(1)中,得出:
由公式(8)可知,rp1除了由底物B和S1确定外,还由总酶量e1,0、酶的催化常数kcat、底物对酶的亲和性常数Kms1、Kmb以及酶反应的速率常数k1,-1确定;因此,在给定底物浓度s1、b以及产物的生成速率rp1时,能够得出e1,0,kcat,Kms1,Kmb,andk1,-1的关系式,进而从计算机的数据库内提取出满足e1,0,kcat,Kms1,Kmb,andk1,-1的关系式酶反应动力学参数;
由公式(3)可知,在kcat求得后,通过任意调整k1,3和k1,4的值,使其满足rp1的定值即可;
将公式(4)除以公式(3),得出:
由公式(9)可知,在满足Kms1和kcat一定的条件下,通过任意调整k1,-2和k1,3的值,使其满足rp1的定值即可;
因此,由公式(1)至(9)可知,对于双底物反应,在得到米氏动力学常数kcat、 Kms1、Kmb,能够求得k1,1、k1,-1、k1,-2、k1,-2、k1,3、k1,-3、k1,4和k1,-4;同理,对于单底物或者三底物反应,当给定的动力学参数的个数少于需要的速率常数的个数时,能够按照公式(1)至(9)的方法,找到数据库参数和速率常数之间的联系,然后求得各个速率常数k,再完成步骤二的步骤。
优选的,在上述体外再现体内反应的方法中,步骤二的具体步骤为:
(A)胞内代谢网络与基元反应式:
(B)胞外代谢网络与基元反应式:
其中,(A)胞内代谢网络及与它相关的基元反应中,B和Q为一组辅因子代谢物CI,他们在网络中不断消耗和产生,且不会耗尽也不会被排出网络,因此B+Q是一个含量很低的常数;(B)胞外代谢网络及与它相关的基元反应中,(1)中的B和S1为底物,Q和P1为相应的产物;(2)中的Q和S2为底物,B和P2为产物;并且(1)中产生的Q不能被(2)作为底物使用,(2)中产生的B也不能被(1)做底物使用;
为了让(A)胞内代谢网络和(B)胞外代谢网络具有可比性,将S1-P1和S2-P2的反应线路设计成完全相同的;即:S1=S2,P1=P2,E1=E2,B=Q,k1,j=k2,j,k1,-j=k2,-j,j=1,2,3,4;因此对于胞内反应来说,根据步骤一中获得的动力学参数进行筛选,得出最小代谢网络的相关动力学参数表格;
在阐述胞内代谢网络相比于胞外代谢网络的优势之前,首先将引入四个定义:反应速率、产物得率、相对产物得率和相对反应速率;其中,反应速率分为产物反应速率rpi和底物反应速率rsi,基于CI的产物得率为qCI(s-1),胞内反应系统的基于CI的产物得率为胞外反应系统的基于CI的产物得率为基于CI的相对产物得率为Rq/q,Rq/q是客观衡量胞内反应相对于胞外反应的效率的标准;CI的相对反应速率为Rr/r,Rr/r是客观衡量胞内反应相对与胞外反应的速率的标准,具体计算公式如下:
其中,in-vivo表示胞内反应,in-vitro表示胞外反应,Δt表示反应时间,Δpi表示从产物积累的角度衡量反应的进程,Δsi表示从底物消耗的角度衡量反应的进程,bss是B在稳态时的值,b0表示B与S1同时耗尽;
根据最小代谢网络的相关动力学参数表格中的数值以及反应物、产物、辅因子代谢中间物CI与时间的变化曲线进行比较,得出如下结论:(1)胞内和胞外代谢网络的反应时间Δt、底物的消耗量ΔS和产物的生成量ΔP完全相同,因此胞内反应相比于胞外反应的基于CI的产物得率根据上述反应结果,表明反应速率中的快反应和慢反应,虽然会影响反应时间,但不会影响基于CI的相对产物得率Rq/q,而底物的初始值会对相对反应速率有影响;(2)胞内和胞外代谢网络中B和Q随着S1和S2的消耗而不断的消耗,因此在这样的胞外系统中用B和Q的初始值来计算基于CI的产物得率更加准确合理;尽管这个最小代谢反应不是真实的,但是它展示了辅因子代谢中间物CI一定是发挥着全局和局部的作用,而且他们在胞内位于低浓度的稳态,胞内代谢相比与胞外代谢的产物得率至少为100-1000倍。
优选的,在上述体外再现体内反应的系统中,步骤三的具体步骤为:
(C)含有NCI、无抑制胞内代谢网络反应过程为:
(D)含有NCI、无抑制胞外代谢网络反应过程为:
(C)和(D)中含有NCI、无抑制胞内和胞外代谢网络的基元反应式:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此在胞外代谢网络中B和Q不具有辅因子代谢物CI的全局性作用,全局性是指(1)或(2)产生的B或Q都作为底物再次参加(2)或(1)的反应,即循环反应;
(E)含有NCI和产物抑制的胞内代谢网络的反应过程为:
(F)含有NCI和产物抑制的胞外代谢网络的反应过程为:
(E)和(F)中含有NCI和产物抑制的胞内和胞外代谢网络的基元反应为:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此胞外系统中B和Q,不具有辅因子代谢物CI的全局性作用,全局性是指(1)或(2)产生的B或Q都作为底物再次参加(2)或(1)的反应,即循环反应;
对于胞内反应来说,根据步骤一中获得的动力学参数进行筛选,得出含有NCI的相关动力学参数表格;
将步骤二的最小代谢网络与步骤三中的代谢网络相比较,得知:在最小代谢网络中,由于没有NCI,总是有ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro,然而在目前复杂的模型中不能使ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro同时成立,否则此时胞内反应中B和Q的摩尔含量将高达30mM左右;在真实的细胞内,渗透压的原因以及细胞摩尔含量的可调性决定了CI的摩尔含量不会太高,大肠内CI的平均含量是1.5mM;因此我们采取了降低反应速度或者降低CI的摩尔含量的方法,从而延长反应时间的方法来评价Rq/q;
由含有NCI的相关动力学参数表格得出,产物抑制一方面减小了产物得率,另一方面也减小了反应的速度;然而,当减小辅因子代谢中间物CI在代谢网络中的含量时,基于CI的相对产物得率却明显上升,虽然这其中伴随值反应速率的略微减小,但是这一损失对细胞来说是划算的;较低的CI摩尔含量提高了代谢的效率,增强了代谢物摩尔含量对代谢流的敏感性,降低了反应速率从而保证其他每一个代谢中间物的摩尔含量都被控制在较低的水平;总的来说这一结果与我们在最小代谢网络中的结果完全一致,从而进一步巩固了CI对与细胞维持高效代谢活动的全局性作用;
进一步地,根据参数绘制在胞内以及胞外代谢网络中非辅因子代谢中间物NCI、反应物和产物的对应物质的量随时间的变化曲线图;
优选的,在上述体外再现体内反应的方法中,步骤四的具体步骤为:在没有抑制的条件下,非辅因子代谢中间物NCI有的维持稳态,有的没有维持稳态;但是当引入抑制后,所有胞内反应的非辅因子代谢中间物NCI都能维持稳态;稳态水平的高低则由抑制强度的高低来决定;
比较胞内反应和胞外反应,虽然胞外反应与对应的胞内反应有相同的抑制机制、抑制水平,但是并不是所有的NCI都有稳态,考虑此时胞内网络和胞外反应唯一的区别是前者包含CI而后者虽然包含CI,但是不具有CI的性质,相当于胞外网络不具有CI,因此唯一的区别在于胞内含有CI,因此得出结论:辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素,而且CI摩尔含量越低,NCI稳态时的摩尔含量水平越低,综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可。
优选的,在上述体外再现体内反应的系统中,步骤五的具体步骤为:在高、中、低三种抑制水平下以CI的浓度大小为自变量,得出CI浓度的改变对胞内系统中其他NCI稳态水平的影响;得出如下结论:NCI总体对CI的改变是很敏感的,当CI浓度很低时,NCI对CI的敏感度最强,因此能够通过降低CI的浓度,使得NCI水平降低,从而提高产物得率。
优选的,在上述体外再现体内反应的方法中,步骤六的具体步骤为:
绘制单底物和双底物中酶反应的酶含量对胞内网络中非辅因子代谢中间物NCI的摩尔含量随着时间的变化曲线图,得出在胞内网路中,酶量的影响是全局性的;增加酶量会使得前面的非辅因子代谢中间物NCI的稳态浓度下降,后面的非辅因子代谢中间物NCI的稳态浓度上升;同时由于B+Q是常数,降低B会使Q升高,从而将酶量的影响由局部扩大到全局;
酶量对反应速率会产生影响,酶量越大,反应速率越大,但是反应速率对不同酶的敏感程度不一样,因此,酶量不管是反应速率还是对非辅因子代谢中间物NCI的稳态浓度的影响都是全局性的。
优选的,在上述体外再现体内反应的方法中,还包括:步骤八,根据所述步骤七中得出的胞外代谢反应系统,将数据库内的相关体内代谢的动力学参数代入胞外代谢反应系统内,在计算机上进行大量的干反应,最终得出一组非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量,使其性能最好,然后在实验室进行湿实验,对胞外代谢反应系统进行微调,得出最优的非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量值,即为最优的胞外代谢反应系统。
本发明的有益效果:
1、本发明是一项既有着深刻的理论影响又有着新颖的实用意义的创造性工作,提出了一种在合成生物学和代谢网络方面的全新的思维模式,从而独创性的构建了一个可以复制并放大的的体外再现体内反应的系统而且本发明必将在基因工程,代谢工程和合成生物学等领域中发挥着重要作用,帮助人类解决日益面临的能源和环境问题,弥补了行业内的空白,具有重要的理论和实践意义。
2、本发明通过数值模拟显示了(a)CI类代谢物全内循环式代谢流的全局性效应,(b)NCI类代谢物等各种反馈抑制的局域性效应,(c)低CI水平更有利于提高胞内代谢反应的产物的得率;这些结果同时验证了本发明的正确性。
3、本发明还通过数值模拟显示了酶是如何全局性地影响代谢系统中所有代谢物和代谢通量的,为计算机模拟基因工程操作提高一个非常有效的方法。
4、应用本发明可以构建出完整的细胞代谢系统,通过调试或更改系统,具体为小规模或大规模的基因变更,来寻求最优目标条件或目标产品,为使用活细胞实验室工作的开展节约了大量的时间,降低了劳动量,节约了成本,使实验的安排更加科学。
5、本发明适用于计算机应用模拟基因或基因组变更以及合成生物学研究等领域,应用领域广泛,提高湿实验的效率,更高效的得到改良的菌种,人类需要的各种生物制品,以解决人类面临的各种能源、疾病和环境等问题;在合成in-silico细胞基础上,与基因组,转录组等网络相匹配,最终达到合成人类所需的某些特定的生命体或类生命体,具有更加广泛的适用性。
附图说明
下面结合附图和具体实施方式来详细说明本发明:
图1为体外再现体内反应的方法的流程图。
图2-5为在最小代谢网络系统中,辅因子代谢中间物CI、底物以及产物所对应物质的量随时间的变化曲线图。
图6为在含有辅因子代谢物的胞内代谢网络系统中,在高,中,低抑制水平下,CI的浓度对相对产物得率Rq/q和相对反应速率Rr/r的影响曲线图。
图7为在含有辅因子代谢物的代谢网络系统中,在零抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图。
图8为在含有辅因子代谢物的代谢网络系统中,在低抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图。
图9为在含有辅因子代谢物的代谢网络系统中,在中等抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图。
图10为在含有辅因子代谢物的代谢网络系统中,在高抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图。
图11为在含有辅因子代谢物的胞内代谢网络系统中,在高,中和低抑制水平下,胞内系统中辅因子代谢物B、Q的浓度对非辅因子代谢物Y1、Y2、Y3、Y4、Y5、Y6稳态浓度的影响曲线图。
图12为在含有辅因子代谢物的胞内代谢网络系统中,模拟基因敲入实验中逐渐增大E2和E3的含量对NCI、B、Q稳态水平的影响曲线图。
图13为在含有辅因子代谢物的胞内代谢网络系统中,模拟基因敲入实验,逐渐增大E2和E3的含量对流过2条代谢途径的代谢流的影响曲线图。
具体实施方式
为了使本发明技术实现的措施、创作特征、达成目的与功效易于明白了解,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
图1为体外再现体内反应的方法的流程图。
如图1所示,本发明的一种体外再现体内反应的方法,包括如下步骤:
步骤一,打开计算机,调取数据库内存储的相关酶反应动力学参数,然后根据体内反应机理,在计算机中实现体内代谢网络反应;
步骤二,针对只含有辅因子代谢中间物CI、酶、底物和产物的最小代谢网络系统进行胞内和胞外的对比干实验,并且绘制出辅因子代谢中间物CI、底物以及产物其对应的浓度随时间的变化曲线图,得出辅因子代谢中间物CI在体内代谢反应中发挥着全局性的作用,而且辅因子代谢中间物CI在胞内的摩尔含量位于不太高的稳态,胞内代谢与胞外代谢相比的产物得率约为100-1000倍;
步骤三,针对含有非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的代谢网络系统在零抑制、低抑制、中抑制和高抑制四种抑制水平下进行胞内和胞外的对比干实验,并且绘制出非辅因子代谢中间物NCI、辅因子代谢中间物CI、底物以及产物所对应浓度随时间的变化曲线图;
步骤四,通过将步骤三条件下得出的多幅浓度随时间的变化曲线图进行比较,得出如下结论:(1)引入抑制后的胞内反应的非辅因子代谢中间物NCI都能维持稳态,稳态水平的高低由抑制强度的高低决定;(2)辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素;综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可;
步骤五,在高、中、低三种抑制水平下,通过调整CI在系统中的摩尔含量,得出CI的浓度高低对NCI稳态水平的影响,得出CI浓度越低,NCI稳态时的水平越低;
步骤六,针对步骤三的实验,通过调整该反应过程中酶的含量来模拟基因敲入、敲出实验,基因敲出实验能够通过设置酶量为零来模拟,得出酶量对非辅因子代谢中间物NCI的稳态摩尔含量和反应速率具有影响,而且能够定量的计算出代谢流对于不同酶的敏感程度;
步骤七,根据步骤二、步骤四至步骤六得出的结论,得出与胞内反应相一致的胞外代谢反应系统,具体如下:
其中,S1、S2代表底物,P1、P2代谢产物,B和Q代表辅因子代谢物CI,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7。
实施例2:
图1是体外再现体内反应的方法的流程图。
如图1所示,体外再现体内反应的方法,包括如下步骤:
步骤一,打开计算机,调取数据库内存储的相关酶反应动力学参数,然后根据体内反应机理,在计算机中实现体内代谢网络反应;
已知在体内代谢网络中,双底物反应的比例为70-80%,第一、第二底物都过量,那么在体内代谢网络的反应中能够得到2个底物对酶的亲和性常数Km、1个酶的催化常数kcat以及最大反应速度Vmax;
Vmax=e1,0×kcat;
双底物的不可逆反应式为:S1+B→P1+Q,该反应中包含的基元反应为:
其中,k1,1,k1,-1,k1,2,k1,-2,k1,3,k1,-3k1,4,k1,-4均为酶反应的速率常数,S1和B表示反应的底物,P1和Q表示反应的产物,E1表示酶,E1B、E1BS1、E1Q表示酶复合物;
根据酶定量守恒定量,总酶量e1,0满足如下公式:
e1,0=e1+e1b+e1bs1+e1q
e1代表E1的摩尔量,e1b代表E1B的摩尔量,e1bs1代表E1BS1的摩尔量,e1q代表E1Q的摩尔量;
由于上述双底物反应式中和为不可逆反应,因此,k1,-3=k1,-4=0;
P1的生成速率的具体计算公式为:
其中,rp1代表P1的生成速率,s1代表底物S1的摩尔量,b代表底物B的摩尔量,由公式(1)可知,rp1除了由底物B和S1确定外,还由k1,,k1,-1,k1,2,k1,-2,k1,3,k1,4和e1,0确定;
其中,Kmb表示底物B对酶的亲和性常数,并且Kmb越大,代表底物B对酶的亲和性越低;
其中,kcat表示酶的催化常数,并且kcat越大,代表没得催化效率越高;
其中,Kms1表示底物S1对酶的亲和性常数,并且Kms1越大,代表底物S1对酶的亲和性越低;
首先,假设s1>>b,同时将公式(2)和(3)代入公式公式(1)中,得出:
然后,将公式(3)除以公式(2),得出k1,1,具体如下:
由公式(6)可知,k1,1为固定值;
进一步地,假设b>>s1,同时将公式(2)和(3)代入公式公式(1)中,得出:
进一步地,将公式(2)、(3)和(4)代入公式(1)中,得出:
由公式(8)可知,rp1除了由底物B和S1确定外,还由总酶量e1,0、酶的催化常数kcat、底物对酶的亲和性常数Kms1、Kmb以及酶反应的速率常数k1,-1确定;因此,在给定底物浓度s1、b以及产物的生成速率rp1时,能够得出e1,0,kcat,Kms1,Kmb,andk1,-1的关系式,进而从计算机的数据库内提取出满足e1,0,kcat,Kms1,Kmb,andk1,-1的关系式酶反应动力学参数;
由公式(3)可知,在kcat求得后,通过任意调整k1,3和k1,4的值,使其满足rp1的定值即可;
将公式(4)除以公式(3),得出:
由公式(9)可知,在满足Kms1和kcat一定的条件下,通过任意调整k1,-2和k1,3的值,使其满足rp1的定值即可;
因此,由公式(1)至(9)可知,对于双底物反应,在得到米氏动力学常数kcat、Kms1、Kmb,能够求得k1,1、k1,-1、k1,2、k1,-2、k1,3、k1,-3、k1,4和k1,-4;同理,对于单底物或者三底物反应,当给定的动力学参数的个数少于需要的速率常数的个数时,能够按照公式(1)至(9)的方法,找到数据库参数和速率常数之间的联系,然后求得各个速率常数k,再完成步骤二的步骤;
步骤二,针对只含有辅因子代谢中间物CI、酶、底物和产物的最小代谢网络系统进行胞内和胞外的对比干实验,并且绘制出辅因子代谢中间物CI、底物以及产物其对应的浓度随时间的变化曲线图,图2-5为在最小代谢网络系统中,辅因子代谢中间物CI、底物以及产物所对应物质的量随时间的变化曲线图,如图2-5所示,得出辅因子代谢中间物CI在体内代谢反应中发挥着全局性的作用,而且辅因子代谢中间物CI在胞内的摩尔含量位于不太高的稳态,胞内代谢与胞外代谢相比的产物得率至少100-1000倍;具体如下:
(A)胞内代谢网络与基元反应式:
(B)胞外代谢网络与基元反应式:
其中,(A)胞内代谢网络及与它相关的基元反应中,B和Q为一组辅因子代谢物CI,他们在网络中不断消耗和产生,且不会耗尽也不会被排出网络,因此B+Q是一个含量很低的常数;(B)胞外代谢网络及与它相关的基元反应中,(1)中的B和S1为底物,Q和P1为相应的产物;(2)中的Q和S2为底物,B和P2为产物;并且(1)中产生的Q不能被(2)作为底物使用,(2)中产生的B也不能被(1)做底物使用;
为了让(A)胞内代谢网络和(B)胞外代谢网络具有可比性,将S1-P1和S2-P2的反应线路设计成完全相同的;即:S1=S2,P1=P2,E1=E2,B=Q,k1,j=k2,j,k1,-j=k2,-j,j=1,2,3,4;因此对于胞内反应来说,根据步骤一中获得的动力学参数进行筛选,得出如表一中所呈现出的最小代谢网络的相关动力学参数,具体数据见下表:
表一 最小代谢模型反应的相关动力学参数和初始条件
动力学参数/参数条件 | 慢反应 | 快反应 |
k1,1=k2,1(mMs)-1 | 105 | 105 |
k1,-1=k2,-1(s-1) | 6.67×103 | 6.67×103 |
k1,2=k2,2(mMs)-1 | 1.33×105 | 1.33×105 |
k1,-2=k2,-2(s-1) | 4.56×104 | 4.56×104 |
k1,3=k2,3(s-1) | 133 | 1330 |
k1,-3=k2,-3 | 0 | 0 |
k1,4=k2,4(s-1) | 100 | 1000 |
k1,-4=k2,-4 | 0 | 0 |
S1=S2(mM) | 100或500 | 100或500 |
P1=P2 | 0 | 0 |
E1=E2(mM) | 0.013 | 0.013 |
酶反应中间体 | 0 | 0 |
在阐述胞内代谢网络相比于胞外代谢网络的优势之前,首先将引入四个定义:反应速率、基于CI的产物得率、基于CI的相对产物得率和相对反应速率;其中,反应速率分为产物反应速率rpi和底物反应速率基于CI的产物得率为qCI(s-1),胞内反应系统的基于CI的产物得率为胞外反应系统的基于CI的产物得率为基于CI的相对产物得率为Rq/q,Rq/q是客观衡量胞内反应相对于胞外反应的效率的标准;CI的相
对反应速率为Rr/r,Rr/r是客观衡量胞内反应相对与胞外反应的速率的标准,具体计算公式如下:
其中,in-vivo表示胞内反应,in-vitro表示胞外反应,Δt表示反应时间,Δpi表示从产物积累的角度衡量反应的进程,Δsi表示从底物消耗的角度衡量反应的进程,bss为B在胞内反应系统中稳态时的值,b0为胞外反应中B的初始值,一般情况下应该保证加入的b0与底物S1同时耗尽;
我们在多种不同的外界条件比较了相对产物得率包括:低底物摩尔含量S=100mM和高底物摩尔含量S=500mM,慢反应k1,3=133(s-1)和k1,4=100(s-1),快反应k2,3=1330(s-1)和k2,4=1000(s-1),在最小代谢网络中,和几乎相同,且Δtin vivo和Δtin vitro几乎相同,当S0=100时,Rq/q=133,S0=500时,Rq/q=666,具体数值如表二;
表二 最小胞内代谢网络和相应的胞外代谢网络计算的初始条件和结果
根据表二中的数值以及图1中的反应物、产物、辅因子代谢中间物CI与时间的变化曲线进行比较,得出如下结论:(1)A(b),B(b),C(b),D(b)表明胞内和胞外代谢网络的反应时间Δt、底物的消耗量ΔS和产物的生成量ΔP完全相同,因此胞内反应相比于胞外反应的基于CI的产物得率根据上述反应结果,表明反应速率中的快反应和慢反应,虽然会影响反应时间,但不会影响基于CI的相对产物得率Rq/q,而底物的初始值会对相对反应速率有影响;(2)根据A(a),B(a),C(a),D(a)曲线图,表明胞内和胞外代谢网络中B和Q随着S1和S2的消耗而不断的消耗,因此在这样的胞外系统中用B和Q的初始值来计算基于CI的产物得率更加准确合理;尽管这个最小代谢反应不是真实的,但是它展示了辅因子代谢中间物CI一定是发挥着全局的作用,而且他们在胞内位于低浓度的稳态,如表二和图2-5所示,当S0=100mM时,Rq/q=133,快反应反应时间为13.6s,慢反应为136s,S0=500mM时,Rq/q=666,快反应反应时间为68s,慢反应为680s;而且根据实际反应时间来看,底物浓度几乎不能在降低,但如果底物浓度S0继续升高至1000mM,将会使得胞外反应所需要的b0同时升高至1000mM升高,这将使得基于CI的产物得率反应Rq/q由666升高到1332。因此保守估计胞内代谢相比与胞外代谢的产物得率至少为100-1000倍;
步骤三,针对含有非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的代谢网络系统在零抑制、低抑制、中抑制和高抑制四种抑制水平下进行胞内和胞外的对比干实验,并且绘制出非辅因子代谢中间物NCI、辅因子代谢中间物CI、底物以及产物其对应浓度随时间的变化曲线图,如图6-8所示,图6为在含有辅因子代谢物的胞内代谢网络系统中,在高,中,低抑制水平下,CI的浓度对相对产物得率Rq/q和相对反应速率Rr/r的影响曲线图;图7为在含有辅因子代谢物的代谢网络系统中,在零抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图;图8为在含有辅因子代谢物的代谢网络系统中,在低抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图;
(C)含有NCI、无抑制胞内代谢网络反应过程为:
(D)含有NCI、无抑制胞外代谢网络反应过程为:
(C)和(D)中含有NCI、无抑制胞内和胞外反系统的基元反应为:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此胞外系统中B和Q,不具有辅因子代谢物CI的全局性作用,全局性是指(1)或(2)产生的B或Q都作为底物再次参加(2)或(1)的反应,即循环反应;
(E)含有NCI和产物抑制的胞内代谢网络的反应过程为:
(F)含有NCI和产物抑制的胞外代谢网络的反应过程为:
(E)和(F)中含有NCI和产物抑制的胞内和胞外代谢网络的基元反应为:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此胞外系统中B和Q,不具有辅因子代谢物CI的全局性作用,全局性是指(1)或(2)产生的B或Q都作为底物再次参加(2)或(1)的反应,即循环反应;
对于胞内和胞外网络反应来说,根据步骤一中获得的动力学参数进行筛选,得出如表三中所呈现出含有NCI的相关动力学参数及初始条件,具体数据见下表三:
表三 含有NCI的相关动力学参数和初始条件
将步骤二中的最小代谢网络与步骤三中的代谢网络相比较,得知:在最小代谢网络中,由于没有NCI,总是有ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro,然而在目前复杂的模型中不能使ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro同时成立,否则此时胞内反应中B和Q的摩尔含量将高达30mM左右;在真实的细胞内,渗透压的原因以及细胞摩尔含量的可调性决定了CI的摩尔含量不会太高,大肠内CI的平均含量是1.5mM;因此我们采取了降低反应速度或者降低CI的摩尔含量的方法,从而延长反应时间的方法来评价Rq/q;
当s=500mM时,使用表三中的动力学参数和初始条件,在高、中、低三个不同抑制水平的条件下绘制基于CI的相对产物得率和相对速率变化曲线,如图2-5所示;
根据图2-5可知,低抑制水平时,当B=1mM、0.1mM和0.05mM时,Rq/q分别为:7.5×102,6.1×103和1.1×104,而此时Rr/r分别为0.97,0.8,0.7;高抑制水平时,如果B=1mM、0.1mM和0.05mM时,Rq/q分别为:3.5×102,1.4×103和2.1×103,Rr/r分别为0.35,0.14,0.14;因此,得出产物抑制一方面减小了产物得率,另一方面也减小了反应的速度;然而,当减小辅因子代谢中间物CI在代谢网络中的含量时,基于CI的相对产物得率却明显上升,虽然这其中伴随着应速率的略微减小,但是这一损失对细胞来说是划算的;较低的CI摩尔含量提高了代谢的效率,增强了代谢流对代谢物摩尔含量对的敏感性,降低了反应速率从而保证其他每一个代谢中间物的摩尔含量都被控制在较低的水平;总的来说这一结果与我们在最小代谢网络中的结果完全一致,从而进一步巩固了CI对与细胞维持高效代谢活动的全局性作用;
进一步地,根据表三中的参数绘制在胞内以及胞外代谢网络中非辅因子代谢中间物NCI、反应物和产物的对应物质的量随时间的变化曲线图,如图6-10,图9为在含有辅因子代谢物的代谢网络系统中,在中等抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图,图10为在含有辅因子代谢物的代谢网络系统中,在高抑制水平,胞内系统和胞外系统中Y1、Y2、Y3、Y4、Y5、Y6、B、Q、底物S1、S2和产物P1、P2的反应进程图;
步骤四,通过将步骤三条件下得出的多幅浓度随时间的变化曲线图进行比较,得出如下结论:(1)引入抑制后的胞内反应的非辅因子代谢中间物NCI都能维持稳态,稳态水平的高低由抑制强度的高低决定;(2)辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素;综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可,具体解释如下;
根据图6和表三可知,在没有抑制的条件下,非辅因子代谢中间物NCI除了Y1和Y4没有维持稳态外,Y2、Y3、Y5和Y6都能维持稳态;但是当引入抑制后,根据图8-10所示,所有胞内反应的非辅因子代谢中间物NCI都能维持的稳态;至于稳态水平的高低则由抑制强度的高低来决定;
比较图7-10中的胞内反应和胞外反应,虽然胞外反应与对应的胞内反应有相同的抑制机制、抑制水平,但是并不是所有的NCI都有稳态,图7-10胞外反应系统中Y1和Y4仍然没有稳态,考虑此时胞内网络和胞外反应唯一的区别是前者包含CI而后者虽然包含CI,但是不具有CI的性质,相当于胞外网络不具有CI,因此唯一的区别在于胞内含有CI,因此得出结论:辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素,而且CI摩尔含量越低,NCI稳态时的摩尔含量水平越低,综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可;
步骤五,在高、中、低三种抑制水平下,通过调整CI在系统中的摩尔含量,得出CI的浓度高低对NCI稳态水平的影响,得出CI浓度越低,NCI稳态时的水平越低;
在高、中、低三种抑制水平下以CI的浓度大小为自变量,得出CI浓度的改变对胞内系统中其他NCI稳态水平的影响,得出如下结论:NCI总体对CI的改变是很敏感的,当CI浓度很低时,NCI对CI的敏感度最强,因此能够通过降低CI的浓度,使得NCI水平降低,从而提高产物得率;
步骤六,针对步骤三的实验,通过调整该反应过程中酶的含量来模拟基因敲入、敲出实验,基因敲出实验能够通过设置酶量为零来模拟,得出酶量对非辅因子代谢中间物NCI的稳态摩尔含量和反应速率具有影响,而且能够定量的计算出代谢流对于不同酶的敏感程度,具体解释如下:
图11为在含有辅因子代谢物的胞内代谢网络系统中,在高,中和低抑制水平下,胞内系统中辅因子代谢物B、Q的浓度对非辅因子代谢物Y1、Y2、Y3、Y4、Y5、Y6稳态浓度的影响曲线图,图12为在含有辅因子代谢物的胞内代谢网络系统中,模拟基因敲入实验中逐渐增大E2和E3的含量对NCI、B、Q稳态水平的影响曲线图,图13为在含有辅因子代谢物的胞内代谢网络系统中,模拟基因敲入实验,逐渐增大E2和E3的含量对流过2条代谢途径的代谢流的影响曲线图;
绘制Y1单底物和Y2双底物酶反应中酶含量对胞内网络中非辅因子代谢中间物NCI的摩尔含量随着时间的变化曲线图,如图11所示;根据图11可知,在胞内网络中,酶量的影响是全局性的,而人们之前基本认为酶量的作用是局部性的;增加E2的量使得Y2、Y3、Y4、Y5、Y6和Q的稳态浓度上升,而Y1和B的稳态浓度下降;根据图12所示,增加E3的量使得Y3、Y4、Y5、Y6和Q的稳态浓度上升,而Y1、Y2和B的稳态浓度下降;因此得出增加酶量会使得前面的非辅因子代谢中间物NCI的稳态浓度下降,后面的非辅因子代谢中间物NCI的稳态浓度上升;同时由于B+Q是常数,降低B会使Q升高,从而将酶量的影响由局部扩大到全局;根据图13可知,酶量对反应速率会产生影响,酶量越大,反应速率越大,但是反应速率对不同酶的敏感程度不一样,反应速率对E2比E3更敏感;计算代谢流控制系数flux controlcoefficients,简称FCC,可知:当e2,0=0.006和0.012mM时,FCC分别为0.295和0.151,当e3,0=0.06和0.12mM,FCC分别为0.106和0.073;因此,酶量不管是反应速率还是对非辅因子代谢中间物NCI的稳态浓度的影响都是全局性的;而且证明通过计算机模拟基因敲入或敲出是可行的,这一做法必将极大的提高基因工程湿实验的效率;
步骤七,根据步骤二、步骤四至步骤六得出的结论,得出与胞内反应相一致的胞外代谢反应系统,具体如下:
其中,S1、S2代表底物,P1、P2代谢产物,B和Q代表辅因子代谢物CI,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7;
步骤八,根据步骤七中得出的胞外代谢反应系统,将数据库内的相关体内代谢的动力学参数代入胞外代谢反应系统内,在计算机上进行大量的干反应,最终得出一组非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量,使其性能最好,然后在实验室进行湿实验,对胞外代谢反应系统进行微调,得出最优的非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量值,即为最优的胞外代谢反应系统。
本发明是一项既有着深刻的理论影响又有着新颖的实用意义的创造性工作,提出了一种在合成生物学和代谢网络方面的全新的思维模式,从而独创性的构建了一个可以复制并放大的的体外再现体内反应的系统而且本发明必将在基因工程,代谢工程和合成生物学等领域中发挥着重要作用,帮助人类解决日益面临的能源和环境问题,弥补了行业内的空白,具有重要的理论和实践意义。
本发明通过数值模拟显示了(a):CI类代谢物全内循环式代谢流的全局性效应,(b):NCI类代谢物等各种反馈抑制的局域性效应,(c):低CI水平更有利于提高胞内代谢反应的产物的得率;这些结果同时验证了本发明的正确性。
本发明还通过数值模拟显示了酶是如何全局性地影响代谢系统中所有代谢物和代谢通量的,为计算机模拟基因工程操作提高一个非常有效的方法。
应用本发明可以构建出完整的细胞代谢系统,通过调试或更改系统,具体为小规模或大规模的基因变更,来寻求最优目标条件或目标产品,为使用活细胞实验室工作的开展节约了大量的时间,降低了劳动量,节约了成本,使实验的安排更加科学。
本发明适用于计算机应用模拟基因或基因组变更以及合成生物学研究等领域,应用领域广泛,提高湿实验的效率,更高效的得到改良的菌种,人类需要的各种生物制品,以解决人类面临的各种能源、疾病和环境等问题;在合成in-silico细胞基础上,与基因组,转录组等网络相匹配,最终达到合成人类所需的某些特定的生命体或类生命体,具有更加广泛的适用性。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点.本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内.本发明要求保护范围由所附的权利要求书及其等同物界定。
Claims (7)
1.体外再现体内反应的方法,其特征在于,包括如下步骤:
步骤一,打开计算机,调取数据库内存储的相关酶反应动力学参数,然后根据体内反应机理,在计算机中实现体内代谢网络反应;
步骤二,针对只含有辅因子代谢中间物CI、酶、底物和产物的最小代谢网络系统进行胞内和胞外的对比干实验,并且绘制出辅因子代谢中间物CI、底物以及产物其对应的浓度随时间的变化曲线图,得出辅因子代谢中间物CI在体内代谢反应中发挥着全局性的作用,而且辅因子代谢中间物CI在胞内的摩尔含量位于不太高的稳态,胞内代谢与胞外代谢相比的产物得率至少100-1000倍;
步骤三,针对含有非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的代谢网络系统在零抑制、低抑制、中抑制和高抑制四种抑制水平下进行胞内和胞外的对比干实验,并且绘制出非辅因子代谢中间物NCI、辅因子代谢中间物CI、底物以及产物其对应浓度随时间的变化曲线图;
步骤四,通过将步骤三条件下得出的多幅浓度随时间的变化曲线图进行比较,得出如下结论:(1)引入抑制后的胞内反应的非辅因子代谢中间物NCI都能维持稳态,稳态水平的高低由抑制强度的高低决定;(2)辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素;综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可;
步骤五,在高、中、低三种抑制水平下,通过调整CI在系统中的摩尔含量,得出CI的浓度高低对NCI稳态水平的影响,得出CI浓度越低,NCI稳态时的水平越低;
步骤六,针对步骤三的实验,通过调整该反应过程中酶的含量来模拟基因敲入、敲出实验,基因敲出实验能够通过设置酶量为零来模拟,得出酶量对非辅因子代谢中间物NCI的稳态摩尔含量和反应速率具有影响,而且能够定量的计算出代谢流对于不同酶的敏感程度;
步骤七,根据步骤二、步骤四至步骤六得出的结论,得出与胞内反应相一致的胞外代谢反应系统,具体如下:
其中,S1、S2代表底物,P1、P2代谢产物,B和Q代表辅因子代谢物CI,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7;
步骤八,根据所述步骤七中得出的胞外代谢反应系统,将数据库内的相关体内代谢的动力学参数代入胞外代谢反应系统内,在计算机上进行大量的干反应,最终得出一组非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量,使其性能最好,然后在实验室进行湿实验,对胞外代谢反应系统进行微调,得出最优的非辅因子代谢中间物NCI、辅因子代谢中间物CI、酶、底物和产物的加入量值,即为最优的胞外代谢反应系统。
2.根据权利要求1所述的体外再现体内反应的方法,其特征在于,所述步骤一的具体步骤为:
已知在体内代谢网络中,双底物反应的比例为70-80%,第一、第二底物都过量,那么在体内代谢网络的反应中能够得到2个底物对酶的亲和性常数Km、1个酶的催化常数kcat以及最大反应速度Vmax;
Vmax=e1,0×kcat;
双底物的不可逆反应式为:S1+B→P1+Q,该反应中包含的基元反应为:
其中,k1,1,k1,-1,k1,2,k1,-2,k1,3,k1,-3k1,4,k1,-4均为酶反应的速率常数,S1和B表示反应的底物,P1和Q表示反应的产物,E1表示酶,E1B、E1BS1、E1Q表示酶复合物;
根据酶定量守恒定量,总酶量e1,0满足如下公式:
e1,0=e1+e1b+e1bs1+e1q
其中,e1代表E1的摩尔量,e1b代表E1B的摩尔量,e1bs1代表E1BS1的摩尔量,e1q代表E1Q的摩尔量;
由于上述双底物反应式中和为不可逆反应,因此,k1,-3=k1,-4=0;
P1的生成速率的具体计算公式为:
<mrow>
<msub>
<mi>r</mi>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>b</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mi>b</mi>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>bs</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,代表P1的生成速率,s1代表底物S1的摩尔量,b代表底物B的摩尔量,由公式(1)可知,除了由底物B和S1确定外,还由k1,1,k1,-1,k1,2,k1,-2,k1,3,k1,4和e1,0确定;
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Kmb表示底物B对酶的亲和性常数,并且Kmb越大,代表底物B对酶的亲和性越低;
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,kcat表示酶的催化常数,并且kcat越大,代表没得催化效率越高;
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Kms1表示底物S1对酶的亲和性常数,并且Kms1越大,代表底物S1对酶的亲和性越低;
首先,假设s1>>b,同时将公式(2)和(3)代入公式公式(1)中,得出:
<mrow>
<msub>
<mi>r</mi>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&cong;</mo>
<mfrac>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mi>b</mi>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>b</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
然后,将公式(3)除以公式(2),得出k1,1,具体如下:
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
由公式(6)可知,k1,1为固定值;
进一步地,假设b>>s1,同时将公式(2)和(3)代入公式公式(1)中,得出:
<mrow>
<msub>
<mi>r</mi>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&cong;</mo>
<mfrac>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
进一步地,将公式(2)、(3)和(4)代入公式(1)中,得出:
<mrow>
<msub>
<mi>r</mi>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>b</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>b</mi>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>bs</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
由公式(8)可知,除了由底物B和S1确定外,还由总酶量e1,0、酶的催化常数kcat、底物对酶的亲和性常数Kms1、Kmb以及酶反应的速率常数k1,-1确定;因此,在给定底物浓度s1、b以及产物的生成速率时,能够得出e1,0,kcat,Kms1,Kmb,andk1,-1的关系式,进而从计算机的数据库内提取出满足e1,0,kcat,Kms1,Kmb,andk1,-1的关系式酶反应动力学参数;
由公式(3)可知,在kcat求得后,通过任意调整k1,3和k1,4的值,使其满足的定值即可;
将公式(4)除以公式(3),得出:
<mrow>
<mfrac>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
由公式(9)可知,在满足Kms1和kcat一定的条件下,通过任意调整k1,-2和k1,3的值,使其满足的定值即可;
因此,由公式(1)至(9)可知,对于双底物反应,在得到米氏动力学常数kcat、 Kms1、Kmb,能够求得k1,1、k1,-1、k1,2、k1,-2、k1,3、k1,-3、k1,4和k1,-4;同理,对于单底物或者三底物反应,当给定的动力学参数的个数少于需要的速率常数的个数时,能够按照公式(1)至(9)的方法,找到数据库参数和速率常数之间的联系,然后求得各个速率常数k,再完成所述步骤二的步骤。
3.根据权利要求2所述的体外再现体内反应的方法,其特征在于,所述步骤二的具体步骤为:
(A)胞内代谢网络与基元反应式:
(B)胞外代谢网络与基元反应式:
其中,(A)胞内代谢网络及与它相关的基元反应中,B和Q为一组辅因子代谢物CI,他们在网络中不断消耗和产生,且不会耗尽也不会被排出网络,因此B+Q是一个含量很低的常数;(B)胞外代谢网络及与它相关的基元反应中,(1)中的B和S1为底物,Q和P1为相应的产物;(2)中的Q和S2为底物,B和P2为产物;并且(1)中产生的Q不能被(2)作为底物使用,(2)中产生的B也不能被(1)做底物使用;
为了让(A)胞内代谢网络和(B)胞外代谢网络具有可比性,将S1-P1和S2-P2的反应线路设计成完全相同的;即:S1=S2,P1=P2,E1=E2,B=Q,k1,j=k2,j,k1,-j=k2,-j,j=1,2,3,4;因此对于胞内反应来说,根据步骤一中获得的动力学参数进行筛选,得出最小代谢网络的相关动力学参数表格;
在阐述胞内代谢网络相比于胞外代谢网络的优势之前,首先将引入四个定义:反应速率、产物得率、相对产物得率和相对反应速率;其中,反应速率分为产物反应速率和底物反应速率基于CI的产物得率为qCI(s-1),胞内反应系统的基于CI的产物得率为胞外反应系统的基于CI的产物得率为基于CI的相对产物得率为Rq/q,Rq/q是客观衡量胞内反应相对于胞外反应的效率的标准;CI的相对反应速率为Rr/r,Rr/r是客观衡量胞内反应相对与胞外反应的速率的标准,具体计算公式如下:
<mrow>
<msub>
<mi>r</mi>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>r</mi>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;s</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>;</mo>
</mrow>
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Delta;p</mi>
<mi>i</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>&Delta;t</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Delta;p</mi>
<mi>i</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>&Delta;t</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>b</mi>
<mn>0</mn>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>q</mi>
<mo>/</mo>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>q</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msubsup>
<msubsup>
<mi>q</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mo>/</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>r</mi>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>v</mi>
<mi>o</mi>
</mrow>
</msubsup>
<msubsup>
<mi>r</mi>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>-</mo>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mfrac>
<mo>;</mo>
</mrow>
其中,in-vivo表示胞内反应,in-vitro表示胞外反应,Δt表示反应时间,Δpi表示从产物积累的角度衡量反应的进程,Δsi表示从底物消耗的角度衡量反应的进程,bss是B在稳态时的值,b0表示B与S1同时耗尽;
根据所述最小代谢网络的相关动力学参数表格中的数值以及反应物、产物、辅因子代谢中间物CI与时间的变化曲线进行比较,得出如下结论:(1)胞内和胞外代谢网络的反应时间Δt、底物的消耗量ΔS和产物的生成量ΔP完全相同,因此胞内反应相比于胞外反应的基于CI的产物得率根据上述反应结果,表明反应速率中的快反应和慢反应,虽然会影响反应时间,但不会影响基于CI的相对产物得率Rq/q,而底物的初始值会对相对反应速率有影响;(2)胞内和胞外代谢网络中B和Q随着S1和S2的消耗而不断的消耗,因此在这样的胞外系统中用B和Q的初始值来计算基于CI的产物得率更加准确合理;尽管这个最小代谢反应不是真实的,但是它展示了辅因子代谢中间物CI一定是发挥着全局和局部的作用,而且他们在胞内位于低浓度的稳态,胞内代谢相比与胞外代谢的产物得率至少为100-1000倍。
4.根据权利要求2所述的体外再现体内反应的方法,其特征在于,所述步骤三的具体步骤为:
(C)含有NCI、无抑制胞内代谢网络反应过程为:
(D)含有NCI、无抑制胞外代谢网络反应过程为:
(C)和(D)中含有NCI、无抑制胞内和胞外代谢网络的基元反应式:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此在胞外代谢网络中B和Q仅仅作为中间物出现,不影响S1、S2、P1、P2的反应和性质,不具有辅因子代谢物CI的全局性作用;
(E)含有NCI和产物抑制的胞内代谢网络的反应过程为:
(F)含有NCI和产物抑制的胞外代谢网络的反应过程为:
(E)和(F)中含有NCI和产物抑制的胞内和胞外代谢网络的基元反应为:
其中,S1、S2代表底物,P1、P2代谢产物,Y1、Y2、Y3、Y4、Y5、Y6代表非辅因子代谢物NCI,E1、E2、E3、E4、E5、E6、E7、E8为相应的酶,Y1抑制E1,Y2抑制E2,Y3抑制E3,Y4抑制E5,Y5抑制E6,Y6抑制E7,B和Q代表辅因子代谢物CI,但是在胞外代谢网络(1)产生的Q不能被(2)作为底物使用,(2)产生的B也不能被(1)作为底物使用,因此胞外系统中B和Q,不具有辅因子代谢物CI的全局性作用;
对于胞内反应来说,根据步骤一中获得的动力学参数进行筛选,得出含有NCI的相关动力学参数表格;
将所述步骤二的最小代谢网络与步骤三中的代谢网络相比较,得知:在最小代谢网络中,由于没有NCI,总是有ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro,然而在目前复杂的模型中不能使ΔPin vivo=ΔPin vitro且Δtin vivo=Δtin vitro同时成立,否则此时胞内反应中B和Q的摩尔含量将高达30mM左右;在真实的细胞内,渗透压的原因以及细胞摩尔含量的可调性决定了CI的摩尔含量不会太高,大肠内CI的平均含量是1.5mM;因此我们采取了降低反应速度或者降低CI的摩尔含量的方法,从而延长反应时间的方法来评价Rq/q;
由含有NCI的相关动力学参数表格得出,产物抑制一方面减小了产物得率,另一方面也减小了反应的速度;然而,当减小辅因子代谢中间物CI在代谢网络中的含量时,基于CI的相对产物得率却明显上升,虽然这其中伴随值反应速率的略微减小,但是这一损失对细胞来说是划算的;较低的CI摩尔含量提高了代谢的效率,增强了代谢物摩尔含量对代谢流的敏感性,降低了反应速率从而保证其他每一个代谢中间物的摩尔含量都被控制在较低的水平;总的来说这一结果与我们在最小代谢网络中的结果完全一致,从而进一步巩固了CI对与细胞维持高效代谢活动的全局性作用;
进一步地,根据参数绘制在胞内以及胞外代谢网络中非辅因子代谢中间物NCI、反应物和产物的对应物质的量随时间的变化曲线图;
5.根据权利要求1所述的体外再现体内反应的方法,其特征在于,所述步骤四的具体步骤为:
在没有抑制的条件下,非辅因子代谢中间物NCI有的维持稳态,有的没有维持稳态;但是当引入抑制后,所有胞内反应的非辅因子代谢中间物NCI都能维持稳态;稳态水平的高低则由抑制强度的高低来决定;
比较胞内反应和胞外反应,虽然胞外反应与对应的胞内反应有相同的抑制机制、抑制水平,但是并不是所有的NCI都有稳态,考虑此时胞内网络和胞外反应唯一的区别是前者包含CI而后者虽然包含CI,但是不具有CI的性质,相当于胞外网络不具有CI,因此唯一的区别在于胞内含有CI,因此得出结论:辅因子代谢中间物CI的存在是非辅因子代谢中间物NCI能否维持稳态的不可或缺的因素,而且CI摩尔含量越低,NCI稳态时的摩尔含量水平越低,综上所述,抑制机制和CI的全局性作用是胞内网络维持稳态的两个重要因素,两者缺一不可。
6.根据权利要求1所述的体外再现体内反应的方法,其特征在于,所述步骤五的具体步骤为:在高、中、低三种抑制水平下以CI的浓度大小为自变量,得出CI浓度的改变对胞内系统中其他NCI稳态水平的影响;得出如下结论:NCI总体对CI的改变是很敏感的,当CI浓度很低时,NCI对CI的敏感度最强,因此能够通过降低CI的浓度,使得NCI水平降低,从而提高产物得率。
7.根据权利要求1所述的体外再现体内反应的方法,其特征在于,所述步骤六的具体步骤为:
绘制单底物和双底物中酶反应的酶含量对胞内网络中非辅因子代谢中间物NCI的摩尔含量随着时间的变化曲线图,得出在胞内网路中,酶量的影响是全局性的;增加酶量会使得前面的非辅因子代谢中间物NCI的稳态浓度下降,后面的非辅因子代谢中间物NCI的稳态浓度上升;同时由于B+Q是常数,降低B会使Q升高,从而将酶量的影响由局部扩大到全局;
酶量对反应速率会产生影响,酶量越大,反应速率越大,但是反应速率对不同酶的敏感程度不一样,因此,酶量不管是反应速率还是对非辅因子代谢中间物NCI的稳态浓度的影响都是全局性的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510926583.5A CN105574358B (zh) | 2015-12-14 | 2015-12-14 | 体外再现体内反应的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510926583.5A CN105574358B (zh) | 2015-12-14 | 2015-12-14 | 体外再现体内反应的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105574358A CN105574358A (zh) | 2016-05-11 |
CN105574358B true CN105574358B (zh) | 2018-01-23 |
Family
ID=55884483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510926583.5A Expired - Fee Related CN105574358B (zh) | 2015-12-14 | 2015-12-14 | 体外再现体内反应的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105574358B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021787B (zh) * | 2017-12-08 | 2021-10-15 | 华东理工大学 | 一种根据基元反应拓扑结构构建动态代谢网络的方法 |
AU2020254517A1 (en) * | 2019-03-29 | 2021-09-30 | Amgen Inc. | Predicting cell culture performance in bioreactors |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1419472A2 (en) * | 2001-08-16 | 2004-05-19 | Biotech Research Ventures Pte Ltd | Method for modelling biochemical pathways |
-
2015
- 2015-12-14 CN CN201510926583.5A patent/CN105574358B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
Bioengineering novel in vitro metabolic pathways using synthetic biology;Andreas Meyer et.al;《Current Opinion in Microbiology》;20071231;第2007年卷;全文 * |
Biohydrogenation from Biomass Sugar Mediated by In Vitro Synthetic Enzymatic Pathways;Yiran Wang et.al;《Chemistry & Biology》;20110325;第2011年卷;全文 * |
Joseph Rollin et.al.High-yield hydrogen production from biomass by in vitro metabolic engineering: Mixed sugars coutilization and kinetic modeling.《Proceedings of the National Academy of Sciences》.2015,第2015年卷 * |
Metabolic Control through Reflexive Enzyme Action;L.KOCH;《J.Theoret.Biol.》;19671231;第1967年卷(第15期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105574358A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Boynton et al. | Ions, cell proliferation, and cancer | |
Boitard et al. | Monitoring single-cell bioenergetics via the coarsening of emulsion droplets | |
Westerhoff et al. | Systems biology towards life in silico: mathematics of the control of living cells | |
Guo et al. | Consistent isomiR expression patterns and 3′ addition events in miRNA gene clusters and families implicate functional and evolutionary relationships | |
CN105574358B (zh) | 体外再现体内反应的方法 | |
Westerhoff et al. | On the control of gene expression | |
Kostenbader et al. | Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme | |
Di Caprio | A fattening factor to quantify the accumulation ability of microorganisms under N-starvation | |
Kabimoldayev et al. | Basics of genome-scale metabolic modeling and applications on C1-utilization | |
Zhang et al. | General theory for integrated analysis of growth, gene, and protein expression in biofilms | |
Ji et al. | Growth characteristics of freeze-tolerant baker’s yeast Saccharomyces cerevisiae AFY in aerobic batch culture | |
Lu et al. | Combined analysis reveals a core set of cycling genes | |
Machado et al. | A kinetic model for curcumin production in Escherichia coli | |
Trausch et al. | Design of modular “plug-and-play” expression platforms derived from natural riboswitches for engineering novel genetically encodable RNA regulatory devices | |
Byrne et al. | Characterization and adaptation of Caldicellulosiruptor strains to higher sugar concentrations, targeting enhanced hydrogen production from lignocellulosic hydrolysates | |
Zhu et al. | Use of flux balance analysis to promote lipid productivity in Chlorella sorokiniana | |
Khoodeeram et al. | An Ockham Razor model of energy metabolism | |
Jariwala et al. | Transport equations in an enzymatic glucose fuel cell | |
Wilkins | Factors influencing microbial growth and its kinetics | |
Arnau et al. | A feedback mechanism controls rDNA copy number evolution in yeast independently of natural selection | |
CN102629304B (zh) | 基于基因组尺度代谢网络模型的代谢工程设计预测方法 | |
Kaplan et al. | Modeling and simulation of metabolic networks for estimation of biomass accumulation parameters | |
Zakharov et al. | Synchronization of circadian rhythms at scale of gene, cell and whole organism | |
Ramakrishnan et al. | Redescription mining and applications in bioinformatics | |
Thi Doan et al. | A Monte‐Carlo approach for the simulation of microbial population dynamics in an heterogeneous scale‐down bioreactor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180123 Termination date: 20191214 |