Bối cảnh và vấn đề nghiên cứu

Bệnh lao là một bệnh truyền nhiễm qua đường hô hấp do trực khuẩn Mycobacterium tuberculosis gây ra và là một trong những nguyên nhân gây tử vong hàng đầu trên thế giới do tác nhân truyền nhiễm. Theo số liệu của Tổ chức Y tế Thế giới (WHO), vào năm 2020 ước tính toàn cầu có khoảng 9,9 triệu người mắc lao mới và khoảng 1,5 triệu người tử vong (trong đó có 1,3 triệu người âm tính với HIV và 214.000 người dương tính với HIV). Việt Nam là quốc gia có gánh nặng bệnh lao cao, đứng thứ 11 trong số 30 quốc gia có tỷ lệ mắc lao cao nhất toàn cầu và đồng thời đứng thứ 11 trong số 30 quốc gia có gánh nặng lao đa kháng thuốc. Mặc dù các mục tiêu thuộc chiến lược chấm dứt bệnh lao (End TB Strategy) của WHO giai đoạn 2020–2035 đặt ra yêu cầu giảm 90% tỷ lệ mắc mới và giảm 95% tỷ lệ tử vong, công tác điều trị vẫn gặp nhiều rào cản do phác đồ kéo dài, độc tính lớn, chi phí cao và nguy cơ gia tăng của vi khuẩn kháng thuốc.

Ethambutol (viết tắt là E hoặc EMB) là thuốc chống lao tổng hợp được sử dụng từ năm 1960 và là thành phần cốt lõi trong phác đồ 4 thuốc hàng 1 điều trị lao nhạy cảm (2RHZE/4RHE). Tuy nhiên, thuốc có khoảng điều trị hẹp và gây độc tính nghiêm trọng trên thị giác (viêm dây thần kinh thị giác), đồng thời sự hấp thu và chuyển hóa thuốc có biến thiên lớn giữa các cá thể. Tại Việt Nam, trước thời điểm đề tài thực hiện, chưa có nghiên cứu nào công bố về dược động học quần thể (Population Pharmacokinetics - PopPK) cũng như phân tích dược động học/dược lực học (PK/PD) của ethambutol trên bệnh nhân lao phổi.

Khóa luận được thực hiện nhằm giải quyết các mục tiêu nghiên cứu cụ thể sau:

  1. Xây dựng mô hình dược động học quần thể của Ethambutol ở bệnh nhân lao phổi.
  2. Mô phỏng Dược động học/Dược lực học (PK/PD) của Ethambutol trên bệnh nhân lao phổi mới và lao tái trị.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: 136 bệnh nhân mắc lao phổi mới hoặc lao phổi tái trị (sau thất bại điều trị hoặc tái phát) đang được điều trị bằng phác đồ chứa các thuốc chống lao thiết yếu hàng 1. Dữ liệu lâm sàng và nồng độ thuốc được trích xuất từ nhiệm vụ khoa học và công nghệ cấp nhà nước (mã số HNQT/SPĐP/01.16).
  • Không gian nghiên cứu: Nghiên cứu thu thập dữ liệu tại 3 cơ sở y tế chuyên khoa: Bệnh viện Phổi Hà Nội, Bệnh viện Phổi Trung ương và Bệnh viện 74 Trung ương.
  • Thời gian nghiên cứu: Dữ liệu bệnh nhân được thu nhận trong giai đoạn từ tháng 3/2017 đến tháng 3/2018; khóa luận hoàn thành xử lý dữ liệu và nghiệm thu vào năm 2022.

Cơ sở lý thuyết và phương pháp

Khóa luận áp dụng phương pháp nghiên cứu định lượng kết hợp mô hình hóa toán sinh học và mô phỏng lâm sàng.

