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

Hội chứng chuyển hóa (Metabolic syndrome - METS) là một tập hợp các rối loạn chuyển hóa phức tạp, bao gồm tình trạng béo bụng, tăng huyết áp, tăng đường huyết lúc đói, kháng insulin hoặc không dung nạp đường, tình trạng tiền đông máu (tăng nồng độ fibrinogen và chất ức chế plasminogen hoạt hóa PAI-1 trong máu), và tình trạng tiền viêm (tăng protein phản ứng C - CRP). METS làm gia tăng trực tiếp nguy cơ mắc bệnh đái tháo đường type 2, bệnh tim mạch do xơ vữa động mạch và đột quỵ. Theo thống kê của Liên đoàn Đái tháo đường Quốc tế (IDF), có khoảng 25% dân số toàn cầu mắc hội chứng này. Tại Việt Nam, nghiên cứu của Viện Dinh dưỡng (2000) trên 620 đối tượng từ 25-64 tuổi ghi nhận tỷ lệ mắc METS là 13%, riêng tại khu vực Hà Nội và Thành phố Hồ Chí Minh tỷ lệ này lên đến 18%.

Địa liền đen (Kaempferia parviflora Wall. ex Baker, họ Gừng - Zingiberaceae; trong tiêu đề bìa khóa luận ghi Kaempferia galanga L.) là loài cây thuốc phân bố phổ biến ở Đông Nam Á, tại Việt Nam được phát hiện ở các tỉnh Đắk Lắk, Gia Lai và Thanh Hóa. Thân rễ của cây được sử dụng trong y học dân gian để hỗ trợ điều trị béo phì, đái tháo đường và cải thiện sức khỏe chuyển hóa nhờ chứa hàm lượng cao các hợp chất polymethoxyflavon và tinh dầu. Tuy nhiên, cơ chế phân tử và mối liên hệ tổng thể giữa các thành phần hóa học của Địa liền đen với các đích sinh học điều trị hội chứng chuyển hóa chưa được làm sáng tỏ một cách có hệ thống. Việc ứng dụng các công cụ tin sinh học và dược lý mạng (network pharmacology) kết hợp mô phỏng tương tác phân tử (molecular docking) mang lại hướng tiếp cận toàn diện nhằm giải mã cơ chế đa thành phần - đa mục tiêu của dược liệu này.

Khóa luận xác định các mục tiêu và nhiệm vụ nghiên cứu cụ thể bao gồm:

  1. Xây dựng mạng lưới hợp chất - đích tác dụng tiềm năng và mạng lưới tương tác protein - protein (PPI), từ đó dự đoán các hợp chất và các đích tác dụng tiềm năng trên hội chứng chuyển hóa của Địa liền đen.
  2. Dự đoán các quá trình và con đường sinh học liên quan đến các đích tác dụng tiềm năng trên hội chứng chuyển hóa của Địa liền đen thông qua phân tích bản thể học gen (Gene Ontology - GO) và phân tích làm giàu con đường KEGG (Kyoto Encyclopedia of Genes and Genomes).
  3. Mô phỏng tương tác phân tử (molecular docking) của các hợp chất tiềm năng có trong Địa liền đen với các protein thụ thể chủ chốt của hội chứng chuyển hóa.

Đối tượng nghiên cứu của khóa luận là các thành phần hóa học (nhóm polymethoxyflavon và tinh dầu) phân lập hoặc xác định từ thân rễ cây Địa liền đen cùng hệ thống các gen và protein đích liên quan đến bệnh sinh hội chứng chuyển hóa ở người (Homo sapiens). Phạm vi nghiên cứu được thực hiện hoàn toàn trên không gian dữ liệu số hóa (in silico) thông qua các cơ sở dữ liệu y sinh học quốc tế (GenCards, PubChem, SwissTargetPrediction, STRING, PDB) và phần mềm tính toán chuyên dụng, hoàn thành năm 2024 tại Bộ môn Dược liệu - Trường Đại học Dược 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 khung lý thuyết của Dược lý mạng (ra đời từ những năm 20 của thế kỷ XXI), kết hợp giữa sinh học phân tử, hóa sinh, tin sinh học và lý thuyết đồ thị mạng lưới để chuyển đổi phương thức nghiên cứu dược liệu từ thực nghiệm đơn lẻ sang nghiên cứu hệ thống dựa trên bằng chứng dữ liệu lớn. Các thông số hình học không gian của mạng lưới được sử dụng để đánh giá bao gồm nút (node đại diện cho gen, protein, hợp chất), cạnh (edge đại diện cho tương tác), bậc (degree - DC), và bậc trung tâm trung gian (betweenness centrality - BC) để xác định các nút trung tâm (hub nodes).

