Bối cảnh và vấn đề nghiên cứu

Bệnh Parkinson (Parkinson’s Disease - PD) là rối loạn thoái hóa thần kinh phổ biến thứ hai trên thế giới sau bệnh Alzheimer, đặc trưng bởi sự mất dần các tế bào thần kinh dopaminergic ở liềm đen và sự xuất hiện bất thường của các thể Lewy chứa protein $\alpha$-synuclein. Theo thống kê dịch tễ được trích dẫn trong văn bản, hiện có hơn 6 triệu người trên toàn cầu đang sống chung với bệnh Parkinson; tỷ lệ mắc dao động từ 100 đến 200 trên 100.000 dân và tăng lên 1–3% ở nhóm dân số trên 65 tuổi. Bệnh khởi phát chủ yếu ở độ tuổi 65–70, nam giới có tỷ lệ mắc cao hơn nữ giới. Về mặt lâm sàng, bệnh tiến triển qua 5 giai đoạn từ tổn thương một bên đến mất hoàn toàn khả năng vận động tự chủ, gây ra cả triệu chứng vận động (run khi nghỉ, cứng cơ, vận động chậm, mất thăng bằng) lẫn triệu chứng không vận động (suy giảm khứu giác, trầm cảm, rối loạn giấc ngủ, sa sút trí tuệ).

Sau hơn hai thế kỷ kể từ khi được bác sĩ James Parkinson mô tả vào năm 1817, y học hiện đại vẫn chưa có biện pháp điều trị triệt để nguyên nhân. Liệu pháp bổ sung tiền chất dopamine bằng Levodopa (L-dopa) vẫn là lựa chọn hàng đầu nhưng việc sử dụng kéo dài thường gây các biến chứng vận động nghiêm trọng như loạn vận động, hạ huyết áp, buồn nôn, đồng thời không giải quyết được các triệu chứng phi vận động. Các nhóm thuốc bổ trợ hiện hành như chất ức chế Catechol-O-methyl-transferase (COMT) hay ức chế Monoamine oxidase B (MAO-B) chỉ làm chậm quá trình giáng hóa dopamine. Vì vậy, việc tìm kiếm các liệu pháp tác động đa đích, vừa cải thiện triệu chứng vừa bảo vệ tế bào thần kinh, là một yêu cầu cấp thiết.

Theo y học cổ truyền phương Đông, dược liệu Câu đằng (Uncaria rhynchophylla, họ Cà phê Rubiaceae) từ lâu đã được sử dụng để điều trị chứng kinh phong, co giật, run rẩy và các bệnh lý mạch máu não. Các bằng chứng thực nghiệm hiện đại xác nhận dịch chiết và các alkaloid từ Câu đằng thể hiện khả năng bảo vệ tế bào thần kinh chống độc tính glutamate, giảm stress oxy hóa và điều hòa tín hiệu nội bào. Đồng thời, các nghiên cứu sinh học phân tử đã xác định thụ thể $A_{2A}$ Adenosine ($A_{2A}AR$) và enzym Monoamine oxidase B (MAO-B) là hai đích phân tử then chốt trong bệnh sinh Parkinson. Chất đối kháng $A_{2A}AR$ giúp tăng cường dẫn truyền dopaminergic thể vân và bảo vệ tế bào thần kinh; trong khi ức chế MAO-B ngăn chặn sự phân hủy dopamine nội sinh và giảm sinh các gốc tự do có hại. Việc ứng dụng phương pháp in silico (docking phân tử kết hợp dự đoán ADMET) cho phép sàng lọc đồng thời trên hai đích này một cách nhanh chóng, tiết kiệm chi phí và định hướng chính xác cho các thử nghiệm tiền lâm sàng tiếp theo.

