Abstract
The human gut is home to a complex array of microorganisms interacting with the host and each other, forming a community known as the microbiome. This community has been linked to human health and disease, but understanding the underlying interactions is still challenging for researchers. Standard studies typically use high-throughput sequencing to analyze microbiome distribution in patient samples. Recent advancements in meta-omic data analysis have enabled computational modeling strategies to integrate this information into an in silico model. However, there is a need for improved parameter fitting and data integration features in microbial community modeling. This study proposes a novel alternative strategy utilizing state-of-the-art dynamic flux balance analysis (dFBA) to provide a simple protocol enabling accurate replication of abundance data composition through dynamic parameter estimation and integration of metagenomic data. We used a recurrent optimization algorithm to replicate community distributions from three different sources: mock, in vitro, and clinical microbiome. Our results show an accuracy of 98% and 96% when using in vitro and clinical bacterial abundance distributions, respectively. The proposed modeling scheme allowed us to observe the evolution of metabolites. It could provide a deeper understanding of metabolic interactions while taking advantage of the high contextualization features of GEM schemes to fit the study case. The proposed modeling scheme could improve the approach in cases where external factors determine specific bacterial distributions, such as drug intake.