In silico modeling of the staphylococcal bacteriophage-derived peptidase CHAPK

doi:10.4161/bact.1.4.18245

Bacteriophage. 2011 July 1; 1(4): 198–206.

doi: 10.4161/bact.1.4.18245

PMCID: PMC3448105

In silico modeling of the staphylococcal bacteriophage-derived peptidase CHAP_K

Mark Fenton,¹ Jakki C. Cooney,² R. Paul Ross,³ Roy D. Sleator,¹ Olivia McAuliffe,³ Jim O'Mahony,¹ and Aidan Coffey¹^*

¹Department of Biological Sciences; Cork Institute of Technology; Cork, Ireland

²Department of Life Sciences and Materials and Surface Science Institute; University of Limerick; Limerick, Ireland

³Teagasc; Moorepark Food Research Centre; County Cork, Ireland

^*Correspondence to: Aidan Coffey, Email: aidan.coffey/at/cit.ie

This is an open-access article licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported License. The article may be redistributed, reproduced, and reused for non-commercial purposes, provided the original source is properly cited.

Abstract

The aim of this study was to use comparative modeling to predict the three-dimensional structure of the CHAP_K protein (cysteine, histidine-dependent amidohydrolase/peptidase domain of the LysK endolysin, derived from bacteriophage K). Iterative PSI-BLAST searches against the Protein Data Bank (PDB) and nonredundant (nr) databases were used to populate a multiple alignment for analysis using the T-Coffee Expresso server. A consensus Maximum Parsimony phylogenetic tree with a bootstrap analysis setting of 1,000 replicates was constructed using MEGA4. Structural templates relevant to our target (CHAP_K) were identified, processed in Expresso and used to generate a 3D model in the alignment mode of SWISS-MODEL. These templates were also processed in the I-TASSER web server. A Staphylococcus saprophyticus CHAP domain protein, 2K3A, was identified as the structural template in both servers. The I-TASSER server generated the CHAP_K model with the best bond geometries when analyzed using PROCHECK and the most logical organization of the structure. The predicted 3D model indicates that CHAP_K has a papain-like fold. Circular dichroism spectropolarimetry also indicated that CHAP_K has an αβ fold, which is consistent with the model presented. The putative active site maintained a highly conserved Cys54-His117-Glu134 charge relay and an oxyanion hole residue Asn136. The residue triplet, Cys-His-Glu, is known to be a viable proteolytic triad in which we predict the Cys residue is used in a nucleophilic attack on peptide bonds at a specific site in the pentaglycine cross bridge of staphylococcal cell wall peptidoglycan. Use of comparative modeling has allowed approximation of the 3D structure of CHAP_K giving information on the structure and an insight into the binding and active site of the catalytic domain. This may facilitate its development as an alternative antibacterial agent.

Keywords: bacteriophage, CHAP, endolysin, in silico, peptidase, staphylococcus

Introduction

Staphylococcus aureus is a common pathogen that plays a major role in various human and animal diseases ranging from skin and soft tissue infections to more serious cases of pneumonia, endocarditis, meningitis and osteomyelitis.¹ Treatment of these infections has become increasingly difficult due to the worldwide prevalence of multidrug-resistant strains including methicillin resistant S. aureus (MRSA) which is a frequent cause of serious nosocomial infections.² As a consequence it is critical to develop new and effective antibacterials with the potential of eliminating such infections irrespective of antibiotic sensitivity.

Over the past decade numerous studies have focused on developing recombinant bacteriophage (phage)-encoded cell wall hydrolases, termed endolysins (lysins), as novel antibacterial agents as recently reviewed by Loessner (2005),³ Fischetti (2008)⁴ and Fenton et al. (2010).⁵ When applied exogenously as purified recombinant proteins to Gram-positive bacteria, lysins bring about rapid cell lysis and bacterial death.³^,⁶^,⁷ It is this ability to rapidly lyse pathogenic Gram-positive cells upon direct contact with peptidoglycan, also termed “lysis from without,” that has laid the foundation for their exploitation as novel therapeutics.⁶

The majority of lysins display a modular structure, usually comprising of at least one N-terminal catalytic domain which attacks bacterial cell wall peptidoglycan, combined with a C-terminal cell wall binding domain which directs the lytic domain to its site of action.⁸^,⁹ In the case of staphylococcal lysins the presence of three domains, comprising an N-terminal cysteine, histidine-dependent amidohydrolase/peptidase (CHAP) domain and an amidase-2 domain linked to a C-terminal SH3b cell wall binding domain, is a common feature. This organization has been observed in LysK,¹⁰ LysWMY,¹¹ Φ11 lysin,¹² MV-L,¹³ LysH5,¹⁴ LysGH15,¹⁵ and SAL-1.¹⁶

