Tổng quan nghiên cứu

Ung thư vú và sarcoma cơ vân đang là những thách thức y tế nghiêm trọng trên phạm vi toàn cầu lẫn tại Việt Nam. Theo thống kê của tổ chức GLOBOCAN năm 2020, ung thư vú là loại ung thư phổ biến nhất ở nữ giới với tỷ lệ mắc mới chiếm 11,7% (tương đương khoảng 2,3 triệu ca mắc mới trên toàn thế giới). Tại Việt Nam, tỷ lệ ung thư vú chiếm 11,8%, đứng thứ ba trong số các bệnh ung thư phổ biến nhất sau ung thư gan (14,5%) và ung thư phổi (14,4%). Song song đó, sarcoma cơ vân là dạng ung thư mô mềm ác tính thường gặp nhất ở trẻ em dưới 15 tuổi; ước tính tại Hoa Kỳ có khoảng 400 đến 500 ca mắc mới mỗi năm (chiếm khoảng 3% ung thư nhi), trong khi tại Việt Nam căn bệnh này chiếm tới 50% các trường hợp sarcoma mô mềm ở trẻ em.

Mặc dù hợp chất tự nhiên curcumin chiết xuất từ củ nghệ vàng sở hữu nhiều hoạt tính sinh học quý giá, việc ứng dụng lâm sàng bị hạn chế đáng kể do cấu trúc β-diketon kém bền, độ tan thấp, tốc độ chuyển hóa nhanh và sinh khả dụng đường uống rất thấp. Nhằm giải quyết vấn đề này, đề tài nghiên cứu tập trung tìm kiếm các dẫn chất tương đồng curcumin kiểu monocarbonyl kết hợp dị vòng 1H-pyrazol nhằm gia tăng độ ổn định hóa học và nâng cao hiệu lực ức chế tế bào ác tính.

Mục tiêu cụ thể của luận văn bao gồm: tổng hợp hóa học 9 dẫn chất tương đồng curcumin kiểu monocarbonyl không đối xứng chứa dị vòng 1H-pyrazol; đánh giá hoạt tính độc tế bào in vitro trên dòng tế bào ung thư vú người MCF-7; sàng lọc in silico ái lực gắn kết của 111 cấu trúc dẫn chất trên thụ thể Estrogen Alpha và thụ thể IGF-1R; đồng thời thực hiện mô phỏng động lực học phân tử 100 ns nhằm xác định các phân tử tiềm năng có tương tác bền vững. Nghiên cứu được triển khai trong giai đoạn 2019-2021 tại Khoa Dược, Đại học Y Dược Thành phố Hồ Chí Minh, đóng góp cơ sở dữ liệu quan trọng để phát triển các hoạt chất điều trị ung thư thế hệ mới có hiệu lực cao gấp 10 đến 20 lần so với curcumin tự nhiên.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết của nghiên cứu dựa trên các cơ chế bệnh sinh phân tử và hóa dược cấu trúc:

Thụ thể Estrogen Alpha (ERα) là yếu tố điều hòa phiên mã quan trọng thuộc nhóm thụ thể nhân, được mã hóa bởi gen ESR1 gồm 595 acid amin với trọng lượng phân tử 66 kDa. Trong các khối u vú, ERα biểu hiện quá mức với tỷ lệ từ 50% đến 80%, so với mức dưới 10% ở biểu mô tuyến vú bình thường. Sự gắn kết của các chất đối kháng tại khoang liên kết phối tử làm biến đổi vị trí của vòng xoắn 12, đặc biệt là hình thành liên kết hydro then chốt với acid amin Arg394, từ đó vô hiệu hóa tín hiệu tăng sinh tế bào ung thư.

Thụ thể yếu tố tăng trưởng giống Insulin 1 (IGF-1R) là một thụ thể tyrosin kinase xuyên màng hiện diện ở hơn 80% trường hợp sarcoma cơ vân. Khi được kích hoạt, IGF-1R kích thích chuỗi truyền tín hiệu qua trục PI3K-AKT-mTOR, ức chế quá trình chết theo chương trình và kích hoạt tăng sinh khối u. Việc ức chế cạnh tranh tại túi liên kết ATP thông qua tương tác với các acid amin bản lề như Met1082 và Glu1080 đóng vai trò quyết định trong việc ngăn chặn đường truyền tín hiệu sinh khối u.

