- 1Biomedical Informatics Research Laboratory, Department of Biology, Syed Babar Ali School of Science and Engineering, Lahore University of Management Sciences, Lahore, Pakistan
- 2Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, United States
Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.
1 Introduction
In 1971, President Richard Nixon declared his euphemistic “war on cancer” through the promulgation of the National Cancer Act (1). Five decades later, despite ground-breaking discoveries and advancements in the field of cancer systems biology, a definitive and affordable cure for all types of cancer still evades humankind (2). Numerous “breakthrough” treatments have also gone on to exhibit adverse side effects (3, 4) that lower patients’ quality of life (QoL) or have reported degrading efficacies (5). At the heart of this problem lies our limited understanding of the bewildering multifactorial biomolecular complexity as well as patient-centricity of cancer.
Recent advances in biomolecular cancer research have helped factor system-level oncological manifestations into mutations across genetic, transcriptomic, proteomic, and metabolomic scales (6–9) that also act in concert (10, 11). Crosstalk between multi-scale pathways comprising of these oncogenic mutations can further exacerbate the etiology of the disease (7, 12–14). The combination of mutational diversity and interplay between the constituent pathways adds genetic heterogeneity and phenotypic plasticity in cancer cells (15, 16). Hanahan and Weinberg (17, 18) summarized this heterogeneity and plasticity into “Hallmarks of Cancer” – a set of progressively acquired traits during the development of cancer.
Experimental techniques such as high-throughput next-generation sequencing, and mass spectrometry-based proteomics are now providing specific spatiotemporal cues on patient-specific biomolecular aberrations involved in cancer development and growth. The voluminous high-throughput patient data coupled with the remarkable complexity of the disease has given impetus to data integrative in silico cancer modeling and therapeutic evaluation approaches (19). Specifically, scale-specific molecular insights into key regulators underpinning each hallmark of cancer are now helping unravel the complex dynamics of the disease (20) besides creating avenues for personalized therapeutics (21, 22). In this review, we will evaluate the emergence, evolution, and integration of multiscale cancer data towards building coherent and biologically plausible in silico models and their integrative analysis for employment in personalized cancer treatment in clinical settings. The review concludes by highlighting the need of integrating and modeling multi-omics data and associated software pipelines for employment in developing personalized therapeutics.
2 Scale-Specific Biomolecular Data and Its Applications in Cancer
Rapid advancements in molecular biology research, particularly in high-throughput genomics (23), transcriptomics (24), and proteomics (25) have resulted in the generation of big data on spatiotemporal measurements of scale-specific biomolecules (Figure 1A) in physiological as well as pathological contexts (26, 27). This vast and complex spatiotemporal data is expected to exceed 40 exabytes by 2025 (28) and is currently populating several online databases and repositories. These databases can be broadly categorized into seven salient database sub-types: biomolecules (29–31), pathways (32–34), networks (35–37), cellular environment (38, 39), cell lines (40–43), histopathological images (44–46), and mutations, and drug (47–52) databases, which are discussed below.
Figure 1 Overview of complex biomolecular regulation in cancer and scale-specific databases. (A) The complexity between genomic, transcriptomic, proteomic, metabolomic, cell-level, and environmental levels in a cancerous cell. Four examples of biomolecular signaling pathways are listed e.g., Hedgehog, Notch, (Wingless) Wnt, TGFβ, and AKT pathway. Stimuli from the extracellular environment signal the downstream pathway activation, in the cell, towards alternating the regulations in the proteomic, metabolomic, transcriptomic, and genomic scales, bringing out a system-level outcome in cancers. Lists (B) biomolecule (genes, transcripts, proteins, and metabolites) databases such as GenBank, GEO, TCGA, HPP, HMDB, etc. (C) Pathways databases such as PathDB, KEGG, STRING, etc. (D) Networks databases such as BioGRID, DIP, BIND, etc. (E) Environment databases e.g, ExoCarta, MatrixDB, MatrisomeDB, etc. (F) Cell lines databases such as CCLE, CLDB, Cellosaurus, etc. (G) Histopathological image database, for instance, TCIA, GTEx, TMAD, etc., and (H) Mutation and drug databases such as DrugBank, KEGG Drug, OncoKB, etc.
2.1 Biomolecular Databases
Biomolecular and clinical data generated from large-scale omics approaches for cancer research can be divided into four sub-categories: (i) genome, (ii) transcriptome, (iii) proteome, and (iv) metabolome (53).
2.1.1 Genome-Scale Databases
The foremost endeavor to collect and organize large-scale genomics data into coherent and accessible repositories led to the establishment of GenBank in 1986 (54) (Figure 1B, Table S1). This open-access resource now forms one of the largest public databases for nucleotide sequences from large-scale sequencing projects comprising over 300,000 species (55, 56). In a salient study employing GenBank, Diez et al. (57) screened breast and ovarian cancer families with mutations in BRCA1 and BRCA2 genes and its distribution in the Spanish population. Medrek et al. (58) employed microarray profile sets from GenBank to analyze gene levels for CD163 and CD68 in different breast cancer patient groups. The study established the need for localization of tumor-associated macrophages as a prognostic marker for breast cancer patients. To date, GenBank remains a comprehensive nucleotide database; however, its data heterogeneity poses a significant challenge in its employment in the development of personalized cancer therapeutics. Towards an improved data stratification and retrieval of genome-scale data, in 2002, Hubbard et al. (59) launched the Ensembl genome database. Ensembl provides a comprehensive resource for human genome sequences capable of automatic annotation and organization of large-scale sequencing data. Amongst various genome-wide studies utilizing Ensembl, Easton et al. (60) used this database to extract human sequence information to identify novel breast cancer susceptibility loci. Patient-specificity (10, 11) and mutational diversity (19) in cancer can manifest across spatiotemporal scales. Hence, the availability of patient-specific data for each type of cancer can furnish valuable insights into the biomolecular foundation of the disease. In an attempt to provide cancer type-specific mutation data, Wellcome Trust’s Sanger Institute developed Catalogue of Somatic Mutations in Cancer (COSMIC) (61) database. COSMIC comprises of 10,000 somatic mutations from 66,634 clinical samples. Schubbert et al. (62) employed COSMIC’s mutation data to investigate Ras activity in cancers as well as developmental disorders. The study concluded that the duration, as well as the strength of hyperactive Ras signaling controls the probability of tumorigenesis. Similarly, Weir et al. (63) utilized COSMIC data on tumor-suppressor and proto-oncogenes in the study to characterize the genome of lung adenocarcinomas. The systematic copy-number analysis with SNP data indicated that several lung cancer genes remain to be elucidated and characterized.
2.1.2 Transcriptome-Scale Databases
Gene-level information can facilitate the development of personalized cancer models; however, gene expression may vary from cell to cell and across cancer patients. As a result, cancer patients have divergent genetic signatures and transcript-level information. Hence, high-throughput transcriptomic data has the potential to provide valuable insights into the transcriptomic complexity in cancer cells and can be useful in investigating cell state, physiology, and relevant biological events (64) (Figure 1B and Table S1). Towards developing such a transcriptional information resource, in 2000, Edgar et al. (31) launched the Gene Expression Omnibus (GEO) initiative. GEO acts as a tertiary resource providing coherent high-throughput transcriptomic and functional genomics data. The platform now hosts over 3800 datasets and is expanding exponentially. GEO was employed by Chakraborty et al. (65) for annotation of chemo-resistant cell line models which helped investigate chemoresistance and glycolysis in ovarian cancers. The study identified Mitochondrial Calcium Uptake 1 (MICU1) as an important component for cancer metabolism that influences aerobic glycolysis and chemoresistance and can have a potential role in cancer therapeutics. The curation of patient-specific gene and protein expression data led to the development of The Cancer Genome Atlas (TCGA) (66). TCGA also captures the copy number variations and DNA methylation profiles for different cancer subtypes. TCGA’s potential (49, 67) is well exhibited by Leiserson et al.’s (68) pan-cancer analysis which helped identify 14 significantly mutated subnetworks containing numerous genes with rare somatic mutations across many cancers types. Davis et al. (69) further evaluated the genomic landscape of chromophobe renal cell carcinomas (ChRCCs) to elicit molecular patterns as clues for determining the origin of cancer cells. To facilitate in data management across different cancer projects as well as to ensure data uniformity towards developing data-driven models, the International Cancer Genome Consortium (ICGC) was launched in 2010 (70). ICGC adopts a federated data storage architecture that enables it to host a collection of scale-specific data from TCGA and 24 other projects (71). Burn et al. (72) estimated the distribution of cytosine in liver tumor data using ICGC. The study reported APOBEC3B (A3B) catalyzed deamination as a chronic source of DNA damage in breast cancer which also explains tumor cell evolution and heterogeneity. Supek et al. (73) compared mutation rates between different human cancers and reported the influence of “silent” mutations as a frequent contributor to cancer. Numerous databases have been established to store large-scale genomic data, however, insights from an integrated analysis of genomic data across databases have the potential to provide precise biomolecular cues into complex processes and evolution in cancer cells. This was enabled by cBio Cancer Genomics Portal (cBioPortal) (74) in 2012, with multidimensional dataset retrieval, and exploration from multiple databases. The platform additionally provides data visualization tools, pathway exploration, statistical analysis, and selective data download features for seamless utilization of large-scale genomics data across genes, patient samples, projects, and databases (75). Numerous studies have effectively employed cBioPortal (76–78); in particular, Jiao et al. (79) evaluated the prognostic value of TP53 and its correlation with EGFR mutations in advanced non-small-cell lung cancer (NSCLC). The study established that TP53 coupled with EGFR mutation can lead to the more accurate prognosis of advanced NSCLC. Hou et al. (80) also used cBioPortal to deduce targetable genotypes which are present in young patients with lung adenocarcinomas and revealed that young patients with lung adenocarcinoma were more likely to harbor targetable genotype.
2.1.3 Proteome-Scale Databases
Transcriptomic data remains limited in providing a deterministic proteome profile (81–83). Particularly, transcripts produced in a cell can be degraded, translated inefficiently, or modified due to post-translational modification (84, 85) resulting in no or a very small amount of functional protein (64). This relatively low correlation between transcriptome and proteome data was highlighted in 2019, by Bathke et al. (86) where it was shown that an increase in transcript synthesis cannot be directly associated with an increase in functional response in a cell. To facilitate functional analysis, there is a need to utilize proteomic-level data, which can help to capture a more accurate quantitative assessment of complex biomolecular regulation for functional studies (Figure 1B and Table S1). Following the successful completion of the Human Genome Project (HGP) (1998), in 2003, a group of Swedish researchers reported the Human Protein Atlas (HPA) (87, 88) with an aim to map the entire set of human proteins in cells, tissues, and organs for normal as well as cancerous state (89). HPA employs large-scale omics-based technologies to localize and quantify protein expression patterns. The database has successfully managed to host comprehensive information on human proteins from cells, tissues, pathology, brain, and blood region-related studies. HPA data can be employed for various purposes such as investigating the spatial distribution of proteins in different tissue and comparing normal and cancerous protein expression patterns across samples, etc (87). In a salient study employing HPA, Gámez-Pozo et al. (90) studied the localized expression pattern of proteins to help profile human lung cancer subtypes. The study reported a combination of peptide-level expressions which can distinguish between non-small lung cancer samples and normal lung cancer in different histological subtypes. Imberg-Kazdan et al. (91) employed HPA to identify novel regulators of androgen receptor (AR) function in prostate cancer towards therapy. Another significant stride towards generation of proteome-level information came with the establishment of the Human Proteome Project (HPP) (92, 93) in 2008 (94) by the Human Proteome Organization (HUPO). HPP consolidated mass spectrometry-based proteomics data, and bioinformatics pipelines, with the aim to organize and map the entire human proteome. To date, numerous studies have utilized HPP towards identifying the complex protein machinery involved in cancer cell fate outcomes (95–100). Amongst the earliest attempts, in 2001, Sebastian et al. (101) employed the HPP platform to deduce the complex regulatory region of the human CYP19 gene (‘armatose’), one of the contributors of breast cancer regulation. HPP project was later segmented into “biology and disease-oriented HPP” (B/D HPP) (102) and chromosome-centric HPP (C-HPP) (103). Specifically, Gupta et al. (104) carried out an extensive analysis of existing experimental and bioinformatics databases to annotate and decipher proteins associated with glioma on chromosome 12, while, Wang et al. (95) performed a qualitative and quantitative assessment of human chromosome 20 genes in cancer tissue and cells using C-HPP resources. The study revealed that several cancer-associated proteins on chromosome 20 were tissue or cell-type specific.
2.1.4 Metabolome-Scale Databases
Metabolic reprogramming is one of the earliest manifestations during tumorigenesis (105) and therefore, potentiates the usefulness of identifying metabolic biomarkers involved in cancer onset, its prognosis, as well as treatment. Large-scale efforts to collect metabolomics data have led to the development of several online databases (106–108) (Figure 1B and Table S1) including the Golm Metabolome Database (GMD) (106), in 2004. GMD provides a comprehensive resource on metabolic profiles, customized mass spectral libraries, along spectral information for use in metabolite identification. GMD was employed in 2011 by Wedge et al. (109) to identify and compare metabolic profiles in serum and plasma samples for small-cell lung cancer patients towards determining optimal agent for onwards analysis. The study showed that the discriminatory ability of both serum and plasma was equivalent. In 2013, Pasikanti et al. (110) utilized GMD to identify biomarker metabolites present in bladder cancer. The study proposed a potential role of kynurenine in malignancy and therapeutic intervention in bladder cancer. To allow for large-scale metabolic data stratification and retrieval, in 2007, Wishart et al. published the Human Metabolome Database (HMDB) (107, 111). HMDB contains organism-specific information on metabolites across various biospecimens and their accompanying environments. It is now the world’s largest metabolomics database with around 114,100 metabolites that have been characterized and annotated. HMDB was employed by Sugimoto et al. (112) to identify environmental compounds specific to oral, breast, and pancreatic cancer profiles. The study identified 57 principal metabolites to help predict disease susceptibility, besides being potential markers for medical screening. Agren et al. (113) employed metabolomic data from HMDB to construct metabolic network models for 69 human cells and 16 cancer types. The study’s comprehensive metabolic network analysis between disease and normal cell types has the potential to provide avenues for the identification of cancer-specific metabolic targets for therapeutic interventions. The HMDB supports data deposition and dissemination, however, integrated exploratory analysis is not available. The Metabolomics Workbench (108), reported in 2016, provides information on metabolomics metadata and experimental data across species, along with an integrated set of exploratory analysis tools. The platform also acts as a resource to integrate, deposit, track, analyze, as well as disseminate large-scale heterogeneous metabolomics data from a variety of studies. In a case study built using this platform, Hattori et al. (114) studied the aberrant BCAA (branched-chain amino acids) metabolism activation by MSI2 (Musashi2)-BCAT1 axis which they reported to drive myeloid leukemia progression.
2.2 Biomolecular Pathway Databases
Investigations restricted to single biomolecular scales have limited translational potential as cancer dysregulation is driven by tightly coupled biomolecular pathways constituted by biomolecules from a variety of spatiotemporal scales (discussed above). Such biomolecular pathways represent organized cascades of interactions integrating different spatiotemporal scales towards reaching specific phenotypic cell fate outcomes. Numerous scale-specific and multi-omics biomolecular pathway databases now exist to help retrieve, store and analyze existing pathway information towards understanding cellular communication in light of complex cancer regulation (32, 34, 115). One of the earliest attempts at integrating genomics data for pathway construction came in 1995, with the establishment of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (32) (Figure 1C and Table S2). Over time, KEGG has significantly expanded to include high-throughput multi-omics data (116). As a result, the resource is divided into fifteen sub-groups including KEGG Genome (for genome-level pathways), KEGG Compound (for small molecules level pathways), KEGG Gene (for gene and protein pathways), KEGG Reaction (for biochemical reaction and metabolic pathways), etc. Li et al. (117) used the KEGG database to perform pathway enrichment analysis for predicting the function of circular RNA (circRNA) dysregulation in colorectal cancer (CRC). The study highlighted circDDX17 potential role as a tumor suppressor and biomarker for CRC. While employing KEGG, Feng et al. (118) identified four up-regulated differentially expressed genes associated with poor prognosis in ovarian cancer. Further, to furnish information on pathways for high-throughput functional analysis studies, PANTHER (protein annotation through evolutionary relationship) database was established in 2010 (34, 115). PANTHER hosts information on ontological gene and protein-protein interaction pathways by leveraging GenBank and Human Gene Mutation Database (HGMD) (119), etc. Turcan et al. (120) employed PANTHER to perform network pathway enrichment for biological processes in differentially expressed genes, especially to investigate IDH1 mutations in glioma hypermethylation phenotype. The study provided a framework for understanding gliomas and the interplay between genomic as well as epigenomic regulation in cancer. To store metabolic pathway information in a cell, Karp et al. (121) developed MetaCyc, a comprehensive reference database comprising of metabolic pathways. MetaCyc is currently available as a web-based resource with metabolic pathway information which can be employed to investigate metabolic reengineering in cancers, carry out biochemistry-based studies, and explore cancer cell metabolism, etc. Miller et al. (122) demonstrated the utility of MetaCyc database by evaluating plasma metabolomic profiles after limonene intervention in breast cancer patients. The study employed MetaCyc to perform pathway-based interpretations and revealed that such alterations were related with tissue-level cyclin D1 expression changes.
2.3 Biomolecular Network Databases
The regulatory complexity of cancer is compounded by the crosstalk between numerous multi-scale biomolecular pathways resulting in the formation of complex interaction networks (Figure 1D and Table S3). One of the earliest biomolecular network databases, the Biomolecular Interaction Network Database (BIND) (123) was established in 2001, with an aim to organize biomolecular interactions between genes, transcripts, proteins, metabolites, as well as small molecules. Chen et al. (124) employed BIND to construct a biological interaction network (BIN) towards investigating tyrosine kinase regulation in breast cancer development. The study identified SLC4A7 and TOLLIP as novel tyrosine kinase substrates which are also linked to tumorigenesis. BIND provides a comprehensive resource of predefined interacting pathways; however, it does not contain ‘indirect’ interaction information. In contrast, the Molecular INTeraction database (MINT) (37), developed in 2002, curates existing literature to develop networks with both direct as well as indirect interactions from large-scale projects with information from genes, transcripts, proteins, promoter regions, etc. MINT can store data on “functional” interaction such as enzymatic properties and modifications present in biomolecular regulatory networks. The database was employed by Vinayagam et al. (125) to construct a human immunodeficiency virus (HIV) network that helped identify novel cancer genes across genomic datasets. The Database of Interacting Proteins (DIP) (126) was developed to mine existing literature and experimental studies on biomolecules and their pathways to construct protein interaction networks (127–129). Goh et al. (127) constructed a protein-protein interaction network for investigating liver cancer using DIP. The study revealed that hepatocellular carcinoma (HCC) at moderate stage is enriched in proteins that are involved in the immune response. Similarly, Zhao et al. (128) identified autophagic pathways in plants lectin-treated cancer cells which are regulated by microRNAs. The study showed that plant lectin has the potential to block sugar-containing receptor EGFR-mediated pathways thereby leading to autophagic cell death. To further consolidate and integrate protein interaction data across pathways as well as organisms, the Search Tool for the Retrieval of Interacting Genes/Proteins – STRING database was developed in 2005 (130). STRING provides a comprehensive text-mining and computational prediction platform which is accessible through an intuitive web interface (131). STRING database also provides information on interaction weights for edges between biomolecules to show an estimated likelihood for each interaction in the network (131). Mlecnik et al. (132) employed STRING database to study T-cells homing factors in colorectal cancer and demonstrated that specific chemokines and adhesion molecules had high densities of T-cell subsets in cancer.
2.4 Cellular Environment Databases
Each pathway within a biomolecular network requires input cues from the extracellular environment for onward downstream signal transduction (133–135). In the case of cancer, the biomolecular milieu constituting the tumor microenvironment (TME) acts as a niche for tumor development, metastasis, as well as therapy response (136) (Figure 1E and Table S4). Efforts to curate information from the environmental factors such as metabolites, matrisome, and other environmental compounds led to the development of MatrixDB (38) in 2011, which hosts matrix-based information on interactions between extracellular proteins and polysaccharides. MatrixDB additionally links databases with information on genes encoding extracellular proteins such as Human Protein Atlas (137) and UniGene (138) as well as host transcripts information. Celik et al. (139) employed MatrixDB data to evaluate epithelial-mesenchymal transition (EMT) inducers in the environment using nine ovarian cancer datasets and discovered a novel biomarker, HOPX, with the potential to drive tumor-associated stroma. To host studies on extracellular matrix (ECM) proteins from normal as well as disease-inflicted tissue samples, MatrisomeDB (39) was established in 2020. The database curates 17 studies on 15 physiologically healthy murine and human tissue as well as 6 cancer types from different stages (including breast, colon, ovarian, and lung cancer) along with other diseases. Levi-Galibov et al. (140) employed MatrisomeDB to investigate the progression of chronic intestinal inflammation in colon cancer. The study revealed the critical role of heat shock factor 1 (HSF1) during early changes in extracellular matrix structure as well as its composition.
2.5 Cell Line Databases
In vitro cell lines derived from cancer patients have become an essential tool for clinical and translational research (141). These cell lines are defined based on gene expression profiles and morphological features which have been cataloged in various databases such as the Cancer Cell Line Encyclopedia (CCLE) (42) (Figure 1F and Table S5). CCLE contains mutation data on 947 different human cancer cell lines coupled with pharmacological profiles of 24 anti-cancer drugs (42) for evaluating therapeutic effectiveness and sensitivity. Li et al. (142) employed CCLE data to investigate cancer cell line metabolism. The study showed that the mutated asparagine synthetase (ASNS) hypermethylation can cause gastric as well as hepatic cancers to sensitized asparaginase therapy. Hanniford et al. (143) demonstrated epigenetic silencing of RNA during invasion and metastasis in melanoma using the CCLE database. Other cell line databases include Cell Line Data Base (CLDB) (43), and The COSMIC Cell Lines Project (40), and CellMinerCDB (41). The CellMinerCDB (2018) curates data from National Cancer Institute (NCI) (144), BROAD institute (145), Sanger institute (146), and Massachusetts General Hospital (MGH) (147) and provides a platform for pharmacological and genomic analysis.
2.6 Histopathological Image Databases
Additionally, histopathological image datasets derived from the microscopic examination of tumor biopsy samples furnish information on cellular structure, function, chemistry, morphology, etc. Numerous histopathological image-based databases have been developed to store, manage, and retrieve such information (Figure 1G and Table S6). Amongst these databases, The Cancer Imaging Archive (TCIA) (46), reported in 2013, provides a multi-component architecture with various types of images including region-specific (e.g., Breast), cancer-type specific (e.g., TCGA-GBM, TCGA-BRCA), radiology, and anatomy images (e.g., Prostate-MRI). The cancer image collection in TCIA has been captured using a variety of modalities including radiation treatment, X-ray, mammography, and computed tomography (CT), etc (148). Li et al. (149) exploited TCIA by using radiomics data in predicting the risk for breast cancer recurrence, while Sun et al. (150) employed image data to perform a cohort study to validate a radiomics-based biomarker in cancer patients. The study developed a radiomic signature for CD8 cells using the TCGA dataset for predicting the immune phenotype of tumors and deduce clinical outcomes. In 2013, image data from TCIA was integrated with The Cancer Digital Slide Archive (CDSA) (44). The CDSA hosts imaging as well as histopathological data and provides more than 20,000 whole-slide images of 22 different cancer types. The whole-slide images of individual patients present in CDSA help in linking tumor morphology with the patient’s genomic and clinical data. Khosravi et al. (151) performed a deep convolution study using CDSA, to distinguish heterogeneous digital pathology images across different types of cancers. To associate patient’s genetic information and histology images, Genotype-Tissue Expression (GTEx) (45) was reported in 2014, and was curated using datasets from non-disease tissues of 1000 individuals. Patel et al. (152) employed GTEx to investigate intratumoral heterogeneity present in glioblastoma and concluded that glioblastoma subtype classifiers are variably expressed in individual cells.
2.7 Mutation and Drug Databases
Pharmacological investigations have elucidated the mechanism as well as efficacies of numerous cancer drugs, in clinical and preclinical studies (153, 154). Databases with drug-target information can be employed in precision oncology towards designing efficacious patient-centric therapeutic strategies (Figure 1H, Table S7). These databases include DrugBank (155), which was established in 2006 and contains information from over 4100 drug entries, 800 FDA-approved small molecules, and 14,000 protein or drug target sequences. DrugBank combines drug data with drug-target information thus enabling applications in cancer biology including in silico drug target discovery, drug design, drug interaction prediction, etc. In a study employing DrugBank, Augustyn et al. (156) evaluated potential therapeutic targets of achaete-scute homolog 1 (ASCL1) genes in lung cancers and reported unique molecular vulnerabilities for potential therapeutics, while Han et al. (157) determined synergistic combinations of drug targets in K562 chronic myeloid leukemia (CML) cells including BCL2L1 and MCL1 combination. Further to evaluate drugs in light of the patient’s genomic signature, PanDrugs (52) database was established in 2019 and currently hosts data from 24 sources and 56297 drug-target pairs along with 9092 unique compounds and 4804 genes. Using PanDrugs, Fernández-Navarro et al. (158) prioritized personalized drug treatments using PanDrugs, for T-cell acute lymphoblastic (T-ALL) patients.
Altogether, the availability of voluminous high-resolution biomolecular data has enabled the development of a quantitative understanding of aberrant mechanisms underpinning hallmarks of cancer as well as create avenues for personalized therapeutic insights. Recently, Karimi et al. (159) systematically evaluated the current multi-omics data generation approaches, as well as their associated analysis pipelines for employment in cancer research.
3 Data-Driven Integrative Modeling in Cancer Systems Biology
The need to prognosticate system-level outcomes in light of oncogenic dysregulation (160, 161) has led researchers to develop integrative data-driven computational models (162–168). Such models can help decode emergent mechanisms underpinning tumorigenesis as well as aid in the development of patient-centered therapeutic strategies (162–168). These in silico models can be broadly grouped into four salient sub-scales as biomolecular (169–171), tumor environment (172–174), cell level (175, 176), and multi-scale integrative cancer models (177–179) (Figure 2).
Figure 2 Evolution timeline of in silico scale-specific and multi-scale data-drive cancer models. Timeline of salient in silico scale-specific and multi-scale cancer models, along with PubMed yearly report (1990-2020) to display the evolutionary trends seen in the development of (A) genome-scale cancer models, (B) Transcript-level cancer models, (C) Proteome-scale models, (D) Metabolome scale models, (E) Environment-based models, (F) Cell-level models, and (G) Multi-scale cancer models.
3.1 Biomolecular-Scale Models
In silico biomolecular models of cancer can be classified into (i) genome-scale, (ii) transcriptome-scale, (iii) proteome-scale, (iv) metabolome-scale models.
3.1.1 Genome-Scale Models
Amongst initial attempts at developing cancer gene regulation models, in 1987, Leppert et al. (169) reported a computational model to study the genetic locus of familial polyposis coli and its involvement in colonic polyposis and colorectal cancer. The study was also validated by Mehl et al. (170) in 1991, which further elucidated the formation and development of familial polyposis coli genes in colorectal cancer patients. In 2014, Stratmann et al. (171) developed a personalized genome-scale 3D lung cancer model to study epithelial-mesenchymal transition (EMT) by TGFβ-based stimulation, while in 2016, Margolin et al. (180) developed a blood-based diagnostic model to help detect DNA hypermethylation of potential pan-cancer marker ZNF154. In 2017, Jahangiri et al. (181) employed an in silico pipeline to evaluate Staphylococcal Enterotoxin B for DNA-based vaccine for cancer therapy. A similar genome-scale model of DNA damage and repair was proposed by Smith et al. (182) to evaluate proton treatment in cancer (Figure 2A).
3.1.2 Transcriptome-Scale Models
In silico models developed using transcriptomic expression data can assist in comparing gene expression patterns in cancer for investigating genetic heterogeneity and cancer development as well as towards precision therapy (Figure 2B). In 2003, Huang et al. (183) presents a mathematical model using a large-scale transcriptional dataset of breast cancer patients to elucidate patterns of metagenes for nodal metastases and relapses. In another large-scale cancer transcriptomics study, 3,000 patient samples from nine different cancer types were used to decode the genomic evolution of cancer by Cheng et al. (184). In 2013, Chen et al. (185) employed an in silico pipeline which helped identify 183 new tumor-associated gene candidates with the potential to be involved in the development of hepatocellular carcinoma (HCC), while, in 2014, Agren et al. (186) developed a personalized transcriptomic data-based model to identify anticancer drugs for HCC. In 2019, Béal et al. (187) reported a logical network modeling pipeline for personalized cancer medicine using individual breast cancer patients’ data. The pipeline was validated in 2021 (188), using in silico personalized logical models for melanomas and colorectal cancers samples in response to BRAF treatments. In a similar study conducted in 2019, Rodriguez et al. (189) developed a mathematical model for breast cancer using transcriptional regulation data to predict hypervariability in a large dynamic dataset which revealed the basis of expression heterogeneity in breast cancer.
3.1.3 Proteome-Scale Models
To capture the quantitative aspects of biomolecular regulations and functional studies (82, 83), data-driven proteomic-based cancer models are essential (Figure 2C). Such models can be particularly helpful in diagnostic as well as prognostic purposes as well as for monitoring response to treatment (82, 83). In a study employing proteome-level information, in 2011, Baloria et al. (190) carried out an in silico proteome-based characterization of the human epidermal growth factor receptor 2 (HER-2) to evaluate its immunogenicity in an in silico DNA vaccine. Akhoon et al. (191) simplified this approach with the development of a new prophylactic in silico DNA vaccine using IL-12 as an adjuvant. In 2017, Fang et al. (192) employed proteome level data towards predicting in silico drug-target interactions for applications in targeted cancer therapeutics, while in 2018, Azevedo et al. (193) designed novel glycobiomarkers in bladder cancer. Recently, in 2020, Lee et al. (194) reported an integrated proteome model of macrophage migration in a complex tumor microenvironment. However, proteome-level models are limited in their ability to provide a complete analysis of the biomolecules present in a cell since they lack information on low molecular weight biomolecular compounds such as metabolites (105).
3.1.4 Metabolome-Scale Models
Metabolic data-driven models can be especially useful in understanding cancer cell metabolism, mitochondrial dysfunction, metabolic pathway alteration, etc (Figure 2D). In 2007, Ma et al. (195) developed the Edinburgh Human Metabolic Network (EHMN) model with more than 3000 metabolic reactions alongside 2000 metabolic genes for employment in metabolite-related studies and functional analysis. In 2011, Folger et al.’s (196) employed the EHMN model to propose a large-scale flux balance analysis (FBA) model for investigating metabolic alterations in different cancer types and for predicting potential drug targets. In 2014, Aurich et al. (197) reported a workflow to characterize cellular metabolic traits using extracellular metabolic data from lymphoblastic leukemia cell lines (Molt-4) towards investigating cancer cell metabolism. Yurkovich et al. (198) augmented this workflow in 2017 and reported eight biomarkers for accurately predicting quantitative metabolite concentration in human red blood cells. Alakwaa et al. (199) employed a mathematical model to predict the status of Estrogen Receptor in breast cancer metabolomics dataset, while in 2018, Azadi et al. (200) used an integrative in silico pipeline to evaluate the anti-cancerous effects of Syzygium aromaticum employing data from the Human Metabolome Database (107).
3.2 Tumor Environment Models
Integrative mathematical models of environmental cues and extracellular matrix can help researchers abstract tumor microenvironments. Such data-driven models can be used to study angiogenesis (201), cell adhesion (202), and vasculature (203), etc (Figure 2E). Amongst initial attempts at developing cancer environment-based in silico models, in 1972, Greenpan et al. (172) designed a solid carcinoma in silico model to evaluate cancer cell behavior in limited diffusion settings. In 1976, the model was expanded (173) to investigate tumor growth in asymmetric conditions. In 1996, Chaplain developed a mathematical model to elucidate avascular growth, angiogenesis, and vascular growth in solid tumors (174). Anderson and Chaplain (204), in 1998, expanded this strategy and reported a continuous and discrete model for tumor-induced angiogenesis. This modeling approach was further augmented in 2005 by Anderson (202) to a hybrid mathematical model of a solid tumor to study cellular adhesion in tumor cell invasion. Organ-specific metastases and associated survival ratios in small cell lung cancer patients have been modeled and evaluated (205) using similar models.
3.3 Cell Models
To model cell population-level behavior in cancer, researchers are increasingly developing innovative cell lines in silico models which can complement in vivo wet-lab experiments, while overcoming wet-lab limitations (206) (Figure 2F). Such models are employed to investigate cell-to-cell interactions as well as evaluate physical features of the synthetic extracellular matrix (ECM) (206), etc (175, 207). Amongst initial attempts at developing in silico cell line models, in 1989, Shackney et al. (175) proposed an in silico cancer cell model to study tumor evolution. Results showed an association between discrete aneuploidy peaks with the activation of growth-promoting genes. In 2007, Aubert et al. (207) developed an in silico glioma cell migration model and validated cell migration preferences for homotype and heterotypic gap junctions with experimental results. Gerlee and Nelander (176) expanded this work, in 2012, to investigate the effect of phenotype switching in glioblastoma growth. In conclusion, cell-based cancer models help provide scale-specific insights into cancer, however, they remain limited in investigating the spatiotemporal tissue diversity and heterogeneity in cancer patients.
3.4 Multi-Scale Models
Recently, data-driven multi-scale models are becoming increasingly popular in cancer (208) (Figure 2G). One of the earliest attempts at developing multi-scale cancer models was in 1985 when Balding developed a mathematical model to demonstrate tumor-induced capillary growth (201). In 2000, Swanson et al. (209) proposed a quantitative model to investigate glioma cells. In a similar study, Zhang et al. (210) generated a 3D, multi-scale agent-based model of the brain to study the role (EGFR)-mediated activation of signaling protein phospholipase role in a cell’s decision to either proliferate or migrate. In 2010, Wang et al. (178) also took a multi-scale agent-based modeling approach to identify therapeutic targets in concurrent EGFR-TGFβ signaling pathway in non-small cell lung cancer (NSCLC). Later in 2011 (179), they employed the approach to identify critical molecular components in NSCLC. Similarly, Perfahl et al. (211) formulated a multi-scale vascular tumor growth model to investigate spatiotemporal regulations in cancer and response to therapy. In 2007, Anderson et al. (212) proposed a mathematical model for studying cancer growth, evolution, and invasion. This model was later built upon by Chaudhary et al. (165, 213) with a multi-scale modeling strategy to investigate tumorigenesis induced by mitochondrial incapacitation in cell death, in 2011. In 2017, Vavourakis et al. (214) developed a multi-scale model to investigate tumor angiogenesis and growth. In 2017, Norton et al. (215) used an agent-based computational model of triple-negative breast cancer to study the effects of migration in CCR5+ cancer cells, stem cell proliferation, and hypoxia on the system. They later (216) reported an agent-based and hybrid model to investigate tumor immune microenvironment, in 2019. In the same year, Karolak et al. (217) modeled in silico breast cancer organoid morphologies (218) to help elucidate efficacies amongst drug treatment based on the morphophenotypic classification. Similarly, Berrouet et al. (219) employed a multi-scale mathematical model to evaluate the effect of drug concentration on monolayers and spheroid cultures.
Summarily, data-integrative computational models have now assumed the forefront in decoding the complex biomolecular regulations involved in cancer and are increasingly been employed for the development of personalized preclinical models as well as therapeutics design (220).
4 Software Platforms for Modeling in Cancer Systems Biology
Over the past decade, in silico modeling of cancer has gained significant popularity in systems biology research (162–168). In particular, data-drive computational models are now acting as an enabling technology for precision medicine and personalized treatment of cancer. To date, several single and multi-scale software platforms have been reported to model biomolecules (171, 221–228), cellular environments (229), as well as cell-level (230), and multi-scale (231) information into coherent in silico cancer models (Figure 3).
Figure 3 Feature-by-feature comparison of networks, environments, cell lines, and multi-scale modeling software in chronological order.
4.1 Biomolecular Modeling Platforms
Software platforms aimed at modeling biomolecular entities, abstract information from published literature as well as high-throughput technologies, and model them using a Boolean modeling approach or a differential equations modeling strategy.
4.1.1 Boolean Modeling Software
Boolean network modeling technique was first introduced by Kauffman (232, 233), in 1969. This approach has been widely adopted as a tool to model gene, transcript, protein, and metabolite regulatory networks. Numerous mathematical and computational cancer models have been developed using this representation (171, 221–225). To facilitate the Boolean model development and analysis process, several platforms have been devised (234–236) (Table S8). The applicability of these platforms can be further categorized into qualitative or quantitative Boolean network modeling.
4.1.1.1 Qualitative Boolean Modeling
Biomolecular qualitative Boolean models are a widely employed approach in cancer systems biology research to cater for cases where there is insufficient quantitative information, and/or lack of mechanistic understanding. Thus far, numerous platforms have been reported to help researchers develop qualitative Boolean network models. Amongst them, FluxAnalyzer (237), reported by Klamt et al. in 2003, was developed to undertake metabolic pathway construction, flux optimization, topological feature detection, flux analysis, etc. To expand the scope of the platform and include cell signaling, gene as well as protein regulatory networks, in 2007, Klamt et al. expanded FluxAnalyzer and reported CellNetAnalyzer (226). Tian et al. (238) employed CellNetAnalyzer to develop a p53 network model for evaluating DNA damage in cancer, while Hetmanski et al. (239) designed a MAPK-driven feedback loop in Rho-A-driven cancer cell invasion. Although CellNetAnalyzer remains a widely used logical modeling software, its programmability and MATLAB dependency hinders its clinical employment for developing personalized cancer models. Towards addressing this challenge, in 2008, Albert et al. published BooleanNet (235), an open-source, freely available Boolean modeling software for large-scale simulations of dynamic biological systems. Saadatpoort et al. (240) employed BooleanNet’s general asynchronous (GA) method to deduce therapeutic targets for granular lymphocyte leukemia. Similarly, in 2008, Kachalo et al. presented NET-SYNTHESIS (227); a platform for undertaking network synthesis, inference, and simplification. Steinway et al. (241) employed both BooleanNet and NET-SYNTHESIS platforms to model epithelial-to-mesenchymal transition (EMT) in light of TFGβ cell signaling, in hepatocellular carcinoma patients towards elucidating potential therapeutic targets. BooleanNet was used to undertake model simulation and NET-SYNTHESIS for carrying out network interference and simplification. Another logical modeling platform, GIMsim (242), published by Naldi et al., in 2009, also employed asynchronous state transition graphs to perform qualitative logical modeling which is especially useful for networks with large state space. This platform was employed by Flobak et al. (243) to map cell fate decisions in gastric adenocarcinoma cell-line towards evaluating drug synergies for treatment purposes, while Remy et al. (244) studied mutually exclusive and co-occurring genetic alterations in bladder cancer. GIMsim also employed multi-valued logical functions, useful in simulating qualitative dynamical behavior in cancer research. However, the platform was unable to program automatic theoretical predictions, moreover, it only employed qualitative analysis approaches and could not be used to accurately map cell fates based on quantitative biomolecular expression data. Taken together, classical qualitative Boolean modeling approaches remain limited in developing predictive cancer models that could leverage quantitative biomolecular expression data generated from next-generation proteomics and related-sequencing projects.
4.1.1.2 Quantitative Boolean Modeling
Platforms aimed to integrate quantitative expression data from existing literature and databases towards carrying out network annotation and onwards analysis can be particularly useful in developing personalized cancer models. One such platform, the Markovian Boolean Stochastic Simulator (MaBoss), was established by Stoll et al. (245) in 2017, for stochastic and semi-quantitative Boolean network model development, mutations, and drug evaluation, sensitivity analysis based on experimental data, and eliciting model predictions. In 2019, Béal et al. (187) employed MaBoss to develop a logical model to evaluate breast cancer in light of individual patients’ genomic signature for personalized cancer medicine. This model was later expanded, in 2021, to investigate BRAF treatments in melanomas and colorectal cancer patients (188). Similarly, Kondratova et al. (246) used MaBoss to model an immune checkpoint network to evaluate the synergistic effects of combined checkpoint inhibitors in different types of cancers. In a similar attempt at developing quantitative Boolean networks, BoolNet (247) was developed in 2010 by Müssel and Kestler. BoolNet allows its users to reconstruct networks from time-series data, perform robustness and perturbation analysis and visualize the resultant cell fates attractor. BoolNet was employed by Steinway et al. (248) to construct a metabolic network model towards evaluating gut microbiome in normal and disease conditions, whereas, Cohen et al. (249) studied tumor cell invasion and migration. BoolNet, however, lacks a graphical user interface, and results from the analysis cannot be visualized interactively, which hindered its employment. To address this issue, Shah et al. (250), in 2018, developed an Attractor Landscape Analysis Toolbox for Cell Fate Discovery and Reprogramming (ATLANTIS). ATLANTIS has an intuitive graphical user interface and interactive result visualization feature, for ease in utilization. The platform can be employed to perform deterministic as well as probabilistic analysis and was validated through literature-based case studies on the yeast cell cycle (251), breast cancer (252), and colorectal cancer (253).
4.1.2 Differential Equations Modeling Software
Boolean models have proven to be a powerful tool in modeling complex biomolecular signaling networks, however, these models are unable to describe continuous concentration, and cannot be used to quantify the time-dependent behavior of biological systems, necessitating the need to switch to quantitative differential equations (254). As a result, numerous stand-alone and web-based tools have been developed to build continuous network models to help describe the temporal evolution of biomolecules towards elucidating more accurate cell fate outcomes from quantitative expression data (Table S8). Amongst initial attempts at developing such software, GEPASI (GEneral PAthway Simulator) was designed by Pedro Mendes et al. (228), in 1993. GEPASI is a stand-alone simulator that facilitates formulating mathematical models of biochemical reaction networks. GEPASI can also be used to perform parametric sensitivity analysis using an automatic pipeline that evaluates networks in light of exhaustive combinatorial input parameters. Ricci et al. (255) employed GEPASI to investigate the mechanism of action of anticancer drugs, while Marín-Hernández et al. (256) constructed kinetic models of glycolysis in cancer. In 2006, Hoops et al. reported a successor of GEPASI; COPASI (COmplex PAthway SImulator) (257), a user-friendly independent biochemical simulator that can handles larger networks for faster simulation results, through parallel computing. Orton et al. (258) employed COPASI to model cancerous mutations in EGFR/ERK pathway, while cellular senescence was evaluated by Pezze et al. (259) for targeted therapeutic interventions. Towards establishing a user-friendly software with an intuitive graphical user interface (GUI), another desktop application, CellDesigner was published by Funahashi et al. (260). CellDesigner application can be extended to include various simulation and analysis packages through integration with systems biology workbench (SBW) (261). In a case study using CellDesigner, Calzone et al. (262) developed a network of retinoblastoma protein (RB/RB1) and evaluated its influence in cell cycle, while Grieco et al. (263) investigated the impact of Mitogen-Activated Protein Kinase (MAPK) network on cancer cell fate outcomes.
4.2 Cell Environment Modeling Software
To facilitate the development of environmental models that can help investigate inter-, intra-, and extracellular interactions between cellular network models and their dynamic environment, several software and platforms have been reported (229, 264). These platforms employ discrete, continuous, and hybrid approaches to develop models of cellular microenvironments towards setting up specific biological contexts such as normoxia, hypoxia, Warburg effect, etc (Table S9). In 2014, Starruß et al. (229) published Morpheus, a platform for modeling complex tumor microenvironment. Morpheus leverages a cellular potts modeling approach to integrate and stimulate cell-based biomolecular systems for modeling intra- and extra cellular dynamics. In a case study using Morpheus, Felix et al. (265) evaluated pancreatic ductal adenocarcinoma’s adaptive and innate immune response levels, while Meyer et al. (266) investigated the dynamics of biliary fluid in the liver lobule. Morpheus is a widely employed modeling software, however, its diffusion solver is limited in its capacity to model large 3D domains. Towards modeling fast simulations for larger cellular systems, the Finite Volume Method for biological problems (BioFVM) software (264) was reported by Ghaffarizadeh et al., in 2016. BioFVM is an efficient transport solver for single as well as multi-cell biological problems such as excretion, decomposition, diffusion, and consumption of substrates, etc (267). In a case study using BioFVM, Ozik et al. (268) evaluated tumor-immune interactions, while Wang et al. (269) elucidated the impact of tumor-parenchyma on the progression of liver metastasis. BioFVM, however, relies on its users to have programming-based knowledge to develop their models, which limits its translational potential. Towards minimizing programming requirements, SALSA (ScAffoLd SimulAtor) was developed by Cortesi et al. (206) in 2020. SALSA is general-purpose software that employs a minimum programming requirement, a significant advantage over its predecessors. The platform has been useful in studying cellular diffusion in 3D cultures. This recent tool was validated in 2021, with a case study that evaluated and predicted therapeutic agents in 3D cell cultures (270).
4.3 Cell-Level Modeling Software
Towards modeling cancer cell-specific behaviors such as cellular adhesion, membrane transport, loss of cell polarity, etc several software have been reported which can help develop in silico cancer cell models (Table S10). The foremost endeavor to develop software for cell-level modeling and simulation, led to the establishment of E-Cell (230), in 1999. E-Cell can be employed to model biochemical regulations and genetic processes using biomolecular regulatory networks in cells. Edwards et al. (271) employed E-cell to predict the metabolic capabilities of Escherichia coli and validated the results using existing literature, while Orton et al. (272) modeled the receptor-tyrosine-kinase-activated MAPK pathway. This software was expanded in 2001, with the development of Virtual Cell (V-Cell) (273), a web-based general-purpose modeling platform. V-Cell has an intuitive graphical and mathematical interface that allows ease in the design and simulation of whole cells, along with sub-cellular biomolecular networks and the external environment. Neves et al. (274) employed V-Cell to investigate the flow of spatial information in cAMP/PKA/B-Raf/MAPK1,2 networks, while Calebiro et al. (275) modeled cell signaling by internalized G-protein–coupled receptors. Similarly, to increase the efficiency in the design and modeling of synthetic regulatory networks in cells, in 2009, Chandran et al. reported TinkerCell (276) with computer-aided design (CAD) functionality which enabled faster simulation and associated analysis. The platform employed a modular approach for constructing networks and provides built-in features for ease in network construction, robustness analysis, and evaluating networks using existing databases. The evolutionary trend for TinkerCell’s platform adaptability and flexibility is highlighted in Figure 4, which shows a gradual shift from being a model-specific platform to a domain modeling one. In a study employing TinkerCell, Renicke et al. (277) constructed a generic photosensitive degron (psd) model to investigate protein degradation and cellular function, while Chandran et al. (278) reported computer-aided biological circuits.
Figure 4 Evolution of scale-specific and multi-scale software. Evolution of multi-scale modeling software for abstracting and simulating the spatiotemporal biomolecular complexity. Highlighting the need for a generic, data-driven, zero-code software requirement.
4.4 Multi-Scale Modeling Software
Multi-scale cancer modeling approaches bring together scale-specific information towards undertaking an integrative analysis of heterogeneous experimental data by building coherent and biologically plausible models (279–281). Several multi-scale modeling platforms have been reported to help develop multi-level cancer models (231) (Table S11). Amongst these, the REcursive Porous Agent Simulation Toolkit (Repast) (282) published in 2003, provides a free and open-source tool for modeling and simulating agent-based models, with high-performance computing (HPC) capability. Repast toolkit was employed in 2007, by Folcik et al. (283) to develop an agent-based model which was then used to study interactions between cells and the immune system. Similarly, Mehdizadeh et al. (284) used Repast to model angiogenesis in porous biomaterial scaffolds. Although Repast can be used for simulating several types of evolutionary trends between agents, there is no established guideline for selecting a mechanism to model such trends, hindering its use by naïve users. Moreover, Repast does not have a GUI or a software development kit (SDK) interface for implementing subcellular mechanisms e.g., gene, protein, and metabolic networks. In contrast, CompuCell (279), published in 2004 by Izaguirre et al., provides an elaborative GUI to model cell-scale or tissue-scale simulations by integrating biomolecular networks, intra- and extracellular environment, and cell to environment interactions. Mahoney et al. (285) employed CompuCell to develop an angiogenesis-based model in cancer for investigating novel cancer therapies. Although CompuCell provides an intuitive framework modeling paradigm, the platform core is not conducive to multi-scale cancer modeling. The focus of the software is primarily multi-agent simulations rather than multi-scale cancer modeling. This poses a significant challenge in the utilization of the software. In 2013, CHASTE (280) (Cancer Heart and Soft Tissue Environment) was launched, which provides a computational simulation pipeline for the mathematical modeling of complex multi-scale models. Users can employ CHASTE for a wide range of problems involving on and off-lattice modeling workflows. CHASTE has also previously been employed to model colorectal cancer crypts. Nonetheless, CHASTE does not have a GUI and can only be executed by command line text commands. Furthermore, it requires recompilation on part of the modeler to use the code updates performed by the group. To further improve the multi-scale modeling approach for investigating cancer, in 2013, Chaudhary et al. (281) published ELECANS (Electronic Cancer System). ELECANS had an intuitive but rigid GUI along with a programmable SDK besides the lack of a high-performance simulation engine (286, 287). ELECANS was employed to model the mitochondrial processes in cancer towards elucidating the hidden mechanisms involved in cell death (165). ELECANS provided a feature-rich environment for constructing multi-scale models, however, the platform lacked a biomolecular network integration pipeline, and also placed a heavy programming requirement on its users. In contrast, in 2018, PhysiCell (288) was reported for 2D and 3D multi-cell off-lattice agent-based simulations. PhysiCell is coupled with BioFVM (264)’s finite volume method to model multi-scale cancer systems (268). In a salient example, Wang et al. (269) employed PhysiCell to model liver metastatic progression. In 2019, PhysiCell’s agent-based modeling features and MaBoss’s Boolean cell signaling network feature were coupled together to develop an integrated platform, PhysiBoss (289). As a result, PhysiBoss provided an agent-based modeling environment to study physical dimension and cell signaling networks in a cancer model. In 2020, Colin et al. (290) employed PhysiBoss’s source code to model diffusion in oocytes during prophase 1 and meiosis 1, while Getz et al. (291) proposed a framework using PhysiBoss to develop a multi-scale model of SARS-CoV-2 dynamics in lung tissue. Recently, in 2021, Gondal et al. (292) reported Theatre for in silico Systems Oncology (TISON), a web-based multi-scale “zero-code” modeling and simulation platform for in silico oncology. TISON aims to develop single or multi-scale models for designing personalized cancer therapeutics. To exemplify the use case for TISON, Gondal et al. employed TISON to model colorectal tumorigenesis in Drosophila melanogaster’s midgut towards evaluating efficacious combinatorial therapies for individual colorectal cancer patients (225).
Summarily, multi-scale modeling software has enabled the development of biologically plausible cancer models to varying degrees. These platforms, however, fall short of providing a generic and high-throughput environment that could be conveniently translated into clinical settings.
5 Pipelining Panomics Data towards In Silico ClinicalSystems Oncology
In silico multi-scale cancer models, annotated with patient-specific biomolecular and clinical data, can help decode complex mechanisms underpinning tumorigenesis and assist in clinical decision-making. Clinically driven in silico multi-scale cancer models simulate in vivo tumor growth and response to therapies across biocomplexity scales, within a clinical environment, towards evaluating efficacious treatment combinations. To facilitate the development of multi-scale cancer models, several large-scale program projects have been launched (Table S12) such as Advancing Clinico-Genomic Trials on Cancer (ACGT) (293), Clinically Oriented Translational Cancer Multilevel Modeling (ContraCancrum) (294), Personalized Medicine (p-medicine) (295), Transatlantic Tumor Model Repositories (TUMOR) (296), and Computational Horizons In Cancer (CHIC) (297), amongst others (Figure 5). Here, we review and evaluate five salient projects for multi-scale cancer modeling towards their clinical deployment.
Figure 5 Salient projects pipelining multi-scale panomics data into clinical settings – a timeline. Timeline highlighting salient project platforms for developing realistic and clinically-driven multi-scale cancer models, along with their associated leading case studies.
5.1 ACGT Project
The ACGT project (293), launched in 2007, proposed to develop Clinico-Genomic infrastructure for organizing clinical and genomic data towards investigating personalized therapeutics regimens for an individual cancer patient. The ACGT platform provides an open-source and open-access infrastructure designed to support the development of “oncosimulators” to help clinicians accurately compare results from different clinical trials and enhance their efficiency towards optimizing cancer treatment. The ACGT framework employs molecular and clinical data generated from different sources including whole genome sequencing, histopathological, imaging, molecular, and clinical data, etc. to develop simulators for mimicking clinical trials. Personalized panomics data employed to develop oncosimulators in ACGT is extracted from real patients which enables the oncosimulators to be clinically relevant for predictive purposes. Additionally, ACGT provides data retrieval, storage, integrative, anonymization, and analysis as well as results presentation capabilities. Using the platform, in 2004, a personalized, spatiotemporal oncosimulator model of breast cancer was developed to mimic a clinical trial based on protocols outlined in the Trial of Principle (TOP), towards evaluating the model response to chemotherapeutic treatment in neoadjuvant settings (298). Similarly, in 2009, Graf et al. (299) modeled nephroblastoma oncosimulator, a childhood cancer of the kidney, based on a clinical trial run by the International Society of Paediatric Oncology (SIOP) for simulating tumor response to therapeutic regimens in clinical trials. The results generated from the TOP and SIOP trials enabled the ACGT oncosimulators to adapt in light of real clinical conditions and the software to be validated against multi-scale patient data. The focus of ACGT oncosimulators, however, is limited to existing clinical trials for predicting efficacious treatment combinations.
5.2 The ContraCancrum Project
Towards establishing a platform for the development of composite multi-scale models for simulating malignant tumor models, in 2008, ContraCancrum project (294) was initiated. ContraCancrum aimed to construct a multi-scale computational framework for translating personalized cancer models into clinical settings towards simulating malignant tumor development and response to therapeutic regimens. For that, the Individualized MediciNe Simulation Environment (IMENSE) platform (300, 301) was established, to undertake the oncosimulator development process. The platform was employed for molecular and clinical data storage, retrieval, integration, and analysis. Oncosimulators, developed under the ContraCancrum project, employ patient data across biologically and clinically relevant scales including molecular, environmental, cell, and tissues level. These oncosimulators can be used to optimize personalized cancer therapeutics for assisting clinician decision-making process. To date, several oncosimulators have been reported, under the ContraCancrum project (27, 294, 299, 301–305). In particular, the initial validation of the ContraCancrum workflow was performed using two case studies; glioblastoma multiforme (GBM) (294, 305) and non-small cell lung cancer (NSCLC) (294, 304). In 2010, Folarin and Stamatakos (306) designed a glioblastoma oncosimulator using personalized molecular patient data to evaluate treatment response under the effect of a chemotherapeutic drug (temozolomide) (300). Similarly, Roniotis et al. (305) developed a multi-scale finite elements workflow to model glioblastoma growth, while Giatili et al. (307) outlined explicit boundary condition treatment in glioblastoma using an in silico tumor growth model. The results from the case studies were validated by comparing in silico prediction with pre- and post-operative imaging and clinical data (294, 302). Moreover, The ContraCancrum project hosts more than 100 lung patients’ tumor and blood samples (294). This data is employed to develop clinically validated in silico multilevel cancer models for NSCLC using patient-specific data (304). In 2010, using a biochemical oncosimulator, Wan et al. (304) investigated the binding affinities for AEE788 and Gefitinib tyrosine kinase inhibitor against mutated epidermal growth factor receptor (EGFR) for NSCLC treatment. Similarly, Wang et al. (308) developed a 2D agent-based NSCLC model to investigate proliferation markers and evaluated ERK as a suitable target for targeted therapy. Although ContraCancrum’s project has created numerous avenues for multilevel cancer model development, the platform is limited to only specific types of cancer for which data is internally available on the platform.
5.3 p-Medicine Project
Towards improving clinical deployment capabilities of oncosimulators, another pilot project: the personalized medicine (p-medicine) project was launched in 2011 (295). The main aim of p-medicine was to create biomedical tools facilitating translation of current medicine to P5 medicine (predictive, personalized, preventive, participatory, and psycho-cognitive) (309). For that, the p-medicine portal provides a web-based environment that hosts specifically-purpose tools for personalized panomics data integration, management, and model development. The portal has an intuitive graphical user interface (GUI) with an integrated workbench application for integrating information from clinical practices, histopathological imaging, treatment, and omics data, etc. Computational models developed under p-medicine workflow, are quantitatively adapted to clinical settings since they are derived using real multi-scale data. Several multi-scale cancer simulation models (oncosimulators) (295, 309) have been devised, using the p-medicine workflow. Amongst these, in 2012, Georgiadi et al. (310) developed a four-dimensional nephroblastoma treatment model and evaluated its employment in clinical decisions making. Towards evaluating personalized therapeutic combinations, in 2014, Blazewicz et al. (311) reported a p-medicine parallelized oncosimulator which evaluated nephroblastoma tumor response to therapy. The parallelization enhanced model usability and accuracy for eventual translation into clinical settings for supporting clinical decisions. In 2016, Argyri et al. (312) developed a breast cancer oncosimulator to evaluate vascular tumor growth in light of single-agent bevacizumab therapy (anti-angiogenic treatment), while in 2014, Stamatakos et al. (313) evaluated breast cancer treatment under an anthracycline drug for chemotherapy (epirubicin). In 2012, Ouzounoglou et al. (314) designed a personalized acute lymphoblastic leukemia (ALL) oncosimulator for evaluating the efficacy of prednisolone (a steroid medication). This study was further augmented, in 2015, by Kolokotroni et al. (315) to investigate the potential cytotoxic side effects of prednisolone. Later in 2017, Ouzounoglou et al. (316) expanded the ALL oncosimulator model to a hybrid oncosimulator for predicting pre-phase treatments for ALL patients. The validation of these models was undertaken using clinical trial data in pre and post-treatment (317). Although p-medicine presents a state-of-the-art in silico multi-scale cancer modeling environment, the project is limited in its application for determining individual biomarkers for potential novel therapy identification. Moreover, the project is limited to its niche set of tools and models that cannot be integrated with other similar model repositories and platforms.
5.4 TUMOR Project
In 2012, a transatlantic USA-EU partnership was initiated with the launch of Transatlantic Tumor Model Repositories (TUMOR) (296). TUMOR provides an integrated, interoperable transatlantic research environment for developing a clinically driven cancer model repository. The repository aimed to integrate computational cancer models developed by different research groups. TUMOR’s transatlantic aim was to couple models from ACGT and ContraCancrum projects with those at the Center for the Development of a Virtual Tumor (CViT) (318) as well as other relevant centers. TUMOR is predicted to serve as an international clinical translation, interoperable and validation platform for in silico oncology hosting multi-scale cancer models from different cancer model repositories and platforms such as E-Cell (230), CellML (319), FieldML (320), and BioModels (35, 321). To support model data interoperability between platforms, TumorML (Tumor model repositories Markup Language) (321) was developed towards facilitating inter-data operability in the TUMOR project. The TUMOR environment also offers a wide range of additional services supporting predictive oncology and individualized optimization of cancer treatment. For example, the platform allows remote access of predictive cancer models in hospitals and to clinicians for the development of quantitative cancer research and personalized cancer therapy model. TUMOR environment also incorporates deterministic as well as stochastic models through COPASI simulator (257). An automatic validation pipeline is also embedded for the execution and deployment of these models in clinical settings (296). TUMOR’s ability to couple and integrate models from different scales, approaches, and repositories towards increasing model accuracy in predictive oncology was exemplified with Wang et al.’s (322) NSCLC model. In 2007, Wang et al. developed a 2D multi-scale NSCLC model for evaluating tumor expansion dynamics in NSCLC patients, in CViT. This model was exported in TumorML and made available in the TUMOR repository to help evaluate growth factor influence in aggressive cancer (321). Similarly, in 2014 Sakkalis et al. (323) employed the TUMOR platform to interlink and couple three independent glioblastoma-specific cancer models (EGFR signaling (324), cancer metabolism (325), Oncosimulator (323), reported by different research groups. The resultant model was used to investigate the impact of radiation and temozolomide (chemotherapy) on glioblastoma multiforme to evaluate treatment effectiveness.
5.5 CHIC Project
Another transatlantic project Computational Horizons In Cancer (CHIC) (297) launched in 2013, aimed to provide an oncosimulator modeling platform for in silico oncology. CHIC was initiated to develop and implement predictive oncology and individualized multi-scale cancer modeling tools towards assisting quantitative cancer research and personalized cancer therapy. The workflow and tools established under the ambient of the CHIC project allowed the development of robust, interoperable, and collaborative in silico models in cancer and normal conditions. The CHIC project also proposed a pipeline to translate the model towards supporting clinicians to make optimal personalized treatment plans for individual patients. To this end, several models are created using CHIC such as non-small cell lung cancer (326), glioblastoma (327), leukemia model (328), etc. These models furnish a quantitative understanding of tumorigenesis towards providing avenues for promoting individual cancer patient treatment combinations. One such model reported by Kolokotroni et al. (326) evaluated the efficacy of cisplatin-based therapy for NSCLC patients using in silico multi-scale cancer model. While, in 2015, Antonopoulos and Stamatakos (327) modeled the infiltration of glioblastoma cells in normal brain regions using a novel treatment. In 2017, Stamatakos and Giatili (329) extended glioblastoma oncosimulator for modeling tumor growth using reaction-diffusion numerical handling based on multi-scale panomics data. They further proposed a clinical pipeline to translate the model into clinical settings towards supporting clinicians to make optimal personalized treatment plans for individual patients. Ouzounoglou et al. (328) developed an in silico multi-scale leukemia oncosimulator model towards modeling deregulations in the G1/S pathway to investigate the altered function of retinoblastoma in ALL patients. However, a clinical translation of these models is currently in the works (297).
6 Discussion
In this work, we have evaluated the use of data-driven multi-scale cancer models in deciphering complex biomolecular underpinnings in cancer research towards developing personalized therapeutics interventions for clinical decision-making. Specifically, we have discussed the chronological evolution of online cancer data repositories that host high-resolution datasets from multiple spatiotemporal scales. Next, we evaluate how this data drives single- and multi-scale systems biology models towards decoding complex cancer regulation in patients. We then track the development of various modeling software and associated applications in enhancing the translational role of cancer systems biology models in clinics.
We conclude that the contemporary multi-scale modeling software line-up remains limited in their clinical employment due to the lack of a generic, zero-code, panomic-based framework for translating research models into clinical settings. Such a framework would help annotate in silico cancer models developed using single and multi-scale databases (45, 61, 330–334). The framework should also provide an environment for developing the extra-cellular matrix of a cancer cell which can then be integrated into cellular models. Existing environment (38, 39) and cell line databases (41–43) need to be integrated to design environmental models along with biologically plausible cell line structures. These cell lines could then be assembled into three-dimensional geometries to create multi-scale in silico organoids. The pipeline should also furnish capabilities such as a convenient import workflow for clinical data integration through histopathological image data repositories (44–46) for designing biologically accurate organoid structures based on each cancer patient’s underlying cellular morphology. Once the personalized multi-scale model has been constructed, the pipeline would allow investigation into the temporal evolution of the multi-scale organoid under personalized inputs and user-designed biomolecular entities. Data generated from the multi-scale model simulation can be analyzed to elicit biomolecular cues for each cancer patient as well as determine its role. Due to the heterogeneity and complexity of biomolecular data a coherent panomics-based pipeline can be challenging to develop and will require a collaborative effort by various research groups, through close collaborations and data standardization.
Taken together, a translational in silico systems oncology pipeline is the need of the hour and will help develop and deliver personalized treatments of cancer as well as substantively inform clinical decision-making processes.
Author Contributions
SC and MG carried out the literature review and drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National ICT-R&D Fund (SRG-209), RF-NCBC-015, NGIRI-2020-4771, HEC (21-30SRGP/R&D/HEC/2014, 20-2269/NRPU/R&D/HEC/12/4792 and 20-3629/NRPU/R&D/HEC/14/585), TWAS (RG 14-319 RG/ITC/AS_C) and LUMS (STG-BIO-1008, FIF-BIO-2052, FIF-BIO-0255, SRP-185-BIO, SRP-058-BIO and FIF-477-1819-BIO) grants.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
This manuscript has been released as a pre-print at bioRxiv (335).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.712505/full#supplementary-material
References
1. DeVita VT. The “War on Cancer” and its Impact. Nat Clin Pract Oncol (2004) 1(2):55. doi: 10.1038/ncponc0036
2. Wong AH-H, Deng C-X. Precision Medicine for Personalized Cancer Therapy. Int J Biol Sci (2015) 11(12):1410–2. doi: 10.7150/ijbs.14154
3. Giacomini KM, Yee SW, Ratain MJ, Weinshilboum RM, Kamatani N, Nakamura Y. Pharmacogenomics and Patient Care: One Size Does Not Fit All. Sci Transl Med (2012) 4(153):153ps18. doi: 10.1126/scitranslmed.3003471
4. Krzyszczyk P, Acevedo A, Davidoff EJ, Timmins LM, Marrero-Berrios I, Patel M, et al. The Growing Role of Precision and Personalized Medicine for Cancer Treatment. Technol (Singap World Sci) (2018) 6(3–4):79–100. doi: 10.1142/s2339547818300020
5. Binotto M, Reinert T, Werutsky G, Zaffaroni F, Schwartsmann G. Health-Related Quality of Life Before and During Chemotherapy in Patients With Early-Stage Breast Cancer. Ecancermedicalscience (2020) 14:1007. doi: 10.3332/ecancer.2020.1007
6. Adler AS, Lin M, Horlings H, Nuyten DSA, van de Vijver MJ, Chang HY. Genetic Regulators of Large-Scale Transcriptional Signatures in Cancer. Nat Genet (2006) 38(4):421–30. doi: 10.1038/ng1752
7. Shi Y, Wei J, Chen Z, Yuan Y, Li X, Zhang Y, et al. Integrative Analysis Reveals Comprehensive Altered Metabolic Genes Linking With Tumor Epigenetics Modification in Pan-Cancer. BioMed Res Int (2019) 2019:1–18. doi: 10.1155/2019/6706354
8. Chen Y, Ni J, Gao Y, Zhang J, Liu X, Chen Y, et al. Integrated Proteomics and Metabolomics Reveals the Comprehensive Characterization of Antitumor Mechanism Underlying Shikonin on Colon Cancer Patient-Derived Xenograft Model. Sci Rep (2020) 10(1):1–14. doi: 10.1038/s41598-020-71116-5
9. Angermueller C, Clark SJ, Lee HJ, Macaulay IC, Teng MJ, Hu TX, et al. Parallel Single-Cell Sequencing Links Transcriptional and Epigenetic Heterogeneity. Nat Methods (2016) 13(3):229–32. doi: 10.1038/nmeth.3728
10. Yoshida A, Kohyama S, Fujiwara K, Nishikawa S, Doi N. Regulation of Spatiotemporal Patterning in Artificial Cells by a Defined Protein Expression System. Chem Sci (2019) 10(48):11064–72. doi: 10.1039/c9sc02441g
11. Rejniak KA, Anderson ARA. Hybrid Models of Tumor Growth. Wiley Interdiscip Rev Syst Biol Med (2011) 3(1):115–25. doi: 10.1002/wsbm.102
12. Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, et al. Oncogenic Signaling Pathways in The Cancer Genome Atlas. Cell (2018) 173(2):321–37. doi: 10.1016/j.cell.2018.03.035
13. Zehir A, Benayed R, Shah RH, Syed A, Middha S, Kim HR, et al. Mutational Landscape of Metastatic Cancer Revealed From Prospective Clinical Sequencing of 10,000 Patients. Nat Med (2017) 23(6):703–13. doi: 10.1038/nm.4333
14. Al Hajri Q, Dash S, Feng W, Garner HR, Anandakrishnan R. Identifying Multi-Hit Carcinogenic Gene Combinations: Scaling Up a Weighted Set Cover Algorithm Using Compressed Binary Matrix Representation on a GPU. Sci Rep (2020) 10(1):1–18. doi: 10.1038/s41598-020-58785-y
15. Roerink SF, Sasaki N, Lee-Six H, Young MD, Alexandrov LB, Behjati S, et al. Intra-Tumour Diversification in Colorectal Cancer at the Single-Cell Level. Nature (2018) 556(7702):457–62. doi: 10.1038/s41586-018-0024-3
16. Song Q, Merajver SD, Li JZ. Cancer Classification in the Genomic Era: Five Contemporary Problems. Hum Genomics (2015) 9(1):1–8. doi: 10.1186/s40246-015-0049-8
17. Hanahan D, Weinberg RA. The Hallmarks of Cancer. Cell (2000) 100(1):57–70. doi: 10.1016/s0092-8674(00)81683-9
18. Hanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013
19. Saha R, Chowdhury A, Maranas CD. Recent Advances in the Reconstruction of Metabolic Models and Integration of Omics Data. Curr Opin Biotechnol (2014) 29:39–45. doi: 10.1016/j.copbio.2014.02.011
20. Deisboeck TS, Berens ME, Kansal AR, Torquato S, Stemmer-Rachamimov AO, Chiocca EA. Pattern of Self-Organization in Tumour Systems: Complex Growth Dynamics in a Novel Brain Tumour Spheroid Model. Cell Prolif (2001) 34(2):115–34. doi: 10.1046/j.1365-2184.2001.00202.x
21. Lorenzo G, Scott MA, Tew K, Hughes TJR, Zhang YJ, Liu L, et al. Tissue-Scale, Personalized Modeling and Simulation of Prostate Cancer Growth. Proc Natl Acad Sci (2016) 113(48):E7663–71. doi: 10.1073/pnas.1615791113
22. Pingle SC, Sultana Z, Pastorino S, Jiang P, Mukthavaram R, Chao Y, et al. In Silico Modeling Predicts Drug Sensitivity of Patient-Derived Cancer Cells. J Transl Med (2014) 12(1):1–13. doi: 10.1186/1479-5876-12-128
23. Zhang J, Chiodini R, Badr A, Zhang G. The Impact of Next-Generation Sequencing on Genomics. J Genet Genomics (2011) 38(3):95–109. doi: 10.1016/j.jgg.2011.02.003
24. Velculescu VE, Zhang L, Zhou W, Vogelstein J, Basrai MA, Bassett DEJ, et al. Characterization of the Yeast Transcriptome. Cell (1997) 88(2):243–51. doi: 10.1016/s0092-8674(00)81845-0
25. Anderson NL, Anderson NG. Proteome and Proteomics: New Technologies, New Concepts, and New Words. Electrophoresis (1998) 19(11):1853–61. doi: 10.1002/elps.1150191103
26. Manzoni C, Kia DA, Vandrovcova J, Hardy J, Wood NW, Lewis PA, et al. Genome, Transcriptome and Proteome: The Rise of Omics Data and Their Integration in Biomedical Sciences. Brief Bioinform (2016) 19(2):286–302. doi: 10.1093/bib/bbw114
27. Stamatakos GS, Dionysiou DD, Graf NM, Sofra NA, Desmedt C, Hoppe A, et al. The “Oncosimulator”: A Multilevel, Clinically Oriented Simulation System of Tumor Growth and Organism Response to Therapeutic Schemes. Towards the Clinical Evaluation of In Silico Oncology. Conf Proc IEEE Eng Med Biol Soc (2007), 6629–32. doi: 10.1109/IEMBS.2007.4353879
28. Stephens ZD, Lee SY, Faghri F, Campbell RH, Zhai C, Efron MJ, et al. Big Data: Astronomical or Genomical? PloS Biol (2015) 13(7):e1002195. doi: 10.1371/journal.pbio.1002195
29. Krupp M, Marquardt JU, Sahin U, Galle PR, Castle J, Teufel A. RNA-Seq Atlas–a Reference Database for Gene Expression Profiling in Normal Tissue by Next-Generation Sequencing. Bioinformatics (2012) 28(8):1184–5. doi: 10.1093/bioinformatics/bts084
30. Park S-J, Yoon B-H, Kim S-K, Kim S-Y. GENT2: An Updated Gene Expression Database for Normal and Tumor Tissues. BMC Med Genomics (2019) 12(5):1–8. doi: 10.1186/s12920-019-0514-7
31. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res (2002) 30(1):207–10. doi: 10.1093/nar/30.1.207
32. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res (1999) 27(1):29–34. doi: 10.1093/nar/27.1.29
33. Nikitin A, Egorov S, Daraselia N, Mazo I. Pathway Studio–the Analysis and Navigation of Molecular Networks. Bioinformatics (2003) 19(16):2155–7. doi: 10.1093/bioinformatics/btg290
34. Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: A Library of Protein Families and Subfamilies Indexed by Function. Genome Res (2003) 13(9):2129–41. doi: 10.1101/gr.772403
35. Le Novère N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, et al. BioModels Database: A Free, Centralized Database of Curated, Published, Quantitative Kinetic Models of Biochemical and Cellular Systems. Nucleic Acids Res (2006) 34(suppl_1):D689–91. doi: 10.1093/nar/gkj092
36. Hermjakob H, Montecchi-Palazzi L, Lewington C, Mudali S, Kerrien S, Orchard S, et al. IntAct: An Open Source Molecular Interaction Database. Nucleic Acids Res (2004) 32(suppl_1):D452–5. doi: 10.1093/nar/gkh052
37. Zanzoni A, Montecchi-Palazzi L, Quondam M, Ausiello G, Helmer-Citterich M, Cesareni G. MINT: A Molecular INTeraction Database. FEBS Lett (2002) 513(1):135–40. doi: 10.1016/s0014-5793(01)03293-8
38. Chautard E, Fatoux-Ardore M, Ballut L, Thierry-Mieg N, Ricard-Blum S. MatrixDB, the Extracellular Matrix Interaction Database. Nucleic Acids Res (2011) 39(suppl_1):D235–40. doi: 10.1093/nar/gkq830
39. Shao X, Taha IN, Clauser KR, Gao YT, Naba A. MatrisomeDB: The ECM-Protein Knowledge Database. Nucleic Acids Res (2020) 48(D1):D1136–44. doi: 10.1093/nar/gkz849
40. Tate JG, Bamford S, Jubb HC, Sondka Z, Beare DM, Bindal N, et al. COSMIC: The Catalogue Of Somatic Mutations In Cancer. Nucleic Acids Res (2019) 47(D1):D941–7. doi: 10.1093/nar/gky1015
41. Rajapakse VN, Luna A, Yamade M, Loman L, Varma S, Sunshine M, et al. CellMinerCDB for Integrative Cross-Database Genomics and Pharmacogenomics Analyses of Cancer Cell Lines. iScience (2018) 10:247–64. doi: 10.1016/j.isci.2018.11.029
42. Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia Enables Predictive Modelling of Anticancer Drug Sensitivity. Nature (2012) 483(7391):603–7. doi: 10.1038/nature11003
43. Romano P, Manniello A, Aresu O, Armento M, Cesaro M, Parodi B. Cell Line Data Base: Structure and Recent Improvements Towards Molecular Authentication of Human Cell Lines. Nucleic Acids Res (2009) 37(suppl_1):D925–32. doi: 10.1093/nar/gkn730
44. Gutman DA, Cobb J, Somanna D, Park Y, Wang F, Kurc T, et al. Cancer Digital Slide Archive: An Informatics Resource to Support Integrated In Silico Analysis of TCGA Pathology Data. J Am Med Inform Assoc (2013) 20(6):1091–8. doi: 10.1136/amiajnl-2012-001469
45. Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The Genotype-Tissue Expression (GTEx) Project. Nat Genet (2013) 45(6):580–5. doi: 10.1038/ng.2653
46. Clark K, Vendt B, Smith K, Freymann J, Kirby J, Koppel P, et al. The Cancer Imaging Archive (TCIA): Maintaining and Operating a Public Information Repository. J Digit Imaging (2013) 26(6):1045–57. doi: 10.1007/s10278-013-9622-7
47. Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. KEGG for Representation and Analysis of Molecular Networks Involving Diseases and Drugs. Nucleic Acids Res (2010) 38(suppl_1):D355–60. doi: 10.1093/nar/gkp896
48. Hamosh A, Scott AF, Amberger J, Bocchini C, Valle D, McKusick VA. Online Mendelian Inheritance in Man (OMIM), a Knowledgebase of Human Genes and Genetic Disorders. Nucleic Acids Res (2002) 30(1):52–5. doi: 10.1093/nar/30.1.52
49. Chakravarty D, Gao J, Phillips SM, Kundra R, Zhang H, Wang J, et al. OncoKB: A Precision Oncology Knowledge Base. JCO Precis Oncol (2017) 1:1–16. doi: 10.1200/PO.17.00011
50. Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, et al. DrugBank: A Knowledgebase for Drugs, Drug Actions and Drug Targets. Nucleic Acids Res (2008) 36(suppl_1):D901–6. doi: 10.1093/nar/gkm958
51. Smirnov P, Kofia V, Maru A, Freeman M, Ho C, El-Hachem N, et al. PharmacoDB: An Integrative Database for Mining In Vitro Anticancer Drug Screening Studies. Nucleic Acids Res (2018) 46(D1):D994–1002. doi: 10.1093/nar/gkx911
52. Piñeiro-Yáñez E, Reboiro-Jato M, Gómez-López G, Perales-Patón J, Troulé K, Rodríguez JM, et al. PanDrugs: A Novel Method to Prioritize Anticancer Drug Treatments According to Individual Genomic Data. Genome Med (2018) 10(1):1–11. doi: 10.1186/s13073-018-0546-1
53. Pinu FR, Beale DJ, Paten AM, Kouremenos K, Swarup S, Schirra HJ, et al. Systems Biology and Multi-Omics Integration: Viewpoints From the Metabolomics Research Community. Metabolites (2019) 9(4):1–31. doi: 10.3390/metabo9040076
54. Bilofsky HS, Burks C, Fickett JW, Goad WB, Wayne FIL, Rindone P, et al. The GenBank Genetic Sequence Databank. Nucleic Acids Res (1986) 14(1):1–4. doi: 10.1093/nar/14.1.1
55. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res (2012) 41(D1):D36–42. doi: 10.1093/nar/gks1195
56. Benson DA, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res (2015) 43(D1):D30–5. doi: 10.1093/nar/gku1216
57. Díez O, Osorio A, Durán M, Martinez-Ferrandis JI, de la Hoya M, Salazar R, et al. Analysis of BRCA1 and BRCA2 Genes in Spanish Breast/Ovarian Cancer Patients: A High Proportion of Mutations Unique to Spain and Evidence of Founder Effects. Hum Mutat (2003) 22(4):301–12. doi: 10.1002/humu.10260
58. Medrek C, Pontén F, Jirström K, Leandersson K. The Presence of Tumor Associated Macrophages in Tumor Stroma as a Prognostic Marker for Breast Cancer Patients. BMC Cancer (2012) 12(1):1–9. doi: 10.1186/1471-2407-12-306
59. Hubbard T, Barker D, Birney E, Cameron G, Chen Y, Clark L, et al. The Ensembl Genome Database Project. Nucleic Acids Res (2002) 30(1):38–41. doi: 10.1093/nar/30.1.38
60. Easton DF, Pooley KA, Dunning AM, Pharoah PDP, Thompson D, Ballinger DG, et al. Genome-Wide Association Study Identifies Novel Breast Cancer Susceptibility Loci. Nature (2007) 447(7148):1087–93. doi: 10.1038/nature05887
61. Bamford S, Dawson E, Forbes S, Clements J, Pettett R, Dogan A, et al. The COSMIC (Catalogue of Somatic Mutations in Cancer) Database and Website. Br J Cancer (2004) 91:355–8. doi: 10.1038/sj.bjc.6601894
62. Schubbert S, Shannon K, Bollag G. Hyperactive Ras in Developmental Disorders and Cancer. Nat Rev Cancer (2007) 7(4):295–308. doi: 10.1038/nrc2109
63. Weir BA, Woo MS, Getz G, Perner S, Ding L, Beroukhim R, et al. Characterizing the Cancer Genome in Lung Adenocarcinoma. Nature (2007) 450(7171):893–8. doi: 10.1038/nature06358
64. Zhao D, Tang Y, Xia X, Sun J, Meng J, Shang J, et al. Integration of Transcriptome, Proteome, and Metabolome Provides Insights Into How Calcium Enhances the Mechanical Strength of Herbaceous Peony Inflorescence Stems. Cells (2019) 8(2):1–22. doi: 10.3390/cells8020102
65. Chakraborty PK, Mustafi SB, Xiong X, Dwivedi SKD, Nesin V, Saha S, et al. MICU1 Drives Glycolysis and Chemoresistance in Ovarian Cancer. Nat Commun (2017) 8(1):1–16. doi: 10.1038/ncomms14634
66. Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): An Immeasurable Source of Knowledge. Contemp Oncol (Pozn) (2015) 19(1A):A68–77. doi: 10.5114/wo.2014.47136
67. Tamborero D, Gonzalez-Perez A, Perez-Llamas C, Deu-Pons J, Kandoth C, Reimand J, et al. Comprehensive Identification of Mutational Cancer Driver Genes Across 12 Tumor Types. Sci Rep (2013) 3(1):1–10. doi: 10.1038/srep02650
68. Leiserson MDM, Vandin F, Wu H-T, Dobson JR, Eldridge JV, Thomas JL, et al. Pan-Cancer Network Analysis Identifies Combinations of Rare Somatic Mutations Across Pathways and Protein Complexes. Nat Genet (2015) 47(2):106–14. doi: 10.1038/ng.3168
69. Davis CF, Ricketts CJ, Wang M, Yang L, Cherniack AD, Shen H, et al. The Somatic Genomic Landscape of Chromophobe Renal Cell Carcinoma. Cancer Cell (2014) 26(3):319–30. doi: 10.1016/j.ccr.2014.07.014
70. Zhang J, Baran J, Cros A, Guberman JM, Haider S, Hsu J, et al. International Cancer Genome Consortium Data Portal—a One-Stop Shop for Cancer Genomics Data. Database (2011) 2011:1–10. doi: 10.1093/database/bar026
71. Jones DTW, Jäger N, Kool M, Zichner T, Hutter B, Sultan M, et al. Dissecting the Genomic Complexity Underlying Medulloblastoma. Nature (2012) 488(7409):100–5. doi: 10.1038/nature11284
72. Burns MB, Lackey L, Carpenter MA, Rathore A, Land AM, Leonard B, et al. APOBEC3B Is an Enzymatic Source of Mutation in Breast Cancer. Nature (2013) 494(7437):366–70. doi: 10.1038/nature11881
73. Supek F, Miñana B, Valcárcel J, Gabaldón T, Lehner B. Synonymous Mutations Frequently Act as Driver Mutations in Human Cancers. Cell (2014) 156(6):1324–35. doi: 10.1016/j.cell.2014.01.051
74. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The Cbio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discovery (2012) 2(5):401–4. doi: 10.1158/2159-8290.CD-12-0095
75. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the Cbioportal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088
76. Robinson D, Allen EMV, Wu Y-M, Schultz N, Lonigro RJ, Mosquera J-M, et al. Integrative Clinical Genomics of Advanced Prostate Cancer. Cell (2015) 161(5):1215–28. doi: 10.1016/j.cell.2015.05.001
77. Zaretsky JM, Garcia-Diaz A, Shin DS, Escuin-Ordinas H, Hugo W, Hu-Lieskovan S, et al. Mutations Associated With Acquired Resistance to PD-1 Blockade in Melanoma. N Engl J Med (2016) 375(9):819–29. doi: 10.1056/NEJMoa1604958
78. Bailey P, Chang DK, Nones K, Johns AL, Patch A-M, Gingras M-C, et al. Genomic Analyses Identify Molecular Subtypes of Pancreatic Cancer. Nature (2016) 531(7592):47–52. doi: 10.1038/nature16965
79. Jiao X-D, Qin B-D, You P, Cai J, Zang Y-S. The Prognostic Value of TP53 and its Correlation With EGFR Mutation in Advanced non-Small Cell Lung Cancer, an Analysis Based on Cbioportal Data Base. Lung Cancer (2018) 123:70–5. doi: 10.1016/j.lungcan.2018.07.003
80. Hou H, Zhang C, Qi X, Zhou L, Liu D, Lv H, et al. Distinctive Targetable Genotypes of Younger Patients With Lung Adenocarcinoma: A Cbioportal for Cancer Genomics Data Base Analysis. Cancer Biol Ther (2020) 21(1):26–33. doi: 10.1080/15384047.2019.1665392
81. Chen G, Gharib TG, Huang C-C, Taylor JMG, Misek DE, Kardia SLR, et al. Discordant Protein and mRNA Expression in Lung Adenocarcinomas. Mol Cell Proteomics (2002) 1(4):304–13. doi: 10.1074/mcp.m200008-mcp200
82. Cox J, Mann M. Is Proteomics the New Genomics? Cell (2007) 130(3):395–8. doi: 10.1016/j.cell.2007.07.032
83. Aslam B, Basit M, Nisar MA, Khurshid M, Rasool MH. Proteomics: Technologies and Their Applications. J Chromatogr Sci (2017) 55(2):182–96. doi: 10.1093/chromsci/bmw167
84. Williams CW, Elmendorf HG. Identification and Analysis of the RNA Degrading Complexes and Machinery of Giardia Lamblia Using an In Silico Approach. BMC Genomics (2011) 12(1):1–15. doi: 10.1186/1471-2164-12-586
85. Gill G. Post-Translational Modification by the Small Ubiquitin-Related Modifier SUMO has Big Effects on Transcription Factor Activity. Curr Opin Genet Dev (2003) 13(2):108–13. doi: 10.1016/S0959-437X(03)00021-2
86. Bathke J, Konzer A, Remes B, McIntosh M, Klug G. Comparative Analyses of the Variation of the Transcriptome and Proteome of Rhodobacter Sphaeroides Throughout Growth. BMC Genomics (2019) 20(1):1–3. doi: 10.1186/s12864-019-5749-3
88. Persson A, Hober S, Uhlén M. A Human Protein Atlas Based on Antibody Proteomics. Curr Opin Mol Ther (2006) 8(3):185–90.
89. Uhlén M, Björling E, Agaton C, Szigyarto CAK, Amini B, Andersen E, et al. A Human Protein Atlas for Normal and Cancer Tissues Based on Antibody Proteomics. Mol Cell Proteomics (2005) 4(12):1920–32. doi: 10.1074/mcp.M500279-MCP200
90. Gámez-Pozo A, Sánchez-Navarro I, Nistal M, Calvo E, Madero R, Díaz E, et al. MALDI Profiling of Human Lung Cancer Subtypes. PloS One (2009) 4(11):e7731. doi: 10.1371/journal.pone.0007731
91. Imberg-Kazdan K, Ha S, Greenfield A, Poultney CS, Bonneau R, Logan SK, et al. A Genome-Wide RNA Interference Screen Identifies New Regulators of Androgen Receptor Function in Prostate Cancer Cells. Genome Res (2013) 23(4):581–91. doi: 10.1101/gr.144774.112
92. Collins FS. New Goals for the U.S. Human Genome Project: 1998-2003. Science (1998) 282(5389):682–9. doi: 10.1126/science.282.5389.682
93. Legrain P, Aebersold R, Archakov A, Bairoch A, Bala K, Beretta L, et al. The Human Proteome Project: Current State and Future Direction. Mol Cell Proteomics (2011) 10(7):M111.009993. doi: 10.1074/mcp.M111.009993
94. Cottingham K. Hupo’s Human Proteome Project: The Next Big Thing? J Proteome Res (2008) 7: (6):2192. doi: 10.1021/pr0837441
95. Wang Q, Wen B, Yan G, Wei J, Xie L, Xu S, et al. Qualitative and Quantitative Expression Status of the Human Chromosome 20 Genes in Cancer Tissues and the Representative Cell Lines. J Proteome Res (2013) 12(1):151–61. doi: 10.1021/pr3008336
96. Hüttenhain R, Soste M, Selevsek N, Röst H, Sethi A, Carapito C, et al. Reproducible Quantification of Cancer-Associated Proteins in Body Fluids Using Targeted Proteomics. Sci Transl Med (2012) 4(142):142ra94–142ra94. doi: 10.1126/scitranslmed.3003989
97. Zhang B, Wang J, Wang X, Zhu J, Liu Q, Shi Z, et al. Proteogenomic Characterization of Human Colon and Rectal Cancer. Nature (2014) 513(7518):382–7. doi: 10.1038/nature13438
98. Huang Z, Ma L, Huang C, Li Q, Nice EC. Proteomic Profiling of Human Plasma for Cancer Biomarker Discovery. Proteomics (2017) 17(6):1–13. doi: 10.1002/pmic.201600240
99. Mertins P, Mani DR, Ruggles KV, Gillette MA, Clauser KR, Wang P, et al. Proteogenomics Connects Somatic Mutations to Signalling in Breast Cancer. Nature (2016) 534(7605):55–62. doi: 10.1038/nature18003
100. Zhang H, Liu T, Zhang Z, Payne SH, Zhang B, McDermott JE, et al. Integrated Proteogenomic Characterization of Human High-Grade Serous Ovarian Cancer. Cell (2016) 166(3):755–65. doi: 10.1016/j.cell.2016.05.069
101. Sebastian S, Bulun SE. A Highly Complex Organization of the Regulatory Region of the Human CYP19 (Aromatase) Gene Revealed by the Human Genome Project. J Clin Endocrinol Metab (2001) 86(10):4600–2. doi: 10.1210/jcem.86.10.7947
102. Aebersold R, Bader GD, Edwards AM, Van Eyk JE, Kussmann M, Qin J, et al. The Biology/Disease-Driven Human Proteome Project (B/D-HPP): Enabling Protein Research for the Life Sciences Community. J Proteome Res (2013) 12(1):23–7. doi: 10.1021/pr301151m
103. Paik Y-K, Jeong S-K, Omenn GS, Uhlen M, Hanash S, Cho SY, et al. The Chromosome-Centric Human Proteome Project for Cataloging Proteins Encoded in the Genome. Nat Biotechnol (2012) 30(3):221–3. doi: 10.1038/nbt.2152
104. Gupta MK, Jayaram S, Madugundu AK, Chavan S, Advani J, Pandey A, et al. Chromosome-Centric Human Proteome Project: Deciphering Proteins Associated With Glioma and Neurodegenerative Disorders on Chromosome 12. J Proteome Res (2014) 13(7):3178–90. doi: 10.1021/pr500023p
105. Fan W, Ge G, Liu Y, Wang W, Liu L, Jia Y. Proteomics Integrated With Metabolomics: Analysis of the Internal Causes of Nutrient Changes in Alfalfa at Different Growth Stages. BMC Plant Biol (2018) 18(1):1–15. doi: 10.1186/s12870-018-1291-8
106. Kopka J, Schauer N, Krueger S, Birkemeyer C, Usadel B, Bergmüller E, et al. GMD@CSB.DB: The Golm Metabolome Database. Bioinformatics (2005) 21(8):1635–8. doi: 10.1093/bioinformatics/bti236
107. Wishart DS, Feunang YD, Marcu A, Guo AC, Liang K, Vázquez-Fresno R, et al. HMDB 4.0: The Human Metabolome Database for 2018. Nucleic Acids Res (2018) 46(D1):D608–17. doi: 10.1093/nar/gkx1089
108. Sud M, Fahy E, Cotter D, Azam K, Vadivelu I, Burant C, et al. Metabolomics Workbench: An International Repository for Metabolomics Data and Metadata, Metabolite Standards, Protocols, Tutorials and Training, and Analysis Tools. Nucleic Acids Res (2016) 44(D1):D463–70. doi: 10.1093/nar/gkv1042
109. Wedge DC, Allwood JW, Dunn W, Vaughan AA, Simpson K, Brown M, et al. Is Serum or Plasma More Appropriate for Intersubject Comparisons in Metabolomic Studies? An Assessment in Patients With Small-Cell Lung Cancer. Anal Chem (2011) 83(17):6689–97. doi: 10.1021/ac2012224
110. Pasikanti KK, Esuvaranathan K, Hong Y, Ho PC, Mahendran R, Raman Nee Mani L, et al. Urinary Metabotyping of Bladder Cancer Using Two-Dimensional Gas Chromatography Time-of-Flight Mass Spectrometry. J Proteome Res (2013) 12(9):3865–73. doi: 10.1021/pr4000448
111. Wishart DS, Tzur D, Knox C, Eisner R, Guo AC, Young N, et al. HMDB: The Human Metabolome Database. Nucleic Acids Res (2007) 35(suppl_1):D521–6. doi: 10.1093/nar/gkl923
112. Sugimoto M, Wong DT, Hirayama A, Soga T, Tomita M. Capillary Electrophoresis Mass Spectrometry-Based Saliva Metabolomics Identified Oral, Breast and Pancreatic Cancer-Specific Profiles. Metabolomics (2010) 6(1):78–95. doi: 10.1007/s11306-009-0178-y
113. Agren R, Bordel S, Mardinoglu A, Pornputtapong N, Nookaew I, Nielsen J. Reconstruction of Genome-Scale Active Metabolic Networks for 69 Human Cell Types and 16 Cancer Types Using INIT. PloS Comput Biol (2012) 8(5):e1002518. doi: 10.1371/journal.pcbi.1002518
114. Hattori A, Tsunoda M, Konuma T, Kobayashi M, Nagy T, Glushka J, et al. Cancer Progression by Reprogrammed BCAA Metabolism in Myeloid Leukaemia. Nature (2017) 545(7655):500–4. doi: 10.1038/nature22314
115. Mi H, Muruganujan A, Casagrande JT, Thomas PD. Large-Scale Gene Function Analysis With the PANTHER Classification System. Nat Protoc (2013) 8(8):1551–66. doi: 10.1038/nprot.2013.092
116. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New Perspectives on Genomes, Pathways, Diseases and Drugs. Nucleic Acids Res (2017) 45(D1):D353–61. doi: 10.1093/nar/gkw1092
117. Li X-N, Wang Z-J, Ye C-X, Zhao B-C, Li Z-L, Yang Y. RNA Sequencing Reveals the Expression Profiles of circRNA and Indicates That Circddx17 Acts as a Tumor Suppressor in Colorectal Cancer. J Exp Clin Cancer Res (2018) 37(1):1–14. doi: 10.1186/s13046-018-1006-x
118. Feng H, Gu Z-Y, Li Q, Liu Q-H, Yang X-Y, Zhang J-J. Identification of Significant Genes With Poor Prognosis in Ovarian Cancer via Bioinformatical Analysis. J Ovarian Res (2019) 12(1):1–9. doi: 10.1186/s13048-019-0508-2
119. Stenson PD, Ball EV, Mort M, Phillips AD, Shaw K, Cooper DN. The Human Gene Mutation Database (HGMD) and its Exploitation in the Fields of Personalized Genomics and Molecular Evolution. Curr Protoc Bioinforma (2012) 39(1):1–13. doi: 10.1002/0471250953.bi0113s39
120. Turcan S, Rohle D, Goenka A, Walsh LA, Fang F, Yilmaz E, et al. IDH1 Mutation is Sufficient to Establish the Glioma Hypermethylator Phenotype. Nature (2012) 483(7390):479–83. doi: 10.1038/nature10866
121. Karp PD, Riley M, Paley SM, Pellegrini-Toole A. The MetaCyc Database. Nucleic Acids Res (2002) 30(1):59–61. doi: 10.1093/nar/30.1.59
122. Miller JA, Pappan K, Thompson PA, Want EJ, Siskos AP, Keun HC, et al. Plasma Metabolomic Profiles of Breast Cancer Patients After Short-Term Limonene Intervention. Cancer Prev Res (Phila) (2015) 8(1):86–93. doi: 10.1158/1940-6207.CAPR-14-0100
123. Bader GD, Donaldson I, Wolting C, Ouellette BF, Pawson T, Hogue CW. BIND–The Biomolecular Interaction Network Database. Nucleic Acids Res (2001) 29(1):242–5. doi: 10.1093/nar/29.1.242
124. Chen Y, Choong L-Y, Lin Q, Philp R, Wong C-H, Ang B-K, et al. Differential Expression of Novel Tyrosine Kinase Substrates During Breast Cancer Development. Mol Cell Proteomics (2007) 6(12):2072–87. doi: 10.1074/mcp.M700395-MCP200
125. Vinayagam A, Gibson TE, Lee HJ, Yilmazel B, Roesel C, Hu Y, et al. Controllability Analysis of the Directed Human Protein Interaction Network Identifies Disease Genes and Drug Targets. Proc Natl Acad Sci USA (2016) 113(18):4976–81. doi: 10.1073/pnas.1603992113
126. Xenarios I, Rice DW, Salwinski L, Baron MK, Marcotte EM, Eisenberg D. DIP: The Database of Interacting Proteins. Nucleic Acids Res (2000) 28(1):289–91. doi: 10.1093/nar/28.1.289
127. Goh WWB, Lee YH, Zubaidah RM, Jin J, Dong D, Lin Q, et al. Network-Based Pipeline for Analyzing MS Data: An Application Toward Liver Cancer. J Proteome Res (2011) 10(5):2261–72. doi: 10.1021/pr1010845
128. Fu LL, Zhao X, Xu HL, Wen X, Wang SY, Liu B, et al. Identification of microRNA-Regulated Autophagic Pathways in Plant Lectin-Induced Cancer Cell Death. Cell Prolif (2012) 45(5):477–85. doi: 10.1111/j.1365-2184.2012.00840.x
129. Graeber TG, Eisenberg D. Bioinformatic Identification of Potential Autocrine Signaling Loops in Cancers From Gene Expression Profiles. Nat Genet (2001) 29(3):295–300. doi: 10.1038/ng755
130. Von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, et al. STRING: Known and Predicted Protein–Protein Associations, Integrated and Transferred Across Organisms. Nucleic Acids Res (2005) 33(suppl_1):D433–7. doi: 10.1093/nar/gki005
131. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING Database in 2017: Quality-Controlled Protein-Protein Association Networks, Made Broadly Accessible. Nucleic Acids Res (2017) 45(D1):D362–8. doi: 10.1093/nar/gkw937
132. Mlecnik B, Tosolini M, Charoentong P, Kirilovsky A, Bindea G, Berger A, et al. Biomolecular Network Reconstruction Identifies T-Cell Homing Factors Associated With Survival in Colorectal Cancer. Gastroenterology (2010) 138(4):1429–40. doi: 10.1053/j.gastro.2009.10.057
133. Azpeitia E, Balanzario EP, Wagner A. Signaling Pathways Have an Inherent Need for Noise to Acquire Information. BMC Bioinf (2020) 21(1):1–21. doi: 10.1186/s12859-020-03778-x
134. Peng R, Xu L, Wang H, Lyu Y, Wang D, Bi C, et al. DNA-Based Artificial Molecular Signaling System That Mimics Basic Elements of Reception and Response. Nat Commun (2020) 11(1):1–10. doi: 10.1038/s41467-020-14739-6
135. Suderman R, Deeds EJ. Intrinsic Limits of Information Transmission in Biochemical Signalling Motifs. Interface Focus (2018) 8(6):1–15. doi: 10.1098/rsfs.2018.0039
136. Quail DF, Joyce JA. Microenvironmental Regulation of Tumor Progression and Metastasis. Nat Med (2013) 19(11):1423–37. doi: 10.1038/nm.3394
137. Thul PJ, Lindskog C. The Human Protein Atlas: A Spatial Map of the Human Proteome. Protein Sci (2018) 27(1):233–44. doi: 10.1002/pro.3307
138. Sehuler GD. Pieces of Use Puzzle: Expressed Sequence Tags and the Catalog of Human Genes. J Mol Med (1997) 75(10):694–8. doi: 10.1007/s001090050155
139. Celik S, Logsdon BA, Battle S, Drescher CW, Rendi M, Hawkins RD, et al. Extracting a Low-Dimensional Description of Multiple Gene Expression Datasets Reveals a Potential Driver for Tumor-Associated Stroma in Ovarian Cancer. Genome Med (2016) 8(1):1–31. doi: 10.1186/s13073-016-0319-7
140. Levi-Galibov O, Lavon H, Wassermann-Dozorets R, Pevsner-Fischer M, Mayer S, Wershof E, et al. Heat Shock Factor 1-Dependent Extracellular Matrix Remodeling Mediates the Transition From Chronic Intestinal Inflammation to Colon Cancer. Nat Commun (2020) 11(1):1–19. doi: 10.1038/s41467-020-20054-x
141. Mirabelli P, Coppola L, Salvatore M. Cancer Cell Lines Are Useful Model Systems for Medical Research. Cancers (Basel) (2019) 11(8):1–18. doi: 10.3390/cancers11081098
142. Li H, Ning S, Ghandi M, Kryukov GV, Gopal S, Deik A, et al. The Landscape of Cancer Cell Line Metabolism. Nat Med (2019) 25(5):850–60. doi: 10.1038/s41591-019-0404-8
143. Hanniford D, Ulloa-Morales A, Karz A, Berzoti-Coelho MG, Moubarak RS, Sánchez-Sendra B, et al. Epigenetic Silencing of CDR1as Drives IGF2BP3-Mediated Melanoma Invasion and Metastasis. Cancer Cell (2020) 37(1):55–70.e15. doi: 10.1016/j.ccell.2019.12.007
144. Birkmeyer NJO, Goodney PP, Stukel TA, Hillner BE, Birkmeyer JD. Do Cancer Centers Designated by the National Cancer Institute Have Better Surgical Outcomes? Cancer (2005) 103(3):435–41. doi: 10.1002/cncr.20785
145. Deng M, Brägelmann J, Kryukov I, Saraiva-Agostinho N, Perner S. FirebrowseR: An R Client to the Broad Institute’s Firehose Pipeline. Database (2017) 2017:1–6. doi: 10.1093/database/baw160
146. Holroyd N, Sanchez-Flores A. Producing Parasitic Helminth Reference and Draft Genomes at the Wellcome Trust Sanger Institute. Parasit Immunol (2012) 34(2–3):100–7. doi: 10.1111/j.1365-3024.2011.01311.x
147. O’Sullivan RL, Keuthen NJ, Hayday CF, Ricciardi JN, Buttolph ML, Jenike MA, et al. The Massachusetts General Hospital (MGH) Hairpulling Scale: 2. Reliability and Validity. Psychother Psychosom (1995) 64(3–4):146–8. doi: 10.1159/000289004
148. Prior F, Smith K, Sharma A, Kirby J, Tarbox L, Clark K, et al. The Public Cancer Radiology Imaging Collections of The Cancer Imaging Archive. Sci Data (2017) 4(1):1–7. doi: 10.1038/sdata.2017.124
149. Li H, Zhu Y, Burnside ES, Drukker K, Hoadley KA, Fan C, et al. MR Imaging Radiomics Signatures for Predicting the Risk of Breast Cancer Recurrence as Given by Research Versions of MammaPrint, Oncotype DX, and PAM50 Gene Assays. Radiology (2016) 281(2):382–91. doi: 10.1148/radiol.2016152110
150. Sun R, Limkin EJ, Vakalopoulou M, Dercle L, Champiat S, Han SR, et al. A Radiomics Approach to Assess Tumour-Infiltrating CD8 Cells and Response to Anti-PD-1 or Anti-PD-L1 Immunotherapy: An Imaging Biomarker, Retrospective Multicohort Study. Lancet Oncol (2018) 19(9):1180–91. doi: 10.1016/S1470-2045(18)30413-3
151. Khosravi P, Kazemi E, Imielinski M, Elemento O, Hajirasouliha I. Deep Convolutional Neural Networks Enable Discrimination of Heterogeneous Digital Pathology Images. EBioMedicine (2018) 27:317–28. doi: 10.1016/j.ebiom.2017.12.026
152. Nahed BV, Curry WT, Martuza RL. Single-Cell RNA-Seq Highlights Intratumoral Heterogeneity in Primary Glioblastoma. Sci (80- ) (2014) 344(6190):1396–401. doi: 10.1126/science.1254257
153. Guo S, Jiang X, Mao B, Li Q-X. The Design, Analysis and Application of Mouse Clinical Trials in Oncology Drug Development. BMC Cancer (2019) 19(1):1–14. doi: 10.1186/s12885-019-5907-7
154. Strickland JC, Lile JA, Rush CR, Stoops WW. Comparing Exponential and Exponentiated Models of Drug Demand in Cocaine Users. Exp Clin Psychopharmacol (2016) 24(6):447–55. doi: 10.1037/pha0000096
155. Wishart DS, Knox C, Guo AC, Shrivastava S, Hassanali M, Stothard P, et al. DrugBank: A Comprehensive Resource for In Silico Drug Discovery and Exploration. Nucleic Acids Res (2006) 34(suppl_1):D668–72. doi: 10.1093/nar/gkj067
156. Augustyn A, Borromeo M, Wang T, Fujimoto J, Shao C, Dospoy PD, et al. ASCL1 Is a Lineage Oncogene Providing Therapeutic Targets for High-Grade Neuroendocrine Lung Cancers. Proc Natl Acad Sci USA (2014) 111(41):14788–93. doi: 10.1073/pnas.1410419111
157. Han K, Jeng EE, Hess GT, Morgens DW, Li A, Bassik MC. Synergistic Drug Combinations for Cancer Identified in a CRISPR Screen for Pairwise Genetic Interactions. Nat Biotechnol (2017) 35(5):463–74. doi: 10.1038/nbt.3834
158. Fernández-Navarro P, López-Nieva P, Piñeiro-Yañez E, Carreño-Tarragona G, Martinez-López J, Sánchez Pérez R, et al. The Use of PanDrugs to Prioritize Anticancer Drug Treatments in a Case of T-ALL Based on Individual Genomic Data. BMC Cancer (2019) 19(1):1–10. doi: 10.1186/s12885-019-6209-9
159. Karimi MR, Karimi AH, Abolmaali S, Sadeghi M, Schmitz U. Prospects and Challenges of Cancer Systems Medicine: From Genes to Disease Networks. Brief Bioinform (2021), 1–31. doi: 10.1093/bib/bbab343
160. Powathil GG, Swat M, Chaplain MAJ. Systems Oncology: Towards Patient-Specific Treatment Regimes Informed by Multiscale Mathematical Modelling. Semin Cancer Biol (2015) 30:13–20. doi: 10.1016/j.semcancer.2014.02.003
161. Valladares-Ayerbes M, Haz-Conde M, Blanco-Calvo M. Systems Oncology: Toward the Clinical Application of Cancer Systems Biology. Futur Oncol (2015) 11(4):553–5. doi: 10.2217/fon.14.255
162. Anderson ARA, Rejniak KA, Gerlee P, Quaranta V. Microenvironment Driven Invasion: A Multiscale Multimodel Investigation. J Math Biol (2009) 58(4–5):579–624. doi: 10.1007/s00285-008-0210-2
163. Stolarska MA, Yangjin KIM, Othmer HG. Multi-Scale Models of Cell and Tissue Dynamics. Philos Trans R Soc A Math Phys Eng Sci (2009) 367(1902):3525–53. doi: 10.1098/rsta.2009.0095
164. Silva A, Anderson A, Gatenby R. A Multiscale Model of the Bone Marrow and Hematopoiesis. Math Biosci Eng (2011) 8(2):643–58. doi: 10.3934/mbe.2011.8.643
165. Chaudhary SU, Shin S, Won J, Cho K, Member S. Multiscale Modeling of Tumorigenesis Induced by Mitochondrial Incapacitation in Cell Death. IEEE Trans BioMed Eng (2011) 58(2):3028–32. doi: 10.1109/TBME.2011.2159713
166. Szabó A, Merks RMH. Blood Vessel Tortuosity Selects Against Evolution of Aggressive Tumor Cells in Confined Tissue Environments: A Modeling Approach. PloS Comput Biol (2017) 13(7):1–32. doi: 10.1371/journal.pcbi.1005635
167. Kumar S, Das A, Sen S. Multicompartment Cell-Based Modeling of Confined Migration: Regulation by Cell Intrinsic and Extrinsic Factors. Mol Biol Cell (2018) 29(13):1599–610. doi: 10.1091/mbc.E17-05-0313
168. Unni P, Seshaiyer P. Mathematical Modeling, Analysis, and Simulation of Tumor Dynamics With Drug Interventions. Comput and Math Methods Med (2019) 2019:1–14. doi: 10.1155/2019/4079298
169. Leppert M, Dobbs M, Scambler P, O’Connell P, Nakamura Y, Stauffer D, et al. The Gene for Familial Polyposis Coli Maps to the Long Arm of Chromosome 5. Science (1987) 238(4832):1411–3. doi: 10.1126/science.3479843
170. Mehl LE. A Mathematical Computer Simulation Model for the Development of Colonic Polyps and Colon Cancer. J Surg Oncol (1991) 47(4):243–52. doi: 10.1002/jso.2930470409
171. Stratmann AT, Fecher D, Wangorsch G, Göttlich C, Walles T, Walles H, et al. Establishment of a Human 3D Lung Cancer Model Based on a Biological Tissue Matrix Combined With a Boolean In Silico Model. Mol Oncol (2014) 8(2):351–65. doi: 10.1016/j.molonc.2013.11.009
172. Greenspan HP. Models for the Growth of a Solid Tumor by Diffusion. Stud Appl Math (1972) 51(4):317–40. doi: 10.1002/sapm1972514317
173. Greenspan HP. On the Growth and Stability of Cell Cultures and Solid Tumors. J Theor Biol (1976) 56(1):229–42. doi: 10.1016/S0022-5193(76)80054-9
174. Chaplain MAJ. Avascular Growth, Angiogenesis and Vascular Growth in Solid Tumours: The Mathematical Modelling of the Stages of Tumour Development. Math Comput Model (1996) 23(6):47–87. doi: 10.1016/0895-7177(96)00019-2
175. Shackney SE, Smith CA, Miller BW, Burholt DR, Murtha K, Pollice AA, et al. Model for the Genetic Evolution of Human Solid Tumors. Cancer Res (1989) 49(12):3344–54.
176. Gerlee P, Nelander S. The Impact of Phenotypic Switching on Glioblastoma Growth and Invasion. PloS Comput Biol (2012) 8(6):e1002556. doi: 10.1371/journal.pcbi.1002556
177. Setty Y. In-Silico Models of Stem Cell and Developmental Systems. Theor Biol Med Model (2014) 11:1–12. doi: 10.1186/1742-4682-11-1
178. Wang Z, Bordas V, Sagotsky J, Deisboeck TS. Identifying Therapeutic Targets in a Combined EGFR–Tgfβr Signalling Cascade Using a Multiscale Agent-Based Cancer Model. Math Med Biol J IMA (2012) 29(1):95–108. doi: 10.1093/imammb/dqq023
179. Wang Z, Bordas V, Deisboeck T. Identification of Critical Molecular Components in a Multiscale Cancer Model Based on the Integration of Monte Carlo, Resampling, and ANOVA. Front Physiol (2011) 2:35. doi: 10.3389/fphys.2011.00035
180. Margolin G, Petrykowska HM, Jameel N, Bell DW, Young AC, Elnitski L. Robust Detection of DNA Hypermethylation of ZNF154 as a Pan-Cancer Locus With In Silico Modeling for Blood-Based Diagnostic Development. J Mol Diagn (2016) 18(2):283–98. doi: 10.1016/j.jmoldx.2015.11.004
181. Jahangiri A, Amani J, Halabian R, Imani fooladi AA. In Silico Analyses of Staphylococcal Enterotoxin B as a DNA Vaccine for Cancer Therapy. Int J Pept Res Ther (2018) 24(1):131–42. doi: 10.1007/s10989-017-9595-3
182. Smith EAK, Henthorn NT, Warmenhoven JW, Ingram SP, Aitkenhead AH, Richardson JC, et al. In Silico Models of DNA Damage and Repair in Proton Treatment Planning: A Proof of Concept. Sci Rep (2019) 9(1):1–10. doi: 10.1038/s41598-019-56258-5
183. Huang E, Cheng SH, Dressman H, Pittman J, Tsou MH, Horng CF, et al. Gene Expression Predictors of Breast Cancer Outcomes. Lancet (2003) 361(9369):1590–6. doi: 10.1016/S0140-6736(03)13308-9
184. Cheng F, Liu C, Lin CC, Zhao J, Jia P, Li WH, et al. A Gene Gravity Model for the Evolution of Cancer Genomes: A Study of 3,000 Cancer Genomes Across 9 Cancer Types. PloS Comput Biol (2015) 11(9):1–25. doi: 10.1371/journal.pcbi.1004497
185. Chen J-S, Hung W-S, Chan H-H, Tsai S-J, Sun HS. In Silico Identification of Oncogenic Potential of Fyn-Related Kinase in Hepatocellular Carcinoma. Bioinformatics (2013) 29(4):420–7. doi: 10.1093/bioinformatics/bts715
186. Agren R, Mardinoglu A, Asplund A, Kampf C, Uhlen M, Nielsen J. Identification of Anticancer Drugs for Hepatocellular Carcinoma Through Personalized Genome-Scale Metabolic Modeling. Mol Syst Biol (2014) 10(3):721. doi: 10.1002/msb.145122
187. Béal J, Montagud A, Traynard P, Barillot E, Calzone L. Personalization of Logical Models With Multi-Omics Data Allows Clinical Stratification of Patients. Front Physiol (2019) 9:1965. doi: 10.3389/fphys.2018.01965
188. Béal J, Pantolini L, Noël V, Barillot E, Calzone L. Personalized Logical Models to Investigate Cancer Response to BRAF Treatments in Melanomas and Colorectal Cancers. PloS Comput Biol (2021) 17(1):e1007900. doi: 10.1371/journal.pcbi.1007900
189. Rodriguez J, Ren G, Day CR, Zhao K, Chow CC, Larson DR. Intrinsic Dynamics of a Human Gene Reveal the Basis of Expression Heterogeneity. Cell (2019) 176(1–2):213–226.e18. doi: 10.1016/j.cell.2018.11.026
190. Baloria U, Akhoon BA, Gupta SK, Sharma S, Verma V. In Silico Proteomic Characterization of Human Epidermal Growth Factor Receptor 2 (HER-2) for the Mapping of High Affinity Antigenic Determinants Against Breast Cancer. Amino Acids (2012) 42(4):1349–60. doi: 10.1007/s00726-010-0830-x
191. Akhoon BA, Slathia PS, Sharma P, Gupta SK, Verma V. In Silico Identification of Novel Protective VSG Antigens Expressed by Trypanosoma Brucei and an Effort for Designing a Highly Immunogenic DNA Vaccine Using IL-12 as Adjuvant. Microb Pathog (2011) 51(1–2):77–87. doi: 10.1016/j.micpath.2011.01.011
192. Fang J, Wu Z, Cai C, Wang Q, Tang Y, Cheng F. Quantitative and Systems Pharmacology. 1. In Silico Prediction of Drug–Target Interactions of Natural Products Enables New Targeted Cancer Therapy. J Chem Inf Model (2017) 57(11):2657–71. doi: 10.1021/acs.jcim.7b00216
193. Azevedo R, Silva AMN, Reis CA, Santos LL, Ferreira JA. In Silico Approaches for Unveiling Novel Glycobiomarkers in Cancer. J Proteomics (2018) 171:95–106. doi: 10.1016/j.jprot.2017.08.004
194. Lee SWL, Seager RJ, Litvak F, Spill F, Sieow JL, Leong PH, et al. Integrated In Silico and 3D In Vitro Model of Macrophage Migration in Response to Physical and Chemical Factors in the Tumor Microenvironment. Integr Biol (2020) 12(4):90–108. doi: 10.1093/intbio/zyaa007
195. Ma H, Sorokin A, Mazein A, Selkov A, Selkov E, Demin O, et al. The Edinburgh Human Metabolic Network Reconstruction and its Functional Analysis. Mol Syst Biol (2007) 3(1):135. doi: 10.1038/msb4100177
196. Folger O, Jerby L, Frezza C, Gottlieb E, Ruppin E, Shlomi T. Predicting Selective Drug Targets in Cancer Through Metabolic Networks. Mol Syst Biol (2011) 7(1):501. doi: 10.1038/msb.2011.35
197. Aurich MK, Paglia G, Rolfsson Ó, Hrafnsdóttir S, Magnúsdóttir M, Stefaniak MM, et al. Prediction of Intracellular Metabolic States From Extracellular Metabolomic Data. Metabolomics (2015) 11(3):603–19. doi: 10.1007/s11306-014-0721-3
198. Yurkovich JT, Yang L, Palsson BO. Biomarkers Are Used to Predict Quantitative Metabolite Concentration Profiles in Human Red Blood Cells. PloS Comput Biol (2017) 13(3):e1005424. doi: 10.1371/journal.pcbi.1005424
199. Alakwaa FM, Chaudhary K, Garmire LX. Deep Learning Accurately Predicts Estrogen Receptor Status in Breast Cancer Metabolomics Data. J Proteome Res (2018) 17(1):337–47. doi: 10.1021/acs.jproteome.7b00595
200. Azadi M, Alemi F, Sadeghi S, Mohammadi M, Rahimi NA, Mirzaie S, et al. An Integrative In Silico Mathematical Modelling Study of the Anti-Cancer Effect of Clove Extract (Syzygium Aromaticum) Combined With In Vitro Metabolomics Study Using 1HNMR Spectroscopy. Iran J Biotechnol (2020) 18(3):45–54. doi: 10.30498/ijb.2020.141102.2336
201. Balding D, McElwain DL. A Mathematical Model of Tumour-Induced Capillary Growth. J Theor Biol (1985) 114(1):53–73. doi: 10.1016/s0022-5193(85)80255-1
202. Anderson ARA. A Hybrid Mathematical Model of Solid Tumour Invasion: The Importance of Cell Adhesion. Math Med Biol A J IMA (2005) 22(2):163–86. doi: 10.1093/imammb/dqi005
203. McDougall SR, Anderson ARA, Chaplain MAJ, Sherratt JA. Mathematical Modelling of Flow Through Vascular Networks: Implications for Tumour-Induced Angiogenesis and Chemotherapy Strategies. Bull Math Biol (2002) 64(4):673–702. doi: 10.1006/bulm.2002.0293
204. Anderson ARA, Chaplain MAJ. Continuous and Discrete Mathematical Models of Tumor-Induced Angiogenesis. Bull Math Biol (1998) 60(5):857–99. doi: 10.1006/bulm.1998.0042
205. Nakazawa K, Kurishima K, Tamura T, Kagohashi K, Ishikawa H, Satoh H, et al. Specific Organ Metastases and Survival in Small Cell Lung Cancer. Oncol Lett (2012) 4(4):617–20. doi: 10.3892/ol.2012.792
206. Cortesi M, Liverani C, Mercatali L, Ibrahim T, Giordano E. An In-Silico Study of Cancer Cell Survival and Spatial Distribution Within a 3D Microenvironment. Sci Rep (2020) 10(1):1–14. doi: 10.1038/s41598-020-69862-7
207. Aubert M, Badoual M, Christov C, Grammaticos B. A Model for Glioma Cell Migration on Collagen and Astrocytes. J R Soc Interface (2008) 5(18):75–83. doi: 10.1098/rsif.2007.1070
208. Harris LA, Beik S, Ozawa PMM, Jimenez L, Weaver AM. Modeling Heterogeneous Tumor Growth Dynamics and Cell–Cell Interactions at Single-Cell and Cell-Population Resolution. Curr Opin Syst Biol (2019) 17:24–34. doi: 10.1016/j.coisb.2019.09.005
209. Swanson KR, Alvord ECJ, Murray JD. A Quantitative Model for Differential Motility of Gliomas in Grey and White Matter. Cell Prolif (2000) 33(5):317–29. doi: 10.1046/j.1365-2184.2000.00177.x
210. Zhang L, Athale CA, Deisboeck TS. Development of a Three-Dimensional Multiscale Agent-Based Tumor Model: Simulating Gene-Protein Interaction Profiles, Cell Phenotypes and Multicellular Patterns in Brain Cancer. J Theor Biol (2007) 244(1):96–107. doi: 10.1016/j.jtbi.2006.06.034
211. Perfahl H, Byrne HM, Chen T, Estrella V, Alarcón T, Lapin A, et al. Multiscale Modelling of Vascular Tumour Growth in 3D: The Roles of Domain Size and Boundary Conditions. PloS One (2011) 6(4):e14790. doi: 10.1371/journal.pone.0014790
212. Anderson ARA, Rejniak KA, Gerlee P, Quaranta V. Modelling of Cancer Growth, Evolution and Invasion: Bridging Scales and Models. Math Model Nat Phenom (2007) 2(3):1–29. doi: 10.1051/mmnp:2007001
213. Chaudhary SU, Sormark JE, Won J-K, Shin SY, Cho K-H. Multi-Scale Modeling and Game-Theoretic Analysis of Mitochondrial Process Elucidates the Hidden Mechanism of Warburg Effect in Tumorigenesis. In: 11th International Conference on Systems Biology (Icsb2010) (2010).
214. Vavourakis V, Wijeratne PA, Shipley R, Loizidou M, Stylianopoulos T, Hawkes DJ. A Validated Multiscale in-Silico Model for Mechano-Sensitive Tumour Angiogenesis and Growth. PloS Comput Biol (2017) 13(1):e1005259. doi: 10.1371/journal.pcbi.1005259
215. Norton K-A, Wallace T, Pandey NB, Popel AS. An Agent-Based Model of Triple-Negative Breast Cancer: The Interplay Between Chemokine Receptor CCR5 Expression, Cancer Stem Cells, and Hypoxia. BMC Syst Biol (2017) 11(1):1–15. doi: 10.1186/s12918-017-0445-x
216. Norton KA, Gong C, Jamalian S, Popel AS. Multiscale Agent-Based and Hybrid Modeling of the Tumor Immune Microenvironment. Processes (2019) 7(1):1–23. doi: 10.3390/pr7010037
217. Karolak A, Poonja S, Rejniak KA. Morphophenotypic Classification of Tumor Organoids as an Indicator of Drug Exposure and Penetration Potential. PloS Comput Biol (2019) 15(7):e1007214. doi: 10.1371/journal.pcbi.1007214
218. Kenny PA, Lee GY, Myers CA, Neve RM, Semeiks JR, Spellman PT, et al. The Morphologies of Breast Cancer Cell Lines in Three-Dimensional Assays Correlate With Their Profiles of Gene Expression. Mol Oncol (2007) 1(1):84–96. doi: 10.1016/j.molonc.2007.02.004
219. Berrouet C, Dorilas N, Rejniak KA, Tuncer N. Comparison of Drug Inhibitory Effects (IC50) in Monolayer and Spheroid Cultures. Bull Math Biol (2020) 82(6):68. doi: 10.1007/s11538-020-00746-7
220. Stamatakos GS, Graf N, Radhakrishnan R. Multiscale Cancer Modeling and In Silico Oncology: Emerging Computational Frontiers in Basic and Translational Cancer Research. J Bioeng BioMed Sci (2013) 3(2):1–6. doi: 10.4172/2155-9538.1000e114
221. Nagaraj SH. Reverter A. A Boolean-Based Systems Biology Approach to Predict Novel Genes Associated With Cancer: Application to Colorectal Cancer. BMC Syst Biol (2011) 5(1):1–15. doi: 10.1186/1752-0509-5-35
222. Fumiã HF, Martins ML. Boolean Network Model for Cancer Pathways: Predicting Carcinogenesis and Targeted Therapy Outcomes. PloS One (2013) 8(7):e69008. doi: 10.1371/journal.pone.0069008
223. Yang L, Meng Y, Bao C, Liu W, Ma C, Li A, et al. Robustness and Backbone Motif of a Cancer Network Regulated by miR-17-92 Cluster During the G1/S Transition. PloS One (2013) 8(3):e57009. doi: 10.1371/journal.pone.0057009
224. der Heyde SV, Bender C, Henjes F, Sonntag J, Korf U, Beißbarth T. Boolean ErbB Network Reconstructions and Perturbation Simulations Reveal Individual Drug Response in Different Breast Cancer Cell Lines. BMC Syst Biol (2014) 8(1):1–22. doi: 10.1186/1752-0509-8-75
225. Gondal MN, Butt RN, Shah OS, Sultan MU, Mustafa G, Nasir Z, et al. A Personalized Therapeutics Approach Using an In Silico Drosophila Patient Model Reveals Optimal Chemo- and Targeted Therapy Combinations for Colorectal Cancer. Front Oncol (2021) 11:692592. doi: 10.3389/fonc.2021.692592
226. Klamt S, Saez-Rodriguez J, Gilles ED. Structural and Functional Analysis of Cellular Networks With CellNetAnalyzer. BMC Syst Biol (2007) 1(1):1–3. doi: 10.1186/1752-0509-1-2
227. Kachalo S, Zhang R, Sontag E, Albert R, DasGupta B. NET-SYNTHESIS: A Software for Synthesis, Inference and Simplification of Signal Transduction Networks. Bioinformatics (2008) 24(2):293–5. doi: 10.1093/bioinformatics/btm571
228. Mendes P. GEPASI: A Software Package for Modelling the Dynamics, Steady States and Control of Biochemical and Other Systems. Bioinformatics (1993) 9(5):563–71. doi: 10.1093/bioinformatics/9.5.563
229. Starruß J, De Back W, Brusch L, Deutsch A. Morpheus: A User-Friendly Modeling Environment for Multiscale and Multicellular Systems Biology. Bioinformatics (2014) 30(9):1331–2. doi: 10.1093/bioinformatics/btt772
230. Tomita M, Hashimoto K, Takahashi K, Shimizu TS, Matsuzaki Y, Miyoshi F, et al. E-CELL: Software Environment for Whole-Cell Simulation. Bioinformatics (1999) 15(1):72–84. doi: 10.1093/bioinformatics/15.1.72
231. Eissing T. A Computational Systems Biology Software Platform for Multiscale Modeling and Simulation: Integrating Whole-Body Physiology, Disease Biology, and Molecular Reaction Networks. Front Physiol (2011) 2:4. doi: 10.3389/fphys.2011.00004
232. Kauffman SA. Metabolic Stability and Epigenesis in Randomly Constructed Genetic Nets. J Theor Biol (1969) 22(3):437–67. doi: 10.1016/0022-5193(69)90015-0
233. Glass L, Kauffman SA. The Logical Analysis of Continuous, non-Linear Biochemical Control Networks. J Theor Biol (1973) 39(1):103–29. doi: 10.1016/0022-5193(73)90208-7
234. Wegner K, Knabe J, Robinson M, Egri-Nagy A, Schilstra M. The Netbuilder’ Project: Development of a Tool for Constructing, Simulating, Evolving, and Analysing Complex Regulatory Networks. BMC Syst Biol (2007) 1(1):1–2. doi: 10.1186/1752-0509-1-S1-P72
235. Albert I, Thakar J, Li S, Zhang R, Albert R. Boolean Network Simulations for Life Scientists. Source Code Biol Med (2008) 3(1):1–8. doi: 10.1186/1751-0473-3-16
236. Zheng J, Zhang D, Przytycki PF, Zielinski R, Capala J, Przytycka TM. Simboolnet-A Cytoscape Plugin for Dynamic Simulation of Signaling Networks. Bioinformatics (2009) 26(1):141–2. doi: 10.1093/bioinformatics/btp617
237. Klamt S, Stelling J, Ginkel M, Gilles ED. FluxAnalyzer: Exploring Structure, Pathways, and Flux Distributions in Metabolic Networks on Interactive Flux Maps. Bioinformatics (2003) 19(2):261–9. doi: 10.1093/bioinformatics/19.2.261
238. Tian K, Rajendran R, Doddananjaiah M, Krstic-Demonacos M, Schwartz J-M. Dynamics of DNA Damage Induced Pathways to Cancer. PloS One (2013) 8(9):e72303. doi: 10.1371/journal.pone.0072303
239. Hetmanski JHR, Zindy E, Schwartz J-M. Caswell PT. A MAPK-Driven Feedback Loop Suppresses Rac Activity to Promote RhoA-Driven Cancer Cell Invasion. PloS Comput Biol (2016) 12(5):e1004909. doi: 10.1371/journal.pcbi.1004909
240. Saadatpour A, Wang R-S, Liao A, Liu X, Loughran TP, Albert I, et al. Dynamical and Structural Analysis of a T Cell Survival Network Identifies Novel Candidate Therapeutic Targets for Large Granular Lymphocyte Leukemia. PloS Comput Biol (2011) 7(11):e1002267. doi: 10.1371/journal.pcbi.1002267
241. Steinway SN, Zañudo JGT, Ding W, Rountree CB, Feith DJ, Loughran TPJ, et al. Network Modeling of Tgfβ Signaling in Hepatocellular Carcinoma Epithelial-to-Mesenchymal Transition Reveals Joint Sonic Hedgehog and Wnt Pathway Activation. Cancer Res (2014) 74(21):5963–77. doi: 10.1158/0008-5472.CAN-14-0225
242. Naldi A, Berenguier D, Fauré A, Lopez F, Thieffry D, Chaouiya C. Logical Modelling of Regulatory Networks With GINsim 2.3. Biosystems (2009) 97(2):134–9. doi: 10.1016/j.biosystems.2009.04.008
243. Flobak Å, Baudot A, Remy E, Thommesen L, Thieffry D, Kuiper M, et al. Discovery of Drug Synergies in Gastric Cancer Cells Predicted by Logical Modeling. PloS Comput Biol (2015) 11(8):e1004426. doi: 10.1371/journal.pcbi.1004426
244. Remy E, Rebouissou S, Chaouiya C, Zinovyev A, Radvanyi F, Calzone L. A Modeling Approach to Explain Mutually Exclusive and Co-Occurring Genetic Alterations in Bladder Tumorigenesis. Cancer Res (2015) 75(19):4042–52. doi: 10.1158/0008-5472.CAN-15-0602
245. Stoll G, Caron B, Viara E, Dugourd A, Zinovyev A, Naldi A, et al. MaBoSS 2.0: An Environment for Stochastic Boolean Modeling. Bioinformatics (2017) 33(14):2226–8. doi: 10.1093/bioinformatics/btx123
246. Kondratova M, Barillot E, Zinovyev A, Calzone L. Modelling of Immune Checkpoint Network Explains Synergistic Effects of Combined Immune Checkpoint Inhibitor Therapy and the Impact of Cytokines in Patient Response. Cancers (Basel) (2020) 12(12):1–19. doi: 10.3390/cancers12123600
247. Müssel C, Hopfensitz M, Kestler HA. BoolNet—an R Package for Generation, Reconstruction and Analysis of Boolean Networks. Bioinformatics (2010) 26(10):1378–80. doi: 10.1093/bioinformatics/btq124
248. Steinway SN, Biggs MB, Loughran TP Jr., Papin JA, Albert R. Inference of Network Dynamics and Metabolic Interactions in the Gut Microbiome. PloS Comput Biol (2015) 11(6):e1004338. doi: 10.1371/journal.pcbi.1004338
249. Cohen DPA, Martignetti L, Robine S, Barillot E, Zinovyev A, Calzone L. Mathematical Modelling of Molecular Pathways Enabling Tumour Cell Invasion and Migration. PloS Comput Biol (2015) 11(11):e1004571. doi: 10.1371/journal.pcbi.1004571
250. Shah OS, Chaudhary MFA, Awan HA, Fatima F, Arshad Z, Amina B, et al. ATLANTIS - Attractor Landscape Analysis Toolbox for Cell Fate Discovery and Reprogramming. Sci Rep (2018) 8(1):3554. doi: 10.1038/s41598-018-22031-3
251. Han B, Wang J. Quantifying Robustness and Dissipation Cost of Yeast Cell Cycle Network: The Funneled Energy Landscape Perspectives. Biophys J (2007) 92(11):3755–63. doi: 10.1529/biophysj.106.094821
252. Choi M, Shi J, Jung SH, Chen X, Cho K-H. Attractor Landscape Analysis Reveals Feedback Loops in the P53 Network That Control the Cellular Response to DNA Damage. Sci Signal (2012) 5(251):1–14. doi: 10.1126/scisignal.2003363
253. Cho S-H, Park S-M, Lee H-S, Lee H-Y, Cho K-H. Attractor Landscape Analysis of Colorectal Tumorigenesis and its Reversion. BMC Syst Biol (2016) 10(1):1–13. doi: 10.1186/s12918-016-0341-9
254. Wittmann DM, Krumsiek J, Saez-Rodriguez J, Lauffenburger DA, Klamt S, Theis FJ. Transforming Boolean Models to Continuous Models: Methodology and Application to T-Cell Receptor Signaling. BMC Syst Biol (2009) 3(1):1–21. doi: 10.1186/1752-0509-3-98
255. Ricci G, De Maria F, Antonini G, Turella P, Bullo A, Stella L, et al. 7-Nitro-2,1,3-Benzoxadiazole Derivatives, a New Class of Suicide Inhibitors for Glutathione S-Transferases: Mechanism of Action of Potential Anticancer Drugs. J Biol Chem (2005) 280(28):26397–405. doi: 10.1074/jbc.M503295200
256. Marín-Hernández A, Gallardo-Pérez JC, Rodríguez-Enríquez S, Encalada R, Moreno-Sánchez R, Saavedra E. Modeling Cancer Glycolysis. Biochim Biophys Acta - Bioenerg (2011) 1807(6):755–67. doi: 10.1016/j.bbabio.2010.11.006
257. Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, et al. COPASI—A COmplex PAthway SImulator. Bioinformatics (2006) 22(24):3067–74. doi: 10.1093/bioinformatics/btl485
258. Orton RJ, Adriaens ME, Gormand A, Sturm OE, Kolch W, Gilbert DR. Computational Modelling of Cancerous Mutations in the EGFR/ERK Signalling Pathway. BMC Syst Biol (2009) 3(1):1–17. doi: 10.1186/1752-0509-3-100
259. Dalle Pezze P, Nelson G, Otten EG, Korolchuk VI, Kirkwood TBL, von Zglinicki T, et al. Dynamic Modelling of Pathways to Cellular Senescence Reveals Strategies for Targeted Interventions. PloS Comput Biol (2014) 10(8):e1003728. doi: 10.1371/journal.pcbi.1003728
260. Funahashi A, Morohashi M, Kitano H, Tanimura N. CellDesigner: A Process Diagram Editor for Gene-Regulatory and Biochemical Networks. Biosilico (2003) 1(5):159–62. doi: 10.1016/S1478-5382(03)02370-9
261. Sauro HM, Hucka M, Finney A, Wellock C, Bolouri H, Doyle J, et al. Next Generation Simulation Tools: The Systems Biology Workbench and BioSPICE Integration. OMICS (2003) 7(4):355–72. doi: 10.1089/153623103322637670
262. Calzone L, Gelay A, Zinovyev A, Radvanyi F, Barillot E. A Comprehensive Modular Map of Molecular Interactions in RB/E2F Pathway. Mol Syst Biol (2008) 4(1):0174. doi: 10.1038/msb.2008.7
263. Grieco L, Calzone L, Bernard-Pierrot I, Radvanyi F, Kahn-Perlès B, Thieffry D. Integrative Modelling of the Influence of MAPK Network on Cancer Cell Fate Decision. PloS Comput Biol (2013) 9(10):e1003286. doi: 10.1371/journal.pcbi.1003286
264. Ghaffarizadeh A, Friedman SH, MacKlin P. BioFVM: An Efficient, Parallelized Diffusive Transport Solver for 3-D Biological Simulations. Bioinformatics (2016) 32(8):1256–8. doi: 10.1093/bioinformatics/btv730
265. Felix TF, Lopez Lapa RM, de Carvalho M, Bertoni N, Tokar T, Oliveira RA, et al. MicroRNA Modulated Networks of Adaptive and Innate Immune Response in Pancreatic Ductal Adenocarcinoma. PloS One (2019) 14(5):e0217421. doi: 10.1371/journal.pone.0217421
266. Meyer K, Ostrenko O, Bourantas G, Morales-Navarrete H, Porat-Shliom N, Segovia-Miranda F, et al. A Predictive 3D Multi-Scale Model of Biliary Fluid Dynamics in the Liver Lobule. Cell Syst (2017) 4(3):277–290.e9. doi: 10.1016/j.cels.2017.02.008
267. Ozik J, Collier N, Heiland R, An G, Macklin P. Learning-Accelerated Discovery of Immune-Tumour Interactions. Mol Syst Des Eng (2019) 4(4):747–60. doi: 10.1039/c9me00036d
268. Ozik J, Collier N, Wozniak JM, Macal C, Cockrell C, Friedman SH, et al. High-Throughput Cancer Hypothesis Testing With an Integrated PhysiCell-EMEWS Workflow. BMC Bioinf (2018) 19(18):81–97. doi: 10.1186/s12859-018-2510-x
269. Wang Y, Brodin E, Nishii K, Frieboes HB, Mumenthaler SM, Sparks JL, et al. Impact of Tumor-Parenchyma Biomechanics on Liver Metastatic Progression: A Multi-Model Approach. Sci Rep (2021) 11(1):1–20. doi: 10.1038/s41598-020-78780-7
270. Cortesi M, Liverani C, Mercatali L, Ibrahim T, Giordano E. Development and Validation of an In-Silico Tool for the Study of Therapeutic Agents in 3D Cell Cultures. Comput Biol Med (2021) 130:1–9. doi: 10.1016/j.compbiomed.2021.104211
271. Edwards JS, Ibarra RU, Palsson BO. In Silico Predictions of Escherichia Coli Metabolic Capabilities Are Consistent With Experimental Data. Nat Biotechnol (2001) 19(2):125–30. doi: 10.1038/84379
272. Orton RJ, Sturm OE, Vyshemirsky V, Calder M, Gilbert DR, Kolch W. Computational Modelling of the Receptor-Tyrosine-Kinase-Activated MAPK Pathway. Biochem J (2005) 392(2):249–61. doi: 10.1042/BJ20050908
273. Loew LM, Schaff JC. The Virtual Cell: A Software Environment for Computational Cell Biology. Trends Biotechnol (2001) 19(10):401–6. doi: 10.1016/S0167-7799(01)01740-1
274. Neves SR, Tsokas P, Sarkar A, Grace EA, Rangamani P, Taubenfeld SM, et al. Cell Shape and Negative Links in Regulatory Motifs Together Control Spatial Information Flow in Signaling Networks. Cell (2008) 133(4):666–80. doi: 10.1016/j.cell.2008.04.025
275. Calebiro D, Nikolaev VO, Gagliani MC, de Filippis T, Dees C, Tacchetti C, et al. Persistent cAMP-Signals Triggered by Internalized G-Protein–Coupled Receptors. PloS Biol (2009) 7(8):e1000172. doi: 10.1371/journal.pbio.1000172
276. Chandran D, Bergmann FT, Sauro HM. TinkerCell: Modular CAD Tool for Synthetic Biology. J Biol Eng (2009) 3(1):1–17. doi: 10.1186/1754-1611-3-19
277. Renicke C, Schuster D, Usherenko S, Essen L-O, Taxis C. A LOV2 Domain-Based Optogenetic Tool to Control Protein Degradation and Cellular Function. Chem Biol (2013) 20(4):619–26. doi: 10.1016/j.chembiol.2013.03.005
278. Chandran D, Bergmann FT, Sauro HM. Computer-Aided Design of Biological Circuits Using TinkerCell. Bioeng Bugs (2010) 1(4):276–83. doi: 10.4161/bbug.1.4.12506
279. Izaguirre JA, Chaturvedi R, Huang C, Cickovski T, Coffland J, Thomas G, et al. CompuCell, a Multi-Model Framework for Simulation of Morphogenesis. Bioinformatics (2004) 20(7):1129–37. doi: 10.1093/bioinformatics/bth050
280. Mirams GR, Arthurs CJ, Bernabeu MO, Bordas R, Cooper J, Corrias A, et al. Chaste: An Open Source C++ Library for Computational Physiology and Biology. PloS Comput Biol (2013) 9(3):1–8. doi: 10.1371/journal.pcbi.1002970
281. Chaudhary SU, Shin S-Y, Lee D, Song J-H, Cho K-H. ELECANS–an Integrated Model Development Environment for Multiscale Cancer Systems Biology. Bioinformatics (2013) 29(7):957–9. doi: 10.1093/bioinformatics/btt063
282. North MJ, Collier NT, Ozik J, Tatara ER, Macal CM, Bragen M, et al. Complex Adaptive Systems Modeling With Repast Simphony. Complex Adapt Syst Model (2013) 1(1):1–26. doi: 10.1186/2194-3206-1-3
283. Folcik VA, An GC, Orosz CG. The Basic Immune Simulator: An Agent-Based Model to Study the Interactions Between Innate and Adaptive Immunity. Theor Biol Med Model (2007) 4(1):1–18. doi: 10.1186/1742-4682-4-39
284. Mehdizadeh H, Sumo S, Bayrak ES, Brey EM, Cinar A. Three-Dimensional Modeling of Angiogenesis in Porous Biomaterial Scaffolds. Biomaterials (2013) 34(12):2875–87. doi: 10.1016/j.biomaterials.2012.12.047
285. Mahoney AW, Podgorski GJ, Flann NS. Multiobjective Optimization Based-Approach for Discovering Novel Cancer Therapies. IEEE/ACM Trans Comput Biol Bioinforma (2012) 9(1):169–84. doi: 10.1109/TCBB.2010.39
286. Cho K-H, Chaudhary SU, Lee D Biosimulation Method and Computing Device With High Expandability. Korea Patent No. 10-2013-0033839. Korea: Korea Patent Office (2013).
287. Cho KH, Chaudhary SU, Lee D Biosimulation Method and Computing Device Using Compiler Embedded There Within. Korea Patent Office; 10-2013–0033837. South Korea: Korea Patent Office (2013).
288. Ghaffarizadeh A, Heiland R, Friedman SH, Mumenthaler SM, Macklin P. PhysiCell: An Open Source Physics-Based Cell Simulator for 3-D Multicellular Systems. PloS Comput Biol (2018) 14(2):e1005991. doi: 10.1371/journal.pcbi.1005991
289. Letort G, Montagud A, Stoll G, Heiland R, Barillot E, MacKlin P, et al. PhysiBoSS: A Multi-Scale Agent-Based Modelling Framework Integrating Physical Dimension and Cell Signalling. Bioinformatics (2019) 35(7):1188–96. doi: 10.1093/bioinformatics/bty766
290. Colin A, Letort G, Razin N, Almonacid M, Ahmed W, Betz T, et al. Active Diffusion in Oocytes Nonspecifically Centers Large Objects During Prophase I and Meiosis I. J Cell Biol (2020) 219(3):1–21. doi: 10.1083/jcb.201908195
291. Getz M, Wang Y, An G, Becker A, Cockrell C, Collier N. Rapid Community-Driven Development of a SARS- CoV-2 Tissue Simulator. bioRxiv (2020), 1–55. doi: 10.1101/2020.04.02.019075
292. Gondal MN, Sultan MU, Arif A, Rehman A, Awan HA, Arshad Z, et al. TISON: A Next-Generation Multi-Scale Modeling Theatre for In Silico Systems Oncology. bioRxiv (2021). doi: 10.1101/2021.05.04.442539
293. Martin L, Anguita A, Graf N, Tsiknakis M, Brochhausen M, Rüping S, et al. ACGT: Advancing Clinico-Genomic Trials on Cancer–Four Years of Experience. Stud Heal Technol Inf (2011) 169:734–8. doi: 10.3233/978-1-60750-806-9-734
294. Marias K, Sakkalis V, Roniotis A, Farmaki C, Stamatakos G, Dionysiou D, et al. Clinically Oriented Translational Cancer Multilevel Modeling: The ContraCancrum Project. World Congress Medical Physics Biomed Eng (2009) 2124–7. doi: 10.1007/978-3-642-03882-2_564
295. Rossi S, Christ-Neumann M, Rüping S, Buffa F, Wegener D, McVie G, et al. P-Medicine: From Data Sharing and Integration via VPH Models to Personalized Medicine. Ecancermedicalscience (2011) 5:1–6. doi: 10.3332/ecancer.2011.218
296. Johnson D, McKeever S, Stamatakos G, Dionysiou D, Graf N, Sakkalis V, et al. Dealing With Diversity in Computational Cancer Modeling. Cancer Inform (2013) 12:1–10. doi: 10.4137/CIN.S11583
297. Stamatakos G, Graf N. Computational Horizons in Cancer (CHIC). Clin Ther (2017) 39(8):e107–8. doi: 10.1016/j.clinthera.2017.05.333
298. Durbecq V, Paesmans M, Cardoso F, Desmedt C, Di Leo A, Chan S, et al. Topoisomerase-Iiα Expression as a Predictive Marker in a Population of Advanced Breast Cancer Patients Randomly Treated Either With Single-Agent Doxorubicin or Single-Agent Docetaxel. Mol Cancer Ther (2004) 3(10):1207–14.
299. Graf N, Hoppe A, Georgiadi E, Belleman R, Desmedt C, Dionysiou D, et al. ‘In Silico’ Oncology for Clinical Decision Making in the Context of Nephroblastoma. Klin Paediatr (2009) 221(03):141–9. doi: 10.1055/s-0029-1216368
300. Zasada SJ, Wang T, Haidar A, Liu E, Graf N, Clapworthy G, et al. IMENSE: An E-Infrastructure Environment for Patient Specific Multiscale Data Integration, Modelling and Clinical Treatment. J Comput Sci (2012) 3(5):314–27. doi: 10.1016/j.jocs.2011.07.001
301. Marias K, Dionysiou D, Sakkalis V, Graf N, Bohle RM, Coveney PV, et al. Clinically Driven Design of Multi-Scale Cancer Models: The ContraCancrum Project Paradigm. Interface Focus (2011) 1(3):450–61. doi: 10.1098/rsfs.2010.0037
302. Li XL, Oduola WO, Qian L, Dougherty ER. Integrating Multiscale Modeling With Drug Effects for Cancer Treatment. Cancer Inform (2015) 14(Suppl 5):21–31. doi: 10.4137/CIN.S30797
303. Stamatakos GS, Georgiadi EC, Graf N, Kolokotroni EA, Dionysiou DD. Exploiting Clinical Trial Data Drastically Narrows the Window of Possible Solutions to the Problem of Clinical Adaptation of a Multiscale Cancer Model. PloS One (2011) 6(3):e17594. doi: 10.1371/journal.pone.0017594
304. Wan S, Coveney PV. Rapid and Accurate Ranking of Binding Affinities of Epidermal Growth Factor Receptor Sequences With Selected Lung Cancer Drugs. J R Soc Interface (2011) 8(61):1114–27. doi: 10.1098/rsif.2010.0609
305. Roniotis A, Panourgias K, Ekaterinaris J, Marias K, Sakkalis V. Approximating the Diffusion–Reaction Equation for Developing Glioma Models for the ContraCancrum Project: A Showcase. In: 4th Int Adv Res Work Silico Oncol Cancer Investig (4th IARWISOCI) – ContraCancrum Work. p. 80–2.
306. Folarin AA, Stamatakos GS. Molecular Personalization of Cancer Treatment via a Multiscale Simulation Model of Tumor Response to Therapy. The Paradigm of Glioblastoma Treated With Temozolomide. In: 4th Int Adv Res Work Silico Oncol Cancer Investig (2010). p. 23–5.
307. Giatili SG, Uzunoglu NK, Stamatakos GS. An Explicit Boundary Condition Treatment of a Diffusion–Based Glioblastoma Tumor Growth Model. 4th Int Adv Res Work Silico Oncol Cancer Investig (4th IARWISOCI) – ContraCancrum Work (2010), 86–8.
308. Wang Z, Deisboeck TS. Application of ANOVA-Based Global Sensitivity Analysis to a Multiscale Cancer Model. In: 4th Int Adv Res Work Silico Oncol Cancer Investig (4th IARWISOCI) – ContraCancrum Work. p. 27–9.
309. Schera F, Weiler G, Neri E, Kiefer S, Graf N. The P-Medicine Portal—A Collaboration Platform for Research in Personalised Medicine. Ecancermedicalscience (2014) 8:1–18. doi: 10.3332/ecancer.2014.398
310. Georgiadi EC, Dionysiou DD, Graf N, Stamatakos GS. Towards In Silico Oncology: Adapting a Four Dimensional Nephroblastoma Treatment Model to a Clinical Trial Case Based on Multi-Method Sensitivity Analysis. Comput Biol Med (2012) 42(11):1064–78. doi: 10.1016/j.compbiomed.2012.08.008
311. Blazewicz M, Georgiadi EC, Pukacki J, Stamatakos GS. Development of the P-Medicine Oncosimulator as a Parallel Treatment Support System. In: 6th Int Adv Res Work Silico Oncol Cancer Investig - CHIC Proj Work. Athens, Greece: Proceedings of the 2014 6th International Advanced Research Workshop on In Silico Oncology and Cancer Investigation - The CHIC Project Workshop (IARWISOCI) (2014). p. 1–5. doi: 10.1109/IARWISOCI.2014.7034641
312. Argyri KD, Dionysiou DD, Misichroni FD, Stamatakos GS. Numerical Simulation of Vascular Tumour Growth Under Antiangiogenic Treatment: Addressing the Paradigm of Single-Agent Bevacizumab Therapy With the Use of Experimental Data. Biol Direct (2016) 11(1):1–31. doi: 10.1186/s13062-016-0114-9
313. Stamatakos G, Dionysiou D, Lunzer A, Belleman R, Kolokotroni E, Georgiadi E, et al. The Technologically Integrated Oncosimulator: Combining Multiscale Cancer Modeling With Information Technology in the In Silico Oncology Context. IEEE J BioMed Heal Informatics. Athens, Greece: The TUMOR Project Workshop (IARWISOCI) (2014) 18(3):840–54. doi: 10.1109/JBHI.2013.2284276
314. Ouzounoglou EN, Dionysiou DD, Stanulla M, Stamatakos GS. Towards Patient Personalization of an Acute Lymphoblastic Leukemia Model During the Oral Administration of Prednisone in Children: Initiating the ALL Oncosimulator. In: 5th Int Adv Res Work Silico Oncol Cancer Investig (2012). p. 1–4.
315. Kolokotroni E, Ouzounoglou E, Systems C, Stanulla M, Dionysiou DD, Systems C. In Silico Oncology: Developing and Clinically Adapting the Acute Lymphoblastic Leukemia (ALL) Oncosimulator by Exploiting Pathway Based Gene Expression Analysis in the Context of the ALL-BFM 2000 Clinical Study. Virtual Physiol Hum Conf (2015). doi: 10.13140/RG.2.1.1430.3123
316. Ouzounoglou E, Kolokotroni E, Stanulla M, Stamatakos GS. A Study on the Predictability of Acute Lymphoblastic Leukaemia Response to Treatment Using a Hybrid Oncosimulator. Interface Focus (2018) 8(1):1–14. doi: 10.1098/rsfs.2016.0163
317. Bucur A, van Leeuwen J, Christodoulou N, Sigdel K, Argyri K, Koumakis L, et al. Workflow-Driven Clinical Decision Support for Personalized Oncology. BMC Med Inform Decis Mak (2016) 16(2):151–62. doi: 10.1186/s12911-016-0314-3
318. Deisboeck TS, Zhang L, Martin S. Advancing Cancer Systems Biology: Introducing the Center for the Development of a Virtual Tumor, CViT. Cancer Inform (2007) 5:1–8. doi: 10.1177/117693510700500001
319. Hedley WJ, Nelson MR, Bellivant DP, Nielsen PF. A Short Introduction to CellML. Philos Trans R Soc London Ser A Math Phys Eng Sci (2001) 359(1783):1073–89. doi: 10.1098/rsta.2001.0817
320. Christie GR, Nielsen PMF, Blackett SA, Bradley CP, Hunter PJ. FieldML: Concepts and Implementation. Philos Trans R Soc A Math Phys Eng Sci (2009) 367(1895):1869–84. doi: 10.1098/rsta.2009.0025
321. Johnson D, McKeever S, Deisboeck TS, Wang Z. Connecting Digital Cancer Model Repositories With Markup: Introducing TumorML Version 1.0. ACM SIGBioinform Rec (2013) 3(3):5–11. doi: 10.1145/2544063.2544064
322. Wang Z, Zhang L, Sagotsky J, Deisboeck TS. Simulating non-Small Cell Lung Cancer With a Multiscale Agent-Based Model. Theor Biol Med Model (2007) 4(1):1–14. doi: 10.1186/1742-4682-4-50
323. Sakkalis V, Sfakianakis S, Tzamali E, Marias K, Stamatakos G, Misichroni F, et al. Web-Based Workflow Planning Platform Supporting the Design and Execution of Complex Multiscale Cancer Models. IEEE J BioMed Heal Informatics (2014) 18(3):824–31. doi: 10.1109/JBHI.2013.2297167
324. Athale C, Mansury Y, Deisboeck TS. Simulating the Impact of a Molecular ‘Decision-Process’ on Cellular Phenotype and Multicellular Patterns in Brain Tumors. J Theor Biol (2005) 233(4):469–81. doi: 10.1016/j.jtbi.2004.10.019
325. Shlomi T, Benyamini T, Gottlieb E, Sharan R, Ruppin E. Genome-Scale Metabolic Modeling Elucidates the Role of Proliferative Adaptation in Causing the Warburg Effect. PloS Comput Biol (2011) 7(3):e1002018. doi: 10.1371/journal.pcbi.1002018
326. Kolokotroni E, Dionysiou D, Veith C, Kim Y-J, Sabczynski J, Franz A, et al. In Silico Oncology: Quantification of the In Vivo Antitumor Efficacy of Cisplatin-Based Doublet Therapy in non-Small Cell Lung Cancer (NSCLC) Through a Multiscale Mechanistic Model. PloS Comput Biol (2016) 12(9):e1005093. doi: 10.1371/journal.pcbi.1005093
327. Antonopoulos M, Stamatakos G. In Silico Neuro-Oncology: Brownian Motion-Based Mathematical Treatment as a Potential Platform for Modeling the Infiltration of Glioma Cells Into Normal Brain Tissue. Cancer Inform (2015) 14:33–40. doi: 10.4137/CIN.S19341
328. Ouzounoglou E, Dionysiou D, Stamatakos GS. Differentiation Resistance Through Altered Retinoblastoma Protein Function in Acute Lymphoblastic Leukemia: In Silico Modeling of the Deregulations in the G1/S Restriction Point Pathway. BMC Syst Biol (2016) 10(1):1–22. doi: 10.1186/s12918-016-0264-5
329. Stamatakos GS, Giatili SG. A Numerical Handling of the Boundary Conditions Imposed by the Skull on an Inhomogeneous Diffusion-Reaction Model of Glioblastoma Invasion Into the Brain: Clinical Validation Aspects. Cancer Inform (2017) 16:1–16. doi: 10.1177/1176935116684824
330. Lee HJ, Palm J, Grimes SM, Ji HP. The Cancer Genome Atlas Clinical Explorer: A Web and Mobile Interface for Identifying Clinical-Genomic Driver Associations. Genome Med (2015) 7(1):1–14. doi: 10.1186/s13073-015-0226-3
331. Shaul YD, Yuan B, Thiru P, Nutter-Upham A, McCallum S, Lanzkron C, et al. MERAV: A Tool for Comparing Gene Expression Across Human Tissues and Cell Types. Nucleic Acids Res (2016) 44(D1):D560–6. doi: 10.1093/nar/gkv1337
332. Jensen MA, Ferretti V, Grossman RL, Staudt LM. The NCI Genomic Data Commons as an Engine for Precision Medicine. Blood (2017) 130(4):453–9. doi: 10.1182/blood-2017-03-735654
333. Li J, Lu Y, Akbani R, Ju Z, Roebuck PL, Liu W, et al. TCPA: A Resource for Cancer Functional Proteomics Data. Nat Methods (2013) 10(11):1046–7. doi: 10.1038/nmeth.2650
334. Pontén F, Jirström K, Uhlen M. The Human Protein Atlas—a Tool for Pathology. J Pathol (2008) 216(4):387–93. doi: 10.1002/path.2440
Keywords: multi-scale cancer modeling, personalized cancer therapeutics, cancer systems biology, data-driven oncology, in silico cancer systems oncology, predictive systems oncology
Citation: Gondal MN and Chaudhary SU (2021) Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics. Front. Oncol. 11:712505. doi: 10.3389/fonc.2021.712505
Received: 20 May 2021; Accepted: 26 October 2021;
Published: 24 November 2021.
Edited by:
George Mattheolabakis, University of Louisiana at Monroe, United StatesReviewed by:
Andrei Leitao, University of São Paulo, BrazilUlf Schmitz, James Cook University, Australia
Copyright © 2021 Gondal and Chaudhary. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Safee Ullah Chaudhary, safee.ullah.chaudhary@gmail.com