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

Bệnh béo phì đang trở thành một thách thức y tế công cộng mang tính toàn cầu do tốc độ gia tăng nhanh chóng cùng các biến chứng chuyển hóa nghiêm trọng. Theo số liệu của Tổ chức Y tế Thế giới (WHO), tỉ lệ béo phì toàn cầu đã tăng gần 3 lần trong giai đoạn từ năm 1975 đến 2016. Tính đến năm 2016, toàn cầu có hơn 1,9 tỷ người trưởng thành (tương đương 39% dân số từ 18 tuổi trở lên) bị thừa cân, trong đó có hơn 650 triệu người (chiếm 13%) mắc bệnh béo phì. Tình trạng này cũng gia tăng nhanh ở nhóm trẻ em và thanh thiếu niên: năm 2016 có khoảng 340 triệu trẻ từ 5 đến 19 tuổi thừa cân hoặc béo phì (tăng từ 4% năm 1975 lên hơn 18% năm 2016), và đến năm 2020 có khoảng 39 triệu trẻ em dưới 5 tuổi gặp tình trạng này. Sự tích tụ mô mỡ bất thường bắt nguồn từ sự mất cân bằng kéo dài giữa năng lượng hấp thu và năng lượng tiêu hao, dẫn đến hàng loạt bệnh lý mạn tính đi kèm như kháng insulin, đái tháo đường týp 2, tăng huyết áp, gan nhiễm mỡ liên quan rối loạn chuyển hóa (viêm gan nhiễm mỡ không do rượu), bệnh tim mạch và một số bệnh ung thư.

Trong cơ chế sinh lý bệnh, phần lớn triglyceride chuỗi dài trong chế độ ăn được nhũ tương hóa và thủy phân bởi enzym lipase tuyến tụy (Human Pancreatic Lipase - HPL), chịu trách nhiệm phân cắt 50–70% lượng chất béo trong lòng ruột non thành monoglyceride và axit béo tự do trước khi hấp thu. Do đó, ức chế HPL là một hướng tiếp cận dược lý quan trọng giúp giảm hấp thu năng lượng tại đường tiêu hóa. Tuy nhiên, tính đến thời điểm tác giả thực hiện nghiên cứu, Orlistat là thuốc không kê đơn duy nhất được Cục Quản lý Thực phẩm và Dược phẩm Hoa Kỳ (FDA) phê duyệt tác động theo cơ chế ức chế lipase tuyến tụy. Dù có nguồn gốc bán tổng hợp từ lipstatin tự nhiên, Orlistat 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, phân mỡ, đi ngoài không tự chủ, cùng các nguy cơ sỏi thận, nhiễm độc thận và viêm tụy. Nhóm hợp chất coumarin (1-benzopyran) là các chất chuyển hóa thứ cấp phong phú trong tự nhiên với hơn 1.300 dẫn xuất phân bố ở gần 30 họ thực vật, thể hiện phổ hoạt tính sinh học rộng và có độ an toàn tương đối cao, mở ra tiềm năng tìm kiếm các chất ức chế HPL mới thay thế.

Khóa luận tốt nghiệp ngành Dược học của sinh viên Nguyễn Thị Hồng Hạnh (Trường Đại học Y Dược – Đại học Quốc gia Hà Nội, năm 2023) xác định 2 mục tiêu nghiên cứu cụ thể:

  1. Sàng lọc các hợp chất Coumarins có tác dụng ức chế enzym lipase tuyến tụy 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 là 50 hợp chất thuộc khung cấu trúc coumarin tự nhiên được phân lập hoặc định danh từ cơ sở dữ liệu hóa học. Thụ thể đích sinh học là cấu trúc tinh thể 3D của enzym lipase tuyến tụy người phức hợp với chất ức chế C11 alkyl phosphonate. 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 sinh học điện toán (in silico), sử dụng các công cụ phần mềm chuyên dụng và cơ sở dữ liệu cấu trúc trực tuyến, hoàn thành vào tháng 5 năm 2023 tại Hà Nội.

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

