Enucleated incompressible red blood cells in shear flow: theoretical analysis of shape instabilities

이 논문은 과잉 막 면적을 고려하기 위해 기존 모델을 확장하여, 형상 불안정성과 그로 인한 스토마토사이트(stomatocyte) 및 트리로브(trilobe) 형태가 막 장력, 재료 탄성 계수 및 점도 비에 어떻게 결정적으로 의존하는지를 밝히는 전핵된 비압축성 적혈세포의 전단 유동에 관한 섭동 이론적 프레임워크를 제시한다.

핵심

개요: 모양이 변하는 적혈구

적혈구(RBC)를 딱딱한 동전이 아니라, 신축성 있고 탄력 있는 피부로 만들어진 작은 물풍선이라고 상상해 보세요. 적혈구의 임무는 좁은 모세혈관을 통과하며 산소를 전달하는 것입니다. 보통 이 풍선은 납작한 도넛 모양(디스코사이트, discocyte)을 하고 있습니다.

하지만 이 세포들이 혈관을 통해 흐를 때, 유체가 세포를 밀어내는 힘이 작용합니다. 특정 속도에 도달하면, 이 세포들은 단순히 늘어나거나 회전하는 대신, 마치 여러 갈래의 꽃잎이나 비대칭적인 덩어리처럼 기이하고 복잡한 모양으로 갑자기 변형됩니다. 과학자들은 이런 현상을 관찰해 왔지만, 왜 이 풍선이 그렇게 급격하게 모양을 바꾸기로 결정했는지 그 이유를 완전히 이해하지 못했습니다.

이 논문은 이러한 갑작스러운 모양 변화 뒤에 숨겨진 규칙을 밝혀내기 위한 이론적 조사(수학적 시뮬레이션)입니다.

세포의 생애: 전과 후

이 불안정성을 이해하기 위해, 저자들은 특정한 생물학적 사건인 탈핵(enucleation) 과정을 고려해야 했습니다.

• 비유: 완벽하게 둥글고 팽팽한 스트레스 볼을 상상해 보세요. 이제 그 공 안의 공기를 빼내되, 고무 피부는 줄어들지 않게 한다고 가정해 봅시다. 그러면 내부의 공은 작아졌지만 피부는 여전히 크기 때문에, 피부가 너무 커져서 접히거나 주름이 생기게 됩니다.
• 과학적 설명: 적혈구가 골수에서 성숙할 때, 세포는 핵(명령 센터)을 밖으로 내보냅니다. 이 과정에서 세포 내부 크기는 작아지지만, 외부 피부의 크기는 그대로 유지됩니다. 이로 인해 "여분의 면적(excess area)" 즉, 약간의 느슨한 여분의 피부가 생겨납니다. 이 논문은 이 "여분(slack)"이 세포를 불안정하게 만들고, 혈류의 압력을 받을 때 모양을 변하게 만드는 핵심 요소라고 주장합니다.

실험: 풍선 밀기

연구진은 이 세포들이 전단 유동(shear flow)(두 개의 미끄러지는 유리판 사이에 세포가 끼어 있는 상황을 상상해 보세요)에 노출될 때 어떤 일이 일어나는지 시뮬레이션하기 위해 수학적 모델을 만들었습니다.

그들은 서로 맞서 싸우는 세 가지 주요 힘을 살펴보았습니다:

1. 유동(The Flow): 세포를 늘리고 회전시키려는 물의 힘.
2. 피부의 탄성(The Skin's Elasticity): 원래 모양으로 되돌아가려는 세포의 성질.
3. 굽힘 저항(The Bending Resistance): 피부를 급격하게 휘게 만드는 데 필요한 에너지.

"장력"의 줄다리기

이 논문의 핵심 발견은 **막 장력(membrane tension)**을 중심으로 이루어집니다. 이것을 풍선 피부의 팽팽함 정도로 생각하면 됩니다.

• 설정: 세포가 그냥 떠 있을 때는 장력이 균형을 이룹니다.
• 트리거(계기): 유속이 빨라짐에 따라 세포는 늘어납니다. 앞서 언급한 "여분의 면적(extra skin)" 때문에 피부의 장력은 떨어집니다.
• 스냅(급격한 변화): 장력이 너무 낮아지면, 피부는 더 이상 매끄러운 타원체 형태를 유지할 수 없습니다. 이는 마치 너무 느슨해진 드럼 가죽이 통제 불능으로 흔들리는 것과 같습니다.

저자들은 이 불안정성이 "여분(slack)"이 많고, 유체의 속도가 빠르며, 세포의 피부가 움직임에 저항할 만큼 딱딱하지 않고 충분히 유연할 때 발생한다는 것을 발견했습니다.

결과: 매끄러운 모양에서 뾰족한 모양으로

불안정성이 작동하면, 세포는 단순히 늘어나는 것이 아니라 "여분의 면적"을 재배치합니다.

• 비유: 매끄럽고 둥근 풍선을 상상해 보세요. 만약 고무가 너무 많은 상태에서 강하게 누른다면, 풍선은 단순히 납작해지는 것이 아니라 갑자기 몇 개의 돌기나 엽(lobe) 모양으로 튀어나올 수 있습니다.
• 발견: 수학적 모델은 세포가 단순한 타원형(럭비공 모양)에서 복잡한 다엽 구조(trilobe 또는 입 모양을 닮은 stomatocyte)로 변형됨을 보여줍니다. "여분의 면적"이 세포의 본체에서 새로운 날카로운 엽(lobe) 부분으로 이동하는 것입니다.

방향의 중요성

논문은 또한 세포가 처음에 어떤 방향으로 놓여 있는지가 중요하다는 점을 언급합니다.

• 만약 세포가 유동에 수직으로 정렬되어 있다면, 타이어처럼 구르는 동작(tank-treading)을 할 수 있습니다.
• 반대로 유동에 평행하게 정렬되어 있다면, 불안정성이 발생하기 더 쉬워져서 세포가 무질서하게 뒤집히거나 회전(tumbling/flipping)하게 됩니다.

결론

이 논문은 적혈구가 언제 매끄러운 모양을 잃고 복잡한 다엽 구조의 덩어리로 변하는지에 대한 수학적 "레시피"를 제공합니다. 핵심 재료는 다음과 같습니다:

1. 여분(The Slack): 세포가 성숙할 때(탈핵) 남겨진 여분의 피부.
2. 압박(The Push): 혈류의 속도.
3. 장력(The Tension): 피부가 모양을 유지하기에 너무 느슨해지는 순간.

저자들은 이 "느슨한 피부"가 적혈구가 극적인 모양 변화를 겪게 만드는 비밀 소스이며, 이는 세포의 기능에 필수적이면서도 동시에 역동적인 불안정성의 신호가 된다고 결론지었습니다.

기술 요약

기술 요약: 전단 유동 내 탈핵된 비압축성 적혈구: 형상 불안정성에 대한 이론적 분석

문제 정의
적혈구(RBC)는 산소 운반을 위해 필수적이며, 이는 상당한 변형을 일으킬 수 있는 능력에 의존한다. 강체 모델(예: 제프리 궤도)은 낮은 전단율에서의 적혈구 역학(텀블링, 탱크 트레딩, 스윙잉)을 성공적으로 설명하지만, 중간 전단율($100\text{--}1000 \text{ s}^{-1}$)에서 관찰되는 복잡한 다엽(multi-lobed) 형상 및 대칭성이 깨진 형태(스토마토사이트와 같은)의 출현을 설명하는 데는 실패한다. 이러한 형상 불안정성을 구동하는 물리적 메커니즘, 즉 막 탄성, 굽힘 저항, 막 장력, 과잉 면적 및 점도 비 사이의 상호작용은 여전히 논쟁의 대상이다. 더욱이, 기존의 이론적 프레임워크는 핵이 배출되는 성숙 과정으로 인한 기하학적 변화, 즉 응력-자유 면적과 갇힌 부피 사이의 관계 변화를 간과하고 적혈구를 정적인 타원체나 캡슐로 취급하는 경향이 있다.

방법론
저자들은 뉴턴 유체 전단 유동 내에 부유하는 폐쇄된 유체 충전 비압축성 막으로 적혈구를 모델링하는 섭동 이론적 프레임워크를 사용한다. 이 접근 방식은 Vlahovska 등[26]에 의해 이전에 확립된 선형 비선형 프레임워크를 확장하여, 두 가지 구별되는 과잉 면적을 명시적으로 포함한다:

1. $\Delta_0$ (응력-자유 과잉 면적): 탈핵 전 적혈구의 타원체 형상에 대응한다.
2. $\Delta$ (현재 과잉 면적): 탈핵 후 디스코이트(discocyte) 형상에 대응하며, 여기서 부피는 감소($V < V_0$)하지만 표면적은 보존되어 내부 응력이 고정된다.

적혈구 형상은 구면 조화 함수($Y_{l,m}$)로 전개되며, 현재의 구성 $\rho$와 응력-자유 구성 $\rho_0$를 분리한다. 운동 방정식은 유체 응력(스토크스 유동)과 막의 힘(탄성(Skalak 모델), 굽힘(Helfrich 에너지), 그리고 막 장력 기여분으로 분해됨) 사이의 균형을 맞춤으로써 유도된다. 막 장력($\Sigma$)은 국소 면적을 일정하게 유지하기 위한 라그랑주 승수로서 작용하며 서로 다른 공간 모드들을 결합시킨다.

분석은 두 단계로 진행된다:

1. 축소된 시스템 ($l=2$): $\Delta$와 $\Delta_0$ 사이의 불일치를 명시적으로 고려하여 텀블링, 탱크 트레딩, 스윙잉의 역학을 재검토하기 위해 시스템을 타원체 모드로 먼저 축소한다.
2. 전체 섭동 시스템 ($l \geq 2$): 선형 불안정성 기준을 도출하기 위해 추가적인 공간 모드를 도입한다. 시스템의 안정성은 무차원 매개변수인 굽힘 수($B$), 모세관 수($C$), 점도 비($\lambda$), 그리고 과잉 면적의 함수로서 고차 모드($l > 2$)의 성장률을 조사함으로써 분석된다.

주요 기여 및 결과

• 이론적 프레임워크의 확장: 본 연구는 탈핵에 의해 유도된 기하학적 및 응력 상태의 변화를 명시적으로 고려하여 기존의 섭동적 접근 방식을 일반화한다. 응력-자유 형상($\Delta_0$)과 현재 형상($\Delta$)을 구분함으로써, 이 모델은 유체 베시클 모델에는 없는 세포의 성숙 과정에 따른 "기억"을 포착한다.
• 축소된 역학에서의 과잉 면적의 역할: $l=2$ 하부 공간에서, 저자들은 $\Delta$와 $\Delta_0$ 사이의 불일치가 유효 탄성 상수($S$)를 재규격화함을 입증한다. 이를 통해 고정된 전단율에서도 참조 형상의 과잉 면적 변화에 의해 구동되는 역학적 전이(예: 스윙잉에서 텀블링으로)를 포착할 수 있다. 초기 방향(면내 또는 면외)은 결과적인 궤적에 상당한 영향을 미친다.
• 불안정성 기준 및 막 장력: 형상 전이를 위한 선형 불안정성 기준의 도출이 주요 기여이다. 분석 결과, 형상 불안정성(다엽 형상의 출현)은 막 장력 $\sigma(t)$에 결정적으로 의존함이 밝혀졌다. 불안정성은 막 장력이 임계값($\sigma < \sigma^*$) 아래로 떨어질 때 발생하며, 이 임계값은 전단 탄성과 굽힘 강성의 비율($\alpha = B/C$)에 의해 결정된다.
  • 불안정성은 굽힘 및 탄성 반응이 외부 유동 및 내부 응력에 대항하기에 불충분할 때 시작된다.
  • 가장 불안정한 모드($l^*$)와 임계 장력은 $\alpha$, 점도 비 $\lambda$, 그리고 과잉 면적에 따라 달라진다.
• 포화 및 재분배: 불안정성 발생 이후의 장기 거동을 분석한다. 불안정한 모드의 성장은 막 장력을 증가시키며, 이는 결국 시스템을 안정화시킨다. 이 과정은 과잉 면적을 지배적인 $l=2$ 모드에서 고차 모드(예: $l=3, 4$)로 재분배하여, 새로운 안정적인 비타원체 형상 구성을 초래한다.
• 위상도: 수치 적분을 통해 분석적 예측을 검증한다. 저자들은 $(B, C)$ 평면에서의 위상도를 구축하여 특정 모드($l=3, 4, \dots$)가 불안정해지는 영역을 식별한다. 이 도표들은 점도 비 $\lambda$와 과잉 면적 $\Delta$를 증가시키는 것이 높은 모세관 수(약한 탄성)에서 적혈구를 불안정하게 만들어 복잡한 형상의 출현을 촉진함을 보여준다.

의의
본 논문은 외부 전단 유동 하에서 스토마토사이트 및 트리로브(trilobe) 형상과 같은 복잡한 적혈구 형태가 나타나는 것에 대한 이론적 설명을 제공한다. 적혈구의 성숙 과정(응력-자유 과잉 면적과 현재 과잉 면적 사이의 불일치)에 따른 구체적인 기계적 결과와 이 불안정성을 연결함으로써, 본 연구는 적혈구가 단순한 타원체 역학에서 복잡한 다엽 상태로 전이되는 메커니즘에 대한 이해를 제공한다. 이 프레임워크는 막 장력이 공간 모드들을 결합하는 핵심 동력임을 식별하며, 굽힘-전단 계수 비와 점도가 형상 불안정성의 발생을 결정하는 데 있어 결정적인 역할을 한다는 점을 강조한다. 이러한 이론적 그림은 강체 모델과 전체 수치 시뮬레이션 사이의 간극을 메우며, 적혈구가 타원체 대칭성을 잃게 되는 조건에 대한 분석적 통찰을 제공한다.