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

Bệnh Parkinson (Parkinson’s Disease - PD) là bệnh lý rối loạn thoái hóa thần kinh mạn tính tiến triển phổ biến thứ hai trên thế giới, chỉ đứng sau bệnh Alzheimer. Bệnh sinh của Parkinson đặc trưng bởi sự mất mát chọn lọc các tế bào thần kinh giải phóng dopamine tại vùng chất đen thuộc thân não, đi kèm với sự tích tụ bất thường của các thể Lewy chứa protein $\alpha$-synuclein. Theo các số liệu dịch tễ học được tác giả trích dẫn, toàn cầu hiện có hơn 6 triệu người đang sống chung với căn bệnh này. Tỷ lệ mắc bệnh dao động trong khoảng 100 đến 200 ca trên 100.000 dân, với tỷ lệ mắc mới hàng năm là 15 trên 100.000 dân; độ tuổi khởi phát trung bình từ 65 đến 70 tuổi (chiếm 1-3% dân số trên 65 tuổi) và nam giới có tỷ lệ mắc cao hơn nữ giới. Bệnh tiến triển qua 5 giai đoạn suy giảm chức năng, từ mức độ nhẹ xuất hiện triệu chứng một bên (giai đoạn 1) đến giai đoạn tàn phế hoàn toàn phải phụ thuộc vào xe lăn hoặc nằm liệt giường (giai đoạn 5).

Biểu hiện lâm sàng của bệnh bao gồm hai nhóm triệu chứng chính:

  • Nhóm triệu chứng vận động (Motor Symptoms - MS): run khi nghỉ, co cứng cơ, vận động chậm chạp, mất thăng bằng và rối loạn tư thế.
  • Nhóm triệu chứng phi vận động (Non-Motor Symptoms - NMS): suy giảm khứu giác, rối loạn giấc ngủ, trầm cảm, lo âu, táo bón và suy giảm nhận thức hoặc sa sút trí tuệ.

