Time-resolved in silico modeling of fine-tuned cAMP signaling in platelets: feedback loops, titrated phosphorylations and pharmacological modulation
© Wangorsch et al; licensee BioMed Central Ltd. 2011
Received: 3 May 2011
Accepted: 28 October 2011
Published: 28 October 2011
Hemostasis is a critical and active function of the blood mediated by platelets. Therefore, the prevention of pathological platelet aggregation is of great importance as well as of pharmaceutical and medical interest. Endogenous platelet inhibition is predominantly based on cyclic nucleotides (cAMP, cGMP) elevation and subsequent cyclic nucleotide-dependent protein kinase (PKA, PKG) activation. In turn, platelet phosphodiesterases (PDEs) and protein phosphatases counterbalance their activity. This main inhibitory pathway in human platelets is crucial for countervailing unwanted platelet activation. Consequently, the regulators of cyclic nucleotide signaling are of particular interest to pharmacology and therapeutics of atherothrombosis. Modeling of pharmacodynamics allows understanding this intricate signaling and supports the precise description of these pivotal targets for pharmacological modulation.
We modeled dynamically concentration-dependent responses of pathway effectors (inhibitors, activators, drug combinations) to cyclic nucleotide signaling as well as to downstream signaling events and verified resulting model predictions by experimental data. Experiments with various cAMP affecting compounds including anti-platelet drugs and their combinations revealed a high fidelity, fine-tuned cAMP signaling in platelets without cross-talk to the cGMP pathway. The model and the data provide evidence for two independent feedback loops: PKA, which is activated by elevated cAMP levels in the platelet, subsequently inhibits adenylyl cyclase (AC) but as well activates PDE3. By multi-experiment fitting, we established a comprehensive dynamic model with one predictive, optimized and validated set of parameters. Different pharmacological conditions (inhibition, activation, drug combinations, permanent and transient perturbations) are successfully tested and simulated, including statistical validation and sensitivity analysis. Downstream cyclic nucleotide signaling events target different phosphorylation sites for cAMP- and cGMP-dependent protein kinases (PKA, PKG) in the vasodilator-stimulated phosphoprotein (VASP). VASP phosphorylation as well as cAMP levels resulting from different drug strengths and combined stimulants were quantitatively modeled. These predictions were again experimentally validated. High sensitivity of the signaling pathway at low concentrations is involved in a fine-tuned balance as well as stable activation of this inhibitory cyclic nucleotide pathway.
On the basis of experimental data, literature mining and database screening we established a dynamic in silico model of cyclic nucleotide signaling and probed its signaling sensitivity. Thoroughly validated, it successfully predicts drug combination effects on platelet function, including synergism, antagonism and regulatory loops.
Cyclic nucleotide signaling is the main inhibitory pathway in platelets in vivo that balances platelet activation. The complex regulation of this pathway includes endothelium released factors activating platelet nucleotide cyclases and in consequence cyclic nucleotide-dependent protein kinases that in turn phosphorylate major components of the platelet activation pathways thus preventing platelet aggregation . Platelet phosphodiesterases (PDEs) counterbalance the action of nucleotide cyclases [2–5]. Regulators of cAMP levels are of strong pharmacological interest and have clinical and therapeutic implications .
Computational modeling can advance understanding of cellular signaling in thrombosis and hemostasis and elucidates furthermore the fine-tuned cAMP signaling in platelets. Great effort has been devoted to model endothelial layer maintenance and senescence  as well as thrombus development . Here we model time-resolved cyclic nucleotide pathway signaling in human platelets involving specific parameters for metabolism, regulation and different effector strengths by means of a data-based in silico modeling approach.
Investigation of signaling pathways by kinetic modeling is commonly applied [8–11] and stimulates new biological insights into the modeled system. Combined with statistical validation based on experimental data, dynamic models can serve as a platform for testing hypotheses e.g. on pharmacological interventions on the system. Furthermore, the predictive features of kinetic models allow leading and directing biochemical experiments and elicit novel and unexpected biological findings.
In the following, our quantitative dynamic model, based on a system of ordinary differential equations (ODEs) is set up according to our philosophy (see below): Using a data-driven modeling approach, it enables us to correctly model the behavior of measured pathway components over time as well as to test the effects of drug combinations. It estimates individual effects of PDE-specific inhibitors and activators of adenylyl cyclase (AC) at different strengths and combinations and allows the prediction of feedback loops. We investigate cAMP accumulation after stimulation of adenylyl cyclase by Iloprost and Forskolin as well as the effects of inhibition of cAMP degrading PDEs by Cilostamide or Milrinone using individual concentration parameters for each platelet PDE isoform (PDE2, PDE3 and PDE5). The reliability of the model is tested and validated with experimental data of human platelets. It allows the prediction of effects of platelet inhibiting drugs and drug combinations on platelet cAMP pathway taking drug interactions, time-scales and regulatory loops into account.
Although the main players of inhibitory signal transduction in human platelets are well investigated, their interplay is still not fully understood. Since platelet inhibition and its up- and downstream signaling events influence this fine-tuned balance around platelet activation and quiescence, we here focus on a pharmacological motivated modeling approach, driven by data-based model validation.
Set-up of a basal model of unstimulated platelets
Model and experimental data indicate low effective PDE concentrations in resting human platelets
The central components of the platelet cyclic nucleotide pathways are the PDEs. Although, data on enzyme kinetics for many isoforms (including those expressed in the platelet) are available, no data on enzyme concentration in human platelets for any isoform are currently accessible. Therefore, we measured intracellular concentrations of the major phosphodiesterases in human platelets. After calibration with recombinant expressed PDEs by Western blot this yields concentrations of 63.46 mg/l (3.3 ng/107 platelets) for PDE2, 225 mg/l (11.7 ng/107 platelets) for PDE3 and 1359 mg/l (70.7 ng/107 platelets) for PDE5. Total molar concentrations are calculated based on a platelet volume of 5.2 fl as reported in .
Using the experimentally determined PDE concentrations and assuming that 80% of total cAMP turnover is provided by PDE3 activity  simulations from the basic model of the unstimulated platelet indicate that basal cyclic nucleotide levels cannot be maintained with these levels of active PDEs but are rapidly diminished (Additional file 1, Figure S2.1). In fact, mathematical calculations yielded a cAMP hydrolysis activity of 1846.1 μM/min for the PDE3 isoform (75% of total activity) as well as 642.9 μM/min for PDE2 (~ 25% of total activity). This would give total cAMP hydrolysis rate of 2.5 mM/min corresponding to a turnover rate of the entire platelet cAMP pool of 625 times per minute. These results suggest that not the entire amount of PDE is enzymatically active but rather the majority of the enzyme remains inactive in resting platelets. Parameter optimization constraining on a constant basal cAMP level of 4 μM yielded an enzymatically active PDE concentration of 0.05 mg/l for PDE2 (only 0.1% of total PDE2 concentration) and 2.3 mg/l for PDE3 (just 1% of total PDE2 concentration), respectively. This precisely reproduced cyclic nucleotide levels under resting conditions (Additional file 1, Figure S2.1). Consequently, this predicts an in vivo activity of 1.5 μM/min for PDE2 and of 6.5 μM/min for PDE3, yielding a total cAMP hydrolysis activity of 8 μM/min. This agrees with apparently low available concentrations suggested by electron microscopy (Additional file 1, Part III), though high total PDE concentrations were measured in vitro. The low active concentration for the antibody stain may be due to e.g. enzyme sequestration or inactivation in a complex (epitope masking) . Thus, either substantial amounts of cAMP degrading PDEs are not catalytically active and/or the in vivo enzyme activity is significantly lower than under in vitro conditions. Furthermore, there is experimental evidence indicating the compartmentalization of PDE3 and PDE5 [20–22]. To incorporate all experimental findings, we thus extended the model to include equilibrium between an active or inactive PDE form (Figure 1, Figure 3).
Model validation using experimental time series data from cAMP elevating compounds
After the development of the basal model of cyclic nucleotide signaling under resting conditions, we experimentally analyzed and modeled pathway effects of drugs causing elevated cyclic nucleotide levels in platelets thus inhibiting platelet aggregation (anti-platelet effect).
Model selection: Multi-experiment fitting and testing hypotheses on cAMP feedback loops
To incorporate all experimental findings, we iteratively expanded and refined the established basal model by fitting the resulting models reflecting PDE inhibition and AC stimulation simultaneously to all collected data sets:
Here, u i represents the concentration of the inhibitor (Milrinone, Cilostamide), whereas k i denotes the respective inhibition constant. These drug-specific constants adjust the apparent Km-value of the competitively inhibited enzyme.
respectively, mimicking the increased cAMP formation by AC (Figure 3C). This rate has been estimated locally for each concentration of the AC activators Forskolin and Iloprost.
Evaluation of integrated feedback loops within the dynamic model.
Factor kf1 thereby weights the cAMP feedback resulting in an activation of PDE3.
Feedback constant kf2 adjusts the strength of this negative feedback loop.
Hypothesis testing and model selection.
H1 versus H0
H2 versus H0
H3 versus H1
H3 versus H2
The validated in silico model reflects processes of platelet activation and models correctly the experimental measurements by considering low effective PDE concentrations, activation of PDE3 by PKA at increased cAMP levels as well as the negative feedback loop of cAMP-dependent PKA inhibiting AC activity.
We analyzed several characteristics of the platelet cAMP signal, indicating the state of inhibition, to investigate the influence of model parameters on this inhibitory signal. Results exemplified for one PDE3 inhibitor (Cilostamide) and one AC stimulator (Iloprost) indicate that in case of PDE3 inhibition (low and high dose; Additional file 2, Figure S1.1A, B), the cAMP signal is mainly controlled by the basal formation of cAMP. The impact of drug-specific and PDE-specific constants augments with the increase of the inhibitory drug dose. A similar sensitivity profile for parameters can be observed for low-dose AC stimulation (Additional file 2, Figure S1.1C) whereas at high doses, particularly the feedback constant kf2 (AC inhibition) gains importance and outweighs the influence of basal cAMP formation (Additional file 2, Figure S1.1D). This corroborates the pivotal role of feedback regulation of the cyclic nucleotide level in platelets.
To highlight the effect of the feedback constants kf1 (equation (3)) and kf2 (equation (4)) on the model predictions for cAMP signal, constants were perturbed (multiplication with factors ranging from 0.1 to 10). Constant kf1 (feedback activation of PDE3) resulted in cAMP signal spans ranging from 3 to 5 μM (Additional file 2, Figure S1.2A, B), however, this range expanded up to 300 μM (500 μM Forskolin). The cAMP signal is highly sensitive to AC feedback inhibition (constant kf2) even at low concentration (1 nM) of Iloprost (Additional file 2, Figure S1.2C). This is also reflected in the model selection procedure (Table 2) again indicating a subtle feedback regulation of cAMP levels.
Combination effects of different drugs
As synergistic effects of PDE inhibitors in combination with Iloprost have been described in human platelets [30–32] we embedded an additional cAMP-dependent constant k to capture this over-additive effect (Additional file 1, Section 8). In adaptation to quantitative data of experiments combining Cilostamide and Iloprost, we find the calculated cAMP level surface in close accordance to measured cAMP concentrations (Figure 5C), resulting in k = 1.443 ± 0.004. This way, the predicted cAMP levels describe perfectly the subsequent experimental measurements in contrast to assuming an additive effect of drug combinations (p < 0.05). This indicates a continuous synergistic effect of Cilostamide and Iloprost interaction potentiating the respective single drug stimulation (Figure 5A).
Furthermore we analyzed the combined effects of Forskolin and Milrinone (Figure 5D) by predicting cAMP levels of combining Forskolin drug doses (10, 30, 100, 200, 500 μM) with Milrinone doses of 10, 50 and 200 μM. The surface of cAMP level interpolates calculated cAMP level plateaus reached at the distinct combinations of drug. Similarly, we predict cAMP levels for the simultaneous inhibition of AC and PDE3, respectively. Figure 5E shows the successive decrease of elevated cAMP level due to Milrinone (1-100 μM) by combining this PDE3 inhibition with several doses of an AC inhibitor (e.g. 2'5'Dideoxyadenosine; non-synergistic, k = 1). In addition, we give a summary and classification of typical platelet drugs (Additional file 1, Table S8.1) allowing computer based modeling of drug (combination) effects.
Pathway integration: Phosphorylation of individual VASP phosphorylation sites
Cross-talk between cyclic nucleotide signaling pathways
Probing the signaling network sensitivity
To investigate the network sensitivity, we considered different perturbations and network cross-linking: We analyzed the model performance by probing drug stimuli (inhibition of PDE2 and PDE3) of different time-scales (transient vs. long-tem; concurrent vs. successive).
The network sensitively responds to transient stimuli even at low stimuli doses (Additional file 2, Figure S2.1). The same holds for transient prostaglandin receptor activation . Physiologically, this receptor is transiently activated by prostacyclin (half-life of less than 5 minutes ), produced by the endothelium. Thus, in case of an injury of the vessel wall parts of the endothelium remain without the production of prostacyclin so that only an insufficient signal is invoked and the platelet lacks inhibition. Similarly, patho-physiological conditions like prostacyclin receptor mutations [37, 38] contribute to a differential platelet inhibition decreasing the protection against unwanted platelet activation and aggregation. Modified platelet activation is crucial for thrombus self-organization and the formation of thrombi. We showed this gradual transient platelet inhibition via different prostacyclin doses as model response (cAMP) as well as experimentally, by considering cAMP-dependent, time-resolved phosphorylation events on VASP (Ser157) (Additional file 2, Figure S2.2).
Furthermore, we showed that an elevation of platelet cAMP level has no effect on the cGMP pathway, however the converse cross-talk has been described [1, 39, 40] and mainly mediated though cross-linking of the platelet phosphodiesterases. We therefore speculated on the effect of cGMP stimuli by incorporating identified PDE interconnections into the model structure (Additional file 2, Figure S2.3).
Mathematical modeling of platelet signaling
ODE-based dynamic models have emerged as fruitful tool for the modeling of biological systems and signaling pathways [8, 11, 41, 42]. In particular, they allow data-based kinetic modeling and the detailed investigation of pharmacological effects on cellular signaling cascades. In general, this works best for moderate-sized networks where sufficient detailed data are available to estimate model parameters. In this sense, the cyclic nucleotide signaling pathway in human platelets is especially well-suited for such an approach as there is enough information on the kinetics available and the number of key components is rather moderate.
Here, we present a mathematical model, based on existing biological knowledge and an extensive set of experimental measurements. Using the powerful approach of multi-experiment fitting, we optimized the model parameters in the context of several in vitro data sets. Since this provides a higher number of data points for the estimation of parameters, it permits not only the estimation of predictive model parameters but also statistically verifies the validity of a given model. Moreover, this allows for discrimination between competing model hypotheses [43, 44]. With this approach we could use the entire set of N = 155 experimental data points (triplicates) to obtain reliable estimates of the basic model parameters by globally fitting them to the dataset.
Pathway simulations and integration
Models on thrombosis or endothelial function have been published before [6, 7], but so far none investigated the cAMP mediated signaling in the platelet in detail. In this study, we established a time-resolved in silico model of fine-tuned cAMP signaling in platelets. The data and model can serve as a basis to gain a deeper insight into the basic effects of platelet cyclic nucleotide signaling. Predictions based on simulations of the basal in silico model indicate that only a very small fraction of the PDEs is enzymatically active under basal conditions which might be a result of PDE compartmentalization as has been reported for adipocytes (PDE3), HEK293 cells and cardiac myocytes (PDE4) and platelets (PDE5) [20–22, 45].
Elevated levels of cyclic nucleotides potently contribute to platelet inhibition targeting cyclic nucleotide-dependent protein kinases that in turn phosphorylate specific substrates. Via cAMP degrading phosphodiesterases, cGMP can regulate the level of cAMP in cardiac cells including platelets [1, 46]. Within this context, we show system state transitions from the non-activated cAMP pathway in platelets to an activated state as well as activatory and inhibitory cross-linking directed from the cGMP and the cAMP pathway (Additional file 1, Figure S1.2). Key feedback regulations involve the resting state of the platelet and the activation of both the cAMP pathway and cGMP pathway. Interestingly, the cAMP pathway is highly specific: Even at considerably high cAMP levels no cross-talk to the cGMP pathway could be observed. This is validated here including experimental measurements.
Furthermore, we investigated how the two branches of cyclic nucleotide signaling in human platelets might be integrated. Therefore, we focused on VASP, a prominent PKA and PKG substrate and highly connected hub protein. Signal integration by VASP protein relies on the activation of the downstream cytoskeletal regulating signaling cascade [35, 47]. Taking into account the two major known PKA/PKG phosphorylation sites, Ser157 and Ser239, we accurately predict and experimentally validate resulting differential activation in quantitative terms of specific PKA mediated phosphorylations on VASP.
Effects of drug combinations
Platelets play a critical role in thrombosis and hemostasis. Therefore, it is important that no accidental activation of platelets occurs and the intricate balance between activation and inhibition is maintained. This has implications for pharmacological fine-tuning of the system. Various pharmacological substances affecting cAMP signaling in platelets are in clinical use. The most prominent mechanisms of anti-platelet drugs are the inhibition of P2Y12 receptors (Clopridogrel, Prasugrel, Ticlopidine), of COX1 and thromboxane (Acetylsalicyclic acid) and of phosphodiesterases (Dipyridamole, Cilostazol, Milrinone) as well as antagonistic effects on integrin αIIbβ3 (Tirofiban, Eptifibatide, Abciximab) [48, 49]. Other drugs such as antagonists of platelet prostaglandin receptor EP3 for prostaglandin E2 targeting platelet G-protein-coupled receptors (GPCRs) or adenylyl cyclases as potential drug targets are under development [4, 50]. Furthermore, anti-platelet therapy benefits of the potential over-additive effect of platelet drugs applied in combination [48, 51] and of the investigation of drug interactions  potentially leading to a decrease in drug dosage and undesirable side effects like bleeding. Hence, a more detailed understanding of the interactions and effects of these substances on platelet activation and inhibition is not only of scientific interest but also of clinical importance. Having calibrated the model with experimental data of different drugs, the final model allows predictions not only about the effect of single drugs but also about combinations thereof. A major challenge of drug interaction modeling is to account for synergistic, additive or antagonistic drug interactions [53, 54]. This can be achieved by individual parameters for each specific combination of drug doses resulting in different types of interactions becoming apparent on the cAMP level. Concerning the modeled drug interactions, we succeeded in modeling these different modes of drug interaction with just one unique cAMP dependent parameter. The assumed simplification led to reliable and robust predictions compatible with the observed experimental data. However, we are aware that more complex effects and different modes of action have to be introduced if a broader range of drugs is considered and larger data sets are available.
Feedback regulation of cAMP level
Moreover, the system analysis based on experiments of platelet stimulation with anti-platelet drugs points to an important role of AC and PDE feedback primarily from PDE3. Focusing on major pathway components it becomes clear that these two feedback loops are to be involved to achieve optimal signal strengths. Induced by elevated cAMP concentration, activated PKA inhibits AC but as well activates PDE3 [25–27]. Short-term feedback regulation of cyclic nucleotide concentrations resulting from activation of cAMP-dependent PKA have been indicated in various cell types (rat hepatocytes, adipocytes, myoblasts, smooth muscle cells, osteoclasts, U937 cells) including platelets [24, 26, 55–59] showing the ubiquitous nature of this mechanism. Here, we focused on the activation of PDE3 as well as the inhibition of AC activity, both in PKA-dependent manner.
The activation of PDE3A and PDE4A by PKA has been observed in smooth muscle cells . Moreover, PDE3A is phosphorylated by PKA at Ser312 in a response to Forskolin  resulting in PDE3A activation in human platelets . Our model results suggest a more prominent role for AC inhibition by PKA phosphorylation. This PKA-dependent inhibitory effect has been described for adenylyl cyclase subtypes 5 (AC5) and 6 (AC6) in smooth muscle cells . As platelet proteomics data indicate the expression AC5 and AC6 besides AC subtype 3 and 7 in human platelets, direct phosphorylation of AC5 or AC6 by PKA might contribute to the feedback inhibition of cAMP synthesis [25, 28].
Our data suggest that both feedback effects are crucial for regulating cAMP levels in human platelets where the inhibitory effect on AC outweighs the positive feedback effect on PDE3 (Table 2).
We established a comprehensive dynamic model using multi-experiment fitting procedures, which enables us to evaluate and test questions of pharmacological interest. Our in silico model integrates central cAMP regulating mechanisms for a comprehensive description of cAMP signaling in human platelets. It has been tested and successfully applied to compare synergistic and non-synergistic as well as activatory and inhibitory drug combination effects on this fine-tuned platelet signaling pathway through which signals are sensitively propagated. The inhibitory cAMP pathway is well balanced: Our modeling approach reveals low active PDE concentrations compared to those experimentally measured in the platelet as well as two feedback loops allowing highly reliable signal transmission.
High sensitivity for low effector concentrations, threshold behavior and stable platelet inhibition are resulting. Perturbations such as prostacyclin receptor mutations and varying time scales of drug stimulation (transient vs. constant, successive vs. concordant) are modeled as well to probe the network sensitivity. In contrast to cGMP, which is known to interact with the cAMP pathway , converse cross-talk of cAMP to the cGMP pathway was neither predicted nor observed. However, direct integration over signals from cAMP and cGMP pathways is achieved via the integrator protein VASP and its phosphorylation sites, thereby also monitoring the effects of anti-platelet drugs.
Future work will extend the model to cover more complex scenarios such as stimulation of Gi-coupled receptors (e.g. P2Y12) which are known to decrease cAMP levels in human platelets thereby leading to platelet activation. Further model development will concern the integration of cGMP pathway effects on platelet cAMP levels upon stimulation of GCs.
In perspective, the in silico modeling approach of platelet cAMP as well as cGMP regulation will support future drug development as well as strategies for anti-aggregatory treatment and provides a unique tool for experimental design of pharmacological studies of platelets.
Washed human platelets  were prepared as described . Briefly, whole human blood was drawn after informed consent by venipuncture from volunteers who had not received any platelet or cyclic nucleotide affecting drugs in the past two weeks. Blood was collected in weakly acidic citrate/dextrose buffer (4:1 v/v). Whole blood was centrifuged at 300 xg for 10 min at 20°C. The platelet rich plasma obtained was centrifuged (380 xg, 20 min at 20°C). The resulting platelet pellet was resuspended in HEPES buffer (145 mM NaCl, 5 mM KCl, 1 mM MgCl2, 10 mM HEPES, 10 mM glucose, pH7.4) to a final cell density of 108/ml. Experiments were carried out by incubation of the platelet suspension at 37°C with the respective compounds and for the time. All experiments were reproduced with at least three different preparations.
Forskolin (Sigma, Munich, Germany), Iloprost (Schering, Berlin, Germany), Milrinone (Sigma, Munich, Germany) and Cilostamide (Tocris Cookson, Avonmouth, UK) were used as indicated in the results. PDE2A (C-terminal fragment aa1-207 deletion, C-terminal His6-tag), PDE3A (aa675-1164 fragment, N-terminal GST-tag) and PDE5A (full length, C-terminal His6) were a kind gift by J. A. Beavo, University of Seattle, Washington, USA.
Measurements of cAMP/cGMP levels
The stimulation was stopped by addition of twice the sample volume of ice cold ethanol (100%). The precipitate was removed by centrifugation, washed again with 100 μl ethanol; the supernatants were combined and dried in the vacuum. Dried platelet samples were dissolved, acetylated and analyzed with commercial cAMP or cGMP enzyme-immuno assays (ENZO Life Sciences, Loerrach, Germany) according to manufacturer's instructions.
Measurements of PDE concentrations
PDE2 and PDE5 were determined from calibration curves with lysates from recombinant expression systems . In brief, E. coli BL21-T1 strain was transformed with plasmids containing inserts coding for PDE2A, PDE3A or PDE5A and lysed after IPTG induction. The lysates were then purified by affinity chromatography. The GST tagged PDE3A fragment was bound to a GST-Sepharose column (Glutathione Sepharose 4B, Amersham Biosciences) and eluted according to the manufacturers' protocol. The His6-tagged PDE2A and PDE5A were purified on a Ni-NTA column (Ni-NTA Superflow Kit, Quiagen) and eluted according to the manufacturers' protocol. The protein content of the eluents was determined by the BCA assay (bicinchoninic acid, BCA protein assay, Thermo Scientific) and dilution series prepared. Platelet lysates and the dilution series of recombinant protein were then submitted to acrylamide gel-electrophoresis, transferred on a nitrocellulose membrane and detected with an appropriate antibody and chemoluminescence reaction (ECL plus, Amersham). The intensity of the protein bands was quantified after scanning with ImageJ software (v1.43, NIH) and the platelet expression determined from the calibration curves obtained from the dilution series. PDE2 and PDE5 were determined from calibration curves with lysates in over-expression systems . PDE3 amount was taken from , kinetic and binding data from [34, 63].
VASP phosphorylation was determined for both phosphorylation sites according to  with the phosphosite specific antibodies 16C2 (pS239) or 5C6 (pS157) and the total VASP antibody IE273. Stimulated washed platelet suspensions were stopped by addition of lysis buffer (20 mM Tris-HCl (pH 7.4), 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 0.5% NP-40, 10 mM ß-glycerolphosphate, 10 mM NaF) and frozen. The frozen samples were thawed and diluted with PBS buffer (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 2 mM KH2PO4, pH = 7.4) by 1:10. As control for background and non-specific binding 5% bovine serum albumin (BSA) solved in lysis buffer was used. Each sample was measured in triplicate for each antibody. VASP was immobilized from the sample by binding on a zyxin coated microtiter plate and incubation for 1 h at room temperature under shaking. The microtiter plates were washed three times with 300 μl/well PBS-T (1% Triton X-100 supplemented PBS) and the primary antibodies were added (5 μg/ml for the phosphospecific antibodies and 1 μg/ml for IE273) in 100 μl portions and again incubated for 1 h and washed 3 times with PBS-T. Detection of antibody binding was achieved with horse radish peroxidase (HRP) coupled goat anti-mouse IgG as secondary antibody and 100 μl of the HRP substrate ABTS (2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt) dissolved in 4 ml ABTS buffer (Roche) and diluted 1:10 with aqua dest. The absorbance of the samples was measured in the microtiter plates with a Wallac Victor 1420TM (Perkin-Elmer) plate reader at 405 nm each for 1 second. From the absorbance of each sample, the absorbance of the background control sample for the respective antibody was subtracted and the data expressed as relative phosphorylation calculated by dividing the absorbance signal obtained with the respective phosphospecific antibody by the signal obtained with the IE273 antibody.
In silico modeling
A set of ordinary differential equations (ODEs) represent platelet cAMP and cGMP signaling pathways (Additional file 1, Part II). Rate constants taken from literature all arise from studies on human platelets [12–20, 24, 27], inhibition constants are from different systems to delimit their fitting ranges. The basal model incorporated Mass Action and Michaelis Menten kinetics (Additional file 1, Part II, Section 2, Appendix). It was next optimized by fitting it simultaneously and experiment-specific to data of PDE inhibition experiments (Additional file 1, Part II, Section 3), AC stimulation experiments (Additional file 1, Part II, Section 4) as well as to phosphorylation data of VASP phosphosites (Additional file 1, Part II, Section 5, 6). Corresponding differential equations were implemented and further analyzed using the MATLAB software (The Mathworks, Inc. Natick, MA).
Estimating model parameters which optimize the χ2 -merit function and set the model statistically compliant with the data is a crucial problem . For integrating the dynamic model and optimal parameter estimation we used the MATLAB toolbox PottersWheel  including multi-experiment fitting. To fit the model to data, we optimize the χ2 -merit function ,with y i being data point i with standard deviation σ i and describes the model value at time point i for a set of parameter values . In case of normally distributed measurement errors, this corresponds to a maximum likelihood estimation. For optimizing this function, we used the trust region algorithm , a powerful deterministic least-square optimizer. This did yield optimal parameters for modeling e.g. the platelet effector experiments, minimizing the distance between model trajectories and time series data.
Model selection as hypothesis testing
For selecting an adequate model structure, being the most crucial part of the modeling process, we conduct the following forward strategy: We start with the most parsimonious reasonable model and refine it iteratively and directed by biochemical knowledge until subsequent refinement does not significantly improve the model fitting process. Therefore, we conducted a commonly used method for model comparison, the likelihood ratio test (LRT) comparing pairs of nested models characterized by a different number of parameters . Assuming a more complex model M C with p C parameters competes with the simpler model M S having p S inherent parameters. Then, the likelihood ratio (twice the difference of the log likelihoods L) is distributed as 2 (L(M C )-L(M S ))~χ2 with p C -p S degrees of freedom, which penalizes overparameterization. Given a chosen level of significance, this enables to determine whether the data supports the simple model (null hypothesis) or if the more complex model offers a sufficient improvement.
Appendix: Basic model - Mathematical formalization
Here, we discuss the basic model and its features in detail. This model, consisting of a system of ordinary differential equations is constructed to simulate the basal concentrations of cyclic nucleotides assuming resting platelet conditions.
Unstimulated platelets hold specific basal levels of the cyclic nucleotides cAMP and cGMP . There is a constant afflux of cyclic nucleotides arising from the corresponding adenylyl and guanylyl cyclase, respectively [17–19]. Active phosphodiesterases PDE2, PDE3 and the cGMP-specific PDE5 counterbalance their activity by enzymatically degrading cyclic nucleotides [2, 16]. The kinetics of these enzyme reactions is well understood as well as the cyclase activities. Thus, this prior biochemical knowledge frames the basis for setting up the basal signaling network as well as the underlying reactions as indicated in Figure 3A.
Dynamic variables and constants
Dynamic variables and constants of basal cyclic nucleotide signaling.
Experimentally: 63.46 mg/l
Experimentally: 225 mg/l
Experimentally: 1359 mg/l
(63.46 -x3) mg/l
(225 - x4) mg/l
(1359 - x5) mg/l
k1: Vmax PDE2
Basal influx of cAMP
Basal influx of cGMP
k13-k18: = 0
for simulation of
Kinetics and reaction rates
Modeled reactions and rates of basal cyclic nucleotide signaling.
Basal AC influx of cAMP (r1):
υ1 = k11;
Basal GC influx of cGMP (r2):
Υ2 = k12;
cAMP turnover via PDE2 (r3):
cAMP turnover via PDE3 (r4):
υ4 = k3·x4·x1/(k4+x1);
cGMP turnover via PDE2 (r5):
υ5 = k5·x3·x2/(k6+x2);
cGMP turnover via PDE3 (r6):
υ6 = k7·x4·x2/(k8+x2);
cGMP turnover via PDE5 (r7):
υ7 = k9·x5·x2/(k10+x2);
(De)activation of PDE (r8- r13):
υ8 = k14·x6;
υ9 = k13·x3;
υ10 = k16·x7;
υ11 = k15·x4;
υ12 = k18·x8;
υ13 = k17·x5;
This leads to the final system of differential equations of the dynamic variables x1-x10 as described in Table 3. A Systems Biology Markup Language file of this basal model is provided (Additional file 3) as well as its extension to the overall model (Additional file 4).
List of abbreviations
cyclic adenosine monophosphate
cyclic guanosine monophosphate
cAMP-dependent protein kinase
cGMP-dependent protein kinase
vasodilator stimulated phosphoprotein
ordinary differential equation
likelihood ratio test
Akaike information criterion
standard error of the mean.
The authors thank DFG (SFB688, TP A2, http://www.sfb688.uni-wuerzburg.de) and BMBF (SARA, FKZ 0315395D, FKZ0315395C, http://sara.informatik.uni-tuebingen.de) for support as well as Prof. Georg Krohne, University of Wuerzburg, for his scientific input into the electron microscopy analysis and Prof. Ulrich Walter, University hospital of Wuerzburg, for helpful and critical discussions on platelet signaling. This publication was funded by the German Research Foundation (DFG) and the University of Wuerzburg in the funding program Open Access Publishing).
- Schwarz UR, Walter U, Eigenthaler M: Taming platelets with cyclic nucleotides. Biochemical pharmacology. 2001, 62: 1153-1161. 10.1016/S0006-2952(01)00760-2.View ArticlePubMedGoogle Scholar
- Bender AT, Beavo JA: Cyclic nucleotide phosphodiesterases: molecular regulation to clinical use. Pharmacol Rev. 2006, 58: 488-520. 10.1124/pr.58.3.5.View ArticlePubMedGoogle Scholar
- Feijge MA, Ansink K, Vanschoonbeek K, Heemskerk JW: Control of platelet activation by cyclic AMP turnover and cyclic nucleotide phosphodiesterase type-3. Biochem Pharmacol. 2004, 67: 1559-1567. 10.1016/j.bcp.2003.12.028.View ArticlePubMedGoogle Scholar
- Noe L, Peeters K, Izzi B, Van Geet C, Freson K: Regulators of platelet cAMP levels: clinical and therapeutic implications. Curr Med Chem. 2010, 17: 2897-2905. 10.2174/092986710792065018.View ArticlePubMedGoogle Scholar
- Geiger J: Inhibitors of platelet signal transduction as anti-aggregatory drugs. Expert Opin Investig Drugs. 2001, 10: 865-890. 10.1517/135437126.96.36.1995.View ArticlePubMedGoogle Scholar
- Wang Y, Aguda BD, Friedman A: A continuum mathematical model of endothelial layer maintenance and senescence. Theor Biol Med Model. 2007, 4: 30-10.1186/1742-4682-4-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu Z, Chen N, Kamocka MM, Rosen ED, Alber M: A multiscale model of thrombus development. J R Soc Interface. 2008, 5: 705-722. 10.1098/rsif.2007.1202.PubMed CentralView ArticlePubMedGoogle Scholar
- Swameye I, Muller TG, Timmer J, Sandra O, Klingmuller U: Identification of nucleocytoplasmic cycling as a remote sensor in cellular signaling by databased modeling. Proc Natl Acad Sci USA. 2003, 100: 1028-1033. 10.1073/pnas.0237333100.PubMed CentralView ArticlePubMedGoogle Scholar
- Schilling M, Maiwald T, Hengl S, Winter D, Kreutz C, Kolch W, Lehmann WD, Timmer J, Klingmuller U: Theoretical and experimental analysis links isoform-specific ERK signalling to cell fate decisions. Mol Syst Biol. 2009, 5: 334-PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh MY, Yang S, Raymond-Stinz MA, Edwards JS, Wilson BS: Spatio-temporal modeling of signaling protein recruitment to EGFR. BMC Syst Biol. 2010, 4: 57-10.1186/1752-0509-4-57.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffmann A, Levchenko A, Scott ML, Baltimore D: The IkappaB-NF-kappaB signaling module: temporal control and selective gene activation. Science. 2002, 298: 1241-1245. 10.1126/science.1071914.View ArticlePubMedGoogle Scholar
- Dittrich M, Birschmann I, Mietner S, Sickmann A, Walter U, Dandekar T: Platelet protein interactions: map, signaling components, and phosphorylation groundstate. Arterioscler Thromb Vasc Biol. 2008, 28: 1326-1331. 10.1161/ATVBAHA.107.161000.View ArticlePubMedGoogle Scholar
- Dittrich M, Birschmann I, Pfrang J, Herterich S, Smolenski A, Walter U, Dandekar T: Analysis of SAGE data in human platelets: features of the transcriptome in an anucleate cell. Thromb Haemost. 2006, 95: 643-651.PubMedGoogle Scholar
- Eigenthaler M, Nolte C, Halbrugge M, Walter U: Concentration and regulation of cyclic nucleotides, cyclic-nucleotide-dependent protein kinases and one of their major substrates in human platelets. Estimating the rate of cAMP-regulated and cGMP-regulated protein phosphorylation in intact cells. Eur J Biochem. 1992, 205: 471-481. 10.1111/j.1432-1033.1992.tb16803.x.View ArticlePubMedGoogle Scholar
- Matsumoto T, Kobayashi T, Kamata K: Phosphodiesterases in the vascular system. J Smooth Muscle Res. 2003, 39: 67-86. 10.1540/jsmr.39.67.View ArticlePubMedGoogle Scholar
- Butt E, Walter U: Platelet phosphodiesterases. Handbook of experimental pharmacology. 1997, 126: 219-230. 10.1007/978-3-642-60639-7_10.View ArticleGoogle Scholar
- Hines LM, Tabakoff B: Platelet adenylyl cyclase activity: a biological marker for major depression and recent drug use. Biol Psychiatry. 2005, 58: 955-962. 10.1016/j.biopsych.2005.05.040.View ArticlePubMedGoogle Scholar
- Parsian A, Todd RD, Cloninger CR, Hoffman PL, Ovchinnikova L, Ikeda H, Tabakoff B: Platelet adenylyl cyclase activity in alcoholics and subtypes of alcoholics. WHO/ISBRA Study Clinical Centers. Alcohol Clin Exp Res. 1996, 20: 745-751. 10.1111/j.1530-0277.1996.tb01681.x.View ArticlePubMedGoogle Scholar
- Juska A, Farndale RW: Inhibition of human platelet adenylate cyclase activity by adrenaline, thrombin and collagen: analysis and reinterpretation of experimental data. Biochem J. 1999, 340 (Pt 1): 245-253.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson LS, Elbatarny HS, Crawley SW, Bennett BM, Maurice DH: Compartmentation and compartment-specific regulation of PDE5 by protein kinase G allows selective cGMP-mediated regulation of platelet functions. Proc Natl Acad Sci USA. 2008, 105: 13650-13655. 10.1073/pnas.0804738105.PubMed CentralView ArticlePubMedGoogle Scholar
- Baillie GS: Compartmentalized signalling: spatial regulation of cAMP by the action of compartmentalized phosphodiesterases. FEBS J. 2009, 276: 1790-1799. 10.1111/j.1742-4658.2009.06926.x.View ArticlePubMedGoogle Scholar
- Houslay MD: Underpinning compartmentalised cAMP signalling through targeted cAMP breakdown. Trends Biochem Sci. 2010, 35: 91-100. 10.1016/j.tibs.2009.09.007.View ArticlePubMedGoogle Scholar
- Floreani M, Fossa P, Gessi S, Mosti L, Borea PA, Dorigo P: New milrinone analogues: in vitro study of structure-activity relationships for positive inotropic effect, antagonism towards endogenous adenosine, and inhibition of cardiac type III phosphodiesterase. Naunyn Schmiedebergs Arch Pharmacol. 2003, 367: 109-118. 10.1007/s00210-002-0675-2.View ArticlePubMedGoogle Scholar
- Macphee CH, Reifsnyder DH, Moore TA, Lerea KM, Beavo JA: Phosphorylation results in activation of a cAMP phosphodiesterase in human platelets. J Biol Chem. 1988, 263: 10353-10358.PubMedGoogle Scholar
- Iwami G, Kawabe J, Ebina T, Cannon PJ, Homcy CJ, Ishikawa Y: Regulation of adenylyl cyclase by protein kinase A. J Biol Chem. 1995, 270: 12481-12484. 10.1074/jbc.270.21.12481.View ArticlePubMedGoogle Scholar
- Murthy KS, Zhou H, Makhlouf GM: PKA-dependent activation of PDE3A and PDE4 and inhibition of adenylyl cyclase V/VI in smooth muscle. Am J Physiol Cell Physiol. 2002, 282: C508-C517.View ArticlePubMedGoogle Scholar
- Hunter RW, Mackintosh C, Hers I: Protein kinase C-mediated phosphorylation and activation of PDE3A regulate cAMP levels in human platelets. J Biol Chem. 2009, 284: 12339-12348. 10.1074/jbc.M807536200.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Harry A, Li J, Smit MJ, Bai X, Magnusson R, Pieroni JP, Weng G, Iyengar R: Adenylyl cyclase 6 is selectively regulated by protein kinase A phosphorylation in a region involved in Galphas stimulation. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94: 14100-14104. 10.1073/pnas.94.25.14100.PubMed CentralView ArticlePubMedGoogle Scholar
- Cox DR, Hinkley DV: Theoretical Statistics. 1994, New York: Chapman and HallGoogle Scholar
- Nolte C, Eigenthaler M, Horstrup K, Honig-Liedl P, Walter U: Synergistic phosphorylation of the focal adhesion-associated vasodilator-stimulated phosphoprotein in intact human platelets in response to cGMP- and cAMP-elevating platelet inhibitors. Biochem Pharmacol. 1994, 48: 1569-1575. 10.1016/0006-2952(94)90201-1.View ArticlePubMedGoogle Scholar
- Hanson MS, Stephenson AH, Bowles EA, Sridharan M, Adderley S, Sprague RS: Phosphodiesterase 3 is present in rabbit and human erythrocytes and its inhibition potentiates iloprost-induced increases in cAMP. Am J Physiol Heart Circ Physiol. 2008, 295: H786-793. 10.1152/ajpheart.00349.2008.PubMed CentralView ArticlePubMedGoogle Scholar
- Netherton SJ, Jimmo SL, Palmer D, Tilley DG, Dunkerley HA, Raymond DR, Russell JC, Absher PM, Sage EH, Vernon RB, et al, et al.: Altered phosphodiesterase 3-mediated cAMP hydrolysis contributes to a hypermotile phenotype in obese JCR:LA-cp rat aortic vascular smooth muscle cells: implications for diabetes-associated cardiovascular disease. Diabetes. 2002, 51: 1194-1200. 10.2337/diabetes.51.4.1194.View ArticlePubMedGoogle Scholar
- Smolenski A, Bachmann C, Reinhard K, Honig-Liedl P, Jarchau T, Hoschuetzky H, Walter U: Analysis and regulation of vasodilator-stimulated phosphoprotein serine 239 phosphorylation in vitro and in intact cells using a phosphospecific monoclonal antibody. J Biol Chem. 1998, 273: 20029-20035. 10.1074/jbc.273.32.20029.View ArticlePubMedGoogle Scholar
- Butt E, Abel K, Krieger M, Palm D, Hoppe V, Hoppe J, Walter U: cAMP- and cGMP-dependent protein kinase phosphorylation sites of the focal adhesion vasodilator-stimulated phosphoprotein (VASP) in vitro and in intact human platelets. J Biol Chem. 1994, 269: 14509-14517.PubMedGoogle Scholar
- Geiger J, Brandmann T, Hubertus K, Tjahjadi B, Schinzel R, Walter U: A protein phosphorylation-based assay for screening and monitoring of drugs modulating cyclic nucleotide pathways. Anal Biochem. 2010, 407: 261-269. 10.1016/j.ab.2010.08.006.View ArticlePubMedGoogle Scholar
- Gomberg-Maitland M, Olschewski H: Prostacyclin therapies for the treatment of pulmonary arterial hypertension. Eur Respir J. 2008, 31: 891-901. 10.1183/09031936.00097107.View ArticlePubMedGoogle Scholar
- Patrignani P, Di Febbo C, Tacconelli S, Douville K, Guglielmi MD, Horvath RJ, Ding M, Sierra K, Stitham J, Gleim S, et al, et al.: Differential association between human prostacyclin receptor polymorphisms and the development of venous thrombosis and intimal hyperplasia: a clinical biomarker study. Pharmacogenet Genomics. 2008, 18: 611-620. 10.1097/FPC.0b013e328301a774.View ArticlePubMedGoogle Scholar
- Arehart E, Stitham J, Asselbergs FW, Douville K, MacKenzie T, Fetalvero KM, Gleim S, Kasza Z, Rao Y, Martel L, et al, et al.: Acceleration of cardiovascular disease by a dysfunctional prostacyclin receptor mutation: potential implications for cyclooxygenase-2 inhibition. Circulation research. 2008, 102: 986-993. 10.1161/CIRCRESAHA.107.165936.PubMed CentralView ArticlePubMedGoogle Scholar
- Omori K, Kotera J: Overview of PDEs and their regulation. Circulation research. 2007, 100: 309-327. 10.1161/01.RES.0000256354.95791.f1.View ArticlePubMedGoogle Scholar
- Aktas B, Honig-Liedl P, Walter U, Geiger J: Inhibition of platelet P2Y12 and alpha2A receptor signaling by cGMP-dependent protein kinase. Biochemical pharmacology. 2002, 64: 433-439. 10.1016/S0006-2952(02)01113-9.View ArticlePubMedGoogle Scholar
- Witt J, Barisic S, Schumann E, Allgower F, Sawodny O, Sauter T, Kulms D: Mechanism of PP2A-mediated IKK beta dephosphorylation: a systems biological approach. BMC Syst Biol. 2009, 3: 71-10.1186/1752-0509-3-71.PubMed CentralView ArticlePubMedGoogle Scholar
- van Riel NA: Dynamic modelling and analysis of biochemical networks: mechanism-based models and model-based experiments. Brief Bioinform. 2006, 7: 364-374. 10.1093/bib/bbl040.View ArticlePubMedGoogle Scholar
- Maiwald T, Timmer J: Dynamical modeling and multi-experiment fitting with PottersWheel. Bionformatics. 2008, 24: 2037-2043. 10.1093/bioinformatics/btn350.View ArticleGoogle Scholar
- Coleman TF, Li AY: An Interior Trust Region Approach for Nonlinear Minimization Subject to Bounds. SIAM Journal on Optimization. 1996, 6: 418-455. 10.1137/0806023.View ArticleGoogle Scholar
- Shakur Y, Takeda K, Kenan Y, Yu ZX, Rena G, Brandt D, Houslay MD, Degerman E, Ferrans VJ, Manganiello VC: Membrane localization of cyclic nucleotide phosphodiesterase 3 (PDE3). Two N-terminal domains are required for the efficient targeting to, and association of, PDE3 with endoplasmic reticulum. The Journal of biological chemistry. 2000, 275: 38749-38761.View ArticlePubMedGoogle Scholar
- Zaccolo M, Movsesian MA: cAMP and cGMP signaling cross-talk: role of phosphodiesterases and implications for cardiac pathophysiology. Circ Res. 2007, 100: 1569-1578. 10.1161/CIRCRESAHA.106.144501.View ArticlePubMedGoogle Scholar
- Dittrich M, Strassberger V, Fackler M, Tas P, Lewandrowski U, Sickmann A, Walter U, Dandekar T, Birschmann I: Characterization of a novel interaction between vasodilator-stimulated phosphoprotein and Abelson interactor 1 in human platelets: a concerted computational and experimental approach. Arterioscler Thromb Vasc Biol. 2010, 30: 843-850. 10.1161/ATVBAHA.109.200683.View ArticlePubMedGoogle Scholar
- Michelson AD: Antiplatelet therapies for the treatment of cardiovascular disease. Nat Rev Drug Discov. 2010, 9: 154-169. 10.1038/nrd2957.View ArticlePubMedGoogle Scholar
- Rao YJ, Xi L: Pivotal effects of phosphodiesterase inhibitors on myocyte contractility and viability in normal and ischemic hearts. Acta Pharmacol Sin. 2009, 30: 1-24. 10.1038/aps.2008.1.PubMed CentralView ArticlePubMedGoogle Scholar
- Pierre S, Eschenhagen T, Geisslinger G, Scholich K: Capturing adenylyl cyclases as potential drug targets. Nat Rev Drug Discov. 2009, 8: 321-335. 10.1038/nrd2827.View ArticlePubMedGoogle Scholar
- Lobato EB, Beaver T, Muehlschlegel J, Kirby DS, Klodell C, Sidi A: Treatment with phosphodiesterase inhibitors type III and V: milrinone and sildenafil is an effective combination during thromboxane-induced acute pulmonary hypertension. Br J Anaesth. 2006, 96: 317-322. 10.1093/bja/ael009.View ArticlePubMedGoogle Scholar
- Mackenzie IS, Coughtrie MW, MacDonald TM, Wei L: Antiplatelet drug interactions. J Intern Med. 2010, 268: 516-529. 10.1111/j.1365-2796.2010.02299.x.View ArticlePubMedGoogle Scholar
- Yan H, Zhang B, Li S, Zhao Q: A formal model for analyzing drug combination effects and its application in TNF-alpha-induced NFkappaB pathway. BMC Syst Biol. 2010, 4: 50-10.1186/1752-0509-4-50.PubMed CentralView ArticlePubMedGoogle Scholar
- Geva-Zatorsky N, Dekel E, Cohen AA, Danon T, Cohen L, Alon U: Protein dynamics in drug combinations: a linear superposition of individual-drug responses. Cell. 2010, 140: 643-651. 10.1016/j.cell.2010.02.011.View ArticlePubMedGoogle Scholar
- Corbin JD, Beebe SJ, Blackmore PF: cAMP-dependent protein kinase activation lowers hepatocyte cAMP. The Journal of biological chemistry. 1985, 260: 8731-8735.PubMedGoogle Scholar
- Gettys TW, Vine AJ, Simonds MF, Corbin JD: Activation of the particulate low Km phosphodiesterase of adipocytes by addition of cAMP-dependent protein kinase. The Journal of biological chemistry. 1988, 263: 10359-10363.PubMedGoogle Scholar
- Ball EH, Seth PK, Sanwal BD: Regulatory mechanisms involved in the control of cyclic adenosine 3':5'-monophosphate phosphodiesterases in myoblasts. The Journal of biological chemistry. 1980, 255: 2962-2968.PubMedGoogle Scholar
- Alvarez R, Sette C, Yang D, Eglen RM, Wilhelm R, Shelton ER, Conti M: Activation and selective inhibition of a cyclic AMP-specific phosphodiesterase, PDE-4D3. Mol Pharmacol. 1995, 48: 616-622.PubMedGoogle Scholar
- Ahlstrom M, Lamberg-Allardt C: Rapid protein kinase A--mediated activation of cyclic AMP-phosphodiesterase by parathyroid hormone in UMR-106 osteoblast-like cells. J Bone Miner Res. 1997, 12: 172-178. 10.1359/jbmr.19188.8.131.52.View ArticlePubMedGoogle Scholar
- Pozuelo Rubio M, Campbell DG, Morrice NA, Mackintosh C: Phosphodiesterase 3A binds to 14-3-3 proteins in response to PMA-induced phosphorylation of Ser428. Biochem J. 2005, 392: 163-172. 10.1042/BJ20051103.PubMed CentralView ArticlePubMedGoogle Scholar
- Aktas B, Utz A, Hoenig-Liedl P, Walter U, Geiger J: Dipyridamole enhances NO/cGMP-mediated vasodilator-stimulated phosphoprotein phosphorylation and signaling in human platelets: in vitro and in vivo/ex vivo studies. Stroke. 2003, 34: 764-769. 10.1161/01.STR.0000056527.34434.59.View ArticlePubMedGoogle Scholar
- Geiger J, Honig-Liedl P, Schanzenbacher P, Walter U: Ligand specificity and ticlopidine effects distinguish three human platelet ADP receptors. Eur J Pharmacol. 1998, 351: 235-246. 10.1016/S0014-2999(98)00305-7.View ArticlePubMedGoogle Scholar
- Poppe H, Rybalkin SD, Rehmann H, Hinds TR, Tang XB, Christensen AE, Schwede F, Genieser HG, Bos JL, Doskeland SO, et al, et al.: Cyclic nucleotide analogs as probes of signaling pathways. Nat Methods. 2008, 5: 277-278. 10.1038/nmeth0408-277.View ArticlePubMedGoogle Scholar
- Schmidt K, Schrammel A, Koesling D, Mayer B: Molecular mechanisms involved in the synergistic activation of soluble guanylyl cyclase by YC-1 and nitric oxide in endothelial cells. Mol Pharmacol. 2001, 59: 220-224.PubMedGoogle Scholar
- Wu CC, Ko FN, Kuo SC, Lee FY, Teng CM: YC-1 inhibited human platelet aggregation through NO-independent activation of soluble guanylate cyclase. Br J Pharmacol. 1995, 116: 1973-1978.PubMed CentralView ArticlePubMedGoogle Scholar