Tradeoffs and optimality in the evolution of gene regulation.

Cellular regulation is believed to evolve in response to environmental variability. However, this has been difficult to test directly. Here, we show that a gene regulation system evolves to the optimal regulatory response when challenged with variable environments. We engineered a genetic module subject to regulation by the lac repressor (LacI) in E. coli, whose expression is beneficial in one environmental condition and detrimental in another. Measured tradeoffs in fitness between environments predict the competition between regulatory phenotypes. We show that regulatory evolution in adverse environments is delayed at specific boundaries in the phenotype space of the regulatory LacI protein. Once this constraint is relieved by mutation, adaptation proceeds toward the optimum, yielding LacI with an altered allosteric mechanism that enables an opposite response to its regulatory ligand IPTG. Our results indicate that regulatory evolution can be understood in terms of tradeoff optimization theory.

DNA surface exploration and operator bypassing during target search.

Many proteins that bind specific DNA sequences search the genome by combining three-dimensional diffusion with one-dimensional sliding on nonspecific DNA1-5. Here we combine resonance energy transfer and fluorescence correlation measurements to characterize how individual lac repressor (LacI) molecules explore the DNA surface during the one-dimensional phase of target search. To track the rotation of sliding LacI molecules on the microsecond timescale, we use real-time single-molecule confocal laser tracking combined with fluorescence correlation spectroscopy (SMCT-FCS). The fluctuations in fluorescence signal are accurately described by rotation-coupled sliding, in which LacI traverses about 40 base pairs (bp) per revolution. This distance substantially exceeds the 10.5-bp helical pitch of DNA; this suggests that the sliding protein frequently hops out of the DNA groove, which would result in the frequent bypassing of target sequences. We directly observe such bypassing using single-molecule fluorescence resonance energy transfer (smFRET). A combined analysis of the smFRET and SMCT-FCS data shows that LacI hops one or two grooves (10-20 bp) every 200-700 µs. Our data suggest a trade-off between speed and accuracy during sliding: the weak nature of nonspecific protein-DNA interactions underlies operator bypassing, but also speeds up sliding. We anticipate that SMCT-FCS, which monitors rotational diffusion on the microsecond timescale while tracking individual molecules with millisecond resolution, will be applicable to the real-time investigation of many other biological interactions and will effectively extend the accessible time regime for observing these interactions by two orders of magnitude.

Engineered transcription activator-like effector dimer proteins confer DNA loop-dependent gene repression comparable to Lac repressor.

Natural prokaryotic gene repression systems often exploit DNA looping to increase the local concentration of gene repressor proteins at a regulated promoter via contributions from repressor proteins bound at distant sites. Using principles from the Escherichia coli lac operon we design analogous repression systems based on target sequence-programmable Transcription Activator-Like Effector dimer (TALED) proteins. Such engineered switches may be valuable for synthetic biology and therapeutic applications. Previous TALEDs with inducible non-covalent dimerization showed detectable, but limited, DNA loop-based repression due to the repressor protein dimerization equilibrium. Here, we show robust DNA loop-dependent bacterial promoter repression by covalent TALEDs and verify that DNA looping dramatically enhances promoter repression in E. coli. We characterize repression using a thermodynamic model that quantitates this favorable contribution of DNA looping. This analysis unequivocally and quantitatively demonstrates that optimized TALED proteins can drive loop-dependent promoter repression in E. coli comparable to the natural LacI repressor system. This work elucidates key design principles that set the stage for wide application of TALED-dependent DNA loop-based repression of target genes.

OptoLacI: optogenetically engineered lactose operon repressor LacI responsive to light instead of IPTG.

Optogenetics' advancement has made light induction attractive for controlling biological processes due to its advantages of fine-tunability, reversibility, and low toxicity. The lactose operon induction system, commonly used in Escherichia coli, relies on the binding of lactose or isopropyl ß-d-1-thiogalactopyranoside (IPTG) to the lactose repressor protein LacI, playing a pivotal role in controlling the lactose operon. Here, we harnessed the light-responsive light-oxygen-voltage 2 (LOV2) domain from Avena sativa phototropin 1 as a tool for light control and engineered LacI into two light-responsive variants, OptoLacIL and OptoLacID. These variants exhibit direct responsiveness to light and darkness, respectively, eliminating the need for IPTG. Building upon OptoLacI, we constructed two light-controlled E. coli gene expression systems, OptoE.coliLight system and OptoE.coliDark system. These systems enable bifunctional gene expression regulation in E. coli through light manipulation and show superior controllability compared to IPTG-induced systems. We applied the OptoE.coliDark system to protein production and metabolic flux control. Protein production levels are comparable to those induced by IPTG. Notably, the titers of dark-induced production of 1,3-propanediol (1,3-PDO) and ergothioneine exceeded 110% and 60% of those induced by IPTG, respectively. The development of OptoLacI will contribute to the advancement of the field of optogenetic protein engineering, holding substantial potential applications across various fields.

