Modular endolysin of Burkholderia AP3 phage has the largest lysozyme-like catalytic subunit discovered to date and no catalytic aspartate residue

Endolysins are peptidoglycan-degrading enzymes utilized by bacteriophages to release the progeny from bacterial cells. The lytic properties of phage endolysins make them potential antibacterial agents for medical and industrial applications. Here, we present a comprehensive characterization of phage AP3 modular endolysin (AP3gp15) containing cell wall binding domain and an enzymatic domain (DUF3380 by BLASTP), both widespread and conservative. Our structural analysis demonstrates the low similarity of an enzymatic domain to known lysozymes and an unusual catalytic centre characterized by only a single glutamic acid residue and no aspartic acid. Thus, our findings suggest distinguishing a novel class of muralytic enzymes having the activity and catalytic centre organization of DUF3380. The lack of amino acid sequence homology between AP3gp15 and other known muralytic enzymes may reflect the evolutionary convergence of analogous glycosidases. Moreover, the broad antibacterial spectrum, lack of cytotoxic effect on human cells and the stability characteristics of AP3 endolysin advocate for its future application development.


Results
In silico analysis of AP3 endolysin. Phage AP3, a representative of Gram-negatives infecting phages, utilizes a complex of four proteins which ensures efficient progeny release from an infected host cell. Its endolysin (AP3gp15), together with antiholin (AP3gp13), holin (AP3gp14) and bimolecular spanin (AP3gp16 and AP3gp17) form the lysis cassette 22 . AP3gp15 is a 266-amino acid protein with molecular mass of 28.9 kDa and theoretical pI of 8.82. The structure of endolysin, predicted in silico in BLASTP, is modular, consisting of CBD (10-65 aa; PG_binding_1 domain; pfam01471) at the N-terminus, and DUF3380 domain (pfam11860) at the C-terminus (90-262 aa) (Fig. 1A). There are at least 73 known significant homologues to modularly organized AP3gp15 and 37 homologues with reversely arranged domains. There are also 2428 and 233 homologues (E-value < 0.0001) of N-terminal CBD and C-terminal EAD, respectively, if compared separately. Domains similar to AP3gp15 CBD have been already well characterized in terms of the structure and function 23,24 , however, structural information concerning DUF3380 is non-existent. Recently, the gp110 endolysin derived from Salmonella phage 10, possessing DUF3380 has been analysed regarding its cleavage specificity and stability features 15 . However, the amino acid identity between phage 10 gp110 and AP3gp15 is less than 50% (46% for the whole protein sequence and 49% for CBD and 45% for EAD, if compared separately) (Fig. 1A). Moreover, both proteins differ significantly in terms of activity strength and stability features. Referring to above differences, the determination of the cleavage specificity and stability of AP3gp15 had to be examined (see below).
The modular structure of AP3 endolysin. To provide the structural basis of the catalytic activity of DUF3380, the X-ray crystal structure of AP3gp15 has been evaluated ( Fig. 1B; Table 1). The overall structure comprises of two clearly distinguished domains, a smaller N-terminal CBD domain (Lys3 -Ser66), and a larger EAD domain (Asp77 -Ala264) located towards the C-terminus (Fig. 1B). CBD domain is composed of a bundle of three antiparallel helices (α1-α3) and is joined with EAD via a hinge region spanning amino acids (Ala67 -Thr76). The EAD domain is organized into two sub-domains: first encompassing helices α4-α5 (Asp77 -Glu101) and helices α10 to α15 (Gly185 -Ala266), and the second containing β1 (Ile115 -Tyr117), β2-β3 (Thr176 -Met184) sheet intervened by helices α6 to α9 (Arg119 -Ala175). The structure of the EAD domain of AP3gp15 (black) and its equivalent Thr (T) in 10gp110 (red) are framed. (B) Two asymmetrical domains are visible, CBD (N-terminal bundle of three helices representing the putative PG binding domain) encompassing residues Lys3 -Ser66; the catalytic domain EAD (residues Asp77 -Ala264); the hinge region (residues Ala67 to Thr76) connects both subunits.  26 and Pseudomonas aeruginosa phage phiKZ gp144 endolysin (PDB ID 3BKH; RMSD = 1.1 for 65 C α atoms (94%)) 24 (Fig. 2). Despite a relatively low homology of amino acid sequence (31% and 39% identity, respectively), a high structural homology of the entire AP3gp15 N-terminal part to the domains of clearly identified function strongly suggests its role in PG interaction. The surface residue analysis of CBD, similarly to homologues, did not reveal any distinct negatively or positively charged patches. However, surface amino acids (Arg6, Arg10, Asp13, Asp31, and Asp54) were spatially homologous to Arg5, Arg9, Asp11, Asp31 and Asp 54 of Zn-dependent amidase, respectively. Also, the residues of AP3gp15: Arg6, Asp13, Asp31, Asp40, Asp54 have their equivalents in phage phiKZ endolysin residues: Lys58, Asp65, Asp83, Asp91, and Asp106. The most common fingerprint between those three proteins comprises of Arg6, Asp31, and Asp54 of AP3gp15 CBD. The conservative character of residues likely suggests their role in PG recognition and binding.
Structural homology suggests an unusual catalytic mechanism of the AP3gp15 enzymatic domain. The structure alignment of AP3 endolysin to HEWL, GEWL, T4 lysozyme and lambda lytic transglycosylase suggests the probable catalytic center elements interacting with the substrate (Fig. 3B,D,G,H). Glu35 and Asp52 of HEWL, Glu73 of GEWL, Glu11 and Asp20 of T4L, and Glu19 of lambda endolysin have been identified as the key residues involved in the catalytic process 29,30 . AP3gp15 contains only Glu equivalent at the presumed catalytic centre (position 101). No Asp equivalent, as the second catalytic residue, could have been identified in the active center of AP3gp15.The substrate is most likely recognized by residues located in α5, α10, a12, α14 helices and a loop region (Ser102 -Val114) of AP3gp15. In particular, Asp221 and His187 as homologues of Asp101 and Trp62-63 (involved in recognition and binding of PG sugar residues by HEWL 27 ) were identified as taking part in the binding of the oligosaccharide substrate (Fig. 4).
AP3 endolysin has lysozyme specificity. The enzymatic specificity of AP3 endolysin was examined in vitro using purified PG as a substrate with previously described methodology 31 . The chromatograms (Fig. 5A) demonstrated AP3gp15 muropeptide profiles similar to those obtained after lysozyme treatment suggesting similar specificity of both enzymes. Further confirmation by MS analysis of all major peaks present in both samples ( Fig. 5B) identified identical products. Thus, comparable to lysozyme, AP3gp15 cleaves β-1,4-glycosidic bond of PG and releases GlcNAc and MurNAc. These results confirm that AP3 endolysin is, in fact, a lysozyme, not a lytic transglycosylase.
Activity and cytotoxicity analysis of recombinant AP3 endolysin. The muralytic activity of AP3 endolysin was assessed on permeabilized Gram-negative strains according to a standardized assay methodology 32 . The AP3gp15 turned out to be two times stronger than the commercially available lysozyme (Table 2). This was consistent for all tested strains, including E. coli, K. pneumoniae, P. aeruginosa, B. cenocepacia and S. enterica Typhimurium. No lytic activity was detected against Gram-positive Staphylococcus aureus and S. epidermidis strains. To verify that Glu101 is a key residue of AP3gp15 catalytic centre, the site-directed mutagenesis was performed and the Glu101 has been converted to an alanine (Glu101Ala). The muralytic activity of modified enzyme tested on P. aeruginosa PAO1 permeabilized cells, dropped down from 14,710 U/mg for wild-type version to 90 U/mg for protein mutant, reducing the lytic activity of 99.4%.
The endolysin showed to be relatively stable at low temperatures but heat-sensitive. It retained > 95% activity after 1 week storage at −20 °C or 4 °C, at the concentration of 0.5 mg/ml in PBS buffer (pH 7.4). After 1 week at RT, the protein started to precipitate but still retained 70% residual activity. Ten minutes incubation at higher temperatures: 30 °C, 40 °C, 50 °C, 60 °C resulted in a decrease of around 2%, 25%, 80% and 92% of its hydrolytic activity, respectively. The incubation for 10 minutes at 80 °C caused a complete inactivation. The pH optimum was in the range of 7-9. The activity dropped at pH 6 and 10 by 10% and 25%, respectively. Low pH (pH = 5) reduced the activity to 58%. The cytotoxicity of AP3gp15 on mammalian cells was also evaluated and 50 µg/ml of the enzyme exhibited no adverse effect on A549 and THP-1 cells. After 48 h of incubation, no changes in cell viability and number were observed in comparison to endolysin free buffer control (data not shown).

