SMART, a simple modular architecture research tool: Identification of signaling domains

Accurate multiple alignments of 86 domains that occur in signaling proteins have been constructed and used to provide a Web-based tool (SMART: simple modular architecture research tool) that allows rapid identification and annotation of signaling domain sequences. The majority of signaling proteins are multidomain in character with a considerable variety of domain combinations known. Comparison with established databases showed that 25% of our domain set could not be deduced from SwissProt and 41% could not be annotated by Pfam. SMART is able to determine the modular architectures of single sequences or genomes; application to the entire yeast genome revealed that at least 6.7% of its genes contain one or more signaling domains, approximately 350 greater than previously annotated. The process of constructing SMART predicted (i) novel domain homologues in unexpected locations such as band 4.1-homologous domains in focal adhesion kinases; (ii) previously unknown domain families, including a citron-homology domain; (iii) putative functions of domain families after identification of additional family members, for example, a ubiquitin-binding role for ubiquitin-associated domains (UBA); (iv) cellular roles for proteins, such predicted DEATH domains in netrin receptors further implicating these molecules in axonal guidance; (v) signaling domains in known disease genes such as SPRY domains in both marenostrin/pyrin and Midline 1; (vi) domains in unexpected phylogenetic contexts such as diacylglycerol kinase homologues in yeast and bacteria; and (vii) likely protein misclassifications exemplified by a predicted pleckstrin homology domain in a Candida albicans protein, previously described as an integrin.

A single-headed dimer of Escherichia coli ribosomal protein L7/L12 supports protein synthesis

During protein synthesis, the two elongation factors Tu and G alternately bind to the 50S ribosomal subunit at a site of which the protein L7/L12 is an essential component. L7/L12 is present in each 50S subunit in four copies organized as two dimers. Each dimer consists of distinct domains: a single N-terminal (“tail”) domain that is responsible for both dimerization and binding to the ribosome via interaction with the protein L10 and two independent globular C-terminal domains (“heads”) that are required for binding of elongation factors to ribosomes. The two heads are connected by flexible hinge sequences to the N-terminal domain. Important questions concerning the mechanism by which L7/L12 interacts with elongation factors are posed by us in response to the presence of two dimers, two heads per dimer, and their dynamic, mobile properties. In an attempt to answer these questions, we constructed a single-headed dimer of L7/L12 by using recombinant DNA techniques and chemical cross-linking. This chimeric molecule was added to inactive core particles lacking wild-type L7/L12 and shown to restore activity to a level approaching that of wild-type two-headed L7/L12.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.

The active digestion of uniparental chloroplast DNA in a single zygote of Chlamydomonas reinhardtii is revealed by using the optical tweezer

The non-Mendelian inheritance of organelle genes is a phenomenon common to almost all eukaryotes, and in the isogamous alga Chlamydomonas reinhardtii, chloroplast (cp) genes are transmitted from the mating type positive (mt+) parent. In this study, the preferential disappearance of the fluorescent cp nucleoids of the mating type negative (mt−) parent was observed in living young zygotes. To study the change in cpDNA molecules during the preferential disappearance, the cpDNA of mt+ or mt− origin was labeled separately with bacterial aadA gene sequences. Then, a single zygote with or without cp nucleoids was isolated under direct observation by using optical tweezers and investigated by nested PCR analysis of the aadA sequences. This demonstrated that cpDNA molecules are digested completely during the preferential disappearance of mt− cp nucleoids within 10 min, whereas mt+ cpDNA and mitochondrial DNA are protected from the digestion. These results indicate that the non-Mendelian transmission pattern of organelle genes is determined immediately after zygote formation.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Nonequilibrium mechanism of transcription termination from observations of single RNA polymerase molecules

Cessation of transcription at specific terminator DNA sequences is used by viruses, bacteria, and eukaryotes to regulate the expression of downstream genes, but the mechanisms of transcription termination are poorly characterized. To elucidate the kinetic mechanism of termination at the intrinsic terminators of enteric bacteria, we observed, by using single-molecule light microscopy techniques, the behavior of surface-immobilized Escherichia coli RNA polymerase (RNAP) molecules in vitro. An RNAP molecule remains at a canonical intrinsic terminator for ≈64 s before releasing DNA, implying the formation of an elongation-incompetent (paused) intermediate by transcription complexes that terminate but not by those that read through the terminator. Analysis of pause lifetimes establishes a complete minimal mechanism of termination in which paused intermediate formation is both necessary and sufficient to induce release of RNAP at the terminator. The data suggest that intrinsic terminators function by a nonequilibrium process in which terminator effectiveness is determined by the relative rates of nucleotide addition and paused state entry by the transcription complex.

