|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Department of Biochemistry and 2 Howard Hughes Medical Institute, University of Texas Southwestern Medical Center, Dallas, Texas 75390-9038, USA
Reprint requests to: Jimin Pei, Department of Biochemistry, University of Texas Southwestern Medical Center, 5323 Harry Hines Boulevard, Dallas, TX 75390-9038, USA; e-mail: jpei{at}chop.swmed.edu; fax: (214) 648-9099.
(RECEIVED November 19, 2004; FINAL REVISION January 19, 2005; ACCEPTED January 20, 2005)
| Abstract |
|---|
|
|
|---|
, C-type, G-type lysozymes, and lytic transglycosylases. Since previous biochemical experiments with P5 of phi-6 have indicated that the purified enzyme possesses endopeptidase activity and not glycosidase activity, our results point to the possibility of a newly evolved molecular function and call for further experimental characterization of this unusual P5 protein. Keywords: structure prediction; bacteriophage phi-6; peptidase; lytic transglycosylase; lysozyme; comparative phage genomics
Article published online ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.041250005.
| Introduction |
|---|
|
|
|---|
In MEROPS, the U40 family contains a single sequence, the P5 protein from bacteriophage phi-6. Bacteriophage phi-6 is a double-stranded RNA virus that infects Pseudomonae (Semancik et al. 1973). The P5 protein has been shown to be a lytic enzyme that cleaves the bacterial peptidoglycan (murein) layer to facilitate membrane fusion during infection and is also responsible for cell lysis in late infection (Mindich and Lehman 1979; Bamford and Palva 1980). Caldentey and Bamford (1992) purified and characterized P5 for substrate specificity. The protein was shown to be active against cell walls of various Gram-negative species. The samples of cell walls treated with P5 reacted positively with 2,4-dinitrophenol, indicating that amino groups were liberated by the enzyme. In contrast, neither muramicitol nor glucosaminitol was detected. Based on these observations, it was suggested that P5 cleaves a peptide bond and is not a glycosidase (Caldentey and Bamford 1992). The investigators of the MEROPS database have thus classified the P5 protein as a peptidase and assigned a function of "murein endopeptidase." However, the catalytic mechanism and active site residues have not been identified for this protein.
CDD (Marchler-Bauer et al. 2002), PFAM (Bateman et al. 2004), and SMART (Letunic et al. 2004) database searches using the P5 protein sequence did not yield significant hits to known domains. A PSI-BLAST (Altschul et al. 1997) search with the P5 protein (gene identification [gi] no. 20330562, 220 residues) did not detect any other sequences with significant e-values (<0.02). The 3D-JURY fold recognition META-server (Ginalski et al. 2003) also did not yield any significant hits to existing structures. The P5 protein is apparently a singleton sequence (Pei and Grishin 2002; Siew et al. 2004) without close homologs in the current protein database (the nr database, October 2004, 2,082,196 sequences; 699,810,385 total letters). However, a sequence (gene identification no. 14422162, local alignment from 10185 out of a length of 247 residues) from bacteriophage phi-12 was found as the best BLAST hit of P5 with an e-value of 0.045.
Bacteriophage phi-12 (Gottlieb et al. 2002a,b) belongs to the same genus (Cystovirus) as bacteriophage phi-6, and they are considered to be closely related evolutionarily. Both phages have three segments of chromosomal double-stranded RNA. This weak hit from bacteriophage phi-12, annotated as a muramidase, has the same chromosomal location (Gottlieb et al. 2002b) as P5 in bacteriophage phi-6, and therefore is likely to be a homolog of P5 protein. A PSI-BLAST search using this muramidase as a query finds the P5 protein with an e-value of 0.017 in the first round as the best hit, and further iterations find significant hits to lytic transglycosylases (Koraimann 2003) that belong to the lysozyme superfamily. For example, the lytic transglycosylase from Bartonella henselae was found with an e-value of 3e - 11 in the fourth iteration. To further verify the homology relationship between the P5 proteins and lysozymes, we started transitive PSI-BLAST searches from the lysozyme domain of the lytic transglycosylase with known structure (Protein Data Bank [PDB] ID 1SLY [PDB] , residues 451618; Thunnissen et al. 1994) (see Materials and Methods). Indeed, P5 proteins from both phages phi-6 and phi-12 were found during the course of these extensive searches. For instance, PSI-BLAST with the query gi|48869064 (a soluble lytic murein transglycosylase) found the P5 protein from bacteriophage phi-6 with an e-value of 0.01 in the third iteration.
Additional evidence for the existence of viral homologs is provided by the P5 protein of bacteriophage phi-13, which was found in transitive PSI-BLAST iterations and contained a domain closely related to lytic transglycosylases, as indicated by CDD or PFAM searches. Bacteriophage phi-13 (Qiao et al. 2000) also belongs to the genus of Cystovirus with three segments of double-stranded RNA and has the same chromosomal location of the P5 protein. A complete survey of existing Cystovirus proteins at the National Center for Biotechnology Information (NCBI) Web site did not reveal other proteins with a lytic transglycosylase domain. Complete genome sequence is available for only one other Cystovirus member: bacteriophage phi-8 (Hoogstraten et al. 2000). However, no phi-8 proteins were found to be close homologs of the P5 proteins from phi-6, phi-12, and phi-13. Unlike bacteriophages phi-12 and phi-13, the bacteriophage phi-8 is very distantly related to phi-6, displaying marked differences in gene structure and sequence (Hoogstraten et al. 2000).
The lytic transglycosylase (PFAM: SLT; COG0741) has the same fold as the well-studied hen egg-white lysozyme (HEWL) and T4 lysozyme (Strynadka and James 1996). In the SCOP database (version 1.65) (Murzin et al. 1995), they all belong to the "Lysozyme-like" fold, which consists of one superfamily and seven families. Extensive PSI-BLAST iterations for the lytic transglycosylases (bacterial muramidase, catalytic domain in SCOP) found ~2000 homologs and linked the eukaryotic families of G-type (goose-type) lysozymes and C-type (chicken-type) lysozymes. Other lysozyme families in the SCOP database, such as the families of "Phage T4 lysozymes" and "Chitosanase," are not found using the lytic transglycosylases as queries. A multiple sequence alignment was constructed by using program PCMA (Pei et al. 2003) for the representative sequences of the lytic transglycosylases, G-type lysozymes, C-type lysozymes, and the P5 proteins from bacteriophages phi-6, phi-12, and phi-13. The alignment was manually adjusted guided by the structural superposition of known structures. To reflect similarities among the sequences, a distance diagram was made by using the Euclidean distance mapping method (Fig. 1
; Grishin and Grishin 2002). Close homologs of lytic transglycosylases, G-type lysozymes, and C-type lysozymes form three well-separated clusters. The P5 protein from bacteriophage phi-13 belongs to the cluster of lytic transglycosylases, although it is relatively far from the cluster center. The P5 proteins from bacteriophages phi-6 and phi-12 are far from the rest of the sequences, suggesting an elevated substitution rate and rapid evolution of these bacteriophage sequences. Consistent with the difficulty of homology inference, the phi-6 P5 protein is the most distant from the rest of the transglycosylase sequences.
|
-helix that contains the conserved catalytic residue glutamate. This
-helix is mostly buried and consists of mainly hydrophobic residues. The catalytic glutamate is situated at the end of the
-helix. The position right after the catalytic residue is also highly conserved and usually contains a serine residue that makes a critical hydrogen bond to the
-sheet in the structure. In the P5 protein from bacteriophage phi-6, this position is occupied by an asparigine. The second conserved motif has a sequence signature of GXXQ, where X is often a hydrophobic residue. This motif corresponds to the second turn of a three-stranded
-sheet in the active site. The glycine is highly conserved because of its unique backbone conformation. The glutamine makes critical interactions with a few backbone polar atoms. The third conserved motif has a sequence signature of
p, where
is an aromatic residue and p stands for a polar residue. The motif is located at the end of an
-helix that points toward the active site in lytic transglycosylases and G-type lysozymes (
2 in Fig. 2B
2 helix is almost deteriorated and the polar residue is not conserved.
|
It is well known that many viruses evolve rapidly (Iyer et al. 2001) and new molecular functions could appear as a result of significant sequence divergence (Todd et al. 2001, 2002; Bartlett et al. 2002). The new peptidase active site could have emerged at a location different from the transglycosylase active site, for instance, in the peptide-binding groove. The experimental results suggested that the peptidase activity was directed at the peptide bond between diaminopimelic acid and D-alanine, which is quite far from the glycosidic bond targeted by lysozymes. This scenario has been suggested for the bifunctional invertebrate lysozyme from T. japonica, for which experimental studies have indicated that the isopeptidase activity and chitinase activity are housed at different active sites (Takeshita et al. 2003). Alternatively, the experimental results based on the phi-6 P5 protein purified to apparent homogeneity from disrupted viral particles (Caldentey and Bamford 1992) might not reflect its in vivo activity; or there could be a slight amount of contamination by other endopeptidases in the purified P5 protein. In addition, some membrane-bound proteases that maintain bacterial cell wall might have not been removed during substrate preparation in Caldentey and Bamfords experiments. Preservation of conserved motifs in the P5 protein strongly suggests that it maintains transglycosylase activity. Further experimental studies are required to clarify the in vivo enzymatic activity of the P5 protein from bacteriophage phi-6 and other similar bacteriophages. Our analysis and predictions structurally annotate the U40 peptidase family and offer testable hypotheses about its potential active site residues.
| Materials and methods |
|---|
|
|
|---|
A multiple sequence alignment was constructed by using the PCMA program (Pei et al. 2003) for representative sequences of C-type lysozymes, G-type lysozymes, bacterial lytic transglycosylases, and the P5 proteins. Manual adjustment of the multiple sequence alignment was made with guidance from available structures. Euclidean distance plot was made from the alignment by using the EESG program (Grishin and Grishin 2002).
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Bachali, S., Jager, M., Hassanin, A., Schoentgen, F., Jolles, P., Fiala-Medioni, A., and Deutsch, J.S. 2002. Phylogenetic analysis of invertebrate lysozymes and the evolution of lysozyme function. J. Mol. Evol. 54: 652664.[CrossRef][Medline]
Bamford, D.H. and Palva, E.T. 1980. Structure of the lipid-containing bacteriophage phi 6: Disruption by Triton X-100 treatment. Biochim. Biophys. Acta 601: 245259.[Medline]
Bartlett, G.J., Porter, C.T., Borkakoti, N., and Thornton, J.M. 2002. Analysis of catalytic residues in enzyme active sites. J. Mol. Biol. 324: 105121.[CrossRef][Medline]
Baskova, I.P., Zavalova, L.L., Basanova, A.V., and Sass, A.V. 2001. Separation of monomerizing and lysozyme activities of destabilase from medicinal leech salivary gland secretion. Biochemistry 66: 13681373.[CrossRef][Medline]
Bateman, A., Coin, L., Durbin, R., Finn, R.D., Hollich, V., Griffiths-Jones, S., Khanna, A., Marshall, M., Moxon, S., Sonnhammer, E.L., et al. 2004. The Pfam protein families database. Nucleic Acids Res. 32: D138D141.
Berman, H.M., Battistuz, T., Bhat, T.N., Bluhm, W.F., Bourne, P.E., Burkhardt, K., Feng, Z., Gilliland, G.L., Iype, L., Jain, S., et al. 2002. The Protein Data Bank. Acta Crystallogr. D Biol. Crystallogr. 58: 899907.[CrossRef][Medline]
Caldentey, J. and Bamford, D.H. 1992. The lytic enzyme of the Pseudomonas phage phi 6: Purification and biochemical characterization. Biochim. Biophys. Acta 1159: 4450.[CrossRef][Medline]
Cheng, H., Shen, N., Pei, J., and Grishin, N.V. 2004. Double-stranded DNA bacteriophage prohead protease is homologous to herpesvirus protease. Protein Sci. 13: 22602269.
Esnouf, R.M. 1997. An extensively modified version of MolScript that includes greatly enhanced coloring capabilities. J. Mol. Graph. Model. 15: 133138.
Ginalski, K., Elofsson, A., Fischer, D., and Rychlewski, L. 2003. 3D-Jury: A simple approach to improve protein structure predictions. Bioinformatics 19: 10151018.
Ginalski, K., Kinch, L., Rychlewski, L., and Grishin, N.V. 2004a. BTLCP proteins: A novel family of bacterial transglutaminase-like cysteine proteinases. Trends Biochem. Sci. 29: 392395.[CrossRef][Medline]
. 2004b. Raptor protein contains a caspase-like domain. Trends Biochem. Sci. 29: 522524.[CrossRef][Medline]
Gottlieb, P., Potgieter, C., Wei, H., and Toporovsky, I. 2002a. Characterization of phi12, a bacteriophage related to phi6: Nucleotide sequence of the large double-stranded RNA. Virology 295: 266271.[CrossRef][Medline]
Gottlieb, P., Wei, H., Potgieter, C., and Toporovsky, I. 2002b. Characterization of phi 12, a bacteriophage related to phi 6: Nucleotide sequence of the small and middle double-stranded RNA. Virology 293: 118124.[CrossRef][Medline]
Grishin, V.N. and Grishin, N.V. 2002. Euclidian space and grouping of biological objects. Bioinformatics 18: 15231534.
Holtje, J.V., Mirelman, D., Sharon, N., and Schwarz, U. 1975. Novel type of murein transglycosylase in Escherichia coli. J. Bacteriol. 124: 10671076.
Hoogstraten, D., Qiao, X., Sun, Y., Hu, A., Onodera, S., and Mindich, L. 2000. Characterization of phi8, a bacteriophage containing three double-stranded RNA genomic segments and distantly related to Phi6. Virology 272: 218224.[CrossRef][Medline]
Iyer, L.M., Aravind, L., and Koonin, E.V. 2001. Common origin of four diverse families of large eukaryotic DNA viruses. J. Virol. 75: 1172011734.
Koraimann, G. 2003. Lytic transglycosylases in macromolecular transport systems of Gram-negative bacteria. Cell. Mol. Life Sci. 60: 23712388.[CrossRef][Medline]
Letunic, I., Copley, R.R., Schmidt, S., Ciccarelli, F.D., Doerks, T., Schultz, J., Ponting, C.P., and Bork, P. 2004. SMART 4.0: Towards genomic data integration. Nucleic Acids Res. 32: D142D144.
Marchler-Bauer, A., Panchenko, A.R., Shoemaker, B.A., Thiessen, P.A., Geer, L.Y., and Bryant, S.H. 2002. CDD: A database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res. 30: 281283.
Mindich, L. and Lehman, J. 1979. Cell wall lysin as a component of the bacteriophage phi 6 virion. J. Virol. 30: 489496.
Murzin, A.G., Brenner, S.E., Hubbard, T., and Chothia, C. 1995. SCOP: A structural classification of proteins database for the investigation of sequences and structures. J. Mol. Biol. 247: 536540.[CrossRef][Medline]
Mushegian, A.R., Fullner, K.J., Koonin, E.V., and Nester, E.W. 1996. A family of lysozyme-like virulence factors in bacterial pathogens of plants and animals. Proc. Natl. Acad. Sci. 93: 73217326.
Pei, J. and Grishin, N.V. 2001. Type II CAAX prenyl endopeptidases belong to a novel superfamily of putative membrane-bound metalloproteases. Trends Biochem. Sci. 26: 275277.[CrossRef][Medline]
. 2002. Breaking the singleton of germination protease. Protein Sci. 11: 691697.
. 2003. Peptidase family U34 belongs to the superfamily of N-terminal nucleophile hydrolases. Protein Sci. 12: 11311135.
Pei, J., Sadreyev, R., and Grishin, N.V. 2003. PCMA: Fast and accurate multiple sequence alignment based on profile consistency. Bioinformatics 19: 427428.
Qiao, X., Qiao, J., Onodera, S., and Mindich, L. 2000. Characterization of phi 13, a bacteriophage related to phi 6 and containing three dsRNA genomic segments. Virology 275: 218224.[CrossRef][Medline]
Rawlings, N.D., Tolle, D.P., and Barrett, A.J. 2004. MEROPS: The peptidase database. Nucleic Acids Res. 32: D160D164.
Semancik, J.S., Vidaver, A.K., and Van Etten, J.L. 1973. Characterization of segmented double-helical RNA from bacteriophage phi6. J. Mol. Biol. 78: 617625.[CrossRef][Medline]
Siew, N., Azaria, Y., and Fischer, D. 2004. The ORFanage: An ORFan database. Nucleic Acids Res. 32: D281D283.
Strynadka, N.C. and James, M.N. 1996. Lysozyme: A model enzyme in protein crystallography. Exs 75: 185222.[Medline]
Takeshita, K., Hashimoto, Y., Ueda, T., and Imoto, T. 2003. A small chimerically bifunctional monomeric protein: Tapes japonica lysozyme. Cell. Mol. Life Sci. 60: 19441951.[CrossRef][Medline]
Thunnissen, A.M., Dijkstra, A.J., Kalk, K.H., Rozeboom, H.J., Engel, H., Keck, W., and Dijkstra, B.W. 1994. Doughnut-shaped structure of a bacterial muramidase revealed by X-ray crystallography. Nature 367: 750753.[CrossRef][Medline]
Thunnissen, A.M., Isaacs, N.W., and Dijkstra, B.W. 1995. The catalytic domain of a bacterial lytic transglycosylase defines a novel class of lysozymes. Proteins 22: 245258.[CrossRef][Medline]
Todd, A.E., Orengo, C.A., and Thornton, J.M. 2001. Evolution of function in protein superfamilies, from a structural perspective. J. Mol. Biol. 307: 11131143.[CrossRef][Medline]
. 2002. Plasticity of enzyme active sites. Trends Biochem. Sci. 27: 419426.[CrossRef][Medline]
Walker, D.R. and Koonin, E.V. 1997. SEALS: A system for easy analysis of lots of sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 5: 333339.[Medline]
Zavalova, L., Lukyanov, S., Baskova, I., Snezhkov, E., Akopov, S., Berezhnoy, S., Bogdanova, E., Barsova, E., and Sverdlov, E.D. 1996. Genes from the medicinal leech (Hirudo medicinalis) coding for unusual enzymes that specifically cleave endo-epsilon (
-Glu)-Lys isopeptide bonds and help to dissolve blood clots. Mol. Gen. Genet. 253: 2025.[CrossRef][Medline]
Zavalova, L.L., Baskova, I.P., Lukyanov, S.A., Sass, A.V., Snezhkov, E.V., Akopov, S.B., Artamonova, II, Archipova, V.S., Nesmeyanov, V.A., Kozlov, D.G., et al. 2000. Destabilase from the medicinal leech is a representative of a novel family of lysozymes. Biochim. Biophys. Acta 1478: 6977.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
J. Soding, M. Remmert, A. Biegert, and A. N. Lupas HHsenser: exhaustive transitive profile search using HMM-HMM comparison. Nucleic Acids Res., July 1, 2006; 34(Web Server issue): W374 - W378. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Pei and N. V. Grishin COG3926 and COG5526: A tale of two new lysozyme-like protein families Protein Sci., October 1, 2005; 14(10): 2574 - 2581. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |