Tổng quan nghiên cứu

Stress oxy hóa sinh ra do sự mất cân bằng giữa các gốc tự do phản ứng chứa oxy (ROS) và hệ thống chống oxy hóa nội sinh là nguyên nhân gốc rễ dẫn đến hơn 50 loại bệnh lý mãn tính nguy hiểm, bao gồm ung thư, đái tháo đường, suy giảm thần kinh và bệnh tim mạch. Trong tự nhiên, nhóm sắc tố thực vật anthocyanin thuộc họ flavonoid bao gồm hơn 600 hợp chất khác nhau, nổi bật nhất là Cyanidin-3-O-β-Glucoside (C3G) với cấu trúc khung carbon C6–C3–C6 sở hữu hoạt tính chống oxy hóa vượt trội hơn nhiều hợp chất cùng nhóm. Mặc dù C3G được ứng dụng rộng rãi trong dược phẩm dinh dưỡng và công nghiệp thực phẩm, nhưng việc xác định chính xác các dạng tồn tại cũng như cơ chế dập tắt gốc tự do ở mức độ phân tử bằng phương pháp thực nghiệm gặp rất nhiều rào cản do các cân bằng chuyển hóa diễn ra nhanh và phức tạp.

Luận văn thạc sĩ chuyên ngành Hóa lý thuyết và Hóa lý (mã số: 8440119) do học viên Tăng Thị Thúy Vi thực hiện dưới sự hướng dẫn khoa học của TS. Mai Văn Bảy tại Trường Đại học Sư phạm – Đại học Đà Nẵng vào năm 2024 đã giải quyết triệt để bài toán này. Mục tiêu cụ thể của đề tài là khảo sát hệ thống 28 cân bằng điện ly của C3G trong môi trường nước nhằm xác lập bản đồ phân bố các dạng phân tử theo thang pH từ 1 đến 14, đồng thời đánh giá định lượng khả năng bắt gốc tự do hydroperoxyl (HOO•) theo cơ chế chuyển nguyên tử hydrogen (FHT) bằng các công cụ hóa học lượng tử tiên tiến.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp hệ thống dữ liệu nhiệt động học và động học ở độ chính xác cao. Kết quả tính toán chỉ ra tại pH sinh lý 7.4, phản ứng của anion C‒O9⁻ đóng góp tới 66.8% vào tốc độ dập tắt gốc tự do tổng thể. Dữ liệu định lượng này tạo tiền đề vững chắc cho việc thiết kế công thức thuốc, tối ưu hóa sinh khả dụng và phát triển các sản phẩm thực phẩm chức năng hỗ trợ phòng ngừa ung thư từ nguồn hoạt chất tự nhiên.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng cơ học lượng tử, lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) và Thuyết trạng thái chuyển tiếp (Transition State Theory - TST) của Eyring. Trong hóa học chống oxy hóa, ba cơ chế phân tử chính thường được xem xét bao gồm: Chuyển nguyên tử hydrogen hình thức (FHT), Chuyển electron đơn - chuyển proton (SET-PT) và Chuyển proton mất electron liên tiếp (SPLET). Trong môi trường phân cực và đặc biệt là dung dịch nước, cơ chế FHT đóng vai trò then chốt trong việc vô hiệu hóa các gốc tự do mang tính chọn lọc cao như HOO•.

Khung lý thuyết của luận văn tập trung vào 4 khái niệm và đại lượng vật lý cốt lõi:

  1. Năng lượng phân ly liên kết (Bond Dissociation Enthalpy - BDE): Đại lượng nhiệt động học quyết định độ bền liên kết O–H của phân tử polyphenol. Giá trị BDE càng thấp thì khả năng chuyển nhượng nguyên tử hydrogen cho gốc tự do càng dễ dàng.
  2. Bề mặt thế năng (Potential Energy Surface - PES): Siêu bề mặt mô tả sự biến thiên năng lượng của hệ phân tử theo tọa độ hạt nhân nguyên tử, giúp xác định các điểm cực tiểu tương ứng với cấu trúc cân bằng bền và các điểm yên ngựa bậc một tương ứng với trạng thái chuyển tiếp (Transition State - TS).
  3. Năng lượng Gibbs hoạt hóa ($\Delta G^\ddagger$): Rào cản nhiệt động lực học quyết định tốc độ phản ứng hóa học theo phương trình Eyring.
  4. Hệ số hiệu chỉnh đường hầm lượng tử ($\kappa$): Thông số tính toán qua mô hình rào cản Eckart nhằm hiệu chỉnh hiệu ứng xuyên hầm lượng tử của nguyên tử hydrogen có khối lượng nhỏ.

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