Về mặt trị liệu, sau hơn 200 năm kể từ khi được bác sĩ James Parkinson mô tả lần đầu vào năm 1817, y học hiện đại vẫn chưa tìm ra phương pháp chữa khỏi hoàn toàn bệnh Parkinson. Phác đồ điều trị nội khoa hiện tại chủ yếu tập trung vào việc bù đắp lượng dopamine thiếu hụt bằng tiền chất Levodopa (L-dopa) hoặc các thuốc chủ vận dopamine, thuốc ức chế enzym phân hủy ngoại vi và trung ương như chất ức chế catechol-O-methyl-transferase (COMT) và monoamine oxidase B (MAO-B). Mặc dù L-dopa mang lại hiệu quả kiểm soát triệu chứng vận động rõ rệt trong giai đoạn đầu, việc sử dụng kéo dài thường dẫn đến các biến chứng vận động nghiêm trọng (như hiện tượng dao động vận động, loạ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 và không làm chậm được quá trình thoái hóa nơ-ron. Do đó, việc tìm kiếm các liệu pháp trị liệu mới tác động lên các đích sinh học phi dopaminergic hoặc đa đích là yêu cầu cấp thiết.

Theo Y học cổ truyền, Câu đằng (Uncaria rhynchophylla), một loài thảo dược thuộc họ Cà phê (Rubiaceae) còn được gọi là "Gou-Teng", từ lâu đã được sử dụng để điều trị các chứng kinh phong, co giật, run rẩy và các bệnh lý liên quan đến hệ thần kinh trung ương. Nghiên cứu thực nghiệm lâm sàng tại Trung Quốc cho thấy việc sử dụng dược liệu này giúp cải thiện đáng kể các rối loạn vận động và hỗ trợ các triệu chứng phi vận động như lo lắng, mất ngủ, táo bón ở bệnh nhân Parkinson. Các nghiên cứu hiện đại đã chứng minh các hoạt chất sinh học (đặc biệt là nhóm alkaloid như rhynchophylline, isorhynchophylline, corynoxeine, isocorynoxeine, hirsutine) có tác dụng bảo vệ tế bào thần kinh, chống thiếu máu cục bộ, ngăn chặn giải phóng canxi nội bào và ức chế các thụ thể thần kinh liên quan.

Trên cơ sở đó, việc kết hợp giữa kinh nghiệm y học cổ truyền và công nghệ thiết kế thuốc hỗ trợ bởi máy tính (in silico) thông qua phương pháp mô phỏng gắn kết phân tử (docking phân tử) mở ra hướng sàng lọc nhanh chóng, chính xác và tiết kiệm chi phí các phân tử hoạt tính tiềm năng từ cây Câu đằng. Đề tài được xác định thực hiện với 2 mục tiêu cụ thể:

  1. Sàng lọc các hợp chất trong cây Câu đằng có tác dụng ức chế A2A adenosine receptor / Monoamine oxidase 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 và tính toán các thông số dược động học và độc tính của các hợp chất tốt nhất thu được sau quá trình sàng lọc.

Đối tượng nghiên cứu của khóa luận là 80 hợp chất hóa học có trong cây Câu đằng (Uncaria rhynchophylla) được thu thập từ các tài liệu khoa học công bố, tương tác với 2 đích phân tử thụ thể protein là A2A adenosine receptor (A2AAR) và enzym monoamine oxidase B (MAO-B). Phạm vi nghiên cứu được thực hiện hoàn toàn trên mô hình tính toán in silico, sử dụng cơ sở dữ liệu cấu trúc thực nghiệm tinh thể từ Ngân hàng Dữ liệu Protein (PDB) và Ngân hàng Dữ liệu Hóa học PubChem, triển khai 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 vào năm 2022.

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

Khung lý thuyết và mô hình sinh học phân tử

Nghiên cứu tập trung vào hai đích phân tử trị liệu quan trọng đối với bệnh Parkinson:

  • Thụ thể A2A Adenosine (A2AAR): Là một thụ thể bắt cặp protein G có cấu trúc gồm 7 vòng xoắn xuyên màng, liên kết với nhau qua 3 vòng nội bào (IL1, IL2, IL3) và 3 vòng ngoại bào (EL1, EL2, EL3). A2AAR tương tác đối kháng sinh lý với thụ thể dopamine tại thể vân. Sự phong bế hoặc đối kháng A2AAR giúp tăng cường tín hiệu dẫn truyền dopaminergic, bù đắp lượng dopamine suy giảm và làm giảm các tác dụng phụ loạn động do dùng L-dopa kéo dài. Bên cạnh đó, các chất đối kháng A2AAR còn thể hiện hoạt tính bảo vệ thần kinh thông qua việc giảm stress oxy hóa và hạn chế tạo gốc tự do. Cấu trúc tinh thể 3D sử dụng là phức hợp A2AAR gắn với phối tử đồng kết tinh ZMA (PDB ID: 3EML).
  • Enzym Monoamine Oxidase B (MAO-B): Là enzym liên kết màng ngoài của ty thể, chịu trách nhiệm thoái hóa dopamine tại khe tiếp hợp và trong tế bào thần kinh đệm, chiếm 80% tổng hoạt tính MAO trong não. Quá trình khử amin oxy hóa dopamine do MAO-B xúc tác sản sinh ra hydrogen peroxide ($H_2O_2$) và các gốc oxy hóa phản ứng, gây tổn thương ty thể; đồng thời MAO-B còn chuyển hóa tiền chất MPTP thành độc tố thần kinh MPP+. Do đó, ức chế MAO-B giúp duy trì nồng độ dopamine nội sinh và ngoại sinh, giảm sản sinh các gốc tự do, bảo vệ chức năng màng ty thể và ức chế sự kết tụ của $\alpha$-synuclein. Cấu trúc tinh thể 3D sử dụng là enzym MAO-B người gắn phối tử đồng kết tinh safinamide (SAG) (PDB ID: 2V5Z).
  • Quy tắc 5 của Lipinski (Rules of 5 - RO5): Bộ tiêu chuẩn sàng lọc đặc tính "giống thuốc" (drug-likeness) đối với các hợp chất phát triển theo đường uống. Một phân tử được xem là có tiềm năng hấp thu tốt qua đường tiêu hóa nếu không vi phạm 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) nằm trong khoảng 40–130.
  • Mô hình đánh giá ADMET: Dự đoán các đặc tính dược động học và độ an toàn in silico bao gồm: hấp thu (độ tan LogS, tính thấm qua màng tế bào Caco-2, tỷ lệ hấp thu qua ruột người HIA); phân bố (thể tích phân bố VD, tính thấm qua hàng rào máu não BBB, tỷ lệ gắn protein huyết tương PPB); chuyển hóa (khả năng ức chế các enzym cytochrom 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); độc tính (độc tính trên gan, thử nghiệm đột biến gen AMES, kích ứng da).

Phương pháp nghiên cứu và công cụ kỹ thuật

