Mixed modes in opening of KcsA potassium channel from a targeted molecular dynamics simulation

Potassium channels conduct K⁺ flow selectively across the membrane through a central pore. During a process called gating, the potassium channels undergo a conformational change that opens or closes the ion-conducting pore. The potassium channel KcsA has been structurally determined in its closed state. However, the dynamic mechanism of the gating transition of the KcsA channel is still being investigated. Here, a targeted molecular dynamics simulation up to 150 ns is performed to investigate the detailed opening process of the KcsA channel with an open Kv1.2 structure serving as the target. The channel arrived at a self-determined quasi-stable state within 60 ns. The rigid-body and hinge-bending modes are observed mixed together in the remaining 90 ns long quasi-stable state. The mixed-mode movement seems come from the competition between the helix rigidity and the biased-applied gating force.     

Modeling diverse range of potassium channels with Brownian dynamics

Using the experimentally determined KcsA structure as a template, we propose a plausible explanation for the diversity of potassium channels seen in nature. A simplified model of KcsA is constructed from its atomic resolution structure by smoothing out the protein-water boundary and representing the atoms forming the channel protein as a homogeneous, low dielectric medium. The properties of the simplified and atomic-detail models, deduced from electrostatic calculations and Brownian dynamics simulations, are shown to be qualitatively similar. We then study how the current flowing across the simplified model channel changes as the shape of the intrapore region is modified. This is achieved by increasing the radius of the intracellular pore systematically from 1.5 to 5 A while leaving the dimensions of the selectivity filter and inner chamber unaltered. The strengths of the dipoles located near the entrances of the channel, the carbonyl groups lining the selectivity filter, and the helix macrodipoles are kept constant. The channel conductance increases steadily as the radius of the intracellular pore is increased. The rate-limiting step for both the outward and inward current is the time it takes for an ion to cross the residual energy barrier located in the intrapore region. The current-voltage relationship obtained with symmetrical solutions is linear when the applied potential is less than approximately 100 mV but deviates slightly from Ohm's law at higher applied potentials. The nonlinearity in the current-voltage curve becomes less pronounced as the radius of the intracellular pore is increased. When the strengths of the dipoles near the intracellular entrance are reduced, the channel shows a pronounced inward rectification. Finally, the conductance exhibits the saturation property observed experimentally. We discuss the implications of these findings on the transport of ions across the potassium channels and membrane channels in general.     

Intrinsic aqueduct orifices facilitate K+ channel gating

The ion-conducting pore of potassium channels, which can open and close to regulate ion passage, was at long thought to be a one-dimensional pore structure with a water-filled central cavity. Here, we find four orifices in the KcsA potassium channel, which are perpendicular to the pore and stretch out from the cavity. Equilibrium molecular dynamics simulations show that water molecules can flow between the cavity and orifices. Targeted molecular dynamics simulations show that during the opening process, water molecules can move into the cavity through the orifices to facilitate channel gating, whereas blocking the aqueduct orifices makes the channel difficult to open.     

Detection of the opening of the bundle crossing in KcsA with fluorescence lifetime spectroscopy reveals the existence of two gates for ion conduction

The closed KcsA channel structure revealed a crossing of the cytosolic ends of the transmembrane helices blocking the permeation pathway. It is generally agreed that during channel opening this helical bundle crossing has to widen in order to enable access to the inner cavity. Here, we address the question of whether the opening of the inner gate is sufficient for ion conduction, or if a second gate, located elsewhere, may interrupt the ion flow. We used fluorescence lifetime measurements on KcsA channels labeled with tetramethylrhodamine at residues in the C-terminal end of TM2 to report on the opening of the lower pore region. We found two populations of channels with different fluorescence lifetimes, whose relative distribution agrees with the open probability of the channel. The absolute fraction of channels found with an open bundle crossing is too high to explain the low open probability of the KcsA-WT channel. We found the same distribution as in the WT channel between open and closed bundle crossing for two KcsA mutants, A73E and E71A, which significantly increase open probability at low pH. These two results strongly suggest that a second gate in the ion permeation pathway exists. The location of the mutations A73E and E71A suggests that the second gate may be the selectivity filter, which resides in an inactivated state under steady-state conditions. Since the long closed times observed in KcsA-WT are not present in KcsA-A73E or -E71A, we propose that KcsA-WT remains predominantly in a state with an open bundle crossing but closed (inactivated) second gate, while the mutations A73E and E71A sharply decrease the tendency to enter in the inactivated state, and as a consequence, the second gate is predominantly open at steady state. The ability to monitor the opening of the bundle crossing optically enables the direct recording of the movement of the pore helices while the channel is functioning.     

KvAP-based model of the pore region of shaker potassium channel is consistent with cadmium- and ligand-binding experiments

