A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

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.

THE ROLE OF MAP KINASES ON THE FUNCTIONAL HETEROGENEITY OF HUMAN CD8+ T CELL MATURATION SUBSETS

While prior studies have focused on naïve (CD45RA+CD27+) and early stage memory (CD45RA-CD27+) CD8+ T cells, late memory CD8+ T cells (CD45RA+CD27) have received less interest because this subset of T cells is generally recognized as effectors, which produce IFNγ (but no IL-2) and perforin. However, multiple studies suggest that late memory CD8+ T cells may provide inadequate protection in infectious diseases and cancer models. To better understand the unique function of late memory CD8+ T cells, I optimized multi-color flow cytometry techniques to assess the cytokine production of each human CD8+ T cell maturation subset. I demonstrated that late memory CD8+ T cells are the predominant producer of CC chemokines (e.g. MIP-1β), but rarely produce IL-2; therefore they do not co-produce IL-2/IFNγ (polyfunctionality), which has been shown to be critical for protective immunity against chronic viral infection. These data suggest that late memory CD8+ T cells are not just cytotoxic effectors, but may have unique functional properties. Determining the molecular signature of each CD8+ T cell maturation subset will help characterize the role of late memory CD8+ T cells. Prior studies suggest that ERK1 and ERK2 play a role in cytokine production including IL-2 in T cells. Therefore, I tested whether differential expression of ERK1 and ERK2 in CD8+ T cell maturation subsets contributes to their functional signature by a novel flow cytometry technique. I found that the expression of total ERK1, but not ERK2, is significantly diminished in late memory CD8+ T cells and that ERK1 expression is strongly associated with IL-2 production and CD28 expression. I also found that IL-2 production is increased in late memory CD8+ T cells by over-expressing ERK1. Collectively, these data suggest that ERK1 is required for IL-2 production in human CD8+ T cells. In summary, this dissertation demonstrated that ERK1 is down-regulated in human late memory CD8+ T cells, leading to decreased production of IL-2. The data in this dissertation also suggested that the functional heterogeneity in human CD8+ T cell maturation subsets results from their differential ERK1 expression.

Role of acetyl-phosphate in activation of the Rrp2-RpoN-RpoS pathway in Borrelia burgdorferi.

Borrelia burgdorferi, the Lyme disease spirochete, dramatically alters its transcriptome and proteome as it cycles between the arthropod vector and mammalian host. During this enzootic cycle, a novel regulatory network, the Rrp2-RpoN-RpoS pathway (also known as the σ(54)-σ(S) sigma factor cascade), plays a central role in modulating the differential expression of more than 10% of all B. burgdorferi genes, including the major virulence genes ospA and ospC. However, the mechanism(s) by which the upstream activator and response regulator Rrp2 is activated remains unclear. Here, we show that none of the histidine kinases present in the B. burgdorferi genome are required for the activation of Rrp2. Instead, we present biochemical and genetic evidence that supports the hypothesis that activation of the Rrp2-RpoN-RpoS pathway occurs via the small, high-energy, phosphoryl-donor acetyl phosphate (acetyl∼P), the intermediate of the Ack-Pta (acetate kinase-phosphate acetyltransferase) pathway that converts acetate to acetyl-CoA. Supplementation of the growth medium with acetate induced activation of the Rrp2-RpoN-RpoS pathway in a dose-dependent manner. Conversely, the overexpression of Pta virtually abolished acetate-induced activation of this pathway, suggesting that acetate works through acetyl∼P. Overexpression of Pta also greatly inhibited temperature and cell density-induced activation of RpoS and OspC, suggesting that these environmental cues affect the Rrp2-RpoN-RpoS pathway by influencing acetyl∼P. Finally, overexpression of Pta partially reduced infectivity of B. burgdorferi in mice. Taken together, these findings suggest that acetyl∼P is one of the key activating molecule for the activation of the Rrp2-RpoN-RpoS pathway and support the emerging concept that acetyl∼P can serve as a global signal in bacterial pathogenesis.

Role of the N- and C-lobes of calmodulin in the activation of Ca(2+)/calmodulin-dependent protein kinase II.

Understanding the principles of calmodulin (CaM) activation of target enzymes will help delineate how this seemingly simple molecule can play such a complex role in transducing Ca (2+)-signals to a variety of downstream pathways. In the work reported here, we use biochemical and biophysical tools and a panel of CaM constructs to examine the lobe specific interactions between CaM and CaMKII necessary for the activation and autophosphorylation of the enzyme. Interestingly, the N-terminal lobe of CaM by itself was able to partially activate and allow autophosphorylation of CaMKII while the C-terminal lobe was inactive. When used together, CaMN and CaMC produced maximal CaMKII activation and autophosphorylation. Moreover, CaMNN and CaMCC (chimeras of the two N- or C-terminal lobes) both activated the kinase but with greater K act than for wtCaM. Isothermal titration calorimetry experiments showed the same rank order of affinities of wtCaM > CaMNN > CaMCC as those determined in the activity assay and that the CaM to CaMKII subunit binding ratio was 1:1. Together, our results lead to a proposed sequential mechanism to describe the activation pathway of CaMKII led by binding of the N-lobe followed by the C-lobe. This mechanism contrasts the typical sequential binding mode of CaM with other CaM-dependent enzymes, where the C-lobe of CaM binds first. The consequence of such lobe specific binding mechanisms is discussed in relation to the differential rates of Ca (2+)-binding to each lobe of CaM during intracellular Ca (2+) oscillations.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

