Stochastic theory of ion movement in channels with single-ion occupancy. Application to sodium permeation of gramicidin channels

Solution of the Poisson-Nernst-Planck equations in the cell-substrate interface

Electrogenic cells are able to generate electrical signals which can be measured by various invasive electrophysiological methods such as patch-clamp or sharp microelectrode recordings. Growing cells on the surfaces of e.g. metal microelectrodes or field-effect transistors allows the recording of an extracellular component of these signals. For an understanding of such extracellular signals it is mandatory to get detailed topographical as well as electrical information about the cell-sensor interface. In a first approximation, this interface can be described by a flat disk between cell membrane and sensor surface. For a correct description of the signals, the electrodiffusion of ions in this interface is modeled by using the stationary Poisson-Nernst-Planck equations. We solve the equations analytically, and derive expressions for the potential, the ionic charge densities, and the seal resistance. The results provide a method for determining the distance h between sensor surface and cell membrane. For human embryonic kidney cells, we receive h approximately 70 nm. Comparison with literature shows good agreement.

Kernel representations for flux and concentration in ion channel models with time-varying concentrations

This paper explores stochastic models for the study of ion transport in biological cells. It considers one-dimensional models with time-varying concentrations at the boundaries. The average concentration and flux in the channel are obtained as kernel representations, where the kernel functions have a probabilistic interpretation which contributes to a better understanding of the models. In particular, the kernel representation is given for the flux at a boundary point, providing a correct version of a representation found in the literature. This requires special attention because one of the kernel functions exhibits a singularity. This kernel representation is feasible due to the linearity of the system that arises from the assumed independence between ions.

Steered Molecular Dynamics Simulations of Na+ Permeation across the Gramicidin A Channel

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.

Equivalence of trans paths in ion channels

We explore stochastic models for the study of ion transport in biological cells. Analysis of these models explains and explores an interesting feature of ion transport observed by biophysicists. Namely, the average time it takes ions to cross certain ion channels is the same in either direction, even if there is an electric potential difference across the channels. It is shown for simple single ion models that the distribution of a path (i.e., the history of location versus time) of an ion crossing the channel in one direction has the same distribution as the time-reversed path of an ion crossing the channel in the reverse direction. Therefore, not only is the mean duration of these paths equal, but other measures, such as the variance of passage time or the mean time a path spends within a specified section of the channel, are also the same for both directions of traversal. The feature is also explored for channels with interacting ions. If a system of interacting ions is in reversible equilibrium (net flux is zero), then the equivalence of the left-to-right trans paths with the time-reversed right-to-left trans paths still holds. However, if the system is in equilibrium, but not reversible equilibrium, then such equivalence need not hold.

Three computational methods for studying permeation, selectivity and dynamics in biological ion channels

The cell membrane, confining some ions and molecules on one side and exchanging others with the other side, is the ultimate unit of the physiology of life. The delicate task of regulating the transport of ions across the membrane is carried out by biological nanotubes called 'ion channels'. Recently, there have been enormous strides in our understanding of the structure-function relationships of biological ion channels. The molecular structures of several ion channels have been determined from crystallographic analysis, including potassium channels, mechanosensitive channels, a chloride channel, as well as gramicidin channels and porins. It is expected that the X-ray structures of other ion channels will soon follow these discoveries, ushering in a new era of ion channel studies in which predicting the function of channels from their atomic structures will become the main quest. In parallel to these experimental findings, there have been important advances in computational biophysics. Here we summarize three theoretical approaches that have been utilized to understand the dynamics of ion permeation across bio-nanotubes, highlighting their advantages and shortcomings, and briefly describe some of the salient properties of ion channels uncovered through computational studies.

Brownian Dynamics Simulation for Modeling Ion Permeation Across Bionanotubes

