Tổng quan nghiên cứu

Bệnh Alzheimer hiện đang là một trong những thách thức y tế và kinh tế xã hội nghiêm trọng nhất trên toàn cầu. Theo báo cáo y tế năm 2024, ước tính có khoảng 7,2 triệu người từ 65 tuổi trở lên tại Mỹ đang sống chung với căn bệnh này (chiếm khoảng 11% dân số cao tuổi) và dự kiến sẽ tăng lên 13,8 triệu người vào năm 2060. Trên phạm vi toàn cầu, số lượng người mắc Alzheimer ghi nhận khoảng 50 triệu người vào năm 2020 và dự báo vượt ngưỡng 152 triệu người vào năm 2050 do tốc độ già hóa dân số nhanh chóng.

Trong bối cảnh các liệu pháp miễn dịch loại bỏ mảng bám amyloid mới chỉ làm chậm tiến trình suy giảm nhận thức từ 27% đến 35%, việc can thiệp vào enzyme Acetylcholinesterase (AChE) nhằm duy trì nồng độ acetylcholin tại khe synap thần kinh vẫn là chiến lược điều trị triệu chứng then chốt. Tuy nhiên, quy trình nghiên cứu và phát triển một hoạt chất ức chế AChE mới theo phương pháp thực nghiệm truyền thống thường đòi hỏi kinh phí lên tới hàng trăm triệu USD và kéo dài từ 10 đến 15 năm.

Đề tài luận văn thạc sĩ dược học chuyên ngành Hóa dược tại Trường Đại học Dược Hà Nội (năm 2025) được triển khai nhằm giải quyết bài toán tối ưu hóa quy trình sàng lọc thuốc. Nghiên cứu tập trung vào 3 mục tiêu cụ thể:

  1. Xây dựng và đánh giá các mô hình học máy phân loại hoạt tính ức chế AChE từ tập dữ liệu cấu trúc hóa học.
  2. Ứng dụng mô hình học máy kết hợp sàng lọc ảo và đánh giá dược động học để tìm kiếm các hợp chất tiềm năng từ cơ sở dữ liệu tự nhiên.
  3. Thiết lập quy trình mô phỏng tương tác và động lực học phân tử để khảo sát độ bền phức hợp giữa phối tử và enzyme mục tiêu.

Nghiên cứu mang ý nghĩa thực tiễn lớn khi rút ngắn hơn 50% thời gian sàng lọc ban đầu, định hướng chính xác các khung cấu trúc hóa học tiềm năng trước khi tiến hành tổng hợp hóa dược và thử nghiệm sinh học in vitro.

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 được xây dựng dựa trên các nền tảng lý thuyết sinh học phân tử, hóa tin học và tính toán lượng tử:

  • Cơ chế bệnh sinh Alzheimer và giả thuyết Cholinergic: Sự thoái hóa của các nơ-ron cholinergic dẫn đến thiếu hụt chất dẫn truyền thần kinh acetylcholin. Enzyme AChE xúc tác quá trình thủy phân acetylcholin với tốc độ cực đại lên tới 25.000 phân tử mỗi giây.
  • Cấu trúc sinh học 3D của Acetylcholinesterase: Enzyme có dạng hình elip với kích thước khoảng 45Å x 60Å x 65Å, chứa một hốc sâu khoảng 20Å dẫn tới trung tâm hoạt động. Cấu trúc gồm bộ ba xúc tác Ser200 - His440 - Glu327, vị trí anion với gốc Trp84 liên kết cation-π, và vị trí anion ngoại vi (PAS) ở lối vào hốc gồm Trp279, Tyr70, Asp72, Tyr121, Tyr334.
  • Lý thuyết học máy có giám sát (Supervised Machine Learning): Ứng dụng các thuật toán phân loại gồm Máy vector hỗ trợ (SVM), Rừng ngẫu nhiên (Random Forest) và Tăng cường độ dốc cực đại (XGBoost) dựa trên các bộ mô tả phân tử 2D/3D (fingerprints và molecular descriptors).
  • Mô phỏng tương tác (Molecular Docking) và Động lực học phân tử (Molecular Dynamics - MD): Dựa trên các trường lực cơ học phân tử để đánh giá năng lượng tự do liên kết Gibbs (ΔG), sự biến đổi không gian - thời gian của hệ thống và độ ổn định liên kết qua phương trình chuyển động Newton.

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

