Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling

Modelling the effect of SPION size in a stent assisted magnetic drug targeting system with interparticle interactions

Cancer is a leading cause of death worldwide and it is caused by the interaction of genomic, environmental, and lifestyle factors. Although chemotherapy is one way of treating cancers, it also damages healthy cells and may cause severe side effects. Therefore, it is beneficial in drug delivery in the human body to increase the proportion of the drugs at the target site while limiting its exposure at the rest of body through Magnetic Drug Targeting (MDT). Superparamagnetic iron oxide nanoparticles (SPIONs) are derived from polyol methods and coated with oleic acid and can be used as magnetic drug carrier particles (MDCPs) in an MDT system. Here, we develop a mathematical model for studying the interactions between the MDCPs enriched with three different diameters of SPIONs (6.6, 11.6, and 17.8 nm) in the MDT system with an implanted magnetizable stent using different magnetic field strengths and blood velocities. Our computational analysis allows for the optimal design of the SPIONs enriched MDCPs to be used in clinical applications.

Identifying anti-growth factors for human cancer cell lines through genome-scale metabolic modeling

Human cancer cell lines are used as important model systems to study molecular mechanisms associated with tumor growth, hereunder how genomic and biological heterogeneity found in primary tumors affect cellular phenotypes. We reconstructed Genome scale metabolic models (GEMs) for eleven cell lines based on RNA-Seq data and validated the functionality of these models with data from metabolite profiling. We used cell line-specific GEMs to analyze the differences in the metabolism of cancer cell lines, and to explore the heterogeneous expression of the metabolic subsystems. Furthermore, we predicted 85 antimetabolites that can inhibit growth of, or even kill, any of the cell lines, while at the same time not being toxic for 83 different healthy human cell types. 60 of these antimetabolites were found to inhibit growth in all cell lines. Finally, we experimentally validated one of the predicted antimetabolites using two cell lines with different phenotypic origins, and found that it is effective in inhibiting the growth of these cell lines. Using immunohistochemistry, we also showed high or moderate expression levels of proteins targeted by the validated antimetabolite. Identified anti-growth factors for inhibition of cell growth may provide leads for the development of efficient cancer treatment strategies.

Human genes with a greater number of transcript variants tend to show biological features of housekeeping and essential genes

Human genes with a greater number of transcript variants are more likely to play functionally important roles such as cellular maintenance and survival.

Phenotype-based cell-specific metabolic modeling reveals metabolic liabilities of cancer

Utilizing molecular data to derive functional physiological models tailored for specific cancer cells can facilitate the use of individually tailored therapies. To this end we present an approach termed PRIME for generating cell-specific genome-scale metabolic models (GSMMs) based on molecular and phenotypic data. We build >280 models of normal and cancer cell-lines that successfully predict metabolic phenotypes in an individual manner. We utilize this set of cell-specific models to predict drug targets that selectively inhibit cancerous but not normal cell proliferation. The top predicted target, MLYCD, is experimentally validated and the metabolic effects of MLYCD depletion investigated. Furthermore, we tested cell-specific predicted responses to the inhibition of metabolic enzymes, and successfully inferred the prognosis of cancer patients based on their PRIME-derived individual GSMMs. These results lay a computational basis and a counterpart experimental proof of concept for future personalized metabolic modeling applications, enhancing the search for novel selective anticancer therapies.

Pathogenesis and prevention of hepatitis C virus-induced hepatocellular carcinoma

Hepatitis C virus (HCV) is one of the major aetiologic agents that causes hepatocellular carcinoma (HCC) by generating an inflammatory, fibrogenic, and carcinogenic tissue microenvironment in the liver. HCV-induced HCC is a rational target for cancer preventive intervention because of the clear-cut high-risk condition, cirrhosis, associated with high cancer incidence (1% to 7% per year). Studies have elucidated direct and indirect carcinogenic effects of HCV, which have in turn led to the identification of candidate HCC chemoprevention targets. Selective molecular targeted agents may enable personalized strategies for HCC chemoprevention. In addition, multiple experimental and epidemiological studies suggest the potential value of generic drugs or dietary supplements targeting inflammation, oxidant stress, or metabolic derangements as possible HCC chemopreventive agents. While the successful use of highly effective direct-acting antiviral agents will make important inroads into reducing long-term HCC risk, there will remain an important role for HCC chemoprevention even after viral cure, given the persistence of HCC risk in persons with advanced HCV fibrosis, as shown in recent studies. The successful development of cancer preventive therapies will be more challenging compared to cancer therapeutics because of the requirement for larger and longer clinical trials and the need for a safer toxicity profile given its use as a preventive agent. Molecular biomarkers to selectively identify high-risk population could help mitigate these challenges. Genome-wide, unbiased molecular characterization, high-throughput drug/gene screening, experimental model-based functional analysis, and systems-level in silico modelling are expected to complement each other to facilitate discovery of new HCC chemoprevention targets and therapies.

