ECMSim: A high-performance interactive web application for real-time spatiotemporal simulation of cardiac ECM signaling and diffusion

ECMSim 기술 요약: 심장 ECM 신호 전달 및 확산의 실시간 시공간적 시뮬레이션을 위한 고성능 대화형 웹 애플리케이션

문제 제기
세포외기질 (ECM) 리모델링은 심혈관 병리생리학, 특히 심장 섬유화와 심부전에서 중심적인 역할을 하는 계산 집약적이고 기작적으로 복잡한 과정입니다. ECM 역학을 정확하게 예측하려면 다중 스케일 반응 - 확산 역학, 기계화학적 피드백, 그리고 공간적으로 이질적인 세포 반응을 동시에 통합해야 합니다. 기존 계산 모델은 다음과 같은 중대한 한계에 직면해 있습니다: 많은 모델이 세포 간 상호작용과 조직 조직화를 모델링하는 데 필요한 공간 해상도가 부족합니다; 종종 특수한 소프트웨어 설치와 전문 지식을 요구하여 접근성을 제한합니다; 그리고 일반적으로 매개변수 탐색을 위한 실시간 상호작용이 부족합니다. furthermore, 전통적인 실험적 접근법은 공간적 분자 확산 및 피드백 조절과 결합된 다중 경로 신호 전달 역학의 전체 복잡성을 포착하는 데 어려움을 겪습니다.

방법론
저자들은 심장 흉터 조직에서의 이질적 기질 리모델링을 시뮬레이션하도록 설계된 웹 기반 플랫폼인 ECMSim을 소개합니다. 이 시스템은 이전에 검증된 단일 세포 논리 기반 ODE 모델 (Zeigler 등) 을 공간적으로 분해된 프레임워크로 확장한, 심장 섬유아세포 신호 전달 네트워크에 대한 포괄적인 시스템 수준 모델을 기반으로 구축되었습니다.

• 수학적 프레임워크: ECMSim 의 핵심은 137 만 개의 결합된 상미분방정식 (ODE) 시스템으로 구성됩니다. 이 시스템은 100 × 100 의 공간 격자 (10,000 개 세포) 를 모델링하며, 각 세포는 137 개의 분자 종 (115 개의 세포 내 신호 전달 종과 17 개의 세포 외 ECM 종) 을 포함합니다. 이 네트워크는 10 개의 입력 분자 (예: TGF-β, AngII, IL6, 장력), 막 수용체, 세포 내 캐스케이드 (MAPK, PI3K-Akt-mTOR, Rho GTPase), 전사 인자, 그리고 ECM 출력물 (콜라겐, MMP, TIMP) 을 포함하는 259 개의 조절 상호작용 (엣지) 을 포함합니다.
• 공간 역학: 세포 간 통신은 5 가지 주요 피드백 신호 분자 (TGF-βfb, AngIIfb, IL6fb, ET1fb, tensionfb) 의 확산과 ECM 분자의 제한된 확산을 통해 매개됩니다. 공간적 진화는 세포 격자에서 유한 차분 근사를 사용한 이산 확산 방정식으로 모델링됩니다. 유한 크기 효과를 최소화하기 위해 주기적 경계 조건이 구현됩니다. ECM 분자는 분자 크기가 크고 국소 침착 경향이 있음을 반영하여 피드백 분자의 확산 속도의 20% 속도로 확산됩니다.
• 수치적 적분: 시스템은 적응형 시간 단계와 함께 오일러 전진법으로 해결됩니다. 수치적 안정성과 생리학적 관련성을 보장하기 위해 분자 농도는 정규화된 범위 [0, 1] 내에서 클램핑됩니다.
• 계산 구현: 시뮬레이션 엔진은 최적화된 C++ 로 작성되었으며 Emscripten 툴체인을 사용하여 WebAssembly로 컴파일됩니다. 이를 통해 애플리케이션은 표준 웹 브라우저에서 직접 실행될 수 있으며, 성능은 네이티브 C++ 코드에 근접합니다 (네이티브 속도의 약 20% 이내). 아키텍처는 객체 지향 설계, 미리 할당된 메모리 전략, 그리고 캐시 지역성 최적화를 활용하여 막대한 계산 부하를 처리합니다.
• 사용자 인터페이스: 플랫폼은 조정 가능한 입력 분자 농도를 적용할 세포 부분집합을 선택하기 위한 "브러시 스타일" 도구가 있는 대화형 인터페이스를 제공합니다. 여기에는 9 가지 주요 속도 상수 (예: 입력 처리, 피드백 민감도, 분해, 확산 속도) 의 실시간 조정을 위한 매개변수 슬라이더와 여러 개의 결합된 시각화 (히트맵, 선 그래프, 벡터장, 3D 표면) 가 포함됩니다.

