Skip to main content

A computational model of mutual antagonism in the mechano-signaling network of RhoA and nitric oxide



RhoA is a master regulator of cytoskeletal contractility, while nitric oxide (NO) is a master regulator of relaxation, e.g., vasodilation. There are multiple forms of cross-talk between the RhoA/ROCK pathway and the eNOS/NO/cGMP pathway, but previous work has not studied their interplay at a systems level. Literature review suggests that the majority of their cross-talk interactions are antagonistic, which motivates us to ask whether the RhoA and NO pathways exhibit mutual antagonism in vitro, and if so, to seek the theoretical implications of their mutual antagonism.


Experiments found mutual antagonism between RhoA and NO in epithelial cells. Since mutual antagonism is a common motif for bistability, we sought to explore through theoretical simulations whether the RhoA-NO network is capable of bistability. Qualitative modeling showed that there are parameters that can cause bistable switching in the RhoA-NO network, and that the robustness of the bistability would be increased by positive feedback between RhoA and mechanical tension.


We conclude that the RhoA-NO bistability is robust enough in silico to warrant the investment of further experimental testing. Tension-dependent bistability has the potential to create sharp concentration gradients, which could contribute to the localization and self-organization of signaling domains during cytoskeletal remodeling and cell migration.


Mutual antagonism between two signaling pathways is a motif of higher-order regulation for building circuits [1,2,3,4,5]. When one competing pathway is active, it turns off the other pathway. By turning off the antagonism from the other pathway, it reinforces and stabilizes its own activity. This double-negative loop is a form of positive feedback that can allow a binary decision to be sharpened and stabilized. Sharpening information signals is crucial because signaling information would otherwise get eroded by environmental noise and by transmission through many sequential steps. Dynamical systems with mutual antagonism have been studied in several contexts, often in conjunction with bistability or oscillation.

Mutual antagonism, when combined with ultrasensitivity, is a common motif for inducing bistability [6]. Bistability literally means that a system has two stable steady-states. The stability of a steady-state means that the system re-converges to the same steady-state after minor perturbations [7]. Meanwhile, the existence of a second stable steady-state means the system can reverse its decision and converge to the opposite state, in response to significant changes from the input stimuli. Bistability has been modeled in many signaling pathways including caspase-dependent apoptosis [8], MAPK signaling [9, 10], direction sensing [11] and establishment of cell polarity [12].

There is growing evidence that bistability in the Rho GTPases network can lead to intracellular spatiotemporal patterns, governing important cellular functions like motility. Mathematical modeling suggests that MLC phosphorylation depends upon bistability between Rho kinase and myosin phosphatase [13]. Bistability in a RhoA-Rac1-PAK network was shown to drive actin cytoskeletal dynamics and cell migration switches in a study combining both modeling and experimental approaches [14]. Bistability in the GTPase network was used to qualitatively model the symmetry breaking and polarization behavior of HeLa cells as observed in microfluidic experiments [15]. In another study, spatially modulated bistability was demonstrated as Rho and Cdc42 zones formed and segregated during Xenopus oocyte wound healing [16]. These studies show that Rho bistability can transform initial signals and graded inputs into spatially precise and temporally co-ordinated events in cytoskeletal regulation. Therefore, we seek to explore the influences and implications of Rho dynamics in other signaling networks, such as the nitric oxide system.

RhoA and nitric oxide (NO) regulate cytoskeletal tension in opposite ways. RhoA signaling activates myosin motors [17] and is a master regulator of cellular contractility. In contrast, NO is a master regulator of relaxation and the 1998 Nobel Prize was awarded for major findings in the role of NO as a signaling molecule in vasorelaxation [18]. Many individual mechanisms of crosstalk have been documented between RhoA and NO. In vascular cell types, many mechanisms exist for the RhoA pathway and its signaling effector ROCK (Rho kinase) to antagonize eNOS and NO. eNOS is a constitutively expressed synthase of NO. Nitric Oxide Synthases (NOS) are categorized as inducible (iNOS) or constitutive (cNOS), including endothelial (eNOS) and neuronal (nNOS) [19]. ROCK decreases the stability of eNOS mRNA and thus inhibits eNOS expression [20, 21]. At the post-translational level, eNOS activity is regulated by Akt-induced phosphorylation at Serine-1177, and Akt activity is regulated by ROCK, which inhibits phosphorylation of Akt [22]. The inhibitory effects of Rho/ROCK on eNOS are well-established, as confirmed with diabetic rats in vivo [23]. Another more rapid mechanism for Rho to antagonize NO production has also been documented in experiments showing that Rho inhibition augmented shear-induced production of NO [24] . The opposite direction of effect, from NO toward RhoA, has been best studied with eNOS-NO signaling to cGMP-PKG in the vascular system [25, 26]. NO activates PKG, which is a serine-threonine kinase that phosphorylates a number of proteins regulating calcium levels and the cytoskeleton. PKG-induced phosphorylation of RhoA at Serine-188 antagonizes RhoA by blocking RhoA translocation to the membrane, which prevents its activation, and by promoting RhoA association with RhoGDI [27,28,29]. In mice lacking the three NOS isoforms, Rho kinase activity was elevated [30]. Most of the reported cross-talk between RhoA and NO has been antagonistic, but positive relationship has also been seen [31,32,33,34].

Prior research on the cross-talk between RhoA and NO has focused on specific elements rather than on systems level interplay. Mutual antagonism would be a logical motif for organizing RhoA and NO signaling in vivo because either extreme could be detrimental to the cell [35] or wasteful of resources [36]. However, it is unclear whether both directions of antagonism co-exist in the same system. The uncertainty arises because NO-induced relaxation has been studied primarily in vascular endothelial cell types, which have high levels of NO, while RhoA-induced contraction has been studied primarily in cell types with strong contractility such as fibroblasts. Because of this uncertainty, the first step in our study is to test experimentally for two key effects in two opposite directions, in a cell type that is neither high-NO nor highly contractile (e.g., an epithelial system). This confirms that the previous forms of antagonistic signaling co-exist, and confirms that they may be studied together in a single model instead of requiring separate endothelial and mesenchymal models. To create NO in our study, we have abstained from using NO donors, which could induce unnaturally high levels of NO, and have instead used Hepatocyte Growth Factor (HGF) to induce endogenous NO production.

