The selective exchange of ions across cellular membranes is a vital biological process. Ca2+-mediated signaling is implicated in a broad array of physiological processes in cells, while elevated intracellular concentrations of Ca2+ are cytotoxic. Due to the significance of this cation, strict Ca2+ concentration gradients are maintained across the plasma and organelle membranes. Therefore, Ca2+ signaling relies on permeation through selective ion channels that control the flux of Ca2+ ions. A key family of Ca2+-permeable membrane channels is the polymodal signal-detecting transient receptor potential (TRP) ion channels. TRP channels are activated by a wide variety of cues including temperature, small molecules, transmembrane voltage, and mechanical stimuli. While most members of this family permeate a broad range of cations non-selectively, TRPV5 and TRPV6 are unique due to their strong Ca2+ selectivity. Here, we address the question of how some members of the TRPV subfamily show a high degree of Ca2+ selectivity while others conduct a wider spectrum of cations. We present results from all-atom molecular dynamics simulations of ion permeation through two Ca2+-selective and two non-selective TRPV channels. Using a new method to quantify permeation cooperativity based on mutual information, we show that Ca2+-selective TRPV channel permeation occurs by a three-binding site knock-on mechanism, whereas a two-binding site knock-on mechanism is observed in non-selective TRPV channels. Each of the ion binding sites involved displayed greater affinity for Ca2+ over Na+. As such, our results suggest that coupling to an extra binding site in the Ca2+-selective TRPV channels underpins their increased selectivity for Ca2+ over Na+ ions. Furthermore, analysis of all available TRPV channel structures shows that the selectivity filter entrance region is wider for the non-selective TRPV channels, slightly destabilizing ion binding at this site, which is likely to underlie mechanistic decoupling.

The significance of Ca2+ in cellular function was first recognized by Sydney Ringer in 1883, who demonstrated that minute amounts of calcium were required for the contraction of cardiac muscle (Ringer, 1883). Ca2+ is now recognized as a versatile signaling agent, with cellular Ca2+ concentrations impacting a broad array of physiological processes ranging from cell proliferation to cell suicide (Carafoli, 2002; Berridge et al., 2003; Clapham, 2007; Patel, 2019). However, the cytoplasmic concentration of Ca2+ ions is usually kept low due to cytotoxic consequences (Bagur and Hajnóczky, 2017). Therefore, the controlled opening of channels in cellular and organellar membranes is one of the required mechanisms to allow the influx of this ion from the exoplasm and internal storage compartments into the cytoplasm; this subsequently initiates the Ca2+ signaling cascade. The question of how Ca2+ channels selectively permeate Ca2+ in low concentrations over vastly more abundant Na+ ions and yet conduct them at high rates has been a longstanding matter of fascination for ion channel researchers (Corry et al., 2001; Hille, 2001).

A key example of ion channels that mediate Ca2+ permeation across the cytoplasmic membrane is the transient receptor potential (TRP) channel superfamily. In their open state, these polymodal signal-detecting TRP channels allow the transmembrane flux of cations down their electrochemical gradient, thereby increasing the intracellular Ca2+ and Na+ concentration (Ramsey et al., 2006). The malfunction of TRP channels underlies a wide range of diseases, and they are therefore of immense biomedical importance, serving as drug targets for a variety of existing and candidate drugs (Moran, 2018).

TRP channels assemble primarily as homotetramers to form functional ion channels. A conserved structural feature across all TRP channels is the presence of six transmembrane helices (S1–S6) per subunit, forming two distinct transmembrane domains; a four-helix bundle comprising of helices S1–S4 forming the voltage-sensing like domain; and the pore-forming domain consisting of helices S5 and S6 (Hilton et al., 2019).

A four-residue ion selectivity filter (SF) is located at the entrance of the channel pore (Fig. 1). In addition to this conserved transmembrane architecture, members of the TRP superfamily display highly diverse extramembrane loops and N- and C-terminal domains between the different subfamilies (Van Goor et al., 2020).

TRP channels are gated open by a particularly wide range of stimuli, which include temperature, small molecules, transmembrane voltage changes, and mechanical cues (Caterina et al., 1997; Tominaga et al., 1998; Caterina et al., 2000). This superfamily of genes can be divided into seven main subfamilies: TRPA (ankyrin), TRPC (canonical), TRPM (melastatin), TRPML (mucolipin), TRPN (no mechanoreceptor potential C), TRPP (polycystin), and TRPV (vanilloid). It should be noted, however, that several less well-characterized TRP subfamilies have also been reported, including the TRPY (Palmer et al., 2001; Zhou et al., 2003; Chang et al., 2010), TRPVL (Peng et al., 2015), and TRPS (Himmel et al., 2020) subfamilies. The TRPV channels are the most intensely studied channel subfamily.

While they are all cation selective, most TRP channels electrophysiologically characterized to date show only limited discrimination between cation types, as well as between divalent and monovalent cations. However, the TRPV5 and TRPV6 channels are unique due to their high selectivity for Ca2+ cations over Na+ cations (PCa/PNa ∼100:1 from reversal potential measurements; Vennekens et al., 2000; Yue et al., 2001). Phylogenetic analysis has demonstrated that the TRPV5 and TRPV6 channels of vertebrates originated from an ancestral TRPV5/6 gene, which then diverged to form TRPV5 and TRPV6 from a duplication event after speciation (Flores-Aldama et al., 2000). Both of these channels are constitutively active due to basal levels of phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) in the cellular membrane and play a key role in Ca2+ homeostasis in the body (Van Goor et al., 2017). Despite their characteristic Ca2+ selectivity, both channels have been shown to permeate monovalent cations such as Na+ when divalent cations are absent (Nilius et al., 2000; Vennekens et al., 2000; Voets et al., 2003; Bödding and Flockerzi, 2004). By contrast, the remaining members of the TRPV subfamily, TRPV1–4, permeate both Ca2+ and Na+ cations, even in the presence of high Ca2+ concentrations, although they are still slightly Ca2+ selective, with a permeability ratio PCa/PNa ∼10:1. These channels gate in response to a number of stimuli, including raised temperature—in particular the archetypal member TRPV1 (Caterina et al., 1997; Tominaga et al., 1998; Caterina et al., 2000), which has led to TRPV1–4 being referred to as thermoTRPV channels—as well as endogenous and exogenous ligands.

In recent years, MD simulations have been successfully employed to shed light on ion channel function and the mode of action of channel-acting drugs in atomistic detail, for instance on K+ channels (Köpfer et al., 2014; Kopec et al., 2018), Na+ channels (Ulmschneider et al., 2013; Ke et al., 2014), Cl channels (McKiernan et al., 2020; Chavan et al., 2020), and ligand-gated ion channels (Sauguet et al., 2013). However, in the case of Ca2+-permeating channels, the conventional point-charge models used to describe uncoordinated Ca2+ ions have historically been inaccurate due to the neglected effects of electronic polarization (Kohagen et al., 2014a). This has resulted in an over-estimation of the binding energies between Ca2+ and proteins, hindering an accurate study of Ca2+ permeation in channels. Previous efforts to resolve this overestimation have included polarizable force field approaches (Li et al., 2015) and rescaled Ca2+ charges (Kohagen et al., 2014a; Kohagen et al., 2014b). More recently, Zhang et al. (2020) published a new multisite Ca2+ model specifically optimized for Ca2+–protein interactions (Zhang et al., 2020). This model correctly replicated experimental values for the hydration-free energy and the number of coordinated water molecules in the first solvation shell and showed Ca2+–protein binding energies comparable to the quantum mechanical and polarizable models. The multisite Ca2+ model has previously been used to investigate ion permeation of Ca2+ in a range of channels, including the type-1 ryanodine receptor (Zhang et al., 2020; Liu et al., 2021), α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (Schackert et al., 2022), and recently the E protein of SARS-CoV-2 (Antonides et al., 2022 Preprint).

In the present work, we set out to elucidate the molecular basis of Ca2+ selectivity and permeation in the TRPV channel subfamily. We conducted atomistic molecular dynamics (MD) simulations of TRPV channels under transmembrane voltage and compared the cation permeation mechanism observed in the Ca2+-selective TRPV5 and TRPV6 channels to the permeation mechanism in two exemplar non-selective TRPV channels, TRPV2 and TRPV3. In total, we observed 2,851 full ion traversals from 17.25 µs of MD simulations, allowing us to decipher the permeation mechanisms and principles of ion selectivity in the TRPV family with statistical power. Our findings suggest that ion conduction in TRPV channels proceeds via a cooperative knock-on mechanism involving multiple ion binding sites. The degree of cooperativity in ion permeation, linking the multiple binding sites, determines the degree of ion selectivity in the channels.

TRPV system construction

Truncated TRPV simulation systems consisting of the membrane-domain of the channels were constructed as described in Table 1. The systems were built using the CHARMM-GUI server (Jo et al., 2008). The charged N- and C-terminal residues were neutralized by capping with acetyl (ACE) and N-methylamide (CT3) groups, respectively, and all missing non-terminal residues were modeled (Jo et al., 2014). In the case of the TRPV5 system, the parameters for PI(4,5) P2 were generated using the CHARMM General Force Field (Vanommeslaeghe, 2010) through the ligand reader and modeler in CHARMM-GUI (Kim et al., 2017).

