Conservation among HSP60 sequences in relation to structure, function, and evolution

The chaperonin HSP60 (GroEL) proteins are essential in eubacterial genomes and in eukaryotic organelles. Functional regions inferred from mutation studies and the Escherichia coli GroEL 3D crystal complexes are evaluated in a multiple alignment across 43 diverse HSP60 sequences, centering on ATP/ADP and Mg2+ binding sites, on residues interacting with substrate, on GroES contact positions, on interface regions between monomers and domains, and on residues important in allosteric conformational changes. The most evolutionary conserved residues relate to the ATP/ADP and Mg2+ binding sites. Hydrophobic residues that contribute in substrate binding are also significantly conserved. A large number of charged residues line the central cavity of the GroEL-GroES complex in the substrate-releasing conformation. These span statistically significant intra- and inter-monomer three-dimensional (3D) charge clusters that are highly conserved among sequences and presumably play an important role interacting with the substrate. Unaligned short segments between blocks of alignment are generally exposed at the outside wall of the Anfinsen cage complex. The multiple alignment reveals regions of divergence common to specific evolutionary groups. For example, rickettsial sequences diverge in the ATP/ADP binding domain and gram-positive sequences diverge in the allosteric transition domain. The evolutionary information of the multiple alignment proffers attractive sites for mutational studies.     

Functional overlap between conserved and diverged KH domains in Saccharomyces cerevisiae SCP160

The K homology (KH) domain is a remarkably versatile and highly conserved RNA-binding motif. Classical KH domains include a characteristic pattern of hydrophobic residues, a Gly-X-X-Gly (GXXG) segment, and a variable loop. KH domains typically occur in clusters, with some retaining their GXXG sequence (conserved), while others do not (diverged). As a first step towards addressing whether GXXG is essential for KH-domain function, we explored the roles of conserved and diverged KH domains in Scp160p, a multiple-KH-domain-containing protein in Saccharomyces cerevisiae. We specifically wanted to know (1) whether diverged KH domains were essential for Scp160p function, and (2) whether diverged KH domains could functionally replace conserved KH domains. To address these questions, we deleted and/or interchanged conserved and diverged KH domains of Scp160p and expressed the mutated alleles in yeast. Our results demonstrated that the answer to each question was yes. Both conserved and diverged KH domains are essential for Scp160p function, and diverged KH domains can function in place of conserved KH domains. These findings challenge the prevailing notions about the requisite features of a KH domain and raise the possibility that there may be more functional KH domains in the proteome than previously appreciated.     

Analysis of the nitrous oxide reduction genes, nosZDFYL, of Achromobacter cycloclastes.

The structural gene, nosZ, for the monomeric N2O reductase has been cloned and sequenced from the denitrifying bacterium Achromobacter cycloclastes. The nosZ gene encodes a protein of 642 amino acid residues and the deduced amino acid sequence showed homology to the previously derived sequences for the dimeric N2O reductases. The relevant DNA region of about 3.6 kbp was also sequenced and found to consist of four genes, nosDFYL based on the similarity with the N2O reduction genes of Pseudomonas stutzeri. The gene product of A. cycloclastes nosF (299 amino acid residues) has a consensus ATP-binding sequence, and the nos Y gene encodes a hydrophobic protein (273 residues) with five transmembrane segments, suggesting the similarity with an ATP-binding cassette (ABC) transporter which has two distinct domains of a highly hydrophobic region and ATP-binding sites. The nosL gene encodes a protein of 193 amino acid residues and the derived sequence showed a consensus sequence of lipoprotein modification/processing site. The expression of nosZ gene in Escherichia coli cells and the comparison of the translated sequences of the nosDFYL genes with those of bacterial transport genes for inorganic ions are discussed.     

Transfer of rps19 to the nucleus involves the gain of an RNP-binding motif which may functionally replace RPS13 in Arabidopsis mitochondria.

