skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An efficient finite element iterative method for solving a nonuniform size modified Poisson-Boltzmann ion channel model
In this paper, a nonuniform size modified Poisson-Boltzmann ion channel (nuSMPBIC) model is presented as a nonlinear system of an electrostatic potential and multiple ionic concentrations. It mixes nonlinear algebraic equations with a Poisson boundary value problem involving Dirichlet-Neumann mixed boundary value conditions and a membrane surface charge density to reflect the effects of ion sizes and membrane charges on electrostatics and ionic concentrations. To overcome the difficulties of strong singularities and exponential nonlinearities, it is split into three submodels with a solution of Model 1 collecting all the singular points and Models 2 and 3 much easier to solve numerically than the original nuSMPBIC model. A damped two-block iterative method is then presented to solve Model 3, along with a novel modified Newton iterative scheme for solving each related nonlinear algebraic system. To this end, an effective nuSMPBIC finite element solver is derived and then implemented as a program package that works for an ion channel protein with a three-dimensional molecular structure and a mixture of multiple ionic species. Numerical results for a voltage-dependent anion channel (VDAC) in a mixture of four ionic species demonstrate a fast convergence rate of the damped two-block iterative method, the high performance of the software package, and the importance of considering nonuniform ion sizes. Moreover, the nuSMPBIC model is validated by the anion selectivity property of VDAC.  more » « less
Award ID(s):
2153376
PAR ID:
10413420
Author(s) / Creator(s):
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Journal of Computational Physics
Volume:
470
Issue:
C
ISSN:
0021-9991
Page Range / eLocation ID:
111556
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper presents an efficient finite element iterative method for solving a nonuniform size- modified Poisson-Nernst-Planck ion channel (SMPNPIC) model, along with an SMPNPIC program package that works for an ion channel protein with a three-dimensional crystallographic structure and an ionic solvent with multiple ionic species. In particular, the SMPNPIC model is constructed and then reformulated by novel mathematical techniques so that each iteration of the method only involves linear boundary value problems and nonlinear algebraic systems, circumventing the numerical difficulties caused by the strong nonlinearities, strong asymmetries, and strong differential equation coupling of the SMPNPIC model. To further improve the method’s efficiency, an efficient modified Newton iterative method is adapted to the numerical solution of each related nonlinear algebraic system. In addition, a uniform SMPNPIC model is introduced and solved, numerically, as a special case of the SMPNPIC model. Numerical results for a voltage-dependent anion channel (VDAC) and two mixture solutions of four ionic species, including ATP4− ions, demonstrate the method’s convergence, the package’s high performance, and the importance of considering nonuniform ion size effects. They also partially validate the SMPNPIC model by the anion selectivity property of VDAC. 
    more » « less
  2. ABSTRACT Voltage‐dependent anion channel (VDAC) is the primary conduit for regulated passage of ions and metabolites into and out of a mitochondrion. Calculating the solvation free energy for VDAC is crucial for understanding its stability, function, and interactions within the cellular environment. In this article, numerical schemes for computing the total solvation free energy for VDAC—comprising electrostatic, ideal gas, and excess free energies plus the nonpolar energy—are developed based on a nonuniform size modified Poisson–Boltzmann ion channel (nuSMPBIC) finite element solver along with tetrahedral meshes for VDAC proteins. The current mesh generation package is also updated to improve mesh quality and accelerate mesh generation. A VDAC Solvation Free Energy Calculation (VSFEC) package is then created by integrating these schemes with the updated mesh package, the nuSMPBIC finite element package, the PDB2PQR package, and the OPM database, as well as one uniform SMPBIC finite element package and one Poisson–Boltzmann ion channel (PBIC) finite element package. With the VSFEC package, many numerical experiments are made using six VDAC proteins, eight ionic solutions containing up to four ionic species, including ATP4−and Ca2+, two reference states, different boundary values, and different permittivity constants. The test results underscore the importance of considering nonuniform ionic size effects to explore the varying patterns of the total solvation free energy, and demonstrate the high performance of the VSFEC package for VDAC solvation free energy calculation. 
    more » « less
  3. A single ion channel is a membrane protein with an ion selectivity filter that allows only a single species of ions (such as potassium ions) to pass through in the “open” state. Its selectivity filter also naturally separates a solvent domain into an intracellular domain and an extracellular domain. Such biological and geometrical characteristics of a single ion channel are novelly adopted in the construction of a new kind of dielectric continuum ion channel model, called the Poisson-Nernst-Planck single ion channel (PNPSIC) model, in this paper. An effective PNPSIC finite element solver is then developed and implemented as a software package workable for a single ion channel with a three-dimensional X-ray crystallographic molecular structure and a mixture of multiple ionic species. Numerical results for a potassium channel confirm the convergence and efficiency of the PNPSIC finite element solver and demonstrate the high performance of the software package. Moreover, the PNPSIC model is applied to the calculation of electric current and validated by biophysical experimental data. 
    more » « less
  4. Abstract An integral equation method is presented for the 1D steady-state Poisson-Nernst-Planck equations modeling ion transport through membrane channels. The differential equations are recast as integral equations using Green’s 3rd identity yielding a fixed-point problem for the electric potential gradient and ion concentrations. The integrals are discretized by a combination of midpoint and trapezoid rules, and the resulting algebraic equations are solved by Gummel iteration. Numerical tests for electroneutral and non-electroneutral systems demonstrate the method’s 2nd order accuracy and ability to resolve sharp boundary layers. The method is applied to a 1D model of the K$$^+$$ +  ion channel with a fixed charge density that ensures cation selectivity. In these tests, the proposed integral equation method yields potential and concentration profiles in good agreement with published results. 
    more » « less
  5. Since mitochondria contribute to tumorigenesis and drug resistance in cancer, mitochondrial genetic engineering promises a new direction for cancer therapy. Here, we report the use of the perimitochondrial enzymatic noncovalent synthesis (ENS) of peptides for delivering genes selectively into the mitochondria of cancer cells for mitochondrial genetic engineering. Specifically, the micelles of peptides bind to the voltage-dependent anion channel (VDAC) on mitochondria for the proteolysis by enterokinase (ENTK), generating perimitochondrial nanofibers in cancer cells. This process, facilitating selective delivery of nucleic acid or gene vectors into mitochondria of cancer cells, enables the mitochondrial transgene expression of CRISPR/Cas9, FUNDC1, p53, and fluorescent proteins. Mechanistic investigation indicates that the interaction of the peptide assemblies with the VDAC and mitochondrial membrane potential are necessary for mitochondria targeting. This local enzymatic control of intermolecular noncovalent interactions enables selective mitochondrial genetic engineering, thus providing a strategy for targeting cancer cells. 
    more » « less