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

Đái tháo đường (ĐTĐ) là một hội chứng rối loạn chuyển hóa không đồng nhất đặc trưng bởi tình trạng tăng glucose huyết mạn tính do rối loạn bài tiết insulin, giảm hoạt tính của insulin hoặc kết hợp cả hai yếu tố. Theo thống kê của Hiệp hội Đái tháo đường Quốc tế (IDF) năm 2021, toàn cầu có khoảng 537 triệu người trưởng thành (20–79 tuổi) đang sống chung với bệnh đái tháo đường. Tại Việt Nam, đái tháo đường được dự báo sẽ trở thành một trong bảy nguyên nhân gây tử vong và tàn tật hàng đầu vào năm 2030.

Trong kiểm soát đường huyết sau bữa ăn đối với bệnh nhân đái tháo đường tuýp 2, ức chế enzym α-glucosidase (enzym xúc tác quá trình thủy phân carbohydrate thành glucose tại ruột non) là một trong những đích can thiệp dược lý quan trọng. Mặc dù các hoạt chất ức chế α-glucosidase tổng hợp như acarbose hay miglitol đang được sử dụng trên lâm sàng, chúng thường gây ra nhiều tác dụng không mong muốn trên đường tiêu hóa như đầy hơi, đau bụng, tiêu chảy và buồn nôn. Do đó, xu hướng tìm kiếm các hợp chất ức chế α-glucosidase có nguồn gốc tự nhiên từ dược liệu có độ an toàn cao đang nhận được nhiều sự quan tâm.

Cây Chè vằng (Jasminum subtriplinerve Blume, họ Nhài Oleaceae) là loài thực vật phân bố phổ biến tại Việt Nam và khu vực châu Á. Dược liệu này đã được nghiên cứu về các hoạt tính sinh học như kháng khuẩn, chống oxy hóa và gây độc tế bào. Tuy nhiên, các nghiên cứu thực nghiệm về tác dụng hạ đường huyết và ức chế enzym α-glucosidase của Chè vằng tại Việt Nam còn hạn chế. Nhằm cung cấp dữ liệu khoa học phục vụ định hướng phát triển sản phẩm hỗ trợ điều trị đái tháo đường, khóa luận thực hiện hai nhiệm vụ nghiên cứu chính:

  1. Đánh giá tác dụng ức chế enzym α-glucosidase của các phân đoạn cao chiết từ lá cây Chè vằng trên mô hình thực nghiệm in vitro.
  2. Xác định các hợp chất thành phần trong cây Chè vằng có tiềm năng ức chế enzym α-glucosidase thông qua phương pháp mô phỏng gắn kết phân tử (molecular docking), đánh giá quy tắc Lipinski và dự đoán dược động học – độc tính (in silico).

Phạm vi nghiên cứu:

  • Đối tượng thực vật: Lá cây Chè vằng (Jasminum subtriplinerve Blume) thu hái tại Hà Nội vào tháng 10 năm 2021, lưu giữ tại Trường Đại học Y Dược – Đại học Quốc gia Hà Nội.
  • Dữ liệu cấu trúc hóa học: Cơ sở dữ liệu gồm 39 hợp chất đã được phân lập và xác định cấu trúc từ lá Chè vằng, tổng hợp từ các công bố khoa học trong nước trên các tạp chí quốc tế (Natural Product Research, Journal of Asian Natural Products Research).
  • Thời gian và địa điểm thực hiện: Năm 2022, tại Trường Đại học Y Dược (ĐHQGHN) và Viện Công nghệ sinh học (Viện Hàn lâm Khoa học và Công nghệ Việt Nam).

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

