Morphogenesis driven by nematic defects in active biological networks

기술 요약: 능동적 생물학적 네트워크 내 네마틱 결함에 의한 형태 발생

문제 제기
본 논문은 세포 형태 발생(cellular morphogenesis), 구체적으로 생물학적 조직이 무질서한 상태에서 특정 형태를 가진 질서 정연한 구조로 어떻게 전이되는지를 모델링하는 근본적인 과제를 다룬다. 세포골격 필라멘트가 부분적인 배향 질서를 가진 네마틱 상(nematic phases)으로 조직될 수 있음은 이미 확립되어 있으나, 이러한 상 내의 위상 수학적 결함(topological defects)이 조직 변형과 패턴 형성을 주도하는 정확한 기계적 역할은 여전히 불분명하다. 기존 모델들은 선형 근사나 미리 정의된 모포젠(morphogen) 기울기에 의존하는 경우가 많으며, 이는 Hydra와 같이 재생되는 유기체에서 관찰되는 네마틱 질서, 탄성, 그리고 스트레스 유도 성장 사이의 복잡하고 비선형적인 상호작용을 포착하는 데 실패할 수 있다. 저자들은 형태 형성 패턴을 사전에 규정하지 않고, 네마틱 질서, 탄성 및 적응형 이완(adaptive relaxation)의 결합으로부터 형태 변화가 발현되는 능동적 생물학적 네트워크(ABN)로서 살아있는 조직을 취급하는 연속체 프레임워크를 개발하고자 한다.

방법론
저자들은 네마틱 상을 형성할 수 있는 교차 결합된 필라멘트 사슬으로 구성된 능동적 생물학적 네트워크로서 살아있는 조직을 모델링하는 준정적(quasi-static) 연속체 프레임워크를 제안한다. 핵심 가정 및 수학적 정식화는 다음과 같다:

• 운동학 및 성장: 조직은 국부적으로 비압축성(incompressible)인 것으로 모델링된다. 형태 발생은 평형 상태의 연속으로 기술된다. 성장은 사전에 규정되지 않으며 국부적인 기계적 스트레스에 의해 유도된다. 변형 구배 $\mathbf{F}$는 네마틱 단계 길이 텐서 $\mathbf{L}$, 탄성 왜곡 $\mathbf{F}_e$, 성장 텐서 $\mathbf{F}_g$, 그리고 이완 텐서 $\mathbf{F}_{rel}$의 곱으로 분해된다.
• 에너지 정식화: 총 자유 에너지 $\mathcal{F}$는 두 가지 주요 기여분을 포함한다:
  1. 네마틱 에너지 ($\mathcal{F}_n$): 에릭센(Ericksen)의 연속체 이론에 기반하며, 디렉터 장(director field) $\mathbf{n}$과 스칼라 질서 매개변수 $s$의 공간적 변화를 설명한다. 결정적으로, $s$는 공간적으로 변화할 수 있도록 허용되며, 위상 수학적 결함에서 0이 되어 특이점(singularity)을 정규화하고 국부적인 등방성 상으로의 전이를 가능하게 한다.
  2. 탄성 에너지 ($\mathcal{F}_e$): 성장과 이완을 포함하도록 확장된 신고전적 변형 에너지 함수이다. 이는 필라멘트 네트워크의 왜곡에 대한 반응을 포착하며 탄성 부분의 변형에 의존한다.
• 결합 메커니즘: 모델은 네 가지 메커니즘을 통합한다: (a) 위상 수학적 결함 근처의 기계적 스트레스, (b) 활동 유도 질서 형성(증가하는 $s$), (c) 스트레스 유도 성장(탄성 에너지 밀도에 비례하여 새로운 물질이 추가됨), (d) 구조적 이완(내부 스트레스를 완화하기 위해 스트레스 없는 구성을 업데이트함).
• 수치 구현: 결합된 비선형 문제는 FEniCSx로 구현된 혼합 유한 요소 정식화를 사용하여 해결된다. 네마틱 디렉터의 머리-꼬리 대칭성과 그에 따른 반정수 결함의 분지 컷(branch cut)을 처리하기 위해, 디렉터 장은 대칭 랭크-2 텐서 장 $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$을 사용하여 재정식화된다. 시뮬레이션은 초기 등방성 구형체로부터 성숙한 형태학적 구조로의 진화를 추적하기 위해 라인 서치(line search)와 증분 하중(incremental loading)을 사용하는 뉴턴 방법(Newton method)을 채택한다.

