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.     

A voltage-sensor water pore

Voltage-sensor (VS) domains cause the pore of voltage-gated ion channels to open and close in response to changes in transmembrane potential. Recent experimental studies suggest that VS domains are independent structural units. This independence is revealed dramatically by a voltage-dependent proton-selective channel (Hv), which has a sequence homologous to the VS domains of voltage-gated potassium channels (Kv). Here we show by means of molecular dynamics simulations that the isolated open-state VS domain of the KvAP channel in a lipid membrane has a configuration consistent with a water channel, which we propose as a common feature underlying the conductance of protons, and perhaps other cations, through VS domains.     

Salt bridges in the miniature viral channel Kcv are important for function

The viral potassium channel Kcv comprises only 94 amino acids, which represent the pore module of more complex K⁺ channels. As for Kir-type channels, Kcv also has a short N-terminal helix exposed to the cytoplasm, upstream of the first transmembrane domain. Here we show that this helix is relevant for Kcv function. The presence of charged amino acids, which form dynamic inter- and intra-subunit salt bridges is crucial. Electrophysiological measurements, yeast rescue experiments and molecular dynamics simulations show that mutants in which the critical salt bridge formation is impaired have no or reduced channel activity. We conclude that these salt bridges destabilise the complexation of K⁺ ions by negative charges on the inner transmembrane domain at the entrance into the cavity. This feature facilitates a continuous and coordinated transfer of ions between the cavity and the cytoplasm for channels without the canonical bundle crossing.     

Water and potassium dynamics inside the KcsA K(+) channel

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.     

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.     

Molecular dynamics simulations and KcsA channel gating

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.     

Molecular simulation of the Kv7.4[ΔS269] mutant channel reveals that ion conduction in the cavity is perturbed due to hydrophobic gating

Mutations in the voltage-gated potassium channel Kv7.4 (encoded as KCNQ4) lead to the early onset of non-syndromic hearing loss, which is significant during language acquisition. The deletion of the S269 pore residue (genetic Δ mutation) in Kv7.4 has been reported to be associated with hearing loss. So far, there is no mechanistic understanding of how this mutation modulates channel function. To understand the role of S269 in ion conduction, we performed molecular dynamics simulations for both wild type and ΔS269 mutant channels. Simulations indicate that the ΔS269 mutation suppresses the fluctuations in the neighboring Y269 residue and thereby consolidates the ring formed by I307 and F310 residues in the adjacent S6 helixes in the cavity region. We show that the long side chains of I307 near the entrance to the cavity form a hydrophobic gate. Comparison of the free energy profiles of a cavity ion in Kv7.4 and Kv7.4[ΔS269] channels reveals a sizable energy barrier in the latter case, which suppresses ion conduction. Thus the simulation studies reveal that the hydrophobic gate resulting from the ΔS269 mutation appears to be responsible for sensorineural hearing loss.     

Molecular dynamics studies of the archaeal translocon

The translocon is a protein-conducting channel conserved over all domains of life that serves to translocate proteins across or into membranes. Although this channel has been well studied for many years, the recent discovery of a high-resolution crystal structure opens up new avenues of exploration. Taking advantage of this, we performed molecular dynamics simulations of the translocon in a fully solvated lipid bilayer, examining the translocation abilities of monomeric SecYEbeta by forcing two helices comprised of different amino acid sequences to cross the channel. The simulations revealed that the so-called plug of SecYEbeta swings open during translocation, closing thereafter. Likewise, it was established that the so-called pore ring region of SecYEbeta forms an elastic, yet tight, seal around the translocating oligopeptides. The closed state of the channel was found to block permeation of all ions and water molecules; in the open state, ions were blocked. Our results suggest that the SecYEbeta monomer is capable of forming an active channel.     

Brownian dynamics theory for predicting internal and external blockages of tetraethylammonium in the KcsA potassium channel