The principles underlying Brownian dynamics (BD), its statistical consistency, and algorithms for practical implementation are outlined here. The ability to compute current flow across ion channels confers a distinct advantage to BD simulations compared to other simulation techniques. Thus, two obvious applications of BD ion channels are in calculation of the current-voltage and current-concentration curves, which can be directly compared to the physiological measurements to assess the reliability of the model and predictive power of the method. We illustrate how BD simulations are used to unravel the permeation dynamics in two biological ion channels-the KcsA K+ channel and CIC Cl- channel.

Ionic permeation free energy in gramicidin: a semimicroscopic perspective

Ion permeation through the gramicidin channel is studied using a model that circumvents two major difficulties inherent to standard simulational methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization, accounting for the influence of electronic polarization by embedding the channel in a dielectric milieu that describes this polarization in a mean sense. The explicit mobile moieties are the ion, multipolar waters, and the carbonyls and amides of the peptide backbone. The model treats the influence of aromatic residues and the membrane dipole potential. A new electrical geometry is introduced that treats long-range electrostatics exactly and avoids problems related to periodic boundary conditions. It permits the translocating ion to make a seamless transition from nearby electrolyte to the channel interior. Other degrees of freedom (more distant bulk electrolyte and nonpolar lipid) are treated as dielectric continua. Reasonable permeation free energy profiles are obtained for potassium, rubidium, and cesium; binding wells are shallow and the central barrier is small. Estimated cationic single-channel conductances are smaller than experiment, but only by factors between 2 (rubidium) and 50 (potassium). When applied to chloride the internal barrier is large, with a corresponding miniscule single-channel conductance. The estimated relative single-channel conductances of gramicidin A, B, and C agree well with experiment.

Saturation of conductance in single ion channels: the blocking effect of the near reaction field

The ionic current flowing through a protein channel in the membrane of a biological cell depends on the concentration of the permeant ion, as well as on many other variables. As the concentration increases, the rate of arrival of bath ions to the channel's entrance increases, and typically so does the net current. This concentration dependence is part of traditional diffusion and rate models that predict Michaelis-Menten current-concentration relations for a single ion channel. Such models, however, neglect other effects of bath concentrations on the net current. The net current depends not only on the entrance rate of ions into the channel, but also on forces acting on ions inside the channel. These forces, in turn, depend not only on the applied potential and charge distribution of the channel, but also on the long-range Coulombic interactions with the surrounding bath ions. In this paper, we study the effects of bath concentrations on the average force on an ion in a single ion channel. We show that the force of the reaction field on a discrete ion inside a channel embedded in an uncharged lipid membrane contains a blocking (shielding) term that is proportional to the square root of the ionic bath concentration. We then show that different blocking strengths yield different behavior of the current-concentration and conductance-concentration curves. Our theory shows that at low concentrations, when the blocking force is weak, conductance grows linearly with concentration, as in traditional models, e.g., Michaelis-Menten formulations. As the concentration increases to a range of moderate shielding, conductance grows as the square root of concentration, whereas at high concentrations, with high shielding, conductance may actually decrease with increasing concentrations: the conductance-concentration curve can invert. Therefore, electrostatic interactions between bath ions and the single ion inside the channel can explain the different regimes of conductance-concentration relations observed in experiments.

Modeling of ionic relaxation around a biomembrane disk

A simplified Brownian dynamics model and the corresponding software implementation have been developed for the simulation of electrolyte dynamics on the mesoscopic scale. In addition to direct control simulations, the model system has been verified by a quantitative comparison with the Debye-Hückel theory. As a first application, the model was used to simulate ionic relaxation processes following abrupt intramembrane charge rearrangements in the case of a disk shaped membrane. In addition to its general implications, the obtained properties of the relaxation kinetics confirm the assumptions of the theory of the so-called suspension method, a technique capable of tracing molecular charge motions of membrane proteins in three dimensions.

Voltage-gated proton channels and other proton transfer pathways

