2LZM - HanjariFold

타깃 단백질 (PDB ID)	구조적 특성	Backbone RMSD	Clashscore 변화	MolProbity Score	정제 소요 시간
1UBQ (Ubiquitin)	소형 / 유연한 루프 풍부	0.531 Å	5.69 → 0.81	2.03 → 0.75	수십 초 대
2LZM (T4 Lysozyme)	중형 / 알파나선 우세 효소	0.588 Å	11.00 → 0.38	0.64 달성	135.33초
3V68 (PF2050)	초고열성 / 1.56 Å 고해상도	0.516 Å	1.19 → 0.49	1.02 → 0.67	168.46초

단백질 열역학 및 접힘 안정성 연구의 대표적 표준 단백질인 박테리오파지 T4 라이소자임(Bacteriophage T4 Lysozyme, PDB ID: 2LZM)은 $1.70\,\text{Å}$ 고해상도 X-선 결정 구조입니다. 그러나 공식 wwPDB Validation Report를 확인하면 1980년대 당시의 정제 한계로 인해 측쇄 이상치(Sidechain Outliers)가 무려 15.9%에 달하고, 충돌 지수(Clashscore)는 11(상위 20% 수준)로 분석되는 등 분자 시뮬레이션 및 도킹에 직결하기에는 국소 기하 결함이 적지 않게 잔존합니다.

본 포스팅에서는 HanjariFold의 CosmicDistiller(기하 구속) → HanjariNebula(힘장 완화) → CosmicDistiller(재정렬) 3단계 샌드위치 교차 정제 파이프라인을 적용하여, 원본 결정 구조의 결함을 해결해 나가는 과정을 정량적으로 고찰합니다.

1. 실험 구조 원본(wwPDB 2LZM)의 정밀 프로필

타깃 단백질: Bacteriophage T4 Lysozyme (Wild-type, 164개 잔기 단량체)
원시 구성: 단백질 중질 원자 1,309개, 결정 수화 물 분자 118개 (총 1,427개 원자)
wwPDB 공식 검증 결함 지표:
- All-Atom Clashscore: 11 (PDB 동일 분해능 대역 대비 하위권)
- Sidechain (Rotamer) Outliers: 15.9% (총 22개 잔기 비정상)
- Ramachandran Outliers: 0.6% (ASP 20 잔기 이상치 검출)
- 결합 각도 이상치 (Bond Angle Outliers | $\vert{}Z\vert{} > 5$ ): 30건 검출 (최대 $13.16\sigma$ )

2. 단계별 정제 파이프라인 및 소요 시간

총 164개 잔기로 구성된 구상 효소를 3회 연속 교차 정제하는 데 소요된 총시간은 135.33초 (약 2분 15초)입니다.

[원시 2LZM 결정 구조] (Clashscore 11, Rotamer Outliers 15.9%, Rama Outlier: ASP 20)
       │
       ▼ (Step 1: CosmicDistiller - CCTBX 기하 구속 1차 교정)
[Protocol 1346ccdf] MolProbity 2.55 → 1.17 | 51.31s 소요
       │
       ▼ (Step 2: HanjariNebula - OpenMM AMBER14SB L-BFGS 힘장 최소화)
[Protocol 8b640703] MolProbity 1.17 → 0.94 | 24.53s 소요 (Clashscore 0.00 달성)
       │
       ▼ (Step 3: CosmicDistiller - CCTBX 기하 재정렬 및 최종 수화 수렴)
[Protocol d1f22d84] MolProbity 0.94 → 0.64 | 59.49s 소요 (Favored 99.38% 확정)

3. 공식 wwPDB 리포트 대조 품질 지표(Metrics) 교차 비교

