[카테고리:] 연구노트

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일

[기술 리포트] 58만 원자급 대칭성 나노케이지(PDB 9OM3) 1회차 지오메트리 정제 결과 분석

개요

거대분자 복합체(Macromolecular assembly) 모델링에서는 전체적인 형태 피팅과 국소 결합 기하학(Local covalent geometry) 간의 균형을 맞추는 것이 핵심 과제입니다.

본 리포트는 극저해상도(17.45 Å) Cryo-EM 전자밀도 맵을 기반으로 구축된 55만 원자 규모의 대칭성 나노케이지 구조체(PDB ID: 9OM3)를 대상으로, 수화 셸(Hydration shell) 배치 및 1회차 지오메트리 최소화(Single-pass geometry minimization)를 수행한 전후 검증 지표를 객관적으로 비교·정리한 기술 보고입니다.

1. 대상 시스템 및 계산 환경

대상 구조: PDB ID 9OM3 (360개 체인, 71,640개 잔기, 비수소 원자 553,500개)
입력 데이터 상태: 저해상도 피팅 특성상 주쇄 및 프롤린 인접 잔기 구간에 다수의 결합각 응력이 누적되어 있는 상태
수화 셸 구축: cKDTree 기반 격자 탐색(4.5 Å 간격)을 통해 계면 완충용 물 분자(HOH) 33,414개 동적 배치 (총 4개 독립 체인 분할 수용)
연산 규모: 총 586,914개 원자 (단백질 553,500 + 수화 셸 33,414)
최소화 파이프라인: CCTBX/Phenix 환경 기반 지오메트리 최소화 (내부 Riding Hydrogen 모델 연동, Hybrid-36 전역 인덱싱 적용)
소요 시간: 단일 워크스테이션 노드 기준 총 2,938.59초 (약 48분 58초)

Validation Benchmark

Target: PDB 9OM3

2. 정제 전후 검증 지표 정량 비교

동일 검증 도구(Phenix MolProbity)를 통해 산출된 원본 리포트(molprobity_pdb_2.out) 및 1회차 정제 모델(distilled_final.cif) 수치 대조

검증 항목	정제 전 (PDB 9OM3 원본)	1회 정제 후 (distilled_final)	변화량 / 상태 비고
전체 원자 수	553,500 (수소/물 없음)	586,914 (물 33,414 포함)	수화 셸 네트워크 반영
결합각 이상치 (> 4.0σ)	2,478개 이상 (최대 8.1σ)	0개	All within 4.0σ 달성
결합길이 이상치 (> 4.0σ)	다수 존재 (Max 5.0σ, Z=4.99)	0개	Max Z = 0.786으로 수렴
Bond Angle RMSD	1.434°	0.224° (Z = 0.122)	↓ 84.4% 감소
Bond Length RMSD	0.0171 Å	0.0007 Å (Z = 0.037)	↓ 95.9% 감소
Ramachandran Favored	100.00%	99.99%	70,920개 잔기 전역 선호 영역 유지
Ramachandran Outliers	0.00%	0.00%	이상치 0건 유지
Peptide Plane 비틀림	–	Twisted Proline/General 0.00%	펩타이드 주쇄 평면성 복원
방향족 고리 최대 편차	0.047 Å (TrpEr 224 등)	0.002 Å (TrpBB 84)	평면 오차 95% 이상 소거
All-atom Clashscore	1.35 (비수소 단독)	2.63 (물 3.3만 개 포함)	백분위 상위권 범위 유지
MolProbity Score	0.86	1.05	초고해상도 구조 수준 범위 안착
시작점 대비 좌표 변위	기준점 (0.000 Å)	Mean 0.635 Å / Max 2.766 Å	전체 나노케이지 외곽 Envelope 보존

Single-pass minimization benchmark without parameter overfitting

Geometry Restraints: GeoStd + CDL v1.2

3. 데이터 관점의 기술적 고찰

1) 공유결합 기하학(Covalent Geometry)의 정상화

