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

Ung thư vú là một trong những bệnh lý ác tính có tỷ lệ mắc và tử vong hàng đầu ở phụ nữ trên toàn cầu cũng như tại Việt Nam. Theo thống kê của tổ chức Ghi nhận Ung thư Toàn cầu (Globocan 2020), toàn thế giới ghi nhận khoảng 2,3 triệu ca mắc mới (chiếm 11,7% tổng số ca ung thư các loại) và 685.000 ca tử vong. Tại Việt Nam, ung thư vú đứng đầu về số ca mắc mới ở nữ giới với 21.555 ca (chiếm 25,8% tổng số ca ung thư ở nữ) và là nguyên nhân tử vong do ung thư đứng hàng thứ tư với 9.345 ca.

Thụ thể estrogen alpha (ERα) là một đích sinh học trọng yếu trong điều trị ung thư vú, biểu hiện ở khoảng 50% đến 80% các trường hợp khối u vú và đóng vai trò thúc đẩy sự hình thành, tiến triển của tế bào ác tính. Các thuốc ức chế ERα hiện hành trên lâm sàng như tamoxifen, raloxifene, toremifene hoặc fulvestrant vẫn bộc lộ nhiều hạn chế đáng kể: xuất hiện hiện tượng kháng thuốc sau thời gian dài điều trị, nguy cơ gây ung thư nội mạc tử cung và gia tăng biến cố thuyên tắc huyết khối tĩnh mạch sâu. Do đó, việc tìm kiếm các hợp chất mới từ nguồn gốc tự nhiên có khả năng ức chế chọn lọc ERα với hoạt tính cao và giảm thiểu tác dụng không mong muốn là một nhu cầu cấp thiết.

Nhóm hợp chất alkaloid từ tự nhiên sở hữu hơn 18.000 cấu trúc đã được xác định, mang nhiều hoạt tính sinh học đa dạng. Mặc dù một số thử nghiệm in vivo đã bước đầu ghi nhận khả năng làm giảm mức biểu hiện protein ERα của các alkaloid (như vinca alkaloid sau 12 giờ xử lý), song các dữ liệu nghiên cứu ở cấp độ phân tử nhằm giải thích cơ chế tương tác trực tiếp giữa khung cấu trúc alkaloid và thụ thể ERα vẫn còn hạn chế. Nhằm giải quyết vấn đề này, 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 được thực hiện với 2 mục tiê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. 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: Cấu trúc 3D của thụ thể ER-alpha người (mã PDB: 3ERT) và cơ sở dữ liệu cấu trúc hóa học của 52 hợp chất alkaloid tự nhiên lấy từ cơ sở dữ liệu PubChem.
  • Phạm vi nghiên cứu: Mô phỏng in silico (docking phân tử, đánh giá quy tắc Lipinski và dự đoán ADMET) được tiến hành tại Trường Đại học Y Dược – Đại học Quốc gia Hà Nội, công bố năm 2022.

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

