Nghiên cứu cấu trúc và khả năng chống oxy hóa của polyphenol dẫn xuất trên nền fullerene

Tổng quan về luận án

Quá trình suy thoái tế bào dưới tác động của stress oxy hóa là căn nguyên hàng đầu dẫn đến các bệnh lý mãn tính nguy hiểm ở người như ung thư, xơ vữa động mạch, rối loạn chuyển hóa và các hội chứng thoái hóa thần kinh. Cụ thể, các gốc tự do và tác nhân phản ứng chứa oxy (Reactive Oxygen Species - ROS), điển hình như superoxide ($O_2^{\bullet-}$), hydroxyl ($^\bullet OH$) và peroxyl ($ROO^\bullet$), làm biến tính protein, peroxy hóa màng lipid tế bào, hình thành sản phẩm cộng DNA và làm bất hoạt các enzyme thiết yếu. Nhằm ức chế chuỗi phản ứng hủy hoại này, các chất chống oxy hóa tự nhiên như polyphenol từ thực vật đóng vai trò bẫy gốc tự do đặc biệt quan trọng. Song song đó, vật liệu nano carbon fullerene ($C_{60}$) nổi lên như một "bọt biển quét gốc tự do" (radical sponge) với ái lực electron mạnh mẽ. Tuy nhiên, $C_{60}$ nguyên bản bộc lộ hạn chế cố hữu: cấu trúc không thể hoạt động theo cơ chế ngắt mạch dây chuyền đối với gốc peroxyl ($ROO^\bullet$).

Khoảng trống nghiên cứu (Research Gap) then chốt được xác định: Trong khi các nguồn dược liệu tự nhiên phong phú của Việt Nam như cây sa kê (Artocarpus altilis) và vỏ quả măng cụt (Garcinia mangostana) đã được ghi nhận có hoạt tính sinh học qua các thử nghiệm in vitro/in vivo sơ bộ, thì bản chất cơ chế phản ứng ở mức độ lượng tử, các thông số nhiệt động học và động học của các polyphenol này vẫn chưa được giải thích tường minh. Đồng thời, trên bình diện quốc tế, việc thiết kế lai hóa cấu trúc giữa polyphenol tự nhiên và khung fullerene ($C_{60}$) nhằm tạo ra các phân tử chống oxy hóa đa cơ chế vẫn còn là một địa hạt sơ khai chưa có hệ thống tính toán lượng tử toàn diện.

Luận án của nghiên cứu sinh Nguyễn Minh Thông, dưới sự hướng dẫn khoa học của PGS. TS. Phạm Cẩm Nam và PGS. TS. Trần Dương tại Trường Đại học Khoa học - Đại học Huế (Chuyên ngành Hóa lý thuyết và Hóa lý, Mã số: 62.19), đã giải quyết triệt để khoảng trống này bằng công cụ hóa học tính toán đa quy mô.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác lập:

• Câu hỏi nghiên cứu 1 (RQ1): Các phân tử polyphenol từ cây sa kê và xanthone từ vỏ măng cụt dập tắt gốc tự do ưu tiên theo cơ chế nhiệt động nào (HAT, SET-PT hay SPLET) trong pha khí và môi trường dung môi?
• Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc hình học, trạng thái chuyển tiếp và hàng rào năng lượng của phản ứng dập tắt gốc tự do peroxyl ($CH_3OO^\bullet$) diễn ra như thế nào?
• Câu hỏi nghiên cứu 3 (RQ3): Cơ chế phản ứng cộng đóng vòng Bingel - Hirsch giữa dẫn xuất malonate của polyphenol và khung $C_{60}$ diễn ra theo lộ trình lập thể nào để đạt độ bền nhiệt động tối ưu?
• Câu hỏi nghiên cứu 4 (RQ4): Việc lai hóa polyphenol lên khung $C_{60}$ có thực sự cải thiện vượt bậc khả năng chống oxy hóa so với các phân tử tự nhiên ban đầu hay không?
• Giả thuyết 1 (H1): Phản ứng dập tắt gốc tự do của polyphenol trong pha khí bị chi phối chủ yếu bởi cơ chế chuyển nguyên tử hydro (HAT) với giá trị BDE(O-H) đóng vai trò quyết định, trong khi dung môi phân cực kích hoạt cơ chế mất proton tuần tự chuyển electron (SPLET).
• Giả thuyết 2 (H2): Phản ứng Bingel - Hirsch ưu tiên diễn ra tại liên kết cầu nối vòng (6,6) của khung $C_{60}$ do mật độ electron $\pi$ cục bộ và sức căng góc vòng tối ưu.
• Giả thuyết 3 (H3): Cấu trúc lai $C_{60}$-polyphenol thể hiện hiệu ứng cộng hưởng hiệp đồng: khung $C_{60}$ thu hút gốc carbon trong khi nhóm polyphenol đính kèm bẻ gãy chuỗi peroxyl qua cơ chế HAT.