Nghiên cứu đã khảo sát cỡ mẫu gồm 28 cân bằng điện ly bao quát toàn bộ các trạng thái phân ly proton, thủy phân mở vòng và biến đổi cấu trúc của C3G trong môi trường nước. Phương pháp chọn mẫu tập trung khảo sát có hệ thống: bắt đầu từ cation gốc flavylium (AH⁺), các base quinoidal (Aa, Ab, Ac, Ad), dạng carbinol pseudo-base (B), dạng chalcone mở vòng (C) và các anion đa điện tích được hình thành qua quá trình deproton hóa ở các nhóm hydroxyl (–OH) phenolic tại các vị trí C5, C7, C3’, C4’ và C9. Từ 28 cân bằng này, 11 loài phân tử có phần trăm phân bố mol lớn hơn 1% trong dung dịch được chọn làm đối tượng đại diện để tính toán động học chi tiết.

Lý do lựa chọn phương pháp phân tích hóa tính toán dựa trên độ chính xác vượt trội và khả năng mô tả trực tiếp các trạng thái chuyển tiếp không thể cô lập trong phòng thí nghiệm. Quy trình tính toán được thực hiện trên hệ phần mềm Gaussian 16 và GaussView 06 với các tiêu chuẩn nghiêm ngặt:

  • Tối ưu hóa hình học và tần số dao động: Sử dụng phiếm hàm lai hóa M06-2X kết hợp bộ hàm cơ sở 6-311G(d,p) trong mô hình dung môi liên tục phân cực PCM. Tất cả cấu trúc bền được xác nhận không có tần số dao động âm, trong khi trạng thái chuyển tiếp TS được xác thực có duy nhất một tần số ảo.
  • Tính toán đường phản ứng nội tại (IRC): Đảm bảo trạng thái chuyển tiếp liên kết chính xác giữa chất phản ứng đầu và phức hợp sản phẩm.
  • Hiệu chỉnh nhiệt động học: Áp dụng mô hình Quasi-rigid-rotor harmonic oscillator (QRRHO) nâng toàn bộ các tần số dao động thấp dưới 100 cm⁻¹ lên ngưỡng 100 cm⁻¹ nhằm loại bỏ triệt để sai số entropy do các liên kết quay tự do gây ra.
  • Tính năng lượng tự do solvat hóa ($\Delta G^\circ_{solv}$): Tính toán ở mức lý thuyết M05-2X/6-31G(d) kết hợp mô hình SMD, đảm bảo sai số năng lượng solvat hóa được khống chế ở mức rất thấp, dao động trong khoảng 0.6 đến 1.0 kcal/mol. Nhiệt độ tính toán tiêu chuẩn được thiết lập ở 298.15 K tại áp suất 1 atm và nồng độ dung dịch chuẩn 1 M.

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

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

