Modeling of open, closed, and open-inactivated states of the hERG1 channel: structural mechanisms of the state-dependent drug binding.

The human ether-a-go-go related gene 1 (hERG1) K ion channel is a key element for the rapid component of the delayed rectified potassium current in cardiac myocytes. Since there are no crystal structures for hERG channels, creation and validation of its reliable atomistic models have been key targets in molecular cardiology for the past decade. In this study, we developed and vigorously validated models for open, closed, and open-inactivated states of hERG1 using a multistep protocol. The conserved elements were derived using multiple-template homology modeling utilizing available structures for Kv1.2, Kv1.2/2.1 chimera, and KcsA channels. Then missing elements were modeled with the ROSETTA De Novo protein-designing suite and further refined with all-atom molecular dynamics simulations. The final ensemble of models was evaluated for consistency to the reported experimental data from biochemical, biophysical, and electrophysiological studies. The closed state models were cross-validated against available experimental data on toxin footprinting with protein-protein docking using hERG state-selective toxin BeKm-1. Poisson-Boltzmann calculations were performed to determine gating charge and compare it to electrophysiological measurements. The validated structures offered us a unique chance to assess molecular mechanisms of state-dependent drug binding in three different states of the channel.