Quy trình nghiên cứu gồm các bước xử lý và công cụ cụ thể:

  • Chuẩn bị cấu trúc phân tử: Cấu trúc 3D của 80 hợp chất từ cây Câu đằng được tải từ cơ sở dữ liệu PubChem, sau đó tối ưu hóa hình học không gian và cực tiểu hóa năng lượng bằng phần mềm Avogadro thông qua thuật toán Gradient liên hợp (Conjugate Gradients), trước khi chuyển đổi sang định dạng .pdbqt bằng AutoDock Tools 1.5.6. Cấu trúc 3D của hai protein đích (3EML và 2V5Z) được tải từ Protein Data Bank (RCSB PDB), tiến hành loại bỏ các phân tử nước, các ion đệm và tách phối tử đồng kết tinh bằng Discovery Studio 2021 Client; sau đó bổ sung nguyên tử hydro phân cực, tính điện tích bán thực nghiệm và gán trường lực Kollman bằng AutoDock Tools.
  • Thẩm định quy trình docking (Re-docking): Tiến hành dock lại hai phối tử đồng kết tinh ZMA và SAG vào vị trí hoạt động ban đầu của phức hợp 3EML và 2V5Z. Sử dụng phần mềm UCSF Chimera 1.14 để tính toán độ lệch bình phương trung bình gốc (RMSD). Tiêu chuẩn quy định quy trình docking đạt độ tin cậy khi RMSD < 1,5 Å.
  • Thiết lập thông số hộp lưới tìm kiếm (Grid box):
    • Đích A2AAR: Kích thước hộp là 70 $\times$ 70 $\times$ 80 Å, tọa độ tâm x = -5,968, y = -10,310, z = 53,267.
    • Đích MAO-B: Kích thước hộp là 70 $\times$ 70 $\times$ 70 Å, tọa độ tâm x = 51,661, y = 155,622, z = 27,829.
    • Bước lưới (spacing) của cả hai hệ thống được cố định ở giá trị 0,375 Å.
  • Mô phỏng docking và phân tích tương tác: Quá trình docking được thực hiện bằng phần mềm AutoDock Vina với hàm tính điểm (scoring function) năng lượng tự do liên kết ($\Delta G$, kcal/mol). Phân tích mạng lưới liên kết hydro, tương tác kỵ nước, liên kết ion và lực van der Waals ở không gian 2D và 3D bằng Discovery Studio 2021 Client.
  • Sàng lọc đặc tính giống thuốc và ADMET: Đánh giá Lipinski RO5 thông qua công cụ trực tuyến của Supercomputing Facility for Bioinformatics & Computational Biology (Scfbio). Tính toán các chỉ số ADMET thông qua nền tảng trực tuyến ADMETlab bằng cách nhập cấu trúc dạng SMILES (Simplified Molecular Input Line Entry Specification). Toàn bộ tính toán được chạy trên hệ thống máy tính cá nhân Dell Vostro 3510.

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 toàn diện các y văn liên quan đến bệnh lý Parkinson, các đích tác dụng sinh học phân tử, nguồn dược liệu Câu đằng và phương pháp luận thiết kế thuốc in silico. Tác giả trình bày chi tiết cơ chế bệnh sinh của Parkinson gắn liền với sự thoái hóa nơ-ron dopaminergic và vai trò của thể Lewy chứa $\alpha$-synuclein; phân tích 5 giai đoạn tiến triển lâm sàng và tiêu chuẩn vàng chẩn đoán dựa trên thăm khám dấu hiệu vận động và phi vận động. Phần tổng quan điều trị phân tích rõ cơ chế hoạt động của Levodopa và các thuốc bổ trợ (ức chế COMT, ức chế MAO-B), đồng thời chỉ ra nhược điểm lớn của các phác đồ hiện nay khi chỉ đóng vai trò điều trị triệu chứng vận động mà không ngăn chặn được tiến trình thoái hóa thần kinh.

Chương 1 đi sâu phân tích cơ chế phân tử của thụ thể A2AAR và enzym MAO-B. Đối kháng A2AAR mang lại tác động kép: bù đắp suy giảm chức năng dopamine ở thể vân, giảm tác dụng phụ loạn động của L-dopa và bảo vệ nơ-ron chống lại stress oxy hóa. Trong khi đó, ức chế MAO-B ngăn chặn quá trình khử amin tạo $H_2O_2$, ức chế độc tính của MPP+, cải thiện tính thấm màng ty thể và làm chậm giai đoạn tạo mầm kết tụ $\alpha$-synuclein. Về nguồn dược liệu, chương này tổng kết đặc điểm thực vật, phân bố, kinh nghiệm dùng thuốc y học cổ truyền và các phát hiện hóa thực vật hiện đại trên cây Câu đằng (Uncaria rhynchophylla), đặc biệt là nhóm hoạt chất alkaloid có khả năng can thiệp vào các thụ thể NMDA và 5-HT2A. Cuối cùng, chương 1 giới thiệu nguyên lý của phương pháp docking phân tử, quy tắc Lipinski và các công cụ dự đoán ADMET trong việc sàng lọc ảo các hợp chất tự nhiên.

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

Chương 2 mô tả chi tiết quy trình thực nghiệm trên máy tính, danh mục trang thiết bị, phần mềm và các bước chuẩn bị dữ liệu. Đối tượng nghiên cứu là 80 hợp chất có trong cây Câu đằng cùng hai cấu trúc tinh thể protein đại diện từ RCSB PDB: A2AAR (PDB ID: 3EML, phối tử ZMA) và MAO-B (PDB ID: 2V5Z, phối tử SAG). Thiết bị thực hiện là máy tính Dell Vostro 3510 với các phần mềm: MGL Tools 1.5.6, AutoDock Vina, UCSF Chimera 1.14, Discovery Studio 2021 Client, Avogadro, công cụ Lipinski trực tuyến của Scfbio và nền tảng ADMETlab.