Quy trình nghiên cứu kết hợp đa tầng tính toán in silico được thực hiện qua các bước chuẩn hóa:

  • Thu thập và tiền xử lý dữ liệu: Cỡ mẫu gồm 1.055 phân tử có dữ liệu hoạt tính ức chế AChE in vitro được trích xuất từ cơ sở dữ liệu ChEMBL. Dữ liệu được làm sạch bằng cách loại bỏ các muối vô cơ, ion kim loại nặng, hợp chất trùng lặp và các mảnh phân tử lỗi.
  • Gán nhãn hoạt tính: Ngưỡng phân loại nhị phân được xác lập dựa trên giá trị IC50 của hoạt chất đối chứng Galantamin (2055 nM). Phân tử có IC50 ≤ 2055 nM được gán nhãn hoạt tính (active), chiếm 674 hợp chất (63,88%); phân tử có IC50 > 2055 nM được gán nhãn không hoạt tính (inactive), chiếm 381 hợp chất (36,12%).
  • Phân chia tập dữ liệu và huấn luyện mô hình: Bộ dữ liệu 1.055 chất được chia ngẫu nhiên theo tỷ lệ 80% cho tập huấn luyện (844 phân tử) và 20% cho tập kiểm tra độc lập (211 phân tử) với seed cố định (random_state = 42). Các mô hình SVM, Random Forest và XGBoost được tối ưu hóa siêu tham số bằng Hyperopt và kiểm định chéo 5 lần (5-fold cross-validation).
  • Đánh giá và sàng lọc ảo: Hiệu năng mô hình được kiểm định qua các chỉ số Accuracy, Precision, Recall, ROC-AUC và PRC-AUC. Các mô hình tối ưu được áp dụng để sàng lọc 100 hợp chất từ cơ sở dữ liệu COCONUT (50 chất) và các dược liệu y học cổ truyền (50 chất).
  • Mô phỏng Docking và MD 100 ns: Thực hiện docking phân tử trên phần mềm AutoDock Vina, dự đoán tính chất ADME/BBB qua SwissADME, pkCSM và PreADMET. Cuối cùng, thực hiện mô phỏng động lực học phân tử 100 ns để phân tích chi tiết các thông số RMSD, RMSF và năng lượng liên kết Gibbs trong 20 ns cuối.

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

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

Quá trình phân tích dữ liệu và tính toán mô phỏng mang lại các kết quả nổi bật:

  • Hiệu năng vượt trội của thuật toán XGBoost và Random Forest: Sau quá trình tinh chỉnh siêu tham số, mô hình XGBoost và Random Forest thể hiện khả năng phân loại vượt trội trên tập kiểm tra 211 phân tử. Chỉ số diện tích dưới đường cong ROC-AUC đạt trên 0,87 và PRC-AUC đạt mức tiệm cận 0,90, cao hơn khoảng 5,8% so với mô hình Support Vector Machine (SVM) cơ sở.
  • Sàng lọc thành công 8 hợp chất tự nhiên có ái lực cao: Từ 100 cấu trúc tự nhiên ban đầu, quá trình sàng lọc học máy và docking phân tử đã xác định được 8 hợp chất tiềm năng với mức năng lượng liên kết dao động từ -10,2 kcal/mol đến -12,4 kcal/mol, vượt trội hơn so với chất đối chứng Galantamin (-8,8 kcal/mol đến -9,8 kcal/mol).
  • Đặc tính dược động học và độ thấm qua hàng rào máu não (BBB): Các hợp chất tiêu biểu như dẫn xuất alkaloid Huperzin A, Hookerianin, Galanthin và một số cấu trúc dị vòng pyrimidine đạt hệ số logPS > -3,0 (thấm tốt vào hệ thần kinh trung ương), giá trị logBB dương từ 0,11 đến 0,72 và tỷ lệ liên kết protein huyết tương (Fu) phù hợp cho việc phân bố tại mô não.
  • Độ ổn định động học phân tử trong 100 ns mô phỏng: Phức hợp giữa các hợp chất được lựa chọn và enzyme AChE duy trì cấu hình ổn định với độ lệch bình phương trung bình (RMSD) dao động dưới 2,0 Å trong suốt 20 ns cuối của quá trình mô phỏng, không xuất hiện hiện tượng tách rời phối tử khỏi túi liên kết.

Thảo luận kết quả