Khung phân tích chức năng sinh học sử dụng hệ thống Gene Ontology (GO) chia theo 3 khía cạnh: Quá trình sinh học (BP), Thành phần tế bào (CC), và Chức năng phân tử (MF); kết hợp với cơ sở dữ liệu làm giàu con đường tín hiệu KEGG (Kyoto Encyclopedia of Genes and Genomes). Đánh giá tương tác phân tử dựa trên phương pháp docking phân tử với ngưỡng năng lượng liên kết $\le -7$ kcal/mol để xác định ái lực gắn kết mạnh và bền vững giữa phối tử (ligand) và thụ thể protein (receptor).

Giai đoạn nghiên cứu Công cụ / Phần mềm / Cơ sở dữ liệu Tiêu chí và tham số xử lý dữ liệu
Thu thập gen bệnh METS GeneCards (genecards.org) Từ khóa "metabolic syndrome", loài Homo sapiens
Thu thập thành phần hóa học PubMed, Google Scholar, PubChem, ChemDraw Thu thập danh mục hợp chất và chuỗi định dạng SMILES
Sàng lọc sinh khả dụng ADMETlab 2.0 Xác suất sinh khả dụng đường uống $F(30%) \le 0.3$
Dự đoán đích tác dụng SwissTargetPrediction (swisstargetprediction.ch) Chỉ số xác suất $\text{Probability} > 0$, loài Homo sapiens (tối đa 100 đích/chất)
Xác định đích tiềm năng Công cụ Venny 2.1 Giao thoa giữa tập hợp đích của dược liệu và gen bệnh METS
Mạng lưới tương tác PPI STRING database v11.5 (string-db.org) Đánh giá tương tác protein người, độ tin cậy tương tác
Phân tích & trực quan mạng Cytoscape 3.0 (tích hợp CytoNCA) Tính toán hệ số topo mạng lưới: Degree Centrality (DC), Betweenness Centrality (BC)
Phân tích GO và KEGG ShinyGO v0.75 Đánh giá giá trị sai số FDR (False Discovery Rate) và độ làm giàu Fold Enrichment
Docking phân tử Ngân hàng dữ liệu Protein (RCSB PDB) Cấu trúc 3D protein chuẩn: 7D9O, 3O9M, 5IKR, 2GIU, 3S79; ngưỡng năng lượng $\le -7$ kcal/mol

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

Chương 1: Tổng quan

