이미지 기반 COMSOL 장기발생 시뮬레이션

초록

본 논문은 마우스 배아의 사지낭 발달을 사례로, 현미경 이미지에서 얻은 3‑차원 해부학적 형상을 MATLAB LiveLink를 통해 COMSOL Multiphysics에 직접 불러와 성장하는 조직 위에 신호전달 모델을 구현하는 방법을 제시한다. ALE(Arbitrary Lagrangian‑Eulerian) 기법을 이용해 급격히 변형되는 도메인에서도 반응‑확산 방정식을 안정적으로 풀 수 있음을 보여준다.

상세 분석

이 연구는 이미지 기반 모델링이라는 최신 흐름을 COMSOL 환경에 성공적으로 접목시킨 점에서 큰 의의가 있다. 먼저 저자들은 형광 및 광학 단층 촬영(OCT)으로 획득한 마우스 배아 사지낭의 3‑D 이미지 데이터를 전처리하여 각각의 조직(외배엽, 중배엽, 내배엽 등)별로 서브도메인을 정의한다. 이 과정에서 MATLAB LiveLink를 활용해 이미지 스택을 voxel‑to‑mesh 변환 파이프라인에 연결함으로써, 복잡한 해부학적 경계를 정확히 재현한 비정형 메쉬를 자동 생성한다. 기존에 COMSOL에서 제공하던 단순 기하학적 모델과 달리, 실제 조직의 비대칭성과 불규칙성을 그대로 반영할 수 있다는 점이 핵심이다.

생성된 메쉬 위에 반응‑확산 방정식과 성장 법칙을 결합한 PDE 시스템을 구현한다. 신호분자(예: FGF, BMP)의 확산·반응을 기술하는 표준 반응‑확산 모델에 더해, 세포 증식·이동에 의해 발생하는 도메인 변형을 ALE 프레임워크로 처리한다. ALE는 메쉬가 변형될 때 물리량을 라그랑지안 좌표계와 오일러안 좌표계 사이에서 자유롭게 매핑하도록 설계돼, 급격한 형태 변화에도 수치적 안정성을 유지한다. 특히, 저자들은 성장률을 이미지 기반으로 추출한 시간‑연속적인 변위 필드와 연동시켜, 실제 배아 성장 속도와 패턴을 모델에 직접 반영하였다.

MATLAB 스크립트를 통해 파라미터 스위핑 및 최적화가 가능하도록 설계했으며, 이는 실험 데이터와 시뮬레이션 결과를 정량적으로 비교할 수 있는 기반을 제공한다. 또한, 서브도메인별 물성(확산계수, 반응속도 등)을 개별적으로 지정함으로써, 조직 특이적 신호전달 차이를 정밀하게 모델링한다. 이러한 접근은 기존에 단일 균일 도메인으로 가정하던 모델링의 한계를 크게 극복한다.

결과적으로, 시뮬레이션은 사지낭의 전방-후방 축에서 FGF와 BMP 신호가 형성하는 농도 구배를 재현했으며, 성장에 따른 도메인 변형이 신호 패턴에 미치는 영향을 정량적으로 분석했다. 특히, 도메인 팽창이 신호전달 경로의 시간적 지연을 초래하고, 이는 조직 패턴 형성에 중요한 피드백 메커니즘으로 작용함을 확인하였다.

이 논문은 이미지 기반 형상 입력, 서브도메인 관리, ALE 기반 변형 해석이라는 세 가지 핵심 기술을 통합함으로써, 복잡한 발달 과정의 수치 모델링을 실험과 긴밀히 연결할 수 있는 새로운 워크플로우를 제시한다. 향후 다중 신호망, 세포 수준의 기계적 상호작용, 그리고 고해상도 시계열 이미지와의 연계까지 확장 가능성이 크다.