Amiloride Docking to Acid-sensing Ion Channel-1

Amiloride is a small molecule diuretic, which has been used to dissect sodium transport pathways in many different systems. This drug is known to interact with the epithelial sodium channel and acid-sensing ion channel proteins, as well as sodium/hydrogen antiporters and sodium/calcium exchangers. The exact structural basis for these interactions has not been elucidated as crystal structures of these proteins have been challenging to obtain, though some involved residues and domains have been mapped. This work examines the interaction of amiloride with acid-sensing ion channel-1, a protein whose structure is available using computational and experimental techniques. Using molecular docking software, amiloride and related molecules were docked to model structures of homomeric human ASIC-1 to generate potential interaction sites and predict which analogs would be more or less potent than amiloride. The predictions made were experimentally tested using whole-cell patch clamp. Drugs previously classified as NCX or NHE inhibitors are shown to also inhibit hASIC-1. Potential docking sites were re-examined against experimental data to remove spurious interaction sites. The voltage sensitivity of inhibitors was also examined. Using the aggregated data from these computational and experimental experiments, putative interaction sites for amiloride and hASIC-1 have been defined. Future work will experimentally verify these interaction sites, but at present this should allow for virtual screening of drug libraries at these putative interaction sites.     

Piroxicam confer neuroprotection in Cerebral Ischemia by inhibiting Cyclooxygenases,Acid- Sensing Ion Channel-1a and Aquaporin-4: an in silico comparison with Aspirin and Nimesulide

Cerebral ischemia (CI), caused by the deprivation of oxygen and glucose to the brain, is the leading cause of permanent disability. Neuronal demise in CI has been linked to several pathways which include cyclooxygenases (COX) − mediated production of prostaglandins (PGs) and subsequently reactive oxygen species (ROS), aquaporin-4 (AQ-4) − mediated brain edema and acidsensing ion channel-1a (ASIC-1a) − mediated acidotoxicity, matrix remodeling, in addition to others. Several non-steroidal antiinflammatory drugs (NSAIDs) are presently in use to prevent these pathways. However, owing to the large number of processes involved, there is high drug load. So, identifying drugs with multimodal role has always been a frequently sought venture. The present in silico study has been performed to find out the relative efficacy of three different NSAIDs (Piroxicam, Aspirin and Nimesulide) in preventing neurodegeneration in CI, with respect to their inhibitory potential on COXs, AQ-4 and ASIC-1a. We find that piroxicam is the most potent inhibitor of these receptors as compared to the NSAIDs under investigation. Since piroxicam has already been reported to inhibit N-methyl-D-aspartate (NMDA) receptor and matrix metalloproteinases (MMPs), which are also linked to CI-induced neurodegeneration, we hereby propose piroxicam to be a gold-standard drug in preventing neurodegeneration in CI.     

Identification of structures for ion channel kinetic models

Markov models of ion channel dynamics have evolved as experimental advances have improved our understanding of channel function. Past studies have examined limited sets of various topologies for Markov models of channel dynamics. We present a systematic method for identification of all possible Markov model topologies using experimental data for two types of native voltage-gated ion channel currents: mouse atrial sodium currents and human left ventricular fast transient outward potassium currents. Successful models identified with this approach have certain characteristics in common, suggesting that aspects of the model topology are determined by the experimental data. Incorporating these channel models into cell and tissue simulations to assess model performance within protocols that were not used for training provided validation and further narrowing of the number of acceptable models. The success of this approach suggests a channel model creation pipeline may be feasible where the structure of the model is not specified a priori.     

The Schistosoma mansoni tegumental allergen protein,SmTAL1: Binding to an IQ-motif from a voltage-gated ion channel and effects of praziquantel