Khung lý thuyết của công trình tích hợp Lý thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT) với phiếm hàm lai B3LYP, phương pháp bán thực nghiệm PM6 và mô hình tích hợp đa lớp ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) do Morokuma phát triển. Đối tượng khảo sát bao gồm 12 hợp chất polyphenol từ sa kê, 14 hợp chất xanthone từ vỏ măng cụt, các flavonoid chuẩn (quercetin, catechin, chalcone, flavone, flavanone) và các hệ lai nano fullerene-polyphenol (như $C_{60}$-altilisin J và $C_{60}$-mangostin).

                      [ BỨC TRANH NGHIÊN CỨU TỔNG QUAN ]
                                      │
            ┌─────────────────────────┴─────────────────────────┐
            ▼                                                   ▼
 ┌───────────────────────┐                           ┌─────────────────────┐
 │ Hợp chất Polyphenol   │                           │ Nanocarbon          │
 │ Tự nhiên (Sa kê,      │                           │ Fullerene (C60)     │
 │ Vỏ măng cụt)          │                           │                     │
 └──────────┬────────────┘                           └──────────┬──────────┘
            │ BDE, IE, PA, ETE, PDE                             │ Sức căng sp2, LUMO
            ▼                                                   ▼
 ┌─────────────────────────────────────────────────────────────────────────┐
 │               Mô phỏng Đa tầng ONIOM (DFT : PM6) & IEFPCM               │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
                                      ▼
 ┌─────────────────────────────────────────────────────────────────────────┐
 │             Phản ứng Đóng vòng Bingel - Hirsch qua TS(6,6)              │
 └────────────────────────────────────┬────────────────────────────────────┘
                                      │
                                      ▼
 ┌─────────────────────────────────────────────────────────────────────────┐
 │  Hệ Dẫn xuất Lai hóa Đột phá: C60-Altilisin J & C60-Mangostin           │
 │  (Tối ưu hóa BDE, IE, EA -> Dập tắt ROS đa cơ chế)                     │
 └─────────────────────────────────────────────────────────────────────────┘

Literature Review và Positioning

Nghiên cứu về cơ chế dập tắt gốc tự do của polyphenol và hoạt tính sinh học của fullerene là giao điểm của hóa học lượng tử, hóa sinh và khoa học vật liệu nano. Trong nước, nhóm tác giả Phạm Cẩm Nam và cộng sự (2012, 2013) đã tiên phong áp dụng mô hình đa lớp ONIOM(ROB3LYP/6-311++G(2df,2p):AM1) để tính toán enthalpy phân ly liên kết BDE(O-H) cho epicatechin (EC) và epigallocatechin gallate (EGCG) từ chè xanh và gừng, khẳng định ưu thế của việc mô phỏng cơ chế HAT trong việc giải thích hoạt tính sinh học.

Trên bình diện quốc tế, các luồng nghiên cứu chính phân hóa thành các nhánh tiếp cận chuyên biệt:

1. Nghiên cứu cấu trúc - hoạt tính polyphenol tự nhiên: Nenadis và cộng sự (2007) khảo sát khả năng dập tắt gốc tự do của dẫn xuất 3',4'-dihydroxy aurones bằng phương pháp DFT/B3LYP thông qua hai tham số BDE(O-H) và năng lượng ion hóa (IE). Mikulski và cộng sự (2010) làm sáng tỏ mối quan hệ giữa cấu trúc lập thể và hoạt tính chống oxy hóa của trans-resveratrol từ đậu phộng dựa trên phân bố mật độ spin và HOMO-LUMO. Xue và cộng sự (2014) khảo sát deoxybenzoins và kết luận cơ chế HAT chiếm ưu thế tuyệt đối trong pha khí, nhưng cơ chế SPLET lại kiểm soát quá trình trong môi trường nước phân cực.
2. Tranh luận khoa học về cơ chế chống oxy hóa: Một cuộc tranh luận kéo dài trong y sinh lý thuyết liên quan đến việc cơ chế nào đóng vai trò cốt lõi trong tế bào sống: HAT (phân cắt đồng ly liên kết O-H) hay SET-PT (chuyển electron đơn sau đó chuyển proton) hay SPLET (tách proton tuần tự rồi chuyển electron). Tranh luận thứ hai liên quan đến hạt nhân $C_{60}$: trong khi Kroto, Smalley và Curl (1985) phát hiện cấu trúc lồng $C_{60}$ với sức căng góc chóp $\theta_p = 101.6^\circ$ đem lại ái lực electron (EA) cực đại để bẫy gốc tự do alkyl, thì các nghiên cứu của Yin, Pratt và cộng sự chỉ ra rằng $C_{60}$ bất hoạt trước các gốc peroxyl ($ROO^\bullet$) trong giai đoạn lan truyền mạch peroxy hóa lipid.

                  [ VỊ TRÍ HỌC THUẬT CỦA LUẬN ÁN ]
                         
  Nghiên cứu Polyphenol Đơn thuần             Nghiên cứu Fullerene C60
  (Nenadis 2007, Mikulski 2010,               (Kroto 1985, Hirsch 1993,
   Pham Cam Nam 2012, Xue 2014)                Yin & Pratt 2000s)
               │                                      │
               │                                      │
               └───────────────────┬──────────────────┘
                                   │
                                   ▼
        ┌────────────────────────────────────────────────────┐
        │        VỊ TRÍ TIÊN PHONG CỦA LUẬN ÁN               │
        │  Thiết kế Hệ Lai Hóa Nano C60 - Polyphenol         │
        │  (Artocarpus altilis & Garcinia mangostana)        │
        │  Đa cơ chế: Bẫy Carbon Radical + Ngắt mạch Peroxyl │
        └────────────────────────────────────────────────────┘

Luận án của Nguyễn Minh Thông định vị chính xác tại điểm hội tụ này: Sử dụng phương pháp hóa tính toán lượng tử hiện đại để phân tích toàn diện 3 cơ chế (HAT, SET-PT, SPLET) trên các hợp chất tự nhiên chiết xuất từ sa kê và măng cụt Việt Nam, đồng thời mô phỏng phản ứng Bingel - Hirsch để gắn kết các polyphenol này lên khung $C_{60}$, khắc phục triệt để nhược điểm trơ với gốc peroxyl của fullerene nguyên bản.

Đó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 và làm phong phú thêm các học thuyết nền tảng của hóa lý lượng tử và động hóa học:

1. Mở rộng Thuyết Phiếm hàm Mật độ (Hohenberg-Kohn-Sham DFT) và Thuyết Orbital Biên tử (Frontier Molecular Orbital Theory - FMO): Chứng minh trực tiếp mối tương quan giữa mức năng lượng HOMO, LUMO, độ biến thiên spin và mật độ phân bố điện tích đối với khả năng cho proton/electron của phân tử.
2. Khẳng định nguyên lý Evans - Polanyi trên phản ứng dập tắt gốc peroxyl: Tìm thấy sự phụ thuộc tuyến tính chặt chẽ giữa enthalpy phân ly liên kết $\text{BDE(O-H)}$ và năng lượng hoạt hóa ($E_a$) của trạng thái chuyển tiếp giữa flavonoid và gốc tự do $CH_3OO^\bullet$.
3. Làm sâu sắc thêm Thuyết Trạng thái Chuyển tiếp (Transition State Theory) trong phản ứng Bingel - Hirsch: Thiết lập cấu hình năng lượng bề mặt thế năng (PES), chứng minh cơ chế đóng vòng cyclopropane diễn ra qua 2 giai đoạn: cộng carbanion nucleophile vào $C_{60}$ tạo anion trung gian, sau đó bứt tách $Br^-$ nội phân tử thông qua trạng thái chuyển tiếp đồng phẳng $TS_{(6,6)b1}$.

Hệ thống mệnh đề khoa học được xây dựng:

• Mệnh đề 1 (P1): Giá trị BDE(O-H) của nhóm hydroxyl là đại lượng chỉ thị nhiệt động quyết định hoạt tính chống oxy hóa trong pha khí; các nhóm hydroxyl có liên kết hydro nội phân tử hoặc cộng hưởng ortho-diphenol (như trong quercetin hoặc dẫn xuất xanthone) sở hữu BDE thấp hơn đáng kể ($< 85\text{ kcal/mol}$).
• Mệnh đề 2 (P2): Sự chuyển dịch cơ chế từ HAT sang SPLET trong dung môi nước được điều khiển bởi ái lực proton (PA) giảm mạnh do hằng số điện môi cao ($\varepsilon = 78.39$) làm bền hóa anion $ArO^-$.
• Mệnh đề 3 (P3): Phản ứng cộng phối vị cyclopropan hóa Bingel - Hirsch có tính chọn lọc vị trí tuyệt đối tại cầu nối liên kết (6,6) do độ dài liên kết ngắn hơn ($1.391\text{ \AA}$) và tính chất giống liên kết đôi so với cầu nối (5,6) ($1.449\text{ \AA}$).
• Mệnh đề 4 (P4): Sự tích hợp dị thể giữa fullerene $C_{60}$ và polyphenol tạo nên hiệu ứng chuyển điện tích nội phân tử (Intramolecular Charge Transfer), gia tăng ái lực electron (EA) và đồng thời hạ thấp thế ion hóa (IE).

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

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết: Lý thuyết Cơ học Lượng tử cấu trúc electron, Lý thuyết Nhiệt động học thống kê và Mô hình Solvat hóa liên tục. Phương pháp phân chia phân tử đa tầng ONIOM 2 lớp được tối ưu hóa đặc thù cho các hệ phân tử sinh học - nano cồng kềnh: $$\text{E}{\text{ONIOM2(high:low)}} = \text{E}{\text{high,model}} + \text{E}{\text{low,real}} - \text{E}{\text{low,model}}$$

Trong đó:

• Lớp mức cao (High Layer - Model System): Chứa trung tâm phản ứng trực tiếp (vòng thơm mang nhóm thế -OH hoặc liên kết cyclopropane đính $C_{60}$) được tính toán chuẩn xác bằng phiếm hàm lai mức cao $\text{ROB3LYP/6-311++G(2df,2p)}$.
• Lớp mức thấp (Low Layer - Real System): Toàn bộ khung phân tử còn lại (bao gồm 60 nguyên tử carbon của $C_{60}$ và các mạch nhánh cồng kềnh) được xử lý bằng phương pháp bán thực nghiệm PM6 nhằm tối ưu hóa tài nguyên tính toán mà vẫn đảm bảo độ chính xác vật lý.
• Điều kiện biên: Nhiệt độ tiêu chuẩn $298.15\text{ K}$, áp suất $1.0\text{ atm}$, các trạng thái dừng được xác nhận không có tần số dao động ảo, trạng thái chuyển tiếp (TS) có duy nhất một tần số dao động âm và được kiểm chứng bằng tọa độ phản ứng nội (IRC).

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng luận (Positivism) kết hợp phương pháp diễn dịch định lượng thông qua cơ học lượng tử tính toán (Computational Quantum Mechanics). Thiết kế nghiên cứu kết hợp đa cấp độ từ bán thực nghiệm (PM6) đến phiếm hàm mật độ chuẩn cao (DFT/B3LYP) và tiếp cận lai đa mức (ONIOM2).

                      [ QUY TRÌNH TÍNH TOÁN ĐA BƯỚC ]
                                      │
            ┌─────────────────────────┴─────────────────────────┐
            ▼                                                   ▼
┌───────────────────────┐                           ┌───────────────────────┐
│ Xây dựng cấu trúc     │                           │ Khảo sát BDE & IE     │
│ 3D trên GaussView 5   │                           │ Chuẩn hóa trên Phenol │
└───────────┬───────────┘                           └───────────┬───────────┘
            │                                                   │
            └─────────────────────────┬─────────────────────────┘
                                      │
                                      ▼
            ┌───────────────────────────────────────────────────┐
            │ Tối ưu hóa Hình học & Tần số dao động             │
            │ (PM6, Opt + Freq, Hệ số ZPE = 1.078)              │
            │ Chuẩn hội tụ: Max Force < 0.00045 Hartree/Bohr    │
            └─────────────────────────┬─────────────────────────┘
                                      │
                                      ▼
            ┌───────────────────────────────────────────────────┐
            │ Xác định Trạng thái Chuyển tiếp & Đường IRC       │
            │ [opt=(ts, calcfc, noeigen), stepsize = 10]        │
            └─────────────────────────┬─────────────────────────┘
                                      │
                                      ▼
            ┌───────────────────────────────────────────────────┐
            │ Tính Năng lượng Điểm đơn Single Point:            │
            │ ONIOM(ROB3LYP/6-311++G(2df,2p):PM6)               │
            │ Mô hình Dung môi IEFPCM (Gas & Water)             │
            └───────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

