Permeation of ions across the potassium channel: brownian dynamics studies

The physical mechanisms underlying the transport of ions across a model potassium channel are described. The shape of the model channel corresponds closely to that deduced from crystallography. From electrostatic calculations, we show that an ion permeating the channel, in the absence of any residual charges, encounters an insurmountable energy barrier arising from induced surface charges. Carbonyl groups along the selectivity filter, helix dipoles near the oval chamber, and mouth dipoles near the channel entrances together transform the energy barrier into a deep energy well. Two ions are attracted to this well, and their presence in the channel permits ions to diffuse across it under the influence of an electric field. Using Brownian dynamics simulations, we determine the magnitude of currents flowing across the channel under various conditions. The conductance increases with increasing dipole strength and reaches its maximum rapidly; a further increase in dipole strength causes a steady decrease in the channel conductance. The current also decreases systematically when the effective dielectric constant of the channel is lowered. The conductance with the optimal choice of dipoles reproduces the experimental value when the dielectric constant of the channel is assumed to be 60. The current-voltage relationship obtained with symmetrical solutions is linear when the applied potential is less than approximately 100 mV but deviates from Ohm's law at a higher applied potential. The reversal potentials obtained with asymmetrical solutions are in agreement with those predicted by the Nernst equation. The conductance exhibits the saturation property observed experimentally. We discuss the implications of these findings for the transport of ions across the potassium channels and membrane channels in general.     

Study of ionic currents across a model membrane channel using Brownian dynamics.

Brownian dynamics simulations have been carried out to study ionic currents flowing across a model membrane channel under various conditions. The model channel we use has a cylindrical transmembrane segment that is joined to a catenary vestibule at each side. Two cylindrical reservoirs connected to the channel contain a fixed number of sodium and chloride ions. Under a driving force of 100 mV, the channel is virtually impermeable to sodium ions, owing to the repulsive dielectric force presented to ions by the vestibular wall. When two rings of dipoles, with their negative poles facing the pore lumen, are placed just above and below the constricted channel segment, sodium ions cross the channel. The conductance increases with increasing dipole strength and reaches its maximum rapidly; a further increase in dipole strength does not increase the channel conductance further. When only those ions that acquire a kinetic energy large enough to surmount a barrier are allowed to enter the narrow transmembrane segment, the channel conductance decreases monotonically with the barrier height. This barrier represents those interactions between an ion, water molecules, and the protein wall in the transmembrane segment that are not treated explicitly in the simulation. The conductance obtained from simulations closely matches that obtained from ACh channels when a step potential barrier of 2-3 kTr is placed at the channel neck. The current-voltage relationship obtained with symmetrical solutions is ohmic in the absence of a barrier. The current-voltage curve becomes nonlinear when the 3 kTr barrier is in place. With asymmetrical solutions, the relationship approximates the Goldman equation, with the reversal potential close to that predicted by the Nernst equation. The conductance first increases linearly with concentration and then begins to rise at a slower rate with higher ionic concentration. We discuss the implications of these findings for the transport of ions across the membrane and the structure of ion channels.     

Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations.

The mechanisms underlying transport of ions across the potassium channel are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We first build open-state configurations of the channel with molecular dynamics simulations, by pulling the transmembrane helices outward until the channel attains the desired interior radius. To gain insights into ion permeation, we construct potential energy profiles experienced by an ion traversing the channel in the presence of other resident ions. These profiles reveal that in the absence of an applied field the channel accommodates three potassium ions in a stable equilibrium, two in the selectivity filter and one in the central cavity. In the presence of a driving potential, this three-ion state becomes unstable, and ion permeation across the channel is observed. These qualitative explanations are confirmed by the results of three-dimensional Brownian dynamics simulations. We find that the channel conducts when the ionizable residues near the extracellular entrance are fully charged and those near the intracellular side are partially charged. The conductance increases steeply as the radius of the intracellular mouth of the channel is increased from 2 A to 5 A. Our simulation results reproduce several experimental observations, including the current-voltage curves, conductance-concentration relationships, and outward rectification of currents.     

Transportation behavior of alkali ions through a cell membrane ion channel. A quantum chemical description of a simplified isolated model

Quantum chemical model calculations were carried out for modeling the ion transport through an isolated ion channel of a cell membrane. An isolated part of a natural ion channel was modeled. The model channel was a calixarene derivative, hydrated sodium and potassium ions were the models of the transported ion. The electrostatic potential of the channel and the energy of the channel-ion system were calculated as a function of the alkali ion position. Both attractive and repulsive ion-channel interactions were found. The calculations - namely the dependence of the system energy and the atomic charges of the water molecules with respect to the position of the alkali ion in the channel - revealed the molecular-structural background of the potassium selectivity of this artificial ion channel. It was concluded that the studied ion channel mimics real biological ion channel quite well.     

