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.     

The pore helix dipole has a minor role in inward rectifier channel function

Ion channels lower the energetic barrier for ion passage across cell membranes and enable the generation of bioelectricity. Electrostatic interactions between permeant ions and channel pore helix dipoles have been proposed as a general mechanism for facilitating ion passage. Here, using genetic selections to probe interactions of an exemplar potassium channel blocker, barium, with the inward rectifier Kir2.1, we identify mutants bearing positively charged residues in the potassium channel signature sequence at the pore helix C terminus. We show that these channels are functional, selective, resistant to barium block, and have minimally altered conductance properties. Both the experimental data and model calculations indicate that barium resistance originates from electrostatics. We demonstrate that potassium channel function is remarkably unperturbed when positive charges occur near the permeant ions at a location that should counteract pore helix electrostatic effects. Thus, contrary to accepted models, the pore helix dipole seems to be a minor factor in potassium channel permeation.     

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.     

Brownian dynamics study of ion transport in the vestibule of membrane channels.

Brownian dynamics simulations have been carried out to study the transport of ions in a vestibular geometry, which offers a more realistic shape for membrane channels than cylindrical tubes. Specifically, we consider a torus-shaped channel, for which the analytical solution of Poisson's equation is possible. The system is composed of the toroidal channel, with length and radius of the constricted region of 80 A and 4 A, respectively, and two reservoirs containing 50 sodium ions and 50 chloride ions. The positions of each of these ions executing Brownian motion under the influence of a stochastic force and a systematic electric force are determined at discrete time steps of 50 fs for up to 2.5 ns. All of the systematic forces acting on an ion due to the other ions, an external electric field, fixed charges in the channel protein, and the image charges induced at the water-protein boundary are explicitly included in the calculations. We find that the repulsive dielectric force arising from the induced surface charges plays a dominant role in channel dynamics. It expels an ion from the vestibule when it is deliberately put in it. Even in the presence of an applied electric potential of 100 mV, an ion cannot overcome this repulsive force and permeate the channel. Only when dipoles of a favorable orientation are placed along the sides of the transmembrane segment can an ion traverse the channel under the influence of a membrane potential. When the strength of the dipoles is further increased, an ion becomes detained in a potential well, and the driving force provided by the applied field is not sufficient to drive the ion out of the well. The trajectory of an ion navigating across the channel mostly remains close to the central axis of the pore lumen. Finally, we discuss the implications of these findings for the transport of ions across the membrane.     

Membrane dipole potential modulates proton conductance through gramicidin channel: movement of negative ionic defects inside the channel.

The effect of membrane dipole potential on gramicidin channel activity in bilayer lipid membranes (BLMs) was studied. Remarkably, it appeared that proton conductance of gramicidin A (gA) channels responded to modulation of the dipole potential oppositely as compared with gA alkali metal cation conductance. In particular, the addition of phloretin, known to reduce the membrane dipole potential, resulted in a decrease in gA proton conductance, on one hand, and an increase in gA alkali metal conductance, on the other hand, whereas 6-ketocholestanol, the agent raising the membrane dipole potential, provoked an increase in gA proton conductance as opposed to a decrease in the alkali metal cation conductance. The peculiarity of the 6-ketocholestanol effect consisted in its dependence on the H(+) concentration. The experiments with the impermeant dipolar compound, phloridzin, showed that the response of proton transport through gramicidin channels to varying the membrane dipole potential did not change qualitatively if the dipole potential of only one monolayer or both monolayers of the BLM was altered. In contrast to gA proton conductance, the single-channel lifetime changed similarly with varying the membrane dipole potential, regardless of the kind of permeant cations (protons or potassium ions). The results of this study could be tentatively accounted for by an assumption that one of the rate-limiting steps of proton conduction through gramicidin channels represents, in fact, movement of negatively charged species (negative ionic defects) across a membrane.     

The predominant role of coordination number in potassium channel selectivity

