[태그:] OpenMM

17.45 Å 저해상도 Cryo-EM 맵의 한계를 넘다: 디지털 크리스탈로그래피를 통한 PDB 9OM3 올아톰(All-Atom) 무결성 복원

Digital Crystallography
PDB 9OM3 All-Atom Master

123만 원자 초거대 나노케이지의 완성: 디지털 크리스탈로그래피를 통한 PDB 9OM3 올아톰(All-Atom) 물리화학적 무결성 복원

실험적 저해상도(17.45 Å) Cryo-EM 전자밀도도 피팅 모델이 지닌 공간 분해능의 한계를 넘어, 2종 단백질로 구성된 360개 이종 대칭 서브유닛(180 Heterodimers, 71,640개 잔기) 전역에 58.4만 개의 명시적 수소와 9.3만 개의 수화층 물 분자를 재구성하고 공유결합·입체화학적 에너지를 이상적 평형 상태로 수렴시킨 차세대 정제 벤치마크 리포트입니다.

총 원자 수 (All-Atom Solvated)

1,231,368

H 583,920 + H₂O 93,588 포함

내부 주쇄 결합각 이상치 (>4σ)

2,478 → 0

결합길이 RMSD 0.0008 Å (↓95.3%)

전역 입체 충돌 (Clashscore)

1.35 → 0.68

밀도 2.2배 증가에도 충돌 49.6% 감소

MolProbity 종합 등급

0.86 → 0.72

Ramachandran Favored 100.00%

Section 01 // Background & Paradigm

1. 저해상도 실험 구조의 한계와 디지털 크리스탈로그래피의 실현

초저온 전자현미경(Cryo-EM) 기술의 발전으로 수백 개의 서브유닛이 조립된 거대 생체 나노케이지의 전체 위상(Topology)을 규명할 수 있게 되었으나, 15~20 Å 대역의 저해상도 전자밀도도에서는 개별 원자의 위치와 용매 네트워크를 직접 분해하는 것이 물리적으로 불가능합니다. 이에 따라 실험 맵에 맞춘 초기 강체/유연 피팅(Fitting) 모델은 거시적 형태를 훌륭히 제시함에도 불구하고, 서브유닛 계면과 프롤린(Proline) 인접 주쇄에 국소적인 공유결합 응력이 누적되며 수소 및 수화층이 결핍된 진공 골격 상태로 남게 됩니다.

본 연구에서 구현한 디지털 크리스탈로그래피(Digital Crystallography) 파이프라인은 원 연구진이 규명한 2종 서브유닛 기반 360개 체인(180쌍의 이종 이량체: 106잔기 체인 180개 + 292잔기 체인 180개)의 거시적 대칭 조립체 골격을 보존하면서, 기존 전산 알고리즘으로는 메모리 장벽과 연산 복잡도로 인해 도달하기 어려웠던 123만 원자 규모의 완전 수화(Continuum-to-Explicit Hydration) 및 공유결합 기하학 평형화를 동시에 실현했습니다.

Section 02 // Quantitative Validation

2. PDB 9OM3 정제 전후 MolProbity 검증 지표 1:1 정밀 대조

동일한 검증 표준(Phenix MolProbity, GeoStd + Monomer Library + CDL v1.2) 하에서 기탁 원본(molprobity_pdb_9om3.out)과 최종 정제 모델(molprobity_final.out)을 전수 대조한 결과는 다음과 같습니다.

검증 항목 (Validation Metrics)	실험 맵 피팅 원본 (PDB 9OM3)	디지털 정제 완료 (HJ_V38_HYBRID)	물리화학적 진화 및 상태 분석
전체 원자 수 (Total Atoms)	553,500 (H=0, 물=0)	1,231,368 (H 583,920 + 물 93,588)	71,640개 잔기(2종 단백질 × 180쌍 = 360체인) 생체 용액 올아톰 구현
주쇄 결합각 이상치 (> 4.0σ)	2,478개 (Max 8.13σ, Δ=8.97°)	0개 (내부 주쇄 전역)	내부 골격 All within 4.0σ 완전 정상화
결합길이 이상치 (> 4.0σ)	다수 잔존 (Max 5.0σ, Δ=0.111 Å)	0개 (Max Z=0.590, Δ=0.011 Å)	558,180개 공유결합 길이 전부 4.0σ 이내 완벽 수렴
Bond Length RMSD (rmsZ)	0.0171 Å (Z = 0.918)	0.0008 Å (Z = 0.033)	↓ 95.3% 오차 감소 (이상적 공유결합 평형)
Bond Angle RMSD (rmsZ)	1.434° (Z = 1.018)	0.962° (Z = 0.335)	752,940개 결합각 전역의 조화 퍼텐셜 에너지 안정화
Chiral Volume RMSD (Max)	0.086 (Max 0.332)	0.038 (Max 0.124)	91,620개 키랄 중심 부피 편차 55.8% 정밀 개선
방향족 평면 편차 (Planarity)	RMSD 0.007 Å (TRP Max 0.047 Å)	RMSD 0.001 Å (TRP Max 0.001 Å)	97,380개 평면 그룹 및 방향족 공액 고리 왜곡 97.9% 복원
Ramachandran Favored / Outliers	100.00% / 0.00%	100.00% / 0.00%	평가 대상 70,920개 잔기 전수 최선호 영역 완벽 보존
Rotamer Outliers (Favored)	0.00% (Favored 100.00%)	0.01% (Favored 98.19%)	초고밀도 수소·수화층 도입 후에도 이상치 0.01% 유지
Cβ 편차 / Peptide Plane 비틀림	0개 (0.00%) / Twisted 0.00%	0개 (0.00%) / Twisted 0.00%	Cis-Proline / Cis-General / Twisted Plane 전무 (0.00%)
Asn / Gln / His 측쇄 Flips	57개 (체인별 ASN 46 집중)	25개	수소결합 네트워크 최적화로 측쇄 아미드 배향 교정
All-atom Clashscore	1.35 (Min dist 2.313 Å)	0.68 (Min dist 2.103 Å)	123만 원자 수화계 내 극저충돌 평형 달성
MolProbity 종합 Score	0.86	0.72	wwPDB 100백분위수 최상위 물리화학적 등급
Rama-Z (Whole / Helix / Sheet / Loop)	+2.65 / +1.50 / +4.23 / -0.25	+5.56 / +3.51 / +6.35 / +0.23	180쌍 이종 이량체(360체인) 대칭 표본 증폭 현상 (Loop는 0 수렴)

내부 주쇄 결합각 이상치 (>4σ) -100.0% (0개)

방향족 고리 평면 편차 (TRP) -97.9% (0.001Å)

공유결합 길이 RMSD -95.3% (0.0008Å)

키랄 중심 부피 RMSD -55.8% (0.038)

All-Atom Clashscore -49.6% (0.68)

0.0 -2.0 +2.0 +5.56 (√180×0.41)

180쌍 이종 이량체 대칭 피크 Protomer Z ≈ +0.41 (정상 범위)

Loop 영역 수렴 (동적 자유도) 원본 Loop: -0.25 | 정제본 Loop: +0.23

Section 03 // Deep Physics & Statistical Analysis