Role of the transcription factor activator protein-2alpha in prostate carcinogenesis

Activator protein 2α (AP-2) is a transcription factor known to play a crucial role in the progression of malignant melanoma, colorectal carcinoma, and breast cancer. Several AP-2 target genes are known to be deregulated in prostate cancer, therefore, we hypothesize that loss AP-2 expression plays a causal role in prostate carcinogenesis. Immunofluorescent staining for AP-2 of 30 radical prostatectomy specimens demonstrated that while AP-2 was highly expressed in normal prostate epithelium, its expression was lost in most cases of high grade prostatic intraepithelial neoplasia (PIN), and all cases of prostate cancer studied. Additional analyses demonstrated that AP-2 was associated with normal luminal differentiation and it was not expressed in the basal cell layer. In cell lines, AP-2 was strongly expressed in immortalized normal prostate epithelial cells, whereas low expression was observed in the LNCaP, LNCaP-LN3, and PC3M-LN4 prostate cancer cell lines. Transfection of the highly tumorigenic and metastatic cell line PC3M-LN4 with the AP-2 gene significantly decreased tumor growth in the prostate of nude mice (p = 0.032) and inhibited metastases to the lymph nodes. Moreover, transfection of the low tumorigenic, low metastatic cell line LNCaP-LN3 with full length AP-2; resulted in complete inhibition of tumor incidence in the AP-2 transfectants (0/19) vs. neo control (10/16). A potential mechanism for this loss of tumorigenicity was the modulation of gene expression in prostate cancer cells that mimicked the normal phenotype. Analysis of differential expression between neo control- and AP-2-transfected cells in vitro and in tumors demonstrated low VEGF expression in AP-2 transfectants. We further demonstrated that AP-2 acted as a transcriptional repressor of the VEGF promoter by binding to a GC-rich region located between −88 and −66. This region contains an AP-2 consensus element overlapping two Sp1 consensus elements. We found that Sp3 and AP-2 bound to this region in a mutually exclusive manner to promote activation or repression. Increased VEGF expression has been observed in high grade PIN and in prostate cancer. Here we provide evidence that this early molecular change could be a result of loss of AP-2 expression in the prostatic epithelium. ^

Glioma-specific alternative RNA splicing: A role for polypyrimidine tract binding protein-1

Alternative RNA splicing plays an integral role in cell fate determination and function, especially in the cells of the brain. Errors in RNA processing contribute to diseases such as cancer, where it leads to the production of oncogenic proteins or the loss of tumor suppressors. In silica mining suggests that hundreds of splice isoforms are misexpressed in the glial cell-derived glioma. However, there is little experimental evidence of the prevalence and contribution of these changes and whether they contribute to the formation and progression of this devastating malignancy. To determine the frequency of these aberrant events, global profiling of alternative RNA splice patterns in glioma and nontumor brain was conducted using an exon array. Most splicing changes were less than 5-fold in magnitude and 14 cassette exon events were validated, including 7 previously published events. To determine the possible causes of missplicing, the differential expression levels of splicing factors in these two tissues were also analyzed. Six RNA splicing factors had greater than 2-fold changes in expression. The highest differentially expressed factor was polypyrimidine tract binding protein-1 (PTB). Evaluation by immunohistochemistry determined that this factor was elevated in both early and late stages of glioma. Glial cell-specific PTB expression in the adult brain led me to examine the role of PTB in gliomagenesis. Downregulation of PTB slowed glioma cell proliferation and migration and enhanced cell adhesion to fibronectin and vitronectin. To determine whether PTB was affecting these processes through splicing, genome-wide exon expression levels were correlated with PTB levels. Surprisingly, previously reported PTB target transcripts were insensitive to changes in PTB levels in both patient samples and PTB-depleted glioma cells. Only one validated glioma-specific splice target, RTN4/Nogo, had a significant PTB-mediated splicing change. Downregulation of PTB enhanced inclusion of its alternative exon 3, which encodes an auxiliary domain within a neurite inhibitor protein. Overexpression of this splice isoform in glioma cells slowed proliferation in a manner similar to that observed in PTB knockdown cells. In summary, aberrant expression of splicing factors such as PTB in glioma may elicit changes in splicing patterns that enhance tumorigenesis. ^

