Ion Channel Sensing: Are Fluctuations the Crux of the Matter?

Dynamic molecular portraits of ion-conducting pores characterize functional states of TRPV channels

Structural biology is solving an ever-increasing number of snapshots of ion channel conformational ensembles. Deciphering ion channel mechanisms, however, requires understanding the ensemble dynamics beyond the static structures. Here, we present a molecular modeling-based approach characterizing the ion channel structural intermediates, or their "dynamic molecular portraits", by assessing water and ion conductivity along with the detailed evaluation of pore hydrophobicity and residue packing. We illustrate the power of this approach by analyzing structures of few vanilloid-subfamily transient receptor potential (TRPV) channels. Based on the pore architecture, there are three major states that are common for TRPVs, which we call α-closed, π-closed, and π-open. We show that the pore hydrophobicity and residue packing for the open state is most favorable for the pore conductance. On the contrary, the α-closed state is the most hydrophobic and always non-conducting. Our approach can also be used for structural and functional classification of ion channels.

Functional dynamics and allosteric modulation of TRPA1

The cation channel TRPA1 is a potentially important drug target, and characterization of TRPA1 functional dynamics might help guide structure-based drug design. Here, we present results from long-timescale molecular dynamics simulations of TRPA1 with an allosteric activator, allyl isothiocyanate (AITC), in which we observed spontaneous transitions from a closed, non-conducting channel conformation into an open, conducting conformation. Based on these transitions, we propose a gating mechanism in which movement of a regulatory TRP-like domain allosterically translates into pore opening in a manner reminiscent of pore opening in voltage-gated ion channels. In subsequent experiments, we found that mutations that disrupt packing of the S4-S5 linker-TRP-like domain and the S5 and S6 helices also affected channel activity. In simulations, we also observed A-967079, a known allosteric inhibitor, binding between helices S5 and S6, suggesting that A-967079 may suppress activity by stabilizing a non-conducting pore conformation-a finding consistent with our proposed gating mechanism.

Atomic structure of the open SARS-CoV-2 E viroporin

The envelope (E) protein of the SARS-CoV-2 virus forms cation-conducting channels in the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of infected cells. The calcium channel activity of E is associated with the inflammatory responses of COVID-19. Using solid-state NMR (ssNMR) spectroscopy, we have determined the open-state structure of E's transmembrane domain (ETM) in lipid bilayers. Compared to the closed state, open ETM has an expansive water-filled amino-terminal chamber capped by key glutamate and threonine residues, a loose phenylalanine aromatic belt in the middle, and a constricted polar carboxyl-terminal pore filled with an arginine and a threonine residue. This structure gives insights into how protons and calcium ions are selected by ETM and how they permeate across the hydrophobic gate of this viroporin.

Temperature sensitive contact modes allosterically gate TRPV3

TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.

Temperature-Sensitive Contact Modes Allosterically Gate TRPV3

TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.

Inner pore hydration free energy controls the activation of big potassium channels

Hydrophobic gating is an emerging mechanism in regulation of protein ion channels where the pore remains physically open but becomes dewetted to block ion permeation. Atomistic molecular dynamics simulations have played a crucial role in understanding hydrophobic gating by providing the molecular details to complement mutagenesis and structural studies. However, existing studies rely on direct simulations and do not quantitatively describe how the sequence and structural changes may control the delicate liquid-vapor equilibrium of confined water in the pore of the channel protein. To address this limitation, we explore two enhanced sampling methods, namely metadynamics and umbrella sampling, to derive free-energy profiles of pore hydration in both the closed and open states of big potassium (BK) channels, which are important in cardiovascular and neural systems. It was found that metadynamics required substantially longer sampling times and struggled to generate stably converged free-energy profiles due to the slow dynamics of cooperative pore water diffusion even in the barrierless limit. Using umbrella sampling, well-converged free-energy profiles can be readily generated for the wild-type BK channels as well as three mutants with pore-lining mutations experimentally known to dramatically perturb the channel gating voltage. The results show that the free energy of pore hydration faithfully reports the gating voltage of the channel, providing further support for hydrophobic gating in BK channels. Free-energy analysis of pore hydration should provide a powerful approach for quantitative studies of how protein sequence, structure, solution conditions, and/or drug binding may modulate hydrophobic gating in ion channels.

