Centro de Engenharia Biológica, Universidade do Minho, Braga, Portugal.
SilicoLife Lda, Braga, Portugal.
BMC Bioinformatics. 2019 Jun 20;20(1):350. doi: 10.1186/s12859-019-2934-y.
Computational strain optimisation methods (CSOMs) have been successfully used to exploit genome-scale metabolic models, yielding strategies useful for allowing compound overproduction in metabolic cell factories. Minimal cut sets are particularly interesting since their definition allows searching for intervention strategies that impose strong growth-coupling phenotypes, and are not subject to optimality bias when compared with simulation-based CSOMs. However, since both types of methods have different underlying principles, they also imply different ways to formulate metabolic engineering problems, posing an obstacle when comparing their outputs.
In this work, we perform an in-depth analysis of potential strategies that can be obtained with both methods, providing a critical comparison of performance, robustness, predicted phenotypes as well as strategy structure and size. To this end, we devised a pipeline including enumeration of strategies from evolutionary algorithms (EA) and minimal cut sets (MCS), filtering and flux analysis of predicted mutants to optimize the production of succinic acid in Saccharomyces cerevisiae. We additionally attempt to generalize problem formulations for MCS enumeration within the context of growth-coupled product synthesis. Strategies from evolutionary algorithms show the best compromise between acceptable growth rates and compound overproduction. However, constrained MCSs lead to a larger variety of phenotypes with several degrees of growth-coupling with production flux. The latter have proven useful in revealing the importance, in silico, of the gamma-aminobutyric acid shunt and manipulation of cofactor pools in growth-coupled designs for succinate production, mechanisms which have also been touted as potentially useful for metabolic engineering.
The two main groups of CSOMs are valuable for finding growth-coupled mutants. Despite the limitations in maximum growth rates and large strategy sizes, MCSs help uncover novel mechanisms for compound overproduction and thus, analyzing outputs from both methods provides a richer overview on strategies that can be potentially carried over in vivo.
计算应变优化方法(CSOMs)已成功用于开发基因组规模的代谢模型,为代谢细胞工厂中的化合物过量生产提供了有用的策略。最小割集特别有趣,因为它们的定义允许搜索施加强生长偶联表型的干预策略,并且与基于模拟的 CSOM 相比,不受最优性偏差的影响。然而,由于这两种方法具有不同的基本原理,它们也意味着制定代谢工程问题的不同方式,当比较它们的输出时,这构成了一个障碍。
在这项工作中,我们对这两种方法都可以获得的潜在策略进行了深入分析,对性能、鲁棒性、预测表型以及策略结构和大小进行了批判性比较。为此,我们设计了一个包括从进化算法(EA)和最小割集(MCS)枚举策略、预测突变体通量分析以优化酿酒酵母琥珀酸生产的流水线。我们还试图在生长偶联产物合成的背景下对 MCS 枚举问题进行形式化处理。进化算法的策略在可接受的生长速率和化合物过量生产之间表现出最佳折衷。然而,受约束的最小割集导致具有多种生长耦联程度和生产通量的表型多样性更大。后者已被证明在揭示γ-氨基丁酸支路的重要性方面是有用的,并且在生长偶联设计中操纵辅因子池对于琥珀酸生产具有潜在的有用性,这些机制也被吹捧为代谢工程的潜在有用性。
两种主要的 CSOM 组对于寻找生长偶联突变体很有价值。尽管存在最大生长速率和大策略尺寸的限制,最小割集有助于揭示化合物过量生产的新机制,因此,分析两种方法的输出提供了对体内潜在策略的更丰富的概述。