On the potential functions used in molecular dynamics simulations of ion channels

The Persistent Question of Potassium Channel Permeation Mechanisms

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K⁺ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K⁺ ions by the channel. However, the molecular mechanism of K⁺ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K⁺ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

An efficient computational approach for modeling protein electrostatic is developed according to static point-charge model distributions based on the linear-scaling EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) quantum mechanical (QM) method. In this approach, the Electrostatic-Potential atomic charges are obtained from ab initio calculation of protein, both polarization and charge transfer effect are taken into consideration. This approach shows a significant improvement in the description of electrostatic potential and solvation energy of proteins comparing with current popular molecular mechanics (MM) force fields. Therefore, it has gorgeous prospect in many applications, including accurate calculations of electric field or vibrational Stark spectroscopy in proteins and predicting protein-ligand binding affinity. It can also be applied in QM/MM calculations or electronic embedding method of ONIOM to provide a better electrostatic environment.

Dynamic water patterns change the stability of the collapsed filter conformation of the KcsA K+ channel

The selectivity filter of the KcsA K+ channel has two typical conformations-the conductive and the collapsed conformations, respectively. The transition from the conductive to the collapsed filter conformation can represent the process of inactivation that depends on many environmental factors. Water molecules permeating behind the filter can influence the collapsed filter stability. Here we perform the molecular dynamics simulations to study the stability of the collapsed filter of the KcsA K+ channel under the different water patterns. We find that the water patterns are dynamic behind the collapsed filter and the filter stability increases with the increasing number of water molecules. In addition, the stability increases significantly when water molecules distribute uniformly behind the four monomeric filter chains, and the stability is compromised if water molecules only cluster behind one or two adjacent filter chains. The altered filter stabilities thus suggest that the collapsed filter can inactivate gradually under the dynamic water patterns. We also demonstrate how the different water patterns affect the filter recovery from the collapsed conformation.

Mechanism of activation at the selectivity filter of the KcsA K+ channel

Potassium channels are opened by ligands and/or membrane potential. In voltage-gated K⁺ channels and the prokaryotic KcsA channel, conduction is believed to result from opening of an intracellular constriction that prevents ion entry into the pore. On the other hand, numerous ligand-gated K⁺ channels lack such gate, suggesting that they may be activated by a change within the selectivity filter, a narrow region at the extracellular side of the pore. Using molecular dynamics simulations and electrophysiology measurements, we show that ligand-induced conformational changes in the KcsA channel removes steric restraints at the selectivity filter, thus resulting in structural fluctuations, reduced K⁺ affinity, and increased ion permeation. Such activation of the selectivity filter may be a universal gating mechanism within K⁺ channels. The occlusion of the pore at the level of the intracellular gate appears to be secondary.

Potassium channels are proteins found in almost all living organisms and are vital for many different biological processes. These proteins contain a pore that allows potassium ions to flow through cell membranes, but only when the channel is open. Most channels have a narrowing at the inward side of the pore, which was proposed to form a gate that controls the flow of ions. This gate only opens when the channel activates. Nearer the outward side of the pore is another narrow region called the selectivity filter. This region interacts selectively with potassium ions as they pass through the channel.

Not all channels form a tight constriction at the inward end of their pore, yet they still only allow potassium ions to flow through when activated. This suggested that these channels might instead use the selectivity filter as a gate. It would also mean that whilst different potassium channels have similar structures, they do not share a common gating mechanism that controls the flow of ions.

Heer et al. studied the bacterial channel called KcsA, which was thought to control the flow of potassium ions via a traditional inward gate. Computer simulations based on this protein’s structure, and experiments with purified KcsA in artificial membranes, showed that the selectivity filter was also involved when KcsA was activated. When the activated channel changed shape, the selectivity filter – which had been constrained by regions forming the pore – could now move and allow potassium ions to flow through the pore. Heer et al. confirmed using a mutant KcsA that the motion of the inward side of the pore upon activation affected the movement of the selectivity filter gate as predicted by the simulation.

