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.

XENOESTROGEN-SPECIFIC MECHANISMS OF DEVELOPMENTAL REPROGRAMMING CORRELATE WITH GENE EXPRESSION AND TUMOR DEVELOPMENT

Environmental exposures during sensitive windows of development can reprogram normal physiological responses and alter disease susceptibility later in life in a process known as developmental reprogramming. We have shown that neonatal exposure to the xenoestrogen diethylstilbestrol (DES) can developmentally reprogram the reproductive tract in genetically susceptible Eker rats giving rise to complete penetrance of uterine leiomyoma. Based on this, we hypothesized that xenoestrogens, including genistein (GEN) and bisphenol A (BPA), reprogram estrogen-responsive gene expression in the myometrium and promote the development of uterine leiomyoma. We proposed the mechanism that is responsible for the developmental reprogramming of gene expression was through estrogen (E2)/ xenoestrogen inducedrapid ER signaling, which modifies the histone methyltransferase Enhancer of Zeste homolog 2 (EZH2) via activation of the PI3K/AKT pathway. We further hypothesized that there is a xenostrogen-specific effect on this pathway altering patterns of histone modification, DNA methylation and gene expression. In addition to our novel finding that E2/DES-induced phosphorylation of EZH2 by AKT reduces the levels of H3K27me3 in vitro and in vivo, this work demonstrates in vivo that a brief neonatal exposure to GEN, in contrast to BPA, activates the PI3K/AKT pathway to regulate EZH2 and decreases H3K27me3 levels in the neonatal uterus. Given that H3K27me3 is a repressive mark that has been shown to result in DNA methylation and gene silencing we investigated the methylation of developmentally reprogrammed genes. In support of this evidence, we show that neonatal DES exposure in comparison to VEH, leads to hypomethylation of the promoter of a developmentally reprogrammed gene, Gria2, that become hyper-responsive to estrogen in the adult myometrium indicating vi that DES exposure alter gene expression via chromatin remodeling and loss of DNA methylation. In the adult uterus, GEN and BPA exposure developmentally reprogrammed expression of estrogen-responsive genes in a manner opposite of one another, correlating with our previous data. Furthermore, the ability of GEN and BPA to developmental reprogram gene expression correlated with tumor incidence and multiplicity. These data show that xenoestrogens have unique effects on the activation of non-genomic signaling in the developing uterus that promotes epigenetic and genetic alterations, which are predictive of developmental reprogramming and correlate with their ability to modulate hormone-dependent tumor development.

Rubisco in complex with Rubisco large subunit methyltransferase.

SET domain protein lysine methyltransferases (PKMT) are a structurally unique class of enzymes that catalyze the specific methylation of lysine residues in a number of different substrates. Especially histone-specific SET domain PKMTs have received widespread attention because of their roles in the regulation of epigenetic gene expression and the development of some cancers. Rubisco large subunit methyltransferase (RLSMT) is a chloroplast-localized SET domain PKMT responsible for the formation of trimethyl-lysine-14 in the large subunit of Rubisco, an essential photosynthetic enzyme. Here, we have used cryoelectron microscopy to produce an 11-A density map of the Rubisco-RLSMT complex. The atomic model of the complex, obtained by fitting crystal structures of Rubisco and RLSMT into the density map, shows that the extensive contact regions between the 2 proteins are mainly mediated by hydrophobic residues and leucine-rich repeats. It further provides insights into potential conformational changes that may occur during substrate binding and catalysis. This study presents the first structural analysis of a SET domain PKMT in complex with its intact polypeptide substrate.

CREB trans-activation of disruptor of telomeric silencing-1 mediates forskolin inhibition of CTGF transcription in mesangial cells.

Connective tissue growth factor (CTGF) participates in diverse fibrotic processes including glomerulosclerosis. The adenylyl cyclase agonist forskolin inhibits CTGF expression in mesangial cells by unclear mechanisms. We recently reported that the histone H3K79 methyltransferase disruptor of telomeric silencing-1 (Dot1) suppresses CTGF gene expression in collecting duct cells (J Clin Invest 117: 773-783, 2007) and HEK 293 cells (J Biol Chem In press). In the present study, we characterized the involvement of Dot1 in mediating the inhibitory effect of forskolin on CTGF transcription in mouse mesangial cells. Overexpression of Dot1 or treatment with forskolin dramatically suppressed basal CTGF mRNA levels and CTGF promoter-luciferase activity, while hypermethylating H3K79 in chromatin associated with the CTGF promoter. siRNA knockdown of Dot1 abrogated the inhibitory effect of forskolin on CTGF mRNA expression. Analysis of the Dot1 promoter sequence identified a CREB response element (CRE) at -384/-380. Overexpression of CREB enhanced forskolin-stimulated Dot1 promoter activity. A constitutively active CREB mutant (CREB-VP16) strongly induced Dot1 promoter-luciferase activity, whereas overexpression of CREBdLZ-VP16, which lacks the CREB DNA-binding domain, abolished this activation. Mutation of the -384/-380 CRE resulted in 70% lower levels of Dot1 promoter activity. ChIP assays confirmed CREB binding to the Dot1 promoter in chromatin. We conclude that forskolin stimulates CREB-mediated trans-activation of the Dot1 gene, which leads to hypermethylation of histone H3K79 at the CTGF promoter, and inhibition of CTGF transcription. These data are the first to describe regulation of the Dot1 gene, and disclose a complex network of genetic and epigenetic controls on CTGF transcription.

