[논문 리뷰] Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

John M. A. Grime, Jesper J. Madsen|arXiv (Cornell University)|2019. 10. 11.

Lipid Membrane Structure and BehaviorBiochemistry, Genetics and Molecular Biology참고 문헌 49인용 수 18

한 줄 요약

이 논문은 다양한 척도와 해상도에서 조절 가능한 지질 집합체를 시뮬레이션하기 위한 계산적으로 효율적인 입자 기반의 군집화 모델을 제시한다. 부드럽고 짧은 범위의 두체 상호작용 잠재력과 모듈식 베드 표현 방식(지질당 2~5개의 베드, 추가로 1.5개의 베드로 구성된 준단층 모델)을 사용하여, 조절 가능한 상 행동을 보이며 밀레르, 이중층, 고세포를 정확하게 형성한다. 이는 기존 모델에 비해 향상된 수치적 안정성과 역학을 보이며, 굽힘과 상 전이와 같은 핵심 세포막 물리 법칙을 잘 캐릭터라이즈한다.

ABSTRACT

We present a minimal model for simulating dynamics of assorted lipid assemblies in a computationally efficient manner. Our model is particle-based and consists of coarse-grained beads put together on a modular platform to give generic molecular lipids with tunable properties. The interaction between coarse-grained beads is governed by soft, short-ranged potentials that account for the renormalized hydrophobic effect in implicit solvent. The model faithfully forms micelles, tubular micelles, or bilayers (periodically infinite, bicelles, or vesicles) depending on the packing ratio of the lipid molecules. Importantly, for the self-assembled bilayer membranes it is straightforward to realize gel and fluid phases, the latter of which is most often of primary interest. We show that the emerging physics over a wide range of scales and resolutions is (with some restrictions) universal. The model, when compared to other popular lipid models, demonstrates improvements in the form of increased numerical stability and boosted dynamics. Possible strategies for customizing the models (e.g., adding chemical specificity) are briefly discussed. An implementation is available for the LAMMPS molecular dynamics simulator [Plimpton. J. Comp. Phys. 117, 1-19. (1995)] including illustrative input examples from the simulations we present.

연구 동기 및 목표

• 다양한 척도와 해상도에서 지질 집합체를 시뮬레이션하기 위한 계산적으로 효율적이고 최소한의 모델을 개발하기 위해.
• 기존의 군집화된 지질 모델, 특히 하드 코어 또는 레너드존 잠재력 기반의 모델에 비해 수치적 안정성과 역학을 향상시키기 위해.
• 조절 가능한 상호작용 매개변수를 통해 유체상과 겔상, 밀도 비율, 굽힘 행동과 같은 조절 가능한 막 특성을 구현하기 위해.
• 맞춤형 화학 특성을 가진 막 관련 생물학적 과정을 시뮬레이션하기 위한 모듈식이고 확장 가능한 프레임워크를 제공하기 위해.
• 모델의 일반성을 입증하여 대규모 형태 전이와 막 물리 법칙을 재현할 수 있도록 하기 위해.

제안 방법

• 모델은 지질당 2~5개의 베드를 사용하는 모듈식 입자 기반 프레임워크를 사용하며, 헤드, 인터페이스, 尾부 영역을 각각 다른 상호작용 매개변수로 표현한다.
• 비결합 상호작용은 부드럽고 짧은 범위의 두체 잠재력으로 제어되며, U_pair(r) = - (A/a) sin(r·a) (r ≤ r₀) 및 - (B/b) sin(π/2 + (r - r_c)·c) (r₀ < r < r_c)로 정의되며, r_c에서 부드럽게 0으로 감쇠된다.
• 내분자력은 조임력( K_b(r - r₀)² )과 각도력( K_θ(θ - θ₀)² )을 통해 하모닉 포텐셜로 모델링되어 지질 사슬의 강성 유지를 도모한다.
• 비극성 효과는 인터페이스 및 尾부 베드 간의 응집 상호작용에 암묵적으로 재규합되며, 헤드 베드는 순수하게 반발 상호작용을 갖는다.
• 모델은 LAMMPS에 구현되었으며, 밀레르, 이중층, 고세포 및 시뮬레이션 기반 당김 실험을 통한 검증이 이루어졌다.
• 고세포의 밀도 조절은 안정성을 높이고 플립플롭 리라크세이션의 필요성을 줄이기 위해 일정한 지질당 부피(VPL) 전략을 사용하여 수행되었다.

실험 결과

연구 질문

• RQ1 부드럽고 짧은 범위의 두체 상호작용 잠재력을 갖는 최소한의 군집화된 지질 모델이 밀레르, 이중층, 고세포와 같은 다양한 형태로 자가 조립될 수 있는가?
• RQ2 이 모델의 역학과 수치적 안정성은 SDK나 레너드존 기반 접근법과 같은 기존 모델에 비해 어떻게 비교되는가?
• RQ3 이 모델은 상 행동(유체상 대비 겔상), 굽힘, 형태 전이와 같은 핵심 막 물리 법칙을 어느 정도 재현할 수 있는가?
• RQ4 이 모델의 조절 가능성은 계산 비용을 줄이며 복잡한 생물학적 과정을 시뮬레이션하는 데 활용될 수 있는가?
• RQ5 지질 밀도 전략의 선택(일정한 APL 대비 VPL)은 고세포의 안정성과 평형화에 어떤 영향을 미치는가?

주요 결과

• 모델은 지질 밀도 비율에 따라 외부 제약 없이 밀레르, 튜불러 밀레르, 이중층(주기적, 비셀라형, 고세포형 포함)을 성공적으로 형성한다.
• 유체 이중층 상은 안정적이며 현실적인 역학을 보이며, 지질의 탄성 및 상호작용 매개변수 조절을 통해 상 행동을 조절할 수 있다.
• 준단층 모델(지질당 1.5개의 베드)은 높은 계산 효율성을 보이며, 고세포에 대한 시뮬레이션 기반 당김 실험을 가능하게 하여 막의 굽힘과 기계적 특성을 드러낸다.
• 특히 매우 군집화된 해상도에서 $1/r^{12}$ 또는 레너드존 잠재력 기반의 모델에 비해 향상된 수치적 안정성과 향상된 역학을 보이며, 이는 안정성 향상에 기여한다.
• 일정한 지질당 부피(VPL) 전략은 일정한 단면적당 부피(APL) 전략보다 더 안정적인 고세포를 생성하며, 플립플롭 리라크세이션의 필요성을 최소화하고 구조적 오류를 감소시킨다.
• 부드럽고 짧은 범위의 잠재력은 잘라내기나 이동 없이 부드럽게 力이 감쇠되므로 안정성을 향상시키며 장시간 시뮬레이션을 가능하게 한다.