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

Theo báo cáo của Cơ quan Nghiên cứu Ung thư Quốc tế (Globocan 2020), ung thư vú là bệnh lý ác tính phổ biến nhất trên phạm vi toàn cầu với khoảng 2,3 triệu ca mắc mới (chiếm 11,7% tổng số ca ung thư) và 685.000 trường hợp tử vong (chiếm 6,9%). Tại Việt Nam, thống kê ghi nhận 21.555 ca mắc mới, đứng đầu trong các loại ung thư ở nữ giới (chiếm 25,8%), và 9.345 ca tử vong. Đa số bệnh nhân ung thư vú tại Việt Nam được chẩn đoán ở giai đoạn II (61,2%), trong khi tỷ lệ phát hiện ở giai đoạn muộn (giai đoạn III và IV) lần lượt là 19,4% và 8,2%.

Về mặt cơ chế phân tử, thụ thể estrogen alpha (ER-alpha hay ERα) là một thành viên thuộc siêu họ thụ thể nhân/steroid, đóng vai trò vừa là đầu dò tín hiệu vừa là yếu tố phiên mã điều hòa biểu hiện gen đích. Thụ thể này chỉ hiện diện ở khoảng dưới 10% tế bào biểu mô tuyến vú bình thường, nhưng lại biểu hiện trong 50% đến 80% các trường hợp khối u vú. Khi phối tử nội sinh 17-β-estradiol (E2) gắn vào vùng liên kết phối tử (LBD), phức hợp sẽ thúc đẩy quá trình phiên mã và tăng sinh tế bào ác tính.

Trong thực hành lâm sàng, liệu pháp nội tiết nhắm vào ERα là lựa chọn điều trị hàng đầu cho bệnh nhân ung thư vú có thụ thể nội tiết dương tính. Các nhóm thuốc chính đang được sử dụng bao gồm thuốc ức chế aromatase (letrozole, anastrozole), thuốc điều hòa chọn lọc thụ thể estrogen (tamoxifen, raloxifen, toremifen) và thuốc làm suy giảm thụ thể estrogen chọn lọc (fulvestrant). Mặc dù mang lại hiệu quả rõ rệt, việc sử dụng các thuốc này trong thời gian dài bộc lộ nhiều tác dụng phụ nghiêm trọng: tamoxifen làm gia tăng nguy cơ ung thư nội mạc tử cung và xuất hiện hiện tượng kháng thuốc sau nhiều năm điều trị; raloxifen làm tăng nguy cơ hình thành cục máu đông, huyết khối tĩnh mạch sâu và thuyên tắc phổi.

Trước thực trạng đó, các hợp chất alkaloid có nguồn gốc tự nhiên với hơn 18.000 cấu trúc đã được phát hiện được xem là nguồn tài nguyên tiềm năng. Một số nghiên cứu in vivo trước đây cho thấy vinca alkaloid có khả năng làm giảm mức protein ERα sau 12 giờ xử lý. Tuy nhiên, các dữ liệu thực nghiệm chưa giải thích đầy đủ các tương tác ở cấp độ phân tử giữa khung cấu trúc alkaloid và thụ thể ERα. Khóa luận tốt nghiệp ngành Dược học của tác giả Bùi Phúc Lộc (2022) tại Trường Đại học Y Dược – Đại học Quốc gia Hà Nội được triển khai nhằm thực hiện hai nhiệm vụ nghiên cứu cụ thể:

  1. Sàng lọc các hợp chất alkaloid có tác dụng ức chế thụ thể ER-alpha bằng phương pháp docking phân tử.
  2. Phân tích 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à cấu trúc không gian ba chiều của thụ thể ER-alpha người và tập hợp 52 hợp chất alkaloid có hoạt tính kháng u đã được ghi nhận. Phạm vi nghiên cứu tập trung vào mô hình tính toán sinh học trên máy tính (in silico), bao gồm mô phỏng docking phân tử, đánh giá quy tắc Lipinski và dự đoán các thông số ADMET.

Cơ sở lý thuyết và phương pháp nghiên cứu

Nghiên cứu được xây dựng dựa trên phương pháp phát minh thuốc có sự hỗ trợ của máy tính (CADD), cụ thể là nhánh phát minh thuốc dựa trên cấu trúc (SBDD). Phương pháp này sử dụng cấu trúc ba chiều của đại phân tử sinh học đích để mô phỏng tương tác liên kết với các phối tử phân tử nhỏ, từ đó tính toán ái lực liên kết thông qua năng lượng tự do Gibbs ($\Delta G$).

