Graph-Theoretic Computation of Degree-Based Topological Indices for Hepatocellular Carcinoma Therapeutics: A Mathematical QSPR Framework for Physicochemical Property Prediction

기술 요약: 간세포암 치료제의 차수 기반 위상 지수를 이용한 그래프 이론적 계산

문제 정의
약물 후보 물질의 물리화학적 성질(끓는점, 인화점, 분극도, 표면 장력, 몰 부피 등)을 실험적으로 결정하는 것은 자원 집약적이며 대규모 화합물 라이브러리에 적용하기 어렵다. 정량적 구조-물성 관계(QSPR) 프레임워크는 분자를 그래프로 모델링함으로써 수학적 대안을 제시하지만, 간세포암(HCC) 치료제에 대해 포괄적인 차수 기반 위상 지수 세트를 적용하는 것과 관련하여 문헌상에 특정 공백이 존재한다. 본 연구는 다섯 가지 임상 승인 HCC 약물의 분자 구조로부터 직접적으로 물리화학적 거동을 예측하기 위한 계산 효율적인 그래프 이론적 프레임워크의 필요성을 다룬다.

연구 방법론
본 연구는 다섯 가지 HCC 약물인 푸티비티닙(Futibatinib), 이필리무맙(Ipilimumab), 펨브롤리주맙(Pembrolizumab), 페미가티닙(Pemigatinib), 레고라페닙(Regorafenib)을 모델링하기 위해 화학 그래프 이론(CGT)을 채택하였다.

• 그래프 구축: 각 분자는 비수소 중원자(non-hydrogen heavy atoms)를 정점($V$)으로 하고 공유 결합을 에지($E$)로 하는 단순 연결 분자 그래프 $G(V, E)$로 표현된다. 수소 원자는 표준 관례에 따라 생략되었다.
• 위상 지수: 여섯 가지 차수 기반 위상 지수가 각 분자의 분자 그래프에 대해 계산되었다. 이 지수들은 차수의 거듭제곱 합을 통해 정점 차수($d_v$)에 대한 민감도가 달라진다:
  • 제1 자그레브 지수 ($M_1$): $(d_u + d_v)$의 합.
  • 제2 자그레브 지수 ($M_2$): $(d_u \times d_v)$의 합.
  • 포가튼 지수 ($F$): $(d_u^2 + d_v^2)$의 합.
  • 예멘 지수 ($Y$): $(d_u^3 + d_v^3)$의 합.
  • 차수 지수 ($D$): $(d_u^6 + d_v^6)$의 합.
  • 합-연결성 지수 ($S$): $(d_u^4 + d_v^4)$의 합.
• 데이터 및 분석: 실험적 물리화학적 데이터(끓는점, 인화점, 분극도, 표면 장력, 몰 부피)는 ChemSpider 데이터베이스에서 추출하였다. 모든 30가지 지수-물성 조합에 대해 SPSS(v25.0)를 사용하여 선형 회귀 모델($T = \alpha + \beta t_i$)을 구축하였다. 모델 성능은 피어슨 상관 계수($r$), 결정 계수($R^2$), 통계적 유의성($p < 0.05$)을 통해 평가되었다.

주요 결과
본 연구는 위상 지수가 서로 다른 물성에 대해 갖는 예측 능력에 관한 뚜렷한 발견을 제시하였다:

• 표면 장력(ST)의 우세성: 표면 장력은 위상 지수와 통계적으로 유의미한 상관관계($p < 0.05$)를 보인 유일한 물리화학적 성질이었다. 여섯 가지 지수 중 다섯 가지가 표면 장력에 대한 유의미한 QSPR 모델을 생성하였으며, $R^2$ 값은 0.907에서 0.978 사이였다.
• 지수 성능 계층: 예멘 지수($Y(G)$)가 가장 높은 예측 충실도($R^2 = 0.978, p = 0.011$)를 달성하였으며, 합-연결성 지수($S(G), R^2 = 0.965$)와 포가튼 지수($F(G), R^2 = 0.967$)가 그 뒤를 바짝 쫓았다. 제1 및 제2 자그레브 지수 또한 유의성을 보였으나 $R^2$ 값은 더 낮았다.
• 비유의적 물성: 끓는점, 인화점, 분극도, 몰 부피에 대한 모델은 통계적 유의성을 확보하는 데 실패하였다. 저자들은 이를 적은 표본 크기(페미가티닙의 데이터 누락으로 인해 끓는점과 인화점의 경우 $N=4$로 감소)와, 차수 기반 기술자만으로는 완전히 포착할 수 없는 전자적 및 입체적 요인에 의존하는 물성의 다요인적 특성 때문이라고 설명한다.
• "차수 지수"의 이상 현상: 6제곱 합을 사용하는 차수 지수($D(G)$)는 표면 장력과 높은 상관관계($r=0.933$)를 보였으나 통계적 유의성($p=0.067$)에는 도달하지 못했다. 저자들은 이를 $D(G)$ 값의 극단적인 동적 범위(거의 두 자릿수 차이의 범위를 가짐)로 인한 소표본 선형 회귀에서의 분산 팽창 현상으로 설명한다.

의의 및 주장
본 논문은 오직 분자 그래프 구조만을 기반으로 HCC 치료제의 표면 활성 물리화학적 거동을 예측할 수 있는 수학적으로 다루기 쉬운 프레임워크를 확립하였다.

• 기전적 통찰: 저자들은 차수 기반 지수(특히 $Y(G)$ 및 $F(G)$와 같은 3제곱 및 2제곱 형태)가 표면 장력을 예측하는 데 성공한 이유가 이러한 지수들이 고차수 정점(예: 방향족 고리 접합부)의 밀도와 일치하기 때문이라고 제안한다. 이러한 구조적 특징은 액체상 응집력을 지배하는 런던 분산 상호작용과 상관관계가 있다.
• 실질적 적용: 실험적 측정 없이 표면 장력을 계산적으로 예측할 수 있다는 본 연구의 결과는, 특히 HCC 약물 전달 시스템을 위한 습윤성, 퍼짐성 및 유화 최적화 측면에서 제약 제형 설계에 직접적인 시사점을 준다.
• 범위의 한 한계: 저자들은 차수 전용 기술자가 3D, 전자적 또는 거리 기반 기술자 없이는 다른 물성(BP, FP, P, MV)을 설명하기에 불충분하다는 점을 언급하며, 결과의 일반화 가능성에 대해 신중한 태도를 유지하였다. 이들은 향후 연구에서 화합물 집단을 확장하고 이러한 한계를 해결하기 위해 비선형 머신러닝 접근법을 통합할 것을 제안한다.