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

Ung thư vú (UTV) là bệnh lý ác tính phát sinh từ sự tăng sinh mất kiểm soát của tế bào biểu mô tuyến vú, chủ yếu bắt nguồn từ tế bào biểu mô ống dẫn (khoảng 85%) hoặc tiểu thùy (khoảng 15%). Theo số liệu dịch tễ học từ báo cáo Globocan 2020, ung thư vú là căn bệnh ung thư 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à gây ra 685.000 ca tử vong (chiếm 6,9%, đứng hàng thứ năm trong các nguyên nhân tử vong do ung thư). Tại Việt Nam, thống kê năm 2020 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à đứng thứ ba tính trên cả hai giới; đồng thời ghi nhận 9.345 ca tử vong. Bệnh nhân ung thư vú tại các quốc gia đang phát triển thường được phát hiện ở giai đoạn muộn; nghiên cứu tại Việt Nam cho thấy tỷ lệ bệnh nhân được chẩn đoán ở giai đoạn II chiếm 61,2%, giai đoạn III chiếm 19,4% và giai đoạn IV chiếm 8,2%.

Về mặt cơ chế phân tử, thụ thể estrogen alpha (ER-alpha hay ERα) là một yếu tố phiên mã thuộc siêu họ thụ thể nhân (nuclear receptors), đóng vai trò then chốt trong quá trình tăng sinh tế bào tuyến vú. Dù chỉ biểu hiện ở dưới 10% tế bào biểu mô tuyến vú bình thường, ERα lại xuất hiện trong khoảng 50% đến 80% các trường hợp khối u vú ác tính. Khi liên kết với phối tử nội sinh 17-β-estradiol (E2), ERα thúc đẩy quá trình phiên mã các gen đích chi phối sự phân chia và tiến triển của tế bào ung thư, đồng thời chịu sự điều hòa của các protein liên kết như MTA1, SRC-1 và PELP1.

Trong thực hành lâm sàng, liệu pháp nội tiết nhắm đích ERα là chỉ định 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 phổ biến hiện nay bao gồm chất điều biến thụ thể estrogen chọn lọc (tamoxifen, raloxifen, toremifen), chất điều hòa thụ thể estrogen chọn lọc làm suy thoái thụ thể (fulvestrant), hoặc thuốc ức chế men aromatase (letrozole, anastrozole). Tuy nhiên, việc điều trị kéo dài thường phát sinh tình trạng kháng thuốc và gây ra các biến cố bất lợi nghiêm trọng: tamoxifen làm tăng nguy cơ ung thư nội mạc tử cung do chuyển đổi cơ chế phản xạ sinh học từ đối kháng sang chủ vận; raloxifen làm tăng tỷ lệ huyết khối tĩnh mạch sâu và thuyên tắc phổi; toremifen làm tăng nguy cơ ung thư tử cung.

Nhằm khắc phục những hạn chế trên, hướng nghiên cứu tìm kiếm hoạt chất mới từ thiên nhiên nhận được sự quan tâm lớn. Alkaloid là nhóm hợp chất tự nhiên đa dạng với hơn 18.000 cấu trúc đã được công bố, sở hữu nhiều hoạt tính sinh học đa dạng như giảm đau, kích thích hoặc ức chế thần kinh trung ương, hạ huyết áp và kháng u. Mặc dù các nghiên cứu in vivo trước đây từng ghi nhận một số vinca alkaloid có khả năng làm giảm biểu hiện protein ERα sau 12 giờ điều trị, các nghiên cứu chuyên sâu ở cấp độ phân tử nhằm giải thích cơ chế tương tác và ái lực của khung alkaloid với thụ thể ERα vẫn còn rất hạn chế.

Từ thực tiễn đó, khóa luận 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 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ấu trúc tinh thể 3D của thụ thể ERα ở người (mã PDB: 3ERT, độ phân giải 1,9 Å) phức hợp với chất đối kháng đồng kết tinh 4-hydroxytamoxifen (OHT); cùng cơ sở dữ liệu gồm 52 hợp chất alkaloid có hoạt tính kháng ung thư thu thập từ thư viện hóa học PubChem.
  • Phạm vi nghiên cứu: Nghiên cứu tính toán mô phỏng in silico, tập trung vào docking phân tử, sàng lọc quy tắc Lipinski và dự đoán các thông số dược động học - độc tính (ADMET). 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, hoàn thành năm 2022.

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