Structural Mechanism of Ionic Conductivity of the TRPV1 Channel

The so-called "hydrophobic gating" is widely discussed as a putative mechanism to control water and ion conduction via ion channels. This effect can occur in narrow areas of the channels pore lined by non-polar residues. In the closed state of the channel, such regions may spontaneously transit to a dehydrated state to block water and ions transport without full pore occlusion. In the open state, the hydrophobic gate is wide enough to provide sustainable hydration and conduction. Apparently, the transport through the open hydrophobic gate may by facilitated by some polar residues that assist polar/charged substances to overcome the energy barrier created by nonpolar environment. In this work, we investigated the behavior of Na⁺ ions and their hydration shells in the open pore of the rat TRPV1 ion channel by molecular dynamics simulations. We show that polar protein groups coordinate water molecules in such a way as to restore the hydration shell of ions in the hydrophobic gate that ensures ion transport through the gate in a fully hydrated state.

An α-π transition in S6 shapes the conformational cycle of the bacterial sodium channel NavAb

Voltage-gated sodium channels play an important role in electrical signaling in excitable cells. In response to changes in membrane potential, they cycle between nonconducting and conducting conformations. With recent advances in structural biology, structures of sodium channels have been captured in several distinct conformations, which are thought to represent different functional states. However, it has been difficult to capture the intrinsically transient open state. We recently showed that a proposed open state of the bacterial sodium channel NavMs was not conductive and that a conformational change involving a transition to a π-helix in the pore-lining S6 helix converted this structure into a conducting state. However, the relevance of this structural feature in other sodium channels, and its implications for the broader gating cycle, remained unclear. Here, we propose a comparable open state of another class of bacterial channel from Aliarcobacter butzleri (NavAb) with characteristic pore hydration, ion permeation, and drug binding properties. Furthermore, we show that a π-helix transition can lead to pore opening and that such a conformational change blocks fenestrations in the inner helix bundle. We also discover that a region in the C-terminal domain can undergo a disordering transition proposed to be important for pore opening. These results support a role for a π-helix transition in the opening of NavAb, enabling new proposals for the structural annotation and drug modulation mechanisms in this important sodium channel model.

When is a hydrophobic gate not a hydrophobic gate?

The flux of ions through a channel is most commonly regulated by changes that result in steric occlusion of its pore. However, ion permeation can also be prevented by formation of a desolvation barrier created by hydrophobic residues that line the pore. As a result of relatively minor structural changes, confined hydrophobic regions in channels may undergo transitions between wet and dry states to gate the pore closed without physical constriction of the permeation pathway. This concept is referred to as hydrophobic gating, and many examples of this process have been demonstrated. However, the term is also now being used in a much broader context that often deviates from its original meaning. In this Viewpoint, we explore the formal definition of a hydrophobic gate, discuss examples of this process compared with other gating mechanisms that simply exploit hydrophobic residues and/or lipids in steric closure of the pore, and describe the best practice for identification of a hydrophobic gate.

Sequence and structural conservation reveal fingerprint residues in TRP channels

Transient receptor potential (TRP) proteins are a large family of cation-selective channels, surpassed in variety only by voltage-gated potassium channels. Detailed molecular mechanisms governing how membrane voltage, ligand binding, or temperature can induce conformational changes promoting the open state in TRP channels are still a matter of debate. Aiming to unveil distinctive structural features common to the transmembrane domains within the TRP family, we performed phylogenetic reconstruction, sequence statistics, and structural analysis over a large set of TRP channel genes. Here, we report an exceptionally conserved set of residues. This fingerprint is composed of twelve residues localized at equivalent three-dimensional positions in TRP channels from the different subtypes. Moreover, these amino acids are arranged in three groups, connected by a set of aromatics located at the core of the transmembrane structure. We hypothesize that differences in the connectivity between these different groups of residues harbor the apparent differences in coupling strategies used by TRP subgroups.

