Similar amino acid sequences: chance or common ancestry?

Characterization of a new oriental-mustard (Brassica juncea) allergen, Bra j IE: detection of an allergenic epitope

Bra j IE, a major allergen from oriental-mustard (Brassica juncea) seeds, has been isolated and characterized. Its primary structure has been elucidated. This protein is composed of two chains (37 and 92 amino acids) linked by disulphide bridges. The amino acid sequence obtained is closely related to that previously determined for Sin a I, an allergen isolated from yellow mustard (Sinapis alba). A common epitope has been detected in the large chain of both Bra j IE and Sin a I by means of electroblotting and immunodetection with 2B3, which is a monoclonal antibody raised against the yellow-mustard allergen. A histidine residue of the large chain of both mustard allergens has been found to be essential for the recognition by 2B3 antibody. A synthetic multiantigenic peptide containing this His was recognized by 2B3 as well as by sera of mustard-hypersensitive individuals. Therefore this antigenic determinant must be involved in the allergenicity of these proteins.

Structural, functional, and evolutionary relationships among extracellular solute-binding receptors of bacteria

Extracellular solute-binding proteins of bacteria serve as chemoreceptors, recognition constituents of transport systems, and initiators of signal transduction pathways. Over 50 sequenced periplasmic solute-binding proteins of gram-negative bacteria and homologous extracytoplasmic lipoproteins of gram-positive bacteria have been analyzed for sequence similarities, and their degrees of relatedness have been determined. Some of these proteins are homologous to cytoplasmic transcriptional regulatory proteins of bacteria; however, with the sole exception of the vitamin B12-binding protein of Escherichia coli, which is homologous to human glutathione peroxidase, they are not demonstrably homologous to any of the several thousand sequenced eukaryotic proteins. Most of these proteins fall into eight distinct clusters as follows. Cluster 1 solute-binding proteins are specific for malto-oligosaccharides, multiple oligosaccharides, glycerol 3-phosphate, and iron. Cluster 2 proteins are specific for galactose, ribose, arabinose, and multiple monosaccharides, and they are homologous to a number of transcriptional regulatory proteins including the lactose, galactose, and fructose repressors of E. coli. Cluster 3 proteins are specific for histidine, lysine-arginine-ornithine, glutamine, octopine, nopaline, and basic amino acids. Cluster 4 proteins are specific for leucine and leucine-isoleucine-valine, and they are homologous to the aliphatic amidase transcriptional repressor, AmiC, of Pseudomonas aeruginosa. Cluster 5 proteins are specific for dipeptides and oligopeptides as well as nickel. Cluster 6 proteins are specific for sulfate, thiosulfate, and possibly phosphate. Cluster 7 proteins are specific for dicarboxylates and tricarboxylates, but these two proteins exhibit insufficient sequence similarity to establish homology. Finally, cluster 8 proteins are specific for iron complexes and possibly vitamin B12. Members of each cluster of binding proteins exhibit greater sequence conservation in their N-terminal domains than in their C-terminal domains. Signature sequences for these eight protein families are presented. The results reveal that binding proteins specific for the same solute from different bacteria are generally more closely related to each other than are binding proteins specific for different solutes from the same organism, although exceptions exist. They also suggest that a requirement for high-affinity solute binding imposes severe structural constraints on a protein. The occurrence of two distinct classes of bacterial cytoplasmic repressor proteins which are homologous to two different clusters of periplasmic binding proteins suggests that the gene-splicing events which allowed functional conversion of these proteins with retention of domain structure have occurred repeatedly during evolutionary history.(ABSTRACT TRUNCATED AT 400 WORDS)

The mitochondrial carrier family of transport proteins: structural, functional, and evolutionary relationships

Energy transduction in mitochondria requires the transport of many specific metabolites across the inner membrane of this eukaryotic organelle. We have screened the protein sequence database for proteins homologous to the mitochondrial ATP/ADP exchange carrier, and the homologous proteins found were similarly screened to ensure that all currently sequenced members of the mitochondrial carrier family (MCF) had been identified. Thirty-seven proteins were identified, 28 of which were less than 90% identical to any other sequenced member of the MCF, and the latter proteins fell into 10 clusters or subfamilies as follows: (1) ATP/ADP exchangers of mammals, plants, algae, yeast, and fungi (11 members); (2) a bovine oxoglutarate/malate exchanger (one member); (3) mammalian uncoupling carriers (five members); (4) yeast and mammalian phosphate carriers (three members); (5) MRS proteins that suppress mitochondrial splicing defects in Saccharomyces cerevisiae (two members); (6) a putative peroxysomal carrier of Candida boidinii; (7) a putative solute carrier from the protozoan, Oxytricha fallax; (8) a putative solute carrier from S. cerevisiae; (9) a putative solute carrier from Zea mays, and (10) two putative solute carriers from the mammalian thyroid gland. The specificities of proteins in clusters 5 to 10 are not known. A multiple alignment and an evolutionary tree of the 28 selected members of the MCF were constructed, thus defining the conserved residues and the phylogenetic relationships of the proteins. Hydropathy plots of the homologous regions were determined and averaged, and the average hydropathy plots were evaluated for sequence similarity. These analyses revealed that the six transmembrane spanners exhibited varying degrees of sequence conservation and hydrophilicity. These spanners, and immediately adjacent hydrophilic loop regions, were more highly conserved than other regions of these proteins. All members of the MCF appear to consist of a tripartite structure with each of the three repeated segments being about 100 residues in length. Each repeat contains two transmembrane spanners, the first being more hydrophobic with conserved glycyl and prolyl residues, the second, preceded by a highly conserved glycyl residue, being more hydrophilic with largely conserved hydrophilic residues in certain positions. Five of the six spanners are followed by the largely conserved sequence (D/E)-Hy (K/R)[- = any residue; Hy = a hydrophobic residue]. Based on both intracluster and intercluster statistical comparisons, repeats 1, 2, and 3 are homologous, but repeats 1 are more similar to each other than they are to repeats 2 or 3 or repeats 2 or 3 are to each other.(ABSTRACT TRUNCATED AT 400 WORDS)

