In silico design of potent antioxidant compounds using 2D-QSAR Molecular docking ADMET prediction and Molecular dynamic simulation

기술 요약: 강력한 항산화 화합물의 인 실리코(In Silico) 설계

문제 제기
활성 산소종(ROS)과 내인성 항산화 방어 체계 사이의 불균형에 의해 발생하는 산화적 스트레스는 노화 및 암, 심혈관 질환, 당뇨병, 알츠하이머병과 같은 퇴행성 질환의 주요 요인이다. 천연 및 합성 항산화제가 존재하지만, 많은 경우 낮은 용해도, 낮은 함량, 독성 등의 한계를 겪고 있다. 특히, 1,4-벤조옥사진(1,4-benzoxazine) 유도체는 항산화제로서 잠재력을 보여주었으나, 구조적 특성과 생물학적 활성, 결합 모드 및 약동학적 프로파일을 연결하는 포괄적인 계산 평가가 부족한 실정이다. 본 연구는 다단계 인 실리코 워크플로우를 사용하여 보다 강력하고 안정적이며 안전한 벤조옥사진 기반 항산화제를 합리적으로 설계하는 것을 목적으로 한다.

방법론
본 연구는 다음의 순차적인 단계를 포함하는 리간드 기반 및 구조 기반 약물 설계 접근법을 채택하였다:

1. 데이터셋 준비: 보고된 DPPH 라디칼 소거 활성(IC50)을 가진 28개의 알려진 1,4-벤조옥사진 유도체 데이터셋을 수집하였다. 구조는 3D로 변환되었으며, 분자 역학(MMFF)을 통해 에너지 최소화를 수행한 후, B3LYP/6-31G* 수준의 밀도 범함수 이론(DFT)을 통해 최적화되었다.
2. 2D-QSAR 모델링: PADEL-Descriptor를 사용하여 분자 기술자(descriptor)를 생성하였다. 데이터셋은 Kennard-Stone 알고리즘을 사용하여 훈련 세트(75%)와 테스트 세트(25%)로 나누었다. Material Studio v8.0에서 유전자 함수 알고리즘(GFA)과 다중 선형 회귀(MLR)를 사용하여 4개의 모델을 구축하였다. 모델 검증에는 내부($R^2_{trn}$, $Q^2_{LOO}$), 외부($R^2_{ext}$), Y-무작위화($cR^2_p$) 테스트가 포함되었다. 적용 범위는 Williams plot을 통해 평가되었다.
3. 분자 도킹: 최적의 QSAR 모델을 사용하여 새로운 유사체를 설계하는 가이드로 사용하였다. 28개의 기존 화합물과 5개의 신규 설계 유도체 모두를 미토콘드리아 전자 전달 및 세포 사멸 조절의 핵심 효소인 사이토크롬 c 퍼옥시다제(Cytochrome c Peroxidase, PDB ID: 2XO8)에 대해 Molegro Virtual Docker(MVD)를 사용하여 도킹하였다. 아스코르브산이 참조 표준으로 사용되었다.
4. 리간드 기반 설계: 가장 활성이 높은 템플릿인 화합물 20(pIC50 = 7.817)에 QSAR 모델의 기술자 기여도에 기초하여 특정 위치(D, X, Y)에 작용기(히드록실, 메톡시, 에톡시, 에틸, 클로레이트, 알데하이드)를 추가하여 변형하였다.
5. ADMET 및 약동학: SwissADME 및 pkCSM 서버를 사용하여 약물 유사성(Lipinski의 Rule of 5, Veber의 규칙)과 ADMET 특성(흡수, 분포, 대사, 배설, 독성)을 예측하였다.
6. 분자 역학(MD) 시뮬레이션: 상위 순위의 설계된 복합체(D1 및 D5)와 수용체 간의 안정성을 Schrödinger의 Desmond 소프트웨어를 사용하여 100 ns 시뮬레이션을 통해 검증하였다. RMSD, RMSF, 회전 반경(Rg), 수소 결합 등의 파라미터를 분석하였다.