Khóa luận kết hợp đồng thời phương pháp thực nghiệm sinh hóa (in vitro) và các công cụ tin sinh dược học (in silico) dựa trên các khung lý thuyết sau:

  • Khung lý thuyết bệnh sinh và đích tác dụng: Áp dụng tiêu chuẩn chẩn đoán và phân loại của Hiệp hội Đái tháo đường Hoa Kỳ (ADA, 2022); cơ chế xúc tác của enzym α-glucosidase (thuộc nhóm glycoside hydrolase GH13 và GH31) trong việc phân cắt liên kết α-1,4-glycoside giải phóng phân tử α-D-glucose; cơ chế ức chế cạnh tranh của các chất ức chế α-glucosidase (AGIs).
  • Khung đánh giá độ giống thuốc và dược động học: Bộ 5 quy tắc Lipinski (Rule of 5 - RO5) nhằm sàng lọc đặc tính giống thuốc đường uống (Khối lượng phân tử MW < 500 Da; LogP < 5; số nhóm cho liên kết hydro HBD ≤ 5; số nhóm nhận liên kết hydro HBA ≤ 10; độ khúc xạ phân tử MR từ 40–130) và mô hình phân tích đặc tính ADMET (Hấp thu, Phân bố, Chuyển hóa, Thải trừ và Độc tính).

Phương pháp nghiên cứu thực nghiệm (in vitro):

  • Chiết xuất và phân đoạn: 500 g bột lá khô Chè vằng được chiết ngấm kiệt bằng dung môi ethanol (EtOH) 70% ở nhiệt độ phòng (tỷ lệ 1:10 kg/L, 3 lần × 3 ngày). Dịch chiết gộp được cô thu hồi dung môi dưới áp suất giảm thu cao toàn phần EtOH. Cao toàn phần được phân tán trong nước nóng và chiết lỏng – lỏng phân đoạn lần lượt với các dung môi có độ phân cực tăng dần: n-hexane, ethyl acetate (EtOAc), n-butanol (BuOH) và cắn nước.
  • Thử nghiệm ức chế enzym: Phản ứng ức chế enzym Yeast α-glucosidase (0,5 U/mL) sử dụng cơ chất p-nitrophenyl-α-D-glucopyranoside (pNPG) trong đệm phosphat pH 6,8. Phản ứng được dừng bằng dung dịch $\text{Na}_2\text{CO}3$ 0,2 M và đo mật độ quang (OD) tại bước sóng 405 nm trên máy đo ELISA Plate Reader (Biotek). Chất đối chứng dương là acarbose. Giá trị $\text{IC}{50}$ được tính toán bằng phần mềm TableCurve 2Dv4.

Phương pháp mô phỏng sàng lọc ảo (in silico):

  • Chuẩn bị cấu trúc protein đích: Cấu trúc tinh thể isomaltase từ nấm men Saccharomyces cerevisiae (mã PDB ID: 3A4A, độ phân giải 1,30 Å, có tính tương đồng cao với α-glucosidase) đồng kết tinh với phối tử α-D-glucopyranose. Protein được chuẩn bị bằng phần mềm MGL AutoDockTools 1.5.6 (loại nước, thêm nguyên tử hydro và điện tích Kollman). Thiết lập hộp lưới (Grid Box) kích thước 28 Å × 28 Å × 28 Å bao trùm trung tâm hoạt động tại tọa độ X = 21,286; Y = -0,781; Z = 18,638.
  • Chuẩn bị phối tử: Cấu trúc 3D của 39 hợp chất và acarbose được tải từ PubChem hoặc dựng bằng ChemOffice 19.0, tối ưu hóa năng lượng dưới trường lực MMFF94 bằng phần mềm Avogadro.
  • Quy trình docking và thẩm định mô hình: Thực hiện bằng phần mềm AutoDock Vina. Thẩm định mô hình docking thông qua phương pháp re-docking phối tử đồng tinh thể α-D-glucopyranose vào thụ thể, đối chiếu độ sai lệch bình phương trung bình gốc (RMSD).
  • Đánh giá ADMET và trực quan hóa: Sử dụng công cụ trực tuyến pkCSM để tính toán các thông số dược động học và độc tính. Phân tích tương tác không gian 2D và 3D của phức hợp protein – ligand bằng phần mềm Discovery Studio Visualizer 2020 và UCSF Chimera 1.15.

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

Chương 1. Tổng quan tài liệu

