Search

The Voltage-Dependent Anion Channel (VDAC) protein is the primary conduit for the regulated passage of ions and metabolites into and out of mitochondria. Calculating its solvation free energy is crucial for understanding its stability, function, and interactions within the cellular environment. In this paper, we introduce a total solvation free energy, $E$, which is the sum of electrostatic, ideal gas, and excess free energies, along with a non-polar energy to yield a zero of $E$ in the absence of charges. We develop numerical schemes for computing $E$ and update the current mesh generation package to accelerate the generation of tetrahedral meshes and improve the quality of meshes for computing $E$. By integrating these schemes and the updated mesh package with our non-uniform size modified Poisson-Boltzmann ion channel (nuSMPBIC), SMPBIC, and PBIC finite element packages, the PDB2PQR package, and the OPM database, we create the VDAC Solvation Free Energy Calculation (VSFEC) package. Using the VSFEC package, we perform comparison tests on the nuSMPBIC, SMPBIC, and PBIC models by using six VDAC proteins and various ionic solutions containing up to four ionic species, including ATP$^{4-}$ and Ca$^{2+}$. We also conduct tests with different neutral voltages and permittivity constants to explore the varying patterns of $E$. Our test results underscore the importance of considering non-uniform ionic size effects and demonstrate the high performance of the VSFEC package in calculating the solvation free energy of VDAC proteins., Comment: 22 pages, 7 figures, 3 tables

This paper presents an efficient finite element iterative method for solving a nonuniform size-modified Poisson-Nernst-Planck ion channel (SMPNPIC) model, along with a 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. Numerical results for a voltage-dependent anion channel (VDAC) and a mixture solution of four ionic species 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., Comment: 24 pages, 5 figures, 3 tables

Pulmonary artery-vein segmentation is crucial for disease diagnosis and surgical planning and is traditionally achieved by Computed Tomography Pulmonary Angiography (CTPA). However, concerns regarding adverse health effects from contrast agents used in CTPA have constrained its clinical utility. In contrast, identifying arteries and veins using non-contrast CT, a conventional and low-cost clinical examination routine, has long been considered impossible. Here we propose a High-abundant Pulmonary Artery-vein Segmentation (HiPaS) framework achieving accurate artery-vein segmentation on both non-contrast CT and CTPA across various spatial resolutions. HiPaS first performs spatial normalization on raw CT volumes via a super-resolution module, and then iteratively achieves segmentation results at different branch levels by utilizing the lower-level vessel segmentation as a prior for higher-level vessel segmentation. We trained and validated HiPaS on our established multi-centric dataset comprising 1,073 CT volumes with meticulous manual annotations. Both quantitative experiments and clinical evaluation demonstrated the superior performance of HiPaS, achieving an average dice score of 91.8% and a sensitivity of 98.0%. Further experiments showed the non-inferiority of HiPaS segmentation on non-contrast CT compared to segmentation on CTPA. Employing HiPaS, we have conducted an anatomical study of pulmonary vasculature on 11,784 participants in China (six sites), discovering a new association of pulmonary vessel anatomy with sex, age, and disease states: vessel abundance suggests a significantly higher association with females than males with slightly decreasing with age, and is also influenced by certain diseases, under the controlling of lung volumes.

In this paper, a deep learning method for solving an improved one-dimensional Poisson-Nernst-Planck ion channel (PNPic) model, called the PNPic deep learning solver, is presented. In particular, it combines a novel local neural network scheme with an effective PNPic finite element solver. Since the input data of the neural network scheme only involves a small local patch of coarse grid solutions, which the finite element solver can quickly produce, the PNPic deep learning solver can be trained much faster than any corresponding conventional global neural network solvers. After properly trained, it can output a predicted PNPic solution in a much higher degree of accuracy than the low cost coarse grid solutions and can reflect different perturbation cases on the parameters, ion channel subregions, and interface and boundary values, etc. Consequently, the PNPic deep learning solver can generate a numerical solution with high accuracy for a family of PNPic models. As an initial study, two types of numerical tests were done by perturbing one and two parameters of the PNPic model, respectively, as well as the tests done by using a few perturbed interface positions of the model as training samples. These tests demonstrate that the PNPic deep learning solver can generate highly accurate PNPic numerical solutions., Comment: 17 pages, 4 figures, 5 tables