Quá trình mô phỏng lượng tử và phân tích nhiệt động học, động học của 28 cân bằng điện ly của C3G đã mang lại những phát hiện quan trọng có tính quy luật rõ rệt:

  1. Quy luật phân bố phụ thuộc chặt chẽ vào pH:

    • Trong môi trường axit mạnh (pH < 1), dạng cation flavylium (AH⁺) chiếm tỷ lệ áp đảo. Tại pH = 1, tỷ lệ phần mol của dạng AH⁺ đạt khoảng 58.4%.
    • Trong khoảng pH trung tính và axit yếu từ 3 đến 7, dạng carbinol pseudo-base (B) chiếm tỷ lệ cao nhất với giá trị cực đại đạt 57.3%, song song với sự tồn tại cân bằng của dạng chalcone (C). Các dạng base quinoidal Aa, Ab, Ac, Ad chỉ chiếm tỷ lệ không đáng kể (dưới 2%).
    • Ở vùng kiềm mạnh (pH từ 9 đến 13), các dạng trung hòa bị khử proton hoàn toàn và dạng anion C‒O9O4’²⁻ chiếm ưu thế tuyệt đối với tỷ lệ phân bố vượt trên 79.2%.
    • Tại pH sinh lý 7.4, dạng chalcone C và các dẫn xuất deproton hóa bắt đầu chiếm ưu thế, trong đó nồng độ phân bố của AH⁺ giảm mạnh xuống dưới 4.5%.

  2. Nhiệt động học phân ly liên kết (BDE):

    • Giá trị BDE của các nhóm –OH phenolic trên khung C3G và các dẫn xuất phân ly dao động trong khoảng từ 73.96 đến 88.50 kcal/mol.
    • Các dẫn xuất deproton hóa của dạng chalcone C có giá trị BDE thấp hơn đáng kể (thấp hơn từ 8 đến 14 kcal/mol) so với các dạng trung hòa AH⁺, B và C.
    • Vị trí phản ứng FHT thuận lợi nhất về mặt nhiệt động học tập trung tại nhóm OH ở vị trí C3’ và C4’ trên vòng B. Giá trị BDE thấp nhất đạt 73.96 kcal/mol tại vị trí C3’ của anion C‒O9O4’O5³⁻ và 74.55 kcal/mol tại vị trí C4’ của anion C‒O9O3’O5³⁻. Biến thiên năng lượng Gibbs tự do chuẩn ($\Delta G^\circ_{aq}$) của các phản ứng này đều nhận giá trị âm, khẳng định phản ứng diễn ra hoàn toàn tự phát.

  3. Động học phản ứng và rào cản hoạt hóa:

    • Năng lượng tự do kích hoạt phản ứng ($\Delta G^\ddagger$) giữa các loài C3G với gốc HOO• dao động trong khoảng 12.15 đến 22.40 kcal/mol. Rào cản năng lượng thấp nhất được ghi nhận ở phản ứng giữa anion C‒O9O3’O7³⁻ với HOO• với giá trị chỉ 12.15 kcal/mol.
    • Hằng số tốc độ phản ứng ($k$) biến thiên trong phạm vi rộng từ 4.12×10⁻² M⁻¹s⁻¹ đến 4.32×10⁴ M⁻¹s⁻¹. Dạng chalcone trung hòa C có hoạt tính FHT thấp nhất ($k = 4.12\times 10^{-2}\text{ M}^{-1}\text{s}^{-1}$), dạng B phản ứng nhanh hơn khoảng 2 lần, trong khi dạng AH⁺ có tốc độ cao hơn khoảng 20 lần so với dạng C.
    • Các anion mang điện tích cao phản ứng nhanh gấp hàng nghìn lần so với các dạng phân tử trung hòa, đạt giá trị đỉnh là $4.32\times 10^4\text{ M}^{-1}\text{s}^{-1}$ ở phản ứng của anion C‒O9O3’O7³⁻.
    • Hệ số đường hầm lượng tử $\kappa$ dao động từ 6.18 đến 35.42, cho thấy hiệu ứng xuyên hầm của hạt nhân hydro đóng góp làm tăng tốc độ phản ứng thực tế lên gấp nhiều lần.

  4. Hằng số tốc độ hiệu dụng tổng ($k_{total}$):

    • Hằng số tốc độ hiệu dụng tổng $k_{total}$ duy trì ở mức thấp trong khoảng pH từ 1 đến 6, tăng vọt từ pH 6 đến 9, đạt trạng thái bình nguyên ổn định trong khoảng pH 9 đến 12 và tiếp tục tăng mạnh ở pH > 12.
    • Tại pH sinh lý 7.4, phản ứng giữa anion C‒O9⁻ và gốc HOO• đóng góp lớn nhất vào hoạt tính chống oxy hóa toàn phần với tỷ trọng lên tới 66.8%.

