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.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

PAT proteins, an ancient family of lipid droplet proteins that regulate cellular lipid stores.

The PAT family of lipid droplet proteins includes 5 members in mammals: perilipin, adipose differentiation-related protein (ADRP), tail-interacting protein of 47 kDa (TIP47), S3-12, and OXPAT. Members of this family are also present in evolutionarily distant organisms, including insects, slime molds and fungi. All PAT proteins share sequence similarity and the ability to bind intracellular lipid droplets, either constitutively or in response to metabolic stimuli, such as increased lipid flux into or out of lipid droplets. Positioned at the lipid droplet surface, PAT proteins manage access of other proteins (lipases) to the lipid esters within the lipid droplet core and can interact with cellular machinery important for lipid droplet biogenesis. Genetic variations in the gene for the best-characterized of the mammalian PAT proteins, perilipin, have been associated with metabolic phenotypes, including type 2 diabetes mellitus and obesity. In this review, we discuss how the PAT proteins regulate cellular lipid metabolism both in mammals and in model organisms.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Characterization of healing tissue in a tooth extraction socket in a rabbit model

The histology of healing in a tooth extraction socket has been described in many studies. The focus of research in bone biology and healing is now centered on molecular events that regulate repair of injured tissue. Rapid progress in cellular and molecular biology has resulted in identification of many signaling molecules (growth factors and cytokines) associated with formation and repair of skeletal tissues. Some of these include members of the transforming growth factor-β superfamily (including the bone morphogenetic proteins), fibroblast growth factors, platelet derived growth factors and insulin like growth factors. ^ Healing of a tooth extraction socket is a complex process involving tissue repair and regeneration. It involves chemotaxis of appropriate cells into the wound, transformation of undifferentiated mesenchymal cells to osteoprogenitor cells, proliferation and differentiation of committed bone forming cells, extracellular matrix synthesis, mineralization of osteoid, maturation and remodeling of bone. Current data suggests that these cellular events are precisely controlled and regulated by specific signaling molecules. A plethora of cytokines; have been identified and studied in the past two decades. Some of these like transforming growth factor beta (TGF-β), vascular endothelial growth factor (VEGF), platelet derived growth factor (PDGF) and fibroblast growth factors (FGFs) are well conserved proteins involved in the initial response to injury and repair in soft and hard tissue. ^ The purpose of this study was to characterize the spatial and temporal localization of TGF-βl, VEGF, PDGF-A, FGF-2 and BMP-2, and secretory IgA in a tooth extraction socket model, and evaluate correlation of spatial and temporal changes of these growth factors to histological events. The results of this study showed positive correlation of histological events to spatial and temporal localization of TGF-β1, BMP-2, FGF-2, PDGF-A, and VEGF in a rabbit tooth extraction model. ^

Requirements for Staphylococcus aureus in vitro cellular internalization

Staphylcoccus aureus is a prokaryotic organism capable of causing numerous superficial and severe human infections. Adhesion of S. aureus to host tissues or cells is believed to be a crucial event in S. aureus infections. Subsequently, S. aureus can seed into the bloodstream resulting in metastasis of the infection. Several reports show that S. aureus can be internalized by non-professional phagocytes, a process which has been proposed to be important in S. aureus dissemination. An intracellular residence has also been proposed to provide safe harbor to reservoirs of dormant bacteria contributing to the persistence of infection. This dissertation describes an investigation into the molecular mechanisms of S. aureus internalization into both fibroblast and epithelial cells. Bacterial requirements for internalization were found to be limited to expression of proteins that bind the extracellular matrix protein fibronectin. A previously unknown fibronectin-binding region in the S. aureus fibronectin-binding protein A was discovered after showing competitive inhibition of S. aureus internalization. This novel fibronectin-binding activity is characterized. Internalization also required cell-based factors. The presence of fibronectin and cell surface receptors of the β1 integrin class, which are known to bind and internalize fibronectin, were found to be necessary for optimal internalization of S. aureus. These results led to the conclusion that fibronectin acts as a bridge between the bacterium and integrins on the host cells. The internalization process exhibits features characteristic of integrin-mediated cell migration on fibronectin-coated surfaces. Both processes involved an active form of the β1 integrin subunit and the protein tyrosine kinase Src. Finally, a Src inhibitor previously shown to be effective in reducing osteoporosis in an in vivo rat model is capable of greatly reducing S. aureus internalization. ^

