Open-access Glycosylation of β1 subunit plays a pivotal role in the toxin sensitivity and activation of BK channels

Abstract

Background:  The accessory β1 subunits, regulating the pharmacological and biophysical properties of BK channels, always undergo post-translational modifications, especially glycosylation. To date, it remains elusive whether the glycosylation contributes to the regulation of BK channels by β1 subunits.

Methods:  Herein, we combined the electrophysiological approach with molecular mutations and biochemical manipulation to investigate the function roles of N-glycosylation in β1 subunits.

Results:  The results show that deglycosylation of β1 subunits through double-site mutations (β1 N80A/N142A or β1 N80Q/N142Q) could significantly increase the inhibitory potency of iberiotoxin, a specific BK channel blocker. The deglycosylated channels also have a different sensitivity to martentoxin, another BK channel modulator with some remarkable effects as reported before. On the contrary to enhancing effects of martentoxin on glycosylated BK channels under the presence of cytoplasmic Ca2+, deglycosylated channels were not affected by the toxin. However, the deglycosylated channels were surprisingly inhibited by martentoxin under the absence of cytoplasmic Ca2+, while the glycosylated channels were not inhibited under this same condition. In addition, wild type BK (α+β1) channels treated with PNGase F also showed the same trend of pharmacological results to the mutants. Similar to this modulation of glycosylation on BK channel pharmacology, the deglycosylated forms of the channels were activated at a faster speed than the glycosylated ones. However, the V1/2 and slope were not changed by the glycosylation.

Conclusion:  The present study reveals that glycosylation is an indispensable determinant of the modulation of β1-subunit on BK channel pharmacology and its activation. The loss of glycosylation of β1 subunits could lead to the dysfunction of BK channel, resulting in a pathological state.

Keywords: BK channels; Glycosylation; β1-subunit; Toxin sensitivity; Kinetic property

Background

BK channels (voltage-dependent large-conductance Ca2+-activated K+ channels),also known as Maxi-K channels, distributed in both excitable and non-excitable cells and are considered as key participants in a variety of physiological functions, including regulating smooth muscle tone[1-3], neuronal firing [4, 5], endocrine cell secretion [6, 7], cell proliferation [8, 9] and migration [10, 11]. Functional BK channels are a tetramer of four pore-forming α subunits encoded by a single gene Slowpoke (Slo) [12]. Different from the close homology of voltage-gated K+ (Kv) channel, the α subunit of BK channels possesses additional hydrophobic segments including a transmembrane helix (S0) which places the N terminus on the extracellular side of the plasma membrane [13, 14] and a long cytosolic C-terminal (S7-S10) where putative Ca2+-binding sites reside [14, 15]. Owing to the various tissue-specific β subunits, BK channels possess a rather complex diversity of subtype family. Four β subunit types have been characterized so far, β1 from smooth muscle[16], β2 from rat chromaffin cells[17, 18], β3 from testis [19, 20], and β4 from brain [21]. The regulatory β subunits are mainly responsible for modulating the BK channel kinetic behavior, Ca2+ sensitivity and pharmacological responses to their specific modulators, such as scorpion toxins [22-25].

Generally, extracellular loops of most membrane proteins, including ion channels, are N-glycosylated. Oligosaccharides attached to the nascent polypeptide by recognizing an iconic amino acid sequon (Asn-X-Ser/Thr). Glycosylation has been demonstrated to be crucial to the activity, intracellular trafficking and targeting, and cell surface expression of ion channels [26-28]. Inhibition of N-linked glycosylation altered the gating properties of Kv1.1[29], Kv1.2 [30, 31] and KvLQT/minK (IsK) [32]. N-linked glycosylation also increases the expression of Shaker potassium channels [33] and HERG channels on the cell surface [34]. Most significantly, some genetic disorders, represented by long QT (LQT) syndrome was found to be induced by the mutations of glycosylation sites in HERG channels [35]. A point mutation of glycosylation site (N604T) changed the inactivation and ion selectivity rather than plasma membrane expression of TRPV1 [36]. In addition, glycosylation could alter the conductance properties of sodium channels [37], the gating properties of rNav1.3 [38] and steady-state inactivation of Nav1.9 endogenously expressed in dorsal root ganglion neurons [39]. Such phenomenon also presents in auxiliary subunits of ion channels. Glycosylation of the sodium channel β4 subunit is developmentally regulated and participates in neuritic degeneration [40]. The inhibitory effect of iberiotoxin on BK channels was enhanced by deglycosylation on β4 subunits [41]. Deglycosylation of the β1-subunit changes the open probability (Po) of BK channels in murine colonic smooth muscle cells [42]. It was also reported that there are 3 N-glycosylated sites (N-X-T/S at N88, N96 and N119) located at the β2 subunits of BK channels. Deglycosylation of these sites could regulate gating kinetics, outward rectification and toxin sensitivity of BK (α+β2) channels [43].

Given the clues mentioned above, glycosylation, either through α or β subunit, may probably participate in regulating the overall performances of channels in response to the dynamic and complicated changes in physiological conditions. Yet, the subtle but significant variations in channel conformation and extracellular environment may bring out uncertainties to elucidate the role of glycosylation in BK channel function. On the other hand, there remain questions on how specific ligands of BK channels, such as iberotoxin and martentoxin behave on deglycosylated BK channels, which needs to be considered in realizing the molecular interactions between toxins and channels. It also needs to be clarified how oligosaccharide chains on β1 influence the kinetic properties of BK channels. To these ends, this study will provide a possible explanation on these questions.

Methods

Construction of channel mutants

The plasmids containing hSloα (U23767) and β1 (KCNMB1; U25138) were gifts from N.W. Davies (University of Leicester) and J.D. Lippiat (Leeds university) [44]. We used the KOD mutagenesis kit (Toyobo, Japan) to make point mutations in the β1 subunit. In brief, PCR reactions were performed by using the wild-type β1 as the template and a pair of complementary mutagenesis primers. The primers were designed with Primer 5.0 (PremierBiosoft, USA). The PCR products were used to transform competent bacterial cells (DH5α, Tiangen, China) to amplify the mutant plasmid of hSloα+β1. All mutant constructs were verified by sequencing (Lifetechnologies, China).

Cell culture and transfection

