Reconstrucción, modelamiento y simulación de la red metabólica y de Quorum-Sensing implicadas en la regulación de un fenotipo específico en Pseudomonas aeruginosa

Abstract

Una meta importante en el periodo posterior al desarrollo de las técnicas de secuenciación de nueva generación es la de relacionar las secuencias anotadas de los genomas con las funciones fisiológicas de una célula. Es por esto que la Biología de Sistemas ha venido trabajando en el diseño de nuevas metodologías computacionales para la reconstrucción de redes metabólicas a escala genómica y para el modelamiento y la simulación dinámica de estos sistemas biológicos en busca de estudiar la regulación de los mecanismos biológicos para la expresión de fenotipos. Debido a la facilidad relativa para la obtención de datos, el tamaño de sus genomas, los costos asociados y el interés clínico entre otras razones, los microorganismos son el grupo en el cual mayor cantidad de redes metabólicas se han reconstruido. Entre los mecanismos que controlan la expresión genética, el Quorum-Sensing es relevante no solo desde la perspectiva de las ciencias básicas, sino que se considera un eslabón importante para lograr avances en el área biotecnológica. El fenómeno conocido como Quorum-Sensing (QS) se basa en la comunicación celular mediada por moléculas de señalización y se encarga de sincronizar la expresión de fenotipos en una comunidad bacteriana, por ejemplo, para controlar su metabolismo y sus funciones como comunidad. Tomando como modelo la red metabólica regulada por QS en Pseudomonas aeruginosa (PAO1) para la producción de un factor de virulencia, de la familia de los sideróforos, conocido como pioverdina (PVD), este trabajo buscó generar la reconstrucción, modelamiento y simulación dinámica de esta red. Para esto se empleó el análisis dinámico de balance de flujo (DFBA) en el modelaje de la red metabólica asociada al fenómeno de QS como una estrategia importante para dar solución al interrogante que se plantea en este trabajo: evidenciar en el modelo la influencia del QS de PAO1 en la síntesis de PVD. Para dar cumplimiento a los objetivos que aborda este trabajo, la metodología propuesta comprende tres grandes etapas: (i) Reconstrucción, modelamiento y validación de la red génica de QS que regula la síntesis de PVD en PAO1. (ii) Construcción, curación y modelamiento bajo la aproximación de FBA de la red metabólica de Pseudomonas aeruginosa. Y (iii) unión, modelamiento bajo la aproximación de DFBA y validación experimental in vitro de la red génica de Quorum-Sensing acoplada con la red metabólica de PAO1 para la síntesis de PVD. La red génica de QS que regula la síntesis de PVD en PAO1, se construyó sobre el estándar SBML, consta de 114 especies químicas y biológicas y 103 reacciones. La red de QS fue modelada como un sistema determinista siguiendo los parámetros de la ley de acción de masas. Los resultados mostraron que a medida que aumenta el crecimiento poblacional, aumenta la producción de moléculas señal de QS en el espacio extracelular emulando así el comportamiento natural de un cultivo bacteriano de PAO1. La reconstrucción de la red metabólica se realizó con base en el modelo iMO1056, la anotación del genoma PAO1 y la vía metabólica para la biosíntesis PVD. El modelo metabólico involucra las reacciones de biosíntesis y de transporte e intercambio de PVD y de las moléculas señal de QS. Se realizó la curación de la red y posteriormente se modeló bajo la aproximación de DFBA, empleando como función objetivo la maximización de biomasa. Se seleccionaron nueve reacciones compartidas por la red de QS y la red metabólica para la fusión de ambas redes. Los flujos de estas reacciones de la red de QS, fueron fijados en el sistema metabólico como restricciones del problema de optimización. Utilizando el DFBA, se realizaron las simulaciones del sistema para obtener (i) los perfiles de flujo para cada reacción, (ii) el perfil de crecimiento, (iii) el perfil de biomasa y (iv) los perfiles de concentración de metabolitos de interés tales como moléculas de señalización de QS, glucosa y PVD. La red metabólica propuesta consta de 1124 reacciones y 881 metabolitos (modelo CCBM1737). El modelamiento dinámico de la red metabólica acoplada a la red de QS de PAO1, permitió evidenciar que el fenómeno de QS ejerce una influencia directa sobre la expresión de diferentes fenotipos metabólicos de acuerdo con el cambio de la intensidad de la señal de QS. Este trabajo es el primer reporte de un modelo in silico de la red génica que comprende todos los sistemas de Quorum-Sensing acoplada con la red metabólica de Pseudomonas aeruginosa. (Texto tomado de la fuente)