Potassium channels are exquisitely selective, allowing K+ to pass across cell membranes while blocking other ion types. Here we demonstrate that the number of carbonyl oxygen atoms that surround permeating ions is the most important factor in determining ion selectivity rather than the size of the pore or the strength of the coordinating dipoles. Although the electrostatic properties of the coordinating ligands can lead to Na+ or K+ selectivity at some values of the dipole moment, no significant selectivity arises at the specific value of the dipole moment for carbonyl groups found in potassium channels when the ligands have complete freedom. Rather, we show that the main contribution to selectivity arises from slight constraints on the conformational freedom of the channel protein that limit the number of carbonyl oxygen atoms to a value better suited to K+ than Na+, despite the pore being flexible. This mechanism provides an example of a general framework for explaining ion discrimination in a range of natural and synthetic macromolecules in which selectivity is controlled by the number of coordinating ligands in addition to their dipole moment.     

Energy barriers for passage of ions through channels. Exact solution of two electrostatic problems

Exact solutions are given to two electrostatic problems relevant to ion permeation through pores in membranes. The first assesses the importance of the pore forming molecule as a dielectric shield. It is shown on the basis of structural and dielectric considerations alone (neglecting effects attributable to possible charge distribution at the interior surface of the pre-former) that the minimum electrostatic barrier for monovalent ion passage through a gramicidin-like channel is 11 kT. It is further shown that given favorable circumstances, dielectric shielding might dramatically reduce the barrier to ion passage through potassium channels. The second problem considers the error introduced by treating ions as point charges. It is shown that for structureless pores the point charge approximation introduces no meaningful error, even if the ratio of ion radius to pore radius is as great as 0.95.     

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.     

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.     

Langevin Poisson-Boltzmann equation: point-like ions and water dipoles near a charged surface

Water ordering near a charged membrane surface is important for many biological processes such as binding of ligands to a membrane or transport of ions across it. In this work, the mean-field Poisson-Boltzmann theory for point-like ions, describing an electrolyte solution in contact with a planar charged surface, is modified by including the orientational ordering of water. Water molecules are considered as Langevin dipoles, while the number density of water is assumed to be constant everywhere in the electrolyte solution. It is shown that the dielectric permittivity of an electrolyte close to a charged surface is decreased due to the increased orientational ordering of water dipoles. The dielectric permittivity close to the charged surface is additionally decreased due to the finite size of ions and dipoles.     

Batrachotoxin-modified sodium channels in planar lipid bilayers. Ion permeation and block

Batrachotoxin-modified, voltage-dependent sodium channels from canine forebrain were incorporated into planar lipid bilayers. Single-channel conductances were studied for [Na+] ranging between 0.02 and 3.5 M. Typically, the single-channel currents exhibited a simple two-state behavior, with transitions between closed and fully open states. Two other conductance states were observed: a subconductance state, usually seen at [NaCl] greater than or equal to 0.5 M, and a flickery state, usually seen at [NaCl] less than or equal to 0.5 M. The flickery state became more frequent as [NaCl] was decreased below 0.5 M. The K+/Na+ permeability ratio was approximately 0.16 in 0.5 and 2.5 M salt, independent of the Na+ mole fraction, which indicates that there are no interactions among permeant ions in the channels. Impermeant and permeant blocking ions (tetraethylammonium, Ca++, Zn++, and K+) have different effects when added to the extracellular and intracellular solutions, which indicates that the channel is asymmetrical and has at least two cation-binding sites. The conductance vs. [Na+] relation saturated at high concentrations, but could not be described by a Langmuir isotherm, as the conductance at low [NaCl] is higher than predicted from the data at [NaCl] greater than or equal to 1.0 M. At low [NaCl] (less than or equal to 0.1 M), increasing the ionic strength by additions of impermeant monovalent and divalent cations reduced the conductance, as if the magnitude of negative electrostatic potentials at the channel entrances were reduced. The conductances were comparable for channels in bilayers that carry a net negative charge and bilayers that carry no net charge. Together, these results lead to the conclusion that negative charges on the channel protein near the channel entrances increase the conductance, while lipid surface charges are less important.     

Two Distinct K+ Channels in Lamprey (Lampetra fluviatilis) Erythrocyte Membrane Characterized by Single Channel Patch Clamp

