Designing Biological Circuits: Synthetic Biology Within the Operon Model and Beyond.

In 1961, Jacob and Monod proposed the operon model of gene regulation. At the model's core was the modular assembly of regulators, operators, and structural genes. To illustrate the composability of these elements, Jacob and Monod linked phenotypic diversity to the architectures of regulatory circuits. In this review, we examine how the circuit blueprints imagined by Jacob and Monod laid the foundation for the first synthetic gene networks that launched the field of synthetic biology in 2000. We discuss the influences of the operon model and its broader theoretical framework on the first generation of synthetic biological circuits, which were predominantly transcriptional and posttranscriptional circuits. We also describe how recent advances in molecular biology beyond the operon model-namely, programmable DNA- and RNA-binding molecules as well as models of epigenetic and posttranslational regulation-are expanding the synthetic biology toolkit and enabling the design of more complex biological circuits.

A natural single-guide RNA repurposes Cas9 to autoregulate CRISPR-Cas expression.

CRISPR-Cas systems provide prokaryotes with acquired immunity against viruses and plasmids, but how these systems are regulated to prevent autoimmunity is poorly understood. Here, we show that in the S. pyogenes CRISPR-Cas system, a long-form transactivating CRISPR RNA (tracr-L) folds into a natural single guide that directs Cas9 to transcriptionally repress its own promoter (Pcas). Further, we demonstrate that Pcas serves as a critical regulatory node. De-repression causes a dramatic 3,000-fold increase in immunization rates against viruses; however, heightened immunity comes at the cost of increased autoimmune toxicity. Using bioinformatic analyses, we provide evidence that tracrRNA-mediated autoregulation is widespread in type II-A CRISPR-Cas systems. Collectively, we unveil a new paradigm for the intrinsic regulation of CRISPR-Cas systems by natural single guides, which may facilitate the frequent horizontal transfer of these systems into new hosts that have not yet evolved their own regulatory strategies.

A Metabolic Pathway for Activation of Dietary Glucosinolates by a Human Gut Symbiont.

Consumption of glucosinolates, pro-drug-like metabolites abundant in Brassica vegetables, has been associated with decreased risk of certain cancers. Gut microbiota have the ability to metabolize glucosinolates, generating chemopreventive isothiocyanates. Here, we identify a genetic and biochemical basis for activation of glucosinolates to isothiocyanates by Bacteroides thetaiotaomicron, a prominent gut commensal species. Using a genome-wide transposon insertion screen, we identified an operon required for glucosinolate metabolism in B. thetaiotaomicron. Expression of BT2159-BT2156 in a non-metabolizing relative, Bacteroides fragilis, resulted in gain of glucosinolate metabolism. We show that isothiocyanate formation requires the action of BT2158 and either BT2156 or BT2157 in vitro. Monocolonization of mice with mutant BtΔ2157 showed reduced isothiocyanate production in the gastrointestinal tract. These data provide insight into the mechanisms by which a common gut bacterium processes an important dietary nutrient.

Some Like it HOT: Horizontal Operon Transfer.

While horizontal gene transfer (HGT) is well documented in bacteria, the role and frequency of HGT across eukaryotes remains poorly understood. Kominek et al. identified a horizontal operon transfer (HOT) event, with clear evidence for selection to facilitate gene expression, that has allowed a group of yeasts to scavenge iron using bacterially derived genes.

Eukaryotic Acquisition of a Bacterial Operon.

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Chalkophores.

Copper-binding metallophores, or chalkophores, play a role in microbial copper homeostasis that is analogous to that of siderophores in iron homeostasis. The best-studied chalkophores are members of the methanobactin (Mbn) family-ribosomally produced, posttranslationally modified natural products first identified as copper chelators responsible for copper uptake in methane-oxidizing bacteria. To date, Mbns have been characterized exclusively in those species, but there is genomic evidence for their production in a much wider range of bacteria. This review addresses the current state of knowledge regarding the function, biosynthesis, transport, and regulation of Mbns. While the roles of several proteins in these processes are supported by substantial genetic and biochemical evidence, key aspects of Mbn manufacture, handling, and regulation remain unclear. In addition, other natural products that have been proposed to mediate copper uptake as well as metallophores that have biologically relevant roles involving copper binding, but not copper uptake, are discussed.

