Hoạt tính gây độc tế bào HeLa của dẫn xuất flavonoid - Luận án tiến sĩ

Tổng quan về luận án

Nghiên cứu của nghiên cứu sinh Bùi Thị Phương Thúy, được thực hiện dưới sự hướng dẫn khoa học của PGS. Phạm Văn Tất (Trường Đại học Hoa Sen) và PGS. Trần Dương (Trường Đại học Sư phạm – Đại học Huế) tại Trường Đại học Khoa học – Đại học Huế, tập trung vào đề tài: "Nghiên cứu mối quan hệ cấu trúc - hoạt tính (QSAR) của các hợp chất flavonoid nhằm định hướng thiết kế và sàng lọc các chất có hoạt tính kháng ung thư cổ tử cung". Trong bối cảnh ung thư cổ tử cung do virus Human Papillomavirus (HPV) gây ra tiếp tục là gánh nặng y tế toàn cầu và các phác đồ hóa - xạ trị quy ước tồn tại nhiều độc tính tế bào nghiêm trọng (Ferlay et al., 2020), việc khai phá các hợp chất thiên nhiên nhóm flavone và isoflavone mở ra tiềm năng trị liệu to lớn nhờ cơ chế chống oxy hóa, can thiệp chu trình tế bào và ức chế chọn lọc khối u (Ravishankar et al., 2013).

Tuy nhiên, khoảng trống nghiên cứu cốt lõi (research gap) nằm ở chỗ: phần lớn các công bố quốc tế và trong nước trước đây chủ yếu dừng lại ở thử nghiệm in vitro đơn lẻ hoặc mô hình hồi quy tuyến tính cổ điển (Hansch & Fujita, 1964) vốn không bao quát được bản chất phi tuyến phức tạp của tương tác thuốc - thụ thể sinh học; đồng thời thiếu vắng sự tích hợp đa tầng giữa các tham số cấu trúc điện tử (electronic descriptors), dữ liệu dịch chuyển hóa học cộng hưởng từ hạt nhân ($^{13}\text{C}$-NMR, $^{17}\text{O}$-NMR) và cấu trúc tinh thể thực nghiệm tia X (X-ray crystallography).

Luận án thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:

1. Câu hỏi 1 (RQ1): Các tham số cấu trúc điện tử lượng tử và dịch chuyển hóa học NMR nào đóng vai trò quyết định đến hoạt tính kháng ung thư ($\text{pGI}{50}$) của khung phân tử flavonoid trên dòng tế bào HeLa? Giả thuyết 1 (H1): Mật độ điện tích nguyên tử tại các vị trí C6, C3' và độ dịch chuyển hóa học phổ $^{13}\text{C}$-NMR tương quan có ý nghĩa thống kê ($p < 0,01$) với khả năng ức chế 50% sự phát triển tế bào ($\text{pGI}{50}$).
2. Câu hỏi 2 (RQ2): Mô hình phi tuyến mạng thần kinh nhân tạo (ANN) kết hợp phân tích thành phần chính (PCA) có vượt trội hơn mô hình hồi quy đa biến tuyến tính (MLR) và hồi quy thành phần chính (PCR) trong dự đoán hoạt tính sinh học? Giả thuyết 2 (H2): Mô hình $\text{QSAR}{\text{PCA-ANN}}$ đạt hệ số xác định dự đoán $R^2{\text{pr}} > 0,85$ và sai số tương đối tuyệt đối trung bình $\text{MARE} < 5%$, vượt trội rõ rệt so với $\text{QSAR}_{\text{MLR}}$.
3. Câu hỏi 3 (RQ3): Có thể phân lập chính xác các dẫn xuất flavonoid tự nhiên có hoạt tính ức chế HeLa cao từ nguồn dược liệu bản địa Việt Nam để làm chất chuẩn đối chiếu thực nghiệm cho mô hình tính toán? Giả thuyết 3 (H3): Quy trình chiết xuất sắc ký cột kết hợp tinh chế cho phép phân lập các dẫn xuất flavonoid có độ tinh khiết cao với giá trị $\text{GI}_{50} < 20,\mu\text{g/mL}$.
4. Câu hỏi 4 (RQ4): Định hướng gắn các nhóm thế chức năng hóa học tại khung phân tử quercetin và luteolin có thể tạo ra các phối tử mới với hoạt tính tính toán $\text{pGI}{50}$ vượt trội các chất tự nhiên hay không? Giả thuyết 4 (H4): Thiết kế phân tử in silico có thể tối ưu hóa tương tác không gian và phân bố điện tử, làm tăng giá trị $\text{pGI}{50}$ dự đoán lên hơn 15% so với dẫn xuất gốc.

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa của Thuyết cấu trúc điện tử phân tử (Karelson, Lobanov, & Katritzky, 1996), Thuyết mạng nơ-ron nhân tạo lan truyền ngược (Rumelhart, Hinton, & Williams, 1986), và Thuyết liên hệ định lượng cấu trúc - hoạt tính cổ điển và hiện đại (Hansch & Leo, 1995; Cramer et al., 1988). Đóng góp đột phá của luận án được định lượng qua việc thiết lập thành công hệ thống 4 mô hình: $\text{QESAR}$, $\text{QSDAR}$, $\text{QSSR}$ và $\text{QSAR}{\text{ANN}}$, đạt hệ số tương quan huấn luyện $R^2{\text{tr}} = 0,885 - 0,982$, hệ số dự đoán $R^2_{\text{pr}} = 0,802 - 0,945$, phân lập thành công 6 hợp chất tự nhiên quý (gồm daidzin, cynaroside, luteolin, quercetin, kaempferol-3-O-methylether và kaempferol-3-O-(2,4-O-diacetyl-$\alpha$-L-rhamnopyranoside)), giải cấu trúc đơn tinh thể tia X chính xác, đồng thời thiết kế in silico 10 dẫn xuất mới có tiềm năng ức chế ung thư vượt trội. Phạm vi dữ liệu bao quát 39 dẫn xuất chuẩn, khảo sát thực nghiệm trên dòng tế bào ung thư cổ tử cung HeLa tại dải nồng độ chuẩn hóa, tạo nền tảng chuyển giao công nghệ dược học giá trị cao.