Differential Activity of the KRAS Oncogene by Method of Activation: Implications for Signaling and Therapeutic Intervention

Despite having been identified over thirty years ago and definitively established as having a critical role in driving tumor growth and predicting for resistance to therapy, the KRAS oncogene remains a target in cancer for which there is no effective treatment. KRas is activated b y mutations at a few sites, primarily amino acid substitutions at codon 12 which promote a constitutively active state. I have found that different amino acid substitutions at codon 12 can activate different KRas downstream signaling pathways, determine clonogenic growth potential and determine patient response to molecularly targeted therapies. Computer modeling of the KRas structure shows that different amino acids substituted at the codon 12 position influences how KRas interacts with its effecters. In the absence of a direct inhibitor of mutant KRas several agents have recently entered clinical trials alone and in combination directly targeting two of the common downstream effecter pathways of KRas, namely the Mapk pathway and the Akt pathway. These inhibitors were evaluated for efficacy against different KRAS activating mutations. An isogenic panel of colorectal cells with wild type KRas replaced with KRas G12C, G12D, or G12V at the endogenous loci differed in sensitivity to Mek and Akt inhibition. In contrast, screening was performed in a broad panel of lung cell lines alone and no correlation was seen between types of activating KRAS mutation due to concurrent oncogenic lesions. To find a new method to inhibit KRAS driven tumors, siRNA screens were performed in isogenic lines with and without active KRas. The knockdown of CNKSR1 (CNK1) showed selective growth inhibition in cells with an oncogenic KRAS. The deletion of CNK1 reduces expression of mitotic cell cycle proteins and arrests cells with active KRas in the G1 phase of the cell cycle similar to the deletion of an activated KRas regardless of activating substitution. CNK1 has a PH domain responsible for localizing it to membrane lipids making KRas potentially amenable to inhibition with small molecules. The work has identified a series of small molecules capable of binding to this PH domain and inhibiting CNK1 facilitated KRas signaling.

Targeted Deletion Of Functionally Validated Enhancers Defines Their Role in Mouse Limb Development

Transcriptional enhancers are genomic DNA sequences that contain clustered transcription factor (TF) binding sites. When combinations of TFs bind to enhancer sequences they act together with basal transcriptional machinery to regulate the timing, location and quantity of gene transcription. Elucidating the genetic mechanisms responsible for differential gene expression, including the role of enhancers, during embryological and postnatal development is essential to an understanding of evolutionary processes and disease etiology. Numerous methods are in use to identify and characterize enhancers. Several high-throughput methods generate large datasets of enhancer sequences with putative roles in embryonic development. However, few enhancers have been deleted from the genome to determine their roles in the development of specific structures, such as the limb. Manipulation of enhancers at their endogenous loci, such as the deletion of such elements, leads to a better understanding of the regulatory interactions, rules and complexities that contribute to faithful and variant gene transcription – the molecular genetic substrate of evolution and disease. To understand the endogenous roles of two distinct enhancers known to be active in the mouse embryo limb bud we deleted them from the mouse genome. I hypothesized that deletion of these enhancers would lead to aberrant limb development. The enhancers were selected because of their association with p300, a protein associated with active transcription, and because the human enhancer sequences drive distinct lacZ expression patterns in limb buds of embryonic day (E) 11.5 transgenic mice. To confirm that the orthologous mouse enhancers, mouse 280 and 1442 (M280 and M1442, respectively), regulate expression in the developing limb we generated stable transgenic lines, and examined lacZ expression. In M280-lacZ mice, expression was detected in E11.5 fore- and hindlimbs in a region that corresponds to digits II-IV. M1442-lacZ mice exhibited lacZ expression in posterior and anterior margins of the fore- and hindlimbs that overlapped with digits I and V and several wrist bones. We generated mice lacking the M280 and M1442 enhancers by gene targeting. Intercrosses between M280 -/+ and M1442 -/+, respectively, generated M280 and M1442 null mice, which are born at expected Mendelian ratios and manifest no gross limb malformations. Quantitative real-time PCR of mutant E11.5 limb buds indicated that significant changes in transcriptional output of enhancer-proximal genes accompanied the deletion of both M280 and M1442. In neonatal null mice we observed that all limb bones are present in their expected positions, an observation also confirmed by histology of E18.5 distal limbs. Fine-scale measurement of E18.5 digit bone lengths found no differences between mutant and control embryos. Furthermore, when the developmental progression of cartilaginous elements was analyzed in M280 and M1442 embryos from E13.5-E15.5, transient development defects were not detected. These results demonstrate that M280 and M1442 are not required for mouse limb development. Though M280 is not required for embryonic limb development it is required for the development and/or maintenance of body size – adult M280 mice are significantly smaller than control littermates. These studies highlight the importance of experiments that manipulate enhancers in situ to understand their contribution to development.