Synergistic inhibition of tumor growth in a murine mammary adenocarcinoma model by combinational gene therapy using IL-12, pro-IL-18, and IL-1β converting enzyme cDNA

We report here that a cancer gene therapy protocol using a combination of IL-12, pro-IL-18, and IL-1β converting enzyme (ICE) cDNA expression vectors simultaneously delivered via gene gun can significantly augment antitumor effects, evidently by generating increased levels of bioactive IL-18 and consequently IFN-γ. First, we compared the levels of IFN-γ secreted by mouse splenocytes stimulated with tumor cells transfected with various test genes, including IL-12 alone; pro-IL-18 alone; pro-IL-18 and ICE; IL-12 and pro-IL-18; and IL-12, pro-IL-18, and ICE. Among these treatments, the combination of IL-12, pro-IL-18, and ICE cDNA resulted in the highest level of IFN-γ production from splenocytes in vitro, and similar results were obtained when these same treatments were delivered to the skin of a mouse by gene gun and IFN-γ levels were measured at the skin transfection site in vivo. Furthermore, the triple gene combinatorial gene therapy protocol was the most effective among all tested groups at suppressing the growth of TS/A (murine mammary adenocarcinoma) tumors previously implanted intradermally at the skin site receiving DNA transfer by gene gun on days 6, 8, 10, and 12 after tumor implantation. Fifty percent of mice treated with the combined three-gene protocol underwent complete tumor regression. In vivo depletion experiments showed that this antitumor effect was CD8+ T cell-mediated and partially IFN-γ-dependent. These results suggest that a combinatorial gene therapy protocol using a mixture of IL-12, pro-IL-18, and ICE cDNAs can confer potent antitumor activities against established TS/A tumors via cytotoxic CD8+ T cells and IFN-γ-dependent pathways.

Gene2EST: a BLAST2 server for searching expressed sequence tag (EST) databases with eukaryotic gene-sized queries

