Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

FoxM1B transcriptionally regulates vascular endothelial growth factor expression and promotes the angiogenesis and growth of glioma cells.

We previously found that FoxM1B is overexpressed in human glioblastomas and that forced FoxM1B expression in anaplastic astrocytoma cells leads to the formation of highly angiogenic glioblastoma in nude mice. However, the molecular mechanisms by which FoxM1B enhances glioma angiogenesis are currently unknown. In this study, we found that vascular endothelial growth factor (VEGF) is a direct transcriptional target of FoxM1B. FoxM1B overexpression increased VEGF expression, whereas blockade of FoxM1 expression suppressed VEGF expression in glioma cells. Transfection of FoxM1 into glioma cells directly activated the VEGF promoter, and inhibition of FoxM1 expression by FoxM1 siRNA suppressed VEGF promoter activation. We identified two FoxM1-binding sites in the VEGF promoter that specifically bound to the FoxM1 protein. Mutation of these FoxM1-binding sites significantly attenuated VEGF promoter activity. Furthermore, FoxM1 overexpression increased and inhibition of FoxM1 expression suppressed the angiogenic ability of glioma cells. Finally, an immunohistochemical analysis of 59 human glioblastoma specimens also showed a significant correlation between FoxM1 overexpression and elevated VEGF expression. Our findings provide both clinical and mechanistic evidence that FoxM1 contributes to glioma progression by enhancing VEGF gene transcription and thus tumor angiogenesis.

The functional role of the tumor suppressor MMAC /PTEN in glioblastoma multiforme and prostate cancer

Investigations into the molecular basis of glioblastoma multiforme led to the identification of a putative tumor suppressor gene, MMAC/ PTEN. Initial studies implicated MMAC/PTEN in many different tumor types, and identified a protein phosphatase motif in its sequence. This project aimed to identify the biological and biochemical functions of MMAC/PTEN by transiently expressing the gene in cancer cells that lack a functional gene product. ^ Expression of MMAC/PTEN mildly suppressed the growth of U251 human glioma cells and abrogated the growth advantage mediated by overexpression of the epidermal growth factor receptor (EGFR). Immunoblotting demonstrated that MMAC/PTEN expression did not affect the phosphorylation of the EGFR itself, or the intermediates of several downstream signaling pathways. However, MMAC/PTEN expression significantly reduced the phosphorylation and catalytic activity of the proto-oncogene Akt/PKB. While Akt/PKB regulates the survival of many cell types, expression of MMAC/PTEN did not induce apoptosis in adherent U251 cells. Instead, MMAC/PTEN expression sensitized the cells to apoptosis when maintained in suspension (anoikis). As the survival of suspended cells is one of the hallmarks leading to metastasis, MMAC/PTEN expression was examined in a system in which metastasis is more clinically relevant, prostate cancer. ^ Expression of MMAC/PTEN in both LNCaP and PC3-P human prostate cancer cells specifically inhibited Akt/PKB phosphorylation. MMAC/PTEN expression in LNCaP cells resulted in a profound inhibition of growth that was significantly greater than that achieved with expression of p53. Expression of MMAC/PTEN in PC3-P cells resulted in greater growth inhibition than was observed in U251 glioma cells, but less than was observed in LNCaP cells, or upon p53 expression. To determine if MMAC/PTEN could function as a tumor suppressor in vivo, the effects of MMAC/PTEN expression on PC3-P cells implanted orthotopically in nude mice were examined. The ex-vivo expression of MMAC/PTEN did not decrease tumor incidence, but it did significantly decrease tumor size and metastasis. In-vivo expression of MMAC/PTEN in pre-established PC3-P tumors did not significantly inhibit tumor incidence or size, but did inhibit metastasis formation. ^ These studies demonstrate that MMAC/PTEN is a novel and important tumor suppressor gene, which functions to downregulate an important cell survival signaling pathway. ^

Characterization and optimization of antigen-specific T cell responses during ex vivo expansion of melanoma tumor-infiltrating lymphocytes