Các khung lý thuyết và bộ quy tắc được áp dụng trong khóa luận gồm:

  • Nguyên lý docking phân tử: Sử dụng hàm tính điểm để ước tính ái lực liên kết dựa trên năng lượng nội tại của phối tử, năng lượng tự do xoắn, tương tác Van der Waals, liên kết hydro, năng lượng desolvat và tương tác tĩnh điện.
  • Quy tắc Lipinski (Rule of 5): Bộ tiêu chuẩn đánh giá tính chất giống thuốc đối với hợp chất dùng đường uống, quy định hợp chất không được vi phạm quá 1 trong 4 chỉ số: 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 ($\text{LogP}$) < 5.
  • Mô hình dự đoán ADMET: Nền tảng pkCSM sử dụng đồ thị đặc tính phân tử để dự đoán các thông số hấp thu ở ruột người (HIA), tính thấm qua màng Caco2, thể tích phân bố ($V_{\text{dss}}$), tính thấm qua hàng rào máu não (BBB), tính thấm vào hệ thần kinh trung ương (CNS), khả năng tương tác với hệ enzym Cytochrome P450 (CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4), độ thanh thải toàn phần và các chỉ số độc tính (AMES, ức chế kênh kali hERG I/II, độc tính gan, liều chết trung bình $\text{LD}_{50}$ đường uống ở chuột).

Nguồn dữ liệu và công cụ phần mềm phục vụ nghiên cứu bao gồm:

  • Cấu trúc thụ thể đích: Cấu trúc tinh thể tia X của thụ thể ER-alpha người được tải về từ Ngân hàng Dữ liệu Protein (Protein Data Bank - PDB) với mã định danh PDB ID là 3ERT (độ phân giải 1,9 Å), chứa phối tử đồng kết tinh là chất đối kháng chọn lọc 4-hydroxytamoxifen (OHT).
  • Cấu trúc phối tử: Dữ liệu cấu trúc hóa học của 52 hợp chất alkaloid thu thập từ cơ sở dữ liệu PubChem dưới định dạng chuẩn.
  • Hệ thống phần mềm và thiết bị: Sử dụng máy tính cá nhân trên nền tảng hệ điều hành Windows 10; phần mềm Autodock Vina 1.1.2 và MGL Tools 1.5.6 phục vụ tính toán docking; phần mềm Avogadro tối ưu hóa năng lượng phối tử bằng phương pháp Gradient liên hợp (Conjugate Gradients); phần mềm UCSF Chimera 1.14 tính toán sai lệch bình phương trung bình gốc (RMSD); phần mềm Discovery Studio 2021 Client xử lý đại phân tử và phân tích tương tác 2D/3D; công cụ trực tuyến pkCSM phục vụ đánh giá tính chất giống thuốc và đặc tính ADMET.

Quy trình nghiên cứu thực tế diễn ra qua 3 bước liên hoàn:

  • Bước 1: Sàng lọc bằng docking phân tử. Thẩm định quy trình docking bằng phương pháp re-dock ligand đồng kết tinh 4-hydroxytamoxifen vào túi gắn của thụ thể 3ERT (yêu cầu RMSD < 1,5 Å). Sau đó, tiến hành dock đồng loạt 52 phân tử alkaloid vào hộp lưới kích thước $40 \times 40 \times 40 \text{ \AA}$ (khoảng cách lưới 0,375 Å, tọa độ tâm $x = 30,282$; $y = -1,913$; $z = 24,207$). Lọc ra các chất có năng lượng liên kết $\Delta G \le -9,8 \text{ kcal/mol}$ (ngưỡng của 4-hydroxytamoxifen).
  • Bước 2: Sàng lọc đặc tính giống thuốc. Đưa mã SMILES của các phối tử vượt qua vòng docking vào công cụ pkCSM để kiểm tra sự tuân thủ quy tắc 5 Lipinski.
  • Bước 3: Dự đoán dược động học và an toàn ADMET. Tính toán các chỉ số hấp thu, phân bố, chuyển hóa, thải trừ và độc tính nhằm lựa chọn ứng viên phù hợp nhất cho định hướng phát triển thuốc.

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 các tài liệu khoa học liên quan đến bốn nhóm vấn đề chính:

  • Bệnh ung thư vú: Trình bày định nghĩa, cơ chế xâm lấn từ biểu mô ống (85%) và tiểu thùy (15%), các yếu tố nguy cơ không thể thay đổi (tuổi tác, tiền sử gia đình, đột biến gen) và có thể thay đổi (lối sống, béo phì, rượu bia), cùng các phương pháp can thiệp ngoại khoa, xạ trị (chiếu ngoài, áp sát, xạ trị trong mổ IORT) và điều trị toàn thân.
  • Thụ thể ER-alpha: Phân tích cấu trúc phân tử của ERα gồm 3 miền chức năng chính (NTD, DBD, LBD), các miền kích hoạt phiên mã AF1 và AF2. Cấu trúc miền gắn phối tử (LBD) dạng cầu chứa 11 vòng xoắn alpha (H1–H12), trong đó khoang gắn kết kỵ nước được bao bọc bởi 22 gốc acid amin quan trọng (như Glu353, Arg394, His524).
  • Thuốc điều trị ung thư vú nội tiết: Tổng thuật cơ chế tác dụng và các tác dụng không mong muốn của tamoxifen, raloxifen, toremifen, fulvestrant, letrozole và anastrozole.
  • Hợp chất alkaloid và nghiên cứu in silico: Khái quát tính đa dạng hóa học của alkaloid tự nhiên, vai trò của nguyên tử nitơ và các nhóm chức cho/nhận proton trong việc tạo liên kết sinh học. Trình bày tổng quan các thuật toán docking phân tử, hàm tính điểm, quy tắc Lipinski và vai trò của các công cụ ADMET trực tuyến trong việc giảm tỷ lệ thất bại ở các giai đoạn thử nghiệm lâm sàng.

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

