Tổng quan về luận án

Nghiên cứu của tác giả Trần Quế Hương thuộc chuyên ngành Hóa dược (Mã số: 9720203) tại Đại học Y Dược TP. Hồ Chí Minh với đề tài "Nghiên cứu sàng lọc các cấu trúc phân tử nhỏ có khả năng gắn kết IL-6, IL-6R và đánh giá khả năng gắn kết bằng các phương pháp lý sinh" đại diện cho một bước tiến tiên phong trong việc phát triển chất ức chế tương tác protein-protein (Small-Molecule Protein-Protein Interaction Inhibitors - SMPPIIs). Trong bối cảnh y sinh học hiện đại, interleukin-6 (IL-6) đóng vai trò trung tâm trong quá trình chuyển tiếp từ viêm cấp tính sang viêm mạn tính, là nguyên nhân then chốt dẫn đến hội chứng bão cytokine trong COVID-19, đồng thời kích hoạt bất thường trục truyền tín hiệu IL-6/JAK/STAT3 thúc đẩy tăng sinh, xâm lấn và di căn tế bào ung thư. Mặc dù các kháng thể đơn dòng sinh học như Tocilizumab, Sarilumab và Olokizumab đã được ứng dụng lâm sàng, chúng tồn tại nhiều nhược điểm nghiêm trọng: chi phí sản xuất đắt đỏ, bắt buộc sử dụng đường tiêm, thời gian bán hủy kéo dài và nguy cơ sinh miễn dịch cao.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: "Trên đường truyền tín hiệu của IL6, các thế hệ thuốc cấu trúc phân tử nhỏ ức chế IL6 thông qua ức chế các janus kinase (JAKs) như tofacitinib, ruxolitinib, baricitinib hay ức chế gp130 như bazedoxifen, raloxifen đã được phát triển, nhưng vẫn chưa có chất phân tử nhỏ nào ức chế thông qua gắn kết trực tiếp lên IL6 hoặc IL6R được công bố." Các nghiên cứu quốc tế trước đây của Li (2014), Zhou (2014) và Sato (2022) chủ yếu dừng lại ở mô hình docking tĩnh, thiếu phân tích động lực học phân tử (Molecular Dynamics - MD) chuyên sâu, hoặc bỏ qua việc đánh giá toàn diện các thông số dược động học ADMET đường uống.

Luận án thiết lập 3 câu hỏi nghiên cứu và 3 giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Những đặc trưng không gian và điểm nóng tương tác (hot-spots) nào tại vị trí I (IL-6/IL-6R) và vị trí IIIa (IL-6/gp130) quyết định ái lực gắn kết của các phối tử phân tử nhỏ?
  • Câu hỏi nghiên cứu 2 (RQ2): Liệu các phân tử nhỏ được sàng lọc từ ngân hàng dữ liệu thông qua tích hợp 3D-pharmacophore, docking phân tử và mô phỏng MD có thể duy trì độ bền liên kết và đạt năng lượng tự do gắn kết âm ($\Delta G_{bind} < 0$) hay không?
  • Câu hỏi nghiên cứu 3 (RQ3): Các phối tử tiềm năng in silico có khả năng ức chế trực tiếp sự tạo phức IL-6/IL-6R trên thử nghiệm lý sinh in vitro và ngăn chặn sự tăng sinh tế bào ung thư phụ thuộc IL-6 ngoại sinh hay không?

Tương ứng là các giả thuyết khoa học:

  • Giả thuyết 1 (H1): Việc nhắm trúng các acid amin trọng yếu tại vị trí I (Phe74, Phe78, Leu174, Leu178 trên IL-6; Phe229, Tyr230, Phe279 trên IL-6R) và vị trí IIIa (Leu57, Trp157, Asp160 trên IL-6) sẽ phá vỡ liên kết không gian lập thể giữa cytokin và thụ thể.
  • Giả thuyết 2 (H2): Phối tử thỏa mãn quy tắc 5 Lipinski và các chỉ số ADMET sàng lọc qua MD 100 ns sẽ thể hiện năng lượng tự do gắn kết MM-PBSA vượt trội so với các phối tử không chọn lọc.
  • Giả thuyết 3 (H3): Các chất ức chế chọn lọc sẽ ức chế tỷ lệ gắn kết IL-6/IL-6R trên kit ELISA cạnh tranh và làm giảm tỷ lệ sống sót của các dòng tế bào ung thư đại tràng HT-29, SW480 và ung thư gan HepG2, HuH-7 dưới sự hoạt hóa của IL-6 ngoại sinh.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Lý thuyết Nhận diện Phân tử (Molecular Recognition Theory), Thuyết Điểm nóng Tương tác Protein-Protein (PPI Hot-Spots Theory) và Động lực học Tín hiệu Cytokin (Classic vs Trans-signaling Framework). Phạm vi nghiên cứu bao quát 3 cơ sở dữ liệu lớn (DrugBank, ZINC15, NCI), mô phỏng động lực học phân tử 100 ns cho hàng chục phức hợp, đánh giá trên kit thử nghiệm sinh hóa lý sinh RayBio® IL-6/IL-6R và 5 dòng tế bào ung thư người.

