Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Targeted cancer diagnosis with radiolabeled endostatin

Nuclear imaging is used for non-invasive detection, staging and therapeutic monitoring of tumors through the use of radiolabeled probes. Generally, these probes are used for applications in which they provide passive, non-specific information about the target. Therefore, there is a significant need for actively-targeted radioactive probes to provide functional information about the site of interest. This study examined endostatin, an endogenous inhibitor of tumor angiogenesis, which has affinity for tumor vasculature. The major objective of this study was to develop radiolabeled analogues of endostatin through novel chemical and radiochemical syntheses, and to determine their usefulness for tumor imaging using in vitro and in vivo models of vascular, mammary and prostate tumor cells. I hypothesize that this binding will allow for a non-invasive approach to detection of tumor angiogenesis, and such detection can be used for therapeutic monitoring to determine the efficacy of anti-angiogenic therapy. ^ The data showed that endostatin could be successfully conjugated to the bifunctional chelator ethylenedicysteine (EC), and radiolabeled with technetium-99m and gallium-68, providing a unique opportunity to use a single precursor for both nuclear imaging modalities: 99mTc for single photon emission computed tomography and 68Ga for positron emission tomography, respectively. Both radiolabeled analogues showed increased binding as a function of time in human umbilical vein endothelial cells and mammary and prostate tumor cells. Binding could be blocked in a dose-dependent manner by unlabeled endostatin implying the presence of endostatin receptors on both vascular and tumor cells. Animal biodistribution studies demonstrated that both analogues were stable in vivo, showed typical reticuloendothelial and renal excretion and produced favorable absorbed organ doses for application in humans. The imaging data provide evidence that the compounds quantitate tumor volumes with clinically-useful tumor-to-nontumor ratios, and can be used for treatment follow-up to depict changes occurring at the vascular and cellular levels. ^ Two novel endostatin analogues were developed and demonstrated interaction with vascular and tumor cells. Both can be incorporated into existing nuclear imaging platforms allowing for potential wide-spread clinical benefit as well as serving as a diagnostic tool for elucidation of the mechanism of action of endostatin. ^

Novel Imaging-Based Techniques Reveal a Role for PD-1/PD-L1 in Tumor Immune Surveillance in the Lung

The binding of immune inhibitory receptor Programmed Death 1 (PD-1) on T cells to its ligand PD-L1 has been implicated as a major contributor to tumor induced immune suppression. Clinical trials of PD-L1 blockade have proven effective in unleashing therapeutic anti-tumor immune responses in a subset of patients with advanced melanoma, yet current response rates are low for reasons that remain unclear. Hypothesizing that the PD-1/PD-L1 pathway regulates T cell surveillance within the tumor microenvironment, we employed intravital microscopy to investigate the in vivo impact of PD-L1 blocking antibody upon tumor-associated immune cell migration. However, current analytical methods of intravital dynamic microscopy data lack the ability to identify cellular targets of T cell interactions in vivo, a crucial means for discovering which interactions are modulated by therapeutic intervention. By developing novel imaging techniques that allowed us to better analyze tumor progression and T cell dynamics in the microenvironment; we were able to explore the impact of PD-L1 blockade upon the migratory properties of tumor-associated immune cells, including T cells and antigen presenting cells, in lung tumor progression. Our results demonstrate that early changes in tumor morphology may be indicative of responsiveness to anti-PD-L1 therapy. We show that immune cells in the tumor microenvironment as well as tumors themselves express PD-L1, but immune phenotype alone is not a predictive marker of effective anti-tumor responses. Through a novel method in which we quantify T cell interactions, we show that T cells are largely engaged in interactions with dendritic cells in the tumor microenvironment. Additionally, we show that during PD-L1 blockade, non-activated T cells are recruited in greater numbers into the tumor microenvironment and engage more preferentially with dendritic cells. We further show that during PD-L1 blockade, activated T cells engage in more confined, immune synapse-like interactions with dendritic cells, as opposed to more dynamic, kinapse-like interactions with dendritic cells when PD-L1 is free to bind its receptor. By advancing the contextual analysis of anti-tumor immune surveillance in vivo, this study implicates the interaction between T cells and tumor-associated dendritic cells as a possible modulator in targeting PD-L1 for anti-tumor immunotherapy.

Dynamic contrast agent -enhanced magnetic resonance imaging assessment of tumor microvasculature

