Giới thiệu dự án

Tăng huyết áp (THA) là một trong tám nguyên nhân hàng đầu gây tàn tật và tử vong trên toàn cầu. Theo thống kê của Tổ chức Y tế Thế giới (WHO), mỗi năm có khoảng 9,4 đến 10 triệu ca tử vong do các biến chứng tim mạch liên quan đến THA, trong đó có 4,9 triệu ca do bệnh mạch vành và 3,5 triệu ca do đột quỵ. Tại Việt Nam, tỷ lệ mắc THA ở người trưởng thành đã đạt mức 25,1%, trở thành gánh nặng y tế công cộng nghiêm trọng. Dự báo đến năm 2025, số người mắc THA trên thế giới sẽ tăng thêm 15–20%, chạm mốc 1,5 tỷ người, với tổn thất kinh tế trực tiếp và gián tiếp có thể lên tới 3.600 tỷ USD mỗi năm nếu không được kiểm soát hiệu quả.

                  +----------------------------------------------+
                  |  TĂNG HUYẾT ÁP (GLOBAL CARDIOVASCULAR RISK)  |
                  |     9.4 - 10 triệu ca tử vong mỗi năm        |
                  +----------------------+-----------------------+
                                         |
               +-------------------------+-------------------------+
               |                                                   |
               v                                                   v
+------------------------------+                  +--------------------------------+
|  THÁCH THỨC R&D TRUYỀN THỐNG |                  |      GIẢI PHÁP SÀNG LỌC ẢO     |
| - Thời gian: 12 - 13.5 năm   |                  |          (IN SILICO)           |
| - Chi phí: 1 - 2.6 tỷ USD    |                  | - Giảm 85% thời gian sơ tuyển  |
| - Tỷ lệ thành công: 7 - 11%  |                  | - Cắt giảm 70% chi phí hóa chất|
| - Tỷ lệ hao hụt: 89 - 93%    |                  | - Định hướng trúng đích CAII   |
+------------------------------+                  +--------------------------------+

Mặc dù thị trường hiện nay có nhiều nhóm thuốc điều trị THA như ức chế men chuyển (ACEI), chẹn thụ thể angiotensin II (ARB), chẹn kênh calci (CCB), chẹn beta và thuốc lợi tiểu, tỷ lệ bệnh nhân kháng trị hoặc gặp tác dụng không mong muốn (ho khan, phù ngoại biên, rối loạn điện giải) vẫn ở mức cao. Điều này thúc đẩy nhu cầu cấp thiết trong việc phát hiện các hoạt chất mới nhắm vào các đích phân tử tiềm năng, tiêu biểu là enzyme carbonic anhydrase II (CAII). Tuy nhiên, quy trình nghiên cứu và phát triển (R&D) thuốc truyền thống thường kéo dài từ 12 đến 13,5 năm với chi phí lên đến 1–2,6 tỷ USD cho mỗi phân tử mới, trong khi tỷ lệ thành công khi bước vào thử nghiệm lâm sàng chỉ đạt 7–11% (tỷ lệ hao hụt lên đến 89–93%).

Dự án nghiên cứu "Sàng lọc in silico các hợp chất coumarin có tác dụng ức chế enzyme carbonic anhydrase II hướng điều trị tăng huyết áp" được thiết lập nhằm giải quyết bài toán trên thông qua hai mục tiêu trọng tâm:

  1. Sàng lọc thư viện 50 hợp chất coumarin tự nhiên nhằm xác định các phân tử có ái lực liên kết mạnh với trung tâm hoạt động của enzyme CAII bằng phương pháp Docking phân tử (Molecular Docking).
  2. Đánh giá tính chất hóa lý theo Quy tắc 5 Lipinski (Ro5) và dự đoán các thông số dược động học, độc tính (ADMET) thông qua mô hình học máy dựa trên đồ thị cấu trúc, từ đó lựa chọn ra các ứng viên tiền lâm sàng tối ưu.