SmTAL1 is a calcium binding protein from the parasitic worm, Schistosoma mansoni. Structurally it is comprised of two domains – an N-terminal EF-hand domain and a C-terminal dynein light chain (DLC)-like domain. The protein has previously been shown to interact with the anti-schistosomal drug, praziquantel (PZQ). Here, we demonstrated that both EF-hands in the N-terminal domain are functional calcium ion binding sites. The second EF-hand appears to be more important in dictating affinity and mediating the conformational changes which occur on calcium ion binding. There is positive cooperativity between the four calcium ion binding sites in the dimeric form of SmTAL1. Both the EF-hand domain and the DLC-domain dimerise independently suggesting that both play a role in forming the SmTAL1 dimer. SmTAL1 binds non-cooperatively to PZQ and cooperatively to an IQ-motif from SmCa_v1B, a voltage-gated calcium channel. PZQ tends to strengthen this interaction, although the relationship is complex. These data suggest the hypothesis that SmTAL1 regulates at least one voltage-gated calcium channel and PZQ interferes with this process. This may be important in the molecular mechanism of this drug. It also suggests that compounds which bind SmTAL1, such as six from the Medicines for Malaria Box identified in this work, may represent possible leads for the discovery of novel antagonists.     

Constraint-based, homology model of the extracellular domain of the epithelial Na+ channel α subunit reveals a mechanism of channel activation by proteases

The epithelial Na(+) channel (ENaC) mediates Na(+) transport across high resistance epithelia. This channel is assembled from three homologous subunits with the majority of the protein's mass found in the extracellular domains. Acid-sensing ion channel 1 (ASIC1) is homologous to ENaC, but a key functional domain is highly divergent. Here we present molecular models of the extracellular region of α ENaC based on a large data set of mutations that attenuate inhibitory peptide binding in combination with comparative modeling based on the resolved structure of ASIC1. The models successfully rationalized the data from the peptide binding screen. We engineered new mutants that had not been tested based on the models and successfully predict sites where mutations affected peptide binding. Thus, we were able to confirm the overall general fold of our structural models. Further analysis suggested that the α subunit-derived inhibitory peptide affects channel gating by constraining motions within two major domains in the extracellular region, the thumb and finger domains.     

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K_v1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K_v1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) α-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K_v open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K_v1.5 channel is unable to close when bupivacaine is bound.     

A KcsA/MloK1 Chimeric Ion Channel Has Lipid-dependent Ligand-binding Energetics

Cyclic nucleotide-modulated ion channels play crucial roles in signal transduction in eukaryotes. The molecular mechanism by which ligand binding leads to channel opening remains poorly understood, due in part to the lack of a robust method for preparing sufficient amounts of purified, stable protein required for structural and biochemical characterization. To overcome this limitation, we designed a stable, highly expressed chimeric ion channel consisting of the transmembrane domains of the well characterized potassium channel KcsA and the cyclic nucleotide-binding domains of the prokaryotic cyclic nucleotide-modulated channel MloK1. This chimera demonstrates KcsA-like pH-sensitive activity which is modulated by cAMP, reminiscent of the dual modulation in hyperpolarization-activated and cyclic nucleotide-gated channels that display voltage-dependent activity that is also modulated by cAMP. Using this chimeric construct, we were able to measure for the first time the binding thermodynamics of cAMP to an intact cyclic nucleotide-modulated ion channel using isothermal titration calorimetry. The energetics of ligand binding to channels reconstituted in lipid bilayers are substantially different from those observed in detergent micelles, suggesting that the conformation of the chimera''s transmembrane domain is sensitive to its (lipid or lipid-mimetic) environment and that ligand binding induces conformational changes in the transmembrane domain. Nevertheless, because cAMP on its own does not activate these chimeric channels, cAMP binding likely has a smaller energetic contribution to gating than proton binding suggesting that there is only a small difference in cAMP binding energy between the open and closed states of the channel.     

Localization of the PIP2 sensor of TRPV1 ion channels

