Refine
Has Fulltext
- yes (121)
Is part of the Bibliography
- yes (121)
Year of publication
Document Type
- Journal article (121) (remove)
Keywords
- metabolism (8)
- apoptosis (7)
- SARS-CoV-2 (5)
- cytokinins (5)
- COVID-19 (4)
- leukemic cells (4)
- regulation (4)
- transcriptome (4)
- Candida albicans (3)
- bioinformatics (3)
Institute
- Theodor-Boveri-Institut für Biowissenschaften (119)
- Klinik und Poliklinik für Anästhesiologie (ab 2004) (6)
- Lehrstuhl für Tissue Engineering und Regenerative Medizin (6)
- Institut für Molekulare Infektionsbiologie (5)
- Center for Computational and Theoretical Biology (4)
- Institut für Experimentelle Biomedizin (4)
- Institut für Pharmakologie und Toxikologie (4)
- Institut für Virologie und Immunbiologie (4)
- Kinderklinik und Poliklinik (3)
- Institut für Humangenetik (2)
Sonstige beteiligte Institutionen
EU-Project number / Contract (GA) number
- 031A408B (1)
- CoG 721016–HERPES (1)
- ESF-ZDEX 4.0 (1)
Fungal infections are a major global health burden where Candida albicans is among the most common fungal pathogen in humans and is a common cause of invasive candidiasis. Fungal phenotypes, such as those related to morphology, proliferation and virulence are mainly driven by gene expression, which is primarily regulated by kinase signaling cascades. Serine-arginine (SR) protein kinases are highly conserved among eukaryotes and are involved in major transcriptional processes in human and S. cerevisiae. Candida albicans harbors two SR protein kinases, while Sky2 is important for metabolic adaptation, Sky1 has similar functions as in S. cerevisiae. To investigate the role of these SR kinases for the regulation of transcriptional responses in C. albicans, we performed RNA sequencing of sky1Δ and sky2Δ and integrated a comprehensive phosphoproteome dataset of these mutants. Using a Systems Biology approach, we study transcriptional regulation in the context of kinase signaling networks. Transcriptomic enrichment analysis indicates that pathways involved in the regulation of gene expression are downregulated and mitochondrial processes are upregulated in sky1Δ. In sky2Δ, primarily metabolic processes are affected, especially for arginine, and we observed that arginine-induced hyphae formation is impaired in sky2Δ. In addition, our analysis identifies several transcription factors as potential drivers of the transcriptional response. Among these, a core set is shared between both kinase knockouts, but it appears to regulate different subsets of target genes. To elucidate these diverse regulatory patterns, we created network modules by integrating the data of site-specific protein phosphorylation and gene expression with kinase-substrate predictions and protein-protein interactions. These integrated signaling modules reveal shared parts but also highlight specific patterns characteristic for each kinase. Interestingly, the modules contain many proteins involved in fungal morphogenesis and stress response. Accordingly, experimental phenotyping shows a higher resistance to Hygromycin B for sky1Δ. Thus, our study demonstrates that a combination of computational approaches with integration of experimental data can offer a new systems biological perspective on the complex network of signaling and transcription. With that, the investigation of the interface between signaling and transcriptional regulation in C. albicans provides a deeper insight into how cellular mechanisms can shape the phenotype.
The ITS2 Database
(2012)
The internal transcribed spacer 2 (ITS2) has been used as a phylogenetic marker for more than two decades. As ITS2 research mainly focused on the very variable ITS2 sequence, it confined this marker to low-level phylogenetics only. However, the combination of the ITS2 sequence and its highly conserved secondary structure improves the phylogenetic resolution1 and allows phylogenetic inference at multiple taxonomic ranks, including species delimitation.
The ITS2 Database presents an exhaustive dataset of internal transcribed spacer 2 sequences from NCBI GenBank accurately reannotated. Following an annotation by profile Hidden Markov Models (HMMs), the secondary structure of each sequence is predicted. First, it is tested whether a minimum energy based fold (direct fold) results in a correct, four helix conformation. If this is not the case, the structure is predicted by homology modeling. In homology modeling, an already known secondary structure is transferred to another ITS2 sequence, whose secondary structure was not able to fold correctly in a direct fold.
The ITS2 Database is not only a database for storage and retrieval of ITS2 sequence-structures. It also provides several tools to process your own ITS2 sequences, including annotation, structural prediction, motif detection and BLAST search on the combined sequence-structure information. Moreover, it integrates trimmed versions of 4SALE and ProfDistS for multiple sequence-structure alignment calculation and Neighbor Joining tree reconstruction. Together they form a coherent analysis pipeline from an initial set of sequences to a phylogeny based on sequence and secondary structure.
In a nutshell, this workbench simplifies first phylogenetic analyses to only a few mouse-clicks, while additionally providing tools and data for comprehensive large-scale analyses.
Plant hormones involving salicylic acid (SA), jasmonic acid (JA), ethylene (Et), and auxin, gibberellins, and abscisic acid (ABA) are known to regulate host immune responses. However, plant hormone cytokinin has the potential to modulate defense signaling including SA and JA. It promotes plant pathogen and herbivore resistance; underlying mechanisms are still unknown. Using systems biology approaches, we unravel hub points of immune interaction mediated by cytokinin signaling in Arabidopsis. High-confidence Arabidopsis protein-protein interactions (PPI) are coupled to changes in cytokinin-mediated gene expression. Nodes of the cellular interactome that are enriched in immune functions also reconstitute sub-networks. Topological analyses and their specific immunological relevance lead to the identification of functional hubs in cellular interactome. We discuss our identified immune hubs in light of an emerging model of cytokinin-mediated immune defense against pathogen infection in plants.
The growing tips of plants grow sterile; therefore, disease-free plants can be generated from them. How plants safeguard growing apices from pathogen infection is still a mystery. The shoot apical meristem (SAM) is one of the three stem cells niches that give rise to the above ground plant organs. This is very well explored; however, how signaling networks orchestrate immune responses against pathogen infections in the SAM remains unclear. To reconstruct a transcriptional framework of the differentially expressed genes (DEGs) pertaining to various SAM cellular populations, we acquired large-scale transcriptome datasets from the public repository Gene Expression Omnibus (GEO). We identify here distinct sets of genes for various SAM cellular populations that are enriched in immune functions, such as immune defense, pathogen infection, biotic stress, and response to salicylic acid and jasmonic acid and their biosynthetic pathways in the SAM. We further linked those immune genes to their respective proteins and identify interactions among them by mapping a transcriptome-guided SAM-interactome. Furthermore, we compared stem-cells regulated transcriptome with innate immune responses in plants showing transcriptional separation among their DEGs in Arabidopsis. Besides unleashing a repertoire of immune-related genes in the SAM, our analysis provides a SAM-interactome that will help the community in designing functional experiments to study the specific defense dynamics of the SAM-cellular populations. Moreover, our study promotes the essence of large-scale omics data re-analysis, allowing a fresh look at the SAM-cellular transcriptome repurposing data-sets for new questions.
Metabolism and signaling of cytokinins was first established in plants, followed by cytokinin discoveries in all kingdoms of life. However, understanding of their role in mammalian cells is still scarce. Kinetin is a cytokinin that mitigates the effects of oxidative stress in mammalian cells. The effective concentrations of exogenously applied kinetin in invoking various cellular responses are not well standardized. Likewise, the metabolism of kinetin and its cellular targets within the mammalian cells are still not well studied. Applying vitality tests as well as comet assays under normal and hyper-oxidative states, our analysis suggests that kinetin concentrations of 500 nM and above cause cytotoxicity as well as genotoxicity in various cell types. However, concentrations below 100 nM do not cause any toxicity, rather in this range kinetin counteracts oxidative burst and cytotoxicity. We focus here on these effects. To get insights into the cellular targets of kinetin mediating these pro-survival functions and protective effects we applied structural and computational approaches on two previously testified targets for these effects. Our analysis deciphers vital residues in adenine phosphoribosyltransferase (APRT) and adenosine receptor (A2A-R) that facilitate the binding of kinetin to these two important human cellular proteins. We finally discuss how the therapeutic potential of kinetin against oxidative stress helps in various pathophysiological conditions.
Intricate mechanisms discriminate between friends and foes in plants. Plant organs deploy overlapping and distinct protection strategies. Despite vulnerability to a plethora of pathogens, the growing tips of plants grow bacteria free. The shoot apical meristem (SAM) is among three stem cells niches, a self-renewable reservoir for the future organogenesis of leaf, stem, and flowers. How plants safeguard this high value growth target from infections was not known until now. Recent reports find the stem cell secreted 12-amino acid peptide CLV3p (CLAVATA3 peptide) is perceived by FLS2 (FLAGELLIN SENSING 2) receptor and activates the transcription of immunity and defense marker genes. No infection in the SAM of wild type plants and bacterial infection in clv3 and fls2 mutants illustrate this natural protection against infections. Cytokinins (CKs) are enriched in the SAM and regulate meristem activities by their involvement in stem cell signaling networks. Auxin mediates plant susceptibility to pathogen infections while CKs boost plant immunity. Here, in addition to the stem-cell-triggered immunity we also highlight a potential link between CK signaling and CLV3p mediated immune response in the SAM.
Introduction
Pro-thrombotic events are one of the prevalent causes of intensive care unit (ICU) admissions among COVID-19 patients, although the signaling events in the stimulated platelets are still unclear.
Methods
We conducted a comparative analysis of platelet transcriptome data from healthy donors, ICU, and non-ICU COVID-19 patients to elucidate these mechanisms. To surpass previous analyses, we constructed models of involved networks and control cascades by integrating a global human signaling network with transcriptome data. We investigated the control of platelet hyperactivation and the specific proteins involved.
Results
Our study revealed that control of the platelet network in ICU patients is significantly higher than in non-ICU patients. Non-ICU patients require control over fewer proteins for managing platelet hyperactivity compared to ICU patients. Identification of indispensable proteins highlighted key subnetworks, that are targetable for system control in COVID-19-related platelet hyperactivity. We scrutinized FDA-approved drugs targeting indispensable proteins and identified fostamatinib as a potent candidate for preventing thrombosis in COVID-19 patients.
Discussion
Our findings shed light on how SARS-CoV-2 efficiently affects host platelets by targeting indispensable and critical proteins involved in the control of platelet activity. We evaluated several drugs for specific control of platelet hyperactivity in ICU patients suffering from platelet hyperactivation. The focus of our approach is repurposing existing drugs for optimal control over the signaling network responsible for platelet hyperactivity in COVID-19 patients. Our study offers specific pharmacological recommendations, with drug prioritization tailored to the distinct network states observed in each patient condition. Interactive networks and detailed results can be accessed at https://fostamatinib.bioinfo-wuerz.eu/.
Synthetically designed alternative photorespiratory pathways increase the biomass of tobacco and rice plants. Likewise, some in planta–tested synthetic carbon-concentrating cycles (CCCs) hold promise to increase plant biomass while diminishing atmospheric carbon dioxide burden. Taking these individual contributions into account, we hypothesize that the integration of bypasses and CCCs will further increase plant productivity. To test this in silico, we reconstructed a metabolic model by integrating photorespiration and photosynthesis with the synthetically designed alternative pathway 3 (AP3) enzymes and transporters. We calculated fluxes of the native plant system and those of AP3 combined with the inhibition of the glycolate/glycerate transporter by using the YANAsquare package. The activity values corresponding to each enzyme in photosynthesis, photorespiration, and for synthetically designed alternative pathways were estimated. Next, we modeled the effect of the crotonyl-CoA/ethylmalonyl-CoA/hydroxybutyryl-CoA cycle (CETCH), which is a set of natural and synthetically designed enzymes that fix CO₂ manifold more than the native Calvin–Benson–Bassham (CBB) cycle. We compared estimated fluxes across various pathways in the native model and under an introduced CETCH cycle. Moreover, we combined CETCH and AP3-w/plgg1RNAi, and calculated the fluxes. We anticipate higher carbon dioxide–harvesting potential in plants with an AP3 bypass and CETCH–AP3 combination. We discuss the in vivo implementation of these strategies for the improvement of C3 plants and in natural high carbon harvesters.
The present study reports the synthesis of new purine bioisosteres comprising a pyrazolo[3,4-d]pyrimidine scaffold linked to mono-, di-, and trimethoxy benzylidene moieties through hydrazine linkages. First, in silico docking experiments of the synthesized compounds against Bax, Bcl-2, Caspase-3, Ki67, p21, and p53 were performed in a trial to rationalize the observed cytotoxic activity for the tested compounds. The anticancer activity of these compounds was evaluated in vitro against Caco-2, A549, HT1080, and Hela cell lines. Results revealed that two (5 and 7) of the three synthesized compounds (5, 6, and 7) showed high cytotoxic activity against all tested cell lines with IC50 values in the micro molar concentration. Our in vitro results show that there is no significant apoptotic effect for the treatment with the experimental compounds on the viability of cells against A549 cells. Ki67 expression was found to decrease significantly following the treatment of cells with the most promising candidate: drug 7. The overall results indicate that these pyrazolopyrimidine derivatives possess anticancer activity at varying doses. The suggested mechanism of action involves the inhibition of the proliferation of cancer cells.
Plant hormones are small regulatory molecules that exert pharmacological actions in mammalian cells such as anti-oxidative and pro-metabolic effects. Kinetin belongs to the group of plant hormones cytokinin and has been associated with modulatory functions in mammalian cells. The mammalian adenosine receptor (A2a-R) is known to modulate multiple physiological responses in animal cells. Here, we describe that kinetin binds to the adenosine receptor (A2a-R) through the Asn253 residue in an adenosine dependent manner. To harness the beneficial effects of kinetin for future human use, we assess its acute toxicity by analyzing different biochemical and histological markers in rats. Kinetin at a dose below 1 mg/kg had no adverse effects on the serum level of glucose or on the activity of serum alanine transaminase (ALT) or aspartate aminotransferase (AST) enzymes in the kinetin treated rats. Whereas, creatinine levels increased after a kinetin treatment at a dose of 0.5 mg/kg. Furthermore, 5 mg/kg treated kinetin rats showed normal renal corpuscles, but a mild degeneration was observed in the renal glomeruli and renal tubules, as well as few degenerated hepatocytes were also observed in the liver. Kinetin doses below 5 mg/kg did not show any localized toxicity in the liver and kidney tissues. In addition to unraveling the binding interaction between kinetin and A2a-R, our findings suggest safe dose limits for the future use of kinetin as a therapeutic and modulatory agent against various pathophysiological conditions.
Modulating key dynamics of plant growth and development, the effects of the plant hormone cytokinin on animal cells gained much attention recently. Most previous studies on cytokinin effects on mammalian cells have been conducted with elevated cytokinin concentration (in the μM range). However, to examine physiologically relevant dose effects of cytokinins on animal cells, we systematically analyzed the impact of kinetin in cultured cells at low and high concentrations (1nM-10μM) and examined cytotoxic and genotoxic conditions. We furthermore measured the intrinsic antioxidant activity of kinetin in a cell-free system using the Ferric Reducing Antioxidant Power assay and in cells using the dihydroethidium staining method. Monitoring viability, we looked at kinetin effects in mammalian cells such as HL60 cells, HaCaT human keratinocyte cells, NRK rat epithelial kidney cells and human peripheral lymphocytes. Kinetin manifests no antioxidant activity in the cell free system and high doses of kinetin (500 nM and higher) reduce cell viability and mediate DNA damage in vitro. In contrast, low doses (concentrations up to 100 nM) of kinetin confer protection in cells against oxidative stress. Moreover, our results show that pretreatment of the cells with kinetin significantly reduces 4-nitroquinoline 1-oxide mediated reactive oxygen species production. Also, pretreatment with kinetin retains cellular GSH levels when they are also treated with the GSH-depleting agent patulin. Our results explicitly show that low kinetin doses reduce apoptosis and protect cells from oxidative stress mediated cell death. Future studies on the interaction between cytokinins and human cellular pathway targets will be intriguing.
Background: Tumor necrosis factor-like weak inducer of apoptosis (TWEAK) and its receptor fibroblast growth factorinducible 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.
Epithelial-to-mesenchymal transition (EMT) is discussed to be centrally involved in invasion, stemness, and drug resistance. Experimental models to evaluate this process in its biological complexity are limited. To shed light on EMT impact and test drug response more reliably, we use a lung tumor test system based on a decellularized intestinal matrix showing more in vivo-like proliferation levels and enhanced expression of clinical markers and carcinogenesis-related genes. In our models, we found evidence for a correlation of EMT with drug resistance in primary and secondary resistant cells harboring KRAS\(^{G12C}\) or EGFR mutations, which was simulated in silico based on an optimized signaling network topology. Notably, drug resistance did not correlate with EMT status in KRAS-mutated patient-derived xenograft (PDX) cell lines, and drug efficacy was not affected by EMT induction via TGF-β. To investigate further determinants of drug response, we tested several drugs in combination with a KRAS\(^{G12C}\) inhibitor in KRAS\(^{G12C}\) mutant HCC44 models, which, besides EMT, display mutations in P53, LKB1, KEAP1, and high c-MYC expression. We identified an aurora-kinase A (AURKA) inhibitor as the most promising candidate. In our network, AURKA is a centrally linked hub to EMT, proliferation, apoptosis, LKB1, and c-MYC. This exemplifies our systemic analysis approach for clinical translation of biomarker signatures.
For SARS-CoV-2, R0 calculations in the range of 2–3 dominate the literature, but much higher estimates have also been published. Because capacity for RT-PCR testing increased greatly in the early phase of the Covid-19 pandemic, R0 determinations based on these incidence values are subject to strong bias. We propose to use Covid-19-induced excess mortality to determine R0 regardless of RT-PCR testing capacity. We used data from the Robert Koch Institute (RKI) on the incidence of Covid cases, Covid-related deaths, number of RT-PCR tests performed, and excess mortality calculated from data from the Federal Statistical Office in Germany. We determined R0 using exponential growth estimates with a serial interval of 4.7 days. We used only datasets that were not yet under the influence of policy measures (e.g., lockdowns or school closures). The uncorrected R0 value for the spread of SARS-CoV-2 based on RT-PCR incidence data was 2.56 (95% CI 2.52–2.60) for Covid-19 cases and 2.03 (95% CI 1.96–2.10) for Covid-19-related deaths. However, because the number of RT-PCR tests increased by a growth factor of 1.381 during the same period, these R0 values must be corrected accordingly (R0corrected = R0uncorrected/1.381), yielding 1.86 for Covid-19 cases and 1.47 for Covid-19 deaths. The R0 value based on excess deaths was calculated to be 1.34 (95% CI 1.32–1.37). A sine-function-based adjustment for seasonal effects of 40% corresponds to a maximum value of R0January = 1.68 and a minimum value of R0July = 1.01. Our calculations show an R0 that is much lower than previously thought. This relatively low range of R0 fits very well with the observed seasonal pattern of infection across Europe in 2020 and 2021, including the emergence of more contagious escape variants such as delta or omicron. In general, our study shows that excess mortality can be used as a reliable surrogate to determine the R0 in pandemic situations.
The variable regions (V1–V9) of the 18S rDNA are routinely used in barcoding and phylogenetics. In handling these data for trypanosomes, we have noticed a misunderstanding that has apparently taken a life of its own in the literature over the years. In particular, in recent years, when studying the phylogenetic relationship of trypanosomes, the use of V7/V8 was systematically established. However, considering the current numbering system for all other organisms (including other Euglenozoa), V7/V8 was never used. In Maia da Silva et al. [Parasitology 2004, 129, 549–561], V7/V8 was promoted for the first time for trypanosome phylogenetics, and since then, more than 70 publications have replicated this nomenclature and even discussed the benefits of the use of this region in comparison to V4. However, the primers used to amplify the variable region of trypanosomes have actually amplified V4 (concerning the current 18S rDNA numbering system).
Molecular characterization of antimicrobial peptide genes of the carpenter ant Camponotus floridanus
(2012)
The production of antimicrobial peptides (AMPs) is a major defense mechanism against pathogen infestation and of particular importance for insects relying exclusively on an innate immune system. Here, we report on the characterization of three AMPs from the carpenter ant Camponotus floridanus. Due to sequence similarities and amino acid composition these peptides can be classified into the cysteine-rich (e.g. defensin) and glycine-rich (e.g. hymenoptaecin) AMP groups, respectively. The gene and cDNA sequences of these AMPs were established and their expression was shown to be induced by microbial challenge. We characterized two different defensin genes. The defensin-2 gene has a single intron, whereas the defensin-1 gene has two introns. The deduced amino acid sequence of the C. floridanus defensins is very similar to other known ant defensins with the exception of a short C-terminal extension of defensin-1. The hymenoptaecin gene has a single intron and a very peculiar domain structure. The corresponding precursor protein consists of a signal- and a pro-sequence followed by a hymenoptaecin-like domain and six directly repeated hymenoptaecin domains. Each of the hymenoptaecin domains is flanked by an EAEP-spacer sequence and a RR-site known to be a proteolytic processing site. Thus, proteolytic processing of the multipeptide precursor may generate several mature AMPs leading to an amplification of the immune response. Bioinformatical analyses revealed the presence of hymenoptaecin genes with similar multipeptide precursor structure in genomes of other ant species suggesting an evolutionary conserved important role of this gene in ant immunity.
Fungal microorganisms frequently lead to life-threatening infections. Within this group of pathogens, the commensal Candida albicans and the filamentous fungus Aspergillus fumigatus are by far the most important causes of invasive mycoses in Europe. A key capability for host invasion and immune response evasion are specific molecular interactions between the fungal pathogen and its human host. Experimentally validated knowledge about these crucial interactions is rare in literature and even specialized host pathogen databases mainly focus on bacterial and viral interactions whereas information on fungi is still sparse. To establish large-scale host fungi interaction networks on a systems biology scale, we develop an extended inference approach based on protein orthology and data on gene functions. Using human and yeast intraspecies networks as template, we derive a large network of pathogen host interactions (PHI). Rigorous filtering and refinement steps based on cellular localization and pathogenicity information of predicted interactors yield a primary scaffold of fungi human and fungi mouse interaction networks. Specific enrichment of known pathogenicity-relevant genes indicates the biological relevance of the predicted PHI. A detailed inspection of functionally relevant subnetworks reveals novel host fungal interaction candidates such as the Candida virulence factor PLB1 and the anti-fungal host protein APP. Our results demonstrate the applicability of interolog-based prediction methods for host fungi interactions and underline the importance of filtering and refinement steps to attain biologically more relevant interactions. This integrated network framework can serve as a basis for future analyses of high-throughput host fungi transcriptome and proteome data.
PRO-Simat is a simulation tool for analysing protein interaction networks, their dynamic change and pathway engineering. It provides GO enrichment, KEGG pathway analyses, and network visualisation from an integrated database of more than 8 million protein-protein interactions across 32 model organisms and the human proteome. We integrated dynamical network simulation using the Jimena framework, which quickly and efficiently simulates Boolean genetic regulatory networks. It enables simulation outputs with in-depth analysis of the type, strength, duration and pathway of the protein interactions on the website. Furthermore, the user can efficiently edit and analyse the effect of network modifications and engineering experiments. In case studies, applications of PRO-Simat are demonstrated: (i) understanding mutually exclusive differentiation pathways in Bacillus subtilis, (ii) making Vaccinia virus oncolytic by switching on its viral replication mainly in cancer cells and triggering cancer cell apoptosis and (iii) optogenetic control of nucleotide processing protein networks to operate DNA storage. Multilevel communication between components is critical for efficient network switching, as demonstrated by a general census on prokaryotic and eukaryotic networks and comparing design with synthetic networks using PRO-Simat. The tool is available at https://prosimat.heinzelab.de/ as a web-based query server.
An epidemic of avian type H7N9 influenza virus, which took place in China in 2013, was enhanced by a naturally occurring R294K mutation resistant against Oseltamivir at the catalytic site of the neuraminidase. To cope with such drug-resistant neuraminidase mutations, we applied the molecular docking technique to evaluate the fitness of the available drugs such as Oseltamivir, Zanamivir, Peramivir, Laninamivir, L-Arginine and Benserazide hydrochloride concerning the N9 enzyme with single (R294K, R119K, R372K), double (R119_294K, R119_372K, R294_372K) and triple (R119_294_372K) mutations in the pocket. We found that the drugs Peramivir and Zanamivir score best amongst the studied compounds, demonstrating their high binding potential towards the pockets with the considered mutations. Despite the fact that mutations changed the shape of the pocket and reduced the binding strength for all drugs, Peramivir was the only drug that formed interactions with the key residues at positions 119, 294 and 372 in the pocket of the triple N9 mutant, while Zanamivir demonstrated the lowest RMSD value (0.7 Å) with respect to the reference structure.