In Silico Screening and Molecular Dynamics Validation of Bryophyllum Pinnatum Flavonoids as Pi3kγ Inhibitors For Rheumatoid Arthritis

기술 요약: 류마티스 관절염 치료를 위한 Bryophyllum Pinnatum 플라보노이드의 PI3Kγ 억제제로서의 인실리코(In Silico) 스크리닝 및 분자 동역학 검증

문제 제기
류마티스 관절염(RA)은 면역 세포 및 활막 섬유아세포에서 인지질 이노시톨 3-키나아제 감마(PI3Kγ) 이형체의 비정상적인 활성화에 의해 유발되는 진행성 관절 파괴를 특징으로 하는 만성 자가면역 질환이다. PI3Kγ는 염증 세포(호중구 및 대식세포)의 침윤을 촉진하고 PI3K/AKT/mTOR 경로를 악화시켜, 염증성 사이토카인(TNF-α, IL-1β, IL-6)의 증가와 활막 섬유아세포의 과성장을 초래한다. 기존의 치료법들이 존재하지만, 새로운 식물 기반 치료제에 대한 필요성이 존재한다. 아프리카 민간요법에서 염증 질환에 널리 사용되는 약용 식물인 Bryophyllum pinnatum은 생리활성 플라보노이드를 함유하고 있으나, PI3Kγ에 대한 구체적인 분자적 기전은 계산 과학적 검증을 통해 아직 밝혀지지 않았다.

방법론
본 연구는 B. pinnatum 유래 플라보노이드를 잠재적인 PI3Kγ 억제제(PDB ID: 2A4Z)로서 평가하기 위해 종합적인 in silico 워크플로우를 채택하였다. 연구 방법론은 다음 단계로 진행되었다:

1. 단백질 및 리간드 준비: PI3Kγ 구조는 OPLS4 포스 필드를 사용하여 전처리(결합 차수, 수소 추가, 루프 채우기)되었다. B. pinnatum 화합물 라이브러리는 PubChem에서 검색되었으며, 3D 형태로 변환된 후 LigPrep을 사용하여 pH 7.0 ± 0.2에서 이온화되었다.
2. 파마코포어 모델링 및 스크리닝: 공결정화된 리간드 복합체를 기반으로 E-파마코포어 모델이 생성되었다. 7개의 소수성 그룹과 10개의 방향족 고리를 특징으로 하는 이 모델은 최소 두 개의 핵심 특징을 가진 화합물만을 선별하기 위해 사용되었다.
3. 분자 도킹: 선별된 화합물(총 25개)은 Glide Extra Precision (XP) 알고리즘을 사용하여 활성 부위에 도킹되었다. 결합 부위는 공결정화된 리간드를 중심으로 하는 그리드 박스로 정의되었다.
4. ADMET 프로파일링: 상위 히트 화합물들은 SwissADME 서버를 사용하여 흡수, 분포, 대사, 배설 및 독성(ADME)을 평가하였으며, 특히 리핀스키의 법칙(Lipinski's Rule of Five), 극성 표면적(TPSA), 피부 투과성(logKp) 준수 여부를 확인하였다.
5. 결합 자유 에너지 계산: 도킹 점수를 검증하기 위해 Prime MM-GBSA를 사용하여 상위 복합체들의 결합 자유 에너지(ΔGbind)를 계산하였다.
6. QSAR 모델링: 리드 화합물들의 pIC₅₀ 값을 예측하기 위해 ChEMBL 데이터베이스의 화합물을 이용한 AutoQSAR 모델을 개발하였다.
7. 분자 동역학(MD) 시뮬레이션: Desmond 소프트웨어(Schrödinger suite)를 사용하여 100-나노초(ns) MD 시뮬레이션을 수행하였으며, OPLS4 포스 필드와 TIP3P 용매 모델을 적용하였다. 단백질-리간드 복합체의 안정성은 RMSD(제곱평균제곱근 편차)와 RMSF(제곱평균제곱근 변동)를 통해 평가되었으며, 상호작용(수소 결합, 소수성, 이온 및 물 가교)을 모니터링하였다.

주요 결과

• 도킹 점수 및 상호작용: 세 가지 플라보노이드인 루테올린(Luteolin), 퀘르세틴(Quercetin), 캠페롤(Kaempferol)은 공결정화된 리간드(-9.628 kcal/mol) 및 표준 약물인 토파시티닙(Tofacitinib, -6.205 kcal/mol)보다 우수한 억제 잠재력을 보였다.
  • 루테올린: -11.044 kcal/mol (LYS 833, VAL 882, GLU 880와의 수소 결합).
  • 퀘르세틴: -10.723 kcal/mol (SER 806, VAL 882, LYS 833와의 수소 결합).
  • 캠페롤: -10.553 kcal/mol (SER 806, LYS 833, VAL 882와의 수소 결합; TYR 867과의 PI-PI 적층).
• ADMET 특성: 모든 리드 화합물은 리핀스키의 법칙(위반 사항 0건)을 충족하였고, 유리한 분자량을 나타냈으며, TPSA 값이 140 Ų 미만으로 나타나 양호한 경구 생체 이용률을 보였다.
• 결합 자유 에너지 (MM-GBSA): 리드 화합물들은 토파시티닙보다 현저히 강한 결합 친화도를 보였다.
  • 루테올린: -45.07 kcal/mol
  • 캠페롤: -44.79 kcal/mol
  • 퀘르세틴: -44.03 kcal/mol
  • (비교: 공결정화 리간드 -44.57 kcal/mol; 토파시티닙 -21.04 kcal/mol).
• QSAR 예측: AutoQSAR 모델(R² = 0.6646, Q² = 0.6625)은 높은 억제 효능을 예측하였으며, 루테올린이 가장 높은 pIC₅₀ 값(7.053)을, 퀘르세틴(6.952), 캠페롤(6.627)이 그 뒤를 이었다.
• MD 시뮬레이션 안정성: 100 ns 동안 단백질 백본(backbone) RMSD는 일반적으로 3.6 Å 미만을 유지하여 구조적 무결성을 나타냈다.
  • 루테올린 복합체는 안정화되기 전 초기 형태 변화를 보였다.
  • 퀘르세틴 복합체는 20 ns 후에 안정화되었다.
  • 캠페롤 복합체는 80–100 ns 사이에 안정화되었다.
  • 토파시티닙은 안정화까지 약 40 ns가 소요되었다.
  • 궤적 전반에 걸쳐 복합체들은 LYS 833, VAL 882, GLU 880과 같은 잔기들과의 수소 결합, 소수성 접촉, 물 가교를 포함한 핵심 상호작용을 유지하였다.

의의 및 주장
저자들은 본 연구가 염증성 질환 치료를 위한 Bryophyllum pinnatum의 전통적 사용에 대한 강력한 계산 과학적 근거를 제공한다고 주장한다. 연구 결과는 루테올린, 퀘르세틴, 캠페롤이 결합 친화도 및 상호작용 안정성 측면에서 표준 약물인 토파시티닙보다 뛰어날 가능성이 있는, 강력하고 안정적이며 약물로서 적합한 PI3Kγ 억제제임을 시사한다.

본 논문은 이 화합물들이 RA에서 호중구/대식세포 침윤, 사이토카인 과잉 생산 및 활막 손상을 일으키는 PI3Kγ 매개 신호 전달 경로를 효과적으로 방해할 수 있다고 상정한다. 연구는 이러한 in silico 결과가 B. pinnatum 구성 성분의 항-RA 특성에 대한 강력한 과학적 근거를 제공한다고 결론짓는다. 다만, 저자들은 이러한 결과가 예비 단계이며 임상 적용 전 in vitro 키나아제 분석, 활막 염증 세포 모델, 그리고 콜라겐 유도 관절염 동물 모델을 통한 검증이 반드시 필요함을 명시하였다.