Skip to main content
  • Methodology article
  • Open access
  • Published:

OptCircuit: An optimization based method for computational design of genetic circuits

Abstract

Background

Recent years has witnessed an increasing number of studies on constructing simple synthetic genetic circuits that exhibit desired properties such as oscillatory behavior, inducer specific activation/repression, etc. It has been widely acknowledged that that task of building circuits to meet multiple inducer-specific requirements is a challenging one. This is because of the incomplete description of component interactions compounded by the fact that the number of ways in which one can chose and interconnect components, increases exponentially with the number of components.

Results

In this paper we introduce OptCircuit, an optimization based framework that automatically identifies the circuit components from a list and connectivity that brings about the desired functionality. Multiple literature sources are used to compile a comprehensive compilation of kinetic descriptions of promoter-protein pairs. The dynamics that govern the interactions between the elements of the genetic circuit are currently modeled using deterministic ordinary differential equations but the framework is general enough to accommodate stochastic simulations. The desired circuit response is abstracted as the maximization/minimization of an appropriately constructed objective function. Computational results for a toggle switch example demonstrate the ability of the framework to generate the complete list of circuit designs of varying complexity that exhibit the desired response. Designs identified for a genetic decoder highlight the ability of OptCircuit to suggest circuit configurations that go beyond the ones compatible with digital logic-based design principles. Finally, the results obtained from the concentration band detector example demonstrate the ability of OptCircuit to design circuits whose responses are contingent on the level of external inducer as well as pinpoint parameters for modification to rectify an existing (non-functional) biological circuit and restore functionality.

Conclusion

Our results demonstrate that OptCircuit framework can serve as a design platform to aid in the construction and finetuning of integrated biological circuits.

Background

Recent years has witnessed an increasing number of studies on constructing simple synthetic genetic circuits that exhibit desired properties such as oscillatory behavior, inducer specific activation/repression, etc. The hope is that these simple circuits are the vanguards of more complex ones with far ranging implications to biotechnology and medicine bringing to fruition the promise of synthetic biology. It has been widely acknowledged that that task of building circuits to meet multiple inducer-specific requirements is a challenging one [1–5]. This is because of the incomplete description of component interactions compounded by the fact that the number of ways in which one can chose and interconnect components, increases exponentially with the number of components. To meet these emerging challenges, in this paper we introduce an optimization based framework that, given an underlying quantitative description, automatically identifies the circuit components from a list and connectivity that brings about the desired functionality.

To date, several small synthetic gene networks that accomplish a specific functionality have been constructed. For example, Gardner and co-workers have developed a genetic toggle switch- a synthetic gene regulatory network that exhibits bistability [6]. Similarly, Elowitz and Leibler have constructed a synthetic circuit termed as repressilator that was designed to produce an oscillatory response [7]. Subsequently, researchers have extended the repressilator circuit design to induce synchronous oscillations [8], design a synthetic gene-metabolic oscillator [9] and many others [10–14]. Several researchers have employed synthetic circuits to investigate the dynamics and inner workings of more complex natural genetic networks. For example, Hooshangi et al. have constructed synthetic transcriptional cascades to investigate the ultrasensitivity and noise propagation in genetic networks [15]. Mangan et al. have investigated the structure and dynamics of the widely occurring feed forward loop motif [16, 17]. Similarly, Becskei and Serrano designed simple gene circuits to examine the effects of autoregulation in gene networks [18].

In addition to uncovering the design principles of natural genetic networks, synthetic genetic networks are now increasingly finding roles in applications ranging from biotechnology, medicine and bio-sensing. For example, Martin et al. have successfully expressed enzymes from plants, yeast and Escherichia coli to produce amorphadine, a precursor to an anti-malarial drug artemisinin [19] and Anderson et al. have engineered the interaction between bacteria and cancer cells to depend on heterologous environmental signals [20]. Similarly, Levskaya et al. have devised a synthetic circuit that switches between different states in response to red light [21]. These ever-expanding applications have spurred the interest for the development of efficient experimental, database and computational techniques to support these efforts [22].

In response to these developments the research community has been rapidly moving towards standardization by creating the Registry of Standard Biological parts [23]. This registry provides a comprehensive compilation of well defined elements of a genetic circuit such as promoters, ribosome binding sites, transcripts, inducer molecules, terminator sites and plasmids among others. The impetus is that these spare parts registries will help usher the development of more rational engineering approaches for designing such circuits. The potential of using modeling and computations to better understand the function of these circuits has already been recognized and mathematical models have been proposed to describe the interactions between genetic elements [24–27].

The recent availability of well-defined spare parts lists and their interactions brings at the forefront the need to develop procedures to design and optimize genetic circuits that exhibit a desired functionality. Previous efforts in this direction include electrical circuit inspired designs proposed by Basu and Weiss [28, 29]. By constructing a library of cellular gates the authors have implemented simple logical functions such as OR, NOT and AND. Similarly, Mason et al. have investigated the behavior of an electronic model of a gene circuit to produce oscillatory behavior [30]. Other efforts include the combinatorial synthesis approach employed by Guet et al. [31]. In this work the authors varied the connectivity of genes and their corresponding promoters thus generating an ensemble of responses from the resulting genetic circuits. This approach, however, becomes intractable for circuits involving a large number of components [2]. Another important consideration associated with the design and fabrication of genetic circuits is the proper matching of kinetic rates of individual elements of the circuits. Several studies have reported that failure to generate the correct response is often due to improper assembly of the basic elements. For example, simulations conducted by Tuttle et al. have confirmed that repressillator circuits constructed by using wild-type promoters do not result in oscillations [32]. Similarly, studies conducted by Hoosangni et al. have revealed that the behavior of a transcriptional cascade depends on the promoter leakiness and expression levels at the previous stage [15]. Several researchers have stressed the need for optimizing the kinetic parameters to ensure functionality and both experimental [33] and computational approaches [2, 34, 35] have been proposed to this end.

To address these questions, in this work we introduce OptCircuit (see Figure 1), an optimization-based framework that (i) automatically identifies the circuit components from a list and connectivity that brings about the desired functionality; (ii) Rectify or redesign an existing (non-functional) biological circuit and restore functionality by modifying an existing component (e.g., through changes in kinetic parameters) and/or identifying additional components to append to the circuit; Multiple literature sources are used to compile a set of kinetic descriptions of promoter-protein, protein-protein and protein-inducer pairs. The dynamics that govern the interactions between the elements of the genetic circuit are currently modeled using deterministic rate equations but the framework is general enough to accommodate stochastic simulations. The desired circuit response is abstracted as the maximization/minimization of an appropriately constructed objective function. Subsequently, an iterative procedure is implemented within our framework to identify an ensemble of circuits that exhibit the desired response. OptCircuit has been applied to a variety of applications ranging from the design of circuits that discriminate between inducer molecules; circuits that detect the combination of inducer molecules (i.e., 2 to 4 genetic decoder) and finally circuits whose responses are dependent on the concentration of the external inducer (concentration band detector).

Figure 1
figure 1

A pictorial illustration of the OptCircuit framework. The three key components of the framework are the basic genetic elements (promoters, transcripts, inducers); the underlying kinetic mechanisms that drive the circuit response and finally the desired behavior of the circuit under construction. These three components are integrated by OptCircuit using an optimization based formulation.

Results

In this section we highlight the capabilities of the OptCircuit framework to design circuits of varying stimulus and complexity. We first examine the design of simple circuit(s) against known architectures that discriminate between inducer molecules. Next, we dialup the complexity of the desired circuit response by seeking circuit configurations that can detect which combination of inducer molecules are present/absent (genetic decoder example). Finally, we test the ability of the framework to identify circuits whose responses are not only dependent on the presence/absence but also on the level of external inducers (concentration band detector example). In addition to identifying optimal configuration of the design variables, in the last example we also explore whether optimizing the kinetic parameters of specific elements within a given circuit can improve its performance.

Inferring circuits with inducer-specific responses

Here we test OptCircuit by generating circuit designs whose responses are contingent on the presence/absence of different inducer molecules and compare the results with known designs [6]. Specifically, in the presence of anhydrotetracyclin (aTc) the desired circuit must express only protein lacI while in response to inducer IPTG the circuit must express only protein tetR.

The desired circuit response is imposed by maximizing the scaled difference between the expression of the desired minus the undesired florescent protein in response to the two different inducers in line with the description provided in the methods section.

