Conserved sequence motifs among bacterial,eukaryotic, and archaeal phosphatases that define a new phosphohydrolase superfamily.

Members of a new molecular family of bacterial nonspecific acid phosphatases (NSAPs), indicated as class C, were found to share significant sequence similarities to bacterial class B NSAPs and to some plant acid phosphatases, representing the first example of a family of bacterial NSAPs that has a relatively close eukaryotic counterpart. Despite the lack of an overall similarity, conserved sequence motifs were also identified among the above enzyme families (class B and class C bacterial NSAPs, and related plant phosphatases) and several other families of phosphohydrolases, including bacterial phosphoglycolate phosphatases, histidinol-phosphatase domains of the bacterial bifunctional enzymes imidazole-glycerolphosphate dehydratases, and bacterial, eukaryotic, and archaeal phosphoserine phosphatases and threalose-6-phosphatases. These conserved motifs are clustered within two domains, separated by a variable spacer region, according to the pattern [FILMAVT]-D-[ILFRMVY]-D-[GSNDE]-[TV]-[ILVAM]-[AT S VILMC]-X-¿YFWHKR)-X-¿YFWHNQ¿-X( 102,191)-¿KRHNQ¿-G-D-¿FYWHILVMC¿-¿QNH¿-¿FWYGP¿-D -¿PSNQYW¿. The dephosphorylating activity common to all these proteins supports the definition of this phosphatase motif and the inclusion of these enzymes into a superfamily of phosphohydrolases that we propose to indicate as "DDDD" after the presence of the four invariant aspartate residues. Database searches retrieved various hypothetical proteins of unknown function containing this or similar motifs, for which a phosphohydrolase activity could be hypothesized.     

Constant relative rate of protein evolution and detection of functional diversification among bacterial, archaeal and eukaryotic proteins

Background  

Detection of changes in a protein's evolutionary rate may reveal cases of change in that protein's function. We developed and implemented a simple relative rates test in an attempt to assess the rate constancy of protein evolution and to detect cases of functional diversification between orthologous proteins. The test was performed on clusters of orthologous protein sequences from complete bacterial genomes (Chlamydia trachomatis, C. muridarum and Chlamydophila pneumoniae), complete archaeal genomes (Pyrococcus horikoshii, P. abyssi and P. furiosus) and partially sequenced mammalian genomes (human, mouse and rat).     

B. subtilis ykuD protein at 2.0 A resolution: insights into the structure and function of a novel, ubiquitous family of bacterial enzymes

The crystal structure of the product of the Bacillus subtilis ykuD gene was solved by the multiwavelength anomalous dispersion (MAD) method and refined using data to 2.0 A resolution. The ykuD protein is a representative of a distinctly prokaryotic and ubiquitous family found among both pathogenic and nonpathogenic Gram-positive and Gram-negative bacteria. The deduced amino acid sequence reveals the presence of an N-terminal LysM domain, which occurs among enzymes involved in cell wall metabolism, and a novel, putative catalytic domain with a highly conserved His/Cys-containing motif of hitherto unknown structure. As the wild-type protein did not crystallize, a double mutant was designed (Lys117Ala/Gln118Ala) to reduce excess surface conformational entropy. As expected, the structure of the LysM domain is similar to the NMR structure reported for an analogous domain from Escherichia coli murein transglycosylase MltD. The molecular model also shows that the 112-residue-long C-terminal domain has a novel tertiary fold consisting of a beta-sandwich with two mixed sheets, one containing five strands and the other, six strands. The two beta-sheets form a cradle capped by an alpha-helix. This domain contains a putative catalytic site with a tetrad of invariant His123, Gly124, Cys139, and Arg141. The stereochemistry of this active site shows similarities to peptidotransferases and sortases, and suggests that the enzymes of the ykuD family may play an important role in cell wall biology.     

Prediction of EF-hand calcium-binding proteins and analysis of bacterial EF-hand proteins

The EF-hand protein with a helix-loop-helix Ca(2+) binding motif constitutes one of the largest protein families and is involved in numerous biological processes. To facilitate the understanding of the role of Ca(2+) in biological systems using genomic information, we report, herein, our improvement on the pattern search method for the identification of EF-hand and EF-like Ca(2+)-binding proteins. The canonical EF-hand patterns are modified to cater to different flanking structural elements. In addition, on the basis of the conserved sequence of both the N- and C-terminal EF-hands within S100 and S100-like proteins, a new signature profile has been established to allow for the identification of pseudo EF-hand and S100 proteins from genomic information. The new patterns have a positive predictive value of 99% and a sensitivity of 96% for pseudo EF-hands. Furthermore, using the developed patterns, we have identified zero pseudo EF-hand motif and 467 canonical EF-hand Ca(2+) binding motifs with diverse cellular functions in the bacteria genome. The prediction results imply that pseudo EF-hand motifs are phylogenetically younger than canonical EF-hand motifs. Our prediction of Ca(2+) binding motifs provides not only an insight into the role of Ca(2+) and Ca(2+)-binding proteins in bacterial systems, but also a way to explore and define the role of Ca(2+) in other biological systems (calciomics).     

