<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="3.0" xml:lang="EN"><front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><!--===== Grouping journal title elements =====--><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">08-PONE-RA-03523R1</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0002372</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline"><subject>Biotechnology</subject><subject>Molecular Biology</subject><subject>Computational Biology/Synthetic Biology</subject><subject>Computational Biology/Transcriptional Regulation</subject></subj-group></article-categories><title-group><article-title>Bimodal and Hysteretic Expression in Mammalian Cells from a Synthetic Gene Circuit</article-title><alt-title alt-title-type="running-head">Synthetic Gene Circuits</alt-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>May</surname><given-names>Tobias</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Eccleston</surname><given-names>Lee</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Herrmann</surname><given-names>Sabrina</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hauser</surname><given-names>Hansjörg</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Goncalves</surname><given-names>Jorge</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wirth</surname><given-names>Dagmar</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Gene Regulation and Differentiation, Helmholtz Centre for Infection Research, Braunschweig, Germany</addr-line>       </aff><aff id="aff2"><label>2</label><addr-line>Department of Engineering, University of Cambridge, Cambridge, United Kingdom</addr-line>       </aff><contrib-group><contrib contrib-type="editor" xlink:type="simple"><name name-style="western"><surname>Califano</surname><given-names>Andrea</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"/></contrib></contrib-group><aff id="edit1">Columbia University, United States of America</aff><author-notes><corresp id="cor1">* E-mail: <email xlink:type="simple">dagmar.wirth@helmholtz-hzi.de</email></corresp><fn fn-type="con"><p>Conceived and designed the experiments: DW TM LE JG. Performed the experiments: TM LE SH. Analyzed the data: DW TM HH LE JG. Wrote the paper: DW TM HH LE JG.</p></fn><fn fn-type="conflict"><p>The authors have declared that no competing interests exist.</p></fn></author-notes><pub-date pub-type="collection"><year>2008</year></pub-date><pub-date pub-type="epub"><day>4</day><month>6</month><year>2008</year></pub-date><volume>3</volume><issue>6</issue><elocation-id>e2372</elocation-id><history><date date-type="received"><day>5</day><month>2</month><year>2008</year></date><date date-type="accepted"><day>7</day><month>4</month><year>2008</year></date></history><!--===== Grouping copyright info into permissions =====--><permissions><copyright-year>2008</copyright-year><copyright-holder>May et al</copyright-holder><license><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><abstract><p>In order to establish cells and organisms with predictable properties, synthetic biology makes use of controllable, synthetic genetic devices. These devices are used to replace or to interfere with natural pathways. Alternatively, they may be interlinked with endogenous pathways to create artificial networks of higher complexity. While these approaches have been already successful in prokaryotes and lower eukaryotes, the implementation of such synthetic cassettes in mammalian systems and even animals is still a major obstacle. This is mainly due to the lack of methods that reliably and efficiently transduce synthetic modules without compromising their regulation properties.</p><p>To pave the way for implementation of synthetic regulation modules in mammalian systems we utilized lentiviral transduction of synthetic modules. A synthetic positive feedback loop, based on the Tetracycline regulation system was implemented in a lentiviral vector system and stably integrated in mammalian cells. This gene regulation circuit yields a bimodal expression response. Based on experimental data a mathematical model based on stochasticity was developed which matched and described the experimental findings. Modelling predicted a hysteretic expression responsewhich was verified experimentally. Thereby supporting the idea that the system is driven by stochasticity.</p><p>The results presented here highlight that the combination of three independent tools/methodologies facilitate the reliable installation of synthetic gene circuits with predictable expression characteristics in mammalian cells and organisms.</p></abstract><funding-group><funding-statement>Lee Eccleston would like to thank The MathWorks for their funding through The Math-Works Scholarship. Jorge Goncalves would like to thank the Nuffield Foundation and EPSRC for their support and funding. The work of Tobias May, Hansjörg Hauser and Dagmar Wirth was supported by the Bundesministerium für Bildung und Forschung (BMBF FKZ 0313940), the DFG (WI2648) and the Helmholtz Association (SBCancer). Further support was obtained from EU (Clinigene and FlpFlex (513769)).</funding-statement></funding-group><counts><page-count count="7"/></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>To date, most of the synthetic gene networks have been created for bacteria or for lower eukaryotes <xref ref-type="bibr" rid="pone.0002372-Becskei1">[1]</xref>–<xref ref-type="bibr" rid="pone.0002372-Becskei2">[7]</xref>. They were successfully employed to investigate different regulatory networks within these biological systems. A wealth of information could be provided by these studies. To facilitate systems biology based approaches for the elucidation of complex regulatory networks involved in immunology, developmental biology or infection, the development of appropriate synthetic modules as well as their establishment in mammalian cells and animals is indispensable. A thorough understanding and mathematical description of these modules in complex organisms would further pave the way to exploit such circuits within systems biology and also for applications in gene or cell therapy <xref ref-type="bibr" rid="pone.0002372-Weber1">[8]</xref>.</p><p>Pioneering studies demonstrated that the creation of synthetic gene circuits is feasible also in mammalian systems. The designed devices encompassed a toggle switch <xref ref-type="bibr" rid="pone.0002372-Kramer1">[9]</xref>, a hysteretic switch <xref ref-type="bibr" rid="pone.0002372-Kramer2">[10]</xref> and a time-delay circuit <xref ref-type="bibr" rid="pone.0002372-Weber2">[11]</xref>. A major limitation is that the manipulation of mammalian cells is tedious and time-consuming. Especially stable integration of synthetic gene circuits remains a challenge and is up to now only achieved in immortalized cell lines which usually do not or only partially reflect relevant properties of the cells they have been derived of. This can be partially attributed to the genetic changes which are required to provide unlimited cell growth. To implement synthetic cassettes in primary cells novel tools have to be developed.</p><p>Viruses have evolved strategies to efficiently transduce their genetic information to mammalian cells. In particular, adenoviruses, adeno-associated viruses and retroviruses were exploited for their potential to genetically modify mammalian cells <xref ref-type="bibr" rid="pone.0002372-Kay1">[12]</xref>, For this purpose, non-replicating viral vectors have been developed that can be transduced with the help of packaging cells <xref ref-type="bibr" rid="pone.0002372-Loew1">[e.g. 13]</xref>, <xref ref-type="bibr" rid="pone.0002372-Schucht1">[14]</xref>. While non-replicating extrachromosomal vectors (e.g. derived from adenoviruses) remain in an extrachromosomal state they lead to a transient expression in proliferating cells. Retroviral and especially lentiviral vectors are the optimal tools to provide stable integration of expression cassettes into the host's genome and thus guarantee the maintenance of the cassettes in replicating cells. Pseudotyping of lentiviral vectors allow to to transduce a broad range of different cell types from different species. Importantly, they are also capable of infecting slowly or non-proliferating cells and thus are preferential tools for many primary cells including stem cells.</p><p>Gene expression patterns in eukaryotes can be classified into two distinct types: graded and bimodal <xref ref-type="bibr" rid="pone.0002372-Fiering1">[15]</xref>. In the graded pattern, cells respond uniformly to increasing levels of the stimulant. In contrast, in the bimodal pattern, the level of the stimulant affects only the probability that a gene will be expressed in a given cell; this expression is always “maximal”. Bistable expression of metabolic genes such as those for galactose utilization in yeast <xref ref-type="bibr" rid="pone.0002372-Acar1">[16]</xref> and lactose utilization in E.coli <xref ref-type="bibr" rid="pone.0002372-Ozbudak1">[17]</xref> have been shown to result from positive feedback. One prominent example for bistability in eukaryotes is the cell-cycle oscillator, a biochemical system that has been shown to promote irreversible transitions between distinct mitotic and interphase states <xref ref-type="bibr" rid="pone.0002372-Pomerening1">[18]</xref>–<xref ref-type="bibr" rid="pone.0002372-Sha1">[20]</xref>. But also in the immune system the T helper activation is governed by an bimodal expression of the transcription factor NFATc2 which integrates graded T cell receptor activation into a bimodal IL2 expression <xref ref-type="bibr" rid="pone.0002372-Serfling1">[21]</xref>, <xref ref-type="bibr" rid="pone.0002372-Podtschaske1">[22]</xref>.</p><p>The bimodal expression mode may lead to multimodal variation in the phenotype of an isogenic cellular culture and in organisms (i.e. differentiation). In genetic networks this confers a memory to past events or environmental changes. Thus, stochasticity is an important feature of many biological processes in mammalian cells. Several means have been described to achieve a bimodal expression pattern. Amongst them are positive feedback loops and mutual repression <xref ref-type="bibr" rid="pone.0002372-Ferrell1">[23]</xref>. As an example, the latter can be accomplished using two interlinked negative feedback loops which create a toggle switch <xref ref-type="bibr" rid="pone.0002372-Gardner1">[4]</xref>, <xref ref-type="bibr" rid="pone.0002372-Kramer1">[9]</xref>. This toggle switch setup was employed in another study as a memory subsystem. This was implemented and used for recording the outputs of a larger DNA-damage sensor network <xref ref-type="bibr" rid="pone.0002372-Kobayashi1">[24]</xref>. Bimodal expression was further achieved by employing two opposing transcriptional regulators <xref ref-type="bibr" rid="pone.0002372-Rossi1">[25]</xref>. Becskei et al. (2001) <xref ref-type="bibr" rid="pone.0002372-Becskei3">[26]</xref> finally proved that in yeast even a less complex network based on the tet-system can result in a bimodal expression pattern upon gradual increase of the inducer tetracycline.</p><p>The aim of this study was to provide and evaluate a strategy to efficiently implement a bimodal synthetic gene circuit in mammalian cells. For this purpose, three different methodologies were exploited. A synthetic positive feedback system resulting in binary expression response was transduced via lentiviral transduction. This combination ensures that the synthetic gene circuit can be installed a broad range of mammalian cells including primary cells. After stable integration into the host cell the synthetic gene circuit displayed a bimodal expression pattern and the resulting experimental data were used to construct a mathematical model. The model followed closely the observed in vivo data and predicted a strong hysteretic response of the positive feedback loop. The in silico data were experimentally verified demonstrating the value of this concerted approach.</p></sec><sec id="s2"><title>Results</title><sec id="s2a"><title>Bimodal expression from a lentiviral transduced positive feedback loop</title><p>Efficient transduction and long-term stable expression of transgenes is achieved with the help of lentiviral transfer vectors both for basic research and gene therapeutic approaches <xref ref-type="bibr" rid="pone.0002372-Wiznerowicz1">[27]</xref>. Several generations of lentiviral vectors have been developed and allowed to successively reduce viral sequences <xref ref-type="bibr" rid="pone.0002372-Delenda1">[28]</xref>. These modifications led to higher safety. Importantly, through the deletion of the viral cis-acting regulatory elements (viral promoter, enhancer) the expression pattern is determined by the nature of the transduced expression cassette.</p><p>We chose a 3<sup>rd</sup> generation lentiviral vector devoid of any viral promoter elements for transduction of the synthetic expression module. This excluded any interference in between the viral vector backbone and the synthetic gene regulation circuit. A transcriptionally regulated gene circuit based on the Tetracycline(tet)-system was implemented in the lentiviral vector. This synthetic module was designed to generate a positive feedback response. Initiation of the positive feedback loop is achieved by a low basal level of the transactivator rtTA (<xref ref-type="fig" rid="pone-0002372-g001">Fig. 1A</xref>). This configuration meets the requirements for a synthetic gene circuit as its regulation is simple and homogenously achievable without any adverse/side effects mediated by the inducer. Secondly, it facilitates rapid induction and repression kinetics as well as multiple cycles of induction and repression. To precisely follow the induction/repression properties of this gene circuit a destabilized enhanced green fluorescent protein (eGFP) mutant <xref ref-type="bibr" rid="pone.0002372-Li1">[29]</xref> is coregulated, thereby reflecting the transcription of the transactivator. eGFP can be monitored in viable cells and facilitates the monitoring of expression on single cell level e.g. by flow cytometry.</p><fig id="pone-0002372-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002372.g001</object-id><label>Figure 1</label><caption><title>Setup and expression pattern of the positive feedback loop.</title><p>(a) Vector map of pJSARGFP transducing the positive feedback loop. (b) NIH3T3 cells were stably transduced through lentiviral infection with pJSARGFP. A representative clone was cultivated for four days with the indicated Dox concentrations (in µg/ml) and analyzed for GFP expression using flow cytometry. The numbers in brackets represent the percentage of cells which fall in the respective marker regions. The bimodal response is clearly visible from this data.</p></caption><graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0002372.g001" xlink:type="simple"/></fig><fig id="pone-0002372-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002372.g002</object-id><label>Figure 2</label><caption><title>Reaction scheme for the feedback system.</title><p>Each node corresponds to a reaction governed by a transition propensity, with lines and arrows representing reactants and products, respectively. A dotted line signifies that the ‘reactant’ is not consumed.</p></caption><graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0002372.g002" xlink:type="simple"/></fig><fig id="pone-0002372-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002372.g003</object-id><label>Figure 3</label><caption><title>Simulations of the expression behaviour of the positive feedback loop.</title><p>The simulations based on the Gillespie algorithm were run as described in the text. Increasing concentrations of Dox (0, 0.5, 1, 1.5 and 2.5 µg/mL) were applied as indicated. The resulting transactivator molecule counts were monitored. The population histograms show the probability of a respective transactivator level. The system parameters were k<sub>on</sub> = 20, k<sub>off</sub> = 20, sA = 50, sR = 5, sP = 0.2, δ<sub>M</sub> = 0.1, δ<sub>P</sub> = 0.05, k+ = 20 and k− = 10.</p></caption><graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0002372.g003" xlink:type="simple"/></fig><fig id="pone-0002372-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002372.g004</object-id><label>Figure 4</label><caption><title>Computer simulations predict hysteretic behaviour.</title><p>Based on the simulations in <xref ref-type="fig" rid="pone-0002372-g003">Figure 3</xref> and corresponding simulations of the on→off transitions. The lower curve represents the percentage of cells in the ‘HIGH’ expression region as the concentration of Dox is gradually increased. The upper curve shows the percentage of cells in the ‘HIGH’ expression region after a transition from a saturating Dox level to that shown.</p></caption><graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0002372.g004" xlink:type="simple"/></fig><fig id="pone-0002372-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002372.g005</object-id><label>Figure 5</label><caption><title>Hysteretic behaviour of cells carrying the positive feedback loop.</title><p>(a) Cells carrying the autoregulated cassette (<xref ref-type="fig" rid="pone-0002372-g001">Fig. 1A</xref>) were cultivated with 1 µg/mL of Dox for 4 days and then sorted. The eGFP positive cell population is designated with ‘high’, the eGFP negative cell population is designated with ‘low’. After sorting the cells were cultivated for additional 4 days with the indicated Dox concentrations. (b) The number of GFP positive cells is plotted against the applied Dox concentration. The blue line reflects the off→on experiment, the red line the corresponding on→off experiment.</p></caption><graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0002372.g005" xlink:type="simple"/></fig><p>We implemented the autoregulatory expression module in NIH3T3 cells using lentiviral gene transfer. Cell clones were established and expression was monitored on the single cell level by flow cytometry (Data from a representative clone is shown.). Increasing concentrations of the inducer were applied to investigate the pattern of expression response (<xref ref-type="fig" rid="pone-0002372-g001">Fig. 1B</xref>). Saturating inducer concentrations (2.5 µg/mL Doxycycline (Dox)) caused homogeneous activation and all cells showed eGFP expression. At non-saturating concentrations of the inducer (0.5–1.5 µg/mL) the population was split and yielded a binary population distribution. By varying the inducer concentration, the ratio of expressing to non-expressing cells could be adjusted. This is a specific property of the positive feedback loop since constant high level expression of the transactivator yielded a graded transcriptional response (data not shown and <xref ref-type="bibr" rid="pone.0002372-Kringstein1">[30]</xref>). Further, constitutive co-expression of the transactivator in these cells turned the response to a graded mode [data not shown].</p></sec><sec id="s2b"><title>Mathematical modelling of the binary expression module</title><p>A mathematical description of the synthetic module was produced. Standard stochastic techniques were employed to try and recreate the binary response <xref ref-type="bibr" rid="pone.0002372-McAdams1">[31]</xref>, <xref ref-type="bibr" rid="pone.0002372-Gillespie1">[32]</xref>. For this purpose, the particular regulatory network was modelled as a discrete, multi-valued, continuous time Markov process. This mathematical simulation was based on the reaction scheme proposed in <xref ref-type="fig" rid="pone-0002372-g002">Fig. 2</xref>. For the sake of simplicity, the eGFP loop that is simultaneously activated and that mirrors tTA expression was excluded. In this reaction scheme, the ‘allowable’ events are synthesis and degradation of mRNA (at two rates, ON and OFF, respectively), synthesis and degradation of proteins, the binding of Dox to the rtTA to form the activator complex, gene switching upon binding of the activator complex to the promoter and a term which reflects competition for Dox within the cellular population. The system was simulated using the Gillespie algorithm <xref ref-type="bibr" rid="pone.0002372-Gillespie1">[32]</xref>. Since exact values for the state transition propensities are not available, we used typical values from Kaern et al. (2005)<xref ref-type="bibr" rid="pone.0002372-Kaern1">[33]</xref>, which were then uniformly scaled to calibrate with the expected time scales.</p><p>The simulations were run using different initial concentrations of the inducer Dox (off→on simulation). <xref ref-type="fig" rid="pone-0002372-g003">Fig. 3</xref> summarises the results obtained from 20,000 independent simulations per given Dox concentration. These results demonstrate two distinct steady state levels of the transactivator, a ‘HIGH’ and a ‘LOW’ fraction. With increasing concentration of Dox the probability for the ‘HIGH’ fraction increases. Thus, the simulation reflects the bimodal distribution across the population as observed in the <italic>in vivo</italic> setting.</p></sec><sec id="s2c"><title>Hysteretic response of the positive-feedback</title><p>Strong positive feedback loops are often associated with hysteresis <xref ref-type="bibr" rid="pone.0002372-Xiong1">[34]</xref>. We tested if the mathematical algorithms as applied here would predict such behaviour; for this purpose, the simulations were run with inducer concentrations of 0, 0.5, 1.0, 1.5 and 2.5 µg/mL, but applying initial condition vectors that were taken from the end of the earlier simulations at saturating inducer levels (on→off simulation). The mathematical simulation showed that for the on→off situation a much lower concentration of the inducer was needed to switch off the system and decrease the probability of a given cell being in the ‘HIGH’ fraction. A summary of the on→off simulations is provided in <xref ref-type="fig" rid="pone-0002372-g004">Figure 4</xref>; for comparison, the corresponding graph for the off→on simulation is also indicated. Thus, the history of the system clearly defines the distribution within the population: for example, while a shift from 0 to 1 µg/mL of Dox gives rise to a clear binary pattern with 70% of cells in the ‘HIGH’ fraction, the downward shift from 2.5 µg/mL to 1 µg/mL results in more than 80% of the cells remaining ‘HIGH’. Importantly, for all applied concentrations of Dox, the final distribution favoured the ‘HIGH’ expression level for the on→off simulation.</p><p>We asked if this mathematical prediction is supported by experimental data and transferred the conditions used for simulation to isogenic cells with different cultivation histories. For this purpose, cells were cultivated for four days with non-saturating inducer concentrations (1 µg/ml). This treatment generates two populations (see <xref ref-type="fig" rid="pone-0002372-g001">Fig. 1B</xref>). Both cell populations were sorted and subsequently again cultivated with different inducer concentrations (0, 0,1, 0,5, 0,75, 1,0, 1,5, 2,0 and 2,5 µg/mL Dox). As expected, both populations behaved uniformly for 0 Dox and 2,5 µg/mL Dox: both ceased expression of eGFP upon cultivation without any Dox, and showed fully activated expression under saturating conditions (2.5 µg/mL). However, the cell pools displayed an inhomogeneous expression pattern at non-saturating levels of 0.1–2.0 µg/mL: while cells from the ‘low’ population provided a bimodal distribution, the ‘high’ population cells remained fully induced at inducer levels of 0.5 µg/mL and above (<xref ref-type="fig" rid="pone-0002372-g005">Fig. 5A</xref>). Only cells, cultivated with low concentrations of the inducer (0.1 µg/mL and below) also showed a bimodal expression pattern resulting in two populations. In (<xref ref-type="fig" rid="pone-0002372-g005">Fig. 5B</xref>) the obtained expression are plotted against the Dox concentration. Together, these data confirm the hysteretic behaviour of the cells as predicted from the mathematical model.</p></sec></sec><sec id="s3"><title>Discussion</title><p>A goal of synthetic biology is to describe and understand naturally occurring networks by constructing synthetic analogues. Mathematical modelling of such synthetic modules allows to describe and evaluate the installed networks. Further, predictions can be made which are subsequently experimentally verified. However, the installation of synthetic gene circuits in cells of higher eukaryotes is still not sufficiently solved. In particular, primary cells or even animals are so far not available for systems biology approaches. This is attributed not only to the genetic and phenotypic complexity of these natural systems but also to the lack of reliable methods to implement synthetic cassettes. Here, we show that expression cassettes conferring bimodal expression can be efficiently transduced with the help of lentiviral vectors. The results presented here were obtained with an immortalized cell line (NiH3T3 cells). However, the real strength of lentiviral transduction is that not only cell lines but also primary cells of different origin and from different species can be transduced. Preliminary data demonstrate that transduction of the synthetic gene regulation circuit is feasible and efficient also in primary cells, Importantly, the expression characteristics of the positive feedback loop are comparable to those observed in the NiH3T3 cells (data not shown). Establishing a mathematical model based on stochasticity we further describe and predict expression patterns.</p><p>Stochasticity is frequently encountered in biological systems. In particular, they contribute to transcriptional and translational processes <xref ref-type="bibr" rid="pone.0002372-Blake1">[2]</xref>, <xref ref-type="bibr" rid="pone.0002372-Ozbudak2">[35]</xref>–<xref ref-type="bibr" rid="pone.0002372-Raser1">[38]</xref>. Generally, it would be expected that stochasticity would compromise the precise control of gene expression levels eventually with detrimental consequences. Indeed, this has been described for several situations <xref ref-type="bibr" rid="pone.0002372-Becskei1">[1]</xref>, <xref ref-type="bibr" rid="pone.0002372-Acar1">[16]</xref>, <xref ref-type="bibr" rid="pone.0002372-Fraser1">[39]</xref>. However, stochasticity also provides a mechanism to generate certain diversification of individual cells within a population–and thus is advantageous under certain conditions, e.g. when the cell encounters stress conditions <xref ref-type="bibr" rid="pone.0002372-Blake2">[40]</xref>. Different natural networks in which stochastic expression leads to a phenotypic heterogeneity of isogenic cells have been described <xref ref-type="bibr" rid="pone.0002372-Kaern1">[reviewed in 33]</xref>. This includes the lysogeny decision in phage lambda; here, an initially homogenous population is split due to molecular fluctuation <xref ref-type="bibr" rid="pone.0002372-Arkin1">[41]</xref>. Also, the HIV-1 Tat expression module gives rise to fluctuations that create distinct phenotypes from a single proviral integration <xref ref-type="bibr" rid="pone.0002372-Weinberger1">[42]</xref>.</p><p>In mammalian systems, stochasticity seems to play a particularly important role in development and differentiation processes such as hematopoiesis <xref ref-type="bibr" rid="pone.0002372-Hoang1">[43]</xref>. Such decisions are also encountered in networks governing apoptosis <xref ref-type="bibr" rid="pone.0002372-Legewie1">[44]</xref> and the cell cycle <xref ref-type="bibr" rid="pone.0002372-Xiong1">[34]</xref>. Also, stochasticity can be involved in diseases, in particular in those which arise from transcription factor haploinsufficiency; upon stochastic simulation evidence was given that haploinsufficient systems are more susceptible to disturbances in expression <xref ref-type="bibr" rid="pone.0002372-Hat1">[45]</xref>. Experimental evidence confirmed this result, demonstrating a decreased probability of gene activation in haploinsufficient cells <xref ref-type="bibr" rid="pone.0002372-Magee1">[46]</xref>. In summary, it seems that stochastic fluctuations, whilst problematic in a homeostatic context, may prove to be essential within the context of development and differentiation of cellular populations.</p><p>In previous studies it was demonstrated that a bimodal expression pattern was realized when a repressor and a transactivator compete for the same DNA binding site <xref ref-type="bibr" rid="pone.0002372-Rossi1">[25]</xref> by two interlinked negative feedback loops <xref ref-type="bibr" rid="pone.0002372-Gardner1">[4]</xref>, <xref ref-type="bibr" rid="pone.0002372-Kramer1">[9]</xref> or by positive feedback and mutual expression <xref ref-type="bibr" rid="pone.0002372-Ferrell1">[23]</xref>. Here, we show that a less complex network consisting of a transactivator and its responsive promoter is sufficient for provoking bimodal expression in a mammalian system. In this respect, the described positive feedback modules represent ideal building blocks to artificially provoke all-or-nothing-answers in mammalian cells.</p><p>The mathematical model predicts a hysteretic expression pattern from the synthetic positive feedback loop. This prediction was experimentally verified. Such a strong a strong hysteretic response is expected to occur only from strong positive feedback loops. The hysteretic response of feedback modules can be tuned by the strength of the positive feedback <xref ref-type="bibr" rid="pone.0002372-Xiong1">[34]</xref>. To achieve such a response with the synthetic modules presented here, transactivator mutants eliciting activation at different inducer levels <xref ref-type="bibr" rid="pone.0002372-Baron1">[47]</xref> might be of benefit. Furthermore, the hysteretic response can be modulated to different strengths through the connection of a positive feedback loop with a transcriptional repressor <xref ref-type="bibr" rid="pone.0002372-Kramer2">[10]</xref>.</p><p>This latter approach pinpoints that for the establishment of complex multigene networks, several non-interacting switches (transactivators) or relais (promoters) are needed. Progress has been made in the creation of independent gene regulation for mammalian systems <xref ref-type="bibr" rid="pone.0002372-May1">[for review see 48]</xref>. The combination of different systems within single cells or animals would allow to differentially control distinct circuits. As exemplified by this study the lentiviral transduction of a synthetic positive feedback circuit proves useful. Exploiting this tool for transduction of independently regulated synthetic modules will allow to reliably install more complex networks in mammalian cells including primary cells. Thereby, it provides an important step to synthetically modify these systems for systems biology approaches.</p></sec><sec id="s4"><title>Materials and Methods</title><sec id="s4a"><title>Cell culture</title><p>Murine fibroblast NIH3T3 cells (ATCC CRL-1658) and 293T cells (human embryonal kidney cells transduced with the simian virus 40 large T antigen) were maintained in Dulbecco's modified Eagle's medium (DMEM). Media were supplemented with 10% fetal calf serum, 2 mM L-glutamine and penicillin (10 U/ml), and streptomycin (100 µg/mL). Doxycycline (Dox) was dissolved in ethanol, filtered and stored at −20°C. The genetic insert, pJSARGFP (synonym: TREAutoR3) was described previously <xref ref-type="bibr" rid="pone.0002372-Markusic1">[49]</xref>.</p></sec><sec id="s4b"><title>Gene transduction</title><p>Lentiviral production was achieved through calcium phosphate transfection of four plasmids. For this purpose, 2×10<sup>5</sup> (293T) cells were plated per 9 cm<sup>2</sup> one day prior to transfection. The lentiviral expression module was transfected into 2×10<sup>5</sup> 293T cells, along with the helper plasmids pLP1 (gag-pol), pLP2 (Rev) and pLP/VSV-g (VSV g) (all obtained from Invitrogen). For transfection, 1,3 µg of lentiviral expression plasmid with 0,9 µg of pLP1, 0,3 µg of LP2, and 0,5 µg of pLP/VSVg were mixed with 15 µL 2.5 M CaCl<sub>2</sub>, and adjusted with H<sub>2</sub>O to give a final volume of 150 µL. The DNA solution was dropwise added to 150 µL HEBS (280 mM NaCl, 50 mM HEPES, 1.5 mM Na<sub>2</sub>HPO<sub>4</sub>, pH 7.1) under vigorous vortexing and incubated at room temperature for 10 min. After 12 h the transfection media was renewed and the virus containing supernatant was collected two and three days after transfection. The supernatant was filtered (pore size 0.45 µm) and frozen at −70°C. For infection 1*10<sup>5</sup> NIH3T3 were plated per 9 cm<sup>2</sup> dish a day prior to infection. At the day of infection the viral supernatant was supplemented with Polybrene (8 µg/mL) and added to the cells for 8 h. The multiplicity of infection was set to 1 to obtain single copy integration of the synthetic gene circuit.</p></sec><sec id="s4c"><title>Reporter Gene Expression</title><p>To determine eGFP expression 1×10<sup>5</sup> cells were cultivated per 6 well plate with the indicated Dox concentrations for the indicated time spans. The cells were washed with PBS, trypsinized and centrifuged at 200 g for five minutes. Cells were stained with Propidium iodide (50 µg/mL) to exclude dead cells from the analysis. The cells were analysed for eGFP expression by flow cytometry using a FACSCalibur (Becton Dickinson).</p></sec><sec id="s4d"><title>Mathematical modelling of the binary expression module</title><p>A mathematical description of the synthetic module was developed. This mathematical simulation was based on the reaction scheme proposed in <xref ref-type="fig" rid="pone-0002372-g002">Fig. 2</xref>. The particular regulatory network was modelled as a discrete, multi-valued, continuous time Markov process which is amenable to standard stochastic techniques. The actual simulations were performed using the Gillespie algorithm as implemented in the SimBiology package, produced by The MathWorks. <xref ref-type="bibr" rid="pone.0002372-McAdams1">[31]</xref>, <xref ref-type="bibr" rid="pone.0002372-Gillespie1">[32]</xref>. In essence, the Gillespie algorithm maintains a count of the number of molecules of each species currently in the cell. These counts are updated at random times based on a set of propensities which are defined and calculated using a given reaction scheme.</p></sec></sec></body><back><ref-list><title>References</title><ref id="pone.0002372-Becskei1"><label>1</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Becskei</surname><given-names>A</given-names></name><name name-style="western"><surname>Serrano</surname><given-names>L</given-names></name></person-group>             <year>2000</year>             <article-title>Engineering stability in gene networks by autoregulation.</article-title>             <source>Nature</source>             <volume>405</volume>             <fpage>590</fpage>             <lpage>593</lpage>          </element-citation></ref><ref id="pone.0002372-Blake1"><label>2</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Blake</surname><given-names>WJ</given-names></name><name name-style="western"><surname>Kaern</surname><given-names>M</given-names></name><name name-style="western"><surname>Cantor</surname><given-names>CR</given-names></name><name name-style="western"><surname>Collins</surname><given-names>JJ</given-names></name></person-group>             <year>2003</year>             <article-title>Noise in eukaryotic gene expression.</article-title>             <source>Nature</source>             <volume>422</volume>             <fpage>633</fpage>             <lpage>637</lpage>          </element-citation></ref><ref id="pone.0002372-Elowitz1"><label>3</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Elowitz</surname><given-names>MB</given-names></name><name name-style="western"><surname>Leibler</surname><given-names>S</given-names></name></person-group>             <year>2000</year>             <article-title>A synthetic oscillatory network of transcriptional regulators.</article-title>             <source>Nature</source>             <volume>403</volume>             <fpage>335</fpage>             <lpage>338</lpage>          </element-citation></ref><ref id="pone.0002372-Gardner1"><label>4</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Gardner</surname><given-names>TS</given-names></name><name name-style="western"><surname>Cantor</surname><given-names>CR</given-names></name><name name-style="western"><surname>Collins</surname><given-names>JJ</given-names></name></person-group>             <year>2000</year>             <article-title>Construction of a genetic toggle switch in Escherichia coli.</article-title>             <source>Nature</source>             <volume>403</volume>             <fpage>339</fpage>             <lpage>342</lpage>          </element-citation></ref><ref id="pone.0002372-Guet1"><label>5</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Guet</surname><given-names>CC</given-names></name><name name-style="western"><surname>Elowitz</surname><given-names>MB</given-names></name><name name-style="western"><surname>Hsing</surname><given-names>W</given-names></name><name name-style="western"><surname>Leibler</surname><given-names>S</given-names></name></person-group>             <year>2002</year>             <article-title>Combinatorial synthesis of genetic networks.</article-title>             <source>Science</source>             <volume>296</volume>             <fpage>1466</fpage>             <lpage>1470</lpage>          </element-citation></ref><ref id="pone.0002372-Hasty1"><label>6</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Hasty</surname><given-names>J</given-names></name><name name-style="western"><surname>McMillen</surname><given-names>D</given-names></name><name name-style="western"><surname>Collins</surname><given-names>JJ</given-names></name></person-group>             <year>2002</year>             <article-title>Engineered gene circuits.</article-title>             <source>Nature</source>             <volume>420</volume>             <fpage>224</fpage>             <lpage>230</lpage>          </element-citation></ref><ref id="pone.0002372-Becskei2"><label>7</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Becskei</surname><given-names>A</given-names></name><name name-style="western"><surname>Boselli</surname><given-names>MG</given-names></name><name name-style="western"><surname>van Oudenaarden</surname><given-names>A</given-names></name></person-group>             <year>2004</year>             <article-title>Amplitude control of cell-cycle waves by nuclear import.</article-title>             <source>Nat Cell Biol</source>             <volume>6</volume>             <fpage>451</fpage>             <lpage>457</lpage>          </element-citation></ref><ref id="pone.0002372-Weber1"><label>8</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Weber</surname><given-names>W</given-names></name><name name-style="western"><surname>Fussenegger</surname><given-names>M</given-names></name></person-group>             <year>2006</year>             <article-title>Pharmacologic transgene control systems for gene therapy.</article-title>             <source>J Gene Med</source>             <volume>8</volume>             <fpage>535</fpage>             <lpage>556</lpage>          </element-citation></ref><ref id="pone.0002372-Kramer1"><label>9</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kramer</surname><given-names>BP</given-names></name><name name-style="western"><surname>Viretta</surname><given-names>AU</given-names></name><name name-style="western"><surname>Daoud-El-Baba</surname><given-names>M</given-names></name><name name-style="western"><surname>Aubel</surname><given-names>D</given-names></name><name name-style="western"><surname>Weber</surname><given-names>W</given-names></name><etal/></person-group>             <year>2004</year>             <article-title>An engineered epigenetic transgene switch in mammalian cells.</article-title>             <source>Nat Biotechnol</source>             <volume>22</volume>             <fpage>867</fpage>             <lpage>870</lpage>          </element-citation></ref><ref id="pone.0002372-Kramer2"><label>10</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kramer</surname><given-names>BP</given-names></name><name name-style="western"><surname>Fussenegger</surname><given-names>M</given-names></name></person-group>             <year>2005</year>             <article-title>Hysteresis in a synthetic mammalian gene network.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>102</volume>             <fpage>9517</fpage>             <lpage>9522</lpage>          </element-citation></ref><ref id="pone.0002372-Weber2"><label>11</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Weber</surname><given-names>W</given-names></name><name name-style="western"><surname>Stelling</surname><given-names>J</given-names></name><name name-style="western"><surname>Rimann</surname><given-names>M</given-names></name><name name-style="western"><surname>Keller</surname><given-names>B</given-names></name><name name-style="western"><surname>Daoud-El Baba</surname><given-names>M</given-names></name><etal/></person-group>             <year>2007</year>             <article-title>A synthetic time-delay circuit in mammalian cells and mice.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>104</volume>             <fpage>2643</fpage>             <lpage>2648</lpage>          </element-citation></ref><ref id="pone.0002372-Kay1"><label>12</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kay</surname><given-names>MA</given-names></name><name name-style="western"><surname>Glorioso</surname><given-names>JC</given-names></name><name name-style="western"><surname>Naldini</surname><given-names>L</given-names></name></person-group>             <year>2001</year>             <article-title>Viral vectors for gene therapy: the art of turning infectious agents into vehicles of therapeutics.</article-title>             <source>Nat Med</source>             <volume>7</volume>             <fpage>33</fpage>             <lpage>40</lpage>          </element-citation></ref><ref id="pone.0002372-Loew1"><label>13</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Loew</surname><given-names>R</given-names></name><name name-style="western"><surname>Selevsek</surname><given-names>N</given-names></name><name name-style="western"><surname>Fehse</surname><given-names>B</given-names></name><name name-style="western"><surname>von Laer</surname><given-names>D</given-names></name><name name-style="western"><surname>Baum</surname><given-names>C</given-names></name><etal/></person-group>             <year>2004</year>             <article-title>Simplified generation of high-titer retrovirus producer cells for clinically relevant retroviral vectors by reversible inclusion of a lox-P-flanked marker gene.</article-title>             <source>Mol Ther</source>             <volume>9</volume>             <fpage>738</fpage>             <lpage>746</lpage>          </element-citation></ref><ref id="pone.0002372-Schucht1"><label>14</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Schucht</surname><given-names>R</given-names></name><name name-style="western"><surname>Coroadinha</surname><given-names>AS</given-names></name><name name-style="western"><surname>Zanta-Boussif</surname><given-names>MA</given-names></name><name name-style="western"><surname>Verhoeyen</surname><given-names>E</given-names></name><etal/></person-group>             <year>2006</year>             <article-title>A new generation of retroviral producer cells: predictable and stable virus production by Flp-mediated site-specific integration of retroviral vectors.</article-title>             <source>Mol Ther</source>             <volume>14</volume>             <fpage>285</fpage>             <lpage>292</lpage>          </element-citation></ref><ref id="pone.0002372-Fiering1"><label>15</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Fiering</surname><given-names>S</given-names></name><name name-style="western"><surname>Whitelaw</surname><given-names>E</given-names></name><name name-style="western"><surname>Martin</surname><given-names>DI</given-names></name></person-group>             <year>2000</year>             <article-title>To be or not to be active: the stochastic nature of enhancer action.</article-title>             <source>Bioessays</source>             <volume>22</volume>             <fpage>381</fpage>             <lpage>387</lpage>          </element-citation></ref><ref id="pone.0002372-Acar1"><label>16</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Acar</surname><given-names>M</given-names></name><name name-style="western"><surname>Becskei</surname><given-names>A</given-names></name><name name-style="western"><surname>van Oudenaarden</surname><given-names>A</given-names></name></person-group>             <year>2005</year>             <article-title>Enhancement of cellular memory by reducing stochastic transitions.</article-title>             <source>Nature</source>             <volume>435</volume>             <fpage>228</fpage>             <lpage>232</lpage>          </element-citation></ref><ref id="pone.0002372-Ozbudak1"><label>17</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Ozbudak</surname><given-names>EM</given-names></name><name name-style="western"><surname>Thattai</surname><given-names>M</given-names></name><name name-style="western"><surname>Lim</surname><given-names>HN</given-names></name><name name-style="western"><surname>Shraiman</surname><given-names>BI</given-names></name><name name-style="western"><surname>van Oudenaarden</surname><given-names>A</given-names></name></person-group>             <year>2004</year>             <article-title>Multistability in the lactose utilization network of Escherichia coli.</article-title>             <source>Nature</source>             <volume>427</volume>             <fpage>737</fpage>             <lpage>740</lpage>          </element-citation></ref><ref id="pone.0002372-Pomerening1"><label>18</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Pomerening</surname><given-names>JR</given-names></name><name name-style="western"><surname>Sontag</surname><given-names>ED</given-names></name><name name-style="western"><surname>Ferrell</surname><given-names>JE</given-names><suffix>Jr</suffix></name></person-group>             <year>2003</year>             <article-title>Building a cell cycle oscillator: hysteresis and bistability in the activation of Cdc2.</article-title>             <source>Nat Cell Biol</source>             <volume>5</volume>             <fpage>346</fpage>             <lpage>351</lpage>          </element-citation></ref><ref id="pone.0002372-Pomerening2"><label>19</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Pomerening</surname><given-names>JR</given-names></name><name name-style="western"><surname>Kim</surname><given-names>SY</given-names></name><name name-style="western"><surname>Ferrell</surname><given-names>JE</given-names><suffix>Jr</suffix></name></person-group>             <year>2005</year>             <article-title>Systems-level dissection of the cell-cycle oscillator: bypassing positive feedback produces damped oscillations.</article-title>             <source>Cell</source>             <volume>122</volume>             <fpage>565</fpage>             <lpage>578</lpage>          </element-citation></ref><ref id="pone.0002372-Sha1"><label>20</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Sha</surname><given-names>W</given-names></name><name name-style="western"><surname>Moore</surname><given-names>J</given-names></name><name name-style="western"><surname>Chen</surname><given-names>K</given-names></name><name name-style="western"><surname>Lassaletta</surname><given-names>AD</given-names></name><name name-style="western"><surname>Yi</surname><given-names>CS</given-names></name><etal/></person-group>             <year>2003</year>             <article-title>Hysteresis drives cell-cycle transitions in Xenopus laevis egg extracts.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>100</volume>             <fpage>975</fpage>             <lpage>980</lpage>          </element-citation></ref><ref id="pone.0002372-Serfling1"><label>21</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Serfling</surname><given-names>E</given-names></name><name name-style="western"><surname>Chuvpilo</surname><given-names>S</given-names></name><name name-style="western"><surname>Liu</surname><given-names>J</given-names></name><name name-style="western"><surname>Hofer</surname><given-names>T</given-names></name><name name-style="western"><surname>Palmetshofer</surname><given-names>A</given-names></name></person-group>             <year>2006</year>             <article-title>NFATc1 autoregulation: a crucial step for cell-fate determination.</article-title>             <source>Trends Immunol</source>             <volume>27</volume>             <fpage>461</fpage>             <lpage>469</lpage>          </element-citation></ref><ref id="pone.0002372-Podtschaske1"><label>22</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Podtschaske</surname><given-names>M</given-names></name><name name-style="western"><surname>Benary</surname><given-names>U</given-names></name><name name-style="western"><surname>Zwinger</surname><given-names>S</given-names></name><name name-style="western"><surname>Hofer</surname><given-names>T</given-names></name><name name-style="western"><surname>Radbruch</surname><given-names>A</given-names></name><etal/></person-group>             <year>2007</year>             <article-title>Digital NFATc2 Activation per Cell Transforms Graded T Cell Receptor Activation into an All-or-None IL-2 Expression.</article-title>             <source>PLoS ONE</source>             <volume>2</volume>             <fpage>e935</fpage>          </element-citation></ref><ref id="pone.0002372-Ferrell1"><label>23</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Ferrell</surname><given-names>JE</given-names><suffix>Jr</suffix></name></person-group>             <year>2002</year>             <article-title>Self-perpetuating states in signal transduction: positive feedback, double-negative feedback and bistability.</article-title>             <source>Curr Opin Cell Biol</source>             <volume>14</volume>             <fpage>140</fpage>             <lpage>148</lpage>          </element-citation></ref><ref id="pone.0002372-Kobayashi1"><label>24</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kobayashi</surname><given-names>H</given-names></name><name name-style="western"><surname>Kaern</surname><given-names>M</given-names></name><name name-style="western"><surname>Araki</surname><given-names>M</given-names></name><name name-style="western"><surname>Chung</surname><given-names>K</given-names></name><name name-style="western"><surname>Gardner</surname><given-names>TS</given-names></name><etal/></person-group>             <year>2004</year>             <article-title>Programmable cells: interfacing natural and engineered gene networks.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>101</volume>             <fpage>8414</fpage>             <lpage>8419</lpage>          </element-citation></ref><ref id="pone.0002372-Rossi1"><label>25</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Rossi</surname><given-names>FM</given-names></name><name name-style="western"><surname>Kringstein</surname><given-names>AM</given-names></name><name name-style="western"><surname>Spicher</surname><given-names>A</given-names></name><name name-style="western"><surname>Guicherit</surname><given-names>OM</given-names></name><name name-style="western"><surname>Blau</surname><given-names>HM</given-names></name></person-group>             <year>2000</year>             <article-title>Transcriptional control: rheostat converted to on/off switch.</article-title>             <source>Mol Cell</source>             <volume>6</volume>             <fpage>723</fpage>             <lpage>728</lpage>          </element-citation></ref><ref id="pone.0002372-Becskei3"><label>26</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Becskei</surname><given-names>A</given-names></name><name name-style="western"><surname>Seraphin</surname><given-names>B</given-names></name><name name-style="western"><surname>Serrano</surname><given-names>L</given-names></name></person-group>             <year>2001</year>             <article-title>Positive feedback in eukaryotic gene networks: cell differentiation by graded to binary response conversion.</article-title>             <source>EMBO J</source>             <volume>20</volume>             <fpage>2528</fpage>             <lpage>2535</lpage>          </element-citation></ref><ref id="pone.0002372-Wiznerowicz1"><label>27</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Wiznerowicz</surname><given-names>M</given-names></name><name name-style="western"><surname>Trono</surname><given-names>D</given-names></name></person-group>             <year>2005</year>             <article-title>Harnessing HIV for therapy, basic research and biotechnology.</article-title>             <source>Trends Biotechnol</source>             <volume>23</volume>             <fpage>42</fpage>             <lpage>47</lpage>          </element-citation></ref><ref id="pone.0002372-Delenda1"><label>28</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Delenda</surname><given-names>C</given-names></name></person-group>             <year>2004</year>             <article-title>Lentiviral vectors: optimization of packaging, transduction and gene expression.</article-title>             <source>J Gene Med</source>             <volume>6</volume>             <supplement>Suppl 1</supplement><fpage>S125</fpage>             <lpage>S138</lpage>          </element-citation></ref><ref id="pone.0002372-Li1"><label>29</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Li</surname><given-names>X</given-names></name><name name-style="western"><surname>Zhao</surname><given-names>X</given-names></name><name name-style="western"><surname>Fang</surname><given-names>Y</given-names></name><name name-style="western"><surname>Jiang</surname><given-names>X</given-names></name><name name-style="western"><surname>Duong</surname><given-names>T</given-names></name><etal/></person-group>             <year>1998</year>             <article-title>Generation of destabilized green fluorescent protein as a transcription reporter.</article-title>             <source>J Biol Chem</source>             <volume>273</volume>             <fpage>34970</fpage>             <lpage>34975</lpage>          </element-citation></ref><ref id="pone.0002372-Kringstein1"><label>30</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kringstein</surname><given-names>AM</given-names></name><name name-style="western"><surname>Rossi</surname><given-names>FM</given-names></name><name name-style="western"><surname>Hofmann</surname><given-names>A</given-names></name><name name-style="western"><surname>Blau</surname><given-names>HM</given-names></name></person-group>             <year>1998</year>             <article-title>Graded transcriptional response to different concentrations of a single transactivator.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>95</volume>             <fpage>13670</fpage>             <lpage>13675</lpage>          </element-citation></ref><ref id="pone.0002372-McAdams1"><label>31</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>McAdams</surname><given-names>HH</given-names></name><name name-style="western"><surname>Arkin</surname><given-names>A</given-names></name></person-group>             <year>1997</year>             <article-title>Stochastic mechanisms in gene expression.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>94</volume>             <fpage>814</fpage>             <lpage>819</lpage>          </element-citation></ref><ref id="pone.0002372-Gillespie1"><label>32</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Gillespie</surname><given-names>DT</given-names></name></person-group>             <year>1976</year>             <article-title>A general method for numerically simulating the stochastic time evolution of coupled chemical reactions.</article-title>             <source>Journal of Computational Physics</source>             <volume>22</volume>             <fpage>403</fpage>             <lpage>434</lpage>          </element-citation></ref><ref id="pone.0002372-Kaern1"><label>33</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Kaern</surname><given-names>M</given-names></name><name name-style="western"><surname>Elston</surname><given-names>TC</given-names></name><name name-style="western"><surname>Blake</surname><given-names>WJ</given-names></name><name name-style="western"><surname>Collins</surname><given-names>JJ</given-names></name></person-group>             <year>2005</year>             <article-title>Stochasticity in gene expression: from theories to phenotypes.</article-title>             <source>Nat Rev Genet</source>             <volume>6</volume>             <fpage>451</fpage>             <lpage>464</lpage>          </element-citation></ref><ref id="pone.0002372-Xiong1"><label>34</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Xiong</surname><given-names>W</given-names></name><name name-style="western"><surname>Ferrell</surname><given-names>JE</given-names><suffix>Jr</suffix></name></person-group>             <year>2003</year>             <article-title>A positive-feedback-based bistable ‘memory module’ that governs a cell fate decision.</article-title>             <source>Nature</source>             <volume>426</volume>             <fpage>460</fpage>             <lpage>465</lpage>          </element-citation></ref><ref id="pone.0002372-Ozbudak2"><label>35</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Ozbudak</surname><given-names>EM</given-names></name><name name-style="western"><surname>Thattai</surname><given-names>M</given-names></name><name name-style="western"><surname>Kurtser</surname><given-names>I</given-names></name><name name-style="western"><surname>Grossman</surname><given-names>AD</given-names></name><name name-style="western"><surname>van Oudenaarden</surname><given-names>A</given-names></name></person-group>             <year>2002</year>             <article-title>Regulation of noise in the expression of a single gene.</article-title>             <source>Nat Genet</source>             <volume>31</volume>             <fpage>69</fpage>             <lpage>73</lpage>          </element-citation></ref><ref id="pone.0002372-Elowitz2"><label>36</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Elowitz</surname><given-names>MB</given-names></name><name name-style="western"><surname>Levine</surname><given-names>AJ</given-names></name><name name-style="western"><surname>Siggia</surname><given-names>ED</given-names></name><name name-style="western"><surname>Swain</surname><given-names>PS</given-names></name></person-group>             <year>2002</year>             <article-title>Stochastic gene expression in a single cell.</article-title>             <source>Science</source>             <volume>297</volume>             <fpage>1183</fpage>             <lpage>1186</lpage>          </element-citation></ref><ref id="pone.0002372-Golding1"><label>37</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Golding</surname><given-names>I</given-names></name><name name-style="western"><surname>Paulsson</surname><given-names>J</given-names></name><name name-style="western"><surname>Zawilski</surname><given-names>SM</given-names></name><name name-style="western"><surname>Cox</surname><given-names>EC</given-names></name></person-group>             <year>2005</year>             <article-title>Real-time kinetics of gene activity in individual bacteria.</article-title>             <source>Cell</source>             <volume>123</volume>             <fpage>1025</fpage>             <lpage>1036</lpage>          </element-citation></ref><ref id="pone.0002372-Raser1"><label>38</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Raser</surname><given-names>JM</given-names></name><name name-style="western"><surname>O'Shea</surname><given-names>EK</given-names></name></person-group>             <year>2005</year>             <article-title>Noise in gene expression: origins, consequences, and control.</article-title>             <source>Science</source>             <volume>309</volume>             <fpage>2010</fpage>             <lpage>2013</lpage>          </element-citation></ref><ref id="pone.0002372-Fraser1"><label>39</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Fraser</surname><given-names>HB</given-names></name><name name-style="western"><surname>Hirsh</surname><given-names>AE</given-names></name><name name-style="western"><surname>Giaever</surname><given-names>G</given-names></name><name name-style="western"><surname>Kumm</surname><given-names>J</given-names></name><name name-style="western"><surname>Eisen</surname><given-names>MB</given-names></name></person-group>             <year>2004</year>             <article-title>Noise minimization in eukaryotic gene expression.</article-title>             <source>PLoS Biol</source>             <volume>2</volume>             <fpage>e137</fpage>          </element-citation></ref><ref id="pone.0002372-Blake2"><label>40</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Blake</surname><given-names>WJ</given-names></name><name name-style="western"><surname>Balazsi</surname><given-names>G</given-names></name><name name-style="western"><surname>Kohanski</surname><given-names>MA</given-names></name><name name-style="western"><surname>Isaacs</surname><given-names>FJ</given-names></name><name name-style="western"><surname>Murphy</surname><given-names>KF</given-names></name><name name-style="western"><surname>Kuang</surname><given-names>Y</given-names></name><name name-style="western"><surname>Cantor</surname><given-names>CR</given-names></name><name name-style="western"><surname>Walt</surname><given-names>DR</given-names></name><name name-style="western"><surname>Collins</surname><given-names>JJ</given-names></name></person-group>             <year>2006</year>             <article-title>Phenotypic consequences of promoter-mediated transcriptional noise.</article-title>             <source>Mol Cell</source>             <volume>24</volume>             <fpage>853</fpage>             <lpage>865</lpage>          </element-citation></ref><ref id="pone.0002372-Arkin1"><label>41</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Arkin</surname><given-names>A</given-names></name><name name-style="western"><surname>Ross</surname><given-names>J</given-names></name><name name-style="western"><surname>McAdams</surname><given-names>HH</given-names></name></person-group>             <year>1998</year>             <article-title>Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells.</article-title>             <source>Genetics</source>             <volume>149</volume>             <fpage>1633</fpage>             <lpage>1648</lpage>          </element-citation></ref><ref id="pone.0002372-Weinberger1"><label>42</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Weinberger</surname><given-names>LS</given-names></name><name name-style="western"><surname>Burnett</surname><given-names>JC</given-names></name><name name-style="western"><surname>Toettcher</surname><given-names>JE</given-names></name><name name-style="western"><surname>Arkin</surname><given-names>AP</given-names></name><name name-style="western"><surname>Schaffer</surname><given-names>DV</given-names></name></person-group>             <year>2005</year>             <article-title>Stochastic gene expression in a lentiviral positive-feedback loop: HIV-1 Tat fluctuations drive phenotypic diversity.</article-title>             <source>Cell</source>             <volume>122</volume>             <fpage>169</fpage>             <lpage>182</lpage>          </element-citation></ref><ref id="pone.0002372-Hoang1"><label>43</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Hoang</surname><given-names>T</given-names></name></person-group>             <year>2004</year>             <article-title>The origin of hematopoietic cell type diversity.</article-title>             <source>Oncogene</source>             <volume>23</volume>             <fpage>7188</fpage>             <lpage>7198</lpage>          </element-citation></ref><ref id="pone.0002372-Legewie1"><label>44</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Legewie</surname><given-names>S</given-names></name><name name-style="western"><surname>Bluthgen</surname><given-names>N</given-names></name><name name-style="western"><surname>Herzel</surname><given-names>H</given-names></name></person-group>             <year>2006</year>             <article-title>Mathematical modeling identifies inhibitors of apoptosis as mediators of positive feedback and bistability.</article-title>             <source>PLoS Comput Biol</source>             <volume>2</volume>             <fpage>e120</fpage>          </element-citation></ref><ref id="pone.0002372-Hat1"><label>45</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Hat</surname><given-names>B</given-names></name><name name-style="western"><surname>Paszek</surname><given-names>P</given-names></name><name name-style="western"><surname>Kimmel</surname><given-names>M</given-names></name><name name-style="western"><surname>Piechor</surname><given-names>K</given-names></name><name name-style="western"><surname>Lipniacki</surname><given-names>T</given-names></name></person-group>             <year>2007</year>             <article-title>How the number of alleles influences gene expression.</article-title>             <source>Journal of Statistical Physics</source>             <volume>128 (1–2)</volume>             <fpage>511</fpage>             <lpage>533</lpage>          </element-citation></ref><ref id="pone.0002372-Magee1"><label>46</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Magee</surname><given-names>JA</given-names></name><name name-style="western"><surname>Abdulkadir</surname><given-names>SA</given-names></name><name name-style="western"><surname>Milbrandt</surname><given-names>J</given-names></name></person-group>             <year>2003</year>             <article-title>Haploinsufficiency at the Nkx3.1 locus. A paradigm for stochastic, dosage-sensitive gene regulation during tumor initiation.</article-title>             <source>Cancer Cell</source>             <volume>3</volume>             <fpage>273</fpage>             <lpage>283</lpage>          </element-citation></ref><ref id="pone.0002372-Baron1"><label>47</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Baron</surname><given-names>U</given-names></name><name name-style="western"><surname>Schnappinger</surname><given-names>D</given-names></name><name name-style="western"><surname>Helbl</surname><given-names>V</given-names></name><name name-style="western"><surname>Gossen</surname><given-names>M</given-names></name><name name-style="western"><surname>Hillen</surname><given-names>W</given-names></name><etal/></person-group>             <year>1999</year>             <article-title>Generation of conditional mutants in higher eukaryotes by switching between the expression of two genes.</article-title>             <source>Proc Natl Acad Sci U S A</source>             <volume>96</volume>             <fpage>1013</fpage>             <lpage>1018</lpage>          </element-citation></ref><ref id="pone.0002372-May1"><label>48</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>May</surname><given-names>T</given-names></name><name name-style="western"><surname>Hauser</surname><given-names>H</given-names></name><name name-style="western"><surname>Wirth</surname><given-names>D</given-names></name></person-group>             <year>2006</year>             <article-title>Current status of transcriptional regulation systems.</article-title>             <source>Cytotechnology</source>             <volume>50 (1–3)</volume>             <fpage>109</fpage>             <lpage>119</lpage>          </element-citation></ref><ref id="pone.0002372-Markusic1"><label>49</label><element-citation publication-type="journal" xlink:type="simple">             <person-group person-group-type="author"><name name-style="western"><surname>Markusic</surname><given-names>D</given-names></name><name name-style="western"><surname>Oude-Elferink</surname><given-names>R</given-names></name><name name-style="western"><surname>Das</surname><given-names>AT</given-names></name><name name-style="western"><surname>Berkhout</surname><given-names>B</given-names></name><name name-style="western"><surname>Seppen</surname><given-names>J</given-names></name></person-group>             <year>2005</year>             <article-title>Comparison of single regulated lentiviral vectors with rtTA expression driven by an autoregulatory loop or a constitutive promoter.</article-title>             <source>Nucleic Acids Res</source>             <volume>33</volume>             <fpage>e63</fpage>          </element-citation></ref></ref-list></back></article>