All experiments were performed on HEK 293T cell lines. HEK 293T cells were obtained from Shanghai cell bank of Chinese Academy of Science. The cells were both cultured in Dulbecco’s modified Eagle medium (DMEM; Life Technologies, Grand Island, NY) supplemented with 10% heat-inactivated fetal bovine serum (FBS; Gibco, Grand Island, NY). Culture dishes were incubated at 37 °C in a humidified atmosphere containing 5% CO2, and subcultured approximately every 2-3 days. One day before transfection, HEK 293T cells were transferred to 24 well plates. At 90% confluence, cells were transiently transfected using Lipofectamine-2000 (Invitrogen, USA) at a ratio of 2 µL reagent with 1 µg total plasmid per well. Electrophysiological experiments and western blot analysis were performed at 1-2 days after transfection.

Western blot analysis

Cells were harvested 48 hours after transfection, washed twice with ice-cold phosphate-buffered saline (PBS) and lysed with cell lysis buffer on ice for 30 min. The cell lysates were centrifuged for 15 min at 14,000×g at 4 °C and the protein concentration of the supernatant was determined. Equivalent amounts of total protein (40 µg/lane) were separated by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and then transferred to poly (vinylidene fluoride) membranes (milipore,USA) using the submerged method. The membranes were incubated with rabbit anti-human IgG at room temperature for 2 h following the incubation of the Maxi-K beta antibody (diluted 1:200, Alomone Labs Ltd., Israel) overnight then washed with phosphate buffered saline plus tween-20 (PBST, tween-20 0.05%), and incubated for 1 h at room temperature with a horseradish peroxidase-conjugated secondary antibody (goat anti-rabbit IgG, 1:10000, Abcam, UK). Blots were visualized with Chemi DocTM XRS+ Imaging System (Bio-RAD, U.S.A).

Electrophysiological recordings

Whole-cell voltage-clamp experiments were performed following the procedures described previously [45], using an EPC-9 amplifier (HEKA Eletronik, Germany) at room temperature (21-25 °C). Patch pipettes were fabricated from glass capillary tubes by PC-10 Puller (Narishige, Japan) with the resistance of 2-3 MΩ. Data acquisition and stimulation protocols were controlled by a Pentium III computer (Legend, Beijing, China) equipped with Pulse/PulseFit 8.3 software (HEKA Eletronik, Germany). Capacitance transients were cancelled. Cells with a seal resistance (Rseal) below 1 GΩ were omitted. Series resistance (Rs) was compensated (80-85%) to minimize voltage errors, and cells with an uncompensated series resistance (Rs) above 10 MΩ were omitted. Leak subtraction was performed using P/6 protocol. Data were low-passed at 10 kHz. The rate of solution exchange was studied using solutions with different KCl concentrations and found to be about 95% complete within 20 s. Unless stated specially, For HEK 293T cells, the holding potential was -80 mV. All the recordings were performed with the pulse of +100 mV.

Solutions and drugs