Nội dung nghiên cứu được triển khai mạch lạc qua 3 bước:

  • Bước 1: Chuẩn bị cấu trúc phân tử protein và phối tử, tiến hành thẩm định quy trình re-docking với ZMA và SAG để kiểm tra độ tin cậy của thuật toán qua hệ số RMSD (< 1,5 Å); sau đó tiến hành docking sàng lọc 80 hợp chất vào khoang gắn cơ chất của A2AAR và MAO-B theo các thông số grid box đã thiết lập.
  • Bước 2: Lựa chọn các hợp chất có điểm số năng lượng liên kết thấp hơn chất chứng dương trên cả hai đích để đưa vào đánh giá đặc tính "giống thuốc" theo 5 tiêu chí của Lipinski.
  • Bước 3: Nạp chuỗi SMILES của các hợp chất đạt chuẩn Lipinski vào hệ thống ADMETlab để dự đoán các thông số hấp thu, phân bố, chuyển hóa, thải trừ và độc tính; kết hợp phân tích tương tác không gian 3D tại trung tâm hoạt động để lựa chọn các ứng viên tiềm năng nhất.

Chương 3: Kết quả

Chương này trình bày toàn bộ số liệu thu được từ các giai đoạn sàng lọc tính toán. Kết quả re-dock cho thấy cấu dạng của phối tử ZMA và SAG sau khi dock lại chồng khít tốt lên cấu trúc tinh thể ban đầu với giá trị RMSD lần lượt là 1,24 Å và 1,382 Å (đều < 1,5 Å), cùng mức năng lượng liên kết tự do là -7,7 kcal/mol đối với A2AAR và -9,8 kcal/mol đối với MAO-B. Điều này chứng minh quy trình và các thông số cài đặt docking đạt độ chính xác cao.

Kết quả mô phỏng docking trên 80 hợp chất ghi nhận:

  • 38 hợp chất có năng lượng liên kết trên A2AAR âm hơn chất chứng ZMA (-7,7 kcal/mol).
  • 11 hợp chất có năng lượng liên kết trên MAO-B âm hơn hoặc bằng chất chứng SAG (-9,8 kcal/mol).
  • 9 hợp chất thể hiện tiềm năng ức chế kép nổi trội nhất trên cả 2 đích phân tử (năng lượng liên kết âm hơn cả hai chất chứng).
STT Tên hợp chất PubChem ID Năng lượng liên kết A2AAR (3EML) (kcal/mol) Năng lượng liên kết MAO-B (2V5Z) (kcal/mol)
-- Chứng dương ZMA -- -7,7 --
-- 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

Đánh giá đặc tính "giống thuốc" theo quy tắc 5 Lipinski cho 9 hợp chất tiềm năng cho thấy có 6/9 hợp chất thỏa mãn các tiêu chí (Angustidine, Cinchonain Id, Acid cryptochlorogenic, Acid hyptatic A, Acid uncarinic B, Acid ursolic lactone). Ba hợp chất bị loại gồm: 3$\alpha$-Dihydrocadambine (MW = 546 Da, HBD = 6, HBA = 11, MR = 133,74), Saponin A (MW = 957 Da, HBD = 11, HBA = 18, MR = 231) và Timosaponin A-III (MW = 741 Da, HBD = 7, HBA = 13, MR = 183). Đáng chú ý, Angustidine thỏa mãn cả 5/5 tiêu chí của Lipinski.

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

Dự đoán thông số ADMET của 6 hợp chất đạt quy tắc Lipinski qua ADMETlab được tổng hợp cụ thể:

Thông số Angustidine Cinchonain Id Acid cryptochlorogenic Acid hyptatic A Acid uncarinic B Acid ursolic lactone
HẤP THU
Độ hòa tan LogS ($\log,\mu\text{g/mL}$) 26,246 56,825 3849,33 0,673 2,774 0,126
Tính thấm màng 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
PHÂN BỐ
Thể tích phân bố ($\log,\text{L/kg}$) 0,698 -0,800 -1,279 -0,574 -0,490 -0,343
Tính thấm BBB ($\log,\text{BB}$) +++ ++ --- + --- ++
Liên kết protein huyết tương (%) 81,347 85,633 50,151 88,088 85,122 82,175
CHUYỂN HÓA
Ức chế CYP1A2 ++ --- --- --- --- ---
Ức chế CYP3A4 - --- --- --- + ---
Ức chế CYP2C9 --- --- --- --- - ---
Ức chế CYP2C19 --- - --- --- --- ---
Ức chế CYP2D6 --- - --- --- - ---
THẢI TRỪ
$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
Độc tính gan ++ - + --- - ---
Độc tính AMES + - --- --- --- ---
Kích ứng da --- --- - --- --- -