A heat shock protein localized to chloroplasts is a member of a eukaryotic superfamily of heat shock proteins.

We have isolated cDNA clones from soybean and pea that specify nuclear-encoded heat shock proteins (HSPs) which localize to chloroplasts. The mRNAs for these HSPs are undetectable at control temperatures, but increase approximately 150-fold during a 2-h heat shock. Hybridization-selection followed by in vitro translation demonstrates that these HSPs are synthesized as precursor proteins which are processed by the removal of 5-6.5 kd during import into isolated chloroplasts. The nucleotide sequence of the cDNAs shows the derived amino acid sequences of the mature pea and soybean proteins are 79% identical. While the predicted transit peptide encoded by the pea cDNA has some characteristics typical of transit sequences, including high Ser content, multiple basic residues and no acidic residues, it lacks two domains proposed to be important for import and maturation of other chloroplast proteins. The carboxy-terminal region of the chloroplast HSP has significant homology to cytoplasmic HSPs from soybean and other eukaryotes. We hypothesize that the chloroplast HSP shares a common structural and functional domain with low mol. wt HSPs which localize to other parts of the cell, and may have evolved from a nuclear gene.     

Linkers of secondary structures in proteins.

Linkers that connect repeating secondary structures fall into conformational classes based on distance and main-chain torsion clustering. A data set of 300 unique protein chains with low pairwise sequence identity was clustered into only a few groups representing the preferred motifs. The linkers of two to eight residues for the nonredundant data set are designated H-Ln-H, H-Ln-E, E-Ln-H, E-Ln-E, where n is the length, H stands for alpha-helices, and E for beta-strands. Most of the clusters identified here corroborate earlier findings. However, 19 new clusters are identified in this paper, with many of them having seven and eight residue linkers. In our first analysis, the secondary structures flanking the linkers are both interacting and noninteracting and there is no precise angle of orientation between them. A second analysis was performed on a set of proteins with restricted orientations for the flanking elements, namely, mainly alpha class of proteins with orthogonal architecture. Two definite clusters are identified, one corresponding to linkers of orthogonal helices and the other to linkers of antiparallel helices. Loops forming binding sites or involved in catalytic activity are important determinants of the function of proteins. Although the structural conservation of the residues around the catalytic triad of serine proteases has been studied widely, there has not been a systematic analysis of the conformation of the loops that contain them. Residues of the catalytic triad reside in the linkers of beta-strands, with varying lengths of more than eight residues. Here, we analyze the structural conservation of such linkers by superposition, and observe a conserved structural feature of the linkers incorporating each of the three residues of the catalytic triad.     

Homologs of eukaryotic Ras superfamily proteins in prokaryotes and their novel phylogenetic correlation with their eukaryotic analogs

Ras superfamily proteins are key regulators in a wide variety of cellular processes. Previously, they were considered to be specific to eukaryotes, and MglA, a group of obviously different prokaryotic proteins, were recognized as their only prokaryotic analogs or even ancestors. Here, taking advantage of quite a current accumulation of prokaryotic genomic databases, we have investigated the existence and taxonomic distribution of Ras superfamily protein homologs in a much wider prokaryotic range, and analyzed their phylogenetic correlation with their eukaryotic analogs. Thirteen unambiguous prokaryotic homologs, which possess the GDP/GTP-binding domain with all the five characteristic motifs of their eukaryotic analogs, were identified in 12 eubacteria and one archaebacterium, respectively. In some other archaebacteria, including four methanogenic archaebacteria and three Thermoplasmales, homologs were also found, but with the GDP/GTP-binding domains not containing all the five characteristic motifs. Many more MglA orthologs were identified than in previous studies mainly in delta-proteobacteria, and all were shown to have common unique features distinct from the Ras superfamily proteins. Our phylogenetic analysis indicated eukaryotic Rab, Ran, Ras, and Rho families have the closest phylogenetic correlation with the 13 unambiguous prokaryotic homologs, whereas the other three eukaryotic protein families (SRbeta, Sar1, and Arf) branch separately from them, but have a relatively close relationship with the methanogenic archaebacterial homologs and MglA. Although homologs were identified in a relative minority of prokaryotes with genomic databases, their presence in a relatively wide variety of lineages, their unique sequence characters distinct from those of eukaryotic analogs, and the topology of our phylogenetic tree altogether do not support their origin from eukaryotes as a result of lateral gene transfer. Therefore, we argue that Ras superfamily proteins might have already emerged at least in some prokaryotic lineages, and that the seven eukaryotic protein families of the Ras superfamily may have two independent prokaryotic origins, probably reflecting the 'fusion' evolutionary history of the eukaryotic cell.     