Literature Review và Positioning

Nghiên cứu về IL-6 bắt đầu từ năm 1986 khi Hirano và cộng sự phân lập BSF-2 (B-cell stimulatory factor 2). Đến năm 1988, Yamasaki cùng nhóm nghiên cứu đã giải mã cấu trúc thụ thể mIL-6R, và năm 1990 Kishimoto phát hiện đồng thụ thể truyền tín hiệu chung gp130. Các công trình kinh điển của Heinrich (2003) và Rose-John (2006, 2018) đã phân định ranh giới giữa con đường tín hiệu chính thống (classic-signaling qua mIL-6R) chịu trách nhiệm cân bằng nội môi, tạo máu, tái tạo gan, và con đường tín hiệu chuyển đổi (trans-signaling qua sIL-6R) kích hoạt viêm mạn tính và phát triển khối u.

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc về cơ chế "tín hiệu cụm" (cluster-signaling) và vai trò của sgp130. Heink và cộng sự (2017) cho rằng phức hợp IL-6/sIL-6R trên tế bào chuyển chỉ liên kết gp130 trên tế bào nhận và sgp130 chỉ đối kháng chọn lọc tín hiệu chuyển đổi mà không ức chế tín hiệu chính thống. Ngược lại, Lamertz và cộng sự (2018) báo cáo bằng chứng thực nghiệm cho thấy sgp130 có thể liên kết trực tiếp với IL-6/mIL-6R trên màng tế bào, vô hiệu hóa cả trục tín hiệu chính thống. Sự mâu thuẫn này đặt ra yêu cầu phải tìm kiếm các phân tử nhỏ ức chế chọn lọc tại vị trí I (IL-6/IL-6R) nhằm can thiệp trực tiếp từ ngọn nguồn tương tác phối tử - thụ thể sơ cấp.

Luận án định vị chính xác khoảng trống nghiên cứu khi so sánh đối chiếu với các công trình quốc tế tiêu biểu:

  1. So sánh với nghiên cứu của Li và cộng sự (2014): Li sử dụng AutoDock 4 xác định Bazedoxifene (-8,2 kcal/mol; $IC_{50} = 24,9\ \mu M$) và Raloxifene (-8,8 kcal/mol; $IC_{50} = 43,2\ \mu M$) gắn kết gp130 tại vị trí IIIa. Tuy nhiên, nghiên cứu của Li chỉ dừng lại ở docking tĩnh, không khảo sát tính biến thiên không gian theo thời gian bằng MD và không tối ưu hóa các thông số dược động học dùng đường uống.
  2. So sánh với nghiên cứu của Sato và cộng sự (2022): Sato công bố dẫn chất A196 (F2432-0054) tương tác với IL-6 tại vị trí I (Phe74, Phe78, Cys73, Ser176) với ái lực $K_D = 1,5\ \mu M$, ức chế tăng sinh tế bào 7TD1 ở nồng độ 5–10 $\mu M$. Dẫu vậy, Sato chưa thiết lập quy trình tính toán năng lượng tự do gắn kết toàn phần MM-PBSA và chưa đánh giá cơ chế chéo trên nhiều dòng tế bào biểu mô ung thư người.
  3. So sánh với nghiên cứu của Sharma và cộng sự (2022): Sharma đánh giá Genistein và Quercetin gắn kết IL-6/IL-6R qua AutoDock và MD 50 ns, ghi nhận $IC_{50}$ của Quercetin đạt $38\ \mu M$ sau 24 giờ. Tuy nhiên, nghiên cứu này không tích hợp mô hình 3D-pharmacophore dựa trên PPIs, không sử dụng bộ lọc ADMET toàn diện và thiếu các thử nghiệm lý sinh cạnh tranh trực tiếp xác nhận cơ chế ức chế tương tác protein.