Cellular mechanisms of operant and classical conditioning

The ability to associate a predictive stimulus with a subsequent salient event (i.e., classical conditioning) and the ability to associate an expressed behavior with the consequences (i.e., operant conditioning) allow for a predictive understanding of a changing environment. Although they are operationally distinct, there has been considerable debate whether at some fundamental level classical and operant conditioning are mechanistically distinct or similar. Feeding behavior of Aplysia (i.e., biting) was chosen as the model system and was successfully conditioned with appetitive forms of both operant and classical conditioning. The neuronal circuitry responsible for feeding is well understood and is suitable for cellular analyses, thus providing for a mechanistic comparison between these two forms of associative learning. ^ Neuron B51 is part of the feeding circuitry of Aplysia and is critical for the expression of ingestive behaviors. B51 also is a locus of plasticity following both operant and classical conditioning. Both in vivo and in vitro operant conditioning increased the input resistance and the excitability of B51. No pairing-specific changes in the input resistance were observed following both in vivo and in vitro classical conditioning. However, classical conditioning decreased the excitability of B51. Thus, both operant and classical conditioning modified the threshold level for activation of neuron B51, but in opposite directions, revealing key differences in the cellular mechanisms underlying these two forms of associative learning. ^ Next, the cellular mechanisms underlying operant conditioning were investigated in more detail using a single-cell analogue. The single-cell analogue successfully recapitulated the previous in vivo and in vitro operant conditioning results by increasing the input resistance and the excitability of B51. Both PKA and PKC were necessary for operant conditioning. Dopamine appears to be the transmitter mediating the reinforcement signal in this form of conditioning. A D1 dopamine receptor antibody revealed that the D1receptor localizes to the axon hillock, which is also the region that gives the strongest response when iontophoresing dopamine. ^ The studies presented herein, thus, provide for a greater understanding of the mechanisms underlying both of these forms of associative learning and demonstrate that they likely operate through distinct cellular mechanisms. ^

SFRP4 regulation of Wnt7a signaling and cellular proliferation in the endometrium

In the endometrium, hormonal effects on epithelial cells are often elicited through stromal hormone receptors via unknown paracrine mechanisms. Several lines of evidence support the hypothesis that Wnts participate in stromal-epithelial cell communication and thus mediate hormone action. Characterization of specific Wnt signaling components in the endometrium was performed using cellular localization studies and evaluating hormone effects in a rat model. Wnt7a was expressed in the luminal epithelium, whereas the extracellular Wnt modulator, SFRP4, was localized to the endometrial stroma. SFRP4 expression is significantly decreased in endometrial carcinoma and aberrant Wnt7a signaling has been shown to cause uterine defects and contribute to the onset of disease. The specific Fzds and SFRPs that bind Wnt7a and the particular signal transduction pathway each Wnt7a-Fzd pair activates have not been identified. Additionally, the function of Wnt7a and SFRP4 in the endometrium has not been addressed. A survey of all Wnt signaling proteins expressed in the endometrium was conducted and Fzd5 and Fzd10 were identified as two receptors capable of transducing the Wnt7a signal. Biologically active recombinant Wnt7a and SFRP4 proteins were purified for quantitative biochemical studies. In Ishikawa cells, Wnt7a binding to Fzd5 activated β-catenin/canonical Wnt signaling and increased cellular proliferation. Wnt7a signaling mediated by Fzd10 induced a non-canonical/JNK-responsive pathway. SFRP4 suppressed Wnt7a action in both an autocrine and paracrine manner. Treatment with SFRP4 protein and overexpression of SFRP4 inhibited endometrial cancer cell growth and induced apoptosis in vitro. A split-eGFP complementation assay was developed to visually detect Wnt7a-Fzd interactions and subsequent pathway activation in cells. By employing a unique ELISA-based protein-protein binding technique, it was demonstrated that Wnt7a binds to SFRP4 and Fzd5 with equal nanomolar affinity. The development of these novel biological tools could lead to a better understanding of Wnt-protein interactions and the identification of new modulators of Wnt signaling. This study supports a mechanism by which the nature of the Wnt7a signal in the endometrium is dependent upon the Fzd repertoire of the cell and can be regulated by SFRP4. The potential tumor suppressor function of SFRP4 suggests it may serve as a therapeutic target for endometrial carcinoma. ^