Five major subfamilies of mariner transposable elements in insects, including the Mediterranean fruit fly, and related arthropods

We have used a PCR assay to screen 404 insects and related arthropods for mariner elements using primers corresponding to amino acids conserved between the mariner elements of Drosophila mauritiana and a moth, Hyalophora cecropia. Potential mariner elements were detected in sixty-three species, representing ten insect orders as well as a centipede and a mite. Phylogenetic analysis of the PCR fragment sequences from thirty species identifies five major subfamilies of mariners. Many species have representatives of multiple subfamilies in their genomes, and the Medfly is an extreme example with representatives of four subfamilies. Two instances of recent horizontal transfer of mariner elements include at least three species each. The widespread but sporadic distribution of mariner elements suggests they are excellent candidates for development as transformation vectors for non-drosophilids.

Analysis of the duplicated human C4/P450c21/X gene cluster

Gene duplications, deletions and rearrangements occur with an unusually high frequency in the region of the P450c21 genes encoding 21-hydroxylase. In the human genome, the locus contains at least 6 genes, oriented 5' C4A, P450c21A, XA, C4B, P450c21B, XB 3'. Sequence analysis of the XA gene, of the 5' flanking DNA of the C4A gene, and of part of the XB gene revealed that this gene cluster was duplicated by nonhomologous recombination at a CAAG tetranucleotide. The location of this duplication suggests that it may have occurred after mammalian speciation. The XA gene is abundantly expressed in the human adrenal as a stable 2.6 kb RNA, but it is not known if that RNA serves a biological function. Knowledge of the anatomy of the XA gene facilitates genetic analysis of disease-causing lesions in the P450c21B gene. Southern blotting data show that about 76% of disordered P450c21B alleles bear gene microconversions that resemble point mutations; the remaining alleles are equally distributed between gene deletions and large gene conversions.

A new subfamily of bacterial ABC-type transport systems catalyzing export of drugs and carbohydrates

Sequence comparison studies revealed that the drug resistance transporter of Streptomyces peucetius (DrrAB) and two nodulation gene products (NodIJ) of Rhizobium leguminosarum are homologous to proteins encoded by three sets of genes that comprise capsular polysaccharide export systems in gram-negative bacteria: KpsTM of Escherichia coli, BexABC of Haemophilus influenzae, and CtrDCB of Neisseria meningitidis. These five systems comprise a new subfamily within the family of ATP binding cassette (ABC)-type transporters. We have termed this subfamily the ABC-2 subfamily. For three of the systems comprising this subfamily (Drr, Nod, and Kps) only one integral membrane constituent has been identified, whereas for the other two systems (Bex and Ctr) two dis-similar integral membrane constituents have been found. This observation suggests that the transmembrane channels of ABC-2-type transporters can be formed of homo- or heterooligomers as is true of several other classes of transport systems.

The inference of evolutionary trees from molecular data

1. Procedures for multiple alignment of sequence data, subsequent phylogenetic inference, and testing of the trees derived are presented. 2. The assumptions underlying different approaches and the extent to which they are valid are discussed.

Exhaustive matching of the entire protein sequence database

The entire protein sequence database has been exhaustively matched. Definitive mutation matrices and models for scoring gaps were obtained from the matching and used to organize the sequence database as sets of evolutionarily connected components. The methods developed are general and can be used to manage sequence data generated by major genome sequencing projects. The alignments made possible by the exhaustive matching are the starting point for successful de novo prediction of the folded structures of proteins, for reconstructing sequences of ancient proteins and metabolisms in ancient organisms, and for obtaining new perspectives in structural biochemistry.

A proposed link between nitrogen and carbon metabolism involving protein phosphorylation in bacteria