Khóa luận áp dụng phương pháp thiết kế thuốc dựa trên cấu trúc (Structure-Based Drug Design - SBDD) thuộc nhánh thiết kế thuốc có sự hỗ trợ của máy tính (Computer-Aided Drug Design - CADD). Cơ sở lý thuyết của mô hình bao gồm:

  • Cấu trúc sinh học mục tiêu: Lipase tuyến tụy người (HPL) là một serine esterase có bộ ba xúc tác Ser152 - Asp176 - His263 nằm ở miền N (các gốc 247–258). Vùng hoạt động được che phủ bởi một miền nắp (lid domain) tạo bởi cầu disulfua giữa Cys237 và Cys261. Ở trạng thái mở, nắp dịch chuyển làm lộ hốc xúc tác để liên kết cơ chất lipid.
  • Mô hình tính toán Docking phân tử: Dự đoán cấu hình liên kết tối ưu giữa phối tử (ligand) và thụ thể (receptor) dựa trên việc cực tiểu hóa năng lượng tự do Gibbs ($\Delta G$) và tính toán hàm tính điểm (scoring function) tích hợp các thành phần năng lượng Van der Waals, liên kết hydro, tương tác tĩnh điện, desolvat hóa và năng lượng xoắn nội tại.
  • Quy tắc 5 Lipinski (Rule of Five - RO5): Bộ tiêu chí đánh giá tính giống thuốc dùng đường uống của hợp chất dựa trên 5 thông số hóa lý: 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 ($\log P$) < 5; và độ khúc xạ phân tử (Molar Refractivity - MR) trong khoảng 40–130. Phân tử đạt tối thiểu 2 trong 5 tiêu chí được xem là có tiềm năng sinh khả dụng đường uống.
  • Mô hình ADMET điện toán: Đánh giá các thuộc tính hấp thu (Absorption), phân bố (Distribution), chuyển hóa (Metabolism), thải trừ (Excretion) và độc tính (Toxicity) thông qua thuật toán đồ thị và học máy trên nền tảng pkCSM.

Quy trình thực nghiệm tính toán bao gồm các bước:

  • Nguồn dữ liệu cấu trúc: Cấu trúc 3D tinh thể HPL ở độ phân giải 2,46 Å được truy xuất từ Ngân hàng Dữ liệu Protein (RCSB Protein Data Bank) với mã định danh PDB: 1LPB. Cấu trúc phức hợp này chứa sẵn phân tử đồng kết tinh là chất ức chế Methoxy undecyl phosphinic acid (MUP). Cấu trúc 2D/3D của 50 dẫn xuất coumarin và chất đối chứng Orlistat được thu thập từ cơ sở dữ liệu PubChem.
  • Chuẩn bị dữ liệu: Cấu trúc protein 1LPB được xử lý trên Discovery Studio 2021 Client để loại bỏ các phân tử nước, ion, chất nền và ligand MUP; sau đó thêm hydro phân cực, gán trường lực Kollman và chuyển sang định dạng .pdbqt bằng AutoDock Tools 1.5.7. Các phân tử coumarin được tối ưu hóa năng lượng hình học bằng phần mềm Avogadro thông qua thuật toán Gradient liên hợp (Conjugate Gradients) rồi chuyển đổi sang file .pdbqt.
  • Thẩm định quy trình (Re-docking): Ligand nội tại MUP sau khi tách khỏi hốc liên kết được dock ngược trở lại thụ thể bằng AutoDock Vina. Tính hợp lý của quy trình được xác thực khi độ lệch bình phương trung bình gốc (RMSD) giữa cấu hình mô phỏng và cấu hình tinh thể thực nghiệm đạt giá trị nhỏ hơn 1,5 Å (kết quả thực tế đo được qua Chimera là 1,271 Å).
  • Mô phỏng Docking phân tử: Sử dụng phần mềm AutoDock Vina với vùng tìm kiếm (grid box) có kích thước $70 \times 70 \times 70$ Å, khoảng cách mắt lưới (grid spacing) 0,375 Å, tọa độ tâm đặt tại $x = 8,431$; $y = 24,417$; $z = 52,623$ bao trọn bộ ba xúc tác và các acid amin lân cận.
  • Sàng lọc tính giống thuốc và ADMET: 10 chất có năng lượng liên kết thấp nhất và tương tác tốt với acid amin khóa được kiểm tra qua công cụ trực tuyến Quy tắc Lipinski của IIT Delhi. 8 chất đạt chuẩn tiếp tục được nạp mã SMILES vào nền tảng pkCSM để tính toán hơn 30 thông số dược động học và an toàn sinh học.

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

