Giới thiệu dự án

Bệnh Alzheimer (Alzheimer's Disease – AD) là dạng sa sút trí tuệ phổ biến nhất thế giới, chiếm 60–80% tổng số ca mất trí nhớ ở người cao tuổi. Theo Tổ chức Y tế Thế giới (WHO), hiện có khoảng 50 triệu bệnh nhân trên toàn cầu, dự báo tăng lên 152 triệu vào năm 2050, tương đương tăng gấp đôi mỗi 5 năm. Tại Việt Nam, tỷ lệ mắc là 5% ở nhóm tuổi 65 và 20% ở nhóm trên 85 tuổi; riêng ảnh hưởng từ phơi nhiễm chất độc chiến tranh ước tính khiến 140.000 cựu chiến binh mắc bệnh vào năm 2020.

Problem statement: Các thuốc đang được FDA chấp thuận — Donepezil, Rivastigmine, Galantamine và Memantine — chỉ giải quyết triệu chứng mà không tác động vào cơ chế bệnh sinh sâu xa. Không có lý thuyết cơ chế bệnh sinh thống nhất, và các ứng viên thuốc như solanezumab, verubecestat hay atabecestat đều thất bại ở giai đoạn thử nghiệm lâm sàng. Khoảng trống nghiên cứu rõ ràng: cần một hướng tiếp cận đa đích (multi-target) nhắm đồng thời vào NMDA, BACE-1, MAO A và AChE.

Mục tiêu nghiên cứu:

  1. Sàng lọc 50 hợp chất từ cây Rau đắng biển (Bacopa monnieri (L.) Wettst) có khả năng ức chế đồng thời ≥3 trong 4 đích phân tử: NMDA, BACE-1, MAO A, AChE bằng phương pháp docking phân tử (Molecular Docking)
  2. Đánh giá đặc điểm giống thuốc theo Quy tắc 5 của Lipinski (Rule of 5 – RO5) cho các ứng viên hàng đầu
  3. Dự đoán các thông số dược động học và độc tính (ADMET) bằng mô hình in silico

Kết quả mong đợi: Xác định ≥3 hợp chất từ Bacopa monnieri thỏa mãn RO5, không có độc tính AMES, qua được hàng rào máu não (Blood-Brain Barrier – BBB) và có năng lượng liên kết thấp hơn chất đối chứng đã được phê duyệt.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

So sánh giải pháp hiện có

Tiêu chí Thuốc đơn đích hiện tại Phương pháp in vitro/in vivo Sàng lọc ảo (in silico) – Dự án này
Chi phí Thấp (đã có) Rất cao Thấp
Thời gian Ngắn Dài (năm) Ngắn (tuần)
Phạm vi khảo sát 1 đích/nghiên cứu Hạn chế 50 hợp chất × 4 đích
Độ độc với người Đã biết Cần thử nghiệm Dự đoán in silico
Tiếp cận đa đích Không Hạn chế

Phân tích yêu cầu theo MoSCoW

  • Must have: Năng lượng liên kết (ΔG) ≤ chất đối chứng; ức chế ≥3 đích; RMSD ≤ 1.5 Å trong re-dock
  • Should have: Thỏa mãn ≥2 tiêu chí RO5; không có độc tính AMES; log BB > −1
  • Could have: Qua BBB tốt (log BB > 0.3); không ức chế CYP3A4/CYP2D6
  • Won't have (lần này): Thử nghiệm in vitro, in vivo hoặc lâm sàng

Thiết kế hệ thống

Technology Stack

Công cụ/Phần mềm Phiên bản Mục đích
AutoDock Vina 1.1.2 Engine docking phân tử chính
AutoDock Tools 1.5.6 Chuẩn bị protein/phối tử (.pdbqt)
Discovery Studio 2021 Client Biểu diễn tương tác 2D/3D protein–ligand
Avogadro Latest Tối ưu hóa cấu trúc 3D phối tử (Conjugate Gradients)
UCSF Chimera 1.x Tính RMSD, đánh giá chồng khít cấu trúc
pkCSM (online) Dự đoán thông số ADMET
SCFBio (online) Tính toán các thông số RO5
PubChem Nguồn cấu trúc phân tử (SMILES, .sdf)

Các đích phân tử và PDB ID

Đích PDB ID Phối tử đồng kết tinh Vai trò trong Alzheimer
NMDA receptor 1PBQ 5,7-dichlorokynurenic acid Ngăn kích thích thần kinh do glutamate
BACE-1 (β-secretase) 4X7I LY2886721 Ức chế tạo Aβ peptide
MAO A 2Z5X Harmine Giảm ROS, tăng dẫn truyền thần kinh
AChE 4EY7 Donepezil Bảo tồn acetylcholine (ACh) trong khớp thần kinh

Workflow tổng thể

PubChem → 50 hợp chất Bacopa + 4 chất đối chứng
    ↓
Chuẩn bị protein (loại H₂O, gắn H, trường lực Kollman)
    ↓
Re-dock phối tử đồng kết tinh → Xác nhận RMSD ≤ 1.5 Å
    ↓
Docking 50 hợp chất × 4 protein đích (AutoDock Vina)
    ↓
Lọc: ΔG ≤ ΔG(đối chứng) + ức chế ≥3 đích → 15 hợp chất
    ↓
Sàng lọc RO5 (SCFBio) → 3 hợp chất giống thuốc
    ↓
Phân tích ADMET (pkCSM) → 3 ứng viên tiềm năng

Methodology

  • Phương pháp phát triển: Nghiên cứu sàng lọc ảo (Virtual Screening) theo quy trình 3 bước tuần tự (Waterfall): docking → RO5 → ADMET
  • Thuật toán tìm kiếm: Gradient Descent kết hợp với Iterated Local Search trong AutoDock Vina; số vòng lặp mặc định = 8; khoảng cách lưới (grid spacing) = 0.375 Å
  • Hàm tính điểm (Scoring function): Hàm thực nghiệm mặc định của AutoDock Vina ước tính năng lượng tự do Gibbs (ΔG, kCal/mol)

Implementation và kết quả

Development process

Pha 1 – Chuẩn bị và xác thực quy trình: Tải 4 cấu trúc protein từ Protein Data Bank (PDB). Loại bỏ phân tử nước và các tinh thể tạp; bổ sung hydro thiếu; gắn trường lực Kollman bằng AutoDock Tools. Thiết lập grid box dựa trên vị trí phối tử đồng kết tinh:

NMDA (1PBQ):  grid size = 88×82×? Å,   center x=34.09
BACE-1(4X7I): grid size = 80×54×62 Å,  center x=27.58
MAO A (2Z5X): grid size = 55×?×12 Å,   center x=2.65
AChE (4EY7):  grid size = 59×?×51 Å,   center x=-2.86

Pha 2 – Re-dock và xác nhận RMSD:

Protein   | Phối tử đồng kết tinh          | ΔG (kCal/mol) | RMSD (Å)
----------|-------------------------------|---------------|----------
NMDA      | 5,7-dichlorokynurenic acid     | -5.687        | < 1.5
BACE-1    | LY2886721                     | -9.261        | < 1.5
MAO A     | Harmine                       | -6.x          | < 1.5
AChE      | Donepezil                     | -12.990       | < 1.5

Tất cả RMSD < 1.5 Å → quy trình docking được xác nhận là đáng tin cậy.

Pha 3 – Sàng lọc 50 hợp chất: Tối ưu hóa cấu trúc 3D phối tử bằng phương pháp Conjugate Gradients trong Avogadro, chuyển sang định dạng .pdbqt, thực hiện docking với AutoDock Vina.

Testing và validation

Tiêu chí lọc pha 1 (docking):

  • ΔG ≤ ΔG(chất đối chứng) tương ứng
  • Có khả năng ức chế ≥3 trong 4 đích

Kết quả: 15/50 hợp chất đáp ứng tiêu chí.

Tiêu chí lọc pha 2 (RO5):

Hợp chất MW (Da) HBD HBA logP MR Giống thuốc
Apigenin 270 3 5 2.195 ✓(40–130) Có (5/5)
Luteolin ~286 4 6 ~1.9 Có (5/5)
Ebelin lactone 454 1 3 7.479 Có (vi phạm logP)
Bacopaside II 928 10 18 1.077 Không
Bacopaside III 846 7 16 2.713 Không

Kết quả: 3/15 hợp chất vượt qua sàng lọc RO5.

Kết quả ADMET

Thông số Apigenin Ebelin lactone Luteolin
Caco2 (log Papp) 1.13 ✓ cao
HIA Cao (>80%) Cao Cao
log BB (BBB) -0.907 Trong ngưỡng
log PS (CNS) -2.251 -1.x (tốt nhất)
log VDss 0.153
CYP3A4 (cơ chất) Không Không
CYP2D6 (ức chế) Không Không
AMES toxicity Không Không Không
Độc tính gan Không Không Không
hERG inhibition Không Không Không

Kết quả đạt được

  • 3/3 hợp chất tiềm năng được xác định (Apigenin, Luteolin, Ebelin lactone)
  • Apigenin và Luteolin thỏa mãn 5/5 tiêu chí RO5 — hiệu suất 100% so với quy tắc giống thuốc
  • Không có hợp chất nào có độc tính AMES, độc tính gan hay ức chế kênh kali tim (hERG)
  • Quy trình sàng lọc giảm từ 50 xuống còn 3 ứng viên, tỷ lệ lọc đạt 94%

Đổi mới và đóng góp

Đổi mới kỹ thuật

1. Tiếp cận đa đích (Multi-target Drug Design): Thay vì sàng lọc theo một đích duy nhất như phần lớn nghiên cứu trước, đề tài này đánh giá đồng thời 4 đích: NMDA, BACE-1, MAO A và AChE trong một quy trình sàng lọc liên tục. Apigenin được phát hiện tạo tương tác liên kết đặc trưng với cả 4 đích:

  • NMDA: liên kết hydro với Arg131; liên kết π–anion với Asp224; liên kết π–σ với Ser180
  • BACE-1: liên kết π–π với Tyr71; liên kết hydro với Thr329
  • MAO A: liên kết alkyl/π–alkyl với Leu259, Lys280, Ile273
  • AChE: liên kết π–π xếp chồng với Trp86, His447, Trp286

2. Ứng dụng thực vật học Ayurvedic vào sàng lọc ảo hiện đại: Bacopa monnieri đã được sử dụng trong y học Ayurvedic hơn 3.000 năm, nhưng đây là một trong các nghiên cứu đầu tiên tại Việt Nam hệ thống hóa 50 hợp chất của cây bằng in silico multi-target docking.

3. Phát hiện Ebelin lactone như ứng viên mới: Ebelin lactone (MW = 454 Da, logP = 7.479) không thỏa mãn hoàn toàn RO5 nhưng có năng lượng liên kết cạnh tranh và ADMET không độc — mở ra hướng tối ưu hóa cấu trúc (structural optimization) để giảm logP về ngưỡng chấp nhận (<5).

So sánh với prior art

Nghiên cứu Đích Phương pháp Số hợp chất Đổi mới của đề tài
Bệnh nhân Alzheimer điều trị bằng Donepezil (FDA, 1996) AChE duy nhất Lâm sàng 1 Sàng lọc 50 hợp chất, 4 đích
Nghiên cứu BACE-1 inhibitor (LY2886721, lâm sàng pha 2) BACE-1 duy nhất Synthetic chemistry 1 Multi-target; nguyên liệu tự nhiên
Đề tài này 4 đích đồng thời In silico docking + RO5 + ADMET 50 Tiết kiệm chi phí; định hướng đa đích

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng

  • Phòng thí nghiệm dược học: Apigenin và Luteolin có thể được đưa vào thử nghiệm in vitro ức chế AChE và MAO A ngay với chi phí thấp vì cả hai flavonoid này đều có sẵn thương mại.
  • Công nghiệp dược phẩm: Ebelin lactone có thể là khung cấu trúc (scaffold) để tổng hợp dẫn xuất mới với logP được điều chỉnh.
  • Nghiên cứu dược liệu học: Kết quả định hướng chiết xuất chuẩn hóa Bacopa monnieri tập trung vào phân đoạn giàu Apigenin/Luteolin.

Yêu cầu triển khai thực tế

  • Thử nghiệm in vitro: Xác định IC₅₀ của Apigenin và Luteolin trên enzyme AChE tái tổ hợp (recombinant AChE)
  • Thử nghiệm in vivo: Mô hình chuột chuyển gen APP/PS1 để đánh giá giảm mảng Aβ
  • Chi phí ước tính giai đoạn tiếp theo: Thử nghiệm in vitro ~50–100 triệu VNĐ; in vivo pha tiền lâm sàng ~500–1.000 triệu VNĐ

Tiềm năng thị trường

Thị trường thuốc Alzheimer toàn cầu ước đạt 13,7 tỷ USD vào 2025. Một thuốc có cơ chế đa đích từ nguồn thảo dược có thể định vị trong phân khúc "functional nutraceutical" hoặc ứng viên thuốc kê đơn mới với lợi thế nguyên liệu sẵn có ở Việt Nam.


Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Phương pháp docking phân tử là tĩnh — không mô phỏng sự linh hoạt của protein (protein flexibility) hay hiệu ứng dung môi (solvation effect) đầy đủ. Cần bổ sung bằng mô phỏng Molecular Dynamics (MD) để đánh giá độ bền phức hợp protein–ligand theo thời gian.
  • Các thông số ADMET dự đoán in silico (pkCSM) có sai số so với thực nghiệm, đặc biệt đối với thấm BBB — cần xác nhận bằng mô hình tế bào MDCK hoặc PAMPA.
  • Ebelin lactone vi phạm tiêu chí logP (7.479 > 5) → cần tối ưu hóa cấu trúc.

Hướng phát triển

  1. MD Simulation: Chạy mô phỏng động lực học phân tử 100 ns cho Apigenin–AChE và Luteolin–BACE-1 để xác nhận độ bền liên kết
  2. QSAR Modeling: Xây dựng mô hình định lượng cấu trúc–hoạt tính (Quantitative Structure–Activity Relationship) từ 50 hợp chất đã docking
  3. Thực nghiệm xác nhận: Đo IC₅₀ của Apigenin, Luteolin với enzyme AChE và MAO A
  4. Tối ưu hóa Ebelin lactone: Thiết kế dẫn xuất giảm logP bằng cách thêm nhóm hydroxyl hoặc amine phân cực

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

  • Sinh viên Dược/Y: Tài liệu thực hành sàng lọc ảo với quy trình đầy đủ từ PDB → docking → RO5 → ADMET; sử dụng phần mềm miễn phí hoặc chi phí thấp
  • Nhà nghiên cứu dược lý: Framework chuẩn để mở rộng sang sàng lọc 200–500 hợp chất từ các thảo dược Việt Nam khác; có thể tái sử dụng toàn bộ pipeline AutoDock Vina
  • Doanh nghiệp dược phẩm: Định hướng R&D từ Bacopa monnieri — ước tính rút ngắn giai đoạn hit identification từ 2–3 năm xuống còn 3–6 tháng nhờ sàng lọc ảo
  • Nhà hoạch định y tế: Bằng chứng khoa học bước đầu hỗ trợ chuẩn hóa và phát triển sản phẩm từ cây Rau đắng biển trong danh mục thuốc quốc gia

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

1. Yêu cầu kỹ thuật để tái hiện nghiên cứu? Máy tính Windows 10 (ví dụ: Dell Latitude E7450 trở lên), AutoDock Vina 1.1.2, AutoDock Tools 1.5.6, Discovery Studio 2021 Client (free), Avogadro, Chimera 1.x, kết nối internet để truy cập PubChem, pkCSM, SCFBio.

2. Tại sao chọn AutoDock Vina thay vì Glide hay GOLD? AutoDock Vina là phần mềm mã nguồn mở, miễn phí, có tốc độ tính toán nhanh hơn AutoDock 4 khoảng 100 lần, phù hợp với quy mô nghiên cứu sinh viên và đã được xác nhận rộng rãi trong cộng đồng docking phân tử.

3. Kết quả có thể mở rộng sang bao nhiêu hợp chất? Pipeline hiện tại có thể mở rộng lên hàng nghìn hợp chất bằng cách sử dụng scripting Python kết hợp AutoDock Vina batch mode và cluster computing.

4. Làm thế nào tích hợp với hệ thống nghiên cứu lâm sàng hiện tại? Kết quả in silico cần được xác nhận qua chuỗi: in vitro enzyme assay → in vitro cell assay → in vivo animal model → Phase I clinical trial. Đây là bước hit identification trong chuỗi Drug Discovery Pipeline.

5. Chi phí và ROI so với phương pháp HTS (High-Throughput Screening)? Sàng lọc ảo in silico chi phí gần như bằng 0 (phần mềm miễn phí, điện và thời gian). HTS truyền thống cho 50 hợp chất × 4 đích ước tính 200–400 triệu VNĐ. Tiết kiệm >95% chi phí ở giai đoạn hit identification.


Kết luận

Khóa luận đã thực hiện thành công quy trình sàng lọc ảo in silico ba bước cho 50 hợp chất từ Bacopa monnieri trên 4 đích phân tử liên quan đến bệnh Alzheimer. Kết quả nổi bật:

  • ApigeninLuteolin — hai flavonoid thỏa mãn toàn bộ 5 tiêu chí RO5, không có độc tính AMES/gan/tim, qua được hàng rào máu não, và có năng lượng liên kết cạnh tranh so với thuốc đối chứng FDA
  • Ebelin lactone — ứng viên tiềm năng với ADMET tốt, cần tối ưu hóa cấu trúc để giảm logP
  • Quy trình sàng lọc đa đích giúp tăng hiệu quả phát hiện ứng viên thuốc, đặc biệt có giá trị trong bối cảnh nguồn nguyên liệu thảo dược phong phú của Việt Nam

Nghiên cứu khẳng định tiềm năng của cây Rau đắng biển trong điều trị bệnh Alzheimer và cung cấp nền tảng khoa học để tiếp tục đánh giá thực nghiệm in vitroin vivo, từng bước hướng tới phát triển thuốc từ thực vật có tác dụng đa đích trong điều trị bệnh thoái hóa thần kinh.