1. Xây dựng cấu trúc và Tối ưu hóa hình học: Các phân tử ban đầu được khởi tạo trên phần mềm GaussView 5. Tối ưu hóa cấu trúc hình học ở trạng thái cơ bản và tính toán tần số dao động được thực hiện tại mức lý thuyết PM6 trong phần mềm Gaussian 09W.
2. Tiêu chuẩn hội tụ khắt khe (Default Convergence Criteria):
   • Giá trị lực cực đại tác dụng lên hạt nhân ($\text{Max Force}$) $< 0.00045\text{ Hartree/Bohr}$.
   • Độ lệch quân phương của lực ($\text{RMS Force}$) $< 0.00030\text{ Hartree/Bohr}$.
   • Độ dịch chuyển cực đại của nguyên tử ($\text{Max Displacement}$) $< 0.00180\text{ Bohr}$.
   • Độ dịch chuyển quân phương ($\text{RMS Displacement}$) $< 0.00120\text{ Bohr}$.
3. Xác định Trạng thái chuyển tiếp (Transition State - TS) và Tọa độ phản ứng nội (IRC):
   • Tìm kiếm điểm yên ngựa bậc một (first-order saddle point) bằng từ khóa opt=(ts, calcfc, noeigen) kết hợp kiểm tra ma trận lực Hessen.
   • Kiểm tra tần số dao động chuẩn: Trạng thái chuyển tiếp được xác thực khi và chỉ khi có duy nhất một tần số dao động âm (tần số ảo - imaginary frequency).
   • Phân tích đường phản ứng nội: Thực hiện tính toán IRC với thông số maxpoints = 10 đến 20 và stepsize = 10 (tương đương $0.1\text{ amu}^{1/2}\text{ Bohr}$) nhằm nối liền hai cực tiểu năng lượng của chất phản ứng và sản phẩm.
4. Hiệu chỉnh nhiệt động học: Năng lượng dao động điểm không (Zero Point Energy - ZPE) ở mức lý thuyết PM6 được nhân hệ số hiệu chỉnh thực nghiệm là $1.078$. Các đại lượng Enthalpy ($H$) và Năng lượng tự do Gibbs ($G$) được tính toán ở điều kiện chuẩn $298.15\text{ K}$ và $1\text{ atm}$ thông qua hiệu chỉnh TCE (Thermal Correction to Enthalpy) và TCGFE (Thermal Correction to Gibbs Free Energy).

Data và phân tích các đại lượng nhiệt động

Khả năng chống oxy hóa được định lượng thông qua 5 thông số nhiệt động cơ bản theo các phương trình toán học lượng tử:

• Cơ chế chuyển nguyên tử hydro (HAT): $$\text{BDE(O-H)} = H(\text{ArO}^\bullet) + H(\text{H}^\bullet) - H(\text{ArOH})$$
• Cơ chế chuyển một electron chuyển proton (SET-PT): $$\text{IE} = H(\text{ArOH}^{+\bullet}) + H(e^-) - H(\text{ArOH})$$ $$\text{PDE} = H(\text{ArO}^\bullet) + H(\text{H}^+) - H(\text{ArOH}^{+\bullet})$$
• Cơ chế tuần tự mất proton chuyển electron (SPLET): $$\text{PA} = H(\text{ArO}^-) + H(\text{H}^+) - H(\text{ArOH})$$ $$\text{ETE} = H(\text{ArO}^\bullet) + H(e^-) - H(\text{ArO}^-)$$

Ảnh hưởng của môi trường dung môi nước được mô phỏng thông qua mô hình phân cực liên tục dựa trên phương trình tích phân IEFPCM (Integral Equation Formalism Polarizable Continuum Model).

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

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

1. Sàng lọc các chất chống oxy hóa tự nhiên vượt trội từ cây Sa kê và Măng cụt:
   • Trong 12 hợp chất polyphenol từ sa kê, hợp chất S12 thể hiện hoạt tính bẻ gãy gốc tự do mạnh nhất với giá trị BDE(O-H) trong pha khí đạt mức kỷ lục thấp: $77.3\text{ kcal/mol}$.
   • Trong 14 hợp chất xanthone từ vỏ quả măng cụt, hợp chất M10 và M11 thể hiện hoạt tính chống oxy hóa xuất sắc với giá trị BDE(O-H) lần lượt là $82.3\text{ kcal/mol}$ và $82.8\text{ kcal/mol}$. Các giá trị này vượt trội hoặc tương đương với các chất chống oxy hóa chuẩn như $\alpha$-tocopherol ($78.2\text{ kcal/mol}$) và quercetin ($81.5\text{ kcal/mol}$).