In silico detection of control signals: mRNA 3′-end-processing sequences in diverse species

We have investigated mRNA 3′-end-processing signals in each of six eukaryotic species (yeast, rice, arabidopsis, fruitfly, mouse, and human) through the analysis of more than 20,000 3′-expressed sequence tags. The use and conservation of the canonical AAUAAA element vary widely among the six species and are especially weak in plants and yeast. Even in the animal species, the AAUAAA signal does not appear to be as universal as indicated by previous studies. The abundance of single-base variants of AAUAAA correlates with their measured processing efficiencies. As found previously, the plant polyadenylation signals are more similar to those of yeast than to those of animals, with both common content and arrangement of the signal elements. In all species examined, the complete polyadenylation signal appears to consist of an aggregate of multiple elements. In light of these and previous results, we present a broadened concept of 3′-end-processing signals in which no single exact sequence element is universally required for processing. Rather, the total efficiency is a function of all elements and, importantly, an inefficient word in one element can be compensated for by strong words in other elements. These complex patterns indicate that effective tools to identify 3′-end-processing signals will require more than consensus sequence identification.

Direct repeat sequences in the Streptomyces chitinase-63 promoter direct both glucose repression and chitin induction

The chi63 promoter directs glucose-sensitive, chitin-dependent transcription of a gene involved in the utilization of chitin as carbon source. Analysis of 5′ and 3′ deletions of the promoter region revealed that a 350-bp segment is sufficient for wild-type levels of expression and regulation. The analysis of single base changes throughout the promoter region, introduced by random and site-directed mutagenesis, identified several sequences to be important for activity and regulation. Single base changes at −10, −12, −32, −33, −35, and −37 upstream of the transcription start site resulted in loss of activity from the promoter, suggesting that bases in these positions are important for RNA polymerase interaction. The sequences centered around −10 (TATTCT) and −35 (TTGACC) in this promoter are, in fact, prototypical of eubacterial promoters. Overlapping the RNA polymerase binding site is a perfect 12-bp direct repeat sequence. Some base changes within this direct repeat resulted in constitutive expression, suggesting that this sequence is an operator for negative regulation. Other base changes resulted in loss of glucose repression while retaining the requirement for chitin induction, suggesting that this sequence is also involved in glucose repression. The fact that cis-acting mutations resulted in glucose resistance but not inducer independence rules out the possibility that glucose repression acts exclusively by inducer exclusion. The fact that mutations that affect glucose repression and chitin induction fall within the same direct repeat sequence module suggests that the direct repeat sequence facilitates both chitin induction and glucose repression.

Recombinant immunotoxins specific for a mutant epidermal growth factor receptor: Targeting with a single chain antibody variable domain isolated by phage display

EGFRvIII is a mutant epidermal growth factor receptor found in glioblastoma, and in carcinoma of the breast, ovary, and lung. The mutant receptor has a deletion in its extracellular domain that results in the formation of a new, tumor-specific extracellular sequence. Mice were immunized with a synthetic peptide corresponding to this sequence and purified EGFRvIII. A single chain antibody variable domain (scFv) phage display library of 8 × 106 members was made from the spleen of one immunized mouse. A scFv specific for EGFRvIII was isolated from this library by panning with successively decreasing amounts of synthetic peptide. This was used to make an immunotoxin by fusing the scFv DNA sequence to sequences coding for domains II and III of Pseudomonas exotoxin A. Purified immunotoxin had a Kd of 22 nM for peptide and a Kd of 11 nM for cell-surface EGFRvIII. The immunotoxin was very cytotoxic to cells expressing EGFRvIII, with an IC50 of 1 ng/ml (16 pM) on mouse fibroblasts transfected with EGFRvIII and an IC50 of 7–10 ng/ml (110–160 pM) on transfected glioblastoma cells. There was no cytotoxic activity at 1000 ng/ml on the untransfected parent glioblastoma cell line. The immunotoxin was completely stable upon incubation at 37°C for 24 h in human serum. The combination of good affinity, cytotoxicity and stability make this immunotoxin a candidate for further preclinical evaluation.

