Transcriptomic meta-analysis to identify potential antifungal targets in Candida albicans.

BACKGROUND: Candida albicans is a fungal pathogen causing human infections. Here we investigated differential gene expression patterns and functional enrichment in C. albicans strains grown under different conditions. METHODS: A systematic GEO database search identified 239 "Candida albicans" datasets, of which 14 were selected after rigorous criteria application. Retrieval of raw sequencing data from the ENA database was accompanied by essential metadata extraction from dataset descriptions and original articles. Pre-processing via the tailored nf-core pipeline for C. albicans involved alignment, gene/transcript quantification, and diverse quality control measures. Quality assessment via PCA and DESeq2 identified significant genes (FDR < = 0.05, log2-fold change > = 1 or <= -1), while topGO conducted GO term enrichment analysis. Exclusions were made based on data quality and strain relevance, resulting in the selection of seven datasets from the SC5314 strain background for in-depth investigation. RESULTS: The meta-analysis of seven selected studies unveiled a substantial number of genes exhibiting significant up-regulation (24,689) and down-regulation (18,074). These differentially expressed genes were further categorized into 2,497 significantly up-regulated and 2,573 significantly down-regulated Gene Ontology (GO) IDs. GO term enrichment analysis clustered these terms into distinct groups, providing insights into the functional implications. Three target gene lists were compiled based on previous studies, focusing on central metabolism, ion homeostasis, and pathogenicity. Frequency analysis revealed genes with higher occurrence within the identified GO clusters, suggesting their potential as antifungal targets. Notably, the genes TPS2, TPS1, RIM21, PRA1, SAP4, and SAP6 exhibited higher frequencies within the clusters. Through frequency analysis within the GO clusters, several key genes emerged as potential targets for antifungal therapies. These include RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101 which exhibited higher occurrence within the identified clusters. CONCLUSION: This comprehensive study significantly advances our understanding of the dynamic nature of gene expression in C. albicans. The identification of genes with enhanced potential as antifungal drug targets underpins their value for future interventions. The highlighted genes, including TPS2, TPS1, RIM21, PRA1, SAP4, SAP6, RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101, hold promise for the development of targeted antifungal therapies.

Improving recombinant protein production by yeast through genome-scale modeling using proteome constraints.

Eukaryotic cells are used as cell factories to produce and secrete multitudes of recombinant pharmaceutical proteins, including several of the current top-selling drugs. Due to the essential role and complexity of the secretory pathway, improvement for recombinant protein production through metabolic engineering has traditionally been relatively ad-hoc; and a more systematic approach is required to generate novel design principles. Here, we present the proteome-constrained genome-scale protein secretory model of yeast Saccharomyces cerevisiae (pcSecYeast), which enables us to simulate and explain phenotypes caused by limited secretory capacity. We further apply the pcSecYeast model to predict overexpression targets for the production of several recombinant proteins. We experimentally validate many of the predicted targets for α-amylase production to demonstrate pcSecYeast application as a computational tool in guiding yeast engineering and improving recombinant protein production.

Whole-cell modeling in yeast predicts compartment-specific proteome constraints that drive metabolic strategies.

When conditions change, unicellular organisms rewire their metabolism to sustain cell maintenance and cellular growth. Such rewiring may be understood as resource re-allocation under cellular constraints. Eukaryal cells contain metabolically active organelles such as mitochondria, competing for cytosolic space and resources, and the nature of the relevant cellular constraints remain to be determined for such cells. Here, we present a comprehensive metabolic model of the yeast cell, based on its full metabolic reaction network extended with protein synthesis and degradation reactions. The model predicts metabolic fluxes and corresponding protein expression by constraining compartment-specific protein pools and maximising growth rate. Comparing model predictions with quantitative experimental data suggests that under glucose limitation, a mitochondrial constraint limits growth at the onset of ethanol formation-known as the Crabtree effect. Under sugar excess, however, a constraint on total cytosolic volume dictates overflow metabolism. Our comprehensive model thus identifies condition-dependent and compartment-specific constraints that can explain metabolic strategies and protein expression profiles from growth rate optimisation, providing a framework to understand metabolic adaptation in eukaryal cells.

Constraint-based models for dominating protein interaction networks.

