Tổng quan nghiên cứu

Năm 2020, Cơ quan Nghiên cứu Ung thư Quốc tế (IARC) ghi nhận khoảng 19,3 triệu ca mắc mới và gần 10 triệu ca tử vong do ung thư trên toàn cầu. Riêng tại Việt Nam, số ca mắc mới trong cùng năm ước tính lên đến 182.563 ca, với 122.690 ca tử vong — trong đó chỉ khoảng 20–30% bệnh nhân được phát hiện ở giai đoạn đầu và tỷ lệ tử vong trong vòng 1 năm chiếm tới 25%. Thực trạng đáng báo động này thúc đẩy nhu cầu cấp thiết phải tìm kiếm các hướng điều trị mới, ít tác dụng phụ hơn so với hóa trị và xạ trị truyền thống.

Thụ thể yếu tố tăng trưởng biểu bì (EGFR) được xác định là một trong những đích phân tử quan trọng nhất trong điều trị ung thư, liên quan trực tiếp đến ung thư phổi không tế bào nhỏ (NSCLC), ung thư gan, ung thư đại trực tràng và ung thư nội mạc tử cung. Ức chế miền tyrosine kinase của EGFR (EGFR-TK) bằng các phân tử nhỏ đã được chứng minh là chiến lược điều trị nhắm đích hiệu quả, song các thuốc tổng hợp hiện hành như gefitinib, erlotinib và lapatinib vẫn gây nhiều tác dụng phụ không mong muốn.

Luận văn thạc sĩ này, thực hiện tại Phòng Thí nghiệm Hóa học Tính toán và Mô phỏng (LCCM), Trường Đại học Quy Nhơn, tập trung nghiên cứu tương tác của các hợp chất thiên nhiên có trong cây Xạ đen (Celastrus hindsii Benth) và cây An xoa (Helicteres hirsuta Loureiro) với miền EGFR-TK (PDB ID: 3POZ) bằng phương pháp in silico. Mục tiêu kép của công trình là: (1) mô phỏng và phân tích định lượng khả năng tương tác phân tử; và (2) đề xuất các hợp chất tiềm năng phục vụ các nghiên cứu thực nghiệm và lâm sàng tiếp theo trong điều trị ung thư liên quan đến EGFR.


Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu xây dựng trên hai trụ cột lý thuyết chính trong hóa học tính toán hiện đại.

Thuyết phiếm hàm mật độ (Density Functional Theory – DFT) là nền tảng tính toán cấu trúc điện tử. Xuất phát từ hai định lý Hohenberg–Kohn (1964), DFT xác lập nguyên lý rằng mật độ electron ρ(r) đủ để xác định toàn bộ tính chất trạng thái cơ bản của hệ. Phương pháp Kohn–Sham cho phép giải bài toán N electron thông qua hệ orbital một electron, trong đó năng lượng trao đổi–tương quan E_xc[ρ] là đại lượng trung tâm cần xấp xỉ. Nghiên cứu áp dụng phiếm hàm hỗn hợp B3LYP với bộ hàm cơ sở mở rộng 6–31++G(d,p) — lựa chọn cân bằng tốt giữa độ chính xác và chi phí tính toán, đặc biệt phù hợp với các phân tử polyphenolic nhiều vòng thơm như flavonoid và lignan. B3LYP kết hợp ba thông số Becke (a = 0,2; b = 0,7; c = 0,8) với phiếm hàm tương quan LYP, cho phép mô tả tốt các tương tác liên kết yếu thường gặp trong phức protein–phối tử.

Phương pháp docking phân tử là công cụ trung tâm để đánh giá ái lực liên kết. Năm khái niệm chuyên ngành cốt lõi gồm: (1) năng lượng liên kết (binding energy – BE) đo lường mức độ bền của phức; (2) hàm đánh giá (scoring functions – SF) dựa vào trường lực AMBER để lượng hóa tương tác Van der Waals, tĩnh điện và liên kết hydrogen; (3) thuật toán di truyền Lamarckian (LGA) để tối ưu cấu dạng phối tử trong không gian thụ thể; (4) điện tích Gasteiger và Kollman cho protein và phối tử; và (5) lưới AutoGrid xác định vùng không gian khảo sát tương tác.

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