After the experimental confirmation of mutual antagonism between RhoA and NO, our next goal is theoretical analysis for the functional implications of the mutual antagonism. We built a model of RhoA-NO signaling and explored whether the system is capable of bistable behavior. The model of RhoA-NO signaling includes 31 reactions of production/consumption (listed in Table S1), which are integrated into 11 governing equations (listed in Table S2). After building the biochemical network model and assessing the potential for bistability, we extend the conceptual scope of the system to include cytokeletal tension, resulting in 12 governing equations. The extended model with tension allows us to assess the influence of mechanical feedback on the qualitative behavior of the system.


Our first objective was to confirm that both directions of RhoA-NO mutual antagonism occur in the same cells. For this verification we selected an epithelial cell line (neither vascular nor highly contractile): MDCK cells stimulated with Hepatocyte Growth Factor (HGF).

NO is required for growth factor-induced phosphorylation of RhoA

To test if NO antagonizes RhoA in epithelial cells, we selected serine-188 phosphorylation of RhoA as a proxy for the ability of NO to antagonize RhoA [25,26,27]. The effect of NO suppression on RhoA phosphorylation was measured in MDCK epithelial cells (Fig. 1A). RhoA is unphosphorylated in resting cells, and HGF triggers phosphorylation of RhoA at Serine 188. MDCK epithelial cells were treated with L-NAME, an inhibitor of cNOS, prior to treatment with HGF. We found that when NO was suppressed by L-NAME pre-treatment, then HGF stimulation was unable to phosphorylate RhoA at Serine 188. An additional side-effect of L-NAME was seen most clearly in cells without growth factor, where L-NAME treatment caused unexpected stimulation of phosphorylation. Such response has been observed in earlier studies where chronic administration of low dose of L-NAME suppressed the negative feedback regulation of NO synthesis by NO [37, 38] itself. For the purposes of confirming previously documented effects, in an epithelial cell type, we can conclude that NO was necessary for RhoA phosphorylation, suggesting that NO-induced antagonism of RhoA occurs significantly in this cell type.

Fig. 1
figure 1

Experimental verification of mutual antagonism between NO and RhoA in epithelial (MDCK) cells. (A and B) NO mediates HGF-induced RhoA-Ser188 phosphorylation. MDCK cells were induced with HGF 1 ng/mL for the specified durations with or without a pre-treatment for 20 min with 50 μM L-NAME, an inhibitor of constitutive NOS. A Western blot analysis with pSer188-RhoA and B Quantification of band intensity. The phospho-RhoA-Ser188 bands in the western blot were quantified using ImageJ and normalized with respect to total RhoA intensity. C, D, E and F The Rho effector, ROCK inhibits p-Akt/p-eNOS signaling. MDCK cells were induced with HGF 1 ng/mL for the specified durations with or without 20 min pre-treatment with 20 uM of ROCK inhibitor Y-27632. C Western blot analysis for shorter durations of HGF treatment. Additional timepoints appear in Additional file 1. Commercial antibodies do not detect the canine isoform of eNOS but phospho-eNOS antibodies are reactive to the canine isoform. Quantification of the bands in the western blot: D phospho-Akt-Ser473 normalized against total Akt, E phospho-Akt-Thr308 normalized against total Akt, and F phospho-eNOS normalized against GAPDH. The white line between lanes denotes that a single gel was cropped to remove other treatment conditions in the middle lanes. (Full gel appears in Additional file 1). **p < 0.005

ROCK inhibits p-Akt/p-eNOS signaling

To test if the antagonistic effect of Rho/ROCK on eNOS is robust and extends to epithelial cells also, we examined whether the Serine 1177 phosphorylation of eNOS is regulated by RhoA/ROCK in MDCK cells. Western blot results (Fig. 1C) showed that treatment with ROCK inhibitor, Y-27632 increased Akt and eNOS phosphorylation. This effect was seen both after initial stimulation with HGF and in unstimulated cells. Additional file 1 Figure S8 (see Addition file 1) shows that long durations of HGF treatment induced an opposite effect, which might be due to a counter-balancing phosphatase [39]. In summary, inhibition of ROCK caused increased Akt and eNOS activation 1 to 2 min after HGF induction, indicating that ROCK inhibited phospho-Akt/phospho-eNOS signaling. We have established that both NO pathway inhibition of RhoA and RhoA pathway inhibition of NO can occur in MDCK cells.

Relationship between mutual antagonism and Bistability

Mutual inhibition between two species is not sufficient by itself to induce bistability, To illustrate this, we simulated a minimal 2-species model of mutual antagonism (double-negative feedback) with mass action kinetics in Fig. 2A. Computational testing of rate parameters (using a multidimensional grid of possible combinations of values) showed that the system was monostable for all combinations (see Additional file 1, Additional file 1 Text S7). A double-negative system can exhibit bistability if one of the inhibitory reactions is ultrasensitive or if there are other forms of ultrasensitivity in the network [6], such as positive feedback. Therefore, the presence of mutual antagonism in the RhoA-NO system is not sufficient to indicate whether the system would be bistable.

Fig. 2
figure 2