Genetic switching by the Lac repressor is based on two-state Monod-Wyman-Changeux allostery.

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod-Wyman-Changeux model.

DNA polymerase ε-dependent modulation of the pausing property of the CMG helicase at the barrier.

The proper pausing of replication forks at barriers on chromosomes is important for genome integrity. However, the detailed mechanism underlying this process has not been well elucidated. Here, we successfully reconstituted fork-pausing reactions from purified yeast proteins on templates that had binding sites for the LacI, LexA, and/or Fob1 proteins; the forks paused specifically at the protein-bound sites. Moreover, although the replicative helicase Cdc45-Mcm2-7-GINS (CMG) complex alone unwound the protein-bound templates, the unwinding of the LacI-bound site was impeded by the presence of a main leading strand DNA polymerase: polymerase ε (Polε). This suggests that Polε modulates CMG to pause at these sites.

Robustness of DNA looping across multiple cell divisions in individual bacteria.

DNA looping has emerged as a central paradigm of transcriptional regulation, as it is shared across many living systems. One core property of DNA looping-based regulation is its ability to greatly enhance repression or activation of genes with only a few copies of transcriptional regulators. However, this property based on a small number of proteins raises the question of the robustness of such a mechanism with respect to the large intracellular perturbations taking place during growth and division of the cell. Here we address the issue of sensitivity to variations of intracellular parameters of gene regulation by DNA looping. We use the lac system as a prototype to experimentally identify the key features of the robustness of DNA looping in growing Escherichia coli cells. Surprisingly, we observe time intervals of tight repression spanning across division events, which can sometimes exceed 10 generations. Remarkably, the distribution of such long time intervals exhibits memoryless statistics that is mostly insensitive to repressor concentration, cell division events, and the number of distinct loops accessible to the system. By contrast, gene regulation becomes highly sensitive to these perturbations when DNA looping is absent. Using stochastic simulations, we propose that the observed robustness to division emerges from the competition between fast, multiple rebinding events of repressors and slow initiation rate of the RNA polymerase. We argue that fast rebinding events are a direct consequence of DNA looping that ensures robust gene repression across a range of intracellular perturbations.

Creating Single-Copy Genetic Circuits.

Synthetic biology is increasingly used to develop sophisticated living devices for basic and applied research. Many of these genetic devices are engineered using multi-copy plasmids, but as the field progresses from proof-of-principle demonstrations to practical applications, it is important to develop single-copy synthetic modules that minimize consumption of cellular resources and can be stably maintained as genomic integrants. Here we use empirical design, mathematical modeling, and iterative construction and testing to build single-copy, bistable toggle switches with improved performance and reduced metabolic load that can be stably integrated into the host genome. Deterministic and stochastic models led us to focus on basal transcription to optimize circuit performance and helped to explain the resulting circuit robustness across a large range of component expression levels. The design parameters developed here provide important guidance for future efforts to convert functional multi-copy gene circuits into optimized single-copy circuits for practical, real-world use.

Cause and Consequence of Tethering a SubTAD to Different Nuclear Compartments.

Detailed genomic contact maps have revealed that chromosomes are structurally organized in megabase-sized topologically associated domains (TADs) that encompass smaller subTADs. These domains segregate in the nuclear space to form active and inactive nuclear compartments, but cause and consequence of compartmentalization are largely unknown. Here, we combined lacO/lacR binding platforms with allele-specific 4C technologies to track their precise position in the three-dimensional genome upon recruitment of NANOG, SUV39H1, or EZH2. We observed locked genomic loci resistant to spatial repositioning and unlocked loci that could be repositioned to different nuclear subcompartments with distinct chromatin signatures. Focal protein recruitment caused the entire subTAD, but not surrounding regions, to engage in new genomic contacts. Compartment switching was found uncoupled from transcription changes, and the enzymatic modification of histones per se was insufficient for repositioning. Collectively, this suggests that trans-associated factors influence three-dimensional compartmentalization independent of their cis effect on local chromatin composition and activity.

Positive supercoiling favors transcription elongation through lac repressor-mediated DNA loops.