Nghiên cứu ứng dụng phương pháp phát minh thuốc có sự hỗ trợ của máy tính (CADD), cụ thể là tiếp cận phát minh thuốc dựa trên cấu trúc (SBDD), kết hợp sàng lọc ảo (virtual screening), đánh giá đặc tính giống thuốc và dự đoán dược động học – độc tính:

  • Khung lý thuyết về đích sinh học: Thụ thể ERα thuộc siêu họ thụ thể nhân, hoạt động như một yếu tố phiên mã điều hòa biểu hiện gen mục tiêu. Vùng liên kết phối tử (LBD) của ERα có cấu trúc hình cầu gồm 11 vòng xoắn alpha (H1 – H12), tạo nên một khoang kỵ nước ellipsoid. Phối tử tự nhiên 17β-estradiol tương tác thông qua các liên kết hydro đặc hiệu với các acid amin Glu353 (trên xoắn H3), Arg394 (trên xoắn H5) và His524 (trên xoắn H11).
  • Phương pháp Docking phân tử: Dự đoán cấu hình liên kết tối ưu và tính toán năng lượng tự do liên kết ($\Delta G$, tính bằng kcal/mol) dựa trên thuật toán tìm kiếm và hàm tính điểm năng lượng (xét tương tác Van der Waals, liên kết hydro, tương tác tĩnh điện và solvat hóa). Quy trình được thẩm định tính chính xác bằng kỹ thuật re-dock phối tử đồng kết tinh 4-hydroxytamoxifen (OHT) với ngưỡng sai số bình phương trung bình gốc RMSD < 1,5 Å.
  • Quy tắc 5 Lipinski (Rule of 5 - RO5): Đánh giá tiềm năng hấp thu đường uống của các hợp chất dựa trên 4 tiêu chí 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) và hệ số phân bố octanol/nước (LogP < 5). Hợp chất không vi phạm quá 1 tiêu chí được xem là đạt đặc tính giống thuốc.
  • Dự đoán ADMET in silico: Sử dụng nền tảng pkCSM dựa trên đồ thị phân tử để dự đoán các chỉ số dược động học: độ tan, tính thấm qua màng tế bào Caco-2, tỷ lệ hấp thu qua ruột người (HIA), 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), khả năng chuyển hóa qua hệ enzyme Cytochrome P450 (CYP), độ thanh thải toàn phần và các chỉ số độc tính (thử nghiệm đột biến gen Ames, độc tính trên gan, ức chế kênh hERG I/II, liều gây chết trung bình LD50 đường uống trên chuột).
  • Nguồn dữ liệu và công cụ phần mềm:
    • Dữ liệu: Ngân hàng Dữ liệu Protein (Protein Data Bank - RCSB PDB) lấy cấu trúc 3ERT; Cơ sở dữ liệu hóa học PubChem lấy cấu trúc các alkaloid và phối tử đối chứng.
    • Phần mềm & Nền tảng: AutoDock / AutoDock Tools (MGL Tools 1.5.6), Autodock Vina, UCSF Chimera 1.14, Discovery Studio 2021 Client, Avogadro, nền tảng trực tuyến pkCSM, phần cứng máy tính Dell Vostro 15 (Windows 10).

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 y văn trong và ngoài nước về bệnh học ung thư vú, phân tích chi tiết dịch tễ học toàn cầu và tại Việt Nam dựa trên báo cáo Globocan 2020, mô tả cấu trúc giải phẫu tuyến vú, cơ chế bệnh sinh liên quan đến vi môi trường khối u và tế bào gốc ung thư (CSC). Tác giả hệ thống hóa các phương pháp điều trị ung thư vú hiện hành gồm phẫu thuật (cắt tuyến vú vét hạch, phẫu thuật bảo tồn, tái tạo), xạ trị (chiếu ngoài, áp sát, xạ trị trong mổ - IORT), hóa trị liệu và liệu pháp nội tiết.

Về mặt hóa dược và sinh học phân tử, chương này làm rõ cấu trúc không gian của thụ thể ERα gồm 3 miền chính: NTD (chứa vùng AF1), DBD (miền gắn DNA) và LBD (chứa vùng AF2 và khoang gắn phối tử). Tác giả phân loại các nhóm thuốc tác động lên ERα (chất ức chế Aromatase, SERD, SERM) và nêu rõ các hạn chế lâm sàng của tamoxifen, raloxifene. Đồng thời, chương 1 cung cấp cái nhìn tổng quan về hóa học alkaloid (phân loại alkaloid dị vòng và không dị vòng, đặc tính tạo liên kết hydro qua nguyên tử nitơ), các cơ sở toán học của nghiên cứu in silico, các thuật toán docking (so khớp, xây dựng tăng dần, Monte Carlo, thuật toán di truyền GA, động lực học phân tử MD), ý nghĩa của hàm tính điểm và các quy tắc dự đoán ADMET/Lipinski.

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

Chương 2 trình bày chi tiết quy trình thực nghiệm in silico được chia thành 3 bước logic:

  1. Chuẩn bị cấu trúc và thẩm định docking: Tải cấu trúc 3ERT (độ phân giải 1,9 Å) từ PDB, sử dụng Discovery Studio 2021 Client để tách phối tử đồng kết tinh 4-hydroxytamoxifen, loại bỏ phân tử nước và ion nền. Cấu trúc protein được bổ sung nguyên tử hydro, tối ưu hóa hydro phân cực, gán trường lực Kollman và chuyển sang định dạng .pdbqt bằng AutoDock Tools. Cấu trúc 52 phối tử alkaloid từ PubChem được tối ưu hóa năng lượng hình học bằng phần mềm Avogadro thông qua phương pháp Gradient liên hợp (Conjugate Gradients). Quy trình re-dock được thực hiện trên AutoDock với hộp tìm kiếm (grid box) có kích thước 40 × 40 × 40 Å, khoảng cách mắt lưới 0,375 Å, tọa độ tâm đặt tại $x = 30,282$; $y = -1,907$; $z = 24,207$ (trung tâm hoạt động của ERα).
  2. Sàng lọc Lipinski: Sử dụng chuỗi mã hóa SMILES của các hợp chất đạt điểm docking cao để phân tích 4 thông số của quy tắc Lipinski trên máy chủ pkCSM.
  3. Phân tích ADMET: Đánh giá định lượng và định tính các thông số dược động học và độc tính qua pkCSM nhằm chọn lọc các ứng viên tối ưu.