Luận án của NCS. Trần Quế Hương đã khắc phục triệt để các hạn chế trên bằng cách thiết lập một chu trình nghiên cứu khép kín từ sàng lọc đa tầng in silico đến thực nghiệm lý sinh và sinh học tế bào.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng Lý thuyết Điểm nóng Tương tác Protein-Protein (Hot-Spot PPI Theory do Clackson & Wells đề xướng) sang hệ thống cytokine - thụ thể phức tạp. Nghiên cứu chứng minh rằng mặc dù bề mặt tiếp xúc IL-6/IL-6R rất rộng và phẳng, phần lớn năng lượng gắn kết tự do lại tập trung vào các cụm acid amin cục bộ:

  • Tại vị trí I trên IL-6: Cụm kỵ nước - vòng thơm Phe74, Phe78 (thuộc vòng nối AB loop), Leu174, Leu178 kết hợp với các gốc tích điện Arg179, Arg182 (chuỗi xoắn D).
  • Tại vị trí I trên IL-6R: Phe229, Tyr230, Phe279 lọt vào các túi kỵ nước tương ứng và tương tác cầu muối/liên kết hydro của Glu277, Glu278, Arg250, Asp272.
  • Tại vị trí IIIa trên IL-6: Leu57, Trp157 và Asp160 tương tác trực tiếp với miền D1 của gp130 (Ile83, Phe36, Tyr94, Asn92).
[Mô hình Khái niệm Tương tác & Ức chế Trục IL-6]

+-------------------------------------------------------------------------+
|                              CYTOKINE IL-6                              |
|   Vị trí I: Phe74, Phe78, Leu174, Leu178, Arg179, Arg182                |
|   Vị trí IIIa: Leu57, Trp157, Asp160                                    |
+------------------------------------+------------------------------------+
                                     |
              (Phong tỏa bằng Phân tử nhỏ SMPPIIs / Flavonoid)
                                     |
                                     v
+------------------------------------+------------------------------------+
|                         THỤ THỂ ĐÍCH                                    |
|   IL-6R (Vị trí I): Phe229, Tyr230, Phe279, Glu277, Glu278              |
|   gp130 (Vị trí IIIa): Ile83, Phe36, Tyr94, Asn92                       |
+------------------------------------+------------------------------------+
                                     |
                      (Ngăn chặn hình thành phức hợp)
                                     |
                                     v
+-------------------------------------------------------------------------+
|                  ỨC CHẾ ĐƯỜNG TRUYỀN TÍN HIỆU HẠ LƯU                    |
|   Giảm phosphoryl hóa JAK -> Bất hoạt STAT3 / MAPK / PI3K-Akt           |
|   Ức chế biểu hiện gen viêm (CRP, SAA, hepcidin) & Tăng sinh ung thư    |
+-------------------------------------------------------------------------+