The discovery of disrupted rps19 genes in Arabidopsis mitochondria prompted speculation about the transfer to the nuclear compartment. We here describe the functional gene transfer of rps19 into the nucleus of Arabidopsis. Molecular cloning and sequence analysis of rps19 show that the nuclear gene encodes a long N-terminal extension. Import studies of the precursor protein indicate that only a small part of this extension is cleaved off during import. The larger part of the extension, which shows high similarity to conserved RNA-binding domains of the RNP-CS type, became part of the S19 protein. In the Escherichia coli ribosome S19 forms an RNA-binding complex as heterodimer with S13. By using immuno-analysis and import studies we show that a eubacterial-like S13 protein is absent from Arabidopsis mitochondria, and is not substituted by either a chloroplastic or a cytosolic homologue of this ribosomal protein. We therefore propose that either a highly diverged or missing RPS13 has been functionally replaced by an RNP domain that most likely derived from a glycine-rich RNA-binding protein. These results represent the first case of a functional replacement of a ribosomal protein by a common RNA-binding domain and offer a new view on the flexibility of biological systems in using well-adapted functional domains for different jobs.     

Cloning, characterization, and inactivation of the Bacillus brevis lon gene.

A gene of Bacillus brevis HPD31 analogous to the Escherichia coli lon gene has been cloned and characterized. The cloned gene (B. brevis lon gene) encodes a polypeptide of 779 amino acids with a molecular weight of 87,400 which resembles E. coli protease La, the lon gene product. Fifty-two percent of the amino acid residues of the two polypeptides were identical. The ATP-binding sequences found in E. coli protease La were highly conserved. The promoter of the B. brevis lon gene resembled that recognized by the major RNA polymerase of Bacillus subtilis and did not contain sequences homologous to the E. coli heat shock promoters. The B. brevis lon gene was inactivated by insertion of the neomycin resistance gene. A mutant B. brevis carrying the inactivated lon gene showed diminished ability for the degradation of abnormal polypeptides synthesized in the presence of puromycin.     

1H, 15N and 13C resonance assignments and secondary structure determination of the RNA-binding domain of E. coli rho protein

Summary Protein fragments containing the RNA-binding domain of Escherichia coli rho protein have been over-expressed in E. coli. NMR spectra of the fragment containing residues 1–116 of rho protein (Rho116) show that a region of this protein is unfolded in solution. Addition of (dC)₁₀ to this fragment stabilizes the folded form of the protein. The fragment comprising residues 1–130 of rho protein (Rho130) is found to be stably folded, both in the absence and presence of nucleic acid. NMR studies of the complex of Rho 130 with RNA and DNA oligonucleotides indicate that the binding-site size, affinity, and specificity of Rho 130 are similar to those of intact rho protein; therefore, Rho 130 is a suitable model of the RNA-binding domain of rho protein. NMR line widths as well as titration experiments of Rho130 complexed with oligonucleotides of various lengths suggest that Rho130 forms oligomers in the presence of longer oligonucleotides. ¹H, ¹⁵N and ¹³C resonance assignments were facilitated by the utilization of two pulse sequences, CN-NOESY and CCH-TOCSY. The secondary structure of unliganded Rho130 has been determined by NMR techniques, and it is clear that the RNA-binding domain of rho is more structurally similar to the cold shock domain than to the RNA recognition motif.Abbreviations Rho116, Rho130 protein containing the first 116 (130) residues of rho - CSD cold shock domain - RRM RNA recognition motif - RBD RNA-binding domain - IPTG isopropyl -D-thiogalactopyranoside - EDTA ethylenediaminetetraacetic acid - NOE nuclear Overhauser enhancement     

crp genes of Shigella flexneri, Salmonella typhimurium, and Escherichia coli.

