Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

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.     

Molecular restraints in the permeation pathway of ion channels

Ion channels assist and control the diffusion of ions through biological membranes. The conduction process depends on the structural characteristics of these nanopores, among which are the hydrophobicity and the afforded diameter of the conduction pathway. In this contribution, we use full atomistic free-energy molecular dynamics simulations to estimate the effect of such characteristics on the energetics of ion conduction through the activation gate of voltage-gated potassium (Kv) channels. We consider specifically the ionic translocation through three different permeation pathways, corresponding to the activation gate of an atomistic model of Shaker channels in closed and partially opened conformations, and that of the open conformation of the Kv1.2 channel. In agreement with experiments, we find that the region of Val(478) constitutes the main gate. The conduction is unfavorable through this gate when the constriction is smaller than an estimated threshold of 4.5-5.0 A, mainly due to incomplete coordination-hydration of the ion. Above this critical size, e.g., for the Kv1.2, the valine gate is wide enough to allow fully coordination of the ion and therefore its diffusion on a flat energy surface. Similar to other ion channels, Kv channels appear therefore to regulate diffusion by constricting hydrophobic regions of the permeation pathway.     

CHAP: A Versatile Tool for the Structural and Functional Annotation of Ion Channel Pores

The control of ion channel permeation requires the modulation of energetic barriers or “gates” within their pores. However, such barriers are often simply identified from the physical dimensions of the pore. Such approaches have worked well in the past, but there is now evidence that the unusual behavior of water within narrow hydrophobic pores can produce an energetic barrier to permeation without requiring steric occlusion of the pathway. Many different ion channels have now been shown to exploit “hydrophobic gating” to regulate ion flow, and it is clear that new tools are required for more accurate functional annotation of the increasing number of ion channel structures becoming available. We have previously shown how molecular dynamics simulations of water can be used as a proxy to predict hydrophobic gates, and we now present a new and highly versatile computational tool, the Channel Annotation Package (CHAP) that implements this methodology.     

Water dynamics and dewetting transitions in the small mechanosensitive channel MscS

The dynamics of confined water in capillaries and nanotubes suggests that gating of ion channels may involve not only changes of the pore geometry, but also transitions between water-filled and empty states in certain locations. The recently solved heptameric structure of the small mechanosensitive channel of Escherichia coli, MscS, has revealed a relatively wide (7-15 A) yet highly hydrophobic transmembrane pore. Continuum estimations based on the properties of pore surface suggest low conductance and a thermodynamic possibility of dewetting. To test the predictions we performed molecular dynamics simulations of MscS filled with flexible TIP3P water. Irrespective to the initial conditions, several independent 6-ns simulations converged to the same stable state with the pore water-filled in the wider part, but predominantly empty in the narrow hydrophobic part, displaying intermittent vapor-liquid transitions. The polar gain-of-function substitution L109S in the constriction resulted in a stable hydration of the entire pore. Steered passages of Cl(-) ions through the narrow part of the pore consistently produced partial ion dehydration and required a force of 200-400 pN to overcome an estimated barrier of 10-20 kcal/mole, implying negligibly low conductance. We conclude that the crystal structure of MscS does not represent an open state. We infer that MscS gate, which is similar to that of the nicotinic ACh receptor, involves a vapor-lock mechanism where limited changes of geometry or surface polarity can locally switch the regime between water-filled (conducting) and empty (nonconducting) states.     

A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.     

The influence of geometry, surface character, and flexibility on the permeation of ions and water through biological pores

A hydrophobic constriction site can act as an efficient barrier to ion and water permeation if its diameter is less than the diameter of an ion's first hydration shell. This hydrophobic gating mechanism is thought to operate in a number of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive channels (MscL and MscS) and perhaps in some potassium channels (e.g. KcsA, MthK and KvAP). Simplified pore models allow one to investigate the primary characteristics of a conduction pathway, namely its geometry (shape, pore length, and radius), the chemical character of the pore wall surface, and its local flexibility and surface roughness. Our extended (about 0.1 micros) molecular dynamic simulations show that a short hydrophobic pore is closed to water for radii smaller than 0.45 nm. By increasing the polarity of the pore wall (and thus reducing its hydrophobicity) the transition radius can be decreased until for hydrophilic pores liquid water is stable down to a radius comparable to a water molecule's radius. Ions behave similarly but the transition from conducting to non-conducting pores is even steeper and occurs at a radius of 0.65 nm for hydrophobic pores. The presence of water vapour in a constriction zone indicates a barrier for ion permeation. A thermodynamic model can explain the behaviour of water in nanopores in terms of the surface tensions, which leads to a simple measure of 'hydrophobicity' in this context. Furthermore, increased local flexibility decreases the permeability of polar species. An increase in temperature has the same effect, and we hypothesize that both effects can be explained by a decrease in the effective solvent-surface attraction which in turn leads to an increase in the solvent-wall surface free energy.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC 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.     

Ion passage pathways and thermodynamics of the amphotericin B membrane channel

Amphotericin B is a polyene macrolide antibiotic used to treat systemic fungal infections. Amphotericin B's chemotherapeutic action requires the formation of transmembrane channels, which are known to transmit monovalent ions. We have investigated the ion passage pathways through the pore of a realistic model structure of the channel and computed the associated thermodynamic properties. Our calculations combined the free energy computations using the Poisson equation with a continuum solvent model and the molecular simulations in which solvent molecules were present explicitly. It was found that there are no substantial structural barriers to a single sodium or chloride ion passage. Thermodynamic free energy calculations showed that the path along which the ions prefer to move is off center from the channel's central axis. In accordance with experiments, Monte Carlo molecular simulations established that sodium ions can pass through the pore. When it encounters a chloride anion in the channel, the sodium cation prefers to form a solvent-bridged pair configuration with the anion.     

Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.     

Ion leakage through transient water pores in protein-free lipid membranes driven by transmembrane ionic charge imbalance

We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na(+) and Cl(-) ions leak through a membrane at approximately the same ratio despite the fact that Na(+) ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore "walls," considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K(+) and Cl(-) ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, "water-defect-mediated" pathway without actual formation of a pore. The "water-defect-mediated" leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.     

Structural determinants of proton blockage in aquaporins

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.     

Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA

Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg²⁺ into cells, we examine early events in the relaxation of Mg²⁺ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg²⁺ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.     

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K⁺ ions in various configurations. For both channels, we find that a pair of K⁺ ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K⁺ ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.     

Role of acetylcholine receptor domains in ion selectivity

The nicotinic acetylcholine receptor (nAChR) is a ligand gated ion channel protein, composed of three domains: a transmembrane domain (TM-domain), extracellular domain (EC-domain), and intracellular domain (IC-domain). Due to its biological importance, much experimental and theoretical research has been carried out to explore its mechanisms of gating and selectivity, but there are still many unresolved issues, especially on the ion selectivity. Moreover, most of the previous theoretical work has concentrated on the TM-domain or EC-domain of nAChR, which may be insufficient to understand the entire structure–function relation. In this work, we perform molecular dynamics, Brownian dynamics simulations and continuum electrostatic calculations to investigate the role of different nAChR domains in ion conduction and selectivity. The results show that although both the EC and IC domains contain strong negative charges that create large cation concentrations at either end of the pore, this alone is not sufficient to create the observed cation selectivity and may play a greater role in determining the channel conductance. The presence of cations in the wide regions of the pore can screen out the protein charge allowing anions to enter, meaning that local regions of the TM-domain are most likely responsible for discriminating between ions. These new results complement our understanding about the ion conduction and selectivity mechanism of nAChR.