We demonstrate that certain phosphoryl transfer proteins of the bacterial phosphotransferase system (PTS), the fructose- and mannitol-specific IIA proteins or domains, are homologous to a class of proteins, one of which is known to affect transcription of some of the nitrogen-regulatory sigma 54-dependent operons in Klebsiella pneumoniae. The phosphorylatable histidyl residue in the homologous PTS proteins and the consensus sequence in the vicinity of the active-site histidine are fully conserved in all members that comprise this family of proteins. A phylogenetic tree of the eight protein members of this family was constructed, and a "signature" sequence that can serve for the identification of new protein members of this family is proposed. These observations suggest that PTS-catalyzed protein phosphorylation may provide a regulatory link between carbon and nitrogen assimilation in bacteria.

Evolutionary divergence plots of homologous proteins

A simple and efficient method is described for analyzing quantitatively multiple protein sequence alignments and finding the most conserved blocks as well as the maxima of divergence within the set of aligned sequences. It consists of calculating the mean distance and the root-mean-square distance in each column of the multiple alignment, averaging the values in a window of defined length and plotting the results as a function of the position of the window. Due attention is paid to the presence of gaps in the columns. Several examples are provided, using the sequences of several cytochromes c, serine proteases, lysozymes and globins. Two distance matrices are compared, namely the matrix derived by Gribskov and Burgess from the Dayhoff matrix, and the Risler Structural Superposition Matrix. In each case, the divergence plots effectively point to the specific residues which are known to be essential for the catalytic activity of the proteins. In addition, the regions of maximum divergence are clearly delineated. Interestingly, they are generally observed in positions immediately flanking the most conserved blocks. The method should therefore be useful for delineating the peptide segments which will be good candidates for site-directed mutagenesis and for visualizing the evolutionary constraints along homologous polypeptide chains.

Comparative analysis of the beta transducin family with identification of several new members including PWP1, a nonessential gene of Saccharomyces cerevisiae that is divergently transcribed from NMT1

While investigating the expression of the Saccharomyces cerevisiae myristoyl-CoA:protein N-myristoyltransferase gene (NMT: E.C. 2.3.1.97) by Northern blot analysis, we observed another RNA transcript whose expression resembled that of NMT1 during meiosis and was derived from a gene located less than 1 kb immediately upstream of NMT1. This new gene, designated PWP1 (for periodic tryptophan protein), is divergently transcribed from NMT1 and encodes a 576-residue protein. Null mutants of PWP1 are viable, but their growth is severely retarded and steady-state levels of several cellular proteins (including at least two proteins that label with exogenous [3H]myristic acid) are drastically reduced. New methods for database searching and assessing the statistical significance of sequence similarities identify PWP1 as a member of the beta-transducin protein superfamily. Two other previously unrecognized beta-transducin-like proteins (S. cerevisiae MAK11 and D. discoideum AAC3) were also identified, and an unexpectedly high degree of sequence homology was found between a Chlamydomonas beta-like polypeptide and the C12.3 gene of chickens. A systematic and quantitative comparative analysis resulted in classifying all beta-transducin-like sequences into 11 nonorthologous families. Based on specific sequence attributes, however, not all beta-transducin-like sequences are expected to be functionally similar, and quantitative criteria for inferring functional analogies are discussed. Possible roles of repetitive tryptophan residues in proteins are also considered.

Phylogenetic continuum indicates "galaxies" in the protein universe: preliminary results on the natural group structures of proteins

The markedly nonuniform, even systematic distribution of sequences in the protein "universe" has been analyzed by methods of protein taxonomy. Mapping of the natural hierarchical system of proteins has revealed some dense cores, i.e., well-defined clusterings of proteins that seem to be natural structural groupings, possibly seeds for a future protein taxonomy. The aim was not to force proteins into more or less man-made categories by discriminant analysis, but to find structurally similar groups, possibly of common evolutionary origin. Single-valued distance measures between pairs of superfamilies from the Protein Identification Resource were defined by two chi 2-like methods on tripeptide frequencies and the variable-length subsequence identity method derived from dot-matrix comparisons. Distance matrices were processed by several methods of cluster analysis to detect phylogenetic continuum between highly divergent proteins. Only well-defined clusters characterized by relatively unique structural, intracellular environmental, organismal, and functional attribute states were selected as major protein groups, including subsets of viral and Escherichia coli proteins, hormones, inhibitors, plant, ribosomal, serum and structural proteins, amino acid synthases, and clusters dominated by certain oxidoreductases and apolar and DNA-associated enzymes. The limited repertoire of functional patterns due to small genome size, the high rate of recombination, specific features of the bacterial membranes, or of the virus cycle canalize certain proteins of viruses and Gram-negative bacteria, respectively, to organismal groups.

Evolution of the alcohol dehydrogenase (ADH) genes in yeast: characterization of a fourth ADH in Kluyveromyces lactis