Network diagram and simulations of the initial model of HGF-activated RhoA-NO network, including known biochemical effects of mutual antagonism. A Two-node model to demonstrate that mutual antagonism is not sufficient to induce bistability. A and B are two species with zeroth order production and first order degradation rates, ksynthA, kdegA, ksynthB, kdegB for A and B respectively. A and B antagonize each other with rate constants kiA and kiB, according to the mass action equations in Additional file 1 Text S7. B The red arrows represent the inhibitory reactions. Complete specifications of the reactions appear in Additional file 1 Tables S1 and S2. C Model simulation with 100 random initial concentrations. The model was simulated, starting with random levels of all species. The resulting trajectories of the RhoA concentration over time (plotted with a random color for each trajectory) shows convergence to two different steady states. D and E Bifurcation analysis of the initial model. The steady state NO concentration is plotted against variation in the parameter D ki1 (rate constant for ROCK inhibition of Akt phosphorylation) and E ki2 (rate constant for RhoA inhibition of eNOS mRNA). F Robustness analysis of the bistability in the initial model. Single-variable robustness analysis was performed (see Procedure 4 in Additional file 1 Methods). Each of the parameters, from k1 to k15 was varied by the specified percentage, and the system was checked for bistability. White boxes denote the presence of bistability while black boxes denote the absence of bistability. The red line highlights the unperturbed parameter set. Note that bistability was retained in some columns, particularly for parameters involved in the degradation or inactivation of NO, eNOS, and cGMP

Pathway model of RhoA-NO mutual antagonism

To address whether the system is capable of bistability, we performed computational modeling of qualitative network behaviors. A signaling network model was constructed (Fig. 2B) using known biochemical effects and using a previous model of cGMP dynamics [40]. The model simulates changes in concentration over time. We applied a constant input of HGF stimulation to prevent the species from decaying to zero. The reactions in the model are described in Additional file 1 Text S3 (see Additional file 1). Although we simulate RhoA and NO concentrations in the physiological range, the model results should not be used for estimating absolute magnitudes; rather it will be used for generating testable predictions about qualitative behaviors. Even if each reaction could be isolated and quantified, the living cell has unknown effects from scaffolds and competitors.

Bistability in the model

We searched the parameter space of the model for bistability (See Additional file 1, Additional file 1 Text S6, Procedure 2). Next we simulated the model using random initial configurations (See Additional file 1, Additional file 1 Text S6, Procedure 3), and found that the model consistently converged to one of two steady-state concentrations (Fig. 2C). The absolute magnitudes are not predictive, but the crucial observation is that the two states are significantly different from each other and that when in each state, the system remains at the same concentration indefinitely. This result implies that the system is capable of bistability.

To study the nature and extent of the bistable regime, we determined the steady-state(s) of the system for a range of values for the antagonism parameters ki1 and ki2. Figure 2D and E show the relationship between the parameter values and the NO steady-state levels. When the antagonism is very strong or very weak, the system has only one steady-state (monostable), but there exists an intermediate range for which the model has multiple steady-states. This is the characteristic S-shaped curve of bistability, which has two stable steady-states (bistable), plus an unstable steady-state that does not occur in practice. For further information about bifurcation diagrams and the fixed points of dynamical systems, please consult [7].

Next we studied whether bistability was robust to changes in the model parameters. Because this network is large and contains poorly studied reactions, robustness would be crucial for the biological plausibility of bistability. Single parameter perturbations (Fig. 2F) showed that for the majority of parameters, even a small percent deviation caused the system to lose bistability. We conclude that the bistability of this model is brittle, not robust. We are not aware of any systems that lacked robustness in theory but displayed robustness in experiments, so we infer that bistability would be an unlikely behavior for this system alone.

Including mechanical tension in the model

Thus far, we have neglected mechanical tension, but NO and RhoA have well-studied molecular mechanisms for exerting opposite effects on cytoskeletal tension. Indeed, their effects on tension are the reason we have studied both together. We next built an extended model (Fig. 3A) including a simple overall measure of cellular tension as if tension were a species in the reaction network. The NO/cGMP pathway has multiple well-studied mechanisms of decreasing cytoskeletal tension [25, 26, 41], which are represented by Reaction 19 in the extended model. Reaction 15 summarizes the mechanisms by which ROCK increases tension [42]. Finally we include the ability of mechanical tension to activate RhoA [43, 44] in Reaction 16, which is a form of positive feedback with ROCK activity [42]. Another potential source of positive feedback from ROCK to RhoA is mediated by phosphorylation of p190A-RhoGAP [45]. Rho-GEFs and Rho-GAPs are not included in this model, but the results we obtain about tension-mediated feedback are potentially relevant to p190A-RhoGAP feedback as well. Full specification of our reactions appears in Additional file 1 Text S3 (see Additional file 1). In summary, our extended model describes known cross-talk between the biochemical network and a first-order approximation of cellular tension.

Fig. 3
figure 3

Network diagram and simulations of the extended model of HGF-activated RhoA-NO network, including known effects of ‘tension’ on NO and RhoA. A Red arrows represent the inhibitory reactions. Complete specification of the extended model appears in Additional file 1 Tables S1 and S2. B Simulation of the extended model, starting with 100 random levels of all species. Trajectories of the RhoA concentration show convergence to two different steady states. Note faster convergence than in Fig. 2. C and D Bifurcation analysis of the extended model. The steady state NO concentration is plotted against change in parameters C ki1 (rate constant for ROCK inhibition of Akt phosphorylation) and D ki2 (rate constant for RhoA inhibition of eNOS mRNA). Note that bistability occurs over a wider range of ki1 and ki2, compared with Fig. 2. E Robustness analysis (as in Fig. 2) applied to the extended model. White boxes denote the presence of bistability while black boxes denote the absence of bistability. The red line highlights the unperturbed parameter set