Aldosterone-induced Sgk1 relieves Dot1a-Af9-mediated transcriptional repression of epithelial Na+ channel alpha.

Aldosterone plays a major role in the regulation of salt balance and the pathophysiology of cardiovascular and renal diseases. Many aldosterone-regulated genes--including that encoding the epithelial Na+ channel (ENaC), a key arbiter of Na+ transport in the kidney and other epithelia--have been identified, but the mechanisms by which the hormone modifies chromatin structure and thus transcription remain unknown. We previously described the basal repression of ENaCalpha by a complex containing the histone H3 Lys79 methyltransferase disruptor of telomeric silencing alternative splice variant a (Dot1a) and the putative transcription factor ALL1-fused gene from chromosome 9 (Af9) as well as the release of this repression by aldosterone treatment. Here we provide evidence from renal collecting duct cells and serum- and glucocorticoid-induced kinase-1 (Sgk1) WT and knockout mice that Sgk1 phosphorylated Af9, thereby impairing the Dot1a-Af9 interaction and leading to targeted histone H3 Lys79 hypomethylation at the ENaCalpha promoter and derepression of ENaCalpha transcription. Thus, Af9 is a physiologic target of Sgk1, and Sgk1 negatively regulates the Dot1a-Af9 repressor complex that controls transcription of ENaCalpha and likely other aldosterone-induced genes.

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. ^

Characterization of the human tissue transglutaminase gene promoter

Transglutaminases are a family of calcium-dependent enzymes, that catalyze the covalent cross-linking of proteins by forming $\varepsilon(\gamma$-glutamyl)lysine isopeptide bonds. In order to investigate the molecular mechanisms regulating the expression of the tissue transglutaminase gene and to determine its biological functions, the goal of this research has been to clone and characterize the human tissue transglutaminase promoter. Thirteen clones of the tissue transglutaminase gene were obtained from the screening of a human placental genomic DNA library. A 1.74 Kb fragment derived from DNA located immediately upstream of the translation start site was subcloned and sequenced. Sequence analysis of this DNA fragment revealed that it contains a TATA box (TATAA), a CAAT box (GGACAAT), and a series of potential transcription factor binding sites and hormone response elements. Four regions of significant homology, a GC-rich region, a TG-rich region, an AG-rich region, and HR1, were identified by aligning 1.8 Kb of DNA flanking the human, mouse, and guinea pig tissue transglutaminase genes.^ To measure promoter activity, we subcloned the 1.74 Kb fragment of the tissue transglutaminase gene into a luciferase reporter vector to generate transglutaminase promoter/luciferase reporter constructs. Transfection experiments showed that this DNA segment includes a functional promoter with high constitutive activity. Deletion analysis revealed that the SP1 sites or corresponding sequences contribute to this activity. We investigated the role of DNA methylation in regulating the activity of the promoter and found that in vitro methylation of tissue transglutaminase promoter/luciferase reporter constructs suppressed their basal activity. Methylation of the promoter is inversely correlated with the expression of the tissue transglutaminase gene in vivo. These results suggest that DNA methylation may be one of the mechanisms regulating the expression of the gene. The tumor suppressor gene product p53 was also shown to inhibit the activity of the promoter, suggesting that induction of the tissue transglutaminase gene is not involved in the p53-dependent programmed cell death pathway. Although retinoids regulate the expression of the tissue transglutaminase gene in vivo, retinoid-inducible activity can not be identified in 3.7 Kb of DNA 5$\sp\prime$ to the tissue transglutaminase gene.^ The structure of the 5$\sp\prime$ end of the tissue transglutaminase gene was mapped. Alignment analysis of the human tissue transglutaminase gene with other human transglutaminases showed that tissue transglutaminase is the simplest member of transglutaminase superfamily. Transglutaminase genes show a conserved core of exons and introns but diverse N-terminuses and promoters. These observations suggest that key regulatory sequences and promoter elements have been appended upstream of the core transglutaminase gene to generate the diversity of regulated expression and regulated activity characteristic of the transglutaminase gene family. ^

Human TNF receptor p75/80 gene structure and promoter characterization

