[논문 리뷰] Small ionic radii limit time step in Martini 3 molecular dynamics simulations

Balázs Fábián, Sebastian Thallmair외 1명|arXiv (Cornell University)|2022. 04. 11.

Protein Structure and Dynamics참고 문헌 46인용 수 16

한 줄 요약

이 논문은 마티니 3에서 Na⁺ 및 Cl⁻ 이온에 대해 '티니'(TQ5) 구슬을 사용할 경우 이온-물 및 이온-이온 충돌 불안정성으로 인해 최대 안정 시간 단계가 25 fs로 제한됨을 규명한다. 더 긴 시간 단계를 가능하게 하기 위해 저자들은 이온 질량을 증가시키거나 구슬 크기를 '스몰'(SQ5) 또는 '레귤러'(RQ5)로 업그레이드하는 것을 제안하며, 이는 구조적 및 동적 성질에 미치는 영향을 최소화하면서도 30–40 fs에서 안정적인 시뮬레이션을 가능하게 한다. 이는 샘플링 효율을 크게 향상시킨다.

ABSTRACT

Among other improvements, the Martini 3 coarse-grained force field provides a more accurate description of the solvation of protein pockets and channels through the consistent use of various bead types and sizes. Here, we show that the representation of Na$^+$ and Cl$^-$ ions as "tiny" (TQ5) beads limits the accessible time step to 25 fs. By contrast, with Martini 2, time steps of 30-40 fs were possible for lipid bilayer systems without proteins. This limitation is relevant for, e.g., phase separating lipid mixtures that require long equilibration times. We derive a quantitative kinetic model of time-integration instabilities in molecular dynamics (MD) as a function of time step, ion concentration and mass, system size, and simulation time. With this model, we demonstrate that ion-water interactions are the main source of instability at physiological conditions, followed closely by ion-ion interactions. We show that increasing the ionic masses makes it possible to use time steps up to 40 fs with minimal impact on static equilibrium properties and on dynamical quantities such as lipid and ion diffusion coefficients. Increasing the size of the bead representing the ions (and thus changing their hydration) also permits longer time steps. The use of larger time steps in Martini 3 simulations results in a more efficient exploration of configuration space. The kinetic model of MD simulation crashes can be used to determine the maximum allowed time step whenever sampling efficiency is critical.

연구 동기 및 목표

• 마티니 3 시뮬레이션에서 Na⁺ 및 Cl⁻ 이온을 사용할 경우 시간 단계 불안정성의 원인을 규명하기 위해.
• 마티니 3에서 25 fs 시간 단계로 제한되는 문제를 해결하여 장시간 스케일 시뮬레이션의 효율적 샘플링을 가능하게 하기 위해.
• 시간 단계, 이온 농도 및 이온 질량의 함수로 MD 장애 빈도를 정량적으로 모델링하기 위해.
• 더 긴 안정적인 시간 단계를 확보하면서 시스템 성질에 영향을 최소화하는 두 가지 해결책—이온 질량 증가 및 구슬 유형 변경—을 평가하기 위해.
• 고속도 MD 시뮬레이션에서 최적의 시간 단계를 결정하는 일반적인 프레임워크를 제공하기 위해.

제안 방법

• MD 장애 빈도 kcrash(cion, ∆t)의 동역학 모델을 유도하고, 시뮬레이션 장애 통계에 대해 전역적으로 피팅하여 불안정성 원인을 정량화하였다.
• 모델을 사용해 이온-물 및 이온-이온 상호작용이 시간 적분 불안정성에 기여하는 정도를 분석하였다.
• 고전적 분할 함수를 유지하면서도 동역학에 미치는 영향을 최소화하기 위해 이온 질량을 균일하게 확대(m′ = 2m 또는 4m)하였다.
• 구슬 유형을 '티니'(TQ5)에서 '스몰'(SQ5) 또는 '레귤러'(RQ5)로 변경하여 마티니 3의 다중 해상도 철학에 부합하는 수화 구조를 재정의하였다.
• 다양한 시간 단계와 이온 매개변수에서 확산 계수, 오더 파라미터, 이온 분포 및 세포막 유동성 등의 시스템 성질을 모니터링하였다.
• 모델은 나트륨 클로라이드 용액, DPPC 이중층 및 PI-(4,5)P2 함유 막을 포함한 여러 시스템에서 검증되었으며, 일관된 결과를 보였다.

실험 결과

연구 질문

• RQ1 Na⁺ 및 Cl⁻ 이온에 대해 '티니'(TQ5) 구슬을 사용할 경우 마티니 3 시뮬레이션에서 시간 단계가 무엇에 의해 제한되는가?
• RQ2 이온-물 및 이온-이온 상호작용은 MD 시간 적분 불안정성에 어떻게 기여하는가?
• RQ3 이온 질량을 증가시키면 25 fs를 초과하는 시간 단계에서 안정적인 시뮬레이션을 수행할 수 있으며, 평형 및 동적 성질에 미치는 영향은 최소한인가?
• RQ4 이온의 구슬 크기를 TQ5에서 SQ5 또는 RQ5로 변경하면 시스템 거동을 유지하면서 더 긴 시간 단계를 가능하게 하는가?
• RQ5 유도된 장애 빈도 동역학 모델을 사용해 고속도 MD 캠프에 최적의 시간 단계를 결정할 수 있는가?

주요 결과

• TQ5 이온을 사용하는 마티니 3 시뮬레이션에서 최대 안정 시간 단계는 이온-물 및 이온-이온 충돌로 인한 불안정성으로 인해 25 fs로 제한된다.
• 이온 질량을 두 배로 늘리면 30 fs에서 안정적인 시뮬레이션을 수행할 수 있고, 네 배로 늘리면 40 fs에서도 안정적이며, 확산 계수 및 구조적 성질에 미치는 영향은 미미하다.
• TQ5에서 SQ5 또는 RQ5로 구슬 유형을 변경하면 마티니 3의 다중 해상도 설계 원칙에 부합하면서도 30–40 fs에서 안정적인 시뮬레이션을 가능하게 한다.
• 동역학적 장애 모델은 불안정성 추세를 정확히 예측할 수 있으며, 고속도 시뮬레이션에 최적의 시간 단계를 결정하는 데 사용될 수 있다.
• 지질의 확산 계수 및 오더 파라미터는 시간 단계와 이온 질량 또는 구슬 유형의 변경에 관계없이 거의 변화하지 않으며, 이는 최소한의 왜곡이 있음을 시사한다.
• 40 fs 시간 단계는 RQ5 이온이 있는 시스템에서만 안정적이며, PI-(4,5)P2 함유 막과 같은 복잡한 시스템에는 권장되지 않는다.