Về mặt hóa dược học, mô hình biến đổi cấu trúc curcuminoid chuyển từ khung β-diketon sang khung monocarbonyl giúp loại bỏ vị trí chuyển hóa nhanh, trong khi việc gắn dị vòng 1H-pyrazol giàu điện tử giúp tăng cường mật độ liên kết hydro và tương tác kỵ nước với khoang thụ thể đích.

+-----------------------------------------------------------------------------------+
|                        KHUNG LÝ THUYẾT & ĐÍCH TÁC ĐỘNG                            |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|   [ Cấu trúc Curcumin ] ---> [ Loại bỏ β-diketon ] ---> [ Gắn dị vòng 1H-Pyrazol ] |
|                                                                |                  |
|                                                                v                  |
|                         [ Dẫn chất Monocarbonyl Curcumin (MACs) ]                 |
|                                        |                                          |
|                    +-------------------+-------------------+                      |
|                    |                                       |                      |
|                    v                                       v                      |
|        [ Kháng Ung Thư Vú ]                       [ Kháng Sarcoma Cơ Vân ]        |
|        - Đích tác động: ERα                       - Đích tác động: IGF-1R         |
|        - Tương tác chìa khóa: Arg394              - Tương tác chìa khóa: Met1082  |
|        - Dòng tế bào: MCF-7                       - Cơ chế: Ức chế Tyrosin Kinase |
|                                                                                   |
+-----------------------------------------------------------------------------------+

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

Nghiên cứu áp dụng quy trình tích hợp giữa hóa học thực nghiệm, thử nghiệm sinh học in vitro và mô phỏng tính toán in silico:

Cỡ mẫu và thiết kế thực nghiệm: Quy trình tổng hợp thực hiện trên 9 dẫn chất không đối xứng thông qua phản ứng ngưng tụ aldol giữa chất trung gian A và chất trung gian B với tỷ lệ mol 1:1,5, sử dụng tiền chất methyl acetoacetat và 4-clorophenylhydrazin qua phản ứng Vilsmeier-Haack. Cấu trúc của 9 hợp chất tổng hợp được kiểm nghiệm và xác nhận độ tinh khiết thông qua các thông số nhiệt độ nóng chảy, sắc ký lớp mỏng (TLC), quang phổ hồng ngoại (IR), khối phổ (MS), cùng phổ cộng hưởng từ hạt nhân 1H-NMR và 13C-NMR.

Đánh giá độc tính tế bào in vitro: Thử nghiệm MTT được tiến hành trên 36 dẫn chất đại diện trong thư viện 111 chất để xác định nồng độ ức chế 50% sự tăng sinh tế bào (IC50) trên dòng tế bào ung thư vú MCF-7 ở các dải nồng độ khảo sát từ 1 µM đến 100 µM.

Gắn kết phân tử và mô phỏng động lực học: Toàn bộ thư viện 111 dẫn chất được gắn kết phân tử vào cấu trúc tinh thể của ERα (mã PDB: 3ERT) và IGF-1R (mã PDB: 3O23). Phương pháp redocking phối tử đồng kết tinh được thực hiện để kiểm chứng quy trình với độ lệch bình phương trung bình nhỏ hơn 2,0 Å. Tiếp theo, 7 phức hợp phối tử - protein tiềm năng nhất (gồm 3 phức hợp trên ERα và 4 phức hợp trên IGF-1R) được đưa vào mô phỏng động lực học phân tử trong thời gian 100 ns trên phần mềm GROMACS dưới điều kiện đẳng nhiệt - đẳng áp (nhiệt độ 300 K, áp suất 1 bar, mô hình dung môi nước TIP3P). Các thông số dao động RMSD, độ linh động acid amin RMSF, bán kính quay Rg, diện tích bề mặt tiếp xúc dung môi SASA và năng lượng gắn kết tự do theo mô hình MM/PBSA được tính toán chi tiết.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Tổng hợp hóa học và hiệu suất phản ứng: Nghiên cứu đã tổng hợp thành công 9 dẫn chất monocarbonyl không đối xứng chứa dị vòng 1H-pyrazol. Hiệu suất tổng hợp chất trung gian A dao động từ 35% đến 65%, trong khi hiệu suất tổng hợp các dẫn chất MACs cuối cùng đạt từ 25% đến 77%. Các dữ liệu phổ cộng hưởng từ hạt nhân và khối phổ xác thực chính xác cấu trúc hóa học dự kiến của toàn bộ 9 dẫn chất.

