⚛️ biophysics

A mechanical model of multicellular remodelling in epithelial monolayers

이 논문은 상피 단층의 세포 내 재구성을 시뮬레이션하기 위해 오프-래티스(off-lattice) 셀-센터 모델 내에 동적 참조 프레임(Dynamic Reference Frame, DRF)을 도입하며, 이를 통해 다중 스케일의 기계적 반응을 재현하고 변형 이력이 조직 회복에 미치는 영향을 분석할 수 있게 한다.

원저자: Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

게시일 2026-09-14

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ⚕️ 이것은 동료 심사를 거치지 않은 프리프린트의 AI 생성 설명입니다. 의학적 조언이 아닙니다. 이 내용을 바탕으로 건강 관련 결정을 내리지 마세요. 전체 면책 조항 읽기

우리 몸의 외층과 내부 장기의 내벽은 상피 조직으로 이루어져 있으며, 이는 신체의 보호 장벽이자 문지기 역할을 하는 세포 시트입니다. 이러한 조직은 정적인 상태가 아닙니다. 장기가 성장하거나 치유되거나 환경에 반응함에 따라 이들은 끊임없이 형태를 바꾸고, 늘어나고, 이완합니다. 과학자들은 오랫동안 이 시트들이 당겨지거나 밀릴 때 단순히 한 번에 튕겨 돌아오거나 일정한 속도로 늘어나는 것이 아니라는 사실을 알고 있었습니다. 대신, 이들은 두 가지 뚜렷한 단계로 반응합니다. 즉, 몇 초 안에 일어나는 빠른 탄성적 튕김(elastic snap)과, 몇 분 또는 몇 시간까지 걸릴 수 있는 훨씬 느린 크리프(creep, 서서히 늘어남) 현상입니다. 이 느린 변화는 세포의 내부 기제, 구체적으로는 세포들을 결합하고 있는 단백질 골격과 끈끈한 결합들이 새로운 형태에 적응하기 위해 재구성되는 리모델링(remodeling)의 한 형태입니다. 이러한 이중 속도 행동을 이해하는 것은 조직이 스트레스 하에서 어떻게 무결성을 유지하는지, 그리고 질병 상황에서 어떻게 실패할 수 있는지를 밝혀내는 데 매우 중요하지만, 기존의 컴퓨터 모델들은 지나치게 복잡해지지 않으면서 이 복잡성을 포착하는 데 어려움을 겪어 왔습니다.

연구팀은 이제 이 행동을 시뮬레이션할 수 있는 새로운 방법을 개발하여, 상피 시트의 빠른 튕김과 느린 크리프를 모두 포착하는 디지털 모델을 만들어냈습니다. 각 세포 내부의 모든 미세한 단백질을 모델링하려고 하면 엄청난 컴퓨팅 능력이 필요하기 때문에, 과학자들은 "동적 참조 프레임(dynamic reference frame)"이라는 개념을 도입했습니다. 조직을 세포들이 스프링으로 연결된 집합체라고 상상해 보십시오. 기존 모델에서 이 스프링의 길이는 고정되어 있었습니다. 즉, 조직을 잡아당기면 스프링이 늘어난 후 세포들이 서로 떨어지지 않는 한 그 늘어난 길이를 영원히 유지했습니다. 이 새로운 모델은 세포의 형태에 대한 내부 기억을 나타내는 두 번째의 보이지 않는 스프링 층을 추가합니다. 조직이 늘어날 때 눈에 보이는 세포들은 빠르게 움직이지만, 이 보이지 않는 층은 새로운 현실에 맞추기 위해 자신의 목표 길이를 천천히 변화시킵니다. 이러한 변화는 세포가 느끼는 스트레스에 적응하기 위해 내부 구조를 재조직하는 생물학적 과정을 나타냅니다.

