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

Riccardo Sacco*, Paolo Airoldi, Aurelio G. Mauri, Joseph W. Jerome

*Corresponding author for this work

Research output: Contribution to journalArticle

3 Citations (Scopus)

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.

Original languageEnglish (US)
Pages (from-to)221-251
Number of pages31
JournalApplied Mathematical Modelling
Volume43
DOIs
StatePublished - Mar 1 2017

Fingerprint

Ion Channels
Charge
Pump
Fluid
Mechanical Stress
Three-dimensional
Fluids
Functional Response
Ions
Potassium
Sodium
Simulation
Pumps
Surface charge
Galerkin Finite Element Method
Membrane Potential
Formulation
Multiphysics
Conformation
Computational Model

Keywords

  • Finite element method
  • Fluid
  • Interface conditions
  • Ionic transport
  • Mechanics
  • Numerical simulation

ASJC Scopus subject areas

  • Modeling and Simulation
  • Applied Mathematics

Cite this

@article{8825afc4685c454a884eaf0417e52bef,
title = "Three-dimensional simulation of biological ion channels under mechanical, thermal and fluid forces",
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.",
keywords = "Finite element method, Fluid, Interface conditions, Ionic transport, Mechanics, Numerical simulation",
author = "Riccardo Sacco and Paolo Airoldi and Mauri, {Aurelio G.} and Jerome, {Joseph W.}",
year = "2017",
month = "3",
day = "1",
doi = "10.1016/j.apm.2016.10.053",
language = "English (US)",
volume = "43",
pages = "221--251",
journal = "Applied Mathematical Modelling",
issn = "0307-904X",
publisher = "Elsevier Inc.",

}

Three-dimensional simulation of biological ion channels under mechanical, thermal and fluid forces. / Sacco, Riccardo; Airoldi, Paolo; Mauri, Aurelio G.; Jerome, Joseph W.

In: Applied Mathematical Modelling, Vol. 43, 01.03.2017, p. 221-251.

Research output: Contribution to journalArticle

TY - JOUR

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

AU - Sacco, Riccardo

AU - Airoldi, Paolo

AU - Mauri, Aurelio G.

AU - Jerome, Joseph W.

PY - 2017/3/1

Y1 - 2017/3/1

N2 - 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.

AB - 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.

KW - Finite element method

KW - Fluid

KW - Interface conditions

KW - Ionic transport

KW - Mechanics

KW - Numerical simulation

UR - http://www.scopus.com/inward/record.url?scp=85006988859&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85006988859&partnerID=8YFLogxK

U2 - 10.1016/j.apm.2016.10.053

DO - 10.1016/j.apm.2016.10.053

M3 - Article

AN - SCOPUS:85006988859

VL - 43

SP - 221

EP - 251

JO - Applied Mathematical Modelling

JF - Applied Mathematical Modelling

SN - 0307-904X

ER -