Systematic approach to Escherichia coli cell population control using a genetic lysis circuit
© Hsu et al.; licensee BioMed Central Ltd. 2014
Published: 12 December 2014
Cell population control allows for the maintenance of a specific cell population density. In this study, we use lysis gene BBa_K117000 from the Registry of Standard Biological Parts, formed by MIT, to lyse Escherichia coli (E. coli). The lysis gene is regulated by a synthetic genetic lysis circuit, using an inducer-regulated promoter-RBS component. To make the design more easily, it is necessary to provide a systematic approach for a genetic lysis circuit to achieve control of cell population density.
Firstly, the lytic ability of the constructed genetic lysis circuit is described by the relationship between the promoter-RBS components and inducer concentration in a steady state model. Then, three types of promoter-RBS libraries are established. Finally, according to design specifications, a systematic design approach is proposed to provide synthetic biologists with a prescribed I/O response by selecting proper promoter-RBS component set in combination with suitable inducer concentrations, within a feasible range.
This study provides an important systematic design method for the development of next-generation synthetic gene circuits, from component library construction to genetic circuit assembly. In future, when libraries are more complete, more precise cell density control can be achieved.
Cell population control is a method of regulating cell density to maintain a self-sustaining population. Self-regulating mechanisms of cell population control have existed in nature for a long time. For example, the spore-forming bacterium Bacillus subtilis delays sporulation under nutrient-limited conditions by killing non-sporulating siblings and feeding on the dead cells to support spore formation . Colicin-producing bacteria produce bacteriocins to destroy nearby competitors through amino-acid starvation or DNA damage . Recently, a programmed control of cell population is proposed, using an endogenous genetic regulatory circuit . However, to make genetic lysis circuit design more easily, it is necessary to provide a systematic approach to achieve control of cell population density.
Naturally occurring lytic systems have the ability to trigger host cell lysis with specific proteins under certain circumstances, such as in the presence of antibiotics or competitors, or under conditions of amino-acid starvation. Some species of bacteria and bacteriophage have lytic protein that induces cell lysis. For example, expression of cloned T4 phage genes e or t can be used to the disrupt E. coli cells [4, 5]. Lysis by the T4 phage usually requires two gene products, transcribed by genes e and t, respectively. Gene e encodes a lysozyme, Gpe; gene t encodes a holin, Gpt. The expression of gene e weakens E. coli cell wall but does not lead to cell disruption, while the expression of gene t enables to induce cell lysis. Another example of a lytic protein is colicin E7 (ColE7), which is encoded by E. coli. Expression of ColE7 is regulated by an SOS response operon. The SOS response operon is composed of the ceaE7, ceiE7 and celE7 genes, the products of which are ColE7, ImE7 and LysE7, respectively . ColE7 can be neutralized through the formation of a protein complex with the immunity protein ImE7. The ColE7-ImE7 complex is exported from the colicin-producing cells by the lysis protein LysE7, which induces lysis of the host cell through break down of the cell membrane [7–12].
Lysis protein is able to activate the outer membrane phospholipase A (OMPLA), which allows colicin to cross the cell envelope and enter into the medium. Colicins are plasmid-encoded bacteriocins produced by E. coli that have antibiotic-like activity against closely related bacteria [12, 13]. Release of lysis protein can therefore control the population density of E. coli. In this study, lysis protein expression level in the genetic lysis circuit relies on tunable genetic components and inducer concentrations. The degree of cell lysis could therefore be fine-tuned by changing external inducer signals. Different inducer concentrations could result in different fates of the cells, for example severe death, modulated death, or slow growth. To make the design more easily, the lytic ability of the constructed genetic lysis circuit is described by the relationship between the promoter-RBS components and inducer concentrations in a steady state model. According to user-oriented specifications, the genetic lysis circuit can be constructed via selecting adequate promoter-RBS components in combination with a feasible range of inducer concentrations. In general, it requires long computation times when component libraries become large. Hence, a genetic algorithm (GA) search method is proposed to save time in evaluating and selecting promoter-RBS components.
This study provides an important systematic design method for the development of next-generation synthetic biology, from component library construction to genetic circuit assembly. When libraries are more complete, more precise cell density control for genetic lysis circuit can be achieved. The mechanism of cell lysis can be used to release useful macromolecules that cannot pass through cell membranes [4, 14]. For example, an operon encoding bcsA, bcsB, bcsC, and bcsD is required for bacterial cellulose synthesis (bcs) in Acetobacter xylinum . The expression of the operon can transform redundant glucose into cellulose to maintain intestinal peristalsis. In addition, an appropriate amount of cellulose has a proven role in preventing obesity. Cellulose, however, is a macromolecule that normally cannot pass through the cell membrane. A genetic lysis circuit can be combined with a bacterial cellulose synthesis system to promote cellulose release. In future, genetic lysis circuits may apply to development in, for example, drug discovery, metabolic control, and therapeutic treatment, with the help of the proposed design methodology.
The contributions of this paper are threefold: (1) Based on promoter-RBS kinetic strengths, we establish three kinds of promoter-RBS libraries. (2) Inducible promoter-RBS components are used to construct genetic lysis circuits with different lytic abilities for cell population control. (3) The proposed GA-based systematic searching method could provide synthetic biologists with a useful tool to design synthetic genetic lysis circuits.
Construction of genetic lysis circuit
Dynamic model of genetic lysis circuit
where x 1 and γ N denote the concentration of the lysis protein and the lysis rate of the lysis protein, respectively. Under normal growth conditions, cells do not suffer from the toxicity of the lysis protein. We used the endogenous gene (BBa_K117000) encoded in E. coli to control cell population density. The lysis protein can activate OMPLA to cause cell lysis and reduce cell density [12, 13]. From equation (2), it can be seen that the concentration of lysis protein x 1 is the negative regulator of cell population density. Because lysis ability is related to the concentration of lysis protein x 1 and the toxicity γ N of the lysis protein, we construct a dynamic model of the genetic lysis circuit to regulate the concentration of lysis protein x 1 , and thus lysis ability.
To construct the dynamic model, promoter-RBS regulation must be introduced. We first define the promoter-RBS regulation function P(P u , P 1 , TF, I), in which P u and P 1 denote the maximum and minimum promoter-RBS strengths, respectively, TF denotes transcription factor concentration, and I denotes inducer concentration. This function describes the biochemical aspect of the transcription and translation process. The details of the function (the mathematical model) can be found in Additional file 1. The dynamic model of the genetic lysis circuit with the repressor-regulated promoter-RBS component in Figure 1(a) thus is described as follows:
where (i, j) denote the i th constitutive promoter-RBS component and the j th repressor-regulated promoter-RBS component, respectively. Pc1 (P u,j ,0,0,0) denotes the promoter-RBS regulation activity of the first stage of the constitutive promoter-RBS components. P u,i denotes the maximum promoter-RBS strength of the i th constitutive promoter-RBS component. Pr2(Pu,j,Pi,j,xr1,I1) denotes the promoter-RBS regulation activity of the second stage of the repressor-regulated components. P u,j and P i,j denote the maximum and minimum promoter-RBS strengths of the j th repressor-regulated promoter-RBS component, respectively. I1 denotes the concentration, and xr1 denotes the concentration of the repressor. Kr2 and nr2 denote the binding affinity and binding cooperativity between the repressor and the corresponding promoter-RBS component in the second stage, respectively. KI1 denotes the dissociation rate between the inducer I1 and the repressor xr1. γ N is the lysis rate of the lysis protein. ω1(t) and ω2(t) denote cellular noise in the transcriptional and translational processes, respectively. Finally, ω3(t) denotes cellular noise in the cell population density.
From equation (3), the change in concentration of repressor protein x r1 is due to the difference between the protein generation rate Pc1(P u,j ,0,0,0) and the protein degradation rate γxr1, in combination with the dilution rate k. The repressor protein xr1 binds to the second stage of repressor-regulated components and reduces the regulation activity Pr2 (P u,j ,P i,j ,x r1 ,I1). The inducer I1 can remove the inhibitory effect of the repressor and increase the regulation activity Pr2 (Pu,j,Pi,j,xr1,I1). Similarly, the change in concentration of lysis protein xl2 is due to the difference between the protein generation rate Pr2 (Pu,j,Pi,j,xr1,I1) and the protein degradation rate γxl2, in combination with the dilution rate k. Cell population density is regulated by the concentration of lysis protein xl2. Lytic ability is therefore controlled by the regulation activity Pr2(Pu,j,Pi,j,xr1,I1), and the four regulated factors Pu,j, Pi,j, xr1, and I1 can be used to control cell population density by selecting appropriate promoter-RBS components and inducer concentrations.
where (i, m) denote the i th constitutive promoter-RBS component and the m th activator-regulated promoter-RBS component, respectively. I2 denotes the concentration of the inducer, xa1 denotes the concentration of the activator protein, and Ka2 and na2 denote the binding affinity and the binding cooperativity between activator and the corresponding promoter-RBS component in the second stage, respectively. Finally, KI2 denotes the dissociation rate between inducer I2 and activator xa1.
for the genetic lysis circuit with activator-regulated promoter RBS component in Figure 1(b). In these equations, the Gaussian noise parameter v i , i = 1, 2, with a zero mean and variance of σ i 2, denotes cellular noise for both the transcriptional and translational gene expression processes in the steady state. v3 denotes cellular noise in cell population density in the steady state. In (5) and (6), the steady state population density N SS changes depending on inducer concentration, constitutive promoter-RBS components, repressor-regulated components, and activator-regulated promoter-RBS components, i.e. i = 1,...,C, j = 1,...,R or i = 1,...,C, m = 1,...,A. In this study, we select values for these components from promoter-RBS libraries (see Additional file 1) to achieve a specified I/O reponse with the genetic lysis circuit. The details of the construction procedure for promoter-RBS component libraries can see in Additional file 1.
where ΔPu,i, ΔPu,j, ΔPl,j, ΔPu,m, ΔPl,m, Δγxr1 Δγxa1, Δγxl2, ΔγN, and Δk denote the standard deviations of the corresponding stochastic parameters, and n i (t), i = 1,2,3 denote Gaussian noise, have zero mean and unit variance, and account for random fluctuation sources. Thus, ΔPu,i, ΔPu,j, ΔPl,j, ΔPu,m, ΔPl,m, Δγxr1, Δγxa1, Δγxl2, ΔγN, and Δk denote the deterministic aspects of parameter variation, and ni(t), i = 1,2,3 denote the random fluctuation sources. The kinetic parameters in the steady state model in equations (5) and (6) are replaced by the parameter perturbations shown in (7) for robust design of the genetic circuit. These parameter fluctuations must be considered in the design procedure so that the synthetic genetic circuit can tolerate fluctuations in vivo .
Design specifications for the genetic lysis circuit
Desired I/O response N ref (I) of the genetic lysis circuit.
Well-characterized promoter-RBS component libraries and a feasible range of inducer concentrations.
Standard derivations of parameter fluctuations and environmental disturbances to be tolerated in vivo.
where S denotes the set of promoter-RBS components i and j (or i and m) selected from the corresponding component libraries, i.e. S = (i, j) for constitutive promoter-RBS components and repressor-regulated promoter-RBS components in (5), and S = (i, m) for constitutive promoter-RBS components and activator-regulated promoter-RBS components in (6). The inducer concentration I for I1 is S = (i, j) and I2 is S = (i, m).
Summary of design procedure for the genetic lysis circuit
Choose a species for cell population density control. In this study, we used E. coli.
Provide user-defined design specifications for the genetic lysis circuit.
Select an initial set S of promoter-RBS components and inducer concentrations.
Calculate the fitness function F(S, I) in equation (9) for each set S of promoter-RBS components and inducer concentrations.
Create an offspring set S using GA operators such as reproduction, crossover, and mutation.
Make copies of possible solutions, on basis of their fitness.
Swap values between two possible solutions.
Randomly alter the value in a possible solution.
Calculate the cost value of the new set S obtained by natural selection. Stop when the design goal is achieved or an acceptable solution is obtained. Otherwise, create the next generation and return to step 5.
Genetic lysis circuit design example for cell population density control in silico and verification via experiment in vivo
In this study, we focus on the design of synthetic genetic lysis circuit to achieve density control of cell population. A mathematical model is introduced to describe the dynamic and steady state regulatory behavior of the genetic lysis circuit. From the steady state mathematical models in (5) or (6), we find that if we want to control cell population density, we need to control lysis ability, which is related to transcriptional and translational processes. The promoter allows the RNA-polymerase molecules to latch onto a DNA strand and initialize the transcription of a downstream gene into mRNA, and the RBS allows the ribosome to bind and translate the mRNA. Hence, in this study, the promoter combined with the RBS is viewed as a genetic component for regulating lytic ability.
For convenient measurement of promoter-RBS components, we use GFP as a reporter to identify characterizations of promoter-RBS components. We can therefore use the parameters identified in Additional file 1 as component libraries to design the genetic lysis circuit. We also provide the design specifications for genetic lysis circuit for use by synthetic biologists. With a systematic approach, we can control the desired cell population density by selecting appropriate promoter-RBS component set in combination with a feasible range of inducer concentrations. In this study, for validation of the proposed design method, we perform some experiments with proper promoter-RBS components and inducer concentrations. The experimental results show that the simulation can robustly predict actual cell population density.
The precise control of cell population density with protein release is useful in medicine to cure disease, because changing medicine dosages may not achieve the desired effects. Additionally, evidence from systems biology indicates that apoptosis is involved in many pathways causing cell death [19, 20]. Each pathway depends on expression levels and can trigger cell death. For example, we embed a genetic lysis circuit into E. coli. Through the control of different expression levels of repressor-regulated or activator-regulated promoter-RBS components, a simple genetic lysis circuit can regulate cell death. This bottom-up design approach can potentially be extended to other complicated processes involving programmed cell death, which could help us to understand systematic phenomena that have always existed in nature.
The construction procedure for cell population density control is very important for engineering a more complex synthetic genetic lysis circuit. We characterize the genetic circuit and provide the desired cell population density beforehand. Using GA, we search for an appropriate set of promoter-RBS components in combinations with a feasible range of inducer concentrations to achieve the desired cell population density. The GA provides a useful tool in the construction of genetic lysis circuit for control of cell population density. In future, when libraries are more complete, more effective and efficient design methods for synthetic genetic lysis circuits may aid developments in drug discovery, metabolic control, and therapeutic treatment.
In this study, we engineer a genetic lysis circuit to control cell population density. Inducible promoter-RBS components are used to construct genetic lysis circuits with different lytic abilities for cell population control. Moreover, we provide four design specifications for cell population density control, allowing designers to select proper promoter-RBS components from libraries for the synthetic genetic lysis circuits. The design problem can be transformed by selecting proper promoter-RBS components in combination with a feasible range of inducer concentrations to achieve a desired I/O response. The proposed GA-based systematic searching methodology could provide synthetic biologists with a useful tool to design synthetic genetic lysis circuits. From the experimental results, we find that the data are close to the prescribed cell densities. Therefore, we believe that the proposed user-oriented design method for cell population density control will provide a useful guide in the rapidly growing field of synthetic biology.
The work was supported by the National Science Council of Taiwan under grant NSC 102-2745-E-007-001-ASP and NSC 103-2745-E-007-001-ASP.
The full funding for the publication fee came from National Science Council of Taiwan.
This article has been published as part of BMC systems Biology Volume 8 Supplement 5, 2014: Proceedings of the 25th International Conference on Genome Informatics (GIW/ISCB-Asia): Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/8/S5.
- Ellermeier CD, Hobbs EC, Gonzalez-Pastor JE, Losick R: A three-protein signaling pathway governing immunity to a bacterial cannibalism toxin. Cell. 2006, 124 (3): 549-559. 10.1016/j.cell.2005.11.041.View ArticlePubMedGoogle Scholar
- Gardner A, West SA, Buckling A: Bacteriocins, spite and virulence. Proceedings of the Royal Society of London Series B: Biological Sciences. 2004, 271 (1547): 1529-1535. 10.1098/rspb.2004.2756.PubMed CentralView ArticlePubMedGoogle Scholar
- You L, Cox RS, Weiss R, Arnold FH: Programmed population control by cell-cell communication and regulated killing. Nature. 2004, 428 (6985): 868-871. 10.1038/nature02491.View ArticlePubMedGoogle Scholar
- Morita M, Asami K, Tanji Y, Unno H: Programmed Escherichia coli cell lysis by expression of cloned T4 phage lysis genes. Biotechnology progress. 2001, 17 (3): 573-576. 10.1021/bp010018t.View ArticlePubMedGoogle Scholar
- Pasotti L, Zucca S, Lupotto M, Cusella De Angelis M, Magni P: Characterization of a synthetic bacterial self-destruction device for programmed cell death and for recombinant proteins release. Journal of biological engineering. 2011, 5 (8):Google Scholar
- Chak KF, Kuo WS, James R: Cloning and characterization of the ColE7 plasmid. Journal of general microbiology. 1991, 137 (1): 91-100. 10.1099/00221287-137-1-91.View ArticlePubMedGoogle Scholar
- Cavard D, Baty D, Howard S, Verheij H, Lazdunski C: Lipoprotein nature of the colicin A lysis protein: effect of amino acid substitutions at the site of modification and processing. Journal of bacteriology. 1987, 169 (5): 2187-2194.PubMed CentralPubMedGoogle Scholar
- CHAK K-F, JAMES R: Localization and characterization of a gene on the ColE3-CA38 plasmid that confers immunity to colicin E8. Journal of general microbiology. 1984, 130 (3): 701-710.PubMedGoogle Scholar
- Inukai M, Ghrayeb J, Nakamura K, Inouye M: Apolipoprotein, an intermediate in the processing of the major lipoprotein of the Escherichia coli outer membrane. Journal of Biological Chemistry. 1984, 259 (2): 757-760.PubMedGoogle Scholar
- PUGSLEY AP, COLE ST: An unmodified form of the ColE2 lysis protein, an envelope lipoprotein, retains reduced ability to promote colicin E2 release and lysis of producing cells. Journal of general microbiology. 1987, 133 (9): 2411-2420.PubMedGoogle Scholar
- Cavard D: Synthesis and functioning of the colicin E1 lysis protein: comparison with the colicin A lysis protein. Journal of bacteriology. 1991, 173 (1): 191-196.PubMed CentralPubMedGoogle Scholar
- Ellis T, Wang X, Collins JJ: Diversity-based, model-guided construction of synthetic gene networks with predicted functions. Nature biotechnology. 2009, 27 (5): 465-471. 10.1038/nbt.1536.PubMed CentralView ArticlePubMedGoogle Scholar
- Cascales E, Buchanan SK, Duché D, Kleanthous C, Lloubes R, Postle K, Riley M, Slatin S, Cavard D: Colicin biology. Microbiology and Molecular Biology Reviews. 2007, 71 (1): 158-229. 10.1128/MMBR.00036-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Lo TM, Tan MH, Hwang IY, Chang MW: Designing a synthetic genetic circuit that enables cell density-dependent auto-regulatory lysis for macromolecule release. Chemical Engineering Science. 2013Google Scholar
- Ross P, Mayer R, Benziman M: Cellulose biosynthesis and function in bacteria. Microbiological Reviews. 1991, 55 (1): 35-58.PubMed CentralPubMedGoogle Scholar
- Alon U: An introduction to systems biology: design principles of biological circuits. 2007, 10: Chapman & Hall/CRCGoogle Scholar
- Chen BS, Chang CH, Wang YC, Wu CH, Lee HC: Robust model matching design methodology for a stochastic synthetic gene network. Mathematical Biosciences. 2011, 230 (1): 23-36. 10.1016/j.mbs.2010.12.007.View ArticlePubMedGoogle Scholar
- Saeidi N, Wong CK, Lo TM, Nguyen HX, Ling H, Leong SSJ, Poh CL, Chang MW: Engineering microbes to sense and eradicate Pseudomonas aeruginosa, a human pathogen. Molecular systems biology. 2011, 7 (1):Google Scholar
- Elmore S: Apoptosis: a review of programmed cell death. Toxicologic pathology. 2007, 35 (4): 495-516. 10.1080/01926230701320337.PubMed CentralView ArticlePubMedGoogle Scholar
- Green DR: Means to an end: apoptosis and other cell death mechanisms. 2011, Cold Spring Harbor Laboratory Press Cold Spring Harbor, NYGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.