Disrupting the ArcA Regulatory Network Amplifies the Fitness Cost of Tetracycline Resistance in Escherichia coli.

There is an urgent need for strategies to discover secondary drugs to prevent or disrupt antimicrobial resistance (AMR), which is causing >700,000 deaths annually. Here, we demonstrate that tetracycline-resistant (TetR) Escherichia coli undergoes global transcriptional and metabolic remodeling, including downregulation of tricarboxylic acid cycle and disruption of redox homeostasis, to support consumption of the proton motive force for tetracycline efflux. Using a pooled genome-wide library of single-gene deletion strains, at least 308 genes, including four transcriptional regulators identified by our network analysis, were confirmed as essential for restoring the fitness of TetR E. coli during treatment with tetracycline. Targeted knockout of ArcA, identified by network analysis as a master regulator of this new compensatory physiological state, significantly compromised fitness of TetR E. coli during tetracycline treatment. A drug, sertraline, which generated a similar metabolome profile as the arcA knockout strain, also resensitized TetR E. coli to tetracycline. We discovered that the potentiating effect of sertraline was eliminated upon knocking out arcA, demonstrating that the mechanism of potential synergy was through action of sertraline on the tetracycline-induced ArcA network in the TetR strain. Our findings demonstrate that therapies that target mechanistic drivers of compensatory physiological states could resensitize AMR pathogens to lost antibiotics. IMPORTANCE Antimicrobial resistance (AMR) is projected to be the cause of >10 million deaths annually by 2050. While efforts to find new potent antibiotics are effective, they are expensive and outpaced by the rate at which new resistant strains emerge. There is desperate need for a rational approach to accelerate the discovery of drugs and drug combinations that effectively clear AMR pathogens and even prevent the emergence of new resistant strains. Using tetracycline-resistant (TetR) Escherichia coli, we demonstrate that gaining resistance is accompanied by loss of fitness, which is restored by compensatory physiological changes. We demonstrate that transcriptional regulators of the compensatory physiologic state are promising drug targets because their disruption increases the susceptibility of TetR E. coli to tetracycline. Thus, we describe a generalizable systems biology approach to identify new vulnerabilities within AMR strains to rationally accelerate the discovery of therapeutics that extend the life span of existing antibiotics.

Reconsidering the in vivo functions of Clostridial Stickland amino acid fermentations.

Stickland amino acid fermentations occur primarily among species of Clostridia. An ancient form of metabolism, Stickland fermentations use amino acids as electron acceptors in the absence of stronger oxidizing agents and provide metabolic capabilities to support growth when other fermentable substrates, such as carbohydrates, are lacking. The reactions were originally described as paired fermentations of amino acid electron donors, such as the branched-chain amino acids, with recipients that include proline and glycine. We present a redox-focused view of Stickland metabolism following electron flow through metabolically diverse oxidative reactions and the defined-substrate reductase systems, including for proline and glycine, and the role of dual redox pathways for substrates such as leucine and ornithine. Genetic studies and Environment and Gene Regulatory Interaction Network (EGRIN) models for the pathogen Clostridioides difficile have improved our understanding of the regulation and metabolic recruitment of these systems, and their functions in modulating inter-species interactions within host-pathogen-commensal systems and uses in industrial and environmental applications.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME.

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

Predictive regulatory and metabolic network models for systems analysis of Clostridioides difficile.

We present predictive models for comprehensive systems analysis of Clostridioides difficile, the etiology of pseudomembranous colitis. By leveraging 151 published transcriptomes, we generated an EGRIN model that organizes 90% of C. difficile genes into a transcriptional regulatory network of 297 co-regulated modules, implicating genes in sporulation, carbohydrate transport, and metabolism. By advancing a metabolic model through addition and curation of metabolic reactions including nutrient uptake, we discovered 14 amino acids, diverse carbohydrates, and 10 metabolic genes as essential for C. difficile growth in the intestinal environment. Finally, we developed a PRIME model to uncover how EGRIN-inferred combinatorial gene regulation by transcription factors, such as CcpA and CodY, modulates essential metabolic processes to enable C. difficile growth relative to commensal colonization. The C. difficile interactive web portal provides access to these model resources to support collaborative systems-level studies of context-specific virulence mechanisms in C. difficile.

In vivo commensal control of Clostridioides difficile virulence.