Hoạt tính gây độc tế bào in vitro: Kết quả sàng lọc 36 dẫn chất trên dòng tế bào ung thư vú MCF-7 ghi nhận 5 hợp chất tiềm năng thể hiện hoạt tính gây độc tế bào mạnh với giá trị IC50 dao động trong khoảng từ 2,55 µM đến 13,49 µM. Đáng chú ý, các dẫn chất này thể hiện hoạt tính ức chế tế bào ung thư cao hơn gấp nhiều lần so với curcumin tự nhiên khi thử nghiệm ở cùng điều kiện nồng độ.

Gắn kết phân tử in silico: Trên thụ thể ERα, 9 dẫn chất đạt điểm số docking âm rất sâu, dao động từ -11,1 kcal/mol đến -9,8 kcal/mol, hình thành liên kết hydro bền vững với acid amin Arg394 và tương tác kỵ nước với Glu353, Leu346, Ala350, Leu387, Met421. Trên thụ thể IGF-1R, 7 dẫn chất đạt năng lượng gắn kết từ -9,6 kcal/mol đến -9,0 kcal/mol, tạo liên kết hydro then chốt với Met1082 và mạng lưới tương tác kỵ nước với Glu1080, Met1156 và Leu1005.

Độ ổn định động lực học phân tử 100 ns: Hợp chất CUR.5 thể hiện độ ổn định vượt trội khi tạo phức với ERα, duy trì giá trị RMSD dao động dưới ngưỡng 2,0 Å và bảo toàn liên kết hydro với Arg394 trong suốt quá trình chạy 100 ns. Tương tự, hợp chất CUR.1 gắn kết chặt chẽ vào túi xúc tác của IGF-1R với độ biến thiên RMSD tối thiểu và năng lượng gắn kết tự do MM/PBSA đạt giá trị âm cao.

+---------------------------------------------------------------------------------------+
|                BẢNG TỔNG HỢP CÁC KẾT QUẢ THỰC NGHIỆM & IN SILICO CHÍNH                |
+---------------------+-------------------+---------------------+-----------------------+
| Hạng mục thông số   | Giá trị đạt được  | Mục tiêu đối chiếu  | Đánh giá ý nghĩa      |
+---------------------+-------------------+---------------------+-----------------------+
| Hiệu suất trung gian| 35% - 65%         | > 30%               | Phản ứng ổn định      |
| Hiệu suất MACs      | 25% - 77%         | > 20%               | Đạt chuẩn tổng hợp    |
| IC50 trên MCF-7     | 2,55 - 13,49 µM   | < 25 µM             | 5 dẫn chất vượt trội  |
| Điểm Docking ERα    | -11,1 đến -9,8    | < -8,0 kcal/mol     | 9 chất gắn kết mạnh   |
| Điểm Docking IGF-1R | -9,6 đến -9,0     | < -8,0 kcal/mol     | 7 chất ức chế cao     |
| Thời gian chạy MDs  | 100 ns            | Chuẩn mô phỏng      | CUR.5 & CUR.1 ổn định |
+---------------------+-------------------+---------------------+-----------------------+

Thảo luận kết quả

Các kết quả thực nghiệm và tính toán cấu trúc giải thích rõ mối quan hệ giữa cấu trúc và hoạt tính sinh học (SAR). Việc thay thế khung β-diketon bằng cầu nối monocarbonyl giúp phân tử duy trì cấu hình phẳng tối ưu, giảm sự cồng kềnh không gian và ngăn chặn sự phân hủy thủy phân trong môi trường sinh lý. Đồng thời, sự hiện diện của dị vòng 1H-pyrazol với các nhóm thế thơm mang halogen làm tăng mật độ điện tử, tạo điều kiện hình thành các liên kết π-alkyl và tương tác phân tán kỵ nước mạnh mẽ bên trong khoang liên kết của protein.

Khi so sánh với các nghiên cứu tương tự trong tài liệu quốc tế về dẫn chất curcuminoid, mức năng lượng gắn kết -11,1 kcal/mol trên thụ thể ERα là một trong những giá trị ái lực cao nhất từng được ghi nhận cho nhóm hợp chất khung lai pyrazol-monocarbonyl. Dữ liệu mô phỏng quỹ đạo động lực học phân tử có thể được biểu diễn trực quan qua biểu đồ đường biến thiên RMSD và RMSF theo thời gian 100 ns, cho thấy vùng liên kết của thụ thể không bị biến dạng bất thường, minh chứng rằng các dẫn chất đóng vai trò là chất điều hòa chọn lọc và ức chế cạnh tranh có độ bền nhiệt động học cao.

  Biến thiên RMSD (Å) trong mô phỏng MDs 100 ns:
  3.0 |                                                     
  2.5 |                                                     
  2.0 |------------------ Phức hợp CUR.5 - ERα (RMSD ~ 1.8 Å ổn định)
  1.5 |.................. Phức hợp CUR.1 - IGF-1R (RMSD ~ 1.6 Å ổn định)
  1.0 |                                                     
  0.5 |                                                     
  0.0 +----------------------------------------------------> Thời gian (ns)
      0       20       40       60       80       100       