These findings show that the KcsA channel opens and closes at the selectivity filter. The changes at the inward side of the pore, previously believed to be the gate, firstly enable the selectivity filter to serve as a gate while also forming a secondary gate. Heer et al. propose that most if not all potassium channels may also use this mechanism. These findings illustrate that molecular simulations can be powerful for predicting how changes in the structure of a protein will affect its behavior. This and future studies of other potassium channels will help scientists to better understand the subtle differences between the diverse range of channels found across different organisms.

Three pairs of weak interactions precisely regulate the G-loop gate of Kir2.1 channel

Kir2.1 (also known as IRK1) plays key roles in regulation of resting membrane potential and cell excitability. To achieve its physiological roles, Kir2.1 performs a series of conformational transition, named as gating. However, the structural basis of gating is still obscure. Here, we combined site-directed mutation, two-electrode voltage clamp with molecular dynamics simulations and determined that H221 regulates the gating process of Kir2.1 by involving a weak interaction network. Our data show that the H221R mutant accelerates the rundown kinetics and decelerates the reactivation kinetics of Kir2.1. Compared with the WT channel, the H221R mutation strengthens the interaction between the CD- and G-loops (E303-R221) which stabilizes the close state of the G-loop gate and weakens the interactions between C-linker and CD-loop (R221-R189) and the adjacent G-loops (E303-R312) which destabilizes the open state of G-loop gate. Our data indicate that the three pairs of interactions (E303-H221, H221-R189 and E303-R312) precisely regulate the G-loop gate by controlling the conformation of G-loop. Proteins 2016; 84:1929-1937. © 2016 Wiley Periodicals, Inc.

Modeling Electronic Polarizability Changes in the Course of a Magnesium Ion Water Ligand Exchange Process

This paper introduces explicit dependence of atomic polarizabilities on intermolecular interactions within the framework of a polarizable force field AMOEBA. Polarizable models used in biomolecular simulations often poorly describe molecular electrostatic induction in condensed phase, in part, due to neglect of a strong dependency of molecular electronic polarizability on intermolecular interactions at short distances. Our variable polarizability model parameters are derived from quantum chemical calculations of small clusters of atoms and molecules, and can be applied in simulations in condensed phase without additional scaling factors. The variable polarizability model is applied to simulate a ligand exchange reaction for a Mg(2+) ion solvated in water. Explicit dependence of water polarizability on a distance between a water oxygen and Mg(2+) is derived from in vacuum MP2 calculations of Mg(2+)-water dimer. The simulations yield a consistent description of the energetics of the Mg(2+)-water clusters of different size. Simulations also reproduce thermodynamics of ion solvation as well as kinetics of a water ligand exchange reaction. In contrast, simulations that used the additive force field or that used the constant polarizability models were not able to consistently and quantitatively describe the properties of the solvated Mg(2+) ion.

Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides

Chondroitin sulfate (CS) is one of several glycosaminoglycans that are major components of proteoglycans. A linear polymer consisting of repeats of the disaccharide -4GlcAβ1-3GalNAcβ1-, CS undergoes differential sulfation resulting in five unique sulfation patterns. Because of the dimer repeat, the CS glycosidic "backbone" has two distinct sets of conformational degrees of freedom defined by pairs of dihedral angles: (ϕ1, ψ1) about the β1-3 glycosidic linkage and (ϕ2, ψ2) about the β1-4 glycosidic linkage. Differential sulfation and the possibility of cation binding, combined with the conformational flexibility and biological diversity of CS, complicate experimental efforts to understand CS three-dimensional structures at atomic resolution. Therefore, all-atom explicit-solvent molecular dynamics simulations with Adaptive Biasing Force sampling of the CS backbone were applied to obtain high-resolution, high-precision free energies of CS disaccharides as a function of all possible backbone geometries. All 10 disaccharides (β1-3 vs β1-4 linkage × five different sulfation patterns) were studied; additionally, ion effects were investigated by considering each disaccharide in the presence of either neutralizing sodium or calcium cations. GlcAβ1-3GalNAc disaccharides have a single, broad, thermodynamically important free-energy minimum, whereas GalNAcβ1-4GlcA disaccharides have two such minima. Calcium cations but not sodium cations bind to the disaccharides, and binding is primarily to the GlcA -COO(-) moiety as opposed to sulfate groups. This binding alters the glycan backbone thermodynamics in instances where a calcium cation bound to -COO(-) can act to bridge and stabilize an interaction with an adjacent sulfate group, whereas, in the absence of this cation, the proximity of a sulfate group to -COO(-) results in two like charges being both desolvated and placed adjacent to each other and is found to be destabilizing. In addition to providing information on sulfation and cation effects, the present results can be applied to building models of CS polymers and as a point of comparison in studies of CS polymer backbone dynamics and thermodynamics.

Ion permeation in K⁺ channels occurs by direct Coulomb knock-on

Potassium channels selectively conduct K(+) ions across cellular membranes with extraordinary efficiency. Their selectivity filter exhibits four binding sites with approximately equal electron density in crystal structures with high K(+) concentrations, previously thought to reflect a superposition of alternating ion- and water-occupied states. Consequently, cotranslocation of ions with water has become a widely accepted ion conduction mechanism for potassium channels. By analyzing more than 1300 permeation events from molecular dynamics simulations at physiological voltages, we observed instead that permeation occurs via ion-ion contacts between neighboring K(+) ions. Coulomb repulsion between adjacent ions is found to be the key to high-efficiency K(+) conduction. Crystallographic data are consistent with directly neighboring K(+) ions in the selectivity filter, and our model offers an intuitive explanation for the high throughput rates of K(+) channels.

Initial steps of inactivation at the K+ channel selectivity filter

K(+) efflux through K(+) channels can be controlled by C-type inactivation, which is thought to arise from a conformational change near the channel's selectivity filter. Inactivation is modulated by ion binding near the selectivity filter; however, the molecular forces that initiate inactivation remain unclear. We probe these driving forces by electrophysiology and molecular simulation of MthK, a prototypical K(+) channel. Either Mg(2+) or Ca(2+) can reduce K(+) efflux through MthK channels. However, Ca(2+), but not Mg(2+), can enhance entry to the inactivated state. Molecular simulations illustrate that, in the MthK pore, Ca(2+) ions can partially dehydrate, enabling selective accessibility of Ca(2+) to a site at the entry to the selectivity filter. Ca(2+) binding at the site interacts with K(+) ions in the selectivity filter, facilitating a conformational change within the filter and subsequent inactivation. These results support an ionic mechanism that precedes changes in channel conformation to initiate inactivation.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Anharmonic vibrations of N-H in Cl(-)(N-methylacetamide)1(H2O)(0-2)Ar2 cluster ions. Combined IRPD experiments and BOMD simulations

Infrared Predissociation (IRPD) spectra of Cl(-)(NMA)1(H2O)0-2Ar2 combined with Born-Oppenheimer Molecular Dynamics (BOMD) IR spectra have been acquired, providing the structure and dynamics of these systems. We show that the chloride ion is bound to the hydrogen of the amide N-H group, forming a strong ionic hydrogen bond, weakening the N-H stretch, and shifting it to lower frequency. The presence of water molecules enhances the ionic hydrogen bond by binding to the amide carbonyl oxygen of NMA and shifts the N-H stretch further to lower frequency. The BOMD IR spectra can recapture all, but about 100 cm(-1), of the 600 to 700 cm(-1) shifts due to the strong N-H stretch anharmonicities observed in experiments. This residual error was found to be due to the lack of zero point energy in the classical treatment of motion in the BOMD method.

Simulations of anionic lipid membranes: development of interaction-specific ion parameters and validation using NMR data