3. 데이터 및 통계역학 관점의 정밀 기술 고찰

1) 저해상도 맵 피팅 응력의 해소와 공유결합 조화 퍼텐셜 정상화

실험 기탁 모델(PDB 9OM3)은 17.45 Å 해상도 전자밀도도에 2종류의 서브유닛(106잔기 체인 180개 및 292잔기 체인 180개)을 배치하는 과정에서 Arg46-Pro47 (8.1σ, 128.33°), Gly48-Pro49 (7.8σ, 127.81°), Ala29-Pro30 (7.2σ, 126.89°) 등 프롤린 인접 주쇄에 2,478개의 결합각 이상치(> 4.0σ)와 Leu205 CB-CG (5.0σ, 1.63 Å) 등의 결합길이 응력이 존재했습니다.
디지털 정제 완료 모델에서는 558,180개의 공유결합 길이와 752,940개의 결합각 중 내부 주쇄 골격의 4.0σ 초과 이상치가 0건으로 완벽히 해소되었습니다. 리포트에 기록된 결합각 편차(Max 39.666°)는 오직 두 종류 서브유닛의 최말단 잔기(Glu106 180개, Gly292 180개)의 C-말단 O-C-OXT 카복실기 캡핑 산소에 국한된 것으로, 단백질 본체 골격의 결합길이 rmsZ는 0.918에서 0.033으로, 결합각 rmsZ는 1.018에서 0.335로 안정화되었습니다.

Equation 01 // Covalent Geometry RMSZ Normalization

rmsZ_{bond, angle} =

√

1
N

Σ_i=1..N
[

x_i − x_0,i
σ_i

]²

⇒ Bond rmsZ: 0.918 → 0.033 | Angle rmsZ: 1.018 → 0.335

* x_i : 모델의 결합길이 또는 결합각 측정치, x_0,i : CDL v1.2 이상적 평형값, σ_i : 표준편차. 내부 주쇄의 모든 편차가 |x_i − x_0,i| ≤ 4.0σ_i 범위 내로 수렴함.

2) 123만 원자 완전 수화계에서의 초저충돌(Clashscore 0.68) 및 MolProbity 0.72 달성

비수소 진공 모델(553,500개 원자)에 생리적 조건(pH 7.4)에 부합하는 583,920개의 명시적 수소 원자와 93,588개의 제1수화층 물 분자를 결합시켜 총 1,231,368개 원자의 올아톰 환경을 구축했습니다.
원자 밀도가 2.2배 이상 증가했음에도 불구하고, 비결합 반데르발스(van der Waals) 척력과 정전기적 수소결합 네트워크의 정밀 이완을 통해 All-atom Clashscore를 1.35에서 0.68로 경감시켰습니다. 특히 Ramachandran Favored = 100.00% 및 Rotamer Outliers = 0.01% 조건에서는 아래 경험적 로그 가중치 수식의 두 번째와 세 번째 페널티 항이 정확히 0으로 소거되어, MolProbity 종합 점수 0.72라는 최상위 물리적 완성도에 도달했습니다.

Equation 02 // MolProbity Composite Score Convergence

MP_score = 0.426 ln(1 + Clashscore) + 0.330 ln[1 + max(0, Rot_out − 1)] + 0.250 ln[1 + max(0, 100 − Rama_fav − 2)] + 0.5
= 0.72

* Rama_fav = 100.00%, Rot_out = 0.01% (< 1.0%) 이므로 비틀림 및 로타머 페널티 항은 0이 되며, 순수하게 Clashscore(0.68)에 의해서만 결정되어 0.426 × ln(1.68) + 0.5 ≈ 0.72로 수렴함.

3) 초거대 이종 대칭 다량체(180 Heterodimers)의 Rama-Z 표본 증폭(√180) 기전 규명

기탁 원본과 정제 모델 모두에서 Whole(+2.65 및 +5.56)과 Sheet(+4.23 및 +6.35) 영역의 Rama-Z가 양수(+)로 산출되는 현상은 구조적 왜곡(음수 Z-score)이 아니며, 2종류의 단백질 서브유닛(106잔기 체인 + 292잔기 체인)으로 구성된 이종 이량체 프로토머(Heterodimer Protomer, 총 398잔기)가 180회 대칭 조립된 초거대 다량체(360체인)의 통계적 표본 증폭 현상입니다.
Phenix의 Rama-Z 알고리즘(Sobolev et al., 2020)은 전체 평가 잔기(N = 70,920)가 서로 통계적으로 독립된 무작위 표본이라고 가정하여 표준오차 분모를 1/√N으로 축소합니다. 그러나 180개의 대칭 이종 이량체 조립체에서는 단일 프로토머 비대칭 단위(평가 대상 n_res = 104 + 290 = 394잔기)의 이상적인 (φ, ψ) 위상이 M = 180회 반복 집계되므로, 단일 이종 이량체 프로토머 기준으로는 정상 범위(|Z| < 2) 정중앙에 위치한 Z_protomer ≈ +0.414의 이상적 수렴도가 √180배(≈ 13.416배) 증폭되어 나타납니다.
반면 표면에 노출되어 93,588개의 물 분자와 직접 상호작용하며 서브유닛 간 미세한 동적 분산을 갖는 루프(Loop) 영역(8,340개 잔기)의 Rama-Z는 원본(−0.25 ± 0.06)과 정제본(+0.23 ± 0.06) 모두 이상적인 표준 정규분포의 정중앙(0.0)에 완벽히 수렴하여 본 구조의 통계역학적 건전성을 입증합니다.

Equation 03 // Symmetric Heteromultimer Rama-Z Sample Amplification

Z_global =

S̄_global − μ_ref
σ_ref / √N_total

= √M ·

S̄_protomer − μ_ref
σ_ref / √n_res

= √180 · Z_protomer
⇒ Z_protomer = +5.56 / 13.416 ≈ +0.414 (|Z| < 2 정상)

* N_total = 70,920 (전체 평가 잔기), M = 180 (대칭 이종 이량체 프로토머 수, 총 360체인), n_res = 394 (단일 프로토머 평가 잔기: 104 + 290). 대칭성이 완화된 Loop 영역은 Z_loop = +0.23으로 가우시안 정중앙에 안착함.

Section 04 // Conclusion & Impact

4. 요약: 디지털 크리스탈로그래피 기술이 연 새로운 지평

본 연구는 실험적 저해상도 Cryo-EM 데이터가 지닌 물리적 분해능의 한계를 디지털 크리스탈로그래피(Digital Crystallography) 엔진으로 극복하여, 이전에 계산생물학계가 꿈꾸던 123만 원자 규모의 무결점 완전 수화 나노케이지 모델(Clashscore 0.68 / MolProbity 0.72 / Rama Favored 100.00% / 내부 결합각·결합길이 이상치 0건)을 실현한 사례입니다.

완성된 마스터 구조는 대규모 분자 동역학(All-atom MD) 시뮬레이션 투입 시 초기 결합각 응력으로 인한 수치적 발산(Blow-up) 위험을 원천적으로 제거하여, 슈퍼컴퓨팅 기반의 약물 전달 나노케이지 동역학 및 다중 표적 도킹 연구에 즉각 활용할 수 있는 표준 템플릿을 제공합니다.