Proton channels exist in a wide variety of membrane proteins where they transport protons rapidly and efficiently. Usually the proton pathway is formed mainly by water molecules present in the protein, but its function is regulated by titratable groups on critical amino acid residues in the pathway. All proton channels conduct protons by a hydrogen-bonded chain mechanism in which the proton hops from one water or titratable group to the next. Voltage-gated proton channels represent a specific subset of proton channels that have voltage- and time-dependent gating like other ion channels. However, they differ from most ion channels in their extraordinarily high selectivity, tiny conductance, strong temperature and deuterium isotope effects on conductance and gating kinetics, and insensitivity to block by steric occlusion. Gating of H(+) channels is regulated tightly by pH and voltage, ensuring that they open only when the electrochemical gradient is outward. Thus they function to extrude acid from cells. H(+) channels are expressed in many cells. During the respiratory burst in phagocytes, H(+) current compensates for electron extrusion by NADPH oxidase. Most evidence indicates that the H(+) channel is not part of the NADPH oxidase complex, but rather is a distinct and as yet unidentified molecule.

Recent advances in ion channel research

The field of ion channels has entered into a rapid phase of development in the last few years, partly due to the breakthroughs in determination of the crystal structures of membrane proteins and advances in computer simulations of biomolecules. These advances have finally enabled the long-dreamed goal of relating function of a channel to its underlying molecular structure. Here we present simplified accounts of the competing permeation theories and then discuss their application to the potassium, gramicidin A and calcium channels.

Continuum electrostatics fails to describe ion permeation in the gramicidin channel

We investigate the validity of continuum electrostatics in the gramicidin A channel using a recently determined high-resolution structure. The potential and electric field acting on ions in and around the channel are computed by solving Poisson's equation. These are then used in Brownian dynamics simulations to obtain concentration profiles and the current passing through the channel. We show that regardless of the effective dielectric constant used for water in the channel or the channel protein, it is not possible to reproduce all the experimental data on gramicidin A; thus, continuum electrostatics cannot provide a valid framework for the description of ion dynamics in gramicidin channels. Using experimental data and molecular dynamics simulations as guides, we have constructed potential energy profiles that can satisfactorily describe the available physiological data. These profiles provide useful benchmarks for future potential of mean force calculations of permeating ions from molecular dynamics simulations of gramicidin A. They also offer a convenient starting point for studying structure-function relationships in modified gramicidin channels.

Computational studies of the gramicidin channel

Ion channels are highly specific membrane-spanning protein structures which serve to facilitate the passage of selected ions across the lipid barrier. In the past decade, molecular dynamics simulations based on atomic models and realistic microscopic interactions with explicit solvent and membrane lipids have been used to gain insight into the function of these complex systems. These calculations have considerably expanded our view of ion permeation at the microscopic level. This Account will mainly focus on computational studies of the gramicidin A channel, one of the simplest and best characterized molecular pore.

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.

Reservoir boundaries in Brownian dynamics simulations of ion channels

Brownian dynamics (BD) simulations provide a practical method for the calculation of ion channel conductance from a given structure. There has been much debate about the implementation of reservoir boundaries in BD simulations in recent years, with claims that the use of improper boundaries could have large effects on the calculated conductance values. Here we compare the simple stochastic boundary that we have been using in our BD simulations with the recently proposed grand canonical Monte Carlo method. We also compare different methods of creating transmembrane potentials. Our results confirm that the treatment of the reservoir boundaries is mostly irrelevant to the conductance properties of an ion channel as long as the reservoirs are large enough.

Cellular microtransport processes: intercellular, intracellular, and aggregate behavior

Ionic and molecular transfer among cells occurs by a variety of transport processes operative at different length scales. Cell membrane permeability and electrical conductance derive from channel proteins producing pores at the molecular (ultrastructural) scale. Intracellular mobility involves the dynamics of motion through the complex ultrastructure of the cytoplasm. These phenomena unite in the larger-scale (microscopic) process of gross intercellular transfer. When such movement occurs among sufficiently many cells, it in turn begins to reflect their average collective (macroscopic) behavior as bulk tissue. This article surveys selected aspects of intercellular and intracellular transport, with emphasis on detailed mechanistic theory, experimental probes of cellular permeability, and systematic transcendence from small to large length scales.