연구진은 이 아이디어를 테스트하기 위해 재료가 힘을 받을 때 어떻게 행동하는지 테스트하는 표준적인 방법인 크리프 및 응력 완화(stress relaxation) 실험의 컴퓨터 시뮬레이션을 실행했습니다. 크리프 실험에서는 디지털 조직 시트의 가장자리에 일정한 힘을 가하고, 연구진은 시간이 지남에 따라 얼마나 늘어나는지를 관찰했습니다. 응력 완화 실험에서는 조직을 고정된 너비로 늘려 놓은 상태에서 그 너비를 유지하기 위해 필요한 힘을 측정했습니다. 모델에 빠른 탄성 반응만을 포함했을 때, 조직은 단순한 고무줄처럼 행동했습니다. 즉, 빠르게 늘어난 후 멈추거나, 힘이 빠르게 떨어졌다가 일정하게 유지되었습니다. 그러나 연구진이 세포의 하부 구조 리모델링을 허용하도록 동적 참조 프레임을 활성화하자, 시뮬레이션은 극적으로 변했습니다. 조직은 처음에는 여전히 빠르게 반응했지만, 그 후에도 시간이 지남에 따라 천천히 계속 늘어났거나 힘이 계속해서 떨어졌으며, 이는 실제 생물학적 실험에서 보이는 이중 단계 행동을 완벽하게 모방했습니다.

이 새로운 모델의 설정을 조정함으로써, 연구팀은 이러한 느린 리모델링의 속도가 세포가 스트레스에 얼마나 민지하게 반응하는지와 얼마나 빨리 내부 "기억"을 업데이트하는지에 달려 있다는 것을 발견했습니다. 만약 세포가 적응하는 속도가 느리다면, 조직은 주로 단순한 탄성 물질처럼 행동한다는 것을 발견했습니다. 하지만 세포가 민감하고 시간이 지남에 따라 내부 구조를 업데이트한다면, 조직은 자연에서 관찰되는 복잡한 다중 속도 반응을 나타냅니다. 결정적으로, 모델은 스트레스의 이력이 중요하다는 것을 보여주었습니다. 만약 조직이 오랜 시간 동안 인장력을 받는다면, 내부 구조는 상당히 적응합니다. 이때 긴장이 해제되면 조직은 즉시 원래대로 돌아가지 않고 천천히 회복되며, 이 회복 속도는 조직이 얼마나 오랫동안 당겨져 있었는지에 따라 달라집니다. 이는 많은 기존 컴퓨터 모델들이 조직이 얼마나 오래 늘려졌는지와 상관없이 회복이 동일할 것이라고 예측했다는 점을 고려할 때 중요한 발견인데, 기존 모델의 결과는 실제 관찰 결과와 모순됩니다.

이 연구는 세포의 형태가 변하는 느린 크리프 변화가 세포의 분열이나 위치 이동이 아니라, 세포 내부 구성 요소의 지속적인 재조직에 의해 발생한다는 것을 확인해 줍니다. 연구진은 세포 내부의 모든 분자를 시뮬레이션할 필요 없이 이러한 행동을 포착할 수 있음을 입증하였으며, 이를 통해 수천 개의 세포를 포함하는 거대한 조직 시트를 연구하는 것이 가능해졌습니다. 그들의 연구는 조직이 과거의 스트레스를 기억하고 적응하는 능력이 그 기계적 설계의 근본적인 특성임을 시사합니다. 세포의 미시적 재조직과 조직의 거시적 행동을 연결하는 프레임을 제공함으로써, 이 새로운 모델은 장기가 어떻게 발달하고, 어떻게 치유되며, 질병의 기계적 힘에 어떻게 반응할 수 있는지 이해하는 데 더 명확한 경로를 제시합니다. 이 발견을 생성하는 데 사용된 코드와 데이터는 다른 과학자들이 사용할 수 있도록 공개되어, 학계가 더욱 상세하고 정확한 살아있는 조직 모델을 구축할 수 있도록 돕고 있습니다.

기술 요약: 상피 단층의 다세포 리모델링에 관한 기계적 모델

문제 제기
많은 포유류 기관의 기초를 형성하는 상피 단층은 외부 힘에 대해 복잡한 기계적 반응을 나타낸다. 실험적 증거(예: Harris 등 [6])는 이러한 조직이 여러 시간 척도에 걸쳐 반응함을 보여준다: 즉, 빠른 탄성 반응(초 단위)에 이어 느린 이완 또는 크리프(분~시간 단위)가 뒤따른다. 기존의 수학적 모델들은 세포 증식이나 재배열을 도입하지 않고는 이러한 거동을 재현하는 데 어려움을 겪는데, 이러한 과정들은 짧은 시간 척도에서는 미미하기 때문이다.