원본 모델의 경우 Clashscore는 1.35로 양호하게 기록되었으나, 내부적으로는 Arg46-Pro47( $8.1\sigma$ , $128.33^\circ$ ), Gly48-Pro49( $7.8\sigma$ , $127.81^\circ$ ) 등 프롤린 인접 주쇄 및 곁사슬에서 2,478개 이상의 결합각 이상치( $>4\sigma$ )가 관측되었습니다.
1회차 최소화 연산 결과, 전역 결합각 RMSD가 1.434°에서 0.224°로 대폭 감소하였으며, $4.0\sigma$ 를 초과하는 결합각 및 결합길이 이상치가 단 한 건도 남지 않고 전량 해소되었습니다.

2) 입체 충돌(Clashscore) 변화에 대한 물리적 해석

정제 모델의 Clashscore는 1.35에서 2.63으로 소폭 조정되었습니다.
이는 비정상적으로 벌어져 있던 2,400여 개의 결합각이 표준 물리 화학적 각도(Ideal geometry)로 복원되고, 33,414개의 물 분자 셸 및 비결합 상호작용 계산을 위한 가상 라이딩 수소(Riding H) 환경이 포함되면서 원자 밀도가 크게 증가함에 따라 나타난 물리적 수렴 결과입니다.
58만 개 원자가 밀집된 상태에서도 Clashscore 2.63과 종합 MolProbity Score 1.05를 기록하여, wwPDB 기준 상위 백분위의 고품질 기하학적 완성도를 안정적으로 유지하고 있음을 확인했습니다.

3) 거대 복합체 Envelope의 안정성

정제 전후의 평균 원자 좌표 변위(Mean shift)는 0.635 Å로, Cryo-EM 전자밀도 맵이 정의하는 나노케이지 고유의 3차원 형태학적 골격을 무너뜨리지 않는 범위 내에서 국소적인 물리적 왜곡만을 외과적으로 교정했음을 보여줍니다.

4. 요약 및 향후 과제

본 실증은 $17.45\text{ \AA}$ 저해상도 맵에서 비롯된 초대형 조립체의 체계적 기하학 결함을 단 1회의 물리 기반 정제(Single-pass minimization)만으로도 유의미하게 안정화할 수 있음을 보여줍니다.

향후에는 잔존하는 Asn/Gln/His 측쇄의 180° 배향 플립 최적화 및 수소결합 네트워크의 추가 미세 조율을 진행하여 접촉 계면의 물리적 완성도를 더욱 고도화할 예정입니다.

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

2026년 10월 11일

