Bối cảnh và vấn đề nghiên cứu

Trong bối cảnh số lượng thuốc mới được phát triển ngày càng giảm trong khi các bệnh nan y, đặc biệt là ung thư, không ngừng gia tăng, các phương pháp sàng lọc ảo (in silico) nổi lên như một chiến lược nghiên cứu hiệu quả, giúp tiết kiệm thời gian và chi phí trong quá trình phát triển thuốc. Theo số liệu từ GLOBOCAN năm 2020, Việt Nam ghi nhận 182.563 ca mắc mới và 122.690 ca tử vong do ung thư, xếp thứ 91/185 quốc gia về tỷ suất mắc mới và thứ 50/185 về tỷ suất tử vong trên 100.000 dân — tăng so với thứ hạng 99/185 và 56/185 của năm 2018. Dữ liệu này cho thấy tình hình ung thư tại Việt Nam đang diễn biến theo chiều hướng xấu và cần có thêm các hướng tiếp cận nghiên cứu thuốc mới.

Mục tiêu nghiên cứu của khóa luận được tác giả xác định gồm hai nhiệm vụ:

  1. Sử dụng docking phân tử để xác định các tương tác giữa các dẫn xuất 2'-hydroxychalconoid, flavanone và flavonol với protein Tubulin (4O2B) và enzyme ACLY (6O0H).
  2. Xác định các pharmacophore trên cấu trúc của các dẫn xuất này, từ đó đề xuất cấu trúc chất dẫn đầu có khả năng ức chế đồng thời hai mục tiêu Tubulin (4O2B) và ACLY (6O0H).

Đối tượng nghiên cứu gồm ba nhóm hợp chất flavonoid: 10 dẫn xuất 2'-hydroxychalconoid (HC1–HC10), 10 dẫn xuất flavanone (HC11–HC20) và 10 dẫn xuất flavonol (HC21–HC30), tổng cộng 30 hợp chất được tổng hợp từ nhóm nghiên cứu của PGS.TS Hoàng Minh Hảo, cùng hai protein đích là Tubulin (mã PDB: 4O2B) và enzyme ATP citrate lyase (mã PDB: 6O0H). Nghiên cứu được thực hiện hoàn toàn in silico, không giới hạn theo không gian địa lý, sử dụng dữ liệu protein từ ngân hàng dữ liệu PDB và cấu trúc hợp chất đã được xác nhận bằng phổ NMR và khối phổ (MS).


Cơ sở lý thuyết và phương pháp

Khung lý thuyết và khái niệm

Khóa luận xây dựng trên nền tảng lý thuyết về sàng lọc ảo (virtual screening — in silico), bao gồm hai hướng tiếp cận chính:

  • Sàng lọc dựa trên phối tử: Xây dựng mô hình pharmacophore hoặc QSAR (2D/3D) từ tập phối tử đã biết để dự đoán khả năng gắn kết của phân tử mới.
  • Sàng lọc dựa trên cấu trúc mục tiêu: Xác định khoang gắn kết của protein rồi áp dụng phương pháp docking để đánh giá tương tác phối tử–đích.

Nghiên cứu này thuộc hướng thứ hai. Nền tảng sinh học phân tử được sử dụng bao gồm:

  • Tubulin: Thành phần chính của vi ống (microtubule), gồm α-tubulin và β-tubulin (khối lượng ~55.000 Dalton, tương ứng 450 amino acid mỗi loại). β-tubulin là đích gắn kết của toàn bộ các thuốc liên kết tubulin hiện nay như taxane, vinca alkaloids và epothilone. Động học vi ống gồm hai kiểu: dynamic instability và treadmilling.
  • ACLY (ATP citrate lyase): Enzyme đồng nhất (homotetramer) gồm bốn chuỗi giống hệt nhau, khối lượng phân tử ~480 kD, mỗi chuỗi chứa 1.101 gốc amino acid. ACLY đóng vai trò liên kết giữa chuyển hóa glucose và tổng hợp acid béo thông qua sản xuất acetyl-CoA, được biểu hiện quá mức trong nhiều loại tế bào ung thư.
  • Flavonoid: Nhóm polyphenol có khung chung C6-C3-C6, trong đó chalcone là tiền chất của flavonol, flavanone, flavone và aurone.

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

Phương pháp chính là docking phân tử (molecular docking) — kỹ thuật mô hình hóa dự đoán vị trí và cấu dạng gắn kết của phân tử ligand vào khoang gắn kết của protein. Quy trình docking gồm ba bước: (1) đặt cấu dạng ligand vào khoang gắn kết bằng thuật toán tam giác phân, (2) chấm điểm gắn kết, (3) xếp hạng dựa trên năng lượng tự do (bao gồm lực tĩnh điện, Van der Waals, entropy và năng lượng solvat hóa). Phần mềm sử dụng là MOE 2008 (Molecular Operating Environment). Cấu trúc 2D các hợp chất được xây dựng bằng ChemBioDraw Ultra 12, sau đó chuyển đổi sang định dạng SMILES và tải vào MOE để tạo cấu trúc 3D.