Potassium channels play fundamental roles in excitable cells. X-ray structures of bacterial potassium channels show that the pore-lining inner helices obstruct the cytoplasmic entrance to the closed channel KcsA, but diverge in widely open channels MthK and KvAP, suggesting a gating-hinge role for a conserved Gly in the inner helix. A different location of the gating hinge and a narrower open pore were proposed for voltage-gated Shaker potassium channels that have the Pro-473-Val-Pro motif. Two major observations back the proposal: cadmium ions lock mutant Val-476-Cys in the open state by bridging Cys-476 and His-486 in adjacent helices, and cadmium blocks the locked-open double mutant Val-474-Cys/Val-476-Cys by binding to Cys-474 residues. Here we used molecular modeling to show that the open Shaker should be as wide as KvAP to accommodate an open-channel blocker, correolide. We further built KvAP-, MthK-, and KcsA-based models of the Shaker mutants and Monte-Carlo-minimized them with constraints Cys-476-Cd(2+)-His-486. The latter were consistent with the KvAP-based model, causing a small-bend N-terminal to the Pro-473-Val-Pro motif. The constraints significantly distorted the MthK-based structure, making it similar to KvAP. The KcsA structure resisted the constraints. Two Cd(2+) ions easily block the locked-open KvAP-based model at Cys-474 residues, whereas constraining a single cadmium ion to four Cys-474 caused large conformational changes and electrostatic imbalance. Although mutual disposition of the voltage-sensor and pore domains in the KvAP x-ray structure is currently disputed, our results suggest that the pore-region domain retains a nativelike conformation in the crystal.     

Basis for allosteric open-state stabilization of voltage-gated potassium channels by intracellular cations

The open state of voltage-gated potassium (Kv) channels is associated with an increased stability relative to the pre-open closed states and is reflected by a slowing of OFF gating currents after channel opening. The basis for this stabilization is usually assigned to intrinsic structural features of the open pore. We have studied the gating currents of Kv1.2 channels and found that the stabilization of the open state is instead conferred largely by the presence of cations occupying the inner cavity of the channel. Large impermeant intracellular cations such as N-methyl-d-glucamine (NMG⁺) and tetraethylammonium cause severe slowing of channel closure and gating currents, whereas the smaller cation, Cs⁺, displays a more moderate effect on voltage sensor return. A nonconducting mutant also displays significant open state stabilization in the presence of intracellular K⁺, suggesting that K⁺ ions in the intracellular cavity also slow pore closure. A mutation in the S6 segment used previously to enlarge the inner cavity (Kv1.2-I402C) relieves the slowing of OFF gating currents in the presence of the large NMG⁺ ion, suggesting that the interaction site for stabilizing ions resides within the inner cavity and creates an energetic barrier to pore closure. The physiological significance of ionic occupation of the inner cavity is underscored by the threefold slowing of ionic current deactivation in the wild-type channel compared with Kv1.2-I402C. The data suggest that internal ions, including physiological concentrations of K⁺, allosterically regulate the deactivation kinetics of the Kv1.2 channel by impairing pore closure and limiting the return of voltage sensors. This may represent a primary mechanism by which Kv channel deactivation kinetics is linked to ion permeation and reveals a novel role for channel inner cavity residues to indirectly regulate voltage sensor dynamics.     

Interaction of local anesthetics with the K+ channel pore domain

Local anesthetics and related drugs block ionic currents of Na⁺, K⁺ and Ca²⁺ conducted across the cell membrane by voltage-dependent ion channels. Many of these drugs bind in the permeation pathway, occlude the pore and stop ion movement. However channel-blocking drugs have also been associated with decreased membrane stability of certain tetrameric K⁺ channels, similar to the destabilization of channel function observed at low extracellular K⁺ concentration. Such drug-dependent stability may result from electrostatic repulsion of K⁺ from the selectivity filter by a cationic drug molecule bound in the central cavity of the channel. In this study we used the pore domain of the KcsA K⁺ channel protein to test this hypothesis experimentally with a biochemical assay of tetramer stability and theoretically by computational simulation of local anesthetic docking to the central cavity. We find that two common local anesthetics, lidocaine and tetracaine, promote thermal dissociation of the KcsA tetramer in a K⁺-dependent fashion. Docking simulations of these drugs with open, open-inactivated and closed crystal structures of KcsA yield many energetically favorable drug-channel complexes characterized by nonbonded attraction to pore-lining residues and electrostatic repulsion of K⁺. The results suggest that binding of cationic drugs to the inner cavity can reduce tetramer stability of K⁺ channels.     

Hierarchical approach to predicting permeation in ion channels

A 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 approximately 110 pS under symmetric 250 mM 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.     