Chương 3: Kết quả

Chương 3 trình bày các số liệu thực nghiệm in silico thu được:

  • Thẩm định quy trình docking: Kết quả re-dock phối tử 4-hydroxytamoxifen cho giá trị RMSD giữa cấu hình kết tinh và cấu hình tái mô phỏng là 1,05 Å (đạt tiêu chuẩn < 1,5 Å), với năng lượng tự do liên kết $\Delta G = -9,8\text{ kcal/mol}$. Phối tử tái lập các tương tác then chốt gồm liên kết hydro với Glu353, Arg394 và tương tác alkyl/$\pi$-alkyl với Ala350, Leu387, Leu391, Leu525, Met421. Mức năng lượng $-9,8\text{ kcal/mol}$ được xác lập làm ngưỡng chuẩn để sàng lọc.
  • Kết quả sàng lọc docking 52 alkaloid: Đã xác định được 4 hợp chất có năng lượng liên kết $\Delta G \le -9,8\text{ kcal/mol}$:
STT Tên hợp chất PubChem CID Năng lượng liên kết ($\Delta G$, kcal/mol)
1 Tubocuraran 9548771 -10,3
2 Penitrem B 3084087 -10,1
3 Penitrem A 6610243 -9,8
4 Pityriacitrin B 56649543 -9,8
Đối chứng 4-hydroxytamoxifen 53788270 -9,8
  • Kết quả sàng lọc đặc điểm giống thuốc (Lipinski Rule of 5):
STT Tên hợp chất Khối lượng phân tử (MW, Da) Số nhóm cho H (HBD) Số nhóm nhận H (HBA) LogP Số tiêu chí vi phạm Đánh giá
1 Penitrem A 634,20 3 6 5,44 2 Loại
2 Penitrem B 584,20 2 5 6,10 2 Loại
3 Pityriacitrin B 376,37 2 4 2,87 0 Đạt
4 Tubocuraran 582,74 1 6 4,44 1 Đạt

Hai hợp chất Penitrem A và Penitrem B bị loại do vi phạm 2 tiêu chí (MW > 500 và LogP > 5). Pityriacitrin B và Tubocuraran được lựa chọn để phân tích ADMET.

  • Kết quả dự đoán thông số ADMET qua pkCSM:
Thông số ADMET Pityriacitrin B Tubocuraran
Hấp thu
Độ tan trong nước ($\log\text{ mol/L}$) -3,407 -2,641
Tính thấm Caco-2 ($\log P_{\text{app}}$, $10^{-6}\text{ cm/s}$) 0,828 0,214
Hấp thu ở ruột người - HIA (%) 93,684 72,002
Tính thấm qua da ($\log K_p$) -2,735 -2,850
Cơ chất 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\text{ L/kg}$) -1,267 -0,442
Tính thấm hàng rào máu não ($\log BB$) 0,052 -0,029
Tính thấm hệ thần kinh trung ương ($\log PS$) -1,388 -2,125
Chuyển hóa
Cơ chất CYP2D6 / CYP3A4 Không / Không Không / Có
Ức chế CYP1A2 Không Có
Ức chế CYP2C19 / CYP2C9 Không / Không 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\text{ mL/min/kg}$) 0,439 0,612
Độ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}$ ($\text{mol/kg}$) 2,502 2,136
Độc tính trên gan (Hepatotoxicity) Không Có

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