Adenosine induced pulmonary fibrosis and the contribution of the adenosine A2B receptor to the induction of osteopontin in a mouse model of pulmonary fibrosis

Pulmonary fibrosis is a devastating and lethal lung disease with no current cure. Research into cellular signaling pathways able to modulate aspects of pulmonary inflammation and fibrosis will aid in the development of effective therapies for its treatment. Our laboratory has generated a transgenic/knockout mouse with systemic elevations in adenosine due to the partial lack of its metabolic enzyme, adenosine deaminase (ADA). These mice spontaneously develop progressive lung inflammation and severe pulmonary fibrosis suggesting that aberrant adenosine signaling is influencing the development and/or progression of the disease in these animals. These mice also show marked increases in the pro-fibrotic mediator, osteopontin (OPN), which are reversed through ADA therapy that serves to lower lung adenosine levels and ameliorate aspects of the disease. OPN is known to be regulated by intracellular signaling pathways that can be accessed through adenosine receptors, particularly the low affinity A2BR receptor, suggesting that adenosine receptor signaling may be responsible for the induction of OPN in our model. In-vitro, adenosine and the broad spectrum adenosine receptor agonist, NECA, were able to induce a 2.5-fold increase in OPN transcripts in primary alveolar macrophages. This induction was blocked through antagonism of the A2BR receptor pharmacologically, and through the deletion of the receptor subtype in these cells genetically, supporting the hypothesis that the A2BR receptor was responsible for the induction of OPN in our model. These findings demonstrate for the first time that adenosine signaling is an important modulator of pulmonary fibrosis in ADA-deficient mice and that this is in part due to signaling through the A2BR receptor which leads to the induction of the pro-fibrotic molecule, otseopontin. ^

CHARACTERIZATION AND TREATMENT OF A NOVEL MOUSE MODEL OF TSC-ASSOCIATED AUTISM

Tuberous sclerosis complex (TSC) is a dominant tumor suppressor disorder caused by mutations in either TSC1 or TSC2. The proteins of these genes form a complex to inhibit the mammalian target of rapamycin complex 1 (mTORC1), which controls protein translation and cell growth. TSC causes substantial neuropathology, often leading to autism spectrum disorders (ASDs) in up to 60% of patients. The anatomic and neurophysiologic links between these two disorders are not well understood. However, both disorders share cerebellar abnormalities. Therefore, we have characterized a novel mouse model in which the Tsc2 gene was selectively deleted from cerebellar Purkinje cells (Tsc2f/-;Cre). These mice exhibit progressive Purkinje cell degeneration. Since loss of Purkinje cells is a well-reported postmortem finding in patients with ASD, we conducted a series of behavior tests to assess if Tsc2f/-;Cre mice displayed autistic-like deficits. Using the three chambered social choice assay, we found that Tsc2f/-;Cre mice showed behavioral deficits, exhibiting no preference between a stranger mouse and an inanimate object, or between a novel and a familiar mouse. Tsc2f/-;Cre mice also demonstrated increased repetitive behavior as assessed with marble burying activity. Altogether, these results demonstrate that loss of Tsc2 in Purkinje cells in a haploinsufficient background lead to behavioral deficits that are characteristic of human autism. Therefore, Purkinje cells loss and/or dysfunction may be an important link between TSC and ASD. Additionally, we have examined some of the cellular mechanisms resulting from mutations in Tsc2 leading to Purkinje cell death. Loss of Tsc2 led to upregulation of mTORC1 and increased cell size. As a consequence of increased protein synthesis, several cellular stress pathways were upregulated. Principally, these included altered calcium signaling, oxidative stress, and ER stress. Likely as a consequence of ER stress, there was also upregulation of ubiquitin and autophagy. Excitingly, treatment with an mTORC1 inhibitor, rapamycin attenuated mTORC1 activity and prevented Purkinje cell death by reducing of calcium signaling, the ER stress response, and ubiquitin. Remarkably, rapamycin treatment also reversed the social behavior deficits, thus providing a promising potential therapy for TSC-associated ASD.

