MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEv​al/downloads

A predicted physicochemically distinct sub-proteome associated with the intracellular organelle of the anammox bacterium Kuenenia stuttgartiensis

Background: Anaerobic ammonium-oxidizing (anammox) bacteria perform a key step in global nitrogen cycling. These bacteria make use of an organelle to oxidize ammonia anaerobically to nitrogen (N2) and so contribute approximately 50% of the nitrogen in the atmosphere. It is currently unknown which proteins constitute the organellar proteome and how anammox bacteria are able to specifically target organellar and cell-envelope proteins to their correct final destinations. Experimental approaches are complicated by the absence of pure cultures and genetic accessibility. However, the genome of the anammox bacterium Candidatus "Kuenenia stuttgartiensis" has recently been sequenced. Here, we make use of these genome data to predict the organellar sub-proteome and address the molecular basis of protein sorting in anammox bacteria. Results: Two training sets representing organellar (30 proteins) and cell envelope (59 proteins) proteins were constructed based on previous experimental evidence and comparative genomics. Random forest (RF) classifiers trained on these two sets could differentiate between organellar and cell envelope proteins with ~89% accuracy using 400 features consisting of frequencies of two adjacent amino acid combinations. A physicochemically distinct organellar sub-proteome containing 562 proteins was predicted with the best RF classifier. This set included almost all catabolic and respiratory factors encoded in the genome. Apparently, the cytoplasmic membrane performs no catabolic functions. We predict that the Tat-translocation system is located exclusively in the organellar membrane, whereas the Sec-translocation system is located on both the organellar and cytoplasmic membranes. Canonical signal peptides were predicted and validated experimentally, but a specific (N- or C-terminal) signal that could be used for protein targeting to the organelle remained elusive. Conclusions: A physicochemically distinct organellar sub-proteome was predicted from the genome of the anammox bacterium K. stuttgartiensis. This result provides strong in silico support for the existing experimental evidence for the existence of an organelle in this bacterium, and is an important step forward in unravelling a geochemically relevant case of cytoplasmic differentiation in bacteria. The predicted dual location of the Sec-translocation system and the apparent absence of a specific N- or C-terminal signal in the organellar proteins suggests that additional chaperones may be necessary that act on an as-yet unknown property of the targeted proteins

Перспективи розвитку експортоорієнтованої стратегії підприємств

Рассмотрен вопрос стратегического развития экспортноориентрованной политики предприятий. Раскрыты перспективы развития международных торговых отношений Украины.Розглянуто питання стратегічного розвитку експортноорієнтовної політики підприємств. Розкрито перспективи розвитку міжнародних торгівельних відносин України

RNF168 cooperates with RNF8 to mediate FOXM1 ubiquitination and degradation in breast cancer epirubicin treatment

The forkhead box M1 (FOXM1) transcription factor has a central role in genotoxic agent response in breast cancer. FOXM1 is regulated at the post-translational level upon DNA damage, but the key mechanism involved remained enigmatic. RNF168 is a ubiquitination E3-ligase involved in DNA damage response. Western blot and gene promoter-reporter analyses showed that the expression level and transcriptional activity of FOXM1 reduced upon RNF168 overexpression and increased with RNF168 depletion by siRNA, suggesting that RNF168 negatively regulates FOXM1 expression. Co-immunoprecipitation studies in MCF-7 cells revealed that RNF168 interacted with FOXM1 and that upon epirubicin treatment FOXM1 downregulation was associated with an increase in RNF168 binding and conjugation to the protein degradation-associated K48-linked polyubiquitin chains. Consistently, RNF168 overexpression resulted in an increase in turnover of FOXM1 in MCF-7 cells treated with the protein synthesis inhibitor cycloheximide. Conversely, RNF168, knockdown significantly enhanced the half-life of FOXM1 in both absence and presence of epirubicin. Using a SUMOylation-defective FOXM1-5x(K>R) mutant, we demonstrated that SUMOylation is required for the recruitment of RNF168 to mediate FOXM1 degradation. In addition, clonogenic assays also showed that RNF168 mediates epirubicin action through targeting FOXM1, as RNF168 could synergise with epirubicin to repress clonal formation in wild-type but not in FOXM1-deficient mouse embryo fibroblasts (MEFs). The physiological relevance of RNF168-mediated FOXM1 repression is further emphasized by the significant inverse correlation between FOXM1 and RNF168 expression in breast cancer patient samples. Moreover, we also obtained evidence that RNF8 recruits RNF168 to FOXM1 upon epirubicin treatment and cooperates with RNF168 to catalyse FOXM1 ubiquitination and degradation. Collectively, these data suggest that RNF168 cooperates with RNF8 to mediate the ubiquitination and degradation of SUMOylated FOXM1 in breast cancer genotoxic response.published_or_final_versio

OTUB1 inhibits the ubiquitination and degradation of FOXM1 in breast cancer and epirubicin resistance