We proceeded to test whether the extended model exhibited bistability. Figure 3B shows that the system converged stably to two different steady-states, RhoA-high and RhoA-low, depending on the initial concentrations given to the model. Like the purely biochemical model, the extended model was capable of bistability. Computing bifurcation diagrams for the extended model (Fig. 3C and D) showed that bistability was maintained over a wider range (five-fold increase) of the antagonism parameters. The difference in NO concentration between the two steady-states also increased 5-fold. Single parameter robustness analysis confirmed that the model was capable of bistability for a wider range of the parameters. The system retained bistabliity for 10 to 15% perturbation in 33 out of 36 parameters (Fig. 3E). The bistability was still sensitive to three parameters: Kmcgmp (the rate constant for cGMP in the Hill equation), k12 and kdeg4 (the rate constants for synthesis and degradation of ROCK). These three parameters belong to the “RhoA-side” of the model, suggesting that the effects of tension on RhoA, and generally the feedback loop between RhoA and tension, would be particularly important topics for future study, as they have disproportionate impact on the feasibility of bistability. Robustness is a relative term, and previous simulations of bistability in comparable biochemical signaling networks have shown several cases where 20 to 40% robustness in theory [1, 46, 47] led to cell culture experiments confirming bistability in practice. Because our system retained bistabliity for 30% perturbation in 23 out of 36 parameters, we conclude that bistability is credible enough to begin experimental testing.

Discussion and conclusions

With the use of qualitative modeling, the primary outcome of our work is to assess whether the network topology is consistent with bistability and to find a range of parameter values and concentration levels that could cause bistability (see Additional file 1, Additional file 1 Table S1). Figure 4 shows the flow of logic for this study. An initial ‘tension-free’ model was found capable of bistability for NO levels ranging between 1 to 8 nM (Fig. 2D, E), but the bistability was easily disrupted by parameter perturbations. In the extended model with tension, the bistable regime could range from 1 to 30 nM of NO (Fig. 3C, D). This is within an order of magnitude of the 220 nM NO, previously measured in cerebellar cells [48]. Since the extended model exhibited bistability for moderate variation in most of the parameters, we conclude that the RhoA-NO bistability in silico is sufficiently robust to warrant the investment of testing for bistability in vitro.

Fig. 4
figure 4

Flowchart of prior work and current work. Theoretical work appears in rounded boxes and experimental work appears in rectangles. Current contributions are in black type on white background. Guiding motivations are shaded black. Prior work is shaded grey

Testing for bistability will be challenging because real-time sensors for RhoA activity can provide poor signal-to-noise ratio [49], and tools for dynamic detection of NO have many limitations, for example DAF has non-specificity [50], irreversible binding to NO, and photo-bleaching. Innovations in biosensor technology may enable these limitations to be overcome [51,52,53]. The size of the swing between two bistable states, and the ease of detecting the switch experimentally, would depend on the cell type. Based on our model, the concentration fold-change between high-NO and low-NO states would be highly dependent on parameters such as ki1 and ki2 (Fig. 3C, D) which would depend on cell type; endothelial cells would be expected to exhibit higher levels of NO in the high-NO steady-state, consistent with low values of ki1 and ki2 in Fig. 3C, D, while epithelial cells might be expected to have lower tolerance for NO, consistent with high values of ki1 and ki2.

Several important phenomena have been omitted from our present model. One is the ability of mechanical tension to regulate NO levels. Mechanosensors like the glycocalyx [54] mediate shear stress-induced NO production in endothelial cells. In addition, softening of the actin cortex may induce NO production [55]. If relaxation does have positive feedback, that could further increase the robustness of bistability. Another consideration is the effect of post-translational modification on protein stability. Ser188 phosphorylation of RhoA protects RhoA from degradation [32], so that NO/cGMP/PKG signaling can up-regulate RhoA protein expression, even though it down-regulates RhoA activity. For example, basal NO release from endothelial cells was found to be necessary for maintaining long-term RhoA levels in vascular smooth muscle cells [31]. Other post-translational effects we have omitted include phosphorylation of p190ARhoGAP [45], S-nitrosation of RhoA [56], and tyrosine nitration of RhoA [57] and p190Rho-GAP [58]. In the temporal realm, we have not yet studied dynamic stimuli such as pulsatile tension or oscillatory RhoA [49].

The importance of our work comes from the potential importance of bistability for converting noisy and gradual inputs into distinct transitions and regulatory decisions in the cell. Bistability may have physiological importance for the ability of cells to establish distinct spatial domains of effect, if nearby regions of space exhibit opposite steady states. In the cytoskeleton, bistability might sharpen the boundary between tension and relaxation, which would be useful for cell migration, direction sensing, and establishment of polarity, which are bistable decisions at the highest scale of cellular decision-making [11, 12]. Regional control of NO is reasonable, even though NO has fast diffusion, because NO exhibits narrow spatial targeting in vivo [59]. Bistability might explain how NO can be so narrowly confined, because there could be some regions having a self-reinforcing state of high NO production and other regions blocking NO production. A sharp concentration gradient would then occur between the two regions. A sharp concentration gradient can be interpreted as a spatiotemporal boundary for localizing tension within the cytosol.

In some examples of bistability, a “digital” yes/no decision at the molecular level might appear gradual or graded when viewed at a higher scale [8]. The assembly and disassembly of focal adhesions appears gradual, but it is oppositely regulated by RhoA and NO [60, 61], and is likely to be a discrete decision at the molecular scale [62,63,64]. Regions of greater mechanical tension would have greater positive feedback toward Rho activation, pushing them to a high-Rho state. The high-Rho state would promote stress fiber growth and focal adhesion assembly, consistent with observations of where adhesions form [65]. Regions with less tension would converge toward a high-NO state, and NO would promote disassembly of filamentous actin [66,67,68] and disassembly of adhesions. Because bistability aids in establishing sharp concentration gradients over space, RhoA and NO might have mutually exclusive regions of local up-regulation. Mutually exclusive regions of RhoA and NO could help explain the simultaneous assembly and disassembly of single focal adhesions, observed in heel-toe remodeling of adhesion structure [69].