Dựa trên kết quả dự đoán ADMET, Angustidine bộc lộ độc tính gan (++) và độc tính gây đột biến AMES (+); Acid cryptochlorogenic có tỷ lệ hấp thu qua ruột thấp (32,5%), thời gian bán thải rất ngắn ($T_{1/2} = 0,473$ giờ), không qua được hàng rào máu não (BBB ---) và dự đoán có độc tính gan (+). Do đó, tác giả đã sàng lọc được 4 hợp chất tối ưu nhất có dược động học khả quan và hồ sơ an toàn tốt: Cinchonain Id, Acid hyptatic A, Acid uncarinic B và Acid ursolic lactone.

Phân tích tương tác 3D tại khoang hoạt động của 4 hợp chất này cho thấy:

  • Trên đích A2AAR: Các hợp chất tạo liên kết tương đồng với phối tử đồng kết tinh ZMA, tương tác với các acid amin quan trọng như Leu-167, Ala-63, Ile-274, Asp-170, Phe-168, Met-270, Tyr-271, Glu-169, Leu-267, Trp-246, Leu-249. Cinchonain Id tạo nhiều liên kết acid amin nhất với điểm docking âm nhất (-10,2 kcal/mol). Acid hyptatic A có điểm số -7,9 kcal/mol do tạo một tương tác không thuận lợi với Gly-5 tại nhóm hydroxyl. Acid ursolic lactone chủ yếu tạo tương tác van der Waals do chỉ có một nhóm hydroxyl nhưng đạt mức năng lượng -9,3 kcal/mol.
  • Trên đích MAO-B: Cả 4 hợp chất đều tạo các liên kết hydro, $\pi$-$\pi$, $\pi$-alkyl và van der Waals với các acid amin thiết yếu trong khoang gắn phối tử tương tự như chất chứng SAG (gồm Pro-102, Leu-171, Cys-172, Ile-198, Ile-199, Gln-206, Ile-316, Tyr-326, Phe-343, Tyr-398, Tyr-435 và vòng isoalloxazine của FAD). Mức năng lượng liên kết của 4 hợp chất tương đương với SAG (-9,8 kcal/mol), trong đó Acid ursolic lactone đạt -9,9 kcal/mol.

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

Chương này tập trung thảo luận sâu về cấu trúc hóa học, mối liên hệ cấu trúc - hoạt tính, tiềm năng dược lý và triển vọng phát triển thành thuốc của 4 hợp chất được lựa chọn:

  • Cinchonain Id: Là một hợp chất thuộc nhóm flavonoid phenolic, xuất hiện trong cây Câu đằng và một số thực vật khác như Trichilia catigua, Eriobotrya japonica. Nghiên cứu trước đây của Wanxia Tang và cộng sự cho thấy hợp chất này có hoạt tính chống oxy hóa mạnh với giá trị $IC_{50} = 2,3,\mu\text{M}$. Đây là công trình đầu tiên ghi nhận Cinchonain Id có khả năng gắn kết mạnh lên cả A2AAR (-10,2 kcal/mol) và MAO-B (-9,8 kcal/mol). Với 6 nhóm hydroxyl liên kết trên các vòng thơm, phân tử hình thành mạng lưới liên kết hydro và tương tác $\pi$ phong phú tại tâm hoạt động. Hợp chất đạt 4/5 tiêu chí Lipinski, có khả năng xuyên qua hàng rào máu não (BBB ++), không gây độc gan, không gây đột biến gen AMES, không kích ứng da và không ức chế hệ enzym CYP. Hạn chế duy nhất là tỷ lệ hấp thu qua ruột ở mức trung bình (45,8%), đòi hỏi các nghiên cứu cải thiện sinh khả dụng.
  • Acid hyptatic A: Là hợp chất triterpenoid từng được phân lập từ Hyptis capitata, Premna microphylla và Uncaria rhynchophylla. Điểm docking trên A2AAR và MAO-B lần lượt là -7,9 và -9,8 kcal/mol. Hợp chất thể hiện tính thấm màng Caco-2 tốt (-5,008 logPapp cm/s), hấp thu qua ruột đạt 57,9%, có khả năng qua hàng rào máu não (BBB +), không bị chuyển hóa qua gan, không ức chế các enzym CYP và cho kết quả âm tính đối với toàn bộ các phép thử độc tính.
  • Acid uncarinic B: Là hợp chất tự nhiên nhóm triterpenoid phân lập từ gai Câu đằng và Hibiscus taiwanensis. Trước đây, nghiên cứu của Lee và cộng sự ghi nhận chất này ức chế phospholipase $C\gamma1$ với $IC_{50} = 44,55,\mu\text{M}$, nhưng chưa từng được đánh giá trên các đích thần kinh. Điểm docking trên A2AAR đạt -8,4 kcal/mol và trên MAO-B đạt -9,8 kcal/mol. Phân tử có tỷ lệ hấp thu qua ruột khá cao (69,2%) và không độc tính. Tuy nhiên, nhược điểm lớn của Acid uncarinic B là ức chế enzym CYP3A4 (+) và không có khả năng vượt qua hàng rào máu não (BBB ---), do đó cần nghiên cứu tối ưu hóa cấu trúc hoặc sử dụng hệ dẫn thuốc nano để đưa hoạt chất vào hệ thần kinh trung ương.
  • Acid ursolic lactone: Là một dẫn chất triterpenoid có trong Eucalyptus camaldulensis, Eucalyptus tereticornis, Tecoma stans, từng được nghiên cứu in vitro và in vivo về tiềm năng điều trị đái tháo đường và béo phì. Đây là lần đầu tiên hợp chất được khám phá tác dụng trên bệnh Parkinson với điểm gắn kết rất thuận lợi (-9,3 kcal/mol trên A2AAR và -9,9 kcal/mol trên MAO-B). Hồ sơ ADMET của hợp chất rất tối ưu: tỷ lệ hấp thu qua ruột cao (80,4%), tính thấm Caco-2 tốt (-4,745 logPapp cm/s), thấm tốt qua hàng rào máu não (BBB ++), không ức chế các isoenzym CYP và không thể hiện độc tính gan, kích ứng da hay độc tính AMES.

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