Overbinding of ions to lipid head groups is a potentially serious artifact in simulations of charged lipid bilayers. In this study, the Lennard-Jones radii in the CHARMM force field for interactions of Na(+) and lipid oxygen atoms of carboxyl, phosphate, and ester groups were revised to match osmotic pressure data on sodium acetate and electrophoresis data on palmitoyloleoyl phosphatidylcholine (POPC) vesicles. The new parameters were then validated by successfully reproducing previously published experimental NMR deuterium order parameters for dimyristoyl phosphatidylglycerol (DMPG) and newly obtained values for palmitoyloleoyl phosphatidylserine (POPS). Although the increases in Lennard-Jones diameters are only 0.02-0.12 Å, they are sufficient to reduce Na+ binding, and thereby increase surface areas per lipid by 5-10% compared with the unmodified parameters.

Atomic-level simulation of current-voltage relationships in single-file ion channels

The difficulty in characterizing ion conduction through membrane channels at the level of individual permeation events has made it challenging to elucidate the mechanistic principles underpinning this fundamental physiological process. Using long, all-atom simulations enabled by special-purpose hardware, we studied K⁺ permeation across the K_V1.2/2.1 voltage-gated potassium channel. At experimentally accessible voltages, which include the physiological range, the simulated permeation rate was substantially lower than the experimentally observed rate. The current–voltage relationship was also nonlinear but became linear at much higher voltages. We observed permeation consistent with a “knock-on” mechanism at all voltages. At high voltages, the permeation rate was in accordance with our previously reported K_V1.2 pore-only simulations, after the simulated voltages from the previous study were recalculated using the correct method, new insight into which is provided here. Including the voltage-sensing domains in the simulated channel brought the linear current–voltage regime closer to the experimentally accessible voltages. The simulated permeation rate, however, still underestimated the experimental rate, because formation of the knock-on intermediate occurred too infrequently. Reducing the interaction strength between the ion and the selectivity filter did not increase conductance. In complementary simulations of gramicidin A, similar changes in interaction strength did increase the observed permeation rate. Permeation nevertheless remained substantially below the experimental value, largely because of infrequent ion recruitment into the pore lumen. Despite the need to apply large voltages to simulate the permeation process, the apparent voltage insensitivity of the permeation mechanism suggests that the direct simulation of permeation at the single-ion level can provide fundamental physiological insight into ion channel function. Notably, our simulations suggest that the knock-on permeation mechanisms in K_V1.2 and KcsA may be different.

Self-Learning Adaptive Umbrella Sampling Method for the Determination of Free Energy Landscapes in Multiple Dimensions

The potential of mean force describing conformational changes of biomolecules is a central quantity that determines the function of biomolecular systems. Calculating an energy landscape of a process that depends on three or more reaction coordinates might require a lot of computational power, making some of multidimensional calculations practically impossible. Here, we present an efficient automatized umbrella sampling strategy for calculating multidimensional potential of mean force. The method progressively learns by itself, through a feedback mechanism, which regions of a multidimensional space are worth exploring and automatically generates a set of umbrella sampling windows that is adapted to the system. The self-learning adaptive umbrella sampling method is first explained with illustrative examples based on simplified reduced model systems, and then applied to two non-trivial situations: the conformational equilibrium of the pentapeptide Met-enkephalin in solution and ion permeation in the KcsA potassium channel. With this method, it is demonstrated that a significant smaller number of umbrella windows needs to be employed to characterize the free energy landscape over the most relevant regions without any loss in accuracy.

Nonselective conduction in a mutated NaK channel with three cation-binding sites