Tumor necrosis factor receptor p75/80 ((TNF-R p75/80) is a 75 kDa type 1 transmembrane protein expressed predominately on cells of hematopoietic lineage. TNF-R p75/80 belongs to the TNF receptor superfamily characterized by cysteine-rich extracellular regions composed of three to six disulfide-linked domains. In the present report, we have characterized, for the first time, the complete gene structure for human TNF-R p75/80 which spans approximately 43 kbp. The gene consists of 10 exons (ranging from 34 bp to 2.5 kbp) and 9 introns (343 bp to 19 kbp). Consensus elements for transcription factors involved in T cell development and activation were noted in the 5$\sp\prime$ flanking region including TCF-1, Ikaros, AP-1, CK-2, IL-6RE, ISRE, GAS, NF-$\kappa$B and SP1, as well as an unusually high GC content and CpG frequency that appears characteristic of some TNF-R family members. The unusual (GATA)$\sb{\rm n}$ and (GAA)(GGA) repeats found within intron 1 may prove useful for further genome analysis within the 1p36 chromosomal locus. The human TNF-R p75/80 gene structure will permit further assessment of its involvement in normal hematopoietic cell development and function, autoimmune disease, and non-random translocations in hematopoietic malignancies. The region 1.8 kb 5$\sp\prime$ of the ATG was able to drive luciferase expression when transfected into cell lines expressing TNF-R p75/80. Further characterization of the 5$\sp\prime$-regulatory region will aid in determining factors and signal transduction pathways involved in regulating TNF-R p75/80 expression. ^

Role of transducer proteins in photoaxis signaling of Halobacterium salinarum

In Halobacterium salinarum phototaxis is mediated by the visual pigment-like photoreceptors sensory rhodopsin I (SRI) and II (SRII). SRI is a receptor for attractant orange and repellent UV-blue light, and SRII is a receptor for repellent blue-green light, and transmit signals through the membrane-bound transducer proteins HtrI and HtrII, respectively. ^ The primary sequences of HtrI and HtrII predict 2 transmembrane helices (TM1 and TM2) followed by a hydrophilic cytoplasmic domain. HtrII shows an additional large periplasmic domain for chemotactic ligand binding. The cytoplasmic regions are homologous to the adaptation and signaling domains of eubacterial chemotaxis receptors and, like their eubacterial homologs, modulate the transfer of phosphate groups from the histidine protein kinase CheA to the response regulator CheY that in turn controls flagellar motor rotation and the cell's swimming behavior. HtrII and Htrl are dimeric proteins which were predicted to contain carboxylmethylation sites in a 4-helix bundle in their cytoplasmic regions, like eubacterial chemotaxis receptors. ^ The phototaxis transducers of H. salinarum have provided a model for studying receptor/tranducer interaction, adaptation in sensory systems, and the role of membrane molecular complexes in signal transduction. ^ Interaction between the transducer HtrI and the photoreceptor SRI was explored by creating six deletion constructs of HtrI, with progressively shorter cytoplasmic domains. This study confirmed a putative chaperone-like function of HtrI, facilitating membrane insertion or stability of the SRI protein, a phenomenon previously observed in the laboratory, and identified the smallest HtrI fragment containing interaction sites for both the chaperone-like function and SRI photocycle control. The active fragment consisted of the N-terminal 147 residues of the 536-residue HtrI protein, a portion of the molecule predicted to contain the two transmembrane helices and the first ∼20% of the cytoplasmic portion of the protein. ^ Phototaxis and chemotaxis sensory systems adapt to stimuli, thereby signaling only in response to changes in environmental conditions. Observations made in our and in other laboratories and homologies between the halobacterial transducers with the chemoreceptors of enteric bacteria anticipated a role for methylation in adaptation to chemo- and photostimuli. By site directed mutagenesis we identified the methylation sites to be the glutamate pairs E265–E266 in HtrI and E513–E514 in HtrII. Cells containing the unmethylatable transducers are still able to perform phototaxis and adapt to light stimuli. By pulse-chase analysis we found that methanol production from carboxylmethyl group hydrolysis occurs upon specific photo stimulation of unmethylatable HtrI and HtrII and is due to turnover of methyl groups on other transducers. We demonstrated that the turnover in wild-type H. salinarum cells that follows a positive stimulus is CheY-dependent. The CheY-feedback pathway does not require the stimulated transducer to be methylatable and operates globally on other transducers present in the cell. ^ Assembly of signaling molecules into architecturally defined complexes is considered essential in transmission of the signals. The spectroscopic characteristics of SRI were exploited to study the stoichiometric composition in the phototaxis complex SRI-HtrI. A molar ratio of 2.1 HtrI: 1 SRI was obtained, suggesting that only 1 SRI binding site is occupied on the HtrI homodimer. We used gold-immunoelectron microscopy and light fluorescence microscopy to investigate the structural organization and the distribution of other halobacterial transducers. We detected clusters of transducers, usually near the cell's poles, providing a ultrastructural basis for the global effects and intertransducer communication we observe. ^