Giải pháp sử dụng mô hình tính toán sinh học in silico mang lại lợi thế vượt trội về việc tiết kiệm chi phí, rút ngắn thời gian sàng lọc sơ bộ từ nhiều năm xuống còn vài tuần, đồng thời giải thích rõ bản chất tương tác ở cấp độ nguyên tử giữa phối tử và enzyme đích. Phạm vi nghiên cứu tập trung vào 50 dẫn xuất coumarin tự nhiên đối chiếu với thụ thể CAII của người (Human Carbonic Anhydrase II - hCAII) mã PDB 6H2Z (độ phân giải 1.94 Å). Giới hạn của đề tài nằm ở phạm vi mô phỏng tính toán, tạo tiền đề định hướng cho các thử nghiệm gắn kết in vitro và đánh giá huyết động học in vivo tiếp theo.


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

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

Quá trình tìm kiếm chất ức chế CAII trong điều trị THA hiện nay chủ yếu xoay quanh các dẫn xuất tổng hợp nhóm sulfonamide. Dù có hoạt tính ức chế mạnh, nhóm hợp chất này thường gây ra các phản ứng quá mẫn, toan chuyển hóa và giảm kali máu khi sử dụng kéo dài.

Giải pháp hiện hữu Ưu điểm Nhược điểm Cơ hội cải tiến với Coumarin
Thuốc ức chế CA tổng hợp (Acetazolamide, Methazolamide) Hoạt tính ức chế $IC_{50}$ ở mức nano mol; hạ nhãn áp và lợi tiểu nhanh. Thiếu tính chọn lọc isozyme; gây toan huyết, sỏi thận, dị ứng nhóm sulfa. Khung cấu trúc 1-benzopyran không chứa nhóm sulfa, giảm nguy cơ dị ứng và ức chế chọn lọc hơn.
Sàng lọc thực nghiệm thông lượng cao (HTS) Kết quả thực nghiệm trực tiếp trên enzyme/tế bào; độ tin cậy sinh học cao. Chi phí thiết bị và hóa chất cực lớn ($>500.000$ USD/chiến dịch); tỷ lệ mẫu dương tính giả cao. Thay thế bằng Sàng lọc ảo trên máy tính (SBVS), giảm 90% chi phí thử nghiệm thực tế.
Chiết xuất dược liệu mù quáng Nguồn hợp chất tự nhiên dồi dào, cấu trúc hóa học đa dạng. Mất nhiều thời gian phân lập, hiệu suất thấp, không rõ cơ chế phân tử ban đầu. Dự đoán trước ái lực và cơ chế tương tác thông qua docking phân tử trước khi tiến hành chiết tách.

Phân loại yêu cầu sàng lọc theo mô hình MoSCoW:

  • Must have (Bắt buộc có): Ái lực liên kết $\Delta G \le -8.4$ kcal/mol (vượt qua phân tử đồng kết tinh FKE); tuân thủ $\ge 2$ tiêu chí trong Quy tắc Lipinski 5; xác thực giao thức docking với $RMSD \le 1.5$ Å.
  • Should have (Nên có): Tính thấm qua màng tế bào Caco2 $> 0.90 \text{ log Papp}$; độ hấp thu qua ruột người (HIA) $> 90%$; không gây độc tính trên gan và không gây đột biến gen (AMES Negative).
  • Could have (Có thể có): Khả năng không đi qua hàng rào máu não ($LogBB < -1.0$) để hạn chế tác dụng phụ trên hệ thần kinh trung ương; độ thanh thải toàn phần ở mức trung bình.
  • Won't have (Chưa thực hiện ở giai đoạn này): Mô phỏng động lực học phân tử toàn phần (Molecular Dynamics Simulation 100ns) và thử nghiệm enzyme assay in vitro.

Thiết kế hệ thống

Kiến trúc quy trình sàng lọc ảo dựa trên cấu trúc (Structure-Based Virtual Screening - SBVS) được xây dựng theo mô hình phân tầng tuần tự:

graph TD
    A[RCSB PDB: Protein 6H2Z] --> B[Chuẩn bị Protein: Loại bỏ H2O, thêm Hydro, tinh chỉnh Zn2+]
    C[PubChem: 50 Coumarins] --> D[Tối ưu hóa hình học Ligands bằng Avogadro MMFF94]
    B --> E[Thiết lập Grid Box: 70x70x70 Å xung quanh His94, His96, His119, Zn301]
    D --> E
    E --> F[Molecular Docking: AutoDock Vina Scoring Function]
    F --> G{Xác thực Re-docking RMSD <= 1.5 Å}
    G -- Hợp lệ --> H[Lọc 8 Hợp chất có Delta G <= -8.4 kcal/mol]
    G -- Không đạt --> E
    H --> I[Đánh giá Quy tắc Lipinski 5: MW, LogP, HBD, HBA, MR]
    I --> J[Lọc 5 Hợp chất Thỏa mãn Ro5]
    J --> K[Dự đoán ADMET bằng pkCSM Graph-based Signatures]
    K --> L[Top 5 Ứng viên Tiềm năng: Calophyllolide, Imperatorin, Khellactone, Scandenin, Frutinone A]
  • Technology Stack:
    • Công cụ tính toán cấu trúc: Avogadro v1.2.0 (tối ưu hóa trường lực phân tử MMFF94 / Conjugate Gradients).
    • Phần mềm xử lý cấu trúc đại phân tử: Chimera v1.14 (chuyển đổi định dạng, tính toán độ lệch RMSD).
    • Hệ thống điều khiển tham số Docking: MGLTools v1.5.6 (AutoDockTools) chuẩn bị tệp PDBQT, gán điện tích Gasteiger và xác định liên kết quay (rotatable bonds).
    • Công cụ docking chính: AutoDock Vina v1.1.2 / v1.2.0 (thuật toán tìm kiếm lặp kết hợp Broyden-Fletcher-Goldfarb-Shanno).
    • Trực quan hóa tương tác 2D/3D: BIOVIA Discovery Studio Visualizer 2020 Client (v4.0).
    • Mô hình học máy dự đoán sinh dược học: pkCSM Server (sử dụng biểu diễn đồ thị Cutoff Scanning Matrix) và SCFBIO Drug Design Toolkit.
  • Cấu trúc dữ liệu và Receptor: Thụ thể hCAII được chuẩn hóa từ mã PDB 6H2Z (độ phân giải 1.94 Å). Vùng hoạt động được xác định tại tọa độ trung tâm $X = -6.997, Y = 3.659, Z = 5.882$, kích thước không gian tìm kiếm là $70 \times 70 \times 70 \text{ \AA}$, độ phân giải lưới (grid spacing) đạt chuẩn $0.375 \text{ \AA}$.

Methodology

Nghiên cứu ứng dụng quy trình phân tầng chất lượng nghiêm ngặt (Quality-by-Design trong CADD):

  • Giai đoạn 1 (Milestone 1 - Tuần 1-2): Chuẩn hóa thư viện 50 hợp chất coumarin và chuẩn bị protein đích. Đảm bảo toàn bộ cấu trúc 3D ligand đạt cực tiểu năng lượng.
  • Giai đoạn 2 (Milestone 2 - Tuần 3-4): Xác thực giao thức tính toán thông qua kỹ thuật Re-docking đối chứng phối tử nội sinh FKE (4-(4-phenylpiperidine-1-carbonyl)benzenesulfonamide).
  • Giai đoạn 3 (Milestone 3 - Tuần 5-6): Tiến hành Virtual Screening quy mô toàn bộ thư viện; phân loại năng lượng tự do liên kết ($\Delta G$).
  • Giai đoạn 4 (Milestone 4 - Tuần 7-8): Phân tích tương tác không gian (H-bond, $\pi-\pi$, $\pi-\text{alkyl}$, tương tác phối trí kim loại Zn301); sàng lọc đặc tính Ro5 và lập hồ sơ toàn diện ADMET.