Kết quả chính

  • Sàng lọc ảo toàn diện 80 hợp chất từ cây Câu đằng (Uncaria rhynchophylla), xác định được 38 hợp chất có ái lực gắn kết với thụ thể A2AAR cao hơn chất chứng dương ZMA (-7,7 kcal/mol) và 11 hợp chất có ái lực với enzym MAO-B tốt hơn hoặc bằng chất chứng dương SAG (-9,8 kcal/mol).
  • Phát hiện 9 hợp chất có tiềm năng ức chế kép trên cả hai đích A2AAR và MAO-B với năng lượng tự do liên kết âm hơn cả hai phối tử đồng kết tinh, trong đó Timosaponin A-III cho năng lượng liên kết âm nhất trên A2AAR (-14,9 kcal/mol) và Saponin A cho năng lượng liên kết âm nhất trên MAO-B (-14,2 kcal/mol).
  • Đánh giá đặc tính "giống thuốc" theo quy tắc 5 Lipinski đã chọn lọc được 6/9 hợp chất đạt tiêu chuẩn (Angustidine, Cinchonain Id, Acid cryptochlorogenic, Acid hyptatic A, Acid uncarinic B, Acid ursolic lactone) và loại bỏ 3 hợp chất vi phạm nghiêm trọng về khối lượng phân tử và liên kết hydro (3$\alpha$-Dihydrocadambine, Saponin A, Timosaponin A-III).
  • Tính toán đầy đủ hồ sơ dược động học và an toàn ADMET trên 6 hợp chất giống thuốc, từ đó loại bỏ 2 hợp chất có độc tính hoặc phân bố kém (Angustidine gây độc gan và đột biến gen; Acid cryptochlorogenic hấp thu thấp và không qua BBB), xác định được 4 hợp chất tiềm năng nhất để phát triển thành thuốc: Cinchonain Id, Acid hyptatic A, Acid uncarinic B và Acid ursolic lactone.

Đóng góp so với các nghiên cứu trước

  • Đây là nghiên cứu đầu tiên áp dụng phương pháp docking phân tử để đánh giá một cách hệ thống tác dụng ức chế đồng thời hai đích phân tử then chốt trong điều trị Parkinson (A2AAR và MAO-B) đối với các hợp chất trong cây Câu đằng.
  • Mở rộng phạm vi tìm kiếm hoạt chất chống Parkinson từ cây Câu đằng: trong khi các y văn trước đây chủ yếu tập trung vào nhóm alkaloid (rhynchophylline, isorhynchophylline), công trình này đã chỉ ra các phân tử thuộc nhóm flavonoid (Cinchonain Id) và triterpenoid (Acid hyptatic A, Acid uncarinic B, Acid ursolic lactone) mới là những chất thể hiện ái lực gắn kết kép mạnh mẽ, mở ra hướng tiếp cận hóa thực vật mới.
  • Cung cấp bộ dữ liệu tính toán hoàn chỉnh về thông số docking, mạng lưới tương tác acid amin, các chỉ số hóa lý Lipinski và các đặc tính ADMET làm cơ sở định hướng cho các thử nghiệm bán tổng hợp và chiết tách dược liệu.

Giải pháp và kiến nghị đề xuất của tác giả

  • Tiến hành chiết xuất, phân lập các hợp chất tiềm năng từ cây Câu đằng để triển khai các thử nghiệm sinh học thực nghiệm in vitro (đo hoạt tính ức chế enzym MAO-B, đo ái lực gắn kết thụ thể A2AAR) và in vivo (trên mô hình động vật Parkinson gây độc bởi MPTP hoặc 6-OHDA) nhằm xác minh lại kết quả tính toán mô phỏng.
  • Thực hiện các nghiên cứu về hóa dược và bào chế nhằm tối ưu hóa cấu trúc: cải thiện độ tan, nâng cao tỷ lệ hấp thu đường tiêu hóa và đặc biệt là cải thiện tính thấm qua hàng rào máu não đối với Acid uncarinic B.
  • Mở rộng mô hình sàng lọc in silico trên các đích phân tử trị liệu khác liên quan đến cơ chế bệnh sinh Parkinson, cũng như khảo sát thêm các nhóm hợp chất và cây thuốc tiềm năng khác.

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