Treatment of metastatic melanoma with tumor reactive T cells (adoptive T cell therapy, ACT) is a promising approach associated with a high clinical response rate. However, further optimization of this treatment modality is required to increase the clinical response after this therapy. ACT in melanoma involves an initial phase (pre-REP) of tumor-infiltrating lymphocyte (TIL) expansion ex vivo from tumor isolates followed by a second phase, “rapid expansion protocol” (REP) generating the billions of cells used as the TIL infusion product. The main question addressed in this thesis was how the currently used REP affected the responsiveness of the CD8+ T cells to defined melanoma antigens. We hypothesized that the REP drives the TIL to further differentiate and become hyporesponsive to antigen restimulation, therefore, proper cytokine treatment or other ways to expand TIL is required to improve upon this outcome. We evaluated the response of CD8+ TIL to melanoma antigen restimulation using MART-1 peptide-pulsed mature DC in vitro. Post-REP TILs were mostly hypo-responsive with poor proliferation and higher apoptosis. Phenotypic analysis revealed that the expression of CD28 was significantly reduced in post-REP TILs. By sorting experiment and microarray analysis, we confirmed that the few CD28+ post-REP TILs had superior survival capacity and proliferated after restimulation. We then went on to investigate methods to maintain CD28 expression during the REP and improve TIL responsiveness. Firstly, IL-15 and IL-21 were found to synergize in maintaining TIL CD28 expression and antigenic responsiveness during REP. Secondly, we found IL-15 was superior as compared to IL-2 in supporting the long-term expansion of antigen-specific CD8+ TIL after restimulation. These results suggest that current expansion protocols used for adoptive T-cell therapy in melanoma yield largely hyporesponsive products containing CD8+ T cells unable to respond in vivo to re-stimulation with antigen. A modification of our current approaches by using IL-15+IL-21 as supporting cytokines in the REP, or/and administration of IL-15 instead of IL-2 after TIL infusion, may enhance the anti-tumor efficacy and long-term persistence of infused T cells in vivo.

THE ROLE OF TYROSINE PHOSPHORYLATION IN THE FUNCTIONS OF THE TUMOR SUPPRESSOR GPRC5A

The retinoic acid inducible G protein coupled receptor family C group 5 type A (GPRC5A) is expressed preferentially in normal lung tissue but its expression is suppressed in the majority of human non-small cell lung cancer cell lines and tissues. This differential expression has led to the idea that GPRC5A is a potential tumor suppressor. This notion was supported by the finding that mice with a deletion of the Gprc5a gene develop spontaneous lung tumors. However, there are various tumor cell lines and tissue samples, including lung, that exhibit higher GPRC5A expression than normal tissues and some reports by other groups that GPRC5A transfection increased cell growth and colony formation. Obviously, GPRC5A has failed to suppress the development of the tumors and the growth of the cell lines where its expression is not suppressed. Since no mutations were detected in the coding sequence of GPRC5A in 20 NSCLC cell lines, it’s possible that GPRC5A acts as a tumor suppressor in the context of some cells but not in others. Alternatively, we raised the hypothesis that the GPRC5A protein may be inactivated by posttranslational modification(s) such as phosphorylation. It is well established that Serine/Threonine phosphorylation of G protein coupled receptors leads to their desensitization and in a few cases Tyrosine phosphorylation of GPCRs has been linked to internalization. Others reported that GPRC5A can undergo tyrosine phosphorylation in the cytoplasmic domain after treatment of normal human mammary epithelial cells (HMECs) with epidermal growth factor (EGF) or Heregulin. This suggested that GPRC5A is a substrate of EGFR. Therefore, we hypothesized that tyrosine phosphorylation of GPRC5A by activation of EGFR signaling may lead to its inactivation. To test this hypothesis, we transfected human embryo kidney (HEK) 293 cells with GPRC5A and EGFR expression vectors and confirmed that GPRC5A can be tyrosine phosphorylated after activation of EGFR by EGF. Further, we found that EGFR and GPRC5A can interact either directly or through other proteins and that inhibition of the EGFR kinase activity decreased the phosphorylation of GPRA5A and the interaction between GPRC5A and EGFR. In c-terminal of GPRC5A, There are four tyrosine residues Y317, Y320, Y347, Y350. We prepared GPRC5A mutants in which all four tyrosine residues had been replaced by phenylalanine (mutant 4F) or each individual Tyr residue was replaced by Phe and found that Y317 is the major site for EGFR mediated phosphorylation in the HEK293T cell line. We also found that EGF can induce GPRC5A internalization both in H1792 transient and stable cell lines. EGF also partially inactivates the suppressive function of GPRC5A on cell invasion activity and anchorage-independent growth ability of H1792 stable cell lines. These finding support our hypothesis that GPRC5A may be inactivated by posttranslational modification- tyrosine phosphorylation.