[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일

[연구 노트] 데 노보(De novo) 바인더 설계와 재조합 단백질 발현 시스템 최적화 가이드
인공지능 기반 단백질 설계(RFdiffusion, ProteinMPNN 등)를 통해 인실리코(In silico) 상에서 결합 친화도와 접힘 신뢰도가 검증된 바인더(Binder)라도, 실제 생체 외(In vitro) 검증 단계에서 발현이 실패하거나 불용성 응집체(Inclusion bodies)를 형성하는 사례가 빈번하게 발생합니다.

단백질의 물리화학적 특성, 2차/3차 구조적 제약, 번역 후 변형(PTMs)의 필요성은 컴퓨터 디자인 단계부터 고려되어야 하며, 이에 부합하는 숙주 발현 시스템(Expression host)을 선제적으로 매핑해야 실험적 검증의 성공률을 극대화할 수 있습니다.

본 글에서는 단백질 발현 표준 프로토콜(STAR Protocols, 2024; doi: 10.1016/j.xpro.2024.103094)의 체계적 의사결정 프레임워크를 기반으로, 바인더 디자인 단계에서 점검해야 할 분자 수준의 변수와 주요 유전자 발현 시스템의 특성을 정리합니다.

1. 바인더 인실리코 설계 단계에서 고려해야 할 생산성 파라미터

컴퓨팅 파이프라인에서 시퀀스를 확정하기 전, 실제 발현 가능성(Developability)과 관련하여 다음 세 가지 요소를 사전에 분석해야 합니다.
1. 이황화 결합(Disulfide bond)의 개수 및 위치
  - 시스테인 잔기 간의 공유 결합은 단백질 열역학적 안정성을 높이지만, 환원 환경인 세포질(Cytoplasm)에서는 자연적으로 형성되기 어렵습니다.
  - 단일 도메인 바인더나 나노바디 골격에 다수의 이황화 결합이 필수적인 구조라면 원핵생물보다는 진핵생물 또는 산화 환경이 조성된 특수 균주 선정이 강제됩니다.
2. 소수성 패킹과 표면 전하 분포 (Aggregation Propensity)
  - 표적 단백질과의 결합 자유 에너지를 높이기 위해 계면에 과도한 소수성 패킹을 도입할 경우, 발현 과정에서 비특이적 응집(Non-specific aggregation)을 유발해 수율이 급감할 수 있습니다.
3. 글리코실화(Glycosylation) 부위 모티프 검증
  - 진핵 숙주 시스템을 사용할 계획이라면, 바인더 결합 에피토프 인근에 의도치 않은 $N$ -연결 당화 모티프( $\text{Asn-X-Ser/Thr}$ )가 위치하지 않도록 인실리코 서열 필터링이 선행되어야 합니다. 당화 사슬의 입체 장애가 표적 결합을 차단할 수 있기 때문입니다.
2. 주요 유전자 발현 시스템별 특성 및 바인더 적용 기준

1) 대장균 (Escherichia coli) 발현계
- 시스템 특성: 세균 기반 시스템으로 증식이 빠르고 배양 비용이 저렴하며, 유전자 조작이 용이하여 고수율(High-yield) 단백질 생산의 1차 표준으로 쓰입니다.
- 장점: 최소 장비(진탕 배양기 등)로 수일 내 그램(g) 단위의 단백질 확보가 가능합니다.
- 한계: 진핵생물 고유의 복잡한 번역 후 변형(당화 등)이 불가능하며, 다중 이황화 결합을 가진 단백질 발현 시 불용성 응집체(봉입체)로 석출되기 쉽습니다.
- 바인더 설계 적용 기준: 분자량 $100\,\text{kDa}$ 이하, 이황화 결합이 없거나 1~3개 이하로 제한된 미니프로테인(Miniprotein), 알파 헬릭스 번들 계열의 데 노보 바인더에 가장 적합합니다. 봉입체 형성 억제를 위해 분자 샤페론(GroEL/ES, DnaK 등) 공발현이나 가용화 태그(MBP, SUMO) 융합 전략이 유효합니다.
2) 효모 (Yeast; Pichia pastoris, Saccharomyces cerevisiae)
- 시스템 특성: 단세포 진핵생물로 원핵생물의 고밀도 배양 편의성과 진핵생물의 단백질 접힘 경로(소포체/골지체)를 모두 갖추고 있습니다.
- 장점: 내독소(LPS) 오염 위험이 없고 분비 생산이 가능하여 다운스트림 정제 공정이 단순화됩니다. 사이토카인, 성장 인자, 단일 사슬 항체 단편(scFv) 생산에 널리 활용됩니다.
- 한계: 포유류와 상이한 고만노오스형(Hypermannose) 당화 패턴을 부가하므로, 면역원성 유발 가능성이 있습니다.
- 바인더 설계 적용 기준: 이황화 결합 형성이 요구되나 포유류 특이적 복합 당화가 결합에 영향을 미치지 않는 스캐폴드 단백질 생산에 추천됩니다.
3) 곤충 세포 (Insect Cells / BEVS)
- 시스템 특성: 바큘로바이러스 발현 벡터 시스템(Baculovirus Expression Vector System)을 기반으로 진핵 단백질의 접힘을 정밀하게 모사합니다.
- 장점: 바이러스 유전체 용량이 커 다중 소단위체(Multi-subunit) 복합체(최대 17개 서브유닛 동시 발현 보고) 및 구조적으로 복잡한 막 단백질(IMPs) 생산에 탁월합니다.
- 한계: 바이러스 스톡 제작 및 감염 단계로 인해 생산 소요 기간(수 주)이 길고, 포유류의 말단 시알산(Sialic acid) 당화 구조는 결여되어 있습니다.
- 바인더 설계 적용 기준: 표적 복합체가 다량체이거나, 막 수용체를 표적으로 삼는 복잡한 구조체의 파트너 발현 및 검증에 최적화되어 있습니다.
4) 포유류 세포 (Mammalian Cells; HEK293, CHO)
- 시스템 특성: 인간 및 동물 유래 세포주로, 천연 생체 조건과 동일한 단백질 접힘, 샤페론 지원, 복합 당화( $N$ -/ $O$ -glycosylation) 환경을 제공합니다.
- 장점: 완전한 형태의 항체(IgG), Fc 융합 단백질, 복합 막 단백질의 구조적 활성을 100% 보존할 수 있습니다.
- 한계: 무균 설비와 $\text{CO}_2$ 인큐베이터 등 고가의 인프라가 필요하며 배양 비용이 높습니다.
- 바인더 설계 적용 기준: 전장 면역글로불린(Full-length IgG) 형태의 항체 디자인이나, 당화 패턴 자체가 결합 에피토프 인식에 직접 관여하는 정밀 치료용 바인더 생산에 필수적입니다.
5) 무세포 단백질 발현계 (Cell-Free Protein Synthesis, CFPS)
- 시스템 특성: 세포막을 파쇄한 라이세이트(E. coli, Wheat germ 등)에 기질과 DNA/RNA를 직접 투입하여 시험관 내에서 번역을 유도하는 시스템입니다.
- 장점: 세포 생존율에 치명적인 세포 독성 단백질 생산이 가능하며, 비천연 아미노산(Non-canonical amino acids) 도입 및 미세 유체 기반 초고속 스크리닝이 가능합니다.
- 한계: 생산 단가가 높고 배치당 회수 수율(μ $\mu\text{g}\sim\text{mg}/\text{mL}$ )이 제한적입니다.
3. 성공적인 바인더 생산을 위한 단계별 결정 트리
[바인더 시퀀스 분석] │ ├─ 1. PTM(당화 등) 불필요 / 단순 도메인 / 이황화 결합 ≤ 3개 │ └─► 대장균 (E. coli) 시스템 선택 (BL21, SHuffle 균주 등) │ ├─ 2. 분비 생산 필요 / 이황화 결합 다수 / 단순 진핵 당화 허용 │ └─► 효모 (Pichia pastoris) 시스템 선택 │ ├─ 3. 대형 다중 복합체 / 막 단백질 표적 │ └─► 곤충 세포 (Baculovirus / BEVS) 시스템 선택 │ └─ 4. 전장 항체(IgG) 구조 / 인간형 복합 당화 필수 └─► 포유류 세포 (HEK293 일시적 발현 또는 CHO 안정 세포주)
컴퓨터 단백질 디자인은 시퀀스 생성 단계에서 끝나는 것이 아니라, 해당 단백질이 물리적으로 접히고 분비될 수 있는 세포 환경을 정의하는 단계까지 아우를 때 실제 기능성 바이오 물질로 완성될 수 있습니다.