In the patch-clamp recordings, the standard bath solution for HEK 293T cells was consisted of the following components (in mM): NaCl 135, KCl 5, MgCl2 1.2, CdCl2 2.5, HEPES 5, glucose 10 (pH 7.4 with NaOH). Pipette solutions for HEK 293T cells were composed of the following components (in mM): NaCl 10, KCl 117, MgSO4 2, HEPES 10, MgATP 2, EGTA 1 (pH 7.2 with KOH). The total Ca2+ to be added to give the desired free concentration was calculated using the program Maxchelator (http://www.stanford.edu/%7Ecpatton/maxc.html).

In the western blot analysis, the phosphate-buffered saline (PBS) was consisted of the following components (in mM): NaCl 135, KCl 4.7, Na2HPO4 10, NaH2PO4 2 (pH = 7.4).

BmK venom collected by electric stimulation was purchased from a scorpion farm in Zhengzhou, Henan province, China. Martentoxin was purified according to the method described by Ji [46].

The toxins were dissolved in the bath solution, supplemented with 1 mg/ml bovine serum albumin (BSA) in order to prevent adherence of the toxin to the vials and the perfusion apparatus. Application of 1 mg/ml BSA alone did not alter BK channel function. Unless otherwise stated, all reagents were purchased from Sigma.

Data analysis

Data were analyzed by PulseFit 8.5 (HEKA Eletronik, Germany) and Origin 7.5 (Northampton, Massachusetts, USA). Results of data analysis were expressed as mean ± S.E.M. and n represents the number of the cells examined.

The Statistical significance was determined using the unpaired Student’s t-Test or one-way ANOVA, and an asterisk denotes P<0.05 unless otherwise stated. The degree of toxin effect was calculated by expressing the remaining current after each drug exposure as a fraction of the current magnitude of the patch prior to the first drug exposure (i.e., fractional current remaining, If).

BK channel (α+β1) currents were elicited by the step pulses ranging from -50 to +120 mV for 200 ms with the increments of 10 mV. The holding potentials were held at -80 mV for BK channel (α+β1). Current density calculation formula (pA/pF), where pA represents the current of BK channel (α+β1) and pF represents the membrane area of measured cell. For determining the voltage dependence of activation, the conductance was calculated using the formula: G(V)= I(V)/(V-ErK), where I(V) is the current of BK channel (α+β1) at the command voltage V, and ErK is the reversal potential. The conductance was normalized to the maximal value and the voltage dependence for activation of BK channel (α+β1) fitted to a Boltzmann equation: f(x) = -1/(1+exp((x-V1/2)/k))+1, where V1/2 is the voltage at which half-maximal activation occurs, and k describes the slope of the fit.

The activation kinetics was analyzed by fitting the decay course of the BK channel currents to a single exponential function: f(t)= C+A×exp[-(t-t0)/τ], where t is time, t0 is the time when the currents were just starting to exponentially increase, A represents the amplitudes of the channels activating with the time constant τ and C is the steady-state asymptote approximating to the maximum value of the current amplitude (Imax).

Results

Glycosylation of β1 co-expressed with hSloα subunit in HEK293T cells

We designed 4 pairs of complementary mutagenesis primers to construct the N80A/N142A mutant and N80Q/N142Q mutant of the β1 subunit of the BK channel, which were identified by sequencing (Additional file 1, 1A-1C). To confirm whether β1-subunit expressed in the HEK 293T cells was glycosylated, the western blotting was used to detect the molecular weight of β1 mutations (N80A/N142A and N80Q/N142Q, Figure 1C) before the whole experiment. Both lysates from wild type (β1)-expressed cells and double mutations (β1 N80A/N142A and β1 N80Q/N142Q)-expressed cells were assayed with antibody specific for BK β1-subunit. A ~25 kDa band observed here represents the wild-type β1-subunit (Figure 1C, first lane), which is consistent with the previous report [42]. The double mutants of glycosylation sites (N80A/N142A) were detected with a smaller molecular mass of ~20 kDa (Figure 1C, second lane). Similar result was also observed from the analysis for another double mutant (N80Q/N142Q) (Figure 1C, third lane). Different mutations (N80A/N142A and N80Q/N142Q) were constructed to rule out the significant alteration of β1 molecular weight and confirm the removal of oligosaccharide from the β1 glycosylation sites. Furthermore, tunicamycin, a potent inhibitor of N-linked glycosylation, was employed [38, 41] to make sure that deglycosylation on β1-subunit was achieved successfully by site-directed mutation. 30 min after HEK293T cells were transfected with the hSloα+β1 (wild-type), tunicamycin at the final concentration of 10 μM was added and the cells were grown under the presence of tunicamycin for 15 h. Cell lysates were assayed by western blot. Wild-type β1-subunit from tunicamycin-treated cells also has a same molecular weight of ~20 kDa (Figure 1C, forth lane) as the double mutant (β1 N80A/N142A or β1 N80Q/N142Q). This indicated that β1 is glycosylated in HEK 293T cells. N80 and N142 were the correct N-linked glycosylation sites in β1-subunit. Western blotting and electrophysiological assay showed that there were no significant difference in protein relative expression level and current density between the wild-type β1 subunit, the two double mutants and the tunicamycin treated group (Figure 1D and Additional file 1D). In other words, β1 subunit mutants or the use of glycation inhibitors did not lead to changes in expression levels.

Figure 1.
Molecular properties of BK channel β1-subunit. (A) The primary structure of a subunit protein, including N-glycosylation site (white boxes) and transmembrane segment (gray boxes). (B) Model of BK (α+β1) channels (two of the four subunits are shown), including segments of α subunits, extracellular loop of β1 and oligosaccharide chains on β1 (dotted line). (C) Western blot showing β1-subunit. Molecular mass of the ~25 kDa immunoreactive band was wild-type β1-subunit (lane 1). Immunoreactive protein, with a molecular mass of ~20 kDa was double glycosylation site mutant β1 N80A/N142A (lane 2) and β1 N80Q/N142Q (lane 3). The immunoreactive band with a molecular mass of ~20 kDa on lane 4 was the β1-subunit by treatment with tunicamycin. There was no significant difference among the wild-type β1 subunit, N80A/N142A mutant, N80Q/N142Q mutant and the tunicamycin treated group in protein relative expression level (p > 0.05).

Deglycosylation of β1 changing the toxin sensitivity of BK channel (α+β1)

To assess the pharmacological sensitivity of deglycosylated BK channels, wild type and mutant BK (α+β1) channels were expressed in HEK293T cells. The cells transfected with wild type and mutant BK (α+β1) channels were pre-treated with 200 nM thapsigargin, an irreversible inhibitor of the sarcoplasmic reticulum Ca2+-ATPase pump for 30 min at 37 °C in order to deplete intracellular Ca2+ stores [47] before performing whole cell patch recordings. When the free Ca2+ concentration in the pipette solution was buffered to 300 nM, the currents of BK (α+β1) channel evoked by +100 mV pulse could not be effectively blocked by iberiotoxin, a specific and potent blocker of BK channels to hSloα [21], even at the concentration of 100 nM (If: 0.92 ± 0.01, n = 7, Figure 2A). However, the double mutant (β1 N80A/N142A) does not protect the channel against iberiotoxin at this concentration (If: 0.46 ± 0.04, n = 7, Figure 2B, p < 0.001). In addition, another double glycosylation site mutant β1 (β1 N80Q/N142Q) also showed high sensitivity to iberiotoxin (If: 0.46 ± 0.04, n = 6, Figure 2C, p < 0.001). As summarized in Figure 2D, both double mutants N80A/N142A and N80Q/N142Q had significant difference of the sensitivity to the inhibition of 100 nM iberiotoxin from wild-type BK (α+β1) channels. It should be noted here that two double mutants don’t have any obvious difference of the sensitivity to iberiotoxin.

Figure 2.
The inhibition effect of iberiotoxin on BK channels (a+β1) could be enhanced after deglycosylation of β1 subunit. (A) Representative whole cell current traces from HEK 293T cells expressing wild-type BK channels (a+β1) before and after the application of IbTx 100 nM. The holding voltage was -80 mV and the currents were elicited by a pulse of +100 mV in the presence of 300 nM free Ca2+ in the pipette solution. (B) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80A/N142A before and after the application of IbTx 100 nM. (C) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80Q/N142Q before and after the application of IbTx 100 nM. (D) Statistical analysis of pharmacological modulation of wild-type BK channels (a+β1) (n = 7) and N80A/N142A (n = 7, ***p < 0.001) or wild-type BK channels (a+β1) (n = 7) and N80Q/N142Q (n = 6, ###p < 0.001) by IbTx 100 nM. The inhibition ratios of IbTx on the currents of N80A/N142A and N80Q/N142Q were not significantly different (p > 0.05).

To investigate this pharmacological modulation of BK channels by deglycosylation further, martentoxin, a specific modulator of BK channel subtypes was also used here. Martentoxin can significantly increase the activity of the BK channel (α+β1) subtype with an EC50 of 495 nM [48]. Under the same condition as above, the evoked currents of BK (α+β1) channels were increased by martentoxin at the concentration of 400 nM (If: 1.13 ± 0.01, n = 6, Figure 3A). In contrast to this enhancement of glycosylated channel activity by martentoxin, both the double mutants N80A/N142A and N80Q/N142Q were insensitive to 400 nM martentoxin (If: 0.99 ± 0.02, n = 5, p < 0.001, Figure 3B; If: 1.04 ± 0.02,n = 6, p < 0.001, Figure 3C), And there were no significant differences between N80A/N142A and N80Q/N142Q of their sensitivity to martentoxin (Figure 3D).

Figure 3.
Deglycosylation of β1 subunit removed the enhancement effect of martentoxin on BK channels (α+β1). (A) Representative whole cell current traces from HEK 293T cells expressing wild-type BK channels (α+β1) before and after the application of martentoxin 400 nM. The holding voltage was -80 mV and the currents were elicited by a pulse of +100 mV in the presence of 300 nM free Ca2+ in the pipette solution. (B) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80A/N142A before and after the application of martentoxin 400 nM. (C) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80Q/N142Q before and after the application of martentoxin 400 nM. (D) Statistical analysis of pharmacological modulation of wild-type BK channels (α+β1) (n = 6) and N80A/N142A (n = 5, ***p < 0.001) or wild-type BK channels (α+β1) (n = 6) and N80Q/N142Q (n = 6, ###p < 0.001) by martentoxin 400 nM. The effect ratios of martentoxin on the currents of N80A/N142A and N80Q/N142Q were not significantly different (p > 0.05).

When the pipette solution contained 10 mM EGTA without any Ca2+, the channels could only be activated by voltage. The currents of glycosylated BK (α+β1) channels could not be modulated by martentoxin at 400 nM (If: 0.97 ± 0.018, n = 3, Figure 4A). In contrast, the currents of double mutants N80A/N142A and N80Q/N142Q could be potently inhibited by 400 nM martentoxin (N80A/N142A If: 0.71 ± 0.064, n = 3, p < 0.001; N80Q/N142Q If: 0.72 ± 0.099, n = 4, p < 0.001, Figures 4B and 4C). The inhibition ratios of martentoxin on the currents of N80A/N142A and N80Q/N142Q were not significantly different (Figure 4D).

Figure 4.
Deglycosylation of β1 subunit altered the pharmacological effect of martentoxin on BK channels (α+β1) in the absence of cytoplasmic Ca2+. (A) Representative whole cell current traces from HEK 293T cells expressing wild-type BK channels (α+β1) before and after the application of martentoxin 400 nM. The holding voltage was -80 mV and the currents were elicited by a pulse of +100 mV. (B) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80A/N142A before and after the application of martentoxin 400 nM. (C) Representative whole cell current traces from HEK 293T cells expressing mutant BK channels N80Q/N142Q before and after the application of martentoxin 400 nM. (D) Statistical analysis of pharmacological modulation of wild-type BK channels (α+β1) (n = 3) and N80A/N142A (n = 3, **p < 0.01) or wild-type BK channels (α+β1) (n = 3) and N80Q/N142Q (n = 4, ##p < 0.01) by martentoxin 400 nM. The inhibition ratios of martentoxin on the currents of N80A/N142A and N80Q/N142Q were not significantly different (p > 0.05).

Deglycosylation of β1 modifying the activation kinetics of BK channel (α+β1) in HEK293T cells

It has been well-established that the gating properties of BK channel depend largely on the interactions between α and β subunits [19, 21, 49-51]. Therefore, the subtle alterations in structure caused by glycosylation would have potential impact on the gating kinetics of BK channel (α+β1). To test this hypothesis, the voltage-dependent activation of the deglycosylated β1 mutants (β1 N80A/N142A and N80Q/N142Q) was examined. HEK293T cells transfected with wild-type or mutant BK (α+β1) channels were also pre-treated with 200 nM thapsigargin for 30 min at 37 °C. The free Ca2+ concentration in the pipette solution was also controlled to 300 nM in the patch clamp recording. The BK (α+β1) currents were elicited by the step pulses ranging from -50 to +120 mV for 200 ms with the increments of 10 mV. The activation curve of the double glycosylation site β1 mutant (β1 N80Q/N142Q) as well as the half-maximal voltage (V1/2) of activation was not significantly shifted compared to the wild-type BK (α+β1) channels (WT: V1/2 = 66.52 ± 1.47,n = 9, p > 0.05; β1 N80Q/N142Q: V1/2 = 65.35 ± 1.48, n = 10, p > 0.05, Figure 5B and Table 1).

Figure 5.
Deglycosylation of β1 subunit did not shift steady-state activation of BK channels (α+β1). (A) Representative whole cell current traces from HEK 293T cells expressing wild type BK channels (α+β1) (left), N80Q/N142Q (right). The channels were stepped from -50 mV to +120 mV in 10 mV steps with the holding potential at -80 mV. (B) The plots of the normalized conductance fit well with Boltzmann function. The voltage-dependent activation curves were insignificantly shifted before and after the deglycosylation. (C) Time constants are shown for wild-type BK channels (α+β1) (n = 5, black circle), N80Q/N142Q (n = 7, gray circle).

Table 1.
Effect of deglycosylation on V1/2 and slope of conductance vs. voltage curves.

To evaluate the influence of deglycosylation on activation rate, we compared the activation time constants between wild-type and mutant channels. Time constants of activation at each potential were determined by fitting each current with a single exponential equation (in Materials and Methods, Data analysis). When the whole cell current tracings were evoked at each potential from +50 mV to +120 mV for wild-type BK (α+β1) channels and mutant channels, it was observed that the activation time constants was remarkably decreased in β1 N80Q/N142Q mutants compared to that of wild-type channels (n = 7, p<0.001 from +50 to +70 mV, p<0.01 from +80 mV to +90 mV, p<0.05 from +100 mV to +120 mV, Figure 5C).

The inhibitory effects of IbTX and MarTX on BK channels (α+β1) were detected by exerting glycosidase

In order to exclude the structural changes induced by a single or double amino acid mutation in β1 subunits, which affect the activity of channels, wild type of BK (α+β1) channels were expressed in HEK293T cell. The cells transfected with wild type BK (α+β1) channels were treated with PNGase F (1000 units) for 1 hour at 37 °C so that PNGase F could remove the N-glycosylation of β1 subunits before performing whole cell patch recordings. When the free Ca2+ concentration in the pipette solution was buffered to 300 nM, the currents of BK (α+β1) channel, pre-treated with PNGase F, evoked by +100 mV pulse could be effectively blocked by iberiotoxin at the concentration of 100 nM and 400 nM (Figure 6A and 6B). Statistical analysis showed that pre-treated with PNGase F, the inhibition of both concentrations of IbTX on BK (α+β1) channel had been significantly enhanced (100 nM IbTX CTRL group: If = 0.89 ± 0.01vs. 100 nM IbTX PNGase F group: If = 0.54 ± 0.05; 400 nM IbTX CTRL group: If = 0.65 ± 0.03 vs. 400 nM IbTX PNGase F group: If = 0.43 ± 0.06, Figure 6C). When the pipette solution without any Ca2+, the BK channel (α+β1) is activated only by voltage. The currents of BK (α+β1) channels could be obviously inhibited by MarTX (Figure 7). Statistical analysis showed that inhibitory ratio of the PNGase F group was more pronounced than that of control group (100 nM MarTX CTRL group: If = 1.00 ± 0.03vs. 100 nM MarTX PNGase F group: If = 0.83 ± 0.03; 400 nM MarTX CTRL group: If = 1.00 ± 0.03 vs. 400 nM MarTX PNGase F group: If = 0.63 ± 0.02, Figure 7C).

Figure 6.
The application of PNGase F enhances the inhibitory effects of IbTX on BK (α+β1) channel. (A-B) Whole-cell current traces from deglycosylation BK channels before and after the application of IbTx (A) 100 nM and (B) 400 nM without treatment (left, CTRL) or pretreatment with PNGase F (right, PNGase F). The holding voltage was -80 mV and the currents were elicited by a pulse of +100 mV. (C) Statistical analysis of the inhibitory ratios of IbTX on BK (α+β1) channel without treatment (left, CTRL) or pretreatment with PNGase F.

Figure 7.
The application of PNGase F enhances the inhibitory effects of MarTX on BK (α+β1) channel. (A-B) Whole-cell current traces from deglycosylation BK channels before and after the application of MarTX (A) 100 nM and (B) 400 nM without treatment (left, CTRL) or pretreatment with PNGase F (right, PNGase F). The holding voltage was -80 mV and the currents were elicited by a pulse of +100 mV. (C) Statistical analysis of the inhibitory ratios of MarTX on BK (α+β1) channels without treatment (left, CTRL) or pretreatment with PNGase F.

Discussion

BK channels (α+β1) are mainly distributed in smooth muscle, hair cells and some of the neurons, which are often considered as the direct target of many natural products. Researches on the role of N-linked glycosylation of β subunit may contribute to the understanding of the molecular mechanisms underlying the interaction of BK channels with specific ligands or auxiliary subunits as well as further reveal physiological functions of BK channels in cells and organs.

Regulation on the pharmacological characterization of BK channel subtypes by deglycosylation of β-subunit

β-subunit usually changes the pharmacological property of BK channels, which is different from that of other potassium channels. β1-subunit co-expressed with the BK channel could increase the sensitivity to charybdotoxin, one of blockers for BK channel subtypes, whereas the BK channels formed by α and β4 subunits (BK (α+β4) channels) become insensitive to charybdotoxin [21, 23]. Another BK channel blocker, iberiotoxin, could potently inhibit the currents of BK channels which consist of α subunit alone. Co-expressed with either β1 or β4 subunits, BK channel subtypes become insensitive to iberiotoxin, respectively [21]. Similar phenomenon could also be seen for the BK channels modulated by slotoxin, a polypeptide purified from the venom of Centruroides noxius [52].

As an extracellular part of β subunits, it is worthy of our consideration whether glycosylation is important to adjust the toxin sensitivity of BK channel subtypes, especially participating in influencing the interaction between channels and toxins. In the present study, BK channels co-expressed with deglycosylated β1 mutants (N80Q/N142Q) were found to be more sensitive to iberiotoxin than that of wild type β1. Actually, such phenomenon is also found in deglycosylated β4 subunit. Deglycosylated β4 becomes less effective than wild type in protecting the channel against the block by iberiotoxin [41]. The results illuminated that oligosaccharide chains on the extracellular loop of β subunit might act as a shield to prevent toxins binding to the BK channel pore. However, some evidences supported another possibility that the removal of the oligosaccharides caused a conformational change in the interaction between α and β subunit. When BK (α+β4) channels were pretreated with martentoxin, the subsequent iberiotoxin acted more potently to inhibit the BK currents than applied it alone [53]. When three basic residues (Lys-120, Arg-121, and Lys-125) of β4 subunit were mutated, the affinity of charybdotoxin on BK channels could be significantly enhanced [54]. Oligosaccharides were not removed from BK channels in the reports above. It may allow us to infer that the deglycosylation on β subunits could alter the conformation of BK channels, the same effects as the pretreated toxin or site mutation. The conformational change increased binding affinities between BK channels and toxins.

As an effective tool to discriminate different BK channel subtypes, martentoxin was tailored to study the regulatory effects of glycosylation on BK channels. The current amplitude of wild-type BK (α+β1) channels could be increased by martentoxin. Interestingly, when co-expressed with non-glycosylated β1 (N80Q/N142Q or N80A/N412A), channels were found to be insensitive to martentoxin. In contrast, martentoxin has no significant effect on BK (α+β1) channels when intracellular free calcium was deprived. When two glycosylation sites (N80Q/N142Q or N80A/N412A) were mutated in β1 subunit, BK (α+β1) channels could be inhibited by martentoxin. The results suggest that the inhibitory sites of the BK (α+β1) channels for martentoxin were prevented by oligosaccharides.

Modification on the kinetic properties of BK channels by deglycosylation of β1-subunit

β1 subunit increases the Ca2+ sensitivity of BK channels by altering the apparent affinity of the channels’ Ca2+ binding sites and slows down the activation kinetics of these channels. BK channels containing the α-subunits co-expressed with β1-subunit produce sustained currents that do not inactivate [49, 55].

Oligosaccharide chains, the ‘branched chain’ of β1-subunit, may regulate the kinetic properties of the BK channels. In the present study, the double glycosylation site mutants (β1 N80Q/N142Q) could remarkably accelerate the activation of BK channel compared to that of wild-type channels. In contrast, the activation time constants of single mutation for either β1 N80Q or N142Q mutant had no significant difference with wild-type channels. It could be infered that these two glycosylation sites on β1 subunit were synergically coupled in β1. The alteration of activation kinetics resulted from completely deglycosylation on β1 subunit. In another report, the open probability (Po) and mean open time of BK (α+β1) channels were both increased by deglycosylation of β1 [42]. The results demonstrated that deglycosylation on β1 could change the activation of BK channel, which might further influence the channel function directly.

It is worth mentioning that the regulation of glycosylation on channel kinetics might be influenced by other environmental factors. In a previous report, deglycosylation of β1 could shift the voltage of steady-state activation toward more positive potentials, when Ca2+ concentration was controlled lower than 1 μM [42], which was different from our data (Figure 5B). In fact, the results are not contradictory possibly because of the following reasons: 1) Depend on calcium sensitivity of BK channels. The calcium sensitivity of BK channels formed by hsloα+β1 might be higher than the BK channels endogenously expressed in murine smooth muscle cells. The hsloα+β1 channels at low calcium concentration (300 nM) show the similar effects of murine BK channels at high calcium concentration (>1 μM) [42]. 2) Depend on membrane environment of BK channels. Unlike previous report, the current of BK (α+β1) channels was recorded by different patch clamping mode in present study. Consequently, the resulted midpoint voltage (V1/2) was somewhat different compared to the previous one, which could be explained by the subtle changes in mechanical state of the lipid membrane, which was closely related to gating properties [56].

Physiological and pharmacological significance of glycosylation on β1 subunit of BK channels

BK (α+β1) channels activated by local Ca2+ release could regulate the membrane potential of arterial smooth muscle cells. Loss of the β1 subunit leads to hypertension, whereas a gain in β1 subunit function is associated with protection against hypertension in humans [57]. In present study, glycosylation on β1 subunits effectively protect BK channels against the regulation of modulators, which possibly implied an important role of glycosylation in maintaining blood pressure and the tone of vascular under physiological condition. On the other hand, deglycosylation on β1 subunits might have great significance in pathological conditions due to accelerating the activation as well as increasing open probability of BK channels, which makes the cell membrane tend to hyperpolarization and smooth muscles more relaxed. Deglycosylation will deserve to be considered as a target in the treatment of channelopathies and drug design. Fusion with glycosylase, specific peptides targeting BK channel might be selected as a therapeutic agent to mitigate the symptom of hypertension. Anyway, the realization of these ideas and resolution of mechanisms still need a large number of experiments in the future.

In summary, deglycosylation on β1 subunits increased the sensitivity to iberotoxin and martentoxin as well as accelerate the activation of BK (α+β1) channels. It would have important implications on the diversity of structure and function of BK channels. Especially at molecular level, in addition to key residues or important region on channel proteins, glycosylation become another indispensable factor for regulating the pharmacological properties of BK channel subtypes.

Conclusion

The present study investigated that deglycosylation of β1 subunits through double-site mutations (β1 N80A/N142A or β1 N80Q/N142Q) could significantly increase the inhibitory potency of iberiotoxin. Meanwhile, the deglycosylated channels also have a different sensitivity to martentoxin. On the contrary to the enhancing effects of martentoxin on glycosylated BK channels under the presence of cytoplasmic Ca2+, deglycosylated channels were not affected by the toxin. Surprisingly, the deglycosylated channels were inhibited by martentoxin under the absence of cytoplasmic Ca2+, while the glycosylated channels were not inhibited under this same condition. In addition, wild type BK (α+β1) channels were treated with PNGase F also showed the same trend of pharmacological results to the mutants. Similar to this modulation of glycosylation on BK channel pharmacology, the deglycosylated forms of the channels were activated at a faster speed than the glycosylated ones though the V1/2 and slope were not changed by the glycosylation. In a word, the present work reveals that glycosylation is an indispensable determinant of the modulation of β1-subunit on BK channel pharmacology and its activation. The loss of glycosylation of β1 subunits could lead to the dysfunction of BK channel, resulting in a pathological state.

Abbreviations

BK channels: voltage-dependent large-conductance Ca2+-activated K+ channels; Kv channel: voltage-gated K+ channel; LQT: long QT syndrome; Rseal: seal resistance; Rs: series resistance; If: fractional current remaining; Imax: the maximum value of the current amplitude; V1/2: the half-maximal voltage; Po: the open probability.

Acknowledgments

We are grateful to Dr. J. D. Lippiat and N. W. Davies (University of Leicester, UK) for providing the plasmid comprising genes of hSloα and β1.

References

  • 1. Barnes EA, Lee L, Barnes SL, Brenner R, Alvira CM, Cornfield DN. β1-Subunit of the calcium-sensitive potassium channel modulates the pulmonary vascular smooth muscle cell response to hypoxia. Am J Physiol Lung Cell Mol Physiol. 2018 Aug;315(2):L265-75.
  • 2. Semenov I, Wang B, Herlihy JT, Brenner R. BK channel beta1-subunit regulation of calcium handling and constriction in tracheal smooth muscle. Am J Physiol Lung Cell Mol Physiol. 2006 Oct;291(4):L802-10.
  • 3. Parajuli SP, Zheng YM, Levin R, Wang YX. Big-conductance Ca(2+)-activated K(+) channels in physiological and pathophysiological urinary bladder smooth muscle cells. Channels (Austin). 2016 Sep 2;10(5):355-64.
  • 4. Du W, Bautista JF, Yang H, Diez-Sampedro A, You SA, Wang L, et al. Calcium-sensitive potassium channelopathy in human epilepsy and paroxysmal movement disorder. Nat Genet. 2005 Jul;37(7):733-8.
  • 5. Castillo K, Contreras GF, Pupo A, Torres YP, Neely A, González C, et al. Molecular mechanism underlying β1 regulation in voltage- and calcium-activated potassium (BK) channels. Proc Natl Acad Sci U S A. 2015 Apr 14;112(15):4809-14.
  • 6. Marty A. The physiological role of calcium-dependent channels. Trends Neurosci. 1989 Nov;12(11):420-4.
  • 7. Duncan PJ, Shipston MJ. BK channels and the control of the pituitary. Int Rev Neurobiol. 2016 Mar;128:343-68.
  • 8. Maqoud F, Curci A, Scala R, Pannunzio A, Campanella F, Coluccia M, et al. Cell cycle regulation by Ca(2+)-activated K⁺ (BK) channels modulators in SH-SY5Y neuroblastoma cells. Int J Mol Sci. 2018 Aug 18;19(8):2442.
  • 9. Weaver AK, Liu X, Sontheimer H. Role for calcium-activated potassium channels (BK) in growth control of human malignant glioma cells. J Neurosci Res. 2004 Oct 15;78(2):224-34.
  • 10. Cheng YY, Wright CM, Kirschner MB, Williams M, Sarun KH, Sytnyk V, et al. KCa1.1, a calcium-activated potassium channel subunit alpha 1, is targeted by miR-17-5p and modulates cell migration in malignant pleural mesothelioma. Mol Cancer. 2016 Jun 1;15(1):44.
  • 11. Kraft R, Krause P, Jung S, Basrai D, Liebmann L, Bolz J, et al. BK channel openers inhibit migration of human glioma cells. Pflugers Arch. 2003 May;446(2):248-55.
  • 12. Atkinson NS, Robertson GA, Ganetzky B. A component of calcium-activated potassium channels encoded by the Drosophila slo locus. Science. 1991 Aug 2;253(5019):551-5.
  • 13. Wallner M, Meera P, Toro L. Determinant for beta-subunit regulation in high-conductance voltage-activated and Ca(2+)-sensitive K+ channels: an additional transmembrane region at the N terminus. Proc Natl Acad Sci U S A. 1996 Dec 10;93(25):14922-7.
  • 14. Almássy J, Nánási PP. Brief structural insight into the allosteric gating mechanism of BK (Slo1) channel (1). Can J Physiol Pharmacol. 2019 Jun;97(6):498-502.
  • 15. Schreiber M, Salkoff L. A novel calcium-sensing domain in the BK channel. Biophys J. 1997 Sep;73(3):1355-63.
  • 16. Knaus HG, Garcia-Calvo M, Kaczorowski GJ, Garcia ML. Subunit composition of the high conductance calcium-activated potassium channel from smooth muscle, a representative of the mSlo and slowpoke family of potassium channels. J Biol Chem. 1994 Feb 11;269(6):3921-4.
  • 17. Martinez-Espinosa PL, Yang C, Gonzalez-Perez V, Xia XM, Lingle CJ. Knockout of the BK β2 subunit abolishes inactivation of BK currents in mouse adrenal chromaffin cells and results in slow-wave burst activity. J Gen Physiol. 2014 Oct;144(4):275-95.
  • 18. Solaro CR, Prakriya M, Ding JP, Lingle CJ. Inactivating and noninactivating Ca(2+)- and voltage-dependent K+ current in rat adrenal chromaffin cells. J Neurosci. 1995 Sep;15(9):6110-23.
  • 19. Brenner R, Jegla TJ, Wickenden A, Liu Y, Aldrich RW. Cloning and functional characterization of novel large conductance calcium-activated potassium channel beta subunits, hKCNMB3 and hKCNMB4. J Biol Chem. 2000 Mar 3;275(9):6453-61.
  • 20. Poulsen AN, Wulf H, Hay-Schmidt A, Jansen-Olesen I, Olesen J, Klaerke DA. Differential expression of BK channel isoforms and beta-subunits in rat neuro-vascular tissues. Biochim Biophys Acta. 2009 Feb;1788(2):380-9.
  • 21. Meera P, Wallner M, Toro L. A neuronal beta subunit (KCNMB4) makes the large conductance, voltage- and Ca2+-activated K+ channel resistant to charybdotoxin and iberiotoxin. Proc Natl Acad Sci U S A. 2000 May 9;97(10):5562-7.
  • 22. Tao J, Shi J, Yan L, Chen Y, Duan YH, Ye P, et al. Enhancement effects of martentoxin on glioma BK channel and BK channel (alpha+beta1) subtypes. PLoS One. 2011 Mar 18;6(3):e15896.
  • 23. Hanner M, Vianna-Jorge R, Kamassah A, Schmalhofer WA, Knaus HG, Kaczorowski GJ, et al. The beta subunit of the high conductance calcium-activated potassium channel. Identification of residues involved in charybdotoxin binding. J Biol Chem. 1998 Jun 26;273(26):16289-96.
  • 24. Gonzalez-Perez V, Lingle CJ. Regulation of BK Channels by Beta and Gamma Subunits. Annu Rev Physiol. 2019 Feb 10;81:113-37.
  • 25. Torres YP, Granados ST, Latorre S. Pharmacological consequences of the coexpression of BK channel α and auxiliary β subunits. Front Physiol. 2014 Oct 10;5:383.
  • 26. Lazniewska J, Weiss N. Glycosylation of voltage-gated calcium channels in health and disease. Biochim Biophys Acta Biomembr. 2017 May;1859(5):662-8.
  • 27. Lazniewska J, Weiss N. The "sweet" side of ion channels. Rev Physiol Biochem Pharmacol. 2014;167:67-114.
  • 28. Wiedmann F, Schlund D, Faustino F, Kraft M, Ratte A, Thomas D, et al., N-Glycosylation of TREK-1/hK(2P)2.1 Two-Pore-Domain Potassium (K(2P)) Channels. Int J Mol Sci. 2019 Oct 20;20(20):5193.
  • 29. Sutachan JJ, Watanabe I, Zhu J, Gottschalk A, Recio-Pinto E, Thornhill WB. Effects of Kv1.1 channel glycosylation on C-type inactivation and simulated action potentials. Brain Res. 2005 Oct 5;1058(1-2):30-43.
  • 30. Thayer DA, Yang SB, Jan YN, Jan LY. N-linked glycosylation of Kv1.2 voltage-gated potassium channel facilitates cell surface expression and enhances the stability of internalized channels. J Physiol. 2016 Nov 15;594(22):6701-13.
  • 31. Watanabe I, Zhu J, Sutachan JJ, Gottschalk A, Recio-Pinto E, Thornhill WB. The glycosylation state of Kv1.2 potassium channels affects trafficking, gating, and simulated action potentials. Brain Res. 2007 May 4;1144:1-18.
  • 32. Freeman LC, Lippold JJ, Mitchel KE. Glycosylation influences gating and pH sensitivity of I(sK). J Membr Biol. 2000 Sep 1;177(1):65-79.
  • 33. Khanna R, Myers MP, Lainé M, Papazian DM. Glycosylation increases potassium channel stability and surface expression in mammalian cells. J Biol Chem. 2001 Sep 7;276(36):34028-34.
  • 34. Lamothe SM, Hulbert M, Guo J, Li W, Yang T, Zhang S. Glycosylation stabilizes hERG channels on the plasma membrane by decreasing proteolytic susceptibility. Faseb J. 2018 Apr;32(4):1933-43.
  • 35. Satler CA, Vesely MR, Duggal P, Ginsburg GS, Beggs AH. Multiple different missense mutations in the pore region of HERG in patients with long QT syndrome. Hum Genet. 1998 Mar;102(3):265-72.
  • 36. Veldhuis NA, Lew MJ, Abogadie FC, Poole DP, Jennings EA, Ivanusic JJ, et al. N-glycosylation determines ionic permeability and desensitization of the TRPV1 capsaicin receptor. J Biol Chem. 2012 Jun 22;287(26):21765-72.
  • 37. Cronin NB, O'Reilly A, Duclohier H, Wallace BA. Effects of deglycosylation of sodium channels on their structure and function. Biochemistry. 2005 Jan 18;44(2):441-9.
  • 38. Xu Q, Cheng HW, He HQ, Liu ZR, He M, Yang HT, et al. Deglycosylation altered the gating properties of rNav1.3: glycosylation/deglycosylation homeostasis probably complicates the functional regulation of voltage-gated sodium channel. Neurosci Bull. 2008 Oct;24(5):283-7.
  • 39. Tyrrell L, Renganathan M, Dib-Hajj SG, Waxman SG. Glycosylation alters steady-state inactivation of sodium channel Nav1.9/NaN in dorsal root ganglion neurons and is developmentally regulated. J Neurosci. 2001 Dec 15;21(24):9629-37.
  • 40. Zhou TT, Zhang ZW, Liu J, Zhang JP, Jiao BH. Glycosylation of the sodium channel beta4 subunit is developmentally regulated and involves in neuritic degeneration. Int J Biol Sci. 2012;8(5):630-9.
  • 41. Jin P, Weiger TM, Levitan IB. Reciprocal modulation between the alpha and beta 4 subunits of hSlo calcium-dependent potassium channels. J Biol Chem. 2002 Nov 15;277(46):43724-9.
  • 42. Hagen BM, Sanders KM. Deglycosylation of the beta1-subunit of the BK channel changes its biophysical properties. Am J Physiol Cell Physiol. 2006 Oct;291(4):C750-6.
  • 43. Huang ZG, Liu HW, Yan ZZ, Wang S, Wang LY, Ding JP. The glycosylation of the extracellular loop of β2 subunits diversifies functional phenotypes of BK Channels. Channels (Austin). 2017 Mar 4;11(2):156-66.
  • 44. Lippiat JD, Standen NB, Harrow ID, Phillips SC, Davies NW. Properties of BK(Ca) channels formed by bicistronic expression of hSloalpha and beta1-4 subunits in HEK293 cells. J Membr Biol. 2003 Mar 15;192(2):141-8.
  • 45. Hamill OP, Marty A, Neher E, Sakmann B, Sigworth FJ. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflugers Arch. 1981 Aug;391(2):85-100.
  • 46. Ji YH, Wang WX, Ye JG, He LL, Li YJ, Yan YP, et al. Martentoxin, a novel K+-channel-blocking peptide: purification, cDNA and genomic cloning, and electrophysiological and pharmacological characterization. J Neurochem. 2003 Jan;84(2):325-35.
  • 47. Weaver AK, Olsen ML, McFerrin MB, Sontheimer H. BK channels are linked to inositol 1,4,5-triphosphate receptors via lipid rafts: a novel mechanism for coupling [Ca(2+)](i) to ion channel activation. J Biol Chem. 2007 Oct 26;282(43):31558-68.
  • 48. Tao J, Shi J, Yan L, Chen Y, Duan YH, Ye P, et al. Enhancement effects of martentoxin on glioma BK channel and BK channel (α+β1) subtypes. PLoS One. 2011 Mar 18;6(3):e15896.
  • 49. Orio P, Rojas P, Ferreira G, Latorre R. New disguises for an old channel: MaxiK channel beta-subunits. News Physiol Sci. 2002 Aug;17:156-61.
  • 50. Xia XM, Ding JP, Lingle CJ. Molecular basis for the inactivation of Ca2+- and voltage-dependent BK channels in adrenal chromaffin cells and rat insulinoma tumor cells. J Neurosci. 1999 Jul 1;19(13):5255-64.
  • 51. Wallner M, Meera P, Toro L. Molecular basis of fast inactivation in voltage and Ca2+-activated K+ channels: a transmembrane beta-subunit homolog. Proc Natl Acad Sci U S A. 1999 Mar 30;96(7):4137-42.
  • 52. Garcia-Valdes J, Zamudio FZ, Toro L, Possani LD. Slotoxin, alphaKTx1.11, a new scorpion peptide blocker of MaxiK channels that differentiates between alpha and alpha+beta (beta1 or beta4) complexes. FEBS Lett. 2001 Sep 21;505(3):369-73.
  • 53. Shi J, He HQ, Zhao R, Duan YH, Chen J, Chen Y, et al. Inhibition of martentoxin on neuronal BK channel subtype (alpha+beta4): implications for a novel interaction model. Biophys J. 2008 May 1;94(9):3706-13.
  • 54. Gan G, Yi H, Chen M, Sun L, Li W, Wu Y, et al. Structural basis for toxin resistance of beta4-associated calcium-activated potassium (BK) channels. J Biol Chem. 2008 Aug 29;283(35):24177-84.
  • 55. Cox DH, Aldrich RW. Role of the beta1 subunit in large-conductance Ca(2+)-activated K(+) channel gating energetics. Mechanisms of enhanced Ca(2+) sensitivity. J Gen Physiol. 2000 Sep;116(3):411-32.
  • 56. Schmidt D, MacKinnon R. Voltage-dependent K+ channel gating and voltage sensor toxin sensitivity depend on the mechanical state of the lipid membrane. Proc Natl Acad Sci U S A. 2008 Dec 9;105(49):19276-81.
  • 57. Nelson MT, Bonev AD. The beta1 subunit of the Ca2+-sensitive K+ channel protects against hypertension. J Clin Invest. 2004 Apr;113(7):955-7.
  • Availability of data and materials
    All data generated or analyzed during this study are included in this article.
  • Funding
    This study was supported by the National Natural Science Foundation of China (n. 82074162 and 81903995), Youth Talent Promotion Project of China Association of Chinese Medicine (n. CACM-2019-QNRC2-C10), Shanghai Municipal Commission of Health and Family Planning Fund (n. 20184Y0086 and 2018JQ003), Project for Capacity Promotion of Putuo District Clinical Special Disease (n. 2019tszb02), Science and Technology Innovation Project of Putuo District Health System (n. ptkwws201902, ptkwws202009 and ptkwws201908).
  • Ethics approval
    Not applicable.
  • Consent for publication
    Not applicable.

Publication Dates

  • Publication in this collection
    02 June 2021
  • Date of issue
    2021

History

  • Received
    08 Dec 2020
  • Accepted
    01 Mar 2021
location_on
Centro de Estudos de Venenos e Animais Peçonhentos (CEVAP/UNESP) Av. Universitária, 3780, Fazenda Lageado, Botucatu, SP, CEP 18610-034, Brasil, Tel.: +55 14 3880-7693 - Botucatu - SP - Brazil
E-mail: editorial.jvatitd@unesp.br
rss_feed Acompanhe os números deste periódico no seu leitor de RSS
Acessibilidade / Reportar erro