The complete nucleotide sequences of the Salmonella typhimurium LT2 and Shigella flexneri 2B crp genes were determined and compared with those of the Escherichia coli K-12 crp gene. The Shigella flexneri gene was almost like the E. coli crp gene, with only four silent base pair changes. The S. typhimurium and E. coli crp genes presented a higher degree of divergence in their nucleotide sequence with 77 changes, but the corresponding amino acid sequences presented only one amino acid difference. The nucleotide sequences of the crp genes diverged to the same extent as in the other genes, trp, ompA, metJ, and araC, which are structural or regulatory genes. An analysis of the amino acid divergence, however, revealed that the catabolite gene activator protein, the crp gene product, is the most conserved protein observed so far. Comparison of codon usage in S. typhimurium and E. coli for all genes sequenced in both organisms showed that their patterns were similar. Comparison of the regulatory regions of the S. typhimurium and E. coli crp genes showed that the most conserved sequences were those known to be essential for the expression of E. coli crp.     

The two homeodomains of the ZmHox2a gene from maize originated as an internal gene duplication and have evolved different target site specificities.

The maize ZmHox2a gene encodes two homeodomains which originated by a 699 bp duplication within an ancestral precursor. The sequences of the two ZmHox2a homeodomains are highly diverged in the N-terminal arm, while residues in the helical part have mostly been conserved. We show here that both ZmHox2a homeodomains are functional DNA-binding motifs but exhibit different target site specificities. CASTing experiments reveal a TCCT motif recognized by HD1 but a GATC tetranucleotide as the recognition sequence of HD2. Mutation of the central nucleotides in both tetranucleotide core motifs abolishes DNA binding. A domain swap experiment indicates that target site specificity is achieved in a combinatorial manner by the contributions of the diverged N-terminal arms together with the slightly different recognition helices. Computer modelling suggests that K47 and H54 in the recognition helices preferentially contact the bases at the 3'-terminus of the tetranucleotide target sequences.     

Analysis of sequence-specific binding of RNA to Hsp70 and its various homologs indicates the involvement of N- and C-terminal interactions.

Members of the 70-kDa family of molecular chaperones assist in a number of molecular interactions that are essential under both normal and stress conditions. These functions require ATP and co-chaperone molecules and are associated with a cyclic transition of intramolecular conformational changes. As a new putative function, we have previously shown that mammalian Hsp/Hsc70 as well as a distant relative, Hsp110, selectively bind certain RNA sequences via their N-terminal ATP-binding domain. To investigate this phenomenon in more detail, here we examined RNA-binding affinity and specificity of various deletion mutants of human Hsp70. We demonstrate, that, although the N-terminal ATPase domain alone is sufficient for RNA binding, its binding affinity is considerably reduced when compared to that of the full-length protein. Additionally, we provide evidence that binding of RNA to a membrane-immobilized protein partner results in complete loss of RNA sequence specificity. Using various Hsp70 homologs, we show distinct RNA-binding properties of these proteins judged by sequence specificity, ribopolymer sensitivity, and northwestern analysis. Finally, we present data disclosing that RNA binding by DnaK, the Escherichia coli homolog, is influenced by the activity of its co-chaperones, DnaJ and GrpE. We conclude that the RNA-binding capability of this class of molecular chaperones is a conserved feature and it is strongly influenced by the structural and conformational properties. Furthermore, the notion that RNA binding of some Hsp70 family members is influenced by co-chaperones suggests an RNA-binding cycle resembling the protein-binding property of the chaperones.     

The Essential Functions of Human Rad51 Are Independent of ATP Hydrolysis