Three alcohol dehydrogenase (ADH) genes have recently been characterized in the yeast Kluyveromyces lactis. We report on a fourth ADH in K. lactis (KADH II: KADH2* gene) which is highly similar to other ADHs in K. lactis and Saccharomyces cerevisiae. KADH II appears to be a cytoplasmic enzyme, and after expression of KADH2 in S. cerevisiae enzyme activity comigrated with a K. lactis ADH present in cells grown in glucose or in ethanol. KADH I was also expressed in S. cerevisiae and it comigrated with a major ADH species expressed under glucose growth conditions in K. lactis. The substrate specificities for KADH I and KADH II were shown to be more similar to that of SADH II than to SADH I. SADH I cannot efficiently utilize long chain alcohols, in contrast to other cytoplasmic yeast ADHs, presumably because of the presence of a methionine (residue 271) in its substrate binding cleft. A comparison of the DNA sequences of ADHs among K. lactis, S. cerevisiae and Schizosaccharomyces pombe suggests that the ancestral yeast species contained one cytoplasmic ADH. After divergence from S. pombe, the ADH in the ancestor to K. lactis and S. cerevisiae was duplicated, and one ADH became localized to the mitochondrion, presumably for the oxidative use of ethanol. Following the speciation of S. cerevisiae and K. lactis, the gene encoding the cytoplasmic ADH in S. cerevisiae duplicated, which resulted in the development of the SADH II protein as the primary oxidative enzyme in place of SADH III. In contrast, the K. lactis mitochondrial ADH duplicated to give rise to the highly expressed KADH3 and KADH4 genes, both of which may still play primary roles in oxidative metabolism. These data suggest that K. lactis and S. cerevisiae use different compartments for their metabolism of ethanol. Our results also indicate that the complex regulatory circuits controlling the glucose-repressible SADH2 in S. cerevisiae are a recent acquisition from regulatory networks used for the control of genes other than SADH2.

Unification of the ferritin family of proteins

Ferritin is the iron-storage protein of eukaryotic organisms. The nucleotide sequence encoding Azotobacter vinelandii bacterioferritin, a hemoprotein, was determined. The deduced amino acid sequence reveals a high degree of identity with Escherichia coli bacterioferritin and a striking similarity to eukaryotic ferritins. Moreover, derivation of a global alignment shows that virtually all key residues specifying the unique structural motifs of eukaryotic ferritin are conserved or conservatively substituted in the A. vinelandii sequence. The alignment suggests specific methionine residues as heme-binding ligands in bacterioferritins. The overall sequence similarity with conservation of key structural residues implies that all ferritins form a unified family of proteins. The results implicate ferritins as proteins potentially common to all aerobic organisms and as such useful in taxonomic classification, evolutionary analysis, and environmental monitoring.

The primary structure of Providencia rettgeri penicillin G amidase gene and its relationship to other gram negative amidases

The nucleotide sequence of Penicillin G amidase (PA,E.C.3.5.1.11) of Providencia rettgeri was determined. We aligned our P. rettgeri PA with other known Gram negative periplasmically located beta-lactam amidases. The analysis revealed a high homology with other Enterobacteric amidases (60%-65%), while with similar Pseudomonas sp. amidases the homology exceeded 25%. These homologies indicate their common ancestry.

Partial gene duplication as a cause of human disease

Tandem duplication of large regions of DNA, including duplication of whole genes, provides a substrate for genetic evolution. Tandem duplication of smaller regions involving parts of genes is now recognized as a contributor to the mutation spectrum that results in genetic disease. In this review, more than 30 unrelated partial gene duplications that have been implicated in the genesis of human genetic disease are presented and the pathogenic effects and frequency of such duplications are summarized. The mechanisms of duplication formation are analyzed with special emphasis on the molecular details of the nucleotide sequences at the duplication junctions. Evidence to date suggests that duplication may arise from either homologous (Alu-Alu) recombination or nonhomologous recombination, the latter possibly mediated by topoisomerases. For the dystrophin gene, in which most duplications have been identified, these recombination events are intrachromosomal, suggesting that unequal sister chromatid exchange is the major mechanism.

Dynamic programming algorithms for biological sequence comparison

Efficient dynamic programming algorithms are available for a broad class of protein and DNA sequence comparison problems. These algorithms require computer time proportional to the product of the lengths of the two sequences being compared [O(N2)] but require memory space proportional only to the sum of these lengths [O(N)]. Although the requirement for O(N2) time limits use of the algorithms to the largest computers when searching protein and DNA sequence databases, many other applications of these algorithms, such as calculation of distances for evolutionary trees and comparison of a new sequence to a library of sequence profiles, are well within the capabilities of desktop computers. In particular, the results of library searches with rapid searching programs, such as FASTA or BLAST, should be confirmed by performing a rigorous optimal alignment. Whereas rapid methods do not overlook significant sequence similarities, FASTA limits the number of gaps that can be inserted into an alignment, so that a rigorous alignment may extend the alignment substantially in some cases. BLAST does not allow gaps in the local regions that it reports; a calculation that allows gaps is very likely to extend the alignment substantially. Although a Monte Carlo evaluation of the statistical significance of a similarity score with a rigorous algorithm is much slower than the heuristic approach used by the RDF2 program, the dynamic programming approach should take less than 1 hr on a 386-based PC or desktop Unix workstation. For descriptive purposes, we have limited our discussion to methods for calculating similarity scores and distances that use gap penalties of the form g = rk. Nevertheless, programs for the more general case (g = q+rk) are readily available. Versions of these programs that run either on Unix workstations, IBM-PC class computers, or the Macintosh can be obtained from either of the authors.