Nghiên cứu cụ thể hóa 4 mệnh đề lý thuyết (Theoretical Propositions):

  • Mệnh đề 1 (P1): Sự phong tỏa đồng thời túi kỵ nước tạo bởi Phe279 (IL-6R) và liên kết ion tại Glu277/Glu278 bởi một phân tử nhỏ sẽ làm suy giảm ái lực gắn kết IL-6/IL-6R trên 50%.
  • Mệnh đề 2 (P2): Các phân tử nhỏ có cấu trúc polyphenol/flavone phẳng sở hữu các nhóm hydroxyl phenolic phân bố đối xứng là cấu trúc tối ưu để bắt chước tương tác $\pi-\pi$ stacking của Phe74 và Phe78.
  • Mệnh đề 3 (P3): Động lực học liên kết tại vị trí IIIa phụ thuộc vào sự chiếm giữ không gian của vòng indole trên Trp157; phối tử chèn vào vị trí này sẽ ngăn chặn sự tiếp cận của miền D1 gp130.
  • Mệnh đề 4 (P4): Thuốc phân tử nhỏ đạt sinh khả dụng đường uống cao đối với đích PPI cytokine phải tuân thủ nghiêm ngặt diện tích bề mặt phân cực $TPSA < 140\ \text{Å}^2$ và số liên kết xoay $\le 10$.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn 5 giai đoạn:

  1. Thiết kế Pharmacophore 3D dựa trên PPIs: Xây dựng các mô hình truy vấn đặc tính hóa học (Ph-IL6-1a, Ph-IL6-1b, Ph-IL6R, Ph-IL6-3a1, Ph-IL6-3a2) từ cấu trúc phức hợp tinh thể X-ray với kháng thể 6F82 và Olokizumab.
  2. Docking phân tử và phân tích dấu vân tay tương tác (PLIF): Thực hiện docking chéo (cross-docking) trên các mô hình D-IL6-1a, D-IL6-1b, D-IL6R, D-IL6-3a; sàng lọc các phối tử có điểm số docking cao và tỷ lệ bắt cặp PLIF trên 70% với các điểm nóng.
  3. Bộ lọc ADMET đa tham số in silico: Dự đoán độ hấp thu qua ruột người (HIA), tính thấm hàng rào máu não (BBB), ái lực ức chế CYP450 và điểm số khả năng tổng hợp (SAscore < 5).
  4. Mô phỏng Động lực học phân tử (MD 100 ns) & Tính năng lượng MM-PBSA: Khảo sát độ biến thiên cấu hình (RMSD), độ linh động acid amin (RMSF), bán kính hồi chuyển (Rg), diện tích tiếp cận dung môi (SASA) và năng lượng tự do liên kết $\Delta G_{bind}$.
  5. Thực nghiệm Lý sinh & Tế bào học: Định lượng ái lực ức chế qua ELISA cạnh tranh (RayBio®) và kiểm chứng ức chế tăng sinh tế bào MTT trên HepG2, HuH-7, LNCaP, SW480, HT-29.

Điều kiện biên (boundary conditions): Khung phân tích áp dụng cho các cấu trúc phân tử nhỏ hữu cơ có khối lượng phân tử $\le 800\ \text{Da}$, nhắm đích ngoại bào tại vị trí tiếp xúc tương tác sơ cấp I và IIIa.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức luận thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp luận định lượng đa tầng (Multi-level Quantitative Design). Thiết kế kết hợp chặt chẽ giữa tính toán lý thuyết in silico với kiểm chứng thực nghiệm lý sinh in vitro nhằm triệt tiêu các kết quả dương tính giả (false-positive hits).

[Quy trình Nghiên cứu Đa tầng Thực chứng]