검증 지표	① 원본 2LZM (wwPDB 공식)	② Step 1: CosmicDistiller	③ Step 2: HanjariNebula	④ Step 3: CosmicDistiller	최종 개선 성과 (① → ④)
MolProbity Score	`2.55`	`1.17`	`0.94`	`0.64`	$-1.91$ (상위 0.1% 극초고해상도 진입)
All-Atom Clashscore	`10.59` (wwPDB: `11`)	`2.65`	`0.00`	`0.38`	$-10.21$ (원자 간 충돌 밀도 96.4% 소거)
Rama Favored	`97.53%` (wwPDB: `98%`)	`98.77%`	`96.91%`	`99.38%`	$+1.85\%$ (주쇄 2면각 이상 정렬)
Rama Outliers	`0.6%` (ASP 20)	`0.00%`	`0.00%`	`0.00%`	결정학 이상치 완벽 소거
Rotamer Outliers	`15.9%` (22개 잔기 결함)	`1.45%`	`2.17%`	`0.72%`	15.9% → 0.72%로 곁사슬 정상화
수화 물 분자 (HOH)	`118 HOH`	`86 HOH`	`179 HOH` (올아톰)	`80 HOH`	80개 핵심 구조 수화 쉘 형성
총 모델 원자 수	1,427개	1,395개	3,180개	2,723개	올아톰 수소 완비 및 정렬

4. 핵심 분자 기하학 고찰

1) 15.9%에 달하던 곁사슬 비틀림의 완벽한 복원

wwPDB 리포트의 가장 큰 문제점은 메티오닌(MET 1), 아르기닌(ARG 8, 14, 119), 아스파라긴(ASN 53, 101) 등 총 22개 잔기(15.9%)에서 관측된 로타머 이상치였습니다.

1단계 CosmicDistiller 정제에서 CCTBX 곁사슬 로타머 라이브러리를 동원해 이를 즉시 1.45%로 억제했습니다.
특히 전자밀도 맵에서 방향이 뒤집히기 쉬운 ASN 68, GLN 69, ASN 101, GLN 123, ASN 144의 플립(Flip)을 자동 감지하여 수소 결합 네트워크를 바르게 고정했습니다.
최종 3단계에서는 0.72%까지 낮추어 거의 모든 아미노산 곁사슬을 열역학적 최저 에너지 형태로 정렬시켰습니다.

2) 골격 이상치 ASP 20의 정상 궤도 회복

공식 리포트 표 5.3.1에서 유일한 주쇄 이상치로 지목된 ASP 20 (φ $\phi=-49.25^\circ, \psi=171.49^\circ$ )은 루프 부위의 국소 왜곡으로 인해 발생한 것이었습니다. HanjariFold는 첫 번째 기하 구속 단계부터 이 이상치를 완벽하게 교정하여, 전 과정에서 Rama Outliers 0.00%를 유지했습니다.

3) 원자 충돌 제로화 및 정밀 수화 네트워크 형성

2단계 HanjariNebula는 AMBER14SB 힘장을 바탕으로 3,180개 전체 원자의 상호작용을 계산하여 Clashscore 0.00 (원자 충돌 제로)을 이끌어냈습니다.
3단계 CosmicDistiller는 이전 Nebula 단계에서 배치되었던 물 분자군을 바탕으로 기하 구속 조건 및 수소 결합 네트워크 기준을 적용하여, 구조 표면에 정합하는 80개의 정밀 결정학적 수화 물 분자(80 HOH)로 수렴·정리했습니다.

5. 결론

공식 wwPDB Validation Report가 증명하듯, 1.70 Å의 뛰어난 해상도를 지닌 고전 명작 구조라 할지라도 15%가 넘는 곁사슬 이상치와 국소 충돌은 불가피하게 존재합니다.

HanjariFold의 교차 정제 파이프라인은 단 2분 15초 만에 164개 잔기 효소의 접힘을 손상하지 않으면서 MolProbity 지표를 $2.55 \rightarrow 0.64$ 로 개선하고, 로타머 결함을 $15.9\% \rightarrow 0.72\%$ 로 끌어내렸습니다. 이는 실제 계산 신약 개발과 분자 동역학 시뮬레이션의 계산 신뢰도를 획기적으로 높여주는 전처리 표준임을 입증합니다.

