A novel directed evolution method specific to the field of thermostability engineering was developed and consequently validated for bacteriolytic enzymes. After only one round of random mutagenesis, an evolved bacteriolytic enzyme, PlyC 29C3, displayed greater than twice the residual activity when compared to the wild-type protein after elevated temperature incubation.
Directed evolution is defined as a method to harness natural selection in order to engineer proteins to acquire particular properties that are not associated with the protein in nature. Literature has provided numerous examples regarding the implementation of directed evolution to successfully alter molecular specificity and catalysis1. The primary advantage of utilizing directed evolution instead of more rational-based approaches for molecular engineering relates to the volume and diversity of variants that can be screened2. One possible application of directed evolution involves improving structural stability of bacteriolytic enzymes, such as endolysins. Bacteriophage encode and express endolysins to hydrolyze a critical covalent bond in the peptidoglycan (i.e. cell wall) of bacteria, resulting in host cell lysis and liberation of progeny virions. Notably, these enzymes possess the ability to extrinsically induce lysis to susceptible bacteria in the absence of phage and furthermore have been validated both in vitro and in vivo for their therapeutic potential3-5. The subject of our directed evolution study involves the PlyC endolysin, which is composed of PlyCA and PlyCB subunits6. When purified and added extrinsically, the PlyC holoenzyme lyses group A streptococci (GAS) as well as other streptococcal groups in a matter of seconds and furthermore has been validated in vivo against GAS7. Significantly, monitoring residual enzyme kinetics after elevated temperature incubation provides distinct evidence that PlyC loses lytic activity abruptly at 45 °C, suggesting a short therapeutic shelf life, which may limit additional development of this enzyme. Further studies reveal the lack of thermal stability is only observed for the PlyCA subunit, whereas the PlyCB subunit is stable up to ~90 °C (unpublished observation). In addition to PlyC, there are several examples in literature that describe the thermolabile nature of endolysins. For example, the Staphylococcus aureus endolysin LysK and Streptococcus pneumoniae endolysins Cpl-1 and Pal lose activity spontaneously at 42 °C, 43.5 °C and 50.2 °C, respectively8-10. According to the Arrhenius equation, which relates the rate of a chemical reaction to the temperature present in the particular system, an increase in thermostability will correlate with an increase in shelf life expectancy11. Toward this end, directed evolution has been shown to be a useful tool for altering the thermal activity of various molecules in nature, but never has this particular technology been exploited successfully for the study of bacteriolytic enzymes. Likewise, successful accounts of progressing the structural stability of this particular class of antimicrobials altogether are nonexistent. In this video, we employ a novel methodology that uses an error-prone DNA polymerase followed by an optimized screening process using a 96 well microtiter plate format to identify mutations to the PlyCA subunit of the PlyC streptococcal endolysin that correlate to an increase in enzyme kinetic stability (Figure 1). Results after just one round of random mutagenesis suggest the methodology is generating PlyC variants that retain more than twice the residual activity when compared to wild-type (WT) PlyC after elevated temperature treatment.
1. Determining Optimal Heating Conditions
First, one must experimentally determine the optimal incubation temperature and time to use for the heating step in the assay. For our PlyC model, it is important to note that E. coli co-transformed with plyCA and plyCB genes on separate expression plasmids has been shown to form a fully functional PlyC holoenzyme6. The 96 well microtiter plate preparation as well as cell growth conditions and the subsequent replica plating technique were adapted from examples provided in the literature12-16. An incubation time of 30 min will be suitable for this assay as results from previous thermal inactivation experiments dictate a rapid loss in activity during short-term incubation in environments exceeding physiological temperature (unpublished observation). The optimal incubation temperature for PlyC was elucidated by the following steps:
After a broad temperature screen is performed to identify the non-permissive temperature, the above steps can be repeated over a more narrow range of temperatures (5 °-10 °C) to elucidate the precise non-permissive temperature for the enzyme of interest. Note, the non-permissive temperature identified in this assay may be different than the melting temperature elucidated by other means due to differences in volume, concentration, etc.
2. Generating the Mutant Library
The creation of the mutant library to be screened involves randomly incorporating mutations by an error-prone DNA polymerase with minimal nucleotide bias in the plyCA gene using the GeneMorph II Random Mutagenesis Kit as follows:
3. 96 Well Microtiter Plate Preparation and Cell Growth Conditions
4. Replica Plating, Protein Expression Induction and Lysate Preparation
5. Soluble Lysate Heat Treatment
At the conclusion of the first round of random mutagenesis, over 6,000 PlyC mutants were screened and a total of 35 mutants with potentially increased thermal behavior were identified, selected and sequenced. Genomic analysis, summarized in Table I, suggests that of the 35 candidates, 7 of the constructs contained WT PlyCA sequences at the level of translation, corresponding to false positives identified by the assay. Of the remaining 28 candidates, the mutation range was from 1 to 6 nucleotide mutations with an average mutation rate of 2.75 nucleotides per plyCA gene, which was in the 2-3 nucleotide mutation range we were targeting. At the translational level, this particular nucleotide mutation range and frequency yielded an amino acid mutation range of 1 to 5 amino acids, with an average mutation rate of 1.9 amino acid mutations per PlyCA polypeptide.
Of the 28 candidates with at least one amino acid mutation, four of these mutant constructs were randomly chosen for further characterization to validate that the extensive screening process of the directed evolution methodology was indeed functioning properly. The mutant PlyC enzymes were purified to > 95% homogeneity based on SDS-PAGE analysis as previously described6-7. Enzyme kinetics of WT PlyC and each of the four PlyC mutants were characterized at equal molar concentrations after incubating the purified enzymes at various elevated temperatures. Activity was monitored after the addition of D471 GAS by measuring the optical density at 600 nm every 15 sec for 20 min. Activity was defined as the residual maximum velocity of the enzyme after heat incubation. Of the four candidates randomly selected for further characterization, mutant 29C3 showed the most enhanced thermal behavior. The thermal behavior of WT PlyC and 29C3 was investigated at different temperatures while incubating and assaying for activity in PBS pH 7.2 at 80 nM and 40 nM concentrations, respectively. The 45-50 °C incubation experiments (Figure 3) were performed in a thermocycler where the samples were incubated in a thin-walled 96 well thermocycler plate in a total volume of 120 μl. The 35 °C, 40 °C and 45 °C incubation experiments (Figure 4 and 5) were performed in a heat block where the samples were incubated in a 1.5 ml microcentrifuge tube in a total volume of 1.3 ml. All experiments were performed in triplicates.
WT PlyC and 29C3 showed no significant difference in kinetic stability at room temperature (25 °C) as depicted in the first set of bars in Figure 3. However, after a short-term incubation for 30 min from temperatures ranging from 45-50 °C, the activity of 29C3 was substantially greater than the activity specific to the WT construct at each temperature point. For example, WT PlyC displayed a 44% loss in activity at 45.2 °C whereas 29C3 showed only a 2% loss in activity at the same temperature.
Long-term incubation studies comparing the residual activity of both WT PlyC and 29C3 were additionally performed at 35 °C and 40 °C involving the measurement of residual activity at 24 and 48 hr time points. At 35 °C, 29C3 displayed 41% and 176% higher activity than WT PlyC at 24 and 48 hr incubation time points, respectively (Figure 4a). At 40 °C, 29C3 displayed 28% and 107% higher activity than WT PlyC at 24 and 48 hr incubation time points, respectively (Figure 4b).
The residual activity of WT PlyC and 29C3 were also monitored every 20 min for a total of 3 hr at 45 °C. WT PlyC was only able to retain 21% activity after a 3 hr incubation at this temperature whereas 29C3 was able to retain 46% activity. The half-life (t1/2) for WT PlyC and 29C3 were 67 and 147 min, respectively, suggesting that 29C3 has a 2.2 fold increase in kinetic stability at 45 °C (Figure 5).
Total Round 1 Candidates | 35 |
Candidates with WT PlyCA Sequence | 7 |
Candidates with ≥1 Amino Acid Mutation | 28 |
— Average Nucleotide Mutation Rate (nt/plyCA) | 2.75 |
— Nucleotide Mutation Range (nt) | 1-6 |
— Average Amino Acid Mutation Rate (AA/PlyCA) | 1.9 |
— Amino Acid Mutation Range (AA) | 1-5 |
Table 1. Candidate Pool Genomic Analysis.
Figure 1. Directed evolution assay methodology. In directed evolution, one starts with the lead enzyme, which would be WT PlyC for the first round. A library of PlyC mutants containing random unbiased nucleotide mutations to the plyCA gene is then constructed by an error-prone DNA polymerase, cloned into the expression vector pBAD24 and transformed into the E. coli strain DH5α already containing pBAD33:plyCB. Individual transformants are inoculated into their own specific well of a 96 well microtiter plate. Through an extensive screening process, individual PlyC mutants that are catalytically active after incubation at a non-permissible temperature are classified as mutants with enhanced kinetic stability. Theconstruct displaying the most progressed thermal behavior will consequently becomes the lead enzyme for the next round of random mutagenesis. Overall, there are three complete rounds of random mutagenesis followed by DNA shuffling resulting in a bacteriolytic molecule with evolved thermal behavior. Click here to view larger figure.
Figure 2. 96 well microtiter plate templates for directed evolution assay. The microtiter plate schematic during the determination of the optimal heating conditions (Figure 2a) consists of inoculating a single parental WT PlyC clone into each row of the microtiter plate. The microtiter plate schematic during mutant screening (Figure 2b) consists of inoculating each well in column 1 with a unique parental WT PlyC clone as well as inoculating each well in columns 2-12 with a distinct PlyC mutant clone. Wells A1-D1 are designated for the positive parental controls, which consist of WT PlyC constructs not exposed to non-permissible temperature incubation. Wells E1-H1 are designated for the negative parental controls, which consist of WT PlyC constructs that are exposed to non-permissible temperature incubation.
Figure 3. Residual activity kinetic analysis comparing WT PlyC and 29C3. Enzymes were purified to homogeneity and incubated for 30 min in a thermocycler at equal molar concentrations at either room temperature or at a temperature gradient ranging from 45-50 °C. Enzymatic activity correlates with the maximum velocity displayed after specific temperature incubation. With the exception of 25 °C and 47.7 °C, the variation in activity between WT and 29C3 at each temperature correlated to a p-value < 0.05. All data are reported as the mean ±SEM of three independent experiments.
Figure 4. Residual activity kinetic analysis comparing WT PlyC to 29C3 at 35 °C and 40 °C. Equal molar concentrations of purified WT PlyC and 29C3 were incubated in a heat block at 35 °C (Figure 4a) or 40 °C (Figure 4b). The residual enzyme activity was measured at 24 and 48 hr time points. The activity displayed by each construct was normalized to the maximum velocity displayed at time point zero. The variation in activity between WT and 29C3 at each temperature and time point correlated to a p-value < 0.05. All data are reported as the mean ±SEM of three independent experiments.
Figure 5. Residual activity kinetic analysis comparing WT PlyC to 29C3 at 45 °C. Equal molar concentrations of purified WT PlyC and 29C3 were incubated in a heat block at 45°C and the residual enzyme activity was measured every 20 min for 3 hr. The activity displayed by each construct was normalized to the maximum velocity displayed at time point zero. With the exception of the data points at 20 min, the variation in activity between WT and 29C3 at each time point correlated to a p-value < 0.05. All data are reported as the mean ±SEM of three independent experiments.
This protocol, outlined in Figure 1, presents a 96 well microtiter plate methodology that allows one to utilize directed evolution to increase the thermostability of any bacteriolytic enzyme. Through the use of an error-prone DNA polymerase, one can introduce random mutations that increase the overall kinetic stability of the translated bacteriolytic molecule of interest, which is typically due to molecular reorganizations consisting of increasing electrostatic, disulphide bridge and hydrophobic interactions that generate improved molecular packing, enhanced modifications of surface charge networks or reinforcement of a higher oligomerization state17-18. After the introduction of nucleotide mutations, an extensive screening procedure was then used to identify mutations that are thermally beneficial. The representative results presented here were based on one round of random mutagenesis. In reality, it has been suggested that at least three successive rounds of random mutagenesis, each using the most robust mutant as the lead enzyme, followed by DNA shuffling is appropriate for these types of directed evolution thermal behavior studies2.
It is important to understand that while this protocol identifies mutations that augment kinetic stability, these variants do not necessarily have increased thermostability. A molecule is considered thermostable when it has the ability to retain its structure at a high temperature. However, in the assay presented, the residual activity of the mutant lysates are not assayed at the same temperature that they were initially incubated at during the heat treatment step and thus one may be selecting for mutations that promote refolding rather than enhanced thermostability19. While we are extrapolating thermal behavior as an indication of thermal stability, true thermal stability must be measured empirically by biophysical methods such as circular dichroism (CD), differential scanning calorimetry (DSC) and differential scanning fluorimetry (DSF). Preliminary data from CD and DSF experiments suggest that several candidates identified from the methodology presented do indeed display progressed thermostability (data not shown).
Although the methodology presented here is specific to implementing increased thermal behavior to PlyC, this same methodology can additionally be employed to enhance thermal activity to other bacteriolytic enzymes with some minor alterations to variables such as buffers, mutation rates, heating conditions and expression systems. One critical variable associated with this assay relates to the nucleotide mutation rate of the error-prone DNA polymerase. We chose to use the error-prone Mutazyme II DNA polymerase in our directed evolution assay due to experimental evidence suggesting this particular enzyme displays the least amount of bias with regards to which nucleotides are randomly incorporated into the gene of interest when compared to other random mutagenesis techniques such as use of hydroxylamine, mutator E. coli strains and other error-prone DNA polyermases20. In general, lower mutation rates tend to be more desirable for two reasons. First, engineering thermostability to proteins typically consists of relatively few amino acid substitutions. High mutation rates can cause dramatic structural alterations to particular molecules that can conclusively result in structural misfolding and functional discrepancies. For example, the incorporation of glycine and proline residues can disrupt alpha helical secondary structures21. Second, high nucleotide mutation rates increase the chances of incorporating premature stop codons, resulting in truncated molecules that are biologically inactive. In general, lower mutation rates are preferred because lower error rates result in the accumulation of adaptive mutations while higher mutation rates generate neutral or deleterious mutations22.
For each round of random mutagenesis, another critical variable that one must optimize involves the incubation temperature used during the screening process. Selection of the experimental non-permissive temperature is relatively challenging. We initially used an incubation temperature that was 10 °C higher than the non-permissive temperature of PlyC, however after screening 3,000 mutants, we were unable to identify any advantageous mutations. Keeping in mind directed evolution methodologies consist of sequentially identifying beneficial amino acid substitutions that have additive effects, it was determined that a less stringent temperature should be used during the screening process even if it increased the chance of generating false-positives. As such, we used an incubation temperature that was only a few degrees above the non-permissive temperature and rescreened selected candidates to rule out the false-positives.
We decided to utilize an incubation temperature that was not too temperate (≤ 1 °C higher than lowest non-permissive temperature) and conversely not too harsh (≥ 5 °C higher than lowest non-permissive temperature). The ideal temperature we decided to use in this particular assay was a moderate 2 °C higher than the experimentally determined non-permissive temperature. Choosing a temperature that is too mild will result in a much higher false positive identification rate than the ~20% we observed (Table I) when using a moderate incubation temperature. Furthermore, using an incubation temperature that is too harsh could ultimately result in the inability to identify mutants with significantly enhanced thermal behavior. For example, using an incubation temperature of ≥ 5 °C would have resulted in the inability to identify the promising 29C3 mutant as well as the majority of the other candidates selected during the first round of screening.
In summation, we provide a protocol for engineering bacteriolytic enzymes to acquire enhanced kinetic stability. Moreover, this same assay, without the heating steps, can be used to screen for mutations that increase catalytic activity. Augmenting both catalytic activity and thermostability is an important developmental hurdle facing any therapeutic enzyme. While the details of this protocol are specific to the endolysin PlyC, the methodology can be adapted to any bacteriolytic enzyme with only a few alterations. We presented preliminary results from only the first round of mutagenesis that validates that functionality of the assay, resulting in the implementation and identification of mutations that generate an increase in molecular thermostability of significant magnitude.
The authors have nothing to disclose.
The authors wish to thank Emilija Renke and Janet Yu for technical assistance. DCN is supported by grants from the United States Department of Defense (DR080205, DM102823, OR09055, and OR090059). RDH is supported by a grant from the Maryland Agriculture Experiment Station and an NIH Training Grant in Cell and Molecular Biology.
Name of the reagent | Company | Catalogue number |
B-PER II Bacterial Protein Extraction Reagent | Thermo Scientific | 78260 |
GeneMorph II Random Mutagenesis Kit | Agilent Technologies | 200500 |
Clear 96 Well, Flat-bottomed Microtiter Plate With Lid | BD Vacutainer Labware Medical | 353072 |
96 Well Nonskirted PCR Plates | Fisher Scientific | 14-230-232 |
Swing-bucket Rotor | Eppendorf | A-2-DWP |
Ampicillin Sodium Salt | Fisher Scientific | BP1760 |
Chloramphenicol | Fisher Scientific | BP904 |
Arabinose | Fisher Scientific | BP2504 |
pBAD24 Expression Vector | ATCC | 87399 |
pBAD33 Expression Vector | ATCC | 87402 |