Search tips
Search criteria 


Logo of plospathPLoS PathogensSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
PLoS Pathog. 2017 February; 13(2): e1006169.
Published online 2017 February 13. doi:  10.1371/journal.ppat.1006169
PMCID: PMC5325619

Functional diversity of secreted cestode Kunitz proteins: Inhibition of serine peptidases and blockade of cation channels

Paul J. Brindley, Editor


We previously reported a multigene family of monodomain Kunitz proteins from Echinococcus granulosus (EgKU-1-EgKU-8), and provided evidence that some EgKUs are secreted by larval worms to the host interface. In addition, functional studies and homology modeling suggested that, similar to monodomain Kunitz families present in animal venoms, the E. granulosus family could include peptidase inhibitors as well as channel blockers. Using enzyme kinetics and whole-cell patch-clamp, we now demonstrate that the EgKUs are indeed functionally diverse. In fact, most of them behaved as high affinity inhibitors of either chymotrypsin (EgKU-2-EgKU-3) or trypsin (EgKU-5-EgKU-8). In contrast, the close paralogs EgKU-1 and EgKU-4 blocked voltage-dependent potassium channels (Kv); and also pH-dependent sodium channels (ASICs), while showing null (EgKU-1) or marginal (EgKU-4) peptidase inhibitory activity. We also confirmed the presence of EgKUs in secretions from other parasite stages, notably from adult worms and metacestodes. Interestingly, data from genome projects reveal that at least eight additional monodomain Kunitz proteins are encoded in the genome; that particular EgKUs are up-regulated in various stages; and that analogous Kunitz families exist in other medically important cestodes, but not in trematodes. Members of this expanded family of secreted cestode proteins thus have the potential to block, through high affinity interactions, the function of host counterparts (either peptidases or cation channels) and contribute to the establishment and persistence of infection. From a more general perspective, our results confirm that multigene families of Kunitz inhibitors from parasite secretions and animal venoms display a similar functional diversity and thus, that host-parasite co-evolution may also drive the emergence of a new function associated with the Kunitz scaffold.

Author summary

Parasite secretions are key players at host-parasite interfaces: parasite establishment and persistence rely, to a great extent, on interactions between these molecules and their host counterparts. We present the functional characterization of a multigene family of secreted Kunitz proteins from the cestode Echinococcus granulosus. Kunitz proteins are a class of metazoan high affinity serine peptidase inhibitors. In addition, families of Kunitz proteins are frequent components of animal venoms; besides inhibiting peptidases, some of these “Kunitz toxins” block cation channels and thus, provide a remarkable example of protein evolution where natural selection has driven the emergence of a new function associated with the same molecular scaffold. Using enzyme kinetics and electrophysiological assays, we demonstrated that the E. granulosus Kunitz family includes peptidase inhibitors as well as channel blockers. This diversity highlights an interesting similarity between animal toxins and parasite secretions that had not been previously described. Furthermore, the presence of analogous families of Kunitz proteins appears to be a distinctive trait of cestode genomes. We thus propose that these molecules contribute to a successful infection acting at the parasite-host interface. In addition, because they bear a strong specificity towards their targets, they are uniquely suited for the development of pharmaceuticals.


Cestodes are a neglected group of platyhelminth parasites, despite causing chronic infections to humans and domestic animals worldwide [1]. Together with other researchers around the world [2], we have been using Echinococcus granulosus as a model to study the molecular basis of the host-parasite cross-talk during cestode infections [3,4,5]. E. granulosus is the agent of cystic echinococcosis, a medically and economically important worldwide zoonosis, with endemic foci in Central Asia, China, South America and Africa [6]. Like all taeniid cestodes, it has a life cycle involving two mammals: a non carnivore intermediate host (harboring the larva) and a carnivore definitive host (harboring the hermaphroditic adult). Intermediate hosts (ungulates such as sheep, cattle and pigs; and, accidentally, also humans) become infected by ingestion of eggs containing oncospheres that develop at visceral sites into bladder-like metacestodes (hydatid cysts). These latter are bounded by a wall whose inner germinal layer gives rise to larval worms (protoscoleces) by asexual budding; protoscoleces are bathed in hydatid fluid that includes host plasmatic proteins and parasite secretions. Infection in the definitive host (always a canid, most often dogs) arises from ingestion of protoscoleces that, upon activation by contact with stomach acid, enzymes and bile acids, evaginate and attach to the mucosa of the duodenum, where they develop into adult tapeworms that can reside in the gut for long periods without causing any apparent damage [7]. Specific anatomical structures allow such a close contact at the canid-worm interface that E. granulosus has been regarded as both a tissue and a luminal parasite [8]. The molecular mechanisms underlying its successful establishment and persistence in the hostile environment of the dog duodenum are unknown.