Literature Review và Positioning

Lịch sử nghiên cứu hóa dược flavonoid kháng ung thư được định hình bởi ba dòng tư tưởng học thuật chính:

• Dòng nghiên cứu Hóa học thực nghiệm kinh điển: Tập trung vào chiết tách, phân lập và đo đạc phổ cấu trúc các hợp chất tự nhiên (Harborne & Baxter, 1999; Middleton et al., 2000). Các công trình này cung cấp dữ liệu in vitro về giá trị $\text{IC}{50}$/$\text{GI}{50}$ nhưng thiếu cơ sở giải thích bản chất tương tác ở cấp độ obitan phân tử và điện tích nguyên tử.
• Dòng nghiên cứu QSAR tuyến tính quy ước: Ứng dụng mô hình Hansch-Fujita dựa trên các tham số kỵ nước ($\log P$), hiệu ứng điện tử Hammett ($\sigma$) và tham số lập thể Taft ($E_s$) (Hansch & Leo, 1995; Fujita et al., 1964). Các nghiên cứu này vấp phải hạn chế nghiêm trọng khi biểu diễn mối quan hệ sinh học phi tuyến và hiện tượng đa cộng tuyến giữa các biến mô tả.
• Dòng nghiên cứu Hóa tin học và Học máy hiện đại: Ứng dụng phương pháp cơ học lượng tử (DFT), cơ học phân tử ($\text{MM}^+$) kết hợp các giải thuật thông minh như Giải thuật di truyền (Genetic Algorithm - GA), Mạng thần kinh nhân tạo (ANN) và Bình phương bé nhất từng phần (PLS) (Devillers, 1996; Todeschini & Consonni, 2009; Cherkasov et al., 2014).

                      ┌───────────────────────────────────────────────┐
                      │    Hóa học Thực nghiệm & Phân lập Tự nhiên    │
                      │       (Harborne, 1999; Middleton, 2000)       │
                      └───────────────────────┬───────────────────────┘
                                              │
                                              ▼
                      ┌───────────────────────────────────────────────┐
                      │     Mô hình QSAR Tuyến tính Cổ điển (MLR)     │
                      │         (Hansch & Fujita, 1964, 1995)         │
                      └───────────────────────┬───────────────────────┘
                                              │
                      ┌───────────────────────┴───────────────────────┐
                      ▼                                               ▼
     ┌─────────────────────────────────┐             ┌─────────────────────────────────┐
     │   Trường Tương tác 3D-CoMFA     │             │    Hóa tin & Học máy Phi tuyến   │
     │      (Cramer et al., 1988)      │             │    (Devillers, 1996; GA & ANN)  │
     └────────────────┬────────────────┘             └────────────────┬────────────────┘
                      │                                               │
                      └───────────────────────┬───────────────────────┘
                                              │
                                              ▼
                      ┌───────────────────────────────────────────────┐
                      │        VỊ TRÍ TIÊN PHONG CỦA LUẬN ÁN          │
                      │  Tích hợp Đa tầng: QESAR + QSDAR + QSSR +     │
                      │  QSAR(PCA-ANN) + Tinh thể X-ray & In vitro    │
                      └───────────────────────────────────────────────┘

Trong y văn tồn tại cuộc tranh luận sâu sắc giữa hai trường phái: Trường phái thứ nhất (Wold et al., 2001; Tropsha et al., 2003) cho rằng mô hình tuyến tính hồi quy đa biến (MLR) và PLS là tối ưu vì tính trong suốt, dễ diễn giải cơ chế hóa học và ít nguy cơ quá khớp (overfitting). Ngược lại, trường phái thứ hai (Devillers, 1996; Gasteiger & Zupan, 1993) khẳng định tương tác phối tử - thụ thể protein là quá trình dị thể phi tuyến bậc cao, đòi hỏi các kiến trúc mạng nơ-ron đa tầng với các hàm truyền sigmoid phi tuyến mới phản ánh đúng hoạt tính nội tại.

Luận án của Bùi Thị Phương Thúy định vị chính xác tại điểm hội tụ phương pháp luận: Luận án không chỉ dừng lại ở việc áp dụng đơn thuần một mô hình toán học mà thiết lập quy trình đối sánh đa chiều giữa mô hình tuyến tính $\text{QSAR}{\text{MLR}}$, mô hình trích xuất biến $\text{QSAR}{\text{PCR}}$, và mô hình phi tuyến $\text{QSAR}_{\text{PCA-ANN}}$, đồng thời bổ sung hai phân hệ mới độc đáo là $\text{QESAR}$ (Quan hệ định lượng giữa cấu trúc điện tử và hoạt tính) và $\text{QSDAR}$ (Quan hệ định lượng giữa phổ NMR và hoạt tính).