Chương 1: Tổng quan

Chương 1 tổng hợp các cơ sở lý luận về bệnh lý béo phì, cấu trúc sinh học enzym HPL, hóa học nhóm coumarin và phương pháp luận in silico:

  • Trình bày định nghĩa béo phì theo WHO, bảng phân loại cân nặng theo chỉ số khối cơ thể (BMI từ thiếu cân $< 18,5$ đến béo phì độ III $\ge 40$ $\text{kg/m}^2$), phân tích nguyên nhân đa yếu tố (chế độ dinh dưỡng dư thừa calo, lối sống tĩnh tại, môi trường kinh tế - xã hội, các gen di truyền trội gây béo phì, tổn thương vùng dưới đồi, hội chứng Cushing, Prader-Willi, Alstrom, Bardet-Biedl, Carpenter, Cohen và tăng cân do thuốc chống loạn thần, corticosteroid, chẹn beta).
  • Hệ thống hóa các phương pháp điều trị béo phì: can thiệp lối sống (giảm 500 kcal/ngày, tập thể dục hiếu khí $\ge 150$ phút/tuần), liệu pháp ngoại khoa, và danh mục thuốc được FDA phê duyệt gồm Orlistat, Liraglutide, Phentermine-Topiramate, Naltrexone-Bupropion, Gelesis (Plenity), Setmelanotide (chất chủ vận MC4R) và Phentermine.
  • Phân tích chi tiết chu trình tiêu hóa mỡ và vai trò phân cắt 50–70% lipid của lipase tụy; làm rõ cấu trúc tiểu phần HPL với bộ ba xúc tác Ser152 - Asp176 - His263 cùng miền nắp chuyển đổi giữa hai trạng thái đóng và mở.
  • Tổng quan về hóa thực vật nhóm coumarin (2H-1-benzopyran-2-one), nêu rõ lịch sử phân lập từ hạt đậu Tonka năm 1820, phân bố tại các họ thực vật chính (Rutaceae, Umbelliferae, Fabaceae...), và dẫn chứng một số coumarin có hoạt tính kháng béo phì đã biết như esculetin từ vỏ cây Fraxinus rhynchophylla (ức chế biệt hóa mỡ qua con đường phụ thuộc PPAR$\gamma$) hoặc các hợp chất chiết từ Peucedanum japonicum Thunb và nấm Grifola frondosa.
  • Trình bày nguyên lý thiết kế thuốc nhờ máy tính (CADD/SBDD), thuật toán tìm kiếm cấu hình và các hàm tính điểm nhiệt động học trong docking, quy tắc 5 Lipinski và các công cụ dự đoán ADMET trực tuyến.

Chương 2: Nguyên liệu, thiết bị, phương pháp nghiên cứu

Chương 2 mô tả chi tiết hạ tầng kỹ thuật, đối tượng dữ liệu và các bước tiến hành thực nghiệm tính toán:

  • Phần cứng: Máy tính xách tay HP Elitebook Folio 9470m chạy hệ điều hành Windows 10.
  • Phần mềm: MGL Tools, AutoDock Vina 1.1.2, UCSF Chimera, Avogadro, Discovery Studio 2021 Client, công cụ trực tuyến Lipinski 5 của Supercomputing Facility for Bioinformatics & Computational Biology (IIT Delhi) và máy chủ web pkCSM.
  • Trình bày quy trình 5 bước thẩm định re-docking ligand MUP để xác định sai số RMSD; quy trình khử nước, khử ion và tối ưu trường lực Kollman cho protein 1LPB; kỹ thuật tối ưu hóa hình học cho 50 cấu trúc coumarin bằng thuật toán Conjugate Gradients trong Avogadro.
  • Thiết lập hộp tọa độ lưới grid box ($70 \times 70 \times 70$ Å, tâm $x = 8,431$; $y = 24,417$; $z = 52,623$) tại trung tâm hoạt động HPL; xác lập 3 tiêu chí chọn lọc sau docking: điểm năng lượng tự do thấp hơn Orlistat, RMSD thấp, và tạo liên kết trực tiếp với các acid amin thiết yếu Ser152, Asp176, His263 và Phe77.
  • Xác định phương thức đánh giá đặc tính giống thuốc theo quy tắc 5 Lipinski và đo lường các nhóm chỉ số ADMET (tính thấm Caco-2, hấp thu ruột, thể tích phân bố VDss, hàng rào máu não BBB, hệ thần kinh trung ương CNS, chuyển hóa qua enzym CYP, thanh thải thận và độc tính tế bào/cơ quan).

