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

Thụ thể yếu tố tăng trưởng biểu bì (Epidermal Growth Factor Receptor - EGFR), một thụ thể tyrosine kinase thuộc họ ErbB, đóng vai trò trung tâm trong việc điều hòa các quá trình sinh học như tăng sinh, biệt hóa và sống sót của tế bào. Sự biểu hiện quá mức hoặc đột biến bất thường của EGFR kích hoạt liên tục các con đường truyền tín hiệu hạ lưu (Ras/MAPK, PI3K/AKT, JAK/STAT), dẫn đến sự hình thành và tiến triển của nhiều khối u ác tính, đặc biệt là ung thư phổi không tế bào nhỏ (NSCLC), ung thư vú và ung thư đại trực tràng.

Các thuốc ức chế tyrosine kinase (TKIs) nhắm vào vị trí liên kết ATP trong miền kinase nội bào của EGFR đã mang lại bước tiến lớn trong điều trị ung thư. Tuy nhiên, sự xuất hiện của các đột biến kháng thuốc tại vị trí gắn ATP—tiêu biểu như đột biến "người gác cổng" T790M ở thế hệ thứ nhất (Gefitinib, Erlotinib) hay đột biến C797S đối với các thuốc thế hệ thứ ba (Osimertinib)—đã làm suy giảm nghiêm trọng hiệu quả điều trị. Do đó, việc tìm kiếm các chiến lược điều trị mới trở nên cấp thiết. Tiếp cận vị trí điều hòa dị lập thể (allosteric site) trên EGFR là một hướng đi triển vọng, cho phép điều chỉnh cấu hình không gian và ức chế hoạt tính enzym mà không cạnh tranh trực tiếp với ATP, đồng thời tạo cơ hội hiệp đồng tác dụng khi kết hợp với các TKI cạnh tranh ATP hiện có.

Về nguồn hoạt chất tự nhiên, nhóm triterpenoid khung lanostane mang vòng lacton 6 hoặc 7 cạnh (được gọi là colossolactone) phân lập từ loài nấm Ganoderma colossum đã được các thử nghiệm in vitro chứng minh có khả năng tăng cường hiệu lực của các thuốc nhắm trúng đích như gefitinib và sorafenib. Mặc dù tiềm năng hiệp đồng đã được ghi nhận, cơ chế tương tác ở cấp độ nguyên tử giữa các dẫn xuất colossolactone với EGFR cũng như bản chất hiệp đồng với gefitinib vẫn chưa được làm sáng tỏ. Khóa luận tốt nghiệp Dược sĩ với đề tài "Unraveling the synergistic mechanisms of colossolactone and gefitinib in cancer therapy via molecular dynamics simulations" được thực hiện nhằm giải quyết khoảng trống tri thức này thông qua phương pháp mô phỏng tính toán động lực học phân tử.

Khóa luận xác định 3 mục tiêu nghiên cứu cụ thể:

  1. Đánh giá tương tác liên kết và độ ổn định cấu trúc của 16 dẫn xuất colossolactone trên EGFR thông qua tích hợp docking phân tử và mô phỏng động lực học phân tử (MD).
  2. Xác định ái lực liên kết của các colossolactone với EGFR bằng phương pháp mô phỏng động lực học phân tử có định hướng (SMD), lấy mẫu dù (US) và tính toán thế năng lực trung bình (PMF).
  3. Phân tích vai trò của dẫn xuất colossolactone tiềm năng trong việc tăng cường ái lực liên kết của gefitinib với EGFR.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng: Cấu trúc không gian 3 chiều của miền tyrosine kinase EGFR ở hai trạng thái hoạt hóa (active state - mã PDB: 2GS6) và bất hoạt (inactive state - mã PDB: 2GS7); phân tử gefitinib; và 16 dẫn xuất colossolactone gồm nhóm số La Mã (I đến VIII) và nhóm chữ cái (A đến H).
  • Phạm vi nghiên cứu: Nghiên cứu tính toán in silico, tập trung vào ba vị trí gắn trên EGFR (túi ATP ở trạng thái hoạt hóa, túi ATP ở trạng thái bất hoạt và túi allosteric ở trạng thái bất hoạt), phân tích tương tác tĩnh, động học phân ly và mô phỏng động lực học phân tử dài hạn (lên tới 40 μs tổng thời gian mô phỏng) để chứng minh cơ chế điều hòa dị lập thể. Đề tài được thực hiện tại Bộ môn Hóa Dược, Khoa Công nghệ Hóa Dược, Trường Đại học Dược Hà Nội, hoàn thành vào tháng 6 năm 2025.

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