M a x i m i z e  Z = ( ( P l a c I a T C − P l a c I I P T G P l a c I a T C ) + ( P t e t R I P T G − P t e t R a T c P t e t R I P T G ) ) / 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyta0KaemyyaeMaemiEaGNaemyAaKMaemyBa0MaemyAaKMaemOEaONaemyzauMaeeiiaaIaeeOwaOLaeyypa0ZaaeWaaeaadaqadaqcfayaamaalaaabaGaemiuaa1aa0baaeaacqWGSbaBcqWGHbqycqWGJbWycqWGjbqsaeaacqWGHbqycqWGubavcqWGdbWqaaGaeyOeI0Iaemiuaa1aa0baaeaacqWGSbaBcqWGHbqycqWGJbWycqWGjbqsaeaacqWGjbqscqWGqbaucqWGubavcqWGhbWraaaabaGaemiuaa1aa0baaeaacqWGSbaBcqWGHbqycqWGJbWycqWGjbqsaeaacqWGHbqycqWGubavcqWGdbWqaaaaaaGccaGLOaGaayzkaaGaey4kaSYaaeWaaKqbagaadaWcaaqaaiabdcfaqnaaDaaabaGaemiDaqNaemyzauMaemiDaqNaemOuaifabaGaemysaKKaemiuaaLaemivaqLaem4raCeaaiabgkHiTiabdcfaqnaaDaaabaGaemiDaqNaemyzauMaemiDaqNaemOuaifabaGaemyyaeMaemivaqLaem4yamgaaaqaaiabdcfaqnaaDaaabaGaemiDaqNaemyzauMaemiDaqNaemOuaifabaGaemysaKKaemiuaaLaemivaqLaem4raCeaaaaaaOGaayjkaiaawMcaaaGaayjkaiaawMcaaiabc+caViabikdaYaaa@8429@
(1.1)

In Eq 1.1, P l a c I a T C , P t e t R a T C MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aa0baaSqaaiabdYgaSjabdggaHjabdogaJjabdMeajbqaaiabdggaHjabdsfaujabdoeadbaakiabcYcaSiabdcfaqnaaDaaaleaacqWG0baDcqWGLbqzcqWG0baDcqWGsbGuaeaacqWGHbqycqWGubavcqWGdbWqaaaaaa@40F9@ represent the levels of transcripts lacI and tetR in presence of inducer aTc and similarly, P l a c I I P T G , P t e t R I P T G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aa0baaSqaaiabdYgaSjabdggaHjabdogaJjabdMeajbqaaiabdMeajjabdcfaqjabdsfaujabdEeahbaakiabcYcaSiabdcfaqnaaDaaaleaacqWG0baDcqWGLbqzcqWG0baDcqWGsbGuaeaacqWGjbqscqWGqbaucqWGubavcqWGhbWraaaaaa@42FB@ represent the levels of lacI and tetR in presence of inducer IPTG respectively.

Using OptCircuit we identify multiple circuits with up to two promoter transcript pairs. The circuit configuration for the best solution is shown in Figure 2A. Interestingly, the configuration is reminiscent of the architecture of the well-studied genetic toggle switch [6]. Briefly, in line with the construction of a genetic toggle switch, in presence of aTC, the activity of protein tetR is suppressed (see Figure 2B). This in turn leads to the expression of protein lacI from Ptet2 promoter (since tetR suppresses expression from Ptet2 promoter) as shown in the Figure 2B. On the other hand, in presence of inducer IPTG, the activity of protein lacI is suppressed which enables expression tetR from Plac1 promoter (see Figure 2C). In terms of computational requirements, a total of 50 iterations required a total of 200 minutes when run on a 16 node LINUX cluster running dual Intel 3.4 GHz Xeon processors. After performing an exhaustive search on circuits having two-promoter transcript pairs, the effect of dialing up the complexity of the designed circuits by allowing for as many as three to four promoter transcript pairs is examined. Our results shown in Figure 3, indicate that in addition to relatively simple circuit designs, akin to known ones, OptCircuit suggests non-intuitive designs with added complexity affording more opportunities for kinetic parameter tuning.

Figure 2
figure 2

Shown in (A) is a relatively simple circuit identified by OptCircuit. Interestingly the circuit is reminiscent of the genetic toggle switch. Similar to the behavior of the toggle switch, in presence of aTc, the activity of tetR protein is suppressed. This enables the expression of lacI as shown in B. Finally, as shown in (C) in presence of IPTG, activity of lacI is suppressed which enables the expression tetR from the P lac promoter. The triangles with open circles at the vertices represent the promoter elements.

Figure 3
figure 3

Alternative circuit configurations proposed by OptCircuit for the first example. OptCircuit is able to identify more complex architectures to realize a particular outcome. The triangles with open circles at the vertices represent the promoter elements.

Design of genetic decoder

In this section, we use OptCircuit to design for more complex responses by constructing a genetic circuit equivalent of a 2–4 bit decoder. A digital decoder is a multiple-input, multiple-output logic circuit that converts coded inputs into coded outputs. Figure 4(a) illustrates the block diagram of a digital decoder and the corresponding truth table is shown in Figure 4(b). In the context of genetic circuits, we seek the design of a circuit architecture that produces four different responses dependent on the presence and/or absence of the sugars glucose (X in Figure 4) and L-arabinose (Y in Figure 4) respectively. Specifically, we would like the circuit to express (i) YFP (F2 in Figure 4) in response to the presence of L-arabinose and absence of glucose, (ii) RFP (F0 in Figure 4) in response to the absence of both glucose and L-arabinose, (iii) BFP (F3 in Figure 4) when both L-arabinose and glucose are present, (iv) GFP (F1 in Figure 4) when L-arabinose is absent but glucose is present. Note that since absence of glucose is known to induce the expression of cAMP, in this work we assume that absence of glucose implies presence of cAMP and equivalently, presence of glucose implies absence of cAMP.

Figure 4
figure 4

A pictorial illustration of the 2–4 bit decoder. In (A), a block diagram for a digital decoder is shown. In (B), the corresponding truth table associated with the decoder shown is illustrated. For the genetic decoder, X and Y represent glucose L-arabinose respectively. Further, F2, F0, F3 and F1 represent YFP, RFP, BFP and GFP responses respectively.

Given N different promoter elements and M transcripts, the total number of design configurations with upto K promoter-transcript pairs is given by (NM)K. This implies that the search space characterizing all circuit configurations is enormous even for relatively modest values for N and M thus preventing its exhaustive navigation. To alleviate this problem, we implemented the OptCircuit framework in a sequential fashion where successive elements are appended to the genetic circuit to meet, one at a time, the four desired responses (see Figure 5). As shown in the Figure, at each step, the objective function values of the ten best circuit architectures are recorded and the circuit producing the best objective value is retained for the next step. The first step shown in Figure 5, involves the expression of YFP under the (-/+) (i.e. absence of glucose and presence of L-arabinose) condition. To this end, we borrowed the circuit configuration from the well studied, feed-forward loop architecture [16, 17]. CRP and AraC are expressed from the constitutive promoters, P c o n s 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aa0baaSqaaiabdogaJjabd+gaVjabd6gaUjabdohaZbqaaiabigdaXaaaaaa@33A5@ and P c o n s 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aa0baaSqaaiabdogaJjabd+gaVjabd6gaUjabdohaZbqaaiabikdaYaaaaaa@33A7@ respectively and YFP is placed under the control of the P BAD promoter. Using this as the seed, the OptCircuit framework is employed to sequentially identify additional components by following the step-wise procedure shown in Figure 6. After the second step (i.e., (-/-) response), our framework identifies the expression of lacI from the P BAD promoter and the expression of RFP from the P lac promoter. In the third step (i.e., (+,+) response) the best objective value was realized following the expression of protein tetR from P BAD and P lac promoters and expression of BFP from the P tet promoter which is repressible by protein tetR (Figure 6 step 3). Finally after the last step ((+,-) response), by allowing for expression of GFP, the additional elements appended to design the decoder include, the expression of protein cI from the P ara and P lac promoters and the expression of proteins tetR and GFP from the P λ promoter.

Figure 5
figure 5

Step wise procedure to design genetic decoder. In the first step, circuit designs to produce YFP in absence of glucose and presence of L-arabinose are identified. Next, in step 2, additional elements are appended to the design from the previous step to produce YFP under (-/+) condition and RFP under (-/-) condition. Similarly at step 3 BFP response in introduced. Finally, additional elements are appended to the best designed identified by OptCircuit at step 3 to accomplish the response shown in step 4.

Figure 6
figure 6

The circuit configurations identified by OptCircuit for the genetic decoder example. Note that the final step (step 4) is shown as a combination of two architectures.

Figure 7(a–d) provides a network illustration describing the workings at each step of the genetic decoder. As shown in the Figure, in the final step, under (-/+) conditions, expression of YFP along with lacI and tetR is induced from P BAD promoter. Further, araC-L-arabinose complex activate expression of cI from P ara promoter which serves to shut OFF expression of GFP. Under -/- conditions, RFP is expression is induced as before; However, observe that P lac promoter induces expression of cI to shut down expression of GFP (see Figure 7d). Under (+/+) condition, BFP expression from P tet promoter is induced as in the previous steps and production of cI from P ara promoter prevents expression of GFP from P λ promoter. Finally, under (+/-) condition, P BAD , P ara and P lac promoters are shut OFF which in turn enables the expression of GFP from P λ promoter.

