Collaborative Research: Improved Boundary Element Methods for Electrostatics of Interacting Proteins in Solvent

蛋白质分子是所有生命形式的基本组成部分，静电效应在决定蛋白质结构和功能方面起着关键作用。然而，研究这些影响的实验是困难的，有很大的兴趣，使用计算机模拟，以推进这一领域的生物化学研究。该项目将通过开发新的数学算法来提高这些模拟的准确性和效率。潜在的生物医学应用包括（1）蛋白质-配体结合，这对合成药物设计很重要，以及（2）蛋白质折叠，这对理解错误折叠如何导致疾病很重要。研究人员将与生物科学家合作，该项目开发的软件将以开源格式提供给科学界使用。项目成果将通过会议和研讨会介绍以及在科学期刊上发表文章加以传播。该项目将通过培养博士后、研究生和本科生级别的初级人员来发展计算数学方面的国家劳动力。带电粒子系统的生物化学模拟通常受到计算远程静电相互作用的成本的限制。为了解决这个问题，PI开发了树码加速边界积分（TABI）方法来求解溶剂化蛋白质静电势的泊松-玻尔兹曼（PB）方程。在这个项目中，PI将把TABI求解器扩展到具有挑战性的溶剂中蛋白质相互作用的情况下，专注于微妙的介电和离子结构效应。本论文的主要工作包括：（1）将TABI程序扩展到溶剂化蛋白质-蛋白质复合物结合能的计算中; (2)树码将被调整，以加速计算的总相关函数的参考相互作用网站模型（RISM），它提供了一个物理上准确的描述溶质/溶剂界面附近的水结构。(3)将进行显式离子模拟以验证PB和RISM模拟，并研究有限离子尺寸对界面电荷结构的影响。(4)树码将通过实现更有效的多极验收标准进行改进，并将开发树码的并行CPU和GPU实现。(5)一种新的方法来计算静电PB力将测试在分子动力学模拟相互作用的蛋白质在溶剂中的应用。(6)该奖项反映了NSF的法定使命，并已被认为值得通过使用基金会的智力价值和更广泛的影响审查标准进行评估的支持。

Protein molecules are basic components in all forms of life, and electrostatic effects play a key role in determining protein structure and function. However studying these effects experimentally is difficult, and there is great interest in using computer simulations to advance this area of biochemical research. This project will contribute by developing new mathematical algorithms to improve the accuracy and efficiency of these simulations. Potential biomedical applications include (1) protein-ligand binding, which is important for synthetic drug design, and (2) protein folding, which is important for understanding how misfolding causes disease. The investigators will collaborate with bioscientists, and the software developed in the project will be made available in open source format for use by the scientific community. The project results will be disseminated through conference and seminar presentations, and publications in scientific journals. The project will develop the national workforce in computational mathematics by training junior personnel at the postdoctoral, graduate student, and undergraduate student levels. Biochemical simulations of charged particle systems are often limited by the cost of computing long-range electrostatic interactions. To address this issue the PIs have developed a treecode-accelerated boundary integral (TABI) method to solve the Poisson-Boltzmann (PB) equation for the electrostatic potential of a solvated protein. In this project the PIs will extend the TABI solver to the challenging case of interacting proteins in solvent, focusing on subtle dielectric and ionic structural effects. The work entails advances in modeling, improvements in treecode and boundary element methods, and applications in molecular biochemistry including the following topics: (1) The TABI code will be extended to compute the binding energy of a solvated protein-protein complex. (2) The treecode will be adapted to accelerate computations of the total correlation function in the Reference Interaction Site Model (RISM) which provides a physically accurate description of the water structure near the solute/solvent interface. (3) Explicit ion simulations will be carried out to validate the PB and RISM simulations, and to investigate the effect of finite ion size on the interfacial charge structure. (4) The treecode will be improved by implementing a more effective multipole acceptance criterion, and parallel CPU and GPU implementations of the treecode will be developed. (5) A new method for computing electrostatic PB forces will be tested for application in molecular dynamics simulations of interacting proteins in solvent. (6) Efficient computation of pH-dependent effects in solvated proteins will be demonstrated using TABI.This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.