참고 자료 (References)
1. STAR Protocols (2024). A decision tree protocol for selecting recombinant protein expression systems. STAR Protocols, 5, 103094. https://doi.org/10.1016/j.xpro.2024.103094
2. Watson, J. L., et al. (2023). De novo design of protein structure and function with RFdiffusion. Nature, 620(7976), 1089–1100. https://doi.org/10.1038/s41586-023-06415-8
2026년 10월 07일
[연구 노트] 단백질 구조 벤치마크 참여 현황 및 향후 분석 계획
HanjariFold 개발팀은 최근 진행된 단백질 3차원 구조 예측 및 정제 분야의 국제 벤치마크 프로그램(CASP)에 참여하여, 공개된 대상 타깃들에 대한 계산 구조 제출 작업을 모두 완료하였습니다.

현재는 공식 제출 기간이 마감되어 주최 측 및 평가위원회의 심층 평가가 진행되는 단계에 접어들었습니다.

연구 개요 및 적용 파이프라인

이번 작업에서는 최신 딥러닝 기반 모델링 기법을 통해 확보된 초기 구조 좌표를 바탕으로, HanjariFold의 독자적인 정제 및 물리 최적화 알고리즘을 연계 적용하였습니다.
- 물리적 힘장 최소화 (HanjariNebula Engine):
  
  원자 간 비물리적 입체 충돌(Steric Clash)을 해소하고 수소결합 및 표면 수화 껍질을 연속적으로 배치하여 위치 에너지를 안정화하였습니다.
