BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

Molecular characteristics of fibroblast growth factor–fibroblast growth factor receptor–heparin-like glycosaminoglycan complex

Fibroblast growth factor (FGF) family plays key roles in development, wound healing, and angiogenesis. Understanding of the molecular nature of interactions of FGFs with their receptors (FGFRs) has been seriously limited by the absence of structural information on FGFR or FGF–FGFR complex. In this study, based on an exhaustive analysis of the primary sequences of the FGF family, we determined that the residues that constitute the primary receptor-binding site of FGF-2 are conserved throughout the FGF family, whereas those of the secondary receptor binding site of FGF-2 are not. We propose that the FGF–FGFR interaction mediated by the ‘conserved’ primary site interactions is likely to be similar if not identical for the entire FGF family, whereas the ‘variable’ secondary sites, on both FGF as well as FGFR mediates specificity of a given FGF to a given FGFR isoform. Furthermore, as the pro-inflammatory cytokine interleukin 1 (IL-1) and FGF-2 share the same structural scaffold, we find that the spatial orientation of the primary receptor-binding site of FGF-2 coincides structurally with the IL-1β receptor-binding site when the two molecules are superimposed. The structural similarities between the IL-1 and the FGF system provided a framework to elucidate molecular principles of FGF–FGFR interactions. In the FGF–FGFR model proposed here, the two domains of a single FGFR wrap around a single FGF-2 molecule such that one domain of FGFR binds to the primary receptor-binding site of the FGF molecule, while the second domain of the same FGFR binds to the secondary receptor-binding site of the same FGF molecule. Finally, the proposed model is able to accommodate not only heparin-like glycosaminoglycan (HLGAG) interactions with FGF and FGFR but also FGF dimerization or oligomerization mediated by HLGAG.

A cobalt complex that selectively disrupts the structure and function of zinc fingers

Zinc finger domains are structures that mediate sequence recognition for a large number of DNA-binding proteins. These domains consist of sequences of amino acids containing cysteine and histidine residues tetrahedrally coordinated to a zinc ion. In this report, we present a means to selectively inhibit a zinc finger transcription factor with cobalt(III) Schiff-base complexes. 1H NMR spectroscopy confirmed that the structure of a zinc finger peptide is disrupted by axial ligation of the cobalt(III) complex to the nitrogen of the imidazole ring of a histidine residue. Fluorescence studies reveal that the zinc ion is displaced from the model zinc finger peptide in the presence of the cobalt complex. In addition, gel-shift and filter-binding assays reveal that cobalt complexes inhibit binding of a complete zinc finger protein, human transcription factor Sp1, to its consensus sequence. Finally, a DNA-coupled conjugate of the cobalt complexes selectively inhibited Sp1 in the presence of several other transcription factors.

Effective inhibition of influenza virus production in cultured cells by external guide sequences and ribonuclease P

The polymerase (PB2) and nucleocapsid (NP) genes encoded by the genome of influenza virus are essential for replication of the virus. When synthetic genes that express RNAs for external guide sequences targeted to the mRNAs of the PB2 and NP genes are stably incorporated into mouse cells in tissue culture, infection of these cells with influenza virus is nonproductive. Endogenous RNase P cleaves the targeted influenza virus mRNAs when they are in a complex with the external guide sequences. Targeting two different mRNAs simultaneously inhibits viral particle production more efficiently than does targeting only one mRNA.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd′ of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd′ suggests that the inversion may be the result of the transposable element’s activity. Characteristics of Odysseus’ terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

Diversification, expression, and γδ T cell recognition of evolutionarily distant members of the MIC family of major histocompatibility complex class I-related molecules

