- Open Access
Computational analysis reveals the coupling between bistability and the sign of a feedback loop in a TGF-β1 activation model
BMC Systems Biology volume 11, Article number: 136 (2017)
Bistable behaviors are prevalent in cell signaling and can be modeled by ordinary differential equations (ODEs) with kinetic parameters. A bistable switch has recently been found to regulate the activation of transforming growth factor-β1 (TGF-β1) in the context of liver fibrosis, and an ordinary differential equation (ODE) model was published showing that the net activation of TGF-β1 depends on the balance between two antagonistic sub-pathways.
Through modeling the effects of perturbations that affect both sub-pathways, we revealed that bistability is coupled with the signs of feedback loops in the model. We extended the model to include calcium and Krüppel-like factor 2 (KLF2), both regulators of Thrombospondin-1 (TSP1) and Plasmin (PLS). Increased levels of extracellular calcium, which alters the TSP1-PLS balance, would cause high levels of TGF-β1, resembling a fibrotic state. KLF2, which suppresses production of TSP1 and plasminogen activator inhibitor-1 (PAI1), would eradicate bistability and preclude the fibrotic steady-state. Finally, the loop PLS − TGF-β1 − PAI1 had previously been reported as negative feedback, but the model suggested a stronger indirect effect of PLS down-regulating PAI1 to produce positive (double-negative) feedback in a fibrotic state. Further simulations showed that activation of KLF2 was able to restore negative feedback in the PLS − TGF-β1 − PAI1 loop.
Using the TGF-β1 activation model as a case study, we showed that external factors such as calcium or KLF2 can induce or eradicate bistability, accompanied by a switch in the sign of a feedback loop (PLS − TGF-β1 − PAI1) in the model. The coupling between bistability and positive/negative feedback suggests an alternative way of characterizing a dynamical system and its biological implications.
Bistability has been found in many biological systems, and mathematical models using ordinary differential equations (ODEs) with respect to time provide a good representation for studying these dynamic behaviors [1,2,3,4]. Some bistable systems show a binary behavior at the level of single cells but exhibit a graded response for a population of cells , while other bistable systems cause a binary output even for populations of cells . Different strategies have been developed to study the two different categories of bistable systems. For systems with bistability at the level of single cells, the bistable nature of the output is often well-known, guiding researchers to elucidate the underlying molecular circuit that enables the bistability. For systems with population-level bistability, the bistability may not be obvious by inspection, in which case, laborious experimental measurement is required to establish the phenomenon . In a recent study, we modeled the bistable activation of transforming growth factor-β1 (TGF-β1), which belongs to the category of population-level bistability. The bistability was validated experimentally by showing hysteresis in an in vitro system . TGF-β1 is a cytokine with broad importance for contexts such as cancer, liver cirrhosis, wound healing and regeneration. Our model captures multiple pathways with positive and negative effects towards TGF-β1 and its activators (Fig. 1; Additional file 1: Supplementary Note 1).
In the current study, we explore some possible triggers that could influence the bistable transition, and look for system properties that correlate with bistability. We demonstrate that factors like calcium and Krüppel-like factor 2 (KLF2) (Additional file 1: Supplementary Note 2) level can be modeled implicitly into the reaction rate parameters of the model, allowing us to discuss the on and off transitions of bistability under different conditions. We identified a coupling between the sign of a signaling loop in the model (i.e., whether the signaling loop shows positive feedback or negative feedback) and the presence of bistability in the model. This could suggest alternative ways to identify and validate of systems with population-level bistability.
Calcium and KLF2 have potential influence on the steady state of TGF-β1 activation
Calcium would promote the steady state with high TGF-β1 activation
We built the low and high calcium variants of the model by considering the potential effects of calcium on the PLS-TSP1 interaction (see Methods, Fig. 1 red arrows). Three parameters for the calcium effect were not known quantitatively and were estimated (see Methods). In Fig. 2a-b, we simulated the low calcium model and the high calcium model over time with 27 total initial configurations. These 27 configurations were combinations of 3 initial concentrations for each of TGF-β1, TSP1, and plasmin, the ssT level, ssP level, and the mean level of ssT and ssP (27 = 33). The initial concentrations of other species were set to the average of their two steady state levels, (i.e., 0.5ssT + 0.5ssP). In the low-calcium model, all trajectories converged to ssP with low TGF-β1, but in the high-calcium model, several of the initial configurations converged to ssT with high TGF-β1. To generalize our understanding of this effect, we plotted the boundary (the separatrix, Fig. 2c) between the initial configurations that caused convergence toward ssT (red) and the initial conditions that caused convergence toward ssP (blue). Initial concentrations were constants for all species other than PLS and TSP1. By comparing the separatrix of the low calcium model (dot) and the high calcium model (circle), we observe a shift of the separatrix toward the blue (ssP) region. This means the red (ssT) region is enlarged in high calcium environment (arrow 3). As expected, calcium tips the balance between PLS and TSP1 to achieve a significant effect on steady state of TGF-β1 activation.
KLF2 would eliminate the steady state with high TGF-β1 activation
KLF2 is a transcription factor studied extensively in atherosclerosis and fibrosis, and previous studies of KLF2 signaling showed TSP1 and PAI1 (plasminogen activator inhibitor-1) to be two of its most strongly affected targets [6, 7]. To study how KLF2 would affect bistable activation of TGF-β1, variants of the TGF-β1 activation model were built as described in methods (Fig. 1 blue arrows). We built a model called “100% KLF2” that downregulated the TSP1 production and PAI1 production rates, proportional to the published effects of KLF2 on the mRNA levels of TSP1 (−7.8 fold) and PAI1 (−7.4 fold). This is a strong effect, so we also built models with 90%, 80%, … 10%, and 0% of the KLF2 effects on the TSP1 and the PAI1 production rates. Each model in the series was simulated to obtain the steady state concentrations. When a dynamical system is bistable, its two steady states are commonly obtained by simulating the model twice, once starting from each side of the separatrix boundary (for example, initializing the system with opposite extreme levels of TGF-β1). For our series of models, the steady states obtained after initialization with high TGF-β1 (resembling ssT) were plotted with open red boxes, and the steady states obtained after initialization with low TGF-β1 (resembling ssP) were plotted with solid blue circles (Fig. 2d). Models for each level of KLF2 were plotted in terms of PAI1 and TGF-β1 steady states (with KLF2 levels decreasing from left to right). For models with KLF2 ≥ 40%, the open red boxes fell at the same points as the solid blue circles, indicating they are monostable. For models with KLF2 ≤ 30%, the red open boxes were distinct from the blue circles, indicating two steady states. For KLF2 levels from 0% to 100% (right to left), the low TGF-β1 steady state (blue dots) remained almost constant, while the high TGF-β1 steady state (red boxes) merged with the low TGF-β1 steady state in an ultrasensitive manner when the KLF2 effect increased from 30% to 40%.
Bifurcation analysis studies how parameter change affects the qualitative behavior and the steady states of a system . Bifurcation plot allows us to see all the equilibria of the system and how the equilibria vary with change of KLF2 levels and other rate parameters. We chose one parameter named “keff2” to show the steady state behavior of the low KLF2 and high KLF2 system. “Keff2”, the enzymatic efficiency of plasmin, is one of many rate parameters that affect the overall bistability of the system. Figure 2e shows the bifurcation plot for the 0% KLF2 model and the 100% KLF2 model with respect to the parameter “keff2”. Solids lines represent the stable steady states of the system (ssP or ssT). Dotted lines represent one unstable steady state between the two stable steady states, which is not achievable through simulation. Black circles represent two limit points in the bifurcation curve, which separates the monostable regime and bistable regime. Bifurcation analysis confirmed that the system with 0% KLF2 retained bistability (an “S-shaped” curve in Fig. 2e) while the system with 100% KLF2 was monostable.
The bistability of the system correlates with the sign of the PLS-PAI1 feedback loop
Positive and negative feedback loops are ubiquitous in biological systems , and necessary for many functions [10, 11]. The TGF-β1 activation network is composed of multiple overlapping feedback loops, including two feedback loops between PLS and PAI1. One obvious loop is the negative feedback loop PLS → TGF-β1 → PAI1⊣PLS, which is frequently cited [12,13,14,15,16]. A less obvious loop is PLS⊣TSP1 → TGF-β1 → PAI1⊣PLS, with two inhibitory effects, meaning positive feedback (See Fig. 3. a-b). Interestingly, experiments have already observed two opposite behaviors of PLS towards TGF-β1 and PAI1 [17,18,19], giving indirect evidence for the possibility of both positive and negative feedback loops involving PLS, TGF-β1, and PAI1.
To characterize the feedback between PLS and PAI1 in this network, we plotted feedback behavior in both low KLF2 (0%) and high KLF2 (100%) models (Fig. 3a). Since PAI1 (plasminogen activator inhibitor-1) is antagonistic towards PLS, the sign of the PLS-PAI1 feedback loop is determined by the response of PAI1 to PLS. We used a stepwise PLS input (black curve) to perturb both the low (%0) KLF2 model and the high (100%) KLF2 model and we simulated the response of PAI1. The low KLF2 model showed a decrease in the level of PAI1, which means that the overall PLS-PAI1 feedback is dominated by the PLS⊣TSP1 → TGF-β1 → PAI1⊣PLS double negative (positive) feedback loop. The high KLF2 model showed an increase in the level of PAI1, which means that the overall PLS-PAI1 feedback is dominated by the PLS → TGF-β1 → PAI1⊣PLS negative feedback loop.
It was interesting to observe in simulation that KLF2 not only was able to eliminate one of the steady states and turned the system into monostable, but also able to change the sign of PLS-PAI1 feedback loop. To further characterize the effect of KLF2 on the system and discuss the reason behind it, we plotted the 2d bistable region of the system in the kp1-kp2 phase plane. kp1 is the TSP1 synthesis rate parameter and kp2 is the PAI1 synthesis parameter, both of which affect the bistability of the system. We found the boundaries of the bistable region through equilibrium continuation of kp2 for a series of kp1 values (Fig. 3b). In our model, KLF2 is represented as a combination of fold changes of kp1 and kp2, therefore, KLF2 levels can be represented as a series of points in the kp1-kp2 phase plane (Fig. 3c). It can be seen on Fig. 3c that KLF2 point is moving out of the bistable region when KLF2 increases from 0 to 100%. We then analyzed the sign of PLS-PAI1 feedback in this kp1-kp2 phase plane. Interestingly, there is a large overlap between the bistable region of the system and positive feedback region of PLS-PAI1 feedback loop (Fig. 3d). KLF2 = 0 point lies in the overlapping area of the bistable region and PLS-PAI1 positive feedback region, while KLF2 = 100% lies in the overlapping area of the non-bistable region and PLS-PAI1 negative feedback region. This explains why the change of KLF2 level can have two different effects on the system.
The large overlap between bistable region and PLS-PAI1 positive feedback region is also an interesting property of the system, since it suggests that in reality, a bistable TGF-β1 activation system most likely also has positive PLS-PAI1 feedback. Although negative feedback loop between PLS-PAI1 has been observed repeatedly and is well accepted [17, 18], positive feedback would be novel. We tested the sign of the feedback from PLS to PAI1 using an experimental system known to exhibit TGF-β1 bistability , a cell culture model of liver fibrosis. In this co-culture with primary hepatocytes and HSC-T6 cell lines, we added different levels of PLS and we measured PAI1 mRNA levels using RT-PCR (Fig. 3b). Increasing PLS was found to cause decreased expression of PAI1 in this bistable system, implying that PLS and PAI1 can indeed exhibit positive feedback.
Using a model of TGF-β1 activation, we explored the ability of external factors to switch bistability on and off, and we characterized the correlation between the bistability of the model, and the sign of a feedback loop in the network.
For the first part, we modeled known effects of calcium on the balance between TSP1 and PLS [20,21,22,23,24,25], and known effects of KLF2 on the gene expression of PAI1 and TSP1 . We then used modeling to show how these effects would propagate through the system. Specifically, the model predicted that calcium would significantly promote TGF-β1 activation, shifting the bistable threshold of the system. The calcium-induced increase in TSP1 would lie within the physiological range of TSP1 . Literature search reveals that extracellular calcium may be relatively easy to perturb via biomaterials of bandages, etc. Therefore, the effect of extracellular calcium on TGF-β1 might have important therapeutic implications for fibrotic or inflammatory diseases where abnormal TGF-β1 contributes to disease. For example, fibrotic diseases are driven by high levels of TGF-β1 , and therapeutic studies in animals have achieved significant access by increasing the PLS pathway  or decreasing the TSP1 pathway . In our model, if we take the ssP state to be healthy and the ssT state to be fibrotic, then a fibrotic system with high calcium could transition toward health through an increase of PLS (Fig. 2c, arrow 1), through a decrease of TSP1 (Fig. 2c, arrow 2), or a combination of both (Fig. 2c, arrow 3).
In contrast, KLF2 was simulated to increase PLS activity and decrease the levels of TGF-β1, by suppressing PAI1 and TSP1 expression. This is consistent with previous work with statin drugs on liver fibrosis , where KLF2 upregulation was observed after treatment with simvastatin. Our model predicts that one of the ways KLF2 may contribute to improvement of liver fibrosis may be by decreasing the activation of TGF-β1 through reduction of the TSP1 and PAI feedback effects.
While modeling KLF2 effects, we noticed that loss of bistability also caused a change in the sign of the PLS - PAI1 feedback loop. Without KLF2, the PLS-PAI1 feedback loop was positive (double negative), but with KLF2 (100% KLF2) and with the destruction of bistability, the PLS-PAI1 feedback loop was negative. Additional bifurcation analysis revealed that high KLF2 is a special case of the general observation, that the bistability of the system is correlated with the sign of the PLS-PAI1 feedback loop.
We demonstrated the positive feedback behavior between PLS and PAI1 in the bistable TGF-β1 system, using an in vitro experiment. Previous studies have already revealed two opposite behaviors of PLS towards TGF-β1 and PAI1 [17,18,19], providing indirect evidence for the possibility of both positive and negative feedback loops. Some aspects of the feedback loop are relatively unambiguous. For example, PAI1 is a specific and potent inhibitor of plasmin activation. PAI1 production follows TGF-β1 signaling so closely that, in practice, PAI1 levels are commonly measured as a readout of TGF-β1 activation . The behavior of the feedback loop thus boils down to the behavior of the PLS – TGF-β1 relationship. In isolation, PLS clearly is able to activate TGF-β1. The same effect has frequently been observed in more physiological contexts, and there is considerable published evidence that PLS and/or plasminogen activators can cause an increase in TGF-β1 and/or PAI1 levels [12,13,14,15,16]. This positive effect of PLS on TGF-β1 or PAI1 serves as evidence that the loop between PLS and PAI1 can have negative feedback. Although the activating ability of PLS toward TGF-β1 is well known and accepted, some studies also suggest the opposite effect. For example, PLS caused TGF-β2 levels to decline in breast cysts . Furthermore, one unconventional finding by Seo et al., showed a positive feedback effect between PAI1 and TGF-β1 , suggesting that PLS can cause a decrease in TGF-β1. In sum, we found that there is some support in the published literature for our prediction that the relationship between PLS and PAI can show either negative feedback or positive feedback, depending on context.
Another important consideration in interpreting this model is the redundancy of proteases and matrix factors that play roles similar to PLS or TSP1. Actually, PLS is only one of many proteases (including elastase, MMP-2, MMP-9, ADAMTS1 and others) that can both activate TGF-β1 and cleave TSP1 [33,34,35,36]. Meanwhile TSP1 can inhibit many of these proteases [20, 37,38,39]. Extracellular proteases often function interdependently by activating each other (e.g. PLS activates several MMPs, which activate other MMPs [40, 41]), and some proteases may have partially redundant effects. Likewise TSP1 may represent a larger class of matrix proteins and mechanical factors with redundant roles in this model. Fibrillin and LTBP1 can promote TGF-β1 activation [42,43,44,45], as can factors that create mechanical tension in the matrix [46, 47]. Fibrillin and LTBP1 can be cleaved by PLS [14, 48], and even mechanical tension would be antagonized by PLS cleavage. In other words, PLS and TSP1 are archetypes of two larger classes of effects, a protease category and a matrix category, that may be capable of antagonizing the effect of each other, even as they contribute individually to TGF-β1 activation. The redundancy of the protease-versus-matrix competition suggests that this antagonism may be an organizing principle of TGF-β1 regulation, with evolutionary importance to the organism. On the other hand, this redundancy also creates many complexities that could perturb the phenomena we simulated. For example, the effects we attribute to PLS itself may actually result from the indirect effects of PLS-activated proteases. Thus, we speculate that the insights we drew from this model point to important properties of general TGF-β1 activation regulatory networks in different contexts.
Using TGF-β as a case study, we demonstrated that external factors could influence the bistability of the model, and that these influences can be modeled implicitly using the reaction parameters of the model. Furthermore, we showed that system-level properties like the sign of feedback loops can correlate with the bistability of a complex model. This provides a novel characterization of the transition between bistable and monostable regimes, and provides a non-obvious explanation for seemingly contradictory experimental findings about the contribution of the PAI1 - PLS feedback loop toward TGF-β activation.
TGF-β1 bistable activation model
We used the model built by Venkatraman et al. as the base model of TGF-β1 regulation. We increased the “kothers” parameter relative to the published model, to allow for higher basal activation of TGF-β1 by other activators such as integrins [42, 49, 50]. Simulations were performed using kroneckerbio toolbox  and the ode15s solver in MATLAB (Mathworks, Natick, MA).
Calcium can affect the structure of TSP1 [22,23,24,25], the enzymatic activity of PLS cleaving TSP1 , and the ability of TSP1 to inhibit serine protease activity . These effects were represented by the rate constants k3, k_3 and k4 in the TGF-β1 activation model. We used the original parameter settings as the low calcium settings. To reflect high calcium conditions, we increased k3 by 10 fold, decreased k_3 by 0.1 fold, and k4 by 0.0001 fold, in order to reflect a high level of calcium in the environment. Details of the model can be found in Table S1 (Additional file 2) and Table S2 (Additional file 3).
It has been shown that KLF2 can decrease TSP1 expression by 7.8 fold and PAI1 expression by 7.4 fold [6, 7]. We simulated the TGF-β1 activation model with no change (0% of the KLF2 effect, original parameter settings), with 100% of the KLF2 effect (7.4 fold decrease of PAI1 synthesis parameter kp2 and 7.8 fold decrease of TSP1 synthesis parameter kp1), as well as a series of intermediate models with 10%, 20%, … 90% of the KLF2 effect, causing intermediate levels of decrease in the PAI1 and TSP1 synthesis rates.
Isolation of primary hepatocytes was performed on male Wistar rats (250-300 g), via a two-step collagenase perfusion method as described previously . A co-culture model of primary rat hepatocytes and hepatic stellate cell line T6 (HSC-T6) was established as described in . Briefly, primary rat hepatocytes were first seeded at a density of 2 × 105 cells on 35 mm collagen-coated dishes (IWAKI) using Williams’s E media with 10% FBS. After 4 h, hepatic stellate cell line T6 (HSC-T6) was seeded at a density of 1.4 × 106 cells. The cells were cultured overnight in 35 °C, and 5% CO2 in William’s E media with 2% FBS to facilitate HSC activation. The next day media was changed to Williams’s E without serum, along with different doses of PLS. After 24 h, the cells were collected.
RT-PCR was performed as described in . Briefly, mRNA was isolated from the cells using RNeasy mini kit (Qiagen), and its concentration was quantified using a Nanodrop 2000 UV-Vis Spectrophotometer. One microgram of mRNA from each sample was converted to cDNA (Invitrogen, Superscript Reverse Transcriptase III) and real-time PCR reaction (Roche, Sybr Green Master mix) was carried out for plasminogen activator inhibitor-1 (PAI1) and β-actin, with in-house primers shown in Table S3 (Additional file 4). The gene expression values were determined by the Del-Del CT relative quantitation method; the target CT values were normalized to the endogeneous reference β-actin, and the normalized mRNA was expressed as a fold-change relative to the untreated control.
Plasminogen activator inhibitor-1
PLS predominant steady state
TSP1 predominant steady state
Transforming growth factor β1
Albeck JG, Burke JM, Spencer SL, Lauffenburger DA, Sorger PK. Modeling a snap-action, variable-delay switch controlling extrinsic cell death. PLoS Biol. 2008;6:2831–52.
Shi Y, Mellier G, Huang S, White J, Pervaiz S, Tucker-Kellogg L. Computational modelling of LY303511 and TRAIL-induced apoptosis suggests dynamic regulation of cFLIP. Bioinformatics. 2013;29:347–54.
Chatterjee A, Kaznessis YN, W-S H. Tweaking biological switches through a better understanding of bistability behavior. Curr Opin Biotechnol. 2008;19:475–81.
Sobie EA. Bistability in Biochemical Signaling Models. Sci Signal. 2011;4:tr10.
Venkatraman L, Chia S-M, Narmada BC, White JK, Bhowmick SS, Forbes Dewey C Jr, So PT, Tucker-Kellogg L, Yu H. Plasmin triggers a switch-like decrease in Thrombospondin-dependent activation of TGF-β1. Biophys J. 2012;103:1060–8.
Dekker RJ. KLF2 provokes a gene expression pattern that establishes functional quiescent differentiation of the endothelium. Blood. 2006;107:4354–63.
Boon RA, Fledderus JO, Volger OL, van Wanrooij EJA, Pardali E, Weesie F, Kuiper J, Pannekoek H, Dijke Ten P, AJG H. KLF2 suppresses TGF-beta signaling in endothelium through induction of Smad7 and inhibition of AP-1. Arterioscler Thromb Vasc Biol. 2007;27:532–9.
Strogatz SH. Nonlinear dynamics and chaos. Westview Press. 2008;352
Neves S, Iyengar R. Modeling of signaling networks - Neves - 2002 - BioEssays - Wiley online library. BioEssays. 2002;
Ratushny AV, Saleem RA, Sitko K, Ramsey SA, Aitchison JD. Asymmetric positive feedback loops reliably control biological responses. Mol Syst Biol. 2012;8:577.
Tsai TY-C, Choi YS, Ma W, Pomerening JR, Tang C, Ferrell JE. Robust, tunable biological oscillations from interlinked positive and negative feedback loops. Science. 2008;321:126–9.
Grainger DJ, Wakefield L, Bethell HW, Farndale RW, Metcalfe JC. Release and activation of platelet latent TGF-beta in blood clots during dissolution with plasmin. Nat Med. 1995;1:932–7.
Khalil N, Corne S, Whitman C, Yacyshyn H. Plasmin regulates the activation of cell-associated latent TGF-beta 1 secreted by rat alveolar macrophages after in vivo bleomycin injury. Am J Respir Cell Mol Biol. 1996;15:252–9.
Yee JA, Yan L, Dominguez JC, Allan EH, Martin TJ. Plasminogen-dependent activation of latent transforming growth factor beta (TGF beta) by growing cultures of osteoblast-like cells. J Cell Physiol. 1993;157:528–34.
Godár S, Horejsi V, Weidle UH, Binder BR, Hansmann C, Stockinger H. M6P/IGFII-receptor complexes urokinase receptor and plasminogen for activation of transforming growth factor-beta1. Eur J Immunol. 1999;29:1004–13.
Binder BR, Christ G, Gruber F, Grubic N, Hufnagl P, Krebs M, Mihaly J, Prager GW. Plasminogen activator inhibitor 1: physiological and pathophysiological roles. News Physiol Sci. 2002;17:56–61.
Thornton AJ, Bruzdzinski CJ, Raper SE, Gelehrter TD. Plasminogen activator inhibitor-1 is an immediate early response gene in regenerating rat liver. Cancer Res. 1994;54:1337–43.
Flaumenhaft R, Abe M, Mignatti P, Rifkin DB. Basic fibroblast growth factor-induced activation of latent transforming growth factor beta in endothelial cells: regulation of plasminogen activator activity. J Cell Biol. 1992;118:901–9.
Seo JY, Park J, Yu MR, Kim YS, Ha H, Lee HB. Positive feedback loop between plasminogen activator Inhibitor-1 and transforming growth factor-Beta1 during renal fibrosis in diabetes. Am J Nephrol. 2009;30:481–90.
Hogg PJ, Stenflo J, Mosher DF. Thrombospondin is a slow tight-binding inhibitor of plasmin. Biochemistry. 1992;31:265–9.
Gotis-Graham I, Hogg PJ, McNeil HP. Significant correlation between thrombospondin 1 and serine proteinase expression in rheumatoid synovium. Arthritis Rheum. 1997;40:1780–7.
Calzada MJ, Kuznetsova SA, Sipes JM, Rodrigues RG, Cashel JA, Annis DS, Mosher DF, Roberts DD. Calcium indirectly regulates immunochemical reactivity and functional activities of the N-domain of thrombospondin-1. Matrix Biol. 2008;27:339–51.
Dixit VM, Galvin NJ, O'Rourke KM, Frazier WA. Monoclonal antibodies that recognize calcium-dependent structures of human thrombospondin. Characterization and mapping of their epitopes. J Biol Chem. 1986;261:1962–8.
Misenheimer TM, Mosher DF. Calcium ion binding to thrombospondin 1. J Biol Chem. 1995;270:1729–33.
Slane JM, Mosher DF, Lai CS. Conformational change in thrombospondin induced by removal of bound Ca2+. A spin label approach. FEBS Lett. 1988;229:363–6.
Bergseth G, Lappegård KT, Videm V, Mollnes TE. A novel enzyme immunoassay for plasma thrombospondin. Comparison with beta-thromboglobulin as platelet activation marker in vitro and in vivo. Thromb Res. 2000;99:41–50.
Branton MH, Kopp JB. ScienceDirect.com - microbes and infection - TGF-β and fibrosis. Microbes Infect. 1999;1:1349–65.
Wang H, Zhang Y, Heuckeroth RO. PAI-1 deficiency reduces liver fibrosis after bile duct ligation in mice through activation of tPA. FEBS Lett. 2007;581:3098–104.
Hugo C. The thrombospondin 1-TGF-beta axis in fibrotic renal disease. Nephrol Dial Transplant. 2003;18:1241–5.
Marrone G, Russo L, Rosado E, Hide D, García-Cardeña G, García-Pagán JC, Bosch J, Gracia-Sancho J. The transcription factor KLF2 mediates hepatic endothelial protection and paracrine endothelial-stellate cell deactivation induced by statins. J Hepatol. 2012;
Mazzieri R, Munger JS, Rifkin DB. Measurement of active TGF-beta generated by cultured cells. Methods Mol Biol. 2000;142:13–27.
Erbas H, Erten O, Irfanoglu ME. Breast cyst fluid plasmin activity and its effect on TGF-beta2 activation. Cancer Investig. 2008;26:22–7.
Bourd-Boittin K, Bonnier D, Leyme A, Mari B, Tuffery P, Samson M, Ezan F, Baffet G, Theret N. Protease profiling of liver fibrosis reveals the ADAM metallopeptidase with thrombospondin type 1 motif, 1 as a central activator of transforming growth factor beta. Hepatology. 2011;54:2173–84.
Sato Y, Rifkin DB. Inhibition of endothelial cell movement by pericytes and smooth muscle cells: activation of a latent transforming growth factor-beta 1-like molecule by plasmin during co-culture. J Cell Biol. 1989;109:309–15.
Yu Q, Stamenkovic I. Cell surface-localized matrix metalloproteinase-9 proteolytically activates TGF-beta and promotes tumor invasion and angiogenesis. Genes Dev. 2000;14:163–76.
Prudova A, Auf Dem Keller U, Butler GS. Overall CM: multiplex N-terminome analysis of MMP-2 and MMP-9 substrate degradomes by iTRAQ-TAILS quantitative proteomics. Mol Cell Proteomics. 2010;9:894–911.
Lawler PR, Lawler J. Molecular basis for the regulation of angiogenesis by thrombospondin-1 and -2. Cold Spring Harb Perspect Med. 2012;2:a006627.
Hogg PJ, Owensby DA, Chesterman CN. Thrombospondin 1 is a tight-binding competitive inhibitor of neutrophil cathepsin G. Determination of the kinetic mechanism of inhibition and localization of cathepsin G binding to the thrombospondin 1 type 3 repeats. J Biol Chem. 1993;268:21811–8.
Bein K, Simons M. Thrombospondin type 1 repeats interact with matrix metalloproteinase 2. Regulation of metalloproteinase activity. J Biol Chem. 2000;275:32167–73.
Davis GE, Pintar Allen KA, Salazar R, Maxwell SA. Matrix metalloproteinase-1 and -9 activation by plasmin regulates a novel endothelial cell-mediated mechanism of collagen gel contraction and capillary tube regression in three-dimensional collagen matrices. J Cell Sci. 2001;114:917–30.
Curry TE, Osteen KG. The matrix metalloproteinase system: changes, regulation, and impact throughout the ovarian and uterine reproductive cycle. Endocr Rev. 2003;24:428–65.
Annes JP, Chen Y, Munger JS, Rifkin DB. Integrin alphaVbeta6-mediated activation of latent TGF-beta requires the latent TGF-beta binding protein-1. J Cell Biol. 2004;165:723–34.
Kaartinen V, Warburton D. Fibrillin controls TGF-beta activation. Nat Genet. 2003;33:331–2.
Tritschler I, Gramatzki D, Capper D, Mittelbronn M, Meyermann R, Saharinen J, Wick W, Keski-Oja J, Weller M. Modulation of TGF-beta activity by latent TGF-beta-binding protein 1 in human malignant glioma cells. Int J Cancer. 2009;125:530–40.
Chaudhry SS, Cain SA, Morgan A, Dallas SL, Shuttleworth CA, Kielty CM. Fibrillin-1 regulates the bioavailability of TGFbeta1. J Cell Biol. 2007;176:355–67.
Buscemi L, Ramonet D, Klingberg F, Formey A, Smith-Clerc J, Meister J-J, Hinz B. The single-molecule mechanics of the latent TGF-β1 complex. Curr Biol. 2011;21:2046–54.
Ahamed J, Janczak CA, Wittkowski KM, Coller BS. PLOS ONE: in vitro and in vivo evidence that Thrombospondin-1 (TSP-1) contributes to stirring- and shear-dependent activation of platelet-derived TGF-β1. PLoS One. 2009;4:e6608.
Kirschner R, Hubmacher D, Iyengar G, Kaur J, Fagotto-Kaufmann C, Brömme D, Bartels R, Reinhardt DP. Classical and neonatal Marfan syndrome mutations in fibrillin-1 cause differential protease susceptibilities and protein function. J Biol Chem. 2011;286:32810–23.
Annes JP, Munger JS, Rifkin DB. Making sense of latent TGFbeta activation. J Cell Sci. 2003;116:217–24.
Fontana L, Chen Y, Prijatelj P, Sakai T, Fässler R, Sakai LY, Rifkin DB. Fibronectin is required for integrin alphavbeta6-mediated activation of latent TGF-beta complexes containing LTBP-1. FASEB J. 2005;19:1798–808.
Toettcher JE, Apgar JF, Castillo AR, Tidor B, White J. Simulation and verification of electronic and biological systems. Dordrecht: Springer Netherlands; 2011. p. 115–36.
Seglen PO. Preparation of isolated rat liver cells. Methods Cell Biol. 1976;13:29–83.
Narmada BC, Chia S-M, Tucker-Kellogg L, Yu H. HGF regulates the activation of TGF-β1 in rat hepatocytes and hepatic stellate cells. J Cell Physiol. 2012;228:393–401.
We thank members of the Singapore-MIT Alliance for helpful suggestions. HL and BCN have since moved to the Genome Institute of Singapore, 60 Biopolis Street, 138672 Singapore, and LV to the Department of Immunology, Genetics and Pathology, Dag Hammarskjölds väg 20, Rudbeck Laboratory, 751 85 Uppsala, Sweden.
Publication costs were paid by the Singapore Ministry of Health’s National Medical Research Council (NMRC) under its Open Fund Individual Research Grant scheme (OFIRG15nov062) to LTK. This work was also supported by a Singapore-MIT Alliance Computational and Systems Biology Flagship Project grant; by Singapore-MIT Alliance Computational and Systems Biology Programme grants to LTK, JKW, and HY; by funding from the Mechanobiology Institute, the Institute of Bioengineering and Nanotechnology, and grants from Jansson-Cilag Singapore (R-185-000-182-592) to HY.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
About this supplement
This article has been published as part of BMC Systems Biology Volume 11 Supplement 7, 2017: 16th International Conference on Bioinformatics (InCoB 2017): Systems Biology. The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-11-supplement-6.
Ethics approval and consent to participate
The rat cell isolation procedure was approved by the IACUC of the National University of Singapore.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Notes. (PDF 107 kb)
List of equations and parameters used for model construction. (PDF 280 kb)
Parameters settings for different models. (PDF 47 kb)
List of primer sequences for genes probed on quantitative real time PCR. (PDF 26 kb)
About this article
Cite this article
Li, H., Venkatraman, L., Narmada, B.C. et al. Computational analysis reveals the coupling between bistability and the sign of a feedback loop in a TGF-β1 activation model. BMC Syst Biol 11 (Suppl 7), 136 (2017). https://doi.org/10.1186/s12918-017-0508-z
- Positive feedback
- Computational modelling
- Dynamical systems
- Biochemical network
- Bifurcation analysis