Dissection of antitumor activity by lymphokine -activated killer cells: Role of FcR+ precursors and obligatory role of tumor necrosis factor-$\alpha$ in antibody -dependent cellular cytotoxicity

Human peripheral blood lymphocytes (PBL) cultured for varying lengths of time in IL-2 are able to mediate antibody independent cellular cytotoxicity (AICC) as well as antibody dependent cellular cytotoxicity (ADCC) against a wide range of tumor targets. The objective of our study is to determine the cytotoxic potential of the subset of LAK cells involved in ADCC, the tumor recognition mechanism in ADCC, the kinetics of ADCC mediated by PBL cultured under various conditions and the role of TNF-$\alpha$ in the development and maturation of ADCC effectors in the LAK population.^ The model system in this study for ADCC used a monoclonal antibody 14G2a (IgG2a), that recognizes the GD2 epitope on human melanoma cell line, SK-Mel-1. The target recognition mechanism operative in AICC (traditionally known as lymphokine activated killing or LAK) is an acquired property of these IL-2 activated cells which confers on them the unique ability to distinguish between tumor and normal cells. This recognition probably involves the presence of a trypsin sensitive N-linked glycoprotein epitope on tumor cells. Proteolytic treatment of the tumor cells with trypsin renders them resistant to AICC by PBL cultured in IL-2. However, ADCC is unaffected. This ADCC, mediated by the relatively small population of cells that are positive for the Fc receptor for IgG (FcR), is an indication that this subset of "LAK" cells does not require the trypsin sensitive epitope on tumor cells to mediate killing. Enriching PBL for FcR+ cells markedly enhanced both AICC and ADCC and also reduced the IL-2 requirement of these cells.^ The stoichiometry of Fc receptor (FcR) expression on the cytotoxic effectors does not correlate with ADCC lytic activity. Although FcRs are necessary to mediate ADCC, other factors, appear to regulate the magnitude of cytolytic activity. In order to investigate these putative factors, the kinetics of ADCC development was studied under various conditions (in IL-2 (10u/ml) and 100u/ml), in IL-2(10u/ml) + TNF$\alpha$ (500u/ml) and in TNF-$\alpha$ (500u/ml) alone). Addition of exogenous TNF-$\alpha$ into the four hour cytotoxicity assay did not increase ADCC, nor did anti-TNF antibodies result in inhibition. On the other hand, addition of anti-TNF antibodies to PBL and IL-2 for 24 hours, resulted in a marked inhibition of the ADCC, suggesting that endogenous TNF-$\alpha$ is obligatory for the maturation and differentiation of ADCC effectors. ^

Circumvention of cellular defense responses to DNA -directed nucleoside analogues by the protein kinase inhibitor, UCN-01