Nghiên cứu sử dụng thiết kế tính toán thuần túy (purely computational), không yêu cầu cỡ mẫu sinh học theo nghĩa thông thường. Thay vào đó, đối tượng phân tích là toàn bộ các hợp chất thiên nhiên đã được phân lập và xác định cấu trúc từ cây Xạ đen và cây An xoa, được phân thành 6 nhóm: flavonoid, lignan, steroid, terpenoid, glycoside và các hợp chất khác — bao gồm hàng chục phân tử được mã hóa theo quy ước F-, L-, S-, T-, G-, D- kết hợp chỉ số nguồn gốc (c: Xạ đen; h: An xoa).

Quy trình phân tích gồm 5 bước tuần tự: (1) tải cấu trúc protein 3POZ từ Protein Data Bank; (2) tối ưu hóa hình học phối tử ở mức B3LYP/6–31++G(d,p) bằng phần mềm Gaussian 09; (3) chuẩn bị file .pdbqt cho cả protein và phối tử qua AutoDockTools 1.5; (4) thực thi docking bằng AutoDock 4 với thuật toán LGA; và (5) trực quan hóa và phân tích tương tác bằng Discovery Studio Visualizer và PyMol.

Để đảm bảo tính tin cậy, nghiên cứu tiến hành thử nghiệm hai điều kiện điện tích (Gasteiger và Kollman) và nhiều tham số hộp lưới AutoGrid khác nhau trước khi cố định điều kiện tối ưu cho toàn bộ phân tích. Phương pháp docking bán linh động được chọn vì cho phép phối tử chuyển động tự do trong không gian thụ thể trong khi protein được xem là cứng nhắc — phù hợp với giai đoạn sàng lọc sơ bộ lead compound.


Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phát hiện 1: Xác nhận điều kiện docking tối ưu. Kết quả so sánh hai loại điện tích (Gasteiger và Kollman) cho thấy điện tích Kollman mang lại năng lượng liên kết âm hơn và tái tạo vị trí gắn kết chuẩn tốt hơn đối với các chất đối chiếu (gefitinib, erlotinib, lapatinib). Khung lưới AutoGrid với kích thước tối ưu bao phủ trọn vùng hoạt tính tyrosine kinase của 3POZ được cố định cho tất cả các thí nghiệm docking tiếp theo, đảm bảo tính so sánh của kết quả.

Phát hiện 2: Nhóm flavonoid thể hiện ái lực liên kết vượt trội. Trong số các nhóm hợp chất được khảo sát, các flavonoid — đặc biệt là các hợp chất có mặt trong cả hai loài cây (ký hiệu F-xch) — cho giá trị năng lượng liên kết âm nhất, tương đương hoặc vượt qua mức tham chiếu của các thuốc FDA-approved. Điều này phù hợp với nghiên cứu thực nghiệm của giáo sư Lê Thế Trung (từ năm 1987) xác định flavonoid là nhóm hoạt chất ức chế ung thư chính trong cây Xạ đen, cũng như với đặc tính đã được chứng minh của flavonoid trong việc ức chế protein kinase. Kết quả docking được minh họa rõ ràng qua các hình 3.3 và 3.4, cho thấy mạng lưới tương tác hydrogen và kỵ nước phong phú tại vùng hoạt tính.

Phát hiện 3: Các triterpene khung lupan trong nhóm terpenoid có triển vọng đáng kể. Các hợp chất terpenoid, đặc biệt các triterpene khung lupan (bao gồm lup-20(29)-ene-3β,11β-diol, lup-20(29)-ene-3-one, lup-5,20(29)-diene-3-one đã được phân lập bởi Trần Văn Sung và cộng sự năm 2007), cho thấy năng lượng liên kết với EGFR-TK đứng thứ hai trong các nhóm khảo sát. Biểu đồ cột (hình 3.8) trực quan hóa sự phân bố BE của toàn bộ nhóm terpenoid, làm nổi bật rõ ràng các hợp chất vượt trội so với phần còn lại của nhóm.