Chương 1 tổng hợp các tài liệu học thuật về 4 nhóm nội dung chính:

  • Bệnh đái tháo đường: Trình bày định nghĩa, phân loại theo ADA 2022 (tuýp 1, tuýp 2, thai kỳ, nguyên nhân khác), cơ chế bệnh sinh do tự miễn phá hủy tế bào β đảo tụy (tuýp 1) hoặc do kháng insulin kết hợp suy giảm chức năng tế bào β (tuýp 2), các biến chứng cấp/mạn tính và các nhóm thuốc điều trị hiện hành.
  • Enzym α-glucosidase và chất ức chế: Cấu trúc phân tử của họ glycoside hydrolase (GH13, GH31), cơ chế xúc tác thủy phân liên kết α-1,4-glycoside ở đầu tận không khử của đường đa/đường đôi. Cơ chế tác dụng và hạn chế của các thuốc AGI kinh điển (acarbose, miglitol), đồng thời tổng kết một số nghiên cứu chiết xuất hoạt chất ức chế α-glucosidase tự nhiên từ thực vật (Curcuma longa, Gypsophila oldhamiana, Piper umbellatum).
  • Phương pháp sàng lọc ảo và docking phân tử: Cơ sở toán học của thuật toán tìm kiếm cấu dạng và hàm tính điểm (scoring function), quy trình chuẩn bị protein/phối tử, bộ tiêu chuẩn Lipinski RO5 và các công cụ dự đoán ADMET (pkCSM).
  • Cây Chè vằng (Jasminum subtriplinerve Blume): Phân loại thực vật học, hình thái, vùng phân bố tại Việt Nam (Hòa Bình, Thái Bình, Thanh Hóa, Nghệ An, Hà Tĩnh), tổng quan thành phần hóa học (terpene glycoside, flavonoid glycoside, phenylpropanoid glycoside, triterpenoid) và các tác dụng sinh học đã được ghi nhận (kháng khuẩn, chống oxy hóa bằng DPPH, độc tế bào Hep-G2 và RD).

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

Chương này mô tả chi tiết quy trình thực nghiệm và các thông số kỹ thuật:

  • Thiết bị và hóa chất phục vụ chiết xuất dung môi phân đoạn và thử nghiệm sinh hóa trên phiến 96 giếng.
  • Công thức tính phần trăm ức chế enzym: $$% \text{ ức chế} = 100% - \left(\frac{A_{\text{mẫu thử}}}{A_{\text{đối chứng}}} \times 100\right)$$ trong đó $A = \text{OD}{\text{thí nghiệm}} - \text{OD}{\text{blank}}$.
  • Quy trình sàng lọc ảo 5 bước: (1) Chuẩn bị thụ thể 3A4A; (2) Chuẩn bị 39 phối tử; (3) Chạy docking bằng AutoDock Vina; (4) Thẩm định giá trị RMSD; (5) Xử lý năng lượng liên kết và phân tích liên kết tương tác.
  • Thiết lập ngưỡng đánh giá đặc tính dược động học theo pkCSM: độ tan, tính thấm Caco2 (log Papp), tỷ lệ hấp thu qua ruột (%), thể tích phân bố VDss, tính thấm hàng rào máu não (logBB), tương tác cytochrome P450 (CYP2D6, CYP3A4), độ thanh thải toàn phần và các thử nghiệm độc tính (AMES, độc gan, nhạy cảm da).

Chương 3. Kết quả nghiên cứu

Chương 3 trình bày các số liệu thực nghiệm và dữ liệu mô phỏng thu được:

1. Kết quả chiết xuất và phân đoạn: Từ 500 g bột lá khô Chè vằng, chiết xuất bằng EtOH 70% thu được 71,79 g cao toàn phần (hiệu suất chiết 14,4%). Sau khi phân bố lỏng – lỏng, khối lượng các phân đoạn thu được gồm:

  • Cao phân đoạn n-hexane: 1,04 g
  • Cao phân đoạn ethyl acetate: 14,66 g
  • Cao phân đoạn n-butanol: 21,78 g
  • Phân đoạn cắn nước: 29,81 g

2. Kết quả ức chế enzym α-glucosidase in vitro: Hoạt tính ức chế enzym của các phân đoạn dịch chiết được so sánh trực tiếp với chất đối chứng dương acarbose:

Mẫu nghiên cứu Nồng độ 500 µg/mL (% ức chế) Nồng độ 100 µg/mL (% ức chế) Nồng độ 20 µg/mL (% ức chế) Nồng độ 4 µg/mL (% ức chế) $\text{IC}_{50}$ (µg/mL)
Cao toàn phần EtOH 19,92 ± 1,56 2,74 ± 0,89 2,47 ± 2,89 1,93 ± 1,24 > 500
Phân đoạn n-Hexane 104,60 ± 4,72 99,75 ± 3,26 89,97 ± 2,19 16,66 ± 0,10 7,27 ± 0,71
Phân đoạn Ethyl acetate 102,19 ± 2,33 99,88 ± 1,07 89,54 ± 3,51 15,57 ± 1,07 7,42 ± 0,95
Phân đoạn n-Butanol 17,14 ± 0,37 4,24 ± 0,82 2,74 ± 0,28 1,53 ± 0,16 > 500
Phân đoạn Cắn nước 12,76 ± 1,07 0,97 ± 1,75 1,25 ± 0,19 1,73 ± 1,47 > 500
Acarbose (Chứng dương) 75,20 ± 0,62 48,57 ± 1,96 20,64 ± 0,34 8,32 ± 0,62 127,53 ± 1,73

3. Kết quả thẩm định mô hình và docking phân tử:

  • Thử nghiệm re-docking phân tử α-D-glucopyranose vào vị trí hoạt động của thụ thể 3A4A cho kết quả trùng khớp cấu dạng với sai số $\text{RMSD} = 0,5\text{ \AA}$ (< 1,0 Å), khẳng định thông số docking có độ tin cậy cao.
  • Docking 39 hợp chất của Chè vằng xác định được 8 hợp chất có năng lượng liên kết thấp nhất (ái lực gắn kết mạnh nhất), đều thấp hơn năng lượng liên kết của acarbose (-8,9 kcal/mol):
STT Tên hợp chất Năng lượng liên kết (kcal/mol) Các acid amin tạo liên kết hydro
1 Rutin -10,4 ASP352, ASP307, ARG315, THR310, GLU277
2 Nicotiflorin -10,4 ARG315, TYR158, HIS280, THR310, GLU277
3 Verbascoside -9,7 ASP307, HIS280, PRO312, GLU277, ASP352, TYR158...
4 Isoquercitrin -9,4 ASP307, ASP242, ARG315, GLU277, ASP352, ARG442...
5 Isoverbascoside -9,2 ASP215, ASP242, SER311, SER240, ASN415...
6 Stigmasterol -9,2 PRO312
7 Astragalin -9,1 THR310, ASP352, GLN353, TYR158, ARG315, PRO312
8 Chevangin B (3) -9,1 ARG315, ASP242, GLU277, TYR158, ASP307, ASP352

4. Kết quả lọc quy tắc Lipinski và dự đoán ADMET: Trong 8 hợp chất trên, có 4 hợp chất thỏa mãn các tiêu chí của quy tắc Lipinski để tiến hành đánh giá ADMET: Astragalin, Isoquercitrin, Verbascoside và Stigmasterol.

Thông số dược động học & độc tính Astragalin Isoquercitrin Verbascoside Stigmasterol
Độ tan trong nước (log mol/L) -2,883 -2,883 -2,878 -6,045
Tính thấm Caco2 (log Papp, $10^{-6}\text{ cm/s}$) 0,306 0,242 0,096 1,213
Hấp thu qua đường ruột (%) 48,052 47,842 41,178 94,970
Tính thấm hàng rào máu não (logBB) -1,374 -1,452 -1,863 0,447
Cơ chất CYP2D6 / CYP3A4 Không / Không Không / Không Không / Không Không / Có
Ức chế CYP2D6 / CYP3A4 Không / Không Không / Không Không / Không Không / Không
Độ thanh thải toàn phần (log mL/min/kg) 0,634 0,554 0,618 0,937
Độc tính AMES / Độc tính gan Không / Không Không / Không Không / Không Không / Không

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