Hạn chế của nghiên cứu

  • Toàn bộ kết quả của đề tài mới dừng lại ở mô hình tính toán mô phỏng in silico, chưa được phối hợp kiểm chứng bằng các phương pháp thực nghiệm sinh học in vitro (đo $IC_{50}$) hoặc in vivo trên cơ thể sống.
  • Quá trình docking sử dụng mô hình thụ thể tĩnh (rigid receptor), chưa tính đến các thay đổi về cấu hình không gian linh động của chuỗi polypeptide theo thời gian (vốn cần được đánh giá bổ sung bằng mô phỏng động lực học phân tử - Molecular Dynamics Simulation).
  • Một số hợp chất tiềm năng vẫn tồn tại nhược điểm dược động học nhất định: Acid uncarinic B không có khả năng vượt qua hàng rào máu não (BBB ---) và ức chế enzym CYP3A4, trong khi Cinchonain Id có tỷ lệ hấp thu qua ruột chưa thật sự cao (45,8%).

Hướng nghiên cứu tiếp theo

  • Triển khai quy trình chiết xuất và phân lập hoạt chất từ dược liệu Câu đằng thực tế để tiến hành khảo sát thực nghiệm in vitro và in vivo.
  • Ứng dụng mô phỏng động lực học phân tử (MD Simulation) để đánh giá độ bền nhiệt động học và tính ổn định của phức hợp protein - ligand trong môi trường dung môi mô phỏng sinh học.
  • Nghiên cứu biến đổi cấu trúc phân tử hoặc ứng dụng công nghệ nano, chất mang phân tử để cải thiện sinh khả dụng và tính thấm qua hàng rào máu não cho các hợp chất triterpenoid.

Giá trị tham khảo

Khóa luận là nguồn tài liệu tham khảo chuyên môn hữu ích cho:

  • Giảng viên, sinh viên và học viên sau đại học chuyên ngành Dược học, Hóa dược, Công nghệ sinh học và Y học cổ truyền đang thực hiện các đề tài nghiên cứu về phát triển thuốc từ nguồn dược liệu tự nhiên.
  • Các nhà nghiên cứu quan tâm đến lĩnh vực thiết kế thuốc hỗ trợ bởi máy tính (Computer-Aided Drug Design - CADD), đặc biệt là quy trình sàng lọc ảo kết hợp giữa AutoDock Vina, quy tắc Lipinski và dự đoán ADMET.
  • Các nội dung có giá trị tham khảo ứng dụng cao trong văn bản gồm:
    • Bộ thông số chuẩn hóa tọa độ tâm và kích thước hộp lưới tìm kiếm (grid box) cho hai đích protein thụ thể A2AAR (3EML) và MAO-B (2V5Z).
    • Quy trình thẩm định tính chính xác của thuật toán docking bằng phương pháp re-dock với tiêu chuẩn sai số RMSD < 1,5 Å.
    • Dữ liệu định lượng chi tiết về năng lượng gắn kết, cấu trúc tương tác acid amin và bảng thông số ADMET của 6 hợp chất tự nhiên từ cây Câu đằng.

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

1. Quy trình thẩm định độ tin cậy của mô hình docking phân tử (re-dock) trong khóa luận được thực hiện như thế nào và đạt kết quả ra sao?

Tác giả đã tách hai phối tử đồng kết tinh thể ban đầu là ZMA (trong phức hợp A2AAR, PDB ID: 3EML) và SAG (trong phức hợp MAO-B, PDB ID: 2V5Z) ra khỏi cấu trúc protein, sau đó tiến hành dock ngược lại vào khoang hoạt động bằng AutoDock Vina. Kết quả cho thấy sự chồng khít cấu trúc không gian trước và sau khi dock đạt giá trị RMSD lần lượt là 1,24 Å đối với ZMA và 1,382 Å đối với SAG. Cả hai giá trị này đều nhỏ hơn ngưỡng quy ước 1,5 Å, chứng minh quy trình và các thông số cài đặt hộp lưới tìm kiếm đạt độ tin cậy cao để sàng lọc 80 hợp chất.

2. Dựa vào những tiêu chí nào tác giả đã loại bỏ 3 trong số 9 hợp chất có điểm docking tốt nhất ở bước đánh giá "giống thuốc"?

Tác giả sử dụng 5 tiêu chí của quy tắc Lipinski (MW < 500 Da, HBD $\le$ 5, HBA $\le$ 10, LogP < 5, MR trong khoảng 40–130). Ba hợp chất bị loại bỏ vì vi phạm nghiêm trọng các tiêu chí cấu trúc cho thuốc đường uống:

  • 3$\alpha$-Dihydrocadambine: vi phạm 4 tiêu chí (MW = 546 Da, HBD = 6, HBA = 11, MR = 133,74).
  • Saponin A: vi phạm 4 tiêu chí (MW = 957 Da, HBD = 11, HBA = 18, MR = 231).
  • Timosaponin A-III: vi phạm 3 tiêu chí (MW = 741 Da, HBD = 7, HBA = 13, MR = 183).