The NaK channel is a cation-selective protein with similar permeability for K(+) and Na(+) ions. Crystallographic structures are available for the wild-type and mutated NaK channels with different numbers of cation-binding sites. We have performed a comparison between the potentials of mean force governing the translocation of K(+) ions and mixtures of one Na(+) and three K(+) ions in a mutated NaK channel with only three cation-binding sites (NaK-CNG). Since NaK-CNG is not selective for K(+) over Na(+), analysis of its multi-ion potential energy surfaces can provide clues about how selectivity originates. Comparison of the potentials of mean force of NaK-CNG and K(+)-selective channels yields observations that strongly suggest that the number of contiguous ion binding sites in a single-file mechanism is the key determinant of the channel's selectivity properties, as already proposed by experimental studies. We conclude that the presence of four binding sites in K(+)-selective channels is essential for highly selective and efficient permeation of K(+) ions, and that a key difference between K(+)-selective and nonselective channels is the absence/presence of a binding site for Na(+) ions at the boundary between S2 and S3 in the context of multi-ion permeation events.

A molecular simulation based assessment of binding of metal ions on micelles

An assessment of the ability of a micellar surface to bind different metal ions using molecular simulation is presented in this study. Sodium dodecyl sulfate (SDS) is considered as the anionic surfactant. Various relevant characteristics of SDS-metal ion systems are estimated to quantify preferential binding of metal ions. These are electrostatic energy, total potential energy of the system, radial distribution function, and entropy and free energy change of the system. By examining these parameters, the relative extents of binding of different metal ions to the micellar surface are assessed.

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.

Molecular dynamics simulations of the TrkH membrane protein

TrkH is a transmembrane protein that mediates uptake of K(+) through the cell membrane. Despite the recent determination of its crystallographic structure, the nature of the permeation mechanism is still unknown, that is, whether K(+) ions move across TrkH by active transport or passive diffusion. Here, molecular dynamics simulations and the umbrella sampling technique have been employed to shed light on this question. The existence of binding site S3 and two alternative binding sites have been characterized. Analysis of the coordination number renders values that are almost constant, with a full contribution from the carbonyls of the protein only at S3. This observation contrasts with observations of K(+) channels, where the contribution of the protein to the coordination number is roughly constant in all four binding sites. An intramembrane loop is found immediately after the selectivity filter at the intracellular side of the protein, which obstructs the permeation pathway, and this is reflected in the magnitude of the energy barriers.

Gating at the selectivity filter of ion channels that conduct Na+ and K+ ions

The NaK channel is a cation selective channel with similar permeability for K(+) and Na(+). The available crystallographic structure of wild-type (WT) NaK is usually associated with a conductive state of the channel. Here, potential of mean force for complete conduction events of Na(+) and K(+) ions through NaK show that: i), large energy barriers prevent the passage of ions through the WT NaK structure, ii), the barriers are correlated to the presence of a hydrogen bond between Asp-66 and Asn-68, and iii), the structure of NaK mutated to mimic cyclic nucleotide-gated channels conducts Na(+) and K(+). These results support the hypothesis that the filter of cation selective channels can adopt at least two different structures: a conductive one, represented by the x-ray structures of the NaK-CNG chimeras, and a closed one, represented by the x-ray structures of the WT NaK.

Absence of ion-binding affinity in the putatively inactivated low-[K+] structure of the KcsA potassium channel

Potassium channels are membrane proteins that selectively conduct K(+) across cellular membranes. The narrowest part of their pore, the selectivity filter, is responsible for distinguishing K(+) from Na(+), and can also act as a gate through a mechanism known as C-type inactivation. It has been proposed that a conformation of the KcsA channel obtained by crystallization in presence of low concentration of K(+) (PDB 1K4D) could correspond to the C-type inactivated state. Here, we show using molecular mechanics simulations that such conformation has little ion-binding affinity and that ions do not contribute to its stability. The simulations suggest that, in this conformation, the selectivity filter is mostly occupied by water molecules. Whether such ion-free state of the KcsA channel is physiologically accessible and representative of the inactivated state of eukaryotic channels remains unclear.

A combined experimental and theoretical study of ion solvation in liquid N-methylacetamide

