Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

Intraembryonic hematopoietic cell migration during vertebrate development.

Vertebrate hematopoietic stem cells are derived from vental mesoderm, which is postulated to migrate to both extra- and intraembryonic positions during gastrula and neurula stages. Extraembryonic migration has previously been documented, but the origin and migration of intraembryonic hematopoietic cells have not been visualized. The zebrafish and most other teleosts do not form yolk sac blood islands during early embryogenesis, but instead hematopoiesis occurs solely in a dorsal location known as the intermediate cell mass (IM) or Oellacher. In this report, we have isolated cDNAs encoding zebrafish homologs of the hematopoietic transcription factors GATA-1 and GATA-2 and have used these markers to determine that the IM is formed from mesodermal cells in a posterior-lateral position on the yolk syncytial layer of the gastrula yolk sac. Surprisingly, cells of the IM then migrate anteriorly through most of the body length prior to the onset of active circulation and exit onto the yolk sac. These findings support a hypothesis in which the hematopoietic program of vertebrates is established by variations in homologous migration pathways of extra- and intraembryonic progenitors.

Nramp defines a family of membrane proteins.

Nramp (natural resistance-associated macrophage protein) is a newly identified family of integral membrane proteins whose biochemical function is unknown. We report on the identification of Nramp homologs from the fly Drosophila melanogaster, the plant Oryza sativa, and the yeast Saccharomyces cerevisiae. Optimal alignment of protein sequences required insertion of very few gaps and revealed remarkable sequence identity of 28% (yeast), 40% (plant), and 55% (fly) with the mammalian proteins (46%, 58%, and 73% similarity), as well as a common predicted transmembrane topology. This family is defined by a highly conserved hydrophobic core encoding 10 transmembrane segments. Other features of this hydrophobic core include several invariant charged residues, helical periodicity of sequence conservation suggesting conserved and nonconserved faces for several transmembrane helices, a consensus transport signature on the intracytoplasmic face of the membrane, and structural determinants previously described in ion channels. These characteristics suggest that the Nramp polypeptides form part of a group of transporters or channels that act on as yet unidentified substrates.

Crystal structure of the cell cycle-regulatory protein suc1 reveals a beta-hinge conformational switch.

The Schizosaccharomyces pombe cell cycle-regulatory protein suc1, named as the suppressor of cdc2 temperature-sensitive mutations, is essential for cell cycle progression. To understand suc1 structure-function relationships and to help resolve conflicting interpretations of suc1 function based on genetic studies of suc1 and its functional homologs in both lower and higher eukaryotes, we have determined the crystal structure of the beta-interchanged suc1 dimer. Each domain consists of three alpha-helices and a four-stranded beta-sheet, completed by the interchange of terminal beta-strands between the two subunits. This beta-interchanged suc1 dimer, when compared with the beta-hairpin single-domain folds of suc1, reveals a beta-hinge motif formed by the conserved amino acid sequence HVPEPH. This beta-hinge mediates the subunit conformation and assembly of suc1: closing produces the intrasubunit beta-hairpin and single-domain fold, whereas opening leads to the intersubunit beta-strand interchange and interlocked dimer assembly reported here. This conformational switch markedly changes the surface accessibility of sequence-conserved residues available for recognition of cyclin-dependent kinase, suggesting a structural mechanism for beta-hinge-mediated regulation of suc1 biological function. Thus, suc1 belongs to the family of domain-swapping proteins, consisting of intertwined and dimeric protein structures in which the dual assembly modes regulate their function.

Mutations in the MSH3 gene preferentially lead to deletions within tracts of simple repetitive DNA in Saccharomyces cerevisiae.