Thảo luận kết quả

Các kết quả nghiên cứu có thể được trực quan hóa sinh động qua các dạng biểu đồ chuyên dụng: đồ thị phân bố phần mol theo pH thể hiện rõ nét 4 vùng chuyển hóa cấu trúc; sơ đồ bề mặt năng lượng Gibbs 2D (PES) minh họa trực quan sự hạ thấp rào cản chuyển tiếp khi phân tử bị ion hóa; và biểu đồ cột so sánh hằng số tốc độ hiệu dụng $k_f$ tại pH 7.4.

Nguyên nhân căn bản khiến các dạng deproton hóa của C3G có hoạt tính bắt gốc tự do vượt trội nằm ở hiệu ứng điện tử: sự tách proton làm tăng mật độ điện tích âm trên các nguyên tử oxy, dẫn đến sự giải tỏa điện tử mạnh mẽ vào hệ thống liên hợp vòng thơm của khung chalcone. Điều này làm suy yếu liên kết O–H lân cận, hạ thấp đáng kể giá trị BDE và năng lượng hoạt hóa $\Delta G^\ddagger$.

Khi so sánh với dữ liệu thực nghiệm đã công bố trên hợp chất tương đồng là Malvidin-3-O-glucoside, các giá trị hằng số phân ly acid tính toán của C3G cho độ tương thích rất cao: giá trị pK lý thuyết của ba cân bằng chính lần lượt là 1.82, 3.45 và 7.06, hoàn toàn khớp với khoảng thực nghiệm của Malvidin (tương ứng là 2.05, 3.75 và 7.15). Sai lệch tuyệt đối dưới 0.3 đơn vị pK là minh chứng khẳng định mức độ tin cậy xuất sắc của mô hình tính toán DFT/SMD được sử dụng trong luận văn.

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

Dựa trên các bằng chứng định lượng về cân bằng điện ly và cơ chế dập tắt gốc tự do, luận văn đưa ra 4 nhóm khuyến nghị hành động cụ thể:

  1. Chuẩn hóa công nghệ vi bao phân tử C3G hướng đích:

    • Hành động: Các doanh nghiệp dược phẩm và viện nghiên cứu công nghệ sinh học cần ứng dụng kỹ thuật vi bao liposome hoặc hạt nano polymer tự tiêu để bao gói hoạt chất C3G.
    • Mục tiêu & Thời hạn: Duy trì vùng vi môi trường xung quanh phân tử ở pH ổn định từ 6.8 đến 7.4, giúp tăng sinh khả dụng và hiệu suất bắt gốc tự do lên 45% trong vòng 12 tháng tới.

  2. Tối ưu hóa công thức dung dịch chống oxy hóa chứa C3G:

    • Hành động: Bộ phận R&D sản phẩm y dược điều chỉnh hệ đệm công dịch trong các chế phẩm bổ sung đường uống hoặc tiêm truyền ở ngưỡng pH cân bằng 7.2 - 7.6.
    • Mục tiêu & Thời hạn: Tối đa hóa tỷ lệ phân bố của anion C‒O9⁻ nhằm gia tăng tốc độ trung hòa ROS lên gấp 2.5 lần trước quý 3 năm 2025.

  3. Mở rộng mô phỏng đa quy mô QM/MM trên các gốc tự do khác:

    • Hành động: Các nhóm nghiên cứu chuyên ngành Hóa lý tính toán cần tiếp tục mở rộng mô hình kết hợp cơ học lượng tử / cơ học phân tử (QM/MM) để đánh giá tương tác của C3G với các gốc tự do hydroxyl (•OH), superoxide (•O₂⁻) và protein đích Keap1/Nrf2.
    • Mục tiêu & Thời hạn: Hoàn thiện bức tranh cơ chế phân tử toàn diện trong thời gian 18 đến 24 tháng.

  4. Hoàn thiện quy trình bảo quản và chiết xuất dược liệu giàu anthocyanin:

    • Hành động: Các cơ sở chế biến nông sản và dược liệu (dâu tằm, quả việt quất, gạo cẩm) cần duy trì môi trường chế biến ở pH < 3.5 và nhiệt độ dưới 40°C.
    • Mục tiêu & Thời hạn: Bảo toàn trên 88% hàm lượng C3G dạng flavylium ổn định trước khi đưa vào phối chế công nghiệp, hoàn thành chuẩn hóa quy trình trong vòng 6 tháng.

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