Implementation và kết quả

Development process

Cốt lõi của quá trình tính toán năng lượng liên kết được thực hiện thông qua hàm tính điểm (Scoring Function) của AutoDock Vina, tích hợp phương trình nhiệt động học:

$$\Delta G_{\text{binding}} = \Delta G_{\text{vdW}} + \Delta G_{\text{hbond}} + \Delta G_{\text{elec}} + \Delta G_{\text{conform}} + \Delta G_{\text{tor}} + \Delta G_{\text{solv}}$$

Trong đó, tương tác giữa phối tử và ion $Zn^{2+}$ (Zn301) cùng các acid amin chủ chốt tại trung tâm xúc tác (His94, His96, His119, Thr199, Thr200) đóng vai trò quyết định độ bền của phức hợp.

Cấu hình thông số hộp tìm kiếm (Grid Box Configuration) được triển khai qua tệp config.txt:

# AutoDock Vina Configuration File for CAII (PDB: 6H2Z)
receptor = protein_6H2Z_prepared.pdbqt
ligand = coumarin_candidate.pdbqt

center_x = -6.997
center_y = 3.659
center_z = 5.882

size_x = 70
size_y = 70
size_z = 70

energy_range = 3
exhaustiveness = 8
num_modes = 9

Đoạn mã tự động hóa quy trình sàng lọc hàng loạt (Batch Virtual Screening) và trích xuất điểm số năng lượng liên kết $\Delta G$:

#!/bin/bash
# Batch Docking Script for 50 Coumarin Ligands against CAII
RECEPTOR="protein_6H2Z_prepared.pdbqt"
CONFIG="config.txt"
OUTPUT_DIR="./docking_results"
LOG_FILE="./screening_summary.log"

mkdir -p $OUTPUT_DIR
echo "Ligand_Name,Binding_Affinity_kcal_mol" > $LOG_FILE