Eukaryotic genomes contain tracts of DNA in which a single base or a small number of bases are repeated (microsatellites). Mutations in the yeast DNA mismatch repair genes MSH2, PMS1, and MLH1 increase the frequency of mutations for normal DNA sequences and destabilize microsatellites. Mutations of human homologs of MSH2, PMS1, and MLH1 also cause microsatellite instability and result in certain types of cancer. We find that a mutation in the yeast gene MSH3 that does not substantially affect the rate of spontaneous mutations at several loci increases microsatellite instability about 40-fold, preferentially causing deletions. We suggest that MSH3 has different substrate specificities than the other mismatch repair proteins and that the human MSH3 homolog (MRP1) may be mutated in some tumors with microsatellite instability.

Cloning and expression of the cDNA of chicken cation-independent mannose-6-phosphate receptor.

We cloned and sequenced the 8767-bp full-length cDNA for the chicken cation-independent mannose-6-phosphate receptor (CI-MPR), of interest because, unlike its mammalian homologs, it does not bind insulin-like growth factor II (IGF-II). The cDNA encodes a protein of 2470 aa that includes a putative signal sequence, an extracytoplasmic domain consisting of 15 homologous repeat sequences, a 23-residue transmembrane sequence, and a 161-residue cytoplasmic sequence. Overall, it shows 60% sequence identity with human and bovine CI-MPR homologs, and all but two of 122 cysteine residues are conserved. However, it shows much less homology in the N-terminal signal sequence, in repeat 11, which is proposed to contain the IGF-II-binding site in mammalian CI-MPR homologs, and in the 14-aa residue segment in the cytoplasmic sequence that has been proposed to mediate G-protein-coupled signal transduction in response to IGF-II binding by the human CI-MPR. Transient expression in COS-7 cells produced a functional CI-MPR which exhibited mannose-6-phosphate-inhibitable binding and mediated endocytosis of recombinant human beta-glucuronidase. Expression of the functional chicken CI-MPR in mice lacking the mammalian CI-MPR should clarify the controversy over the physiological role of the IGF-II-binding site in mammalian CI-MPR homologs.

Hydrophobic cluster analysis predicts an amino-terminal domain of varicella-zoster virus open reading frame 10 required for transcriptional activation.

Varicella-zoster virus open reading frame 10 (ORF10) protein, the homolog of the herpes simplex virus protein VP16, can transactivate immediate-early promoters from both viruses. A protein sequence comparison procedure termed hydrophobic cluster analysis was used to identify a motif centered at Phe-28, near the amino terminus of ORF10, that strongly resembles the sequence of the activating domain surrounding Phe-442 of VP16. With a series of GAL4-ORF10 fusion proteins, we mapped the ORF10 transcriptional-activation domain to the amino-terminal region (aa 5-79). Extensive mutagenesis of Phe-28 in GAL4-ORF10 fusion proteins demonstrated the importance of an aromatic or bulky hydrophobic amino acid at this position, as shown previously for Phe-442 of VP16. Transactivation by the native ORF10 protein was abolished when Phe-28 was replaced by Ala. Similar amino-terminal domains were identified in the VP16 homologs of other alphaherpesviruses. Hydrophobic cluster analysis correctly predicted activation domains of ORF10 and VP16 that share critical characteristics of a distinctive subclass of acidic activation domains.

An integral membrane component of coatomer-coated transport vesicles defines a family of proteins involved in budding.

We have isolated a major integral membrane protein from Golgi-derived coatomer-coated vesicles. This 24-kDa protein, p24, defines a family of integral membrane proteins with homologs present in yeast and humans. In addition to sequence similarity, all p24 family members contain a motif with the characteristic heptad repeats found in coiled coils. When the yeast p24 isoform, yp24A, is knocked out in a strain defective for vesicle fusion, a dramatic reduction in the accumulation of transport vesicles is observed. Together, these results indicate a role for this protein family in the budding of coatamer-coated and other species of coated vesicles.

Transcription in archaea: similarity to that in eucarya.