Mycoparasitic interaction relieves binding of the Cre1 carbon catabolite repressor protein to promoter sequences of the ech42 (endochitinase-encoding) gene in Trichoderma harzianum

The fungus Trichoderma harzianum is a potent mycoparasite of various plant pathogenic fungi. We have studied the molecular regulation of mycoparasitism in the host/mycoparasite system Botrytis cinerea/T. harzianum. Protein extracts, prepared from various stages of mycoparasitism, were used in electrophoretic mobility-shift assays (EMSAs) with two promoter fragments of the ech-42 (42-kDa endochitinase-encoding) gene of T. harzianum. This gene was chosen as a model because its expression is triggered during mycoparasitic interaction [Carsolio, C., Gutierrez, A., Jimenez, B., van Montagu, M. & Herrera-Estrella, A. (1994) Proc. Natl. Acad. Sci. USA 91, 10903–10907]. All cell-free extracts formed high-molecular weight protein–DNA complexes, but those obtained from mycelia activated for mycoparasitic attack formed a complex with greater mobility. Competition experiments, using oligonucleotides containing functional and nonfunctional consensus sites for binding of the carbon catabolite repressor Cre1, provided evidence that the complex from nonmycoparasitic mycelia involves the binding of Cre1 to both fragments of the ech-42 promoter. The presence of two and three consensus sites for binding of Cre1 in the two ech-42 promoter fragments used is consistent with these findings. In contrast, the formation of the protein–DNA complex from mycoparasitic mycelia is unaffected by the addition of the competing oligonucleotides and hence does not involve Cre1. Addition of equal amounts of protein of cell-free extracts from nonmycoparasitic mycelia converted the mycoparasitic DNA–protein complex into the nonmycoparasitic complex. The addition of the purified Cre1::glutathione S-transferase protein to mycoparasitic cell-free extracts produced the same effect. These findings suggest that ech-42 expression in T. harzianum is regulated by (i) binding of Cre1 to two single sites in the ech-42 promoter, (ii) binding of a “mycoparasitic” protein–protein complex to the ech-42 promoter in vicinity of the Cre1 binding sites, and (iii) functional inactivation of Cre1 upon mycoparasitic interaction to enable the formation of the mycoparasitic protein–DNA complex.

Protective DNA vaccination against organ-specific autoimmunity is highly specific and discriminates between single amino acid substitutions in the peptide autoantigen

DNA vaccines that encode encephalitogenic sequences in tandem can protect from subsequent experimental autoimmune encephalomyelitis induced with the corresponding peptide. The mechanism for this protection and, in particular, if it is specific for the amino acid sequence encoding the vaccine are not known. We show here that a single amino acid exchange in position 79 from serine (nonself) to threonine (self) in myelin basic protein peptide MBP68–85, which is a major encephalitogenic determinant for Lewis rats, dramatically alters the protection. Moreover, vaccines encoding the encephalitogenic sequence MBP68–85 do not protect against the second encephalitogenic sequence MBP89–101 in Lewis rats and vice versa. Thus, protective immunity conferred by DNA vaccination exquisitely discriminates between peptide target autoantigens. No bystander suppression was observed. The exact underlying mechanisms remain elusive because no simple correlation between impact on ex vivo responses and protection against disease were noted.

Diversity, functionality, and stability of the T cell repertoire derived in vivo from a single human T cell precursor

In this report, we have analyzed the human T cell repertoire derived in vivo from a single T cell precursor. A unique case of X-linked severe combined immunodeficiency in which a reverse mutation occurred in an early T cell precursor was analyzed to this end. It was determined that at least 1,000 T cell clones with unique T cell receptor-β sequences were generated from this precursor. This diversity seems to be stable over time and provides protection from infections in vivo. A similar estimation was obtained in an in vitro murine model of T cell generation from a single cell precursor. Overall, our results document the large diversity potential of T cell precursors and provide a rationale for gene therapy of the block of T cell development.