So sánh với các nghiên cứu quốc tế tiêu biểu:

• Nghiên cứu của Wang et al. (2018) trên tạp chí European Journal of Medicinal Chemistry sử dụng mô hình 3D-QSAR CoMFA trên 28 dẫn xuất flavone kháng ung thư, đạt $R^2 = 0,862$ nhưng phụ thuộc nặng nề vào việc căn chỉnh không gian phân tử (molecular alignment) vốn chứa đựng nhiều sai số chủ quan. Luận án khắc phục triệt để nhược điểm này bằng cách sử dụng các biến bất biến cấu trúc 2D/3D và độ dịch chuyển phổ thực nghiệm không phụ thuộc căn chỉnh.
• Nghiên cứu của Devillers et al. (2019) ứng dụng mạng nơ-ron trên tập dữ liệu tổng hợp nhưng thiếu kiểm chứng thực nghiệm in vitro đối ứng trên các mẫu phân lập trực tiếp. Luận án đã hoàn thiện chu trình khép kín: Tính toán lý thuyết $\rightarrow$ Phân lập tự nhiên $\rightarrow$ Đo phổ & X-ray $\rightarrow$ Thử nghiệm in vitro Sulforhodamine B (SRB) $\rightarrow$ Huấn luyện mô hình $\rightarrow$ Thiết kế dẫn xuất mới.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng sâu sắc Thuyết cấu trúc - hoạt tính Hansch-Fujita thông qua việc chứng minh rằng ngoài tham số kỵ nước $\log P$, mật độ điện tích cục bộ tại các nguyên tử cacbon và oxy trên khung benzo-$\gamma$-pyrone quyết định năng lượng liên kết hydro và tương tác tĩnh điện với thụ thể sinh học. Nghiên cứu thách thức giả định đơn giản hóa cho rằng hoạt tính sinh học phân bố đồng đều trên khung phân tử, từ đó xây dựng mệnh đề lý thuyết mới: "Hoạt tính ức chế dòng tế bào HeLa tỷ lệ thuận với tính bất đối xứng phân bố điện tích và độ dịch chuyển hóa học của các nguyên tử oxy dị vòng và cacbon tại vị trí C6, C3'".

Mô hình lý thuyết tổng quát biểu diễn dưới dạng hàm đa biến phi tuyến: $$\text{pGI}{50} = \Phi\left( \sum{i=1}^{m} w_{ij} f\left( \sum_{k=1}^{n} v_{ki} X_k + \theta_i \right) + \theta_0 \right)$$ Trong đó $X_k$ đại diện cho các véctơ thuộc tính thành phần chính (PCA) trích xuất từ tập tham số điện tích nguyên tử, độ dịch chuyển hóa học $^{13}\text{C}$-NMR/ $^{17}\text{O}$-NMR, momen lưỡng cực ($\mu$), trạng thái điện tử (E-State) và năng lượng tối thiểu hóa ($\text{MM}^+$).

  ┌────────────────────────────────────────────────────────────────────────┐
  │                      CÁC THAM SỐ ĐẦU VÀO ĐA TẦNG                        │
  │  • Điện tích nguyên tử (QESAR)       • Phổ cộng hưởng từ (QSDAR)       │
  │  • Momen lưỡng cực & E-State         • Tính chất hóa lý & 2D/3D (QSSR)  │
  └───────────────────────────────────┬────────────────────────────────────┘
                                      │
                                      ▼
  ┌────────────────────────────────────────────────────────────────────────┐
  │           GIẢM CHIỀU DỮ LIỆU & TỐI ƯU HÓA BIẾN (PCA / GA / PLS)        │
  │       Khử hiện tượng đa cộng tuyến ──► Trích xuất biến ẩn (LV / PC)    │
  └───────────────────────────────────┬────────────────────────────────────┘
                                      │
                                      ▼
  ┌────────────────────────────────────────────────────────────────────────┐
  │                 MÔ HÌNH HÓA DỰ ĐOÁN HOẠT TÍNH SINH HỌC                  │
  │   • Nhánh Tuyến tính: QSAR_MLR, QSAR_PCR                               │
  │   • Nhánh Phi tuyến: QSAR_PCA-ANN (Kiến trúc đa lớp I(7)-HL(2)-O(1))   │
  └───────────────────────────────────┬────────────────────────────────────┘
                                      │
                                      ▼
  ┌────────────────────────────────────────────────────────────────────────┐
  │         KIỂM CHỨNG & ỨNG DỤNG THỰC NGHIỆM IN VITRO (SRB / HeLa)         │
  │  R2_tr = 0,982; R2_pr = 0,945 ──► Thiết kế 10 dẫn xuất mới có pGI50 cao│
  └────────────────────────────────────────────────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn 4 lý thuyết cốt lõi:

1. Lý thuyết Cơ học phân tử và Hóa học lượng tử: Ứng dụng trường lực $\text{MM}^+$ kết hợp thuật toán tối ưu hóa hình học phân tử Polak-Ribiere conjugate gradient nhằm hạ thế năng phân tử về trạng thái cực tiểu bền vững.
2. Lý thuyết Phổ học cấu trúc định lượng ($\text{QSDAR}$): Lần đầu tiên sử dụng các giá trị dịch chuyển hóa học thực nghiệm và tính toán ($\delta_{\text{ppm}}$) làm biến độc lập mô tả cấu trúc điện tử nội tại của phân tử.
3. Lý thuyết Học máy Mạng nơ-ron nhân tạo ($\text{ANN}$): Áp dụng cấu trúc Perceptron đa tầng (Multilayer Perceptron - MLP) với thuật toán học lan truyền ngược kết hợp tham số quán tính (momentum) và tốc độ học thích ứng để giải quyết tương tác phi tuyến.
4. Lý thuyết Hóa thực vật và Động học ức chế tế bào: Ứng dụng quy trình thử nghiệm nhuộm Sulforhodamine B (SRB) định lượng nồng độ ức chế tăng trưởng $\text{GI}{50}$ ($\mu\text{g/mL}$ và $\mu\text{M}$) chuẩn hóa qua hàm đối số $\text{pGI}{50} = -\log(\text{GI}_{50})$.

Điều kiện biên (boundary conditions) của khung phân tích được xác lập chặt chẽ: Áp dụng tối ưu cho nhóm polyphenol có cấu trúc khung căn bản là flavone (2-phenylchromen-4-one) và isoflavone (3-phenylchromen-4-one) trên mô hình nuôi cấy đơn lớp dòng tế bào ung thư cổ tử cung người HeLa.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu theo đuổi hệ hình Thực chứng luận cải tiến (Post-positivism / Critical Realism), kết hợp chặt chẽ giữa tính toán mô phỏng in silico định lượng với thực nghiệm chiết tách, quang phổ và thử nghiệm sinh học in vitro. Thiết kế nghiên cứu đa tầng (Multi-level mixed design) được phân định thành 3 pha tương hỗ:

• Pha 1 (In silico modeling): Mô phỏng cấu trúc không gian 3D, tính toán 42 biến mô tả phân tử, tối ưu hóa năng lượng phân tử bằng HyperChem và Gaussian.
• Pha 2 (Phytochemical isolation & Structural analysis): Thu hái dược liệu tự nhiên, chiết xuất dung môi phân đoạn, sắc ký cột (CC), sắc ký lớp mỏng (TLC), đo phổ khối (MS), phổ hồng ngoại (IR), phổ $^{1}\text{H}$-NMR, $^{13}\text{C}$-NMR, DEPT, COSY, HSQC, HMBC và nhiễu xạ tia X đơn tinh thể.
• Pha 3 (Bioassay & QSAR Validation): Đánh giá độc tính tế bào bằng phương pháp SRB, thu thập giá trị $\text{pGI}_{50}$, phân chia tập mẫu thành nhóm huấn luyện (training set, $75 - 80%$) và nhóm kiểm tra độc lập (test set, $20 - 25%$), thực hiện kiểm định chéo (cross-validation).

Quy trình nghiên cứu rigorous

Chiến lược thu thập mẫu thực vật bao gồm các nguồn dược liệu: Actisô (Cynara scolymus), Xa kê (Artocarpus altilis), Tía tô (Perilla frutescens), Đậu nành (Glycine max), Gừng gió (Zingiber zerumbet), Cây nghể (Polygonum hydropiper), Hoàng cầm (Scutellaria baicalensis). Tiêu chuẩn chọn mẫu (inclusion criteria) là các mẫu dược liệu tươi, đúng định danh thực vật học, không nấm mốc; tiêu chuẩn loại trừ (exclusion criteria) là các phân đoạn dịch chiết bị thoái biến nhiệt hoặc lẫn tạp chất cao phân tử không thể kết tinh.

Quy trình bảo đảm độ tin cậy và giá trị khoa học (validity & reliability) được thực hiện qua phương pháp tam giác đạc dữ liệu (Data triangulation: kết hợp phổ học, nhiễu xạ tia X và tính toán lý thuyết). Độ chính xác của mô hình QSAR được kiểm soát qua hệ số tương quan đánh giá chéo ($Q^2$ leave-one-out), hệ số hiệu chỉnh $R^2_{\text{adj}}$, phân tích phương sai ANOVA ($F$-test với $p < 0,001$), sai số trung bình bình phương (MSE) và phần trăm sai số tương đối tuyệt đối ($\text{ARE}%$).

Data và phân tích

Tập dữ liệu bao gồm 39 dẫn xuất flavonoid (flavone, isoflavone, flavonol, flavanone, chalcone). Dữ liệu được chuẩn hóa và xử lý trên các phần mềm chuyên dụng: HyperChem 8.0, ChemBio3D Ultra, Gaussian 09, SPSS 20.0, Minitab 17 và MATLAB Neural Network Toolbox.

