Abstract
The tumor microenvironment (TME) presents a challenging barrier for effective nanotherapy-mediated drug delivery to solid tumors. In particular for tumors less vascularized than the surrounding normal tissue, as in liver metastases, the structure of the organ itself conjures with cancer-specific behavior to impair drug transport and uptake by cancer cells. Cells and elements in the TME of hypovascularized tumors play a key role in the process of delivery and retention of anti-cancer therapeutics by nanocarriers. This brief review describes the drug transport challenges and how they are being addressed with advanced in vitro 3D tissue models as well as with in silico mathematical modeling. This modeling complements network-oriented techniques, which seek to interpret intra-cellular relevant pathways and signal transduction within cells and with their surrounding microenvironment. With a concerted effort integrating experimental observations with computational analyses spanning from the molecular- to the tissue-scale, the goal of effective nanotherapy customized to patient tumor-specific conditions may be finally realized.
Challenges of the Tumor Microenvironment to Drug Delivery
The Tumor Microenvironment Inhibits Drug Delivery
In order for drug molecules to elicit a pharmacological response, the molecules must arrive in sufficient quantities to the tissue of interest and bind to the molecular target activating or inhibiting particular pathways. To achieve therapeutic responses in solid tumors, drug molecules need to overcome various barriers at different physical scales. The tumor microenvironment (TME) includes several scales: (a) molecular (nano-) scale, including up- and down-regulation of various proteins that can signal for tumor growth or drug-efflux mechanisms; (b) nano- to micro- scale, which incorporates gradients of cell nutrients and oxygen, growth factors, and other means of cell-to-cell communication; (c) micro-scale, in which interactions between cells occur in the acellular stroma compartment of the tumor; (d) micro to macro scale, which incorporates the organ architecture, blood supply, lymphatics, and other physiological factors. While these barriers span several orders of magnitude, they are intricately linked and cross-communicate. As an example, the architecturally/anatomically irregular and functionally impaired tumor neovasculature (micro to macro scale) is characterized by reduced oxygen tension, oscillating flow, constricted blood vessels, and other abnormal features (Jain, 2003; Folkman, 2007; Chung et al., 2010; Dewhirst and Secomb, 2017). Consequently, the TME becomes heterogeneous in terms of gradients of solutes and nutrients (nano- to micro- scale) as well as differences in pH and cell viability due to hypoxia (micro-scale). Hypoxia promotes recruitment of immune cells to the tissue, while prompting release of cytokines and chemokines (molecular (nano-) scale) that affect cell-to-cell interactions (micro-scale).
The heterogeneous TME has a significant effect on therapeutic outcomes. First, TME heterogeneity and three-dimensionality represent a significant barrier to systemically administered therapeutics, including nanotherapeutics. As a result, in vitro efficiencies of anti-cancer drugs (especially those shown in 2D cultures) do not correlate well with potencies observed in vivo, as has been shown in several studies (Agiostratidou et al., 2001; Fruehauf, 2002). The results highlight discrepancies in positive predictive values between clinical efficacies and in vitro therapy selection (50–70%) and negative predictive accuracy (∼90%), demonstrating that enhanced potency of drugs in 2D cultures largely disregards the barriers in the heterogeneous TME. Second, cells in the TME such as endothelial cells, macrophages and other cells of the immune system, and fibroblasts/myofibroblasts actively interact with the tumor cells in most solid tumors and affect cancer cell proliferation, survival, polarity and invasive capacity (Williams et al., 2016; Guo and Deng, 2018; Sarode et al., 2020).
Multiple studies have demonstrated that while the normal cellular microenvironment can inhibit or even prevent the growth of tumor cells, the changes that happen in the TME synergistically support tumor growth. Tumors shape their microenvironment promoting the growth not only malignant cells, but also non-malignant TME or stromal cells. There are many mechanisms that still need to be elucidated in the tumor-stroma interactions, although the importance of an altered TME in the process of tumorigenesis is no longer questioned. Numerous successful cancer therapies targeting the TME have been approved or are being developed, highlighting the importance of the cells and tissue neighboring malignant cells for tumor survival and invasion. As an example, tumor macrophages can be polarized to be tumor-growth supportive (M2 phenotype) or inhibiting (M1 phenotype) (Pollard, 2008; Ruffell et al., 2012; Ruffell and Coussens, 2015; Williams et al., 2016; Kielbassa et al., 2019). The M1/M2 ratio has been shown to be a strong prognostic factor in a variety of solid tumors including liver metastasis (Cui et al., 2013; He et al., 2013; Zhang et al., 2014; Yuan et al., 2017).
The Organ Microenvironment Promotes Metastasis Development
The TME is dependent not only on the origin and the characteristics of tumor cells, but also on the anatomy and physiology of the organ to which tumor cells disseminate. We discuss the liver as one particular example to illustrate this complexity (van den Eynden et al., 2013) and the challenges it poses to drug delivery. Metastatic lesions represent the most common malignancy in the liver and are up to 40 times more frequent in clinical practice than primary liver tumors (Rummeny and Marchal, 1997). The liver is a highly vascularized organ that has a dense network of capillaries, sinusoids, efficiently providing oxygen and soluble nutrients to the innermost cells in the organ. Two physiological factors have been linked to the high incidence of the liver being an organ of choice for distant metastasis: (Chung et al., 2010) increased likelihood of invasion due to dual blood supply from systemic and portal circulation; (Dewhirst and Secomb, 2017) presence of fenestrations in the liver sinusoids that allow for tumor cell invasion from the circulation.
Uniquely, incipient liver metastases preserve the stromal structure of the liver and do not rely on angiogenesis for survival (Stessels et al., 2004). This vascularization pattern, in which tumor cells primarily use existing vasculature in surrounding parenchyma, is unconventional, compared to most solid tumors (Pezzella et al., 1997; Danet et al., 2003a, b; Namasivayam et al., 2007), and significantly limits diffusive transport into the lesions Thus, the most frequent tumor types with liver metastasis, including breast, colon, lung, and gastric carcinomas create hypovascular lesions usually showing perilesional enhancement (Namasivayam et al., 2007; Hazhirkarzar et al., 2020). This characteristic significantly impairs the delivery of systemically administered therapeutics to tumors in the liver.
For instance, it has been shown that breast cancer and pancreatic ductal adenocarcinoma liver metastases are characterized by poor permeation of molecules and are clinically observed as hypo-attenuating spots, which intravenously injected contrast agents do not permeate (Liu et al., 2003). Impaired diffusion is an important factor limiting adequate concentration of therapeutics, and could explain why chemotherapy fails to cure unresectable liver lesions. Poor permeation is especially acute with high molecular weight (HMW) molecules, such as m99Tc microaggregated albumin (Daly et al., 1985).
Experimental Evidence With Hypo-Perfused Tumors
Recent data from in vivo studies in breast cancer liver tumors with low vascularization patterns and high macrophage (Mϕ) content confirm that these lesions lack efficient perfusion. In intravital microscopy (IVM) studies, transport of HMW molecules (fluorescent dextrans) into tumor lesions was impeded compared to healthy tissue (Figure 1). When a HMW drug nAb-PTX, possessing similar hydrodynamic diameter as 40 KDa dextran (∼12 nm), was packaged in a solid multistage porous silicon NV (MSV) (Godin et al., 2010a, b, 2011, 2012; Tanaka et al., 2010; Tasciotti et al., 2011; Srinivasan et al., 2013; Yokoi et al., 2013; Jaganathan et al., 2014; Tanei et al., 2016), the nanotherapeutic was taken up by Mϕ in the proximity of breast cancer liver metastases, thus enabling higher drug concentrations in the lesions and dramatically improving nAb-PTX therapeutic efficacy and animal survival (Figure 2).
FIGURE 1
FIGURE 2
Evidence with other types of hypo-perfused tumors supports these data. The composition of the stroma in preclinical models of orthotopic primary pancreatic based on L3.6pL cells were analyzed in Yokoi et al. (2013). Profound differences in the cellular elements of pancreatic stroma as a result of different treatments were observed (Figure 3, Borsoi et al., 2015). Mice were treated with a combination of gemcitabine and nAb-PTX encapsulated or nAb-PTX packed into an engineered MSV (Tasciotti et al., 2011; Yokoi et al., 2013; Borsoi et al., 2015, 2017). These data demonstrate that cellular elements in the TME of poorly perfused tumors dramatically change during disease progression and in response to therapy.
FIGURE 3
Three-Dimensional in vitro and ex vivo Models of the Cancer Microenvironment
Tools and Considerations for Preclinical Models to Evaluate Drug Delivery and Efficacy
In vitro approaches have been a traditional stronghold for anti-cancer therapeutic screening, with the goal of maximizing predictive potency while closely mimicking tumor physiology and cell-cell interactions (Zoli et al., 2001; Breslin and O’Driscoll, 2016). Importantly, 3D in vitro models complement in vivo rodent models and overcome their limitations, such as high cost, long latency and ethical minimization of the use of animals. Furthermore, microenvironments that recapitulate more precisely human tumor physiology can be created in vitro, enabling evaluation of the effect of various factors on the mechanisms of tumor initiation, progression and response to therapy. The addition of the third spatial dimension to in vitro models has introduced improved cell-cell interaction setting up transport limitations to oxygen, nutrients and potential therapeutics, thereby increasing the utility for therapeutic screening (Griffith and Swartz, 2006). Here, we briefly review considerations in evaluating drug efficacy in terms of adopting the appropriate 3D in vitro/ex vivo models, and consider integrated engineering examples of therapeutic efficiency directed toward hypovascularized lesions such as liver metastases.
A majority of therapeutics suffer from high attrition rates as they move through the oncologic discovery pipeline, due to a lack of translation from the preclinical to the clinical stage (Hay et al., 2014; Wong et al., 2019). Preliminary considerations while identifying the most appropriate in vitro model draw a balance between the complex physiological relevance (cellular complexity, molecular pathology) of 3D models and their ability to scale up for drug screening. While cellular simplicity may be sufficient for evaluating key molecular targets for targeted therapy, models that include elements of tumor stroma have to be considered for more nuanced drug delivery approaches. For example, a nanoparticle-based approach targeting collagenase in hypovascular pancreatic ductal adenocarcinoma had the goal of demonstrating increased drug penetration after nanoparticle treatment (Zinger et al., 2019). Similarly, novel nano-immunotherapies in development target phagocytic cells like macrophages or antigen presenting cells within the TME (Caster et al., 2019).
Another important consideration in preclinical screening is how therapeutic efficacy is defined. Specific to 3D models, IC50 alone may not be the most suitable metric to evaluate therapeutic efficiency in isolation (Zoli et al., 2001). Traditional plate-reader based viability/metabolic activity assays are not designed to discern individual contributions of cellular compartments within tumoroid/organoid models (Horman et al., 2015). However, a combination of IC50 with imaging methodologies that evaluate drug uptake, molecular read outs and morphometric analysis may offer more information of the dynamics of the therapeutic process (Raghavan et al., 2015, 2016; Stock et al., 2016; Santo et al., 2017). Such high content imaging approaches can reliably provide morphological and cellular dynamics information that can be correlated with therapeutic outcomes (Garvey et al., 2016; Ahonen et al., 2017; Bulin et al., 2017).
Integrated Engineered Tumor Model Approaches for Anti-cancer Therapeutic Screening
A fundamental consideration in preclinical in vitro models remains the methodologies utilized to create the actual models – engineering approaches can often dictate resultant drug sensitivity. This section briefly describes some of the most popular approaches utilized currently to manufacture 3D in vitro/ex vivo tumor models, with examples of how they have been utilized in drug delivery to hypovascular lesions such as liver metastases and pancreatic ductal adenocarcinoma.
Spheroids and Organoids
Spheroids represent the most frequently used 3D cancer model in preclinical research. Most cancer cell lines will grow and organize as avascular 3D structures called spheroids, either aided by scaffolds, hydrogels, magnetic levitations, extracellular matrix (ECM) gels or in scaffold-free suspension or liquid overlay cultures. The organoid/tumoroid model has been demonstrated as a promising tool to recapitulate patient response in colorectal (Sato et al., 2011), ovarian (Kopper et al., 2019), pancreatic (Boj et al., 2015; Ware et al., 2016a, b), breast (Jaganathan et al., 2014; Leonard and Godin, 2016), and gastrointestinal cancers (Aberle et al., 2018). Spheroids and organoids derived from cancer stem cells/tumor-initiating cells have been described to capture inter-patient heterogeneity, with the opportunity to rapidly relate phenotype to genotype, and high fidelity drug sensitivity (Sachs and Clevers, 2014; Raghavan et al., 2017). In regards to 3D models of hypovascularized tumors (e.g., liver metastasis of breast, lung and colorectal tumors, pancreatic, etc.) tumor spheroids recapitulate the nutrient and therapeutic supply patterns (mainly from the surrounding tissue).
The biggest advantage of the spheroid/organoid approach is to be able to deconstruct the cellular and environmental elements of the TME, and reverse engineer them as appropriate, to study specific targeting strategies. For example, avascular spheroids of breast cancer cells surrounded by macrophages were constructed using magnetic levitation and bioprinting to recapitulate the microenvironment of liver metastasis (Leonard and Godin, 2016; Leonard et al., 2016, 2020). In another study, hepatocellular carcinoma (HCC) spheroids were created by combining alginate microbead technology with decellularized liver matrix-derived ECM, to encapsulate HCC cells (Sun et al., 2018). To model the anchorage-independent aggregation of pre-metastatic ovarian cancer cells within malignant ascites, spheroids containing ovarian cancer stem cells and macrophages were utilized (Raghavan et al., 2019). Pancreatic tumor-rich spheroids were designed for drug screening with various human pancreatic cancer cells and human stellate cells using a combination of polymer matrix and hanging drop techniques (Ware et al., 2016b).
Ex vivo Explant Cultures/Sheet Models
The establishment of these explant cultures/sheet models relies on precision-cut slicing mechanisms, and access to primary patient-derived tumors. These models have the advantage of preserving complete tissue architecture, as well as, the cellular heterogeneity of the tumor, and can virtually be derived from any accessible surgical tumor specimen. The challenges of preservation and proliferation, however, still remain. Using explant cultures, anticancer drug screening has been demonstrated in rectal cancer liver metastases, which could be replicated in rodent xenograft models (Zhang et al., 2020). Tumor slices have also been utilized to test nanoparticle-based drug delivery of oligonucleotides to ovarian cancer and gliomas, and non-small cell lung cancer (Dong et al., 2011; Ewe et al., 2017).
Decellularized Liver Scaffolds and Matrices
The use of naturally derived scaffolds for solid organs has been well developed for regenerative medicine (Baptista et al., 2009), and has made inroads into use in anticancer therapeutics (Hinderer et al., 2016). The most advantageous aspect of the use of natural ECM-based bio-scaffolds is the ability to retain bioactive molecules like growth factors and other signaling compounds embedded within the scaffolding structure, providing cues to tumor cells that are seeded on them. Metastatic colorectal cancer cells and native HCC cells grown on human decellularized liver scaffolds demonstrated less efficacy to chemotherapeutic regimens that were used at standard 2D determined IC50 (Hussein et al., 2016; D’Angelo et al., 2020). This implies the retaining of bioactivity in decellularized scaffolds, which speaks not only to their utility in re-engineering the architectural and scaffolding components of the TME, but also highlights that the metastatic TME triggers a physiologically more resistant cancer phenotype. Recently, decellularized liver scaffolds have been manufactured by immersion decellularization of chick embryos, which bypasses the slow whole-organ decellularization approach. A dramatic reduction in doxorubicin efficiency against triple negative breast cancer liver metastases was demonstrated (Guller et al., 2020). Importantly, the speed at which the matrices can be produced using this approach lends itself high throughput amenable for anticancer screening. It should be noted that decellularized tissue scaffolds can be combined with other techniques for 3D tumor growth, as described above.
Microfluidic Organ-on-a-Chip
Microfluidic approaches have been implemented to generate tumor spheroids, relying on perfusion and microwells. These models are mainly used for tumors with enhanced degree of angiogenesis and vascularization as compared to the surrounding tissue. The utility of microfluidic approaches in the study of cancer drug delivery is well reviewed (Unger et al., 2014; Ozcelikkale et al., 2017). The advantage of using microfluidic approaches is the inclusion of dynamics of the TME (Agastin et al., 2011; Kim et al., 2012). Microfluidic models of pancreatic ductal adenocarcinoma demonstrated transcriptome level similarity to patient-derived pancreatic stellate cells (Gioeli et al., 2019). To model liver metastasis, gelatin methacryloyl was combined with decellularized liver matrix components including bioactive factors derived from liver matrix, in a dynamic culture system to assess drug-dose responses to acetaminophen and sorafenib (Lu et al., 2018). In a similar metastasis-on-a-chip model of kidney cancer metastasized to the liver, NP-conjugated 5-fluorouracil was shown to be more efficacious than free drug (Wang et al., 2020).
3D Bioprinting
Bioprinting is an emerging engineering strategy to develop 3D tumor models, where 3D printing is exploited to deposit cells and biomaterials in tissue-like structures. The biggest limitation of bioprinting approaches has been the maintenance of cellular viability through the simultaneous layer-by-layer assembly process of cells and bio-inks (Albritton and Miller, 2017). Early attempts using bio-inks like gelatin, alginate and/or fibrinogen demonstrated that heterotypic cell interactions can be maintained in bioprinted tumors, which result in them being more chemoresistant (Zhao et al., 2014; Dai et al., 2016; Zhang et al., 2016; Zhou et al., 2016). This technique has been interestingly utilized to create metastasis models at the tumor-vasculature interface (Meng et al., 2019), as well as, mimicking organ specific metastasis using decellularized liver matrix ECM as a bioink (Lee et al., 2017). This technology has potential for combining organoid technology with rapid production, potentially rendering it high-throughput amenable for personalized medicine and therapeutic screening applications.
Drug-Based Nanotherapy in Cancer
Nanotherapy Can Be More Effective Compared to Free Drug Infusion
Nanoparticles have been used clinically for tumor therapy since the early 1990s. Currently available nanotherapies improve safety and efficacy of chemotherapies. As an example, Doxyl®, the first nano-drug, liposomal doxorubicin, was introduced to reduce the cardiotoxicity of doxorubicin, leveraging the differences in biodistribution of free drug vs. liposomal entity (Barenholz, 2012). Since then, a few dozens of nanotherapeutics have been approved for cancer therapy. Many other nanomedicine-based approaches are under investigation worldwide. Nanotherapy offers the possibility to improve metastatic disease treatment by increasing the concentrations and attaining controlled release of therapeutics in distant lesions (Sanga et al., 2007; Farokhzad Omid and Langer, 2009; Misra et al., 2010; Bourzac, 2012; Parhi et al., 2012; Zamboni et al., 2012; Bertrand et al., 2014; Grodzinski and Farrell, 2014; Fernandes et al., 2015). Nanovectors (NV) can favorably change pharmacokinetics of drugs in plasma and in tissues, prolonging circulation time and enhancing delivery to tumors (Bae et al., 2005, 2007; Bae and Kataoka, 2009; Maniwa et al., 2010; Ponta and Bae, 2014; Khawar et al., 2015; Reichel et al., 2015, 2017; Curtis et al., 2016b; Khalid et al., 2016).
Generally, systemically administered nanotherapies first flow in the circulation and either attach to the elements in the TME (e.g., receptors overexpressed in the tumor-associated endothelium and stroma elements) or extravasate through the gaps in the leaky and disorganized tumor neovasculature. A key mechanism related to the extravasation of nanotherapies and macromolecules in the tumor is called the Enhanced Permeation and Retention (EPR) effect (Greish, 2010). Factors that contribute to the EPR effect include: enhanced vascular permeability and angiogenesis, which sustain rapid growth of tumor on one side, and impaired lymphatic drainage, on the other side (Maeda et al., 2003; Nakamura et al., 2015).
Another mechanism for tumor targeting proposed specifically for nanotherapies includes their interaction with immune cells, such as macrophages, which are abundant in the TME. Macrophages are professional phagocytes and, as such, efficiently take up particles from the circulation (Gustafson et al., 2015). These immune cells can thus serve as a cellular depot for therapeutics in the TME (Choi et al., 2007; Leonard et al., 2016; Tanei et al., 2016; Hu et al., 2019). Additionally, nanotherapies can affect the polarization and the function of the immune cells, reversing their phenotype from pro-tumorigenic M2 to anti-tumorigenic M1 (Ponzoni et al., 2018; Hu et al., 2019; Reichel et al., 2019).
Challenges for Nanotherapy in the Tumor Microenvironment
Clinical benefits from nanotherapy remain controversial because in vivo efficacy varies from tumor to tumor (Kyle et al., 2007; Vicent et al., 2009; Cukierman and Khan, 2010; Fang et al., 2011; Bhatia et al., 2012; Mehta et al., 2012; Cho et al., 2013; Crist et al., 2013; Prabhakar et al., 2013; Venditto and Szoka, 2013; Dawidczyk et al., 2014; Danhier, 2016). Further, it is well known that organ physiology and microenvironment significantly impact the efficacy of cancer therapy (Curtis and Frieboes, 2016). As an example, while well-documented in primary tumors, the EPR effect does not usually apply to organs enriched with blood vessels, such as the liver. In the liver, tumor lesions are less vascularized than the organ itself (e.g., they appear as “white spots” on a “red bed”). Thus, management of cancer metastasized to different organs should ideally account for these variations, providing personalized therapy based on the location of the metastasis.
A number of HMW-based therapeutic strategies has been clinically approved and proposed for the therapy of advanced breast cancer, pancreatic ductal denocarcinoma and other tumors, including albumin-bound drug conjugates, such as nanoalbumin bound paclitaxel, nAb-PTX or Abraxane® (Blomstrand et al., 2019; Li and Kwon, 2019; Untch et al., 2019), monoclonal antibodies (mAb), such as anti-HER2 mAb trastuzumab (Figueroa-Magalhaes et al., 2014) or anti-EGFR mAb, cetuximab (Huang and Buchsbaum, 2009), and genetic materials, including siRNA, mRNA, and aptamers (Kang et al., 2015; Ngamcherdtrakul and Yantasee, 2019). In liver metastatic lesions with low vascularization patterns, these potent therapeutics are unable to be transported deeply enough into the lesions prior to their clearance from circulation. Thus, new approaches to enhance their accumulation in liver metastases are necessary.
Mathematical Modeling of Drug-Based Cancer Nanotherapy
The “One Drug” (a.k.a. “Silver Bullet”) Approach to Cancer Drug Therapy Has Failed
The complexity of the TME coupled with cellular plasticity may preclude any one particular therapeutic from fully succeeding in eradicating a tumor. Yet, historically, the search for the one “silver bullet” that could cure patients has garnered much attention in research and in clinical medicine (e.g., Zempolich and Dubeau, 1998; Koontz, 2017). More recently, approaches that consider information from different physical scales and perspectives, and that are able to integrate the associated information into a coherent picture, have been pursued. These approaches include the combination of experimentation with mathematical modeling and computational simulation (Frieboes et al., 2006, 2012, 2019; Sanga et al., 2007; Sinek et al., 2007; Brocato et al., 2018; Dogra et al., 2019).
A Systems-Level Approach Is Required for Nanotherapy to Succeed
Newer approaches to cancer therapy involving molecular profiling represent a promising avenue, especially as they seek to elucidate the role of the microenvironment in the evolution of acquired drug resistance. In particular, the development of gene analysis tools has offered the opportunity to more quickly and cheaply assess variations in the genetic make-up of tumors (Heath et al., 2016; Miryala et al., 2018). Yet the interpretation of these data for clinical application is non-trivial. Although the detection of genetic variation in individual patient tumors is considered crucial for the success of personalized medicine, for the most part it remains unclear how this variation translates to tumor-scale phenotype. With the exception of a few well-studied genes [e.g., BRCA1 and BRCA2 in breast cancer (Grimmett et al., 2018)], most of the genetic information elucidated from these analyses has yet to be meaningfully interpreted.
Mathematical Modeling as a Systems-Level Approach
A major reason for the challenge to link the molecular to the tissue scale is that the growth and treatment response of many cancers do not solely depend on variation at the genetic scale but rather on the combination of characteristics at multiple scales, including genetic, cellular, and tissue conditions such as the tumor and organ microenvironment. A purely empirical approach would be insurmountable to determine optimal drug therapy, due to the many variables involved. For this task, which requires a systems-level perspective, mathematical modeling and computational simulation are ideally suited. Network-oriented approaches have sought to “connect the dots,” so to speak, to make sense of intra-cellular relevant pathways as well as signal transduction within cells and with their surrounding microenvironment (Kreeger and Lauffenburger, 2010; Bachmann et al., 2012; Lecca and Re, 2017). In addition to the effort to model the molecular-scale, the integration of experimental and computational modeling to enable realistic, predictive evaluation of tumor behavior during therapy has been pursued. The application of mathematical modeling can help to bridge the associated physical scales (from molecules to tissue), and thus lead to improved interpretation of particular tumor characteristics and how they might influence the drug transport in the tumor, and consequently, the drug response.
Mathematical Modeling Provides a Link to the Molecular Scale
Mathematical modeling complements network-oriented techniques (e.g., principal network analysis) and data-based approaches (e.g., statistical and machine learning) by providing for mechanistic insight of tissue behavior in time as well as space. This would be considered important to evaluate the effects of the microenvironment on tumor response to therapy. In particular when evaluating omics strategies, recent work strongly suggests that the best predictive results are obtained via studies accounting for multiple types of datasets (Cho et al., 2014; Ge et al., 2015; Li et al., 2016; Ren et al., 2016; Dong et al., 2017), as a single approach may not suffice for personalized cancer patient treatment. The interdependence of the datasets indicates an integrated approach that links genes to phenotype (Su et al., 2011; Zhang et al., 2013; Bolger et al., 2014; Minton et al., 2015; Murakami et al., 2015). For example, metabolomics provides insight into outcomes of transcription changes, which reflect differential functionality of specific metabolites influencing the tumor response. Attempts to predict therapy response based on metabolomics analysis only (Tian et al., 2018) or gene expression (Geeleher and Cox, 2014) yield results typically relying on statistical analyses that do not necessarily represent any particular tumor (Peng et al., 2018; Mucaki et al., 2019). Ideally, a mathematical-based framework would be capable of recreating particular patients’ tumors for in silico prediction of behavior in time and space prior to treatment, incorporating omics data for patient customization, and thus move toward the goal of predictive personalized treatment.
The ability to predict personalized drug response would help to obviate unnecessary, high-cost and high-morbidity treatments, allowing more efficient patient management. Although a number of theoretical models of tumor drug response have been developed in recent years (e.g., Byrne et al., 2006; Ribba et al., 2006; Stamatakos et al., 2006; Enderling et al., 2007; Roeder and Glauche, 2008; Frieboes et al., 2009; Hinow et al., 2009; Sinek et al., 2009; Gevertz, 2011), few have focused on a multiscale integration of molecular data to evaluate cancer treatment response. By extracting mathematical model parameter values from tumor-specific cell proliferation, apoptosis, and molecular characteristics, and simulating the effects of nanocarrier and drug transport and retention in tissue, it may become possible to predict the drug response customized for individual patients, beyond what would have been expected from sole consideration of any one of these parameters or the intrinsic resistance of the cancer cells themselves.
Modeling of Cancer Nanotherapy Taking Into Account the Microenvironment
Mathematical Modeling Focusing on Drug Delivery
Mathematical modeling and computational analysis are actively being pursued in several aspects of oncology to personalize and improve therapeutic outcomes (e.g., Ibrahim-Hashim et al., 2017). In particular, tissue structure and transport in liver (Rani et al., 2006; Hoehme et al., 2007; Campbell et al., 2008; Hoehme et al., 2010, 2017, 2018; Holzhutter et al., 2012; Drasdo et al., 2014; Dutta-Moscato et al., 2014; Lettmann and Hardtke-Wolenski, 2014; Schliess et al., 2014; Siggers et al., 2014; Bethge et al., 2015; Ricken et al., 2015; Schwen et al., 2015; Nishii et al., 2016; Sluka et al., 2016; White et al., 2016; Friedman and Hao, 2017; Hudson et al., 2017, 2019; Meyer et al., 2017; Fu et al., 2018; Mahlbacher et al., 2018; Clendenon et al., 2019; Van Liedekerke et al., 2020) as well as pancreas (Haeno et al., 2012; Louzoun et al., 2014; Ng and Frieboes, 2017, 2018; Roy and Finley, 2017; Yamamoto et al., 2017; Chen et al., 2020; Dogra et al., 2020b) have been modeled. While numerous studies have simulated tumor growth and angiogenesis [see recent reviews and related work (Cristini et al., 2008; Edelman et al., 2010; Lowengrub et al., 2010; Osborne et al., 2010; Rejniak and McCawley, 2010; Vineis et al., 2010; Andasari et al., 2011; Chaplain, 2011; Deisboeck et al., 2011; Frieboes et al., 2011; Michor et al., 2011; Rejniak and Anderson, 2011; Swanson et al., 2011)] including metastatic conditions (Campbell et al., 2008; Haeno et al., 2012; Bethge et al., 2015; Hudson et al., 2017, 2019), as well as pharmacokinetics/pharmacodynamics (Lee et al., 2013; Pascal et al., 2013a, b; Wang et al., 2015, 2016; Barbolosi et al., 2016; Enriquez-Navas et al., 2016; Curtis et al., 2018) and drug discovery (Wang and Deisboeck, 2014; Zhang and Brusic, 2014), few have focused on HMW-based therapeutics. To address this need, models have been proposed which, coupled with experimentally measured parameters, have paved the way for more realistic multiscale modeling of cancer nanotherapy (Dogra et al., 2019), integrating spatial scales nm to cm- and temporal scales from sub-sec to weeks (as summarized in Table 1).
TABLE 1
| Nanotherapy Focus | References |
| Biodistribution/physiological pharmacokinetics | Li et al., 2010, 2012; Li and Reineke, 2011; Dogra et al., 2020a |
| Transport in avascular tumors | Gao et al., 2013; Curtis et al., 2016a |
| Transport in irregularly vascularized tumors | van de Ven et al., 2013; Wu et al., 2014; Curtis et al., 2016a; Miller and Frieboes, 2019a, b |
| Transport based on nanoparticle physical characteristics | Decuzzi and Ferrari, 2006; Decuzzi et al., 2009; Godin et al., 2010b |
| Binding to tumor vasculature | Frieboes et al., 2013; Curtis et al., 2015; Chamseddine et al., 2018, 2020; |
| Interactions with macrophages | Leonard et al., 2016, 2017, 2020; Mahlbacher et al., 2018 |
| Intracellular pharmacokinetics | Li et al., 2013; Miller and Frieboes, 2019b |
| For tumor detection | Reichel et al., 2015 |
| For hyperthermia applications | Kaddi et al., 2013 |
| For drug delivery | van de Ven et al., 2012; Li et al., 2013; Curtis et al., 2015, 2016a,b; Leonard et al., 2016, 2017, 2020; Chamseddine et al., 2018, 2020; Miller and Frieboes, 2019a, b; |
Overview of recent mathematical modeling to study nanoparticle delivery and efficacy in tumors.
Traditional pharmacokinetics and pharmacodynamics (PK/PD) correlate efficacy of small molecule drugs with time-dependent changes in average bulk drug concentration at the site of action. However, it has been shown that these PK/PD methods often fail to predict efficacy of drug-loaded NV that control drug concentrations in vascular, interstitial and intracellular spaces of tumor (Ponta and Bae, 2014; Ponta et al., 2015; Curtis et al., 2016b). In addition, traditional analyses frequently simplify or rule out tumor physiology changing over time (Khawar et al., 2015; Khalid et al., 2016). To overcome these issues, computational modeling that integrates key parameters influencing NV-based drug delivery into tumor tissue to realistically predict in vivo efficacy of nanotherapy has been proposed.
Physiologically based pharmacokinetic (PBPK) models have been utilized to evaluate drug efficacy, including nanoparticle-delivered drugs (Li et al., 2010). These models evaluate the absorption, distribution, metabolism, and excretion of small molecules such as drugs or nanoparticles by organizing tissues and organs as distinct compartments linked via mass transport. In Li et al. (2013), PBPK modeling was applied to evaluate the intracellular pharmacokinetics of paclitaxel delivered by nanoparticles, with the model parameters set from experimental data with human breast cancer MCF7 cells in vitro. The results revealed that the simulated intracellular pharmacokinetics corresponded with relevant parameters, including nanoparticle PK, drug-release kinetics, and drug dose.
The effect of irregular vascularization on nanoparticle transport and drug release in tumor tissue has been evaluated in several studies (van de Ven et al., 2012; Li et al., 2013; Curtis et al., 2015, 2016a,b; Leonard et al., 2016, 2017, 2020; Chamseddine et al., 2018, 2020; Miller and Frieboes, 2019a, b). The role of macrophages in nanotherapeutic transport and effect on hypo-vascularized tumor lesions such as breast cancer liver metastases was evaluated in Leonard et al. (2016, 2017). In particular, the response due to repetitive therapy with MSV-nAb-PTX and nAb-PTX, showed that encapsulation of the drug in multi-stage vectors (MSV-nAb-PTX) could maximize the tumor regression (Figure 4).
FIGURE 4
Case Study: Modeling of Liver Metastasis Nanotherapy Response
Shifting macrophage polarization from an anti-inflammatory and tumorigenic M2 phenotype to a pro-inflammatory and anti-cancerous M1 phenotype has recently garnered increased focus (Pyonteck et al., 2013; Leonard et al., 2017; Tariq et al., 2017; Poh and Ernst, 2018), with some promising results (Pyonteck et al., 2013). This shift would ideally be combined with standard therapies. Computational modeling has predicted that the tumor response depends non-linearly on the M1:M2 ratio (Leonard et al., 2017). To explore this further, mathematical modeling was recently employed to analyze the effects of the nanotherapy while simulating manipulation of the macrophage phenotype via a hypothetical immunotherapy affecting macrophage polarization (Leonard et al., 2020). Although the role of macrophages in cancer therapy has been evaluated in the past via mathematical modeling (Mahlbacher et al., 2019), the effect of varying macrophage phenotypes on nanotherapy response has not been extensively explored. The simulations indicated that the M2-tumor interaction may have a dual role in the response to MSV-nab-PTX, initially promoting tumor death and subsequently aiding tumor regrowth.
To test this model-derived hypothesis, CRISPR technology was employed in the laboratory to achieve a stable polarization of macrophages and avoid their repolarization in the dynamically changing TME. The experiments showed that the response to MSV-nab-PTX was non-uniform with respect to the M1:M2 ratio. For 72 h exposure, an M1:M2 ratio of 1500:500 reached lower viability than 2000:0 with only M1, demonstrating that the M2 subtype increases the therapeutic efficacy. Similarly, the ratio of 0:2000 with only M2 had lower viability than the 500:1500 ratio. The tumor response to MSV-nab-PTX loaded macrophages predicted by the computational model is in Figure 5. Figure 5A shows a general trend of decreased tumor size when the M1:M2 ratio increases. Simulating the inactivation of M2 macrophages to gauge the M1-only effect (Figure 5B) while maintaining the same number of activated M1 shows that the response is significantly less than when the M2 are active, even for a high M1:M2 ratio of 3.8:1. Hence, a dual action of the M2 macrophages is forecast by the model. Since PTX is a cell-cycle inhibitor, M2 macrophages synergistically augment the drug effect during treatment by promoting cell proliferation, and support tumor recovery after the treatment. By simulating repeated treatment cycles with MSV-nab-PTX (Figures 5C,D), the model showed that the presence of both macrophage subtypes significantly supports tumor regression.
FIGURE 5
These modeling results suggest that immunotherapy strategies primarily dependent on raising the M1:M2 ratio may be less effective than protocols that establish an M1:M2 ratio that maximizes tumor regression during chemotherapeutic exposure, and then tilts this ratio in favor of the M1 macrophages during the tumor recovery phase in order to leverage their cytotoxic effect.
Conclusion
While molecular targets in tumors are currently clinically evaluated and considered in determining therapeutic strategies, the physical and physiological barriers in the TME may not be taken into account. As discussed in this review, in many instances, the resistance to therapy can have physical or physiological origins. In the case of tumor metastasis to the liver or other hypovascularized lesions, the tumor lesion blood supply represents a critical limiting factor. The ability to retain the drug in the proximity of tumor cells, for example, by anchoring it to the cells of the TME, could bring significant therapeutic advantage. Modeling therapeutic responses and the efficiency of advanced tools, such as nanomedicines, for enhancing these responses would be of prime interest to improve outcomes.
Three-dimensional tumor models are being designed and utilized to bridge the gap between 2D cell cultures and the gold-standard animal models. While for some purposes, such as drug delivery from vasculature to the tumor mass in hypovascular tumors, simple, high-throughput and highly reproducible spheroids can be used, organoids may represent an advantageous systems when TME-tumor cell interactions are important to consider (e.g., in the case of immunotherapy evaluation). With the advent of bioprinting technologies and the ability to recreate complex heterotypic cellular interactions within 3D models, it is expected that the use of 3D tumor models in anticancer therapeutic screening will be significantly expanded. Combined with novel nanotherapy-based targeting strategies and integrated with computational predictions of NP behavior within the TME, preclinical testing in 3D tissue models could benefit from in vitro-ex vivo/in silico approaches in the oncologic drug discovery pipeline.
Nanotherapy-tumor interactions are expected to depend non-linearly (i.e., non-additively) on nanotherapy and tumor tissue-specific conditions, including vascularization, hypoxia, and other microenvironment characteristics affecting tumor response. We have illustrated in this review that the analysis of such interactions could benefit from mathematical modeling that provides a capability for system analysis. In order to leverage the power of these models, their parameters need to be based on biologically relevant data, including clinical information. Recent advances in computational power may enable simulations with enhanced biology to more fully capture the complexity of the TME, including fibroblast cells and extra-cellular matrix components. These models could then be integrated with network-oriented approaches to fully link the molecular- to the tissue-scale. Ideally, nanotherapy candidates would be evaluated prior to treatment by informing the parameters of such models with patient tumor-specific characteristics, as can be observed via analysis of biopsy samples, imaging, and omics information. The models could then be used to determine protocols for optimal tumor response.
Statements
Author contributions
All authors jointly wrote and revised the manuscript, and approved the final version.
Funding
BG acknowledges partial support by Cancer Prevention and Research Institute of Texas (CPRIT RP150611). BG and SR acknowledge partial support by the Houston Methodist Cancer Center and the Dr. and Mrs. Alan Kaplan Gynecologic Cancer Research Fund. HF acknowledges partial support by the National Institutes of Health/National Cancer Institute grant R15CA203605. SR acknowledges support from the Texas A&M Engineering Experiment Station and the Texas A&M University Department of Biomedical Engineering.
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.
References
1
AberleM. R.BurkhartR. A.TiriacH.Olde DaminkS. W. M.DejongC. H. C.TuvesonD. A.et al (2018). Patient-derived organoid models help define personalized management of gastrointestinal cancer.Br. J. Surg.105e48–e60. 10.1002/bjs.10726
2
AgastinS.GiangU.-B. T.GengY.DeLouiseL. A.KingM. R. (2011). Continuously perfused microbubble array for 3D tumor spheroid model.Biomicrofluidics5:024110. 10.1063/1.3596530
3
AgiostratidouG.SgourosI.GalaniE.VoulgariA.ChondrogianniN.SamantasE.et al (2001). Correlation of in vitro cytotoxicity and clinical response to chemotherapy in ovarian and breast cancer patients.Anticancer Res.21455–459.
4
AhonenI.AkerfeltM.TorisevaM.OswaldE.SchulerJ.NeesM. (2017). A high-content image analysis approach for quantitative measurements of chemosensitivity in patient-derived tumor microtissues.Sci. Rep.7:6600.
5
AlbrittonJ. L.MillerJ. S. (2017). 3D bioprinting: improving in vitro models of metastasis with heterogeneous tumor microenvironments.Dis. Model. Mech.103–14. 10.1242/dmm.025049
6
AndasariV.GerischA.LolasG.SouthA. P.ChaplainM. A. (2011). Mathematical modeling of cancer cell invasion of tissue: biological insight from mathematical analysis and computational simulation.J. Math. Biol.63141–171. 10.1007/s00285-010-0369-1
7
BachmannJ.RaueA.SchillingM.BeckerV.TimmerJ.KlingmullerU. (2012). Predictive mathematical models of cancer signalling pathways.J. Intern. Med.271155–165. 10.1111/j.1365-2796.2011.02492.x
8
BaeY.KataokaK. (2009). Intelligent polymeric micelles from functional poly(ethylene glycol)-poly(amino acid) block copolymers.Adv. Drug Deliv. Rev.61768–784. 10.1016/j.addr.2009.04.016
9
BaeY.NishiyamaN.FukushimaS.KoyamaH.MatsumuraY.KataokaK. (2005). Preparation and biological characterization of polymeric micelle drug carriers with intracellular pH-triggered drug release property: Tumor permeability, controlled subcellular drug distribution, and enhanced in vivo antitumor efficacy.Bioconjugate Chem.16122–130. 10.1021/bc0498166
10
BaeY.NishiyamaN.KataokaK. (2007). In vivo antitumor activity of the folate-conjugated pH-sensitive polymeric micelle selectively releasing adriamycin in the intracellular acidic compartments.Bioconjugate Chem.181131–1139. 10.1021/bc060401p
11
BaptistaP. M.OrlandoG.Mirmalek-SaniS. H.SiddiquiM.AtalaA.SokerS. (2009). Whole organ decellularization - a tool for bioscaffold fabrication and organ bioengineering.Conf. Proc. IEEE Eng. Med. Biol. Soc.20096526–6529.
12
BarbolosiD.CiccoliniJ.LacarelleB.BarlesiF.AndreN. (2016). Computational oncology–mathematical modelling of drug regimens for precision medicine.Nat. Rev. Clin. Oncol.13242–254. 10.1038/nrclinonc.2015.204
13
BarenholzY. (2012). Doxil(R)–the first FDA-approved nano-drug: lessons learned.J. Control Release160117–134. 10.1016/j.jconrel.2012.03.020
14
BertrandN.WuJ.XuX.KamalyN.FarokhzadO. C. (2014). Cancer nanotechnology: the impact of passive and active targeting in the era of modern cancer biology.Adv. Drug Del. Rev.662–25. 10.1016/j.addr.2013.11.009
15
BethgeA.SchumacherU.WedemannG. (2015). Simulation of metastatic progression using a computer model including chemotherapy and radiation therapy.J. Biomed. Inform.5774–87. 10.1016/j.jbi.2015.07.011
16
BhatiaS.FrangioniJ. V.HoffmanR. M.IafrateA. J.PolyakK. (2012). The challenges posed by cancer heterogeneity.Nat. Biotechnol.30604–610. 10.1038/nbt.2294
17
BlomstrandH.ScheiblingU.BratthallC.GreenH.ElanderN. O. (2019). Real world evidence on gemcitabine and nab-paclitaxel combination chemotherapy in advanced pancreatic cancer.BMC Cancer19:40.
18
BojS. F.HwangC.-I.BakerL. A.ChioI. I. C.EngleD. D.CorboV.et al (2015). Organoid models of human and mouse ductal pancreatic cancer.Cell160324–338.
19
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170
20
BorsoiC.LeonardF.LeeY.ZaidM.ElganainyD.AlexanderJ. F.et al (2017). Gemcitabine enhances the transport of nanovector-albumin-bound paclitaxel in gemcitabine-resistant pancreatic ductal adenocarcinoma.Cancer Lett.403296–304. 10.1016/j.canlet.2017.06.026
21
BorsoiC.YokoiK.LeonardF.FerrariM.GodinB. (2015). “Enhanced therapeutic efficacy of a combination of gemcitabine and albumin-bound paclitaxel in multistage nanovectors in pancreatic ductal adenocarcinoma: Evaluation of transport phenomena,”Proceedings AACR 106th Annual Meeting 2015 (Philadelphia, PA: Cancer Res 2015), 2833.
22
BourzacK. (2012). Nanotechnology: carrying drugs.Nature491S58–S60.
23
BreslinS.O’DriscollL. (2016). The relevance of using 3D cell cultures, in addition to 2D monolayer cultures, when evaluating breast cancer drug sensitivity and resistance.Oncotarget745745–45756. 10.18632/oncotarget.9935
24
BrocatoT. A.CokerE. N.DurfeeP. N.LinY. S.TownsonJ.WyckoffE. F.et al (2018). Understanding the connection between nanoparticle uptake and cancer treatment efficacy using mathematical modeling.Sci. Rep.8:7538.
25
BulinA. L.BroekgaardenM.HasanT. (2017). Comprehensive high-throughput image analysis for therapeutic efficacy of architecturally complex heterotypic organoids.Sci. Rep.7:16645.
26
ByrneH. M.OwenM. R.AlarconT.MurphyJ.MainiP. K. (2006). Modelling the response of vascular tumours to chemotherapy: A multiscale approach.Math. Mod. Meth. Appl. S161219–1241. 10.1142/s0218202506001522
27
CampbellA.SivakumaranT.DavidsonM.LockM.WongE. (2008). Mathematical modeling of liver metastases tumour growth and control with radiotherapy.Phys. Med. Biol.537225–7239. 10.1088/0031-9155/53/24/015
28
CasterJ. M.CallaghanC.SeyedinS. N.HendersonK.SunB.WangA. Z. (2019). Optimizing Advances in Nanoparticle Delivery for Cancer Immunotherapy.Adv. Drug Deliv. Rev.1443–15. 10.1016/j.addr.2019.07.009
29
ChamseddineI. M.FrieboesH. B.KokkolarasM. (2018). Design optimization of tumor vasculature-bound nanoparticles.Sci. Rep.8:17768.
30
ChamseddineI. M.FrieboesH. B.KokkolarasM. (2020). Multi-objective optimization of tumor response to drug release from vasculature-bound nanoparticles.Sci. Rep.10:8294.
31
ChaplainM. A. J. (2011). Multiscale mathematical modelling in biology and medicine.IMA J. Appl. Math.76371–388. 10.1093/imamat/hxr025
32
ChenJ.WeihsD.VermolenF. J. (2020). Computational modeling of therapy on pancreatic cancer in its early stages.Biomech. Model Mechanobiol.19427–444. 10.1007/s10237-019-01219-0
33
ChoE. J.HolbackH.LiuK. C.AbouelmagdS. A.ParkJ.YeoY. (2013). Nanoparticle characterization: state of the art. Challenges, and emerging technologies.Mol. Pharm.102093–2110. 10.1021/mp300697h
34
ChoK.EvansB. S.WoodB. M.KumarR.ErbT. J.WarlickB. P.et al (2014). Integration of untargeted metabolomics with transcriptomics reveals active metabolic pathways.Metabolomics11503–517. 10.1007/s11306-014-0713-3
35
ChoiM. R.Stanton-MaxeyK. J.StanleyJ. K.LevinC. S.BardhanR.AkinD.et al (2007). A cellular Trojan Horse for delivery of therapeutic nanoparticles into tumors.Nano Lett.73759–3765. 10.1021/nl072209h
36
ChungA. S.LeeJ.FerraraN. (2010). Targeting the tumour vasculature: insights from physiological angiogenesis.Nat. Rev. Cancer10505–514. 10.1038/nrc2868
37
ClendenonS. G.FuX.Von HoeneR. A.ClendenonJ. L.SlukaJ. P.WinfreeS.et al (2019). Spatial temporal analysis of fieldwise flow in microvasculature.J. Vis. Exp.153.
38
CristR. M.GrossmanJ. H.PatriA. K.SternS. T.DobrovolskaiaM. A.AdiseshaiahP. P.et al (2013). Common pitfalls in nanotechnology: lessons learned from NCI’s Nanotechnology Characterization Laboratory.Integr. Biol.566–73. 10.1039/c2ib20117h
39
CristiniV.FrieboesH. B.LiX.LowengrubJ.MacklinP.SangaS.et al (2008). Nonlinear Modeling and Simulation of Tumor Growth. Selected Topics in Cancer Modeling.Boston: Birkhäuser.
40
CuiY. L.LiH. K.ZhouH. Y.ZhangT.LiQ. (2013). Correlations of tumor-associated macrophage subtypes with liver metastases of colorectal cancer.Asian Pac. J. Cancer Prev.141003–1007. 10.7314/apjcp.2013.14.2.1003
41
CukiermanE.KhanD. R. (2010). The benefits and challenges associated with the use of drug delivery systems in cancer therapy.Biochem. Pharmacol.80762–770. 10.1016/j.bcp.2010.04.020
42
CurtisL. T.EnglandC. G.WuM.LowengrubJ.FrieboesH. B. (2016a). An interdisciplinary computational/experimental approach to evaluate drug-loaded gold nanoparticle tumor cytotoxicity.Nanomedicine11197–216. 10.2217/nnm.15.195
43
CurtisL. T.FrieboesH. B. (2016). The tumor microenvironment as a barrier to cancer nanotherapy.Adv. Exp. Med. Biol.936165–190. 10.1007/978-3-319-42023-3_9
44
CurtisL. T.RychahouP.BaeY.FrieboesH. B. (2016b). A computational/experimental assessment of antitumor activity of polymer nanoassemblies for ph-controlled drug delivery to primary and metastatic tumors.Pharm. Res.332552–2564.
45
CurtisL. T.van BerkelV. H.FrieboesH. B. (2018). Pharmacokinetic/pharmacodynamic modeling of combination-chemotherapy for lung cancer.J. Theor. Biol.44838–52. 10.1016/j.jtbi.2018.03.035
46
CurtisL. T.WuM.LowengrubJ.DecuzziP.FrieboesH. B. (2015). Computational modeling of tumor response to drug release from vasculature-bound nanoparticles.PLoS One10:e0144888. 10.1371/journal.pone.0144888
47
DaiX.MaC.LanQ.XuT. (2016). 3D bioprinted glioma stem cells for brain tumor model and applications of drug susceptibility.Biofabrication8:045005. 10.1088/1758-5090/8/4/045005
48
DalyJ. M.ButlerJ.KemenyN.YehS. D.RidgeJ. A.BotetJ.et al (1985). Predicting tumor response in patients with colorectal hepatic metastases.Ann. Surg.202384–393. 10.1097/00000658-198509000-00017
49
DanetI. M.SemelkaR. C.LeonardouP.BragaL.VaideanG.WoosleyJ. T.et al (2003a). Spectrum of MRI appearances of untreated metastases of the liver.AJR Am. J. Roentgenol.181809–817. 10.2214/ajr.181.3.1810809
50
DanetI. M.SemelkaR. C.NagaseL. L.WooselyJ. T.LeonardouP.ArmaoD. (2003b). Liver metastases from pancreatic adenocarcinoma: MR imaging characteristics.J. Magn. Reson. Imaging18181–188. 10.1002/jmri.10337
51
D’AngeloE.NatarajanD.SensiF.AjayiO.FassanM.MammanoE.et al (2020). Patient-Derived Scaffolds of Colorectal Cancer Metastases as an Organotypic 3D Model of the Liver Metastatic Microenvironment.Cancers12:364. 10.3390/cancers12020364
52
DanhierF. (2016). To exploit the tumor microenvironment: Since the EPR effect fails in the clinic, what is the future of nanomedicine?J. Control Release244108–121. 10.1016/j.jconrel.2016.11.015
53
DawidczykC. M.KimC.ParkJ. H.RussellL. M.LeeK. H.PomperM. G.et al (2014). State-of-the-art in design rules for drug delivery platforms: lessons learned from FDA-approved nanomedicines.J. Control Release187133–144. 10.1016/j.jconrel.2014.05.036
54
DecuzziP.FerrariM. (2006). The adhesive strength of non-spherical particles mediated by specific interactions.Biomaterials275307–5314. 10.1016/j.biomaterials.2006.05.024
55
DecuzziP.PasqualiniR.ArapW.FerrariM. (2009). Intravascular delivery of particulate systems: does geometry really matter?Pharm. Res.26235–243. 10.1007/s11095-008-9697-x
56
DeisboeckT. S.WangZ.MacklinP.CristiniV. (2011). Multiscale cancer modeling.Annu. Rev. Biomed. Eng.13127–155.
57
DewhirstM. W.SecombT. W. (2017). Transport of drugs from blood vessels to tumour tissue.Nat. Rev. Cancer17738–750. 10.1038/nrc.2017.93
58
DograP.ButnerJ. D.ChuangY. L.CasertaS.GoelS.BrinkerC. J.et al (2019). Mathematical modeling in cancer nanomedicine: a review.Biomed. Microdevices21:40.
59
DograP.ButnerJ. D.Ruiz RamirezJ.ChuangY. L.NoureddineA.Jeffrey BrinkerC.et al (2020a). A mathematical model to predict nanomedicine pharmacokinetics and tumor delivery.Comput. Struct. Biotechnol. J.18518–531. 10.1016/j.csbj.2020.02.014
60
DograP.RamírezJ. P.PeláezM. J.WangZ.CristiniV.ParasherG.et al (2020b). Mathematical modeling to address challenges in pancreatic cancer.Curr. Top. Med. Chem.20367–376. 10.2174/1568026620666200101095641
61
DongM.PhilippiC.LoretzB.NafeeN.SchaeferU. F.FriedelG.et al (2011). Tissue slice model of human lung cancer to investigate telomerase inhibition by nanoparticle delivery of antisense 2’-O-methyl-RNA.Int. J. Pharm.41933–42. 10.1016/j.ijpharm.2011.07.009
62
DongY.HuJ.FanL.ChenQ. R. N. A. - (2017). Seq-based transcriptomic and metabolomic analysis reveal stress responses and programmed cell death induced by acetic acid in Saccharomyces cerevisiae.Sci. Rep.7:42659.
63
DrasdoD.HoehmeS.HengstlerJ. G. (2014). How predictive quantitative modelling of tissue organisation can inform liver disease pathogenesis.J. Hepatol.61951–956. 10.1016/j.jhep.2014.06.013
64
Dutta-MoscatoJ.SolovyevA.MiQ.NishikawaT.Soto-GutierrezA.FoxI. J.et al (2014). A Multiscale Agent-Based in silico Model of Liver Fibrosis Progression.Front. Bioeng. Biotechnol.2:18.
65
EdelmanL. B.EddyJ. A.PriceN. D. (2010). In silico models of cancer.Wiley Interdiscip. Rev. Syst. Biol. Med.2438–459.
66
EnderlingH.ChaplainM. A.AndersonA. R.VaidyaJ. S. (2007). A mathematical model of breast cancer development, local treatment and recurrence.J. Theor. Biol.246245–259. 10.1016/j.jtbi.2006.12.010
67
Enriquez-NavasP. M.KamY.DasT.HassanS.SilvaA.ForoutanP.et al (2016). Exploiting evolutionary principles to prolong tumor control in preclinical models of breast cancer.Sci. Transl. Med.8:327ra24. 10.1126/scitranslmed.aad7842
68
EweA.HobelS.HeineC.MerzL.KallendruschS.BechmannI.et al (2017). Optimized polyethylenimine (PEI)-based nanoparticles for siRNA delivery, analyzed in vitro and in an ex vivo tumor tissue slice culture model.Drug Deliv. Transl. Res.7206–216. 10.1007/s13346-016-0306-y
69
FangJ.NakamuraH.MaedaH. (2011). The EPR effect: Unique features of tumor blood vessels for drug delivery, factors involved, and limitations and augmentation of the effect.Adv. Drug Deliv. Rev.63136–151. 10.1016/j.addr.2010.04.009
70
Farokhzad OmidC.LangerR. (2009). Impact of nanotechnology on drug delivery.ACS Nano316–20.
71
FernandesE.FerreiraJ. A.AndreiaP.LuísL.BarrosoS.SarmentoB.et al (2015). New trends in guided nanotherapies for digestive cancers: A systematic review.J. Control. Release209288–307. 10.1016/j.jconrel.2015.05.003
72
Figueroa-MagalhaesM. C.JelovacD.ConnollyR.WolffA. C. (2014). Treatment of HER2-positive breast cancer.Breast23128–136.
73
FolkmanJ. (2007). Angiogenesis: an organizing principle for drug discovery?Nat. Rev. Drug Discov.6273–286. 10.1038/nrd2115
74
FrieboesH. B.ChaplainM. A.ThompsonA. M.BearerE. L.LowengrubJ. S.CristiniV. (2011). Physical oncology: a bench-to-bedside quantitative and predictive approach.Cancer Res.arch71298–302. 10.1158/0008-5472.can-10-2676
75
FrieboesH. B.EdgertonM. E.FruehaufJ. P.RoseF. R.WorrallL. K.GatenbyR. A.et al (2009). Prediction of drug response in breast cancer using integrative experimental/computational modeling.Cancer Res.694484–4492. 10.1158/0008-5472.can-08-3740
76
FrieboesH. B.SinekJ. P.NalciogluO.FruehaufJ. P.CristiniV. (2006). Nanotechnology in Cancer Drug Therapy: A Biocomputational Approach. Biological and Biomedical Nanotechnology, Vol. 20p. New York, NY: Springer-Verlag, 435–460.
77
FrieboesH. B.WuM.LowengrubJ.DecuzziP.CristiniV. (2013). A computational model for predicting nanoparticle accumulation in tumor vasculature.PLoS One8:e56876. 10.1371/journal.pone.0056876
78
FrieboesH. B.YokoiK.DaveB.HussainF.GodinB. (2012). Modeling the Tumor Microenvironment as a Biobarrier in Cancer Nanotherapeutics, Vol. 20. Singapore: Pan Stanford Publishing Pte. Ltd, 137–184.
79
FrieboesH. B.DecuzziP.SinekJ. P.FerrariM.CristiniV. (2019). Computational Modeling of Tumor Biobarriers: Implications for Delivery of Nano-Based Therapeutics, Vol. 20. New Jersey: World Scientific, 201–244.
80
FriedmanA.HaoW. (2017). Mathematical modeling of liver fibrosis.Math. Biosci. Eng.14143–164. 10.3934/mbe.2017010
81
FruehaufJ. P. (2002). In vitro assay-assisted treatment selection for women with breast or ovarian cancer.Endocr. Relat. Cancer9171–182. 10.1677/erc.0.0090171
82
FuX.SlukaJ. P.ClendenonS. G.DunnK. W.WangZ.KlaunigJ. E.et al (2018). Modeling of xenobiotic transport and metabolism in virtual hepatic lobule models.PLoS One13:e0198060. 10.1371/journal.pone.0198060
83
GaoY.LiM.ChenB.ShenZ.GuoP.WientjesM. G.et al (2013). Predictive models of diffusive nanoparticle transport in 3-dimensional tumor cell spheroids.AAPS J.15816–831. 10.1208/s12248-013-9478-2
84
GarveyC. M.SpillerE.LindsayD.ChiangC. T.ChoiN. C.AgusD. B.et al (2016). A high-content image-based method for quantitatively studying context-dependent cell population dynamics.Sci. Rep.6:29752.
85
GeQ.ZhangY.HuaW. P.WuY. C.JinX. X.SongS. H.et al (2015). Combination of transcriptomic and metabolomic analyses reveals a JAZ repressor in the jasmonate signaling pathway of Salvia miltiorrhiza.Sci Rep514048.
86
GeeleherP.CoxN. J. (2014). S. HR. Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines.Genome Biol.15:R47.
87
GevertzJ. L. (2011). Computational modeling of tumor response to vascular-targeting therapies–part I: validation.Comput. Math. Methods Med.2011:830515.
88
GioeliD.SnowC. J.SimmersM. B.HoangS. A.FiglerR. A.AllendeJ. A.et al (2019). Development of a multicellular pancreatic tumor microenvironment system using patient-derived tumor cells.Lab Chip191193–1204. 10.1039/c8lc00755a
89
GodinB.ChiappiniC.SrinivasanS.AlexanderJ. F.YokoiK.FerrariM.et al (2012). Discoidal porous silicon particles: fabrication and biodistribution in breast cancer bearing mice.Adv. Funct. Mater.224225–4235. 10.1002/adfm.201200869
90
GodinB.DriessenW. H.PronethB.LeeS. Y.SrinivasanS.RumbautR.et al (2010a). An integrated approach for the rational design of nanovectors for biomedical imaging and therapy.Adv. Genet.6931–64. 10.1016/s0065-2660(10)69009-8
91
GodinB.GuJ.SerdaR. E.BhavaneR.TasciottiE.ChiappiniC.et al (2010b). Tailoring the degradation kinetics of mesoporous silicon structures through PEGylation.J. Biomed. Mater. Res. A941236–1243.
92
GodinB.TasciottiE.LiuX.SerdaR. E.FerrariM. (2011). Multistage nanovectors: from concept to novel imaging contrast agents and therapeutics.Acc Chem. Res.44979–989. 10.1021/ar200077p
93
GreishK. (2010). Enhanced permeability and retention (EPR) effect for anticancer nanomedicine drug targeting.Methods Mol. Biol.62425–37. 10.1007/978-1-60761-609-2_3
94
GriffithL. G.SwartzM. A. (2006). Capturing complex 3D tissue physiology in vitro.Nat. Rev. Mol. Cell Biol.7211–224. 10.1038/nrm1858
95
GrimmettC.PickettK.ShepherdJ.WelchK.Recio-SaucedoA.StreitE.et al (2018). Systematic review of the empirical investigation of resources to support decision-making regarding BRCA1 and BRCA2 genetic testing in women with breast cancer.Patient Educ. Couns.101779–788. 10.1016/j.pec.2017.11.016
96
GrodzinskiP.FarrellD. (2014). Future opportunities in cancer nanotechnology—NCI strategic workshop report.Cancer Res.741307–1310. 10.1158/0008-5472.can-13-2787
97
GullerA. E.RozovaV. S.KuschnerusI.KhabirZ.NadortA.Garcia-BennettA. E.et al (2020). Tissue engineered model of hepatic breast cancer micrometastasis shows host-dependent colonization patterns and drug responses.bioRxiv [Preprint]10.1101/2020.01.08.898163v1
98
GuoS.DengC. X. (2018). Effect of stromal cells in tumor microenvironment on metastasis initiation.Int. J. Biol. Sci.142083–2093. 10.7150/ijbs.25720
99
GustafsonH. H.Holt-CasperD.GraingerD. W.GhandehariH. (2015). Nanoparticle uptake: the phagocyte problem.Nano Today10487–510. 10.1016/j.nantod.2015.06.006
100
HaenoH.GonenM.DavisM. B.HermanJ. M.Iacobuzio-DonahueC. A.MichorF. (2012). Computational modeling of pancreatic cancer reveals kinetics of metastasis suggesting optimum treatment strategies.Cell148362–375. 10.1016/j.cell.2011.11.060
101
HayM.ThomasD. W.CraigheadJ. L.EconomidesC.RosenthalJ. (2014). Clinical development success rates for investigational drugs.Nat. Biotechnol.3240–51. 10.1038/nbt.2786
102
HazhirkarzarB.KhoshpouriP.ShaghaghiM.GhasabehM. A.PawlikT. M.KamelI. R. (2020). Current state of the art imaging approaches for colorectal liver metastasis.Hepatobiliary Surg. Nutr.935–48. 10.21037/hbsn.2019.05.11
103
HeY. F.ZhangM. Y.WuX.SunX. J.XuT.HeQ. Z.et al (2013). High MUC2 expression in ovarian cancer is inversely associated with the M1/M2 ratio of tumor-associated macrophages and patient survival time.PLoS One8:e79769. 10.1371/journal.pone.0079769
104
HeathJ. R.RibasA.MischelP. S. (2016). Single-cell analysis tools for drug discovery and development.Nat. Rev. Drug Discov.15204–216. 10.1038/nrd.2015.16
105
HindererS.LaylandS. L.Schenke-LaylandK. (2016). ECM and ECM-like materials - Biomaterials for applications in regenerative medicine and cancer therapy.Adv. Drug Deliv. Rev.97260–269. 10.1016/j.addr.2015.11.019
106
HinowP.GerleeP.McCawleyL. J.QuarantaV.CiobanuM.WangS.et al (2009). A spatial model of tumor-host interaction: application of chemotherapy.Math. Biosci. Eng.6521–546. 10.3934/mbe.2009.6.521
107
HoehmeS.BertauxF.WeensW.Grasl-KrauppB.HengstlerJ. G.DrasdoD. (2018). Model Prediction and Validation of an Order Mechanism Controlling the Spatiotemporal Phenotype of Early Hepatocellular Carcinoma.Bull. Math. Biol.801134–1171. 10.1007/s11538-017-0375-1
108
HoehmeS.BrulportM.BauerA.BedawyE.SchormannW.HermesM.et al (2010). Prediction and validation of cell alignment along microvessels as order principle to restore tissue architecture in liver regeneration.Proc. Natl. Acad. Sci. U.S.A.10710371–10376. 10.1073/pnas.0909374107
109
HoehmeS.FriebelA.HammadS.DrasdoD.HengstlerJ. G. (2017). Creation of three-dimensional liver tissue models from experimental images for systems medicine.Methods Mol. Biol.1506319–362. 10.1007/978-1-4939-6506-9_22
110
HoehmeS.HengstlerJ. G.BrulportM.SchaferM.BauerA.GebhardtR.et al (2007). Mathematical modelling of liver regeneration after intoxication with CCl(4).Chem. Biol. Interact.16874–93. 10.1016/j.cbi.2007.01.010
111
HolzhutterH. G.DrasdoD.PreusserT.LippertJ.HenneyA. M. (2012). The virtual liver: a multidisciplinary, multilevel challenge for systems biology.Wiley Interdiscip. Rev. Syst. Biol. Med.4221–235. 10.1002/wsbm.1158
112
HormanS. R.HoganC.Delos ReyesK.LoF.AntczakC. (2015). Challenges and opportunities toward enabling phenotypic screening of complex and 3D cell models.Future Med. Chem.7513–525. 10.4155/fmc.14.163
113
HuG.GuoM.XuJ.WuF.FanJ.HuangQ.et al (2019). Nanoparticles Targeting Macrophages as Potential Clinical Therapeutic Agents Against Cancer and Inflammation.Front. Immunol.10:1998.
114
HuangZ. Q.BuchsbaumD. J. (2009). Monoclonal antibodies in the treatment of pancreatic cancer.Immunotherapy1223–229.
115
HudsonS. V.DolinC. E.PooleL. G.MasseyV. L.WilkeyD.BeierJ. I.et al (2017). Modeling the kinetics of integrin receptor binding to hepatic extracellular matrix proteins.Sci. Rep.7:12444.
116
HudsonS. V.MillerH. A.MahlbacherG. E.SaforoD.BeverlyL. J.ArteelG. E.et al (2019). Computational/experimental evaluation of liver metastasis post hepatic injury: interactions with macrophages and transitional ECM.Sci. Rep.9:15077.
117
HusseinK. H.ParkK. M.GhimJ. H.YangS. R.WooH. M. (2016). Three dimensional culture of HepG2 liver cells on a rat decellularized liver matrix for pharmacological studies.J. Biomed. Mater. Res. B Appl. Biomater.104263–273. 10.1002/jbm.b.33384
118
Ibrahim-HashimA.Robertson-TessiM.Enriquez-NavasP. M.DamaghiM.BalagurunathanY.WojtkowiakJ. W.et al (2017). Defining cancer subpopulations by adaptive strategies rather than molecular properties provides novel insights into intratumoral evolution.Cancer Res.772242–2254. 10.1158/0008-5472.can-16-2844
119
JaganathanH.GageJ.LeonardF.SrinivasanS.SouzaG. R.DaveB.et al (2014). Three-dimensional in vitro co-culture model of breast tumor using magnetic levitation.Sci. Rep.4:6468.
120
JainR. K. (2003). Molecular regulation of vessel maturation.Nat. Med.9685–693. 10.1038/nm0603-685
121
KaddiC. D.PhanJ. H.WangM. D. (2013). Computational nanomedicine: modeling of nanoparticle-mediated hyperthermal cancer therapy.Nanomedicine81323–1333. 10.2217/nnm.13.117
122
KangS. A.HasanN.MannA. P.ZhengW.ZhaoL.MorrisL.et al (2015). Blocking the adhesion cascade at the premetastatic niche for prevention of breast cancer metastasis.Mol. Ther.231044–1054. 10.1038/mt.2015.45
123
KhalidA.PersanoS.ShenH.ZhaoY.BlancoE.FerrariM.et al (2016). Strategies for improving drug delivery: nanocarriers and microenvironmental priming.Expert. Opin. Drug Del.141–13.
124
KhawarI. A.KimJ. H.KuhH.-J. (2015). Improving drug delivery to solid tumors: priming the tumor microenvironment.J. Control Release20178–89. 10.1016/j.jconrel.2014.12.018
125
KielbassaK.VegnaS.RamirezC.AkkariL. (2019). Understanding the origin and diversity of macrophages to tailor their targeting in solid cancers.Front. Immunol.10:2215.
126
KimC.BangJ. H.KimY. E.LeeS. H.KangJ. Y. (2012). On-chip anticancer drug test of regular tumor spheroids formed in microwells by a distributive microchannel network.Lab Chip124135–4142.
127
KoontzB. F. (2017). Stereotactic body radiation therapy for oligometastatic prostate cancer: the hunt for the silver bullet.Int. J. Radiat. Oncol. Biol. Phys.99761–763. 10.1016/j.ijrobp.2017.05.020
128
KopperO.de WitteC. J.LohmussaarK.Valle-InclanJ. E.HamiN.KesterL.et al (2019). An organoid platform for ovarian cancer captures intra- and interpatient heterogeneity.Nat. Med.25838–849. 10.1038/s41591-019-0422-6
129
KreegerP. K.LauffenburgerD. A. (2010). Cancer systems biology: a network modeling perspective.Carcinogenesis312–8. 10.1093/carcin/bgp261
130
KyleA. H.HuxhamL. A.YeomanD. M.MinchintonA. I. (2007). Limited tissue penetration of taxanes: A mechanism for resistance in solid tumors.Clin. Cancer Res.132804–2810. 10.1158/1078-0432.ccr-06-1941
131
LeccaP.ReA. (2017). Network-oriented approaches to anticancer drug response.Methods Mol. Biol.1513101–117. 10.1007/978-1-4939-6539-7_8
132
LeeH.HanW.KimH.HaD. H.JangJ.KimB. S.et al (2017). Development of liver decellularized extracellular matrix bioink for three-dimensional cell printing-based liver tissue engineering.Biomacromolecules181229–1237. 10.1021/acs.biomac.6b01908
133
LeeJ. J.HuangJ.EnglandC. G.McNallyL. R.FrieboesH. B. (2013). Predictive modeling of in vivo response to gemcitabine in pancreatic cancer.PLoS Comput. Biol.9:e1003231. 10.1371/journal.pcbi.1003231
134
LeonardF.CurtisL. T.HamedA. R.ZhangC.ChauE.SievingD.et al (2020). Nonlinear response to cancer nanotherapy due to macrophage interactions revealed by mathematical modeling and evaluated in a murine model via CRISPR-modulated macrophage polarization.Cancer Immunol. Immunother.69731–744. 10.1007/s00262-020-02504-z
135
LeonardF.CurtisL. T.WareM. J.NosratT.LiuX.YokoiK.et al (2017). Macrophage polarization contributes to the anti-tumoral efficacy of mesoporous nanovectors loaded with albumin-bound paclitaxel.Front. Immunol.8:693.
136
LeonardF.CurtisL. T.YesantharaoP.TaneiT.AlexanderJ. F.WuM.et al (2016). Enhanced performance of macrophage-encapsulated nanoparticle albumin-bound-paclitaxel in hypo-perfused cancer lesions.Nanoscale812544–12552. 10.1039/c5nr07796f
137
LeonardF.GodinB. (2016). 3D In Vitro Model for Breast Cancer Research Using Magnetic Levitation and Bioprinting Method.Methods Mol. Biol.1406239–251. 10.1007/978-1-4939-3444-7_21
138
LettmannK. A.Hardtke-WolenskiM. (2014). The importance of liver microcirculation in promoting autoimmune hepatitis via maintaining an inflammatory cytokine milieu–a mathematical model study.J. Theor. Biol.34833–46. 10.1016/j.jtbi.2014.01.016
139
LiM.Al-JamalK. T.KostarelosK.ReinekeJ. (2010). Physiologically based pharmacokinetic modeling of nanoparticles.ACS Nano46303–6317.
140
LiM.CzyszczonE.ReinekeJ. (2013). Delineating intracellular pharmacokinetics of paclitaxel delivered by PLGA nanoparticles.Drug Deliv. Transl. Res.3551–561. 10.1007/s13346-013-0162-y
141
LiM.PanagiZ.AvgoustakisK.ReinekeJ. (2012). Physiologically based pharmacokinetic modeling of PLGA nanoparticles with varied mPEG content.Int. J. Nanomed.71345–1356.
142
LiM.ReinekeJ. (2011). Mathematical modelling of nanoparticle biodistribution: extrapolation among intravenous, oral and pulmonary administration routes.Int. J. Nano Biomater.3222–223.
143
LiS.TodorA.LuoR. (2016). Blood transcriptomics and metabolomics for personalized medicine.Comput. Struct. Biotechnol. J.141–7. 10.1016/j.csbj.2015.10.005
144
LiX.KwonH. (2019). Efficacy and safety of nanoparticle albumin-bound paclitaxel in elderly patients with metastatic breast cancer: a meta-analysis.J. Clin. Med.8:1689. 10.3390/jcm8101689
145
LiuL. X.ZhangW. H.JiangH. C. (2003). Current treatment for liver metastases from colorectal cancer.World J. Gastroentero9193–200.
146
LouzounY.XueC.LesinskiG. B.FriedmanA. (2014). A mathematical model for pancreatic cancer growth and treatments.J. Theor. Biol.35174–82. 10.1016/j.jtbi.2014.02.028
147
LowengrubJ. S.FrieboesH. B.JinF.ChuangY. L.LiX.MacklinP.et al (2010). Nonlinear modelling of cancer: bridging the gap between cells and tumours.Nonlinearity23R1–R9.
148
LuS.CuzzucoliF.JiangJ.LiangL. G.WangY.KongM.et al (2018). Development of a biomimetic liver tumor-on-a-chip model based on decellularized liver matrix for toxicity testing.Lab Chip183379–3392. 10.1039/c8lc00852c
149
MaedaH.FangJ.InutsukaT.KitamotoY. (2003). Vascular permeability enhancement in solid tumor: various factors, mechanisms involved and its implications.Int. Immunopharmacol.3319–328. 10.1016/s1567-5769(02)00271-0
150
MahlbacherG.CurtisL. T.LowengrubJ.FrieboesH. B. (2018). Mathematical modeling of tumor-associated macrophage interactions with the cancer microenvironment.J. Immunother. Cancer6:10.
151
MahlbacherG. E.ReihmerK. C.FrieboesH. B. (2019). Mathematical modeling of tumor-immune cell interactions.J. Theor. Biol.46947–60. 10.1016/j.jtbi.2019.03.002
152
ManiwaY.YoshimuraM.HashimotoS.TakataM.NishioW. (2010). Chemosensitivity of lung cancer: differences between the primary lesion and lymph node metastasis.Oncol. Lett.1345–349. 10.3892/ol_00000061
153
MehtaG.HsiaoA. Y.IngramM.LukerG. D.TakayamaS. (2012). Opportunities and challenges for use of tumor spheroids as models to test drug delivery and efficacy.J. Control. Release164192–204. 10.1016/j.jconrel.2012.04.045
154
MengF.MeyerC. M.JoungD.ValleraD. A.McAlpineM. C.Panoskaltsis-MortariA. (2019). 3D bioprinted in vitro metastatic models via reconstruction of tumor microenvironments.Adv. Mater.31:e1806899.
155
MeyerK.OstrenkoO.BourantasG.Morales-NavarreteH.Porat-ShliomN.Segovia-MirandaF.et al (2017). A predictive 3D multi-scale model of biliary fluid dynamics in the liver lobule.Cell Syst.4:e9.
156
MichorF.LiphardtJ.FerrariM.WidomJ. (2011). What does physics have to do with cancer?Nat. Rev. Cancer11657–670. 10.1038/nrc3092
157
MillerH. A.FrieboesH. B. (2019a). Evaluation of Drug-Loaded Gold Nanoparticle Cytotoxicity as a Function of Tumor Vasculature-Induced Tissue Heterogeneity.Ann. Biomed. Eng.47257–271. 10.1007/s10439-018-02146-4
158
MillerH. A.FrieboesH. B. (2019b). Pharmacokinetic/Pharmacodynamics Modeling of Drug-Loaded PLGA Nanoparticles Targeting Heterogeneously Vascularized Tumor Tissue.Pharm. Res.36:185.
159
MintonD. R.FuL.ChenQ.RobinsonB. D.GrossS. S.NanusD. M.et al (2015). Analyses of the transcriptome and metabolome demonstrate that HIF1alpha mediates altered tumor metabolism in clear cell renal cell carcinoma.PLoS One10:e0120649. 10.1371/journal.pone.0120649
160
MiryalaS. K.AnbarasuA.RamaiahS. (2018). Discerning molecular interactions: A comprehensive review on biomolecular interaction databases and network analysis tools.Gene64284–94. 10.1016/j.gene.2017.11.028
161
MisraR.AcharyaS.SahooS. K. (2010). Cancer nanotechnology: application of nanotechnology in cancer therapy.Drug Discov. Today15842–850. 10.1016/j.drudis.2010.08.006
162
MucakiE. J.ZhaoJ. Z. L.LizotteD. J.RoganP. K. (2019). Predicting responses to platin chemotherapy agents with biochemically-inspired machine learning.Signal Transduct. Target Ther.4:1.
163
MurakamiY.KuboS.TamoriA.ItamiS.KawamuraE.IwaisakoK.et al (2015). Comprehensive analysis of transcriptome and metabolome analysis in Intrahepatic Cholangiocarcinoma and Hepatocellular Carcinoma.Sci. Rep.5:16294.
164
NakamuraH.FangJ.MaedaH. (2015). Development of next-generation macromolecular drugs based on the EPR effect: challenges and pitfalls.Expert. Opin. Drug Deliv.1253–64. 10.1517/17425247.2014.955011
165
NamasivayamS.MartinD. R.SainiS. (2007). Imaging of liver metastases: MRI.Cancer Imaging72–9. 10.1102/1470-7330.2007.0002
166
NgC. F.FrieboesH. B. (2017). Model of vascular desmoplastic multispecies tumor growth.J. Theor. Biol.430245–282. 10.1016/j.jtbi.2017.05.013
167
NgC. F.FrieboesH. B. (2018). Simulation of multispecies desmoplastic cancer growth via a fully adaptive non-linear full multigrid algorithm.Front. Physiol.9:821.
168
NgamcherdtrakulW.YantaseeW. (2019). siRNA therapeutics for breast cancer: recent efforts in targeting metastasis, drug resistance, and immune evasion.Transl. Res.214105–120. 10.1016/j.trsl.2019.08.005
169
NishiiK.ReeseG.MoranE. C.SparksJ. L. (2016). Multiscale computational model of fluid flow and matrix deformation in decellularized liver.J. Mech. Behav. Biomed. Mater.57201–214. 10.1016/j.jmbbm.2015.11.033
170
OsborneJ. M.WalterA.KershawS. K.MiramsG. R.FletcherA. G.PathmanathanP.et al (2010). A hybrid approach to multi-scale modelling of cancer.Philos. Trans. A Math. Phys. Eng. Sci.3685013–5028.
171
OzcelikkaleA.MoonH. R.LinnesM.HanB. (2017). In vitro microfluidic models of tumor microenvironment to screen transport of drugs and nanoparticles.Wiley Interdiscip. Rev. Nanomed. Nanobiotechnol.9:e1460. 10.1002/wnan.1460
172
ParhiP.MohantyC.SahooS. K. (2012). Nanotechnology-based combinational drug delivery: an emerging approach for cancer therapy.Drug Discov. Today171044–1052. 10.1016/j.drudis.2012.05.010
173
PascalJ.AshleyC. E.WangZ.BrocatoT. A.ButnerJ. D.CarnesE. C.et al (2013a). Mechanistic modeling identifies drug-uptake history as predictor of tumor drug resistance and nano-carrier-mediated response.ACS Nano711174–11182. 10.1021/nn4048974
174
PascalJ.BearerE. L.WangZ.KoayE. J.CurleyS. A.CristiniV. (2013b). Mechanistic patient-specific predictive correlation of tumor drug response with microenvironment and perfusion measurements.Proc. Natl. Acad. Sci. U.S.A.11014266–14271. 10.1073/pnas.1300619110
175
PengF.LiuY.HeC.KongY.OuyangQ.XieX.et al (2018). Prediction of platinum-based chemotherapy efficacy in lung cancer based on LC-MS metabolomics approach.J. Pharm. Biomed. Anal.15495–101. 10.1016/j.jpba.2018.02.051
176
PezzellaF.PastorinoU.TagliabueE.AndreolaS.SozziG.GaspariniG.et al (1997). Non-small-cell lung carcinoma tumor growth without morphological evidence of neo-angiogenesis.Am. J. Pathol.1511417–1423.
177
PohA. R.ErnstM. (2018). Targeting Macrophages in Cancer: From Bench to Bedside.Front. Oncol.8:49.
178
PollardJ. W. (2008). Macrophages define the invasive microenvironment in breast cancer.J. Leukoc Biol.84623–630. 10.1189/jlb.1107762
179
PontaA.BaeY. (2014). Tumor-preferential sustained drug release enhances antitumor activity of block copolymer micelles.J. Drug Target22619–628. 10.3109/1061186x.2014.910793
180
PontaA.FugitK. D.AndersonB. D.BaeY. (2015). Release, partitioning, and conjugation stability of doxorubicin in polymer micelles determined by mechanistic modeling.Pharm. Res.321752–1763. 10.1007/s11095-014-1573-2
181
PonzoniM.PastorinoF.Di PaoloD.PerriP.BrignoleC. (2018). Targeting macrophages as a potential therapeutic intervention: impact on inflammatory diseases and cancer.Int. J. Mol. Sci.19:1953. 10.3390/ijms19071953
182
PrabhakarU.MaedaH.JainR. K.Sevick-MuracaE. M.ZamboniW.FarokhzadO. C.et al (2013). Challenges and key considerations of the enhanced permeability and retention effect for nanomedicine drug delivery in oncology.Cancer Res.732412–2417. 10.1158/0008-5472.can-12-4561
183
PyonteckS. M.AkkariL.SchuhmacherA. J.BowmanR. L.SevenichL.QuailD. F.et al (2013). CSF-1R inhibition alters macrophage polarization and blocks glioma progression.Nat. Med.191264–1272. 10.1038/nm.3337
184
RaghavanS.MehtaP.HorstE. N.WardM. R.RowleyK. R.MehtaG. (2016). Comparative analysis of tumor spheroid generation techniques for differential in vitro drug toxicity.Oncotarget716948–16961. 10.18632/oncotarget.7659
185
RaghavanS.MehtaP.WardM. R.BregenzerM. E.FleckE. M. A.TanL.et al (2017). Personalized medicine-based approach to model patterns of chemoresistance and tumor recurrence using ovarian cancer stem cell spheroids.Clin Cancer Res236934–6945. 10.1158/1078-0432.ccr-17-0133
186
RaghavanS.MehtaP.XieY.LeiY. L.MehtaG. (2019). Ovarian cancer stem cells and macrophages reciprocally interact through the WNT pathway to promote pro-tumoral and malignant phenotypes in 3D engineered microenvironments.J. Immunother. Cancer7:190.
187
RaghavanS.WardM. R.RowleyK. R.WoldR. M.TakayamaS.BuckanovichR. J.et al (2015). Formation of stable small cell number three-dimensional ovarian cancer spheroids using hanging drop arrays for preclinical drug sensitivity assays.Gynecol. Oncol.138181–189. 10.1016/j.ygyno.2015.04.014
188
RaniH. P.SheuT. W.ChangT. M.LiangP. C. (2006). Numerical investigation of non-Newtonian microcirculatory blood flow in hepatic lobule.J. Biomech.39551–563. 10.1016/j.jbiomech.2004.11.029
189
ReichelD.CurtisL. T.EhlmanE.EversB. M.RychahouP.FrieboesH. B.et al (2017). Development of halofluorochromic polymer nanoassemblies for the potential detection of liver metastatic colorectal cancer tumors using experimental and computational approaches.Pharm. Res.342385–2402. 10.1007/s11095-017-2245-9
190
ReichelD.RychahouP.BaeY. (2015). Polymer nanoassemblies with solvato- and halo-fluorochromic properties for real-time monitoring of drug release and pH-enhanced ex vivo imaging of metastatic tumors.Ther. Deliv.61221–1237. 10.4155/tde.15.59
191
ReichelD.TripathiM.PerezJ. M. (2019). Biological effects of nanoparticles on macrophage polarization in the tumor microenvironment.Nanotheranostics366–88. 10.7150/ntno.30052
192
RejniakK. A.AndersonA. R. (2011). Hybrid models of tumor growth.Wiley Interdiscip. Rev. Syst. Biol. Med.3115–125. 10.1002/wsbm.102
193
RejniakK. A.McCawleyL. J. (2010). Current trends in mathematical modeling of tumor-microenvironment interactions: a survey of tools and applications.Exp. Biol. Med.235411–423.
194
RenS.ShaoY.ZhaoX.HongC. S.WangF.LuX.et al (2016). Integration of metabolomics and transcriptomics reveals major metabolic pathways and potential biomarker involved in prostate cancer.Mol. Cell Proteomics15154–163. 10.1074/mcp.m115.052381
195
RibbaB.SautO.ColinT.BreschD.GrenierE.BoisselJ. P. (2006). A multiscale mathematical model of avascular tumor growth to investigate the therapeutic benefit of anti-invasive agents.J. Theor. Biol.243532–541. 10.1016/j.jtbi.2006.07.013
196
RickenT.WernerD.HolzhutterH. G.KonigM.DahmenU.DirschO. (2015). Modeling function-perfusion behavior in liver lobules including tissue, blood, glucose, lactate and glycogen by use of a coupled two-scale PDE-ODE approach.Biomech. Model. Mechanobiol.14515–536. 10.1007/s10237-014-0619-z
197
RoederI.GlaucheI. (2008). Pathogenesis, treatment effects, and resistance dynamics in chronic myeloid leukemia–insights from mathematical model analyses.J. Mol. Med.8617–27. 10.1007/s00109-007-0241-y
198
RoyM.FinleyS. D. (2017). Computational model predicts the effects of targeting cellular metabolism in pancreatic cancer.Front. Physiol.8:217.
199
RuffellB.AffaraN. I.CoussensL. M. (2012). Differential macrophage programming in the tumor microenvironment.Trends Immunol.33119–126. 10.1016/j.it.2011.12.001
200
RuffellB.CoussensL. M. (2015). Macrophages and therapeutic resistance in cancer.Cancer Cell27462–472. 10.1016/j.ccell.2015.02.015
201
RummenyE. J.MarchalG. (1997). Liver imaging. Clinical applications and future perspectives.Acta Radiol.38626–630. 10.1080/02841859709172392
202
SachsN.CleversH. (2014). Organoid cultures for the analysis of cancer phenotypes.Curr. Opin. Genet. Dev.2468–73. 10.1016/j.gde.2013.11.012
203
SangaS.FrieboesH. B.SinekJ. P.CristiniV. A. (2007). “Multiscale Approach for Computational Modeling of Biobarriers to Cancer Chemotherapy via Nanotechnology,” in Cancer Nanotechnology: Nanomaterials for Cancer Diagnosis And Therapy, edsNalwaH. S.WebsterT. J. (Valencia, California: American Scientific), 1–21.
204
SantoV. E.RebeloS. P.EstradaM. F.AlvesP. M.BoghaertE.BritoC. (2017). Drug screening in 3D in vitro tumor models: overcoming current pitfalls of efficacy read-outs.Biotechnol. J.12:1600505. 10.1002/biot.201600505
205
SarodeP.SchaeferM. B.GrimmingerF.SeegerW.SavaiR. (2020). Macrophage and tumor cell cross-talk is fundamental for lung tumor progression: we need to talk.Front. Oncol.10:324.
206
SatoT.StangeD. E.FerranteM.VriesR. G.Van EsJ. H.Van den BrinkS.et al (2011). Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett’s epithelium.Gastroenterology1411762–1772. 10.1053/j.gastro.2011.07.050
207
SchliessF.HoehmeS.HenkelS. G.GhallabA.DrieschD.BottgerJ.et al (2014). Integrated metabolic spatial-temporal model for the prediction of ammonia detoxification during liver damage and regeneration.Hepatology602040–2051. 10.1002/hep.27136
208
SchwenL. O.SchenkA.KreutzC.TimmerJ.Bartolome RodriguezM. M.KuepferL.et al (2015). Representative sinusoids for hepatic four-scale pharmacokinetics simulations.PLoS One10:e0133653. 10.1371/journal.pone.0133653
209
SiggersJ. H.LeungchavaphongseK.HoC. H.RepettoR. (2014). Mathematical model of blood and interstitial flow and lymph production in the liver.Biomech. Model Mechanobiol.13363–378. 10.1007/s10237-013-0516-x
210
SinekJ. P.FrieboesH. B.SivaramanB.SangaS.CristiniV. (2007). Mathematical and Computational Modeling: Towards the Development and Application of Nanodevices for Drug Delivery. Nanotechnologies for the Life Sciences.Hoboken, NJ: Wiley-VCH Verlag GmbH & Co. KGaA.
211
SinekJ. P.SangaS.ZhengX.FrieboesH. B.FerrariM.CristiniV. (2009). Predicting drug pharmacokinetics and effect in vascularized tumors using computer simulation.J. Math. Biol.58485–510. 10.1007/s00285-008-0214-y
212
SlukaJ. P.FuX.SwatM.BelmonteJ. M.CosmanescuA.ClendenonS. G.et al (2016). A Liver-Centric Multiscale Modeling Framework for Xenobiotics.PLoS One11:e0162428. 10.1371/journal.pone.0162428
213
SrinivasanS.AlexanderJ. F.DriessenW. H.LeonardF.YeH.LiuX.et al (2013). Bacteriophage associated silicon particles: design and characterization of a novel theranostic vector with improved payload carrying potential.J. Mater. Chem. B Mater. Biol. Med.1:10.1039/C3TB20595A.
214
StamatakosG. S.AntipasV. P.UzunogluN. K. (2006). A spatiotemporal, patient individualized simulation model of solid tumor response to chemotherapy in vivo: the paradigm of glioblastoma multiforme treated by temozolomide.IEEE Trans. Biomed. Eng.531467–1477. 10.1109/tbme.2006.873761
215
StesselsF.Van Den EyndenG.Van Der AuweraI.SalgadoR.Van Den HeuvelE.HarrisA. L.et al (2004). Breast adenocarcinoma liver metastases, in contrast to colorectal cancer liver metastases, display a non-angiogenic growth pattern that preserves the stroma and lacks hypoxia.Br. J. Cancer901429–1436. 10.1038/sj.bjc.6601727
216
StockK.EstradaM. F.VidicS.GjerdeK.RudischA.SantoV. E.et al (2016). Capturing tumor complexity in vitro: Comparative analysis of 2D and 3D tumor models for drug discovery.Sci. Rep.6:28951.
217
SuG.BurantC. F.BeecherC. W.AtheyB. D.MengF. (2011). Integrated metabolome and transcriptome analysis of the NCI60 dataset.BMC Bioinformatics12(Suppl. 1):S36.
218
SunD.LiuY.WangH.DengF.ZhangY.ZhaoS.et al (2018). Novel decellularized liver matrix-alginate hybrid gel beads for the 3D culture of hepatocellular carcinoma cells.Int. J. Biol. Macromol.1091154–1163. 10.1016/j.ijbiomac.2017.11.103
219
SwansonK. R.RockneR. C.ClaridgeJ.ChaplainM. A.AlvordE. C.Jr.AndersonA. R. (2011). Quantifying the role of angiogenesis in malignant progression of gliomas: in silico modeling integrates imaging and histology.Cancer Res.717366–7375. 10.1158/0008-5472.can-11-1399
220
TanakaT.GodinB.BhavaneR.Nieves-AliceaR.GuJ.LiuX.et al (2010). In vivo evaluation of safety of nanoporous silicon carriers following single and multiple dose intravenous administrations in mice.Int. J. Pharm.402190–197. 10.1016/j.ijpharm.2010.09.015
221
TaneiT.LeonardF.LiuX.AlexanderJ. F.SaitoY.FerrariM.et al (2016). Redirecting transport of nanoparticle albumin-bound paclitaxel to macrophages enhances therapeutic efficacy against liver metastases.Cancer Res.76429–439. 10.1158/0008-5472.can-15-1576
222
TariqM.ZhangJ.LiangG.DingL.HeQ.YangB. (2017). Macrophage polarization: anti-cancer strategies to target tumor-associated macrophage in breast cancer.J. Cell Biochem.1182484–2501. 10.1002/jcb.25895
223
TasciottiE.GodinB.MartinezJ. O.ChiappiniC.BhavaneR.LiuX.et al (2011). Near-infrared imaging method for the in vivo assessment of the biodistribution of nanoporous silicon particles.Mol. Imaging1056–68.
224
TianY.WangZ.LiuX.DuanJ.FengG.YinY.et al (2018). Prediction of chemotherapeutic efficacy in non-small cell lung cancer by serum metabolomic profiling.Clin. Cancer Res.242100–2109. 10.1158/1078-0432.ccr-17-2855
225
UngerC.KramerN.WalzlA.ScherzerM.HengstschlagerM.DolznigH. (2014). Modeling human carcinomas: physiologically relevant 3D models to improve anti-cancer drug development.Adv. Drug Deliv. Rev.79-8050–67. 10.1016/j.addr.2014.10.015
226
UntchM.JackischC.SchneeweissA.SchmatlochS.AktasB.DenkertC.et al (2019). NAB-Paclitaxel Improves Disease-Free Survival in Early Breast Cancer: GBG 69-GeparSepto.J. Clin. Oncol.372226–2234. 10.1200/jco.18.01842
227
van de VenA. L.AbdollahiB.MartinezC. J.BureyL. A.LandisM. D.ChangJ. C.et al (2013). Modeling of nanotherapeutics delivery based on tumor perfusion.New J. Phys.15:55004.
228
van de VenA. L.WuM.LowengrubJ.McDougallS. R.ChaplainM. A.CristiniV.et al (2012). Integrated intravital microscopy and mathematical modeling to optimize nanotherapeutics delivery to tumors.AIP Adv.2:11208.
229
van den EyndenG. G.MajeedA. W.IllemannM.VermeulenP. B.BirdN. C.Høyer-HansenG.et al (2013). The multifaceted role of the microenvironment in liver metastasis: Biology and clinical implications.Cancer Res.732031–2043. 10.1158/0008-5472.can-12-3931
230
Van LiedekerkeP.NeitschJ.JohannT.WarmtE.Gonzalez-ValverdeI.HoehmeS.et al (2020). A quantitative high-resolution computational mechanics cell model for growing and regenerating tissues.Biomech. Model Mechanobiol.19189–220. 10.1007/s10237-019-01204-7
231
VendittoV. J.SzokaF. C. (2013). Cancer nanomedicines: so many papers and so few drugs!Adv. Drug Del. Rev.6580–88. 10.1016/j.addr.2012.09.038
232
VicentM. J.RingsdorfH.DuncanR. (2009). Polymer therapeutics: Clinical applications and challenges for development.Adv. Drug Deliv. Rev.611117–1120. 10.1016/j.addr.2009.08.001
233
VineisP.SchatzkinA.PotterJ. D. (2010). Models of carcinogenesis: an overview.Carcinogenesis311703–1709. 10.1093/carcin/bgq087
234
WangY.WuD.WuG.WuJ.LuS.LoJ.et al (2020). Metastasis-on-a-chip mimicking the progression of kidney cancer in the liver for predicting treatment efficacy.Theranostics10300–311. 10.7150/thno.38736
235
WangZ.ButnerJ. D.CristiniV.DeisboeckT. S. (2015). Integrated PK-PD and agent-based modeling in oncology.J. Pharmacokinet. Pharmacodyn.42179–189. 10.1007/s10928-015-9403-7
236
WangZ.DeisboeckT. S. (2014). Mathematical modeling in cancer drug discovery.Drug Discov. Today19145–150. 10.1016/j.drudis.2013.06.015
237
WangZ.KerkettaR.ChuangY. L.DograP.ButnerJ. D.BrocatoT. A.et al (2016). Theory and experimental validation of a spatio-temporal model of chemotherapy transport to enhance tumor cell kill.PLoS Comput. Biol.12:e1004969. 10.1371/journal.pcbi.1004969
238
WareM. J.ColbertK.KeshishianV.HoJ.CorrS. J.CurleyS. A.et al (2016a). Generation of Homogenous Three-Dimensional Pancreatic Cancer Cell Spheroids Using an Improved Hanging Drop Technique.Tissue Eng. Part C Methods22312–321. 10.1089/ten.tec.2015.0280
239
WareM. J.KeshishianV.LawJ. J.HoJ. C.FavelaC. A.ReesP.et al (2016b). Generation of an in vitro 3D PDAC stroma rich spheroid model.Biomaterials108129–142. 10.1016/j.biomaterials.2016.08.041
240
WhiteD.CoombeD.RezaniaV.TuszynskiJ. (2016). Building a 3D virtual liver: methods for simulating blood flow and hepatic clearance on 3d structures.PLoS One11:e0162215. 10.1371/journal.pone.0162215
241
WilliamsC. B.YehE. S.SoloffA. C. (2016). Tumor-associated macrophages: unwitting accomplices in breast cancer malignancy.NPJ Breast. Cancer2:15025.
242
WongC. H.SiahK. W.LoA. W. (2019). Estimation of clinical trial success rates and related parameters.Biostatistics20273–286. 10.1093/biostatistics/kxx069
243
WuM.FrieboesH. B.ChaplainM. A.McDougallS. R.CristiniV.LowengrubJ. S. (2014). The effect of interstitial pressure on therapeutic agent transport: coupling with the tumor blood and lymphatic vascular systems.J. Theor. Biol.355194–207. 10.1016/j.jtbi.2014.04.012
244
YamamotoK. N.YachidaS.NakamuraA.NiidaA.OshimaM.DeS.et al (2017). Personalized management of pancreatic ductal adenocarcinoma patients through computational modeling.Cancer Res.773325–3335. 10.1158/0008-5472.can-16-1208
245
YokoiK.GodinB.ObornC. J.AlexanderJ. F.LiuX.FidlerI. J.et al (2013). Porous silicon nanocarriers for dual targeting tumor associated endothelial cells and macrophages in stroma of orthotopic human pancreatic cancers.Cancer Lett.334319–327. 10.1016/j.canlet.2012.09.001
246
YuanX.ZhangJ.LiD.MaoY.MoF.DuW.et al (2017). Prognostic significance of tumor-associated macrophages in ovarian cancer: A meta-analysis.Gynecol. Oncol.147181–187. 10.1016/j.ygyno.2017.07.007
247
ZamboniW. C.TorchilinV.PatriA. K.HrkachJ.SternS.LeeR.et al (2012). Best practices in cancer nanotechnology: perspective from NCI nanotechnology alliance.Clin. Cancer Res.183229–3241. 10.1158/1078-0432.ccr-11-2938
248
ZempolichK. A.DubeauL. (1998). Telomerase and cancer management: silver bullet or fool’s gold?Gynecol. Oncol.68143–144. 10.1006/gyno.1998.4948
249
ZhangG.HeP.TanH.BudhuA.GaedckeJ.GhadimiB. M.et al (2013). Integration of metabolomics and transcriptomics revealed a fatty acid network exerting growth inhibitory effects in human pancreatic cancer.Clin. Cancer Res.194983–4993. 10.1158/1078-0432.ccr-13-0209
250
ZhangM.HeY.SunX.LiQ.WangW.ZhaoA.et al (2014). A high M1/M2 ratio of tumor-associated macrophages is associated with extended survival in ovarian cancer patients.J. Ovarian Res.7:19. 10.1186/1757-2215-7-19
251
ZhangP.BrusicV. (2014). Mathematical modeling for novel cancer drug discovery and development.Expert. Opin. Drug Discov.91133–1150. 10.1517/17460441.2014.941351
252
ZhangY.HuangW. J.YangQ. R.ZhangH. D.ZhuX. J.ZengM.et al (2020). Cryopreserved biopsy tissues of rectal cancer liver metastasis for assessment of anticancer drug response in vitro and in vivo.Oncol. Rep.43405–414.
253
ZhangY. S.DuchampM.OkluR.EllisenL. W.LangerR.KhademhosseiniA. (2016). Bioprinting the Cancer Microenvironment.ACS Biomater. Sci. Eng.21710–1721. 10.1021/acsbiomaterials.6b00246
254
ZhaoY.YaoR.OuyangL.DingH.ZhangT.ZhangK.et al (2014). Three-dimensional printing of Hela cells for cervical tumor model in vitro.Biofabrication6:035001. 10.1088/1758-5082/6/3/035001
255
ZhouX.LiuC.ZhaoX.WangX. (2016). A 3D bioprinting liver tumor model for drug screening.World J. Pharm. Pharm. Sci.5196–213.
256
ZingerA.KorenL.AdirO.PoleyM.AlyanM.YaariZ.et al (2019). Collagenase nanoparticles enhance the penetration of drugs into pancreatic tumors.ACS nano1311008–11021. 10.1021/acsnano.9b02395
257
ZoliW.RicottiL.TeseiA.BarzantiF.AmadoriD. (2001). In vitro preclinical models for a rational design of chemotherapy combinations in human tumors.Crit. Rev. Oncol. Hematol.3769–82. 10.1016/s1040-8428(00)00110-4
Summary
Keywords
liver metastasis, nanotherapy, tumor microenvironment, macrophages, mathematical modeling, computational simulation
Citation
Frieboes HB, Raghavan S and Godin B (2020) Modeling of Nanotherapy Response as a Function of the Tumor Microenvironment: Focus on Liver Metastasis. Front. Bioeng. Biotechnol. 8:1011. doi: 10.3389/fbioe.2020.01011
Received
16 June 2020
Accepted
03 August 2020
Published
19 August 2020
Volume
8 - 2020
Edited by
Paola Lecca, Free University of Bozen-Bolzano, Italy
Reviewed by
Katarzyna Anna Rejniak, Moffitt Cancer Center, United States; Viviana Salvatore, University of Bologna, Italy
Updates
Copyright
© 2020 Frieboes, Raghavan and Godin.
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: Biana Godin, BGodin@houstonmethodist.org; bianagodinv@gmail.com
This article was submitted to Nanobiotechnology, a section of the journal Frontiers in Bioengineering and Biotechnology
Disclaimer
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.