Although a large number of ion channels are now believed to be regulated by phosphoinositides, particularly phosphoinositide 4,5-bisphosphate (PIP2), the mechanisms involved in phosphoinositide regulation are unclear. For the TRP superfamily of ion channels, the role and mechanism of PIP2 modulation has been especially difficult to resolve. Outstanding questions include: is PIP2 the endogenous regulatory lipid; does PIP2 potentiate all TRPs or are some TRPs inhibited by PIP2; where does PIP2 interact with TRP channels; and is the mechanism of modulation conserved among disparate subfamilies? We first addressed whether the PIP2 sensor resides within the primary sequence of the channel itself, or, as recently proposed, within an accessory integral membrane protein called Pirt. Here we show that Pirt does not alter the phosphoinositide sensitivity of TRPV1 in HEK-293 cells, that there is no FRET between TRPV1 and Pirt, and that dissociated dorsal root ganglion neurons from Pirt knock-out mice have an apparent affinity for PIP2 indistinguishable from that of their wild-type littermates. We followed by focusing on the role of the C terminus of TRPV1 in sensing PIP2. Here, we show that the distal C-terminal region is not required for PIP2 regulation, as PIP2 activation remains intact in channels in which the distal C-terminal has been truncated. Furthermore, we used a novel in vitro binding assay to demonstrate that the proximal C-terminal region of TRPV1 is sufficient for PIP2 binding. Together, our data suggest that the proximal C-terminal region of TRPV1 can interact directly with PIP2 and may play a key role in PIP2 regulation of the channel.     

Channel Gating Regulation by the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) First Cytosolic Loop

In this study, we present data indicating a robust and specific domain interaction between the cystic fibrosis transmembrane conductance regulator (CFTR) first cytosolic loop (CL1) and nucleotide binding domain 1 (NBD1) that allows ion transport to proceed in a regulated fashion. We used co-precipitation and ELISA to establish the molecular contact and showed that binding kinetics were not altered by the common clinical mutation F508del. Both intrinsic ATPase activity and CFTR channel gating were inhibited severely by CL1 peptide, suggesting that NBD1/CL1 binding is a crucial requirement for ATP hydrolysis and channel function. In addition to cystic fibrosis, CFTR dysregulation has been implicated in the pathogenesis of prevalent diseases such as chronic obstructive pulmonary disease, acquired rhinosinusitis, pancreatitis, and lethal secretory diarrhea (e.g. cholera). On the basis of clinical relevance of the CFTR as a therapeutic target, a cell-free drug screen was established to identify modulators of NBD1/CL1 channel activity independent of F508del CFTR and pharmacologic rescue. Our findings support a targetable mechanism of CFTR regulation in which conformational changes in the NBDs cause reorientation of transmembrane domains via interactions with CL1 and result in channel gating.     

Demonstration of a Direct Interaction between σ-1 Receptors and Acid-Sensing Ion Channels

The σ-1 receptor is a widely expressed protein that interacts with a variety of ion channels, including the acid-sensing ion channel (ASIC) 1a. Here we used atomic force microscopy to determine the architecture of the ASIC1a/σ-1 receptor complex. When isolated His₈-tagged ASIC1a was imaged in complex with anti-His₆ antibodies, the angle between pairs of bound antibodies was 135°, consistent with the known trimeric structure of the channel. When ASIC1a was coexpressed with FLAG/His₆-tagged σ-1 receptor, ASIC1a became decorated with small particles, and pairs of these particles bound at an angle of 131°. When these complexes were incubated with anti-FLAG antibodies, pairs of antibodies bound at an angle of 134°, confirming that the small particles were σ-1 receptors. Of interest, we found that the σ-1 receptor ligand haloperidol caused an ∼50% reduction in ASIC1a/σ-receptor binding, suggesting a way in which σ-1 ligands might modulate channel properties. For the first time, to our knowledge, we have resolved the structure of a complex between the σ-1 receptor and a target ion channel, and demonstrated that the stoichiometry of the interaction is 1 σ-1 receptor/1 ASIC1a subunit.     

Human PIEZO1: Removing Inactivation

