Giới thiệu dự án

Đại dịch COVID-19 do virus SARS-CoV-2 gây ra đã tạo nên một cuộc khủng hoảng y tế toàn cầu chưa từng có tiền lệ. Theo số liệu thống kê từ Tổ chức Y tế Thế giới (WHO), tính đến tháng 4/2022, toàn cầu đã ghi nhận hơn 487 triệu ca nhiễm và trên 6,1 triệu trường hợp tử vong. Tại Việt Nam, số ca nhiễm đã vượt mốc 10 triệu với hơn 42.000 ca tử vong. Mặc dù tỷ lệ bao phủ vaccine đạt mức cao (với 79,9% dân số Việt Nam tiêm đủ liều cơ bản tính đến tháng 3/2022), sự xuất hiện liên tục của các biến chủng mới cùng hội chứng "hậu COVID-19" (Long COVID) đặt ra yêu cầu cấp bách về việc phát triển các liệu pháp kháng virus đặc hiệu.

Quy trình nghiên cứu và phát triển (R&D) một loại thuốc mới theo phương pháp thực nghiệm truyền thống (de novo drug discovery) thường kéo dài từ 10–15 năm với chi phí trung bình ước tính từ 1 đến 2,6 tỷ USD. Trong bối cảnh đại dịch diễn biến nhanh, phương pháp sàng lọc ảo dựa trên cấu trúc (Structure-Based Virtual Screening - SBVS) và mô phỏng in silico trở thành giải pháp đột phá. Phương pháp này cho phép rút ngắn thời gian sàng lọc hàng trăm lần, tối ưu hóa chi phí nghiên cứu tiền lâm sàng và giảm thiểu rủi ro thất bại sinh học nhờ dự đoán chính xác ái lực liên kết thụ thể - phối tử cùng các đặc tính dược động học ADMET.

+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH XÂM NHẬP VÀ SAO CHÉP SARS-COV-2                    |
+-----------------------------------------------------------------------------------+
|  [SARS-CoV-2 Spike] ---> Liên kết ---> [Host ACE2 Receptor]                       |
|           |                                                                       |
|           v (Hòa màng / Thực bào qua TMPRSS2 & Cathepsin L)                       |
|  [Giải phóng viral RNA (+)] ---> Dịch mã Polyproteins (pp1a, pp1ab)               |
|           |                                                                       |
|           v (Cắt phân đoạn nhờ Mpro / 3CLpro & PLpro)                             |
|  [Các protein phi cấu trúc (nsps)] ---> Hình thành phức hợp nsp12-nsp7-nsp8 (RdRp)|
|           |                                                                       |
|           v (Sao chép và phiên mã bộ gen RNA)                                     |
|  [Virion mới lắp ráp và xuất bào]                                                 |
+-----------------------------------------------------------------------------------+

