The aldehyde dehydrogenase, AldA, is essential for L-1,2-propanediol utilization in laboratory-evolved Escherichia coli.

Most Escherichia coli strains are naturally unable to grow on 1,2-propanediol (PDO) as a sole carbon source. Recently, however, a K-12 descendent E. coli strain was evolved to grow on 1,2-PDO, and it was hypothesized that this evolved ability was dependent on the aldehyde dehydrogenase, AldA, which is highly conserved among members of the family Enterobacteriacea. To test this hypothesis, we first performed computational model simulation, which confirmed the essentiality of the aldA gene for 1,2-PDO utilization by the evolved PDO-degrading E. coli. Next, we deleted the aldA gene from the evolved strain, and this deletion was sufficient to abolish the evolved phenotype. On re-introducing the gene on a plasmid, the evolved phenotype was restored. These findings provide experimental evidence for the computationally predicted role of AldA in 1,2-PDO utilization, and represent a good example of E. coli robustness, demonstrated by the bacterial deployment of a generalist enzyme (here AldA) in multiple pathways to survive carbon starvation and to grow on a non-native substrate when no native carbon source is available.

Model-driven discovery of synergistic inhibitors against E. coli and S. enterica serovar Typhimurium targeting a novel synthetic lethal pair, aldA and prpC.

Mathematical models of biochemical networks form a cornerstone of bacterial systems biology. Inconsistencies between simulation output and experimental data point to gaps in knowledge about the fundamental biology of the organism. One such inconsistency centers on the gene aldA in Escherichia coli: it is essential in a computational model of E. coli metabolism, but experimentally it is not. Here, we reconcile this disparity by providing evidence that aldA and prpC form a synthetic lethal pair, as the double knockout could only be created through complementation with a plasmid-borne copy of aldA. Moreover, virtual and biological screening against the two proteins led to a set of compounds that inhibited the growth of E. coli and Salmonella enterica serovar Typhimurium synergistically at 100-200 µM individual concentrations. These results highlight the power of metabolic models to drive basic biological discovery and their potential use to discover new combination antibiotics.

Systems biology-guided identification of synthetic lethal gene pairs and its potential use to discover antibiotic combinations.

Mathematical models of metabolism from bacterial systems biology have proven their utility across multiple fields, for example metabolic engineering, growth phenotype simulation, and biological discovery. The usefulness of the models stems from their ability to compute a link between genotype and phenotype, but their ability to accurately simulate gene-gene interactions has not been investigated extensively. Here we assess how accurately a metabolic model for Escherichia coli computes one particular type of gene-gene interaction, synthetic lethality, and find that the accuracy rate is between 25% and 43%. The most common failure modes were incorrect computation of single gene essentiality and biological information that was missing from the model. Moreover, we performed virtual and biological screening against several synthetic lethal pairs to explore whether two-compound formulations could be found that inhibit the growth of Gram-negative bacteria. One set of molecules was identified that, depending on the concentrations, inhibits E. coli and S. enterica serovar Typhimurium in an additive or antagonistic manner. These findings pinpoint specific ways in which to improve the predictive ability of metabolic models, and highlight one potential application of systems biology to drug discovery and translational medicine.