In endothelial remodeling, shear stress triggers NO production through calcium-dependent pathways not included in our model [70], but the downstream effects of NO might still participate in switch-like behavior with RhoA. Blood flow creates localized regions of high force, inducing adhesions to grow at the upstream interface of the cell [71] with the activation of RhoA [72]. Meanwhile, the remainder of the cell (and its neighboring vascular smooth muscle cells) might relax in response to the shear-induced NO. RhoA-NO bistability might facilitate an orderly juxtapositioning of adhesion and relaxation that benefits cell health, in addition to the known importance of laminar flow versus turbulent flow. There are many possible contexts where tension-dependent bistability between RhoA and NO could contribute to the self-organization of cytoskeletal machinery into discrete spatiotemporal domains.

Materials and methods

Cell culture experiments

Madin-Darby Canine Kidney (MDCK) strain II cells were kindly provided by Walter Hunziker from IMCB-Singapore [73]. Chemicals (L-NAME (L-NG-Nitroarginine Methyl Ester), Y-27362, HGF) were from Sigma. Additional file 1 Texts S4 and S5 (see Additional file 1) describe cell culture methods, sample preparation, Western blotting method and antibody reagents.

Computational procedures and bifurcation analysis

The differential equations were simulated in Matlab (Mathworks, Natick, MA). Procedures for assessing bistability, for parameter search, for random generation of 100 initial conditions, and for robustness analysis appear in Additional file 1 Text S6 (see Additional file 1). Similar methods have been used previously [74, 75]. For bifurcation diagrams, we computed the change in steady-states over a range of parameters, using Govaerts and Kuznetsov’s continuation software Matcont [76], and Matcont parameters were chosen to ensure the bistable region was evaluated.

Statistical analysis

Experiments were performed in triplicate. Error bars display the standard errors of the mean (±SEM). Statistical significance was determined using Student’s t-test.

Availability of data and materials

Not applicable.



Nitric Oxide


Rho Kinase


Inducible Nitric Oxide Synthase


Endothelial Nitric Oxide Synthase


Constitutional Nitric Oxide Synthase


Hepatocyte Growth Factor


Madin-Darby Canine Kidney cells


L-N-Nitro arginine methyl ester


Rho GTPase-activating protein