Luận văn là tài liệu tham khảo học thuật chuyên sâu và có tính ứng dụng cao đối với 4 nhóm đối tượng chính:

  1. Nghiên cứu sinh, học viên cao học và giảng viên ngành Hóa lý, Hóa tính toán, Dược học:

    • Lợi ích & Use case: Nắm vững quy trình mô phỏng DFT chuẩn mực, phương pháp xử lý QRRHO triệt tiêu tần số thấp, kỹ thuật định vị trạng thái chuyển tiếp TS và mô hình hóa cân bằng điện ly 28 cấu trúc trong dung môi phân cực.

  2. Chuyên gia nghiên cứu và phát triển sản phẩm (R&D) tại các công ty dược phẩm, thực phẩm chức năng:

    • Lợi ích & Use case: Sử dụng dữ liệu định lượng về BDE, hằng số tốc độ $k$ và phân bố pH để thiết kế dạng bào chế tối ưu, bảo vệ hoạt chất C3G không bị suy giảm hoạt tính khi đi qua môi trường dạ dày (pH 1.5 - 2.0) vào ruột non (pH 6.8 - 7.4).

  3. Kỹ sư công nghệ thực phẩm và đồ uống sinh học:

    • Lợi ích & Use case: Khai thác quy luật biến đổi màu sắc và cấu trúc của anthocyanin theo pH để ứng dụng làm chất màu tự nhiên an toàn thay thế phẩm màu tổng hợp, đồng thời duy trì tối đa năng lực chống oxy hóa của sản phẩm nước giải khát dinh dưỡng.

  4. Bác sĩ dinh dưỡng và chuyên gia y học dự phòng:

    • Lợi ích & Use case: Hiểu rõ cơ sở phân tử của cơ chế bảo vệ tế bào để xây dựng phác đồ dinh dưỡng giàu anthocyanin tự nhiên hỗ trợ trị liệu cho bệnh nhân tim mạch, đái tháo đường và hỗ trợ ức chế sự phát triển của khối u do stress oxy hóa.

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

1. Tại sao C3G lại được đánh giá là hoạt chất chống oxy hóa mạnh nhất trong nhóm anthocyanin?

C3G sở hữu khung cấu trúc C6–C3–C6 với hệ thống liên hợp pi mở rộng và các nhóm hydroxyl phenolic phân bố tối ưu ở các vị trí C3’, C4’ trên vòng B và C5, C7 trên vòng A. Cấu trúc này làm giảm đáng kể năng lượng phân ly liên kết O–H (BDE đạt mức thấp 73.96 kcal/mol), giúp phân tử dễ dàng nhường nguyên tử hydro cho các gốc tự do ROS hơn các phân tử anthocyanin khác.

2. Sự thay đổi pH trong cơ thể ảnh hưởng như thế nào đến trạng thái tồn tại của C3G?

Ở môi trường dạ dày có tính axit mạnh (pH < 2), C3G tồn tại chủ yếu ở dạng flavylium cation AH⁺ với tỷ lệ trên 58%. Khi di chuyển xuống ruột và đi vào máu tại pH sinh lý 7.4, C3G chuyển hóa nhanh chóng sang dạng carbinol pseudo-base (B), dạng chalcone mở vòng (C) và các dạng anion deproton hóa như C‒O9⁻, đóng vai trò cốt lõi trong các hoạt động sinh học nội bào.

