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

Đái tháo đường (ĐTĐ), đặc biệt là đái tháo đường typ 2 (ĐTĐ2), đang là một trong những thách thức y tế công cộng lớn nhất trên phạm vi toàn cầu. Theo báo cáo năm 2019 của Liên đoàn Đái tháo đường Quốc tế (IDF), toàn thế giới có khoảng 463 triệu người trong độ tuổi từ 20 đến 79 mắc bệnh, dự báo con số này sẽ gia tăng lên 578 triệu người vào năm 2030 và chạm ngưỡng 700 triệu người vào năm 2045. Trong năm 2019, hơn 4 triệu trường hợp tử vong có liên quan trực tiếp đến các biến chứng mạn tính của bệnh như suy thận, tim mạch, mù lòa và đoạn chi. Tại Việt Nam, tình hình dịch tễ ghi nhận xu hướng gia tăng đáng lo ngại: số liệu năm 2017 từ IDF cho thấy có 3,53 triệu người trưởng thành mắc bệnh với ít nhất 80 ca tử vong mỗi ngày do biến chứng; điều tra năm 2012 của Bệnh viện Nội tiết Trung ương ghi nhận tỷ lệ hiện mắc là 5,42% (với 63% chưa được chẩn đoán); và số liệu điều tra STEPwise của Bộ Y tế năm 2015 chỉ ra tỷ lệ mắc ĐTĐ toàn quốc ở nhóm 18–69 tuổi là 4,1%, bên cạnh 3,6% tiền đái tháo đường.

Thị trường thuốc điều trị ĐTĐ đứng thứ hai toàn cầu về doanh số chỉ sau thuốc điều trị ung thư, đạt quy mô 39,2 tỷ USD vào năm 2019 và ước tính chi phí chăm sóc y tế toàn cầu có thể vượt mốc 500 tỷ USD vào năm 2035. Mặc dù nhiều nhóm thuốc hạ đường huyết đường uống đã được đưa vào lâm sàng (như sulfonylurea, biguanide, thiazolidinedione, chất ức chế DPP-4, chất ức chế kênh SGLT-2, chất ức chế α-glucosidase), các liệu pháp này vẫn bộc lộ nhiều hạn chế về độc tính, nguy cơ hạ đường huyết quá mức, tăng cân và tác dụng phụ trên đường tiêu hóa. Quy trình nghiên cứu và phát triển (R&D) một hoạt chất mới theo phương pháp truyền thống tiêu tốn trung bình 1,8 tỷ USD với tỷ lệ đào thải lên tới 96%, trong đó có tới 50% ứng cử viên thất bại ở thử nghiệm lâm sàng do các đặc tính dược động học kém (39%) và độc tính trên động vật (11%). Bên cạnh đó, nghiên cứu tương quan toàn bộ nhiễm sắc thể (GWAS) đã xác định ít nhất 65 biến thể DNA liên quan đến ĐTĐ2, chứng minh tính phức tạp của bệnh sinh và thúc đẩy xu hướng tiếp cận điều trị theo cơ chế đa đích.

Việt Nam sở hữu nguồn dược liệu tự nhiên phong phú với nhiều loài thực vật chứa nhóm hợp chất flavonoid có tiềm năng hạ đường huyết đã được ứng dụng trong y học cổ truyền. Việc ứng dụng công cụ sàng lọc ảo trên máy tính (in silico) dựa trên cấu trúc (Structure-based virtual screening - SBVS) kết hợp đánh giá các thông số dược động học và an toàn phân tử là giải pháp tối ưu giúp rút ngắn thời gian, giảm chi phí và sàng lọc chính xác các phân tử dẫn đường tiềm năng.

Trên cơ sở đó, khóa luận xác định 02 mục tiêu nghiên cứu cụ thể:

  1. Sàng lọc hợp chất flavonoid có tác dụng ức chế đồng thời enzym α-amylase và α-glucosidase hướng điều trị ĐTĐ2 sử dụng mô hình in silico.
  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, độc tính (ADMET) 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 và phạm vi nghiên cứu:

  • Đối tượng phân tử đích: Cấu trúc không gian ba chiều của enzym α-amylase tuyến tụy người (mã PDB: 2QV4, độ phân giải 1,97 Å) và enzym α-glucosidase (mã PDB: 3W37, độ phân giải 1,7 Å).
  • Đối tượng phối tử: Thư viện gồm 73 hợp chất flavonoid tự nhiên được thu thập từ cơ sở dữ liệu PubChem.
  • Đối chứng dương: Thuốc hạ đường huyết chuẩn acarbose (PubChem CID: 41774).
  • Phạm vi nghiên cứu: Nghiên cứu được thực hiện tại Trường Đại học Y Dược, Đại học Quốc gia Hà Nội vào năm 2023, giới hạn trong phạm vi tính toán tin sinh học (in silico).

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

