Flexible controllability of interlocked feedback loops in biological system

JOURNAL OF SCIENCE OF HNUE FIT., 2013, Vol. 58, pp. 88-101 This paper is available online at FLEXIBLE CONTROLLABILITY OF INTERLOCKED FEEDBACK LOOPS IN BIOLOGICAL SYSTEM Nguyen Cuong1∗, Hoang Do Thanh Tung1, Trinh Thi Xuan2 1Institute of Information Technology, VAST; 2Faculty of Information Technology, Hanoi Open University ∗Email: ncuong@ioit.ac.vn Abstract. The positive and negative feedback loops are ubiquitous basic elements in regulatory biological complex networks. Positive feedback loops are responsible for bistability, creating discontinuous output response from continuous input and are regarded as reliable toggle switches for making all-or-none decisions. Negative feedback loops are responsible for homeostasis and oscillations. The feedback loops hook up together and cooperate in biological complex networks rather than work alone. In this paper, we study the behavior of interlocked feedback loops. We find that interlocked positive feedback loops which are classified into two classes based on their dynamic and biological functions could significantly enhance reliability and provide flexible controllability of toggle switches. The interlocking of a positive feedback loop with a negative feedback loop brings up relaxation oscillation with amplitude and frequency that could be controlled freely and independently. Biological implications are also discussed. Keywords: Interlocked feedback loop, flexible controllability, tristability. 1. Introduction The positive and negative feedback loops are ubiquitous basic elements in the regulatory biological complex networks such as transcriptional regulation networks [1], signaling transduction pathways [2], cell cycle regulatory networks [3-5] and circadian rhythm regulation networks [6]. The positive feedback loops are responsible for bistability. This means that the positive feedback loop creates discontinuous output response from continuous input. They are regarded as toggle switches for making all-or-none decisions [7] and self-perpetuating states [2]. In certain cases, the positive feedback loop could make a one way switch or irreversible switch, and once the output turns on, it never turns off [7-9]. However, the positive feedback loop alone does not guarantee bistability and high nonlinearity within positive feedback loops is required [9]. 88 Flexible controllability of interlocked feedback loops in biological system In negative feedback loops (NFL), output P somehow suppresses its own production through a feedback loop. Typically, in negative feedback loops, output P activates its inhibitor X through a signal pathway including several components which introduce a time delay in a feedback loop. The delayed negative feedback loop (dNFL) could be able to generate oscillation [7, 10]. dNFLs are found in many biological systems such as signaling pathway [11], cell cycle regulation networks [3-5] and circadian clocks [6, 12]. However, in many biological systems, positive and negative feedbacks hook up and cooperate to work out the desired functions of the systems. For instance, the Mitotic trigger and regulation of the cell cycle of various organisms is regulated by interlocked PFL and dNFL related to the Mitotic Promoting Factor (MPF) and its friends and enemies (APC) [3, 5, 13]. The interlocked positive and delay negative feedback loops related to the CLOCK genes are also found in circadian regulation networks in a broad range of species [6]. Some studies have been made to address the role of interlocked feedback loops from perspectives such as robustness [14] [15] and noisy signal coping [16]. However, the link between biological structure and biological function has not yet been discovered. The question is, are there any new features of interlocked systems that did not exist as features of single systems? And, what are the biological advantages and disadvantages of the new features, if any? In this paper, we show that the interlocking of feedback loops makes the dynamic behaviors of single feedback loop more robust and brings new dynamic behaviors, such as tri-stability, controllability of amplitude and frequency of oscillation. Those new features are utilized in real biological processes. 2. Method 2.1. The wired diagrams In a positive feedback loop (PFL), the output somehow promotes its own production. There are two type of positive feedback loops, mutual activation (Figure 1A) and mutual inhibition (Figure 1B). In mutual activation, product P promotes its helper TF which in turn promotes the production rate of output P , a so-called self-enhancing positive feedback loop (ePFL). In mutual inhibition, output P removes its inhibitors E to release itself from suppression, a so-called self-recovering positive feedback loop (rPFL). In a negative feedback loop (NFL), the output somehow suppresses itself by activating an inhibitor X which in turn accelerates the degradation rate of output P (Figure 1C). 89 Nguyen Cuong, Hoang Do Thanh Tung, Trinh Thi Xuan Figure 1. Wired diagrams of (A) a Self-enhancing positive feedback loop, (B) a Self-recovering positive feedback loop and (C) a negative feedback loop. Solid and dashed arrows represent the reactions and regulatory effect of components, respectively 2.2. The mathematical model Output P is, in general, synthesized proportional to input S and the transcription factor TF and is degraded in proportion to inhibitor E and X and the output P itself. The general ODE for P would be dP dt = ksS + keP − (kd + krE + knX)P (2.1) Where ks and kd are synthesis and degradation rate constants. The constants ke, kr and kn represent the effect of TF, E, and X on output P and is referred to as positive and negative feedback strength. By setting kr = kn = 0, kr = kn = 0, ke = kn = 0 or ke = kr = 0, we get the equation for output P with ePFL, rPFL and NFL, respectively. The activation and inhibition of helper TF and inhibitors E and X are governed by Michaelis-Menten kinetics in order to utilize the nonlinearity[17]. However, any other reaction types possessing high nonlinearity such as the hill-coefficient equation would do. The ODE for transcription factor TF, enzyme E and inhibitorX would be dTF dt = katfP (1− TF ) Jatf + 1− TF − kitfTF Jitf + TF (2.2) dE dt = kae(1− E) Jae + 1− E − kieP ·E Jie + E (2.3) dX dt = ( kaxP (1−X) Jax + 1−X − kixX Jix +X ) 1 τ (2.4) 90 Flexible controllability of interlocked feedback loops in biological system The constants of ka and ki are maximal activation and inhibition rates, and those of Ja and Ji are Michaelis constants (relative to the total concentrations of TF). The transcription factor TF is more than half of maximum when P > θe = kitf/katf , and vice versa. The ratio θe = kitf/katf is called the activation level of PFL. The inhibitor E is more than half of maximum when P > θr = kie/kae, and vice versa. The ratio θn = kie/kae is called the activation level of PFL. The constant τ provides a time scale for the activation and inhibition of enzyme E. The larger τ is, the slower the activation and inhibition of enzyme E is. The parameter values are listed below: ks = 0.4, kp = 1, kd = 2, kn = 2, katf = 1, kitf = 0.5, kae = 1, kie = 1, kax = 0.25, kix = 0.1, Jatf = Jitf = 0.01, Jae = Jie = 0.05, Jax = Jix = 0.005, τ = 1. Changed parameters are specified in the text. 3. Results In positive feedback loops, output P somehow promotes its own production. The output can elevate itself by promoting a helper which in turn elevates the output level or releases itself from suppression by suppressing its inhibitor. 3.1. Self-enhancing positive feedback loop: signal amplifier The typical example of a self-enhancing positive feedback loop is autocatalytic regulation in which output P promotes its own production by activating activator TF, Figure 2A As shown in Figure 2B, the input-output response curves with different feedback strength are shown. As long as input S increases, product P will increase monotonically. When output P reaches the activation level of ePFL, P > re, it is able to turn on its helper, the TF. However, if there is no feedback strength, ke = 0, there is no help from TF to affect the rate of output P , thus PFL would continue to monotonically increase with input S, shown as a black, dashed-dotted line. By increasing feedback strength ke, the transcription factor TF in turn enhances the synthesis rate of product P . If feedback strength is weak, the enhancement of synthesis rate is small and thus output P monotonically increases and is enhanced due to simple regulation (red line). At stronger feedback strength, enhancement of the synthesis rate is larger and thus output P discontinuously increases, creating bistability or hysteresis (blue line). Bistability means that output P could be in an either low of high stable steady state (solid lines) at one input S. Those stable steady states are separated by an unstable steady state, represented by a dashed line. Output P abruptly changes from one stable state to other one at two critical values, eSN1 and eSN2, called saddle node bifurcation. The region bound by SNs, eSN2 < S < eSN1 is called a bistable region. If feedback strength gets too strong, enhancement of the synthesis rate is strong thus making the bistability become irreversible bistability (magenta line). That means that once output P is enhanced, it never goes down even when input S is washed out, S = 0. 91 Nguyen Cuong, Hoang Do Thanh Tung, Trinh Thi Xuan The summarization is shown in Figure 2C, in which the phase diagram of input S and feedback strength ke is plotted. The two saddle node bifurcations, eSN1 and eSN2, occur at smaller input S with stronger feedback strength. When feedback strength is smaller than a critical value (cusp point), ke < kec, there will be no bistable region with any input S. And the bistable region becomes irreversible when the feedback strength is larger than the critical value kei, ke > kei , when eSN2 goes to a negative part of input S. This ePFL is referred to as an amplification module. Input signal S is passed through and amplified. The amount of amplification is proportional to the strength of transcription factor TF, keTF. The stronger the feedback strength is, the larger the amplification is. Figure 2. Self-enhancing and Self-recovering positive feedback loop (ePFL) (A, C) Input-response curve of ePFL and rPFL with none (dashed-dotted), weak (red), strong (blue) and very strong (magenta) feedback strength, respectively: The solid and dashed black lines represent stable and unstable steady state of output P where eSN1 and eSN2 are saddle node bifurcation of ePFL. (B, D) The Input-feedback strength phase diagrams show regions of monostability and bistability. 3.2. Self-recovering Positive feedback loop: signal buffer Typical examples of a self-recovering positive feedback loop are mutual inhibition and double negative feedback loops, in which output P removes its suppression by inhibiting suppressor enzyme E, Figure 2D. In this case, output P is synthesized proportional to input S. The degradation rate is proportional to output P and the effect of suppressor E, (kd + krE)P , where 92 Flexible controllability of interlocked feedback loops in biological system kd is degradation rate, kr represents the effect of enzyme E on the degradation rate, the so-called feedback strength. Note that transcription factor E is suppressed when P > rr = kar/kir, and vice versa. See Figure 2E. The ratio rr = kar/kir is called the deactivation level of rPFL. The input-output response curves with different feedback strength kr are shown in Figure 2E. As long as input S increases, product P is initially suppressed due to the suppression of enzyme E, which is proportional to feedback strength kr. If there is no feedback strength, kr = 0, there is no suppression and thus output P monotonically increases as that of simple regulation, shown as a dashed-dotted line. With the presence of feedback strength, the output increases but is suppressed. When output is strong enough to remove its suppression, P > rr, the enzyme is removed and the output is released. If feedback strength is weak, the suppression on P is small and the recovery is continuous, shown as a red line. Output P would be discontinuously and abruptly recurred (bistability) due to the simple regulation once the feedback is stronger, shown as a blue line. Obviously, feedback strength represents the suppression of P from its inhibitor. The output is more suppressed with stronger feedback strength and therefore more input S is needed to produce output P . The ratio rr = kar/kir determines when the inhibitor is removed. A summarization is shown in Figure 2F, in which the phase diagram of input S and feedback strength kr is plotted. As long as feedback strength increases, it is hard to remove the suppression on inhibiting enzyme E and thus more input S is needed to remove E. The two saddle node bifurcations move to a larger input S. Note that when feedback strength is smaller than the critical value (cusp point), kr < krc, there will be no bistable region with any input S. This ePFL works as a buffer. When the input signal increases, the output response is increased but buffered in an inactive state until it overrides its suppressor. Thus, the output response is delayed from input S. The delay strongly depends on suppressor capacity (feedback strength) kr. The stronger the suppressor is, the more output is delayed by P. 3.3. The interlocking of two positive feedback loops: monostability, bistability, and tristability The interlocking of two positive feedback loops which show monostability and bistability would generate monostability, bistability and multistability. Obviously, each positive feedback loop has two important control parameters - feedback strength and the activation level. Feedback strength determines how much the output is enhanced (ke in ePFL) or suppressed (kr in rPFL). The activation levels, re and rr, determine the working region of ePFL and rPFL with level of output P , respectively. Therefore, there are two possibilities that the feedback loop in ePFL and rPFL will be activated far from each other, (i) re rr, or (iii) the feedback loop in ePFL and rPFL can be cooperatively activated, re ≈ rr, 93 Nguyen Cuong, Hoang Do Thanh Tung, Trinh Thi Xuan (i) First, consider the case that the ePFL activate at a smaller range of output P than that in rPFL, re < rr. (See Figure 3A1-Figure 3A4) Output P in iPFL (shown in black) is initially suppressed from simple regulation (shown as a dashed-dotted line) under the effect of inhibiting enzyme E, thus output P in iPFL initially follows that of rPFL (shown in blue). When outputP in iPFL reaches the activation level of ePFL, P > re, the TF turns on to enhance output P . The output P in iPFL continues to increase due to the activation level of rPFL, P > rr, at which time the inhibiting enzyme E is removed to release output P . Therefore, there are two transitions of output P corresponding to the different working region of ePFL and rPFL. The transition might be continuous (monostability) or discontinuous (bistability) depending on the feedback strength. Both transitions are continuous when feedback strengths ke and kr are weak (See Figure 3A1) One of them will be weak and the other will be strong and cause a continuous and a discontinuous transition. (See Figure 3A2,3) If feedback strengths ke and kr are strong, there will be two discontinuous transitions. (See Figure 3A4) The two discontinuous transitions might or might not overlap. Once they overlap, mutltistability emerges. Mutltistability implies that at one input S there are more than two stable steady states separated by two unstable steady states. Here, three stable steady states, a low, middle and high level of P , are separated by two unstable steady states. (ii) Consider the case that the ePFL activate at a larger range of output P than in rPFL, re > rm. (See Figure 3C1-Figure 3C4) This means that as long as output P in iPFL (shown in black) increases and reaches an activation level of rPFL, P > rm, the inhibiting enzyme E is removed and output P reaccurs as simple regulation (shown as a black, dashed-dotted line). Output P in iPFL keep increasing until it reaches the activation level of ePFL, P > re„ and the TF turns on to enhance output P from simple regulation. Therefore, there are two transitions of output P with corresponding to different working region of ePFL and rPFL. Depending the feedback strengths ke and kr, the input-output response of P might have two continuous transitions (See Figure 3C1), a continuous and a discontinuous transitions (See Figure 3C2,3), or two discontinuous transitions (See Figure 3C4). Here we show a case where the two discontinuous transitions do not overlap and there is thus no mutltistability. However, mutltistability could be obtained easily by tuning the activation level or feedback strength. (iii) When the ePFL and rPFL activate at the same range of output P , re ≈ rm, activation of ePFL causes activation of rPFL, and vice versa. Thus they might cooperate with each other. An interlocking of two PFLs with weak feedback strengths which show monostability could generate bistability (See Figure 3B1) With an interlocking of two PFL with weak feedback strength and strong feedback strength, the bistable region of iPFL is moved and extended from that of a single PFL. For instance, as seen in Figure 3B2, the bistable region of rPFL is moved and extended to the left under the effect of ePFL. And the bistable region of ePFL is moved and extended to the right under the effect of rPFL (See Figure 3B3) An interlocking of two PFL with strong feedback strengths in the bistable region 94 Flexible controllability of interlocked feedback loops in biological system of iPFL is extended in both the left and right directions from that of a single PFL. (See Figure 3C4) Figure 3. (A1-A4, B1-B4, C1-C4) input-output response of output P at different activation levels of ePFL and rPFL with different feedback strengths as shown Red, blue and black lines represent the input-output response of a single ePFL, rPFL and the interlocking of ePFL and rPFL (iPFL). The dashed-dotted line represents the input-output response curve of simple regulation. For ePFL feedback strength, ke = 0.3 is weak and 1 is strong. For rPFL feedback strength kr = 0.5 is weak while 1.5 is strong. For the A’s: re = 0.5, rr = 0.8; B′s : re = rr = 0.5; C ′s : re = 1, rr = 0.5. 3.4. The interlocking of a positive and a negative feedback loop: Relaxation oscillation with flexible controllability of amplitude and frequency Oscillation comes from the cooperation of positive and negative feedback loops. The negative feedback loop drives output P back and forth between low and high steady states which are from a positive feedback loop. (See Figure 4A) Therefore, oscillation properties such as amplitude and frequency strongly depend on the feedback strengths of and the embedded time delay in feedback loops. Indeed, a too weak or too strong positive feedback strength pushes the system toward a stable fixed point. In Figure 4C, D increases feedback strength kp and elevates the right branch of a reverted-N shape null-cline of output P (solid line) causing the oscillatory system to move toward the stable steady state. If the positive feedback strength 95 Nguyen Cuong, Hoang Do Thanh Tung, Trinh Thi Xuan is too low, the iPNFL turns out to be a single NFL and thus the iPNFL oscillation will be shut down. The bifurcation diagram in Figure 6A shows that with either a too weak (ke < ke1) or a too strong (ke2 <