A modular systems biological modelling framework studies cyclic nucleotide signaling in platelets.

BACKGROUND: The cyclic nucleotides cAMP and cGMP inhibit platelet activation. Different platelet signaling modules work together. We develop here a modelling framework to integrate different signaling modules and apply it to platelets. RESULTS: We introduce a novel standardized bilinear coupling mechanism allowing sub model debugging and standardization of coupling with optimal data driven modelling by methods from optimization. Besides cAMP signaling our model considers specific cGMP effects including external stimuli by drugs. Moreover, the output of the cGMP module serves as input for a modular model of VASP phosphorylation and for the activity of cAMP and cGMP pathways in platelets. Experimental data driven modeling allows us to design models with quantitative output. We use the condensed information about involved regulation and system responses for modeling drug effects and obtaining optimal experimental settings. Stepwise further validation of our model is given by direct experimental data. CONCLUSIONS: We present a general framework for model integration using modules and their stimulus responses. We demonstrate it by a multi-modular model for platelet signaling focusing on cGMP and VASP phosphorylation. Moreover, this allows to estimate drug action on any of the inhibitory cyclic nucleotide pathways (cGMP, cAMP) and is supported by experimental data.

A systems-biology model of the tumor necrosis factor (TNF) interactions with TNF receptor 1 and 2.

MOTIVATION: Clustering enables TNF receptors to stimulate intracellular signaling. The differential soluble ligand-induced clustering behavior of TNF receptor 1 (TNFR1) and TNFR2 was modeled. A structured, rule-based model implemented ligand-independent pre-ligand binding assembly domain (PLAD)-mediated homotypic low affinity interactions of unliganded and liganded TNF receptors. RESULTS: Soluble TNF initiates TNFR1 signaling but not TNFR2 signaling despite receptor binding unless it is secondarily oligomerized. We consider high affinity binding of TNF to signaling-incompetent pre-assembled dimeric TNFR1 and TNFR2 molecules and secondary clustering of liganded dimers to signaling competent ligand-receptor clusters. Published receptor numbers, affinities and measured different activities of clustered receptors validated model simulations for a large range of receptor and ligand concentrations. Different PLAD-PLAD affinities and different activities of receptor clusters explain the observed differences in the TNF receptor stimulating activities of soluble TNF. AVAILABILITY AND IMPLEMENTATION: All scripts and data are in manuscript and supplement at Bioinformatics online. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Aluminium toxicokinetics after intramuscular, subcutaneous, and intravenous injection of Al citrate solution in rats.

Knowledge of dose linearity, plasma clearance, rate and extent of subcutaneous (SC) and intramuscular (IM) absorption of soluble aluminium (Al) citrate is considered a prerequisite for evaluation of toxicokinetic data obtained from SC or IM administration of Al adjuvants in medicinal products. Therefore, total Al plasma kinetics was investigated after SC, IM, and IV administration of single Al doses (36 and 360 µg/kg IM or SC; 30 and 300 µg/kg IV) given as citrate solution in rats. Control groups receiving vehicle (saline) were run in parallel to monitor background plasma Al levels over time resulting from dietary intake. Evaluation of Al plasma profiles was done by both non-compartmental analysis of baseline-corrected data and simultaneous model fitting to the raw data using a population kinetics approach. High and dose-independent total plasma clearance (6.6 mL/min/kg) was observed after IV administration corresponding to 60-82% of normal rat GFR. This supports the previous assumptions that parenterally administered Al citrate is more rapidly cleared from plasma than other Al species (e.g., chloride or lactate). Furthermore, plasma exposure of Al (Cmax and AUC0-inf) increased dose-proportionally at all administration routes. Fast and complete absorption of Al was observed at each dose level after both SC and IM administration (bioavailability estimates: 88 and 110%). Estimates for the first-order absorption rate constant ka correspond to absorption half-lives of 36 min (SC) and ≤ 13 min (IM). There was no increase in tissue Al content (whole bone and brain) after 36 µg/kg IM compared to control rats.

Quantitative single-molecule localization microscopy combined with rule-based modeling reveals ligand-induced TNF-R1 reorganization toward higher-order oligomers.

We report on the assembly of tumor necrosis factor receptor 1 (TNF-R1) prior to ligand activation and its ligand-induced reorganization at the cell membrane. We apply single-molecule localization microscopy to obtain quantitative information on receptor cluster sizes and copy numbers. Our data suggest a dimeric pre-assembly of TNF-R1, as well as receptor reorganization toward higher oligomeric states with stable populations comprising three to six TNF-R1. Our experimental results directly serve as input parameters for computational modeling of the ligand-receptor interaction. Simulations corroborate the experimental finding of higher-order oligomeric states. This work is a first demonstration how quantitative, super-resolution and advanced microscopy can be used for systems biology approaches at the single-molecule and single-cell level.

