Collaborative Research: Advances in Nonlocal Dielectric Modeling and Free Energy Calculation for Protein in Ionic Solvent

非局域介电模型考虑了水分子间的极化相关性，对经典的泊松介电模型进行了改进。然而，由于离子溶剂中出现的建模和算法复杂性，目前对该方法的研究大多仅限于水溶剂。目前确实存在的离子模型未能纳入关键的非局部介电效应。最近的发展也表明，熵的变化，由于蛋白质配体协会是至关重要的了解结合亲和力。然而，熵的计算仍然是一项非常困难的任务。受这些挑战的启发，本项目旨在开发新的非局部连续静电模型和新的数值求积，用于直接计算蛋白质在离子溶剂中的熵和自由能。新的非局部模型将构造从一个新的组合的非局部介电方法与硬球混合流体的基本测量理论下的约束功能优化协议。与经典的Poisson-Boltzmann方程相比，它们有望显着提高静电势计算的准确性，因为它们反映了水分子之间的离子尺寸效应和极化相关性。新的数值求积将开发通过使用一个特殊的棱柱形单元插值构造从棱柱形网格的有界状态区域附近的势能最小点。通过采用一种新的非局部模型计算所涉及的静电势，进一步降低了计算复杂度。最后，将开发新的快速数值算法和程序包来求解新的电介质模型和实现新的数值求积。 蛋白质在离子溶剂中的静电势能、熵和自由能的计算是生物分子模拟中的一项基本任务。新的非局部介电模型，计算熵和自由能的数值求积，以及伴随的有效的数值算法和程序包，从这个项目产生的将是一个相当大的贡献，数学生物学，计算生物化学，计算数学和计算机科学领域。它们将在离子通道研究、合理药物设计和其他生物工程应用中发挥重要作用，并将提高我们对关键生理过程、细胞能量学和蛋白质-配体结合亲和力以及一般健康和疾病的定量理解。 该项目的发现预计将对数学，计算机科学，生物化学和生物工程的发展产生重大影响。 由于最近人们对解决高维问题的算法产生了浓厚的兴趣，因此该项目中所取得的任何进展都将在与具有指数衰减的被积函数的高维积分相关的各种领域产生广泛的潜在影响。

The nonlocal dielectric approach can significantly enhance the classic Poisson dielectric model by considering the polarization correlations among water molecules. However, current studies on the approach are mostly restricted to the water solvent, due to modeling and algorithmic complications that arise in the case of ionic solvents. The current ionic models that do exist fail to incorporate crucial nonlocal dielectric effects. Recent developments also indicate that entropic changes due to protein-ligand association are critical to understanding binding affinity. The computation of entropy, however, remains a very difficult task. Motivated by these challenges, this project aims to develop new nonlocal continuum electrostatic models and new numerical quadratures for the direct calculation of entropy and free energy for protein in ionic solvent. The new nonlocal models will be constructed from a novel combination of the nonlocal dielectric approach with the fundamental measure theory of hard-sphere mixture fluids under the constrained functional optimization protocol. They are expected to significantly improve the accuracy of electrostatic potential calculations in comparison to the classic Poisson-Boltzmann equation, since they reflect both ionic size effects and polarization correlations among water molecules. The new numerical quadratures will be developed by using a special prismatic element interpolation constructed from a prismatic mesh of a bounded state region near a potential energy minimum point. The computing complexity will be further reduced through using a new nonlocal model for computing involved electrostatic potential. Lastly, new fast numerical algorithms and program packages will be developed for solving the new dielectric models and for implementing the new numerical quadratures. Calculation of electrostatic potential energy, entropy, and free energy for protein in ionic solvent is a fundamental task in biomolecular simulations. The new nonlocal dielectric models, numerical quadratures for computing entropy and free energy, and the accompanying efficient numerical algorithms and program packages produced from this project will be a considerable contribution to the fields of mathematical biology, computational biochemistry, computational mathematics, and computer science. They will play important roles in ion channel studies, rational drug design, and other bioengineering applications, and will improve our quantitative understanding of critical physiological processes, cellular energetics, and affinity in protein-ligand binding, and of health and disease in general. The findings from this project are expected to have a significant impact on the development of mathematics, computer science, biochemistry, and bioengineering. Because there has been substantial interest recently in algorithms for solving high-dimensional problems, any advances made in this project will have broad potential impact in a variety of areas that are related to high-dimensional integrals with integrands that decay exponentially.