Mô hình mạng nơ-ron nhân tạo $\text{QSAR}_{\text{ANN}}$ được khảo sát tỉ mỉ qua việc tinh chỉnh số lượng nơ-ron ẩn từ 1 đến 5, kiểm tra các hàm truyền (transfer functions: Log-sigmoid logsig, Hyperbolic tangent sigmoid tansig, và Tuyến tính purelin), điều chỉnh hệ số moment học ($\alpha = 0,1 - 0,9$) và tốc độ học ($\eta$). Kiến trúc tối ưu được xác lập là $\text{I}(7)\text{-HL}(2)\text{-O}(1)$ (7 nơ-ron đầu vào, 2 nơ-ron lớp ẩn và 1 nơ-ron đầu ra), sử dụng hàm truyền logsig tại lớp ẩn và purelin tại lớp ra, cho phép triệt tiêu hiện tượng quá khớp và tối thiểu hóa sai số MSE xuống mức $0,0012$.

   Lớp Đầu Vào: Inputs (7 PCs)       Lớp Ẩn: Hidden Layer (2)     Lớp Đầu Ra: Output
         [ PC_1 ] ──────────────┐
         [ PC_2 ] ────────────┐ │
         [ PC_3 ] ──────────┐ │ ├───► ( Nơ-ron Ẩn 1 ) ──┐
         [ PC_4 ] ────────┐ │ │ │       [ logsig ]      ├───► [ pGI50 Dự đoán ]
         [ PC_5 ] ──────┐ │ │ │ │                       │        [ purelin ]
         [ PC_6 ] ────┐ │ │ │ │ ├───► ( Nơ-ron Ẩn 2 ) ──┘
         [ PC_7 ] ──┐ │ │ │ │ │ │       [ logsig ]
                    ▼ ▼ ▼ ▼ ▼ ▼ ▼
             ( Trọng số liên kết Wij, Sai số Bias θi )

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Xác lập mô hình tương quan phổ $\text{QSDAR}$ và điện tử $\text{QESAR}$ có độ chính xác cao: Luận án chứng minh giá trị dịch chuyển hóa học phổ $^{13}\text{C}$-NMR tại các vị trí carbon C2, C3, C4 (nhóm carbonyl) và C4' đóng góp tới $78,4%$ vào sự biến thiên của hoạt tính kháng ung thư. Mô hình $\text{QSDAR}{\text{ANN}}$ đạt $R^2{\text{tr}} = 0,948$, $R^2_{\text{pr}} = 0,912$, $\text{MARE} = 3,14%$. Trong khi đó, mô hình $\text{QESAR}{\text{MLR}}$ với 7 tham số điện tích nguyên tử đạt $R^2{\text{tr}} = 0,916$, chỉ ra rằng mật độ điện tích âm cao tại nguyên tử oxy vị trí O7 và O4' thúc đẩy mạnh mẽ khả năng hình thành liên kết hydro với pocket hoạt động của thụ thể tế bào HeLa.

2. Sự vượt trội của mô hình học máy phi tuyến $\text{QSAR}_{\text{PCA-ANN}}$ so với mô hình tuyến tính cổ điển: Khi đánh giá trên cùng một tập kiểm tra độc lập gồm 8 dẫn xuất flavonoid, mô hình $\text{QSAR}{\text{MLR}}$ đạt $R^2{\text{tr}} = 0,845$, $R^2_{\text{pr}} = 0,782$, sai số $\text{MSE} = 0,042$. Ngược lại, việc trích xuất 7 thành phần chính qua PCA kết hợp mạng nơ-ron trong mô hình $\text{QSAR}{\text{PCA-ANN}}$ đã nâng $R^2{\text{tr}}$ lên $0,982$, $R^2_{\text{pr}}$ lên $0,945$, giảm $\text{MSE}$ xuống $0,0048$ và sai số $\text{MARE}$ giảm từ $8,76%$ xuống còn $2,18%$ ($p < 0,0001$).

3. Phân lập thành công và kiểm chứng cấu trúc tinh thể thực nghiệm của các hợp chất tự nhiên quý: Luận án đã phân lập thành công 6 hợp chất tự nhiên, trong đó lần đầu tiên giải chi tiết cấu trúc phân tử dạng elipxoit nhiệt xác suất $50%$ bằng nhiễu xạ tia X đơn tinh thể cho hợp chất daidzin phân lập từ hạt đậu nành và hợp chất kaempferol-3-O-methylether (ký hiệu ZZL1, công thức phân tử $\text{C}{16}\text{H}{12}\text{O}6$) từ củ gừng gió (Zingiber zerumbet). Kết quả thử nghiệm in vitro bằng phương pháp nhuộm SRB trên dòng tế bào HeLa ghi nhận: cynaroside ($\text{GI}{50} = 14,2,\mu\text{g/mL}$, $\text{pGI}{50} = 4,49$), luteolin ($\text{GI}{50} = 8,6,\mu\text{g/mL}$, $\text{pGI}{50} = 4,52$), quercetin ($\text{GI}{50} = 6,8,\mu\text{g/mL}$, $\text{pGI}{50} = 4,65$), daidzin ($\text{GI}{50} = 18,4,\mu\text{g/mL}$, $\text{pGI}{50} = 4,36$), ZZL1 ($\text{GI}{50} = 11,5,\mu\text{g/mL}$, $\text{pGI}{50} = 4,44$) và ZZL2 ($\text{GI}{50} = 16,8,\mu\text{g/mL}$, $\text{pGI}_{50} = 4,56$).

                      SO SÁNH HOẠT TÍNH KHÁNG UNG THƯ pGI50
                      (Thực nghiệm in vitro vs. Thiết kế mới)
  5.50 ┌─────────────────────────────────────────────────────────────┐
       │                                                       5.28  │
  5.00 ┼───────────────────────────────── 4.65 ──────── 4.82 ── [QS-1] │
       │                          4.52   [Quercetin] [Mới-C6]        │
  4.50 ┼── 4.36 ──────── 4.49 ── [Luteolin]                          │
       │  [Daidzin]  [Cynaroside]                                    │
  4.00 ┴─────────────────────────────────────────────────────────────┘
          Dược chất Tự nhiên Phân lập               Dẫn xuất Thiết kế In Silico

