Next Article in Journal
IoT for Monitoring Fungal Growth and Ochratoxin A Development in Grapes Solar Drying in Tunnel and in Open Air
Next Article in Special Issue
Snake Venom: A Promising Source of Neurotoxins Targeting Voltage-Gated Potassium Channels
Previous Article in Journal
The Role of Heterocysts in Cyanotoxin Production during Nitrogen Limitation
Previous Article in Special Issue
Recombinant Production, NMR Solution Structure, and Membrane Interaction of the Phα1β Toxin, a TRPA1 Modulator from the Brazilian Armed Spider Phoneutria nigriventer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 15
Issue 10
10.3390/toxins15100612
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells

by
Ekaterina N. Lyukmanova
1,2,3,4,*,†,
Maxim M. Zaigraev
2,3,†,
Dmitrii S. Kulbatskii
2,†,
Aizek B. Isaev
2,3,
Ilya D. Kukushkin
2,3,
Maxim L. Bychkov
2,
Mikhail A. Shulepko
1,
Anton O. Chugunov
2,3 and
Mikhail P. Kirpichnikov
2,4
1
Faculty of Biology, MSU-BIT Shenzhen University, Shenzhen 518172, China
2
Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Miklukho-Maklaya 16/10, Moscow 117997, Russia
3
Phystech School of Biological and Medical Physics, Moscow Institute of Physics and Technology (National Research University), Institutsky Lane 9, Dolgoprudny, Moscow 141701, Russia
4
Interdisciplinary Scientific and Educational School of Moscow University «Molecular Technologies of the Living Systems and Synthetic Biology», Faculty of Biology, Lomonosov Moscow State University, Leninskie Gory, Moscow 119234, Russia
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Toxins 2023, 15(10), 612; https://doi.org/10.3390/toxins15100612
Submission received: 8 September 2023 / Revised: 30 September 2023 / Accepted: 7 October 2023 / Published: 13 October 2023
(This article belongs to the Special Issue Ion Channels, Venom, and Toxins)
Browse Figures
Versions Notes

Abstract

:
Cancer progression is characterized by microenvironmental acidification. Tumor cells adapt to low environmental pH by activating acid-sensing trimeric ion channels of the DEG/ENaC family. The α-ENaC/ASIC1a/γ-ENaC heterotrimeric channel is a tumor-specific acid-sensing channel, and its targeting can be considered a new strategy for cancer therapy. Mambalgin-2 from the Dendroaspis polylepis venom inhibits the α-ENaC/ASIC1a/γ-ENaC heterotrimer more effectively than the homotrimeric ASIC1a channel, initially proposed as the target of mambalgin-2. Although the molecular basis of such mambalgin selectivity remained unclear. Here, we built the models of the complexes of mambalgin-2 with the α-ENaC/ASIC1a/γ-ENaC and ASIC1a channels, performed MD and predicted the difference in the binding modes. The importance of the ‘head’ loop region of mambalgin-2 for the interaction with the hetero-, but not with the homotrimeric channel was confirmed by site-directed mutagenesis and electrophysiology. A new mode of allosteric regulation of the ENaC channels by linking the thumb domain of the ASIC1a subunit with the palm domain of the γ-ENaC subunit was proposed. The data obtained provide new insights into the regulation of various types of acid-sensing ion channels and the development of new strategies for cancer treatment.
Keywords:
three-finger toxins; acid-sensing ion channels; cancer; DEG/ENaC; ASIC1a; α-ENaC; γ-ENaC; mambalgin-2
Key Contribution: The acidification of tumor milieu drives cancer cell adaptation via tumor-specific acid-sensing heterotrimeric channels composed of the ASIC1a; α-ENaC, and γ-ENaC subunits from the DEG/ENaC family. Three-finger protein mambalgin-2 from Dendroaspis polylepis venom inhibits the heterotrimeric α-ENaC/ASIC1a/γ-ENaC more effectively than the homotrimeric ASIC1a channel due to the additional contacts of the mambalgin-2 ‘head’ loop with the γ-ENaC subunit. Targeting of the α-ENaC/ASIC1a/γ-ENaC heterotrimer by mambalgin-2 may become perspective strategy for cancer treatment.

Graphical Abstract

1. Introduction

The intensive metabolism is characteristic of tumors. Due to high cell density and insufficient oxygenation, metabolic pathways are activated in tumor cells, leading to the formation of acidic products that are removed from tumor cells to maintain a neutral intracellular pH [1,2,3]. This decreases the pH of the tumor microenvironment from ~7.4 to ~5.5–6.5 [4,5,6,7,8]. Tumor cells adapt to low environmental pH by activating pH sensors such as proton-sensitive cation channels of the degenerin/epithelial sodium channel (DEG/ENaC) family, which are trimers consisting of homologous subunits ASIC1a, ASIC1b, ASIC2a, ASIC2b, ASIC3, ASIC4, and α–δ ENaC [9,10,11,12,13]. The implication of these pH sensors in maintaining different (gastric, breast, hepatic, and glioblastoma) tumor cell growth is well documented [14,15,16,17,18,19]. The subunits that form these channels determine their pH sensitivity and pharmacological properties [20]. The most sensitive to a pH drop are the ASIC1a channels, which are expressed in the central nervous system and mediate synaptic plasticity and pain perception. Channels formed by ASIC1a subunits are expressed in non-neuronal cells, where they regulate gas exchange, the secretion of antibodies by lymphocytes, and the differentiation of keratinocytes in the skin [11]. ASIC1a channels are activated by a pH drop to ~6.4–6.6 and provide an entry of sodium ions into the cell. The inward sodium current regulates important cellular functions such as calcium homeostasis, volume and shape maintenance, activity of mitogenic signaling pathways based on activation of Wnt, ERK, or p38 MAP kinase, gene transcription (via β-catenin, the transcription factor NFκB, or by regulation of the maturation of microRNA miR-350) [11,21].
The ENaC channels are also DEG/ENaC family members and are activated in tumor cells by protease cleavage [22]. Due to their similar structure as well as the fundamental similarity of selective filters, the ASIC and ENaC subunits can form heterocomplexes [23]. The interaction between the α-ENaC and ASIC1a subunits, resulting in a cationic current in alveolar cells, has been shown previously in a hybrid system [24]. The heterocomplex formed by the α-ENaC, ASIC1a, and γ-ENaC subunits is characteristic of late-stage glioma cells but not of normal astrocytes [12,13] or lung adenocarcinoma cells [25].
Activated DEG/ENaC channels regulate the important signaling pathways that promote the proliferation and migration of tumor cells of different origins [16,26,27,28,29,30,31], thus targeting these channels is a perspective anti-tumor strategy. For example, inhibitors of the ASIC1a-containing channels—mambalgin-2 (Mamb-2), found in the black mamba (Dendroaspis polylepis) venom [32],—suppresses the growth of glioma, leukemia, melanoma, and lung adenocarcinoma cells [25,33,34,35]. Recently, it was demonstrated that Mamb-2 interacts with the tumor-specific heterotrimeric channel composed of the α-ENaC, ASIC1a, and γ-ENaC subunits in melanoma and lung adenocarcinoma cells [25,35]. Moreover, Mamb-2 inhibits the heterotrimeric α-ENaC/ASIC1a/γ-ENaC channel with a higher efficiency than the homotrimeric ASIC1a channel [25]. Although the molecular basis of such a toxin’s selectivity remains unknown.
Here, we studied the Mamb-2 interaction with the homotrimeric ASIC1a and heterotrimeric α-ENaC/ASIC1a/γ-ENaC channels by MD simulations and revealed that the toxin forms additional ionic and hydrophobic contacts with the complementary γ-ENaC(−) subunit. Site-directed mutagenesis of the residues of Mamb-2 and the γ-ENaC(−) subunit from the interacting interface confirmed the discovered mode of the toxin/channel interaction. Our data provide new information about the regulation of acid-sensing ion channels and the demand for rational design of new anti-cancer drugs targeting the tumor-specific α-ENaC/ASIC1a/γ-ENaC heterotrimer.

2. Results

Previously, significantly stronger inhibitory activity of Mamb-2 at the α-ENaC/ASIC1a/γ-ENaC heterotrimeric channels in comparison to the ASIC1 homotrimeric channels was reported [25]. Here, we aimed to reveal the molecular and structural basis for such a toxin’s selectivity. To achieve that, we built the models for hetero- and homotrimeric channels in complex with Mamb-2 and ran molecular dynamics (MD) simulations for them in the membrane environment. To prove the initially postulated modeling hypothesis, we performed site-directed mutagenesis and electrophysiology.

2.1. Modeling of the Mambalgin-2 Complexes with the ASIC1a and α-ENaC/ASIC1/γ-ENaC Channels

