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

Ung thư là một trong những nguyên nhân gây tử vong hàng đầu trên thế giới, ghi nhận gần 10 triệu ca tử vong vào năm 2020, trong đó ung thư vú chiếm khoảng 2,3 triệu ca. Theo số liệu của Tổ chức Y tế Thế giới (WHO) năm 2019, ung thư là nguyên nhân gây tử vong trước 70 tuổi đứng thứ nhất hoặc thứ hai tại 112 trên 182 quốc gia. Tại Việt Nam, thống kê của GLOBOCAN năm 2020 chỉ ra cả nước có 182.563 ca mắc mới và 122.690 ca tử vong do ung thư; riêng ung thư vú có 21.555 ca mắc mới (chiếm 11,8% tổng số ca ung thư) và 9.345 ca tử vong. Ở nữ giới Việt Nam, ung thư vú chiếm tỷ lệ mắc cao nhất với 25,8%.

Về mặt sinh học phân tử, họ thụ thể yếu tố tăng trưởng biểu bì ErbB (gồm EGFR/HER1, ErbB2/HER2, ErbB3/HER3, ErbB4/HER4) đóng vai trò điều hòa các con đường dẫn truyền tín hiệu tăng sinh và sống sót của tế bào. Sự biểu hiện quá mức hoặc đột biến của các thụ thể ErbB, đặc biệt là HER2 (chiếm khoảng 20-30% ở các trường hợp ung thư vú giai đoạn đầu), kích hoạt mạnh mẽ các trục tín hiệu MAPK (RAS-RAF-MEK-ERK) và PI3K/AKT/mTOR. Tình trạng này liên quan chặt chẽ đến tiên lượng xấu, khả năng di căn, kháng thuốc và làm giảm tỷ lệ sống sót của bệnh nhân. Do đó, việc phát triển các tác nhân ức chế chọn lọc thụ thể HER2 là một hướng đi trọng tâm trong điều trị đích ung thư vú.

Từ góc độ dịch tễ học, các nghiên cứu ghi nhận phụ nữ tại các quốc gia châu Á như Nhật Bản và Trung Quốc có tỷ lệ mắc ung thư vú thấp hơn đáng kể so với các nước phương Tây nhờ chế độ ăn giàu đậu nành chứa các hợp chất isoflavone. Tổng hợp 23 nghiên cứu trước đây trên 330.826 người cho thấy việc tiêu thụ isoflavone tỷ lệ nghịch với tỷ lệ tử vong do ung thư; cụ thể, bổ sung 10 mg isoflavone mỗi ngày giúp giảm 7-9% nguy cơ tử vong do mọi bệnh ung thư. Isoflavone có khả năng ức chế tăng sinh tế bào, ngăn chặn hình thành mạch, kích hoạt quá trình chết theo chương trình (apoptosis) và cạnh tranh với estrogen nội sinh.

Tuy nhiên, quy trình tổng hợp hóa học và thử nghiệm hoạt tính sinh học truyền thống đòi hỏi nguồn kinh phí lớn và thời gian kéo dài. Nhằm tối ưu hóa quá trình tìm kiếm thuốc mới, phương pháp sàng lọc ảo trên máy tính (in silico), đặc biệt là kỹ thuật docking phân tử, được ứng dụng để mô phỏng tương tác và dự đoán ái lực liên kết giữa phối tử với protein đích.

Khóa luận được thực hiện nhằm giải quyết 2 mục tiêu nghiên cứu cụ thể:

  1. Sàng lọc các hợp chất isoflavone có tác dụng ức chế enzym HER2 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 (ADMET) 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 enzym HER2 (mã PDB: 3PP0) và tập dữ liệu 60 hợp chất isoflavone tự nhiên được thu thập từ y văn, đối chiếu với các thuốc đối chứng dương đã được FDA phê duyệt gồm Trastuzumab, Lapatinib và Tucatinib.
  • Phạm vi nghiên cứu: Mô phỏng tính toán in silico, sàng lọc các đặc tính hóa lý, dược động học và an toàn tế bào dựa trên mô hình thuật toán máy tính; 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 vào tháng 5 năm 2022.

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

