Abstract:In this article we address the three-dimensional modeling and simulation of biological ion channels using a continuum-based approach. Our multi-physics formulation self-consistently combines, to the best of our knowledge for the first time, ion electrodiffusion, channel fluid motion, thermal self-heating and mechanical deformation. The resulting system of nonlinearly coupled partial differential equations in conservation form is discretized using the Galerkin Finite Element Method. The validation of the proposed computational model is carried out with the simulation of two ion nanochannels. The first is a voltage operated channel with K+ and Na+ ions, and the second is the sodium-potassium pump in which also chlorine (Cl−) and bicarbonate (HCO3−) ions are considered. In the first case study, we investigate the coupling between electrochemical and fluid-dynamical effects. Then, we enrich the modeling picture by investigating the influence of a thermal gradient. Finally, we add a mechanical stress responsible for channel deformation and investigate its effect on the functional response of the channel. Results show that fluid and thermal fields have no influence in absence of mechanical deformation whereas ion distributions and channel functional response are significantly modified if mechanical stress is included in the model. These predictions agree with biophysical conjectures on the importance of protein conformation in the modulation of channel electrochemical properties. In the second case study, we exploit our multiphysical mathematical formulation to investigate the effect of permanent surface charge on the function of the sodium-potassium pump. In particular, we consider the biophysical application in which the pump actively participates to the process of aqueous humor production across the transepithelial membrane in the ciliary body of the eye. In this study we are motivated by the fact that several data are available to verify the accuracy of model predictions and that the role of surface charge has not yet been mathematically analyzed in the specific context at hand. Results show that the model is able to predict in a very accurate manner the correct aqueous humor flow direction, the magnitude of aqueous velocity and the experimentally measured transepithelial membrane potential only if the (negative) surface permanent charge is larger than a limiting value, whereas if the charge is below this value, fluid inversion occurs.

Computational Modeling of Three-Dimensional Electrodiffusion in Biological Systems: Application to the Node of Ranvier

Electrodiffusion Models of Axon and Extracellular Space Using the Poisson-Nernst-Planck Equations

Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries

Finite element simulation of ionicelectrodiffusion in cellular geometries

Modeling and Simulation of Thermo-Fluid-Electrochemical Ion Flow in Biological Channels

Modeling ionic flow between small targets: insights from diffusion and electro-diffusion theory

Spatio-temporal modeling of saltatory conduction in neurons using Poisson-Nernst-Planck treatment and estimation of conduction velocity

A Primitive Model for Predicting Membrane Currents in Excitable Cells Based Only on Ion Diffusion Coefficients

Modeling the voltage distribution in a non-locally but globally electroneutral confined electrolyte medium: applications for nanophysiology

Neural bidomain model for multidimensional ephaptic coupling of neural spike propagation along myelinated fiber bundles with the nodes of Ranvier

A Cell-Based Model for Ionic Electrodiffusion in Excitable Tissue

Three-dimensional simulation of biological ion channels under mechanical, thermal and fluid forces

Ephaptic conduction in a cardiac strand model with 3D electrodiffusion

An electrodiffusive network model with multicompartmental neurons and synaptic connections

Accurate numerical simulation of electrodiffusion and water movement in brain tissue

Modelling in vivo action potential propagation along a giant axon

Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions

Electro-diffusive modeling and the role of spine geometry on action potential propagation in neurons

Electrodiffusion models of synaptic potentials in dendritic spines

Sensitivity of the electrical response of a node of Ranvier model to alterations of the myelin sheath geometry