Chương 1 tổng thuật các kiến thức đại cương về hội chứng chuyển hóa, đặc điểm thực vật, thành phần hóa học, tác dụng sinh học của Địa liền đen cùng nền tảng lý thuyết các phương pháp in silico:

  • Đại cương về Hội chứng chuyển hóa: Tổng hợp các định nghĩa, yếu tố nguy cơ (tuổi tác, chủng tộc, béo phì với BMI > 23 và phân bố mỡ dạng quả táo, tiền sử gia đình đái tháo đường, hội chứng buồng trứng đa nang) và 5 bộ tiêu chuẩn chẩn đoán phổ biến (NCEP, WHO, IDF, AACE, EGIR). Trong đó tiêu chuẩn của WHO xác định khi có từ 3 yếu tố: vòng bụng $\ge 90$ cm (nam) hoặc $\ge 80$ cm (nữ); triglycerid $\ge 150$ mg/dl; HDL-C $< 40$ mg/dl (nam) hoặc $< 50$ mg/dl (nữ); huyết áp $\ge 130/85$ mmHg; glucose máu lúc đói $\ge 100$ mg/dl. Tỷ lệ rối loạn lipid máu ở người trưởng thành Việt Nam có tăng triglycerid từ 15,8% đến 46,6% và HDL-C thấp từ 15,7% đến 52,3%.
  • Thành phần hóa học và tác dụng của Địa liền đen: Thống kê 16 hợp chất thuộc nhóm polymethoxyflavon (như 3,5,7,3',4'-pentamethoxyflavon, 5,7,4'-trimethoxyflavon, 3,5,7-trimethoxyflavon...) và 28 thành phần tinh dầu (như $\alpha$-pinen, camphen, eucalyptol, linalool, bornyl acetat...). Tổng quan các nghiên cứu thực nghiệm chứng minh tác dụng hạ glucose máu (trên chuột béo phì do HFD), chống béo phì (hoạt hóa ATGL và HSL độc lập PPAR$\gamma$, thử nghiệm lâm sàng cho thấy 12 mg/ngày polymethoxyflavon làm giảm mỡ nội tạng ở người có BMI từ 23–30 $\text{kg/m}^2$), giãn mạch hạ huyết áp (qua con đường eNOS và NO), chống viêm/giảm stress oxy hóa và tăng cường chuyển hóa năng lượng cơ bắp thông qua hoạt hóa AMPK, tăng biểu hiện GLUT4 và MCT1.
  • Nền tảng phương pháp luận: Giới thiệu các khái niệm, quy trình 3 giai đoạn của dược lý mạng (xây dựng mạng, phân tích mạng, thẩm định), lý thuyết phân tích bản thể gen GO, phân tích con đường KEGG và mô phỏng tương tác phân tử docking.

Chương 2: Nguyên liệu và phương pháp nghiên cứu

Chương 2 trình bày chi tiết trang thiết bị (máy tính Macbook Pro), phần mềm (Cytoscape 3.0, ShinyGO v0.75, Venny 2.1) cùng quy trình thực hiện 6 nội dung nghiên cứu theo sơ đồ thiết kế:

  1. Tìm kiếm và thu thập gen liên quan đến METS trên cơ sở dữ liệu GenCards.
  2. Thu thập thành phần hóa học Địa liền đen qua PubMed, Google Scholar, PubChem, sau đó sàng lọc sinh khả dụng đường uống bằng ADMETlab 2.0 ($F(30%) \le 0.3$), dự đoán đích tác dụng tương ứng qua SwissTargetPrediction.
  3. Xác định tập hợp đích tác dụng tiềm năng chung bằng biểu đồ Venny và xây dựng mạng lưới "Hợp chất - Đích tác dụng tiềm năng".
  4. Thiết lập mạng lưới tương tác protein - protein (PPI) trên STRING v11.5 và tính toán thông số mạng bằng CytoNCA trên phần mềm Cytoscape.
  5. Thực hiện phân tích GO (quá trình sinh học, thành phần tế bào, chức năng phân tử) và con đường KEGG bằng ShinyGO v0.75 trên loài Homo sapiens.
  6. Mô phỏng docking phân tử giữa các hợp chất then chốt và các thụ thể protein đại diện cho các đích tác dụng trung tâm.

Chương 3: Thực nghiệm, kết quả và bàn luận

Chương 3 công bố toàn bộ dữ liệu tính toán, phân tích topo mạng lưới và đánh giá cơ chế tác dụng:

  • Thu thập gen bệnh và thành phần hóa học: Thu được 2243 gen liên quan đến METS từ GenCards. Tổng hợp được 44 hợp chất từ thân rễ Địa liền đen, qua sàng lọc thông số sinh khả dụng đường uống $F(30%) \le 0.3$ đã chọn ra 26 hợp chất đạt tiêu chuẩn (gồm 17 hợp chất tinh dầu/terpenoid và 9 hợp chất polymethoxyflavon).
  • Xác định đích tác dụng tiềm năng: 26 hợp chất đạt sinh khả dụng được dự đoán có 252 đích tác dụng sinh học không trùng lặp. Phân tích giao thoa giữa 252 đích này và 2243 gen bệnh METS thu được 248 đích tác dụng tiềm năng chung.
  • Mạng lưới Hợp chất - Đích tác dụng: Có 22 hợp chất có đích liên quan đến METS. Mạng lưới cô đọng xác định 16 nút cốt lõi gồm 11 hợp chất tiềm năng và 5 đích tác dụng trung tâm với 55 cạnh liên kết. Trong đó, 7 hợp chất tác động đến nhiều mục tiêu nhất đều thuộc nhóm polymethoxyflavon: 3,5,7-trimethoxyflavon tác động đến 108 đích; 6 hợp chất gồm 5-hydroxy-3,7,3',4'-tetramethoxyflavon, 5-hydroxy-3,7,4'-trimethoxyflavon, 5-hydroxy-3,7-dimethoxyflavon, 3,5-dihydroxy-7,3',4'-trimethoxyflavon, 5,3'-dihydroxy-3,7,4'-trimethoxyflavon và 3,5,7,3',4'-pentamethoxyflavon mỗi chất tác động đến 104 đích.
STT Tên hợp chất quan trọng Nhóm cấu trúc Số đích tác dụng trong METS
1 3,5,7-trimethoxyflavon Polymethoxyflavon 108
2 5-hydroxy-3,7,3',4'-tetramethoxyflavon Polymethoxyflavon 104
3 5-hydroxy-3,7,4'-trimethoxyflavon Polymethoxyflavon 104
4 5-hydroxy-3,7-dimethoxyflavon Polymethoxyflavon 104
5 3,5-dihydroxy-7,3',4'-trimethoxyflavon Polymethoxyflavon 104
6 5,3'-dihydroxy-3,7,4'-trimethoxyflavon Polymethoxyflavon 104
7 3,5,7,3',4'-pentamethoxyflavon Polymethoxyflavon 104
8 5,7,3',4'-tetramethoxyflavon Polymethoxyflavon -
9 L-$\alpha$-Terpineol Tinh dầu / Terpenoid -
10 Bornyl acetat Tinh dầu / Terpenoid -
11 $\alpha$-Cadinol Tinh dầu / Terpenoid -
  • Mạng lưới PPI và 5 đích trung tâm: Phân tích mạng lưới tương tác protein - protein của 248 đích bằng CytoNCA xác định PTGS2 là nút trung tâm quan trọng nhất với hệ số BC và DC cao nhất. Năm đích trung tâm hàng đầu được chọn lọc gồm: PTGS2, CYP19A1, ACHE, ESR2, và BCHE.
  • Phân tích GO và KEGG: Xác định 20 quá trình sinh học, 20 thành phần tế bào, 20 chức năng phân tử và 20 con đường KEGG có ý nghĩa thống kê cao nhất. Các con đường chuyển hóa tiêu biểu gồm: con đường dẫn truyền tín hiệu PI3K-Akt, con đường EGFR kháng thuốc ức chế tyrosin kinase, sinh tổng hợp hormon steroid, các quá trình trao đổi chất và đề kháng insulin.
  • Mô phỏng Docking phân tử: Thực hiện mô phỏng tương tác giữa 11 hợp chất quan trọng với 5 thụ thể protein trung tâm từ cơ sở dữ liệu PDB.
Đích tác dụng Tên Protein thụ thể Mã PDB cấu trúc 3D Đường dẫn nguồn dữ liệu
CYP19A1 Aromatase / Cytochrome P450 19A1 3S79 rcsb.org/structure/3S79
ACHE Acetylcholinesterase 7D9O rcsb.org/structure/7D9O
PTGS2 Prostaglandin-endoperoxide synthase 2 (COX-2) 5IKR rcsb.org/structure/5IKR
ESR2 Estrogen receptor 2 2GIU rcsb.org/structure/2GIU
BCHE Butyrylcholinesterase 3O9M rcsb.org/structure/3O9M

Kết quả docking cho thấy toàn bộ 11 hợp chất đều có năng lượng liên kết với các protein đạt ngưỡng $\le -7$ kcal/mol (riêng L-$\alpha$-Terpineol có giá trị $-6.3$ kcal/mol trên một số protein). Đích ACHE có số lượng hợp chất tương tác nhiều nhất và lực liên kết mạnh nhất, tiếp theo là BCHE, CYP19A1, PTGS2 và ESR2. Hai hợp chất có năng lượng liên kết $\le -7$ kcal/mol trên toàn bộ 5 đích thụ thể là 5-hydroxy-3,7-dimethoxyflavon và 3,5,7-trimethoxyflavon.

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

  • Hệ thống hóa dữ liệu và sàng lọc hoạt chất: Xác định được 26/44 hợp chất trong thân rễ Địa liền đen đạt yêu cầu về sinh khả dụng đường uống ($F(30%) \le 0.3$), dự đoán 248 đích tác dụng tiềm năng can thiệp vào cơ chế bệnh sinh của hội chứng chuyển hóa.
  • Làm rõ 11 hợp chất và 5 đích protein trung tâm: Thu hẹp và chứng minh vai trò chủ đạo của 11 hợp chất hoạt tính cao (tiêu biểu là 8 polymethoxyflavon và 3 thành phần tinh dầu) tác động trực tiếp lên 5 đích enzyme/thụ thể trung tâm gồm ACHE, BCHE, PTGS2, ESR2 và CYP19A1.
  • Khám phá hai phân tử tiềm năng nhất: Qua docking phân tử, xác định 5-hydroxy-3,7-dimethoxyflavon và 3,5,7-trimethoxyflavon là hai phân tử có khả năng gắn kết mạnh mẽ và đồng thời lên toàn bộ 5 đích thụ thể chuyển hóa với mức năng lượng liên kết $\le -7$ kcal/mol.
  • Đóng góp khoa học: Đây là công trình đầu tiên ứng dụng tích hợp phương pháp Dược lý mạng và Docking phân tử trên Địa liền đen để giải thích cơ chế tác dụng hiệp đồng đa thành phần, đa mục tiêu, đa con đường (PI3K-Akt, chuyển hóa lipid, ức chế enzyme cholinesterase và kháng viêm) trong điều trị hội chứng chuyển hóa.
  • Đề xuất kiến nghị: Tác giả kiến nghị tiến hành các thử nghiệm thực nghiệm sinh học in vitro và in vivo tiếp theo nhằm thẩm định tính chính xác của các kết quả dự đoán in silico.

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

  • Hạn chế: Toàn bộ dữ liệu của khóa luận được triển khai trên các mô hình tính toán và thuật toán mô phỏng in silico dựa vào cơ sở dữ liệu mở; chưa tiến hành thẩm định thực nghiệm sinh học in vitro (đo hoạt tính ức chế enzyme AChE, BChE, COX-2) hoặc in vivo (trên mô hình động vật mắc hội chứng chuyển hóa) trong khuôn khổ khóa luận.
  • Hướng nghiên cứu tiếp: Thực hiện chiết xuất, phân lập hoặc tổng hợp 2 hợp chất tiềm năng (5-hydroxy-3,7-dimethoxyflavon và 3,5,7-trimethoxyflavon) để tiến hành các thử nghiệm dược lý thực nghiệm, kiểm chứng ái lực liên kết thực tế trên các dòng tế bào và mô hình động vật tiền lâm sàng.

Giá trị tham khảo

  • Đối tượng thụ hưởng: Tài liệu hữu ích cho sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Dược học, Dược liệu học, Hóa dược, Tin sinh học và Dược học tính toán.
  • Nội dung tham khảo giá trị:
    • Quy trình từng bước chuẩn hóa trong nghiên cứu Dược lý mạng kết hợp Docking phân tử đối với cây thuốc có thành phần hóa học phức tạp.
    • Bộ cơ sở dữ liệu sàng lọc về 26 hợp chất đạt sinh khả dụng đường uống và 248 đích tác dụng liên quan đến hội chứng chuyển hóa của Địa liền đen.
    • Các tham số cấu trúc phân tử, mã PDB protein thụ thể và giá trị năng lượng liên kết làm tiền đề cho việc thiết kế, phát triển các chế phẩm hoặc dẫn chất mới hỗ trợ điều trị hội chứng chuyển hóa.

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

1. Khóa luận sử dụng tiêu chuẩn và công cụ nào để sàng lọc các hợp chất có sinh khả dụng đường uống đạt yêu cầu?
Nghiên cứu sử dụng công cụ ADMETlab 2.0 để tính toán thông số xác suất sinh khả dụng đường uống $F(30%)$ từ định dạng chuỗi SMILES của hợp chất. Các hợp chất có giá trị $F(30%) \le 0.3$ (tương ứng với xác suất sinh khả dụng đường uống $\ge 30%$) được lựa chọn; kết quả thu được 26 hợp chất đạt yêu cầu từ 44 hợp chất ban đầu.

2. Năm đích tác dụng trung tâm của Địa liền đen trên hội chứng chuyển hóa được xác định qua phân tích mạng lưới PPI là những đích nào?
Năm đích trung tâm có hệ số bậc trung tâm trực tiếp (DC) và bậc trung tâm trung gian (BC) cao nhất trong mạng lưới tương tác protein - protein gồm có: PTGS2 (Prostaglandin-endoperoxide synthase 2 / COX-2, là nút có BC và DC cao nhất), CYP19A1 (Aromatase), ACHE (Acetylcholinesterase), ESR2 (Estrogen receptor 2), và BCHE (Butyrylcholinesterase).

3. Mã PDB của các cấu trúc 3D protein thụ thể được sử dụng trong mô phỏng Docking phân tử là gì?
Cấu trúc 3D của 5 protein thụ thể được tải từ Ngân hàng Dữ liệu Protein (RCSB PDB) bao gồm: ACHE (mã PDB: 7D9O), BCHE (mã PDB: 3O9M), PTGS2 (mã PDB: 5IKR), ESR2 (mã PDB: 2GIU), và CYP19A1 (mã PDB: 3S79).

4. Hai hợp chất nào trong Địa liền đen được đánh giá là tiềm năng nhất qua mô phỏng tương tác phân tử?
Hai hợp chất tiềm năng nhất được lựa chọn là 5-hydroxy-3,7-dimethoxyflavon và 3,5,7-trimethoxyflavon (đều thuộc nhóm polymethoxyflavon). Cả hai hợp chất này đều đạt mức năng lượng liên kết $\le -7$ kcal/mol với toàn bộ 5 đích protein thụ thể nghiên cứu, thể hiện ái lực tương tác mạnh và đồng thời trên nhiều cơ chế chuyển hóa.

5. Những con đường tín hiệu và chuyển hóa nào nổi bật nhất qua phân tích KEGG?
Phân tích KEGG xác định 20 con đường làm giàu có ý nghĩa, trong đó các con đường quan trọng liên quan trực tiếp đến cơ chế bệnh sinh của hội chứng chuyển hóa bao gồm: con đường dẫn truyền tín hiệu PI3K-Akt, con đường EGFR kháng thuốc ức chế tyrosin kinase, sinh tổng hợp hormon steroid, các quá trình chuyển hóa nội sinh và điều hòa kháng insulin.

Kết luận

Khóa luận tốt nghiệp của tác giả Ngô Hằng Nga đã ứng dụng thành công phương pháp tiếp cận in silico tích hợp Dược lý mạng, phân tích làm giàu GO/KEGG và Docking phân tử để làm sáng tỏ mối liên quan giữa các thành phần hóa học của Địa liền đen (Kaempferia parviflora) và tác dụng trên hội chứng chuyển hóa. Nghiên cứu đã xác định được 248 đích tác dụng tiềm năng, làm nổi bật vai trò trung tâm của 11 hợp chất và 5 đích protein (PTGS2, CYP19A1, ACHE, ESR2, BCHE), đồng thời phát hiện 5-hydroxy-3,7-dimethoxyflavon cùng 3,5,7-trimethoxyflavon là hai phân tử có tiềm năng tác dụng hiệp đồng đa mục tiêu mạnh nhất. Kết quả cung cấp cơ sở dữ liệu khoa học tin cậy và định hướng then chốt cho các nghiên cứu thực nghiệm in vitro, in vivo tiếp theo trong phát triển dược phẩm từ Địa liền đen.