Khung lý thuyết và cơ sở sinh học phân tử

Cơ chế điều hòa đường huyết sau ăn phụ thuộc chặt chẽ vào hoạt động của hai enzym thủy phân carbohydrate chính:

  • Enzym α-amylase (AA): Là một endoenzym giàu tyrosine, acid glutamic và acid aspartic (hai acid amin này chiếm khoảng 1/4 tổng số acid amin), có phân tử lượng trung bình khoảng 50 kDa. Enzym này xúc tác phân cắt ngẫu nhiên các liên kết α-1,4-glucosidic nằm sâu bên trong mạch phân tử tinh bột và glycogen, chuyển hóa cơ chất qua ba giai đoạn: dextrin hóa tạo α-dextrin phân tử thấp, đường hóa tạo tetra- và trimaltose, và cuối cùng tạo thành maltotetrose, maltotriose và maltose.
  • Enzym α-glucosidase (AG): Là một exoenzym khu trú tại diềm bàn chải tế bào niêm mạc ruột non (chủ yếu thuộc các họ glycoside hydrolase GH13 và GH31), chịu trách nhiệm thủy phân chặng cuối cùng các disaccharide và oligosaccharide thành monosaccharide glucose tự do để hấp thu qua màng ruột vào hệ tuần hoàn.

Các chất ức chế enzym α-amylase (AAIs) và α-glucosidase (AGIs) đóng vai trò làm chậm quá trình tiêu hóa và hấp thu carbohydrate, ngăn chặn đỉnh tăng glucose huyết sau ăn mà không kích thích tế bào β tụy tiết insulin, từ đó giảm thiểu nguy cơ hạ đường huyết quá mức. Tuy nhiên, các thuốc AGIs tổng hợp đang lưu hành trên thị trường (acarbose, miglitol, voglibose, DNJ) thường gây tác dụng không mong muốn trên đường tiêu hóa như đầy hơi, co thắt ruột và tiêu chảy.

Nhóm hợp chất flavonoid là các chất chuyển hóa thứ cấp nhóm polyphenol gồm hơn 4.000 cấu trúc, có khung cơ bản là 2-phenyl benzopyrone (C6-C3-C6) liên kết với các nhóm hydroxyl có khả năng dọn gốc tự do, chống oxy hóa, kháng viêm và ức chế các enzym chuyển hóa carbohydrate. Flavonoid được phân thành ba nhóm chính: euflavonoid (2-phenylbenzopyrans), isoflavonoid (3-benzopyrans) và neoflavonoid (4-benzopyrans).

Phương pháp nghiên cứu và quy trình sàng lọc