The structures were aligned in the membrane using the PPM server (Lomize et al., 2012), inserted into a 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine bilayer of 150 × 150 Å size using the CHARMM-GUI membrane builder (Jo et al., 2007; Wu et al., 2014), and then solvated. Ions were added using GROMACS 2020.2 (Abraham et al, 2015; Lindahl et al, 2020) to neutralize any system charges and add ions to a concentration of either 150 mM NaCl, 150 mM CaCl2, or a mixture of 75 mM NaCl and 75 mM CaCl2. In the case of simulations containing Ca2+, the standard CHARMM36m Ca2+ ions were then replaced with the multisite Ca2+ model of Zhang et al. (2020). Hydrogen mass repartitioning of the system was used to allow the use of 4 fs time steps in simulations of NaCl solutions. The multisite Ca2+ model used for simulations of CaCl2, however, is incompatible with a 4 fs time step, and therefore any simulations including Ca2+ cations were performed with hydrogen mass repartitioning but at a time step of 2 fs. The protein was restrained in the open-state by applying harmonic restraints on the α-carbon atoms of the lower gate residues (see Table 1).

MD simulations details

All simulations were performed using GROMACS 2020.2 (Abraham et al, 2015; Lindahl et al, 2020) and the CHARMM36m force field for the proteins, lipids, and ions (Huang et al., 2017). The TIP3P water model was used to model solvent molecules (Jorgensen et al., 1983). The system was minimized and equilibrated using the suggested equilibration inputs from CHARMM-GUI (Lee et al., 2016). In brief, the system was equilibrated using the isothermal-isobaric (NPT) ensemble for a total time of 1.85 ns with the force constraints on the system components being gradually released over six equilibration steps. The systems were then further equilibrated by performing a 15 ns simulation with no electric field applied. To prevent closing of the lower gate of the pore, harmonic restraints were applied to maintain the distance between the α-carbon atoms of the lower gate residues of each respective chain (Table 1). To drive ion permeation, an external electric field was applied by using the method of Aksimentiev and Schulten (2005) to production simulations with an E0 of −0.03 V nm−1; this resulted in a transmembrane voltage of ∼410 mV with negative polarity in the intracellular region. The temperature was maintained at 310 K using the Nosé-Hoover thermostat (Evans and Holian, 1985) and the pressure was maintained semi-isotropically at 1 bar using the Parrinello-Rahman barostat (Parrinello and Rahman, 1981). Periodic boundary conditions were used throughout the simulations. Long-range electrostatic interactions were modeled using the particle-mesh Ewald method (Darden et al., 1993) with a cutoff of 12 Å. The linear constraint solver algorithm (Hess et al., 1997) was used to constrain bonds with hydrogen atoms. All individual simulations were 250 ns long and repeated five times for each system, as summarized in Table 2. Additional details for all simulations of Ca2+-selective and non-selective TRPV channels are described in Tables S1 and S2, respectively, and details for additional control simulations are described in Table S3.

Simulation analysis

Analysis of MD trajectory data was performed using in-house written Python scripts, utilizing GROMACS modules (Abraham et al, 2015; Lindahl et al, 2020), the SciPy library of tools (Oliphant, 2007; Pérez and Granger, 2007; Millman and Aivazis, 2011; Van der Walt et al., 2011), and MDAnalysis (Michaud-Agrawal et al., 2011; Gowers et al., 2016). Analysis of the pore architecture was performed using the Channel Annotation Package (CHAP; Klesse et al, 2019). All plots were generated in Python using Matplotlib (Hunter, 2007) and Seaborn (Waskom et al., 2018).

Calculating conductance and selectivity from in silico electrophysiology experiments

The conductance of the channels (Gchannel) was calculated according to Eq. 1, where Np is the number of permeation events, Qion is the charge of the permeating ion in Coulomb, ttraj is the length of the trajectory, and Vm is the transmembrane voltage. The mean conductance and standard error were calculated from overlapping 50 ns windows of the trajectory.
Gchannel=Np×Qionttraj×Vm.
(1)

The selectivity (PCa/PNa) was calculated as the ratio between the total sum of Ca2+ permeation events and the total sum of Na+ permeation events from fivefold replicated 250 ns simulations in dicationic solutions with 75 mM CaCl2 and 75 mM NaCl.

Identification of cation binding sites from MD simulations of TRPV channels

Cation binding sites were identified by plotting timeseries of each permeating ion with respect to their position along the pore axis. To further validate this, the Featurizer function of PENSA (Vögele et al., 2021; Vögele et al., 2022) was used to identify the 12 most occupied ion binding sites, as determined by 3-D density maxima within 7 Å of the protein. This analysis was performed on a trajectory of concatenated fivefold replicated 250 ns simulations with a 200 ps time step from monocationic simulations.

Calculating ion occupancy probabilities and residence times in the identified cation binding sites

The ion occupancy (Oion) of the identified cation binding sites was calculated by dividing the number of frames (Noccupied), in which an ion’s center of geometry is within 3.5 Å of the center of geometry of the ion coordinating binding site atoms by the number of frames in the time window (Nframes). These atoms for the respective binding sites were: the carboxylate oxygen atoms of the α-position residue of the SF for binding site A, the carbonyl oxygen atoms of the β- and γ-position residues of the SF for binding site B, and the terminal carbon atoms of the hydrocarbon side chain of the isoleucine of the hydrophobic lower gate (I575 in TRPV5) and the amide oxygen atoms of neighboring asparagine residue (N572 in TRPV5). A cut-off distance of 3.5 Å was chosen based on the maximum reported distance for calcium–oxygen interactions (Nayal and Di Cera, 1994; Dudev and Lim, 2003; Deng et al., 2006). The mean ion occupancies and standard error were calculated from non-overlapping 50 ns windows of the fivefold replicated 250 ns simulation trajectories with a 20 ps time step.

The ion residence times (tr) were calculated by averaging the amount of time an individual ion was located within 3.5 Å of the center of geometry of the ion coordinating binding site atoms. The mean tr and standard error were calculated from fivefold replicated 250 ns simulation trajectories with a 20 ps time step.

Characterizing permeation cooperativity through mutual information using state-specific information (SSI) from PENSA

To characterize the level of cooperativity in the knock-on permeation mechanisms in TRPV channels, we used PENSA to calculate the SSI shared between discrete state transitions in the occupancy distributions of each binding site (Thomson et al., 2021,Preprint; Vögele et al., 2021; Vögele et al., 2022). A timeseries distribution with a time step of 20 ps for each binding site was obtained, whereby for each frame, if an ion occupied the binding site, then this ion’s atom ID number was recorded, whereas if the binding site was unoccupied, an ID of −1 was recorded. The ID numbers were discrete, and changes between ID numbers in each binding site, therefore, represent discrete state transitions. By quantifying the mutual information shared between changes to the ID numbers in each site, we were able to determine whether ion transitions at one site were coupled to transitions at another during a 20 ps time interval. From this, we concluded whether cations are “knocking” each other, or dissociation occurred independently from one another. The time interval was iteratively optimized to keep noise and finite sampling effects to a minimum (see below). We found that both were smallest when we used an interval of 20 ps.

Similar to McClendon et al. (2009), we observed that finite sampling resulted in independent distributions sharing mutual information (McClendon et al., 2009; Pethel and Hahs, 2014). To overcome this, we calculated a statistical threshold for each simulation via randomly permuted copies of the original data. Random permutations of the original data maintained marginal probabilities for binding site occupation in each simulation while at the same time quantifying the effect of finite sampling on the measurement of SSI. SSI was then calculated between two independently permuted versions of the occupancy distribution for the minimum entropy binding site. Since the upper bound of mutual information between two variables is equal to the lowest entropy of those variables, we used the binding site corresponding to the lowest entropy for obtaining the threshold. This ensured that the portion of SSI which could be attributed to random noise between any two binding sites was always less than or equal to the SSI. This measurement was repeated 1,000 times to resolve a Gaussian distribution from which we obtained the 99% confidence threshold. We subtracted this threshold from the measured values to resolve excess mutual information, or excess SSI (exSSI), shared in discrete state transitions. As it is not possible to transfer negative information, negative exSSI values were corrected to a value of 0.

We also derived a maximum SSI value representing a theoretical upper limit for the information that can be shared between two binding sites, where exSSImax is given by subtracting the random threshold from the minimum entropy of the two binding sites in question.
exSSI(A,B)max=min(H(A),H(B))threshold(A,B).
(2)
To quantify the interdependence of all three ion binding sites within the TRPV pores, the total correlation (TotCorr) was obtained using Eq. 3, where H(A), H(B), and H(C) represent the entropy of binding sites A, B, and C, respectively, and H(A, B, C) the joint entropy of binding sites A, B, and C.
TotCorr=H(A)+H(B)+H(C)H(A,B,C).
(3)

Characterizing the architecture of the SF of TRPV channels

To determine the area formed between residues in the SF, the area of the quadrilateral between the adjacent chains was calculated on the X and Y axes. For this, we used the carboxylate oxygen atoms of the adjacent chains for the α-position residue, and the carbonyl oxygen atoms for the β-, γ-, and δ-position residues. To quantify the SF areas from our MD simulations, the mean and standard error were calculated from non-overlapping 50 ns windows of the fivefold replicated 250 ns trajectories with a 200 ps time step. Furthermore, the SF areas of all TRPV structures deposited in the Protein Data Bank (PDB), available as of February 4, 2022, were determined. A total of 101 structures were analyzed, with non-tetrameric structures or structures without all the atoms of interest modeled not included. The mean and SEM were calculated for all the available structures for any particular TRPV channel.

Online supplemental material