Phát hiện 4: Nhóm steroid bị giới hạn bởi khung cấu trúc cứng nhắc. Các hợp chất steroid, dù có khung cyclopenta[a]phenantren cố định (hình 3.7), cho năng lượng liên kết kém ổn định hơn. Nguyên nhân được quy cho sự thiếu linh hoạt cấu dạng — bộ khung phẳng, cứng của steroid hạn chế khả năng hình thành các tiếp xúc bổ sung với các vị trí không gian hẹp trong túi liên kết của EGFR-TK. Tương tác chủ đạo được ghi nhận là loại kỵ nước, với rất ít liên kết hydrogen so với nhóm flavonoid.

Thảo luận kết quả

Sự vượt trội của flavonoid trong tương tác với EGFR-TK có thể giải thích qua ba cơ chế: (1) cấu trúc vòng thơm đa nhân tạo điều kiện thuận lợi cho các tương tác π–π với phenylalanine, tyrosine trong vùng hoạt tính; (2) nhóm hydroxyl phong phú trên khung flavone cung cấp nhiều điểm hình thành liên kết hydrogen với các amino acid quan trọng như Thr790, Met793 trong EGFR-TK; và (3) khối lượng phân tử vừa phải cho phép phối tử xâm nhập sâu vào túi ATP-binding. So sánh với các thuốc tổng hợp tham chiếu (khối lượng phân tử từ 393,1 Da), một số hợp chất flavonoid và terpenoid trong nghiên cứu này thỏa mãn tương tự quy tắc Lipinski, gợi mở tiềm năng phát triển thành thuốc uống.

Kết quả cũng tương đồng với báo cáo của Nguyễn Huy Cường (luận án tiến sĩ, 2008), trong đó một số triterpenoid cây Xạ đen thể hiện hoạt tính gây độc tế bào tốt trên dòng Hep-G2 và Lu-1 — khẳng định sự nhất quán giữa kết quả in silico và in vitro, tăng cường độ tin cậy của phương pháp docking áp dụng.


Đề xuất và khuyến nghị

Từ kết quả docking phân tử, luận văn đưa ra bốn nhóm giải pháp và khuyến nghị cụ thể:

1. Ưu tiên nghiên cứu thực nghiệm cho các hợp chất đầu bảng. Các nhà nghiên cứu thực nghiệm nên tập trung trước tiên vào các flavonoid (nhóm F-xch) và triterpene khung lupan có giá trị BE âm nhất trong nghiên cứu này. Các thí nghiệm MTT assay, đánh giá IC₅₀ trên các dòng tế bào ung thư NSCLC (A549), ung thư gan (Hep-G2) và ung thư đại trực tràng (HCT116) nên được ưu tiên triển khai trong giai đoạn 2025–2027. Mục tiêu đạt giá trị IC₅₀ ở mức dưới 10 μM để đảm bảo ý nghĩa dược lý.

2. Áp dụng phương pháp động lực học phân tử (MD) để xác nhận độ bền phức. Các phức protein–phối tử tiềm năng được xác định qua docking cần được kiểm chứng thêm bằng mô phỏng MD trong vòng ít nhất 100 ns, nhằm đánh giá tính ổn định nhiệt động lực học của liên kết theo thời gian. Phòng LCCM và các đơn vị hợp tác tại Việt Nam có thể triển khai bước này trong năm 2025–2026 với phần mềm GROMACS hoặc AMBER.

3. Mở rộng khảo sát sang các đột biến EGFR kháng thuốc. Đột biến T790M của EGFR là cơ chế kháng thuốc phổ biến nhất với gefitinib và erlotinib. Nghiên cứu tiếp theo nên tiến hành docking các hợp chất tiềm năng lên cấu trúc EGFR mang đột biến T790M và L858R để đánh giá khả năng vượt qua cơ chế kháng thuốc — một hướng mà các thuốc thiên nhiên thường có lợi thế nhờ cơ chế tác động đa đích.