Chương 3: Kết quả

Chương 3 trình bày các phát hiện và hệ thống dữ liệu định lượng thu được qua các giai đoạn sàng lọc:

  • Thẩm định docking: Quá trình re-dock MUP vào thụ thể 1LPB đạt giá trị RMSD là 1,271 Å ($< 1,5$ Å), xác nhận mô hình mô phỏng đạt độ tin cậy. MUP tạo liên kết hydro với Ser152, Phe77, His263 và tương tác alkyl/$\pi$-alkyl với Tyr114, Pro180, Phe215, Ala178, Ile209, Val210.
  • Sàng lọc 50 hợp chất coumarin: Toàn bộ 50 hợp chất đều có năng lượng liên kết âm hơn thuốc đối chứng Orlistat ($-6,7$ kcal/mol), dao động từ $-7,0$ đến $-10,9$ kcal/mol. Phân bố theo 4 nhóm cấu trúc:
    • Coumarin đơn (14 chất): từ $-7,0$ kcal/mol (Coumarin, Herniarin) đến $-8,6$ kcal/mol (Cichoriin, Novobioxin).
    • Furanocoumarin (17 chất): từ $-7,5$ kcal/mol (Sphondin, Isobergapten) đến $-9,6$ kcal/mol (Ammijin).
    • Pyranocoumarin (12 chất): từ $-8,0$ kcal/mol (Seselin, Alloxanthoxyletin) đến $-9,2$ kcal/mol (Xanthalin).
    • Nhóm coumarin khác (7 chất): từ $-8,7$ kcal/mol (Scandenin) đến $-10,9$ kcal/mol (Frutinon A).
  • Lựa chọn 10 hợp chất có năng lượng liên kết tốt nhất và tương tác chọn lọc với Ser152, Phe77, His263 và Phe215.
STT Tên hợp chất Năng lượng liên kết (kcal/mol) Liên kết Hydro Tương tác $\pi$-$\pi$
1 Frutinon A -10,9 SER152, PHE77, HIS263 TYR114, PHE215
2 Chartreusin -10,5 SER152, PHE215 PHE77
3 3''-Demetylchartreusin -10,2 SER152 TYR114, PHE77, PHE215
4 Ammijin -9,6 SER152, TYR114 PHE77
5 Hydrangenol -9,6 SER152, HIS263 TYR114, PHE215
6 Peucenidin -9,5 SER152, PHE77, HIS263, ARG256 TYR114, PHE215
7 Calophyllolid -9,5 SER152 PHE77, TYR114, PHE215
8 Phyllocoumarin -9,5 HIS263 PHE77, TYR114, PHE215
9 Edultin -9,4 SER152, HIS263, PHE77 TYR114, PHE215
10 Xanthalin -9,2 SER152, PHE215 TYR114
  • Đánh giá quy tắc 5 Lipinski cho 10 hợp chất: Chartreusin (MW = 640,594 Da; HBD = 5; HBA = 14; MR = 154,03) và 3''-Demetylchartreusin (MW = 626,567 Da; HBD = 6; HBA = 14; MR = 149,239) vi phạm nhiều tiêu chí nên bị loại bỏ. Tám hợp chất còn lại (Frutinon A, Ammijin, Hydrangenol, Peucenidin, Calophyllolid, Phyllocoumarin, Edultin, Xanthalin) đáp ứng đầy đủ hoặc vi phạm không quá 1 tiêu chí, đủ điều kiện làm thuốc đường uống.