Studies have shown that of the two lytic domains contained within staphylococcal lysins, the CHAP domain confers the principle enzymatic activity of the protein, whereas the amidase domain contributes minimal detectable activity.¹⁷^-²⁰ Furthermore, a recent study in our laboratory showed that the activity of the single-domain truncated CHAP (later designated CHAP_K), was 2-fold higher than the native three-domain LysK protein.²⁰ Studies have also demonstrated that the CHAP domain of staphylococcal lysins acts as a _D-Ala-Gly endopeptidase, specifically cleaving the peptide bond between _D-alanine and the first glycine in the pentaglycine cross-bridge of staphylococcal cell wall peptidoglycan¹⁷^,¹⁸ (Fig. 1).

Figure 1. CHAP_K cleavage site of staphylococcal cell wall peptidoglycan.

CHAP proteins contain three highly conserved amino acid residues, two of which are an invariant cysteine (Cys) and histidine (His) along with a third, polar residue such as asparagine (Asn), aspartic acid (Asp) or glutamic acid (Glu).²¹ These residues are principally involved in catalytic activity, forming part of the active site of the enzyme as well as being equivalent to the catalytic triad of papain-like thiol proteases.²¹^,²² Site-directed mutagenesis studies on these conserved residues within the CHAP domain of the LysWMY staphylococcal lysin resulted in reduced activity, suggesting that a Cys-His-Asn catalytic triad is crucial for enzymatic function.¹¹

The CHAP domain is a member of the NlpC/P60 family of peptidases and can be found in proteins from bacteria, archaea and eukaryotes of the Trypanosomidae family.²² However, very little structural information is available in relation to the CHAP domains of phage lysins. To date the crystal structure of five lysins has been elucidated, Cpl-1,²³ PlyL,²⁴ PlyPSA,²⁵ PlyB,²⁶ and Ply500,²⁷ none of which are derived from staphylococcal phage and none of which contain a CHAP domain.

In the current study a comparative modeling approach is employed to predict the three dimensional (3D) structure of the 162-amino-acid CHAP domain protein, CHAP_K, derived from the staphylococcal phage K lysin; LysK.¹⁰ The potential of CHAP_K as an alternative antibacterial agent which displays rapid lytic activity against pathogenic staphylococci including MRSA, both in vitro and in vivo, has previously been demonstrated by our group.²⁰^,²⁸

Comparative modeling allows accurate prediction of the 3D structure of a target protein (CHAP_K) based on its similarity to a related template protein whose structure has been experimentally resolved by either X-ray crystallography or Nuclear Magnetic Resonance (NMR) spectroscopy.²⁹ The basis of this approach is that evolutionary related proteins share similar 3D structures if they have a statistically significant sequence similarity (generally > 25% identity).³⁰ The model prediction process consists of four consecutive steps starting with template selection, followed by target-template alignment, model building and finally model evaluation.³¹

Resolution of the 3D structure of CHAP_K by comparative modeling can provide a valuable insight into the structural basis of interaction with its cell wall peptidoglycan substrate and mechanism of catalytic activity. This knowledge may subsequently be used in the design of novel lysins with altered and improved catalytic domains, or the creation of designer lysins with unique capabilities (e.g., enhanced turnover rate, altered specificity) furthering the development of these unique enzymes as alternative antibacterial agents in the fight against multi-drug resistant bacterial infections.

Results and Discussion

Template selection

The first step in the comparative modeling of the target protein, CHAP_K, was to identify potential templates whose structure has been resolved by either X-ray crystallography or NMR. The two main criteria for selecting suitable templates for the target were sequence length (choosing proteins of similar length to CHAP_K) and percentage identity between the target and template sequence of > 25% over the entire length of the sequence. Generally, proteins with > 25% sequence identity can be labeled homologous and share at least 80% of their structure.³⁰^,³² With the aim of collecting the most suitable template sequences to form a strong multiple alignment in the second step of comparative modeling, a PSI-BLAST (Position-Specific Iterated BLAST) search was performed against the Protein Data Bank (PDB) over a number of iterations. This bank is a repository of the 3D structural data of proteins and nucleic acids which has been experimentally obtained by X-ray crystallography or NMR. The most suitable template obtained from the PDB search was the enzyme 2K3A; a 155-amino-acid secreted CHAP protein from Staphylococcus saprophyticus³³ with 28% sequence identity to CHAP_K. This template returned the lowest (E)-value (0.01) and largest % coverage (58%). The BLAST search converged after three iterations and no known structure received a significantly high score, over a range of scoring matrix and gap penalty settings. This suggested that not many CHAP domain structures with similarity to phage-derived CHAP domains have been experimentally resolved to date.