The geometry of domain combination in proteins.

Most proteins in genomes are the result of the recombination of two or more domains. It has been found that if proteins are formed by a combination of domains from superfamilies A and B, then the domains may occur in the sequential order AB or BA but only in about 2% of cases do they occur in both sequential orders. The classical Rossmann domains of known structure are combined with catalytic domains from seven different superfamilies. In addition, there are eight cases where structures with both AB and BA domain combinations are known. For these two sets of structures, we analysed: (i) the relative orientation of the domains; (ii) the type of domain connection; (iii) the structure of the interdomain links; and (iv) domain function. The results of this analysis indicate that in most cases domain order is conserved because recombination of the domains has only occurred once during the course of evolution. Functional reasons become important when the domain connections are short. In seven out of the eight known cases where domains are combined in the AB and BA sequential orders they have different geometrical relationships that give them different functional properties.     

An archaeal protein homologous to mammalian SRP54 and bacterial Ffh recognizes a highly conserved region of SRP RNA.

The gene encoding the 54 kDa protein of signal recognition particle (SRP54) in the hyperthermophilic archaeon Pyrococcus furiosus has been cloned and sequenced. Recombinant P. furiosus SRP54 (pf-SRP54) and the N-terminal G-domain and C-terminal M-domain (pf-SRP54M) of pf-SRP54 with an amino-terminal addition of six histidine residues were expressed in Escherichia coli and subjected to binding experiments for SRP RNA, non-conserved 213-nucleotide RNA (helices 1, 2, 3, 4 and 5) and conserved 107-nucleotide RNA (helices 6 and 8) from SRP RNA. The RNA binding properties of the purified protein were determined by filter binding assays. The histidine-tagged pf-SRP54M bound specifically to the conserved 107-nucleotide RNA in the absence of pf-SRP19, unlike the eukaryotic homologue, with an apparent binding constant (K) of 18 nM.     

Gneg-mPLoc: A top-down strategy to enhance the quality of predicting subcellular localization of Gram-negative bacterial proteins

By incorporating the information of gene ontology, functional domain, and sequential evolution, a new predictor called Gneg-mPLoc was developed. It can be used to identify Gram-negative bacterial proteins among the following eight locations: (1) cytoplasm, (2) extracellular, (3) fimbrium, (4) flagellum, (5) inner membrane, (6) nucleoid, (7) outer membrane, and (8) periplasm. It can also be used to deal with the case when a query protein may simultaneously exist in more than one location. Compared with the original predictor called Gneg-PLoc, the new predictor is much more powerful and flexible. For a newly constructed stringent benchmark dataset in which none of proteins included has ≥25% pairwise sequence identity to any other in a same subset (location), the overall jackknife success rate achieved by Gneg-mPLoc was 85.5%, which was more than 14% higher than the corresponding rate by the Gneg-PLoc. As a user friendly web-server, Gneg-mPLoc is freely accessible at http://www.csbio.sjtu.edu.cn/bioinf/Gneg-multi/.     

A structure-based proposal for the catalytic mechanism of the bacterial acid phosphatase AphA belonging to the DDDD superfamily of phosphohydrolases

The Escherichia coli gene aphA codes for a periplasmic acid phosphatase called AphA, belonging to class B bacterial phosphatases, which is part of the DDDD superfamily of phosphohydrolases. After our first report about its crystal structure, we have started a series of crystallographic studies aimed at understanding of the catalytic mechanism of the enzyme. Here, we report three crystal structures of the AphA enzyme in complex with the hydrolysis products of nucleoside monophosphate substrates and a fourth with a proposed intermediate analogue that appears to be covalently bound to the enzyme. Comparison with the native enzyme structure and with the available X-ray structures of different phosphatases provides clues about the enzyme chemistry and allows us to propose a catalytic mechanism for AphA, and to discuss it with respect to the mechanism of other bacterial and human phosphatases.     