Tên hợp chất Phân tử khối (Da) Nhóm cho H (HBD) Nhóm nhận H (HBA) $\log P$ Khúc xạ phân tử (MR)
Frutinon A 264,236 0 4 2,592 70,397
Chartreusin 640,594 5 14 0,948 154,030
3''-Demetylchartreusin 626,567 6 14 0,294 149,239
Ammijin 408,403 4 9 -0,483 98,162
Hydrangenol 256,257 2 4 2,552 68,074
Peucenidin 386,400 0 7 3,272 99,636
Calophyllolid 416,473 0 5 5,378 119,652
Phyllocoumarin 342,303 4 7 1,773 85,627
Edultin 386,400 0 7 3,272 99,636
Xanthalin 426,465 0 7 4,218 113,393
  • Đánh giá ADMET cho 8 hợp chất tiềm năng:
    • Hấp thu: 6 hợp chất có tính thấm màng Caco-2 cao ($> 0,9 \times 10^{-6}\text{ cm/s}$) và tỉ lệ hấp thu qua ruột ở người đạt mức cao: Frutinon A ($100%$), Peucenidin ($100%$), Edultin ($100%$), Xanthalin ($100%$), Calophyllolid ($98,5%$), Hydrangenol ($93,86%$). Ammijin ($45,0%$) và Phyllocoumarin ($79,1%$) có mức hấp thu thấp hơn.
    • Phân bố: Calophyllolid có thể tích phân bố cao nhất ($\log VDss = 0,47$). Khả năng qua hàng rào máu não ($\log BB > 0,3$) ghi nhận ở Frutinon A ($0,387$) và Xanthalin ($0,866$). Tuy nhiên tính thấm hệ thần kinh trung ương ($\log PS$) của Xanthalin là $-2,894$ ($< -2$) nên không thấm vào mô não, trong khi Frutinon A và Calophyllolid có $\log PS > -2$.
    • Chuyển hóa và Thải trừ: 5 hợp chất (Frutinon A, Peucenidin, Calophyllolid, Edultin, Xanthalin) là cơ chất của CYP3A4; không chất nào là cơ chất của CYP2D6 hay OCT2 ở thận. Bốn chất (Peucenidin, Calophyllolid, Edultin, Xanthalin) có tính chất ức chế CYP3A4. Tốc độ thanh thải toàn phần cao nhất ở Xanthalin ($1,047 \text{ log ml/min/kg}$) và Edultin ($0,971 \text{ log ml/min/kg}$).
    • Độc tính: Frutinon A và Hydrangenol dương tính với thử nghiệm gây đột biến AMES. Calophyllolid và Phyllocoumarin có nguy cơ ức chế kênh kali hERG II gây độc tính trên tim. Bốn chất Peucenidin, Edultin, Xanthalin và Ammijin có cảnh báo độc tính trên gan. Toàn bộ 8 hợp chất không gây kích ứng da và không ức chế hERG I.

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 thu được và đối chiếu với cơ chế sinh học:

  • Bàn luận về mối quan hệ giữa ái lực liên kết và vị trí gắn kết: Việc cả 10 hợp chất dẫn đầu đều tạo liên kết hydro trực tiếp với Ser152 và tương tác với các gốc His263, Phe77 chứng minh khả năng phong tỏa bộ ba xúc tác trung tâm của HPL, tương tự cơ chế ức chế của phối tử nội tại MUP nhưng vượt trội hơn Orlistat về mặt năng lượng tính toán.
  • Đánh giá toàn diện các thông số ADMET: Mặc dù Frutinon A có ái lực mạnh nhất ($-10,9$ kcal/mol), hợp chất này lại dương tính với thử nghiệm đột biến gen AMES và có khả năng thấm vào thần kinh trung ương. Calophyllolid có năng lượng tốt ($-9,5$ kcal/mol) nhưng lại mang nguy cơ ức chế kênh hERG II gây độc tính tim mạch và $\log P > 5$.
  • Tổng hợp sàng lọc đa tiêu chí xác định 3 ứng viên có tiềm năng phát triển thuốc cao nhất gồm Peucenidin, Edultin và Xanthalin. Cả 3 chất đều có ái lực liên kết mạnh ($-9,2$ đến $-9,5$ kcal/mol), âm tính với thử nghiệm AMES, không ức chế hERG I/II, độ hấp thu ruột đạt $100%$ và tính thấm Caco-2 cao.
  • Nêu rõ các cảnh báo cần khắc phục: Cả 3 hợp chất được chọn đều được dự đoán có khả năng gây độc tính trên gan và đóng vai trò vừa là cơ chất vừa là chất ức chế enzym CYP3A4, dẫn đến nguy cơ tương tác thuốc bất lợi nếu phối hợp điều trị với các dược chất chuyển hóa qua cùng con đường.

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