The forkhead transcription factor FOXM1 has a key role in DNA damage response, and its deregulated overexpression is associated with genotoxic drug resistance in breast cancer. However, little is known about the posttranslational mechanisms by which FOXM1 expression is regulated by genotoxic agents and how they are deregulated in resistant cells. Initial co-immunoprecipitation studies verified previous proteomic analysis finding that the OTUB1 is a novel FOXM1-interacting protein. Western blot analysis showed that both OTUB1 and FOXM1 expression reduced upon genotoxic agent treatment in MCF-7 cells, but remained relatively constant in resistant cells. FOXM1 expression reduced upon OTUB1 depletion by siRNA and increased with OTUB1 overexpression in MCF-7 cells, arguing that OTUB1 positively regulates FOXM1 expression. In agreement, co-immunoprecipitation experiments demonstrated that FOXM1 expression is associated with OTUB1 binding but inversely correlates with conjugation to the protein degradation-associated Lys-48-linked ubiquitin-chains. Overexpression of wild-type (WT) OTUB1, but not the OTUB1(C91S) mutant, disrupted the formation of Lys48-linked ubiquitin-conjugates on FOXM1. Importantly, knockdown of OTUB1 by siRNA resulted in an increase in turnover of FOXM1 in MCF-7 cells treated with the protein synthesis inhibitor cycloheximide, whereas overexpression of WT OTUB1, but not the OTUB1(C91S) mutant, significantly enhances the half-life of FOXM1. In addition, proliferative and clonogenic assays also show that OTUB1 can enhance the proliferative rate and epirubicin resistance through targeting FOXM1, as OTUB1 has little effect on FOXM1-deficient cells. The physiological relevance of the regulation of FOXM1 by OTUB1 is further underscored by the significant correlations between FOXM1 and OTUB1 expression in breast cancer patient samples. Cox-regression survival analysis indicates that OTUB1 overexpression is linked to poorer outcome in particular in patients treated with chemotherapy. Collectively, these data suggest that OTUB1 limits the ubiquitination and degradation of FOXM1 in breast cancer and has a key role in genotoxic agent resistance

Indexing the Pseudomonas specialized metabolome enabled the discovery of poaeamide B and the bananamides

Pseudomonads are cosmopolitan microorganisms able to produce a wide array of specialized metabolites. These molecules allow Pseudomonas to scavenge nutrients, sense population density and enhance or inhibit growth of competing microorganisms. However, these valuable metabolites are typically characterized one-molecule–one-microbe at a time, instead of being inventoried in large numbers. To index and map the diversity of molecules detected from these organisms, 260 strains of ecologically diverse origins were subjected to mass-spectrometry-based molecular networking. Molecular networking not only enables dereplication of molecules, but also sheds light on their structural relationships. Moreover, it accelerates the discovery of new molecules. Here, by indexing the Pseudomonas specialized metabolome, we report the molecular-networking-based discovery of four molecules and their evolutionary relationships: a poaeamide analogue and a molecular subfamily of cyclic lipopeptides, bananamides 1, 2 and 3. Analysis of their biosynthetic gene cluster shows that it constitutes a distinct evolutionary branch of the Pseudomonas cyclic lipopeptides. Through analysis of an additional 370 extracts of wheat-associated Pseudomonas, we demonstrate how the detailed knowledge from our reference index can be efficiently propagated to annotate complex metabolomic data from other studies, akin to the way in which newly generated genomic information can be compared to data from public databases

Dissecting Disease-Suppressive Rhizosphere Microbiomes by Functional Amplicon Sequencing and 10x Metagenomics

Disease-suppressive soils protect plants against soilborne fungal pathogens that would otherwise cause root infections. Soil suppressiveness is, in most cases, mediated by the antagonistic activity of the microbial community associated with the plant roots. Considering the enormous taxonomic and functional diversity of the root-associated microbiome, identification of the microbial genera and mechanisms underlying this phenotype is challenging. One approach to unravel the underlying mechanisms is to identify metabolic pathways enriched in the disease-suppressive microbial community, in particular, pathways that harbor natural products with antifungal properties. An important class of these natural products includes peptides produced by nonribosomal peptide synthetases (NRPSs). Here, we applied functional amplicon sequencing of NRPS-associated adenylation domains (A domains) to a collection of eight soils that are suppressive or nonsuppressive (i.e., conducive) to Fusarium culmorum, a fungal root pathogen of wheat. To identify functional elements in the root-associated bacterial community, we developed an open-source pipeline, referred to as dom2BGC, for amplicon annotation and putative gene cluster reconstruction through analyzing A domain co-occurrence across samples. We applied this pipeline to rhizosphere communities from four disease-suppressive and four conducive soils and found significant similarities in NRPS repertoires between suppressive soils. Specifically, several siderophore biosynthetic gene clusters were consistently associated with suppressive soils, hinting at competition for iron as a potential mechanism of suppression. Finally, to validate dom2BGC and to allow more unbiased functional metagenomics, we performed 10× metagenomic sequencing of one suppressive soil, leading to the identification of multiple gene clusters potentially associated with the disease-suppressive phenotyp

The tomato cytochrome P450 CYP712G1 catalyzes the double oxidation of orobanchol <i>en route</i> to the rhizosphere signaling strigolactone, solanacol

Strigolactones (SLs) are rhizosphere signalling molecules and phytohormones. The biosynthetic pathway of SLs in tomato has been partially elucidated, but the structural diversity in tomato SLs predicts that additional biosynthetic steps are required. Here, root RNA-seq data and co-expression analysis were used for SL biosynthetic gene discovery. This strategy resulted in a candidate gene list containing several cytochrome P450s. Heterologous expression in Nicotiana benthamiana and yeast showed that one of these, CYP712G1, can catalyse the double oxidation of orobanchol, resulting in the formation of three didehydro-orobanchol (DDH) isomers. Virus-induced gene silencing and heterologous expression in yeast showed that one of these DDH isomers is converted to solanacol, one of the most abundant SLs in tomato root exudate. Protein modelling and substrate docking analysis suggest that hydroxy-orbanchol is the likely intermediate in the conversion from orobanchol to the DDH isomers. Phylogenetic analysis demonstrated the occurrence of CYP712G1 homologues in the Eudicots only, which fits with the reports on DDH isomers in that clade. Protein modelling and orobanchol docking of the putative tobacco CYP712G1 homologue suggest that it can convert orobanchol to similar DDH isomers as tomato