Model channel ion currents in NaCl-extended simple point charge water solution with applied-field molecular dynamics

Using periodic boundary conditions and a constant applied field, we have simulated current flow through an 8.125-A internal diameter, rigid, atomistic channel with polar walls in a rigid membrane using explicit ions and extended simple point charge water. Channel and bath currents were computed from 10 10-ns trajectories for each of 10 different conditions of concentration and applied voltage. An electric field was applied uniformly throughout the system to all mobile atoms. On average, the resultant net electric field falls primarily across the membrane channel, as expected for two conductive baths separated by a membrane capacitance. The channel is rarely occupied by more than one ion. Current-voltage relations are concentration dependent and superlinear at high concentrations.

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.

Noncontact dipole effects on channel permeation. V. Computed potentials for fluorinated gramicidin

Experimental and theoretical calculations indicate that the dipole moment of the four Trp side chains in gramicidin A (gA) channels modify channel conductance through long-range electrostatic interactions. Electrostatic ion/side-chain interaction energies along the channel were computed with CHARMM using ab initio atom charges for native and 4-, 5-, or 6-fluorinated Trp side chains. The bulk water reaction to the polar side chains was included using the method of images as implemented by, and channel waters in idealized structures were included. Ion/Trp interaction energies were approximately -0.6 kcal/mol throughout the channel for all four of the native Trp pairs. Channel waters produced a modest reduction in the magnitude of interactions, essentially offsetting images representing the bulk water outside the channel. The effects of side-chain fluorination depended on ring position and, to a lesser extent, residue number. Compared with native Trp, 5-fluorination reduces the translocation barrier with minor effects on the exit barrier. In contrast, 6-fluorination primarily reduces exit barrier. 4-Fluorination produces a more complex double-well energy profile. Effects of measured side-chain movements resulting from fluorination or change in lipid bilayer were negligible whereas thermal side chain librations cause large effects, especially in the region of the ion-binding sites.

Framework model for single proton conduction through gramicidin

This paper describes a framework model for proton conduction through gramicidin; a model designed to incorporate information from molecular dynamics and use this to predict conductance properties. The state diagram describes both motion of an excess proton within the pore as well as the reorientation of waters within the pore in the absence of an excess proton. The model is constructed as the diffusion limit of a random walk, allowing control over the boundary behavior of trajectories. Simple assumptions about the boundary behavior are made, which allow an analytical solution for the proton current and conductance. This is compared with corresponding expressions from statistical mechanics. The random walk construction allows diffusing trajectories underlying the model to be simulated in a simple way. Details of the numerical algorithm are described.

A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels

A computational algorithm based on Grand Canonical Monte Carlo (GCMC) and Brownian Dynamics (BD) is described to simulate the movement of ions in membrane channels. The proposed algorithm, GCMC/BD, allows the simulation of ion channels with a realistic implementation of boundary conditions of concentration and transmembrane potential. The method is consistent with a statistical mechanical formulation of the equilibrium properties of ion channels (; Biophys. J. 77:139-153). The GCMC/BD algorithm is illustrated with simulations of simple test systems and of the OmpF porin of Escherichia coli. The approach provides a framework for simulating ion permeation in the context of detailed microscopic models.

Tests of continuum theories as models of ion channels. II. Poisson-Nernst-Planck theory versus brownian dynamics

