Matthew Hoyles

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 Ohms 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.

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 Poissons 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.

Solutions of Poisson's equation in channel-like geometries

Electric forces play a key role in the conductance of ions in biological channels. Therefore, their correct treatment is very important in making physical models of ion channels. Here, we present FORTRAN 90 codes for solution of Poissons equation satisfying the Dirichlet boundary conditions in realistic channel geometries that can be used in studies of ion channels. For a general channel shape, we discuss a numerical solution of Poissons equation based on an iterative technique. We also provide an analytical solution of Poissons equation in toroidal coordinates and its numerical implementation. A torus shaped channel is closer to reality than a cylindrical one, hence it could serve as a useful test model.

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.

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 Poissons 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 Poissons 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.

Molecular and Brownian dynamics study of ion selectivity and conductivity in the potassium channel

Abstract We employ recently revealed structural information for the potassium channel in molecular and Brownian dynamics simulations to investigate the physical mechanisms involved in the transport of ions across this channel. We show that ion selectivity arises from the ability of the channel protein to completely solvate potassium ions but not the smaller sodium ions. From energy and free energy perturbation profiles, we estimate the size of the energy barrier experienced by a sodium ion. Brownian dynamics simulations are carried out to determine conductance properties of this channel under various conditions.