With the aim of identifying molecules participating in the E. granulosus–dog cross-talk, we surveyed the genes expressed by protoscoleces and pepsin/H+-treated protoscoleces. We thus identified a multigene family of Kunitz-type inhibitors (EgKUs). These molecules were associated mostly with treated protoscoleces, suggesting that they play roles at the initial phases of infection [3]. Kunitz inhibitors are a class of metazoan serine peptidase inhibitors, whose prototype is the bovine pancreatic inhibitor of trypsin (BPTI; family I2 of the MEROPS database; [9]. They are competitive inhibitors acting in a substrate-like manner, that form very stable complexes of 1:1 stoichiometry with their target enzymes, devoid of activity. The interaction between the enzyme and the inhibitor is highly dependent on the residue located at the position P1 of the antipeptidase loop (position 15 of mature BPTI) [10]. In addition, families of Kunitz inhibitors are frequent components of the saliva and secretions from hematophagous animals and also of animal venoms. These “Kunitz-type toxins” have been described in the venoms from snakes [11], sea anemones [12,13], cone snails [14], spiders [15], scorpions [16,17] as well as in the saliva of blood-sucking arthropods [18,19] and in the secretions of hematophagous nematodes [20]. Interestingly, besides inhibiting peptidases, some Kunitz toxins, until now described only in animal venoms, block various types of cation channels. Furthermore, some act solely as channel blockers. A set of neurotoxins present in the venoms of mamba snakes (‘‘dendrotoxins”), whose function is to paralyze the prey, is the best known example [21].

We previously reported the molecular features of eight EgKUs (that we named EgKU-1-EgKU-8) and provided evidence that some of them (notably, EgKU-3 and EgKU-8) are secreted by protoscoleces. Although diverse, these EgKUs were found to group into three pairs of close paralogs (EgKU-1/EgKU-4; EgKU-3/EgKU-8; EgKU-6/EgKU-7), which would be the products of recent gene duplications. In addition, we carried out detailed kinetic studies with native EgKU-1 and EgKU-8 purified from protoscoleces that revealed their possible functionalities. EgKU-8 behaved as a slow, tight-binding inhibitor of trypsins, with global inhibition constants (KI*) in the 10−11 M range, and interacted with enzymes through a mechanism involving two reversible steps; an initial relatively fast formation of an enzyme-inhibitor complex followed by a slow transition to a tight complex. In sharp contrast, EgKU-1 did not inhibit any of the assayed peptidases. Interestingly, molecular modeling revealed that structural elements associated with activity in Kunitz cation-channel blockers are also present in EgKU-1. Indeed, α-dendrotoxin (α-DTX), a well characterized blocker of specific voltage-activated K+-channels (Kv) [21], was—at the time—the best overall template of EgKU-1; and several amino acids important for toxin activity were found to be conserved in the consensus model of the parasite molecule, supporting the notion that it is a putative cation channel blocker. Presumed orthologs of the EgKUs (peptidase inhibitors as well as channel blockers) were also found to be present in the transcriptomes from the other medically important cestodes (notably, E. multilocularis and Taenia solium, the agents of alveolar echinococcosis and cysticercosis, respectively), indicating that families of monodomain Kunitz inhibitors are also present in closely related organisms [3].

In this article, we characterize the activity of EgKU-1–EgKU-8 using enzyme kinetics and whole-cell patch clamp assays. We thus demonstrate that the E. granulosus Kunitz family is indeed functionally diverse. On the one hand, we show that all but EgKU-1 and EgKU-4 behave as high affinity inhibitors of either chymotrypsin or trypsin. On the other hand, patch-clamp assays on rat dorsal root ganglion (DRG) neurons confirmed that EgKU-1, and also its close paralog EgKU-4, block Kv. Furthermore, the two proteins also block pH-dependent sodium channels (acid sensing ion channels, ASICs), a previously unreported activity for Kunitz inhibitors, that we recently described for α-DTX [22]. In addition, we provide further evidence of the presence of EgKUs in parasite secretions. We discuss the significance of these results taking into account available genomic and transcriptomic data from E. granulosus and related cestodes.


Except for EgKU-1 and EgKU-4, the EgKUs are high affinity inhibitors of serine peptidases

In our previous study, working with native EgKU-1 and EgKU-8, we demonstrated that EgKU-8 is a high affinity inhibitor of trypsins, whereas EgKU-1 did not inhibit any of the assayed peptidases [3]. To further advance in the functional characterization of the family, we prepared recombinants of the eight EgKUs. We carried out a preliminary screening of the serine peptidase inhibition activity of recombinant EgKU-2–EgKU-7. Of note, in the case of EgKU-8 both the native inhibitor and the recombinant protein behaved similarly (KI* 60 ± 13 versus 50 ± 10 pM, for native and recombinant EgKU-8, respectively). Using pancreatic enzymes, we analyzed whether the EgKUs showed the inhibition profiles that may be predicted from the respective amino acid at position P1: EgKU-2 (Trp in P1) and EgKU-3 (Leu in P1) are candidate chymotrypsin inhibitors, whereas EgKU-4-EgKU-7 (Arg in P1) are predicted to inhibit trypsin. All six EgKUs showed the expected activities. We subsequently performed titration assays to analyze whether they behaved as high affinity inhibitors. These studies indicated that, except for EgKU-4, the close paralog of EgKU-1, the EgKUs are high affinity inhibitors of bovine chymotrypsin (EgKU-2 and EgKU-3) or trypsin (EgKU-5-EgKU-8) (Table 1 and Fig 1). In view of these results, we further characterized the inhibitory activity of EgKU-3, the closest paralog of EgKU-8.

Fig 1
Titration assays of recombinant EgKUs: results for EgKU-3 and EgKU-4.
Table 1
Screening of serine peptidase inhibitory activity of EgKUs.

EgKU-3 inhibits chymotrypsin with high affinity through a two-step mechanism

We carried out kinetic studies using bovine chymotrypsin A and also chymotrypsin purified from dog pancreas, i. e. chymotrypsin B (chymotrypsin A is absent from dogs, see S01.001 at MEROPS— EgKU-3 inhibited with high affinity both peptidases; Fig 2A shows a representative experiment with the bovine enzyme and Table 2 the global inhibition constants calculated for the two chymotrypsins. The values of KI* were of the same order (53 ± 19 and 84 ± 49 pM for the bovine and canine enzymes, respectively) indicating no bias in specificity towards any of them. EgKU-3 also inhibited elastase, although with substantially lower affinity than chymotrypsins (KI* of 5 ± 2 nM).

Fig 2
Inhibition studies with EgKU-3: results for bovine chymotrypsin A.
Table 2
Global inhibition constants (KI*) of EgKU-3, EgKU-4 and EgKU-8 acting on pancreatic serine peptidases.

In order to study the inhibition mechanism of EgKU-3 towards chymotrypsins, we carried out time course experiments with chymotrypsin A. The progress curves for the inhibition (Fig 2B) indicated that the enzyme-inhibitor complex reaches equilibrium in a time scale of minutes and that EgKU-3 is a slow-binding inhibitor as defined by Morrison [23]. The interaction of EgKU-3 with chymotrypsin was reversible, since progress curves reached appreciable slopes even at higher than stoichiometric inhibitor concentrations. This is the expected behavior for Kunitz-type inhibitors [10] and the one observed for EgKU-8 [3]. Similarly, the plot of the apparent rate constant (kobs) versus EgKU-3 concentration was hyperbolical (Fig 2C), in accordance with a mechanism involving two steps, a fast initial binding of the inhibitor to the target enzyme followed by a slow transition [24]. The kinetic constants of EgKU-3 binding to chymotrypsin obtained from analyses of the progress curves are shown in Table 3. Note that the value of KI* calculated from the kinetic constants compared very well with the value obtained through the fit of steady-state rate versus inhibitor concentration data to the Morrison equation (Table 2).

Table 3
Inhibitory kinetics of EgKU-3 on bovine chymotrypsin A.

EgKU-1 and EgKU-4 block voltage-activated potassium channels (Kv)

As already mentioned, our results indicate that the paralogs EgKU-1/EgKU-4 do not show the typical serine peptidase inhibitory activity of Kunitz-type inhibitors. In fact, EgKU-1 did not inhibit any assayed peptidase [3]; whereas EgKU-4 inhibited trypsin albeit with low affinity, with a KI* of 47 ± 2 nM, i. e. 1000-fold higher than the KI* of EgKU-3 and EgKU-8 versus their target enzymes (Table 2). In view of these results and taking into account the structural similarity between EgKU-1 and α-DTX [3], we analyzed whether EgKU-1 and EgKU-4 acted on Kv using whole-cell patch-clamp assays on neurons isolated from DRG. Both EgKUs inhibited Kv; Fig 3A illustrates the effect of recombinant EgKU-1. The blockade was more pronounced over the steady-state component of the current than over the peak current: 25 ± 11% versus 20 ± 14% and 27 ± 12% versus 23 ± 7% for 200 nM of recombinant EgKU-1 and EgKU-4, respectively (n = 7). We also tested the activity of native EgKU-1 and verified that the recombinant inhibitor reproduced reasonably well the behavior of the native inhibitor (100 nM of native EgKU-1 blocked the steady-state current by 19% and the peak current by 10%; n = 4). Thus, the effect appeared to be stronger on the currents at the end of the pulse (accounting for non inactivating -delayed-rectifier- K+ currents, IKDR) than on those at the beginning of the pulse (corresponding to fast -transient A-type- K+ currents, IKA). The effect was only partially reversible because about 60% persisted 3 min after washing (Fig 3B). In addition, it was clearly observed on the currents elicited over -40 mV, as highlighted by the activity profile of the EgKUs over the K+ currents activated by different voltages (Fig 3D–3F). In contrast, the perfusion of 1 μM EgKU-8 (n = 7) produced no significant changes in the peak amplitude (2.6 ± 3.2%, P = 0.20) or the steady-state current (1.1 ± 3.5%, P = 0.38), whereas 1 μM of EgKU-3 (n = 10) produced a slight non-significant reduction of the peak current (6.3 ± 3.8%, P = 0.13) and had no effect on the steady-state current (4.1 ± 4.5%, P = 0.13) (Fig 4). We also analyzed the effect of EgKU-1 and EgKU-4 on voltage-activated sodium channels (Nav) and observed no effect (S1 Fig).

Fig 3
Inhibition studies with EgKU-1 and EgKU-4: results for Kv from DRG neurons.
Fig 4
Studies with EgKU-3 and EgKU-8 on total K+ currents from DRG neurons.

Subsequently, we further characterized the effect of EgKU-1 and EgKU-4 on Kv. For this purpose, we recorded the currents after a pre-pulse of -120 mV to activate all voltage-dependent K+ currents, transient IKA as well as slow-inactivating IKDR; and also those remaining when the pre-pulse was of -45 mV, voltage at which IKA is inactivated. Thus, the second recording corresponded to IKDR; whereas IKA could be deduced by subtracting IKDR from the first recording. A representative experiment with recombinant EgKU-1 is shown in Fig 5.This setup allowed us to analyze the effect of the inhibitors on total K+ currents (Fig 5A and 5C) as well as on both types of isolated K+ currents, IKDR (Fig 5B and 5D) and IKA (Fig 5E). As anticipated by the previous experiment, EgKU-1 principally affected IKDR, with virtually no effect on IKA (Fig 5D versus 5E).5E). Finally, we studied the concentration-response relationship for native EgKU-1 and estimated an IC50 of about 200 nM when acting on all K+ currents activated by a pulse of -100 to 0 mV (Fig 6). Although we did not determine the IC50 for EgKU-4, its behavior was similar to the one of EgKU-1.

Fig 5
Inhibition studies with EgKU-1 on isolated K+ currents from DRG neurons.
Fig 6
Concentration-response analysis of native EgKU-1 on total K+ currents from DRG neurons.

EgKU-1 and EgKU-4 block Acid Sensing Ion Channels (ASICs)

Taking into account the recently described activity of α-DTX on ASIC currents in DRG neurons [22], we also analyzed the effect of EgKU-1 and EgKU-4 on pH-dependent Na+ currents. The sustained application of both EgKUs blocked the ASIC currents elicited by a pH change from 7.4 to 6.1. The blocking effect was on the peak amplitude (Ipeak) and no significant effect was observed on the desensitization time course (τdes). The effect on the current amplitude was fully reversible after 1 min washing (Fig 7A–7C). We similarly analyzed the effect of EgKU-3 and EgKU-8: EgKU-3 (n = 8) produced a slight but significant decrement of the Ipeak (6.4 ± 3.0%; P = 0.02), whereas EgKU-8 (n = 7) had no effect (6.2 ± 5.1%; P = 0.11) (Fig 7D–7E). Finally, we studied the concentration-response relationship for native EgKU-1 (Fig 8A); the estimated IC50 value (about 8 nM) was 25-fold lower than the one determined for Kv (about 200 nM), suggesting higher selectivity of EgKU-1 for ASICs than Kv.

Fig 7
Inhibition studies with EgKU-1 and EgKU-4: results for ASIC currents from DRG neurons.
Fig 8
Concentration-response analysis of native EgKU-1 on ASIC currents from DRG neurons.

The EgKUs may be detected in metacestode and adult worm secretions

In our previous study, using mass spectrometry analysis, we showed that members of the Kunitz family, notably EgKU-3 and EgKU-8, would be present in protoscolex secretions from untreated and pepsin/H+-treated larval worms [3]. To further approach the question of whether Kunitz inhibitors are secreted to the parasite-host interface, we similarly analyzed hydatid fluid and adult worm secretions. Fig 9A shows a representative MALDI-TOF MS profile (5000–10000 Da) of hydatid fluid from bovine cysts. Peaks of m/z 6407.5 and 6519.6, matching the predicted MH+ value for EgKU-3 (6406.4 Da) and EgKU-8 (6520.4 Da), respectively, were observed. Furthermore, the intensity of the signals putatively corresponding to both EgKUs was substantially increased in the chymotrypsin A-affinity purified fraction from the same sample (Fig 9B). In addition, the MS profile of an analogous fraction from the supernatant of in vitro cultured immature adults (Fig 9C) also showed peaks matching the predicted MH+ value for EgKU-3 and EgKU-8 (m/z of 6409.8 and 6521.3, respectively). Peptide mass fingerprinting of the components purified from cyst fluid allowed the detection of signals that could be assigned to tryptic peptides of these EgKUs. In particular, we detected signals with m/z values 1074.54 and 1491.63 that corresponded to EgKU-8 sequences 7LPLDPGFCR15 and 21WGFHQESGECVR32 with S-carboxymethylated cysteines (theoretical m/z values 1074.53 and 1491.64, respectively; see Fig 10C). In addition, a signal corresponding to EgKU-3 sequence 49–57 was detected in the same spectrum. MS/MS analysis of the corresponding ions further corroborated the amino acidic sequences (S1 Dataset). Note that, although lower than towards trypsins, the affinity of EgKU-8 towards chymotrypsin A (KI* 10−9 M; Table 2) was high enough to allow its purification from the secretions. In contrast, this approach is not suitable to purify EgKU-1 and EgKU-4 that have no affinity towards chymotrypsin; thus, even if they had been present in the original sample, we would not have detected them.

Fig 9
Detection of EgKU-3 and EgKU-8 in parasite secretions.
Fig 10
Structural analyses of EgKU-1/EgKU-4 and EgKU-3/EgKU-8.


In the present study, we described the functional characterization of eight members of the E. granulosus family of secreted Kunitz inhibitors (EgKUs). Using recombinant forms of EgKU-1-EgKU-8 and native EgKU-1, we demonstrated that six EgKUs behave as high affinity inhibitors of either chymotrypsin (EgKU-2 and EgKU-3) or trypsin (EgKU-5-EgKU-8), whereas the close paralogs EgKU-1/EgKU-4 act as cation channel blockers (of Kv as well as ASICs), while showing either null (EgKU-1, our previous study [3]) or marginal (EgKU-4) serine peptidase inhibition activity. This degree of functional diversity, commonly observed in animal venoms, had not been previously described for Kunitz inhibitors present in parasite secretions.

Regarding serine peptidase inhibition, detailed kinetic studies showed that the interaction of EgKU-3 with chymotrypsins mimics the one of the close paralog EgKU-8 with trypsins: it is slow, of very high affinity and involves two steps. EgKU-3 strongly inhibited isoforms A and B of chymotrypsin with KI* in the 10−11 M range. Notably, according to MEROPS, dogs have two chymotrypsins B (with > 95% identity, encoded by CTRB1 and CTRB2 genes) and lack chymotrypsin A. The values of KI* are among the smaller registered for chymotrypsin inhibitors. In fact, only two peptides have been reported to have similar affinity, both towards bovine chymotrypsin A (see S01.001 in MEROPS) [25,26]. No high affinity inhibitors of chymotrypsin B have been described so far, most likely because very few studies have been carried out with this isoform (see S01.152 in MEROPS). Therefore, EgKU-3 appears as an interesting titration reagent for chymotrypsins A and B, especially for the latter because adequate titration reagents are currently unavailable. The stability of the EgKU-3-chymotrypsin A complex is similar to that of EgKU-8 complexes with trypsins, with k2/KI, the apparent second order rate constant for complex formation (kon), in the 106 M-1 s-1 range; and the dissociation rate constant (k-2) in the 10−4 s-1 range (Table 3 and our previous study, [3]). These values of k2/KI are in good agreement with reports for other members of the family, including BPTI with bovine trypsin [27], whereas those of k-2 are several orders faster than the one reported for BPTI (10−8 s-1 with bovine trypsin [27]).

The activity of EgKU-3 as a strong tight-binding inhibitor of chymotrypsins and a less potent inhibitor of pancreatic elastase, as well as its lack of activity towards trypsin (Table 2) are consistent with the presence of a Leu in P1. In turn, similar to BPTI (Lys in P1), EgKU-8 (Arg in P1) strongly inhibits trypsins, less potently chymotrypsin A and does not inhibit elastase (Table 2). Notably, neither EgKU-8 nor BPTI inhibit chymotrypsin B. The antipeptidase loops of EgKU-3 and EgKU-8 differ by 50% (the corresponding mature polypeptides differ by only 37%); differences are mainly on the P side of the loop (residues 10 to 15), and involve some non-conservative substitutions (notably Asp10 instead of Lys10 in P6; see Fig 10C). Because the loop contributes to the inhibition specificity primarily determined by the P1 site, some of these residues are likely involved in the interaction with chymotrypsin B. In any case, consistent with the fact that isoforms A and B have similar affinities towards substrates with Phe as P1 residue [28], the values of KM for the substrate we used were of the same order.

Regarding cation channel inhibition, patch-clamp studies carried out on rat DRG neurons showed that EgKU-1 and EgKU-4 block voltage-activated potassium currents. The effect was voltage-dependent and, as described for dendrotoxins, it was not totally reversible (Fig 3; [29]). The detailed characterization of EgKU-1 activity on isolated K+ currents indicated that it preferentially blocks IKDR, as compared to IKA. This behavior differs from the one of α–DTX and resembles the one of δ–DTX [29]. In any case, the IC50 determined for EgKU-1 was two orders of magnitude higher than those of dendrotoxins assayed on DRG neurons (10−7 versus 10−9 M [29]). This could be due, at least in part, to the fact that, although EgKU-1 shares with dendrotoxins several residues that would participate in channel interaction (notably Leu7, Lys26 and Lys27), it lacks the Lys5 that has been described as a primary determinant of activity ([21]; this is also the case for EgKU-4 that shares Leu7 and Lys27 with α–DTX and EgKU-1; see Fig 10C). Because we studied the effect of the EgKUs on total Kv currents of primary cultures, we cannot comment on their activity over specific Kv. Nevertheless, our results indicate that EgKU-1 could be more active over some Kv than others. Indeed, the dose-response curve does not start from zero (Fig 6), as if a specific Kv was highly blocked at low concentrations of the inhibitor.

EgKU-1 and EgKU-4 also showed a potent dose-dependent blocking effect on the ASIC currents in DRG neurons, which was totally reversible after one minute washing (Fig 7). These neurons express at least two subpopulations of transient ASIC currents as judged by their inactivation constants [30]. One of them derives from channels of ASIC1a, ASIC1b and ASIC3 subunits; the other from channels of ASIC2a subunits (reviewed by [31]), which are the least expressed in DRG [32]. Although our experimental setup does not allow us to conclude which channels are sensitive to the EgKUs, the dispersion of the values in the dose-response study with native EgKU-1 (Fig 8A) points to some variability of the blocking effect among different cells, suggesting that the effect could be stronger for some channel(s). EgKU-1 could thus mimic the performance of other peptide blockers of ASICs, such as APETx2 (reviewed by [31]).

We recently reported that α-DTX, the well-known blocker of voltage activated K+ channels, also inhibits ASIC currents in rat DRG, although with significantly less potency than Kv (IC50 ~ 10−7 M [22] versus 10−9 M [29], respectively; see also Figs Figs44 and and8B).8B). This result indicates that the Kunitz domain is yet another structural scaffold for ASIC-blocking polypeptides. Interestingly, an exposed basic-aromatic cluster identified in structurally different ASIC blocking peptides [33] was also found to be present in the structure of α-DTX [22]. Notably, this feature is observable towards one side of the model structures of EgKU-1 and EgKU-4 and not in those of EgKU-3/EgKU-8 (Fig 10A). In any case, functionally distinct EgKUs differ mainly in surface charge distribution (Fig 10B). The relatively low selectivity of Kunitz inhibitors towards cation channels contrasts with the high specificity of their interaction with serine peptidases. Not surprisingly, structure-activity analyses focused at identifying the “channel-blocking site” of Kunitz proteins have usually highlighted regions on their surface involved in channel interaction but not a defined structural motif comparable to the antipeptidase loop involved in serine peptidase interaction. Furthermore, key residues for channel blockade are frequently located in the N- or C-terminal extensions of the Kunitz domain [17,34]