We test the validity of the mean-field approximation in Poisson-Nernst-Planck theory by contrasting its predictions with those of Brownian dynamics simulations in schematic cylindrical channels and in a realistic potassium channel. Equivalence of the two theories in bulk situations is demonstrated in a control study. In simple cylindrical channels, considerable differences are found between the two theories with regard to the concentration profiles in the channel and its conductance properties. These differences are at a maximum in narrow channels with a radius smaller than the Debye length and diminish with increasing radius. Convergence occurs when the channel radius is over 2 Debye lengths. These tests unequivocally demonstrate that the mean-field approximation in the Poisson-Nernst-Planck theory breaks down in narrow ion channels that have radii smaller than the Debye length.

Molecular dynamics of synthetic leucine-serine ion channels in a phospholipid membrane

Molecular dynamics calculations were carried out on models of two synthetic leucine-serine ion channels: a tetrameric bundle with sequence (LSLLLSL)(3)NH(2) and a hexameric bundle with sequence (LSSLLSL)(3)NH(2). Each protein bundle is inserted in a palmitoyloleoylphosphatidylcholine bilayer membrane and solvated by simple point charge water molecules inside the pore and at both mouths. Both systems appear to be stable in the absence of an electric field during the 4 ns of molecular dynamics simulation. The water motion in the narrow pore of the four-helix bundle is highly restricted and may provide suitable conditions for proton transfer via a water wire mechanism. In the wider hexameric pore, the water diffuses much more slowly than in bulk but is still mobile. This, along with the dimensions of the pore, supports the observation that this peptide is selective for monovalent cations. Reasonable agreement of predicted conductances with experimentally determined values lends support to the validity of the simulations.

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.

Effects of nonequilibrium fluctuations on ionic transport through biomembranes

We investigate the effects of nonequilibrium fluctuations on ionic transport through ion channels in membranes using the concept of localized ratchet. Due to the localization, the ionic population in the binding site can be enhanced or suppressed depending upon ionic potential and its fluctuations, affecting the gating kinetics of the channel. The localized dichotomic fluctuations of ionic potential are shown to give rise to a current reversal differing from the results of periodic ratchets. It is also found that strong correlations between binding energy and membrane potential fluctuations induce resonancelike behaviors in ionic current as the fluctuating rate varies.

Effective diffusion coefficients of K+ and Cl- ions in ion channel models

We have used molecular dynamics simulations, corresponding to a total simulation time of 11 ns, to investigate the effective short-time local diffusion coefficient of potassium and chloride ions in a series of model ion channels. These models, which include channels formed by the fungal peptide alamethicin, by a synthetic leucine-serine peptide, and by the pore-lining M2 helix bundle of the nicotinic acetylcholine receptor, have a range of different secondary structures, diameters and hydrophobicities. We find that the diffusion coefficients of both ions are appreciably reduced in the narrower channels, the extent of the reduction being similar for both the anionic and cationic species. This suggests that a difference in mobility cannot be the source of the ion selectivity exhibited by some of the channels (for example, the leucine-serine peptide). We find no evidence for a reduction in mobility of either ion in the nAChR model. These results are broadly in line with a previous similar study of Na+ ions, and may be useful in Poisson-Nernst-Planck, Eyring rate theory or Brownian dynamics calculations of channel conductance.

Dynamic properties of Na+ ions in models of ion channels: a molecular dynamics study

We present simulation results for the effective diffusion coefficients of a sodium ion in a series of model ion channels of different diameters and hydrophobicities, including models of alamethicin, a leucine-serine peptide, and the M2 helix bundle of the nicotinic acetylcholine receptor. The diffusion coefficient, which in the simulations has a value of 0.15(2) A2ps-1 in bulk water, is found to be reduced to as little as 0.02(1) A2ps-1 in the narrower channels, and to about 0.10(5) A2ps-1 in wider channels such as the nicotinic acetylcholine receptor. It is anticipated that this work will be useful in connection with calculations of channel conductivity using such techniques as the Poisson-Nernst-Planck equation, Eyring rate theory, or Brownian dynamics.

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.

Using theory and simulation to understand permeation and selectivity in ion channels

