An Integrated Mathematical Model of Thrombin-, Histamine-and VEGF-Mediated Signalling in Endothelial Permeability
© Wei et al; licensee BioMed Central Ltd. 2011
Received: 16 February 2011
Accepted: 15 July 2011
Published: 15 July 2011
Endothelial permeability is involved in injury, inflammation, diabetes and cancer. It is partly regulated by the thrombin-, histamine-, and VEGF-mediated myosin-light-chain (MLC) activation pathways. While these pathways have been investigated, questions such as temporal effects and the dynamics of multi-mediator regulation remain to be fully studied. Mathematical modeling of these pathways facilitates such studies. Based on the published ordinary differential equation models of the pathway components, we developed an integrated model of thrombin-, histamine-, and VEGF-mediated MLC activation pathways.
Our model was validated against experimental data for calcium release and thrombin-, histamine-, and VEGF-mediated MLC activation. The simulated effects of PAR-1, Rho GTPase, ROCK, VEGF and VEGFR2 over-expression on MLC activation, and the collective modulation by thrombin and histamine are consistent with experimental findings. Our model was used to predict enhanced MLC activation by CPI-17 over-expression and by synergistic action of thrombin and VEGF at low mediator levels. These may have impact in endothelial permeability and metastasis in cancer patients with blood coagulation.
Our model was validated against a number of experimental findings and the observed synergistic effects of low concentrations of thrombin and histamine in mediating the activation of MLC. It can be used to predict the effects of altered pathway components, collective actions of multiple mediators and the potential impact to various diseases. Similar to the published models of other pathways, our model can potentially be used to identify important disease genes through sensitivity analysis of signalling components.
The endothelium is a semi-permeable barrier that regulates the flux of liquid and solutes between the blood and surrounding tissues. Endothelial permeability increases paracellular leakage of plasma fluid and proteins to surrounding tissues, and intravasation of tissue-released contents to the blood in the development of multiple diseases related to injury (such as edema, trauma, ischaemia-reperfusion injury, respiratory distress syndrome, and thrombosis), inflammation (such as atherosclerosis and sepsis), diabetes, and cancer [1–4]. The level of endothelial permeability is regulated individually or in combination by multiple mediators, particularly thrombin, histamine, and vascular endothelial growth factor (VEGF), under various disease conditions .
The proinflammatory and vasoactive factors thrombin, generated in thrombosis and inflammatory diseases, and histamine, produced in acute inflammatory responses to trauma, burns, allergy, and infection, induce transient endothelial permeability to link inflammation, tissue injury and vascular leakage to cellular responses and symptoms [5–7]. VEGF, released in diabetic retinopathy, I-R injury, vasculogenesis, angiogenesis, and tumor development and metastasis, causes endothelial permeability to enable extravasation of fluids and solutes and intravasation of tumor cells [8–10]. These three key mediators stimulate their respective receptors on endothelial cells to individually and collectively activate Ca2+, Rho GTPase/ROCK, and Myosin light chain kinase (MLCK) signalling pathways that subsequently activate myosin light chain (MLC) to induce cytoskeleton contraction in endothelial cells and dissociation of cell-cell junctions, resulting in endothelial hyper-permeability [4, 11].
Significant progress has been made in understanding the molecular mechanism and dynamics of the relevant signalling events [4, 7, 9, 11, 12] and the roles of different regulators [13, 14]. Nonetheless, some puzzles still remain to be elucidated. For instance, it is unclear what contributes to the different temporal effects and permeability recovery rates by histamine, thrombin, and VEGF mediated signalling, given that they share similar signalling cascades in triggering endothelial permeability. Another question is how multiple mediators under certain complicated inflammatory conditions collectively reduce the effectiveness of antagonizing agents directed at individual mediator-mediated signalling .
As part of the efforts for solving these puzzles and for quantitative and mechanistic study of the relevant signalling events, mathematical models have been developed for analyzing the relevant signalling and regulation processes [15–20]. In particular, ordinary differential equation (ODE) based mathematical models of thrombin, Ca2+-calmodulin (CaM), and Rho activation have been developed for investigating the thrombin-mediated activation of MLC , and Ca2+-CaM, MLCK and Myosin Light chain phosphatase (MYCP) on MLC activation [15, 16, 21]. To enable more comprehensive analysis of signalling in endothelial permeability, there is a need to develop an expanded ODE model that covers the signalling mediated by multiple mediators, particularly thrombin, histamine and VEGF.
Thrombin-, Histamine-and VEGF-Mediated Signalling Cascades in Endothelial Permeability
Thrombin mediated GPCR activation
Thrombin regulates endothelial permeability, inflammation and other events via activation of thrombin receptors such as PAR-1 by proteolytically cleaving the N-terminus of these receptors . PAR-1 is the main receptor in the regulation of endothelial permeability (Additional File 1, Figure S1). It interacts with Gq to increase the concentration of Ca2+ and activate protein kinase C, inositol 1, 4, 5-triphosphate and diacylglycerol . It is also linked to G12/13  to activate the small G-protein Rho .
VEGF mediated ERK activation
VEGF regulates angiogenesis, cancer and microvascular permeability under various physiological and pathological conditions by activating transmembrane tyrosine kinase receptors VEGFR-2 and Flt-1, which promotes mitogenic, chemotactic, and prosurvival signalling and activates phospholipase C (PLC), intracellular Ca2+, and various protein kinase C (PKC) isoforms. In particular, VEGF activates ERK-1/2 via the Raf-MEK-ERK cascade . Accumulative evidences suggest that ERK-MLCK-meditated cytoskeletal responses contribute to VEGF-elicited microvascular hyperpermeability. For instance, MLCK has been found to contain multiple MAP kinase consensus phosphorylation sites (P-x-S-P) that can be directly phosphorylated by MAP kinase , which is supported by additional experimental evidence indicating MLCK as a substrate for MAP kinase .
Histamine, VEGF mediated NO activation
Cytosolic Ca2+ elevation is a typical initial response of endothelial cells to hormonal and chemical signal and to changes in physical parameters, and many endothelial functions are dependent on changes in Ca2+ concentration . For instance, the activity of endothelial nitric oxide synthase (eNOS) in producing nitric oxide in endothelial cells absolutely requires CaM  and it appears to also require Ca2+ to sustain elevated level of activity .
Nitric oxide plays a critical role in the endothelial cell proliferation, migration, and tube formation, as well as increased vascular permeability, hypotension, and angiogenesis in vivo [35, 36]. VEGF- and histamine-induced microvascular hyperpermeability are both mediated by a signalling cascade triggered by receptor binding and transduced by a serial activation of intracellular enzymes, including PLC, eNOS, soluble guanylate cyclase (sGC), and protein kinase G (PKG). Subsequently, the VEGF-activated NO-PKG pathway was linked to ERK1/2-mediated proliferation of cultured endothelial cells via phosphorylation and activation of the upstream p42/44 MAPK cascade component RAF by PKG .
Thrombin, VEGF and Histamine mediated MLC activation
MLC of myosin II plays a critical role in controlling actomyosin contractility in both smooth muscle and nonmuscle cells [38–40]. MLC phosphorylation is regulated by the balance of two enzymatic activities, i.e., MLCK and myosin phosphatase (MYCP). MLCK is regulated by Ca2+/calmodulin and is believed to be a major kinase in both smooth muscle and nonmuscle cells. In addition, Rho-kinase (ROCK) can directly phosphorylate MLC in vitro . MYCP is a holoenzyme composed of three subunits: a catalytic subunit of 38 kDa that was identified as protein phosphatase 1 (PP1) catalytic subunit δ-isoform (PP1Cδ)  and two noncatalytic subunits of 21 and 110-130 kDa [43, 44]. The larger one, called myosin phosphatase targeting subunit 1 (MYPT1), binds to the catalytic subunit and targets it to MLC, providing substrate specificity . ROCK and PKC have been proposed to mediate the inhibition of smooth muscle MYCP, leading to increased MLC phosphorylation in response to various agonists. Phosphorylation of the MYPT1 regulatory site (Thr695 in chicken MYPT1) by ROCK induces inhibition of MYCP activity . Some experimental findings suggest that CPI-17, a soluble globular protein, is involved in PKC-dependent inhibition of MYCP and it has thus been considered as a specific inhibitor for MYCP .
Detail description of signaling cascades used in this model was provided the Additional File 3.
Results and Discussion
Model validation with experimental studies of the regulation of MLC activation, calcium release, and Rho activation by thrombin
Elevation in cytosolic Ca2+ concentration () is a common initial response of endothelial cells to various changes such as the exposure to hormonal and inflammatory stimuli and variation of physical conditions . Jeng et al.  have shown that the binding of thrombin and PAR-1 induces rapid calcium mobilization and increase of [Ca2+]i, with [Ca2+]i peaking at 30-40 s followed by a rapid drop. The simulated calcium release profile in Figure 2 (Right) exhibits a peak concentration at 38 s followed by a rapid decay, consistent with Jeng et al's experiment results. The increased intracellular Ca2+ influx is expected to enhance the binding of Ca2+ to CaM, which subsequently activates MLCK to phosphorylate the MLC of myosin II. To evaluate which signalling event is primarily responsible for the large transient increase in the level of MLC activation (the first peak at ~2.5 min in the left Figure 2), we systematically varied the strength of protein-protein interactions upstream of MLC. As shown in Additional File 1, Figure S5, the first peak disappears when the Ca2+-dependent MLC activation (Reaction 73-86) was switched off, while that peak remains largely intact when the ROCK-dependent MLC activation and CPI-17-MYPT interactions were switched off (Reactions 57-58, 63-70, 99-102), Therefore, our analysis suggests that this Ca2+-dependent mechanism was primarily responsible for the large transient increase of the levels of MLC activation.
Thrombin induces a prolonged increase of endothelial permeability lasting for 1-1.5 hr. This prolonged elevated permeability is attributed to the activation of the small Rho GTPase and Rho kinase [54, 55]. It has been found that Rho GTPase activation can be observed after 2 min and the elevated activation is maintained up to 60 min after thrombin stimulation, and the time course of Rho GTPase activation correlates well with the time course of MLC activation increase by experiment observation . Our simulation result presented in Additional File 1, Figure S6 is consistent with this observation, which shows that the simulated Rho GTPase activation was maintained for 60 min. Rho GTPase activation induces MLC activation via both direct and indirect routes. Rho GTPase and ROCK directly activate MLC to subsequently induce the contraction of the non-muscle cell systems [41, 56]. In the indirect route, ROCK inhibits myosin phosphatase activity by phosphorylating the myosin binding subunit (MBS) of myosin phosphatase , which increases the activation level of MLC, actomyosin interaction, stress fiber formation, and subsequent endothelial permeability. We studied whether these direct and indirect Rho GTPase -dependent mechanisms are primarily responsible for the sustained levels of MLC activation (the main peak at ~30 min in the left Figure 2) by systematically varying the protein-protein interactions upstream of MLC. As shown in Additional File 1, Figure S7, this peak remains largely intact when the Ca2+ -dependent MLC activation and P-CPI-17-MYPT interaction (Reaction 57-58, 71-76) were switched off, but disappeared when the ROCK-dependent MLC activation (Reactions 63-70, 99-102) were switched off. Therefore, our analysis suggests that both the direct and indirect Rho GTPase -dependent mechanisms play an important role for the sustained levels of MLC activation.
Model validation with experimental studies of MLC activation and ERK activation by VEGF
Model validation with experimental studies of MLC activation by histamine
The model was further validated by determining whether the simulation results are consistent with experimentally observed regulation of MLC activation by the third individual mediator histamine. This histamine-mediated process was simulated by using our model with thrombin and VEGF switched off by setting their initial concentrations to zero values. The simulation results in Additional File 1, Figure S9 indicated that histamine causes a transient increase of MLC activation that peaked at 2.5 min, which is consistent with the experimental finding that histamine induces a transient endothelial permeability peaked at 2-5 min . Further investigation showed that this peak is primarily induced by Ca2+-dependent mechanism and the contribution from the NO-dependent ERK activation path is very small, as shown in Additional File 1, Figure S10 by switching off each individual path. Moreover, the contribution from the NO-dependent ERK activation path is much weaker compared with Ras-dependent ERK activation and MLC activation by the individual mediator VEGF.
Comparison of the simulated thrombin-mediated IP3 and Ca2+ release with that of an existing model
Simulation of the effects of thrombin receptor PAR-1 over-expression on thrombin-mediated MLC activation
Simulation of the effects of Rho GTPase and ROCK over-expression on thrombin-mediated MLC activation
Simulation of effects of VEGF and VEGFR2 over-expression on VEGF-mediated MLC activation
VEGFR2 is recognized as the principal mediator of physiological and pathological effects of VEGF on endothelial cells, which include proliferation, migration, survival, and permeability . The expression of VEGF and VEGFR2 in endothelial cells has been found to be elevated in oxidative stress , type 1 leprosy reaction , and during diabetes to induce microvascular complications, especially diabetic retinopathy . Over-expression of VEGF and VEGFR2 has been shown to correlate with increased risk of metastatic disease and overall poor prognosis in different carcinomas . Apart from their primary functions in angiogenesis, the roles of VEGF and VEGFR2 in metastasis likely involve the regulation of endothelial permeability to facilitate cell transmigration and plasma accumulation in the matrix in support of new vessel formation . The effects of VEGF and VEGFR2 over-expression on VEGF-mediated MLC activation were quantitatively evaluated by using our model with thrombin and histamine switched off.
Simulation of synergistic activation of MLC by thrombin and histamine
Comparison of the areas with respect to different time ranges in Figure 8
MLC activation curve
Area under MLC activation curve with respect to different time range (time unit:min)
Curve 1:Histamine + Thrombin induced activation
Curve 2:Histamine-mediated activation
Curve 3:Thrombin-mediated activation
Simple sum of curve 2 and 3
Area difference between curve 1 and simple sum of curve 2 and 3
As shown in Figure 8, the synergistic effect at low concentrations of thrombin and histamine only occur during the time range from 10 min to 50 min. Before and after this time range, the level of MLC activation by thrombin + histamine is less than the simple sum of that by thrombin and histamine alone. The less than additive effect during the first 10 min is primarily due to the time-dependent behavior of MLC activation by the Ca2+-dependent signalling cascade. The transient MLC activation curve by the Ca2+-dependent cascade is largely the same for the thrombin, histamine, and thrombin + histamine mediated processes (Additional File 1, Figure S11, S12 and S13, solid line). It is thus not difficult to understand that the simple sum of the level of MLC activation by thrombin and histamine alone is superficially larger than that by thrombin + histamine. The less than additive effect after 50 min is primarily due to the variation of time-dependent behavior of MLC activation by the ROCK-dependent signalling cascade. The level of MLC activation slowly rises to significant levels without decay in the presence of thrombin alone for up to 100 min (Additional File 1, Figure S11, dotted and dash-dotted line). On the other hand, the MLC activation level rises slowly to moderate levels without decay in the presence of histamine alone (Additional File 1, Figure S12). In contrast, the MLC activation level quickly rises to high levels and rapidly decays to low levels after 50 min in the presence of thrombin + histamine, the signalling strength thus becomes less than additive after 50 min (Additional File 1, Figure S13, dash-dotted line ).
The underling mechanism of the significant synergistic effect during 10-50 min time period can be elucidated from the perspective of network regulation. MLC activation is regulated by at least four signalling cascades Ca2+-dependent, CPI-17-dependent, NO-dependent, and ROCK-dependent cascades (Figure 1). As shown in Additional File 1, Figure S11, S12, and S13, The MLC activation curve induced by the Ca2+-dependent cascade is roughly the same for the thrombin, histamine and thrombin + histamine mediated processes. The level of MLC activation induced by the CPI-17-dependent, NO-dependent cascade in the presence of thrombin + histamine is close to the simple sum of that in the presence of thrombin and histamine alone. While the MLC activation induced by the ROCK-dependent cascade in the presence of thrombin + histamine is at significantly higher levels and shows more transient pattern than that in the presence of thrombin and histamine alone. These differences in signalling behavior lead to synergistic effect within 10 min to 50 min time range.
The different signalling behavior of the Rho-ROCK signalling stimulated by different mediators or mediator combinations primarily arises from the dynamics of ROCK activation . The kinase activity of ROCK is off when ROCK is intra-molecularly folded. ROCK can be activated only after it is unfolded by the binding of Rho GTPase to its Rho-binding domain to disrupt the auto-inhibitory interaction, which subsequently allows such proteins as Rho GTPase and PKG to activate ROCK at phosphorylation site. Hence, in the presence of thrombin + histamine, thrombin-activated Rho GTPase unfolds ROCK to allow histamine-activated PKG to activate ROCK thereby enhancing the level of ROCK activation in combination with thrombin-mediated Rho GTPase activation of ROCK. When stimulated by histamine or VEGF alone, ROCK is in the inactive state and does not contribute to MLC activation. When stimulated by thrombin alone, ROCK is activated only by Rho GTPase without the contribution from PKG, leading to a slower increase and lower peak strength of MLC activation than that in the presence of thrombin + histamine. Such an integrated communication network is expected to enable fine tuning of the strength and duration of MLC activation, thereby enabling fine regulation of physiological responses, including synergistic or more complex effects. Network models have suggested that partial inhibition of a surprisingly small number of targets can be more efficient than the complete inhibition of a single target . Experimental and simulation studies of synergistic effects of thrombin and histamine on endothelial monolayer permeability may provide useful information for developing multi-target drugs against endothelial permeability and related diseases .
Prediction of the collective regulation of MLC activation by thrombin and VEGF
Our simulation model was further used to study the collective regulation of MLC activation by thrombin and VEGF, with a particular focus on whether or not thrombin and VEGF synergistically activate MLC in certain time ranges. Systemic activation of blood coagulation is often present in cancer patients, and thrombin generated during thrombosis can augment malignant phenotypes by activating tumor cell adhesion to platelets and endothelial cells, enhancing tumor cell growth and metastasis, and stimulating tumor cell angiogenesis . Moreover, thrombin promotes VEGF secretion and platelet activation, thus causing a mutual stimulation between endothelial cells and cancer cells [76, 77]. Therefore, the collective effect of thrombin and VEGF on MLC activation and subsequently endothelial hyperpermability may have substantial influence on the tumor growth and metastasis process in cancer patients with blood coagulation near and at the tumor sites .
Comparison of the areas with respect to different time ranges in Figure 9
MLC activation curve
Area under MLC activation curve with respect to different time range (time unit:min)
Curve 1:Thrombin + VEGF induced activation
Curve 2:VEGF-mediated activation
Curve 3:Thrombin-mediated activation
Simple sum of curve 2 and 3
Area difference between curve 1 and simple sum of curve 2 and 3
Prediction of the effect of CPI-17 over-expression on MLC activation in the presence of lower concentration of thrombin, histamine and VEGF
Thrombin, VEGF, and histamine are hallmarks of endothelial hyper-permeability, which perform their regulatory roles individually and collectively under different disease conditions, and with different dynamic profiles. Thrombin and VEGF can increase microvascular permeability ~50,000 times more potently than histamine . Thrombin, VEGF, and histamine induce prolonged (1-1.5 hr), intermediate (15-20 min) and transient (~5 min) increases of endothelial permeability, respectively. An integrated simulation model that includes the signalling of all these hallmark mediators enables more comprehensive analysis of the signalling processes involved in different disease processes and regulated by different combinations of these mediators.
Based on published models of relevant signalling, we developed an integrated mathematical model including the signalling pathways of all three of these mediators. Simulation results from our model were consistent with available experimental data of signalling mediated by both individual mediators and combinations of two mediators, and could be used to interpret the sustained and transient phases of MLC activation. Our model was able to predict the effects of altered pathway components and synergistic combination of multiple mediators, some of which are consistent with experimental findings . Similar to the published models of other pathways, our model can potentially be used to identify important disease genes through sensitivity analysis of signalling components . Our model may also be extended to emphasize other components to facilitate further investigation of the effects of different mediators, cascades, and cross-talk on endothelial permeability and related diseases.
One of the most commonly used approaches to model biological systems is that of ODEs. In general, a differential equation can be used to describe the chemical reaction rate that depends on the change of participating species over time. The temporal dynamic behavior of molecular species in the biological signalling pathway network can be captured by a set of coupled ODEs. Our pathway model is illustrated in Figure 1. Thrombin, VEGF and histamine induced MLC activation, as well as Ca2+-dependent and ROCK-dependent activation of MLC, were included in the model. The constituent molecular interactions, their kinetic constants and molecular concentrations are described in detail in Additional File 2. The ODEs for these interactions were derived based on mass-action laws with interaction rate constants defined by the forward and reverse rate constants Kf and Kb or turnover number Kcat for enzymatic reactions derived from published models [88, 89] and other literature. Our simulation model contains 200 equations and interactions and 185 distinct molecular species, characterized by 319 kinetic parameters and 48 initial molecular concentrations. These ODEs were then solved by using the Dormand-Prince pair based Ode45 solver of Matlab with the absolute tolerance of 1.0E-6 and relative tolerance 0.0010. A Systems Biology Markup Language (SBML) version of our model is provided at http://bidd.nus.edu.sg/group/Supplement.htm, and uploaded into the BioModels  and KDBI  databases (Additional File 4).
Collection and estimation of kinetic parameters
The types of parameters used in our model are parameters governing protein-protein interactions and catalytic activities. The published simulation studies have found that most parameters are robust and moderate changes do not significantly alter the overall pathway behavior [15, 16, 18]. Apart from the use of the parameters of the published simulation models, additional parameters were obtained from the literature based on the widely used assumption that parameters measured in vitro and in some cell lines are generally applicable in most cases. For those protein-protein interactions without available parameters, their parameters were putatively estimated from the known parameters of the relevant interacting domain profile pairs [92, 93] or other interacting protein pairs of similar sequences. As a biological network is believed to be robust, and protein-protein binding interactions for proteins in similar families that mediate similar types of biochemical reactions (such as Ras and Rho) often differ within a 10-fold range, the values of kinetic parameters obtained from previous models were optimized within this range.
The parameters of the protein-pairs not available from previous models were obtained by the following procedure: The first step in finding the parameters of a specific protein-pair is to search protein-pairs that are both with available parameters and with each individual protein similar in sequence with the respective protein of the studied protein-pair. If one or more such protein-pairs are found, then the average values of the parameters of these protein-pairs are used as the initial parameters of the studied protein-pair, which are further optimized in ± 10-fold range with respect to experimental data. For instance, the parameters for the Rho activation cycle were obtained from the Ras activation cycle and were further optimized within a 10-fold range. The cycle of optimization and validation was repeated in order to obtain simulated results that agreed well with known experimental trends. If no such protein-pair is found, we proceed to the second step to search protein-pairs that are both with available parameters and with each individual protein belonging to the same domain family of the respective protein of the studied protein-pair. If one or more such protein-pairs are found, then the average values of the parameters of these protein-pairs are used as the initial parameters of the studied protein-pair, which are further optimized in ± 10-fold range with respect to experimental data.
The parameters of RhoGAP and PKC related protein-pairs were determined by the first and second step. The parameters of 14 protein-pairs could not be determined by the first and second step due to lack of experimental data and relevant protein-pairs with known parameters. These parameters were determined by using the trust regions algorithm  to fit the simulation results to the experimental data of RAS, ERK, MYPT and CPI-17 activation curves [95–97]. Additional File 1, Figure S14 shows the fitting curve against experimental RAS activation data. The level of fitting is based on the least-squares method and the fitting process proceeds in iterations until the R-square value is >0.6. In each iteration, the parameter values derived from the previous iteration were used as the starting parameters for further optimization.
Model Optimization, Validation and Parameter Sensitivity Analysis
Mathematical models developed by 'one-set-fits-all' generic parameters need not reproduce quantitative behavior in all systems, but may be able to reproduce the behavior or trend for specific systems. For instance, a mathematical model developed for a biological pathway from parameters obtained experimentally from one cell type can behave slightly different in another cell types. Differences in model behavior between cell types can be due to the presence or absence of crosstalk (i.e., differences in model topology) and variation in effective values of kinetic parameters. Hence, in this study, we developed a generic model of the thrombin-, VEGF-, and histamine-MLC signalling pathway to investigate the role of these three main mediators in regulating MLC activation. The simulated results are represented as trajectories of concentrations of chemical species with respect to time that are validated against available experimental data. If the trend or dynamic behavior of a particular reactant or product behaves in such a way that is consistent with the experimental data, then the model is said to be reasonable and can be used to analyze and predict unknown biological phenomena within some difficulty to define range of conditions. If the simulation results are not in agreement with known experimental facts, then the model has to be revisited to examine possible errors, such as incorrect interaction kinetics or values of kinetic parameters. Optimized parameters obtained from previous mathematical models are not necessarily optimized in a new study as the scope of these models can be different. The cycle of optimization and validation are repeated in order to obtain simulated results that agree well with known experimental trends.
The sensitivity of the simulation results with respect to the optimized and other parameters need to be systematically analyzed to determine if the model is sufficiently robust to be able to analyze and predict the true dynamic behavior of biological networks without the artifact of parameters. Differential analysis of parameter sensitivity, also referred to as the direct method, was utilized to compute the time-dependent sensitivities of all the species states with respect to each parameter values in the model . Complex-step derivative approximation  was used to calculate numerical derivatives of the reactions in the model to achieve near analytical accuracy, robustness and easy implementation. We used sensitivity analysis function of Matlab to conduct sensitivity analysis. The sensitivity value of ppMLC with respect to all parameters in the model was provided in Additional File 5 and Additional File 1, Figure S15. As shown in Additional File 1, Figure S15, only 14 (4%) kinetic parameters including CPI-17, PKC and ROCK related reactions showed some sensitivity in affecting the output. The majority of the parameters are insensitive in affecting the output. Thus, our model can be considered as sufficiently robust.
This work is supported by the Singapore-MIT Alliance (SMA), National University of Singapore, Grant (2007AA02Z160, 2009ZX09501-004) from the National Natural Science Foundation of China (NSFC), the Singapore-MIT Alliance for Research and Technology (SMART) through funding from the Singapore National Research Foundation. Ministry of Science and Technology Singapore and Key Special Project Grant 2009ZX09501-00 of China.
- Nierodzik ML, Karpatkin S: Thrombin induces tumor growth, metastasis, and angiogenesis: Evidence for a thrombin-regulated dormant tumor phenotype. Cancer Cell. 2006, 10: 355-362. 10.1016/j.ccr.2006.10.002View ArticlePubMedGoogle Scholar
- Martorell L, Martinez-Gonzalez J, Rodriguez C, Gentile M, Calvayrac O, Badimon L: Thrombin and protease-activated receptors (PARs) in atherothrombosis. Thromb Haemost. 2008, 99: 305-315.PubMedGoogle Scholar
- Finigan JH: The coagulation system and pulmonary endothelial function in acute lung injury. Microvasc Res. 2009, 77: 35-38. 10.1016/j.mvr.2008.09.002View ArticlePubMedGoogle Scholar
- Kumar P, Shen Q, Pivetti CD, Lee ES, Wu MH, Yuan SY: Molecular mechanisms of endothelial hyperpermeability: implications in inflammation. Expert Rev Mol Med. 2009, 11: e19-PubMed CentralView ArticlePubMedGoogle Scholar
- Jutel M, Blaser K, Akdis CA: The role of histamine in regulation of immune responses. Chem Immunol Allergy. 2006, 91: 174-187.View ArticlePubMedGoogle Scholar
- Coughlin SR: Thrombin signalling and protease-activated receptors. Nature. 2000, 407: 258-264. 10.1038/35025229View ArticlePubMedGoogle Scholar
- van Nieuw Amerongen GP, Draijer R, Vermeer MA, van Hinsbergh VW: Transient and prolonged increase in endothelial permeability induced by histamine and thrombin: role of protein kinases, calcium, and RhoA. Circ Res. 1998, 83: 1115-1123.View ArticlePubMedGoogle Scholar
- Keck PJ, Hauser SD, Krivi G, Sanzo K, Warren T, Feder J, Connolly DT: Vascular permeability factor, an endothelial cell mitogen related to PDGF. Science. 1989, 246: 1309-1312. 10.1126/science.2479987View ArticlePubMedGoogle Scholar
- Sun H, Breslin JW, Zhu J, Yuan SY, Wu MH: Rho and ROCK signaling in VEGF-induced microvascular endothelial hyperpermeability. Microcirculation. 2006, 13: 237-247. 10.1080/10739680600556944View ArticlePubMedGoogle Scholar
- Langley RR, Fidler IJ: Tumor cell-organ microenvironment interactions in the pathogenesis of cancer metastasis. Endocr Rev. 2007, 28: 297-321. 10.1210/er.2006-0027View ArticlePubMedGoogle Scholar
- Vandenbroucke E, Mehta D, Minshall R, Malik AB: Regulation of endothelial junctional permeability. Ann N Y Acad Sci. 2008, 1123: 134-145. 10.1196/annals.1420.016View ArticlePubMedGoogle Scholar
- Hirano K, Derkach DN, Hirano M, Nishimura J, Kanaide H: Protein kinase network in the regulation of phosphorylation and dephosphorylation of smooth muscle myosin light chain. Mol Cell Biochem. 2003, 248: 105-114. 10.1023/A:1024180101032View ArticlePubMedGoogle Scholar
- Wang L, Dudek SM: Regulation of vascular permeability by sphingosine 1-phosphate. Microvasc Res. 2009, 77: 39-45. 10.1016/j.mvr.2008.09.005PubMed CentralView ArticlePubMedGoogle Scholar
- Hu G, Minshall RD: Regulation of transendothelial permeability by Src kinase. Microvasc Res. 2009, 77: 21-25. 10.1016/j.mvr.2008.10.002View ArticlePubMedGoogle Scholar
- Lukas TJ: A signal transduction pathway model prototype I: From agonist to cellular endpoint. Biophys J. 2004, 87: 1406-1416. 10.1529/biophysj.103.035253PubMed CentralView ArticlePubMedGoogle Scholar
- Lukas TJ: A signal transduction pathway model prototype II: Application to Ca2+-calmodulin signaling and myosin light chain phosphorylation. Biophys J. 2004, 87: 1417-1425. 10.1529/biophysj.104.042721PubMed CentralView ArticlePubMedGoogle Scholar
- Moraru II, Loew LM: Intracellular signaling: spatial and temporal control. Physiology (Bethesda). 2005, 20: 169-179.View ArticleGoogle Scholar
- Maeda A, Ozaki Y, Sivakumaran S, Akiyama T, Urakubo H, Usami A, Sato M, Kaibuchi K, Kuroda S: Ca2+ -independent phospholipase A2-dependent sustained Rho-kinase activation exhibits all-or-none response. Genes Cells. 2006, 11: 1071-1083. 10.1111/j.1365-2443.2006.01001.xView ArticlePubMedGoogle Scholar
- Viswanathan GA, Seto J, Patil S, Nudelman G, Sealfon SC: Getting started in biological pathway construction and analysis. PLoS Comput Biol. 2008, 4: e16- 10.1371/journal.pcbi.0040016PubMed CentralView ArticlePubMedGoogle Scholar
- van Nieuw Amerongen GP, van Hinsbergh VW: Targets for pharmacological intervention of endothelial hyperpermeability and barrier function. Vascul Pharmacol. 2002, 39: 257-272. 10.1016/S1537-1891(03)00014-4View ArticlePubMedGoogle Scholar
- Fajmut A, Dobovisek A, Brumen M: Mathematical modeling of the relation between myosin phosphorylation and stress development in smooth muscles. J Chem Inf Model. 2005, 45: 1610-1615. 10.1021/ci050178aView ArticlePubMedGoogle Scholar
- Caunt M, Hu L, Tang T, Brooks PC, Ibrahim S, Karpatkin S: Growth-regulated oncogene is pivotal in thrombin-induced angiogenesis. Cancer Res. 2006, 66: 4125-4132. 10.1158/0008-5472.CAN-05-2570View ArticlePubMedGoogle Scholar
- Zania P, Papaconstantinou M, Flordellis CS, Maragoudakis ME, Tsopanoglou NE: Thrombin mediates mitogenesis and survival of human endothelial cells through distinct mechanisms. Am J Physiol Cell Physiol. 2008, 294: C1215-1226. 10.1152/ajpcell.00452.2007View ArticlePubMedGoogle Scholar
- Grand RJ, Turnell AS, Grabham PW: Cellular consequences of thrombin-receptor activation. Biochem J. 1996, 313 (Pt 2): 353-368.PubMed CentralView ArticlePubMedGoogle Scholar
- Vouret-Craviari V, Van Obberghen-Schilling E, Rasmussen UB, Pavirani A, Lecocq JP, Pouyssegur J: Synthetic alpha-thrombin receptor peptides activate G protein-coupled signaling pathways but are unable to induce mitogenesis. Mol Biol Cell. 1992, 3: 95-102.PubMed CentralView ArticlePubMedGoogle Scholar
- Parry MA, Myles T, Tschopp J, Stone SR: Cleavage of the thrombin receptor: identification of potential activators and inactivators. Biochem J. 1996, 320 (Pt 1): 335-341.PubMed CentralView ArticlePubMedGoogle Scholar
- Buhl AM, Johnson NL, Dhanasekaran N, Johnson GL: G alpha 12 and G alpha 13 stimulate Rho-dependent stress fiber formation and focal adhesion assembly. J Biol Chem. 1995, 270: 24631-24634. 10.1074/jbc.270.42.24631View ArticlePubMedGoogle Scholar
- Cobb MH: MAP kinase pathways. Prog Biophys Mol Biol. 1999, 71: 479-500. 10.1016/S0079-6107(98)00056-XView ArticlePubMedGoogle Scholar
- Aixinjueluo W, Furukawa K, Zhang Q, Hamamura K, Tokuda N, Yoshida S, Ueda R: Mechanisms for the apoptosis of small cell lung cancer cells induced by anti-GD2 monoclonal antibodies: roles of anoikis. J Biol Chem. 2005, 280: 29828-29836. 10.1074/jbc.M414041200View ArticlePubMedGoogle Scholar
- Shoemaker MO, Lau W, Shattuck RL, Kwiatkowski AP, Matrisian PE, Guerra-Santos L, Wilson E, Lukas TJ, Van Eldik LJ, Watterson DM: Use of DNA sequence and mutant analyses and antisense oligodeoxynucleotides to examine the molecular basis of nonmuscle myosin light chain kinase autoinhibition, calmodulin recognition, and activity. J Cell Biol. 1990, 111: 1107-1125. 10.1083/jcb.111.3.1107View ArticlePubMedGoogle Scholar
- Klemke RL, Cai S, Giannini AL, Gallagher PJ, de Lanerolle P, Cheresh DA: Regulation of cell motility by mitogen-activated protein kinase. J Cell Biol. 1997, 137: 481-492. 10.1083/jcb.137.2.481PubMed CentralView ArticlePubMedGoogle Scholar
- Bottazzi B, Vouret-Craviari V, Bastone A, De Gioia L, Matteucci C, Peri G, Spreafico F, Pausa M, D'Ettorre C, Gianazza E, Tagliabue A, Salmona M, Tedesco F, Introna M, Mantovani A: Multimer formation and ligand recognition by the long pentraxin PTX3. Similarities and differences with the short pentraxins C-reactive protein and serum amyloid P component. J Biol Chem. 1997, 272: 32817-32823. 10.1074/jbc.272.52.32817View ArticlePubMedGoogle Scholar
- Bredt DS, Snyder SH: Isolation of nitric oxide synthetase, a calmodulin-requiring enzyme. Proc Natl Acad Sci USA. 1990, 87: 682-685. 10.1073/pnas.87.2.682PubMed CentralView ArticlePubMedGoogle Scholar
- Lin S, Fagan KA, Li KX, Shaul PW, Cooper DM, Rodman DM: Sustained endothelial nitric-oxide synthase activation requires capacitative Ca2+ entry. J Biol Chem. 2000, 275: 17979-17985. 10.1074/jbc.275.24.17979View ArticlePubMedGoogle Scholar
- Noiri E, Lee E, Testa J, Quigley J, Colflesh D, Keese CR, Giaever I, Goligorsky MS: Podokinesis in endothelial cell migration: role of nitric oxide. Am J Physiol. 1998, 274: C236-244.PubMedGoogle Scholar
- Ziche M, Morbidelli L, Choudhuri R, Zhang HT, Donnini S, Granger HJ, Bicknell R: Nitric oxide synthase lies downstream from vascular endothelial growth factor-induced but not basic fibroblast growth factor-induced angiogenesis. J Clin Invest. 1997, 99: 2625-2634. 10.1172/JCI119451PubMed CentralView ArticlePubMedGoogle Scholar
- Hood J, Granger HJ: Protein kinase G mediates vascular endothelial growth factor-induced Raf-1 activation and proliferation in human endothelial cells. J Biol Chem. 1998, 273: 23504-23508. 10.1074/jbc.273.36.23504View ArticlePubMedGoogle Scholar
- Kamm KE, Stull JT: The function of myosin and myosin light chain kinase phosphorylation in smooth muscle. Annu Rev Pharmacol Toxicol. 1985, 25: 593-620. 10.1146/annurev.pa.25.040185.003113View ArticlePubMedGoogle Scholar
- Moussavi RS, Kelley CA, Adelstein RS: Phosphorylation of vertebrate nonmuscle and smooth muscle myosin heavy chains and light chains. Mol Cell Biochem. 1993, 127-128: 219-227. 10.1007/BF01076773View ArticlePubMedGoogle Scholar
- Somlyo AP, Somlyo AV: Signal transduction and regulation in smooth muscle. Nature. 1994, 372: 231-236. 10.1038/372231a0View ArticlePubMedGoogle Scholar
- Amano M, Ito M, Kimura K, Fukata Y, Chihara K, Nakano T, Matsuura Y, Kaibuchi K: Phosphorylation and activation of myosin by Rho-associated kinase (Rho-kinase). J Biol Chem. 1996, 271: 20246-20249. 10.1074/jbc.271.34.20246View ArticlePubMedGoogle Scholar
- Hartshorne DJ, Ito M, Erdodi F: Myosin light chain phosphatase: subunit composition, interactions and regulation. J Muscle Res Cell Motil. 1998, 19: 325-341. 10.1023/A:1005385302064View ArticlePubMedGoogle Scholar
- Alessi D, MacDougall LK, Sola MM, Ikebe M, Cohen P: The control of protein phosphatase-1 by targetting subunits. The major myosin phosphatase in avian smooth muscle is a novel form of protein phosphatase-1. Eur J Biochem. 1992, 210: 1023-1035. 10.1111/j.1432-1033.1992.tb17508.xView ArticlePubMedGoogle Scholar
- Shirazi A, Iizuka K, Fadden P, Mosse C, Somlyo AP, Somlyo AV, Haystead TA: Purification and characterization of the mammalian myosin light chain phosphatase holoenzyme. The differential effects of the holoenzyme and its subunits on smooth muscle. J Biol Chem. 1994, 269: 31598-31606.PubMedGoogle Scholar
- Johnson D, Cohen P, Chen MX, Chen YH, Cohen PT: Identification of the regions on the M110 subunit of protein phosphatase 1 M that interact with the M21 subunit and with myosin. Eur J Biochem. 1997, 244: 931-939. 10.1111/j.1432-1033.1997.00931.xView ArticlePubMedGoogle Scholar
- Kimura K, Ito M, Amano M, Chihara K, Fukata Y, Nakafuku M, Yamamori B, Feng J, Nakano T, Okawa K, Iwamatsu A, Kaibuchi K: Regulation of myosin phosphatase by Rho and Rho-associated kinase (Rho-kinase). Science. 1996, 273: 245-248. 10.1126/science.273.5272.245View ArticlePubMedGoogle Scholar
- Eto M, Ohmori T, Suzuki M, Furuya K, Morita F: A novel protein phosphatase-1 inhibitory protein potentiated by protein kinase C. Isolation from porcine aorta media and characterization. J Biochem. 1995, 118: 1104-1107.PubMedGoogle Scholar
- Goeckeler ZM, Wysolmerski RB: Myosin light chain kinase-regulated endothelial cell contraction: the relationship between isometric tension, actin polymerization, and myosin phosphorylation. J Cell Biol. 1995, 130: 613-627. 10.1083/jcb.130.3.613View ArticlePubMedGoogle Scholar
- Kolodney MS, Elson EL: Correlation of myosin light chain phosphorylation with isometric contraction of fibroblasts. J Biol Chem. 1993, 268: 23850-23855.PubMedGoogle Scholar
- Zhi G, Ryder JW, Huang J, Ding P, Chen Y, Zhao Y, Kamm KE, Stull JT: Myosin light chain kinase and myosin phosphorylation effect frequency-dependent potentiation of skeletal muscle contraction. Proc Natl Acad Sci USA. 2005, 102: 17519-17524. 10.1073/pnas.0506846102PubMed CentralView ArticlePubMedGoogle Scholar
- Benardeau A, Hatem SN, Rucker-Martin C, Le Grand B, Mace L, Dervanian P, Mercadier JJ, Coraboeuf E: Contribution of Na+/Ca2+ exchange to action potential of human atrial myocytes. Am J Physiol. 1996, 271: H1151-1161.PubMedGoogle Scholar
- Tran QK, Watanabe H: Calcium signalling in the endothelium. Handb Exp Pharmacol. 2006, 145-187.Google Scholar
- Jeng JH, Chan CP, Wu HL, Ho YS, Lee JJ, Liao CH, Chang YK, Chang HH, Chen YJ, Perng PJ, Chang MC: Protease-activated receptor-1-induced calcium signaling in gingival fibroblasts is mediated by sarcoplasmic reticulum calcium release and extracellular calcium influx. Cell Signal. 2004, 16: 731-740. 10.1016/j.cellsig.2003.11.008View ArticlePubMedGoogle Scholar
- Birukova AA, Smurova K, Birukov KG, Kaibuchi K, Garcia JG, Verin AD: Role of Rho GTPases in thrombin-induced lung vascular endothelial cells barrier dysfunction. Microvasc Res. 2004, 67: 64-77. 10.1016/j.mvr.2003.09.007View ArticlePubMedGoogle Scholar
- Yazaki A, Tamaru S, Sasaki Y, Komatsu N, Wada H, Shiku H, Nishikawa M: Inhibition by Rho-kinase and protein kinase C of myosin phosphatase is involved in thrombin-induced shape change of megakaryocytic leukemia cell line UT-7/TPO. Cell Signal. 2005, 17: 321-330. 10.1016/j.cellsig.2004.07.009View ArticlePubMedGoogle Scholar
- Kureishi Y, Kobayashi S, Amano M, Kimura K, Kanaide H, Nakano T, Kaibuchi K, Ito M: Rho-associated kinase directly induces smooth muscle contraction through myosin light chain phosphorylation. J Biol Chem. 1997, 272: 12257-12260. 10.1074/jbc.272.19.12257View ArticlePubMedGoogle Scholar
- Senger DR, Galli SJ, Dvorak AM, Perruzzi CA, Harvey VS, Dvorak HF: Tumor cells secrete a vascular permeability factor that promotes accumulation of ascites fluid. Science. 1983, 219: 983-985. 10.1126/science.6823562View ArticlePubMedGoogle Scholar
- Breslin JW, Pappas PJ, Cerveira JJ, Hobson RW, Duran WN: VEGF increases endothelial permeability by separate signaling pathways involving ERK-1/2 and nitric oxide. Am J Physiol Heart Circ Physiol. 2003, 284: H92-H100.View ArticlePubMedGoogle Scholar
- Nelken NA, Soifer SJ, O'Keefe J, Vu TK, Charo IF, Coughlin SR: Thrombin receptor expression in normal and atherosclerotic human arteries. J Clin Invest. 1992, 90: 1614-1621. 10.1172/JCI116031PubMed CentralView ArticlePubMedGoogle Scholar
- Tellez C, Bar-Eli M: Role and regulation of the thrombin receptor (PAR-1) in human melanoma. Oncogene. 2003, 22: 3130-3137. 10.1038/sj.onc.1206453View ArticlePubMedGoogle Scholar
- Even-Ram S, Uziely B, Cohen P, Grisaru-Granovsky S, Maoz M, Ginzburg Y, Reich R, Vlodavsky I, Bar-Shavit R: Thrombin receptor overexpression in malignant and physiological invasion processes. Nat Med. 1998, 4: 909-914. 10.1038/nm0898-909View ArticlePubMedGoogle Scholar
- Jin HG, Yamashita H, Nagano Y, Fukuba H, Hiji M, Ohtsuki T, Takahashi T, Kohriyama T, Kaibuchi K, Matsumoto M: Hypoxia-induced upregulation of endothelial small G protein RhoA and Rho-kinase/ROCK2 inhibits eNOS expression. Neurosci Lett. 2006, 408: 62-67. 10.1016/j.neulet.2006.08.038View ArticlePubMedGoogle Scholar
- Li B, Zhao WD, Tan ZM, Fang WG, Zhu L, Chen YH: Involvement of Rho/ROCK signalling in small cell lung cancer migration through human brain microvascular endothelial cells. FEBS Lett. 2006, 580: 4252-4260. 10.1016/j.febslet.2006.06.056View ArticlePubMedGoogle Scholar
- Price JT, Bonovich MT, Kohn EC: The biochemistry of cancer dissemination. Crit Rev Biochem Mol Biol. 1997, 32: 175-253. 10.3109/10409239709082573View ArticlePubMedGoogle Scholar
- Worthylake RA, Lemoine S, Watson JM, Burridge K: RhoA is required for monocyte tail retraction during transendothelial migration. J Cell Biol. 2001, 154: 147-160. 10.1083/jcb.200103048PubMed CentralView ArticlePubMedGoogle Scholar
- Adamson P, Etienne S, Couraud PO, Calder V, Greenwood J: Lymphocyte migration through brain endothelial cell monolayers involves signaling through endothelial ICAM-1 via a rho-dependent pathway. J Immunol. 1999, 162: 2964-2973.PubMedGoogle Scholar
- Ferrara N, Gerber HP, LeCouter J: The biology of VEGF and its receptors. Nat Med. 2003, 9: 669-676. 10.1038/nm0603-669View ArticlePubMedGoogle Scholar
- Chua CC, Hamdy RC, Chua BH: Upregulation of vascular endothelial growth factor by H2O2 in rat heart endothelial cells. Free Radic Biol Med. 1998, 25: 891-897. 10.1016/S0891-5849(98)00115-4View ArticlePubMedGoogle Scholar
- Fiallo P, Clapasson A, Favre A, Pesce C: Overexpression of vascular endothelial growth factor and its endothelial cell receptor KDR in type 1 leprosy reaction. Am J Trop Med Hyg. 2002, 66: 180-185.PubMedGoogle Scholar
- Caldwell RB, Bartoli M, Behzadian MA, El-Remessy AE, Al-Shabrawey M, Platt DH, Liou GI, Caldwell RW: Vascular endothelial growth factor and diabetic retinopathy: role of oxidative stress. Curr Drug Targets. 2005, 6: 511-524. 10.2174/1389450054021981View ArticlePubMedGoogle Scholar
- Ferrara N, Alitalo K: Clinical applications of angiogenic growth factors and their inhibitors. Nat Med. 1999, 5: 1359-1364. 10.1038/70928View ArticlePubMedGoogle Scholar
- Padro T, Bieker R, Ruiz S, Steins M, Retzlaff S, Burger H, Buchner T, Kessler T, Herrera F, Kienast J, Müller-Tidow C, Serve H, Berdel WE, Mesters RM: Overexpression of vascular endothelial growth factor (VEGF) and its cellular receptor KDR (VEGFR-2) in the bone marrow of patients with acute myeloid leukemia. Leukemia. 2002, 16: 1302-1310. 10.1038/sj.leu.2402534View ArticlePubMedGoogle Scholar
- Beynon HL, Haskard DO, Davies KA, Haroutunian R, Walport MJ: Combinations of low concentrations of cytokines and acute agonists synergize in increasing the permeability of endothelial monolayers. Clin Exp Immunol. 1993, 91: 314-319.PubMed CentralView ArticlePubMedGoogle Scholar
- Riento K, Ridley AJ: Rocks: multifunctional kinases in cell behaviour. Nat Rev Mol Cell Biol. 2003, 4: 446-456. 10.1038/nrm1128View ArticlePubMedGoogle Scholar
- Csermely P, Agoston V, Pongor S: The efficiency of multi-target drugs: the network approach might help drug design. Trends Pharmacol Sci. 2005, 26: 178-182. 10.1016/j.tips.2005.02.007View ArticlePubMedGoogle Scholar
- Maragoudakis ME, Tsopanoglou NE, Andriopoulou P, Maragoudakis MM: Effects of thrombin/thrombosis in angiogenesis and tumour progression. Matrix Biol. 2000, 19: 345-351. 10.1016/S0945-053X(00)00079-2View ArticlePubMedGoogle Scholar
- Roselli M, Mineo TC, Basili S, Mariotti S, Martini F, Bellotti A, Ambrogi V, Spila A, D'Alessandro R, Gazzaniga PP, Guadagni F, Ferroni P: Vascular endothelial growth factor (VEGF-A) plasma levels in non-small cell lung cancer: relationship with coagulation and platelet activation markers. Thromb Haemost. 2003, 89: 177-184.PubMedGoogle Scholar
- Gieseler F, Luhr I, Kunze T, Mundhenke C, Maass N, Erhart T, Denker M, Beckmann D, Tiemann M, Schulte C, Dohrmann P, Cavaillé F, Godeau F, Gespach C: Activated coagulation factors in human malignant effusions and their contribution to cancer cell metastasis and therapy. Thromb Haemost. 2007, 97: 1023-1030.PubMedGoogle Scholar
- MacDonald JA, Eto M, Borman MA, Brautigan DL, Haystead TA: Dual Ser and Thr phosphorylation of CPI-17, an inhibitor of myosin phosphatase, by MYPT-associated kinase. FEBS Lett. 2001, 493: 91-94. 10.1016/S0014-5793(01)02277-3View ArticlePubMedGoogle Scholar
- Ohama T, Hori M, Sato K, Ozaki H, Karaki H: Chronic treatment with interleukin-1beta attenuates contractions by decreasing the activities of CPI-17 and MYPT-1 in intestinal smooth muscle. J Biol Chem. 2003, 278: 48794-48804. 10.1074/jbc.M310166200View ArticlePubMedGoogle Scholar
- Sakai H, Chiba Y, Hirano T, Misawa M: Possible involvement of CPI-17 in augmented bronchial smooth muscle contraction in antigen-induced airway hyper-responsive rats. Mol Pharmacol. 2005, 68: 145-151.PubMedGoogle Scholar
- Dakshinamurti S, Mellow L, Stephens NL: Regulation of pulmonary arterial myosin phosphatase activity in neonatal circulatory transition and in hypoxic pulmonary hypertension: a role for CPI-17. Pediatr Pulmonol. 2005, 40: 398-407. 10.1002/ppul.20290View ArticlePubMedGoogle Scholar
- Chang S, Hypolite JA, DiSanto ME, Changolkar A, Wein AJ, Chacko S: Increased basal phosphorylation of detrusor smooth muscle myosin in alloxan-induced diabetic rabbit is mediated by upregulation of Rho-kinase beta and CPI-17. Am J Physiol Renal Physiol. 2006, 290: F650-656.View ArticlePubMedGoogle Scholar
- Woodsome TP, Eto M, Everett A, Brautigan DL, Kitazawa T: Expression of CPI-17 and myosin phosphatase correlates with Ca(2+) sensitivity of protein kinase C-induced contraction in rabbit smooth muscle. J Physiol. 2001, 535: 553-564. 10.1111/j.1469-7793.2001.t01-1-00553.xPubMed CentralView ArticlePubMedGoogle Scholar
- Aslam M, Hartel FV, Arshad M, Gunduz D, Abdallah Y, Sauer H, Piper HM, Noll T: cAMP/PKA antagonizes thrombin-induced inactivation of endothelial myosin light chain phosphatase: role of CPI-17. Cardiovasc Res. 2010, 87: 375-384. 10.1093/cvr/cvq065View ArticlePubMedGoogle Scholar
- Wang W, Merrill MJ, Borchardt RT: Vascular endothelial growth factor affects permeability of brain microvessel endothelial cells in vitro. Am J Physiol. 1996, 271: C1973-1980.PubMedGoogle Scholar
- Braun P, Tasan M, Dreze M, Barrios-Rodiles M, Lemmens I, Yu H, Sahalie JM, Murray RR, Roncari L, de Smet AS, Venkatesan K, Rual JF, Vandenhaute J, Cusick ME, Pawson T, Hill DE, Tavernier J, Wrana JL, Roth FP, Vidal M: An experimentally derived confidence score for binary protein-protein interactions. Nat Methods. 2009, 6: 91-97. 10.1038/nmeth.1281PubMed CentralView ArticlePubMedGoogle Scholar
- Kholodenko BN, Demin OV, Moehren G, Hoek JB: Quantification of short term signaling by the epidermal growth factor receptor. J Biol Chem. 1999, 274: 30169-30181. 10.1074/jbc.274.42.30169View ArticlePubMedGoogle Scholar
- Schoeberl B, Eichler-Jonsson C, Gilles ED, Muller G: Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors. Nat Biotechnol. 2002, 20: 370-375. 10.1038/nbt0402-370View ArticlePubMedGoogle Scholar
- Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, Chelliah V, Li L, He E, Henry A, Stefan MI, Snoep JL, Hucka M, Novère NL, Laibe C: BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Syst Biol. 2010, 4: 92- 10.1186/1752-0509-4-92PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar P, Han BC, Shi Z, Jia J, Wang YP, Zhang YT, Liang L, Liu QF, Ji ZL, Chen YZ: Update of KDBI: Kinetic Data of Bio-molecular Interaction database. Nucleic Acids Res. 2009, 37: D636-641. 10.1093/nar/gkn839PubMed CentralView ArticlePubMedGoogle Scholar
- Singhal M, Resat H: A domain-based approach to predict protein-protein interactions. BMC Bioinformatics. 2007, 8: 199- 10.1186/1471-2105-8-199PubMed CentralView ArticlePubMedGoogle Scholar
- Wojcik J, Schachter V: Protein-protein interaction map inference using interacting domain profile pairs. Bioinformatics. 2001, 17 (Suppl 1): S296-305. 10.1093/bioinformatics/17.suppl_1.S296View ArticlePubMedGoogle Scholar
- Coleman TF, Li YY: An interior trust region approach for nonlinear minimization subject to bounds. Siam J Optimiz. 1996, 6: 418-445. 10.1137/0806023.View ArticleGoogle Scholar
- Kawano Y, Fukata Y, Oshiro N, Amano M, Nakamura T, Ito M, Matsumura F, Inagaki M, Kaibuchi K: Phosphorylation of myosin-binding subunit (MBS) of myosin phosphatase by Rho-kinase in vivo. J Cell Biol. 1999, 147: 1023-1038. 10.1083/jcb.147.5.1023PubMed CentralView ArticlePubMedGoogle Scholar
- Sasagawa S, Ozaki Y, Fujita K, Kuroda S: Prediction and validation of the distinct dynamics of transient and sustained ERK activation. Nat Cell Biol. 2005, 7: 365-373. 10.1038/ncb1233View ArticlePubMedGoogle Scholar
- Dimopoulos GJ, Semba S, Kitazawa K, Eto M, Kitazawa T: Ca2+-dependent rapid Ca2+ sensitization of contraction in arterial smooth muscle. Circ Res. 2007, 100: 121-129. 10.1161/01.RES.0000253902.90489.dfPubMed CentralView ArticlePubMedGoogle Scholar
- Cameron AC, Windmeijer FAG: An R-squared measure of goodness of fit for some common nonlinear regression models. J Econometrics. 1997, 77: 329-342. 10.1016/S0304-4076(96)01818-0.View ArticleGoogle Scholar
- Hamby DM: A Review of Techniques for Parameter Sensitivity Analysis of Environmental-Models. Environ Monit Assess. 1994, 32: 135-154. 10.1007/BF00547132.View ArticlePubMedGoogle Scholar
- Martins JRRA, Sturdza P, Alonso JJ: The complex-step derivative approximation. Acm T Math Software. 2003, 29: 245-262. 10.1145/838250.838251.View ArticleGoogle Scholar