Dynamic Analysis of Interlocked Positive Feedback Loops
|
|
- Joan Gilbert
- 6 years ago
- Views:
Transcription
1 Journal of the Korean Physical Society, Vol. 53, No. 6, December 2008, pp Dynamic Analysis of Interlocked Positive Feedback Loops Cuong Nguyen and Seung Kee Han Physics Department, Chungbuk National University, Cheongju (Received 3 September 2008) Positive feedback loops, which are responsible for reliable switches in biological systems, not only work alone but are also interlocked together. In this paper, we analyzed the possible dynamics of the interlocking of two dierent types of positive feedback loops, self-enhanced and self-recovering positive feedback loops. Using a bifurcation analysis, we showed that the interlocked positive feedback loop was able to generate diverse dynamical behaviors. They are classied as monostability, single bistability, non-overlapping, overlapping, nested and merged bistability with or without irreversibility. The diversity of dynamical behaviors implies the exibility of biological systems in response to external stimuli. PACS numbers: Aa, Yc, Xa, d Keywords: Positive feedback loop, interlocked feedback loops, Bifurcation analysis, Dynamical diversity, Multistability I. INTRODUCTION Biological complex networks possess many simple modules with specic functions [1, 2]. These modules are built up from many building blocks with specic dynamical properties [3]. Positive feedback loops (PFLs), which are one of the building blocks, are repeatedly reported in many biological systems, such as cell cycles [3{6] and circadian rhythms [7]. PFLs have potential to exhibit bistability and irreversibility and play an important role in biological systems by converting a graded signal into a discontinuous response with hysteresis or reliable switches. The PFLs work in an interlocked manner, rather than alone. Interlocked PFLs can be found widely in many biological complex networks, such as mammalian cell cycle regulatory networks [6, 8, 9], transcriptional networks [10], gene regulation networks, etc. For example, in the budding yeast cell cycle model [4], the coupling of S and M phases is governed by interlocking of two PFLs of mutual inhibition. The inhibitor CKI mutually inhibits Clb5 6 and Clb1 2; Clb5 6 and Clb1 2 can, in turn, degrade CKI through uibiquitination [4]. Other interlocking modules can also be found in the ssion yeast cell cycle model [5], in which the whole cell cycle model is divided into three modules: G1=S, G2=M and Finish (M=G1). The cyclin-dependent kinase complex MP F is a unique regulator of those modules. Module G1=S is composed by interlocking two positive feedback loops of mutual inhibition (between M BF and its inhibitor Rum1) and mutual phosphorylation (between skhan@chungbuk.ac.kr; Fax: MBF and Ste9); G2=M is composed by interlocking two PFLs, MP F with Cdc25 and MP F with W ee1. Other examples can be found in Ref. 9. Interlocking two or more PFLs together reinforces the properties of a single PFL and might bring new features that are not present in the single PFL case. In this study, we identied all possible dynamics of the interlocking of two types of PFL by performing a bifurcation analysis [11] in the one and the two parameters space of the stimulus amplitude and coupling strength. Finally, we constructed and compared phase diagrams in two coupling strengths. The diagrams show the dynamical structure of interlocking of the two types of PFLs. II. TWO DIFFERENT TYPES OF POSITIVE FEEDBACK LOOPS Here, we examined two types of positive feedback loop. One is the so-called self-enhancing PFL (epfl) in which the output P somehow increases its own synthesis rate; the other is the so-called self-recovering PFL (rpfl) in which the output P somehow removes its suppression. The two typical models of epfl and rpfl are auto-regulation (Figure 1(a)) and mutual phosphorylation (Figure 1(b)), respectively. er 1. Self-enhancing Positive Feedback: an Ampli- In the auto-regulation (Figure 1(a)), the output P positively regulates its own upstream regulator, T F, which,
2 Dynamic Analysis of Interlocked Positive Feedback Loops { Cuong Nguyen and Seung Kee Han Fig. 1. (a) Autoregulation: output P positively regulates its own upstream regulator, T F, which, in turn, increases the synthesis rate of output P. (b) Mutual phosphorylation: the output P (red) removes its inhibited form P p (grey) by activating an enzyme E. Here, [P ] and [T F ] denote the concentrations of the output P and the regulator T F, respectively. The output P synthesis is proportional to stimulus S and is regulated by a transcription factor T F. Its degradation is proportional to itself. The activation and the deactivation of the regulator T F follow the Michaelis-Menten kinetic. The k a [P ] and the k itf are the maximal rates of activation and deactivation, respectively. The J atf and the J itf are Michaelis constants. The activation of the T F is induced by P with a coupling strength constant k a, which is the eect of P on T F. The chosen parameters are k s1 = 0.4, k s2 = 1, k d = 2, k a = 1, k itf = 0.5 and J atf = J itf = The parameters are dimensionless. The small Michaelis constants, J's, are to ensure that the autoregulation possesses non-linearity, which is required to exhibit bistability [13]. Taking the stimulus S as the primary parameter, we created a one-parameter bifurcation diagram and steadystate values of the output P as a function of the stimulus S for dierent positive coupling strengths k a as in Figure 2(a)-(c). Fig. 2. Auto-regulation bifurcation diagram of a oneparameter bifurcation diagram: (a) monostability (k a = 0.01), (b) bistability (k a = 1), (c) irreversible bistability (k a = 1.5), (d) two-parameter bifurcation diagram of input S and coupling strength k a. Solid and dashed lines denote stable and unstable steady states, respectively. Dashed-dotted lines are a bifurcation diagram of auto-regulation without a feedback loop (k a = 0). in turn, increases the synthesis rate of the output P. Using the law of mass reaction and Michaelis-Menten kinetic [12], we can write the dierential equations governing the output P and the transcription factor T F of P as d[p ] = k s1 S + k s2 [T F ] k d [P ]; (1) d[t F ] = k a[p ](1 [T F ]) k itf [T F ] J atf + 1 [T F ] J itf + 1 [T F ] : (2) A. Monostability for weak coupling strength With a very small coupling strength k a, autoregulation exhibits monostability, as depicted by a monoincreasing line because there is no feedback eect in Figure 2(a). B. Bistability for strong coupling strength As the coupling strength k a gets stronger, the autoregulation exhibits bistability, which is the coexistence of two stable xed points (solid lines) separated by an unstable one (dashed line, Figure 2(b)). As long as the input S increases, the output P increases and reaches the threshold P 1 and then turns on the helper T F. The helper T F, in turn, increases the output level abruptly by increasing the output synthesis rate to a higher state (upper branch). If the input S decreases, the output P then decreases and falls below another threshold P 2 at which the output P is not strong enough to maintain the T F at the ON state. Therefore, the output P dramatically drops back to a lower state (lower branch). The two transitions occur at two dierent input S, called two saddle node bifurcation points, SN1 auto and SN2 auto, where the feedback loop turns on and o, respectively. To turn the feedback loop on, activating
3 Journal of the Korean Physical Society, Vol. 53, No. 6, December 2008 transcription factor T F, the output P has to exceed the threshold P 1. Similarly, P 2 is the threshold for turning the feedback o. The thresholds P 1 and P 2 are called feedback on/o levels. Also, the region of SN2 auto < S < SN1 auto is called the bistable region in which there are two stable and one unstable steady states. Here, the bistability of auto-regulation is said to be reversible because the output P can be ipped up or down by driving the stimulus S beyond SN1 auto or below SN2 auto. C. Irreversibility for very strong coupling strength As the coupling strength gets extremely large, the auto-regulation still exhibits bistability, but is irreversible (Figure 2(c)). This means that, once P ips up at S > SN1 auto, it never ips down, not even when all stimuli S are washed out (S = 0). The epfl could exhibit dierent dynamical behaviors depending on the coupling strength k a. As shown in the Figure 2(d), along the coupling axis, k a, two saddle node bifurcation points SN1 auto and SN2 auto emerge at k a = ka 1 (ka 1 = ); then move along to the left, the SN2 auto crosses zero at k a = ka 2 (ka 2 = ). Thus, there are three regions along the k a -axis: monostability k a < ka 1, bistability ka 1 < k a < ka 2 and irreversible bistability ka 2 < k a. We note that, once the helper T F turns on, it remarkably enhances the output P (solid line) in comparison with that of the simple regulation (dashed dotted line). This kind of positive feedback loop is, therefore, called a self-enhancing PFL (epfl). The epfl works as an amplier that amplies the output P. 2. Self-recovering Positive Feedback: the Buer In the PFL of mutual phosphorylation, the output P (red) removes the inhibited form P p (gray) by activating the enzyme E (Figure 1(b)). The active output P and its inactive form P p are thus antagonistic. The mathematical model of mutual phosphorylation is written in the following equations: d[p ] T = k s S k d [P ] T ; (3) d[p ] = k s S + k ap[e]([p ] T [P ]) J ap + [P ] T [P ] k ip [P ] k d [P ]; (4) J ip + [P ] d[e] = k b[p ](1 [E]) J ae + 1 [E] k ie [E] J ie + [E] : (5) Here, [P ] T, [P ] and [E] denote the total output P, the output P and the enzyme E, respectively. We note that output P can be in a free form P or in a phosphorylated Fig. 3. Mutual phosphorylation: one-parameter bifurcation diagram of P (solid and dashed lines): (a) monobistability, (b) bistability and (c) two-parameter bifurcation diagram of input S and coupling strength kb. Solid and dashed lines denote stable and unstable steady states, respectively. Dashdotted lines denote a one-parameter bifurcation diagram of total P. form P p ([P ] + [P ] p = [P ] T, where P T denotes the total P ). Here, the coupling strength k b is the aectivity of P on activation of the enzyme E. The parameters are k s = 0.4, k d = 2, k ap = 1, k ip = 0.5, k b = 1.3, k ie = 0.25 and J ap = J ip = J ae = J ie = The parameters are dimensionless. The two-parameter bifurcation diagram of mutual phosphorylation (Figure 3(c)) indicates that, for a coupling strength k b smaller than a certain value, k b < kb1 (kb1 = ), mutual phosphorylation exhibits monostability (Figure 3(a)). It exhibits bistability (Figure 3(b)) when k b gets stronger, kb1 < k b. It does not, however, exhibit irreversibility, as shown in the twoparameter bifurcation diagram. Instead of crossing the vertical axis, the SN2 m phos line approaches the vertical axis while the stimulus S reaches 0. It is remarkable to note that, with strong coupling strength, k b > kb1 (Figure 3(b)), the output (solid line) is initially suppressed from simple regulation (dashed dotted line) by an inactive form of P p until the input S is strong enough, S > SN1 m phos, to activate the enzyme E, which recovers the output P by removing the inactive form. Therefore, the inactive form P p works as a buer to \store" the output P until the buer is fullled and then broken by the output P. This kind of PFL is called a self-recovering PFL (rpfl). III. INTERLOCKED POSITIVE FEEDBACK LOOPS (IPFLS) In a biological complex network, two or more PFLs could share common components to work together [3{
4 Dynamic Analysis of Interlocked Positive Feedback Loops { Cuong Nguyen and Seung Kee Han Fig. 4. Interlocking of two positive feedback loops: autoregulation and mutual phosphorylation. The output P exists in two forms, active P (left) and inactive form (right, phosphorylated form). The active output P could activate its helper T F and enzyme E. In turn, the helper T F enhances the synthesis rate of P, while the enzyme E dephosphorylates the inactive form of the output P to become an active form. k a and k b are the coupling strengths of autoregulation and mutual phosphorylation, respectively. 7,10,13]. As shown in Figure 4, an epfl and an rpfl share a common output P composed of interlocked PFLs. In this ipfl, the output P is cooperatively regulated by auto-regulation and mutual phosphorylation. The stimulus S and the transcription factor T F regulate P together; and the output P can be dephosphorylated or phosphorylated with the help of enzyme E. Both T F and E are activated by P with coupling strengths k a and k b, respectively (Figure 4). The mathematical model and parameters are inherited from individual PFLs as in the following equations: d[p ] T = k s S k d [P ] T ; (6) d[p ] = k s S + k ap[e]([p ] T [P ]) k ip [P ] k d [P ]; J ap + [P ] T [P ] J ip + [P ] (7) d[t F ] d[e] = k a[p ](1 [T F ]) J atf + 1 [T F ] = k b[p ](1 [E]) J ae + 1 [E] k itf [T F ] J itf + [T F ] ; (8) k ie [E] J ie + [E] : (9) All the parameter values here are the same as those for single PFLs. 1. Bifurcation Analysis in Two-dimensional Parameter Space Let us start with the case in which the epfl coupling strength k a is xed in a bistable region, ka 1 < k a < ka 2 Fig. 5. Two-parameter bifurcation diagram of the interlocked PFLs (black) on (S, k b ) parameter space with xed k a (k a = 1), the auto-regulation (dashed blue lines) with coupling strength k a = 1 and mutual phosphorylation (dashed red lines). The value of the k's are k 1 = 0.063, k 2 = 0.221, k 3 = 0.224, k 4 = 0.315, k 5 = 0.393, k 6 = 0.582, k 7 = 1.031; k 8 = 2.8 and k 9 = (Figure 2(d)) while the rpfl coupling strength k b is varying. We performed a two-parameter bifurcation analysis of input S and rpfl coupling strength k b. The two-parameter bifurcation diagram of ipfl (black) is shown in Figure 5. Along the k b axis, two saddle node bifurcation points, SN1e and SN2e, which come from a single epfl, always exist. Meanwhile, other two saddle node bifurcation points, SN2m and SN1m, which correspond to those of a single rpfl, emerge at k b = k 1 and then move to the left and interfere with the SN1e and SN2e. The SN2m crosses the SN1e and SN2e at k b = k 2 and k b = k 4, respectively. The SN1m crosses the SN1e and then collapses with SN2e at k b = k 3 and k b = k 5, respectively. In the region of k 5 < k b < k 6, SN2m turns to be SN2e and SN1e turns to be SN1m due to the on/o level swapping between the epfl and the rpfl (refer to the next section for more details). The creation of SN2m and SN1e occurs at k b = k 6. Then, the SN1m crosses SN1e and SN2e at k b = k 7 and k b = k 9, respectively. The SN2m crosses SN2e at k b = k 8. Based on the interference between those saddle node bifurcation points, the k b axis could be divided into ten regions of dierent dynamical behaviors of ipfls as specied by the k's on the right axis. The detailed dynamical behavior of each region is explained in the next section. 2. Bifurcation Analysis in One-dimensional Parameter Space The ten regions along the k b axis are represented by one-parameter bifurcation diagrams, as shown in Figure 6(a)-(j), respectively. In Figure 6(a)-(j), the epfl coupling strength k a is xed so that the single epfl could
5 Journal of the Korean Physical Society, Vol. 53, No. 6, December 2008 Fig. 6. One-parameter bifurcation diagrams of epfl (blue), rpfl (red) and ipfl (black) with k a xed (k a = 1) for various values of k b : 0.05, 0.2, 0.23, 0.3, 0.35, 0.5, 0.8, 2, 5, 26 in correspondence to (a), (b), (j). yield bistability; meanwhile. the rpfl coupling strength k b is chosen according to the regions specied by k's in Figure 5. SN2e < S < SN1e and SN2m < S < SN1m, are separated from each other. Consequently, this conguration of the ipfl bifurcation diagram is called non-overlapping bistability (NoB). A. k b < k 1 In this region, a single rpfl exhibits suppressed monostability (red line in Figure 6(a)). Meanwhile, a single epfl could yield bistability with bistable regions of SN2 auto < S < SN1 auto (blue line in Figure 6(a)). Therefore, the ipfls exhibits a single bistable region of SN1e < S < SN2e (Figure 6(a)). However, the bistable region of ipfl is at larger stimulus S than that of epfl, SN2 auto < S < SN1 auto, due to suppression in rpfl. B. k 1 < k b < k 2 In the region of k 1 < k b < k 2, the single rpfl exhibits bistability, SN2m phos < S < SN 1m phos, at very large stimulus S (Figure 6(b)). Therefore, ipfl exhibits an additional bistable region, SN 2m < S < SN 1m, coming from rpfl, beside the bistable region coming from epfl, SN2e < S < SN1e. However, the bistable region of ipfl, SN 1m < S < SN 2m is at smaller stimulus S than that of rpfl, SN2m phos < S < SN1m phos due to the enhancement of epfl, which is activated rst. Here, the two bistable regions of ipfl, C. k 2 < k b < k 3 By increasing k b, the two bifurcation points SN1m and SN2m move to the left and down and then interfere with SN1e and SN2e. In the region of k 2 < k b < k 3, the two bistable regions overlap, SN 2e < SN 2m < SN1e < SN1m (Figure 6(c)), creating a tristability region, SN2m < S < SN1e. Tristability means that three stable xed points coexist and are separated by two unstable xed points. Therefore, this conguration is called the overlapping bistable region (OLB). D. k 3 < k b < k 4 In this region, the bistable region of SN2m < S < SN1m is located inside the bistable region of SN2e < S < SN 1e, making a nesting bistability (NB) (Figure 6(d)). In this NB, switching from the lowest stable state to the highest one requires only one transition. Meanwhile, switching from the highest state to the lowest one requires two transitions. E. k 2 < k b < k 3
6 Dynamic Analysis of Interlocked Positive Feedback Loops { Cuong Nguyen and Seung Kee Han Here, two bistable regions merge together to create a large bistable region with a single transition between the two lowest and the highest stable states, SN2m < S < SN 1e (Figure 6(e)). Though inaccessible, there is a small region of tristability. Therefore, this is called merge bistability (MB'). F. k 5 < k b < k 6 In this region, two bistable regions emerge and the SN 1m and SN 2e collide and disappear, creating a bistable region without tristability inside (Figure 6(f)), the so-called merged bistability (MB). The merged bistable region is much larger than the two single bistable regions. This would enhance the reliability of biological toggle switches. G. k 6 < k b < k 7 Here, the ipfl exhibits another merged bistability with tristability (MB') again (Figure 6(g)). However, the rpfl on/o level is less than that of epfl in this conguration. Therefore, the bifurcation points of ipfl, SN1e, SN2e, SN1m and SN2m, occur at those of a single epfl and rpfl, respectively. H. k 7 < k b < k 8 Other nesting bistability (NB) is shown (Figure 6(h)). I. k 8 < k b < k 9 Other overlapping bistability (OLB) is shown (Figure 6(i)). J. k 8 < k b < k 9 Other non-overlapping bistability (NoB) is shown (Figure 6(j)). Thus, since k b > k 1, a single rpfl exhibits bistability and its two loci of saddle node bifurcations move to the left when k b is increasing (Figure 3(c)). As depicted in Figure 6(b)-(j), when k b is increased, the bistable region of a single rpfl, SN 2m-phos<S<SN 1mphos (red), moves to the left and down while those of a single epfl, SN1e and SN2e (blue), are xed. As a consequence, the ipfls exhibit an additional bistable region, SN2m < S < SN1m to SN2e < S < SN1e. Also, SN 2m < S < SN 1m moves to the left and down and then interferes with former bistable regions (SN2e < S < SN1e), creating diverse types of dynamical behaviors of ipfls. 3. Diversity of Dynamical Behaviors of ipfls To sum up, depending on the two coupling strengths k a and k b, interlocked system shows diverse dynamical behaviors. 1. Single bistability, SB, in which usually comes out when either auto-regulation or mutual activation exhibits bistability while the other does not (Figure 6(a)). The output turns on when the input exceeds the bifurcation point SN1e and the output turns o when the input S decreases below SN2e. 2. Non-overlapping bistability, NoB, in which bistable regions of each positive feedback loop are separated from each other (Figure 6(b) and (j)). To completely turn on the output P, the input S must reach SN1e to make the rst transition and then reach SN1m to make the second transition. As to turning o the output P, there are two transitions at SN2m and SN2e. 3. Overlapping bistability, OLB, in which two stable regions of an interlocked system overlap each other, creating tristability. Tristability means that three stable xed points coexist and are separated by two unstable xed points (Figure 6(c) and (i)). Turning on or o the output requires two transitions of output P at bifurcations. 4. Nesting bistability, NB, in which one bistable region nests in the other (Figure 6(d) and (h)). In this type, the transition number is asymmetric. For example, turning on the output P requires only one transition while two transitions are required to turn o the output P (Figure 6(d)). 5. Merging bistability, MB and MB', in which two individual bistable regions become a large merged one (Figure 6(e)-(g)). In this region, two single PFLs, epfl and rpfl, cooperate with each other to create a large bistable region. We also searched for all possible dynamical behaviors of an ipfl for a rpfl coupling strength k b when the epfl coupling strength k a is in the irreversible bistable region, k a > ka 2 (Figure 2(d)). All possible one-parameter bifurcation diagrams of a single epfl (blue), a single rpfl (red) and an ipfl (black) are shown in Figure 7(a)-(j). Because the epfl bistable region (blue) is irreversible, the ipfl bistable region (black) is irreversible. Despite that, the dynamical behavior of the ipfl is similar to that described in the case of ka 1 < k a < ka 2 (Figure 6(a)-(j)), except that their bistability is irreversible. Therefore, these dynamical behaviors could be classied into dierent groups as follows: 1. Irreversible single bistability (isb), Figure 7(a). There is one irreversible bistable region.
7 Journal of the Korean Physical Society, Vol. 53, No. 6, December 2008 Fig. 7. One-parameter bifurcation diagrams of epfl with strong coupling strength (blue), rpfl (red) and ipfl (black). The epfl coupling strength is k a = 3.5. The rpfl coupling strength is k b = [0.05, 0.29, 0.34, 0.5, 0.58, 0.7, 0.9, 2, 10.0, 30.0] in correspondence to (a), (b), (j). The ipfl coupling strengths k a and k b are same as those of a single epfl and rpfl. one of them is irreversible. 3. Irreversible overlapping bistability (iolb), Figure 7(c). Two bistable regions overlap creating tristability; one of them is irreversible. 4. Irreversible nesting bistability (inb). One bistable region nests in the other; one of them is irreversible (Figure 7(d) and (j)). 5. Double irreversible nesting bistability (i 2 NB). One bistable region nests in the other; both of them are irreversible (Figure 7(e)). Fig. 8. Phase diagram on the coupling-strength plane of (k a, k b ). Mono: monostability; SB: single bistability; NoB: non-overlapped bistability; OLB: overlapped bistability; NB: nested bistability; and MB: merged bistability. \i" and \i 2 " denote irreversibility and double irreversibility. Regions with the same color are of the same type. See the text for more details. 2. Irreversible non-overlapping bistability (inob), Figure 7(b). Two bistable regions are separated; 6. Irreversible merged bistability with or without tristability inside (imb and imb'), Figure 7(f)- (i). Two individual bistable regions become a large merged one, which is irreversible. and k b 4. Phase Diagram of the coupling strength k a The ipfls are composed of two PFLs with two different positive coupling strengths, k a and k b. To understand more how ipfls depend on their two coupling strengths, we built up a parameter-phase diagram in the (k a, k b ) plane (Figure 8). The whole phase plane
8 Dynamic Analysis of Interlocked Positive Feedback Loops { Cuong Nguyen and Seung Kee Han of (k a, k b ) is divided into seven regions with dierent colors, as shown in Figure 8. The monostability region (Mono) at small k a and k b (k a < ka 1, k b < k 1 ) is where both epfl and rpfl exhibit monostability; so the ipfls also exhibit monostability. The single bistability regions with or without irreversibility (SB and isb) are where either epfl or rpfl exhibits bistability, (ka 1 < k a and k b < k 1 ) or (k a < ka 1 and k 1 < k b ). The region of (ka 1 < k a, k 1 < k b ) is divided into four small regions and show a diversity of dynamical behaviors of ipfl, such as non-overlapping bistability with or without irreversibility (NoB and inob), overlapping bistability with or without irreversibility (OLB and iolb), nesting bistability with or without irreversibility (NB, inb and i 2 NB) and merged bistability with or without irreversibility (MB, MB', imb and imb'). IV. CONCLUSION In complex biological systems, epfl and rpfl cooperate and interlock together rather than working alone [9]. Our dynamical analysis indicated that, by interlocking two PFLs, the ipfl would show a diversity of dynamical behaviors, such as single bistability (SB), non-overlapping bistability (NoB), overlapping bistability (OLB), nesting bistability (NB) and merging bistability (MB) (Figure 6), or irreversible single bistability (isb), irreversible non-overlapping bistability (inob), double irreversible non-overlapping bistability (i 2 NoB), irreversible overlapping bistability (iolb), irreversible nesting bistability (inb), double irreversible Nesting bistability (i 2 NB) and irreversible Merging bistability (imb) (Figure 7). Those dynamical behaviors could be easily transformed to any one of the others by changing one of the coupling strengths, k a or k b. This provides exibility for biological systems in response to various external stimuli. Remarkable enough, the bistable region of ipfl in the MB type is greatly enhanced in comparison to that of a single one (Figure 6(e)-(g) and Figure 7(f)-(i)). In this region, the epfl and the rpfl cooperate to greatly improve the reliability of the decisive switch. The cooperation of the two PFLs creates a more reliable toggle switch, which is important in biological decisive systems. In this work, positive feedback loops were classied into two groups, self-enhancing (epfl) and self-recovery (rpfl) positive feedback loops. Then the epfl and rpfl were interlocked together, the interlocked positive feedback loop (ipfl) showed diverse dynamical behaviors with respect to the relative coupling strengths of the feedback loops. Interestingly, when the epfl and the rpfl cooperate with comparable coupling strengths, the two bistable regions of the single epfl and rpfl merge together, resulting in a much larger bistable region. Therefore, cooperativity of the epfl and the rpfl greatly enhances the reliability of switches. ACKNOWLEDGMENTS This work was supported by a Korea Science and Engineering Foundation grant (R ) funded by the Ministry of Science of Technology, Korea and also by the BK21 program at Chungbuk National University. REFERENCES [1] L. H. Hartwell, J. J. Hopeld, S. Leibler and A. W. Murray, Nature 402, C47 (1999). [2] C. Nguyen, C. N. Yoon, S. K. Han, D. L. Pham and H. Y. Kim, J. Korean Phys. Soc. 50, 332 (2007). [3] J. J. Tyson, K. C. Chen and B. Novak, Curr. Opin. Cell. Biol. 15, 221 (2003). [4] K. C. Chen, A. Csikasz-Nagy, B. Gyory, J. Val, B. Novak and J. J. Tyson, Mol. Biol. Cell. 11, 369 (2000); C. N. Yoon, S. K. Han and H. Y. Kim, J. Korean Phys. Soc. 44, 638 (2004). [5] B. Novak, Z. Pataki, A. Ciliberto and J. J. Tyson, Chaos 11, 277 (2001). [6] Z. Qu, W. R. MacLellan and J. N. Weiss, Biophys. J. 85, 3600 (2003). [7] P. Smolen, D. A. Baxter and J. H. Byrne, J. Neurosci. 21, 6644 (2001). [8] M. Swat, A. Kel and H. Herzel, Bioinformatics 20, 1506 (2004). [9] O. Brandman, J. E. Ferrell, Jr., R. Li and T. Meyer, Science 310, 496 (2005). [10] T. I. Lee, N. J. Rinaldi, F. Robert, D. T. Odom, Z. Bar- Joseph, G. K. Gerber, N. M. Hannett, C. T. Harbison, C. M. Thompson, I. Simon, J. Zeitlinger, E. G. Jennings, H. L. Murray, D. B. Gordon, B. Ren, J. J. Wyrick, J. B. Tagne, T. L. Volkert, E. Fraenkel, D. K. Giord and R. A. Young, Science 298, 799 (2002). [11] S. G. Lee, J. Korean Phys. Soc. 52, 11 (2008). [12] G. E. Briggs and J. B. Haldane, Biochem J. 19, 338 (1925). [13] J. E. Ferrell and W. Xiong, Chaos 11, 227 (2001).
Cell cycle regulation in the budding yeast
Cell cycle regulation in the budding yeast Bởi: TS. Nguyen Cuong Introduction The cell cycle is the sequence of events by which a growing cell duplicates all its components and then divides into two daughter
More informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 Modelling of Signalling Pathways Curr. Op. Cell Biol. 15 (2003) 221 1) How do the magnitudes of signal output and signal duration depend on the kinetic
More informationBioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"
Bioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014" Modelling of Signalling Pathways" Curr. Op. Cell Biol. 15 (2003) 221" 1) How do the magnitudes of signal output and signal duration depend on the
More informationnutrients growth & division repellants movement
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Outline 1. Cell Signaling: Physiology 2. Cell Signaling: Molecular Biology 3. Chemical
More informationNetwork Dynamics and Cell Physiology. John J. Tyson Dept. Biological Sciences Virginia Tech
Network Dynamics and Cell hysiology John J. Tyson Dept. Biological Sciences Virginia Tech Collaborators Budapest Univ. Techn.. & Econ. Bela Novak Attila Csikasz-Nagy Andrea Ciliberto Virginia Tech Kathy
More informationAnalysis and Simulation of Biological Systems
Analysis and Simulation of Biological Systems Dr. Carlo Cosentino School of Computer and Biomedical Engineering Department of Experimental and Clinical Medicine Università degli Studi Magna Graecia Catanzaro,
More informationA simple model for the eukaryotic cell cycle. Andrea Ciliberto
A simple model for the eukaryotic cell cycle Andrea Ciliberto The cell division cycle G1 cell division Start S (DNA Replication) Finish M (mitosis) G2/M G2 Kohn, Mol. Biol. Cell., 1999 How did we get to
More informationNetwork Dynamics and Cell Physiology. John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Signal Transduction Network Hanahan & Weinberg (2000) Gene Expression Signal-Response
More informationA minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle
A minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle K. Sriram, G. Bernot and F. Képès Abstract: A novel topology of regulatory networks abstracted from the
More informationProblem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)
Problem Set 2 1 Competitive and uncompetitive inhibition (12 points) a. Reversible enzyme inhibitors can bind enzymes reversibly, and slowing down or halting enzymatic reactions. If an inhibitor occupies
More informationMultistability in the lactose utilization network of Escherichia coli
Multistability in the lactose utilization network of Escherichia coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Agenda Motivation Intro to multistability Purpose
More informationDNA replication. M (mitosis)
Introduction to cell cycle modeling Attila Csikász-Nagy G1 The cell cycle is the sequence of events whereby a growing cell replicates all its components and divides them more-or-less evenly between two
More informationRegulation of the mammalian cell cycle: a model of the G 1 -to-s transition
Am J Physiol Cell Physiol 284: C349 C364, 2003. First published October 9, 2002; 10.1152/ajpcell.00066.2002. Regulation of the mammalian cell cycle: a model of the G 1 -to-s transition ZHILIN QU, JAMES
More informationANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides
ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY Nael H El-Farra, Adiwinata Gani & Panagiotis D Christofides Department of Chemical Engineering University of California, Los Angeles 2003 AIChE
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationModelling the cell cycle regulatory network
Chapter 3 Modelling the cell cycle regulatory network 3.1 Dynamical modelling Dynamical properties of a system arise from the interaction of its components. In the case of the cell division cycle, cellular
More informationarxiv: v1 [q-bio.mn] 7 Nov 2018
Role of self-loop in cell-cycle network of budding yeast Shu-ichi Kinoshita a, Hiroaki S. Yamada b a Department of Mathematical Engineering, Faculty of Engeneering, Musashino University, -- Ariake Koutou-ku,
More informationStochastic Models for Circadian Oscillations: Emergence of a Biological Rhythm
Stochastic Models for Circadian Oscillations: Emergence of a Biological Rhythm DIDIER GONZE, JOSÉ HALLOY, ALBERT GOLDBETER Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles,
More informationMultistability in the lactose utilization network of E. coli. Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J.
Multistability in the lactose utilization network of E. coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Motivation Understanding biological switches in the
More informationBioControl - Week 2, Lecture 1
BioControl - Week 2, Lecture Goals of this lecture: CDS tools for bio-molecular systems Study of equilibria Dynamic performance Periodic behaviors Main eample: Negative autoregulation Suggested readings
More informationCore genetic module: The mixed feedback loop
PHYSICAL REVIEW E 72, 031908 2005 Core genetic module: The mixed feedback loop Paul François and Vincent Hakim Laboratoire de Physique Statistique, * CNRS-UMR 8550, Ecole Normale Supérieure, 24, rue Lhomond
More informationDr. Fred Cross, Rockefeller (KITP Bio Networks 3/26/2003) 1
Outline Cell growth as the driver for cell cycle (in microbes): coordination of growth and division A basic principle organizing cell cycle control: why cyclin-dependent kinase activity must oscillate
More informationSuperstability of the yeast cell-cycle dynamics: Ensuring causality in the presence of biochemical stochasticity
Journal of Theoretical Biology 245 (27) 638 643 www.elsevier.com/locate/yjtbi Superstability of the yeast cell-cycle dynamics: Ensuring causality in the presence of biochemical stochasticity Stefan Braunewell,
More informationPlant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter
Plant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter 9/10/2008 1 Learning Objectives Explain why a cell cycle was selected for during evolution
More information56:198:582 Biological Networks Lecture 10
56:198:582 Biological Networks Lecture 10 Temporal Programs and the Global Structure The single-input module (SIM) network motif The network motifs we have studied so far all had a defined number of nodes.
More informationNotes, Mathematical Cell Biology Course. Leah Edelstein-Keshet
Notes, Mathematical Cell Biology Course Leah Edelstein-Keshet May 2, 2012 ii Leah Edelstein-Keshet Contents 1 Qualitative behaviour of simple ODEs and bifurcations 1 1.0.1 Cubic kinetics.......................
More informationAlternating Oscillations and Chaos in a Model of Two Coupled Biochemical Oscillators Driving Successive Phases of the Cell Cycle
Alternating Oscillations and Chaos in a Model of Two Coupled Biochemical Oscillators Driving Successive Phases of the Cell Cycle PIERRE-CHARLES ROMOND, a MAURO RUSTICI, b DIDIER GONZE, AND ALBERT GOLDBETER
More informationResearch Article Robust Cell Size Checkpoint from Spatiotemporal Positive Feedback Loop in Fission Yeast
BioMed Research International Volume 213, Article ID 91941, 9 pages http://dx.doi.org/1.1155/213/91941 Research Article Robust Cell Size Checkpoint from Spatiotemporal Positive Feedback Loop in Fission
More informationBiophysical Journal Volume 92 May
Biophysical Journal Volume 92 May 2007 3407 3424 3407 Dynamics of a Minimal Model of Interlocked Positive and Negative Feedback Loops of Transcriptional Regulation by camp-response Element Binding Proteins
More informationIntroduction to Biological Modeling
Introduction to Biological odeling Lecture : odeling dynamics Sept. 9, 010 Steve Andrews Brent lab, Basic Sciences Division, FHCRC Last week Why model biology? Example: E. coli chemotaxis Typical modeling
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process https://www.youtube.com/
More informationMATHEMATICAL MODELING OF MITOSIS
BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS FACULTY OF CHEMICAL TECHNOLOGY AND BIOTECHNOLOGY GYÖRGY OLÁH DOCTORAL SCHOOL MATHEMATICAL MODELING OF MITOSIS Summary of the thesis AUTHOR: Attila Tóth SUPERVISOR:
More informationMeasuring TF-DNA interactions
Measuring TF-DNA interactions How is Biological Complexity Achieved? Mediated by Transcription Factors (TFs) 2 Regulation of Gene Expression by Transcription Factors TF trans-acting factors TF TF TF TF
More informationAnalysis of Mode Transitions in Biological Networks
Analysis of Mode Transitions in Biological Networks Nael H. El-Farra, Adiwinata Gani, and Panagiotis D. Christofides Dept. of Chemical Engineering, University of California, Los Angeles, CA 90095 DOI 10.1002/aic.10499
More informationIntroduction to Mathematical Physiology I - Biochemical Reactions
Introduction to Mathematical Physiology I - Biochemical Reactions J. P. Keener Mathematics Department Math Physiology p.1/28 Introduction The Dilemma of Modern Biology The amount of data being collected
More information56:198:582 Biological Networks Lecture 11
56:198:582 Biological Networks Lecture 11 Network Motifs in Signal Transduction Networks Signal transduction networks Signal transduction networks are composed of interactions between signaling proteins.
More informationMathematical model of the cell division cycle of fission yeast
CHAOS VOLUME 11, NUMBER 1 MARCH 2001 Mathematical model of the cell division cycle of fission yeast Bela Novak a) and Zsuzsa Pataki Department of Agricultural Chemical Technology, Budapest University of
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More informationNetwork Dynamics and Cell Physiology
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Outline 1. The cell is an information processing system. 2. Information is processed
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationSystem-level feedbacks control cell cycle progression
FEBS Letters 583 (2009) 3992 3998 journal homepage: www.febsletters.org Review System-level feedbacks control cell cycle progression Orsolya Kapuy a, Enuo He a, Sandra López-Avilés b, Frank Uhlmann b,
More informationModelling the fission yeast cell cycle Akos Sveiczer, John J. Tyson and Bela Novak Date received (in revised form): 20th October, 2003
Akos Sveiczer is a Reader in the Department of Agricultural Chemical Technology at Budapest University of Technology and Economics. John J. Tyson is a University Distinguished Professor in the Department
More informationStochastic simulations
Stochastic simulations Application to molecular networks Literature overview Noise in genetic networks Origins How to measure and distinguish between the two types of noise (intrinsic vs extrinsic)? What
More informationProceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007
Proceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007 A FEEDBACK APPROACH TO BIFURCATION ANALYSIS IN BIOCHEMICAL NETWORKS WITH MANY PARAMETERS Steffen Waldherr 1 and Frank Allgöwer Institute
More informationUNIVERSITY OF CALIFORNIA SANTA BARBARA DEPARTMENT OF CHEMICAL ENGINEERING. CHE 154: Engineering Approaches to Systems Biology Spring Quarter 2004
UNIVERSITY OF CALIFORNIA SANTA BARBARA DEPARTMENT OF CHEMICAL ENGINEERING CHE 154: Engineering Approaches to Systems Biology Spring Quarter 2004 Lecture: Tue/Thu 9:30-10:45am Engr-II Room 3301 Instructor
More informationIs chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?
Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value
More informationIntroduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.
Introduction Dagmar Iber Jörg Stelling joerg.stelling@bsse.ethz.ch CSB Deterministic, SS 2015, 1 Origins of Systems Biology On this assumption of the passage of blood, made as a basis for argument, and
More informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
More informationFrom quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle
ORIGINAL RESEARCH ARTICLE published: 02 November 2012 doi: 10.3389/fphys.2012.00413 From quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle Claude Gérard and Albert Goldbeter*
More informationBistability in cell signaling: How to make continuous processes discontinuous, and reversible processes irreversible
CHAOS VOLUME 11, NUMBER 1 MARCH 2001 Bistability in cell signaling: How to make continuous processes discontinuous, and reversible processes irreversible James E. Ferrell, Jr. a) and Wen Xiong Department
More informationStem Cell Reprogramming
Stem Cell Reprogramming Colin McDonnell 1 Introduction Since the demonstration of adult cell reprogramming by Yamanaka in 2006, there has been immense research interest in modelling and understanding the
More informationBifurcation and Singularity Analysis of a Molecular Network for the Induction of Long-Term Memory
Biophysical Journal Volume 90 April 2006 2309 2325 2309 Bifurcation and Singularity Analysis of a Molecular Network for the Induction of Long-Term Memory Hao Song, Paul Smolen, Evyatar Av-Ron, Douglas
More informationMonomeric Bistability and the Role of Autoloops in Gene Regulation
Monomeric Bistability and the Role of Autoloops in Gene Regulation Stefanie Widder Javier Macía Ricard V. Solé SFI WORKING PAPER: 2008-10-044 SFI Working Papers contain accounts of scientific work of the
More informationSimulation of Gene Regulatory Networks
Simulation of Gene Regulatory Networks Overview I have been assisting Professor Jacques Cohen at Brandeis University to explore and compare the the many available representations and interpretations of
More informationCS-E5880 Modeling biological networks Gene regulatory networks
CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene
More informationModeling Multiple Steady States in Genetic Regulatory Networks. Khang Tran. problem.
Modeling Multiple Steady States in Genetic Regulatory Networks Khang Tran From networks of simple regulatory elements, scientists have shown some simple circuits such as the circadian oscillator 1 or the
More informationRobust Network Topologies for Generating Switch-Like Cellular Responses
University of Pennsylvania ScholarlyCommons Departmental Papers (BE) Department of Bioengineering 6-22-2011 Robust Network Topologies for Generating Switch-Like Cellular Responses Najaf A. Shah University
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process DNA codes for DNA replication
More informationSystems biology and complexity research
Systems biology and complexity research Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Interdisciplinary Challenges for
More informationDynamic modeling and analysis of cancer cellular network motifs
SUPPLEMENTARY MATERIAL 1: Dynamic modeling and analysis of cancer cellular network motifs Mathieu Cloutier 1 and Edwin Wang 1,2* 1. Computational Chemistry and Bioinformatics Group, Biotechnology Research
More information56:198:582 Biological Networks Lecture 9
56:198:582 Biological Networks Lecture 9 The Feed-Forward Loop Network Motif Subgraphs in random networks We have discussed the simplest network motif, self-regulation, a pattern with one node We now consider
More informationBistability from double phosphorylation in signal transduction
Bistability from double phosphorylation in signal transduction Kinetic and structural requirements Fernando Ortega 1,2,3, José L. Garcés 1,2, Francesc Mas 1,2, Boris N. Kholodenko 4 and Marta Cascante
More informationAn Introduction to the spls Package, Version 1.0
An Introduction to the spls Package, Version 1.0 Dongjun Chung 1, Hyonho Chun 1 and Sündüz Keleş 1,2 1 Department of Statistics, University of Wisconsin Madison, WI 53706. 2 Department of Biostatistics
More informationLecture 8: Temporal programs and the global structure of transcription networks. Chap 5 of Alon. 5.1 Introduction
Lecture 8: Temporal programs and the global structure of transcription networks Chap 5 of Alon 5. Introduction We will see in this chapter that sensory transcription networks are largely made of just four
More informationDynamics of the Cell Cycle: Checkpoints, Sizers, and Timers
3600 Biophysical Journal Volume 85 December 2003 3600 3611 Dynamics of the Cell Cycle: Checkpoints, Sizers, and Timers Zhilin Qu, W. Robb MacLellan, and James N. Weiss Cardiovascular Research Laboratory,
More informationA synthetic oscillatory network of transcriptional regulators
A synthetic oscillatory network of transcriptional regulators Michael B. Elowitz & Stanislas Leibler, Nature, 403, 2000 igem Team Heidelberg 2008 Journal Club Andreas Kühne Introduction Networks of interacting
More informationA Model of Periodic Oscillation for Genetic Regulatory Systems
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL. 49, NO. 10, OCTOBER 2002 1429 A Model of Periodic Oscillation for Genetic Regulatory Systems Luonan Chen, Senior Member,
More informationFUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models
FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models Markus W. Covert Stanford University 0 CRC Press Taylor & Francis Group Boca Raton London New York Contents /... Preface, xi
More informationOverview of the cell cycle
Chapter 2 Overview of the cell cycle 2.1 The organisation of cell cycle in eukaryotes During the cell cycle, the typical eukaryotic cell goes through a series of well defined phases, to divide into two
More informationMonotonicity and Bistability of Calcium/Calmodulin-Dependent Protein Kinase-Phosphatase Activation
2010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 02, 2010 FrB13.2 Monotonicity and Bistability of Calcium/Calmodulin-Dependent Protein Kinase-Phosphatase Activation
More informationNo oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism
No oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism Luna Bozeman 1 and Adriana Morales 2 1 Department of Mathematical Sciences Clemson
More informationCOMPUTER SIMULATION OF DIFFERENTIAL KINETICS OF MAPK ACTIVATION UPON EGF RECEPTOR OVEREXPRESSION
COMPUTER SIMULATION OF DIFFERENTIAL KINETICS OF MAPK ACTIVATION UPON EGF RECEPTOR OVEREXPRESSION I. Aksan 1, M. Sen 2, M. K. Araz 3, and M. L. Kurnaz 3 1 School of Biological Sciences, University of Manchester,
More informationDynamics of the Mixed Feedback Loop Integrated with MicroRNA
The Second International Symposium on Optimization and Systems Biology (OSB 08) Lijiang, China, October 31 November 3, 2008 Copyright 2008 ORSC & APORC, pp. 174 181 Dynamics of the Mixed Feedback Loop
More informationDistinctive Topologies of Partner-switching Signaling Networks Correlate with their Physiological Roles
doi:10.1016/j.jmb.2007.04.021 J. Mol. Biol. (2007) 369, 1333 1352 Distinctive Topologies of Partner-switching Signaling Networks Correlate with their Physiological Roles Oleg A. Igoshin 1, Margaret S.
More informationChapter 15 Active Reading Guide Regulation of Gene Expression
Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,
More informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
More informationBi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016
Bi 8 Lecture 11 Quantitative aspects of transcription factor binding and gene regulatory circuit design Ellen Rothenberg 9 February 2016 Major take-home messages from λ phage system that apply to many
More information1 Periodic stimulations of the Incoherent Feedforward Loop network
1 Periodic stimulations of the Incoherent Feedforward Loop network In this Additional file, we give more details about the mathematical analysis of the periodic activation of the IFFL network by a train
More informationBiophysical Chemistry 72 (1998)
Biophysical Chemistry 72 (1998) 153 167 Link between fertilization-induced Ca 2+ oscillations and relief from metaphase II arrest in mammalian eggs: a model based on calmodulin-dependent kinase II activation
More informationModelling Biochemical Reaction Networks. Lecture 17: Modeling the cell cycle, Part I
Modelling Biochemical Reaction Networks Lecture 17: Modeling the cell cycle, Part I Marc R. Roussel Department of Chemistry and Biochemistry Recommended reading Fall, Marland, Wagner and Tyson, sections
More informationLecture 10: Cyclins, cyclin kinases and cell division
Chem*3560 Lecture 10: Cyclins, cyclin kinases and cell division The eukaryotic cell cycle Actively growing mammalian cells divide roughly every 24 hours, and follow a precise sequence of events know as
More informationA local sensitivity analysis of Ca 2+ -calmodulin binding and its influence over PP1 activity
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 A local sensitivity analysis of Ca 2+ -calmodulin binding and its influence
More information56:198:582 Biological Networks Lecture 8
56:198:582 Biological Networks Lecture 8 Course organization Two complementary approaches to modeling and understanding biological networks Constraint-based modeling (Palsson) System-wide Metabolism Steady-state
More informationCHAPTER 8. An Introduction to Metabolism
CHAPTER 8 An Introduction to Metabolism WHAT YOU NEED TO KNOW: Examples of endergonic and exergonic reactions. The key role of ATP in energy coupling. That enzymes work by lowering the energy of activation.
More informationButterfly color patterns are built as finely tiles mosaics of colored scales, each the product
CHAPTER 8 Stochastic Gene Expression: Dominance, Thresholds and Boundaries H. Frederik Nijhout Abstract Butterfly color patterns are built as finely tiles mosaics of colored scales, each the product of
More informationThe dynamics of cell cycle regulation
The dynamics of cell cycle regulation John J. Tyson, 1 * Attila Csikasz-Nagy, 2,3 and Bela Novak 2 Summary Major events of the cell cycle DNA synthesis, mitosis and cell division are regulated by a complex
More informationProblem Set # 3
20.320 Problem Set # 3 October 1 st, 2010 Due on October 8 th, 2010 at 11:59am. No extensions, no electronic submissions. General Instructions: 1. You are expected to state all your assumptions and provide
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationProgram for the rest of the course
Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling
More informationRotational Number Approach to a Damped Pendulum under Parametric Forcing
Journal of the Korean Physical Society, Vol. 44, No. 3, March 2004, pp. 518 522 Rotational Number Approach to a Damped Pendulum under Parametric Forcing Eun-Ah Kim and K.-C. Lee Department of Physics,
More informationCompact Modeling of Allosteric Multisite Proteins: Application to a Cell Size Checkpoint
Compact Modeling of Allosteric Multisite Proteins: Application to a Cell Size Checkpoint Germán Enciso 1 *, Douglas R. Kellogg 2, Arturo Vargas 3 1 Department of Mathematics, Department of Developmental
More informationStochastic Oscillator Death in Globally Coupled Neural Systems
Journal of the Korean Physical Society, Vol. 52, No. 6, June 2008, pp. 19131917 Stochastic Oscillator Death in Globally Coupled Neural Systems Woochang Lim and Sang-Yoon Kim y Department of Physics, Kangwon
More informationCooperative Effects of Noise and Coupling on Stochastic Dynamics of a Membrane-Bulk Coupling Model
Commun. Theor. Phys. (Beijing, China) 51 (2009) pp. 455 459 c Chinese Physical Society and IOP Publishing Ltd Vol. 51, No. 3, March 15, 2009 Cooperative Effects of Noise and Coupling on Stochastic Dynamics
More informationSrc Family Kinases and Receptors: Analysis of Three Activation Mechanisms by Dynamic Systems Modeling
Biophysical Journal Volume 94 March 2008 1995 2006 1995 Src Family Kinases and Receptors: Analysis of Three Activation Mechanisms by Dynamic Systems Modeling Hendrik Fuß, Werner Dubitzky, C. Stephen Downes,
More information2 Dilution of Proteins Due to Cell Growth
Problem Set 1 1 Transcription and Translation Consider the following set of reactions describing the process of maing a protein out of a gene: G β g G + M M α m M β m M + X X + S 1+ 1 X S 2+ X S X S 2
More informationSupplementary Figure 1. Confirmation of REF52-hE2F1p::4NLS-d4Venus reporter cells and characterization of E2F dynamics. (a) Alignment of E2F dynamics
Supplementary Figure 1. Confirmation of REF52-hE2F1p::4NLS-d4Venus reporter cells and characterization of E2F dynamics. (a) Alignment of E2F dynamics trajectories to endogenous E2F1 mrna expression. Gray
More informationSeminar in Bioinformatics, Winter Network Motifs. An Introduction to Systems Biology Uri Alon Chapters 5-6. Meirav Zehavi
Seminar in Bioinformatics, Winter 2010 Network Motifs An Introduction to Systems Biology Uri Alon Chapters 5-6 Meirav Zehavi Contents Network Motifs in - Sensory Transcription Networks - SIM - Multi-Output
More informationIn: Proc. BENELEARN-98, 8th Belgian-Dutch Conference on Machine Learning, pp 9-46, 998 Linear Quadratic Regulation using Reinforcement Learning Stephan ten Hagen? and Ben Krose Department of Mathematics,
More informationOscillatory Turing Patterns in a Simple Reaction-Diffusion System
Journal of the Korean Physical Society, Vol. 50, No. 1, January 2007, pp. 234 238 Oscillatory Turing Patterns in a Simple Reaction-Diffusion System Ruey-Tarng Liu and Sy-Sang Liaw Department of Physics,
More information