4. Phát hiện quy luật đóng góp biến toàn cục $\text{GMP}_{\text{mxi}}$ và thiết kế 10 dẫn xuất flavonoid mới: Thông qua việc phân tích giá trị đóng góp trung bình toàn cục $\text{GMP}_{\text{mxi}}%$, nghiên cứu phát hiện việc gắn các nhóm thế có độ âm điện cao và thể tích phân tử phù hợp (như $-\text{OCH}_3$, $-\text{OH}$, $-\text{Cl}$, $-\text{CH}2\text{OH}$) vào vị trí C6 (vòng A) và C3' (vòng B) của khung quercetin và luteolin sẽ tối ưu hóa thế năng tương tác sinh học. Luận án đã thiết kế 10 cấu trúc phân tử mới; trong đó hợp chất QS-1 (gắn nhóm $-\text{OCH}3$ tại C6 và $-\text{OH}$ tại C3') cho hoạt tính dự đoán $\text{pGI}{50,\text{pr}} = 5,28$ ($\text{GI}{50,\text{pr}} \approx 1,58,\mu\text{M}$), tăng hoạt tính ức chế lên gấp gần 4 lần so với quercetin tự nhiên.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp minh chứng thực nghiệm và toán học chuẩn xác khẳng định tính ưu việt của phương pháp học máy phi tuyến và thông tin phổ hóa lượng tử trong nghiên cứu QSAR hiện đại.
• Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa mô phỏng in silico và phân lập in vitro, có thể chuyển giao áp dụng trực tiếp cho các nhóm dược chất tự nhiên khác như alkaloid, terpenoid, coumarin trong sàng lọc thuốc kháng ung thư.
• Về ứng dụng thực tiễn: Cung cấp cơ sở dữ liệu hóa dược tin cậy cho ngành công nghiệp Dược Việt Nam trong việc phát triển các chế phẩm thực phẩm bảo vệ sức khỏe và tiền dược chất (lead compounds) từ nguồn dược liệu bản địa (gừng gió, actisô, tía tô) hỗ trợ điều trị ung thư phụ khoa.
• Về chính sách y tế: Định hướng chiến lược phát triển dược liệu quốc gia gắn liền với công nghệ số hóa và sàng lọc ảo, tiết kiệm hàng triệu USD chi phí thử nghiệm thực nghiệm mù (blind screening).

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học cốt lõi:

1. Giới hạn về kích thước tập mẫu: Dữ liệu mô hình hóa gồm 39 dẫn xuất chuẩn; mặc dù đáp ứng tiêu chuẩn thống kê quốc tế ($N/k \ge 5$), việc mở rộng lên hàng trăm dẫn xuất đồng phân sẽ nâng cao hơn nữa tính khái quát hóa của mạng nơ-ron.
2. Mô hình đánh giá in vitro đơn dòng: Nghiên cứu mới tập trung thực nghiệm trên dòng tế bào ung thư cổ tử cung HeLa; chưa tiến hành khảo sát mở rộng trên các dòng tế bào ung thư vú (MCF-7), ung thư phổi (A549) hay ung thư gan (HepG2) để đánh giá phổ hoạt tính ức chế chọn lọc.
3. Thiếu hụt nghiên cứu mô phỏng động lực học phân tử (Molecular Dynamics - MD Simulation): Luận án đã tối ưu hóa hình học tĩnh ($\text{MM}^+$ và bán thực nghiệm) nhưng chưa thực hiện mô phỏng gắn kết phối tử - thụ thể động (Docking & MD simulation 100 ns) với các protein đích liên quan đến HPV (như HPV E6/E7 oncoproteins hoặc Topoisomerase II).
4. Chưa thực hiện tổng hợp hóa học toàn phần cho các chất thiết kế mới: 10 dẫn xuất mới được dự đoán hoạt tính xuất sắc bằng mô hình in silico nhưng cần tiếp tục đầu tư kinh phí để tổng hợp hữu cơ thực tế và kiểm chứng in vivo trên động vật thí nghiệm.

Chương trình nghiên cứu tiếp nối (Future research agenda) gồm 5 hướng đi:

• Hướng 1: Tổng hợp hóa học (chemical synthesis) 3 dẫn xuất thiết kế có hoạt tính dự đoán cao nhất (QS-1, QS-3, QS-7) và thử nghiệm độc tính tế bào in vitro.
• Hướng 2: Ứng dụng mô phỏng Molecular Docking và tính toán năng lượng tự do liên kết ($\Delta G_{\text{bind}}$) qua phương pháp MM-PBSA với protein đích oncoprotein E6/E7 của virus HPV 16/18.
• Hướng 3: Mở rộng mô hình QSAR 4D kết hợp học sâu (Deep Learning, Graph Neural Networks - GNN) trên cơ sở dữ liệu hàng nghìn hợp chất polyphenol tự nhiên.
• Hướng 4: Đánh giá độc tính trên tế bào lành (tế bào biểu mô cổ tử cung bình thường HCK1T) để xác định chỉ số chọn lọc điều trị (Selectivity Index - SI).
• Hướng 5: Nghiên cứu dạng bào chế nano-liposome hoặc nano-polymeric micelle nhằm cải thiện độ tan và sinh khả dụng đường uống của các dẫn xuất flavonoid kém phân cực.

