Accelerating metabolic models evaluation with statistical metamodels: application to Salmonella infection models

Clémence Frioux (1) , Sylvie Huet (2) , Simon Labarthe (1, 3, 2) , Julien Martinelli (4, 5) , Thibault Malou (6) , David James Sherman (1) , Marie-Luce Taupin (7) , Pablo Ugalde-Salas (1)

1 PLEIADE - Pleiade, from patterns to models in computational biodiversity and biotechnology
2 MaIAGE - Mathématiques et Informatique Appliquées du Génome à l'Environnement [Jouy-En-Josas]
3 BioGeCo - Biodiversité, Gènes & Communautés
4 CBIO - Centre de Bioinformatique
5 Lifeware - Computational systems biology and optimization
6 INSA Toulouse - Institut National des Sciences Appliquées - Toulouse
7 LaMME - Laboratoire de Mathématiques et Modélisation d'Evry

Mathematical and numerical models are increasingly used in microbial ecology to model the fate of microbial communities in their ecosystem. These models allow to connect in a mechanistic framework species-level informations, such as the microbial genomes, with macro-scale features, such as species spatial distributions or metabolite gradients. Numerous models are built upon specieslevel metabolic models that predict the metabolic behaviour of a microbe by solving an optimization problem knowing its genome and its nutritional environment. However, screening the community dynamics with these metabolic models implies to solve such an optimization problem by species at each time step, leading to a significant computational load further increased by several orders of magnitude when spatial dimensions are added. In this paper, we propose a statistical framework based on Reproducing Kernel Hilbert Space (RKHS) metamodels that are used to provide fast approximations of the original metabolic model. The metamodel can replace the optimization step in the system dynamics, providing comparable outputs at a much lower computational cost. We will first build a system dynamics model of a simplified gut microbiota composed of a unique commensal bacterial strain in interaction with the host and challenged by a Salmonella infection. Then, the machine learning method will be introduced, and particularly the ANOVA-RKHS that will be exploited to achieve variable selection and model parsimony. A training dataset will be constructed with the original system dynamics model and hyper-parameters will be carefully chosen to provide fast and accurate approximations of the original model. Finally, the accuracy of the trained metamodels will be assessed, in particular by comparing the system dynamics outputs when the original model is replaced by its metamodel. The metamodel allows an overall relative error of 4.71% but reducing the computational load by a speed-up factor higher than 45, while correctly reproducing the complex behaviour occurring during Salmonella infection. These results provide a proof-of-concept of the potentiality of machine learning methods to give fast approximations of metabolic model outputs and pave the way towards PDEbased spatio-temporal models of microbial communities and host-microbiota-pathogen interactions.

Domaines: Modélisation et simulation

Clémence Frioux, Sylvie Huet, Simon Labarthe, Julien Martinelli, Thibault Malou, et al.. Accelerating metabolic models evaluation with statistical metamodels: application to Salmonella infection models. 2022. ⟨hal-03635862v1⟩