4. Xây dựng cơ sở dữ liệu hợp chất dược liệu Việt Nam phục vụ sàng lọc ảo. Việt Nam có gần 4.000 loài thực vật dùng làm thuốc trong y học cổ truyền, song hầu hết chưa được nghiên cứu ở cấp độ phân tử. Bộ Giáo dục và Đào tạo cùng Bộ Y tế nên phối hợp đầu tư xây dựng cơ sở dữ liệu hợp chất thiên nhiên Việt Nam chuẩn hóa theo định dạng SMILES/SDF, phục vụ sàng lọc ảo quy mô lớn (virtual screening) với chi phí thấp hơn thực nghiệm ít nhất 10 lần trong giai đoạn sơ bộ.


Đối tượng nên tham khảo luận văn

Nhà nghiên cứu hóa học dược liệu và hóa học tính toán sẽ tìm thấy trong luận văn một quy trình docking phân tử chuẩn hóa, từ tối ưu hóa cấu trúc B3LYP/6–31++G(d,p) đến phân tích tương tác phối tử–protein bằng Discovery Studio. Bộ dữ liệu năng lượng liên kết của hàng chục hợp chất thuộc 6 nhóm chất từ hai loài cây Xạ đen và An xoa là tài nguyên tham chiếu có giá trị cho các nghiên cứu sàng lọc dược liệu tiếp theo.

Sinh viên và học viên cao học ngành Hóa lý thuyết, Hóa sinh và Dược học có thể sử dụng luận văn như một hướng dẫn thực hành in silico toàn diện: từ lý thuyết DFT, nguyên lý docking phân tử, cách chọn lựa và so sánh hàm đánh giá (force field-based, empirical, knowledge-based), đến kỹ năng sử dụng Gaussian 09, AutoDock 4 và PyMol trong thực tế.

Dược sĩ, bác sĩ và chuyên gia y học cổ truyền quan tâm đến cơ sở khoa học phân tử của các bài thuốc dân gian sẽ tìm thấy bằng chứng định lượng lý giải tại sao cây Xạ đen và cây An xoa có tiềm năng kháng ung thư, đặc biệt qua dữ liệu tương tác của các flavonoid và terpenoid với thụ thể EGFR — làm nền tảng cho tư vấn sử dụng dược liệu có căn cứ khoa học.

Nhà hoạch định chính sách và quản lý nghiên cứu khoa học trong lĩnh vực y tế và dược phẩm sẽ thấy đây là minh chứng cụ thể cho tiềm năng khai thác nguồn tài nguyên thực vật phong phú của Việt Nam (gần 14.000 loài thực vật bậc cao, khoảng 4.000 loài dùng làm thuốc) thông qua công nghệ tính toán chi phí thấp, rút ngắn đáng kể giai đoạn phát hiện thuốc từ nhiều thập kỷ xuống còn vài năm.


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

1. Phương pháp in silico và docking phân tử là gì, và tại sao lại áp dụng trong nghiên cứu này?

Phương pháp in silico là các phương pháp mô phỏng và tính toán lý thuyết thực hiện trên máy tính. Docking phân tử cụ thể là kỹ thuật dự đoán cách một phân tử nhỏ (phối tử) gắn vào protein đích và ước tính năng lượng liên kết. Ứng dụng trong nghiên cứu này giúp sàng lọc hàng chục hợp chất thiên nhiên một cách nhanh chóng và kinh tế, giảm thiểu thực nghiệm thử–sai tốn kém, trước khi chuyển sang kiểm chứng in vitro cho các ứng viên hàng đầu.

2. Cây Xạ đen và cây An xoa chứa những nhóm hợp chất nào có tiềm năng kháng ung thư?

Cả hai loài cây đều giàu các hợp chất polyphenolic tự nhiên gồm: flavonoid, lignan, steroid, terpenoid, glycoside và một số hợp chất đặc thù khác. Trong đó, flavonoid và triterpene khung lupan được xác định là hai nhóm có ái lực liên kết cao nhất với EGFR-TK qua kết quả docking, phù hợp với các nghiên cứu thực nghiệm in vitro trên các dòng tế bào ung thư gan Hep-G2 và ung thư phổi Lu-1 đã được công bố trước đó.

3. Tại sao thụ thể EGFR là đích phân tử quan trọng trong điều trị ung thư?