Dynamic contrast agent-enhanced magnetic resonance imaging (DCE MRI) data, when analyzed with the appropriate pharmacokinetic models, have been shown to provide quantitative estimates of microvascular parameters important in characterizing the angiogenic activity of malignant tissue. These parameters consist of the whole blood volume per unit volume of tissue, v b, transport constant from the plasma to the extravascular, extracellular space (EES), k1 and the transport constant from the EES to the plasma, k2. Parameters vb and k1 are expected to correlate with microvascular density (MVD) and vascular permeability, respectively, which have been suggested to serve as surrogate markers for angiogenesis. In addition to being a marker for angiogenesis, vascular permeability is also useful in estimating tumor penetration potential of chemotherapeutic agents. ^ Histological measurements of the intratumoral microvascular environment are limited by their invasiveness and susceptibility to sampling errors. Also, MVD and vascular permeability, while useful for characterizing tumors at a single time point, have shown less utility in longitudinal studies, particularly when used to monitor the efficacy of antiangiogenic and traditional chemotherapeutic agents. These limitations led to a search for a non-invasive means of characterizing the microvascular environment of an entire tumor. ^ The overall goal of this project was to determine the utility of DCE MRI for monitoring the effect of antiangiogenic agents. Further applications of a validated DCE MRI technique include in vivo measurements of tumor microvascular characteristics to aid in determining prognosis at presentation and in estimating drug penetration. DCE MRI data were generated using single- and dual-tracer pharmacokinetic models with different molecular-weight contrast agents. The resulting pharmacokinetic parameters were compared to immunohistochemical measurements. The model and contrast agent combination yielding the best correlation between the pharmacokinetic parameters and histological measures was further evaluated in a longitudinal study to evaluate the efficacy of DCE MRI in monitoring the intratumoral microvascular environment following antiangiogenic treatment. ^

Patient Characteristics and Outcomes of Gastrointestinal Stromal Tumor Patients Undergoing Imatinib Plasma Level Testing

Although gastrointestinal stromal tumor (GIST) is effectively treated with imatinib, there are a number of clinical challenges in the optimal treatment of these patients. The plasma steady-state trough level of imatinib has been proposed to correlate with clinical outcome. Plasma imatinib level may be affected by a number of patient characteristics. Additionally, the ideal plasma trough concentration of imatinib is likely to vary based on the KIT genotype (genotype determines imatinib binding affinity) of the individual patient. Patients’ genotype or plasma imatinib level may influence the type and duration of response that is appreciable by clinical evaluation. The objectives of this study were to determine effects of genotype on the type of response appreciable by current imaging criteria, to determine the distribution of plasma imatinib levels in patients with GIST, to determine factors that correlate with plasma imatinib level, to determine the incremental effects of imatinib dose escalation; and to explore the median plasma levels and outcomes of patients with various KIT mutations. We therefore obtained KIT mutation information and analyzed CT response for size and density measurement of GISTs at baseline and within the first four moths of imatinib treatment. In 126 patients with metastatic/unresectable disease, the KIT genotype of patients’ tumor was significantly associated with unique response characteristics measurable by CT. Furthermore, hepatic and peritoneal metastases differed in their response characteristics. A subgroup of patients with KIT exon 9 mutation, who received higher doses of imatinib and experienced higher trough imatinib levels, experienced improved progression-free survival similar to that of KIT exon 11 patients. Therefore, we have found that imatinib plasma levels were higher in patients with elevated Aspartate amino transferase, were women, were older, or were being treated concomitantly with CYP450 substrate drugs. As expected, CYP450 inducers correlated with a lower plasma imatinib levels in GIST patients. Renal metabolism of imatinib accounts for <10%, so it was not included in the analysis but may affect covariates. Interestingly, there was a trend for low imatinib levels and inferior progression-free survival in patients who had undergone complete gastrectomy. Patients with KIT exon 9 mutation in our cohort received higher imatinib doses, experienced higher trough imatinib levels, and experienced a PFS similar to that of KIT exon 11 patients. In conclusion, imatinib plasma levels are influenced by a number of patient characteristics. The optimal imatinib plasma level for individual patients is not known but is an area of intense investigation. Our study confirms patients with KIT exon 9 mutations benefit from high-dose imatinib and higher trough imatinib levels.

NOVEL AMINO ACID TRANSPORTER-TARGETED RADIOTRACERS FOR BREAST CANCER IMAGING

