@article{KunzLiangNillaetal.2016, author = {Kunz, Meik and Liang, Chunguang and Nilla, Santosh and Cecil, Alexander and Dandekar, Thomas}, title = {The drug-minded protein interaction database (DrumPID) for efficient target analysis and drug development}, series = {Database}, volume = {2016}, journal = {Database}, doi = {10.1093/database/baw041}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-147369}, pages = {baw041}, year = {2016}, abstract = {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.}, language = {en} } @article{KaltdorfSrivastavaGuptaetal.2016, author = {Kaltdorf, Martin and Srivastava, Mugdha and Gupta, Shishir K. and Liang, Chunguang and Binder, Jasmin and Dietl, Anna-Maria and Meir, Zohar and Haas, Hubertus and Osherov, Nir and Krappmann, Sven and Dandekar, Thomas}, title = {Systematic Identification of Anti-Fungal Drug Targets by a Metabolic Network Approach}, series = {Frontiers in Molecular Bioscience}, volume = {3}, journal = {Frontiers in Molecular Bioscience}, doi = {10.3389/fmolb.2016.00022}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-147396}, pages = {22}, year = {2016}, abstract = {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.}, language = {en} } @article{FoersterBeisserGrohmeetal.2012, author = {F{\"o}rster, Frank and Beisser, Daniela and Grohme, Markus A. and Liang, Chunguang and Mali, Brahim and Siegl, Alexander Matthias and Engelmann, Julia C. and Shkumatov, Alexander V. and Schokraie, Elham and M{\"u}ller, Tobias and Schn{\"o}lzer, Martina and Schill, Ralph O. and Frohme, Marcus and Dandekar, Thomas}, title = {Transcriptome analysis in tardigrade species reveals specific molecular pathways for stress adaptations}, series = {Bioinformatics and biology insights}, volume = {6}, journal = {Bioinformatics and biology insights}, doi = {10.4137/BBI.S9150}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-123089}, pages = {69-96}, year = {2012}, abstract = {Tardigrades have unique stress-adaptations that allow them to survive extremes of cold, heat, radiation and vacuum. To study this, encoded protein clusters and pathways from an ongoing transcriptome study on the tardigrade \(Milnesium\) \(tardigradum\) were analyzed using bioinformatics tools and compared to expressed sequence tags (ESTs) from \(Hypsibius\) \(dujardini\), revealing major pathways involved in resistance against extreme environmental conditions. ESTs are available on the Tardigrade Workbench along with software and databank updates. Our analysis reveals that RNA stability motifs for \(M.\) \(tardigradum\) are different from typical motifs known from higher animals. \(M.\) \(tardigradum\) and \(H.\) \(dujardini\) protein clusters and conserved domains imply metabolic storage pathways for glycogen, glycolipids and specific secondary metabolism as well as stress response pathways (including heat shock proteins, bmh2, and specific repair pathways). Redox-, DNA-, stress- and protein protection pathways complement specific repair capabilities to achieve the strong robustness of \(M.\) \(tardigradum\). These pathways are partly conserved in other animals and their manipulation could boost stress adaptation even in human cells. However, the unique combination of resistance and repair pathways make tardigrades and \(M.\) \(tardigradum\) in particular so highly stress resistant.}, language = {en} } @article{DandekarLiangKrueger2013, author = {Dandekar, Thomas and Liang, Chunguang and Kr{\"u}ger, Beate}, title = {GoSynthetic database tool to analyse natural and engineered molecular processes}, series = {Database}, journal = {Database}, doi = {10.1093/database/bat043}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-97023}, year = {2013}, abstract = {An essential topic for synthetic biologists is to understand the structure and function of biological processes and involved proteins and plan experiments accordingly. Remarkable progress has been made in recent years towards this goal. However, efforts to collect and present all information on processes and functions are still cumbersome. The database tool GoSynthetic provides a new, simple and fast way to analyse biological processes applying a hierarchical database. Four different search modes are implemented. Furthermore, protein interaction data, cross-links to organism-specific databases (17 organisms including six model organisms and their interactions), COG/KOG, GO and IntAct are warehoused. The built in connection to technical and engineering terms enables a simple switching between biological concepts and concepts from engineering, electronics and synthetic biology. The current version of GoSynthetic covers more than one million processes, proteins, COGs and GOs. It is illustrated by various application examples probing process differences and designing modifications.}, language = {en} } @article{RatzkaFoersterLiangetal.2012, author = {Ratzka, Carolin and F{\"o}rster, Frank and Liang, Chunguang and Kupper, Maria and Dandekar, Thomas and Feldhaar, Heike and Gross, Roy}, title = {Molecular characterization of antimicrobial peptide genes of the carpenter ant Camponotus floridanus}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-75985}, year = {2012}, abstract = {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.}, subject = {Biologie}, language = {en} } @phdthesis{Liang2009, author = {Liang, Chunguang}, title = {Tools for functional genomics applied to Staphylococci, Listeriae, Vaccinia virus and other organisms}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-48051}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2009}, abstract = {Genome sequence analysis A combination of genome analysis application has been established here during this project. This offers an efficient platform to interactively compare similar genome regions and reveal loci differences. The genes and operons can be rapidly analyzed and local collinear blocks (LCBs) categorized according to their function. The features of interests are parsed, recognized, and clustered into reports. Phylogenetic relationships can be readily examined such as the evolution of critical factors or a certain highly-conserved region. The resulting platform-independent software packages (GENOVA and inGeno), have been proven to be efficient and easy to handle in a number of projects. The capabilities of the software allowed the investigation of virulence factors, e.g., rsbU, strains' biological design, and in particular pathogenicity feature storage and management. We have successfully investigated the genomes of Staphylococcus aureus strains (COL, N315, 8325, RN1HG, Newman), Listeria spp. (welshimeri, innocua and monocytogenes), E.coli strains (O157:H7 and MG1655) and Vaccinia strains (WR, Copenhagen, Lister, LIVP, GLV-1h68 and parental strains). Metabolic network analysis Our YANAsquare package offers a workbench to rapidly establish the metabolic network of such as Staphylococcous aureus bacteria in genome-scale size as well as metabolic networks of interest such as the murine phagosome lipid signalling network. YANAsquare recruits reactions from online databases using an integrated KEGG browser. This reduces the efforts in building large metabolic networks. The involved calculation routines (METATOOL-derived wrapper or native Java implementation) readily obtain all possible flux modes (EM/EP) for metabolite fluxes within the network. Advanced layout algorithms visualize the topological structure of the network. In addition, the generated structure can be dynamically modified in the graphic interface. The generated network as well as the manipulated layout can be validated and stored (XML file: scheme of SBML level-2). This format can be further parsed and analyzed by other systems biology software, such as CellDesigner. Moreover, the integrated robustness-evaluation routine is able to examine the synthesis rates affected by each single mutation throughout the whole network. We have successfully applied the method to simulate single and multiple gene knockouts, and the affected fluxes are comprehensively revealed. Recently we applied the method to proteomic data and extra-cellular metabolite data of Staphylococci, the physiological changes regarding the flux distribution are studied. Calculations at different time points, including different conditions such as hypoxia or stress, show a good fit to experimental data. Moreover, using the proteomic data (enzyme amounts) calculated from 2D-Gel-EP experiments our study provides a way to compare the fluxome and the enzyme expression. Oncolytic vaccinia virus (VACV) We investigated the genetic differences between the de novo sequence of the recombinant oncolytic GLV-1h68 and other related VACVs, including function predictions for all found genome differences. Our phylogenetic analysis indicates that GLV-1h68 is closest to Lister strains but has lost several ORFs present in its parental LIVP strain, including genes encoding CrmE and a viral Golgi anti-apoptotic protein, v-GAAP. Functions of viral genes were either strain-specific, tissue-specific or host-specific comparing viral genes in the Lister, WR and COP strains. This helps to rationally design more optimized oncolytic virus strains to benefit cancer therapy in human patients. Identified differences from the comparison in open reading frames (ORFs) include genes for host-range selection, virulence and immune modulation proteins, e.g. ankyrin-like proteins, serine proteinase inhibitor SPI-2/CrmA, tumor necrosis factor (TNF) receptor homolog CrmC, semaphorin-like and interleukin-1 receptor homolog proteins. The contribution of foreign gene expression cassettes in the therapeutic and oncolytic virus GLV-1h68 was studied, including the F14.5L, J2R and A56R loci. The contribution of F14.5L inactivation to the reduced virulence is demonstrated by comparing the virulence data of GLV-1h68 with its F14.5L-null and revertant viruses. The comparison suggests that insertion of a foreign gene expression cassette in a nonessential locus in the viral genome is a practical way to attenuate VACVs, especially if the nonessential locus itself contains a virulence gene. This reduces the virulence of the virus without compromising too much the replication competency of the virus, the key to its oncolytic activity. The reduced pathogenicity of GLV-1h68 was confirmed by our experimental collaboration partners in male mice bearing C6 rat glioma and in immunocompetent mice bearing B16-F10 murine melanoma. In conclusion, bioinformatics and experimental data show that GLV-1h68 is a promising engineered VACV variant for anticancer therapy with tumor-specific replication, reduced pathogenicity and benign tissue tropism.}, subject = {Genanalyse}, language = {en} } @article{KaltdorfBreitenbachKarletal.2023, author = {Kaltdorf, Martin and Breitenbach, Tim and Karl, Stefan and Fuchs, Maximilian and Kessie, David Komla and Psota, Eric and Prelog, Martina and Sarukhanyan, Edita and Ebert, Regina and Jakob, Franz and Dandekar, Gudrun and Naseem, Muhammad and Liang, Chunguang and Dandekar, Thomas}, title = {Software JimenaE allows efficient dynamic simulations of Boolean networks, centrality and system state analysis}, series = {Scientific Reports}, volume = {13}, journal = {Scientific Reports}, doi = {10.1038/s41598-022-27098-7}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-313303}, year = {2023}, abstract = {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.}, language = {en} } @article{LiangRiosMiguelJaricketal.2021, author = {Liang, Chunguang and Rios-Miguel, Ana B. and Jarick, Marcel and Neurgaonkar, Priya and Girard, Myriam and Fran{\c{c}}ois, Patrice and Schrenzel, Jacques and Ibrahim, Eslam S. and Ohlsen, Knut and Dandekar, Thomas}, title = {Staphylococcus aureus transcriptome data and metabolic modelling investigate the interplay of Ser/Thr kinase PknB, its phosphatase Stp, the glmR/yvcK regulon and the cdaA operon for metabolic adaptation}, series = {Microorganisms}, volume = {9}, journal = {Microorganisms}, number = {10}, issn = {2076-2607}, doi = {10.3390/microorganisms9102148}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-248459}, year = {2021}, abstract = {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.}, language = {en} } @article{AmpattuHagmannLiangetal.2017, author = {Ampattu, Biju Joseph and Hagmann, Laura and Liang, Chunguang and Dittrich, Marcus and Schl{\"u}ter, Andreas and Blom, Jochen and Krol, Elizaveta and Goesmann, Alexander and Becker, Anke and Dandekar, Thomas and M{\"u}ller, Tobias and Schoen, Christoph}, title = {Transcriptomic buffering of cryptic genetic variation contributes to meningococcal virulence}, series = {BMC Genomics}, volume = {18}, journal = {BMC Genomics}, number = {282}, doi = {10.1186/s12864-017-3616-7}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-157534}, year = {2017}, abstract = {Background: Commensal bacteria like Neisseria meningitidis sometimes cause serious disease. However, genomic comparison of hyperinvasive and apathogenic lineages did not reveal unambiguous hints towards indispensable virulence factors. Here, in a systems biological approach we compared gene expression of the invasive strain MC58 and the carriage strain α522 under different ex vivo conditions mimicking commensal and virulence compartments to assess the strain-specific impact of gene regulation on meningococcal virulence. Results: Despite indistinguishable ex vivo phenotypes, both strains differed in the expression of over 500 genes under infection mimicking conditions. These differences comprised in particular metabolic and information processing genes as well as genes known to be involved in host-damage such as the nitrite reductase and numerous LOS biosynthesis genes. A model based analysis of the transcriptomic differences in human blood suggested ensuing metabolic flux differences in energy, glutamine and cysteine metabolic pathways along with differences in the activation of the stringent response in both strains. In support of the computational findings, experimental analyses revealed differences in cysteine and glutamine auxotrophy in both strains as well as a strain and condition dependent essentiality of the (p)ppGpp synthetase gene relA and of a short non-coding AT-rich repeat element in its promoter region. Conclusions: Our data suggest that meningococcal virulence is linked to transcriptional buffering of cryptic genetic variation in metabolic genes including global stress responses. They further highlight the role of regulatory elements for bacterial virulence and the limitations of model strain approaches when studying such genetically diverse species as N. meningitidis.}, language = {en} } @article{PradaMaagSiegmundetal.2022, author = {Prada, Juan Pablo and Maag, Luca Estelle and Siegmund, Laura and Bencurova, Elena and Liang, Chunguang and Koutsilieri, Eleni and Dandekar, Thomas and Scheller, Carsten}, title = {Estimation of R0 for the spread of SARS-CoV-2 in Germany from excess mortality}, series = {Scientific Reports}, volume = {12}, journal = {Scientific Reports}, number = {1}, doi = {10.1038/s41598-022-22101-7}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-301415}, year = {2022}, abstract = {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.}, language = {en} }