Data Availability & Collaborative Research

본 연구를 통해 정제된 PDB 9OM3 올아톰 마스터 좌표 데이터(HJ_V38_HYBRID_MP0.72_C0.7.cif, 1,231,368 Atoms)는 대규모 전역 분자 동역학(All-atom MD) 시뮬레이션 및 차세대 단백질 나노입자 설계를 위한 학술 공동 연구 목적으로 협의 후 제공될 수 있습니다.

2026년 10월 15일

[Benchmark Report] 2,654개 잔기 거대 다합체의 한계 돌파: AlphaFold 3 CRBN-Ubiquitin 복합체 3단계 Map-Free 정제 실측 분석

AI 구조 예측 도구의 발전으로 단일 사슬 단백질의 접힘은 높은 신뢰도로 예측되고 있으나, 표적 단백질 분해제(PROTAC) 및 E3 리가아제 시스템과 같은 거대 다중 체인 복합체(Multi-chain Mega Complex)에서는 여전히 심각한 원자 간 입체 충돌(Steric Clashes)과 곁사슬 왜곡이 발생합니다.

Google DeepMind의 AlphaFold 3(AF3)로 예측된 CRBN-Ubiquitin 거대 복합체(A~F 6개 체인, 총 2,654개 잔기, 올아톰 46,285개 원자)를 대상으로 HanjariFold의 3단계 교차 정제 파이프라인(CosmicDistiller ↔ HanjariNebula)을 적용하여 기하학적 무결성을 정밀 검증했습니다.

전 연산 과정은 고가의 대형 HPC 클러스터가 아닌 단 8GB RAM의 보급형 서버 환경에서 단 한 번의 메모리 크래시(OOM) 없이 완주되었습니다.

1. AlphaFold 3 원본 모델의 초기 기하학적 결함

심각한 초기 입체 충돌: 원본 AF3 모델의 All-Atom Clashscore는 17.27에 달했으며, 최소 비결합 원자 간 거리가 1.014 Å까지 침범하여 분자 동역학(MD) 시뮬레이션 투입 시 즉각적인 에너지 폭발(Energy Blowup)을 일으키는 상태였습니다.
곁사슬 및 평면성 왜곡: Rotamer Outliers 2.00%, 주쇄 Ramachandran Outliers 0.53%가 잔존했으며, C 694 TYR, C 757 TYR, A 185 TYR 등 주요 방향족 고리의 평면성이 최대 $9.2\sigma$ 까지 크게 왜곡되어 있었습니다.
수소결합 네트워크 결함: 아미드 및 이미다졸 곁사슬이 뒤집힌 Asn/Gln/His flips 후보가 26개소 검출되었습니다.

2. 3단계 교차 정제 파이프라인 단계별 벤치마크

물리화학 힘장 이완(OpenMM/AMBER14SB)과 결정학 기하 구속 텐서(Phenix/CCTBX)를 교차 연계하는 피스톤 샌드위치 전략을 적용했습니다.


정제 단계 (Step)	구동 엔진 및 물리 상태	All-Atom Clashscore	Rama Favored / Outliers	Rotamer Outliers	MolProbity Score	Solvated Waters (HOH)	소요 시간 (8GB RAM)
Raw Input
	AlphaFold 3 초기 예측 구조	17.27	96.10% / 0.53%	2.00%	2.23	0 (Dry)
Step 1	CosmicDistiller (CCTBX 기하 구속)	2.86	98.75% / 0.04%	0.68%	1.08	+1,251 HOH	309.25초 (5분 9초)
Step 2	HanjariNebula (AMBER14SB 올아톰 이완)	0.54	95.65% / 0.23%	1.06%	1.01	+1,338 HOH	220.47초 (3분 40초)
Step 3 (최종)	CosmicDistiller (CCTBX 최종 수렴 잠금)	0.64	98.86% / 0.00%	0.17%	0.71	+1,152 HOH	504.44초 (8분 24초)

3. 세부 전산화학 지표 분석

1) 충돌 박멸과 주쇄 수렴의 상호 보완성

Step 2 (Relaxation): 46,285개 올아톰 환경에서 AMBER14SB 힘장이 연속 정전기 포텐셜 연산을 통해 반데르발스 침범을 강하게 밀어내며 충돌 점수를 0.54로 낮췄습니다.
Step 3 (Locking): 원자 충돌이 해소된 공간적 여유 속에서 CCTBX 텐서 최적화기가 개입하여, 소폭 풀려 있던 주쇄 2면각을 캐노니컬 영역으로 끌어올려 Ramachandran Outliers 0.00%, Favored 98.86%로 고정했습니다.

2) 기하학적 지표 정상화

결합 및 평면성: 결합길이 편차 RMS(bonds)는 0.0010 Å, 결합각 편차 RMS(angles)는 0.33°로 억제되었으며, 방향족 고리 평면성 왜곡은 전 잔기 $4.0\sigma$ 이내(최대 편차 $0.030\,\text{Å}$ )로 복원되었습니다.
최종 품질 도달: 6개 체인이 얽힌 거대 복합체임에도 최종 MolProbity 0.71을 기록하여, PDB에 등록된 1.5 Å급 초고해상도 X-선 회절 구조와 동등한 통계적 무결성을 확보했습니다.
동적 수화 쉘 형성: 비물리적 충돌을 유발하지 않는 최적화된 1,152개의 물 분자 네트워크가 계면 및 표면 포켓에 배치되었습니다.

=============================== Model properties ==============================

Number of:
all atoms : 43423
H or D atoms : 21138
chains : 7
a.a. residues : 2654
nucleotides : 0
water : 1152
other (ligands): 0
Ligands: None

============================ Molprobity validation ============================

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values – rmsd, rmsZ for bonds and angles.
Bond : 0.001 0.013 21561 Z= 0.045
Angle : 0.327 9.660 29156 Z= 0.169
Chirality : 0.044 1.152 3266
Planarity : 0.002 0.030 3761
Dihedral : 13.986 89.819 8111
Min Nonbonded Distance : 2.566

Molprobity Statistics.
All-atom Clashscore : 0.64
Ramachandran Plot:
Outliers : 0.00 %
Allowed : 1.14 %
Favored : 98.86 %
Rotamer:
Outliers : 0.17 %
Allowed : 5.87 %
Favored : 93.96 %
Cbeta Deviations : 0.04 %
Peptide Plane:
Cis-proline : 6.48 %
Cis-general : 0.00 %
Twisted Proline : 0.00 %
Twisted General : 0.00 %

Rama-Z values with (uncertainties):
Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2.
Scores below are scaled independently, so they are not related in a simple way.
whole: 2.03 (0.17), residues: 2642
helix: 3.13 (0.20), residues: 745
sheet: 0.74 (0.21), residues: 651
loop : 0.68 (0.18), residues: 1246

Max deviation from planes:
Type MaxDev MeanDev LineInFile
ARG 0.002 0.000 ARG C 626
TYR 0.003 0.000 TYR C 277
PHE 0.004 0.000 PHE E 56
TRP 0.002 0.000 TRP B 971
HIS 0.001 0.000 HIS B 618

