Selectivity of the KcsA potassium channel: Analysis and computation

Mathematical Analysis on Current-Voltage Relations via Classical Poisson-Nernst-Planck Systems with Nonzero Permanent Charges under Relaxed Electroneutrality Boundary Conditions

We focus on a quasi-one-dimensional Poisson-Nernst-Planck model with small permanent charges for ionic flows of two oppositely charged ion species through an ion channel. Of particular interest is to examine the dynamics of ionic flows in terms of I-V (current-voltage) relations with boundary layers due to the relaxation of neutral conditions on boundary concentrations. This is achieved by employing the regular perturbation analysis on the solutions established through geometric singular perturbation analysis. Rich dynamics are observed, particularly, the nonlinear interplays among different physical parameters are characterized. Critical potentials are identified, which play critical roles in the study of ionic flows and can be estimated experimentally. Numerical simulations are performed to further illustrate and provide more intuitive understandings of our analytical results.

Investigating ion transport inside the pentameric ion channel encoded in COVID-19 E protein

Ion flow inside an ion channel can be described through continuum based Born-Poisson-Nernst-Planck (BPNP) equations in conjunction with the Lennard-Jones potential. Keeping in mind the ongoing pandemic, in this study, an attempt has been made to understand the selectivity and the current voltage relation of the COVID-19 E protein pentameric ion channel. Two ionic species, namely Na^{+} and Cl^{-}, have been considered here. E protein is one of the smallest structural protein which is embedded in the outer membrane of the virus. Once the virus is inside the host cell, this protein is expressed abundantly and is responsible for activities such as replication and budding of the virus. In the literature, we can find a few experimental studies focusing on understanding the activity of the channel formed by E proteins of different viruses. Here, we attempt the same study for the COVID-19 E protein ion channel through mathematical modeling. The channel geometry is calculated from the protein data bank file which was provided by NARLabs, Taiwan, using the hole program. Further, it was used to obtain the charge distribution using the pdbtopqr online program. The immersed boundary-lattice Boltzmann method (IB-LBM) has been implemented to numerically solve the system of equations in the channel generated by the protein data bank file. Further, an in-house code which operates on multiple GPUs and uses the cuda platform has been developed to achieve the goal of performing the current investigation.

Electric discharge of electrocytes: Modelling, analysis and simulation

In this paper, we investigate the electric discharge of electrocytes by extending our previous work on the generation of electric potential. We first give a complete formulation of a single cell unit consisting of an electrocyte and a resistor, based on a Poisson-Nernst-Planck (PNP) system with various membrane currents as interfacial conditions for the electrocyte and a Maxwell's model for the resistor. Our previous work can be treated as a special case with an infinite resistor (or open circuit). Using asymptotic analysis, we simplify our PNP system and reduce it to an ordinary differential equation (ODE) based model. Unlike the case of an infinite resistor, our numerical simulations of the new model reveal several distinct features. A finite current is generated, which leads to non-constant electric potentials in the bulk of intracellular and extracellular regions. Furthermore, the current induces an additional action potential (AP) at the non-innervated membrane, contrary to the case of an open circuit where an AP is generated only at the innervated membrane. The voltage drop inside the electrocyte is caused by an internal resistance due to mobile ions. We show that our single cell model can be used as the basis for a system with stacked electrocytes and the total current during the discharge of an electric eel can be estimated by using our model.

Review and Modification of Entropy Modeling for Steric Effects in the Poisson-Boltzmann Equation

The classical Poisson-Boltzmann model can only work when ion concentrations are very dilute, which often does not match the experimental conditions. Researchers have been working on the modification of the model to include the steric effect of ions, which is non-negligible when the ion concentrations are not dilute. Generally the steric effect was modeled to correct the Helmholtz free energy either through its internal energy or entropy, and an overview is given here. The Bikerman model, based on adding solvent entropy to the free energy through the concept of volume exclusion, is a rather popular steric-effect model nowadays. However, ion sizes are treated as identical in the Bikerman model, making an extension of the Bikerman model to include specific ion sizes desirable. Directly replacing the ions of non-specific size by specific ones in the model seems natural and has been accepted by many researchers in this field. However, this straightforward modification does not have a free energy formula to support it. Here modifications of the Bikerman model to include specific ion sizes have been developed iteratively, and such a model is achieved with a guarantee that: (1) it can approach Boltzmann distribution at diluteness; (2) it can reach saturation limit as the reciprocal of specific ion size under extreme electrostatic conditions; (3) its entropy can be derived by mean-field lattice gas model.

Electric potential generation of electrocytes: Modelling, analysis, and computation

In this paper, we developed a one-dimensional model for electric potential generation of electrocytes in electric eels. The model is based on the Poisson-Nernst-Planck system for ion transport coupled with membrane fluxes including the Hodgkin-Huxley type. Using asymptotic analysis, we derived a simplified zero-dimensional model, which we denote as the membrane model in this paper, as a leading order approximation. Our analysis provides justification for the assumption in membrane models that electric potential is constant in the intracellular space. This is essential to explain the superposition of two membrane potentials that leads to a significant transcellular potential. Numerical simulations are also carried out to support our analytical findings.