Expressed sequence tags (ESTs) are randomly sequenced cDNA clones. Currently, nearly 3 million human and 2 million mouse ESTs provide valuable resources that enable researchers to investigate the products of gene expression. The EST databases have proven to be useful tools for detecting homologous genes, for exon mapping, revealing differential splicing, etc. With the increasing availability of large amounts of poorly characterised eukaryotic (notably human) genomic sequence, ESTs have now become a vital tool for gene identification, sometimes yielding the only unambiguous evidence for the existence of a gene expression product. However, BLAST-based Web servers available to the general user have not kept pace with these developments and do not provide appropriate tools for querying EST databases with large highly spliced genes, often spanning 50 000–100 000 bases or more. Here we describe Gene2EST (http://woody.embl-heidelberg.de/gene2est/), a server that brings together a set of tools enabling efficient retrieval of ESTs matching large DNA queries and their subsequent analysis. RepeatMasker is used to mask dispersed repetitive sequences (such as Alu elements) in the query, BLAST2 for searching EST databases and Artemis for graphical display of the findings. Gene2EST combines these components into a Web resource targeted at the researcher who wishes to study one or a few genes to a high level of detail.

Combinatorial association and abundance of components of interferon-stimulated gene factor 3 dictate the selectivity of interferon responses.

Genes containing the interferon-stimulated response element (ISRE) enhancer have been characterized as transcriptionally responsive primarily to type I interferons (IFN alpha/beta). Induction is due to activation of a multimeric transcription factor, interferon-stimulated gene factor 3 (ISGF3), which is activated by IFN alpha/beta but not by IFN gamma. We found that ISRE-containing genes were induced by IFN gamma as well as by IFN alpha in Vero cells. The IFN gamma response was dependent on the ISRE and was accentuated by preexposure of cells to IFN alpha, a treatment that increases the abundance of ISGF3 components. Overexpression of ISGF3 polypeptides showed that the IFN gamma response depended on the DNA-binding protein ISGF3 gamma (p48) as well as on the 91-kDa protein STAT91 (Stat1 alpha). The transcriptional response to IFN alpha required the 113-kDa protein STAT113 (Stat2) in addition to STAT91 and p48. Mutant fibrosarcoma cells deficient in each component of ISGF3 were used to confirm that IFN gamma induction of an ISRE reporter required p48 and STAT91, but not STAT113. A complex containing p48 and phosphorylated STAT91 but lacking STAT113 bound the ISRE in vitro. IFN gamma-induced activation of this complex, preferentially formed at high concentrations of p48 and STAT91, may explain some of the overlapping responses to IFN alpha and IFN gamma.

Inhibition of gene expression in human cells through small molecule-RNA interactions

Small molecules that bind their biological receptors with high affinity and selectivity can be isolated from randomized pools of combinatorial libraries. RNA-protein interactions are important in many cellular functions, including transcription, RNA splicing, and translation. One example of such interactions is the mechanism of trans-activation of HIV-1 gene expression that requires the interaction of Tat protein with the trans-activation responsive region (TAR) RNA, a 59-base stem-loop structure located at the 5′ end of all nascent HIV-1 transcripts. Here we demonstrate the isolation of small TAR RNA-binding molecules from an encoded combinatorial library. We have made an encoded combinatorial tripeptide library of 24,389 possible members from d-and l-alpha amino acids on TentaGel resin. Using on-bead screening we have identified a small family of mostly heterochiral tripeptides capable of structure-specific binding to the bulge loop of TAR RNA. In vitro binding studies reveal stereospecific discrimination when the best tripeptide ligand is compared with diastereomeric peptide sequences. In addition, the most strongly binding tripeptide was shown to suppress transcriptional activation by Tat protein in human cells with an IC50 of ≈50 nM. Our results indicate that tripeptide RNA ligands are cell permeable, nontoxic to cells, and capable of inhibiting expression of specific genes by interfering with RNA-protein interactions.

An artificial cell-cycle inhibitor isolated from a combinatorial library

Understanding the genetic networks that operate inside cells will require the dissection of interactions among network members. Here we describe a peptide aptamer isolated from a combinatorial library that distinguishes among such interactions. This aptamer binds to cyclin-dependent kinase 2 (Cdk2) and inhibits its kinase activity. In contrast to naturally occurring inhibitors, such as p21Cip1, which inhibit the activity of Cdk2 on all its substrates, inhibition by pep8 has distinct substrate specificity. We show that the aptamer binds to Cdk2 at or near its active site and that its mode of inhibition is competitive. Expression of pep8 in human cells retards their progression through the G1 phase of the cell cycle. Our results suggest that the aptamer inhibits cell-cycle progression by blocking the activity of Cdk2 on substrates needed for the G1-to-S transition. This work demonstrates the feasibility of selection of artificial proteins to perform functions not developed during evolution. The ability to select proteins that block interactions between a gene product and some partners but not others should make sophisticated genetic manipulations possible in human cells and other currently intractable systems.

The DNMT3B DNA methyltransferase gene is mutated in the ICF immunodeficiency syndrome

DNA methylation is an important regulator of genetic information in species ranging from bacteria to humans. DNA methylation appears to be critical for mammalian development because mice nullizygous for a targeted disruption of the DNMT1 DNA methyltransferase die at an early embryonic stage. No DNA methyltransferase mutations have been reported in humans until now. We describe here the first example of naturally occurring mutations in a mammalian DNA methyltransferase gene. These mutations occur in patients with a rare autosomal recessive disorder, which is termed the ICF syndrome, for immunodeficiency, centromeric instability, and facial anomalies. Centromeric instability of chromosomes 1, 9, and 16 is associated with abnormal hypomethylation of CpG sites in their pericentromeric satellite regions. We are able to complement this hypomethylation defect by somatic cell fusion to Chinese hamster ovary cells, suggesting that the ICF gene is conserved in the hamster and promotes de novo methylation. ICF has been localized to a 9-centimorgan region of chromosome 20 by homozygosity mapping. By searching for homologies to known DNA methyltransferases, we identified a genomic sequence in the ICF region that contains the homologue of the mouse Dnmt3b methyltransferase gene. Using the human sequence to screen ICF kindreds, we discovered mutations in four patients from three families. Mutations include two missense substitutions and a 3-aa insertion resulting from the creation of a novel 3′ splice acceptor. None of the mutations were found in over 200 normal chromosomes. We conclude that mutations in the DNMT3B are responsible for the ICF syndrome.

Making artificial antibodies: A format for phage display of combinatorial heterodimeric arrays

The gene VII protein (pVII) and gene IX protein (pIX) are associated closely on the surface of filamentous bacteriophage that is opposite of the end harboring the widely exploited pIII protein. We developed a phagemid format wherein antibody heavy- and light-chain variable regions were fused to the amino termini of pVII and pIX, respectively. Significantly, the fusion proteins interacted to form a functional Fv-binding domain on the phage surface. Our approach will be applicable to the display of generic peptide and protein libraries that can form combinatorial heterodimeric arrays. Consequently, it represents a first step toward artificial antibodies and the selection of novel biological activities.

Combinatorial roles for pRB, p107, and p130 in E2F-mediated cell cycle control

Numerous studies have implicated the pRB family of nuclear proteins in the control of cell cycle progression. Although over-expression experiments have revealed that each of these proteins, pRB, p107, and p130, can induce a G1 cell cycle arrest, mouse knockouts demonstrated distinct developmental requirements for these proteins, as well as partial functional redundancy between family members. To study the mechanism by which the closely related pRB family proteins contribute to cell cycle progression, we generated 3T3 fibroblasts derived from embryos that lack one or more of these proteins (pRB−/−, p107−/−, p130−/−, pRB−/−/p107−/−, pRB−/−/p130−/−, and p107−/−/p130−/−). By comparing the growth and cell cycle characteristics of these cells, we have observed clear differences in the manner in which they transit through the G1 and S phases as well as exit from the cell cycle. Deletion of Rb, or more than one of the family members, results in a shortening of G1 and a lengthening of S phase, as well as a reduction in growth factor requirements. In addition, the individual cell lines showed differential regulation of a subset of E2F-dependent gene promoters, as well as differences in cell cycle-dependent kinase activity. Taken together, these observations suggest that the closely related pRB family proteins affect cell cycle progression through distinct biochemical mechanisms and that their coordinated action may contribute to their diverse functions in various physiological settings.

Selection of peptides that functionally replace a zinc finger in the Sp1 transcription factor by using a yeast combinatorial library

We have developed a strategy for the identification of peptides able to functionally replace a zinc finger domain in a transcription factor. This strategy could have important ramifications for basic research on gene regulation and for the development of therapeutic agents. In this study in yeast, we expressed chimeric proteins that included a random peptide combinatorial library in association with two zinc finger domains and a transactivating domain. The library was screened for chimeric proteins capable of activating transcription from a target sequence in the upstream regulatory regions of selectable or reporter genes. In a screen of approximately 1.5 × 107 transformants we identified 30 chimeric proteins that exhibited transcriptional activation, some of which were able to discriminate between wild-type and mutant DNA targets. Chimeric library proteins expressed as glutathione S-transferase fusions bound to double-stranded oligonucleotides containing the target sequence, suggesting that the chimeras bind directly to DNA. Surprisingly, none of the peptides identified resembled a zinc finger or other well-known transcription factor DNA binding domain.

Analysis of endoplasmic reticulum trafficking signals by combinatorial screening in mammalian cells

To improve the accuracy of predicting membrane protein sorting signals, we developed a general methodology for defining trafficking signal consensus sequences in the environment of the living cell. Our approach uses retroviral gene transfer to create combinatorial expression libraries of trafficking signal variants in mammalian cells, flow cytometry to sort cells based on trafficking phenotype, and quantitative trafficking assays to measure the efficacy of individual signals. Using this strategy to analyze arginine- and lysine-based endoplasmic reticulum localization signals, we demonstrate that small changes in the local sequence context dramatically alter signal strength, generating a broad spectrum of trafficking phenotypes. Finally, using sequences from our screen, we found that the potency of di-lysine, but not di-arginine, mediated endoplasmic reticulum localization was correlated with the strength of interaction with α-COP.

Identification and dissection of an enhancer controlling epithelial gene expression in skin

Keratins 14 and 5 are the structural hallmarks of the basal keratinocytes of the epidermis and outer root sheath (ORS) of the hair follicle. Their genes are controlled in a tissue-specific manner and thus serve as useful tools to elucidate the regulatory mechanisms involved in keratinocyte-specific transcription. Previously we identified several keratinocyte-specific DNase I hypersensitive sites (HSs) in the 5′ regulatory sequences of the K14 gene and showed that a 700-bp regulatory domain encompassing HSs II and III can confer epidermal and ORS-specific gene expression in transgenic mice in vivo. Although HS II harbored much of the transactivation activity in vitro, it was not sufficient to restrict expression to keratinocytes in vivo. We now explore the HS III regulatory element. Surprisingly, this element on its own confers gene expression to the keratinocytes of the inner root sheath (IRS) of the hair follicle, whereas a 275-bp DNA fragment containing both HSs II and III shifts the expression from the IRS to the basal keratinocytes and ORS in vivo. Electrophoretic mobility-shift assays and mutational studies of HSs III reveal a role for CACCC-box binding proteins, Sp1 family members, and other factors adding to the list of previously described factors that are involved in keratinocyte-specific gene expression. These studies highlight a cooperative interaction of the two HSs domains and strengthen the importance of combinatorial play of transcription factors that govern keratinocyte-specific gene regulation.

Pseudoknots in prion protein mRNAs confirmed by comparative sequence analysis and pattern searching

The human prion gene contains five copies of a 24 nt repeat that is highly conserved among species. An analysis of folding free energies of the human prion mRNA, in particular in the repeat region, suggested biased codon selection and the presence of RNA patterns. In particular, pseudoknots, similar to the one predicted by Wills in the human prion mRNA, were identified in the repeat region of all available prion mRNAs available in GenBank, but not those of birds and the red slider turtle. An alignment of these mRNAs, which share low sequence homology, shows several co-variations that maintain the pseudoknot pattern. The presence of pseudoknots in yeast Sup35p and Rnq1 suggests acquisition in the prokaryotic era. Computer generated three-dimensional structures of the human prion pseudoknot highlight protein and RNA interaction domains, which suggest a possible effect in prion protein translation. The role of pseudoknots in prion diseases is discussed as individuals with extra copies of the 24 nt repeat develop the familial form of Creutzfeldt–Jakob disease.

Tbx19, a tissue-selective regulator of POMC gene expression

Pituitary cell types arise in a temporally and spatially specific fashion, in response to combinatorial actions of transcription factors induced by transient signaling gradients. The critical transcriptional determinants of the two pituitary cell types that express the pro-opiomelanocortin (POMC) gene, the anterior lobe corticotropes, producing adrenocorticotropin, and the intermediate lobe melanotropes, producing melanocyte-stimulating hormone (MSHα), have remained unknown. Here, we report that a member of the T-box gene family, Tbx19, which is expressed only in the rostral ventral diencephalon and pituitary gland, commencing on e11.5, marks pituitary cells that will subsequently express the POMC gene and is capable of altering progression of ventral cell types and inducing adrenocorticotropin in rostral tip cells. It is suggested that Tbx19, depending on the presence of synergizing transcription factors, can activate POMC gene expression and repress the α glycoprotein subunit and thyroid-stimulating hormone β promoters.

Combinatorial control of muscle development by basic helix-loop-helix and MADS-box transcription factors.

Members of the MyoD family of muscle-specific basic helix-loop-helix (bHLH) proteins function within a genetic pathway to control skeletal muscle development. Mutational analyses of these factors suggested that their DNA binding domains mediated interaction with a coregulator required for activation of muscle-specific transcription. Members of the myocyte enhancer binding factor 2 (MEF2) family of MADS-box proteins are expressed at high levels in muscle and neural cells and at lower levels in several other cell types. MEF2 factors are unable to activate muscle gene expression alone, but they potentiate the transcriptional activity of myogenic bHLH proteins. This potentiation appears to be mediated by direct interactions between the DNA binding domains of these different types of transcription factors. Biochemical and genetic evidence suggests that MEF2 factors are the coregulators for myogenic bHLH proteins. The presence of MEF2 and cell-specific bHLH proteins in other cell types raises the possibility that these proteins may also cooperate to regulate other programs of cell-specific gene expression. We present a model to account for such cooperative interactions.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.