The availability of the E. granulosus genome [5,35] has allowed us to identify genes coding for at least eight additional monodomain Kunitz proteins with the same molecular architecture as EgKU-1–EgKU-8, i. e. a signal peptide followed by a single Kunitz domain. Similar to the rest of the family, the newly-identified members are diverse and include several pairs of close paralogs, consistent with an accelerated evolution of the family. Fig 11 shows an unrooted phylogenetic tree of the Kunitz domains from the sixteen EgKUs together with eleven close paralogs from T. solium and five from functionally characterized monodomain Kunitz proteins from Lophotrochozoa, including four from trematodes. A true phylogenetic tree is not intended, as the signal might be blurred by homoplasy. Rather, the tree is aimed to mirror functional groupings of the sequences in an approximate evolutionary context. Not surprisingly, the sequences from T. solium pair with their close E. granulosus paralogs. The groupings roughly correlate with functional features, whereas EgKU-2 (and a putative T. solium ortholog) appears very distant from the rest. The red sub-clade includes several serine peptidase inhibitors: in addition to EgKU-3/EgKU-8 and EgKU-5, EGR_07242 (EgKI-2 in [36]) and the schistosome proteins SjKI-1 [37] and SmKI-1 [38]. EGR_07242 (Arg in P1) was recently found to inhibit trypsin, although with relatively low affinity (KI ~ 10−9 M; [36]), probably due to the lack of Cys14, i. e. the one forming the disulphide bond that stabilizes the antipeptidase loop. SjKI-1 and SmKI-1 (both with Arg in P1) also inhibit trypsin with IC50 in the 10−10 and 10−8 M range, respectively [37,38]. The green sub-clade appears to group a different set of serine peptidase inhibitors (EgKU-6/EgKU-7 and closely related proteins from T. solium). In turn, the blue sub-clade includes the channel blockers EgKU-1/EgKU-4 together with another pair of close E. granulosus paralogs (EgrG001136600/EgrG001137000), and two T. solium proteins (TsM_000410200 and TsM_000513000). Although it is difficult to predict their function without further data, these proteins could also act as channel blockers because, similar to EgKU-1/EgKU-4, they feature the conserved Leu7 and a positively charged β–turn that form the Kv-blocking site of α-DTX and related toxins [34,39,40,41]. As to the other groupings, EgrG_001136500 (Leu in P1) was recently found to be a potent inhibitor of neutrophil elastase (KI ~ 10−11 M) and cathepsin G (KI ~ 10−10 M) and, interestingly, to reduce neutrophil infiltration in a local inflammation model (EgKI-2 in [36]); thus, its close paralog (EgrG_00113800; Arg in P1) could also be a serine peptidase inhibitor. Finally, the sequences from Fasciola hepatica (FhKTM [42] and FhKT1 [43], both with Leu in P1, whose Kunitz domains differ in 3/51 amino acids) define a basal, separate sub-clade, that could also reflect functional diversity: FhKTM was found to be a marginal inhibitor of trypsin with virtually no effect over chymotrypsin [42] but, notably, FhKT1 was recently characterized as an inhibitor of cysteine peptidases, including the major parasite cathepsin L secreted peptidases and related human peptidases [43]).

Fig 11
Expanded view of the E. granulosus Kunitz family.

Another interesting finding of our work refers to the demonstration of the presence of some EgKUs (notably, EgKU-3 and EgKU-8) in cyst fluid (from bovine cysts) and secretions from immature adult worms, which complement our previous results with secretions from protoscoleces and pepsin/H+-treated protoscoleces [3]. EgKU-8 was also detected by proteomic analyses in fertile cyst fluids from ovine and human infections, but not in infertile cysts from infected cattle [44]. Members of the Kunitz family could thus be secreted to the E. granulosus-dog interface not only at the initial stages of infection (as indicated by their presence in larval worm secretions) but also at late stages, and contribute to the establishment and persistence of dog echinococcosis. In turn, their presence in cyst fluid would point to a role at the onset of infection in dogs and/or during the chronic stage of infection in intermediate hosts.

In addition, available RNASeq data [5,35] indicate that members of the family are expressed in all the analyzed stages (immature adult, activated oncosphere, cyst, protoscolex, pepsin/H+-treated protoscolex); interestingly, most EgKUs are expressed in adults. Furthermore, several of them appear up-regulated in this stage, notably, EgKU-3, EgKU-7, EgKU-8, and also EGR_07242 and EgrG_0001137000 [35]. In addition, the orthologs of EgKU-1 and EgKU-4 were found to be highly up-regulated in E. multilocularis gravid adults (as compared to non-gravid worms) [5], whereas EgrG_001136500 and a close paralog (EmuJ_001136900) were among the transcripts with higher expression in oncospheres from E. granulosus [35] and E. multilocularis [45], respectively. Together with our results highlighting the presence of EgKUs in parasite secretions, these data support the concept that Kunitz proteins are involved in parasite interaction with definitive and intermediate hosts, and indicate that specific members of the family would be engaged in particular moments of the life-cycle.