Guanine nucleotide exchange factor


  1. Venkatraman L, Chia SM, Narmada BC, White JK, Bhowmick SS, Forbes Dewey C Jr, et al. Plasmin triggers a switch-like decrease in thrombospondin-dependent activation of TGF-beta1. Biophys J. 2012;103(5):1060–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Kim SY, Ferrell JE Jr. Substrate competition as a source of ultrasensitivity in the inactivation of Wee1. Cell. 2007;128(6):1133–45.

    Article  CAS  PubMed  Google Scholar 

  3. Bartos M, Vida I, Jonas P. Synaptic mechanisms of synchronized gamma oscillations in inhibitory interneuron networks. Nat Rev Neurosci. 2007;8(1):45–56.

    Article  CAS  PubMed  Google Scholar 

  4. Chickarmane V, Enver T, Peterson C. Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility. PLoS Comput Biol. 2009;5(1):e1000268.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Kurosawa G, Aihara K, Iwasa Y. A model for the circadian rhythm of cyanobacteria that maintains oscillation without gene expression. Biophys J. 2006;91(6):2015–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Zhang Q, Bhattacharya S, Andersen ME. Ultrasensitive response motifs: basic amplifiers in molecular signalling networks. Open Biol. 2013;3(4):130031.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Strogatz SH. Nonlinear dynamics and chaos : with applications to physics, biology, chemistry, and engineering. Cambridge: Westview Press; 2000.

    Google Scholar 

  8. Eissing T, Conzelmann H, Gilles ED, Allgower F, Bullinger E, Scheurich P. Bistability analyses of a caspase activation model for receptor-induced apoptosis. J Biol Chem. 2004;279(35):36892–7.

    Article  CAS  PubMed  Google Scholar 

  9. Ferrell JE Jr, Machleder EM. The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science. 1998;280(5365):895–8.

    Article  CAS  PubMed  Google Scholar 

  10. Bagowski CP, Ferrell JE Jr. Bistability in the JNK cascade. Curr Biol. 2001;11(15):1176–82.

    Article  CAS  PubMed  Google Scholar 

  11. Beta C, Amselem G, Bodenschatz E. A bistable mechanism for directional sensing. New J Phys. 2008;10(8):083015.

    Article  CAS  Google Scholar 

  12. Semplice M, Veglio A, Naldi G, Serini G, Gamba A. A bistable model of cell polarity. PLoS One. 2012;7(2):e30977.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kaneko-Kawano T, Takasu F, Naoki H, Sakumura Y, Ishii S, Ueba T, et al. Dynamic regulation of myosin light chain phosphorylation by rho-kinase. PLoS One. 2012;7(6):e39269.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Byrne KM, Monsefi N, Dawson JC, Degasperi A, Bukowski-Wills JC, Volinsky N, et al. Bistability in the Rac1, PAK, and RhoA signaling network drives actin cytoskeleton dynamics and cell motility switches. Cell systems. 2016;2(1):38–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Holmes WR, Lin B, Levchenko A, Edelstein-Keshet L. Modelling cell polarization driven by synthetic spatially graded Rac activation. PLoS Comput Biol. 2012;8(6):e1002366.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Simon CM, Vaughan EM, Bement WM, Edelstein-Keshet L. Pattern formation of rho GTPases in single cell wound healing. Mol Biol Cell. 2013;24(3):421–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Amano M, Ito M, Kimura K, Fukata Y, Chihara K, Nakano T, et al. Phosphorylation and activation of myosin by rho-associated kinase (rho-kinase). J Biol Chem. 1996;271(34):20246–9.

    Article  CAS  PubMed  Google Scholar 

  18. SoRelle R. Nobel prize awarded to scientists for nitric oxide discoveries. Circulation. 1998;98(22):2365–6.

    Article  CAS  PubMed  Google Scholar 

  19. Alderton WK, Cooper CE, Knowles RG. Nitric oxide synthases: structure, function and inhibition. Biochem J. 2001;357(Pt 3):593–615.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Eto M, Barandier C, Rathgeb L, Kozai T, Joch H, Yang Z, et al. Thrombin suppresses endothelial nitric oxide synthase and upregulates endothelin-converting enzyme-1 expression by distinct pathways: role of rho/ROCK and mitogen-activated protein kinase. Circ Res. 2001;89(7):583–90.

    Article  CAS  PubMed  Google Scholar 

  21. Laufs U, Liao JK. Post-transcriptional regulation of endothelial nitric oxide synthase mRNA stability by rho GTPase. J Biol Chem. 1998;273(37):24266–71.

    Article  CAS  PubMed  Google Scholar 

  22. Ming XF, Viswambharan H, Barandier C, Ruffieux J, Kaibuchi K, Rusconi S, et al. Rho GTPase/rho kinase negatively regulates endothelial nitric oxide synthase phosphorylation through the inhibition of protein kinase B/Akt in human endothelial cells. Mol Cell Biol. 2002;22(24):8467–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bivalacqua TJ, Champion HC, Usta MF, Cellek S, Chitaley K, Webb RC, et al. RhoA/rho-kinase suppresses endothelial nitric oxide synthase in the penis: a mechanism for diabetes-associated erectile dysfunction. Proc Natl Acad Sci U S A. 2004;101(24):9121–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kolluru GK, Majumder S, Chatterjee S. Rho-kinase as a therapeutic target in vascular diseases: striking nitric oxide signaling. Nitric Oxide. 2014;43:45–54.

    Article  CAS  PubMed  Google Scholar 

  25. Lee MR, Li L, Kitazawa T. Cyclic GMP causes Ca2+ desensitization in vascular smooth muscle by activating the myosin light chain phosphatase. J Biol Chem. 1997;272(8):5063–8.

    Article  CAS  PubMed  Google Scholar 

  26. Carvajal JA, Germain AM, Huidobro-Toro JP, Weiner CP. Molecular mechanism of cGMP-mediated smooth muscle relaxation. J Cell Physiol. 2000;184(3):409–20.<409::AID-JCP16>3.0.CO;2-K.

  27. Chitaley K, Webb RC. Nitric oxide induces dilation of rat aorta via inhibition of rho-kinase signaling. Hypertension. 2002;39(2 Pt 2):438–42.

    Article  CAS  PubMed  Google Scholar 

  28. Sawada N, Itoh H, Yamashita J, Doi K, Inoue M, Masatsugu K, et al. cGMP-dependent protein kinase phosphorylates and inactivates RhoA. Biochem Biophys Res Commun. 2001;280(3):798–805.

    Article  CAS  PubMed  Google Scholar 

  29. Ellerbroek SM, Wennerberg K, Burridge K. Serine phosphorylation negatively regulates RhoA in vivo. J Biol Chem. 2003;278(21):19023–31.

    Article  CAS  PubMed  Google Scholar 

  30. Noda K, Godo S, Saito H, Tsutsui M, Shimokawa H. Opposing roles of nitric oxide and rho-kinase in lipid metabolism in mice. Tohoku J Exp Med. 2015;235(3):171–83.

    Article  CAS  PubMed  Google Scholar 

  31. Sauzeau V, Rolli-Derkinderen M, Marionneau C, Loirand G, Pacaud P. RhoA expression is controlled by nitric oxide through cGMP-dependent protein kinase activation. J Biol Chem. 2003;278(11):9472–80.

    Article  CAS  PubMed  Google Scholar 

  32. Rolli-Derkinderen M, Sauzeau V, Boyer L, Lemichez E, Baron C, Henrion D, et al. Phosphorylation of serine 188 protects RhoA from ubiquitin/proteasome-mediated degradation in vascular smooth muscle cells. Circ Res. 2005;96(11):1152–60.

    Article  CAS  PubMed  Google Scholar 

  33. Zhan R, Yang S, He W, Wang F, Tan J, Zhou J, et al. Nitric oxide enhances keratinocyte cell migration by regulating rho GTPase via cGMP-PKG signalling. PLoS One. 2015;10(3):e0121551.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Zhan R, He W, Wang F, Yao Z, Tan J, Xu R, et al. Nitric oxide promotes epidermal stem cell migration via cGMP-rho GTPase signalling. Sci Rep. 2016;6(1):30687.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Ridnour LA, Thomas DD, Mancardi D, Espey MG, Miranda KM, Paolocci N, et al. The chemistry of nitrosative stress induced by nitric oxide and reactive nitrogen oxide species. Putting perspective on stressful biological situations. Biol Chem. 2004;385(1):1–10.

    Article  CAS  PubMed  Google Scholar 

  36. Bernstein BW, Bamburg JR. Actin-ATP hydrolysis is a major energy drain for neurons. J Neurosci. 2003;23(1):1–6.

    Article  CAS  Google Scholar 

  37. Bernatova I, Kopincova J, Puzserova A, Janega P, Babal P. Chronic low-dose L-NAME treatment increases nitric oxide production and vasorelaxation in normotensive rats. Physiol Res / Academia Scientiarum Bohemoslovaca. 2007;56(Suppl 2):S17–24.

    CAS  Google Scholar 

  38. Kopincova J, Puzserova A, Bernatova I. L-NAME in the cardiovascular system - nitric oxide synthase activator? Pharmacol Rep. 2012;64(3):511–20.

    Article  CAS  PubMed  Google Scholar 

  39. Greif DM, Kou R, Michel T. Site-specific dephosphorylation of endothelial nitric oxide synthase by protein phosphatase 2A: evidence for crosstalk between phosphorylation sites. Biochemistry. 2002;41(52):15845–53.

    Article  CAS  PubMed  Google Scholar 

  40. Yang J, Clark JW, Bryan RM, Robertson CS. Mathematical modeling of the nitric oxide/cGMP pathway in the vascular smooth muscle cell. Am J Physiol Heart Circ Physiol. 2005;289(2):H886–97.

    Article  CAS  PubMed  Google Scholar 

  41. Yang B, Gwozdz T, Dutko-Gwozdz J, Bolotina VM. Orai1 and Ca2+−independent phospholipase A2 are required for store-operated Icat-SOC current, Ca2+ entry, and proliferation of primary vascular smooth muscle cells. Am J Physiol Cell Physiol. 2012;302(5):C748–56.

    Article  CAS  PubMed  Google Scholar 

  42. Bhadriraju K, Yang M, Alom Ruiz S, Pirone D, Tan J, Chen CS. Activation of ROCK by RhoA is regulated by cell adhesion, shape, and cytoskeletal tension. Exp Cell Res. 2007;313(16):3616–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zhao XH, Laschinger C, Arora P, Szaszi K, Kapus A, McCulloch CA. Force activates smooth muscle alpha-actin promoter activity through the rho signaling pathway. J Cell Sci. 2007;120(Pt 10):1801–9.

    Article  CAS  PubMed  Google Scholar 

  44. Lessey EC, Guilluy C, Burridge K. From mechanical force to RhoA activation. Biochemistry. 2012;51(38):7420–32.

    Article  CAS  PubMed  Google Scholar 

  45. Mori K, Amano M, Takefuji M, Kato K, Morita Y, Nishioka T, et al. Rho-kinase contributes to sustained RhoA activation through phosphorylation of p190A RhoGAP. J Biol Chem. 2009;284(8):5067–76.

    Article  CAS  PubMed  Google Scholar 

  46. Venkatraman L, Li H, Dewey CF Jr, White JK, Bhowmick SS, Yu H, et al. Steady states and dynamics of urokinase-mediated plasmin activation in silico and in vitro. Biophys J. 2011;101(8):1825–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Chen C, Cui J, Lu H, Wang R, Zhang S, Shen P. Modeling of the role of a Bax-activation switch in the mitochondrial apoptosis decision. Biophys J. 2007;92(12):4304–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Bellamy TC, Wood J, Goodwin DA, Garthwaite J. Rapid desensitization of the nitric oxide receptor, soluble guanylyl cyclase, underlies diversity of cellular cGMP responses. Proc Natl Acad Sci U S A. 2000;97(6):2928–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Machacek M, Hodgson L, Welch C, Elliott H, Pertz O, Nalbant P, et al. Coordination of rho GTPase activities during cell protrusion. Nature. 2009;461(7260):99–103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Zhang X, Kim WS, Hatcher N, Potgieter K, Moroz LL, Gillette R, et al. Interfering with nitric oxide measurements. 4,5-diaminofluorescein reacts with dehydroascorbic acid and ascorbic acid. J Biol Chem. 2002;277(50):48472–8.

    Article  CAS  PubMed  Google Scholar 

  51. Eroglu E, Gottschalk B, Charoensin S, Blass S, Bischof H, Rost R, et al. Development of novel FP-based probes for live-cell imaging of nitric oxide dynamics. Nat Commun. 2016;7(1):10623.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Shcherbakova DM, Cox Cammer N, Huisman TM, Verkhusha VV, Hodgson L. Direct multiplex imaging and optogenetics of rho GTPases enabled by near-infrared FRET. Nat Chem Biol. 2018;14(6):591–600.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Marston DJ, Vilela M, Huh J, Ren J, Azoitei ML, Glekas G, et al. Multiplexed GTPase and GEF biosensor imaging enables network connectivity analysis. Nat Chem Biol. 2020;16(8):826–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Yao Y, Rabodzey A, Dewey CF Jr. Glycocalyx modulates the motility and proliferative response of vascular endothelium to fluid shear stress. Am J Physiol Heart Circ Physiol. 2007;293(2):H1023–30.

    Article  CAS  PubMed  Google Scholar 

  55. Fels J, Jeggle P, Kusche-Vihrog K, Oberleithner H. Cortical actin nanodynamics determines nitric oxide release in vascular endothelium. PLoS One. 2012;7(7):e41520.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Zuckerbraun BS, Stoyanovsky DA, Sengupta R, Shapiro RA, Ozanich BA, Rao J, et al. Nitric oxide-induced inhibition of smooth muscle cell proliferation involves S-nitrosation and inactivation of RhoA. Am J Physiol Cell Physiol. 2007;292(2):C824–31.

    Article  CAS  PubMed  Google Scholar 

  57. Rafikov R, Dimitropoulou C, Aggarwal S, Kangath A, Gross C, Pardo D, et al. Lipopolysaccharide-induced lung injury involves the nitration-mediated activation of RhoA. J Biol Chem. 2014;289(8):4710–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Siddiqui MR, Komarova YA, Vogel SM, Gao X, Bonini MG, Rajasingh J, et al. Caveolin-1-eNOS signaling promotes p190RhoGAP-A nitration and endothelial permeability. J Cell Biol. 2011;193(5):841–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Gladwin MT, Kim-Shapiro DB. Vascular biology: nitric oxide caught in traffic. Nature. 2012;491(7424):344–5.

    Article  CAS  PubMed  Google Scholar 

  60. Goligorsky MS, Abedi H, Noiri E, Takhtajan A, Lense S, Romanov V, et al. Nitric oxide modulation of focal adhesions in endothelial cells. Am J Phys. 1999;276(6 Pt 1):C1271–81.

    Article  CAS  Google Scholar 

  61. Ridley AJ, Hall A. The small GTP-binding protein rho regulates the assembly of focal adhesions and actin stress fibers in response to growth factors. Cell. 1992;70(3):389–99.

    Article  CAS  PubMed  Google Scholar 

  62. Changede R, Cai H, Wind SJ, Sheetz MP. Integrin nanoclusters can bridge thin matrix fibres to form cell-matrix adhesions. Nat Mater. 2019;18(12):1366–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Kalappurakkal JM, Anilkumar AA, Patra C, van Zanten TS, Sheetz MP, Mayor S. Integrin Mechano-chemical Signaling Generates Plasma Membrane Nanodomains that Promote Cell Spreading. Cell. 2019;177(7):1738–56.e23.

    Article  CAS  Google Scholar 

  64. Yang B, Wolfenson H, Chung VY, Nakazawa N, Liu S, Hu J, et al. Stopping transformed cancer cell growth by rigidity sensing. Nat Mater. 2020;19(2):239–50.

    Article  CAS  PubMed  Google Scholar 

  65. Shemesh T, Verkhovsky AB, Svitkina TM, Bershadsky AD, Kozlov MM. Role of focal adhesions and mechanical stresses in the formation and progression of the lamellipodium-lamellum interface [corrected]. Biophys J. 2009;97(5):1254–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Gorodeski GI. NO increases permeability of cultured human cervical epithelia by cGMP-mediated increase in G-actin. Am J Physiol Cell Physiol. 2000;278(5):C942–52.

    Article  CAS  PubMed  Google Scholar 

  67. Lu J, Katano T, Okuda-Ashitaka E, Oishi Y, Urade Y, Ito S. Involvement of S-nitrosylation of actin in inhibition of neurotransmitter release by nitric oxide. Mol Pain. 2009;5:58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Clancy R, Leszczynska J, Amin A, Levartovsky D, Abramson SB. Nitric oxide stimulates ADP ribosylation of actin in association with the inhibition of actin polymerization in human neutrophils. J Leukoc Biol. 1995;58(2):196–202.

    Article  CAS  PubMed  Google Scholar 

  69. Wolfenson H, Henis YI, Geiger B, Bershadsky AD. The heel and toe of the cell’s foot: a multifaceted approach for understanding the structure and dynamics of focal adhesions. Cell Motil Cytoskeleton. 2009;66(11):1017–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Koo A, Nordsletten D, Umeton R, Yankama B, Ayyadurai S, Garcia-Cardena G, et al. In silico modeling of shear-stress-induced nitric oxide production in endothelial cells through systems biology. Biophys J. 2013;104(10):2295–306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Girard PR, Nerem RM. Shear stress modulates endothelial cell morphology and F-actin organization through the regulation of focal adhesion-associated proteins. J Cell Physiol. 1995;163(1):179–93.

    Article  CAS  PubMed  Google Scholar 

  72. Wojciak-Stothard B, Ridley AJ. Shear stress-induced endothelial cell polarization is mediated by rho and Rac but not Cdc42 or PI 3-kinases. J Cell Biol. 2003;161(2):429–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Phua DC, Humbert PO, Hunziker W. Vimentin regulates scribble activity by protecting it from proteasomal degradation. Mol Biol Cell. 2009;20(12):2841–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Ferrell JE, Tsai TY-C, Yang Q. Modeling the cell cycle: why do certain circuits oscillate? Cell. 2011;144(6):874–85.

    Article  CAS  PubMed  Google Scholar 

  75. Gérard C, Gonze D, Goldbeter A. Revisiting a skeleton model for the mammalian cell cycle: from bistability to Cdk oscillations and cellular heterogeneity. J Theor Biol. 2019;461:276–90.

    Article  CAS  PubMed  Google Scholar 

  76. Dhooge A, Govaerts W, Kuznetsov YA. MATCONT: a MATLAB package for numerical bifurcation analysis of ODEs. ACM Trans Math Softw. 2003;29(2):141–64.

    Article  Google Scholar 

