세포가 만든 주름, 힘과 형태의 비밀

초록

본 논문은 세포가 연성 기판 위에 가하는 수축성 트랙션 힘을 얇은 막의 푸플-폰 카르만(FvK) 이론과 2차원 활성 네마틱 위상장 모델로 결합해, 주름 형성의 역학을 정량적으로 설명한다. 세포의 수축 강도, 공간 분포, 대칭성에 따라 방사형·원주형·비등방성 주름이 전이되는 메커니즘을 제시하고, 파라미터에 따른 주름 파장·진폭의 스케일링 관계를 도출한다. 실험 데이터와 비교해 모델이 정성·정량적으로 일치함을 확인함으로써, 관찰된 주름 패턴을 세포 힘의 지표로 활용할 수 있는 최소 모델 프레임워크를 제공한다.

상세 분석

이 연구는 세포‑기판 상호작용을 물리적으로 해석하기 위해 두 개의 핵심 이론을 통합한다. 첫 번째는 얇은 탄성 막의 변형을 기술하는 푸플‑폰 카르만(FvK) 방정식이다. 여기서는 얇은 상부 필름(탄성계수 μ_f, 포아송 비 ν_f)이 두꺼운 기판(점성 η_s, 두께 H) 위에 놓인 구조를 가정하고, 평면 응력 σ_ij와 변형 ϵ_ij을 연결한다. 세포가 기판에 가하는 인-플레인 수축 응력 σ_cell_ij는 FvK 식에 외부 하중으로 추가되어, 외부 힘에 의한 굽힘 변위 ω_z와 인-플레인 변위 ω_i의 시간 진화를 제어한다.

두 번째는 세포 자체를 2차원 활성 네마틱 위상장 모델로 묘사한다. 스칼라 필드 φ(0~1)는 세포 내부와 외부를 구분하고, Cahn‑Hilliard 자유에너지로 경계의 두께 λ와 표면 장력 γ를 정의한다. 네마틱 텐서 Q_ij=2(n_i n_j−δ_ij/2)와 Beris‑Edwards 방정식으로 동역학을 기술하며, Navier‑Stokes와 결합해 유동·점성 항을 포함한다. 가장 중요한 활성 응력 σ_active_ij=−ζ(δ_ij+Q_ij)에서 ζ는 이중극성 활성 계수로, 양(수축)·음(팽창) 구분을 가능하게 한다. 특히, Q_ij가 비등방성을 부여해 세포의 장축과 단축에서 서로 다른 트랙션을 생성한다는 점이 주름의 비대칭성을 설명한다.

수치 구현은 OpenFOAM 기반의 유한체적법으로, φ와 Q의 동역학을 고정된 형태(정해진 타원)로 초기화한 뒤, σ_cell_ij를 계산해 FvK 방정식에 입력한다. 시간 스케일 τ_f와 τ_ω를 구분해 안정성을 확보하고, 격자 해상도 Δx=0.5, 시스템 크기 256×256을 사용한다. 파라미터 탐색에서는 μ_f=0.1, ν_s=0.001, η_s=10 등 실험적 PDMS 기판에 근접한 값을 채택했으며, 세포 활동 ζ를 0.005~0.1 범위로 변화시켜 주름 전이를 조사한다.

결과적으로, 세포가 원형이든 타원형이든 상관없이 응력 집중은 세포 극점(긴 축 끝)에서 발생한다. 이는 네마틱 결함(+½)이 두 개 형성되어 극점으로 이동하면서 발생하는 탄성 에너지 최소화 현상과 일치한다. 응력 집중 영역에서 발생한 압축력은 FvK 식에 의해 표면 파동을 유발하고, 파동은 세포 경계에서 최대 진폭을 보이며 멀리까지 전파된다. 파장 λ_w는 기판 전단계수 μ_f와 두께 h_f에 대한 스케일링 λ_w∝h_f(μ_f/μ_s)^{1/3}와 같은 기존 이론과 일치하지만, ζ가 증가하면 λ_w이 감소하고 진폭 A_w가 비선형적으로 증가한다는 새로운 관계를 제시한다.

기판 강성(μ_s) 변화에 따라 전이 현상이 관찰된다. 부드러운 기판에서는 방사형 주름이 주를 이루고, 강성 기판에서는 원주형·비등방성 주름이 우세해진다. 또한, 세포 종횡비 α가 0.2~1 사이에서 증가할수록 주름이 긴 축을 따라 정렬되는 경향이 뚜렷해져, 세포 형태가 주름 방향성을 결정한다는 결론을 얻는다. 이러한 예측은 기존의 Wrinkle‑Force Microscopy 실험 데이터와 정량적으로 일치하며, 실험에서 측정된 파장·진폭 분포를 모델이 재현함을 보여준다.

전반적으로, 이 논문은 (1) 세포 내부 활성 네마틱 스트레스를 물리적으로 정량화, (2) 얇은 막의 비선형 굽힘 이론과 결합, (3) 파라미터에 따른 주름 형태·스케일 전이를 예측하는 통합 모델을 제시함으로써, 세포‑기판 상호작용을 역학적으로 해석하고, 관찰된 주름 패턴을 세포 힘의 비침투적 측정법으로 활용할 수 있는 이론적 기반을 제공한다.