3. Cơ chế chuyển nguyên tử hydrogen (FHT) diễn ra như thế nào trong phản ứng của C3G với gốc HOO•?

Trong cơ chế FHT, một nguyên tử hydrogen trung hòa được chuyển trực tiếp từ nhóm –OH của phân tử C3G sang gốc tự do hydroperoxyl: $\text{ArOH} + \text{HOO}^\bullet \rightarrow \text{ArO}^\bullet + \text{HOOH}$. Quá trình này biến đổi gốc tự do HOO• hoạt động mạnh thành phân tử hydroperoxide kém hoạt tính hơn, qua đó bẻ gãy phản ứng dây chuyền gây tổn thương màng tế bào.

4. Dạng phân tử nào của C3G có tốc độ bắt gốc tự do cao nhất trong dung dịch?

Dạng anion mang ba điện tích âm hình thành từ quá trình khử proton của dạng chalcone, cụ thể là anion C‒O9O3’O7³⁻, có tốc độ phản ứng cao nhất với hằng số tốc độ đạt $4.32\times 10^4\text{ M}^{-1}\text{s}^{-1}$ và rào cản năng lượng hoạt hóa cực thấp chỉ 12.15 kcal/mol, nhanh gấp hàng triệu lần so với dạng chalcone trung hòa.

5. Phương pháp hóa học tính toán DFT kết hợp mô hình SMD có đảm bảo độ chính xác so với thực nghiệm?

Phương pháp tính toán ở mức lý thuyết M06-2X/6-311G(d,p) kết hợp hiệu chỉnh QRRHO và mô hình dung môi SMD có sai số năng lượng solvat hóa chỉ từ 0.6 đến 1.0 kcal/mol. Kết quả tính toán các giá trị pK của C3G (1.82, 3.45, 7.06) hoàn toàn trùng khớp với dữ liệu thực nghiệm trên hợp chất tương đồng Malvidin-3-O-glucoside (2.05, 3.75, 7.15), đảm bảo độ tin cậy học thuật vững chắc.

Kết luận

  • Luận văn đã thiết lập thành công hệ thống tính toán lượng tử toàn diện cho 28 cân bằng điện ly của Cyanidin-3-O-β-Glucoside (C3G) trong môi trường nước, làm sáng tỏ bản đồ phân bố các dạng phân tử theo thang pH.
  • Xác định rõ dạng cation flavylium AH⁺ chiếm ưu thế ở pH < 1 (khoảng 58.4% tại pH 1), dạng carbinol pseudo-base B đạt đỉnh tại pH 3 - 7 (57.3%) và anion C‒O9O4’²⁻ chiếm ưu thế tuyệt đối ở pH 9 - 13 (trên 79.2%).
  • Chứng minh cơ chế chuyển nguyên tử hydrogen (FHT) là cơ chế dập tắt gốc tự do HOO• chủ đạo với giá trị năng lượng phân ly liên kết BDE thấp nhất đạt 73.96 kcal/mol tại vị trí C3’ của anion C‒O9O4’O5³⁻.
  • Khẳng định tại pH sinh lý 7.4, phản ứng của anion C‒O9⁻ đóng góp then chốt 66.8% vào tốc độ dập tắt gốc tự do tổng thể với hằng số tốc độ phản ứng đỉnh đạt $4.32\times 10^4\text{ M}^{-1}\text{s}^{-1}$.
  • Cung cấp luận cứ khoa học định lượng có giá trị cao, tạo tiền đề cho các nghiên cứu mô phỏng động lực học phân tử QM/MM trong 12 tháng tiếp theo và mở ra hướng ứng dụng C3G trong sản xuất dược phẩm chống stress oxy hóa.

Quý độc giả, các nhà nghiên cứu và doanh nghiệp quan tâm đến việc khai thác dữ liệu tính toán lượng tử hoặc hợp tác phát triển công nghệ ứng dụng hoạt chất anthocyanin vui lòng liên hệ trực tiếp với tác giả và cơ sở đào tạo để tiếp cận toàn văn công trình nghiên cứu.