Solid-to-solid transition in dense assemblies of elongated cells

기술 요약: 조밀한 신장된 세포 집합체에서의 고체-고체 전이

문제 제기
합류된 생물학적 조직에서 세포 형태는 등방성 상피 세포 형태부터 고도로 신장된 내피 세포 형태까지 다양하게 나타난다. 표준 버텍스 모델(vertex models)은 "목표 형상 지수(target shape index)" $p_0$(세포의 선호 둘레와 선호 면적의 제곱근 사이의 비율)를 통해 조직의 강성을 제어함으로써 조직 역학을 설명한다. 이러한 표준 모델에서는 $p_0$가 증가함에 따라 이웃 교환(neighbor exchange)에 대한 에너지 장벽이 낮아지며, 결국 조직의 강성이 사라지는 고체-유체 전이(solid-to-fluid transition)를 유도한다. 그러나 많은 생물학적 세포 유형(예: 근아세포, 신경 전구세포, 내피 세포)은 면적-둘레 경쟁보다는 세포 골격 조직이나 극성에 의해 구동되는 고유하고 수동적인 신장된 형태에 대한 선호도를 가진다. 본 연구가 다루는 핵심 질문은, 이러한 자율적인 세포 신장이 예측된 대로 형상 지수에 의한 유체화를 일으키는지, 아니면 강성으로 가는 별개의 역학적 경로를 구축하는지 여부이다.

방법론
저자들은 "세포 신장 탄성(Cell Elongation Elasticity, CEE)" 항을 도입하여 표준적인 합류 세포 시트의 버텍스 모델을 확장하였다.

• 모델 프레임워크: 조직은 꼭짓점 위치 $\mathbf{r}_i$를 자유도로 갖는 다각형 타일링으로 모델링된다. 전체 역학 에너지 $E$는 표준 버텍스 모델 에너지($E_{svm}$)와 CEE 에너지($E_{cee}$)의 합이다.
• 에너지 항:
  • $E_{svm}$은 목표 면적 $A_0$ 및 둘레 $P_0$로부터의 편차에 대해 페널티를 부여한다.
  • $E_{cee} = \sum_{\alpha} \frac{1}{2}K_Q [\text{tr}(\mathbf{Q}_\alpha^2) - \frac{1}{2}q_0^2]^2$는 선호하는 비등방성 형상 $q_0$로부터의 편차에 대해 페널티를 부여한다. 여기서 $K_Q$는 신장 강성이고, $\mathbf{Q}_\alpha$는 세포 변(cell edge)에 기반한 2차 텐서로 세포 형상의 비등방성을 기술한다.
• 역학 및 시뮬레이션: 시스템은 마찰력과 형상 완화력(shape-relaxation forces) 사이의 균형을 맞추는 과감쇠 역학(overdamped dynamics)을 통해 진화하며, 힘은 총 에너지로부터 유도된다. 시뮬레이션은 주기적 경계 조건을 가진 정규 육각형 격자(및 무작위 보로노이 격자에서 테스트됨) 위에 $N_c = 10^3$개의 세포를 초기화하여 수행되었다. 파라미터들은 $\sqrt{A_0}$ 및 응력 척도에 의해 재조정되었다.
• 분석: 저자들은 아핀 변환(affine transformations) 하에서 역학적 에너지, 전단 탄성률($G_{xy}$), 그리고 항복 응력($\sigma_{yield}$)에 대한 해석적 표현식을 유도한다. 이들은 이론적 예측을 다양한 $K_Q$와 $P_0$에 대해 수치 시뮬레이션과 비교한다.

주요 결과

1. 고체-고체 전이: 신장이 조직을 유체화할 것이라는 표준적인 예상과는 반대로, CEE의 도입은 고체-고체 전이를 유도한다. 신장 강성 $K_Q$가 증가함에 따라, 조직은 질서 정연한 등방성 고체 상태에서 무질서한 비등방성 고체 상태로 전이된다.
2. 강성의 지속성: 비등방성 고체 상태(높은 $K_Q$)에서 조직은 유한한 전단 강성과 유한한 항복 응력을 유지한다. 장기 전단 탄성률 $G_{xy}^{(LT)}$는 임계 강성에서 0으로 감소한 후 다시 상승하여 고원(plateau)을 형성하는데, 이는 강성의 회복을 나타낸다.
3. 경화 메커니즘:
   • 작은 $K_Q$ (등방성 영역)에서, CEE는 전단 탄성률에 음의 기여를 한다.
   • 큰 $K_Q$ (비등방성 영역)에서, 세포 신장은 세포 간 장력($T_{\alpha\beta}$)의 증가를 유도한다. 이 장력은 세포 둘레 강성 $K_P$에 의해 지배되는 선형 변형 경화(strain-stiffening) 영역을 이끌어내어, 양(+)의 유한한 전단 탄성률을 복원한다.
4. 항복 응력 거동: 항복 응력은 단기 전단 탄성률이 사라지는 지점을 넘어서도 유한하게 유지된다. 시스템은 $K_Q$가 더 증가함에 따라 취성 항복(brittle yielding, 급격한 응력 강하)에서 연성 항복(ductile yielding, 부드러운 전이)으로의 전이를 보인다.
5. 상태도:
   • $P_0 < P_0^*$ (표준 고체 영역)의 경우, $K_Q$를 증가시키면 고체-고체 전이가 발생한다.
   • $P_0 > P_0^*$ (표준 유체 영역)의 경우, 미소한 $K_Q$만으로도 조직을 고체화할 수 있으며, 임계점이 $K_Q \to 0^+$로 이동한다.
6. 견고성: 이러한 결과는 정규 육각형 및 무작위 보로노이 초기 패턴 모두에서 성립하며, 선호하는 신장량 $q_0$의 변화에도 견고하다.

의의 및 주장
본 논문은 세포 형상과 조직의 유변학(rheology)을 분리함으로써, 신장된 세포 형태가 반드시 유체와 같은 역학을 의미한다는 추론에 도전한다.

• 새로운 제어 파라미터: 저자들은 목표 형상 지수와 구별되는, 조직 역학을 위한 독립적인 제어 파라미터로서 CEE를 확립하였다.
• 미시적 메커니즘: 본 연구는 세포가 어떻게 자율적으로 자신의 형상을 설정하는지와 조직이 변형에 저항하는 방식을 연결함으로써, 생물학적 조직에서의 고체-고체 전이에 대한 미시적 메커니즘을 제공한다.
• 능동 모델과의 대비: 자율적인 신장이 자발적인 흐름과 유체화를 유발하는 능동적 극성 또는 네마틱 모델과 달리, 이 수동적 역학 평형 모델은 자율적 신장이 조직의 강성을 유지하거나 복원할 수 있음을 보여준다.
• 형태 형성의 함의: 결과는 신장된 세포들로 구성된 조직(예: C2C12 근아세포 또는 HUVEC)이 상당한 기계적 저항력을 가질 수 있음을 시사하며, 이는 배아 접힘(embryo folding)이나 상처 폐쇄와 같은 형태 형성 과정을 모델링할 때 반드시 고려되어야 한다.

저자들은 강하게 신장된 세포 배양물의 강성을 측정하는 것이 CEE와 조직 유변학 사이의 연결 고리를 직접적으로 테스트하는 방법이 될 것이라고 결론지으며, 세포 형태만을 토대로 조직의 유동성을 추론하는 것에 대해 경고한다.