[논문 리뷰] A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics

Jaydeep P. Bardhan, D. A. Tejani외 3명|arXiv (Cornell University)|2015. 05. 24.

Nanopore and Nanochannel Transport Studies참고 문헌 34인용 수 3

한 줄 요약

이 논문은 연속체 전기역학 모델에 대해 전하 기호 비대칭성—용매 구조와 유한한 용매 크기에서 기인하는 핵심 물리적 효과—을 반영하기 위해 표준 맥스웰 경계 조건을 수정한 비선형 경계 조건(NLBC)을 제안한다. 검증된 NLBC 모델은 전체 원자 분자 동역학 시뮬레이션과 비교하여 표면 전하에서의 비대칭 전기적 자유 에너지를 정확히 재현하며, 특히 용매-용질 인터페이스 근처에서 높은 정확도를 보이고, 용매 영역에서 가우스 법칙을 준수한다.

ABSTRACT

Understanding the behavior of biomolecules such as proteins requires understanding the critical influence of the surrounding fluid (solvent) environment--water with mobile salt ions such as sodium. Unfortunately, for many studies, fully atomistic simulations of biomolecules, surrounded by thousands of water molecules and ions are too computationally slow. Continuum solvent models based on macroscopic dielectric theory (e.g. the Poisson equation) are popular alternatives, but their simplicity fails to capture well-known phenomena of functional significance. For example, standard theories predict that electrostatic response is symmetric with respect to the sign of an atomic charge, even though response is in fact strongly asymmetric if the charge is near the biomolecule surface. In this work, we present an asymmetric continuum theory that captures the essential physical mechanism--the finite size of solvent atoms--using a nonlinear boundary condition (NLBC) at the dielectric interface between the biomolecule and solvent. Numerical calculations using boundary-integral methods demonstrate that the new NLBC model reproduces a wide range of results computed by more realistic, and expensive, all-atom molecular-dynamics (MD) simulations in explicit water. We discuss model extensions such as modeling dilute-electrolyte solvents with Debye-Huckel theory (the linearized Poisson-Boltzmann equation) and opportunities for the electromagnetics community to contribute to research in this important area of molecular nanoscience and engineering.

연구 동기 및 목표

• 생물분자의 표면 근처에서 전하 기호 비대칭성을 포착하지 못하는 연속체 전기역학 모델의 오랜 한계를 해결한다.
• 표준 푸아송 및 푸아송-볼츠만 이론이 예측하는 대칭적 반응이 실험 및 원자 수준의 시뮬레이션 관측과 다름을 해결한다.
• 용매 분자의 유한한 크기와 계면 수분 구조 효과를 포함하는 물리적으로 타당한 비선형 경계 조건을 개발한다.
• 용매 영역에서 가우스 법칙을 준수함으로써 기본 전기역학과의 일관성을 확보한다.
• 특히 표면에 위치한 이온과 아미노산 잔기의 경우에 대해 정확하고 계산 효율적인 전기적 자유 에너지 예측을 가능하게 한다.

제안 방법

• 용질-용매 인터페이스에서 표준 맥스웰 경계 조건(정규 전기력선의 연속성)을 비선형 경계 조건(NLBC)으로 대체한다.
• 복잡한 분자 기하구조를 위한 푸아송 방정식을 해결하기 위해 경계 적분 방정식(BIE)을 사용하여 전기적 문제를 수립한다.
• 전하 크기의 조각별 선형 함수로 용매 반응을 모델링하며, 영점에서의 불연속성은 계면 수분 구조 효과를 반영한다.
• 희석 전해질에서의 선형화된 푸아송-볼츠만 방정식(LPBE)을 통해 스크리닝을 고려하여 그린 함수를 수정함으로써 외부 용매 영역에서 가우스 법칙을 준수한다.
• 재현 가능성을 위해 MATLAB에 경계 요소 방법 솔버를 구현하고, 공개된 소스 코드와 데이터를 제공한다.
• 기준 시스템 테스트를 위해 단원자 이온과 아미노산을 위한 구형 및 타원형 기하구조를 근사치로 사용한다.

실험 결과

연구 질문

• RQ1 비선형 경계 조건이 전기적 자유 에너지의 전하 기호 비대칭성을 전자기적 분자 동역학 시뮬레이션에서 관측된 바와 같이 정확히 재현할 수 있는가?
• RQ2 용매 영역에서 가우스 법칙을 준수하는 것이 전하를 띤 생물분자의 연속체 전기역학 모델 정확도에 어떤 영향을 미치는가?
• RQ3 전하 기호 비대칭성의 크기가 용질-용매 인터페이스에서의 거리에 따라 어느 정도 의존하는가?
• RQ4 NLBC 모델이 1가 이온과 아미노산 잔기에 대해 고정밀도와 계산 효율성으로 충전 자유 에너지를 정량적으로 재현할 수 있는가?
• RQ5 용매 스크리닝(데바이-후크스/선형화된 푸아송-볼츠만)은 NLBC 모델의 성능과 일관성에 어떤 영향을 미치는가?

주요 결과

• 새로운 NLBC 모델은 원래 대칭 모델과 표준 푸아송 이론보다도 단일가 이온의 충전 자유 에너지를 훨씬 더 정확하게 재현한다.
• 모델은 표면 전하에서 상당한 전하 기호 비대칭성을 포착하며, 단백질 크기의 구에서조차도 양전하와 음전하 간의 에너지 차이가 크게 유지된다.
• 비대칭성의 크기는 용질-용매 인터페이스에서의 거리가 증가함에 따라 급격히 감소하며, 원자 수준의 시뮬레이션과 물리적 직관과 일치한다.
• 용매 영역에서 가우스 법칙을 준수하는 것은 모델의 일관성을 향상시키며, 특히 선형화된 푸아송-볼츠만 방정식을 사용해 전해질 용액을 모델링할 경우에 매우 중요하다.
• 아미노산 잔기의 경우, 전하가 단백질 표면에 가까울 때 비대칭 전기적 반응이 가장 두드러지며, 깊이가 깊어질수록 감소한다.
• NLBC 모델은 물리적으로 타당한 경계 조건을 통해 계면 수분 효과를 정확히 포착하는 첫 번째 비대칭 푸아송 기반 연속체 이론이며, 고정밀 전자기적 분자 동역학 시뮬레이션을 통해 검증되었다.