Đề tài khóa luận đặt ra 2 mục tiêu nghiên cứu cụ thể:

  1. Sàng lọc các hợp chất trong cây Câu đằng (Uncaria rhynchophylla) có tác dụng ức chế đồng thời hai đích phân tử $A_{2A}$ adenosine receptor ($A_{2A}AR$) và Monoamine oxidase B (MAO-B) bằng phương pháp docking phân tử.
  2. Nghiên cứu các đặc điểm "giống thuốc" (drug-likeness) theo quy tắc Lipinski và tính toán các thông số dược động học, độc tính (ADMET) của các hợp chất tiềm năng nhất thu được sau quá trình sàng lọc.
  • Đối tượng nghiên cứu: 80 hợp chất hóa học có nguồn gốc tự nhiên từ cây Câu đằng (Uncaria rhynchophylla).
  • Phạm vi nghiên cứu: Mô phỏng tính toán tin sinh học (in silico); sử dụng cấu trúc tinh thể ba chiều của hai protein đích gồm $A_{2A}AR$ (PDB ID: 3EML) và MAO-B (PDB ID: 2V5Z) lấy từ Ngân hàng Dữ liệu Protein (RCSB PDB); dữ liệu cấu trúc các phối tử lấy từ PubChem; thời gian và địa điểm hoàn thành khóa luận là tháng 6 năm 2022 tại Bộ môn Dược lý, Trường Đại học Y Dược, Đại học Quốc gia Hà Nội.

Cơ sở lý thuyết và phương pháp

Khóa luận được xây dựng dựa trên cơ sở sinh học phân tử của hai đích tác dụng thần kinh trung ương:

  • Thụ thể $A_{2A}$ Adenosine ($A_{2A}AR$): Là thụ thể bắt cặp protein G có cấu trúc gồm 7 chuỗi xoắn xuyên màng (mỗi chuỗi 25–27 acid amin) nối bởi 3 vòng nội bào (IL1–IL3) và 3 vòng ngoại bào (EL1–EL3). Đối kháng $A_{2A}AR$ giúp tạo tương tác bù trừ sinh lý với thụ thể dopamine tại thể vân, giảm các rối loạn vận động thứ phát do dùng L-dopa và bảo vệ tế bào thần kinh chống lại stress oxy hóa.
  • Enzym Monoamine oxidase B (MAO-B): Là enzym liên kết ở màng ngoài ty thể của tế bào thần kinh đệm và tiểu cầu, chiếm 80% hoạt tính MAO trong não. MAO-B xúc tác khử amin oxy hóa dopamine tạo ra hydrogen peroxide ($H_2O_2$) và chuyển hóa tiền độc chất MPTP thành ion $MPP^+$ gây thoái hóa tế bào dopaminergic. Ức chế MAO-B giúp bảo tồn nồng độ dopamine tại khe synap, ức chế dòng calci nội bào qua ty thể và hạn chế quá trình kết tập protein $\alpha$-synuclein.