Defining the human gallbladder proteome by transcriptomics and affinity proteomics

Global protein analysis of human gallbladder tissue is vital for identification of molecular regulators and effectors of its physiological activity. Here, we employed a genome-wide deep RNA sequencing analysis in 28 human tissues to identify the genes overrepresented in the gallbladder and complemented it with antibody-based immunohistochemistry in 48 human tissues. We characterized human gallbladder proteins and identified 140 gallbladder-specific proteins with an elevated expression in the gallbladder as compared to the other analyzed tissues. Five genes were categorized as enriched, with at least fivefold higher levels in gallbladder, 60 genes were categorized as group enriched with elevated transcript levels in gallbladder shared with at least one other tissue and 75 genes were categorized as enhanced with higher expression than the average expression in other tissues. We explored the localization of the genes within the gallbladder through cell-type specific antibody-based protein profiling and the subcellular localization of the genes through immunofluorescent-based profiling. Finally, we revealed the biological processes and metabolic functions carried out by these genes through the use of GO, KEGG Pathway, and HMR2.0 that is compilation of the human metabolic reactions. We demonstrated the results of the combined analysis of the transcriptomics and affinity proteomics.

Defining the human adipose tissue proteome to reveal metabolic alterations in obesity

White adipose tissue (WAT) has a major role in the progression of obesity. Here, we combined data from RNA-Seq and antibody-based immunohistochemistry to describe the normal physiology of human WAT obtained from three female subjects and explored WAT-specific genes by comparing WAT to 26 other major human tissues. Using the protein evidence in WAT, we validated the content of a genome-scale metabolic model for adipocytes. We employed this high-quality model for the analysis of subcutaneous adipose tissue (SAT) gene expression data obtained from subjects included in the Swedish Obese Subjects Sib Pair study to reveal molecular differences between lean and obese individuals. We integrated SAT gene expression and plasma metabolomics data, investigated the contribution of the metabolic differences in the mitochondria of SAT to the occurrence of obesity, and eventually identified cytosolic branched-chain amino acid (BCAA) transaminase 1 as a potential target that can be used for drug development. We observed decreased glutaminolysis and alterations in the BCAAs metabolism in SAT of obese subjects compared to lean subjects. We also provided mechanistic explanations for the changes in the plasma level of BCAAs, glutamate, pyruvate, and α-ketoglutarate in obese subjects. Finally, we validated a subset of our model-based predictions in 20 SAT samples obtained from 10 lean and 10 obese male and female subjects.

A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration

Over the last decade, the field of cancer metabolism has mainly focused on studying the role of tumorigenic metabolic rewiring in supporting cancer proliferation. Here, we perform the first genome-scale computational study of the metabolic underpinnings of cancer migration. We build genome-scale metabolic models of the NCI-60 cell lines that capture the Warburg effect (aerobic glycolysis) typically occurring in cancer cells. The extent of the Warburg effect in each of these cell line models is quantified by the ratio of glycolytic to oxidative ATP flux (AFR), which is found to be highly positively associated with cancer cell migration. We hence predicted that targeting genes that mitigate the Warburg effect by reducing the AFR may specifically inhibit cancer migration. By testing the anti-migratory effects of silencing such 17 top predicted genes in four breast and lung cancer cell lines, we find that up to 13 of these novel predictions significantly attenuate cell migration either in all or one cell line only, while having almost no effect on cell proliferation. Furthermore, in accordance with the predictions, a significant reduction is observed in the ratio between experimentally measured ECAR and OCR levels following these perturbations. Inhibiting anti-migratory targets is a promising future avenue in treating cancer since it may decrease cytotoxic-related side effects that plague current anti-proliferative treatments. Furthermore, it may reduce cytotoxic-related clonal selection of more aggressive cancer cells and the likelihood of emerging resistance.

Hepatocellular carcinoma: genome-scale metabolic models for hepatocellular carcinoma

A new study proposes a modelling strategy to identify reactions, genes and metabolites relevant in hepatocellular carcinoma using in silico and in vivo analyses. The proposed genome-scale metabolic model integrates genomic and proteomic information, and points to statins, among others, as potential chemopreventive and anticancer drugs.

Elucidating the interactions between the human gut microbiota and its host through metabolic modeling

