세포 극성·운동을 연결하는 기계·생화학 통합 모델

초록

이 논문은 이머스드 바운더리 방법을 이용해, 탄성 막으로 둘러싸인 가변형 세포 경계 안에서 반응‑확산 기반 극성 신호망과 돌출·탄성 힘을 결합한 2차원 컴퓨팅 모델을 제시한다. 단순한 두‑물질 반응‑확산 시스템이 세포 앞쪽에 돌출력을 발생시키고, 탄성 막이 형태를 유지·복원한다. 모델을 통해 순수 기계적·화학적 효과와 완전 결합 시스템을 각각 검증했으며, 자발적 극성 형성과 지속적 이동을 재현함으로써 실제 크롤링 세포의 핵심 행동을 설명한다.

상세 분석

본 연구는 세포극성·운동을 동시에 기술하기 위해 두 가지 핵심 요소를 통합하였다. 첫째, 세포 외곽을 탄성 섬유망으로 모델링하고, 이 섬을 이머스드 바운더리(IB) 프레임워크에 삽입함으로써 변형 가능한 비정상 영역을 자연스럽게 다룰 수 있었다. IB 방법은 라그랑지안 형태의 막을 유클리드 격자에 ‘침투’시켜, 유체·물질 흐름과 구조 변형을 하나의 연산 루프에서 해결한다는 장점이 있다. 여기서는 2D 평면에서 막을 이산화된 스프링 네트워크로 구현했으며, 스프링 상수와 휴식 길이를 통해 탄성 복원력과 전단 저항을 조절했다.

둘째, 내부 생화학적 신호는 전형적인 두‑물질 반응‑확산 시스템(활성형 A와 억제형 B)으로 설정하였다. A는 자기촉진(양의 피드백)과 B에 의한 억제(음의 피드백)를 동시에 갖으며, 확산 계수 D_A<D_B 로 설정해 국소화된 고농도 ‘극성점’이 형성되도록 했다. 반응 항은 FitzHugh‑Nagumo 형태를 차용했으며, 파라미터 조정을 통해 안정적인 전파 파동 대신 정적인 패턴(극성)이나 진동을 유도할 수 있었다.

기계·화학 결합 메커니즘은 A 농도가 높은 영역에 비례하는 돌출력 f_protrude를 정의하고, 이를 막의 법선 방향으로 적용하였다. 이 힘은 탄성 복원력과 균형을 이루어 세포 앞쪽이 팽창하고 뒤쪽이 수축되는 비대칭 형태를 만든다. 시간적 스케일 차이를 고려해, 화학 반응·확산은 초 단위, 막 변형·운동은 분 단위로 통합하였다.

시뮬레이션 결과는 세 가지 실험적 변형을 제시한다. (1) 순수 기계적 경우, 외부 힘 없이도 스프링 네트워크의 초기 불균형이 점진적으로 평탄화되며 원형에 수렴한다. (2) 순수 화학적 경우, 고정된 원형 도메인 안에서 A‑B 패턴이 자발적으로 한쪽에 집중돼 극성점이 형성되고, 이는 경계 조건에 무관하게 유지된다. (3) 완전 결합 경우, 극성점이 형성된 위치에 돌출력이 집중되어 세포가 해당 방향으로 지속적으로 이동한다. 이동 속도와 형태는 스프링 상수, 돌출력 계수, 확산 비율 등에 민감하게 반응한다. 특히, 탄성 복원력이 너무 강하면 돌출력이 억제돼 이동이 멈추고, 반대로 복원력이 약하면 세포가 과도하게 늘어나 비현실적인 형태를 만든다.

이 모델은 기존에 별도로 다루어졌던 기계적 변형과 화학적 극성 형성을 하나의 수치 프레임워크에 통합함으로써, 세포가 내부 신호에 따라 형태를 재구성하고 이동하는 과정을 정량적으로 탐색할 수 있게 한다. 또한, 파라미터 스페이스 탐색을 통해 다양한 세포 유형(예: 래디칼 세포, 섬모 세포)과 환경(경도, 외부 화학 구배)에서 나타나는 행동 양식을 재현할 가능성을 시사한다. 향후 3D 확장, 다중 신호 네트워크, 그리고 실제 실험 데이터와의 정량적 검증이 진행된다면, 이 접근법은 세포 수준의 물리‑생물학적 모델링에 중요한 도구가 될 것이다.