Chương này mô tả chi tiết thông số kỹ thuật của các đối tượng nghiên cứu và thiết lập quy trình tính toán:

  • Chuẩn bị dữ liệu: Quy trình xử lý cấu trúc thụ thể 3ERT (loại bỏ phân tử nước, tách phối tử gốc, bổ sung nguyên tử hydro phân cực, gán điện tích Kollman và xuất định dạng pdbqt bằng Autodock Tools). Tối ưu hóa cấu trúc 52 phối tử alkaloid bằng Avogadro và chuyển đổi cấu trúc phân tử sang định dạng SMILES.
  • Quy trình thẩm định và docking: Thiết lập thông số hộp tìm kiếm (grid box) tại trung tâm hoạt động của ERα; tiêu chuẩn thẩm định phương pháp re-dock với độ sai lệch RMSD dưới 1,5 Å.
  • Bộ tiêu chí đánh giá: Xác định các mốc tham chiếu cho quy tắc Lipinski và các ngưỡng đánh giá ADMET trên pkCSM (hấp thu tốt khi HIA > 70%, tính thấm màng Caco2 cao khi log Papp > 0,9, thấm hàng rào máu não khi logBB > 0,3, thấm thần kinh trung ương khi logPS > -2).

Chương 3: Kết quả

Chương kết quả thể hiện các dữ liệu thực nghiệm tính toán thu được:

  • Thẩm định quy trình: Re-dock 4-hydroxytamoxifen vào thụ thể 3ERT đạt giá trị RMSD = 1,5 Å, chứng minh độ tin cậy của thuật toán. Năng lượng liên kết của 4-hydroxytamoxifen đạt $\Delta G = -9,8 \text{ kcal/mol}$, tạo liên kết hydro với Glu353, Arg394 và tương tác kỵ nước alkyl/$\pi$-alkyl với Ala350, Leu387, Leu391, Leu525, Met421.
  • Sàng lọc docking 52 alkaloid: Xác định được 4 hợp chất có điểm số năng lượng liên kết $\Delta G \le -9,8 \text{ kcal/mol}$ gồm: Penitrem A (-9,8 kcal/mol), Pityriacitrin B (-9,8 kcal/mol), Penitrem B (-10,1 kcal/mol) và Tubocuraran (-10,8 kcal/mol theo Bảng 3.4 hoặc -10,2 kcal/mol tại phần thảo luận).
  • Kiểm tra quy tắc Lipinski: Pityriacitrin B đạt đầy đủ 4/4 tiêu chí; Tubocuraran đạt 3/4 tiêu chí (chấp nhận được theo quy ước không vi phạm quá 1 tiêu chí). Penitrem A và Penitrem B vi phạm 2 tiêu chí (trọng lượng phân tử MW và hệ số LogP) nên bị loại khỏi các bước phân tích tiếp theo.
  • Phân tích ADMET và tương tác phối tử: Pityriacitrin B có tỷ lệ hấp thu ruột đạt 72,002% (xấp xỉ 73%), không qua hàng rào máu não (logBB = -0,029), không ức chế và không phải cơ chất của CYP3A4, CYP2D6, không ức chế hERG I/II, nhưng có dự đoán độc tính gan. Tubocuraran có tỷ lệ hấp thu ruột 95%, là cơ chất CYP3A4, ức chế CYP2D6 và ức chế kênh hERG II (nguy cơ độc tính tim). Về tương tác, Pityriacitrin B tạo liên kết hydro với His356, Val355 và liên kết với Ile358, Trp360.

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

