A geometry-dependent, force balance-driven model of Staphylococcus epidermidis biofilm cell cluster detachment

기술 요약: 스테파로코커스 에피더미디스 (Staphylococcus epidermidis) 바이오필름 세포 군집 박리의 기하학적 의존성 및 힘 균형 기반 모델

문제 제기
병원 내 지속성 감염의 주요 원인인 Staphylococcus epidermidis에 의해 형성된 박테리아 바이오필름은, 박테리아 군집이 박리되어 순환하고 새로운 감염 부위를 점령할 수 있는 분해 및 분산 단계를 포함하는 생명주기를 가집니다. 이러한 박리된 군집의 생물학적 운명, 즉 면역 제거를 회피하고 혈관 벽에 부착하며 항생제에 저항하는 능력은 그 크기와 형태에 크게 영향을 받습니다. 기존 바이오필름 박리 수학적 모델들은 종종 이 과정을 바이오필름 두께, 전단 응력, 또는 세포외 고분자 물질 (EPS) 밀도의 함수로 단순화합니다. 이러한 프레임워크들은 종종 예측된 생물학적 결과에 결정적인 기하학적 구조와 공간적 배열과 같은 박리된 군집의 특정 특성을 추적하지 못합니다. 또한, 세포 수준에서 시뮬레이션된 바이오필름을 검증하기 위해 실험적 미세구조 데이터를 통합하는 모델은 거의 없습니다.

방법론
저자들은 바이오필름 성장 시뮬레이션과 새로운 기하학적 의존성 박리 알고리즘을 결합한 통합 모델링 프레임워크를 개발했습니다. 이 접근법은 다음 세 가지 주요 단계로 진행됩니다:

1. 바이오필름 성장 (iDynoMiCS): 저자들은 오픈소스 에이전트 기반 플랫폼인 iDynoMiCS 를 활용하여 24 시간 된 S. epidermidis 바이오필름의 성장을 시뮬레이션했습니다. 이 모델은 정지된 구형 박테리아 에이전트와 EPS 에이전트를 추적하며, 이들의 성장과 분열을 모노드 (Monod) 동역학을 통해 용질 농도 (산소와 탄소) 와 결합시킵니다. 생물학적 관련성을 보장하기 위해 시뮬레이션 매개변수는 바이오필름 높이, 세포 수 밀도, 세포 조직 (군집화), 그리고 EPS 바이오매스에 관한 실험 데이터와 비교 검증되었습니다. 침식률과 EPS 수율을 변화시키기 위해 네 가지 서로 다른 매개변수 세트가 생성되어 서로 다른 구조적 특성을 가진 바이오필름을 만들었습니다.

2. 군집 형성 알고리즘: 잠재적 박리 사건을 식별하기 위해 저자들은 4 단계 알고리즘 (태깅, 추가, 선별, 군집화) 을 도입했습니다.

   • 태깅: 에이전트 (박테리아 및 EPS) 는 국소 높이, 이웃 밀도, 그리고 "점착성" 매개변수 ($\beta$) 를 기반으로 박리 확률이 할당됩니다. 이 점착성은 효소적 분해 또는 화학적 교란제로 인한 손상된 EPS 바이오매스의 비율을 나타내는 교란 비율 ($r_{H^*}$) 에 의해 조절됩니다.
   • 추가 및 선별: 태그가 지정된 이웃으로 둘러싸인 태그가 지정되지 않은 에이전트는 집합에 추가되는 반면, 바이오필름 깊숙이 위치하여 충분한 태그가 지정된 이웃이 없는 에이전트는 제거됩니다.
   • 군집화: 남은 태그가 지정된 에이전트들은 DBSCAN 알고리즘을 사용하여 군집으로 그룹화됩니다. 군집 내 에이전트의 기하학적 배열은 타원 (2 차원) 또는 타원체 (3 차원) 를 적합시켜 종횡비와 방향 각도를 결정함으로써 특징지어집니다.

