In silico analysis of design of experiment methods for metabolic pathway optimization.

Microbial cell factories allow the production of chemicals presenting an alternative to traditional fossil fuel-dependent production. However, finding the optimal expression of production pathway genes is crucial for the development of efficient production strains. Unlike sequential experimentation, combinatorial optimization captures the relationships between pathway genes and production, albeit at the cost of conducting multiple experiments. Fractional factorial designs followed by linear modeling and statistical analysis reduce the experimental workload while maximizing the information gained during experimentation. Although tools to perform and analyze these designs are available, guidelines for selecting appropriate factorial designs for pathway optimization are missing. In this study, we leverage a kinetic model of a seven-genes pathway to simulate the performance of a full factorial strain library. We compare this approach to resolution V, IV, III, and Plackett Burman (PB) designs. Additionally, we evaluate the performance of these designs as training sets for a random forest algorithm aimed at identifying best-producing strains. Evaluating the robustness of these designs to noise and missing data, traits inherent to biological datasets, we find that while resolution V designs capture most information present in full factorial data, they necessitate the construction of a large number of strains. On the other hand, resolution III and PB designs fall short in identifying optimal strains and miss relevant information. Besides, given the small number of experiments required for the optimization of a pathway with seven genes, linear models outperform random forest. Consequently, we propose the use of resolution IV designs followed by linear modeling in Design-Build-Test-Learn (DBTL) cycles targeting the screening of multiple factors. These designs enable the identification of optimal strains and provide valuable guidance for subsequent optimization cycles.

Combinatorial optimization of pathway, process and media for the production of p-coumaric acid by Saccharomyces cerevisiae.

Microbial cell factories are instrumental in transitioning towards a sustainable bio-based economy, offering alternatives to conventional chemical processes. However, fulfilling their potential requires simultaneous screening for optimal media composition, process and genetic factors, acknowledging the complex interplay between the organism's genotype and its environment. This study employs statistical design of experiments to systematically explore these relationships and optimize the production of p-coumaric acid (pCA) in Saccharomyces cerevisiae. Two rounds of fractional factorial designs were used to identify factors with a significant effect on pCA production, which resulted in a 168-fold variation in pCA titre. Moreover, a significant interaction between the culture temperature and expression of ARO4 highlighted the importance of simultaneous process and strain optimization. The presented approach leverages the strengths of experimental design and statistical analysis and could be systematically applied during strain and bioprocess design efforts to unlock the full potential of microbial cell factories.

Machine Learning-Guided Optimization of p-Coumaric Acid Production in Yeast.

Industrial biotechnology uses Design-Build-Test-Learn (DBTL) cycles to accelerate the development of microbial cell factories, required for the transition to a biobased economy. To use them effectively, appropriate connections between the phases of the cycle are crucial. Using p-coumaric acid (pCA) production in Saccharomyces cerevisiae as a case study, we propose the use of one-pot library generation, random screening, targeted sequencing, and machine learning (ML) as links during DBTL cycles. We showed that the robustness and flexibility of the ML models strongly enable pathway optimization and propose feature importance and Shapley additive explanation values as a guide to expand the design space of original libraries. This approach allowed a 68% increased production of pCA within two DBTL cycles, leading to a 0.52 g/L titer and a 0.03 g/g yield on glucose.