It is clear that the function of ion channels must flow from their structure. With recent advances in computational power and methodology, it appears feasible to correlate structure to ion channel permeation at an atomistically detailed level of description. The overall strategy is to structure the calculations in a hierarchy, ranging from coarse-grained thermodynamic and kinetic descriptions to fine-grained molecular dynamics descriptions with atomic detail. Each level of description is connected to the others by appropriate statistical mechanical theory. The coarse-grained descriptions can be correlated directly with electrophysiological experiment. The fine-grained descriptions are used to parameterize the coarse-grained descriptions and to describe the permeation process at the most detailed level. This strategy has so far had varying degrees of success. It has successfully described water permeation through lipid bilayers and gramicidin channels. It has revealed the essential events of ion permeation through gramicidin channels at an atomistically detailed level. The role of channel protein motions in permeation has been elucidated. However, it appears that force fields used to describe molecular dynamics must be refined further to achieve completely accurate predictions of the permeation of such small ions as sodium. Channels with more complex structure and more multiion occupancy than gramicidin pose major computational challenges with respect to sampling protein conformations and ion distributions involved in the permeation process. Possible approaches to meeting these challenges are discussed.

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.

Boundary conditions for- single-ion diffusion

We have constructed a theory for diffusion through the pore of a single-ion channel by taking a limit of a random walk around a cycle of states. Similar to Levitt's theory of single-ion diffusion, one obtains boundary conditions for the Nernst-Planck equation that guarantee that the pore is occupied by at most one ion. Two of the terms in the boundary conditions are identical to those given by Levitt. However, the construction gives rise to a third term not found in Levitt's theory. With this term, the channel spends exponentially distributed intervals in the empty state. Ion sample paths have been simulated to help visualize trajectories near the channel entrances, with and without the new term. We use the modified Levitt theory to fit several potential profiles to the conductance data of Russell et al. In particular, we have analyzed the profile for Na+ in gramicidin calculated by Roux and Karplus. The peak-to-peak amplitude of their result must be reduced to at most 35% of its original value to fit the data. But with this reduction, excellent fits are obtained.

A semi-microscopic Monte Carlo study of permeation energetics in a gramicidin-like channel: the origin of cation selectivity

The influence of a gramicidin-like channel former on ion free energy barriers is studied using Monte Carlo simulation. The model explicitly describes the ion, the water dipoles, and the peptide carbonyls; the remaining degrees of freedom, bulk electrolyte, non-polar lipid and peptide regions, and electronic (high frequency) permittivity, are treated in continuum terms. Contributions of the channel waters and peptide COs are studied both separately and collectively. We found that if constrained to their original orientations, the COs substantially increase the cationic permeation free energy; with or without water present, CO reorientation is crucial for ion-CO interaction to lower cation free energy barriers; the translocation free energy profiles for potassium-, rubidium-, and cesium-like cations exhibit no broad barriers; the lipid-bound peptide interacts more effectively with anions than cations; anionic translocation free energy profiles exhibit well defined maxima. Using experimental data to estimate transfer free energies of ions and water from bulk electrolyte to a non-polar dielectric (continuum lipid), we found reasonable ion permeation profiles; cations bind and permeate, whereas anions cannot enter the channel. Cation selectivity arises because, for ions of the same size and charge, anions bind hydration water more strongly.

Extended dipolar chain model for ion channels: electrostriction effects and the translocational energy barrier

We reinvestigate the dipolar chain model for an ion channel. Our goal is to account for the influence that ion-induced electrostriction of channel water has on the translocational energy barriers experienced by different ions in the channel. For this purpose, we refine our former model by relaxing the positional constraint on the ion and the water dipoles and by including Lennard-Jones contributions in addition to the electrostatic interactions. The positions of the ion and the waters are established by minimization of the free energy. As before, interaction with the external medium is described via the image forces. Application to alkali cations show that the short range interactions modulate the free energy profiles leading to a selectivity sequence for translocation. We study the influence of some structural parameters on this sequence and compare our theoretical predictions with observed results for gramicidin.