Plant progesterone 5beta-reductase is not homologous to the animal enzyme. Molecular evolutionary characterization of P5betaR from Digitalis purpurea

Plants of the genus Digitalis produce cardiac glycosides, i.e. digoxin, which are widely used for congestive heart failure. Progesterone 5beta-reductase (P5betaR) is a key enzyme in the biosynthesis of these natural products. Here, we have carried out the purification and partial amino acid sequencing of the native P5betaR from foxglove (Digitalis purpurea), and isolated a cDNA encoding this enzyme. Similarly to other steroid 5beta-reductases, the recombinant P5betaR catalyzes the stereospecific reduction of the Delta(4)-double bond of several steroids with a 3-oxo,Delta(4,5) structure. The gene encoding P5betaR is expressed in all plant organs, and maximally transcribed in leaves and mature flowers. P5betaR belongs to the short-chain dehydrogenase/reductase (SDR) superfamily, bearing no structural homology to its mammalian counterpart, which is a member of the aldo-keto reductase (AKR) superfamily. A similar situation occurs with 3beta-hydroxy-Delta(5)-steroid dehydrogenase (3betaHSD), the gene immediately preceding P5betaR in the cardenolide pathway, which suggests that the entire route has evolved independently in animals and plants. P5betaR is retained only in plants, where it is ubiquitous, and a few distantly related bacterial lineages after its diversification from the last universal common ancestor. Evolutionary conserved changes in its putative active site suggest that plant P5betaR is a member of a novel subfamily of extended SDRs, or a new SDR family.     

Isolation,screening, and optimization of bacterial strains for novel transglutaminase production

Transglutaminases are a class of transferases known to form isopeptide bond between glutamine and lysine residues in a protein molecule. Increasing demand for transglutaminase in food and other industries and its low productivity have compelled researchers to isolate and screen micro-organisms with potential to produce it. In the present investigation around 200 isolates were screened for extracellular secretion of microbial transglutaminase (MTGase). Isolate B4 showed enzyme activity of 1.71?±?0.2?U/mL followed by isolate C2 which showed 1.61?±?0.17?U/mL activity, comparable with the activity of industrially used microbial strains. Biochemical analysis along with 16S r-RNA sequencing revealed these isolates (B4 and C2) to be Bacillus nakamurai and a variant of Bacillus subtilis, respectively. Amongst the various production media screened, a medium containing starch and peptone was found best for MTGase production. Correlation between growth, enzyme production, and sugar utilization was also studied and maximum enzyme production was obtained after 48 to 60?hr. Highest MTGase titer (3.95?±?0.03?U/mL for B4 and 2.65?±?0.17?U/mL for C2) was obtained by optimization of parameters. The enzyme was characterized for temperature and pH optima, pH and thermal stability, and effect of metal ions, suggesting its potential use in future applications.     

Merlin, a new superfamily of DNA transposons identified in diverse animal genomes and related to bacterial IS1016 insertion sequences

Several new families of DNA transposons were identified by computer-assisted searches in a wide range of animal species that includes nematodes, flat worms, mosquitoes, sea squirt, zebrafish, and humans. Many of these elements have coding capacity for transposases, which are related to each other and to those encoded by the IS1016 group of bacterial insertion sequences. Although these transposases display a motif similar to the DDE motif found in many transposases and integrases, they cannot be directly allied to any of the previously described eukaryotic transposases. Other common features of the new eukaryotic and bacterial transposons include similarities in their terminal inverted repeats and 8-bp or 9-bp target-site duplications. Together, these data indicate that these elements belong to a new superfamily of DNA transposons, called Merlin/IS1016, which is common in many eubacterial and animal genomes. We also present evidence that these transposons have been recently active in several animal species. This evidence is particularly strong in the parasitic blood fluke Schistosoma mansoni, in which Merlin is also the first described DNA transposon family.     

Cloning, expression, and crystallization of a hyperthermophilic protein that is homologous to the eukaryotic translation initiation factor, eIF5A.

A gene coding for a protein homologous to a translation initiation factor of eukaryotes, eIF5A, was cloned from Methanococcus jannaschii, a hyperthermophile with an optimum growth temperature of 85 degrees C. The protein was overexpressed, purified and crystallized. The crystals were obtained by vapor diffusion method with 8% PEG 4000 as precipitant and belong to space group P4(1)22 with unit cell dimensions a = b = 45.52 A and c = 155.59 A. These crystals diffract to at least 2.2 A resolution.     

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.