Evolution of the SARS-CoV-2 main protease under fitness constraints imposed by folding stability and activity

기술 요약: 적응도 제약 조건 하에서의 SARS-CoV-2 메인 프로테아제 진화

문제 정의
백신의 성공에도 불구하고, SARS-CoV-2는 보완적인 항바이러스 전략을 필요로 하는 면역 회피 변이를 생성하며 계속해서 진화하고 있다. 메인 프로테아제(Mpro, 또는 3CLpro 또는 nsp5로 알려짐)는 바이러스 복제에서의 필수적인 역할, Coronaviridae 패밀리 전반에 걸친 높은 보존성, 그리고 유사한 인간 상동체의 부재로 인해 중요한 치료 표적이다. 그러나 Mpro 변이가 기질 인식, 접힘 안정성(folding stability) 및 잠재적인 약물 내성에 미치는 진화적 결과는 아직 충분히 탐구되지 않았다. Mpro는 고도로 보존되어 있지만, 특정 바이러스 계통(예: Omicron 관련 P132H)을 중심으로 아미노산 치환이 축적된다. 본 연구가 다루는 핵심 질문은 이러한 널리 퍼진 돌연변이들이 효소의 구조적 안정성, 촉매 기능 또는 기질 결합 친화력을 저해하는지 여부와, 진화적 제약이 단백질의 적응도 지형(fitness landscape)을 어떻게 형성하는가이다.

연구 방법론
본 연구는 진화 유전학과 구조 생물물리학을 통합한 다중 스케일 계산 접근 방식을 채택하였다:

1. 데이터 수집 및 변이 특성 분석: 저자들은 NCBI에서 다운로드한 101,000개의 SARS-CoV-2 ORF1ab 아미노산 서열을 분석하였다. 서열은 길이와 동일성(Wuhan 참조 서열 대비 >75%)을 기준으로 필터링되었으며, MAFFT를 사용하여 정렬되었다. Mpro 서열은 추출된 후 돌연변이 프로필과 시간적 빈도(2019–2025)에 따라 그룹화되었다.
2. 구조 모델링: 야생형(WT) Mpro와 가장 빈번하게 나타나는 10가지 변이체를 세 가지 천연 기질(nsp4-5, nsp5-6, nsp10-11)과 결합된 상태로 AlphaFold3 (v3.0.1)를 사용하여 모델링하였다.
3. 분자 동역학(MD) 시뮬레이션: 모든 원자 모델 MD 시뮬레이션은 AMBER22와 ff14SB 포스 필드, TIP3P 용매를 사용하여 수행되었다. 시스템은 에너지 최소화, 평형화(NVT 및 NPT), 그리고 10 ns의 생성 실행(production run) 과정을 거쳤다. 궤적은 백본 RMSD와 수소 결합(HB) 점유율을 분석하는 데 사용되었다.
4. 결합 에너지 분석: 결합 자유 에너지($\Delta G$)는 MM-GBSA(Molecular Mechanics-Generalized Born Surface Area) 방법을 사용하여 추정되었다. 주요 잔기가 기질 결합에 기여하는 정도를 식별하기 위해 잔기별 에너지 분해(per-residue energy decomposition)를 수행하였다.
5. 컨포메이션 역학 분석: C$\alpha$ 원자 좌표에 주성분 분석(PCA)을 적용하여 지배적인 집단 운동을 특성화하였다. 컨포메이션 랜드스케이프는 첫 두 주성분(PC1 및 PC2)에 대한 2차원 커널 밀도 추정(KDE)을 통해 시각화되었다.
6. 정량적 발산 측정: 쿨백-라이블러(KL) 발산을 계산하여 변이체들과 WT 사이의 컨포메이션 샘플링 확률 분포를 정량적으로 비교함으로써, 컨포메이션 엔트로피의 변화를 평가하였다.

주요 결과