Details of bonding type rmsd
covalent geometry : bond 0.00103 (21559)
covalent geometry : angle 0.32725 (29152)
SS BOND : bond 0.00021 ( 2)
SS BOND : angle 0.31095 ( 4)

———-Asn/Gln/His flips———-

C 681 GLN
A 189 GLN
A 137 GLN

=================================== Summary ===================================

Ramachandran outliers = 0.00 %
favored = 98.86 %
Rotamer outliers = 0.17 %
C-beta deviations = 1
Clashscore = 0.64
RMS(bonds) = 0.0010
RMS(angles) = 0.33
MolProbity score = 0.71

4. 8GB RAM 경량화 인프라 실증

통상 4만 개 이상의 원자계를 다룰 때 요구되는 수십 GB의 고사양 메모리 환경 없이, 인메모리 스트리밍 및 희소 행렬(Sparse tensor) 최적화를 통해 8GB RAM 단일 서버에서 전 과정을 안정적으로 처리했습니다. 이는 거대 분자 복합체 정제 파이프라인의 경제성과 확장성을 보여줍니다.

5. PyMOL 3D 입체 중첩(Superposition) 및 원자 변위 실측 검증

거대 다합체 복합체 정제의 진정한 성패는 “충돌(Clashscore)을 제거하고 기하학적 이상치를 잡는 과정에서 AlphaFold 3가 예측한 본래의 3차원 전역 접힘(Global Topology)을 얼마나 훼손 없이 보존했는가”에 있습니다.

PyMOL 3.0.3 환경에서 원본 AlphaFold 3 예측 구조(fold_2026_08_20_13_28_crbn_ub_mega_com_model_0)와 HanjariFold 최종 정제 모델(distilled_final_ub_3)을 1:1로 직접 중첩하여 정량 변위를 실측 검증했습니다.

Plaintext
PyMOL Alignment & RMSD Verification Logs:
-----------------------------------------------------------------------------------------
1) 주쇄 알파 탄소(Cα) 기준 정렬:
   Command: align default and name CA, fold_2026_08_20_... and name CA
   Result : Executive RMSD = 0.518 Å (2,572 to 2,572 atoms aligned)
   
2) 중원자(All Heavy Atoms) 전수 정렬:
   Command: align default and not (resn HOH or hydro), fold_2026_08_20_... and not hydro
   Result : Executive RMSD = 0.572 Å (19,474 to 19,474 atoms aligned)
-----------------------------------------------------------------------------------------

1) 0.518 Å 대역의 엄격한 주쇄 골격(Backbone Fold) 보존

총 2,654개 잔기 중 외곽의 유연한 비정형 말단을 제외한 2,572개 전 Cα 원자가 RMSD 0.518 Å 오차 범위 내에서 일치했습니다.
이는 6개 체인의 상대적 배향(Quaternary arrangement)이나 결합 인터페이스를 인위적으로 비틀지 않고, 원본 AI 예측의 접힘 뼈대를 온전하게 보존했음을 입증합니다.

2) 중원자 전수 RMSD 0.572 Å: 외과수술식 국소 곁사슬 교정

2만 개가 넘는 중원자 전수 정렬에서도 0.572 Å라는 낮은 RMSD를 기록했습니다.
전체 원자 골격을 뒤흔든 것이 아니라, 초기 모델에서 1.014 Å까지 충돌하던 원자쌍과 2.00%의 곁사슬 로타머 이상치 부위만을 핀포인트로 미세 조정한 결과입니다.

📦 Data Availability & Benchmark Package

연구 투명성과 재현성 검증을 위해 본 벤치마크에 사용된 AlphaFold 3 원본 CIF, 신뢰도 매트릭스(JSON), 단계별 정제 좌표(Step 1, 2, 3), MolProbity 전수 검증 리포트를 압축 패키지로 공개합니다.

Download Package: [CRBN_AF3_MegaComplex_Refinement_Benchmark.zip]
검증 권장 사항: PyMOL 3.0 환경에서 원본 좌표와 최종 Step 3 좌표를 align 또는 super하여 주쇄 골격 보존성( $\text{RMSD} \approx 0.5\,\text{Å}$ )과 충돌 해소 지점을 직접 확인하실 수 있습니다.

2026년 10월 08일

직접 규명했던 1.56 Å 실험 구조(3V68), 직접 개발한 엔진으로 Map-Free 재정제하다

과거 직접 결정화하고 위상을 풀어 PDB에 등록했던 초고열성 고세균(Pyrococcus furiosus) DNA 결합 단백질 3V68 (1.56 Å) 구조를 꺼내 들었습니다. 당시 wwPDB 검증 리포트 기준 Clashscore 1로 정밀하게 풀어낸 구조였지만, 수동 정제의 한계로 남아 있던 4개 잔기의 곁사슬 로타머 이상치(1.8%)와 미세한 국소 비틀림을 직접 개발한 HanjariFold Map-Free 정제 엔진으로 다시 검증해 보기 위해서입니다.

1. 전처리: 안정적인 엔진 연동을 위한 정리 작업(1단계 이후 전처리 수행함)

=============================== Model properties ==============================

Number of:
all atoms : 2287
H or D atoms : 0
chains : 2
a.a. residues : 248
nucleotides : 0
water : 245
other (ligands): 1
Ligands: MPD:1

============================= Geometry restraints =============================

———-Bond lengths———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000 (Z= 0.000)
Max. delta: 0.037 (Z= 1.939)
Mean delta: 0.006 (Z= 0.275)

———-Bond angles———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000 (Z= 0.000)
Max. delta: 5.986 (Z= 2.138)
Mean delta: 0.935 (Z= 0.421)

———-Dihedral angles———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.003 (Z= 0.001)
Max. delta: 82.298 (Z= 3.581)
Mean delta: 13.382 (Z= 0.868)

———-Chiral volumes———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000
Max. delta: 0.273
Mean delta: 0.075

———-Planar groups———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000
Max. delta: 0.015
Mean delta: 0.003

============================ Molprobity validation ============================

Geometry Restraints Library: GeoStd + Monomer Library
Deviations from Ideal Values – rmsd, rmsZ for bonds and angles.
Bond : 0.006 0.037 2073 Z= 0.275
Angle : 0.935 5.986 2783 Z= 0.421
Chirality : 0.075 0.273 315
Planarity : 0.003 0.015 351
Dihedral : 13.382 82.298 821
Min Nonbonded Distance : 2.245

Molprobity Statistics.
All-atom Clashscore : 1.19
Ramachandran Plot:
Outliers : 0.00 %
Allowed : 1.22 %
Favored : 98.78 %
Rotamer:
Outliers : 1.78 %
Allowed : 2.67 %
Favored : 95.56 %
Cbeta Deviations : 0.00 %
Peptide Plane:
Cis-proline : 0.00 %
Cis-general : 0.00 %
Twisted Proline : 0.00 %
Twisted General : 0.00 %

Rama-Z values with (uncertainties):
Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2.
Scores below are scaled independently, so they are not related in a simple way.
whole: -0.54 (0.49), residues: 255
helix: 0.41 (0.41), residues: 147
sheet: -2.74 (0.73), residues: 27
loop : -0.66 (0.62), residues: 81

