Polycyclic aromatic hydrocarbon metabolic network in Mycobacterium vanbaalenii PYR-1

This study investigated a metabolic network (MN) from Mycobacterium vanbaalenii PYR-1 for polycyclic aromatic hydrocarbons (PAHs) from the perspective of structure, behavior, and evolution, in which multilayer omics data are integrated. Initially, we utilized a high-throughput proteomic analysis to assess the protein expression response of M. vanbaalenii PYR-1 to seven different aromatic compounds. A total of 3,431 proteins (57.38% of the genome-predicted proteins) were identified, which included 160 proteins that seemed to be involved in the degradation of aromatic hydrocarbons. Based on the proteomic data and the previous metabolic, biochemical, physiological, and genomic information, we reconstructed an experiment-based system-level PAH-MN. The structure of PAH-MN, with 183 metabolic compounds and 224 chemical reactions, has a typical scale-free nature. The behavior and evolution of the PAH-MN reveals a hierarchical modularity with funnel effects in structure/function and intimate association with evolutionary modules of the functional modules, which are the ring cleavage process (RCP), side chain process (SCP), and central aromatic process (CAP). The 189 commonly upregulated proteins in all aromatic hydrocarbon treatments provide insights into the global adaptation to facilitate the PAH metabolism. Taken together, the findings of our study provide the hierarchical viewpoint from genes/proteins/metabolites to the network via functional modules of the PAH-MN equipped with the engineering-driven approaches of modularization and rationalization, which may expand our understanding of the metabolic potential of M. vanbaalenii PYR-1 for bioremediation applications.