Establishment of a human 3D lung cancer model based on a biological tissue matrix combined with a Boolean in silico model.

For the development of new treatment strategies against cancer, understanding signaling networks and their changes upon drug response is a promising approach to identify new drug targets and biomarker profiles. Pre-requisites are tumor models with multiple read-out options that accurately reflect the clinical situation. Tissue engineering technologies offer the integration of components of the tumor microenvironment which are known to impair drug response of cancer cells. We established three-dimensional (3D) lung carcinoma models on a decellularized tissue matrix, providing a complex microenvironment for cell growth. For model generation, we used two cell lines with (HCC827) or without (A549) an activating mutation of the epidermal growth factor receptor (EGFR), exhibiting different sensitivities to the EGFR inhibitor gefitinib. EGFR activation in HCC827 was inhibited by gefitinib, resulting in a significant reduction of proliferation (Ki-67 proliferation index) and in the induction of apoptosis (TUNEL staining, M30-ELISA). No significant effect was observed in conventional cell culture. Results from the 3D model correlated with the results of an in silico model that integrates the EGFR signaling network according to clinical data. The application of TGFß1 induced tumor cell invasion, accompanied by epithelial-mesenchymal transition (EMT) both in vitro and in silico. This was confirmed in the 3D model by acquisition of mesenchymal cell morphology and modified expression of fibronectin, E-cadherin, ß-catenin and mucin-1. Quantitative read-outs for proliferation, apoptosis and invasion were established in the complex 3D tumor model. The combined in vitro and in silico model represents a powerful tool for systems analysis.

Exogenous administration of a recombinant variant of TWEAK impairs healing after myocardial infarction by aggravation of inflammation.

BACKGROUND: Tumor necrosis factor-like weak inducer of apoptosis (TWEAK) and its receptor fibroblast growth factor-inducible 14 (Fn14) are upregulated after myocardial infarction (MI) in both humans and mice. They modulate inflammation and the extracellular matrix, and could therefore be important for healing and remodeling after MI. However, the function of TWEAK after MI remains poorly defined. METHODS AND RESULTS: Following ligation of the left coronary artery, mice were injected twice per week with a recombinant human serum albumin conjugated variant of TWEAK (HSA-Flag-TWEAK), mimicking the activity of soluble TWEAK. Treatment with HSA-Flag-TWEAK resulted in significantly increased mortality in comparison to the placebo group due to myocardial rupture. Infarct size, extracellular matrix remodeling, and apoptosis rates were not different after MI. However, HSA-Flag-TWEAK treatment increased infiltration of proinflammatory cells into the myocardium. Accordingly, depletion of neutrophils prevented cardiac ruptures without modulating all-cause mortality. CONCLUSION: Treatment of mice with HSA-Flag-TWEAK induces myocardial healing defects after experimental MI. This is mediated by an exaggerated neutrophil infiltration into the myocardium.

A Boolean view separates platelet activatory and inhibitory signalling as verified by phosphorylation monitoring including threshold behaviour and integrin modulation.

Platelets are critical for haemostasis and blood clotting. However, since under normal circumstances blood should flow without clotting, its function is regulated via a complex interplay of activating and inhibiting signal transduction pathways. Understanding this network is crucial for treatment of cardiovascular and bleeding diseases. Detailed protein interaction and phosphorylation data are explored to establish a simplified Boolean model of the central platelet cascades. We implemented the model by means of CellNetAnalyzer and showed how different signalling events coalesce into a fully activated system state. Furthermore, we examined the networks' inherent threshold behaviour using the semi-quantitative modelling software SQUAD. Finally, predictions are verified monitoring phosphorylations which mark different activation phases as modelled. The model can also be applied to simulate different pharmacological conditions as they modify node activity (aspirin, clopidogrel, milrinon, iloprost, combination) and is available for further studies. It agrees well with observations. Activatory pathways are diversified to cope with complex environmental conditions. Platelet activation needs several activation steps to integrate over different network subsets, as they are formed by the interplay of activating kinases, calcium mobilization, and the inhibiting cAMP-PKA system. System stability analysis shows two phases: a sub-threshold behaviour, characterized by integration over different activatory and inhibitory conditions, and a beyond threshold phase, represented by competition and shutting down of counter-regulatory pathways. The integrin network and Akt-protein are critical for stable effector response. Dynamic threshold-analysis reveals a dependency of the relative activating input strength necessary to irreversibly engage the system from the absolute inhibitory signal strength.

Time-resolved in silico modeling of fine-tuned cAMP signaling in platelets: feedback loops, titrated phosphorylations and pharmacological modulation.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.