In order to populate the multiple alignment in step two, iterative PSI-BLAST searches were also performed in the non-redundant (nr) database. Although the inclusion of homologous sequences from the nr database does not contribute structural information, it aids in constructing a better quality multiple alignment, which is generally regarded as more reliable than a pair-wise structural alignment alone.³¹

Target-template alignment and model building

In the second step of comparative modeling, the combination of selected template structures from the PDB and sequences from the nr-database, collected in step one, were uploaded to the T-Coffee Expresso server together with the amino acid sequence of the target (CHAP_K). The Expresso suite is a multiple sequence alignment server which accurately aligns sequences that include structural information.³⁴ The resulting multiple alignment of the CHAP_K target and template proteins was used to construct a phylogenetic tree and identify templates which were most closely related to the target. A consensus maximum parsimony phylogenetic tree with bootstrap analysis inferred from 1,000 replicates was constructed using the MEGA4 software suite³⁵ and taken to represent the evolutionary history of the protein sequences analyzed (Fig. 2A). Upon examination of the tree, proteins in close proximity to the target, including the protein 2K3A (with experimentally derived 3D structure) were selected and processed again using the Expresso alignment algorithm (Fig. 2B).

Figure 2. (A) Phylogenetic tree showing evolutionary relationship of 40 CHAP_K homologs (15 from PDB BLAST and 25 from nr BLAST databases). The maximum parsimony tree was obtained using the Close-Neighbor-Interchange algorithm level 3 in MEGA4 (more ...)

For the third step in the process (model building), the above alignment was submitted to the Alignment mode of the web-based homology modeling server, SWISS-MODEL.³⁶ This program facilitates the construction of 3D models of a target protein using experimentally determined structures of related family members.³⁷ Models were constructed using the 2K3A protein and each structurally resolved protein from the clade containing CHAP_K in the phylogenetic tree (Fig. 2A). Inspection of predicted 3D CHAP_K models in relation to structural compatibility between the target and template, sharing of conserved amino acids and the location of these residues in the overall structure, resulted in the identification of the 2K3A protein as the most promising candidate for accurate modeling of CHAP_K.

In model building, accurate alignment of the target sequence with selected template is a crucial step and poor alignments are commonly responsible for inaccurate models.³⁷ In order to examine the quality of the target-template alignment and the models created in SWISS-MODEL Alignment mode, CHAP_K was submitted to an additional web-based server, I-TASSER.³⁸ The results from I-TASSER independently identified the 2K3A protein as the most suitable template for constructing an accurate 3D model of CHAP_K. Upon inspection and comparison of the predicted models from both I-TASSER (Fig. 3A) and SWISS-MODEL (Fig. 3B) servers, it was obvious that I-TASSER had generated a better quality alignment as its predicted model had a fold and catalytic residue position that was consistent with enzymatic function³³ (Fig. 3A).

Figure 3. Predicted 3D structure of CHAP_K by I-TASSER (A) and SwissModel (B) servers. Highly conserved cysteine residue (Cys54) is highlighted in orange and conserved histidine residue (His117) is highlighted in blue for both predicted structures. (more ...)

Model Evaluation

In the fourth step, the quality of the I-TASSER generated CHAP_K model was evaluated using PROCHECK³⁹ and visualized using PyMOL.⁴⁰ The predicted model had two short α-helices packed against a β-sheet saddle, which is a typical feature of NplC/P60 papain-related peptidases²¹ (Fig. 3A). Model refinement was performed using CNS_SOLVE,⁴¹ and the resultant structure had good geometries with 54.3% of residues falling in the most favored region of the Ramachandran plot, and only 2.9% outliers. A comparison between the newly constructed CHAP_K model and the 2K3A CHAP protein template was also assessed. Superposition of CHAP_K on the template structure used LSQMAN.⁴² The calculated root mean square distance (r.m.s.d) value of 1.9 Ångstroms (Å) on 143 C-α atoms indicated the expected degree of agreement between the two structures. Inspection of the superimposed structures in PyMOL showed good agreement between CHAP_K and the S. saprophyticus CHAP domain throughout the structures (Fig. 4).