Distant relatives of major histocompatibility complex (MHC) class I molecules, human MICA and MICB, function as stress-induced antigens that are broadly recognized by intestinal epithelial γδ T cells. They may thus play a central role in the immune surveillance of damaged, infected, or otherwise stressed intestinal epithelial cells. However, the generality of this system in evolution and the mode of recognition of MICA and MICB are undefined. Analysis of cDNA sequences from various primate species defined translation products that are homologous to MICA and MICB. All of the MIC polypeptides have common characteristics, although they are extraordinarily diverse. The most notable alterations are several deletions and frequent amino acid substitutions in the putative α-helical regions of the α1α2 domains. However, the primate MIC molecules were expressed on the surfaces of normal and transfected cells. Moreover, despite their sharing of relatively few identical amino acids in potentially accessible regions of their α1α2 domains, they were recognized by diverse human intestinal epithelial γδ T cells that are restricted by MICA and MICB. Thus, MIC molecules represent a family of MHC proteins that are structurally diverse yet appear to be functionally conserved. The promiscuous mode of γδ T cell recognition of these antigens may be explained by their sharing of a single conserved interaction site.

Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination

One-third of humans are infected with Mycobacterium tuberculosis, the causative agent of tuberculosis. Sequence analysis of two megabases in 26 structural genes or loci in strains recovered globally discovered a striking reduction of silent nucleotide substitutions compared with other human bacterial pathogens. The lack of neutral mutations in structural genes indicates that M. tuberculosis is evolutionarily young and has recently spread globally. Species diversity is largely caused by rapidly evolving insertion sequences, which means that mobile element movement is a fundamental process generating genomic variation in this pathogen. Three genetic groups of M. tuberculosis were identified based on two polymorphisms that occur at high frequency in the genes encoding catalase-peroxidase and the A subunit of gyrase. Group 1 organisms are evolutionarily old and allied with M. bovis, the cause of bovine tuberculosis. A subset of several distinct insertion sequence IS6110 subtypes of this genetic group have IS6110 integrated at the identical chromosomal insertion site, located between dnaA and dnaN in the region containing the origin of replication. Remarkably, study of ≈6,000 isolates from patients in Houston and the New York City area discovered that 47 of 48 relatively large case clusters were caused by genotypic group 1 and 2 but not group 3 organisms. The observation that the newly emergent group 3 organisms are associated with sporadic rather than clustered cases suggests that the pathogen is evolving toward a state of reduced transmissability or virulence.

A peptide derived from an extracellular domain selectively inhibits receptor internalization: Target sequences on insulin and insulin-like growth factor 1 receptors

Certain peptides derived from the α1 domain of the major histocompatibility class I antigen complex (MHC-I) inhibit receptor internalization, increasing the steady-state number of active receptors on the cell surface and thereby enhancing the sensitivity to hormones and other agonists. These peptides self-assemble, and they also bind to MHC-I at the same site from which they are derived, suggesting that they could bind to receptor sites with significant sequence similarity. Receptors affected by MHC-I peptides do, indeed, have such sequence similarity, as illustrated here by insulin receptor (IR) and insulin-like growth factor-1 receptor. A synthetic peptide with sequence identical to a certain extracellular receptor domain binds to that receptor in a ligand-dependent manner and inhibits receptor internalization. Moreover, each such peptide is selective for its cognate receptor. An antibody to the IR peptide not only binds to IR and competes with the peptide but also inhibits insulin-dependent internalization of IR. These observations, and binding studies with deletion mutants of IR, indicate that the sequence QILKELEESSF encoded by exon 10 plays a key role in IR internalization. Our results illustrate a principle for identifying receptor-specific sites of importance for receptor internalization, and for enhancing sensitivity to hormones and other agonists.

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.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

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.

Subunit structure of the mammalian exocyst complex

The exocyst is a protein complex required for the late stages of secretion in yeast. Unlike the SNAREs (SNAP receptors), important secretory proteins that are broadly distributed on the target membrane, the exocyst is specifically located at sites of vesicle fusion. We have isolated cDNAs encoding the rexo70, rsec5, and rsec15 subunits of the mammalian complex. The amino acid sequences encoded by these genes are between 21% and 24% identical to their yeast homologs. All three genes are broadly expressed and multiple transcripts are observed for rexo70 and rsec15. Characterization of cDNAs encoding the 84-kDa subunit of the mammalian complex revealed a novel protein. mAbs were generated to the mammalian rsec6 subunit of the exocyst complex. rsec6 immunoreactivity is found in a punctate distribution at terminals of PC12 cell processes at or near sites of granule exocytosis.

Evolutionary instability of the major histocompatibility complex class I loci in New World primates

Homologues of the human major histocompatibility complex (MHC) HLA-A, -B, -E, -F, and -G loci are present in all the Catarrhini (Old World primates, apes, and humans), and some of their allelic lineages have survived several speciation events. Analysis of 26 MHC class I cDNAs from seven different genera of New World primates revealed that the Callitrichinae (tamarins and marmosets) are an exception to these rules of MHC stability. In gene trees of primate MHC class I genes, sequences from the Callitrichinae cluster in a genus-specific fashion, whereas in the other genera of New World primates, as in the Catarrhini, they cluster in a transgeneric way. The genus-specific clustering of the Callitrichinae cDNAs indicates that there is no orthology between MHC class I loci in genera of this phyletic group. Additionally, the Callitrichinae genera exhibit limited variability of their MHC class I genes, in contrast to the high variability displayed by all other primates. Each Callitrichinae genus, therefore, expresses its own set of MHC class I genes, suggesting that an unusually high rate of turnover of loci occurs in this subfamily. The limited variability of MHC class I genes in the Callitrichinae is likely the result of the recent origin of these loci.

Signal Recognition Particle-dependent Targeting of Ribosomes to the Rough Endoplasmic Reticulum in the Absence and Presence of the Nascent Polypeptide-associated Complex

Proteins with RER-specific signal sequences are cotranslationally translocated across the rough endoplasmic reticulum through a proteinaceous channel composed of oligomers of the Sec61 complex. The Sec61 complex also binds ribosomes with high affinity. The dual function of the Sec61 complex necessitates a mechanism to prevent signal sequence-independent binding of ribosomes to the translocation channel. We have examined the hypothesis that the signal recognition particle (SRP) and the nascent polypeptide-associated complex (NAC), respectively, act as positive and negative regulatory factors to mediate the signal sequence-specific attachment of the ribosome-nascent chain complex (RNC) to the translocation channel. Here, SRP-independent translocation of a nascent secretory polypeptide was shown to occur in the presence of endogenous wheat germ or rabbit reticulocyte NAC. Furthermore, SRP markedly enhanced RNC binding to the translocation channel irrespective of the presence of NAC. Binding of RNCs, but not SRP-RNCs, to the Sec61 complex is competitively inhibited by 80S ribosomes. Thus, the SRP-dependent targeting pathway provides a mechanism for delivery of RNCs to the translocation channel that is not inhibited by the nonselective interaction between the ribosome and the Sec61 complex.

Nuclear Localization of Schizosaccharomyces pombe Mcm2/Cdc19p Requires MCM Complex Assembly

The minichromosome maintenance (MCM) proteins MCM2–MCM7 are conserved eukaryotic replication factors that assemble in a heterohexameric complex. In fission yeast, these proteins are nuclear throughout the cell cycle. In studying the mechanism that regulates assembly of the MCM complex, we analyzed the cis and trans elements required for nuclear localization of a single subunit, Mcm2p. Mutation of any single mcm gene leads to redistribution of wild-type MCM subunits to the cytoplasm, and this redistribution depends on an active nuclear export system. We identified the nuclear localization signal sequences of Mcm2p and showed that these are required for nuclear targeting of other MCM subunits. In turn, Mcm2p must associate with other MCM proteins for its proper localization; nuclear localization of MCM proteins thus requires assembly of MCM proteins in a complex. We suggest that coupling complex assembly to nuclear targeting and retention ensures that only intact heterohexameric MCM complexes remain nuclear.