Thermodynamic coupling between activation and inactivation gating in potassium channels revealed by free energy molecular dynamics simulations

The amount of ionic current flowing through K(+) channels is determined by the interplay between two separate time-dependent processes: activation and inactivation gating. Activation is concerned with the stimulus-dependent opening of the main intracellular gate, whereas inactivation is a spontaneous conformational transition of the selectivity filter toward a nonconductive state occurring on a variety of timescales. A recent analysis of multiple x-ray structures of open and partially open KcsA channels revealed the mechanism by which movements of the inner activation gate, formed by the inner helices from the four subunits of the pore domain, bias the conformational changes at the selectivity filter toward a nonconductive inactivated state. This analysis highlighted the important role of Phe103, a residue located along the inner helix, near the hinge position associated with the opening of the intracellular gate. In the present study, we use free energy perturbation molecular dynamics simulations (FEP/MD) to quantitatively elucidate the thermodynamic basis for the coupling between the intracellular gate and the selectivity filter. The results of the FEP/MD calculations are in good agreement with experiments, and further analysis of the repulsive, van der Waals dispersive, and electrostatic free energy contributions reveals that the energetic basis underlying the absence of inactivation in the F103A mutation in KcsA is the absence of the unfavorable steric interaction occurring with the large Ile100 side chain in a neighboring subunit when the intracellular gate is open and the selectivity filter is in a conductive conformation. Macroscopic current analysis shows that the I100A mutant indeed relieves inactivation in KcsA, but to a lesser extent than the F103A mutant.     

Structures of Gating Intermediates in a K+ channel

Regulation of ion conduction through the pore of a K⁺ channel takes place through the coordinated action of the activation gate at the bundle crossing of the inner helices and the inactivation gate located at the selectivity filter. The mechanism of allosteric coupling of these gates is of key interest. Here we report new insights into this allosteric coupling mechanism from studies on a W67F mutant of the KcsA channel. W67 is in the pore helix and is highly conserved in K⁺ channels. The KcsA W67F channel shows severely reduced inactivation and an enhanced rate of activation. We use continuous wave EPR spectroscopy to establish that the KcsA W67F channel shows an altered pH dependence of activation. Structural studies on the W67F channel provide the structures of two intermediate states: a pre- open state and a pre-inactivated state of the KcsA channel. These structures highlight key nodes in the allosteric pathway. The structure of the KcsA W67F channel with the activation gate open shows altered ion occupancy at the second ion binding site (S2) in the selectivity filter. This finding in combination with previous studies strongly support a requirement for ion occupancy at the S2 site for the channel to inactivate.     

Open-state models of a potassium channel

The structure of the bacterial potassium channel, KcsA, corresponds to the channel in a closed state. Two lines of evidence suggest that the channel must widen its intracellular mouth when in an open state: 1) internal block by a series of tetraalkylammonium ions and 2) spin labeling experiments. Thus it is known that the protein moves in this region, but it is unclear by how much and the mechanisms that are involved. To address this issue we have applied a novel approach to generate plausible open-state models of KcsA. The approach can be thought of as placing a balloon inside the channel and gradually inflating it. Only the protein sees the balloon, and so water is free to move in and out of the channel. The balloon is a van der Waals sphere whose parameters change by a small amount at each time step, an approach similar to methods used in free energy perturbation calculations. We show that positioning of this balloon at various positions along the pore axis generates similar open-state models, thus indicating that there may be a preferred pathway to an open state. We also show that the resulting structures from this process are conformationally unstable and need to undergo a relaxation process for up to 4 ns. We show that the channel can relax into a new state that has a larger pore radius at the region of the intracellular mouth. The resulting models may be useful in exploring models of the channel in the context of ion permeation and blocking agents.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Molecular docking of the scorpion toxin Tc1 to the structural model of the voltage-gated potassium channel Kv1.1 from human Homo sapiens

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K+ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K_v1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K_v1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) α-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K_v open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K_v1.5 channel is unable to close when bupivacaine is bound.     

Binding of kappa-conotoxin PVIIA to Shaker K+ channels reveals different K+ and Rb+ occupancies within the ion channel pore