Đề xuất và khuyến nghị

Dựa trên những phát hiện khoa học từ đề tài, các giải pháp cụ thể được đề xuất nhằm tiếp tục phát triển hướng nghiên cứu:

Tối ưu hóa quy trình tổng hợp hóa học: Bộ môn Hóa Dược và Công nghệ Dược phẩm cần tiến hành chuẩn hóa điều kiện phản ứng xúc tác base và dung môi hữu cơ thân thiện môi trường trong vòng 6 đến 12 tháng tới, nhằm mục tiêu nâng hiệu suất tổng hợp toàn phần của các dẫn chất MACs từ mức bình quân 25-77% lên đạt trên 85%.

Mở rộng thử nghiệm độc tế bào in vitro: Nhóm nghiên cứu thuộc Bộ môn Dược lực - Dược lâm sàng nên thực hiện mở rộng thử nghiệm đánh giá IC50 trên các dòng tế bào sarcoma cơ vân đặc thù như RD, RH30, cùng các dòng ung thư vú kháng nội tiết như MDA-MB-231 và T47D trong lộ trình 12 tháng, nhằm xác định phổ kháng u rộng với mục tiêu IC50 duy trì dưới 5 µM.

Đánh giá dược động học và an toàn in vivo: Các trung tâm nghiên cứu phát triển thuốc và Viện Kiểm nghiệm thuốc cần chủ trì thử nghiệm tiền lâm sàng trên mô hình động vật gặm nhấm trong thời gian từ 18 đến 24 tháng, tập trung đo lường thời gian bán thải sinh học, độc tính cấp, độc tính bán trường diễn và thể tích phân bố mô để kiểm chứng tính an toàn của dẫn chất CUR.5 và CUR.1.

Phát triển hệ phân phối thuốc dạng nano: Các chuyên gia công nghệ bào chế cần nghiên cứu thiết kế các hệ mang thuốc kích thước nano như hạt nano lipid tự phân tán (SEDDS), liposome hoặc micelle polyme trong vòng 12 đến 18 tháng, nhằm nâng cao độ tan trong nước và tăng tỷ lệ hấp thu sinh khả dụng của hoạt chất lên trên 90%.

Đối tượng nên tham khảo luận văn

Nhà nghiên cứu Hóa Dược và Hóa Hữu cơ: Tiếp cận quy trình tổng hợp chi tiết các dị vòng pyrazol kết hợp khung curcuminoid không đối xứng, cùng bộ dữ liệu thực nghiệm phổ 1H-NMR, 13C-NMR, MS và IR để phục vụ việc thiết kế và tổng hợp các dẫn chất kháng ung thư mới.

Chuyên viên mô hình hóa phân tử và thiết kế thuốc tính toán (CADD): Ứng dụng quy trình sàng lọc in silico chuẩn hóa từ bước chuẩn bị cấu trúc, docking phân tử trên AutoDock/Vina cho đến phân tích động lực học phân tử 100 ns trên GROMACS với các chỉ số đánh giá nhiệt động học chuyên sâu.

Dược sĩ và nhà phát triển công nghệ bào chế: Tham khảo các đặc tính lý hóa, độ kỵ nước và cấu trúc phân tử của dẫn chất monocarbonyl curcumin để định hướng nghiên cứu công thức bào chế các dạng thuốc giải phóng có kiểm soát và hệ dẫn thuốc hướng đích.

Học viên sau đại học và sinh viên ngành Dược học/Y sinh: Sử dụng làm tài liệu tham khảo học thuật giá trị về phương pháp luận kết hợp liên ngành giữa tổng hợp hóa dược, thử nghiệm sinh học tế bào và mô phỏng tin sinh học trong quy trình khám phá thuốc hiện đại.

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

Dẫn chất monocarbonyl curcumin (MACs) có ưu điểm gì vượt trội so với curcumin tự nhiên? Curcumin tự nhiên sở hữu cấu trúc β-diketon rất dễ bị thủy phân ở pH sinh lý, chuyển hóa nhanh tại gan và có sinh khả dụng đường uống dưới 1%. Việc chuyển đổi sang cấu trúc monocarbonyl giúp loại bỏ vị trí kém bền, gia tăng độ ổn định hóa học trong môi trường sinh học và tăng hoạt tính gây độc tế bào ung thư lên gấp 10 đến 20 lần.