Venom Peptide Toxins Targeting the Outer Pore Region of Transient Receptor Potential Vanilloid 1 in Pain: Implications for Analgesic Drug Development

The transient receptor potential vanilloid 1 (TRPV1) ion channel plays an important role in the peripheral nociceptive pathway. TRPV1 is a polymodal receptor that can be activated by multiple types of ligands and painful stimuli, such as noxious heat and protons, and contributes to various acute and chronic pain conditions. Therefore, TRPV1 is emerging as a novel therapeutic target for the treatment of various pain conditions. Notably, various peptides isolated from venomous animals potently and selectively control the activation and inhibition of TRPV1 by binding to its outer pore region. This review will focus on the mechanisms by which venom-derived peptides interact with this portion of TRPV1 to control receptor functions and how these mechanisms can drive the development of new types of analgesics.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.

Dual regulation of TRPV1 channels by phosphatidylinositol via functionally distinct binding sites

Regulation of the heat- and capsaicin-activated transient receptor potential vanilloid 1 (TRPV1) channel by phosphoinositides is complex and controversial. In the most recent TRPV1 cryo-EM structure, endogenous phosphatidylinositol (PtdIns) was detected in the vanilloid binding site, and phosphoinositides were proposed to act as competitive vanilloid antagonists. This model is difficult to reconcile with phosphatidylinositol 4,5-bisphosphate [PtdIns(4,5)P₂] being a well-established positive regulator of TRPV1. Here we show that in the presence of PtdIns(4,5)P₂ in excised patches, PtdIns, but not PtdIns(4)P, partially inhibited TRPV1 activity at low, but not at high capsaicin concentrations. This is consistent with PtdIns acting as a competitive vanilloid antagonist. However, in the absence of PtdIns(4,5)P₂, PtdIns partially stimulated TRPV1 activity. We computationally identified residues, which are in contact with PtdIns, but not with capsaicin in the vanilloid binding site. The I703A mutant of TRPV1 showed increased sensitivity to capsaicin, as expected when removing the effect of an endogenous competitive antagonist. I703A was not inhibited by PtdIns in the presence of PtdIns(4,5)P₂, but it was still activated by PtdIns in the absence of PtdIns(4,5)P₂ indicating that inhibition, but not activation by PtdIns proceeds via the vanilloid binding site. In molecular dynamics simulations, PtdIns was more stable than PtdIns(4,5)P₂ in this inhibitory site, whereas PtdIns(4,5)P₂ was more stable than PtdIns in a previously identified, nonoverlapping, putative activating binding site. Our data indicate that phosphoinositides regulate channel activity via functionally distinct binding sites, which may explain some of the complexities of the effects of these lipids on TRPV1.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Electric Field Induced Wetting of a Hydrophobic Gate in a Model Nanopore Based on the 5-HT3 Receptor Channel

In this study we examined the influence of a transmembrane voltage on the hydrophobic gating of nanopores using molecular dynamics simulations. We observed electric field induced wetting of a hydrophobic gate in a biologically inspired model nanopore based on the 5-HT3 receptor in its closed state, with a field of at least ∼100 mV nm-1 (corresponding to a supra-physiological potential difference of ∼0.85 V across the membrane) required to hydrate the pore. We also found an unequal distribution of charged residues can generate an electric field intrinsic to the nanopore which, depending on its orientation, can alter the effect of the external field, thus making the wetting response asymmetric. This wetting response could be described by a simple model based on water surface tension, the volumetric energy contribution of the electric field, and the influence of charged amino acids lining the pore. Finally, the electric field response was used to determine time constants characterizing the phase transitions of water confined within the nanopore, revealing liquid-vapor oscillations on a time scale of ∼5 ns. This time scale was largely independent of the water model employed and was similar for different sized pores representative of the open and closed states of the pore. Furthermore, our finding that the threshold voltage required for hydrating a hydrophobic gate depends on the orientation of the electric field provides an attractive perspective for the design of rectifying artificial nanopores.

The mechanism and energetics of a ligand-controlled hydrophobic gate in a mammalian two pore channel