연속체 모델(예: 표준 선형 고체 모델)은 일반적으로 단일 특성 시간 척도를 가정한다.
복합 점탄성 모델(예: 다중 맥스웰 요소)은 여러 시간 척도를 생성할 수 있지만, 종종 현상학적이며 파라미터를 특정 하위 세포 메커니즘과 연결하기 어렵다.
세포 기반 모델(예: 버텍스 역학, 중첩 구체 모델)은 일반적으로 고정된 상호작용 법칙에 의존한다. 증식과 재배열을 무시할 경우, 이들 또한 단일 시간 척도 반응을 보이는 경향이 있다. 하위 세포 요소(Subcellular Element, SE) 모델은 내부 재구성을 통해 다중 시간 척도 거동을 포착할 수 있지만, 높은 계산 비용으로 인해 매우 작은 시스템(<10개 세포)에만 국한되어 조직 규모의 시뮬레이션에는 부적합하다.

방법론
저자들은 특히 보로노이 테셀레이션(Voronoi Tessellation, VT) 프레임워크를 활용하는 오프-래티스(off-lattice) 세포 중심 모델의 새로운 확장을 제안한다. 핵심 혁신은 **동적 참조 프레임(Dynamic Reference Frame, DRF)**의 도입이다.

모델 아키텍처: 조직은 세 가지 세포 중심 집합으로 표현된다:
1. 물리적 세포 ( $r_i$ ): 실제로 변형되는 세포의 위치.
2. 참조 세포 ( $\rho_i$ ): 시간의 흐름에 따라 진화하며 조직의 "리모델링된" 상태를 나타내는 동적 구성.
3. 정적 세포 ( $s_i$ ): 기준 역할을 하는 고정된 초기 구성.
지배 방정식:
- 물리적 운동: 물리적 세포는 점성 항력, 외부 힘, 그리고 물리적 세포와 참조 세포를 연결하는 스프링 힘(식 2) 사이의 힘 균형에 따라 움직인다. 스프링 힘은 물리적 거리( $||r_{ij}||$ )와 참조 거리( $||\rho_{ij}||$ )의 차이에 의존한다.
- 참조 진화: 참조 세포는 별도의 운동 방정식(식 3)에 따라 진화한다. 이들은 두 가지 경쟁하는 힘에 의해 구동된다:
  1. 정적 구성( $s_i$ )으로 끌어당기는 복원력.
  2. 메모리 윈도우 $T_M$ 동안의 시간 평균적 물리적 구성( $\bar{r}_{ij}$ )으로 끌어당기는 리모델링 힘.
무차원화: 시스템은 다음 주요 파라미터들을 정의하기 위해 무차원화된다:
- $\alpha$ (리모델링 비율): 복원력 대비 리모델링 힘의 상대적 강도.
- $\eta$ (리모델링 속도): 물리적 프레임에 상대적인 참조 프레임의 진화 속도.
- $T_M$ (리모델링 메모리): 참조 프레임이 물리적 변형을 평균화하는 시간 창.
시뮬레이션 설정: 모델은 두 가지 실험 프로토콜 하에서 2D 상피 단층(105개 세포)을 시뮬레이션한다:
1. 크리프 이완(Creep Relaxation): 일정한 힘을 가하고 변형률을 기록함.
2. 응력 이완(Stress Relaxation): 변형률을 최대치까지 램프(ramp) 형태로 높인 후 유지하며 응력을 기록함.
분석: 시ما이션 시계열 데이터는 지수 함수의 합으로 피팅된다. 시간 척도의 개수는 단일 지수 함수가 간결한 적합성(오차 $E_1 < 0.05$ )을 제공하는지, 아니면 다중 항이 필요한지에 따라 결정된다.

주요 기여