+-------------------------------------------------------------------------+
|                  NGÂN HÀNG CƠ SỞ DỮ LIỆU PHÂN TỬ VÀ PROTEIN              |
|   - RCSB PDB: IL-6, IL-6R, Phức hợp 6F82, Olokizumab                    |
|   - Thư viện: DrugBank (DB), ZINC15, NCI (> 100.000 hợp chất)           |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
|                  SÀNG LỌC ẢO IN SILICO ĐA THẦM                          |
|   - 3D-Pharmacophore (MOE) -> Molecular Docking -> Phân tích PLIF       |
|   - Sàng lọc ADMET (SwissADME, pkCSM): Lipinski, SAscore < 5            |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
|             MÔ PHỎNG ĐỘNG LỰC HỌC PHÂN TỬ (MD 100 ns)                   |
|   - GROMACS, trường lực CHARMM36/AMBER                                  |
|   - Đánh giá: RMSD < 2 Å, RMSF, Rg, SASA, H-bond persistence            |
|   - Năng lượng tự do gắn kết MM-PBSA / MM-GBSA                          |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
|                  KIỂM CHỨNG LÝ SINH & IN VITRO                          |
|   - Thử nghiệm gắn kết cạnh tranh: Kit ELISA RayBio® IL-6/IL-6R         |
|   - Đo mật độ quang OD 450 nm, tính toán tỷ lệ % ức chế và IC50         |
|   - Thử nghiệm MTT: Tăng sinh tế bào HT-29, SW480, HepG2, HuH-7, LNCaP  |
+-------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

  • Thu thập dữ liệu cấu trúc: Trích xuất tọa độ không gian 3D của phức hợp IL-6/IL-6R/gp130 từ Protein Data Bank (RCSB PDB), xử lý proton hóa, loại bỏ phân tử nước không tham gia liên kết cấu trúc và tối thiểu hóa năng lượng với phần mềm MOE (Molecular Operating Environment).
  • Sàng lọc Pharmacophore và Docking: Xây dựng mô hình pharmacophore gồm 4 đến 6 đặc tính không gian (vùng kỵ nước, nhân thơm, nhóm cho/nhận liên kết hydro). Thực hiện docking bán linh động (semi-flexible) trên AutoDock Vina và MOE.
  • Tiêu chuẩn chọn/loại trừ ADMET:
    • Tiêu chuẩn chọn: Khối lượng phân tử $\le 500\ \text{Da}$; $LogP \le 5$; số nhóm cho hydro $HBD \le 5$; số nhóm nhận hydro $HBA \le 10$; độ hấp thu qua ruột $HIA > 30%$; $SAscore < 5$.
    • Tiêu chuẩn loại trừ: Hợp chất vi phạm trên 1 quy tắc Lipinski, chứa nhóm cấu trúc PAINS (Pan-Assay Interference Compounds), hoặc có độc tính dự đoán qua pkCSM cao (đột biến Ames, ức chế kênh hERG).
  • Giao thức mô phỏng Động lực học phân tử: Thiết lập hộp mô phỏng điều hòa với dung môi nước TIP3P, trung hòa điện tích bằng ion $Na^+/Cl^-$, cân bằng hệ thống ở điều kiện nhiệt độ không đổi $NVT$ ($300\ \text{K}$) và áp suất không đổi $NPT$ ($1\ \text{bar}$) trong $100\ \text{ps}$. Thời gian chạy quỹ đạo MD đạt $100\ \text{ns}$ ($100.000\ \text{ps}$) trên phần mềm GROMACS.
  • Giao thức thử nghiệm sinh hóa ELISA: Sử dụng kit RayBio® Human IL-6/IL-6R Binding Assay. Protein IL-6R được cố định trên giếng 96, ủ đồng thời với phối tử phân tử nhỏ và IL-6 tái tổ hợp ở $37^\circ\text{C}$ trong 4 giờ. Sử dụng kháng thể kháng IL-6 gắn HRP kết hợp cơ chất TMB, dừng phản ứng bằng dung dịch acid và đo mật độ quang quang phổ tại bước sóng $\lambda = 450\ \text{nm}$.
  • Giao thức MTT sinh học tế bào: Nuôi cấy tế bào ung thư gan (HepG2, HuH-7), ung thư tuyến tiền liệt (LNCaP) và ung thư đại tràng (SW480, HT-29) trong môi trường chuyên dụng. Xử lý tế bào bằng IL-6 ngoại sinh kết hợp hoặc không kết hợp với chất ức chế thử nghiệm ở các nồng độ pha loãng bậc 2 trong 24, 48, 72 và 96 giờ. Đo tinh thể formazan hòa tan tại $\lambda = 540\ \text{nm}$ bằng máy BioTek.

