Solutions of Poisson's equation in channel-like geometries

doi:10.1016/S0010-4655(98)00090-3

Computer Physics Communications

Volume 115, Issue 1, 1 December 1998, Pages 45-68

https://doi.org/10.1016/S0010-4655(98)00090-3 Get rights and content

Abstract

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 Poisson's 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 Poisson's equation based on an iterative technique. We also provide an analytical solution of Poisson's 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.

References (13)

M. Hoyles et al.
Biophys. J.
(1996)
S. Kuyucak et al.
Biophys. J.
(1998)
S.C. Li et al.
Biophys. J.
(1998)
D.G. Levitt
Biophys. J.
(1978)
B. Hille
Ionic Channels of Excitable Membranes
(1992)
P.M. Morse et al.

There are more references available in the full text version of this article.

Cited by (55)

Comparison of efficient techniques for the simulation of dielectric objects in electrolytes
2015, Journal of Computational Physics
We review two recently developed efficient approaches for the numerical evaluation of the electrostatic polarization potential in particle-based simulations. The first is an image-charge method that can be applied to systems of spherical dielectric objects and provides a closed-form solution of Poisson's equation through multiple image-charge reflections and numerical evaluation of the resulting line integrals. The second is a boundary-element method that computes the discretized surface bound charge through a combination of the generalized minimal residual method (GMRES) and a fast Ewald solver. We compare the accuracy and efficiency of both approaches as a function of the pertinent numerical parameters. We demonstrate use of the image-charge method in a Monte Carlo simulation using the Barnes–Hut octree algorithm and the boundary-element method in a molecular dynamics simulation using the Particle–Particle Particle–Mesh (PPPM) Ewald method, and present numerical results for the ensemble-averaged induced force between two spherical colloids immersed in an electrolyte.
Mechanism of tetrodotoxin block and resistance in sodium channels
2014, Biochemical and Biophysical Research Communications
Citation Excerpt :
The position and velocity of each ion evolves according to a stochastic dynamical system. The electrostatic forces experienced by the ions are derived from pre-calculated lookup tables containing the solutions to Poisson’s equation [16]. The adaptive Poisson–Boltzmann Solver [17] is used to derive the electric field generated by the partial charges in the channel protein.
Tetrodotoxin (TTX) has been used for many decades to characterize the structure and function of biological ion channels. Yet, the precise mechanism by which TTX blocks voltage-gated sodium (Na_V) channels is not fully understood. Here molecular dynamics simulations are used to elucidate how TTX blocks mammalian voltage-gated sodium (Nav) channels and why it fails to be effective for the bacterial sodium channel, Na_VAb. We find that, in Na_VAb, a sodium ion competes with TTX for the binding site at the extracellular end of the filter, thus reducing the blocking efficacy of TTX. Using a model of the skeletal muscle channel, Na_V1.4, we show that the conduction properties of the channel observed experimentally are faithfully reproduced. We find that TTX occludes the entrance of Na_V1.4 by forming a network of hydrogen-bonds at the outer lumen of the selectivity filter. The guanidine group of TTX adopts a lateral orientation, rather than pointing into the filter as proposed previously. The acidic residues just above the selectivity filter are important in stabilizing the hydrogen-bond network between TTX and Na_V1.4. The effect of two single mutations of a critical tyrosine residue in the filter of Na_V1.4 on TTX binding observed experimentally is reproduced using computational mutagenesis.
Monte carlo simulation of electrolyte solutions in biology: In and out of equilibrium
2014, Annual Reports in Computational Chemistry
A concise account is given about Monte Carlo (MC) simulation techniques for homogeneous and inhomogeneous systems of electrolyte solutions at concentrations characteristic of biological situations. While the techniques are quite general, the focus here is on reduced models of ion channels where the electrolyte is modeled in the implicit solvent framework, while the channel and the membrane are modeled using minimal structural information. The ion channels under consideration are the L-type calcium channel (present in muscle cell membranes) and the ryanodine receptor calcium release channel (present in the membrane of the sarcoplasmic reticulum). Although MC simulations are known to be appropriate to simulate equilibrium systems, these examples show how to use MC techniques to simulate nonequilibrium ionic transport. The Local Equilibrium Monte Carlo method is shown to be especially useful to simulate globally nonequilibrium systems by coupling MC to the Nernst–Planck transport equation.
Molecular dynamics simulations of scorpion toxin recognition by the Ca <sup>2+</sup>-activated potassium channel K<inf>Ca</inf>3.1
2013, Biophysical Journal
The Ca²⁺-activated channel of intermediate-conductance (K_Ca3.1) is a target for antisickling and immunosuppressant agents. Many small peptides isolated from animal venoms inhibit K_Ca3.1 with nanomolar affinities and are promising drug scaffolds. Although the inhibitory effect of peptide toxins on K_Ca3.1 has been examined extensively, the structural basis of toxin-channel recognition has not been understood in detail. Here, the binding modes of two selected scorpion toxins, charybdotoxin (ChTx) and OSK1, to human K_Ca3.1 are examined in atomic detail using molecular dynamics (MD) simulations. Employing a homology model of K_Ca3.1, we first determine conduction properties of the channel using Brownian dynamics and ascertain that the simulated results are in accord with experiment. The model structures of ChTx-K_Ca3.1 and OSK1-K_Ca3.1 complexes are then constructed using MD simulations biased with distance restraints. The ChTx-K_Ca3.1 complex predicted from biased MD is consistent with the crystal structure of ChTx bound to a voltage-gated K⁺ channel. The dissociation constants (K_d) for the binding of both ChTx and OSK1 to K_Ca3.1 determined experimentally are reproduced within fivefold using potential of mean force calculations. Making use of the knowledge we gained by studying the ChTx-K_Ca3.1 complex, we attempt to enhance the binding affinity of the toxin by carrying out a theoretical mutagenesis. A mutant toxin, in which the positions of two amino acid residues are interchanged, exhibits a 35-fold lower K_d value for K_Ca3.1 than that of the wild-type. This study provides insight into the key molecular determinants for the high-affinity binding of peptide toxins to K_Ca3.1, and demonstrates the power of computational methods in the design of novel toxins.
Ion conduction in ligand-gated ion channels- Brownian dynamics studies of four recent crystal structures
2010, Biophysical Journal
Four x-ray crystal structures of prokaryotic homologs of ligand-gated ion channels have recently been determined: ELIC from Erwinia chrysanthemi, two structures of a proton-activated channel from Gloebacter violaceus (GLIC1 and GLIC2) and that of the E221A mutant (GLIC1M). The availability of numerous structures of channels in this family allows for aspects of channel gating and ion conduction to be examined. Here, we determine the likely conduction states of the four structures as well as IV curves, ion selectivity, and steps involved in ion permeation by performing extensive Brownian dynamics simulations. Our results show that the ELIC structure is indeed nonconductive, but that GLIC1 and GLIC1M are both conductive of ions with properties different from those seen in experimental studies of the channel. GLIC2 appears to reflect an open state of the channel with a predicted conductance of 10.8–12.4 pS in 140 mM NaCl solution, which is comparable to the experimental value 8 ± 2 pS. The extracellular domain of the channel is shown to have an important influence on the channel current, but a less significant role in ion selectivity.
Role of acetylcholine receptor domains in ion selectivity
2009, Biochimica et Biophysica Acta - Biomembranes
Citation Excerpt :
The total force acting on each and every ion in the assembly is calculated and then new positions are determined for the ions a short time later. Electrostatic forces are calculated by assigning dielectric constants of 2 to the protein and 60 to the water in the channel and solving Poisson's equation using an iterative method [45]. It should be noted that while the dielectric constant of bulk water is close to 80, this is likely to be reduced in the confined space inside the pore.
The nicotinic acetylcholine receptor (nAChR) is a ligand gated ion channel protein, composed of three domains: a transmembrane domain (TM-domain), extracellular domain (EC-domain), and intracellular domain (IC-domain). Due to its biological importance, much experimental and theoretical research has been carried out to explore its mechanisms of gating and selectivity, but there are still many unresolved issues, especially on the ion selectivity. Moreover, most of the previous theoretical work has concentrated on the TM-domain or EC-domain of nAChR, which may be insufficient to understand the entire structure–function relation. In this work, we perform molecular dynamics, Brownian dynamics simulations and continuum electrostatic calculations to investigate the role of different nAChR domains in ion conduction and selectivity. The results show that although both the EC and IC domains contain strong negative charges that create large cation concentrations at either end of the pore, this alone is not sufficient to create the observed cation selectivity and may play a greater role in determining the channel conductance. The presence of cations in the wide regions of the pore can screen out the protein charge allowing anions to enter, meaning that local regions of the TM-domain are most likely responsible for discriminating between ions. These new results complement our understanding about the ion conduction and selectivity mechanism of nAChR.

View all citing articles on Scopus

View full text

Solutions of Poisson's equation in channel-like geometries

Abstract

Biophys. J.

Biophys. J.

Biophys. J.

Biophys. J.

Ionic Channels of Excitable Membranes