Khung lý thuyết và mô hình sinh học phân tử

Nghiên cứu dựa trên cấu trúc phân tử của miền tyrosine kinase EGFR (khoảng 540 acid amin), gồm hai thùy N-lobe và C-lobe liên kết qua vùng bản lề (hinge region). Túi liên kết ATP nằm giữa hai thùy với acid amin Thr790 đóng vai trò gác cổng. Túi allosteric kế cận chỉ xuất hiện khi vòng xoắn $\alpha\text{-C}$ dịch chuyển ra phía ngoài, đưa EGFR về trạng thái bất hoạt, liên quan mật thiết đến cầu muối Lys745–Glu762 và motif DFG (Asp855, Phe856, Gly857).

Các hợp chất colossolactone là nhóm triterpenoid có khung oxa-A-homo-steroidal hoặc oxa-A,B-dihomo-steroidal, mang nhiều nhóm chức phân cực và khung steroid kỵ nước linh hoạt, tạo điều kiện thuận lợi cho việc hình thành các liên kết hydro và tương tác Van der Waals với các hốc liên kết trên thụ thể.

Phương pháp nghiên cứu và công cụ tính toán

Khóa luận thiết lập một quy trình tính toán đa bước từ mô phỏng tĩnh đến mô phỏng cân bằng và không cân bằng:

  1. Chuẩn bị cấu trúc và hóa học lượng tử:

    • Cấu trúc protein EGFR (2GS6, 2GS7) được xử lý loại bỏ nước tinh thể và bổ sung các đoạn acid amin khuyết thiếu bằng máy chủ MODELLER.
    • Cấu trúc 2D của gefitinib và 16 dẫn xuất colossolactone (thu thập từ cơ sở dữ liệu PubChem) được tối ưu hóa hình học và tính toán phân bố điện tích nguyên tử RESP (Restrained Electrostatic Potential) bằng lý thuyết phiếm hàm mật độ (DFT) ở mức lý thuyết B3LYP/6-31+G(d,p) trên phần mềm Gaussian 16. Tham số topology phối tử được trích xuất bằng mô-đun antechamber của AmberTools23. Quỹ đạo phân tử biên (HOMO/LUMO) cũng được tính toán chi tiết.
  2. Docking phân tử (Molecular Docking):

    • Thực hiện trên phần mềm AutoDock Vina phiên bản 1.0 thông qua giao diện AutoDockTools (ADT) thuộc bộ MGLTools.
    • Hộp lưới tính toán (grid box) có kích thước $80 \times 80 \times 80\text{ \AA}$ bao phủ 3 vị trí chức năng: túi ATP trạng thái hoạt hóa, túi ATP trạng thái bất hoạt và túi allosteric trạng thái bất hoạt với độ thấu triệt (exhaustiveness) thiết lập ở mức 200 và bước lưới $0{,}375\text{ \AA}$.
  3. Mô phỏng Động lực học phân tử cổ điển (MD):

    • Thực hiện trên phần mềm GROMACS phiên bản 2023 sử dụng trường lực Amber99SB-ILDN và mô hình nước TIP3P.
    • Hệ thống được hòa tan trong hộp chữ nhật với khoảng cách tối thiểu đến biên $1{,}0\text{ nm}$, chứa hơn 96.000 nguyên tử, trung hòa điện tích bằng ion ngược.
    • Sau khi tối ưu hóa năng lượng (đến khi lực cực đại $< 100\text{ kJ/mol/nm}$), hệ được cân bằng NVT (100 ps) và NPT (100 ps), sau đó chạy mô phỏng sản xuất (production run) 100 ns ở nhiệt độ 300 K và áp suất 1 atm.
  4. Mô phỏng Động lực học phân tử có định hướng (SMD):

    • Áp dụng lực kéo cơ học không cân bằng lên khối tâm (COM) của phối tử với tốc độ kéo không đổi $1\text{ m/s}$ và hằng số lò xo $600\text{ kJ/mol/nm}^2$.
    • Hướng kéo được căn chỉnh dọc trục z của hộp mô phỏng dựa trên phân tích đường hầm thoát phối tử từ phần mềm CAVER. Nguyên tử Cα của protein được cố định bằng thế điều hòa ($100\text{ kJ/mol/nm}^2$).
    • Với mỗi phức hợp, thực hiện 50 quỹ đạo SMD độc lập (mỗi quỹ đạo 3 ns) để ghi nhận lực bứt phá cực đại ($F_m$), công kéo ($W_p$) và rào cản năng lượng tự do phân ly ($\Delta G_u^\ddagger$) theo đẳng thức Jarzynski cải biên bởi Hummer và Szabo.
  5. Lấy mẫu dù (Umbrella Sampling - US) và tính toán PMF:

    • Dùng để xác định năng lượng tự do liên kết trạng thái cân bằng. Hơn 120 cửa sổ (windows) được trích xuất từ quỹ đạo SMD kéo chậm ($0{,}1\text{ m/s}$), mỗi cửa sổ cách nhau $0{,}1\text{ nm}$.
    • Mỗi cửa sổ được cân bằng và mô phỏng MD 10 ns. Đường cong thế năng lực trung bình (PMF) được tái cấu trúc bằng phương pháp phân tích lược đồ trọng số (WHAM).
  6. Mô phỏng cơ chế hiệp đồng quy mô lớn:

    • Xây dựng 2 hệ mô phỏng: hệ đơn phối tử (chỉ có Gefitinib tại túi ATP) và hệ kép (Gefitinib tại túi ATP kết hợp Colossolactone H tại túi allosteric).
    • Thực hiện 10 quỹ đạo MD độc lập kéo dài $2\text{ }\mu\text{s}$ cho mỗi hệ, tạo ra tổng thời lượng mô phỏng lên tới $40\text{ }\mu\text{s}$ để đánh giá định lượng độ ổn định liên kết thông qua năng lượng tương tác (IE), khoảng cách không gian và số lượng tiếp xúc nguyên tử với 10 acid amin cốt lõi của túi ATP.

