A Holistic Approach from Systems Biology Reveals the Direct Influence of the Quorum-Sensing Phenomenon on Pseudomonas aeruginosa Metabolism to Pyoverdine Biosynthesis

Computational modeling and simulation of biological systems have become valuable tools for understanding and predicting cellular performance and phenotype generation. This work aimed to construct, model, and dynamically simulate the virulence factor pyoverdine (PVD) biosynthesis in through a systemic approach, considering that the metabolic pathway of PVD synthesis is regulated by the quorum-sensing (QS) phenomenon. The methodology comprised three main stages: (i) Construction, modeling, and validation of the QS gene regulatory network that controls PVD synthesis in strain PAO1; (ii) construction, curating, and modeling of the metabolic network of using the flux balance analysis (FBA) approach; (iii) integration and modeling of these two networks into an integrative model using the dynamic flux balance analysis (DFBA) approximation, followed, finally, by an in vitro validation of the integrated model for PVD synthesis in as a function of QS signaling. The QS gene network, constructed using the standard System Biology Markup Language, comprised 114 chemical species and 103 reactions and was modeled as a deterministic system following the kinetic based on mass action law. This model showed that the higher the bacterial growth, the higher the extracellular concentration of QS signal molecules, thus emulating the natural behavior of PAO1. The metabolic network model was constructed based on the iMO1056 model, the PAO1 strain genomic annotation, and the metabolic pathway of PVD synthesis. The metabolic network model included the PVD synthesis, transport, exchange reactions, and the QS signal molecules. This metabolic network model was curated and then modeled under the FBA approximation, using biomass maximization as the objective function (optimization problem, a term borrowed from the engineering field). Next, chemical reactions shared by both network models were chosen to combine them into an integrative model. To this end, the fluxes of these reactions, obtained from the QS network model, were fixed in the metabolic network model as constraints of the optimization problem using the DFBA approximation. Finally, simulations of the integrative model (CCBM1146, comprising 1123 reactions and 880 metabolites) were run using the DFBA approximation to get (i) the flux profile for each reaction, (ii) the bacterial growth profile, (iii) the biomass profile, and (iv) the concentration profiles of metabolites of interest such as glucose, PVD, and QS signal molecules. T