We present homologies between archaeal and eucaryal DNA-dependent RNA polymerase (RNAP) subunits and transcription factors. The sequences of the Sulfolobus acidocaldarius subunits D, E, and N and alignments with eucaryal homologs are presented here. The similarities between archaeal transcription factors and their eucaryal homologs TFIIB and TBP have been established in other laboratories. The archaeal RNAP subunits H, K, and N, respectively, show high sequence similarity to ABC27, ABC23, and ABC10 beta (found in all three eucaryal RNAPs); subunit D, to AC40 (common to polymerase II and polymerase III) and B44 (polymerase II); and subunit L, to AC19 and B12.5. The similarity of subunit D and its eucaryal homologs to bacterial alpha is limited to the "alpha-motif," which is also present in subunit L and its eucaryal homologs. Genes encoding homologs of the related eucaryal RNAP subunits A12.2/B12.6 and also homologs of eucaryal transcription elongation factors of the TFIIS family have been detected in Sulfolobus acidocaldarius and Thermococcus celer. In archaea, the protein is not an RNAP subunit. Together with the sequence similarities between archaeal box A-containing and eucaryal TATA box-containing promoters, this shows that the archaeal and eucaryal transcription systems are truly homologous and that they differ structurally and functionally from the bacterial transcription machinery. In contrast, however, a number of genes for the archaeal transcription apparatus are organized in clusters resembling the clusters of transcription-associated genes in Bacteria.

Cloning of a receptor for amphibian [Phe13]bombesin distinct from the receptor for gastrin-releasing peptide: identification of a fourth bombesin receptor subtype (BB4).

Bombesin is a tetradecapeptide originally isolated from frog skin and demonstrated to have a wide range of actions in mammals. Based on structural homology and similar biological activities, gastrin-releasing peptide (GRP) has been considered the mammalian equivalent of bombesin. We previously reported that frogs have both GRP and bombesin, which therefore are distinct peptides. We now report the cloning of a bombesin receptor subtype (BB4) that has higher affinity for bombesin than GRP. PCR was used to amplify cDNAs related to the known bombesin receptors from frog brain. Sequence analysis of the amplified cDNAs revealed 3 classes of receptor subtypes. Based on amino acid homology, two classes were clearly the amphibian homologs of the GRP and neuromedin B receptors. The third class was unusual and a full-length clone was isolated from a Bombina orientalis brain cDNA library. Expression of the receptor in Xenopus oocytes demonstrated that the receptor responded to picomolar concentrations of [Phe13]-bombesin, the form of bombesin most prevalent in frog brain. The relative rank potency of bombesin-like peptides for this receptor was [Phe13]bombesin > [Leu13]bombesin > GRP > neuromedin B. In contrast, the rank potency for the GRP receptor is GRP > [Leu13]bombesin > [Phe13]bombesin > neuromedin B. Transient expression in CHOP cells gave a Ki for [Phe13]bombesin of 0.2 nM versus a Ki of 2.1 nM for GRP. Distribution analysis showed that this receptor was expressed only in brain, consistent with the distribution of [Phe13]-bombesin. Thus, based on distribution and affinity, this bombesin receptor is the receptor for [Phe13]bombesin. Phylogenetic analysis suggests that this receptor separated prior to separation of the GRP and neuromedin B receptors; thus, BB4 receptors and their cognate ligands may also exist in mammals.

Immunolocalization of the mercurial-insensitive water channel and glycerol intrinsic protein in epithelial cell plasma membranes.