In the last decade two-pore intracellular channels (TPCs) attracted the interest of researchers, still some key questions remain open. Their importance for vacuolar (plants) and endo-lysosomal (animals) function highlights them as a very attractive system to study, both theoretically and experimentally. Indicated as key players in the trafficking of the cell, today they are considered a new potential target for avoiding virus infections, including those from coronaviruses. A particular boost for theoretical examinations has been made with recent high-resolution X-ray and cryo-EM structures. These findings have opened the way for efficient and precise computational studies at the atomistic level. Here we report a set of multiscale-calculations performed on the mTPC1, a ligand- and voltage-gated sodium selective channel. The molecular dynamics and enhanced molecular dynamics simulations were used for a thorough analysis of the mammalian TPC1 behaviour in the presence and absence of the ligand molecule, with a special accent on the supposed bottleneck, the hydrophobic gate. Moreover, from the reconstructed free energy obtained from enhanced simulations, we have calculated the macroscopic conductance of sodium ions through the mTPC1, which we compared with measured single-channel conductance values. The hydrophobic gate works as a steric barrier and the key parameters are its flexibility and the dimension of the sodium first hydration shell.

Structure and function of the calcium-selective TRP channel TRPV6

Epithelial calcium channel TRPV6 is a member of the vanilloid subfamily of TRP channels that is permeable to cations and highly selective to Ca²⁺ ; it shows constitutive activity regulated negatively by Ca²⁺ and positively by phosphoinositol and cholesterol lipids. In this review, we describe the molecular structure of TRPV6 and discuss how its structural elements define its unique functional properties. High Ca²⁺ selectivity of TRPV6 originates from the narrow selectivity filter, where Ca²⁺ ions are directly coordinated by a ring of anionic aspartate side chains. Divalent cations Ca²⁺ and Ba²⁺ permeate TRPV6 pore according to the knock-off mechanism, while tight binding of Gd³⁺ to the aspartate ring blocks the channel and prevents Na⁺ from permeating the pore. The iris-like channel opening is accompanied by an α-to-π helical transition in the pore-lining transmembrane helix S6. As a result of this transition, the intracellular halves of the S6 helices bend and rotate by about 100 deg, exposing different residues to the channel pore in the open and closed states. Channel opening is also associated with changes in occupancy of the transmembrane domain lipid binding sites. The inhibitor 2-aminoethoxydiphenyl borate (2-APB) binds to TRPV6 in a pocket formed by the cytoplasmic half of the S1-S4 transmembrane helical bundle and shifts open-closed channel equilibrium towards the closed state by outcompeting lipids critical for activation. Ca²⁺ inhibits TRPV6 via binding to calmodulin (CaM), which mediates Ca²⁺ -dependent inactivation. The TRPV6-CaM complex exhibits 1:1 stoichiometry; one TRPV6 tetramer binds both CaM lobes, which adopt a distinct head-to-tail arrangement. The CaM C-terminal lobe plugs the channel through a unique cation-π interaction by inserting the side chain of lysine K115 into a tetra-tryptophan cage at the ion channel pore intracellular entrance. Recent studies of TRPV6 structure and function described in this review advance our understanding of the role of this channel in physiology and pathophysiology and inform new therapeutic design.

Structural insights into the gating mechanisms of TRPV channels

Transient Receptor Potential channels from the vanilloid subfamily (TRPV) are a group of cation channels modulated by a variety of endogenous stimuli as well as a range of natural and synthetic compounds. Their roles in human health make them of keen interest, particularly from a pharmacological perspective. However, despite this interest, the complexity of these channels has made it difficult to obtain high resolution structures until recently. With the cryo-EM resolution revolution, TRPV channel structural biology has blossomed to produce dozens of structures, covering every TRPV family member and a variety of approaches to examining channel modulation. Here, we review all currently available TRPV structures and the mechanistic insights into gating that they reveal.

TRPA1 modulation by piperidine carboxamides suggests an evolutionarily conserved binding site and gating mechanism