Most current biomolecular simulations are based on potential energy functions that treat the electrostatic energy as a sum of pairwise Coulombic interactions between effective fixed atomic charges. This approximation, in which many-body induced polarization effects are included in an average way, is expected to be satisfactory for a wide range of systems, but less accurate for processes involving the transfer and partition of ions among heterogeneous environments. The limitations of these potential energy functions are perhaps most obvious in studies of ion permeation through membrane channels. In many cases, the pore is so narrow that the permeating ion must shed most of its surrounding water molecules and the large energetic loss due to dehydration must be compensated by coordination with protein atoms. Interactions of cations with protein backbone carbonyl oxygens, in particular, play a critical role in several important biological channels. As a first step toward meeting the challenge of developing an accurate explicit accounting for induced polarization effects, the present work combines experiments and computation to characterize the interactions of alkali and halide ions with N-methylacetamide chosen to represent the peptide bond. From solubility measurements, we extract the solvation free energies of KCl and NaCl in liquid N-methylacetamide. Polarizable models based on the Drude oscillator are then developed and compared with available experimental and ab initio data. The good agreement for a range of structural and thermodynamic properties in the gas and condensed phases suggests that the polarizable models provide an accurate representation of ion-amide interactions in biological systems.

Ion selectivity of the KcsA channel: a perspective from multi-ion free energy landscapes

Potassium (K(+)) channels are specialized membrane proteins that are able to facilitate and regulate the conduction of K(+) through cell membranes. Comprising five specific cation binding sites (S(0)-S(4)) formed by the backbone carbonyl groups of conserved residues common to all K(+) channels, the narrow selectivity filter allows fast conduction of K(+) while being highly selective for K(+) over Na(+). To extend our knowledge of the microscopic mechanism underlying selectivity in K(+) channels, we characterize the free energy landscapes governing the entry and translocation of a Na(+) or a K(+) from the extracellular side into the selectivity filter of KcsA. The entry process of an extracellular ion is examined in the presence of two additional K(+) in the pore, and the three-ion potential of mean force is computed using extensive all-atom umbrella sampling molecular dynamics simulations. A comparison of the potentials of mean force yields a number of important results. First, the free energy minima corresponding to configurations with extracellular K(+) or Na(+) in binding site S(0) or S(1) are similar in depth, suggesting that the thermodynamic selectivity governed by the free energy minima for those two binding sites is insignificant. Second, the free energy barriers between stable multi-ion configurations are generally higher for Na(+) than for K(+), implying that the kinetics of ion conduction is slower when a Na(+) enters the pore. Third, the region corresponding to binding site S(2) near the center of the narrow pore emerges as the most selective for K(+) over Na(+). In particular, while there is a stable minimum for K(+) in site S(2), Na(+) faces a steep free energy increase with no local free energy well in this region. Lastly, analysis shows that selectivity is not correlated with the overall coordination number of the ion entering the pore, but is predominantly affected by changes in the type of coordinating ligands (carbonyls versus water molecules). These results further highlight the importance of the central region near binding site S(2) in the selectivity filter of K(+) channels.

On the utilization of energy minimization to the study of ion selectivity

The major pitfalls in studying ion selectivity in binding site models using energy minimization based methods are examined and discussed. It is shown that the properties derived from energy minimization are strongly configuration-dependent and that the results should be interpreted with caution. It is concluded that computational studies of ion selectivity must include thermal fluctuations and entropic effects.

Atypical mechanism of conduction in potassium channels

Potassium channels can conduct passively K+ ions with rates of up to approximately 10(8) ions per second at physiological conditions, and they are selective to these species by a factor of 10(4) over Na+ ions. Ion conduction has been proposed to involve transitions between 2 main states, with 2 or 3 K+ ions occupying the selectivity filter separated by an intervening water molecule. The largest free energy barrier of such a process was reported to be of the order of 2-3 kcal mol(-1). Here, we present an alternative mechanism for conduction of K+ in potassium channels where site vacancies are involved, and we propose that coexistence of several ion permeation mechanisms is energetically possible. Conduction can be described as a more anarchic phenomenon than previously characterized by the concerted translocations of K+-water-K+.

