Giới thiệu dự án
Trong ngành công nghiệp dược phẩm và kiểm nghiệm hóa phân tích hiện đại, việc kiểm soát chất lượng các chế phẩm đa thành phần đóng vai trò then chốt nhằm đảm bảo hiệu quả điều trị và an toàn sinh học. Các dạng thuốc phối hợp kháng khuẩn, kháng nấm và kháng viêm – điển hình như viên nén phụ khoa Mycogynax – chứa đồng thời ba hoạt chất có hoạt tính dược lý mạnh gồm Metronidazole (200 mg), Chloramphenicol (80 mg) và Dexamethasone acetate (0,5 mg).
+-------------------------------------------------------------+
| THÁCH THỨC ĐỊNH LƯỢNG HỖN HỢP |
+-------------------------------------------------------------+
|
+------------------------------+------------------------------+
| |
v v
+---------------------------------------------------------+ +---------------------------------------------------------+
| SẮC KÝ LỎNG HPLC | | QUANG PHỔ UV-VIS TRUYỀN THỐNG |
| • Chi phí thiết bị và dung môi pha động rất đắt đỏ | | • Phổ hấp thụ của 3 hoạt chất chồng lấn hoàn toàn |
| • Thời gian chuẩn bị mẫu và chạy sắc ký kéo dài | | • Định luật Beer cổ điển không giải được hệ phương trình|
| • Tạo ra lượng lớn chất thải hóa chất độc hại | | • Yêu cầu chiết tách hóa học phức tạp, dễ gây sai số |
+---------------------------------------------------------+ +---------------------------------------------------------+
| |
+------------------------------+------------------------------+
|
v
+-------------------------------------------------------------+
| GIẢI PHÁP ĐỀ TÀI |
| QUANG PHỔ HẤP THỤ UV-VIS + THUẬT TOÁN PLS2 |
| • Không cần tách pha – Xử lý phổ chồng lấn bằng Chemometrics|
| • Chi phí vận hành thấp – Thời gian phân tích < 10 phút |
| • Thuật toán NIPALS tối ưu hóa trên nền tảng C++ |
+-------------------------------------------------------------+
Vấn đề nghiên cứu và Điểm nghẽn kỹ thuật (Problem Statement)
Phương pháp tiêu chuẩn hiện nay tại các trung tâm kiểm nghiệm là sắc ký lỏng hiệu năng cao (HPLC). Mặc dù HPLC đạt độ chọn lọc cao, phương pháp này đòi hỏi chi phí đầu tư thiết bị rất lớn, tiêu tốn dung môi hữu cơ tinh khiết cao (Methanol, Acetonitrile pha động) và quy trình tiền xử lý mẫu phức tạp. Ngược lại, phương pháp quang phổ hấp thụ phân tử UV-Vis là công cụ phổ biến, thao tác nhanh, chi phí vận hành thấp nhưng gặp phải rào cản kỹ thuật nghiêm trọng: phổ hấp thụ tử ngoại của Metronidazole, Chloramphenicol và Dexamethasone acetate chồng lấn sâu sắc trong vùng bước sóng 220 – 350 nm. Khi đó, việc áp dụng hệ phương trình cộng tính quang phổ thông thường qua định luật Bouguer – Lambert – Beer gặp hiện tượng đa cộng tuyến, sai số tích lũy lớn và không thể định lượng độc lập từng chất nếu không qua bước tách chiết hóa học độc hại.
Mục tiêu của dự án
- Khảo sát toàn diện các điều kiện hóa lý tối ưu (hệ dung môi, pH đệm, độ bền quang phổ theo thời gian, khoảng nồng độ tuyến tính tuân theo định luật Beer) cho cả 3 hoạt chất riêng rẽ và hỗn hợp.
- Xây dựng thuật toán hồi quy đa biến bình phương tối thiểu riêng phần đa biến phụ thuộc (Partial Least Squares 2 – PLS2) dựa trên quy trình lặp NIPALS (Non-linear Iterative Partial Least Squares).
- Lập trình module phần mềm chuyên dụng bằng ngôn ngữ C++ trên môi trường Microsoft Visual C++ 6.0 để xử lý trực tiếp ma trận phổ đa bước sóng ($m \times 27$).
- Kiểm định mô hình thông qua phương pháp dự báo hợp thức chéo (Leave-One-Out Cross-Validation) với chỉ số thống kê PRESS và RMS.
- Ứng dụng quy trình định lượng đồng thời 3 hoạt chất trong chế phẩm thực tế (viên nén Mycogynax) và đối chứng với kết quả công bố của nhà sản xuất cũng như Trung tâm Dịch vụ Phân tích Thí nghiệm TP.HCM (CASE).
Phương pháp tiếp cận và Biện giải
Dự án tích hợp phương pháp quang phổ UV-Vis với công cụ hóa tin học (Chemometrics). Thuật toán PLS2 đồng thời giải mã cấu trúc phương sai của ma trận độ hấp thụ quang $\mathbf{X}$ và ma trận nồng độ $\mathbf{Y}$, chiếu các biến độc lập ban đầu lên các biến ẩn (latent variables/scores) trực giao. Cách tiếp cận này loại bỏ hoàn toàn hiện tượng đa cộng tuyến giữa 27 bước sóng khảo sát, lọc bỏ nhiễu nền phổ và giữ lại tối đa thông tin tương quan có ý nghĩa định lượng.
Kết quả kỳ vọng và Phạm vi áp dụng
- Chỉ tiêu đo lường: Tỷ lệ thu hồi hoạt chất đạt $98,0% - 102,0%$; sai số tương đối giữa nồng độ thực nghiệm và nồng độ lý thuyết $< 3,0%$; chỉ số sai số dự báo $RMS < 0,05$.
- Phạm vi nghiên cứu: Định lượng đồng thời hệ 3 thành phần Metronidazole, Chloramphenicol, Dexamethasone acetate trong môi trường dung dịch đệm cồn-nước ở bước sóng 230 – 350 nm. Không bao gồm định lượng thành phần tá dược trơ và Nystatin (hoạt chất kháng nấm được loại bỏ trong pha xử lý sơ bộ).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
| Phương pháp phân tích |
Ưu điểm kỹ thuật |
Nhược điểm / Rào cản |
Chi phí & Thời gian |
Độ chính xác |
| Sắc ký lỏng HPLC |
Tách triệt để các cấu tử, độ nhạy cao, độ đặc hiệu chuẩn tuyệt đối |
Cần cột tách C18 đắt tiền, tiêu tốn dung môi gradient, bảo trì phức tạp |
Rất cao (~500.000 VNĐ/mẫu); 30 - 45 phút/mẫu |
Rất cao (Sai số $< 1%$) |
| Trắc quang đơn bước sóng cổ điển |
Thao tác đơn giản, thiết bị đo UV-Vis phổ biến, giá thành rẻ |
Hoàn toàn bất khả thi khi các đỉnh hấp thụ quang phổ chồng lấn nhau |
Thấp (~20.000 VNĐ/mẫu); 5 phút/mẫu |
Sai số không chấp nhận được ($> 25%$) |
| Quang phổ đạo hàm bậc 1 - bậc 2 |
Triệt tiêu được nhiễu nền, tách được một số đỉnh phổ xen phủ |
Độ nhạy suy giảm ở đạo hàm bậc cao, khuếch đại nhiễu tín hiệu quang học |
Thấp; 10 - 15 phút/mẫu |
Trung bình (Sai số $3 - 7%$) |
| Hồi quy PLS2 kết hợp UV-Vis (Đề xuất) |
Không cần tách pha, khai thác dữ liệu toàn dải phổ (27 bước sóng), khử nhiễu đa biến |
Mô hình toán học phức tạp, đòi hỏi phần mềm tính toán ma trận ma trận tải |
Rất thấp (~25.000 VNĐ/mẫu); 8 - 10 phút/mẫu |
Cao (Sai số $< 2,5%$) |
Phân tích yêu cầu hệ thống theo mô hình MoSCoW
- Must have (Bắt buộc):
- Khảo sát được khoảng nồng độ tuyến tính tuân theo định luật Beer: Metronidazole ($2 - 26\ \mu\text{g/mL}$), Chloramphenicol ($2 - 42\ \mu\text{g/mL}$), Dexamethasone acetate ($1 - 22\ \mu\text{g/mL}$).
- Giải thuật NIPALS cho bài toán PLS2 với ma trận $Y$ đa biến (3 cột nồng độ).
- Thuật toán tự động tìm số cấu tử chính tối ưu $n_c$ thông qua cực tiểu hóa hàm PRESS.
- Should have (Nên có):
- Module C++ chuẩn hóa dữ liệu tự động (Mean Centering và Auto-scaling).
- Xuất bảng tương quan nồng độ thực tế ($C_{TT}$) và nồng độ dự báo ($C_{DB}$).
- Could have (Có thể có):
- Giao diện đồ họa biểu diễn đồ thị hồi quy trực tiếp trên Visual C++.
- Won't have (Chưa thực hiện đợt này):
- Phân tích phi tuyến tính bằng mạng nơ-ron nhân tạo (ANN-PLS) cho hỗn hợp trên 5 cấu tử.
Thiết kế hệ thống
Kiến trúc ngăn xếp công nghệ (Technology Stack)
- Hệ quang phổ: Máy quang phổ UV-Vis quét chùm tia kép (Double Beam Spectrophotometer), cuvet thạch anh quang học $l = 1\text{ cm}$, dải bước sóng 200 – 400 nm.
- Môi trường phát triển: Microsoft Visual C++ 6.0 Standard Edition, trình biên dịch C++ Native x86.
- Cấu trúc dữ liệu & Thiết kế ma trận:
Matrix_A[m][27]: Lưu trữ giá trị mật độ quang học $A$ tại 27 bước sóng phân tích.Matrix_C[m][3]: Lưu trữ nồng độ 3 cấu tử trong $m$ mẫu chuẩn calibration.Vector_W[27], Vector_T[m], Vector_P[27], Vector_Q[3]: Các vector trọng số, vector vết và vector tải trong vòng lặp NIPALS.
Phương pháp luận và Quy trình triển khai
Quy trình nghiên cứu thực nghiệm
Giai đoạn 1: Khảo sát hóa lý & Tối ưu hóa điều kiện đo quang (Tuần 1 - 4)
Giai đoạn 2: Xây dựng thuật toán & Lập trình Module PLS2/C++ (Tuần 5 - 8)
Giai đoạn 3: Thực nghiệm trên mẫu tự pha & Hiệu chỉnh mô hình (Tuần 9 - 11)
Giai đoạn 4: Ứng dụng trên mẫu thuốc thật Mycogynax & Đối chuẩn (Tuần 12 - 14)
Cơ sở toán học của thuật toán PLS2 NIPALS
Mô hình toán học phân rã đồng thời hai ma trận dữ liệu $\mathbf{A}{m \times p}$ (ma trận phổ) và $\mathbf{C}{m \times n}$ (ma trận nồng độ):
$$\mathbf{A} = \mathbf{T}\mathbf{P}^T + \mathbf{E} = \sum_{h=1}^{a} \mathbf{t}_h \mathbf{p}_h^T + \mathbf{E}$$
$$\mathbf{C} = \mathbf{U}\mathbf{Q}^T + \mathbf{F} = \sum_{h=1}^{a} \mathbf{u}_h \mathbf{q}_h^T + \mathbf{F}$$
Mối liên hệ nội tại giữa biến ẩn phổ $\mathbf{t}_h$ và biến ẩn nồng độ $\mathbf{u}_h$ được thiết lập qua phương trình hồi quy tuyến tính:
$$\mathbf{u}_h = b_h \mathbf{t}_h \quad \text{với} \quad b_h = \frac{\mathbf{u}_h^T \mathbf{t}_h}{\mathbf{t}_h^T \mathbf{t}_h}$$
Chỉ số tổng bình phương sai số dự báo (PRESS - Prediction Error Sum of Squares) được xác định theo công thức:
$$\text{PRESS} = \sum_{i=1}^{m} \sum_{j=1}^{n} (c_{ij} - \hat{c}_{ij})^2$$
Trong đó $c_{ij}$ là nồng độ thực tế của cấu tử $j$ trong mẫu $i$, và $\hat{c}_{ij}$ là nồng độ dự báo khi loại bỏ chính mẫu $i$ khỏi quá trình lập chuẩn.
Hiện thực hóa và Kết quả thực nghiệm
Đoạn mã thực thi thuật toán lõi PLS2 NIPALS trên C++
Dưới đây là đoạn mã nguồn C++ hiện thực hóa vòng lặp NIPALS trích xuất từ chương trình tính toán của đề tài:
#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
const int MAX_SAMPLES = 50;
const int MAX_WAVELENGTHS = 27;
const int MAX_COMPONENTS = 3;
const double EPSILON = 1e-7;
// Ham thuc thi vong lap NIPALS PLS2 cho mot cau tu an (Latent Variable)
bool CalculatePLS2Component(
int m, int p, int n,
double A[MAX_SAMPLES][MAX_WAVELENGTHS],
double C[MAX_SAMPLES][MAX_COMPONENTS],
double w[MAX_WAVELENGTHS],
double t[MAX_SAMPLES],
double p_load[MAX_WAVELENGTHS],
double q_load[MAX_COMPONENTS],
double &b_coeff
) {
double u[MAX_SAMPLES];
double t_old[MAX_SAMPLES];
// Khoi tao vector u bang cot dau tien cua ma tran C
for (int i = 0; i < m; ++i) {
u[i] = C[i][0];
}
int iter = 0;
while (iter < 500) {
// 1. Tinh vector trong so w = A^T * u / (u^T * u)
double u_norm_sq = 0.0;
for (int i = 0; i < m; ++i) u_norm_sq += u[i] * u[i];
for (int j = 0; j < p; ++j) {
w[j] = 0.0;
for (int i = 0; i < m; ++i) {
w[j] += A[i][j] * u[i];
}
w[j] /= (u_norm_sq + 1e-12);
}
// 2. Chuan hoa vector w: w = w / ||w||
double w_norm = 0.0;
for (int j = 0; j < p; ++j) w_norm += w[j] * w[j];
w_norm = sqrt(w_norm);
for (int j = 0; j < p; ++j) w[j] /= w_norm;
// 3. Tinh vector vet t = A * w / (w^T * w) -> w^T * w = 1
for (int i = 0; i < m; ++i) {
t_old[i] = t[i];
t[i] = 0.0;
for (int j = 0; j < p; ++j) {
t[i] += A[i][j] * w[j];
}
}
// 4. Tinh vector tai q = C^T * t / (t^T * t)
double t_norm_sq = 0.0;
for (int i = 0; i < m; ++i) t_norm_sq += t[i] * t[i];
for (int k = 0; k < n; ++k) {
q_load[k] = 0.0;
for (int i = 0; i < m; ++i) {
q_load[k] += C[i][k] * t[i];
}
q_load[k] /= (t_norm_sq + 1e-12);
}
// 5. Tinh vector u moi: u = C * q / (q^T * q)
double q_norm_sq = 0.0;
for (int k = 0; k < n; ++k) q_norm_sq += q_load[k] * q_load[k];
for (int i = 0; i < m; ++i) {
u[i] = 0.0;
for (int k = 0; k < n; ++k) {
u[i] += C[i][k] * q_load[k];
}
u[i] /= (q_norm_sq + 1e-12);
}
// Kiem tra hoi tu dua tren su thay doi cua vector t
double diff = 0.0;
for (int i = 0; i < m; ++i) diff += pow(t[i] - t_old[i], 2);
if (iter > 0 && sqrt(diff) < EPSILON) break;
iter++;
}
// 6. Tinh vector tai p = A^T * t / (t^T * t)
double t_sq = 0.0;
for (int i = 0; i < m; ++i) t_sq += t[i] * t[i];
for (int j = 0; j < p; ++j) {
p_load[j] = 0.0;
for (int i = 0; i < m; ++i) p_load[j] += A[i][j] * t[i];
p_load[j] /= (t_sq + 1e-12);
}
// 7. He so hoi quy noi tai b = (u^T * t) / (t^T * t)
double ut_dot = 0.0;
for (int i = 0; i < m; ++i) ut_dot += u[i] * t[i];
b_coeff = ut_dot / (t_sq + 1e-12);
return true;
}
Kết quả kiểm định trên mẫu tự pha (Synthetic Mixtures)
Quá trình hợp thức chéo (Cross-validation) trên 16 mẫu hỗn hợp tự pha bậc 3 với các tỷ lệ nồng độ thay đổi ngẫu nhiên cho thấy giá trị PRESS thay đổi rõ rệt theo số lượng cấu tử chính $n_c$:
Giá trị PRESS theo số cấu tử chính (nc):
nc = 1: PRESS = 142.85 (Underfitting - Thiếu biến ẩn)
nc = 2: PRESS = 18.42 (Mô hình bắt đầu ổn định)
nc = 3: PRESS = 2.15 (CỰC TIỂU TỐI ƯU - PRESS Min)
nc = 4: PRESS = 3.89 (Overfitting - Bắt đầu thu nạp nhiễu)
nc = 5: PRESS = 7.64 (Nhiễu nền quang phổ khuếch đại)
Bảng kết quả định lượng trên tập mẫu tự pha ($n_c = 3$)
| Mẫu |
$C_{LT}$ Metro ($\mu\text{g/mL}$) |
$C_{TN}$ Metro ($\mu\text{g/mL}$) |
Sai số (%) |
$C_{LT}$ Chlor ($\mu\text{g/mL}$) |
$C_{TN}$ Chlor ($\mu\text{g/mL}$) |
Sai số (%) |
$C_{LT}$ Dexa ($\mu\text{g/mL}$) |
$C_{TN}$ Dexa ($\mu\text{g/mL}$) |
Sai số (%) |
| M1 |
10.00 |
10.08 |
+0.80 |
4.00 |
3.96 |
-1.00 |
0.025 |
0.0252 |
+0.80 |
| M2 |
15.00 |
14.85 |
-1.00 |
6.00 |
6.07 |
+1.17 |
0.037 |
0.0366 |
-1.08 |
| M3 |
20.00 |
20.14 |
+0.70 |
8.00 |
7.91 |
-1.12 |
0.050 |
0.0504 |
+0.80 |
| M4 |
5.00 |
5.06 |
+1.20 |
10.00 |
9.88 |
-1.20 |
0.020 |
0.0198 |
-1.00 |
| M5 |
12.00 |
11.91 |
-0.75 |
5.00 |
5.04 |
+0.80 |
0.030 |
0.0303 |
+1.00 |
| TB |
— |
— |
0.89 |
— |
— |
1.06 |
— |
— |
0.90 |
Ghi chú: $C_{LT}$: Nồng độ lý thuyết pha chế; $C_{TN}$: Nồng độ tính toán từ mô hình PLS2.
Kết quả định lượng trên mẫu thuốc thật Mycogynax
Tiến hành phân tích 10 lô mẫu chế phẩm viên nén phụ khoa Mycogynax trên thị trường. Kết quả định lượng được quy đổi về hàm lượng hoạt chất trong 1 viên nén và so sánh với phương pháp HPLC tiêu chuẩn tại Trung tâm Phân tích Thí nghiệm TP.HCM (CASE):
| Hoạt chất phân tích |
Hàm lượng nhãn (mg/viên) |
Phương pháp PLS2 (mg/viên) |
Phương pháp HPLC / CASE (mg/viên) |
Độ thu hồi PLS2 (%) |
Độ lệch so với HPLC (%) |
| Metronidazole |
200,0 |
$198,65 \pm 1,42$ |
$199,10 \pm 0,85$ |
$99,32%$ |
$-0,23%$ |
| Chloramphenicol |
80,0 |
$79,12 \pm 0,68$ |
$79,50 \pm 0,42$ |
$98,90%$ |
$-0,48%$ |
| Dexamethasone acetate |
0,50 |
$0,494 \pm 0,008$ |
$0,498 \pm 0,004$ |
$98,80%$ |
$-0,80%$ |
Kiểm định thống kê Student ($t$-test) tại độ tin cậy $95%$ ($\alpha = 0,05$) cho thấy giá trị $t_{thực\ nghiệm} < t_{bảng}\ (2,262)$, khẳng định không có sự khác biệt có ý nghĩa thống kê giữa phương pháp quang phổ UV-Vis PLS2 và phương pháp HPLC đối chuẩn.
Đổi mới và Đóng góp khoa học
- Giải quyết triệt để bài toán phổ hấp thụ chồng lấn cao: Khác với phương pháp cổ điển OLS (Ordinary Least Squares) dễ bị suy biến ma trận nghịch đảo khi các dải phổ phụ thuộc tuyến tính, thuật toán PLS2 trích xuất các biến ẩn tối đa hóa hiệp phương sai giữa phổ và nồng độ, cho phép định lượng chính xác 3 hoạt chất chênh lệch nồng độ tới 400 lần (Metronidazole 200 mg so với Dexamethasone acetate 0,5 mg).
- So sánh ưu thế với các thuật toán hồi quy kinh điển:
- So với PCR (Principal Component Regression): PCR chỉ tối đa hóa phương sai trên ma trận $\mathbf{A}$ độc lập, dễ loại bỏ nhầm các thành phần phổ có cường độ nhỏ nhưng mang thông tin nồng độ cốt lõi (như của Dexamethasone acetate). PLS2 tối ưu hóa đồng thời cả $\mathbf{A}$ và $\mathbf{C}$, giữ trọn vẹn tín hiệu hữu ích.
- So với ILS (Inverse Least Squares): ILS giới hạn số bước sóng phải nhỏ hơn số mẫu thí nghiệm ($p < m$). PLS2 cho phép $p = 27 \gg m = 16$ mà không gặp hiện tượng quá khớp (overfitting).
- Hiệu quả kinh tế và Môi trường (Green Chemistry): Giảm $85%$ lượng dung môi độc hại thải ra môi trường so với HPLC, rút ngắn thời gian phân tích từ 45 phút xuống dưới 10 phút/mẫu ($77,7%$ thời gian), tiết kiệm hơn $90%$ chi phí hóa chất kiểm nghiệm.
Ứng dụng thực tế và Triển khai công nghiệp
Kịch bản ứng dụng thực tế
- Kiểm soát chất lượng tại nhà máy dược phẩm (In-process QA/QC): Kiểm tra độ đồng đều hàm lượng của viên nén ngay trên dây chuyền sản xuất thuốc phụ khoa (Merynal, Megyfar, Mycogynax, Megyna) mà không cần chờ kết quả sắc ký chậm trễ.
- Phân viện kiểm nghiệm và Thanh tra Dược: Triển khai các trạm kiểm nghiệm lưu động nhanh, phát hiện thuốc giả, thuốc kém chất lượng với chi phí đầu tư thiết bị thấp.
+-------------------------------------------------------------+
| QUY TRÌNH TRIỂN KHAI THỰC TẾ |
+-------------------------------------------------------------+
|
+------------------------------+------------------------------+
| |
v v
+---------------------------------------------------------+ +---------------------------------------------------------+
| BƯỚC 1: TIỀN XỬ LÝ MẪU | | BƯỚC 2: QUÉT PHỔ UV-VIS |
| • Cân 20 viên nén, tính khối lượng trung bình | | • Chuyển dung dịch vào cuvet thạch anh l = 1cm |
| • Nghiền mịn, hòa tan bằng Ethanol 96% | | • Quét phổ hấp thụ từ 230 nm đến 350 nm |
| • Lắc siêu âm 15 phút, định mức bằng đệm pH 9,0 | | • Xuất file dữ liệu mật độ quang (.TXT / .CSV) |
| • Lọc qua màng micropore 0.45 um loại bỏ tá dược trơ | | • Tự động trích xuất 27 bước sóng phân tích |
+---------------------------------------------------------+ +---------------------------------------------------------+
|
v
+---------------------------------------------------------+
| BƯỚC 3: DỰ BÁO QUA SOFTWARE PLS2 |
| • Nạp ma trận phổ vào module C++ thực thi |
| • Tự động nhân ma trận hồi quy chuẩn hóa B_pls |
| • Xuất kết quả hàm lượng (mg/viên) sau 3 giây |
+---------------------------------------------------------+
Yêu cầu triển khai hệ thống
- Phần cứng: Máy quang phổ UV-Vis có cổng kết nối RS-232/USB; Máy vi tính văn phòng (tối thiểu CPU Pentium III, 128MB RAM, HĐH Windows 98/2000/XP/7/10/Linux).
- Phần mềm: File thực thi standalone
PLS2_Analyzer.exe (không yêu cầu cài đặt môi trường MATLAB cồng kềnh, dung lượng $< 2\text{ MB}$).
Hạn chế và Hướng phát triển
Hạn chế kỹ thuật
- Nhiễu tán xạ quang học do tá dược không tan: Mặc dù tá dược trơ (Tinh bột ngô, Microcrystalline cellulose, Magnesium stearate) không tan trong cồn nhưng nếu quy trình lọc qua màng $0,45\ \mu\text{m}$ không kỹ, hiện tượng tán xạ ánh sáng sẽ làm trượt đường nền phổ (baseline drift).
- Giới hạn nồng độ định luật Beer: Khi nồng độ Metronidazole vượt quá $30\ \mu\text{g/mL}$, hiện tượng bất đối xứng phân tử và tương tác quang hóa bắt đầu xuất hiện, làm mất tính tuyến tính của vector vết $\mathbf{t}$.
Hướng nâng cấp đề xuất
- Tích hợp giải thuật mạng nơ-ron phi tuyến: Phát triển module lai ANN-PLS hoặc Kernel-PLS (K-PLS) để mô hình hóa chính xác cả các vùng nồng độ cao vượt ra ngoài định luật Beer.
- Hiện đại hóa giao diện người dùng: Chuyển đổi mã nguồn C++ sang nền tảng C# .NET Core WPF hoặc Python PyQt6, hỗ trợ vẽ biểu đồ 3D tương quan phổ - nồng độ trực quan và tích hợp trực tiếp giao thức truyền dữ liệu chuẩn LIMS (Laboratory Information Management System).
Đối tượng hưởng lợi từ công trình
- Sinh viên & Học viên chuyên ngành Hóa học/Dược học: Tiếp cận tài liệu thực chứng về Chemometrics, hiểu rõ cách kết nối giữa lý thuyết toán đại số tuyến tính ma trận với các bài toán phân tích thực tế.
- Kỹ sư phần mềm & Nhà phát triển Chemometrics: Nắm bắt thuật toán NIPALS PLS2 viết bằng ngôn ngữ lập trình bậc thấp C++, dễ dàng tối ưu hóa tính toán ma trận mà không bị phụ thuộc vào các thư viện thương mại độc quyền.
- Doanh nghiệp sản xuất dược phẩm: Ứng dụng quy trình kiểm nghiệm thay thế HPLC cho các dòng sản phẩm viên đặt phụ khoa đa thành phần, giảm chi phí vận hành hàng năm từ 150 – 200 triệu VNĐ trên mỗi dây chuyền sản xuất.
- Cơ quan nghiên cứu khoa học: Cung cấp bộ dữ liệu thực nghiệm chuẩn và phương pháp luận mở để mở rộng định lượng cho các hệ phức tạp 4 – 5 hoạt chất (như bổ sung thêm Nystatin, Neomycin).
Câu hỏi thường gặp (FAQ)
1. Yêu cầu kỹ thuật tối thiểu để triển khai phần mềm PLS2/C++ là gì?
Phần mềm được biên dịch dưới dạng mã máy nhị phân C++ gốc độc lập nên cực kỳ nhẹ. Hệ thống chỉ yêu cầu máy tính chạy hệ điều hành Windows XP trở lên, tương thích với mọi máy quang phổ UV-Vis có khả năng xuất dữ liệu dạng bảng điểm (ASCII/CSV). Không cần cài đặt MATLAB hay bất kỳ Runtime Engine bổ trợ nào.
2. Giới hạn số cấu tử tối đa mà mô hình PLS2 có thể định lượng đồng thời là bao nhiêu?
Về mặt toán học, thuật toán PLS2 không giới hạn số lượng cấu tử $n$. Tuy nhiên, trong thực nghiệm quang phổ hấp thụ, giới hạn tối ưu thường là từ $3 - 6$ cấu tử. Khi số cấu tử vượt quá 6, hiện tượng xen phủ quang phổ quá dày đặc kết hợp với sai số thiết bị sẽ làm tăng nhanh chỉ số RMS, đòi hỏi phải tăng số lượng mẫu lập chuẩn calibration ($m > 50$).
3. Phương pháp này xử lý như thế nào đối với các tá dược có trong viên thuốc?
Các tá dược thông thường trong viên nén Mycogynax (như lactose, tinh bột ngô, magnesium stearate) phần lớn không tan trong ethanol $96%$ hoặc không có nhóm mang màu (chromophore) hấp thụ bức xạ trong vùng $230 - 350\text{ nm}$. Các hạt lơ lửng gây tán xạ ánh sáng được loại bỏ triệt để qua màng lọc micro $0,45\ \mu\text{m}$. Ngoài ra, bước trừ phổ mẫu trắng (Blank Correction) trong thuật toán giúp triệt tiêu hoàn toàn tín hiệu nền của tá dược tan.
4. Chi phí bảo trì và kiểm định định kỳ mô hình PLS2 như thế nào?
Chi phí bảo trì phần mềm gần như bằng 0. Về mặt hóa học, mô hình hiệu chuẩn (Calibration Matrix) cần được kiểm chứng lại sau mỗi 6 tháng hoặc khi thay đổi lô dung môi, thay đèn nguồn Deuterium/Tungsten của máy đo quang bằng cách chạy lại tập 5 mẫu kiểm định chuẩn (Validation set).
5. Thời gian hoàn vốn (ROI) khi chuyển đổi từ HPLC sang UV-Vis PLS2 là bao lâu?
Một phòng thí nghiệm kiểm nghiệm trung bình thực hiện 1.000 mẫu/năm. Chi phí phân tích HPLC ước tính 500.000 VNĐ/mẫu (bao gồm khấu hao cột sắc ký, dung môi pha động HPLC grade, điện năng). Chi phí UV-Vis PLS2 chỉ khoảng 25.000 VNĐ/mẫu. Mức tiết kiệm hàng năm đạt gần 475.000.000 VNĐ, giúp doanh nghiệp thu hồi toàn bộ chi phí nghiên cứu và thiết lập quy trình chỉ sau 3 - 5 tháng vận hành.
Kết luận
Đề tài nghiên cứu đã chứng minh tính khả thi, độ tin cậy và hiệu quả vượt trội của việc kết hợp phương pháp trắc quang hấp thụ phân tử UV-Vis với thuật toán hồi quy đa biến PLS2 trong phân tích dược phẩm. Bằng việc xây dựng thành công module tính toán độc lập trên nền tảng C++ (NIPALS Algorithm), công trình đã giải quyết trọn vẹn điểm nghẽn phổ chồng lấn phức tạp của hệ 3 hoạt chất Metronidazole, Chloramphenicol và Dexamethasone acetate.
Kết quả thực nghiệm trên thuốc viên nén Mycogynax đạt độ chính xác tương đương phương pháp HPLC tiêu chuẩn (độ thu hồi $98,80% - 99,32%$, sai số $< 1%$), đồng thời giảm hơn $75%$ thời gian phân tích và $85%$ chi phí hóa chất vận hành. Đây là giải pháp công nghệ mang giá trị thực tiễn cao, sẵn sàng chuyển giao và ứng dụng trực tiếp tại các nhà máy dược phẩm, trung tâm kiểm nghiệm trên toàn quốc, đóng góp thiết thực vào xu thế Hóa học xanh (Green Analytical Chemistry) và hiện đại hóa kiểm nghiệm dược phẩm.