Influence of a channel-forming peptide on energy barriers to ion permeation, viewed from a continuum dielectric perspective

The continuum three-dielectric model for an aqueous ion channel pore-forming peptide-membrane system is extended to account for the finite length of the channel. We focus on the electrostatic influence that a channel-forming peptide may exert on energy barriers to ion permeation. The nonlinear dielectric behavior of channel water caused by dielectric saturation in the presence of an ion is explicitly modeled by assigning channel water a mean dielectric constant much less than that of bulk water. An exact solution of the continuum problem is formulated by approximating the dielectric behavior of bulk water, assigning it a dielectric constant of infinity. The validity of this approximation is verified by comparison with a Poisson-Boltzmann description of the electrolyte. The formal equivalence of high ionic strength and high electrolyte dielectric constant is demonstrated. We estimate limits on the reduction of the electrostatic free energy caused by ionic interaction with the channel-forming peptide. We find that even assigning this region an epsilon of 100, its influence is insufficient to lower permeation free energy barriers to values consistent with observed channel conductances. We provide estimates of the effective dielectric constant of this highly polarizable region, by comparing energy barriers computed using the continuum approach with those found from a semi-microscopic analysis of a simplified model of a gramicidin-like charge distribution. Possible ways of improving both models are discussed.

Brownian dynamics study of a multiply-occupied cation channel: application to understanding permeation in potassium channels

The behavior of a multiply-occupied cation-selective channel has been computed by Brownian dynamics. The length, cross-section, ion-ion repulsion force, and ionic mobility within the channel are all estimated from data and physical reasoning. The only free parameter is a partition energy at the mouth of the channel, defining the free energy of an ion in the channel compared to the bath. It is presumed that this partition energy is associated with the energetics of exchanging a bulk hydration environment for a channel hydration environment. Varying the partition energy alone, keeping all other parameters fixed, gives approximately the full range of magnitudes of single channel conductances seen experimentally for K channels. Setting the partition energy at -11 kT makes the computed channel look similar to a squid axon K channel with respect to magnitude of conductance, shape of the I-V curve, non-unity of Ussing flux ratio exponents, decrease of current and increase of conductance with extracellular ion accumulation, and saturation at high ion concentration in the bathing solution. The model includes no preferred binding sites (local free energy minima) for ions in the channel. Therefore it follows that none of the above-mentioned properties of K channels are strong evidence for the existence of such sites. The model does not show supersaturation of current at very high bathing concentrations nor any pronounced voltage-dependence of the Ussing flux ratio exponent, suggesting that these features would require additional details not included in the model presented herein.

Permeation in ionic channels: a statistical rate theory approach

A novel way to model permeation through ionic channels is formulated. Our method does not require that equilibrium exists in the channel or at the channel interfaces. In addition, the potential profile does not need to be specified and the assumption of constant field across the membrane does not need to be made. Our formulation relies on statistical rate theory for its development and uses a form of the electrochemical potential which assumes that the ions are in solution. We show that the conductance and the degree of nonlinearity are dependent on the relative equilibrium exchange rates in the channel and at the interfaces. Nonlinear current-voltage plots can be obtained in symmetric solutions as well as a nonunity exponent for the Ussing flux ratio. Due to the dependence of the partition coefficient on solubility, it is highly unlikely that the intracellular and extracellular partition coefficients are the same. A manifestation of unequal partition coefficients is a current reversal at a membrane voltage that is different from the Nernst potential of the current-carrying ionic species.

The permeation properties of small organic cations in gramicidin A channels