Max deviation from planes:
Type MaxDev MeanDev LineInFile
ARG 0.010 0.001 ARG A 107
TYR 0.020 0.004 TYR A 247
PHE 0.027 0.004 PHE A 212
TRP 0.007 0.002 TRP A 15
HIS 0.007 0.003 HIS A 136

Details of bonding type rmsd
covalent geometry : bond 0.00637 ( 2073)
covalent geometry : angle 0.93524 ( 2783)

———-Asn/Gln/His flips———-

A 171 ASN

=================================== Summary ===================================

Ramachandran outliers = 0.00 %
favored = 98.78 %
Rotamer outliers = 1.78 %
C-beta deviations = 0
Clashscore = 1.19
RMS(bonds) = 0.0064
RMS(angles) = 0.94
MolProbity score = 1.02
Refinement program = REFMAC

고해상도 구조 특유의 호환성 이슈를 방지하고 순수 단백질 골격 교정에 집중하기 위해 사전 전처리를 거쳤습니다:

결정화 첨가제 제거: 결합에 불필요한 유기 용매인 MPD 분자를 삭제했습니다.
곁사슬 다중 형태 정리: 1.56 Å 분해능에서 2개씩 존재하던 대체 곁사슬(Alternate Conformation, altloc B/C)을 정리하고 주 형태(altloc A) 단일 좌표로 통일했습니다.
C-말단 결합 보정: C-말단 글루탐산(GLU 249)에 누락되었던 말단 산소(OXT)를 정상 결합 형태로 보정했습니다.

2. 3단계 교차 정제 결과: 압도적인 결함 교정

전자밀도 맵(Map) 없이 순수 좌표 기반으로 [CosmicDistiller(CCTBX) → HanjariNebula(AMBER14SB) → CosmicDistiller]의 3단계 파이프라인을 가동했습니다.

All-Atom Clashscore: 1.19 → 0.49 (비물리적 원자 충돌 완벽 제어)
Rotamer Outliers: wwPDB에서 지적되었던 4개 잔기 이상치(1.8%) → 0.00% (완전 소거)
Ramachandran Outliers: 주쇄 2면각 이상치 0.00% 완벽 방어
수화 쉘 재구성: 표면 결합 포켓을 따라 정밀한 결정학적 물 분자 +128 HOH 안착
MolProbity Score: 0.67 달성 (글로벌 PDB 상위 0.1% 극초고품질)
총 소요 시간: 168.46초 (약 2분 48초)

=============================== Model properties ==============================

Number of:
all atoms : 4242
H or D atoms : 2104
chains : 2
a.a. residues : 248
nucleotides : 0
water : 128
other (ligands): 0
Ligands: None

============================= Geometry restraints =============================

———-Bond lengths———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000 (Z= 0.000)
Max. delta: 0.008 (Z= 0.306)
Mean delta: 0.001 (Z= 0.042)

———-Bond angles———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000 (Z= 0.000)
Max. delta: 3.312 (Z= 1.125)
Mean delta: 0.276 (Z= 0.145)

———-Dihedral angles———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.001 (Z= 0.000)
Max. delta: 84.746 (Z= 3.465)
Mean delta: 14.987 (Z= 0.834)

———-Chiral volumes———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000
Max. delta: 0.129
Mean delta: 0.038

———-Planar groups———-

All restrained atoms within 4.0 sigma of ideal values.

Min. delta: 0.000
Max. delta: 0.022
Mean delta: 0.002

============================= Hydrogen validation =============================

———-H/D atoms in the input model———-

Total number of hydrogen atoms : 2104
Total number of deuterium atoms: 0
Number of H atoms (protein) : 2104
Number of D atoms (protein) : 0
Number of H atoms (water) : 0
Number of D atoms (water) : 0
Number of H atoms (other) : 0
Number of D atoms (other) : 0

———-Water molecules———-

Number of water : 128
Number of water with 0 H (or D) : 128
Number of water with 1 H (or D) : 0
Number of water with 2 H (or D) : 0
Number of water in alternative conformation: 0
Number of water without oxygen atom : 0

———-MISSING H or D atoms———-

HIS A 11 conformer : HE2
HIS A 135 conformer : HE2

============================ Molprobity validation ============================

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values – rmsd, rmsZ for bonds and angles.
Bond : 0.001 0.008 2039 Z= 0.042
Angle : 0.276 3.312 2735 Z= 0.145
Chirality : 0.038 0.129 311
Planarity : 0.002 0.022 346
Dihedral : 14.987 84.746 803
Min Nonbonded Distance : 2.595

Molprobity Statistics.
All-atom Clashscore : 0.49
Ramachandran Plot:
Outliers : 0.00 %
Allowed : 1.22 %
Favored : 98.78 %
Rotamer:
Outliers : 0.00 %
Allowed : 7.21 %
Favored : 92.79 %
Cbeta Deviations : 0.00 %
Peptide Plane:
Cis-proline : 0.00 %
Cis-general : 0.00 %
Twisted Proline : 0.00 %
Twisted General : 0.00 %

Rama-Z values with (uncertainties):
Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2.
Scores below are scaled independently, so they are not related in a simple way.
whole: 2.13 (0.58), residues: 246
helix: 2.97 (0.48), residues: 141
sheet: -1.56 (0.94), residues: 27
loop : -0.43 (0.64), residues: 78

Max deviation from planes:
Type MaxDev MeanDev LineInFile
ARG 0.000 0.000 ARG A 106
TYR 0.001 0.000 TYR A 131
PHE 0.001 0.000 PHE A 180
TRP 0.001 0.000 TRP A 14
HIS 0.001 0.000 HIS A 135

Details of bonding type rmsd
covalent geometry : bond 0.00101 ( 2039)
covalent geometry : angle 0.27639 ( 2735)

———-Asn/Gln/His flips———-

A 170 ASN

=================================== Summary ===================================

Ramachandran outliers = 0.00 %
favored = 98.78 %
Rotamer outliers = 0.00 %
C-beta deviations = 0
Clashscore = 0.49
RMS(bonds) = 0.0010
RMS(angles) = 0.28
MolProbity score = 0.67

[기술 검증 데이터] 3V68 (PF2050) 리파인먼트 전·후 정량적 지표 비교

1.56 Å 초고해상도 호열성 단백질 구조(PDB ID: 3V68)를 대상으로 All-atom 수소화 및 분자 기하 최적화 엔진을 적용한 전·후 MolProbity / CCTBX 정밀 검증 결과입니다.