Genetic recombination and the repair of double-strand DNA breaks in Saccharomyces cerevisiae require Rad51, a homologue of the Escherichia coli RecA protein. In vitro, Rad51 binds DNA to form an extended nucleoprotein filament and catalyzes the ATP-dependent exchange of DNA between molecules with homologous sequences. Vertebrate Rad51 is essential for cell proliferation. Using site-directed mutagenesis of highly conserved residues of human Rad51 (hRad51) and gene targeting of the RAD51 locus in chicken DT40 cells, we examined the importance of Rad51's highly conserved ATP-binding domain. Mutant hRad51 incapable of ATP hydrolysis (hRad51K-133R) binds DNA less efficiently than the wild type but catalyzes strand exchange between homologous DNAs. hRad51 does not need to hydrolyze ATP to allow vertebrate cell proliferation, form nuclear foci, or repair radiation-induced DNA damage. However, cells expressing hRad51K-133R show greatly reduced targeted integration frequencies. These findings show that ATP hydrolysis is involved in DNA binding by hRad51 and suggest that the extent of DNA complexed with hRad51 in nucleoprotein influences the efficiency of recombination.     

Covariance analysis of RNA recognition motifs identifies functionally linked amino acids.

The RNA recognition motif (RRM) is one of the most common eukaryotic protein motifs. RRM sequences form a conserved globular structure known as the RNA-binding domain (RBD) or the ribonucleoprotein domain. Many proteins that contain RRM sequences bind RNA in a sequence-specific manner. To investigate the basis for the RNA-binding specificity of RRMs, we subjected 330 aligned RRM sequences to covariance analysis. The analysis revealed a single network of covariant amino acid pairs comprising the buried core of the RBD and a surface patch. Structural studies have implicated a subset of these residues in RNA binding. The covariance linkages identify a larger set of amino acid residues, including some not directly in contact with bound RNA, that may influence RNA-binding specificity.     

Sensory transducers of E. coli are composed of discrete structural and functional domains

The tar and tsr genes of E. coli encode functionally analogous transducer proteins that mediate two distinct classes of chemotactic response. The tap gene lies adjacent to tar, and is thought to encode another transducer protein. We present here the complete nucleotide sequence of the tar-tap region of the E. coli genome, together with a comparative analysis of the sequences of the Tar, Tap, and Tsr proteins. The proteins appear to have a simple transmembrane structure consisting of an extracytoplasmic amino-terminal domain, a membrane-spanning domain, and an intracellular carboxy-terminal domain. The carboxy-terminal domains of three proteins possess highly homologous sequences and contain sites of methylation involved in sensory adaptation, while the amino-terminal sequences are only distantly related to one another, consistent with their serving as chemoreceptor domains that have diverged functionally.     

Structure of a cDNA for the pro alpha 2 chain of human type I procollagen. Comparison with chick cDNA for pro alpha 2(I) identifies structurally conserved features of the protein and the gene

Nucleotide sequences were determined for cloned cDNAs encoding for more than half of the pro alpha 2 chain of type I procollagen from man. Comparisons with previously published data on homologous cDNAs from chick embryos made it possible to examine evolution of the gene in two species which have diverged for 250-300 million years. The amino acid sequence of the alpha-chain domain supported previous indications that there is a strong selective pressure to maintain glycine as every third amino acid and to maintain a prescribed distribution of charged amino acids. However, there is little apparent selective pressure on other amino acids. The amino acid sequence of the C-propeptide domain showed less divergence than the alpha-chain domain. The 5' end or N terminus of the human C-propeptide, however, contained an insert of 12 bases coding for 4 amino acids not found in the chick C-propeptide. About 100 amino acid residues from the N terminus, two residues found in the chick sequence were missing from the human. In the second half of the C-propeptide, there was complete conservation of a 37 amino acid sequence and conservation of 50 out of 51 amino acids in the same region, an observation which suggested that the region serves some special purpose such as directing the association of one pro alpha 2(I) C-propeptide with two pro alpha 1(I) C-propeptides so as to produce the heteropolymeric structure of type I procollagen. In addition, comparison of human and chick DNAs for pro alpha 2(I) revealed three different classes of conservation of nucleotide sequence which have no apparent effect on the structure of the protein: a preference for U on the third base position of codons for glycine, proline, and alanine; a high degree of nucleotide conservation in the 51 amino acid highly conserved region of the C-propeptide; a high degree of nucleotide conservation in the 3'-noncoding region. These three classes of nucleotide conservation may reflect unusual features of collagen genes, such as their high GC content or their highly repetitive coding sequences.     