Data và phân tích

  • Mẫu dữ liệu: Hơn 100.000 cấu trúc phân tử nhỏ từ 3 cơ sở dữ liệu (DrugBank, ZINC15, NCI).
  • Phân tích động lực học: Tính toán độ lệch RMSD nguyên tử protein và phối tử qua $100\ \text{ns}$ (ngưỡng ổn định $< 2\ \text{Å}$), phân tích dao động RMSF tại từng vị trí acid amin, diện tích SASA và năng lượng tự do MM-PBSA với công cụ gmx_MMPBSA.
  • Phân tích thống kê sinh học: Giá trị phần trăm ức chế (%) và nồng độ ức chế $50%$ ($IC_{50}$) được tính toán từ các đường cong đáp ứng nồng độ (nồng độ thử nghiệm từ 0,1 $\mu M$ đến 100 $\mu M$). So sánh thống kê giữa các nhóm xử lý và nhóm đối chứng bằng phân tích phương sai ANOVA một chiều và phép kiểm định sau ANOVA Dunnett/Tukey, mức ý nghĩa được xác lập tại $^P<0,05$, $^{}P<0,01$, $^{}P<0,001$, $^{****}P<0,0001$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Giải mã bản đồ năng lượng và hình thái liên kết tại vị trí I và IIIa: Khẳng định Phe74, Phe78, Leu174, Leu178 (IL-6) và Phe229, Tyr230, Phe279 (IL-6R) là các điểm nóng không thể thay thế. Phân tích PLIF trên mô hình D-IL6R cho thấy các phối tử đạt ái lực cao nhất đều tạo liên kết hydro bền vững với Glu277 và Glu278 với tần suất duy trì $> 60%$ trong suốt $100\ \text{ns}$ mô phỏng.
  2. Xác định các phối tử phân tử nhỏ tiềm năng ức chế trực tiếp tương tác IL-6/IL-6R: Đã phát hiện và kiểm chứng thực nghiệm 6 hợp chất hàng đầu gồm Quercetin, Dabigatran, Dacomitinib, Imatinib, Entospletinib và Rhamnetin. Trong đó:
    • Thử nghiệm ELISA cạnh tranh sơ bộ ở nồng độ $100\ \mu M$ ghi nhận Quercetin, Dabigatran và Dacomitinib ức chế mạnh sự hình thành phức hợp IL-6/IL-6R với tỷ lệ ức chế tương ứng từ $45%$ đến $> 70%$.
    • Quercetin và Dacomitinib duy trì cấu hình siêu ổn định trong khoang gắn kết IL-6R với giá trị dao động RMSD trung bình $< 1,8\ \text{Å}$ và năng lượng tự do gắn kết MM-PBSA đạt mức âm sâu ($\Delta G_{bind} < -25\ \text{kcal/mol}$).
  3. Phát hiện tái định vị thuốc (Drug Repurposing) đột phá đối với nhóm thuốc "tinib": Dacomitinib (chất ức chế EGFR kinase) và Entospletinib (chất ức chế SYK kinase) thể hiện khả năng gắn kết kép độc đáo: vừa ức chế cạnh tranh trực tiếp tại vị trí I của IL-6/IL-6R, vừa có tiềm năng ngăn chặn tầng tín hiệu nội bào hạ lưu.
  4. Hiệu quả ức chế tăng sinh tế bào ung thư đại tràng qua trung gian IL-6:
    • Trên dòng tế bào HT-29 và SW480, IL-6 ngoại sinh thúc đẩy tăng sinh tế bào rõ rệt sau 24 và 48 giờ ($P < 0,001$).
    • Sự hiện diện của 6 chất ức chế tiềm năng (Imatinib, Dacomitinib, Entospletinib, Dabigatran, Quercetin, Rhamnetin) ở nồng độ $0,1 - 10\ \mu M$ đã đảo ngược hoàn toàn tác động kích thích tăng sinh của IL-6 ngoại sinh, làm giảm tỷ lệ tế bào sống sót một cách có ý nghĩa thống kê ($P < 0,0001$).
  5. Khám phá các hợp chất mới nhắm vị trí IIIa: Sàng lọc từ thư viện ZINC15 và NCI đã tìm ra các cấu trúc ZINC01183742 và ZINC13348421 có khả năng chèn sâu vào rãnh tiếp xúc giữa chuỗi D của IL-6 và gp130, tạo liên kết hydro ổn định với Trp157 và Asp160.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp minh chứng thực nghiệm xác thực mô hình ức chế dị lập thể và cạnh tranh trực tiếp PPIs của cytokine; mở rộng hiểu biết về việc phân tử nhỏ có thể vô hiệu hóa tín hiệu chuyển đổi (trans-signaling) mà không cần dùng kháng thể protein cồng kềnh.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình CADD tích hợp từ 3D-pharmacophore PPIs $\rightarrow$ Docking $\rightarrow$ ADMET $\rightarrow$ MD 100 ns $\rightarrow$ MM-PBSA $\rightarrow$ ELISA $\rightarrow$ MTT. Quy trình này có tính khả chuyển cao, dễ dàng áp dụng cho các cytokine khác thuộc họ 4 vòng xoắn như IL-11, LIF, OSM, CNTF.
  • Về mặt ứng dụng lâm sàng và thực tiễn: Khám phá cấu trúc khung polyphenol (Quercetin, Rhamnetin, Chrysin, Baicalin) và khả năng tái định vị Dacomitinib mở ra hướng bào chế các thuốc điều trị phối hợp dạng uống, giúp hạ giá thành điều trị cho bệnh nhân viêm khớp dạng thấp, bão cytokine và ung thư đại trực tràng kháng thuốc.
  • Khuyến nghị chính sách y tế: Thiết lập lộ trình phát triển tiền lâm sàng cho các dẫn chất phân tử nhỏ ức chế IL-6/IL-6R, giảm sự phụ thuộc vào các chế phẩm sinh học nhập khẩu đắt tiền.