Độ tin cậy của phương pháp được kiểm định qua re-docking (docking lặp lại): ligand đồng kết tinh được tách ra và docking lại vào khoang gắn kết; kết quả được đánh giá qua chỉ số RMSD (Root Mean Square Deviation) — nếu RMSD ≤ 2 Å thì phương pháp được xem là đáng tin cậy.

Nguồn dữ liệu

  • Cấu trúc tinh thể protein tải từ Protein Data Bank (PDB): 4O2B (phức Tubulin–Colchicine) và 6O0H (ATP citrate lyase).
  • 30 hợp chất flavonoid do nhóm nghiên cứu của PGS.TS Hoàng Minh Hảo tổng hợp, cấu trúc xác nhận bằng NMR và MS.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương này cung cấp nền tảng lý thuyết trên bốn mảng lớn. Thứ nhất, tổng quan về ung thư dựa trên số liệu GLOBOCAN 2020 — trong đó nhấn mạnh tình hình ung thư tại Việt Nam ngày càng nghiêm trọng hơn. Ung thư được định nghĩa là quá trình tế bào thoát khỏi kiểm soát sinh lý, xâm lấn và di căn đến các cơ quan khác. Tác giả phân biệt khối u lành tính và khối u ác tính, cũng như ba cơ chế gây tử vong do ung thư.

Thứ hai, phần về Tubulin trình bày siêu họ protein tubulin (α, β, γ, δ, ε, ζ-tubulin ở người), đặc tính hình thành vi ống qua dị thể (heterodimer) α/β-tubulin, hai kiểu động học dynamic instability và treadmilling của vi ống, và ứng dụng ức chế tubulin trong điều trị ung thư thông qua các thuốc taxane, vinca alkaloids và epothilone.

Thứ ba, phần về enzyme ACLY phân tích cấu trúc homotetramer, vai trò trung gian giữa chuyển hóa glucose và tổng hợp lipid, sự biểu hiện quá mức trong tế bào ung thư và tiềm năng làm đích thuốc kháng ung thư. Các thí nghiệm đã phát hiện rằng ức chế ACLY khiến tế bào ung thư ngừng phát triển và phân chia.

Thứ tư, chương tổng quan về flavonoid — cấu trúc khung C6-C3-C6, các phân nhóm (chalcone, flavonol, flavanone, flavone, isoflavone...) — và lý thuyết sàng lọc ảo (dựa trên phối tử và dựa trên cấu trúc mục tiêu), trong đó phương pháp docking được trình bày chi tiết.

Chương 2: Phương pháp nghiên cứu

Chương này mô tả đầy đủ quy trình thực nghiệm in silico. Phần chuẩn bị mục tiêu tác động trình bày quy trình xử lý protein bằng công cụ LigX của MOE 2008, gồm các bước: xác định ligand đồng kết tinh, proton hóa các amino acid, tối thiểu hóa năng lượng (phương pháp gradient liên hợp — Conj Grad, điểm dừng < 0,0001 kcal/mol, số bước lặp tối đa 10, điện tích Gasteiger-Huckel), xóa phân tử nước không liên kết, và lưu file protein.

Phần cơ sở dữ liệu hợp chất giới thiệu 30 phân tử: 10 dẫn xuất 2'-hydroxychalconoid (HC1–HC10, thay đổi nhóm thế trên vòng B), 10 dẫn xuất flavanone (HC11–HC20) và 10 dẫn xuất flavonol (HC21–HC30). Mỗi nhóm được trình bày kèm bảng cấu trúc và tên gọi IUPAC đầy đủ.

Quy trình chuẩn bị cấu trúc phân tử gồm ba bước: (1) xây dựng cấu trúc 2D bằng ChemBioDraw Ultra 12, chuyển sang SMILES rồi tạo cấu trúc 3D trong MOE; (2) tối thiểu hóa năng lượng; (3) lưu file và tối thiểu hóa lần hai trong Database Viewer. Thông số docking được cài đặt: vị trí docking là khoang gắn kết của ligand đồng kết tinh, phương pháp đặt là tam giác phân, giữ lại 10 cấu dạng tốt nhất mỗi phân tử.

Chương 3: Kết quả và bàn luận

Chương trình bày kết quả theo hai phần: re-docking và docking toàn bộ thư viện hợp chất.

