Molecular dynamics simulation of TMEM16A channel: Linking structure with gating

Abstract

TMEM16A, the calcium-activated chloride channel, is broadly expressed and plays pivotal roles in diverse physiological processes. To understand the structural and functional relationships of TMEM16A, it is necessary to fully clarify the structural basis of the gating of the TMEM16A channel. Herein, we performed the protein electrostatic analysis and molecular dynamics simulation on the TMEM16A in the presence and absence of Ca2+. Data showed that the separation of TM4 and TM6 causes pore expansion, and Q646 may be a key residue for the formation of π-helix in the middle segment of TM6. Moreover, E705 was found to form a group of H-bond interactions with D554/K588/K645 below the hydrophobic gate to stabilize the closed conformation of the pore in the Ca2+-free state. Interestingly, in the Ca2+ bound state, the E705 side chain swings 100o to serve as Ca2+-binding coordination and released K645. K645 is closer to the hydrophobic gate in the calcium-bound state, which facilitates the provision of electrostatic forces for chloride ions as the ions pass through the hydrophobic gate. Our findings provide the structural-based insights to understanding the mechanisms of gating of TMEM16A.