Phương pháp nghiên cứu thực hiện hoàn toàn trên máy tính (in silico), sử dụng máy tính Dell Vostro 3510 với các công cụ chuyên dụng:

  • Mô phỏng Docking phân tử: Sử dụng phần mềm Autodock Vina với hàm tính điểm năng lượng tự do liên kết ($\Delta G$, kcal/mol). Quy trình chuẩn bị protein đích được thực hiện bằng Discovery Studio 2021 Client (loại nước, loại phối tử đồng kết tinh ZMA/SAG) và AutoDockTools 1.5.6 (thêm hydro phân cực, gán điện tích Kollman, xuất định dạng pdbqt). Cấu trúc 80 phối tử từ cơ sở dữ liệu PubChem được tối ưu hóa cấu hình không gian 3 chiều bằng thuật toán Gradient liên hợp (Conjugate Gradients) trong phần mềm Avogadro.
  • Thẩm định quy trình (Re-docking): Thực hiện tái gắn kết phối tử đồng kết tinh vào vùng hoạt động của protein mục tiêu. Tính toán độ lệch bình phương trung bình gốc (RMSD) bằng phần mềm UCSF Chimera. Tiêu chuẩn quy trình đạt độ tin cậy là giá trị RMSD < 1,5 Å. Phối tử đồng kết tinh ZMA (ức chế $A_{2A}AR$) và SAG (Safinamide, ức chế MAO-B) được sử dụng làm chất chứng dương đối chứng năng lượng liên kết.
  • Không gian hộp tìm kiếm (Grid box): Đích $A_{2A}AR$ (PDB ID: 3EML) có kích thước 70 × 70 × 80 Å, tọa độ tâm x = -5,968, y = -10,31, z = 53,267. Đích MAO-B (PDB ID: 2V5Z) có kích thước 70 × 70 × 70 Å, tọa độ tâm x = 51,661, y = 155,622, z = 27,829; khoảng cách mắt lưới (spacing) của cả hai hệ thống là 0,375 Å.
  • Đánh giá đặc điểm giống thuốc: Áp dụng bộ Quy tắc số 5 của Lipinski (Rule of 5 - RO5) thông qua công cụ trực tuyến của Scfbio. Hợp chất được xem là có tiềm năng hấp thu tốt qua đường uống khi vi phạm không quá 2 trong 5 tiêu chí: Trọng lượng phân tử (MW) < 500 Da; 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; Hệ số phân bố octanol/nước (LogP) < 5; Độ khúc xạ mol (MR) từ 40 đến 130.
  • Dự đoán ADMET: Sử dụng nền tảng ADMETlab với đầu vào là chuỗi ký hiệu phân tử SMILES để sàng lọc các đặc tính hấp thu (độ tan LogS, tính thấm tế bào Caco-2, hấp thu ruột người), phân bố (thể tích phân bố VD, tính thấm hàng rào máu não BBB, liên kết protein huyết tương PPB), chuyển hóa (khả năng ức chế 5 isoenzym Cytochrome P450: CYP1A2, CYP3A4, CYP2C9, CYP2C19, CYP2D6), thải trừ (thời gian bán thải $T_{1/2}$, độ thanh thải toàn phần CL) và độc tính (độc tính gan, độc tính gây đột biến AMES, kích ứng da).

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương này hệ thống hóa các tài liệu khoa học về bệnh sinh, chẩn đoán lâm sàng và phân loại 5 giai đoạn tiến triển của bệnh Parkinson theo thang điểm vận động. Tác giả phân tích ưu - nhược điểm của các phác đồ điều trị dùng thuốc (Levodopa, thuốc chủ vận dopamine, thuốc ức chế COMT, thuốc ức chế MAO-B) và không dùng thuốc (kích thích não sâu, phẫu thuật, phục hồi chức năng). Tiếp theo, chương trình bày chi tiết về cấu trúc sinh học phân tử, cơ chế tác dụng của hai đích thụ thể $A_{2A}AR$ và enzym MAO-B trong việc duy trì nồng độ dopamine và bảo vệ tế bào thần kinh. Phần dược liệu mô tả đặc điểm thực vật, phân bố vùng thượng du Việt Nam, bộ phận dùng (đoạn thân/cành có gai móc câu sấy khô), tính vị quy kinh theo Đông y và các nhóm hoạt chất sinh học (alkaloid oxindol như rhynchophylline, isorhynchophylline, corynoxeine) trong cây Câu đằng (Uncaria rhynchophylla). Cuối cùng, chương giới thiệu nguyên lý của phương pháp tính toán in silico, cơ chế hoạt động của thuật toán docking phân tử, quy chuẩn đánh giá RMSD, quy tắc Lipinski 5 đối với thuốc đường uống và mô hình phân tích dược động học/độc tính ADMET.

Chương 2: Đối tượng và phương pháp nghiên cứu

