Previous Article |
Table of Contents
| Next Article
Genetics
Programmable cells: Interfacing natural and engineered gene networks



Department of Biomedical Engineering, Center for BioDynamics, and Center for Advanced Biotechnology, Boston University, 44 Cummington Street, Boston, MA 02215
Contributed by Charles R. Cantor, April 26, 2004
| Abstract |
|---|
|
|
|---|
Novel cellular behaviors and characteristics can be obtained by coupling engineered gene networks to the cell's natural regulatory circuitry through appropriately designed input and output interfaces. Here, we demonstrate how an engineered genetic circuit can be used to construct cells that respond to biological signals in a predetermined and programmable fashion. We employ a modular design strategy to create Escherichia coli strains where a genetic toggle switch is interfaced with: (i) the SOS signaling pathway responding to DNA damage, and (ii) a transgenic quorum sensing signaling pathway from Vibrio fischeri. The genetic toggle switch endows these strains with binary response dynamics and an epigenetic inheritance that supports a persistent phenotypic alteration in response to transient signals. These features are exploited to engineer cells that form biofilms in response to DNA-damaging agents and cells that activate protein synthesis when the cell population reaches a critical density. Our work represents a step toward the development of "plug-and-play" genetic circuitry that can be used to create cells with programmable behaviors.
heterologous gene expression | synthetic biology | Escherichia coli
Many cell regulatory systems are organized as modules (23-25) and a similar design strategy may allow the construction of cells with desired behaviors and characteristics. We envision that engineered gene networks can be used as regulatory modules and interfaced with the cell's genetic circuitry as "plug-and-play" devices to execute specific programs in response to particular biological signals. The simplest programmable cell obtained with this design strategy would be comprised of three distinct modules (Fig. 1): (i) a signaling pathway (the biosensor module) that detects relevant signals and interfaces these signals to a regulatory circuit, (ii) an artificial genetic module (the regulatory circuit) capable of responding to the signals transmitted by the biosensor module, and directing output signals according to its engineered properties, and (iii) an output interface that converts the signals transmitted by the regulatory circuit into a biological response. The behavior of the programmed cell is then determined by the dynamical and logical properties of the regulatory module and by the signaling pathways that are used as input and output interfaces.
|
|
| Experimental Procedures |
|---|
|
|
|---|
Fluorescence Measurements. GFP expression was quantified by using a FACSCalibur flow cytometer. Samples were prepared by pelleting cells from 1 ml of culture followed by resuspension in phosphate-buffered saline.
DNA Damage. Cells were grown aerobically in LB medium containing the appropriate antibiotics at 37°C and 300 rpm. Colonies were picked from selective plates and grown for 17-24 h, followed by an additional 16 h in medium containing 2 mM isopropyl-
-thiogalactopyranoside (IPTG). DNA damage was induced with mitomycin C (MMC) or UV irradiation. In the experiments with MMC treatment, the IPTG-containing culture was used to inoculate fresh LB medium with different MMC concentrations and grown for 15 h. The MMC-treated cells were grown for 3-56 h with dilutions every 12 h to keep the cells in the logarithmic growth phase. In the experiments with UV treatment, cells were plated and incubated for 2 h at 30°C before being exposed to irradiation (Stratalinker 2400) for 1-10 s. Cells were subsequently collected and grown in fresh medium for 4 h before being filtered (0.22-µm Millipore Millex-GV membrane filter) and assayed.
Biofilm Formation. Cells were grown aerobically in M63 minimal medium [1.052 g/liter KH2PO4/5.613 g/liter K2HPO4/2.0 g/liter (NH4)2SO4/0.50 mg FeSO4(H2O)7/1.0 mmol MgSO4, pH 7.2] containing 0.2% glucose and appropriate antibiotics at 37°C and 300 rpm. After exposure to MMC or UV irradiation, a small number of cells were used to inoculate 100-µl fresh M63 loaded into 96-well polystyrene plates. The plates were incubated for 24 h before the level of biofilm was quantified by using a crystal violet staining assay (47). Absorbance at 600 nm was measured by using a TECAN SPECTRAfluor Plus plate reader. Microfermentor experiments were carried out by using 20-ml continuous-flow fermentors (flow rate, 13 ml/h), stirred by aeration with sterile air and containing submerged Corning glass plates as the substratum for the biofilm. The fermentors were inoculated with 10 µl of culture treated with MMC as described above. Digital pictures were taken 48 h later.
AHL-Dependent Expression. All experiments involving the strains B1 and B2 were carried out in LB medium at 30°C unless otherwise stated. Cells were kept in the logarithmic growth phase by dilutions at appropriate intervals. AHL used to induce strain B1 [N-(
-ketocaproyl)-L-homoserine lactone] was obtained from Sigma. Cells with high and low initial GFP expression were obtained by growth in medium containing 2 mM IPTG for 12 h and growth at 42°C for 12 h, respectively. The cells were subsequently washed and used to inoculate fresh medium. The density-dependent expression experiment was carried out by growing the transformed cells on selective plates containing 2 mM IPTG, followed by growth at very low cell densities for 8 h in LB containing 2 mM IPTG. Cells were subsequently pelleted, washed three times, and used to inoculate batch cultures at various initial cell densities. The absorbance (cell density) of the cultures at 600 nm (A600) was determined with a SPECTRAfluor Plus plate reader.
| Results |
|---|
|
|
|---|
cI, that encode the transcriptional regulator proteins, LacR and
CI. The lacI gene is expressed from a modified PL promoter, PL*, which is repressed by
CI. The
cI gene is expressed from a promoter, Ptrc, which is repressed by LacR. This design endows cells with two distinct phenotypic states (4): one where the
CI activity is high and the expression of lacI is low, and one where the activity of LacR is high and the expression of
cI is low (Fig. 2A).
|
CI/low LacR state to the high LacR/low
CI state. When
CI activity is decreased, lacI expression is derepressed and LacR activity increases. This represses
cI expression, which decreases
CI activity and further increases LacR activity. The same result can be achieved with a perturbation that increases the activity of LacR. In both cases, a transition from one stable state to the other occurs if the perturbation is sufficiently large to bring the system across a certain threshold (see supporting information). Transitions from one stable state to the other can be induced by high-amplitude random fluctuations, referred to as noise-induced transitions (48), or by signals that temporarily change the parameters of the system. In bistable gene circuits, noise-induced transitions can cause individual cells to change expression state at random (49). The result is the emergence of a mixed population consisting of cells in different expression states, which appears as a bimodal population distribution when protein levels are measured in single cells (7, 13, 50).
The genetic toggle switch is a robust bistable system, and noise-induced transitions are rare (4). In such systems, transitions from one stable state to the other can be induced by a signal that temporarily brings the system out of the region of bistability. A mathematical analysis (see supporting information) indicates that transitions from the high
CI state to the high LacR state can be induced by a signal that temporarily increases (i) the
CI decay rate or (ii) the LacR basal synthesis rate. The simulated response of a single cell to such signals is shown in Fig. 2B. It illustrates how a cell initially in the high
CI state switches to the high LacR state as a result of a transient increase in
CI proteolysis. Increasing the basal LacR synthesis rate gives a similar response (see supporting information). In both cases, a transition to the high LacR state occurs when the signal reaches a threshold value where the high
CI state is destabilized. Because individual cells have slightly different threshold values, due, for instance, to variability in plasmid copy number, and because the probability of a noise-induced transition increases as the bifurcation parameter approaches the threshold value (48), it is expected that intermediate signals will give rise to bimodal population distributions.
Guided by the mathematical analysis, we interfaced the toggle switch with a natural signaling pathway that increases the rate of
CI decay and an engineered signaling pathway that increases the rate of LacR synthesis, respectively. The signaling pathway that degrades
CI in strains A1 and A2 (Table 1) is the SOS-response pathway, where the RecA coprotease is activated in the presence of single-stranded DNA (24). Activated RecA cleaves the
CI repressor protein, causing derepression of the PL promoter (51). The signaling pathway that increases the basal expression of the lacI gene in strains B1 and B2 (Table 1) is based on the quorum sensing pathway V. fischeri (29-31). In this pathway, the regulator protein of the lux operon, LuxR, is induced by AHL, and the induced LuxR protein activates expression from the lux promoter, PluxI. By placing the lacI gene downstream of PluxI, the rate of LacR synthesis is increased when AHL molecules are present in the environment.
Strain A1: Interfacing the SOS Pathway. Interfacing the genetic toggle switch (the regulatory circuit) with the SOS network (the biosensor module) required a series of alterations of the original pTAK plasmid (4). The toggle switch plasmid (pTSMa, see Fig. 3A) was made by replacing the cI857 gene, which encodes a
CI variant that is cleaved inefficiently by RecA (52), with wild-type
cI, and by changing the origin of replication to decrease the plasmid copy number. This was required to achieve compatibility between the biosensor module and the regulatory circuit (see below). As the output interface, we used a medium-copy number reporter plasmid (pCIRa), carrying a fusion of PL* and the gfp gene (Fig. 3A).
|
CI concentrations giving rise to variability in induction threshold.
The A1 strain can also detect brief exposures (<10 s) to UV irradiation (Fig. 3C). As in the experiment with MMC, UV irradiation at intermediate intensities induces a binary cellular response, resulting in bimodal population distributions. In both MMC- and UV-treated cells, the feedback architecture of the toggle switch module prevents expression of the
cI gene, even after the damaged DNA has been repaired and cells resume their pretreatment activities (see Fig. 2B). This allows cells to retain memory of DNA damage over successive generations, as demonstrated by the high expression state >48 h (corresponding to 50-60 generations) after the removal of MMC (see Fig. 3B).
Detecting DNA Damage with Strain A1. Strain A1 is a highly sensitive sensor of DNA damaging agents. Treatment with 1 ng/ml and 10 ng/ml MMC (Fig. 3B) gave a 1.9-fold and 19-fold increase in the population-averaged fluorescence signal (geometric mean), respectively. For comparison, the two sensor strains developed by Vollmer et al. (36) showed a 1.8-fold and 5.0-fold increase in the detected signal in response to 10 ng/ml MMC, whereas Kostrzynska et al. (37) reported a minimum detection limit of 4 ng/ml MMC (0.012 µM). In addition, the response of the A1 strain to UV irradiation at 6 J/m2 and 12 J/m2 was a 44-fold and 250-fold increase in average fluorescence (Fig. 3C). This represents a significant improvement in yield compared to previous reports of 4- to 5-fold increases in signal intensity at 10 J/m2 (37, 38).
To evaluate how the architecture of the regulatory circuit affects the ability of the A1 strain to detect DNA damage, we tested the response to MMC treatment of a strain that contains a regulatory circuit identical to the pTSMa toggle switch, except that it lacks the lacI feedback gene (plasmid pCIE). Fluorescence could not be detected after 15-h treatments at concentrations <1,000 ng/ml. A relatively weak fluorescence signal was detected when pCIE/pCIRa cells were assayed 30-60 min after the removal of MMC at concentrations between 1,000 and 4,000 ng/ml. The poor sensitivity and yield are probably due to the cellular activity of RecA being unable to cleave
CI at a sufficient rate (see supporting information for further discussion). However, GFP expression could not be detected in cells assayed 3 h after the removal of MMC. This indicates that the PL* promoter is active only for a limited time period after DNA damage in the circuit lacking the lacI gene. Comparing these results with those obtained from the A1 strain demonstrates that the feedback architecture of the genetic toggle switch provides at least a 1,000-fold improvement in sensitivity and enables readout of a detection event long after the DNA-damaging agent is removed. The latter could significantly improve the signal-to-noise ratio, because this feature allows for long signal integration. The disadvantages of a toggle switch-based biosensor include a loss of temporal information and a requirement of resetting, i.e., application of IPTG (4), between detection events.
Strain A2: Permanent Phenotypic Alteration. The above experiments indicate that the epigenetic inheritance capabilities of the genetic toggle switch might enable a permanent phenotypic change in response to a transient signal. To demonstrate this feature, we transferred the control of biofilm formation from the cell's natural circuitry to the genetic toggle switch in strain A2. This was done by deleting the traA gene (53) from the genome of the host strain and by constructing a biofilm-forming output plasmid (pBFR) where the expression of the traA gene is controlled by the PL* promoter. The engineered regulatory circuits of the A2 strain are illustrated in Fig. 4A. In this strain, the traA gene is constitutively expressed when the cells are in the high LacR/low
CI state. As a result, the strain is programmed to produce biofilm only when it has been subjected to DNA damage.
|
Strain B1: Interfacing General Input Signals. The experiments described above demonstrate that a natural signaling pathway can be interfaced with an engineered gene network. However, those studies exploit a preexisting molecular compatibility: the
CI protein is naturally cleaved upon the activation of the SOS pathway. As indicated by the mathematical analysis (supporting information), a transition between stable expression states in the genetic toggle switch can also occur if the expression of the repressed transcription factor protein is increased in response to an incoming signal. Thus, in principle, any cellular signal that activates the expression from a bacterial promoter might be used to couple the genetic toggle switch to natural regulatory circuits.
To demonstrate the generality of the input interface, and the plug-and-play features of the design strategy, we created a strain (B1) where a biosensor of AHL molecules interacts with the genetic toggle switch via the lacI gene. The engineered regulatory circuitry in the B1 strain (Fig. 5A) consists of a low-copy number AHL sensor plasmid (pAHLa), carrying a fusion of the lacI gene and the luxR-PluxI fragment from the V. fischeri lux operon, and a medium-copy number toggle switch plasmid (pTSMb1). In this strain, the toggle switch plasmid carries a copy of the gfp gene, such that cells fluoresce in the high
CI/low LacR state.
|
CI is the dominant repressor, but LacR still has some basal activity. The stability of the distinct expression states was confirmed in a separate control experiment where stable expression was observed for up to 50 h (corresponding to 50-60 generations) after the removal of the inducing factor (see supporting information).
To evaluate the switching dynamics and the sensitivity of the B1 strain, we conducted a series of experiments where cells initially in the high or low GFP expression states were exposed to AHL at various concentrations for 24 h (Fig. 5C). Regardless of the concentration of AHL, cells that were initially in the high LacR state (low GFP expression, open circles in Fig. 5C) remained in this state. Cells initially in the high
CI state (high GFP expression, closed circles in Fig. 5C) remained in that state at AHL concentrations <20 nM. All cells switched to the low GFP state when treated with AHL at 50 nM concentration or higher. Bimodal population distributions (Fig. 5C Inset) were observed at AHL concentrations between 20 and 50 nM. It is clear from Fig. 5 B and C that the long-term stability of the two expression states and the switching properties of the A1 strain (see Fig. 3B) are preserved in the B1 strain.
Strain B2: Density-Dependent Gene Activation. AHL is a natural biological signal secreted by Gram-negative bacteria as a means of coordinating cellular activity with the cell population density (29-31). To enable the E. coli population to measure its own density through AHL, we created the plasmid pAHLb where the luxI gene from V. fischeri is expressed polycistronically with the luxR gene and lacI is expressed from the PluxI promoter (Fig. 6A). The protein encoded by luxI is a synthetase that converts common precursor metabolites into AHL signaling molecules (29-31), and the extracellular concentration of AHL correlates with the cell density in cultures of cells that carry the luxI gene. As a result, LuxR should be activated, and lacI expression from the pAHLb plasmid increased, when the cell density increases.
|
Because of the modular design of this system, density-dependent synthesis of any protein can be achieved simply by replacing the gfp gene on the high-copy number reporter plasmid with a gene of interest. For example, programmed population control could be achieved by replacing gfp with a killer gene, as it was recently shown (34), by fusing the ccdB gene to the PluxI promoter and synthesizing LuxR and LuxI constitutively inside E. coli cells. Moreover, the sharp switching threshold of our system might be useful in industrial-scale production of proteins that inhibit cell growth because the target protein is synthesized only when the population has reached a high density.
| Discussion |
|---|
|
|
|---|
Our investigations also revealed some of the current challenges in constructing artificial gene circuits with sophisticated dynamical and computational properties. Interfacing these circuits with natural signaling pathways (or with each other) requires that the signals (e.g., activating or repressing transcription factors) are appropriately adjusted to allow effective information transmission between circuit modules while, at the same time, maintaining the proper function of the system as a whole. In many cases, the properties of the system must be optimized rather than those of the individual components (see supplemental information). In this respect, the modular design strategy could benefit significantly from the development of directed evolution technologies (54, 55) that can select for nontrivial dynamical behaviors. Moreover, more complex gene regulatory modules and interfaces need to be constructed to fully realize the capabilities of modular genetic control circuits. This could enable sophisticated processing capabilities, including event counting and signal integration.
| Acknowledgements |
|---|
| Footnotes |
|---|
-thiogalactopyranoside; MMC, mitomycin C; AHL, acylhomoserine lactone.
H.K. and M.K. contributed equally to this work. ![]()
To whom correspondence should be addressed. E-mail: jcollins{at}bu.edu.
© 2004 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
and Higher Organisms (Cell Press & Blackwell, Oxford).
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
E. L. Haseltine and F. H. Arnold Implications of Rewiring Bacterial Quorum Sensing Appl. Envir. Microbiol., January 15, 2008; 74(2): 437 - 445. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Brenner, D. K. Karig, R. Weiss, and F. H. Arnold Engineered bidirectional communication mediates a consensus in a microbial biofilm consortium PNAS, October 30, 2007; 104(44): 17300 - 17304. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Cai, B. Hartnett, C. Gustafsson, and J. Peccoud A syntactic model to design and verify synthetic genetic constructs derived from standard biological parts Bioinformatics, October 15, 2007; 23(20): 2760 - 2767. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Batt, B. Yordanov, R. Weiss, and C. Belta Robustness analysis and tuning of synthetic gene networks Bioinformatics, September 15, 2007; 23(18): 2415 - 2422. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. N. Win and C. D. Smolke From the Cover: A modular and extensible RNA-based gene-regulatory platform for engineering cellular function PNAS, September 4, 2007; 104(36): 14283 - 14288. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. F. Murphy, G. Balazsi, and J. J. Collins Combinatorial promoter design for engineering noisy gene expression PNAS, July 31, 2007; 104(31): 12726 - 12731. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Weber, M. Daoud-El Baba, and M. Fussenegger From the Cover: Synthetic ecosystems based on airborne inter- and intrakingdom communication PNAS, June 19, 2007; 104(25): 10435 - 10440. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Song, P. Smolen, E. Av-Ron, D. A. Baxter, and J. H. Byrne Dynamics of a Minimal Model of Interlocked Positive and Negative Feedback Loops of Transcriptional Regulation by cAMP-Response Element Binding Proteins Biophys. J., May 15, 2007; 92(10): 3407 - 3424. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Weber, J. Stelling, M. Rimann, B. Keller, M. Daoud-El Baba, C. C. Weber, D. Aubel, and M. Fussenegger A synthetic time-delay circuit in mammalian cells and mice PNAS, February 20, 2007; 104(8): 2643 - 2648. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Drubin, J. C. Way, and P. A. Silver Designing biological systems Genes & Dev., February 1, 2007; 21(3): 242 - 254. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Tian, S. Xu, J. Gao, and K. Burrage Simulated maximum likelihood method for estimating kinetic rates in gene expression Bioinformatics, January 1, 2007; 23(1): 84 - 91. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Burrage, L. Hood, and M. A. Ragan Advanced computing for systems biology Brief Bioinform, December 1, 2006; 7(4): 390 - 398. [Abstract] [Full Text] [PDF] |
||||
![]() |
|