Evolutionary Convergence of Pathway-Specific Enzyme Expression Stoichiometry.

Coexpression of proteins in response to pathway-inducing signals is the founding paradigm of gene regulation. However, it remains unexplored whether the relative abundance of co-regulated proteins requires precise tuning. Here, we present large-scale analyses of protein stoichiometry and corresponding regulatory strategies for 21 pathways and 67-224 operons in divergent bacteria separated by 0.6-2 billion years. Using end-enriched RNA-sequencing (Rend-seq) with single-nucleotide resolution, we found that many bacterial gene clusters encoding conserved pathways have undergone massive divergence in transcript abundance and architectures via remodeling of internal promoters and terminators. Remarkably, these evolutionary changes are compensated post-transcriptionally to maintain preferred stoichiometry of protein synthesis rates. Even more strikingly, in eukaryotic budding yeast, functionally analogous proteins that arose independently from bacterial counterparts also evolved to convergent in-pathway expression. The broad requirement for exact protein stoichiometries despite regulatory divergence provides an unexpected principle for building biological pathways both in nature and for synthetic activities.

Transcription-replication interactions reveal bacterial genome regulation.

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

A eukaryotic-like ubiquitination system in bacterial antiviral defence.

Ubiquitination pathways have crucial roles in protein homeostasis, signalling and innate immunity1-3. In these pathways, an enzymatic cascade of E1, E2 and E3 proteins conjugates ubiquitin or a ubiquitin-like protein (Ubl) to target-protein lysine residues4. Bacteria encode ancient relatives of E1 and Ubl proteins involved in sulfur metabolism5,6, but these proteins do not mediate Ubl-target conjugation, leaving open the question of whether bacteria can perform ubiquitination-like protein conjugation. Here we demonstrate that a bacterial operon associated with phage defence islands encodes a complete ubiquitination pathway. Two structures of a bacterial E1-E2-Ubl complex reveal striking architectural parallels with canonical eukaryotic ubiquitination machinery. The bacterial E1 possesses an amino-terminal inactive adenylation domain and a carboxy-terminal active adenylation domain with a mobile α-helical insertion containing the catalytic cysteine (CYS domain). One structure reveals a pre-reaction state with the bacterial Ubl C terminus positioned for adenylation, and a second structure mimics an E1-to-E2 transthioesterification state with the E1 CYS domain adjacent to the bound E2. We show that a deubiquitinase in the same pathway preprocesses the bacterial Ubl, exposing its C-terminal glycine for adenylation. Finally, we show that the bacterial E1 and E2 collaborate to conjugate Ubl to target-protein lysine residues. Together, these data reveal that bacteria possess bona fide ubiquitination systems with strong mechanistic and architectural parallels to canonical eukaryotic ubiquitination pathways, suggesting that these pathways arose first in bacteria.

Antiviral signalling by a cyclic nucleotide activated CRISPR protease.

CRISPR defence systems such as the well-known DNA-targeting Cas9 and the RNA-targeting type III systems are widespread in prokaryotes1,2. The latter orchestrates a complex antiviral response that is initiated through the synthesis of cyclic oligoadenylates after recognition of foreign RNA3-5. Among the large set of proteins that are linked to type III systems and predicted to bind cyclic oligoadenylates6,7, a CRISPR-associated Lon protease (CalpL) stood out to us. CalpL contains a sensor domain of the SAVED family7 fused to a Lon protease effector domain. However, the mode of action of this effector is unknown. Here we report the structure and function of CalpL and show that this soluble protein forms a stable tripartite complex with two other proteins, CalpT and CalpS, that are encoded on the same operon. After activation by cyclic tetra-adenylate (cA4), CalpL oligomerizes and specifically cleaves the MazF homologue CalpT, which releases the extracytoplasmic function σ factor CalpS from the complex. Our data provide a direct connection between CRISPR-based detection of foreign nucleic acids and transcriptional regulation. Furthermore, the presence of a SAVED domain that binds cyclic tetra-adenylate in a CRISPR effector reveals a link to the cyclic-oligonucleotide-based antiphage signalling system.

Small RNA-mediated activation of sugar phosphatase mRNA regulates glucose homeostasis.