The conductance properties of organic cations in single gramicidin A channels were studied using planar lipid bilayers. From measurements at 10 mM and at 27 mV the overall selectivity sequence was found to be NH4+ > K+ > hydrazinium > formamidinium > Na+ > methylammonium, which corresponds to Eisenman polyatomic cation sequence X'. Methylammonium and formamidinium exhibit self block, suggesting multiple occupancy and single filing. Formamidinium has an apparent dissociation constant (which is similar to those of alkali metal cations) for the first ion being 22 mM from the Eadie-Hofstee plot (G0 vs. G0/C), 12 mM from the rate constants of a three-step kinetic model. The rate-limiting step for formamidinium is translocation judging from supralinear I-V relations at low concentrations. 1 M formamidinium solutions yields exceptionally long single channel lifetimes, 20-fold longer than methylammonium, which yields lifetimes similar to those found with alkali metal cations. The average lifetime in formamidinium solution significantly decreases with increasing voltage up to 100 mV but is relatively voltage independent between 100 and 200 mV. At lower voltages (< or = 100 mV), the temperature and concentration dependences of the average lifetime of formamidinium were steep. At very low salt concentrations (0.01 M, 100 mV), there was no significant difference in average lifetime from that formed with 0.01 M methylammonium or hydrazinium. We conclude that formamidinium very effectively stabilizes the dimeric channel while inside the channel and speculate that it does so by affecting tryptophan-reorientation or tryptophan-lipid interactions at binding sites.

The nature of ion and water barrier crossings in a simulated ion channel

Using a combination of techniques, including molecular dynamics, time-correlation analysis, stochastic dynamics, and fitting of continuum diffusion theory to electrophysiological data, a characterization is made of thermally driven sodium, water, and D2O motion within the gramicidin A channel. Since the channel contents are constrained to move in a single-file fashion, the motion that corresponds to experimentally measurable rates of permeation of the membrane is the motion of the center of mass of the channel contents. We therefore emphasize channel contents center-of-mass motion in our analysis of molecular dynamics computations. The usual free energy calculation techniques would be of questionable validity when applied to such motion. As an alternative to those techniques, we postulate a periodic sinusoidal free energy profile (related to the periodic structure of the helical channel) and deduce the fluid dynamic diffusion coefficient and the height and spacing of the free energy barriers from the form of the mean-square-deviation function, using stochastic computations. The fluid dynamic friction in each case appears similar to that for aqueous solution. However, the diffusive motions are modulated by a spatially periodic free energy profile with a periodicity characteristic of an L-D pair of amino acids in the gramicidin helix, approximately 1.7 A in the model we use. The barrier height depends on which substance is moving in the channel, but in each case is several times thermal energy. For barriers of this width and height, the motion is intermediate between the low-friction (transition-state) and high-friction (Brownian) limits. Thus, neither of these formalisms that have been used commonly to describe membrane permeation gives an accurate picture of the underlying physical process (although the Brownian description seems closer to correct). The non-Markovian Langevin equation must be solved to describe properly the statistics of the process. The "channel state of matter" characteristic of the channel contents appears to have some properties typical of the solid and some typical of the liquid state. The magnitude of the local friction and nature of the ion solvation are similar to the liquid state, but the periodicities of structure, free energy, and dynamics are somewhat solid-like. The alignment of water dipoles in the channel bears some resemblance to the orientational ordering of a nematic liquid crystal, but unlike a nematic liquid crystal, the waters have a degree of translational order as well. Thus, the "channel state" is not adequately described by analogy to either the solid or liquid states or to liquid crystals but must be dealt with as its own characteristic type of condensed matter.

Conformational model for ion permeation in membrane channels: a comparison with multi-ion models and applications to calcium channel permeability

The permeation properties of ion channels existing in several conductive states were analyzed. Each state was represented by the one-ion model. A special emphasis was placed on features, assumed to be indicative of a multi-ion mode of channel occupancy such as a deviation of concentration dependence of channel conductance from the Michaelis-Menten equation, an anomalous mole fraction effect, a strong voltage dependence of ion block and coupling of unidirectional fluxes (anomalous Ussing flux ratio). The conformational model was shown to have all these properties. The ion permeation through voltage-sensitive calcium channels fulfilled all the characteristics of the model proposed.