Generalized protein tertiary structure recognition using associative memory Hamiltonians

In previous papers, a method of protein tertiary structure recognition was described based on the construction of an associative memory Hamiltonian, which encoded the amino acid sequence and the C alpha co-ordinates of a set of database proteins. Using molecular dynamics with simulated annealing, the ability of the Hamiltonian to successfully recall the structure of a protein in the memory database was successfully demonstrated, as long as the total number of database proteins did not exceed a characteristic value, called the capacity of the Hamiltonian, equal to 0.5N to 0.7N, where N is the number of amino acid residues in the protein to be recalled. In this paper, we describe the development of additional methods to increase the capacity of the Hamiltonian, including use of a more complete representation of the protein backbone and the incorporation of contextual information into the Hamiltonian through the use of secondary structure prediction. In addition, we further extend the ability of associative memory models to predict the tertiary structures of proteins not present in the protein data set, by making the Hamiltonian invariant with respect to biological symmetries that represent site mutations and insertions and deletions. The ability of the Hamiltonian to generalize from homologous proteins to an unknown protein in the presence of other unrelated proteins in the data set is demonstrated.

A new distinct group of 2 S albumins from rapeseed. Amino acid sequence of two low molecular weight napins

Two napins (nIa and nIb), isolated from Brassica napus (rapeseed) seeds, have been sequenced. The two proteins show the common structural pattern of the 2 S albumins, since they are composed of two disulfide-linked chains of different size, yet they exhibit an atypical low molecular weight (12.5 kDa vs. 14.5 kDa of the major napins). High sequence similarity has been found between these 2 proteins, but only 54% similarity can be estimated from their comparison with the 14.5 kDa major napins. Thus, nIa and nIb are considered representatives of a new distinct group of rapeseed napins since all the previously known napins exhibit 95% sequence similarity. Unexpectedly, the similarity increases when compared with the 2 S proteins from other species.

Studies of translocation catalysis

There is a symbiotic relationship between the evolution of fundamental theory and the winning of experimentally-based knowledge. The impact of the General Chemiosmotic Theory on our understanding of the nature of membrane transport processes is described and discussed. The history of experimental studies on transport catalysed by ionophore antibiotics and the membrane proteins of mitochondria and bacteria are used to illustrate the evolution of knowledge and theory. Recent experimental approaches to understanding the lactose-H+ symport protein of Escherichia coli and other sugar porters are described to show that the lack of experimental knowledge of the three-dimensional structures of the proteins currently limits the development of theories about their molecular mechanism of translocation catalysis.

Sequence and analysis of a portion of the genomes of Shope fibroma virus and malignant rabbit fibroma virus that is important for viral replication in lymphocytes

The 10.7-kb BamHI "C" restriction fragment of malignant rabbit fibroma virus (MV) contains genes that are important for its immunosuppressive activity. When this fragment is transferred to a related avirulent leporipoxvirus, Shope fibroma virus (SFV), recombinant viruses show clinical features characteristic of MV: they replicate in lymphocytes and alter immune function in vitro, induce disseminated tumors in recipient rabbits, and are immunosuppressive in vivo. The 10.7-kb BamHI "C" restriction fragment of MV was sequenced in its entirety. Its DNA sequence and the 14 ORF's derived from analyzing this sequence are discussed. Analysis of known open reading frames to which the ORF's from MV's Bam "C" fragment show homology permits us to identify some MV ORF's showing high degrees of similarity to known and postulated proteins produced by vaccinia virus. Functions for some of these vaccinia proteins are known, while functions for others are hypothetical or unknown. Further analysis of genetic determinants of MV's virulence has indicated that two overlapping restriction subfragments of the BamHI "C" fragment can transfer MV's virulent behavior to SFV. The 0.7-kb region in which these two subfragments overlap includes the C-terminus of MV orf C-7 and the N terminus of MV orf C-8. These correspond to the C- and N-termini, respectively, of SFV orf's D-9 and D-10 and to vaccinia orf's D-6 (early transcription factor) and D-7 (subunit of RNA polymerase). We sequenced the region of SFV's BamHI "D" fragment in this area and illustrate here the comparative sequences of this portion of SFV's genome and orf's. On the basis of comparisons between MV, SFV, and vaccinia in this area we discuss the potential significance of these observations.