Chương 4 tập trung phân tích ý nghĩa của các phát hiện:

  • Về hoạt tính ức chế in vitro: Phân đoạn n-hexane ($\text{IC}{50} = 7,27\text{ }\mu\text{g/mL}$) và ethyl acetate ($\text{IC}{50} = 7,42\text{ }\mu\text{g/mL}$) thể hiện khả năng ức chế α-glucosidase mạnh gấp hơn 17 lần so với chất đối chứng acarbose ($\text{IC}_{50} = 127,53\text{ }\mu\text{g/mL}$). Điều này chứng minh các hoạt chất có hoạt tính tập trung chủ yếu ở nhóm phân cực trung bình và kém phân cực trong lá Chè vằng.
  • Về mô hình sàng lọc in silico và tương tác phân tử: 4 hợp chất tiềm năng (Astragalin, Isoquercitrin, Verbascoside, Stigmasterol) tương tác chặt chẽ với các acid amin xúc tác quan trọng tại trung tâm hoạt động của enzym như GLU277, ASP352, ARG442, ARG446, GLN279 và HIS280. Trong đó, verbascoside và stigmasterol chia sẻ lần lượt 15 và 17 acid amin tương tác chung với thuốc chuẩn acarbose.
  • Về đặc tính ADMET: Stigmasterol có tỷ lệ hấp thu đường ruột cao (94,97%), trong khi các flavonoid glycoside (astragalin, isoquercitrin) và verbascoside có mức hấp thu trung bình (41–48%). Cả 4 chất đều không thể hiện độc tính gây đột biến gen (AMES) và không gây độc tế bào gan trên mô hình dự đoán.

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

  • Dữ liệu thực nghiệm mới về hoạt tính sinh học: Khóa luận chứng minh dịch chiết phân đoạn n-hexane và ethyl acetate từ lá cây Chè vằng có tác dụng ức chế mạnh enzym α-glucosidase in vitro với $\text{IC}_{50}$ lần lượt là $7,27 \pm 0,71\text{ }\mu\text{g/mL}$ và $7,42 \pm 0,95\text{ }\mu\text{g/mL}$, vượt trội so với chất đối chứng acarbose ($127,53 \pm 1,73\text{ }\mu\text{g/mL}$).
  • Xây dựng cơ sở dữ liệu và xác định các phân tử tiềm năng: Tổng hợp cơ sở dữ liệu cấu trúc của 39 hợp chất có trong cây Chè vằng, sàng lọc thành công 4 hợp chất có tiềm năng phát triển thuốc/sản phẩm hỗ trợ điều trị đái tháo đường gồm: Astragalin, Isoquercitrin, Verbascoside và Stigmasterol. Cả 4 chất đều đáp ứng quy tắc Lipinski, có năng lượng liên kết với thụ thể thấp hơn acarbose và sở hữu hồ sơ an toàn ADMET thuận lợi.
  • Đóng góp học thuật: Cung cấp cơ sở khoa học hiện đại kết hợp giữa sinh hóa thực nghiệm và hóa tin học cho việc sử dụng Chè vằng trong y học cổ truyền, đồng thời định hướng cho việc phân lập các hoạt chất hạ đường huyết từ hai phân đoạn dịch chiết tiềm năng.

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

  • Hạn chế: Nghiên cứu thực nghiệm in vitro mới thực hiện trên các phân đoạn cao chiết thô và bán tinh khiết mà chưa tiến hành phân lập riêng lẻ từng đơn chất tinh khiết trong phân đoạn n-hexane và ethyl acetate để thử hoạt tính trực tiếp. Ngoài ra, kết quả mô phỏng in silico là dự đoán trên mô hình tính toán máy tính, cần kiểm chứng thực nghiệm sinh học chuyên sâu.
  • Hướng nghiên cứu tiếp được đề xuất:
    1. Tiến hành phân lập và xác định cấu trúc hóa học chi tiết của các đơn chất có trong phân đoạn n-hexane và ethyl acetate từ lá cây Chè vằng nhằm kiểm chứng hoạt tính ức chế enzym α-glucosidase của từng chất.
    2. Đánh giá tác dụng hạ glucose huyết và cơ chế tác dụng của các phân đoạn n-hexane và ethyl acetate trên mô hình động vật thực nghiệm (in vivo).