Dị vòng 1H-pyrazol đóng vai trò then chốt như thế nào trong hoạt tính kháng ung thư? Dị vòng 1H-pyrazol là một khung cấu trúc mang tính sinh học cao, cung cấp các nguyên tử nitơ giàu điện tử để tạo liên kết hydro bền vững với acid amin trong khoang thụ thể. Cấu trúc này không chỉ mở rộng khả năng kẹp giữ trong túi liên kết của ERα và IGF-1R mà còn góp phần ức chế tăng sinh mạch máu và kích hoạt con đường apoptosis của tế bào ác tính.

Tại sao nghiên cứu lại lựa chọn hai thụ thể ERα và IGF-1R làm đích tác động chính? Thụ thể Estrogen Alpha biểu hiện quá mức trong 50% đến 80% trường hợp ung thư vú, đóng vai trò sống còn trong điều hòa phiên mã tăng sinh khối u. Trong khi đó, IGF-1R hiện diện ở hơn 80% ca bệnh sarcoma cơ vân, chịu trách nhiệm truyền tín hiệu sinh tồn tế bào. Tấn công có định hướng vào hai thụ thể này đem lại hiệu quả điều trị trúng đích cao.

Ý nghĩa thực tiễn của giai đoạn mô phỏng động lực học phân tử (MDs) 100 ns là gì? Quá trình gắn kết phân tử tĩnh chỉ dự đoán tư thế gắn tại một thời điểm cố định. Mô phỏng động lực học 100 ns trên GROMACS cho phép khảo sát sự chuyển động của phức hợp trong môi trường nước sinh lý thực tế, chứng minh phức hợp giữa dẫn chất CUR.5 với ERα và CUR.1 với IGF-1R có dao động RMSD dưới 2,0 Å, khẳng định liên kết duy trì cực kỳ ổn định theo thời gian.

Phương pháp thử nghiệm MTT đánh giá hoạt tính độc tế bào có độ tin cậy như thế nào? Thử nghiệm MTT là phương pháp chuẩn quốc tế đo lường hoạt tính enzym dehydrogenase trong ty thể của các tế bào sống, giúp chuyển hóa muối tetrazolium màu vàng thành tinh thể formazan màu tím. Độ hấp thụ quang học tỷ lệ thuận với số lượng tế bào sống sót, cho phép xác định chính xác giá trị IC50 của 36 hợp chất trên dòng tế bào MCF-7 với độ lặp lại cao và sai số thực nghiệm thấp.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu thông qua 5 đóng góp khoa học nổi bật:

  • Tổng hợp và kiểm chứng thành công cấu trúc hóa học của 9 dẫn chất monocarbonyl curcumin không đối xứng chứa dị vòng 1H-pyrazol với hiệu suất phản ứng đạt tới 77%.
  • Xác định được 5 dẫn chất có hoạt tính độc tế bào vượt trội trên dòng tế bào ung thư vú MCF-7 với nồng độ ức chế IC50 ấn tượng trong khoảng 2,55 µM đến 13,49 µM.
  • Hoàn thành sàng lọc in silico trên 111 cấu trúc, phát hiện 9 dẫn chất gắn kết mạnh trên thụ thể ERα (-11,1 đến -9,8 kcal/mol) và 7 dẫn chất gắn kết tốt trên thụ thể IGF-1R (-9,6 đến -9,0 kcal/mol).
  • Chứng minh độ bền liên kết và tính ổn định nhiệt động học của hai hợp chất tiềm năng CUR.5 (trên ERα) và CUR.1 (trên IGF-1R) qua mô phỏng động lực học phân tử 100 ns.
  • Thiết lập quy trình nghiên cứu tích hợp hoàn chỉnh và hiệu quả giữa hóa dược tổng hợp, thử nghiệm sinh học in vitro và mô hình hóa phân tử in silico.

Lộ trình nghiên cứu tiếp theo trong giai đoạn 12 đến 24 tháng tới sẽ tập trung vào thử nghiệm độc tính in vivo và xây dựng hệ phân phối thuốc nano. Quý độc giả và các nhà khoa học quan tâm được khuyến khích tham khảo chi tiết toàn văn luận văn thạc sĩ để khai thác toàn bộ dữ liệu phổ cấu trúc và thông số mô phỏng phân tử phục vụ cho công tác nghiên cứu phát triển thuốc.