Breast cancer is the most common malignancy among women in the world. Its 5-year survival rate ranges from 23.4% in patients with stage IV to 98% in stage I disease, highlighting the importance of early detection and diagnosis. 18F-2-Fluoro-2-deoxy-glucose (18F-FDG), using positron emission tomography (PET), is the most common functional imaging tool for breast cancer diagnosis currently. Unfortunately, 18F-FDG-PET has several limitations such as poorly differentiating tumor tissues from inflammatory and normal brain tissues. Therefore, 18F-labeled amino acid-based radiotracers have been reported as an alternative, which is based on the fact that tumor cells uptake and consume more amino acids to sustain their uncontrolled growth. Among those radiotracers, 18F-labeled tyrosine and its derivatives have shown high tumor uptake and great ability to differentiate tumor tissue from inflammatory sites in brain tumors and squamous cell carcinoma. They enter the tumor cells via L-type amino acid transporters (LAT), which were reported to be highly expressed in many cancer cell lines and correlate positively with tumor growth. Nevertheless, the low radiosynthesis yield and demand of an on-site cyclotron limit the use of 18F-labeled tyrosine analogues. In this study, four Technetium-99m (99mTc) labeled tyrosine/ AMT (α-methyl tyrosine)-based radiotracers were successfully synthesized and evaluated for their potentials in breast cancer imaging. In order to radiolabel tyrosine and AMT, the chelators N,N’-ethylene-di-L-cysteine (EC) and 1,4,8,11-tetra-azacyclotetradecane (N4 cyclam) were selected to coordinate 99mTc. These chelators have been reported to provide stable chelation ability with 99mTc. By using the chelator technology, the same target ligand could be labeled with different radioisotopes for various imaging modalities for tumor diagnosis, or for internal radionuclide therapy in future. Based on the in vitro and in vivo evaluation using the rat mammary tumor models, 99mTc-EC-AMT is considered as the most suitable radiotracer for breast cancer imaging overall, however, 99mTc-EC-Tyrosine will be more preferred for differential diagnosis of tumor from inflammation.

APPLICATIONS OF EPHB4 RECEPTOR SPECIFIC PEPTIDES IN TARGETED CANCER IMAGING AND THERAPY

EphB4 receptors, a member of the largest family of receptor tyrosine kinases, are found over-expressed in a variety of tumors cells including glioma cells as well as angiogenic blood vessels. Noninvasive imaging of EphB4 could potentially increase early detection rates, monitor response to therapy directed against EphB4, and improve patient outcomes. Targeted delivery of EphB4 receptor specific peptide conjugated hollow gold nanoshells (HAuNS) into tumors has great potential in cancer imaging and photothermal therapy. In this study, we developed an EphB4 specific peptide named TNYL-RAW and labeled with radioisotope 64Cu and Cy5.5 dye. We also conjugate this specific peptide with hollow gold nanoshells (HAuNS) to evaluate targeted photothermal therapy of cancers. In vitro, 64Cu-DOTA-TNYL- RAW specifically bind to CT26 and PC-3M cells but not to A549 cells. In vivo, Small-animal PET/CT clearly showed the significant uptake of 64Cu-DOTA-TNYL-RAW in CT26 and PC-3M tumors but not in A549 tumors. Furthermore, µPET/CT and near-infrared optical imaging clearly showed the uptake of the dual labeled TNYL-RAW peptide in both U251 and U87 tumors in the brains of nude mice. In U251 tumors, Cy5.5-labeled peptide can bind to EphB4-expressing tumor blood vessels and tumors cells. But in U87 models, dual labeled peptide only could bind to tumor associated blood vessels. Also, Irradiation of PC-3M and CT-26 cell treated with TNYL-PEG-HAuNS nanopatilces with near-infrared (NIR) laser resulted in selective destruction of these cells in vitro. EphB4 targeted TNYL-PEG-HAuNS showed more photothermal killing effect on CT26 tumor model than PEG-HAuNS did. In summary, tumors with overexpression of EphB4 receptors can be noninvasively visualized by micro PET/CT with 64Cu labeled or dual labeled TNYL-RAW peptide. Targeted delivery of TNYL-RAW conjugated HAuNS into tumors can greatly improve the treatment effect of photothermal therapy. The information acquired with this study should be advantageous in improving diagnostics and future applications in photothermal ablation therapy in clinical.

Elucidating the role of mesenchymal stem cell participation in the tumor microenvironment