The minimum dominating set (MDSet) comprises the smallest number of graph nodes, where other graph nodes are connected with at least one MDSet node. The MDSet has been successfully applied to extract proteins that control protein-protein interaction (PPI) networks and to reveal the correlation between structural analysis and biological functions. Although the PPI network contains many MDSets, the identification of multiple MDSets is an NP-complete problem, and it is difficult to determine the best MDSets, enriched with biological functions. Therefore, the MDSet model needs to be further expanded and validated to find constrained solutions that differ from those generated by the traditional models. Moreover, by identifying the critical set of the network, the set of nodes common to all MDSets can be time-consuming. Herein, the authors adopted the minimisation of metabolic adjustment (MOMA) algorithm to develop a new framework, called maximisation of interaction adjustment (MOIA). In MOIA, they provide three models; the first one generates two MDSets with a minimum number of shared proteins, the second model generates constrained multiple MDSets ( k -MDSets), and the third model generates user-defined MDSets, containing the maximum number of essential genes and/or other important genes of the PPI network. In practice, these models significantly reduce the cost of finding the critical set and classifying the graph nodes. Herein, the authors termed the critical set as the k -critical set, where k is the number of MDSets generated by the proposed model. Then, they defined a new set of proteins called the ( k - 1 ) -critical set, where each node belongs to ( k - 1 ) MDSets. This set has been shown to be as important as the k -critical set and contains many essential genes, transcription factors, and protein kinases as the k -critical set. The ( k - 1 ) -critical set can be used to extend the search for drug target proteins. Based on the performance of the MOIA models, the authors believe the proposed methods contribute to answering key questions about the MDSets of PPI networks, and their results and analysis can be extended to other network types.

Absolute Quantification of Protein and mRNA Abundances Demonstrate Variability in Gene-Specific Translation Efficiency in Yeast.

Protein synthesis is the most energy-consuming process in a proliferating cell, and understanding what controls protein abundances represents a key question in biology and biotechnology. We quantified absolute abundances of 5,354 mRNAs and 2,198 proteins in Saccharomyces cerevisiae under ten environmental conditions and protein turnover for 1,384 proteins under a reference condition. The overall correlation between mRNA and protein abundances across all conditions was low (0.46), but for differentially expressed proteins (n = 202), the median mRNA-protein correlation was 0.88. We used these data to model translation efficiencies and found that they vary more than 400-fold between genes. Non-linear regression analysis detected that mRNA abundance and translation elongation were the dominant factors controlling protein synthesis, explaining 61% and 15% of its variance. Metabolic flux balance analysis further showed that only mitochondrial fluxes were positively associated with changes at the transcript level. The present dataset represents a crucial expansion to the current resources for future studies on yeast physiology.

Systems biology analysis of hepatitis C virus infection reveals the role of copy number increases in regions of chromosome 1q in hepatocellular carcinoma metabolism.

Hepatitis C virus (HCV) infection is a worldwide healthcare problem; however, traditional treatment methods have failed to cure all patients, and HCV has developed resistance to new drugs. Systems biology-based analyses could play an important role in the holistic analysis of the impact of HCV on hepatocellular metabolism. Here, we integrated HCV assembly reactions with a genome-scale hepatocyte metabolic model to identify metabolic targets for HCV assembly and metabolic alterations that occur between different HCV progression states (cirrhosis, dysplastic nodule, and early and advanced hepatocellular carcinoma (HCC)) and healthy liver tissue. We found that diacylglycerolipids were essential for HCV assembly. In addition, the metabolism of keratan sulfate and chondroitin sulfate was significantly changed in the cirrhosis stage, whereas the metabolism of acyl-carnitine was significantly changed in the dysplastic nodule and early HCC stages. Our results explained the role of the upregulated expression of BCAT1, PLOD3 and six other methyltransferase genes involved in carnitine biosynthesis and S-adenosylmethionine metabolism in the early and advanced HCC stages. Moreover, GNPAT and BCAP31 expression was upregulated in the early and advanced HCC stages and could lead to increased acyl-CoA consumption. By integrating our results with copy number variation analyses, we observed that GNPAT, PPOX and five of the methyltransferase genes (ASH1L, METTL13, SMYD2, TARBP1 and SMYD3), which are all located on chromosome 1q, had increased copy numbers in the cancer samples relative to the normal samples. Finally, we confirmed our predictions with the results of metabolomics studies and proposed that inhibiting the identified targets has the potential to provide an effective treatment strategy for HCV-associated liver disorders.