PIEZO1 is an inactivating eukaryotic cation-selective mechanosensitive ion channel. Two sites have been located in the channel that when individually mutated lead to xerocytotic anemia by slowing inactivation. By introducing mutations at two sites, one associated with xerocytosis and the other artificial, we were able to remove inactivation. The double mutant (DhPIEZO1) has a substitution of arginine for methionine (M2225R) and lysine for arginine (R2456K). The loss of inactivation was accompanied by ∼30-mmHg shift of the activation curve to lower pressures and slower rates of deactivation. The slope sensitivity of gating was the same for wild-type and mutants, indicating that the dimensional changes between the closed and open state are unaffected by the mutations. The unitary channel conductance was unchanged by mutations, so these sites are not associated with pore. DhPIEZO1 was reversibly inhibited by the peptide GsMTx4 that acted as a gating modifier. The channel kinetics were solved using complex stimulus waveforms and the data fit to a three-state loop in detailed balance. The reaction had two pressure-dependent rates, closed to open and inactivated to closed. Pressure sensitivity of the opening rate with no sensitivity of the closing rate means that the energy barrier between them is located near the open state. Mutant cycle analysis of inactivation showed that the two sites interacted strongly, even though they are postulated to be on opposite sides of the membrane.     

Human PIEZO1: Removing Inactivation

PIEZO1 is an inactivating eukaryotic cation-selective mechanosensitive ion channel. Two sites have been located in the channel that when individually mutated lead to xerocytotic anemia by slowing inactivation. By introducing mutations at two sites, one associated with xerocytosis and the other artificial, we were able to remove inactivation. The double mutant (DhPIEZO1) has a substitution of arginine for methionine (M2225R) and lysine for arginine (R2456K). The loss of inactivation was accompanied by ∼30-mmHg shift of the activation curve to lower pressures and slower rates of deactivation. The slope sensitivity of gating was the same for wild-type and mutants, indicating that the dimensional changes between the closed and open state are unaffected by the mutations. The unitary channel conductance was unchanged by mutations, so these sites are not associated with pore. DhPIEZO1 was reversibly inhibited by the peptide GsMTx4 that acted as a gating modifier. The channel kinetics were solved using complex stimulus waveforms and the data fit to a three-state loop in detailed balance. The reaction had two pressure-dependent rates, closed to open and inactivated to closed. Pressure sensitivity of the opening rate with no sensitivity of the closing rate means that the energy barrier between them is located near the open state. Mutant cycle analysis of inactivation showed that the two sites interacted strongly, even though they are postulated to be on opposite sides of the membrane.     

Lysophosphatidic acid directly activates TRPV1 through a C-terminal binding site

Since 1992, there has been growing evidence that the bioactive phospholipid lysophosphatidic acid (LPA), whose amounts are increased upon tissue injury, activates primary nociceptors resulting in neuropathic pain. The TRPV1 ion channel is expressed in primary afferent nociceptors and is activated by physical and chemical stimuli. Here we show that in control mice LPA produces acute pain-like behaviors, which are substantially reduced in Trpv1-null animals. Our data also demonstrate that LPA activates TRPV1 through a unique mechanism that is independent of G protein-coupled receptors, contrary to what has been widely shown for other ion channels, by directly interacting with the C terminus of the channel. We conclude that TRPV1 is a direct molecular target of the pain-producing molecule LPA and that this constitutes, to our knowledge, the first example of LPA binding directly to an ion channel to acutely regulate its function.     

Kinetics of Ca2+ binding to the SR Ca-ATPase in the E1 state

The time-resolved kinetics of Ca2+ binding to the SR Ca-ATPase in the E1 state was investigated by Ca(2+)-concentration jump experiments. Ca2+ was released by an ultraviolet-light flash from caged calcium, and charge movements in the membrane domain of the ion pumps were detected by the fluorescent styryl dye 2BITC. The partial reaction (H3E1 <-->) E1 <--> CaE1 <--> Ca2E1 can be characterized by two time constants, tau1 and tau2, both of which are not significantly Ca(2+)-concentration-dependent and only weakly pH-dependent at pH < 7.5. Both time constants differ by a factor of approximately 50 (4.7 vs. 200 ms). The weak substrate-dependence indicates that the rate-limiting process is not related to Ca2+ migration through the access channel and ion binding to the binding sites but to conformational rearrangements preceding the ion movements. The high activation energy obtained for both processes, 42.3 kJ mol(-1) and 60.3 kJ mol(-1) at pH 7.2, support this concept. Transient binding of Ca ions to the loop L67 and a movement of the Ca-loaded loop are discussed as a mechanism that facilitates the entrance of both Ca ions into the access channel to the ion-binding sites.     