The nucleotide sequence and genome organization of the RNA-1 segment in two bromoviruses: broad bean mottle virus and cowpea chlorotic mottle virus

The complete nucleotide sequences of the RNA-1 segments in broad bean mottle virus (BBMV) and cowpea chlorotic mottle virus (CCMV) were determined. BBMV RNA-1 consists of 3158 nucleotides and CCMV RNA-1 has 3171 nucleotides. Both BBMV and CCMV RNA-1 are capped at the 5' end but, unlike in other tricornaviruses, BBMV RNA-1 initiates with an A residue. Both BBMV and CCMV RNA-1 are monocistronic encoding for highly homologous 1a proteins of 966 and 958 amino acids, respectively. The highest homologies are clustered within two domains: the N-domain that aligns with the nsP1 Sindbis virus protein, a putative methyl transferase, and the C-domain which has a conserved nucleotide binding motif. Previous sequence comparisons suggest that the C-terminal domain may function as an NTP-dependent RNA helicase. In addition, we find that the C-domain has patterns similar to those of the reovirus and vaccinia virus guanylyl transferases. All this implies that 1a protein is a multifunctional polypeptide involved in both RNA capping and RNA polymerization processes.

Rat sulfite oxidase antibodies cross-react with two gene family-related proteins: albumin and vitamin D-binding protein

Screening lambda cDNA libraries from rat liver with antibody to native rat liver sulfite oxidase (RLSO) showed cross-reaction with two proteins that belong to the same gene family: serum albumin and vitamin D-binding protein. Antibodies raised against native RLSO or sodium dodecyl sulfate-denatured protein cross-reacted with these proteins by Western blot analysis. The relative effectiveness of RLSO antibody binding was estimated to be 1/5 for rat serum albumin and 1/10 for rat vitamin D-binding protein. This result was not caused by contaminating proteins in the RLSO used for immunization as the RLSO preparation did not react with rat serum albumin antibody. RLSO antibodies, selected for their ability to bind rat serum albumin immobilized on nitrocellulose, recognized both rat serum albumin and RLSO. RLSO antibody, with albumin-reactive antibody removed, still recognized vitamin D-binding protein, suggesting that multiple determinants specific to each protein are involved in the cross-reaction. Comparison of RLSO antibody binding to the rat and human proteins indicated that the determinants were species-specific. cDNA clones identified by screening cDNA libraries with RLSO antibody demonstrated that these determinants reside in the C-terminal domain of these proteins. These results suggest that these proteins contain some common immunological features and may be evolutionarily related.

Evolution of enzymatic regulation of prostaglandin action: novel connections to regulation of human sex and adrenal function, antibiotic synthesis and nitrogen fixation

The recent determination of the amino acid sequences of enzymes that metabolize prostaglandins and steroids has revealed interesting connections between some of these enzymes. Human placental 15-hydroxyprostaglandin dehydrogenase, which catalyzes the oxidation of the C15 alcohol on prostaglandins E2 and F2 alpha, is homologous to 11 beta-hydroxysteroid, 17 beta-hydroxysteroid, and 3 alpha, 20 beta-hydroxysteroid dehydrogenases. That is, these four enzymes are derived from a common ancestor. Moreover, enzymes important in synthesis of antibiotics and proteins synthesized by soil bacteria that form nitrogen-fixing nodules in alfalfa and soybeans are homologous to 15-hydroxyprostaglandin dehydrogenase. These homologies provide important insights into the origins of intercellular communication that is mediated by prostaglandins, steroids, and fatty acids.

Searching protein sequence libraries: comparison of the sensitivity and selectivity of the Smith-Waterman and FASTA algorithms

The sensitivity and selectivity of the FASTA and the Smith-Waterman protein sequence comparison algorithms were evaluated using the superfamily classification provided in the National Biomedical Research Foundation/Protein Identification Resource (PIR) protein sequence database. Sequences from each of the 34 superfamilies in the PIR database with 20 or more members were compared against the protein sequence database. The similarity scores of the related and unrelated sequences were determined using either the FASTA program or the Smith-Waterman local similarity algorithm. These two sets of similarity scores were used to evaluate the ability of the two comparison algorithms to identify distantly related protein sequences. The FASTA program using the ktup = 2 sensitivity setting performed as well as the Smith-Waterman algorithm for 19 of the 34 superfamilies. Increasing the sensitivity by setting ktup = 1 allowed FASTA to perform as well as Smith-Waterman on an additional 7 superfamilies. The rigorous Smith-Waterman method performed better than FASTA with ktup = 1 on 8 superfamilies, including the globins, immunoglobulin variable regions, calmodulins, and plastocyanins. Several strategies for improving the sensitivity of FASTA were examined. The greatest improvement in sensitivity was achieved by optimizing a band around the best initial region found for every library sequence. For every superfamily except the globins and immunoglobulin variable regions, this strategy was as sensitive as a full Smith-Waterman. For some sequences, additional sensitivity was achieved by including conserved but nonidentical residues in the lookup table used to identify the initial region.