3D visualization of HIV virions by cryoelectron tomography.

The structure of the human immunodeficiency virus (HIV) and some of its components have been difficult to study in three-dimensions (3D) primarily because of their intrinsic structural variability. Recent advances in cryoelectron tomography (cryo-ET) have provided a new approach for determining the 3D structures of the intact virus, the HIV capsid, and the envelope glycoproteins located on the viral surface. A number of cryo-ET procedures related to specimen preservation, data collection, and image processing are presented in this chapter. The techniques described herein are well suited for determining the ultrastructure of bacterial and viral pathogens and their associated molecular machines in situ at nanometer resolution.

Structural equation modeling of the medical outcome study: Short Form 36 (SF-36)

The factorial validity of the SF-36 was evaluated using confirmatory factor analysis (CFA) methods, structural equation modeling (SEM), and multigroup structural equation modeling (MSEM). First, the measurement and structural model of the hypothesized SF-36 was explicated. Second, the model was tested for the validity of a second-order factorial structure, upon evidence of model misfit, determined the best-fitting model, and tested the validity of the best-fitting model on a second random sample from the same population. Third, the best-fitting model was tested for invariance of the factorial structure across race, age, and educational subgroups using MSEM.^ The findings support the second-order factorial structure of the SF-36 as proposed by Ware and Sherbourne (1992). However, the results suggest that: (a) Mental Health and Physical Health covary; (b) general mental health cross-loads onto Physical Health; (c) general health perception loads onto Mental Health instead of Physical Health; (d) many of the error terms are correlated; and (e) the physical function scale is not reliable across these two samples. This hierarchical factor pattern was replicated across both samples of health care workers, suggesting that the post hoc model fitting was not data specific. Subgroup analysis suggests that the physical function scale is not reliable across the "age" or "education" subgroups and that the general mental health scale path from Mental Health is not reliable across the "white/nonwhite" or "education" subgroups.^ The importance of this study is in the use of SEM and MSEM in evaluating sample data from the use of the SF-36. These methods are uniquely suited to the analysis of latent variable structures and are widely used in other fields. The use of latent variable models for self reported outcome measures has become widespread, and should now be applied to medical outcomes research. Invariance testing is superior to mean scores or summary scores when evaluating differences between groups. From a practical, as well as, psychometric perspective, it seems imperative that construct validity research related to the SF-36 establish whether this same hierarchical structure and invariance holds for other populations.^ This project is presented as three articles to be submitted for publication. ^

BIOLOGICAL AND BIOCHEMICAL CHARACTERIZATION OF MURINE-TUMOR SPECIFIC TRANSPLANTATION ANTIGENS (ISOTACHOPHORESIS, IMMUNOGENICITY)

Tumor-specific transplantation antigens (TSTA) are individually distinct neoantigens expressed on the cells of chemically-induced neoplasms. TSTA are operationally defined by immunization of syngeneic mice against challenge with viable tumor cells. Immunization with cell surface or extracted TSTA induces specific resistance to transplanted tumor cells. The biological and biochemical nature of TSTA was investigated in the 3-methylcholanthrene-induced fibrosarcomas of female C3H/HeJ mice, MCA-F and MCA-D. Tumor cell suspensions were extracted by treatment with 3M KCl or 2.5% butanol solutions and the TSTA was partially purified by preparative isoelectric focusing. The isoelectric pH of TSTA purified from 3M KCl extracts was 5.8-6.0, and from butanol extracts was 6.4-6.6. Whereas immunization with 10('5) and 10('6) irradiated tumor cells induces complete rejection of tumor cell challenge over a two-fold-log dose range, immunization with ug quantities within a one-fold-log dose range of extracted TSTA induces only partial resistance to tumor challenge. Reduced immunogenicity of extracted TSTA is hypothesized to result from immunization of mice with insufficiently purified TSTA preparations. The hypothesis predicts that immunization with highly purified TSTA, free from interfering substances, induces complete rejection of tumor challenge over a broad dose range. To test the hypothesis preparative isotachophoresis (pITP) was used to purify TSTA from electrofocused TSTA fractions. Significant purification was achieved, as immunization with 15 pg to 1.5 ug (5 logs) of pITP-purified TSTA extracted from the MCA-F, or with 1 pg to 10 ng (4 logs) of TSTA from the MCA-D tumor induced specific resistance to tumor challenge. Despite 50,000 fold purification of TSTA, immunization induced partial, not complete, rejection of transplanted tumor cells. This suggests a clear dissociation of the immunogenicity and purification of extracted TSTA, indicating that the induction of partial immunity to tumor challenge is an intrinsic property of extracted TSTA.^