Limitations và Future Research

  • Hạn chế 1: Luận án sử dụng cấu trúc tinh thể X-ray tĩnh làm điểm khởi đầu cho in silico; các biến đổi cấu hình lớn (large-scale conformational changes) của phức hợp thụ thể màng gp130 khi có mặt lipid kép chưa được mô phỏng toàn phần.
  • Hạn chế 2: Thử nghiệm in vitro mới dừng lại ở mức tế bào dòng (cell lines) 2D và kit sinh hóa; chưa đánh giá trên mô hình nuôi cấy 3D organoid hoặc vi môi trường khối u có sự tương tác miễn dịch phức tạp.
  • Hạn chế 3: Chưa tiến hành các thử nghiệm in vivo trên động vật biến đổi gen (mô hình chuột viêm khớp collagen-induced arthritis hoặc chuột mang khối u dị ghép xenograft) để kiểm tra dược động học và độc tính trường diễn thực tế.

Chương trình nghiên cứu tương lai:

  1. Tổng hợp hóa học các dẫn xuất tối ưu hóa cấu trúc (SAR optimization) từ khung phân tử Quercetin và ZINC01183742 để nâng cao độ hòa tan và ái lực gắn kết ($K_D$ ở mức nanomolar).
  2. Kiểm chứng cấu trúc phức hợp tinh thể phân tử nhỏ - protein bằng phương pháp nhiễu xạ tia X (X-ray crystallography) hoặc kính hiển vi điện tử đông lạnh (Cryo-EM).
  3. Đánh giá tác động chọn lọc trên con đường tín hiệu chính thống so với tín hiệu chuyển đổi thông qua đo lường mức độ phosphoryl hóa STAT3 (p-STAT3) và các protein pha cấp tính.
  4. Triển khai thử nghiệm tiền lâm sàng trên mô hình động vật viêm loét đại tràng (DSS-induced colitis) và ung thư biểu mô gan.

Tác động và ảnh hưởng

  • Tác động học thuật: Đóng góp một cơ sở dữ liệu hoàn chỉnh về các phân tử nhỏ gắn kết IL-6 và IL-6R tại vị trí I và IIIa; công bố các bài báo trên các tạp chí quốc tế chuyên ngành Hóa dược uy tín thuộc danh mục Scopus/SCIE; dự kiến tạo nguồn trích dẫn học thuật lớn cho các nghiên cứu về SMPPIIs.
  • Chuyển đổi ngành công nghiệp dược phẩm: Mở ra phân khúc nghiên cứu phát triển (R&D) thuốc phân tử nhỏ dùng đường uống kháng IL-6, thay thế hoặc bổ trợ cho các kháng thể đơn dòng giá hàng ngàn USD mỗi liều điều trị.
  • Lợi ích xã hội và người bệnh: Thuốc phân tử nhỏ có chi phí tổng hợp thấp, bảo quản ở điều kiện nhiệt độ phòng thông thường, tiện lợi khi sử dụng đường uống, giúp bệnh nhân mắc các bệnh viêm mạn tính và ung thư dễ dàng tuân thủ phác đồ điều trị lâu dài.
  • Tầm vóc quốc tế: Đưa năng lực nghiên cứu thiết kế thuốc in silico kết hợp lý sinh của Việt Nam tiếp cận ngang tầm với các nhóm nghiên cứu tiên tiến trên thế giới.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận quy trình phương pháp luận chuẩn mực về sàng lọc ảo dựa trên cấu trúc kết hợp mô phỏng MD chuyên sâu và phân tích MM-PBSA.
  • Các nhà khoa học Hóa dược & Dược lý: Sở hữu các cấu trúc khởi nguồn (hit compounds) đã được kiểm chứng hoạt tính để tiếp tục phát triển dẫn chất và thử nghiệm cơ chế phân tử.
  • Doanh nghiệp Dược phẩm R&D: Rút ngắn thời gian và tiết kiệm hàng triệu USD chi phí sàng lọc ngẫu nhiên nhờ áp dụng phễu lọc in silico chính xác cao.
  • Bác sĩ lâm sàng và Bệnh nhân: Có triển vọng tiếp cận các liệu pháp ức chế cytokine dạng uống hiệu quả, an toàn và chi phí hợp lý trong tương lai.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Luận án đã mở rộng Lý thuyết Điểm nóng Tương tác Protein-Protein (Hot-Spot PPI Theory) bằng cách xác định chính xác rằng việc ức chế tương tác IL-6/IL-6R tại vị trí I không đòi hỏi phải che phủ toàn bộ diện tích bề mặt tiếp xúc rộng lớn mà chỉ cần một phân tử nhỏ khóa chặt túi kỵ nước tại Phe279 kết hợp tạo liên kết hydro bền vững với cụm acid amin tích điện âm Glu277/Glu278 trên IL-6R, hoặc tạo tương tác $\pi-\pi$ xếp chồng với Phe74/Phe78 trên IL-6.