Chương 2 mô tả chi tiết nguồn dữ liệu và thiết lập các bước thực nghiệm tính toán:

  • Thu thập danh mục 80 hợp chất từ các công trình công bố về thành phần hóa thực vật của cây Câu đằng; tải cấu trúc không gian dạng SDF từ PubChem.
  • Lựa chọn cấu trúc protein độ phân giải cao từ RCSB PDB: thụ thể $A_{2A}AR$ người phức hợp với chất đối kháng ZMA (PDB ID: 3EML) và enzym MAO-B người gắn chất ức chế chọn lọc Safinamide SAG (PDB ID: 2V5Z).
  • Chuẩn hóa quy trình re-docking: tách phối tử đồng kết tinh, tối ưu hóa phối tử và protein, thiết lập tọa độ hộp lưới bao phủ trung tâm hoạt động và thẩm định giá trị RMSD.
  • Thiết lập quy trình 3 bước nghiên cứu: (1) Sàng lọc docking ảo 80 chất trên cả hai đích protein; (2) Lọc các hợp chất thỏa mãn quy tắc Lipinski 5; (3) Phân tích chuyên sâu 5 nhóm chỉ số ADMET bằng công cụ ADMETlab để chọn lọc các phân tử tiềm năng nhất.

Chương 3: Kết quả

Chương 3 trình bày các số liệu tính toán cụ thể thu được qua từng bước:

  • Thẩm định phương pháp: Tái gắn kết hai phối tử đồng kết tinh cho thấy cấu dạng trước và sau khi dock chồng khít rất tốt; giá trị RMSD của ZMA trên $A_{2A}AR$ là 1,24 Å và của SAG trên MAO-B là 1,382 Å (đều < 1,5 Å). Điểm số năng lượng tự do liên kết của ZMA là -7,7 kcal/mol và SAG là -9,8 kcal/mol, được cố định làm ngưỡng đối chứng dương.
  • Sàng lọc docking 80 hợp chất: Có 38/80 hợp chất cho năng lượng gắn kết trên $A_{2A}AR$ âm hơn chứng dương ZMA; 11/80 hợp chất trên MAO-B âm hơn hoặc bằng chứng dương SAG. Có 9 hợp chất thể hiện ái lực vượt trội trên đồng thời cả hai đích: 3$\alpha$-dihydrocadambine, angustidine, cinchonain Id, acid cryptochlorogenic, acid hyptatic A, saponin A, timosaponin A-III, acid uncarinic B và acid ursolic lactone. Trong đó, timosaponin A-III đạt điểm docking âm nhất trên $A_{2A}AR$ (-14,9 kcal/mol) và saponin A đạt điểm âm nhất trên MAO-B (-14,2 kcal/mol).
  • Sàng lọc Lipinski: Đánh giá 9 hợp chất cho thấy 6/9 chất đạt tiêu chuẩn "giống thuốc" (vi phạm $\le$ 2 tiêu chí), riêng angustidine thỏa mãn cả 5/5 tiêu chí. Ba hợp chất bị loại do khối lượng phân tử quá lớn và vi phạm nhiều thông số là 3$\alpha$-dihydrocadambine, saponin A và timosaponin A-III.
  • Phân tích ADMET và tương tác không gian: Đánh giá 6 hợp chất giống thuốc trên ADMETlab cho thấy angustidine có nguy cơ gây độc gan và độc tính AMES; acid cryptochlorogenic có tỷ lệ hấp thu qua ruột thấp (32,5%) và không qua được hàng rào máu não (BBB). Từ đó, tác giả lựa chọn được 4 hợp chất an toàn và có dược động học khả quan nhất gồm: cinchonain Id, acid hyptatic A, acid uncarinic B và acid ursolic lactone. Mô phỏng tương tác 3D chỉ ra 4 phân tử này chiếm trọn khoang gắn kết của hai đích protein, thiết lập các liên kết hydro, liên kết $\pi$-$\pi$, $\pi$-alkyl và tương tác van der Waals với các acid amin chủ chốt (Leu-167, Ala-63, Ile-274, Phe-168... trên $A_{2A}AR$ và Pro-102, Cys-172, Tyr-326, Phe-343, Tyr-398... trên MAO-B).