주요 결과

• QSAR 모델 성능: 통계적 견고성($R^2_{trn} = 0.7889$, $R^2_{ext} = 0.7449$, $Q^2_{LOO} = 0.6889$, $cR^2_p = 0.6015$)을 기준으로 모델 4가 최적의 모델로 선정되었다. 이 모델은 GATS2m, GATS3m, SpMin6_Bhv 및 IC2 기술자를 활용하였다. Y-무작위화 테스트는 모델이 우연한 상관관계가 아님을 확인하였다.
• 설계된 화합물: 5개의 새로운 벤조옥사진 유도체(D1–D5)가 설계되었다. 이들은 템플릿(화합물 20, pIC50 = 7.817)과 아스코르브산(pIC50 = 7.817)을 능가하는 pIC50 값 범위(8.002 ~ 8.532)의 우수한 예측 항산화 활성을 나타냈다.
• 도킹 친화도: 설계된 화합물들은 강한 결합 친화도를 보였다. 화합물 D5가 가장 높은 MolDock 점수(-135.37 kcal/mol)를 기록하였고, D1(-134.38 kcal/mol), D4(-131.14 kcal/mol)가 그 뒤를 이었다. 이 점수들은 참조 물질인 아스코르브산(-82.772 kcal/mol)보다 유의미하게 높았다. 상호작용 분석 결과, ARG 48, TRP 51, HIS 181과 같은 잔기들과의 수소 결합, pi-cation, pi-pi 적층 및 소수성 상호작용을 포함한 다양한 결합 모드가 확인되었다.
• ADMET 프로파일: 모든 설계된 화합물은 Lipinski와 Veber의 규칙을 준수하여 양호한 경구 생체 이용률을 나타냈다. 이들은 높은 위장관 흡수율(90–94%), 낮은 혈액-뇌 장벽 투과성(log BB < 0.3), 유리한 대사 프로파일(CYP3A4 기질, CYP2D6 비억제제)을 보였다. 독성 예측 결과, D1, D2, D5는 비돌연변이성인 반면, D3와 D4는 돌연변이 유발 가능성이 있는 것으로 예측되었다.
• MD 시뮬레이션: 100 ns MD 시뮬레이션은 D1 및 D5 복합체의 구조적 안정성을 확인하였다. 단백질 백본 RMSD는 2.2–2.6 Å 사이에서 안정화되었으며, 리간드는 활성 부위 내에 안정적으로 결합되어 있었고, RMSD 값은 약 0.9 Å에서 안정화되었다. 수소 결합 분석 및 접촉 맵은 시뮬레이션 전반에 걸쳐 지속적이고 강력한 상호작용을 나타냈다.

의의 및 주장
저자들은 본 연구가 신규 항산화제를 위한 합리적 설계를 위해 통합된 인 실리코 파이프라인(QSAR, 도킹, MD, ADMET)의 유용성을 성공적으로 입증했다고 주장한다. 본 연구는 1,4-벤조옥사진 골격에 항산화 효능 및 사이토크롬 c 퍼옥시다제에 대한 결합 친화도를 향상시키는 특정 구조적 변형을 식별하였다.

본 논문은 설계된 벤조옥사진 유도체, 특히 D1, D2, D5가 산화적 스트레스 관련 질환에 대한 치료제로서 개발될 수 있는 유망한 후보 물질임을 결론짓는다. 그러나 저자들은 이러한 결과가 "추가적인 인 비트로(in vitro) 및 인 비보(in vivo) 실험의 대상이 되는 유망한 항산화제 역할을 할 수 있음"을 명시하며, 계산적 예측이라는 점을 밝혀 신중한 입장을 유지하고 있다. 본 연구는 미래의 합성 및 생물학적 테스트를 위한 검증된 구조적 틀을 제공함으로써 실험 연구와 실제적 치료 적용 사이의 간극을 메우는 것을 목표로 한다.