Transcription coupled repair of 8-oxoguanine in murine cells: The Ogg1 protein is required for repair in nontranscribed sequences but not in transcribed sequences

To assess the role of the Ogg1 DNA glycosylase in the transcription-coupled repair (TCR) of the mutagenic lesion, 7,8-dihydro-8oxoguanine (8-OxoG), we have investigated the removal of this lesion in wild-type and ogg1−/− null mouse embryo fibroblast (MEF) cell lines. We used nonreplicating plasmids containing a single 8-OxoG·C base pair in a different assay that allowed us to study the removal of 8-OxoG located in a transcribed sequence (TS) or in a nontranscribed sequence (NTS). The results show that the removal of 8-OxoG in a wild-type MEF cell line is faster in the TS than in the NTS, indicating TCR of 8-OxoG in murine cells. In the homozygous ogg1−/− MEF cell line, 8-OxoG was not removed from the NTS whereas there was still efficient 8-OxoG repair in the TS. Expression of the mouse Ogg1 protein in the homozygous ogg1−/− cell line restored the ability to remove 8-OxoG in the NTS. Therefore, we have demonstrated that Ogg1 is essential for the repair of 8-OxoG in the NTS but is not required in the TS. These results indicate the existence of an Ogg1-independent pathway for the TCR of 8-OxoG in vivo.

Determination of gene organization in individual haplotypes by analyzing single DNA fragments from single spermatozoa

To determine human Ig heavy chain variable region (VH) gene segment organization on individual homologous chromosomes, an efficient approach has been developed. Single spermatozoa were used as subjects for the study. Upon sperm lysis, VH regions in each sperm were randomly sheared into fragments by the random Brownian force. The fragments were separated from each other by aliquoting the lysate into a certain number of tubes. The gene segments in the VH1 and VH4 families in each tube were identified by denaturing gradient gel electrophoresis after PCR amplification. The polymorphic VH sequences were used to determine the parental origins of the analyzed sperm. VH segment organization in the parental haplotypes was determined by aligning the overlapping fragments from the spermatozoa with the corresponding haplotypes. Based on this comparison between the resulting haplotype maps and the composite map reported previously, the VH region on chromosome 14 could be subdivided into four portions. The numbers and compositions of the VH gene segments differ considerably among the maps in two portions, but are highly conserved in the other two. The data also indicate that the VH region on chromosome 15 may contain a large duplicated block with copy number varying among haplotypes. The approach used in the present study may be used to construct high-resolution haplotype maps without molecular cloning.

Bidirectional imprinting of a single gene: GNAS1 encodes maternally, paternally, and biallelically derived proteins

The GNAS1 gene encodes the α subunit of the guanine nucleotide-binding protein Gs, which couples signaling through peptide hormone receptors to cAMP generation. GNAS1 mutations underlie the hormone resistance syndrome pseudohypoparathyroidism type Ia (PHP-Ia), so the maternal inheritance displayed by PHP-Ia has raised suspicions that GNAS1 is imprinted. Despite this suggestion, in most tissues Gsα is biallelically encoded. In contrast, the large G protein XLαs, also encoded by GNAS1, is paternally derived. Because the inheritance of PHP-Ia predicts the existence of maternally, rather than paternally, expressed transcripts, we have investigated the allelic origin of other mRNAs derived from GNAS1. We find this gene to be remarkable in the complexity of its allele-specific regulation. Two upstream promoters, each associated with a large coding exon, lie only 11 kb apart, yet show opposite patterns of allele-specific methylation and monoallelic transcription. The more 5′ of these exons encodes the neuroendocrine secretory protein NESP55, which is expressed exclusively from the maternal allele. The NESP55 exon is 11 kb 5′ to the paternally expressed XLαs exon. The transcripts from these two promoters both splice onto GNAS1 exon 2, yet share no coding sequences. Despite their structural unrelatedness, the encoded proteins, of opposite allelic origin, both have been implicated in regulated secretion in neuroendocrine tissues. Remarkably, maternally (NESP55), paternally (XLαs), and biallelically (Gsα) derived proteins all are produced by different patterns of promoter use and alternative splicing of GNAS1, a gene showing simultaneous imprinting in both the paternal and maternal directions.