Fig. S1 shows voltage-dependence of occupancy and ion selectivity. Fig. S2 shows additional cation binding sites in TRPV5. Fig. S3 shows pore radii and hydrophobicity in simulations of TRPV2, TRPV3, TRPV5, and TRPV6. Fig. S4 shows occupancy and residence times at postulated W583 ion binding site. Fig. S5 shows permeation traces at low ion concentration. Fig. S6 shows permeation traces in low voltage simulations. Fig. S7 shows impact of binding site affinity differences on exSSI. Fig. S8 shows total correlation of cation permeation between all binding sites in TRPV2, TRPV3, TRPV5, and TRPV6. Fig. S9 shows SF backbone dihedral distributions. Fig. S10 shows multiple sequence alignment of TRPV SF sequences. Fig. S11 shows effect of PI(4,5)P2 on pore radii. Tables S1, S2, and S3 show MD simulation details. Table S4 shows permeation times. Table S5 shows SSI details. Table S6 shows selectivity in dicationic solutions. Table S7 shows root-mean-square fluctuation (RMSF) of the SF backbone. Supplemental text at the end of the PDF provides additional information.

Continuous permeation of Ca2+ and Na+ in open-state TRPV5 and TRPV6 channels

We performed MD simulations of the pore domain of open-state TRPV5 (Hughes et al., 2018) and TRPV6 (McGoldrick et al., 2018) channels embedded in 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine lipid bilayers under transmembrane voltage (approximately −410 mV). The aqueous solutions contained either 150 mM CaCl2 or 150 mM NaCl (herein referred to as monocationic solutions) or a mixture consisting of 75 mM CaCl2 and 75 mM NaCl (herein referred to as dicationic solutions). All simulations performed with Ca2+ ions utilized the multisite Ca2+ model developed by Zhang et al. (2020) unless otherwise stated. In both the monocationic and the dicationic solutions, the applied voltage drove a continuous flow of permeating ions through all investigated open-state TRPV channels. Overall, we recorded 433 complete inward channel crossings for Ca2+ and 417 for Na+ in simulations of the Ca2+-selective TRPV channels.

In TRPV5 and TRPV6, Ca2+ ions traversed the entire pore within average time spans of 28.4 ± 3.9 ns (TRPV5) and 12.0 ± 1.0 ns (TRPV6; Table S4). The calculated Ca2+ and Na+ conductances from our simulations are shown in Table 3. The considerable Na+ conductances we observed agree with the experimental finding that the highly Ca2+-selective TRPV channels conduct Na+ well in the absence of Ca2+ (Nilius et al., 2000; Vennekens et al., 2000; Voets et al., 2003; Bödding and Flockerzi, 2004). Notably, these conductances are in quantitative agreement with published values measured for Na+ in vitro (Yue et al., 2001; Cha et al., 2007). By contrast, control simulations of TRPV5 using the default CHARMM36m force field parameters for Ca2+, but otherwise identical conditions, did not exhibit ion permeation; instead, the Ca2+ ions remained tightly bound to the protein ion binding sites for the entire course of the simulations. This observation is reflective of the shortcomings of standard parameters for divalent cations in fixed-point charge force fields and highlights the improved accuracy of multisite Ca2+ models in simulating divalent cation permeation and reproducing in vitro conductances. In addition, no Cl anions were observed to permeate TRPV channels in any of our simulations.

We note that the PCa/PNa values obtained from our simulations overall show lower Ca2+ selectivity than the reported literature values (Table S6). We surmised that this might be, at least partially, due to the higher voltages used in our simulations to enhance the sampling rate. Supplementary simulations performed at a lower voltage demonstrated that, indeed, the selectivity for Ca2+ increases with lower voltages across the membrane (Fig. S1). Below a voltage threshold of 205 mV, however, the sampling of permeation events in the simulations became very poor, such that we were not able to reliably probe the precise voltage range used in the experiments.

We can, additionally, not rule out the contribution of force field inaccuracies. Our simulations with Na+ show a remarkable agreement between experimentally recorded and simulated conductance. Even though the Ca2+ model we used has been carefully parameterized (Zhang et al., 2020), modeling divalent cations is a far from trivial task, and this is the first multisite Ca2+ model with which simulations of ion channel current have become possible. It can therefore not be excluded that further iterations of model refinement may eventually be required to not only reflect experimental solvation-free energies and protein affinity (Zhang et al., 2020) but also accurately reproduce experimental conductance values. This includes kinetic factors such as the correct characteristics of ion (electro-)diffusion in bulk solvent.

Pore cation binding sites and their preference for Ca2+ binding

Prior to the determination of the atomic structures of Ca2+ channels and the development of channel-permeable models for Ca2+ ions, it had been suggested from experimental observations that Ca2+ channels may obtain their selectivity through competition, i.e., by divalent cations, such as Ca2+, binding more tightly to their ion binding sites than monovalent cations, such as Na+ (Corry et al., 2001; Hille, 2001).

By analyzing the individual traces of permeating Ca2+ cations along the pore axis z of TRPV5 and TRPV6 over time, we identified three cation-binding sites inside the channels (Fig. 2). We refer to these cation binding sites as sites A, B, and C, viewed from the extracellular entrance of the channel SF to the hydrophobic lower gate. The three cation binding sites were further confirmed by 3-D density analysis using PENSA (Vögele et al., 2021; Vögele et al., 2022; Fig. S2). The PENSA analysis also identified a number of cation-binding sites outside of the pore within the extracellular loops of both TRPV5 and TRPV6 (Fig. S2), in line with the previous suggestion that TRPV6 contains negatively charged “recruitment sites” that funnel cations toward the entrance of the pore (Saotome et al., 2016; Sakipov et al., 2018).

Of the three binding sites we observed, binding site A is formed by the carboxylate oxygen atoms of the ring of acidic residues at the SF entrance (referred to here as the SF α-position); binding site B is formed by the carbonyl oxygen atoms of the bottom two SF residues (SF γ- and δ-positions); and binding site C is formed jointly by the hydrophobic gate consisting of a ring of isoleucine residues (I575 in TRPV5) and the amide oxygen atoms of the neighboring asparagine residues (N572 in TRPV5) near the cytoplasmic exit of the pore (Fig. 2). The location of these binding sites coincides with constrictions in the pore profile, as determined using CHAP (Klesse et al, 2019; Fig. S3; see Fig. S11). The distance between binding sites A and B is ∼5 Å, and that between binding sites B and C is ∼14 Å. We note that Hughes et al. (2018) reported a further constriction below the hydrophobic gate (binding site C) formed by W583 in the TRPV5 structure, and an analogous constriction at W583 can be observed in TRPV6. However, our simulations do not suggest that the side chains of W583 constitute a functionally important ion binding site, as shown in Fig. S4.

In monocationic Ca2+ solutions, the three binding sites showed Ca2+ occupancy probabilities of 0.69 ± 0.05, 0.67 ± 0.04, and 0.57 ± 0.03 in TRPV5, and 0.43 ± 0.04, 0.54 ± 0.05, and 0.29 ± 0.02 in TRPV6 (from A to C, respectively; Fig. 3). In monocationic Na + solutions, similar occupancies were observed. However, the Na+ residence times (tr) at the three binding sites were markedly lower than those observed for Ca2+, with ratios of tr (Ca2+):tr (Na+) varying between ∼35:1 and ∼3:1 (Fig. 3). These residence times suggest that Ca2+ ions have a greater affinity for these binding sites than Na+.

This observation was further substantiated when the occupancy of binding sites in dicationic solutions was analyzed, in which Ca2+ and Na+ cations are competing for the binding sites. Binding sites A, B, and C in the Ca2+-selective channels all showed high occupancies with Ca2+ in the mixed cationic solutions (Fig. 3). Across all the replicate simulations, we recorded average Ca2+ occupancies of 0.51 ± 0.05, 0.48 ± 0.06, and 0.37 ± 0.03 for binding sites A, B, and C in TRPV5, respectively, and of 0.68 ± 0.04, 0.61 ± 0.04, and 0.40 ± 0.04 for binding sites A, B, and C in TRPV6 (Fig. 3). By contrast, the Na+ occupancy of each of the three binding sites under these conditions was found to be below 0.07, both in TRPV5 and TRPV6; that is, the ratio between Ca2+ and Na+ occupancy varies between ∼85:1 and ∼7:1 (Fig. 3). These values indicate a free energy difference of between 11.5 and 5.0 kJ mol−1 for the preferential binding of Ca2+ over Na+.

A highly cooperative knock-on mechanism between three cation binding sites underpins selective Ca2+ permeation in TRPV channels

The observed increased affinity for Ca2+ cations at the pore binding sites compared to Na+ means that, in a mixed solution, Ca2+ will preferentially occupy these binding sites; however, this also implies that Ca2+ ions face a greater energy barrier when they dissociate from the binding sites. In monocationic solutions, this would result in a greatly reduced Ca2+ conductance with respect to Na+. For instance, based on our observed residence times in monocationic solutions, we would expect an ∼12-fold reduced Ca2+ unbinding rate compared to Na+ for binding site A in TRPV5. However, a much-reduced Ca2+ conductance is neither observed in our simulations nor the experimental literature. Due to the divalent charge of Ca2+, increasing the affinity to cation binding sites, this dichotomy had previously been suggested to exist, and it was hypothesized that this paradox could be resolved by assuming cooperativity between successive unbinding events such as in a knock-on mechanism (Corry et al., 2001; Hille, 2001).

In the classic knock-on mechanism, which for example underpins K+ channel function, ions transition into and out of multiple ion-binding sites in a highly correlated fashion (Armstrong, 1971; Hille and Schwarz, 1978; Neyton and Miller, 1988; Morais-Cabral et al., 2001; Köpfer et al., 2014). For example, early experiments by Hodgkin and Keynes and later flux-ratio measurements established that 3–3.4 K+ ions moved in lockstep with each other during permeation through K+ channels (Hodgkin and Keynes, 1955; Stampe and Begenisich, 1996).

