[논문 리뷰] Bayesian Estimation of the Hydroxyl Radical Diffusion Coefficient at Low Temperature and High Pressure from Atomistic Molecular Dynamics

Carter T. Butts, Rachel W. Martin|arXiv (Cornell University)|2021. 08. 25.

Spectroscopy and Quantum Chemical Studies참고 문헌 60인용 수 7

한 줄 요약

이 연구는 TIP4P/2005 수분 모델을 사용한 원자적 분자 동역학(MD) 시뮬레이션을 통해 저온·고압 조건에서 수용액 내 하이드록실 라디칼(·OH)의 확산 계수를 추정하기 위해 계층적 베이지안 추론 프레임워크를 제시한다. 저자들은 298 K 및 1 atm에서 실험적 확산 계수 0.23 Ų/ps를 재현하도록 Lennard-Jones 파라미터를 최적화하여 CHARMM 호환 ·OH 힘장(Force field)을 개발하였으며, 유한한 크기 효과를 Yeh와 Hummer의 분석 보정식을 통해 보정한 log-log 다항식 모델을 유도하여 다양한 T 및 P 조건에서 D·OH를 정확하게 예측하였다.

ABSTRACT

The hydroxyl radical is the primary reactive oxygen species produced by the radiolysis of water, and is a significant source of radiation damage to living organisms. Mobility of the hydroxyl radical at low temperatures and/or high pressures is hence a potentially important factor in determining the challenges facing psychrophilic and/or barophilic organisms in high-radiation environments (e.g., ice-interface or undersea environments in which radiative heating is a potential heat and energy source). Here, we estimate the diffusion coefficient for the hydroxyl radical in aqueous solution, using a hierarchical Bayesian model based on atomistic molecular dynamics trajectories in TIP4P/2005 water over a range of temperatures and pressures.

연구 동기 및 목표

• 실험적 측정이 불가능한 저온·고압 조건에서 수용액 내 하이드록실 라디칼(·OH)의 확산 계수를 추정하기 위해.
• TIP4P/2005 수분에서 298 K 및 1 atm 조건에서 알려진 확산 계수 0.23 Ų/ps를 정확히 재현할 수 있는 CHARMM 호환 ·OH 힘장을 개발하기 위해.
• 분자 동역학 시뮬레이션에서의 유한한 크기 효과를 Yeh와 Hummer의 분석 보정식을 사용하여 보정하기 위해.
• 베이지안 사후 추론 기반으로 온도 및 압력에 대한 함수로 D·OH를 예측할 수 있는 log-log 다항식 모델을 구축하기 위해.
• 극한 환경(iced 또는 저해도 해양 환경)에 서식하는 극한성 생물체에서 방사선 손상에 대한 ·OH 이동성의 역할을 정량화하여 분석하기 위해.

제안 방법

• 다양한 시뮬레이션 조건에서의 불확실성을 통합하여 원자적 MD 궤적에서 ·OH 확산 계수를 추론하기 위해 계층적 베이지안 모델을 사용하였다.
• O 및 H 원자에 대해 Lennard-Jones 우물 깊이(ε) 및 반경 절반(rmin/2) 파라미터를 최적화하여, 298 K 및 1 atm 조건에서 실험적 D·OH = 0.23 Ų/ps를 재현하는 CHARMM 기반 ·OH 힘장을 설정하였다.
• 4차원 파라미터 공간(εO, εH, rminO/2, rminH/2)을 탐색하기 위해 250개의 비정규적 Halton 수열 샘플링을 사용하였다.
• 각 시뮬레이션은 NpT 캐노니컬 앙상블을 사용한 1 ns 생산 주기와 온도 제어를 위한 랭지에인 동역학, 그리고 TIP4P/2005 수분에서의 압력 조절을 위한 랭지에인-노제-후우어 피스톤을 적용하였다.
• ·OH 궤적에서 DMD를 계산하기 위해 Bullerjahn 등이 제안한 공분산 기반 추정기법을 적용한 후, Yeh와 Hummer의 공식을 사용하여 유한 크기로 보정된 D∞로 보정: D∞ = DMD + (ζkBT)/(2πηL).
• 수소 결합 에너지 분석을 기반으로 한 파라미터 선택 기준을 적용하여, O–H 기여자 및 O–O 수용체 기하학적 배열 간의 상대적 불안정성 구조를 정확히 반영하는 구성에 우선순위를 두었으며, 물리적 현실성 확보를 위해 설계되었다.

실험 결과

연구 질문

• RQ1 실험 데이터가 없는 저온·고압 조건에서 수분 내 하이드록실 라디칼의 확산 계수는 얼마인가?
• RQ2 기본 조건에서 알려진 확산 계수를 재현할 수 있는 물리적으로 정확한 CHARMM 호환 ·OH 힘장을 어떻게 최적화할 수 있는가?
• RQ3 MD 시뮬레이션에서의 유한한 크기 효과가 추정된 확산 계수에 얼마나 큰 영향을 미치며, 이를 어떻게 보정할 수 있는가?
• RQ4 ·OH 확산 계수는 온도 및 압력에 따라 어떻게 변화하며, 시뮬레이션 데이터 기반으로 예측 모델을 유도할 수 있는가?
• RQ5 ·OH 이동성은 저온성 및 고압성 생물체의 생체분자에 대한 방사선 손상에 어떤 역할을 하는가?

주요 결과

• 최적의 ·OH 힘장 파라미터는 εO = -0.30 kcal/mol, rminO/2 = 1.5 Å, εH = -0.15 kcal/mol, rminH/2 = 1.0 Å이며, 이는 298 K 및 1 atm 조건에서 실험적 D·OH = 0.23 Ų/ps를 정확히 재현한다.
• Yeh와 Hummer의 보정은 평균적으로 추정된 확산 계수를 약 10% 감소시키며, 특히 작은 시뮬레이션 박스에서 보정 폭이 가장 크다.
• 최종 사후 분포 D·OH는 로그-로그 다항식 모델로 가장 잘 기술되며: log10(D·OH) = a + b·log10(T) + c·log10(P) + d·log10(T)·log10(P), 계수는 베이지안 추론를 통해 유도되었다.
• 모델은 온도 감소 및 압력 증가에 따라 D·OH가 크게 감소함을 예측하며, 이는 용매 점도 증가 및 분자 이동성 감소와 일치한다.
• D·OH 추정의 절대 상대 오차(ARE)를 최소화하는 파라미터 조합은 1.3–1.8 Å 범위 내에서 O–O 상호작용의 반발력을 강하게 선호하며, 이는 확산에 대한 공간적 배제 원리의 중요성을 시사한다.
• 수소 결합 분석 결과, 선택된 파라미터 조합은 O–O 대면 구조를 O–H 기여자 구조보다 정확히 불안정화함을 확인하였으며, 이는 현실적인 용매 환경 행동을 보장한다.