Some proteins, like the lac repressor (LacI), mediate long-range loops that alter DNA topology and create torsional barriers. During transcription, RNA polymerase generates supercoiling that may facilitate passage through such barriers. We monitored E. coli RNA polymerase progress along templates in conditions that prevented, or favored, 400 bp LacI-mediated DNA looping. Tethered particle motion measurements revealed that RNA polymerase paused longer at unlooped LacI obstacles or those barring entry to a loop than those barring exit from the loop. Enhanced dissociation of a LacI roadblock by the positive supercoiling generated ahead of a transcribing RNA polymerase within a torsion-constrained DNA loop may be responsible for this reduction in pause time. In support of this idea, RNA polymerase transcribed 6-fold more slowly through looped DNA and paused at LacI obstacles for 66% less time on positively supercoiled compared to relaxed templates, especially under increased tension (torque). Positive supercoiling propagating ahead of polymerase facilitated elongation along topologically complex, protein-coated templates.

Deep representation learning improves prediction of LacI-mediated transcriptional repression.

Recent progress in DNA synthesis and sequencing technology has enabled systematic studies of protein function at a massive scale. We explore a deep mutational scanning study that measured the transcriptional repression function of 43,669 variants of the Escherichia coli LacI protein. We analyze structural and evolutionary aspects that relate to how the function of this protein is maintained, including an in-depth look at the C-terminal domain. We develop a deep neural network to predict transcriptional repression mediated by the lac repressor of Escherichia coli using experimental measurements of variant function. When measured across 10 separate training and validation splits using 5,009 single mutations of the lac repressor, our best-performing model achieved a median Pearson correlation of 0.79, exceeding any previous model. We demonstrate that deep representation learning approaches, first trained in an unsupervised manner across millions of diverse proteins, can be fine-tuned in a supervised fashion using lac repressor experimental datasets to more effectively predict a variant's effect on repression. These findings suggest a deep representation learning model may improve the prediction of other important properties of proteins.

Programmable gene regulation for metabolic engineering using decoy transcription factor binding sites.

Transcription factor decoy binding sites are short DNA sequences that can titrate a transcription factor away from its natural binding site, therefore regulating gene expression. In this study, we harness synthetic transcription factor decoy systems to regulate gene expression for metabolic pathways in Escherichia coli. We show that transcription factor decoys can effectively regulate expression of native and heterologous genes. Tunability of the decoy can be engineered via changes in copy number or modifications to the DNA decoy site sequence. Using arginine biosynthesis as a showcase, we observed a 16-fold increase in arginine production when we introduced the decoy system to steer metabolic flux towards increased arginine biosynthesis, with negligible growth differences compared to the wild type strain. The decoy-based production strain retains high genetic integrity; in contrast to a gene knock-out approach where mutations were common, we detected no mutations in the production system using the decoy-based strain. We further show that transcription factor decoys are amenable to multiplexed library screening by demonstrating enhanced tolerance to pinene with a combinatorial decoy library. Our study shows that transcription factor decoy binding sites are a powerful and compact tool for metabolic engineering.

Switching off: The phenotypic transition to the uninduced state of the lactose uptake pathway.

The lactose uptake pathway of E. coli is a paradigmatic example of multistability in gene regulatory circuits. In the induced state of the lac pathway, the genes comprising the lac operon are transcribed, leading to the production of proteins that import and metabolize lactose. In the uninduced state, a stable repressor-DNA loop frequently blocks the transcription of the lac genes. Transitions from one phenotypic state to the other are driven by fluctuations, which arise from the random timing of the binding of ligands and proteins. This stochasticity affects transcription and translation, and ultimately molecular copy numbers. Our aim is to understand the transition from the induced to the uninduced state of the lac operon. We use a detailed computational model to show that repressor-operator binding and unbinding, fluctuations in the total number of repressors, and inducer-repressor binding and unbinding all play a role in this transition. Based on the timescales on which these processes operate, we construct a minimal model of the transition to the uninduced state and compare the results with simulations and experimental observations. The induced state turns out to be very stable, with a transition rate to the uninduced state lower than 2×10-9 per minute. In contrast to the transition to the induced state, the transition to the uninduced state is well described in terms of a 2D diffusive system crossing a barrier, with the diffusion rates emerging from a model of repressor unbinding.

The genotype-phenotype landscape of an allosteric protein.