2. Quy luật động học phản ứng dập tắt gốc peroxyl ($CH_3OO^\bullet$): Khảo sát phản ứng dập tắt gốc $CH_3OO^\bullet$ bởi quercetin và các polyphenol cho thấy phản ứng diễn ra qua phức trung gian tiền phản ứng ($\text{Int-1}$), vượt qua trạng thái chuyển tiếp (TS) với liên kết $O\cdots H\cdots O$ kéo dài, tạo phức sau phản ứng ($\text{Int-2}$) và phân ly thành sản phẩm. Tồn tại một sự tương quan tuyến tính chặt chẽ giữa BDE(O-H) và hàng rào thế năng hoạt hóa ($E_a$), chứng minh BDE là thông số quyết định tốc độ phản ứng dập tắt gốc tự do.
3. Bản chất nhiệt động và lập thể của phản ứng Bingel - Hirsch trên nền $C_{60}$:
   • Nghiên cứu phản ứng giữa ion âm dimethyl bromomalonate $[CBr(COOCH_3)2]^-$ và $C{60}$ chứng minh trong 4 đường phản ứng khả dĩ, lộ trình đi qua trạng thái chuyển tiếp $TS_{(6,6)b1}$ để tạo vòng cyclopropane tại vị trí liên kết (6,6) có hàng rào năng lượng thấp nhất và sản phẩm có năng lượng tự do Gibbs $\Delta G$ âm nhất, hoàn toàn phù hợp với thực nghiệm tổng hợp hữu cơ.
   • Khoảng cách liên kết carbon-carbon tại cầu nối (6,6) sau khi đóng vòng giãn ra từ $1.391\text{ \AA}$ lên xấp xỉ $1.58\text{ \AA}$, giải tỏa sức căng hình học của lồng fullerene.
4. Thiết kế thành công hệ vật liệu lai nano $C_{60}$-Polyphenol thế hệ mới: Các dẫn xuất lai hóa giữa fullerene và altilisin J ($C_{60}\text{-altilisin J}$) cũng như mangostin ($C_{60}\text{-mangostin-1}$, $C_{60}\text{-mangostin-2}$) đã được thiết kế và tính toán cấu trúc thành công. Hệ lai nano này không chỉ bảo tồn nguyên vẹn giá trị BDE(O-H) thấp của polyphenol mà còn tối ưu hóa ái lực electron (EA) và năng lượng ion hóa (IE), biến phân tử lai thành chất chống oxy hóa lưỡng tính cực mạnh: vừa bẫy gốc alkyl vừa dập tắt gốc peroxyl.

Implications đa chiều

• Về mặt học thuật và lý thuyết: Giải quyết dứt điểm các tranh luận về cơ chế ưu tiên của polyphenol trong các môi trường điện môi khác nhau; bổ sung cơ sở dữ liệu lượng tử chuẩn xác cho các hợp chất tự nhiên đặc hữu Đông Nam Á.
• Về mặt phương pháp luận: Chuẩn hóa quy trình tính toán ONIOM 2 lớp $\text{ROB3LYP/6-311++G(2df,2p):PM6}$ như một công cụ chuẩn mực để nghiên cứu các hệ lai vật liệu nano - phân tử hữu cơ lớn, giúp tiết kiệm hàng ngàn giờ tính toán trên siêu máy tính mà vẫn giữ sai số năng lượng $< 1.5\text{ kcal/mol}$ so với thực nghiệm.
• Về mặt ứng dụng thực tiễn và công nghiệp: Cung cấp bản thiết kế phân tử chính xác cho ngành dược phẩm và mỹ phẩm để định hướng bán tổng hợp các hoạt chất chống lão hóa, bảo vệ tế bào thần kinh và hỗ trợ điều trị ung thư từ nguồn phụ phẩm nông nghiệp (vỏ măng cụt, lá sa kê).

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận một số giới hạn phương pháp luận:

1. Mô hình dung môi liên tục IEFPCM: Việc sử dụng mô hình điện môi liên tục chưa mô tả được chi tiết các liên kết hydro đặc hiệu giữa phân tử nước riêng lẻ với nhóm -OH của polyphenol, điều có thể làm dịch chuyển nhẹ giá trị năng lượng phân ly proton (PDE).
2. Phương pháp bán thực nghiệm PM6 cho khung $C_{60}$: Dù PM6 tái hiện tốt cấu trúc hình học, việc tính toán năng lượng ion hóa (IE) ở mức PM6 đôi khi có độ lệch hệ thống nhỏ so với các phép tính phiếm hàm thuần túy với hàm cơ sở mở rộng khuếch tán.
3. Thiếu vắng dữ liệu động học mô phỏng động lực học phân tử (MD): Nghiên cứu dừng lại ở mức khảo sát bề mặt thế năng tĩnh (static PES) ở $0\text{ K}$ và hiệu chỉnh nhiệt động học ở $298.15\text{ K}$, chưa khảo sát hành vi phân tử trong màng lipid kép phospholipid qua mô phỏng MD.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