The theory of Brownian dynamics is used to model permeation and the blocking of KcsA potassium channels by tetraethylammonium (TEA). A novel Brownian dynamics simulation algorithm is implemented that comprises two free energy profiles; one profile is seen by the potassium ions and the other by the TEA molecules whose shape is approximated by a sphere. Our simulations reveal that internally applied TEA blocks the passage of K⁺ ions by physically occluding the pore. A TEA molecule in the external reservoir encounters an attractive energy-well created by four tyrosine residues at position 82, in addition to all other attractive and repulsive forces impinging on it. Using Brownian dynamics, we investigate how deep the energy-well needs to be to reproduce the experimentally determined inhibitory constant k_i for the TEA blockade of KcsA or the mutant Shaker T449Y. The one-dimensional free energy profile obtained from molecular dynamics is first converted into a one-dimensional potential energy profile, and is then transformed into a three-dimensional free energy profile in Brownian dynamics by adding the short-range potential from the channel walls. When converted, the free energy profile calculated from molecular dynamics gives a well-depth of ∼10 kT. We systematically alter the depths of the profiles, and then use Brownian dynamics simulations to numerically determine the current versus TEA-concentration curves. We show that the sequence of binding and unbinding events of the TEA molecule to the binding pocket can be modeled by a first-order Markov process. The Brownian dynamics simulations also reveal that the probability of a TEA molecule binding to the binding pocket in KcsA potassium channels increases exponentially with TEA concentration and depends also on the applied potential and the K⁺ concentration in the simulation assembly.     

Free energy of a potassium ion in a model of the channel formed by an amphipathic leucine-serine peptide

We use molecular dynamics simulations to investigate the position-dependent free energy of a potassium ion in a model of an ion channel formed by the synthetic amphipathic leucine-serine peptide, LS3. The channel model is a parallel bundle of six LS3 helices around which are packed 146 methane-like spheres in order to mimic a membrane. At either end of and within the channel are 1051 water molecules, plus four ions (two potassium and two chloride). The free energy of a potassium ion in the channel was estimated using the weighted histogram analysis (WHAM) method. This is the first time to our knowledge that such a calculation has been carried out as a function of the position of an ion in three dimensions within a channel. The results indicate that for this channel, which is lined by hydrophilic serine sidechains, there is a relatively weak dependence of the free energy on the axial/off-axial position of the ion. There are some off-axis local minima, especially in the C-terminal half of the channel. Using the free energy results, a single channel current-voltage curve was estimated using a one-dimensional Nernst-Planck equation. Although reasonable agreement with experiment is achieved for K(+) ions flowing from the N-terminal to the C-terminal mouth, in the opposite direction the current is underestimated. This underestimation may be a consequence of under-sampling of the conformational dynamics of the channel. We suggest that our simulations may have captured, for example, a sub-conductance level (i.e. an incompletely open state) of the LS3 channel.     

Dynamics, hydration, and motional averaging of a loop-gated artificial protein cavity: the W191G mutant of cytochrome c peroxidase in water as revealed by molecular dynamics simulations

Five molecular dynamics simulations of the W191G cavity mutant of cytochrome c peroxidase in explicit water reveal distinct dynamic and hydration behavior depending on the closed or open state of the flexible loop gating the cavity, the binding of (K+ or small molecule) cations, and the system temperature. The conformational spaces sampled by the loop region and by the cavity significantly reduce upon binding. The largest ordering factor on water dynamics is the presence of the K+ ion occupying the gated cavity. Considerable water exchange occurs for the open-gate cavity when no ligand or cation is bound. In all cases, good correspondence is found between the calculated (ensemble-averaged) location of water molecules and the water sites determined by X-ray crystallography experiments. However, our simulations suggest that these sites do not necessarily correspond to the presence of bound water molecules. In fact, individual water molecules may repeatedly exchange within the cavity volume yet occupy on average these water sites. Four major conclusions emerge. First, it seems misleading to interpret the conformation of protein loop regions in terms of single dominant structures. Second, our simulations support the general picture of Pro 190 cis-trans isomerization as a determinant of the loop-opening mechanism. Third, receptor flexibility is fundamental for ligand binding and molecular recognition, and our results suggest its importance for the docking of small compounds to the artificial cavity. Fourth, after validation against the available experimental data, molecular dynamics simulations can be used to characterize the dynamics and exchange of water molecules and ions, providing atomic level and time-dependent information otherwise inaccessible to experiments.     

The fast gating mechanism in ClC-0 channels

We investigate and then modify the hypothesis that a glutamate side chain acts as the fast gate in ClC-0 channels. We first create a putative open-state configuration of the prokaryotic ClC Cl- channel using its crystallographic structure as a basis. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to ClC-0 by replacing all the nonconserved polar and charged residues. Using this open-state channel model, we carry out molecular dynamics simulations to study how the glutamate side chain can move between open and closed configurations. When the side chain extends toward the extracellular end of the channel, it presents an electrostatic barrier to Cl- conduction. However, external Cl- ions can push the side chain into a more central position where, pressed against the channel wall, it does not impede the motion of Cl- ions. Additionally, a proton from a low-pH external solution can neutralize the extended glutamate side chain, which also removes the barrier to conduction. Finally, we use Brownian dynamics simulations to demonstrate the influence of membrane potential and external Cl- concentration on channel open probability.     