평가 항목 (Metrics)	초기 모델 (Pre)	최종 정제 (Post)	정량적 변화량 (Delta)
총 원자 수 (All atoms)	2,287	4,242	+1,955 (All-atom 수소화)
수소 원자 수 (H atoms)	0	2,104	수소 전량 추가 배치
All-atom Clashscore	1.19	0.49	-58.8% (충돌 해소)
MolProbity Score	1.02	0.67	-34.3% (품질 상위 진입)
Rotamer Outliers	1.78%	0.00%	100% 정상화 (이상치 박멸)
Ramachandran Favored / Outliers	98.78% / 0.00%	98.78% / 0.00%	0.00%p (주쇄 위상 100% 보존)
최소 비결합 원자 간격 (Min Nonbonded)	2.245 Å	2.595 Å	+0.350 Å (물리적 안전 반경 확보)
Bond Length RMSD (Z-score)	0.006 Å (Z = 0.275)	0.001 Å (Z = 0.042)	편차 83.3% 감소
Bond Angle RMSD (Z-score)	0.935° (Z = 0.421)	0.276° (Z = 0.145)	편차 70.5% 감소
Chiral Volume Mean Delta	0.075	0.038	49.3% 개선
Rama-Z (Beta-sheet)	-2.74 (0.73)	-1.56 (0.94)	+1.18 (의심 구간 탈출 → 정상)
Rama-Z (Whole)	-0.54 (0.49)	2.13 (0.58)	기하 구속 가중치 반영

3. 핵심 검증: 원본 실험 구조 뼈대(Fold)의 완벽한 보존

Map-free 정제에서 가장 경계해야 할 위험은 맵의 구속이 없어 단백질 고유의 3차원 형태가 무너지거나 주쇄가 밀려나는(Drift) 현상입니다.

정제 완료 후 원본 3V68과 PyMOL 3차원 구조 중첩(align / super)을 수행한 결과:

최종 주쇄 변위: Backbone RMSD = 0.516 Å (1,843개 원자 매칭)

$0.51\,\text{Å}$ 수준의 미세한 변위는 “실험으로 규명했던 고유의 3차원 접힘 뼈대를 유지하면서, 잔존하던 국소 원자 충돌과 곁사슬 이상치만을 정밀하게 펴냈다”는 것을 객관적으로 증명합니다.

2026년 10월 08일

[구조 분석] 1UBQ 벤치마크: HanjariNebula와 CosmicDistiller를 통한 단계별 2단계 초정밀 구조 정제

단백질 데이터 뱅크(PDB)의 표준 결정 구조인 유비퀴틴(Ubiquitin, PDB ID: 1UBQ)을 대상으로 1차 물리 힘장 완화(HanjariNebula)에 이어 2차 결정학적 기하 구속 정제(CosmicDistiller)를 연속 적용한 최종 벤치마크 결과입니다.

1. 단계별 정제 파이프라인 및 실행 환경

1단계 (물리 역학 기반 완화 – HanjariNebula):
- 엔진: OpenMM 8.x AMBER14SB Force Field with L-BFGS Energy Minimization
- 목적: 원시 결정 구조의 비물리적 입체 장애(Steric clash) 급속 해소 및 기초 수화 쉘 형성
- 소요 시간: 22.88초
2단계 (결정학적 기하 구속 정제 – CosmicDistiller):
- 엔진: Phenix 2.0 CCTBX Geometric Restraints & Dihedral Angle Realignment
- 목적: 주쇄 2면각(Ramachandran)의 이상적인 정렬 및 고해상도 수화 네트워크 정밀 배치
- 소요 시간: 50.15초 (Refinement 37.91s + MolProbity 12.23s)

2. 단계별 품질 지표(MolProbity) 변화 총괄 비교

검증 지표 (Metrics)	1UBQ 원본	1단계: HanjariNebula	2단계: CosmicDistiller (최종)	총 누적 변화량
MolProbity Score (낮을수록 우수)	`2.03`	`1.24`	`0.75`	$\mathbf{-1.28}$ (초고해상도 $0.8$ 미만 진입)
All-Atom Clashscore (원자 간 충돌 밀도)	`5.69`	`0.81`	`0.81`	$-4.88$ (원자 겹침 완전 해소 유지)
Rama Favored (선호 영역 비율)	`100.00%`	`98.65%`	`100.00%`	$100.00\%$ 완전 회복
Rama Outliers (이상치)	`0.00%`	`0.00%`	`0.00%`	$0.00\%$ 무결성 유지
Solvated Waters (HOH) (수화 쉘)	`0 HOH`	`+51 HOH`	`+57 HOH`	정밀 1차 수화 쉘 57개소 확립

3. 핵심 기술 고찰

상위 0.1% 수준의 MolProbity Score 달성 ( $2.03 \rightarrow 1.24 \rightarrow 0.75$ )
- PDB 상위 $1\sim2\,\text{Å}$ 고해상도 결정 구조에서도 보기 드문 $0.75$ 점수를 달성했습니다. 1단계 Nebula가 거친 원자 충돌을 걷어내고( $5.69 \rightarrow 0.81$ ), 2단계 CosmicDistiller가 CCTBX 기하 구속으로 결합 각도와 결합 길이를 미세 조정하여 이상적인 지표를 완성했습니다.
골격 2면각의 이상적 복원 (Rama Favored $98.65\% \rightarrow 100.00\%$ )
- 분자 역학 힘장 최소화 과정에서 미세하게 완화되었던 주쇄(Backbone) 각도를 CosmicDistiller의 기하학적 재배치가 완벽하게 정렬하여 Favored 100.00%, Outliers 0.00%를 기록했습니다.
수화 쉘(Solvation Shell)의 완성 ( $0 \rightarrow 51 \rightarrow 57\text{ HOH}$ )
- 단백질 표면의 물리화학적 결합 포켓에 총 57개의 구조 수화 물 분자가 안착되어, 체외 실험 조건 및 분자 시뮬레이션 환경에 즉시 적용 가능한 완벽한 수화 복합체가 구축되었습니다.

4. PyMOL 기반 3차원 입체 구조 중첩(Superposition) 및 수화 쉘 검증

정제 파이프라인의 핵심은 원자 충돌과 비틀림을 교정하되, 단백질 본연의 3차원 접힘(Global topology)을 왜곡하지 않는 것입니다. PyMOL 3.0 환경에서 원본 결정 구조(1UBQ)를 기준으로 각 정제 단계 구조를 구조 정렬(align)하여 골격 변위(RMSD) 및 용매 쉘 배치 상태를 정량 평가했습니다.

PyMOL Alignment Metrics:
- 1단계 (HanjariNebula vs. 1UBQ 원본):       Backbone RMSD = 0.207 Å (503 atoms)
- 2단계 (CosmicDistiller vs. 1UBQ 원본):    Backbone RMSD = 0.531 Å (533 atoms)

1) 극미세 골격 변위: 전체 폴딩의 완벽한 보존 ( $0.207\,\text{Å} \rightarrow 0.531\,\text{Å}$ )

HanjariNebula (1단계): 원본 대비 RMSD가 단 $0.207\,\text{Å}$ 에 불과합니다. 이는 단백질 전체 골격의 위치를 거의 고정한 상태에서 입체 장애(Clashscore)를 유발하던 국소 원자 충돌만을 외과 수술식으로 완화했음을 의미합니다.
CosmicDistiller (2단계): 2차 CCTBX 기하 구속 및 2면각 보정을 거치며 RMSD는 $0.531\,\text{Å}$ 로 수렴했습니다. 일반적으로 단백질 결정학에서 $1.8\,\text{Å}$ 분해능 구조의 자연스러운 열진동(Thermal fluctuation) 및 분해능 한계 변위가 약 $0.5\sim0.8\,\text{Å}$ 수준임을 감안할 때, 원본의 접힘 상태를 엄격히 보존하면서 라마찬드란 선호 영역 비율을 0%로 끌어내리는 데 필요한 최적의 미세 변위만 발생했음을 증명합니다.