Khóa luận vận dụng hệ thống lý thuyết về cấu trúc sinh học phân tử của thụ thể ERα và các phương pháp tính toán trong phát minh thuốc dựa trên cấu trúc (Structure-Based Drug Design - SBDD):

  • Cấu trúc thụ thể ERα: Thụ thể có 3 miền chính gồm miền đầu N (NTD chứa vùng hoạt hóa AF1 độc lập với hormon), miền liên kết DNA (DBD), và miền liên kết phối tử (LBD chứa vùng hoạt hóa AF2 phụ thuộc hormon). Cấu trúc miền LBD hình cầu bao gồm 11 vòng xoắn alpha (H1–H12), tạo nên một khoang kỵ nước hình elip được bao quanh bởi 22 gốc acid amin. Các acid amin chủ chốt định hướng phối tử trong túi liên kết gồm Glu353 (trên xoắn H3), Arg394 (trên xoắn H5) và His524 (trên xoắn H11).
  • Quy tắc Lipinski (Rule of 5): Bộ tiêu chuẩn đánh giá đặc tính giống thuốc dùng đường uống, quy định một hợp chất không được vi phạm quá 1 trong 4 tiêu chí: 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.
  • Hệ thống đánh giá ADMET: Dự đoán các đặc tính dược động học gồm hấp thu (tính tan trong nước, tính thấm màng Caco-2, mức độ hấp thu ở ruột người HIA, cơ chất/chất ức chế P-glycoprotein), phân bố (thể tích phân bố VDss, tính thấm qua hàng rào máu não BBB, tính thấm hệ thần kinh trung ương CNS), chuyển hóa (vai trò cơ chất và khả năng ức chế các isozym cytochrome P450: CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4), thải trừ (độ thanh thải toàn phần), và độc tính (độc tính gây đột biến gen AMES, ức chế kênh kali hERG I/II, liều chết trung bình cấp tính đường uống LD50 ở chuột, độc tính bán trường diễn LOAEL, độc tính trên gan và kích ứng da).

Phương pháp nghiên cứu thực nghiệm trên máy tính được triển khai 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 (Re-docking): Tách phối tử đồng kết tinh 4-hydroxytamoxifen ra khỏi phức hợp 3ERT bằng BIOVIA Discovery Studio 2021 Client. Chuẩn bị protein bằng Autodock Tools (loại bỏ phân tử nước, thêm nguyên tử hydro phân cực, gán trường lực Kollman, xuất định dạng pdbqt). Chuẩn bị ligand và tối ưu hóa hình học năng lượng bằng Avogadro (thuật toán Gradient liên hợp - Conjugate Gradients). Thực hiện docking lại phối tử vào protein bằng Autodock. Tính toán độ lệch bình phương trung bình gốc (RMSD) giữa cấu dạng sau dock và tinh thể gốc bằng UCSF Chimera; ngưỡng RMSD < 1,5 Å xác nhận tính tin cậy của mô hình. Năng lượng tự do liên kết ($\Delta G$) của 4-hydroxytamoxifen được lấy làm điểm chuẩn sàng lọc.
    • Mô phỏng sàng lọc 52 alkaloid: Thiết lập kích thước hộp tìm kiếm (grid box) là $40 \times 40 \times 40$ Å, khoảng cách các ô lưới là 0,375 Å, tọa độ tâm đặt tại $x = 30,282$; $y = -1,913$; $z = 24,207$ tương ứng với vị trí khoang hoạt động của ERα. Sàng lọc các phối tử có $\Delta G \le -9,8$ kcal/mol, RMSD thấp và tạo liên kết tốt với các acid amin đích.
  • Bước 2: Sàng lọc đặc điểm giống thuốc theo quy tắc Lipinski
    • Nhập mã cấu trúc phân tử SMILES của các hợp chất vượt qua bước docking vào công cụ máy chủ tính toán pkCSM để trích xuất các thông số MW, HBD, HBA và LogP.
  • Bước 3: Dự đoán các thông số ADMET
    • Sử dụng nền tảng tính toán pkCSM (dựa trên mô hình học máy và chữ ký đồ thị phân tử) để trích xuất toàn diện các chỉ số dược động học và an toàn sinh học của các hợp chất đạt tiêu chuẩn Lipinski. Phân tích tương tác 2D và 3D tại vị trí gắn kết bằng BIOVIA Discovery Studio 2021 Client.