Nghiên cứu mang lại các kết quả định lượng cụ thể:

  • Xác lập quy trình docking đạt chuẩn trên thụ thể HPL (PDB: 1LPB) với chỉ số RMSD kiểm định là 1,271 Å.
  • Xây dựng bộ dữ liệu năng lượng liên kết của 50 dẫn chất coumarin tự nhiên đối với enzym HPL, chứng minh $100%$ hợp chất thử nghiệm có năng lượng liên kết âm hơn Orlistat ($-6,7$ kcal/mol).
  • Phân loại và chỉ ra 8 hợp chất thỏa mãn các tiêu chí hóa lý đường uống của quy tắc Lipinski (loại trừ Chartreusin và 3''-Demetylchartreusin do phân tử khối $> 600$ Da và số nhóm nhận/cho hydro quá lớn).
  • Phân tích chi tiết hồ sơ an toàn và dược động học của 8 hợp chất, từ đó lựa chọn được 3 ứng viên coumarin tiềm năng nhất gồm: Peucenidin ($-9,5$ kcal/mol), Edultin ($-9,4$ kcal/mol) và Xanthalin ($-9,2$ kcal/mol). Cả 3 hợp chất đều gắn chọn lọc vào Ser152, hấp thu qua đường tiêu hóa $100%$, thấm màng Caco-2 tốt và không có độc tính gây đột biến hay độc tim mạch.
  • Về mặt đóng góp: Đề tài đã bước đầu chứng minh tiềm năng ức chế lipase tuyến tụy của khung cấu trúc coumarin tự nhiên thông qua sàng lọc điện toán, cung cấp cơ sở dữ liệu cấu trúc - tương tác chi tiết giúp định hướng cho các nghiên cứu tổng hợp, bán tổng hợp và tối ưu hóa phân tử nhằm tạo ra thuốc điều trị béo phì mới có hiệu quả ức chế HPL chọn lọc.

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

Khóa luận tự chỉ ra các hạn chế cố hữu trong mô hình nghiên cứu:

  • Toàn bộ kết quả đều dựa trên các mô hình mô phỏng tính toán và thuật toán dự đoán (in silico), do đó chưa phản ánh hoàn toàn động học phân tử linh động trong môi trường sinh học thực tế và chưa được thẩm định bằng thực nghiệm hóa sinh.
  • Các ứng viên tiềm năng được lựa chọn (Peucenidin, Edultin, Xanthalin) vẫn còn tồn tại hạn chế về mặt dược lý học: có nguy cơ gây độc tính trên gan và khả năng ức chế CYP3A4 gây tương tác thuốc.

Hướng nghiên cứu tiếp theo được đề xuất:

  • Tiến hành các thử nghiệm ức chế enzym lipase tuyến tụy in vitro để xác định giá trị nồng độ ức chế bán tối đa ($\text{IC}_{50}$) thực tế của Peucenidin, Edultin và Xanthalin.
  • Thực hiện các nghiên cứu in vivo trên mô hình động vật béo phì do chế độ ăn giàu chất béo nhằm đánh giá tác dụng hạ lipid và độc tính thực tế.
  • Tối ưu hóa cấu trúc hóa học của các phân tử dẫn đường để làm giảm độc tính trên gan và giảm ái lực ức chế enzym chuyển hóa CYP3A4.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên 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 và Tin sinh học:

  • Cung cấp quy trình chi tiết từng bước về thực hành sàng lọc ảo (virtual screening) kết hợp giữa AutoDock Vina, UCSF Chimera, Avogadro và pkCSM trên protein đích dạng tinh thể PDB.
  • Dữ liệu tương tác giữa các acid amin thuộc bộ ba xúc tác Ser152 - Asp176 - His263 với các nhóm coumarin là tài liệu tham khảo giá trị cho các nghiên cứu thiết kế phân tử ức chế HPL.
  • Bảng tổng hợp các thông số hóa lý, năng lượng liên kết và đặc tính ADMET của 50 dẫn chất coumarin là nguồn dữ liệu thô phục vụ cho các nghiên cứu tối ưu hóa cấu trúc dẫn chất tự nhiên hướng điều trị rối loạn chuyển hóa lipid.

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