• Hướng 1: Ứng dụng các hàm phiếm hàm có hiệu chỉnh phân tán tầm xa (DFT-D3, $\omega\text{B97X-D}$) kết hợp mô phỏng solvat hóa tường minh (explicit solvent with water clusters).
• Hướng 2: Thực hiện mô phỏng Động lực học Phân tử (Molecular Dynamics - MD) để khảo sát sự xâm nhập và định hướng của $C_{60}$-mangostin xuyên qua màng tế bào ung thư.
• Hướng 3: Triển khai thực nghiệm tổng hợp hóa học (Organic Synthesis) dẫn xuất $C_{60}$-altilisin J theo quy trình Bingel - Hirsch và thử nghiệm độc tính tế bào in vitro.
• Hướng 4: Mở rộng nghiên cứu khả năng bẫy các gốc tự do nguy hiểm khác như gốc tự do nitrogen hoạt tính (Reactive Nitrogen Species - RNS, điển hình là $NO^\bullet$, $ONOO^-$).

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

Các kết quả nghiên cứu của luận án đã được công bố trên 5 bài báo quốc tế thuộc danh mục ISI uy tín (Journal of Molecular Modeling, Chemical Physics Letters, Springer International Publishing) và 3 bài báo trên các tạp chí khoa học chuyên ngành quốc gia (Tạp chí KH&CN Viện Hàn lâm KH&CN Việt Nam, ĐH Đà Nẵng, ĐH Huế).

• Tác động học thuật: Dự kiến tạo ra làn sóng trích dẫn cao trong lĩnh vực hóa sinh tính toán và hóa học nano; mở ra mô hình chuẩn cho việc sàng lọc ảo (virtual screening) dược liệu nhiệt đới.
• Tác động kinh tế - xã hội: Thúc đẩy kinh tế tuần hoàn và công nghệ xanh tại Việt Nam thông qua việc nâng cao chuỗi giá trị của cây sa kê và tái sử dụng vỏ quả măng cụt - vốn là phế phẩm nông nghiệp sau chế biến - thành nguyên liệu thô để chiết xuất các xanthone chống oxy hóa cao cấp.
• Chiến lược công nghệ y dược: Mở đường cho các liệu pháp điều trị hướng đích sử dụng dẫn xuất fullerene hòa tan trong nước làm chất ức chế enzyme thoái hóa thần kinh (Alzheimer, Parkinson) và chất cảm quang trong xạ trị khối u.

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

• Nghiên cứu sinh và Giới học thuật lý thuyết: Tiếp cận quy trình tính toán đa lớp ONIOM chuẩn hóa, bộ thông số nhiệt động lượng tử chi tiết và phương pháp phân tích trạng thái chuyển tiếp phản ứng gốc tự do.
• Các nhà Hóa dược và Hóa tổng hợp hữu cơ: Nhận được các chỉ dẫn lập thể và điều kiện nhiệt động chính xác để tối ưu hóa hiệu suất phản ứng Bingel - Hirsch khi tổng hợp dẫn xuất fullerene.
• Doanh nghiệp Dược - Mỹ phẩm: Sở hữu cơ sở khoa học vững chắc để phát triển các dòng sản phẩm thực phẩm chức năng chống oxy hóa, mỹ phẩm chống lão hóa từ chiết xuất chuẩn hóa sa kê và măng cụt.
• Các nhà hoạch định chính sách KH&CN: Có thêm bằng chứng thực tiễn về hiệu quả đầu tư nghiên cứu cơ bản định hướng ứng dụng từ Quỹ Phát triển Khoa học và Công nghệ Quốc gia (NAFOSTED).

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ì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc tích hợp thành công Thuyết Phiếm hàm Mật độ (DFT) với Mô hình Đa tầng ONIOM để giải mã trọn vẹn 3 cơ chế cạnh tranh (HAT, SET-PT, SPLET) trên hệ lai phân tử nano $C_{60}$-polyphenol. Công trình mở rộng Thuyết Trạng thái Chuyển tiếp của Eyring và nguyên lý Evans - Polanyi vào lĩnh vực động học dập tắt gốc tự do peroxyl trên bề mặt lồng carbon cầu.