Figure 7
figure 7

A network illustration describing the workings of the decoder at each step of the stepwise procedure. Activating interactions are represented by pointed arrows and inhibiting interactions are represented by arrows with flat heads. In (a), CRP-cAMP complex together with L-arabinose-araC complex activate expression of YFP from P BAD promoter under (-/+) condition. In (b),(c), (d), the networks at step 2, step 3 and step 4 are illustrated.

The results shown in the previous paragraph also sheds light onto the design principles for the construction of a genetic decoder. Figure 8, illustrates the binary logic schematics for the genetic circuits characterizing each step shown in Figure 6 and Figure 7. Observe that at each step, the OptCircuit framework allows for the addition of components that are activated only when the corresponding inducer conditions are met. In addition, by expressing appropriate repressor molecules, the OptCircuit framework ensures the repression of all the other promoters expressing fluorescent proteins. For example, after Step 1 the AND gate expressing YFP is active only when glucose is absent of L-arabinose is present. Subsequently, a NOT gate logic is introduced after the second step to turn OFF the P lac promoter when YFP is expressed. After the third step, an OR gate with two inputs followed by an NOT gate is introduced. The OR gate combines the indirect repressive effect of tetR that turns OFF the production of BFP if either YFP or RFP are expressed. Finally, after the last step the OR gate with three inputs of tetR and and an OR gate with two inputs of cI followed by a NOT gate are introduced. This ensures that if either one of BFP, RFP or YFP are expressed then GFP is turned OFF and conversely GFP is expressed only if none of the other three reporters are expressed.

Figure 8
figure 8

Binary logic representation of the circuits obtained in Figure 5. The additional elements appended at each step are indicated using black colored logic gates. Observe that at each step, the OptCircuit framework appends an OR gate followed by a NOT gate to generate the required response. The OR gate serves to integrate the indirect effect of the expression under corresponding inducer conditions. A legend describing the representation of the logic gates is presented in supplementary material (see Additional file 1, Figure S2).

The identified circuit design (Figure 6) happens to be consistent with a purely binary logic viewpoint of regulation. This is not the case with all identified designs. For example, one such circuit configuration involves the expression of protein lacI from P ara promoter instead of the expression of protein tetR from the P BAD promoter (see Figure 9) leading to a behavior that is inconsistent with Boolean-only regulation. To illustrate this, consider the truth table of the design shown in Figure 9. Under -/+ condition YFP is expressed and RFP and GFP are shut-off. However, unlike the circuit described in previous paragraph, expression of YFP is not accompanied by expression of tetR and hence the P tet promoter is free to express the fluorescent protein BFP. Nevertheless, OptCircuit identified this circuit configuration as an optimal architecture for a genetic decoder because the employed kinetic description accounts for not only the presence but also the level of each participating molecule needed to activate transcription. Figure 10, provides a comparison of the steady-state levels of proteins tetR and BFP for the circuits described in Figure 6 (step 4) and Figure 9. In circuit (A), the level of tetR is relatively high (~60 nm) which in turn strongly represses the expression from the P tet promoter. This is expected since, in circuit A, expression of YFP is accompanied by expression of tetR. In contrast, in circuit B, even though the level of protein tetR is relatively low (~10 nm), examination of the level of protein BFP suggests that even low levels of protein tetR are able to effectively repress the expression of BFP from the P tet promoter. The low level of tetR is a manifestation of the leaky repression exerted on the P lac promoter by the lacI protein. This observation is further substantiated by expression of protein RFP, albeit at low levels. These results indicate that by taking into account the underlying kinetic description of the interactions, the OptCircuit framework is able to expand upon possible circuit designs by drawing from architectures that may not be valid based on digital logic viewpoint though adequately meet the imposed requirement due to the careful matching of kinetic parameters as often observed in nature.

Figure 9
figure 9

Alternative circuit design predicted by the OptCircuit framework. The corresponding binary logic diagram is also provided for comparison. The colored arrows indicate the expression of corresponding colored fluorescent proteins. The logic gates are represented using established conventions. The gates include, AND, the triangles represent NOT gates or inverters and the crescent shaped gates are the OR gates. A legend describing the representation of the logic gates is presented in supplementary material (see Additional file Figure S2).

Figure 10
figure 10

The steady-state levels of proteins for the circuits shown in Figure 6(step 4) and Figure 9. Circuit (A), represents the circuit shown in final step of Figure 6. (B) represents the circuit in Figure 9.

Design of Concentration Band Detector

With this example, we explore whether OptCircuit can pinpoint design configurations whose responses are dependent not only on the presence/absence of external inducers but also on their concentrations. We use the concentration band detector example [36] to demonstrate the OptCircuit application. Briefly, this circuit expresses high levels of a reporter protein only when the concentration of the external inducer (i.e. L-arabinose) is within a specific range [36] (i.e., neither too high or too low) as shown in Figure 11.

Figure 11
figure 11

The desired response from the concentration band detector example. The solid line represents the desired response of the inducer molecule (GFP). The dotted line represents the desired response of the protein tetR. In (b), the desired response of protein tetR in the LTD circuit is shown, and in (c), the corresponding response form the HTD circuit is illustrated.

In line with the design proposed by Basu and coworkers, OptCircuit first places the reporter protein under the control of a repressible promoter (i.e., P tet promoter) which is repressed by protein tetR (dotted line in Figure 11(A)). Subsequently, we use OptCircuit to design two circuits, a low threshold detector (LTD) and a high threshold detector (HTD). The LTD circuit expresses high levels of tetR at low levels of L-arabinose and low levels of tetR at high levels of L-arabinose (see Figure 11(B)). In contrast, the HTD circuit is designed to express low levels of tetR at low levels of inducer and high levels of tetR at high levels of L-arabinose (see Figure 11(C)). Finally, the LTD and HTD circuits are fused together to obtain an inverted bell shaped response for protein tetR.

The best circuit configurations proposed by OptCircuit are shown in Figure 12. The only difference between the LTD and HTD is that while tetR is expressed from the Plac promoter in the LTD, it is expressed from PBAD promoter in the HTD. Examination of the circuit behavior reveals that at low levels of L-arabinose, the PBAD promoter is not sufficiently activated ensuring low levels of protein lacI. This in turn implies that the Plac promoter is free to express tetR from the LTD circuit (see Figure 12 and Figure 13(A)). As the amount of L-arabinose accumulates in the system, the transcriptional expression from the PBAD promoter is enhanced leading to expression of lacI from LTD and tetR from HTD (see Figure 12 and 13(B)). Finally, expression of lacI from HTD turns off expression of tetR from HTD. The final OptCircuit design enables the expression of protein tetR from Plac and PBAD promoters, lacI from PBAD promoter and reporter protein GFP from Ptet promoter. The level of protein tetR as a function of level of L-arabinose is shown in Figure 13(C). As shown in Figure 13(C), we find that the circuit response deviates significantly from the desired response implying that by simply reshuffling existing components the desired response is not attainable.

Figure 12
figure 12

The best circuit configurations proposed by OptCircuit. In (A), the LTD circuit is shown; In (B), HTD circuit is illustrated and finally the combination of both the circuit designs is shown in (C).

Figure 13
figure 13

In (A), the response of the LTD circuit is shown. In (B), the response obtained from HTD is illustrated. In (C), the response obtained prior to optimization of kinetic parameters is shown. Note that this response deviates significantly from the desired response. In (D), the response obtained after optimizing kinetic parameters is shown. Specifically the efficiency of constitutive promoter expressing CRP is decreased 13 fold. In all the figures cubic splines were used to interpolate the observed steady-state expression levels.

To address this remaining challenge we next explore whether modifying any existing component in the circuit will shift the circuit response closer to the desired response. Specifically, the circuit described in the previous paragraph is "fixed" and subsequently starting from the current parameter values as an initial guess we optimize the kinetic parameter values using Eq 1.13. Results indicate that a considerable improvement in circuit response is obtained when the transcriptional efficiency of the constitutive promoter expressing protein CRP is decreased 13 fold. This resulted in a 18.69% decrease in the objective value (8.1071 → 6.5976). The effect of this parameter modifications are quantified in Figure 13(D) demonstrating that OptCircuit can be used to pinpoint kinetic parameter modifications improving its functionality.

Discussion