• 진화 패턴: WT Mpro는 2021년 말까지 우세했다. Omicron 관련 P132H 변이는 2022년 초부터 2025년까지 우세하게 나타났다. 다른 돌연변이들(예: L89F, K90R)은 일시적이거나 낮은 빈도의 증가를 보였다. 결정적으로, 활성 부위 및 이량체 계면(dimerization interface)에서의 돌연변이는 극히 드물거나 존재하지 않았는데(빈도 일반적으로 <0.001), 이는 이 기능적 영역들에서 강력한 정화 선택(purifying selection)이 작용하고 있음을 나타낸다.
• 기질 결합 친화도: MM-GBSA 계산 결과, 널리 퍼진 치환들은 세 가지 천연 기질 모두에 대해 전체적인 기질 결합 친화력에 실질적인 변화를 일으키지 않았다. 특정 변이-기질 쌍(예: V73I와 nsp4-5)에서 통계적으로 유의미한 차이가 발견되었으나, 이러한 효과는 국소적이었으며 의미 있는 적응도 차이를 시사하기에는 불충분했다.
• 상호작용 네트워크: 수소 결합 및 잔기별 기여도에 대한 상세 분석 결과, 전반적인 결합 패턴과 HB 네트워크는 대체로 보존되어 있었다. 특정 국소적 섭동(예: V73I 변이와 nsp4-5 V4 사이의 상호작용 감소)이 발생했음에도 불구하고, 핵심 활성 부위 잔기(G23, T24, T25, T26, H41, C145 등)를 포함하는 핵심 상호작용 네트워크는 안정적으로 유지되었다.
• 컨포메이션 역학: PCA 결과, WT와 변이체 모두 내재적인 단백질 유연성을 반영하는 여러 메타스테이블(metastable) 컨포메이션 하부 상태를 샘플링함을 보여주었다.
  • PC1 (글로벌 운동): 모든 변이체에서 PC1에 대한 KL 발산이 0으로 나타났으며, 이는 지배적이고 큰 진폭의 글로벌 집단 운동이 모든 변이체에서 엄격하게 보존됨을 의미한다.
  • PC2 (이차 운동): PC2에서 완만한 KL 발산(약 0.08에서 0.76 사이)이 관찰되었으며, 이는 특정 치환이 전체적인 역학 구조를 무너뜨리지 않으면서 이차 운동 및 국소 역학을 미세하게 조정함을 시사한다.

주요 기여

1. 포괄적 적응도 평가: 본 연구는 가장 널리 퍼진 SARS-CoV-2 Mpro 변이들이 접힘 안정성과 기질 결합에 어떻게 영향을 미치는지에 대한 체계적인 평가를 제공하며, 현재 순환하는 변이들이 높은 기능적 적응도를 유지하고 있음을 확인하였다.
2. 역학적 보존성: PCA와 KL 발산을 결와하여, 저자들은 Mpro의 진화가 제약되어 있으면서도 기능에 필수적인 글로벌 역학적 행동(PC1)을 보존하는 동시에, 이차 운동(PC2)에서만 미세한 조정을 거친다는 것을 입증하였다.
3. 에피스타시스(Epistasis) 고려: 분석에는 다중 돌연변이(예: K90R_P132H)가 포함되어, 단일 돌연변이 연구에서는 놓칠 수 있는 잠재적인 에피스타시스 효과를 탐지할 수 있었다.
4. 방법론적 프레임워크: 본 논문은 진화적 감시와 생물물리학적 모델링(AlphaFold3, MD, MM-GBSA, PCA)을 결합하여 바이러스 단백질의 적응을 모니터링하는 워크플로우를 확립하였다.

의의 및 주장
본 논문은 SARS-CoV-2 Mpro의 진화가 구조적 및 기능적 요구 사항, 특히 기질 인식 및 촉매 능력을 유지해야 한다는 필요성에 의해 강력하게 제약받고 있다고 주장한다. 연구 결과는 Mpro가 견고한 항바이러스 표적로서 적합하다는 점을 뒷받침하며, 바이러스는 아직 천연 기질과의 상호작용이나 글로벌 컨포메이션 역학을 유의미하게 저해하는 변이를 진화시키지 못했다.

저자들은 현재의 변이들이 Mpro 억제제의 효능을 위협하지는 않지만, 미래의 내성 경로를 예측하기 위해 취약할 수 있는 특정 국소 접촉 부위를 식별하는 것이 중요하다고 주장한다. 저자들은 진화 역학을 생물물리학적 적응도 평가(안정성, 결합 에너지 및 컨포메이션 앙상블)와 통합하는 것이 더 나은 장기적 견고성을 갖춘 치료 전략을 설계하고 바이러스 진화를 모니터링하는 실질적인 프레임워크를 제공한다고 제안한다. 이 연구는 내성이 임박했다고 주장하는 것이 아니라, 컨포메이션 앙상블이나 결합 친화력을 약물 회피에 유리한 방향으로 변화시킬 수 있는 변이를 탐지하기 위한 지속적인 감시의 필요성을 강조한다.