Arabidopsis proteins containing similarity to the universal stress protein domain of bacteria

We have collected a set of 44 Arabidopsis proteins with similarity to the USPA (universal stress protein A of Escherichia coli) domain of bacteria. The USPA domain is found either in small proteins, or it makes up the N-terminal portion of a larger protein, usually a protein kinase. Phylogenetic tree analysis based upon a multiple sequence alignment of the USPA domains shows that these domains of protein kinases 1.3.1 and 1.3.2 form distinct groups, as do the protein kinases 1.4.1. This indicates that their USPA domain structures have diverged appreciably and suggests that they may subserve distinct cellular functions. Two USPA fold classes have been proposed: one based on Methanococcus jannaschii MJ0577 (1MJH) that binds ATP, and the other based on the Haemophilus influenzae universal stress protein (1JMV), highly similar to E. coli UspA, which does not bind ATP. A set of common residues involved in ATP binding in 1MJH and conserved in similar bacterial sequences is also found in a distinct cluster of Arabidopsis sequences. Threading analysis, which examines aspects of secondary and tertiary structure, confirms this Arabidopsis sequence cluster as highly similar to 1MJH. This structural approach can distinguish between the characteristic fold differences of 1MJH-like and 1JMV-like bacterial proteins and was used to assign the complete set of candidate Arabidopsis proteins to one of these fold classes. It is clear that all the plant sequences have arisen from a 1MJH-like ancestor.     

Molecular cloning of human XPAC gene homologs from chicken, Xenopus laevis and Drosophila melanogaster.

We cloned homologs of the human Xeroderma Pigmentosum Group A complementing (XPAC) gene from chicken, Xenopus laevis and Drosophila melanogaster. A comparison of the amino acid sequences of these homologs with that of the human XPAC protein revealed that in the NH2-terminal domain there are only two conserved regions, one of which is presumed to function as the nuclear localization signal, whereas the COOH-terminal domain is highly conserved, the frequency of identical amino acids in all four XPAC proteins being 50%, and the four cysteine residues predicted to form a zinc-finger motif, and three other cysteine residues are all conserved. These results strongly suggest that the COOH-terminal domain containing a zinc-finger motif plays an important role in the function of these proteins.     

Residues Required for Activity in Escherichia coli o-Succinylbenzoate Synthase (OSBS) Are Not Conserved in All OSBS Enzymes

Understanding how enzyme specificity evolves will provide guiding principles for protein engineering and function prediction. The o-succinylbenzoate synthase (OSBS) family is an excellent model system for elucidating these principles because it has many highly divergent amino acid sequences that are <20% identical, and some members have evolved a second function. The OSBS family belongs to the enolase superfamily, members of which use a set of conserved residues to catalyze a wide variety of reactions. These residues are the only conserved residues in the OSBS family, so they are not sufficient to determine reaction specificity. Some enzymes in the OSBS family catalyze another reaction, N-succinylamino acid racemization (NSAR). NSARs cannot be segregated into a separate family because their sequences are highly similar to those of known OSBSs, and many of them have both OSBS and NSAR activities. To determine how such divergent enzymes can catalyze the same reaction and how NSAR activity evolved, we divided the OSBS family into subfamilies and compared the divergence of their active site residues. Correlating sequence conservation with the effects of mutations in Escherichia coli OSBS identified two nonconserved residues (R159 and G288) at which mutations decrease efficiency ≥200-fold. These residues are not conserved in the subfamily that includes NSAR enzymes. The OSBS/NSAR subfamily binds the substrate in a different orientation, eliminating selective pressure to retain arginine and glycine at these positions. This supports the hypothesis that specificity-determining residues have diverged in the OSBS family and provides insight into the sequence changes required for the evolution of NSAR activity.