Hydration number, topological control, and ion selectivity

The topological control hypothesis presented by Bostick and Brooks [Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 9260] has sought to explain binding selectivity in potassium channels based on the premise that a universal measure of ion solvation in different environments is provided by its average coordination structure in bulk water. This leads to the view that ion selectivity is predominantly controlled by the number of ligands coordinating the ion and that the chemical type of those ligands plays a minor role. The significance of the topological control hypothesis and its ability to predict ion selectivity in protein binding sites are examined. It is shown that the framework encounters increasing difficulties when different protein binding sites with similar coordination numbers are considered.

Dynamics, energetics, and selectivity of the low-K+ KcsA channel structure

Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure.

Many-body effects and simulations of potassium channels

The electronic polarizability of an ion or a molecule is a measure of the relative tendency of its electron cloud to be distorted from its normal shape by an electric field. On the molecular scale, in a condensed phase, any species sits in an electric field due to its neighbours, and the resulting polarization is an important contribution to the total interaction energy. Electrostatic interactions are crucial for determining the majority of chemical–physical properties of the system and electronic polarization is a fundamental component of these interactions. Thus, polarization effects should be taken into account if accurate descriptions are desired. In classical computer simulations, the forces required to drive the system are typically based on interatomic interaction potentials derived in part from electronic structure calculations or from experimental data. Owing to the difficulties in including polarization effects in classical force fields, most of them are based just on pairwise additive interaction potentials. At present, major efforts are underway to develop polarizable interaction potentials for biomolecular simulations. In this review, various ways of introducing explicit polarizability into biomolecular models and force fields are reviewed, and the progress that might be achieved in applying such methods to study potassium channels is described.

Continuum polarizable force field within the Poisson-Boltzmann framework

We have developed and tested a complete set of nonbonded parameters for a continuum polarizable force field. Our analysis shows that the new continuum polarizable model is consistent with B3LYP/cc-pVTZ in modeling electronic response upon variation of dielectric environment. Comparison with experiment also shows that the new continuum polarizable model is reasonable, with accuracy similar to that of B3LYP/cc-pVTZ in reproduction of dipole moments of selected organic molecules in the gas phase. We have further tested the validity to interchange the Amber van der Waals parameters between the explicit and continuum polarizable force fields with a series of dimers. It can be found that the continuum polarizable model agrees well with MP2/cc-pVTZ, with deviations in dimer binding energies less than 0.9 kcal/mol in the aqueous dielectric environment. Finally, we have optimized atomic cavity radii with respect to experimental solvation free energies of 177 training molecules. To validate the optimized cavity radii, we have tested these parameters against 176 test molecules. It is found that the optimized Poisson-Boltzmann atomic cavity radii transfer well from the training set to the test set, with an overall root-mean-square deviation of 1.30 kcal/mol, an unsigned average error of 1.07 kcal/mol, and a correlation coefficient of 92% for all 353 molecules in both the training and test sets. Given the development documented here, the next natural step is the construction of a full protein/nucleic acid force field within the new continuum polarization framework.

Molecular dynamics simulations of potassium channels

Despite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.

A theoretical study of aqueous solvation of K comparing ab initio, polarizable, and fixed-charge models