The localization and identification of candidate tumor suppressor genes involved in prostate cancer

Frequent loss of heterozygosity (LOH) at specific chromosomal regions are highly associated with the inactivation of tumor suppressor genes (TSGs) (Weinberg, 1991; Bishop, 1989). Chromosome 8p is the most frequently reported site of LOH (∼60%) in prostate cancer (PC), suggesting that there may be inactivated TSG(s) involved in PC on chromosome 8p. (Bergerheim et. al., 1991; Kagan et. al., 1995). In order to identify the smallest common regions of frequent LOH (SCLs) on chromosome 8, we screened 52 PC patient/tumor samples with 39 polymorphic markers in successive screenings. In the course of refining the SCLs, we identified 3 tumors with >6 Mb homozygous deletions (HZDs) at 8p22 and 8p21, suggesting the presence of candidate TSGs at both loci. These HZDs spanned the two SCLs at 8p22 (46%) and 8p21 (45%). The SCLs were narrowed to 3.2 cM at 8p22 and less than 3 cM at 8p21. ^ In order to identify candidate TSGs within the SCLs on 8p, two approaches were used. In the candidate gene approach, thirty genes that mapped to the SCLs were evaluated for expression in normal prostate and in PC cell lines. One of the candidate genes, Clusterin, showed decreased expression in 4/7 (57%) prostate cancer cell lines by Northern blot analysis. Clusterin will be further examined as a candidate TSG. ^ The second approach involved utilizing subtractive hybridization and hybrid affinity capture to generate pools of expressed sequence tags (ESTs) enriched for genes that are downregulated or deleted in PC and that map to specific regions of interest. We took advantage of a prostate cancer cell line (PC3) with a known HZD of a candidate TSG, CTNNA1 on 5q31, to develop and validate a model system. We then developed subtracted libraries enriched for 8p22 and 8p21 ESTs by this method, using two cell lines, MDAPCa-2b and PC3. The ESTs were cloned, and 40 were sequenced and evaluated for expression in normal prostate and PC cell lines. Three ESTs from the subtracted libraries, C2, C17 and F12, showed decreased expression in 29–57% of the prostate tumor cell lines studied, and will be further examined as candidate TSGs. ^

Towards the identification of a tumor suppressor gene on chromosome 3p12: Physical mapping and candidate gene screening

Renal cell carcinoma (RCC) is the most common malignant tumor of the kidney. Characterization of RCC tumors indicates that the most frequent genetic event associated with the initiation of tumor formation involves a loss of heterozygosity or cytogenetic aberration on the short arm of human chromosome 3. A tumor suppressor locus Nonpapillary Renal Carcinoma-1 (NRC-1, OMIM ID 604442) has been previously mapped to a 5–7 cM region on chromosome 3p12 and shown to induce rapid tumor cell death in vivo, as demonstrated by functional complementation experiments. ^ To identify the gene that accounts for the tumor suppressor activities of NRC-1, fine-scale physical mapping was conducted with a novel real-time quantitative PCR based method developed in this study. As a result, NRC-1 was mapped within a 4.6-Mb region defined by two unique sequences within UniGene clusters Hs.41407 and Hs.371835 (78,545Kb–83,172Kb in the NCBI build 31 physical map). The involvement of a putative tumor suppressor gene Robo1/Dutt1 was excluded as a candidate for NRC-1. Furthermore, a transcript map containing eleven candidate genes was established for the 4.6-Mb region. Analyses of gene expression patterns with real-time quantitative RT-PCR assays showed that one of the eleven candidate genes in the interval (TSGc28) is down-regulated in 15 out of 20 tumor samples compared with matched normal samples. Three exons of this gene have been identified by RACE experiments, although additional exon(s) seem to exist. Further gene characterization and functional studies are required to confirm the gene as a true tumor suppressor gene. ^ To study the cellular functions of NRC-1, gene expression profiles of three tumor suppressive microcell hybrids, each containing a functional copy of NRC-1, were compared with those of the corresponding parental tumor cell lines using 16K oligonucleotide microarrays. Differentially expressed genes were identified. Analyses based on the Gene Ontology showed that introduction of NRC-1 into tumor cell lines activates genes in multiple cellular pathways, including cell cycle, signal transduction, cytokines and stress response. NRC-1 is likely to induce cell growth arrest indirectly through WEE1. ^

