TRPV6 is a calcium-selective ion channel implicated in epithelial Ca2+ uptake. TRPV6 inhibitors are needed for the treatment of a broad range of diseases associated with disturbed calcium homeostasis, including cancers. Here we combine cryo-EM, calcium imaging, and mutagenesis to explore molecular bases of human TRPV6 inhibition by the antifungal drug econazole and the universal ion channel blocker ruthenium red (RR). Econazole binds to an allosteric site at the channel’s periphery, where it replaces a lipid. In contrast, RR inhibits TRPV6 by binding in the middle of the ion channel’s selectivity filter and plugging its pore like a bottle cork. Despite different binding site locations, both inhibitors induce similar conformational changes in the channel resulting in closure of the gate formed by S6 helices bundle crossing. The uncovered molecular mechanisms of TRPV6 inhibition can guide the design of a new generation of clinically useful inhibitors.
TRPV6 is a representative of the vanilloid subfamily of transient receptor potential (TRP) channels that is considered the principal epithelial calcium channel for Ca2+ absorption in intestines and many other organs1. TRPV6 is highly selective to Ca2+ (PCa/PNa > 100)2,3,4,5 and its constitutive activity is regulated positively by membrane lipids, such as phosphatidylinositol 4,5-bisphosphate (PIP2)6,7,8, and negatively by Mg2+, which is at least partially responsible for the strong inward rectification displayed by these channels9,10, and calmodulin (CaM) that causes Ca2+-dependent inactivation11,12,13 by plugging the channel through a unique cation-π interaction by inserting the side chain of lysine K115 into a tetra-tryptophan cage at the pore’s intracellular entrance14. Dysregulation of TRPV6 leads to disturbed calcium homeostasis. Thus, knockout TRPV6–/– mice exhibit defective intestinal Ca2+ absorption, increased urinary Ca2+ excretion, decreased femoral bone mineral density, lower body weight, alopecia, dermatitis, and severely impaired male fertility15,16,17,18,19, whereas changes in TRPV6 expression have been shown in several mouse models of human diseases20,21,22,23,24. Decreased mRNA and protein levels of TRPV6 have been found in placentas of women suffering from pre-eclampsia25. Several human TRPV6 gene mutations have been recently linked to transient neonatal hyperparathyroidism, skeletal under-mineralization and dysplasia, and chronic pancreatitis26,27,28,29,30,31,32. Because of the central role of calcium in cancer development33, TRPV6 has also been classified as an oncochannel that is involved in both increased cell proliferation and inhibition of apoptosis34, and overexpressed in various types of human cancer, including breast, prostate, colon, ovarian, thyroid, endometrial cancers, and leukemia33,35. TRPV6 inhibitors can therefore be used for the treatment of a broad range of diseases associated with disturbed calcium homeostasis, including TRPV6-rich tumors36,37,38,39.
Several small-molecule inhibitors of TRPV6 have been identified, including 2-aminoethoxydiphenyl borate (2-APB) and TH-1177, their derivatives, highly potent and selective derivatives of (4-phenylcyclohexyl)piperazine (PCHPDs), trinuclear ruthenium amine ruthenium red (RR), and antifungal drugs econazole and miconazole2,40,41,42,43,44,45,46,47. Structures of TRPV6 identified a 2-APB-binding site at the cytoplasmic end of the S1–S4 transmembrane helical bundle and suggested an allosteric mechanism of channel closure through modulation of protein–lipid interactions48. Several structures of TRPV6 in complex with PCHPDs identified their main binding site at the intracellular entrance to the ion channel pore, which overlaps with the binding site of CaM49. Therefore, PCHPDs act as TRPV6 open-channel blockers that convert the channel into a non-conducting state, mimicking inactivation by CaM14,49.
Here we use cryo-EM, calcium imaging, and mutagenesis to study structural mechanisms of TRPV6 inhibition by econazole and RR. RR, which is used in many biological applications50,51,52, inhibits a broad spectrum of ion channels, including ryanodine receptors, mitochondrial calcium uniporters, calcium, K2P, TRP, CALHM, and Piezo channels53,54. RR appears to act as a channel blocker of TRPV6 but compared with PCHPDs, its binding site is located on the opposite side of the ion channel pore, in the selectivity filter. In contrast, econazole acts as an allosteric inhibitor of TRPV6, similar to 2-APB. However, econazole binds between S1-S4 bundle of one TRPV6 subunit and S5 of the neighboring subunit, at the transmembrane domain (TMD) location that is different from the 2-APB-binding site48. Nevertheless, econazole replaces a lipid that is otherwise present at the binding site and induces a similar closure of the ion channel gate as do RR and 2-APB. Our results expand the spectrum of inhibitor binding sites and molecular mechanisms of inhibition that can be explored for the design of new drugs targeting TRPV6.
Structure of human TRPV6 in the open state
To study molecular details of TRPV6 inhibition, we first optimized our protein preparation. As a control, we aimed to determine the open (apo) state structure of human TRPV6 (TRPV6Open) at a higher resolution than that of the previously solved structures49,55. To avoid possible interference of inhibition with CaM-mediated inactivation, we used the C-terminally truncated functional human TRPV6 construct (hTRPV6-CtD) that lacks the CaM-binding site49. Initially, we improved the structural resolution for hTRPV6-CtD purified in glyco-diosgenin (GDN) detergent from 3.26 Å49 to 2.69 Å (Fig. 1a, Supplementary Fig. 1, Supplementary Table 1). As an alternative to the detergent, we reconstituted hTRPV6-CtD into cNW11 nanodiscs. This further improved the resolution to 2.54 Å (Supplementary Fig. 1, Supplementary Table 1), resulting in a structure that is essentially identical to the one in detergent (Supplementary Fig. 2), although with slightly fewer and less well resolved annular lipids. These two hTRPV6-CtD structures were also very similar to TRPV6 structures in amphipols and MSP2N2 nanodiscs (Supplementary Fig. 2), suggesting that TRPV6 structure determination does not depend on the type of membrane mimic or the presence or absence of the C-terminus and yields similar open-state conformations, albeit with a slightly better resolution in cNW11 nanodiscs.
The new open-state structures have the same architecture as previously published TRPV6 structures (Fig. 1b, Supplementary Fig. 2)14,48,55,56,57. The channel is assembled of four subunits and contains two main components: a TMD with a central ion channel pore and an intracellular skirt that is mostly built of ankyrin repeat domains connected to each other by three-stranded β-sheets, N-terminal helices, and C-terminal hooks. Amphipathic TRP helices run nearly parallel to the membrane and interact with both the TMD and the skirt. The TMD is composed of six transmembrane helices S1–S6 and a re-entrant pore loop (P-loop) between S5 and S6. A bundle of the first four transmembrane helices represents the S1–S4 domain. A domain homologous to the S1–S4 domain acts as a voltage sensor in voltage-gated ion channels58. The pore domain of each subunit includes S5, P-loop, and S6, and is packed against the S1–S4 domain of the neighboring subunit in a domain-swapped arrangement55,57.
Structures of TRPV6 in complex with RR and econazole
As cNW11 nanodiscs yielded the open-state structure at better resolution than in detergent, amphipols, or MSP2N2 nanodiscs, especially for the protein part, we used cNW11 nanodiscs to determine TRPV6 structures in complex with the inhibitors. Indeed, the resolution of structures obtained in the presence of RR (TRPV6RR) and econazole (TRPV6Eco) was about the same as for the open-state structure in cNW11 nanodiscs (Fig. 2, Supplementary Fig. 1, Supplementary Table 1). Compared with the open-state structure, the structure solved in the presence of RR showed a distinct non-protein density in the middle of the ion channel pore that perfectly matched the shape of RR (Fig. 2a, c). On the contrary, the structure of TRPV6 determined in the presence of econazole did not show densities in the ion channel pore that would resemble a molecule of econazole. Instead, it had four densities (one per subunit) in the TMD region peripheral to the ion channel pore with the shape that closely matched an econazole molecule (Fig. 2b, d). The location of these densities coincides with the location of lipid in other structures (Supplementary Fig. 3).
Binding site of RR in TRPV6 channel pore
Different parts of channel pores serve as binding sites for ion channel blockers. In TRPV6, two types of ion channel blockers have been characterized structurally, gadolinium56,57 and PCHPDs49. Although PCHPDs bind at the intracellular entrance to the ion channel pore, the gadolinium site is located at the extracellular entrance and formed by aspartates D542. The same site is shared by permeable ions56,57,59 and calcium binds at this site in the open state (Fig. 3a, d). RR binding occurs close to the extracellular pore entrance and its site partially overlaps with the binding site for ions (Fig. 3b, e). High-resolution reconstruction of TRPV6RR allows precise placement of the inhibitor inside the ion channel pore (Fig. 3b). On the extracellular side, the RR-binding site is flanked by aspartates D542 with a spherical density on the top that might represent a calcium ion typically bound at D542 in the absence of RR. The inhibitor occupies the entire selectivity filter of TRPV6, which is lined by carboxyl groups of D542, backbone carbonyl oxygens of I541, I540, and T539 as well as the hydroxyl group of T539, and spans the upper half of the pore’s central cavity. The surface of the selectivity filter is highly electronegative (Fig. 3e), creating a favorable environment for the positively charged RR that carries a 6+ total charge. Similar to TRPV6Open structure (Fig. 3a, d), the structure of TRPV6Eco shows no RR-like density in the selectivity filter but a strong density at the D542 site, presumably representing Ca2+ (Fig. 3c, f). There is also a noticeable density in the central cavity of TRPV6Eco pore. The shape of this density did not match the shape of econazole, even after multiple rounds of focused classification and symmetry relaxation. Nevertheless, as TRPV6Open and TRPV6RR structures do not contain any density in the central cavity above the noise level, there is a possibility that the observed density in the TRPV6Eco structure represents a low-affinity or non-specific binding of the inhibitor in random orientations that cannot be revealed by cryo-EM-processing tools.
Allosteric binding site of econazole
The econazole-binding site is located at the interface between S5 of one subunit and S4 of the neighboring subunit (Fig. 4a, c). Econazole is sandwiched between W495 of S5 and F472 of S4 and is also surrounded by hydrophobic side chains of L496 and V499 of S5 as well as M466, A469, M474, and L475 of S4. It is not surprising that the hydrophobic character of the econazole-binding site creates a favorable environment for lipids, which bind at this site in the absence of econazole (Supplementary Fig. 3; Fig. 4b, d).
We probed the econazole-binding site by mutating residues that contribute to this site and compared the concentration dependence of econazole inhibition of the corresponding mutants with wild-type channels (Fig. 4e). Econazole inhibited Ca2+ uptake through wild-type TRPV6 with IC50 of 4.4 ± 0.3 µM (n = 9). Substantial weakening of inhibition observed for L475A (IC50 = 151 ± 20 µM, n = 3), F472A (3.14 ± 0.24 mM, n = 6) and W495A (IC50 = 3.43 ± 0.52 mM, n = 6) mutants strongly supported the identified binding site of econazole. Mutation T539A in the ion channel pore, which is distal to the econazole-binding site, did not change the efficacy of calcium uptake inhibition by econazole (IC50 = 3.1 ± 0.2 µM, n = 3), confirming that the effects of L475A, F472A, and W495A mutations were specific. The strong effects of W495A and F472A mutations on econazole inhibition also indicate that the binding site at the S4–S5 interface is the main inhibitory site of econazole. Therefore, if density observed in the pore’s central cavity of TRPV6Eco (Fig. 3c) represents econazole molecule(s), it can only justify non-specific binding or low-affinity inhibition.
Effects of RR and econazole on TRPV6 gating
To explore the possible effects of RR and econazole on TRPV6 gating, we compared the ion channel pore in TRPV6Open with the pore in TRPV6RR and TRPV6Eco (Fig. 5). Consistent with previous findings, the pore of the open-state TRPV6 was in the conducting state (Fig. 5a, d). Typical for the open conformation of TRPV6, S6 was not entirely α-helical and contained a π-bulge resulting from an α-to-π transition in its middle (Figs. 5a, 3a). In contrast, the pores in TRPV6RR and TRPV6Eco were apparently closed due to a hydrophobic seal imposed by the gate residues L574 and M578 (Figs. 5b–d, 3e, f). Correspondingly, S6 in TRPV6RR and TRPV6Eco was entirely α-helical, typical for the closed conformation of TRPV6. The open-to-closed state transformation of the lower pore in TRPV6Open compared with TRPV6RR and TRPV6Eco included a ~100° rotation and bending away from the pore axis of the intracellular part of S6, as it was described previously48,55. This conformational change in S6 started below the gating hinge alanine A566 and did not cause alteration in the upper pore of TRPV6, leaving the selectivity filter apparently intact (Figs. 3 and 5).
By studying structural bases of TRPV6 inhibition by RR and econazole we identified two novel TRPV6 inhibitory binding sites: in the ion channel selectivity filter (Fig. 3) and at the interface between transmembrane helices S4 and S5 (Fig. 4). When added to the previously discovered inhibitory sites of gadolinium formed by the side chains of aspartates D542 at the pore extracellular entrance56,57, 2-APB in the pocket formed by the cytoplasmic half of the S1–S4 transmembrane helix bundle48, and PCHPDs at the pore intracellular entrance49, they increase the number of target sites for potential therapeutic interventions as well as synergistic treatment approaches. The number of distinct ligand binding sites in the vanilloid subfamily of TRP channels (TRPVs) therefore reaches twelve (Supplementary Fig. 4). Given the similarity in amino-acid sequences (Supplementary Fig. 5) and structural architecture of TRPVs, there is a good chance that all these sites can be targeted in every member of the subfamily. Some of these sites, like the vanilloid binding site initially identified in TRPV160, are already found to bind ligands in different TRPVs61,62. Other sites might only bind ligands in certain TRPVs and their specificity will require additional studies.
Previously, binding sites of RR were identified structurally in CALHM channels54 and K2P potassium channels53. In undecameric CALHM channels, RR binds near the pore-lining helix S1 and was proposed to stabilize this helix in a lifted conformation, giving rise to channel inhibition54. In K2P channels, RR binds under the CAP domain archway above the channel pore and was proposed to act according to a “finger in the dam” mechanism53. The two proposed mechanisms are fundamentally different from the mechanism of TRPV6 inhibition, where RR acts like a cork in the bottle by plugging the pore in the narrow region of the selectivity filter (Fig. 6). This region of the pore is likely optimal for RR block owing to its unique geometry as well as the overall negative surface charge that complements the positive charge of the blocker (Fig. 3).
Previous structural determination of econazole-binding sites in ion channels is limited to TRPV5 channels, which are highly homologous to TRPV6 (Supplementary Fig. 5). Based on the TRPV5 structure determined in the presence of econazole62, econazole was proposed to bind to the vanilloid site, which in TRPV1 accommodates vanilloid ligands like agonists resiniferatoxin and capsaicin, competitive antagonist capsazepine, or phosphatidylinositol (PI) lipid60,63,64. The vanilloid binding site is different from the econazole-binding site in TRPV6 identified in the present study. All residues coordinating econazole at its binding site in TRPV6 are conserved between TRPV5 and TRPV6, except A469, which is a threonine in TRPV5 (Supplementary Fig. 5). It is therefore unclear why econazole binding in these two highly homologous channels would be different. In all high-resolution structures of TRPV6, the vanilloid site is clearly occupied by cholesteryl hemisuccinate (CHS), which is routinely used during TRPV6 purification49 (Supplementary Fig. 3), and likely accommodates cholesterol in cellular conditions. It is therefore possible that in a relatively low, 4.8-Å resolution reconstruction of TRPV5, the blob of density assigned to the molecule of econazole62 represents a lipid. Higher-resolution structural information about TRPV5 inhibition by econazole would be helpful to better understand the difference from TRPV6. Of other inhibitory binding sites in TRP channels, the location of the econazole-binding site in TRPV6 (Fig. 4) has the highest similarity with the location of a piperlongumine-binding site in TRPV265. Compared with the econazole-binding site in TRPV6, the piperlongumine-binding site in TRPV2 is wedged deeper between S4 and S5 and makes contact with S6.
Compared with RR, which directly occludes the channel pore and prevents ionic conductance, econazole acts as an allosteric inhibitor. Econazole binds to the site, which is occupied by a lipid in the apo-state structure (Fig. 4). Because of the smaller size of econazole compared to the lipid, econazole binding creates a void that allows Q473 and M474 side chains to move closer to F472 and W495 and correspondingly away from R589 (Fig. 4f). Separation of Q473 and R589 results in breakage of the salt bridge connecting these residues in the open state55. The loss of this salt bridge, which energetically compensates for the unfavorable α-to-π transition in S6, reverses the transition. Concurrently, the lower portions of S6 helices rotate by ~100°, leading to the separation of D489 and T581 and the loss of the open-state-stabilizing hydrogen bond between them55. The rotated lower portions of S6 expose the side chains of L574 and M578 towards the pore center and hydrophobically seal the pore to prevent ion conductance, thus converting the channel into the closed, non-conducting state (Figs. 3, 5, and 6).
What remains a puzzle is why RR, which binds at the selectivity filter and physically occludes this region of the pore for ion conductance, also causes the gating transition in the lower pore similar to econazole. The selectivity filter conformations in TRPV6RR, TRPV6Eco, and TRPV6Open are nearly identical, as well as of the P-helix that directly contacts S6. We hypothesize that the positive charge of RR creates an electric field inside the pore’s central cavity that interacts with the electric dipole of the S6 helix (Fig. 6a). This interaction causes repulsion of the lower portion of S6 away from RR, which results in partial S6 rotation and helps S6 to become completely α-helical, with side chains of L574 and M578 sealing the pore (Fig. 6a, d). Interestingly, all currently known TRPV6 inhibitors either directly block the ion channel pore (Ga3+, calmodulin, PCHPDs, RR) or produce an allosteric effect by displacing lipids (2-APB, econazole). More work and high-resolution structures will be necessary to reveal novel druggable sites on the surface of TRPV6 and better understand the mechanisms and energetics of TRPV6 inhibition.
With RR and econazole joining PCHPDs and 2-APB in the list of TRPV6 inhibitors with structurally identified binding sites, there is a stronger foundation for the structure-based design and development of new therapeutic agents potentially capable of targeting a broad spectrum of TRPV6-linked diseases, including various types of human cancer. Importantly, since all of these compounds occupy distinct binding sites, they can likely bind simultaneously, resulting in synergistic inhibition of TRPV6. Thus, co-administration of molecules representing these different classes of inhibitors might be a future strategy for a more efficient, disease-dependent, and fine-tuned personalized therapy approach. The close homology and sequence similarity between TRPV6 and TRPV5 (Supplementary Fig. 5) suggest that these approaches can also be extended to TRPV5-linked diseases, including hypercalciuria, nephrolithiasis, and bone disorders66,67.
C-terminally truncated human TRPV6 (hTRPV6-CtD, residues 1–666 of wild-type channel)49 used for cryo-EM was cloned into a pEG BacMam vector68 with a C-terminal thrombin cleavage site followed by a streptavidin affinity tag (WSHPQFEK). For Fura-2 AM measurements, point mutations in wild-type human TRPV6 were introduced using standard molecular biology techniques (see the list of primers in Supplementary Table 2).
Expression and purification
hTRPV6-CtD was expressed and purified based on our previously developed protocol49. Bacmids and baculoviruses were produced using standard procedures68. Baculovirus was made in Sf9 cells for ~72 hours (Thermo Fisher Scientific, mycoplasma test negative, GIBCO #12659017) and was added to suspension-adapted HEK 293 cells lacking N-acetyl-glucosaminyltransferase I (GnTI–, mycoplasma test negative, ATCC #CRL-3022) that were maintained in Freestyle 293 media (Gibco-Life Technologies #12338-018) supplemented with 2% FBS at 37 °C and 5% CO2. Twenty-four hours after transduction, 10 mM sodium butyrate was added to enhance protein expression, and the temperature was reduced to 30 °C. Seventy-two hours after transduction, cells were harvested by centrifugation at 5471 × g for 15 min using a Sorvall Evolution RC centrifuge (Thermo Fisher Scientific), washed in phosphate-buffered saline pH 8.0, and pelleted by centrifugation at 3202 × g for 10 min using an Eppendorf 5810 centrifuge. The cell pellet was resuspended in ice-cold lysis buffer containing 20 mM Tris-Cl pH 8.0, 150 mM NaCl, 0.8 μM aprotinin, 4.3 μM leupeptin, 2 μM pepstatin A, 1 mM phenylmethylsulfonyl fluoride, and 1 mM β-mercaptoethanol (βME). Cells were subsequently lysed using a Misonix Sonicator with a preset program (six cycles of 15 s “on” at the amplitude of eight followed by 15 s “off”; this program was repeated three times for optimal cell lysis) under constant stirring on ice. Unbroken cells and cell debris were pelleted in an Eppendorf 5810 centrifuge at 3202 × g and 4 °C for 10 min. The supernatant was subjected to ultracentrifugation in a Beckman Coulter ultracentrifuge using a Beckman Coulter Type 45Ti rotor at 186,000 g and 4 °C for 1 hour to pellet the membranes. The membrane pellet was mechanically homogenized and solubilized in the lysis buffer supplemented with 1% (w/v) n-dodecyl β-d-maltoside and 0.1% (w/v) CHS under stirring at 4 °C for 2 hours. Insoluble material was removed by ultracentrifugation for 40 min in a Beckman Coulter Type 45Ti rotor at 186,000 × g and the supernatant was added to strep resin and rotated for 14–16 hours at 4 °C. The resin was washed with 10 column volumes of wash buffer containing 20 mM Tris-HCl pH 8.0, 150 mM NaCl, 1 mM βME, 0.01% (w/v) GDN, and 0.001% (w/v) CHS, and the protein was eluted with the same buffer supplemented with 2.5 mM D-desthiobiotin. The eluted protein was concentrated using a 100 kDa NMWL centrifugal filter (MilliporeSigma™ Amicon™) to 0.5 ml and then centrifuged in a Sorvall MTX 150 Micro-Ultracentrifuge (Thermo Fisher Scientific) using an S100AT4 rotor for 30 min at 66,000 × g and 4 °C before injection into a size-exclusion chromatography (SEC) column. All constructs were further purified using a Superose™ 6 10/300 GL SEC column attached to an AKTA FPLC (GE Healthcare) and equilibrated in 150 mM NaCl, 20 mM Tris-HCl pH 8.0, 1 mM βME, 0.01% GDN, and 0.001% CHS. The tetrameric peak fractions were pooled and concentrated using 100 kDa NMWL centrifugal filter to 2–4 mg/ml.
Circularized cNW11 nanodiscs were prepared according to the standard protocol69 and stored (~2–3 mg/ml) before usage at –80 °C in 20 mM Tris (pH 8.0) and 150 mM NaCl. For nanodisc reconstitution, the purified protein was mixed with circularized cNW11 nanodiscs and soybean lipids (Soy polar extract, Avanti Polar Lipids, USA) at a molar ratio of 1:3:166 (TRPV6 monomer:cNW11:lipid). The lipids were dissolved to a concentration of 100 mg/ml in 150 mM NaCl and 20 mM Tris (pH 8), and subjected to 5–10 cycles of freezing in liquid nitrogen and thawing in a water bath sonicator. For TRPV6RR, 1 mM RR was added to the TRPV6 monomer:cNW11:lipid mixture. The nanodisc mixture (500 µl) was rocked at room temperature for 1 h. Subsequently, 20 mg of Bio-beads SM2 (Bio-Rad), pre-wet in a buffer containing 20 mM Tris pH 8, 150 mM NaCl and 1 mM βME, were added to the nanodisc mixture, which was then rotated for one hour at 4 °C. Additional 20 mg of Bio-beads SM2 were added and the resulting mixture was rotated at 4 °C for another ~14–20 hours. Bio-beads SM2 were then removed by pipetting and nanodisc-reconstituted hTRPV6-CtD was purified from empty nanodiscs by SEC using Superose™ 6 10/300 GL SEC column equilibrated in 150 mM NaCl, 20 mM Tris (pH 8.0), and 1 mM βME. Fractions of nanodisc-reconstituted hTRPV6-CtD were pooled and concentrated using 100 kDa NMWL centrifugal filter to 0.6–3.0 mg/ml.
RR was dissolved in water to a final concentration of 50 mM. Econazole was dissolved in dimethyl sulfoxide to a final concentration of 100 mM. These stocks were added to the cNW11-reconstituted protein at the final concentrations of 0.5 mM econazole or 1 mM RR, and the resulting samples were incubated for two hours at room temperature before grid preparation.
Cryo-EM sample preparation and data collection
UltrAuFoil R 1.2/1.3, Au 300 grids were used for plunge-freezing. Prior to sample application, grids were plasma treated in a PELCO easiGlow glow discharge cleaning system (0.39 mBar, 15 mA, “glow” 25 s, “hold” 10 s). A Mark IV Vitrobot (Thermo Fisher Scientific) set to 100% humidity at 4 °C was used to plunge-freeze the grids in liquid ethane after applying 3 µl of protein sample to their gold-coated side using a blot time of 5 s, a blot force of 5, and a wait time of 15 s. To achieve optimal particle distribution for TRPV6RR, 3 µl of the protein sample was applied twice. The grids were stored in liquid nitrogen before imaging.
Images of frozen-hydrated particles of TRPV6Open in GDN and CHS were collected on a Titan Krios transmission electron microscope (TEM) (Thermo Fisher Scientific) operating at 300 kV and equipped with a post-column GIF Quantum energy filter and a Gatan K3 Summit direct electron detection (DED) camera (Gatan, Pleasanton, CA, USA) using EPU software (Thermo Fisher). 7904 micrographs were collected in counting mode with an image pixel size of 0.855 Å across a defocus range of –1.0 to –2.0 µm. The total dose of ~58 e− Å−2 was attained by using a dose rate of ~14.1 e− pixel−1 s−1 across 50 frames for a 3.0-s total exposure time.
Images of frozen-hydrated particles of cNW11-reconstituted TRPV6Open were collected on a Titan Krios TEM operating at 300 kV with a post-column GIF Quantum energy filter of 20 eV and a Gatan K3 Summit DED camera using SerialEM70. 5706 micrographs were collected in super-resolution mode with an image pixel size of 0.413 Å across a defocus range of –0.8 to –2.0 µm. The total dose of ~60 e− Å−2 was attained by using a dose rate of ~16 e− pixel−1 s−1 across 50 frames for a 2.0-s total exposure time.
Images of frozen-hydrated particles of cNW11-reconstituted TRPV6RR were collected on a Titan Krios TEM operating at 300 kV with a post-column GIF Quantum energy filter of 20 eV and a Gatan K3 Summit DED camera using SerialEM. 4431 micrographs were collected in super-resolution mode with an image pixel size of 0.4 Å across a defocus range of –0.5 to –2.0 µm. The total dose of ~58 e− Å−2 was attained by using a dose rate of ~16 e− pixel−1 s−1 across 50 frames for a 2.0-s total exposure time.
Images of frozen-hydrated particles of cNW11-reconstituted TRPV6Eco were collected on a Titan Krios TEM operating at 300 kV with a post-column GIF Quantum energy filter of 20 eV and a Gatan K3 Summit DED camera using SerialEM. 3950 micrographs were collected in counting mode with an image pixel size of 0.86 Å across a defocus range of –0.8 to –2.0 µm. The total dose of ~60 e− Å−2 was attained by using a dose rate of ~16 e− pixel−1 s−1 across 50 frames for a 2.2-s total exposure time.
Image processing and 3D reconstruction
All data sets were processed in RELION71 and cryoSPARC72 (see Supplementary Table 1 for details). Movies were motion-corrected with MotionCor273 algorithm implemented in RELION. Contrast transfer function (CTF) estimation was performed on non-dose-weighted micrographs using Gctf74. Subsequent data processing was done on dose-weighted micrographs. Following CTF estimation, micrographs were manually inspected and those with outliers in defocus values, ice thickness, and astigmatism as well as micrographs with lower predicted CTF-correlated resolution (>6 Å) were excluded from the rest of the processing pipeline (individually assessed for each parameter relative to overall distribution; no set threshold). Particles were picked using a Laplacian-of-Gaussian approach or with the help of unpublished TRPV6 cryo-EM density templates in RELION and 2D-classified in RELION or cryoSPARC. A selection of 2D classes was used to generate templates to be used for template-based picking. Picked particles were further 2D- and 3D-classified in iterative classification and selection rounds. Selected 3D classes were used for 3D template generation and refinement. The reported resolutions of the final maps were estimated using the gold-standard Fourier shell correlation. The local resolution predictions were made in RELION, with the resolution range estimated using the FSC = 0.143 criterion71. EM density visualization was done in UCSF Chimera75 and UCSF ChimeraX76.
As a representative of the image processing workflow, data for TRPV6RR were processed as follows. Initially, 997,105 particles were automatically picked using Laplacian-of-Gaussian approach in RELION. These particles were subjected to several rounds of 2D classification and 2D-class selection to generate 2D templates for template-based picking in RELION. After template-based picking, 584,265 particle images (2x-binned, 0.8 Å/pix) were imported into cryoSPARC for several iterative rounds of 2D classification, 2D-class selection, and 3D classification (heterogeneous refinement) to eliminate distorted particles. Particle images from the final 3D classification (329,613 particles) were imported back to RELION and subjected to 3D auto-refinement with C1-symmetry without masking, and with the previously published low-pass filtered to 40 Å map of hTRPV6-CtD in the apo-state (EMD-22622) as a reference. The resulting 3.53-Å reconstruction was used as a template for further processing. The refined particles were iteratively CTF-corrected for anisotropic magnification, beamtilt, and defocus, subjected to Bayesian polishing, and refined to 2.97-Å resolution with C1-symmetry and a soft mask. Particles from the last step were subjected to 3D classification into 10 classes using a soft mask and a 2.97-Å map as a reference. Particles from the two highest resolution classes (154,831) were refined to 2.43-Å resolution with C4 symmetry and a soft mask. The map from the final 3D reconstruction was automatically sharpened using RELION postprocessing and further used for model building.
To build TRPV6 models in Coot77, we used the previously published crystal or cryo-EM structures of TRPV6 as guides49. The models were tested for overfitting by shifting their coordinates by 0.5 Å (using shake) in Phenix78, refining each shaken model against a corresponding unfiltered half map, and generating densities from the resulting models in Chimera. Structures were visualized and figures were prepared in UCSF Chimera, UCSF ChimeraX, and Pymol79. The pore radius was calculated using HOLE80.
Fura-2 AM measurements
Wild-type or mutant human TRPV6 was expressed in HEK cells as described above. 48–72 hours after transduction, cells were harvested by centrifugation at 550 × g for 5 mins. The cells were resuspended in prewarmed HEPES-buffered saline (HBS: 118 mM NaCl, 4.8 mM KCl, 1 mM MgCl2, 5 mM d-glucose, 10 mM HEPES pH 7.4) containing 5 µg/ml of Fura-2 AM (Life Technologies) and incubated at 37 °C for 45 minutes. The loaded cells were then centrifuged for 5 mins at 550 × g, resuspended again in prewarmed HBS, and incubated at 37 °C for 20–30 mins in the dark. The cells were subsequently pelleted and washed twice, then resuspended in HBS for experiments. The cells were kept on ice in the dark for a maximum of ~2 hours before fluorescence measurements, which were conducted using spectrofluorometer QuantaMasterTM 40 (Photon Technology International) at room temperature in a quartz cuvette under constant stirring. Intracellular Ca2+ was measured by taking the ratio of fluorescence measurements at two excitation wavelengths (340 and 380 nm) and one emission wavelength (510 nm). The excitation wavelength was switched at 200-ms intervals.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Information. Cryo-EM density maps have been deposited to the Electron Microscopy Data Bank (EMDB) under the accession codes EMD-24890 [https://www.emdataresource.org/EMD-24890] (TRPV6open, GDN/CHS), EMD-24891 [https://www.emdataresource.org/EMD-24891] (TRPV6open, cNW11), EMD-24892 [https://www.emdataresource.org/EMD-24892] (TRPV6RR, cNW11), and EMD-24893 [https://www.emdataresource.org/EMD-24893] (TRPV6Eco, cNW11) (see Supplementary Table 1). The corresponding model coordinates have been deposited to the PDB under accession numbers 7S88 [https://doi.org/10.2210/pdb7S88/pdb] (TRPV6open, GDN/CHS), 7S89 [https://doi.org/10.2210/pdb7S89/pdb] (TRPV6open, cNW11), 7S8B [https://doi.org/10.2210/pdb7S8B/pdb] (TRPV6RR, cNW11), and 7S8C [https://doi.org/10.2210/pdb7S8C/pdb] (TRPV6Eco, cNW11) (see Supplementary Table 1). All other data are available from the corresponding author upon request. Source data are provided with this paper.
Yelshanskaya, M. V., Nadezhdin, K. D., Kurnikova, M. G. & Sobolevsky, A. I. Structure and function of the calcium-selective TRP channel TRPV6. J. Physiol. 599, 2673–2697 (2021).
Hoenderop, J. G. et al. Function and expression of the epithelial Ca2+ channel family: comparison of mammalian ECaC1 and 2. J. Physiol. 537, 747–761 (2001).
Vennekens, R., Voets, T., Bindels, R. J., Droogmans, G. & Nilius, B. Current understanding of mammalian TRP homologues. Cell Calcium 31, 253–264 (2002).
den Dekker, E., Hoenderop, J. G., Nilius, B. & Bindels, R. J. The epithelial calcium channels, TRPV5 & TRPV6: from identification towards regulation. Cell Calcium 33, 497–507 (2003).
Fecher-Trost, C., Weissgerber, P. & Wissenbach, U. TRPV6 channels. Handb. Exp. Pharm. 222, 359–384 (2014).
Zakharian, E., Cao, C. & Rohacs, T. Intracellular ATP supports TRPV6 activity via lipid kinases and the generation of PtdIns(4,5)P2. FASEB J. 25, 3915–3928 (2011).
Cai, R. et al. Autoinhibition of TRPV6 channel and regulation by PIP2. iScience 23, 101444 (2020).
Thyagarajan, B., Lukacs, V. & Rohacs, T. Hydrolysis of phosphatidylinositol 4,5-bisphosphate mediates calcium-induced inactivation of TRPV6 channels. J. Biol. Chem. 283, 14980–14987 (2008).
Voets, T., Janssens, A., Prenen, J., Droogmans, G. & Nilius, B. Mg2+-dependent gating and strong inward rectification of the cation channel TRPV6. J. Gen. Physiol. 121, 245–260 (2003).
Owsianik, G., Talavera, K., Voets, T. & Nilius, B. Permeation and selectivity of TRP channels. Annu. Rev. Physiol. 68, 685–717 (2006).
Lambers, T. T., Weidema, A. F., Nilius, B., Hoenderop, J. G. & Bindels, R. J. Regulation of the mouse epithelial Ca2+ channel TRPV6 by the Ca2+-sensor calmodulin. J. Biol. Chem. 279, 28855–28861 (2004).
Derler, I. et al. Dynamic but not constitutive association of calmodulin with rat TRPV6 channels enables fine tuning of Ca2+-dependent inactivation. J. Physiol. 577, 31–44 (2006).
Nilius, B. et al. Fast and slow inactivation kinetics of the Ca2+ channels ECaC1 and ECaC2 (TRPV5 and TRPV6). Role of the intracellular loop located between transmembrane segments 2 and 3. J. Biol. Chem. 277, 30852–30858 (2002).
Singh, A. K., McGoldrick, L. L., Twomey, E. C. & Sobolevsky, A. I. Mechanism of calmodulin inactivation of the calcium-selective TRP channel TRPV6. Sci. Adv. 4, eaau6088 (2018).
Bianco, S. D. et al. Marked disturbance of calcium homeostasis in mice with targeted disruption of the Trpv6 calcium channel gene. J. Bone Miner. Res. 22, 274–285 (2007).
Suzuki, Y. et al. Calcium channel TRPV6 is involved in murine maternal-fetal calcium transport. J. Bone Miner. Res 23, 1249–1256 (2008).
Weissgerber, P. et al. Male fertility depends on Ca2+ absorption by TRPV6 in epididymal epithelia. Sci. Signal. 4, ra27 (2011).
Weissgerber, P. et al. Excision of Trpv6 gene leads to severe defects in epididymal Ca2+ absorption and male fertility much like single D541A pore mutation. J. Biol. Chem. 287, 17930–17941 (2012).
Lieben, L. et al. Trpv6 mediates intestinal calcium absorption during calcium restriction and contributes to bone homeostasis. Bone 47, 301–308 (2010).
Wangemann, P. et al. Loss of cochlear HCO3- secretion causes deafness via endolymphatic acidification and inhibition of Ca2+ reabsorption in a Pendred syndrome mouse model. Am. J. Physiol. Ren. Physiol. 292, F1345–F1353 (2007).
Huybers, S. et al. Murine TNF(DeltaARE) Crohn’s disease model displays diminished expression of intestinal Ca2+ transporters. Inflamm. Bowel Dis. 14, 803–811 (2008).
Wu, G. et al. Suppression of intestinal calcium entry channel TRPV6 by OCRL, a lipid phosphatase associated with Lowe syndrome and Dent disease. Am. J. Physiol. Cell Physiol. 302, C1479–C1491 (2012).
Yang, S. S. et al. Generation and analysis of the thiazide-sensitive Na+-Cl- cotransporter (Ncc/Slc12a3) Ser707X knockin mouse as a model of Gitelman syndrome. Hum. Mutat. 31, 1304–1315 (2010).
Frick, K. K. et al. Increased biological response to 1,25(OH)2D3 in genetic hypercalciuric stone-forming rats. Am. J. Physiol. Ren. Physiol. 304, F718–F726 (2013).
Hache, S. et al. Alteration of calcium homeostasis in primary preeclamptic syncytiotrophoblasts: effect on calcium exchange in placenta. J. Cell Mol. Med. 15, 654–667 (2011).
Suzuki, Y. et al. TRPV6 variants interfere with maternal-fetal calcium transport through the placenta and cause transient neonatal hyperparathyroidism. Am. J. Hum. Genet 102, 1104–1114 (2018).
Burren, C. P. et al. TRPV6 compound heterozygous variants result in impaired placental calcium transport and severe undermineralization and dysplasia of the fetal skeleton. Am. J. Med. Genet. A 176, 1950–1955 (2018).
Nett, V., Erhardt, N., Wyatt, A. & Wissenbach, U. Human TRPV6-pathies caused by gene mutations. Biochim. Biophys. Acta Gen. Subj. 1865, 129873 (2021).
Masamune, A. et al. Variants that affect function of calcium channel TRPV6 are associated with early-onset chronic pancreatitis. Gastroenterology 158, 1626–1641 (2020).
Suzuki, Y. et al. Novel TRPV6 mutations in the spectrum of transient neonatal hyperparathyroidism. J. Physiol. Sci. 70, 33 (2020).
Yamashita, S., Mizumoto, H., Sawada, H., Suzuki, Y. & Hata, D. TRPV6 gene mutation in a dizygous twin with transient neonatal hyperparathyroidism. J. Endocr. Soc. 3, 602–606 (2019).
Zou, W. B. et al. TRPV6 variants confer susceptibility to chronic pancreatitis in the Chinese population. Hum. Mutat. 41, 1351–1357 (2020).
Stewart, J. M. TRPV6 as a target for cancer therapy. J. Cancer 11, 374–387 (2020).
Huber, S. M. Oncochannels. Cell Calcium 53, 241–255 (2013).
Peng, J. B., Suzuki, Y., Gyimesi, G. & Hediger, M. A. in Calcium Entry Channels in Non-Excitable Cells (eds J. A. Kozak & J. W. Putney, Jr.) 241–274 (2018).
Schwarz, E. C. et al. TRPV6 potentiates calcium-dependent cell proliferation. Cell Calcium 39, 163–173 (2006).
Lehen’kyi, V., Flourakis, M., Skryma, R. & Prevarskaya, N. TRPV6 channel controls prostate cancer cell proliferation via Ca2+/NFAT-dependent pathways. Oncogene 26, 7380–7385 (2007).
Bolanz, K. A., Hediger, M. A. & Landowski, C. P. The role of TRPV6 in breast carcinogenesis. Mol. Cancer Ther. 7, 271–279 (2008).
Bowen, C. V. et al. In vivo detection of human TRPV6-rich tumors with anti-cancer peptides derived from soricidin. PLoS ONE 8, e58866 (2013).
Haverstick, D. M., Heady, T. N., Macdonald, T. L. & Gray, L. S. Inhibition of human prostate cancer proliferation in vitro and in a mouse model by a compound synthesized to block Ca2+ entry. Cancer Res. 60, 1002–1008 (2000).
Landowski, C. P., Bolanz, K. A., Suzuki, Y. & Hediger, M. A. Chemical inhibitors of the calcium entry channel TRPV6. Pharm. Res. 28, 322–330 (2011).
Kovacs, G. et al. Inhibition of the human epithelial calcium channel TRPV6 by 2-aminoethoxydiphenyl borate (2-APB). Cell Calcium 52, 468–480 (2012).
Hofer, A. et al. Design, synthesis and pharmacological characterization of analogs of 2-aminoethyl diphenylborinate (2-APB), a known store-operated calcium channel blocker, for inhibition of TRPV6-mediated calcium transport. Bioorg. Med. Chem. 21, 3202–3213 (2013).
Simonin, C. et al. Optimization of TRPV6 calcium channel inhibitors using a 3D ligand-based virtual screening method. Angew. Chem. 54, 14748–14752 (2015).
Cunha, M. R. et al. Photoswitchable inhibitor of the calcium channel TRPV6. ACS Medicinal Chem. Lett. 10, 1341–1345 (2019).
Cunha, M. R. et al. Natural product inspired optimization of a selective TRPV6 calcium channel inhibitor. RSC Med. Chem. 11, 1032–1040 (2020).
Nilius, B. et al. Pharmacological modulation of monovalent cation currents through the epithelial Ca2+ channel ECaC1. Br. J. Pharm. 134, 453–462 (2001).
Singh, A. K., Saotome, K., McGoldrick, L. L. & Sobolevsky, A. I. Structural bases of TRP channel TRPV6 allosteric modulation by 2-APB. Nat. Commun. 9, 2465 (2018).
Bhardwaj, R. et al. Inactivation-mimicking block of the epithelial calcium channel TRPV6. Sci. Adv. 6, abe1508 (2020).
Clarke, M. J. Ruthenium metallopharmaceuticals. Coord. Chem. Rev. 232, 69–93 (2002).
Tapia, R. & Velasco, I. Ruthenium red as a tool to study calcium channels, neuronal death and the function of neural pathways. Neurochem. Int. 30, 137–147 (1997).
Yamada, K. Dual staining of some sulfated mucopolysaccharides with alcian blue (pH 1.0) and ruthenium red (pH 2.5). Histochemie 23, 13–20 (1970).
Pope, L., Lolicato, M. & Minor, D. L. Jr. Polynuclear ruthenium amines inhibit K2P channels via a “Finger in the Dam” mechanism. Cell Chem. Biol. 27, 511–524 (2020).
Choi, W., Clemente, N., Sun, W., Du, J. & Lu, W. The structures and gating mechanism of human calcium homeostasis modulator 2. Nature 576, 163–167 (2019).
McGoldrick, L. L. et al. Opening of the human epithelial calcium channel TRPV6. Nature 553, 233–237 (2018).
Saotome, K., Singh, A. K., Yelshanskaya, M. V. & Sobolevsky, A. I. Crystal structure of the epithelial calcium channel TRPV6. Nature 534, 506–511 (2016).
Singh, A. K., Saotome, K. & Sobolevsky, A. I. Swapping of transmembrane domains in the epithelial calcium channel TRPV6. Sci. Rep. 7, 10669 (2017).
Long, S. B., Campbell, E. B. & Mackinnon, R. Crystal structure of a mammalian voltage-dependent Shaker family K+ channel. Science 309, 897–903 (2005).
Sakipov, S., Sobolevsky, A. I. & Kurnikova, M. G. Ion permeation mechanism in epithelial calcium channel TRVP6. Sci. Rep. 8, 5715 (2018).
Cao, E., Liao, M., Cheng, Y. & Julius, D. TRPV1 structures in distinct conformations reveal activation mechanisms. Nature 504, 113–118 (2013).
Zubcevic, L., Le, S., Yang, H. & Lee, S. Y. Conformational plasticity in the selectivity filter of the TRPV2 ion channel. Nat. Struct. Mol. Biol. 25, 405–415 (2018).
Hughes, T. E. T. et al. Structural basis of TRPV5 channel inhibition by econazole revealed by cryo-EM. Nat. Struct. Mol. Biol. 25, 53–60 (2018).
Gao, Y., Cao, E., Julius, D. & Cheng, Y. TRPV1 structures in nanodiscs reveal mechanisms of ligand and lipid action. Nature 534, 347–351 (2016).
Nadezhdin, K. D. et al. Extracellular cap domain is an essential component of the TRPV1 gating mechanism. Nat. Commun. 12, 2154 (2021).
Conde, J. et al. Allosteric antagonist modulation of TRPV2 by piperlongumine impairs glioblastoma progression. ACS Cent. Sci. 7, 868–881 (2021).
van Goor, M. K. C., Hoenderop, J. G. J. & van der Wijst, J. TRP channels in calcium homeostasis: from hormonal control to structure-function relationship of TRPV5 and TRPV6. Biochim Biophys. Acta Mol. Cell Res. 1864, 883–893 (2017).
Na, T. & Peng, J. B. TRPV5: a Ca2+ channel for the fine-tuning of Ca2+ reabsorption. Handb. Exp. Pharm. 222, 321–357 (2014).
Goehring, A. et al. Screening and large-scale expression of membrane proteins in mammalian cells for structural studies. Nat. Protoc. 9, 2574–2585 (2014).
Nasr, M. L. et al. Covalently circularized nanodiscs for studying membrane proteins and viral entry. Nat. Methods 14, 49–52 (2017).
Suloway, C. et al. Automated molecular microscopy: the new Leginon system. J. Struct. Biol. 151, 41–60 (2005).
Scheres, S. H. RELION: implementation of a Bayesian approach to cryo-EM structure determination. J. Struct. Biol. 180, 519–530 (2012).
Punjani, A., Rubinstein, J. L., Fleet, D. J. & Brubaker, M. A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat. Methods 14, 290 (2017).
Zheng, S. Q. et al. MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat. Methods 14, 331–332 (2017).
Zhang, K. Gctf: real-time CTF determination and correction. J. Struct. Biol. 193, 1–12 (2016).
Pettersen, E. F. et al. UCSF Chimera–a visualization system for exploratory research and analysis. J. Comput Chem. 25, 1605–1612 (2004).
Pettersen, E. F. et al. UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci. 30, 70–82 (2021).
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr. D. Biol. Crystallogr. 66, 486–501 (2010).
Afonine, P. V. et al. Towards automated crystallographic structure refinement with phenix.refine. Acta Crystallogr. D. 68, 352–367 (2012).
The PyMOL Molecular Graphics System (DeLano Scientific, San Carlos, CA, USA, 2002).
Smart, O. S., Neduvelil, J. G., Wang, X., Wallace, B. A. & Samsom, M. S. HOLE: a program for the analysis of the pore dimensions of ion channel structural models. J. Mol. Graph. 14, 354–360 (1996).
We thank Harry Scott, Omar Davulcu (Pacific Northwest Center for Cryo-EM), Corey Hecksel (Stanford-SLAC), Ulrich Baxa, and Thomas Edwards (National Cancer Institute, The Frederick National Laboratory) for help with microscope operation and data collection. This research was, in part, supported by the National Cancer Institute’s National Cryo-EM Facility at the Frederick National Laboratory for Cancer Research under contract HSSN261200800001E. Some of this work was performed at the Stanford-SLAC Cryo-EM Center (S2C2) supported by the NIH Common Fund Transformative High-Resolution Cryo-Electron Microscopy program (U24 GM129541). A portion of this research was supported by NIH grant U24GM129547 and performed at the PNCC at OHSU and accessed through EMSL (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research. A.N. is a Walter Benjamin Fellow funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)–464295817. A.I.S. was supported by the NIH (R01 CA206573, R01 NS083660, R01 NS107253) and NSF (1818086).
The authors declare no competing interests.
Peer review information Nature Communications thanks Juan Du and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Neuberger, A., Nadezhdin, K.D. & Sobolevsky, A.I. Structural mechanisms of TRPV6 inhibition by ruthenium red and econazole. Nat Commun 12, 6284 (2021). https://doi.org/10.1038/s41467-021-26608-x