Cơ sở lý thuyết và công cụ nghiên cứu

Nghiên cứu dựa trên cơ sở cấu trúc không gian ba chiều của miền tyrosine kinase thuộc thụ thể HER2 (ErbB2). Cấu trúc tinh thể HER2 được lấy từ Ngân hàng Dữ liệu Protein (RCSB Protein Data Bank) với mã định danh 3PP0, có độ phân giải 2,25 Å. Cấu trúc này chứa sẵn phối tử đồng kết tinh là 03Q (công thức hóa học: 2-{2-[4-({5-chloro-6-[3-(trifluoromethyl)phenoxy]pyridin-3-yl}amino)-5H-pyrrolo[3,2-d]pyrimidin-5-yl]ethoxy}ethanol).

Thư viện phối tử gồm 60 hợp chất isoflavone có nguồn gốc tự nhiên được xác định từ các công bố khoa học trước đó và tải về từ cơ sở dữ liệu PubChem dưới định dạng SDF. Ba hoạt chất kháng HER2 đã được Cục Quản lý Thực phẩm và Dược phẩm Hoa Kỳ (FDA) phê duyệt gồm Trastuzumab, Lapatinib và Tucatinib được sử dụng làm chất đối chứng.

Hệ thống thiết bị và phần mềm phục vụ nghiên cứu gồm:

  • Thiết bị phần cứng: Máy tính để bàn CPU Intel Core i5-10500, GPU NVIDIA GeForce GTX 1660, hệ điều hành Windows 11.
  • Phần mềm chuyên dụng: MGL Tools 1.5.6 (AutoDockTools), AutoDock Vina, UCSF Chimera 1.14, Discovery Studio 2021 Client, Avogadro và công cụ tính toán trực tuyến pkCSM.

Phương pháp nghiên cứu

