A BAYESIAN APPROACH TO DOSE-RESPONSE ASSESSMENT AND DRUG-DRUG INTERACTION ANALYSIS: APPLICATION TO IN VITRO STUDIES

The considerable search for synergistic agents in cancer research is motivated by the therapeutic benefits achieved by combining anti-cancer agents. Synergistic agents make it possible to reduce dosage while maintaining or enhancing a desired effect. Other favorable outcomes of synergistic agents include reduction in toxicity and minimizing or delaying drug resistance. Dose-response assessment and drug-drug interaction analysis play an important part in the drug discovery process, however analysis are often poorly done. This dissertation is an effort to notably improve dose-response assessment and drug-drug interaction analysis. The most commonly used method in published analysis is the Median-Effect Principle/Combination Index method (Chou and Talalay, 1984). The Median-Effect Principle/Combination Index method leads to inefficiency by ignoring important sources of variation inherent in dose-response data and discarding data points that do not fit the Median-Effect Principle. Previous work has shown that the conventional method yields a high rate of false positives (Boik, Boik, Newman, 2008; Hennessey, Rosner, Bast, Chen, 2010) and, in some cases, low power to detect synergy. There is a great need for improving the current methodology. We developed a Bayesian framework for dose-response modeling and drug-drug interaction analysis. First, we developed a hierarchical meta-regression dose-response model that accounts for various sources of variation and uncertainty and allows one to incorporate knowledge from prior studies into the current analysis, thus offering a more efficient and reliable inference. Second, in the case that parametric dose-response models do not fit the data, we developed a practical and flexible nonparametric regression method for meta-analysis of independently repeated dose-response experiments. Third, and lastly, we developed a method, based on Loewe additivity that allows one to quantitatively assess interaction between two agents combined at a fixed dose ratio. The proposed method makes a comprehensive and honest account of uncertainty within drug interaction assessment. Extensive simulation studies show that the novel methodology improves the screening process of effective/synergistic agents and reduces the incidence of type I error. We consider an ovarian cancer cell line study that investigates the combined effect of DNA methylation inhibitors and histone deacetylation inhibitors in human ovarian cancer cell lines. The hypothesis is that the combination of DNA methylation inhibitors and histone deacetylation inhibitors will enhance antiproliferative activity in human ovarian cancer cell lines compared to treatment with each inhibitor alone. By applying the proposed Bayesian methodology, in vitro synergy was declared for DNA methylation inhibitor, 5-AZA-2'-deoxycytidine combined with one histone deacetylation inhibitor, suberoylanilide hydroxamic acid or trichostatin A in the cell lines HEY and SKOV3. This suggests potential new epigenetic therapies in cell growth inhibition of ovarian cancer cells.

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

Analysis of a human placental lactogen gene promoter

Human placental lactogen (hPL) and human growth hormone (hGH) comprise a multigene family that share $>$90% nucleic acid sequence homology including 500 bp of 5$\sp\prime$ flanking sequence. Despite these similarities, hGH is produced in the anterior pituitary while hPL is expressed in the placenta. For most genes studied to date, regulation of expression occurs by alterations at the level of transcriptional initiation. Nuclear proteins bind specific DNA sequences in the promoter to regulate gene expression. In this study, the hPL$\sb3$ promoter was analyzed for DNA sequences that contribute to its expression. The interaction between the hPL$\sb3$ promoter and nuclear proteins was examined using nuclear extracts from placental and non-placental cells.^ To identify regulatory elements in the promoter of the hPL$\sb3$ gene, 5$\sp\prime$ deletion mutants were constructed by cleaving 1200 bp of upstream sequence with various restriction enzymes. These DNA fragments were ligated 5$\sp\prime$ to a promoterless bacterial gene chloramphenicol acetyltransferase (CAT) and transfected into JEG-3 cells, a human placental choriocarcinoma cell line. The level of CAT activity reflects the ability of the promoter mutants to activate transcription. Deletion of the sequence between $-$142 bp and $-$129 bp, relative to the start of transcription, resulted in an 8-fold decrease in CAT activity. Nuclear proteins from JEG-3, HeLa, and HepG2 (human liver cells), formed specific binding complexes with this region of the hPL$\sb3$ promoter, as shown by gel mobility shift assay. The $-$142 bp to $-$129 bp region contains a sequence similar to that of a variant binding site for the transcription factor Sp1. Sp1-like proteins were identified by DNA binding assay, in the nuclear extracts of the three cell lines. A series of G nucleotides in the hPL$\sb3$ promoter regulatory region were identified by methylation interference assay to interact with the DNA-binding proteins and the pattern obtained is similar to that for other Sp1 binding sites that have been studied. This suggests that hPL$\sb3$ may be transcriptionally regulated by Sp1 or a Sp1-like transacting factor. ^

DNA methylation inhibits p53 mediated survivin repression in cancer cells

Survivin (BIRC5) is a member of the Inhibitor of Apoptosis (IAP) gene family and functions as a chromosomal passenger protein as well as a mediator of cell survival. Survivin is widely expressed during embryonic development then becomes transcriptionally silent in most highly differentiated adult tissues. It is also overexpressed in virtually every type of tumor. The survivin promoter contains a canonical CpG island that has been described as epigenetically regulated by DNA methylation. We observed that survivin is overexpressed in high grade, poorly differentiated endometrial tumors, and we hypothesized that DNA hypomethylation could explain this expression pattern. Surprisingly, methylation specific PCR and bisulfite pyrosequencing analysis showed that survivin was hypermethylated in endometrial tumors and that this hypermethylation correlated with increased survivin expression. We proposed that methylation could activate survivin expression by inhibit the binding of a transcriptional repressor. ^ The tumor suppressor protein p53 is a well documented transcriptional repressor of survivin and examination of the survivin promoter showed that the p53 binding site contains 3 CpG sites which often become methylated in endometrial tumors. To determine if methylation regulates survivin expression, we treated HCT116 cells with decitabine, a demethylation agent, and observed that survivin transcript and protein levels were significantly repressed following demethylation in a p53 dependent manner. Subsequent binding studies confirmed that DNA methylation inhibited the binding of p53 protein to its binding site in the survivin promoter. ^ We are the first to report this novel mechanism of epigenetic regulation of survivin. We also conducted microarray analysis which showed that many other cancer relevant genes may also be regulated in this manner. While demethylation agents are traditionally thought to inhibit cancer cell growth by reactivating tumor suppressors, our results indicate that an additional important mechanism is to decrease the expression of oncogenes. ^