Chương 4 tập trung phân tích ý nghĩa của các kết quả tính toán:

  • Về ái lực liên kết: Cả Tubocuraran ($-10,3\text{ kcal/mol}$) và Pityriacitrin B ($-9,8\text{ kcal/mol}$) đều thể hiện năng lượng liên kết tương đương hoặc vượt trội so với chất ức chế chuẩn 4-hydroxytamoxifen ($-9,8\text{ kcal/mol}$). Phân tích tương tác không gian cho thấy các hợp chất này gắn kết chặt chẽ vào túi gắn LBD của ERα thông qua mạng lưới liên kết hydro và tương tác kỵ nước với các acid amin thiết yếu.
  • Về tính chất dược động học và độ an toàn: Pityriacitrin B biểu hiện các ưu điểm vượt trội về mặt hóa lý với khả năng hấp thu tại ruột người rất cao (93,684%), không ức chế các enzyme chuyển hóa chính của gan (CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4), không gây độc tính trên gan, không gây đột biến gen và không ức chế kênh kali tim hERG I/II. Ngược lại, Tubocuraran dù có ái lực gắn kết cao hơn nhưng lại vi phạm 1 quy tắc Lipinski (MW = 582,74 Da), có khả năng gây độc gan, ức chế kênh hERG II và ức chế các enzyme CYP1A2, CYP2D6, đòi hỏi các nghiên cứu tối ưu hóa cấu trúc sâu hơn.
  • Về phương pháp và hạn chế: Tác giả thảo luận về bản chất gần đúng của hàm tính điểm trong docking phân tử, tính chất tĩnh của cấu trúc tinh thể protein PDB và nhấn mạnh sự cần thiết của các kiểm chứng thực nghiệm sinh học.

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

  • Kết quả chính:
    1. Thẩm định thành công quy trình docking phân tử trên thụ thể ERα người (PDB: 3ERT) với giá trị RMSD re-dock đạt 1,05 Å và năng lượng liên kết đối chứng của 4-hydroxytamoxifen đạt $-9,8\text{ kcal/mol}$.
    2. Sàng lọc 52 hợp chất alkaloid tự nhiên, phát hiện 4 hợp chất có ái lực gắn kết cao ($\Delta G \le -9,8\text{ kcal/mol}$), trong đó Tubocuraran đạt $-10,3\text{ kcal/mol}$ và Pityriacitrin B đạt $-9,8\text{ kcal/mol}$.
    3. Lọc thành công 2 hợp chất thỏa mãn quy tắc 5 Lipinski và xây dựng bảng hồ sơ dược động học, độc tính (ADMET) chi tiết cho Pityriacitrin B và Tubocuraran.
  • Giải pháp và đề xuất của tác giả: Đề xuất Pityriacitrin B là ứng viên tiềm năng nhất cần được ưu tiên đưa vào các giai đoạn nghiên cứu thực nghiệm in vitro và in vivo nhằm phát triển thuốc kháng ERα mới hỗ trợ điều trị ung thư vú.
  • Đóng góp học thuật: Nghiên cứu cung cấp dữ liệu định lượng ở cấp độ phân tử về cơ chế gắn kết của các alkaloid tự nhiên lên thụ thể ERα, chứng minh tính khả thi của mô hình sàng lọc ảo in silico trong việc rút ngắn thời gian và tối ưu hóa chi phí nghiên cứu phát triển thuốc điều trị ung thư từ nguồn dược liệu.

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

  • Hạn chế được tác giả nêu rõ:
    • Nghiên cứu dừng lại ở mức độ mô phỏng tính toán in silico, chưa tiến hành thực nghiệm hóa sinh 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 như MCF-7) hoặc in vivo để kiểm chứng hoạt tính ức chế thực tế.
    • Docking phân tử sử dụng mô hình thụ thể tĩnh, chưa khảo sát được sự biến đổi cấu hình động của phức hợp protein-phối tử theo thời gian trong môi trường dung môi sinh lý.
  • Hướng nghiên cứu tiếp theo:
    • Thực hiện mô phỏng động lực học phân tử (Molecular Dynamics simulation) để đánh giá độ bền nhiệt động và tính ổn định không gian của phức hợp ERα – Pityriacitrin B.
    • Tiến hành chiết tách, bán tổng hợp hoặc tổng hợp toàn phần hợp chất Pityriacitrin B để thử nghiệm hoạt tính sinh học in vitro và in vivo trên các mô hình ung thư vú.

Giá trị tham khảo

  • Đối tượng thụ hưởng:
    • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Dược học, Hóa dược, Công nghệ Sinh học và Tin sinh học (Bioinformatics).
    • Các nhà nghiên cứu đang thực hiện đề tài sàng lọc thuốc điều trị ung thư, đặc biệt là ung thư phụ thuộc nội tiết tố (ER-positive breast cancer).
  • Nội dung giá trị nhất để tham khảo:
    • Quy trình từng bước chuẩn hóa mô phỏng docking phân tử trên AutoDock và AutoDock Tools kết hợp kiểm định RMSD.
    • Thiết lập thông số grid box và tọa độ không gian chính xác tại vùng LBD của thụ thể ERα (3ERT).
    • Phương pháp luận kết hợp sàng lọc ái lực liên kết, quy tắc Lipinski và dự đoán ADMET bằng công cụ pkCSM để thu hẹp thư viện hợp chất tự nhiên.

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

