Kumar Vinay Satish, Maranas Costas D
Department of Industrial and Manufacturing Engineering, The Pennsylvania State University, University Park, Pennsylvania, United States of America.
PLoS Comput Biol. 2009 Mar;5(3):e1000308. doi: 10.1371/journal.pcbi.1000308. Epub 2009 Mar 13.
Genome-scale metabolic reconstructions are typically validated by comparing in silico growth predictions across different mutants utilizing different carbon sources with in vivo growth data. This comparison results in two types of model-prediction inconsistencies; either the model predicts growth when no growth is observed in the experiment (GNG inconsistencies) or the model predicts no growth when the experiment reveals growth (NGG inconsistencies). Here we propose an optimization-based framework, GrowMatch, to automatically reconcile GNG predictions (by suppressing functionalities in the model) and NGG predictions (by adding functionalities to the model). We use GrowMatch to resolve inconsistencies between the predictions of the latest in silico Escherichia coli (iAF1260) model and the in vivo data available in the Keio collection and improved the consistency of in silico with in vivo predictions from 90.6% to 96.7%. Specifically, we were able to suggest consistency-restoring hypotheses for 56/72 GNG mutants and 13/38 NGG mutants. GrowMatch resolved 18 GNG inconsistencies by suggesting suppressions in the mutant metabolic networks. Fifteen inconsistencies were resolved by suppressing isozymes in the metabolic network, and the remaining 23 GNG mutants corresponding to blocked genes were resolved by suitably modifying the biomass equation of iAF1260. GrowMatch suggested consistency-restoring hypotheses for five NGG mutants by adding functionalities to the model whereas the remaining eight inconsistencies were resolved by pinpointing possible alternate genes that carry out the function of the deleted gene. For many cases, GrowMatch identified fairly nonintuitive model modification hypotheses that would have been difficult to pinpoint through inspection alone. In addition, GrowMatch can be used during the construction phase of new, as opposed to existing, genome-scale metabolic models, leading to more expedient and accurate reconstructions.
基因组规模的代谢重建通常通过将利用不同碳源的不同突变体的计算机模拟生长预测与体内生长数据进行比较来验证。这种比较会产生两种类型的模型预测不一致情况;要么模型预测生长但实验中未观察到生长(生长非生长不一致,即GNG不一致),要么模型预测不生长但实验显示生长(非生长生长不一致,即NGG不一致)。在此,我们提出了一个基于优化的框架GrowMatch,以自动协调GNG预测(通过抑制模型中的功能)和NGG预测(通过向模型中添加功能)。我们使用GrowMatch来解决最新的大肠杆菌(iAF1260)计算机模拟模型的预测与Keio文库中的体内数据之间的不一致,并将计算机模拟与体内预测的一致性从90.6%提高到96.7%。具体而言,我们能够为56/72个GNG突变体和13/38个NGG突变体提出恢复一致性的假设。GrowMatch通过在突变体代谢网络中提出抑制作用解决了18个GNG不一致情况。通过抑制代谢网络中的同工酶解决了15个不一致情况,其余23个与基因阻断对应的GNG突变体通过适当修改iAF1260的生物量方程得以解决。GrowMatch通过向模型中添加功能为5个NGG突变体提出了恢复一致性的假设,而其余8个不一致情况通过找出可能执行缺失基因功能的替代基因得以解决。在许多情况下,GrowMatch识别出了相当非直观的模型修改假设,这些假设仅凭检查很难确定。此外,GrowMatch可用于新的(而非现有的)基因组规模代谢模型的构建阶段,从而实现更便捷、准确的重建。