Re-docking trên 4O2B: Ligand đồng kết tinh LOC (N-[(7S)-1,2,3,10-tetramethoxy-9-oxo-6,7-dihydro-5H-benzo[d]heptalen-7-yl]ethanamide) được tách ra và docking lại. Giá trị RMSD thu được thấp hơn 2,00 Å, xác nhận phương pháp và thông số MOE phù hợp.

Re-docking trên 6O0H: Ligand đồng kết tinh LBG (methyl 3-chloro-5-[(4,6-difluoro[1,1'-biphenyl]-3-yl)sulfamoyl]-4-hydroxybenzoate) được tách ra và docking lại; điểm docking (DS) đạt -14.1067, RMSD < 2,00 Å.

Kết quả docking toàn bộ được phân tích dựa trên điểm số docking (đơn vị kcal/mol, âm hơn là tốt hơn), biểu đồ tần suất tương tác amino acid tại khoang gắn kết và sơ đồ tương tác 2D của các hợp chất tiềm năng. Các loại tương tác được phân tích gồm: liên kết hydro, tương tác cation-π, tương tác π-π, tương tác ion và tương tác Van der Waals.


Kết quả và đóng góp

Kết quả chính

Tất cả 30 hợp chất đều được docking thành công vào khoang gắn kết của cả Tubulin (4O2B) và ACLY (6O0H). Nghiên cứu xác định 04 hợp chất tiềm năng dẫn đầu mỗi nhóm trên từng đích:

Hợp chất DS trên 4O2B (kcal/mol) DS trên 6O0H (kcal/mol)
HC9 -11,2206 -12,3810
HC19 -10,8180 -11,5916
HC30 (có trong văn bản, giá trị đầy đủ không hiển thị hoàn chỉnh trong fulltext) (có trong văn bản)

Kết quả re-docking xác nhận: RMSD < 2,00 Å trên cả 4O2B và 6O0H, khẳng định độ tin cậy của phương pháp và thông số MOE đã cài đặt. Điểm số docking âm của ligand đồng kết tinh LBG trên 6O0H đạt -14.1067 kcal/mol.

Giải pháp và kiến nghị

Tác giả đề xuất các hợp chất tiềm năng (trong đó có HC9, HC19, HC30) làm chất dẫn đầu để tiến hành các thử nghiệm thực nghiệm tiếp theo, bao gồm in vitro (thử nghiệm trong ống nghiệm) nhằm xác nhận hoạt tính kháng ung thư gan thực tế. Ngoài ra, tác giả đề xuất nghiên cứu xác định pharmacophore từ các cấu trúc tiềm năng để định hướng thiết kế và tổng hợp các hợp chất thế hệ tiếp theo.

Đóng góp so với nghiên cứu trước

Các dẫn xuất 2'-hydroxychalconoid, flavanone và flavonol sử dụng trong nghiên cứu được tổng hợp từ nhóm nghiên cứu của PGS.TS Hoàng Minh Hảo và chưa được đánh giá tương tác với Tubulin (4O2B) và ACLY (6O0H) thông qua docking phân tử trước đó. Nghiên cứu cung cấp dữ liệu in silico về các amino acid quan trọng tại khoang gắn kết của hai đích này khi tương tác với nhóm flavonoid — thông tin có giá trị cho việc tối ưu hóa cấu trúc trong các nghiên cứu tiếp theo.


Hạn chế và hướng nghiên cứu tiếp

Nghiên cứu được thực hiện hoàn toàn in silico bằng phương pháp docking phân tử, do đó kết quả mang tính dự đoán và cần được kiểm chứng bằng thực nghiệm in vitro và in vivo. Bản thân tác giả nêu rõ các hợp chất tiềm năng cần "tiếp tục nghiên cứu thực nghiệm kháng ung thư gan". Phạm vi nghiên cứu giới hạn ở 30 hợp chất thuộc ba nhóm flavonoid cụ thể và hai đích phân tử — chưa mở rộng sang các protein đích khác hoặc các dạng đột biến của Tubulin và ACLY vốn liên quan đến kháng thuốc. Phương pháp docking tĩnh (không tích hợp mô phỏng động lực học phân tử) cũng là giới hạn của phương pháp luận được sử dụng trong phạm vi đề tài này.


Giá trị tham khảo

Khóa luận phù hợp làm tài liệu tham khảo cho:

  • Sinh viên ngành Công nghệ Kỹ thuật Hóa học, Dược học, Hóa dược và Sinh học phân tử muốn tìm hiểu quy trình sàng lọc ảo in silico theo hướng tiếp cận dựa trên cấu trúc mục tiêu.
  • Người thực hiện đề tài docking phân tử bằng phần mềm MOE: phần Chương 2 cung cấp mô tả đầy đủ và có hệ thống về quy trình chuẩn bị protein (LigX), chuẩn bị phối tử (ChemBioDraw Ultra 12, tối thiểu hóa năng lượng) và cài đặt thông số docking.
  • Sinh viên nghiên cứu về flavonoid có thể tham khảo phần phân loại và cấu trúc các nhóm hợp chất (chalcone, flavanone, flavonol) cùng bảng danh mục cấu trúc IUPAC chi tiết của 30 dẫn xuất.
  • Người nghiên cứu về đích phân tử Tubulin và ACLY trong ung thư gan: phần Chương 1 và Chương 3 cung cấp thông tin về sinh học phân tử của hai đích này và dữ liệu tương tác amino acid tại khoang gắn kết.

Câu hỏi thường gặp

1. Tiêu chí nào được dùng để đánh giá độ tin cậy của kết quả docking trong nghiên cứu này? Nghiên cứu sử dụng phương pháp re-docking: tách ligand đồng kết tinh ra khỏi phức hợp tinh thể rồi docking lại vào khoang gắn kết. Kết quả được chấp nhận khi RMSD ≤ 2 Å — giá trị đo độ lệch giữa cấu dạng sau docking và cấu dạng có sẵn trong tinh thể. Cả 4O2B và 6O0H đều đạt tiêu chí này.

2. Tại sao nghiên cứu chọn Tubulin (4O2B) và ACLY (6O0H) làm đích phân tử, không phải các protein ung thư khác? Tubulin (β-tubulin) là thành phần không thể thiếu cho sự hình thành thoi phân bào trong quá trình nguyên phân của tế bào ung thư; ức chế Tubulin phá vỡ phân bào và kích hoạt apoptosis. ACLY là enzyme cuối cùng của con đường đường phân đồng thời là điểm khởi đầu cho tổng hợp lipid — hai quá trình trao đổi chất đặc trưng bị thay đổi trong tế bào ung thư. Cả hai đều là mục tiêu đã được các nghiên cứu trước xác nhận là có tiềm năng trong điều trị ung thư, trong đó có ung thư gan.

3. Ligand đồng kết tinh của 4O2B và 6O0H là gì, và chúng có vai trò gì trong nghiên cứu? Ligand đồng kết tinh của 4O2B là LOC (N-[(7S)-1,2,3,10-tetramethoxy-9-oxo-6,7-dihydro-5H-benzo[d]heptalen-7-yl]ethanamide) — dẫn xuất của colchicine. Ligand đồng kết tinh của 6O0H là LBG (methyl 3-chloro-5-[(4,6-difluoro[1,1'-biphenyl]-3-yl)sulfamoyl]-4-hydroxybenzoate). Cả hai được sử dụng để: (1) xác định vị trí khoang gắn kết cho docking, và (2) thực hiện re-docking để kiểm định phương pháp.

4. Trong số 30 hợp chất, những chất nào cho điểm số docking tốt nhất trên cả hai đích? HC9 và HC19 nổi bật trên cả hai đích: HC9 đạt -11,2206 kcal/mol trên Tubulin (4O2B) và -12,3810 kcal/mol trên ACLY (6O0H); HC19 đạt -10,8180 kcal/mol trên 4O2B và -11,5916 kcal/mol trên 6O0H. Đây đều là dẫn xuất của 2'-hydroxychalconoid (thuộc nhóm HC1–HC10).

5. Các tương tác phân tử nào được phân tích trong kết quả docking? Phần mềm MOE 2008 được dùng để phân tích và hiển thị các loại tương tác giữa hợp chất và amino acid tại khoang gắn kết, bao gồm: liên kết hydro, tương tác cation-π, tương tác π-π, tương tác ion và tương tác bề mặt Van der Waals (qua sự tiếp xúc bề mặt thân nước/kỵ nước). Kết quả được trình bày dưới dạng sơ đồ tương tác 2D và biểu đồ tần suất tương tác của từng amino acid tại khoang gắn kết của 4O2B và 6O0H.


Kết luận

Khóa luận thực hiện sàng lọc ảo in silico ba nhóm dẫn xuất flavonoid (2'-hydroxychalconoid, flavanone và flavonol) trên hai đích phân tử kháng ung thư là Tubulin (4O2B) và ACLY (6O0H) bằng phần mềm MOE 2008, với quy trình được kiểm định độ tin cậy qua re-docking (RMSD < 2 Å trên cả hai đích). Từ 30 hợp chất khảo sát, nghiên cứu xác định được các chất dẫn đầu tiềm năng — trong đó HC9 và HC19 cho điểm số docking âm sâu trên cả hai đích — và cung cấp dữ liệu về các amino acid chủ chốt tại khoang gắn kết có thể được khai thác trong thiết kế pharmacophore. Các kết quả in silico này tạo cơ sở định hướng cho các thử nghiệm in vitro và in vivo tiếp theo nhằm xác nhận hoạt tính kháng ung thư gan thực tế của các hợp chất tiềm năng.