Hierarchical Approach to Predicting Permeation in Ion Channels

AbstractA hierarchical computational strategy combining molecular modeling, electrostatics calculations, molecular dynamics, and Brownian dynamics simulations is developed and implemented to compute electrophysiologically measurable properties of the KcsA potassium channel. Models for a series of channels with different pore sizes are developed from the known x-ray structure, using insights into the gating conformational changes as suggested by a variety of published experiments. Information on the pH dependence of the channel gating is incorporated into the calculation of potential profiles for K+ ions inside the channel, which are then combined with K+ ion mobilities inside the channel, as computed by molecular dynamics simulations, to provide inputs into Brownian dynamics simulations for computing ion fluxes. The open model structure has a conductance of ∼110 pS under symmetric 250mM K+ conditions, in reasonable agreement with experiments for the largest conducting substate. The dimensions of this channel are consistent with electrophysiologically determined size dependence of quaternary ammonium ion blocking from the intracellular end of this channel as well as with direct structural evidence that tetrabutylammonium ions can enter into the interior cavity of the channel. Realistic values of Ussing flux ratio exponents, distribution of ions within the channel, and shapes of the current-voltage and current-concentration curves are obtained. The Brownian dynamics calculations suggest passage of ions through the selectivity filter proceeds by a “knock-off” mechanism involving three ions, as has been previously inferred from functional and structural studies of barium ion blocking. These results suggest that the present calculations capture the essential nature of K+ ion permeation in the KcsA channel and provide a proof-of-concept for the integrated microscopic/mesoscopic multitiered approach for predicting ion channel function from structure, which can be applied to other channel structures