🧬 biology

Constitutive parameter inference using physics-based data-driven modeling in full volume datasets of intact and torn rotator cuff tendons

본 연구는 전체 부피 데이터를 기반으로 변분 시스템 식별 및 편미분방정식 제약 최적화를 적용하여 건전한 회전근개 힘줄과 파열된 회전근개 힘줄의 구성 매개변수를 추정하였으며, 단순화된 모델이 주요 변형 경향을 포착할 수 있지만 임상적 예측 시뮬레이션을 위한 복잡한 내부 역학을 완전히 재현하려면 현재 구성 방정식들이 여전히 정교화가 필요함을 입증하였다.

원저자: Carla Nathaly Villacís Núñez, Siddhartha Srivastava, Ulrich Scheven, Asheesh Bedi, Krishna Garikipati, Ellen M. Arruda

게시일 2026-05-07

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Carla Nathaly Villacís Núñez, Siddhartha Srivastava, Ulrich Scheven, Asheesh Bedi, Krishna Garikipati, Ellen M. Arruda

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

회전근개 힘줄을 단순히 로프가 아니라, 미세한 섬유로 이루어진 복잡하고 층상 구조를 가진 3 차원 스펀지로 상상해 보십시오. 이 스펀지를 당기면 단순히 곧게 늘어나는 것이 아니라, 특히 찢어진 상태라면 비틀어지고, 전단되며, 복잡한 방식으로 찌그러집니다.

이 논문은 마치 고기술 탐정 이야기와 같습니다. 연구자들은 특수한 MRI 카메라와 매우 정교한 수학을 사용하여 이 스펀지의 거동을 지배하는 '비밀 레시피'(수학적 규칙) 를 찾아내려 했습니다.

다음은 그들의 모험에 대한 요약입니다:

1. 문제: 표면만 보는 것과 전체 그림을 보는 것

일반적으로 과학자들이 힘줄이 어떻게 파열되는지 연구할 때는 자동차의 도장 상태를 보고 엔진이 어떻게 작동하는지 추측하듯, 표면만 관찰합니다. 그들은 카메라를 사용하여 표면이 어떻게 늘어나는지 지켜봅니다. 하지만 저자들은 이것이 땅에 떨어진 나뭇잎만 지켜보며 토네이도를 이해하려는 것과 같다고 깨달았습니다. 내부의 소용돌이치는 바람을 놓치게 됩니다.

그들은 힘줄이 실제로 어떻게 움직이고 찢어지는지 이해하기 위해 깊고 숨겨진 층을 포함한 힘줄의 전체 3 차원 부피를 보고 싶어 했습니다.

2. 실험: '당기고 끊기' MRI

연구진은 인간 어깨와 매우 유사한 양의 어깨를 특수 MRI 기계 안에 넣었습니다.

준비: 그들은 힘줄을 맞춤형 3D 프린팅 홀더에 고정하고 부드럽게 당긴 후, 더 강하게 당겼습니다.
비틀기: 그들은 두 가지 유형의 샘플로 이를 수행했습니다. 하나는 건강한 (완전한) 샘플이고, 다른 하나는 수술적으로 찢어진 (75% 박리된) 샘플이었습니다.
마법 카메라: 단순히 사진을 찍는 대신, MRI 는 힘줄 내부의 모든 미세한 픽셀 (voxel) 을 위한 '모션 트래커'처럼 작동했습니다. 3 차원 블록 내부의 모든 단일 점이 어떻게 움직이고, 비틀리고, 전단되는지 정확하게 매핑했습니다.

3. 탐정 작업: 변분 시스템 식별 (VSI)

이제 그들은 힘줄이 어떻게 움직였는지 보여주는 방대한 양의 데이터를 가지고 있었지만, 그 움직임을 일으킨 '규칙'(수학 공식) 을 알지 못했습니다.

그들은 **변분 시스템 식별 (Variational System Identification, VSI)**이라는 방법을 사용했습니다. 이는 '규칙 맞추기' 게임과 같습니다.