Increased understanding of the interactions between the gut microbiota, diet and environmental effects may allow us to design efficient treatment strategies for addressing global health problems. Existence of symbiotic microorganisms in the human gut provides different functions for the host such as conversion of nutrients, training of the immune system, and resistance to pathogens. The gut microbiome also plays an influential role in maintaining human health, and it is a potential target for prevention and treatment of common disorders including obesity, type 2 diabetes, and atherosclerosis. Due to the extreme complexity of such disorders, it is necessary to develop mathematical models for deciphering the role of its individual elements as well as the entire system and such models may assist in better understanding of the interactions between the bacteria in the human gut and the host by use of genome-scale metabolic models (GEMs). Recently, GEMs have been employed to explore the interactions between predominant bacteria in the gut ecosystems. Additionally, these models enabled analysis of the contribution of each species to the overall metabolism of the microbiota through the integration of omics data. The outcome of these studies can be used for proposing optimal conditions for desired microbiome phenotypes. Here, we review the recent progress and challenges for elucidating the interactions between the human gut microbiota and host through metabolic modeling. We discuss how these models may provide scaffolds for analyzing high-throughput data, developing probiotics and prebiotics, evaluating the effects of probiotics and prebiotics and eventually designing clinical interventions.

The human liver-specific proteome defined by transcriptomics and antibody-based profiling

Human liver physiology and the genetic etiology of the liver diseases can potentially be elucidated through the identification of proteins with enriched expression in the liver. Here, we combined data from RNA sequencing (RNA-Seq) and antibody-based immunohistochemistry across all major human tissues to explore the human liver proteome with enriched expression, as well as the cell type-enriched expression in hepatocyte and bile duct cells. We identified in total 477 protein-coding genes with elevated expression in the liver: 179 genes have higher expression as compared to all the other analyzed tissues; 164 genes have elevated transcript levels in the liver shared with at least one other tissue type; and an additional 134 genes have a mild level of increased expression in the liver. We identified the precise localization of these proteins through antibody-based protein profiling and the subcellular localization of these proteins through immunofluorescent-based profiling. We also identified the biological processes and metabolic functions associated with these proteins, investigated their contribution in the occurrence of liver diseases, and identified potential targets for their treatment. Our study demonstrates the use of RNA-Seq and antibody-based immunohistochemistry for characterizing the human liver proteome, as well as the use of tissue-specific proteins in identification of novel drug targets and discovery of biomarkers.-Kampf, C., Mardinoglu, A., Fagerberg, L., Hallström, B. M., Edlund, K., Lundberg, E., Pontén, F., Nielsen, J., Uhlen, M. The human liver-specific proteome defined by transcriptomics and antibody-based profiling.

Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling

Genome-scale metabolic models (GEMs) have proven useful as scaffolds for the integration of omics data for understanding the genotype-phenotype relationship in a mechanistic manner. Here, we evaluated the presence/absence of proteins encoded by 15,841 genes in 27 hepatocellular carcinoma (HCC) patients using immunohistochemistry. We used this information to reconstruct personalized GEMs for six HCC patients based on the proteomics data, HMR 2.0, and a task-driven model reconstruction algorithm (tINIT). The personalized GEMs were employed to identify anticancer drugs using the concept of antimetabolites; i.e., drugs that are structural analogs to metabolites. The toxicity of each antimetabolite was predicted by assessing the in silico functionality of 83 healthy cell type-specific GEMs, which were also reconstructed with the tINIT algorithm. We predicted 101 antimetabolites that could be effective in preventing tumor growth in all HCC patients, and 46 antimetabolites which were specific to individual patients. Twenty-two of the 101 predicted antimetabolites have already been used in different cancer treatment strategies, while the remaining antimetabolites represent new potential drugs. Finally, one of the identified targets was validated experimentally, and it was confirmed to attenuate growth of the HepG2 cell line.

Generalized framework for context-specific metabolic model extraction methods

Genome-scale metabolic models (GEMs) are increasingly applied to investigate the physiology not only of simple prokaryotes, but also eukaryotes, such as plants, characterized with compartmentalized cells of multiple types. While genome-scale models aim at including the entirety of known metabolic reactions, mounting evidence has indicated that only a subset of these reactions is active in a given context, including: developmental stage, cell type, or environment. As a result, several methods have been proposed to reconstruct context-specific models from existing genome-scale models by integrating various types of high-throughput data. Here we present a mathematical framework that puts all existing methods under one umbrella and provides the means to better understand their functioning, highlight similarities and differences, and to help users in selecting a most suitable method for an application.