Nguồn dữ liệu và trang thiết bị

  • Cơ sở dữ liệu: RCSB PDB (mã 2GS6, 2GS7) và PubChem (CID của 16 colossolactone và gefitinib).
  • Hệ thống tính toán: 10 máy trạm hiệu năng cao, mỗi máy trang bị Dual CPU Intel Xeon Platinum 8180 (tổng 56 nhân, 112 luồng), 256GB RAM DDR4-3200 ECC, card đồ họa NVIDIA GeForce RTX 4080 SUPER 16GB và ổ cứng Samsung 990 PRO NVMe SSD.

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

Introduction (Tổng quan tài liệu và cơ sở nghiên cứu)

Chương mở đầu hệ thống hóa các kiến thức nền tảng về cấu trúc phân tử và con đường dẫn truyền tín hiệu của thụ thể EGFR trong sinh học ung thư. Tác giả phân tích chi tiết bốn thế hệ thuốc ức chế EGFR-TKI:

  • Thế hệ 1 (Gefitinib, Erlotinib, Lapatinib, Icotinib): Ức chế thuận nghịch cạnh tranh ATP.
  • Thế hệ 2: Gắn cộng hóa trị không thuận nghịch với Cys797.
  • Thế hệ 3 (Osimertinib, Rociletinib, Avitinib...): Nhắm chọn lọc đột biến T790M.
  • Thế hệ 4 (EAI-045, BLU-945, TQB-3804): Các chất ức chế dị lập thể không cạnh tranh ATP đang trong giai đoạn thử nghiệm.

Chương này cũng tổng hợp nguồn gốc hóa thực vật của chi nấm Ganoderma (đặc biệt là loài Ganoderma colossum), cấu trúc hóa học của 16 dẫn xuất colossolactone và các hoạt tính sinh học đã công bố (kháng HIV-1 protease, kháng khuẩn, kháng sốt rét, gây độc tế bào ung thư và hiệp đồng với 5-FU, gemcitabine, sorafenib, gefitinib). Phần cuối chương tổng kết các nguyên lý tính toán in silico gồm docking phân tử, mô phỏng MD cổ điển, SMD, US/PMF cùng những ưu thế và hạn chế của từng phương pháp.

Chapter 2: Materials and Methods (Nguyên liệu và phương pháp nghiên cứu)

