[논문 리뷰] Coupling particle-based reaction-diffusion simulations with reservoirs mediated by reaction-diffusion PDEs

Margarita Kostré, Christof Schütte외 2명|arXiv (Cornell University)|2020. 05. 29.

Block Copolymer Self-Assembly참고 문헌 53인용 수 14

한 줄 요약

이 논문은 시간에 따라 변화하고 공간적으로 비균일한 저장소를 나타내는 입자 기반 반응-확산(PBRD) 시뮬레이션과 매크로스코픽 반응-확산 PDE 간의 수학적으로 일관된 결합 방법을 제시한다. 개방계에 대한 평균장 이론과 이차 반응을 유도함으로써, 입자 수가 보존되지 않는 개방 생화학계를 정확하고 확장 가능한 방식으로 시뮬레이션할 수 있으며, 스케일 간에 적절한 유량 교환을 보장한다.

ABSTRACT

Open biochemical systems of interacting molecules are ubiquitous in life-related processes. However, established computational methodologies, like molecular dynamics, are still mostly constrained to closed systems and timescales too small to be relevant for life processes. Alternatively, particle-based reaction-diffusion models are currently the most accurate and computationally feasible approach at these scales. Their efficiency lies in modeling entire molecules as particles that can diffuse and interact with each other. In this work, we develop modeling and numerical schemes for particle-based reaction-diffusion in an open setting, where the reservoirs are mediated by reaction-diffusion PDEs. We derive two important theoretical results. The first one is the mean-field for open systems of diffusing particles; the second one is the mean-field for a particle-based reaction-diffusion system with second-order reactions. We employ these two results to develop a numerical scheme that consistently couples particle-based reaction-diffusion processes with reaction-diffusion PDEs. This allows modeling open biochemical systems in contact with reservoirs that are time-dependent and spatially inhomogeneous, as in many relevant real-world applications.

연구 동기 및 목표

• 입자 기반 반응-확산(PBRD) 시뮬레이션과 저장소를 나타내는 매크로스코픽 반응-확산 PDE 간의 일관된 수학적 프레임워크를 개발하기 위해.
• 입자 수가 보존되지 않는 개방계를 모델링하는 데 도전하는 바, 그로 인해 그랜드 캐노니컬 집단 접근이 필요함을 고려하기 위해.
• 특히 이차 반응에 대해 미세구조 반응 속도와 매크로스코픽 PDE 매개변수 간의 정확한 관계를 수립하기 위해.
• PBRD 시뮬레이션과 PDE 저장소 간의 수치적 일관성을 확보하여, 복잡한 생화학계의 정확한 다스케일 모델링을 가능하게 하기 위해.

제안 방법

• 확산 입자로 이루어진 개방계에 대한 평균장 이론을 유도하여 저장소와의 입자 교환을 모델링하기 위해.
• 이차 반응을 포함한 PBRD 시스템에 대한 평균장 수식을 개발하여, 미세구조 반응 속도 상수와 매크로스코픽 속도 상수 간의 관계를 설정하기 위해.
• 유도된 평균장 관계를 바탕으로 PBRD 시뮬레이션과 PDE 저장소 간의 일관된 결합 기법을 구성하기 위해.
• PDE를 해결하기 위해 크랭크-니콜슨 유한차분 방법을 구현하여 공간 및 시간에 대해 두 번째 차수의 정확도를 확보하기 위해.
• PDE의 반응 및 확산 항을 처리하기 위해 스트랑 분할(Strang splitting)을 적용하여 수치적 안정성과 정확도를 유지하기 위해.
• PDE 저장소가 PBRD 시뮬레이션에서 입자 유입 및 유출를 올바르게 매개함을 확인함으로써 방법을 검증하여 열역학적 일관성을 유지하기 위해.

실험 결과

연구 질문

• RQ1 시간 및 공간에 따라 변화하는 저장소를 나타내는 매크로스코픽 반응-확산 PDE와 입자 기반 반응-확산 시뮬레이션을 어떻게 일관되게 결합할 수 있는가?
• RQ2 특히 이차 반응에 대해 PBRD의 미세구조 반응 속도와 PDE의 매크로스코픽 속도 상수 사이의 정확한 평균장 관계는 무엇인가?
• RQ3 그랜드 캐노니컬 프레임워크 내에서 개방 PBRD 시스템의 입자 수 변동성을 어떻게 적절히 모델링할 수 있는가?
• RQ4 PBRD 시뮬레이션과 PDE 저장소를 다단계 스케일에서 결합할 때 일관성과 안정성을 보장하는 수치적 기법은 무엇인가?
• RQ5 공간적으로 비균일하고 시간에 따라 변화하는 저장소를 포함할 경우, PBRD 시뮬레이션의 정확도와 동역학에 어떤 영향을 미치는가?

주요 결과

• 저자들은 확산 입자로 이루어진 개방계에 대한 평균장 이론을 유도하여, 그랜드 캐노니컬 집단에서 저장소와의 입자 교환을 일관되게 모델링할 수 있도록 하였다.
• PBRD의 미세구조 이차 반응 속도와 매크로스코픽 PDE 속도 상수 간의 정확한 평균장 관계를 수립하였으며, 이는 확산 계수가 속도 관계에 포함되지 않는다는 핵심 수정 사항을 보여주었다 — 이는 이전의 가정에 대한 중요한 수정이다.
• 제안된 결합 기법은 PBRD 시뮬레이션과 PDE 저장소 간의 수학적 일관성을 확보하여, 공간적·시간적 변화가 있는 경계 조건을 갖는 개방 생화학계의 정확한 시뮬레이션을 가능하게 하였다.
• 수치적 검증을 통해 크랭크-니콜슨 방법과 스트랑 분할이 공간 및 시간에 대해 두 번째 차수의 정확도를 제공함을 확인하였으며, 안정적이고 신뢰할 수 있는 PDE 해를 보장하였다.
• 이 방법은 빈도가 공간과 시간에 따라 변화하는 세포 외부 환경을 갖는 세포 신호 전달과 같은 복잡하고 현실적인 생물학적 시나리오를 시뮬레이션할 수 있게 하였다.
• 구현 코드는 GitHub에 공개되어 있으며 MIT 라이선스 하에 제공되어 재현 가능성을 높이고 다스케일 모델링 워크플로우에의 통합을 촉진한다.