The EcoCyc database: reflecting new knowledge about Escherichia coli K-12

A single rare σ70 variant establishes a unique gene expression pattern in the E. coli pathobiont LF82

LF82, an adherent-invasive Escherichia coli (AIEC) pathobiont, is associated with Crohn's disease, an inflammatory bowel disease of unknown etiology. Although AIEC phenotypes differ from those of 'commensal' or pathogenic E. coli, work has failed to identify genetic features accounting for these differences. We have investigated a natural, but rare, single nucleotide polymorphism (SNP) in LF82 present within the highly conserved rpoD gene, encoding σ70 [primary sigma factor, RNA polymerase (RNAP)]. We demonstrate that σ70 D445V results in transcriptomic and phenotypic changes consistent with LF82 phenotypes, including increased antibiotic resistance and biofilm formation and increased capacity for methionine biosynthesis. RNA-seq analyses comparing σ70 V445 versus σ70 D445 identified 24 genes upregulated by σ70 V445 in both LF82 and the laboratory E. coli K-12 strain MG1655. Using in vitro transcription, we demonstrate that σ70 D445V directly increases transcription from promoters for several of the up-regulated genes and that the presence of a 16 bp spacer and -14 G:C is associated with this increase. The position of D445V within RNAP suggests that it could affect RNAP/spacer interaction. Our work represents the first identification of a distinguishing SNP for this pathobiont and suggests an underrecognized mechanism by which pathobionts and strain variants can emerge.

Carbapenem-resistant Acinetobacter baumannii (CRAB): metabolic adaptation and transcriptional response to human urine (HU)

Carbapenem-resistant Acinetobacter baumannii (CRAB) is a major human pathogen and a research priority for developing new antimicrobial agents. CRAB is a causative agent of a variety of infections in different body sites. One of the manifestations is catheter-associated urinary tract infection, which exposes the bacteria to the host's urine, creating a particular environment. Exposure of two CRAB clinical isolates, AB5075 and AMA40, to human urine (HU) resulted in the differential expression levels of 264 and 455 genes, respectively, of which 112 were common to both strains. Genes within this group play roles in metabolic pathways such as phenylacetic acid (PAA) catabolism, the Hut system, the tricarboxylic acid (TCA) cycle, and other processes like quorum sensing and biofilm formation. These results indicate that the presence of HU induces numerous adaptive changes in gene expression of the infecting bacteria. These changes presumably help bacteria establish and thrive in the hostile conditions in the urinary tract. These analyses advance our understanding of CRAB's metabolic adaptations to human fluids, as well as expand knowledge on bacterial responses to distinct human fluids containing different concentrations of human serum albumin (HSA).

An efficient in vivo-inducible CRISPR interference system for group A Streptococcus genetic analysis and pathogenesis studies

While genome-wide transposon mutagenesis screens have identified numerous essential genes in the significant human pathogen Streptococcus pyogenes (group A Streptococcus or GAS), many of their functions remain elusive. This knowledge gap is attributed in part to the limited molecular toolbox for controlling GAS gene expression and the bacterium's poor genetic transformability. CRISPR interference (CRISPRi), using catalytically inactive GAS Cas9 (dCas9), is a powerful approach to specifically repress gene expression in both bacteria and eukaryotes, but ironically, it has never been harnessed for controlled gene expression in GAS. In this study, we present a highly transformable and fully virulent serotype M1T1 GAS strain and introduce a doxycycline-inducible CRISPRi system for efficient repression of bacterial gene expression. We demonstrate highly efficient, oligo-based single guide RNA cloning directly to GAS, enabling the construction of a gene knockdown strain in just 2 days, in contrast to the several weeks typically required. The system is shown to be titratable and functional both in vitro and in vivo using a murine model of GAS infection. Furthermore, we provide direct in vivo evidence that the expression of the conserved cell division gene ftsZ is essential for GAS virulence, highlighting its promise as a target for emerging FtsZ inhibitors. Finally, we introduce SpyBrowse (https://veeninglab.com/SpyBrowse), a comprehensive and user-friendly online resource for visually inspecting and exploring GAS genetic features. The tools and methodologies described in this work are poised to facilitate fundamental research in GAS, contribute to vaccine development, and aid in the discovery of antibiotic targets.

IMPORTANCE