Chương 2 mô tả toàn bộ thông số kỹ thuật, quy trình chuẩn bị dữ liệu và phương pháp tính toán:

  • Quy trình tối ưu hóa hóa học lượng tử DFT/B3LYP và gán điện tích RESP cho 16 phối tử colossolactone và gefitinib.
  • Tọa độ tâm và kích thước hộp lưới docking cho 3 vị trí trên EGFR (ATP hoạt hóa, ATP bất hoạt và allosteric bất hoạt).
  • Các bước thiết lập mô phỏng động lực học phân tử trên GROMACS 2023 với trường lực Amber99SB-ILDN.
  • Thiết lập thông số lực kéo SMD (vận tốc $1\text{ m/s}$, hằng số đàn hồi $600\text{ kJ/mol/nm}^2$, 50 quỹ đạo/hệ) và phương pháp xác định rào cản phân ly Hummer–Szabo.
  • Thiết lập lấy mẫu dù US (>120 cửa sổ, bước $0{,}1\text{ nm}$, thời gian $10\text{ ns/cửa sổ}$) và tích hợp dữ liệu bằng thuật toán WHAM.
  • Thiết kế hệ mô phỏng đánh giá tính hiệp đồng giữa Gefitinib và Colossolactone H với 10 quỹ đạo lặp lại, mỗi quỹ đạo $2\text{ }\mu\text{s}$ (tổng thời gian mô phỏng $40\text{ }\mu\text{s}$).

Chapter 3: Results and Discussion (Kết quả và bàn luận)

1. Đánh giá tương tác và độ ổn định của 16 dẫn xuất colossolactone trên EGFR

  • Docking phân tử: Cả 16 dẫn xuất colossolactone đều cho điểm năng lượng docking âm tại cả ba vị trí gắn.
    • Tại túi ATP hoạt hóa, colossolactone A tạo nhiều liên kết hydro nhất (5 liên kết với Leu718, Glu762, Phe795, Thr854, Asp855), trong khi colossolactone I, III, VII, F có năng lượng docking thấp hơn $-9{,}0\text{ kcal/mol}$.
    • Tại túi ATP bất hoạt, các hợp chất I, II, III, D, E, G đạt điểm năng lượng dưới $-10\text{ kcal/mol}$.
    • Tại túi allosteric bất hoạt, năng lượng docking dao động từ $-6{,}7$ đến $-10{,}3\text{ kcal/mol}$, với Phe856, Lys745 và Thr790 là các gốc cho liên kết hydro chính. Kết quả cho thấy số lượng liên kết hydro không tỷ lệ thuận tuyến tính với năng lượng docking.
  • Mô phỏng MD 100 ns: Toàn bộ 16 dẫn xuất duy trì trạng thái liên kết ổn định trong các túi gắn. Năng lượng tương tác Coulombic và Van der Waals đều mang giá trị âm, trong đó tương tác Van der Waals chiếm vai trò chủ đạo (hệ số tương quan Pearson giữa tổng số tương tác và năng lượng VDW đạt 0,94). Khoảng cách trung bình giữa phối tử và các acid amin cốt lõi (Gly721, Lys745, Thr790, Gln791, Leu792, Met793, Cys797, Arg841) đa số duy trì dưới $0{,}5\text{ nm}$.
  • Cấu hình đại diện (Representative Configurations): Phân tích cảnh quan năng lượng tự do Gibbs (trích xuất từ 50 ns cuối) cho thấy các liên kết hydro tại túi ATP biến động tương đối lớn trong quá trình mô phỏng, trong khi liên kết tại túi allosteric thể hiện tính ổn định cao và bền vững hơn.

2. Đánh giá ái lực liên kết qua mô phỏng SMD và US/PMF

  • Mô phỏng SMD 50 quỹ đạo:
    • Tại túi ATP hoạt hóa: Colossolactone F và VIII có lực bứt phá lớn nhất ($F_m = 1190\text{ pN}$ và $1130\text{ pN}$), công kéo $W_p$ đạt $120\text{ kcal/mol}$ và $119\text{ kcal/mol}$. Ngược lại, colossolactone H và G có lực bứt phá yếu nhất ($F_m = 509\text{ pN}$ và $596\text{ pN}$; $W_p = 45\text{ kcal/mol}$ và $47\text{ kcal/mol}$). Ái lực yếu tại túi ATP hoạt hóa là một đặc tính thuận lợi, tránh cạnh tranh trực tiếp với các thuốc ức chế ATP.
    • Tại túi ATP bất hoạt và túi allosteric: Colossolactone H thể hiện ưu thế vượt trội hoàn toàn. Tại túi allosteric, colossolactone H đạt lực bứt phá cao nhất ($F_m = 876\text{ pN}$), công kéo lớn nhất ($W_p = 87\text{ kcal/mol}$) và rào cản phân ly cao nhất ($\Delta G_u^\ddagger = 61\text{ kcal/mol}$).