Two water channel homologs were cloned recently from rat kidney, mercurial-insensitive water channel (MIWC) and glycerol intrinsic protein (GLIP). Polyclonal antibodies were raised against synthetic C-terminal peptides and purified by affinity chromatography. MIWC and GLIP antibodies recognized proteins in rat kidney with an apparent molecular mass of 30 and 27 kDa, respectively, and did not cross-react. By immunofluorescence, MIWC and GLIP were expressed together on the basolateral plasma membrane of collecting duct principal cells in kidney. By immunohistochemistry, MIWC and GLIP were expressed on tracheal epithelial cells with greater expression of GLIP on the basal plasma membrane and MIWC on the lateral membrane; only MIWC was expressed in bronchial epithelia. In eye, GLIP was expressed in conjunctival epithelium, whereas MIWC was found in iris, ciliary body, and neural cell layers in retina. MIWC and GLIP colocalized on the basolateral membrane of villus epithelial cells in colon and brain ependymal cells. Expression of MIWC and GLIP was not detected in small intestine, liver, spleen, endothelia, and cells that express water channels CHIP28 or WCH-CD. These studies suggest water/solute transporting roles for MIWC and GLIP in the urinary concentrating mechanism, cerebrospinal fluid absorption, ocular fluid balance, fecal dehydration, and airway humidification. The unexpected membrane colocalization of MIWC and GLIP in several tissues suggests an interaction at the molecular and/or functional levels.

Anterior axis duplication in Xenopus induced by the over-expression of the cadherin-binding protein plakoglobin.

Plakoglobin interacts with both classical and desmosomal cadherins. It is closely related to Drosophila aramadillo (arm) gene product; arm acts in the wingless (wg)-signaling pathway to establish segment polarity. In Xenopus, homologs of wg--i.e., wnts, can produce anterior axis duplications by inducing dorsal mesoderm. Studies in Drosophila suggest that wnt acts by increasing the level of cytoplasmic armadillo protein (arm). To test whether simply increasing the level of plakoglobin mimics the effects of exogenous wnts in Xenopus, we injected fertilized eggs with RNA encoding an epitope-tagged form of plakoglobin; this induced both early radial gastrulation and anterior axis duplication. Exogenous plakoglobin accumulates in the nuclei of embryonic cells. Plakoglobin binds to the tail domain of the desmosomal cadherin desmoglein 1. When RNA encoding the tail domain of desmoglein was coinjected with plakoglobin RNA, both the dorsalizing effect and nuclear accumulation of plakoglobin were suppressed. Mutational analysis indicates that the central arm repeat region of plakoglobin is sufficient to induce axis duplication and that this polypeptide accumulates in the nuclei of embryonic cells. These data show that increased plakoglobin levels can, by themselves, generate the intracellular signals involved in the specification of dorsal mesoderm.

Delivery of antisense oligodeoxyribonucleotides against the human epidermal growth factor receptor into cultured KB cells with liposomes conjugated to folate via polyethylene glycol.

Antisense oligodeoxyribonucleotides targeted to the epidermal growth factor (EGF) receptor were encapsulated into liposomes linked to folate via a polyethylene glycol spacer (folate-PEG-liposomes) and efficiently delivered into cultured KB cells via folate receptor-mediated endocytosis. The oligonucleotides were a phosphodiester 15-mer antisense to the EGF receptor (EGFR) gene stop codon (AEGFR2), the same sequence with three phosphorothioate linkages at each terminus (AEGFR2S), a randomized 15-mer control of similar base composition to AEGFR2 (RC15), a 14-mer control derived from a symmetrized Escherichia coli lac operator (LACM), and the 5'-fluorescein-labeled homologs of several of the above. Cellular uptake of AEGFR2 encapsulated in folate-PEG-liposomes was nine times higher than AEGFR2 encapsulated in nontargeted liposomes and 16 times higher than unencapsulated AEGFR2. Treatment of KB cells with AEGFR2 in folate-PEG-liposomes resulted in growth inhibition and significant morphological changes. Curiously, AEGFR2 and AEGFR2S encapsulated in folate-PEG-liposomes exhibited virtually identical growth inhibitory effects, reducing KB cell proliferation by > 90% 48 hr after the cells were treated for 4 hr with 3 microM oligonucleotide. Free AEGFR2 caused almost no growth inhibition, whereas free AEGFR2S was only one-fifth as potent as the folate-PEG-liposome-encapsulated oligonucleotide. Growth inhibition of the oligonucleotide-treated cells was probably due to reduced EGFR expression because indirect immunofluorescence staining of the cells with a monoclonal antibody against the EGFR showed an almost quantitative reduction of the EGFR in cells treated with folate-PEG-liposome-entrapped AEGFR2. These results suggest that antisense oligonucleotide encapsulation in folate-PEG-liposomes promise efficient and tumor-specific delivery and that phosphorothioate oligonucleotides appear to offer no major advantage over native phosphodiester DNA when delivered by this route.