2. Đổi mới phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế trước đây? Trả lời: So với nghiên cứu của Nenadis et al. (2007) chỉ dừng ở mức tính DFT thuần túy cho phân tử nhỏ aurones hay Xue et al. (2014) khảo sát deoxybenzoins, luận án đã tiên phong ứng dụng mô hình 2 lớp $\text{ONIOM(ROB3LYP/6-311++G(2df,2p):PM6)}$ kết hợp hiệu chỉnh IEFPCM và ZPE scale ($1.078$). Cách tiếp cận này cho phép khảo sát chính xác các cấu trúc nano khổng lồ ($> 100$ nguyên tử) với độ chính xác tiệm cận thực nghiệm nhưng giảm $85%$ thời gian tính toán so với DFT toàn phần.

3. Phát hiện bất ngờ nhất từ dữ liệu tính toán lượng tử là gì? Trả lời: Việc gắn khung $C_{60}$ không làm suy giảm hoạt tính bẻ gãy liên kết O-H của nhóm polyphenol (BDE chỉ dao động trong khoảng sai số $\pm 0.8\text{ kcal/mol}$), nhưng lại làm giảm mạnh thế ion hóa (IE giảm từ $8.01\text{ eV}$ xuống $7.22\text{ eV}$) và tăng ái lực electron (EA), biến hệ lai $C_{60}$-mangostin thành một siêu chất chống oxy hóa hoạt động đồng thời trên cả pha kỵ nước lẫn màng tế bào.

4. Quy trình lặp lại (Replication Protocol) có được cung cấp minh bạch không? Trả lời: Luận án cung cấp chi tiết tuyệt đối các từ khóa tính toán Gaussian 09W (opt, freq, opt=(ts, calcfc, noeigen), IRC, scrf=iefpcm), tọa độ không gian nguyên tử, phương pháp chọn mô hình phân chia lớp ONIOM, và hệ số hiệu chỉnh dao động nhiệt, đảm bảo khả năng tái lập độc lập $100%$ trên bất kỳ hệ thống tính toán hiệu năng cao nào.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào? Trả lời: Chương trình 10 năm tập trung vào chuỗi giá trị khép kín: Mô phỏng động lực học phân tử (MD) $\rightarrow$ Bán tổng hợp hóa học dẫn xuất $C_{60}$-xanthone $\rightarrow$ Đánh giá hoạt tính kháng viêm/chống ung thư in vitro $\rightarrow$ Thử nghiệm tiền lâm sàng trên mô hình động vật in vivo.

Kết luận

1. Đã thiết lập thành công quy trình mô phỏng hóa học tính toán đa quy mô ONIOM(DFT:PM6) với độ tin cậy và độ chính xác cao để khảo sát các hệ sinh học nano phức tạp.
2. Xác định cơ chế chuyển nguyên tử hydro (HAT) giữ vai trò chủ đạo trong pha khí, trong khi cơ chế mất proton tuần tự chuyển electron (SPLET) chiếm ưu thế trong môi trường nước đối với các polyphenol từ sa kê và măng cụt.
3. Tìm ra các hợp chất chống oxy hóa tự nhiên tiềm năng vượt bậc gồm S12 ($\text{BDE} = 77.3\text{ kcal/mol}$), M10 ($\text{BDE} = 82.3\text{ kcal/mol}$) và M11 ($\text{BDE} = 82.8\text{ kcal/mol}$).
4. Làm sáng tỏ cơ chế phản ứng Bingel - Hirsch ở mức độ lượng tử, chứng minh tính chọn lọc vị trí tuyệt đối tại cầu nối liên kết (6,6) của lồng fullerene $C_{60}$ qua trạng thái chuyển tiếp $TS_{(6,6)b1}$.
5. Thiết kế thành công cấu trúc phân tử của các hệ vật liệu lai nano $C_{60}$-altilisin J và $C_{60}$-mangostin, mở ra thế hệ chất chống oxy hóa đa năng khắc phục triệt để nhược điểm của $C_{60}$ nguyên bản.
6. Mở ra 3 nhánh nghiên cứu liên ngành mới: Hóa học lượng tử dược liệu Việt Nam, Thiết kế vật liệu nano y sinh thông minh, và Công nghệ tận thu phụ phẩm nông nghiệp tạo dược phẩm cao cấp.