Hiệu quả ức chế mạnh mẽ của các hợp chất tiềm năng bắt nguồn từ cơ chế gắn kết đa điểm (dual-site binding). Thay vì chỉ tương tác đơn thuần tại trung tâm xúc tác đáy hốc (Ser200, His440, Trp84), các phối tử được chọn còn hình thành liên kết π-π stacking và liên kết hydro vững chắc với các acid amin thuộc vị trí anion ngoại vi PAS (đặc biệt là Trp279 và Tyr70). Tương tác kép này không chỉ ngăn cản cơ chất acetylcholin tiếp cận bộ ba xúc tác mà còn ức chế quá trình kết tập mảng amyloid qua trung gian PAS của enzyme AChE.

Khi so sánh với các nghiên cứu quốc tế công bố trong giai đoạn 2023-2024 (với độ chính xác mô hình học máy đạt từ 85% đến 91%), kết quả của luận văn thể hiện sự tương đồng cao về độ tin cậy. Việc phân tích dữ liệu phân bố hoạt tính thông qua biểu đồ histogram log10(IC50) và biểu đồ phân tán năng lượng tự do Gibbs giúp chứng minh rằng quá trình tiền xử lý 1.055 mẫu dữ liệu đã triệt tiêu hiệu quả các yếu tố nhiễu thực nghiệm, mang lại giá trị dự đoán chuẩn xác.

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

Dựa trên các kết quả tính toán in silico đã đạt được, luận văn đưa ra 4 khuyến nghị và định hướng hành động cụ thể:

  1. Tiến hành tổng hợp hóa học và thử nghiệm sinh học in vitro: Ưu tiên thiết kế quy trình tổng hợp và thử nghiệm hoạt tính ức chế AChE trên mô hình enzyme tái tổ hợp đối với 3 dẫn xuất tự nhiên tiềm năng nhất (gồm khung cấu trúc pyrimidine chọn lọc và dẫn xuất alkaloid tương tự Huperzin A). Mục tiêu hướng đến việc đạt giá trị IC50 thực nghiệm < 100 nM, do nhóm nghiên cứu Bộ môn Hóa dược thực hiện trong lộ trình 6 đến 12 tháng.
  2. Mở rộng quy mô cơ sở dữ liệu học máy: Nâng số lượng mẫu huấn luyện từ 1.055 hợp chất lên hơn 5.000 hợp chất thông qua việc tích hợp thêm dữ liệu từ BindingDB và PubChem. Đồng thời, bổ sung các bộ mô tả không gian 3D nâng cao (như WHIM, GETAWAY, 3D-MoRSE) nhằm nâng cao độ chính xác dự đoán tổng thể lên trên 92% trước quý 4 năm 2026.
  3. Kéo dài thời gian mô phỏng động lực học phân tử: Tận dụng hệ thống máy chủ tính toán hiệu năng cao (HPC) để gia tăng thời gian mô phỏng MD từ 100 ns lên 500 ns hoặc 1.000 ns (1 µs) đối với các phức hợp hàng đầu. Điều này giúp tính toán chính xác năng lượng liên kết tự do bằng phương pháp MM-PBSA/MM-GBSA và khảo sát các thay đổi cấu hình allosteric dài hạn.
  4. Đánh giá độc tính tế bào và kiểm tra dược động học tiền lâm sàng: Phối hợp cùng các đơn vị dược lý thực hiện các thử nghiệm độc tính tế bào thần kinh (neurotoxicity) và kiểm tra khả năng vượt hàng rào máu não trên mô hình động vật gặm nhấm, hoàn thành trong giai đoạn 18 đến 24 tháng tới.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng chính:

  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Dược học, Hóa dược: Nguồn tư liệu tham khảo giá trị về phương pháp luận kết hợp giữa Hóa tính toán (Computational Chemistry), Hóa tin học (Cheminformatics) và mô hình hóa sinh học phân tử.
  • Kỹ sư tin sinh học và chuyên gia khoa học dữ liệu y sinh: Tham khảo chi tiết quy trình xử lý dữ liệu cấu trúc hóa học, tối ưu hóa thuật toán XGBoost/Random Forest và phương pháp đánh giá đường cong Precision-Recall trên tập dữ liệu mất cân bằng.
  • Các viện nghiên cứu và trung tâm R&D tại các doanh nghiệp dược phẩm: Ứng dụng quy trình sàng lọc ảo tích hợp nhằm rút ngắn thời gian và cắt giảm từ 40% đến 60% chi phí thử nghiệm thực nghiệm ban đầu trong các dự án phát triển thuốc hướng thần kinh.
  • Dược sĩ lâm sàng và bác sĩ chuyên khoa thần kinh: Cung cấp cái nhìn chuyên sâu ở cấp độ nguyên tử về cơ chế tác động của các thuốc ức chế AChE, giúp hiểu rõ động học tương tác của các hoạt chất điều trị Alzheimer thế hệ mới.

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