To meet the requirements for rapid tumor growth, a complex array of non-neoplastic vascular, fibroblastic, and immune cells are recruited to the tumor microenvironment. Understanding the origin, composition, and mechanism(s) for recruitment of these stromal components will help identify areas for therapeutic intervention. Previous findings have suggested that ex-vivo expanded bone marrow-derived MSC home to the sites of tumor development, responding to inflammatory signals and can serve as effective drug delivery vehicles. Therefore, we first sought to fully assess conditions under which MSC migrate to and incorporate into inflammatory microenvironments and the consequences of modulated inflammation. MSC delivered to animals bearing inflammatory insults were monitored by bioluminescence imaging and displayed specific tropism and selective incorporation into all tumor and wound sites. These findings were consistent across routes of tumor establishment, MSC administration, and immunocompetence. MSC were then used as drug delivery vehicles, transporting Interferon β to sites of pancreatic tumors. This therapy was effective at inhibiting pancreatic tumor growth under homeostatic conditions, but inhibition was lost when inflammation was decreased with CDDO-Me combination treatment. Next, to examine the endogenous tumor microenvironment, a series of tissue transplant experiments were carried out in which tissues were genetically labeled and engrafted in recipients prior to tumor establishment. Tumors were then analyzed for markers of tumor associated fibroblasts (TAF): α-smooth muscle actin (α-SMA), nerve glia antigen 2 (NG2), fibroblast activation protein (FAP), and fibroblast specific protein (FSP) as well as endothelial marker CD31 and macrophage marker F4/80. We determined the majority of α-SMA+, NG2+ and CD31+ cells were non-bone marrow derived, while most FAP+, FSP+, and F4/80+ cells were recruited from the bone marrow. In accord, transplants of prospectively isolated BM MSC prior to tumor development indicated that these cells were recruited to the tumor microenvironment and co-expressed FAP and FSP. In contrast, fat transplant experiments revealed recruited fat derived cells co-expressed α-SMA, NG2, and CD31. These results indicate TAF are a heterogeneous population composed of subpopulations with distinct tissues of origin. These models have provided a platform upon which further investigation into tumor microenvironment composition and tests for candidate drugs can be performed. ^

Radiomics of NSCLC: Quantitative CT Image Feature Characterization and Tumor Shrinkage Prediction

Radiomics is the high-throughput extraction and analysis of quantitative image features. For non-small cell lung cancer (NSCLC) patients, radiomics can be applied to standard of care computed tomography (CT) images to improve tumor diagnosis, staging, and response assessment. The first objective of this work was to show that CT image features extracted from pre-treatment NSCLC tumors could be used to predict tumor shrinkage in response to therapy. This is important since tumor shrinkage is an important cancer treatment endpoint that is correlated with probability of disease progression and overall survival. Accurate prediction of tumor shrinkage could also lead to individually customized treatment plans. To accomplish this objective, 64 stage NSCLC patients with similar treatments were all imaged using the same CT scanner and protocol. Quantitative image features were extracted and principal component regression with simulated annealing subset selection was used to predict shrinkage. Cross validation and permutation tests were used to validate the results. The optimal model gave a strong correlation between the observed and predicted shrinkages with . The second objective of this work was to identify sets of NSCLC CT image features that are reproducible, non-redundant, and informative across multiple machines. Feature sets with these qualities are needed for NSCLC radiomics models to be robust to machine variation and spurious correlation. To accomplish this objective, test-retest CT image pairs were obtained from 56 NSCLC patients imaged on three CT machines from two institutions. For each machine, quantitative image features with concordance correlation coefficient values greater than 0.90 were considered reproducible. Multi-machine reproducible feature sets were created by taking the intersection of individual machine reproducible feature sets. Redundant features were removed through hierarchical clustering. The findings showed that image feature reproducibility and redundancy depended on both the CT machine and the CT image type (average cine 4D-CT imaging vs. end-exhale cine 4D-CT imaging vs. helical inspiratory breath-hold 3D CT). For each image type, a set of cross-machine reproducible, non-redundant, and informative image features was identified. Compared to end-exhale 4D-CT and breath-hold 3D-CT, average 4D-CT derived image features showed superior multi-machine reproducibility and are the best candidates for clinical correlation.

Development of HIF-1α/HIF-1β heterodimerization inhibitors using a novel bioluminescence reporter assay system for in vitro high throughput screening and in vivo imaging

Tumor growth often outpaces its vascularization, leading to development of a hypoxic tumor microenvironment. In response, an intracellular hypoxia survival pathway is initiated by heterodimerization of hypoxia-inducible factor (HIF)-1α and HIF-1β, which subsequently upregulates the expression of several hypoxia-inducible genes, promotes cell survival and stimulates angiogenesis in the oxygen-deprived environment. Hypoxic tumor regions are often associated with resistance to various classes of radio- or chemotherapeutic agents. Therefore, development of HIF-1α/β heterodimerization inhibitors may provide a novel approach to anti-cancer therapy. To this end, a novel approach for imaging HIF-1α/β heterodimerization in vitro and in vivo was developed in this study. Using this screening platform, we identified a promising lead candidate and further chemically derivatized the lead candidate to assess the structure-activity relationship (SAR). The most effective first generation drug inhibitors were selected and their pharmacodynamics and anti-tumor efficacy in vivo were verified by bioluminescence imaging (BLI) of HIF-1α/β heterodimerization in the xenograft tumor model. Furthermore, the first generation drug inhibitors, M-TMCP and D-TMCP, demonstrated efficacy as monotherapies, resulting in tumor growth inhibition via disruption of HIF-1 signaling-mediated tumor stromal neoangiogenesis.