Khung lý thuyết và công thức ứng dụng

  • Mô hình hóa ảnh hưởng hỗn hợp phi tuyến (Non-linear Mixed Effects Modeling - NLMEM): Phân tích đồng thời các thông số cố định (fixed effects - đại diện cho giá trị trung bình quần thể như độ thanh thải, thể tích phân bố) và thông số ngẫu nhiên (random effects - biến thiên giữa các cá thể $\omega$ và sai số thặng dư $\varepsilon$).
  • Đích Dược động học/Dược lực học (PK/PD): Hai chỉ số hiệu quả được áp dụng gồm $f\text{AUC}{0-24}/\text{MIC} > 11{,}9$ và $fC{\max}/\text{MIC} > 0{,}48$ dựa trên các nghiên cứu in vitro và lâm sàng quốc tế. Ngưỡng xác suất đạt đích (Probability of Target Attainment - PTA) tối ưu được quy chuẩn ở mức $\ge 90%$. Dải nồng độ ức chế tối thiểu (MIC) được khảo sát là 0,5; 1; 2; 4 và 8 mg/L.
  • Hiệu chỉnh kích thước cơ thể:
    • Cân nặng trừ mỡ (Fat-Free Mass - FFM) tính theo công thức Janmahasatian (2005):
      • Nam giới: $\text{FFM} = \frac{9{,}27 \times 10^3 \times \text{BWt}}{6{,}68 \times 10^3 + 216 \times \text{BMI}}$
      • Nữ giới: $\text{FFM} = \frac{9{,}27 \times 10^3 \times \text{BWt}}{8{,}78 \times 10^3 + 244 \times \text{BMI}}$
    • Tốc độ lọc cầu thận ước tính (eGFR) tính theo công thức MDRD: $$\text{eGFR} = 175 \times (\text{Scr})^{-1{,}154} \times (\text{Tuổi})^{-0{,}203} \times (0{,}742\text{ nếu là nữ})$$