Given the activity profile of EgKU-1-EgKU-8 and the fact that they are mostly expressed and secreted by larval and adult worms [3,35], we can speculate about their potential counterparts in the dog duodenum. The apical end of the scolex contains a gland (the rostellar gland) whose secretion is believed to play a key role in host-parasite cross-talk, due to the very intimate contact of the scolex with the mucosa (reviewed by [2]). Interestingly, seminal studies demonstrated that the rostellar gland secretion is cystine-rich [46]; the gland could thus be the site of synthesis and concentration of the EgKUs. Pancreatic enzymes appear as clear targets of the EgKUs acting as serine peptidase inhibitors: those present in cyst fluid could initially protect the larval worms from digestion; whereas those secreted could protect the scolex (whose glycocalix is thin) after the parasite has attached to the mucosa. The EgKUs could also inhibit other serine peptidases, such as those secreted by immune cells (see for example [47]) and membrane peptidases from epithelial cells (see for example [48,49]); in turn, this effect could prevent the activation of proteinase-activated receptors (PARs; [50]). As to putative targets of EgKU-1/EgKU-4, very little is known about the expression and functional properties of Kv and ASICs in the gut. In any case, both types of channels participate in the physiology of epithelial cells and afferent neurons [51,52,53]. In addition, Kv and ASICs are involved in the activation and maturation of dendritic cells and macrophages. In particular, Kv1.3 and Kv1.5 modulate the Ca2+-dependent functions of these cells and their blockade down-regulates their activation [54]. Regarding ASICs, dendritic cells express ASIC1, ASIC2 and ASIC3, and extracellular acidosis induces currents that are blocked by ASIC inhibitors. In addition, acidosis triggers the activation of dendritic cells and macrophages and ASIC inhibitors block these effects [54,55]. Taking into account that extracellular acidosis is a hallmark of inflammation, the blockade of ASICs may be crucial to weaken the induction of innate immunity and to favor the development of a chronic infection. In this context, it is pertinent to mention that excretory-secretory (E/S) products from E. granulosus adults have recently been found to impair dendritic cell function and induce the development of regulatory T cells [56]. Furthermore, a similar result was observed with the Kunitz protein FhKTM from the trematode F. hepatica, which is known to be present in parasite E/S products [57].As already mentioned, FhKTM showed a marginal serine peptidase inhibitory activity [42] but, notably, the very closely related FhKT1 (> 90% overall identity with FhKTM) was recently found to inhibit cysteine peptidases [43].

In our previous study, we included data from an extensive survey of platyhelminth ESTs available at the time, indicating that the expression of families of monodomain Kunitz proteins would be a distinctive trait of cestodes [3]. Genomic and transcriptomic data currently available (accessible from WormBase ParaSite:; [58]) confirm and extend these initial observations, in particular, for the other medically important cestodes, E. multilocularis and T. solium, and indicate that this family is expanded in cestodes. Putative orthologs of virtually all the EgKUs are present in the E. multilocularis genome, and several are also predicted for T. solium (Fig 11). In addition, recent genome-wide analyses of E/S proteins showed that the Kunitz domain is either the most (in E. multilocularis, 17/673 predicted E/S proteins; [59]) or the third most (in T. solium, 14/838 predicted E/S proteins; [60]) represented domain in the predicted secretome of these parasites, whereas manual inspection of the putative secreted Kunitz proteins indicates that a majority of them contain a single Kunitz domain. In contrast, parasitic trematodes express only a few monodomain Kunitz inhibitors (see for example [61]).

The secretion of monodomain Kunitz proteins thus appears to be a strategy evolved by cestodes to block, through high affinity interactions, the function of host proteins (either serine peptidases or cation channels) and contribute to the establishment and persistence of infection. The putative immunomodulatory role of these molecules opens the way to further studies of their involvement in immunoevasion, acting as single molecules as well as synergistically. From a more general perspective, the data confirm that multigene families of Kunitz inhibitors from parasite secretions and animal venoms display a similar functional diversity. As we had previously mentioned, because the genes coding for parasite secretions and predator toxins arise from an arms race between different organisms, it is interesting to consider that both sets of molecules display analogous evolutionary patterns. Finally, the strong target specificity of some of these molecules makes them uniquely suited as tools for the characterization of biological processes as well as for the development of pharmaceuticals.

Materials and methods

Preparation of EgKUs

Native EgKU-1 and EgKU-8 were purified to homogeneity from a protoscolex lysate by cation exchange followed by reverse-phase chromatography, as previously described [3].

EgKU-1–EgKU-8 [3] were overexpressed as amino-terminal His6-tagged fusion proteins in Escherichia coli strain BL21(DE3) using pET28a recombinant plasmids prepared according to standard procedures. The expression constructs included the His6 leader sequence followed by the cDNA sequence encoding the corresponding predicted full-length mature EgKU. These were amplified from E. granulosus pepsin/H+-activated protoscolex cDNA [62] using Vent DNA polymerase (New England Biolabs) and specific primers containing restriction sites to allow directional cloning into the pET28a vector. BamHI and HindIII sites were used except for EgKU-6 and EgKU-7 whose coding sequences are cleaved by BamHI; EcoRI was used instead. The integrity of the expression constructs was checked by sequencing. The His6-tagged EgKU fusion proteins were expressed in transformed E. coli grown in LB containing 10 mg/L of kanamycin and induced at 37°C with 0.1 mM isopropyl thiogalactopyranoside. Induction of expression was at late-log-phase (A600 0.6–1.0) during 4 h, in the case of the pairs of paralogs EgKU-1/EgKU-4 and EgKU-3/EgKU-8 that yielded good amounts of soluble recombinant peptides [63]. Expression of EgKU-2, EgKU-5 and the paralogs EgKU-6/EgKU-7, whose recombinants are recovered mostly as inclusion bodies [63], was induced earlier (A600 0.2–0.3) to maximize the yield of the soluble proteins. In all cases, the induced cells were harvested by centrifugation, the pellet was suspended in “lysis” buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM imidazole), and the cells were lyzed by sonication. The lysates were centrifuged (20,000 g for 30 min at 4°C) and the supernatants used to purify the His-tagged fusion proteins using a Ni2+-charged affinity matrix (Ni-NTA, Invitrogen), following the manufacturer’s instructions. The soluble fraction of the bacterial lysates was loaded onto the column equilibrated with lysis buffer, washed with equilibration buffer including 30 mM imidazole, and the recombinant EgKUs eluted with the same buffer containing 250 mM imidazole, that was subsequently dialyzed.