그들은 연조직의 거동을 설명할 수 있는 세 가지 다른 수학적 '레시피'(모델) 라이브러리를 가지고 있었습니다.
- 레시피 A (네오 - 후크안): 매우 단순하고 기본적인 레시피입니다. 복잡한 케이크를 밀가루와 물만으로 설명하려는 것과 같습니다.
- 레시피 B (수정된 HGO): 섬유 방향과 고차 수학을 포함한 특정 항을 포함하는 정교하고 복잡한 레시피입니다.
- 레시피 C (다항식): 항들을 혼합한 중간 수준의 레시피입니다.
그들은 MRI 운동 데이터를 컴퓨터 프로그램에 입력하여 어떤 레시피가 힘줄의 운동을 가장 정확하게 '재연'할 수 있는지 확인했습니다.

4. 주요 발견

수학이 그들에게 알려준 바는 다음과 같습니다:

단순한 레시피의 실패: 기본적인 '네오 - 후크안' 레시피 (레시피 A) 는 너무 단순했습니다. 일반적인 늘어남은 추측할 수 있었지만, 찢어진 힘줄 내부에서 발생한 비틀림과 전단을 완전히 예측하지 못했습니다. 이는 장난감 자동차 모델만 사용하여 실제 자동차 사고를 예측하려는 것과 같아서, 실제 물리학을 놓쳤습니다.
정교함 대 단순함의 놀라움: 연구자들은 가장 복잡한 레시피 (레시피 B) 가 승리할 것으로 예상했습니다. 놀랍게도, 다항식 레시피 (레시피 C) 의 단순화된 버전이 정교한 레시피만큼 잘 수행되었지만 훨씬 적은 성분을 사용했습니다. 일을 처리하기 위해 가장 복잡한 수학이 필요한 것이 아니라, 올바른 성분이 필요하다는 것이 밝혀졌습니다.
'섬유' 성분이 핵심: 그들이 배운 가장 중요한 점은 수학에 반드시 섬유의 방향을 포함해야 한다는 것입니다. 섬유가 어느 방향을 가리키는지 무시하면 모델이 무너집니다. 그들이 힘줄 내부의 미세한 실들을 위한 GPS 와 같은 섬유 방향 지도를 추가했을 때, 모델들은 갑자기 물리적으로 현실적이 되었습니다.
추측 정교화: 첫 번째 추측 (VSI) 은 좋았지만, 그들은 두 번째 단계 (PDE 제약 최적화) 를 사용하여 숫자를 '미세 조정'했습니다. 이는 대략적인 스케치를 선이 완벽해질 때까지 다듬는 것과 같습니다. 이 단계는 컴퓨터가 힘줄이 뻣뻣하고 단단해야 (비압축성) 함에도 불구하고 찌그러진다고 (압축성) 생각했던 주요 문제를 해결했습니다.

5. '아하!' 순간: 보이지 않는 것 보기

가장 흥미로운 부분은 **전단대 (shear bands)**를 보는 것이었습니다.

찢어진 힘줄에서 MRI 는 힘줄의 층들이 찢어짐의 가장자리에서 서로 미끄러지는 것 (박리) 을 보여주었습니다.
단순한 모델은 이를 볼 수 없었습니다.
섬유 방향과 올바른 수학을 포함한 모델들은 이 미끄러짐을 볼 수 있었습니다. 그들은 힘줄이 찢어질 때 발생하는 '내부 마찰'을 성공적으로 재현했는데, 이는 찢어짐이 왜 악화되는지 이해하는 데 중요합니다.

6. 결론

이 논문은 다음과 같이 결론 내립니다:

힘줄을 작은 조각으로 잘라내어 연구할 필요가 없습니다. 전체 부피 MRI 를 사용하여 한 조각으로 전체를 연구할 수 있습니다.
가능한 가장 복잡한 수학이 필요하지 않습니다. 섬유 방향을 포함한다면 더 간단한 모델이 작동합니다.
현재 모델들은 힘줄이 어떻게 찢어지는지 예측하는 데 점점 더 나아지고 있지만, 완벽해지기 위해서는 여전히 작업이 필요합니다.

간단히 말해: 그들은 찢어진 힘줄의 3 차원 디지털 트윈을 구축하고, 그것이 움직이게 만드는 수학 규칙을 찾아냈으며, 힘줄이 어떻게 파열되는지 진정으로 이해하려면 힘줄의 내부를 보는 것이 유일한 방법임을 증명했습니다.