In this work, we introduced an optimization-based approach termed OptCircuit that (i) automatically identifies the circuit components from a list and connectivity that brings about the desired functionality; (ii) Rectify or redesign an existing (non-functional) biological circuit and restore functionality by modifying an existing component and/or identifying additional components to append to the circuit. The dynamics that govern the interactions between the elements of the genetic circuit were modeled using deterministic rate equations and the desired circuit response is abstracted as the maximization/minimization of an appropriately constructed objective function. Subsequently, an iterative procedure was implemented within our framework to identify an ensemble of circuits that exhibit the desired response. The capabilities of the developed tool were investigated by synthesizing circuits that exhibit a wide array of responses. The genetic toggle switch example demonstrated the ability of the framework to suggest simple or more complex circuit configurations capable of discriminating between inducer molecules. The 2 to 4 genetic decoder example led to complex circuit designs consisting of as many as 13 promoter-transcript pairs that may or may not be identifiable through a digital logic based design procedure. Finally, the concentration band detector example illustrated how OptCircuit can be used to design not the architecture of the synthetic circuit but also suggest modifications on its kinetic parameters for optimized performance. OptCircuit can also be employed in tandem with existing computational methods for fine-tuning circuit performance by providing initial configurations. For example, Feng and co-workers developed a global sensitivity analysis based approach to identify the optimal parameter configuration [2]. In their approach, they start from a representative circuit configuration and then proceed to identify the optimal parameter set by estimating the sensitivity of the parameter variation on the circuit response.

It is important to emphasize that all the kinetic parameter modifications suggested by OptCircuit can be realized using a host of experimental strategies. For example in the construction of repressilator circuit, the authors control the rate of protein degradation by ssrA tagging whereby an amino acid sequence is introduced into the proteins which makes them a target for all proteases [37]. Similarly Yokobayashi and co-workers have used directed evolution to restore the performance of an unoptimized circuit [33]. Specifically, by focusing on the cI gene and its corresponding ribosome binding site, the authors report mutations that potentially reduce the translational efficiency or reduce the ribosome binding affinity. Other promising strategies include the approaches developed by Lutz and Bujard [38] to control the promoter activity and repression for the Ptet and Plac promoters.

While our results demonstrate the ability of the framework to design circuits that accomplish a wide variety of responses, it is important to emphasize here that the framework does not take into account potential interactions with other biological components and processes present in the cellular environment. Further, the deterministic nature of OptCircuit may render the circuit deign sensitive to stochastic fluctuations and hence may fail to function properly in wake of noise. This observation assumes further significance in light of multiple modeling studies that have demonstrated the strong influence of noise and stochastic events on circuit performance [15, 39–43]. For example, examination of the decoder circuit response to varying level of input signal reveals that the while the circuit is optimal with respect to connectivity; it is sensitive to changes in the input signal levels (see Additional file 1, Figure S1). This implies that the circuit behavior can be realized in only a very narrow range of concentrations [44] motivating the need to design circuits that are inherently robust to noise and leakiness of specific components. Key features that confer robustness are redundancy, modularity and the ability to decouple perturbations [45]. While most of the current literature regarding biological robustness has focused on elucidating the architectural and mechanistic features of a network, much less effort has been devoted to developing quantitative and qualitative criteria for quantifying robustness. Efforts in this direction include the works of Cherry and Adler [46] and Tomshine and Kaznessis [41] who have proposed that large separation between steady states is likely to render the biological switch immune to stochastic fluctuations. Subsequently, H. El-Farra et al. [47] employed these performance measures to develop optimization problems to identify parameters that confer robustness. The above observations highlight one potential way of improving the robustness of the circuits using OptCircuit. Specifically, performance measures such as separation between steady-states could be imposed as appropriate formulated objective functions to systematically synthesize circuits that are likely to be robust to stochastic fluctuations. Similarly, enforcing alternative ways of realizing an outcome can incorporate other qualitative metrics such as redundancy. We are currently exploring ways of protecting against component failure and incorporating robustness into the circuits we developed in this manuscript.

Another important limitation of the current approach lies in the computational requirements to accomplish a specific circuit design. We find that while OptCircuit was able to readily design circuits with relatively simple architecture (genetic toggle switch), the design procedure for more complex responses (decoder and band detector) entails expensive computational resources. However, our results indicate that the computational burden may be alleviated to a certain extent by exploiting the inherent decomposable structure built into the circuit responses. Specifically, both the decoder and band detector circuits were constructed by assembling smaller less complex building blocks (e.g. LTD and HTD in concentration bad detector). These observations suggest that by carefully identifying key decomposable structures within the desired response, it is possible to overcome this important limitation of the framework.

Conclusion

In recent years, researchers have deposited several standard and interchangeable biological parts in the registry of standard biological parts (e.g. composite parts such as Isoamyl alcohol generating device (BBa_J45400), Elowitz repressilator (BBa_I5610)).

Currently efforts are underway to specify the functionality of these parts in terms of parameter estimates and behavior. Our results and those proposed by other researchers conclusively demonstrate that proper parameter compatibility is essential to ensure funtionality. As the characterization of these parts is moving at a fast pace, the OptCircuit framework could serve as a design platform to aid in the construction and finetuning of integrated biological circuits.

Methods

Modeling framework

The basic elements constituting a genetic circuit include promoter elements, protein/transcript molecules and inducers. Briefly, promoters are regions of DNA where RNA polymerases bind to initiate transcription. Transcripts referred to ORFs which upon transcription and translation produce proteins which in synthetic circuits act as transcriptional regulators repressing or activating a promoter's strength. Finally, inducers are small molecules (e.g., aTC, IPTG) which by directly interacting with transcription factors can block, enable or simply modulate a transcriptional regulation event. The quantitative description of the mechanistic detail underlying the interactions embedded in the genetic circuitry requires the definition of the following sets and variables.

