TY - JOUR A1 - Wangorsch, Gaby A1 - Butt, Elke A1 - Mark, Regina A1 - Hubertus, Katharina A1 - Geiger, Jörg A1 - Dandekar, Thomas A1 - Dittrich, Marcus T1 - Time-resolved in silico modeling of fine-tuned cAMP signaling in platelets: feedback loops, titrated phosphorylations and pharmacological modulation N2 - 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 antiplatelet drugs and their combinations revealed a high fidelity, fine-tuned cAMP signaling in platelets without crosstalk 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. KW - Vasodilatator-stimuliertes Phosphoprotein KW - VASP KW - cyclic nucleotide signaling KW - silico model Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-69145 ER - TY - JOUR A1 - Beisser, Daniela A1 - Grohme, Markus A. A1 - Kopka, Joachim A1 - Frohme, Marcus A1 - Schill, Ralph O. A1 - Hengherr, Steffen A1 - Dandekar, Thomas A1 - Klau, Gunnar W. A1 - Dittrich, Marcus A1 - Müller, Tobias T1 - Integrated pathway modules using time-course metabolic profiles and EST data from Milnesium tardigradum N2 - Background: Tardigrades are multicellular organisms, resistant to extreme environmental changes such as heat, drought, radiation and freezing. They outlast these conditions in an inactive form (tun) to escape damage to cellular structures and cell death. Tardigrades are apparently able to prevent or repair such damage and are therefore a crucial model organism for stress tolerance. Cultures of the tardigrade Milnesium tardigradum were dehydrated by removing the surrounding water to induce tun formation. During this process and the subsequent rehydration, metabolites were measured in a time series by GC-MS. Additionally expressed sequence tags are available, especially libraries generated from the active and inactive state. The aim of this integrated analysis is to trace changes in tardigrade metabolism and identify pathways responsible for their extreme resistance against physical stress. Results: In this study we propose a novel integrative approach for the analysis of metabolic networks to identify modules of joint shifts on the transcriptomic and metabolic levels. We derive a tardigrade-specific metabolic network represented as an undirected graph with 3,658 nodes (metabolites) and 4,378 edges (reactions). Time course metabolite profiles are used to score the network nodes showing a significant change over time. The edges are scored according to information on enzymes from the EST data. Using this combined information, we identify a key subnetwork (functional module) of concerted changes in metabolic pathways, specific for de- and rehydration. The module is enriched in reactions showing significant changes in metabolite levels and enzyme abundance during the transition. It resembles the cessation of a measurablemetabolism (e.g. glycolysis and amino acid anabolism) during the tun formation, the production of storage metabolites and bioprotectants, such as DNA stabilizers, and the generation of amino acids and cellular components from monosaccharides as carbon and energy source during rehydration. Conclusions: The functional module identifies relationships among changed metabolites (e.g. spermidine) and reactions and provides first insights into important altered metabolic pathways. With sparse and diverse data available, the presented integrated metabolite network approach is suitable to integrate all existing data and analyse it in a combined manner. KW - Milnesium tardigradum KW - Integrated network analysis KW - Functional modules KW - Metabolic profiles KW - Metabolic pathways KW - Trend test Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-75241 ER - TY - JOUR A1 - Staiger, Christine A1 - Cadot, Sidney A1 - Kooter, Raul A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Klau, Gunnar W. A1 - Wessels, Lodewyk F. A. T1 - A Critical Evaluation of Network and Pathway-Based Classifiers for Outcome Prediction in Breast Cancer JF - PLoS One N2 - Recently, several classifiers that combine primary tumor data, like gene expression data, and secondary data sources, such as protein-protein interaction networks, have been proposed for predicting outcome in breast cancer. In these approaches, new composite features are typically constructed by aggregating the expression levels of several genes. The secondary data sources are employed to guide this aggregation. Although many studies claim that these approaches improve classification performance over single genes classifiers, the gain in performance is difficult to assess. This stems mainly from the fact that different breast cancer data sets and validation procedures are employed to assess the performance. Here we address these issues by employing a large cohort of six breast cancer data sets as benchmark set and by performing an unbiased evaluation of the classification accuracies of the different approaches. Contrary to previous claims, we find that composite feature classifiers do not outperform simple single genes classifiers. We investigate the effect of (1) the number of selected features; (2) the specific gene set from which features are selected; (3) the size of the training set and (4) the heterogeneity of the data set on the performance of composite feature and single genes classifiers. Strikingly, we find that randomization of secondary data sources, which destroys all biological information in these sources, does not result in a deterioration in performance of composite feature classifiers. Finally, we show that when a proper correction for gene set size is performed, the stability of single genes sets is similar to the stability of composite feature sets. Based on these results there is currently no reason to prefer prognostic classifiers based on composite features over single genes classifiers for predicting outcome in breast cancer. KW - modules KW - protein-interaction networks KW - expression signature KW - classification KW - set KW - metastasis KW - stability KW - survival KW - database KW - markers Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-131323 VL - 7 IS - 4 ER - TY - JOUR A1 - Floren, Andreas A1 - Mupepele, Anne-Christine A1 - Müller, Tobias A1 - Dittrich, Marcus T1 - Are Temperate Canopy Spiders Tree-Species Specific? N2 - Arboreal spiders in deciduous and coniferous trees were investigated on their distribution and diversity. Insecticidal knock-down was used to comprehensively sample spiders from 175 trees from 2001 to 2003 in the Białowieża forest and three remote forests in Poland. We identified 140 species from 9273 adult spiders. Spider communities were distinguished between deciduous and coniferous trees. The richest fauna was collected from Quercus where beta diversity was also highest. A tree-species-specific pattern was clearly observed for Alnus, Carpinus, Picea and Pinus trees and also for those tree species that were fogged in only four or three replicates, namely Betula and Populus. This hitherto unrecognised association was mainly due to the community composition of common species identified in a Dufrene-Legendre indicator species analysis. It was not caused by spatial or temporal autocorrelation. Explaining tree-species specificity for generalist predators like spiders is difficult and has to involve physical and ecological tree parameters like linkage with the abundance of prey species. However, neither did we find a consistent correlation of prey group abundances with spiders nor could differences in spider guild composition explain the observed pattern. Our results hint towards the importance of deterministic mechanisms structuring communities of generalist canopy spiders although the casual relationship is not yet understood. KW - trees KW - spiders KW - conifers KW - forests KW - predation KW - oaks KW - community structures KW - pines Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-111413 ER - TY - JOUR A1 - Morton, Charles Oliver A1 - Fliesser, Mirjam A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Bauer, Ruth A1 - Kneitz, Susanne A1 - Hope, William A1 - Rogers, Thomas Richard A1 - Einsele, Hermann A1 - Löffler, Jürgen T1 - Gene Expression Profiles of Human Dendritic Cells Interacting with Aspergillus fumigatus in a Bilayer Model of the Alveolar Epithelium/Endothelium Interface N2 - The initial stages of the interaction between the host and Aspergillus fumigatus at the alveolar surface of the human lung are critical in the establishment of aspergillosis. Using an in vitro bilayer model of the alveolus, including both the epithelium (human lung adenocarcinoma epithelial cell line, A549) and endothelium (human pulmonary artery epithelial cells, HPAEC) on transwell membranes, it was possible to closely replicate the in vivo conditions. Two distinct sub-groups of dendritic cells (DC), monocyte-derived DC (moDC) and myeloid DC (mDC), were included in the model to examine immune responses to fungal infection at the alveolar surface. RNA in high quantity and quality was extracted from the cell layers on the transwell membrane to allow gene expression analysis using tailored custom-made microarrays, containing probes for 117 immune-relevant genes. This microarray data indicated minimal induction of immune gene expression in A549 alveolar epithelial cells in response to germ tubes of A. fumigatus. In contrast, the addition of DC to the system greatly increased the number of differentially expressed immune genes. moDC exhibited increased expression of genes including CLEC7A, CD209 and CCL18 in the absence of A. fumigatus compared to mDC. In the presence of A. fumigatus, both DC subgroups exhibited up-regulation of genes identified in previous studies as being associated with the exposure of DC to A. fumigatus and exhibiting chemotactic properties for neutrophils, including CXCL2, CXCL5, CCL20, and IL1B. This model closely approximated the human alveolus allowing for an analysis of the host pathogen interface that complements existing animal models of IA. KW - aspergillus fumigatus KW - gene expression KW - immune receptors KW - immune response KW - denritic cells KW - B cell receptors KW - gene regulation KW - RNA extraction Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-112893 ER - TY - JOUR A1 - Karulin, Alexey Y. A1 - Caspell, Richard A1 - Dittrich, Marcus A1 - Lehmann, Paul V. T1 - Normal distribution of CD8+ T-cell-derived ELISPOT counts within replicates justifies the reliance on parametric statistics for identifying positive responses JF - Cells N2 - Accurate assessment of positive ELISPOT responses for low frequencies of antigen-specific T-cells is controversial. In particular, it is still unknown whether ELISPOT counts within replicate wells follow a theoretical distribution function, and thus whether high power parametric statistics can be used to discriminate between positive and negative wells. We studied experimental distributions of spot counts for up to 120 replicate wells of IFN-γ production by CD8+ T-cell responding to EBV LMP2A (426 – 434) peptide in human PBMC. The cells were tested in serial dilutions covering a wide range of average spot counts per condition, from just a few to hundreds of spots per well. Statistical analysis of the data using diagnostic Q-Q plots and the Shapiro-Wilk normality test showed that in the entire dynamic range of ELISPOT spot counts within replicate wells followed a normal distribution. This result implies that the Student t-Test and ANOVA are suited to identify positive responses. We also show experimentally that borderline responses can be reliably detected by involving more replicate wells, plating higher numbers of PBMC, addition of IL-7, or a combination of these. Furthermore, we have experimentally verified that the number of replicates needed for detection of weak responses can be calculated using parametric statistics. KW - ELISPOT KW - statistics KW - t-Test KW - ANOVA KW - T-cells KW - normal distribution Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-149968 VL - 4 IS - 1 ER - TY - JOUR A1 - Kuhn, Joachim A1 - Gripp, Tatjana A1 - Flieder, Tobias A1 - Dittrich, Marcus A1 - Hendig, Doris A1 - Busse, Jessica A1 - Knabbe, Cornelius A1 - Birschmann, Ingvild T1 - UPLC-MRM Mass Spectrometry Method for Measurement of the Coagulation Inhibitors Dabigatran and Rivaroxaban in Human Plasma and Its Comparison with Functional Assays JF - PLOS ONE N2 - Introduction The fast, precise, and accurate measurement of the new generation of oral anticoagulants such as dabigatran and rivaroxaban in patients' plasma my provide important information in different clinical circumstances such as in the case of suspicion of overdose, when patients switch from existing oral anticoagulant, in patients with hepatic or renal impairment, by concomitant use of interaction drugs, or to assess anticoagulant concentration in patients' blood before major surgery. Methods Here, we describe a quick and precise method to measure the coagulation inhibitors dabigatran and rivaroxaban using ultra-performance liquid chromatography electrospray ionization-tandem mass spectrometry in multiple reactions monitoring (MRM) mode (UPLC-MRM MS). Internal standards (ISs) were added to the sample and after protein precipitation; the sample was separated on a reverse phase column. After ionization of the analytes the ions were detected using electrospray ionization-tandem mass spectrometry. Run time was 2.5 minutes per injection. Ion suppression was characterized by means of post-column infusion. Results The calibration curves of dabigatran and rivaroxaban were linear over the working range between 0.8 and 800 mu g/L (r > 0.99). Limits of detection (LOD) in the plasma matrix were 0.21 mu g/L for dabigatran and 0.34 mu g/L for rivaroxaban, and lower limits of quantification (LLOQ) in the plasma matrix were 0.46 mu g/L for dabigatran and 0.54 mu g/L for rivaroxaban. The intraassay coefficients of variation (CVs) for dabigatran and rivaroxaban were < 4% and 6%; respectively, the interassay CVs were < 6% for dabigatran and < 9% for rivaroxaban. Inaccuracy was < 5% for both substances. The mean recovery was 104.5% (range 83.8-113.0%) for dabigatran and 87.0%(range 73.6-105.4%) for rivaroxaban. No significant ion suppressions were detected at the elution times of dabigatran or rivaroxaban. Both coagulation inhibitors were stable in citrate plasma at -20 degrees C, 4 degrees C and even at RT for at least one week. A method comparison between our UPLC-MRM MS method, the commercially available automated Direct Thrombin Inhibitor assay (DTI assay) for dabigatran measurement from CoaChrom Diagnostica, as well as the automated anti-Xa assay for rivaroxaban measurement from Chromogenix both performed by ACL-TOP showed a high degree of correlation. However, UPLC-MRM MS measurement of dabigatran and rivaroxaban has a much better selectivity than classical functional assays measuring activities of various coagulation factors which are susceptible to interference by other coagulant drugs. Conclusions Overall, we developed and validated a sensitive and specific UPLC-MRM MS assay for the quick and specific measurement of dabigatran and rivaroxaban in human plasma. KW - LC-MS/MS KW - validation KW - serum KW - quantification KW - apixaban KW - diagnostic accuracy KW - performance liquid chromatography KW - factor XA inhibitor KW - direct oral anticoagulants KW - direct thrombin inhibitor Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-136023 VL - 10 IS - 12 ER - TY - JOUR A1 - Floren, Andreas A1 - Krüger, Dirk A1 - Müller, Tobias A1 - Dittrich, Marcus A1 - Rudloff, Renate A1 - Hoppe, Björn A1 - Linsenmair, Karl Eduard T1 - Diversity and interactions of wood-inhabiting fungi and beetles after deadwood enrichment JF - PLoS ONE N2 - Freshly cut beech deadwood was enriched in the canopy and on the ground in three cultural landscapes in Germany (Swabian Alb, Hainich-Dun, Schorfheide-Chorin) in order to analyse the diversity, distribution and interaction of wood-inhabiting fungi and beetles. After two years of wood decay 83 MOTUs (Molecular Operational Taxonomic Units) from 28 wood samples were identified. Flight Interception Traps (FITs) installed adjacent to the deadwood enrichments captured 29.465 beetles which were sorted to 566 species. Geographical 'region' was the main factor determining both beetle and fungal assemblages. The proportions of species occurring in all regions were low. Statistic models suggest that assemblages of both taxa differed between stratum and management praxis but their strength varied among regions. Fungal assemblages in Hainich-Dun, for which the data was most comprehensive, discriminated unmanaged from extensively managed and age-class forests (even-aged timber management) while canopy communities differed not from those near the ground. In contrast, the beetle assemblages at the same sites showed the opposite pattern. We pursued an approach in the search for fungus-beetle associations by computing cross correlations and visualize significant links in a network graph. These correlations can be used to formulate hypotheses on mutualistic relationships for example in respect to beetles acting as vectors of fungal spores. KW - european beech forests KW - bark beetles KW - management KW - decay KW - ecology KW - norway spruce KW - substrate quality KW - communities KW - rare Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-145129 VL - 10 IS - 11 ER - TY - JOUR A1 - Karulin, Alexey Y. A1 - Karacsony, Kinga A1 - Zhang, Wenji A1 - Targoni, Oleg S. A1 - Moldova, Ioana A1 - Dittrich, Marcus A1 - Sundararaman, Srividya A1 - Lehmann, Paul V. T1 - ELISPOTs produced by CD8 and CD4 cells follow Log Normal size distribution permitting objective counting JF - Cells N2 - Each positive well in ELISPOT assays contains spots of variable sizes that can range from tens of micrometers up to a millimeter in diameter. Therefore, when it comes to counting these spots the decision on setting the lower and the upper spot size thresholds to discriminate between non-specific background noise, spots produced by individual T cells, and spots formed by T cell clusters is critical. If the spot sizes follow a known statistical distribution, precise predictions on minimal and maximal spot sizes, belonging to a given T cell population, can be made. We studied the size distributional properties of IFN-γ, IL-2, IL-4, IL-5 and IL-17 spots elicited in ELISPOT assays with PBMC from 172 healthy donors, upon stimulation with 32 individual viral peptides representing defined HLA Class I-restricted epitopes for CD8 cells, and with protein antigens of CMV and EBV activating CD4 cells. A total of 334 CD8 and 80 CD4 positive T cell responses were analyzed. In 99.7% of the test cases, spot size distributions followed Log Normal function. These data formally demonstrate that it is possible to establish objective, statistically validated parameters for counting T cell ELISPOTs. KW - ELISPOT KW - software KW - IFN-γ KW - IL-17 KW - T cells KW - Normal Distribution KW - spot size KW - gating KW - cytokines KW - IL-2 KW - IL-4 KW - IL-5 KW - CD8 KW - CD4 Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-149648 VL - 4 IS - 1 ER - TY - JOUR A1 - Appel, Mirjam A1 - Scholz, Claus-Jürgen A1 - Müller, Tobias A1 - Dittrich, Marcus A1 - König, Christian A1 - Bockstaller, Marie A1 - Oguz, Tuba A1 - Khalili, Afshin A1 - Antwi-Adjei, Emmanuel A1 - Schauer, Tamas A1 - Margulies, Carla A1 - Tanimoto, Hiromu A1 - Yarali, Ayse T1 - Genome-Wide Association Analyses Point to Candidate Genes for Electric Shock Avoidance in Drosophila melanogaster JF - PLoS ONE N2 - Electric shock is a common stimulus for nociception-research and the most widely used reinforcement in aversive associative learning experiments. Yet, nothing is known about the mechanisms it recruits at the periphery. To help fill this gap, we undertook a genome-wide association analysis using 38 inbred Drosophila melanogaster strains, which avoided shock to varying extents. We identified 514 genes whose expression levels and/or sequences covaried with shock avoidance scores. We independently scrutinized 14 of these genes using mutants, validating the effect of 7 of them on shock avoidance. This emphasizes the value of our candidate gene list as a guide for follow-up research. In addition, by integrating our association results with external protein-protein interaction data we obtained a shock avoidance- associated network of 38 genes. Both this network and the original candidate list contained a substantial number of genes that affect mechanosensory bristles, which are hairlike organs distributed across the fly's body. These results may point to a potential role for mechanosensory bristles in shock sensation. Thus, we not only provide a first list of candidate genes for shock avoidance, but also point to an interesting new hypothesis on nociceptive mechanisms. KW - functional analysis KW - disruption project KW - natural variation KW - complex traits KW - networks KW - behavior KW - flies KW - temperature KW - genetics KW - painful Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-152006 VL - 10 IS - 5 ER -