Nonindependent K+ Movement through the Pore in IRK1 Potassium Channels

We measured unidirectional K⁺ in- and efflux through an inward rectifier K channel (IRK1) expressed in Xenopus oocytes. The ratio of these unidirectional fluxes differed significantly from expectations based on independent ion movement. In an extracellular solution with a K⁺ concentration of 25 mM, the data were described by a Ussing flux-ratio exponent, n′, of ∼2.2 and was constant over a voltage range from −50 to −25 mV. This result indicates that the pore of IRK1 channels may be simultaneously occupied by at least three ions. The IRK1 n′ value of 2.2 is significantly smaller than the value of 3.5 obtained for Shaker K channels under identical conditions. To determine if other permeation properties that reflect multi-ion behavior differed between these two channel types, we measured the conductance (at 0 mV) of single IRK1 channels as a function of symmetrical K⁺ concentration. The conductance could be fit by a saturating hyperbola with a half-saturation K⁺ activity of 40 mM, substantially less than the reported value of 300 mM for Shaker K channels. We investigated the ability of simple permeation models based on absolute reaction rate theory to simulate IRK1 current–voltage, conductance, and flux-ratio data. Certain classes of four-barrier, three-site permeation models are inconsistent with the data, but models with high lateral barriers and a deep central well were able to account for the flux-ratio and single channel data. We conclude that while the pore in IRK1 and Shaker channels share important similarities, including K⁺ selectivity and multi-ion occupancy, they differ in other properties, including the sensitivity of pore conductance to K⁺ concentration, and may differ in the number of K⁺ ions that can simultaneously occupy the pore: IRK1 channels may contain three ions, but the pore in Shaker channels can accommodate four or more ions.     

Atomistic Detailed Mechanism and Weak Cation-Conducting Activity of HIV-1 Vpu Revealed by Free Energy Calculations

The viral protein U (Vpu) encoded by HIV-1 has been shown to assist in the detachment of virion particles from infected cells. Vpu forms cation-specific ion channels in host cells, and has been proposed as a potential drug target. An understanding of the mechanism of ion transport through Vpu is desirable, but remains limited because of the unavailability of an experimental structure of the channel. Using a structure of the pentameric form of Vpu – modeled and validated based on available experimental data – umbrella sampling molecular dynamics simulations (cumulative simulation time of more than 0.4 µs) were employed to elucidate the energetics and the molecular mechanism of ion transport in Vpu. Free energy profiles corresponding to the permeation of Na⁺ and K⁺ were found to be similar to each other indicating lack of ion selection, consistent with previous experimental studies. The Ser23 residue is shown to enhance ion transport via two mechanisms: creating a weak binding site, and increasing the effective hydrophilic length of the channel, both of which have previously been hypothesized in experiments. A two-dimensional free energy landscape has been computed to model multiple ion permeation, based on which a mechanism for ion conduction is proposed. It is shown that only one ion can pass through the channel at a time. This, along with a stretch of hydrophobic residues in the transmembrane domain of Vpu, explains the slow kinetics of ion conduction. The results are consistent with previous conductance studies that showed Vpu to be a weakly conducting ion channel.     

KCNE1 binds to the KCNQ1 pore to regulate potassium channel activity

Potassium channels control the resting membrane potential and excitability of biological tissues. Many voltage-gated potassium channels are controlled through interactions with accessory subunits of the KCNE family through mechanisms still not known. Gating of mammalian channel KCNQ1 is dramatically regulated by KCNE subunits. We have found that multiple segments of the channel pore structure bind to the accessory protein KCNE1. The sites that confer KCNE1 binding are necessary for the functional interaction, and all sites must be present in the channel together for proper regulation by the accessory subunit. Specific gating control is localized to a single site of interaction between the ion channel and accessory subunit. Thus, direct physical interaction with the ion channel pore is the basis of KCNE1 regulation of K+ channels.     

Novel Mechanisms of Trafficking Defect Caused by KCNQ1 Mutations Found in Long QT Syndrome

