Tổng quan nghiên cứu

Glycoside tim là nhóm dược chất steroid giữ vai trò quan trọng trong điều trị suy tim sung huyết và chống loạn nhịp tim với hơn 300 hợp chất tự nhiên được ghi nhận từ hơn 45 loài thực vật. Trong số các hoạt chất đó, digoxin và digitoxin là hai loại thuốc được chỉ định phổ biến nhất trong lâm sàng nhưng lại thuộc nhóm chất độc bảng A với cửa sổ điều trị cực kỳ hẹp. Hoạt chất digoxin có thời gian bán thải dài, tích lũy sâu trong các mô tim, gan, thận và dễ dẫn đến ngộ độc nghiêm trọng nếu hàm lượng đưa vào cơ thể vượt ngưỡng kiểm soát. Do đó, việc kiểm soát chính xác hàm lượng hoạt chất trong các chế phẩm viên nén là yêu cầu bắt buộc và cấp thiết trong công tác kiểm nghiệm dược phẩm.

Thách thức lớn nhất hiện nay là trong các công thức thuốc viên thương mại, tá dược độn như lactose, bột talc và magie stearat chiếm tới gần 98% tổng khối lượng, trong khi hàm lượng hoạt chất chỉ dao động từ 0,0% đến 0,3%. Các phương pháp phân tích tiêu chuẩn hiện hành như sắc ký lỏng hiệu năng cao (HPLC) đòi hỏi từ 15 đến 45 phút cho mỗi mẫu phân tích, sử dụng nhiều dung môi hữu cơ đắt tiền như axetonitrin và metanol, đồng thời quy trình chiết tách pha rắn phức tạp làm kéo dài thời gian kiểm nghiệm. Trong khi đó, phương pháp quang phổ hồng ngoại (FTIR) thông thường lại gặp hiện tượng xen phủ tín hiệu hấp thụ giữa tá dược và hoạt chất trong vùng dao động của khung steroid.

Nhằm giải quyết bài toán này, đề tài tập trung xây dựng quy trình phân tích định lượng đồng thời hai hoạt chất digoxin và digitoxin bằng phương pháp quang phổ hồng ngoại kết hợp kỹ thuật thống kê đa biến (Chemometrics). Nghiên cứu được thực hiện trong khuôn khổ chương trình hợp tác quốc tế giữa Việt Nam và Pháp, thực hiện thực nghiệm tại Bộ môn Hóa phân tích với hệ thống máy quang phổ Agilent Cary 600 FTIR và phần mềm tính toán MATLAB.

Ý nghĩa thực tiễn của công trình thể hiện rõ qua các chỉ số tối ưu hóa: cắt giảm thời gian phân tích từ 15 phút xuống dưới 3 phút cho mỗi phép đo, tiết kiệm 100% chi phí dung môi hữu cơ độc hại trong khâu xử lý mẫu, và nâng cao hiệu suất giám định chất lượng thuốc lưu hành tại các cơ sở kiểm nghiệm lên gấp 5 lầ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 sự giao thoa giữa hai trụ cột lý thuyết chính: lý thuyết quang phổ hấp thụ hồng ngoại phân tử và lý thuyết hóa học tính toán đa biến (Chemometrics).

Kỹ thuật quang phổ hồng ngoại biến đổi Fourier (FTIR) dựa trên định luật Lambert-Beer mở rộng, phản ánh dao động hóa trị và biến dạng của các liên kết hóa học cơ bản (C-H, C-C, C-O, O-H) và vòng lacton chưa bão hòa tại vị trí C-17 của khung steroid. Do digoxin và digitoxin có cấu trúc gần như đồng nhất (chỉ khác nhau một nhóm -OH tự do tại vị trí C-12), các dải hấp thụ trong vùng 3000 cm-1 đến 2800 cm-1 xảy ra hiện tượng chồng phổ sâu sắc.

Để giải quyết bài toán định lượng đa thành phần không cần tách chất, lý thuyết hồi quy đa biến tuyến tính được đưa vào ứng dụng. Trọng tâm của nghiên cứu là phương pháp Hồi quy cấu tử chính (Principal Component Regression - PCR), kết hợp giữa Phân tích thành phần chính (PCA) nhằm nén không gian dữ liệu phổ đa chiều và phương pháp Bình phương tối thiểu nghịch đảo (ILS) để xây dựng phương trình dự báo hàm lượng. Mô hình này được đối sánh trực tiếp với ba mô hình hồi quy khác gồm: Bình phương tối thiểu thông thường (CLS), Bình phương tối thiểu nghịch đảo thuần túy (ILS) và Bình phương tối thiểu từng phần (PLS).

Ba khái niệm cốt lõi được áp dụng xuyên suốt mô hình gồm: tính chất cộng tính của mật độ quang học hỗn hợp, ma trận điểm số và ma trận trọng số (Score và Loading matrices) giúp cô lập các biến ẩn, cùng phương sai tích lũy từ 1 đến 4 cấu tử chính (PC) nhằm loại trừ tối đa phương sai nhiễu ngẫu nhiên.

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

Nghiên cứu thiết kế một hệ thống thực nghiệm chuẩn mực với các tiêu chí cụ thể:

  • Cỡ mẫu và nguồn dữ liệu: Chuẩn bị 40 mẫu hỗn hợp chuẩn nhân tạo với tổng khối lượng 3,0 g mỗi mẫu, chia thành 30 mẫu cho tập chuẩn hóa ma trận (Calibration set) và 10 mẫu độc lập cho tập kiểm tra đánh giá sai số (Validation set). Dữ liệu đối chứng thực tế được thu thập trên 20 viên nén dược phẩm thương mại đang lưu hành trên thị trường.

  • Phương pháp chọn mẫu: Ma trận hàm lượng được thiết kế có chủ đích nhằm bao phủ dải nồng độ hoạt chất từ 0,0% đến 0,3% khối lượng, kết hợp tỷ lệ các tá dược điển hình gồm 60% đến 70% lactose, 23% đến 30% talc và 7% đến 9% magie stearat, mô phỏng chính xác cấu trúc nền mẫu thực tế của các nhà sản xuất.

  • Phương pháp phân tích và lý do lựa chọn: Phép đo quang phổ hồng ngoại truyền qua thu nhận đồng thời tín hiệu tại 106 điểm số sóng trong vùng 3000-2800 cm-1 trên thiết bị Agilent Cary 600 FTIR. Thuật toán PCR trên phần mềm MATLAB được lựa chọn làm phương pháp phân tích cốt lõi vì khả năng chuyển đổi ma trận tín hiệu 30x106 thành các trục trực giao mới, cô lập hoàn toàn các sai số ngẫu nhiên do lực nén ép viên và hiện tượng tán xạ quang học trên nền chất mang KBr.

  • Timeline nghiên cứu: Quá trình thực nghiệm được triển khai có hệ thống qua 4 giai đoạn trong 12 tháng: tối ưu hóa thông số ép viên KBr, thu thập ma trận phổ 106 bước sóng, lập trình thuật toán xử lý dữ liệu trên MATLAB, và kiểm định chéo với sắc ký lỏng hiệu năng cao HPLC Shimadzu sử dụng detector PDA tại bước sóng 225 nm với tốc độ pha động 0,7 ml/phút ở nhiệt độ cột 30 độ C.

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

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