기술적 요약: 회전근개 힘줄의 구성 매개변수 추론

문제 제기
회전근개 힘줄 파열은 어깨 질환의 주요 원인이나, 파열 진행을 예측할 수 있는 역학적 프레임워크의 부재로 인해 부분 두께 파열의 임상적 관리에 대한 논쟁이 지속되고 있습니다. 전통적인 생역학적 특성 분석은 표면 변형률 측정 (예: 디지털 이미지 상관법) 이나 절제된 시료로부터의 응력 - 변형률 곡선에 의존합니다. 이러한 접근법은 파열 성장을 이해하는 데 필수적인 전단 변형률과 박리 메커니즘과 같은 복잡하고 3 차원적인 내부 역학을 포착하지 못합니다. 또한, 기존 유한 요소 모델은 종종 힘줄의 본래적이고 연속적인 역학적 거동을 반영하지 못할 수 있는 단순화된 경계 조건 또는 지역적 분할에 의존합니다. 따라서 완전한 부피의 비파괴 데이터를 활용하여 손상되지 않은 조직과 손상된 조직 상태를 모두 정확하게 나타내는 구성 매개변수를 추론할 수 있는 모델링 프레임워크가 필요합니다.

방법론
저자들은 완전 부피의 자기 공명 영상 (MRI) 데이터를 변분 시스템 식별 (VSI) 과 편미분 방정식 (PDE) 제약 최적화와 통합한 물리 기반 데이터 주도 모델링 프레임워크를 개발했습니다.

실험 데이터 획득: 실험은 인간 상완골 극상근 힘줄의 대리로 검증된 양의 극하근 힘줄에서 수행되었습니다. 시료는 1mm 와 2mm 의 두 가지 변위 수준에서 준일축 인장 하에 완전한 상태와 파열 상태 (75% 활액낭 부착면 박리) 로 테스트되었습니다. 맞춤형 변위 부호화 MRI 시퀀스 (APGSTEi) 가 전체 부피 변위장을 포착하여 체적, 인장, 전단 메커니즘을 포함한 내부 변형률 거동을 밝혔습니다.
수학적 프레임워크: 본 연구는 응력 평형 방정식의 약형 (weak form) 을 활용하는 VSI 를 사용했습니다. 이 접근법은 후보 변형 에너지 밀도 함수 ( $W$ ) 에서 유도된 기저 연산자를 관측된 변위장에 연결합니다. 변형 에너지는 오른쪽 카uchy-그린 텐서의 불변량 ( $I_1, I_2, I_3, I_4$ ) 에 의존하는 체적, 등체적, 이방성 기여로 분해되었습니다.
구성 모델: 매개변수 추론을 위해 세 가지 후보 모델이 평가되었습니다.
- 네오 - 후크 (NH): 단순화된 등방성 모델.
- 수정된 홀차펠 - 가서 - 오그덴 (m-HGO): 섬유 분산과 불변량 $I_1$ 및 $I_2$ 의 고차 항을 포함하는 이방성 모델.
- 다항식 (Pol): 불변량 $I_1, I_2$ 및 섬유 관련 $I_4$ 에 기반한 축소된 다항식 모델.
추론 및 최적화:
- 단계별 회귀를 통한 VSI: 초기 매개변수는 평형 잔차에 기반한 손실 함수를 최소화하여 추론되었습니다. F-검정을 사용한 단계별 회귀 전략이 비활성 변형 메커니즘을 억제하여 변형 에너지 함수의 간결한 표현을 도출했습니다.
- PDE 제약 최적화: 역문제 (특히 거의 비압축성 재료에서 체적 탄성률과 관련하여) 의 부적절성을 해결하고 매개변수를 정제하기 위해 어드저인트 기반 최적화 프레임워크가 적용되었습니다. 이는 평형 제약 조건 하에서 전방 유한 요소 해와 실험적 변위장 간의 오차를 최소화했습니다.
모델 설정: 이 프레임워크는 상완골 두 주위의 힘줄 곡률을 고려하여 고해상도 해부학적 이미지에서 유도된 볼트 단위 섬유 방향 맵을 통합했습니다. 경계 조건은 접촉 시뮬레이션의 복잡성을 피하기 위해 실험적 측정값에서 직접 부여되었습니다.