Discussion
Phage endolysins have a large, but yet unexplored potential as antimicrobials, due to poor structural and biochemical characterization. The annotation of endolysin encoding genes and their qualification to a specific enzymatic group are generally based on the amino acid sequence similarity to known proteins. Lack of a large reference library results, among others, in a false classification of newly characterized enzymes. One of the examples is the phage lambda endolysin which operates in the databases and in numerous publications as lysozyme, but in fact, has lytic transglycosylase activity 29,33 . The experimental verification of structural features and enzymatic specificity is a crucial step for prospective enzyme modification and further application. The thorough characterization of representative members is a key element in the case of conservative domains of unknown function (DUFs) which represent 20% of all currently annotated protein domains in Pfam database. The data obtained from the first member usually defines the primary characteristics of the entire protein family allowing a better understanding of the biology/function of many organisms.
In the present study, we describe a modular endolysin (rare for Gram-negative infecting phages) originating from B. cenocepacia AP3 phage. The domains arrangement of AP3gp15 (CBD at the N-terminus and EAD at the C-terminus) is characteristic for all eight endolysins originating from Gram-negatives specific phages with confirmed modular structure 10,[12][13][14][15]34 . The AP3gp15 CBD element revealed a high structural similarity to CBDs of Zn-dependent amidase from C. acetobutylicum and Pseudomonas phage phiKZ endolysin. The latter domain has a high-affinity to Gram-negatives PG and is able to increase an enzymatic activity of phiKZ endolysin 23 . It should Homologies between HEWL and AP3gp15 residues that are possibly involved in the substrate recognition and other interactions. Only the homologous catalytic Glu35 of HEWL is highlighted as it is represented by its counterpart residue of Glu101 in AP3gp15. Asp101 of HEWL plays a role in substrate positioning and its putative homologue is Asp221 of AP3gp15. The closest homologue of similarly functioning Trp62 and Trp63 of HEWL is represented by His187 that positions between these two amino acids in the structural alignment.
be noted that genomes of both phages: phiKZ (giant myovirus, 280,334 bp) and AP3 (Peduovirinae; 36,499 bp) belong to entirely different taxonomic lineages. Moreover, CBD of AP3 endolysin shows significant similarity to almost two and a half thousand homologues (phages and bacteria) in genome databases including a recently described endolysin from Salmonella phage 10 15 . Therefore, AP3gp15-like CBDs seem to be universal in general, irrespectively to the origin. It stays in contrast to endolysins encoded by Gram-positives specific phages, where CBDs demonstrate a large sequence diversity, and genus or even species specificity 4,35 . The structural alignment highlights the similarity of DUF3380 of AP3gp15 to proteins having an activity of lysozyme and lytic transglycosylase (phage lambda endolysin). The catalytic mechanism of a vast majority of lysozymes relies on two residues: Glu and Asp. Catalytic centers of canonical lysozymes derived from T4 and P21 phages are based on the Asp residue separated in a sequence by eight residues relative to catalytic Glu. Although, the AP3gp15 also contains Asp at analogous sequence position (Asp110), at the structural level this residue is not equivalent. In T4 and P21 lysozymes, the residue forms the active site whereas in AP3gp15 it is positioned away from the catalytic center. The crystal structure analysis suggests that only a single catalytic residue (Glu101) is present at the active site of AP3 endolysin. Moreover, the comparative in silico analysis of AP3gp15 EAD homologues revealed that Asp is not a conserved residue in this protein family (present only in 123 out of 233 analyzed DUF3380 domains) what again might suggest an insignificant role in the catalysis. Instead of given Asp, other replacing amino acids: asparagine, threonine, serine, histidine, proline, tyrosine or alanine, can be found in sequenced DUF3380 homologues. Of functionally characterized member, a DUF3380 domain of Salmonella phage 10 endolysin has threonine residue at a position equivalent to 110 in AP3gp15 (Fig. 1A) and its active site most probably consists only of Glu101 residue, although no structure is available. Experimental data of AP3gp15 mutant Glu101Ala with a decline in activity by 99.4%, confirmed of Glu101 involvement in the catalytic reaction. However, based on the putative catalytic residues identified on other lysozymes by sequence and structure similarity, we were not able to identify a second catalytic residue. To date, only one phage endolysin (lysozyme Lys68) with Glu as a putative single catalytic residue, has been described, but no sequence homology to AP3gp15 is seen 10 . A single glutamic acid is characteristics for the catalytic centre of GEWL group of lysozymes 20 and lytic transglycosylases 36 . Although the structural comparison of AP3gp15 EAD suggests a closer similarity to phage lambda endolysin (lytic transglycosylase) rather than GEWL, the biochemical analysis showed lysozyme-like cleavage products of PG digestion by AP3gp15. Similar results showing a lysozyme cleavage specificity were obtained for a related modular endolysin derived from Salmonella phage 10 by Rodríguez-Rubio and co-workers 15 . This indicates that DUF3380 has a HEWL-like fold, but a novel type of catalytic center basing on one residue as in GEWL and lytic transglycosylase. The structural similarity of the catalytic domain AP3gp15 to other glycosidases allows speculating on the catalytic mechanism. The substrate binding and related distortion of the cleaved bond, the formation of a covalent intermediate and the electrophilic movement of C-1 carbon along the reaction coordinate are the canonical basis 30 . HEWL was the first enzyme for which the three-dimensional structure was determined by X-ray diffraction 37 and currently it is the best structurally characterized glycosidase in terms of the catalytic mechanism and substrate binding. The cleavage is preceded by the binding of PG hexasaccharide unit with concomitant distortion within MurNAc. MurNAc unit in the -1 (Glu) subsite adopts a half-chair conformation, thereby forming a kink in the oligosaccharide between sites -1 and + 1 (Glu and Asp). Glu35 of HEWL transfers a proton to the O1 position while the negative charge on Asp52 stabilizes the positively charged oxonium ion intermediate 30,38,39 . Despite the catalytic domain of AP3 endolysin being much larger than HEWL and the core of AP3gp15 aligns well with HEWL, the core of AP3gp15 contains additional protrusions which function is unknown. Although, the significantly different binding mode to HEWL is unlikely, the substrate binding surface of AP3gp15 may be extended in comparison to HEWL. If so, EAD of AP3 endolysin (DUF3380) could be the largest catalytic subunit described to date among lysozymes. The broad distribution and conservative character of DUF3380 domains in phages and bacteria advocate its important function. Moreover, the lack of amino acid sequence homology between AP3 endolysin and other glycosidases may reflect evolutionary convergence. That would imply that DUF3380-like and other classes of muralytic glycosidases had aroused in nature independently in different ancestral phages/bacteria as a result of an evolutionary trend in the development of peptidoglycan-degrading proteins. Particularly interesting is the fact, that despite differences in the known structures of lysozymes, the enzymatic specificity remains the same.
The AP3gp15 exhibits almost twice higher muralytic activity comparing to commercially available HEWL lysozyme, both tested on previously permeabilized bacterial cells. Compared to other endolysins derived from Gram-negatives specific phages, the AP3gp15 was classified as a medium effective enzyme, having 65-fold lower than the strongest 10gp110, to 6-fold higher activity comparing to the weakest KP32gp15 10,[13][14][15]31,34 . Similarly to others, the AP3 endolysin shows a broad spectrum of activity against Gram-negative strains, but no PG-degrading properties with respect to Gram-positives 10,13,14,31,40,41 . This is consistent with the differences in PG structure. The thick PG of Gram-positive bacteria varies significantly in the peptide composition, crosslinks, and modifications of the glycan chain, whereas the PG of Gram-negatives is conserved having 1-3 layers of A1γ chemotype 42 . Like the vast majority of endolysins derived from Gram-negative specific phages, the AP3gp15 is not able to lyse intact  cells and for this purpose requires the combination with the outer membrane permeabilizers or the binding to specific cationic peptide elements [9][10][11] . It is worth to mention that despite amino acid homology between AP3 and Salmonella phage endolysins, they differ in terms of bactericidal potential and thermostability. The AP3gp15 has 65-fold lower PG-hydrolysing efficacy than Salmonella phage endolysin (958,720 U/mg). Moreover, in contrast to the latter one, the AP3gp15 is unstable above 60 °C. The stability of Salmonella phage endolysin is probably related to the presence of disulphide bridges (DiANNA 1.1 web server prediction) 43 which is absent in AP3gp15. The accurate determination of protein thermostability is an important issue in protein characterization and allows for the rational optimization of bioprocesses and for prospective protein modification to increase their thermal resistance [44][45][46] . Possibly, the introduction of disulphide bridges guided by Salmonella phage endolysin into AP3gp15 would increase its stability and thus the utility for industrial application. No adverse effect on eukaryotic cell lines (lack of cytotoxicity) promotes the AP3gp15 potential as an antimicrobial agent, as well. Summarizing, in silico comparative analysis of AP3gp15 CBD and EAD (DUF3380) domains demonstrated unexpectedly a widespread occurrence within phages and bacteria and to a lesser extent also among archaea. The conserved character and wide distribution of both domains suggest their important function and evolutionary ancient origin. Therefore, the first report of DUF3380 domain structure and experimental characterization of its enzymatic activity presented in this study provides a significant contribution to the understanding of this widespread protein family with muralytic properties.