For each permeating ion in a simulation of TRPV5, Fig. 4 shows the association and dissociation of Ca2+ and Na+ ions at binding sites A, B, and C from top to bottom as color code (bound to A, red; bound to B, orange; bound to C, yellow; transiting within the pore but not bound to a binding site, blue; located in extracellular solvent, dark gray; and located in intracellular solvent, light gray). As can be seen for Ca2+ in TRPV5 for example (Fig. 4 left), the plot demonstrates that (i) permeating Ca2+ ions spend the vast majority of their time within the pore at the three binding sites (reflected in the scarcity of blue boxes vs. red, orange, and yellow), (ii) dual and triple occupancy of the three sites, A, B, and C, with Ca2+ frequently observed (horizontal slices across plot: triple occupancy is observed in 27.2% of the simulation frames, dual occupancy in 49.7%), and (iii) transitions between states show a high degree of correlation, i.e., the ions frequently move in concert into and out of their respective binding sites (horizontal slices; binding state transitions). By contrast, during Na+ permeation (Fig. 4 center), the ions are predominantly transiting across the pore without occupying particular binding sites for extended time spans (blue, on average 53% of the traversal time for each ion).

To go beyond visual inspection of the trajectories and to assess the cooperativity of ion permeation in a quantitative way, we developed a new approach based on mutual information, taking into account the “state” of each ion binding site. To achieve this, we assigned a specific binding state (unoccupied or occupied with a specific ion) to binding sites A, B, and C and used our recently developed approach, SSI (Thomson et al., 2021,Preprint), on pairs of adjacent sites to quantify the degree of coupling between ion binding transitions at each of these sites (see Materials and methods). This analysis yields a coefficient quantifying the cooperativity between ion binding and unbinding at neighboring or more distant binding sites, where a greater coefficient signifies a higher degree of coupling. This coupling suggests that when an ion transitions from one site it is more likely that there is a transition at the other. To correct for the non-zero mutual information that samples of two completely independent variables can display due to finite-size effects, we followed the approach of McClendon et al. (2009) and Pethel and Hahs (2014) to yield excess mutual information, or exSSI. We also determined a theoretical upper limit for the maximum mutual information that can be shared between two binding sites by using the minimum state entropy among the two sites. Note that this quantity represents an absolute upper limit; reaching it would require both binding sites to exhibit idealized simultaneous states and state transitions throughout the entire simulated time.

The SSI analysis showed that in the Ca2+-selective TRPV channels under the simulated conditions, TRPV5 and TRPV6, a high level of information above noise is shared between the transition of ions into and out of binding sites A and B, respectively, both for the permeation of Ca2+ and Na+ (exSSI between 0.8 and 1.6 bits; Table S5 and Fig. 5). This suggests that the ion binding and unbinding processes at each of these binding sites are coupled to one another, constituting a knock-on mechanism at relatively short range. We observed three to four water molecules on average between cations bound at binding sites A and B during knock-on, demonstrating a “soft” knock-on mechanism to be in place, as opposed to the “direct” knock-on mechanism between dehydrated K+ ions in K+-selective cation channels (Köpfer et al., 2014). As detailed further below, our simulations indicate that only a moderate level of ion desolvation occurs in the SF of the studied TRPV channels, such that the hydration shell of the permeating ions remains largely intact during knock-on.

Similarly, the transitions of ions into and out of binding sites B and C show a large degree of correlation for both Ca2+ and Na+ (Fig. 5). In the case of binding sites B and C, however, this requires a remote knock-on mechanism to be in operation, since these sites are ∼14 Å apart. The concept of a remote knock-on event was first proposed by Tindjong et al. (2013) based upon Brownian dynamics simulations and observed by Zhang et al. (2020) in atomistic MD simulations of Ca2+ permeation in the RyR1 channel. Our SSI analysis suggests that the degree of cooperativity in the remote knock-on mechanism between binding sites B and C (Fig. 5 right) is slightly smaller than the cooperativity in the closer knock-on mechanism between binding sites A and B (Fig. 5 left). Simulations conducted at a lower CaCl2 concentration (Fig. S5) and reduced transmembrane voltage (Fig. S6) confirmed that the knock-on mechanism observed occurs also under these milder conditions; however, we did not achieve enough sampling to recalculate the exSSI values at these conditions.

Cation permeation in non-selective TRPV channels shows a lower degree of cooperativity

To determine if the remaining, non-selective TRPV channels showed a different permeation mechanism, we next performed simulations of the open-state TRPV2 (Dosey et al., 2019) and TRPV3 (Singh et al., 2019) channels using the same simulation approach as described for the Ca2+-selective TRPV channels. These simulations of non-selective TRPV channels also showed continuous ion permeation, with cation conductances, again, in good agreement with published conductance values measured in vitro (Table 4). Overall, we recorded 706 complete inward channel crossings for Ca2+ and 1,176 for Na+ from simulations of the non-selective TRPV channels.

The occupancy of binding sites B and C in the TRPV2 and TRPV3 systems showed no clear difference to the Ca2+-selective channels. By contrast, the occupancy of binding site A was reduced by ∼50% for both Na+ and Ca2+ ions in the monocationic solutions (Fig. 3). This suggests that cations are less well coordinated at binding site A in the non-selective TRPV channels, leading to lower affinity binding in the monocationic solutions. All binding sites, however, exhibited a preference for binding Ca2+ in the dicationic solutions. Note that, in the non-selective channels TRPV2 and TRPV3, this is coupled with a particularly low residence time for Ca2+ at binding site A, again suggesting higher exchange rates and weaker binding, despite its occupancy. That is, many Ca2+ ions are observed to diffuse back from A to bulk solution in TRPV2 and TRPV3.

We were therefore curious if the three-site knock-on mechanism described previously for Ca2+-selective TRPV channels is also at play in the non-selective TRPV channels. Our SSI analysis confirmed that the cooperativity between binding sites B and C in the non-selective TRPV channels was comparable with those calculated for the Ca2+-selective TRPV channels (Table S5 and Fig. 5). However, the correlation between ion binding transitions at binding sites A and B was substantially reduced in both of the non-selective TRPV systems (Fig. 5). The nearly complete absence of cooperativity from binding sites A and B demonstrates that a knock-on mechanism is not occurring between these two sites in the non-selective TRPV channels. Instead, our findings suggest that cation permeation in the non-selective TRPV channels occurs via a two-site knock-on mechanism between binding sites B and C.

Since the ion occupancy observed at binding site A is reduced in the case of the non-selective TRPV channels, it is plausible that this lower affinity also impacts the coupling between transitions at binding sites A and B. To test this notion further, the relationship between affinity differences of a pair of binding sites and the knock-on co-operativity was tested systematically by using a toy model with two energy wells (binding sites) possessing a range of different depths (affinities). As shown in Fig. S7, there is a linear relationship between the affinity difference and the observed exSSI. This demonstrates that the diminished affinity of binding site A is likely to be the major reason for the loss of cooperativity in the SF of the non-selective TRPV channels. Our results show that similar binding affinity is a necessary but not sufficient condition for a high degree of cooperativity between two cation binding sites.

Based on our SSI approach to quantify mutual information in binding and unbinding events at different ion binding sites and using the concept of total correlation to evaluate the overall cooperativity in a system across all coupled events, we next calculated the total correlation of ion permeation for all the TRP channels investigated. The reduction in the number of correlated knock-on sites within the non-selective TRP channels can be distinctly seen when comparing the total correlation for each of the non-selective and Ca2+-selective TRPV channels (Fig. S8).

The ion binding sites in non-selective TRPV channels preferentially bind Ca2+ over Na+, as in Ca2+-selective TRPV channels, which explains why the non-selective TRPV channels in fact show slight Ca2+ selectivity (PCa/PNa ∼10:1; Owsianik et al., 2006). However, our data suggest that the reduced level of coordination at binding site A, and especially the effect on the three-site cooperativity it imparts together with sites B and C, reduces the Ca2+ selectivity from PCa/PNa ∼100:1 seen in TRPV5 and TRPV6, and in this way ultimately determines the difference between Ca2+-selective and non-selective permeation.

Structural features distinguishing Ca2+-selective from non-selective permeation

To understand why cation coordination at binding site A is weakened in non-selective TRPV channels, decoupling its cooperativity, we investigated the area formed between the four subunits of the TRPV channel for each residue in the selectivity filter. The cross-sectional area formed by the carboxylate oxygen atoms at the SF α-position is larger in the non-selective TRPV channels than in the Ca2+-selective TRPV channels (Fig. 6). Interestingly, despite the increased average area formed by the carboxylate oxygen atoms at the α-position residue in the selectivity filter, the average area formed by the carbonyl oxygen atoms at both the γ- and δ-positions are smaller in non-selective TRPV channels than in Ca2+-selective TRPV channels (Fig. 6). These differences in selectivity architecture were further confirmed using the pore profile calculated using CHAP (Klesse et al, 2019; Fig. S3). Our simulations showed no differences in the flexibility of selectivity filters between Ca2+-selective and non-selective TRPV channels as determined by RMSF calculations of the backbone atoms (Table S7); we did however observe a small difference in the backbone dihedral angle distribution of the β-position residue, which we predict to be due to differences in stabilizing hydrophobic interactions of the residue side-chain (Fig. S9, S10, and S11).