Long QT syndrome (LQTS) is a hereditary arrhythmia caused by mutations in genes for cardiac ion channels, including a potassium channel, KvLQT1. Inheritance of LQTS is usually autosomal-dominant, but autosomal-recessive inheritance can be observed in patients with LQTS accompanied by hearing loss. In this study, we investigated the functional alterations caused by KCNQ1 mutations, a deletion (delV595) and a frameshift (P631fs/19), which were identified in compound heterozygous state in two patients with autosomal-recessive LQTS not accompanied by hearing loss. Functional analyses showed that both mutations impaired cell surface expression due to trafficking defects. The mutations severely affected outward potassium currents without apparent dominant negative effects. It was found that delV595 impaired subunit binding, whereas P631fs/19 was retained in endoplasmic reticulum due to the newly added 19-amino acid sequence containing two retention motifs (R⁶³³GR and R⁶⁴⁶LR). This is the first report of novel mechanisms for trafficking abnormality of cardiac ion channels, providing us new insights into the molecular mechanisms of LQTS.     

Coupling between Residues on S4 and S1 Defines the Voltage-Sensor Resting Conformation in NaChBac

The voltage sensor is a four-transmembrane helix bundle (S1-S4) that couples changes in membrane potential to conformational alterations in voltage-gated ion channels leading to pore opening and ion conductance. Although the structure of the voltage sensor in activated potassium channels is available, the conformation of the voltage sensor at rest is still obscure, limiting our understanding of the voltage-sensing mechanism. By employing a heterologously expressed Bacillus halodurans sodium channel (NaChBac), we defined constraints that affect the positioning and depolarization-induced outward motion of the S4 segment. We compared macroscopic currents mediated by NaChBac and mutants in which E43 on the S1 segment and the two outermost arginines (R1 and R2) on S4 were substituted. Neutralization of the negatively charged E43 (E43C) had a significant effect on channel gating. A double-mutant cycle analysis of E43 and R1 or R2 suggested changes in pairing during channel activation, implying that the interaction of E43 with R1 stabilizes the voltage sensor in its closed/available state, whereas interaction of E43 with R2 stabilizes the channel open/unavailable state. These constraints on S4 dynamics that define its stepwise movement upon channel activation and positioning at rest are novel, to the best of our knowledge, and compatible with the helical-screw and electrostatic models of S4 motion.     

Effects of activation of skin ion channels TRPM8, TRPV1, and TRPA1 on the immune response. Comparison with effects of cold and heat exposure

The effects of pharmacological stimulation of skin ion channels TRPA1, TRPM8, TRPV1 on the immune response are presented. These effects are compared with the effects of different types of temperature exposures - skin cooling, deep cooling, and deep heating. This analysis allows us to clear the differences in the influence on the immune response of thermosensitive ion channels localized in the skin; (2) whether the changes in the immune response under temperature exposures are due to these thermosensitive ion channels. Experiments were performed on Wistar rats. For stimulation of TRPM8 ion channel, an application to the skin of 1% menthol was used, for TRPA1 - 0.04% allylisotiocianate, and for TRPV1 - capsaicin in a concentration of 0.001.The antigen binding in the spleen was two-times stimulated by activation of the cold-sensitive ion channel TRPM8 and much weaker by activation of warm-sensitive TRPV1 (by 15%), and another cold-sensitive ion channel TRPA1 (by 40%). Only the stimulation of TRPA1 significantly (by 140%) increased antibody formation in the spleen, while TRPM8 had practically no effect on this process, and activation of TRPV1 significantly (by 60%) inhibited antibody formation. Stimulation of the TRPM8 ion channel significantly (by 60%) reduced the level of IgG in the blood, which is believed to control of infectious diseases.The obtained results show that pharmacological activation of the skin TRPA1, TRPM8, TRPV1 ion channels can differently affect the immune system. At the epicenter of changes there were the antigen binding and antibody formation in the spleen, as well as the level of IgG in the blood. Exactly stimulation of the TRPM8 ion channel determines the changes in the immune response when only the skin is cooling, while at deep body heating, the changes in the immune response are mostly determined by the activation of the skin TRPV1 ion channel.