Vị trí gắn trên EGFR Thông số SMD của Colossolactone H Giá trị ghi nhận
Túi ATP (trạng thái hoạt hóa) Lực bứt phá cực đại ($F_m$)
Công kéo ($W_p$)
Rào cản phân ly ($\Delta G_u^\ddagger$)
$509\text{ pN}$
$45\text{ kcal/mol}$
$28\text{ kcal/mol}$
Túi ATP (trạng thái bất hoạt) Lực bứt phá cực đại ($F_m$)
Công kéo ($W_p$)
Rào cản phân ly ($\Delta G_u^\ddagger$)
$892\text{ pN}$
$93\text{ kcal/mol}$
$68\text{ kcal/mol}$
Túi Allosteric (trạng thái bất hoạt) Lực bứt phá cực đại ($F_m$)
Công kéo ($W_p$)
Rào cản phân ly ($\Delta G_u^\ddagger$)
$876\text{ pN}$
$87\text{ kcal/mol}$
$61\text{ kcal/mol}$
  • Tính toán năng lượng tự do liên kết trạng thái cân bằng bằng US/PMF:
    • Năng lượng tự do liên kết ($\Delta G$) của Colossolactone H tại túi ATP hoạt hóa là $-5\text{ kcal/mol}$.
    • Tại túi ATP bất hoạt, $\Delta G = -13\text{ kcal/mol}$.
    • Tại túi allosteric bất hoạt, $\Delta G = -16\text{ kcal/mol}$ (cao hơn gấp 3 lần so với túi ATP hoạt hóa).
    • Kết quả khẳng định Colossolactone H có tính chọn lọc và ái lực liên kết nhiệt động học cao nhất tại vị trí allosteric ở trạng thái bất hoạt của EGFR.

3. Đặc điểm tương tác của Colossolactone H tại túi allosteric

  • Tương tác giữa Colossolactone H và miền kinase EGFR tại vị trí allosteric chủ yếu do lực Van der Waals chi phối (đóng góp $-52\text{ kcal/mol}$ so với Coulombic chỉ $-6\text{ kcal/mol}$).
  • Phân tích tương tác với 40 acid amin xung quanh cho thấy 12 acid amin có tần suất tiếp xúc cao nhất ($> 100$ tương tác/snapshot), dẫn đầu là Leu788 (311 tiếp xúc), Ile759 (304 tiếp xúc), Met766 (259 tiếp xúc), Leu747 (264 tiếp xúc) và Leu777 (231 tiếp xúc).
  • Nhóm 8 acid amin Leucine đóng góp 982 tiếp xúc (chiếm 30% tổng số 3270 tiếp xúc của toàn bộ phân tử) với tổng năng lượng VDW lên tới $-14\text{ kcal/mol}$.
  • Ba acid amin tạo tương tác VDW mạnh nhất là Ile759 ($-4{,}8\text{ kcal/mol}$), Met766 ($-4{,}5\text{ kcal/mol}$) và Asp855 ($-4{,}4\text{ kcal/mol}$).
  • Phân tích theo 3 phân đoạn cấu trúc của Colossolactone H (G1, G2, G3) ghi nhận năng lượng VDW lần lượt là $-14$, $-19$ và $-19\text{ kcal/mol}$, cùng một số tương tác tĩnh điện cục bộ đáng kể (G1–Ile759, G3–Met766, G3–Leu777, G2–Asp855).

4. Vai trò điều hòa dị lập thể và cơ chế hiệp đồng với Gefitinib

Trong mô phỏng đối chứng $2\text{ }\mu\text{s}$ lặp lại 10 quỹ đạo độc lập ($40\text{ }\mu\text{s}$ tổng cộng):

  • Hệ Gefitinib đơn độc (không có Colossolactone H): Phân tử gefitinib có xu hướng mất ổn định, phân ly khỏi túi gắn ATP, khuếch tán ra dung môi và tái gắn không đặc hiệu lên bề mặt protein. Chỉ có 1/10 quỹ đạo (quỹ đạo 6) duy trì được liên kết chức năng ổn định tại túi ATP (thỏa mãn đồng thời tiêu chí: năng lượng tương tác $\text{IE} < -46\text{ kcal/mol}$, khoảng cách $< 0{,}5\text{ nm}$, số tiếp xúc $> 188$).
  • Hệ phối hợp Gefitinib và Colossolactone H: Có tới 9/10 quỹ đạo duy trì ổn định đồng thời cả ba tiêu chí liên kết chức năng của gefitinib tại túi ATP.
