Tổng quan về luận án
Nghiên cứu của tiến sĩ Nguyễn Minh Quang với đề tài "Thiết kế, sàng lọc và tổng hợp một số dẫn xuất thiosemicarbazone và phức chất dựa trên các tính toán hóa lượng tử kết hợp phương pháp mô hình hóa QSPR" thuộc chuyên ngành Hóa lý thuyết và Hóa lý (Mã số: 944.19) tại Trường Đại học Khoa học – Đại học Huế phối hợp cùng Trường Đại học Công nghiệp Thành phố Hồ Chí Minh dưới sự hướng dẫn khoa học của PGS. Phạm Văn Tất và TS. Trần Xuân Mậu, là công trình tiên phong xác lập cầu nối tích hợp giữa hóa học tính toán đa tầng và thực nghiệm phân tích - tổng hợp vô cơ hiện đại.
Bối cảnh khoa học của công trình khởi phát từ nhu cầu cấp bách trong việc quan trắc môi trường và phát triển dược chất. Các ion kim loại chuyển tiếp nặng như $\text{Cd}^{2+}$, $\text{Ni}^{2+}$, $\text{Cu}^{2+}$, $\text{Zn}^{2+}$, $\text{Hg}^{2+}$, $\text{Pb}^{2+}$ là tác nhân tích lũy sinh học nguy hiểm gây ung thư và dị tật bẩm sinh, đòi hỏi các phương pháp định lượng trắc quang có độ nhạy cao. Dẫn xuất thiosemicarbazone (TSC) mang hệ dị nguyên tử nitơ - lưu huỳnh ($N, S$) sở hữu khả năng phối trí linh động và hoạt tính kháng sinh, kháng nấm vượt trội. Tuy nhiên, một research gap then chốt trong y văn quốc tế đã được tác giả xác định trực tiếp: "chưa có công bố nào phát triển mô hình QSPR trên đối tượng phức chất giữa thiosemicarbazone và các ion kim loại với đại lượng đặc trưng là hằng số bền của phức". Phần lớn nghiên cứu trước đây (như công trình của Domagk [39], Reddy [113], hay Singh [118]) thuần túy tiến hành thử nghiệm - sai (trial-and-error) hoặc khảo sát lý thuyết đơn lập mà thiếu hẳn mô hình dự báo định lượng khả năng tạo phức trước khi tổng hợp.
Để giải quyết khoảng trống này, luận án thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:
- RQ1 / H1: Liệu các tham số cấu trúc 2D/3D và chỉ số hóa lượng tử bán thực nghiệm (PM7) có thể giải thích định lượng độ biến thiên hằng số bền ($\log\beta_{11}, \log\beta_{12}$) của phức kim loại - thiosemicarbazone hay không?
- RQ2 / H2: Sự kết hợp giữa giải thuật di truyền (GA) với các kỹ thuật học máy hồi quy tuyến tính bội (MLR), máy học vectơ hỗ trợ (SVR) và mạng thần kinh nhân tạo (ANN) có tạo ra được mô hình dự báo thỏa mãn nghiêm ngặt 5 nguyên tắc OECD hay không?
- RQ3 / H3: Tính toán bề mặt thế năng (PES) và tối ưu hóa cấu hình lập thể có dự đoán chính xác cấu dạng bền nhất và cơ chế phối trí hai càng (bidentate) qua nguyên tử lưu huỳnh thione/thiolate và nitơ azomethine trước khi chuyển sang pha thực nghiệm hay không?
- RQ4 / H4: Các dẫn xuất phối tử mới thiết kế mang khung dị vòng phenothiazine (BEPT) và carbazole (BECT) có thể tổng hợp thành công và hình thành phức chất màu bền với hằng số bền thực nghiệm khớp với giá trị dự báo từ mô hình QSPR trong sai số cho phép ($MSE < 0,05$)?
Khung lý thuyết của nghiên cứu tích hợp lý thuyết quan hệ định lượng cấu trúc - tính chất (QSPR) khởi nguồn từ nguyên lý năng lượng tự do tuyến tính (LFER) của Hammett - Taft, mô hình đa thông số Hansch, lý thuyết đồ thị hóa học Wiener - Kier - Hall, cơ học lượng tử bán thực nghiệm Stewart (PM7) và phương trình quang phổ Benesi-Hildebrand. Luận án khai thác tập mẫu gồm 54 phối tử thiosemicarbazone thực nghiệm, xây dựng bộ 230 mô tả phân tử 2D/3D, phân loại và sàng lọc dữ liệu phức chất đơn nhân phân tử bậc 1 vị trí phối trí ($ML$) và bậc 2 ($ML_2$). Từ đó, tác giả thiết kế 44 phối tử thiosemicarbazone mới, dự báo hằng số bền cho 176 phức chất kim loại giả định, lựa chọn và tổng hợp thành công 2 phối tử mới gồm BEPT (2-((7-bromo-10-ethyl-10H-phenothiazin-3-yl)methylene)hydrazine-1-carbothioamide) và BECT (2-((6-bromo-9-ethyl-9H-carbazol-3-yl)methylene)hydrazine-1-carbothioamide), cùng 4 phức chất mới $\text{Cd(II)-BEPT}$, $\text{Ni(II)-BEPT}$, $\text{Cu(II)-BECT}$, $\text{Zn(II)-BECT}$. Kết quả kiểm định cho thấy mô hình GA-MLR đạt $R^2_{train} = 0,9446$ và hệ số đánh giá chéo loại trừ từng cá thể $Q^2_{LOO} = 0,9262$, mang lại giải pháp định lượng đột phá giúp rút ngắn hơn 70% thời gian sàng lọc và chi phí hóa chất thực nghiệm.
Literature Review và Positioning
Nghiên cứu về mối quan hệ giữa cấu trúc và tính chất/hoạt tính sinh học (QSAR/QSPR) đã trải qua hơn 150 năm tiến hóa lý thuyết với các cột mốc nền tảng:
- Crum Brown và Fraser (1868) [14] lần đầu tiên thiết lập mối tương quan toán học giữa cấu trúc hóa học và đáp ứng sinh học trên các dẫn xuất strychnine.
- Richet (1893) [86] và Meyer - Overton (1899) chứng minh độc tính và hoạt tính gây mê tương quan trực tiếp với hệ số phân bố dầu-nước ($\log P$), đặt nền móng cho tham số kỵ nước.
- Hammett (1935) và Taft (1952) [86] phát triển mô hình LFER, chuẩn hóa hằng số nhóm thế $\sigma$, hiệu ứng không gian lập thể ($E_s$) và hiệu ứng cộng hưởng.
- Hansch và Fujita (1962) [86] hợp nhất các tham số điện tử, lập thể và kỵ nước trong mô hình hồi quy đa biến đa thông số.
- Wiener (1947), Platt (1947), Randić (1975), Kier và Hall (1977) [86] cách mạng hóa lĩnh vực bằng lý thuyết đồ thị hóa học, giới thiệu các chỉ số kết nối phân tử (molecular connectivity indices) và trạng thái điện - cấu trúc topo (electrotopological state).
Trong lĩnh vực hóa học phối trí của thiosemicarbazone, Domagk và cộng sự (1946) [39] mở đầu cho làn sóng tổng hợp phức chất kim loại chuyển tiếp nhằm tăng cường hoạt tính kháng khuẩn. Tiếp đó, hàng loạt nghiên cứu quốc tế đã khai thác tính chất tạo phức của dẫn xuất TSC:
- Reddy và cộng sự (2013) [113] ứng dụng dẫn xuất 5-bromosalicylaldehyde thiosemicarbazone để định lượng quang phổ Cu(II).
- Wang và cộng sự (2014) [156] phát triển mô hình 3D-QSAR (CoMFA/CoMSIA) trên các dẫn xuất benzaldehyde thiosemicarbazone ức chế enzyme phenoloxidase.
- Singh và cộng sự (2011, 2013) [118], [120] sử dụng 2-acetylpyridine 4-methyl-3-thiosemicarbazone (APMT) và 3-APT làm thuốc thử trắc quang phân tích vi lượng Cu(II) và Ni(II).
- Al-Amiery và cộng sự (2012) [138] kết hợp tính toán DFT/B3LYP và AIM để giải thích cấu trúc phân tử của dẫn xuất thiosemicarbazone dị vòng pyrrole.
Tuy nhiên, tồn tại hai luồng quan điểm đối nghịch trong y văn: Một trường phái cho rằng các mô hình hồi quy tuyến tính cổ điển (MLR) là đủ để biểu diễn mối tương quan cấu trúc - hằng số bền dựa trên các hiệu ứng điện tử định xứ; trường phái thứ hai khẳng định tương tác phối trí kim loại - phối tử chịu chi phối bởi hiệu ứng bất đối xứng trường tinh thể và tương tác quỹ đạo phi tuyến, bắt buộc phải dùng các thuật toán hộp đen như mạng thần kinh nhân tạo (ANN) hay máy học vectơ hỗ trợ (SVM).
Luận án của NCS. Nguyễn Minh Quang định vị chính xác ở điểm giao thoa này bằng cách triển khai phân tích so sánh song song giữa các mô hình tuyến tính (MLR, PLSR, PCR) và phi tuyến (ANN, SVR) kết hợp tối ưu hóa biến bằng giải thuật di truyền (GA). Điểm đột phá vượt trội so với các công trình quốc tế nằm ở chỗ: thay vì chỉ dừng lại ở mô hình hóa in silico (như Wang [156]) hoặc tổng hợp thực nghiệm đơn thuần (như Reddy [113], Singh [118]), luận án xây dựng một chu trình khép kín: Thu thập dữ liệu thực nghiệm $\rightarrow$ Mô hình hóa QSPR $\rightarrow$ Sàng lọc & Thiết kế phối tử ảo mới $\rightarrow$ Tính toán cấu dạng năng lượng lượng tử $\rightarrow$ Tổng hợp thực nghiệm $\rightarrow$ Kiểm định quang phổ & Tái thẩm định mô hình.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đóng góp trực tiếp vào việc mở rộng thuyết liên kết phối trí hiện đại và lý thuyết năng lượng tự do tuyến tính (LFER):
- Mở rộng lý thuyết tương quan LFER của Hammett-Taft cho hệ phức chất kim loại: Chứng minh rằng hằng số bền $\log\beta_{pq}$ không chỉ phụ thuộc vào mật độ điện tích tại tâm phối trí mà còn bị chi phối mạnh mẽ bởi cấu trúc topo không gian và độ bất đối xứng phân tử.
- Làm sáng tỏ cân bằng hỗ biến Thione - Thiol trong pha dung dịch: Luận án trích dẫn cơ chế: "Dạng thione hoạt động như một phối tử hai càng (hai phối trí) trung hòa, trong khi sự mất đi proton ở vị trí thio (lưu huỳnh) từ thiol tạo ra phối tử hai càng mang điện tích". Tính toán hóa lượng tử bán thực nghiệm PM7 đã cung cấp bằng chứng định lượng về hàng rào năng lượng hỗ biến và thế năng bề mặt (PES), chứng minh dạng thiolate phối trí ưu tiên qua hệ liên hợp liên phân tử đối với các ion kim loại chuyển tiếp hóa trị hai ($\text{Ni}^{2+}, \text{Cd}^{2+}, \text{Cu}^{2+}, \text{Zn}^{2+}$).
- Xác lập hệ biến số mô tả cốt lõi (Core Descriptors): Bằng thuật toán GA-MLR với 44 tổ hợp huấn luyện, luận án chứng minh sự đóng góp phần trăm trung bình (APC%) của các tham số: chỉ số liên kết bậc cao $xp5$, $xp3$, số vòng thơm ($nrings$) và độ elip phân tử ($ovality$) là những biến số quyết định độ bền nhiệt động học của phức $ML$ và $ML_2$.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp liên ngành 3 trục lý thuyết trụ cột:
- Trục Hóa học lượng tử (Quantum Chemistry): Khảo sát trường lực phân tử MM+ kết hợp cơ học lượng tử bán thực nghiệm PM7 trên phần mềm MOPAC2016 (phiên bản Version 17 do GS. James Stewart cung cấp) để tối ưu hóa hình học, tính toán năng lượng orbital phân tử HOMO/LUMO, nhiệt tạo thành ($\Delta H_f$), năng lượng liên kết nhân-nhân, diện tích và thể tích bề mặt phân cực COSMO.
- Trục Hóa trắc học & Học máy (Chemometrics & Machine Learning): Áp dụng thuật toán phân cụm k-means và phân cụm phân cấp tập hợp (AHC) để phân chia khách quan tập dữ liệu thành tập huấn luyện (training set) và tập kiểm định ngoại (external validation set), kết hợp GA-MLR, GA-SVR (với hàm nhân RBF tối ưu hóa 5 cấp qua lưới tham số $\gamma$ và $C$), cùng kiến trúc mạng nơ-ron truyền thẳng lan truyền ngược sai số QSPR-ANN dạng $I(11)-HL(n)-O(1)$ và $I(6)-HL(6)-O(1)$.
- Trục Thực nghiệm Phổ nghiệm & Trắc quang (Spectroscopy & Photometry): Khung phân tích cân bằng dung dịch thông qua định luật Lambert-Beer và phương trình Benesi-Hildebrand tuyến tính hóa:
$$\frac{[\text{M}]0}{A} = \frac{1}{\varepsilon \cdot l \cdot \beta{1n} \cdot [\text{L}]_0^n} + \frac{1}{\varepsilon \cdot l}$$
Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ: áp dụng cho các phức chất kim loại chuyển tiếp phân nhóm $d^{8}, d^{9}, d^{10}$ với phối tử thiosemicarbazone mang khung thơm dị vòng đa nhân trong hệ dung môi phân cực (DMF/nước), lực ion $I$ được kiểm soát bởi dung dịch đệm và nhiệt độ phản ứng tiêu chuẩn $298 - 303\text{ K}$.
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ủ triết lý nghiên cứu thực chứng nghiêm ngặt (positivism paradigm), vận dụng thiết kế thực nghiệm kết hợp tính toán đa cấp (multi-level computational-experimental design).
flowchart TD
A[Bộ dữ liệu 54 Ligand & Phức chất thực nghiệm] --> B[Tối ưu hóa cấu trúc 3D MM+ & PM7 trên MOPAC2016]
B --> C[Trích xuất 230 Bộ mô tả phân tử 2D/3D/Lượng tử]
C --> D[Phân cụm dữ liệu k-means & AHC: Training & Test Sets]
D --> E[Lựa chọn biến bằng GA: GA-MLR, GA-SVR, GA-ANN]
E --> F[Đánh giá chéo LOO & Miền ứng dụng Williams Plot / Grubb Test]
F --> G[Sàng lọc ảo 44 TSC mới & Dự báo log Beta cho 176 phức chất]
G --> H[Khảo sát cấu dạng bền PES & Năng lượng quay góc nhị diện]
H --> I[Tổng hợp hóa học đa giai đoạn: BEPT & BECT]
I --> J[Tổng hợp phức Ni-BEPT, Cd-BEPT, Cu-BECT, Zn-BECT]
J --> K[Xác định cấu trúc: FT-IR, 1D/2D-NMR, HR-MS, EDX-SEM]
K --> L[Xác định thành phần Job/Molar Ratio & Hằng số bền Benesi-Hildebrand]
L --> M[So sánh Thực nghiệm vs Dự báo: Đánh giá độ chuẩn xác mô hình]
Cấu trúc mẫu nghiên cứu được chuẩn hóa từ 54 phối tử ban đầu, phân tách thành các tập dữ liệu phức $ML$ và $ML_2$ qua 9 nhóm dữ liệu phân tầng nhằm đảm bảo tính đại diện thống kê và tránh hiện tượng quá khớp (overfitting).
Quy trình nghiên cứu rigorous
- Sàng lọc cấu trúc và tiền xử lý dữ liệu: Cấu trúc không gian của các phối tử và phức chất được vẽ và tối ưu hóa năng lượng ban đầu bằng trường lực cơ học phân tử MM+ trên phần mềm HyperChem, sau đó tối ưu hóa toàn phần ở mức bán thực nghiệm PM7 trên phần mềm MOPAC2016. Các cấu trúc dị thường hình học hoặc không hội tụ gradient năng lượng bị loại bỏ sau kiểm định ngoại lệ Grubb ở mức tin cậy 95%.
- Trích xuất và chuẩn hóa biến: 230 bộ mô tả phân tử thuộc các nhóm 2D (chỉ số topo, liên kết, điện tích nguyên tử riêng phần) và 3D (diện tích bề mặt phân tử, thể tích tương tác, momen lưỡng cực, ovality, HOMO, LUMO) được tính toán. Loại bỏ các biến bất biến (zero variance) và các biến có hệ số tương quan chéo liên biến $r > 0,90$.
- Chiến lược lấy mẫu và kiểm định: Sử dụng thuật toán phân cụm phân cấp tập hợp (AHC) với khoảng cách bình phương Euclide để chia tập mẫu thành Training Set (huấn luyện mô hình) và Test Set (kiểm định ngoại) theo tỷ lệ tối ưu 4:1.
- Quy trình tam giác hóa phương pháp (Triangulation):
- Tam giác hóa phương pháp mô hình hóa: So sánh độc lập kết quả giữa MLR, PLSR, SVR và ANN.
- Tam giác hóa phổ nghiệm cấu trúc: Kết hợp đồng thời phổ hồng ngoại biến đổi Fourier (FT-IR), phổ cộng hưởng từ hạt nhân 1 chiều ($^1\text{H}$-NMR, $^{13}\text{C}$-NMR, DEPT-90, DEPT-135, CPD), phổ 2 chiều ($^1\text{H}-^{13}\text{C}$ HSQC, $^{1}\text{H}-^{13}\text{C}$ HMBC), phổ khối lượng phân giải cao (HR-MS với nguồn ion hóa phun mù điện tử ESI), phổ phân tán năng lượng tia X (EDX) và kính hiển vi điện tử quét (SEM).
- Tam giác hóa xác định công thức phức: Kết hợp phương pháp biến đổi liên tục Job và phương pháp tỷ lệ mol (molar ratio method) trên hệ phổ hấp thụ tử ngoại - khả kiến (UV-Vis).
Data và phân tích
Quá trình phân tích dữ liệu lượng tử và học máy được xử lý trên các phần mềm chuyên dụng: HyperChem 8.0, MOPAC2016 (v17), Gaussian 09, SPSS, QSARIS và MATLAB:
- Đánh giá thống kê mô hình: Đánh giá qua các chỉ số: hệ số xác định hồi quy $R^2$, hệ số hiệu chỉnh $R^2_{adj}$, tổng bình phương phần dư dự đoán $PRESS$, căn bậc hai sai số trung bình phương nội ($RMSE_{train}, RMSE_{CV}$) và ngoại ($RMSE_{test}$), sai số tuyệt đối trung bình ($AUE$), sai số phần trăm tuyệt đối trung bình ($MAPE, MARE$).
- Kiểm định tính vững (Robustness & Validation): Thực hiện đánh giá chéo nội bộ loại trừ từng cá thể (Leave-One-Out, LOO-CV) đạt $Q^2_{LOO} > 0,90$, đánh giá chéo loại trừ nhiều cá thể (Leave-Many-Out, LMO-CV), và kiểm định ngẫu nhiên hoán vị biến phụ thuộc ($Y$-randomization test) để loại trừ tương quan ngẫu nhiên.
- Xác định miền ứng dụng (Applicability Domain - AD): Sử dụng phân tích đòn bẩy (leverage approach) kết hợp biểu đồ Williams plot để phát hiện các mẫu ngoại lai (outliers) vượt ngưỡng cảnh báo đòn bẩy $h^* = 3(k+1)/n$, đảm bảo tất cả các dẫn xuất thiết kế mới đều nằm trọn vẹn trong không gian mô tả của miền ứng dụng.
- Phân tích phương sai (ANOVA) và kiểm định Student $t$-test ở mức ý nghĩa $p < 0,001$ khẳng định sự khác biệt có ý nghĩa thống kê vượt trội của mô hình QSPR tối ưu so với các mô hình ngẫu nhiên.
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Hiệu năng vượt trội của mô hình GA-MLR và GA-SVR trong dự báo hằng số bền:
Đối với phức đơn nhân $ML$, mô hình QSPR-MLR tối ưu hóa bằng giải thuật di truyền lựa chọn tổ hợp biến đặc trưng đạt hệ số tương quan huấn luyện $R^2_{train} = 0,9446$, sai số chuẩn $s = 0,214$, và hệ số đánh giá chéo $Q^2_{LOO} = 0,9262$. Kiểm định Fisher $F$-test đạt giá trị vượt ngưỡng giới hạn ở độ tin cậy 99% ($p < 0,0001$). Mô hình phi tuyến GA-SVR và GA-ANN ($I(6)-HL(6)-O(1)$) cho thấy sự tương thích hoàn hảo giữa giá trị thực nghiệm và tính toán với $MSE < 0,028$.
-
Giải mã bản chất hóa lượng tử chi phối khả năng tạo phức:
Phân tích đóng góp biến số chứng minh rằng hằng số bền $\log\beta_{11}$ phụ thuộc thuận nghịch có hệ thống vào 4 biến cấu trúc chủ đạo:
- Chỉ số kết nối topo $xp5$ và $xp3$ (phản ánh độ phân nhánh và chiều dài chuỗi liên hợp $C=N-NH-C(=S)-$).
- Số lượng vòng thơm ngưng tụ ($nrings$).
- Độ elip phân tử ($ovality$) biểu thị hình thái khối không gian của phối tử khi tiếp cận ion kim loại trung tâm.
Đặc biệt, năng lượng orbital phân tử $E_{HOMO}$ của phối tử tỷ lệ thuận với hằng số bền phức chất, phản ánh khả năng cho cặp electron không phân chia từ dị nguyên tử $S$ và $N$ vào các orbital $d$ trống của ion kim loại chuyển tiếp.
-
Khám phá rào cản năng lượng quay và cấu dạng không gian tối ưu:
Phân tích quét thế năng bề mặt (PES) ở dải nhiệt độ $298 - 473\text{ K}$ cho thấy hai phối tử BEPT và BECT có 4 góc nhị diện xoay linh động ($\alpha_1: \text{H-N1-C2-N3}$, $\alpha_2: \text{N1-C2-N3-N4}$, $\alpha_3: \text{C2-N3-N4-C5}$, $\alpha_4: \text{N4-C5-C6-C7}$). Hàng rào năng lượng quay xác định cấu dạng phẳng bền vững nhất của BEPT và BECT có năng lượng thấp nhất, tạo điều kiện thuận lợi tuyệt đối cho việc đóng vòng chelat 5 cạnh bền vững với ion kim loại.
-
Tổng hợp thành công 2 phối tử mới và 4 phức chất mới với độ tinh khiết cao:
- Phối tử BEPT được tổng hợp qua 4 giai đoạn từ phenothiazine (ethyl hóa $\rightarrow$ formyl hóa carbonyl $\rightarrow$ brom hóa bằng NBS $\rightarrow$ ngưng tụ với thiosemicarbazide).
- Phối tử BECT được tổng hợp qua 4 giai đoạn từ carbazole tương tự.
- 4 phức chất $\text{Cd(II)-BEPT}$, $\text{Ni(II)-BEPT}$, $\text{Cu(II)-BECT}$, $\text{Zn(II)-BECT}$ được cô lập ở dạng tinh thể màu đặc trưng, bền vững trong không khí và giữ độ bền màu quang phổ vượt quá 120 phút trong dung dịch.
-
Sự trùng khớp xuất sắc giữa hằng số bền thực nghiệm và dự báo lý thuyết:
Bằng phương pháp trắc quang Benesi-Hildebrand kết hợp phương pháp Job và tỷ lệ mol:
- Phức chất $\text{Cd(II)-BEPT}$ và $\text{Ni(II)-BEPT}$ xác định tỷ lệ phối trí $M:L = 1:2$ ($\text{Cd(BEPT)}_2$ và $\text{Ni(BEPT)}2$), hằng số bền thực nghiệm $\log\beta{12}$ lần lượt đạt $11,84 \pm 0,06$ và $12,15 \pm 0,05$.
- Phức chất $\text{Cu(II)-BECT}$ và $\text{Zn(II)-BECT}$ xác định tỷ lệ phối trí $M:L = 1:2$ ($\text{Cu(BECT)}_2$ và $\text{Zn(BECT)}2$), hằng số bền thực nghiệm $\log\beta{12}$ lần lượt đạt $12,42 \pm 0,04$ và $11,68 \pm 0,05$.
Sai số tương đối tuyệt đối ($ARE$) giữa giá trị thực nghiệm và giá trị dự báo từ mô hình QSPR đều nằm dưới $3,8%$, chứng minh độ tin cậy tuyệt đối của quy trình mô hình hóa.
graph LR
A[Mô hình QSPR GA-MLR/SVR] -->|Dự báo log Beta| B[Phức chất mới giả định]
C[Thực nghiệm Trắc quang UV-Vis] -->|Đo log Beta Benesi-Hildebrand| D[Phức chất mới tổng hợp]
B --> E{Đánh giá sai số ARE < 3.8%}
D --> E
E -->|Khẳng định tính đúng đắn| F[Quy trình chuẩn hóa In Silico - In Vitro]
Implications đa chiều
- Ý nghĩa lý thuyết: Cung cấp hệ phương pháp luận hoàn chỉnh kết hợp hóa học lượng tử bán thực nghiệm (PM7) với các thuật toán học máy hiện đại, mở rộng biên giới của lý thuyết liên kết phối trí trong dung dịch.
- Đổi mới phương pháp luận: Thiết lập quy trình chuẩn 6 bước từ sàng lọc ảo đến tổng hợp có định hướng, có thể chuyển giao và áp dụng trực tiếp cho các họ phối tử khác như hydrazone, Schiff base, oxime hoặc porphyrin.
- Ứng dụng thực tiễn trong phân tích trắc quang: Các phối tử BEPT và BECT là các thuốc thử hữu cơ có độ nhạy và độ chọn lọc cực cao, cho phép phát hiện định lượng hàm lượng vết các kim loại nặng độc hại ($\text{Cd}^{2+}, \text{Ni}^{2+}, \text{Cu}^{2+}, \text{Zn}^{2+}$) trong các mẫu nước thải công nghiệp, mẫu thực phẩm và bệnh phẩm y tế bằng thiết bị trắc quang UV-Vis phổ thông giá thành rẻ.
- Ứng dụng trong hóa dược: Mở ra triển vọng phát triển các phức chất kim loại chuyển tiếp mang hoạt tính ức chế tế bào ung thư và kháng khuẩn đa kháng thuốc dựa trên bộ khung dị vòng phenothiazine và carbazole liên hợp thiosemicarbazone.
Limitations và Future Research
Nhằm duy trì tính khách quan và chuẩn mực học thuật, luận án chỉ rõ các giới hạn nội tại:
- Giới hạn tập dữ liệu ban đầu: Bộ dữ liệu thực nghiệm gồm 54 phối tử chủ yếu là các dẫn xuất thiosemicarbazone thơm đơn nhân và nhị nhân; số lượng phối tử ba càng ($ONS, NNS$) và các phức kim loại quý phân nhóm $4d, 5d$ ($\text{Pt}^{2+}, \text{Pd}^{2+}, \text{Ru}^{3+}$) còn khiêm tốn.
- Hiệu ứng dung môi trong tính toán lượng tử: Các mô phỏng lượng tử sử dụng mô hình điện môi liên tục COSMO ở mức bán thực nghiệm PM7; chưa xét đến vi cấu trúc solvat hóa riêng phần của các phân tử dung môi phối trí cạnh tranh trong dung dịch nước thực tế.
- Đặc trưng cấu trúc tinh thể đơn tinh thể: Do hạn chế về nuôi tinh thể đơn, việc xác minh cấu trúc phức chất dựa trên tổ hợp đa phổ ($^1\text{H}/^{13}\text{C}$-NMR, DEPT, 2D-NMR HSQC-HMBC, HR-MS, FT-IR, EDX-SEM) thay vì nhiễu xạ tia X đơn tinh thể (SC-XRD).
Chương trình nghiên cứu tương lai (Future Research Agenda) 5-10 năm tới bao gồm:
- Mở rộng mô hình hóa QSPR sang hệ 3D/4D-QSPR có tính đến động lực học phân tử (Molecular Dynamics - MD) trong dung môi tường minh.
- Ứng dụng lý thuyết phiếm hàm mật độ phụ thuộc thời gian (TD-DFT) ở mức B3LYP/def2-TZVP để mô phỏng chính xác phổ kích thích điện tử UV-Vis của phức chất.
- Khảo sát hoạt tính sinh học in vitro và in vivo (kháng tế bào ung thư dòng MCF-7, HepG2, kháng nấm Candida albicans) của 4 phức chất mới tổng hợp.
- Phát triển các bộ kit trắc quang thương mại hóa và cảm biến quang học (optical chemosensors) cố định trên màng polymer nano để quan trắc ion kim loại nặng hiện trường.
Tác động và ảnh hưởng
- Tác động học thuật (Academic Impact): Đặt nền tảng xuất bản cho hàng loạt công trình trên các tạp chí quốc tế chuyên ngành (ISI/Scopus), định hình một hướng tiếp cận kết hợp QSPR - Lượng tử - Thực nghiệm tại Việt Nam với tiềm năng trích dẫn cao trong lĩnh vực Hóa học tính toán và Hóa phân tích vô cơ.
- Chuyển dịch công nghệ R&D công nghiệp (Industry Transformation): Cung cấp giải pháp cho các viện nghiên cứu và phòng thí nghiệm kiểm nghiệm môi trường quy trình tổng hợp thuốc thử phân tích hữu cơ hiệu năng cao với chi phí giảm 60-80% so với thuốc thử nhập ngoại.
- Chính sách và Lợi ích xã hội: Hỗ trợ các cơ quan quản lý tài nguyên và môi trường (Tổng cục Môi trường, các Trung tâm Quan trắc) phương pháp kiểm soát dư lượng kim loại nặng trong nước thải công nghiệp với độ tin cậy đạt chuẩn quốc tế, góp phần bảo vệ nguồn nước và sức khỏe cộng đồng trước nguy cơ ngộ độc kim loại tích lũy.
- Vị thế quốc tế: Khẳng định năng lực nghiên cứu hóa học tính toán độc lập của các nhà khoa học Việt Nam trên bản đồ học thuật thế giới thông qua việc làm chủ phần mềm và phương pháp luận tiên tiến.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Giới nghiên cứu Hóa học tính toán: Tiếp cận quy trình từng bước xây dựng mô hình QSPR chuẩn mực theo OECD, phương pháp tối ưu hóa hình học bán thực nghiệm PM7 và giải thuật di truyền chọn lọc biến số.
- Giảng viên & Nhà khoa học Hóa vô cơ - Phân tích: Nắm bắt phương pháp luận tích hợp giữa mô phỏng in silico và phổ nghiệm đa chiều hiện đại trong xác định cơ chế tạo phức và tính toán hằng số bền.
- Kỹ sư R&D Hóa dược & Môi trường: Sở hữu quy trình tổng hợp và công thức điều chế 2 phối tử mới (BEPT, BECT) cùng các phức chất kim loại tương ứng ứng dụng làm chất kháng khuẩn và thuốc thử trắc quang vi lượng.
- Cơ quan Kiểm nghiệm & Quan trắc Môi trường: Tiếp nhận phương pháp phân tích trắc quang đơn giản, rẻ tiền, nhanh chóng với độ lặp lại cao để quan trắc ion $\text{Cd}^{2+}, \text{Ni}^{2+}, \text{Cu}^{2+}, \text{Zn}^{2+}$.
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 mở rộng Thuyết tương quan năng lượng tự do tuyến tính (LFER) của Hammett-Taft và thuyết phối trí vào việc lượng hóa hằng số bền phức chất kim loại - thiosemicarbazone. Luận án chứng minh rằng năng lượng tạo phức không chỉ do điện tích hình thức quyết định mà là hàm đa biến đồng thời của chỉ số kết nối topo không gian ($xp3, xp5$), độ elip ($ovality$) và năng lượng orbital biên ($E_{HOMO}$), giải quyết triệt để vấn đề dự báo ái lực liên kết trước thực nghiệm.
2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?
Trả lời: So với nghiên cứu của Wang và cộng sự (2014) [156] (chỉ dừng ở mô hình 3D-QSAR tĩnh không có kiểm chứng tạo phức thực nghiệm) và nghiên cứu của Reddy và cộng sự (2013) [113] (chỉ tổng hợp thực nghiệm thuốc thử trắc quang thuần túy bằng phương pháp thử - sai), luận án tạo ra bước đột phá về phương pháp luận khi thiết lập chu trình "Khám phá ảo $\rightarrow$ Dự báo lượng tử $\rightarrow$ Thực nghiệm hóa học $\rightarrow$ Hồi quy phản hồi". Đây là nghiên cứu đầu tiên trên thế giới xây dựng thành công mô hình QSPR cho hằng số bền phức kim loại - thiosemicarbazone đáp ứng toàn diện 5 nguyên tắc OECD.
3. Phát hiện bất ngờ nhất (surprising finding) thu được từ dữ liệu là gì?
Trả lời: Phát hiện bất ngờ nhất là ảnh hưởng vượt trội của chỉ số lập thể $ovality$ (độ elip phân tử) và số vòng thơm ngưng tụ ($nrings$) so với các tham số điện tử thuần túy trong việc quyết định hằng số bền của phức bậc hai $ML_2$. Mặc dù mật độ điện tích tại tâm $S$ và $N$ quyết định liên kết phối trí ban đầu, nhưng chính độ cồng kềnh và hình học không gian của phối tử mới là yếu tố quyết định hàng rào năng lượng cản trở lập thể khi phân tử phối tử thứ hai tiếp cận tạo cấu hình bát diện hoặc vuông phẳng.
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Trả lời: Luận án cung cấp quy trình tái lập chi tiết ở mức độ tuyệt đối: từ danh mục 230 biến số toán học, thông số cài đặt giải thuật di truyền (GA) và tham số mạng nơ-ron ANN, đến quy trình tổng hợp hóa học chi tiết từng gam hóa chất, nhiệt độ, dung môi, chất xúc tác, cùng toàn bộ dữ liệu phổ nghiệm gốc ($^1\text{H}$-NMR, $^{13}\text{C}$-NMR, DEPT, HSQC, HMBC, HR-MS, FT-IR, UV-Vis) của BEPT, BECT và 4 phức chất tương ứng.
5. Lộ trình nghiên cứu 10 năm (10-year research agenda) được phác thảo như thế nào?
Trả lời: Lộ trình 10 năm định hướng: (1) Nâng cấp mô hình QSPR lên nền tảng trí tuệ nhân tạo học sâu (Deep Learning/Graph Neural Networks) kết hợp mô phỏng động lực học lượng tử đa quy mô QM/MM; (2) Mở rộng thư viện cấu trúc sang phức chất của các nguyên tố phóng xạ và kim loại nhóm platin ứng dụng trong xạ trị và hóa trị ung thư; (3) Tự động hóa quy trình tổng hợp dòng chảy liên tục (continuous flow synthesis) dựa trên dự báo in silico.
Kết luận
Công trình luận án tiến sĩ của NCS. Nguyễn Minh Quang đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp khoa học cốt lõi:
- Xây dựng thành công hệ thống mô hình QSPR dự báo hằng số bền phức chất giữa thiosemicarbazone và các ion kim loại chuyển tiếp với độ chính xác cao ($R^2_{train} = 0,9446; Q^2_{LOO} = 0,9262; MSE < 0,028$), thỏa mãn đầy đủ các tiêu chuẩn quốc tế nghiêm ngặt của OECD.
- Lượng hóa và giải mã thành công các yếu tố cấu trúc và hóa lượng tử bán thực nghiệm PM7 ($xp5, xp3, nrings, ovality, E_{HOMO}$) đóng vai trò then chốt chi phối độ bền nhiệt động học của phức chất dạng $ML$ và $ML_2$.
- Sàng lọc ảo 44 dẫn xuất thiosemicarbazone thiết kế mới, mô phỏng chính xác bề mặt thế năng (PES) và các hàng rào năng lượng quay góc nhị diện, xác lập cấu dạng phẳng bền vững nhất của các phối tử dị vòng mới.
- Tổng hợp thành công bằng con đường hóa học 2 phối tử mới chứa dị vòng thơm đa nhân gồm BEPT và BECT qua quy trình tổng hợp hữu cơ 4 giai đoạn với hiệu suất cao.
- Tổng hợp và cô lập thành công 4 phức chất kim loại mới gồm $\text{Cd(II)-BEPT}$, $\text{Ni(II)-BEPT}$, $\text{Cu(II)-BECT}$, $\text{Zn(II)-BECT}$; giải đoán sáng tỏ cấu trúc không gian và cơ chế tạo phức hai càng ($N, S$) thông qua hệ thống phổ nghiệm hiện đại 1D/2D-NMR, HR-MS, FT-IR, EDX và SEM.
- Xác định chính xác thành phần và hằng số bền quang phổ ($\log\beta_{12}$ từ $11,68$ đến $12,42$) bằng phương pháp Benesi-Hildebrand, Job và tỷ lệ mol, qua đó kiểm chứng thực nghiệm tính đúng đắn và độ tin cậy vượt bậc của các mô hình QSPR đã phát triển với sai số $ARE < 3,8%$.
Luận án khẳng định bước tiến vượt bậc trong việc hiện đại hóa nghiên cứu hóa vô cơ - hóa lý thuyết tại Việt Nam, mở ra 3 hướng nghiên cứu mới về cảm biến hóa học phân tích môi trường, vật liệu phức chất quang hoạt và dược chất phức kim loại hướng đích trong tương lai.
