[Bioinformatics] PyMOL을 활용한 소분자 화합물(Ligand) 구조 중첩 및 분석 가이드

단백질-리간드 복합체 연구에서 여러 유도체의 결합 모드를 비교할 때, 구조 중첩(Superposition)은 필수적인 과정입니다. 흔히 사용하는 ssm 명령어는 단백질 2차 구조 정렬에 최적화되어 있으므로, 소분자 화합물을 정교하게 겹치기 위해서는 리간드 특성에 맞는 별도의 방법이 필요합니다.

1. PyMOL 리간드 정렬 방법 (Structural Alignment)

PDB 형식으로 생성된 화합물 모델들을 PyMOL에 불러온 뒤, 공통된 화학적 골격(Core Scaffold)을 기준으로 정렬하면 각 유도체의 치환기가 포켓 외부로 뻗어 나가는 방향과 상호작용의 차이를 한눈에 비교할 수 있습니다.

방법 A: align 또는 super 명령어 (가장 간편한 자동 정렬)

PyMOL이 두 화합물 간의 유사한 원자들을 자동으로 식별하여 최적의 위치로 겹쳐주는 방식입니다.

• 사용 사례: 전체적인 결합 유사성을 빠르게 확인하고 싶을 때 적합합니다.

• 명령어:

  Python

  # 'target_ligand'를 'reference_ligand'에 맞춰 자동 정렬
  align target_ligand, reference_ligand
  

방법 B: pair_fit 명령어 (공통 골격 중심의 정밀 정렬)

사용자가 지정한 **특정 원자들(예: 핵심 페놀 고리의 탄소들)**만을 수학적으로 정확하게 일치시키는 방식입니다.

• 사용 사례: 공통 골격은 고정하고 치환기의 변화에 따른 구조적 변위만을 정밀하게 분석할 때 가장 추천되는 방법입니다.

• 실행 순서:

  1. 중첩의 기준이 될 두 화합물의 공통 원자들을 각각 선택(Selection)합니다.

  2. 선택한 원자 번호를 매칭하여 명령어를 실행합니다.

  Python

  # 선택한 지점들을 기준으로 정밀하게 중첩
  pair_fit target_selection, reference_selection
  

2. 구조 중첩 작업이 중요한 이유

리간드들을 하나로 겹쳐보는 작업은 단순한 시각화를 넘어, 신약 설계의 핵심적인 판단 근거를 제공합니다.

• SAR(Structure-Activity Relationship) 분석: 핵심 골격을 고정한 채 유도체들을 중첩해 보면, 특정 작용기 도입 시 화합물이 포켓의 물리적 공간을 얼마나 점유하고 루프 구조를 어떻게 밀어내는지 시각적으로 즉각 확인할 수 있습니다.

• Clash(원자 충돌) 확인: 정밀화된 단백질 구조 모델과 중첩했을 때, 유도체의 새로운 작용기가 단백질 원자와 입체적 충돌을 일으키지 않는지 판단하는 명확한 기준이 됩니다.

💡 실무를 위한 팁

화합물 데이터를 다룰 때 발생할 수 있는 시행착오를 줄이기 위한 몇 가지 팁을 공유합니다.

1. 원자 이름 및 순서 통일: PDB 변환 시 원자 이름(Atom Name)이 제각각이면 PyMOL이 대응 원자를 찾는 데 혼란을 겪을 수 있습니다. RDKit 등의 도구를 사용해 PDB를 생성할 때 원자 순서를 가급적 통일하는 것이 좋습니다.

2. PyMOL Wizard 활용: 다수의 오브젝트를 관리하기 번거롭다면 Wizard -> Copy-to-object 기능을 활용하여 하나의 오브젝트 안에서 여러 컨포머(Conformer)를 비교하는 것도 효율적입니다.

3. 3D 좌표 생성의 중요성: SMILES로부터 PDB를 생성할 때는 반드시 수소 원자를 추가하고 에너지 최적화 과정을 거쳐 신뢰할 수 있는 3차원 좌표를 확보해야 합니다.

이 포스팅이 리간드 기반 신약 설계 및 구조 최적화 연구를 진행하시는 분들께 도움이 되기를 바랍니다.