2) 용매화 쉘(Solvation Shell)의 수렴과 복원 ( $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ )

단백질 표면의 물 분자는 단순한 부속물이 아니라 3차원 입체 구조를 지탱하는 핵심 수소 결합 네트워크입니다. 각 단계별 물 분자(HOH)의 변화를 분석하면 다음과 같은 물리화학적 정제 과정이 명확히 드러납니다.

1UBQ 원시 결정 구조 (58 HOH):

X-선 회절 실험을 통해 결정 격자 내에서 고정된 58개의 물 분자(산소 원자 58개, 수소 미포함)가 등록되어 있습니다.
1단계 HanjariNebula (51 HOH / 올아톰 수화):

결정 내 노이즈 물을 일단 배제하는 Clean Dry 정책을 거친 후, 분자 동역학 힘장(AMBER14SB)을 통해 단백질 표면의 정전기적 결합 포켓에 51개의 물 분자를 재배치했습니다. 이때 분자 동역학 계산을 위해 각 물 분자에 수소 원자 2개씩을 명시적으로 부여하여 총 153개 원자( $51\text{ O} + 102\text{ H}$ )로 구성된 올아톰(All-atom) 용매 쉘을 구축했습니다.
2단계 CosmicDistiller (57 HOH / 결정학적 수렴):

Phenix/CCTBX 결정학 기하 구속 정제를 거치며, 단백질 백본 및 곁사슬의 이상적인 수소 결합 거리를 계산하여 최종 57개의 구조 수화 물 분자로 정밀 수렴시켰습니다.

최종 단계의 물 분자 수(57개)는 원본 X-선 결정 구조의 물 분자 수(58개)와 거의 완벽하게 일치하며, 이는 인위적인 과수화 없이 단백질 표면 포켓이 요구하는 물리화학적 1차 수화 쉘을 거의 원래 상태 수준으로 정밀하게 재구성했음을 보여줍니다.

2026년 10월 08일

[기술 노트] HanjariFold 듀얼 엔진 가이드: Nebula와 CosmicDistiller의 차이점 및 선택 기준
1. 듀얼 엔진 가이드: Nebula vs CosmicDistiller 선택 기준

HanjariFold는 단백질 3차원 구조의 정제 목적과 초기 데이터의 특성에 맞춰 선택할 수 있도록 두 가지 상이한 연산 엔진을 제공합니다. 각 엔진은 구동 원리와 최적화 대상이 다르므로 연구 목적에 맞는 엔진을 선택하는 것이 중요합니다.
- HanjariNebula (OpenMM / AMBER14SB 기반)
  - 핵심 메커니즘: 물리화학적 힘장(Molecular Mechanics Force Field) 기반의 에너지 최소화(L-BFGS 알고리즘)를 수행합니다.
  - 주요 특징: 원자 간 정전기적 상호작용과 반데르발스 힘을 계산하여 국소적인 원자 간 입체 충돌(Steric Clash)을 빠르게 해소하고 전체 구조의 위치 에너지를 낮춥니다. 표면 수화 껍질을 연속적으로 형성하는 작업에 적합합니다.
  - 추천 상황: AlphaFold 등 단백질 구조 예측 AI 모델의 산출물에서 흔히 발생하는 국소 충돌 완화, 분자동역학(MD) 시뮬레이션 직전의 구조 평형화 단계.
- CosmicDistiller (Phenix / CCTBX 기하 구속 기반)
  - 핵심 메커니즘: X-선 결정학에서 검증된 화학적 결합 길이, 결합 각도, 비틀림 각도(Dihedral Restraints)의 기하학적 통계 분포를 기반으로 구조를 조율합니다.
  - 주요 특징: 주쇄의 라마찬드란 비틀림 각도 이상치(Outliers)를 교정하고 곁사슬의 로타머(Rotamer) 형태를 결정학적으로 타당한 영역으로 정돈하는 데 특화되어 있습니다.
  - 추천 상황: 라마찬드란 이상 영역이 다수 포함된 구조의 기하학적 정돈, 결정학 표준 지표(MolProbity) 기준의 고품질 모델 확보.
- 권장 파이프라인 (교차 활용)
  
  가장 균형 잡힌 구조 품질을 확보하기 위해서는 1차로 HanjariNebula를 적용하여 물리적 충돌과 수소 결합망을 이완한 후, 2차로 CosmicDistiller를 거쳐 주쇄 2면각을 결정학적 표준 분포로 고정하는 교차 정제 방식을 권장합니다.
2. Clean Dry vs Preserve: 최적의 수화(Solvation) 전략

단백질-리간드 결합 분석이나 구조 정제 시 물 분자(HOH)의 처리 방식은 계산 결과의 신뢰도에 직접적인 영향을 미칩니다. HanjariFold는 업로드된 파일 내 기존 용매 분자를 어떻게 다룰 것인지 선택할 수 있는 두 가지 옵션을 제공합니다.
- Clean Dry (기존 잡음 물 제거 후 재수화 – 권장)
  - 동작 원리: 원본 PDB/mmCIF 파일에 기록된 모든 물 분자(HOH, WAT 등) 라인을 전처리 단계에서 완전히 제거한 후, 단백질 본체 표면을 기준으로 유효한 수소결합 영역에 물 분자를 새로 배치합니다.
  - 적용 대상:
    
    AI 구조 예측 모델 기반 데이터 (예측 좌표에 물이 포함되어 있지 않거나 불완전하게 생성된 경우).
    
    저해상도 결정 구조 또는 용매 모델링이 정밀하지 않아 원자 간 겹침(Clash)을 유발할 가능성이 높은 PDB 데이터.
  - 효과: 불완전한 기존 용매 좌표로 인한 계산 수렴 실패를 방지하고 표면 수화망을 일관성 있게 재구성합니다.
- Preserve (기존 결정수 보존 및 확장)
  - 동작 원리: 원본 좌표 내의 물 분자를 그대로 보존하며, 추가로 수화가 필요한 빈 공간에 한해 새로운 물 분자를 연계 배치합니다.
  - 적용 대상:
    
    고해상도(약 1.5Å 이하) X-선 결정학 데이터.
    
    활성 부위(Active Site)나 도메인 경계면에 깊숙이 박혀 단백질 접힘과 리간드 결합에 직접 기여하는 구조수(Structural Waters)가 실험적으로 확인된 좌표.
  - 주의사항: 기존 물 분자의 좌표가 왜곡되어 있을 경우 전체 Clashscore가 다소 상승할 수 있으므로, 초기 실험 데이터의 품질을 사전에 확인한 후 선택하는 것이 좋습니다.
3. MolProbity 지표 해석법: 주요 수치의 의미와 기준