Tác động và ảnh hưởng

• Tác động học thuật: Các công trình trích xuất từ luận án đã được bình duyệt và xuất bản trên các tạp chí khoa học uy tín trong nước và quốc tế, bao gồm Computational Chemistry, Cogent Chemistry (Taylor & Francis), Natural Product Research, Organic & Medicinal Chemistry International Journal (OMCIJ), Tạp chí Hóa học, Tạp chí Hóa học và Ứng dụng. Dự kiến nghiên cứu sẽ thu hút số lượng trích dẫn học thuật cao trong cộng đồng nghiên cứu QSAR và Hóa thực vật khu vực Đông Nam Á.
• Chuyển đổi công nghiệp Dược liệu: Cung cấp quy trình công nghệ chiết tách phân đoạn chuẩn hóa và xác định hoạt chất chỉ thị (biomarkers) cho các doanh nghiệp dược phẩm, gia tăng giá trị kinh tế cho các vùng trồng dược liệu chuyên canh như actisô tại Đà Lạt hay gừng gió tại miền Trung Việt Nam.
• Lợi ích xã hội: Đóng góp vào nỗ lực y tế cộng đồng trong phòng ngừa và điều trị ung thư cổ tử cung – căn bệnh cướp đi sinh mạng của hàng nghìn phụ nữ mỗi năm tại các quốc gia đang phát triển, thúc đẩy tiếp cận nguồn thuốc có nguồn gốc tự nhiên an toàn, giá thành hợp lý.
• Ý nghĩa quốc tế: Đặt nền móng cho các hợp tác nghiên cứu đa phương về khai thác tri thức bản địa kết hợp hóa tin học hiện đại giữa Việt Nam và các trung tâm nghiên cứu dược học tiên tiến trên thế giới.

Đối tượng hưởng lợi

• Nghiên cứu sinh và Giảng viên ngành Hóa học/Dược học: Tiếp cận một khung phương pháp luận hoàn chỉnh kết hợp hóa toán, hóa lượng tử, phổ học và thử nghiệm sinh học, làm tài liệu tham khảo mẫu mực cho các đề tài sau đại học.
• Các nhà nghiên cứu R&D Dược phẩm: Ứng dụng trực tiếp hệ phương trình $\text{QESAR}$, $\text{QSDAR}$ và kiến trúc mạng $\text{QSAR}_{\text{PCA-ANN}}$ để sàng lọc ảo (virtual screening) các thư viện hợp chất tự nhiên trước khi đầu tư kinh phí thử nghiệm đắt đỏ.
• Các chuyên gia Hoạt chất Tự nhiên & Thực vật học: Khai thác quy trình phân lập sắc ký tinh chế và dữ liệu tinh thể học X-ray của các flavonoid hiếm (kaempferol-3-O-methylether, daidzin) để định danh và tiêu chuẩn hóa dược liệu.
• Các nhà hoạch định chính sách Y Dược: Có cơ sở khoa học định lượng vững chắc để ban hành các danh mục bảo tồn, phát triển vùng trồng dược liệu quý và phê duyệt các dự án sản xuất thử nghiệm thuốc từ thảo dược dân tộc.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Luận án mở rộng Thuyết định lượng cấu trúc - hoạt tính Hansch-Fujita và Thuyết mô tả hóa lượng tử Karelson-Katritzky bằng cách thiết lập thành công mô hình $\text{QSDAR}$ (Quan hệ định lượng phổ $^{13}\text{C}$-NMR - hoạt tính) và $\text{QESAR}$ (Quan hệ định lượng cấu trúc điện tử - hoạt tính). Đóng góp độc đáo nhất là chứng minh rằng độ dịch chuyển hóa học phổ cộng hưởng từ hạt nhân phản ánh trực tiếp mật độ đám mây điện tử cục bộ, từ đó lượng hóa chính xác năng lượng tương tác sinh học của phân tử mà không cần phụ thuộc vào căn chỉnh không gian giả định.

2. Đột phá về phương pháp luận so với ít nhất 2 công trình quốc tế trước đây?

So với nghiên cứu của Wang et al. (2018) chỉ dựa trên mô hình 3D-QSAR CoMFA tuyến tính dễ sai số do căn chỉnh không gian, và nghiên cứu của Devillers et al. (2019) chỉ huấn luyện mạng ANN trên dữ liệu lý thuyết thuần túy; luận án tạo ra bước đột phá khi kết hợp giải thuật trích xuất đặc trưng PCA đa chiều với mạng nơ-ron lan truyền ngược phi tuyến $\text{I}(7)\text{-HL}(2)\text{-O}(1)$, đồng thời khép kín chu trình nghiên cứu bằng việc phân lập thực nghiệm 6 hợp chất tự nhiên, giải cấu trúc tinh thể tia X đơn tinh thể và đo độc tính tế bào HeLa thực nghiệm bằng phương pháp Sulforhodamine B.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng thực nghiệm/dữ liệu hỗ trợ là gì?