Giá trị tham khảo

  • Đối tượng sử dụng: Khóa luận là tài liệu tham khảo cho sinh viên, học viên cao học, nghiên cứu sinh và các nhà nghiên cứu thuộc các chuyên ngành Dược học, Dược liệu học, Dược lý học, Hóa dược và Công nghệ sinh học.
  • Nội dung giá trị: Quy trình chiết xuất phân đoạn dung môi dược liệu; phương pháp thực nghiệm đánh giá ức chế enzym α-glucosidase sử dụng cơ chất pNPG trên phiến 96 giếng; quy trình tích hợp mô phỏng sàng lọc ảo AutoDock Vina, thẩm định mô hình bằng chỉ số RMSD và phân tích dược động học/độc tính bằng pkCSM.

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

1. Phân đoạn chiết nào từ lá cây Chè vằng thể hiện tác dụng ức chế α-glucosidase mạnh nhất?
Phân đoạn n-hexane và phân đoạn ethyl acetate thể hiện hoạt tính ức chế mạnh nhất với giá trị $\text{IC}{50}$ lần lượt là $7,27 \pm 0,71\text{ }\mu\text{g/mL}$ và $7,42 \pm 0,95\text{ }\mu\text{g/mL}$, trong khi cao toàn phần ethanol, phân đoạn n-butanol và phân đoạn nước đều có $\text{IC}{50} > 500\text{ }\mu\text{g/mL}$.

2. Mô hình docking phân tử sử dụng cấu trúc protein đích nào và tiêu chuẩn thẩm định ra sao?
Nghiên cứu sử dụng cấu trúc tinh thể isomaltase từ nấm men Saccharomyces cerevisiae (PDB ID: 3A4A, độ phân giải 1,30 Å). Mô hình được thẩm định bằng cách re-docking phân tử α-D-glucopyranose vào vị trí hoạt động; kết quả đạt giá trị $\text{RMSD} = 0,5\text{ \AA}$ (đạt chuẩn cho phép < 1,0 Å), chứng minh phương pháp tính toán có độ chính xác cao.

3. Trong số 39 hợp chất của Chè vằng, những chất nào có ái lực liên kết mạnh nhất và chất nào được chọn sau sàng lọc?
Có 8 hợp chất đạt năng lượng liên kết thấp nhất từ -9,1 đến -10,4 kcal/mol (mạnh hơn acarbose ở mức -8,9 kcal/mol) gồm: Rutin, Nicotiflorin, Verbascoside, Isoquercitrin, Isoverbascoside, Stigmasterol, Astragalin và Chevangin B. Sau khi sàng lọc qua quy tắc 5 tiêu chí của Lipinski, 4 hợp chất thỏa mãn điều kiện giống thuốc được lựa chọn là Astragalin, Isoquercitrin, Verbascoside và Stigmasterol.

4. Kết quả dự đoán dược động học và độc tính (ADMET) của 4 hợp chất tiềm năng cho thấy điều gì?
Cả 4 hợp chất (Astragalin, Isoquercitrin, Verbascoside, Stigmasterol) đều không có độc tính gây đột biến gen (thử nghiệm AMES âm tính), không gây độc cho gan và không gây kích ứng da. Về hấp thu, Stigmasterol có khả năng hấp thu qua ruột cao nhất (94,97%), trong khi 3 chất còn lại đạt mức từ 41,18% đến 48,05%.


Kết luận

Khóa luận tốt nghiệp của tác giả Lê Minh Ngọc đã hoàn thành việc đánh giá hoạt tính ức chế enzym α-glucosidase của cây Chè vằng (Jasminum subtriplinerve Blume) thông qua sự kết hợp giữa mô hình thực nghiệm sinh hóa in vitro và mô phỏng gắn kết phân tử in silico. Kết quả thực nghiệm xác định hai phân đoạn chiết n-hexane và ethyl acetate có hoạt tính ức chế α-glucosidase mạnh vượt trội so với thuốc chuẩn acarbose. Đồng thời, mô hình sàng lọc ảo từ 39 hợp chất đã xác định được 4 phân tử tiềm năng (Astragalin, Isoquercitrin, Verbascoside, Stigmasterol) thỏa mãn các tiêu chí giống thuốc Lipinski và an toàn về mặt dược động học. Công trình cung cấp cơ sở khoa học thực chứng quan trọng cho việc tiếp tục nghiên cứu phân lập hoạt chất và thử nghiệm in vivo trên mô hình điều trị đái tháo đường từ nguồn dược liệu Chè vằng.