S e t s : V a r i a b l e s : I = { i } = set of promoters J = { j } = set of transcripts K = { k } = set of inducers T = { t } = time P j ( t ) = protein level of transcript  j  at time  t Y i j = { 1  if transcript  j  is expressed from promoter  i 0  otherwise MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiGaaaqaaiabdofatjabdwgaLjabdsha0jabdohaZjabcQda6aqaaiabdAfawjabdggaHjabdkhaYjabdMgaPjabdggaHjabdkgaIjabdYgaSjabdwgaLjabdohaZjabcQda6aabaeqabaGaemysaKKaeyypa0Jaei4EaSNaemyAaKMaeiyFa0Naeyypa0Jaee4CamNaeeyzauMaeeiDaqNaeeiiaaIaee4Ba8MaeeOzayMaeeiiaaIaeeiCaaNaeeOCaiNaee4Ba8MaeeyBa0Maee4Ba8MaeeiDaqNaeeyzauMaeeOCaiNaee4CamhabaGaemOsaOKaeyypa0Jaei4EaSNaemOAaOMaeiyFa0Naeyypa0Jaee4CamNaeeyzauMaeeiDaqNaeeiiaaIaee4Ba8MaeeOzayMaeeiiaaIaeeiDaqNaeeOCaiNaeeyyaeMaeeOBa4Maee4CamNaee4yamMaeeOCaiNaeeyAaKMaeeiCaaNaeeiDaqNaee4CamhabaGaem4saSKaeyypa0Jaei4EaSNaem4AaSMaeiyFa0Naeyypa0Jaee4CamNaeeyzauMaeeiDaqNaeeiiaaIaee4Ba8MaeeOzayMaeeiiaaIaeeyAaKMaeeOBa4MaeeizaqMaeeyDauNaee4yamMaeeyzauMaeeOCaiNaee4CamhabaGaemivaqLaeyypa0Jaei4EaSNaeeiDaqNaeiyFa0Naeyypa0JaeeiDaqNaeeyAaKMaeeyBa0MaeeyzaugaaqaabeqaaiabdcfaqnaaBaaaleaacqWGQbGAaeqaaOGaeiikaGIaemiDaqNaeiykaKIaeyypa0JaeeiCaaNaeeOCaiNaee4Ba8MaeeiDaqNaeeyzauMaeeyAaKMaeeOBa4MaeeiiaaIaeeiBaWMaeeyzauMaeeODayNaeeyzauMaeeiBaWMaeeiiaaIaee4Ba8MaeeOzayMaeeiiaaIaeeiDaqNaeeOCaiNaeeyyaeMaeeOBa4Maee4CamNaee4yamMaeeOCaiNaeeyAaKMaeeiCaaNaeeiDaqNaeeiiaaIaemOAaOMaeeiiaaIaeeyyaeMaeeiDaqNaeeiiaaIaeeiDaqNaeeyAaKMaeeyBa0MaeeyzauMaeeiiaaIaemiDaqhabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpdaGabaqcfayaauaabaqaceaaaeaacqaIXaqmcqqGGaaicqqGPbqAcqqGMbGzcqqGGaaicqqG0baDcqqGYbGCcqqGHbqycqqGUbGBcqqGZbWCcqqGJbWycqqGYbGCcqqGPbqAcqqGWbaCcqqG0baDcqqGGaaicqWGQbGAcqqGGaaicqqGPbqAcqqGZbWCcqqGGaaicqqGLbqzcqqG4baEcqqGWbaCcqqGYbGCcqqGLbqzcqqGZbWCcqqGZbWCcqqGLbqzcqqGKbazcqqGGaaicqqGMbGzcqqGYbGCcqqGVbWBcqqGTbqBcqqGGaaicqqGWbaCcqqGYbGCcqqGVbWBcqqGTbqBcqqGVbWBcqqG0baDcqqGLbqzcqqGYbGCcqqGGaaicqWGPbqAaeaacqaIWaamcqqGGaaicqqGVbWBcqqG0baDcqqGObaAcqqGLbqzcqqGYbGCcqqG3bWDcqqGPbqAcqqGZbWCcqqGLbqzaaaakiaawUhaaaaaaaaa@2B3E@

The set I represents all the promoter elements investigated in this study. J represents the set of transcripts and finally K is the set of all inducer molecules. Model variables encode the structure of the synthetic circuit and quantify the protein levels. Specifically, the binary variable Y ij determines which transcript j is expressed from a promoter i and P j (t) quantifies the level of protein j at a given time t.

Kinetic description of interactions

Genetic circuits are characterized by a number of interactions including protein-promoter and protein-inducer and protein-protein interactions. For example, protein lacI, in its tetramer form functions as a repressor for Plac promoter while inducer molecule aTc suppresses the activity of protein tetR. In genetic circuits, unlike digital or binary logic based circuits, the presence/absence of a particular set of interactions alone is insufficient to accurately predict correctly all possible responses. In fact, several studies have reported that in addition to interactions, the kinetic rates of individual elements have to accurately match in order to ensure function. To this end, the kinetic description of each element of a genetic circuit is embedded into the OptCircuit framework.

Specifically, for every transcript j, the set of ordinary differential equations (ODE's) that govern the time evolution of the protein is given by Eq 1.2.

d P j d t = ∑ i Y i j [ Rate of Production of j from i ] − K d e c a y j P j ( t ) ∀ j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGqbaudaWgaaqaaiabdQgaQbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaabuaeaacqWGzbqwdaWgaaWcbaGaemyAaKMaemOAaOgabeaakiabcUfaBjabbkfasjabbggaHjabbsha0jabbwgaLjabbccaGiabb+gaVjabbAgaMjabbccaGiabbcfaqjabbkhaYjabb+gaVjabbsgaKjabbwha1jabbogaJjabbsha0jabbMgaPjabb+gaVjabb6gaUjabbccaGiabb+gaVjabbAgaMjabbccaGiabbQgaQjabbccaGiabbAgaMjabbkhaYjabb+gaVjabb2gaTjabbccaGiabbMgaPjabc2faDjabgkHiTiabdUealnaaDaaaleaacqWGKbazcqWGLbqzcqWGJbWycqWGHbqycqWG5bqEaeaacqWGQbGAaaaabaGaemyAaKgabeqdcqGHris5aOGaemiuaa1aaSbaaSqaaiabdQgaQbqabaGccqGGOaakcqWG0baDcqGGPaqkcqGHaiIicqWGQbGAaaa@7626@
(1.2)

The first term in Eq 1.2, accounts for the cumulative rate of production of a particular protein j from the promoter elements and the second term represents the first order decay of the protein. Also observe that the production of a protein j from a promoter i is turned ON if and only if the corresponding binary variable Y ij is equal to one.

OptCircuit accounts for the activating and repressing effects on every promoter i within the framework by using the modeling formulation proposed by Hasty et al. [25]. Briefly, all biochemical reactions characterizing the interactions affecting a particular promoter are listed and divided into fast and slow steps. The fast reaction set typically includes protein dimerization and protein promoter binding while transcription and degradation steps compose the slow reaction set. The dynamics governing the promoter kinetics are derived using mass action kinetics with fast reactions that have rate constants in the order of seconds, assumed to be in equilibrium [25]. The modeling environment in OptCircuit is versatile enough to incorporate finer levels of mechanistic detail whenever available (e.g., modeling of mRNA [7]). The complete list of kinetic formulations for the promoter elements and the corresponding kinetic parameters employed in this study are provided in the supplementary material (see Additional file 1).

Objective Function Modeling

The reliance on an optimization framework for designing synthetic circuits with a desired response implies that the objective function must be carefully chosen so as its maximization or minimization is a good surrogate of the desired response(s). The type of desired responses is partitioned into inducer-free and inducer-dependent ones. Inducer-free responses translate into the design of circuits whose response is consistent with a targeted time-course. This response may be oscillatory, constant or ramping up/down. For all these case the objective function, Z minimizes the sum of the squared departures from the targeted responses at all time points:

M i n i m i z e  Z = ∑ t ( P j ∗ ( t ) − P j ∗ exp ( t ) ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyta0KaemyAaKMaemOBa4MaemyAaKMaemyBa0MaemyAaKMaemOEaONaemyzauMaeeiiaaIaeeOwaOLaeyypa0ZaaabuaeaacqGGOaakcqWGqbaudaWgaaWcbaGaemOAaO2aaWbaaWqabeaacqGHxiIkaaaaleqaaOGaeiikaGIaemiDaqNaeiykaKIaeyOeI0Iaemiuaa1aa0baaSqaaiabdQgaQnaaCaaameqabaGaey4fIOcaaaWcbaGagiyzauMaeiiEaGNaeiiCaahaaOGaeiikaGIaemiDaqNaeiykaKIaeiykaKYaaWbaaSqabeaacqaIYaGmaaaabaGaemiDaqhabeqdcqGHris5aaaa@5381@
(1.3)

In Equation 1.3, P j ∗ exp ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiuaa1aa0baaSqaaiabdQgaQnaaCaaameqabaGaey4fIOcaaaWcbaGagiyzauMaeiiEaGNaeiiCaahaaOGaeiikaGIaemiDaqNaeiykaKcaaa@3711@ denotes the experimentally observed profile. Inducer-dependent responses require a clear distinction between states corresponding to presence/absence of multiple inducers. Specifically, to accomplish this, an objective function is constructed as follows that maximizes the scaled separation between the inducer-present/absent responses:

M a x i m i z e  Z = ∑ j ∈ R ∑ k ∈ K P j ( T ) | k − ∑ j ′ ≠ j ∈ R P j ′ ( T ) | k P j ( T ) | k MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyta0KaemyyaeMaemiEaGNaemyAaKMaemyBa0MaemyAaKMaemOEaONaemyzauMaeeiiaaIaeeOwaOLaeyypa0ZaaabuaeaadaaeqbqcfayaamaalaaabaWaaqGaaeaacqWGqbaudaWgaaqaaiabdQgaQbqabaGaeiikaGIaemivaqLaeiykaKcacaGLiWoadaWgaaqaaiabdUgaRbqabaGaeyOeI0YaaabuaeaadaabcaqaaiabdcfaqnaaBaaabaGafmOAaOMbauaaaeqaaiabcIcaOiabdsfaujabcMcaPaGaayjcSdWaaSbaaeaacqWGRbWAaeqaaaqaaiqbdQgaQzaafaGaeyiyIKRaemOAaOMaeyicI4SaemOuaifabeGaeyyeIuoaaeaacqWGqbaudaWgaaqaaiabdQgaQbqabaGaeiikaGIaemivaqLaeiykaKIaeiiFaW3aaSbaaeaacqWGRbWAaeqaaaaaaSqaaiabdUgaRjabgIGiolabdUealbqab0GaeyyeIuoaaSqaaiabbQgaQjabgIGiolabbkfasbqab0GaeyyeIuoaaaa@6AF7@
(1.4)

In Equation 1.4, K represents the set of inducer molecules present in the system, R represents the set of reporter proteins (e.g. GFP, YFP etc) and P j (T)| k represents the steady state levels of transcript j in presence of inducer k. Alternatively, if the circuit response must be inducer concentration dependent then the objective function can be formulated again as a the minimization of a least squares sum by considering multiple inducer concentrations.

M i n i m i z e  Z = ∑ r ( P j | k , r − P j exp | k , r ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyta0KaemyAaKMaemOBa4MaemyAaKMaemyBa0MaemyAaKMaemOEaONaemyzauMaeeiiaaIaeeOwaOLaeyypa0ZaaabuaeaadaqadaqaamaaeiaabaGaemiuaa1aaSbaaSqaaiabdQgaQbqabaaakiaawIa7amaaBaaaleaacqWGRbWAcqGGSaalcqWGYbGCaeqaaOGaeyOeI0YaaqGaaeaacqWGqbaudaqhaaWcbaGaemOAaOgabaGagiyzauMaeiiEaGNaeiiCaahaaaGccaGLiWoadaWgaaWcbaGaem4AaSMaeiilaWIaemOCaihabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaaqaaiabdkhaYbqab0GaeyyeIuoaaaa@55AE@
(1.5)

In Eq 1.5, r represents the discretizations levels for the inducer concentration. P j |k, rand P j exp | k , r MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaaqGaaeaacqWGqbaudaqhaaWcbaGaemOAaOgabaGagiyzauMaeiiEaGNaeiiCaahaaaGccaGLiWoadaWgaaWcbaGaem4AaSMaeiilaWIaemOCaihabeaaaaa@3834@ represent the simulated and the desired steady-state levels of reporter transcript j at inducer discretization level r. These are only some examples of desired circuit responses. Using this optimization-based framework even more complicated responses can be modeled limited only by the imagination of the circuit designer. In the examples investigated in this study Eq 1.4 represents the objective function imposed for the genetic decoder example and Eq 1.5, for the concentration band detector example.

Optimization model

Using the notation listed above, the problem of designing a genetic circuit that exhibits a desired response is formulated as the following mixed integer dynamic optimization problem (MIDO) [48–52].

Min/Max Z = f (P j (t))

s.t.

d P j d t = ∑ i Y i j [ Rate of Production of j from i ] − K d e c a y j P j ( t ) ∀ j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGqbaudaWgaaqaaiabdQgaQbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaabuaeaacqWGzbqwdaWgaaWcbaGaemyAaKMaemOAaOgabeaakiabcUfaBjabbkfasjabbggaHjabbsha0jabbwgaLjabbccaGiabb+gaVjabbAgaMjabbccaGiabbcfaqjabbkhaYjabb+gaVjabbsgaKjabbwha1jabbogaJjabbsha0jabbMgaPjabb+gaVjabb6gaUjabbccaGiabb+gaVjabbAgaMjabbccaGiabbQgaQjabbccaGiabbAgaMjabbkhaYjabb+gaVjabb2gaTjabbccaGiabbMgaPjabc2faDjabgkHiTiabdUealnaaDaaaleaacqWGKbazcqWGLbqzcqWGJbWycqWGHbqycqWG5bqEaeaacqWGQbGAaaaabaGaemyAaKgabeqdcqGHris5aOGaemiuaa1aaSbaaSqaaiabdQgaQbqabaGccqGGOaakcqWG0baDcqGGPaqkcqGHaiIicqWGQbGAaaa@7626@
(1.7)
∑ j Y i j ≤ P max ∀ i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGqbaudaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGQbGAaeqaniabggHiLdaakeaacqGHaiIicqWGPbqAaaaaaa@3D34@
(1.8)
∑ i Y i j ≤ T max ∀ j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGubavdaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGPbqAaeqaniabggHiLdaakeaacqGHaiIicqWGQbGAaaaaaa@3D3C@
(1.9)
∑ i ∑ j Y i j ≤ M M a x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeqaaaqaamaaqafabaWaaabuaeaacqWGzbqwdaWgaaWcbaGaemyAaKMaemOAaOgabeaakiabgsMiJkabd2eannaaBaaaleaacqWGnbqtcqWGHbqycqWG4baEaeqaaaqaaiabdQgaQbqab0GaeyyeIuoaaSqaaiabdMgaPbqab0GaeyyeIuoaaaaaaa@3E35@
(1.10)

The objective function in Equation 1.6 models the circuit response imposed by the circuit designer. Equation 1.7, describes the time evolution of protein levels as a set of ordinary differential equations as described in the previous section. Equation 1.8 imposes an upper limit on the number of transcripts a particular promoter i can express. Similarly, Equation 1.9 imposes a limit on the number of times a particular transcript j can be expressed from different promoters. Finally, Equation 1.10 imposes a limit on the total number of promoter-transcript pairs in the designed genetic circuit.

The boolean constraints (Equations 1.8–1.10) offer the flexibility to incorporate the design of an existing biological circuit and probe its behavior. This can be accomplished by incorporating constraints of the form

Y ij = 1 ∀ ({i, j} ∈ EX)

where, the set EX contains the connectivity information of the circuit. This feature confers upon us the ability to readily extend the framework to rectify or redesign an existing (non functional) biological circuit by identifying additional components to append to the circuit to ensure its functionality.

The solution procedure for the MIDO class of optimization problems is difficult [52] due to the simultaneous presence of binary variables Y ij and constraints in form of ODE's. Reliable solution methodologies that guarantee a global optimal solution for this class of problems are still in infancy [53]. Therefore, in this paper we rely on a decomposition procedure to bracket an optimal solution. The basic idea of proposed approach is to generate a converging sequence of upper and lower bounds to the original problem. The solution procedure is listed in a step-wise manner below.

Step 1: Initialize iteration counter, iter = 1; SET upper bound UB = ∞; SET lower bound LB = -∞; Generate an initial guess for a feasible circuit design Y i j = Y i j i t e r MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpcqWGzbqwdaqhaaWcbaGaemyAaKMaemOAaOgabaGaemyAaKMaemiDaqNaemyzauMaemOCaihaaaaa@3AB0@ . In this implementation of the framework, the initial guess is generated by simply solving the following optimization problem with a objective function set to zero.

Min Z = 0

s.t.

∑ j Y i j ≤ P max ∀ i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGqbaudaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGQbGAaeqaniabggHiLdaakeaacqGHaiIicqWGPbqAaaaaaa@3D34@
(1.8a)
∑ i Y i j ≤ T max ∀ j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGubavdaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGPbqAaeqaniabggHiLdaakeaacqGHaiIicqWGQbGAaaaaaa@3D3C@
(1.9a)
∑ i ∑ j Y i j ≤ M M a x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeqaaaqaamaaqafabaWaaabuaeaacqWGzbqwdaWgaaWcbaGaemyAaKMaemOAaOgabeaakiabgsMiJkabd2eannaaBaaaleaacqWGnbqtcqWGHbqycqWG4baEaeqaaaqaaiabdQgaQbqab0GaeyyeIuoaaSqaaiabdMgaPbqab0GaeyyeIuoaaaaaaa@3E35@
(1.10a)

After convergence to local optima, this initial guess is excluded by using integer cuts and the above problem is used to generate the next starting point.

Step 2: Integrate the system of ordinary differential equations (1.2) for fixed values of the design variables Y i j = Y i j i t e r MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpcqWGzbqwdaqhaaWcbaGaemyAaKMaemOAaOgabaGaemyAaKMaemiDaqNaemyzauMaemOCaihaaaaa@3AB0@ to obtain the objective function value Z; SET UB = min(UB,Z). Store the solution corresponding to the best upper bound.

Step 3: In step2, we compute an upper bound for the objective function value. In step3 we compute a lower bound. The main ideas involved in the computation of the lower bound are presented in the development of the outer approximation procedure for solving non-linear and mixed integer non-linear optimization algorithms [54]. Briefly, the Master problem described in step 3 of the optimization procedure constructs a relaxation of the original feasible space by constructing supporting hyper planes at the point of interest. It has been shown previously that the solution to the Master problem yields a lower bound to the objective function value [54]. In our case the supporting hyper planes are constructed at the integer solution(s) obtained from step 2 and the partial derivates are computed by finite difference method. The master problem yields a lower bound on the objective function value and a new integer solution. A brief description of the main ideas governing the outer approximation procedure is provided in supplementary material (see Additional file 1).

Construct the master (lower bounding) problem as follows.

min imize μ

s.t.

μ ≥ Z i t e r + ∑ i ∑ j ( Y i j − Y i j k ) ( ∂ Z ∂ Y i j ) Y i j = Y i j k ∀ k = 1 , 2 , ... i t e r MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaaiabeY7aTjabgwMiZkabdQfaAnaaCaaaleqabaGaemyAaKMaemiDaqNaemyzauMaemOCaihaaOGaey4kaSYaaabuaeaadaaeqbqaamaabmaabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHsislcqWGzbqwdaqhaaWcbaGaemyAaKMaemOAaOgabaGaem4AaSgaaaGccaGLOaGaayzkaaWaaeWaaKqbagaadaWcaaqaaiabgkGi2kabdQfaAbqaaiabgkGi2kabdMfaznaaBaaabaGaemyAaKMaemOAaOgabeaaaaaakiaawIcacaGLPaaadaWgaaWcbaGaemywaK1aaSbaaWqaaiabdMgaPjabdQgaQbqabaWccqGH9aqpcqWGzbqwdaqhaaadbaGaemyAaKMaemOAaOgabaGaem4AaSgaaaWcbeaaaeaacqWGQbGAaeqaniabggHiLdaaleaacqWGPbqAaeqaniabggHiLdaakeaacqGHaiIicqWGRbWAcqGH9aqpcqaIXaqmcqGGSaalcqaIYaGmcqGGSaalcqGGUaGlcqGGUaGlcqGGUaGlcqWGPbqAcqWG0baDcqWGLbqzcqWGYbGCaaaaaa@6EF8@
(1.12)
∑ j Y i j ≤ P max ∀ i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGqbaudaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGQbGAaeqaniabggHiLdaakeaacqGHaiIicqWGPbqAaaaaaa@3D34@
(1.8b)
∑ i Y i j ≤ T max ∀ j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaamaaqafabaGaemywaK1aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGHKjYOcqWGubavdaahaaWcbeqaaiGbc2gaTjabcggaHjabcIha4baaaeaacqWGPbqAaeqaniabggHiLdaakeaacqGHaiIicqWGQbGAaaaaaa@3D3C@
(1.9b)
∑ i ∑ j Y i j ≤ M M a x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeqaaaqaamaaqafabaWaaabuaeaacqWGzbqwdaWgaaWcbaGaemyAaKMaemOAaOgabeaakiabgsMiJkabd2eannaaBaaaleaacqWGnbqtcqWGHbqycqWG4baEaeqaaaqaaiabdQgaQbqab0GaeyyeIuoaaSqaaiabdMgaPbqab0GaeyyeIuoaaaaaaa@3E35@
(1.10b)

The partial derivates are computed using finite difference methods.

Solve the master problem to obtain the objective function value μ* and integer solution, Y i j ∗ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemywaK1aa0baaSqaaiabdMgaPjabdQgaQbqaaiabgEHiQaaaaaa@30E4@ .

SET LB = μ*

Step 4: If LB ≥ UB, then STOP (crossover). Otherwise, increase iteration counter iter → iter + 1. Y i j i t e r = Y i j ∗ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemywaK1aa0baaSqaaiabdMgaPjabdQgaQbqaaiabdMgaPjabdsha0jabdwgaLjabdkhaYbaakiabg2da9iabdMfaznaaDaaaleaacqWGPbqAcqWGQbGAaeaacqGHxiIkaaaaaa@3BA0@ . Return to Step 2.

In addition to identifying the optimal configuration of design variables ( Y i j ∗ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemywaK1aa0baaSqaaiabdMgaPjabdQgaQbqaaiabgEHiQaaaaaa@30E4@ ), OptCircuit can also be employed to optimize kinetic parameters of specific elements within the genetic circuit. For example, given a genetic circuit, the task of determining the optimal promoter strength of a particular promoter i* can be achieved by replacing Equation 1.12 with

μ ≥ Z i t e r + ( α i − α i k ) ( ∂ Z ∂ α i ) α i = α j k ∀ k = 1 , 2 , ... i t e r MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaaiabeY7aTjabgwMiZkabdQfaAnaaCaaaleqabaGaemyAaKMaemiDaqNaemyzauMaemOCaihaaOGaey4kaSYaaeWaaeaacqaHXoqydaWgaaWcbaGaemyAaKgabeaakiabgkHiTiabeg7aHnaaDaaaleaacqWGPbqAaeaacqWGRbWAaaaakiaawIcacaGLPaaadaqadaqcfayaamaalaaabaGaeyOaIyRaemOwaOfabaGaeyOaIyRaeqySde2aaSbaaeaacqWGPbqAaeqaaaaaaOGaayjkaiaawMcaamaaBaaaleaacqaHXoqydaWgaaadbaGaemyAaKgabeaaliabg2da9iabeg7aHnaaDaaameaacqWGQbGAaeaacqWGRbWAaaaaleqaaaGcbaGaeyiaIiIaem4AaSMaeyypa0JaeGymaeJaeiilaWIaeGOmaiJaeiilaWIaeiOla4IaeiOla4IaeiOla4IaemyAaKMaemiDaqNaemyzauMaemOCaihaaaaa@632A@
(1.13)

Note that given the nonlinear nature of the problem under investigation the above procedure is carried out multiple times starting for several starting initial guesses and the local optimum solution identified at each iteration is stored along with a sorted list of the best circuit configurations.

References

  1. Sprinzak D, Elowitz MB: Reconstruction of genetic circuits. Nature. 2005, 438: 443-448. 10.1038/nature04335

    Article  CAS  PubMed  Google Scholar 

  2. Feng XJ, Hooshangi S, Chen D, Li G, Weiss R: Optimizing genetic circuits by global sensitivity analysis. Biophys J. 2004, 87: 2195-2202. 10.1529/biophysj.104.044131

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Hasty J, McMillen D, Collins JJ: Engineered gene circuits. Nature. 2002, 420: 224-230. 10.1038/nature01257

    Article  CAS  PubMed  Google Scholar 

  4. Endy D: Foundations for engineering biology. Nature. 2005, 438: 449-453. 10.1038/nature04342

    Article  CAS  PubMed  Google Scholar 

  5. Alon U: Biological networks: the tinkerer as an engineer. Science. 2003, 301: 1866-1867. 10.1126/science.1089072

    Article  CAS  PubMed  Google Scholar 

  6. Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403: 339-342. 10.1038/35002131

    Article  CAS  PubMed  Google Scholar 

  7. Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403: 335-338. 10.1038/35002125

    Article  CAS  PubMed  Google Scholar 

  8. McMillen D, Kopell N, Hasty J, Collins JJ: Synchronizing genetic relaxation oscillators by intercell signaling. Proc Natl Acad Sci USA. 2002, 99: 679-684. 10.1073/pnas.022642299

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Fung E, Wong WW, Suen JK, Bulter T, Lee SG: A synthetic gene-metabolic oscillator. Nature. 2005, 435: 118-122. 10.1038/nature03508

    Article  CAS  PubMed  Google Scholar 

  10. Atkinson MR, Savageau MA, Myers JT, Ninfa AJ: Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli. Cell. 2003, 113: 597-607. 10.1016/S0092-8674(03)00346-5

    Article  CAS  PubMed  Google Scholar 

  11. Basu S, Gerchman Y, Collins CH, Arnold FH, Weiss R: A synthetic multicellular system for programmed pattern formation. Nature. 2005, 434: 1130-1134. 10.1038/nature03461

    Article  CAS  PubMed  Google Scholar 

  12. Hasty J, Pradines J, Dolnik M, Collins JJ: Noise-based switches and amplifiers for gene expression. Proc Natl Acad Sci USA. 2000, 97: 2075-2080. 10.1073/pnas.040411297

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Judd EM, Laub MT, McAdams HH: Toggles and oscillators: new genetic circuit designs. Bioessays. 2000, 22: 507-509. 10.1002/(SICI)1521-1878(200006)22:6<507::AID-BIES3>3.0.CO;2-0

    Article  CAS  PubMed  Google Scholar 

  14. Kobayashi H, Kaern M, Araki M, Chung K, Gardner TS: Programmable cells: interfacing natural and engineered gene networks. Proc Natl Acad Sci USA. 2004, 101: 8414-8419. 10.1073/pnas.0402940101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Hooshangi S, Thiberge S, Weiss R: Ultrasensitivity and noise propagation in a synthetic transcriptional cascade. Proc Natl Acad Sci USA. 2005, 102: 3581-3586. 10.1073/pnas.0408507102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Mangan S, Zaslaver A, Alon U: The coherent feedforward loop serves as a sign-sensitive delay element in transcription networks. J Mol Biol. 2003, 334: 197-204. 10.1016/j.jmb.2003.09.049

    Article  CAS  PubMed  Google Scholar 

  17. Mangan S, Alon U: Structure and function of the feed-forward loop network motif. Proc Natl Acad Sci USA. 2003, 100: 11980-11985. 10.1073/pnas.2133841100

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Becskei A, Serrano L: Engineering stability in gene networks by autoregulation. Nature. 2000, 405: 590-593. 10.1038/35014651

    Article  CAS  PubMed  Google Scholar 

  19. Martin VJ, Pitera DJ, Withers ST, Newman JD, Keasling JD: Engineering a mevalonate pathway in Escherichia coli for production of terpenoids. Nat Biotechnol. 2003, 21: 796-802. 10.1038/nbt833

    Article  CAS  PubMed  Google Scholar 

  20. Anderson JC, Clarke EJ, Arkin AP, Voigt CA: Environmentally controlled invasion of cancer cells by engineered bacteria. J Mol Biol. 2006, 355: 619-627. 10.1016/j.jmb.2005.10.076

    Article  CAS  PubMed  Google Scholar 

  21. Levskaya A, Chevalier AA, Tabor JJ, Simpson ZB, Lavery LA: Synthetic biology: engineering Escherichia coli to see light. Nature. 2005, 438: 441-442. 10.1038/nature04405

    Article  CAS  PubMed  Google Scholar 

  22. Tian J, Gong H, Sheng N, Zhou X, Gulari E: Accurate multiplex gene synthesis from programmable DNA microchips. Nature. 2004, 432: 1050-1054. 10.1038/nature03151

    Article  CAS  PubMed  Google Scholar 

  23. Website title., http://parts.mit.edu

  24. Tyson JJ, Chen KC, Novak B: Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr Opin Cell Biol. 2003, 15: 221-231. 10.1016/S0955-0674(03)00017-6

    Article  CAS  PubMed  Google Scholar 

  25. Hasty J, Isaacs F, Dolnik M, McMillen D, Collins JJ: Designer gene networks: Towards fundamental cellular control. Chaos. 2001, 11: 207-220. 10.1063/1.1345702

    Article  CAS  PubMed  Google Scholar 

  26. Gilman A, Arkin AP: Genetic "code": representations and dynamical models of genetic components and networks. Annu Rev Genomics Hum Genet. 2002, 3: 341-369. 10.1146/annurev.genom.3.030502.111004

    Article  CAS  PubMed  Google Scholar 

  27. Glass L, Perkins TJ, Mason J, Siegelmann HT, Edwards R: Chaotic Dynamics in an Electronic Model of a Genetic Network. Journal of Statistical Physics. 2005, 121: 969-994. 10.1007/s10955-005-7009-y.

    Article  Google Scholar 

  28. Basu S, Weiss R: The Device physics of Cellular Logic Gates. 2002, 54-61.

    Google Scholar 

  29. Weiss R, Homay GE, Knight TF: Toward in vivo Digital Circuits. 1999

    Google Scholar 

  30. Mason J, Linsay PS, Collins JJ, Glass L: Evolving complex dynamics in electronic models of genetic networks. Chaos. 2004, 14: 707-715. 10.1063/1.1786683

    Article  PubMed  Google Scholar 

  31. Guet CC, Elowitz MB, Hsing W, Leibler S: Combinatorial synthesis of genetic networks. Science. 2002, 296: 1466-1470. 10.1126/science.1067407

    Article  CAS  PubMed  Google Scholar 

  32. Tuttle LM, Salis H, Tomshine J, Kaznessis YN: Model-driven designs of an oscillating gene network. Biophys J. 2005, 89: 3873-3883. 10.1529/biophysj.105.064204

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Yokobayashi Y, Weiss R, Arnold FH: Directed evolution of a genetic circuit. Proc Natl Acad Sci USA. 2002, 99: 16587-16591. 10.1073/pnas.252535999

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Francois P, Hakim V: Design of genetic networks with specified functions by evolution in silico. Proc Natl Acad Sci USA. 2004, 101: 580-585. 10.1073/pnas.0304532101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Battogtokh D, Asch DK, Case ME, Arnold J, Schuttler HB: An ensemble method for idenitfying regulatory circuits with special reference to the qa gene cluster of Neurospora crassa. Proc Natl Acad Sci USA. 2002, 99: 16904-16909. 10.1073/pnas.262658899

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Basu S, Karig D, Weiss R: Engineering signal processing in cells: Towards molecular concentration band detection. Natural Computing. 2003, 2: 463-478. 10.1023/B:NACO.0000006774.27778.f0.

    Article  CAS  Google Scholar 

  37. Keiler KC, Waller PR, Sauer RT: Role of a peptide tagging system in degradation of proteins synthesized from damaged messenger RNA. Science. 1996, 271: 990-993. 10.1126/science.271.5251.990

    Article  CAS  PubMed  Google Scholar 

  38. Lutz R, Bujard H: Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1–I2 regulatory elements 10.1093/nar/25.6.1203. Nucl Acids Res. 1997, 25: 1203-1210.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Hooshangi S, Weiss R: The effect of negative feedback on noise propagation in transcriptional gene networks. Chaos. 2006, 16: 26108-26108. 10.1063/1.2208927.

    Article  Google Scholar 

  40. Salis H, Kaznessis YN: Computer-aided design of modular protein devices: Boolean AND gene activation. Phys Biol. 2006, 3: 295-310. 10.1088/1478-3975/3/4/007

    Article  CAS  PubMed  Google Scholar 

  41. Tomshine J, Kaznessis YN: Optimization of a stochastically simulated gene network model via simulated annealing. Biophys J. 2006, 91: 3196-3205. 10.1529/biophysj.106.083485

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Sotiropoulos V, Kaznessis YN: Synthetic tetracycline-inducible regulatory networks: computer-aided design of dynamic phenotypes. BMC Syst Biol. 2007, 1: 7- 10.1186/1752-0509-1-7

    Article  PubMed Central  PubMed  Google Scholar 

  43. Zak DE, Gonye GE, Schwaber JS, Doyle FJ: Importance of input perturbations and stochastic gene expression in the reverse engineering of genetic regulatory networks: insights from an identifiability analysis of an in silico network. Genome Res. 2003, 13: 2396-2405. 10.1101/gr.1198103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Kramer BP, Weber W, Fussenegger M: Artificial regulatory networks and cascades for discrete multilevel transgene control in mammalian cells. Biotechnol Bioeng. 2003, 83: 810-820. 10.1002/bit.10731

    Article  CAS  PubMed  Google Scholar 

  45. Kitano H: Biological robustness. Nature Reviews Genetics. 2004, 5: 826-837. 10.1038/nrg1471

    Article  CAS  PubMed  Google Scholar 

  46. Cherry JL, Adler FR: How to make a biological switch. J Theor Biol. 2000, 203: 117-133. 10.1006/jtbi.2000.1068

    Article  CAS  PubMed  Google Scholar 

  47. El-Farra NH: An Optimization-Based Method for the Design of Robust Synthetic Switches in Biological Networks. Cincinnati. 2005

    Google Scholar 

  48. Mohideen MJ, Perkins JD, Pistikopoulos EN: Towards an Efficient Numerical Procedure for Mixed Integer Optimal Control. Computers and Chemical Engineering. 1997, 21: S457-S462.

    Article  CAS  Google Scholar 

  49. Sirdeshpande AR, Ierapetritou MG, Androulakis IP: Design of Flexible Reduced Kinetic Mechanisms. Process Systems Engineering. 2001, 47: 2461-2473.

    CAS  Google Scholar 

  50. Tlacuahuac AF, Beigler LT: A Robust and Efficient Mixed-Integer Non-Linear Dynamic Optimization Approach for Simultaneous Design and Control. 2004

    Google Scholar 

  51. Bansal V, Perkins JD, Pistikopoulos EN: A Case Study in Simultaneous Design and Control Using Rigorous Mixed-Integer Dynamic Optimization Models. IndEngChemRes. 2002, 41: 760-778.

    CAS  Google Scholar 

  52. Bansal V, Sakizlis V, Ross R, Perkins JD, Pistikopoulos EN: New algorithms for mixed-integer dynamic optimization. Computers and Chemical Engineering. 2003, 27: 647-688. 10.1016/S0098-1354(02)00261-2.

    Article  CAS  Google Scholar 

  53. Chachuat B, Singer AB, Barton PI: Global Mixed-Integer Dynamic optimization. AIChe Journal. 2005, 51: 2235-2253. 10.1002/aic.10494.

    Article  CAS  Google Scholar 

  54. Floudas CA: Nonlinear and Mixed-Integer Optimization. 1995, Fundamentals and Applications: Oxford University Press

    Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Patrick Cirino for helpful discussions and carefully reading the manuscript. The authors gratefully acknowledge funding from the department of energy grant DOE: DE-FG02-05ER25684.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Costas D Maranas.

Additional information

Authors' contributions

MSD and CDM designed and developed the computational algorithms. MSD performed the experiments. MSD and CDM wrote the manuscript. Both the authors have read and approved the final manuscript.

Electronic supplementary material

12918_2007_181_MOESM1_ESM.doc

Additional file 1: Supplementary material. This file provides the following information. i) List of promoters, transcripts and inducer molecules employed in this study. ii) A description of the activating and inhibiting interactions employed in this work. iii) The set of ordinary differential equations used for the genetic toggle switch example. In addition a brief description of the mechanistic detail embedded in these equations is provided. iv) List of nominal parameter values that were used for the genetic toggle switch example (1st example). v) Equations describing the production terms for the promoters used for the genetic decoder and concentration band detector examples. vi) List of nominal parameter values that were used for the genetic decoder and concentration band detector examples (2nd/3rd examples). vii) The sensitivity of the reporter proteins to varying levels of input signals in the genetic decoder example. viii) Legend describing the representation adopted for the logic gates. ix) A brief description of the main ideas behind the outer approximation procedure. (DOC 336 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Dasika, M.S., Maranas, C.D. OptCircuit: An optimization based method for computational design of genetic circuits. BMC Syst Biol 2, 24 (2008). https://doi.org/10.1186/1752-0509-2-24

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1752-0509-2-24

Keywords