Chương 4: Bàn luận

Tác giả phân tích sâu về đặc điểm cấu trúc hóa học, mối quan hệ cấu trúc - hoạt tính và các bằng chứng thực nghiệm đã công bố của 4 hợp chất tiềm năng:

  • Cinchonain Id: Là một flavonoid phenylpropanoid-substituted epicatechin. Hợp chất chứa 6 nhóm hydroxyl liên kết với các vòng thơm, tạo điều kiện thuận lợi hình thành mạng lưới liên kết hydro dày đặc với túi gắn kết, cho năng lượng docking rất tốt (-10,2 kcal/mol trên $A_{2A}AR$ và -9,8 kcal/mol trên MAO-B). Mặc dù tỷ lệ hấp thu ruột ở mức trung bình (45,8%), chất này có khả năng qua hàng rào máu não, không ức chế các isoenzym CYP và không thể hiện độc tính gan hay đột biến gen.
  • Acid hyptatic A: Thuộc nhóm triterpene, có tính thấm Caco-2 tối ưu (-5,008 logPapp cm/s), hấp thu ruột đạt 57,9%, qua được hàng rào máu não, âm tính với các phép thử độc tính và không ức chế hệ enzym chuyển hóa CYP.
  • Acid uncarinic B: Là một triterpenoid tự nhiên có hoạt tính ức chế phospholipase C$\gamma$1 đã biết ($IC_{50}$ = 44,55 $\mu$M). Năng lượng docking trên hai đích lần lượt là -8,4 và -9,8 kcal/mol, tỷ lệ hấp thu ruột đạt 69,2%. Tuy nhiên, chất này bị dự đoán ức chế CYP3A4 và không thấm qua hàng rào máu não (BBB), đòi hỏi các nghiên cứu cải biến cấu trúc hoặc hệ dẫn thuốc để vượt qua màng não.
  • Acid ursolic lactone: Đạt năng lượng liên kết cao (-9,3 kcal/mol trên $A_{2A}AR$ và -9,9 kcal/mol trên MAO-B). Phân tử có tính thấm Caco-2 rất tốt (-4,745 logPapp cm/s), tỷ lệ hấp thu qua ruột người cao nhất trong nhóm triterpen (80,4%), qua được hàng rào máu não và không gây độc tế bào gan hay đột biến gen.

Kết quả và đóng góp

Nghiên cứu đã xử lý và phân tích có hệ thống 80 hợp chất từ cây Câu đằng trên mô hình docking hai đích thụ thể. Dữ liệu mô phỏng năng lượng liên kết và các đặc tính hóa lý được tổng hợp cụ thể trong các bảng dưới đây:

STT Tên hợp chất PubChem ID Năng lượng gắn kết $A_{2A}AR$ (3EML) (kcal/mol) Năng lượng gắn kết MAO-B (2V5Z) (kcal/mol)
- Chất chứng dương ZMA - -7,7 -
- Chất chứng dương SAG - - -9,8
1 3$\alpha$-Dihydrocadambine 162138 -12,6 -11,7
2 Angustidine 3084770 -10,6 -9,8
3 Cinchonain Id 21676382 -10,2 -9,8
4 Acid cryptochlorogenic 9798666 -8,3 -10,4
5 Acid hyptatic A 12305222 -7,9 -9,8
6 Saponin A 85144169 -11,2 -14,2
7 Timosaponin A-III 15953793 -14,9 -13,4
8 Acid uncarinic B 100967916 -8,4 -9,8
9 Acid ursolic lactone 21669105 -9,3 -9,9