To build a model of the α-ENaC/ASIC1a/γ-ENaC channel, we considered the following points: (1) These channels are assembled counterclockwise; such assembly is provided by inter-subunit complementarity due to a network of hydrogen bonds [36]. Therefore, the α-ENaC, ASIC1a, and γ-ENaC subunits in the heterochannel should also be arranged accordingly. (2) Mamb-2 inhibits currents through the homotrimeric ASIC1a and heterotrimeric α-ENaC/ASIC1a/γ-ENaC channels upon a pH decrease [25]; thus, we proposed that Mamb-2 in complex with the heterotrimeric channel interacts at least with the ASIC1a subunit. (3) According to the cryo-EM structure of the homotrimeric ASIC1a channel with mambalgin-1 (Mamb-1), the closest Mamb-2 homologue (differing by a single substitution Y4F at the N-terminus), mambalgins interact mainly with the ASIC1a subunit located clockwise relative to Mamb-2 (view from the extracellular domain, in the direction of the channel pore), but not with the ASIC1a subunit located counterclockwise [37,38]. We labeled the subunit interacting with Mamb-2 in the hetero- and homocomplexes as the primary (+), and those located counterclockwise as the complementary (−) ones.
The models of the hetero- and homomeric channels were built in PyMOL (Figure 1) by superimposition of the individual subunit models on the experimental cryo-EM structure of the homomeric ASIC1a channel in complex with three Mamb-1 molecules whose positions were not altered (PDB code: 7CFT [37]). Structural models of the human ASIC1a (13-470), α-ENaC (50-593), and γ-ENaC subunits (20-579) with Uniprot codes P78348, P37088, and P51170, respectively, were obtained from the AlphaFold database [39,40]. To avoid unwanted artifacts, the disordered N- and C-termini residues were removed from the structures (the final length of the subunits used in the work is shown in the brackets in the previous sentence). The Mamb-2 structure was obtained from the Mamb-1 structure (pdb code 7CFT) by the Y4F substitution using the mutagenesis tool in PyMOL.
The models of the Mamb-2 complex with the heterotrimeric α-ENaC/ASIC1a/γ-ENaC and homotrimeric ASIC1a channels were built in a lipid bilayer composed of palmitoyloleoylphosphatidylcholine (POPC) molecules similar to [41] using the CHARMM-GUI [42] Membrane Builder tool [43]. To represent N-glycans at sites 368/395 and 209/497 at the ASIC1a and γ-ENac subunits, respectively, a classic pentasaccharide core composed of two N-Acetylglucosamine (GlcNAc) and three Mannose (Man) residues was used. N-glycans were introduced into the built models of the channels using the Glycan Reader and Modeller tool [44]. The resulting models are shown in Figure 1; further details are provided in Methods.
Molecular dynamics (MD) simulations for both modeled complexes were performed for 500 ns. In both models, Mamb-2 was bound at the adjacent subunits’ interfaces: ASIC1a/γ-ENaC or ASIC1a/ASIC1a (Figure 1). Due to the different number of toxin binding sites, one and three toxin molecules were bound to the heterotrimeric and homotrimeric channels, respectively. Similarly to the cryo-EM structure of ASIC1a/Mamb-1 (PDB code: 7CFT [37]), Mamb-2 interacted mainly with the thumb domain of the ASIC1a(+) subunits and with the lower palm domain of the complementary (−) subunits without contacting the membrane (Figure 1). No considerable contacts of Mamb-2 with N-glycans were observed. The general position of Mamb-2 relative to the subunits and cell membrane in both complexes, as well as the stability of the complexes, did not change during MD.

2.2. Mambalgin-2 Interacts More Extensively with the Complementary γ-ENaC(−)-Subunit as Compared with ASIC1(−)

The MD results with the subsequent intermolecular contact analysis showed that the toxin’s interaction sites at the primary ASIC1a(+) subunit were almost identical in the hetero- and homotrimeric channels. Mamb-2 interacted with the primary ASIC1a(+) subunit in both complexes mainly by the central (II) loop, although some contacts formed by the toxin’s loop I and ‘head’ also were revealed (Figure 2 and Figure 3c, and Tables S1 and S2). Significant differences in the interacting interfaces Mamb-2/channel were observed for the complementary γ-ENaC(−) and ASIC1a(−) subunits of the hetero- and homotrimeric channels, respectively (Figure 3a,b). Mamb-2 interacted with ASIC1a(−) very sparsely (via few hydrophobic contacts), while its contacts with γ-ENaC(−) were more numerous and not limited to just hydrophobic interactions but also included the ionic bridge D15–K91 (Figure 3a,b and Table 1).
In both cases, Mamb-2 is bound to the complementary subunit by just a tip of its “head”, involving mainly hydrophobic interactions. Although the number of such contacts was significantly higher for the α-ENaC/ASIC1a/γ-ENaC channel, this complements the pure hydrophobic interactions with the ionic/hydrogen bond between the toxin’s D15 and the channel’s K91 (Figure 3, Table 1). Notably, in both cases, the toxin’s binding surface was formed by the ‘head’ residues H13–M16 (Figure 3c). No significant staking or π-cation interactions were found in both complexes. An unusual feature arises in the Mamb-2/γ-ENaC interaction (Table 1): the ligand’s R14 interacts (although hydrophobically) with the receptor’s R90, which was weird given the presumed repulsion of these positively charged residues. Indeed, guanidine groups seemed to repel in MD, while aliphatic chains interacted. At the same time, such unusual Arg–Arg short-range interactions, stabilized by water and contributing to intermolecular recognition, were reported previously [45]. We suggest that in the Mamb-2/γ-ENaC interaction, a similar case can be realized.

2.3. Interaction with γ-ENaC(−) Determines the High Inhibitory Activity of Mamb-2 at α-ENaC/ASIC1a/γ-ENaC Heterotrimers

To validate the modeling results, we obtained the mutant variants of the γ-ENaC subunit with F89A, R90A, and K91A substitutions. The influence of these mutations on the activity of Mamb-2 at the α-ENaC/ASIC1a/γ-ENaC channel was analyzed in electrophysiological experiments in X. laevis oocytes, which express wild-type (WT) and mutant variants of the channel. Notably, the introduced mutations didn’t affect the mean expression level of the channel, as was shown by the comparison of peak current distributions (see Figure S1). We found that 370 nM of Mamb-2 almost completely inhibited the current through the WT α-ENaC/ASIC1a/γ-ENaC heterotrimer. However, in the case of the α-ENaC/ASIC1a/γ-ENaC mutants, Mamb-2 demonstrated lower activity (Figure 4a).
The analysis of “dose-response” curves showed that the mutations R90A and K91A in the γ-ENaC subunit significantly diminished the activity of Mamb-2 at the α-ENaC/ASIC1a/γ-ENaC heterotrimeric channel (Figure 4b, Table 2). The curves for R90A and K91A channel mutants moved right regarding the curve for the WT channel, demonstrating a significant increase in IC50 from ~130 nM of the WT heterotrimeric channel to ~250 nM and 190 nM, respectively. Simultaneously, the changes in the bottom of the curves (the value of maximal inhibitory response) for the R90A and K91A mutants of the α-ENaC/ASIC1a/γ-ENaC channel were observed, pointing to possible allosteric regulation by Mamb-2.

2.4. Mamb-2 “Head” Mutations Diminish Its Inhibitory Activity at the α-ENaC/ASIC1a/γ-ENaC Channel but Not at ASIC1a Channel

To prove the participation of the Mamb-2 ‘head’ in the interaction with the complementary γ-ENaC subunit revealed by the MD study (Figure 2, Table 1) and the greater role of this interaction in α-ENaC/ASIC1a/γ-ENaC as compared with the ASIC1a channel, we produced the toxin mutants with substitutions H13A, R14A, D15A, and R14A/D15A. Mutations did not alter the overall Mamb-2 structure, according to 1D 1H-NMR analysis. The NMR spectra of the mutants were very similar to the spectrum of WT Mamb-2, suggesting that mutant toxins’ structures did not change (Figure S2).
We found that no mutant toxins were able to fully inhibit the currents through the WT α-ENaC/ASIC1a/γ-ENaC channel at the concentration of 370 nM, while WT Mamb-2 (Figure 5a). For each Mamb-2 mutant, we obtained “dose–response” curves and compared them to the curve of the WT toxin. It appeared that every single mutation resulted in significantly diminished toxin activity at the α-ENaC/ASIC1a/γ-ENaC heterochannel (Figure 5b). IC50 for the mutants H13A, R14A, D15A, and R14A/D15A became ~350, 200, 230, and 220 nM, respectively, while the IC50 of WT Mamb-2 was found to be ~130 nM. Notably, the diminished activity of the double mutant R14A/D15A was obvious but didn’t reach statistical significance (Table 3). No significant changes in the bottom of the curves for any of the Mamb-2 mutants were revealed (Table 3).
At the same time, the ’head’ mutations H13A, R14A, and D15A caused no influence on the Mamb-2 activity at the ASIC1a homotrimeric channel (Figure 6a,b). In contrast to the situation with the α-ENaC/ASIC1a/γ-ENaC channel, the double mutation R14A/D15A and single mutation H13A resulted in the increased toxin’s activity at the ASIC1a channel with significantly enhanced maximal inhibitory response, although without significant change in IC50 (Figure 6a,b, Table 4). Thus, Mamb-2 WT and its mutants H13A, R14A, and D15A demonstrated only partial inhibitory activity at the ASIC1a channel even at high concentrations (10 µM), while the double mutation R14A/D15A resulted in nearly full inhibition of the ASIC1a channel at this concentration (Figure 6b, Table 4).