To expand the geometric analysis to all available TRPV channel structures, we also calculated the average area formed by selectivity filter residues from all available TRPV structures deposited in the PDB (Fig. 6). Analysis of these static structures showed the same two trends observed within our MD simulations: (1) Ca2+-selective TRPV channels clearly have a smaller average area formed by the carboxylate oxygen atoms of the SF α-position; and (2) non-selective TRPV channels have a slightly narrower constriction at the SF γ- and δ-position residues. Our MD simulations suggested that the wider opening at the α-position leads to a weaker cation binding interaction at binding site A and cation coordination, which, in turn, decouples binding site A from the co-operative knock-on mechanism that underpins permeation in the Ca2+-selective TRPV channels.

The non-selective channels possess a narrower constriction formed by the side chains and carbonyl groups of the γ- and δ-position residues forming binding site B, whereas the Ca2+-selective TRPV channels exhibit an even narrower constriction at binding site A. Our simulations suggest that the greater occupancy of binding site A is a result of the more confined geometry of this charged binding site in the Ca2+-selective TRPV channels. This, in turn, leads to a higher occupancy of the uncharged binding site B in the Ca2+-selective TRPV channels (see Fig. 3), as binding site B receives cations from the adjacent binding site A via the knock-on mechanism, rather than having to bind them from bulk solution, as would be the case in the non-selective channels. Thus, despite the constriction not being as narrow, binding site B has a higher occupancy in the Ca2+-selective than in the non-selective TRPV channels.

Ca2+-selective permeation is not strongly linked to the solvation states of permeating cations

A previously reported mechanism of how cation selectivity can be achieved is by desolvation of permeating cations. Differences in desolvation energies between cationic species provide a thermodynamic penalty that can be more favorable for the permeation of one cationic species over another. Such a mechanism has been reported to underpin K+ selectivity over Na+ in K+ channels (Köpfer et al., 2014; Kopec et al., 2018), for example. Alternative mechanisms suggested to yield K+ selectivity are based on protein flexibility, especially the plasticity of the K+ channel SF (Noskov et al., 2004), the number of stacked ion binding sites in the SF (Derebe et al., 2011), where a reduction from four to three has been shown to abolish K+ selectivity, the kinetic model of selectivity (Thompson et al., 2009; Kim and Allen, 2011), and the coordination model (Bostick and Brooks, 2007; Varma et al., 2008). By contrast, Na+ selectivity has been suggested to rely chiefly on a “snug fit” coordination of the Na+ ion in the SFs of eukaryotic NaV channels and the preservation of its solvation shell while permeating bacterial NaV channels (Dudev and Lim, 2014).

To investigate whether a high degree of desolvation was a dominant factor in Ca2+ selectivity in TRPV channels, we determined the number of oxygen atoms within a 3 Å radius of the cations, representing their first solvation shell (Fig. 7). In the bulk solution, both Ca2+ and Na+ cations showed their expected water coordination number of 7 and 5.6, respectively. As the cations entered the pore, we saw a small degree of partial dehydration of permeating cations at the SF (Fig. 7). In particular, the carboxylate oxygen atoms of the acidic residue at the entrance of the SF coordinated an incoming cation, with these displacing up to approximately two coordinated water molecules from the first solvation shell of the cation. We also observe some desolvation around, or below, binding site C. However, no major differences in the solvation shell of permeating Ca2+ or Na+ cations, or indeed between the Ca2+-selective and non-selective TRPV channels were observed. The finding that the degree of desolvation does not differ substantially between Ca2+-selective and non-selective channels indicates that their selectivity is not exclusively based on a mechanism of ion dehydration. The largest level of desolvation is seen at ion binding site A in the TRPV5 channel for Ca2+. Since the overall desolvation penalty for Ca2+ is larger than for Na+ (Marcus, 1991), one would expect this to lead to the preferential binding of Na+ at this site, whereas the opposite is observed (Fig. 3), further arguing against dehydration as a major selectivity mechanism in TRPV channels.

Our simulations showed three main cation binding sites in the permeation pathway of TRPV channels, which we term binding sites A, B, and C. Binding sites A and B are formed by the carboxylate oxygen atoms of the α-position residue of the SF and by the carbonyl oxygen atoms of the γ- and δ-position residues of the SF, respectively. Binding site C is located just above the hydrophobic lower gate of the pore and is formed by the isoleucine residues of the lower gate and amide oxygen atoms of the neighboring asparagine residues.

The identification of these cation-binding sites in our MD simulations is in agreement with previously published TRPV structures and other MD simulations. In their crystal structure of the Rattus norvegicus TRPV6, Saotome et al. (2016) identified three regions of electron density within the channel pore which they interpreted as cation-binding sites. It should be noted that there are some minor differences in the residues forming binding site C to the binding site reported here, likely due to the structure of TRPV6 of Saotome et al. (2016) being in the closed-state, rather than the open-state structure of McGoldrick et al. (2018) used in the present study. During the transition from a closed state to the open state, TRP channels undergo a rotation about the pore-forming S6 helix, which changes the pore-facing residues (McGoldrick et al., 2018).

Moreover, numerous structures of TRPV channels deposited within the PDB are resolved with cations bound at one of the cation-binding sites identified in our MD simulations. These include structures from several orthologs of the TRPV channels simulated in this study, as well as of the TRPV1 and TRPV4 channels that were not simulated in this study. For example, structures of the closed-state TRPV1 channel of R. norvegicus (Kwon et al., 2021) and the open-state TRPV4 channel of Homo sapiens in complex with 4α-PDD (Botte et al., 2020 Preprint), both model cations bound at binding site B. This observation further validates the existence of the identified cation binding sites, as well as their conservation among TRPV channels and perhaps the wider TRP superfamily.

In addition to structural data, Sakipov et al. (2018) performed equilibrium MD simulations of Ca2+ movement in the closed-state structure of R. norvegicus TRPV6 and identified cation-binding sites in the SF (Sakipov et al., 2018). The binding sites from this work are also generally in agreement with our simulation results and the aforementioned structural data. However, Sakipov et al. (2018) reported that their simulation data identified two Ca2+ cations residing at binding site A, with adjacent chains each occupying an ion. In accordance with previous x-ray crystallographic data, two ions associated with binding site A are not seen in a major population of our simulation ensembles, which we attribute to the use of an optimized multisite model for Ca2+ ions in the present study (Zhang et al., 2020) and increased sampling for improved statistical analysis.

In monocationic solutions, we observed a high probability of at least two of the three binding sites (A, B, and C) being simultaneously occupied with either Na+ or Ca2+, respectively, with mostly insignificant differences between the occupancy values for Na+ and Ca2+ at each individual binding site. However, the residence times were markedly reduced at all sites for Na+ ions, and overall a tendency toward weaker binding for either ion at binding site A at the extracellular SF entrance was observed. In mixed, dicationic solutions of Na+ and Ca2+, by contrast, Ca2+ ions strongly outcompeted Na+ ions for association at all pore-binding sites, both in the Ca2+-selective and non-selective TRPV channels. Therefore, according to our data, all the binding sites display a much greater affinity for Ca2+. This gave rise to the question of how permeation efficiency for conducting Ca2+ ions is achieved in these channels and how this relates to their varying degrees of selectivity since higher affinity binding is usually expected to lead to slower permeation rates.

By ensuring cooperativity between binding/unbinding events at multiple binding sites, permeation rates can be enhanced (Corry et al., 2001; Hille, 2001). We hypothesized that the level of cooperativity between sites A, B, and C in the pore could underpin the difference between highly and less Ca2+-selective TRPV channels. We, therefore, developed a novel method to quantify co-operativity during ion permeation in pores with multiple ion binding sites based on mutual information and total correlation measures using the SSI approach (Thomson et al., 2021 Preprint). We anticipate that this method will be similarly useful for the study of permeation mechanisms and the basis of selectivity in other channels. Our analysis showed that there is a substantial degree of co-operativity for Ca2+ permeation between binding sites B and C across all TRPV channels, whereas a clear distinction exists between the co-operativity between binding sites A and B in Ca2+-selective and non-selective TRPV channels. In the non-selective TRPV channels, binding site A is decoupled from binding site B. By contrast, binding sites A and B are even more strongly coupled than binding sites B and C in the case of the Ca2+-selective channels. We suggest that this marked difference in co-operativity mechanistically explains the different levels of Ca2+ selectivity in TRPV channels.

Conclusion

We have characterized the cation permeation mechanisms in four members of the TRPV channel family. We identified three cation binding sites within the pore, each of which displayed greater affinity for Ca2+ binding over Na+. A novel application of mutual information between ion binding and unbinding at consecutive binding sites (SSI) enabled us to quantify the degree of knock-on taking place in ion permeation. This showed that the level of Ca2+ selectivity in TRPV channels is determined by the co-operativity, or coupling, between the transitions at three pore cation-binding sites. The Ca2+-selective TRPV channels display a highly correlated three-site knock-on, whereas the cation binding site at the extracellular entrance is decoupled from the mechanism in the non-selective TRPV channels, which reduces the overall preference for Ca2+ permeation.

MD simulation inputs and analysis scripts used for this study are deposited in a public GitHub repository, available at: https://github.com/cmives/Ca_selectivity_mechanism_of_TRPV_channels.

Crina M. Nimigean served as editor.

We thank Professor Chen Song of Peking University for providing parameters and discussion on the multisite Ca2+ model. We also thank the University of Dundee I.T. services for maintenance of the School of Life Sciences high-performance computing cluster that is utilized in this research.

C.M. Ives was supported by the Medical Research Council (grant number MR/N013735/1), and N.J. Thomson was supported by the UK Research and Innovation Biotechnology and Biological Sciences Research Council (grant number BB/M010996/1). Open Access funding provided by the University of Dundee.

The authors declare no competing financial interests.