Hệ mô phỏng (10 quỹ đạo $\times 2\text{ }\mu\text{s}$) Số quỹ đạo duy trì liên kết chức năng ổn định Năng lượng tương tác tổng với 10 acid amin túi ATP ($\sum \text{IE}$) Tương tác đặc hiệu tại Arg841 ($\text{IE}_{\text{Arg841}}$)
Gefitinib đơn độc (G) 1 / 10 quỹ đạo (chỉ quỹ đạo 6 đạt yêu cầu) Dao động từ $-0{,}1$ đến $-43{,}0\text{ kcal/mol}$ (9/10 quỹ đạo bị mất tương tác) Yếu hoặc bằng 0 ở hầu hết các quỹ đạo (riêng traj6 đạt $-43{,}0\text{ kcal/mol}$)
Gefitinib + Colossolactone H (G&C) 9 / 10 quỹ đạo (trừ traj9) Đạt từ $-46{,}1$ đến $-108{,}0\text{ kcal/mol}$ ở 9 quỹ đạo Rất mạnh và ổn định: $-49{,}1$ đến $-112{,}2\text{ kcal/mol}$ ở 9 quỹ đạo
  • Vai trò then chốt của gốc Arg841: Trong hệ phối hợp, Arg841 duy trì năng lượng tương tác âm rất lớn với gefitinib. Về mặt sinh học, Arg841 điều hòa quá trình tự phosphoryl hóa và hoạt hóa qua tín hiệu oxy hóa khử (sulfenyl hóa Cys797). Sự có mặt của Colossolactone H tại túi allosteric tạo hiệu ứng cản trở không gian hợp lý, định hướng nguyên tử nitơ N1 trên vòng quinazoline của gefitinib tạo 2 liên kết hydro bền vững với Arg841 ngay từ thời điểm $1\text{ ns}$, khóa chặt gefitinib trong túi gắn ATP và ngăn chặn sự phân ly.

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

Kết quả chính

  • Xác định toàn diện cấu hình liên kết và độ ổn định của 16 dẫn xuất colossolactone trên miền tyrosine kinase của EGFR ở cả trạng thái hoạt hóa và bất hoạt.
  • Chứng minh bằng nhiệt động học và động học phân ly rằng Colossolactone H là dẫn xuất tối ưu nhất để ức chế dị lập thể, với ái lực liên kết cân bằng tại túi allosteric đạt $\Delta G = -16\text{ kcal/mol}$, lực bứt phá $F_m = 876\text{ pN}$ và rào cản phân ly $\Delta G_u^\ddagger = 61\text{ kcal/mol}$.
  • Giải mã mạng lưới tương tác kỵ nước dày đặc giữa Colossolactone H và 40 acid amin tại túi allosteric, trong đó các gốc Leucine (chiếm 30% số tiếp xúc) cùng với Ile759, Met766 và Asp855 đóng vai trò mỏ neo cấu trúc.
  • Cung cấp bằng chứng định lượng qua $40\text{ }\mu\text{s}$ mô phỏng MD chứng minh Colossolactone H nâng tỷ lệ duy trì liên kết chức năng của Gefitinib từ 10% (1/10 quỹ đạo) lên 90% (9/10 quỹ đạo) thông qua cơ chế cố định không gian và tăng cường tương tác hydro với gốc Arg841.

Giải pháp và kiến nghị tác giả đề xuất

  • Tiến hành các mô phỏng MD mở rộng trên 16 dẫn xuất colossolactone đối với các thể đột biến kháng thuốc phổ biến của EGFR như T790M và C797S để đánh giá khả năng khắc phục tình trạng kháng TKI.
  • Ứng dụng quy trình mô phỏng in silico tích hợp này để nghiên cứu cặp phối hợp giữa Colossolactone G và Sorafenib nhằm tìm hiểu cơ chế hiệp đồng tương tự trên ung thư biểu mô tế bào gan.
  • Dựa trên khung cấu trúc và phân bố tương tác của Colossolactone H, tiến hành thiết kế hợp lý (rational drug design) các dẫn xuất mới nhằm tối ưu hóa ái lực và độ chọn lọc cho túi allosteric, phát triển các chất điều hòa dị lập thể thế hệ mới.

Đóng góp khoa học

