RESUMO
Genomics-based metabolic models of microorganisms currently have no easy way of corroborating predicted biomass with the actual metabolites being produced. This study uses untargeted mass spectrometry-based metabolomics data to generate a list of accurate metabolite masses produced from the human commensal bacteria Citrobacter sedlakii grown in the presence of a simple glucose carbon source. A genomics-based flux balance metabolic model of this bacterium was previously generated using the bioinformatics tool PyFBA and phenotypic growth curve data. The high-resolution mass spectrometry data obtained through timed metabolic extractions were integrated with the predicted metabolic model through a program called MS_FBA. This program correlated untargeted metabolomics features from C. sedlakii with 218 of the 699 metabolites in the model using an exact mass match, with 51 metabolites further confirmed using predicted isotope ratios. Over 1400 metabolites were matched with additional metabolites in the ModelSEED database, indicating the need to incorporate more specific gene annotations into the predictive model through metabolomics-guided gap filling.
RESUMO
WHO reported 10.4 million new tuberculosis (TB) cases and 1.8 million deaths in 2015, making M. tuberculosis the most successful human pathogen with highest mortality among infectious diseases [1,2]. Drug-resistant TB is a major threat to global TB control [2,3]. Recently Torres et al. [4] identified 14 novel substitutions in M. tuberculosis-KatG (the enzyme associated with resistance to isoniazid-an important first-line anti-TB drug) and demonstrated that 12 of the 14 can cause INH-resistance in M. smegmatis. This study presents an in silico structure-based analysis of these 14 amino acid substitutions using homology models and x-ray crystal structures (when available) in M. tuberculosis. Our models demonstrate that several of these mutations cluster around three openings in the KatG tertiary structure which appear to initiate channels to the heme group at the catalytic center of the enzyme. We studied the effects of these mutations on the tertiary structure of KatG, focusing on conformational changes in the three channels in the protein structure. Our results suggest that the 14 novel mutations sufficiently restrict one or more of these access channels, thus potentially preventing INH from reaching the catalytic heme. These observations provide valuable insights into the structure-based origins of INH resistance and provide testable hypotheses for future experimental studies.