K+ channels: gating mechanisms and lipid interactions

Thesis

Computational methods, including homology modelling, in-silico dockings, and molecular dynamics simulations have been used to study the functional dynamics and interactions of K⁺ channels. Molecular models were built of the inwardly rectifying K⁺ channel Kir2.2, the bacterial homolog K⁺ channel KirBac3.1, and the twin pore (K2P) K⁺ channels TREK-1 and TRESK. To investigate the electrostatic energy profile of K⁺ permeating through these homology models, continuum electrostatic calculations were performed. The primary mechanism of KirBac3.1 gating is believed to involve an opening at the helix bundle crossing (HBC). However, simulations of Kir channels have not yet revealed opening at the HBC. Here, in simulations of the new KirBac3.1-S129R X-ray crystal structure, in which the HBC was trapped open by the S129R mutation in the inner pore-lining helix (TM2), the HBC was found to exhibit considerable mobility. In a simulation of the new KirBac3.1-S129R-S205L double mutant structure, if the S129R and the S205L mutations were converted back to the wild-type serine, the HBC would close faster than in the simulations of the KirBac3.1-S129R single mutant structure. The double mutant structure KirBac3.1-S129R-S205L therefore likely represents a higher-energy state than the single mutant KirBac3.1-S129R structure, and these simulations indicate a staged pathway of gating in KirBac channels. Molecular modelling and MD simulations of the Kir2.2 channel structure demonstrated that the HBC would tend to open if the C-linker between the transmembrane and cytoplasmic domain was modelled helical. The electrostatic energy barrier for K⁺ permeation at the helix bundle crossing was found to be sensitive to subtle structural changes in the C-linker. Charge neutralization or charge reversal of the PIP2-binding residue R186 on the C-linker decreased the electrostatic barrier for K⁺ permeation through the HBC, suggesting an electrostatic contribution to the PIP2-dependent gating mechanism. Multi-scale simulations determined the PIP2 binding site in Kir2.2, in good agreement with crystallographic predictions. A TREK-1 homology model was built, based on the TRAAK structure. Two PIP2 binding sites were found in this TREK-1 model, at the C-terminal end, in line with existing functional data, and between transmembrane helices TM2 and TM3. The TM2-TM3 site is in reasonably good agreement with electron density attributed to an acyl tail in a recently deposited TREK-2 structure.

Authors: Schmidt, M

• Publication date: 2013
• DOI: 10.5287/ora-yqazwnazm
• Type of award: DPhil
• Level of award: Doctoral
• Awarding institution: University of Oxford
• Language: English
• Keywords: molecular dynamics simulations; phosphatidyl inositol phosphates; homology modelling; ion channel gating; small molecule docking; potassium ion channel; complex bilayers; membrane proteins; Ion channels
• Subjects: Computational biochemistry; Molecular biophysics (biochemistry); Biochemistry; Life Sciences