Quy trình nghiên cứu được triển khai theo 3 bước tuần tự:

  1. Thẩm định quy trình docking (Re-docking): Tách phối tử đồng kết tinh 03Q ra khỏi phức hợp 3PP0, loại bỏ các phân tử nước, bổ sung nguyên tử hydro phân cực và gán trường lực Kollman cho protein bằng Discovery Studio và AutoDockTools. Tiến hành docking lại phối tử 03Q vào vùng hoạt động của enzym với hộp lưới (grid box) có kích thước 30 Å × 30 Å × 30 Å, tọa độ trung tâm (x, y, z) = (34, 46, -12), số vòng lặp tính toán là 8. Tính toán giá trị sai lệch bình phương trung bình gốc (RMSD) giữa cấu dạng sau khi re-dock và cấu dạng gốc bằng phần mềm UCSF Chimera; ngưỡng chấp nhận độ tin cậy của quy trình là RMSD ≤ 1,5 Å.
  2. Mô phỏng docking phân tử sàng lọc: Tối ưu hóa cấu trúc không gian 3D của 60 hợp chất isoflavone bằng thuật toán Gradient liên hợp (Conjugate Gradients) trên phần mềm Avogadro, chuyển đổi sang định dạng .pdbqt. Tiến hành docking vào vị trí hoạt động của enzym 3PP0 bằng AutoDock Vina với cùng thông số hộp lưới. Mỗi hợp chất được docking lặp lại 3 lần để lấy giá trị năng lượng tự do liên kết ($\Delta G$, tính bằng kcal/mol) trung bình. Tiêu chí sàng lọc là năng lượng liên kết phải nhỏ hơn hoặc bằng năng lượng của chất chứng có ái lực tốt nhất (Trastuzumab đạt -9,4 kcal/mol).
  3. Đánh giá đặc tính giống thuốc và dự đoán ADMET:
    • Đánh giá quy tắc 5 của Lipinski (Rule of 5 - RO5) gồm 4 tiêu chí: Khối lượng phân tử (MW) < 500 Da; hệ số phân bố octanol/nước (LogP) < 5; số nhóm nhận liên kết hydro (HBA) < 10; số nhóm cho liên kết hydro (HBD) < 5. Hợp chất vi phạm không quá 1 tiêu chí được xem là có đặc tính giống thuốc dùng đường uống.
    • Dự đoán các chỉ số Dược động học và Độc tính (ADMET) qua công cụ pkCSM bằng chuỗi mã hóa SMILES: Hấp thu (tính thấm màng Caco-2, tỷ lệ hấp thu qua ruột người HIA); Phân bố (thể tích phân bố $V_{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); Chuyển hóa (khả năng ức chế enzym cytochrom P450 gồm CYP2D6 và CYP3A4); Thải trừ (thanh thải toàn phần và cơ chất của chất vận chuyển cation hữu cơ OCT2 ở thận); Độc tính (đột biến gen qua thử nghiệm AMES, ức chế kênh kali hERG gây độc tim, độc tính trên gan, kích ứng da, liều gây chết trung bình $LD_{50}$ và liều dung nạp tối đa).

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

Đặt vấn đề và Tổng quan tài liệu

Phần tổng quan hệ thống hóa toàn diện các kiến thức nền tảng về bệnh học ung thư vú, các phân loại mô bệnh học và phân tử, cùng các phương thức can thiệp lâm sàng tương ứng. Tác giả trình bày bảng phân loại giai đoạn bệnh và tỷ lệ sống sót sau 5 năm: giai đoạn tại chỗ (giai đoạn 0) đạt 100%; giai đoạn xâm lấn sớm (I, IIa, IIb) đạt 98,0% đối với khu trú tại chỗ và 83,6% tại vùng; giai đoạn tiến triển cục bộ (IIIa, IIIb, IIIc) đạt 57%; giai đoạn di căn xa (IV) giảm xuống còn 23,3%.

Về phân loại phân tử, ung thư vú được chia thành các nhóm: Luminal A (ER+, PR+, HER2-, tiên lượng tốt), Luminal B (ER+, PR+, HER2 âm tính hoặc dương tính, tiên lượng khá), nhóm HER2 dương tính (ER-, PR-, HER2+, tiên lượng xấu) và ung thư vú thể tam âm (ER-, PR-, HER2-, tiên lượng xấu). Phương pháp điều trị biến thiên theo từng giai đoạn, kết hợp giữa phẫu thuật (bảo tồn hoặc cắt bỏ triệt căn Patey), xạ trị hậu phẫu, hóa trị bổ trợ/tân bổ trợ, nội tiết trị liệu (tamoxifen, exemestane qua các thử nghiệm TEXT, SOFT) và các kháng thể đơn dòng điều trị đích như Trastuzumab, Pertuzumab.

Chương này cũng phân tích chi tiết cơ chế phân tử của thụ thể HER2. Khi HER2 tạo dimer dị phân (đặc biệt là cặp dị phân HER2/HER3), các dòng tín hiệu nội bào được khuếch đại mạnh mẽ nhất, thúc đẩy tế bào ác tính tăng sinh. Khung cấu trúc hóa học diphenylpropan ($C_6-C_3-C_6$) của isoflavone và sản phẩm chuyển hóa đường ruột S-equol được phân tích về khả năng cạnh tranh thụ thể estrogen và ức chế tyrosine kinase. Đồng thời, chương 1 giới thiệu lý thuyết về các thuật toán docking (Matching, Incremental Construction, Monte Carlo, Genetic Algorithm, Molecular Dynamics) và hàm tính điểm năng lượng tự do Gibbs ($\Delta G$).

Nguyên liệu, thiết bị, nội dung và phương pháp nghiên cứu

Chương 2 mô tả chi tiết quy trình chuẩn bị dữ liệu đầu vào và các bước thực nghiệm tính toán. Cấu trúc tinh thể 3PP0 được chuẩn hóa bằng cách loại bỏ các phân tử nước và phối tử nội sinh 03Q trên Discovery Studio Client, bổ sung hydro phân cực và tính toán điện tích bán thực nghiệm Kollman trên AutoDockTools.

Hộp lưới mô phỏng được giới hạn chính xác tại khoang gắn phối tử tự nhiên của HER2. Quy trình re-docking phối tử 03Q được thực hiện nhằm kiểm định tính chính xác của thuật toán AutoDock Vina. Danh sách 60 phối tử isoflavone được chuyển đổi từ tệp SDF sang PDB bằng Chimera, sau đó tối ưu hóa hình học năng lượng tối thiểu bằng Avogadro trước khi chuyển thành tệp .pdbqt. Các tiêu chuẩn đánh giá được lượng hóa bằng các ngưỡng toán học và dược lý cụ thể: tiêu chí Lipinski RO5, ngưỡng thấm Caco-2 (> 0,9), ngưỡng hấp thu đường ruột HIA (> 70% là tốt, < 30% là kém), ngưỡng phân bố thể tích $V_{Dss}$ (> 0,45 là cao, < -0,15 là thấp), ngưỡng an toàn thần kinh ($LogBB < -1$ và $LogPS < -3$), cùng các chỉ số an toàn độc tính.

Kết quả và bàn luận

Kết quả re-dock phối tử 03Q vào enzym 3PP0 đạt giá trị RMSD là 1,5 Å, năng lượng liên kết $\Delta G = -10,1$ kcal/mol. Phối tử re-dock tái lập được các tương tác với các acid amin cốt lõi ở khoang hoạt động tương tự như cấu trúc tinh thể nguyên bản, gồm liên kết $\pi$-alkyl với LYS753, VAL734, ALA751, MET774, LEU785, LEU796; liên kết $\pi$-$\sigma$ với LEU852, LEU726; liên kết hydro với GLN799, MET801; liên kết halogen với GLU770; tương tác với PHE1004 và Van der Waals với ASP863 (chỉ thiếu tương tác với ASN850 so với cấu trúc gốc). Kết quả này chứng minh tính hợp lệ của hệ thống mô phỏng.

Trong bước sàng lọc docking 60 hợp chất isoflavone, chất đối chứng Trastuzumab có mức năng lượng liên kết là -9,4 kcal/mol (thấp nhất trong 3 chất chứng), được dùng làm mốc so sánh. Kết quả xác định được 35 hợp chất có năng lượng liên kết $\Delta G \le -9,4$ kcal/mol.

Tiến hành đánh giá quy tắc 5 Lipinski cho 35 hợp chất này, tác giả đã sàng lọc được 20 hợp chất thỏa mãn điều kiện giống thuốc dùng đường uống (không vi phạm quá 1 tiêu chuẩn). Tiếp tục phân tích hồ sơ ADMET trên 20 hợp chất này, nghiên cứu ghi nhận:

  • Có 4 hợp chất đạt tính thấm tốt qua màng Caco-2 (giá trị xấp xỉ 0,897 - 0,9) và tỷ lệ hấp thu qua ruột người rất cao (93% - 96%), nổi bật là Genistein (93,387%).
  • Đánh giá phân bố cho thấy các chất Puerarin, 6''-O-Acetyldaidzin an toàn với hệ thần kinh do không qua hàng rào máu não ($LogBB < -1$) và không thấm vào hệ thần kinh trung ương ($LogPS < -3$).
  • Toàn bộ 20 hợp chất đều không ức chế CYP2D6 và không phải là cơ chất của chất vận chuyển OCT2 ở thận.
  • Đánh giá độc tính phát hiện nhiều hợp chất có nguy cơ gây độc: Isolupalbigenin và Kraussianone 3 có nguy cơ gây đột biến gen (AMES) và độc tim (ức chế hERG II); Alpinum isoflavone, Isoderrone gây độc tim và gan; Daidzin, Puerarin, Erypoegin K, Angustone A gây độc tim; Ononin gây độc gan. Có 8 hợp chất hoàn toàn không thể hiện độc tính.

Tổng hợp toàn diện các thông số, tác giả đã chọn lọc được 2 hợp chất ưu việt nhất là Genistein (PubChem ID: 5280961) và Biochanin A (PubChem ID: 5280373). Genistein có năng lượng liên kết -9,5 kcal/mol, tương tác với 7/16 acid amin vùng hoạt động (liên kết $\pi$-alkyl với LYS753, VAL734, ALA751, LEU796, MET774; liên kết $\pi$-$\sigma$ với LEU785; Van der Waals với ASP863, LEU852); hấp thu ở ruột đạt 93,387%; $V_{Dss}$ đạt 0,45; không độc tim, gan, da, không gây đột biến gen; $LD_{50}$ đạt 2,268 mol/kg và liều tối đa 0,478 log mg/kg/ngày (mặc dù có ức chế CYP3A4). Biochanin A đạt năng lượng liên kết -9,7 kcal/mol, tạo liên kết với 10/14 acid amin cốt lõi (LYS753, VAL734, ALA751, LEU852, LEU726, MET774, LEU796, LEU785, ASP863, GLU770); không ức chế CYP2D6, CYP3A4; không gây độc tế bào hay đột biến gen; $LD_{50}$ đạt 1,476 mol/kg.

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

Khóa luận đã hoàn thành quy trình sàng lọc ảo đa tầng từ thư viện phân tử tự nhiên đến các cấu trúc tiềm năng có khả năng ức chế enzym HER2. Bảng dưới đây tổng hợp các dữ liệu định lượng quan trọng nhất của 2 hợp chất tiềm năng so với phối tử đồng kết tinh và chất chứng Trastuzumab:

Tiêu chí phân tích Phối tử 03Q (Đồng kết tinh) Trastuzumab (Chất chứng) Genistein (ID: 5280961) Biochanin A (ID: 5280373)
Năng lượng liên kết ($\Delta G$) -10,1 kcal/mol -9,4 kcal/mol -9,5 kcal/mol -9,7 kcal/mol
Giá trị RMSD 1,5 Å - - -
Tương tác acid amin cốt lõi 13 acid amin (LYS753, VAL734, ALA751, MET774, LEU785...) Tương tác nhóm LYS753, VAL734, ALA751, GLN799... 7 acid amin cốt lõi (LYS753, VAL734, ALA751, LEU796, MET774...) 10 acid amin cốt lõi (LYS753, VAL734, ALA751, LEU852, LEU726...)
Tuân thủ quy tắc 5 Lipinski Không đánh giá Không áp dụng Đạt chuẩn (0 vi phạm) Đạt chuẩn (0 vi phạm)
Hấp thu ở ruột người (HIA) - - 93,387% ~93% - 96%
Tính thấm Caco-2 - - 0,897 (xấp xỉ 0,9) Khá (gần bằng Genistein)
Thể tích phân bố ($V_{Dss}$) - - 0,45 log L/kg (phân bố cao) Phân bố kém (< -0,15)
Chuyển hóa qua CYP - - Ức chế CYP3A4, không ức chế CYP2D6 Không ức chế CYP2D6, không ức chế CYP3A4
Cơ chất thải trừ OCT2 - - Không Không
Độc tính AMES & hERG - - Không gây độc Không gây độc
Độc tính gan & kích ứng da - - Không Không
Độc tính cấp $LD_{50}$ (chuột) - - 2,268 mol/kg (an toàn > 2) 1,476 mol/kg

Nghiên cứu đóng góp một quy trình sàng lọc in silico chuẩn hóa, kết hợp giữa mô phỏng liên kết không gian (docking bằng AutoDock Vina) với các bộ lọc hóa lý (Lipinski) và sinh dược học (pkCSM). Kết quả cung cấp bằng chứng điện toán chứng minh Genistein và Biochanin A có ái lực liên kết mạnh hơn chất đối chứng Trastuzumab trên miền xúc tác của HER2, đồng thời xác lập vị trí gắn kết và các liên kết phân tử chi tiết với các acid amin quan trọng của enzym.

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

Khóa luận đã thẳng thắn chỉ ra các hạn chế cố hữu của nghiên cứu:

  • Kết quả nghiên cứu hoàn toàn dựa trên các mô hình mô phỏng tính toán in silico; tồn tại sai số nhất định do việc sử dụng nhiều phễu lọc và phần mềm thuật toán học máy khác nhau.
  • Có khoảng cách giữa cấu trúc mô hình tĩnh trên máy tính và điều kiện sinh lý động học thực tế phức tạp trong cơ thể người.
  • Hợp chất Biochanin A vẫn còn một số điểm hạn chế về mặt dược động học như thể tích phân bố thấp ($V_{Dss}$), tính an toàn đối với hệ thần kinh chưa thật sự rõ ràng và chỉ số $LD_{50}$ chưa đạt ngưỡng tối ưu (> 2 mol/kg). Genistein có đặc tính ức chế enzym gan CYP3A4, tiềm ẩn nguy cơ tương tác thuốc khi sử dụng phối hợp.

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

  1. Triển khai các thử nghiệm thực nghiệm in vitro và in vivo đối với Genistein và Biochanin A nhằm xác định chính xác hoạt tính ức chế enzym HER2, đánh giá mức độ an toàn tế bào ở các dải nồng độ khác nhau và nghiên cứu cải thiện độ tan, tính thấm của hợp chất.
  2. Tiếp tục ứng dụng mô hình quy trình sàng lọc in silico này để tìm kiếm các hoạt chất mới tác động lên các đích phân tử khác trong điều trị ung thư vú hoặc mở rộng sang các nhóm hợp chất tự nhiên khác ngoài isoflavone.

Giá trị tham khảo

Tài liệu khóa luận có giá trị tham khảo hữu ích cho:

  • Sinh viên và học viên ngành Dược học, Hóa dược, Công nghệ sinh học: Cung cấp tài liệu mẫu mực về phương pháp thiết kế thuốc nhờ máy tính (CADD), đặc biệt là kỹ thuật docking phân tử và phương pháp dự đoán ADMET in silico.
  • Người nghiên cứu về ung thư và các hợp chất tự nhiên: Cung cấp dữ liệu chi tiết về cơ chế bệnh học của thụ thể HER2 trong ung thư vú, cấu trúc phân tử của nhóm isoflavone, và danh mục các tương tác phối tử - protein tại khoang xúc tác 3PP0.
  • Phần đáng tham khảo nhất:
    • Quy trình re-docking phối tử nội sinh để kiểm định tính hợp lệ của mô phỏng (tính toán RMSD và đối chiếu tương tác acid amin 2D/3D).
    • Khung phương pháp luận kết hợp phễu lọc 3 tầng: Năng lượng docking $\rightarrow$ Quy tắc Lipinski $\rightarrow$ Dự đoán ADMET toàn diện (Hấp thu, Phân bố, Chuyển hóa qua CYP, Đào thải qua OCT2 và các chỉ số độc tính hERG, AMES, $LD_{50}$).

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

1. Tiêu chí nào được sử dụng để chứng minh quy trình docking trong nghiên cứu là đáng tin cậy?

Quy trình docking được kiểm định thông qua việc re-dock phối tử đồng kết tinh 03Q vào cấu trúc tinh thể protein 3PP0. Tiêu chuẩn đánh giá là giá trị độ lệch bình phương trung bình gốc (RMSD) giữa phối tử sau khi re-dock so với phối tử gốc phải $\le 1,5$ Å. Kết quả thực nghiệm tính toán đạt RMSD chính xác là 1,5 Å và phối tử tái lập được các liên kết với hầu hết các acid amin thiết yếu trong khoang hoạt động (LYS753, VAL734, ALA751, MET774, LEU785, LEU796, LEU852, LEU726, GLN799, MET801, GLU770, PHE1004, ASP863), khẳng định quy trình docking đạt độ tin cậy khoa học.

2. Chất đối chứng nào được sử dụng làm ngưỡng năng lượng để sàng lọc 35 hợp chất tiềm năng?

Tác giả đã đưa vào mô phỏng 3 thuốc ức chế HER2 được FDA phê duyệt gồm Trastuzumab, Lapatinib và Tucatinib. Trong đó, Trastuzumab có năng lượng tự do liên kết âm nhất ($\Delta G = -9,4$ kcal/mol). Do đó, mức năng lượng $\Delta G \le -9,4$ kcal/mol của Trastuzumab được lựa chọn làm mốc tiêu chuẩn để sàng lọc các hợp chất isoflavone; kết quả có 35 trên tổng số 60 chất đạt tiêu chuẩn này.

3. Vì sao các hợp chất có năng lượng docking tốt như Isolupalbigenin hay Kraussianone 3 lại bị loại khỏi danh sách tối ưu?

Mặc dù có ái lực liên kết cao với protein đích, các hợp chất này bị loại ở bước sàng lọc an toàn sinh học thông qua công cụ pkCSM. Kết quả phân tích chỉ ra Isolupalbigenin và Kraussianone 3 đều có nguy cơ gây đột biến gen (dương tính với thử nghiệm AMES) và gây độc tính trên tim thông qua ức chế kênh kali hERG II; riêng Isolupalbigenin còn thể hiện độc tính trên gan, không đảm bảo yêu cầu an toàn để phát triển thành dược phẩm.

4. Hai hợp chất tiềm năng nhất được lựa chọn là gì và các chỉ số năng lượng liên kết cụ thể của chúng ra sao?

Hai hợp chất ưu việt nhất được chọn lọc là Genistein (PubChem ID: 5280961) và Biochanin A (PubChem ID: 5280373). Năng lượng tự do liên kết $\Delta G$ của Genistein đạt -9,5 kcal/mol (tương tác với 7 acid amin cốt lõi) và của Biochanin A đạt -9,7 kcal/mol (tương tác với 10 acid amin cốt lõi). Cả hai hợp chất đều có năng lượng liên kết tốt hơn chất đối chứng Trastuzumab (-9,4 kcal/mol).

5. Dữ liệu ADMET của Genistein có những ưu điểm và hạn chế nào được ghi nhận trong nghiên cứu?

  • Ưu điểm: Genistein có khả năng hòa tan tốt trong nước, tỷ lệ hấp thu qua ruột non người (HIA) rất cao đạt 93,387%, tính thấm Caco-2 đạt 0,897 (tiệm cận mức chuẩn 0,9), thể tích phân bố tối ưu ($V_{Dss} = 0,45$), an toàn đối với hệ thần kinh trung ương (khó qua hàng rào máu não), không gây độc tim (không ức chế hERG), không độc gan, không gây đột biến AMES và có mức liều gây chết trung bình an toàn ($LD_{50} = 2,268$ mol/kg).
  • Hạn chế: Hợp chất này có khả năng ức chế enzym cytochrom gan CYP3A4, cần được theo dõi nguy cơ tương tác chuyển hóa thuốc trong các nghiên cứu lâm sàng.

Kết luận

Khóa luận tốt nghiệp của sinh viên Đinh Xuân Tôn Tài đã hoàn thành xuất sắc việc xây dựng và thẩm định quy trình sàng lọc ảo in silico nhằm đánh giá tác dụng ức chế enzym HER2 của 60 hợp chất isoflavone tự nhiên định hướng điều trị ung thư vú. Thông qua phễu lọc khoa học kết hợp giữa mô phỏng docking phân tử trên AutoDock Vina, quy tắc 5 Lipinski và hệ thống dự đoán ADMET bằng pkCSM, nghiên cứu đã loại bỏ các chất kém hấp thu hoặc mang độc tính tim, gan, đột biến gen. Kết quả đã xác định được hai hoạt chất sáng giá nhất là Genistein và Biochanin A với ái lực liên kết mạnh hơn Trastuzumab và sở hữu các đặc tính dược động học đường uống khả quan. Đề tài cung cấp cơ sở dữ liệu lý thuyết và cấu trúc phân tử vững chắc, mở ra hướng nghiên cứu thực nghiệm in vitro và in vivo tiếp theo trong phát triển thuốc điều trị đích ung thư vú từ nguồn dược liệu tự nhiên.