Nguồn dữ liệu và trang thiết bị thực hiện:

  • Dữ liệu cấu trúc thụ thể: Ngân hàng dữ liệu protein quốc tế Protein Data Bank (PDB ID: 3ERT).
  • Dữ liệu cấu trúc phối tử: Thư viện hóa chất PubChem (NCBI, Hoa Kỳ).
  • Thiết bị phần cứng: Máy tính Dell Vostro 15, chạy hệ điều hành Windows 10.
  • Hệ thống phần mềm: MGL Tools 1.5.6/1.5.7, Autodock Vina 1.1.2 / Autodock, UCSF Chimera 1.14, Avogadro 1.2.0, BIOVIA Discovery Studio 2021 Client, máy chủ trực tuyến pkCSM, Microsoft Office 2020.

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ác kiến thức nền tảng về bệnh học ung thư vú, cấu trúc sinh học phân tử của thụ thể ERα, các thuốc điều trị hiện hành, nhóm hợp chất tự nhiên alkaloid và các phương pháp nghiên cứu in silico. Tác giả phân tích giải phẫu bệnh ung thư vú, nêu rõ 85% bắt nguồn từ biểu mô ống và 15% từ tiểu thùy, cùng các yếu tố nguy cơ không thể điều chỉnh (giới tính, tuổi tác, tiền sử gia đình, đột biến gen) và có thể điều chỉnh (béo phì, lối sống tĩnh tại, rượu bia). Về mặt thụ thể, ERα được trình bày chi tiết về cấu trúc từ đầu N đến đầu C gồm 3 miền chính (NTD, DBD, LBD), vị trí các vùng AF1, AF2 và miền bản lề D. Tác giả phân loại cơ chế tác dụng của các thuốc nội tiết hiện có (SERMs, SERDs, AIs) cùng các hạn chế lâm sàng tương ứng như nguy cơ ung thư nội mạc tử cung của tamoxifen và toremifen, hoặc biến cố huyết khối tĩnh mạch của raloxifen. Đối với alkaloid, chương tổng quan phân loại theo cấu trúc (dị vòng và không dị vòng) và nguồn gốc tự nhiên, đồng thời nhấn mạnh hoạt tính của các vinca alkaloid trên thụ thể ERα. Phần cuối chương trình bày nguyên lý phát minh thuốc nhờ máy tính (CADD/SBDD), cơ chế tính điểm năng lượng tự do Gibbs trong docking phân tử, định nghĩa 4 chỉ số của quy tắc Lipinski và vai trò của việc sàng lọc sớm ADMET nhằm giảm thiểu tỷ lệ thất bại của thuốc trong 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 2 mô tả chi tiết nguồn nguyên liệu, công cụ phần mềm, thiết bị phần cứng và quy trình thực nghiệm mô phỏng. Đối tượng gồm cấu trúc thụ thể ERα người gắn 4-hydroxytamoxifen (PDB ID: 3ERT) và 52 phân tử alkaloid từ PubChem. Tác giả trình bày phương pháp chuẩn bị protein (loại bỏ nước và ligand, bổ sung hydro phân cực, gán điện tích Kollman), phương pháp chuẩn bị phối tử (chuyển đổi cấu trúc, tối ưu hóa năng lượng trường lực bằng Avogadro), và các thông số thiết lập cho hộp lưới tính toán docking (kích thước $40 \times 40 \times 40$ Å, khoảng cách mắt lưới 0,375 Å tại tâm hoạt động). Quy trình thẩm định re-docking được chuẩn hóa với tiêu chí RMSD < 1,5 Å. Chương này cũng quy định rõ các tiêu chuẩn loại trừ và chọn lọc: các phối tử phải có điểm năng lượng liên kết $\Delta G \le -9,8$ kcal/mol (ngang bằng hoặc thấp hơn chất đối chứng 4-hydroxytamoxifen); vượt qua bộ lọc Lipinski (không vi phạm quá 1 tiêu chí); và đạt các chỉ số dược động học tối ưu trên mô hình pkCSM (HIA > 70%, Caco-2 > 0,9, logBB < 0,3, không gây độc tính AMES hay ức chế kênh hERG).