Microscopic Kinetics of DNA Translocation through synthetic nanopores

We have previously demonstrated that a nanometer-diameter pore in a nanometer-thick metal-oxide-semiconductor-compatible membrane can be used as a molecular sensor for detecting DNA. The prospects for using this type of device for sequencing DNA are avidly being pursued. The key attribute of the sensor is the electric field-induced (voltage-driven) translocation of the DNA molecule in an electrolytic solution across the membrane through the nanopore. To complement ongoing experimental studies developing such pores and measuring signals in response to the presence of DNA, we conducted molecular dynamics simulations of DNA translocation through the nanopore. A typical simulated system included a patch of a silicon nitride membrane dividing water solution of potassium chloride into two compartments connected by the nanopore. External electrical fields induced capturing of the DNA molecules by the pore from the solution and subsequent translocation. Molecular dynamics simulations suggest that 20-basepair segments of double-stranded DNA can transit a nanopore of 2.2 x 2.6 nm(2) cross section in a few microseconds at typical electrical fields. Hydrophobic interactions between DNA bases and the pore surface can slow down translocation of single-stranded DNA and might favor unzipping of double-stranded DNA inside the pore. DNA occluding the pore mouth blocks the electrolytic current through the pore; these current blockades were found to have the same magnitude as the blockade observed when DNA transits the pore. The feasibility of using molecular dynamics simulations to relate the level of the blocked ionic current to the sequence of DNA was investigated.     

Electrostatics of the intracellular vestibule of K+ channels

Previous calculations using continuum electrostatic calculations showed that a fully hydrated monovalent cation is electrostatically stabilized at the center of the cavity of the KcsA potassium channel. Further analysis demonstrated that this cavity stabilization was controlled by a balance between the unfavorable reaction field due to the finite size of the cavity and the favorable electrostatic field arising from the pore helices. In the present study, continuum electrostatic calculations are used to investigate how the stability of an ion in the intracellular vestibular cavity common to known potassium channels is affected as the inner channel gate opens and the cavity becomes larger and contiguous with the intracellular solution. The X-ray structure of the calcium-activated potassium channel MthK, which was crystallized in the open state, is used to construct models of the KcsA channel in the open state. It is found that, as the channel opens, the barrier at the helix bundle crossing decreases to approximately 0 kcal/mol, but that the ion in the cavity is also significantly destabilized. The results are compared and contrasted with additional calculations performed on the KvAP (voltage-activated) and KirBac1.1 (inward rectifier) channels, as well as models of the pore domain of Shaker in the open and closed state. In conclusion, electrostatic factors give rise to energetic constraints on ion permeation that have important functional consequences on the various K+ channels, and partly explain the presence or absence of charged residues near the inner vestibular entry.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

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.     

Molecular dynamics study of the KcsA channel at 2.0-A resolution: stability and concerted motions within the pore

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-A resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K+ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K+ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K+ ions and water molecules displacements shows concerted file motions on the simulated time scale (approximately 1 ns), which could act as precursor to the diffusion of K+ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.     

Identification of the Molecular Site of Ivabradine Binding to HCN4 Channels

Ivabradine is a specific heart rate-reducing agent approved as a treatment of chronic stable angina. Its mode of action involves a selective and specific block of HCN channels, the molecular components of sinoatrial "funny" (f)-channels. Different studies suggest that the binding site of ivabradine is located in the inner vestibule of HCN channels, but the molecular details of ivabradine binding are unknown. We thus sought to investigate by mutagenesis and in silico analysis which residues of the HCN4 channel, the HCN isoform expressed in the sinoatrial node, are involved in the binding of ivabradine. Using homology modeling, we verified the presence of an inner cavity below the channel pore and identified residues lining the cavity; these residues were replaced with alanine (or valine) either alone or in combination, and WT and mutant channels were expressed in HEK293 cells. Comparison of the block efficiency of mutant vs WT channels, measured by patch-clamp, revealed that residues Y506, F509 and I510 are involved in ivabradine binding. For each mutant channel, docking simulations correctly explain the reduced block efficiency in terms of proportionally reduced affinity for ivabradine binding. In summary our study shows that ivabradine occupies a cavity below the channel pore, and identifies specific residues facing this cavity that interact and stabilize the ivabradine molecule. This study provides an interpretation of known properties of f/HCN4 channel block by ivabradine such as the “open channel block”, the current-dependence of block and the property of "trapping" of drug molecules in the closed configuration.