Imaging aerosolized viruses with an X-ray free-electron laser using single-particle rotational invariants

기술 요약: 회전 불변량을 이용한 X선 자유 전자 레이저 기반 에어로졸화 바이러스 이미징

문제 제기
X선 자유 전자 레이저(XFEL)를 이용한 단일 입자 이미징(SPI)은 비결정성 생체 입자(예: 바이러스)를 생물학적으로 유사한 조건 하에서 "파괴 전 회절(diffraction-before-destruction)" 측정을 가능하게 한다. 그러나 비결정성 표적의 경우 병진 불변성(translational invariance)의 결여로 인해 구조 결정 과정이 매우 복잡해진다. 기존의 접근 방식은 직접 접근 가능한 관측량인 산란 강도 $I(q)$에 의존하며, 실공간 전자 밀도를 재구성하기 위해 손실된 위상을 복구해야 한다. 반복적 위상 결정 알고리즘이 존재하지만, 실험적 해상도가 제한적이거나 대칭성 제약이 필요한 경우 어려움을 겪을 수 있으며, 이는 고유한 구조적 비대칭성과 역동성을 가릴 위험이 있다. 또한, 에어로졸화된 입자 앙상블 내의 구조적 변이를 특성화하는 것은 무작위 방향 정렬의 복잡성과 환경적 이질성으로 인해 여전히 어려운 과제이다.

방법론
본 연구는 에어로졸화된 박테리오파지 PR772의 모델 가이드 및 ab initio 3차원(3D) 구조 결정을 용이하게 하기 위해 X선 회절 측정으로부터 유도된 **회전 불변량(rotational invariants)**을 활용한다.

1. 데이터 획득: 에어로데믹 렌즈 스택을 통해 에어로졸화된 PR72 입자를 사용하여 Linac Coherent Light Source (LCLS)에서 펨토초 XFEL 펄스(1.6 및 1.7 keV)로 회절 패턴을 기록하였다(amo06516, amo86615, amo11416의 세 가지 독립적인 실험).
2. 회전 불변량을 통한 XCCA: 저자들은 개별 회절 패턴에만 의존하는 대신, 각도 X선 교차 상관 분석(XCCA)을 활용하였다. 신호 대 잡음비를 높이고 배경 노이즈를 억제하기 위해 평균 차이 패턴 이점 각도 교차 상관 함수(CCF)인 $C_{diff}(q, q', \Delta)$를 계산하였다.
3. 불변량 추출: CCF는 $C_n(q, q')$로 푸리에 성분 분해되며, 이는 회전 불변량 $B_l(q, q')$의 선형 조합이다. 회전 불변량은 구면 조화 함수 확장 계수의 곱으로 정의되며, 입자의 방향과 무관한 고유한 구조 정보를 인코딩한다.
4. 구조 결정:
   • Ab Initio 재구성: 실험적으로 결정된 불변량 $B_l(q, q')$로부터 3D 전자 밀도를 재구성하기 위해 xFrame 소프트웨어 제품군에 구현된 다단계 반복 위상 결정(MTIP) 알고리즘을 사용하였다. 이 과정은 대칭성 제약을 부과하지 않고 미지의 유니터리 행렬(방향 결정)과 위상을 동시에 해결한다.
   • 모델 가이드 분석: 다양한 구조 모델(예: 완벽한 정이십면체, 왜곡된 캡시드, 코어-쉘 구조)로부터의 회절을 시뮬레이션하기 위해 비드 모델링(bead-modeling) 접근법을 채택하였다. 이러한 시뮬레이션은 이론적 불변량을 생성하여 실험 데이터와 비교함으로써 밀도 불균일성이나 형상 왜곡과 같은 특정 구조적 특징을 검증할 수 있게 한다.
5. 정규화: 노이즈와 수치적 오류를 완화하고 MTIP 알고리즘에 최적화된 데이터를 제공하기 위해 Ruiz 평형화(equilibration) 및 저계수 근사(low-rank approximation)를 포함한 고급 정규화 기술을 적용하였다.

주요 결과

• 구조적 비대칭성: 완만한 실험적 해상도(데이터셋 및 메트릭에 따라 6 ~ 13 nm 범위)에도 불구하고, ab initio 재구성은 PR772 캡시드가 완벽한 정이십면체 대칭에서 벗어나 있음을 일관되게 보여주었다. 이러한 왜곡은 세 가지 독립적인 실험 모두에서 관찰되었으며, 이는 특정 실험 조건의 산물이 아닌 공통적인 구조적 특성임을 시사한다.
• 내부 밀도 변이: 재구성과 모델 가이드 분석은 바이러스 입자 내의 비균일한 밀도 분포를 나타냈다. 구체적으로, amo06516 데이터셋은 amo86615 데이터셋에 비해 상당한 내부 밀도 이질성을 보였으며, 후자는 더 균일한 것으로 나타났다.
• 꼭짓점 확장: amo06516 서브셋에서 정이십면체 꼭짓점 중 하나에서 솟아나온 희미한 밀도 확장이 감지되었다. 이 특징은 관련 바이러스인 PRD1에서 게놈 전달 시 관찰되는 나노튜브 돌출부와 유사하며, 이는 구조적 변형이나 샘플의 이질성을 시사한다.
• 불변량의 견고성: 회전 불변량 $C_n(q, q')$은 높은 견고성을 보였으며, 회절 패턴의 수가 달라지더라도 서로 다른 실험 실행 및 입자 크기 서브셋 간에 강한 유사성을 나타냈다.
• 모델 검증: 시뮬레이션 결과 완벽한 정이십면체 모델은 실험적 불변량 형태를 재현하는 데 실패했다. 반면, 반경 방향 단축 압축(radial uniaxial compression)과 복합 왜곡(반경 및 접선 방향)을 도입했을 때, 특히 낮은 차수의 불변량($n=2, 4$)에 대해 실험적 위상 맵과의 일치도가 크게 향상되었다. 코어-쉘 모델은 내부 코어 밀도에 대한 민감도를 드러내며 적합도를 더욱 정밀하게 개선하였다.

의의 및 주장
저자들은 회전 불변량이 XFEL을 사용하여 나노 및 생체 입자의 구조와 시간적 진화를 조사하는 강력한 도구 역할을 한다고 주장한다. 회전 자유도를 제거함으로써, 이 불변량들은 순방향 모델링(forward modeling)과 역문제 해결(inverse problem solving)에 동시에 사용할 수 있는 광범위한 구조 정보를 제공한다.

본 논문은 ab initio 재구성과 모델 가이드 분석을 결합한 접근 방식이 실험적 해상도가 제한될 때 구조 분석의 충실도를 높인다는 점을 강조한다. 결정적으로, 이 방법은 전통적인 결정학이나 저온 전자 현미경(cryo-EM) 분석에서 가려질 수 있는 구조적 특징(예: 비대칭 캡시드 왜곡 및 밀도 불균일성)을 탐지할 수 있게 한다. 결과는 회전 불변량이 측정된 앙상블 내의 유사성과 구조적 변이를 효과적으로 특성화할 수 있음을 보여주며, 이는 근접 자연 상태(near-native conditions)에서 바이러스 전파 경로와 구조적 역동성을 연구하는 경로를 제시한다. 본 연구는 회전 불변량이 상관 관계 기반 이미징 기술을 통해 바이러스의 고유한 비대칭성과 동적 구조 변이(잠재적으로 에어로졸화 과정 중의 환경 변화 또는 고유한 평형 역동성에 의해 유도된 것)를 특성화할 수 있음을 입증한다.