Authors’ contributions: C.M. Ives and U. Zachariae conceived the idea and designed the computational study, C.M. Ives conducted the simulations, C.M. Ives analyzed the simulation data, C.M. Ives and N.J. Thomson developed the use of SSI for knock-on permeation, N.J. Thomson analyzed the SSI data, U. Zachariae supervised the work, C.M. Ives and U. Zachariae wrote the manuscript with contributions from N.J. Thomson, and all authors edited and reviewed the manuscript.

Abraham
,
M.J.
,
T.
Murtola
,
R.
Schulz
,
S.
Páll
,
J.C.
Smith
,
B.
Hess
, and
E.
Lindahl
.
2015
.
GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers
.
SoftwareX
.
1
:
19
25
.
Aksimentiev
,
A.
, and
K.
Schulten
.
2005
.
Imaging alpha-hemolysin with molecular dynamics: Ionic conductance, osmotic permeability, and the electrostatic potential map
.
Biophys. J.
88
:
3745
3761
.
Antonides
,
L.H.
,
Q.W.
Hurst
,
C.M.
Ives
,
K.
Ramberg
,
N.
Ostrovitsa
,
E.
Scanlan
,
M.
Caffrey
,
S.J.
Pitt
, and
U.
Zachariae
.
2022
.
The SARS-CoV-2 envelope (E) protein forms a calcium- and voltage-activated calcium channel
.
bioRxiv
.
(Preprint posted October 11, 2022)
.
Armstrong
,
C.M.
1971
.
Interaction of tetraethylammonium ion derivatives with the potassium channels of giant axons
.
J. Gen. Physiol.
58
:
413
437
.
Bagur
,
R.
, and
G.
Hajnóczky
.
2017
.
Intracellular Ca2+ sensing: Its role in calcium homeostasis and signaling
.
Mol. Cell
.
66
:
780
788
.
Berridge
,
M.J.
,
M.D.
Bootman
, and
H.L.
Roderick
.
2003
.
Calcium signalling: Dynamics, homeostasis and remodelling
.
Nat. Rev. Mol. Cell Biol.
4
:
517
529
.
Bödding
,
M.
, and
V.
Flockerzi
.
2004
.
Ca2+ dependence of the Ca2+-selective TRPV6 channel
.
J. Biol. Chem.
279
:
36546
36552
.
Bostick
,
D.L.
, and
C.L.
Brooks
III
.
2007
.
Selectivity in K+ channels is due to topological control of the permeant ion’s coordinated state
.
Proc. Natl. Acad. Sci. USA
.
104
:
9260
9265
.
Botte
,
M.
,
A.K.C.
Ulrich
,
R.
Adaixo
,
D.
Gnutt
,
A.
Brockmann
,
D.
Bucher
,
M.
Chami
,
N.
Bocquet
,
U.
Ebbinghaus-Kintscher
,
V.
Puetter
, et al
.
2020
.
Cryo-EM structural studies of the agonist complexed human TRPV4 ion-channel reveals novel structural rearrangements resulting in an open-conformation
.
bioRxiv
.
(Preprint posted October 13, 2020)
.
Carafoli
,
E.
2002
.
Calcium signaling: A tale for all seasons
.
Proc. Natl. Acad. Sci. USA
.
99
:
1115
1122
.
Caterina
,
M.J.
,
M.A.
Schumacher
,
M.
Tominaga
,
T.A.
Rosen
,
J.D.
Levine
, and
D.
Julius
.
1997
.
The capsaicin receptor: A heat-activated ion channel in the pain pathway
.
Nature
.
389
:
816
824
.
Caterina
,
M.J.
,
A.
Leffler
,
A.B.
Malmberg
,
W.J.
Martin
,
J.
Trafton
,
K.R.
Petersen-Zeitz
,
M.
Koltzenburg
,
A.I.
Basbaum
, and
D.
Julius
.
2000
.
Impaired nociception and pain sensation in mice lacking the capsaicin receptor
.
Science
.
288
:
306
313
.
Cha
,
S.K.
,
W.
Jabbar
,
J.
Xie
, and
C.L.
Huang
.
2007
.
Regulation of TRPV5 single-channel activity by intracellular pH
.
J. Membr. Biol.
220
:
79
85
.
Chang
,
Y.
,
G.
Schlenstedt
,
V.
Flockerzi
, and
A.
Beck
.
2010
.
Properties of the intracellular transient receptor potential (TRP) channel in yeast, Yvc1
.
FEBS Lett.
584
:
2028
2032
.
Chavan
,
T.S.
,
R.C.
Cheng
,
T.
Jiang
,
I.I.
Mathews
,
R.A.
Stein
,
A.
Koehl
,
H.S.
Mchaourab
,
E.
Tajkhorshid
, and
M.
Maduke
.
2020
.
A CLC-ec1 mutant reveals global conformational change and suggests a unifying mechanism for the CLC Cl/H+ transport cycle
.
Elife
.
9
:e53479.
Chung
,
M.-K.
,
H.
Lee
,
A.
Mizuno
,
M.
Suzuki
, and
M.J.
Caterina
.
2004
.
2-aminoethoxydiphenyl borate activates and sensitizes the heat-gated ion channel TRPV3
.
J. Neurosci.
24
:
5177
5182
.
Clapham
,
D.E.
2007
.
Calcium signaling
.
Cell
.
131
:
1047
1058
.
Corry
,
B.
,
T.W.
Allen
,
S.
Kuyucak
, and
S.-H.
Chung
.
2001
.
Mechanisms of permeation and selectivity in calcium channels
.
Biophys. J.
80
:
195
214
.
Darden
,
T.
,
D.
York
, and
L.
Pedersen
.
1993
.
Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems
.
J. Chem. Phys.
98
:
10089
10092
.
Deng
,
H.
,
G.
Chen
,
W.
Yang
, and
J.J.
Yang
.
2006
.
Predicting calcium-binding sites in proteins - a graph theory and geometry approach
.
Proteins
.
64
:
34
42
.
Derebe
,
M.G.
,
D.B.
Sauer
,
W.
Zeng
,
A.
Alam
,
N.
Shi
, and
Y.
Jiang
.
2011
.
Tuning the ion selectivity of tetrameric cation channels by changing the number of ion binding sites
.
Proc. Natl. Acad. Sci. USA
.
108
:
598
602
.
Dosey
,
T.L.
,
Z.
Wang
,
G.
Fan
,
Z.
Zhang
,
I.I.
Serysheva
,
W.
Chiu
, and
T.G.
Wensel
.
2019
.
Structures of TRPV2 in distinct conformations provide insight into role of the pore turret
.
Nat. Struct. Mol. Biol.
26
:
40
49
.
Dudev
,
T.
, and
C.
Lim
.
2014
.
Ion selectivity strategies of sodium channel selectivity filters
.
Acc. Chem. Res.
47
:
3580
3587
.
Evans
,
D.J.
, and
B.L.
Holian
.
1985
.
The Nose-Hoover thermostat
.
J. Chem. Phys.
83
:
4069
4074
.
Flores-Aldama
,
L.
,
M.W.
Vandewege
,
K.
Zavala
,
C.K.
Colenso
,
W.
Gonzalez
,
S.E.
Brauchi
, and
J.C.
Opazo
.
2020
.
Evolutionary analyses reveal independent origins of gene repertoires and structural motifs associated to fast inactivation in calcium-selective TRPV channels
.
Sci. Rep.
10
:
8684
.
Gowers
,
R.
,
M.
Linke
,
J.
Barnoud
,
T.
Reddy
,
M.
Melo
,
S.
Seyler
,
J.
Domański
,
D.
Dotson
,
S.
Buchoux
,
I.
Kenney
, and
O.
Beckstein
.
2016
.
MDAnalysis: A python package for the rapid analysis of molecular dynamics simulations
.
Proc.15th Python Sci. conf.
98
105
.
Hess
,
B.
,
H.
Bekker
,
H.J.
Berendsen
, and
J.G.
Fraaije
.
1997
.
LINCS: A linear constraint solver for molecular simulations
.
J. Comput. Chem.
18
:
1463
1472
.
Hille
,
B.
2001
.
Ion Channels of Excitable Membranes
. Third edition.
Sinauer
,
Sunderland, MA
.
Hille
,
B.
, and
W.
Schwarz
.
1978
.
Potassium channels as multi-ion single-file pores
.
J. Gen. Physiol.
72
:
409
442
.
Hilton
,
J.K.
,
M.
Kim
, and
W.D.
Van Horn
.
2019
.
Structural and evolutionary insights point to allosteric regulation of TRP ion channels
.
Acc. Chem. Res.
52
:
1643
1652
.
Himmel
,
N.J.
,
T.R.
Gray
, and
D.N.
Cox
.
2020
.
Phylogenetics identifies two eumetazoan TRPM clades and an eighth TRP family, TRP soromelastatin (TRPS)
.
Mol. Biol. Evol.
37
:
2034
2044
.
Hoenderop
,
J.G.J.
,
R.
Vennekens
,
D.
Müller
,
J.
Prenen
,
G.
Droogmans
,
R.J.M.
Bindels
, and
B.
Nilius
.
2001
.
Function and expression of the epithelial Ca2+ channel family: comparison of mammalian ECaC1 and 2
.
J. Physiol
.
537
:
747
761
.
Hodgkin
,
A.L.
, and
R.D.
Keynes
.
1955
.
The potassium permeability of a giant nerve fibre
.
J. Physiol.
128
:
61
88
.
Huang
,
J.
,
S.
Rauscher
,
G.
Nawrocki
,
T.
Ran
,
M.
Feig
,
B.L.
de Groot
,
H.
Grubmüller
, and
A.D.
MacKerell
Jr
.
2017
.
CHARMM36m: An improved force field for folded and intrinsically disordered proteins
.
Nat. Methods
.
14
:
71
73
.
Hughes
,
T.E.T.
,
R.A.
Pumroy
,
A.T.
Yazici
,
M.A.
Kasimova
,
E.C.
Fluck
,
K.W.
Huynh
,
A.
Samanta
,
S.K.
Molugu
,
Z.H.
Zhou
,
V.
Carnevale
, et al
.
2018
.
Structural insights on TRPV5 gating by endogenous modulators
.
Nat. Commun.
9
:
4198
4209
.
Hunter
,
J.D.
2007
.
Matplotlib: A 2D graphics environment
.
Comput. Sci. Eng.
9
:
99
104
.
Jo
,
S.
,
T.
Kim
, and
W.
Im
.
2007
.
Automated builder and database of protein/membrane complexes for molecular dynamics simulations
.
PLoS One
.
2
:e880.
Jo
,
S.
,
T.
Kim
,
V.G.
Iyer
, and
W.
Im
.
2008
.
CHARMM-GUI: A web-based graphical user interface for CHARMM
.
J. Comput. Chem.
29
:
1859
1865
.
Jo
,
S.
,
X.
Cheng
,
S.M.
Islam
,
L.
Huang
,
H.
Rui
,
A.
Zhu
,
H.S.
Lee
,
Y.
Qi
,
W.
Han
,
K.
Vanommeslaeghe
, et al
.
2014
.
CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues
.
Adv. Protein Chem. Struct. Biol.
96
:
235
265
.
Jorgensen
,
W.L.
,
J.
Chandrasekhar
,
J.D.
Madura
,
R.W.
Impey
, and
M.L.
Klein
.
1983
.
Comparison of simple potential functions for simulating liquid water
.
J. Chem. Phys.
79
:
926
935
.
Ke
,
S.
,
E.N.
Timin
, and
A.
Stary-Weinzinger
.
2014
.
Different inward and outward conduction mechanisms in NaVMs suggested by molecular dynamics simulations
.
PLOS Comput. Biol.
10
:e1003746.
Kim
,
I.
, and
T.W.
Allen
.
2011
.
On the selective ion binding hypothesis for potassium channels
.
Proc. Natl. Acad. Sci. USA
.
108
:
17963
17968
.
Kim
,
S.
,
J.
Lee
,
S.
Jo
,
C.L.
Brooks
III
,
H.S.
Lee
, and
W.
Im
.
2017
.
CHARMM-GUI ligand reader and modeler for CHARMM force field generation of small molecules
.
J. Comput. Chem.
38
:
1879
1886
.
Klesse
,
G.
,
S.
Rao
,
M.S.P.
Sansom
, and
S.J.
Tucker
.
2019
.
CHAP: A versatile tool for the structural and functional annotation of ion channel pores
.
J. Mol. Biol
.
431
:
3353
3365
.
Köpfer
,
D.A.
,
C.
Song
,
T.
Gruene
,
G.M.
Sheldrick
,
U.
Zachariae
, and
B.L.
de Groot
.
2014
.
Ion permeation in K⁺ channels occurs by direct Coulomb knock-on
.
Science
.
346
:
352
355
.
Kohagen
,
M.
,
M.
Lepšík
, and
P.
Jungwirth
.
2014a
.
Calcium binding to calmodulin by molecular dynamics with effective polarization
.
J. Phys. Chem. Lett.
5
:
3964
3969
.
Kohagen
,
M.
,
P.E.
Mason
, and
P.
Jungwirth
.
2014b
.
Accurate description of calcium solvation in concentrated aqueous solutions
.
J. Phys. Chem. B
.
118
:
7902
7909
.
Kopec
,
W.
,
D.A.
Köpfer
,
O.N.
Vickery
,
A.S.
Bondarenko
,
T.L.C.
Jansen
,
B.L.
de Groot
, and
U.
Zachariae
.
2018
.
Direct knock-on of desolvated ions governs strict ion selectivity in K+ channels
.
Nat. Chem.
10
:
813
820
.
Kwon
,
D.H.
,
F.
Zhang
,
Y.
Suo
,
J.
Bouvette
,
M.J.
Borgnia
, and
S.Y.
Lee
.
2021
.
Heat-dependent opening of TRPV1 in the presence of capsaicin
.
Nat. Struct. Mol. Biol.
28
:
554
563
.
Lee
,
J.
,
X.
Cheng
,
J.M.
Swails
,
M.S.
Yeom
,
P.K.
Eastman
,
J.A.
Lemkul
,
S.
Wei
,
J.
Buckner
,
J.C.
Jeong
,
Y.
Qi
, et al
.
2016
.
CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field
.
J. Chem. Theor. Comput.
12
:
405
413
.
Li
,
H.
,
V.
Ngo
,
M.C.
Da Silva
,
D.R.
Salahub
,
K.
Callahan
,
B.
Roux
, and
S.Y.
Noskov
.
2015
.
Representation of ion-protein interactions using the drude polarizable force-field
.
J. Phys. Chem. B
.
119
:
9401
9416
.
Lindahl
,
E.
,
M.J.
Abraham
,
B.
Hess
, and
D.
van der Spoel
.
2020
.
GROMACS 2020.2 manual
.
Liu
,
C.
,
A.
Zhang
,
N.
Yan
, and
C.
Song
.
2021
.
Atomistic details of charge/space competition in the Ca2+ selectivity of ryanodine receptors
.
J. Phys. Chem. Lett.
12
:
4286
4291
.
Lomize
,
M.A.
,
I.D.
Pogozheva
,
H.
Joo
,
H.I.
Mosberg
, and
A.L.
Lomize
.
2012
.
OPM database and PPM web server: Resources for positioning of proteins in membranes
.
Nucleic Acids Res.
40
:
D370
D376
.
Marcus
,
Y.
1991
.
Thermodynamics of solvation of ions. Part 5—Gibbs free energy of hydration at 298.15 K
.
J. Chem. Soc. Faraday Trans.
87
:
2995
2999
.
McClendon
,
C.L.
,
G.
Friedland
,
D.L.
Mobley
,
H.
Amirkhani
, and
M.P.
Jacobson
.
2009
.
Quantifying correlations between allosteric sites in thermodynamic ensembles
.
J. Chem. Theor. Comput.
5
:
2486
2502
.
McGoldrick
,
L.L.
,
A.K.
Singh
,
K.
Saotome
,
M.V.
Yelshanskaya
,
E.C.
Twomey
,
R.A.
Grassucci
, and
A.I.
Sobolevsky
.
2018
.
Opening of the human epithelial calcium channel TRPV6
.
Nature
.
553
:
233
237
.
McKiernan
,
K.A.
,
A.K.
Koster
,
M.
Maduke
, and
V.S.
Pande
.
2020
.
Dynamical model of the CLC-2 ion channel reveals conformational changes associated with selectivity-filter gating
.
PLOS Comput. Biol.
16
:e1007530.
Michaud-Agrawal
,
N.
,
E.J.
Denning
,
T.B.
Woolf
, and
O.
Beckstein
.
2011
.
MDAnalysis: A toolkit for the analysis of molecular dynamics simulations
.
J. Comput. Chem.
32
:
2319
2327
.
Millman
,
K.J.
, and
M.
Aivazis
.
2011
.
Python for scientists and engineers
.
Comput. Sci. Eng.
13
:
9
12
.
Morais-Cabral
,
J.H.
,
Y.
Zhou
, and
R.
MacKinnon
.
2001
.
Energetic optimization of ion conduction rate by the K+ selectivity filter
.
Nature
.
414
:
37
42
.
Moran
,
M.M.
2018
.
TRP channels as potential drug targets
.
Annu. Rev. Pharmacol. Toxicol.
58
:
309
330
.
Nayal
,
M.
, and
E.
Di Cera
.
1994
.
Predicting Ca2+-binding sites in proteins
.
Proc. Natl. Acad. Sci. USA
.
91
:
817
821
.
Neyton
,
J.
, and
C.
Miller
.
1988
.
Discrete Ba2+ block as a probe of ion occupancy and pore structure in the high-conductance Ca2+ -activated K+ channel
.
J. Gen. Physiol.
92
:
569
586
.
Nilius
,
B.
,
R.
Vennekens
,
J.
Prenen
,
J.G.J.
Hoenderop
,
R.J.M.
Bindels
, and
G.
Droogmans
.
2000
.
Whole-cell and single channel monovalent cation currents through the novel rabbit epithelial Ca2+ channel ECaC
.
J. Physiol.
527
:
239
248
.
Noskov
,
S.Y.
,
S.
Bernèche
, and
B.
Roux
.
2004
.
Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands
.
Nature
.
431
:
830
834
.
Oliphant
,
T.E.
2007
.
Python for scientific computing
.
Comput. Sci. Eng.
9
:
10
20
.
Owsianik
,
G.
,
K.
Talavera
,
T.
Voets
, and
B.
Nilius
.
2006
.
Permeation and selectivity of TRP channels
.
Annu. Rev. Physiol.
68
:
685
717
.
Palmer
,
C.P.
,
X.-L.
Zhou
,
J.
Lin
,
S.H.
Loukin
,
C.
Kung
, and
Y.
Saimi
.
2001
.
A TRP homolog in Saccharomyces cerevisiae forms an intracellular Ca2+-permeable channel in the yeast vacuolar membrane
.
Proc. Natl. Acad. Sci. USA
.
98
:
7801
7805
.
Parrinello
,
M.
, and
A.
Rahman
.
1981
.
Polymorphic transitions in single crystals: A new molecular dynamics method
.
J. Appl. Phys.
52
:
7182
7190
.
Patel
,
S.
2019
.
The secret life of calcium in cell signaling
.
Biochemist
.
41
:
34
37
.
Pérez
,
F.
, and
B.E.
Granger
.
2007
.
IPython: A system for interactive scientific computing
.
Comput. Sci. Eng.
9
:
21
29
.
Peng
,
G.
,
X.
Shi
, and
T.
Kadowaki
.
2015
.
Evolution of TRP channels inferred by their classification in diverse animal species
.
Mol. Phylogenet. Evol.
84
:
145
157
.
Pethel
,
S.D.
, and
D.W.
Hahs
.
2014
.
Exact test of independence using mutual information
.
Entropy
.
16
:
2839
2849
.
Ramsey
,
I.S.
,
M.
Delling
, and
D.E.
Clapham
.
2006
.
An introduction to TRP channels
.
Annu. Rev. Physiol.
68
:
619
647
.
Ringer
,
S.
1883
.
A further contribution regarding the influence of the different constituents of the blood on the contraction of the heart
.
J. Physiol.
4
:
29
42: 3
.
Sakipov
,
S.
,
A.I.
Sobolevsky
, and
M.G.
Kurnikova
.
2018
.
Ion permeation mechanism in epithelial calcium channel TRVP6
.
Sci. Rep.
8
:
5715
.
Saotome
,
K.
,
A.K.
Singh
,
M.V.
Yelshanskaya
, and
A.I.
Sobolevsky
.
2016
.
Crystal structure of the epithelial calcium channel TRPV6
.
Nature
.
534
:
506
511
.
Sauguet
,
L.
,
F.
Poitevin
,
S.
Murail
,
C.
Van Renterghem
,
G.
Moraga-Cid
,
L.
Malherbe
,
A.W.
Thompson
,
P.
Koehl
,
P.-J.
Corringer
,
M.
Baaden
, and
M.
Delarue
.
2013
.
Structural basis for ion permeation mechanism in pentameric ligand-gated ion channels
.
EMBO J.
32
:
728
741
.
Schackert
,
F.K.
,
J.
Biedermann
,
S.
Abdolvand
,
S.
Minniberger
,
C.
Song
,
A.J.R.
Plested
,
P.
Carloni
, and
H.
Sun
.
2022
.
Mechanism of calcium permeation in a glutamate receptor ion channel
.
ChemRxiv
.
Singh
,
A.K.
,
L.L.
McGoldrick
,
L.
Demirkhanyan
,
M.
Leslie
,
E.
Zakharian
, and
A.I.
Sobolevsky
.
2019
.
Structural basis of temperature sensation by the TRP channel TRPV3
.
Nat. Struct. Mol. Biol.
26
:
994
998
.
Stampe
,
P.
, and
T.
Begenisich
.
1996
.
Unidirectional K+ fluxes through recombinant Shaker potassium channels expressed in single Xenopus oocytes
.
J. Gen. Physiol.
107
:
449
457
.
Thompson
,
A.N.
,
I.
Kim
,
T.D.
Panosian
,
T.M.
Iverson
,
T.W.
Allen
, and
C.M.
Nimigean
.
2009
.
Mechanism of potassium-channel selectivity revealed by Na+ and Li+ binding sites within the KcsA pore
.
Nat. Struct. Mol. Biol
.
16
:
1317
1324
.
Thomson
,
N.J.
,
O.N.
Vickery
,
C.M.
Ives
, and
U.
Zachariae
.
2021
.
Ion-water coupling controls class A GPCR signal transduction pathways
.
bioRxiv
. (Preprint posted August 28, 2020).
Tindjong
,
R.
,
I.
Kaufman
,
D.G.
Luchinsky
,
P.V.E.
McClintock
,
I.A.
Khovanov
, and
R.S.
Eisenberg
.
2013
.
Stochastic dynamics of remote knock-on permeation in biological ion channels
.
22nd Int. Conf. Noise Fluctuations
:
1
4
.
Tominaga
,
M.
,
M.J.
Caterina
,
A.B.
Malmberg
,
T.A.
Rosen
,
H.
Gilbert
,
K.
Skinner
,
B.E.
Raumann
,
A.I.
Basbaum
, and
D.
Julius
.
1998
.
The cloned capsaicin receptor integrates multiple pain-producing stimuli
.
Neuron
.
21
:
531
543
.
Ulmschneider
,
M.B.
,
C.
Bagnéris
,
E.C.
McCusker
,
P.G.
Decaen
,
M.
Delling
,
D.E.
Clapham
,
J.P.
Ulmschneider
, and
B.A.
Wallace
.
2013
.
Molecular dynamics of ion transport through the open conformation of a bacterial voltage-gated sodium channel
.
Proc. Natl. Acad. Sci. USA
.
110
:
6364
6369
.
Van Der Walt
,
S.
,
S.C.
Colbert
, and
G.
Varoquaux
.
2011
.
The NumPy array: A structure for efficient numerical computation
.
Comput. Sci. Eng.
13
:
22
30
.
van Goor
,
M.K.C.
,
J.G.J.
Hoenderop
, and
J.
van der Wijst
.
2017
.
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
.
Van Goor
,
M.K.
,
L.
de Jager
,
Y.
Cheng
, and
J.
van der Wijst
.
2020
.
High-resolution structures of TRPV channels: Unveiling a functionally diverse group of ion channels
.
Protein Sci.
29
:
1569
1580
.
Vanommeslaeghe
,
K.
,
E.
Hatcher
,
C.
Acharya
,
S.
Kundu
,
S.
Zhong
,
J.
Shim
,
E.
Darian
,
O.
Guvench
,
P.
Lopes
,
I.
Vorobyov
, and
A.D.
Mackerell
Jr
.
2010
.
CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields
.
J. Comput. Chem.
31
:
671
690
Varma
,
S.
,
D.
Sabo
, and
S.B.K.
Rempe
.
2008
.
K+/Na+ selectivity in K channels and valinomycin: Over-coordination versus cavity-size constraints
.
J. Mol. Biol.
376
:
13
22
.
Vennekens
,
R.
,
J.G.
Hoenderop
,
J.
Prenen
,
M.
Stuiver
,
P.H.
Willems
,
G.
Droogmans
,
B.
Nilius
, and
R.J.
Bindels
.
2000
.
Permeation and gating properties of the novel epithelial Ca2+ channel
.
J. Biol. Chem
.
275
:
3963
3969
.
Vögele
,
M.
,
N.J.
Thomson
,
S.T.
Truong
, and
J.
McAvity
.
2021
. PENSA 0.2.7 version v0.2.7. https://github.com/drorlab/pensa
Vögele
,
M.
,
N.J.
Thomson
,
S.T.
Truong
,
J.
McAvity
,
U.
Zachariae
, and
R.O.
Dror
.
2022
.
Systematic analysis of biomolecular conformational ensembles with PENSA
.
arXiv
.
Voets
,
T.
,
A.
Janssens
,
J.
Prenen
,
G.
Droogmans
, and
B.
Nilius
.
2003
.
Mg2+-dependent gating and strong inward rectification of the cation channel TRPV6
.
J. Gen. Physiol.
121
:
245
260
.
Waskom
,
M.
,
O.
Botvinnik
,
D.
O’Kane
,
P.
Hobson
,
J.
Ostblom
,
S.
Lukauskas
,
D.C.
Gemperline
,
T.
Augspurger
,
Y.
Halchenko
,
J.B.
Cole
, et al
.
2018
.
mwaskom/seaborn: v0.9.0 (July 2018)
.
Watanabe
,
H.
,
J.
Vriens
,
A.
Janssens
,
R.
Wondergem
,
G.
Droogmans
, and
B.
Nilius
.
2003
.
Modulation of TRPV4 gating by intra- and extracellular Ca2+
.
Cell Calcium
.
33
:
489
495
.
Wu
,
E.L.
,
X.
Cheng
,
S.
Jo
,
H.
Rui
,
K.C.
Song
,
E.M.
Davila-Contreras
,
Y.
Qi
,
J.
Lee
,
V.
Monje-Galvan
,
R.M.
Venable
, et al
.
2014
.
CHARMM-GUI membrane builder toward realistic biological membrane simulations
.
J. Comput. Chem.
35
:
1997
2004
.
Yue
,
L.
,
J.B.
Peng
,
M.A.
Hediger
, and
D.E.
Clapham
.
2001
.
CaT1 manifests the pore properties of the calcium-release-activated calcium channel
.
Nature
.
410
:
705
709
.
Zhang
,
F.
,
S.M.
Hanson
,
A.
Jara-Oseguera
,
D.
Krepkiy
,
C.
Bae
,
L.V.
Pearce
,
P.M.
Blumberg
,
S.
Newstead
, and
K.J.
Swartz
.
2016
.
Engineering vanilloid-sensitivity into the rat TRPV2 channel
.
eLife
.
5
. e16409.
Zhang
,
A.
,
H.
Yu
,
C.
Liu
, and
C.
Song
.
2020
.
The Ca2+ permeation mechanism of the ryanodine receptor revealed by a multi-site ion model
.
Nat. Commun
.
11
:
922
.
Zhou
,
X.-L.
,
A.F.
Batiza
,
S.H.
Loukin
,
C.P.
Palmer
,
C.
Kung
, and
Y.
Saimi
.
2003
.
The transient receptor potential channel on the yeast vacuole is mechanosensitive
.
Proc. Natl. Acad. Sci. USA
.
100
:
7105
7110
.
This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).