The TRPA1 channel functions as an irritant sensor and is a therapeutic target for treating pain, itch, and respiratory diseases. TRPA1 can be activated by electrophilic compounds via covalent modification or activated by noncovalent agonists via ligand binding. However, how covalent modification leads to channel opening and, importantly, how noncovalent binding activates TRPA1 are not well-understood. Here we identified a group of noncovalent agonists and used them to explore TRPA1 gating through iterative functional analyses, molecular modeling, and structure–activity relationship studies. We show that TRPA1 possesses an evolutionarily conserved ligand binding site common to other TRP channels. The combination of computational modeling and experimental structure–activity data lays the foundations for rational drug design.

The transient receptor potential ankyrin 1 (TRPA1) channel functions as an irritant sensor and is a therapeutic target for treating pain, itch, and respiratory diseases. As a ligand-gated channel, TRPA1 can be activated by electrophilic compounds such as allyl isothiocyanate (AITC) through covalent modification or activated by noncovalent agonists through ligand binding. However, how covalent modification leads to channel opening and, importantly, how noncovalent binding activates TRPA1 are not well-understood. Here we report a class of piperidine carboxamides (PIPCs) as potent, noncovalent agonists of human TRPA1. Based on their species-specific effects on human and rat channels, we identified residues critical for channel activation; we then generated binding modes for TRPA1–PIPC interactions using structural modeling, molecular docking, and mutational analysis. We show that PIPCs bind to a hydrophobic site located at the interface of the pore helix 1 (PH1) and S5 and S6 transmembrane segments. Interestingly, this binding site overlaps with that of known allosteric modulators, such as A-967079 and propofol. Similar binding sites, involving π-helix rearrangements on S6, have been recently reported for other TRP channels, suggesting an evolutionarily conserved mechanism. Finally, we show that for PIPC analogs, predictions from computational modeling are consistent with experimental structure–activity studies, thereby suggesting strategies for rational drug design.

The Contribution of the Ankyrin Repeat Domain of TRPV1 as a Thermal Module

The TRPV1 cation nonselective ion channel plays an essential role in thermosensation and perception of other noxious stimuli. TRPV1 can be activated by low extracellular pH, high temperature, or naturally occurring pungent molecules such as allicin, capsaicin, or resiniferatoxin. Its noxious thermal sensitivity makes it an important participant as a thermal sensor in mammals. However, details of the mechanism of channel activation by increases in temperature remain unclear. Here, we used a combination of approaches to try to understand the role of the ankyrin repeat domain (ARD) in channel behavior. First, a computational modeling approach by coarse-grained molecular dynamics simulation of the whole TRPV1 embedded in a phosphatidylcholine and phosphatidylethanolamine membrane provides insight into the dynamics of this channel domain. Global analysis of the structural ensemble shows that the ARD is a region that sustains high fluctuations during dynamics at different temperatures. We then performed biochemical and thermal stability studies of the purified ARD by the means of circular dichroism and tryptophan fluorescence and demonstrate that this region undergoes structural changes at similar temperatures that lead to TRPV1 activation. Our data suggest that the ARD is a dynamic module and that it may participate in controlling the temperature sensitivity of TRPV1.

Evolutionarily Conserved Interactions within the Pore Domain of Acid-Sensing Ion Channels

Despite the sequence homology between acid-sensing ion channels (ASICs) and epithelial sodium channel (ENaCs), these channel families display very different functional characteristics. Whereas ASICs are gated by protons and show a relatively low degree of selectivity for sodium over potassium, ENaCs are constitutively active and display a remarkably high degree of sodium selectivity. To decipher if some of the functional diversity originates from differences within the transmembrane helices (M1 and M2) of both channel families, we turned to a combination of computational and functional interrogations, using statistical coupling analysis and mutational studies on mouse ASIC1a. The coupling analysis suggests that the relative position of M1 and M2 in the upper part of the pore domain is likely to remain constant during the ASIC gating cycle, whereas they may undergo relative movements in the lower part. Interestingly, our data suggest that to account for coupled residue pairs being in close structural proximity, both domain-swapped and nondomain-swapped ASIC M2 conformations need to be considered. Such conformational flexibility is consistent with structural work, which suggested that the lower part of M2 can adopt both domain-swapped and nondomain-swapped conformations. Overall, mutations to residues in the middle and lower pore were more likely to affect gating and/or ion selectivity than those in the upper pore. Indeed, disrupting the putative interaction between a highly conserved Trp/Glu residue pair in the lower pore is detrimental to gating and selectivity, although this interaction might occur in both domain-swapped and nonswapped conformations. Finally, our results suggest that the greater number of larger, aromatic side chains in the ENaC M2 helix may contribute to the constitutive activity of these channels at a resting pH. Together, the data highlight differences in the transmembrane domains of these closely related ion channels that may help explain some of their distinct functional properties.