Chương 3: Kết quả

Chương 3 trình bày toàn bộ số liệu thu được qua các bước sàng lọc tính toán:

  • Kết quả thẩm định mô hình: Quá trình re-dock phân tử đối chứng 4-hydroxytamoxifen vào thụ thể 3ERT đạt giá trị RMSD = 1,15 Å (nhỏ hơn ngưỡng 1,5 Å), cấu dạng tái docking trùng khít với cấu dạng trong tinh thể thực nghiệm. Năng lượng liên kết đạt $\Delta G = -9,8$ kcal/mol. Phối tử đối chứng hình thành liên kết hydro với Glu353, Arg394 và tương tác alkyl/$\pi$-alkyl với Ala350, Leu387, Leu391, Leu525, Met421.
  • Kết quả sàng lọc docking 52 alkaloid: Năng lượng liên kết của 52 chất dao động từ -6,3 kcal/mol đến -10,8 kcal/mol. Có 4 hợp chất đạt năng lượng liên kết $\Delta G \le -9,8$ 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,2 kcal/mol, trong bảng phân tích tương tác ghi nhận mức năng lượng liên kết cực tiểu đạt -10,8 kcal/mol).
  • Kết quả sàng lọc quy tắc Lipinski: Được biểu diễn chi tiết qua bảng số liệu thực nghiệm:
STT Tên hợp chất PubChem ID Phân tử khối (MW, Da) Số nhóm cho H (HBD) Số nhóm nhận H (HBA) LogP Đánh giá Lipinski
1 Penitrem A 6610243 634,44 3 6 4,54 Vi phạm 1 tiêu chí (MW > 500), đạt 2/4 chỉ số chung (bị loại)
2 Pityriacitrin B 56649543 374,39 2 4 3,60 Đạt 4/4 tiêu chí (thỏa mãn)
3 Penitrem B 3084087 583,76 2 5 5,21 Vi phạm 2 tiêu chí (MW > 500, LogP > 5) (bị loại)
4 Tubocuraran 9548771 580,71 2 6 4,81 Vi phạm 1 tiêu chí (MW > 500), đạt 3/4 chỉ số (thỏa mãn)
  • Kết quả phân tích ADMET của Pityriacitrin B và Tubocuraran:
Thuộc tính ADMET Thông số cụ thể Pityriacitrin B Tubocuraran
Hấp thu Độ tan trong nước (log mol/L) -2,641 -3,812
Tính thấm màng Caco-2 ($\log P_{\text{app}}$, $10^{-6}$ cm/s) 0,749 0,828
Hấp thu ở ruột người HIA (%) 72,002 95,142
Tính thấm qua da (log Kp) -2,735 -2,751
Cơ chất của P-glycoprotein Có Có
Ức chế P-glycoprotein I / II Không / Không Có / Có
Phân bố Thể tích phân bố VDss (log L/kg) -1,267 -0,421
Tính thấm hàng rào máu não (log BB) -0,029 -0,512
Tính thấm hệ TKTW (log PS) -1,882 -1,388
Chuyển hóa Cơ chất của CYP2D6 / CYP3A4 Không / Không Không / Có
Ức chế CYP1A2 / CYP2C19 / CYP2C9 Không / Không / Không Có / Không / Không
Ức chế CYP2D6 / CYP3A4 Không / Không Có / Không
Thải trừ Độ thanh thải toàn phần (log mL/min/kg) 0,780 0,439
Độc tính Độc tính đột biến gen AMES Không Không
Ức chế kênh hERG I / hERG II Không / Không Không / Có
Độc tính cấp đường uống chuột $LD_{50}$ (mol/kg) 2,136 2,458
Độc tính gan / Kích ứng da Có / Không Có / Không
  • Kết quả tương tác phối tử - thụ thể: Pityriacitrin B tạo liên kết hydro với His356, Val355 và tương tác cộng hóa trị/kỵ nước với Ile358, Trp360. Tubocuraran liên kết với các gốc Glu353, Arg394, Ala350, Leu387, Leu391, Leu525, Met421 và tạo tương tác $\pi$-$\pi$ với Tyr526.

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