Figure 4. Structural comparison by superposition of the experimentally resolved 2K3A protein (cyan) and predicted CHAP_K structure (magenta). Highly conserved active site residues are identified and labeled for both proteins. CHAP_K Cys54 (orange) (more ...)

Further analysis of secondary structure in the CHAP_K model was performed using the software package Digital Shape Sampling and Processing (DSSP).⁴³^,⁴⁴ The predicted structure differed from the parental template model 2K3A in that the region between residues 86 and 102 contains a β-strand in 2K3A (Fig. 5B). However, the CHAP_K model has this region somewhat more distant from the sheet (Fig. 5A). As a result these residues are not being categorized as a β-strand in the CHAP_K model by DSSP⁴³^,⁴⁴ (Fig. 5A).

Figure 5. Topology maps for CHAP_K and 2K3A proteins. Secondary structure was determined using DSSP. The topology maps show α-helices in red and the strands of the β-sheets in yellow. (A) shows the topology map for CHAP_K while (B) (more ...)

Circular dichroism (CD) spectropolarimetry was performed to investigate the presence of an αβ fold, predicted in the CHAP_K model. The CD spectrum obtained for CHAP_K showed characteristic minima in the range 205–220 nm, typical of a protein containing α-helical and β-sheet structure (Fig. 6). Deconvolution of the spectrum using CDPro⁴⁵^,⁴⁶ and a test set comprising 43 soluble and 13 membrane proteins, indicated that CHAP_K contains an αβ fold, which supports the model. The exact α-helix and β-sheet content determined from the spectrum varied depending on which algorithm was used in the CDPro⁴⁵^,⁴⁶ suite of programs. These data are shown in Table 1. In the case of SELCON and CONTINLL algorithms, the content of α-helix is close to that found in the model (13.9%) while the CDSSTR algorithm is under-estimating the α-helical content significantly. For all components of CDPro⁴⁵^,⁴⁶ the estimated β-sheet structure is at variance with the amount of β-sheet seen in the model (33.9%, 35.2% and 39% for SELCON, CONTINLL and CDSSTR respectively, compared with an estimated 18.7% in the model). This could be explained by additional secondary structure in the protein that is not well modeled. The region from residues 86–102 (Fig. 5A) in the predicted CHAP_K model is an area where there is potentially an additional β-strand. This would bring the CHAP_K model closer in structure to the 2K3A structure. In addition, the N-terminal region of the 2K3A protein is an area of extended coil in the NMR structure, this may become more ordered and contribute additional β-sheet structure. Regardless of the precise secondary structure composition, CD spectropolarimetry indicated that the CHAP_K protein has an αβ fold, which is consistent with the predicted 3D model.

Figure 6. CD spectropolarimetry spectrum for CHAP_K.

Table 1. Secondary structure composition based on CD spectropolarimetry

Both target and template proteins have highly-conserved cysteine and histidine residues with Cys54-His117 residues of CHAP_K (magenta in Fig. 7) located approximately in the same position as the Cys57-His109 residues in the active site of 2K3A (cyan in Fig. 7). The distance between the Cys57 and His109 catalytic residues of 2K3A was found to be 3.3 Ångströms (Å) (Fig. 7), which is within hydrogen bond distance.³³ This distance is a critical feature in imparting enzymatic activity to CHAP proteins. Typically, the presence of a substrate in the Cys-His active site, triggers the deprotonation of cysteine by the histidine residue, resulting in a nucleophilic attack by the deprotonated cysteine on a substrate acyl carbonyl carbon, releasing a fragment of the substrate with an N-terminus. This leads to the formation of a tetrahedral ionic intermediate, restoration of the histidine residue to its deprotonated form and finally the release of the carboxylic acid moiety from the cysteine residue via hydrolysis. This results in the regeneration of the Cys-His active site.⁴⁷^,⁴⁸ The distance of 3.7 Å between the corresponding Cys54-His117 active site residues in the CHAP_K model are comparable with the active site of the 2K3A template (Fig. 7). It is therefore plausible that the CHAP_K Cys54-His117 highly conserved residues have a similar catalytic mode of action. When localized to their staphylococcal peptidoglycan substrate, they initiate a rapid nucleophilic attack on the peptide bond between _D-alanine and the first glycine in the pentaglycine cross-bridge, resulting in cleavage of the bacterial cell wall and rapid cell death. To further investigate the proposed role of the conserved Cys54-His117 active-site residues, a QuikChange II Site-Directed Mutagenesis Kit (Agilent Technologies Irealand Ltd) was used to create point mutations where a change from Cys54 to alanine (Ala) or His117 to Ala resulted in loss of enzymatic activity (data not shown). This strongly indicates that these two residues play a crucial role in CHAP_K catalytic activity.