주요 기여

1. ABN을 위한 연속체 프레임워크: 본 논문은 네마틱 질서, 탄성, 그리고 스트레스 유도 성장/이완을 결합하는 수학적 모델을 도입하여, 선형 쉘 이론을 넘어 비선형 형태 발생 결과를 포착한다.
2. 결함 주도 형태 발생: 연구는 네마틱 결함의 위상 수학적 성질(특히 전하와 공간적 배치)이 형태 진화를 주도하는 주요 동인임을 입증한다. 모델은 결함이 단순한 특이점이 아니라 스트레스를 국부화하고 성장을 직접 유도하는 조직화 중심 역할을 한다는 것을 보여준다.
3. 복잡한 형태학적 구조의 재현: 초기 결함 구성을 변화시킴으로써, 모델은 Hydra 재생에서 관찰되는 세 가지 뚜렷한 형태 발생 결과를 성공적으로 재현한다:
   • 단축(Head-Foot) 재생: 서로 반대 방향에 위치한 두 개의 $+1$ 결함이 신장(elongation)을 유도하고 입(apical rupture)과 발(foot)의 형성을 이끈다.
   • 두두(Two-Head) 발달: 상반구에 위치한 두 개의 $+1$ 결함이 $-1/2$ 결함들에 의해 균형을 이루며, 두 개의 별도 머리 형성을 이끈다.
   • 촉수 형성: 촉수 끝의 $+1$ 결함과 촉수 기저부의 $-1/2$ 결함 구성이 돌출 성장을 유도한다.
4. 스트레스 유도 성장 메커니즘: 시뮬레이션은 성장이 결함 근처의 높은 탄성 스트레스 지역에 국부화됨을 명시적으로 보여주는 한편, 이완은 스트레스가 기하학적으로 해소될 수 있는 영역(입 부위)에서 구멍을 확장시키는 반면, 발이나 촉수 끝에서는 지속적인 스트레스를 유지함을 보여준다.

결과
Hydra 재생에 대한 유한 요소 시뮬레이션은 다음을 밝혀냈다:

• 스트레스 국부화: 네마틱 디렉터 장의 왜곡은 위상 수학적 결함 근처에서 높은 기계적 스트레스를 생성한다. 이러한 스트레스는 새로운 물질의 축적(성장)을 유도한다.
• 비대칭적 결과: 단축 구성에서, 스트레스 이완은 정단부(입)가 파열되고 확장되도록 허용하는 반면, 기저부(발)는 닫힌 상태로 구조적 스트레스를 유지하게 한다.
• 결함의 역할: $+1$ 결함은 돌출 형성을 유도하고 조직 파열을 일으킨다. $-1/2$ 결함은 초기 단계에서 주로 구조적 안정제 또는 "위상 수학적 균형자" 역할을 하지만, 촉수 구성에서는 스트레스 축적이 결국 촉수 기저부에서의 국부적 성장을 촉발한다.
• 형태학적 진화: 모델은 네마믹 질서 매개변수와 그에 따른 스트레스 유도 성장의 진화만을 통해, 구형의 등방성 조직 파편에서 특정 특징(구강 구멍, 촉수)을 가진 복잡하고 신장된 구조로의 전이를 포착한다.

의의 및 주장
저자들은 본 연구가 사전에 지정된 성장 패턴이나 외부 모포젠 기울기 없이도 어떻게 위상 수학적 결함이 능동적 네마틱 네트워크 내에서 복잡한 형태 발생 과정을 주도할 수 있는지에 대한 견고한 이론적 토대를 제공한다고 주장한다. 조직을 네마틱 왜곡에 의해 생성된 기계적 스트레스에 직접 반응하는 능동적 생물학적 네트워크로 취급함으로써, 모델은 미세한 세포골격 조직과 거시적인 조직 형태 사이의 간극을 메운다. 결과는 Hydra의 실험적 관찰과 일치하며, 신체 축 극성의 "기억"과 특정 구조의 출현이 네마틱 결함의 위상 수학적 전하와 배치에 인코딩되어 있음을 시사한다. 본 논문은 이 프레 framework를 순수 생화학적 신호 전달 모델과는 구별되는, 살아있는 조직의 패턴 형성을 이해하기 위한 내재적 기계적 동인으로서의 단계로 위치시킨다.