While group A Streptococcus (GAS) remains a predominant cause of bacterial infections worldwide, there are limited genetic tools available to study its basic cell biology. Here, we bridge this gap by creating a highly transformable, fully virulent M1T1 GAS strain. In addition, we established a tight and titratable doxycycline-inducible system and developed CRISPR interference (CRISPRi) for controlled gene expression in GAS. We show that CRISPRi is functional in vivo in a mouse infection model. Additionally, we present SpyBrowse, an intuitive and accessible genome browser (https://veeninglab.com/SpyBrowse). Overall, this work overcomes significant technical challenges of working with GAS and, together with SpyBrowse, represents a valuable resource for researchers in the GAS field.

Stress response in Escherichia coli following sublethal phenalene-1-one mediated antimicrobial photodynamic therapy: an RNA-Seq study

Since the molecular mechanisms behind adaptation and the bacterial stress response toward antimicrobial photodynamic therapy (aPDT) are not entirely clear yet, the aim of the present study was to investigate the transcriptomic stress response in Escherichia coli after sublethal treatment with aPDT using RNA sequencing (RNA-Seq). Planktonic cultures of stationary phase E. coli were treated with aPDT using a sublethal dose of the photosensitizer SAPYR. After treatment, RNA was extracted, and RNA-Seq was performed on the Illumina NextSeq 500. Differentially expressed genes were analyzed and validated by qRT-PCR. Furthermore, expression of specific stress response proteins was investigated using Western blot analysis.The analysis of the differential gene expression following pathway enrichment analysis revealed a considerable number of genes and pathways significantly up- or down-regulated in E. coli after sublethal treatment with aPDT. Expression of 1018 genes was up-regulated and of 648 genes was down-regulated after sublethal treatment with aPDT as compared to irradiated controls. Analysis of differentially expressed genes and significantly de-regulated pathways showed regulation of genes involved in oxidative stress response and bacterial membrane damage. In conclusion, the results show a transcriptomic stress response in E. coli upon exposure to aPDT using SAPYR and give an insight into potential molecular mechanisms that may result in development of adaptation.

Host-microbe interactions rewire metabolism in a C. elegans model of leucine breakdown deficiency

In humans, defects in leucine catabolism cause a variety of inborn errors in metabolism. Here, we use Caenorhabditis elegans to investigate the impact of mutations in mccc-1, an enzyme that functions in leucine breakdown. Through untargeted metabolomic and transcriptomic analyses we find extensive metabolic rewiring that helps to detoxify leucine breakdown intermediates via conversion into previously undescribed metabolites and to synthesize mevalonate, an essential metabolite. We also find that the leucine breakdown product 3,3-hydroxymethylbutyrate (HMB), commonly used as a human muscle-building supplement, is toxic to C. elegans and that bacteria modulate this toxicity. Unbiased genetic screens revealed interactions between the host and microbe, where components of bacterial pyrimidine biosynthesis mitigate HMB toxicity. Finally, upregulated ketone body metabolism genes in mccc-1 mutants provide an alternative route for biosynthesis of the mevalonate precursor 3-hydroxy-3-methylglutaryl-CoA. Our work demonstrates that a complex host-bacteria interplay rewires metabolism to allow host survival when leucine catabolism is perturbed.

Global regulation via modulation of ribosome pausing by the ABC-F protein EttA

Having multiple rounds of translation of the same mRNA creates dynamic complexities along with opportunities for regulation related to ribosome pausing and stalling at specific sequences. Yet, mechanisms controlling these critical processes and the principles guiding their evolution remain poorly understood. Through genetic, genomic, physiological, and biochemical approaches, we demonstrate that regulating ribosome pausing at specific amino acid sequences can produce ~2-fold changes in protein expression levels which strongly influence cell growth and therefore evolutionary fitness. We demonstrate, both in vivo and in vitro, that the ABC-F protein EttA directly controls the translation of mRNAs coding for a subset of enzymes in the tricarboxylic acid (TCA) cycle and its glyoxylate shunt, which modulates growth in some chemical environments. EttA also modulates expression of specific proteins involved in metabolically related physiological and stress-response pathways. These regulatory activities are mediated by EttA rescuing ribosomes paused at specific patterns of negatively charged residues within the first 30 amino acids of nascent proteins. We thus establish a unique global regulatory paradigm based on sequence-specific modulation of translational pausing.

FLiPPR: A Processor for Limited Proteolysis (LiP) Mass Spectrometry Data Sets Built on FragPipe

Here, we present FLiPPR, or FragPipe LiP (limited proteolysis) Processor, a tool that facilitates the analysis of data from limited proteolysis mass spectrometry (LiP-MS) experiments following primary search and quantification in FragPipe. LiP-MS has emerged as a method that can provide proteome-wide information on protein structure and has been applied to a range of biological and biophysical questions. Although LiP-MS can be carried out with standard laboratory reagents and mass spectrometers, analyzing the data can be slow and poses unique challenges compared to typical quantitative proteomics workflows. To address this, we leverage FragPipe and then process its output in FLiPPR. FLiPPR formalizes a specific data imputation heuristic that carefully uses missing data in LiP-MS experiments to report on the most significant structural changes. Moreover, FLiPPR introduces a data merging scheme and a protein-centric multiple hypothesis correction scheme, enabling processed LiP-MS data sets to be more robust and less redundant. These improvements strengthen statistical trends when previously published data are reanalyzed with the FragPipe/FLiPPR workflow. We hope that FLiPPR will lower the barrier for more users to adopt LiP-MS, standardize statistical procedures for LiP-MS data analysis, and systematize output to facilitate eventual larger-scale integration of LiP-MS data.

Carbapenem-resistant Acinetobacter baumannii (CRAB): metabolic adaptation and transcriptional response to human urine (HU)

Carbapenem-resistant Acinetobacter baumannii (CRAB) is a major human pathogen and a research priority for developing new antimicrobial agents. CRAB is a causative agent of a variety of infections in different body sites. One of the manifestations is catheter-associated urinary tract infection, which exposes the bacteria to the host's urine, creating a particular environment. Exposure of two CRAB clinical isolates, AB5075 and AMA40, to human urine (HU) resulted in the differential expression levels of 264 and 455 genes, respectively, of which 112 were common to both strains. Genes within this group play roles in metabolic pathways such as phenylacetic acid (PAA) catabolism, the Hut system, the tricarboxylic acid (TCA) cycle, and other processes like quorum sensing and biofilm formation. These results indicate that the presence of HU induces numerous adaptive changes in gene expression of the infecting bacteria. These modifications presumably help bacteria establish and thrive in the hostile conditions in the urinary tract. These analyses advance our understanding of CRAB's metabolic adaptations to human fluids, as well as expanding knowledge on bacterial responses to distinct human fluids containing different concentrations of human serum albumin (HSA).

ZEPPI: Proteome-scale sequence-based evaluation of protein-protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence coevolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI can be implemented on a proteome-wide scale and is applied here to millions of structural models of dimeric complexes in the Escherichia coli and human interactomes found in the PrePPI database. PrePPI's scoring function is based primarily on the evaluation of protein-protein interfaces, and ZEPPI adds a new feature to this analysis through the incorporation of evolutionary information. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. As we discuss, the standard CAPRI scores used to evaluate docking models are based on model quality and not on the ability to give yes/no answers as to whether two proteins interact. ZEPPI is able to detect weak signals from PPI models that the CAPRI scores define as incorrect and, similarly, to identify potential PPIs defined as low confidence by the current PrePPI scoring function. A number of examples that illustrate how the combination of PrePPI and ZEPPI can yield functional hypotheses are provided.

Two synonymous single-nucleotide polymorphisms promoting fluoroquinolone resistance of Escherichia coli in the environment

Single-nucleotide polymorphism (SNP) is one of the core mechanisms that respond to antibiotic resistance of Escherichia coli (E. coli), which is a major issue in environmental pollution. A specific type of SNPs, synonymous SNPs, have been generally considered as the "silent" SNPs since they do not change the encoded amino acid. However, the impact of synonymous SNPs on mRNA splicing, nucleo-cytoplasmic export, stability, and translation was gradually discovered in the last decades. Figuring out the mechanism of synonymous SNPs in regulating antibiotic resistance is critical to improve antimicrobial therapy strategies in clinics and biological treatment strategies of antibiotic-resistant E. coli-polluted materials. With our newly designed antibiotic resistant SNPs prediction algorithm, Multilocus Sequence Type based Identification for Phenotype-single nucleotide polymorphism Analysis (MIPHA), and in vivo validation, we identified 2 important synonymous SNPs 522 G>A and 972 C>T, located at hisD gene, which was previously predicted as a fluoroquinolone resistance-related gene without a detailed mechanism in the E. coli samples with environmental backgrounds. We first discovered that hisD causes gyrA mutation via the upregulation of sbmC and its downstream gene umuD. Moreover, those 2 synonymous SNPs of hisD cause its own translational slowdown and further reduce the expression levels of sbmC and its downstream gene umuD, making the fluoroquinolone resistance determining region of gyrA remains unmutated, ultimately causing the bacteria to lose their ability to resist drugs. This study provided valuable insight into the role of synonymous SNPs in mediating antibiotic resistance of bacteria and a new perspective for the treatment of environmental pollution caused by drug-resistant bacteria.

Combined functional genomic and metabolomic approaches identify new genes required for growth in human urine by multidrug-resistant Escherichia coli ST131

Urinary tract infections (UTIs) are one of the most common bacterial infections in humans, with ~400 million cases across the globe each year. Uropathogenic Escherichia coli (UPEC) is the major cause of UTI and increasingly associated with antibiotic resistance. This scenario has been worsened by the emergence and spread of pandemic UPEC sequence type 131 (ST131), a multidrug-resistant clone associated with extraordinarily high rates of infection. Here, we employed transposon-directed insertion site sequencing in combination with metabolomic profiling to identify genes and biochemical pathways required for growth and survival of the UPEC ST131 reference strain EC958 in human urine (HU). We identified 24 genes required for growth in HU, which mapped to diverse pathways involving small peptide, amino acid and nucleotide metabolism, the stringent response pathway, and lipopolysaccharide biosynthesis. We also discovered a role for UPEC resistance to fluoride during growth in HU, most likely associated with fluoridation of drinking water. Complementary nuclear magnetic resonance (NMR)-based metabolomics identified changes in a range of HU metabolites following UPEC growth, the most pronounced being L-lactate, which was utilized as a carbon source via the L-lactate dehydrogenase LldD. Using a mouse UTI model with mixed competitive infection experiments, we demonstrated a role for nucleotide metabolism and the stringent response in UPEC colonization of the mouse bladder. Together, our application of two omics technologies combined with different infection-relevant settings has uncovered new factors required for UPEC growth in HU, thus enhancing our understanding of this pivotal step in the UPEC infection pathway.

IMPORTANCE

Uropathogenic Escherichia coli (UPEC) cause ~80% of all urinary tract infections (UTIs), with increasing rates of antibiotic resistance presenting an urgent threat to effective treatment. To cause infection, UPEC must grow efficiently in human urine (HU), necessitating a need to understand mechanisms that promote its adaptation and survival in this nutrient-limited environment. Here, we used a combination of functional genomic and metabolomic techniques and identified roles for the metabolism of small peptides, amino acids, nucleotides, and L-lactate, as well as the stringent response pathway, lipopolysaccharide biosynthesis, and fluoride resistance, for UPEC growth in HU. We further demonstrated that pathways involving nucleotide metabolism and the stringent response are required for UPEC colonization of the mouse bladder. The UPEC genes and metabolic pathways identified in this study represent targets for the development of innovative therapeutics to prevent UPEC growth during human UTI, an urgent need given the rapidly rising rates of global antibiotic resistance.

Genome-Wide Screens Identify Genes Responsible for Intrinsic Boric Acid Resistance in Escherichia coli

Boric acid (BA) has antimicrobial properties and is used to combat bacterial infections, including Enterobacteria. However, the molecular mechanisms and cellular responses to BA are still unknown. This genomics study aims to provide new information on the genes and molecular mechanisms related to the antimicrobial effect of BA in Escherichia coli. The Keio collection of E. coli was used to screen 3985 single-gene knockout strains in order to identify mutant strains that were sensitive or hypersensitive to BA at certain concentrations. The mutant strains were exposed to different concentrations of BA ranging from 0 to 120 mM in LB media. Through genome-wide screens, 92 mutants were identified that were relatively sensitive to BA at least at one concentration tested. The related biological processes in the particular cellular system were listed. This study demonstrates that intrinsic BA resistance is the result of various mechanisms acting together. Additionally, we identified eighteen out of ninety-two mutant strains (Delta_aceF, aroK, cheZ, dinJ, galS, garP, glxK, nohA, talB, torR, trmU, trpR, yddE, yfeS, ygaV, ylaC, yoaC, yohN) that exhibited sensitivity using other methods. To increase sensitivity to BA, we constructed double and triple knockout mutants of the selected sensitive mutants. In certain instances, engineered double and triple mutants exhibited significantly amplified effects. Overall, our analysis of these findings offers further understanding of the mechanisms behind BA toxicity and intrinsic resistance in E. coli.

Engineering Escherichia coli for increased Und-P availability leads to material improvements in glycan expression technology

BACKGROUND

Bacterial surface glycans are assembled by glycosyltransferases (GTs) that transfer sugar monomers to long-chained lipid carriers. Most bacteria employ the 55-carbon chain undecaprenyl phosphate (Und-P) to scaffold glycan assembly. The amount of Und-P available for glycan synthesis is thought to be limited by the rate of Und-P synthesis and by competition for Und-P between phosphoglycosyl transferases (PGTs) and GTs that prime glycan assembly (which we collectively refer to as PGT/GTs). While decreasing Und-P availability disrupts glycan synthesis and promotes cell death, less is known about the effects of increased Und-P availability.

RESULTS

To determine if cells can maintain higher Und-P levels, we first reduced intracellular competition for Und-P by deleting all known non-essential PGT/GTs in the Gram-negative bacterium Escherichia coli (hereafter called ΔPGT/GT cells). We then increased the rate of Und-P synthesis in ΔPGT/GT cells by overexpressing the Und-P(P) synthase uppS from a plasmid (puppS). Und-P quantitation revealed that ΔPGT/GT/puppS cells can be induced to maintain 3-fold more Und-P than wild type cells. Next, we determined how increasing Und-P availability affects glycan expression. Interestingly, increasing Und-P availability increased endogenous and recombinant glycan expression. In particular, ΔPGT/GT/puppS cells could be induced to express 7-fold more capsule from Streptococcus pneumoniae serotype 4 than traditional E. coli cells used to express recombinant glycans.

CONCLUSIONS

We demonstrate that the biotechnology standard bacterium E. coli can be engineered to maintain higher levels of Und-P. The results also strongly suggest that Und-P pathways can be engineered to increase the expression of potentially any Und-P-dependent polymer. Given that many bacterial glycans are central to the production of vaccines, diagnostics, and therapeutics, increasing Und-P availability should be a foremost consideration when designing bacterial glycan expression systems.

Creating new-to-nature carbon fixation: A guide

Synthetic biology aims at designing new biological functions from first principles. These new designs allow to expand the natural solution space and overcome the limitations of naturally evolved systems. One example is synthetic CO2-fixation pathways that promise to provide more efficient ways for the capture and conversion of CO2 than natural pathways, such as the Calvin Benson Bassham (CBB) cycle of photosynthesis. In this review, we provide a practical guideline for the design and realization of such new-to-nature CO2-fixation pathways. We introduce the concept of "synthetic CO2-fixation", and give a general overview over the enzymology and topology of synthetic pathways, before we derive general principles for their design from their eight naturally evolved analogs. We provide a comprehensive summary of synthetic carbon-assimilation pathways and derive a step-by-step, practical guide from the theoretical design to their practical implementation, before ending with an outlook on new developments in the field.

Bacterial RNase III: Targets and physiology

Bacteria can rapidly adapt to changes in their environment thanks to the innate flexibility of their genetic expression. The high turnover rate of RNAs, in particular messenger and regulatory RNAs, provides an important contribution to this dynamic adjustment. Recycling of RNAs is ensured by ribonucleases, among which RNase III is the focus of this review. RNase III enzymes are highly conserved from prokaryotes to eukaryotes and have the specific ability to cleave double-stranded RNAs. The role of RNase III in bacterial physiology has remained poorly explored for a long time. However, transcriptomic approaches recently uncovered a large impact of RNase III in gene expression in a wide range of bacteria, generating renewed interest in the physiological role of RNase III. In this review, we first describe the RNase III targets identified from global approaches in 8 bacterial species within 4 Phyla. We then present the conserved and unique functions of bacterial RNase III focusing on growth, resistance to stress, biofilm formation, motility and virulence. Altogether, this review highlights the underestimated impact of RNase III in bacterial adaptation.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

The EcoCyc Database (2023)

EcoCyc is a bioinformatics database available online at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on the regulation of gene expression, E. coli gene essentiality, and nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for the analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed online. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. Data generated from a whole-cell model that is parameterized from the latest data on EcoCyc are also available. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

FLiPPR: A Processor for Limited Proteolysis (LiP) Mass Spectrometry Datasets Built on FragPipe

Here, we present FLiPPR, or FragPipe LiP (limited proteolysis) Processor, a tool that facilitates the analysis of data from limited proteolysis mass spectrometry (LiP-MS) experiments following primary search and quantification in FragPipe. LiP-MS has emerged as a method that can provide proteome-wide information on protein structure and has been applied to a range of biological and biophysical questions. Although LiP-MS can be carried out with standard laboratory reagents and mass spectrometers, analyzing the data can be slow and poses unique challenges compared to typical quantitative proteomics workflows. To address this, we leverage the fast, sensitive, and accurate search and label-free quantification algorithms in FragPipe and then process its output in FLiPPR. FLiPPR formalizes a specific data imputation heuristic that carefully uses missing data in LiP-MS experiments to report on the most significant structural changes. Moreover, FLiPPR introduces a new data merging scheme (from ions to cut-sites) and a protein-centric multiple hypothesis correction scheme, collectively enabling processed LiP-MS datasets to be more robust and less redundant. These improvements substantially strengthen statistical trends when previously published data are reanalyzed with the FragPipe/FLiPPR workflow. As a final feature, FLiPPR facilitates the collection of structural metadata to identify correlations between experiments and structural features. We hope that FLiPPR will lower the barrier for more users to adopt LiP-MS, standardize statistical procedures for LiP-MS data analysis, and systematize output to facilitate eventual larger-scale integration of LiP-MS data.

Multi-scale models of whole cells: progress and challenges

Whole-cell modeling is "the ultimate goal" of computational systems biology and "a grand challenge for 21st century" (Tomita, Trends in Biotechnology, 2001, 19(6), 205-10). These complex, highly detailed models account for the activity of every molecule in a cell and serve as comprehensive knowledgebases for the modeled system. Their scope and utility far surpass those of other systems models. In fact, whole-cell models (WCMs) are an amalgam of several types of "system" models. The models are simulated using a hybrid modeling method where the appropriate mathematical methods for each biological process are used to simulate their behavior. Given the complexity of the models, the process of developing and curating these models is labor-intensive and to date only a handful of these models have been developed. While whole-cell models provide valuable and novel biological insights, and to date have identified some novel biological phenomena, their most important contribution has been to highlight the discrepancy between available data and observations that are used for the parametrization and validation of complex biological models. Another realization has been that current whole-cell modeling simulators are slow and to run models that mimic more complex (e.g., multi-cellular) biosystems, those need to be executed in an accelerated fashion on high-performance computing platforms. In this manuscript, we review the progress of whole-cell modeling to date and discuss some of the ways that they can be improved.

Proteome efficiency of metabolic pathways in Escherichia coli increases along the nutrient flow

IMPORTANCE

Protein translation is the most expensive cellular process in fast-growing bacteria, and efficient proteome usage should thus be under strong natural selection. However, recent studies show that a considerable part of the proteome is unneeded for instantaneous cell growth in Escherichia coli. We still lack a systematic understanding of how this excess proteome is distributed across different pathways as a function of the growth conditions. We estimated the minimal required proteome across growth conditions in E. coli and compared the predictions with experimental data. We found that the proteome allocated to the most expensive internal pathways, including translation and the synthesis of amino acids and cofactors, is near the minimally required levels. In contrast, transporters and central carbon metabolism show much higher proteome levels than the predicted minimal abundance. Our analyses show that the proteome fraction unneeded for instantaneous cell growth decreases along the nutrient flow in E. coli.

ZEPPI: proteome-scale sequence-based evaluation of protein-protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence co-evolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. ZEPPI can be implemented on a proteome-wide scale as evidenced by calculations on millions of structural models of dimeric complexes in the E. coli and human interactomes found in the PrePPI database. A number of examples that illustrate how these tools can yield novel functional hypotheses are provided.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

Progress toward Proteome-Wide Photo-Cross-Linking to Enable Residue-Level Visualization of Protein Structures and Networks In Vivo

Cross-linking mass spectrometry (XL-MS) is emerging as a method at the crossroads of structural and cellular biology, uniquely capable of identifying protein-protein interactions with residue-level resolution and on the proteome-wide scale. With the development of cross-linkers that can form linkages inside cells and easily cleave during fragmentation on the mass spectrometer (MS-cleavable cross-links), it has become increasingly facile to identify contacts between any two proteins in complex samples, including in live cells or tissues. Photo-cross-linkers possess the advantages of high temporal resolution and high reactivity, thereby engaging all residue-types (rather than just lysine); nevertheless, photo-cross-linkers have not enjoyed widespread use and are yet to be employed for proteome-wide studies because their products are challenging to identify. Here, we demonstrate the synthesis and application of two heterobifunctional photo-cross-linkers that feature diazirines and N-hydroxy-succinimidyl carbamate groups, the latter of which unveil doubly fissile MS-cleavable linkages upon acyl transfer to protein targets. Moreover, these cross-linkers demonstrate high water-solubility and cell-permeability. Using these compounds, we demonstrate the feasibility of proteome-wide photo-cross-linking in cellulo. These studies elucidate a small portion of Escherichia coli's interaction network, albeit with residue-level resolution. With further optimization, these methods will enable the detection of protein quinary interaction networks in their native environment at residue-level resolution, and we expect that they will prove useful toward the effort to explore the molecular sociology of the cell.

The glutathione import system satisfies the Staphylococcus aureus nutrient sulfur requirement and promotes interspecies competition

Sulfur is an indispensable element for bacterial proliferation. Prior studies demonstrated that the human pathogen Staphylococcus aureus utilizes glutathione (GSH) as a source of nutrient sulfur; however, mechanisms of GSH acquisition are not defined. Here, we identify a five-gene locus comprising a putative ABC-transporter and predicted γ-glutamyl transpeptidase (ggt) that promotes S. aureus proliferation in medium supplemented with either reduced or oxidized GSH (GSSG) as the sole source of nutrient sulfur. Based on these phenotypes, we name this transporter operon the glutathione import system (gisABCD). Ggt is encoded within the gisBCD operon, and we show that the enzyme is capable of liberating glutamate using either GSH or GSSG as substrates, demonstrating it is a bona fide γ-glutamyl transpeptidase. We also determine that Ggt is expressed in the cytoplasm, representing only the second example of cytoplasmic Ggt localization, the other being Neisseria meningitidis. Bioinformatic analyses revealed that Staphylococcus species closely related to S. aureus encode GisABCD-Ggt homologs. However, homologous systems were not detected in Staphylococcus epidermidis. Consequently, we establish that GisABCD-Ggt provides a competitive advantage for S. aureus over S. epidermidis in a GSH- and GSSG-dependent manner. Overall, this study describes the discovery of a nutrient sulfur acquisition system in S. aureus that targets GSSG in addition to GSH and promotes competition against other staphylococci commonly associated with the human microbiota.

DksA is a conserved master regulator of stress response in Acinetobacter baumannii

Coordination of bacterial stress response mechanisms is critical for long-term survival in harsh environments for successful host infection. The general and specific stress responses of well-studied Gram-negative pathogens like Escherichia coli are controlled by alternative sigma factors, archetypically RpoS. The deadly hospital pathogen Acinetobacter baumannii is notoriously resistant to environmental stresses, yet it lacks RpoS, and the molecular mechanisms driving this incredible stress tolerance remain poorly defined. Here, using functional genomics, we identified the transcriptional regulator DksA as a master regulator for broad stress protection and virulence in A. baumannii. Transcriptomics, phenomics and in vivo animal studies revealed that DksA controls ribosomal protein expression, metabolism, mutation rates, desiccation, antibiotic resistance, and host colonization in a niche-specific manner. Phylogenetically, DksA was highly conserved and well-distributed across Gammaproteobacteria, with 96.6% containing DksA, spanning 88 families. This study lays the groundwork for understanding DksA as a major regulator of general stress response and virulence in this important pathogen.

The Unique N-Terminal Domain of Chlamydial Bactofilin Mediates Its Membrane Localization and Ring-Forming Properties

Chlamydia trachomatis is an obligate intracellular bacterial pathogen. In evolving to the intracellular niche, Chlamydia has reduced its genome size compared to other bacteria and, as a consequence, has a number of unique features. For example, Chlamydia engages the actin-like protein MreB, rather than the tubulin-like protein FtsZ, to direct peptidoglycan (PG) synthesis exclusively at the septum of cells undergoing polarized cell division. Interestingly, Chlamydia possesses another cytoskeletal element-a bactofilin ortholog, BacA. Recently, we reported BacA is a cell size-determining protein that forms dynamic membrane-associated ring structures in Chlamydia that have not been observed in other bacteria with bactofilins. Chlamydial BacA possesses a unique N-terminal domain, and we hypothesized this domain imparts the membrane-binding and ring-forming properties of BacA. We show that different truncations of the N terminus result in distinct phenotypes: removal of the first 50 amino acids (ΔN50) results in large ring structures at the membrane whereas removal of the first 81 amino acids (ΔN81) results in an inability to form filaments and rings and a loss of membrane association. Overexpression of the ΔN50 isoform altered cell size, similar to loss of BacA, suggesting that the dynamic properties of BacA are essential for the regulation of cell size. We further show that the region from amino acid 51 to 81 imparts membrane association as appending it to green fluorescent protein (GFP) resulted in the relocalization of GFP from the cytosol to the membrane. Overall, our findings suggest two important functions for the unique N-terminal domain of BacA and help explain its role as a cell size determinant. IMPORTANCE Bacteria use a variety of filament-forming cytoskeletal proteins to regulate and control various aspects of their physiology. For example, the tubulin-like FtsZ recruits division proteins to the septum whereas the actin-like MreB recruits peptidoglycan (PG) synthases to generate the cell wall in rod-shaped bacteria. Recently, a third class of cytoskeletal protein has been identified in bacteria-bactofilins. These proteins have been primarily linked to spatially localized PG synthesis. Interestingly, Chlamydia, an obligate intracellular bacterium, does not have PG in its cell wall and yet possesses a bactofilin ortholog. In this study, we characterize a unique N-terminal domain of chlamydial bactofilin and show that this domain controls two important functions that affect cell size: its ring-forming and membrane-associating properties.

Revisiting Fur Regulon Leads to a Comprehensive Understanding of Iron and Fur Regulation

Iron is an essential element because it functions as a cofactor of many enzymes, but excess iron causes cell damage. Iron hemostasis in Escherichia coli was transcriptionally maintained by the ferric uptake regulator (Fur). Despite having been studied extensively, the comprehensive physiological roles and mechanisms of Fur-coordinated iron metabolism still remain obscure. In this work, by integrating a high-resolution transcriptomic study of the Fur wild-type and knockout Escherichia coli K-12 strains in the presence or absence of iron with high-throughput ChIP-seq assay and physiological studies, we revisited the regulatory roles of iron and Fur systematically and discovered several intriguing features of Fur regulation. The size of the Fur regulon was expanded greatly, and significant discrepancies were observed to exist between the regulations of Fur on the genes under its direct repression and activation. Fur showed stronger binding strength to the genes under its repression, and genes that were repressed by Fur were more sensitive to Fur and iron regulation as compared to the genes that were activated by Fur. Finally, we found that Fur linked iron metabolism to many essential processes, and the systemic regulations of Fur on carbon metabolism, respiration, and motility were further validated or discussed. These results highlight how Fur and Fur-controlled iron metabolism affect many cellular processes in a systematic way.

The Gene Ontology knowledgebase in 2023

The Gene Ontology (GO) knowledgebase (http://geneontology.org) is a comprehensive resource concerning the functions of genes and gene products (proteins and noncoding RNAs). GO annotations cover genes from organisms across the tree of life as well as viruses, though most gene function knowledge currently derives from experiments carried out in a relatively small number of model organisms. Here, we provide an updated overview of the GO knowledgebase, as well as the efforts of the broad, international consortium of scientists that develops, maintains, and updates the GO knowledgebase. The GO knowledgebase consists of three components: (1) the GO-a computational knowledge structure describing the functional characteristics of genes; (2) GO annotations-evidence-supported statements asserting that a specific gene product has a particular functional characteristic; and (3) GO Causal Activity Models (GO-CAMs)-mechanistic models of molecular "pathways" (GO biological processes) created by linking multiple GO annotations using defined relations. Each of these components is continually expanded, revised, and updated in response to newly published discoveries and receives extensive QA checks, reviews, and user feedback. For each of these components, we provide a description of the current contents, recent developments to keep the knowledgebase up to date with new discoveries, and guidance on how users can best make use of the data that we provide. We conclude with future directions for the project.

Inhibiting fatty acid synthesis overcomes colistin resistance

Treating multidrug-resistant infections has increasingly relied on last-resort antibiotics, including polymyxins, for example colistin. As polymyxins are given routinely, the prevalence of their resistance is on the rise and increases mortality rates of sepsis patients. The global dissemination of plasmid-borne colistin resistance, driven by the emergence of mcr-1, threatens to diminish the therapeutic utility of polymyxins from an already shrinking antibiotic arsenal. Restoring sensitivity to polymyxins using combination therapy with sensitizing drugs is a promising approach to reviving its clinical utility. Here we describe the ability of the biotin biosynthesis inhibitor, MAC13772, to synergize with colistin exclusively against colistin-resistant bacteria. MAC13772 indirectly disrupts fatty acid synthesis (FAS) and restores sensitivity to the last-resort antibiotic, colistin. Accordingly, we found that combinations of colistin and other FAS inhibitors, cerulenin, triclosan and Debio1452-NH₃, had broad potential against both chromosomal and plasmid-mediated colistin resistance in chequerboard and lysis assays. Furthermore, combination therapy with colistin and the clinically relevant FabI inhibitor, Debio1452-NH₃, showed efficacy against mcr-1 positive Klebsiella pneumoniae and colistin-resistant Escherichia coli systemic infections in mice. Using chemical genomics, lipidomics and transcriptomics, we explored the mechanism of the interaction. We propose that inhibiting FAS restores colistin sensitivity by depleting lipid synthesis, leading to changes in phospholipid composition. In all, this work reveals a surprising link between FAS and colistin resistance.

A bacterial genome assembly and annotation laboratory using a virtual machine

With the global increase of infections caused by antibiotic-resistant bacterial strains, there is an urgent need for new methods of tackling the issue. Genomic analysis of bacterial strains can help to understand their virulence and antibiotic resistance profile. Bioinformatic skills are in great demand across the biological sciences. We designed a workshop that allows university students to learn the process of genome assembly using command-line tools within a virtual machine on a Linux operating system. We use Illumina and Nanopore short and long-read raw sequences to reveal the advantages and disadvantages of short, long, and hybrid assembly methods. The workshop teaches how to assess read and assembly quality, perform genome annotation, and analyze pathogenicity, antibiotic and phage resistance. The workshop is intended for a five-week teaching period and is concluded by a student poster presentation assessment.

Quantitative model for genome-wide cyclic AMP receptor protein binding site identification and characteristic analysis

Cyclic AMP receptor proteins (CRPs) are important transcription regulators in many species. The prediction of CRP-binding sites was mainly based on position-weighted matrixes (PWMs). Traditional prediction methods only considered known binding motifs, and their ability to discover inflexible binding patterns was limited. Thus, a novel CRP-binding site prediction model called CRPBSFinder was developed in this research, which combined the hidden Markov model, knowledge-based PWMs and structure-based binding affinity matrixes. We trained this model using validated CRP-binding data from Escherichia coli and evaluated it with computational and experimental methods. The result shows that the model not only can provide higher prediction performance than a classic method but also quantitatively indicates the binding affinity of transcription factor binding sites by prediction scores. The prediction result included not only the most knowns regulated genes but also 1089 novel CRP-regulated genes. The major regulatory roles of CRPs were divided into four classes: carbohydrate metabolism, organic acid metabolism, nitrogen compound metabolism and cellular transport. Several novel functions were also discovered, including heterocycle metabolic and response to stimulus. Based on the functional similarity of homologous CRPs, we applied the model to 35 other species. The prediction tool and the prediction results are online and are available at: https://awi.cuhk.edu.cn/∼CRPBSFinder.

Genome-wide identification of Kanamycin B binding RNA in Escherichia coli

BACKGROUND

The aminoglycosides are established antibiotics that inhibit bacterial protein synthesis by binding to ribosomal RNA. Additional non-antibiotic aminoglycoside cellular functions have also been identified through aminoglycoside interactions with cellular RNAs. The full extent, however, of genome-wide aminoglycoside RNA interactions in Escherichia coli has not been determined. Here, we report genome-wide identification and verification of the aminoglycoside Kanamycin B binding to Escherichia coli RNAs. Immobilized Kanamycin B beads in pull-down assays were used for transcriptome-profiling analysis (RNA-seq).

RESULTS

Over two hundred Kanamycin B binding RNAs were identified. Functional classification analysis of the RNA sequence related genes revealed a wide range of cellular functions. Small RNA fragments (ncRNA, tRNA and rRNA) or small mRNA was used to verify the binding with Kanamycin B in vitro. Kanamycin B and ibsC mRNA was analysed by chemical probing.

CONCLUSIONS

The results will provide biochemical evidence and understanding of potential extra-antibiotic cellular functions of aminoglycosides in Escherichia coli.

Molecular Signatures of the Eagle Effect Induced by the Artificial Siderophore Conjugate LP-600 in E. coli

Achieving cellular uptake is a central challenge for novel antibiotics targeting Gram-negative bacterial pathogens. One strategy is to hijack the bacterial iron transport system by siderophore-antibiotic conjugates that are actively imported into the cell. This was realized with the MECAM-ampicillin conjugate LP-600 we recently reported that was highly active against E. coli. In the present study, we investigate a paradoxical regrowth of E. coli upon treatment of LP-600 at concentrations 16-32 times above the minimum inhibitory concentration (MIC). The phenomenon, coined "Eagle-effect" in other systems, was not due to resistance formation, and it occurred for the siderophore conjugate but not for free ampicillin. To investigate the molecular imprint of the Eagle effect, a combined transcriptome and untargeted metabolome analysis was conducted. LP-600 induced the expression of genes involved in iron acquisition, SOS response, and the e14 prophage upon regrowth conditions. The Eagle effect was diminished in the presence of sulbactam, which we ascribe to a putative synergistic antibiotic action but not to β-lactamase inhibition. The study highlights the relevance of the Eagle effect for siderophore conjugates. Through the first systematic -omics investigations, it also demonstrates that the Eagle effect manifests not only in a paradoxical growth but also in unique gene expression and metabolite profiles.

Genes Vary Greatly in Their Propensity for Collateral Fitness Effects of Mutations

Mutations can have deleterious fitness effects when they decrease protein specific activity or decrease active protein abundance. Mutations will also be deleterious when they cause misfolding or misinteractions that are toxic to the cell (i.e., independent of whether the mutations affect specific activity and abundance). The extent to which protein evolution is shaped by these and other collateral fitness effects is unclear in part because little is known of their frequency and magnitude. Using deep mutational scanning (DMS), we previously found at least 42% of missense mutations in the TEM-1 β-lactamase antibiotic resistance gene cause deleterious collateral fitness effects. Here, we used DMS to comprehensively determine the collateral fitness effects of missense mutations in three genes encoding the antibiotic resistance proteins New Delhi metallo-β-lactamase (NDM-1), chloramphenicol acetyltransferase I (CAT-I), and 2″-aminoglycoside nucleotidyltransferase (AadB). AadB (20%), CAT-I (0.9%), and NDM-1 (0.2%) were less susceptible to deleterious collateral fitness effects than TEM-1 (42%) indicating that genes have different propensities for these effects. As was observed with TEM-1, all the studied deleterious aadB mutants increased aggregation. However, aggregation did not correlate with collateral fitness effects for many of the deleterious mutants of CAT-I and NDM-1. Select deleterious mutants caused unexpected phenotypes to emerge. The introduction of internal start codons in CAT-1 caused loss of the episome and a mutation in aadB made its cognate antibiotic essential for growth. Our study illustrates how the complexity of the cell provides a rich environment for collateral fitness effects and new phenotypes to emerge.

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.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

Evolved bacterial resistance to the chemotherapy gemcitabine modulates its efficacy in co-cultured cancer cells

Drug metabolism by the microbiome can influence anticancer treatment success. We previously suggested that chemotherapies with antimicrobial activity can select for adaptations in bacterial drug metabolism that can inadvertently influence the host's chemoresistance. We demonstrated that evolved resistance against fluoropyrimidine chemotherapy lowered its efficacy in worms feeding on drug-evolved bacteria (Rosener et al., 2020). Here, we examine a model system that captures local interactions that can occur in the tumor microenvironment. Gammaproteobacteria-colonizing pancreatic tumors can degrade the nucleoside-analog chemotherapy gemcitabine and, in doing so, can increase the tumor's chemoresistance. Using a genetic screen in Escherichia coli, we mapped all loss-of-function mutations conferring gemcitabine resistance. Surprisingly, we infer that one third of top resistance mutations increase or decrease bacterial drug breakdown and therefore can either lower or raise the gemcitabine load in the local environment. Experiments in three E. coli strains revealed that evolved adaptation converged to inactivation of the nucleoside permease NupC, an adaptation that increased the drug burden on co-cultured cancer cells. The two studies provide complementary insights on the potential impact of microbiome adaptation to chemotherapy by showing that bacteria-drug interactions can have local and systemic influence on drug activity.

Enzyme expression kinetics by Escherichia coli during transition from rich to minimal media depends on proteome reserves

Bacterial fitness depends on adaptability to changing environments. In rich growth medium, which is replete with amino acids, Escherichia coli primarily expresses protein synthesis machineries, which comprise ~40% of cellular proteins and are required for rapid growth. Upon transition to minimal medium, which lacks amino acids, biosynthetic enzymes are synthesized, eventually reaching ~15% of cellular proteins when growth fully resumes. We applied quantitative proteomics to analyse the timing of enzyme expression during such transitions, and established a simple positive relation between the onset time of enzyme synthesis and the fractional enzyme 'reserve' maintained by E. coli while growing in rich media. We devised and validated a coarse-grained kinetic model that quantitatively captures the enzyme recovery kinetics in different pathways, solely on the basis of proteomes immediately preceding the transition and well after its completion. Our model enables us to infer regulatory strategies underlying the 'as-needed' gene expression programme adopted by E. coli.

Promoter–motif extraction from co-regulated genes and their relevance to co-expression using E. coli as a model

Gene expression varies due to the intrinsic stochasticity of transcription or as a reaction to external perturbations that generate cellular mutations. Co-regulation, co-expression and functional similarity of substances have been employed for indoctrinating the process of the transcriptional paradigm. The difficult process of analysing complicated proteomes and biological switches has been made easier by technical improvements, and microarray technology has flourished as a viable platform. Therefore, this research enables Microarray to cluster genes that are co-expressed and co-regulated into specific segments. Copious search algorithms have been employed to ascertain diacritic motifs or a combination of motifs that are performing regular expression, and their relevant information corresponding to the gene patterns is also documented. The associated genes co-expression and relevant cis-elements are further explored by engaging Escherichia coli as a model organism. Various clustering algorithms have also been used to generate classes of genes with similar expression profiles. A promoter database ‘EcoPromDB’ has been developed by referring RegulonDB database; this promoter database is freely available at www.ecopromdb.eminentbio.com and is divided into two sub-groups, depending upon the results of co-expression and co-regulation analyses.

Nutrient stress is a target for new antibiotics

Novel approaches are required to address the looming threat of pan-resistant Gram-negative pathogens and forestall the rise of untreatable infections. Unconventional targets that are uniquely important during infection and tractable to high-throughput drug discovery methods hold high potential for innovation in antibiotic discovery programs. In this context, inhibitors of bacterial nutrient stress are particularly exciting candidates for future antibiotic development. Amino acid, nucleotide, and vitamin biosynthesis pathways are critical for bacterial growth in nutrient-limiting conditions in the laboratory and the host. Although historically dismissed as dispensable for pathogens, a wealth of transposon mutagenesis and single-mutant studies have emerged which demonstrate that several such pathways are critical for infection. Indeed, high-throughput screens of diverse synthetic compounds and natural products have uncovered inhibitors of nutrient biosynthesis. Herein, we review bacterial nutrient biosynthesis and its role during host infection. Further, we explore screening platforms developed to search for inhibitors of these targets and highlight successes among these. Finally, we feature important and sometimes surprising connections between bacterial nutrient biosynthesis, antibiotic activity, and antibiotic resistance.

Creating De Novo Overlapped Genes

Future applications of synthetic biology will rely on deploying engineered cells outside of lab environments for long periods of time. Currently, a significant roadblock to this application is the potential for deactivating mutations in engineered genes. A recently developed method to protect engineered coding sequences from mutation is called Constraining Adaptive Mutations using Engineered Overlapping Sequences (CAMEOS). In this chapter we provide a workflow for utilizing CAMEOS to create synthetic overlaps between two genes, one essential (infA) and one non-essential (aroB), to protect the non-essential gene from mutation and loss of protein function. In this workflow we detail the methods to collect large numbers of related protein sequences, produce multiple sequence alignments (MSAs), use the MSAs to generate hidden Markov models and Markov random field models, and finally generate a library of overlapping coding sequences through CAMEOS scripts. To assist practitioners with basic coding skills to try out the CAMEOS method, we have created a virtual machine containing all the required packages already installed that can be downloaded and run locally.

Predicting stress response and improved protein overproduction in Bacillus subtilis

Bacillus subtilis is a well-characterized microorganism and a model for the study of Gram-positive bacteria. The bacterium can produce proteins at high densities and yields, which has made it valuable for industrial bioproduction. Like other cell factories, metabolic modeling of B. subtilis has discovered ways to optimize its metabolism toward various applications. The first genome-scale metabolic model (M-model) of B. subtilis was published more than a decade ago and has been applied extensively to understand metabolism, to predict growth phenotypes, and served as a template to reconstruct models for other Gram-positive bacteria. However, M-models are ill-suited to simulate the production and secretion of proteins as well as their proteomic response to stress. Thus, a new generation of metabolic models, known as metabolism and gene expression models (ME-models), has been initiated. Here, we describe the reconstruction and validation of a ME model of B. subtilis, iJT964-ME. This model achieved higher performance scores on the prediction of gene essentiality as compared to the M-model. We successfully validated the model by integrating physiological and omics data associated with gene expression responses to ethanol and salt stress. The model further identified the mechanism by which tryptophan synthesis is upregulated under ethanol stress. Further, we employed iJT964-ME to predict amylase production rates under two different growth conditions. We analyzed these flux distributions and identified key metabolic pathways that permitted the increase in amylase production. Models like iJT964-ME enable the study of proteomic response to stress and the illustrate the potential for optimizing protein production in bacteria.

Identification of Attenuators of Transcriptional Termination: Implications for RNA Regulation in Escherichia coli

The regulatory function of many bacterial small RNAs (sRNAs) requires the binding of the RNA chaperone Hfq to the 3' portion of the sRNA intrinsic terminator, and therefore sRNA signaling might be regulated by modulating its terminator. Here, using a multicopy screen developed with the terminator of sRNA SgrS, we identified an sRNA gene (cyaR) and three protein-coding genes (cspD, ygjH, and rof) that attenuate SgrS termination in Escherichia coli. Analyses of CyaR and YgjH, a putative tRNA binding protein, suggested that the CyaR activity was indirect and the effect of YgjH was moderate. Overproduction of the protein attenuators CspD and Rof resulted in more frequent readthrough at terminators of SgrS and two other sRNAs, and regulation by SgrS of target mRNAs was reduced. The effect of Rof, a known inhibitor of Rho, was mimicked by bicyclomycin or by a rho mutant, suggesting an unexpected role for Rho in sRNA termination. CspD, a member of the cold shock protein family, bound both terminated and readthrough transcripts, stabilizing them and attenuating termination. By RNA sequencing analysis of the CspD overexpression strain, we found global effects of CspD on gene expression across some termination sites. We further demonstrated effects of endogenous CspD under slow growth conditions where cspD is highly expressed. These findings provided evidence of changes in the efficiency of intrinsic termination, confirming this as an additional layer of the regulation of sRNA signaling. IMPORTANCE Growing evidence suggests that the modulation of intrinsic termination and readthrough of transcription is more widespread than previously appreciated. For small RNAs, proper termination plays a critical role in their regulatory function. Here, we present a multicopy screen approach to identify factors that attenuate small RNA termination and therefore abrogate signaling dependent on the small RNA. This study highlights a new aspect of regulation of small RNA signaling as well as the modulation of intrinsic termination.

Dynamic feedback regulation for efficient membrane protein production using a small RNA-based genetic circuit in Escherichia coli

BACKGROUND

Membrane proteins (MPs) are an important class of molecules with a wide array of cellular functions and are part of many metabolic pathways. Despite their great potential-as therapeutic drug targets or in microbial cell factory optimization-many challenges remain for efficient and functional expression in a host such as Escherichia coli.

RESULTS

A dynamically regulated small RNA-based circuit was developed to counter membrane stress caused by overexpression of different MPs. The best performing small RNAs were able to enhance the maximum specific growth rate with 123%. On culture level, the total MP production was increased two-to three-fold compared to a system without dynamic control. This strategy not only improved cell growth and production of the studied MPs, it also suggested the potential use for countering metabolic burden in general.

CONCLUSIONS

A dynamically regulated feedback circuit was developed that can sense metabolic stress caused by, in casu, the overexpression of an MP and responds to it by balancing the metabolic state of the cell and more specifically by downregulating the expression of the MP of interest. This negative feedback mechanism was established by implementing and optimizing simple-to-use genetic control elements based on post-transcriptional regulation: small non-coding RNAs. In addition to membrane-related stress when the MP accumulated in the cytoplasm as aggregates, the sRNA-based feedback control system was still effective for improving cell growth but resulted in a decreased total protein production. This result suggests promiscuity of the MP sensor for more than solely membrane stress.

Toxic effect and inability of L-homoserine to be a nitrogen source for growth of Escherichia coli resolved by a combination of in vivo evolution engineering and omics analyses

L-homoserine is a pivotal intermediate in the carbon and nitrogen metabolism of E. coli. However, this non-canonical amino acid cannot be used as a nitrogen source for growth. Furthermore, growth of this bacterium in a synthetic media is potently inhibited by L-homoserine. To understand this dual effect, an adapted laboratory evolution (ALE) was applied, which allowed the isolation of a strain able to grow with L-homoserine as the nitrogen source and was, at the same time, desensitized to growth inhibition by this amino acid. Sequencing of this evolved strain identified only four genomic modifications, including a 49 bp truncation starting from the stop codon of thrL. This mutation resulted in a modified thrL locus carrying a thrL* allele encoding a polypeptide 9 amino acids longer than the thrL encoded leader peptide. Remarkably, the replacement of thrL with thrL* in the original strain MG1655 alleviated L-homoserine inhibition to the same extent as strain 4E, but did not allow growth with this amino acid as a nitrogen source. The loss of L-homoserine toxic effect could be explained by the rapid conversion of L-homoserine into threonine via the thrL*-dependent transcriptional activation of the threonine operon thrABC. On the other hand, the growth of E. coli on a mineral medium with L-homoserine required an activation of the threonine degradation pathway II and glycine cleavage system, resulting in the release of ammonium ions that were likely recaptured by NAD(P)-dependent glutamate dehydrogenase. To infer about the direct molecular targets of L-homoserine toxicity, a transcriptomic analysis of wild-type MG1655 in the presence of 10 mM L-homoserine was performed, which notably identified a potent repression of locomotion-motility-chemotaxis process and of branched-chain amino acids synthesis. Since the magnitude of these effects was lower in a ΔthrL mutant, concomitant with a twofold lower sensitivity of this mutant to L-homoserine, it could be argued that growth inhibition by L-homoserine is due to the repression of these biological processes. In addition, L-homoserine induced a strong upregulation of genes in the sulfate reductive assimilation pathway, including those encoding its transport. How this non-canonical amino acid triggers these transcriptomic changes is discussed.

ChEAP: ChIP-exo analysis pipeline and the investigation of Escherichia coli RpoN protein-DNA interactions

Genome-scale studies of the bacterial regulatory network have been leveraged by declining sequencing cost and advances in ChIP (chromatin immunoprecipitation) methods. Of which, ChIP-exo has proven competent with its near-single base-pair resolution. While several algorithms and programs have been developed for different analytical steps in ChIP-exo data processing, there is a lack of effort in incorporating them into a convenient bioinformatics pipeline that is intuitive and publicly available. In this paper, we developed ChIP-exo Analysis Pipeline (ChEAP) that executes the one-step process, starting from trimming and aligning raw sequencing reads to visualization of ChIP-exo results. The pipeline was implemented on the interactive web-based Python development environment - Jupyter Notebook, which is compatible with the Google Colab cloud platform to facilitate the sharing of codes and collaboration among researchers. Additionally, users could exploit the free GPU and CPU resources allocated by Colab to carry out computing tasks regardless of the performance of their local machines. The utility of ChEAP was demonstrated with the ChIP-exo datasets of RpoN sigma factor in E. coli K-12 MG1655. To analyze two raw data files, ChEAP runtime was 2 min and 25 s. Subsequent analyses identified 113 RpoN binding sites showing a conserved RpoN binding pattern in the motif search. ChEAP application in ChIP-exo data analysis is extensive and flexible for the parallel processing of data from various organisms.

MlaC belongs to a unique class of non-canonical substrate-binding proteins and follows a novel phospholipid-binding mechanism

The outer membrane (OM) of Gram-negative bacteria acts as a formidable barrier against a plethora of detrimental compounds owing to its asymmetric nature. This is because the OM possesses lipopolysaccharides (LPSs) in the outer leaflet and phospholipids (PLs) in the inner leaflet. The maintenance of lipid asymmetry (Mla) system is involved in preserving the distribution of PLs in OM. The periplasmic component of the system MlaC serves as the substrate-binding protein (SBP) that shuttles PLs between the inner and outer membranes. However, an in-depth report highlighting its mechanism of ligand binding is still lacking. This study reports the crystal structure of MlaC from Escherichia coli (EcMlaC) at a resolution of 2.5 Å in a quasi-open state, complexed with PL. The structural analysis reveals that EcMlaC and orthologs comprise two major domains, viz. nuclear transport factor 2-like (NTF2-like) and phospholipid-binding protein (PBP). Each domain can be further divided into two subdomains arranged in a discontinuous fashion. This study further reveals that EcMlaC is polyspecific in nature and follows a reverse mechanism of the opening of the substrate-binding site during the ligand binding. Furthermore, MlaC can bind two PLs by forming subsites in the binding pocket. These findings, altogether, have led to the proposition of the unique "segmented domain movement" mechanism of PL binding, not reported for any known SBP to date. Further, unlike typical SBPs, MlaC has originated from a cystatin-like fold. Overall, this study establishes MlaC to be a non-canonical SBP with a unique ligand-binding mechanism.

Iron availability enhances the cellular energetics of aerobic Escherichia coli cultures while upregulating anaerobic respiratory chains

Aerobic Escherichia coli growth at restricted iron concentrations (≤ 1.75 ± 0.04 μM) is characterized by lower biomass yield, higher acetate accumulation and higher activation of the siderophore iron-acquisition systems. Although iron homeostasis in E. coli has been studied intensively, previous studies focused only on understanding the regulation of the iron import systems and the iron-requiring enzymes. Here, the effect of iron availability on the energy metabolism of E. coli has been investigated. It was established that aerobic cultures growing under limiting iron conditions showed lower ATP yield per glucose, lower growth rate and lower TCA cycle activity and respiration, at the same time as increased glucose consumption, acetate and pyruvate accumulation, practically mimicking microaerobic growth. However, at excess iron, independent of oxygen availability, the cultures showed high cellular energetics (5.8 ATP/mol of glucose) by using pathways requiring iron-rich complex proteins found in the TCA cycle and respiratory chain. In conditions of iron excess, some iron-requiring terminal reductases of the respiratory chain, that were thought to function only under anaerobiosis, were used by the E. coli, when in aerobic conditions, to maintain high respiratory activity. This allowed it to produce more biomass and more reactive oxygen species that were controlled by the higher activity of the antioxidant defenses (SOD, peroxidase and catalase) and the iron-sulfur cluster repair systems.

Application of comparative genomics of Acetobacter species facilitates genome-scale metabolic reconstruction of the Acetobacter ghanensis LMG 23848T and Acetobacter senegalensis 108B cocoa strains

Acetobacter species play an import role during cocoa fermentation. However, Acetobacter ghanensis and Acetobacter senegalensis are outcompeted during fermentation of the cocoa pulp-bean mass, whereas Acetobacter pasteurianus prevails. In this paper, an in silico approach aimed at delivering some insights into the possible metabolic adaptations of A. ghanensis LMG 23848^T and A. senegalensis 108B, two candidate starter culture strains for cocoa fermentation processes, by reconstructing genome-scale metabolic models (GEMs). Therefore, genome sequence data of a selection of strains of Acetobacter species were used to perform a comparative genomic analysis. Combining the predicted orthologous groups of protein-encoding genes from the Acetobacter genomes with gene-reaction rules of GEMs from two reference bacteria, namely a previously manually curated model of A. pasteurianus 386B (iAp386B454) and two manually curated models of Escherichia coli (EcoCyc and iJO1366), allowed to predict the set of reactions present in A. ghanensis LMG 23848^T and A. senegalensis 108B. The predicted metabolic network was manually curated using genome re-annotation data, followed by the reconstruction of species-specific GEMs. This approach additionally revealed possible differences concerning the carbon core metabolism and redox metabolism among Acetobacter species, pointing to a hitherto unexplored metabolic diversity. More specifically, the presence or absence of reactions related to citrate catabolism and the glyoxylate cycle for assimilation of C2 compounds provided not only new insights into cocoa fermentation but also interesting guidelines for future research. In general, the A. ghanensis LMG 23848^T and A. senegalensis 108B GEMs, reconstructed in a semi-automated way, provided a proof-of-concept toward accelerated formation of GEMs of candidate functional starter cultures for food fermentation processes.