염증 반응과 과립 조직 특성이 인간 골절 치유에 미치는 영향: 계산 모델링 연구

초록

본 연구는 유한요소(FE) 기반의 기계생물학 모델을 이용해 골절 초기 단계인 염증 반응·과립 조직 형성·초기 골막(콜러스) 형성이 치유에 미치는 영향을 정량화하였다. MSC 확산계수, 과립 조직 탄성계수, 콜러스 두께, 그리고 골절 간극 크기를 변수로 설정하고 파라메트릭 시뮬레이션을 수행하였다. MSC 이동이 빠르고, 과립 조직이 더 단단하며, 콜러스가 두껍고 간극이 작을수록 치유가 촉진되지만, 일정 수준을 넘으면 포화 현상이 나타난다. 또한, 골절면이 매우 가깝다면 콜러스 형성이 필요 없으며, 이는 1차·2차 치유 개념과 일치한다.

상세 분석

이 논문은 기존 기계생물학 모델이 간과해 온 골절 초기 단계, 즉 염증 반응과 과립 조직(그라뉼레이션 티슈) 형성을 정밀하게 재현하려는 시도이다. 저자들은 Lacroix & Prendergast(2002)의 모델을 기반으로, 4노드 사각형 변위·공극압 요소를 사용해 원통형 골절 모델을 구축하였다. 주요 변수는 네 가지로 설정했는데, (1) MSC(중간엽줄기세포) 확산계수(D), (2) 과립 조직의 탄성계수(Eg), (3) 초기 콜러스 두께(d), (4) 절단면 절반 간극(h)이다. 각 변수는 실험적 보고와 물리적 한계를 고려해 7단계씩 변동시켰으며, 특히 D는 0.001100 mm²/일, Eg는 0.012 MPa, d는 18 mm, h는 0.54 mm 범위로 설정했다.

시뮬레이션은 하루 단위로 120일까지 진행했으며, 각 단계에서 옥타헥사드 전단 변형률과 유체 흐름을 계산해 조직 분화 규칙(Prendergast et al., 1997; Sapotnick & Nackenhorst, 2015)을 적용했다. 조직이 섬유, 연골, 골로 전환되는 임계값을 정하고, ‘연골 콜러스(CC)’, ‘골연결(BB)’, ‘완전 골콜러스(BC)’가 형성되는 시점을 치유 지표로 삼았다.

결과는 네 변수 모두 치유 속도에 비선형적인 영향을 미침을 보여준다. MSC 확산계수가 0.11 mm²/일 구간에서는 정상적인 치유가 이루어졌으며, 0.01 mm²/일 이하에서는 세포 밀도가 120일 내에 50 %에 도달하지 못해 비연합이 발생했다. 확산계수가 10 mm²/일 이상으로 증가해도 치유 속도는 크게 변하지 않아 포화 현상이 확인되었다. 과립 조직의 탄성계수는 0.2 MPa 이하에서는 치유 시점에 차이가 미미했지만, 0.22 MPa 구간에서는 연골 콜러스가 평균 10일, 골연결이 16일, 완전 골콜러스가 24일 빨라졌다. 콜러스 두께는 1 mm일 때 섬유성 콜러스와 비연합을 초래했으며, 2 mm에서는 4개월 내에 완전 골콜러스를 형성했다. 36 mm 범위에서는 치유가 최적화되어 610주 내에 완전 골콜러스가 형성되었고, 7~8 mm에서는 1개월 이내에 치유가 완료되었다. 간극(h)은 작을수록(0.5 mm) 콜러스 형성이 불필요해 1차 치유와 일치했으며, 간극이 커질수록 연골 및 골 형성에 필요한 콜러스 두께가 증가했다.

이러한 파라미터 간 상호작용은 ‘포화’와 ‘임계값’ 개념을 제시한다. 즉, MSC 이동 속도와 과립 조직 강성, 콜러스 두께가 일정 수준을 초과하면 추가적인 개선 효과가 없으며, 오히려 과도한 조직 형성은 비효율적인 기계적 환경을 초래할 수 있다. 연구는 또한 초기 콜러스가 ‘이상적인’ 두께와 강성을 가져야 최적의 치유 시간을 달성한다는 가설을 실험적으로 뒷받침한다.

임상적 함의로는, 초기 염증 억제제 투여가 MSC 이동을 저해하거나 과립 조직 형성을 방해하면 치유 지연이 발생한다는 점을 강조한다. 또한, 수술적 고정(플레이트·스크류 등)에서 간극을 최소화하면 콜러스 형성을 억제하고 1차 치유를 유도할 수 있다. 모델은 파라미터 조정을 통해 약물·기계적 개입의 최적 시점을 예측할 수 있어, 개인 맞춤형 치료 설계에 활용 가능하다. 다만, 모델이 선형 탄성·다공성 가정에 머물러 있고, 혈관신생·면역세포 동역학을 포함하지 않은 점은 향후 보완이 필요하다.