Tác giả tiến hành so sánh, biện luận chuyên sâu về ái lực liên kết và hồ sơ an toàn sinh học của hai ứng viên hàng đầu:

  • Pityriacitrin B: Là một alkaloid chứa nhân indole có khả năng hấp thụ tia cực tím, phân lập ban đầu từ nấm men Malassezia furfur. Dù các nghiên cứu trước đây chỉ ghi nhận hoạt tính chống tăng sinh yếu trên dòng tế bào ung thư vú và tuyến tiền liệt, kết quả docking cho thấy chất này gắn kết tốt với ERα ($\Delta G = -9,8$ kcal/mol). Về dược động học, chất có khả năng hấp thu tại ruột tốt (72-73%), tính thấm qua hàng rào máu não kém (logBB = -0,029 < 0,3) giúp hạn chế tác dụng phụ trên thần kinh trung ương, thải trừ nhanh qua thận, không phải là cơ chất và không ức chế hai enzym chuyển hóa chủ lực CYP3A4 và CYP2D6 (tránh tương tác thuốc). Tuy nhiên, nhược điểm là tính thấm màng Caco-2 chưa đạt mức tối ưu (< 0,9) và mô hình cảnh báo nguy cơ độc tính trên gan.
  • Tubocuraran: Sở hữu ái lực liên kết rất mạnh ($\Delta G = -10,2$ đến -10,8 kcal/mol), hấp thu ruột cao (95%), nhưng bộc lộ nhiều hạn chế về an toàn: khối lượng phân tử lớn (580,71 Da), ức chế kênh hERG II gây nguy cơ kéo dài khoảng QT và độc tính trên tim, chuyển hóa qua CYP3A4 và ức chế CYP2D6 (dễ gây tương tác chuyển hóa làm giảm hiệu lực thuốc phối hợp), đồng thời có dự đoán độc tính trên gan.
  • Biện luận phương pháp: Tác giả chỉ ra ưu thế của sàng lọc ảo trong việc tiết kiệm chi phí và thời gian, nhưng cũng thảo luận thẳng thắn về hạn chế của phương pháp docking tĩnh khi chưa phản ánh được sự biến đổi cấu hình linh động của protein trong môi trường sinh lý thực tế.

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

Khóa luận đã trích xuất và định lượng được các kết quả cụ thể:

  • Thẩm định thành công quy trình docking phân tử trên thụ thể ERα (PDB ID: 3ERT) với giá trị sai số RMSD = 1,15 Å, xác lập mốc năng lượng liên kết chuẩn của 4-hydroxytamoxifen là -9,8 kcal/mol.
  • Sàng lọc trong số 52 hợp chất alkaloid từ PubChem, xác định được 4 hợp chất có năng lượng liên kết vượt trội hoặc tương đương đối chứng: Tubocuraran (-10,2 đến -10,8 kcal/mol), Penitrem B (-10,1 kcal/mol), Penitrem A (-9,8 kcal/mol) và Pityriacitrin B (-9,8 kcal/mol).
  • Đánh giá quy tắc 5 Lipinski loại trừ Penitrem A và Penitrem B do vi phạm kích thước phân tử và tính kỵ nước; xác định Pityriacitrin B thỏa mãn toàn bộ 4/4 tiêu chí và Tubocuraran đạt 3/4 tiêu chí.
  • Thiết lập hồ sơ ADMET chi tiết, chỉ ra Pityriacitrin B là ứng cử viên tiềm năng nhất nhờ hấp thu đường ruột tốt (72,002%), thanh thải nhanh (0,78 log mL/min/kg), không ức chế enzym gan CYP2D6/CYP3A4 và không mang độc tính tim mạch qua kênh hERG.
  • Tác giả đề xuất định hướng nghiên cứu tối ưu hóa cấu trúc hóa học của dẫn xuất Pityriacitrin B nhằm cải thiện tính thấm tế bào (tăng chỉ số Caco-2) và loại bỏ cảnh báo độc tính trên gan.
  • Về mặt học thuật, công trình đóng góp dữ liệu mô phỏng tương tác phân tử đầu tiên của Pityriacitrin B trên thụ thể ERα, cung cấp cơ sở cấu trúc giải thích khả năng ức chế thụ thể này của các alkaloid tự nhiên, bổ sung bằng chứng phân tử cho các quan sát in vivo trước đây về tác dụng kháng ung thư vú của nhóm alkaloid.

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

