Structural Determinants of the Transient Receptor Potential 1 (TRPV1) Channel Activation by Phospholipid Analogs*

Background: The TRPV1 ion channel can be regulated by negatively charged lipids. Results: TRPV1 shows specificity for LPA analogs containing monounsaturated hydrocarbon chains with a negatively charged phosphate, cyclic phosphate, and thiophosphate headgroup. Conclusion: TRPV1 activation is highly restricted to natural lipids with oleyl or oleoyl side chains. Significance: Production of endogenous C18:1 lysophospholipids can selectively trigger activation of TRPV1 and nociceptive neuronal responses. The transient receptor potential vanilloid 1 (TRPV1) ion channel is a polymodal protein that responds to various stimuli, including capsaicin (the pungent compound found in chili peppers), extracellular acid, and basic intracellular pH, temperatures close to 42 °C, and several lipids. Lysophosphatidic acid (LPA), an endogenous lipid widely associated with neuropathic pain, is an agonist of the TRPV1 channel found in primary afferent nociceptors and is activated by other noxious stimuli. Agonists or antagonists of lipid and other chemical natures are known to possess specific structural requirements for producing functional effects on their targets. To better understand how LPA and other lipid analogs might interact and affect the function of TRPV1, we set out to determine the structural features of these lipids that result in the activation of TRPV1. By changing the acyl chain length, saturation, and headgroup of these LPA analogs, we established strict requirements for activation of TRPV1. Among the natural LPA analogs, we found that only LPA 18:1, alkylglycerophosphate 18:1, and cyclic phosphatidic acid 18:1, all with a monounsaturated C18 hydrocarbon chain activate TRPV1, whereas polyunsaturated and saturated analogs do not. Thus, TRPV1 shows a more restricted ligand specificity compared with LPA G-protein-coupled receptors. We synthesized fatty alcohol phosphates and thiophosphates and found that many of them with a single double bond in position Δ9, 10, or 11 and Δ9 cyclopropyl group can activate TRPV1 with efficacy similar to capsaicin. Finally, we developed a pharmacophore and proposed a mechanistic model for how these lipids could induce a conformational change that activates TRPV1.

The transient receptor potential vanilloid 1 (TRPV1) ion channel is a polymodal protein that responds to various stimuli, including capsaicin (the pungent compound found in chili peppers), extracellular acid, and basic intracellular pH, temperatures close to 42°C, and several lipids. Lysophosphatidic acid (LPA), an endogenous lipid widely associated with neuropathic pain, is an agonist of the TRPV1 channel found in primary afferent nociceptors and is activated by other noxious stimuli. Agonists or antagonists of lipid and other chemical natures are known to possess specific structural requirements for producing functional effects on their targets. To better understand how LPA and other lipid analogs might interact and affect the function of TRPV1, we set out to determine the structural features of these lipids that result in the activation of TRPV1. By changing the acyl chain length, saturation, and headgroup of these LPA analogs, we established strict requirements for activation of TRPV1. Among the natural LPA analogs, we found that only LPA 18:1, alkylglycerophosphate 18:1, and cyclic phosphatidic acid 18:1, all with a monounsaturated C18 hydrocarbon chain activate TRPV1, whereas polyunsaturated and saturated analogs do not. Thus, TRPV1 shows a more restricted ligand specificity compared with LPA G-protein-coupled receptors. We synthesized fatty alcohol phosphates and thiophosphates and found that many of them with a single double bond in position ⌬9, 10, or 11 and ⌬9 cyclopropyl group can activate TRPV1 with efficacy similar to capsaicin. Finally, we developed a pharmacophore and proposed a mechanistic model for how these lipids could induce a conformational change that activates TRPV1.
The transient receptor potential vanilloid 1 (TRPV1) 3 is a member of the TRP channel family of ion channels that has been widely associated with painful processes such as angina pectoris, arthritis, and cancer, as well as inflammatory events (1)(2)(3)(4)(5). The activation of this receptor in the peripheral sensory cells where it is expressed leads to pain perception (6,7).
This polymodal ion channel is activated by various stimuli that include low extracellular and high intracellular pH, pungent compounds found in some plants, high temperatures (8) and lipids such as lysophosphatidic acid 18:1 (LPA) (9).
LPA is found in serum and plasma, and it can reach up to micromolar concentrations in biological fluids (10,11). LPA has been importantly linked to neuropathic pain, and for several years it was thought to produce it solely through its interaction with specific GPCRs (LPA 1-6 ) (12,13), with the LPA 1 receptor being particularly important in this process (14). Recently, our group demonstrated a direct and physiologically relevant interaction of TRPV1 with LPA. By using a combination of electrophysiological and biochemical techniques, we showed that LPA 18:1 and the bromophosphonate analog of LPA (BrP-LPA), which is an antagonist of several of the LPA receptors, activate TRPV1 through a direct interaction with the proximal C-terminal region of the ion channel, a site that is also occupied by phosphatidylinositol 4,5-bisphosphate (PIP 2 ) (9).
Finally, we were able to observe nocifensive behavior generation when we injected LPA into the paws of mice. This response was diminished in TRPV1 knock-out animals, indicating that LPA can produce acute pain through the activation of TRPV1 (9).
Mechanistically, we found that the neutralization of the positively charged residue Lys-710 in the C terminus of the protein rendered the mutant TRPV1-K710Q channels less sensitive to LPA, whereas when the charge was reversed (K710D), the channels were almost completely insensitive to the lipid, indicating that activation of the channel was partly dependent on an electrostatic interaction (9).
Another lipid that regulates the activity of TRPV1 is PIP 2 (15)(16)(17)(18). The role of PIP 2 has been controversial, with one research group suggesting that PIP 2 is not necessary for recombinant purified TRPV1 activation in artificial liposomes (19). In fact, it has been shown that PIP 2 can inhibit TRPV1 function but under conditions where it is incorporated to the extracellular side of the plasma membrane (18). However, it is accepted that intracellular PIP 2 is necessary for the activation of TRPV1 by capsaicin in biological membranes (20,21). Moreover, similar to what we had previously shown for LPA (9), several negatively charged lipids have been shown to positively regulate the activity of TRPV1 (16). We set out to determine the structural requirements for the activation of TRPV1 by lipids similar to LPA, some of which are naturally occurring lysophospholipids.
For other channel types and lipids, it has been shown that the length of the acyl chain, the presence or absence of unsaturations, as well as the nature of the headgroup charge linked with the acyl chain can influence the effect of these molecules on ion channel function. For example, M-type potassium channels are activated by a large variety of lipid phosphates with minimum requirements for acyl chain and the presence of a phosphate headgroup (22,23).
In this study, we have evaluated the impact of the LPA structural motifs, the hydrocarbon chain, the glycerol backbone, and the phosphate headgroup, on the structure-activity relationship of TRPV1 activation, and we have delineated the minimum requirements for agonist activity in this family of compounds.
Evoked currents from TRPV1-expressing HEK293 cells were low-pass filtered at 2.9 kHz and sampled at 10 kHz with an EPC 10 amplifier (HEKA Elektronik GmbH, Pfalz, Germany). Inside-out patches were held at 0 mV, and voltage was then stepped from 0 to Ϫ120 to 120 mV for 100 ms in 10-mV increments and then returned to 0 mV. All recordings were performed at room temperature (19 Ϯ 3°C). Currents obtained in the presence of the different lipids were normalized to the maximal current obtained with a saturating (4 M) capsaicin concentration (I/I max ). To perform dose-response curves, a given concentration of the lipid was applied to a TRPV1-expressing membrane patch, and currents were measured and normalized to the maximal current obtained in the presence of 4 M capsaicin. A single concentration of lipid was used to obtain the average for activation from several different membrane patches for each point in the dose-response curves, and then the data were fitted to the Hill equation, as described previously (26). Data were acquired and analyzed with PULSE software (HEKA Elektronik) and were plotted and analyzed with programs written using Igor Pro (WaveMetrics Inc., Portland, OR).
In Vitro Interaction Assays-Surface proteins were obtained from HEK293 cells transiently expressing TRPV1 channels using the cell surface isolation kit (Pierce) according to the manufacturer's instructions. For the overlay assay, lipids were dissolved in DMEM (Invitrogen) with 1% fatty acid-free BSA and prepared as 10 mM stocks and spotted (100 nmol/spot) onto a nitrocellulose membrane (GE Healthcare). The membrane was air-dried for 1 h at room temperature and blocked with PBS containing 10% nonfat dry milk and 1% fatty acid-free BSA for 2 h with gentle shaking. After washing twice with PBS, 0.05% Tween 20 (PBS-T), the membrane was incubated with the surface protein solution containing 10 -20 g of total surface protein in 50 mM Tris-HCl (pH 8), 150 mM NaCl, 0.5% IGEPAL (Sigma), 1ϫ complete protease inhibitor mixture (Roche Applied Science), and 0.5% BSA at 4°C for 48 h with shaking. After three washes with PBS-T, the membranes were incubated with anti-TRPV1 antibody (1:500; sc-12498, P-19; Santa Cruz Biotechnology, Santa Cruz, CA) in PBS-T with 5% nonfat dry milk overnight at 4°C. Finally, following three more washes with PBS-T, the membranes were incubated with horseradish peroxidase-conjugated secondary anti-goat antibody, and the binding of lipids to TRPV1 was visualized by chemiluminescence by exposing the blot for 15 min (Amersham Biosciences). For assays with LPA-coated beads, plasma membrane protein extracts were incubated with LPA beads or with control beads (Echelon Biosciences) for 2 h at 4°C and washed five times with binding buffer containing 0.5% IGEPAL (Sigma). For competition assays, protein extracts were previously incubated with 10 M FAP-4 lipid for 2 h. Pro-teins pulled down from each sample were eluted from the beads by adding an equal volume of 2ϫ Laemmli sample buffer to the beads and heating to 95°C for 5 min. Eluted proteins were separated by SDS-PAGE and analyzed by immunoblot assays using an anti-TRPV1 (sc-12498, P-19; Santa Cruz Biotechnology) antibody. In every experiment, a small sample of the total membrane protein extract was loaded to the gel as a positive control for TRPV1 protein expression (lane 1 in Fig. 3C). Immunoblots were developed by ECL reagent, and densitometric analysis was done using ImageJ software (National Institutes of Health) and expressed as relative protein levels versus the amount of TRPV1 bound to LPA beads (9). Densitometric analysis for the overlay assays was performed also by using the ImageJ software, although these assays are not strictly quantitative.
Ligand-based TRPV1 Pharmacophore Modeling-The flexible alignment models were created with the Molecular Operating Environment (MOE) 2011.10 from Chemical Computing Group (Montreal, Canada). Active compounds were aligned with the flexible alignment module implemented in MOE using default parameters. The alignment with the lowest alignment score S and average strain energy score U were selected. The active alignment model was then fixed, and inactive compounds were separately aligned with the flexible alignment module using the same settings. The final conformation for each compound was selected on the basis of the compound's lowest energy strain and overlap of molecular features.
Pharmacophore modeling was performed using the pharmacophore query editor function in MOE. Models were manually created using the lowest energy conformation of the flexibly aligned models. The pharmacophore points (one anionic group and two hydrophobic groups) were identified as those structural features of the active compounds also sharing a common volume.
Molecular Modeling-Monte Carlo conformational analysis of lipids was carried out using the Multiple Minimum program in the MacroModel suite of software Version 9.7 (Schrödinger Inc., Surrey, UK). Coordinates of LPA 18:0, 18:1, and 18:2 were obtained from sdf files found at the Pub-Chem database (CID: 9547179, 5311263, and 53478601). A 10,000-step conformational search was performed on the three molecules using the OPLS-2005 force field and the water solvation model (27). Conformers were selected using a 3.0 kcal/mol energy cutoff and pooled into geometrically similar families based on the root mean square difference between corresponding torsion angles in pairs of structures using XCluster.
Modeling of TRPV1-Lipid Interactions-The structures of monomers of the TRPV1 channel in either the closed (PDB code 3J5P-B) or open (PDB code 3J5Q-D) state were used as templates for docking the most abundant structural conformations of 18-carbon LPAs. Residues 361-719 were used in both cases because the binding of PIP 2 (17,18,21) and LPA (9) has been found to occur on the proximal C terminus, which was also supported by molecular modeling experiments (28). Structures were prepared in Autodock Tools, and all the dockings were carried out with Autodock 4.2 (29) with a Lamarckian  (30).
Statistical Analysis-Statistical comparisons were made with a one-way analysis of variance test. p Ͻ 0.05 was considered statistically significant. Group data are reported as the mean Ϯ S.E.

TRPV1 Channel Is Activated by Lipids Similar to LPA-Our
approach to a comprehensive study of the structural requirements of the lipid ligands necessary to activate TRPV1 was based on the designation of five motifs in the LPA scaffold as follows: 1) the glycerol backbone and its substitutions; 2) the  To begin to determine the nature of the requirements for activation of TRPV1 by LPA, we investigated what effect the lack of a double bond or the presence of more than one double bond in the hydrocarbon chain would have on TRPV1 activation. For this purpose, we used commercially available saturated LPA 18:0 (5 M) and LPA 18:2 (5 M). To our surprise, although structurally similar to LPA 18:1, both LPA 18:0 and LPA 18:2 failed to activate the channel (Fig. 1, B and C, red traces, and E). Moreover, FAP-7, a synthetic analog with an 18:2 fatty alcohol chain also failed to activate TRPV1 (Fig. 1, D, red trace, and E), consistent with the result we obtained with LPA 18:2. Additionally, both LPA 18:0 and LPA 18:2 failed to activate TRPV1 even at concentrations as high as 25 M (Table 1), the highest concentration we could use without breaking the membrane patches. These data suggested that the absence of an unsaturation or the presence of more than one unsaturation leads to a lack of activation of TRPV1 by lipids with these structural features.
Using an overlay assay, we detected an interaction of LPA 18:1, LPA 18:2, and the fatty alcohol thiophosphate FAP-7 can also interact with TRPV1 (Fig. 1, F and G). However, this inter-action does not ensure changes in the function of TRPV1, supporting the idea that lipid-protein interactions that lead to activation of the channel depend tightly on the structural features of the lipid.
We then tested another commercially available saturated lipid, CPA 16:0, and two naturally occurring unsaturated lipids, CPA 18:1 that is found in serum (31) and dioleoylglycerophosphate (DGPP) that is found in plants and yeast (32)(33)(34). When CPA 16:0 (5 M) was applied to inside-out excised membrane patches from HEK293 cells expressing TRPV1, it did not produce TRPV1 current activation (Fig. 2, A and D). However, CPA 18:1 (5 M) produced currents of a similar magnitude at a saturating capsaicin concentration (4 M) (Fig. 2, B and D), and the dose-response to this lipid yielded a K1 ⁄ 2 of 2.7 M. However, when DGPP 18:1 with two oleyl chains was applied to these TRPV1-expressing membrane patches, no activation was observed (Fig. 2, C and D). Moreover, using an overlay assay, we also determined an interaction of the activator CPA 18:1 with TRPV1 (Fig. 2, E and F). The overlay assays also showed that DGPP 18:1 interacted with TRPV1. However, even though DGPP 18:1 bound to TRPV1, this interaction did not result in channel activation.
Thus, the data shown in Figs. 1 and 2 demonstrated that unsaturated mono-oleoylglycerophosphates were capable of activating TRPV1, whereas a dioleoylglycerophosphate was ineffective in promoting the opening of TRPV1.   Synthetic Monounsaturated Fatty Alcohol Phosphates Activate TRPV1-To evaluate the role of the glycerol backbone in TRPV1 activation, we synthesized 16 fatty alcohol phosphate and thiophosphate compounds (FAPs) ( Table 1) that shared long hydrocarbon chains (18 -20 carbons) with one double bond located at position 9, 10, or 11. Fig. 3 shows that FAP-3, FAP-4, FAP-5, and FAP-6 all activated TRPV1 with an efficiency similar to that of 4 M capsaicin as measured by I lipid / I max capsaicin (Fig. 3B). Interestingly, compound FAP-8, which presents a cyclopropyl group at position 9, also activated TRPV1 with similar efficiency as 4 M capsaicin (Fig. 3B). FAP-10 with a branched hydrocarbon side chain was inactive (data not shown; see Table 1 for structure). These results point to a possible role of a region that can provide a certain degree of torsion to the molecule (i.e. unsaturation and a cyclopropyl group) (35) in the specificity for activation of TRPV1 by these lipids.
We had previously shown that LPA interacts with the Lys-710 residue in the C terminus of TRPV1 (9). We used the K710D mutant TRPV1 channel to determine whether these compounds interact with the same site in the C terminus as LPA 18:1 does to activate the channel. As shown in Fig. 3B, all of the compounds tested showed a clear decrease in the activation of the mutant TRPV1 channel as compared with the WT TRPV1 channel.
Moreover, as demonstrated previously (9), we found that the TRPV1 channel interacts with LPA-coated beads in pulldown assays (Fig. 3C, lane 2). In addition, when TRPV1 was preincubated with one of the fatty alcohol phosphates described above, FAP-4, the amount of channel pulled down by the LPA-coated beads was greatly diminished (Fig. 3C, lane 3), also indicating that lipid FAP-4 competes for the same site that LPA 18:1 binds to in the proximal C terminus of TRPV1.
Lipid Analogs with Dodecanoyl or Shorter Hydrocarbon Chains Fail to Activate TRPV1-We have previously shown that LPA binds to the proximal region of the C terminus of TRPV1 through a mechanism that partly involves electrostatic interactions, presumably with the charged phosphate group of LPA and the charged Lys-710 residue of TRPV1. Under this scenario, the hydrophobic acyl chain of LPA would be inserted in the membrane or buried in the membrane-protein interface with its polar headgroup facing the cytoplasmic side interacting with Lys-710. Hence, we hypothesized that changes in the length of the acyl chain of LPA may curtail the effect of this lipid on TRPV1 activity.
To test the latter possibility, we purchased a short-chain LPA and synthesized short-chain FAPs that ranged from aliphatic chain lengths of 6 -12 carbons. In this case, neither LPA 6:0, FAP-15, FAP-16, nor FAP-17 (see Table 1 and Fig. 4) produced activation of TRPV1 currents (Fig. 4, A-D, red traces, and E) as compared with activation by capsaicin. Moreover, we tested whether these lipids could bind to TRPV1, and the representative biochemical assay shown in Fig. 4, F and G, is indicative of either a weak or a complete lack of interaction of these shortchain lipids with the channel. However, the aforementioned lipids are not only molecules with short chains but these chains are also saturated, a trait that was also distinctive from those lipids that activate TRPV1 as shown in Fig. 1.
Because double bonds provide a considerable degree of torsion and motion to the fatty acid chains (35), which may be important for proper physical interaction with a site in the ion channel, we synthesized the ⌬4-monounsaturated hexyl chain FAP-14 and tested whether it could activate TRPV1 by applying this lipid to the intracellular face of the channel. The results in Fig. 4, D (red trace) and E, show that FAP-14 failed to induce TRPV1-mediated currents, as the currents in the presence of this lipid were indistinguishable from those obtained initially in the absence of chemical ligands. Moreover, we performed overlay assays to determine whether there was an interaction of FAP-14 with TRPV1 and found that in fact it did bind to the channel (Fig. 4, F and G). Thus, these data indicated that the presence of an unsaturation at position ⌬4 in a short-chain lipid did not render it functionally relevant as an agonist of TRPV1, although it could interact with the protein.
Other Structural Requirements for the Activation of TRPV1 by Long-chain Lipids-To this point, using lipids similar to LPA, we have shown that a long chain and the presence of one unsaturation or a torsion-providing cyclopropyl group are necessary for the activation of TRPV1. We now sought to determine whether lipids with an ether link instead of an ester link between the hydrocarbon chain and the glycerol backbone could also promote TRPV1 activation. For this purpose, we also used two commercially available alkyl lipids (AGP 18:1 and AGP 18:0), and we synthesized another two similar ether- linked compounds, FAP-11 and FAP-12, that differed from each other in the position of the ether with respect to the thiophosphate headgroup. When TRPV1-expressing membrane patches from HEK293 cells were exposed to AGP 18:1 (Fig. 5, A and E), we found that the channel activated similarly to what we observed with LPA 18:1 (Fig. 1A). In contrast, TRPV1 currents could not be activated by the application of AGP 18:0 (Fig. 5, B  and E), reinforcing the hypothesis that the presence of an unsaturation is required for activation of the channel. Moreover, we once again found that a lipid such as AGP 18:0 that does not activate the channel can bind to TRPV1 (Fig. 5F).
Interestingly, although the FAP-12 compound clearly activated TRPV1 (Fig. 5, C and E), the FAP-11 compound, at a concentration of 5 M, failed to induce TRPV1-mediated currents (Fig. 5, D and E), indicating that the position of the ether bond with respect to the negatively charged headgroup is also important. However, when we tested FAP-11 at a concentration of 25 M, we were able to observe nearly 80% of activation of the TRPV1 channel as compared with a saturating concentration of capsaicin (Table 1), indicating that although less effective than FAP-12, FAP-11 can induce channel activation at high concentrations.
Then we decided to test whether the glycerol phosphate (which contains the charged headgroup; Fig. 6A) oleic acid (an unsaturated 18-carbon acyl chain without a charged head-group; Fig. 6B) or whether 1-O-hexadecyl glycerol (a saturated 16-carbon alkyl chain ether without a charged headgroup; Fig.  6C, HG) could influence TRPV1 activity on their own. In summary, we found that none of these compounds could promote TRPV1 function (Fig. 6, A-C, red traces, and F), demonstrating that these different component parts of lipids alone are not structurally sufficient for activation of TRPV1.
To further probe the role of the negatively charged headgroup, we synthesized the fatty alcohol analog FAP-9 with a boron phosphate and a hydrocarbon chain of 18 carbons with monounsaturation at the ⌬9 position. This compound was ineffective in activating the channel (Fig. 6, D, red trace, and F), suggesting that a variation in the electronegativity of the charged headgroup is responsible for this observation. Moreover, a compound similar to FAP-9 but without the boron phosphate group, FAP-1, was able to activate TRPV1 (Fig. 6, E  and F).
Finally, we had previously shown that the LPA G-proteincoupled receptor (GPCR) pan-antagonist BrP-LPA (9, 36) is an activator of TRPV1, despite the fact that it has a saturated acyl chain. This observation is contrary to the findings we presented above for a host of compounds demonstrating a strict requirement for a single double bond in the hydrocarbon chain of active ligands. To resolve this apparent contradiction between our present results showing that ligands only with a single double bond are effective activators of TRPV1, all active compounds were superimposed using the "flexible alignment" procedure implemented in MOE (MOE 2011.08 version, Chemical Computing Group Inc., Montreal, Quebec, Canada). A ligand-based pharmacophore consisting of an anionic headgroup and two hydrophobic groups (Fig. 7, red and green dots) highlights key features shared by all of the active compounds. All active compounds were able to adopt a conformation that the inactive compounds do not. In addition to the negatively charged headgroup, the two hydrophobic points where the active compounds intersect appear to be critical for activity. In Fig. 7A, an overlay of all active compounds is shown. As an example, LPA 18:1 (Fig. 7B) shares all three features defined by the overlay of all active ligands in Fig. 7A. In contrast, in the case of the inactive LPA 18:0, whose superimposition with LPA 18:1 is shown in Fig. 7C, the hydrocarbon chain bends away from the middle point of the pharmacophore. Surprisingly, BrP-LPA (Fig. 7D), although saturated, is predicted to cross the middle hydrophobic feature similarly to the unsaturated active compounds. The bromophosphonate headgroup appears to push its hydrophobic feature in the middle closer to the end of the hydrocarbon tail. However, BrP-LPA curls from left to right instead of the predicted right to left curl in LPA 18:1 (Fig. 7, B and D). All of the above support the hypothesis that the hydrophobic feature in the middle of the pharmacophore and the nature of the negatively charged headgroup together play an important role in defining activity.

DISCUSSION
We have shown previously that the TRPV1 channel is activated by LPA 18:1. This was the first description of a direct interaction of LPA, a molecule implicated in neuropathic pain, with a TRP channel that importantly mediates noxious signals (9).
In light of these previous findings and other recent studies showing that negatively charged lipids can positively regulate the activity of TRPV1 channels (16 -18, 20), we sought to determine the structural requirements that lipids and compounds similar to LPA have to meet to elicit TRPV1 activation. For this purpose, we used a variety of naturally occurring and synthetic analogs to discern their requirements for TRPV1 activation.
By varying the chain lengths, unsaturations, their positions, and their headgroups, we were able to identify several lipids that activate TRPV1 (AGP 18:1, CPA 18:1, and FAP-1-9) in a similar fashion to LPA 18:1 (see Table 1). Moreover, we tested (37,38) whether the endogenously produced lipid CPA 18:1 could bind to the protein in biochemical assays, and we found that it does bind to the TRPV1 channel. The effects of long acyl chain lipids on ion channel regulation is not unique for the TRPV1 channels (16), but it has also been described for K ATP ion channels (39), large conductance Ca 2ϩ potassium (BK) channels (40), and inwardly rectifying potassium (Kir) channels (41). However, short-chain LPAs, saturated or monounsaturated, did not activate TRPV1. Based on these observations, we hypothesized that the chain length might be important because the acyl chain could interact directly with the core of the protein, whereas the phosphate group interacts with charged residues in the association domain, including Lys-710, a residue previously shown to form part of a binding site for LPA (9).
Our results using AGP and the FAP series indicate that the chemical bond linking the headgroup to the hydrocarbon chain can be either an ester or an ether bond, and TRPV1 activation is still achieved. We also determined that the glycerol backbone is not essential for TRPV1 activity. The negatively charged headgroup can be a phosphate, thiophosphate, or cyclic phosphate (Table 1). However, the boron phosphate does not elicit TRPV1 activation. Single double bonds in positions ⌬9, 10 or 11 of the hydrocarbon chain, in either cis or trans configuration or a cyclopropyl group in the ⌬9 position, were all tolerated and active. The ligand-based pharmacophore identified three features shared by all active compounds. Our data support the ionic interaction between the negatively charged headgroup with Lys-710 in the proximal C terminus, as shown by the decreased effect of the FAP compounds with the K710D TRPV1 FIGURE 7. Ligand-based pharmacophore of TRPV1 activation reveals three key features. A, three-dimensional alignment of the active compounds (yellow). Three-point pharmacophore query modeled on the bioactive conformation of active compounds: the anionic group with two hydrophobic features. B, LPA 18:1 shown as representative of the ability of unsaturated compounds to adopt an active conformation. C, superimposition of LPA 18:1 (yellow) and LPA 18:0 (stick) reveals that LPA 18:1 occupies the three pharmacophore features, whereas LPA 18:0 is unable to extend to the critical hydrophobic intersection. D, flexible alignment of LPA 18:1 (yellow) and BrP-LPA (stick) indicates that a similar low energy conformation can be adopted, although the two compounds curl in opposite ways. mutant channel. Nonetheless, it is possible that they have a complex binding site, and they also interact with other residues in the ion channel because some of these compounds displayed more current activation of the TRPV1-K710D mutant than what we had previously reported for LPA 18:1 (9). Moreover, we do not know which part of the channel protein interacts with the hydrocarbon chain. Identification of this interaction site awaits further investigation. We also find that the structural determinants for activation of the TRPV1 channel by the lipids studied here are far more restricted for LPA species than that of the LPA GPCRs. LPA GPCRs are activated by several analogs that were ineffective at activating TRPV1, including shortchain LPAs with C10 -12 carbons and saturated or polyunsaturated LPA species with chain lengths Ն18.
In the acetylcholine receptor, it has been shown that the position of the double bond in monounsaturated free fatty acids is pivotal for the inhibitory effects of these lipids on the receptor channel (42). By studying the effects of having a double bond at different positions, Perillo et al. (42) determined that only those in positions 6 and 9 could exert changes in the open probability of the receptor channel. From these results, it was concluded that the position and isomerism of the torsion angle of unsaturated free fatty acids are important for their role as inhibitors of the acetylcholine receptor.
Here, we have tested long-chain molecules with none, one, or two double bonds on the activation of TRPV1. We find that neither LPA 18:0 nor LPA 18:2 (not even at concentrations as high as 25 M) could activate the channel. Previous studies have shown that, for the case of C18 N-acylethanolamines, the presence of two unsaturations renders these molecules more effective as agonists of TRPV1 (27). The unsaturation of LPA 18:1 provides an increase in the rate of motion and a decrease in the order of the acyl chain (35), probably allowing it to gain a conformation suitable for interaction with amino acids in a hydrophobic pocket in the transmembrane region of the channel. However, because DGPP 18:1 could not activate TRPV1, it is possible that the hydrophobic pocket with which LPA 18:1 interacts is unable to accommodate two fatty acid chains.
Lipids that possess different degrees of unsaturation acquire different structures. As can be seen from the modeling of the lipids in Fig. 8A, the LPA 18:0 structure is linear, whereas the LPA 18:1 and LPA 18:2 molecules acquire different degrees of torsion. It is also possible that the presence of an unsaturation in the aliphatic chain of a lipid not only provides the right degree of torsion but also provides the ability to establish a -electron bond with an aromatic residue (43) located in the transmembrane region of the channel, stabilizing the lipid-protein interaction.
It has been suggested that the polyunsaturated LPA species 18:2 and 20:4 and the monounsaturated 18:1 constitute the majority of the LPA in serum with 38, 39, and 9%, respectively (44). Using an animal model where nerve injury was experimentally induced, it was previously determined that the levels of 18:1, 16:0, and 18:0 LPA in the spinal dorsal horn are increased (45). In this model, it was shown that LPA 18:1 is the major species of LPA in quantity and the most relevant species in pain generation (45).
For other channels, lipid-protein interactions have been proposed. For example, for the Kir 2.2 channel there is an x-ray crystal structure of the protein in the presence of PIP 2 . Hansen et al. (46) found that PIP 2 binds at an interface between the transmembrane domain and a nonspecific phospholipid-binding region in the cytoplasmic domain. When PIP 2 binds, there is contraction of a flexible expansion linker that changes this region into a compact helical structure (46). This could be a similar scenario to what happens with LPA and TRPV1. Fig. 8B shows an example of a lipid-protein interaction between LPA 18:0 (that does not activate the channel) and TRPV1. In this case, the hypothesis is that the charged phosphate headgroup of LPA 18:0 is not correctly oriented and thus is incapable of interacting with Lys-710.
In a recently obtained cryostructure of the TRPV1 channel, the Lys-710 residue in the proximal C terminus of the channel (located near the TRP box) moves parallel to the axis of the membrane when the channel is in an activated state (Fig. 8C) with either capsaicin or resiniferatoxin and double-knot toxin as compared with the apo-state (47). Here, we propose that long-chain LPA-like molecules of 18 carbons with one double bond at the pharmacophore feature (Fig. 7), which provides a "kink" to the lipids, adopt a conformation that promotes that the positively charged Lys-710 residue moves parallel to the membrane to be able to interact with the charged headgroup, causing a conformational change that can render TRPV1 active (Fig. 8C, right).
In our biochemical assays, we found that lipids that did not activate TRPV1 could still bind to the channel (e.g. FAP-14 and DGPP 18:1), probably through their interaction with other charged residues in the protein. However, the fact that these lipids could not affect TRPV1 function shows that not all interactions of lipids with the channel lead to a change in its function, similar to what happens with the Kir2.2 channel (48). We interpret these findings as indicating that binding to and activation of the channel represent two distinct events in the lipidprotein interaction of TRPV1.
We propose that the binding of the 18:1 fatty acid of LPA in a hydrophobic pocket of the channel and the electrostatic attraction of Lys-710 produce a distortion/conformational change in the association domain that is transduced to the pore segments and the gating machinery, triggering the channel to open. Our data are in accordance with previous findings that, in excised membrane patches, lipids positively regulate the activity of TRPV1.