Mechanism of anti-tumor activity of adenovirus type 5 E1A in ovarian cancer

The adenovirus type 5 E1A gene was originally developed as a gene therapy to inhibit tumorigenicity of HER-2-overexpressing cells by transcriptional downregulation of HER-2. Our goal is to improve the overall efficacy of E1A gene therapy. To achieve this goal, we have conducted two preclinical experiments. ^ First, we hypothesized that Bcl-2 overexpressing ovarian cancer is resistant to E1A gene therapy. This hypothesis is based on that the 19 kDa protein product of the adenoviral E1B gene which is homologous to Bcl-2 inhibits E1A-induced apoptosis. Treating high Bcl-2-xpressing cells with E1A in combination with an antisense oligonucleotide to Bcl-2 (Bcl-2-ASO) resulted in a significant decrease in cell viability due to an increased rate of apoptosis relative to cells treated with E1A alone. In an ovarian cancer xenograft model, mice implanted with low HER-2, high Bcl-2 cells, treated with E1A plus Bcl-2-ASO led to prolonged survival. Bcl-2 thus may serve as a predictive molecular marker enabling us to select patients with ovarian cancer who will benefit significantly from E1A gene therapy. ^ Second, we elucidated the molecular mechanism governing the anti-tumor effect of E1A in ovarian cancer to identify a more potent tumor suppressor gene. We identified PEA-15 (phospho-protein enriched in astrocytes) upregulated in E1A transfected low HER-2-expressing OVCAR-3 ovarian cancer cell, which showed decreased cell proliferation. PEA-15 moved ERK from the nucleus to the cytoplasm and inhibited ERK-dependent transcription and proliferation. Using small interfering RNA to knock down PEA-15 expression in OVCAR-3 cells made to constitutively express E1A resulted in accumulation of phosphoERK in the nucleus, an increase in Elk-1 activity, DNA synthesis, and anchorage-independent growth. PEA-15 also independently suppressed colony formation in some breast and ovarian cancer cell lines in which E1A is known to have anti-tumor activity. We conclude that the anti-tumor activity of E1A depends on PEA-15. ^ In summary, (1) Bcl-2 may serve as a predictive molecular marker of E1A gene therapy, allowing us to select patients and improve efficacy of E1A gene therapy. (2) PEA-15 was identified as a component of the molecular mechanism governing the anti-tumor activity of E1A in ovarian cancer, (3) PEA-15 may be developed as a novel therapeutic gene. ^

The role of functional tumor suppressor 15-lipoxygenase 2 (15-LOX2) in normal human prostate epithelial cells and prostate cancer development

15-Lipoxygenase 2 (15-LOX2) is a recently cloned human lipoxygenase that shows tissue-restricted expression in prostate, lung, skin, and cornea. The protein level and enzymatic activity of 15-LOX2 have been shown to be down-regulated in prostate cancers compared with normal and benign prostate tissues. We report the cloning and functional characterization of 15-LOX2 and its three splice variants (termed 15-LOX2sv-a, 15-LOX2sv-b, and 15-LOX2sv-c) from primary prostate epithelial (NHP) cells. Western blotting with multiple NHP cell strains and prostate cancer (PCa) cell lines reveals that the expression of 15-LOX2 is lost in all PCa cell lines, accompanied by decreased enzymatic activity. 15-LOX2 is expressed at multiple subcellular locations, including cytoplasm, cytoskeleton, cell-cell border, and nucleus. Surprisingly, the three splice variants of 15-LOX2 are mostly excluded from the nucleus. To elucidate the relationship between nuclear localization, enzymatic activity, and tumor suppressive functions, we established PCa cell clones stably expressing 15-LOX2 or 15-LOX2sv-b. The 15-LOX2 clones express 15-LOX2 in the nuclei and possess robust enzymatic activity, whereas 15-LOX2sv-b clones show neither nuclear protein localization nor arachidonic acid-metabolizing activity. Interestingly, both 15-LOX2- and 15-LOX2sv-b-stable clones proliferate much slower in vitro when compared with control clones. When orthotopically implanted in nude mouse prostate, both 15-LOX2 and 15-LOX2sv-b suppress PC3 tumor growth in vivo. Finally, cultured NHP cells lose the expression of putative stem/progenitor cell markers, slow down in proliferation, and enter senescence. Several pieces of evidence implicate 15-LOX2 plays a role in replicative senescence of NHP cells: (1) promoter activity and the mRNA and protein levels of 15-LOX2 and its splice variants are upregulated in serially passaged NHP cells, which precede replicative senescence and occur in a cell-autonomous manner; (2) PCa cells stably expressing 15-LOX2 or 15-LOX2sv-b show a passage-related senescence-like phenotype; (3) enforced expression of 15-LOX2 or 15-LOX2sv-b in young NHP cells induce partial cell-cycle arrest and senescence-like phenotypes. Together, these results suggest that 15-LOX2 suppress prostate tumor development and do not necessarily depend on arachidonic acid-metabolizing activity and nuclear localization. Also, 15-LOX2 may serve as an endogenous prostate senescence gene and its tumor-suppressing functions might be associated with its ability to induce cell senescence. ^