Quá trình nghiên cứu thực nghiệm đã ghi nhận 4 phát hiện quan trọng có giá trị khoa học và ứng dụng cao:

  1. Tối ưu hóa điều kiện đo mẫu rắn trên nền KBr: Nghiên cứu xác định tỷ lệ khối lượng phối trộn tối ưu giữa mẫu và bột KBr là 3:97 (tương ứng 3 mg hỗn hợp dược chất nghiền cùng 97 mg KBr). Với lượng mẫu ép viên 25 mg và lực nén 15 lần, mật độ quang hấp thụ tại hai số sóng đặc trưng 2900 cm-1 và 2932 cm-1 đạt giá trị ổn định lần lượt là 1,1020 và 0,8421. Mức độ lặp lại của phép đo đạt độ lệch chuẩn tương đối RSD chỉ từ 0,02% đến 0,03% khi nhiệt độ môi trường dao động từ 25 đến 35 độ C.

  2. Xác định số cấu tử chính tối ưu cho mô hình PCR: Thuật toán PCR áp dụng phép phân rã giá trị số ít (SVD) trên ma trận tín hiệu 30x106 đã chứng minh rằng khi chọn số cấu tử chính từ 3 đến 4 (PC = 3-4), mô hình giải thích được hơn 99% tổng phương sai của toàn bộ tập dữ liệu phổ. Mô hình bóc tách thành công tín hiệu riêng biệt của digoxin và digitoxin bất chấp sự tương đồng cao về công thức phân tử.

  3. Độ nhạy và độ tin cậy của mô hình định lượng: Giới hạn phát hiện (LOD) và giới hạn định lượng (LOQ) của phương pháp PCR đạt mức microgram, đáp ứng hoàn hảo tiêu chuẩn kiểm nghiệm vi lượng. Đánh giá trên mẫu thêm chuẩn cho thấy hiệu suất thu hồi của phương pháp đạt dải từ 98,2% đến 101,5%, với độ lệch chuẩn tương đối RSD luôn duy trì dưới 2,5%.

  4. Kiểm chứng trên mẫu dược phẩm thương mại: Khi phân tích các mẫu thuốc lưu hành trên thị trường, phương pháp Chemometrics-FTIR cho kết quả hàm lượng hoạt chất tương đương với phép đo sắc ký HPLC chuẩn (sai số tương đối dưới 3,5%), nhưng rút ngắn thời gian phân tích hơn 80% (từ 15 phút của HPLC xuống dưới 3 phút cho một mẫu).

Thảo luận kết quả

Hiện tượng xen phủ phổ nghiêm trọng giữa tá dược (chiếm 98% khối lượng) và hai glycoside tim đã khiến các phương pháp quang phổ đơn biến truyền thống hoàn toàn bất lực. Sự khác biệt duy nhất giữa digoxin (C41H64O14) và digitoxin (C41H64O13) là một nhóm -OH tại vị trí C-12, khiến các đỉnh hấp thụ tại vùng 3600-2800 cm-1 gần như hòa lẫn vào nhau.

Khi áp dụng phương pháp hồi quy cấu tử chính PCR, các biến độc lập tương quan chặt chẽ được chuyển đổi thành các biến trực giao không tương quan. Nhờ đó, các biến động ngẫu nhiên về độ dày viên nén, kích thước hạt và độ hấp thụ của KBr được chuyển vào sai số dư và loại bỏ khỏi phương trình hồi quy. Ngược lại, mô hình CLS cho thấy sai số lớn vì phụ thuộc vào độ hấp thụ tuyệt đối của phổ thô, còn mô hình PLS và ILS gặp khó khăn khi chọn lọc số bước sóng đặc trưng cho hệ 3 cấu tử tá dược phức tạp.

Dữ liệu phân tích trong nghiên cứu có thể được trực quan hóa hiệu quả thông qua biểu đồ phân bố điểm số không gian ba chiều (3D Score Plot) thể hiện sự phân cụm rõ nét của các nhóm mẫu chuẩn và mẫu kiểm tra. Đồng thời, bảng tương quan tuyến tính giữa nồng độ thực tế và nồng độ dự báo từ mô hình PCR phản ánh hệ số tương quan R2 đạt mức 0,9985, khẳng định độ chuẩn xác tuyệt đối của phương pháp mới. So với quy trình HPLC đòi hỏi tiêu tốn nhiều dung môi axetonitrin, phương pháp FTIR-PCR hoàn toàn là một giải pháp xanh và có tính thời sự cao.

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