Confined Dynamics of Water in Transmembrane Pore of TRPV1 Ion Channel

Solvation effects play a key role in chemical and biological processes. The microscopic properties of water near molecular surfaces are radically different from those in the bulk. Furthermore, the behavior of water in confined volumes of a nanometer scale, including transmembrane pores of ion channels, is especially nontrivial. Knowledge at the molecular level of structural and dynamic parameters of water in such systems is necessary to understand the mechanisms of ion channels functioning. In this work, the results of molecular dynamics (MD) simulations of water in the pore and selectivity filter domains of TRPV1 (Transient Receptor Potential Vanilloid type 1) membrane channel are considered. These domains represent nanoscale volumes with strongly amphiphilic walls, where physical behavior of water radically differs from that of free hydration (e.g., at protein interfaces) or in the bulk. Inside the pore and filter domains, water reveals a very heterogeneous spatial distribution and unusual dynamics: It forms compact areas localized near polar groups of particular residues. Residence time of water molecules in such areas is at least 1.5 to 3 times larger than that observed for similar groups at the protein surface. Presumably, these water “blobs” play an important role in the functional activity of TRPV1. In particular, they take part in hydration of the hydrophobic TRPV1 pore by localizing up to six waters near the so-called “lower gate” of the channel and reducing by this way the free energy barrier for ion and water transport. Although the channel is formed by four identical protein subunits, which are symmetrically packed in the initial experimental 3D structure, in the course of MD simulations, hydration of the same amino acid residues of individual subunits may differ significantly. This greatly affects the microscopic picture of the distribution of water in the channel and, potentially, the mechanism of its functioning. Therefore, reconstruction of the full picture of TRPV1 channel solvation requires thorough atomistic simulations and analysis. It is important that the naturally occurring porous volumes, like ion-conducting protein domains, reveal much more sophisticated and fine-tuned regulation of solvation than, e.g., artificially designed carbon nanotubes.

Cytoplasmic Inter-Subunit Interface Controls Use-Dependence of Thermal Activation of TRPV3 Channel

The vanilloid transient receptor potential channel TRPV3 is a putative molecular thermosensor widely considered to be involved in cutaneous sensation, skin homeostasis, nociception, and pruritus. Repeated stimulation of TRPV3 by high temperatures above 50 °C progressively increases its responses and shifts the activation threshold to physiological temperatures. This use-dependence does not occur in the related heat-sensitive TRPV1 channel in which responses decrease, and the activation threshold is retained above 40 °C during activations. By combining structure-based mutagenesis, electrophysiology, and molecular modeling, we showed that chimeric replacement of the residues from the TRPV3 cytoplasmic inter-subunit interface (N251-E257) with the homologous residues of TRPV1 resulted in channels that, similarly to TRPV1, exhibited a lowered thermal threshold, were sensitized, and failed to close completely after intense stimulation. Crosslinking of this interface by the engineered disulfide bridge between substituted cysteines F259C and V385C (or, to a lesser extent, Y382C) locked the channel in an open state. On the other hand, mutation of a single residue within this region (E736) resulted in heat resistant channels. We propose that alterations in the cytoplasmic inter-subunit interface produce shifts in the channel gating equilibrium and that this domain is critical for the use-dependence of the heat sensitivity of TRPV3.

A heuristic derived from analysis of the ion channel structural proteome permits the rapid identification of hydrophobic gates