Bayesian joint modeling of longitudinal and survival data

The joint modeling of longitudinal and survival data is a new approach to many applications such as HIV, cancer vaccine trials and quality of life studies. There are recent developments of the methodologies with respect to each of the components of the joint model as well as statistical processes that link them together. Among these, second order polynomial random effect models and linear mixed effects models are the most commonly used for the longitudinal trajectory function. In this study, we first relax the parametric constraints for polynomial random effect models by using Dirichlet process priors, then three longitudinal markers rather than only one marker are considered in one joint model. Second, we use a linear mixed effect model for the longitudinal process in a joint model analyzing the three markers. In this research these methods were applied to the Primary Biliary Cirrhosis sequential data, which were collected from a clinical trial of primary biliary cirrhosis (PBC) of the liver. This trial was conducted between 1974 and 1984 at the Mayo Clinic. The effects of three longitudinal markers (1) Total Serum Bilirubin, (2) Serum Albumin and (3) Serum Glutamic-Oxaloacetic transaminase (SGOT) on patients' survival were investigated. Proportion of treatment effect will also be studied using the proposed joint modeling approaches. ^ Based on the results, we conclude that the proposed modeling approaches yield better fit to the data and give less biased parameter estimates for these trajectory functions than previous methods. Model fit is also improved after considering three longitudinal markers instead of one marker only. The results from analysis of proportion of treatment effects from these joint models indicate same conclusion as that from the final model of Fleming and Harrington (1991), which is Bilirubin and Albumin together has stronger impact in predicting patients' survival and as a surrogate endpoints for treatment. ^

Preclinical modeling of multi-drug cancer therapies

Anticancer drugs typically are administered in the clinic in the form of mixtures, sometimes called combinations. Only in rare cases, however, are mixtures approved as drugs. Rather, research on mixtures tends to occur after single drugs have been approved. The goal of this research project was to develop modeling approaches that would encourage rational preclinical mixture design. To this end, a series of models were developed. First, several QSAR classification models were constructed to predict the cytotoxicity, oral clearance, and acute systemic toxicity of drugs. The QSAR models were applied to a set of over 115,000 natural compounds in order to identify promising ones for testing in mixtures. Second, an improved method was developed to assess synergistic, antagonistic, and additive effects between drugs in a mixture. This method, dubbed the MixLow method, is similar to the Median-Effect method, the de facto standard for assessing drug interactions. The primary difference between the two is that the MixLow method uses a nonlinear mixed-effects model to estimate parameters of concentration-effect curves, rather than an ordinary least squares procedure. Parameter estimators produced by the MixLow method were more precise than those produced by the Median-Effect Method, and coverage of Loewe index confidence intervals was superior. Third, a model was developed to predict drug interactions based on scores obtained from virtual docking experiments. This represents a novel approach for modeling drug mixtures and was more useful for the data modeled here than competing approaches. The model was applied to cytotoxicity data for 45 mixtures, each composed of up to 10 selected drugs. One drug, doxorubicin, was a standard chemotherapy agent and the others were well-known natural compounds including curcumin, EGCG, quercetin, and rhein. Predictions of synergism/antagonism were made for all possible fixed-ratio mixtures, cytotoxicities of the 10 best-scoring mixtures were tested, and drug interactions were assessed. Predicted and observed responses were highly correlated (r2 = 0.83). Results suggested that some mixtures allowed up to an 11-fold reduction of doxorubicin concentrations without sacrificing efficacy. Taken together, the models developed in this project present a general approach to rational design of mixtures during preclinical drug development. ^