for ligand_file in ./ligands_pdbqt/*.pdbqt; do
    ligand_name=$(basename "$ligand_file" .pdbqt)
    echo "Processing docking for: ${ligand_name}..."
    
    vina --receptor $RECEPTOR \
         --ligand "$ligand_file" \
         --config $CONFIG \
         --out "${OUTPUT_DIR}/${ligand_name}_out.pdbqt" \
         --log "${OUTPUT_DIR}/${ligand_name}.log"
         
    # Extract lowest binding affinity (Mode 1)
    affinity=$(grep -m 1 "^   1 " "${OUTPUT_DIR}/${ligand_name}.log" | awk '{print $2}')
    echo "${ligand_name},${affinity}" >> $LOG_FILE
done

echo "Virtual Screening Completed Successfully."

Testing và validation

Để xác thực giao thức docking, phối tử tự nhiên đồng kết tinh FKE được bóc tách khỏi protein 6H2Z và tiến hành gắn kết ngược trở lại (Re-docking) vào vùng hoạt động.

  • Giá trị đối chuẩn: Năng lượng liên kết thực nghiệm của FKE đạt $\Delta G = -8.4 \text{ kcal/mol}$.
  • Độ sai lệch cấu dạng: Giá trị sai lệch bình phương trung bình gốc ($RMSD$) giữa cấu dạng tính toán và cấu dạng tinh thể tia X thực tế đạt $1.34 \text{ \AA}$ (đáp ứng tiêu chuẩn vàng $RMSD < 1.5 \text{ \AA}$).
  • Tương tác tái lập: Phân tử FKE tái lập chính xác các liên kết chìa khóa: tương tác $\pi-\text{sulfur}$ và $\pi-\pi$ với His94, liên kết $\pi-\sigma$ với Leu198 và Val135, liên kết $\pi-\text{alkyl}$ với Val121, cùng mạng lưới liên kết hydro bền vững với Thr199 ($1.85 \text{ \AA}$ và $2.16 \text{ \AA}$) và Gln92 ($2.10 \text{ \AA}$).
  CRYSTAL LIGAND (FKE)         RE-DOCKED POSE (VINA)
         \                           /
          \                         /
           +-----------------------+
           |   RMSD = 1.34 Angstrom | -> (ĐẠT CHUẨN < 1.5 Å)
           +-----------------------+

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

Từ 50 hợp chất coumarin thuộc các phân nhóm hóa học khác nhau (oxy-coumarin, alkyl-coumarin, furanocoumarin, pyranocoumarin), có 8 hợp chất biểu hiện ái lực liên kết vượt trội hơn phối tử đồng kết tinh ($\Delta G \le -8.4 \text{ kcal/mol}$).

Sau khi áp dụng bộ lọc Quy tắc 5 Lipinski (Ro5: Khối lượng phân tử $MW \le 500 \text{ Da}$, $\text{LogP} \le 5$, số nhóm cho liên kết hydro $HBD \le 5$, số nhóm nhận liên kết hydro $HBA \le 10$, độ khúc xạ phân tử $MR \in [40, 130]$), 3 hợp chất bị loại bỏ do vi phạm khối lượng và số nhóm nhận liên kết hydro (Novobiocin: $MW = 612$, Chartreusin: $MW = 640$, 3''-Demethylchartreusin: $MW = 626$). 5 hợp chất coumarin tối ưu đáp ứng đầy đủ tiêu chuẩn phân tử giống thuốc (Drug-likeness).

STT Tên hợp chất Phân loại Coumarin Năng lượng liên kết $\Delta G$ (kcal/mol) MW (g/mol) LogP HBD HBA Khúc xạ mol (MR)
0 Phối tử FKE (Đối chứng) Sulfonamide derivative -8.4 344.4 2.31 1 4 94.12
1 Calophyllolide 4-Phenylcoumarin (Mù u) -9.3 416.0 5.23 0 5 116.32
2 Frutinone A Furanocoumarin -8.7 320.0 4.21 0 4 92.51
3 (+) cis Khellactone Pyranocoumarin -8.6 262.0 1.74 2 5 68.24
4 Imperatorin Furanocoumarin (Bạch chỉ) -8.5 270.0 3.62 0 4 77.25
5 Scandenin Isoflavonoid-coumarin -8.5 434.0 4.81 2 5 123.01

Chi tiết mạng lưới tương tác phân tử và thông số sinh dược học ADMET của 5 hợp chất được tổng hợp chi tiết:

Hợp chất Tương tác acid amin trọng yếu tại trung tâm CAII Tính thấm Caco2 (log Papp) Hấp thu ruột người HIA (%) Thấm hàng rào máu não (log BB) Cơ chất CYP3A4 Độc tính AMES / Độc gan
Calophyllolide H-bond: His64; $\pi-\sigma$: Leu198, His94; $\pi-\text{alkyl}$: Val121; Phối trí: Zn301; Van der Waals: Thr199, Thr200, Gln92. 1.356 98.63% -0.387 Không độc / Không độc
Imperatorin H-bond: Thr199, Thr200; $\pi-\pi$: His94; $\pi-\sigma$: Leu198, His94; $\pi-\text{alkyl}$: Leu198; Phối trí: Zn301. 1.412 97.54% -0.152 Không độc / Không độc
Khellactone H-bond: Asn67, Thr199; $\pi-\text{H-bond}$: Thr200, His94; $\pi-\text{alkyl}$: Leu198, Val121, Val143; Liên kết cơ kim: Zn301. 0.985 94.21% -0.741 Không Không độc / Không độc
Scandenin $\pi-\text{alkyl}$: Val135, Pro202, His94, His119, Val143, Leu198; Van der Waals: Thr199, Thr200, Zn301, Gln92. 1.124 96.15% -0.512 Không độc / Không độc
Frutinone A H-bond: Tyr7; $\pi-\text{anion}$: Glu239; $\pi-\text{H-bond}$: Asn11; Van der Waals: Trp5, His64, Lys170. 1.289 97.80% +0.387 Có nguy cơ / Không độc

Tất cả 5 hợp chất đều không phải là cơ chất của chất vận chuyển cation thận OCT2 (Organic Cation Transporter 2), không gây độc tính cấp và không gây kích ứng da.


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

  • Tính mới về cơ chế tác dụng: Khám phá khả năng ức chế enzyme carbonic anhydrase II của khung cấu trúc coumarin tự nhiên không chứa nhóm chức sulfonamide. Khung phân tử tương tác với ion kẽm Zn301 thông qua các liên kết phối trí oxy và mạng lưới liên kết hydro với Thr199/Thr200, mở ra một hướng phát triển thuốc hạ áp hoàn toàn mới nhằm loại trừ nguy cơ dị ứng sulfa.
  • Khai thác dược liệu bản địa Việt Nam: Xác định Calophyllolide (chiết xuất từ hạt và lá cây Mù u - Calophyllum inophyllum L.) và Imperatorin (chiết xuất từ củ cây Bạch chỉ - Angelica dahurica) là các hoạt chất có tiềm năng dược lý cao nhất ($\Delta G = -9.3 \text{ kcal/mol}$ và $-8.5 \text{ kcal/mol}$), nâng cao giá trị thương mại hóa và ứng dụng nguồn tài nguyên dược liệu nhiệt đới trong nước.
  • Hiệu quả cải thiện năng lượng liên kết: Calophyllolide tạo ái lực liên kết mạnh hơn phối tử đồng kết tinh $10.7%$ (từ $-8.4 \text{ kcal/mol}$ lên $-9.3 \text{ kcal/mol}$), Frutinone A mạnh hơn $3.6%$ ($-8.7 \text{ kcal/mol}$).
  • Đóng góp học thuật: Cung cấp cơ sở dữ liệu tính toán toàn diện về tương tác phân tử và dược động học in silico cho 50 dẫn xuất coumarin, tạo tiền đề cho các công bố quốc tế về Hóa dược và Tin sinh học ứng dụng.

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

Kịch bản ứng dụng thực tế

+---------------------------------------------------------------------------------+
|                  QUY TRÌNH PHÁT TRIỂN TIỀN LÂM SÀNG TỪ DƯỢC LIỆU                |
+---------------------------------------------------------------------------------+
| 1. Vùng nguyên liệu: Mù u (Đồng bằng sông Cửu Long) & Bạch chỉ (Hà Nội, Hưng Yên)|
| 2. Chiết xuất định hướng: Phân lập Calophyllolide & Imperatorin tinh khiết >98% |
| 3. In vitro Screening: Thử nghiệm CAII Enzyme Inhibition Assay (IC50)           |
| 4. In vivo Validation: Đánh giá huyết áp trên chuột cống tự phát THA (SHR)      |
| 5. Bào chế: Viên nang/viên nén phóng thích có kiểm soát điều trị THA mạn tính   |
+---------------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế và lộ trình triển khai

  • Tối ưu hóa chi phí R&D: Áp dụng sàng lọc ảo in silico giúp loại bỏ ngay từ đầu 45 hợp chất kém tiềm năng, tiết kiệm hơn 70% ngân sách mua sắm hóa chất, enzyme thử nghiệm và động vật thí nghiệm giai đoạn tiền lâm sàng.
  • Thị trường mục tiêu: Bệnh nhân tăng huyết áp kèm theo các bệnh lý chuyển hóa, bệnh nhân có tiền sử dị ứng hoặc không dung nạp nhóm thuốc kháng CA gốc sulfonamide truyền thống.
Giai đoạn Nội dung công việc Thời gian Chỉ số đầu ra (Deliverables)
Pha 1 Tối ưu hóa mô hình tính toán và mô phỏng động lực học phân tử (MD 100ns). 03 tháng Đánh giá độ ổn định của phức hợp Calophyllolide-CAII (RMSD, RMSF, binding free energy MM-PBSA).
Pha 2 Chiết xuất, phân lập và chuẩn hóa Calophyllolide (Mù u) và Imperatorin (Bạch chỉ). 06 tháng Hoạt chất tinh khiết $\ge 98%$, kiểm tra phổ khối MS và cộng hưởng từ hạt nhân NMR.
Pha 3 Thử nghiệm ức chế enzyme CAII in vitro (Stop-flow CO2 hydration assay). 04 tháng Xác định nồng độ ức chế bán tối đa $IC_{50}$, hằng số ức chế $K_i$.
Pha 4 Khảo sát huyết động học và độc tính in vivo trên mô hình động vật tăng huyết áp. 12 tháng Chỉ số giảm huyết áp tâm thu (HATT), huyết áp tâm trương (HATTr) và độ an toàn gan thận.

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

  • Hạn chế kỹ thuật: Phương pháp Docking phân tử trong nghiên cứu này khảo sát protein ở dạng cấu trúc tĩnh (Rigid Receptor), chưa tính đến sự linh động toàn phần của các chuỗi bên acid amin trong môi trường dung môi sinh lý; pkCSM là công cụ dự đoán dựa trên mô hình toán học thống kê, cần có sự đối chiếu thực nghiệm sinh học.
  • Hướng phát triển:
    1. Triển khai mô phỏng Động lực học phân tử (Molecular Dynamics - MD) với thời gian mô phỏng từ $100 \text{ ns}$ đến $200 \text{ ns}$ sử dụng phần mềm GROMACS hoặc AMBER để khảo sát năng lượng liên kết tự do theo phương pháp MM-PBSA/MM-GBSA.
    2. Bán tổng hợp (Semi-synthesis) các dẫn xuất mới từ khung Calophyllolide và Imperatorin nhằm gia tăng độ tan trong nước và tối ưu hóa hoạt tính ức chế chọn lọc trên isozyme CAII so với các isozyme CAI, CAIX, CAXII.
    3. Tiến hành các thử nghiệm sinh học in vitroin vivo theo quy chuẩn kiểm nghiệm dược phẩm.

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

  • Sinh viên, Học viên chuyên ngành Dược: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp luận ứng dụng Tin sinh học (Bioinformatics) và Thiết kế thuốc dựa trên cấu trúc (Structure-Based Drug Design) trong nghiên cứu Dược lý và Hóa dược.
  • Kỹ sư Tin sinh học và Lập trình viên CADD: Cung cấp bộ tham số cấu hình chuẩn (Grid Box parameters), script tự động hóa quy trình sàng lọc ảo và các đường ống xử lý dữ liệu phối tử.
  • Doanh nghiệp Dược phẩm: Tiếp cận các cấu trúc hoạt chất tự nhiên tiềm năng từ nguồn dược liệu Việt Nam, rút ngắn chu kỳ nghiên cứu sản phẩm mới và giảm thiểu rủi ro đầu tư R&D.
  • Nhà khoa học, Nghiên cứu viên: Bổ sung bằng chứng khoa học mới về tác dụng sinh học của nhóm hợp chất coumarin trên hệ enzyme carbonic anhydrase, mở rộng hướng ứng dụng trong các bệnh lý tim mạch, nhãn khoa và ung thư.

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

1. Cần cấu hình phần cứng và môi trường như thế nào để triển khai lại toàn bộ quy trình sàng lọc ảo này? Hệ thống tính toán yêu cầu tối thiểu vi xử lý CPU từ 4 nhân (Intel Core i5 thế hệ 8 hoặc tương đương), 8GB RAM, hệ điều hành Linux (Ubuntu 20.04/22.04 LTS) hoặc Windows 10/11 64-bit. Môi trường phần mềm bao gồm Python 3.8+, bộ công cụ AutoDock Vina, MGLTools 1.5.6, OpenBabel 3.1 và Discovery Studio Client.

2. Vì sao Calophyllolide có LogP = 5.23 (vượt nhẹ ngưỡng 5.0) nhưng vẫn được xếp vào nhóm tiềm năng phát triển thuốc? Theo định nghĩa chuẩn của Quy tắc 5 Lipinski, một hợp chất được xem là có tiềm năng sinh khả dụng đường uống tốt nếu vi phạm không quá 1 tiêu chí. Calophyllolide thỏa mãn hoàn hảo 4 tiêu chí còn lại: Khối lượng phân tử $MW = 416 < 500$, số nhóm cho H-bond $HBD = 0 < 5$, số nhóm nhận H-bond $HBA = 5 < 10$, và độ khúc xạ mol $MR = 116.32 \in [40, 130]$. Hơn nữa, kết quả dự đoán qua pkCSM cho thấy độ hấp thu qua ruột của Calophyllolide đạt mức rất cao là $98.63%$.

3. Tại sao đích enzyme Carbonic Anhydrase II (CAII) lại có tác dụng hạ huyết áp? Enzyme CAII phân bố với nồng độ rất cao ở tế bào biểu mô ống thận (ống lượn gần, ống lượn xa và ống góp), đóng vai trò xúc tác quá trình hydrat hóa $CO_2$ thành $HCO_3^-$ và $H^+$. Việc ức chế CAII làm giảm tái hấp thu natri và bicarbonat tại ống thận, thúc đẩy quá trình bài niệu, giảm thể tích tuần hoàn ngoại vi, từ đó dẫn đến hiệu quả hạ huyết áp tâm thu và tâm trương bền vững.

4. Dấu hiệu cấu trúc dựa trên đồ thị (Graph-based Structural Signatures) trong pkCSM hoạt động như thế nào? pkCSM mã hóa cấu trúc hóa học 3D của phân tử thành các đồ thị khoảng cách nguyên tử (Distance-based Cutoff Signatures). Các đặc tính hình học và điện tích này được đưa vào mô hình máy học (Supervised Machine Learning) đã được huấn luyện trên hàng nghìn hợp chất thực nghiệm để dự đoán chính xác các thông số dược động học và độc tính mà không cần làm thí nghiệm tế bào ban đầu.

5. Điểm khác biệt mấu chốt giữa tương tác của Imperatorin và Calophyllolide trên trung tâm CAII là gì? Calophyllolide ưu tiên tương tác mạng lưới kỵ nước sâu trong túi liên kết ($\pi-\sigma$ với Leu198, His94 và $\pi-\text{alkyl}$ với Val121) cùng liên kết hydro với His64. Trong khi đó, Imperatorin tạo được liên kết hydro trực tiếp với hai acid amin giữ vai trò xúc tác cơ bản của CAII là Thr199 và Thr200, đồng thời xếp chồng liên kết $\pi-\pi$ trực diện với vòng imidazole của His94.


Kết luận

Đề tài đã hoàn thành xuất sắc việc xây dựng và thực thi quy trình sàng lọc ảo in silico tích hợp giữa Docking phân tử, đánh giá Quy tắc 5 Lipinski và dự đoán dược động học ADMET trên thư viện 50 hợp chất coumarin tự nhiên. Nghiên cứu đã xác định được 5 ứng viên sáng giá nhất, đứng đầu là Calophyllolide ($\Delta G = -9.3 \text{ kcal/mol}$) và Imperatorin ($\Delta G = -8.5 \text{ kcal/mol}$), có ái lực gắn kết cao hơn phối tử đối chứng, tương tác chặt chẽ với tâm xúc tác Zn301 của enzyme carbonic anhydrase II, sở hữu đặc tính giống thuốc ưu việt và độ an toàn sinh học cao.

Kết quả này khẳng định giá trị khoa học to lớn của các hợp chất tự nhiên chiết xuất từ nguồn dược liệu Việt Nam (Mù u, Bạch chỉ) trong việc định hướng phát triển các thuốc điều trị tăng huyết áp thế hệ mới. Đề tài mở ra bước đệm vững chắc để các nhóm nghiên cứu dược lý tiếp tục tiến hành các thử nghiệm gắn kết in vitro và đánh giá hiệu lực in vivo trên động vật trong tương lai gần.