3. Discussion

Acid-sensing trimeric ion channels of the DEG/ENaC family are activated by external acidification (pH drop) and are important regulators of different essential processes such as mechanosensitivity, pain perception, synaptic plasticity, anxiety response, cancer progression, etc., [10,46]. The composition of the channel’s trimers determines its sensitivity to the pH drop or other physiological stimuli and shapes the adequate response to an external signal [11]. The formation of channels with altered composition and response patterns in comparison to physiologically normal channels may serve as a mechanism of adaptation of the cells to different pathological conditions. For example, during cancer progression, the formation of the tumor-specific heterotrimeric α-ENaC/ASIC1a/γ-ENaC channels was shown [12,13]. The activation of such channels may drive the adaptation of tumor cells to external acidification and activate different mitogenic signaling pathways [10]. Thus, targeting such cancer-specific heterotrimeric channels may be a perspective strategy for the treatment of different oncopathologies.
Previously, it was shown that mambalgin-2, the three-finger protein (TFP, [47]) from the Mamba venom initially described as the inhibitor of the ASIC1a channels [48], interacts with the heterotrimeric α-ENaC/ASIC1a/γ-ENaC channels in various cancer cells, and this interaction results in control of the homeostasis of cancer cells [25,35]. Moreover, mambalgin-2 inhibits currents evoked by the pH drop at the α-ENaC/ASIC1a/γ-ENaC channels with higher efficiency than at the ASIC1a channels (Table 3 and Table 4, and [25]). Here, we aimed to elucidate the molecular basis underlying the mambalgin-2 selectivity regarding the heterotrimeric α-ENaC/ASIC1a/γ-ENaC channels.
We performed computer modeling of the interaction of mambalgin-2 with two types of channels: heterotrimeric α-ENaC/ASIC1a/γ-ENaC and homotrimeric ASIC1a. We found that the toxin forms a similar interaction pattern with the primary ASIC1a(+) subunit in both cases, mostly by the central (II) loop (Figure 2), which agrees well with the previous mutagenesis data [48]. However, the interaction mode of the toxin with the complementary (−) subunits was different. Almost all residues forming the interaction interface Mamb-2/γ-ENaC(−) in complex with the heterotrimeric channel were located in the ‘head’ of the toxin’s molecule (Figure 3a). The ‘head’ region of Mamb-2 formed the set of hydrophobic and ionic contacts with the γ-ENaC(−) subunit, while the number of the ‘head’ contacts with the ASIC1a(−) subunit in complex with the homotrimeric channel was obviously lower (Figure 3a,b and Table 1, Tables S1 and S2). Site-directed mutagenesis with the subsequent electrophysiological studies confirmed the importance of the ‘head’ residues for the toxin’s interaction just with the α-ENaC/ASIC1a/γ-ENaC channels, but not with the ASIC1a channels. Thus, the stronger activity of Mamb-2 at the α-ENaC/ASIC1a/γ-ENaC channels can be explained by the existence of additional contacts of the toxin’s ‘head’ with the complementary γ-ENaC(−) subunit.
The main functional epitopes of TFPs are located in the disordered loops I-III [49,50,51,52]. The ‘Head’ region of TFPs also has a disordered structure (Figure 3c) and can be considered an auxiliary binding site regulating TFP’s interaction with other supplementary partners necessary for the function of the ‘main’ receptor. These supplementary partners can be either the cell membrane surrounding the nicotinic acetylcholine receptor [53,54] or other receptors, for example, receptor-tyrosine kinases in the case of nicotinic acetylcholine receptors and the human SLURP-1 protein in cancer cells [51]. Here, on the example of mambalgin-2, we revealed an additional ‘application’ of the ‘head’ region of TFPs: the regulation of the interaction selectivity between various types of acid-sensing channels. Probably, the ‘head’ moiety can be used for transplanting one or another activity to TFPs and to drug design with directed specificity.
Surprisingly, wild-type Mamb-2 fully inhibited the α-ENaC/ASIC1a/γ-ENaC channels and only partially inhibited the ASIC1a ones (Figure 4 and Figure 5), pointing to possible allosteric modulation of the homotrimeric channel by the toxin. Mutations R90A and K91A of the heterotrimeric channel led to increased IC50 (from 130 ± 20 to 250 ± 60 and 190 ± 30 nM, respectively) and the diminished maximal inhibitory activity of Mamb-2 (from 100% to ~70%, Figure 4, Table 2), indicating the appearance of the allosteric effects with switching of the pharmacology of the α-ENaC/ASIC1a/γ-ENaC channel to the properties of the homotrimeric ASIC1a channel. On the other hand, the introduction of the double mutation R14A/D15A into the Mamb-2 molecule also changed the mode of the toxin interaction with the homotrimeric channel in the direction of the heterotrimeric channel (increase in the maximal inhibitory activity of Mamb-2 from ~60% to ~85% with slight decrease in IC50, Figure 6, Table 4). Thus, the changes in the binding interface both from the toxin’s and channel’s sides could shift the pharmacology of heterotrimeric and homotrimeric acid-sensing channels to each other.
Previously, indirect data indicating that Mamb-2 inhibits and stabilizes the desensitization state of the channels containing the ASIC1a subunit in cancer cells were obtained [33,34]. Here, we built the models of the complexes of Mamb-2 with the α-ENaC/ASIC1a/γ-ENaC and ASIC1a channels, and these complexes were stable at MD during the 500 ns trajectory. According to these models, Mamb-2 interacts with the α4 and α5 helices of the thumb domain of the primary ASIC1a(+) subunit (Figure 1), which has an allosteric effect on the acid (proton-binding) pocket, keeping it in an expanded state with a simultaneously closed ion channel [38]. This channel’s state is characterized by a noticeable weakening of pH-sensing and a decrease in the current through the channel [38], which is in line with our present data (Figure 4, Figure 5 and Figure 6). On the other hand, Mamb-2 interacts with the residues from the β1-β2 linker of the palm domain of the complementary (−) subunits (F89, R90, K91 in the case of the γ-ENaC (−) subunit and V80, A81, A82, and S83 in the case of ASIC1a (−) subunit). Previously, in the example of the homotrimeric ASIC1a channel, the important role of reorientation of the β1-β2 linker in switching the channel from the low-pH open state to the low-pH desensitized state was shown [55]. This means that the interaction of Mamb-2 with the β1-β2 linker can stabilize either the open or desensitized states of the channels. Keeping in mind that no activation effect from Mamb-2 was observed (Figure 4, Figure 5 and Figure 6), this indicates that the toxin stabilizes the desensitization state of the α-ENaC/ASIC1a/γ-ENaC and ASIC1a channels.
The function of the α-ENaC/β-ENaC/γ-ENaC channels can be regulated by different molecules and factors such as extracellular Na+, Cl−, H+, metals, proteases, amiloride, and sheer stress [56]. Presently, several extracellular regulatory sites in the heterotrimeric channel were described: the extracellular sites located at the finger domains of the α-ENaC and γ-ENaC subunits, which can affect Na+ self-inhibition [22,57], and the amiloride binding pore domain [58]. Moreover, cross-linking of the palm domains of α-ENaC or β-ENaC with the thumb domains of γ-ENaC or α-ENaC, respectively, resulted in strong inhibition of the α-ENaC/β-ENaC/γ-ENaC channel activity [56]. Here, we substituted the β-ENaC subunit with the ASIC1a subunit and showed that linking the thumb domain of ASIC1a with the palm domain of γ-ENaC within the α-ENaC/ASIC1a/γ-ENaC channel by the Mamb-2 molecule (Figure 3a) also resulted in the channel inhibition (Figure 4). Thus, a new mode of regulation of the heterotrimeric ASIC1a/ENaC channels was revealed.

4. Conclusions

Thus, by computer modeling, we predicted the difference in the binding modes of mambalgin-2 with the heterotrimeric α-ENaC/ASIC1a/γ-ENaC and homotrimeric ASIC1a channels. The importance of the ‘head’ region of the toxin for the interaction with the heterotrimeric channel was confirmed by site-directed mutagenesis and electrophysiology. A new mode of allosteric regulation of ENaC channels was revealed. The data obtained provide new insights into the regulation of various types of acid-sensing ion channels and the development of new strategies for cancer treatment.

5. Materials and Methods

5.1. Computer Modeling Methods