2. Tính đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây?

So với nghiên cứu của Li (2014) chỉ sử dụng docking tĩnh trên gp130, và nghiên cứu của Sato (2022) thiếu tính toán năng lượng liên kết nhiệt động lực học, luận án đã tiên phong tích hợp mô phỏng động lực học phân tử GROMACS $100\ \text{ns}$, giải mã biến thiên cấu hình thời gian thực (RMSD, RMSF, Rg, SASA), tính năng lượng tự do MM-PBSA, kết hợp màng lọc ADMET dự đoán đường uống và kiểm chứng bằng kit ELISA cạnh tranh RayBio® chuẩn hóa.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là các thuốc ức chế kinase phân tử nhỏ thế hệ mới như Dacomitinib và Entospletinib có khả năng gắn kết chéo (cross-binding) với ái lực rất cao vào vị trí tương tác I của IL-6/IL-6R, ức chế trực tiếp sự tạo phức protein ngoại bào với hiệu quả tương đương hoặc vượt trội hơn các phối tử tự nhiên, làm giảm mạnh sự tăng sinh tế bào ung thư HT-29 ở nồng độ cực thấp $0,1\ \mu M$.

4. Luận án có cung cấp giao thức tái lập (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp chi tiết toàn bộ mã PDB, cấu hình phần mềm MOE, tham số trường lực CHARMM36 trong GROMACS, các bước tối thiểu hóa năng lượng, điều kiện cân bằng $NVT/NPT$, nồng độ thử nghiệm ELISA cạnh tranh, bước sóng đo quang phổ ($\lambda = 450\ \text{nm}$ cho ELISA và $\lambda = 540\ \text{nm}$ cho MTT), cùng quy trình nuôi cấy 5 dòng tế bào ung thư đảm bảo tính tái lập 100% tại các phòng thí nghiệm tiêu chuẩn.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm gồm 3 giai đoạn: (1) Tối ưu hóa cấu trúc hóa dược tạo thư viện dẫn xuất có $K_D < 100\ \text{nM}$; (2) Xác định cấu trúc tinh thể tia X đồng kết tinh của phức hợp protein - phối tử và đánh giá dược động học in vivo trên chuột; (3) Thử nghiệm tiền lâm sàng pha I/II trên các mô hình bệnh viêm ruột mạn tính và ung thư đại tràng kháng hóa trị.

Kết luận

Luận án tiến sĩ của NCS. Trần Quế Hương mang lại 5 đóng góp học thuật cốt lõi:

  1. Xây dựng thành công chu trình sàng lọc ảo tích hợp kết hợp 3D-pharmacophore PPIs, docking phân tử, màng lọc ADMET và mô phỏng động lực học phân tử $100\ \text{ns}$ có độ chính xác cao.
  2. Lập bản đồ tương tác chi tiết các điểm nóng tại vị trí I (IL-6/IL-6R) và vị trí IIIa (IL-6/gp130), làm sáng tỏ cơ chế nhận diện phân tử của các phối tử phân tử nhỏ.
  3. Phát hiện và chứng minh thực nghiệm 6 phân tử nhỏ hàng đầu (đặc biệt là Quercetin, Dacomitinib, Dabigatran) có khả năng ức chế trực tiếp sự tạo phức IL-6/IL-6R trên kit ELISA cạnh tranh với tỷ lệ ức chế $> 45 - 70%$ ở nồng độ $100\ \mu M$.
  4. Kiểm chứng hoạt tính sinh học tế bào, chứng minh các chất ức chế tiềm năng ngăn chặn hiệu quả sự tăng sinh của các dòng tế bào ung thư HT-29, SW480, HepG2 được hoạt hóa bởi IL-6 ngoại sinh ($P < 0,0001$).
  5. Mở ra hướng nghiên cứu mới về tái định vị thuốc và thiết kế thuốc SMPPIIs đường uống nhắm vào cytokine, tạo tiền đề vững chắc cho việc phát triển các liệu pháp điều trị trúng đích bệnh lý viêm mạn tính và ung thư trong tương lai.