Allostery is a fundamental biophysical mechanism that underlies cellular sensing, signaling, and metabolism. Yet a quantitative understanding of allosteric genotype-phenotype relationships remains elusive. Here, we report the large-scale measurement of the genotype-phenotype landscape for an allosteric protein: the lac repressor from Escherichia coli, LacI. Using a method that combines long-read and short-read DNA sequencing, we quantitatively measure the dose-response curves for nearly 105 variants of the LacI genetic sensor. The resulting data provide a quantitative map of the effect of amino acid substitutions on LacI allostery and reveal systematic sequence-structure-function relationships. We find that in many cases, allosteric phenotypes can be quantitatively predicted with additive or neural-network models, but unpredictable changes also occur. For example, we were surprised to discover a new band-stop phenotype that challenges conventional models of allostery and that emerges from combinations of nearly silent amino acid substitutions.

T7-lac promoter vectors spontaneous derepression caused by plant-derived growth media may lead to serious expression problems: a systematic evaluation.

BACKGROUND: The widespread usage of protein expression systems in Escherichia coli (E. coli) is a workhorse of molecular biology research that has practical applications in biotechnology industry, including the production of pharmaceutical drugs. Various factors can strongly affect the successful construction and stable maintenance of clones and the resulting biosynthesis levels. These include an appropriate selection of recombinant hosts, expression systems, regulation of promoters, the repression level at an uninduced state, growth temperature, codon usage, codon context, mRNA secondary structure, translation kinetics, the presence/absence of chaperons and others. However, optimization of the growth medium's composition is often overlooked. We systematically evaluate this factor, which can have a dramatic effect on the expression of recombinant proteins, especially those which are toxic to a recombinant host. RESULTS: Commonly used animal tissue- and plant-based media were evaluated using a series of clones in pET vector, containing expressed Open Reading Frames (ORFs) with a wide spectrum of toxicity to the recombinant E. coli: (i) gfpuv (nontoxic); (ii) tp84_28-which codes for thermophilic endolysin (moderately toxic); and (iii) tthHB27IRM-which codes for thermophilic restriction endonuclease-methyltransferase (REase-MTase)-RM.TthHB27I (very toxic). The use of plant-derived peptones (soy peptone and malt extract) in a culture medium causes the T7-lac expression system to leak. We show that the presence of raffinose and stachyose (galactoside derivatives) in those peptones causes premature and uncontrolled induction of gene expression, which affects the course of the culture, the stability of clones and biosynthesis levels. CONCLUSIONS: The use of plant-derived peptones in a culture medium when using T7-lac hybrid promoter expression systems, such as Tabor-Studier, can lead to uncontrolled production of a recombinant protein. These conclusions also extend to other, lac operator-controlled promoters. In the case of proteins which are toxic to a recombinant host, this can result in mutations or deletions in the expression vector and/or cloned gene, the death of the host or highly decreased expression levels. This phenomenon is caused by the content of certain saccharides in plant peptones, some of which (galactosides) may act as T7-lac promoter inducer by interacting with a Lac repressor. Thus, when attempting to overexpress toxic proteins, it is recommended to either not use plant-derived media or to use them with caution and perform a pilot-scale evaluation of the derepression effect on a case-by-case basis.

Engineered signal-coupled inducible promoters: measuring the apparent RNA-polymerase resource budget.

Inducible promoters are a central regulatory component in synthetic biology, metabolic engineering, and protein production for laboratory and commercial uses. Many of these applications utilize two or more exogenous promoters, imposing a currently unquantifiable metabolic burden on the living system. Here, we engineered a collection of inducible promoters (regulated by LacI-based transcription factors) that maximize the free-state of endogenous RNA polymerase (RNAP). We leveraged this collection of inducible promotors to construct simple two-channel logical controls that enabled us to measure metabolic burden - as it relates to RNAP resource partitioning. The two-channel genetic circuits utilized sets of signal-coupled transcription factors that regulate cognate inducible promoters in a coordinated logical fashion. With this fundamental genetic architecture, we evaluated the performance of each inducible promoter as discrete operations, and as coupled systems to evaluate and quantify the effects of resource partitioning. Obtaining the ability to systematically and accurately measure the apparent RNA-polymerase resource budget will enable researchers to design more robust genetic circuits, with significantly higher fidelity. Moreover, this study presents a workflow that can be used to better understand how living systems adapt RNAP resources, via the complementary pairing of constitutive and regulated promoters that vary in strength.

Predictive shifts in free energy couple mutations to their phenotypic consequences.