5.1.1. Molecular Dynamics Calculations

MD trajectories for the heterotrimeric α-ENaC/ASIC1a/γ-ENaC and homotrimeric ASIC1a channels in complex with Mamb-2 were calculated in the all-atom CHARMM36m forcefield [59] and TIP3P water model via GROMACS [60]. For a more accurate description of π-cation interactions, an additional set of CHARMM36-WYF parameters [61] was used. MD calculations included the following stages: energy minimization (steepest descent algorithm), NVT relaxation (250 ps), NPT relaxation with C-rescale barostat [62] (2000 ps), and production MD (500 ns).
The MD calculations for an equilibrated system were performed in an NPT ensemble at a temperature of 310 K with a V-rescale thermostat [63] and a Parinello-Raman barostat [64] with a time step of 2 fs. No position restraints were applied to any molecules during the production MD phase, so the structures of the proteins and N-glycans were allowed to relax and change.

5.1.2. Contacts Analysis

For analysis of intermolecular contacts between Mamb-2 and the heterotrimeric α-ENaC/ASIC1a/γ-ENaC and homotrimeric ASIC1a channels, we utilized our in-house IMPULSE software package (cont_stat.js procedure) [65]. Ionic, ion-dipole, π-cationic, staking interactions, hydrogen bonds, and hydrophobic contacts were analyzed. Hydrophobic contacts were determined according to the concept of molecular hydrophobic potential (MHP) [66]. In the case of the homotrimeric ASIC1a channel data, the Mamb-2 interaction lifetime was averaged between three toxin molecules.
A complete list of all found contacts is available in Tables S1 and S2 for the hetero- and homotrimer channels, respectively. The values in the tables are the relative lifetimes of interactions between the Mamb-2 and channels’ primary (+) and complementary (−) subunit residues; one can analyze them using the sorting and filtering tools of the table processing program. Relative lifetime is the fraction of MD time when the contact is observed, from absence (0) to presence throughout the whole trajectory (1).
Using the obtained data on the contacts lifetimes from the Impulse software package, we built maps of the interaction intensity of Mamb-2 with the primary ASIC1a(+) subunit in the ASIC1a(+)/ASIC1a(+) and ASIC1a(+)/γ-ENaC(−) interfaces in the homotrimeric and heterotrimeric channels, respectively. To do that, we performed a summation of all the interaction lifetimes of different types of contacts between each residue of Mamb-2 and each residue of ASIC1a(+) in both complexes. The data on the contacts lifetime used for summation are represented in Tables S1 and S2, as are the initial csv files obtained from the Impulse software package after MD trajectory analysis. The summation of lifetimes was performed using our own Python script. The following types of contacts lifetimes were summarized: ionic bridges, hydrogen bonds, π-cation, stacking, and hydrophobic interactions. The resulting heatmaps were prepared via Matplotlib and Seaborn library utilities. Developed scripts for data analysis and contacts maps preparation are available at the Zenodo repository.
The following software versions were used: GROMACS 2022.4; CHARMM-GUI 3.8; IMPULSE 21.09; and PyMOL 2.5.5.

5.2. Recombinant Mambalgin-2 and Its Mutants Production

Recombinant Mamb-2 and its mutant variants with the H13A, R14A, D15A, and R14A/D15A substitutions were produced in E. coli as described previously [25]. The homogeneity and purity of the mambalgin-2 were confirmed by SDS-PAGE, HPLC, and MALDI-MS in the reaction, and Ellman’s reagent (Sigma-Aldrich, St. Louis, MO, USA) was used to confirm the disulfide bond formation. The proper spatial structure was confirmed by 1D 1H NMR.

5.3. Electrophysiological Recordings in X. laevis Oocytes

For expression of the human homotrimeric ASIC1a, heterotrimeric α-ENaC/ASIC1a/γ-ENaC, or heteromeric channel with F89A, R90A, and K91A mutations in the γ-ENAC subunit, Xenopus laevis oocytes were harvested and injected with synthesized mRNA coding the corresponding genes as described in [67]. Two-electrode voltage-clamp recordings were done using the TEC-03X amplifier (npi electronic GmbH, Tamm, Germany) at a holding potential of −50 mV. The output signal was filtered at 20 Hz by a four-pole Bessel filter and digitized at 1 kHz using the National Instruments BNC-6251 card, controlled by WinWCP 5.2.7. Microelectrodes were fabricated with 1.0–1.5 MOhm resistance and filled with 3 M KCl. The external perfusion solution for oocytes was ND-96 medium (96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 1 mM MgCl2, 10 mM HEPES, pH 7.4) for both ASIC1a and α-ENaC/ASIC1a/γ-ENaC, and it was fed to a perfusion chamber (volume ~50 µL) using a gravity-flow system at speed 2.5 mL/min. The currents were stimulated by a pH drop to 5.0 for both ASIC1a and α-ENaC/ASIC1a/γ-ENaC.
The pH drop was performed by fast perfusion solution exchange for 0.1 s to ND-96 medium with 10 mM MES instead of HEPES, adjusted to the target pH, and kept for 7 s; after that, the solution was exchanged to ND-96 with a pH of 7.4. The treatment of oocytes with the Mamb-2-containing solution was performed by solution exchange in the recording chamber for 15 s before the current stimulation. Lyophilized Mamb-2 was dissolved in 100% DMSO in a 10-mM stock solution and then diluted to a target concentration with ND-96 medium before the experiment. Amiloride (Sigma-Aldrich) was diluted in ND-96 medium (pH 7.4) from the 500-mM stock in 100% DMSO. The DMSO concentration in the recording solution did not exceed 1%. The exchange of extracellular solutions was performed using a computer-controlled valve system.
Mamb-2 concentrations for the dose–response curves were 0–10 μM. Each oocyte was tested with all concentrations of Mamb-2. Dose–response curves were fitted with the Hill equation: Y = 1/(1 + 10((LogIC50 − X)×nH)), where Y is the normalized current amplitude (independent normalization for each oocyte), X is the lg of concentration of the ligand concentration, and nH is the Hill coefficient (slope factor). Fit parameters were compared with one-way ANOVA and Dunnett’s multiple comparison test; for IC50, the comparison of LogIC50 was done instead of IC50 itself.

5.4. Statistical Analysis

Data are presented as mean ± SEM. The number of samples (n) is shown in the figure legends. The absence of outliers in each dataset was confirmed by the Grubbs’ test (alpha = 0.05). The data were tested for normality (Shapiro–Wilk test, at p = 0.05) and for the homogeneity of variances (Levene’s test, at p = 0.05). The data were analyzed using the F-test, as indicated in the figure legends. Differences in the data were considered statistically significant at p < 0.05. Analysis was performed using the GraphPad Prism 8.0 software (GraphPad Software).

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/toxins15100612/s1, Figure S1. Scatter plot of peak current amplitudes for α-ENaC/ASIC1a/γ-ENaC oocytes with WT, F89A, R90A, K91A γ-ENaC subunits. No statistically significant differences between means were detected by one-way ANOVA (F(3, 21) = 0.4985, p = 0.6873), Figure S2. Amide regions of 1H 1D NMR spectra of Mamb-2 variants (700 MHz, 30 °C, pH 4.0); Supplementary Tables S1 and S2 for contacts lifetime data, MD trajectories, Python scripts and CSV data can be downloaded at Zenodo repository https://doi.org/10.5281/zenodo.8303194 (accessed on 8 October 2023).

Author Contributions

