570 Biowissenschaften; Biologie
Refine
Has Fulltext
- yes (147)
Is part of the Bibliography
- yes (147)
Year of publication
- 2019 (147) (remove)
Document Type
- Journal article (86)
- Doctoral Thesis (60)
- Preprint (1)
Keywords
- Tissue Engineering (7)
- cancer (4)
- Candida albicans (3)
- Genexpression (3)
- In vitro (3)
- Signaltransduktion (3)
- leukemic cells (3)
- metagenomics (3)
- microbiome (3)
- 3D Tumormodell (2)
Institute
- Theodor-Boveri-Institut für Biowissenschaften (74)
- Graduate School of Life Sciences (33)
- Julius-von-Sachs-Institut für Biowissenschaften (10)
- Rudolf-Virchow-Zentrum (10)
- Lehrstuhl für Tissue Engineering und Regenerative Medizin (9)
- Fakultät für Biologie (4)
- Institut für Hygiene und Mikrobiologie (4)
- Institut für Molekulare Infektionsbiologie (4)
- Institut für Pharmakologie und Toxikologie (4)
- Medizinische Klinik und Poliklinik II (4)
- Institut für Virologie und Immunbiologie (3)
- Lehrstuhl für Biochemie (3)
- Abteilung für Molekulare Innere Medizin (in der Medizinischen Klinik und Poliklinik II) (2)
- Center for Computational and Theoretical Biology (2)
- Deutsches Zentrum für Herzinsuffizienz (DZHI) (2)
- Institut für Anatomie und Zellbiologie (2)
- Institut für Humangenetik (2)
- Institut für Medizinische Strahlenkunde und Zellforschung (2)
- Klinik und Poliklinik für Psychiatrie, Psychosomatik und Psychotherapie (2)
- Pathologisches Institut (2)
- Abteilung für Funktionswerkstoffe der Medizin und der Zahnheilkunde (1)
- Comprehensive Cancer Center Mainfranken (1)
- Institut für Experimentelle Biomedizin (1)
- Institut für Geographie und Geologie (1)
- Institut für Klinische Epidemiologie und Biometrie (1)
- Institut für Klinische Neurobiologie (1)
- Institut für Pharmazie und Lebensmittelchemie (1)
- Institut für Systemimmunologie (1)
- Klinik und Poliklinik für Dermatologie, Venerologie und Allergologie (1)
- Klinik und Poliklinik für Kinder- und Jugendpsychiatrie, Psychosomatik und Psychotherapie (1)
- Klinik und Poliklinik für Nuklearmedizin (1)
- Lehrstuhl für Molekulare Psychiatrie (1)
- Medizinische Klinik und Poliklinik I (1)
- Urologische Klinik und Poliklinik (1)
Sonstige beteiligte Institutionen
Metagenomic sequencing has greatly improved our ability to profile the composition of environmental and host-associated microbial communities. However, the dependency of most methods on reference genomes, which are currently unavailable for a substantial fraction of microbial species, introduces estimation biases. We present an updated and functionally extended tool based on universal (i.e., reference-independent), phylogenetic marker gene (MG)-based operational taxonomic units (mOTUs) enabling the profiling of >7700 microbial species. As more than 30% of them could not previously be quantified at this taxonomic resolution, relative abundance estimates based on mOTUs are more accurate compared to other methods. As a new feature, we show that mOTUs, which are based on essential housekeeping genes, are demonstrably well-suited for quantification of basal transcriptional activity of community members. Furthermore, single nucleotide variation profiles estimated using mOTUs reflect those from whole genomes, which allows for comparing microbial strain populations (e.g., across different human body sites).
Thermal melanism theory states that dark-colored ectotherm organisms are at an advantage at low temperature due to increased warming. This theory is generally supported for ectotherm animals, however, the function of colors in the fungal kingdom is largely unknown. Here, we test whether the color lightness of mushroom assemblages is related to climate using a dataset of 3.2 million observations of 3,054 species across Europe. Consistent with the thermal melanism theory, mushroom assemblages are significantly darker in areas with cold climates. We further show differences in color phenotype between fungal lifestyles and a lifestyle differentiated response to seasonality. These results indicate a more complex ecological role of mushroom colors and suggest functions beyond thermal adaption. Because fungi play a crucial role in terrestrial carbon and nutrient cycles, understanding the links between the thermal environment, functional coloration and species’ geographical distributions will be critical in predicting ecosystem responses to global warming.
How insects promote crop pollination remains poorly understood in terms of the contribution of functional trait differences between species. We used meta-analyses to test for correlations between community abundance, species richness and functional trait metrics with oilseed rape yield, a globally important crop. While overall abundance is consistently important in predicting yield, functional divergence between species traits also showed a positive correlation. This result supports the complementarity hypothesis that pollination function is maintained by non-overlapping trait distributions. In artificially constructed communities (mesocosms), species richness is positively correlated with yield, although this effect is not seen under field conditions. As traits of the dominant species do not predict yield above that attributed to the effect of abundance alone, we find no evidence in support of the mass ratio hypothesis. Management practices increasing not just pollinator abundance, but also functional divergence, could benefit oilseed rape agriculture.
The cytoskeletal crosslinking protein MACF1 is dispensable for thrombus formation and hemostasis
(2019)
Coordinated reorganization of cytoskeletal structures is critical for key aspects of platelet physiology. While several studies have addressed the role of microtubules and filamentous actin in platelet production and function, the significance of their crosstalk in these processes has been poorly investigated. The microtubule-actin cross-linking factor 1 (MACF1; synonym: Actin cross-linking factor 7, ACF7) is a member of the spectraplakin family, and one of the few proteins expressed in platelets, which possess actin and microtubule binding domains thereby facilitating actin-microtubule interaction and regulation. We used megakaryocyte- and platelet-specific Macf1 knockout (Macf1fl/fl, Pf4-Cre) mice to study the role of MACF1 in platelet production and function. MACF1 deficient mice displayed comparable platelet counts to control mice. Analysis of the platelet cytoskeletal ultrastructure revealed a normal marginal band and actin network. Platelet spreading on fibrinogen was slightly delayed but platelet activation and clot traction was unaffected. Ex vivo thrombus formation and mouse tail bleeding responses were similar between control and mutant mice. These results suggest that MACF1 is dispensable for thrombopoiesis, platelet activation, thrombus formation and the hemostatic function in mice.
Immunotherapy with chimeric antigen receptor-engineered T-cells (CAR-T) is under investigation in multiple myeloma. There are reports of myeloma remission after CD19 CAR-T therapy, although CD19 is hardly detectable on myeloma cells by flow cytometry (FC). We apply single molecule-sensitive direct stochastic optical reconstruction microscopy (dSTORM), and demonstrate CD19 expression on a fraction of myeloma cells (10.3–80%) in 10 out of 14 patients (density: 13–5,000 molecules per cell). In contrast, FC detects CD19 in only 2 of these 10 patients, on a smaller fraction of cells. Treatment with CD19 CAR-T in vitro results in elimination of CD19-positive myeloma cells, including those with <100 CD19 molecules per cell. Similar data are obtained by dSTORM analyses of CD20 expression on myeloma cells and CD20 CAR-T. These data establish a sensitivity threshold for CAR-T and illustrate how super-resolution microscopy can guide patient selection in immunotherapy to exploit ultra-low density antigens.
Hsp90 is a dimeric molecular chaperone that is essential for the folding and activation of hundreds of client proteins. Co-chaperone proteins regulate the ATP-driven Hsp90 client activation cycle. Aha-type co-chaperones are the most potent stimulators of the Hsp90 ATPase activity but the relationship between ATPase regulation and in vivo activity is poorly understood. We report here that the most strongly conserved region of Aha-type co-chaperones, the N terminal NxNNWHW motif, modulates the apparent affinity of Hsp90 for nucleotide substrates. The ability of yeast Aha-type co-chaperones to act in vivo is ablated when the N terminal NxNNWHW motif is removed. This work suggests that nucleotide exchange during the Hsp90 functional cycle may be more important than rate of catalysis.
Periodontal treatment prevents arthritis in mice and methotrexate ameliorates periodontal bone loss
(2019)
Recent studies indicate a causal relationship between the periodontal pathogen P. gingivalis and rheumatoid arthritis involving the production of autoantibodies against citrullinated peptides. We therefore postulated that therapeutic eradication P. gingivalis may ameliorate rheumatoid arthritis development and here turned to a mouse model in order to challenge our hypothesis. F1 (DBA/1 x B10.Q) mice were orally inoculated with P. gingivalis before collagen-induced arthritis was provoked. Chlorhexidine or metronidazole were orally administered either before or during the induction phase of arthritis and their effects on arthritis progression and alveolar bone loss were compared to intraperitoneally injected methotrexate. Arthritis incidence and severity were macroscopically scored and alveolar bone loss was evaluated via microcomputed tomography. Serum antibody titres against P. gingivalis were quantified by ELISA and microbial dysbiosis following oral inoculation was monitored in stool samples via microbiome analyses. Both, oral chlorhexidine and metronidazole reduced the incidence and ameliorated the severity of collagen-induced arthritis comparable to methotrexate. Likewise, all three therapies attenuated alveolar bone loss. Relative abundance of Porphyromonadaceae was increased after oral inoculation with P. gingivalis and decreased after treatment. This is the first study to describe beneficial effects of non-surgical periodontal treatment on collagen-induced arthritis in mice and suggests that mouthwash with chlorhexidine or metronidazole may also be beneficial for patients with rheumatoid arthritis and a coexisting periodontitis. Methotrexate ameliorated periodontitis in mice, further raising the possibility that methotrexate may also positively impact on the tooth supporting tissues of patients with rheumatoid arthritis.
Enormous efforts have been made to target metabolic dependencies of cancer cells for developing new therapies. However, the therapeutic efficacy of glycolysis inhibitors is limited due to their inability to elicit cell death. Hexokinase 2 (HK2), via its mitochondrial localization, functions as a central nexus integrating glycolysis activation and apoptosis resilience. Here we identify that K63-linked ubiquitination by HectH9 regulates the mitochondrial localization and function of HK2. Through stable isotope tracer approach and functional metabolic analyses, we show that HectH9 deficiency impedes tumor glucose metabolism and growth by HK2 inhibition. The HectH9/HK2 pathway regulates cancer stem cell (CSC) expansion and CSC-associated chemoresistance. Histological analyses show that HectH9 expression is upregulated and correlated with disease progression in prostate cancer. This work uncovers that HectH9 is a novel regulator of HK2 and cancer metabolism. Targeting HectH9 represents an effective strategy to achieve long-term tumor remission by concomitantly disrupting glycolysis and inducing apoptosis.
In most organisms, ribosomal RNA (rRNA) contributes to >85% of total RNA. Thus, to obtain useful information from RNA-sequencing (RNA-seq) analyses at reasonable sequencing depth, typically, mature polyadenylated transcripts are enriched or rRNA molecules are depleted. Targeted depletion of rRNA is particularly useful when studying transcripts lacking a poly(A) tail, such as some non-coding RNAs (ncRNAs), most bacterial RNAs and partially degraded or immature transcripts. While several commercially available kits allow effective rRNA depletion, their efficiency relies on a high degree of sequence homology between oligonucleotide probes and the target RNA. This restricts the use of such kits to a limited number of organisms with conserved rRNA sequences. In this study we describe the use of biotinylated oligos and streptavidin-coated paramagnetic beads for the efficient and specific depletion of trypanosomal rRNA. Our approach reduces the levels of the most abundant rRNA transcripts to less than 5% with minimal off-target effects. By adjusting the sequence of the oligonucleotide probes, our approach can be used to deplete rRNAs or other abundant transcripts independent of species. Thus, our protocol provides a useful alternative for rRNA removal where enrichment of polyadenylated transcripts is not an option and commercial kits for rRNA are not available.
Icefishes (suborder Notothenioidei; family Channichthyidae) are the only vertebrates that lack functional haemoglobin genes and red blood cells. Here, we report a high-quality genome assembly and linkage map for the Antarctic blackfin icefish Chaenocephalus aceratus, highlighting evolved genomic features for its unique physiology. Phylogenomic analysis revealed that Antarctic fish of the teleost suborder Notothenioidei, including icefishes, diverged from the stickleback lineage about 77 million years ago and subsequently evolved cold-adapted phenotypes as the Southern Ocean cooled to sub-zero temperatures. Our results show that genes involved in protection from ice damage, including genes encoding antifreeze glycoprotein and zona pellucida proteins, are highly expanded in the icefish genome. Furthermore, genes that encode enzymes that help to control cellular redox state, including members of the sod3 and nqo1 gene families, are expanded, probably as evolutionary adaptations to the relatively high concentration of oxygen dissolved in cold Antarctic waters. In contrast, some crucial regulators of circadian homeostasis (cry and per genes) are absent from the icefish genome, suggesting compromised control of biological rhythms in the polar light environment. The availability of the icefish genome sequence will accelerate our understanding of adaptation to extreme Antarctic environments.
Dishevelled (DVL) is the key component of the Wnt signaling pathway. Currently, DVL conformational dynamics under native conditions is unknown. To overcome this limitation, we develop the Fluorescein Arsenical Hairpin Binder- (FlAsH-) based FRET in vivo approach to study DVL conformation in living cells. Using this single-cell FRET approach, we demonstrate that (i) Wnt ligands induce open DVL conformation, (ii) DVL variants that are predominantly open, show more even subcellular localization and more efficient membrane recruitment by Frizzled (FZD) and (iii) Casein kinase 1 ɛ (CK1ɛ) has a key regulatory function in DVL conformational dynamics. In silico modeling and in vitro biophysical methods explain how CK1ɛ-specific phosphorylation events control DVL conformations via modulation of the PDZ domain and its interaction with DVL C-terminus. In summary, our study describes an experimental tool for DVL conformational sampling in living cells and elucidates the essential regulatory role of CK1ɛ in DVL conformational dynamics.
Zinc (Zn2+) can modulate platelet and coagulation activation pathways, including fibrin formation. Here, we studied the (patho)physiological consequences of abnormal platelet Zn2+ storage and release. To visualize Zn2+ storage in human and mouse platelets, the Zn2+ specific fluorescent dye FluoZin3 was used. In resting platelets, the dye transiently accumulated into distinct cytosolic puncta, which were lost upon platelet activation. Platelets isolated from Unc13d−/− mice, characterized by combined defects of α/δ granular release, showed a markedly impaired Zn2+ release upon activation. Platelets from Nbeal2−/− mice mimicking Gray platelet syndrome (GPS), characterized by primarily loss of the α-granule content, had strongly reduced Zn2+ levels, which was also confirmed in primary megakaryocytes. In human platelets isolated from patients with GPS, Hermansky-Pudlak Syndrome (HPS) and Storage Pool Disease (SPD) altered Zn2+ homeostasis was detected. In turbidity and flow based assays, platelet-dependent fibrin formation was impaired in both Nbeal2−/− and Unc13d−/− mice, and the impairment could be partially restored by extracellular Zn2+. Altogether, we conclude that the release of ionic Zn2+ store from secretory granules upon platelet activation contributes to the procoagulant role of Zn2+ in platelet-dependent fibrin formation.
The NEDD8-activating enzyme (NAE) inhibitor MLN4924 inhibits cullin-RING ubiquitin ligase complexes including the SKP1-cullin-F-box E3 ligase βTrCP. MLN4924 therefore inhibits also the βTrCP-dependent activation of the classical and the alternative NFĸB pathway. In this work, we found that a subgroup of multiple myeloma cell lines (e.g., RPMI-8226, MM.1S, KMS-12BM) and about half of the primary myeloma samples tested are sensitized to TNF-induced cell death by MLN4924. This correlated with MLN4924-mediated inhibition of TNF-induced activation of the classical NFκB pathway and reduced the efficacy of TNF-induced TNFR1 signaling complex formation. Interestingly, binding studies revealed a straightforward correlation between cell surface TNFR1 expression in multiple myeloma cell lines and their sensitivity for MLN4924/TNF-induced cell death. The cell surface expression levels of TNFR1 in the investigated MM cell lines largely correlated with TNFR1 mRNA expression. This suggests that the variable levels of cell surface expression of TNFR1 in myeloma cell lines are decisive for TNF/MLN4924 sensitivity. Indeed, introduction of TNFR1 into TNFR1-negative TNF/MLN4924-resistant KMS-11BM cells, was sufficient to sensitize this cell line for TNF/MLN4924-induced cell death. Thus, MLN4924 might be especially effective in myeloma patients with TNFR1+ myeloma cells and a TNFhigh tumor microenvironment.
Fanconi anemia (FA) is a genetically heterogeneous disorder with 22 disease-causing genes reported to date. In some FA genes, monoallelic mutations have been found to be associated with breast cancer risk, while the risk associations of others remain unknown. The gene for FA type C, FANCC, has been proposed as a breast cancer susceptibility gene based on epidemiological and sequencing studies. We used the Oncoarray project to genotype two truncating FANCC variants (p.R185X and p.R548X) in 64,760 breast cancer cases and 49,793 controls of European descent. FANCC mutations were observed in 25 cases (14 with p.R185X, 11 with p.R548X) and 26 controls (18 with p.R185X, 8 with p.R548X). There was no evidence of an association with the risk of breast cancer, neither overall (odds ratio 0.77, 95%CI 0.44–1.33, p = 0.4) nor by histology, hormone receptor status, age or family history. We conclude that the breast cancer risk association of these two FANCC variants, if any, is much smaller than for BRCA1, BRCA2 or PALB2 mutations. If this applies to all truncating variants in FANCC it would suggest there are differences between FA genes in their roles on breast cancer risk and demonstrates the merit of large consortia for clarifying risk associations of rare variants.
MYC paralogs are frequently activated in small cell lung cancer (SCLC) but represent poor drug targets. Thus, a detailed mapping of MYC-paralog-specific vulnerabilities may help to develop effective therapies for SCLC patients. Using a unique cellular CRISPR activation model, we uncover that, in contrast to MYCN and MYCL, MYC represses BCL2 transcription via interaction with MIZ1 and DNMT3a. The resulting lack of BCL2 expression promotes sensitivity to cell cycle control inhibition and dependency on MCL1. Furthermore, MYC activation leads to heightened apoptotic priming, intrinsic genotoxic stress and susceptibility to DNA damage checkpoint inhibitors. Finally, combined AURK and CHK1 inhibition substantially prolongs the survival of mice bearing MYC-driven SCLC beyond that of combination chemotherapy. These analyses uncover MYC-paralog-specific regulation of the apoptotic machinery with implications for genotype-based selection of targeted therapeutics in SCLC patients.
Animals must slow or halt locomotion to integrate sensory inputs or to change direction. In Caenorhabditis elegans, the GABAergic and peptidergic neuron RIS mediates developmentally timed quiescence. Here, we show RIS functions additionally as a locomotion stop neuron. RIS optogenetic stimulation caused acute and persistent inhibition of locomotion and pharyngeal pumping, phenotypes requiring FLP-11 neuropeptides and GABA. RIS photoactivation allows the animal to maintain its body posture by sustaining muscle tone, yet inactivating motor neuron oscillatory activity. During locomotion, RIS axonal Ca2+ signals revealed functional compartmentalization: Activity in the nerve ring process correlated with locomotion stop, while activity in a branch correlated with induced reversals. GABA was required to induce, and FLP-11 neuropeptides were required to sustain locomotion stop. RIS attenuates neuronal activity and inhibits movement, possibly enabling sensory integration and decision making, and exemplifies dual use of one cell across development in a compact nervous system.
Coordinated regulation of the lysosomal and autophagic systems ensures basal catabolism and normal cell physiology, and failure of either system causes disease. Here we describe an epigenetic rheostat orchestrated by c-MYC and histone deacetylases that inhibits lysosomal and autophagic biogenesis by concomitantly repressing the expression of the transcription factors MiT/TFE and FOXH1, and that of lysosomal and autophagy genes. Inhibition of histone deacetylases abates c-MYC binding to the promoters of lysosomal and autophagy genes, granting promoter occupancy to the MiT/TFE members, TFEB and TFE3, and/or the autophagy regulator FOXH1. In pluripotent stem cells and cancer, suppression of lysosomal and autophagic function is directly downstream of c-MYC overexpression and may represent a hallmark of malignant transformation. We propose that, by determining the fate of these catabolic systems, this hierarchical switch regulates the adaptive response of cells to pathological and physiological cues that could be exploited therapeutically.
Background
Shotgun metagenomes contain a sample of all the genomic material in an environment, allowing for the characterization of a microbial community. In order to understand these communities, bioinformatics methods are crucial. A common first step in processing metagenomes is to compute abundance estimates of different taxonomic or functional groups from the raw sequencing data.
Given the breadth of the field, computational solutions need to be flexible and extensible, enabling the combination of different tools into a larger pipeline.
Results
We present NGLess and NG-meta-profiler. NGLess is a domain specific language for describing next-generation sequence processing pipelines. It was developed with the goal of enabling user-friendly computational reproducibility. It provides built-in support for many common operations on sequencing data and is extensible with external tools with configuration files.
Using this framework, we developed NG-meta-profiler, a fast profiler for metagenomes which performs sequence preprocessing, mapping to bundled databases, filtering of the mapping results, and profiling (taxonomic and functional). It is significantly faster than either MOCAT2 or htseq-count and (as it builds on NGLess) its results are perfectly reproducible.
Conclusions
NG-meta-profiler is a high-performance solution for metagenomics processing built on NGLess. It can be used as-is to execute standard analyses or serve as the starting point for customization in a perfectly reproducible fashion.
NGLess and NG-meta-profiler are open source software (under the liberal MIT license) and can be downloaded from https://ngless.embl.de or installed through bioconda.
Zn\(^{2+}\) deficiency in the human population is frequent in underdeveloped countries. Worldwide, approximatively 2 billion people consume Zn\(^{2+}\)-deficient diets, accounting for 1–4% of deaths each year, mainly in infants with a compromised immune system. Depending on the severity of Zn\(^{2+}\) deficiency, clinical symptoms are associated with impaired wound healing, alopecia, diarrhea, poor growth, dysfunction of the immune and nervous system with congenital abnormalities and bleeding disorders. Poor nutritional Zn\(^{2+}\) status in patients with metastatic squamous cell carcinoma or with advanced non-Hodgkin lymphoma, was accompanied by cutaneous bleeding and platelet dysfunction. Forcing Zn\(^{2+}\) uptake in the gut using different nutritional supplementation of Zn\(^{2+}\) could ameliorate many of these pathological symptoms in humans. Feeding adult rodents with a low Zn\(^{2+}\) diet caused poor platelet aggregation and increased bleeding tendency, thereby attracting great scientific interest in investigating the role of Zn\(^{2+}\) in hemostasis. Storage protein metallothionein maintains or releases Zn\(^{2+}\) in the cytoplasm, and the dynamic change of this cytoplasmic Zn\(^{2+}\) pool is regulated by the redox status of the cell. An increase of labile Zn\(^{2+}\) pool can be toxic for the cells, and therefore cytoplasmic Zn\(^{2+}\) levels are tightly regulated by several Zn\(^{2+}\) transporters located on the cell surface and also on the intracellular membrane of Zn\(^{2+}\) storage organelles, such as secretory vesicles, endoplasmic reticulum or Golgi apparatus. Although Zn\(^{2+}\) is a critical cofactor for more than 2000 transcription factors and 300 enzymes, regulating cell differentiation, proliferation, and basic metabolic functions of the cells, the molecular mechanisms of Zn\(^{2+}\) transport and the physiological role of Zn\(^{2+}\) store in megakaryocyte and platelet function remain elusive. In this review, we summarize the contribution of extracellular or intracellular Zn\(^{2+}\) to megakaryocyte and platelet function and discuss the consequences of dysregulated Zn\(^{2+}\) homeostasis in platelet-related diseases by focusing on thrombosis, ischemic stroke and storage pool diseases.
Mathematical optimization framework allows the identification of certain nodes within a signaling network. In this work, we analyzed the complex extracellular-signal-regulated kinase 1 and 2 (ERK1/2) cascade in cardiomyocytes using the framework to find efficient adjustment screws for this cascade that is important for cardiomyocyte survival and maladaptive heart muscle growth. We modeled optimal pharmacological intervention points that are beneficial for the heart, but avoid the occurrence of a maladaptive ERK1/2 modification, the autophosphorylation of ERK at threonine 188 (ERK\(^{Thr188}\) phosphorylation), which causes cardiac hypertrophy. For this purpose, a network of a cardiomyocyte that was fitted to experimental data was equipped with external stimuli that model the pharmacological intervention points. Specifically, two situations were considered. In the first one, the cardiomyocyte was driven to a desired expression level with different treatment strategies. These strategies were quantified with respect to beneficial effects and maleficent side effects and then which one is the best treatment strategy was evaluated. In the second situation, it was shown how to model constitutively activated pathways and how to identify drug targets to obtain a desired activity level that is associated with a healthy state and in contrast to the maleficent expression pattern caused by the constitutively activated pathway. An implementation of the algorithms used for the calculations is also presented in this paper, which simplifies the application of the presented framework for drug targeting, optimal drug combinations and the systematic and automatic search for pharmacological intervention points. The codes were designed such that they can be combined with any mathematical model given by ordinary differential equations.