Leukemia inhibitory factor receptor is structurally related to the IL-6 signal transducer, gp130

Leukemia inhibitory factor (LIF) is a cytokine with a broad range of activities that in many cases parallel those of interleukin-6 (IL-6) although LIF and IL-6 appear to be structurally unrelated. A cDNA clone encoding the human LIF receptor was isolated by expression screening of a human placental cDNA library. The LIF receptor is related to the gp130 'signal-transducing' component of the IL-6 receptor and to the G-CSF receptor, with the transmembrane and cytoplasmic regions of the LIF receptor and gp130 being most closely related. This relationship suggests a common signal transduction pathway for the two receptors and may help to explain similar biological effects of the two ligands. Murine cDNAs encoding soluble LIF receptors were isolated by cross-hybridization and share 70% amino acid sequence identity to the human sequence.

Sequence similarity between Photosystems I and II. Identification of a Photosystem I reaction center transmembrane helix that is similar to transmembrane helix IV of the D2 subunit of Photosystem II and the M subunit of the non-sulfur purple and flexible green bacteria

There are basic structural similarities between plant PS II and bacterial RCs of the Chloroflexaceae and Rhodospirillaceae. These RCs are referred to as PS II-type RCs. A similar relationship of PS I RC to PS II-type RCs has not been established. Although plant PS I and PS II RCs show structural and functional differences, they also share similarities. Therefore, the A and B polypeptides of PS I were searched for PS II D1 and D2 polypeptide-like sequences. An alignment without gaps was found between PS II-type D2/M helix IV and PS I B helix X, as well as a weaker alignment of PS II-type D1/L with PS I B helix X. No comparable alignment with PS I A was found. In the M/D2 alignment there were eight identities and some conservative substitutions in twenty nine residues. PS I B helix X appeared to contain a modified chlorophyll dimer and monomer binding site and a modified non-heme iron-quinone binding site. The conserved residue sequence was found only in RC polypeptides. The proposed chlorophyll dimer-monomer binding site was located transmembrane from the iron-sulfur cluster X binding site. The conserved residues generally are those that interact with prosthetic groups. Half of the conserved residues are located on the same side of the helix. Thus, although there are impediments to concluding firmly that PS I B helix X has a functional and evolutionary relatedness to the D2 PS II and bacterial M RC polypeptides, our analysis gives reasonable support to the idea.

Generation and analysis of nonhomologous RNA-RNA recombinants in brome mosaic virus: sequence complementarities at crossover sites

All three single-stranded RNAs of the brome mosaic virus (BMV) genome contain a highly conserved, 193-base 3' noncoding region. To study the recombination between individual BMV RNA components, barley plants were infected with a mixture of in vitro-transcribed wild-type BMV RNAs 1 and 2 and an RNA3 mutant that carried a deletion near the 3' end. This generated a population of both homologous and nonhomologous 3' recombinant BMV RNA3 variants. Sequencing revealed that these recombinants were derived by either single or double crossovers with BMV RNA1 or RNA2. The primary sequences at recombinant junctions did not show any similarity. However, they could be aligned to form double-stranded heteroduplexes. This suggested that local hybridizations among BMV RNAs may support intermolecular exchanges.

Mitochondrial gene sequences show fungal homology for Pneumocystis carinii

A 6.8 kilobase fragment of mitochondrial DNA from Pneumocystis carinii encodes for apocytochrome b, NADH dehydrogenase subunits 1, 2, 3, and 6, cytochrome oxidase subunit II, and the small subunit of ribosomal RNA. Comparative sequence analysis with a series of organisms representative of the fungal and protozoan groups shows that P. carinii has, consistently, an average similarity of 60% with the fungi but only 20% with the protozoa. The data indicate homology with the fungi for this opportunistic pathogen.

Resonant recognition model and protein topography. Model studies with myoglobin, hemoglobin and lysozyme

This study describes the further extension of the resonant recognition model for the analysis and prediction of protein--protein and protein--DNA structure/function dependencies. The model is based on the significant correlation between spectra of numerical presentations of the amino acid or nucleotide sequences of proteins and their coded biological activity. According to this physico-mathematical method, it is possible to define amino acids in the sequence which are predicted to be the most critical for protein function. Using sperm whale myoglobin, human hemoglobin and hen egg white lysozyme as model protein examples, sets of predicted amino acids, or so-called 'hot spots', have been identified within the tertiary structure. It was found for each protein that the predicted 'hot spots', which are distributed along the primary sequence, are spatially grouped in a dome-like arrangement over the active site. The identified amino acids did not correspond to the amino acid residues which are involved in the chemical reaction site of these proteins. It is thus proposed that the resonant recognition model helps to identify amino acid residues which are important for the creation of the molecular structure around the catalytic active site and also the associated physical field conditions required for biorecognition, docking of the specific substrate and full biological activity.