The sulfur controller-2 negative regulatory gene of Neurospora crassa encodes a protein with beta-transducin repeats.

The sulfur regulatory system of Neurospora crassa is composed of a set of structural genes involved in sulfur catabolism controlled by a genetically defined set of trans-acting regulatory genes. These sulfur regulatory genes include cys-3+, which encodes a basic region-leucine zipper transcriptional activator, and the negative regulatory gene scon-2+. We report here that the scon-2+ gene encodes a polypeptide of 650 amino acids belonging to the expanding beta-transducin family of eukaryotic regulatory proteins. Specifically, SCON2 protein contains six repeated G beta-homologous domains spanning the C-terminal half of the protein. SCON2 represents the initial filamentous fungal protein identified in the beta-transducin group. Additionally, SCON2 exhibits a specific amino-terminal domain that potentially defines another subfamily of beta-transducin homologs. Expression of the scon-2+ gene has been examined using RNA hybridization and gel mobility-shift analysis. The dependence of scon-2+ expression on CYS3 function and the binding of CYS3 to the scon-2+ promoter indicate the presence of an important control loop within the N. crassa sulfur regulatory circuit involving CYS3 activation of scon-2+ expression. On the basis of the presence of beta-transducin repeats, the crucial role of SCON2 in the signal-response pathway triggered by sulfur limitation may be mediated by protein-protein interactions.

Root of the universal tree of life based on ancient aminoacyl-tRNA synthetase gene duplications.

Universal trees based on sequences of single gene homologs cannot be rooted. Iwabe et al. [Iwabe, N., Kuma, K.-I., Hasegawa, M., Osawa, S. & Miyata, T. (1989) Proc. Natl. Acad. Sci. USA 86, 9355-9359] circumvented this problem by using ancient gene duplications that predated the last common ancestor of all living things. Their separate, reciprocally rooted gene trees for elongation factors and ATPase subunits showed Bacteria (eubacteria) as branching first from the universal tree with Archaea (archaebacteria) and Eucarya (eukaryotes) as sister groups. Given its topical importance to evolutionary biology and concerns about the appropriateness of the ATPase data set, an evaluation of the universal tree root using other ancient gene duplications is essential. In this study, we derive a rooting for the universal tree using aminoacyl-tRNA synthetase genes, an extensive multigene family whose divergence likely preceded that of prokaryotes and eukaryotes. An approximately 1600-bp conserved region was sequenced from the isoleucyl-tRNA synthetases of several species representing deep evolutionary branches of eukaryotes (Nosema locustae), Bacteria (Aquifex pyrophilus and Thermotoga maritima) and Archaea (Pyrococcus furiosus and Sulfolobus acidocaldarius). In addition, a new valyl-tRNA synthetase was characterized from the protist Trichomonas vaginalis. Different phylogenetic methods were used to generate trees of isoleucyl-tRNA synthetases rooted by valyl- and leucyl-tRNA synthetases. All isoleucyl-tRNA synthetase trees showed Archaea and Eucarya as sister groups, providing strong confirmation for the universal tree rooting reported by Iwabe et al. As well, there was strong support for the monophyly (sensu Hennig) of Archaea. The valyl-tRNA synthetase gene from Tr. vaginalis clustered with other eukaryotic ValRS genes, which may have been transferred from the mitochondrial genome to the nuclear genome, suggesting that this amitochondrial trichomonad once harbored an endosymbiotic bacterium.