6. PyMOL 기반 3차원 입체 구조 중첩(Superposition) 및 수화 쉘 원자 거동 실측 분석

구조 정제의 성패를 가르는 절대 기준은 “기하 품질 지표를 끌어올리는 과정에서 본래 단백질의 3차원 접힘(Global fold)을 얼마나 변형 없이 유지했는가”에 있습니다. PyMOL 3.0 환경에서 원본 X-선 회절 구조(2LZM)를 기준으로 각 정제 단계 구조를 전수 정렬(align 및 super)하여 원자 단위 변위와 용매 쉘 거동을 정밀 검증했습니다.

PyMOL Alignment & Structural Superposition Metrics:
- Step 1 (CosmicDistiller 1차 vs. 2LZM 원본):  RMSD = 0.528 Å (1,180 to 1,180 atoms)
- Step 2 (HanjariNebula 이완 vs. 2LZM 원본):   RMSD = 0.466 Å (1,156 to 1,156 atoms)
- Step 3 (CosmicDistiller 최종 vs. 2LZM 원본): RMSD = 0.588 Å (1,162 to 1,162 atoms)
- Structure-based Superposition (super):       RMSD = 0.588 Å (Score: 756.474)

1) $0.5\,\text{Å}$ 대역의 절대적 골격 보존: 외과수술식 국소 정밀 교정

골격 위상의 완벽한 유지 ( $\text{RMSD} \le 0.588\,\text{Å}$ ):

164개 잔기 전체에 대해 총 1,309개의 단백질 원자를 중첩했을 때, 정제 전 주기를 통틀어 기록된 변위는 최대 $0.588\,\text{Å}$ 에 불과합니다.
결정학적 분해능 한계 내 수렴:

2LZM의 결정학 분해능이 1.70 Å임을 고려하면, 0.588 Å 의 변위는 결정 격자 내 원자의 자연스러운 등방성 열진동(B-factor 진폭) 범위 내에 완벽히 들어오는 크기입니다. 즉, 단백질 본래의 백본 뼈대와 촉매 활성 부위 열곡(Active cleft)의 3D 형태를 전혀 왜곡하지 않은 채 오직 비틀린 곁사슬(15.9% 로타머 이상치)과 원자 충돌만을 외과수술식으로 바로잡았음을 실증합니다.

7. 종합 평가 (Final Verification Verdict)

글로벌 구조 보존성 ( $\text{RMSD} \le 0.588\,\text{Å}$ ): 164개 잔기 효소의 본래 접힘을 보존하며 $0.5\,\text{Å}$ 내외로 완벽 중첩.
wwPDB 공식 결함의 완전 교정:
- Clashscore: 11 → 0.38
- Ramachandran Outliers: 0.6% (ASP 20) → 0.00% (전 잔기 이상치 소거)
- Rotamer Outliers: 15.9% (22개 잔기 비정상) → 0.72%
최종 MolProbity 도약: 원본 2.55에서 0.64로 수렴하여 전 세계 결정 구조 상위 0.1% 수준의 초고해상도 기하 품질 확보.

PyMOL의 실측 RMSD 지표와 원자 카운트는 HanjariFold의 3단계 교차 정제 파이프라인이 1980년대 레거시 실험 구조의 기하 결함을 현대적 기준에 부합하도록 교정해 내는 표준 파이프라인임을 보여줍니다.

[태그:] 2LZM

[Benchmark Summary] 리파인먼트인가, 리모델링인가? 3종 대표 실험 구조로 입증한 골격 보존의 한계치

1. 대표 3종 실험 구조 벤치마크 총괄 요약

2. 핵심 시사점: 일관된 기하학적 안정성

[구조 정밀화] 2LZM 벤치마크: 공식 wwPDB 리포트 대조를 통해 검증한 T4 라이소자임 3단계 교차 정제 (MolProbity 0.64 달성)