Khóa luận nêu rõ các giới hạn nội tại bắt nguồn từ mô hình tính toán và phạm vi nghiên cứu:

  • Toàn bộ kết quả đều được xây dựng trên các mô hình dự đoán in silico, hoàn toàn chưa có thử nghiệm thực nghiệm sinh học (in vitro trên dòng tế bào hoặc enzyme tinh chế, in vivo trên động vật) để kiểm chứng ái lực liên kết thực tế và nồng độ ức chế $IC_{50}$.
  • Mô hình docking phân tử sử dụng thụ thể dạng tĩnh (rigid receptor), chưa tính đến độ linh động cấu hình và dung môi nước đầy đủ như trong các mô phỏng động lực học phân tử (Molecular Dynamics simulation).
  • Hợp chất Pityriacitrin B vẫn còn tồn tại các cảnh báo về tính thấm màng Caco-2 thấp hơn ngưỡng tiêu chuẩn (0,749 so với mức khuyến nghị > 0,9) và nguy cơ độc tính trên gan.
  • Hướng nghiên cứu tiếp theo được tác giả đề xuất: Thực hiện mô phỏng động lực học phân tử (MD) để khảo sát độ bền phức hợp theo thời gian; tiến hành bán tổng hợp hoặc biến đổi cấu trúc phân tử Pityriacitrin B nhằm giảm độc tính tế bào gan; triển khai các thử nghiệm sinh học in vitro đánh giá khả năng ức chế ERα trên dòng tế bào ung thư vú (như MCF-7).

Giá trị tham khảo

Khóa luận là tài liệu học thuật có giá trị tham khảo 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:

  • Cung cấp quy trình từng bước chuẩn hóa về sàng lọc ảo dựa trên cấu trúc (SBDD), từ cách xử lý file tinh thể PDB, quy trình re-docking thẩm định thông số, đến phương pháp sử dụng chuỗi công cụ phần mềm Autodock, UCSF Chimera, Avogadro và BIOVIA Discovery Studio.
  • Cung cấp phương pháp luận ứng dụng máy chủ pkCSM để đánh giá quy tắc Lipinski và dự đoán các thông số dược động học - độc tính ADMET, phục vụ định hướng sàng lọc các thư viện hoạt chất tự nhiên.
  • Nội dung Chương 3 và Chương 4 là nguồn dữ liệu tham khảo chuyên biệt cho các đề tài nghiên cứu hoạt tính kháng ung thư vú hướng đích thụ thể estrogen alpha từ nguồn tài nguyên alkaloid thực vật và vi nấm.

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

1. Thụ thể ER-alpha được lựa chọn trong nghiên cứu lấy từ nguồn nào và ligand đối chứng là gì?
Cấu trúc tinh thể 3 chiều của thụ thể ERα được tải về từ Ngân hàng dữ liệu protein (Protein Data Bank) với mã PDB là 3ERT, ở độ phân giải 1,9 Å. Cấu trúc này mô tả thụ thể ERα ở người được gắn kết với ligand đối kháng nội tiết chọn lọc đồng kết tinh là 4-hydroxytamoxifen (OHT).