Chương bàn luận phân tích chuyên sâu về ý nghĩa sinh học và hóa dược của các kết quả:

  • Phân tích tiềm năng của Pityriacitrin B: Đây là một indole alkaloid có nguồn gốc từ nấm men Malassezia furfur, trước đây mới chỉ ghi nhận hoạt tính chống tăng sinh yếu trên ung thư vú và tuyến tiền liệt. Nghiên cứu này làm sáng tỏ khả năng tương tác ổn định của Pityriacitrin B tại túi gắn LBD của ERα với năng lượng tự do ngang bằng phối tử đối kháng chuẩn.
  • So sánh với Tubocuraran: Mặc dù Tubocuraran có điểm số docking cao hơn (-10,8 kcal/mol so với -9,8 kcal/mol), việc hợp chất này ức chế hERG II, ức chế enzym chuyển hóa CYP2D6 và là cơ chất của CYP3A4 tạo ra nguy cơ cao về độc tính trên tim cũng như tương tác thuốc phức tạp.
  • Hạn chế của mô hình in silico: Phân tích tính chất tĩnh của cấu trúc tinh thể tia X, sự phụ thuộc vào các thuật toán ước tính năng lượng và sự cần thiết phải kiểm chứng bằng mô phỏng động lực học phân tử (MD) cũng như thực nghiệm sinh học in vitro/in vivo.

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

Nghiên cứu đã định lượng hóa và sàng lọc thành công các hợp chất alkaloid tiềm năng từ tập dữ liệu ban đầu thông qua chuỗi tiêu chí khắt khe về năng lượng gắn kết, tính chất lý hóa và dược động học.

Tên hợp chất PubChem ID Năng lượng liên kết ($\Delta G$, kcal/mol) Trọng lượng phân tử (MW, Da) Số nhóm cho H (HBD) Số nhóm nhận H (HBA) LogP Kết luận Lipinski
4-Hydroxytamoxifen (Đối chứng) 449459 -9,8 387,5 1 2 5,5 Đạt (vi phạm 1)
Penitrem A 6610243 -9,8 634,4 3 6 5,32 Loại (vi phạm 2)
Penitrem B 3084087 -10,1 > 500 - - > 5 Loại (vi phạm 2)
Pityriacitrin B 56649543 -9,8 < 500 < 5 < 10 < 5 Đạt (4/4 tiêu chí)
Tubocuraran 9548771 -10,8 - - - - Đạt (3/4 tiêu chí)

Đối với hai hợp chất vượt qua vòng sàng lọc Lipinski, các thông số dược động học và an toàn ADMET được tổng hợp và đối chiếu trực tiếp:

Thông số ADMET Pityriacitrin B Tubocuraran Ý nghĩa tham chiếu
Độ tan trong nước (log mol/L) -2,641 - Chỉ số hòa tan
Tính thấm Caco2 (log Papp, $10^{-6}$ cm/s) 0,749 0,828 Ngưỡng cao khi > 0,9
Hấp thu ở ruột người (HIA) (%) 72,002% (~73%) 95,0% Hấp thu tốt khi > 70%
Thể tích phân bố ($\log V_{\text{dss}}$) -1,267 - Phân bố mô kém khi < -0,15
Thấm hàng rào máu não (logBB) -0,029 < 0,3 Không qua BBB (ngưỡng > 0,3)
Thấm hệ thần kinh TW (logPS) -1,388 (hoặc -1,882) > -2 Thấm kém qua CNS (ngưỡng > -2)
Chuyển hóa qua CYP3A4 Không (cơ chất/ức chế) Cơ chất Nguy cơ tương tác chuyển hóa
Chuyển hóa qua CYP2D6 Không (cơ chất/ức chế) Ức chế Nguy cơ tương tác thuốc
Ức chế kênh hERG I / II Không / Không Không / Có Độc tính tim mạch (hERG II)
Độc tính gan Có Có Cần tối ưu hóa cấu trúc
Đột biến gen (AMES test) Không Không Nguy cơ gây ung thư/đột biến
Kích ứng da Không Không Phản ứng tại chỗ
Độ thanh thải toàn phần (log ml/min/kg) 0,78 (hoặc 0,439) 0,439 Tốc độ đào thải qua thận

