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.

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

Molecular analysis of a neurovirulent murine leukemia virus, {\it ts\/}1

ts1 is a neurovirulent spontaneous temperature-sensitive mutant of Moloney murine leukemia virus TB (MoMuLV-TB). MoMuLV-TB causes T-cell lymphoma or lymphoid leukemia in mice after a long latency period whereas ts1 causes a progressive hindlimb paralytic disease after a much shorter latency period. In previous studies, it had been shown that the temperature-sensitive defect resided in the $env$ gene. At the restrictive temperature, the envelope precursor polyprotein, gPr80$\sp{env}$, is inefficiently processed intracellularly into a heterodimer consisting of two cleavage products, gp70 and Prp15E. This inefficient processing is correlated with neurovirulence. In this study, the nucleotide sequences of the env genes for both ts1 and MoMuLV-TB were determined, and the encoded amino acid sequences were deduced from the DNA sequences. There were four unique amino acid substitutions in the gPr80$\sp{env}$ of ts1. In order to determine which unique amino acid was responsible for the phenotypic characteristics of ts1, a set of hybrid genomes was constructed by exchanging restriction fragments between ts1 and MoMuLV-TB. NIH 3T3 cells were transfected with the hybrid genomes to obtain infectious hybrid viruses. Assays of the hybrid viruses showed that a Val-25$\to$Ile substitution in gPr80$\sp{env}$ was responsible for the temperature sensitivity, inefficient processing, and neurovirulence of ts1. In further studies, the Ile-25 in gPr80$\sp{env}$ was substituted with Thr, Ala, Leu, Gly, and Glu by site-directed mutagenesis to generate a new set of mutant viruses, i.e., ts1-T, -A, -L, -G, and -E, respectively. The rank order of the mutants for temperature sensitivity was: ts1-E $>$ ts1-G $>$ ts1-L $>$ ts1-A $>$ ts1 $>$ ts1-T. The degree of temperature sensitivity of each of the mutants also correlated with the degree of inefficient processing of gPr80$\sp{env}$. The mutant viruses were assayed for neurovirulence. ts1-T caused whole body tremor, ts1-A caused hindlimb paralysis, ts1-L caused paraparesis, but ts1-G and -E were not neurovirulent. These results show that inefficient processing of gPr80$\sp{env}$ is correlated with neurovirulence, but if processing of gPr80$\sp{env}$ is too inefficient there is no neurovirulence. Furthermore, the disease profile of each of the neurovirulent viruses depends on the degree of inefficient processing of gPr80$\sp{env}$. ^

Constitutively active cAMP receptor mutants dominantly inhibit receptor-mediated signaling pathways in Dictyostelium discoideum

Dictyostelium, a soil amoeba, is able to develop from free-living cells to multicellular fruiting bodies upon starvation using extracellular cAMP to mediate cell-cell communication, chemotaxis and developmental gene expression. The seven transmembrane G protein-coupled cAMP receptor-1 (cAR1) mediated responses, such as the activation of adenylyl cyclase and guanylyl cyclase, are transient, due to the existence of poorly understood adaptation mechanisms. For this dissertation, the powerful genetics of the Dictyostelium system was employed to study the adaptation mechanism of cAR1-mediated cAMP signaling as well as mechanisms intrinsic to cAR1 that regulate its activation. ^ We proposed that constitutively active cAR1 would cause constant adaptation, thus inhibiting downstream pathways that are essential for aggregation and development. Therefore, a screen for dominant negative cAR1 mutants was undertaken to identify constitutively active receptor mutants. Three dominant negative cAR1 mutants were identified. All appear to be constitutively active receptor mutants because they are constitutively phosphorylated and possess high affinity for cAMP. Biochemical studies showed that these mutant receptors prevented the activation of downstream effectors, including adenylyl and guanylyl cyclases. In addition, these cells also were defective in cAMP chemotaxis and cAR1-mediated gene expression. These findings suggest that the mutant receptors block development by constantly activating multiple adaptation pathways. ^ Sequence analysis revealed that these mutations (I104N, L100H) are clustered in a conserved region of the third transmembrane helix (TM3) of cAR1. To investigate the role of this region in receptor activation, one of these residues, I104, was mutated to all the other 19 possible amino acids. We found that all but the most conservative substitutions increase the receptor's affinity about 20- to 70-fold. However, only highly polar substitutions of I104, particularly basic residues, resulted in receptors that are constitutively phosphorylated and dominantly inhibit development, suggesting that highly polar substitutions not only disrupt an interaction constraining the receptor in its low-affinity, inactive state but also promote an additional conformational change that resembles the ligand-bound conformation. Our findings suggest that I104 plays a specific role in constraining the receptor in its inactive state and that substituting it with highly polar residues results in constitutive activation. ^