1. Phối tử đối chứng nào được sử dụng để thẩm định quy trình docking và kết quả đạt được là bao nhiêu?
Nghiên cứu sử dụng chất ức chế chọn lọc 4-hydroxytamoxifen (OHT) được tách trực tiếp từ phức hợp tinh thể thụ thể ERα (mã PDB: 3ERT). Kết quả re-dock cho giá trị RMSD là 1,05 Å (đạt tiêu chuẩn tin cậy < 1,5 Å) và năng lượng tự do liên kết $\Delta G = -9,8\text{ kcal/mol}$.

2. Trong số 52 hợp chất alkaloid khảo sát, những hợp chất nào có năng lượng liên kết tốt hơn hoặc bằng đối chứng?
Có 4 hợp chất đạt $\Delta G \le -9,8\text{ kcal/mol}$, bao gồm: Tubocuraran ($-10,3\text{ kcal/mol}$), Penitrem B ($-10,1\text{ kcal/mol}$), Penitrem A ($-9,8\text{ kcal/mol}$) và Pityriacitrin B ($-9,8\text{ kcal/mol}$).

3. Vì sao Penitrem A và Penitrem B bị loại bỏ ở giai đoạn sàng lọc đặc tính giống thuốc?
Cả hai hợp chất này đều vi phạm 2 tiêu chuẩn trong Quy tắc 5 Lipinski: Penitrem A có khối lượng phân tử $\text{MW} = 634,20\text{ Da}$ và $\text{LogP} = 5,44$; Penitrem B có khối lượng phân tử $\text{MW} = 584,20\text{ Da}$ và $\text{LogP} = 6,10$ (vượt ngưỡng cho phép $\text{MW} < 500\text{ Da}$ và $\text{LogP} < 5$). Theo quy định không được vi phạm quá 1 tiêu chí, hai chất này không đạt yêu cầu về đặc tính giống thuốc đường uống.

4. Hợp chất Pityriacitrin B thể hiện những ưu điểm dược động học và độ an toàn nào nổi bật?
Pityriacitrin B đáp ứng cả 4 tiêu chuẩn Lipinski, có tỷ lệ hấp thu qua ruột người rất cao (93,684%), độ thanh thải $0,439\log\text{ mL/min/kg}$, không gây độc tính trên gan, không gây đột biến gen (Ames âm tính), không ức chế kênh tim mạch hERG I/II và không ức chế bất kỳ enzyme chuyển hóa nào thuộc họ Cytochrome P450 (CYP1A2, 2C19, 2C9, 2D6, 3A4).

5. Tọa độ và kích thước hộp tìm kiếm (grid box) trong mô phỏng docking được thiết lập như thế nào?
Hộp tìm kiếm được cấu hình với kích thước không gian $40 \times 40 \times 40\text{ \AA}$, khoảng cách giữa các điểm lưới là $0,375\text{ \AA}$, tâm vùng tìm kiếm đặt tại tọa độ $x = 30,282$; $y = -1,907$; $z = 24,207$, bao trọn khoang liên kết phối tử LBD của cấu trúc thụ thể ERα.

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 trọn vẹn hai mục tiêu đề ra thông qua việc xây dựng và thực thi quy trình sàng lọc in silico chuẩn hóa trên thụ thể ER-alpha (PDB: 3ERT). Từ thư viện 52 hợp chất alkaloid tự nhiên, nghiên cứu đã xác định được Pityriacitrin B là ứng viên tiềm năng hàng đầu với năng lượng liên kết mạnh ($-9,8\text{ kcal/mol}$), thỏa mãn hoàn toàn quy tắc 5 Lipinski và sở hữu hồ sơ ADMET an toàn, thuận lợi cho hấp thu đường uống. Kết quả này cung cấp cơ sở khoa học và dữ liệu phân tử vững chắc cho các nghiên cứu thực nghiệm tiếp theo nhằm phát triển các thuốc ức chế ER-alpha mới có nguồn gốc từ tự nhiên phục vụ điều trị ung thư vú.