Glucose homeostasis is strictly controlled in all domains of life. Bacteria that are unable to balance intracellular sugar levels and deal with potentially toxic phosphosugars cease growth and risk being outcompeted. Here, we identify the conserved haloacid dehalogenase (HAD)-like enzyme YigL as the previously hypothesized phosphatase for detoxification of phosphosugars and reveal that its synthesis is activated by an Hfq-dependent small RNA in Salmonella typhimurium. We show that the glucose-6-P-responsive small RNA SgrS activates YigL synthesis in a translation-independent fashion by the selective stabilization of a decay intermediate of the dicistronic pldB-yigL messenger RNA (mRNA). Intriguingly, the major endoribonuclease RNase E, previously known to function together with small RNAs to degrade mRNA targets, is also essential for this process of mRNA activation. The exploitation of and targeted interference with regular RNA turnover described here may constitute a general route for small RNAs to rapidly activate both coding and noncoding genes.

Cryo-EM structure of an active bacterial TIR-STING filament complex.

Stimulator of interferon genes (STING) is an antiviral signalling protein that is broadly conserved in both innate immunity in animals and phage defence in prokaryotes1-4. Activation of STING requires its assembly into an oligomeric filament structure through binding of a cyclic dinucleotide4-13, but the molecular basis of STING filament assembly and extension remains unknown. Here we use cryogenic electron microscopy to determine the structure of the active Toll/interleukin-1 receptor (TIR)-STING filament complex from a Sphingobacterium faecium cyclic-oligonucleotide-based antiphage signalling system (CBASS) defence operon. Bacterial TIR-STING filament formation is driven by STING interfaces that become exposed on high-affinity recognition of the cognate cyclic dinucleotide signal c-di-GMP. Repeating dimeric STING units stack laterally head-to-head through surface interfaces, which are also essential for human STING tetramer formation and downstream immune signalling in mammals5. The active bacterial TIR-STING structure reveals further cross-filament contacts that brace the assembly and coordinate packing of the associated TIR NADase effector domains at the base of the filament to drive NAD+ hydrolysis. STING interface and cross-filament contacts are essential for cell growth arrest in vivo and reveal a stepwise mechanism of activation whereby STING filament assembly is required for subsequent effector activation. Our results define the structural basis of STING filament formation in prokaryotic antiviral signalling.

An α helix to ß barrel domain switch transforms the transcription factor RfaH into a translation factor.

NusG homologs regulate transcription and coupled processes in all living organisms. The Escherichia coli (E. coli) two-domain paralogs NusG and RfaH have conformationally identical N-terminal domains (NTDs) but dramatically different carboxy-terminal domains (CTDs), a ß barrel in NusG and an α hairpin in RfaH. Both NTDs interact with elongating RNA polymerase (RNAP) to reduce pausing. In NusG, NTD and CTD are completely independent, and NusG-CTD interacts with termination factor Rho or ribosomal protein S10. In contrast, RfaH-CTD makes extensive contacts with RfaH-NTD to mask an RNAP-binding site therein. Upon RfaH interaction with its DNA target, the operon polarity suppressor (ops) DNA, RfaH-CTD is released, allowing RfaH-NTD to bind to RNAP. Here, we show that the released RfaH-CTD completely refolds from an all-α to an all-ß conformation identical to that of NusG-CTD. As a consequence, RfaH-CTD binding to S10 is enabled and translation of RfaH-controlled operons is strongly potentiated. PAPERFLICK:

Cross-regulation and cross-talk of conserved and accessory two-component regulatory systems orchestrate Pseudomonas copper resistance.

