Kinetic Modeling of Plant Metabolism and Its Predictive Power: Peppermint Essential Oil Biosynthesis as an Example

解决的核心问题

植物代谢网络在细胞、组织及整体水平上受到区室化、细胞类型特化和昼夜节律等影响，处于非稳态，静态建模难以描述其动态变化。本文以薄荷精油生物合成为例，阐述如何构建一个基于实验数据的确定性动力学模型，以定量预测单萜代谢物浓度随时间的动态变化并解析通路调控机制。

研究策略

采用实验与建模迭代整合的策略，将宏观测量（腺毛产量、分布与油体积）与微观酶动力学（体外动力学常数、亚细胞区室体积分数）及发育动态（腺毛数量变化、酶活性变化）相结合，基于质量守恒建立常微分方程组，用MATLAB数值求解模拟单萜积累曲线，并通过统计检验评估模型拟合度。

研究路线图

100%

下载 .drawio 源文件

核心内容

该研究以薄荷精油生物合成途径为例，构建了基于实验数据的确定性动力学模型，用以定量预测单萜代谢物浓度随时间的动态变化并解析通路调控机制。研究将宏观测量（腺毛数量与分布、储油体积）与微观酶动力学参数、亚细胞区室体积分数及发育动态数据整合，基于质量守恒建立常微分方程组，用MATLAB数值求解模拟41天内的单萜积累曲线。模型的关键创新包括用逻辑函数模拟腺毛发育、高斯函数模拟酶活性随叶龄的变化，并引入Comp因子校正区室体积对酶有效浓度的影响，以及z因子和w因子分别校正(+)-薄荷呋喃与(+)-胡薄荷酮在分泌细胞与储油腔间的浓度差异。结果显示，成熟叶片含约10,151个腺毛，整叶精油量约1,535 μg；反馈抑制常数中GPP对IPPI的Ki为96 μM，(+)-薄荷呋喃对PR的Ki为300 μM。模拟的单萜谱与实验数据统计拟合良好，表明该模型能准确捕捉植物代谢网络在非稳态条件下的动态行为，为其他植物代谢途径的动态建模提供了可复用的方法学框架。

创新点

（1）利用逻辑函数模拟薄荷叶片腺毛发育动力学；（2）采用高斯函数（或高斯函数叠加）模拟酶活性随叶龄的动态变化；（3）引入Comp因子校正亚细胞区室体积对酶有效浓度的影响；（4）引入z因子和w因子分别校正(+)-薄荷呋喃及(+)-胡薄荷酮在分泌细胞与储油腔之间的浓度差异（z=100，w=0.05）；（5）提供了可直接复制运行的完整MATLAB代码。

研究对象（18）

