A lattice relaxation algorithm is developed to solve the Poisson-Nernst-Planck (PNP) equations for ion transport through arbitrary three-dimensional volumes. Calculations of systems characterized by simple parallel plate and cylindrical pore geometries are presented in order to calibrate the accuracy of the method. A study of ion transport through gramicidin A dimer is carried out within this PNP framework. Good agreement with experimental measurements is obtained. Strengths and weaknesses of the PNP approach are discussed.

Paracellular ion transport in epithelia is mediated by pores formed by members of the claudin family. The degree of selectivity and the molecular mechanism of ion permeation through claudin pores are poorly understood. By expressing a high-conductance claudin isoform, claudin-2, in high-resistance Madin-Darby canine kidney cells under the control of an inducible promoter, we were able to quantitate claudin pore permeability. Claudin-2 pores were found to be narrow, fluid filled, and cation selective. Charge selectivity was mediated by the electrostatic interaction of partially dehydrated permeating cations with a negatively charged site within the pore that is formed by the side chain carboxyl group of aspartate-65. Thus, paracellular pores use intrapore electrostatic binding sites to achieve a high conductance with a high degree of charge selectivity.

No abstract

Nuclear Pore Complexes (NPCs) are key cellular transporter that control nucleocytoplasmic transport in eukaryotic cells, but its transport mechanism is still not understood. The centerpiece of NPC transport is the assembly of intrinsically disordered polypeptides, known as FG nucleoporins, lining its passageway. Their conformations and collective dynamics during transport are difficult to assess in vivo. In vitro investigations provide partially conflicting results, lending support to different models of transport, which invoke various conformational transitions of the FG nucleoporins induced by the cargo-carrying transport proteins. We show that the spatial organization of FG nucleoporin assemblies with the transport proteins can be understood within a first principles biophysical model with a minimal number of key physical variables, such as the average protein interaction strengths and spatial densities. These results address some of the outstanding controversies and suggest how molecularly divergent NPCs in different species can perform essentially the same function.

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.

A recently introduced real-space lattice methodology for solving the three-dimensional Poisson-Nernst-Planck equations is used to compute current-voltage relations for ion permeation through the gramicidin A ion channel embedded in membranes characterized by surface dipoles and/or surface charge. Comparisons to a variety of experimental results, presented herein, have proven largely successful. Strengths and weaknesses of the method are discussed.

A dynamic lattice Monte Carlo (DLMC) simulation approach to the description of ion transport in dielectric environments is presented. Conventional approaches using periodic boundary conditions are inefficient for nonequilibrium situations in inhomogeneous systems. Instead, the simulated system is embedded in a bigger system that determines the average electrostatic potential and the ionic concentrations at its boundaries. Two issues are of special importance: implementing the given boundary conditions in the treatment of dynamical processes at and near the boundaries, and efficient evaluation of ion-ion interaction in the heterogeneous dielectric medium during the Monte Carlo simulation. The performance of the method is checked by comparing numerical results to exact solutions for simple geometries, and to mean field (Poisson-Nernst-Planck, PNP) theory in a system where the latter should provide a reasonable description. Other examples in which the PNP theory fails in various degrees are shown and discussed. In particular, PNP results deviate considerably from the DLMC dynamics for ion transport through rigid narrow membrane channels with large disparity between the dielectric constants of the protein and the water environments.

The local diffusion constant of K + inside the Gramicidin A (GA) channel has been calculated using four computational methods based on molecular dynamics (MD) simulations, specifically: Mean Square Displacement (MSD), Velocity Autocorrelation Function (VACF), Second Fluctuation Dissipation Theorem (SFDT) and analysis of the Generalized Langevin Equation for a Harmonic Oscillator (GLE-HO). All methods were first tested and compared for K + in bulk water-all predicted the correct diffusion constant. Inside GA, MSD and VACF methods were found to be unreliable because they are biased by the systematic force exerted by the membrane-channel system on the ion. SFDT and GLE-HO techniques properly unbias the influence of the systematic force on the diffusion properties and predicted a similar diffusion constant of K + inside GA, namely, ca. 10 times smaller than in the bulk. It was found that both SFDT and GLE-HO methods require extensive MD sampling on the order of tens of nanoseconds to predict a reliable diffusion constant of K + inside GA.