Bacteria use diverse strategies and molecular machinery to maintain copper homeostasis and to cope with its toxic effects. Some genetic elements providing copper resistance are acquired by horizontal gene transfer; however, little is known about how they are controlled and integrated into the central regulatory network. Here, we studied two copper-responsive systems in a clinical isolate of Pseudomonas paraeruginosa and deciphered the regulatory and cross-regulation mechanisms. To do so, we combined mutagenesis, transcriptional fusion analyses and copper sensitivity phenotypes. Our results showed that the accessory CusRS two-component system (TCS) responds to copper and activates both its own expression and that of the adjacent nine-gene operon (the pcoA2 operon) to provide resistance to elevated levels of extracellular copper. The same locus was also found to be regulated by two core-genome-encoded TCSs-the copper-responsive CopRS and the zinc-responsive CzcRS. Although the target palindromic sequence-ATTCATnnATGTAAT-is the same for the three response regulators, transcriptional outcomes differ. Thus, depending on the operon/regulator pair, binding can result in different activation levels (from none to high), with the systems demonstrating considerable plasticity. Unexpectedly, although the classical CusRS and the noncanonical CopRS TCSs rely on distinct signaling mechanisms (kinase-based vs. phosphatase-based), we discovered cross-talk in the absence of the cognate sensory kinases. This cross-talk occurred between the proteins of these two otherwise independent systems. The cusRS-pcoA2 locus is part of an Integrative and Conjugative Element and was found in other Pseudomonas strains where its expression could provide copper resistance under appropriate conditions. The results presented here illustrate how acquired genetic elements can become part of endogenous regulatory networks, providing a physiological advantage. They also highlight the potential for broader effects of accessory regulatory proteins through interference with core regulatory proteins.

Reciprocal regulation of enterococcal cephalosporin resistance by products of the autoregulated yvcJ-glmR-yvcL operon enhances fitness during cephalosporin exposure.

Enterococci are commensal members of the gastrointestinal tract and also major nosocomial pathogens. They possess both intrinsic and acquired resistance to many antibiotics, including intrinsic resistance to cephalosporins that target bacterial cell wall synthesis. These antimicrobial resistance traits make enterococcal infections challenging to treat. Moreover, prior therapy with antibiotics, including broad-spectrum cephalosporins, promotes enterococcal proliferation in the gut, resulting in dissemination to other sites of the body and subsequent infection. As a result, a better understanding of mechanisms of cephalosporin resistance is needed to enable development of new therapies to treat or prevent enterococcal infections. We previously reported that flow of metabolites through the peptidoglycan biosynthesis pathway is one determinant of enterococcal cephalosporin resistance. One factor that has been implicated in regulating flow of metabolites into cell wall biosynthesis pathways of other Gram-positive bacteria is GlmR. In enterococci, GlmR is encoded as the middle gene of a predicted 3-gene operon along with YvcJ and YvcL, whose functions are poorly understood. Here we use genetics and biochemistry to investigate the function of the enterococcal yvcJ-glmR-yvcL gene cluster. Our results reveal that YvcL is a DNA-binding protein that regulates expression of the yvcJ-glmR-yvcL operon in response to cell wall stress. YvcJ and GlmR bind UDP-GlcNAc and reciprocally regulate cephalosporin resistance in E. faecalis, and binding of UDP-GlcNAc by YvcJ appears essential for its activity. Reciprocal regulation by YvcJ/GlmR is essential for fitness during exposure to cephalosporin stress. Additionally, our results indicate that enterococcal GlmR likely acts by a different mechanism than the previously studied GlmR of Bacillus subtilis, suggesting that the YvcJ/GlmR regulatory module has evolved unique targets in different species of bacteria.

A multi-iron enzyme installs copper-binding oxazolone/thioamide pairs on a nontypeable Haemophilus influenzae virulence factor.

The multinuclear nonheme iron-dependent oxidases (MNIOs) are a rapidly growing family of enzymes involved in the biosynthesis of ribosomally synthesized, posttranslationally modified peptide natural products (RiPPs). Recently, a secreted virulence factor from nontypeable Haemophilus influenzae (NTHi) was found to be expressed from an operon, which we designate the hvf operon, that also encodes an MNIO. Here, we show by Mössbauer spectroscopy that the MNIO HvfB contains a triiron cofactor. We demonstrate that HvfB works together with HvfC [a RiPP recognition element (RRE)-containing partner protein] to perform six posttranslational modifications of cysteine residues on the virulence factor precursor peptide HvfA. Structural characterization by tandem mass spectrometry and NMR shows that these six cysteine residues are converted to oxazolone and thioamide pairs, similar to those found in the RiPP methanobactin. Like methanobactin, the mature virulence factor, which we name oxazolin, uses these modified residues to coordinate Cu(I) ions. Considering the necessity of oxazolin for host cell invasion by NTHi, these findings point to a key role for copper during NTHi infection. Furthermore, oxazolin and its biosynthetic pathway represent a potential therapeutic target for NTHi.