Download references


We are very grateful to Prof. Alpha Yap for helpful scientific discussions, to Dr. Huipeng Li for help with computational methods, and to Dr. Lakshmi Venkatraman and Dr. Narendra Suhas Jagannathan for advice.

About this supplement

This article has been published as part of BMC Molecular and Cell Biology Volume 22 Supplement 1, 2021: 19th International Conference on Bioinformatics 2020 (InCoB2020): molecular and cell biology. The full contents of the supplement are available online at


This research is supported by the Singapore Ministry of Health’s National Medical Research Council under its OFLCG scheme, grant NMRC/OFLCG/002/2018 to LTK, and by the Singapore Ministry of Education Academic Research Fund Tier 2 grant MOE2019-T2–1-138 to LTK. Publication costs were funded by MOE2019-T2–1-138. Funders had no role in the design of the study. Any opinions, findings and conclusions or recommendations expressed in this material are those of the author(s) and do not reflect the views of the funding agencies.

Author information

Authors and Affiliations



AS carried out all the experiments in vitro and in silico. LTK conceived of the project. LBC supervised the in vitro experiments. The manuscript was written by AS and LTK. All authors reviewed and commented on the paper, and all authors approved of the final manuscript.

Corresponding author

Correspondence to Lisa Tucker-Kellogg.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary materials for “A computational model of mutual antagonism in the mechano-signaling network of RhoA and Nitric Oxide” including Table S1 (Reactions, parameters and variables for the initial and the extended models), Table S2 (Ordinary Differential Equations and the corresponding reactions from Table S1), Table S3 (Initial concentrations for extreme initialization of ‘NO-high’ & ‘RhoA-high’ states), Table S4 (Ordinary Differential Equations (ODEs) and parameter values), Texts S3, S4, S5, S6, S7 (Supplementary Materials and Methods), Figure S8 (Raw, unedited gel images of the western blot in Fig. 1A), Figure S9 (Raw, unedited gel images of the western blot in Fig. 1C) and Figure S10. (ROCK inhibits Akt/eNOS phosphorylation at shorter durations of HGF treatment while long durations of HGF treatment induced an opposite effect.)

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Surendran, A., Forbes Dewey, C., Low, B.C. et al. A computational model of mutual antagonism in the mechano-signaling network of RhoA and nitric oxide. BMC Mol and Cell Biol 22 (Suppl 1), 47 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: