Trực khuẩn mủ xanh Pseudomonas aeruginosa là một trong những tác nhân gây nhiễm trùng bệnh viện nguy hiểm và thách thức nhất hiện nay. Trong thực hành lâm sàng, sự kết hợp giữa piperacillin và tazobactam là một trong những vũ khí chủ lực để đối phó với tác nhân này. Tuy nhiên, việc tối ưu hóa liều lượng để đảm bảo hiệu quả diệt khuẩn tối đa vẫn là một bài toán phức tạp đòi hỏi sự hỗ trợ của khoa học dữ liệu.
Vào tháng hai năm 2024, Cơ quan Quản lý Thực phẩm và Dược phẩm Hoa Kỳ FDA đã đưa ra một tuyên bố quan trọng liên quan đến điểm gãy nhạy cảm phụ thuộc liều lượng của trực khuẩn mủ xanh ở mức 16 microgram trên mililit. FDA khuyến cáo nên sử dụng liều 4.5 gram mỗi 6 giờ truyền kéo dài trong 3 giờ thay vì phương pháp truyền tĩnh mạch tiêu chuẩn trong 30 phút, vì phương pháp tiêu chuẩn không đạt được khả năng đạt mục tiêu điều trị trên 90 phần trăm. Bài viết này sẽ sử dụng ngôn ngữ lập trình R và thư viện mrgsolve để xây dựng mô hình dược động học quần thể, tiến hành mô phỏng và phân tích sâu về vấn đề này.
Xây dựng mô hình dược động học quần thể
Để mô phỏng chính xác nồng độ thuốc trong cơ thể, chúng ta cần một mô hình dược động học quần thể tối ưu. Chúng ta sẽ sử dụng một mô hình gộp phức tạp được xây dựng từ nhiều nhóm bệnh nhân khác nhau từ trẻ sơ sinh đến người cao tuổi. Dưới đây là mã nguồn R để thiết lập mô hình mrgsolve tập trung vào nồng độ piperacillin tự do không liên kết với protein, vốn là thành phần quyết định hiệu quả kháng khuẩn.
1library(mrgsolve)
2library(tidyverse)
3mod <- mcode(model = "piptazo_ken", code = '
4$PARAM
5theta_v1_pip = 10.4,
6tbw = 70,
7theta_cl_pip = 10.6,
8pma_year = 35,
9gamma1 = 3.35,
10gamma2 = 1.92,
11scr = 0.82,
12theta_scr = 0.346,
13dec50_pip = 89.1,
14theta_cl_sime = 1,
15theta_v2_sime = 1,
16funb_pip = 0.645,
17fdbs_pip = 1,
18theta_v2_pip = 11.6,
19theta_q2_pip = 15.2,
20mat50 = 54.2
21$CMT CENT_PIP PERI_PIP
22$MAIN
23double fsize = tbw / 70;
24double pma_week = pma_year * 52.17;
25double fmat = pow(pma_week,gamma1) / (pow(pma_week,gamma1) + pow(mat50,gamma1));
26double fdec_pip = 1 - (pow(pma_year, gamma2) / (pow(pma_year, gamma2) + pow(dec50_pip, gamma2)));
27double scrstd = exp(1.42 - (1.17 + 0.203 * log(pma_year / 100.0)) / sqrt(pma_year / 100.0));
28double fscr = exp(-theta_scr*(scr-scrstd));
29double v1_pip = theta_v1_pip * fsize * exp(ETA(1));
30double cl_pip = theta_cl_pip * pow(fsize, 0.75) * fmat * fdec_pip * fscr * theta_cl_sime * exp(ETA(2));
31double v2_pip = theta_v2_pip * fsize * theta_v2_sime * exp(ETA(4));
32double q2_pip = theta_q2_pip * pow(v2_pip/theta_v2_pip, 0.75) * exp(ETA(5));
33$OMEGA
340.1667645
350.1711123
360.1589037
370.547726
380.3579094
39$ODE
40dxdt_CENT_PIP = -(cl_pip/v1_pip)*CENT_PIP - (q2_pip/v1_pip)*CENT_PIP + (q2_pip/v2_pip)*PERI_PIP;
41dxdt_PERI_PIP = -(q2_pip/v2_pip)*PERI_PIP + (q2_pip/v1_pip)*CENT_PIP;
42$TABLE
43double free_pip = (CENT_PIP/v1_pip) * funb_pip;
44$CAPTURE
45free_pip')Mô hình này định nghĩa các tham số cố định và các biến số sinh lý của bệnh nhân như cân nặng, tuổi hiệu chỉnh và nồng độ creatinine huyết thanh. Khác với mô hình lý thuyết đầy đủ, chúng ta chỉ tập trung mô phỏng nồng độ piperacillin tự do vì đây là yếu tố chính tương tác trực tiếp với nồng độ ức chế tối thiểu MIC của vi khuẩn. Việc bỏ qua các sai số ngẫu nhiên dư thừa giúp mô hình gọn nhẹ và chạy nhanh hơn khi thực hiện mô phỏng quần thể lớn.
Phân tích khả năng đạt mục tiêu điều trị trên các quần thể
Chỉ số khả năng đạt mục tiêu điều trị PTA là xác suất để một chế độ liều cụ thể đạt được mục tiêu dược động học dược lực học mong muốn, thường là thời gian nồng độ thuốc tự do vượt quá MIC lớn hơn 50 phần trăm khoảng liều. Chúng ta sẽ tiến hành khảo sát chỉ số này trên nhiều nhóm đối tượng khác nhau để thấy rõ sự khác biệt.
Nhóm bệnh nhân không chăm sóc đặc biệt
Đầu tiên, chúng ta kiểm tra một bệnh nhân tiêu chuẩn nặng 90 kg, 35 tuổi với nồng độ creatinine huyết thanh là 0.83 miligram trên decilit.
Biểu đồ cho thấy chế độ truyền tĩnh mạch tiêu chuẩn trong 30 phút hoàn toàn thất bại trong việc đạt mức PTA 90 phần trăm khi MIC ở mức 16 microgram trên mililit. Chỉ khi áp dụng phương pháp truyền kéo dài trong 3 giờ, khả năng đạt mục tiêu mới vượt qua ngưỡng yêu cầu một cách an toàn.
Ảnh hưởng của tuổi và chức năng thận
Khi chúng ta điều chỉnh mô phỏng cho một quần thể có độ tuổi trung bình cao hơn là 50 tuổi và chức năng thận suy giảm nhẹ với creatinine huyết thanh là 1 miligram trên decilit, kết quả có sự thay đổi rõ rệt.
Nhờ độ thanh thải của thận giảm tự nhiên theo độ tuổi và chức năng thận, nồng độ thuốc trong máu được duy trì lâu hơn. Lúc này, phương pháp truyền ngắt quãng thông thường với liều 4 gram bắt đầu cho thấy kết quả khả quan hơn, tiệm cận được mục tiêu điều trị.
Nhóm bệnh nhân hồi sức cấp cứu nhiễm trùng huyết sớm
Đối với nhóm bệnh nhân nặng trong phòng hồi sức cấp cứu bị nhiễm trùng huyết sớm, các thông số sinh lý thay đổi rất phức tạp.
Kết quả mô phỏng cho thấy nồng độ tự do dường như cao hơn, một phần do giả định tỷ lệ thuốc tự do đạt 100 phần trăm trong một số mô hình nghiên cứu lâm sàng của nhóm bệnh nhân này. Tuy nhiên, trên thực tế lâm sàng, sự tăng lọc cầu thận và tăng thể tích phân bố ở bệnh nhân nhiễm trùng huyết có thể làm tăng tốc độ đào thải thuốc cực kỳ nhanh chóng.
Nhóm bệnh nhân sốt giảm bạch cầu
Đây là nhóm bệnh nhân có sự suy giảm khả năng đạt mục tiêu điều trị nghiêm trọng nhất trong các mô phỏng.
Ngay cả khi áp dụng chế độ truyền kéo dài 3 giờ, hầu hết các mức liều đều gặp khó khăn để vượt qua ngưỡng PTA 90 phần trăm tại mức MIC bằng 16, ngoại trừ liều cao 4 gram truyền kéo dài mỗi 6 giờ. Điều này giải thích tại sao việc điều trị nhiễm trùng trực khuẩn mủ xanh trên bệnh nhân sốt giảm bạch cầu luôn là một thách thức cực kỳ lớn đối với các bác sĩ lâm sàng.
Khảo sát các ngưỡng mục tiêu nghiêm ngặt hơn
Nhiều tài liệu y khoa gần đây cho rằng đích điều trị truyền thống là thời gian nồng độ thuốc tự do vượt quá MIC lớn hơn 50 phần trăm là chưa đủ đối với trực khuẩn mủ xanh, đặc biệt là ở những bệnh nhân nguy kịch. Các chuyên gia đề xuất nâng mục tiêu lên 100 phần trăm thời gian nồng độ thuốc tự do vượt quá MIC, hoặc thậm chí là 100 phần trăm thời gian nồng độ thuốc tự do vượt quá gấp 4 lần MIC.
Khi áp dụng các tiêu chuẩn khắt khe này vào mô hình mô phỏng, biểu đồ PTA sụt giảm rất mạnh. Chế độ truyền ngắt quãng tiêu chuẩn gần như hoàn toàn vô dụng, và ngay cả chế độ truyền kéo dài cũng đòi hỏi những mức liều tối đa mới có thể hy vọng đạt đích điều trị. Điều này phản ánh sự tương đồng lớn với các báo cáo lâm sàng thực tế, nhấn mạnh tầm quan trọng của việc cá thể hóa liều lượng dựa trên dược động học quần thể.
Vai trò thực sự của Tazobactam đối với trực khuẩn mủ xanh
Một câu hỏi thú vị được đặt ra là: Tại sao chúng ta cần thêm tazobactam vào phối hợp thuốc khi bản thân piperacillin đã có hoạt tính chống lại trực khuẩn mủ xanh? Phải chăng là do sự hiện diện của các enzyme beta lactamase kháng thuốc?
Để làm rõ điều này, chúng ta tiến hành truy vấn cơ sở dữ liệu bộ gen của Trung tâm Thông tin Công nghệ sinh học Quốc gia Hoa Kỳ NCBI đối với các chủng trực khuẩn mủ xanh nhạy cảm với piperacillin và tazobactam. Chúng ta tìm kiếm sự hiện diện của các gen mã hóa enzyme beta lactamase phổ biến.
Kết quả phân tích từ 124 chủng trực khuẩn mủ xanh có dữ liệu giải trình tự gen hoàn chỉnh cho thấy có tới 46 phần trăm chủng vi khuẩn mang gen blaOXA-2. Điều này có nghĩa là nếu không có sự hỗ trợ của chất ức chế beta lactamase là tazobactam để bảo vệ cấu trúc, piperacillin sẽ nhanh chóng bị thủy phân và mất hoàn toàn hoạt tính điều trị. Phát hiện này khẳng định vai trò không thể thay thế của tazobactam trong việc duy trì hiệu lực của kháng sinh phối hợp.
✨ Giá trị đắt giá: Việc tích hợp kỹ thuật mô phỏng dược động học quần thể bằng công cụ lập trình R giúp chúng ta không chỉ hiểu rõ hơn về mặt lý thuyết mà còn đưa ra được những quyết định lâm sàng chính xác. Sự kết hợp giữa khoa học dữ liệu và y học cá thể hóa chính là chìa khóa để tối ưu hóa liều lượng thuốc, cải thiện đáng kể tỷ lệ thành công trong điều trị các ca nhiễm trùng nặng do vi khuẩn đa kháng thuốc.
Câu hỏi tư duy hoặc bài tập ứng dụng
Dựa trên mô hình mrgsolve đã được thiết lập ở trên, bạn hãy thực hiện các bước sau:
1. Hãy viết mã nguồn R để thay đổi thông số độ thanh thải creatinine trong mô hình nhằm mô phỏng một bệnh nhân bị suy thận nặng với nồng độ creatinine huyết thanh tăng lên mức 2.5 miligram trên decilit.
2. Vẽ biểu đồ PTA tương ứng và đưa ra nhận xét: Chế độ truyền ngắt quãng 30 phút với liều 4 gram mỗi 6 giờ lúc này có đạt được mục tiêu PTA trên 90 phần trăm tại mức MIC bằng 16 hay không? Giải thích ý nghĩa lâm sàng của kết quả này đối với việc hiệu chỉnh liều trên bệnh nhân suy thận.