동적 참조 프레임 (DRF): 본 논문은 세포의 즉각적인 변형과 그 기계적 참조 구성을 분리하는 메커니즘을 도입한다. 이를 통해 저자들은 하위 세포 구성 요소를 명시적으로 모델링하지 않고도 하위 세포 리모델링(세포골격 재구성, 마이오신 회전, 접착 결합 리모델링 등)을 시뮬레이션할 수 있으며, 이를 통해 SE 모델의 높은 계산 비용을 피할 수 있다.
다중 시간 척도의 발현: 저자들은 DRF 프레임워크가 (증식 없이도) 표준 VT 또는 세포 중심 모델에서 결여되었던 특징인 다중 시간 척도 반응(빠른 탄성 변형 + 느린 이완/크리프)을 자연스럽게 생성함을 입증한다.
파라미터 민감도: 연구는 리모델링 파라미터( $\alpha, \eta, T_M$ )가 단일 및 다중 시간 척도 거동 사이의 전이를 어떻게 제어하는지 정량화한다.
이력 의존성(History Dependence): 모델은 변형 해제 시 조직의 회복이 가해진 변형 지속 시간에 따라 달라진다는 실험적 관찰 결과를 성공적으로 재현한다. 이는 표준 탄성 모델이 포착하지 못하는 현상이다.

결과

크리프 실험:
- 리모델링이 없는 경우( $\alpha=0$ ), 조직은 단일 시간 척도 반응(빠른 탄성 변형 후 평탄화)을 보인다.
- 리모델링이 있는 경우( $\alpha > 0$ ), 참조 프레임이 진화함에 따라 추가적인 변형을 허용하는 두 번째의 더 느린 시간 척도가 나타난다.
- 다중 시간 척도의 발현은 가해진 힘의 크기에 관계없이 견고하게 나타나지만, $\alpha$ (리모델링 강도)와 $\eta$ (리모델링 속도)에는 매우 민감하다. 메모리 윈도우( $T_M$ )를 증가시키는 것 또한 다중 시간 척도 거동을 촉진한다.
응력 이완 실험:
- 크리프와 유사하게, 리모델링은 이차적인 이완 단계를 도입한다.
- 램프 시간( $T_R$ )은 제어 파라미터 역할을 한다. 중간 정도의 램프 시간은 뚜렷한 사후 이완 시간 척도를 드러내지만, 매우 짧거나 매우 긴 램프는 이를 가릴 수 있다.
- $T_M$ 을 증가시키면 다중 시간 척도 반응이 관찰되는 파라미터 공간이 확장된다.
변형으로부터의 회복:
- 리모델링이 없는 모델에서, 변형을 해제한 후의 회복률은 유지 시간( $T_H$ )과 무관하다.
- DRF 모델에서 회복률은 $T_H$ 에 반비례한다. 모델은 더 긴 유지 시간이 참조 프레임을 더 많이 변형시켜(축적된 하위 세포 리모델링을 의미), 해제 시 더 낮은 응력과 더 느린 후속 회복을 초래한다고 예측한다. 이 거동은 변형 지속 시간에 따라 회복이 달라진다는 Harris 등 [6]의 일반적인 실험적 발견과 일치하지만, 구체적인 $1/T_H$ 관계는 제안된 모델의 정량적 결과이다.

의의 및 주장
본 논문은 하위 세포 리모델링을 조직 규모에서 연구하기 위한 계산 효율적인 프레임워크를 제공한다고 주장한다. DRF를 통해 표준 세포 중심 모델을 확장함으로써, 저자들은 과도하게 비싼 하위 세포 시뮬레이션이나 측정 불가능한 파라미터를 가진 현상학적 연속체 모델에 의존하지 않고도 복잡한 다중 시간 척도 기계적 거동을 재현할 수 있었다.

저자들은 이 프레임워크가 "생물학적 조직의 발달과 기능을 설명하고 분석하는 데 더 상세한 묘사를 가능하게 한다"고 기술한다. 그들은 이 모델이 모듈식임을 강조한다: DRF 개념은 다른 오프-래티스 모델(예: 버텍스 역학)에 적용될 수 있으며, 중간 참조 상태를 추가함으로써 3D로 확장하거나 세포 성장 및 분화와 같은 과정을 포함할 수 있다. 이 연구는 변형의 이력이 하위 세포 턴오버를 통해 매개되며, 이는 이 접근법을 사용하여 체계적으로 조사할 수 있는 중요한 요소임을 시사한다.

유사한 논문