Đóng góp của khóa luận thể hiện ở việc cung cấp dữ liệu mô phỏng cấp độ phân tử đầu tiên về cơ chế gắn kết của Pityriacitrin B lên thụ thể ER-alpha. Tác giả đã chứng minh Pityriacitrin B hình thành các liên kết hydro đặc hiệu với His356 và Val355, đồng thời duy trì ái lực liên kết tương đương với chất đối kháng chuẩn 4-hydroxytamoxifen. Kết quả này mở ra một hướng tiếp cận cụ thể cho việc nghiên cứu cấu trúc indole alkaloid tự nhiên trong ức chế thụ thể estrogen alpha.

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

Khóa luận ghi nhận các hạn chế chính xuất phát từ phương pháp và phạm vi nghiên cứu:

  • Hạn chế của mô hình in silico tĩnh: Nghiên cứu chỉ thực hiện docking phân tử trên một cấu trúc tinh thể cố định của thụ thể ERα (PDB ID: 3ERT). Phương pháp này chưa tính đến tính linh động toàn phần của protein trong môi trường sinh lý, điều mà mô phỏng động lực học phân tử (Molecular Dynamics - MD) có thể giải quyết.
  • Quy mô tập dữ liệu phối tử: Quá trình sàng lọc mới chỉ thực hiện trên tập hợp 52 hợp chất alkaloid tải về từ PubChem, chưa bao phủ toàn diện kho tàng phong phú của nhóm hợp chất alkaloid trong tự nhiên.
  • Hồ sơ độc tính của các ứng viên hàng đầu: Cả hai chất được chọn đều được dự đoán có độc tính trên gan; riêng Tubocuraran còn có nguy cơ gây độc tính trên tim thông qua ức chế kênh hERG II và ảnh hưởng đến hệ enzym chuyển hóa CYP.
  • Thiếu thực nghiệm kiểm chứng: Các kết luận về ái lực liên kết và thông số ADMET hoàn toàn dựa trên thuật toán tính toán và dự đoán sinh học, chưa có dữ liệu thực nghiệm in vitro (như thử nghiệm gắn kết thụ thể, thử nghiệm độc tế bào trên các dòng ung thư vú dương tính với ER) và thử nghiệm in vivo đối chứng.

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

  • Mở rộng quy mô sàng lọc ảo trên các cơ sở dữ liệu hóa học lớn hơn đối với các phân lớp alkaloid.
  • Ứng dụng mô phỏng động lực học phân tử (MD simulation) nhằm đánh giá độ bền nhiệt động và sự ổn định theo thời gian của phức hợp Pityriacitrin B – ERα.
  • Tiến hành thiết kế, bán tổng hợp hoặc biến đổi cấu trúc hóa học của Pityriacitrin B nhằm loại bỏ độc tính gan dự đoán, cải thiện tính thấm màng sinh học và tăng cường ái lực liên kết.
  • Triển khai các thử nghiệm sinh học in vitro để xác định nồng độ ức chế tối đa một nửa ($\text{IC}_{50}$) của Pityriacitrin B trên thụ thể ER-alpha và các dòng tế bào ung thư vú nhạy cảm với nội tiết.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo phù hợp 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, Tin sinh học và Công nghệ sinh học dược phẩm:

  • Tham khảo quy trình sàng lọc ảo chuẩn hóa: Toàn bộ các bước từ chuẩn bị cấu trúc thụ thể/phối tử, thẩm định hệ thống bằng re-docking, thiết lập hộp lưới, đến phân tích tương tác không gian bằng phần mềm mở (AutoDock, AutoDock Vina, Chimera, Avogadro, Discovery Studio) được trình bày logic và có thể áp dụng trực tiếp cho các đề tài phát hiện phối tử trên các đích protein khác.
  • Tham khảo phương pháp luận kết hợp: Cách thức kết hợp có hệ thống giữa docking phân tử, bộ quy tắc Lipinski và công cụ dự đoán ADMET (pkCSM) cung cấp một mô hình thực hành mẫu mực cho các nghiên cứu SBDD giai đoạn sớm nhằm tiết kiệm thời gian và kinh phí thực nghiệm.
  • Tài liệu chuyên đề về thụ thể ER-alpha: Phần tổng quan cung cấp hệ thống kiến thức chi tiết về cấu trúc LBD, cơ chế dẫn truyền tín hiệu nội tiết và các tương tác acid amin then chốt trong túi liên kết của thụ thể estrogen alpha.

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