Kết quả sàng lọc đặc điểm "giống thuốc" theo 5 tiêu chí của Lipinski trên 9 hợp chất tiềm năng:

STT Tên hợp chất Khối lượng phân tử (MW, Da) Nhóm cho H (HBD) Nhóm nhận H (HBA) LogP Độ khúc xạ mol (MR) Kết luận giống thuốc
1 3$\alpha$-Dihydrocadambine 546 6 11 -0,71 133,74 Không
2 Angustidine 301 1 3 3,15 88,95 Có
3 Cinchonain Id 452 6 9 2,69 112,82 Có
4 Acid cryptochlorogenic 345 6 9 -0,65 82,52 Có
5 Acid hyptatic A 488 4 5 5,18 133,48 Có
6 Saponin A 957 11 18 -0,12 231,00 Không
7 Timosaponin A-III 741 7 13 1,44 183,00 Không
8 Acid uncarinic B 649 3 6 8,19 181,66 Có
9 Acid ursolic lactone 456 1 3 6,76 130,33 Có

Kết quả phân tích 5 nhóm chỉ số dược động học và an toàn ADMET trên 6 hợp chất đạt quy tắc Lipinski:

Thông số đánh giá Angustidine Cinchonain Id Acid cryptochlorogenic Acid hyptatic A Acid uncarinic B Acid ursolic lactone
Độ tan LogS ($\log \mu\text{g/mL}$) 26,246 56,825 3849,330 0,673 2,774 0,126
Tính thấm Caco-2 ($\text{logPapp cm/s}$) -4,750 -6,682 -6,553 -5,008 -5,246 -4,745
Hấp thu ở ruột (%) 86,3 45,8 32,5 57,9 69,2 80,4
Thể tích phân bố VD ($\log\text{L/kg}$) 0,698 -0,800 -1,279 -0,574 -0,490 -0,343
Tính thấm màng não BBB +++ ++ --- + --- ++
Liên kết protein huyết tương (%) 81,347 85,633 50,151 88,088 85,122 82,175
Ức chế CYP1A2 ++ --- --- --- --- ---
Ức chế CYP3A4 - --- --- --- + ---
Thời gian bán thải $T_{1/2}$ (giờ) 2,035 1,631 0,473 1,860 2,317 1,639
Độ thanh thải toàn phần ($\log\text{ml/phút/kg}$) 2,164 1,138 1,195 1,394 1,307 1,462
Độc tính gan ++ - + --- - ---
Độc tính đột biến gen AMES + - --- --- --- ---
Kích ứng da --- --- - --- --- -