Đây là công trình nghiên cứu tính toán in silico đầu tiên khảo sát toàn diện dòng hợp chất colossolactone trên thụ thể EGFR. Khóa luận đã làm sáng tỏ ở mức độ nguyên tử cơ chế tác động hiệp đồng đã quan sát được trên thực nghiệm giữa triterpenoid tự nhiên từ nấm Ganoderma colossum và thuốc trúng đích gefitinib, mở ra hướng tiếp cận phối hợp thuốc giữa chất điều hòa dị lập thể tự nhiên và TKI cạnh tranh ATP trong điều trị ung thư.


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

Hạn chế của nghiên cứu

  • Do giới hạn về chi phí tài nguyên tính toán, việc tính toán năng lượng tự do liên kết trạng thái cân bằng bằng kỹ thuật lấy mẫu dù (Umbrella Sampling với hơn 120 cửa sổ) chỉ mới thực hiện tập trung trên dẫn xuất tiêu biểu nhất là Colossolactone H, chưa thể áp dụng trên toàn bộ 15 dẫn xuất còn lại.
  • Các mô phỏng MD cổ điển sử dụng trường lực Amber99SB-ILDN với các xấp xỉ cơ học cổ điển, chưa tính đến các hiệu ứng phân cực động, chuyển dịch điện tích lượng tử hay khả năng tạo liên kết cộng hóa trị trong môi trường sinh học phức tạp.
  • Khảo sát cơ chế hiệp đồng mới dừng lại ở cấu trúc EGFR dạng hoang dã (wild-type), chưa thực hiện mô phỏng tương tác trên các cấu trúc EGFR mang đột biến kháng thuốc kép hoặc kháng thuốc bậc cao (T790M/C797S).

Hướng nghiên cứu tiếp theo

  • Mở rộng mô phỏng tính toán trên các dạng đột biến EGFR kháng TKI thế hệ mới (T790M, C797S, L858R).
  • Mở rộng khảo sát cơ chế điều hòa trên hệ phức hợp Colossolactone G – Sorafenib.
  • Thiết kế và tối ưu hóa cấu trúc hóa học các dẫn xuất bán tổng hợp từ Colossolactone H để nâng cao độ gắn kết dị lập thể trước khi tiến hành thử nghiệm sinh học in vitro và in vivo.

Giá trị tham khảo

Tài liệu khóa luận có giá trị tham khảo chuyên môn cao đối với:

  • Sinh viên, học viên cao học và nghiên cứu sinh ngành Dược học, Hóa dược và Hóa tin sinh: Tham khảo một quy trình chuẩn mực, hoàn chỉnh và hiện đại trong thiết kế thuốc hỗ trợ bởi máy tính (CADD/SBDD), tích hợp đồng bộ từ Hóa tính toán lượng tử (DFT/RESP), Docking phân tử, Mô phỏng động lực học phân tử cổ điển (MD 100 ns), Mô phỏng kéo không cân bằng (SMD), Lấy mẫu dù (Umbrella Sampling/PMF) đến Mô phỏng dài hạn đánh giá tính hiệp đồng ($40\text{ }\mu\text{s}$).
  • Các nhà nghiên cứu phát triển thuốc ung thư: Cung cấp dữ liệu cấu trúc chi tiết về túi gắn allosteric của EGFR và cơ chế phân tử của việc phối hợp chất điều hòa allosteric nhằm tăng cường hiệu lực của các thuốc ức chế kinase cạnh tranh ATP.
  • Phần đáng tham khảo nhất: Chương 2 (phương pháp luận thiết lập thông số mô phỏng và công thức tính rào cản phân ly Hummer–Szabo), Mục 3.2 (so sánh năng lượng tự do giữa túi ATP và allosteric qua US/PMF) và Mục 3.4 (phương pháp phân tích 10 quỹ đạo MD ensemble $2\text{ }\mu\text{s}$ để lượng hóa cơ chế hiệp đồng qua acid amin Arg841).

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

1. Tại sao Colossolactone H lại được lựa chọn để phân tích sâu bằng Umbrella Sampling thay vì các dẫn xuất khác?

Trong số 16 dẫn xuất được sàng lọc qua mô phỏng SMD 50 quỹ đạo, Colossolactone H thể hiện các thông số vượt trội nhất tại túi allosteric ở trạng thái bất hoạt của EGFR, bao gồm lực bứt phá cực đại cao nhất ($F_m = 876\text{ pN}$), công kéo lớn nhất ($W_p = 87\text{ kcal/mol}$) và rào cản phân ly cao nhất ($\Delta G_u^\ddagger = 61\text{ kcal/mol}$). Đồng thời, Colossolactone H có ái lực yếu tại túi ATP hoạt hóa ($F_m = 509\text{ pN}$), giúp hạn chế sự cạnh tranh bất lợi với các thuốc ức chế ATP.