1. Tại sao nghiên cứu lại chọn mốc IC50 = 2055 nM để phân loại hoạt tính?
Ngưỡng 2055 nM được thiết lập dựa trên giá trị IC50 thực nghiệm của Galantamin – một thuốc ức chế AChE đã được FDA phê duyệt và sử dụng phổ biến trên lâm sàng. Việc lấy mốc này làm ranh giới giúp mô hình học máy tập trung nhận diện những cấu trúc có hiệu lực sinh học tương đương hoặc vượt trội hơn các loại thuốc hiện hành.

2. Mô hình học máy nào cho hiệu quả sàng lọc chất ức chế AChE tốt nhất trong luận văn?
Thuật toán XGBoost kết hợp kỹ thuật tinh chỉnh siêu tham số và đánh giá chéo 5 lần mang lại hiệu năng cao nhất. Mô hình đạt chỉ số ROC-AUC trên 0,87 và PRC-AUC tiệm cận 0,90 trên tập kiểm tra độc lập gồm 211 phân tử, thể hiện khả năng phân biệt chính xác giữa nhóm hoạt tính và không hoạt tính.

3. Ý nghĩa của việc thực hiện mô phỏng động lực học phân tử (MD) 100 ns là gì?
Molecular Docking chỉ cung cấp cấu hình tĩnh tại một thời điểm, trong khi mô phỏng MD 100 ns cho phép quan sát sự chuyển động của toàn bộ hệ thống protein - phối tử trong môi trường nước và ion sinh lý. Điều này giúp kiểm chứng độ ổn định liên kết thông qua chỉ số RMSD < 2,0 Å và xác nhận liên kết hydro bền vững theo thời gian.

4. Khả năng vượt qua hàng rào máu não của các chất tìm được có đảm bảo không?
Các phân tử tiềm năng được chọn lọc đều trải qua quy trình đánh giá khắt khe bằng công cụ pkCSM và PreADMET. Kết quả cho thấy các hoạt chất đều có chỉ số logPS > -3,0 và logBB dương (từ 0,11 đến 0,72), đáp ứng đầy đủ tiêu chuẩn thấm qua hàng rào máu não để phát huy tác dụng tại hệ thần kinh trung ương.

5. Luận văn đã sàng lọc từ những nguồn dữ liệu hóa học tự nhiên nào?
Nghiên cứu đã tiến hành sàng lọc trên 100 hợp chất nguồn gốc tự nhiên, bao gồm 50 phân tử trích xuất từ cơ sở dữ liệu COCONUT và 50 hợp chất hoạt tính có nguồn gốc từ các cây thuốc y học cổ truyền (như Thông đất, Tỏi, Nghệ, Saffron, Hương thảo).

Kết luận

  • Luận văn xây dựng thành công quy trình sàng lọc in silico đa tầng kết hợp học máy, docking phân tử, dự đoán ADME/BBB và mô phỏng động lực học phân tử 100 ns.
  • Bộ dữ liệu 1.055 hợp chất từ ChEMBL được làm sạch và phân loại chuẩn xác với 674 chất active (63,88%) và 381 chất inactive (36,12%).
  • Các thuật toán XGBoost và Random Forest thể hiện hiệu năng phân loại xuất sắc với ROC-AUC trên 0,87 và PRC-AUC xấp xỉ 0,90.
  • Phát hiện 8 cấu trúc tự nhiên tiềm năng có năng lượng gắn kết vượt trội (-10,2 đến -12,4 kcal/mol) và duy trì độ ổn định cao (RMSD < 2,0 Å) trong mô phỏng MD.
  • Đóng góp giải pháp tính toán khoa học đáng tin cậy, giúp định hướng hiệu quả cho các nghiên cứu tổng hợp hóa dược và thử nghiệm sinh học in vitro trong giai đoạn 6 đến 12 tháng tiếp theo.

Công trình là bước đệm quan trọng thúc đẩy ứng dụng trí tuệ nhân tạo trong hóa dược, mở ra hướng đi triển vọng cho công cuộc nghiên cứu thuốc điều trị bệnh Alzheimer tại Việt Nam.