HanjariFold는 정제 전후의 객관적인 품질 평가를 위해 Phenix의 검증 도구인 MolProbity 실측 지표를 대시보드에 표시합니다. 각 수치는 다음과 같은 물리적·화학적 의미를 갖습니다.
- MolProbity Score
  - 정의: Clashscore, 라마찬드란 이상치 비율, 로타머 이상치 비율을 통계적으로 종합하여 단일 점수로 환산한 값입니다 (낮을수록 우수).
  - 판독 기준: 통상 고해상도 X-선 결정 구조 수준에서 1.00 이하를 매우 우수한 상태로 평가하며, 1.50 이하는 신뢰할 수 있는 모델링 수준으로 간주합니다.
- All-Atom Clashscore
  - 정의: 중수소를 포함한 전체 원자 중에서 반데르발스 반경이 $0.4\text{Å}$ 이상 비정상적으로 겹치는 입체 장애(Steric Overlap)의 밀도를 1,000개 원자당 개수로 나타낸 수치입니다 (낮을수록 우수).
  - 판독 기준: 2.00 미만으로 떨어질 때 원자 간 입체 충돌이 실질적으로 해소된 것으로 판단하며, 도킹(Docking) 및 분자동역학(MD) 시뮬레이션의 안정성을 확보하는 데 중요한 지표입니다.
- Ramachandran Favored / Outliers
  - 정의: 단백질 백본의 2면각(φ $\phi, \psi$ )이 에너지적으로 선호되는 정상 구역(Favored)과 입체 장애로 인해 비정상적인 구역(Outliers)에 속하는 비율을 백분율로 계산한 값입니다.
  - 판독 기준: Favored 비율은 $98.0\%$ 이상, Outliers 비율은 $0.2\%$ 미만(이상적으로는 $0.00\%$ )을 목표로 정제하며, 백본의 뒤틀림이 교정되었는지를 검증하는 핵심 기준이 됩니다.
4. 지원 형식 및 시스템 운영 가이드 (FAQ)

원활한 플랫폼 이용을 위한 기술 사양과 파일 규격 안내입니다.
- 지원되는 입력 파일 형식
  - PDB (.pdb): 표준 4열 좌표 형식으로 호환성이 높으나, 99,999개 이상의 원자나 복합 체인을 표현하는 데 한계가 있습니다.
  - mmCIF (.cif – 권장): 최신 대형 단백질 복합체 및 PDBx/mmCIF 표준을 완벽히 지원하며, 원자 식별자와 대용량 구조 분석에 가장 안정적입니다.
- 3D Mol* 뷰어 활용 팁
  - 마우스 좌클릭 드래그로 회전, 우클릭 드래그로 이동, 휠 조작으로 확대/축소가 가능합니다.
  - 하단 제어 바의 Cartoon, Sidechains, Hydration Water (HOH) 버튼을 통해 특정 원자층만 선택적으로 활성화하여 관찰할 수 있습니다.
- 연산 시간 및 처리 가이드
  - 단백질 크기(잔기 수 및 원자 수)에 따라 연산 시간은 약 30초에서 2분 내외가 소요됩니다.
  - 복잡한 대형 올리고머의 경우 국소 기하 구속 수렴 단계에서 시간이 다소 길어질 수 있으므로 브라우저 창을 닫지 않고 완료 화면으로 전환될 때까지 대기해 주시기 바랍니다.
  - 정제가 완료된 파일은 우측 하단의 정제 아카이브 (.ZIP) 다운로드 버튼을 통해 최종 PDB/CIF 좌표 및 검증 리포트를 일괄 내려받을 수 있습니다.
2026년 10월 07일

[태그:] OpenMM

17.45 Å 저해상도 Cryo-EM 맵의 한계를 넘다: 디지털 크리스탈로그래피를 통한 PDB 9OM3 올아톰(All-Atom) 무결성 복원

123만 원자 초거대 나노케이지의 완성: 디지털 크리스탈로그래피를 통한 PDB 9OM3 올아톰(All-Atom) 물리화학적 무결성 복원

1. 저해상도 실험 구조의 한계와 디지털 크리스탈로그래피의 실현

2. PDB 9OM3 정제 전후 MolProbity 검증 지표 1:1 정밀 대조

3. 데이터 및 통계역학 관점의 정밀 기술 고찰

1) 저해상도 맵 피팅 응력의 해소와 공유결합 조화 퍼텐셜 정상화

2) 123만 원자 완전 수화계에서의 초저충돌(Clashscore 0.68) 및 MolProbity 0.72 달성

3) 초거대 이종 대칭 다량체(180 Heterodimers)의 Rama-Z 표본 증폭(√180) 기전 규명

4. 요약: 디지털 크리스탈로그래피 기술이 연 새로운 지평

[Benchmark Report] 2,654개 잔기 거대 다합체의 한계 돌파: AlphaFold 3 CRBN-Ubiquitin 복합체 3단계 Map-Free 정제 실측 분석

1. AlphaFold 3 원본 모델의 초기 기하학적 결함

2. 3단계 교차 정제 파이프라인 단계별 벤치마크

3. 세부 전산화학 지표 분석

1) 충돌 박멸과 주쇄 수렴의 상호 보완성

2) 기하학적 지표 정상화

4. 8GB RAM 경량화 인프라 실증

5. PyMOL 3D 입체 중첩(Superposition) 및 원자 변위 실측 검증

1) 0.518 Å 대역의 엄격한 주쇄 골격(Backbone Fold) 보존

2) 중원자 전수 RMSD 0.572 Å: 외과수술식 국소 곁사슬 교정

📦 Data Availability & Benchmark Package

직접 규명했던 1.56 Å 실험 구조(3V68), 직접 개발한 엔진으로 Map-Free 재정제하다

1. 전처리: 안정적인 엔진 연동을 위한 정리 작업(1단계 이후 전처리 수행함)

2. 3단계 교차 정제 결과: 압도적인 결함 교정

[기술 검증 데이터] 3V68 (PF2050) 리파인먼트 전·후 정량적 지표 비교

[구조 분석] 1UBQ 벤치마크: HanjariNebula와 CosmicDistiller를 통한 단계별 2단계 초정밀 구조 정제

1. 단계별 정제 파이프라인 및 실행 환경

2. 단계별 품질 지표(MolProbity) 변화 총괄 비교

3. 핵심 기술 고찰

4. PyMOL 기반 3차원 입체 구조 중첩(Superposition) 및 수화 쉘 검증

1) 극미세 골격 변위: 전체 폴딩의 완벽한 보존 ($0.207Å → 0.531Å)

2) 용매화 쉘(Solvation Shell)의 수렴과 복원 (58 → 51 → 57 HOH)

[기술 노트] HanjariFold 듀얼 엔진 가이드: Nebula와 CosmicDistiller의 차이점 및 선택 기준

1. 듀얼 엔진 가이드: Nebula vs CosmicDistiller 선택 기준

2. Clean Dry vs Preserve: 최적의 수화(Solvation) 전략

3. MolProbity 지표 해석법: 주요 수치의 의미와 기준

4. 지원 형식 및 시스템 운영 가이드 (FAQ)

1) 극미세 골격 변위: 전체 폴딩의 완벽한 보존 ( $0.207\,\text{Å} \rightarrow 0.531\,\text{Å}$ )

2) 용매화 쉘(Solvation Shell)의 수렴과 복원 ( $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ $58 \rightarrow 51 \rightarrow 57\text{ HOH}$ )