Simulation of the effect of an external GHz electric field on the potential energy profile of Ca2+ ions in the selectivity filter of the CaVAb channel

Ca_V channels are transmembrane proteins that mediate and regulate ion fluxes across cell membranes, and they are activated in response to action potentials to allow Ca²⁺ influx. Since ion channels are composed of charge or polar groups, an external alternating electric field may affect the ion‐selective membrane transport and the performance of the channel. In this article, we have investigated the effect of an external GHz electric field on the dynamics of calcium ions in the selectivity filter of the Ca_VAb channel. Molecular dynamics (MD) simulations and the potential of mean force (PMF) calculations were carried out, via the umbrella sampling method, to determine the free energy profile of Ca²⁺ ions in the Ca_VAb channels in presence and absence of an external field. Exposing Ca_VAb channel to 1, 2, 3, 4, and 5 GHz electric fields increases the depth of the potential energy well and this may result in an increase in the affinity and strength of Ca²⁺ ions to binding sites in the selectivity filter the channel. This increase of strength of Ca²⁺ ions binding in the selectivity filter may interrupt the mechanism of Ca²⁺ ion conduction, and leads to a reduction of Ca²⁺ ion permeation through the Ca_VAb channel.     

Monte-Carlo Simulation of Molten CsCl Using a ‘Deformation Dipole’ Polarisable Ion Potential

Abstract

We have simulated the structure of molten CsCl at 943 K using a form of interionic potential, based on the deformation dipole model of Hardy and Karo, where ions have both point charges and dipoles. The magnitude of the dipole on an ion is determined by two factors, one being the local electric field and the other the repulsive potential due to overlap with neighbouring ions. Parameters are determined directly from the crystal lattice constant, compressibility and dielectric constants. Good agreement with experimental results is obtained. The potential is compared to the shell model of ionic polarisability which does not produce good agreement with experiment.     

Mechanisms of permeation and selectivity in calcium channels

The mechanisms underlying ion transport and selectivity in calcium channels are examined using electrostatic calculations and Brownian dynamics simulations. We model the channel as a rigid structure with fixed charges in the walls, representing glutamate residues thought to be responsible for ion selectivity. Potential energy profiles obtained from multi-ion electrostatic calculations provide insights into ion permeation and many other observed features of L-type calcium channels. These qualitative explanations are confirmed by the results of Brownian dynamics simulations, which closely reproduce several experimental observations. These include the current-voltage curves, current-concentration relationship, block of monovalent currents by divalent ions, the anomalous mole fraction effect between sodium and calcium ions, attenuation of calcium current by external sodium ions, and the effects of mutating glutamate residues in the amino acid sequence.     

Modeling diverse range of potassium channels with Brownian dynamics

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

Energy barrier presented to ions by the vestibule of the biological membrane channel.

The role of the vestibule in influencing the permeation of ions through biological ion channels is investigated. We derive analytical expressions for the electric potential satisfying Poisson's equation with prolate spheroidal boundary conditions. To allow more realistic geometries we devise an iterative method to calculate the electric potential arising from a fixed charge and an arbitrary dielectric boundary, and confirm that the analytical expressions and iterative method give similar potential values. We then investigate the size of the potential barrier presented to an ion by model vestibules of conical and catenary shapes. The height of the potential barrier increases steeply as an ion enters the vestibule and moves toward the constricted region of the channel. We show that the barrier presented by, for example, a 15 degrees conical vestibule can be canceled by placing dipoles with a total moment of about 50 Debyes near the constricted region of the pore. The selectivity of cations and anions can result from the polarity of charge groups or the orientation of dipoles located near the constricted region of the channel.     

Analytical solutions of Poisson's equation for realistic geometrical shapes of membrane ion channels.

Analytical solutions of Poisson's equations satisfying the Dirichlet boundary conditions for a toroidal dielectric boundary are presented. The electric potential computed anywhere in the toroidal conduit by the analytical method agrees with the value derived from an iterative numerical method. We show that three different channel geometries, namely, bicone, catenary, and toroid, give similar potential profiles as an ion traverses along their central axis. We then examine the effects of dipoles in the toroidal channel wall on the potential profile of ions passing through the channel. The presence of dipoles eliminates the barrier for one polarity of ion, while raising the barrier for ions of the opposite polarity. We also examine how a uniform electric field from an external source is affected by the protein boundary and a mobile charge. The channel distorts the field, reducing it in the vestibules, and enhancing it in the constricted segment. The presence of an ion in one vestibule effectively excludes ions of the same polarity from that vestibule, but has little effect in the other vestibule. Finally, we discuss how the solutions we provide here may be utilized to simulate a system containing a channel and many interacting ions.     

