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)
Aspergillus is an important fungal genus containing economically important species, as well as pathogenic species of animals and plants. Using eighteen fungal species of the genus Aspergillus, we conducted a comprehensive investigation of conserved genes and their evolution. This also allows us to investigate the selection pressure driving the adaptive evolution in the pathogenic species A. fumigatus. Among single-copy orthologs (SCOs) for A. fumigatus and the closely related species A. fischeri, we identified 122 versus 50 positively selected genes (PSGs), respectively. Moreover, twenty conserved genes of unknown function were established to be positively selected and thus important for adaption. A. fumigatus PSGs interacting with human host proteins show over-representation of adaptive, symbiosis-related, immunomodulatory and virulence-related pathways, such as the TGF-β pathway, insulin receptor signaling, IL1 pathway and interfering with phagosomal GTPase signaling. Additionally, among the virulence factor coding genes, secretory and membrane protein-coding genes in multi-copy gene families, 212 genes underwent positive selection and also suggest increased adaptation, such as fungal immune evasion mechanisms (aspf2), siderophore biosynthesis (sidD), fumarylalanine production (sidE), stress tolerance (atfA) and thermotolerance (sodA). These genes presumably contribute to host adaptation strategies. Genes for the biosynthesis of gliotoxin are shared among all the close relatives of A. fumigatus as an ancient defense mechanism. Positive selection plays a crucial role in the adaptive evolution of A. fumigatus. The genome-wide profile of PSGs provides valuable targets for further research on the mechanisms of immune evasion, antimycotic targeting and understanding fundamental virulence processes.
Patient-tailored therapy based on tumor drivers is promising for lung cancer treatment. For this, we combined in vitro tissue models with in silico analyses. Using individual cell lines with specific mutations, we demonstrate a generic and rapid stratification pipeline for targeted tumor therapy. We improve in vitro models of tissue conditions by a biological matrix-based three-dimensional (3D) tissue culture that allows in vitro drug testing: It correctly shows a strong drug response upon gefitinib (Gef) treatment in a cell line harboring an EGFR-activating mutation (HCC827), but no clear drug response upon treatment with the HSP90 inhibitor 17AAG in two cell lines with KRAS mutations (H441, A549). In contrast, 2D testing implies wrongly KRAS as a biomarker for HSP90 inhibitor treatment, although this fails in clinical studies. Signaling analysis by phospho-arrays showed similar effects of EGFR inhibition by Gef in HCC827 cells, under both 2D and 3D conditions. Western blot analysis confirmed that for 3D conditions, HSP90 inhibitor treatment implies different p53 regulation and decreased MET inhibition in HCC827 and H441 cells. Using in vitro data (western, phospho-kinase array, proliferation, and apoptosis), we generated cell line-specific in silico topologies and condition-specific (2D, 3D) simulations of signaling correctly mirroring in vitro treatment responses. Networks predict drug targets considering key interactions and individual cell line mutations using the Human Protein Reference Database and the COSMIC database. A signature of potential biomarkers and matching drugs improve stratification and treatment in KRAS-mutated tumors. In silico screening and dynamic simulation of drug actions resulted in individual therapeutic suggestions, that is, targeting HIF1A in H441 and LKB1 in A549 cells. In conclusion, our in vitro tumor tissue model combined with an in silico tool improves drug effect prediction and patient stratification. Our tool is used in our comprehensive cancer center and is made now publicly available for targeted therapy decisions.
Glycoprotein VI (GPVI) is a platelet-specific receptor for collagen and fibrin, regulating important platelet functions such as platelet adhesion and thrombus growth. Although the blockade of GPVI function is widely recognized as a potent anti-thrombotic approach, there are limited studies focused on site-specific targeting of GPVI. Using computational modeling and bioinformatics, we analyzed collagen- and CRP-binding surfaces of GPVI monomers and dimers, and compared the interacting surfaces with other mammalian GPVI isoforms. We could predict a minimal collagen-binding epitope of GPVI dimer and designed an EA-20 antibody that recognizes a linear epitope of this surface. Using platelets and whole blood samples donated from wild-type and humanized GPVI transgenic mice and also humans, our experimental results show that the EA-20 antibody inhibits platelet adhesion and aggregation in response to collagen and CRP, but not to fibrin. The EA-20 antibody also prevents thrombus formation in whole blood, on the collagen-coated surface, in arterial flow conditions. We also show that EA-20 does not influence GPVI clustering or receptor shedding. Therefore, we propose that blockade of this minimal collagen-binding epitope of GPVI with the EA-20 antibody could represent a new anti-thrombotic approach by inhibiting specific interactions between GPVI and the collagen matrix.
Automatic image reconstruction is critical to cope with steadily increasing data from advanced microscopy. We describe here the Fiji macro 3D ART VeSElecT which we developed to study synaptic vesicles in electron tomograms. We apply this tool to quantify vesicle properties (i) in embryonic Danio rerio 4 and 8 days past fertilization (dpf) and (ii) to compare Caenorhabditis elegans N2 neuromuscular junctions (NMJ) wild-type and its septin mutant (unc-59(e261)). We demonstrate development-specific and mutant-specific changes in synaptic vesicle pools in both models. We confirm the functionality of our macro by applying our 3D ART VeSElecT on zebrafish NMJ showing smaller vesicles in 8 dpf embryos then 4 dpf, which was validated by manual reconstruction of the vesicle pool. Furthermore, we analyze the impact of C. elegans septin mutant unc-59(e261) on vesicle pool formation and vesicle size. Automated vesicle registration and characterization was implemented in Fiji as two macros (registration and measurement). This flexible arrangement allows in particular reducing false positives by an optional manual revision step. Preprocessing and contrast enhancement work on image-stacks of 1nm/pixel in x and y direction. Semi-automated cell selection was integrated. 3D ART VeSElecT removes interfering components, detects vesicles by 3D segmentation and calculates vesicle volume and diameter (spherical approximation, inner/outer diameter). Results are collected in color using the RoiManager plugin including the possibility of manual removal of non-matching confounder vesicles. Detailed evaluation considered performance (detected vesicles) and specificity (true vesicles) as well as precision and recall. We furthermore show gain in segmentation and morphological filtering compared to learning based methods and a large time gain compared to manual segmentation. 3D ART VeSElecT shows small error rates and its speed gain can be up to 68 times faster in comparison to manual annotation. Both automatic and semi-automatic modes are explained including a tutorial.
Synaptic vesicles (SVs) are a key component of neuronal signaling and fulfil different roles depending on their composition. In electron micrograms of neurites, two types of vesicles can be distinguished by morphological criteria, the classical “clear core” vesicles (CCV) and the typically larger “dense core” vesicles (DCV), with differences in electron density due to their diverse cargos. Compared to CCVs, the precise function of DCVs is less defined. DCVs are known to store neuropeptides, which function as neuronal messengers and modulators [1]. In C. elegans, they play a role in locomotion, dauer formation, egg-laying, and mechano- and chemosensation [2]. Another type of DCVs, also referred to as granulated vesicles, are known to transport Bassoon, Piccolo and further constituents of the presynaptic density in the center of the active zone (AZ), and therefore are important for synaptogenesis [3].
To better understand the role of different types of SVs, we present here a new automated approach to classify vesicles. We combine machine learning with an extension of our previously developed vesicle segmentation workflow, the ImageJ macro 3D ART VeSElecT. With that we reliably distinguish CCVs and DCVs in electron tomograms of C. elegans NMJs using image-based features. Analysis of the underlying ground truth data shows an increased fraction of DCVs as well as a higher mean distance between DCVs and AZs in dauer larvae compared to young adult hermaphrodites. Our machine learning based tools are adaptable and can be applied to study properties of different synaptic vesicle pools in electron tomograms of diverse model organisms.
The signal modelling framework JimenaE simulates dynamically Boolean networks. In contrast to SQUAD, there is systematic and not just heuristic calculation of all system states. These specific features are not present in CellNetAnalyzer and BoolNet. JimenaE is an expert extension of Jimena, with new optimized code, network conversion into different formats, rapid convergence both for system state calculation as well as for all three network centralities. It allows higher accuracy in determining network states and allows to dissect networks and identification of network control type and amount for each protein with high accuracy. Biological examples demonstrate this: (i) High plasticity of mesenchymal stromal cells for differentiation into chondrocytes, osteoblasts and adipocytes and differentiation-specific network control focusses on wnt-, TGF-beta and PPAR-gamma signaling. JimenaE allows to study individual proteins, removal or adding interactions (or autocrine loops) and accurately quantifies effects as well as number of system states. (ii) Dynamical modelling of cell–cell interactions of plant Arapidopsis thaliana against Pseudomonas syringae DC3000: We analyze for the first time the pathogen perspective and its interaction with the host. We next provide a detailed analysis on how plant hormonal regulation stimulates specific proteins and who and which protein has which type and amount of network control including a detailed heatmap of the A.thaliana response distinguishing between two states of the immune response. (iii) In an immune response network of dendritic cells confronted with Aspergillus fumigatus, JimenaE calculates now accurately the specific values for centralities and protein-specific network control including chemokine and pattern recognition receptors.
New antimycotic drugs are challenging to find, as potential target proteins may have close human orthologs. We here focus on identifying metabolic targets that are critical for fungal growth and have minimal similarity to targets among human proteins. We compare and combine here: (I) direct metabolic network modeling using elementary mode analysis and flux estimates approximations using expression data, (II) targeting metabolic genes by transcriptome analysis of condition-specific highly expressed enzymes, and (III) analysis of enzyme structure, enzyme interconnectedness (“hubs”), and identification of pathogen-specific enzymes using orthology relations. We have identified 64 targets including metabolic enzymes involved in vitamin synthesis, lipid, and amino acid biosynthesis including 18 targets validated from the literature, two validated and five currently examined in own genetic experiments, and 38 further promising novel target proteins which are non-orthologous to human proteins, involved in metabolism and are highly ranked drug targets from these pipelines.
Background: Chagas disease (CD) is a major burden in Latin America, expanding also to non-endemic countries. A gold standard to detect the CD causing pathogen Trypanosoma cruzi is currently not available. Existing real time polymerase chain reactions (RT-PCRs) lack sensitivity and/or specificity. We present a new, highly specific RT-PCR for the diagnosis and monitoring of CD. Material and Methods: We analyzed 352 serum samples from Indigenous people living in high endemic CD areas of Colombia using three leading RT-PCRs (k-DNA-, TCZ-, 18S rRNA-PCR), the newly developed one (NDO-PCR), a Rapid Test/enzyme-linked immuno sorbent assay (ELISA), and immunofluorescence. Eighty-seven PCR-products were verified by sequence analysis after plasmid vector preparation. Results: The NDO-PCR showed the highest sensitivity (92.3%), specificity (100%), and accuracy (94.3%) for T. cruzi detection in the 87 sequenced samples. Sensitivities and specificities of the kDNA-PCR were 89.2%/22.7%, 20.5%/100% for TCZ-PCR, and 1.5%/100% for the 18S rRNA-PCR. The kDNA-PCR revealed a 77.3% false positive rate, mostly due to cross-reactions with T. rangeli (NDO-PCR 0%). TCZ- and 18S rRNA-PCR showed a false negative rate of 79.5% and 98.5% (NDO-PCR 7.7%), respectively. Conclusions: The NDO-PCR demonstrated the highest specificity, sensitivity, and accuracy compared to leading PCRs. Together with serologic tests, it can be considered as a reliable tool for CD detection and can improve CD management significantly.
Background: Boolean networks capture switching behavior of many naturally occurring regulatory networks. For semi-quantitative modeling, interpolation between ON and OFF states is necessary. The high degree polynomial interpolation of Boolean genetic regulatory networks (GRNs) in cellular processes such as apoptosis or proliferation allows for the modeling of a wider range of node interactions than continuous activator-inhibitor models, but suffers from scaling problems for networks which contain nodes with more than ~10 inputs. Many GRNs from literature or new gene expression experiments exceed those limitations and a new approach was developed.
Results: (i) As a part of our new GRN simulation framework Jimena we introduce and setup Boolean-tree-based data structures; (ii) corresponding algorithms greatly expedite the calculation of the polynomial interpolation in almost all cases, thereby expanding the range of networks which can be simulated by this model in reasonable time. (iii) Stable states for discrete models are efficiently counted and identified using binary decision diagrams. As application example, we show how system states can now be sampled efficiently in small up to large scale hormone disease networks (Arabidopsis thaliana development and immunity, pathogen Pseudomonas syringae and modulation by cytokinins and plant hormones).
Conclusions: Jimena simulates currently available GRNs about 10-100 times faster than the previous implementation of the polynomial interpolation model and even greater gains are achieved for large scale-free networks. This speed-up also facilitates a much more thorough sampling of continuous state spaces which may lead to the identification of new stable states. Mutants of large networks can be constructed and analyzed very quickly enabling new insights into network robustness and behavior.
Control of genetic regulatory networks is challenging to define and quantify. Previous control centrality metrics, which aim to capture the ability of individual nodes to control the system, have been found to suffer from plausibility and applicability problems. Here we present a new approach to control centrality based on network convergence behaviour, implemented as an extension of our genetic regulatory network simulation framework Jimena (http://stefan-karl.de/jimena). We distinguish three types of network control, and show how these mathematical concepts correspond to experimentally verified node functions and signalling pathways in immunity and cell differentiation: Total control centrality quantifies the impact of node mutations and identifies potential pharmacological targets such as genes involved in oncogenesis (e.g. zinc finger protein GLI2 or bone morphogenetic proteins in chondrocytes). Dynamic control centrality describes relaying functions as observed in signalling cascades (e.g. src kinase or Jak/Stat pathways). Value control centrality measures the direct influence of the value of the node on the network (e.g. Indian hedgehog as an essential regulator of proliferation in chondrocytes). Surveying random scale-free networks and biological networks, we find that control of the network resides in few high degree driver nodes and networks can be controlled best if they are sparsely connected.
Background: In several studies, secondary structures of ribosomal genes have been used to improve the quality of phylogenetic reconstructions. An extensive evaluation of the benefits of secondary structure, however, is lacking. Results: This is the first study to counter this deficiency. We inspected the accuracy and robustness of phylogenetics with individual secondary structures by simulation experiments for artificial tree topologies with up to 18 taxa and for divergency levels in the range of typical phylogenetic studies. We chose the internal transcribed spacer 2 of the ribosomal cistron as an exemplary marker region. Simulation integrated the coevolution process of sequences with secondary structures. Additionally, the phylogenetic power of marker size duplication was investigated and compared with sequence and sequence-structure reconstruction methods. The results clearly show that accuracy and robustness of Neighbor Joining trees are largely improved by structural information in contrast to sequence only data, whereas a doubled marker size only accounts for robustness. Conclusions: Individual secondary structures of ribosomal RNA sequences provide a valuable gain of information content that is useful for phylogenetics. Thus, the usage of ITS2 sequence together with secondary structure for taxonomic inferences is recommended. Other reconstruction methods as maximum likelihood, bayesian inference or maximum parsimony may equally profit from secondary structure inclusion. Reviewers: This article was reviewed by Shamil Sunyaev, Andrea Tanzer (nominated by Frank Eisenhaber) and Eugene V. Koonin. Open peer review: Reviewed by Shamil Sunyaev, Andrea Tanzer (nominated by Frank Eisenhaber) and Eugene V. Koonin. For the full reviews, please go to the Reviewers’ comments section.
Cyclin-dependent kinase-like kinases (CLKs) are dual specificity protein kinases that phosphorylate Serine/Arginine-rich (SR) proteins involved in pre-mRNA processing. Four CLKs, termed PfCLK-1-4, can be identified in the human malaria parasite Plasmodium falciparum, which show homology with the yeast SR protein kinase Sky1p. The four PfCLKs are present in the nucleus and cytoplasm of the asexual blood stages and of gametocytes, sexual precursor cells crucial for malaria parasite transmission from humans to mosquitoes. We identified three plasmodial SR proteins, PfSRSF12, PfSFRS4 and PfSF-1, which are predominantly present in the nucleus of blood stage trophozoites, PfSRSF12 and PfSF-1 are further detectable in the nucleus of gametocytes. We found that recombinantly expressed SR proteins comprising the Arginine/Serine (RS)-rich domains were phosphorylated by the four PfCLKs in in vitro kinase assays, while a recombinant PfSF-1 peptide lacking the RS-rich domain was not phosphorylated. Since it was hitherto not possible to knock-out the pfclk genes by conventional gene disruption, we aimed at chemical knock-outs for phenotype analysis. We identified five human CLK inhibitors, belonging to the oxo-beta-carbolines and aminopyrimidines, as well as the antiseptic chlorhexidine as PfCLK-targeting compounds. The six inhibitors block P. falciparum blood stage replication in the low micromolar to nanomolar range by preventing the trophozoite-to-schizont transformation. In addition, the inhibitors impair gametocyte maturation and gametogenesis in in vitro assays. The combined data show that the four PfCLKs are involved in phosphorylation of SR proteins with essential functions for the blood and sexual stages of the malaria parasite, thus pointing to the kinases as promising targets for antimalarial and transmission blocking drugs.
Two-component systems (TCS) are short signalling pathways generally occurring in prokaryotes. They frequently regulate prokaryotic stimulus responses and thus are also of interest for engineering in biotechnology and synthetic biology. The aim of this study is to better understand and describe rewiring of TCS while investigating different evolutionary scenarios. Based on large-scale screens of TCS in different organisms, this study gives detailed data, concrete alignments, and structure analysis on three general modification scenarios, where TCS were rewired for new responses and functions: (i) exchanges in the sequence within single TCS domains, (ii) exchange of whole TCS domains; (iii) addition of new components modulating TCS function. As a result, the replacement of stimulus and promotor cassettes to rewire TCS is well defined exploiting the alignments given here. The diverged TCS examples are non-trivial and the design is challenging. Designed connector proteins may also be useful to modify TCS in selected cases.
MicroRNAs are well-known strong RNA regulators modulating whole functional units in complex signaling networks. Regarding clinical application, they have potential as biomarkers for prognosis, diagnosis, and therapy. In this review, we focus on two microRNAs centrally involved in lung cancer progression. MicroRNA-21 promotes and microRNA-34 inhibits cancer progression. We elucidate here involved pathways and imbed these antagonistic microRNAs in a network of interactions, stressing their cancer microRNA biology, followed by experimental and bioinformatics analysis of such microRNAs and their targets. This background is then illuminated from a clinical perspective on microRNA-21 and microRNA-34 as general examples for the complex microRNA biology in lung cancer and its diagnostic value. Moreover, we discuss the immense potential that microRNAs such as microRNA-21 and microRNA-34 imply by their broad regulatory effects. These should be explored for novel therapeutic strategies in the clinic.
The drug-minded protein interaction database (DrumPID) has been designed to provide fast, tailored information on drugs and their protein networks including indications, protein targets and side-targets. Starting queries include compound, target and protein interactions and organism-specific protein families. Furthermore, drug name, chemical structures and their SMILES notation, affected proteins (potential drug targets), organisms as well as diseases can be queried including various combinations and refinement of searches. Drugs and protein interactions are analyzed in detail with reference to protein structures and catalytic domains, related compound structures as well as potential targets in other organisms. DrumPID considers drug functionality, compound similarity, target structure, interactome analysis and organismic range for a compound, useful for drug development, predicting drug side-effects and structure–activity relationships.
Lung cancer is currently the leading cause of cancer related mortality due to late diagnosis and limited treatment intervention. Non-coding RNAs are not translated into proteins and have emerged as fundamental regulators of gene expression. Recent studies reported that microRNAs and long non-coding RNAs are involved in lung cancer development and progression. Moreover, they appear as new promising non-invasive biomarkers for early lung cancer diagnosis. Here, we highlight their potential as biomarker in lung cancer and present how bioinformatics can contribute to the development of non-invasive diagnostic tools. For this, we discuss several bioinformatics algorithms and software tools for a comprehensive understanding and functional characterization of microRNAs and long non-coding RNAs.
High attrition-rates entailed by drug testing in 2D cell culture and animal models stress the need for improved modeling of human tumor tissues. In previous studies our 3D models on a decellularized tissue matrix have shown better predictivity and higher chemoresistance. A single porcine intestine yields material for 150 3D models of breast, lung, colorectal cancer (CRC) or leukemia. The uniquely preserved structure of the basement membrane enables physiological anchorage of endothelial cells and epithelial-derived carcinoma cells. The matrix provides different niches for cell growth: on top as monolayer, in crypts as aggregates and within deeper layers. Dynamic culture in bioreactors enhances cell growth. Comparing gene expression between 2D and 3D cultures, we observed changes related to proliferation, apoptosis and stemness. For drug target predictions, we utilize tumor-specific sequencing data in our in silico model finding an additive effect of metformin and gefitinib treatment for lung cancer in silico, validated in vitro. To analyze mode-of-action, immune therapies such as trispecific T-cell engagers in leukemia, as well as toxicity on non-cancer cells, the model can be modularly enriched with human endothelial cells (hECs), immune cells and fibroblasts. Upon addition of hECs, transmigration of immune cells through the endothelial barrier can be investigated. In an allogenic CRC model we observe a lower basic apoptosis rate after applying PBMCs in 3D compared to 2D, which offers new options to mirror antigen-specific immunotherapies in vitro. In conclusion, we present modular human 3D tumor models with tissue-like features for preclinical testing to reduce animal experiments.
Whereas the role of calcium ions (Ca\(^{2+}\)) in plant signaling is well studied, the physiological significance of pH‐changes remains largely undefined.
Here we developed CapHensor, an optimized dual‐reporter for simultaneous Ca\(^{2+}\) and pH ratio‐imaging and studied signaling events in pollen tubes (PTs), guard cells (GCs), and mesophyll cells (MCs). Monitoring spatio‐temporal relationships between membrane voltage, Ca\(^{2+}\)‐ and pH‐dynamics revealed interconnections previously not described.
In tobacco PTs, we demonstrated Ca\(^{2+}\)‐dynamics lag behind pH‐dynamics during oscillatory growth, and pH correlates more with growth than Ca\(^{2+}\). In GCs, we demonstrated abscisic acid (ABA) to initiate stomatal closure via rapid cytosolic alkalization followed by Ca2+ elevation. Preventing the alkalization blocked GC ABA‐responses and even opened stomata in the presence of ABA, disclosing an important pH‐dependent GC signaling node. In MCs, a flg22‐induced membrane depolarization preceded Ca2+‐increases and cytosolic acidification by c. 2 min, suggesting a Ca\(^{2+}\)/pH‐independent early pathogen signaling step. Imaging Ca2+ and pH resolved similar cytosol and nuclear signals and demonstrated flg22, but not ABA and hydrogen peroxide to initiate rapid membrane voltage‐, Ca\(^{2+}\)‐ and pH‐responses.
We propose close interrelation in Ca\(^{2+}\)‐ and pH‐signaling that is cell type‐ and stimulus‐specific and the pH having crucial roles in regulating PT growth and stomata movement.
We observed substantial differences in predicted Major Histocompatibility Complex II (MHCII) epitope presentation of SARS-CoV-2 proteins for different populations but only minor differences in predicted MHCI epitope presentation. A comparison of this predicted epitope MHC-coverage revealed for the early phase of infection spread (till day 15 after reaching 128 observed infection cases) highly significant negative correlations with the case fatality rate. Specifically, this was observed in different populations for MHC class II presentation of the viral spike protein (p-value: 0.0733 for linear regression), the envelope protein (p-value: 0.023), and the membrane protein (p-value: 0.00053), indicating that the high case fatality rates of COVID-19 observed in some countries seem to be related with poor MHC class II presentation and hence weak adaptive immune response against these viral envelope proteins. Our results highlight the general importance of the SARS-CoV-2 structural proteins in immunological control in early infection spread looking at a global census in various countries and taking case fatality rate into account. Other factors such as health system and control measures become more important after the early spread. Our study should encourage further studies on MHCII alleles as potential risk factors in COVID-19 including assessment of local populations and specific allele distributions.
Serine/threonine kinase PknB and its corresponding phosphatase Stp are important regulators of many cell functions in the pathogen S. aureus. Genome-scale gene expression data of S. aureus strain NewHG (sigB\(^+\)) elucidated their effect on physiological functions. Moreover, metabolic modelling from these data inferred metabolic adaptations. We compared wild-type to deletion strains lacking pknB, stp or both. Ser/Thr phosphorylation of target proteins by PknB switched amino acid catabolism off and gluconeogenesis on to provide the cell with sufficient components. We revealed a significant impact of PknB and Stp on peptidoglycan, nucleotide and aromatic amino acid synthesis, as well as catabolism involving aspartate transaminase. Moreover, pyrimidine synthesis was dramatically impaired by stp deletion but only slightly by functional loss of PknB. In double knockouts, higher activity concerned genes involved in peptidoglycan, purine and aromatic amino acid synthesis from glucose but lower activity of pyrimidine synthesis from glucose compared to the wild type. A second transcriptome dataset from S. aureus NCTC 8325 (sigB\(^−\)) validated the predictions. For this metabolic adaptation, PknB was found to interact with CdaA and the yvcK/glmR regulon. The involved GlmR structure and the GlmS riboswitch were modelled. Furthermore, PknB phosphorylation lowered the expression of many virulence factors, and the study shed light on S. aureus infection processes.