1. Tiêu chí kiểm định quy trình docking (re-docking) của khóa luận là gì và kết quả đạt được như thế nào?
Khóa luận sử dụng ligand đồng kết tinh MUP (Methoxy undecyl phosphinic acid) có sẵn trong cấu trúc tinh thể 1LPB để tách ra và tiến hành dock ngược lại vào protein. Tiêu chuẩn đặt ra là giá trị RMSD giữa phối tử sau khi re-dock so với vị trí tinh thể ban đầu phải nhỏ hơn 1,5 Å. Kết quả đo được trên phần mềm Chimera đạt RMSD = 1,271 Å, khẳng định quy trình docking và bộ thông số hộp lưới là đáng tin cậy.

2. Điểm docking của 10 hợp chất hàng đầu so với chất đối chứng Orlistat chênh lệch ra sao?
Chất đối chứng dương Orlistat có năng lượng liên kết là $-6,7$ kcal/mol. Cả 10 hợp chất hàng đầu đều có năng lượng liên kết âm hơn đáng kể, dao động từ $-9,2$ đến $-10,9$ kcal/mol. Hợp chất có năng lượng thấp nhất là Frutinon A ($-10,9$ kcal/mol), tiếp theo là Chartreusin ($-10,5$ kcal/mol), 3''-Demetylchartreusin ($-10,2$ kcal/mol), Ammijin và Hydrangenol ($-9,6$ kcal/mol), Peucenidin, Calophyllolid và Phyllocoumarin ($-9,5$ kcal/mol), Edultin ($-9,4$ kcal/mol) và Xanthalin ($-9,2$ kcal/mol).

3. Vì sao hai hợp chất Chartreusin và 3''-Demetylchartreusin có năng lượng liên kết rất tốt nhưng bị loại khỏi danh sách phát triển thuốc?
Dù xếp thứ 2 và thứ 3 về ái lực liên kết ($-10,5$ và $-10,2$ kcal/mol), cả hai hợp chất này đều vi phạm nghiêm trọng quy tắc 5 Lipinski dành cho thuốc đường uống: phân tử khối lớn vượt ngưỡng ($640,594$ Da và $626,567$ Da so với chuẩn $< 500$ Da), số nhóm nhận liên kết hydro quá cao (14 nhóm so với chuẩn $< 10$), số nhóm cho liên kết hydro cao (5 và 6 nhóm so với chuẩn $< 5$) và độ khúc xạ phân tử vượt giới hạn ($154,03$ và $149,24$ so với ngưỡng 40–130).

4. Ba hợp chất tiềm năng nhất được lựa chọn cuối cùng là những chất nào và hồ sơ an toàn của chúng có điểm gì cần lưu ý?
Ba hợp chất được đánh giá tốt nhất là Peucenidin, Edultin và Xanthalin. Cả ba đều có ái lực liên kết mạnh với HPL ($-9,2$ đến $-9,5$ kcal/mol), đạt chuẩn Lipinski, độ hấp thu ruột $100%$, không gây đột biến gen (AMES âm tính) và không ức chế kênh hERG. Tuy nhiên, điểm cần lưu ý là cả 3 chất đều được dự đoán có khả năng gây độc tính trên gan và đều có hoạt tính ức chế enzym CYP3A4, đòi hỏi phải được tối ưu hóa cấu trúc trong các bước nghiên cứu tiếp theo.

Kết luận

Khóa luận tốt nghiệp của sinh viên Nguyễn Thị Hồng Hạnh đã thực hiện thành công quy trình sàng lọc in silico 50 dẫn chất coumarin tự nhiên tác động ức chế enzym lipase tuyến tụy người (1LPB). Thông qua các bước đánh giá tuần tự từ docking phân tử, sàng lọc tính chất giống thuốc theo Lipinski đến phân tích toàn diện hơn 30 thông số ADMET, đề tài đã xác định được 3 ứng viên tiềm năng nhất gồm Peucenidin, Edultin và Xanthalin với ái lực liên kết vượt trội so với Orlistat và khả năng hấp thu đường uống tối ưu. Kết quả nghiên cứu cung cấp cơ sở dữ liệu khoa học tin cậy và mở ra hướng đi triển vọng trong việc phát triển các thuốc điều trị béo phì mới có nguồn gốc từ tự nhiên.