Brownian dynamics simulations of ion transport through the VDAC

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.     

Effect of 3-phenylindole on lipophilic ion and carrier-mediated ion transport across bilayer lipid membranes

The physical effects of 3-phenylindole, an antimicrobial compound which interacts with phospholipids, on ion transport across phosphatidylcholine-cholesterol bilayers have been investigated using three lipophilic ions and one ion-carrier complex. It was found that 3-phenylindole increased membrane electrical conductance of positively charged membrane probes and decreased electrical conductance of negatively charged probes. The enhancement of conductance detected by nonactin-K+ complex and tetraphenylarsonium+ was several orders of magnitude, whereas the suppression of conductance due to tetraphenylborate- and dipicrylamine- was less than a factor of ten. Presence of 3-phenylindole in aqueous phase slightly decreased adsorption of tetraphenylborate- and dipicrylamine- at the membrane surface. From the voltage dependence of the steady-state conductance it was shown that 3-phenylindole induced kinetic limitation of membrane transport of potassium mediated by nonactin. No such limitation was found in the case of tetraphenylarsonium+ transport. These results are shown to be consistent with the present concept of ion diffusion in membranes and the assumption that 3-phenylindole decreases the electric potential in the membrane interior. The asymmetry of the effect of 3-phenylindole on the magnitude of conductance changes for positively and negatively charged membrane permeable ions is also discussed as a reflection of the discreteness of both the absorbed 3-phenylindole and lipid dipoles.     

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.     

Focused electric field across the voltage sensor of potassium channels

Voltage-gated ion channels respond to changes in membrane potential by movement of their voltage sensors across the electric field between cytoplasmic and extracellular solutions. The principal voltage sensors in these proteins are positively charged S4 segments. The absolute magnitude of S4 movement discriminates two competing classes of gating models. In one class, the movement is <10 Angstrom due to the fact that the electric field is focused by aqueous crevices in the channel protein. In an alternative model, based in part on the crystal structure of a potassium channel, the side chains of S4 arginines move their charges across the bilayer's electric field, a distance of >25 Angstrom. Here, using tethered charges attached to an S4 segment, we provide evidence that the electric field falls across a distance of <4 Angstrom, supporting a model in which the relative movement between S4 and the electric field is very small.     

Cell Membrane Tethers Generate Mechanical Force in Response to Electrical Stimulation

Living cells maintain a huge transmembrane electric field across their membranes. This electric field exerts a force on the membrane because the membrane surfaces are highly charged. We have measured electromechanical force generation by cell membranes using optically trapped beads to detach the plasma membrane from the cytoskeleton and form long thin cylinders (tethers). Hyperpolarizing potentials increased and depolarizing potentials decreased the force required to pull a tether. The membrane tether force in response to sinusoidal voltage signals was a function of holding potential, tether diameter, and tether length. Membrane electromechanical force production can occur at speeds exceeding those of ATP-based protein motors. By harnessing the energy in the transmembrane electric field, cell membranes may contribute to processes as diverse as outer hair cell electromotility, ion channel gating, and transport.     

The influence of hydrophobic ions and dipolar molecules on the electrostatic barrier in biomembranes

To calculate the electric field inside a membrane the aqueous phase can be approximated by a conductor since the dielectric constant of water is much larger than that of the membrane. Then, using the method of image charges, ions adsorbed inside the membrane can be considered as dipoles and dipolar molecules adsorbed inside the membrane may similarly be regarded as sets of two similarly oriented dipoles. The microscopic interactions and, therefore, the spatial correlations of the adsorbed species can then be obtained. Together with the Gouy theory for the diffuse double layer these results allow the determination of the adsorbed phase—aqueous phase equilibrium. From the densities and spatial correlations of the adsorbed ions and dipolar species, their influence upon the electrostatic barrier as experienced by an ion translocating the membrane can be calculated. Changes observed in the relaxation time and initial conductance of translocating hydrophobic ions in voltage-pulse experiments on bilayer membranes are predicted using this model of the electrostatic barrier. In addition, an equation giving the surface tension as a function of the (non-ideal) adsorption of hydrophobic ions and dipoles is derived.