Interplay between VSD, pore, and membrane lipids in electromechanical coupling in HCN channels

Figures

• 

  Figure 1.. Position and contact network of W281. The three cryoEM structures of the resting (A; PDB ID 5U6P), intermediate (B; PDB ID 6UQG), and activated (C; PDB ID 6UQF) states of hHCN1, shown in blue, yellow, and red, respectively, were investigated by equilibrium molecular dynamics (MD) simulations. The position of W281 is highlighted in ball/stick, and transparent surface and surrounding residues are shown in ball and stick in the resting (D), intermediate (E), and activated (F) conformations. (G) Contact occurrence between W281 and nearby residues and lipid molecules along 1-μs-long MD simulations.

• 

  Figure 1—figure supplement 1.. S4-S5 and S4-S1 hydrophobic interfaces in the resting, intermediate, and activated states. Probability distributions of (A) the hydrophobic contact areas at the S4-S5 (A, top panel – measured between the innermost parts of S4 [residues 277–289] and S5 [residues 290–305]) and S4-S1 (A, bottom panel – measured between the bottom parts of S4 [res 277–279] and S1 [res 139–150]) interfaces and (B) the W281-N300 distances along the equilibrium simulations of the resting (blue), intermediate (yellow), and activated (red) states.

• 

  Figure 2.. HCN1 activation induces an increase in the lower S4/S5-lipids hydrophobic interfaces. (A, D) First and last snapshots of voltage-sensor domains (VSDs) extracted from activating simulations of intact (panel A, subunit D) or HCND-less system (panel D, subunit A). The gating charges and residues W281, Y277 (S4), N300 (S5), and F143 (S1) are highlighted using ball and stick representation. (B, E) Downward displacement of the gating charges with respect to the charge transfer center F186 during the activation of S4 in subunit D (panel B, intact HCN1 model) and subunit A (panel E, HCND-less HCN1 model) during activating simulations. (C, F) Correlation between the S4 bending angle and the rotation of W281 in subunit D (panel C, intact HCN1 model) and subunit A (panel F, HCND-less HCN1 model) during activating simulations. The time evolution is represented by a color gradient, from white/light pink to dark pink. (G) Evolution of the S4/S5 (brown) and S4/S1 (gray) hydrophobic interfaces (contact areas) along activating simulations of HCN1 in the presence (subunit D, top box) or absence (subunit A, bottom box) of the HCN domain. (H) Evolution of the number of lipid carbon atoms, N, in…

• 

  Figure 2—figure supplement 1.. S4 displacement in simulations of the intact (A, 1 µs) or in the HCND-less system (B, 750 ns). In both panels, the first two rows represent the gating charges (R267, R270, R273, and R276 colored from light gray to black) movement with respect to the charge transfer center (F186, z = 0) for the four protein chains (chains A, B, C, and D). The last rows show the evolution of the bending angle of S4 (defined as the angle between segments encompassing residues 252–272 and 272–289) and the rotation angle of W281 (right panel, defined as the angle between atoms Cɑ-V296 [on S5], Cɑ-W281 and CE2-W281 [on S4]). The four monomers, chains A, B, C, and D, are colored from light gray to black.

• 

  Figure 2—figure supplement 2.. Evolution of the S4/S5 and S4/S1 interfaces in activating molecular dynamics (MD) simulations of intact (A–C) and in HCND-less HCN1 models (D–F). Each column represents a different subunit. First row: last snapshot from Hyperpolarizing EF-simulation showing the position of W281 and Y277 with respect to F143 and N300. Second row: hydrophobic (gray) and electrostatic (green) contact areas measured between the innermost parts of S4 (residues 277–289) and S5 (residues 290–305) as a function of simulation time. Third row: hydrophobic (gray) and electrostatic (green) contact areas measured between the bottom parts of S4 (res 277–279) and S1 (res 139–150) as a function of simulation time.

• 

  Figure 2—figure supplement 3.. Hydrophobic and electrostatic interactions in the lower S4/S5 region are state-dependent. Average contact areas between each residue/subunit from the lower S4 (left) and S5 (right) and surrounding hydrophobic (light gray) and hydrophilic (spring green) protein residues, lipid tails (dark gray) and headgroups (light green), and water molecules (turquoise) in the resting (panel A) and activated (panel B) states equilibrium molecular dynamics (MD) simulations (1 µs). (C) Contact area between lower S4/S5 protein residues and lipid tails averaged over the four subunits and over the simulation time, in the resting (blue) and activated (red) states. Error bars are given as standard deviation calculated over the four subunits and over the simulation time.

• 

  Figure 3.. Structural comparison of the S4/S5 linker and C-linker regions. (A) Structural comparison of the S4-S6 region of the resting (5U6P; blue), intermediate (6UQG; yellow), and activated (6UQF; red) states of hHCN1. Alignment was performed on the entire structure. (B) Structural comparison of the S4/S5 linker and C-linker regions of the resting, intermediate, and activated states of hHCN1. Residue K412 of the C-linker is shown as grey spheres. (C) D290-K412 distance distribution in simulations of the resting, intermediate, and activated states. Dashed lines represent distances in the starting cryoEM structures of the same color. (D) Peak tail current versus voltage plot for whole-cell recordings from HEK cells expressing HCN1 D290C/K412C with (red) or without (black) 100 uM CdCl2 in the patch pipette. Error bars represent SEM for n = 5 independent cells for each condition.

• 

  Figure 3—figure supplement 1.. Functional recordings of the D290C-K412C double mutant. (A) Comparison of the relative peak tail current versus voltage relationships for HCN1-EM with (pink) and without (gray) cadmium chloride and the D290C-K412C double mutant with (red) and without (black) cadmium chloride in the patch pipette. Error bars represent SEM for n = 5 independent cells for each condition. (B) Example of cesium blockage and residual leak current from cells expressing the D290C-K412C mutant. Currents are elicited from the same cell from test pulses to –140 mV with (red) and without (black) 1 mM CsCl present in the bath solution. (C) Same as (B) showing just the test pulse to highlight the lack of time-dependence to the ‘instantaneous’ current component (Iinstant).

• 

  Figure 4.. Pore opening of HCN1. (A) Residues and distances used as collective variables in adiabatic bias molecular dynamics (ABMD) simulations are highlighted by ball and stick and black double arrows, respectively. (B) Inter-subunit distances between opposite V390 residues are represented as a function of time in the resting (blue), activated (red), and ABMD (crimson) simulations starting from the closed activated (A1/A2 in light purple and light orchid) and resting (R1/R2 in dark purple and dark orchid) states. (C) Contact occurrences of R297, D401, and H392 in HCN1 are shown for the resting (blue), activated (red), and ABMD (crimson) simulations starting from the closed activated (A1/A2 in light purple and light orchid). (D) Lower S4/S5, S5/S6, S6/S6(bis) S4/CL(bis) interfaces/interactions. Important/interesting residues from these interfaces are highlighted in ball and stick.

• 

  Figure 4—figure supplement 1.. Open pore models from adiabatic bias molecular dynamics (ABMD) simulations compared to HCN4. Left and middle: Overlay of pore domains from the start (blue) and end (purple) for two ABMD replicates starting from the activated-closed state. Right: overlay of the closed (white) and open (gray) pores from HCN4. Spheres highlight the alpha carbon of G382 for HCN1 or G503 for HCN4. Sticks highlight V390 in HCN1 and the corresponding position I511 in HCN4. Only two opposing subunits are shown for clarity.

• 

  Figure 4—figure supplement 2.. Potassium permeation along the HCN1 pore. Positions of potassium ions along the Z-axis of the simulation box, each shown in a different symbol and color, as a function of the simulation time. The membrane is located between Z = ~90 Å and Z = ~120 Å. Data from six simulations are shown: resting and activated/closed HCN1 structures equilibrium molecular dynamics (MD) simulations, and adiabatic bias molecular dynamics (ABMD) simulation of the transition toward open channel starting from resting/closed (R1 and R2) and activated/closed (A1 and A2) HCN1 structures.

• 

  Figure 4—figure supplement 3.. Key state-dependent interactions. (A) Contact occurrences of H279 and E282 from S4 and R404 and R405 from the C-linker in HCN1 are shown for the resting (blue), activated (red), and adiabatic bias molecular dynamics (ABMD) simulations starting from the closed activated (A1/A2 in light purple and light orchid). (B, C) Probability distributions of H392-Q398 (B) and R404-E282 distances (C) along the equilibrium simulations of the resting (blue), activated (red), and the ABMD simulations (A1 in purple and A2 in orchid).