Modeling the hERG potassium channel in a phospholipid bilayer: Molecular dynamics and drug docking studies

The hERG potassium channel has recently been a matter of extensive studies both at experimental and computational levels, because of its possible involvement in the potentially lethal drug-induced long QT syndrome. In this context, in the absence of an experimentally determined 3D structure, to acquire a detailed description of the channel pore and an understanding of atomic determinants for drug binding is of enormous interest at both academic and industrial levels. To contribute to this aim, we first built the open and closed states of the channel by homology modeling techniques, and then submitted both channel models to ns-time-scale MD simulations in explicit membrane. An in-depth analysis of the dynamical behavior of the channel pore with particular attention to the cavity volume was carried out. Finally, using hERG conformations coming from MD simulations, docking experiments were performed to identify a possible binding mode for the most potent hERG channel blocker so far known, the antihistaminic drug astemizole. We show that the combined use of MD and docking is suitable to identify a possible binding mode of drugs in a fairly good agreement with experiments. Moreover, the exploitation of MD snapshots in the docking experiments allowed us to capture some induced-fit effects related to the side chain conformations of Tyr652 and Phe656, which are residues playing a pivotal role in the hERG drug binding.