DNA-directed nucleoside analogues, such as ara-C, fludarabine, and gemcitabine, are antimetabolites effective in the treatment of a variety of cancers. However, resistance to nucleoside analogue-based chemotherapy in treatments is still a major problem in therapy. Therefore, it is essential to develop rationales for optimizing the use of nucleoside analogues in combination with other anticancer drugs or modalities such as radiation. The present study focuses on establishing mechanism-based combination strategy to overcome resistance to nucleoside analogues. ^ I hypothesized that the cytostatic concentrations of nucleoside analogues may cause S-phase arrest by activating an S-phase checkpoint that consists of a series of kinases. This may allow cells to repair damaged DNA over time and spare cytotoxicity. Thus, the ability of cells to enact an S-phase arrest in response to incorporation of potentially lethal amounts of nucleoside analogue may serve as a mechanism of resistance to S-phase-specific agents. As a corollary, the addition of a kinase inhibitor, such as UCN-01, may dysregulate the checkpoint response and abrogate the survival of S-phase-arrested cells by suppression of the survival signaling pathways. Using gemcitabine as a model of S-phase-specific nucleoside analogues in human acute myelogenous leukemia ML-1 cells, I demonstrated that cells arrested in S-phase in response to cytostatic conditions. Proliferation continued after washing the cells into drug-free medium, suggesting S-phase arrest served as a resistance mechanism of cancer cells to spare cytotoxicity of nucleoside analogues. However, nontoxic concentrations of UCN-01 rapidly killed S-phase-arrested cells by apoptosis. Furthermore, the molecular mechanism for UCN-01-induced apoptosis in S-phase-arrested cells was through inhibition of survival pathways associated with these cells. In this regard, suppression of the PI 3-kinase-Akt-Bad survival pathway as well as the NF-κB signaling pathway were associated with induction of apoptosis in S-phase-arrested cells by UCN-01, whereas the Ras-Raf-MEK-ERK pathway appeared not involved. This study has provided the rationales and strategies for optimizing the design of effective combination therapies to overcome resistance to nucleoside analogues. In fact, a clinical trial of the combination of ara-C with UCN-01 to treat relapsed or refractory AML patients has been initiated at U.T.M.D. Anderson Cancer Center. ^

Binding of the HIV-1 gp120 -V3 to cellular glycosphingolipids is necessary for CXCR4 recruitment and infection

Infection by human immunodeficiency virus type 1 (HIV-1) is a multi-step process, and detailed analyses of the various events critical for productive infection are necessary to clearly understanding the infection process and identifying novel targets for therapeutic interventions. Evidence from this study reveals binding of the viral envelope protein to host cell glycosphingolipids (GSLs) as a novel event necessary for the orderly progression of the host cell-entry and productive infection by HIV-1. Data obtained from co-immunoprecipitation analyses and confocal microscopy showed that the ability of viral envelope to interact with the co-receptor CXCR4 and productive infection of HIV-1 were inhibited in cells rendered GSL-deficient, while both these activities were restored after reconstitution of the cells with specific GSLs like GM3. Furthermore, evidence was obtained using peptide-inhibitors of HIV-1 infection to show that binding of a specific region within the V3-loop of the envelope protein gp120 to the host cell GSLs is the trigger necessary for the CD4-bound gp120 to recruit the CXCR4 co-receptor. Infection-inhibitory activity of the V3 peptides was compromised in GSL-deficient cells, but could be restored by reconstitution of GSLs. Based on these findings, a revised model for HIV-1 infection is proposed that accounts for the established interactions between the viral envelope and host cell receptors while enumerating the importance of the new findings that fill the gap in the current knowledge of the sequential events for the HIV-1 entry. According to this model, post-CD4 binding of the HIV-1 envelope surface protein gp120 to host cell GSLs, mediated by the gp120-V3 region, enables formation of the gp120-CD4-GSL-CXCR4 immune-complex and productive infection. The identification of cellular GSLs as an additional class of co-factors necessary for HIV-1 infection is important for enhancing the basic knowledge of the HIV-1 entry that can be exploited for developing novel antiviral therapeutic strategies. ^