Conceptualization, M.M.Z. and E.N.L.; methodology, M.M.Z., D.S.K. and E.N.L.; software, M.M.Z. and D.S.K.; validation, E.N.L., M.M.Z., A.O.C. and D.S.K.; formal analysis, M.M.Z., M.L.B. and D.S.K.; investigation, M.M.Z., D.S.K., A.B.I., I.D.K. and M.A.S.; resources, E.N.L., M.P.K. and A.O.C.; data curation, E.N.L. and A.O.C.; writing—original draft preparation, E.N.L., M.L.B., M.M.Z. and A.O.C.; writing—review and editing, D.S.K.; visualization, M.M.Z., A.O.C., E.N.L. and D.S.K.; supervision, E.N.L. and M.P.K.; project administration, E.N.L. and M.P.K.; funding acquisition, E.N.L. and M.P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the grant from the Ministry of Science and Higher Education of Russian Federation (project № 075-15-2020-773).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Putnam, R.W. 22—Intracellular pH Regulation. In Cell Physiology Source Book, 3rd ed.; Sperelakis, N., Ed.; Academic Press: San Diego, CA, USA, 2001; pp. 357–372. ISBN 978-0-12-656976-6. [Google Scholar]
  2. Dratkiewicz, E.; Simiczyjew, A.; Mazurkiewicz, J.; Ziętek, M.; Matkowski, R.; Nowak, D. Hypoxia and Extracellular Acidification as Drivers of Melanoma Progression and Drug Resistance. Cells 2021, 10, 862. [Google Scholar] [CrossRef]
  3. Böhme, I.; Bosserhoff, A.K. Acidic Tumor Microenvironment in Human Melanoma. Pigment Cell Melanoma Res. 2016, 29, 508–523. [Google Scholar] [CrossRef] [PubMed]
  4. Gillies, R.J.; Raghunand, N.; Karczmar, G.S.; Bhujwalla, Z.M. MRI of the Tumor Microenvironment. J. Magn. Reson. Imaging 2002, 16, 430–450. [Google Scholar] [CrossRef]
  5. Gerweck, L.E.; Seetharaman, K. Cellular pH Gradient in Tumor versus Normal Tissue: Potential Exploitation for the Treatment of Cancer. Cancer Res. 1996, 56, 1194–1198. [Google Scholar]
  6. Hjelmeland, A.B.; Wu, Q.; Heddleston, J.M.; Choudhary, G.S.; MacSwords, J.; Lathia, J.D.; McLendon, R.; Lindner, D.; Sloan, A.; Rich, J.N. Acidic Stress Promotes a Glioma Stem Cell Phenotype. Cell Death Differ. 2011, 18, 829–840. [Google Scholar] [CrossRef] [PubMed]
  7. Vaupel, P.; Kallinowski, F.; Okunieff, P. Blood Flow, Oxygen and Nutrient Supply, and Metabolic Microenvironment of Human Tumors: A Review. Cancer Res. 1989, 49, 6449–6465. [Google Scholar] [PubMed]
  8. Shi, Q.; Maas, L.; Veith, C.; Van Schooten, F.J.; Godschalk, R.W. Acidic Cellular Microenvironment Modifies Carcinogen-Induced DNA Damage and Repair. Arch. Toxicol. 2017, 91, 2425–2441. [Google Scholar] [CrossRef]
  9. Kellenberger, S.; Schild, L. Epithelial Sodium Channel/Degenerin Family of Ion Channels: A Variety of Functions for a Shared Structure. Physiol. Rev. 2002, 82, 735–767. [Google Scholar] [CrossRef] [PubMed]
  10. Kellenberger, S.; Schild, L. International Union of Basic and Clinical Pharmacology. XCI. Structure, Function, and Pharmacology of Acid-Sensing Ion Channels and the Epithelial Na+ Channel. Pharmacol. Rev. 2015, 67, 1–35. [Google Scholar] [CrossRef] [PubMed]
  11. Wemmie, J.A.; Taugher, R.J.; Kreple, C.J. Acid-Sensing Ion Channels in Pain and Disease. Nat. Rev. Neurosci. 2013, 14, 461–471. [Google Scholar] [CrossRef]
  12. Kapoor, N.; Lee, W.; Clark, E.; Bartoszewski, R.; McNicholas, C.M.; Latham, C.B.; Bebok, Z.; Parpura, V.; Fuller, C.M.; Palmer, C.A.; et al. Interaction of ASIC1 and ENaC Subunits in Human Glioma Cells and Rat Astrocytes. Am. J. Physiol. Cell Physiol. 2011, 300, C1246–C1259. [Google Scholar] [CrossRef]
  13. Kapoor, N.; Bartoszewski, R.; Qadri, Y.J.; Bebok, Z.; Bubien, J.K.; Fuller, C.M.; Benos, D.J. Knockdown of ASIC1 and Epithelial Sodium Channel Subunits Inhibits Glioblastoma Whole Cell Current and Cell Migration. J. Biol. Chem. 2009, 284, 24526–24541. [Google Scholar] [CrossRef] [PubMed]
  14. Sheng, Y.; Wu, B.; Leng, T.; Zhu, L.; Xiong, Z. Acid-Sensing Ion Channel 1 (ASIC1) Mediates Weak Acid-Induced Migration of Human Malignant Glioma Cells. Am. J. Cancer Res. 2021, 11, 997–1008. [Google Scholar]
  15. Jin, C.; Ye, Q.-H.; Yuan, F.-L.; Gu, Y.-L.; Li, J.-P.; Shi, Y.-H.; Shen, X.-M.; Liu, B.; Lin, Z.-H. Involvement of Acid-Sensing Ion Channel 1α in Hepatic Carcinoma Cell Migration and Invasion. Tumour Biol. 2015, 36, 4309–4317. [Google Scholar] [CrossRef] [PubMed]
  16. Zhu, S.; Zhou, H.-Y.; Deng, S.-C.; Deng, S.-J.; He, C.; Li, X.; Chen, J.-Y.; Jin, Y.; Hu, Z.-L.; Wang, F.; et al. ASIC1 and ASIC3 Contribute to Acidity-Induced EMT of Pancreatic Cancer through Activating Ca2+/RhoA Pathway. Cell Death Dis. 2017, 8, e2806. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, X.; Sun, X.; Wang, Z.; Zhou, X.; Xu, L.; Li, F.; Zhang, X.; Pan, J.; Qi, L.; Qian, H.; et al. Involvement of Acid-Sensing Ion Channel 1a in Gastric Carcinoma Cell Migration and Invasion. Acta Biochim. Biophys. Sin. 2018, 50, 440–446. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, B.; Liu, J.; Ho, T.-T.; Ding, X.; Mo, Y.-Y. ERK-Mediated NF-κB Activation through ASIC1 in Response to Acidosis. Oncogenesis 2016, 5, e279. [Google Scholar] [CrossRef]
  19. Harguindey, S.; Stanciu, D.; Devesa, J.; Alfarouk, K.; Cardone, R.A.; Polo Orozco, J.D.; Devesa, P.; Rauch, C.; Orive, G.; Anitua, E.; et al. Cellular Acidification as a New Approach to Cancer Treatment and to the Understanding and Therapeutics of Neurodegenerative Diseases. Semin. Cancer Biol. 2017, 43, 157–179. [Google Scholar] [CrossRef]
  20. Osmakov, D.I.; Khasanov, T.A.; Andreev, Y.A.; Lyukmanova, E.N.; Kozlov, S.A. Animal, Herb, and Microbial Toxins for Structural and Pharmacological Study of Acid-Sensing Ion Channels. Front. Pharmacol. 2020, 11, 991. [Google Scholar] [CrossRef] [PubMed]
  21. Cullinan, M.M.; Klipp, R.C.; Bankston, J.R. Regulation of Acid-Sensing Ion Channels by Protein Binding Partners. Channels 2021, 15, 635–647. [Google Scholar] [CrossRef]
  22. Kleyman, T.R.; Eaton, D.C. Regulating ENaC’s Gate. Am. J. Physiol. Cell Physiol. 2020, 318, C150–C162. [Google Scholar] [CrossRef] [PubMed]
  23. Hanukoglu, I. ASIC and ENaC Type Sodium Channels: Conformational States and the Structures of the Ion Selectivity Filters. FEBS J. 2017, 284, 525–545. [Google Scholar] [CrossRef]
  24. Trac, P.T.; Thai, T.L.; Linck, V.; Zou, L.; Greenlee, M.; Yue, Q.; Al-Khalili, O.; Alli, A.A.; Eaton, A.F.; Eaton, D.C. Alveolar Nonselective Channels Are ASIC1a/α-ENaC Channels and Contribute to AFC. Am. J. Physiol. Lung Cell. Mol. Physiol. 2017, 312, L797–L811. [Google Scholar] [CrossRef] [PubMed]
  25. Sudarikova, A.V.; Bychkov, M.L.; Kulbatskii, D.S.; Chubinskiy-Nadezhdin, V.I.; Shlepova, O.V.; Shulepko, M.A.; Koshelev, S.G.; Kirpichnikov, M.P.; Lyukmanova, E.N. Mambalgin-2 Inhibits Lung Adenocarcinoma Growth and Migration by Selective Interaction with ASIC1/α-ENaC/γ-ENaC Heterotrimer. Front. Oncol. 2022, 12, 904742. [Google Scholar] [CrossRef]
  26. Gupta, S.C.; Singh, R.; Asters, M.; Liu, J.; Zhang, X.; Pabbidi, M.R.; Watabe, K.; Mo, Y.-Y. Regulation of Breast Tumorigenesis through Acid Sensors. Oncogene 2016, 35, 4102–4111. [Google Scholar] [CrossRef] [PubMed]
  27. Wu, Y.; Gao, B.; Xiong, Q.-J.; Wang, Y.-C.; Huang, D.-K.; Wu, W.-N. Acid-Sensing Ion Channels Contribute to the Effect of Extracellular Acidosis on Proliferation and Migration of A549 Cells. Tumour Biol. 2017, 39, 1010428317705750. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, C.; Zhu, L.-L.; Xu, S.-G.; Ji, H.-L.; Li, X.-M. ENaC/DEG in Tumor Development and Progression. J. Cancer 2016, 7, 1888–1891. [Google Scholar] [CrossRef] [PubMed]
  29. Xu, S.; Liu, C.; Ma, Y.; Ji, H.-L.; Li, X. Potential Roles of Amiloride-Sensitive Sodium Channels in Cancer Development. Biomed. Res. Int. 2016, 2016, 2190216. [Google Scholar] [CrossRef]
  30. Zhang, Y.; Cao, N.; Gao, J.; Liang, J.; Liang, Y.; Xie, Y.; Zhou, S.; Tang, X. ASIC1a Stimulates the Resistance of Human Hepatocellular Carcinoma by Promoting EMT via the AKT/GSK3β/Snail Pathway Driven by TGFβ/Smad Signals. J. Cell. Mol. Med. 2022, 26, 2777–2792. [Google Scholar] [CrossRef] [PubMed]
  31. Rooj, A.K.; McNicholas, C.M.; Bartoszewski, R.; Bebok, Z.; Benos, D.J.; Fuller, C.M. Glioma-Specific Cation Conductance Regulates Migration and Cell Cycle Progression. J. Biol. Chem. 2012, 287, 4053–4065. [Google Scholar] [CrossRef]
  32. Diochot, S.; Baron, A.; Salinas, M.; Douguet, D.; Scarzello, S.; Dabert-Gay, A.-S.; Debayle, D.; Friend, V.; Alloui, A.; Lazdunski, M.; et al. Black Mamba Venom Peptides Target Acid-Sensing Ion Channels to Abolish Pain. Nature 2012, 490, 552–555. [Google Scholar] [CrossRef] [PubMed]
  33. Bychkov, M.L.; Shulepko, M.A.; Vasileva, V.Y.; Sudarikova, A.V.; Kirpichnikov, M.P.; Lyukmanova, E.N. ASIC1a Inhibitor Mambalgin-2 Suppresses the Growth of Leukemia Cells by Cell Cycle Arrest. Acta Naturae 2020, 12, 101–116. [Google Scholar] [CrossRef] [PubMed]
  34. Bychkov, M.; Shulepko, M.; Osmakov, D.; Andreev, Y.; Sudarikova, A.; Vasileva, V.; Pavlyukov, M.S.; Latyshev, Y.A.; Potapov, A.A.; Kirpichnikov, M.; et al. Mambalgin-2 Induces Cell Cycle Arrest and Apoptosis in Glioma Cells via Interaction with ASIC1a. Cancers 2020, 12, 1837. [Google Scholar] [CrossRef] [PubMed]
  35. Bychkov, M.L.; Kirichenko, A.V.; Shulepko, M.A.; Mikhaylova, I.N.; Kirpichnikov, M.P.; Lyukmanova, E.N. Mambalgin-2 Inhibits Growth, Migration, and Invasion of Metastatic Melanoma Cells by Targeting the Channels Containing an ASIC1a Subunit Whose Up-Regulation Correlates with Poor Survival Prognosis. Biomedicines 2021, 9, 1324. [Google Scholar] [CrossRef] [PubMed]
  36. Noreng, S.; Bharadwaj, A.; Posert, R.; Yoshioka, C.; Baconguis, I. Structure of the Human Epithelial Sodium Channel by Cryo-Electron Microscopy. Elife 2018, 7, e39340. [Google Scholar] [CrossRef] [PubMed]
  37. Sun, D.; Liu, S.; Li, S.; Zhang, M.; Yang, F.; Wen, M.; Shi, P.; Wang, T.; Pan, M.; Chang, S.; et al. Structural Insights into Human Acid-Sensing Ion Channel 1a Inhibition by Snake Toxin Mambalgin1. Elife 2020, 9, e57096. [Google Scholar] [CrossRef]
  38. Salinas, M.; Kessler, P.; Douguet, D.; Sarraf, D.; Tonali, N.; Thai, R.; Servent, D.; Lingueglia, E. Mambalgin-1 Pain-Relieving Peptide Locks the Hinge between A4 and A5 Helices to Inhibit Rat Acid-Sensing Ion Channel 1a. Neuropharmacology 2021, 185, 108453. [Google Scholar] [CrossRef] [PubMed]
  39. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A.; et al. AlphaFold Protein Structure Database: Massively Expanding the Structural Coverage of Protein-Sequence Space with High-Accuracy Models. Nucleic Acids Res. 2022, 50, D439–D444. [Google Scholar] [CrossRef] [PubMed]
  40. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef]
  41. Mernea, M.; Ulăreanu, R.Ș.; Cucu, D.; Al-Saedi, J.H.; Pop, C.-E.; Fendrihan, S.; Anghelescu, G.D.C.; Mihăilescu, D.F. Epithelial Sodium Channel Inhibition by Amiloride Addressed with THz Spectroscopy and Molecular Modeling. Molecules 2022, 27, 3271. [Google Scholar] [CrossRef]
  42. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef] [PubMed]
  43. Lee, J.; Patel, D.S.; Ståhle, J.; Park, S.-J.; Kern, N.R.; Kim, S.; Lee, J.; Cheng, X.; Valvano, M.A.; Holst, O.; et al. CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans. J. Chem. Theory Comput. 2019, 15, 775–786. [Google Scholar] [CrossRef] [PubMed]
  44. Park, S.-J.; Lee, J.; Qi, Y.; Kern, N.R.; Lee, H.S.; Jo, S.; Joung, I.; Joo, K.; Lee, J.; Im, W. CHARMM-GUI Glycan Modeler for Modeling and Simulation of Carbohydrates and Glycoconjugates. Glycobiology 2019, 29, 320–331. [Google Scholar] [CrossRef] [PubMed]
  45. Magalhaes, A.; Maigret, B.; Hoflack, J.; Gomes, J.N.; Scheraga, H.A. Contribution of Unusual Arginine-Arginine Short-Range Interactions to Stabilization and Recognition in Proteins. J. Protein Chem. 1994, 13, 195–215. [Google Scholar] [CrossRef]
  46. Wemmie, J.A.; Price, M.P.; Welsh, M.J. Acid-Sensing Ion Channels: Advances, Questions and Therapeutic Opportunities. Trends Neurosci. 2006, 29, 578–586. [Google Scholar] [CrossRef] [PubMed]
  47. Vasilyeva, N.A.; Loktyushov, E.V.; Bychkov, M.L.; Shenkarev, Z.O.; Lyukmanova, E.N. Three-Finger Proteins from the Ly6/uPAR Family: Functional Diversity within One Structural Motif. Biochemistry 2017, 82, 1702–1715. [Google Scholar] [CrossRef] [PubMed]
  48. Mourier, G.; Salinas, M.; Kessler, P.; Stura, E.A.; Leblanc, M.; Tepshi, L.; Besson, T.; Diochot, S.; Baron, A.; Douguet, D.; et al. Mambalgin-1 Pain-Relieving Peptide, Stepwise Solid-Phase Synthesis, Crystal Structure, and Functional Domain for Acid-Sensing Ion Channel 1a Inhibition. J. Biol. Chem. 2016, 291, 2616–2629. [Google Scholar] [CrossRef] [PubMed]
  49. Bourne, Y.; Taylor, P.; Marchot, P. Acetylcholinesterase Inhibition by Fasciculin: Crystal Structure of the Complex. Cell 1995, 83, 503–512. [Google Scholar] [CrossRef] [PubMed]
  50. Wu, M.; Liu, C.Z.; Joiner, W.J. Structural Analysis and Deletion Mutagenesis Define Regions of QUIVER/SLEEPLESS That Are Responsible for Interactions with Shaker-Type Potassium Channels and Nicotinic Acetylcholine Receptors. PLoS ONE 2016, 11, e0148215. [Google Scholar] [CrossRef]
  51. Bychkov, M.L.; Shulepko, M.A.; Shlepova, O.V.; Kulbatskii, D.S.; Chulina, I.A.; Paramonov, A.S.; Baidakova, L.K.; Azev, V.N.; Koshelev, S.G.; Kirpichnikov, M.P.; et al. SLURP-1 Controls Growth and Migration of Lung Adenocarcinoma Cells, Forming a Complex With A7-nAChR and PDGFR/EGFR Heterodimer. Front. Cell Dev. Biol. 2021, 9, 739391. [Google Scholar] [CrossRef] [PubMed]
  52. Lyukmanova, E.N.; Shenkarev, Z.O.; Schulga, A.A.; Ermolyuk, Y.S.; Mordvintsev, D.Y.; Utkin, Y.N.; Shoulepko, M.A.; Hogg, R.C.; Bertrand, D.; Dolgikh, D.A.; et al. Bacterial Expression, NMR, and Electrophysiology Analysis of Chimeric Short/Long-Chain α-Neurotoxins Acting on Neuronal Nicotinic Receptors. J. Biol. Chem. 2007, 282, 24784–24791. [Google Scholar] [CrossRef]
  53. Shenkarev, Z.O.; Chesnokov, Y.M.; Zaigraev, M.M.; Chugunov, A.O.; Kulbatskii, D.S.; Kocharovskaya, M.V.; Paramonov, A.S.; Bychkov, M.L.; Shulepko, M.A.; Nolde, D.E.; et al. Membrane-Mediated Interaction of Non-Conventional Snake Three-Finger Toxins with Nicotinic Acetylcholine Receptors. Commun. Biol. 2022, 5, 1344. [Google Scholar] [CrossRef] [PubMed]
  54. Lesovoy, D.M.; Bocharov, E.V.; Lyukmanova, E.N.; Kosinsky, Y.A.; Shulepko, M.A.; Dolgikh, D.A.; Kirpichnikov, M.P.; Efremov, R.G.; Arseniev, A.S. Specific Membrane Binding of Neurotoxin II Can Facilitate Its Delivery to Acetylcholine Receptor. Biophys. J. 2009, 97, 2089–2097. [Google Scholar] [CrossRef]
  55. Yoder, N.; Yoshioka, C.; Gouaux, E. Gating Mechanisms of Acid-Sensing Ion Channels. Nature 2018, 555, 397–401. [Google Scholar] [CrossRef]
  56. Zhang, L.; Wang, X.; Chen, J.; Sheng, S.; Kleyman, T.R. Extracellular Intersubunit Interactions Modulate Epithelial Na+ Channel Gating. J. Biol Chem. 2023, 299, 102914. [Google Scholar] [CrossRef] [PubMed]
  57. Winarski, K.L.; Sheng, N.; Chen, J.; Kleyman, T.R.; Sheng, S. Extracellular Allosteric Regulatory Subdomain within the Gamma Subunit of the Epithelial Na+ Channel. J. Biol. Chem. 2010, 285, 26088–26096. [Google Scholar] [CrossRef]
  58. Kashlan, O.B.; Kleyman, T.R. ENaC Structure and Function in the Wake of a Resolved Structure of a Family Member. Am. J. Physiol. Ren. Physiol. 2011, 301, F684–F696. [Google Scholar] [CrossRef]
  59. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmüller, H.; MacKerell, A.D. CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins. Nat. Methods 2017, 14, 71–73. [Google Scholar] [CrossRef] [PubMed]
  60. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef]
  61. Khan, H.M.; MacKerell, A.D.; Reuter, N. Cation-π Interactions between Methylated Ammonium Groups and Tryptophan in the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2019, 15, 7–12. [Google Scholar] [CrossRef]
  62. Bernetti, M.; Bussi, G. Pressure Control Using Stochastic Cell Rescaling. J. Chem. Phys. 2020, 153, 114107. [Google Scholar] [CrossRef]
  63. Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef]
  64. Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  65. Krylov, N.A.; Efremov, R.G. Libxtc: An Efficient Library for Reading XTC-Compressed MD Trajectory Data. BMC Res. Notes 2021, 14, 124. [Google Scholar] [CrossRef]
  66. Pyrkov, T.V.; Chugunov, A.O.; Krylov, N.A.; Nolde, D.E.; Efremov, R.G. PLATINUM: A Web Tool for Analysis of Hydrophobic/Hydrophilic Organization of Biomolecular Complexes. Bioinformatics 2009, 25, 1201–1202. [Google Scholar] [CrossRef]
  67. Osmakov, D.I.; Koshelev, S.G.; Palikov, V.A.; Palikova, Y.A.; Shaykhutdinova, E.R.; Dyachenko, I.A.; Andreev, Y.A.; Kozlov, S.A. Alkaloid Lindoldhamine Inhibits Acid-Sensing Ion Channel 1a and Reveals Anti-Inflammatory Properties. Toxins 2019, 11, 542. [Google Scholar] [CrossRef]