주요 결과

다중 메커니즘의 활성화: 힘줄의 복잡한 기하학적 구조와 경계 조건으로 인해 단일 준일축 인장 테스트가 체적, 인장, 전단 메커니즘을 동시에 활성화했습니다. 파열된 힘줄은 절단 - 미절단 계면에서 박리를 나타내는 뚜렷한 고전단 변형률 대역을 보여주었습니다.
모델 성능:
- 네오 - 후크 (NH) 모델은 손상된 힘줄의 전단 거동을 포함한 주요 내부 특징을 재현하지 못했으며, 실험적으로 관찰되지 않은 압축 영역과 같은 비현실적인 변형률 분포를 예측했습니다.
- 수정된 HGO (m-HGO) 모델은 합리적인 정확도로 힘줄의 변형 메커니즘을 성공적으로 포착했습니다. 그러나 VSI 는 $I_1$ 과 $I_2$ 의 고차 항이 미미함을 나타내어 축소된 표현이 충분함을 시사했습니다.
- 다항식 (Pol) 모델은 세 항 (네오 - 후크 형태 plus 2 차 섬유 항) 으로 축소되어 m-HGO 모델과 유사한 성능을 발휘했습니다. 이는 단순화된 이방성 공식이 지배적인 역학을 포착할 수 있음을 시사합니다.
매개변수 정제: PDE 제약 최적화는 추론된 매개변수의 물리적 타당성을 크게 향상시켰습니다. VSI 는 종종 체적 탄성률을 하한 경계에서 예측하여 (낮은 식별 가능성을 의미함), 정제된 모델은 거의 비압축성 거동과 일치하는 매개변수를 산출했습니다. 정제는 또한 대부분의 힘줄에서 섬유 분산 ( $\kappa > 0$ ) 을 도입한 반면, VSI 는 종종 완벽한 정렬을 예측했습니다.
오차 감소: 정제된 매개변수를 사용한 전방 예측은 특히 종방향과 박리된 조직의 미끄러짐 거동을 재현하는 데 변위 및 변형률 오차가 감소했습니다. m-HGO 및 Pol 모델은 내부 전단 변형률 국소화를 포착하는 데 NH 모델보다 훨씬 우수했습니다.

의의 및 주장
저자들은 이 작업이 완전 부피 데이터를 사용하여 훈련 및 검증된 최초의 힘줄 전용 모델링 프레임워크이며, 고충실도 추론 체계에 본래의 구조적 특징을 통합했다고 주장합니다.

통합 모델링 접근법: 결과는 회전근개 힘줄이 지역적 분할의 경향에 도전하는 균질한 물성치를 가진 연속 구조로 모델링될 수 있음을 시사합니다. 내부 역학적 이질성과 변형률 경향은 완전 부피 데이터에 기반할 때 전통적인 구성 모델에 의해 효과적으로 포착되었습니다.
간결성: 본 연구는 VSI 가 변형 에너지 밀도 함수의 최소적이지만 충분한 표현을 추출할 수 있음을 보여줍니다. 세 항 다항식 모델은 더 복잡한 m-HGO 공식만큼 효과적이었으며, 예비 조사를 위한 계산 효율적인 대안을 제공했습니다.
임상적 관련성: 조직 절제를 피하고 본래 부착부를 활용함으로써 이 방법론은 in situ에서 물성치를 특성화할 수 있게 합니다. 저자들은 이 프레임워크가 힘줄 손상 및 수복에 대한 개인화된 예측 시뮬레이션의 기초를 마련하여 부분 두께 파열의 치료 선택을 안내할 수 있다고 주장합니다.
한계점: 저자들은 현재 모델이 모든 복잡한 내부 역학을 완전히 재현하지는 않는다고 겸손하게 인정합니다. 한계점에는 최적화 과정의 국소 최소값 가능성, 손상 후 물성 변화에 대한 가정, 그리고 보간된 섬유 방향장에 대한 의존성이 포함됩니다. 향후 연구는 시뮬레이션과 임상 현실 간의 격차를 더 좁히기 위해 위상장 (phase-field) 파괴 모델링과 시간 의존적 거동을 통합할 것을 제안합니다.