The purity of the EgKUs was checked by SDS-PAGE analysis and the protein concentration was determined with the bicinchoninic acid reagent (BCA, Pierce, USA) using bovine serum albumin as standard, or by A280. The quality of the recombinants was further controlled by: i) confirming the presence of three disulphide bonds (i. e. that the proteins were fully oxidized) through determination of the molecular masses of the EgKUs and their reduced and alkylated derivatives by MALDI-TOF MS, as described by Calvete [64] (S2 Fig); ii) checking that recombinant EgKU-8 reproduced the performance of the native inhibitor towards bovine trypsin and thus, that the recombinant was properly folded and the N-terminal extension contributed by the expression vector did not interfere with enzyme interaction (Ki* was 50 ± 10 pM for recombinant EgKU-8 and 60 ± 13 pM for the native inhibitor [3]); iii) during this study, we also verified that recombinant EgKU-1 reproduced reasonably well the performance of the native inhibitor acting on Kv and ASIC currents from DRG neurons (Figs (Figs33 and and7,7, respectively). Usual yields of the EgKUs recovered as soluble recombinants and used for activity assays were as follows: ~5 mg/L of culture for EgKU-1/EgKU-4 and EgKU-3/EgKU-8; ~300 μg/L for EgKU-7 (~30% of the total); and ~5 μg/L for EgKU-2, EgKU-5 and EgKU-6 (~5% of the total). The purity of the proteins used for detailed activity studies (EgKU-1, EgKU-3, EgKU-4 and EgKU-8) was always > 95%.

Assays of peptidase inhibition by EgKUs

The inhibitory activity of recombinant EgKU-1–EgKU-8 was tested against bovine trypsin (EC, bovine and canine chymotrypsins (EC, and porcine elastase (, essentially as previously described [3]. Bovine enzymes and porcine elastase were obtained from Sigma-Aldrich, whereas the canine peptidase was purified from the pancreas of a dog that had passed away due to an accidental cause, following the procedure of Waritani et al. [65]. According to MEROPS, dogs have two chymotrypsins B (with > 95% identity, encoded by CTRB1 and CTRB2 genes); the fraction we isolated from dog pancreas most likely contained a mixture of both enzymes.The following peptidases were thus assayed (MEROPS——identifiers are indicated in brackets; [9]): from Bos taurus, chymotrypsin A (S01.001) and trypsin 1 (cationic, S01.151); from Sus scrofa, elastase (S01.153); from Canis familiaris, chymotrypsin B (S01.152)).

Prior to inhibition studies, proteolytic activity in enzyme preparations was determined with fluorogenic substrates using initial steady-state rate conditions at 37°C and pH 8.0. Assays (200 μl) were performed in black 96-well microplates (Costar, Corning Life Sciences). Enzymes and substrates were dissolved in 50 mM Tris-HCl, pH 8.0 containing 0.01% Triton X-100 (v/v), and reactions were initiated by the addition of enzyme. The changes in fluorescence intensity, corresponding to the formation of the hydrolysis product 7-amino-4-methylcoumarin (AMC), were registered at excitation and emission wavelengths of 390 and 460 nm, respectively, with a microplate fluorescence reader (FLUOstar* OPTIMA, BMG Labtechnologies). For trypsin activity, the artificial substrate N-t-BOC-Ile-Glu-Gly-Arg-AMC was used; for chymotrypsin, Suc-Ala-Ala-Pro-Phe-AMC; and for elastase, Suc-Ala-Ala-Ala-AMC. Calibration curves using AMC were carried out in each experiment. Initial steady-state rates of substrate hydrolysis were calculated from the linear portion of product (AMC) versus time plots when less than 10% of substrate had been consumed. The substrates and AMC were also obtained from Sigma-Aldrich.

Protein concentrations of enzyme preparations were determined with the BCA reagent using bovine serum albumin as standard; and the active site concentration of trypsin and bovine chymotrypsin A by specific titration with the high affinity inhibitor BPTI. Initially, the active site concentration of canine chymotrypsin could not be estimated because, similar to bovine chymotrypsin B [66], it was not inhibited by BPTI. The enzyme was subsequently titrated with EgKU-3 that inhibits chymotrypsins A and B with high affinity (see results).

The kinetic parameters for substrate and enzyme pairs were calculated from the non-linear fitting to the Michaelis-Menten equation. The values determined with the substrates specified above were: KM = 85 ± 9 μM and kCat = 50 ± 6 s-1 for bovine trypsin; KM = 30 ± 2 μM and kCat = 19 ± 2 s-1 for bovine chymotrypsin A; KM = 39 ± 2 μM and kCat = 6 ± 1 s-1 for canine chymotrypsin B.

For inhibition studies, each of the enzymes was incubated with the purified recombinant EgKUs for 15 min at 37°C prior to the addition of the appropriate fluorogenic substrate, to allow for the equilibration of the enzyme-inhibitor complexes. The substrate concentration (5 μM) was chosen so as to be well below the corresponding KM, as specified above.

To check whether the EgKUs behaved as high affinity inhibitors, the purified recombinants were titrated against active-site titrated bovine trypsin and chymotrypsin, as described by Olson et al. [67].

Peptidase inhibition studies with EgKU-3 and EgKU-4

The activities of EgKU-3 and EgKU-4 were further analyzed by characterizing the kinetics of enzyme inhibition, as previously described for EgKU-8 [3]. All experiments were carried out at least two independent times. Within each experiment, measurements were performed in duplicates. The inhibition and rate constants reported are the average ± standard error of independent experiments.

Tight-binding kinetics

In order to determine the global inhibition constants (KI*) of the recombinant inhibitors towards the assayed serine peptidases, the initial steady-state rates of substrate hydrolysis in the presence of increasing concentrations of each inhibitor were measured after pre-incubation of the enzyme with inhibitor. The inhibition constants were calculated by nonlinear fitting to the Morrison equation for tight binding inhibitors [68,69]:


where KI*app is the apparent global dissociation constant of the enzyme-inhibitor complex, vi is the inhibited steady-state rate, v is the uninhibited rate, [I] is the total inhibitor concentration and [E] is the total enzyme concentration. The true inhibition constants, KI*, were corrected from KI*app according to the Eq 2 for competitive inhibitors:


Slow-binding kinetics

The decrease in the rate of product formation during the first minutes after mixing the enzyme with the inhibitor and substrate (5 μM) was studied for increasing inhibitor concentrations. Progress curves were analyzed using the Eq 3 [70] that describes the slow establishment of equilibrium between the enzyme and the inhibitor according to:


where P is the concentration of AMC produced by hydrolysis of the substrate, vo is the initial rate, vi is the rate once the enzyme-inhibitor equilibrium was reached, and kobs represents the apparent first order rate constant. Computer fitting of progress curves estimated values for vo, vi and kobs.

For Kunitz inhibitors that bind to the enzyme rapidly and reversibly forming an initial “loose” complex EI that isomerizes slowly to the final complex EI*, the reaction mechanism can be represented by Eq 4:

E + I ⇌ EI ⇌ EI*

In this mechanism, the value of the apparent rate constant (kobs) is related to the kinetic constants of the second step, k2 and k-2, and to the equilibrium dissociation constant of the initial loose complex KI (KI = k-1/k1), by Eq 5 [24]:


In the case of EgKU-3 that was found to follow the mechanism represented in Eq 4, the constants k2 and KI were determined from plots of kobs versus [I], by computer fitting to Eq 5. Because k-2 was too small to be accurately estimated from these plots, it was determined with Eq 6 [23] using data from situations where the ratio vi/vo was higher than 0.05:


The values of k-2, k2 and KI thus determined allowed to corroborate the inhibition constant KI*, according to Eq 7:


Data analysis

Computer fitting to non-linear equations was performed using the software Origin version 8 (OriginLab).

Patch-clamp experiments with EgKU-1 and EgKU-4 on DRG neurons

The effect of EgKU-1 and EgKU-4 on voltage-gated (Na+ and K+) and ASIC currents was studied using the whole cell patch-clamp technique in primary cultured rat DRG neurons. The effect of EgKU-3 and EgKU-8 on voltage-gated K+ and ASIC currents was similarly analyzed. α-DTX (kindly donated by Dr. Carlos Cerveñansky from the Unidad de Bioquímica y Proteómica Analíticas, Institut Pasteur de Montevideo/Instituto de Investigaciones Biológicas Clemente Estable, Montevideo, Uruguay) and Bovine Serum Albumin (Sigma Chemicals) were used as positive and negative controls, respectively.

Ethics statement

The study was performed in strict accordance with the Guiding Principles from the Committee on Guide for the Care and Use of Laboratory Animals of the National Research Council of the National Academies of the United States of America and with the regulations of the General Law of Health and of Research in Health Matters of the Ministry of Health of Mexico and Technical Specifications for production, use and Care of laboratory animals (NOM-062-ZOO-1999). The animal protocol was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of the Vice-rectory of Research and Postgraduate Studies of the Autonomous University of Puebla (BUAP/VIEP 2014–236, BUAP/VIEP 2015–273 and BUAP/VIEP 2016–265). All efforts were made to minimize animal suffering and to reduce the number of animals used. The animals were provided by the “Claude Bernard” animal facility of the Autonomous University of Puebla.

Cell culture

The DRG neurons were isolated from the vertebral column of Wistar rats at postnatal ages P5 to P9 without sex distinction, and cultured according to the procedure described by Salceda and coworkers [71]. In brief, the neurons were incubated (30 min at 37°C) in Leibovitz L15 medium (Invitrogen, USA) containing 1.25 mg/ml trypsin and 1.25 mg/ml collagenase (both from Sigma-Aldrich). After the enzymatic treatment, the ganglia were washed 3 times with sterile L15. Cells were mechanically dissociated and plated on 12-mm x 10-mm glass coverslips (Corning, USA), pretreated with poly-D-lysine (Sigma-Aldrich), which were placed onto 35-mm culture dishes (Corning). The isolated cells were allowed to settle and adhere to the coverslips during 4 to 8 h in a humid atmosphere (95% air and 5% CO2, at 37°C) using a CO2 water-jacketed incubator (Nuaire, USA). The plating medium was L15 supplemented with 15.7 mM NaHCO3 (Merck, Mexico), 10% fetal bovine serum, 2.5 μg/ml fungizone (both from Invitrogen), 100 U/mL penicillin (Lakeside, Mexico), and 15.8 mM HEPES (Sigma-Aldrich).

Electrophysiological recording

A coverslip with attached cells was transferred to a 500 μl perfusion chamber mounted on the stage of an inverted phase-contrast microscope (TMS, Nikon Co, Japan). Cells were bath-perfused with extracellular solution by means of a peristaltic pump (Masterflex, L/S Easy-Load II; Cole Parmer, USA). Current recording was carried out at room temperature (23–25°C) using an Axopatch 1D amplifier (Molecular Devices, USA). The cells selected for recording were refringent, they were not adhered to other cells, showed no neurite outgrowth, and had a round soma. Command-pulse generation and data sampling were controlled with the pClamp 9.2 software (Molecular Devices) using a 16-bit data-acquisition system (Digidata 1320, Molecular Devices). Signals were low-pass filtered at 5 kHz and digitized at 10 or 20 kHz, depending on the experiment. Patch pipettes were pulled from borosilicate glass capillaries (TW120-3; WPI, USA) using a Flaming-Brown electrode puller (P-80/-PC; Sutter Instruments, USA). They typically had a resistance of 1 to 3 MΩ when filled with intracellular solution. The series resistance was electronically compensated for by ≈ 80%. Seal and series resistance were continuously monitored to guarantee stable recording conditions.

Experimental protocols and data analysis

To study their effects on voltage-gated K+ and Na+ currents, the EgKUs were ejected under pressure using a microinjector (Baby Bee, USA) from a micropipette positioned in the vicinity of the cell under recording. For the recording of total currents, the cells were perfused with an extracellular solution containing (mM): 140 NaCl, 5.4 KCl, 1.8 CaCl2, 1.2 MgCl2, 10 HEPES, adjusted to pH 7.4 with NaOH. For the recording of ASIC currents, the same solution was employed but HEPES was substituted by MES (10 mM) and pH was adjusted to 6.1. For these experiments, the patch pipettes were filled with an intracellular solution containing (mM): 10 NaCl, 125 KCl, 0.1 CaCl2, 10 EGTA, 5 HEPES, 1 NaGTP, 2 MgATP, adjusted to pH 7.2 with KOH.

Potassium transmembrane currents were elicited by a single-step voltage protocol applying a 120 ms pre-pulse to -100 mV before the 0 mV (800 ms) test pulse with an interval between sweeps of 8 s (holding potencial, Vh = -60 mV). For current-voltage (I/V) analyses, the total potassium current was obtained from a Vh of -45 mV through a series of voltage steps (from -65 to 55 mV, every 8 s) preceded by a conditioning step to -120 mV (200 ms). To record the sustained component of the current (IKDR), the conditioning step was set to -45 mV to inactivate the transient component of the current (IKA). Then IKA was calculated by a trace by trace subtraction between the aforementioned protocols. For this group of experiments, the extracellular solution contained (mM): 10 KCl, 1.8 CaCl2, 1.2 MgCl2, 0.3 CdCl2, 130 choline chloride, 10 HEPES, adjusted to pH 7.4; and the intracellular solution (mM): 50 KCl, 60 choline chloride, 0.1 CaCl2, 40 KF, 10 EGTA, 5 HEPES, adjusted to pH 7.2.

The EgKUs were always co-applied with the single-step voltage protocol. I/V analyses were carried out after the effect had stabilized. Two parameters were used to evaluate the effect: the maximum amplitude of the current and the amplitude of the current in the last 10 ms of the test pulse.

Sodium currents were evoked by a 40 ms depolarization step to a membrane potential of -10 mV from a Vh of -100 mV, with an interval between sweeps of 8 s. Three parameters were used to evaluate the effect of the EgKUs: the maximum amplitude of the current (INamax), the time constant of the current inactivation (τh), as derived from an exponential fit, and the ratio between the current amplitude at the end of the voltage pulse and INamax (INaend/INamax), which gives an estimate of the probability for the channels not being inactivated at the end of the voltage pulse. Cells were bath-perfused with an extracellular solution containing (mM): 20 NaCl, 70 choline chloride, 1.8 CaCl2, 1 MgCl2, 10 HEPES, 45 tetraethylammonium chloride, 10 4-aminopyridine adjusted to pH 7.4. Pipettes were filled with the following intracellular solution (mM): 10 NaCl, 30 CsCl, 100 CsF, 5 HEPES, 8 EGTA, 10 tetraethylammonium chloride, adjusted to pH 7.2.

The ASIC currents were generated by a fast (about 100 ms) pH change from 7.4 to 6.1, by shifting one of the three outlets of a fast change perfusion system (SF-77B, Warner Inst, USA) while keeping the cell at a Vh of -60 mV (the extracellular solution of pH 6.1 contained MES, pK = 6.15, instead of HEPES, pK = 7.55). The pH was kept at 6.1 during 5 s. The time between pH changes was 1 min to guarantee that the ASIC current was completely recovered from desensitization. The transient receptor potential vanilloid 1 (TRPV1) antagonist capsazepine (10 μM) (Sigma-Aldrich) was added to the extracellular solution (pH 6.1) to prevent activation of the TRPV1 receptor present in DRG neurons. The ASIC currents were characterized by the maximum peak amplitude (Ipeak), the desensitization time-constant (τdes, determined by fitting the decay phase of the current with a single exponential function), and the amplitude at the steady state (Iend) measured at the last 100 ms of the acid pulse. The EgKUs were applied 20 s before and during the acid pulse (sustained application). At least two control responses were recorded before any experimental manipulation. The pH of the perfusion solution was checked not to be affected by addition of the inhibitors.

Concentration-response data were fitted with the function Y = min + (max-min)/[1+(x/EC50)H], where Y is the effect of the inhibitor, x is the concentration, max and min are the maximum and minimum effects, EC50 is the concentration at which 50% of the effect is obtained and H is the Hill coefficient. Experimental data are presented as the mean ± standard error. To define the statistical significance, a paired Student’s t-test was used and P ≤ 0.05 was considered as significant, when comparing the effects in the presence and absence of inhibitor.

Analysis of EgKUs in parasite secretions

Fresh hydatid fluid was recovered under aseptic conditions from individual fertile cysts of the G1 genotype (E. granulosus sensu stricto), present in the lungs of naturally infected bovines in Uruguay, and kept at -70°C. Cysts were collected during the routine work of local abattoirs in Montevideo. Cyst fluid was analyzed by MALDI-TOF MS using a Voyager DE-PRO spectrometer (Applied Biosystems). The sample was concentrated by vacuum drying, desalted using C18 reverse phase micro-columns (OMIX Pipette tips, Varian) and eluted with matrix solution (α-cyano-4-hydroxycinnamic acid in 0.2% trifluoroacetic acid in 50% v/v acetonitrile-H2O) directly on the MALDI sample plate.

An aliquot of adult worm in vitro secretions was kindly provided by MSc Noelia Morel; the culture was carried out in our department in the context of a project to develop a copro-ELISA kit for canine echinococcosis [72]. Worm secretions were prepared essentially as described by Casaravilla and coworkers [73]. In brief, adult worms were recovered from the intestine of dogs experimentally infected with E. granulosus (euthanized 30 days post-infection), washed with sterile phosphate-buffered saline (PBS) and cultured in RPMI containing 105 UI/L penicillin, 100 mg/L streptomycin and 250 μg/L amphotericin B, at 37°C in 5% CO2.The supernatant was collected every 8 h for two days and kept at -70°C. An aliquot was clarified by centrifugation at 10,000 g and analyzed as described for cyst fluid.

EgKU-3 and EgKU-8 were affinity purified from hydatid fluid and worm secretions using chymotrypsin A agarose (Sigma-Aldrich, USA) and analyzed by MALDI-TOF MS (4800 MALDI TOF/TOF analyzer, ABi Sciex) to confirm their presence in the parasite secretions. The resin containing agarose-bound bovine chymotrypsin A (1.5 mg) was rehydrated and equilibrated with 10 mM Tris-HCl buffer pH 7.5. Cyst fluid (1 ml) was incubated in batch with the resin during 10 min at 20°C; the resin was then washed 3 times with 10 mM Tris-HCl pH 7.5, and the EgKUs were eluted by incubation during 10 min with 25 μl of 0.2% trifluoroacetic acid. After centrifugation, 1 μl of the eluate was applied directly on the MALDI sample plate with 1 μl of the matrix solution (α-cyano-4-hydroxycinnamic acid in 0.1% trifluoroacetic acid in 60% v/v acetonitrile-H2O). Mass spectra were acquired in positive ion linear mode and externally calibrated using protein standards (Applied Biosystems). The EgKUs were purified from adult worm secretions with the same protocol, using 3 ml of the sample previously concentrated by vacuum drying.

EgKU-3 and EgKU-8 purified from cyst fluid were characterized by peptide mass fingerprinting; the eluate was reduced and alkylated with iodoacetamide prior to treatment with trypsin (Sequencing-grade, Promega). The sequence of selected peptides was confirmed by collision-induced dissociation MS/MS experiments, as previously described [3].

Structural modeling

The full-length mature sequences of EgKU-1, EgKU-3, EgKU-4 and EgKU-8 were used to compute structural models using the i-Tasser server [74]. The C-scores for all models were higher than 0 (EgKU-1 = 0.93; EgKU-3 = 1.23; EgKU-4 = 0.05; EgKU-8 = 1.13). Typical C-scores range from [–2, 5], with higher scores meaning more reliable models. α-DTX did not stand amongst the top 10 threading templates used by i-Tasser, which allowed direct comparisons with no circularities in the analyses. Electrostatic properties were calculated at pH 7.4 with the Adaptive Poisson Boltzmann Solver [75] for the best i-Tasser models as well as for the crystal structure of α-DTX (PDB access code 1DTX). Next, per-residue solvent-accessibility was computed with the “areaimol” program from the CCP4 suite [76]. Basic, acidic and aromatic residues with solvent accessibilities above 40 Å2 were displayed in van der Waals representation onto a cartoon backbone of the model, using VMD [77]. Electrostatic molecular surface representations were produced and rendered with PyMol (

Phylogenetic analysis

The Kunitz domains from EgKU-1-EgKU-8 as well as from eight additional monodomain Kunitz proteins identified in the E. granulosus genome sequences (encoded by EgrG_001136500, EgrG_001136600, EgrG_001136800, EgrG_001137000, EgrG_001137200, EgrG_1137300 and EgrG_1137400 from the genome produced at the Wellcome Trust Sanger Centre [5]; and the protein EGR_07242, from the genome produced at the Chinese National Human Genome Center [35]) were multiply aligned with Mafft [78] (L-insi option). Eleven close paralogs (some of which are putative orthologs) of the E. granulosus monodomain Kunitz proteins identified within genomic and transcriptomic data from the T. solium genome project [5] (encoded by genes TsM_000321400, TsM_000410200, TsM_000513000, TsM_000576900, TsM_000647700, TsM_000724900, TsM_001022000, TsM_001027800, TsM_001085400, TsM_001162600 retrieved from GeneDB—; and the EST EL746785 retrieved as TSE0004567 from PartiGeneDB— were added to this set. A group of five functionally annotated sequences from other Lophotrochozoa was also included: FhKTM (UniProt Q9TXD3; [42]) and FhKT1 [43], from F. hepatica; SjKI-1 (Sjp_0020270) from S. japonicum [37]; SmKI-1 (Smp_147730) from S. mansoni [38] and Conkunitzin-S1 (UniProt P0C1X2) from the mollusk Conus striatus [14]. The alignment of these 32 proteins was used as input for MrBayes [79] for a Bayesian phylogenetic reconstruction using the Poisson substitution model in a run of 1,000,000 generations, discarding the first 100,000 for summarizing results (mcmc ngen = 1000000; sump burnin = 100000; sumt burnin = 100000). The short sequence length of the Kunitz domain prevents robust and reliable identification of many branching events. Furthermore, saturation events are guaranteed to occur and may be difficult to pinpoint. The tree is one of the best we can get with current methods. In fact, maximum likelihood reconstruction with 1000 bootstraps provided poorer support for branching events. The final dendrogram was visualized and rendered in FigTree ( It should be noted that the proteins encoded by EgrG_001136500, EgrG_001136600, EgrG_001137000, EgrG_001137200, and EgrG_001137400 [5] correspond to EGR_08721, EGR_08720, EGR_08716, EGR_9006, and EGR_9007 [35], respectively.

Supporting information

S1 Dataset

Peptide mass fingerprinting of chymotrypsin affinity purified proteins from hydatid cyst fluid.

(A) MALDI-TOF mass spectrum of peptides generated by tryptic digestion of cyst fluid proteins recovered after an affinity purification step using chymotrypsin A. The signals that can be assigned to peptides derived from EgKU-8 or EgKU-3 are indicated (refer to Fig 10C). The signal of m/z 1074.54 corresponds to EgKU-8 sequence LPLDPGFcR (theoretical MH+ 1074.53); the signal of m/z 1164.52 to EgKU-3 sequence EQcELLcGR (theoretical MH+ 1164.51); and the signal of m/z 1491.63 to EgKU-8 sequence WGFHQESGEcVR (theoretical MH+ 1491.64); c indicates carbamidomethyl Cys. (B) MS/MS analysis of tryptic peptides from EgKU-3 and EgKU-8. The list of theoretical m/z values of fragment ions is shown for each sequence, and the ions detected in the MS/MS spectra are highlighted in bold. a and b ions correspond to N-terminal fragments, and y ions to C-terminal fragments, according to the accepted nomenclature [80,81].


S1 Fig

Studies with EgKU-1 and EgKU-4 on isolated Na+ currents from DRG neurons.

Representative traces showing that the sustained (25 s) perfusion of 100 nM recombinant EgKU-1 (A) or EgKU-4 (C) does not block voltage-activated sodium channels (Nav). (B) and (D) are the current-voltage relationships of the peak Na+ current from the traces in (A) and (C), respectively. The black traces correspond to control conditions and the gray ones after EgKU perfusion.


S2 Fig

MALDI-TOF MS analysis of recombinant EgKUs: results for EgKU-3.

Mass spectrometry analyses of free thiols and disulphide bonds were carried out as described in Calvete [64]. Three samples were examined for each purified recombinant EgKU: i) untreated (A); ii) denatured with guanidinium hydrochloride and treated with iodoacetamide (IA), to assess the presence of free thiols (B); iii) denatured with guanidinium hydrochloride, reduced with DTT and treated with IA to confirm the presence of 6 Cys residues (C). Predicted MH+ values for EgKU-3 are as follows: untreated EgKU-3 (3 Cys-Cys) = 9819 Da; reduced and alkylated EgKU-3 = 10167 Da.


S1 Table

cDNA and gene sequence data from EgKU-1-EgKU-8.



We thank José Pacheco (Departamento de Patología, Facultad de Veterinaria, UdelaR) for providing us with the dog pancreas; Madelón Portela (Unidad de Bioquímica y Proteómica Analíticas, Institut Pasteur de Montevideo) for help with MS analysis; Héctor Romero (Departamento de Ecología y Evolución, Facultad de Ciencias, UdelaR) for advice on phylogenetic issues; and Carla Palermo (Departamento de Biociencias, Facultad de Química, UdelaR) and Gonzalo Pérez (Instituto de Química Biológica, Facultad de Ciencias, UdelaR) for help with the preparation of recombinant EgKU-1/EgKU-4 and EgKU-5, respectively.

Funding Statement

The work was supported by grants from Comisión Sectorial de Investigación Científica (CSIC-I+D 2008 and CSIC-Grupos I+D 2014, UdelaR); and by fellowships from Agencia Nacional de Investigación e Innovación (ANII, Uruguay) to LP (BE_INI_2008_462), and MF (BE_POS_2009_1206 and BE_POS_2011_3635); and from Comisión Académica de Posgrado (UdelaR) to MF. MF, MG, RD, BA and CF also acknowledge PEDECIBA and ANII (SNI) for support. Work at the laboratory of ES was funded by grants from Consejo Nacional de Ciencia y Tecnología de México (CONACyT, 167052, 169835 and 229866), Vicerrectoría de Investigación y Estudios de Posgrado (BUAP/VIEP 2014-236, BUAP/VIEP 2015-273 and BUAP/VIEP 2016-265) and PROFOCIE 2014. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files. The sequences of the E. granulosus Kunitz proteins whose characterization is described were derived from an EST project and previously published [3]. GenBank accession numbers of the corresponding full-length coding cDNAs are as follows: KY438959 for EgKU-1; KY438960 for EgKU-2; KY438963 for EgKU-3; KY438961 for EgKU-4; KY438965 for EgKU-5; KY438962 for EgKU-6; KY438964 for EgKU-7; and FJ031017 for EgKU-8. Refer to S1 Table for further information on cDNA and gene sequence data related to EgKU-1-EgKU-8.


1. Budke CM, White AC Jr., Garcia HH (2009) Zoonotic larval cestode infections: neglected, neglected tropical diseases? PLoS Negl Trop Dis 3: e319 doi: 10.1371/journal.pntd.0000319 [PMC free article] [PubMed]
2. Thompson RC, Jenkins DJ (2014) Echinococcus as a model system: biology and epidemiology. Int J Parasitol 44: 865–877. doi: 10.1016/j.ijpara.2014.07.005 [PubMed]
3. González S, Fló M, Margenat M, Durán R, González-Sapienza G, et al. (2009) A family of diverse Kunitz inhibitors from Echinococcus granulosus potentially involved in host-parasite cross-talk. PLoS One 4: e7009 doi: 10.1371/journal.pone.0007009 [PMC free article] [PubMed]
4. Parkinson J, Wasmuth J, Salinas G, Bizarro CV, Sanford C, et al. (2012) A transcriptomic analysis of Echinococcus granulosus larval stages: implications for parasite biology and host adaptation. PLoS Negl Trop Dis 6: e1897 doi: 10.1371/journal.pntd.0001897 [PMC free article] [PubMed]
5. Tsai IJ, Zarowiecki M, Holroyd N, Garciarrubio A, Sanchez-Flores A, et al. (2013) The genomes of four tapeworm species reveal adaptations to parasitism. Nature 496: 57–63. doi: 10.1038/nature12031 [PMC free article] [PubMed]
6. Cardona GA, Carmena D (2013) A review of the global prevalence, molecular epidemiology and economics of cystic echinococcosis in production animals. Vet Parasitol 192: 10–32. doi: 10.1016/j.vetpar.2012.09.027 [PubMed]
7. Thompson RCA (1995) Biology and systematics of Echinococcus In: Thompson RCA, Lymbery AJ, editors. Echinococcus and Hydatid Disease. Wallingford: CAB International; pp. 1–50.
8. Smyth JD, McManus DP (1989) The cestodes: general considerations The Physiology and Biochemistry of Cestodes. Cambridge, USA: Cambridge University Press.
9. Rawlings ND, Barrett AJ, Finn RD (2016) Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res 44: D343–350. doi: 10.1093/nar/gkv1118 [PMC free article] [PubMed]
10. Laskowski M Jr., Kato I (1980) Protein inhibitors of proteinases. Annu Rev Biochem 49: 593–626. doi: 10.1146/ [PubMed]
11. Fry BG (2005) From genome to "venome": molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res 15: 403–420. doi: 10.1101/gr.3228405 [PubMed]
12. Schweitz H, Bruhn T, Guillemare E, Moinier D, Lancelin JM, et al. (1995) Kalicludines and kaliseptine. Two different classes of sea anemone toxins for voltage sensitive K+ channels. J Biol Chem 270: 25121–25126. [PubMed]
13. Andreev YA, Kozlov SA, Koshelev SG, Ivanova EA, Monastyrnaya MM, et al. (2008) Analgesic compound from sea anemone Heteractis crispa is the first polypeptide inhibitor of vanilloid receptor 1 (TRPV1). J Biol Chem 283: 23914–23921. doi: 10.1074/jbc.M800776200 [PMC free article] [PubMed]
14. Bayrhuber M, Vijayan V, Ferber M, Graf R, Korukottu J, et al. (2005) Conkunitzin-S1 is the first member of a new Kunitz-type neurotoxin family. Structural and functional characterization. J Biol Chem 280: 23766–23770. doi: 10.1074/jbc.C500064200 [PubMed]
15. Yuan CH, He QY, Peng K, Diao JB, Jiang LP, et al. (2008) Discovery of a distinct superfamily of Kunitz-type toxin (KTT) from tarantulas. PLoS One 3: e3414 doi: 10.1371/journal.pone.0003414 [PMC free article] [PubMed]
16. Chen Z, Luo F, Feng J, Yang W, Zeng D, et al. (2013) Genomic and structural characterization of Kunitz-type peptide LmKTT-1a highlights diversity and evolution of scorpion potassium channel toxins. PLoS One 8: e60201 doi: 10.1371/journal.pone.0060201 [PMC free article] [PubMed]
17. Chen ZY, Hu YT, Yang WS, He YW, Feng J, et al. (2012) Hg1, novel peptide inhibitor specific for Kv1.3 channels from first scorpion Kunitz-type potassium channel toxin family. J Biol Chem 287: 13813–13821. doi: 10.1074/jbc.M112.343996 [PMC free article] [PubMed]
18. Schwarz A, Cabezas-Cruz A, Kopecky J, Valdes JJ (2014) Understanding the evolutionary structural variability and target specificity of tick salivary Kunitz peptides using next generation transcriptome data. BMC Evol Biol 14: 4 doi: 10.1186/1471-2148-14-4 [PMC free article] [PubMed]
19. Dai SX, Zhang AD, Huang JF (2012) Evolution, expansion and expression of the Kunitz/BPTI gene family associated with long-term blood feeding in Ixodes Scapularis. BMC Evol Biol 12: 4 doi: 10.1186/1471-2148-12-4 [PMC free article] [PubMed]
20. Cappello M, Vlasuk GP, Bergum PW, Huang S, Hotez PJ (1995) Ancylostoma caninum anticoagulant peptide: a hookworm-derived inhibitor of human coagulation factor Xa. Proc Natl Acad Sci U S A 92: 6152–6156. [PubMed]
21. Harvey AL (2001) Twenty years of dendrotoxins. Toxicon 39: 15–26. [PubMed]
22. Báez A, Salceda E, Fló M, Graña M, Fernández C, et al. (2015) alpha-Dendrotoxin inhibits the ASIC current in dorsal root ganglion neurons from rat. Neurosci Lett 606: 42–47. doi: 10.1016/j.neulet.2015.08.034 [PubMed]
23. Morrison JF (1982) The slow-binding and slow, tight-binding inhibition of enzyme-catalysed reactions. Trends Biochem Sci 7: 102–105.
24. Morrison JF, Walsh CT (1988) The behavior and significance of slow-binding enzyme inhibitors. Adv Enzymol Relat Areas Mol Biol 61: 201–301. [PubMed]
25. Lima CA, Torquato RJ, Sasaki SD, Justo GZ, Tanaka AS (2010) Biochemical characterization of a Kunitz type inhibitor similar to dendrotoxins produced by Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) hemocytes. Vet Parasitol 167: 279–287. doi: 10.1016/j.vetpar.2009.09.030 [PubMed]
26. Bhattacharyya A, Mazumdar S, Leighton SM, Babu CR (2006) A Kunitz proteinase inhibitor from Archidendron ellipticum seeds: purification, characterization, and kinetic properties. Phytochemistry 67: 232–241. doi: 10.1016/j.phytochem.2005.11.010 [PubMed]
27. Vincent JP, Lazdunski M (1972) Trypsin-pancreatic trypsin inhibitor association. Dynamics of the interaction and role of disulfide bridges. Biochemistry 11: 2967–2977. [PubMed]
28. Hudaky P, Kaslik G, Venekei I, Graf L (1999) The differential specificity of chymotrypsin A and B is determined by amino acid 226. Eur J Biochem 259: 528–533. [PubMed]
29. Hall A, Stow J, Sorensen R, Dolly JO, Owen D (1994) Blockade by dendrotoxin homologues of voltage-dependent K+ currents in cultured sensory neurones from neonatal rats. Br J Pharmacol 113: 959–967. [PMC free article] [PubMed]
30. Escoubas P, De Weille JR, Lecoq A, Diochot S, Waldmann R, et al. (2000) Isolation of a tarantula toxin specific for a class of proton-gated Na+ channels. J Biol Chem 275: 25116–25121. doi: 10.1074/jbc.M003643200 [PubMed]
31. Diochot S, Salinas M, Baron A, Escoubas P, Lazdunski M (2007) Peptide inhibitors of acid-sensing ion channels. Toxicon 49: 271–284. doi: 10.1016/j.toxicon.2006.09.026 [PubMed]
32. Lingueglia E, de Weille JR, Bassilana F, Heurteaux C, Sakai H, et al. (1997) A modulatory subunit of acid sensing ion channels in brain and dorsal root ganglion cells. J Biol Chem 272: 29778–29783. [PubMed]
33. Chagot B, Escoubas P, Diochot S, Bernard C, Lazdunski M, et al. (2005) Solution structure of APETx2, a specific peptide inhibitor of ASIC3 proton-gated channels. Protein Sci 14: 2003–2010. doi: 10.1110/ps.051378905 [PubMed]
34. Gasparini S, Danse JM, Lecoq A, Pinkasfeld S, Zinn-Justin S, et al. (1998) Delineation of the functional site of alpha-dendrotoxin. The functional topographies of dendrotoxins are different but share a conserved core with those of other Kv1 potassium channel-blocking toxins. J Biol Chem 273: 25393–25403. [PubMed]
35. Zheng H, Zhang W, Zhang L, Zhang Z, Li J, et al. (2013) The genome of the hydatid tapeworm Echinococcus granulosus. Nat Genet 45: 1168–1175. doi: 10.1038/ng.2757 [PubMed]
36. Ranasinghe SL, Fischer K, Zhang W, Gobert GN, McManus DP (2015) Cloning and characterization of two potent Kunitz type protease inhibitors from Echinococcus granulosus. PLoS Negl Trop Dis 9: e0004268 doi: 10.1371/journal.pntd.0004268 [PMC free article] [PubMed]
37. Ranasinghe SL, Fischer K, Gobert GN, McManus DP (2015) A novel coagulation inhibitor from Schistosoma japonicum. Parasitology 142: 1663–1672. doi: 10.1017/S0031182015001328 [PubMed]
38. Ranasinghe SL, Fischer K, Gobert GN, McManus DP (2015) Functional expression of a novel Kunitz type protease inhibitor from the human blood fluke Schistosoma mansoni. Parasit Vectors 8: 408 doi: 10.1186/s13071-015-1022-z [PMC free article] [PubMed]
39. Harvey AL, Robertson B (2004) Dendrotoxins: structure-activity relationships and effects on potassium ion channels. Curr Med Chem 11: 3065–3072. [PubMed]
40. Katoh E, Nishio H, Inui T, Nishiuchi Y, Kimura T, et al. (2000) Structural basis for the biological activity of dendrotoxin-I, a potent potassium channel blocker. Biopolymers 54: 44–57. doi: 10.1002/(SICI)1097-0282(200007)54:1<44::AID-BIP50>3.0.CO;2-Z [PubMed]
41. Smith LA, Reid PF, Wang FC, Parcej DN, Schmidt JJ, et al. (1997) Site-directed mutagenesis of dendrotoxin K reveals amino acids critical for its interaction with neuronal K+ channels. Biochemistry 36: 7690–7696. doi: 10.1021/bi963105g [PubMed]
42. Bozas SE, Panaccio M, Creaney J, Dosen M, Parsons JC, et al. (1995) Characterisation of a novel Kunitz-type molecule from the trematode Fasciola hepatica. Mol Biochem Parasitol 74: 19–29. [PubMed]
43. Smith D, Tikhonova IG, Jewhurst HL, Drysdale OC, Dvorak J, et al. (2016) Unexpected activity of a novel Kunitz-type inhibitor: Inhibition of cysteine proteases but not serine proteases. J Biol Chem 291: 19220–19234. doi: 10.1074/jbc.M116.724344 [PMC free article] [PubMed]
44. Aziz A, Zhang W, Li J, Loukas A, McManus DP, et al. (2011) Proteomic characterisation of Echinococcus granulosus hydatid cyst fluid from sheep, cattle and humans. J Proteomics 74: 1560–1572. doi: 10.1016/j.jprot.2011.02.021 [PubMed]
45. Huang F, Dang Z, Suzuki Y, Horiuchi T, Yagi K, et al. (2016) Analysis on gene expression profile in oncospheres and early stage metacestodes from Echinococcus multilocularis. PLoS Negl Trop Dis 10: e0004634 doi: 10.1371/journal.pntd.0004634 [PMC free article] [PubMed]
46. Smyth JD (1964) Observations on the scolex of Echinococcus granulosus, with special reference to the occurrence and cytochemistry of secretory cells in the rostellum. Parasitology 54: 515–526. [PubMed]
47. McDermott JR, Bartram RE, Knight PA, Miller HR, Garrod DR, et al. (2003) Mast cells disrupt epithelial barrier function during enteric nematode infection. Proc Natl Acad Sci U S A 100: 7761–7766. doi: 10.1073/pnas.1231488100 [PubMed]
48. Antalis TM, Buzza MS, Hodge KM, Hooper JD, Netzel-Arnett S (2010) The cutting edge: membrane-anchored serine protease activities in the pericellular microenvironment. Biochem J 428: 325–346. doi: 10.1042/BJ20100046 [PMC free article] [PubMed]
49. Buzza MS, Netzel-Arnett S, Shea-Donohue T, Zhao A, Lin CY, et al. (2010) Membrane-anchored serine protease matriptase regulates epithelial barrier formation and permeability in the intestine. Proc Natl Acad Sci U S A 107: 4200–4205. doi: 10.1073/pnas.0903923107 [PubMed]
50. Gieseler F, Ungefroren H, Settmacher U, Hollenberg MD, Kaufmann R (2013) Proteinase-activated receptors (PARs)—focus on receptor-receptor-interactions and their physiological and pathophysiological impact. Cell Commun Signal 11: 86 doi: 10.1186/1478-811X-11-86 [PMC free article] [PubMed]
51. Heitzmann D, Warth R (2008) Physiology and pathophysiology of potassium channels in gastrointestinal epithelia. Physiol Rev 88: 1119–1182. doi: 10.1152/physrev.00020.2007 [PubMed]
52. Holzer P (2015) Acid-sensing ion channels in gastrointestinal function. Neuropharmacology 94: 72–79. doi: 10.1016/j.neuropharm.2014.12.009 [PMC free article] [PubMed]
53. Warth R, Barhanin J (2003) Function of K+ channels in the intestinal epithelium. J Membr Biol 193: 67–78. doi: 10.1007/s00232-002-2001-9 [PubMed]
54. Matzner N, Zemtsova IM, Nguyen TX, Duszenko M, Shumilina E, et al. (2008) Ion channels modulating mouse dendritic cell functions. J Immunol 181: 6803–6809. [PubMed]
55. Kong X, Tang X, Du W, Tong J, Yan Y, et al. (2013) Extracellular acidosis modulates the endocytosis and maturation of macrophages. Cell Immunol 281: 44–50. doi: 10.1016/j.cellimm.2012.12.009 [PubMed]
56. Wang Y, Zhou H, Shen Y, Wu W, Liu H, et al. (2015) Impairment of dendritic cell function and induction of CD4(+)CD25(+)Foxp3(+) T cells by excretory-secretory products: a potential mechanism of immune evasion adopted by Echinococcus granulosus. BMC Immunol 16: 44 doi: 10.1186/s12865-015-0110-3 [PMC free article] [PubMed]
57. Falcon CR, Masih D, Gatti G, Sanchez MC, Motran CC, et al. (2014) Fasciola hepatica Kunitz type molecule decreases dendritic cell activation and their ability to induce inflammatory responses. PLoS One 9: e114505 doi: 10.1371/journal.pone.0114505 [PMC free article] [PubMed]
58. Howe KL, Bolt BJ, Cain S, Chan J, Chen WJ, et al. (2016) WormBase 2016: expanding to enable helminth genomic research. Nucleic Acids Res 44: D774–780. doi: 10.1093/nar/gkv1217 [PMC free article] [PubMed]
59. Wang S, Wei W, Cai X (2015) Genome-wide analysis of excretory/secretory proteins in Echinococcus multilocularis: insights into functional characteristics of the tapeworm secretome. Parasit Vectors 8: 666 doi: 10.1186/s13071-015-1282-7 [PMC free article] [PubMed]
60. Gomez S, Adalid-Peralta L, Palafox-Fonseca H, Cantu-Robles VA, Soberon X, et al. (2015) Genome analysis of Excretory/Secretory proteins in Taenia solium reveals their Abundance of Antigenic Regions (AAR). Sci Rep 5: 9683 doi: 10.1038/srep09683 [PMC free article] [PubMed]
61. Protasio AV, Tsai IJ, Babbage A, Nichol S, Hunt M, et al. (2012) A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni. PLoS Negl Trop Dis 6: e1455 doi: 10.1371/journal.pntd.0001455 [PMC free article] [PubMed]
62. Fernández C, Gregory WF, Loke P, Maizels RM (2002) Full-length-enriched cDNA libraries from Echinococcus granulosus contain separate populations of oligo-capped and trans-spliced transcripts and a high level of predicted signal peptide sequences. Mol Biochem Parasitol 122: 171–180. [PubMed]
63. Salinas G, Pellizza L, Margenat M, Fló M, Fernández C (2011) Tuned Escherichia coli as a host for the expression of disulfide-rich proteins. Biotechnol J 6: 686–699. doi: 10.1002/biot.201000335 [PubMed]
64. Calvete JJ (2007) Determinación del número de grupos sulfidrilo y de enlaces disulfuro mediante espectrometría de masas. Proteómica 00:21–28. Available from:
65. Waritani T, Okuno Y, Ashida Y, Tsuchiya R, Kobayashi K, et al. (2001) Development and characterization of monoclonal antibodies against canine trypsin. Vet Immunol Immunopathol 80: 333–338. [PubMed]
66. Wu FC, Laskowski M (1955) Action of the naturally occurring trypsin inhibitors against chymotrypsins alpha and beta. J Biol Chem 213: 609–619. [PubMed]
67. Olson ST, Bock PE, Kvassman J, Shore JD, Lawrence DA, et al. (1995) Role of the catalytic serine in the interactions of serine proteinases with protein inhibitors of the serpin family. Contribution of a covalent interaction to the binding energy of serpin-proteinase complexes. J Biol Chem 270: 30007–30017. [PubMed]
68. Morrison JF (1969) Kinetics of the reversible inhibition of enzyme-catalysed reactions by tight-binding inhibitors. Biochim Biophys Acta 185: 269–286. [PubMed]
69. Greco WR, Hakala MT (1979) Evaluation of methods for estimating the dissociation constant of tight binding enzyme inhibitors. J Biol Chem 254: 12104–12109. [PubMed]
70. Williams JW, Morrison JF (1979) The kinetics of reversible tight-binding inhibition. Methods Enzymol 63: 437–467. [PubMed]
71. Salceda E, Garateix A, Soto E (2002) The sea anemone toxins BgII and BgIII prolong the inactivation time course of the tetrodotoxin-sensitive sodium current in rat dorsal root ganglion neurons. J Pharmacol Exp Ther 303: 1067–1074. doi: 10.1124/jpet.102.038570 [PubMed]
72. Morel N, Lassabe G, Elola S, Bondad M, Herrera S, et al. (2013) A monoclonal antibody-based copro-ELISA kit for canine echinococcosis to support the PAHO effort for hydatid disease control in South America. PLoS Negl Trop Dis 7: e1967 doi: 10.1371/journal.pntd.0001967 [PMC free article] [PubMed]
73. Casaravilla C, Malgor R, Rossi A, Sakai H, Nonaka N, et al. (2005) Production and characterization of monoclonal antibodies against excretory/secretory products of adult Echinococcus granulosus, and their application to coproantigen detection. Parasitol Int 54: 43–49. doi: 10.1016/j.parint.2004.08.006 [PubMed]
74. Yang J, Zhang Y (2015) I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res 43: W174–181. doi: 10.1093/nar/gkv342 [PMC free article] [PubMed]
75. Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA (2001) Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci U S A 98: 10037–10041. doi: 10.1073/pnas.181342398 [PubMed]
76. Collaborative Computational Project, Number 4 (1994) The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr 50: 760–763. doi: 10.1107/S0907444994003112 [PubMed]
77. Hsin J, Arkhipov A, Yin Y, Stone JE, Schulten K (2008) Using VMD: an introductory tutorial. Curr Protoc Bioinformatics Chapter 5: Unit 5 7. [PMC free article] [PubMed]
78. Katoh K, Toh H (2008) Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform 9: 286–298. doi: 10.1093/bib/bbn013 [PubMed]
79. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. [PubMed]
80. Roepstorff P, Fohlman J (1984) Biological Mass Spectrometry 11: 601. [PubMed]
81. Johnson RS, Martin SA, Biemann K, Stults JT, Watson JT (1987) Novel fragmentation process of peptides by collision-induced decomposition in a tandem mass spectrometer: differentiation of leucine and isoleucine. Anal Chem 59: 2621–2625. [PubMed]
82. Beitz E (2000) TEXshade: shading and labeling of multiple sequence alignments using LATEX2 epsilon. Bioinformatics 16: 135–139. [PubMed]

Articles from PLoS Pathogens are provided here courtesy of Public Library of Science