Toxins 15 00612 g001
Figure 1. Models of the complexes of Mamb-2 with ASIC1 homotrimer (a,b) and α-ENaC/ASIC1/γ-ENaC heterotrimer (c,d): side (a,c) and top (b,d) views are shown. Both models are immersed in the phospholipid bilayer and solvated to perform MD calculations; the presented structures result from 500 ns MD.
Figure 1. Models of the complexes of Mamb-2 with ASIC1 homotrimer (a,b) and α-ENaC/ASIC1/γ-ENaC heterotrimer (c,d): side (a,c) and top (b,d) views are shown. Both models are immersed in the phospholipid bilayer and solvated to perform MD calculations; the presented structures result from 500 ns MD.
Toxins 15 00612 g001
Toxins 15 00612 g002
Figure 2. Contact maps of the Mamb-2 interaction with the primary ASIC1a(+) subunit in heterotrimeric and homotrimeric complexes. Map cells are colored according to the contact intensity, which corresponds to the summarized lifetime of the contacts of all types between each pair of residues under consideration during the MD simulation of 500 ns. The following types of contacts were considered: ionic bridges, hydrogen bonds, π-cation, stacking, and hydrophobic interactions. ASIC1a residues outside the 280–370 range did not give a significant contribution to the interaction and therefore were not included in the map. The position of the loops’ and ‘head’ regions of Mamb-2 is shown.
Figure 2. Contact maps of the Mamb-2 interaction with the primary ASIC1a(+) subunit in heterotrimeric and homotrimeric complexes. Map cells are colored according to the contact intensity, which corresponds to the summarized lifetime of the contacts of all types between each pair of residues under consideration during the MD simulation of 500 ns. The following types of contacts were considered: ionic bridges, hydrogen bonds, π-cation, stacking, and hydrophobic interactions. ASIC1a residues outside the 280–370 range did not give a significant contribution to the interaction and therefore were not included in the map. The position of the loops’ and ‘head’ regions of Mamb-2 is shown.
Toxins 15 00612 g002
Toxins 15 00612 g003
Figure 3. Amino-acid residues implicated in Mamb-2 interaction with the complementary (−) subunits in the α-ENaC/ASIC1a/γ-ENaC hetero- (a) and ASIC1a homotrimer (b). Note that interactions with (−)-subunits determine the difference between these two complexes, while the constant part of the Mamb-2/ASIC1a(+) binding interface is not shown for clarity (see also Figure 2). The binding sites Mamb-2/(−)-subunit are shown by dotted lines. The contacting residues from the (−)-subunits are shown in italics. (c). Mamb-2 model: loops and ‘head’ are subscribed, along with receptor-interacting residues in the latter.
Figure 3. Amino-acid residues implicated in Mamb-2 interaction with the complementary (−) subunits in the α-ENaC/ASIC1a/γ-ENaC hetero- (a) and ASIC1a homotrimer (b). Note that interactions with (−)-subunits determine the difference between these two complexes, while the constant part of the Mamb-2/ASIC1a(+) binding interface is not shown for clarity (see also Figure 2). The binding sites Mamb-2/(−)-subunit are shown by dotted lines. The contacting residues from the (−)-subunits are shown in italics. (c). Mamb-2 model: loops and ‘head’ are subscribed, along with receptor-interacting residues in the latter.
Toxins 15 00612 g003
Toxins 15 00612 g004
Figure 4. Mamb-2 inhibitory activity at the α-ENaC/ASIC1a/γ-ENaC channel with mutated γ-ENaC(−) subunit (F89A, R90A and K91A). (a) Representative current traces recorded in X. laevis the wild-type (WT) α-ENaC/ASIC1a/γ-ENaC channel or its mutants with substitutions F89A, R90A, and K91A in the γ-ENaC subunit after pH decrease without (black lines) and with (red lines) application of 370 nM of Mamb-2. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the stimulation phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 inhibitory effect at WT and mutated α-ENaC/ASIC1a/γ-ENaC channels. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ± SEM (n = 4–7). The fitted curves are described by Hill’s equation. Parameters describing the curves’ fit are presented in Table 2. (c) Protocol of experiments within the two-electrode patch-clamp technique.
Figure 4. Mamb-2 inhibitory activity at the α-ENaC/ASIC1a/γ-ENaC channel with mutated γ-ENaC(−) subunit (F89A, R90A and K91A). (a) Representative current traces recorded in X. laevis the wild-type (WT) α-ENaC/ASIC1a/γ-ENaC channel or its mutants with substitutions F89A, R90A, and K91A in the γ-ENaC subunit after pH decrease without (black lines) and with (red lines) application of 370 nM of Mamb-2. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the stimulation phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 inhibitory effect at WT and mutated α-ENaC/ASIC1a/γ-ENaC channels. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ± SEM (n = 4–7). The fitted curves are described by Hill’s equation. Parameters describing the curves’ fit are presented in Table 2. (c) Protocol of experiments within the two-electrode patch-clamp technique.
Toxins 15 00612 g004
Toxins 15 00612 g005
Figure 5. Inhibitory activity of Mamb-2 and its mutants H13A, R14A, D15A, and R14A/D15A at the α-ENaC/ASIC1a/γ-ENaC channel. (a) Representative traces of currents recorded in X. laevis oocytes, which express WT α-ENaC/ASIC1a/γ-ENaC channel after pH decrease without (black lines) and with (red lines) application of 370 nM of Mamb-2 and its mutants. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the pH drop phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 and its mutants at the α-ENaC/ASIC1a/γ-ENaC channel. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ±SEM (n = 4–7). The fitted curves are described by Hill’s equation. Parameters describing the curves’ fit are presented in Table 3.
Figure 5. Inhibitory activity of Mamb-2 and its mutants H13A, R14A, D15A, and R14A/D15A at the α-ENaC/ASIC1a/γ-ENaC channel. (a) Representative traces of currents recorded in X. laevis oocytes, which express WT α-ENaC/ASIC1a/γ-ENaC channel after pH decrease without (black lines) and with (red lines) application of 370 nM of Mamb-2 and its mutants. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the pH drop phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 and its mutants at the α-ENaC/ASIC1a/γ-ENaC channel. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ±SEM (n = 4–7). The fitted curves are described by Hill’s equation. Parameters describing the curves’ fit are presented in Table 3.
Toxins 15 00612 g005
Toxins 15 00612 g006
Figure 6. Inhibitory activity of Mamb-2 and its mutants H13A, R14A, D15A, and R14A/D15A at the ASIC1a channel. (a) Representative traces of currents recorded in X. laevis oocytes, which express WT ASIC1a channel after pH decrease without (black lines) and with (red lines) application of 1.1 µM of Mamb-2 and its mutants. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the pH drop phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 and its mutants at the ASIC1a channel. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ±SEM (n = 4–7). The fitted curves are described by the single-component Hill’s equation. Parameters describing the curves’ fit are presented in Table 4.
Figure 6. Inhibitory activity of Mamb-2 and its mutants H13A, R14A, D15A, and R14A/D15A at the ASIC1a channel. (a) Representative traces of currents recorded in X. laevis oocytes, which express WT ASIC1a channel after pH decrease without (black lines) and with (red lines) application of 1.1 µM of Mamb-2 and its mutants. Oocytes were pre-incubated with Mamb-2 for 15 s (red bar above the trace, the beginning of pre-incubation is off the time scale), and the pH drop phase (pH 5.0) was 7 s. Time and current scale are shown by bars near the traces. (b) Dose–response curves for the Mamb-2 and its mutants at the ASIC1a channel. The response current was normalized to the control experiment for each Mamb-2 concentration. Each data point represents an average from independent experiments in different oocytes ±SEM (n = 4–7). The fitted curves are described by the single-component Hill’s equation. Parameters describing the curves’ fit are presented in Table 4.
Toxins 15 00612 g006
Table 1. Mambalgin-2 contacts with the complementary γ-ENaC(–) and ASIC1(–) subunits found in MD. The interaction type is denoted in the brackets: hydrogen bond (H), ionic bridge or ion–dipole interaction (I), and hydrophobic contact (M) (according to “Molecular Hydrophobic Potential”). H and I contacts with a lifetime >10% (of 500 ns MD trajectory) and M contacts with a lifetime >30% are listed. Hydrophobic (M) contacts between apparently polar/charged residues denote the interaction of non-polar fragments of these residues.
Table 1. Mambalgin-2 contacts with the complementary γ-ENaC(–) and ASIC1(–) subunits found in MD. The interaction type is denoted in the brackets: hydrogen bond (H), ionic bridge or ion–dipole interaction (I), and hydrophobic contact (M) (according to “Molecular Hydrophobic Potential”). H and I contacts with a lifetime >10% (of 500 ns MD trajectory) and M contacts with a lifetime >30% are listed. Hydrophobic (M) contacts between apparently polar/charged residues denote the interaction of non-polar fragments of these residues.
Mambalgin-2γ-ENaC(–)ASIC1a(–)
H13F89 (M)A81 (M)
R90 (M)
K91 (M)
R14F89 (M)V80 (M)
R90 (M)A81 (M)
K91 (M)
D15K91 (I, H, M)
M16K91 (M)A81 (M)
S298 (M)A82 (M)
M299 (M)S83 (M)
S39K91 (M)
Table 2. Parameters of Mamb-2 inhibitory activity at the α-ENaC/ASIC1a/γ-ENaC channel with the mutated γ-ENAC subunit. ‘Bottom’ is the value of maximal inhibitory response, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. ** (p < 0.01) and **** (p < 0.0001) indicate significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
Table 2. Parameters of Mamb-2 inhibitory activity at the α-ENaC/ASIC1a/γ-ENaC channel with the mutated γ-ENAC subunit. ‘Bottom’ is the value of maximal inhibitory response, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. ** (p < 0.01) and **** (p < 0.0001) indicate significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
γ-ENaCIC50, nMBottomnH
WT130 ± 200.07 ± 0.041.6 ± 0.4
F89A150 ± 20 n.s.0.14 ± 0.03 n.s.2.5 ± 0.7 n.s.
R90A250 ± 60 ****0.18 ± 0.07 n.s.1.8 ± 0.6 n.s.
K91A190 ± 30 **0.26 ± 0.05 n.s.2.1 ± 0.6 n.s.
Table 3. Parameters of inhibitory activity of the Mamb-2 mutants at the α-ENaC/ASIC1a/γ-ENaC channel. ‘Bottom’ is the value of the maximal inhibitory response, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. * (p < 0.05), **** (p < 0.0001) indicates significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
Table 3. Parameters of inhibitory activity of the Mamb-2 mutants at the α-ENaC/ASIC1a/γ-ENaC channel. ‘Bottom’ is the value of the maximal inhibitory response, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. * (p < 0.05), **** (p < 0.0001) indicates significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
Mambalgin-2 MutationIC50, nMBottomnH
WT130 ± 200.07 ± 0.04−1.6 ± 0.4
H13A350 ± 60 ****0.07 ± 0.05 n.s.−1.4 ± 0.3 n.s.
R14A200 ± 20 ****0.10 ± 0.04 n.s.−2.2 ± 0.3 *
D15A230 ± 30 ****0.09 ± 0.05 n.s.−2.0 ± 0.4 n.s.
R14A/D15A220 ± 40 ****0.03 ± 0.08 n.s.−1.2 ± 0.2 n.s.
Table 4. Parameters of inhibitory activity of the Mamb-2 mutants at the ASIC1a channel. ‘Bottom’ is the value of maximal inhibition, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. * (p < 0.05), ** (p < 0.01), and **** (p < 0.0001) indicate significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
Table 4. Parameters of inhibitory activity of the Mamb-2 mutants at the ASIC1a channel. ‘Bottom’ is the value of maximal inhibition, nH—the Hill coefficient, IC50—concentration of the ligand at which the amplitude of the current through the channel was inhibited by 50%. * (p < 0.05), ** (p < 0.01), and **** (p < 0.0001) indicate significant differences between the data groups and the wild-type group (WT) by one-way ANOVA with Dunnett post-hoc test, n.s. (p > 0.05) indicates the absence of statistical difference.
Mambalgin-2 MutationIC50, nMBottomnH
WT310 ± 600.39 ± 0.042.9 ± 1.5
H13A510 ± 1500 n.s.0.28 ± 0.05 **1.5 ± 0.4 *
R14A250 ± 300 n.s.0.44 ± 0.07 n.s.1.7 ± 0.9 n.s.
D15A300 ± 400 n.s.0.42 ± 0.06 n.s.1.0 ± 0.35 **
R14A/D15A470 ± 70 n.s.0.13 ± 0.05 ****1.9 ± 0.4 n.s.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lyukmanova, E.N.; Zaigraev, M.M.; Kulbatskii, D.S.; Isaev, A.B.; Kukushkin, I.D.; Bychkov, M.L.; Shulepko, M.A.; Chugunov, A.O.; Kirpichnikov, M.P. Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells. Toxins 2023, 15, 612. https://doi.org/10.3390/toxins15100612

AMA Style

Lyukmanova EN, Zaigraev MM, Kulbatskii DS, Isaev AB, Kukushkin ID, Bychkov ML, Shulepko MA, Chugunov AO, Kirpichnikov MP. Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells. Toxins. 2023; 15(10):612. https://doi.org/10.3390/toxins15100612

Chicago/Turabian Style

Lyukmanova, Ekaterina N., Maxim M. Zaigraev, Dmitrii S. Kulbatskii, Aizek B. Isaev, Ilya D. Kukushkin, Maxim L. Bychkov, Mikhail A. Shulepko, Anton O. Chugunov, and Mikhail P. Kirpichnikov. 2023. "Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells" Toxins 15, no. 10: 612. https://doi.org/10.3390/toxins15100612

APA Style

Lyukmanova, E. N., Zaigraev, M. M., Kulbatskii, D. S., Isaev, A. B., Kukushkin, I. D., Bychkov, M. L., Shulepko, M. A., Chugunov, A. O., & Kirpichnikov, M. P. (2023). Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells. Toxins, 15(10), 612. https://doi.org/10.3390/toxins15100612

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop