Force distribution during the push-to-open gating process of NOMPC.

(a) The simulation systems of NOMPC. The NOMPC molecule was divided into two subsystems, denoted by the orange and purple rectangular boxes, for molecular dynamics (MD) simulations. The interface between the linker helices (LH) and the TRP domains is shown in a circular inset (bottom view). (b) The membrane-parallel (purple line) and membrane-normal (yellow line) components of the net forces on the TRP domain exerted by the LH domain. The error bars denote the standard deviations from two compressive trajectories (CI-1/2) and two free trajectories (FI-1/2). (c) The net forces on the TRP domain exerted by the LH domain. The grey quiver on each residue showed the net force exerted by the LH domain. The black showed the resultant force on the residues of the TRP domain. Two major points of force application, E1571 and R1581, were marked. ⊗ represented the inward force component. The data were averaged for the four subunits by rotationally symmetrized around the protein center.

A membrane-parallel torsion force can open NOMPC.

(a-b) The side view (a) and bottom view (b) of the net-FDA of the E1571 and R1581 on the TRP domain exerted from the LH domain, which was decomposed into the membrane-parallel and membrane-normal components. ⊗ represented membrane-normal forces pointing into the membrane. (c-d) The transmembrane pore size evolution in the simulations with membrane-parallel forces (c) or membrane-normal forces (d), three replicates for each condition (trajectories: MPI-1/2/3, MNI-1/2/3, see Methods).

PDBID of resolved cryo-EM structures uesd in Fig. 6.

Torsional mechanical properties of the AR bundle.

(a) The sketch of the steered molecular dynamics simulation. The dotted arrows revealed the direction of twisting. (b) The relationships between the torque on AR1-8 (left Y-axis), the length of the AR domain (right Y-axis), and the rotational angle of AR1-8. This was calculated from the rotating trajectories RII (0-4).

Force transmission along the AR domain.

(a) The quiver plot of net-FDA (XY planar) of each (i + 1)th AR unit exerted from the ith AR unit showed the force transmission from the intra-cellular terminal to the transmembrane side of the AR domain viewing from AR1 toward AR29. ⊙ represented the outward Z-component force. (b) The distribution of the pairwise residue-residue net-FDA within AR9 to AR29 (weak interaction under 5 pN neglected). The inset showed local information over 70 pN. (c) The representative interaction heatmaps (AR13-AR15) showed the magnitude of non-bonded pairwise residue-residue net-FDA (left) and hydrogen bonds occupancy (right). (d) The Pearson correlation coefficient (PCC) between net-FDA and hbonds in AR units. Inset showed the mean value and standard deviation of Intra-AR (diagonal) and Inter-AR (sub-diagonal) PCCs.

A schematic diagram of the twist-to-open model.

The side view (left) and bottom view (right) of the twist-to-open model. Both the compression and clockwise twisting of the AR domain can generate a membrane-parallel twisting force on the TRP domain, which is the key component for gating.

Conformational changes of the TRP domain in various TRP channels with resolved closed- and open-state cryo-EM structures.

(a) The membrane-parallel rotational angle of the TRP helix (a positive value indicates clockwise rotation when viewed from the intracellular side). (b) The tilt angle of the TRP domain relative to the membrane surface (a positive value indicates tilt-up). S1-S4 helices were aligned for analysis. (c) and (d) represent the rotational and tilting conformational changes of the TRP domain upon gating, respectively. In the majority of the studied cases, the TRP domain rotates clockwise (c). The exception of counterclockwise rotation is observed for HsTRPA1. (please refer to Table 1 for the protein structures used in this analysis)

MD simulation systems.

Simulation trajectories of system I.

Simulation trajectories of system II.

Kinematic analysis of the push-to-open trajectory CI-1.

(a) The simulation systems and domain structures. (b) The pore dilation process in the push-to-open trajectory. (c) The motion tendency of the LH domain. The purple arrows represent each α-C, and the orange arrows represent the mass center of the LH domain. The data were averaged for the four subunits by rotationally symmetrized around the protein center. (d) The rotation of the TRP domain, projected in XY plane (light purple line, positive values a clockwise rotation viewing from the intracellular side) and Z direction (dark purple line, the positive value meant a tilt-up rotation from intracellular side to extracellular side), respectively.

The pairwise FDA results of the push-to-open trajectory CI-1 between the AR29 domain, the LH domain, and the TRP domain.

The size and color of each marker denoted the magnitude of the corresponding pairwise FDA results labeled in the two sides. The circle marker and square marker showed the FDA results between AR29-LH and TRP-LH, respectively.

The FDA results of the push-to-open trajectories between the LH domain and the TRP domain on the XY plane (TRP).

A superimposed quiver plot of the FDA results for each residue of the TRP domain, exerted by the LH domain, was calculated from the trajectories FI-1/2 and CI-1/2.

The FDA results of the push-to-open trajectory CI-1 between the TRP and transmembrane (TM) domain.

(a) The scalar FDA of the TM domain exerted by the TRP domain. The size of each marker indicates the magnitude of the corresponding pairwise interaction in FDA. The TM domain is subdivided into four regions with strong interactions, which are illustrated by circles of distinct colors: mint green for S1, violet for the S2-S3 linker, pink for the S4-S5 linker, and yellow for S6. (b) The quiver plot illustrates the membrane-normal component of the net-FDA of the TM domain exerted by the TRP domain. (c) Intracellular view and (d) side view of the structure. The regions of significant interactions are highlighted using the same color code as in (a).

Water molecules around pore region.

(a) A snapshot of trajectory MPI-1 at 800 ns. The pore was dilated so water molecules (red) could spontaneously pass through. (b) The number of water molecules in the gate region around I1554, which was the narrowest part. We counted the water molecules whose oxygen atom was within 4 Å in the Z-direction from the α-C of I1554.

The transmembrane pore size evolution in the simulations with (a) half or (b) one-third magnitude of original membrane-parallel force, three replicates for each condition (trajectories: MPI-h1/2/3, MPI-t1/2/3).

The torque (a) and rotational angle (b) versus time of trajectory RII, respectively.

The black line illustrated the trajectory RII-0 (the second step of the rotation of system II), while the other colored lines starting from a circle depicted the third step of the rotation of system II at every 2.5-degree rotation for a duration of 100 ns in the trajectory RII-1/2/3/4. The solid line corresponded to the data of the AR1-8, whereas the dotted line corresponded to the LH domain (Fig. 3).

Torsional mechanical properties of the AR bundle.

(a) The sketch of the steered molecular dynamics simulation. The dotted arrows revealed the direction of twisting. (b) The relationships between the torque on AR1-8 (left Y-axis), the length of the AR domain (right Y-axis), and the rotational angle of AR1-8. This was calculated from the rotating trajectories RII-fast (0-4). The torque (c) and rotational angle (d) versus time of trajectory RII-fast, respectively. The black line illustrated the trajectory RII-fast-0 (the second step of the rotation of system II), while the other colored lines starting from a circle depicted the third step of the rotation of system II at every 2.5-degree rotation for a duration of 100 ns in the trajectory RII-fast-1/2/3/4. The solid line corresponded to the data of the AR1-8, whereas the dotted line corresponded to the LH domain. The fitted torsion coefficient from RII-fast is (2.31 ± 0.44) × 103 kJ mol−1 rad−1 and the compression–twist coupling coefficient, , is (1.67 ± 0.14) nm rad−1.

Location of AR13-15 in the AR spring.

Rotation of the TRP domain of TRPV1 under torsional or pushing forces.

(a) Structure of squirrel TRPV1. (b) Overlay of TRPV1 conformations at 0 ns (grey) and 500 ns (red) under torsional force, as viewed from the intracellular side. The dark red denotes the AR1-2 units, while the light red indicates the TRP domain. (c) Rotation of AR1-2 units (dark red, rotating around the center of TRPV1) and the TRP domain (light red, self-rotating) under torsional force. (d) Similar to (b), but under pushing force. (e) Rotation of the TRP domain under pushing force (light red, self-rotating).

Secondary structure analysis of the TRP domain in the trajectory MPI-1.

The vertical coordinates denote all residue IDs of the TRP domains across the four chains. The analysis was conducted with DSSP.

Secondary structure analysis of the force-acting AR units in the trajectory RII-0.

The vertical coordinates denote all residue IDs of the AR9-AR11 units of one chain. The analysis was conducted with DSSP.