<?xml version='1.0'?>
<!DOCTYPE art SYSTEM 'http://www.biomedcentral.com/xml/article.dtd'>
<art>
   <ui>1754-1611-3-15</ui>
   <ji>1754-1611</ji>
   <fm>
      <dochead>Research</dochead>
      <bibl>
         <title>
            <p>A switchable light-input, light-output system modelled and constructed in yeast</p>
         </title>
         <aug>
            <au id="A1" ce="yes">
               <snm>Sorokina</snm>
               <fnm>Oxana</fnm>
               <insr iid="I1"/>
               <email>o.sorokina@sms.ed.ac.uk</email>
            </au>
            <au id="A2" ce="yes">
               <snm>Kapus</snm>
               <fnm>Anita</fnm>
               <insr iid="I2"/>
               <email>kapusa@brc.hu</email>
            </au>
            <au id="A3">
               <snm>Terecskei</snm>
               <fnm>Kata</fnm>
               <insr iid="I2"/>
               <email>terecs@brc.hu</email>
            </au>
            <au id="A4">
               <snm>Dixon</snm>
               <mi>E</mi>
               <fnm>Laura</fnm>
               <insr iid="I1"/>
               <insr iid="I3"/>
               <email>Laura.Dixon@ed.ac.uk</email>
            </au>
            <au id="A5">
               <snm>Kozma-Bognar</snm>
               <fnm>Laszlo</fnm>
               <insr iid="I2"/>
               <email>kozmab@brc.hu</email>
            </au>
            <au id="A6">
               <snm>Nagy</snm>
               <fnm>Ferenc</fnm>
               <insr iid="I1"/>
               <insr iid="I2"/>
               <email>nagyf@brc.hu</email>
            </au>
            <au ca="yes" id="A7">
               <snm>Millar</snm>
               <mi>J</mi>
               <fnm>Andrew</fnm>
               <insr iid="I1"/>
               <insr iid="I3"/>
               <email>Andrew.Millar@ed.ac.uk</email>
            </au>
         </aug>
         <insg>
            <ins id="I1">
               <p>Institute of Molecular Plant Sciences, The University of Edinburgh, Kings Buildings, Mayfield Road, Edinburgh EH9 3JH, UK</p>
            </ins>
            <ins id="I2">
               <p>Institute of Plant Biology, Biological Research Center, Temesvari krt. 62, H-6726, Szeged, Hungary</p>
            </ins>
            <ins id="I3">
               <p>Centre for Systems Biology at Edinburgh, C.H. Waddington Building, Kings Buildings, Mayfield Road, Edinburgh EH9 3JD, UK</p>
            </ins>
         </insg>
         <source>Journal of Biological Engineering</source>
         <issn>1754-1611</issn>
         <pubdate>2009</pubdate>
         <volume>3</volume>
         <issue>1</issue>
         <fpage>15</fpage>
         <url>http://www.jbioleng.org/content/3/1/15</url>
         <xrefbib>
            <pubidlist>
               <pubid idtype="pmpid">19761615</pubid>
               <pubid idtype="doi">10.1186/1754-1611-3-15</pubid>
            </pubidlist>
         </xrefbib>
      </bibl>
      <history>
         <rec>
            <date>
               <day>28</day>
               <month>4</month>
               <year>2009</year>
            </date>
         </rec>
         <acc>
            <date>
               <day>17</day>
               <month>9</month>
               <year>2009</year>
            </date>
         </acc>
         <pub>
            <date>
               <day>17</day>
               <month>9</month>
               <year>2009</year>
            </date>
         </pub>
      </history>
      <cpyrt>
         <year>2009</year>
         <collab>Sorokina et al; licensee BioMed Central Ltd.</collab>
         <note>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</note>
      </cpyrt>
      <abs>
         <sec>
            <st>
               <p>Abstract</p>
            </st>
            <sec>
               <st>
                  <p>Background</p>
               </st>
               <p>Advances in synthetic biology will require spatio-temporal regulation of biological processes in heterologous host cells. We develop a light-switchable, two-hybrid interaction in yeast, based upon the Arabidopsis proteins PHYTOCHROME A and FAR-RED ELONGATED HYPOCOTYL 1-LIKE. Light input to this regulatory module allows dynamic control of a light-emitting LUCIFERASE reporter gene, which we detect by real-time imaging of yeast colonies on solid media.</p>
            </sec>
            <sec>
               <st>
                  <p>Results</p>
               </st>
               <p>The reversible activation of the phytochrome by red light, and its inactivation by far-red light, is retained. We use this quantitative readout to construct a mathematical model that matches the system's behaviour and predicts the molecular targets for future manipulation.</p>
            </sec>
            <sec>
               <st>
                  <p>Conclusion</p>
               </st>
               <p>Our model, methods and materials together constitute a novel system for a eukaryotic host with the potential to convert a dynamic pattern of light input into a predictable gene expression response. This system could be applied for the regulation of genetic networks - both known and synthetic.</p>
            </sec>
         </sec>
      </abs>
   </fm>
   <meta>
      <classifications>
         <classification id="endnote" subtype="user_supplied_xml" type="bmc"/>
      </classifications>
   </meta>
   <bdy>
      <sec>
         <st>
            <p>Background</p>
         </st>
         <p>Gene expression systems with both spatial and temporal regulation are key components of engineered and synthetic biological networks. Engineered systems generally use a controlled external stimulus to signal to a specific promoter element, producing a rapid and dose-dependent response <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>. The external stimulus, used at the level of both the whole organism and cell culture, has often been a small, cell permeable molecule, which functions as an activator for the corresponding promoters <abbrgrp><abbr bid="B2">2</abbr><abbr bid="B3">3</abbr><abbr bid="B4">4</abbr></abbrgrp>. Heat shock gene promoter systems can also be utilised for conditional gene expression using heat or irradiation as the stimulus <abbrgrp><abbr bid="B5">5</abbr></abbrgrp>.</p>
         <p>The yeast artificial light switchable promoter system proposed by Shimizu-Sato <it>et al</it>. demonstrates many of the advantages of inducible systems, including low background expression, high inducibility, reversibility and dose-dependence <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. It combines these desirable features with non-toxicity and a lack of pleiotropic and unanticipated effects which are inherent properties of chemically inducible systems. This system is based on the properties of the plant phytochrome B photoreceptor (PhyB), which reversibly changes its conformation in response to red (&#955;<sub>max </sub>= 660 nm) or far-red light (&#955;<sub>max </sub>= 730 nm). The far-red light absorbing conformer (PhyB Pfr) binds to the phytochrome interacting factor 3 (PIF3) protein, whereas interaction between the red light absorbing conformer (PhyB Pr) and PIF3 is much less efficient <abbrgrp><abbr bid="B7">7</abbr></abbrgrp>. In the proposed system, PhyB and PIF3 are expressed as chimeric proteins, fused to the DNA-binding (GBD) or the transcriptional activator (GAD) domain of the GAL4 transcription factor, respectively, giving a typical two - hybrid interaction assay. The <it>cis </it>component of the system is the <it>lacZ </it>reporter gene controlled by a GAL4-responsive artificial promoter. In darkness, PhyB-GBD binds the promoter, but does not induce transcription. Red light illumination converts PhyB into the Pfr form, therefore facilitating PhyB-PIF3 interaction, which recruits PIF3-GAD to the GAL4-dependent promoter resulting in the activation of transcription. Subsequent far-red light illumination coverts PhyB Pfr to Pr and this is followed by the dissociation of the PhyB-GBD - PIF3-GAD complex and abrogation of transcription. The authors demonstrated the dose-dependent response of the system and the dynamics of photoreversible activation of the <it>lacZ </it>reporter gene, derived from quantitative liquid culture assays.</p>
         <p>Recently, another genetically encoded signalling system based on PhyB - PIF3 interaction, with different chimeric proteins, has been successfully used for photoswitching of actin assembly through the Cdc42-WASP-Arp2/3 pathway in <it>E.coli </it><abbrgrp><abbr bid="B8">8</abbr></abbrgrp>.</p>
         <p>All phytochromes (PhyA-E) in the model plant <it>Arabidopsis thaliana </it>are capable of light-dependent conformational changes, but interacting proteins have only been investigated for the two most abundant phytochromes (PhyA and PhyB) <abbrgrp><abbr bid="B7">7</abbr><abbr bid="B9">9</abbr><abbr bid="B10">10</abbr></abbrgrp>. FAR-RED ELONGATED HYPOCOTYL 1 (FHY1) and FHY1 LIKE (FHL) proteins control the nuclear import of PhyA via specific interactions with the Pfr conformer <abbrgrp><abbr bid="B11">11</abbr><abbr bid="B12">12</abbr></abbrgrp>. It follows that, besides the PhyB-PIF3 pair, other phytochrome-interacting protein combinations could be employed as the "light sensing" module of the expression system.</p>
         <p>Functional phytochrome receptors consist of the apoprotein and the covalently linked chromophore called phytochromobilin. Since the chromophore is not synthesised in yeast, an analogous compound, phycocyanobilin (PCB), purified from cyanobacteria, is added to the media. PCB is taken up readily by yeast cells and is autoligated by phytochrome apoproteins resulting in photochemically functional phytochrome photoreceptors <abbrgrp><abbr bid="B13">13</abbr><abbr bid="B14">14</abbr><abbr bid="B15">15</abbr></abbrgrp>. When expressed in yeast with PCB, PhyA behaves like other phytochrome receptors: the Pr &#8596; Pfr conversion is controlled by red and far-red light <abbrgrp><abbr bid="B15">15</abbr><abbr bid="B16">16</abbr><abbr bid="B17">17</abbr></abbrgrp>.</p>
         <p>The light switch described by Shimizu-Sato at al., translates light-dependent protein interactions into transcriptional regulation of a selected gene <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. Beta-galactosidase is the most widely used reporter gene in yeast; however, the protein has a half-life of more than 20 hours in this system, and it can be detected <it>in vitro </it>only <abbrgrp><abbr bid="B18">18</abbr></abbrgrp>. By comparison, the firefly luciferase has a 1.5 hour half-life in yeast, and luciferase activity (luminescence) can be monitored in real-time and <it>in vivo</it>, which makes this reporter a better tool for monitoring dynamic changes in transcription, as has been elegantly demonstrated recently through the monitoring of the cell-cycle and respiratory oscillations monitoring in agitated liquid yeast culture <abbrgrp><abbr bid="B19">19</abbr><abbr bid="B20">20</abbr></abbrgrp>.</p>
         <p>Our aim was to create and mathematically model an inducible gene expression system, based on the principles described above, but containing novel components that provide more stringent regulation and <it>in vivo </it>real-time detection of transcription in yeast colonies on solid media.</p>
      </sec>
      <sec>
         <st>
            <p>Results</p>
         </st>
         <sec>
            <st>
               <p>Selection and testing of components for the light inducible expression system</p>
            </st>
            <p>Detection of promoter induction via beta-galactosidase activity is a well characterised method in <it>S. cerevisiae </it>however it requires time-consuming sampling and <it>in vitro </it>analysis. In order to provide a real-time, <it>in vivo </it>detectable reporter in our system, the GAL4-responsive <it>GAL1 </it>promoter was fused to the firefly luciferase gene (<it>GAL1:LUC</it>) (Fig <figr fid="F1">1</figr>). Figure <figr fid="F1">1B</figr> shows the resulting gene circuit in the community standard Systems Biology Graphic Notation (SBGN) <abbrgrp><abbr bid="B21">21</abbr></abbrgrp>. As a constitutive control, the <it>ADH1:LUC </it>(ALCOCHOL DEHYDROGENASE I) construct was prepared and stably integrated into the genome. Yeast colonies prepared as described in Materials and Methods reached a steady state of luminescence 16-18 hr after luciferin was applied (Additional file <supplr sid="S1">1</supplr>). As expected, <it>ADH1:LUC </it>produced much higher light emission than <it>GAL1:LUC </it>independent of the GBD/GAD fusion proteins expressed (Additional files <supplr sid="S1">1</supplr> and <supplr sid="S2">2</supplr>). Separate set of patches were irradiated with red light (R) or far-red light (FR), or R immediately followed by FR (R/FR), or were kept in darkness. R light induced a rapid increase of luminescence in the case of yeast patches expressing <it>GAL1:LUC</it>, but not in the <it>ADH1:LUC</it>-expressing patches or in patches expressing <it>GAL1:LUC </it>without GBD/GAD fusion proteins (Fig <figr fid="F2">2</figr> and Additional files <supplr sid="S1">1</supplr> and <supplr sid="S2">2</supplr>). Luminescence reached a maximum 14-16 hr after the R light, followed by a slow decrease. In contrast, FR light alone induced very low levels of luciferase activity, which was essentially the same, when R light treatments were followed by FR light immediately (Fig <figr fid="F2">2</figr>). Since the relative (fold) induction was the highest in yeast cells having the <it>GAL1:LUC </it>reporter and expressing PHYA-GBD and FHY1-GAD (Fig <figr fid="F2">2A</figr>), this set of interacting proteins were used in further experiments. These results demonstrate that (i) appropriate LUC markers can be used to report phytochrome photoconversion and light-induced protein-protein interactions in our system; (ii) LUC enzyme activity is unaffected by light in yeast and (iii) yeast patches grown on solid media and treated with luciferin represent stable and reliable experimental material for luminescence imaging.</p>
            <suppl id="S1">
               <title>
                  <p>Additional file 1</p>
               </title>
               <text>
                  <p><b>Light switch and luciferin pre-treatment</b>. The data provided demonstrate the behavior of the light switch system from the time point of the luciferin pretreatment. Yeast cells harboring the <it>GAL1:LUC </it>reporter and expressing PHYA-GBD/GAD-FHY1 (A), PHYA-GBD/GAD-FHL (B), or PHYBNT-GBD/GAD-PIF3 (C) fusion protein-pairs were grown in darkness to form patches (merged colonies) for two days at 30&#176;C, treated with 2.5 mM luciferin and transferred to 22&#176;C for 17.5 h. Separate yeast patches were irradiated at 17.5 h with single red (R), or far-red (FR) light pulses, or with red pulses immediately followed by far-red pulses (R/FR), or were kept in darkness (Dark). Time zero corresponds to the time of luciferin application +1 h. Absolute luminescence levels corrected for the background of the camera are shown.</p>
               </text>
               <file name="1754-1611-3-15-S1.PPT">
                  <p>Click here for file</p>
               </file>
            </suppl>
            <suppl id="S2">
               <title>
                  <p>Additional file 2</p>
               </title>
               <text>
                  <p><b>Basal expression level of the <it>GAL1:LUC </it>reporter</b>. The data provided show the basal expression level of the <it>GAL1:LUC </it>reporter in yeast patches grown with and without added PCB. Yeast cells harboring the <it>GAL1:LUC </it>reporter, but lacking any plasmids expressing PHYA, PHYB, FHY1 or FHL proteins were grown in darkness to form patches for two days at 30&#176;C, treated with 2.5 mM luciferin and transferred to 22&#176;C for 21 h. Separate yeast patches were irradiated at 21 h with single red (R), or far-red (FR) light pulses, or with red pulses immediately followed by far-red pulses (R/FR), or were kept in darkness (Dark). Yeast patches were grown on media with PCB (A) or without PCB (B). Time zero corresponds to the time of luciferin application +1 h. Absolute luminescence levels corrected for the background of the camera are shown.</p>
               </text>
               <file name="1754-1611-3-15-S2.PPT">
                  <p>Click here for file</p>
               </file>
            </suppl>
            <fig id="F1">
               <title>
                  <p>Figure 1</p>
               </title>
               <caption>
                  <p>Light-responsive gene promoter system</p>
               </caption>
               <text>
                  <p><b>Light-responsive gene promoter system</b>. Target gene expression is activated by red light, when photoactivated phytochrome A (Pfr_GBD) interacts with PFL_GAD and recruits it to the target promoter with Gal UAS sites for transcription activation. Expression of the gene could be switched off by far red application, when photoinactivated phytochrome (PR_GBD) dissociates from FHL_GAD. (A) Schematic representation. (B) SBGN representation <abbrgrp><abbr bid="B21">21</abbr></abbrgrp>. Numbers correspond to following reactions: 1 Phytochrome Pr-Pfr photoconversion by red light (R). <b>2 </b>Phytochrome Pfr-Pr photoconversion by far red light (FR). <b>3 </b>Pfr-Pr dark reversion. <b>4 </b>Pfr dissociation from SAP (Sequestrated Areas of Phytochrome). <b>5 </b>Pr association to SAP. <b>6 </b>Phytochrome (Pfr_GBD) association with FHL_GAD. <b>7 </b>Pfr_FHL-Pr_FHL transition by FR <b>8 </b>Intermediate Pr_FHL complex dissociation. <b>9 </b>Transcription of luciferase gene. <b>10 </b>Translation of luciferase mRNA. <b>11 </b>Luciferin- luciferase enzymatic reaction</p>
               </text>
               <graphic file="1754-1611-3-15-1"/>
            </fig>
            <fig id="F2">
               <title>
                  <p>Figure 2</p>
               </title>
               <caption>
                  <p>Light switchable light output</p>
               </caption>
               <text>
                  <p><b>Light switchable light output</b>. Yeast cells harboring the GAL1:LUC reporter and expressing PHYA-GBD/GAD-FHY1 (A), PHYA-GBD/GAD-FHL (B), or PHYBNT-GBD/GAD-PIF3 (C) fusion protein-pairs were grown in darkness to form patches (merged colonies) for two days at 30&#176;C, treated with 2.5 mM luciferin and transferred to 22&#176;C for 17.5 h. Separate yeast patches were irradiated with single red (R), or far-red (FR) light pulses, or with red pulses immediately followed by far-red pulses (R/FR), or were kept in darkness (Dark). Luminescence values normalised to the pre-pulse levels are shown; time 0 h is the start of the light treatment. The luciferin pretreatment is shown in Additional file <supplr sid="S1">1</supplr>. E: Selected luminescent images of yeast patches used to obtain data in panel A. I: red light-induced, NI: non-induced dark control, T0: last images before the light pulse. Consecutive luminescent images taken in every two hours are shown. Pictures are displayed in pseudo-colors: red-white or blue-black colors indicate high or low expression levels, respectively.</p>
               </text>
               <graphic file="1754-1611-3-15-2"/>
            </fig>
         </sec>
         <sec>
            <st>
               <p>Luciferase as a reporter for gene expression in yeast</p>
            </st>
            <p>In order to calculate changes in the rate of transcription from real-time luminescence data, it was necessary to determine the relationship between transcription and enzyme activity.</p>
            <p>The copper inducible <it>CUP1 </it>promoter was fused to luciferase gene, expression was induced and <it>CUP1: LUC </it>RNA and luciferase activity were tested over 7 hours time course. Figure <figr fid="F3">3</figr> shows a 3-4 h delay in the induction of LUC activity relative to <it>LUC </it>mRNA expression. These data contributed to determine the kinetic parameters (mRNA half-life, translation rate) for the Luc reporter model.</p>
            <fig id="F3">
               <title>
                  <p>Figure 3</p>
               </title>
               <caption>
                  <p>Kinetics of induction of luciferase mRNA and luciferase activity from the <it>CUP1:LUC </it>reporter gene</p>
               </caption>
               <text>
                  <p><b>Kinetics of induction of luciferase mRNA and luciferase activity from the <it>CUP1:LUC </it>reporter gene</b>. Yeast cells harboring the <it>CUP1:LUC </it>construct were grown in liquid rich media overnight at 30&#176;C. <it>CUP1:LUC </it>expression was induced by 1.5 mM copper-sulphate (final concentration) at time = 0 and aliquots were harvested hourly. The samples were used to prepare crude protein extracts and to isolate total RNA. Luciferase activity was measured by in vitro assays and luciferase mRNA was determined by qRT-PCR reactions.</p>
               </text>
               <graphic file="1754-1611-3-15-3"/>
            </fig>
         </sec>
         <sec>
            <st>
               <p>An unidentified compound functions as a chromophore Text for this sub-section</p>
            </st>
            <p>Phytochromes are chromoproteins consisting of the apoprotein and a covalently linked, linear tetrapyrrole chromophore, phytochromobilin (P&#934;B) <abbrgrp><abbr bid="B22">22</abbr></abbrgrp>. In the absence of chromophore, phytochromes cannot absorb light, do not show light dependent conformation changes and, therefore, do not function as photoreceptors. Phytochrome apoproteins are synthesised in the Pr form in plants and after autoligation of P&#934;B are capable of light absorption and photoconversion into the Pfr conformer. Cyanobacteria (Synechococcus and Synechocystis sp.) harbour phytochrome-like photoreceptors, which use a chromophore (PCB) with similar structure to that of P&#934;B <abbrgrp><abbr bid="B23">23</abbr></abbrgrp>. Plant phytochromes binding PCB are fully functional photoreceptors and, because of the relative ease of PCB purification it is generally used as an exogenously added chromophore <abbrgrp><abbr bid="B24">24</abbr></abbrgrp>. However, it was unclear whether yeast cultures contained chromophore-like compounds that could serve as the chromophore for plant phytochromes. To test this, yeast cells with the <it>GAL1: LUC </it>reporter and expressing PHYA-GBD and FHY1-GAD were grown on media lacking PCB and the same light treatments were administered as in Fig <figr fid="F2">2A</figr>. To our surprise, significant R induction was detected in the absence of PCB (Fig <figr fid="F4">4</figr>). The fold-induction was reduced to 30% compared to the results with added PCB (Fig <figr fid="F4">4</figr> vs. Fig <figr fid="F2">2A</figr>). Moreover, Figure <figr fid="F4">4</figr> shows that FR light alone, or R followed by FR light also gave qualitatively very similar results to the photoreceptor with PCB. The basal expression level of <it>GAL1:LUC </it>was not affected significantly by the presence or absence of PCB (Additional file <supplr sid="S2">2</supplr>). These results demonstrate that an unidentified compound naturally present in yeast, can serve as a chromophore for phytochromes expressed in this heterologous system.</p>
            <fig id="F4">
               <title>
                  <p>Figure 4</p>
               </title>
               <caption>
                  <p>Switching without PCB</p>
               </caption>
               <text>
                  <p><b>Switching without PCB</b>. Yeast cells harboring the <it>GAL1:LUC </it>reporter and expressing PHYA-GBD and GAD-FHY1 fusion proteins were grown in darkness at 30&#176;C for two days on media without PCB. Luciferin and light treatments and imaging were performed as in Fig. 2A.</p>
               </text>
               <graphic file="1754-1611-3-15-4"/>
            </fig>
         </sec>
         <sec>
            <st>
               <p>Model assumptions and structure</p>
            </st>
            <p>To use our regulatory system for synthetic biology we developed an ordinary differential equation (ODE) model of its function based on kinetic data from the literature and experimentally determined parameter values (Fig <figr fid="F1">1A, B</figr> and Additional file <supplr sid="S3">3</supplr>). The model describes all the known phytochrome properties (e.g. photoconversion, dark reversion, sequestration, etc), using yeast phytochrome data to provide a realistic description of the light-switch function (for the detailed model description and structure, see Methods).</p>
            <p>In summary, the model assumptions are:</p>
            <suppl id="S3">
               <title>
                  <p>Additional file 3</p>
               </title>
               <text>
                  <p><b>Model for PhyA_FHL system in SBML format</b>. Model in SBML</p>
               </text>
               <file name="1754-1611-3-15-S3.XML">
                  <p>Click here for file</p>
               </file>
            </suppl>
            <p>1) Overall concentrations of Phys and PIF3/FHY1/FHL are constant,</p>
            <p>2) Before the light impulse all the phytochromes are in the inactive (Pr) form, and sequestered in a slow acting pool. This might be related to their inclusion in SAPs(Sequestrated Areas of Phytochromes)-like structures similar to those observed by microscopy in cytosol <abbrgrp><abbr bid="B15">15</abbr></abbrgrp> (see Methods for more details).</p>
            <p>3) The dark reversion rate is the same for the free Phy and the Phy-FHY1/FHL complex.</p>
            <p>4) Photoconversion or dark reversion of the Pfr_FHY1/FHL complex creates an intermediate step (Pr_FHY1/FHL) with distinct kinetics, prior to dissociation of the complex.</p>
            <p>5) The luciferase - luciferin subsystem is approximated as a steady state before light treatments.</p>
            <p>6) The initial sharp decrease in luminescence following the application of luciferin is due to the diffusion of luciferin from the point of application on the yeast patch into the agar medium.</p>
            <p>This enables the model to provide a good fit to conceptually similar systems with different interaction partners. For example, the adapted model fits Shimizu-Sato's data with good accuracy using parameter values derived from the literature <abbrgrp><abbr bid="B6">6</abbr><abbr bid="B7">7</abbr><abbr bid="B16">16</abbr><abbr bid="B25">25</abbr></abbrgrp>. This model is simpler, for the main part because the slow LacZ degradation obscures the long-term kinetics (see the model equations in Methods and simulation results in Additional file <supplr sid="S4">4</supplr>).</p>
            <suppl id="S4">
               <title>
                  <p>Additional file 4</p>
               </title>
               <text>
                  <p><b>Simulation of PhyB-PIF3 system</b>. The data represent the results of model fitting to the data from Shimizu-Sato et al. <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. <b>A</b>. Far-red reversal timecourse (long-term): Red light (RL) at 0, Far-Red(FRL) at 0.5 h, 1 h, 2 h. <b>B </b>Far-red reversal timecourse(short-term): Red light at 0, Far-Red at 0.5 h. Data are presented with dots and model simulation with solid lines.</p>
               </text>
               <file name="1754-1611-3-15-S4.PDF">
                  <p>Click here for file</p>
               </file>
            </suppl>
            <p>To account for the initial difference in cell density for each yeast patch that affects luminescence intensity (Fig <figr fid="F5">5</figr>) and for the non-uniformity of the solid media, the initial conditions were set for each patch (each experiment) individually, so that each experimental curve is considered for two regimes, "diffusion" and "phytochrome". The former starts from the application of luciferin with initially decreasing luminescence, which approaches an approximate steady state after 17-18 hours. The latter begins from the light treatment and continues to the end (Fig <figr fid="F6">6A</figr>). The first regime fits separately to the diffusion part of the model and provides the initial luciferin level at the time of light application. Such decomposition of the initial conditions is introduced to describe the temporally changing substrate availability that emerges from the solid culture conditions. Modelling the solid media allows a much wider variety of experimental applications in conditions that are most common in yeast and synthetic biology. The diffusion coefficient in our experiments corresponds to diffusion rate of 3 to 10 mm/h, which is in good agreement with literature data for agar gel <abbrgrp><abbr bid="B26">26</abbr></abbrgrp>.</p>
            <fig id="F5">
               <title>
                  <p>Figure 5</p>
               </title>
               <caption>
                  <p>Effect of cell density on LUC luminescence intensity</p>
               </caption>
               <text>
                  <p><b>Effect of cell density on LUC luminescence intensity</b>. Yeast cells harboring the <it>ADH1:LUC </it>reporter were grown in liquid culture at 30&#176;C for 16 h. Dilution series were prepared from the overnight culture and 20 &#956;l of each dilution were pipetted on solid media and allowed to dry for 15 min at 22&#176;C. The patches were treated with 2.5 mM luciferin and imaged subsequently. Average luminescence values of three replicates for each dilution were plotted against the relative cell density. Graphs with linear (A) or semi-logarithmic (B) scale are shown. Error bars represent standard deviation. The culture with 0.01 relative cell density had OD260 = 0.232 and the corresponding patches contained 8.2&#183;10<sup>4 </sup>cells. Yeast cells were counted by a Burker counting chamber.</p>
               </text>
               <graphic file="1754-1611-3-15-5"/>
            </fig>
            <fig id="F6">
               <title>
                  <p>Figure 6</p>
               </title>
               <caption>
                  <p>Time course of luciferase luminescence intensity in different light conditions</p>
               </caption>
               <text>
                  <p><b>Time course of luciferase luminescence intensity in different light conditions</b>. Experiment results (markers) and simulation (solid lines) are presented at the same scale. (A) Yeast cells containing PhyA_GBD and FHL_GAD were incubated in the dark. Luciferin was added at time zero with application of: no pulse (dark),10 min red pulse (RL), or 10 min red pulse followed immediately by 10 min far red pulse (FRL) at 18 h. (B) On-Off experiment with different time intervals between red and far red pulses. Yeast cells containing both PHYA-GBD andd FHL-GAD were preincubated with PCB in the dark. Luciferin was added at time zero. At time 18 h all the cultures were given the 10 min red pulse (RL) or 10 min red pulse with subsequent 10 min far- red pulse at 0 (RL-FRL), 0.5 (RL-0.5 h-FRL), 1 (RL-1 h-FRL), 3 (RL-3 h-FRL) and 9 h (RL-9 h-FRL) after red pulse. A representative dataset of 5 replicates is shown.</p>
               </text>
               <graphic file="1754-1611-3-15-6"/>
            </fig>
            <p>Parameter values for model equations were obtained from fitting the model to all the timeseries data from luciferase imaging (Fig <figr fid="F6">6A, B</figr>), within the parameres ranges derived from the literature (Table <tblr tid="T1">1</tblr>). The parameter data for luciferase protein degradation rate was supported by additional experiments using cycloheximide treatment of yeast cultures constitutively expressing the luciferase gene (data not shown). The degradation rate constant is estimated to be 0.2-0.8, which corresponds to a 0.8-3 hour half-life. This is similar to the value measured in yeast and mammalian cell cultures and plants <abbrgrp><abbr bid="B27">27</abbr><abbr bid="B28">28</abbr><abbr bid="B29">29</abbr></abbrgrp>.</p>
            <tbl id="T1">
               <title>
                  <p>Table 1</p>
               </title>
               <caption>
                  <p>Final parameter set for PhyA-FHL model</p>
               </caption>
               <tblbdy cols="5">
                  <r>
                     <c ca="left">
                        <p>
                           <b>Parameter description</b>
                        </p>
                     </c>
                     <c ca="center">
                        <p>
                           <b>Parameter name</b>
                        </p>
                     </c>
                     <c ca="center">
                        <p>
                           <b>Parameter value</b>
                        </p>
                     </c>
                     <c ca="left">
                        <p>
                           <b>Dimensions</b>
                        </p>
                     </c>
                     <c ca="left">
                        <p>
                           <b>Literature data</b>
                        </p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Total PhyA concentration</p>
                     </c>
                     <c ca="center">
                        <p>Pool_PhyA</p>
                     </c>
                     <c ca="center">
                        <p>4.93E+01</p>
                     </c>
                     <c ca="left">
                        <p>nM</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Total FHL concentration</p>
                     </c>
                     <c ca="center">
                        <p>Pool_FHL</p>
                     </c>
                     <c ca="center">
                        <p>1.00E+02</p>
                     </c>
                     <c ca="left">
                        <p>nM</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of assocoation for PhyA and FHL</p>
                     </c>
                     <c ca="center">
                        <p>K1</p>
                     </c>
                     <c ca="center">
                        <p>1.50E+00</p>
                     </c>
                     <c ca="left">
                        <p>1/nM*h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of basic (background) transcription for GAL4 promoter</p>
                     </c>
                     <c ca="center">
                        <p>k_base</p>
                     </c>
                     <c ca="center">
                        <p>4.50E-01</p>
                     </c>
                     <c ca="left">
                        <p>nM/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate transition to free pool</p>
                     </c>
                     <c ca="center">
                        <p>K2</p>
                     </c>
                     <c ca="center">
                        <p>1.00E-02</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of sequastration to SAPs</p>
                     </c>
                     <c ca="center">
                        <p>K3</p>
                     </c>
                     <c ca="center">
                        <p>9.90E+00</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of Pr_FHL complex dissociation</p>
                     </c>
                     <c ca="center">
                        <p>K_dis</p>
                     </c>
                     <c ca="center">
                        <p>3.00E-02</p>
                     </c>
                     <c ca="left">
                        <p>1/nM*h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of dark reversion</p>
                     </c>
                     <c ca="center">
                        <p>K_rev</p>
                     </c>
                     <c ca="center">
                        <p>2.00E-01</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c ca="left">
                        <p>0.027-1.2 <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B35">35</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of photoactivation by R light</p>
                     </c>
                     <c ca="center">
                        <p>Ka_R</p>
                     </c>
                     <c ca="center">
                        <p>4.96E-03</p>
                     </c>
                     <c ca="left">
                        <p>m2/um</p>
                     </c>
                     <c ca="left">
                        <p>4.96E-03 <abbrgrp><abbr bid="B30">30</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of photoactivation by FR light</p>
                     </c>
                     <c ca="center">
                        <p>Ka_FR</p>
                     </c>
                     <c ca="center">
                        <p>3.55E-05</p>
                     </c>
                     <c ca="left">
                        <p>m2/um</p>
                     </c>
                     <c ca="left">
                        <p>3.55E-05 <abbrgrp><abbr bid="B30">30</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of photoinactivation by R light</p>
                     </c>
                     <c ca="center">
                        <p>Ki_R</p>
                     </c>
                     <c ca="center">
                        <p>7.44E-04</p>
                     </c>
                     <c ca="left">
                        <p>m2/um</p>
                     </c>
                     <c ca="left">
                        <p>7.44E-04 <abbrgrp><abbr bid="B30">30</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of photoinactivation by FR light</p>
                     </c>
                     <c ca="center">
                        <p>Ki_FR</p>
                     </c>
                     <c ca="center">
                        <p>1.70E-03</p>
                     </c>
                     <c ca="left">
                        <p>m2/um</p>
                     </c>
                     <c ca="left">
                        <p>1.70E-03 <abbrgrp><abbr bid="B30">30</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of luciferase RNA degradatin</p>
                     </c>
                     <c ca="center">
                        <p>m_luc</p>
                     </c>
                     <c ca="center">
                        <p>3 (12 min)</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of luciferase translation</p>
                     </c>
                     <c ca="center">
                        <p>p_luc</p>
                     </c>
                     <c ca="center">
                        <p>1.00E+01</p>
                     </c>
                     <c ca="left">
                        <p>1/h*nM</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of luciferase protein degradation</p>
                     </c>
                     <c ca="center">
                        <p>m_LUC</p>
                     </c>
                     <c ca="center">
                        <p>0.8 (48 min)</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c ca="left">
                        <p>0.23-1.3 <abbrgrp><abbr bid="B28">28</abbr><abbr bid="B29">29</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Rate of the luciferin-dependent luciferase degradation</p>
                     </c>
                     <c ca="center">
                        <p>m_LUC1</p>
                     </c>
                     <c ca="center">
                        <p>0.00E+00</p>
                     </c>
                     <c ca="left">
                        <p>1/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Luciferase activity constant</p>
                     </c>
                     <c ca="center">
                        <p>k_kat</p>
                     </c>
                     <c ca="center">
                        <p>4.99E+02</p>
                     </c>
                     <c ca="left">
                        <p>nM/h</p>
                     </c>
                     <c ca="left">
                        <p>144 <abbrgrp><abbr bid="B28">28</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Michaelis constant for luciferin</p>
                     </c>
                     <c ca="center">
                        <p>Km</p>
                     </c>
                     <c ca="center">
                        <p>1.90E+00</p>
                     </c>
                     <c ca="left">
                        <p>mM</p>
                     </c>
                     <c ca="left">
                        <p>0.1-2 <abbrgrp><abbr bid="B28">28</abbr></abbrgrp></p>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Hill coefficient for transcription</p>
                     </c>
                     <c ca="center">
                        <p>a_luc</p>
                     </c>
                     <c ca="center">
                        <p>2.00E+00</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Maximum transcription rate</p>
                     </c>
                     <c ca="center">
                        <p>n_luc</p>
                     </c>
                     <c ca="center">
                        <p>2.94E+01</p>
                     </c>
                     <c ca="left">
                        <p>nM/h</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Light conversin factor</p>
                     </c>
                     <c ca="center">
                        <p>RLU</p>
                     </c>
                     <c ca="center">
                        <p>1.00E+01</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
                  <r>
                     <c cspan="5">
                        <hr/>
                     </c>
                  </r>
                  <r>
                     <c ca="left">
                        <p>Michaelis constant for transcription</p>
                     </c>
                     <c ca="center">
                        <p>g_luc</p>
                     </c>
                     <c ca="center">
                        <p>1.04E+00</p>
                     </c>
                     <c ca="left">
                        <p>nM</p>
                     </c>
                     <c>
                        <p/>
                     </c>
                  </r>
               </tblbdy>
            </tbl>
         </sec>
         <sec>
            <st>
               <p>Model predictions</p>
            </st>
            <p>The refined model both captures qualitative dynamics and enables a quantitative description of the light switching behaviour. Moreover, it enables us to deduce which parameters would be critical for particular behaviours of the system. By varying these parameters we showed that predicted intermediate state during the photoconversion of Pfr_FHL is crucial to match the slow switching off of the observed LUC expression (Fig <figr fid="F7">7A</figr>). The biochemical nature of this state as well as its experimental measurement is the subject of further experiments. Also, according to the model simulation (Fig <figr fid="F7">7B, C</figr>), shortening of reporter protein half-life does not affect the longevity of reversal of the transcription activation but significantly reduces the intensity of the luminescence.</p>
            <fig id="F7">
               <title>
                  <p>Figure 7</p>
               </title>
               <caption>
                  <p>Investigation of the crucial parameters by simulation of the model</p>
               </caption>
               <text>
                  <p><b>Investigation of the crucial parameters by simulation of the model</b>. <b>A</b>. Simulation of different stability of Pr_FHL intermediate state. Dashed line corresponds to initial value of dissociation constant K_dis, derived from fitting to the experimental curves. Dotted line corresponds to 10 time's faster Pr_FHL complex dissociation. Solid line corresponds to 100-times faster complex dissociation. <b>B</b>. Simulation of different stability of Luc reporter. Dashed line corresponds to initial value of degradation rate constant m_luc, derived from fitting to experimental data. Dotted line corresponds to twice faster degradation rate of luciferase. Solid line corresponds to more than 4-times faster degradation rate. <b>C</b>. Same as B, but in semilog scale. Reporter instability affects the amplitude of the signal but not the overall timing of events</p>
               </text>
               <graphic file="1754-1611-3-15-7"/>
            </fig>
            <fig id="F8">
               <title>
                  <p>Figure 8</p>
               </title>
               <caption>
                  <p>Model simulation and predictions</p>
               </caption>
               <text>
                  <p><b>Model simulation and predictions</b>. (A) Single peaks. Long term On-Off model simulation with far red light application in different time intervals after red pulse. Luciferin was added at time zero. At time 18 h all the cultures were given the 10 min red pulse or the 10 min red pulse with a subsequent 10 min far-red pulse at 0, 0.5, 1, 3 and 9 h after the red pulse. (B) Multiple peaks. Model simulation with different time intervals between red and far red pulses. The green curve corresponds to 10 h between RL pulses and 5 h between RL and FRL pulses, the black curve to 24 h between RL pulses and 5 h between RL and FRL pulses, the blue curve to 10 h between RL pulses and 2 min between RL and FRK pulses, and the red curve to 24 h between RL pulses and 2 min between RL and FRL pulses. (C) Possible profile of transcription activation with RL pulse at time 18 h and subsequent RL-FRL pulses with 2 h between RL and 2 min between RL-FRL until time 133 h-"square" shape. (D) Possible profile of transcription activation with RL and FRL pulses of different frequencies of RL (RL-FRL is 2 min everywhere):5 h-24 h-5 h-"sigma" shape</p>
               </text>
               <graphic file="1754-1611-3-15-8"/>
            </fig>
            <p>The light switch model also gives several predictions about the long-term system behaviour (Fig <figr fid="F8">8A</figr>). In particular, we predicted based on experimental data for 50 h, that in the experimental conditions considered the complete removal of the transcriptional activation effect should take a relatively long time (100 hours), and this has been confirmed by experiments (data not shown). Furthermore, with the given dynamics, we can manipulate subsequent applications of R and FR, to achieve a wide range of desirable profiles of transcription activation (Fig <figr fid="F8">8B, c</figr>, and <figr fid="F8">8D</figr>). Fig <figr fid="F8">8B</figr> illustrates the different types of behaviour of light input, which depend on the interval between R and FR treatments. It is clear from simulations that a small interval (2 min) between R and FR causes the increase of transcription rate and, accordingly, the increase of the luminescence intensity with time. Meanwhile, a longer interval (5 hours) produces the stable base line of input oscillations. On the basis of these simulations one can create a specific protocol of light input, combining the given modes as required, and thus obtaining the "square" shape (Fig <figr fid="F8">8C</figr>) with two modes of light regime, and "sigmoid" shape (Fig <figr fid="F8">8D</figr>) with three modes. Thus, the overall system could be used as a tool for the design of experiments with flexible perturbations of the system, for example, by changing the time intervals between.</p>
         </sec>
      </sec>
      <sec>
         <st>
            <p>Discussion</p>
         </st>
         <p>We developed a photo-regulatory genetic switch for yeast cells that combines several desirable properties. In addition to the widely recognised interacting pair PhyB-PIF3 we have tested other possible protein combinations. We found that the PhyA-FHY1 (and PhyA-FHL) pair provides higher induction level with lower background than that of the PhyB-PIF3 pair in our experimental conditions.</p>
         <p>Previous experiments were carried out using agitated liquid yeast cultures at 30&#176;C <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. We used yeast colonies grown on solid media, for the reason that this setup facilitates light treatments, continuous monitoring of luminescence and potentially allows spatial patterning of light input and biological response. Our experimental system represents a reliable, reproducible and simple set up for investigation of dynamic transcription.</p>
         <p>Phytochrome photoperception in yeast has previously been reported with addition of exogenous chromophore (PCB) <abbrgrp><abbr bid="B6">6</abbr><abbr bid="B11">11</abbr></abbrgrp>. Our system also showed light responsiveness without exogenous chromophore, albeit at a lower level (Fig <figr fid="F4">4</figr>). We propose that phytochromes can employ an unidentified compound from yeast as a chromophore, but the light absorbing efficiency of the constituted receptor is less than that of the holoprotein binding PCB. As a result, R treatment induces a reduced amount of phytochrome Pfr, which results in less efficient induction of transcription. The heterologous chromophore may be specific for some yeast strains, or its weak activating effect could have been difficult to detect using previously-employed reporter genes.</p>
         <p>We developed a mathematical model that describes the system and fits the experimental data with great accuracy. The model incorporated experimental variability arising from the cultures on solid media, <it>via </it>substrate diffusion that corresponds to observations and sets the initial substrate concentration. To our knowledge, this is the first attempt to model the specificity of substrate distribution in solid cell culture. The model results demonstrate good overall accuracy (see Fig <figr fid="F6">6A, B</figr>); nevertheless, the success of the fit varies between experiments. For example, in the On-Off experiments (Fig <figr fid="F6">6B</figr>), model fits better to the longer intervals (at least 1 hour between R and FR treatments), compared with the shorter treatment intervals, when FR is given immediately or 30 min after R. It can be seen from the time series that there is only a small quantitative difference between immediate and 30-min-delayed FR, while after a 1 hour delay the shape of the response resembles that of a single R treatment (Fig <figr fid="F6">6B</figr>), differing only in amplitude. This qualitative shift in the system behaviour requires further analysis, which may shed light on mechanism of transcription activation by R light and inactivation by FR light in the system.</p>
         <p>We found that the kinetics of induction were slower in our conditions compared to previously reported experiments <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. For our tests, yeast patches were grown at 30&#176;C for two days, and then, due to technical issues, the plates were moved to 22&#176;C in the imaging chamber, so the effect of light treatments was investigated at 22&#176;C. We found evidence that the system responds more quickly at 30&#176;C (data not shown), but a complete explanation requires further investigations.</p>
         <p>Our system does not display an instantaneous shutting off of target gene expression. It takes a substantial period of time to completely remove the Luc signal after FR treatment. Modelling suggests that this is not simply due to stability of the Luc reporter (see Fig <figr fid="F7">7B, C</figr>), but rather reflects persistent PhyA activity. It should be noted that far-red exposure does not convert all the active PhyA into Pr form, but by itself produces about 3% of Pfr form <abbrgrp><abbr bid="B30">30</abbr></abbrgrp>. Additionally, we propose a residual physical interaction between the Pr form of PhyA and FHY1/FHL as a possible explanation for the slow kinetics. This was supported with model simulations that correspond to the experimental kinetics. However, the properties of the intermediate state remain to be determined.</p>
      </sec>
      <sec>
         <st>
            <p>Conclusion</p>
         </st>
         <p>The current work initially aimed to create a system to provide well-defined, light-induced perturbations in transcription to a genetic oscillatory circuit, to effect entrainment of the oscillations to a rhythmic light regime. The light switchable system presented here meets the requirements for an entrainment tool; moreover the mathematical model will facilitate the design of any desired entrainment mode. Hence, using the light switch with the corresponding model provides a powerful tool for regular perturbing any gene system of interest with a predictable amplitude and period. Moreover, with spatially-patterned light inputs, such as images, the system would allow spatio-temporal regulation, which could facilitate a greater understanding of biological processes in which inter-cellular communication is involved.</p>
      </sec>
      <sec>
         <st>
            <p>Methods</p>
         </st>
         <sec>
            <st>
               <p>Constructs, yeast strains and growth conditions</p>
            </st>
            <p>Plasmids expressing PHYA-GBD, PHYBNT-GBD, GAD-PIF3, GAD-FHY1 and GAD-FHL fusion proteins have been described <abbrgrp><abbr bid="B6">6</abbr><abbr bid="B11">11</abbr><abbr bid="B12">12</abbr></abbrgrp>. PHYBNT corresponds to an N-terminal fragment of PHYB containing residues 1-621. <it>GAL1</it>, <it>CUP1 </it>and <it>ADH1 </it>promoters containing full 5' un-translated regions and the 3' un-translated region (terminator) of the <it>GAL2 </it>and <it>ADH1 </it>gene were amplified from S. cerevisiae PJ69-4A genomic DNA using the following primers:</p>
            <p>GAL1 Fwd: 5'-AAA<ul>GTCGAC</ul>ATTACCACCATATACATATCC-3'</p>
            <p>GAL1 Rev: 5'-TT<ul>GAATTC</ul>TCTCCTTGACGTTAAAGTAT-3'</p>
            <p>CUP1 Fwd: 5'-CGG<ul>GTCGAC</ul>GCCTTGTTACTAGTTAGAAAAAGACATT-3'</p>
            <p>CUP1 Rev: 5'-CGCC<ul>GAATTC</ul>TTTATGTGATGATTGATTGATTGATT-3'</p>
            <p>ADH1 Fwd: 5'-CTG<ul>GGATCC</ul>GCTGCAGGTCGAGATCCGGGATC-3'</p>
            <p>ADH1 Rev: 5'-CGC<ul>GAATTC</ul>TGGAGTTGATTGTATGCTTGGTATAGCTTGA-3'</p>
            <p>GAL2t Fwd: 5'-ATA<ul>CTGCAG</ul>TGCGTTTGAAGTGAGACGC-3'</p>
            <p>GAL2t Rev: 5'-ATA<ul>CCCGGG</ul>TGGAAGAAAGTCCAGGCAAG-3'</p>
            <p>ADH1t Fwd: 5'-CGC<ul>CTGCAG</ul>AGCTTTGGACTTCTTCGCCAGAGGTTT-3'</p>
            <p>ADH1t Rev: 5'-TAT<ul>CCCGGG</ul>GGCCGGTAGAGGTGTGGTCAATAAGAGC-3'</p>
            <p>The firefly luciferase gene was amplified from plasmid pGL3 (Promega) using the following primers:</p>
            <p>LUC Fwd: 5'-ATGAATTCATGGAAGACGCCAAAAACATA-3'</p>
            <p>LUC Rev: 5'-ATACTGCAGTTACACGGCGATCTTTCC-3'</p>
            <p>The promoter:luciferase-terminator constructs were assembled in pBluescript SK plasmid using the restriction sites designed for the PCR primers (sites are underlined in the sequences of the primers above). All plasmids were transformed into <it>E.coli </it>by the SEM method and cultured under standard conditions <abbrgrp><abbr bid="B31">31</abbr></abbrgrp>. The <it>GAL2 </it>or the <it>ADH1 </it>terminator was used for the <it>GAL1</it>, <it>CUP1:LUC </it>or the <it>ADH1:LUC </it>construct, respectively. The constructs were verified by sequencing and re-cloned in the integrating plasmid p&#948;-UB <abbrgrp><abbr bid="B32">32</abbr></abbrgrp>. The final clones were linearized with XhoI and transformed in yeast strain PJ69-4A by standard LiAC/carrier DNA/PEG protocol. Transformants were plated on synthetic dropout media (SD) without uracil SD(-U). Selected strains carrying the <it>GAL1:LUC </it>construct were co-transformed with plasmids pD153 or pGADT7 (Clontech) expressing GBD- or GAD-fusion proteins, respectively <abbrgrp><abbr bid="B6">6</abbr></abbrgrp>. Transformants were selected and maintained on SD(-LW) plates. Preparation of media and transformation of yeast cells was done according to the Clontech Yeast Protocols Handbook (Clontech).</p>
         </sec>
         <sec>
            <st>
               <p>In vivo luminescence imaging, light treatments</p>
            </st>
            <p>2 ml of selective SD media was inoculated with yeast cells and incubated for 16 hr at 30&#176;C with agitation. 20 &#956;l drops of the cultures were transferred to SD agar plates containing 10 &#956;M PCB, irradiated with far-red light at 70 &#956;molm<sup>-2</sup>s<sup>-1 </sup>fluence rate for 10 min and incubated for 48 hr at 30&#176;C in darkness. All further manipulations were conducted under green safety light. Yeast cells formed merged colonies (or patches) with 5-8 mm diameter. 20 &#956;l of 2.5 mM luciferin solution was pipetted at the center of each patch and the plates were transferred in the imaging chamber at 22&#176;C. Images were taken every 15 minutes using a liquid nitrogen-cooled CCD camera (Visitron Systems GmbH, Munich, Germany). Luminescence was quantified using the Metamorph software (Molecular Devices, Downingtown, PA). Unless stated otherwise, light treatments were administered 17-18 hr after the application of luciferin. The duration of each light treatment was 10 min and fluence rate of light (independent of wavelength) was 70 &#956;molm<sup>-2</sup>s<sup>-1</sup>. Red and far-red light was provided by Snap-Lite LED modules (Quantum Devices, Barneveld, WI).</p>
         </sec>
         <sec>
            <st>
               <p>Induction of CUP1:LUC expression, qRT-PCR and in vitro luciferase assays</p>
            </st>
            <p>Yeast cells carrying the <it>CUP1:LUC </it>construct were inoculated in 10 ml of SD(-U) media and were grown for 16 hr at 30&#176;C with agitation. The starter cultures were diluted to a final volume of 100 ml with fresh SD(-U) media. <it>CUP1:LUC </it>expression was induced by adding CuSO<sub>4 </sub>solution to a final concentration of 1.5 mM. Samples were harvested hourly from induced and non-induced cultures. 2 ml or 100 &#956;l of the cultures were pelleted and frozen for RNA quantification or for luciferase assays, respectively. Total RNA was isolated by using the RNeasy Plant Mini Kit (QIAGEN) according to the manufacturer's instructions. cDNA synthesis and qRT-PCR was performed as described <abbrgrp><abbr bid="B33">33</abbr></abbrgrp>. Primers for qRT-PCR were:</p>
            <p>LUC Fwd: 5'-GGAGCACGGAAAGACGATGACGG-3'</p>
            <p>LUC Rev: 5'-ACAAACACAACTCCTCCGCGCA-3'</p>
            <p>ACT1 Fwd: 5'-CACCAACTGGGACGATATGGA-3'</p>
            <p>ACT1 Rev: 5'-GGCAACTCTCAATTCGTTGTAGAA-3'</p>
            <p>Luciferase-specific signals were normalised to ACTIN 1 (ACT1) levels for each sample. For <it>in vitro </it>luciferase activity measurements, frozen cell pellets were re-suspended in 100 &#956;l of Cell Culture Lysis Buffer (Promega), and vigorously vortexed. After incubation on ice for 5 min, cell debris was pelleted by centrifugation and the supernatant was used as crude protein extract. 20 &#956;l of protein extracts was mixed with 30 &#956;l of the Steady-Glo Luciferase Assay Reagent (Promega) in the wells of a microtiter plate and luminescence was measured in the TopCount NXT luminometer (Perkin-Elmer) for an hour after the addition of the reagent. Counts during monitoring were averaged and normalized to total protein content of the extracts. Protein concentrations were determined by the Bradford assay <abbrgrp><abbr bid="B34">34</abbr></abbrgrp>.</p>
         </sec>
         <sec>
            <st>
               <p>Modelling</p>
            </st>
            <sec>
               <st>
                  <p>1. Model description and structure</p>
               </st>
               <p>Our principal model system (Fig <figr fid="F1">1A, B</figr>) includes two chimeric proteins: phytochrome fused to the GAL4 DNA-binding domain (Phy_GBD), in the active (Pfr) and inactive (Pr) forms, and binding protein PIF3 (or FHL/FHY1) fused to the GAL4 activation domain (FHL_GAD). According to existing experimental data, the recombinant phytochromes are quite stable in yeast. Although the light lability of plant PhyA Pfr is well-described, no detectable difference was observed between the stability of the Pfr and Pr forms of oat PhyA over an 80 hour time period in yeast <abbrgrp><abbr bid="B24">24</abbr></abbrgrp>; moreover, no significant decay in the total PhyA and PhyB amounts over 120 hours was reported <abbrgrp><abbr bid="B35">35</abbr><abbr bid="B15">15</abbr></abbrgrp>. This provides the basis for assuming that our model proteins are present constitutively, so neither production nor degradation occurs in the model system. Two pools, Pool_Phy and Pool_PIF3, fulfil the mass conservation laws for Phy and PIF3.</p>
               <p>In plants, the Pr forms of phytochromes are localized in the cytoplasm in the dark and are translocated to the nucleus in their Pfr form after light absorption <abbrgrp><abbr bid="B36">36</abbr></abbrgrp>. In the yeast system, however, all fusion proteins are constitutively nuclear-localized due to the natural Nuclear Localisation Sequence (NLS) present in the GBD tag or the presence of the SV40 NLS motif fused to the GAD fusion partner. Therefore, in this system the only light-dependent event is the interaction of phytochromes with their corresponding protein partners. Taken together, these details give us reason to locate the interacting proteins and the processes of association and dissociation in the nucleus.</p>
               <p>Not instantaneous kinetics of induction (Fig <figr fid="F2">2</figr>) prompted us to suggest the existence of two phytochrome pools: slow and fast. It has been reported that the sequestration of recombinant PhyA into the cytosolic SAPs (sequestered areas of phytochrome) in yeast has no dependence on light <abbrgrp><abbr bid="B15">15</abbr></abbrgrp>. We, therefore, propose the presence of sequestered and free Phy pools (less and more easy to access, respectively) in the nucleus with a reversible interchange occurring between them. We assume that only the free pool is available for binding to its interaction partner, and, thus, the transition between slow (sequestrated) and fast (free) pool is responsible for the shape of initial light response.</p>
               <p>It is well known that the phytochrome photoconversion cross-section (&#963;) for Pr and Pfr forms depends on the wavelength of light. Red (approximately 660 nm) and far red (approximately 730 nm) light are the most effective for Pr &#8594; Pfr and Pfr &#8594; Pr photoconversions, respectively. Nevertheless, it is evident from the cross-section data that the absorption spectra of the Pr and Pfr forms of Phy significantly overlap <abbrgrp><abbr bid="B30">30</abbr></abbrgrp>. This means that monochromatic light of biologically relevant wavelength (i.e. red) does not convert all the Phy to the Pr or Pfr form, but rather determines the specific distribution ratio of the forms in the total Phy pool. We thus have to account for the activation and inactivation of phytochrome by both red and far-red light, so that:</p>
               <p>
                  <display-formula>
                     <graphic file="1754-1611-3-15-i1.gif"/>
                  </display-formula>
               </p>
               <p>Exact values for Ki and Ka for the different wavelengths were adopted from <abbrgrp><abbr bid="B30">30</abbr></abbrgrp>. In the model Pr &#8596; Pfr transitions are applied for both associated and free form of the phytochromes.</p>
               <p>Dark reversion has been reported for PhyA and PhyB in yeast cultures <abbrgrp><abbr bid="B15">15</abbr><abbr bid="B35">35</abbr><abbr bid="B13">13</abbr></abbrgrp>. According to these data, only a fraction of the total Pfr pool is subject to dark reversion (20-40% of the total amount) with a half-life of 20-40 min. For simplicity in the current model we assume a single Pfr pool that is dark reversible and has a longer half-life than the range suggested by Hennig et al <abbrgrp><abbr bid="B13">13</abbr></abbrgrp>; however, the model is still in good agreement with the overall kinetics described in the literature <abbrgrp><abbr bid="B35">35</abbr><abbr bid="B15">15</abbr></abbrgrp>.</p>
               <p>We assume that dark reversion of the complex Pfr_PIF3 (Pfr_FHY1, Pfr_FHL occurs with the same rate. Therefore, both the photoconversion and dark reversion processes contribute to dissociation of the transcriptional activation complex.</p>
               <p>Finally, for the PhyA_FHY1/FHL complexes, we assume the existence of an additional state, Pr_FHY1/FHL, that has the ability to activate transcription to some extent, as it has been previously demonstrated by <abbrgrp><abbr bid="B11">11</abbr></abbrgrp>. Although the reference above corresponds to PhyA, in our experimental conditions PhyB demonstrated the same kinetics (Fig <figr fid="F2">2C</figr>), so we assume the intermediate state for PhyB-PIf3 as well. According to our hypothesis, this complex is produced as an intermediate product of photoconversion of the Pfr_FHY1/FHL complex after FR exposure. Thus, we propose that Pr proteins that have previously been Pfr can interact with FHY1/FHL and activate transcription.</p>
               <p>Mass Action kinetics were used to describe complex formation and dissociation, translocation, translation, and degradation. Transcription was described with a Hill function and the reporter enzymatic reaction follows Michaelis-Menten kinetics (see Fig <figr fid="F1">1B</figr> for the reactions presented).</p>
               <p>The model equations for the PhyA_FHY1/FHL system are presented below:</p>
               <p>
                  <display-formula id="M1">
                     <graphic file="1754-1611-3-15-i2.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M2">
                     <graphic file="1754-1611-3-15-i3.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M3">
                     <graphic file="1754-1611-3-15-i4.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M4">
                     <graphic file="1754-1611-3-15-i5.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M5">
                     <graphic file="1754-1611-3-15-i6.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M6">
                     <graphic file="1754-1611-3-15-i7.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M7">
                     <graphic file="1754-1611-3-15-i8.gif"/>
                  </display-formula>
               </p>
               <p>The equations (1)-(5) describe changes in concentrations of <it>all </it>the phytochrome components, while (6) and (7) correspond to changes in concentrations of luciferase mRNA and protein, respectively.</p>
               <p>Luminescence level is calculated according to the Michaelis-Menten equation:</p>
               <p>
                  <display-formula id="M8">
                     <graphic file="1754-1611-3-15-i9.gif"/>
                  </display-formula>
               </p>
               <p>Light emission is measured in terms of Relative Light Units (RLU) per second and this corresponds to the rate of the light emission reaction for the colony <abbrgrp><abbr bid="B28">28</abbr></abbrgrp>. The parameter RLU is a conversion factor that translates the number of moles of luciferin reacted into the RLU measurement by the instrument. This also accounts for the discrepancies in colony sizes (Fig <figr fid="F5">5</figr>), growth rate, and instrument characteristics.</p>
            </sec>
         </sec>
         <sec>
            <st>
               <p>'Diffusion' part of the model</p>
            </st>
            <p>Our experimental setup involves the application of a relatively small amount of luciferin substrate (20 &#956;l) to a yeast patch, growing in a 100-mm diameter plate on an agar gel of 5-7 mm thickness. We assumed that the initial decrease in luminescence level just after luciferin application predominantly resulted from the diffusion of substrate through the gel. This was confirmed by an additional experiment (Additional file <supplr sid="S5">5</supplr>). Taking into account that the thickness of the gel is much smaller than the diameter of the plate, we assumed that the diffusion of luciferin could be described with the diffusion equation in polar cylindrical coordinates:</p>
            <suppl id="S5">
               <title>
                  <p>Additional File 5</p>
               </title>
               <text>
                  <p><b>Diffusion rate estimation experiment</b>. The data represent the experiment for the estimation of the diffusion rate. Yeast colonies were grown from 20 &#956;l of OD<sub>600 </sub>0.6 cultures and placed at 1 cm intervals as the pattern indicates in the figure above. Cultures were grown following same conditions as described in Materials and Methods. Luciferin was applied to the central well (orange on figure) and images were taken at 6 min intervals for 26 hours. The rate of diffusion was calculated from the successive time intervals between colonies showing luminescence.</p>
               </text>
               <file name="1754-1611-3-15-S5.PDF">
                  <p>Click here for file</p>
               </file>
            </suppl>
            <p>
               <display-formula>
                  <graphic file="1754-1611-3-15-i10.gif"/>
               </display-formula>
            </p>
            <p>The particular solution of form</p>
            <p>
               <display-formula id="M9">
                  <graphic file="1754-1611-3-15-i11.gif"/>
               </display-formula>
            </p>
            <p>was found to fit experimental data with the best accuracy. Here, S is the cytosolic luciferin concentration, D is the diffusion coefficient, r0 is the effective colony radius, A and B are constants of integration.</p>
            <sec>
               <st>
                  <p>2. PHYA_FHL Model Reactions</p>
               </st>
               <p>
                  <display-formula id="M1a">
                     <graphic file="1754-1611-3-15-i12.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M2a">
                     <graphic file="1754-1611-3-15-i13.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M3a">
                     <graphic file="1754-1611-3-15-i14.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M4a">
                     <graphic file="1754-1611-3-15-i15.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M5a">
                     <graphic file="1754-1611-3-15-i16.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M6a">
                     <graphic file="1754-1611-3-15-i17.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M7a">
                     <graphic file="1754-1611-3-15-i18.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M8a">
                     <graphic file="1754-1611-3-15-i19.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M9a">
                     <graphic file="1754-1611-3-15-i20.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M10">
                     <graphic file="1754-1611-3-15-i21.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M11">
                     <graphic file="1754-1611-3-15-i22.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M12">
                     <graphic file="1754-1611-3-15-i23.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M13">
                     <graphic file="1754-1611-3-15-i24.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M14">
                     <graphic file="1754-1611-3-15-i25.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M15">
                     <graphic file="1754-1611-3-15-i26.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M16">
                     <graphic file="1754-1611-3-15-i27.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M17">
                     <graphic file="1754-1611-3-15-i28.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M18">
                     <graphic file="1754-1611-3-15-i29.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M19">
                     <graphic file="1754-1611-3-15-i30.gif"/>
                  </display-formula>
               </p>
            </sec>
         </sec>
         <sec>
            <st>
               <p>Conservation laws</p>
            </st>
            <p>
               <display-formula>
                  <graphic file="1754-1611-3-15-i31.gif"/>
               </display-formula>
            </p>
            <sec>
               <st>
                  <p>3. PHYB-PIF3 Model Reactions (via Shimizu-Sato's system)</p>
               </st>
               <p>Model for the Shimizu-Sato's system has the similar structure but differs in reporter -LacZ. Model lacks the description of the reporter protein kinetics due to the stability of LacZ protein and the overall relative shortness of the timescale investigated in the paper (2 h) (See Additional file <supplr sid="S4">4</supplr> for the data and model simulation).</p>
               <p>
                  <display-formula id="M1b">
                     <graphic file="1754-1611-3-15-i32.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M2b">
                     <graphic file="1754-1611-3-15-i33.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M3b">
                     <graphic file="1754-1611-3-15-i34.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M4b">
                     <graphic file="1754-1611-3-15-i35.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M5b">
                     <graphic file="1754-1611-3-15-i36.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M6b">
                     <graphic file="1754-1611-3-15-i37.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M7b">
                     <graphic file="1754-1611-3-15-i38.gif"/>
                  </display-formula>
               </p>
               <p>
                  <display-formula id="M8b">
                     <graphic file="1754-1611-3-15-i39.gif"/>
                  </display-formula>
               </p>
            </sec>
            <sec>
               <st>
                  <p>4. Estimation of photoconversion rates</p>
               </st>
               <p>For estimation of photoconversion rates we used the data for the photoconversion cross-section of Pr, Pfr andP and Pfr/P ratios at photoequilibrium of type -I phytochrome <abbrgrp><abbr bid="B30">30</abbr></abbrgrp>.</p>
               <p indent="1">&#8226; Ka for R and FR were estimated according to the wavelength in experiment for 660 nm R and 730 nm FR: Ka_R = 4963e-6 (m2/um); Ka_FR = 35.53e-6(m2/um); Ki_R = 743.9e-6(m2/um); Ki_FR = 1701e-6(m2/um)</p>
               <p indent="1">&#8226; Light intensity for the R from experiment conditions is 70 (umol/m2/s) = >70*3600 = 252000(umol/m2/h)</p>
               <p indent="1">&#8226; For FR -is 80 (umol/m2/s) = >80*3600 = 288000 (umol/m2/h)</p>
            </sec>
            <sec>
               <st>
                  <p>5. Fitting to experimental results</p>
               </st>
               <p>The model was developed in SBTOOLBOX2 for MATLAB and fitted with a particle-swarm optimisation algorithm from the SBPD package in SBTOOLBOX2 <abbrgrp><abbr bid="B37">37</abbr><abbr bid="B38">38</abbr></abbrgrp>. Experiments were designed to cover all possible states of the system that have to be addressed in the model. We started with fitting the model to the simple experimental protocol, including dark conditions and red light application with or without the subsequent immediate far-red application (Fig <figr fid="F6">6A</figr>). Dark experiments taken separately provided us with parameter values for the luciferase system (see Table <tblr tid="T1">1</tblr>), namely the degradation and translation rates, that were fixed during the following optimization procedure. Light response parameter values were estimated from R and R-FR experiments. For that the model was simultaneously fitted to five sets of ON-OFF experiments, each containing seven experiments: R, dark and five combination of R followed by FR with intervals 0 h, 0.5 h, 1 h, 3 h and 9 h (Fig <figr fid="F6">6B</figr>). Thus, a total of 35 timeseries (each of 210-360 timepoints) were fitted simultaneously. Fitting results demonstrate a good accuracy (see Fig <figr fid="F6">6A, B</figr>) with the root mean square deviation of 1.9*10<sup>-3</sup>.</p>
               <p>As we aimed to account for increasing variability arising from solid culture conditions, our model parameters comprise the members which appear specific in each experiment. First of all, this relates to parameters corresponding to 'diffusion' section (D, r0, A and B) as they establish the initial conditions by the time of light treatment. Secondly, parameter RLU that accounts for variability in colony size and growth rate also has to be locally estimated for each experiment. Therefore, besides the parameters identical for all experiments shown in Table <tblr tid="T1">1</tblr>, the system contains 5 local parameters that have to be fitted in SBPD individually for each timeseries to obtain the results in the same scale.</p>
            </sec>
            <sec>
               <st>
                  <p>6. Identifiability analysis</p>
               </st>
               <p>The MOTA tool was used to assess the identifiability of the model and reduce the number of free parameters in the fitting procedure <abbrgrp><abbr bid="B39">39</abbr></abbrgrp>. We started from non-identifiable system with 17 parameters and followed the procedure described in <abbrgrp><abbr bid="B39">39</abbr></abbrgrp>. After we fixed in series the functionally related parameters, we ended up with an identifiable 10-parameter system (Table <tblr tid="T1">1</tblr>).</p>
               <p indent="1">&#8226; 500 parameter set were obtained from independent fittings from 500 random initial points</p>
               <p indent="1">&#8226; Parameter sets were analysed to find statistical dependence between parameters</p>
               <p indent="1">&#8226; Seven parameters were fixed to make Pool_Phy, k_base, m_luc1, K_rev, g_Luc, K_dis, k_kat</p>
               <p indent="1">&#8226; The Final model contained 10 free parameters</p>
            </sec>
         </sec>
      </sec>
      <sec>
         <st>
            <p>Competing interests</p>
         </st>
         <p>The authors declare that they have no competing interests.</p>
      </sec>
      <sec>
         <st>
            <p>Authors' contributions</p>
         </st>
         <p>All authors designed experiments; OS, AK, LED and LKB performed experiments and analysed data; KT constructed materials; OS designed and performed modeling, simulation and analysis; OS, LED, LKB, FN and AJM wrote the paper.</p>
      </sec>
   </bdy>
   <bm>
      <ack>
         <sec>
            <st>
               <p>Acknowledgements</p>
            </st>
            <p>We thank Peter H. Quail, Andreas Hiltbrunner and Nancy A. DaSilva for providing the PHYB-N, PIF3 clones, the PHYA, FHY1 and FHL constructs and the p&#948;-UB plasmid, respectively. This work was supported by grants from the EC (EUCLOCK project, LSHG-CT-2006-018741) to F.N. and A.J.M.; from the Hungarian Science Research Fund (OTKA K-73362 to L.K.B. and OTKA NK-60106 to F.N.) and the Hungarian Academy of Sciences (J&#225;nos B&#243;lyai Research Scholarship) to L.K.B.</p>
            <p>The Centre for Systems Biology at Edinburgh is a Centre for Integrative Systems Biology funded by BBSRC and EPSRC award D019621.</p>
            <p>We thank Anatoly Sorokin for the helpful discussions during the model development and analysis.</p>
         </sec>
      </ack>
      <refgrp>
         <bibl id="B1">
            <title>
               <p>Conditional Control of Protein Function</p>
            </title>
            <aug>
               <au>
                  <snm>Banaszynski</snm>
                  <fnm>LA</fnm>
               </au>
               <au>
                  <snm>Wandless</snm>
                  <fnm>TJ</fnm>
               </au>
            </aug>
            <source>Chemistry &amp; Biology</source>
            <pubdate>2006</pubdate>
            <volume>13</volume>
            <fpage>11</fpage>
            <lpage>21</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmpid" link="fulltext">16426967</pubid>
                  <pubid idtype="doi">10.1016/j.chembiol.2005.10.010</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B2">
            <title>
               <p>P{Switch}, a system for spatial and temporal control of gene expression in Drosophila melanogaster</p>
            </title>
            <aug>
               <au>
                  <snm>Roman</snm>
                  <fnm>G</fnm>
               </au>
               <au>
                  <snm>Endo</snm>
                  <fnm>K</fnm>
               </au>
               <au>
                  <snm>Zong</snm>
                  <fnm>L</fnm>
               </au>
               <au>
                  <snm>Davis</snm>
                  <fnm>RL</fnm>
               </au>
            </aug>
            <source>PNAS</source>
            <pubdate>2001</pubdate>
            <volume>98</volume>
            <fpage>12602</fpage>
            <lpage>12607</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">60100</pubid>
                  <pubid idtype="pmpid" link="fulltext">11675496</pubid>
                  <pubid idtype="doi">10.1073/pnas.221303998</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B3">
            <title>
               <p>Tet B or not tet B: Advances in tetracycline-inducible gene expression</p>
            </title>
            <aug>
               <au>
                  <snm>Blau</snm>
                  <fnm>HM</fnm>
               </au>
               <au>
                  <snm>Rossi</snm>
                  <fnm>FMV</fnm>
               </au>
            </aug>
            <source>PNAS</source>
            <pubdate>1999</pubdate>
            <volume>96</volume>
            <fpage>797</fpage>
            <lpage>799</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">33520</pubid>
                  <pubid idtype="pmpid" link="fulltext">9927646</pubid>
                  <pubid idtype="doi">10.1073/pnas.96.3.797</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B4">
            <title>
               <p>Small-molecule regulation of zebrafish gene expression</p>
            </title>
            <aug>
               <au>
                  <snm>Esengil</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Chang</snm>
                  <fnm>V</fnm>
               </au>
               <au>
                  <snm>Mich</snm>
                  <fnm>JK</fnm>
               </au>
               <au>
                  <snm>Chen</snm>
                  <fnm>JK</fnm>
               </au>
            </aug>
            <source>Nature Chemical Biology</source>
            <pubdate>2007</pubdate>
            <volume>3</volume>
            <fpage>154</fpage>
            <lpage>155</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1038/nchembio858</pubid>
                  <pubid idtype="pmpid" link="fulltext">17237798</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B5">
            <title>
               <p>A technique enabling the controlled gene expression in targeted single cells of C. elegans</p>
            </title>
            <aug>
               <au>
                  <snm>Suzuki</snm>
                  <fnm>M</fnm>
               </au>
               <au>
                  <snm>Kamei</snm>
                  <fnm>Y</fnm>
               </au>
               <au>
                  <snm>Yuba</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Takagi</snm>
                  <fnm>S</fnm>
               </au>
            </aug>
            <source>IEEE Xplore</source>
            <pubdate>2006</pubdate>
            <fpage>1</fpage>
            <lpage>5</lpage>
         </bibl>
         <bibl id="B6">
            <title>
               <p>A light-switchable gene promoter system</p>
            </title>
            <aug>
               <au>
                  <snm>Shimizu-Sato</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Huq</snm>
                  <fnm>E</fnm>
               </au>
               <au>
                  <snm>Tepperman</snm>
                  <fnm>JM</fnm>
               </au>
               <au>
                  <snm>Quail</snm>
                  <fnm>PH</fnm>
               </au>
            </aug>
            <source>Nature Biotechnology</source>
            <pubdate>2002</pubdate>
            <fpage>1041</fpage>
            <lpage>1044</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1038/nbt734</pubid>
                  <pubid idtype="pmpid" link="fulltext">12219076</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B7">
            <title>
               <p>Binding of phytochrome B to its nuclear signalling partner PIF3 is reversibly induced by light</p>
            </title>
            <aug>
               <au>
                  <snm>Ni</snm>
                  <fnm>M</fnm>
               </au>
               <au>
                  <snm>Tepperman</snm>
                  <fnm>JM</fnm>
               </au>
               <au>
                  <snm>Quail</snm>
                  <fnm>PH</fnm>
               </au>
            </aug>
            <source>Nature</source>
            <pubdate>1999</pubdate>
            <fpage>781</fpage>
            <lpage>784</lpage>
            <xrefbib>
               <pubid idtype="pmpid" link="fulltext">10466729</pubid>
            </xrefbib>
         </bibl>
         <bibl id="B8">
            <title>
               <p>Genetically encoded photoswitching of actin assembly through the Cdc42-WASP-Arp2/3 complex pathway</p>
            </title>
            <aug>
               <au>
                  <snm>Leung</snm>
                  <fnm>DW</fnm>
               </au>
               <au>
                  <snm>Otomo</snm>
                  <fnm>C</fnm>
               </au>
               <au>
                  <snm>Chory</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>Rosen</snm>
                  <fnm>MK</fnm>
               </au>
            </aug>
            <source>PNAS</source>
            <pubdate>2008</pubdate>
            <volume>105</volume>
            <fpage>12797</fpage>
            <lpage>12802</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">2525560,2525560</pubid>
                  <pubid idtype="pmpid" link="fulltext">18728185</pubid>
                  <pubid idtype="doi">10.1073/pnas.0801232105</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B9">
            <title>
               <p>Phytochromes: photosensory perception and signal transduction</p>
            </title>
            <aug>
               <au>
                  <snm>Quail</snm>
                  <fnm>P</fnm>
               </au>
               <au>
                  <snm>Boylan</snm>
                  <fnm>M</fnm>
               </au>
               <au>
                  <snm>Parks</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>Short</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Xu</snm>
                  <fnm>Y</fnm>
               </au>
               <au>
                  <snm>Wagner</snm>
                  <fnm>D</fnm>
               </au>
            </aug>
            <source>Science</source>
            <pubdate>1995</pubdate>
            <volume>268</volume>
            <fpage>675</fpage>
            <lpage>680</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1126/science.7732376</pubid>
                  <pubid idtype="pmpid" link="fulltext">7732376</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B10">
            <title>
               <p>Functional Characterization of Phytochrome Interacting Factor 3 in Phytochrome-Mediated Light Signal Transduction</p>
            </title>
            <aug>
               <au>
                  <snm>Kim</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>Yi</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Choi</snm>
                  <fnm>G</fnm>
               </au>
               <au>
                  <snm>Shin</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>Song</snm>
                  <fnm>P-S</fnm>
               </au>
               <au>
                  <snm>Choi</snm>
                  <fnm>G</fnm>
               </au>
            </aug>
            <source>The Plant Cell</source>
            <pubdate>2003</pubdate>
            <volume>15</volume>
            <fpage>2399</fpage>
            <lpage>2407</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">197304</pubid>
                  <pubid idtype="pmpid" link="fulltext">14508006</pubid>
                  <pubid idtype="doi">10.1105/tpc.014498</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B11">
            <title>
               <p>Nuclear accumulation of the phytochrome A photoreceptor requires FHY1</p>
            </title>
            <aug>
               <au>
                  <snm>Hiltbrunner</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Viczian</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Bury</snm>
                  <fnm>E</fnm>
               </au>
               <au>
                  <snm>Tscheuschler</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Kircher</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Toth</snm>
                  <fnm>R</fnm>
               </au>
               <au>
                  <snm>Honsberger</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Nagy</snm>
                  <fnm>F</fnm>
               </au>
               <au>
                  <snm>Fankhauser</snm>
                  <fnm>C</fnm>
               </au>
               <au>
                  <snm>Schafer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>Current biology</source>
            <pubdate>2005</pubdate>
            <volume>15</volume>
            <fpage>2125</fpage>
            <lpage>2130</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1016/j.cub.2005.10.042</pubid>
                  <pubid idtype="pmpid" link="fulltext">16332538</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B12">
            <title>
               <p>FHY1 and FHL act together to mediate nuclear accumulation of the phytochrome A photoreceptor</p>
            </title>
            <aug>
               <au>
                  <snm>Hiltbrunner</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Tscheuschler</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Viczian</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Kunkel</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Kircher</snm>
                  <fnm>RS</fnm>
               </au>
               <au>
                  <snm>Schafer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>Plant and cell physiology</source>
            <pubdate>2006</pubdate>
            <volume>47</volume>
            <fpage>1023</fpage>
            <lpage>1034</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1093/pcp/pcj087</pubid>
                  <pubid idtype="pmpid" link="fulltext">16861711</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B13">
            <title>
               <p>Both subunits of the dimeric plant photoreceptor phytochrome require chomophore for Pfr stability</p>
            </title>
            <aug>
               <au>
                  <snm>Hennig</snm>
                  <fnm>L</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>J Biol Chem. </source>
            <pubdate>2001</pubdate>
            <volume>276</volume>
            <issue>11</issue>
            <fpage>7913</fpage>
            <lpage>7918</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1074/jbc.M009793200</pubid>
                  <pubid idtype="pmpid" link="fulltext">11106666</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B14">
            <title>
               <p>Functional analysis of yeast-derived phytochrome A and B phycocyanobilin adducts</p>
            </title>
            <aug>
               <au>
                  <snm>Kunkel</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Neuhaus</snm>
                  <fnm>G</fnm>
               </au>
               <au>
                  <snm>Batschauer</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Chua</snm>
                  <fnm>N-H</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>The Plant Journal</source>
            <pubdate>1996</pubdate>
            <volume>10</volume>
            <issue>4</issue>
            <fpage>625</fpage>
            <lpage>636</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1046/j.1365-313X.1996.10040625.x</pubid>
                  <pubid idtype="pmpid" link="fulltext">8893541</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B15">
            <title>
               <p>In Vivo Characterization of Phytochrome-Phycocyanobilin Adducts in Yeast</p>
            </title>
            <aug>
               <au>
                  <snm>Kunkel</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Speth</snm>
                  <fnm>V</fnm>
               </au>
               <au>
                  <snm>B&#252;che</snm>
                  <fnm>C</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>Journal of Biological Chemistry</source>
            <pubdate>1995</pubdate>
            <volume>270</volume>
            <fpage>20193</fpage>
            <lpage>20200</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1074/jbc.270.34.20193</pubid>
                  <pubid idtype="pmpid" link="fulltext">7650038</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B16">
            <title>
               <p>Recombinant Type A and B Phytochromes from Potato. Transient Absorption Spectroscopy</p>
            </title>
            <aug>
               <au>
                  <snm>Ruddat</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Schmidt</snm>
                  <fnm>P</fnm>
               </au>
               <au>
                  <snm>Gatz</snm>
                  <fnm>C</fnm>
               </au>
               <au>
                  <snm>Braslavsky</snm>
                  <fnm>SE</fnm>
               </au>
               <au>
                  <snm>G&#228;rtner</snm>
                  <fnm>W</fnm>
               </au>
               <au>
                  <snm>Schaffner</snm>
                  <fnm>K</fnm>
               </au>
            </aug>
            <source>Biochemistry</source>
            <pubdate>1996</pubdate>
            <volume>36</volume>
            <fpage>103</fpage>
            <lpage>111</lpage>
            <xrefbib>
               <pubid idtype="doi">10.1021/bi962012w</pubid>
            </xrefbib>
         </bibl>
         <bibl id="B17">
            <title>
               <p>Recombinant Phytochrome A in Yeast Differs by its Spectroscopic and Photochemical Properties from the Major phyA' and is Close to the Minor phyA": Evidence for Posttranslational Modification of the Pigment in Plants</p>
            </title>
            <aug>
               <au>
                  <snm>Sineshchekov</snm>
                  <fnm>V</fnm>
               </au>
               <au>
                  <snm>Hennig</snm>
                  <fnm>L</fnm>
               </au>
               <au>
                  <snm>Lamparter</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Hughes</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>G&#228;rtner</snm>
                  <fnm>W</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>Photochemistry and Photobiology</source>
            <pubdate>2001</pubdate>
            <volume>73</volume>
            <issue>6</issue>
            <fpage>692</fpage>
            <lpage>696</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1562/0031-8655(2001)073&lt;0692:RPAIYD>2.0.CO;2</pubid>
                  <pubid idtype="pmpid" link="fulltext">11421077</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B18">
            <title>
               <p>In vivo half-life of a protein is a function of its amino-terminal residue</p>
            </title>
            <aug>
               <au>
                  <snm>Bachmair</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Finley</snm>
                  <fnm>D</fnm>
               </au>
               <au>
                  <snm>Varshavsky</snm>
                  <fnm>A</fnm>
               </au>
            </aug>
            <source>Science</source>
            <pubdate>1986</pubdate>
            <volume>234</volume>
            <fpage>179</fpage>
            <lpage>186</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1126/science.3018930</pubid>
                  <pubid idtype="pmpid" link="fulltext">3018930</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B19">
            <title>
               <p>Dual Luciferase Assay System for Rapid Assessment of Gene Expression in Saccharomyces cerevisiae</p>
            </title>
            <aug>
               <au>
                  <snm>McNabb</snm>
                  <fnm>DS</fnm>
               </au>
               <au>
                  <snm>Reed</snm>
                  <fnm>R</fnm>
               </au>
               <au>
                  <snm>Marciniak</snm>
                  <fnm>RA</fnm>
               </au>
            </aug>
            <source>Eukaryotic Cell</source>
            <pubdate>2005</pubdate>
            <volume>4</volume>
            <fpage>1539</fpage>
            <lpage>1549</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">1214206</pubid>
                  <pubid idtype="pmpid" link="fulltext">16151247</pubid>
                  <pubid idtype="doi">10.1128/EC.4.9.1539-1549.2005</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B20">
            <title>
               <p>Real-time luminescence monitoring of cell-cycle and respiratory oscillations in yeast</p>
            </title>
            <aug>
               <au>
                  <snm>Robertson</snm>
                  <fnm>JB</fnm>
               </au>
               <au>
                  <snm>Stowers</snm>
                  <fnm>CC</fnm>
               </au>
               <au>
                  <snm>Boczko</snm>
                  <fnm>E</fnm>
               </au>
               <au>
                  <snm>Johnson</snm>
                  <fnm>CH</fnm>
               </au>
            </aug>
            <source>PNAS</source>
            <pubdate>2008</pubdate>
            <volume>105</volume>
            <fpage>17988</fpage>
            <lpage>17993</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">2584751</pubid>
                  <pubid idtype="pmpid" link="fulltext">19004762</pubid>
                  <pubid idtype="doi">10.1073/pnas.0809482105</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B21">
            <title>
               <p>The Systems Biology Graphical Notation</p>
            </title>
            <aug>
               <au>
                  <snm>Nov&#232;re</snm>
                  <fnm>NL</fnm>
               </au>
               <au>
                  <snm>Hucka</snm>
                  <fnm>M</fnm>
               </au>
               <au>
                  <snm>Mi</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Moodie</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Schreiber</snm>
                  <fnm>F</fnm>
               </au>
               <au>
                  <snm>Sorokin</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Demir</snm>
                  <fnm>E</fnm>
               </au>
               <au>
                  <snm>Wegner</snm>
                  <fnm>K</fnm>
               </au>
               <au>
                  <snm>Aladjem</snm>
                  <fnm>MI</fnm>
               </au>
               <au>
                  <snm>Wimalaratne</snm>
                  <fnm>SM</fnm>
               </au>
               <etal/>
            </aug>
            <source>Nature Biotechnology</source>
            <pubdate>2009</pubdate>
            <volume>27</volume>
            <fpage>735</fpage>
            <lpage>741</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1038/nbt.1558</pubid>
                  <pubid idtype="pmpid" link="fulltext">19668183</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B22">
            <title>
               <p>Multiple Heme Oxygenase Family Members Contribute to the Biosynthesis of the Phytochrome Chromophore in Arabidopsis</p>
            </title>
            <aug>
               <au>
                  <snm>Emborg</snm>
                  <fnm>TJ</fnm>
               </au>
               <au>
                  <snm>Walker</snm>
                  <fnm>JM</fnm>
               </au>
               <au>
                  <snm>Noh</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>Vierstra</snm>
                  <fnm>RD</fnm>
               </au>
            </aug>
            <source>Plant Physiology</source>
            <pubdate>2006</pubdate>
            <volume>140</volume>
            <fpage>856</fpage>
            <lpage>868</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">1400562</pubid>
                  <pubid idtype="pmpid" link="fulltext">16428602</pubid>
                  <pubid idtype="doi">10.1104/pp.105.074211</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B23">
            <title>
               <p>The Chromophore Structure of the Cyanobacterial Phytochrome Cph1 As Predicted by Time-Dependent Density Functional Theory</p>
            </title>
            <aug>
               <au>
                  <snm>Matute</snm>
                  <fnm>RA</fnm>
               </au>
               <au>
                  <snm>Contreras</snm>
                  <fnm>R</fnm>
               </au>
            </aug>
            <source>J Phys Chem B. </source>
            <pubdate>2008</pubdate>
            <volume>112</volume>
            <issue>51</issue>
            <fpage>16253</fpage>
            <lpage>16256</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1021/jp807471e</pubid>
                  <pubid idtype="pmpid">19368024</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B24">
            <title>
               <p>Phytochrome assembly in living cells of the yeast Saccharomyces cerevisiae</p>
            </title>
            <aug>
               <au>
                  <snm>Li</snm>
                  <fnm>L</fnm>
               </au>
               <au>
                  <snm>Lagarias</snm>
                  <fnm>JC</fnm>
               </au>
            </aug>
            <source>Proc Natl Acad Sci U S A. </source>
            <pubdate>1994</pubdate>
            <volume>91</volume>
            <issue>26</issue>
            <fpage>12535</fpage>
            <lpage>12539</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">45473</pubid>
                  <pubid idtype="pmpid">7809073</pubid>
                  <pubid idtype="doi">10.1073/pnas.91.26.12535</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B25">
            <title>
               <p>Phytochrome B binds with greater apparent affinity than phytochrome A to the basic helix-loop-helix factor PIF3 in a reaction requiring the PAS domain of PIF3</p>
            </title>
            <aug>
               <au>
                  <snm>Zhu</snm>
                  <fnm>Y</fnm>
               </au>
               <au>
                  <snm>Tepperman</snm>
                  <fnm>JM</fnm>
               </au>
               <au>
                  <snm>Fairchild</snm>
                  <fnm>CD</fnm>
               </au>
               <au>
                  <snm>Quail</snm>
                  <fnm>PH</fnm>
               </au>
            </aug>
            <source>PNAS</source>
            <pubdate>2000</pubdate>
            <volume>97</volume>
            <fpage>13419</fpage>
            <lpage>13424</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">27239</pubid>
                  <pubid idtype="pmpid" link="fulltext">11069292</pubid>
                  <pubid idtype="doi">10.1073/pnas.230433797</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B26">
            <title>
               <p>Principles of assessing bacterial susceptibility to antibiotics using the agar diffusion method</p>
            </title>
            <aug>
               <au>
                  <snm>Bonev</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>Hooper</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>Parisot</snm>
                  <fnm>J</fnm>
               </au>
            </aug>
            <source>Journal of Antimicrobial Chemotherapy</source>
            <pubdate>2008</pubdate>
            <volume>61</volume>
            <issue>6</issue>
            <fpage>1295</fpage>
            <lpage>1301</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1093/jac/dkn090</pubid>
                  <pubid idtype="pmpid" link="fulltext">18339637</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B27">
            <title>
               <p>Modulation of firefly luciferase stability and impact on studies of gene regulation</p>
            </title>
            <aug>
               <au>
                  <snm>Thompson</snm>
                  <fnm>JF</fnm>
               </au>
               <au>
                  <snm>Hayes</snm>
                  <fnm>LS</fnm>
               </au>
               <au>
                  <snm>Lloyd</snm>
                  <fnm>DB</fnm>
               </au>
            </aug>
            <source>Gene</source>
            <pubdate>1991</pubdate>
            <volume>103</volume>
            <fpage>171</fpage>
            <lpage>177</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1016/0378-1119(91)90270-L</pubid>
                  <pubid idtype="pmpid">1889744</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B28">
            <title>
               <p>Kinetic analysis and modeling of firefly luciferase as a quantitative reporter gene in live mammalian cells</p>
            </title>
            <aug>
               <au>
                  <snm>Ignowski</snm>
                  <fnm>JM</fnm>
               </au>
               <au>
                  <snm>Schaffer</snm>
                  <fnm>DV</fnm>
               </au>
            </aug>
            <source>Biotechnol Bioeng. </source>
            <pubdate>2004</pubdate>
            <volume>86</volume>
            <issue>7</issue>
            <fpage>827</fpage>
            <lpage>834</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1002/bit.20059</pubid>
                  <pubid idtype="pmpid" link="fulltext">15162459</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B29">
            <title>
               <p>Reconstruction of transcriptional dynamics from gene reporter data using differential equations</p>
            </title>
            <aug>
               <au>
                  <snm>Finkenst&#228;dt</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>Heron</snm>
                  <fnm>EA</fnm>
               </au>
               <au>
                  <snm>Komorowski</snm>
                  <fnm>M</fnm>
               </au>
               <au>
                  <snm>Edwards</snm>
                  <fnm>K</fnm>
               </au>
               <au>
                  <snm>Tang</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Harper</snm>
                  <fnm>CV</fnm>
               </au>
               <au>
                  <snm>Davis</snm>
                  <fnm>JR</fnm>
               </au>
               <au>
                  <snm>White</snm>
                  <fnm>MR</fnm>
               </au>
               <au>
                  <snm>Millar</snm>
                  <fnm>AJ</fnm>
               </au>
               <au>
                  <snm>Rand</snm>
                  <fnm>DA</fnm>
               </au>
            </aug>
            <source>Bioinformatics. </source>
            <pubdate>2008</pubdate>
            <volume>24</volume>
            <issue>24</issue>
            <fpage>2901</fpage>
            <lpage>2907</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="pmcid">2639297</pubid>
                  <pubid idtype="pmpid" link="fulltext">18974172</pubid>
                  <pubid idtype="doi">10.1093/bioinformatics/btn562</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B30">
            <title>
               <p>Photomorphogenesis in Plants</p>
            </title>
            <aug>
               <au>
                  <snm>Kendrick</snm>
                  <fnm>RE</fnm>
               </au>
               <au>
                  <snm>Cronenberg</snm>
                  <fnm>GHM</fnm>
               </au>
            </aug>
            <publisher>Dordrecht, Kluwer Academic Publishers</publisher>
            <edition>2</edition>
            <pubdate>1994</pubdate>
         </bibl>
         <bibl id="B31">
            <title>
               <p>High efficiency transformation of Escherichia coli with plasmids</p>
            </title>
            <aug>
               <au>
                  <snm>Inoue</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Nojima</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Okayama</snm>
                  <fnm>H</fnm>
               </au>
            </aug>
            <source>Gene</source>
            <pubdate>1990</pubdate>
            <volume>96</volume>
            <fpage>23</fpage>
            <lpage>28</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1016/0378-1119(90)90336-P</pubid>
                  <pubid idtype="pmpid" link="fulltext">2265755</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B32">
            <title>
               <p>Sequential-Integration for the Regulated Insertion of Cloned Genes in Saccharomyces cerevisiae</p>
            </title>
            <aug>
               <au>
                  <snm>Lee</snm>
                  <fnm>FW</fnm>
               </au>
               <au>
                  <snm>Silva</snm>
                  <fnm>NA</fnm>
               </au>
            </aug>
            <source>Biotechnology Progress</source>
            <pubdate>1997</pubdate>
            <volume>13</volume>
            <fpage>368</fpage>
            <lpage>373</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1021/bp970055d</pubid>
                  <pubid idtype="pmpid" link="fulltext">9265775</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B33">
            <title>
               <p>Arabidopsis thaliana Circadian Clock Is Regulated by the Small GTPase LIP1</p>
            </title>
            <aug>
               <au>
                  <snm>Kevei</snm>
                  <fnm>&#201;</fnm>
               </au>
               <au>
                  <snm>Gyula</snm>
                  <fnm>P</fnm>
               </au>
               <au>
                  <snm>Feh&#233;r</snm>
                  <fnm>B</fnm>
               </au>
               <au>
                  <snm>T&#243;th</snm>
                  <fnm>R</fnm>
               </au>
               <au>
                  <snm>Viczi&#225;n</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Kircher</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Rea</snm>
                  <fnm>D</fnm>
               </au>
               <au>
                  <snm>Dorjgotov</snm>
                  <fnm>D</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
               <au>
                  <snm>Millar</snm>
                  <fnm>AJ</fnm>
               </au>
               <etal/>
            </aug>
            <source>Current Biology</source>
            <pubdate>2007</pubdate>
            <volume>17</volume>
            <fpage>1456</fpage>
            <lpage>1464</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1016/j.cub.2007.07.018</pubid>
                  <pubid idtype="pmpid" link="fulltext">17683937</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B34">
            <title>
               <p>A rapid and sensitive method for quantitation of microgram quantities of protein utilizing the principle of protein-dye-binding</p>
            </title>
            <aug>
               <au>
                  <snm>Bradford</snm>
                  <fnm>MM</fnm>
               </au>
            </aug>
            <source>Anal Biochem</source>
            <pubdate>1976</pubdate>
            <volume>72</volume>
            <fpage>248</fpage>
            <lpage>254</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1016/0003-2697(76)90527-3</pubid>
                  <pubid idtype="pmpid" link="fulltext">942051</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B35">
            <title>
               <p>In Vivo Characterization of Chimeric Phytochromes in Yeast</p>
            </title>
            <aug>
               <au>
                  <snm>Eichenberg</snm>
                  <fnm>K</fnm>
               </au>
               <au>
                  <snm>Kunkel</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Kretsch</snm>
                  <fnm>T</fnm>
               </au>
               <au>
                  <snm>Speth</snm>
                  <fnm>V</fnm>
               </au>
               <au>
                  <snm>Sch&#228;fer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>J Biol Chem</source>
            <pubdate>1999</pubdate>
            <volume>274</volume>
            <fpage>354</fpage>
            <lpage>359</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1074/jbc.274.1.354</pubid>
                  <pubid idtype="pmpid" link="fulltext">9867850</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B36">
            <title>
               <p>Intracellular trafficking of photoreceptors during light-induced signal transduction in plants</p>
            </title>
            <aug>
               <au>
                  <snm>Nagy</snm>
                  <fnm>F</fnm>
               </au>
               <au>
                  <snm>Kircher</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Schafer</snm>
                  <fnm>E</fnm>
               </au>
            </aug>
            <source>Journal of Cell Science</source>
            <pubdate>2001</pubdate>
            <volume>114</volume>
            <fpage>475</fpage>
            <lpage>480</lpage>
            <xrefbib>
               <pubid idtype="pmpid" link="fulltext">11171317</pubid>
            </xrefbib>
         </bibl>
         <bibl id="B37">
            <title>
               <p>Systems Biology Toolbox for MATLAB: a computational platform for research in systems biology</p>
            </title>
            <aug>
               <au>
                  <snm>Schmidt</snm>
                  <fnm>H</fnm>
               </au>
               <au>
                  <snm>Jirstrand</snm>
                  <fnm>M</fnm>
               </au>
            </aug>
            <source>Bioinformatics</source>
            <pubdate>2006</pubdate>
            <volume>22</volume>
            <fpage>514</fpage>
            <lpage>515</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1093/bioinformatics/bti799</pubid>
                  <pubid idtype="pmpid" link="fulltext">16317076</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
         <bibl id="B38">
            <title>
               <p>A review of particle swarm optimization. Part II: hybridisation, combinatorial, multicriteria and constrained optimization, and indicative applications</p>
            </title>
            <aug>
               <au>
                  <snm>Banks</snm>
                  <fnm>A</fnm>
               </au>
               <au>
                  <snm>Vincent</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>Anyakoha</snm>
                  <fnm>C</fnm>
               </au>
            </aug>
            <source>Natural Computing: an international journal</source>
            <pubdate>2008</pubdate>
            <volume>7</volume>
            <issue>(1)x</issue>
            <fpage>109</fpage>
            <lpage>124</lpage>
            <xrefbib>
               <pubid idtype="doi">10.1007/s11047-007-9050-z</pubid>
            </xrefbib>
         </bibl>
         <bibl id="B39">
            <title>
               <p>Data-based identifiability analysis of non-linear dynamical models</p>
            </title>
            <aug>
               <au>
                  <snm>Hengl</snm>
                  <fnm>S</fnm>
               </au>
               <au>
                  <snm>Kreutz</snm>
                  <fnm>C</fnm>
               </au>
               <au>
                  <snm>Timmer</snm>
                  <fnm>J</fnm>
               </au>
               <au>
                  <snm>Maiwald</snm>
                  <fnm>T</fnm>
               </au>
            </aug>
            <source>Bioinformatics</source>
            <pubdate>2007</pubdate>
            <volume>23</volume>
            <fpage>2612</fpage>
            <lpage>2618</lpage>
            <xrefbib>
               <pubidlist>
                  <pubid idtype="doi">10.1093/bioinformatics/btm382</pubid>
                  <pubid idtype="pmpid" link="fulltext">17660526</pubid>
               </pubidlist>
            </xrefbib>
         </bibl>
      </refgrp>
   </bm>
</art>
