The transient receptor potential channel vanilloid type 1 (TRPV1) is activated by a variety of endogenous and exogenous stimuli and is involved in nociception and body temperature regulation. Although the structure of TRPV1 has been experimentally determined in both the closed and open states, very little is known about its activation mechanism. In particular, the conformational changes that occur in the pore domain and result in ionic conduction have not yet been identified. Here we suggest a hypothetical molecular mechanism for TRPV1 activation, which involves rotation of a conserved asparagine in S6 from a position facing the S4–S5 linker toward the pore. This rotation is associated with hydration of the pore and dehydration of the four peripheral cavities located between each S6 and S4–S5 linker. In light of our hypothesis, we perform bioinformatics analyses of TRP and other evolutionary related ion channels, evaluate newly available structures, and reexamine previously reported water accessibility and mutagenesis experiments. These analyses provide several independent lines of evidence to support our hypothesis. Finally, we show that our proposed molecular mechanism is compatible with the prevailing theory that the selectivity filter acts as a secondary gate in TRPV1.
Transient receptor potential (TRP) channel vanilloid type 1 (TRPV1) is a member of the TRP family, which promotes nonselective cation current across cell membranes in response to heat, low pH, and inflammatory agents (Caterina et al., 1997; Tominaga et al., 1998; Hayes et al., 2000). Activation of this channel induces burning pain: the action potential triggered by TRPV1 in sensory ganglia and small sensory C and Aδ fibers is transmitted through the second-order projection neurons in the spinal cord to the thalamus and specific higher order brain areas, where it is perceived as pain (Julius and Basbaum, 2001; Carnevale and Rohacs, 2016).
TRPV1 is a homotetramer with each subunit composed of the transmembrane and cytoplasmic regions (Cao et al., 2013b; Liao et al., 2013; Gao et al., 2016). The transmembrane region has six helical segments (S1–S6); four of them (S1–S4) constitute a domain structurally homologous to the voltage sensor domain of voltage-gated ion channels. The remaining two helices (S5 and S6) assemble into the pore domain. Cryo-EM structures of TRPV1 in different functional states have revealed key elements of its activation mechanism (Cao et al., 2013b; Liao et al., 2013; Gao et al., 2016). For instance, it became apparent that the gate is located at the crossing of the S6 helical bundle. The gate is closed in the absence of an activating stimulus, while at room temperature it opens upon binding of the agonist capsaicin between S4 and S5 (Cao et al., 2013b). The molecular mechanism of this conformational change has been in part clarified by mutagenesis experiments (Yang et al., 2015; Gao et al., 2016): capsaicin forms a hydrogen bond with the S4–S5 linker and pulls it in the outward direction; as for voltage-gated ion channels, this motion relieves the constriction exerted by the linker on the S6 helical bundle.
However, the cryo-EM structure of the capsaicin-bound state is closed at cryogenic temperatures: the gate is narrower than that of the open state (Cao et al., 2013b), and in simulations it is dehydrated (Kasimova et al., 2018). Consistently, Jara-Oseguera et al. (2016) have shown that TRPV1 activation becomes increasingly less probable as the temperature is decreased even when capsaicin is present in large concentrations. Thus, the capsaicin-bound state is closed at cryogenic temperatures and opens as temperature is increased, with the midpoint of activation located approximately at 300°K.
We recently used molecular dynamics and free energy calculations to study the molecular mechanism of the closed-to-open transition in the capsaicin-bound state (Kasimova et al., 2018). Our major finding was that the conserved asparagine N676 in S6 adopts two alternative conformations in the open and closed states, projecting its side chain toward either the pore or four peripheral cavities located between the S6 helix and the S4–S5 linker, respectively. When TRPV1 is in the closed state, these cavities are connected to the intracellular solution and host several water molecules. Dehydration of the peripheral cavities correlates with the rotation of N676 from these cavities toward the pore. In this rotated conformation, N676 promotes pore hydration and thus permeability to ions (Fig. 1).
Here, we will briefly summarize our previous computational findings and discuss several independent pieces of experimental evidence that corroborate our hypothesis. We will show how these lend credibility to the rotation of N676 and the presence of the peripheral cavities. Finally, we will present an analysis of our previously reported simulations suggesting that the selectivity filter acts as a secondary gate.
Materials and methods
Molecular dynamics simulations
The molecular dynamics trajectories of the capsaicin-bound closed and open states were taken from our previous study (Kasimova et al., 2018). Briefly, the structure of the TRPV1 capsaicin-bound state (Cao et al., 2013b) with four capsaicin molecules (Elokely et al., 2016) was embedded in a hydrated 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer and surrounded by 150 mM NaCl solution. Four to six water molecules were added inside each peripheral cavity in one case (this system further relaxed to the closed state), while these cavities were left empty in the other case (this system further relaxed to the open state; Kasimova et al., 2018). The CHARMM36 force field (Mackerell et al., 2004) was used for the protein and the POPC lipids, the parameters derived in (Elokely et al., 2016) were used for capsaicin, and the TIP3P model (Jorgensen et al., 1983) was used for water molecules. The molecular dynamics simulations were performed following a multistep protocol (for the details, see Kasimova et al., 2018) using NAMD 2.10 (Phillips et al., 2005). The pressure and temperature were kept constant at 1 atm and 300°K, respectively, using Langevin dynamics. A cutoff of 11 Å with a switching function between 8 and 11 Å was used for the van der Waals interactions. A cutoff of 11 Å was considered for the short-range component of electrostatic interactions, and the long-range component was computed using particle mesh Ewald (Darden et al., 1993). The overall length of the molecular dynamics trajectory for each state was 750 ns. In addition, we performed the molecular dynamics simulations of the TRPV1 apo state (Liao et al., 2013) using an analogous setup. The overall length of the molecular dynamics trajectory for this state was 200 ns.
Analysis of sequence conservation
The multiple sequence alignments of the TRP and voltage-gated sodium and calcium channels were taken from our previous studies (Palovcak et al., 2015; Kasimova et al., 2016). To analyze the sequence conservation for these families, we computed the Shannon entropy profile using the logos website (Crooks et al., 2004).
Estimation of the accessibility of S6 residues to the intracellular solution
To estimate the accessibility of S6 residues to the intracellular solution, we considered the second halves of the molecular dynamics trajectories. We separated the intracellular solution from the rest applying the following procedure. We first calculated the three-dimensional map of water occupancy using the Volmap tool of VMD (Humphrey et al., 1996). Based on this, we generated a binary map by setting each voxel to either 0 or 1 depending on whether or not the local value of water occupancy overcomes a predefined threshold. After that, we used the gmx cluster of Gromacs 5.0.4 (Abraham et al., 2015) to group the voxels. This allowed us to identify the cluster containing the intracellular solution and all the regions connected to it. Finally, as a measure of accessibility, we calculated the overlap between the occupancy of the intracellular solution and that of each S6 residue.
Results and discussion
The π-bulge in the S6 helix is evolutionarily conserved across the TRP family
We start by recalling a remarkable feature of the capsaicin-bound structure (Cao et al., 2013b): the presence of a π-helical segment between Y671 and N676 (Fig. 2 A). α- and π-helices show an important difference: while in α-helices, hydrogen bonds join main chain groups that are four residues apart (i–i+4), this spacing changes to five residues in π-helices (i–i+5). Therefore, the presence of a π-helical segment in an α-helix implies a mismatch in the pattern of hydrogen bonds, which results in unpaired carbonyl and amino groups. In practice, one residue does not participate in the network of hydrogen bonds and bulges out of the α-helix, forming a so-called π-bulge (Fig. 2 A). The π-bulge is mobile: the bulging residue can “reabsorb” by forming hydrogen bonds at the expense of a neighboring residue, which, as a result, bulges out of the helix.
In one of our previous investigations, we hypothesized, on the basis of the sole TRPV1 structure, that this structural motif is conserved across the entire family of TRP channels (Palovcak et al., 2015). To test this hypothesis, we analyzed patterns of correlations in a large multiple sequence alignment containing ∼3,000 distinct sequences of genes encoding for TRP channels. The presence of a π-helical segment changes the register of the downstream α-helix: it offsets by one position all the subsequent residues, effectively imparting a rotation of approximately 100° to the downstream α-helix (Fig. 2 B). The change of the register results in a specific pattern of contacts between S6 and the neighboring S5. We detected this pattern in the multiple sequence alignment as strong correlations between S5 and S6 residues; importantly, these pairs of residues are in spatial proximity only with the “correct” register of S6, i.e., if the π-bulge is present (Palovcak et al., 2015). The conservation of this specific pattern of contacts between S5 and S6 across the entire TRP family suggests that the π-bulge is conserved as well.
Since then, several experimental structures of TRP channels have confirmed our hypothesis about the structural conservation of the π-bulge. The structures of TRPV1, TRPV3, TRPV5, TRPV6, TRPA1, TRPM2, TRPM4, TRPM7, TRPC4, TRPP2, TRPN-like channel NOMPC, TRPML1, and TRPML3 all show a π-helical segment in S6 (Fig. 2 C and Table 1; Cao et al., 2013b; Liao et al., 2013; Barad et al., 2015; Paulsen et al., 2015; Gao et al., 2016; Shen et al., 2016; Chen et al., 2017; Grieben et al., 2017; Guo et al., 2017; Hirschi et al., 2017; Jin et al., 2017; Schmiege et al., 2017; Winkler et al., 2017; Zhou et al., 2017; Autzen et al., 2018; Duan et al., 2018a,b; Hughes et al., 2018b; McGoldrick et al., 2018; Singh et al., 2018a,b; Su et al., 2018a; Vinayagam et al., 2018; Zhang et al., 2018). Interestingly, structures without the π-bulge have also been determined for some of these TRP channels (TRPV3, TRPV5, TRPV6, and TRPP2; Saotome et al., 2016; Singh et al., 2017, 2018a,b; Hughes et al., 2018a; Hulse et al., 2018; McGoldrick et al., 2018; Su et al., 2018b; Zheng et al., 2018), suggesting that a transition between the π- and α-helical conformations might be part of the functional cycle (Fig. 2 C and Table 1). Finally, only TRPV2, TRPV4, TRPM8, and TRPC3 do not show a π-helical segment in S6 (Huynh et al., 2016; Zubcevic et al., 2016, 2018; Deng et al., 2018; Fan et al., 2018; Hughes et al., 2018a; Yin et al., 2018). Either these TRP channels lack altogether the π-helical conformation or, alternatively, this conformation has not been experimentally captured yet.
The π-bulge causes conformational flexibility of neighboring residues
Analysis of a comprehensive multiple sequence alignment of TRP channels revealed that the pore-lining residues are all hydrophobic except for one that in most of the sequences is an asparagine (N676 in TRPV1, Fig. 3, A and B; Palovcak et al., 2015). This feature is strictly conserved in the two evolutionarily divergent families of TRPP and TRPML channels as well. In TRPML channels, instead of an asparagine, some sequences have a serine; also in these cases, serine is the only hydrophilic pore-lining residue. The conserved asparagine is located only ∼1.5 helical turns away from the π-bulge, which potentially endows this residue with conformational flexibility. Structures of different TRP channels (Cao et al., 2013b; Liao et al., 2013; Barad et al., 2015; Paulsen et al., 2015; Gao et al., 2016; Huynh et al., 2016; Saotome et al., 2016; Shen et al., 2016; Zubcevic et al., 2016, 2018; Chen et al., 2017; Grieben et al., 2017; Guo et al., 2017; Hirschi et al., 2017; Jin et al., 2017; Schmiege et al., 2017; Singh et al., 2017, 2018a,b; Winkler et al., 2017; Zhou et al., 2017; Autzen et al., 2018; Deng et al., 2018; Duan et al., 2018a,b; Fan et al., 2018; Hughes et al., 2018b,a; Hulse et al., 2018; McGoldrick et al., 2018; Su et al., 2018a,b; Vinayagam et al., 2018; Yin et al., 2018; Zhang et al., 2018; Zheng et al., 2018) indeed show that the orientation of N676 (or the corresponding amino acid) with respect to the pore can differ by as much as 180° between two structures (Fig. 3 C and Table 1). For instance, in the TRPV subfamily, the conserved asparagine points toward the pore in TRPV1 and some structures of TRPV3, TRPV5, and TRPV6 (Cao et al., 2013b; Liao et al., 2013; Barad et al., 2015; Gao et al., 2016; McGoldrick et al., 2018; Singh et al., 2018a,b ; Hughes et al., 2018b); points in the opposite direction (toward the S4–S5 linker) in TRPV2, TRPV4, and the remaining TRPV3, TRPV5, and TRPV6 structures (McGoldrick et al., 2018; Singh et al., 2018a,b ; Huynh et al., 2016; Saotome et al., 2016; Zubcevic et al., 2016, 2018; Singh et al., 2017; Deng et al., 2018; Hughes et al., 2018a); and, finally, lies just outside the pore in the TRPV1 capsaicin-bound structure (Cao et al., 2013b). Similar conformations are observed in other TRP families including TRPM, TRPC, and TRPP (Shen et al., 2016; Grieben et al., 2017; Guo et al., 2017; Winkler et al., 2017; Autzen et al., 2018; Duan et al., 2018b,a; Fan et al., 2018; Hulse et al., 2018; Su et al., 2018a,b; Vinayagam et al., 2018; Yin et al., 2018; Zhang et al., 2018; Zheng et al., 2018). Importantly, the structures with the conserved asparagine inside the pore or just outside of it all have the π-bulge; by contrast, those in which this residue points to the S4–S5 linker do not show any π-helical segment in S6 (Table 1).
This structural heterogeneity is in itself striking: evolutionary conservation points to functional relevance, yet the residue lacks a well defined conformation. Is the role of the conserved asparagine related to its flexibility? In our previous study (Kasimova et al., 2018), we have shown that in the TRPV1 capsaicin-bound structure, the rotation of N676 inside or outside the pore significantly affects the polarity of the latter and hence ionic transport: only when N676 is exposed to the pore, the latter is continuously hydrated and the free energy barrier for ionic transport at the level of the gate is ∼5.2 kcal/mol. By contrast, when N676 is interacting with the S4–S5 linker, the free energy barrier is larger than 12 kcal/mol, thereby preventing ionic transport. Consistent with these predictions, recently determined experimental structures of TRPV3, TRPV5, and TRPV6 show that the conserved asparagine faces the pore in the open conformations and the S4–S5 linker in the closed ones (Saotome et al., 2016; Singh et al., 2017, 2018a,b; Hughes et al., 2018a,b; McGoldrick et al., 2018).
It is important to note that the location of the conserved asparagine side chain inside the pore may not be the only determinant of channel opening. The radius of the pore should be, in any case, large enough to permit the passage of hydrated ions. Accordingly, several structures of TRP channels showing the conserved asparagine in the pore-facing configuration and the pore radius smaller than 2 Å are likely representatives of the closed states (Table 1).
Y671 motion affects gating at the selectivity filter and is coupled to N676 rotation
Are there functional data showing that the rotation of S6 underlies TRP channels’ gating? Steinberg et al. (2017) recently reported that TRPV1 activation correlates with a conformational change at the level of the π-bulge where a nonnatural amino acid bearing a coumarin side chain was introduced. Coumarin is a fluorophore with exquisite sensitivity to environmental polarity (Wang et al., 2006; Wagner, 2009); when genetically encoded into proteins, it can be used as reporter of changes in local hydration and conformational rearrangements. In the TRPV1 channel with Y671 substituted for coumarin, activation by capsaicin was found to cause an increase in both the photon counts and optical fluctuations. This effect suggests that coumarin undergoes a conformational transition from a more hydrated state to a less hydrated one (Steinberg et al., 2017).
While the exact conformational change cannot be inferred from these experiments, molecular dynamics simulations performed on the TRPV1 wild type suggested a possible explanation for the observed changes in fluorescence (Steinberg et al., 2017). The simulations showed that the surface area of Y671 accessible to solvent decreases from 30% to 15% on passing from the closed to open states, in agreement with the experimental data for coumarin. The reason for this change is the conformational rearrangement of Y671 coupled to gating through N676. In the TRPV1 open state, Y671 is oriented perpendicularly with respect to the pore axis, and it is shielded from water molecules by the pore helices. The N676 side chain is located inside the pore, where it forms a hydrogen bond with the I672 carbonyl group, which stabilizes the π-bulge on I672 (Fig. 4 A). In the closed state, the N676 side chain moves toward the S4–S5 linker and no longer donates a hydrogen bond to the main chain. As a consequence, the π-bulge is found preferentially on Y671 (Fig. 4 A). The displacement of the π-bulge from I672 to Y671 causes a change of Y671 orientation with respect to the pore: it goes from being perpendicular to the pore axis in the open state to being parallel to it in the closed state (Fig. 4, A and B). In the conformation parallel to the pore axis, Y671 becomes less shielded by the pore helices and hence more accessible to water molecules (Fig. 4 B).
Interestingly, the molecular dynamics simulations also suggest that the observed reorientation of Y671 is coupled to the conformational change of the selectivity filter. In particular, when the Y671 side chain is parallel to the pore axis, it donates a hydrogen bond to the T641 carbonyl group on the pore helix. This interaction displaces the pore helix toward the extracellular solution, resulting in a slight constriction of the selectivity filter. In our previous study, we have shown that this constriction is sufficient to increase the free energy barrier for ionic transport from ∼1.5 in the open state to ∼5.0 kcal/mol in the closed one (Kasimova et al., 2018). Thus our findings suggest that in TRPV1, the selectivity filter can act as a gate whose opening is coupled to the conformation of Y671 and N676. Can this coupling mechanism be conserved in other TRP channels? The analysis of the multiple sequence alignment argues against this hypothesis by showing a lack of conservation at position 671 (Fig. 3 A). Moreover, recent structural data indicate that in TRPML3, the selectivity filter adopts the same conformation in the closed and open states (Hirschi et al., 2017; Zhou et al., 2017).
Are the peripheral cavities the missing piece of the puzzle?
The molecular mechanism suggested here agrees with the one emerging from the TRPV3, TRPV5, and TRPV6 structures (Saotome et al., 2016; Singh et al., 2017, 2018a,b; Hughes et al., 2018a,b; McGoldrick et al., 2018) and complements it with a detailed description of the conformational rearrangements involved in activation: the rotation of the conserved asparagine increases the hydrophilicity of the pore surface, thereby promoting the hydration of this compartment. An interesting finding provided by the molecular dynamics simulations and not anticipated by these structures is the factor stabilizing the closed state, in which the asparagine is located outside the pore. Our simulations have shown that in the closed state of TRPV1, N676 is accommodated inside a small hydrophobic cavity located between the S4–S5 linker and S6 and connected to the intracellular solution (Kasimova et al., 2018). This cavity, which we call peripheral in contrast to the central cavity along the ionic pathway, hosts several water molecules forming hydrogen bonds with N676 carboxamide. Importantly, the hydration of the peripheral cavity and the N676 rotation are correlated (Kasimova et al., 2018): when the peripheral cavity dehydrates, N676 rotates inside the pore. Based on this observation, we suggested that at least for some stimuli, activation of TRPV1 is triggered by the dewetting of the peripheral cavities (Kasimova et al., 2018). Is there any experimental evidence that supports the existence of these peripheral cavities? Overall, we found three independent observations, which can be retrospectively interpreted in light of this suggested molecular mechanism: (i) the experiments probing S6 accessibility to intracellular solution (Salazar et al., 2009), (ii) the site-directed mutagenesis experiments targeting S6 residues (Susankova et al., 2007), and (iii) the structural studies of the evolutionary related voltage-gated sodium and calcium channels (Payandeh et al., 2011, 2012; Zhang et al., 2012; Tsai et al., 2013; Tang et al., 2014, 2016; Ahuja et al., 2015; Wu et al., 2015, 2016; Lenaeus et al., 2017; Shen et al., 2017; Sula et al., 2017; Yan et al., 2017; Irie et al., 2018).
In a milestone paper by Salazar et al. (2009), the authors probed the accessibility of the TRPV1 pore residues to the intracellular solution; they mutated, one by one, all the S6 residues to cysteine and assessed their accessibility to thiol-modifying reagents in the open and closed states. Their results showed that modification of D654-Y671 is more prominent in the open state than in the closed one, and that of M682-G683 is state-independent. Interestingly, the rest of S6, namely the I672-L681 segment, showed different accessibility only when a bulky reagent was used; application of small reagents in the open or closed states, in contrast, resulted in a similar effect. Based on these data, the authors suggested the presence of lower and upper gates located at the levels of L681 and Y671, respectively; both gates are open in the open state, while in the closed state, the upper gate is closed and the lower gate is partially open (Salazar et al., 2009). This elegant explanation, however, turned out to be inconsistent with the structural information subsequently obtained via cryo-EM on the apo state (Liao et al., 2013), in which the gate at the level of the S6 bundle crossing was found to be tightly constricted with a radius of ∼1.0 Å, i.e., impermeable to any solute.
Can the presence of the peripheral cavities solve this apparent contradiction? To answer this question, we characterized the accessibility of the TRPV1 pore residues in the open and two closed states, including capsaicin-bound closed and apo closed (Fig. 5). In agreement with Salazar et al. (2009), we found that all of the S6 residues are accessible to the intracellular solution in the capsaicin-bound open state. Moreover, we found that the same set of residues, namely the I672-G683 segment, is accessible in the capsaicin-bound closed state. Despite the constriction at the level of I679, the residues up to and including I672 do interact with the intracellular solution. However, the interacting water molecules are not located in the pore, but rather in the peripheral cavities. Therefore, considering that the pore residues can be accessible to the intracellular solution through the peripheral cavities, one can resolve the contradiction between the water accessibility experiments (Salazar et al., 2009) and the structural information (Liao et al., 2013).
Interestingly, several pore residues above the gate are also accessible in the apo closed state, and similar to the capsaicin-bound closed state, they interact with the water molecules present in the peripheral cavities. However, in the apo state, the peripheral cavities are significantly smaller and host fewer water molecules. We hypothesize that this difference is the result of the S4–S5 linker movement: upon binding, capsaicin pulls the S4–S5 linker in the outward direction, thus expanding the space between the S4–S5 linker and S6 (Yang et al., 2015; Gao et al., 2016). The smaller size of the peripheral cavities in the apo closed state explains why only some of the I672-G683 residues are accessible to the intracellular solution.
Another piece of evidence lending credibility to our hypothesis comes from the alanine-scanning mutagenesis of the S6 residues performed by Susankova et al. (2007). The authors have identified a clear periodic pattern of the functional effects of mutations, which is consistent with the α-helical structure of S6 (Susankova et al., 2007). Subsequently obtained cryo-EM structures of TRPV1 showed that the residues whose mutation has the most significant effect on activation are located outside the pore and line the surface of the peripheral cavities (Fig. 6). Importantly, the response of the corresponding mutants to capsaicin, voltage, and heat was significantly reduced, indicating stabilization of the closed state versus the open one. A possible interpretation for these data is that the substitution of hydrophobic residues (I672, L673, L674, L678, and M682) with a less hydrophobic one (alanine) promotes the hydration of the peripheral cavities and thus stabilizes the closed state.
Finally, we analyzed the families of voltage-gated sodium and calcium channels (Nav and Cav) that are structurally homologous and evolutionarily related to TRP channels (Payandeh et al., 2011, 2012; Zhang et al., 2012; Tsai et al., 2013; Tang et al., 2014, 2016; Ahuja et al., 2015; Wu et al., 2015, 2016; Lenaeus et al., 2017; Shen et al., 2017; Sula et al., 2017; Yan et al., 2017; Irie et al., 2018). Interestingly, these families also show a predominantly hydrophobic pore with the only exception of a conserved asparagine at the position corresponding to N676 (Fig. 7 A). In NaChBac, mutagenesis of this residue into several other amino acids yielded functional channels only in the case of the N225D mutant, and even in that case, gating of N225D was drastically affected (O’Reilly et al., 2017). Furthermore, the structures of Nav and Cav channels show that the conserved asparagine can adopt both the pore inward-facing (NavPas, all subunits; NavEel, one subunit out of four) and outward-facing conformations (Cav1.1, all subunits; NavAb, NavRh, NavMs, and NavEel, three subunits out of four; Fig. 7, B and C; Payandeh et al., 2011, 2012; Zhang et al., 2012; Tsai et al., 2013; Tang et al., 2014, 2016; Ahuja et al., 2015; Wu et al., 2015, 2016; Lenaeus et al., 2017; Shen et al., 2017; Sula et al., 2017; Yan et al., 2017; Irie et al., 2018). Do these channels have the peripheral cavities? A recent high-resolution structure of the bacterial sodium channel NavMs shows the presence of four cavities, whose location coincides with that of the peripheral ones in TRPV1 (Sula et al., 2017). Moreover, similar to the peripheral cavities, they host the side chain of the conserved asparagine and, importantly, few water molecules interacting with it.
Despite two decades of extensive investigations, the molecular mechanism of TRPV1 activation is only partially understood. Catching an ion channel in the act of gating is, in general, a complex task requiring techniques with both high temporal and high spatial resolution. NMR spectroscopy and fluorescence-/luminescence-based approaches are typically the tools of choice, but they cannot easily probe the motions of an ion channel when embedded in a cell membrane. This limitation is particularly severe in the case of TRPV1, whose activity shows exquisite dependency on numerous environmental factors, such as lipid composition (Chuang et al., 2001; Cao et al., 2013a; Lukacs et al., 2013; Poblete et al., 2015; Rohacs, 2015; Gao et al., 2016), extracellular concentration of Na+ (Jara-Oseguera et al., 2016), and the presence of cholesterol or fatty acids (Liu et al., 2006; Matta et al., 2007; Picazo-Juárez et al., 2011; Morales-Lázaro et al., 2016). Accordingly, advanced experimental techniques, such as cell unroofing (Gordon et al., 2016; Zagotta et al., 2016) and incorporation of novel unnatural amino acids (Zagotta et al., 2016; Steinberg et al., 2017), are being constantly developed to overcome these limitations and introduce only minimal perturbations to the native environment of the channel. Despite the great progress, these techniques have focused so far only on specific aspects of TRPV1 activation. A global picture of this molecular process is still lacking, and it is unlikely that a single experimental technique will provide it in the foreseeable future.
Here, prompted by our previous molecular dynamics simulations (Kasimova et al., 2018) and bioinformatics analyses (Palovcak et al., 2015), we propose a molecular mechanism for TRPV1 activation. More than definite answers, our contribution provides testable predictions and an overall framework to rationalize a wealth of existing experimental data. Indeed, despite being consistent with the results of several independent studies, our hypothesis is still in part speculative. Most importantly, our work was only focused on the capsaicin-bound state, and the transition between that and the apo one was not investigated. It remains unclear how the conformation of the peripheral cavities changes along this transition and how capsaicin binding triggers this change. Based on the previous studies (Yang et al., 2015; Gao et al., 2016), we surmise that capsaicin promotes the expansion of the peripheral cavities between S6 and the S4–S5 linker. Enhanced sampling simulations starting from the apo state could provide answers to this question and also assess possible temperature dependency of the apo to capsaicin-bound transition.
Further research will be required to validate the molecular mechanism proposed here, including design of novel experiments able to discriminate among multiple alternative interpretations. For instance, incorporating unnatural amino acids in TRPV1 may allow for a fine-tuning of the polar properties of the pore and the peripheral cavities without perturbing their geometry too much. The same approach could be used to investigate the coupling between the selectivity filter and the gate. According to our model, this coupling involves the interaction between the side chain of N676 and the S6 helix main chain. The latter could be modified by introducing appropriately designed unnatural amino acids. From the computational point of view, long unbiased simulations could be used to sample gating events and provide direct confirmation of the activation mechanism. The Anton 2 supercomputer developed by D. E. Shaw Research (Shaw et al., 2009, 2014) and available through the Pittsburgh Supercomputing Center provides the opportunity to reach timescales comparable to those required for TRPV1 activation and thus to test the hypothesized chain of events. These simulations ought to involve, ideally, both the apo and capsaicin-bound states.
In closing, it is interesting to comment on an intriguing possibility emerging from our analyses: we found that the key structural features supporting our hypothesized mechanism are present also in the evolutionarily related families of voltage-gated sodium and calcium channels. In particular, the conserved asparagine in the S6 helix, the π-bulge, and the peripheral cavities are seemingly conserved beyond the family of TRP channels. This raises the possibility that voltage-gated sodium and calcium channels might share with TRP channels important aspects of the activation mechanism.
However, the sequence and structure conservation of S6 across TRP and voltage-gated sodium and calcium channels highlights a potentially problematic aspect of our hypothesis: in our previous computational investigation (Kasimova et al., 2018), we suggested that the rotation of N676 underlies temperature sensitivity in TRPV1. The evolutionary conservation of this asparagine cannot be easily reconciled with the fact that several TRP and voltage-gated sodium and calcium channels are insensitive to temperature. Future investigations ought to improve our understanding of the temperature sensitivity of ion channels and thus contribute to solving this conundrum.
This work was partially supported by the National Institutes of Health through grants R01GM093290, S10OD020095, NS055159, and P01GM055876 and the National Science Foundation through grants ACI-1614804 and CNS-1625061.
The authors declare no competing financial interests.
Author contributions: M.A. Kasimova and V. Carnevale designed the study. M.A. Kasimova collected the data and carried out analyses and simulations. M.A. Kasimova and D. Granata developed computational code to analyze molecular dynamics trajectories. M.A. Kasimova, A.T. Yazici, Y. Yudin, D. Granata, M.L. Klein, T. Rohacs, and V. Carnevale analyzed and interpreted the outcome of bioinformatics analyses and molecular dynamics simulations. T. Rohacs and V. Carnevale secured funding for the project. M.A. Kasimova and V. Carnevale wrote the paper with input from the other coauthors.
José D. Faraldo-Gómez served as editor.