Methods
In silico analysis of AP3gp15. In silico analysis of AP3 endolysin included comparative protein sequence homology analysis, conserved domains recognition (NCBI's BLASTP, HMMER, Phyre2, HHpred) and predicted physicochemical parameters (ExPASy ProtParam) using methodology described elsewhere 47-50 . Recombinant endolysin preparation. The AP3 endolysin (AP3gp15) sequence was amplified from AP3 genomic DNA by PCR using Pfu DNA polymerase (Thermo Fischer Scientific) and specific primers: forward: 5′-ATGTATAAAACCCTGCGCCTCGGCG-3′, 10 pM; reverse: 5′-CGCGGCCGCCCGGCTGTAGCGATCG-3 ′, 10 pM. The PCR product was cloned into commercially available pEXP5-CT/TOPO expression vector (Thermo Fischer Scientific) according to the manufacturer's instruction. The construct was verified by sequencing and subsequently transformed into Escherichia coli BL21-AI (Thermo Fischer Scientific) expression strain. The culture was performed in LB medium at 37 °C with shaking (200 rpm) to achieve of mid-exponential-phase (OD600 ~ 0,6). Next, the expression was induced with L-arabinose (to a final concentration of 50.2%) and the incubation was continued at 20 °C for 18 hours. Cells were pelleted at 6000 x g, resuspended in lysis buffer (500 mM NaCl, 20 mM NaH 2 PO 4 , pH 7.7) and disrupted by the combination of freeze-thawing (−80 °C/RT; three times) and sonication (5 cycles of 15 s pulse and 30 s rest, Bandelin Sonopuls). Subsequently, the lysate was clarified at 15000 x g and recombinant protein was purified from filtered supernatant on NGC Medium Pressure Chromatography Systems (Bio-Rad) combined with Bio-Scale Mini Profinity IMAC Cartridges (Bio-Rad). The fractions containing recombinant protein were eluted (500 mM NaCl, 20 mM NaH 2 PO 4 , 500 mM imidazole, pH 7.7) and pooled. Protein  Table 2). In Gram-negative strains the outer membrane was removed by 45 min incubation in chloroform-saturated 50 mM Tris-HCl buffer (pH 7.7), according to Lavigne et al. 52 methodology. The muralytic activity of AP3 phage endolysin was assayed as described by Briers et al. 32 . 270 µl of bacterial culture was combined with 30 µl of endolysin (final concentration range of 50 ng/ml -200 ng/ml). As a control, the same concentrations of commercially available lysozyme from chicken egg white (Sigma Aldrich) was used. The kinetic drop in culture turbidity in the presence of endolysin or commercial lysozyme was measured spectrophotometrically using microplate reader (Asys UVM340). The muralytic activity was quantified using a standardized method developed by Briers et al. 32 , where 1 unit corresponds to the amount of endolysin resulting in an OD 600 decrease in culture turbidity of 0.001/min. Determination of AP3 phage endolysin stability. Storage stability in PBS buffer (pH 7.4) was determined by incubating the endolysin for 1 week at RT and 1 month at 4 °C and −20 °C. The thermostability of enzyme was determined in PBS (pH 7.4) by 10 min incubation at 30 °C, 40 °C, 50 °C, 60 °C and 70 °C. The pH stability was defined in the universal buffer (50 mM KCl, 10 mM KH 2 PO 4 , 10 mM Na-citrate and 10 mM H 3 BO 4 ) adjusted to different pH (from 3 to 10 using) by concentrated NaOH or HCl. In each case, the residual activity was determined by monitoring the PG degrading activity against permeabilized P. aeruginosa PAO.
PG isolation and analysis of endolysin cleavage specificity. E. coli cells from 0.2 l overnight culture were pelleted at 4,500 x g and resuspended in 5 ml of PBS. An equal volume of 10% SDS was added and the sample was incubated in a boiling water bath with vigorous stirring for 4 h. Following the sample was further stirred overnight at RT. The insoluble fraction (PG) was pelleted at 400,000 x g, 15 min, 30 °C (TLA-100.3 rotor; OptimaTM Max ultracentrifuge, Beckman). SDS was washed out and the PG was treated with Pronase E 0.1 mg/ml at 60 °C for 1 h and further boiled in 1% SDS for 2 h to stop the reaction. The sacculus was resuspended in 200 μl of 50 mM sodium phosphate buffer pH 4.9 and digested overnight with 30 μg/ml muramidase (Cellosyl) or with 0.4 mg/ml of AP3gp15 in Tris-HCL 20 mM buffer pH 8.0, 1 mM MgCl 2 and 1 mM ZnCl 2 . Samples were incubated at 37 °C. PG digestion was stopped by 5 min incubation in a boiling water bath. Coagulated protein was removed by centrifugation. The supernatants were mixed with 150 μl 0.5 M sodium borate pH 9.5, and subjected to reduction of muramic acid residues into muramitol by sodium borohydride treatment (10 mg/ml final concentration, 30 min at RT). Samples were adjusted to pH 3.5 with phosphoric acid. Chromatographic analyses of muropeptides were performed on ACQUITY Ultra Performance Liquid Chromatography (UPLC) BEH C18 column (130 Å, 1.7 μm, 2.1 mm by 150 mm; Waters), and peptides were detected at Abs. 204 nm using ACQUITY UPLC UV-Visible Detector. Muropeptides were separated using a linear gradient from buffer A (PBS 50 mM, pH 4.35) to buffer B (PBS 50 mM, pH 4.95 methanol 15% (v/v)) in 20 min, and flow 0.25 ml/min. The peaks from chromatogram that did not show similarity to any of previously described muropeptides were considered as background noise and removed from the chromatogram. Eluting peak samples were analyzed with Agilent 6550 iFunnel Q-TOF LC/MS System (Agilent Technologies) using an organic separation method 53 . Crystallization, X-ray data collection, and structure determination. Native and SeMet-labelled The diffraction data were collected at beamline P11 PETRA III at DESY. A native data set extending to 1.72 Å resolution was collected using 0.9184 Å wavelength, while a data set for SeMet-labelled protein crystals were collected at the peak wavelength (0.9795 Å) and extended to 2.42 Å resolution. Data were integrated and scaled using XDS 54 and SCALA 55 from the CCP4 suite, respectively. The native crystals belonged to space group H32, and the Matthews coefficient suggested the presence of two monomers in the asymmetric unit. SeMet-labelled crystals belonged to space group P1211 with four monomers in the asymmetric unit.
The initial structure was determined by single-wavelength anomalous dispersion (SAD) using AutoSol program included in the PHENIX suite 56 . The model was constructed and partially refined in an automated mode using AutoBuild and Buccaneer 57 from PHENIX and CCP4 suites, respectively. Initial, partially refined structure was used as a search model for molecular replacement using native data and performed with Phaser 58 . Refinement was performed with REFMAC5 59 to complement manual rebuilding using COOT 60 . R free was used to monitor the refinement strategy. Water molecules were added in 2Fobs -Fcalc map at densities above 2σ at reasonable sites upon visual inspection. Data processing and refinement statistics are summarized in Table 1. The geometry of the final model was analyzed with Molprobity 61 . Structural similarity searches were carried out using Dali server 25 and Z scores of > 10 were considered as significant.
Site-directed mutagenesis. To establish the putative catalytic residue within the endolysin sequence, the active site mutation (Glu101Ala) was introduced. Two overlapping specific primers: forward 5′-GCCGTCAATGAGGTTGCATCGAAAGGTGCCGGG-3′ (10 pM) and reverse 5′-CCCGGCACCTTT CGATGCAACCTCATTGACGGC-3′ (10pM) with mutation basepairs underlined, were applied. The previously constructed expression vector (pEXP5-CT/TOPO plasmid containing the AP3gp15 DNA) was mutated using above primers and the GeneArt ® Site-Directed Mutagenesis System from Thermo Fischer Scientific. The mutant was prepared, purified using the same protocol as for the wild type of AP3gp15. The PG muralytic activity of the purified AP3gp15 Glu101Ala mutant was quantified using permeabilized P. aeruginosa PAO1 cells, and protocol described below in "determination of muralytic potential" section.
Determination of endolysin cytotoxicity. Determination of endolysin cytotoxicity against human cells has been chosen as preliminary safety evaluation. Cytotoxicity of AP3gp15 was examined with respect to lung carcinoma epithelial cell line A549 (ATCC CCL-185) and acute monocytic leukaemia cell line THP-1 (ATCC TIB-202) using trypan blue assay. Endolysin suspended in PBS buffer pH 7.4 was added to exponentially growing cells (1 × 10 5 cells/ml) to a final concentration of 50 µg/ml. Cells were incubated for 24 h and 48 h in a suitable medium: DMEM (for A549) or RPMI 1640 (for THP-1), both supplemented with 2 mM glutamax, 1% antibioticantimycotic, and 10% heat-inactivated fetal bovine serum (all media and supplements from Gibco BRL, Thermo Fischer Scientific). Cell viability was determined by trypan blue staining (0.4% in PBS; Sigma Aldrich, Germany) and compared to a number of cells in the control (PBS treatment only). Data availability. The protein sequence of AP3gp15 can be found in GenBank under accession number: