A binding site for phosphoinositides described by multiscale simulations explains their modulation of voltage-gated sodium channels
- Research School of Biology, Australian National University, Australia
Figures
Figure 1
Structure of Naᵥ and modulation by phosphoinositides.
(A) Naᵥ channel topology featuring transmembrane helices (S1–S6), the selectivity filter (SF), and the DIII–IV linker (containing the IFM motif) located between DIII and DIV. (B) Naᵥ1.4 structure (6agf) showing the four domain-swapped voltage-sensing domains (VSDs I–IV), pore, and DIII–IV linker on the intracellular side. (C) Summary of PI(4,5)P2 effects on transitions between Naᵥ channel functional states (Gada et al., 2023). (D) Structure of the PI(4,5)P2 headgroup with the 4’- and 5’-phosphates indicated.
Figure 2 with 3 supplements
Lipid fingerprint and binding of all phosphoinositide (PIP) types to Naᵥ1.4.
(A) Naᵥ1.4 embedded in a 360 Å × 360 Å model mammalian membrane containing 63 lipid species. (B) Lipid depletion enrichment index of lipids around Naᵥ1.4 grouped into 12 headgroup classes. (C) Naᵥ1.4, shown from the intracellular (left) and membrane (right) sides, colored by PIP occupancy (darker purple = greater PIP occupancy). (D) Distribution of PIP binding occupancies (left) and occupancy of lipid species at four residues with the highest PIP occupancy. Error bars show standard error, n = 3.
Figure 2—figure supplement 1
Lipid z-density maps (left) and contact occupancy structures (right) for all 12 lipid classes in the mammalian membrane.
Shows lipid species (LPC), glycosphingolipid (GM), ceramide (CER), phosphoinositide (PIP), phosphatidic acid (PA), and phosphotidylinositol (PI). Density maps for each lipid type are show on the right. Contact occupancy structures are visualized from the intracellular view, the voltage-sensing domain (VSD)-IV side view, the extracellular view and the VSD-II side view (left to right).
Figure 2—figure supplement 2
Lipid z-density maps (left) and contact occupancy structures (right) for all 12 lipid classes in the mammalian membrane.
Shows phosphatidylserine (PS), sphingomyelin (SM), phosphatidylethanolamine (PE), phosphatidylcholine (PC), cholesterol (CL), and diacylglycerol (DG). Density maps for each lipid type are show on the right. Contact occupancy structures are visualized from the intracellular view, the voltage-sensing domain (VSD)-IV side view, the extracellular view and the VSD-II side view (left to right).
Figure 2—figure supplement 3
Contact occupancy distributions (left) and outlying residues (right) for all 12 lipid classes in the mammalian membrane.
Error bars show standard error, n = 3.
Figure 3 with 3 supplements
Binding of different phosphoinositide (PIP) species in enriched PIP simulations.
(A) Enriched PIP simulation system, with Naᵥ1.4 embedded in a POPC membrane (transparent gray) and 5% each of PIP1 (blue), PIP2 (purple), and PIP3 (pink) added to the cytoplasmic leaflet. (B) Representative snapshots from the five longest binding events from different replicates, showing the three different PIP species (PIP1 in blue, PIP2 in purple, and PIP3 in pink) binding to voltage-sensing domain (VSD)-IV and the DIII–IV linker. Naᵥ1.4 is shown in white with interacting residues on the DIV S4–S5 linker and the DIII–IV linker colored in green and orange, respectively. (C) A frequency distribution showing interaction times for each PIP species, defined as the length of a continuous period in which a PIP was within 0.7 nm of two VSD-IV binding site residues. (D) Frequency plots showing number of positive residues interacting with bound PIP in the DIII–IV linker (vertical) and VSD-IV (horizontal). (E) Minimum distance between binding residues on Naᵥ1.4 and bound PIPs lipid across simulation time for the five longest binding events, colored by distance and the type of PIP bound.
Figure 3—figure supplement 1
Minimum distance between binding residues on Naᵥ1.4 and bound phosphoinositides (PIPs) lipid across 31 long duration binding events (>10 µs), colored by distance and the type of PIP bound: PIP1 (blue-green), PIP2 (purple), and PIP3 (pink).
Figure 3—figure supplement 2
Minimum distance between binding residues on Naᵥ1.4 and bound phosphoinositides (PIPs) lipid across 31 long duration binding events (>10 µs), colored by distance and the type of PIP bound: PIP1 (blue-green), PIP2 (purple), and PIP3 (pink).
Figure 3—figure supplement 3
Phosphoinositide (PIP) occupancy at putative binding site residues in WT (wild-type) (brown) and when all eight binding site residues are mutated to leucine (L8, light blue) or glutamate (E8, dark blue).
Error bars show standard error, n = 3.
Figure 4 with 3 supplements
PI(4,5)P2 and PI(4)P binding to Naᵥ1.4 stabilizes the DIII–IV linker in atomistic and flexible coarse-grained simulations.
(A) Representative snapshots of PI(4,5)P2 bound from the voltage-sensing domain (VSD)-I side (purple stick) and PI(4)P bound from the VSD-III side (cyan stick), with six basic residues forming the binding site located on the DIII–IV linker (orange VDW representation) and VSD-IV S4–S5 linker (shown in green VDW representation) visualized from the intracellular face of the protein. (B) Proportion of frames where each of the binding site residues were identified to be within 4.5 Å with the different headgroup regions, P4, P5, and PO4, for PI(4,5)P2 (left) and PI(4)P (right). Error bars show standard error, n = 5. (C) Comparison of RMSF per carbon-alpha for simulations with and without bound PI(4,5)P2, showing residues on the S4–S5 linker and DIV linker with significant differences in mobility (student’s t-test, p-value <0.05). (D) Interaction network plots between the phosphoinositide (PIP) headgroup and basic binding residues on DIII–IV linker (orange) and DIV S4–S5 linker (green), generated by ProLIF – showing the dominant interactions across simulations of PI(4,5)P2 and PI(4)P. (E) Density plots showing differences in the distributions of distance between IFM/IQM motif and its binding pocket in the presence and absence of PIPs for the Naᵥ1.4 wild-type (left) and IFM->IQM mutant (right); with representative snapshots showing the two distinct conformations of the IQM motif in the mutant.
Figure 4—figure supplement 1
Root mean square deviation (RMSD) of the Naᵥ1.4 backbone, voltage-sensing domain (VSD)-IV backbone, DIII–IV linker backbone, and S4–S5 backbone (and PI(4,5)P2/PI(4)P1 headgroup) over 1.5 µs of atomistic simulations without any phosphoinositide (PIP) (top row), with a single PI(4,5)P2 lipid bound (middle row) and with PI(4)P1 bound (bottom row).
Five replicates were simulated – rep1 (blue), rep2 (red), rep3 (green), rep4 (orange), rep5 (purple).
Figure 4—figure supplement 2
Electrostatic interactions between the PI(4,5)P2 headgroup and binding residues.
ProLIF bar charts describing the occurrence of anionic and hydrogen bond interactions between the PI(4,5)P2 headgroup and binding residues; error bars show standard error, n = 5 (top). Per replicate data, showing the number of electrostatic interactions detected between the PI(4,5)P2 headgroup and each of the binding residues (left); atomic interactions summarized in a ProLIF LigPlot (right).
Figure 4—figure supplement 3
Electrostatic interactions between the PI(4)P1 headgroup and binding residues.
ProLIF bar charts describing the occurrence of anionic and hydrogen bond interactions between the PI(4)P1 headgroup and binding residues; error bars show standard error, n = 5 (top). Per replicate data, showing the number of electrostatic interactions detected between the PI(4)P1 headgroup and each of the binding residues (left); atomic interactions summarized in a ProLIF LigPlot (right).
Figure 5
Comparison of the identified phosphoinositide binding site to Naᵥ subtypes and other ion channels.
(A) Binding poses of PI(4,5)P2 (in purple) and PI(4)P (in cyan) aligned with two other tetrameric channels structures Kᵥ7.1 (6v01, in yellow) and Caᵥ2.2 (7mix, in red) that were resolved with PI(4,5)P2 at their respective voltage-sensing domains (VSDs). (B) Sequence alignment of the S4 helix and S4–S5 linker of the four domains of Naᵥ1.4, compared to VSD-II of Caᵥ2.2 and one of the four identical VSDs of Kᵥ7.1; residues colored by amino acid class; purple boxes indicate PI(4,5)P2 binding residues (identified with 5 Å of the headgroup). (C) Sequence alignment of the nine human Naᵥ channel subtypes shows high sequence similarity in the S4 helix, S4–S5 linker, and DIII–IV linker regions.
Figure 6 with 3 supplements
Phosphoinositide (PIP) binding to Naᵥ1.7 with different voltage-sensing domain (VSD) states in coarse-grained simulations.
(A) Atomistic representation of the three different Naᵥ1.7 structures simulated: (1) the inactivated state (blue, PDB ID: 6j8g) with the VSDs all in the activated, up state, (2) the Naᵥ1.7-NaᵥPas chimera (pink, PDB ID: 6nt4) with the Naᵥ1.7 VSD-IV in the deactivated, down state and a bound NaᵥPas C-terminal domain (CTD), (3) a Naᵥ1.7 resting state model (orange, model generation detailed in Materials and methods) with all four VSDs in the deactivated, down state. Panel insets show the different conformations of VSD-IV (left) and VSD-II (right) across different structures. The inactivation switch formed by the CTD and VSD-IV S4–S5 linker proposed by Clairfeuille et al. is shown (middle). (B) Combined occupancy of all PIP species (PIP1, PIP2, PIP3) at binding site residues in the three systems; error bars show standard error, n = 3. (C) Distribution of PIP binding durations at the identified site. (D) Intracellular view of CTD covering the resting state VSD-IV. Representative snapshot of PIP binding at DIV S4–S5 linker in the Naᵥ1.7-NaᵥPas system. The CTD (dark pink) prevents PIP access to DIII–IV linker lysines, K1492 and K1495. (E) Combined PIP occupancy at the bottom three gating charges on VSD I–III in the inactivated (blue) and resting state model (orange) simulations. For VSD-IV, PIP occupancy in the Naᵥ1.7-NaᵥPas system (pink) is also shown. (F) Representative simulation snapshot showing PIP (purple) binding at the gating charges (orange) in the resting state model simulations.
Figure 6—figure supplement 1
Naᵥ1.7 resting state model validation.
Each of the four voltage-sensing domains (VSDs) aligned separately, comparing the inactivated structure (6j8g, in blue) to the resting state model. All the S4 gating charges and the PHE/TYR residue on S2 to indicate the hydrophobic constriction site (HCS) for each VSD (in stick representation). The templates that each corresponding VSD was modeled from are shown in translucent representation. S3 helix hidden for better visualization of gating charge positions.
Figure 6—figure supplement 2
Minimum distance between binding residues on Naᵥ1.7 and bound phosphoinositides (PIPs) shown for long duration binding events (>20 µs) in each system, colored by distance and the type of PIP bound: PIP1 (blue-green), PIP2 (purple), and PIP3 (pink).
Five long-term binding events for the inactivated structure of Naᵥ1.7 (PDB ID: 6j8g) are shown on the left. Naᵥ1.7-NaᵥPas chimera (top right, PDB ID: 6nt4) and Naᵥ1.7 resting state model (bottom right) have fewer longer term interactions with PIP at this site. These interactions do not involve residues belonging to the DIII–IV linker (K1147/1492 and K1150/1495).
Figure 6—figure supplement 3
Phosphoinositide (PIP) binding at S4 gating charges.
(A) Binding duration distributions for PIP binding events at residues belonging to the bottom three gating charges in the S4 helix of voltage-sensing domain (VSD) I–III are shown for the inactivated structure of Naᵥ1.7 (PDB ID: 6j8g, all VSDs up) and resting Naᵥ1.7 model (all VSDs down, C-terminal domain [CTD] not modeled). For VSD-IV, the distribution for the Naᵥ1.7-NaᵥPas chimera (VSD-IV down, CTD bound) is also shown. (B) Minimum distance between terminal three gating charges on each VSD and bound PIP for long duration (>20 µs) binding events shown for the resting Naᵥ1.7 model (top) and inactivated Naᵥ1.7 structure (bottom). (C) Combined PIP occupancy at the bottom three gating charges on VSD I–IV in the inactivated (blue) and resting state model (orange) simulations of a Naᵥ1.4 resting state model; error bars show standard error, n = 3.
Figure 7
Proposed mechanism of phosphoinositide (PIP) effects on the sodium channel functional cycle.
Tables
Table 1
Disease causing point mutations at analogous phosphoinositide (PIP) binding residues in Naᵥ subtypes (described in the UniProt database).
| Naᵥ1.4 residue # | Analogous residue # | Subtype | Disease information; mechanism |
|---|---|---|---|
| K1330 | K1505N | Naᵥ1.5 | Long QT3 syndrome; unknown significance |
| R1463 | K1641N | Naᵥ1.2 | Benign familial infantile seizure; unknown significance |
| R1469 | R1657C | Naᵥ1.1 | Generalized epilepsy with febrile seizures plus; depolarizing shift in voltage dependence of activation, reduced current, accelerated recovery from slow inactivation |
| R1642C | Naᵥ1.3 | Developmental epileptic encephalopathy; accelerated recovery from inactivation | |
| R1644C R1644H | Naᵥ1.5 | Long QT3 syndrome Brugada syndrome |
Table 2
Generation of the resting state model using a combination of multiple templates for different voltage-sensing domains (VSDs) and pore domains (PDs) of Naᵥ1.7.
| Template pdb | Domains used | Res IDs (from template) | Res IDs (human Naᵥ1.7 numbering) | Template info |
|---|---|---|---|---|
| 7xve | VSDI, PDI, PDII | 1–404, 541–650 | 8–411, 864–972 | Naᵥ1.7 mutant (L866F, T870M, and A874F on S5II; V947F, M952F, and V953F on S6II; and V1438I, V1439F, and G1454C on S6III, E156K on S2I and G779R on S2II) VSDI deactivated (VSDII partially deactivated but not used) |
| 7k48 | VSDII, PDII, VSDIII, PDIII, PDIV | 405, 938, 1115–1244 | 728–973, 1175–1462, 1639–1768 | Naᵥ1.7/NaᵥAb chimera, where top half of each VSD is Naᵥ1.7 VSDII; all VSDs deactivated with engineered tarantula toxin m3-Huwentoxin-IV bound |
| 6nt4 | VSDIV, PDIV, PDI | 939–1244, 239–404 | 1463–1768, 246–411 | Naᵥ1.7/NaᵥPas chimera, full Naᵥ1.7 VSDI, and VSDIV deactivated with α-scorpion neurotoxin AaH2 bound |
| 7w9k | Loops between S5 and S6 regions for PDII, PDIII, and PDIV | 563–619, 811–918, 1141–1213 | 886–942, 1335–1442, 1665–1737 | Fast-inactivated Naᵥ1.7 (all VSDs activated) DIII–IV linker bound to pore |
Additional files
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
A binding site for phosphoinositides described by multiscale simulations explains their modulation of voltage-gated sodium channels
eLife 12:RP91218.
https://doi.org/10.7554/eLife.91218.3
