Characterizing structural features of cuticle-degrading proteases from fungi by molecular modeling
© Liu et al; licensee BioMed Central Ltd. 2007
Received: 02 January 2007
Accepted: 18 May 2007
Published: 18 May 2007
Serine proteases secreted by nematode and insect pathogenic fungi are bio-control agents which have commercial potential for developing into effective bio-pesticides. A thorough understanding of the structural and functional features of these proteases would significantly assist with targeting the design of efficient bio-control agents.
Structural models of serine proteases PR1 from entomophagous fungus, Ver112 and VCP1 from nematophagous fungi, have been modeled using the homology modeling technique based on the crystal coordinate of the proteinase K. In combination with multiple sequence alignment, these models suggest one similar calcium-binding site and two common disulfide bridges in the three cuticle-degrading enzymes. In addition, the predicted models of the three cuticle-degrading enzymes present an essentially identical backbone topology and similar geometric properties with the exception of a limited number of sites exhibiting relatively large local conformational differences only in some surface loops and the N-, C termini. However, they differ from each other in the electrostatic surface potential, in hydrophobicity and size of the S4 substrate-binding pocket, and in the number and distribution of hydrogen bonds and salt bridges within regions that are part of or in close proximity to the S2-loop.
These differences likely lead to variations in substrate specificity and catalytic efficiency among the three enzymes. Amino acid polymorphisms in cuticle-degrading enzymes were discussed with respect to functional effects and host preference. It is hoped that these structural models would provide a further basis for exploitation of these serine proteases from pathogenic fungi as effective bio-control agents.
Proteases are enzymes that catalyze the hydrolysis of peptide bonds in other proteins and are often used in biotechnology, industry, and agriculture as bio-control agents against parasites. Serine proteases (EC 3.4.21.-) are present in virtually all organisms and exist as two major families, the trypsin-like (EC 220.127.116.11) and the subtilisin-like (EC 18.104.22.168) families. These two families have independently evolved with a similar catalytic mechanism which has been widely investigated [1–4]. Although the overall fold of various serine proteases may differ, they all follow the same mechanism of action through an identical stereochemistry of the catalytic triad and oxyanion hole. In this mechanism, the Ser functions as the primary nucleophile and the His plays a dual role as the proton donor and acceptor at different steps in the reaction. The role of Asp is believed to bring the His into the correct orientation to facilitate the nucleophilic attack by the Ser. The role of the oxyanion hole is to stabilize the developing negative charge on the oxygen atom of the substrate during the formation of the tetrahedral intermediate [5–7].
Serine proteases secreted by nematophagous and entomophagous fungi are one of the important virulence determinants in the infection processes of invertebrate mycopathogens [8, 9]. They function mainly to degrade insect and nematode cuticles. These cuticles comprise primarily chitin fibrils embedded in a protein matrix, in which the quantity and type of proteins vary between species, tissues and growth stages of insects and nematodes [10, 11]. During the infection process, an initial mechanical pressure may play a role but of major importance is the subsequent enzymatic degradation of cuticles within the immediate vicinity of the infection hyphae. Being accompanied by the enzymatic degradation, the cuticles are penetrated and the prey is invaded and eventually digested by fungi. Associated with the ability to penetrate the complex and variable cuticle structures, the pathogenic fungi produce a number of different cuticle-degrading enzymes such as lipases, chitinases and at least four different classes of proteases [12, 13]. Among these proteases the most extensively studied insect cuticle-degrading protease is a serine protease, PR1, produced by an entomopathogenic fungus Metarhizium anisopliae [14, 15]. PR1 is a crucial factor in the infection of the insect because it is produced more rapidly and is in a higher concentration than the other cuticle-degrading enzymes . Furthermore, it has been observed, using gold-labeled antibodies, that the PR1 is secreted by the Metarhizium anisopliae during penetration of the host cuticle . A similar protease, VCP1 from the nematophagous fungus, Pochonia chlamydosporia (Syn. Verticillium chlamydosporium), has also been observed to play an important role during the infection process of nematode eggs [18, 19]. More recently, an alkaline serine protease, Ver112 from the nematophagous fungus Lecanicillium psalliotae exhibiting cuticle-degrading and nematicidal activity was cloned and sequenced by our laboratory [9, 20].
Cloning and sequencing of genes encoding the three proteases (PR1, VCP1 and Ver112) have shown that they have a high degree of sequence similarity with each other and belong to the proteinase K family of subtilisin-like serine proteases (subtilases) which is a large family of endopeptidases found only in fungi and Gram-negative bacteria . These enzymes show conservation of the Asp-His-Ser catalytic triad and catalytic domain. However, the three-dimensional (3D) structures have not been resolved for any of the cuticle-degrading proteases so far by either X-ray crystallography or nuclear magnetic resonance (NMR) techniques. A thorough understanding of the 3D structures of these enzymes can provide greater insight into the mechanism and specificity of enzyme catalysis which would guide future targeted mutagenesis for improved efficiencies.
In this paper, the 3D structural models of the mature enzymes for PR1, VCP1, and Ver112 were modeled based on the available high resolution crystal structure (at 0.98 Å) of Proteinase K using the homology modeling technique. Based on both the multiple sequence alignment and structural models, predictions were made for these cuticle-degrading enzymes concerning calcium-binding site, disulfide bridge, salt bridge, hydrogen bond, aromatic interaction, substrate-binding site and electrostatic surface potential. The physiochemical and geometrical properties of the three cuticle-degrading enzymes and the template protease K were compared to identity factors that might contribute to substrate specificity and host preference by these fungi.
Comparison of physiochemical properties
Physiochemical properties and optimal reaction conditions of the proteinase K (1IC6) and the three cuticle-degrading enzymes (PR1, VCP1 and Ver112).
T j (°C)
Model description and structural comparison
Pairwise least square superpositions of the four structural models give backbone root mean deviation (RMSD) values of 0.04–0.75 Å. The structural resemblance with regard to the RMSD value is in the following order: PR1-VCP1(0.04 Å) > Ver112-PR1 (0.34 Å) > VCP1-Ver112 (0.35 Å) > Ver112-1IC6 (0.58 Å) > 1IC6-PR1 (0.75 Å) = 1IC6-VCP1 (0.75 Å), indicating that the models of the cuticle-degrading enzymes are more similar to each other than to 1IC6. This roughly reflects the amino acid sequence identities between them. The sequence identities of PR1-VCP1, Ver112-PR1, VCP1-Ver112, Ver112-1IC6, 1IC6-PR1, and 1IC6-VCP1 are 78.6%, 69.5%, 64.3%, 64.1%, 66.4%, and 62.5%, respectively.
Geometrical properties of the proteinase K (1IC6) and the three cuticle-degrading enzymes (PR1, VCP1 and Ver112).
Prediction of Ca2+ binding sites
Multiple sequence alignment and structural models suggest that there are two disulfide bridges occupying the same positions in each of the four proteases. One of the disulfide bridges, Cys34–Cys123, connects β1 carrying the catalytic residue Asp39 and a loop connecting α3 and β4. Although this loop is far away from the catalytic triad residues and substrate-binding regions, it links α3 located next to the substrate-binding segment 100–104 and β4 preceding the substrate-binding segment 132–136. As a result, this loop may function as a "hinge" responsible for modulating the orientation of the substrate-binding segments. Thus the Cys34–Cys123 can be thought to get involved in stabilizing the regions close to the catalytic residue Asp39 and substrate-binding sites. In addition, an absolutely conserved salt bridge, Asp117:Arg121, is also found to participate in stabilizing the "hinge" loop. The other disulfide bridge, Cys178–Cys249, connecting β8 located next to the S1 pocket and α6 situated close to the C-terminus, makes the peripheral α6 anchor to the rest of the molecule and contributes to the stability of the C-terminal region. In addition, Cys178–Cys249 may also influence the orientation of the S1 residues via a loop comprising residues 171–177.
Hydrogen bonds and salt bridges
The numbers of hydrogen bonds/salt bridges in 1IC6, PR1, VCP1 and Ver112 are 223/14, 189/12, 189/12, and 195/13, respectively (Table 2). Hydrogen bonds play an important role in the overall stability of the protein due to their large number and wide distribution. In contrast, the contribution of salt bridges to protein stability is more localized . In this paper, we mainly focus on the salt bridges that contribute to the stability of local structures, and on the number and type of hydrogen bonds and salt bridges that are close to or within the substrate-binding sites.
A shared salt bridge network Arg12:Asp187:Lys/Arg18 is found in the three cuticle-degrading enzymes to participate in stabilizing the N-terminal region which is stabilized by Ca2 in 1IC6. Another common salt bridge network, Arg52:Glu50:Lyr/Arg80, is found in all four enzymes, connecting and stabilizing the surface loops between a 3/10 helix (residues 47–49) and β2 and between α2 and β3. As described above, the absolutely conserved salt bridge, Asp117:Arg121, together with the disulfide bridge Cys34–Cys123, anchors and stabilizes the "hinge" loop. In 1IC6, PR1 and VCP1, there are two successive negatively charged residues, Asp26 and Glu/Asp27, in a loop connecting α1 and β1, and these two residues can form at least two salt bridges with two successive positively charged residues, Lys86 and Lys87 in a loop connecting α2 and β3. In contrast, Ver112 contains just one salt bridge at the equivalent positions due to the replacement of the negatively charged residue at position 27 with a neutral residue, Thr. In addition, a salt bridge Asp112:Arg147 is found in both 1IC6 and Ver112 to bridge α3 and α4, which are directly connected to the substrate-binding segments comprising residues 100–104 and 132–136. Although the negatively charged residues (Asp/Glu112) are conserved among the four proteases, the replacement of Arg147 with Ala in both PR1 and VCP1 disrupts the formation of this salt bridge. In PR1 and VCP1, a salt bridge network, Lys250:Asp254:Lys251, is found to participates in stabilizing the C-terminal helix α6. However, only one salt bridge, Arg250:Asp254, is found in α6 of 1IC6, and no salt bridge is found in α6 of Ver112. Interestingly, a common disulfide bond Cys178–Cys249 is found to participate in stabilizing the α6 in all four investigated protease models.
There are at least one positively charged residue, Arg, and one negatively charged residue, Asp, in a region comprising residues 162–169 which is part of, or in close vicinity to the S1 pocket. In this region, 1IC6 contains a salt bridge, Asp165:Arg167, PR1 and Ver112 share a salt bridge, Asp162:Arg164, and VCP1 contains a salt bridge network comprising three residue, Asp162:Arg168:Asp165. The number of hydrogen bonds in this region is 8, 6, 5, and 4 in 1IC6, PR1, Ver112, and VCP1, respectively. These salt bridges and hydrogen bonds contribute to the structural stability of this region and may cause the P1 residue of the substrate to adopt the correct orientation to facilitate nucleophilic attack by the Ser. Interestingly, the observed reduction in the hydrogen bond number within this region for VCP1 likely is compensated by the salt bridge networks involving the three residues, Asp162:Arg168:Asp165.
The charged residue clusters are found in a loop comprising residues 59–68 and in a region comprising residues 94–98. The loop 59–68 preceding the catalytic residue His69 contains abundant polar residues and resides on the molecular surface, hereafter it is referred to as the polar surface loop (PSL). The loop 95–101 is part of, or in close proximity to the S2 site, and hereafter it is referred to as the S2-loop. In both PR1 and VCP1, there is one salt bridge network, Asp65:Lyr94:Asp97, that bridges the PSL and the S2-loop together and hence stabilizes them. In 1IC6 and Ver112, three salt bridge networks, Asp65:Lyr94:Asp97, Asp65:Lyr94:Asp98 and Arg64:Asp98:Lyr94 are responsible for stabilizing these two loops. Interestingly, an absolutely conserved hydrogen bond network Asp97-[Gly100, S101] which participates in stabilizing the S2-loop is found in all the four enzymes. In addition to Asp65:Lyr94 and Arg64:Asp98, another salt bridge involved in stabilizing the PSL in 1IC6 and Ver112 is Glu43:Arg64. This salt bridge is not found in the PR1 and VCP1 due to the replacement of the Arg64 by a neutral residue, Asn or Thr. Furthermore, there is an absolutely conserved hydrogen network Asp65-[Gly68, Thr71] involved in stabilizing the PSL in the four enzymes, and the number of hydrogen bonds participating in stabilizing the PSL is 10, 9, 7 and 8 in 1IC6, PR1, VCP1 and Ver112, respectively. Taken together, it is thought that the abundant hydrogen bonds in combination with a small quantity of salt bridges provide determinants for the stabilization of the PSL. It has also been observed that the PSL can be stabilized by calcium in two proteinase K-like enzymes, the VPRK from Vibrio sp. PA44 [PDB:1SH7]  and the thermitase from Thermoactinomyces vulgaris [PDB:1THM] .
Aromatic residues and interactions
Clustered aromatic residues with potential interactions between them assist in enhancing the thermal stability of the protein . Aromatic residues that are conserved among the four enzymes are at positions 8(Trp), 23(Tyr), 25(Tyr), 36(Tyr), 49(Phe), 59(Tyr/Phe), 82(Tyr/Trp), 113(Phe/Tyr), 137(Tyr/Phe), 192(Phe), 195(Tyr/Phe), 202(Phe), 212(Trp), 236(Tyr), 274(Tyr/Phe), and 277(Tyr/Phe). Out of them two absolutely conserved aromatic interactions, Tyr23–Tyr25 and Tyr/Phe59-Phe/Tyr113, are found in all four protease models. The aromatic ring stack of Tyr23–Tyr25 can contribute to the stabilization of the N-terminal region, and Tyr/Phe59-Phe/Tyr113 can contribute to the stabilizations of the PSL and α3. The α3 is in close proximity to the substrate-binding segment 100–104. A shared aromatic interaction, Tyr36-Phe/Tyr91, observed in 1IC6, PR1 and Ver112 can potentially contribute to the stabilizations of both β1 and β3. However, this aromatic interaction does not occur in VCP1 due to the replacement of the aromatic residue Phe or Tyr at position 91 by Leu in VCP1. Another aromatic interaction Tyr60-Tyr61 involved in stabilizing PSL is found only in 1IC6.
Substrate-binding region and electrostatic surface potential
The S2' site is a hydrophobic pocket formed primarily by residues 192, 221 and 222. The S1 is a distinct, large and elongated pocket that is primarily formed at the side and bottom by segments 132–135 and 158–161, respectively. In addition, the S1 pocket is also surrounded at the rim by residue 162, at the bottom end by 169, and at the top by segment 222-225 carrying the nucleophilic residue Ser224, respectively. The S2 is a less distinct and smaller pocket than S1. It is bounded at one side by residue 100 and at the other side by the active site residue His69, at the bottom by both the hydrophobic residue 96 and the active site residue Asp39, at the bottom end by residue 40, and at the rim by residue 67, respectively. S3 is not a pocket because the side chain of the residue (Ser101 in all four enzymes) at the S3 site points away from the substrate-binding region/surface, leading to the formation of a relatively protrudent patch on the protein surface (Figure 5). The other residue that may also get involved in the interaction with the P3 side chain is at a close position 100. S4 is a very distinct pocket that is situated between two segments, 100–104 and 132–136. This pocket appears to have two subsites, S4a and S4b. S4a has residues 96 and 107 at the bottom, 133 at the side, and 102 at the rim, respectively. S4b has residues 104 at the side, 141 at the bottom, and 135 and 136 at the rim, respectively. The side chains of these residues determine the size of the S4 pocket.
Among the four investigated enzymes, the amino acid sequences in the substrate-binding sites of S2' and S3 are absolutely conserved but of S1, S2 and S4 are relatively conserved. Residue changes in the S1 pocket are found at positions 162 and 169 that are located at the rim and bottom end of this pocket, respectively. The neutral Asn162 in 1IC6 is replaced by the acidic Asp in the three fungal cuticle-degrading enzymes PR1, VCP1, and Ver112. Although Asn162 and Asp162 are practically identical in conformation, the acidic residue Asp162 would contribute to the more anionic characteristics of the S1 site in the cuticle-degrading enzymes (Figure 5), which in turn may shift the specificity to the basic P1 residues [3, 4]. In addition, it has been observed that the anionic character can increase flexibility of an enzyme , and in particular increase flexibility in the active site region . We therefore suggest that Asp162 in the cuticle-degrading enzymes would increase the flexibility of the S1 site, which in turn may enhance the catalytic efficiency of the cuticle-degrading enzymes and infection virulence by the fungi. Position 169 is occupied by the Tyr in 1IC6 and VCP1 but the Thr in PR1 and Ver112. Mutational studies on subtilisin BPN' have shown that P1 specificity can be dramatically modulated by substitution of residue at this position [29, 30]. Residue changes in the S2 site are only found at position 67 with replacement of the Asn in 1IC6 by the His in PR1, VCP1 and Ver112, implying a preference of the cuticle-degrading enzymes for the negatively charged P2 residues at a low pH. The most obvious differences observed within the substrate-binding sites are located in the S4 pocket at positions 103, 104, 107, 136, and 141, in which all four enzymes contain different amino acids. For example, the Ser103 in the three cuticle-degrading enzymes in place of the Gln103 in 1IC6 would enlarge the size of the S4 pocket due to its relatively smaller side chain than that of the Gln. The Tyr104 in PR1, VCP1 and 1IC6 is observed to be replaced by the Leu in Ver112. Since previous studies have shown that the Tyr104 forms a flexible lid to the S4 pocket [3, 4] of subtilases, we speculate that the Leu104 in Ver112 would broaden the entrance for P4 residues and simultaneously increase the hydrophobic property of the S4 pocket. The Ile107 in 1IC6, PR1 and Ver112 is replaced by the Val in VCP1 and the Val141 in 1IC6, PR1 and VCP1 is replaced by the Leu in Ver112. These two substitutions likely have a minor influence on the size and hydrophobicity of the S4 pocket as Ile/Val107 and Val/Leu141 are largely buried at the S4 pocket bottom and are all hydrophobic residues. The substitution of Pro136 in VCP1 for Gly136 in 1IC6, PR1 and Ver112 results in a larger rim of the S4 pocket in VCP1 than in the other enzymes. In addition, the Pro136 is located at the other side of the Tyr104 with which a gate is formed at the entrance of the S4 pocket, as thus we speculate that the Pro136 substitution may influence the entry of the P4 residues with large side chains. The hydropathy score reveals that the hydrophobic values of the S4 pockets in 1IC6, PR1, VCP1 and Ver112 are 26.0, 28.6, 27.1 and 33.1, respectively, indicating that the S4 pockets in cuticle-degrading enzymes are more hydrophobic than in 1IC6.
A notable difference in the surfaces of the proteins compared here is their different electrostatic surface potentials (Figure 5). Large parts of the surface of PR1 and Ver112 are positively charged, reflecting that these two proteases contain more positively charged residues than 1IC6 and VCP1 (Table 1). Although 1IC6 and VCP1 contain the same number of positive/negative residues (Table 1), the negatively charged surfaces are relatively evenly dispersed on the positively charged or less charged surfaces in 1IC6 (Figure 5A), whereas in VCP1 (Figure 5C), the negatively charged surfaces are mainly concentrated on regions that are part of or in close vicinity to the substrate-binding sites, and the positively charged surfaces are concentrated on the other regions.
The 3D structural models of the three cuticle-degrading enzymes, PR1, VCP1 and Ver112, have been modeled using the homology modeling technique. According to Baker and Sali , high-accuracy homology models with about 1.0 Å RMS error for the main-chain atoms can be obtained given more than 50% sequence identity to their templates, and the calculated models are comparable to the accuracy of a medium-resolution NMR or X-ray structure. Since the sequence identities of the three cuticle-degrading enzymes to the template 1IC6 are all greater than 60% and there is no large insertion or deletion relative to 1IC6, the modeled structures should be fairly accurate and can be used for predicting the Ca2+ binding site, disulfide bridge, hydrogen bond, salt bridge, substrate specificity and catalytic activity. Model assessments using the program PROCHECK  reveal that for structural models of PR1, VCP1, and Ver112, the residues located in the most favored regions of the Ramachandran plot are 85.8%, 86.2%, and 86.0%, and in the additionally allowed regions are 13.8%, 12.9%, and 13.6%, respectively. Only one residue, the Arg126 in Ver112 is located in the disallowed region. Interestingly, this Arg is located in a loop region (the "hinge" loop) between α3 and β4, and it can be replaced by Gly in the isoform of Ver112, resulting in an amino acid polymorphism at position 126 . The amino acid polymorphisms will be discussed below.
Although the total quantities of hydrogen bonds and salt bridges are comparable among the three cuticle-degrading enzymes, the differences in their number, type, and distribution in regions that are part of or in close vicinity to the substrate-binding sites can have a significant effect on the local flexibility, which will influence the substrate binding, catalytic activity and substrate specificity of these enzymes. Helland et al.  has argued that a strong hydrogen bond network, Asn97-[S99, S101], is able to bring the S2-loop of a proteinase K-like enzyme, SPRK from the psychrotroph Serratia species [PDB:2B6N], into a 'hub and spokes' arrangement, which confers more rigidity to the S2-loop and reduces the binding affinity (higher KM value) of SPRK towards the synthetic substrate suc-Ala-Ala-Pro-Phe-nitroanilide. In addition, the tight S2-loop is also able to influence the substrate specificity profile of SPRK . With regard to the three cuticle-degrading enzymes compared here, a conserved hydrogen bond network, Asp97-[Gly100, S101], and a conserved salt bridge network, Asp65:Lyr94:Asp97, are found in the S2-loops. However, the number of hydrogen bonds involved in stabilizing the S2-loop is 4, 4 and 6 in PR1, VCP1 and Ver112, respectively, and the number of salt bridges involved in stabilizing the S2-loop is 1, 1, and 3 in PR1, VCP1, and Ver112, respectively. Therefore the S2-loop appears to be more rigid in Ver112 than in the other two enzymes, and it is not unreasonable to believe that the relatively high rigidity of the S2-loop in Ver112 would lead to the reduction in substrate affinity, which in turn would lower the catalytic efficiency for Ver112 in comparison to PR1 and VCP1.
In terms of substrate specificity, we have identified variable residues in the substrate-binding sites because any variation in specificity of subtilases is most likely to be caused by these variable residues, especially those residues whose side chains interact directly with the P1 and P4 residues of the substrate [3, 4]. As previously mentioned, the substitution of Asp162 in the three cuticle-degrading enzymes for Asn in 1IC6 would not only enhance the catalytic efficiency but also shift the specificity to the basic P1 residues. The position 169 at the bottom end of the S1 pocket is occupied by either the Tyr or Thr in the four enzymes. Although Tyr is larger than Thr, the S1 pocket in 1IC6 and VCP1 remains as wide as that in PR1 and Ver112 as the side chains of these two residues rotate away from the S1 pocket bottom end (Figure 5). In addition, docking data has shown that this Tyr in 1IC6 is not likely to interact directly with the P1 residue of a substrate . For such interactions to occur, the side chain will have to rotate, and hence, cause steric clashes with the enzyme main chain that will require energy to resolve. Taken together, we suggest that the residue substitution observed at position 169 may have a minor effect on substrate specificity for the three cuticle-degrading enzymes. On the contrary, the S4 site may contribute significantly to the substrate specificity because of the five observed variable residues in the S4 pocket. Briefly, the cuticle-degrading enzymes have larger and more hydrophobic S4 pockets than that of the 1IC6. The difference in size is mainly caused by the substitution of the smaller Ser in the cuticle-degrading enzymes for Gln in 1IC6 at position 103. In addition, Ver112 may have a preference for larger and more hydrophobic P4 residues than the other two cuticle-degrading enzymes due to the Leu substitution for Tyr at position 104, and VCP1 may have a preference for relatively smaller P4 residues due to the Pro substitution for Gly at position 136. Furthermore, the differences in conformational flexibility of the substrate-binding regions may also contribute to substrate specificity due to the induced fit process of the substrate binding . It is therefore reasonable to believe that the more rigid S2-loop in Ver112 than in PR1 and VCP1 might cause variation in the specificity profile. It should be pointed out that, although the possible differences in substrate specificity between the three cuticle-degrading enzymes have been predicted from our models, more kinetic data with substrates which mimic the real "target" of the proteolytic attack should be obtained in order to verity these predictions.
The common electrostatic surface feature shared by the three cuticle-degrading enzymes is that their substrate-binding pockets are largely negatively charged while most of the other surfaces are positively charged. As mentioned above, the anionic feature of the substrate-binding regions would increase the local conformational flexibility of them, which in turn would enhance the catalytic efficiency. In addition, the large patches of the positively charged surfaces increase the adsorption of the cuticle-degrading enzymes to cuticles carrying abundant acidic residues [19, 35], and also explain why the cuticle-degrading enzymes studied here are all alkaline proteases.
Some isoforms of PR1, VCP1 and Ver112 carrying the amino acid polymorphism sites have been cloned and sequenced. For example, up to 18 residue variations within PR1 of several isolates from different hosts have been detected . In addition, the number of the amino acid polymorphisms detected for VCP1 and Ver112 is 2 and 1, respectively. However, no polymorphism has been found at the catalytic residue triad of the Asp, His and Ser, indicating a strong selective pressure for the maintenance of the structure/function of this site. Some of the polymorphisms observed near the substrate-binding regions may have effects on the substrate utilization, kinetics or surface property of the enzymes . It has been suggested that the abundant polymorphisms within PR1 may render its isoforms recalcitrant to some of the antiproteolytic hemolymphal peptides and could enable the pathogenic fungus M. anisopliae to exploit a wider range of substrates . For VCP1, its two polymorphism sites are observed at positions 62 and 99 . The polymorphism at position 62 is the replacement of an acidic residue, Glu in the root knot nematode isolates by a neutral residue, Gln in the cyst nematode isolates. This substitution is far away from the substrate-binding regions and may have a negligible influence on the local conformation because the Glu is similar in size to Gln. However, this substitution reduces the negative charge of the VCP1 isoform which in turn enhances enzyme adsorption to the acidic cuticles . The other polymorphism at position 99 is the replacement of a Gly in the root knot nematode isolates by an Ala in the cyst nematode isolates. This replacement is situated in the S2-loop and we speculate that the smaller Gly in VCP1 isoforms from the root knot nematode isolates would increase flexibility of the S2-loop and reduce steric hindrance of the larger P2 and P3 residues with corresponding substrate-binding sites, the S2 and S3. Indeed, substrate utilization assays  reveal that VCP1 isoforms from the root knot nematode isolates have a higher activity towards Suc-Ala-Val-Pro-Phe-pNA than towards Suc-Ala-Ala-Pro-Phe-pNA. As mentioned above, only one polymorphism at position 126 has been found so far in Ver112. This polymorphism is the replacement of the Arg126 in isolates L. psalliotae 112 and 608 by the Gly in isolate L. psalliotae 602 . On the one hand, Arg126 increases the positive charge of VCP1, which in turn enhances the adsorption of the enzyme to the negatively charged nematode cuticle. On the other hand, Arg126 is located in the 'hinge' loop that is responsible for modulation of the orientation of the substrate-binding segments 100–104 and 132–136. Although the "hinge" loop is stabilized by the disulfide bridge Cys34–Cys123 and the salt bridge Asp117:Arg121, Arg126 may make it unstable due to the large size of Arg and no apparent negatively charged complementary residue being found to form the salt bridge. We speculate, therefore, that the Arg126 would make the "hinge" loop more flexible and influence the orientation of the substrate-binding segments (100–104 and 132–136) and the catalytic activity. However, more kinetic data should be obtained to clarify these effects. In addition, the host related amino acid polymorphism in PR1, VCP1 and Ver112 between different isolates may also contribute to host preference by fungi.
The 3D structural models of the three cuticle-degrading proteases, PR1, VCP1 and Ver112, have been modeled using the homology modeling technique. These models reported here are very similar to each other and to the template structure of 1IC6, with restricted conformational differences found only in certain surface loops and the N-, C-termini. These structures share the common disulfide bridges and similar calcium binding sites, physiochemical properties, optimal reaction conditions and high sequence identities. These similarities support the initial hypotheses that the three cuticle-degrading enzymes belong to the proteinase K-like serine protease family. These enzymes secreted by nematophagous and entomophagous fungi can be considered to be general-purpose proteases because they are required either for degrading cuticles or for growth on purified proteinaceous cuticles. However, the variable residues observed within the substrate-binding sites and the different conformational flexibilities in the S2-loop can be expected to cause variation in substrate specificity as well as in catalytic activity. The host related amino acid polymorphisms observed for these cuticle-degrading enzymes are not only of great interest to the functional pathology of fungi, but also to future exploitation of these fungi as biocontrol agents.
The protein sequences for PR1 , VCP1  and Ver112  were obtained from the GenBank database with accession numbers [GenBank:S22387], [GenBank:AJ427460] and [GenBank:AY692148], respectively. The signal peptide sequences and the pro-peptide sequences were removed and the obtained mature proteins consisted of 281, 281, and 280 amino acids for the PR1, VCP1 and Ver112, respectively. The X-ray structure of the serine protease proteinase K  ([PDB:1IC6], hereafter it is referred to as 1IC6 according to its PDB accession number) from Tritirachium album limber at 0.98 Å resolution was used as a template to model the 3D structural models of PR1, VCP1 and Ver112. All modelings were performed using the MODELLER 7v7 software package  running under the Linux operating system. Briefly, to obtain the structural model of PR1, the target sequence was first aligned with the template sequence of 1IC6 using an in-house written script "alignment" for the MODELLER 7v7. Second, 20 PR1 models were generated from the template 1IC6 using an in-house written script "get-model". The insertions and deletions were generated and optimized using a loop modeling sub-routine within the "get-model" script with a 3D_INTERPOLATION algorithm and a thorough optimization protocol (a thorough molecular dynamics/simulated annealing procedure and the MD_LEVEL is set to refine_3). Finally, these 20 models were clustered and an optimized average model was generated using a MODELLER 7v7 script "cluster". The transferred coordinates for a given target atom were the average of the 20 models in the largest cluster of the corresponding atom with the cluster cut off at 1.5 Å. The modeling procedures for VCP1 and Ver112 were identical to that for the PR1. To further optimize their structures, the four proteases (including protease K, 1IC6) were respectively subjected to the steepest descent energy minimization with the tolerance 100 kJ mol-1 nm-1, followed by a second conjugate gradient energy minimization with the tolerance 10 kJ mol-1 nm-1 using the program "mdrun" within the GROMACS software package . For PDB files of the homology models after energy minimizations, [seeAdditional file 1: file1_PR1.pdb], [see Additional file 2: file2_VCP1.pdb] and [see Additional file 3: file3_Ver112.pdb]. The overall geometrical qualities of these models were assessed using the program PROCHECK .
Physicochemical property calculations
The molecular weight, amino acid composition and charge distribution of each mature protease were calculated using the program SAPS . The isoelectric point (pI) and the extinction coefficient (EXTCOEF) at 280 nm of the mature enzymes were calculated using the pI service  and the EXTCOEF  program, respectively.
Structural analyses and comparisons
The geometrical properties of the mature proteases such as the radius of gyration (Rg), number of native contacts (NNC), number of hydrogen bonds (NHB), and number of salt bridges (NSB) were calculated using the programs g_gyrate, g_mindist, g_hbond and g_saltbr within the GROMACS software package , respectively. Pairwise superposition of structures was performed with the g_confrms program within GROMACS. Secondary structure element and solvent accessible surface area (SASA) of the models were analyzed with the DSSP program . The electrostatic surface potential was calculated using the Swiss-pdbViewer software package . The hydrogen network was represented by residue1-[residue2, residue3, residue4...], where residue1 could form hydrogen bonds with residue2, residue3 and residue4, etc. The salt bridge network was represented by residue1:residue2:residue3, where residue2 could form salt bridges with residue1 and residue3.
Amino acid sequence alignment
Multiple sequence alignment was performed using the CLUSTALW program  with the gap open and extension penalties 10.0 and 0.20, respectively. Next, the obtained multiple sequence alignment was improved by manual adjustments according to the known alignments from MODELLER 7v7 and the superimpositions of 3D structural models. The amino acid numbering used throughout this paper corresponds to that of 1IC6. Substrate-binding site residues in the S2', S1, S2, S3, S4 sites of the cuticle-degrading enzymes were assigned based on those of 1IC6 and the crystal structure of proteinase K in complex with an inhibitor Methoxysuccinyl-Ala-Ala-Pro-Ala-chloromethyl ketone [PDB:3PRK] . The nomenclature of the substrate residues and the substrate-binding sites of enzymes is according to Schechter and Berger . The scissile peptide bond is located between the P1 and P1' substrate residues and residues of the substrate are numbered P2', P1', P1, P2, P3, P4, etc. towards the N-terminus. The complementary sites of the substrate-binding region of the enzyme are correspondingly numbered S2', S1', S1, S2, S3, S4, etc.
The authors thank Sara Barton for careful reading and helpful comments, three anonymous reviewers for their valuable comments, and High Performance Computer Center of Yunnan University for computational support. This work was funded by the National Natural Science Foundation of China (approved number 30630003), and partially funded by grants from the Yunnan University (2004Q013B), Yunnan Province (2006C008M, 2006C0071M), open fund from the Laboratory for Conservation and Utilization of Bio-resources, Yunnan University, and Innovation Group Project from Yunnan University.
- Blow DM: Structure and mechanism of chymotrypsin. Acc Chem Res 1976, 9: 145–152. 10.1021/ar50100a004View ArticleGoogle Scholar
- Tsukuda H, Blow D: Structure of alpha-chymotrypsin refined at 1.68 Å resolution. J Mol Biol 1985, 184: 703–711. 10.1016/0022-2836(85)90314-6View ArticleGoogle Scholar
- Siezen RJ, Leunissen JAM: Subtilases: the superfamily of subtilisin-like serine proteases. Protein Sci 1997, 6: 501–523.PubMed CentralView ArticlePubMedGoogle Scholar
- Siezen RJ, de Vos WM, Leunissen JA, Dijkstra BW: Homology modelling and protein engineering strategy of subtilases, the family of subtilisin-like serine proteinases. Protein Eng 1991, 4: 719–737. 10.1093/protein/4.7.719View ArticlePubMedGoogle Scholar
- Russell AJ, Fersht AR: Rational modification of enzyme catalysis by engineering surface charge. Nature 1987, 328: 496–500. 10.1038/328496a0View ArticlePubMedGoogle Scholar
- Dodson G, Wlodawer A: Catalytic triads and their relatives. Trends Biochem Sci 1998, 23: 347–352. 10.1016/S0968-0004(98)01254-7View ArticlePubMedGoogle Scholar
- Birktoft JJ, Blow DM: Structure of alpha-chymotrypsin. J Mol Biol 1972, l68: 187–240. 10.1016/0022-2836(72)90210-0View ArticleGoogle Scholar
- Clarkson JM, Charnley AK: New insights into the mechanisms of fungal pathogenesis in insects. Trends Microbiol 1996, 4: 197–203. 10.1016/0966-842X(96)10022-6View ArticlePubMedGoogle Scholar
- Yang JK, Huang XW, Tian BY, Wang M, Niu QH, Zhang KQ: Isolation and characterization of a serine protease from the nematophagous fungus, Lecanicillium psalliotae, displaying nematicidal activity. Biotechnol Lett 2005, 27: 1123–1128. 10.1007/s10529-005-8461-0View ArticlePubMedGoogle Scholar
- Andersen SO: Biochemistry of insect cuticle. Annu Rev Entomol 1974, 24: 29–61. 10.1146/annurev.en.24.010179.000333View ArticleGoogle Scholar
- Jansson HB, Nordbring-Hertz B: Infection mechanism in the fungus-nematode system. In Diseases of Nematodes. Volume 2. Edited by: Poinar GO Jr, Jansson HB. Boca Raton Fla: CRC Press; 1988:59–72.Google Scholar
- Pinto FGS, Fungaro MHP, Ferreira JM, Valadares-Inglis MC, Furlaneto MC: Genetic variation in the cuticle-degrading protease activity of the entomopathogen Metarhizium flavoviride. Genet Mol Biol 2002, 25: 231–234. 10.1590/S1415-47572002000200018View ArticleGoogle Scholar
- Khan A, Williams KL, Nevalainen HKM: Effects of paecilomyces lilacinus protease and chitinase on the eggshell structures and hatching of Meloidogyne javanica juveniles. Biol Control 2004, 31: 346–352. 10.1016/j.biocontrol.2004.07.011View ArticleGoogle Scholar
- St Leger RJ, Charnley AK, Cooper RM: Characterization of cuticle-degrading proteases produced by the entomopathogen Metarhizium anisopliae . Arch Biochem Biophys 1987, 253: 221–232. 10.1016/0003-9861(87)90655-2View ArticlePubMedGoogle Scholar
- St Leger RJ, Frank DC, Roberts DW, Staples RC: Molecular cloning and regulatory analysis of the cuticle degrading protease structural gene from the entomopathogenic fungus Metarhizium anisopliae . Eur J Biochem 1993, 204: 991–1001. 10.1111/j.1432-1033.1992.tb16721.xView ArticleGoogle Scholar
- St Leger RJ, Cooper RM, Charnley AK: Cuticle-degrading enzymes of entomopathogenic fungi: degradation in vitro by enzymes from entomopathogens. J Invertebr Pathol 1986, 48: 167–177. 10.1016/0022-2011(86)90043-1View ArticleGoogle Scholar
- Goettel MS, St Leger RJ, Rizzo NW, Staples RC, Roberts DW: Ultrastructural localization of a cuticle-degrading protease produced by the entomopathogenic fungus Metarhizium anisopliae during penetration of host (Manduca sexta) cuticle. J Gen Microbiol 1989, 135: 2233–2239.Google Scholar
- Segers R, Butt TM, Kerry BR, Beckett A, Peberdy JF: The role of the proteinase VCP1 produced by the nematophagous Verticillium chlamydosporium in the infection process of nematode eggs. Mycol Res 1996, 100: 421–428.View ArticleGoogle Scholar
- Morton CO, Hirsch PR, Peberdy JP, Kerry BR: Cloning of and genetic variation in protease VCP1 from the nematophagous fungus Pochonia chlamydosporia. Mycol Res 2003, 107: 38–46. 10.1017/S0953756202007050View ArticlePubMedGoogle Scholar
- Yang JK, Huang XW, Tian BY, Sun H, Duan JX, Wu WP, Zhang KQ: Characterization of an extracellular serine protease gene from the nematophagous fungus Lecanicillium psalliotae . Biotechnol Lett 2005, 27: 1329–1334. 10.1007/s10529-005-0482-1View ArticlePubMedGoogle Scholar
- Betzel Ch, Pal GP, Saenger W: Three-dimensional structure of proteinase K at 0.15 nm resolution. Eur J Biochem 1988, 178: 155–171. 10.1111/j.1432-1033.1988.tb14440.xView ArticlePubMedGoogle Scholar
- Betzel Ch, Teplyakov AV, Harutyunyan EH, Saenger W, Wilson KS: Thermitase and proteinase K: a comparison of the refined three-dimensional structures of the native enzymes. Protein Eng 1990, 3: 161–172. 10.1093/protein/3.3.161View ArticlePubMedGoogle Scholar
- Betzel C, Gourinath S, Kumar P, Kaur P, Perbandt M, Eschenburg S, Singh TP: Structure of a serine protease proteinase K from Tritirachium album limber at 0.98 Å resolution. Biochemistry 2001, 40: 3080–3088. 10.1021/bi002538nView ArticlePubMedGoogle Scholar
- Arnorsdottir J, Kristjansson MM, Ficner R: Crystal structure of a subtilisin-like serine proteinase from a psychrotrophic Vibrio species reveals structural aspects of cold adaptation. FEBS J 2005, 272: 832–845. 10.1111/j.1742-4658.2005.04523.xView ArticlePubMedGoogle Scholar
- Teplyakov AV, Kuranova IP, Harutyunyan EH, Vainstein BK, Frommel C, Hohne WE, Wilson KS: Crystal structure of thermitase at 1.4 Å. J Mol Biol 1990, 214: 261–279. 10.1016/0022-2836(90)90160-NView ArticlePubMedGoogle Scholar
- Wolf WM, Bajorath J, Muller A, Raghunathan S, Singh TP, Hinrichs W, Saenger W: Inhibition of proteinase K by methoxysuccinyl-Ala-Ala-Pro-Ala-chloromethyl ketone. An X-ray study at 2.2-Å resolution. J Biol Chem 1991, 266: 17695–17699.PubMedGoogle Scholar
- Pasternak A, Ringe D, Hedstrom L: Comparison of anionic and cationic trypsinogens: the anionic activation domain is more flexible in solution and differs in its mode of BPTI binding in the crystal structure. Protein Sci 1999, 8: 253–258.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar S, Nussinov R: Different roles of electrostatics in heat and in cold: adaptation by citrate synthase. Chembiochem 2004, 5: 280–290. 10.1002/cbic.200300627View ArticlePubMedGoogle Scholar
- Wells JA, Powers DB, Bott RR, Graycar TP, Estell DA: Designing substrate specificity by protein engineering of electrostatic interactions. Proc Natl Acad Sci USA 1987, 84: 1219–1223. 10.1073/pnas.84.5.1219PubMed CentralView ArticlePubMedGoogle Scholar
- Estell DA, Graycar TP, Miller JV, Powers DB, Burnier JP, Ng PG, Wells JA: Probing steric and hydrophobic effects on enzyme interactions by protein engineering. Science 1986, 233: 659–663. 10.1126/science.233.4764.659View ArticlePubMedGoogle Scholar
- Baker D, Sali A: Protein structure prediction and structural genomics. Science 2001, 294: 93–96. 10.1126/science.1065659View ArticlePubMedGoogle Scholar
- Laskowski RA, MacArthur M, Moss DS, Thornton JM: PROCHECK – a program to check the stereochemical quality of protein structures. J Appl Crystallogr 1993, 26: 283–291. 10.1107/S0021889892009944View ArticleGoogle Scholar
- Müller A, Hinrichs W, Wolf WM, Saenger W: Crystal structure of calcium-free proteinase K at 1.5-Å resolution. J Biol Chem 1994, 269: 23108–23111.PubMedGoogle Scholar
- Helland R, Larsen AN, Smalas AO, Willassen NP: The 1.8 Å crystal structure of a proteinase K-like enzyme from a psychrotroph Serratia species. FEBS J 2006, 273: 61–71. 10.1111/j.1742-4658.2005.05040.xView ArticlePubMedGoogle Scholar
- Bidochka MJ, Khachatourians GG: Basic proteases of entomopathogenic fungi differ in their adsorption properties to insect cuticle. J Invertebr Pathol 1994, 64: 26–32. 10.1006/jipa.1994.1064View ArticleGoogle Scholar
- Leal SCM, Bertioli DJ, Butt TM, Carder JH, Burrows PR, Peberdy JF: Amplification and restriction endonuclease digestion of the Pr1 genes for the detection and characterization of Metarhizium strains. Mycol Res 1997, 101: 257–265. 10.1017/S0953756296002560View ArticleGoogle Scholar
- Segers R, Butt TM, Keen JN, Kerry BR, Peberdy JF: The subtilisins of the invertebrate mycopathogens Verticillium chlamydosporium and Metarhizium anisopliae are serologically and functionally related. FEMS Microbiol Lett 1995, 126: 227–231. 10.1111/j.1574-6968.1995.tb07423.xView ArticlePubMedGoogle Scholar
- Bidochka MJ, Melzer MJ: Genetic polymorphisms in three subtilisin-like protease isoforms (Pr1A, Pr1B, and Pr1C) from Metarhizium strains. Can J Microbiol 2000, 46: 1138–1144. 10.1139/cjm-46-12-1138View ArticlePubMedGoogle Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234: 779–815. 10.1006/jmbi.1993.1626View ArticlePubMedGoogle Scholar
- Lindahl E, Hess B, van der Spoel D: GROMACS 3.0: A package for molecular simulation and trajectory analysis. J Mol Mod 2001, 7: 306–317.Google Scholar
- Brendel V, Bucher P, Nourbakhsh I, Blaisdell BE, Karlin S: Methods and algorithms for statistical analysis of protein sequences. Proc Natl Acad Sci USA 1992, 89: 2002–2006. 10.1073/pnas.89.6.2002PubMed CentralView ArticlePubMedGoogle Scholar
- The Isoelectric point service[http://www.embl-heidelberg.de/cgi/pi-wrapper.pl]
- Gill SC, von Hipple P: Calculation of protein extinction coefficients from amino acid sequence data. Anal Biochem 1989, 182: 319–326. 10.1016/0003-2697(89)90602-7View ArticlePubMedGoogle Scholar
- Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22: 2577–2637. 10.1002/bip.360221211View ArticlePubMedGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis 1997, 18: 2714–2723. 10.1002/elps.1150181505View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22: 4673–4680. 10.1093/nar/22.22.4673PubMed CentralView ArticlePubMedGoogle Scholar
- Schechter I, Berger A: On the size of the active site in proteases. I. Papain. Biochem Biophys Res Commun 1967, 27: 157–162. 10.1016/S0006-291X(67)80055-XView ArticlePubMedGoogle Scholar