주요 기여

1. 확장 가능한 공간 모델링: ECMSim 은 심장 섬유아세포 모델링을 단일 세포 표현에서 반응 - 확산 결합을 갖춘 공간적으로 분해된 10,000 개 세포 격자로 확장하여 조직 수준의 조직화와 창발적 공간 패턴 연구를 가능하게 합니다.
2. 실시간 상호작용: WebAssembly 를 활용함으로써 플랫폼은 137 만 개의 변수를 포함하는 시뮬레이션에 대한 실시간 성능을 달성하여, 사용자가 고성능 컴퓨팅 클러스터의 전형적인 지연 시간 없이 매개변수 공간을 대화식으로 탐색하고 동적 반응을 관찰할 수 있게 합니다.
3. 접근성과 민주화: 웹 기반 배포는 특수한 소프트웨어 설치 필요성을 제거하여 다양한 플랫폼의 더 넓은 범위의 연구자와 교육자에게 정교한 시스템 생물학 모델링을 접근 가능하게 만듭니다.
4. 포괄적인 네트워크 통합: 이 모델은 자가분비 및 파라크린 피드백 루프를 고려한 통합 프레임워크에 심장 섬유화와 관련된 주요 신호 전달 경로 (TGF-β/Smad, MAPK, PI3K-Akt-mTOR, 기계적 전달) 를 통합합니다.

결과

• 성능 벤치마크: MacBook Pro (Apple M1 Max) 에서 ECMSim 은 ODE 속도 평가, 오일러 적분 업데이트, 그리고 확산 이웃 계산을 포함하여 시간 단계당 약 484 만 개의 연산을 실행합니다. 시스템은 약 523 MB 의 총 메모리 사용량을 유지하면서 실시간 상호작용을 유지하며, CPU 시간의 99% 이상이 ODE 솔버에 할당됩니다.
• 시공간 역학: 시뮬레이션은 섬유화 반응의 전파를 모델링하는 플랫폼의 능력을 보여줍니다. 국소적 자극 (예: TGF-β를 사용한 X 자형 브러시 선택) 은 ECM 생산으로 이어지는 세포 내 캐스케이드를 유발하며, 이는 이웃 세포로 확산되어 공간적으로 이질적인 농도 구배를 생성합니다. 시스템은 신호 활성화와 ECM 축적 사이의 지연 시간뿐만 아니라 (예: TGF-β + AngII + IL6) 조합 사이토카인 자극의 시너지 효과를 포착합니다.
• 매개변수 민감도: 민감도 분석은 ECM 생산률 ($k_{production}$) 이 ECM 리모델링 결과의 지배적 조절자임을 밝히며, 그 다음으로 활성화 및 분해 속도가 따릅니다. 이 모델은 다른 매개변수에서의 중간 정도 변이에 대해 견고하며, 이는 활성화된 섬유아세포에서 생산률이 기질 축적을 주로 주도한다는 생물학적 이해와 일치합니다.

의의 및 주장
저자들은 ECMSim 이 복잡한 생물학적 시스템을 위한 접근 가능한 계산 모델링에서 중요한 진전을 나타낸다고 주장합니다. 수학적 엄밀함과 사용자 친화적인 웹 기술을 결합함으로써, 이 플랫폼은 현재 심장 섬유화 연구의 중대한 격차, 특히 공간 해상도와 실시간 상호작용의 부재를 해결합니다.

이 논문은 ECMSim 이 전통적 실험의 시간적 및 자원적 제약 없이 병리적 조건, 가설적 시나리오, 그리고 실험적 약물이 표적 수용체에 미치는 영향을 조사하는 데 기여한다고 제시합니다. 이는 가설을 테스트하고 모델을 검증하기 위한 "디지털 실험실" 역할을 합니다. 주로 심장 섬유화에 초점을 맞추고 있지만, 저자들은 근본적인 반응 - 확산 프레임워크가 암 미세환경, 상처 치유, 그리고 재생 생물학을 포함한 다른 응용 분야를 모델링하는 기초를 제공한다고 지적합니다. 이 플랫폼은 분자 수준의 실험적 관찰과 조직 수준의 병리생리학 사이의 간극을 해소하여 학제 간 협력을 장려하고 새로운 치료 전략 탐색을 가속화하는 것을 목표로 합니다.