Intra-tumor and inter-patient heterogeneity are two challenges in developing mathematical models for precision medicine diagnostics. Here we review several techniques that can be used to aid the mathematical modeller in inferring and quantifying both sources of heterogeneity from patient data. These techniques include virtual populations, nonlinear mixed effects modeling, non-parametric estimation, Bayesian techniques, and machine learning. We create simulated virtual populations in this study and then apply the four remaining methods to these datasets to highlight the strengths and weak-nesses of each technique. We provide all code used in this review at https://github.com/jtnardin/Tumor-Heterogeneity/ so that this study may serve as a tutorial for the mathematical modelling community. This review article was a product of a Tumor Heterogeneity Working Group as part of the 2018–2019 Program on Statistical, Mathematical, and Computational Methods for Precision Medicine which took place at the Statistical and Applied Mathematical Sciences Institute.
Citation: Rebecca Everett, Kevin B. Flores, Nick Henscheid, John Lagergren, Kamila Larripa, Ding Li, John T. Nardini, Phuong T. T. Nguyen, E. Bruce Pitman, Erica M. Rutter. A tutorial review of mathematical techniques for quantifying tumor heterogeneity[J]. Mathematical Biosciences and Engineering, 2020, 17(4): 3660-3709. doi: 10.3934/mbe.2020207
[1] | Judit Faus-Garriga, Isabel Novoa, Andrés Ozaita . mTOR signaling in proteostasis and its relevance to autism spectrum disorders. AIMS Biophysics, 2017, 4(1): 63-89. doi: 10.3934/biophy.2017.1.63 |
[2] | Warin Rangubpit, Pornthep Sompornpisut, R.B. Pandey . Thermal-induced unfolding-refolding of a nucleocapsid COVN protein. AIMS Biophysics, 2021, 8(1): 103-110. doi: 10.3934/biophy.2021007 |
[3] | Jonathan Trauth, Johannes Scheffer, Sophia Hasenjäger, Christof Taxis . Strategies to investigate protein turnover with fluorescent protein reporters in eukaryotic organisms. AIMS Biophysics, 2020, 7(2): 90-118. doi: 10.3934/biophy.2020008 |
[4] | Alessandro Didonna, Federico Benetti . Post-translational modifications in neurodegeneration. AIMS Biophysics, 2016, 3(1): 27-49. doi: 10.3934/biophy.2016.1.27 |
[5] | Shan-Yang Lin . Salmon calcitonin: conformational changes and stabilizer effects. AIMS Biophysics, 2015, 2(4): 695-723. doi: 10.3934/biophy.2015.4.695 |
[6] | Jany Dandurand, Angela Ostuni, Maria Francesca Armentano, Maria Antonietta Crudele, Vincenza Dolce, Federica Marra, Valérie Samouillan, Faustino Bisaccia . Calorimetry and FTIR reveal the ability of URG7 protein to modify the aggregation state of both cell lysate and amylogenic α-synuclein. AIMS Biophysics, 2020, 7(3): 189-203. doi: 10.3934/biophy.2020015 |
[7] | Ateeq Al-Zahrani, Natasha Cant, Vassilis Kargas, Tracy Rimington, Luba Aleksandrov, John R. Riordan, Robert C. Ford . Structure of the cystic fibrosis transmembrane conductance regulator in the inward-facing conformation revealed by single particle electron microscopy. AIMS Biophysics, 2015, 2(2): 131-152. doi: 10.3934/biophy.2015.2.131 |
[8] | Davide Sala, Andrea Giachetti, Antonio Rosato . Molecular dynamics simulations of metalloproteins: A folding study of rubredoxin from Pyrococcus furiosus. AIMS Biophysics, 2018, 5(1): 77-96. doi: 10.3934/biophy.2018.1.77 |
[9] | Nicola Gaetano Gatta, Gaetano Cammarota, Vittorio Gentile . Possible roles of transglutaminases in molecular mechanisms responsible for human neurodegenerative diseases. AIMS Biophysics, 2016, 3(4): 529-545. doi: 10.3934/biophy.2016.4.529 |
[10] | Jay L. Brewster . Signaling hubs at ER/mitochondrial membrane associations. AIMS Biophysics, 2017, 4(2): 222-239. doi: 10.3934/biophy.2017.2.222 |
Intra-tumor and inter-patient heterogeneity are two challenges in developing mathematical models for precision medicine diagnostics. Here we review several techniques that can be used to aid the mathematical modeller in inferring and quantifying both sources of heterogeneity from patient data. These techniques include virtual populations, nonlinear mixed effects modeling, non-parametric estimation, Bayesian techniques, and machine learning. We create simulated virtual populations in this study and then apply the four remaining methods to these datasets to highlight the strengths and weak-nesses of each technique. We provide all code used in this review at https://github.com/jtnardin/Tumor-Heterogeneity/ so that this study may serve as a tutorial for the mathematical modelling community. This review article was a product of a Tumor Heterogeneity Working Group as part of the 2018–2019 Program on Statistical, Mathematical, and Computational Methods for Precision Medicine which took place at the Statistical and Applied Mathematical Sciences Institute.
Production of any protein requires that the information for the specific sequence of its amino acids is transcribed from DNA into an mRNA molecule. The mRNA will then be translated by the complex machinery of the ribosome and its co-factors into a polypeptide chain that must be folded, with assistance of chaperones and folding enzymes, in its unique architecture. First thoughts on the flow of genetic information were enounced by Francis Crick in 1956, when he exposed his “Ideas on protein synthesis” [1]. These were later published and partially revisited in a Nature article in 1970 defining what he pompously (and inappropriately, as himself later admitted) defined as “Central dogma of molecular biology” [2] (Figure 1).
Although the (mis) use of the word “dogma” generated some misunderstanding with colleagues, the principles he wanted to highlight were rather simple: once the information, intended as the sequence of amino acids in a protein is transferred from DNA into proteins, it cannot go back. In 1970, Crick proposed a scenario where the flow of information would generally occur from DNA to DNA, DNA to RNA and from RNA to protein (what he called “general transfer”). RNA replication and flow of information from RNA to DNA and from DNA directly to protein were defined as “special transfers” meaning that they could happen under particular conditions. “Unknown transfers” were the ones from protein to protein as well as the reverse flow of information from protein to RNA or DNA. At the time, there were no experimental evidences for the latter and thus Crick hypothesized that these transfers were improbable to exist [2].
As science moved on, the concepts highlighted in Crick’s “central dogma” have dramatically changed and experimental evidences have proven the existence of “unknown transfers” such as reverse transcription (RNA to DNA) by Temin and Baltimore [3,4] and the description of the prion-mediated inheritance [5].
In any case, and independently of what Crick meant to say at the time he posted the “Ideas on protein synthesis” [1] and later on [2], it remains that translation and attainment of the correct ternary and quaternary structure of polypeptides are crucial steps in the flow of genetic information. Proteins are the active molecules that carry out almost every function in a living organism, while DNA stores the information to produce each one of them.
Translation of genetic information (sequences of nucleotides) into a polypeptide chain (sequence of amino acids) is mediated by ribosomes. These are complex structures that act as RNA-based enzymes (ribozymes) to catalyze peptide bond formation. Eukaryotic 80S ribosomes are large (4.3 MDa for the human ribosome) and complex structures composed of a large 60S subunit, which contains 47 proteins and 3 ribosomal (r)RNAs, and a small 40S subunit containing 33 proteins and 1 rRNA molecule [6,7].
Translation can be divided in 3 steps: Initiation, elongation and termination (Figure 2).
i) Initiation: Amino acids reach the ribosome conjugated with a specific transfer (t)RNA molecule (amino-acyl-tRNAs) that possesses a complementary anticodon for the codon encoding the specific amino acid. In eukaryotes, at least eleven different initiation factors are required to properly initiate translation. These ensure that the methionyl-initiator tRNA (Met-tRNAiMet) is brought in the P site of the ribosomal complex to the initiator AUG codon of an engaged mRNA [8] (Figure 2, step 1).
ii) Elongation: Successive amino-acyl-tRNAs are recruited to the ribosome on the A-site (Figure 2, step 2). While the mRNA molecule is pulled through the ribosome, codon-matching amino-acyl-tRNAs are recruited one after the other in order to link the amino acids into a peptide chain. During elongation the tRNA is recruited to the A-site while an adjacent tRNA in the P-site is covalently bound to the growing amino acid chain (Figure 2, step 3).
At this point a peptidyl transferase contained in the large ribosomal subunit breaks the high energy bond between the tRNA in the P-site and its amino acid, and links it to the free amino acid group on the tRNA in the A-site (Figure 2, step 4). This reaction causes a conformational shift of the large ribosomal subunit relative to the mRNA (which is bound by the small ribosomal subunit) leading to the shift of the bound tRNAs from the P-to the E-site and from the A-to the P-site (Figure 2, step 5). Next, the small ribosomal subunit undergoes a conformational change and moves, together with the mRNA, exactly 3 nucleotides. This resets the ribosomal structure leaving an empty A-site that can now accommodate the next complementary amino-acyl-tRNA, which binding will lead to the release of the “empty” tRNA from the E-site (Figure 2, step 6) [9].
iii) Termination: Translation gets terminated when a stop codon (UAA, UAG or UGA) is encountered by the A-site of the ribosome. These codons stop translation by engaging so called release factors that bind the A-site and force the peptidyl-transferase to catalyze the addition of a water molecule to the C-terminal peptidyl-tRNA thereby releasing the carboxyl terminus from its tRNA and leading to the disassembly of the large and small subunits of the ribosome (Figure 2, step 7) [10].
Translation is a rapid process. The catalytic activity of the ribozyme elongates the polypeptide chain at a rate of 3-10 amino acids per second in eukaryotes [11,12] and up to 10-20 amino acids per second in E. coli [13,14,15].
Every step in the biosynthetic process of proteins from DNA replication to protein folding (including transcription, amino-acyl-tRNA synthesis and translation) is error prone. Following chapters will focus on the fidelity of the mentioned processes and on quality control and repair mechanisms that operate in order to minimize the insertion of errors in the flow of genetic information from DNA to the mature protein. We will also discuss how cells maintain proteostasis i.e., how they ensure production of the proteome in correct amount and appropriate quality.
DNA replication is performed with very high fidelity since the correct maintenance of the genetic material is fundamental for the survival not only of organisms but, most importantly, for species. DNA polymerases (the enzymes that catalyze the synthesis of new DNA strands) and their co-factors have evolved in order to minimize the insertion of wrong nucleotides. In fact, DNA replication machineries are estimated to insert a mismatched base every 104-105 nucleotides, but through the combined actions of proofreading and post-replication mismatch repair mechanisms, eukaryotic cells can replicate their genome with remarkable fidelity of less than 1 mutation per genome per cell division (error rate of 10-8-10-10) [16,17,18,19,20] (Figure 1). Although DNA replication error rate is exquisitely low, nucleotide mis-incorporation is fundamental for adaptation to changing environmental conditions and for driving evolution [20].
mRNA is the principal template molecule used for protein synthesis. Genes encoded in genomic DNA are transcribed, in the nucleus, by multi-component transcription machineries, which comprise RNA polymerases (RNA pol) as principal component, into a “transient” single stranded nucleotides chain (mRNA). This is transported into the cytosol where it can be engaged by the ribosomes (for reviews see [21,22,23]). The mRNA molecule has a relatively short half-life (around 9 h) and is re-synthesized constantly (typically 17 copies/gene are present) [24]. The continuous renewal of the mRNA pool implies that errors during transcription have less severe consequences than errors occurring during DNA replication and also fidelity of the process is less stringent (about 1 wrong nucleotide every 104-106 bases both in vitro and in vivo)[25,26,27] (Figure 1).
When it comes to translation of mRNA polynucleotides into amino acid chains, the error rate increases. Generally, the error rate for translation is one order of magnitude higher than that for transcription, laying between a wrong amino acid every 103-105 codons [28,29,30,31,32] (Figure 1). An evolutionary interpretation for the lower error rate in transcription than in translation is that an error in transcription would lead to many erroneous protein copies (on average 17 mRNA/gene code for 50, 000 proteins/gene [24], meaning ca. 2950 proteins/mRNA), whereas an error in translation affects only one polypeptide chain if the error is due to stochastic mis-incorporation (and not a defective machinery). Note that in addition to the mistranslation of mRNAs, the protein synthesis efficiency can also be hampered by the incorrect charging of the tRNAs themselves. Incorporation of the wrong amino acid on a given tRNA has been measured to occur with error rates of 10-4-10-5 [28] (Figure 2).
Protein biogenesis is not terminated upon translation. Most proteins must achieve a specific three dimensional structure, which is determined by its amino acid sequence. The final, or native, conformation of a polypeptide chain is generally the one that minimizes its free energy and can occur naturally in a test tube [33]. In the living cell, however, the polypeptide chain emerging from ribosomes has to undergo an active, chaperone-assisted process of folding (see 4.2), and in some cases of oligomerization, i.e., of pairing with other polypeptide chains, in order to achieve its functional conformation (tertiary and quaternary structures). Only then, proteins can be transported to their final destination. As one can imagine, the process of folding is all but error-free and needs strict quality assessment in order to permit the synthesis of a functional proteome.
The efficiency of polypeptide folding, as well as the rate, varies substantially for different gene products. Intuitively, short, single domain proteins achieve their native structure rapidly and efficiently, whereas long, multidomain and multi-subunit complexes have folding pathways that can be aborted at any time. Moreover, efficiency and rate of folding are substantially affected or even fully blocked by mutations in the polypeptide chain or by perturbation of the environmental and homeostatic conditions of the folding environment (e.g., the ER) making it very difficult to make trustable predictions (Figure 1). As such, although it is generally thought that the folding process is efficient for most proteins [34,35] there are controversial results/opinions. For example a study reported that ca. 30% of the newly synthesized proteome was rapidly degraded by the proteasome as a possible indication of low efficiency [36], whereas another reported diametrically different conclusions [37].
In order to perform their functions, proteins and protein complexes must be transported to the right intra-or extra-cellular place. Protein targeting is as important as any other step in protein synthesis and as such, defect in protein localization have been linked to several medical conditions or diseases [38]. Approximately half of the synthesized cellular proteins need to be transported across or into a membrane [39] and regulation of most trafficking processes is based on specific signal sequences displayed by the polypeptide chains [40]. Protein transport within cells is carried out by both bulk (e.g., unselective ER export via COPII vesicles) and specific, receptor based mechanisms that enrich select cargos into COPII vesicles (e.g., ERGIC53-based export of secretory glycoproteins from the ER) [41,42]. Protein translocation across organelle membranes might happen post-translationally, as for mitochondrial, peroxisomal, secretory and nuclear proteins, or co-translationally, as in the case of proteins of the secretory pathway and integral membrane proteins, which are synthesized directly into the ER. Protein trafficking into organelles involves specific mechanisms for each compartment. As a general rule, proteins carry a specific address tag for the organelle they are sent to. This tag is recognized by a receptor (bound on the organelle membrane or soluble), which engages the necessary import machinery. Protein import is energy driven and occurs through membrane-spanning channels (e.g., the TOM and TIM complexes of the outer and inner mitochondrial membrane [43] or the Sec61/Sec62/Sec63 complex of the ER membrane [44]). In the case of the nuclear envelope, small proteins up to 40 kDa can also traffic by simple diffusion [45]. The signal sequence may, or may not be removed upon translocation by signal peptidases. Studies on the fidelity of the cellular trafficking processes have not been performed, but several factors have been shown to be important and their loss significantly affects transport efficacy and precision and ultimately leads to disease [38]. As an example, two distinct mutations in the ERGIC-53 gene were shown to underlie the autosomal recessive bleeding disorder caused by combined factor V and factor VIII deficiency. Homozygous patients completely lack the expression of coagulation factors V and VIII due to their impaired export from the ER via COPII vesicles [46,47].
Protein folding occurs in cellular compartments with different chimico-physical properties (i.e., cytosol, mitochondria, ER), which contain compartment-specific chaperones and enzymes that maximize protein folding efficiency. Cells have developed dedicated mechanisms that screen the quality of the proteome during its synthesis and intervene when a protein is not able to reach its native conformation by promoting degradation of the erroneous product. Degradation of cellular proteins generally occurs in the cytosol and nucleus via the ubiquitin proteasome system (UPS) [48,49,50,51,52] and additionally by lysosomal pathways upon engulfment of the cytosolic proteins or organelles by the process of autophagy [51,52,53]. Following chapters will describe quality control and degradation machineries that act during the proteins synthesis process, from translation to protein folding and secretion.
The quality of newly synthesized proteins is constantly monitored from the moment they start emerging from the ribosome up to their degradation [54]. The earliest defective ribosomal products (DRiPs) can be selected for degradation when they are still emerging from the ribosome. This mechanism of co-translational degradation takes place when the translation machinery is defective or engages faulty mRNAs, which leads to ribosomal/translational stalling and the formation of incomplete nascent chains (Figure 3).
Recent work in yeast showed that stalled 80S ribosomes carrying unfinished polypeptide chains are engaged by the ribosome quality control complex (RQC) which mediates ubiquitination and degradation of the incomplete nascent chain [55,56,57,58]. Molecular components of the RQC are the E3 ubiquitin ligase Listerin (Ltn1), Tae2, Rqc1, Rqc2 and Cdc48 with its cofactors Ufd1 and Npl14 [55,59]. Importantly, a recent study described a novel mechanism by which the RQC, in association with the 60S ribosomal subunit, adds polypeptide stretches to the stalled nascent chain and eventually ubiquitinates and targets the polypeptide for degradation via the proteasome [59] (Figure 3). In particular, ribosome stalling leads to the dissociation of the 60S and 40S subunits (Figure 3, step 1), followed by association of Ltn1 and Rqc2 with the peptidyl-tRNA-60S species (Figure 3, step 2). Here Rqc2 through binding of Alanyl-and Threonyl-tRNA orchestrates the, nota bene mRNA free and 40S free, incorporation of carboxy-terminal Ala and Thr extension (termed CAT tails) into the incomplete protein, followed by subsequent ubiquitination by Ltn1 (Figure 3, step 3) and extraction through the Cdc48 complex (Figure 3, step 4).
The ER plays a crucial role in the maintenance of a quantitatively and qualitatively correct proteome. In fact, it is the organelle where the synthesis and folding of most membrane bound and secreted proteins (which make up ca. 1/3 of a cell’s proteome [60]) occurs (Figure 4). Most proteins that are synthesized into the ER (Figure 4, step 1) may be subjected to the addition of a complex glycan, composed of three residues of glucose, nine mannoses and two N-acetylglucosamines (Glc3Man9GlcNAc2) on side chains of asparagine residues in the sequon Asn-Xxx-Ser/Thr/Cys (Xxx all but Pro or Asp). N-glycosylation is catalyzed by a multiprotein complex oligosaccharyltransferase (OST) [61,62] (Figure 4, step 2). The attached sugar increases solubility of nascent and folding polypeptides since it shields hydrophobic regions. Moreover, it engages folding, quality control and degradation factors to determine the fate of the folding protein [62,63,64]. For non-glycosylated proteins the processes that regulate folding and degradation are far less known and thus not mentioned in the following chapters.
Early upon addition of the glycan, the outermost glucose is trimmed by the ER enzyme α-glucosidase I. The di-glucosylated protein-bound oligosaccharide engages Malectin [65] (Figure 4, step 3) a type I membrane-bound ER lectin, which was shown to preferentially associate with misfolded proteins [66]. Malectin operates the first quality control checkpoint for sequestration and delivery of misfolded proteins for ER-associated degradation (ERAD, see below). Elimination of the second glucose by α-glucosidase II leads to association of the maturing protein with the lectins Calnexin (CNX) and Calreticulin (CRT), which in turn engage enzymes necessary for proper folding of the bound polypeptide: protein disulfide isomerases (PDIs) like PDI, ERp57, ERp44, TMX1[67] and others, and peptidyl-prolyl isomerases (PPIs) like cyclophilin B (CYPB) (for a review see [62]) (Figure 4, step 4). Trimming of the third glucose by α-glucosidase II releases the substrate from CNX/CRT (Figure 4, step 5). This is a crucial step in protein folding in the ER. In fact, if a protein collapses to its proper conformation after release from CNX/CRT, it will be transported through the Golgi apparatus to reach its final intra-or extra-cellular destination [41,62] (Figure 4, step 6). On the contrary, if a protein is still misfolded after release from Cnx/Crt there are two options: (1) the N-glycans are re-glucosylated by the glucosyltransferase UGT1 and the polypeptide enters another round of folding in association with CNX/CRT (hence CNX/CRT cycle, [62,64,68]) (Figure 4, step 7), or (2) the N-glycans get further trimmed by ER mannosidases such as ER α1, 2-mannosidase I [69], EDEM1 [70], EDEM2 [71] and/or EDEM3 [72,73]. Removal of the mannoses impedes the re-glucosylation by UGT1 [70] and enables the recruitment of OS-9 and XTP3B [74,75], two lectins that deliver misfolded proteins to retro-translocation channels (dislocons) built around E3 ubiquitin ligases such as HRD1 and GP78 embedded in the ER membrane [76] (Figure 4, step 8). Proteins that escape retention-based quality control in the ER might be retained in post-ER compartments. This relies on the intervention of UGT1 and p97 as has recently been shown for proteins displaying ionizable residues in the intramembrane domain [77] (Figure 4, step 9).
Since protein folding is liable to error and its efficiency differs from protein to protein, misfolded products constantly form in the ER. These must be rapidly eliminated to avoid the formation of toxic aggregates and disturbances of ER homeostasis. Misfolded proteins generally undergo selection for degradation via de-mannosylation of their N-glycans by EDEM proteins and associate with the ER lectins OS-9 and XTP3B [78]. As mentioned above, these two ERAD lectins act as shuttles that deliver the misfolded polypeptide to dislocons. Once there, misfolded proteins are retrotranslocated into the cytosol, poly-ubiquitinated by E3 ligases, extracted from the membrane by the AAA-ATPase p97 and its co-factors, and eventually targeted to the 26S proteasome for extensive proteolysis [78,79] (Figure 4, step 8). ERAD activity must be carefully regulated to prevent excessive or inefficient degradation, which can in turn result in the onset of loss-or gain-of-function disorders, respectively. Indeed, the turnover and the intracellular level of some ERAD factors such as EDEM1, OS-9 and SEL1L, is controlled in a post-translational manner through what has been defined as ERAD tuning [79,80,81].
As mentioned above, degradation of misfolded translational products not only happens through the UPS, but also via autophagy. Autophagy is a cellular degradative pathway that removes proteins, aggregates and even entire organelles, by engulfing them into double-membrane vesicles (the autophagosomes) that then fuse with lysosomes containing hydrolases (Figure 4, step 9). There are several types of autophagy (e.g., macroautophagy, microautophagy and chaperone mediated autophagy (CMA), for a review see [82]). They are mediated by the autophagy-related genes (ATG) and their associated enzymes [83,84,85]. Autophagy was first studied as a response to starvation, but several lines of research in the last decades made it clear that autophagy exists in many flavors and plays crucial roles in cell and organelle homeostasis [86] and, as such, is involved in human health and disease [87].
In line with this, the degradation of misfolded proteins through macroautophagy or both UPS and macroautophagy has been reported. Examples of disease-related proteins degraded via the autophagic machinery are mutant Huntingtin (HTT) protein (Huntington’s disease) [88], α-synuclein (Parkinson’s disease) [89], Aβ and tau (Alzheimer’s disease) [90,91]. Other well-characterized autophagy substrate is the Z-variant of α-1-antitrypsin (ATZ) [92] and cystic-fibrosis conductance regulator (CFTR) aggregates [93].
Although autophagic degradation of cytosolic components might happen stochastically as a bulk phenomenon, several types of highly selective autophagy have recently been described. Autophagy receptors are essential for such targeted degradation, which allow the recognition of the specific cargo and its coupling to the autophagic machinery [94]. Specific degradation of parts of the ER (and its content) have been shown to be crucial both during homeostasis [95] and recovery from an unfolded protein response (UPR, see below), during which excess ER chaperones and folding factors are degraded via the lysosome in a Sec62 dependent manner [96].
The maintenance of proteostasis is of fundamental importance for the integrity of cells and organisms. As described above, cells are equipped with quality control and clearance mechanisms that ensure removal of misfolded proteins. These machineries that work efficiently under homeostatic conditions might need to be adjusted to fluctuations in cargo load in order to permit cells to adapt to, and tolerate stressful stimuli. Mutated or defective gene products might lead to loss-of-function or gain-of-toxic-function phenotypes. Fluctuations in protein synthesis due to environmental conditions or pathogen attacks can create conditions of overload of the quality control systems. Pathogens may subvert cellular proteostasis by hijacking the protein synthesis and quality control machineries for the purpose of their reproduction and immune escape [78,97,98]. In addition to quality control and degradation system for monitoring the quality of their proteome, cells have developed responses to various stresses in order to preserve and restore protein homeostasis. If maintenance or recovery of homeostatic conditions is not possible (i.e., in case of pathologies) and cellular viability is threatened, modulation of protein synthesis, protein folding, quality control and/or stress responses via chemical and biological compounds might be an efficient strategy to counteract these diseases.
A striking example of cellular stress response is the unfolded protein response (UPR) elicited by the accumulation of misfolded proteins in the ER (Figure 5). Since the flux of nascent chains into the ER is dynamic and might change due to fluctuating environmental conditions, genetic mutations and attack by pathogens, a cell has to tune the level of its folding, quality control and ERAD components in order to preserve the production of membrane and secreted cargo proteins. Such homeostatic control is achieved by the constant monitoring of “ER load” by ER membrane embedded stress sensors, namely inositol-requiring-protein-1 (IRE1), activating transcription factor-6 (ATF6) and protein kinase RNA (PKR)-like ER kinase (PERK). These sensors are activated upon accumulation of misfolded proteins in the ER and subsequently orchestrate a signaling cascade from the ER to the nucleus [99] (Figure 5).
The IRE1 branch of the UPR is the most conserved (found in yeast, plants and metazoans [100,101]). IRE1 is a type 1 protein with a lumenal N-terminal domain and a cytosolic C-terminal region that possesses kinase and RNase activities. Mammalian cells possess two homologs, namely IRE1α and IRE1β characterized by tissue-specific distribution. IRE1α is activated upon accumulation of misfolded proteins, which leads to its oligomerization, trans-autophosphorylation with consequent activation of its RNase domain. IRE1α signaling occurs via the cleavage of its principal substrate, XBP1 (X-box binding protein-1) mRNA (Figure 5). An unconventional splicing event leads to the excision of an intron and subsequent ligation leads to the formation of the XBP1s mRNA that translates into a potent bZIP transcription factor that upregulates transcription of ER quality control components and enzymes for lipid synthesis, which controls ER membrane expansion [102] (Figure 5). In addition to XBP1 splicing, IRE1 is responsible for the cleavage of a specific subset of mRNAs, mechanism termed Regulated IRE1-Dependent Decay (RIDD), leading to the reduction of ER cargo protein synthesis [103] (Figure 5).
ATF6 is a metazoan specific stress sensor with a type II topology consisting of a luminal carboxyl-terminal stress sensing domain and a cytoplasmic amino-terminal bZIP transcription factor domain. Upon ER stress, ATF6 is transported to the Golgi apparatus, where it is sequentially cleaved by site 1 and site 2 proteases (S1P and S2P) leading to the release of the bZIP domain (ATF6f), which translocates to the nucleus and enhances transcription of ER quality control and ERAD components [102,104] (Figure 5). Interestingly, XBP1s and ATF6 upregulate UPR target genes with similar functions (e.g., chaperones and PDIs) and even share some of their targets, which upregulation mainly aims at enhancing folding and degradation capacity of the ER [102,105]. Additionally XBP1 is also a target gene of the ATF6 pathway, thus tightly linking the IRE1 and the ATF6 branches of the UPR.
The third branch of the UPR differs in its outcome compared to the IRE1 and ATF6 pathways. PERK is a type I membrane protein with a luminal stress sensing domain, similar to the one of IRE1 [106], and a cytosolic kinase domain. Upon activation by ER stress, PERK oligomerizes and phosphorylates itself and other targets. Importantly, PERK phosphorylates the α-subunit of eukaryotic initiation factor 2 (eIF2α) at Ser51, which results in a diminished global protein synthesis, thus reducing burden of ER cargo proteins [107]. Phosphorylation of eIF2α specifically upregulates the translation of activating transcription factor 4 (ATF4) from an open reading frame in its 5’ untranslated region [108,109]. Important among ATF4 target genes are activating transcription factor 3 (ATF3), transcription factor C/EBP homologous protein (CHOP) and DNA damage-inducible 34 (GADD34). ATF3 targets both CHOP and GADD34 [110] (Figure 5). CHOP is a transcription factor responsible for the upregulated transcription of several apoptotic genes, while GADD34, induced late during ER stress, encodes a phosphatase, which counteracts eIF2α phosphorylation (reducing translational inhibition). Thus, the PERK pathway plays a dual role during UPR: induction of apoptosis if the ER stress cannot be resolved and recovery of protein synthesis if eIF2α phosphorylation is reduced via both reduction of ER stress (decreased PERK activation) and GADD34. Additionally, PERK phosphorylates NRF2 which induces genes responsible for antioxidant responses [111] (Figure 5)
Several cellular machineries ensure the synthesis of the proper amount and quality of the proteome. Moreover, quality control mechanisms operate at every stage of protein biogenesis in order to maintain proteostasis and to respond to perturbations of various origins by inducing stress responses.
Mutations in the genome can give rise to ribosomal products that can elude the quality control systems. Amino acid substitutions in a newly synthesized polypeptide might destabilize the entire protein’s structure [112] causing protein aggregation and deposition (gain-of-toxic-function), or resulting in rapid degradation of the aberrant product (loss-of-function) [113,114,115]. Such pathological conditions are referred to as conformational diseases [114,116]. The aggregation or misfolding of a single species of proteins can affect maturation of other cellular proteins that follow or engage similar pathways. This means that defects in a single protein might spread to several other protein species and induce their misfolding and aggregation, thus harnessing the entire proteostasis network [117,118]. If activation of stress response pathways cannot re-establish proteostasis, consequences will be deleterious for single cells, organs and the entire organism.
Therapeutic approaches are available to restore the folding of proteins or to inhibit their aggregation by delivering molecules that facilitate folding of mutant proteins, stabilize their structure or enhance their solubility. Chemical chaperones are membrane permeable molecules that may non-selectively rescue the folding and functionality of mutant misfolded proteins [114]. Pharmacological chaperones, which have the same biological function, specifically bind only to certain proteins. The specificity makes pharmacological chaperones very suitable for clinical use [119,120,121].
The importance of the UPR is underlined by its implication in many diseases. Cancer cells undergo somatic mutations that are not detected by their DNA damage machineries, which enable them to ignore growth control and disable apoptotic signaling. Especially solid-tumor cells, which are poorly vascularized, might also encounter other type of stresses such as hypoxia, nutrient deprivation and pH changes [122]. These conditions can alter ER homeostasis and interfere with the physiological protein folding environment leading to the activation of the UPR. Activation of the IRE1/XBP1s and the ATF6 branches has been shown in several cancer types such as breast tumors [123], hepatocellular carcinomas [124], gastric tumors [125] and several cancer cell lines [126]. ER chaperone levels were shown to correlate with tumor aggressiveness and alterations in reaction to chemical agents. It is not clear if apoptotic signaling induced by the UPR is escaped due to specific activation of the IRE1 and ATF6 branches or to inactivation of the apoptotic program. In any case, cytoprotective branches of the UPR enable cancer cells to produce and fold large numbers of the proteins that they produce. Thus inhibition of the UPR might be a promising way to treat some types of cancer [127].
In contrast, other conditions are linked to an insufficient activation of the UPR, meaning that induction of the UPR, especially the IRE1 and/or ATF6 branches, could increase ER protein folding and export capacity decreasing aberrant ERAD of unstable, mutant proteins, thus ameliorating loss-of-function diseases such as cystic fibrosis or lysosomal storage diseases [128,129,130]. Similarly, increasing ERAD activity could attenuate the secretion of destabilized, aggregation-prone proteins as in the case of human amyloid diseases [113,131,132,133]. UPR signaling pathways are thereby promising targets for pharmacological intervention in a heterogeneous variety of human diseases. Small molecules, which selectively activate or inhibit only one branch of the UPR have been described and successfully used to influence pathological states [134,135,136,137,138]. Most recently, ATF6 specific modulators were described. Ceapins selectively inhibit the ATF6 pathway by blocking ATF6 delivery (and cleavage) to the Golgi apparatus [139,140]. Moreover, several ATF6 activators that enable the reduction of extracellular protein aggregates were described [141]. The palette of available chemical compounds offers the possibility to specifically modulate one or more branches of the UPR (at least in cultured cells) paving the road for the development of new therapeutic options for a variety of medical conditions including cancer and conformational diseases.
Maintenance of proteostasis is of fundamental importance for cellular viability. Through the combined action of a highly interconnected network of signaling pathways, which modulate gene expression, translation, protein folding, trafficking and degradation, cells tightly control the abundance and quality of the produced proteome and adapt it to stressful conditions. Nevertheless, pathological conditions, such as protein conformational diseases, may arise due to defective protein folding and transport.
The social and medical impact of protein misfolding diseases (most often rare genetic conditions that affect children) highlight the importance of understanding basic molecular mechanisms of cellular functions. Common pathways can be engaged by gene products at the basis of completely different pathological phenotypes (e.g. cancer vs lysosomal storage diseases). It will therefore be crucial in future to continue working on mechanisms involved in protein folding, ERAD and UPR in order to describe new targets for pharmaceutical intervention. Several open questions remain to be answered: What are the chimico-physical features of a protein that determine the engagement of specific folding or clearance machineries? How important are certain mechanisms in different cell types or tissues and how differently are they regulated? Do cells experiencing stresses communicate with and influence surrounding cells/tissues? Which are the components involved in the recovery from stresses? Do cells remember experiencing stressful conditions in order to adapt to future insults (epigenetics)? And of course how can these pathways be targeted in order to restore proteostasis?
To the latter question the recently developed clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 (CRISPR associated protein 9) gene editing technology seems very promising to correct gene defects [142]. Gene-editing as a therapeutic option has the advantage of eliminating the problem at its source, but what if the intervention causes unpredictable side effects? In contrast to pharmacologic interventions, whose use can be promptly interrupted in case of unwanted side effects, irreversibility is the major drawback in “genetic surgery”. Thus, basic research aimed at identify new possible targets regulating processes involved in protein biogenesis (in its larger meaning) remains of fundamental importance.
M. M. is supported by Signora Alessandra, the Foundation for Research on Neurodegenerative Diseases, the Swiss National Science Foundation (SNF), the Sinergia program of the SNF and the Comel, Gabriele and Gelu Foundations.
The authors declare no conflict of interest.
[1] | J. Liu, H. Dang, X. W. Wang, The significance of intertumor and intratumor heterogeneity in liver cancer, Exp. Mol. Med., 50 (2018), e416. |
[2] |
T. M. Grzywa, W. Paskal, P. K. Włodarski, Intratumor and intertumor heterogeneity in melanoma, Transl. Oncol., 10 (2017), 956-975. doi: 10.1016/j.tranon.2017.09.007
![]() |
[3] | K. Anderson, C. Lutz, F. W. Van Delft, C. M. Bateman, Y. Guo, S. M. Colman, et al., Genetic variegation of clonal architecture and propagating cells in leukaemia, Nature, 469 (2011), 356-361. |
[4] | S. Li, F. E. Garrett-Bakelman, S. S. Chung, M. A. Sanders, T. Hricik, F. Rapaport, et al., Distinct evolution and dynamics of epigenetic and genetic heterogeneity in acute myeloid leukemia, Nat. Med., 22 (2016), 792-799. |
[5] | G. C. Webb, D. D. Chaplin, Genetic variability at the human tumor necrosis factor loci., J. Immunol., 145 (1990), 1278-1285. |
[6] |
A. Marusyk, K. Polyak, Tumor heterogeneity: Causes and consequences, Biochim. Biophys. Acta, Rev. Cancer, 1805 (2010), 105-117. doi: 10.1016/j.bbcan.2009.11.002
![]() |
[7] | A. A. Alizadeh, V. Aranda, A. Bardelli, C. Blanpain, C. Bock, C. Borowski, et al., Toward understanding and exploiting tumor heterogeneity, Nat. Med., 21 (2015), 846-853. |
[8] | M. Lee, G. T. Chen, E. Puttock, K. Wang, R. A. Edwards, M. L. Waterman, et al., Mathematical modeling links wnt signaling to emergent patterns of metabolism in colon cancer, Mol. Syst. Biol., 13 (2017), 912. |
[9] |
S. Terry, S. Buart, S. Chouaib, Hypoxic stress-induced tumor and immune plasticity, suppression, and impact on tumor heterogeneity, Front. Immunol., 8 (2017), 1625. doi: 10.3389/fimmu.2017.01625
![]() |
[10] | S. Akgul, A.-M. Patch, R. C. J. D'Souza, P. Mukhopadhyay, K. Nones, S. Kempe, et al., ¨ Intratumoural heterogeneity underlies distinct therapy responses and treatment resistance in glioblastoma, Cancers, 11 (2019), 190. |
[11] |
T. D. Laajala, T. Gerke, S. Tyekucheva, J. C. Costello, Modeling genetic heterogeneity of drug response and resistance in cancer, Curr. Opin. Syst. Biol., 17 (2019), 8-14. doi: 10.1016/j.coisb.2019.09.003
![]() |
[12] |
X. Ma, J. Huang, Predicting clinical outcome of therapy-resistant prostate cancer, Proc. Natl. Acad. Sci., 116 (2019), 11090-11092. doi: 10.1073/pnas.1906812116
![]() |
[13] |
A. M. Stein, T. Demuth, D. Mobley, M. Berens, L. M. Sander, A mathematical model of glioblastoma tumor spheroid invasion in a three-dimensional in vitro experiment, Biophys. J., 92 (2007), 356-365. doi: 10.1529/biophysj.106.093468
![]() |
[14] |
T. L. Stepien, E. M. Rutter, Y. Kuang, Traveling waves of a go-or-grow model of glioma growth, SIAM J. Appl. Math., 78 (2018), 1778-1801. doi: 10.1137/17M1146257
![]() |
[15] |
A. Martínez-González, G. F. Calvo, L. A. P. Romasanta, V. M. Pérez-García, Hypoxic cell waves around necrotic cores in glioblastoma: A biomathematical model and its therapeutic implications, Bull. Math. Biol., 74 (2012), 2875-2896. doi: 10.1007/s11538-012-9786-1
![]() |
[16] |
H. Hatzikirou, D. Basanta, M. Simon, K. Schaller, A. Deutsch, 'Go or grow': The key to the emergence of invasion in tumour progression?, Math. Med. Biol.: A J. IMA, 29 (2012), 49-65. doi: 10.1093/imammb/dqq011
![]() |
[17] |
Z. Husain, P. Seth, V. P. Sukhatme, Tumor-derived lactate and myeloid-derived suppressor cells: Linking metabolism to cancer immunology, Oncoimmunology, 2 (2013), e26383. doi: 10.4161/onci.26383
![]() |
[18] |
A. N. Mendler, B. Hu, P. U. Prinz, M. Kreutz, E. Gottfried, E. Noessner, Tumor lactic acidosis suppresses CTL function by inhibition of p38 and JNK/c-Jun activation, Int. J. Cancer, 131 (2012), 633-640. doi: 10.1002/ijc.26410
![]() |
[19] | K. Fischer, P. Hoffmann, S. Voelkl, N. Meidenbauer, J. Ammer, M. Edinger, et al., Inhibitory effect of tumor cell-derived lactic acid on human T cells, Blood, 109 (2007), 3812-3819. |
[20] |
E. C. Fiedler, M. T. Hemann, Aiding and abetting: How the tumor microenvironment protects cancer from chemotherapy, Annu. Rev. Cancer Biol., 3 (2019), 409-428. doi: 10.1146/annurev-cancerbio-030518-055524
![]() |
[21] | J. K. Saggar, M. Yu, Q. Tan, I. F. Tannock, The tumor microenvironment and strategies to improve drug distribution, Front. Oncol., 3 (2013), 154. |
[22] |
A. I. Minchinton, I. F. Tannock, Drug penetration in solid tumours, Nat. Rev. Cancer, 6 (2006), 583-592. doi: 10.1038/nrc1893
![]() |
[23] | A. Baldock, R. Rockne, A. Boone, M. Neal, C. Bridge, L. Guyman, et al., From patient-specific mathematical neuro-oncology to precision medicine, Front. Oncol., 3 (2013), 62. |
[24] | A. L. Baldock, S. Ahn, R. Rockne, S. Johnston, M. Neal, D. Corwin, et al., Patient-specific metrics of invasiveness reveal significant prognostic benefit of resection in a predictable subset of gliomas, PLoS One, 9 (2014), e99057. |
[25] |
S. Barish, M. F. Ochs, E. D. Sontag, J. L. Gevertz, Evaluating optimal therapy robustness by virtual expansion of a sample population, with a case study in cancer immunotherapy, Proc. Nati. Acad. Sci., 114 (2017), E6277-E6286. doi: 10.1073/pnas.1703355114
![]() |
[26] | J. C. L. Alfonso, L. Berk, Modeling the effect of intratumoral heterogeneity of radiosensitivity on tumor response over the course of fractionated radiation therapy, Radiat. Oncol., 14 (2019), 88. |
[27] |
F. Forouzannia, H. Enderling, M. Kohandel, Mathematical modeling of the effects of tumor heterogeneity on the efficiency of radiation treatment schedule, Bull. Math. Biol., 80 (2018), 283-293. doi: 10.1007/s11538-017-0371-5
![]() |
[28] |
C. H. Morrell, J. D. Pearson, H. B. Carter, L. J. Brant, Estimating unknown transition times using a piecewise nonlinear mixed-effects model in men with prostate cancer, J. Am. Stat. Assoc., 90 (1995), 45-53. doi: 10.1080/01621459.1995.10476487
![]() |
[29] | N. Hartung, S. Mollard, D. Barbolosi, A. Benabdallah, G. Chapuisat, G. Henry, et al., Mathematical modeling of tumor growth and metastatic spreading: Validation in tumor-bearing mice, Cancer Res., 74 (2014), 6397-6407. |
[30] |
J. Poleszczuk, R. Walker, E. G. Moros, K. Latifi, J. J. Caudell, H. Enderling, Predicting patientspecific radiotherapy protocols based on mathematical model choice for proliferation-saturation index, Bull. Math. Biol., 80 (2018), 1195-1206. doi: 10.1007/s11538-017-0279-0
![]() |
[31] | A. Rafii, C. Touboul, H. A. Thani, K. Suhre, J. A. Malek, Where cancer genomics should go next: a clinician's perspective, Hum. Mol. Genet., 23 (2014), R69-R75. |
[32] |
D. Zardavas, A. Irrthum, C. Swanton, M. Piccart, Clinical management of breast cancer heterogeneity, Nat. Rev. Clin. Oncol., 12 (2015), 381-394. doi: 10.1038/nrclinonc.2015.73
![]() |
[33] |
A. D. Barker, C. C. Sigman, G. J. Kelloff, N. M. Hylton, D. A. Berry, L. J. Esserman, I-spy 2: an adaptive breast cancer trial design in the setting of neoadjuvant chemotherapy, Clin. Pharmacol. Ther., 86 (2009), 97-100. doi: 10.1038/clpt.2009.68
![]() |
[34] |
A. S. Meyer, L. M. Heiser, Systems biology approaches to measure and model phenotypic heterogeneity in cancer, Curr. Opin. Syst. Biol., 17 (2019), 35-40. doi: 10.1016/j.coisb.2019.09.002
![]() |
[35] | R. A. Gatenby, E. T. Gawlinski, A reaction-diffusion model of cancer invasion, Cancer Res., 56 (1996), 5745-5753. |
[36] |
P. Tracqui, G. C. Cruywagen, D. E. Woodward, G. T. Bartoo, J. D. Murray, E. Alvord Jr, A mathematical model of glioma growth: the effect of chemotherapy on spatio-temporal growth, Cell Proliferation, 28 (1995), 17-31. doi: 10.1111/j.1365-2184.1995.tb00036.x
![]() |
[37] |
M. A. J. Chaplain, Reaction-diffusion prepatterning and its potential role in tumour invasion, J. Biol. Syst., 3 (1995), 929-936. doi: 10.1142/S0218339095000824
![]() |
[38] |
H. L. Harpold, E. C. Alvord Jr, K. R. Swanson, The evolution of mathematical modeling of glioma proliferation and invasion, J. Neuropathol. Exp. Neurol., 66 (2007), 1-9. doi: 10.1097/nen.0b013e31802d9000
![]() |
[39] | R. Rockne, J. K. Rockhill, M. Mrugala, A. M. Spence, I. Kalet, K. Hendrickson, et al., Predicting the efficacy of radiotherapy in individual glioblastoma patients in vivo: a mathematical modeling approach, Phys. Med. Biol., 55 (2010), 3271. |
[40] | T. E. Yankeelov, N. Atuegwu, D. Hormuth, J. A. Weis, S. L. Barnes, M. I. Miga, et al., Clinically relevant modeling of tumor growth and treatment response, Sci. Transl. Med., 5 (2013), 187ps9-187ps9. |
[41] | D. A. Hormuth II, J. A. Weis, S. L. Barnes, M. I. Miga, E. C. Rericha, V. Quaranta, et al., Predicting in vivo glioma growth with the reaction diffusion equation constrained by quantitative magnetic resonance imaging data, Phys. Biol., 12 (2015), 046006. |
[42] |
K. R. Swanson, R. C. Rostomily, E. C. Alvord Jr, A mathematical modelling tool for predicting survival of individual patients following resection of glioblastoma: A proof of principle, Br. J. Cancer, 98 (2008), 113-119. doi: 10.1038/sj.bjc.6604125
![]() |
[43] | S. Prokopiou, E. G. Moros, J. Poleszczuk, J. Caudell, J. F. Torres-Roca, K. Latifi, et al., A proliferation saturation index to predict radiation response and personalize radiotherapy fractionation, Radiat. Oncol., 10 (2015), 159. |
[44] | C. Vaghi, A. Rodallec, R. Fanciullino, J. Ciccolini, J. Mochel, M. Mastri, et al., Population modeling of tumor growth curves, the reduced Gompertz model and prediction of the age of a tumor, in Mathematical and Computational Oncology (eds. G. Bebis, T. Benos, K. Chen, K. Jahn, E. Lima), Lecture Notes in Computer Science, Springer International Publishing, Cham, (2019), 87-97. |
[45] | T. J. Sullivan, Introduction to Uncertainty Quantification, Springer, 2015. |
[46] | G. J. Lord, C. E. Powell, T. Shardlow, Introduction to Computational Stochastic PDEs, Cambridge, 2014. |
[47] | R. Smith, Uncertainty Quantification: Theory, Implementation, and Applications, SIAM, 2014. |
[48] | N. Henscheid, E. Clarkson, K. J. Myers, H. H. Barrett, Physiological random processes in precision cancer therapy, PLoS One, 13 (2018), e0199823. |
[49] |
J. P. O'Connor, C. J. Rose, J. C. Waterton, R. A. Carano, G. J. Parker, A. Jackson, Imaging intratumor heterogeneity: Role in therapy response, resistance, and clinical outcome, Clin. Cancer Res., 21 (2015), 249-257. doi: 10.1158/1078-0432.CCR-14-0990
![]() |
[50] |
R. J. Allen, T. R. Rieger, C. J. Musante, Efficient generation and selection of virtual populations in quantitative systems pharmacology models, CPT: Pharmacometrics Syst. Pharmacol., 5 (2016), 140-146. doi: 10.1002/psp4.12063
![]() |
[51] | T. R. Rieger, R. J. Allen, L. Bystricky, Y. Chen, G. W. Colopy, Y. Cui, et al., Improving the generation and selection of virtual populations in quantitative systems pharmacology models, Prog. Biophys. Mol. Biol., 139 (2018), 15-22. |
[52] |
D. J. Klinke, Integrating epidemiological data into a mechanistic model of type 2 diabetes: Validating the prevalence of virtual patients, Ann. Biomed. Eng., 36 (2008), 321-334. doi: 10.1007/s10439-007-9410-y
![]() |
[53] |
B. J. Schmidt, F. P. Casey, T. Paterson, J. R. Chan, Alternate virtual populations elucidate the type i interferon signature predictive of the response to rituximab in rheumatoid arthritis, BMC Bioinf., 14 (2013), 221. doi: 10.1186/1471-2105-14-221
![]() |
[54] | J. Bassaganya-Riera, Accelerated Path to Cures, Springer, 2018. |
[55] | N. Henscheid, Quantifying Uncertainties in Imaging-Based Precision Medicine, Ph.D thesis, University of Arizona, 2018. |
[56] | J. G. Albeck, J. M. Burke, S. L. Spencer, D. A. Lauffenburger, P. K. Sorger, Modeling a snapaction, variable-delay switch controlling extrinsic cell death, PLoS Biol., 6 (2008), e299. |
[57] |
Q. Wu, S. D. Finley, Predictive model identifies strategies to enhance tsp1-mediated apoptosis signaling, Cell Commun. Signaling, 15 (2017), 53. doi: 10.1186/s12964-017-0207-9
![]() |
[58] | M. R. Birtwistle, J. Rauch, A. Kiyatkin, E. Aksamitiene, M. Dobrzyński, J. B. Hoek, et al., Emergence of bimodal cell population responses from the interplay between analog single-cell signaling and protein expression noise, BMC Syst. Biol., 6 (2012), 109. |
[59] | D. Brown, R. A. Namas, K. Almahmoud, A. Zaaqoq, J. Sarkar, D. A. Barclay, et al., Trauma in silico: Individual-specific mathematical models and virtual clinical populations, Sci. Transl. Med., 7 (2015), 285ra61-285ra61. |
[60] | A. Badano, C. G. Graff, A. Badal, D. Sharma, R. Zeng, F. W. Samuelson, et al., Evaluation of digital breast tomosynthesis as replacement of full-field digital mammography using an in silico imaging trial, JAMA Net. Open, 1 (2018), e185474-e185474. |
[61] | W. Kainz, E. Neufeld, W. E. Bolch, C. G. Graff, C. H. Kim, N. Kuster, et al., Advances in computational human phantoms and their applications in biomedical engineering—a topical review, IEEE Trans. Radiat. Plasma Med. Sci., 3 (2019), 1-23. |
[62] | E. Roelofs, M. Engelsman, C. Rasch, L. Persoon, S. Qamhiyeh, D. De Ruysscher, et al., Results of a multicentric in silico clinical trial (rococo): Comparing radiotherapy with photons and protons for non-small cell lung cancer, J. Thorac. Oncol., 7 (2012), 165-176. |
[63] |
J. K. Birnbaum, F. O. Ademuyiwa, J. J. Carlson, L. Mallinger, M. W. Mason, R. Etzioni, Comparative effectiveness of biomarkers to target cancer treatment: Modeling implications for survival and costs, Med. Decis. Making, 36 (2016), 594-603. doi: 10.1177/0272989X15601998
![]() |
[64] |
D. Li, S. D. Finley, The impact of tumor receptor heterogeneity on the response to anti-angiogenic cancer treatment, Integr. Biol., 10 (2018), 253-269. doi: 10.1039/C8IB00019K
![]() |
[65] |
K. Mustapha, Q. Gilli, J.-M. Frayret, N. Lahrichi, E. Karimi, Agent-based simulation patient model for colon and colorectal cancer care trajectory, Procedia Comput. Sci., 100 (2016), 188-197. doi: 10.1016/j.procs.2016.09.140
![]() |
[66] |
J. P. Rolland, H. H. Barrett, Effect of random background inhomogeneity on observer detection performance, J. Opt. Soc. Am. A, 9 (1992), 649-658. doi: 10.1364/JOSAA.9.000649
![]() |
[67] | H. H. Barrett, K. J. Myers, Foundations of Image Science, NY: John Wiley and Sons, New York, 2004. |
[68] |
M. A. Kupinski, E. Clarkson, J. W. Hoppin, L. Chen, H. H. Barrett, Experimental determination of object statistics from noisy images, JOSA A, 20 (2003), 421-429. doi: 10.1364/JOSAA.20.000421
![]() |
[69] |
F. Bochud, C. Abbey, M. Eckstein, Statistical texture synthesis of mammographic images with clustered lumpy backgrounds, Opt. Express, 4 (1999), 33-42. doi: 10.1364/OE.4.000033
![]() |
[70] | N. Henscheid, Generating patient-specific virtual tumor populations with reaction-diffusion models and molecular imaging data, Math. Biosci. Eng., Submitted. |
[71] | W. Hundsdorfer, J. G. Verwer, Numerical solution of time-dependent advection-diffusion-reaction equations, Springer Science and Business Media, 2013. |
[72] | S. KM, M. Prah, Data from brain-tumor-progression, The Cancer Imaging Archive, 2018. |
[73] |
M. Davidian, D. M. Giltinan, Nonlinear models for repeated measurement data: An overview and update, J. Agric., Biol., and Environ. Stat., 8 (2003), 387-419. doi: 10.1198/1085711032697
![]() |
[74] | K. Clark, B. Vendt, K. Smith, J. Freymann, J. Kirby, P. Koppel, et al., The cancer imaging archive (tcia): maintaining and operating a public information repository, J. Digital Imaging, 26 (2013), 1045-1057. |
[75] | S. Desmée, F. Mentré, C. Veyrat-Follet, J. Guedj, Nonlinear mixed-effect models for prostatespecific antigen kinetics and link with survival in the context of metastatic prostate cancer: A comparison by simulation of two-stage and joint approaches, Am. Assoc. Pharm. Sci. J., 17 (2015), 691-699. |
[76] |
A. Król, C. Tournigand, S. Michiels, V. Rondeau, Multivariate joint frailty model for the analysis of nonlinear tumor kinetics and dynamic predictions of death, Stat. Med., 37 (2018), 2148-2161. doi: 10.1002/sim.7640
![]() |
[77] |
E. Grenier, V. Louvet, P. Vigneaux, Parameter estimation in non-linear mixed effects models with SAEM algorithm: extension from ODE to PDE, ESAIM: Math. Modell. Numer. Anal., 48 (2014), 1303-1329. doi: 10.1051/m2an/2013140
![]() |
[78] | H. Liang, Modeling antitumor activity in xenograft tumor treatment, Biom. J.: J. Math. Methods Biosci., 47 (2005), 358-368. |
[79] |
H. Liang, N. Sha, Modeling antitumor activity by using a non-linear mixed-effects model, Math. Biosci., 189 (2004), 61-73. doi: 10.1016/j.mbs.2004.01.002
![]() |
[80] | A. M. Stein, D. Bottino, V. Modur, S. Branford, J. Kaeda, J. M. Goldman, et al., Bcr-abl transcript dynamics support the hypothesis that leukemic stem cells are reduced during imatinib treatment, Clin. Cancer Res., 17 (2011), 6812-6821. |
[81] | T. Ferenci, J. Sápi, L. Kovács, Modelling tumor growth under angiogenesis inhibition with mixed-effects models, Acta Polytech. Hung., 14 (2017), 221-234. |
[82] | B. Ribba, G. Kaloshi, M. Peyre, D. Ricard, V. Calvez, M. Tod, et al., A tumor growth inhibition model for low-grade glioma treated with chemotherapy or radiotherapy, Clin. Cancer Res., 18 (2012), 5071-5080. |
[83] | B. Ribba, N. H. Holford, P. Magni, I. Trocóniz, I. Gueorguieva, P. Girard, et al., A review of mixed-effects models of tumor growth and effects of anticancer drug treatment used in population analysis, CPT: Pharmacometrics Syst. Pharmacol., 3 (2014), 1-10. |
[84] | N. W. Shock, R. C. Greulich, R. Andres, D. Arenberg, P. T. Costa, E. G. Lakatta, et al., Normal human aging. the baltimore longituidinal study of aging, U.S. Government Printing Office Publication. |
[85] | https://www.mathworks.com/help/stats/nlmefit.html, accessed November 16, 2019. |
[86] | H. T. Banks, V. A. Bokil, S. Hu, A. K. Dhar, R. A. Bullis, C. L. Browdy, et al., Modeling shrimp biomass and viral infection for production of biological countermeasures, Math. Biosci. Eng., 3 (2006), 635-660. |
[87] |
H. T. Banks, J. L. Davis, A comparison of approximation methods for the estimation of probability distributions on parameters, Appl. Numer. Math., 57 (2007), 753-777. doi: 10.1016/j.apnum.2006.07.016
![]() |
[88] | M. Sirlanci, S. E. Luczak, C. E. Fairbairn, D. Kang, R. Pan, X. Yu, et al., Estimating the distribution of random parameters in a diffusion equation forward model for a transdermal alcohol biosensor, Automatica, 106 (2019), 101-109. |
[89] |
E. M. Rutter, H. T. Banks, K. B. Flores, Estimating intratumoral heterogeneity from spatiotemporal data, J. Math. Biol., 77 (2018), 1999-2022. doi: 10.1007/s00285-018-1238-6
![]() |
[90] | H. T. Banks, K. B. Flores, I. G. Rosen, E. M. Rutter, M. Sirlanci, W. C. Thompson, The prohorov metric framework and aggregate data inverse problems for random pdes, Commun. Appl. Anal., 22 (2018), 415-446. |
[91] | H. T. Banks, W. C. Thompson, Existence and consistency of a nonparametric estimator of probability measures in the prohorov metric framework, Int. J. Pure Appl. Math., 103 (2015), 819-843. |
[92] | H. Akaike, A new look at the statistical model identification, in Selected Papers of Hirotugu Akaike (eds. E. Parzen, K. Tanabe, G. Kitagawa), Springer New York, New York, (1998), 215-222. |
[93] | R. Ullrich, H. Backes, H. Li, L. Kracht, H. Miletic, K. Kesper, et al., Glioma proliferation as assessed by 3-fluoro-3-deoxy-l-thymidine positron emission tomography in patients with newly diagnosed high-grade glioma, Clin. Cancer Res., 14 (2008), 2049-2055. |
[94] | E. Stretton, E. Geremia, B. Menze, H. Delingette, N. Ayache, Importance of patient dti's to accurately model glioma growth using the reaction diffusion equation, in 2013 IEEE 10th Int. Symp. Biomed. Imaging, IEEE, 2013, 1142-1145. |
[95] | N. C. Atuegwu, D. C. Colvin, M. E. Loveless, L. Xu, J. C. Gore, T. E. Yankeelov, Incorporation of diffusion-weighted magnetic resonance imaging data into a simple mathematical model of tumor growth, Phys. Med. Biol., 57 (2011), 225-240. |
[96] | J. Lipková, P. Angelikopoulos, S. Wu, E. Alberts, B. Wiestler, C. Diehl, et al., Personalized radiotherapy design for glioblastoma: Integrating mathematical tumor models, multimodal scans, and bayesian inference, IEEE Trans. Med. Imaging, 38 (2019), 1875-1884. |
[97] | J. Kaipio, E. Somersalo, Statistical and Computational Inverse Problems, Springer Science and Business Media, 2006. |
[98] |
A. M. Stuart, Inverse problems: a bayesian perspective, Acta Numer., 19 (2010), 451-559. doi: 10.1017/S0962492910000061
![]() |
[99] |
P. R. Jackson, J. Juliano, A. Hawkins-Daarud, R. C. Rockne, K. R. Swanson, Patient-specific mathematical neuro-oncology: using a simple proliferation and invasion tumor model to inform clinical practice, Bull. Math. Biol., 77 (2015), 846-856. doi: 10.1007/s11538-015-0067-7
![]() |
[100] | J. Oden, Cse 397 Lecture Notes: Foundations of Predictive Computational Science, 2017. Available from: https://www.oden.utexas.edu/media/reports/2017/1701.pdf. |
[101] | J. Collis, A. J. Connor, M. Paczkowski, P. Kannan, J. Pitt-Francis, H. M. Byrne, et al., Bayesian calibration, validation and uncertainty quantification for predictive modelling of tumour growth: A tutorial, Bull. Math. Biol., 79 (2017), 939-974. |
[102] |
E. J. Kostelich, Y. Kuang, J. M. McDaniel, N. Z. Moore, N. L. Martirosyan, M. C. Preul, Accurate state estimation from uncertain data and models: An application of data assimilation to mathematical models of human brain tumors, Biol. Direct, 6 (2011), 64. doi: 10.1186/1745-6150-6-64
![]() |
[103] | A. Hawkins-Daarud, S. Prudhomme, K. G. van der Zee, J. T. Oden, Bayesian calibration, validation, and uncertainty quantification of diffuse interface models of tumor growth, J. Math. Biol., 67 (2013), 1457-1485. |
[104] |
I. Babuška, F. Nobile, R. Tempone, A systematic approach to model validation based on bayesian updates and prediction related rejection criteria, Comput. Methods Appl. Mech. Eng., 197 (2008), 2517-2539. doi: 10.1016/j.cma.2007.08.031
![]() |
[105] | J. T. Oden, E. A. B. F. Lima, R. C. Almeida, Y. Feng, M. N. Rylander, D. Fuentes, et al., Toward predictive multiscale modeling of vascular tumor growth, Arch. Comput. Methods Eng., 23 (2016), 735-779. |
[106] | M. Lê, H. Delingette, J. Kalpathy-Cramer, E. R. Gerstner, T. Batchelor, J. Unkelbach, et al., MRI based bayesian personalization of a tumor growth model, IEEE Trans. Med. Imaging, 35 (2016), 2329-2339. |
[107] | A. Hawkins-Daarud, S. K. Johnston, K. R. Swanson, Quantifying uncertainty and robustness in a biomathematical model-based patient-specific response metric for glioblastoma, JCO Clin. Cancer Inf., 3 (2019), 1-8. |
[108] | I. M. Sobol, Sensitivity estimates for nonlinear mathematical models, Math. Modell. Comput. Exp., 1 (1993), 407-414. |
[109] | P. Constantine, Active Subspaces: Emerging Ideas for Dimension Reduction in Parameter Studies, SIAM, 2015. |
[110] |
M. D. Morris, Factorial sampling plans for preliminary computational experiments, Technometrics, 33 (1991), 161-174. doi: 10.1080/00401706.1991.10484804
![]() |
[111] | C. Rasumssen, C. Williams, Gaussian Processes for Machine Learning, MIT Press, 2006. |
[112] |
X. Zhu, M. Welling, F. Jin, J. Lowengrub, Predicting simulation parameters of biological systems using a gaussian process model, Stat. Anal. Data Mining: ASA Data Sci. J., 5 (2012), 509-522. doi: 10.1002/sam.11163
![]() |
[113] | M. C. Kennedy, A. O'Hagan, Bayesian calibration of computer models, J. R. Stat. Soc.: S. B (Stat. Method., 63 (2001), 425-464. |
[114] | Y. LeCun, L. Bottou, Y. Bengio, P. Haffner, et al., Gradient-based learning applied to document recognition, Proc. IEEE, 86 (1998), 2278-2324. |
[115] | A. Krizhevsky, I. Sutskever, G. E. Hinton, Imagenet classification with deep convolutional neural networks, in Advances in Neural Information Processing Systems (eds. F. Pereira, C. J. C. Burges, L. Bottou, K. Q. Weinberger), Curran Associates, (2012), 1097-1105. |
[116] | O. Ronneberger, P. Fischer, T. Brox, U-net: Convolutional networks for biomedical image segmentation, in Medical Image Computing and Computer-Assisted Intervention (eds. N. Navab, J. Hornegger, W. M. Wells, A. F. Frangi), Springer International Publishing, Cham, (2015), 234-241. |
[117] | A. Esteva, A. Robicquet, B. Ramsundar, V. Kuleshov, M. DePristo, K. Chou, et al., A guide to deep learning in healthcare, Nat. Med., 25 (2019), 24-29. |
[118] | Z. Cao, L. Duan, G. Yang, T. Yue, Q. Chen, An experimental study on breast lesion detection and classification from ultrasound images using deep learning architectures, BMC Med. Imaging, 19 (2019), 51. |
[119] |
D. Ribli, A. Horváth, Z. Unger, P. Pollner, I. Csabai, Detecting and classifying lesions in mammograms with deep learning, Sci. Rep., 8 (2018), 4165. doi: 10.1038/s41598-018-29025-1
![]() |
[120] | E. M. Rutter, J. H. Lagergren, K. B. Flores, Automated object tracing for biomedical image segmentation using a deep convolutional neural network, in Medical Image Computing and Computer Assisted Intervention (eds. A. F. Frangi, J. A. Schnabel, C. Davatzikos, C. Alberola López & G. Fichtinger), Springer International Publishing, Cham, (2018), 686-694. |
[121] | S. Bakas, M. Reyes, A. Jakab, S. Bauer, M. Rempfler, A. Crimi, et al., Identifying the best machine learning algorithms for brain tumor segmentation, progression assessment, and overall survival prediction in the BRATS challenge, preprint, arXiv: 1811.02629. |
[122] | M. Havaei, A. Davy, D. Warde-Farley, A. Biard, A. Courville, Y. Bengio, et al., Brain tumor segmentation with deep neural networks, Med. Image Anal., 35 (2017), 18-31. |
[123] |
Q. Xie, K. Faust, R. Van Ommeren, A. Sheikh, U. Djuric, P. Diamandis, Deep learning for image analysis: Personalizing medicine closer to the point of care, Crit. Rev. Clin. Lab. Sci., 56 (2019), 61-73. doi: 10.1080/10408363.2018.1536111
![]() |
[124] | X. Feng, N. Tustison, C. Meyer, Brain tumor segmentation using an ensemble of 3d U-nets and overall survival prediction using radiomic features, in Brainlesion: Glioma, Multiple Sclerosis, Stroke and Traumatic Brain Injuries (eds. A. Crimi, S. Bakas, H. Kuijf, F. Keyvan, M. Reyes, T. van Walsum), Springer International Publishing, Cham, (2019), 279-288. |
[125] | E. Puybareau, G. Tochon, J. Chazalon, J. Fabrizio, Segmentation of gliomas and prediction of patient overall survival: A simple and fast procedure, in Brainlesion: Glioma, Multiple Sclerosis, Stroke and Traumatic Brain Injuries (eds. A. Crimi, S. Bakas, H. Kuijf, F. Keyvan, M. Reyes & T. van Walsum), Springer International Publishing, Cham, (2019), 199-209. |
[126] | L. Sun, S. Zhang, L. Luo, Tumor segmentation and survival prediction in glioma with deep learning, in Brainlesion: Glioma, Multiple Sclerosis, Stroke and Traumatic Brain Injuries (eds. A. Crimi, S. Bakas, H. Kuijf, F. Keyvan, M. Reyes & T. van Walsum), Springer International Publishing, Cham, (2019), 83-93. |
[127] | B. Lyu, A. Haque, Deep learning based tumor type classification using gene expression data, in Proceedings of the 2018 ACM International Conference on Bioinformatics, Computational Biology, and Health Informatics, ACM, New York, (2018), 89-96. |
[128] | P. Mobadersany, S. Yousefi, M. Amgad, D. A. Gutman, J. S. Barnholtz-Sloan, J. E. V. Vega, et al., Predicting cancer outcomes from histology and genomics using convolutional networks, Proc. Natl. Acad. Sci., 115 (2018), E2970-E2979. |
[129] | J. Zou, M. Huss, A. Abid, P. Mohammadi, A. Torkamani, A. Telenti, A primer on deep learning in genomics, Nat. Genet., 12-18. |
[130] | B. Shickel, P. J. Tighe, A. Bihorac & P. Rashidi, Deep EHR: A survey of recent advances in deep learning techniques for electronic health record (EHR) analysis, IEEE J. Biomed. Health Inf., 22 (2017), 1589-1604. |
[131] | R. A. Taylor, J. R. Pare, A. K. Venkatesh, H. Mowafi, E. R. Melnick, W. Fleischman, et al., Prediction of in-hospital mortality in emergency department patients with sepsis: a local big data-driven, machine learning approach, Academic Emerg. Med., 23 (2016), 269-278. |
[132] | T. Zheng, W. Xie, L. Xu, X. He, Y. Zhang, M. You, et al., A machine learning-based framework to identify type 2 diabetes through electronic health records, Int. J. Med. Inf., 97 (2017), 120-127. |
[133] | H. Xu, Z. Fu, A. Shah, Y. Chen, N. B. Peterson, Q. Chen, et al., Extracting and integrating data from entire electronic health records for detecting colorectal cancer cases, in AMIA Annual Symposium Proceedings, American Medical Informatics Association, (2011), 1564. |
[134] |
D. Zhao, C. Weng, Combining pubmed knowledge and ehr data to develop a weighted bayesian network for pancreatic cancer prediction, J. Biomed. Inf., 44 (2011), 859-868. doi: 10.1016/j.jbi.2011.05.004
![]() |
[135] |
R. Caruana, Multitask learning, Mach. Learn., 28 (1997), 41-75. doi: 10.1023/A:1007379606734
![]() |
[136] | S. Ruder, An overview of multi-task learning in deep neural networks, preprint, arXiv: 1706.05098. |
[137] | N. Jaques, S. Taylor, E. Nosakhare, A. Sano, R. Picard, Multi-task learning for predicting health, stress, and happiness, in NIPS Workshop on Machine Learning for Healthcare, 2016. |
[138] | R. C. Rockne, A. Hawkins-Daarud, K. R. Swanson, J. P. Sluka, J. A. Glazier, P. Macklin, et al., The 2019 mathematical oncology roadmap, Phys. Biol., 16 (2019), 041005. |
[139] |
M. Raissi, P. Perdikaris, G. E. Karniadakis, Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations, J. Comput. Phys., 378 (2019), 686-707. doi: 10.1016/j.jcp.2018.10.045
![]() |
[140] |
S. H. Rudy, S. L. Brunton, J. L. Proctor, J. N. Kutz, Data-driven discovery of partial differential equations, Sci. Adv., 3 (2017), e1602614. doi: 10.1126/sciadv.1602614
![]() |
[141] | J. H. Lagregren, J. T. Nardini, G. M. Lavigne, E. M. Rutter, K. B. Flores, Learning partial differential equations for biological transport models from noisy spatiotemporal data, preprint, arXiv: 1902.04733. |
[142] |
F. Hamilton, A. L. Lloyd, K. B. Flores, Hybrid modeling and prediction of dynamical systems, PLoS Comput. Biol., 13 (2017), e1005655. doi: 10.1371/journal.pcbi.1005655
![]() |
[143] |
J. Lagergren, A. Reeder, F. Hamilton, R. C. Smith, K. B. Flores, Forecasting and uncertainty quantification using a hybrid of mechanistic and non-mechanistic models for an age-structured population model, Bull. Math. Biol., 80 (2018), 1578-1595. doi: 10.1007/s11538-018-0421-7
![]() |
[144] | N. Gaw, A. Hawkins-Daarud, L. S. Hu, H. Yoon, L. Wang, Y. Xu, et al., Integration of machine learning and mechanistic models accurately predicts variation in cell density of glioblastoma using multiparametric MRI, Sci. Rep., 9 (2019), 10063. |
[145] | D. P. Kingma, M. Welling, Auto-encoding variational bayes, preprint,, arXiv: 1312.6114. |
[146] | C. Doersch, Tutorial on variational autoencoders, preprint, arXiv: 1606.05908. |
[147] | A. Myronenko, 3d mri brain tumor segmentation using autoencoder regularization, in Brainlesion: Glioma, Multiple Sclerosis, Stroke and Traumatic Brain Injuries (eds. A. Crimi, S. Bakas, H. Kuijf, F. Keyvan, M. Reyes, T. van Walsum), Springer International Publishing, Cham, (2019), 311-320. |
[148] | J. Xu, L. Xiang, Q. Liu, H. Gilmore, J. Wu, J. Tang, et al., Stacked sparse autoencoder (ssae) for nuclei detection on breast cancer histopathology images, IEEE Trans. Med. Imaging, 35 (2015), 119-130. |
[149] | I. Goodfellow, J. Pouget-Abadie, M. Mirza, B. Xu, D. Warde-Farley, S. Ozair, et al., Generative adversarial nets, in Advances in Neural Information Processing Systems (eds. Z. Ghahramani, M. Welling, C. Cortes, N. D. Lawrence, K. Q. Weinberger), Curran Associates, (2014), 2672-2680. |
[150] | L. Metz, B. Poole, D. Pfau, J. Sohl-Dickstein, Unrolled generative adversarial networks, preprint, arXiv: 1611.02163. |
[151] | M. Arjovsky, L. Bottou, Towards principled methods for training generative adversarial networks, preprint, arXiv: 1701.04862. |
[152] | A. Radford, L. Metz, S. Chintala, Unsupervised representation learning with deep convolutional generative adversarial networks, preprint, arXiv: 1511.06434. |
[153] |
M. Frid-Adar, I. Diamant, E. Klang, M. Amitai, J. Goldberger, H. Greenspan, GANbased synthetic medical image augmentation for increased CNN performance in liver lesion classification, Neurocomputing, 321 (2018), 321-331. doi: 10.1016/j.neucom.2018.09.013
![]() |
[154] | H.-C. Shin, N. A. Tenenholtz, J. K. Rogers, C. G. Schwarz, M. L. Senjem, J. L. Gunter, et al., Medical image synthesis for data augmentation and anonymization using generative adversarial networks, in Simulation and Synthesis in Medical Imaging (eds. A. Gooya, O. Goksel, I. Oguz, N. Burgos), Springer International Publishing, Cham, (2018), 1-11. |
[155] | J.-Y. Zhu, T. Park, P. Isola, A. A. Efros, Unpaired image-to-image translation using cycleconsistent adversarial networks, in The IEEE International Conference on Computer Vision (ICCV), 2017. |
[156] | Y. Zhang, S. Miao, T. Mansi, R. Liao, Task driven generative modeling for unsupervised domain adaptation: Application to x-ray image segmentation, in Medical Image Computing and Computer Assisted Intervention (eds. A. F. Frangi, J. A. Schnabel, C. Davatzikos, C. Alberola-López, G. Fichtinger), Springer International Publishing, Cham, (2018), 599-607. |
[157] | C. E. Gast, A. D. Silk, L. Zarour, L. Riegler, J. G. Burkhart, K. T. Gustafson, et al., Cell fusion potentiates tumor heterogeneity and reveals circulating hybrid cells that correlate with stage and survival, Sci. Adv., 4 (2018), eaat7828. |
[158] | K. Mitamura, Y. Yamamoto, N. Kudomi, Y. Maeda, T. Norikane, K. Miyake, et al., Intratumoral heterogeneity of 18f-FLT uptake predicts proliferation and survival in patients with newly diagnosed gliomas, Ann. Nucl. Med., 31 (2017), 46-52. |
[159] |
D. A. Hormuth II, J. A. Weis, S. L. Barnes, M. I. Miga, V. Quaranta, T. E. Yankeelov, Biophysical modeling of in vivo glioma response after whole-brain radiation therapy in a murine model of brain cancer, Int. J. Radiat. Oncol. Biol. Phys., 100 (2018), 1270-1279. doi: 10.1016/j.ijrobp.2017.12.004
![]() |
[160] | E. M. Rutter, T. L. Stepien, B. J. Anderies, J. D. Plasencia, E. C. Woolf, A. C. Scheck, et al., Mathematical analysis of glioma growth in a murine model, Sci. Rep., 7 (2017), 2508. |
[161] | S. Benzekry, C. Lamont, J. Weremowicz, A. Beheshti, L. Hlatky, P. Hahnfeldt, Tumor growth kinetics of subcutaneously implanted Lewis Lung carcinoma cells, 2019. Available from: https://zenodo.org/record/3572401#.Xqt0TlNKjRZ. |
[162] | M. Mastri, A. Tracz, J. M. L. Ebos, Tumor growth kinetics of human LM2-4LUC+ triple negative breast carcinoma cells, 2019. Available from: https://zenodo.org/record/3574531#.Xqtzv1NKjRZ. |
[163] | A. Rodallec, S. Giocometti, J. Ciccolini, R. Fanciullino, Tumor growth kinetics of human MDA-MB-231 cells transfected with dTomato lentivirus, 2019. Available from: https://zenodo.org/record/3593919#.Xqt0PlNKjRZ. |
1. | Luan Américo-Da-Silva, Jheimmy Diaz, Mario Bustamante, Georthan Mancilla, Ingrid Oyarzún, Hugo E. Verdejo, Clara Quiroga, A new role for HERPUD1 and ERAD activation in osteoblast differentiation and mineralization, 2018, 32, 0892-6638, 4681, 10.1096/fj.201701229RR | |
2. | Vera G. Volpi, Thierry Touvier, Maurizio D'Antonio, Endoplasmic Reticulum Protein Quality Control Failure in Myelin Disorders, 2017, 9, 1662-5099, 10.3389/fnmol.2016.00162 | |
3. | Timothy J. Bergmann, Fiorenza Fumagalli, Marisa Loi, Maurizio Molinari, Role of SEC62 in ER maintenance: A link with ER stress tolerance in SEC62-overexpressing tumors?, 2017, 4, 2372-3556, e1264351, 10.1080/23723556.2016.1264351 | |
4. | Daria Sicari, Aeid Igbaria, Eric Chevet, Control of Protein Homeostasis in the Early Secretory Pathway: Current Status and Challenges, 2019, 8, 2073-4409, 1347, 10.3390/cells8111347 |