Đóng góp mới của tác giả so với các công bố trước:

  • Đây là nghiên cứu đầu tiên áp dụng phương pháp sàng lọc ảo in silico đồng thời trên cả hai đích phân tử $A_{2A}AR$ và MAO-B đối với thành phần hóa học của cây Câu đằng (Uncaria rhynchophylla).
  • Trong khi các nghiên cứu trước đây về Câu đằng chỉ tập trung vào nhóm hoạt chất alkaloid chính (rhynchophylline, isorhynchophylline), khóa luận đã phát hiện 4 hợp chất không thuộc nhóm alkaloid chính (cinchonain Id, acid hyptatic A, acid uncarinic B, acid ursolic lactone) có ái lực gắn kết đồng thời rất mạnh trên cả hai thụ thể đích, mở ra hướng tiếp cận hóa thực vật mới cho nhóm chất triterpene và flavonoid trong điều trị Parkinson.
  • Tác giả đề xuất 2 nhóm kiến nghị cụ thể: (1) Tiến hành thực nghiệm chiết xuất, phân lập các hợp chất này từ mẫu dược liệu Câu đằng thực tế để đánh giá hoạt tính ức chế enzym trên các mô hình in vitro và in vivo; (2) Nghiên cứu bào chế cải thiện độ tan, khả năng hấp thu và tính thấm qua hàng rào máu não đối với các phân tử triterpene, đồng thời mở rộng mô hình docking trên các đích phân tử khác liên quan đến cơ chế bệnh sinh thoái hóa thần kinh.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận được hoàn thành dựa trên phương pháp tính toán mô phỏng in silico, do đó tồn tại một số hạn chế khách quan được tác giả chỉ rõ:

  • Kết quả năng lượng gắn kết và thông số ADMET hoàn toàn là các dự đoán dựa trên thuật toán máy tính, chưa được đối chứng và thẩm định qua các thử nghiệm sinh học thực nghiệm (in vitro hoặc in vivo).
  • Một số phân tử có điểm docking rất cao nhưng lại gặp rào cản dược động học: acid uncarinic B bị cản trở bởi hàng rào máu não (BBB âm tính) và gây ức chế enzym CYP3A4; acid ursolic lactone có độ tan trong nước rất thấp (LogS = 0,126 $\log \mu\text{g/mL}$).
  • Các hướng nghiên cứu tiếp theo bao gồm:
    • Chiết xuất và thử nghiệm hoạt tính ức chế trực tiếp enzym MAO-B và thụ thể $A_{2A}AR$ trên dòng tế bào và mô hình động vật mắc hội chứng Parkinson gây bởi MPTP.
    • Tối ưu hóa cấu trúc hóa học hoặc ứng dụng công nghệ nano, chất mang lipid để khắc phục nhược điểm về độ tan và tăng khả năng vận chuyển qua hàng rào máu não của acid uncarinic B và acid ursolic lactone.
    • Mở rộng quy mô sàng lọc trên các đích thụ thể khác như thụ thể dopamine $D_2$, thụ thể NMDA và quá trình kết tập sợi tơ $\alpha$-synuclein.

Giá trị tham khảo

Công trình khóa luận có giá trị tham khảo chuyên môn thiết thực cho các đối tượng và định hướng nghiên cứu sau:

  • Sinh viên và học viên ngành Dược học, Hóa dược, Dược lý học: Cung cấp tài liệu tổng quan toàn diện về bệnh sinh thoái hóa thần kinh, cơ chế tác dụng của hệ thụ thể adenosine và enzym monoamine oxidase, cùng tổng quan hóa thực vật chi tiết của dược liệu Câu đằng (Uncaria rhynchophylla).
  • Người làm nghiên cứu về thiết kế thuốc tính toán (CADD/In silico): Cung cấp quy trình thực nghiệm chuẩn mực về sàng lọc ảo: từ bước chuẩn bị cấu trúc thụ thể - phối tử, xác thực mô hình bằng phương pháp re-docking tính chỉ số RMSD, cài đặt không gian hộp lưới (grid box), đến việc liên kết kết quả docking với bộ lọc Lipinski và nền tảng ADMETlab.
  • Nghiên cứu sinh phát triển thuốc từ nguồn gốc tự nhiên: Dữ liệu tương tác phân tử 3D và các giá trị năng lượng gắn kết chi tiết của 4 hợp chất (cinchonain Id, acid hyptatic A, acid uncarinic B, acid ursolic lactone) là cơ sở dữ liệu ban đầu đáng tin cậy để định hướng phân lập hoạt chất có mục tiêu từ thảo dược.

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

1. Tiêu chuẩn nào được sử dụng để chứng minh quy trình docking trong khóa luận là đáng tin cậy?
Tác giả đã tiến hành re-docking hai phối tử đồng kết tinh ZMA và SAG vào vị trí hoạt động của cấu trúc tinh thể $A_{2A}AR$ (PDB ID: 3EML) và MAO-B (PDB ID: 2V5Z). Giá trị sai lệch RMSD đo được lần lượt là 1,24 Å và 1,382 Å. Cả hai giá trị này đều nhỏ hơn ngưỡng quy chuẩn 1,5 Å, chứng minh quy trình mô phỏng đạt độ chính xác cao.