Nhằm đưa kết quả nghiên cứu vào ứng dụng thực tiễn trong ngành dược phẩm, các giải pháp cụ thể được đề xuất như sau:

  1. Chuẩn hóa quy trình chế bị mẫu ép viên tự động: Các nhà máy dược phẩm và trung tâm kiểm nghiệm cần đầu tư hệ thống dập viên thủy lực tự động với lực nén cố định 15 chu kỳ và cân phân tích có độ chính xác 0,0001 g. Mục tiêu là đưa độ lệch chuẩn tương đối RSD của thao tác chuẩn bị mẫu xuống dưới 1,5% trong vòng 6 tháng tới, do bộ phận đảm bảo chất lượng (QA) chủ trì thực hiện.

  2. Tích hợp mô-đun phần mềm PCR vào thiết bị phân tích: Doanh nghiệp sản xuất cần triển khai đóng gói các thuật toán hồi quy PCR viết bằng ngôn ngữ MATLAB thành mô-đun phần mềm tích hợp trực tiếp vào giao diện điều khiển máy quang phổ FTIR. Thời gian hoàn thành dự kiến trong 3 tháng, giúp kỹ thuật viên vận hành tự động trả kết quả hàm lượng trong vòng dưới 2 phút mỗi mẫu.

  3. Mở rộng cơ sở dữ liệu phổ cho các glycoside tim khác: Các viện nghiên cứu và trường đại học chuyên ngành dược cần phối hợp xây dựng ma trận phổ chuẩn cho ít nhất 5 hoạt chất glycoside tim tự nhiên khác như lanatoside C, ouabain và acetyldigoxin. Lộ trình thực hiện trong 12 tháng nhằm phục vụ công tác kiểm tra nguyên liệu dược liệu thô nhập khẩu.

  4. Ban hành tiêu chuẩn cơ sở và hướng dẫn thẩm định đa biến: Cục Quản lý Dược và Hội đồng Dược điển cần ban hành văn bản hướng dẫn kỹ thuật công nhận phương pháp Chemometrics-FTIR tương đương với tiêu chuẩn HPLC trong Dược điển Việt Nam. Lộ trình triển khai trong 24 tháng tới nhằm giúp hệ thống kiểm nghiệm toàn quốc tiết kiệm ước tính 40% chi phí hóa chất hàng năm.

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

Luận văn mang lại giá trị thực tiễn và học thuật to lớn cho 4 nhóm đối tượng trọng tâm:

  1. Chuyên viên kiểm nghiệm chất lượng (QC/QA) tại các nhà máy dược phẩm: Nhận được quy trình định lượng nhanh không phá hủy mẫu, áp dụng trực tiếp để kiểm tra độ đồng đều hàm lượng các lô thuốc viên nén digoxin lưu hành với tốc độ nhanh gấp 5 lần so với chạy máy HPLC.

  2. Giảng viên và nhà nghiên cứu chuyên ngành Hóa phân tích, Hóa dược: Tiếp cận tài liệu tham khảo chuyên sâu về phương pháp ứng dụng Chemometrics và xử lý tín hiệu phổ đa chiều trên phần mềm MATLAB, làm cơ sở xây dựng giáo trình giảng dạy sau đại học.

  3. Thanh tra viên và cơ quan quản lý dược phẩm: Nắm bắt công cụ sàng lọc nhanh tại hiện trường giúp giám định thuốc thật, thuốc giả, thuốc kém chất lượng chứa hoạt chất độc bảng A mà không cần vận hành các hệ thống máy sắc ký cồng kềnh.

  4. Học viên cao học và sinh viên ngành Công nghệ Dược: Khảo cứu phương pháp thiết kế ma trận thực nghiệm 30x106, quy trình xử lý dữ liệu phổ toàn phần và kỹ thuật thẩm định phương pháp phân tích phục vụ việc thực hiện các đề tài khóa luận tốt nghiệp.

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

  1. Tại sao không thể sử dụng quang phổ hồng ngoại truyền thống để định lượng digoxin và digitoxin? Cấu trúc phân tử của digoxin và digitoxin gần như giống hệt nhau, chỉ khác biệt duy nhất một nhóm hydroxyl ở vị trí C-12, khiến các đỉnh hấp thụ trong vùng 3600-2800 cm-1 trùng lặp hoàn toàn. Ngoài ra, ba loại tá dược chiếm tới 98% khối lượng viên tạo ra hiện tượng chồng phổ phức tạp, làm cho định luật Lambert-Beer đơn biến không thể phân tách riêng từng hoạt chất.

  2. Thuật toán PCR vượt trội hơn các mô hình PLS và CLS ở những điểm nào? Thuật toán PCR thực hiện phân tích thành phần chính để nén dữ liệu từ 106 bước sóng thành 3 đến 4 cấu tử chính, giúp cô lập toàn bộ nhiễu nền KBr và sai số lực nén vào các thành phần phương sai thấp. Trong khi đó, mô hình CLS bị sai số lớn do phụ thuộc vào độ hấp thụ tuyệt đối, còn PLS dễ bị mất ổn định trước các hệ mẫu có độ tán xạ cao.

  3. Quy trình chuẩn bị mẫu ép viên KBr có những yêu cầu kỹ thuật nghiêm ngặt nào? Kỹ thuật đòi hỏi tỷ lệ trộn chuẩn xác 3 mg mẫu với 97 mg KBr (tỷ lệ 3:97), nghiền mịn đồng nhất trong 10 phút và nén ép 15 lần với lượng bột 25 mg để viên đạt độ trong suốt tối đa. Quá trình thao tác bắt buộc thực hiện trong phòng có độ ẩm dưới 50% nhằm triệt tiêu hoàn toàn đỉnh hấp thụ giả của nước tại vùng 3400 cm-1.

  4. Độ chính xác của phương pháp Chemometrics-FTIR so với HPLC ra sao? Thực nghiệm đối chứng trên 10 mẫu kiểm tra cho thấy phương pháp FTIR-PCR đạt hiệu suất thu hồi từ 98,2% đến 101,5% với độ lệch chuẩn tương đối RSD dưới 2,5%, hoàn toàn tương đương với sắc ký HPLC chuẩn (RSD dưới 2,0%). Tuy nhiên, phương pháp FTIR giúp tiết kiệm 100% dung môi độc hại và rút ngắn thời gian phân tích xuống dưới 3 phút.

  5. Phương pháp này có thể áp dụng cho các dạng bào chế khác ngoài viên nén không? Phương pháp hoàn toàn có thể mở rộng cho dung dịch tiêm, cồn thuốc hoặc cao chiết dược liệu bằng cách áp dụng kỹ thuật đo màng mỏng chất lỏng hoặc đầu đo phản xạ toàn phần suy giảm (ATR-FTIR). Khi thay đổi dạng bào chế, người nghiên cứu chỉ cần chuẩn hóa lại ma trận nồng độ và lựa chọn số cấu tử chính phù hợp trong MATLAB.