Two channels, distinguished by using single-channel patch-clamp, carry out potassium transport across the red cell membrane of lamprey erythrocytes. A small-conductance, inwardly rectifying K⁺-selective channel was observed in both isotonic and hypotonic solutions (osmolarity decreased by 50%). The single-channel conductance was 26 ± 3 pS in isotonic (132 mm K⁺) solutions and 24 ± 2 pS in hypotonic (63 mm K⁺) solutions. No outward conductance was found for this channel, and the channel activity was completely inhibited by barium. Cell swelling activated another inwardly rectifying K⁺ channel with a larger inward conductance of 65 pS and outward conductance of 15 pS in the on-cell configuration. In this channel, rectification was due to the block of outward currents by Mg²⁺ and Ca²⁺ ions, since when both ions were removed from the cytosolic side in inside-out patches the conductance of the channel was nearly ohmic. In contrast to the small-conductance channel, the swelling-activated channel was observed also in the presence of barium in the pipette. Neither type of channel was dependent on the presence of Ca²⁺ ions on the cytosolic side for activity. Received: 18 July 1997/Revised: 30 January 1998     

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.     

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.     

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.     

The role of proteins in a dipole model for steady-state ionic transport through biological membranes.

The steady-state current-voltage characteristics of biological membranes are analyzed for means of an application of the electrodiffusion theory to the passage of ions through "dielectric pores", with orientable dipoles at the pore-water interfaces. A detailed evaluation of the electrostatic potential barrier shows, indeed, that the ions have practically no chance to penetrate into the phospholipid bilayer, but that they can cross the membrane through local protein inclusions, of high dielectric constant. A "gating mechanism" can be provided, moreover, by a change of the potential barrier, resulting from a dipole reorientation at the pore-water interface. Dipole-dipole interactions are opposed to the orienting effect of an applied field, but they can be neglected when the separation between the dipoles exceeds a certain critical value. The high polarizability of the pore material leads to an amplification of the effect of an applied field on the orientable dipoles. It is therefore possible to achieve a satisfactory agreement with the experimental results of Gilbert and Ehrenstein (Biophys. J., 9: 447, 1969) for the squid axon, and, in particular, to account for the width of the negative resistance regions with a relatively small value for the length of the orientable dipoles.     

How does vestibule surface charge affect ion conduction and toxin binding in a sodium channel?

We describe various models for the dielectric geometry and pore mouth charge distribution of a Na channel. The electric potential due to the vestibule charges is then computed on the basis of the nonlinear Possion-Boltzmann equation. The results are used to account for the effect of permeant ion concentration and ionic strength on channel conductance and on toxin association rate constants for Na channels. We find that a single negatively charged group near the entrance to the channel constriction is adequate to account for deviations from Michaelis-Menten conductance kinetics and for the concentration dependence of toxin-binding coefficients. We find further that only a limited range of vestibule geometries and pore mouth charge distributions are consistent with experiment.     

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.     

The different screening of electric charges and dipoles near a dielectric interface

The electric fields within a planar slab of material due to both charges and dipoles within the slab and near to its surface are simply calculated using the method of images. If the slab is immersed in a fluid of high dielectric constant the electric field within the slab due to the charge is always reduced but that of a suitably oriented electric dipole is enhanced by as much as a factor of two.     

A novel mechanism for voltage control of channel conductance

Many channel-formers can exist in conformational states with varying degrees of conductance. When the difference in the energy level between states is voltage dependent the result is a voltage-dependent channel. This voltage-dependence is usually attributable to the movement of charges or alignment of dipoles associate with channel-former. This paper presents a simpler mechanism by which the energy difference can be voltage dependent without the need for charge movement or dipole alignment. The voltage dependent energy difference between the high and low conducting states can arise from a change in electrical potential at a fixed charge located on the walls lining the pore (or field at a dipole fixed within the pore) which occurs as a result of the conformational change. This mechanism takes advantage of structural features normally found in channel formers and local changes resulting from channel opening and closing to generate the energy difference needed for voltage-dependence.