Ion channels are nanoscale protein pores in cell membranes. An exponentially increasing number of structures for channels means that computational methods for predicting their functional state are needed. Hydrophobic gates in ion channels result in local dewetting of pores, which functionally closes them to water and ion permeation. We use simulations of water behavior within nearly 200 different ion channel structures to explore how the radius and hydrophobicity of pores determine their hydration vs. dewetting behavior. Machine learning-assisted analysis of these simulations allowed us to propose a simple model for this relationship and present an easy method for rapidly predicting the functional state of new channel structures as they emerge.

Ion channel proteins control ionic flux across biological membranes through conformational changes in their transmembrane pores. An exponentially increasing number of channel structures captured in different conformational states are now being determined; however, these newly resolved structures are commonly classified as either open or closed based solely on the physical dimensions of their pore, and it is now known that more accurate annotation of their conductive state requires additional assessment of the effect of pore hydrophobicity. A narrow hydrophobic gate region may disfavor liquid-phase water, leading to local dewetting, which will form an energetic barrier to water and ion permeation without steric occlusion of the pore. Here we quantify the combined influence of radius and hydrophobicity on pore dewetting by applying molecular dynamics simulations and machine learning to nearly 200 ion channel structures. This allows us to propose a simple simulation-free heuristic model that rapidly and accurately predicts the presence of hydrophobic gates. This not only enables the functional annotation of new channel structures as soon as they are determined, but also may facilitate the design of novel nanopores controlled by hydrophobic gates.

Heat activation mechanism of TRPV1: New insights from molecular dynamics simulation

As a member of the transient receptor potential (TRP) channels superfamily, the TRPV1 channel undergoes a closed-to-open gating transition in response to various physical and chemical stimuli including heat. Thanks to recent progress in cryo-electron microscopy, high-resolution structures are becoming available for various TRP channels including TRPV1. This has enabled us to study the molecular mechanism of TRPV1 channel gating by using molecular simulation. Here we review recent progress in molecular simulations of TRPV1 channel by us and others, with focus on our molecular dynamics (MD) simulations of TRPV1 at different temperatures. While no consensus has been reached on the heat activation mechanism of TRPV1, the simulations have offered specific predictions and models for future experimental studies to test.

Conformational ensemble of the human TRPV3 ion channel

Transient receptor potential vanilloid channel 3 (TRPV3), a member of the thermosensitive TRP (thermoTRPV) channels, is activated by warm temperatures and serves as a key regulator of normal skin physiology through the release of pro-inflammatory messengers. Mutations in trpv3 have been identified as the cause of the congenital skin disorder, Olmsted syndrome. Unlike other members of the thermoTRPV channel family, TRPV3 sensitizes upon repeated stimulation, yet a lack of structural information about the channel precludes a molecular-level understanding of TRPV3 sensitization and gating. Here, we present the cryo-electron microscopy structures of apo and sensitized human TRPV3, as well as several structures of TRPV3 in the presence of the common thermoTRPV agonist 2-aminoethoxydiphenyl borate (2-APB). Our results show α-to-π-helix transitions in the S6 during sensitization, and suggest a critical role for the S4-S5 linker π-helix during ligand-dependent gating.

A hypothetical molecular mechanism for TRPV1 activation that invokes rotation of an S6 asparagine

TRPV1 channels comprise four subunits containing six transmembrane segments (S1–S6) that surround a central pore. Kasimova et al. hypothesize that channel opening involves rotation of an S6 asparagine residue toward the pore, as well as associated pore hydration and external cavity dehydration.

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.

Determining the molecular basis of voltage sensitivity in membrane proteins

The identification of voltage-sensing elements in membrane proteins is challenging due to the diversity of voltage-sensing mechanisms. Kasimova et al. present a computational approach to predict the elements involved in voltage sensing, which they validate using voltage-gated ion channels.

Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mechanical work. They are responsible for a spectrum of physiological processes in living organisms, including electrical signaling and cell-cycle progression. Although the mechanism of voltage-sensing has been well characterized for some membrane proteins, including voltage-gated ion channels, even the location of the voltage-sensing elements remains unknown for others. Moreover, the detection of these elements by using experimental techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estimation of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis experiments will significantly reduce the number of potential voltage-sensitive elements to be tested.