An important goal in the period following the development of new generation sequencing techniques is to relate the annotated sequences of genomes to the physiological functions of a cell. That is why Systems Biology has been working on the design of new computational methodologies for the reconstruction of metabolic networks at genomic scale and for the modeling and dynamic simulation of these biological systems in order to study the regulation of biological mechanisms for the expression of phenotypes. Due to the relative ease of data collection, genome size, associated costs and clinical interest, among other reasons, microorganisms are the group in which most metabolic networks have been rebuilt. Among the mechanisms that control gene expression, the Quorum-Sensing is relevant not only from the perspective of the basic sciences but is also considered an important link to achieve advances in the biotechnological area. The phenomenon known as Quorum-Sensing (QS) is based on cell communication mediated by signaling molecules and is responsible for synchronizing the expression of phenotypes in a bacterial community, for example, to control its metabolism and functions as a community. Taking as a model the metabolic network regulated by QS in Pseudomonas aeruginosa (PAO1) for the production of a virulence factor, of the siderophore family, known as pyoverdine (or this purpose, the dynamic flow balance analysis (DFBA) was used in the modeling of the metabolic network associated with the QS phenomenon as an important strategy to provide a solution to the question posed in this paper: to highlight in the model the influence of the QS of PAO1 on the synthesis of PVD. In order to meet the objectives of this work, the proposed methodology comprises three main stages: (i) Reconstruction, modeling, and validation of the QS gene network that regulates the synthesis of ODP in ODP1. (ii) Construction, healing, and modeling under the ABF approach of the metabolic network of Pseudomonas aeruginosa. And (iii) binding, modeling under the DFBA approach and in vitro experimental validation of the Quorum-Sensing gene network coupled with the PAO1 metabolic network for PVD synthesis. The QS gene network that regulates the synthesis of PVD in PAO1 was built on the SBML standard, consisting of 114 chemical and biological species and 103 reactions. The QS network was modeled as a deterministic system following the parameters of the law of mass action. The results showed that as population growth increases, the production of QS signal molecules in the extracellular space increases, emulating the natural behavior of a bacterial culture of PAO1. The reconstruction of the metabolic network was carried out based on the iMO1056 model, the annotation of the PAO1 genome and the metabolic pathway for PVD biosynthesis. The metabolic model involves the reactions of biosynthesis and transport and exchange of PVD and QS signal molecules. The network was cured and then modeled under the DFBA approach, using biomass maximization as an objective function. Nine reactions shared by the QS network and the metabolic network were selected to merge the two networks. The flows of these reactions from the QS network were set in the metabolic system as restrictions of the optimization problem. Using DFBA, system simulations were performed to obtain (i) the flow profiles for each reaction, (ii) the growth profile, (iii) the biomass profile and (iv) the concentration profiles of metabolites of interest such as QS, glucose, and PVD signaling molecules. The proposed metabolic network consists of 1124 reactions and 881 metabolites (model CCBM1737). The dynamic modeling of the metabolic network coupled to the QS network of PAO1, allowed to show that the QS phenomenon has a direct influence on the expression of different metabolic phenotypes according to the change in the intensity of the QS signal. This paper is the first report of an in silico model of the gene network comprising all Quorum-Sensing systems coupled with the metabolic network of Pseudomonas aeruginosa.