Dự án "Nghiên cứu tác dụng ức chế các đích phân tử virus SARS-CoV-2 của các hợp chất có trong cây Cam thảo (Glycyrrhiza glabra L.) bằng phương pháp in silico" được thực hiện nhằm giải quyết bài toán cấp thiết trên. Đề tài tập trung vào 4 mục tiêu cốt lõi:

  1. Sàng lọc phân tử đồng thời trên 4 đích sinh học: Men chuyển Angiotensin 2 (ACE2 - PDB: 1R4L), Protein gai (Spike Protein RBD - PDB: 6M0J), Protease chính (Main Protease / $M^{\text{pro}}$ - PDB: 7L0D) và RNA polymerase phụ thuộc RNA (RdRp - PDB: 7BV2).
  2. Thẩm định độ tin cậy quy trình docking (Re-docking validation): Đảm bảo độ lệch bình phương trung bình gốc (Root Mean Square Deviation - RMSD) của phối tử đồng kết tinh đạt giá trị $< 1,5\text{ \AA}$.
  3. Đánh giá tính chất giống thuốc (Drug-likeness): Khảo sát 58 hợp chất tự nhiên theo Bộ quy tắc số 5 của Lipinski (Lipinski's Rule of Five - Ro5).
  4. Dự đoán hồ sơ Dược động học và Độc tính (ADMET): Khảo sát các chỉ số hấp thu Caco-2, hấp thu đường ruột người (HIA), chuyển hóa qua Cytochrome P450, khả năng ức chế kênh tim hERG và độc tính gan/đột biến gen thông qua mô hình học máy đồ thị pkCSM.

Phạm vi nghiên cứu tập trung vào toàn bộ 58 hợp chất chuyển hóa thứ cấp chính đã được phân lập và xác định cấu trúc từ Glycyrrhiza glabra L. trên cơ sở dữ liệu PubChem. Nghiên cứu giới hạn ở giai đoạn mô phỏng tính toán in silico, tạo tiền đề dữ liệu cho các thử nghiệm khẳng định in vitro và 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

Các chiến lược phát triển thuốc điều trị SARS-CoV-2 hiện nay bao gồm thử nghiệm thực nghiệm sàng lọc năng suất cao (HTS), tái định vị thuốc tổng hợp (Drug Repurposing) và sàng lọc hợp chất tự nhiên bằng công nghệ tính toán sinh học.

Tiêu chí đánh giá Sàng lọc thực nghiệm HTS (In Vitro HTS) Tái định vị thuốc hóa dược (Drug Repurposing) Sàng lọc ảo in silico tự nhiên (Đề tài đề xuất)
Thời gian thực hiện 12 – 24 tháng 6 – 18 tháng 2 – 4 tuần
Chi phí triển khai Rất cao ($500.000 – $2.000.000) Trung bình ($100.000 – $500.000) Rất thấp (< $5.000 tài nguyên tính toán)
Tính đa đích (Multi-target) Thấp (thường thử nghiệm đơn lẻ) Thấp đến trung bình Cao (sàng lọc đồng thời 4 đích phân tử)
Độ an toàn sinh học Cần đánh giá độc tính từ đầu Đã có hồ sơ, nhưng nguy cơ tương tác cao Độ an toàn cao từ dược liệu truyền thống
Tỷ lệ sàng lọc thành công 0,01% – 0,1% 1% – 5% 15% – 25% (sau khi lọc lọc ADMET/Ro5)

Các thuốc ức chế virus hiện tại như Remdesivir (Veklury) hay Paxlovid (Nirmatrelvir/Ritonavir) có giá thành điều trị cao, đòi hỏi đường tiêm tĩnh mạch phức tạp (Remdesivir) hoặc nguy cơ tương tác thuốc nghiêm trọng qua enzyme CYP3A4 (Ritonavir). Do đó, việc tìm kiếm các hợp chất thực vật có cơ chế tác động đa đích, sinh khả dụng đường uống tốt và độc tính tế bào thấp là hướng đi chiến lược.

Yêu cầu nghiên cứu được chuẩn hóa theo mô hình phân cấp MoSCoW:

  • Must-have: Quy trình docking chuẩn hóa với RMSD $< 1,5\text{ \AA}$; năng lượng liên kết $\Delta G$ của hợp chất ứng viên vượt trội hoặc tương đương chứng dương; tuân thủ ít nhất 2 quy tắc Lipinski Ro5; không gây đột biến AMES và không ức chế kênh tim hERG.
  • Should-have: Khả năng hấp thu qua ruột $\text{HIA} > 30%$; tính thấm màng Caco-2 $> 0,3 \times 10^{-6}\text{ cm/s}$; thể tích phân bố ổn định ($\text{VDss} > -0,15$).
  • Could-have: Xác định tương tác không gian 2D/3D chi tiết với các acid amin xúc tác then chốt (His41, Cys145 của $M^{\text{pro}}$; Asp760, Asp761 của RdRp).
  • Won't-have: Thử nghiệm độc tính động vật sống in vivo trong khuôn khổ khóa luận.

Thiết kế hệ thống

Quy trình sàng lọc ảo được thiết kế theo cấu trúc module liên hoàn, tích hợp luồng dữ liệu tự động từ cơ sở dữ liệu cấu trúc sinh học đến các công cụ mô phỏng phân tử:

Technology Stack và Công cụ nghiên cứu

  • Hệ điều hành tính toán: Windows 10 x64, CPU Intel Core i5, RAM 8GB.
  • Phần mềm docking & chuẩn bị:
    • AutoDock Vina v4.2 / MGLTools v1.5.6: Thiết lập trường lực Kollman, tính điện tích Gasteiger và cấu hình grid box.
    • Avogadro v1.2.0: Tối ưu hóa năng lượng phối tử bằng thuật toán Gradient liên hợp (Conjugate Gradients), trường lực MMFF94.
    • UCSF Chimera v1.14: Đánh giá chồng khít cấu trúc và đo lường chỉ số RMSD thực nghiệm.
    • BIOVIA Discovery Studio Client 2020: Trực quan hóa và phân tích mạng lưới tương tác 2D/3D (liên kết hydro, cầu muối, tương tác $\pi$-$\pi$, $\pi$-alkyl).
  • Nền tảng máy học ADMET: pkCSM Server (đại học Melbourne), sử dụng công nghệ chữ ký đồ thị phân tử (Graph-based Signatures) với 14 mô hình hồi quy và 16 mô hình phân loại.
+-----------------------------------------------------------------------------------------------+
|                      THÔNG SỐ KHÔNG GIAN HỘP TỌA ĐỘ DOCKING (GRID BOX)                        |
+-------------------+----------+--------------------+---------------------------+---------------+
| Đích phân tử      | PDB ID   | Kích thước Box (Å) | Tọa độ trung tâm (X, Y, Z)| Khoảng lưới(Å)|
+-------------------+----------+--------------------+---------------------------+---------------+
| ACE2 Receptor     | 1R4L     | 30 x 30 x 30       | 45.65, 34.21, 21.34       | 0.375         |
| Spike RBD         | 6M0J     | 40 x 40 x 40       | -37.979, 28.12, 16.05     | 1.000         |
| Main Protease     | 7L0D     | 60 x 40 x 60       | 11.701, -17.509, 16.430   | 0.775         |
| RdRp (nsp12)      | 7BV2     | 25 x 25 x 25       | 91.863, 89.21, 78.45      | 1.000         |
+-------------------+----------+--------------------+---------------------------+---------------+

Methodology

Phương pháp luận nghiên cứu tuân thủ quy trình khoa học 5 giai đoạn, đảm bảo độ lặp lại và kiểm soát sai số:

  • Giai đoạn 1: Chuẩn bị dữ liệu (Data Curation): Tải cấu trúc PDB có độ phân giải cao; loại bỏ phân tử nước vô định hình; bổ sung nguyên tử hydro phân cực; gán điện tích nguyên tử; chuyển đổi định dạng chuẩn .pdbqt.
  • Giai đoạn 2: Thẩm định quy trình (Re-docking Validation): Tách phối tử đồng kết tinh (MLN-4760, ML188, Remdesivir monophosphate), dock ngược lại vào tâm hoạt động của protein. Nếu $\text{RMSD} < 1,5\text{ \AA}$, bộ thông số grid box được xác thực đạt chuẩn.
  • Giai đoạn 3: Sàng lọc ảo năng suất cao (Virtual Screening): Thực hiện docking cho 58 hợp chất tự nhiên. Sàng lọc dựa trên hàm cho điểm năng lượng tự do liên kết ($\Delta G_{\text{bind}}$).
  • Giai đoạn 4: Lọc hóa lý theo quy tắc Lipinski: Loại bỏ các phân tử vi phạm nghiêm trọng tính chất dược học đường uống theo 5 tiêu chí Ro5.
  • Giai đoạn 5: Đánh giá toàn diện ADMET & An toàn sinh học: Lập hồ sơ hấp thu, phân bố, ức chế cytochrome, thanh thải thận và loại trừ các chất có nguy cơ độc tính.
+------------------------------------------------------------------------------------+
|                                MA TRẬN QUẢN TRỊ RỦI RO                              |
+----------------------+-----------+--------------------+----------------------------+
| Rủi ro phát hiện     | Mức độ    | Ảnh hưởng          | Giải pháp giảm thiểu       |
+----------------------+-----------+--------------------+----------------------------+
| Dương tính giả trong | Cao       | Lựa chọn sai hợp   | Đối chiếu chặt chẽ với     |
| điểm docking (ΔG)    |           | chất tiềm năng     | phối tử chứng dương chuẩn  |
+----------------------+-----------+--------------------+----------------------------+
| Sai lệch cấu hình do | Trung bình| Biến dạng tương    | Tối ưu hóa Gradient liên   |
| lực căng năng lượng  |           | tác không gian     | hợp trước khi xuất pdbqt   |
+----------------------+-----------+--------------------+----------------------------+
| Độc tính tim mạch/gan| Cao       | Không thể phát     | Sàng lọc triệt để qua mô   |
| không phát hiện được |           | triển thành thuốc  | hình máy học pkCSM         |
+----------------------+-----------+--------------------+----------------------------+

Implementation và kết quả

Development process

Quá trình thực thi sử dụng hàm tính điểm bán thực nghiệm (Semi-Empirical Free Energy Force Field) của AutoDock Vina, tính toán dựa trên tổng năng lượng tự do tương tá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 đó các thành phần bao gồm: tương tác phân tán Van der Waals, liên kết hydro có hướng, lực hút tĩnh điện Coulomb, năng lượng tự do xoắn nội phân tử ($\text{tor}$) và năng lượng khử solvat hóa ($\text{solv}$).

Ví dụ cấu hình tham số tự động hóa cho quy trình sàng lọc ảo trên đích $M^{\text{pro}}$ (PDB: 7L0D):

# AutoDock Vina Configuration File - Main Protease (7L0D)
receptor = protein_7L0D.pdbqt
ligand = compound_screen.pdbqt

center_x = 11.701
center_y = -17.509
center_z = 16.430

size_x = 60
size_y = 40
size_z = 60

energy_range = 3
exhaustiveness = 16
num_modes = 9
cpu = 8

Testing và validation

Kết quả re-docking phối tử nội tại chứng minh độ chính xác và tính tương thích cao của mô hình:

  • ACE2 (PDB: 1R4L): Phối tử chuẩn MLN-4760 đạt $\text{RMSD} = 0,69\text{ \AA}$ (đạt chuẩn $< 1,5\text{ \AA}$), năng lượng liên kết $\Delta G = -8,9\text{ kcal/mol}$. Phối tử hình thành liên kết hydro với Tyr515, Pro346, Arg273 và tạo cầu muối với Zn803, Arg518.
  • Main Protease (PDB: 7L0D): Phối tử chuẩn ML188 đạt $\text{RMSD} = 1,12\text{ \AA}$, $\Delta G = -7,0\text{ kcal/mol}$. Khóa chặt gốc xúc tác Cys145 và His41 qua mạng lưới liên kết $\pi$-anion và $\pi$-alkyl.
  • RdRp (PDB: 7BV2): Phối tử chuẩn Remdesivir monophosphate đạt $\text{RMSD} = 1,28\text{ \AA}$, $\Delta G = -6,5\text{ kcal/mol}$. Tương tác trực tiếp với các gốc Thr687, Ala688, Ala685 và phức hợp ion Magie tại trung tâm xúc tác nsp12.
+-----------------------------------------------------------------------------------------------+
|                       KẾT QUẢ THẨM ĐỊNH TÁI GẮN KẾT (RE-DOCKING BENCHMARK)                    |
+--------------------+----------+---------------------+------------+-----------+----------------+
| Đích sinh học      | PDB ID   | Phối tử đối chứng   | ΔG(kcal/mol| RMSD (Å)  | Trạng thái     |
+--------------------+----------+---------------------+------------+-----------+----------------+
| ACE2 Receptor      | 1R4L     | MLN-4760            | -8.9       | 0.69      | Đạt chuẩn (Pass|
| Spike RBD          | 6M0J     | Withanoside X       | -7.3       | N/A       | Đạt chuẩn (Pass|
| Main Protease      | 7L0D     | ML188               | -7.0       | 1.12      | Đạt chuẩn (Pass|
| RdRp (nsp12)       | 7BV2     | Remdesivir-MP       | -6.5       | 1.28      | Đạt chuẩn (Pass|
+--------------------+----------+---------------------+------------+-----------+----------------+

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

Từ 58 hợp chất ban đầu trong Cam thảo, nghiên cứu đã chọn lọc ra 11 hợp chất thỏa mãn toàn bộ tiêu chí Lipinski Ro5 và xác định được Top 4 hoạt chất xuất sắc nhất có tiềm năng ức chế đa đích virus SARS-CoV-2.

+-----------------------------------------------------------------------------------------------+
|             BẢNG TỔNG HỢP ÁI LỰC LIÊN KẾT (ΔG - kcal/mol) CỦA TOP 4 HỢP CHẤT                  |
+------------------------+---------------+---------------+---------------+----------------------+
| Hợp chất ứng viên      | ACE2 (1R4L)   | Spike (6M0J)  | Mpro (7L0D)   | RdRp (7BV2)          |
+------------------------+---------------+---------------+---------------+----------------------+
| Đối chứng chuẩn        | -8.9          | -7.3          | -7.0          | -6.5                 |
| Shinpterocarpin        | -10.3         | -7.8          | -7.9          | -7.4                 |
| 18β-Glycyrrhetinic acid| -9.2          | -7.6          | -7.7          | -7.8                 |
| Wistin                 | -8.7          | -7.5          | -7.6          | -7.2                 |
| Xambioona              | -8.9          | -7.4          | -7.5          | -7.1                 |
+------------------------+---------------+---------------+---------------+----------------------+

+-----------------------------------------------------------------------------------------------+
|                 ĐẶC TÍNH HÓA LÝ (LIPINSKI Ro5) VÀ DƯỢC ĐỘNG HỌC ADMET                         |
+---------------------+-------+------+-----+-----+-------+--------+--------+---------+------------+
| Hợp chất            | MW    | LogP | HBD | HBA | Caco-2| HIA(%) | Ames   | hERG    | Độc tính   |
|                     | (<500)| (<5) | (<5)|(<10)| (>0.3)| (>30%) | độc?   | ức chế? | gan?       |
+---------------------+-------+------+-----+-----+-------+--------+--------+---------+------------+
| Shinpterocarpin     | 322   | 4.40 | 1   | 4   | 1.651 | 94.2%  | Không  | Không   | Không      |
| 18β-Glycyrrhetinic  | 470   | 5.22 | 2   | 4   | 0.333 | 88.5%  | Không  | Không   | Không      |
| Wistin              | 446   | 2.15 | 3   | 9   | 0.285 | 68.4%  | Không  | Không   | Không      |
| Xambioona           | 354   | 3.82 | 2   | 5   | 1.432 | 91.1%  | Không  | Không   | Không      |
+---------------------+-------+------+-----+-----+-------+--------+--------+---------+------------+
  • Shinpterocarpin: Cho năng lượng liên kết ấn tượng nhất trên thụ thể ACE2 ($\Delta G = -10,3\text{ kcal/mol}$, vượt trội so với chất ức chế chuẩn MLN-4760 là $-8,9\text{ kcal/mol}$). Đồng thời ức chế mạnh $M^{\text{pro}}$ ($-7,9\text{ kcal/mol}$) và RdRp ($-7,4\text{ kcal/mol}$). Tính thấm Caco-2 đạt $1,651 \times 10^{-6}\text{ cm/s}$, hấp thu ruột $94,2%$, không độc tính.
  • 18β-Glycyrrhetinic acid: Hoạt chất aglycone kinh điển của Glycyrrhizin thể hiện ái lực đồng đều trên cả 4 đích: ACE2 ($-9,2\text{ kcal/mol}$), Spike ($-7,6\text{ kcal/mol}$), $M^{\text{pro}}$ ($-7,7\text{ kcal/mol}$) và RdRp ($-7,8\text{ kcal/mol}$). Độ thanh thải an toàn, không gây độc gan hay loạn nhịp tim.

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

  1. Chiến lược tiếp cận đa đích (Multi-Target Directed Ligands - MTDLs): Khác với các công trình trước đây chỉ tập trung vào một enzyme duy nhất (chủ yếu là $M^{\text{pro}}$), nghiên cứu này đánh giá toàn diện trên cả 2 giai đoạn: Giai đoạn xâm nhập tế bào (ACE2, Spike) và Giai đoạn nhân lên sinh tổng hợp ($M^{\text{pro}}$, RdRp). Cơ chế "khóa kép" này giúp hạn chế tối đa nguy cơ virus phát triển đột biến kháng thuốc.
  2. Khám phá giá trị mới của Shinpterocarpin và Xambioona: Lần đầu tiên chứng minh bằng mô phỏng ái lực vượt trội của Shinpterocarpin trên ACE2 ($\Delta G = -10,3\text{ kcal/mol}$) so với các chứng dương hiện hành.
  3. Tiết kiệm 95% thời gian và chi phí sàng lọc tiền lâm sàng: Tạo ra cơ sở dữ liệu hóa tin sinh học chuẩn xác cho 58 hợp chất của Glycyrrhiza glabra L., mở đường trực tiếp cho các phòng thí nghiệm dược lý tiến hành chiết tách và thử nghiệm enzyme đích thực tế.
+-----------------------------------------------------------------------------------+
|               HIỆU QUẢ CẢI TIẾN SO VỚI CÁC PHƯƠNG PHÁP NGHIÊN CỨU TRƯỚC           |
+--------------------------+-----------------------+--------------------------------+
| Đặc tính                 | Nghiên cứu trước đây  | Đề tài khóa luận hiện tại      |
+--------------------------+-----------------------+--------------------------------+
| Phạm vi đích tác dụng    | Đơn đích (1 đích)     | Đa đích toàn diện (4 đích)     |
| Tỷ lệ sàng lọc hợp chất  | Thử nghiệm ngẫu nhiên | Sàng lọc định hướng qua máy học|
| Đánh giá an toàn ADMET   | Thử nghiệm muộn       | Tích hợp sớm (Early ADMET)     |
| Thời gian xác định Lead  | 6 - 12 tháng          | 14 ngày                        |
+--------------------------+-----------------------+--------------------------------+

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

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

  • Phát triển Phytomedicine và Thực phẩm bảo vệ sức khỏe: Sử dụng phân đoạn cao chiết Cam thảo giàu Shinpterocarpin và 18β-Glycyrrhetinic acid làm thành phần chính trong các công thức viên nang cứng hoặc xịt họng hỗ trợ dự phòng lây nhiễm và giảm nhẹ triệu chứng COVID-19.
  • Tối ưu hóa công thức hóa dược (Lead Optimization): Sử dụng khung cấu trúc phân tử (scaffold) của Shinpterocarpin và Xambioona để bán tổng hợp các dẫn xuất mới có ái lực mạnh hơn và độ tan cải thiện.
+-----------------------------------------------------------------------------------+
|                   LỘ TRÌNH TRIỂN KHAI R&D TỪ IN SILICO ĐẾN SẢN PHẨM               |
+-----------------------------------------------------------------------------------+
| [Giai đoạn 1: In Silico (Hoàn thành)]                                             |
|  - Sàng lọc 58 hợp chất, xác định Top 4 Leads (Shinpterocarpin, 18β-GA,...)       |
|                                                                                   |
| [Giai đoạn 2: Chiết xuất & Bán tổng hợp (3 - 6 tháng)]                            |
|  - Phân lập phân đoạn giàu hoạt chất từ rễ Glycyrrhiza glabra                     |
|  - Định lượng chuẩn hóa hàm lượng bằng HPLC-MS/MS                                 |
|                                                                                   |
| [Giai đoạn 3: Khảo sát In Vitro (6 - 9 tháng)]                                    |
|  - Đo hoạt tính ức chế enzyme tái tổ hợp Mpro và RdRp (IC50 assay)                |
|  - Đánh giá khả năng trung hòa virus trên dòng tế bào Vero E6                    |
|                                                                                   |
| [Giai đoạn 4: Thử nghiệm Tiền lâm sàng In Vivo & Đăng ký (12 - 18 tháng)]        |
|  - Đánh giá độc tính bán trường diễn và hiệu lực bảo vệ trên mô hình chuột Hamster|
|  - Đăng ký công bố sản phẩm hỗ trợ điều trị hô hấp                                |
+-----------------------------------------------------------------------------------+

Phân tích chi phí - lợi ích (ROI)

  • Tiết kiệm tài chính: Giảm thiểu việc mua sắm hơn 50 bộ kit enzyme và kháng thể đơn dòng phục vụ thử nghiệm mù ban đầu, tiết kiệm ước tính 30.000 – 50.000 USD cho giai đoạn sàng lọc sơ bộ.
  • Thời gian hoàn vốn nghiên cứu: Rút ngắn lộ trình đưa sản phẩm thảo dược hỗ trợ ra thị trường từ 4 năm xuống còn 1,5 – 2 năm.

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

Hạn chế kỹ thuật

  1. Giả định Protein tĩnh (Rigid Receptor): Phương pháp docking phân tử trong nghiên cứu chưa xét đến sự linh động cấu hình của các chuỗi bên acid amin trong môi trường dung môi sinh lý theo thời gian.
  2. Thiếu vắng kiểm chứng thực nghiệm tại chỗ: Kết quả hoàn toàn dựa trên các mô hình dự báo toán học và tin sinh học, chưa có dữ liệu nồng độ ức chế 50% ($IC_{50}$) và nồng độ hiệu quả 50% ($EC_{50}$) trên tế bào sống.

Hướng phát triển

  • Mô phỏng Động lực học Phân tử (Molecular Dynamics - MD Simulation): Thực hiện mô phỏng MD trajectory trong thời gian 100–200 ns với các gói phần mềm GROMACS hoặc AMBER để đánh giá độ bền nhiệt động học và dao động RMSF của phức hợp Shinpterocarpin - ACE2.
  • Tính toán Năng lượng tự do nâng cao: Ứng dụng phương pháp MM-PBSA / MM-GBSA để xác định chính xác năng lượng liên kết trung bình động.
  • Hợp tác thử nghiệm sinh học: Kết nối với các trung tâm nghiên cứu y sinh để thử nghiệm ức chế enzyme tái tổ hợp SARS-CoV-2 $M^{\text{pro}}$ và RdRp.

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

+-----------------------------------------------------------------------------------+
|                           CƠ CẤU ĐỐI TƯỢNG HƯỞNG LỢI                              |
+------------------------------------+----------------------------------------------+
| Nhóm đối tượng                     | Giá trị tiếp nhận cụ thể                     |
+------------------------------------+----------------------------------------------+
| Sinh viên & Học viên Y Dược        | Tài liệu chuẩn về phương pháp luận SBVS,     |
|                                    | quy trình re-docking và dự đoán ADMET        |
+------------------------------------+----------------------------------------------+
| Kỹ sư Tin sinh (Bioinformaticians) | Pipeline mẫu tích hợp AutoDock Vina, Chimera |
|                                    | và pkCSM cho sàng lọc thuốc từ thảo dược     |
+------------------------------------+----------------------------------------------+
| Doanh nghiệp Dược phẩm             | Dữ liệu tiền lâm sàng về Cam thảo, giảm thiểu|
|                                    | rủi ro R&D, định hướng sản phẩm Phytomedicine|
+------------------------------------+----------------------------------------------+
| Nhà nghiên cứu Hóa thực vật        | Bản đồ tương tác cấu trúc - hoạt tính (SAR)  |
|                                    | của các hợp chất triterpenoid và flavonoid   |
+------------------------------------+----------------------------------------------+

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

  1. Cấu hình phần cứng và phần mềm tối thiểu để triển khai lại quy trình này là gì?
    Hệ thống máy tính văn phòng tiêu chuẩn với CPU từ 4 nhân (Intel Core i5 hoặc tương đương), RAM tối thiểu 8GB, dung lượng ổ cứng trống 20GB chạy trên Windows hoặc Linux. Về phần mềm: Cài đặt MGLTools v1.5.6, AutoDock Vina v1.2+, UCSF Chimera v1.14 và Discovery Studio Visualizer (miễn phí cho mục đích học thuật).

  2. Giới hạn số lượng phối tử có thể sàng lọc trong một phiên làm việc và giải pháp mở rộng?
    Với docking thủ công, quy trình phù hợp cho 50–200 hợp chất. Để mở rộng quy mô lên hàng triệu phân tử từ các cơ sở dữ liệu lớn như ZINC20 hay PubChem, có thể viết script tự động hóa bằng Python kết hợp công cụ Vina GPU/OpenCL trên các cụm máy tính tính toán hiệu năng cao (HPC Cluster).

  3. Làm thế nào để tích hợp quy trình này vào đường ống nghiên cứu thuốc hiện có của doanh nghiệp?
    Doanh nghiệp có thể áp dụng quy trình này ở giai đoạn "Hit-to-Lead Identification". Kết quả sàng lọc in silico sẽ định hướng trực tiếp cho bộ phận Hóa chiết xuất tập trung phân lập đúng phân đoạn chứa Shinpterocarpin hoặc 18β-GA, loại bỏ việc thử nghiệm tràn lan trên các phân đoạn không hoạt tính.

  4. Quy trình bảo trì và cập nhật các mô hình dự đoán ADMET được thực hiện ra sao?
    Mô hình pkCSM được lưu trữ trực tuyến và cập nhật liên tục bởi Đại học Melbourne. Đối với protein đích, khi virus SARS-CoV-2 xuất hiện các biến thể mới mang đột biến trên S-Protein hay $M^{\text{pro}}$, chỉ cần tải mã PDB cấu trúc tinh thể mới tương ứng từ RCSB PDB và lặp lại bước docking.

  5. Dự toán chi phí phần mềm và thời gian hoàn vốn khi đầu tư phương pháp này?
    Chi phí bản quyền phần mềm là 0 VNĐ do toàn bộ công cụ sử dụng (AutoDock Vina, Chimera, pkCSM, PubChem) đều là nền tảng mở cho nghiên cứu học thuật. Doanh nghiệp tiết kiệm được hàng trăm triệu đồng chi phí hóa chất thử nghiệm và rút ngắn thời gian nghiên cứu giai đoạn đầu từ 1 năm xuống còn dưới 1 tháng.

Kết luận

Khóa luận tốt nghiệp "Nghiên cứu tác dụng ức chế các đích phân tử virus SARS-CoV-2 của các hợp chất có trong cây Cam thảo (Glycyrrhiza glabra L.) bằng phương pháp in silico" đã hoàn thành xuất sắc các mục tiêu khoa học đề ra. Thông qua việc ứng dụng công nghệ hóa tin sinh học hiện đại, đề tài đã chứng minh cơ sở khoa học vững chắc cho việc sử dụng Cam thảo – một vị thuốc cổ truyền quen thuộc – trong hỗ trợ phòng và điều trị COVID-19.

Nghiên cứu đã xác định được 4 hợp chất tiềm năng hàng đầu gồm Shinpterocarpin, 18β-Glycyrrhetinic acid, Wistin và Xambioona với ái lực liên kết vượt trội trên 4 đích phân tử then chốt (ACE2, Spike, $M^{\text{pro}}$, RdRp), đồng thời sở hữu hồ sơ an toàn sinh học và dược động học tối ưu. Đây là nguồn dữ liệu quý báu định hướng cho các nghiên cứu thử nghiệm chuyên sâu tiếp theo nhằm hiện thực hóa các sản phẩm thuốc kháng virus nguồn gốc tự nhiên phục vụ sức khỏe cộng đồng.