2. Năng lượng tự do liên kết ($\Delta G$) của Colossolactone H tại các vị trí trên EGFR chênh lệch như thế nào qua tính toán US/PMF?

Kết quả tái cấu trúc đường cong thế năng lực trung bình (PMF) từ hơn 120 cửa sổ lấy mẫu dù cho thấy:

  • Túi ATP ở trạng thái hoạt hóa: $\Delta G = -5\text{ kcal/mol}$.
  • Túi ATP ở trạng thái bất hoạt: $\Delta G = -13\text{ kcal/mol}$.
  • Túi Allosteric ở trạng thái bất hoạt: $\Delta G = -16\text{ kcal/mol}$. Mức năng lượng $-16\text{ kcal/mol}$ tại túi allosteric cao hơn gấp hơn 3 lần so với túi ATP hoạt hóa, khẳng định tính chọn lọc nhiệt động học vượt trội của Colossolactone H đối với vị trí allosteric khi EGFR ở trạng thái bất hoạt.

3. Tương tác hóa học nào đóng vai trò cốt lõi giúp Colossolactone H gắn ổn định vào túi allosteric của EGFR?

Tương tác Van der Waals (VDW) đóng vai trò áp đảo với mức năng lượng $-52\text{ kcal/mol}$, so với đóng góp tương tác tĩnh điện Coulombic chỉ $-6\text{ kcal/mol}$. Khung kỵ nước của Colossolactone H liên kết chặt chẽ với mạng lưới các gốc kỵ nước, đặc biệt là 8 acid amin Leucine (chiếm 30% tổng số tiếp xúc với năng lượng VDW đạt $-14\text{ kcal/mol}$), cùng với ba acid amin đóng góp lực VDW mạnh nhất là Ile759 ($-4{,}8\text{ kcal/mol}$), Met766 ($-4{,}5\text{ kcal/mol}$) và Asp855 ($-4{,}4\text{ kcal/mol}$).

4. Cơ chế phân tử cụ thể nào chứng minh tính hiệp đồng giữa Colossolactone H và Gefitinib trong mô phỏng $40\text{ }\mu\text{s}$?

Sự có mặt của Colossolactone H tại túi allosteric kế cận tạo ra sự chiếm chỗ không gian hợp lý, ngăn không cho gefitinib đi quá sâu hoặc bị mất định hướng. Nhờ đó, nguyên tử nitơ N1 trên vòng quinazoline của gefitinib được cố định để tạo 2 liên kết hydro bền vững với gốc Arg841 ngay từ $1\text{ ns}$. Tương tác này duy trì năng lượng liên kết rất mạnh ($-49{,}1$ đến $-112{,}2\text{ kcal/mol}$) tại Arg841 trên 9/10 quỹ đạo, giữ gefitinib không bị phân ly khỏi túi ATP và ức chế cấu hình xúc tác phosphoryl hóa của EGFR.


Kết luận

Khóa luận tốt nghiệp Dược sĩ của sinh viên Thái Chính Tâm đã ứng dụng thành công quy trình mô phỏng tính toán in silico chuyên sâu để làm sáng tỏ cơ chế tương tác và tính chất hiệp đồng của nhóm hợp chất tự nhiên colossolactone trên thụ thể EGFR. Nghiên cứu đã xác định Colossolactone H là chất điều hòa dị lập thể tiềm năng nhất với ái lực liên kết nhiệt động học vượt trội tại túi allosteric ($\Delta G = -16\text{ kcal/mol}$) thông qua mạng lưới tương tác Van der Waals bền vững. Đặc biệt, phân tích động lực học phân tử quy mô $40\text{ }\mu\text{s}$ đã chứng minh Colossolactone H nâng cao rõ rệt độ ổn định liên kết của gefitinib tại túi ATP (từ 10% lên 90% số quỹ đạo) nhờ cơ chế cố định tương tác với acid amin then chốt Arg841. Công trình cung cấp cơ sở khoa học mức độ nguyên tử vững chắc, mở ra định hướng phát triển các phác đồ phối hợp thuốc nhắm trúng đích và thiết kế các chất ức chế EGFR allosteric mới trong điều trị ung thư.