Nghiên cứu sử dụng hệ thống phần cứng máy tính cá nhân HP Pavilion 14 (hệ điều hành Windows 11) kết hợp hệ thống phần mềm chuyên dụng và công cụ trực tuyến:

  • Chuẩn bị cấu trúc phân tử: Tinh thể protein 2QV4 và 3W37 được tải về từ Protein Data Bank (PDB). Bằng phần mềm Discovery Studio 2021 Client và AutoDockTools 1.5.6, các phân tử nước và ligand đồng kết tinh được loại bỏ, bổ sung nguyên tử hydro phân cực, gán trường lực Kollman và chuyển đổi sang định dạng .pdbqt. Cấu trúc không gian ba chiều của 73 flavonoid từ PubChem được tối ưu hóa năng lượng hình học bằng phương pháp Gradient liên hợp (Conjugate Gradients) trên phần mềm Avogadro 2.
  • Thẩm định quy trình docking (Re-docking): Tiến hành tách phối tử đồng kết tinh (acarbose/NAG/QV4) và dock ngược lại vào khoang hoạt động của protein tương ứng. Sự tương đồng về mặt cấu dạng trước và sau khi docking được đánh giá thông qua chỉ số độ lệch bình phương trung bình gốc (RMSD) bằng phần mềm UCSF Chimera 1.14. Quy trình được công nhận có giá trị tin cậy khi RMSD < 1,5 Å.
  • Mô phỏng docking phân tử: Thực hiện bằng phần mềm AutoDock với các thông số hộp tìm kiếm (grid box):
    • Đối với α-amylase (2QV4): Kích thước hộp 70 × 70 × 70 Å, khoảng cách mắt lưới 0,375 Å, tọa độ tâm x = 16,11; y = 49,57; z = 25,06.
    • Đối với α-glucosidase (3W37): Kích thước hộp 60 × 60 × 70 Å, khoảng cách mắt lưới 0,375 Å, tọa độ tâm x = 14,103; y = -27,101; z = -42,876.
    • Mỗi hợp chất được docking lặp lại 3 lần để lấy giá trị năng lượng tự do liên kết ($\Delta G$, kcal/mol) trung bình, so sánh với giá trị của acarbose.
  • Đánh giá đặc tính giống thuốc (Lipinski's Rule of 5): Thực hiện qua công cụ trực tuyến để kiểm tra 5 chỉ số hóa lý thiết yếu cho thuốc dùng đường uống: trọng lượng phân tử (MW < 500 Da), số nhóm cho liên kết hydro (HBD < 5), số nhóm nhận liên kết hydro (HBA < 10), hệ số phân bố octanol/nước (LogP < 5) và độ khúc xạ mol (MR trong khoảng 40–130).
  • Dự đoán các thông số ADMET: Sử dụng máy chủ pkCSM (dựa trên cấu trúc chuỗi SMILES) để dự đoán tính thấm màng Caco-2, tỷ lệ hấp thu qua ruột người (HIA %), thể tích phân bố ổn định (VDss), khả năng qua hàng rào máu não (BBB) và hệ thần kinh trung ương (CNS), khả năng ức chế enzym cytochrom P450 (CYP2D6, CYP3A4), độ thanh thải toàn phần (Clr), cơ chất vận chuyển cation hữu cơ tại thận (OCT2), độc tính đột biến gen (AMES), ức chế kênh kali hERG (độc tính tim), liều gây chết trung bình (LD50 trên chuột), độc tính trên gan và kích ứng da.

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa cơ sở lý luận về bệnh sinh đái tháo đường typ 2, các số liệu dịch tễ học trong nước và quốc tế, cùng phân loại 4 nhóm bệnh đái tháo đường lâm sàng theo chuẩn y khoa. Tác giả phân tích ưu, nhược điểm của 9 nhóm thuốc điều trị ĐTĐ2 hiện hành (Bảng 1.1 trong khóa luận) gồm: insulin, sulfonylurea, meglitinide, biguanide, thiazolidinedione, chất ức chế α-glucosidase, chất ức chế DPP-4, đồng vận thụ thể GLP-1 và chất ức chế SGLT-2. Chương này đi sâu vào cơ chế xúc tác enzyme của α-amylase và α-glucosidase, cấu trúc hóa học phân nhóm của flavonoid tự nhiên và nguyên lý vận hành của các mô hình sàng lọc ảo dựa trên cấu trúc (SBVS). Tác giả làm rõ các thuật toán tìm kiếm cấu dạng phối tử (thuật toán so khớp, xây dựng tăng dần, Monte Carlo, thuật toán di truyền GA, động lực học phân tử MD), bản chất của các hàm tính điểm năng lượng tự do Gibbs ($\Delta G$) và tầm quan trọng của các bộ lọc sàng lọc tính chất giống thuốc (Quy tắc số 5 Lipinski) cùng mô hình máy học dự đoán dược động học, độc tính (ADMET).

Chương 2: Thiết bị, nội dung và phương pháp nghiên cứu

Chương 2 mô tả chi tiết danh mục phần cứng, các công cụ tính toán và quy trình thực nghiệm gồm ba bước kế tiếp nhau:

  1. Sàng lọc ái lực liên kết của 73 flavonoid trên hai đích tác dụng enzym bằng phương pháp docking phân tử;
  2. Đánh giá tính chất giống thuốc đường uống theo quy tắc số 5 của Lipinski;
  3. Dự đoán các thông số ADMET của các chất vượt qua vòng lọc để xác định các ứng viên tối ưu.

Quy trình thẩm định phương pháp được chuẩn hóa thông qua thử nghiệm re-docking phối tử nội tại vào cấu trúc tinh thể 2QV4 và 3W37. Tọa độ khoang hoạt động và kích thước hộp lưới tính toán cho từng enzym được cố định chính xác dựa trên các acid amin xúc tác quan trọng tại trung tâm hoạt động. Chương này cũng nêu rõ các ngưỡng đánh giá sinh học: tính thấm Caco-2 được coi là cao khi log Papp > 0,9; tỷ lệ hấp thu qua ruột người (HIA) được phân loại kém (< 20%), trung bình (20–70%) và tốt (70–100%); phân bố mô tốt khi log VDss > 0,45; tính thấm thần kinh an toàn khi log BB < 0,3 và log PS < -2.

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 sàng lọc:

  • Thẩm định quy trình: Kết quả re-docking cho thấy cấu dạng phối tử đồng kết tinh tái tạo có độ chồng khít cao với cấu trúc thực nghiệm trong tinh thể X-quang: chỉ số RMSD đối với phức hợp α-amylase (QV4/2QV4) đạt 0,919 Å và đối với α-glucosidase (NAG/3W37) đạt 0,500 Å. Cả hai giá trị đều nhỏ hơn 1,5 Å, khẳng định tính chính xác của thông số và thuật toán docking.
  • Sàng lọc ái lực liên kết: Đối chứng acarbose thể hiện năng lượng liên kết $\Delta G$ là -7,9 kcal/mol trên α-amylase và -8,2 kcal/mol trên α-glucosidase. Trong số 73 flavonoid sàng lọc:
    • Có 67 hợp chất có ái lực liên kết với α-amylase mạnh hơn hoặc bằng acarbose ($\Delta G \le -7,9$ kcal/mol);
    • Có 48 hợp chất có ái lực liên kết với α-glucosidase mạnh hơn hoặc bằng acarbose ($\Delta G \le -8,2$ kcal/mol);
    • Xác định được 47 hợp chất biểu hiện tác dụng kép mạnh trên cả hai enzym.
  • Sàng lọc đặc tính giống thuốc: Áp dụng bộ quy tắc Lipinski trên 47 hợp chất tiềm năng kép, có 14 hợp chất đáp ứng tiêu chuẩn dùng đường uống gồm: Corylin, Cyanidin chloride, Daidzin, Epicatechin, Galangin, Kaempferol, Liquiritin, Luteolin, Maackiain, Myricetin, Ononin, Quercetin, Tectochrysin và Tectorigenin.
  • Dự đoán ADMET và phân tích tương tác: Phân tích toàn diện 14 hợp chất qua pkCSM đã xác định 4 hợp chất tối ưu nhất về khả năng hấp thu qua ruột non, tính thấm Caco-2 cao và hồ sơ an toàn vượt trội: Tectochrysin, Galangin, Maackiain và Corylin. Mô phỏng tương tác không gian 2D/3D bằng Discovery Studio làm rõ các liên kết hydro, tương tác tĩnh điện $\pi$-cation, $\pi$-anion và tương tác kỵ nước $\pi$-alkyl giữa 4 phối tử này với các acid amin cốt lõi trong khoang xúc tác của hai enzym.

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

Chương 4 tập trung biện giải các kết quả tính toán trên cơ sở hóa học dược dụng và phân tích đối chiếu với các y văn quốc tế:

  • Phân tích mối quan hệ giữa cấu trúc và hoạt tính (SAR): Khung polyphenolic của flavonoid với các vòng thơm tạo tương tác xếp chồng $\pi$-$\pi$ bền vững và hệ thống nhóm hydroxyl định hướng liên kết hydro với khoang xúc tác là yếu tố then chốt tạo nên ái lực gắn kết cao hơn chất chuẩn acarbose.
  • Luận giải 4 hợp chất tiềm năng: Tectochrysin (chiết xuất từ chi Populus) có tác dụng bảo vệ tế bào và kéo dài tuổi thọ; Galangin (từ rễ cây củ riềng Alpinia galanga) là chất chống viêm, chống oxy hóa mạnh và có tiềm năng hạ đường huyết đã được ghi nhận; Maackiain (từ dược liệu Sophora flavescens) có tiềm năng điều chỉnh phản ứng viêm nhưng tương tác gắn kết với enzym còn hạn chế, gợi mở hướng bán tổng hợp hóa dược; Corylin (từ hạt phá cố chỉ Psoralea corylifolia) thể hiện khả năng gắn kết đa điểm vượt trội.
  • Đánh giá phương pháp luận: Khẳng định tính hiệu quả và độ tin cậy của thuật toán AutoDock và công cụ pkCSM, đồng thời chỉ rõ sự cần thiết phải phối hợp mô hình in silico với các thử nghiệm thực nghiệm sinh học (in vitro, in vivo).
  • Thảo luận các rào cản kỹ thuật của nghiên cứu máy tính: Sai số tích lũy qua các bước lọc, tính linh động của protein chưa được phản ánh toàn diện, phụ thuộc vào chất lượng cơ sở dữ liệu và yêu cầu kiểm soát chặt chẽ quá trình tiền xử lý dữ liệu.

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

Nghiên cứu đã xây dựng thành công quy trình sàng lọc ảo tích hợp từ docking phân tử đến dự đoán ADMET. Bảng tổng hợp dưới đây chiết xuất các dữ liệu quan trọng nhất của 4 hợp chất dẫn đường tiềm năng được lựa chọn cuối cùng so sánh với chất đối chứng acarbose:

Tên hoạt chất PubChem CID $\Delta G$ với AA (kcal/mol) $\Delta G$ với AG (kcal/mol) MW (g/mol) LogP HBA / HBD Thấm Caco-2 HIA (%) Thần kinh (BBB / CNS) Nguy cơ độc tính (AMES / hERG / Gan) Acid amin tương tác chính trên AA Acid amin tương tác chính trên AG
Corylin 5316097 -9,0 -9,1 320,34 4,19 4 / 1 1,26 96,87 An toàn (-1,65) AMES (+) / hERG (-) / Gan (-) LEU165, TRP59, ASP300, GLU233 ARG670, ARG699, ARG676, GLU301
Galangin 5281616 -8,4 -8,7 270,24 2,60 5 / 3 0,99 93,98 Không thấm (-2,17) Không độc tính GLN63, TYR62, TRP59 GLU310, ARG699, ARG670, ARG676, ASP666
Maackiain 91510 -8,8 -8,3 284,27 2,73 5 / 1 1,24 96,92 Không thấm (-2,06) AMES (+) / hERG (-) / Gan (-) TYR62 ARG699, ARG676, ILE759, ARG670
Tectochrysin 5281954 -8,1 -8,6 268,27 3,02 4 / 1 1,25 95,23 An toàn (-1,99) Không độc tính GLN63, TYR62, TRP59, HIS305 ARG670, ARG699, ARG814, ARG676, ASP666, GLU792, VAL760, TYR659
Acarbose (Đối chứng) 41774 -7,9 -8,2 — — — — — — — ASP300, TRP59, GLU233, GLY164, HIS299 ASN425, ARG422, ASP423, TYR427, PRO426, LEU428

Các đóng góp khoa học chính mà tác giả đã đạt được:

  • Cung cấp dữ liệu định lượng về ái lực liên kết của 73 hợp chất flavonoid tự nhiên trên hai enzym đích liên quan trực tiếp đến hấp thu carbohydrate trong điều trị ĐTĐ2, xác định được tỷ lệ cao các hợp chất có ái lực mạnh hơn thuốc đối chứng lâm sàng acarbose (67 chất trên α-amylase và 48 chất trên α-glucosidase).
  • Thu hẹp thành công thư viện hóa học từ 73 phân tử xuống 47 hợp chất có hoạt tính ức chế kép, tiếp tục qua phễu lọc cấu trúc giống thuốc Lipinski còn 14 chất và cuối cùng xác lập được 4 hợp chất (Corylin, Galangin, Maackiain, Tectochrysin) sở hữu đặc tính dược động học vượt trội (hấp thu qua ruột > 93%, tính thấm màng tế bào cao > 0,9).
  • Giải trình chi tiết cơ chế tương tác ở cấp độ nguyên tử giữa các phân tử flavonoid tiềm năng với các gốc acid amin tại vị trí hoạt động của α-amylase và α-glucosidase, chứng minh Corylin tạo nhiều liên kết bền vững nhất (hydrogen, alkyl, $\pi$-alkyl, $\pi$-anion).

Kiến nghị từ tác giả:

  • Cần triển khai các thử nghiệm thực nghiệm sinh học in vitro (đo hoạt độ enzym ức chế IC50) và in vivo trên mô hình động vật nhằm đánh giá chính xác hoạt tính hạ đường huyết và khoảng nồng độ an toàn sinh học.
  • Đề xuất áp dụng biến đổi cấu trúc, bán tổng hợp hóa học đối với Maackiain để gia tăng số lượng liên kết bền vững với protein đích.
  • Mở rộng phạm vi ứng dụng của mô hình sàng lọc in silico này sang các nhóm hoạt chất thiên nhiên khác và các đích phân tử mới điều trị ĐTĐ2 (như DPP-4, SGLT-2, PTP1B).

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

Khóa luận chỉ rõ các hạn chế nội tại xuất phát từ bản chất của mô hình tính toán in silico:

  • Kết quả thu được hoàn toàn dựa trên các mô phỏng tĩnh và thuật toán dự đoán của phần mềm, chưa phản ánh được đầy đủ sự thay đổi cấu hình động học của enzym trong môi trường sinh lý thực tế của cơ thể người.
  • Việc sử dụng liên tiếp nhiều công cụ và phần mềm khác nhau (Avogadro, AutoDock, pkCSM) tiềm ẩn nguy cơ tích lũy sai số tính toán; một lỗi nhỏ trong định dạng cấu trúc hoặc xử lý trường lực có thể ảnh hưởng đến kết quả cuối cùng.
  • Cảnh báo an toàn phân tử: Hai trong số bốn hợp chất tiềm năng là Corylin và Maackiain cho kết quả dương tính với phép thử AMES trên mô hình pkCSM, biểu hiện nguy cơ đột biến gen/sinh ung thư cần được kiểm chứng nghiêm ngặt; đồng thời Maackiain tạo ít liên kết với trung tâm hoạt động của enzym nên có năng lượng liên kết chưa tối ưu.

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

  • Thực hiện kiểm chứng hoạt tính ức chế enzym trên mô hình thực nghiệm in vitro đối với 4 hợp chất đã chọn lọc.
  • Đánh giá độc tính tế bào và độc tính cấp, bán trường diễn trên động vật thực nghiệm.
  • Nghiên cứu bán tổng hợp hóa học nhằm tối ưu hóa cấu trúc của Maackiain, bổ sung các nhóm chức để tăng cường ái lực liên kết với protein đích.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo hữu ích cho sinh viên, học viên cao học và nghiên cứu viên thuộc các chuyên ngành Dược học, Hóa dược, Dược lý học, Hóa sinh học và Tin sinh học y dược:

  • Phần đáng tham khảo nhất là quy trình từng bước chuẩn hóa mô phỏng docking phân tử trên hai enzym chuyển hóa carbohydrate (α-amylase mã 2QV4 và α-glucosidase mã 3W37), bao gồm thông số tọa độ tâm lưới và kích thước hộp tìm kiếm đã được kiểm định qua bước re-docking với giá trị RMSD đạt chuẩn (< 1,0 Å).
  • Bảng dữ liệu chi tiết về năng lượng tự do liên kết ($\Delta G$) của 73 dẫn xuất flavonoid và ma trận thông số dược động học, an toàn phân tử ADMET cung cấp nguồn dữ liệu đối chiếu giá trị cho các nghiên cứu sàng lọc thuốc mới từ nguồn dược liệu tự nhiên hướng đích ĐTĐ2.

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

1. Vì sao khóa luận lại lựa chọn ức chế đồng thời cả hai enzym α-amylase và α-glucosidase thay vì chỉ tác động lên một đích đơn lẻ?
Quá trình phân giải carbohydrate từ thức ăn diễn ra theo hai giai đoạn nối tiếp: enzym α-amylase thủy phân tinh bột thành các oligosaccharide mạch ngắn, sau đó enzym α-glucosidase tại niêm mạc ruột tiếp tục phân cắt các phân tử này thành đường đơn glucose để hấp thu vào máu. Việc ức chế đồng thời cả hai enzym tạo nên tác động hiệp đồng, làm chậm triệt để quá trình sinh glucose từ thức ăn, kiểm soát hiệu quả đỉnh tăng đường huyết sau ăn, đồng thời phù hợp với định hướng điều trị bệnh lý đa cơ chế phân tử của đái tháo đường typ 2.

2. Độ tin cậy của quy trình docking phân tử trong nghiên cứu này được thẩm định bằng cách nào và kết quả cụ thể ra sao?
Độ tin cậy được thẩm định bằng phương pháp re-docking: tách các phối tử đồng kết tinh trong cấu trúc tinh thể thực nghiệm của enzym (acarbose/QV4/NAG) ra khỏi protein rồi tiến hành dock ngược lại vào khoang hoạt động theo cùng bộ thông số tính toán. Kết quả tính toán bằng phần mềm Chimera cho thấy độ lệch bình phương trung bình gốc (RMSD) đạt 0,919 Å đối với α-amylase (2QV4) và 0,500 Å đối với α-glucosidase (3W37). Cả hai giá trị đều thấp hơn ngưỡng quy chuẩn 1,5 Å, chứng minh quy trình docking đạt độ chính xác cao.

3. Từ 73 hợp chất flavonoid ban đầu, quy trình sàng lọc đã loại bỏ các phân tử qua các phễu lọc nào để chọn ra 4 chất cuối cùng?
Quá trình sàng lọc diễn ra qua 3 phễu lọc liên tiếp:

  1. Phễu năng lượng liên kết docking: Chọn các chất có ái lực mạnh hơn hoặc bằng đối chứng acarbose trên cả 2 enzym ($\Delta G \le -7,9$ kcal/mol trên AA và $\le -8,2$ kcal/mol trên AG), thu được 47 hợp chất.
  2. Phễu tính chất giống thuốc: Đánh giá theo Quy tắc số 5 của Lipinski, chọn các chất đạt ít nhất 2/5 tiêu chuẩn về độ hòa tan và hấp thu đường uống, thu hẹp còn 14 hợp chất.
  3. Phễu dược động học và độc tính (ADMET): Phân tích qua công cụ pkCSM, chọn các phân tử có độ hấp thu ở ruột cao (HIA > 90%), tính thấm Caco-2 cao (> 0,9), an toàn với hệ thần kinh và ít nguy cơ độc tính tạng, xác định được 4 hợp chất tối ưu: Tectochrysin, Galangin, Maackiain và Corylin.

4. Trong 4 hợp chất tiềm năng nhất, hợp chất nào tạo liên kết bền vững nhất với trung tâm hoạt động của enzym và vì sao?
Corylin (PubChem CID: 5316097) là hợp chất thể hiện khả năng tương tác tốt nhất trên cả hai enzym. Kết quả phân tích tương tác 2D/3D cho thấy Corylin tạo được nhiều liên kết hydro, liên kết alkyl, $\pi$-alkyl và $\pi$-anion với các acid amin quan trọng tại vị trí hoạt động: tương tác với LEU165, TRP59, ASP300, GLU233 trên α-amylase và ARG670, ARG699, ARG676, GLU301 trên α-glucosidase. Số lượng liên kết bền vững phong phú này giúp phức hợp giữa Corylin và enzym đạt trạng thái ổn định với mức năng lượng tự do thấp nhất.

5. Kết quả dự đoán ADMET có ghi nhận nguy cơ độc tính nào đối với 4 hợp chất được lựa chọn hay không?
Có. Trong khi Tectochrysin và Galangin hoàn toàn âm tính với tất cả các chỉ số độc tính sàng lọc (không độc tính AMES, không ức chế hERG gây độc tim, không gây độc gan và da), thì Corylin và Maackiain đều bị cảnh báo dương tính với phép thử độc tính AMES (nguy cơ gây đột biến gen). Đây là dữ liệu cảnh báo quan trọng đòi hỏi phải có các nghiên cứu thực nghiệm độc tính chuyên sâu trước khi định hướng phát triển tiền lâm sàng.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Huyền đã hoàn thành trọn vẹn hai mục tiêu nghiên cứu thông qua việc ứng dụng hiệu quả phương pháp sàng lọc ảo in silico kết hợp đánh giá ADMET trên 73 hợp chất flavonoid tự nhiên. Nghiên cứu đã chứng minh 47 hợp chất có khả năng ức chế kép trên cả α-amylase và α-glucosidase vượt trội so với thuốc chuẩn acarbose, đồng thời xác lập được 4 phối tử tiềm năng nhất gồm Tectochrysin, Galangin, Maackiain và Corylin thỏa mãn quy tắc Lipinski cùng các đặc tính hấp thu đường ruột tối ưu. Kết quả đạt được cung cấp cơ sở dữ liệu phân tử xác thực, đóng vai trò là bước định hướng quan trọng để tiếp tục triển khai các thử nghiệm thực nghiệm in vitro và in vivo nhằm tìm kiếm các thuốc hạ đường huyết mới từ nguồn dược liệu tự nhiên.