3. 힘 균형 박리 모델: 식별된 각 군집에 대해, 항력 모멘트 ($M_d$) 와 접착 모멘트 ($M_a$) 사이의 균형에 따라 박리가 발생하는지 여부가 결정됩니다.

   • 항력 모멘트 ($M_d$): 낮은 레이놀즈 수 유동에 대한 유체 역학 원리를 사용하여 계산되며, 군집의 종횡비, 유동 방향에 대한 방향 각도, 그리고 단면적을 포함합니다.
   • 접착 모멘트 ($M_a$): 손상된 EPS 영역에서의 접착 약화를 고려하여 군집과 주변 에이전트 간의 상호작용 포텐셜을 계산하기 위해 수정된 Gay-Berne 포텐셜 이론을 사용하여 모델링됩니다.
   • 결정 규칙: $M_d > M_a$인 경우 박리가 발생합니다.

이 모델은 군집과 단일 에이전트가 제거됨에 따라 바이오필름 지형이 진화하도록 1 시간 동안 1 분 간격의 시간 단계로 박리를 시뮬레이션합니다.

주요 결과
이 연구는 2 차원 및 3 차원 바이오필름에 대해 다양한 EPS 교란 수준 (0% 에서 95%) 에서 박리를 시뮬레이션했습니다:

• 군집 크기와 빈도: 대부분의 박리된 군집은 작았으며, 직경은 일반적으로 7 µm 미만이었습니다. EPS 교란이 증가함에 따라 박리 사건의 빈도가 증가했고 박리된 군집의 평균 크기가 약간 증가했으나, 대부분 면역 회피와 관련된 5–10 µm 임계값 아래에 머물렀습니다.
• 기하학과 방향: 더 큰 박리된 군집은 더 작고 구형인 군집에 비해 2 차원에서는 더 선형적으로, 3 차원에서는 더 이방성으로 경향했습니다. 3 차원 시뮬레이션에서 대부분의 박리된 군집은 장타원체형이었습니다. 흥미롭게도, 이 모델은 군집이 종종 유동 방향에 대해 큰 방향 각도에서 박리된다고 예측했는데, 이는 단일 입자에 대한 일부 이론적 기대와 대조되지만 바이오필름 응집체의 복잡한 기하학과 일치합니다.
• 교란의 영향: 더 높은 수준의 EPS 교란은 무교란 (0% 교란) 기준선에 비해 바이오필름 높이의 더 빠른 감소와 박리된 군집의 더 높은 비율로 이어졌습니다. 그러나 95% 교란에서도 평균 박리된 군집 크기는 5.5 µm 를 초과하지 않았습니다.
• 매개변수 민감도: 결과는 서로 다른 매개변수 세트 (침식률과 EPS 수율 변화) 에서 견고했으며, 일반적으로 더 높은 바이오필름이 더 큰 군집의 더 높은 비율을 생성했습니다.

의의 및 주장
이 논문은 세 가지 주요 신규성을 주장합니다:

1. 실험적 검증: 이전의 많은 모델들과 달리, 이 프레임워크는 S. epidermidis에 대한 실험 데이터와 비교하여 시뮬레이션된 바이오필름 미세구조 (높이, 밀도, 군집화, EPS 함량) 를 명시적으로 검증하여 초기 조건이 생물학적으로 현실적임을 보장합니다.
2. 기하학적 의존성 역학: 이 모델은 박테리아와 EPS 의 공간적 위치를 국소 접착 강도와 명시적으로 연결하고, 어떤 세포 그룹이 함께 박리될지 결정하기 위해 기하학적 요소 (형태, 방향) 를 활용함으로써 단순화된 박리 함수를 넘어섭니다.
3. 시간적 역학: 이 프레임워크는 군집 빈도, 크기, 그리고 모양이 EPS 교란에 반응하여 시간에 따라 어떻게 변화하는지에 대한 기계론적 통찰력을 제공하며, 이는 실험적으로 얻기 어려운 데이터입니다.

저자들은 이 통합 모델링 프레임워크를 바이오필름 박리 과정의 예측 능력에서 중요한 진전으로 위치시킵니다. 그들은 손상된 EPS 가 군집 형성 및 박리 역학에 미치는 영향을 이해하는 것이 면역 제거에 더 취약하거나 의료 기기에 부착할 가능성이 적은 특정 군집 크기를 생성하기 위한 바이오필름의 표적 교란과 같은 치료 전략에 정보를 제공할 수 있다고 제안합니다. 이 논문은 일정한 EPS 교란에 대한 가정과 군집을 완벽한 타원체로 단순화하는 것 등 한계를 인정하지만, 현재 모델이 바이오필름 분해의 복잡한 역학을 예측하기 위한 실현 가능한 첫걸음을 제공한다고 주장합니다.