[태그:] AMBER14SB

  • [구조 분석] 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 -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. 핵심 기술 고찰

    1. 상위 0.1% 수준의 MolProbity Score 달성 (2.03 → 1.24 → 0.75)

      • PDB 상위 1 ~ 2Å 고해상도 결정 구조에서도 보기 드문 0.75 점수를 달성했습니다. 1단계 Nebula가 거친 원자 충돌을 걷어내고(5.69 → 0.81), 2단계 CosmicDistiller가 CCTBX 기하 구속으로 결합 각도와 결합 길이를 미세 조정하여 이상적인 지표를 완성했습니다.

    2. 골격 2면각의 이상적 복원 (Rama Favored 98.65% → 100.00%)

      • 분자 역학 힘장 최소화 과정에서 미세하게 완화되었던 주쇄(Backbone) 각도를 CosmicDistiller의 기하학적 재배치가 완벽하게 정렬하여 Favored 100.00%, Outliers 0.00%를 기록했습니다.

    3. 수화 쉘(Solvation Shell)의 완성 (0 → 51 → 57 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Å → 0.531Å)

    • HanjariNebula (1단계): 원본 대비 RMSD가 단 0.207Å에 불과합니다. 이는 단백질 전체 골격의 위치를 거의 고정한 상태에서 입체 장애(Clashscore)를 유발하던 국소 원자 충돌만을 외과 수술식으로 완화했음을 의미합니다.

    • CosmicDistiller (2단계): 2차 CCTBX 기하 구속 및 2면각 보정을 거치며 RMSD는 0.531Å로 수렴했습니다. 일반적으로 단백질 결정학에서 1.8Å 분해능 구조의 자연스러운 열진동(Thermal fluctuation) 및 분해능 한계 변위가 약 0.5 ~ 0.8Å 수준임을 감안할 때, 원본의 접힘 상태를 엄격히 보존하면서 라마찬드란 선호 영역 비율을 0%로 끌어내리는 데 필요한 최적의 미세 변위만 발생했음을 증명합니다.

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

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

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

      X-선 회절 실험을 통해 결정 격자 내에서 고정된 58개의 물 분자(산소 원자 58개, 수소 미포함)가 등록되어 있습니다.

    • 1단계 HanjariNebula (51 HOH / 올아톰 수화):

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

    • 2단계 CosmicDistiller (57 HOH / 결정학적 수렴):

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

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

  • [기술 노트] 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Å 이상 비정상적으로 겹치는 입체 장애(Steric Overlap)의 밀도를 1,000개 원자당 개수로 나타낸 수치입니다 (낮을수록 우수).

      • 판독 기준: 2.00 미만으로 떨어질 때 원자 간 입체 충돌이 실질적으로 해소된 것으로 판단하며, 도킹(Docking) 및 분자동역학(MD) 시뮬레이션의 안정성을 확보하는 데 중요한 지표입니다.

    • Ramachandran Favored / Outliers

      • 정의: 단백질 백본의 2면각(φ, ψ)이 에너지적으로 선호되는 정상 구역(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 좌표 및 검증 리포트를 일괄 내려받을 수 있습니다.