Bacteria elicit a phage tolerance response subsequent to infection of their neighbors.

Appearance of plaques on a bacterial lawn is a sign of successive rounds of bacteriophage infection. Yet, mechanisms evolved by bacteria to limit plaque spread have been hardly explored. Here, we investigated the dynamics of plaque development by lytic phages infecting the bacterium Bacillus subtilis. We report that plaque expansion is followed by a constriction phase owing to bacterial growth into the plaque zone. This phenomenon exposed an adaptive process, herein termed "phage tolerance response", elicited by non-infected bacteria upon sensing infection of their neighbors. The temporary phage tolerance is executed by the stress-response RNA polymerase sigma factor σX (SigX). Artificial expression of SigX prior to phage attack largely eliminates infection. SigX tolerance is primarily conferred by activation of the dlt operon, encoding enzymes that catalyze D-alanylation of cell wall teichoic acid polymers, the major attachment sites for phages infecting Gram-positive bacteria. D-alanylation impedes phage binding and hence infection, thus enabling the uninfected bacteria to form a protective shield opposing phage spread.

Metabolic Gene Clusters in Eukaryotes.

In bacteria, more than half of the genes in the genome are organized in operons. In contrast, in eukaryotes, functionally related genes are usually dispersed across the genome. There are, however, numerous examples of functional clusters of nonhomologous genes for metabolic pathways in fungi and plants. Despite superficial similarities with operons (physical clustering, coordinate regulation), these clusters have not usually originated by horizontal gene transfer from bacteria, and (unlike operons) the genes are typically transcribed separately rather than as a single polycistronic message. This clustering phenomenon raises intriguing questions about the origins of clustered metabolic pathways in eukaryotes and the significance of clustering for pathway function. Here we review metabolic gene clusters from fungi and plants, highlight commonalities and differences, and consider how these clusters form and are regulated. We also identify opportunities for future research in the areas of large-scale genomics, synthetic biology, and experimental evolution.

Pseudomonas aeruginosa senses and responds to epithelial potassium flux via Kdp operon to promote biofilm.

Mucosa-associated biofilms are associated with many human disease states, but the host mechanisms promoting biofilm remain unclear. In chronic respiratory diseases like cystic fibrosis (CF), Pseudomonas aeruginosa establishes chronic infection through biofilm formation. P. aeruginosa can be attracted to interspecies biofilms through potassium currents emanating from the biofilms. We hypothesized that P. aeruginosa could, similarly, sense and respond to the potassium efflux from human airway epithelial cells (AECs) to promote biofilm. Using respiratory epithelial co-culture biofilm imaging assays of P. aeruginosa grown in association with CF bronchial epithelial cells (CFBE41o-), we found that P. aeruginosa biofilm was increased by potassium efflux from AECs, as examined by potentiating large conductance potassium channel, BKCa (NS19504) potassium efflux. This phenotype is driven by increased bacterial attachment and increased coalescence of bacteria into aggregates. Conversely, biofilm formation was reduced when AECs were treated with a BKCa blocker (paxilline). Using an agar-based macroscopic chemotaxis assay, we determined that P. aeruginosa chemotaxes toward potassium and screened transposon mutants to discover that disruption of the high-sensitivity potassium transporter, KdpFABC, and the two-component potassium sensing system, KdpDE, reduces P. aeruginosa potassium chemotaxis. In respiratory epithelial co-culture biofilm imaging assays, a KdpFABCDE deficient P. aeruginosa strain demonstrated reduced biofilm growth in association with AECs while maintaining biofilm formation on abiotic surfaces. Furthermore, we determined that the Kdp operon is expressed in vivo in people with CF and the genes are conserved in CF isolates. Collectively, these data suggest that P. aeruginosa biofilm formation can be increased by attracting bacteria to the mucosal surface and enhancing coalescence into microcolonies through aberrant AEC potassium efflux sensed by the KdpFABCDE system. These findings suggest host electrochemical signaling can enhance biofilm, a novel host-pathogen interaction, and potassium flux could be a therapeutic target to prevent chronic infections in diseases with mucosa-associated biofilms, like CF.

Quantitative Control for Stoichiometric Protein Synthesis.

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production.