2. Tiêu chuẩn nào được sử dụng để chứng minh quy trình docking phân tử là hợp lý và đáng tin cậy?
Quy trình docking được thẩm định thông qua việc tái docking (re-dock) phân tử 4-hydroxytamoxifen vào lại vị trí hoạt động của thụ thể ERα. Kết quả tính toán bằng phần mềm UCSF Chimera cho thấy độ lệch bình phương trung bình gốc (RMSD) giữa cấu dạng sau khi dock và cấu dạng thực nghiệm trong tinh thể gốc là 1,15 Å. Do giá trị này nhỏ hơn ngưỡng quy ước 1,5 Å, quy trình docking được xác nhận đạt độ tin cậy khoa học.

3. Bốn hợp chất alkaloid nào đạt điểm năng lượng liên kết tốt nhất trong số 52 chất được sàng lọc?
Bốn hợp chất có năng lượng liên kết $\Delta G \le -9,8$ kcal/mol (ngang bằng hoặc thấp hơn mốc đối chứng của 4-hydroxytamoxifen) gồm có: Penitrem A (PubChem ID: 6610243, $\Delta G = -9,8$ kcal/mol), Pityriacitrin B (PubChem ID: 56649543, $\Delta G = -9,8$ kcal/mol), Penitrem B (PubChem ID: 3084087, $\Delta G = -10,1$ kcal/mol), và Tubocuraran (PubChem ID: 9548771, $\Delta G = -10,2$ kcal/mol, điểm tương tác cực tiểu đạt -10,8 kcal/mol).

4. Vì sao Penitrem A và Penitrem B bị loại bỏ dù có năng lượng liên kết docking rất tốt?
Penitrem A và Penitrem B bị loại sau bước sàng lọc đặc điểm giống thuốc do vi phạm quy tắc Lipinski: Penitrem A có khối lượng phân tử lớn (MW = 634,44 Da > 500 Da) và chỉ đạt 2/4 tiêu chí chung; Penitrem B vi phạm cùng lúc 2 tiêu chí là khối lượng phân tử vượt ngưỡng (MW = 583,76 Da > 500 Da) và hệ số phân bố octanol/nước quá cao (LogP = 5,21 > 5).

5. Vì sao Pityriacitrin B được đánh giá triển vọng hơn Tubocuraran dù Tubocuraran có điểm docking thấp hơn?
Mặc dù Tubocuraran có năng lượng liên kết thấp hơn (-10,2 đến -10,8 kcal/mol so với -9,8 kcal/mol), hợp chất này vi phạm quy tắc Lipinski (MW = 580,71 Da) và mang nhiều nguy cơ dược động học bất lợi: ức chế kênh hERG II gây độc tính tim mạch, là cơ chất của CYP3A4 và ức chế CYP2D6 gây nguy cơ tương tác thuốc. Ngược lại, Pityriacitrin B đạt cả 4 tiêu chuẩn Lipinski, không ức chế các enzym chuyển hóa gan quan trọng, không có độc tính tim mạch hERG và không gây đột biến gen AMES.

Kết luận

Khóa luận tốt nghiệp của tác giả Bùi Phúc Lộc đã hoàn thành việc sàng lọc in silico 52 hợp chất alkaloid tự nhiên trên đích tác dụng thụ thể estrogen alpha (ERα, mã PDB 3ERT) nhằm tìm kiếm hướng điều trị ung thư vú. Thông qua quy trình tích hợp giữa mô phỏng docking phân tử, bộ lọc quy tắc Lipinski và dự đoán dược động học ADMET trên nền tảng pkCSM, nghiên cứu đã loại bỏ các phân tử có nguy cơ độc tính cao và xác định được Pityriacitrin B là hợp chất tiềm năng nhất với năng lượng liên kết $\Delta G = -9,8$ kcal/mol và các chỉ số an toàn sinh học khả quan. Đề tài cung cấp cơ sở dữ liệu cấu trúc và phương pháp luận tính toán có giá trị cho việc tiếp tục tối ưu hóa phân tử và triển khai các thử nghiệm sinh học thực nghiệm tiếp theo.