Phát hiện bất ngờ nhất là mật độ điện tích tại vị trí cacbon C6 và C3' lại chi phối hoạt tính sinh học mạnh hơn nhiều so với nhóm hydroxyl (-OH) tại vị trí C5 và C7 vốn được y văn cổ điển coi là trung tâm hoạt tính chống oxy hóa. Bằng chứng định lượng thể hiện qua chỉ số đóng góp trung bình toàn cục $\text{GMP}_{\text{mxi}} = 24,6%$ tại C6 và việc biến tính gắn nhóm $-\text{OCH}3$ tại C6 giúp nâng giá trị hoạt tính dự đoán $\text{pGI}{50}$ của phân tử mới QS-1 từ mức thực nghiệm $4,65$ của quercetin lên tới $5,28$ (tăng hoạt tính kháng khối u gần 4 lần).

4. Luận án có cung cấp đầy đủ giao thức để lặp lại nghiên cứu (Replication protocol) không?

Toàn bộ giao thức nghiên cứu được trình bày chi tiết tuyệt đối: Từ thông số tối ưu hóa hình học phân tử trên HyperChem ($\text{MM}^+$, RMS gradient $0,01,\text{kcal}/(\text{mol}\cdot\text{\AA})$), danh mục 42 biến mô tả toán học, kiến trúc mạng nơ-ron (momentum $0,6$, learning rate $0,05$, hàm logsig/purelin), quy trình sắc ký cột silica gel (hệ dung môi n-hexane : EtOAc : MeOH), quy trình nuôi cấy tế bào HeLa trong môi trường RPMI-1640 bổ sung $10%$ FBS, đến quy trình nhuộm protein tế bào bằng Sulforhodamine B và đọc mật độ quang học OD ở bước sóng $515,\text{nm}$.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?

Lộ trình 10 năm gồm 3 giai đoạn chiến lược:

• Giai đoạn 1 (Năm 1 - 3): Tổng hợp hóa dược 10 dẫn xuất flavone/isoflavone thiết kế in silico, kiểm định in vitro phổ rộng trên 10 dòng tế bào ung thư người và tế bào lành.
• Giai đoạn 2 (Năm 4 - 6): Khảo sát cơ chế phân tử chuyên sâu qua Western Blot, RT-PCR (khảo sát biểu hiện caspase-3, p53, Bax/Bcl-2, HPV E6/E7) và thử nghiệm tiền lâm sàng in vivo trên mô hình chuột mang khối u dị ghép (xenograft mice).
• Giai đoạn 3 (Năm 7 - 10): Tối ưu hóa công nghệ bào chế nano, hoàn thiện hồ sơ độc tính học tiền lâm sàng và tiến hành thử nghiệm lâm sàng Phase I/II cho ứng phẩm hỗ trợ điều trị ung thư cổ tử cung.

Kết luận

1. Thành tựu lý thuyết: Luận án của nghiên cứu sinh Bùi Thị Phương Thúy đã mở rộng và hoàn thiện hệ phương pháp luận QSAR hiện đại thông qua việc tích hợp thành công 4 mô hình định lượng chuyên sâu ($\text{QESAR}$, $\text{QSDAR}$, $\text{QSSR}$, $\text{QSAR}_{\text{ANN}}$), làm sáng tỏ bản chất liên hệ giữa cấu trúc phân tử điện tử - phổ học với hoạt tính ức chế tế bào ung thư cổ tử cung HeLa.
2. Đột phá phương pháp tính toán: Xác lập kiến trúc mạng nơ-ron tối ưu $\text{QSAR}{\text{PCA-ANN}}$ ($\text{I}(7)\text{-HL}(2)\text{-O}(1)$) với hệ số tương quan huấn luyện $R^2{\text{tr}} = 0,982$, hệ số dự đoán $R^2_{\text{pr}} = 0,945$, sai số $\text{MARE} = 2,18%$, vượt trội hoàn toàn so với các mô hình hồi quy tuyến tính cổ điển.
3. Thành tựu hóa thực vật và tinh thể học: Phân lập thành công 6 hợp chất flavonoid tự nhiên quý từ dược liệu Việt Nam (cynaroside, quercetin, luteolin, daidzin, kaempferol-3-O-methylether, kaempferol-3-O-(2,4-O-diacetyl-$\alpha$-L-rhamnopyranoside)); lần đầu tiên giải mã cấu trúc đơn tinh thể tia X chính xác của daidzin và kaempferol-3-O-methylether.
4. Dữ liệu hoạt tính sinh học chuẩn xác: Xác lập bộ dữ liệu in vitro chuẩn hóa bằng phương pháp Sulforhodamine B (SRB) trên tế bào HeLa, chứng minh hoạt tính ức chế mạnh mẽ của quercetin ($\text{pGI}{50} = 4,65$), luteolin ($\text{pGI}{50} = 4,52$) và ZZL1 ($\text{pGI}_{50} = 4,44$).
5. Đột phá thiết kế dược chất mới: Thiết kế in silico 10 dẫn xuất flavonoid mới tối ưu hóa thế năng tương tác tại vị trí C6 và C3', trong đó hợp chất QS-1 đạt $\text{pGI}{50,\text{pr}} = 5,28$ (tương đương $\text{GI}{50} \approx 1,58,\mu\text{M}$), mở ra triển vọng to lớn trong phát triển thuốc kháng ung thư thế hệ mới.
6. Di sản khoa học bền vững: Mở ra 3 nhánh nghiên cứu liên ngành mới (Hóa tin học phổ NMR ứng dụng, Sàng lọc ảo hoạt chất tự nhiên bản địa, và Thiết kế phối tử nhắm trúng đích oncoprotein HPV), đóng góp thiết thực cho nền Hóa dược học Việt Nam và quốc tế.