Phương pháp nghiên cứu và công cụ xử lý

  • Quy trình thu thập mẫu sinh học: Nồng độ ethambutol trong huyết tương được lấy theo 2 chế độ lấy mẫu thưa (sparse sampling):
    • Chế độ 1: Lấy mẫu tại các thời điểm trước khi dùng thuốc và sau khi uống 1, 4, 6 giờ.
    • Chế độ 2: Lấy mẫu tại các thời điểm trước khi dùng thuốc và sau khi uống 2, 5, 8 giờ.
  • Quy trình xây dựng mô hình PopPK: Thực hiện trên phần mềm MONOLIX phiên bản 2021R1 qua các bước:
    1. Khảo sát mô hình cấu trúc: Kiểm tra mô hình 1 ngăn và 2 ngăn kết hợp các đặc tính hấp thu (bậc 1, thời gian trễ $T_{\text{lag}}$, mô hình ngăn chuyển tiếp - transit compartment, hấp thu tuần tự bậc 0 và bậc 1).
    2. Mô hình sai số thặng dư: Đánh giá sai số cộng, sai số tỉ lệ, sai số hỗn hợp cộng - tỉ lệ, sai số mũ, sai số log-normal.
    3. Sàng lọc yếu tố ảnh hưởng (Covariates): Áp dụng quy trình từng bước gồm thêm lần lượt (forward inclusion, tiêu chuẩn $\Delta -2\text{LL} \ge 3{,}841$, $p < 0{,}05$) và loại trừ lần lượt (backward elimination, tiêu chuẩn $\Delta -2\text{LL} \ge 6{,}635$, $p < 0{,}01$).
    4. Thẩm định mô hình: Sử dụng đồ thị tương quan giữa nồng độ quan sát và dự đoán, đồ thị sai số dự đoán có trọng số cá thể (IWRES), đồ thị kiểm tra chẩn đoán trực quan (Visual Predictive Checks - VPC) dựa trên 10.000 bộ dữ liệu mô phỏng.
  • Mô phỏng PK/PD: Thực hiện bằng mô phỏng Monte Carlo qua gói phần mềm mlxR trên 1.000 bệnh nhân ảo để đánh giá PTA theo cân nặng (25–39 kg, 40–54 kg, $\ge 55\text{ kg}$), theo giới tính và theo các mức liều giả định tính theo FFM (15, 20, 25, 30, 35 mg/kg FFM).
  • Thống kê y sinh: Sử dụng phần mềm SPSS và Microsoft Excel 365. Các biến định lượng được mô tả bằng trung vị (khoảng tứ phân vị) hoặc trung bình $\pm$ độ lệch chuẩn; biến định tính trình bày dưới dạng tần số và tỷ lệ phần trăm.
  • Y đức: Nghiên cứu nằm trong đề tài cấp nhà nước HNQT/SPĐP/01.16 đã được Hội đồng Đạo đức trong nghiên cứu Y sinh học – Khoa Y Dược, Đại học Quốc gia Hà Nội phê duyệt và có văn bản báo cáo Ban đánh giá vấn đề đạo đức của Bộ Y tế.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 hệ thống hóa cơ sở lý luận về bệnh học, dược lý học và phương pháp luận PopPK:

  • Bệnh học và dịch tễ: Trực khuẩn lao M. tuberculosis lây qua các hạt khí dung khi bệnh nhân ho, nói hoặc hắt hơi. Dòng vi khuẩn Bắc Kinh (Beijing genotype) chiếm ưu thế tại miền Bắc Việt Nam, có liên quan chặt chẽ đến khả năng lây lan và tỷ lệ kháng thuốc cao. Tỷ lệ mắc lao ở nam giới tại Việt Nam cao gấp 5,1 lần so với nữ giới.
  • Dược lý học ethambutol: Cấu trúc $(2S,2'S)\text{-}2,2'\text{-(ethylenediimino)di-1-butanol}$, khối lượng phân tử phân bố dưới dạng muối hydroclorid dễ tan trong nước. Cơ chế tác dụng chính là ức chế enzym arabinosyltransferase (mã hóa bởi cụm gen embCAB), ngăn chặn tổng hợp arabinogalactan và lipoarabinomannan của thành tế bào vi khuẩn. Đột biến điểm tại codon embB 306 là cơ chế kháng thuốc phổ biến nhất.
  • Đặc tính Dược động học: Thuốc hấp thu nhanh qua đường tiêu hóa (75–80%), thể tích phân bố lớn ($V_d \approx 1{,}6\text{ L/kg}$), chuyển hóa qua gan khoảng 15%, thải trừ chủ yếu qua thận (đến 80% trong 24 giờ, trong đó 50% ở dạng nguyên vẹn qua lọc cầu thận và bài tiết tích cực tại ống thận; độ thanh thải thận khoảng 7 ml/phút/kg).
  • Tổng quan nghiên cứu quốc tế: Tác giả phân tích 6 công trình PopPK tiêu biểu trên thế giới (chủ yếu thực hiện tại Châu Phi như Nam Phi, Tanzania, Rwanda, Ghana). Hầu hết các nghiên cứu đều chỉ ra mô hình 2 ngăn là mô hình cấu trúc tối ưu; độ thanh thải ($CL$) và thể tích phân bố ($V$) có liên quan mật thiết đến kích thước cơ thể (cân nặng toàn phần, FFM).

Chương 2: Đối tượng và phương pháp nghiên cứu

Trình bày chi tiết dữ liệu 136 bệnh nhân lao phổi điều trị thuốc hàng 1 tại 3 bệnh viện từ tháng 3/2017 đến tháng 3/2018. Tác giả thiết lập quy trình nghiên cứu 5 bước rõ ràng:

  1. Chuẩn hóa dữ liệu nồng độ – thời gian và các yếu tố ảnh hưởng lâm sàng sang định dạng chuẩn của phần mềm MONOLIX.
  2. Thiết lập mô hình PopPK cơ bản (khảo sát số ngăn, kiểu hấp thu, phương trình sai số).
  3. Đánh giá tính cộng tuyến giữa các biến lâm sàng, đưa biến vào mô hình bằng thuật toán step-wise (forward/backward) dựa trên tiêu chí Likelihood Ratio Test.
  4. Thẩm định mô hình qua các biểu đồ chẩn đoán toán học và đồ thị VPC.
  5. Mô phỏng Monte Carlo trên mlxR để tính toán chỉ số PTA ở các ngưỡng MIC khác nhau.

Chương 3: Kết quả nghiên cứu

Chương này cung cấp toàn bộ các số liệu thực nghiệm thu được từ quá trình phân tích dữ liệu và mô hình hóa:

  • Đặc điểm mẫu nghiên cứu ($N = 136$):
    • Tuổi trung vị: 45,5 (tứ phân vị: 33 – 58).
    • Giới tính: Nam giới chiếm ưu thế với 73,53% ($n = 100$).
    • Cân nặng trung vị: 49 kg (44 – 55); BMI trung vị: 18,69 $\text{kg/m}^2$; FFM trung vị: 41,29 kg (35,07 – 46,2).
    • Phân loại bệnh: Lao mới 81 bệnh nhân (59,56%), lao tái trị 55 bệnh nhân (40,44%). Bệnh nhân đồng nhiễm HIV chỉ có 1 trường hợp (0,74%); đái tháo đường 21 trường hợp (15,44%). Độ thanh thải creatinin ước tính theo MDRD đạt trung vị 93,35 $\text{mL/phút/1{,}73 m}^2$.
  • Thực trạng phân liều ethambutol:
    • Liều dùng phổ biến nhất là 800 mg/ngày (72,06%) và 1200 mg/ngày (25,74%).
    • Liều tính theo cân nặng thực tế có trung vị là 18,18 mg/kg (dao động rộng từ 10,53 đến 25,00 mg/kg). Liều tính theo FFM đạt trung vị 22,57 mg/kg FFM (14,24 đến 33,00 mg/kg FFM).
    • Xuất hiện tình trạng lệch liều theo khoảng cân nặng: nhóm bệnh nhân 25–39 kg và $>70\text{ kg}$ có mức liều thực tế theo mg/kg cao hơn ngưỡng khuyến cáo; nhóm 40–54 kg dùng liều sát mức sàn khuyến cáo (trung bình 17,29 mg/kg); có bệnh nhân dùng mức liều chỉ 10,53 mg/kg (dưới mức sàn 15 mg/kg của WHO) và có bệnh nhân dùng tới 25 mg/kg (vượt mức trần 20 mg/kg).
  • Mô hình PopPK cấu trúc và sai số:
    • Khảo sát 10 mô hình cấu trúc cho thấy mô hình 2 ngăn kết hợp quá trình hấp thu tuần tự bậc 0 và bậc 1, có thời gian trễ trước giai đoạn hấp thu bậc 0 (kí hiệu MH 2.5) cho độ khớp tốt nhất với $-2\text{LL} = 650{,}09$; $\text{AIC} = 672{,}09$; $\text{BIC} = 704{,}19$; $\text{BICc} = 676{,}31$.
    • Mô hình sai số tỉ lệ ($C_{\text{obs}} = C_{\text{pred}} + b \times C_{\text{pred}} \times \varepsilon$) cho giá trị tiêu chí thông tin thấp nhất ($-2\text{LL} = 603{,}64$).
    • Khảo sát hiệp phương sai ngẫu nhiên ghi nhận tương quan có ý nghĩa giữa độ thanh thải $Cl$ và thanh thải liên ngăn $Q$ ($\text{corr_Cl_Q} = 0{,}714$), làm giảm $-2\text{LL}$ tới 45,87 đơn vị.
  • Mô hình có yếu tố dự đoán (Covariate Model):
    • Quá trình forward inclusion xác định FFM làm giảm $-2\text{LL}$ mạnh nhất khi liên kết với $Cl$ ($\Delta -2\text{LL} = -6{,}8$). Tiếp tục bổ sung biến Tuổi vào $Cl$ làm giảm thêm 6,76 đơn vị $-2\text{LL}$.
    • Quá trình backward elimination xác nhận khi loại bỏ $Cl\text{-}FFM$, $-2\text{LL}$ tăng 10,79 đơn vị; khi loại $Cl\text{-Tuổi}$, $-2\text{LL}$ tăng 9,87 đơn vị (đều vượt ngưỡng cắt 6,635 ứng với $p < 0{,}01$).
    • Phương trình hồi quy cuối cùng biểu diễn độ thanh thải cá thể: $$\log(Cl) = \log(Cl_{\text{pop}}) + 0{,}346 \times \log\left(\frac{\text{FFM}}{41{,}29}\right) - 0{,}00316 \times \log\left(\frac{\text{Tuổi}}{45{,}5}\right) + \varepsilon_{Cl}$$ Trong đó giá trị độ thanh thải quần thể chuẩn $Cl_{\text{pop}} = 45{,}8\text{ L/h}$ đối với người có FFM = 41,29 kg và 45,5 tuổi.
  • Kết quả mô phỏng PK/PD:
    • Mô phỏng theo nhóm cân nặng: Tại $\text{MIC} = 0{,}5\text{ mg/L}$, PTA của $f\text{AUC}/\text{MIC} > 11{,}9$ đạt tối ưu ($>98%$). Tại $\text{MIC} = 1\text{ mg/L}$, tỷ lệ đạt đích giảm xuống mức trên 50%. Tại $\text{MIC} = 2\text{ mg/L}$, PTA giảm nghiêm trọng xuống dưới 10%; tại $\text{MIC} \ge 4\text{ mg/L}$, không còn bệnh nhân nào đạt đích. Đối với đích $fC_{\max}/\text{MIC} > 0{,}48$, 100% bệnh nhân đạt đích tại $\text{MIC} \le 1\text{ mg/L}$; $>98%$ tại $\text{MIC} = 2\text{ mg/L}$; $>85%$ tại $\text{MIC} = 4\text{ mg/L}$ và giảm xuống tối đa 30,4% tại $\text{MIC} = 8\text{ mg/L}$.
    • Mô phỏng theo giới tính: Bệnh nhân nữ có tỷ lệ đạt cả hai đích PK/PD cao hơn nam giới ở các dải MIC cao. Cụ thể tại $\text{MIC} = 8\text{ mg/L}$, PTA của $fC_{\max}/\text{MIC}$ ở nữ đạt 54,2% so với 34,8% ở nam.
    • Mô phỏng theo liều FFM: Tăng liều theo FFM làm tăng rõ rệt PTA. Ở $\text{MIC} = 4\text{ mg/L}$, PTA của $fC_{\max}/\text{MIC}$ tăng từ 55,3% (ở liều 15 mg/kg FFM) lên 97,9% (ở liều 35 mg/kg FFM). Tại $\text{MIC} = 8\text{ mg/L}$, tỷ lệ đạt đích tăng từ 3,5% (liều 15 mg/kg FFM) lên 68,5% (liều 35 mg/kg FFM).

Chương 4: Bàn luận

Tác giả đối chiếu phát hiện của khóa luận với y văn thế giới:

  • Đặc điểm dịch tễ và liều dùng: Tỷ lệ nam giới chiếm 73,53% phù hợp với các nghiên cứu dịch tễ lao tại Việt Nam (nam giới thường có hành vi nguy cơ cao hơn như hút thuốc, uống rượu). Việc phân chia liều dùng theo từng bậc cân nặng cố định (400, 800, 1200 mg) dẫn đến độ phân tán liều thực tế theo mg/kg quá rộng, làm phát sinh nguy cơ thất bại điều trị ở người bị dưới liều hoặc ngộ độc thị giác ở người thừa liều.
  • Cơ chế mô hình cấu trúc: Quá trình hấp thu phức tạp tuần tự bậc 0 và bậc 1 có thời gian trễ được lý giải do dạng bào chế viên nén cần thời gian rã và hòa tan trong dạ dày trước khi thuốc được vận chuyển xuống ruột non và hấp thu nhanh theo cơ chế thụ động.
  • Ảnh hưởng của FFM và Tuổi: FFM là đại lượng đại diện chính xác hơn tổng cân nặng (WT) vì ethambutol là thuốc thân nước, phân bố chủ yếu vào khối cơ thể không chứa mỡ. Độ thanh thải $Cl$ giảm nhẹ theo tuổi tác, phản ánh sự suy giảm sinh lý chức năng lọc cầu thận và bài tiết qua ống thận ở người lớn tuổi.
  • Ý nghĩa lâm sàng của mô phỏng: Phác đồ liều hiện hành chỉ đảm bảo hiệu quả điều trị tối ưu khi vi khuẩn có mức $\text{MIC} \le 0{,}5\text{ mg/L}$. Đối với các chủng có $\text{MIC} \ge 1\text{ mg/L}$ (vẫn nằm trong khoảng nhạy cảm theo lâm sàng), việc áp dụng liều thông thường tiềm ẩn nguy cơ thất bại điều trị rất cao, chứng minh sự cần thiết phải cá thể hóa liều dựa trên cân nặng trừ mỡ (FFM).

Kết quả và đóng góp

Bảng tổng hợp thông số mô hình PopPK cuối cùng của Ethambutol

Dữ liệu ước tính các thông số dược động học quần thể từ phần mềm MONOLIX được tổng hợp chi tiết trong bảng sau:

Thông số Ý nghĩa dược động học Giá trị quần thể Sai số chuẩn tương đối (RSE %) Biến thiên cá thể ($\omega$) RSE % của $\omega$
Tk01 (h) Khoảng thời gian hấp thu bậc 0 1,81 0,103% - -
ka2 ($\text{h}^{-1}$) Hằng số tốc độ hấp thu bậc 1 0,308 4,06% - -
F1 Tỷ lệ thuốc hấp thu qua pha bậc 0 0,218 7,52% 0,651 10,8%
Tlag1 (h) Thời gian trễ trước khi hấp thu bậc 0 0,051 48,3% 1,68 14,7%
Cl (L/h) Độ thanh thải từ ngăn trung tâm 45,8 7,11% 0,305 7,98%
V1 (L) Thể tích phân bố ngăn trung tâm 6,73 61,1% 2,00 18,5%
Q (L/h) Độ thanh thải liên ngăn 21,2 4,82% 0,382 13,7%
V2 (L) Thể tích phân bố ngăn ngoại vi 1250 29,3% 0,397 38,1%
$\beta_{Cl_tFFM}$ Hệ số ảnh hưởng của FFM lên Cl 0,346 36,5% - -
$\beta_{Cl_Age}$ Hệ số ảnh hưởng của Tuổi lên Cl -0,00316 44,4% - -
b Hệ số sai số tỉ lệ 0,172 6,34% - -
corr_Cl_Q Tương quan hiệp phương sai Cl và Q 0,713 11,0% - -

Kết quả chính và đề xuất

  • Đã xác định được mô hình 2 ngăn kết hợp hấp thu tuần tự bậc 0 và bậc 1 có thời gian trễ là mô hình toán học mô tả tối ưu nhất động thái nồng độ ethambutol trong máu trên bệnh nhân lao phổi người Việt Nam.
  • Đã chứng minh FFM và Tuổi là hai biến số sinh lý có ý nghĩa thống kê quyết định sự biến thiên độ thanh thải ethambutol giữa các cá thể.
  • Đề xuất: Cần xem xét áp dụng chiến lược cá thể hóa chế độ liều ethambutol dựa trên cân nặng trừ mỡ (FFM) thay vì áp dụng thang chia liều theo các khoảng cân nặng cố định hiện hành của Chương trình Chống lao Quốc gia, đặc biệt trên những bệnh nhân nhiễm chủng vi khuẩn có mức MIC cao.

Đóng góp học thuật

  • Đây là nghiên cứu đầu tiên tại Việt Nam thiết lập thành công mô hình PopPK của ethambutol trên quần thể bệnh nhân lao phổi (bao gồm cả thể bệnh lao mới và lao tái trị).
  • Cung cấp bằng chứng định lượng về mối liên hệ giữa nồng độ thuốc và đáp ứng diệt khuẩn thông qua mô phỏng PK/PD Monte Carlo, tạo nền tảng lý thuyết và thực tiễn cho các nghiên cứu tối ưu hóa phác đồ điều trị lao tại Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận ghi nhận các hạn chế cụ thể được tác giả nêu trong văn bản:

  • Đặc điểm mẫu nghiên cứu: Toàn bộ 136 bệnh nhân được thu nhận tại 3 bệnh viện thuộc khu vực miền Bắc, chưa bao phủ được tính đại diện cho bệnh nhân lao trên toàn quốc.
  • Tỷ lệ đồng nhiễm HIV rất thấp: Chỉ ghi nhận 1/136 bệnh nhân nhiễm HIV (0,74%), do đó mô hình chưa thể đánh giá được tác động của tình trạng suy giảm miễn dịch và tương tác thuốc kháng retrovirus (ARV) lên dược động học của ethambutol.
  • Thiết kế lấy mẫu nồng độ: Áp dụng phương pháp lấy mẫu thưa (sparse sampling với 4 mẫu/bệnh nhân) khiến việc ước tính một số thông số phân bố ban đầu có sai số chuẩn tương đối còn cao (chẳng hạn thể tích ngăn trung tâm $V_1$ có RSE = 61,1%).
  • Dữ liệu vi sinh: Đánh giá PK/PD mới dừng lại ở bước mô phỏng toán học dựa trên các dải MIC giả định theo y văn quốc tế, chưa đo lường đồng thời giá trị MIC thực tế trên từng chủng M. tuberculosis phân lập trực tiếp từ các bệnh nhân tham gia nghiên cứu.

Hướng nghiên cứu tiếp theo: Mở rộng nghiên cứu đa trung tâm trên cỡ mẫu lớn hơn ở cả ba miền Bắc - Trung - Nam, bổ sung nhóm bệnh nhân đồng nhiễm HIV và tiến hành nuôi cấy, xác định MIC thực tế của từng chủng vi khuẩn nhằm xây dựng điểm gãy lâm sàng (clinical breakpoints) đặc thù cho bệnh nhân Việt Nam.

Giá trị tham khảo

Tài liệu khóa luận có giá trị ứng dụng và tham khảo khoa học rõ rệt:

  • Đối tượng thụ hưởng: Dược sĩ lâm sàng tại các bệnh viện chuyên khoa lao và bệnh phổi; bác sĩ điều trị bệnh truyền nhiễm; học viên sau đại học và sinh viên chuyên ngành Dược học, Dược lý - Dược lâm sàng hoặc Y đa khoa.
  • Nội dung giá trị nhất để tham khảo:
    • Quy trình và kỹ thuật xây dựng mô hình PopPK bằng phương pháp NLMEM trên phần mềm MONOLIX (từ khâu chuẩn hóa ma trận dữ liệu, thử nghiệm hàm số sai số đến thuật toán sàng lọc biến số ảnh hưởng).
    • Thuật toán và câu lệnh mô phỏng Monte Carlo trên gói phần mềm mlxR để đánh giá xác suất đạt đích PTA theo các dải MIC.
    • Bộ thông số dược động học tham chiếu của ethambutol trên người Việt Nam, đóng vai trò làm dữ liệu đầu vào cho các phần mềm hỗ trợ điều chỉnh liều lâm sàng (Bayesian dosing tools) hoặc theo dõi nồng độ thuốc điều trị (TDM).

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

1. Cấu trúc mô hình PopPK cuối cùng được lựa chọn cho Ethambutol có đặc điểm gì đặc biệt? Mô hình cấu trúc tối ưu là mô hình 2 ngăn có độ thanh thải tuyến tính bậc 1, đặc trưng bởi cơ chế hấp thu kép: hấp thu tuần tự bậc 0 và bậc 1, có kèm theo thời gian trễ ($T_{\text{lag}} = 0{,}051\text{ h}$) trước khi bắt đầu pha hấp thu bậc 0. Mô hình sai số thặng dư phù hợp nhất là mô hình sai số tỉ lệ ($b = 0{,}172$).

2. Những yếu tố lâm sàng nào ảnh hưởng có ý nghĩa thống kê đến độ thanh thải (Cl) của thuốc? Qua hai bước kiểm định Likelihood Ratio (thêm lần lượt $p < 0{,}05$ và loại trừ lần lượt $p < 0{,}01$), hai yếu tố dự đoán có ý nghĩa thống kê được giữ lại trong mô hình là Cân nặng trừ mỡ (FFM) với hệ số hồi quy $\beta = 0{,}346$ và Tuổi của bệnh nhân với hệ số hồi quy $\beta = -0{,}00316$.

3. Liều dùng Ethambutol thực tế trong mẫu nghiên cứu có đáp ứng đúng khuyến cáo của WHO và Bộ Y tế không? Trung vị liều tính theo cân nặng chung của mẫu là 18,18 mg/kg, nằm trong khoảng khuyến cáo (15–20 mg/kg). Tuy nhiên, do thực tế thuốc được chia theo các hàm lượng cố định (800 mg và 1200 mg), sự phân bố liều bị lệch: nhóm bệnh nhân 25–39 kg và $>70\text{ kg}$ có xu hướng dùng mức liều cao hơn khuyến cáo; một số cá thể nhận liều quá thấp (10,53 mg/kg) hoặc quá cao (25,00 mg/kg).

4. Kết quả mô phỏng PK/PD cho thấy nguy cơ thất bại điều trị xuất hiện ở ngưỡng MIC nào? Với chế độ liều hiện hành, khả năng đạt đích $f\text{AUC}/\text{MIC} > 11{,}9$ chỉ đạt tối ưu ($>98%$) khi vi khuẩn có $\text{MIC} \le 0{,}5\text{ mg/L}$. Khi $\text{MIC} = 1\text{ mg/L}$, PTA giảm xuống còn trên 50%; khi $\text{MIC} = 2\text{ mg/L}$, tỷ lệ đạt đích tụt dốc xuống dưới 10% và bằng 0% tại $\text{MIC} \ge 4\text{ mg/L}$.

5. Việc điều chỉnh liều theo cân nặng trừ mỡ (FFM) mang lại ưu thế gì qua mô phỏng? Mô phỏng 5 mức liều theo FFM (15, 20, 25, 30, 35 mg/kg FFM) cho thấy khi tăng liều theo FFM, tỷ lệ đạt đích PK/PD tăng rõ rệt trên tất cả các dải MIC. Tại mức $\text{MIC} = 4\text{ mg/L}$, PTA của đích $fC_{\max}/\text{MIC}$ tăng mạnh từ 55,3% ở mức liều 15 mg/kg FFM lên đến 97,9% ở mức liều 35 mg/kg FFM.

Kết luận

Khóa luận của tác giả Tạ Việt Hà đã xây dựng thành công mô hình dược động học quần thể 2 ngăn của ethambutol trên 136 bệnh nhân lao phổi tại Việt Nam bằng phương pháp NLMEM trên phần mềm MONOLIX. Nghiên cứu đã xác định định lượng hai yếu tố sinh lý có ảnh hưởng mang tính quyết định đến độ thanh thải thuốc là cân nặng trừ mỡ (FFM) và tuổi của bệnh nhân. Kết quả mô phỏng PK/PD Monte Carlo cảnh báo nguy cơ thất bại điều trị cao ở chế độ liều thông thường khi nồng độ ức chế tối thiểu MIC của vi khuẩn từ 1 mg/L trở lên. Công trình cung cấp luận cứ khoa học quan trọng cho việc chuyển đổi phương pháp phân liều từ bậc cân nặng truyền thống sang cá thể hóa liều theo cân nặng trừ mỡ (FFM) nhằm nâng cao hiệu quả điều trị và hạn chế độc tính trên thị giác cho bệnh nhân lao.