- 기하 구속 기반 정밀 조율 (CosmicDistiller Engine):
  
  결정학적 표준에 부합하는 주쇄 2면각(Ramachandran distribution) 및 화학 결합 기하 구조를 정돈하여 구조의 물리화학적 신뢰도를 보정하였습니다.
데이터 공개 원칙 및 엠바고 준수

벤치마크 평가의 독립성과 공정성을 존중하기 위해, 대회 운영 규정과 비공개 실험 데이터 보호 원칙에 따라 다음과 같은 기준을 유지합니다.
1. 타깃별 세부 정보 및 좌표 비공개:
  
  공식 평가가 완료되기 전까지 개별 타깃의 식별자(Target ID) 및 제출된 세부 원자 좌표는 외부에 공개하지 않습니다.
2. 사전 비교 분석 배제:
  
  실험 그룹의 정답 결정 구조(Ground Truth)가 공식적으로 릴리즈되기 전까지 어떠한 형태의 예단이나 사전 정합성 분석도 진행하지 않습니다.
향후 분석 계획 (12월 이후)

주최 측의 공식 평가 결과가 공표되고 실험 구조 데이터가 정식으로 공개되는 12월 이후, 다음과 같은 분석 포스팅을 순차적으로 게재할 예정입니다.
- 실험 정답 구조(Ground Truth)와의 1:1 정합성 분석:
  
  실제 고해상도 X-선/Cryo-EM 구조 대비 HanjariFold 정제 엔진이 달성한 백본 정밀도 및 곁사슬 배치 개선 정도를 객관적 수치로 검증합니다.
- MolProbity 및 물리적 지표 변화 추적:
  
  초기 예측 좌표와 최종 제출 좌표 사이의 원자 충돌 점수(Clashscore), 라마찬드란 적합도, 결합 에너지 변화를 종합 비교합니다.
- 수화 정책(Solvation Policy)에 따른 영향 고찰:
  
  용매 분자 배치 유무가 결합 포켓 및 복합체 안정성에 미친 실제 효과를 데이터 기반으로 공유합니다.
HanjariFold는 앞으로도 순수 학술 연구 목적에 부합하는 투명하고 재현 가능한 단백질 구조 정제 파이프라인을 구축해 나가겠습니다. 12월 공식 결과 발표 이후 심도 있는 데이터 분석 글로 다시 찾아뵙겠습니다.
2026년 10월 07일

[카테고리:] 연구노트

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. 요약: 디지털 크리스탈로그래피 기술이 연 새로운 지평

[기술 리포트] 58만 원자급 대칭성 나노케이지(PDB 9OM3) 1회차 지오메트리 정제 결과 분석

개요

1. 대상 시스템 및 계산 환경

2. 정제 전후 검증 지표 정량 비교

3. 데이터 관점의 기술적 고찰

1) 공유결합 기하학(Covalent Geometry)의 정상화

2) 입체 충돌(Clashscore) 변화에 대한 물리적 해석

3) 거대 복합체 Envelope의 안정성

4. 요약 및 향후 과제

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

[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

[연구 노트] 데 노보(De novo) 바인더 설계와 재조합 단백질 발현 시스템 최적화 가이드

1. 바인더 인실리코 설계 단계에서 고려해야 할 생산성 파라미터

2. 주요 유전자 발현 시스템별 특성 및 바인더 적용 기준

1) 대장균 (Escherichia coli) 발현계

2) 효모 (Yeast; Pichia pastoris, Saccharomyces cerevisiae)

3) 곤충 세포 (Insect Cells / BEVS)

4) 포유류 세포 (Mammalian Cells; HEK293, CHO)

5) 무세포 단백질 발현계 (Cell-Free Protein Synthesis, CFPS)

3. 성공적인 바인더 생산을 위한 단계별 결정 트리

참고 자료 (References)

[연구 노트] 단백질 구조 벤치마크 참여 현황 및 향후 분석 계획

연구 개요 및 적용 파이프라인

데이터 공개 원칙 및 엠바고 준수

향후 분석 계획 (12월 이후)