The hydration of K(+) is studied using a hierarchy of theoretical approaches, including ab initio Born-Oppenheimer molecular dynamics and Car-Parrinello molecular dynamics, a polarizable force field model based on classical Drude oscillators, and a nonpolarizable fixed-charge potential based on the TIP3P water model. While models based more directly on quantum mechanics offer the possibility to account for complex electronic effects, polarizable and fixed-charges force fields allow for simulations of large systems and the calculation of thermodynamic observables with relatively modest computational costs. A particular emphasis is placed on investigating the sensitivity of the polarizable model to reproduce key aspects of aqueous K(+), such as the coordination structure, the bulk hydration free energy, and the self diffusion of K(+). It is generally found that, while the simple functional form of the polarizable Drude model imposes some restrictions on the range of properties that can simultaneously be fitted, the resulting hydration structure for aqueous K(+) agrees well with experiment and with more sophisticated computational models. A counterintuitive result, seen in Car-Parrinello molecular dynamics and in simulations with the Drude polarizable force field, is that the average induced molecular dipole of the water molecules within the first hydration shell around K(+) is slightly smaller than the corresponding value in the bulk. In final analysis, the perspective of K(+) hydration emerging from the various computational models is broadly consistent with experimental data, though at a finer level there remain a number of issues that should be resolved to further our ability in modeling ion hydration accurately.

Polarization effects and charge transfer in the KcsA potassium channel

The electronic structure of the selectivity filter of KcsA K(+) channel is investigated by density functional theory (DFT/BLYP) and QM/MM methods. The quantum part includes the selectivity filter, which is polarized by the electrostatic field of the environment treated with the Amber force field. The details of the electronic structure were investigated using the maximally localized Wannier function centers of charge and Bader's atoms in molecules charge analysis. Our results show that the channel backbone carbonyl groups are largely polarized and that there is a sizeable charge transfer from the backbone to the cations. These effects are expected to be important for an accurate description of the carbonyl groups and the ion-ion electrostatic repulsion, which have been proposed to play a central role for the selectivity mechanism of the channel [S.Y. Noskov, S. Berneche, B. Roux, Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands. Nature 431 (2004) 830-834].

Ion selectivity in potassium channels

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.

Ion permeation through a narrow channel: using gramicidin to ascertain all-atom molecular dynamics potential of mean force methodology and biomolecular force fields

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.

Poisson–Nernst–Planck Theory Approach to the Calculation of Current Through Biological Ion Channels

The Poisson-Nernst-Planck (PNP) theory of electro-diffusion is reviewed. Techniques for numerical solution of the three-dimensional PNP equations are summarized, and several illustrative applications to ion transport through protein channels are presented. Strengths and weaknesses of the theory are discussed, as well as attempts to improve it via increasingly realistic evaluation of the force acting on each ion due to the protein/membrane environment.

Dependence of ion hydration on the sign of the ion’s charge

The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.

Energetics of ion conduction through the gramicidin channel

The free energy governing K(+) conduction through gramicidin A channels is characterized by using over 0.1 micros of all-atom molecular dynamics simulations with explicit solvent and membrane. The results provide encouraging agreement with experiments and insights into the permeation mechanism. The free energy surface of K(+), as a function of both axial and radial coordinates, is calculated. Correcting for simulation artifacts due to periodicity and the lack of hydrocarbon polarizability, the calculated single-channel conductance for K(+) ions is 0.8 pS, closer to experiment than any previous calculation. In addition, the estimated single ion dissociation constants are within the range of experimental determinations. The relatively small free energy barrier to ion translocation arises from a balance of large opposing contributions from protein, single-file water, bulk electrolyte, and membrane. Mean force decomposition reveals a remarkable ability of the single-file water molecules to stabilize K(+) by -40 kcal/mol, roughly half the bulk solvation free energy. The importance of the single-file water confirms the conjecture of Mackay et al. [Mackay, D. H. J., Berens, P. H., Wilson, K. R. & Hagler, A. T. (1984) Biophys. J. 46, 229-248]. Ion association with the channel involves gradual dehydration from approximately six to seven water molecules in the first shell, to just two inside the narrow pore. Ion permeation is influenced by the orientation of the single-file water column, which can present a barrier to conduction and give rise to long-range coupling of ions on either side of the pore. Small changes in the potential function, including contributions from electronic polarization, are likely to be sufficient to obtain quantitative agreement with experiments.