Electroneutral models for dynamic Poisson-Nernst-Planck systems

Cellular communication among smooth muscle cells: The role of membrane potential via connexins

Communication via action potentials among neurons has been extensively studied. However, effective communication without action potentials is ubiquitous in biological systems, yet it has received much less attention in comparison. Multi-cellular communication among smooth muscles is crucial for regulating blood flow, for example. Understanding the mechanism of this non-action potential communication is critical in many cases, like synchronization of cellular activity, under normal and pathological conditions. In this paper, we employ a multi-scale asymptotic method to derive a macroscopic homogenized bidomain model from the microscopic electro-neutral (EN) model. This is achieved by considering different diffusion coefficients and incorporating nonlinear interface conditions. Subsequently, the homogenized macroscopic model is used to investigate communication in multi-cellular tissues. Our computational simulations reveal that the membrane potential of syncytia, formed by interconnected cells via connexins, plays a crucial role in propagating oscillations from one region to another, providing an effective means for fast cellular communication. Statement of Significance: In this study, we investigated cellular communication and ion transport in vascular smooth muscle cells, shedding light on their mechanisms under normal and abnormal conditions. Our research highlights the potential of mathematical models in understanding complex biological systems. We developed effective macroscale electro-neutral bi-domain ion transport models and examined their behavior in response to different stimuli. Our findings revealed the crucial role of connexinmediated membrane potential changes and demonstrated the effectiveness of cellular communication through syncytium membranes. Despite some limitations, our study provides valuable insights into these processes and emphasizes the importance of mathematical modeling in unraveling the complexities of cellular communication and ion transport.

Piezoionic High Performance Hydrogel Generator and Active Protein Absorber via Microscopic Porosity and Phase Blending

Generating electricity in hydrogel is very important but remains difficult. Hydrogel with electricity generation capability is more capable in bio-relevant tasks such as tissue engineering, artificial skin, or medical treatment, because electricity is indispensable in regulating physiological activities. Here, a porous and phase blending hydrogel structure for effective piezoionic electricity generation is developed. Dynamic electric field is generated taking advantage of the difference in streaming speeds of sodium and chloride in the material. Microscopic porosity and hydrophilic-hydrophobic phase blending are the two key factors for prominent piezoionic performance. Voltages as high as 600 mV are first realized in hydrogels in response to medical ultrasound stimulation. The hydrogel structure is also subjective to effective substance exchange and can actively enrich proteins from surroundings under mechanical stimuli. Preliminary applications in neural stimulation, constructing complex spatial-temporal chemical and electric field distribution patterns, mimetic tactile sensor, sample pretreatment in fast detection, and enzyme immobilization are demonstrated.

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.

Bounds for systems of coupled advection-diffusion equations with application to the Poisson-Nernst-Plank equations

We consider coupled systems of advection-diffusion equations with initial and boundary conditions and determine conditions on the advection terms that allow us to obtain solutions that can be explicitly bounded above and below using the initial and boundary conditions. Given the advection terms, using our methodology one can easily check if such bounds can be obtained and then one can construct the necessary nonlinear transformation to allow the bounds to be determined. We apply this technique to determine bounding quantities for a number of examples. In particular, we show that the three-ion electroneutral Poisson-Nernst-Planck system of equations can be transformed into a system, which allows for the use of our techniques and we determine the bounding quantities. In addition, we determine the general form of advection terms that allow these techniques to be applied and show that our method can be applied to a very wide class of advection-diffusion equations.

A tridomain model for potassium clearance in optic nerve of Necturus

Complex fluids flow in complex ways in complex structures. Transport of water and various organic and inorganic molecules in the central nervous system are important in a wide range of biological and medical processes. However, the exact driving mechanisms are often not known. In this work, we investigate flows induced by action potentials in an optic nerve as a prototype of the central nervous system. Different from traditional fluid dynamics problems, flows in biological tissues such as the central nervous system are coupled with ion transport. They are driven by osmosis created by concentration gradient of ionic solutions, which in turn influence the transport of ions. Our mathematical model is based on the known structural and biophysical properties of the experimental system used by the Harvard group Orkand et al. Asymptotic analysis and numerical computation show the significant role of water in convective ion transport. The full model (including water) and the electrodiffusion model (excluding water) are compared in detail to reveal an interesting interplay between water and ion transport. In the full model, convection due to water flow dominates inside the glial domain. This water flow in the glia contributes significantly to the spatial buffering of potassium in the extracellular space. Convection in the extracellular domain does not contribute significantly to spatial buffering. Electrodiffusion is the dominant mechanism for flows confined to the extracellular domain.

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.

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.

Selectivity of the KcsA potassium channel: Analysis and computation

Ion channels regulate the flux of ions through cell membranes and play significant roles in many physiological functions. Most of the existing literature focuses on computational approaches based on molecular dynamics simulation or numerical solution of the modified Poisson-Nernst-Planck (PNP) system. In this paper, we present an analytical and computational study of a mathematical model of the KcsA potassium channel, including the effects of ion size (Bikerman model) and solvation energy (Born model). Under equilibrium conditions, we obtain an analytical solution of our modified PNP system, which is used to explain selectivity of KcsA of various ions (K^{+}, Na^{+}, Cl^{-}, Ca^{2+}, and Ba^{2+}) due to negative permanent charges inside the filter region and the effect of ion sizes. Our results show that K^{+} is always selected over Na^{+}, as smaller Na^{+} ions have larger solvation energy. As the amount of negative charges in the filter exceeds a critical value, divalent ions (Ca^{2+} and Ba^{2+}) can enter the filter region and block the KcsA channel. For the nonequilibrium cases, due to difficulties associated with a pure analytical or numerical approach, we use a hybrid analytical-numerical method to solve the modified PNP system. Our predictions of selectivity of KcsA channels and saturation phenomenon of the current-voltage (I-V) curve agree with experimental observations.

A Bidomain Model for Lens Microcirculation

There exists a large body of research on the lens of the mammalian eye over the past several decades. The objective of this work is to provide a link between the most recent computational models and some of the pioneering work in the 1970s and 80s. We introduce a general nonelectroneutral model to study the microcirculation in the lens of the eye. It describes the steady-state relationships among ion fluxes, between water flow and electric field inside cells, and in the narrow extracellular spaces between cells in the lens. Using asymptotic analysis, we derive a simplified model based on physiological data and compare our results with those in the literature. We show that our simplified model can be reduced further to the first-generation models, whereas our full model is consistent with the most recent computational models. In addition, our simplified model captures in its equations the main features of the full computational models. Our results serve as a useful link intermediate between the computational models and the first-generation analytical models. Simplified models of this sort may be particularly helpful as the roles of similar osmotic pumps of microcirculation are examined in other tissues with narrow extracellular spaces, such as cardiac and skeletal muscle, liver, kidney, epithelia in general, and the narrow extracellular spaces of the central nervous system, the "brain." Simplified models may reveal the general functional plan of these systems before full computational models become feasible and specific.