Figure 7. Close-up stereoview of highly conserved Cys-His active site residues of CHAP_K (magenta) and 2K3A (cyan) superimposed. CHAP_K Cys54 (orange) and His117 (light blue), 2K3A Cys57 (yellow), His 109 (blue). Broken lines with labels indicate (more ...)

In total, eight conserved residues, Cys54, His117, Glu134, Asn 136, Gly116, Gly74, Thr59 and Gln55 at or near the surface of the predicted CHAP_K model may be essential for enzymatic activity. These residues coordinate the formation of a broad and shallow cleft which encompasses the active site of the enzyme (Fig. 8). This was consistent with the 2K3A template structure.³³ In the 2K3A protein, active site residues Cys57-His109-Glu126 are arranged in a viable proteolytic triad both in a sequence order and a topology typical of cysteine peptidases, with Asn128 acting as an additional stabilizing group.³³ Analysis of the CHAP_K model revealed corresponding highly conserved residues constituted by a Cys54-His117-Glu134 charge relay and an oxyanion hole (stabilizing) residue Asn136, near the surface of the active site cleft.

Figure 8. Electrostatic surface representation of CHAP_K with the active site cleft indicated by a black arrow. The protein surface is colored according to its electrostatic potential from regions of red (negatively charged) to blue (positively (more ...)

The exact mechanism of binding of CHAP_K to its peptidoglycan substrate remains unknown. However, examination of the electrostatic surface potential of the active site of CHAP_K revealed an overall electropositive charge (Fig. 8) indicating electrostatic complementarity with the electronegative character of its peptidoglycan substrate. This combination of attractive and repulsive forces may help guide the substrate into the correct position for binding before contacting the active site of CHAP_K. In addition, it was suggested that in the case of 2K3A, two exposed tyrosine aromatic groups, Tyr107 and Tyr129, either side of the active site possibly function as rails to anchor the substrate in position.³³ This may also be the case for CHAP_K as it has two exposed tyrosine residues (Tyr112 and Tyr140) on either side of its active site (Fig. 3A) and whose location corresponds closely to the position of Tyr107 and Tyr129 in the 2K3A template structure.

Conclusion

In the present study, comparative modeling was used to construct a 3D model of the phage derived peptidase, CHAP_K. Analysis of the 3D model provided a structural basis for CHAP_K activity and presented a probable mechanism involving the highly conserved Cys54-His117-Glu134-Asn136 active site residues. These residues have been identified in other lysins as the key catalytic residues involved in the hyrolysis of a _D-Ala-Gly substrate, in the pentaglycine cross-bridge of staphylococcal cell wall peptidoglycan. The proposed nucleophilic-attack mechanism of the conserved Cys54 residue on the peptide bond in an acid-base reaction is thought to be the most likely catalytic mechanism for all members of the CHAP NplC/P60 family of proteases.

Further exploration and investigation of phage lysins and their various domains, particularly at a structural and bioinformatic level, can provide important information for understanding biochemical function and specific substrate interaction mechanisms. This in turn may contribute to specifically engineered designer lysins with diverse and improved applications for development as novel antibacterials. Additional research is currently underway to experimentally resolve the crystal structure of CHAP_K.

Methods

The following web-based servers, suites, molecular analysis and visualization programs were used in the comparative modeling process for CHAP_K structure prediction. Basic Local Alignment Search Tool (BLAST) (www.ncbi.nlm.nih.gov/BLAST/), T-Coffee (www.tcoffee.org), MEGA4,³⁵ SWISS-MODEL (www.swissmodel.expasy.org),³⁶ iterative threading assembly refinement (I-TASSER),³⁸ PyMOL,⁴⁰ PROCHECK,³⁹ CNS_SOLVE⁴¹ LSQMAN⁴² and DSSP.⁴³^,⁴⁴ For CD spectropolarimetry The CHAP_K domain was dialysed extensively against 10 mM potassium phosphate buffer, pH 7.0. Protein was diluted to 0.019 mg/mL in the same buffer for analysis. Data was collected from 260 nm to 180 nm on a Chirascan instrument (Applied Photophysics). Data was converted to Δε values prior to analysis. Determination of the secondary structure from the CD spectrum was performed using the CDPro⁴⁵^,⁴⁶ software package.

Acknowledgments

This work was financed by Science Foundation Ireland Ref 06/RF/BIM004. We thank Dr. Todd Kagawa for his valuable input in the analysis of CD data.

Footnotes

Previously published online: www.landesbioscience.com/journals/bacteriophage/article/18245

References

1. Lowy FD. Staphylococcus aureus infections. N Engl J Med. 1998;339:520–32. doi: 10.1056/NEJM199808203390806. [PubMed] [Cross Ref]

2. European Centre for Disease Prevention and Control. Antimicrobial resistance surveillance in Europe 2009. Antimicrobial Resistance Surveillance Network (EARS-Net). Stockholm: ECDC; 2010. ISBN 978-92-9193-227-6.

3. Loessner MJ. Bacteriophage endolysins–current state of research and applications. Curr Opin Microbiol. 2005;8:480–7. doi: 10.1016/j.mib.2005.06.002. [PubMed] [Cross Ref]

4. Fischetti VA. Bacteriophage lysins as effective antibacterials. Curr Opin Microbiol. 2008;11:393–400. doi: 10.1016/j.mib.2008.09.012. [PMC free article] [PubMed] [Cross Ref]

5. Fenton M, Ross P, McAuliffe O, O'Mahony J, Coffey A. Recombinant bacteriophage lysins as antibacterials. Bioeng Bugs. 2010;1:9–16. doi: 10.4161/bbug.1.1.9818. [PMC free article] [PubMed] [Cross Ref]

6. Fischetti VA. Bacteriophage lytic enzymes: novel anti-infectives. Trends Microbiol. 2005;13:491–6. doi: 10.1016/j.tim.2005.08.007. [PubMed] [Cross Ref]

7. Borysowski J, Weber-Dabrowska B, Gorski A. Bacteriophage endolysins as a novel class of antibacterial agents. Exp Biol Med (Maywood) 2006;231:366–77. [PubMed]

8. García P, Garcia JL, Garcia E, Sanchez-Puelles JM, Lopez R. Modular organization of the lytic enzymes of Streptococcus pneumoniae and its bacteriophages. Gene. 1990;86:81–8. doi: 10.1016/0378-1119(90)90116-9. [PubMed] [Cross Ref]

9. Loessner MJ, Kramer K, Ebel F, Scherer S. C-terminal domains of Listeria monocytogenes bacteriophage murein hydrolases determine specific recognition and high-affinity binding to bacterial cell wall carbohydrates. Mol Microbiol. 2002;44:335–49. doi: 10.1046/j.1365-2958.2002.02889.x. [PubMed] [Cross Ref]

10. O'Flaherty S, Coffey A, Meaney W, Fitzgerald GF, Ross RP. The recombinant phage lysin LysK has a broad spectrum of lytic activity against clinically relevant staphylococci, including methicillin-resistant Staphylococcus aureus. J Bacteriol. 2005;187:7161–4. doi: 10.1128/JB.187.20.7161-7164.2005. [PMC free article] [PubMed] [Cross Ref]

11. Yokoi KJ, Kawahigashi N, Uchida M, Sugahara K, Shinohara M, Kawasaki K, et al. The two-component cell lysis genes holWMY and lysWMY of the Staphylococcus warneri M phage varphiWMY: cloning, sequencing, expression, and mutational analysis in Escherichia coli. Gene. 2005;351:97–108. doi: 10.1016/j.gene.2005.03.006. [PubMed] [Cross Ref]

12. Sass P, Bierbaum G. Lytic activity of recombinant bacteriophage phi11 and phi12 endolysins on whole cells and biofilms of Staphylococcus aureus. Appl Environ Microbiol. 2007;73:347–52. doi: 10.1128/AEM.01616-06. [PMC free article] [PubMed] [Cross Ref]

13. Rashel M, Uchiyama J, Ujihara T, Uehara Y, Kuramoto S, Sugihara S, et al. Efficient elimination of multidrug-resistant Staphylococcus aureus by cloned lysin derived from bacteriophage phi MR11. J Infect Dis. 2007;196:1237–47. doi: 10.1086/521305. [PubMed] [Cross Ref]

14. Obeso JM, Martinez B, Rodriguez A, Garcia P. Lytic activity of the recombinant staphylococcal bacteriophage PhiH5 endolysin active against Staphylococcus aureus in milk. Int J Food Microbiol. 2008;128:212–8. doi: 10.1016/j.ijfoodmicro.2008.08.010. [PubMed] [Cross Ref]

15. Gu J, Xu W, Lei L, Huang J, Feng X, Sun C, et al. A Novel Bacteriophage Lysin, LysGH15, Protects a Murine Bacteremia Model Efficiently against Lethal Methicillin-Resistant Staphylococcus aureus Infection. J Clin Microbiol. 2011;49:111–7. doi: 10.1128/JCM.01144-10. [PMC free article] [PubMed] [Cross Ref]

16. Jun SY, Jung GM, Son JS, Yoon SJ, Choi YJ, Kang SH. Comparison of the Antibacterial Properties of Phage Endolysins SAL-1 and LysK. Antimicrob Agents Chemother. 2011;55:1764–7. doi: 10.1128/AAC.01097-10. [PMC free article] [PubMed] [Cross Ref]

17. Navarre WW, Ton-That H, Faull KF, Schneewind O. Multiple enzymatic activities of the murein hydrolase from staphylococcal phage phi11. Identification of a D-alanyl-glycine endopeptidase activity. J Biol Chem. 1999;274:15847–56. doi: 10.1074/jbc.274.22.15847. [PubMed] [Cross Ref]

18. Becker SC, Dong S, Baker JR, Foster-Frey J, Pritchard DG, Donovan DM. LysK CHAP endopeptidase domain is required for lysis of live staphylococcal cells. FEMS Microbiol Lett. 2009;294:52–60. doi: 10.1111/j.1574-6968.2009.01541.x. [PubMed] [Cross Ref]

19. Rigden DJ, Jedrzejas MJ, Galperin MY. Amidase domains from bacterial and phage autolysins define a family of gamma-D,L-glutamate-specific amidohydrolases. Trends Biochem Sci. 2003;28:230–4. doi: 10.1016/S0968-0004(03)00062-8. [PubMed] [Cross Ref]

20. Horgan M, O'Flynn G, Garry J, Cooney J, Coffey A, Fitzgerald GF, et al. Phage lysin LysK can be truncated to its CHAP domain and retain lytic activity against live antibiotic-resistant staphylococci. Appl Environ Microbiol. 2009;75:872–4. doi: 10.1128/AEM.01831-08. [PMC free article] [PubMed] [Cross Ref]

21. Anantharaman V, Aravind L. Evolutionary history, structural features and biochemical diversity of the NlpC/P60 superfamily of enzymes. Genome Biol. 2003;4:R11. doi: 10.1186/gb-2003-4-2-r11. [PMC free article] [PubMed] [Cross Ref]

22. Bateman A, Rawlings ND. The CHAP domain: a large family of amidases including GSP amidase and peptidoglycan hydrolases. Trends Biochem Sci. 2003;28:234–7. doi: 10.1016/S0968-0004(03)00061-6. [PubMed] [Cross Ref]

23. Hermoso JA, Monterroso B, Albert A, Galan B, Ahrazem O, Garcia P, et al. Structural basis for selective recognition of pneumococcal cell wall by modular endolysin from phage Cp-1. Structure. 2003;11:1239–49. doi: 10.1016/j.str.2003.09.005. [PubMed] [Cross Ref]

24. Low LY, Yang C, Perego M, Osterman A, Liddington RC. Structure and lytic activity of a Bacillus anthracis prophage endolysin. J Biol Chem. 2005;280:35433–9. doi: 10.1074/jbc.M502723200. [PubMed] [Cross Ref]

25. Korndörfer IP, Danzer J, Schmelcher M, Zimmer M, Skerra A, Loessner MJ. The crystal structure of the bacteriophage PSA endolysin reveals a unique fold responsible for specific recognition of Listeria cell walls. J Mol Biol. 2006;364:678–89. doi: 10.1016/j.jmb.2006.08.069. [PubMed] [Cross Ref]

26. Porter CJ, Schuch R, Pelzek AJ, Buckle AM, McGowan S, Wilce MC, et al. The 1.6 A crystal structure of the catalytic domain of PlyB, a bacteriophage lysin active against Bacillus anthracis. J Mol Biol. 2007;366:540–50. doi: 10.1016/j.jmb.2006.11.056. [PubMed] [Cross Ref]

27. Korndörfer IP, Kanitz A, Danzer J, Zimmer M, Loessner MJ, Skerra A. Structural analysis of the L-alanoyl-D-glutamate endopeptidase domain of Listeria bacteriophage endolysin Ply500 reveals a new member of the LAS peptidase family. Acta Crystallogr D Biol Crystallogr. 2008;64:644–50. doi: 10.1107/S0907444908007890. [PubMed] [Cross Ref]

28. Fenton M, Casey PG, Hill C, Gahan CGM, Ross RP, McAuliffe O, et al. The truncated phage lysin CHAPk eliminates Staphylococcus aureus in the nares of mice. Bioeng Bugs. 2010;1:404–7. doi: 10.4161/bbug.1.6.13422. [PMC free article] [PubMed] [Cross Ref]

29. Baker D, Sali A. Protein structure prediction and structural genomics. Science. 2001;294:93–6. doi: 10.1126/science.1065659. [PubMed] [Cross Ref]

30. Claverie J-M, Notredame C. Bioinformatics for dummies 2nd ed.–For dummies. 2007, Hoboken, N.J.: Wiley Pub. xviii, 436 p.

31. Martí-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, Sali A. Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 2000;29:291–325. [PubMed]

32. Ginalski K. Comparative modeling for protein structure prediction. Curr Opin Struct Biol. 2006;16:172–7. doi: 10.1016/j.sbi.2006.02.003. [PubMed] [Cross Ref]

33. Rossi P, Aramini JM, Xiao R, Chen CX, Nwosu C, Owens LA, et al. Structural elucidation of the Cys-His-Glu-Asn proteolytic relay in the secreted CHAP domain enzyme from the human pathogen Staphylococcus saprophyticus. Proteins. 2009;74:515–9. doi: 10.1002/prot.22267. [PMC free article] [PubMed] [Cross Ref]

34. Armougom F, Moretti S, Poirot O, Audic S, Dumas P, Schaeli B, et al. Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee. Nucleic Acids Res. 2006;34(Web Server issue):W604-8. doi: 10.1093/nar/gkl092. [PMC free article] [PubMed] [Cross Ref]

35. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9. doi: 10.1093/molbev/msm092. [PubMed] [Cross Ref]

36. Kiefer F, Arnold K, Kunzli M, Bordoli L, Schwede T. The SWISS-MODEL Repository and associated resources. Nucleic Acids Res. 2009;37(Database issue):D387–92. doi: 10.1093/nar/gkn750. [PMC free article] [PubMed] [Cross Ref]

37. Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T. Protein structure homology modeling using SWISS-MODEL workspace. Nat Protoc. 2009;4:1–13. doi: 10.1038/nprot.2008.197. [PubMed] [Cross Ref]

38. Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40. doi: 10.1186/1471-2105-9-40. [PMC free article] [PubMed] [Cross Ref]

39. Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996;8:477–86. doi: 10.1007/BF00228148. [PubMed] [Cross Ref]

40. DeLano WL. The PyMOL Molecular Graphics System 2002, DeLano Scientific, San Carlos, CA, USA.

41. Brünger AT, Adams PD, Clore GM, DeLano WL, Gros P, Grosse-Kunstleve RW, et al. Crystallography & NMR system: A new software suite for macromolecular structure determination. Acta Crystallogr D Biol Crystallogr. 1998;54:905–21. doi: 10.1107/S0907444998003254. [PubMed] [Cross Ref]

42. Novotny M, Madsen D, Kleywegt GJ. Evaluation of protein fold comparison servers. Proteins. 2004;54:260–70. doi: 10.1002/prot.10553. [PubMed] [Cross Ref]

43. Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983;22:2577–637. doi: 10.1002/bip.360221211. [PubMed] [Cross Ref]

44. Joosten RP, te Beek TA, Krieger E, Hekkelman ML, Hooft RW, Schneider R, et al. A series of PDB related databases for everyday needs. Nucleic Acids Res. 2011;39(Database issue):D411–9. doi: 10.1093/nar/gkq1105. [PMC free article] [PubMed] [Cross Ref]

45. Sreerama N, Woody RW. On the analysis of membrane protein circular dichroism spectra. Protein Sci. 2004;13:100–12. doi: 10.1110/ps.03258404. [PubMed] [Cross Ref]

46. Sreerama N, Woody RW. Computation and analysis of protein circular dichroism spectra. Methods Enzymol. 2004;383:318–51. doi: 10.1016/S0076-6879(04)83013-1. [PubMed] [Cross Ref]

47. Lau EY, Bruice TC. Consequences of breaking the Asp-His hydrogen bond of the catalytic triad: effects on the structure and dynamics of the serine esterase cutinase. Biophys J. 1999;77:85–98. doi: 10.1016/S0006-3495(99)76874-8. [PubMed] [Cross Ref]

48. Lin CH, Kwon DS, Bollinger JM, Jr., Walsh CT. Evidence for a glutathionyl-enzyme intermediate in the amidase activity of the bifunctional glutathionylspermidine synthetase/amidase from Escherichia coli. Biochemistry. 1997;36:14930–8. doi: 10.1021/bi9714464. [PubMed] [Cross Ref]

Articles from Bacteriophage are provided here courtesy of
Landes Bioscience