Viscous vertex model for active epithelial tissues

1. 문제 제기 (Problem)

• 기존 모델의 한계: 기존의 정점 모델 (Vertex Model) 은 주로 세포 간 접합부의 탄성 에너지와 기질 (substrate) 과의 마찰력을 기반으로 조직 역학을 설명합니다. 이 모델에서는 접합부 장력이 탄성 요소로 간주되며, 점성 (viscosity) 기여는 종종 무시되거나 기질 마찰로만 처리됩니다.
• 물리적 현실과의 괴리: 실제 자유 부유형 (free-floating) 배아 조직이나 오르가노이드 (organoids) 와 같은 시스템에서는 기질 마찰력이 내부 점성력 (세포질 및 피질 점성) 에 비해 매우 작거나 무시할 수 있습니다. 특히 세포 피질 (cortex) 의 점성은 세포질보다 수천 배 더 크며, 조직 내 에너지 소산의 주된 원인이 됩니다.
• 수치적 불안정성: 기질 마찰이 0 에 수렴하는 극한에서 기존 모델은 전체 병진 및 회전 운동이 구속되지 않아 운동 방정식 행렬이 특이 (singular) 해지며, 수치적 발산을 일으킵니다.
• 결론: 회전 불변성 (rotational invariance) 을 유지하면서, 기질 마찰이 없는 상태에서도 잘 정의된 (well-posed) 상태가 되며, 세포의 접합부 점성과 벌크 (bulk) 점성을 모두 체계적으로 포함하는 새로운 모델이 필요합니다.

2. 방법론 (Methodology)

저자들은 회전 불변 점성 정점 모델 (Rotationally Invariant Viscous Vertex Model) 을 개발했습니다.

• 점성 소산 메커니즘:
  1. 세포 - 세포 접합부 점성 (Junctional Viscosity, $\eta_s$): 인접한 정점 (vertex) 간의 상대 속도 차이를 접합부 방향에 투영하여 정의합니다. 이는 접합부 길이 변화율 ($\dot{\ell}_{ij}$) 에 비례하는 점성 장력을 생성합니다.
  2. 세포 벌크 점성 (Cell Bulk Viscosity, $\eta_b$): 각 정점과 해당 세포의 기하학적 중심 사이의 상대 운동을 기반으로 정의합니다. 이는 세포 내부의 점성 소산을 모사합니다.
  • 두 점성 모두 회전 불변성을 가지며, 임의의 기준 좌표계에 의존하지 않습니다.
• 운동 방정식 및 라그랑주 승수법:
  • 정점의 힘 평형 방정식은 $C \cdot v = F(t)$ 형태의 행렬식으로 표현됩니다. 여기서 $C$는 기질 마찰, 접합부 점성, 벌크 점성을 모두 포함하는 마찰 - 점성 계수 행렬입니다.
  • 기질 마찰이 0 인 경우 ($\gamma=0$): 행렬 $C$가 특이 행렬이 되는 문제를 해결하기 위해 라그랑주 승수 (Lagrange multipliers) 를 도입하여 전체 운동량과 각운동량 보존 (또는 경계 조건) 을 제약 조건으로 추가합니다. 이를 통해 확장된 안장점 행렬 (extended saddle-point matrix, $C_{ext}$) 을 구성하여 수치적 안정성을 확보합니다.
• 점도 측정 프로토콜:
  • 슬랩 전단 (Slab-shear) 프로토콜: 조직에 전단 변형률 속도를 가하고, 정점 속도와 선 장력 (line tension) 을 기반으로 유효 조직 전단 점도 ($\eta_{tissue}$) 를 추출합니다.
  • 단시간 (Short-time) 및 장시간 (Long-time) 점도: 세포 재배열이 없는 순간적인 응답과 재배열이 포함된 장기 응답을 구분하여 분석합니다.

3. 주요 기여 (Key Contributions)

• 체계적인 점성 정점 모델 정립: 세포의 접합부 점성과 벌크 점성을 동일한 수준에서 처리하고, 회전 불변성을 만족하며 기질 마찰이 0 인 극한에서도 수치적으로 잘 정의된 최초의 모델 중 하나를 제시했습니다.
• 라그랑주 승수 기반 경계 조건 구현: 고정 경계, 지정된 당김 (traction), 자유 슬라이딩 경계 등 다양한 경계 조건을 라그랑주 승수 프레임워크를 통해 통합적으로 구현할 수 있는 방법을 제시했습니다.
• 미시적 점도와 거시적 조직 점도의 연결: 단일 육각형 세포에 대한 해석적 유도 및 무질서한 세포 패킹에 대한 수치 시뮬레이션을 통해, 조직의 유효 점도가 세포의 접합부 점성 ($\eta_s$) 과 벌크 점성 ($\eta_b$) 에 선형적으로 비례함을 증명했습니다.

4. 결과 (Results)

• 점도와 조직 역학:
  • 점도 증가 효과: 접합부 및 벌크 점성이 증가하면 세포 모양이 더 길어지고 (elongated), 공간적으로 상관관계가 있는 형태를 띠게 됩니다.
  • 전단 밴딩 (Shear Banding) 억제: 기질 마찰이 있는 경우 전단 흐름이 경계면 근처에 국소화되는 전단 밴딩 현상이 발생하지만, 점성이 지배적이거나 기질 마찰이 없는 경우 이 현상이 사라지고 뉴턴 유체와 유사한 균일한 흐름이 관찰됩니다.
• 활성 조직 (Active Tissues) 의 거동:
  • 극성 활성 (Polar Activity): 세포의 전방 돌기 (traction) 를 모델링할 때, 점성이 증가하면 난류와 유사한 흐름에서 질서 있는 회전 모드로 전환됩니다.
  • 네마틱 활성 (Nematic Activity): 세포 모양에 의존하는 활성 응력을 모델링할 때, 점성 증가는 위상 결함 (topological defects, $\pm 1/2$) 의 밀도를 줄이고 결함의 조직화를 유도하며, 일관된 흐름 구조를 형성합니다.
  • 결함 역학: 점성이 높을수록 결함의 수가 감소하고, 결함 주변의 응력 및 흐름 패턴의 상관 길이가 길어집니다.

5. 의의 및 중요성 (Significance)

• 자유 부유 조직 및 오르가노이드 연구의 핵심 도구: 기질 마찰이 없는 환경 (예: 배아 발생 초기 단계, 3D 오르가노이드) 에서 조직 역학을 정확하게 모델링할 수 있는 수학적 기반을 제공합니다.
• 세포 기반 모델과 연속체 이론의 연결고리: 이 모델은 세포 수준의 미시적 시뮬레이션과 조직 수준의 연속체 활성 네마틱 (active nematic) 이론 사이의 간극을 메우는 역할을 합니다.
• 실험적 해석 가능성: 향후 생체 상피 조직의 유변학적 측정 데이터를 정량적으로 해석하고, 세포 내부의 점성 소산 메커니즘을 규명하는 데 필수적인 프레임워크를 제공합니다.

이 논문은 활성 물질 물리학 (Active Matter Physics) 분야에서 세포 조직의 복잡한 유변학적 거동을 이해하는 데 있어 점성 효과의 중요성을 부각시키고, 이를 정량적으로 다룰 수 있는 강력한 계산 도구를 제시했다는 점에서 중요한 의의를 가집니다.