Leveraging systems biology approaches, we illustrate how metabolically distinct species of Clostridia protect against or worsen Clostridioides difficile infection in mice by modulating the pathogen's colonization, growth, and virulence to impact host survival. Gnotobiotic mice colonized with the amino acid fermenter Paraclostridium bifermentans survive infection with reduced disease severity, while mice colonized with the butyrate-producer, Clostridium sardiniense, succumb more rapidly. Systematic in vivo analyses revealed how each commensal alters the gut-nutrient environment to modulate the pathogen's metabolism, gene regulatory networks, and toxin production. Oral administration of P. bifermentans rescues conventional, clindamycin-treated mice from lethal C. difficile infection in a manner similar to that of monocolonized animals, thereby supporting the therapeutic potential of this commensal species. Our findings lay the foundation for mechanistically informed therapies to counter C. difficile disease using systems biology approaches to define host-commensal-pathogen interactions in vivo.

History and Future Perspectives on the Discipline of Quantitative Systems Pharmacology Modeling and Its Applications.

Mathematical biology and pharmacology models have a long and rich history in the fields of medicine and physiology, impacting our understanding of disease mechanisms and the development of novel therapeutics. With an increased focus on the pharmacology application of system models and the advances in data science spanning mechanistic and empirical approaches, there is a significant opportunity and promise to leverage these advancements to enhance the development and application of the systems pharmacology field. In this paper, we will review milestones in the evolution of mathematical biology and pharmacology models, highlight some of the gaps and challenges in developing and applying systems pharmacology models, and provide a vision for an integrated strategy that leverages advances in adjacent fields to overcome these challenges.

Integrated genetic and metabolic landscapes predict vulnerabilities of temozolomide resistant glioblastoma cells.

Metabolic reprogramming and its molecular underpinnings are critical to unravel the duality of cancer cell function and chemo-resistance. Here, we use a constraints-based integrated approach to delineate the interplay between metabolism and epigenetics, hardwired in the genome, to shape temozolomide (TMZ) resistance. Differential metabolism was identified in response to TMZ at varying concentrations in both the resistant neurospheroidal (NSP) and the susceptible (U87MG) glioblastoma cell-lines. The genetic basis of this metabolic adaptation was characterized by whole exome sequencing that identified mutations in signaling pathway regulators of growth and energy metabolism. Remarkably, our integrated approach identified rewiring in glycolysis, TCA cycle, malate aspartate shunt, and oxidative phosphorylation pathways. The differential killing of TMZ resistant NSP by Rotenone at low concentrations with an IC50 value of 5 nM, three orders of magnitude lower than for U87MG that exhibited an IC50 value of 1.8 mM was thus identified using our integrated systems-based approach.

Transcriptome signature of cell viability predicts drug response and drug interaction in Mycobacterium tuberculosis.

There is an urgent need for new drug regimens to rapidly cure tuberculosis. Here, we report the development of drug response assayer (DRonA) and "MLSynergy," algorithms to perform rapid drug response assays and predict response of Mycobacterium tuberculosis (Mtb) to drug combinations. Using a transcriptome signature for cell viability, DRonA detects Mtb killing by diverse mechanisms in broth culture, macrophage infection, and patient sputum, providing an efficient and more sensitive alternative to time- and resource-intensive bacteriologic assays. Further, MLSynergy builds on DRonA to predict synergistic and antagonistic multidrug combinations using transcriptomes of Mtb treated with single drugs. Together, DRonA and MLSynergy represent a generalizable framework for rapid monitoring of drug effects in host-relevant contexts and accelerate the discovery of efficacious high-order drug combinations.

Integrated constraints based analysis of an engineered violacein pathway in Escherichia coli.

Strategies towards optimal violacein biosynthesis, a potential drug molecule, need systems level coordination of enzymatic activities of individual genes in a multigene operon vioABCDE. Constraints-based flux balance analysis of an extended iAF1260 model (iAF1260vio) with a reconstructed violacein module predicted growth and violacein yields in Escherichia coli accurately. Shadow price (SP) analysis identified tryptophan metabolism and NADPH as limiting. Increased tryptophan levels in Δpgi & ΔpheA were validated using in silico gene deletion analysis. Phenotypic phase plane (PhPP) analysis highlighted sensitivity between tryptophan and NADPH for violacein synthesis at molar growth yields. A synthetic VioABCDE operon (SYNO) sequence was designed to maximize Codon Adaptive Index (CAI: 0.9) and tune translation initiation rates (TIR: 2-50 fold higher) in E. coli. All pSYN E. coli transformants produced higher violacein, with a maximum six-fold increase in yields. The rational design E. coli: ΔpheA SYN: gave the highest violacein titers (33.8 mg/l). Such integrated approaches targeting multiple molecular hierarchies in the cell can be extended further to increase violacein yields.