- 1Bioinformatics Research Group, Artificial Intelligence Center, SRI International, Menlo Park, CA, United States
- 2Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México, Cuernavaca, México
- 3Department of Molecular Sciences, Macquarie University, Sydney, NSW, Australia
- 4Instituto de Energías Renovables, Universidad Nacional Autónoma de México, Temixco, México
- 5Department of Microbiology and Immunology, Stritch School of Medicine, Loyola University Chicago, Maywood, IL, United States
- 6Department of Biomedical Engineering, Boston University, Boston, MA, United States
The EcoCyc model-organism database collects and summarizes experimental data for Escherichia coli K-12. EcoCyc is regularly updated by the manual curation of individual database entries, such as genes, proteins, and metabolic pathways, and by the programmatic addition of results from select high-throughput analyses. Updates to the Pathway Tools software that supports EcoCyc and to the web interface that enables user access have continuously improved its usability and expanded its functionality. This article highlights recent improvements to the curated data in the areas of metabolism, transport, DNA repair, and regulation of gene expression. New and revised data analysis and visualization tools include an interactive metabolic network explorer, a circular genome viewer, and various improvements to the speed and usability of existing tools.
Introduction
Escherichia coli is the most well-studied bacterial model organism. The scientific literature reports on more than a century of research on E. coli, including paradigm-shifting research on enzyme function, gene regulation and genetic engineering. Knowledge gained about the biology of E. coli is often the basis for assigning gene product functions in less studied organisms, and scientists turn to the body of E. coli research to begin to understand these functions in the context of their organism of interest. However, despite the long history of research, the functions of a surprising number of E. coli gene products remain unknown (Ghatak et al., 2019). Knowledge gaps remain even in areas that have been studied for decades, and the genes of unknown function that are essential for growth in rich media exist.
The EcoCyc database has been manually curated by PhD-level scientists for nearly three decades (Karp and Riley, 1993; Keseler et al., 2017), and its coverage has been expanded from metabolism to the entire genome. Extensive literature searches enable curators to capture both established knowledge and new insights. Perhaps equally important, the curation process can capture a lack of knowledge via the assignment of detailed evidence codes. For example, the participation of an enzyme in a metabolic pathway is often established by assaying its biochemical function in vitro, resulting in an IDA (inferred from direct assay) evidence code. Occasionally, an enzyme’s function within a metabolic pathway is known only by its mutant phenotype, resulting in an IMP (inferred from mutant phenotype) evidence code. Therefore, EcoCyc provides an overview of current knowledge and serves as a resource for the identification of knowledge gaps.
EcoCyc collects research conducted with the laboratory workhorse K-12 strains projected on the genome sequence of the first sequenced E. coli K-12 strain, MG1655. Many other E. coli strains have been sequenced since that first genome sequence. To leverage the EcoCyc curation effort and enhance the quality and usability of all E. coli databases within the BioCyc database collection (of which EcoCyc is a member database), curated gene and protein data have also been propagated from EcoCyc to orthologs in databases for 480 other E. coli strains via a new automated method (Paley et al., 2021). In this article, we highlight and summarize additions to the data content and improvements to search, data-analysis, and visualization tools since our last publication reporting on updates to EcoCyc (Keseler et al., 2017).
Results
Curated Data in EcoCyc
An overview of many of the data types captured in EcoCyc version 24.5, released on January 7, 2021, is shown in Table 1. This section highlights some notable updates since release version 21.1 (Keseler et al., 2017).
Metabolism
EcoCyc integrates historical data with the most recent insights from the published literature. For example, the enzymes involved in the biosynthesis of ubiquinol-8 were genetically identified decades ago. The current representation of this pathway in EcoCyc can be seen by following this link: https://ecocyc.org/ECOLI/NEW-IMAGE?type=PATHWAY&object=PWY-6708&detail-level=2.
For most of the enzymes, curators were unable to find the published reports of biochemical assays of the activities of ubiquinol-8 biosynthesis enzymes, which is likely due to the general difficulty of, lack of interest in, and/or obstacles to publishing negative data. The unavailability of this information highlights the importance of recording the lack of specific types of data, as is being done in EcoCyc: the evidence codes associated with many of the individual enzymatic reactions in this pathway remain at the “inferred by mutant phenotype” level.
This lack of biochemical data seemed surprising, because most sof the enzymes in ubiquinol-8 biosynthesis, like those in menaquinol-8 biosynthesis,1 are located in the cytoplasm. However, unlike menaquinol-8 biosynthesis, where the hydrophobic octaprenyl tail is added late in the pathway by the inner membrane-localized enzyme MenA, mutant phenotype data showed that the octaprenyl tail of ubiquinol-8 is added early in the pathway. Also, two accessory factors with no predicted biochemical function, UbiJ and UbiK, were identified only by their mutant phenotypes (Aussel et al., 2014; Agrawal et al., 2017; Loiseau et al., 2017). The puzzle pieces fell into place in 2019, when Hajj Chehade et al. discovered that most of the ubiquinol-8 biosynthetic enzymes and the two accessory factors form a soluble complex (metabolon) in the cytoplasm. This complex is able to perform the biochemical transformations while shielding the octaprenyl tail from the aqueous environment (Hajj Chehade et al., 2019). However, other questions remain. The UbiB protein is involved in ubiquinol-8 biosynthesis based on a mutant phenotype. It was originally thought to provide a catalytic activity within the pathway (Cox et al., 1969), but is now proposed to function as a regulator (Poon et al., 2000; Hajj Chehade et al., 2013). Each of these pieces of data can be accessed in multiple ways, for example, by hovering over enzyme names to show the evidence codes associated with their functions and by reading the free-text summaries for the pathway and each enzyme.
Transmembrane Transport
Newly characterized transporters reported in the literature remain a focus for curation. Recent highlights include the curation of the pyruvate:proton symporters BtsT (Kristoficova et al., 2018) and CstA (Hwang et al., 2018; Gasperotti et al., 2020), the Zn2+:proton symporter ZntA (Gati et al., 2017), and a guanidinium:proton antiporter Gdx (Kermani et al., 2018). The latter transporter is regulated by a guanidine-II riboswitch predicted to act as a translation “on” switch (Huang et al., 2017; Sherlock et al., 2017). As part of the curation process, the gene names and free-text summaries for these proteins were updated, and transport reactions (Figure 1A) and regulatory information (Figure 1B) were added.
Figure 1. Curation of Gdx in EcoCyc. (A) The guanidinium:proton antiport reaction mediated by Gdx. (B) Gdx regulation by a guanidine II riboswitch.
The guanidinium:proton antiporter Gdx is a member of the small multidrug resistance (SMR) family of proton-dependent drug efflux transporters. EcoCyc currently represents 25 known energy-dependent drug efflux transporters, including representatives from five of the seven major families of efflux transporters (Chitsaz and Brown, 2017). We have reviewed and updated the curation of all the drug efflux transporters in EcoCyc and improved our representation of the specific substrates, both physiological and non-physiological, that are exported by these proteins. Many new reactions and compounds have been added to the database as a result of this update. Readers interested in this area can view a freely available SmartTable of all drug efflux transporters and their reactions at the following link: https://ecocyc.org/group?id=biocyc14-4655-3823813233.
DNA Repair
Significant improvements have been made to the curation of DNA repair enzymes, with a particular focus on the addition of reactions that accurately reflect the catalytic activities of these important proteins. Eleven new reactions were created as part of this process, including those for two newly described enzymes: the genome maintenance protein encoded by yedK (Mohni et al., 2019; Thompson et al., 2019; Wang et al., 2019) and an interstrand DNA crosslink repair glycosylase encoded by ycaQ (Bradley et al., 2020). Figure 2 shows the new reactions assigned to YedK and YcaQ.
Figure 2. Reactions assigned to recently characterized DNA repair proteins in EcoCyc. (A) Genome maintenance protein YedK. (B) Interstrand DNA crosslink repair glycosylase YcaQ.
Lysine Acetylation Sites
Protein Nε-lysine acetylation is a common post-translational modification, resulting from transfer of an acetyl group (CH3CO) to the ε-amino group (N-ε) of lysine residues within a protein. Acetylation increases the side-chain size and neutralizes the positive charge of the lysine residue, potentially altering protein activity (Christensen et al., 2019). Some proteins regulated by Nε-lysine acetylation include the central metabolic enzymes acetyl-CoA synthetase (Starai and Escalante-Semerena, 2004), enolase (Nakayasu et al., 2017), and malate dehydrogenase (Venkat et al., 2017), as well as the transcription factors PhoP (Ren et al., 2019) and CRP (Davis et al., 2018). Nε-lysine acetylation can be catalyzed by lysine acetyltransferases (KATs) using acetyl-CoA as the acetyl donor. The best studied KAT in E. coli is YfiQ (also known as Pat, PatZ, and Pka). Recently, four novel KATs – YjaB, YiaC, RimI, and PhnO – were revealed (Christensen et al., 2018). Nε-lysine acetylation can also occur without the help of a dedicated enzyme; in this case, the acetyl donor is acetyl phosphate, a high energy central metabolic intermediate that accumulates when carbon is in excess (Weinert et al., 2013; Kuhn et al., 2014; Christensen et al., 2017).
We greatly expanded the coverage of lysine acetylation in EcoCyc by importing five acetylome datasets that identify specific lysine positions in proteins that have been subject to acetylation (Kuhn et al., 2014; Schilling et al., 2015; Christensen et al., 2018). The lysine acetylation sites are recorded and displayed as protein features. When visiting a protein page, clicking on the tab “Protein Features” will show the amino acid sequence and a table of annotations that indicate specific sites or regions with evidence for a variety of functional properties including known acetylation sites. Two examples can be found by following these links for proteins AceF and LipA, respectively: https://ecocyc.org/gene?orgid=ECOLI&id=EG10025#tab=FTRS and https://ecocyc.org/gene?orgid=ECOLI&id=EG11306#tab=FTRS.
In summary, 914 proteins were updated by data showing at least one lysine that can be acetylated. Acetylation data were added to 2,065 distinct lysine residues in the proteome.
The preceding protein pages for AceF and LipA illustrate the ability of EcoCyc to capture the functions of substitution mutants in the Protein Features tab. For example, the page for AceF captures the fact that an H to C substitution at position 603 abolishes the catalytic activity of the protein (see the first feature table). A total of 6,792 such “mutagenesis variant” protein features are present in EcoCyc, although there must be additional such information in the experimental literature. EcoCyc contains 40,051 protein features in total (including the preceding 6,792), including, for example, enzyme active sites and metal ion binding sites.
Regulation of Gene Expression
Since 2017, a significant amount of new data related to specific promoters, regulatory interactions (RIs) and transcription units in E. coli K-12 has been published. This increase is reflected in new database objects and in modifications to existing objects as shown in Table 2. The largest number of modifications comes from enriching summaries and adding new evidence to existing objects.
Table 2. Summary of curation of the regulation of gene expression in EcoCyc between releases 21.1 and 24.5.
We have continued expanding the description of transcriptional regulation by including the binding of regulatory molecules directly to RNA polymerase. Examples are the allosteric regulation of RNA polymerase by ppGpp and DksA.
Regulatory Interactions Extracted From High-Throughput Experiments
As a result of the increasing E. coli K-12 literature involving the use of high-throughput technologies (HTs; Santos-Zavaleta et al., 2018), we have increased the number of DNA binding sites and their associated RIs (Table 2). Of the total number of new RIs, over 1,000 come from HT experiments with seven transcription factors. These RIs were identified by the authors through the combination of genome binding and expression profiling experiments, such as variants of chromatin immunoprecipitation (ChIP) and RNA-seq and microarray analyses, respectively (Table 3).
Table 3. Transcription factors and their regulatory interactions (RIs) extracted from high-throughput experiments.
Redefinition of Basic Concepts in Gene Regulation
The conceptual data model used in EcoCyc to organize the knowledge about transcriptional regulation derives from the initial model by Jacob and Monod of the operon concept (Jacob and Monod, 1961). After 60 years of research with many technological advances before and after the explosion of HT methodologies in genomics, it was the time to revise the classic definitions to update them with our current knowledge on the regulation of transcription initiation in bacteria. Based on the consensus view of a group of experts (Mejía-Almonte et al., 2020), we have modified some aspects of modeling this knowledge in EcoCyc. For instance, a single promoter object was previously used to represent transcription start sites (TSSs) for RNA polymerase holoenzymes containing different sigma factors. Now, each of those TSSs belongs to a different promoter because each may be subject to different regulation even if the TSS is at exactly the same genome location (Mejía-Almonte et al., 2020). Conversely, given the known flexibility of RNA polymerase, one promoter may have more than one TSS within a region of five base pairs (Liu and Turnbough, 1994; Walker and Osuna, 2002; Winkelman et al., 2016). This limit is now being used in EcoCyc to add newly identified TSSs to known promoters. In particular, this is the case with experiments identifying TSSs and their associated transcription units from HT experiments (Yan et al., 2018; Ju et al., 2019).
The Escherichia coli K-12 MG1655 GenBank File, U00096.3
EcoCyc has worked together with the original submitter, Dr. Guy Plunkett III, and staff from UniProt and NCBI to update the E. coli GenBank entry U00096.3, with the last update deposited on September 23, 2020. All genome annotation data within this entry, such as gene symbols, gene positions, and updated function names, are drawn directly from EcoCyc. Gene names are updated from the originally assigned “y-names” if a new name was assigned in the experimental literature. We encourage renaming “y-genes” with Demerec-style gene names (Demerec et al., 1966) once a function has been discovered. A brief summary on the history of the sequenced genome and guidelines for new gene names can be accessed on the following website: https://www.genome.wisc.edu/sequencing/k12.htm.
New Tools in EcoCyc
Metabolic Network Explorer
The Metabolic Network Explorer (see website command Tools → Metabolism → Metabolic Network Explorer) is a new tool for interactively exploring the E. coli metabolic network around a metabolite of interest, as shown in Figure 3. The user specifies a starting metabolite, and the software displays that metabolite along with a full list of potential precursor and successor metabolites derived from the complete reaction network in EcoCyc. The tooltip for each potential precursor or successor metabolite lists all the reactions and enzymes that carry out the transformation and any pathways they belong to. After the user selects a precursor or successor metabolite to add it and its connecting reaction to a central path, that metabolite’s potential precursor and successor metabolites are added to the display. The user can continue to expand the central path in either or both directions by selecting metabolites at the start or end or the user can change the central path by selecting metabolites connected to internal metabolites. A list of paths previously generated in the current session is maintained to allow the user to quickly switch among them. The display includes several customization options such as whether to show metabolite structures or pathway names.
Figure 3. Example display of the Metabolic Network Explorer to explore Escherichia coli metabolism starting from the metabolite D-glyceraldehyde 3-phosphate.
Circular Genome Viewer
A new circular genome viewer (Tools → Genome → Circular Genome Viewer) provides a global view of the organization of the chromosome as a set of concentric circles (tracks) containing features (genes, promoters, binding-sites, and other extragenic sites) of interest. A given track can be filtered at the outset to only show features that match certain criteria (the available selection criteria depend on the feature type) or it can include a larger set of features; various selection criteria can be applied after the fact to highlight subsets of features. Possible feature types that can be displayed include genes, pseudogenes, promoters, transcription factor binding sites, REP elements, and others. The set of filtering and highlighting criteria for genes include product type (e.g., RNAs, enzymes, and transporters), name substrings, pathway classes, regulons, GO terms, and gene identifiers from an uploaded file. Figure 4 shows an example display with a variety of feature types and highlights. The circular genome viewer can also combine tracks from multiple strains or related species and highlight the orthologs between them.
Figure 4. An example Circular Genome Viewer display, with tracks that showcase a variety of feature types, filtering and highlighting options, listed in order from the outermost circle inwards.
Revised Tools in EcoCyc
EcoCyc contains extensive web search options including a new command for searching for pseudogenes and different types of RNAs (website command Tools → Search → Genes, Proteins, or Tools → Search → RNAs → Search/Filter by type/subunits). We have also added a web-based search tool for searching for DNA and RNA sites of various types such as attenuators, riboswitches, phage attachment sites, and transposons (website command Tools → Search → Search DNA or mRNA sites).
We have upgraded the multiple-sequence alignment tools available for EcoCyc to use Clustal Omega (Sievers and Higgins, 2021) to compute alignments and MSA Viewer (Yachdav et al., 2016) to display the alignments (website command Analysis → Multiple Sequence Alignment).
The Genome Overview diagram depicts the entire E. coli gene in a single screen (Figure 5 and website command Tools → Genome → Genome Overview). Each gene is shown as a single arrow with an arrowhead style distinguishing protein-coding genes from RNA-coding genes, and arrow direction indicating transcription direction. Adjacent genes drawn in the same color are within the same operon. We recently added the ability to search the diagram for genes by name or by substring (e.g., find all the genes whose name contains “arg”) and to highlight the search results on the diagram.
Figure 5. The Genome Overview diagram captures an entire genome with each arrow representing a single gene; neighboring genes that are part of the same operon are displayed in the same color. Arrow sizes are not to scale.
The Regulatory Overview diagram depicts the E. coli regulatory network, more specifically, transcriptional regulation (including transcription factors and sigma factors), and translational regulation (including small RNAs). The diagram (Figure 6 and website command Tools → Genome → Regulatory Overview) is organized into three concentric ellipses; the inner ellipse depicts global regulatory genes, the middle ellipse depicts other regulatory genes, and the outer ring depicts genes that are not regulators. The diagram supports a variety of operations, including searching for genes by names and highlighting the regulators or regulatory targets of a given gene. A new command enables the user to output either the entire regulatory network or a subnetwork starting at a given gene to an ASCII file whose indentation describes the hierarchy of regulatory relationships.
Figure 6. The Regulatory Overview diagram captures user-specified portions of a regulatory network. The user has specified the display of regulation by FliA (red) and RhyB (yellow).
The Cellular Overview diagram depicts the full E. coli metabolic and transport network (see website command Tools → Metabolism → Cellular Overview). All EcoCyc pathways are included, grouped by class, along with a section for reactions that have not been assigned to pathways. Transporters and other membrane proteins are shown on a schematic of the double membrane, with periplasmic reactions and proteins between the membranes. The diagram supports highlighting operations for genes, proteins, metabolites, reactions, and pathways using a variety of criteria. This diagram is also used by the Omics Viewer, in which omics data, such as transcriptomics or metabolomics data, are overlaid on the cellular overview to illustrate experimental results in a metabolic context. The Omics Viewer has also been substantially revamped to give the user extensive interactive control over the mapping of omics data values to colors, including the ability to selectively hide or show specified data ranges.
All three of the overview diagrams have been re-engineered to use modern, high-quality graphics that draw more rapidly and to provide real-time semantic zooming capabilities.
Discussion
The EcoCyc database is unique in its extensive coverage of E. coli biology captured from a century of research. Ongoing manual curation enables the addition of new gene product functions and other important new research results, while the incorporation of new high-throughput datasets expands the types of data stored in the database. EcoCyc also welcomes user input. The “Provide Feedback” button on each data page can be used to submit information on new publications, to point out errors or omissions, and to suggest other improvements.
Future directions for EcoCyc include integrating EcoCyc with the E. coli whole cell model developed by the laboratory of Prof. M. Covert (Macklin et al., 2020) and improving the EcoCyc search and visualization tools.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: www.ecocyc.org.
Author Contributions
IK, SP, PK, AM, MK, JC-V, and AW: writing of manuscript. IK, AM, AS-Z, SG-C, VT, RC, and WO: EcoCyc curation. LM-R, CB-M, SP, MK, AK, and PM: EcoCyc data import. PS and RB: EcoCyc releases and website. SP, MK, WO, AK, PM, PS, and RB: Pathway Tools software development. PK, JC-V, and IP: guidance and oversight. PK and JC-V: funding. All authors contributed to the article and approved the submitted version.
Funding
This work was funded under awards from the National Institute of General Medical Sciences of the National Institutes of Health GM077678 to PK and RO1GM110597 to JC-V. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The NIH did not play any role in the design of the study; nor in collection, analysis, or interpretation of data; nor in writing the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Agrawal, S., Jaswal, K., Shiver, A. L., Balecha, H., Patra, T., and Chaba, R. (2017). A genome-wide screen in Escherichia coli reveals that ubiquinone is a key antioxidant for metabolism of long-chain fatty acids. J. Biol. Chem. 292, 20086–20099. doi: 10.1074/jbc.M117.806240
Aquino, P., Honda, B., Jaini, S., Lyubetskaya, A., Hosur, K., Chiu, J. G., et al. (2017). Coordinated regulation of acid resistance in Escherichia coli. BMC Syst. Biol. 11:1. doi: 10.1186/s12918-016-0376-y
Aussel, L., Loiseau, L., Hajj Chehade, M., Pocachard, B., Fontecave, M., Pierrel, F., et al. (2014). ubiJ, a new gene required for aerobic growth and proliferation in macrophage, is involved in coenzyme Q biosynthesis in Escherichia coli and Salmonella enterica serovar Typhimurium. J. Bacteriol. 196, 70–79. doi: 10.1128/JB.01065-13
Bradley, N. P., Washburn, L. A., Christov, P. P., Watanabe, C. M. H., and Eichman, B. F. (2020). Escherichia coli YcaQ is a DNA glycosylase that unhooks DNA interstrand crosslinks. Nucleic Acids Res. 48, 7005–7017. doi: 10.1093/nar/gkaa346
Caldara, M., Charlier, D., and Cunin, R. (2006). The arginine regulon of Escherichia coli: whole-system transcriptome analysis discovers new genes and provides an integrated view of arginine regulation. Microbiology 152, 3343–3354. doi: 10.1099/mic.0.29088-0
Chitsaz, M., and Brown, M. H. (2017). The role played by drug efflux pumps in bacterial multidrug resistance. Essays Biochem. 61, 127–139. doi: 10.1042/EBC20160064
Cho, B.-K., Barrett, C. L., Knight, E. M., Park, Y. S., and Palsson, B. Ø. (2008). Genome-scale reconstruction of the Lrp regulatory network in Escherichia coli. Proc. Natl. Acad. Sci. U. S. A. 105, 19462–19467. doi: 10.1073/pnas.0807227105
Cho, S., Cho, Y.-B., Kang, T. J., Kim, S. C., Palsson, B., and Cho, B.-K. (2015). The architecture of ArgR-DNA complexes at the genome-scale in Escherichia coli. Nucleic Acids Res. 43, 3079–3088. doi: 10.1093/nar/gkv150
Christensen, D. G., Baumgartner, J. T., Xie, X., Jew, K. M., Basisty, N., Schilling, B., et al. (2019). Mechanisms, detection, and relevance of protein acetylation in prokaryotes. MBio 10, e02708–e02718. doi: 10.1128/mBio.02708-18
Christensen, D. G., Meyer, J. G., Baumgartner, J. T., D’Souza, A. K., Nelson, W. C., Payne, S. H., et al. (2018). Identification of novel protein lysine acetyltransferases in Escherichia coli. MBio 9, e01905–e01918. doi: 10.1128/mBio.01905-18
Christensen, D. G., Orr, J. S., Rao, C. V., and Wolfe, A. J. (2017). Increasing growth yield and decreasing acetylation in Escherichia coli by optimizing the carbon-to-magnesium ratio in peptide-based media. Appl. Environ. Microbiol. 83, e03034–e03016. doi: 10.1128/AEM.03034-16
Cox, G. B., Young, I. G., McCann, L. M., and Gibson, F. (1969). Biosynthesis of ubiquinone in Escherichia coli K-12: location of genes affecting the metabolism of 3-octaprenyl-4-hydroxybenzoic acid and 2-octaprenylphenol. J. Bacteriol. 99, 450–458. doi: 10.1128/jb.99.2.450-458.1969
Davis, R., Écija-Conesa, A., Gallego-Jara, J., de Diego, T., Filippova, E. V., Kuffel, G., et al. (2018). An acetylatable lysine controls CRP function in E. coli. Mol. Microbiol. 107, 116–131. doi: 10.1111/mmi.13874
Demerec, M., Adelberg, E. A., Clark, A. J., and Hartman, P. E. (1966). A proposal for a uniform nomenclature in bacterial genetics. Genetics 54, 61–76. doi: 10.1093/genetics/54.1.61
Federowicz, S., Kim, D., Ebrahim, A., Lerman, J., Nagarajan, H., Cho, B., et al. (2014). Determining the control circuitry of redox metabolism at the genome-scale. PLoS Genet. 10:e1004264. doi: 10.1371/journal.pgen.1004264
Gasperotti, A., Göing, S., Fajardo-Ruiz, E., Forné, I., and Jung, K. (2020). Function and regulation of the pyruvate transporter CstA in Escherichia coli. Int. J. Mol. Sci. 21:E9068. doi: 10.3390/ijms21239068
Gati, C., Stetsenko, A., Slotboom, D. J., Scheres, S. H. W., and Guskov, A. (2017). The structural basis of proton driven zinc transport by ZntB. Nat. Commun. 8:1313. doi: 10.1038/s41467-017-01483-7
Ghatak, S., King, Z. A., Sastry, A., and Palsson, B. O. (2019). The y-ome defines the 35% of Escherichia coli genes that lack experimental evidence of function. Nucleic Acids Res. 47, 2446–2454. doi: 10.1093/nar/gkz030
Hajj Chehade, M., Loiseau, L., Lombard, M., Pecqueur, L., Ismail, A., Smadja, M., et al. (2013). ubiI, a new gene in Escherichia coli coenzyme Q biosynthesis, is involved in aerobic C5-hydroxylation. J. Biol. Chem. 288, 20085–20092. doi: 10.1074/jbc.M113.480368
Hajj Chehade, M., Pelosi, L., Fyfe, C. D., Loiseau, L., Rascalou, B., Brugière, S., et al. (2019). A soluble metabolon synthesizes the isoprenoid lipid ubiquinone. Cell Chem. Biol. 26, 482.e7–492.e7. doi: 10.1016/j.chembiol.2018.12.001
Huang, L., Wang, J., and Lilley, D. M. J. (2017). The structure of the guanidine-II riboswitch. Cell Chem. Biol. 24, 695.e2–702.e2. doi: 10.1016/j.chembiol.2017.05.014
Hwang, S., Choe, D., Yoo, M., Cho, S., Kim, S. C., Cho, S., et al. (2018). Peptide transporter CstA imports pyruvate in Escherichia coli K-12. J. Bacteriol. 200, e00771–e00717. doi: 10.1128/JB.00771-17
Jacob, F., and Monod, J. (1961). Genetic regulatory mechanisms in the synthesis of proteins. J. Mol. Biol. 3, 318–356. doi: 10.1016/S0022-2836(61)80072-7
Ju, X., Li, D., and Liu, S. (2019). Full-length RNA profiling reveals pervasive bidirectional transcription terminators in bacteria. Nat. Microbiol. 4, 1907–1918. doi: 10.1038/s41564-019-0500-z
Karp, P. D., and Riley, M. (1993). Representations of metabolic knowledge: pathways. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1, 207–215.
Kermani, A. A., Macdonald, C. B., Gundepudi, R., and Stockbridge, R. B. (2018). Guanidinium export is the primal function of SMR family transporters. Proc. Natl. Acad. Sci. U. S. A. 115, 3060–3065. doi: 10.1073/pnas.1719187115
Keseler, I. M., Mackie, A., Santos-Zavaleta, A., Billington, R., Bonavides-Martínez, C., Caspi, R., et al. (2017). The EcoCyc database: reflecting new knowledge about Escherichia coli K-12. Nucleic Acids Res. 45, D543–D550. doi: 10.1093/nar/gkw1003
Kristoficova, I., Vilhena, C., Behr, S., and Jung, K. (2018). BtsT, a novel and specific pyruvate/H+ symporter in Escherichia coli. J. Bacteriol. 200, e00599–e00517. doi: 10.1128/JB.00599-17
Kroner, G. M., Wolfe, M. B., and Freddolino, P. L. (2019). Escherichia coli Lrp regulates one-third of the genome via direct, cooperative, and indirect routes. J. Bacteriol. 201, e00411–e00418. doi: 10.1128/JB.00411-18
Kuhn, M. L., Zemaitaitis, B., Hu, L. I., Sahu, A., Sorensen, D., Minasov, G., et al. (2014). Structural, kinetic and proteomic characterization of acetyl phosphate-dependent bacterial protein acetylation. PLoS One 9:e94816. doi: 10.1371/journal.pone.0094816
Liu, J., and Turnbough, C. L. (1994). Effects of transcriptional start site sequence and position on nucleotide-sensitive selection of alternative start sites at the pyrC promoter in Escherichia coli. J. Bacteriol. 176, 2938–2945. doi: 10.1128/jb.176.10.2938-2945.1994
Loiseau, L., Fyfe, C., Aussel, L., Hajj Chehade, M., Hernández, S. B., Faivre, B., et al. (2017). The UbiK protein is an accessory factor necessary for bacterial ubiquinone (UQ) biosynthesis and forms a complex with the UQ biogenesis factor UbiJ. J. Biol. Chem. 292, 11937–11950. doi: 10.1074/jbc.M117.789164
Macklin, D. N., Ahn-Horst, T. A., Choi, H., Ruggero, N. A., Carrera, J., Mason, J. C., et al. (2020). Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation. Science 369:eaav3751. doi: 10.1126/science.aav3751
Mejía-Almonte, C., Busby, S. J. W., Wade, J. T., van Helden, J., Arkin, A. P., Stormo, G. D., et al. (2020). Redefining fundamental concepts of transcription initiation in bacteria. Nat. Rev. Genet. 21, 699–714. doi: 10.1038/s41576-020-0254-8
Mohni, K. N., Wessel, S. R., Zhao, R., Wojciechowski, A. C., Luzwick, J. W., Layden, H., et al. (2019). HMCES maintains genome integrity by shielding abasic sites in single-strand DNA. Cell 176, 144.e13–153.e13. doi: 10.1016/j.cell.2018.10.055
Nakayasu, E. S., Burnet, M. C., Walukiewicz, H. E., Wilkins, C. S., Shukla, A. K., Brooks, S., et al. (2017). Ancient regulatory role of lysine acetylation in central metabolism. MBio 8, e01894–e01817. doi: 10.1128/mBio.01894-17
Paley, S., Keseler, I. M., Krummenacker, M., and Karp, P. D. (2021). Leveraging curation among Escherichia coli pathway/genome databases using ortholog-based annotation propagation. Front. Microbiol. 12:614355. doi: 10.3389/fmicb.2021.614355
Poon, W. W., Davis, D. E., Ha, H. T., Jonassen, T., Rather, P. N., and Clarke, C. F. (2000). Identification of Escherichia coli ubiB, a gene required for the first monooxygenase step in ubiquinone biosynthesis. J. Bacteriol. 182, 5139–5146. doi: 10.1128/JB.182.18.5139-5146.2000
Ren, J., Sang, Y., Qin, R., Su, Y., Cui, Z., Mang, Z., et al. (2019). Metabolic intermediate acetyl phosphate modulates bacterial virulence via acetylation. Emerg. Microbes Infect. 8, 55–69. doi: 10.1080/22221751.2018.1558963
Santos-Zavaleta, A., Sánchez-Pérez, M., Salgado, H., Velázquez-Ramírez, D. A., Gama-Castro, S., Tierrafría, V. H., et al. (2018). A unified resource for transcriptional regulation in Escherichia coli K-12 incorporating high-throughput-generated binding data into RegulonDB version 10.0. BMC Biol. 16:91. doi: 10.1186/s12915-018-0555-y
Schilling, B., Christensen, D., Davis, R., Sahu, A. K., Hu, L. I., Walker-Peddakotla, A., et al. (2015). Protein acetylation dynamics in response to carbon overflow in Escherichia coli. Mol. Microbiol. 98, 847–863. doi: 10.1111/mmi.13161
Sherlock, M. E., Malkowski, S. N., and Breaker, R. R. (2017). Biochemical validation of a second guanidine riboswitch class in bacteria. Biochemistry 56, 352–358. doi: 10.1021/acs.biochem.6b01270
Shimada, T., Takada, H., Yamamoto, K., and Ishihama, A. (2015). Expanded roles of two-component response regulator OmpR in Escherichia coli: genomic SELEX search for novel regulation targets. Genes Cells Devoted Mol. Cell. Mech. 20, 915–931. doi: 10.1111/gtc.12282
Sievers, F., and Higgins, D. G. (2021). The clustal omega multiple alignment package. Methods Mol. Biol. 2231, 3–16. doi: 10.1007/978-1-0716-1036-7_1
Starai, V. J., and Escalante-Semerena, J. C. (2004). Identification of the protein acetyltransferase (pat) enzyme that acetylates acetyl-CoA synthetase in Salmonella enterica. J. Mol. Biol. 340, 1005–1012. doi: 10.1016/j.jmb.2004.05.010
Thompson, P. S., Amidon, K. M., Mohni, K. N., Cortez, D., and Eichman, B. F. (2019). Protection of abasic sites during DNA replication by a stable thiazolidine protein-DNA cross-link. Nat. Struct. Mol. Biol. 26, 613–618. doi: 10.1038/s41594-019-0255-5
Venkat, S., Gregory, C., Sturges, J., Gan, Q., and Fan, C. (2017). Studying the lysine acetylation of malate dehydrogenase. J. Mol. Biol. 429, 1396–1405. doi: 10.1016/j.jmb.2017.03.027
Walker, K. A., and Osuna, R. (2002). Factors affecting start site selection at the Escherichia coli fis promoter. J. Bacteriol. 184, 4783–4791. doi: 10.1128/JB.184.17.4783-4791.2002
Wang, N., Bao, H., Chen, L., Liu, Y., Li, Y., Wu, B., et al. (2019). Molecular basis of abasic site sensing in single-stranded DNA by the SRAP domain of E. coli yedK. Nucleic Acids Res. 47, 10388–10399. doi: 10.1093/nar/gkz744
Weinert, B. T., Iesmantavicius, V., Wagner, S. A., Schölz, C., Gummesson, B., Beli, P., et al. (2013). Acetyl-phosphate is a critical determinant of lysine acetylation in E. coli. Mol. Cell 51, 265–272. doi: 10.1016/j.molcel.2013.06.003
Winkelman, J. T., Vvedenskaya, I. O., Zhang, Y., Zhang, Y., Bird, J. G., Taylor, D. M., et al. (2016). Multiplexed protein-DNA cross-linking: scrunching in transcription start site selection. Science 351, 1090–1093. doi: 10.1126/science.aad6881
Yachdav, G., Wilzbach, S., Rauscher, B., Sheridan, R., Sillitoe, I., Procter, J., et al. (2016). MSAViewer: interactive JavaScript visualization of multiple sequence alignments. Bioinformatics 32, 3501–3503. doi: 10.1093/bioinformatics/btw474
Keywords: Escherichia coli, EcoCyc, model-organism database, drug efflux transporters, metabolism, gene regulation
Citation: Keseler IM, Gama-Castro S, Mackie A, Billington R, Bonavides-Martínez C, Caspi R, Kothari A, Krummenacker M, Midford PE, Muñiz-Rascado L, Ong WK, Paley S, Santos-Zavaleta A, Subhraveti P, Tierrafría VH, Wolfe AJ, Collado-Vides J, Paulsen IT and Karp PD (2021) The EcoCyc Database in 2021. Front. Microbiol. 12:711077. doi: 10.3389/fmicb.2021.711077
Edited by:
Jörg Stülke, University of Göttingen, GermanyReviewed by:
Jürgen Lassak, Ludwig Maximilian University of Munich, GermanyMichael Y. Galperin, National Center for Biotechnology Information (NLM), United States
Deborah A. Siegele, Texas A&M University, United States
Copyright © 2021 Keseler, Gama-Castro, Mackie, Billington, Bonavides-Martínez, Caspi, Kothari, Krummenacker, Midford, Muñiz-Rascado, Ong, Paley, Santos-Zavaleta, Subhraveti, Tierrafría, Wolfe, Collado-Vides, Paulsen and Karp. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peter D. Karp, cGthcnBAYWkuc3JpLmNvbQ==
†Present address: Wai Kit Ong, X, The Moonshot Factory, Mountain View, CA, United States