EGFR điều tiết các con đường tín hiệu thiết yếu cho tăng trưởng, phân chia và biệt hóa tế bào. Khi bị kích hoạt liên tục do đột biến hoặc biểu hiện quá mức, EGFR thúc đẩy tế bào phân chia mất kiểm soát, gây ra ung thư phổi (NSCLC), ung thư gan, ung thư đại trực tràng và nhiều loại khác. Ức chế miền tyrosine kinase của EGFR cắt đứt chuỗi tín hiệu này, làm chậm hoặc dừng sự phát triển khối u mà ít ảnh hưởng đến tế bào lành hơn hóa trị thông thường.

4. Kết quả docking này có thể ứng dụng ngay vào điều trị lâm sàng không?

Chưa thể ứng dụng trực tiếp. Docking phân tử là bước sàng lọc sơ bộ in silico, cung cấp bằng chứng lý thuyết về tiềm năng tương tác. Để tiến tới ứng dụng lâm sàng, các hợp chất tiềm năng cần trải qua ít nhất 3 giai đoạn tiếp theo: đánh giá hoạt tính in vitro trên dòng tế bào ung thư, thử nghiệm in vivo trên mô hình động vật, và cuối cùng là thử nghiệm lâm sàng trên người — một lộ trình có thể kéo dài 10–15 năm theo tiêu chuẩn FDA.

5. Sự khác biệt giữa gefitinib/erlotinib và các hợp chất thiên nhiên trong nghiên cứu này là gì?

Gefitinib và erlotinib là các anilinoquinazoline tổng hợp (khối lượng từ 393,1 Da), được thiết kế đặc hiệu gắn vào vùng hoạt tính tyrosine kinase của EGFR theo cơ chế cạnh tranh với ATP. Tuy hiệu quả ban đầu cao, chúng gây kháng thuốc sau một thời gian và có tác dụng phụ đáng kể. Các hợp chất thiên nhiên (flavonoid, terpenoid) có cơ chế tác động đa đích tiềm năng, có thể vượt qua một số cơ chế kháng thuốc, đồng thời ít độc hơn do cơ thể người đã có quá trình thích nghi tiến hóa lâu dài với các hợp chất này.


Kết luận

  • Đóng góp phương pháp: Thiết lập quy trình docking phân tử chuẩn hóa với AutoDock 4, kết hợp tối ưu hóa cấu trúc B3LYP/6–31++G(d,p) và so sánh điều kiện điện tích/lưới, áp dụng cho hệ hợp chất thiên nhiên Việt Nam.
  • Đóng góp khoa học: Lần đầu tiên làm rõ bản chất tương tác phân tử của các hợp chất cây Xạ đen và cây An xoa với EGFR-TK (PDB: 3POZ) ở cấp độ nguyên tử.
  • Phát hiện nổi bật: Flavonoid và triterpene khung lupan là hai nhóm hợp chất có ái lực liên kết cao nhất với EGFR-TK, với năng lượng liên kết cạnh tranh so với thuốc tổng hợp tham chiếu.
  • Ý nghĩa thực tiễn: Cung cấp danh sách ứng viên hợp chất tiềm năng có cơ sở dữ liệu in silico, rút ngắn và định hướng giai đoạn nghiên cứu thực nghiệm tiếp theo.
  • Hướng phát triển tiếp theo: Ưu tiên kiểm chứng in vitro cho các flavonoid đầu bảng (2025–2026), tiếp theo là mô phỏng MD và đánh giá trên mô hình đột biến EGFR T790M (2026–2027).

Nghiên cứu góp phần khẳng định tiềm năng to lớn của kho tàng dược liệu Việt Nam trong phát triển thuốc kháng ung thư thế hệ mới — hướng tới liệu pháp nhắm đích EGFR từ nguồn gốc thiên nhiên, an toàn hơn và có khả năng vượt qua cơ chế kháng thuốc hiện tại. Các nhà nghiên cứu và tổ chức quan tâm có thể liên hệ với Phòng Thí nghiệm LCCM, Trường Đại học Quy Nhơn để hợp tác phát triển nghiên cứu tiếp theo.