Analysis of p53 function and regulation

p53 is a tumor suppressor gene that is the most frequent target inactivated in cancers. Overexpression of wild-type p53 in rat embryo fibroblasts suppresses foci formation by other cooperating oncogenes. Introduction of wild-type p53 into cells that lack p53 arrests them at the G1/S boundary and reverses the transformed phenotype of some cells. The function of p53 in normal cells is illustrated by the ability of p53 to arrest cells at G1 phase of the cell cycle upon exposure to DNA-damaging agents including UV-irradiation and biosynthesis inhibitors.^ Since the amino acid sequence of p53 suggested that it may function as a transcription factor, we used GAL4 fusion assays to test that possibility. We found that wild-type p53 could specifically activate transcription when anchored by the GAL4 DNA binding domain. Mutant p53s, which have lost the ability to suppress foci formation by other oncogenes, were not able to activate transcription in this assay. Thus, we established a direct correlation between the tumor suppression and transactivation functions of p53.^ Having learned that p53 was a transcriptional activator, we next sought targets of p53 activation. Because many transcription factors regulate their own expression, we tested whether p53 had this autoregulatory property. Transient expression of wild-type p53 in cells increased the levels of endogenous p53 mRNA. Cotransfection of p53 together with a reporter bearing the p53 promoter confirmed that wild-type p53 specifically activates its own promoter. Deletion analysis from both the 5$\sp\prime$ and 3$\sp\prime$ ends of the promoter minimized the region responsible for p53 autoregulation to 45 bp. Methylation interference identified nucleotides involved in protein-DNA interaction. Mutations within this protected site specifically eliminated the response of the promoter to p53. In addition, multiple copies of this element confer responsiveness to wild-type p53 expression. Thus, we identified a F53 responsive element within the p53 promoter.^ The presence of a consensus NF-$\kappa$B site in the p53 promoter suggested that NF-KB may regulate p53 expression. Gel-shift experiments showed that both the p50 homodimer and the p50/p65 heterodimer bind to the p53 promoter. In addition, the p65 subunit of NF-$\kappa$B activates the p53 promoter in transient transfection experiments. TNF $\alpha$, a natural NF-$\kappa$B inducer, also activates the p53 promoter. Both p65 activation and TNF $\alpha$ induction require an intact NF-$\kappa$B site in the p53 promoter. Since NF-$\kappa$B activation occurs as a response to stress and p53 arrests cells in G1/S, where DNA repair occurs, activation of p53 by NF-$\kappa$B could be a mechanism by which cells recover from stress.^ In conclusion, we provided the first data that wild-type p53 functions as a transcriptional activator, whereas mutant p53 cannot. The correlation between growth suppression and transcriptional activation by p53 implies a pathway of tumor suppression. We have analyzed upstream components of the pathway by the identification of both p53 and NF-$\kappa$B as regulators of the p53 promoter. ^

Molecular and genetic analysis of Myxococcus xanthus early *development

To identify more mutations that can affect the early development of Myxococcus xanthus, the synthetic transposon TnT41 was designed and constructed. By virtue of its special features, it can greatly facilitate the processes of mutation screening/selection, mapping, cloning and DNA sequencing. In addition, it allows for the systematic discovery of genes in regulatory hierarchies using their target promoters. In this study, the minimal regulatory region of the early developmentally regulated gene 4521 was used as a reporter in the TnT41 mutagenesis. Both positive (P) mutations and negative (N) mutations were isolated based on their effects on 4521 expression.^ Four of these mutations, i.e. N1, N2, P52 and P54 were analyzed in detail. Mutations N1 and N2 are insertion mutations in a gene designated sasB. The sasB gene is also identified in this study by genetic and molecular analysis of five UV-generated 4521 suppressor mutations. The sasB gene encodes a protein without meaningful homology in the databases. The sasB gene negatively regulates 4521 expression possibly through the SasS-SasR two component system. A wild-type sasB gene is required for normal M. xanthus fruiting body formation and sporulation.^ Cloning and sequencing analysis of the P52 mutation led to the identification of an operon that encodes the M. xanthus high-affinity branched-chain amino acid transporter system. This liv operon consists of five genes designated livK, livH, livM, livC, and livF, respectively. The Liv proteins are highly similar to their counterparts from other bacteria in both amino acid sequences, functional motifs and predicted secondary structures. This system is required for development since liv null mutations cause abnormality in fruiting body formation and a 100-fold decrease in sporulation efficiency.^ Mutation P54 is a TnT41 insertion in the sscM gene of the ssc chemotaxis system, which has been independently identified by Dr. Shi's lab. The sscM gene encodes a MCP (methyl-accepting chemotaxis protein) homologue. The SscM protein is predicted to contain two transmembrane domains, a signaling domain and at least one putative methylation site. Null mutations of this gene abolish the aggregation of starving cells at a very early stage, though the sporulation levels of the mutant can reach 10% that of wild-type cells. ^