底盘细胞	代谢通路	基因	蛋白	功能	说明
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	1-脱氧-D-木酮糖-5-磷酸合酶（DXS, E1）	催化	质体定位；利用Comp=0.139校正浓度；催化丙酮酸与GAP缩合生成DXP
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	1-脱氧-D-木酮糖-5-磷酸还原异构酶（DXR, E2）	催化	质体定位；可逆反应，正反向动力学常数分别采用KM2f与KM2r
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	2-C-甲基-D-赤藓醇-4-磷酸胞苷转移酶（MCT, E3）	催化	质体定位；催化ME4P与CTP生成CDP-ME
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	4-(胞苷5'-二磷酸)-2-C-甲基-D-赤藓醇激酶（CMK, E4）	催化	质体定位；催化CDP-ME磷酸化生成CDP-ME2P
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	2-C-甲基-D-赤藓醇-2,4-环二磷酸合酶（MECPS, E5）	催化	质体定位；催化CDP-ME2P环化生成MECPP
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	4-羟基-3-甲基丁-2-烯基二磷酸合酶（HDS, E6）	催化	质体定位；催化MECPP开环生成HMBPP
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	4-羟基-3-甲基丁-2-烯基二磷酸还原酶（HDR, E7）	催化	质体定位；产物比例为85% IPP和15% DMAPP，模型中按比例拆分为E7a和E7b
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	异戊烯基二磷酸异构酶（IPPI, E8）	催化	质体定位；受产物GPP竞争性抑制，Ki=96 μM；可逆反应，正反向动力学参数分别为KM8f与KM8r
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	香叶基二磷酸合酶（GPPS, E9）	催化	质体定位；假设随机顺序双底物机制，将DMAPP与IPP缩合生成GPP
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-柠檬烯合酶（LS, E10）	催化	质体定位；催化GPP环化生成(-)-柠檬烯
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-柠檬烯-3-羟化酶（L3H, E11）	催化	内质网膜定位（Comp=0.365）；细胞色素P450酶，催化(-)-柠檬烯羟化生成(-)-反式异胡椒烯醇
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-反式异胡椒烯醇脱氢酶（IsoDH, E12）	催化	线粒体定位（Comp=0.044）；催化(-)-trans-异胡椒烯醇氧化为(-)-异胡椒烯酮
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-异胡椒烯酮还原酶（IsOR, E13）	催化	细胞质定位（Comp=0.204）；催化(-)-异胡椒烯酮还原为(+)-顺式异胡薄荷酮
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(+)-顺式异胡薄荷酮异构酶（IsoI, E14）	催化	细胞质定位（Comp=0.204）；催化(+)-cis-异胡薄荷酮异构化为(+)-胡薄荷酮
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(+)-薄荷呋喃合酶（MFS, E15）	催化	内质网定位（Comp=0.365）；催化(+)-胡薄荷酮生成(+)-薄荷呋喃
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(+)-胡薄荷酮还原酶（PR, E16）	催化	细胞质定位（Comp=0.204）；产物90%（-）-薄荷酮、10%（+）-异薄荷酮；受(+)-薄荷呋喃竞争性抑制（Ki=300 μM）及底物(+)-胡薄荷酮底物抑制（Kis=112 μM），z因子和w因子分别校正相关浓度
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-薄荷酮：(-)-薄荷醇还原酶（MMR, E17）	催化	细胞质定位（Comp=0.204）；底物分别为(-)-薄荷酮（KM17a=3 μM）和(+)-异薄荷酮（KM17b=41 μM），产物分别为(-)-薄荷醇和(+)-新异薄荷醇
薄荷（Mentha × piperita）腺毛分泌细胞	薄荷单萜生物合成途径（MEP途径及下游修饰）	—	(-)-薄荷酮：(+)-新薄荷醇还原酶（MNR, E18）	催化	细胞质定位（Comp=0.204）；可逆反应，正反向动力学参数分别为KM18af/KM18ar；底物包括(-)-薄荷酮和(+)-异薄荷酮，产物包括(+)-新薄荷醇和(-)-异薄荷醇

关键发现

成熟薄荷腺毛的储油体积为2.03×10^-4 μl
成熟叶片含10,151个腺毛，按大小分布为大型3,959个（39%）、中型5,786个（57%）、小型406个（4%）
油密度0.9时整叶精油量约1,535 μg。亚细胞体积分数：质体/亮体0.139、线粒体0.044、细胞质0.204、内质网0.365。E7（HDR）产物比例为85% IPP和15% DMAPP，E16（PR）产物比例为90%（-）-薄荷酮和10%（+）-异薄荷酮。反馈抑制常数：GPP对IPPI的Ki=96 μM，(+)-薄荷呋喃对PR的Ki=300 μM，(+)-胡薄荷酮对PR的底物抑制常数Kis=112 μM
分泌细胞中(+)-薄荷呋喃浓度仅为储油腔的1%（z=100），(+)-胡薄荷酮保留比例为5%（w=0.05）。模型初始底物丙酮酸和3-磷酸甘油醛各为1900 μg/叶。使用ode15s求解器模拟41天（3456000秒）动态变化，模拟单萜谱与实验数据通过统计检验拟合良好。

研究结论

该研究建立的薄荷精油生物合成动力学模型整合了酶动力学、亚细胞区室化、发育动态及底物供给等多层面实验数据，能够准确模拟单萜积累的动态变化，为理解植物代谢网络的结构与调控提供了有力工具，也可为其他植物代谢途径的动态建模提供可复用的方法学框架。