Dynamical Systems-Constrained Metabolic Modeling Enables Forecasting of Host-Microbiome Dynamics

Forecasting how microbiome-host ecosystems evolve through time simultaneously at the compositional and functional level remains a central challenge in biology. While dynamical systems models (DSMs) can infer and predict community composition from longitudinal abundance data, and constraint-based metabolic models (CBMMs) can estimate metabolic fluxes from genome-scale reconstructions, no existing framework unifies these approaches to generate mechanistically grounded, time-resolved forecasts of both microbial abundances and metabolite dynamics from ecological data alone. Here, we introduce the Dynamical Systems Constrained Metabolic Modeling (DySCoMeMo) framework, a new hybrid computational pipeline that integrates ecological DSMs with CBMMs to predict temporal dynamics of biomass and metabolites across microbial communities and hosts. DySCoMeMo leverages parameters inferred from application of DSMs to microbiome time series data to constrain metabolic modeling over time, thereby bridging ecological interaction networks with genome-scale metabolic modeling. DySCoMeMo is able to predict future community and metabolite dynamics in vitro with accuracy that is superior or on-par compared to that achieved with established methods that require actual microbial abundances and/or metabolites data for metabolite network inference or for estimating the per-microbe contribution to the extracellular metabolic pool. DySCoMeMo also generalizes to in vivo data as it is capable of accurately forecasting microbial and metabolite dynamics in response to dietary perturbations even when host metabolism is included. Finally, DySCoMeMo uniquely enables the identification of keystone species by quantifying their contributions to sustaining metabolic environments. Together, our work establishes a generalizable, mechanistically grounded framework for time-resolved forecasting of microbiome-host microbial and metabolic dynamics, bridging ecological interaction inference with genome-scale metabolism of communities.