Search PubMed⌕ Search

Biomedical subjects

S Kuyucak

Publications and source records attributed to S Kuyucak.

At least 19 recordsLinked to original sources
PubMed2001

Mechanisms of permeation and selectivity in calcium channels.

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

Biophysical Phenomena↗
PubMed2000

A model of calcium channels.

We propose a model of calcium channels that can explain most of their observed properties, including the anomalous mole fraction effect and mutation of the glutamate residues. The structure grossly resembles that of the KcsA potassium channel except for the presence of an extracellular vestibule and a shorter selectivity filter containing four glutamate residues. Using this model in electrostatic calculations and Brownian dynamics simulations, we study mechanisms of ion permeation and selectivity in the channel. Potential energy profiles calculated for multiple ions in the channel provide explanations of ion permeation, the block of Na(+) currents by Ca(2+) ions, and many other observed properties. Brownian dynamics simulations provide quantitative predictions for the channel currents which reproduce available experimental data.

Calcium Channels↗
PubMed2000

Molecular dynamics estimates of ion diffusion in model hydrophobic and KcsA potassium channels.

Molecular dynamics simulations are carried out to obtain estimates of diffusion coefficients of biologically important Na+, K+, Ca2+ and Cl- ions in hydrophobic cylindrical channels with varying radii and large reservoirs. Calculations for the cylindrical channels are compared to those for the KcsA potassium channel, for which the protein structure has recently been determined from X-ray diffraction experiments. Our results show that ion diffusion is maintained at reasonably high levels even within narrow channels, and does not support the very small diffusion coefficients used in some continuum models in order to fit experimental data. The present estimates of ion diffusion coefficients are useful in the calculation of channel conductance using the Poisson-Nernst-Planck theory, or Brownian dynamics.

Algorithms↗
PubMed2000

Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics.

Continuum theories of electrolytes are widely used to describe physical processes in various biological systems. Although these are well-established theories in macroscopic situations, it is not clear from the outset that they should work in small systems whose dimensions are comparable to or smaller than the Debye length. Here, we test the validity of the mean-field approximation in Poisson-Boltzmann theory by comparing its predictions with those of Brownian dynamics simulations. For this purpose we use spherical and cylindrical boundaries and a catenary shape similar to that of the acetylcholine receptor channel. The interior region filled with electrolyte is assumed to have a high dielectric constant, and the exterior region representing protein a low one. Comparisons of the force on a test ion obtained with the two methods show that the shielding effect due to counterions is overestimated in Poisson-Boltzmann theory when the ion is within a Debye length of the boundary. As the ion gets closer to the boundary, the discrepancy in force grows rapidly. The implication for membrane channels, whose radii are typically smaller than the Debye length, is that Poisson-Boltzmann theory cannot be used to obtain reliable estimates of the electrostatic potential energy and force on an ion in the channel environment.

Algorithms↗
PubMed2000

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.

Algorithms↗
PubMed1999

Molecular dynamics study of the KcsA potassium channel.

The structural, dynamical, and thermodynamic properties of a model potassium channel are studied using molecular dynamics simulations. We use the recently unveiled protein structure for the KcsA potassium channel from Streptomyces lividans. Total and free energy profiles of potassium and sodium ions reveal a considerable preference for the larger potassium ions. The selectivity of the channel arises from its ability to completely solvate the potassium ions, but not the smaller sodium ions. Self-diffusion of water within the narrow selectivity filter is found to be reduced by an order of magnitude from bulk levels, whereas the wider hydrophobic section of the pore maintains near-bulk self-diffusion. Simulations examining multiple ion configurations suggest a two-ion channel. Ion diffusion is found to be reduced to approximately 1/3 of bulk diffusion within the selectivity filter. The reduced ion mobility does not hinder the passage of ions, as permeation appears to be driven by Coulomb repulsion within this multiple ion channel.

Bacterial Proteins↗
PubMed1999

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.

Diffusion↗
PubMed1998

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.

Algorithms↗
PubMed1998

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

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

Cations↗
PubMed1998

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.

Biophysical Phenomena↗
PubMed1996

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.

Biophysical Phenomena↗
PubMed1995

Changes in the kinetics and conductance of N-methyl-D-aspartate (NMDA)-receptor activated single channels with temperature.

Single channel currents activated by N-methyl-D-aspartate (NMDA) were recorded in inside-out membrane patches from hippocampal pyramidal cells. The effects of temperature on the conductance level and the kinetics of single channels were accurately determined by using a digital signal processing technique based on hidden Markov models. With increasing temperature, the probability of the channel being in the open state increases steeply, as does the amplitude of channel currents. The proportional change in channel conductance with increasing temperature is greater than the corresponding change in bulk conductivity of 150 mM electrolyte solutions. From the temperature dependence of channel conductance and kinetics, we derived several thermodynamical properties of the NMDA-receptor activated channel. Among these are the height of the energy barrier V presented to an ion traversing the pore, the activation energy Ea and the enthalpy delta H of the channel.

Animals↗
PubMed1994

Temperature dependence of conductivity in electrolyte solutions and ionic channels of biological membranes.

Temperature is a key parameter in the description of any physical system. Experimental study of the temperature dependence of conductivity is very valuable in building and testing theoretical models. This fact does not appear to be fully appreciated and exploited in the study of ionic channels of biological membranes owing in part to the lack of an adequate theory for the temperature dependence of conductivity in electrolyte solutions. To redress this imbalance, and to encourage further temperature-dependence studies in ionic channels, we first give explicit expressions for the conductivity of ions in electrolyte solutions in terms of the microscopic parameters of the liquid. We then propose that the dynamics of ion transport in membrane channels are similar to that in bulk electrolyte solutions, except that ions permeating the pore need to surmount a potential barrier, the height of which can be deduced experimentally. Finally, we use our model to analyze the conductance-temperature relationships obtained in two types of single ionic channels.

Cell Membrane↗