3. Vì sao hai hợp chất Angustidine và Acid cryptochlorogenic đạt quy tắc Lipinski nhưng không được chọn vào nhóm 4 hợp chất tối ưu cuối cùng?

Mặc dù Angustidine đạt cả 5/5 tiêu chí Lipinski và Acid cryptochlorogenic đạt 4/5 tiêu chí, kết quả dự đoán ADMETlab cho thấy:

  • Angustidine bộc lộ độc tính gan (++), độc tính đột biến gen AMES (+) và có khả năng ức chế mạnh isoenzym chuyển hóa CYP1A2 (++).
  • Acid cryptochlorogenic có tỷ lệ hấp thu qua ruột thấp nhất trong các chất (32,5%), thời gian bán thải quá ngắn ($T_{1/2} = 0,473$ giờ), không có khả năng qua hàng rào máu não (BBB ---) và bị dự đoán có độc tính trên gan (+).

4. Bốn hợp chất tiềm năng nhất được lựa chọn có những ưu điểm nổi bật nào về mặt gắn kết và dược động học?

  • Cinchonain Id: Năng lượng liên kết A2AAR là -10,2 kcal/mol, MAO-B là -9,8 kcal/mol; có khả năng qua hàng rào máu não (BBB ++), không ức chế enzym chuyển hóa CYP, không độc tính gan/da/AMES và có hoạt tính chống oxy hóa đã biết ($IC_{50} = 2,3,\mu\text{M}$).
  • Acid hyptatic A: Năng lượng liên kết A2AAR là -7,9 kcal/mol, MAO-B là -9,8 kcal/mol; tính thấm Caco-2 tốt (-5,008 logPapp cm/s), hấp thu ruột 57,9%, qua được hàng rào máu não (BBB +) và hoàn toàn không thể hiện độc tính.
  • Acid uncarinic B: Năng lượng liên kết A2AAR là -8,4 kcal/mol, MAO-B là -9,8 kcal/mol; hấp thu ruột đạt 69,2%, không độc tính và từng được biết có hoạt tính ức chế phospholipase $C\gamma1$.
  • Acid ursolic lactone: Năng lượng liên kết A2AAR là -9,3 kcal/mol, MAO-B là -9,9 kcal/mol; hấp thu ruột vượt trội (80,4%), tính thấm Caco-2 cao (-4,745 logPapp cm/s), thấm tốt qua hàng rào máu não (BBB ++), không ức chế các enzym CYP và an toàn tuyệt đối với gan, da và gen.

5. Hai đích phân tử A2AAR và MAO-B giữ vai trò sinh học cụ thể nào trong cơ chế điều trị bệnh Parkinson?

Thụ thể A2AAR phân bố ở thể vân và có tương tác đối kháng sinh lý với thụ thể dopamine; ức chế hoặc đối kháng A2AAR giúp tăng cường dẫn truyền dopaminergic, bù đắp lượng dopamine thiếu hụt và làm giảm biến chứng loạn động do dùng L-dopa kéo dài, đồng thời bảo vệ nơ-ron trước các gốc oxy hóa tự do. Enzym MAO-B chịu trách nhiệm thoái hóa 80% dopamine trong não và tạo ra sản phẩm phụ độc hại $H_2O_2$, đồng thời chuyển hóa MPTP thành độc tố MPP+; do đó, ức chế MAO-B giúp bảo tồn nồng độ dopamine, bảo vệ màng ty thể khỏi stress oxy hóa và làm chậm quá trình kết tụ thể Lewy chứa $\alpha$-synuclein.

Kết luận

Khóa luận tốt nghiệp của tác giả Cao Thị Lan Anh đã ứng dụng thành công phương pháp docking phân tử kết hợp sàng lọc quy tắc 5 Lipinski và mô hình dự đoán ADMET để đánh giá toàn diện 80 hợp chất từ cây Câu đằng (Uncaria rhynchophylla) hướng tới điều trị bệnh Parkinson. Công trình đã thẩm định quy trình tính toán đạt độ tin cậy cao với giá trị RMSD nhỏ hơn 1,5 Å và xác định được 4 hoạt chất tiềm năng nhất có khả năng ức chế kép trên cả hai đích A2AAR và MAO-B cùng hồ sơ an toàn thuận lợi gồm Cinchonain Id, Acid hyptatic A, Acid uncarinic B và Acid ursolic lactone. Kết quả nghiên cứu là tiền đề khoa học có giá trị định hướng cho các nghiên cứu thực nghiệm chiết xuất, thử nghiệm hoạt tính sinh học in vitro, in vivo và phát triển các dạng thuốc mới hỗ trợ điều trị bệnh thoái hóa thần kinh từ nguồn dược liệu truyền thống.