Kết luận

  • Xây dựng thành công quy trình định lượng đồng thời hai hoạt chất độc bảng A là digoxin và digitoxin trong dải số sóng 3000-2800 cm-1.
  • Tối ưu hóa điều kiện chuẩn bị mẫu rắn ép viên với tỷ lệ mẫu/KBr đạt chuẩn 3:97 và khối lượng ép 25 mg, cho độ ổn định cao với RSD dưới 0,03%.
  • Ứng dụng xuất sắc thuật toán PCR trên phần mềm MATLAB với 3-4 cấu tử chính, giải thích trên 99% phương sai dữ liệu và loại bỏ hoàn toàn nhiễu nền tá dược.
  • Thẩm định phương pháp đạt hiệu suất thu hồi 98,2% đến 101,5%, độ lệch chuẩn RSD dưới 2,5%, cho kết quả tương đương phương pháp HPLC tiêu chuẩn.
  • Cắt giảm hơn 80% thời gian phân tích và loại bỏ 100% lượng dung môi hữu cơ độc hại, tạo ra bước tiến quan trọng cho xu hướng phân tích xanh.

Đóng góp lớn nhất của luận văn là đã chuyển giao thành công kỹ thuật Chemometrics từ nghiên cứu lý thuyết sang ứng dụng thực tiễn trong kiểm nghiệm dược phẩm, giải quyết triệt để bài toán phân tích nhanh các thuốc có phạm vi điều trị hẹp.

Trong 6 tháng tới, nhóm nghiên cứu sẽ tiếp tục mở rộng thử nghiệm trên các dạng bào chế viên nang và hoàn thiện phần mềm thương mại hóa.

Quý độc giả, các nhà khoa học và các đơn vị kiểm định chất lượng hãy tham khảo toàn văn luận văn thạc sĩ và liên hệ ngay với nhóm tác giả để tiếp nhận chuyển giao quy trình phân tích hiện đại này.