2. Hợp chất nào đạt năng lượng liên kết âm nhất trên từng đích protein và tại sao chúng không được chọn làm chất tiềm năng phát triển thuốc uống?
Hợp chất đạt năng lượng gắn kết âm nhất trên $A_{2A}AR$ là timosaponin A-III (-14,9 kcal/mol) và trên MAO-B là saponin A (-14,2 kcal/mol). Tuy nhiên, cả hai chất này đều bị loại vì vi phạm nghiêm trọng quy tắc 5 Lipinski: khối lượng phân tử quá lớn (lần lượt là 741 Da và 957 Da), số nhóm cho/nhận liên kết hydro quá cao (saponin A có 11 HBD và 18 HBA; timosaponin A-III có 7 HBD và 13 HBA), không phù hợp để hấp thu qua đường uống.

3. Vì sao hợp chất angustidine thỏa mãn cả 5 tiêu chí của Lipinski nhưng lại không nằm trong nhóm 4 chất được khuyến nghị cuối cùng?
Mặc dù angustidine có thông số "giống thuốc" tối ưu và độ thấm qua hàng rào máu não rất tốt (BBB +++), công cụ ADMETlab đã dự đoán hợp chất này có nguy cơ gây độc tính trên tế bào gan (mức độ ++) và cho kết quả dương tính với thử nghiệm độc tính đột biến gen AMES (+), đồng thời gây ức chế isoenzym CYP1A2.

4. Bốn hợp chất tiềm năng nhất được lựa chọn trong khóa luận là những chất nào?
Bốn hợp chất được chọn lọc gồm: cinchonain Id (flavonoid), acid hyptatic A, acid uncarinic B và acid ursolic lactone (đều thuộc nhóm triterpen/triterpenoid). Cả 4 chất đều có năng lượng liên kết trên hai đích tương đương hoặc âm hơn chứng dương, thỏa mãn quy tắc Lipinski và không thể hiện độc tính gan, kích ứng da hay đột biến gen AMES.

5. Rào cản dược động học chính của acid uncarinic B khi phát triển thành thuốc điều trị bệnh thần kinh trung ương là gì?
Dự đoán ADMET cho thấy acid uncarinic B không có khả năng đi qua hàng rào máu não (kết quả BBB ở mức âm tính "---") và có khả năng ức chế isoenzym CYP3A4 (+), do đó cần phải nghiên cứu các giải pháp bào chế dẫn thuốc chuyên biệt hoặc biến đổi cấu trúc để thuốc có thể thâm nhập vào mô não.


Kết luận

Khóa luận tốt nghiệp của sinh viên Cao Thị Lan Anh đã xây dựng thành công quy trình sàng lọc ảo in silico kết hợp mô phỏng docking phân tử, đánh giá quy tắc Lipinski và phân tích ADMET cho 80 hợp chất từ cây Câu đằng (Uncaria rhynchophylla). Nghiên cứu đã chứng minh tính tin cậy của mô hình tính toán qua phép thử re-docking với giá trị RMSD dưới 1,5 Å và xác định được 9 hợp chất có ái lực gắn kết đồng thời vượt trội trên cả hai đích phân tử then chốt của bệnh Parkinson là $A_{2A}AR$ và MAO-B. Từ 6 hợp chất thỏa mãn tiêu chuẩn "giống thuốc" đường uống, tác giả đã chọn lọc được 4 phân tử tiềm năng nhất gồm cinchonain Id, acid hyptatic A, acid uncarinic B và acid ursolic lactone với đặc tính dược động học khả quan và mức độ an toàn cao. Đây là luận cứ khoa học quan trọng làm tiền đề cho các định hướng chiết xuất phân lập và thử nghiệm tiền lâm sàng tiếp theo nhằm phát triển các thuốc mới điều trị thoái hóa thần kinh từ nguồn dược liệu Việt Nam.