Solution conformations of the N-terminal CNBr fragment of glycogen phosphorylase and its interaction with calmodulin

The CNBr peptides, CBPa and CBPb, corresponding to the N-terminal 1-91 amino acid residues of glycogen-phosphorylase a and b, respectively, were purified and characterized. CD, 31P-NMR and fluorescence spectroscopy were used to assess the structural organization of the cyanogen bromide peptides in solution. The cyanogen bromide peptides yielded 21% of alpha-helical structures by CD compared to a calculated value of 36.3%. These peptides interact with calmodulin which induces measurable alpha-helices in the cyanogen bromide peptides. The helix stabilizing reagent, trifluoroethanol, induces high numbers of alpha-helices in CBP, thereby demonstrating the conformation fluidity of this peptide. The dissociation constants for calmodulin and CBP estimated by fluorescence titrations were 36.0 and 29.9 nM for CBPb in the presence of Ca2+ and EGTA, respectively. The phosphorylated residue in CBPa causes a decrease in binding interactions with calmodulin and corresponding values obtained for CBPa by fluorescence titration are 56.0 and 141.0 nM, respectively. The Ser-P-14 of CBPa is titratable, yielding a pKa = 5.45 and a Hill coefficient of 1.5. A helical wheel analysis using a computer program in PC/GENE of the CBP shows that peptide stretches in the alpha-1 and alpha-2 helices are most basic and fairly amphiphilic and therefore represent the most probable segment for CaM binding. It is this structural character of these segments which presumably confer the ability to bind CaM and facilitate some of the allosteric transitions of glycogen phosphorylase.

Sequence similarity between opioid peptide precursors and DNA-binding proteins

The opioid peptide precursors, preprodynorphin and preproenkephalin show structure similarity with a transcription factor, hunchback and the putative helix-loop-helix DNA-binding proteins, lil-1, tal and twist. Segments with similarity contain the three enkephalin sequences in preprodynorphin and one in preproenkephalin which are present within heptapeptide repeats characteristic of an alpha-helical coiled-coil structure distinctive of an amphipathic helix-loop-helix DNA-binding motif. Hunchback and the opioid prohormones also have cystein-rich regions characteristic of zinc-finger domains in common.

Evolutionary origin of thyroglobulin by duplication of esterase gene

Rat liver microsomal carboxyesterase E1 was found to have homology with five esterases and with the COOH-terminal parts of two thyroglobulins. A phylogenetic tree constructed for these proteins shows that this new superfamily has evolved from a common ancestral gene that encoded a carboxyesterase. The tree also shows that the ancestral gene already existed before the divergence of vertebrates and invertebrates and later its duplicated genes gained various kinds of esterase activity. According to the tree, one of the duplicated genes evolved into the COOH-terminal half of thyroglobulin by a gene fusion with a DNA sequence whose evolutionary origin is unknown.

Exons--original building blocks of proteins?

In a recent paper, Walter Gilbert's group has estimated the number of original exons from which all extant proteins might have been constructed. The approach used is subjected to a critical analysis here. It is shown that there are flawed assumptions about both the mechanism and generality of exon-shuffling and in the sequence comparison procedures employed, the latter failing to distinguish chance similarity from similarity due to common ancestry. These methodological errors lead to the omission of many known cases of exon-shuffling and the inclusion of others which may not be genuine. In consequence, the analysis from the Gilbert group cannot give a reliable estimate of those modules that actually participated in exon-shuffling and provides no information on the number of protein archetypes that did not participate in these processes.

The product of the KIN1 locus in Saccharomyces cerevisiae is a serine/threonine-specific protein kinase

The catalytic domain (30 kDa) of all protein kinases can be aligned for maximum homology, thereby revealing both invariant and highly conserved residues. The KIN1 locus from Saccharomyces cerevisiae was isolated by hybridization to a degenerate oligonucleotide encoding the conserved protein kinase domain, DVWSFG. The predicted amino acid sequence revealed significant homology to the catalytic domain of protein kinases. Using antibodies raised against a bacterial LacZ/KIN1 fusion protein, we have identified by immunoprecipitation the yeast KIN1 gene product as a 145,000 dalton protein (p145KIN1). In exponentially growing yeast cells, the KIN1 protein is phosphorylated primarily on serine residues. The gene product of KIN1 was shown to be a serine/threonine-specific protein kinase in immune complexes, as determined by the transfer of label from [gamma-32P]ATP to either pp145KIN1 or to an exogenously added substrate, alpha-casein. The optimal metal ion concentration in this assay was 20 mM-MnCl2. Subsequent phosphoamino acid analysis of the radiolabelled product, pp145KIN1, demonstrated that this autophosphorylation was specific for serine/threonine residues. There is no apparent difference between wild-type cells and cells containing a disrupted KIN1 gene. The biochemical characterization of protein kinases in simple eukaryotes such as yeast will aid us in determining the role of phosphorylation in cellular growth and physiology.