1. Tiêu chuẩn nào được sử dụng để thẩm định độ tin cậy của quy trình docking phân tử trong nghiên cứu này? Nghiên cứu sử dụng phương pháp re-dock ligand đồng kết tinh 4-hydroxytamoxifen vào cấu trúc thụ thể ER-alpha (PDB ID: 3ERT). Quy trình được xem là đạt độ tin cậy khi độ sai lệch vị trí không gian (RMSD) giữa cấu hình sau khi dock lại và cấu hình tinh thể ban đầu nhỏ hơn 1,5 Å. Kết quả thực nghiệm tính toán của khóa luận đạt RMSD = 1,5 Å và năng lượng liên kết $\Delta G = -9,8 \text{ kcal/mol}$.

2. Tiêu chí lựa chọn các hợp chất alkaloid tiềm năng từ kết quả docking phân tử là gì? Tiêu chí tiên quyết là hợp chất phải có năng lượng liên kết $\Delta G \le -9,8 \text{ kcal/mol}$ (tương đương hoặc thấp hơn năng lượng liên kết của phối tử chuẩn 4-hydroxytamoxifen), đồng thời cấu dạng liên kết phải có chỉ số RMSD thấp nhất và tạo được tương tác với các acid amin quan trọng trong vùng hoạt động của thụ thể.

3. Tại sao hai hợp chất Penitrem A và Penitrem B bị loại dù có năng lượng liên kết rất tốt? Dù Penitrem A có $\Delta G = -9,8 \text{ kcal/mol}$ và Penitrem B có $\Delta G = -10,1 \text{ kcal/mol}$, cả hai chất này đều vi phạm 2 trong số 4 tiêu chuẩn của quy tắc Lipinski (trọng lượng phân tử MW > 500 Da và hệ số phân bố LogP > 5). Theo quy ước thuốc dùng đường uống không được vi phạm quá 1 tiêu chuẩn, hai hợp chất này không đảm bảo đặc tính giống thuốc nên bị loại.

4. Dựa trên cơ sở nào tác giả kết luận Pityriacitrin B có triển vọng phát triển thuốc tốt hơn Tubocuraran? Mặc dù Tubocuraran có điểm số docking tốt hơn (-10,8 kcal/mol so với -9,8 kcal/mol), Pityriacitrin B vượt trội hơn về mặt an toàn và chuyển hóa: Pityriacitrin B đạt đủ 4/4 tiêu chuẩn Lipinski, không tương tác với các enzym gan CYP3A4, CYP2D6, và quan trọng nhất là không ức chế kênh hERG II. Ngược lại, Tubocuraran là cơ chất CYP3A4, ức chế CYP2D6 và ức chế hERG II, dẫn đến nguy cơ tương tác thuốc phức tạp và độc tính trên tim.

5. Các liên kết hóa học cụ thể giữa Pityriacitrin B và thụ thể ER-alpha gồm những tương tác nào? Tại vị trí hoạt động của thụ thể ERα, Pityriacitrin B hình thành các liên kết hydro trực tiếp với gốc acid amin His356 và Val355, đồng thời tạo liên kết tương tác đặc hiệu với Ile358 và Trp360, giúp neo giữ phối tử ổn định trong khoang liên kết.

Kết luận

Khóa luận tốt nghiệp của tác giả Bùi Phúc Lộc đã xây dựng thành công quy trình sàng lọc in silico trên 52 hợp chất alkaloid nhằm tìm kiếm chất ức chế thụ thể ER-alpha hướng điều trị ung thư vú. Trong số các phân tử khảo sát, Pityriacitrin B nổi bật là ứng viên tiềm năng nhất nhờ ái lực liên kết tương đương chất đối kháng chuẩn 4-hydroxytamoxifen (-9,8 kcal/mol), tuân thủ hoàn toàn quy tắc Lipinski và sở hữu các đặc tính dược động học thuận lợi dù cần được tối ưu để giảm độc tính gan. Công trình cung cấp cơ sở dữ liệu phân tử quan trọng, định hình hướng đi tiếp theo cho các nghiên cứu mô phỏng động lực học và thử nghiệm sinh học thực nghiệm trên nhóm hoạt chất alkaloid tự nhiên.