[논문 리뷰] Linking Cellular and Mechanical Processes in Articular Cartilage Lesion Formation: A Mathematical Model

G Kapitanov, Xiayi Wang외 3명|arXiv (Cornell University)|2016. 01. 06.

Osteoarthritis Treatment and Mechanisms참고 문헌 48인용 수 19

한 줄 요약

이 연구는 유한요소 기반 기계적 시뮬레이션과 연령 구조를 가진 세포 및 반응-확산 화학 역학을 통합한 새로운 다스케일 수학 모델을 제시하여 관절 연골에서 외상성 관절염 병변 형성을 시뮬레이션한다. 모델은 초기 기계적 하중이 괴사성 세포 사멸을 유도하며, 이로 인해 ROS, DAMPs, PIC, EPO 등의 염증성 및 항염증성 신호 전달이 조직을 관통해 퍼지는 연쇄 반응을 유도함을 밝혀냈다. 특히 PIC 및 DAMP의 분해율(λF, λM)과 같은 핵심 매개변수들이 14일 간의 병변 진행과 세포 상태 전이에 결정적인 영향을 미친다.

ABSTRACT

A severe application of stress on articular cartilage can initiate a cascade of biochemical reactions that can lead to the development of osteoarthritis. We constructed a multiscale mathematical model of the process with three components: cellular, chemical, and mechanical. The cellular component describes the different chondrocyte states according to the chemicals these cells release. The chemical component models the change in concentrations of those chemicals. The mechanical component contains a simulation of pressure application onto a cartilage explant and the resulting strains that initiate the biochemical processes. The model creates a framework for incorporating explicit mechanics, simulated by finite element analysis, into a theoretical biology framework.

연구 동기 및 목표

• 기계적 손상이 관절 연골의 생화학적 및 세포 반응과 어떻게 연결되는지를 기계적 프레임워크를 개발하는 것.
• 기계적 변형 시뮬레이션과 세포 및 화학 역학의 생물수학 모델을 명시적으로 통합함으로써 기존 모델의 격차를 메우는 것.
• 손상 후 연골세포 사멸과 염증 확산을 이끄는 핵심 임계값과 매개변수 민감도를 규명하는 것.
• 다양한 기계적 및 생화학적 조건에서 치료 효과가 PTOA 발달에 미치는 영향을 평가할 수 있는 예측 도구를 제공하는 것.

제안 방법

• 기계적 구성요소는 Abaqus 유한요소 분석을 사용하여 실린더형 연골 생검 조직에 금속 프로브를 압착함으로써 기계적 변형을 시뮬레이션하고, 변형 분포와 초기 괴사성 세포 밀도를 계산한다.
• 세포 구성요소는 시간에 따라 각 상태에 머무른 시간에 따라 연골세포 상태(CU, CT, CE, ST, SA, DA, DN)를 추적하는 연령 구조 모델을 사용하며, 전이가 화학 신호에 의해 유도된다.
• 화학 구성요소는 ROS, DAMPs, 프로염증세포활성물질(PIC), EPO의 농도를 상태 의존적 생성 및 분해를 고려한 반응-확산 편미분 방정식으로 모델링한다.
• 모델은 기계적 변형을 초기 괴사성 세포 밀도와 연결하며, 이로 인해 DAMPs가 방출되고 염증 및 항염증 반응을 포함한 피드백 루프를 통해 생화학적 연쇄 반응이 시작된다.
• 핵심 매개변수(예: λF, λM, β11)에 대한 민감도 분석을 수행하여 세포 및 화학 역학에 미치는 영향을 평가한다.
• 수치적 해는 유한차분 방법과 암시적 시간 스텝을 사용하여 2차원 축방향-반경 방향 공간에서 손상 후 14일 간의 진행 상황을 시뮬레이션한다.

실험 결과

연구 질문

• RQ1 초기 기계적 하중이 관절 연골에서 괴사성 세포 사멸을 어떻게 유도하고, 그 후 생화학적 신호 전달을 어떻게 유도하는가?
• RQ2 ROS, DAMPs, PIC, EPO와 같은 핵심 신호 분자들이 손상 후 연골세포 사멸과 염증 확산에 어떤 역할을 하는가?
• RQ3 분해율(λF, λM)과 생성율(β11, β13)의 매개변수 변화가 병변 형성 확산에 어떤 영향을 미치는가?
• RQ4 유한요소 기반 변형 데이터 통합이 다스케일 연골 손상 모델의 생물학적 현실성과 예측 능력을 어느 정도 향상시키는가?
• RQ5 염증 연쇄 반응에서 치료적으로 타겟팅할 수 있는 임계값 또는 전환점(티핑 포인트)을 모델이 식별할 수 있는가?

주요 결과

• 모델은 정성적으로 기대되는 행동을 성공적으로 재현하였다: 염증 매개체(PIC, ROS, DAMPs)와 분해성 세포 상태(ST, SA)의 농도 최고치는 압력이 가해진 영역 중심에서 관찰되었으며, 이는 손상 부위의 국소화와 일치한다.
• PIC 분해율 λF를 0.1로 감소시켰을 때 ST, SA 및 ROS 수치가 유의미하게 증가하여 염증세포활성물질의 청소 기능 저하가 병변 진행을 가속화함을 시사한다.
• DAMP 분해율 λM을 0.1로 감소시켰을 경우, 건강한 연골세포(CU)가 거의 완전히 소실되고, 14일째에 염증 화학 물질의 뚜렷한 확산이 관찰되어 DAMP의 지속성이 만성 염증을 이끄는 주요 원인임을 입증한다.
• 민감도 분석은 생물학적으로 타당한 반응을 확인하였다: β11(PIC 생성) 및 β13(EPO 억제)의 변동이 세포 상태 전이 및 사이토카인 농도의 역학에 상당한 영향을 미친다.
• 모델은 EPO 농도가 EPOR 활성 세포(SA)를 건강한 상태(CU)로 되돌리는 데 농도 의존적 효과를 보임을 입증하여, EPO가 내재된 항염증 조절자로서의 역할을 뒷받침한다.
• 유한요소 기반 변형 데이터를 생물수학적 프레임워크에 통합함으로써 생화학적 연쇄 반응의 현실적인 시작이 가능해졌으며, 이는 기계적 자극에 의해 유도된 병리 상태를 모델링할 수 있음을 검증한다.