Mutation is a critical mechanism by which evolution explores the functional landscape of proteins. Despite our ability to experimentally inflict mutations at will, it remains difficult to link sequence-level perturbations to systems-level responses. Here, we present a framework centered on measuring changes in the free energy of the system to link individual mutations in an allosteric transcriptional repressor to the parameters which govern its response. We find that the energetic effects of the mutations can be categorized into several classes which have characteristic curves as a function of the inducer concentration. We experimentally test these diagnostic predictions using the well-characterized LacI repressor of Escherichia coli, probing several mutations in the DNA binding and inducer binding domains. We find that the change in gene expression due to a point mutation can be captured by modifying only the model parameters that describe the respective domain of the wild-type protein. These parameters appear to be insulated, with mutations in the DNA binding domain altering only the DNA affinity and those in the inducer binding domain altering only the allosteric parameters. Changing these subsets of parameters tunes the free energy of the system in a way that is concordant with theoretical expectations. Finally, we show that the induction profiles and resulting free energies associated with pairwise double mutants can be predicted with quantitative accuracy given knowledge of the single mutants, providing an avenue for identifying and quantifying epistatic interactions.

DNA looping by protamine follows a nonuniform spatial distribution.

DNA looping plays an important role in cells in both regulating and protecting the genome. Often, studies of looping focus on looping by prokaryotic transcription factors like lac repressor or by structural maintenance of chromosomes proteins such as condensin. Here, however, we are interested in a different looping method whereby condensing agents (charge ≥+3) such as protamine proteins neutralize the DNA, causing it to form loops and toroids. We considered two previously proposed mechanisms for DNA looping by protamine. In the first mechanism, protamine stabilizes spontaneous DNA fluctuations, forming randomly distributed loops along the DNA. In the second mechanism, protamine binds and bends the DNA to form a loop, creating a distribution of loops that is biased by protamine binding. To differentiate between these mechanisms, we imaged both spontaneous and protamine-induced loops on short-length (≤1 µm) DNA fragments using atomic force microscopy. We then compared the spatial distribution of the loops to several model distributions. A random looping model, which describes the mechanism of spontaneous DNA folding, fit the distribution of spontaneous loops, but it did not fit the distribution of protamine-induced loops. Specifically, it failed to predict a peak in the spatial distribution of loops at an intermediate location along the DNA. An electrostatic multibinding model, which was created to mimic the bind-and-bend mechanism of protamine, was a better fit of the distribution of protamine-induced loops. In this model, multiple protamines bind to the DNA electrostatically within a particular region along the DNA to coordinate the formation of a loop. We speculate that these findings will impact our understanding of protamine's in vivo role for looping DNA into toroids and the mechanism of DNA condensation by condensing agents more broadly.

Evolution of a New Function by Fusion between Phage DNA and a Bacterial Gene.

Mobile genetic elements, such as plasmids, phages, and transposons, are important sources for evolution of novel functions. In this study, we performed a large-scale screening of metagenomic phage libraries for their ability to suppress temperature-sensitivity in Salmonella enterica serovar Typhimurium strain LT2 mutants to examine how phage DNA could confer evolutionary novelty to bacteria. We identified an insert encoding 23 amino acids from a phage that when fused with a bacterial DNA-binding repressor protein (LacI) resulted in the formation of a chimeric protein that localized to the outer membrane. This relocalization of the chimeric protein resulted in increased membrane vesicle formation and an associated suppression of the temperature sensitivity of the bacterium. Both the host LacI protein and the extracellular 23-amino acid stretch are necessary for the generation of the novel phenotype. Furthermore, mutational analysis of the chimeric protein showed that although the native repressor function of the LacI protein is maintained in this chimeric structure, it is not necessary for the new function. Thus, our study demonstrates how a gene fusion between foreign DNA and bacterial DNA can generate novelty without compromising the native function of a given gene.

Effect of Photocaged Isopropyl ß-d-1-thiogalactopyranoside Solubility on the Light Responsiveness of LacI-controlled Expression Systems in Different Bacteria.

Photolabile protecting groups play a significant role in controlling biological functions and cellular processes in living cells and tissues, as light offers high spatiotemporal control, is non-invasive as well as easily tuneable. In the recent past, photo-responsive inducer molecules such as 6-nitropiperonyl-caged IPTG (NP-cIPTG) have been used as optochemical tools for Lac repressor-controlled microbial expression systems. To further expand the applicability of the versatile optochemical on-switch, we have investigated whether the modulation of cIPTG water solubility can improve the light responsiveness of appropriate expression systems in bacteria. To this end, we developed two new cIPTG derivatives with different hydrophobicity and demonstrated both an easy applicability for the light-mediated control of gene expression and a simple transferability of this optochemical toolbox to the biotechnologically relevant bacteria Pseudomonas putida and Bacillus subtilis. Notably, the more water-soluble cIPTG derivative proved to be particularly suitable for light-mediated gene expression in these alternative expression hosts.