The x-ray structure of the KcsA channel at different [K(+)] and [Rb(+)] provided insight into how K(+) channels might achieve high selectivity and high K(+) transit rates and showed marked differences between the occupancies of the two ions within the ion channel pore. In this study, the binding of kappa-conotoxin PVIIA (kappa-PVIIA) to Shaker K(+) channel in the presence of K(+) and Rb(+) was investigated. It is demonstrated that the complex results obtained were largely rationalized by differences in selectivity filter occupancy of this 6TM channels as predicted from the structural work on KcsA. kappa-PVIIA inhibition of the Shaker K(+) channel differs in the closed and open state. When K(+) is the only permeant ion, increasing extracellular [K(+)] decreases kappa-PVIIA affinity for closed channels by decreasing the "on" binding rate, but has no effect on the block of open channels, which is influenced only by the intracellular [K(+)]. In contrast, extracellular [Rb(+)] affects both closed- and open-channel binding. As extracellular [Rb(+)] increases, (a) binding to the closed channel is slightly destabilized and acquires faster kinetics, and (b) open channel block is also destabilized and the lowest block seems to occur when the pore is likely filled only by Rb(+). These results suggest that the nature of the permeant ions determines both the occupancy and the location of the pore site from which they interact with kappa-PVIIA binding. Thus, our results suggest that the permeant ion(s) within a channel pore can determine its functional and pharmacological properties.     

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Sequence-function analysis of the K+-selective family of ion channels using a comprehensive alignment and the KcsA channel structure

Sequence-function analysis of K(+)-selective channels was carried out in the context of the 3.2 A crystal structure of a K(+) channel (KcsA) from Streptomyces lividans (Doyle et al., 1998). The first step was the construction of an alignment of a comprehensive set of K(+)-selective channel sequences forming the putative permeation path. This pathway consists of two transmembrane segments plus an extracellular linker. Included in the alignment are channels from the eight major classes of K(+)-selective channels from a wide variety of species, displaying varied rectification, gating, and activation properties. Segments of the alignment were assigned to structural motifs based on the KcsA structure. The alignment's accuracy was verified by two observations on these motifs: 1), the most variability is shown in the turret region, which functionally is strongly implicated in susceptibility to toxin binding; and 2), the selectivity filter and pore helix are the most highly conserved regions. This alignment combined with the KcsA structure was used to assess whether clusters of contiguous residues linked by hydrophobic or electrostatic interactions in KcsA are conserved in the K(+)-selective channel family. Analysis of sequence conservation patterns in the alignment suggests that a cluster of conserved residues is critical for determining the degree of K(+) selectivity. The alignment also supports the near-universality of the "glycine hinge" mechanism at the center of the inner helix for opening K channels. This mechanism has been suggested by the recent crystallization of a K channel in the open state. Further, the alignment reveals a second highly conserved glycine near the extracellular end of the inner helix, which may be important in minimizing deformation of the extracellular vestibule as the channel opens. These and other sequence-function relationships found in this analysis suggest that much of the permeation path architecture in KcsA is present in most K(+)-selective channels. Because of this finding, the alignment provides a robust starting point for homology modeling of the permeation paths of other K(+)-selective channel classes and elucidation of sequence-function relationships therein. To assay these applications, a homology model of the Shaker A channel permeation path was constructed using the alignment and KcsA as the template, and its structure evaluated in light of established structural criteria.     

Molecular Docking of the Scorpion Toxin Tc1 to the Structural Model of the Voltage-gated Potassium Channel Kv1.1 from Human Homo sapiens

Abstract

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K⁺ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

Open-state conformation of the KcsA K+ channel: Monte Carlo normal mode following simulations

Potassium channels fluctuate between closed and open states. The detailed mechanism of the conformational changes opening the intracellular pore in the K+ channel from Streptomyces lividans (KcsA) is unknown. Applying Monte Carlo normal mode following, we find that gating involves rotation and unwinding of the TM2 bundle, lateral movement of the TM2 helices away from the channel axis, and disappearance of the TM2 bundle. The open-state conformation of KcsA exhibits a very wide inner vestibule, with a radius approximately 5-7 A and inner helices bent at the A98-G99 hinge. Computed conformational changes demonstrate that spin labeling and X-ray experiments illuminate different stages in gating: transition begins with clockwise rotation of the TM2 helices ending at a final state with the TM2 bend hinged near residues A98-G99. The concordance between the computational and experimental results provides atomic-level insights into the structural rearrangements of the channel's inner pore.     

Influence of interionic interactions on functional state and blocker binding of voltage-gated potassium channels

A mechanism of ion conduction of a voltage-gated potassium channel KcsA was investigated in full-atomic approximation at a trajectory length of 100 ns using the Lomonosov supercomputer. Methods of molecular dynamics were employed. A structure of the KcsA channel in the open state obtained by X-ray structure analysis (PDB ID 3fb7) was used. Free energy profiles of the KcsA pore occupied with either one or three potassium ions were calculated. It was shown that, under physiological conditions, ions pass through the channel pore cooperatively and the mechanism most probably includes three ions permeating in concert. Interactions of the mammalian voltage-gated channel Kv1.2 with neurotoxin were investigated. It was demonstrated that the effect of interionic interactions on binding of a blocker is rather insufficient.