Bối cảnh và vấn đề nghiên cứu
Ghép thận là phương pháp điều trị tối ưu cho người bệnh mắc bệnh thận mạn tính giai đoạn cuối. Thành công của quá trình ghép thận phụ thuộc vào liệu pháp ức chế miễn dịch duy trì nhằm ngăn ngừa phản ứng thải ghép. Phác đồ chuẩn thường được áp dụng là sự phối hợp của bộ ba thuốc: chất ức chế calcineurin (CNIs), mycophenolate mofetil (MMF) và corticosteroid, trong đó CNIs đóng vai trò trọng tâm. Hiện nay, tacrolimus đã thay thế phần lớn cyclosporin do hiệu lực ức chế miễn dịch cao và hiệu quả lâm sàng đã được khẳng định.
Tuy nhiên, tacrolimus là thuốc có khoảng điều trị hẹp và ghi nhận sự biến thiên dược động học rất lớn giữa các cá thể cũng như trong cùng một cá thể theo thời gian. Nồng độ thuốc trong máu thấp làm tăng nguy cơ thải ghép cấp tính, trong khi nồng độ cao dẫn đến các độc tính nghiêm trọng trên thận, thần kinh, đái tháo đường sau ghép, rối loạn tiêu hóa và nhiễm trùng. Trong thực hành lâm sàng, việc giám sát nồng độ trị liệu (TDM) tacrolimus chủ yếu dựa vào nồng độ đáy trong máu toàn phần ($C_0$). Cách tiếp cận này bộc lộ nhược điểm là làm chậm trễ việc điều chỉnh liều, đặc biệt ở giai đoạn sớm sau ghép, đồng thời nồng độ đáy không phản ánh hoàn toàn chính xác mức độ phơi nhiễm toàn thân (diện tích dưới đường cong - AUC).
Để giải quyết vấn đề cá thể hóa liều dùng, các mô hình dược động học quần thể (Pop-PK) kết hợp thuật toán ước tính Bayes (Bayesian estimation) đã được nghiên cứu và ứng dụng trên thế giới nhằm hỗ trợ xác định liều khởi đầu tối ưu (tiếp cận a priori) và hiệu chỉnh liều duy trì dựa trên nồng độ đo được (tiếp cận a posteriori). Tại Việt Nam, Bệnh viện Quân y 103 là một trong những trung tâm lớn về ghép thận và đang triển khai TDM tacrolimus thường quy. Dù vậy, chưa có mô hình Pop-PK nào được xây dựng riêng cho người Việt, và khả năng dự đoán của các mô hình công bố quốc tế trên quần thể bệnh nhân trong nước vẫn chưa được kiểm chứng qua phương pháp thẩm định ngoại (external validation).
Đề tài đặt ra hai mục tiêu nghiên cứu cụ thể:
- Xác định mô hình Pop-PK đã được công bố trong y văn phù hợp với quần thể bệnh nhân ghép thận đang được theo dõi và điều trị tại Bệnh viện Quân y 103.
- Đánh giá hiệu suất dự đoán của các mô hình Pop-PK tacrolimus phù hợp với quần thể bệnh nhân ghép thận đang được theo dõi và điều trị tại Bệnh viện Quân y 103.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng lâm sàng: 45 bệnh nhân ghép thận trưởng thành đang được điều trị và theo dõi ngoại trú ổn định bằng phác đồ tacrolimus dạng giải phóng ngay (Prograf), dùng liều mỗi 12 giờ kết hợp MMF và corticosteroid tại Bệnh viện Quân y 103.
- Đối tượng dữ liệu y văn: Các bài báo công bố mô hình Pop-PK tacrolimus trên cơ sở dữ liệu PubMed tính đến tháng 5 năm 2022.
- Phạm vi thời gian theo dõi dữ liệu lâm sàng: Từ ngày 01/12/2021 đến ngày 31/05/2022.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và dược động học quần thể
Tacrolimus là một lacton macrolid phân lập từ Streptomyces tsukubaensis. Thuốc gắn với protein nội bào immunophilin FKBP-12 tạo thành phức hợp ức chế enzym calcineurin phosphatase, ngăn chặn quá trình dephosphoryl hóa yếu tố nhân của tế bào T hoạt hóa (NF-AT), từ đó ức chế phiên mã các gen sinh tổng hợp cytokin (đặc biệt là Interleukin-2), ức chế sự hoạt hóa và tăng sinh tế bào lympho T.
Về mặt dược động học, tacrolimus hấp thu nhanh qua đường tiêu hóa ($T_{max}$ từ 0,5 đến 1 giờ), sinh khả dụng đường uống thấp và dao động lớn (trung bình 25%, khoảng 4–93%) do khả năng thấm kém, bị chuyển hóa bởi enzym CYP3A4 ở ruột và bị đẩy ngược vào lòng ruột bởi bơm tống thuốc P-glycoprotein. Trong máu toàn phần, trên 95% tacrolimus liên kết với hồng cầu, thể tích phân bố trung bình khoảng 1 L/kg; trong huyết tương, 72–77% thuốc gắn với protein (chủ yếu là $\alpha_1$-acid glycoprotein và albumin). Thuốc chuyển hóa chủ yếu qua hệ enzym cytochrom P450 (CYP3A4 và CYP3A5) tại gan và thành ruột thành hơn 15 chất chuyển hóa không hoặc có rất ít hoạt tính, thải trừ trên 95% qua đường mật và dưới 1% qua nước tiểu dưới dạng nguyên vẹn, thời gian bán thải ($t_{1/2}$) trung bình là 11,3 giờ.
Một mô hình Pop-PK hiệu ứng hỗn hợp phi tuyến tính (non-linear mixed effects model) bao gồm 3 mô hình con (sub-models):
- Mô hình cấu trúc (Structural sub-model): Biểu diễn số ngăn (1 ngăn, 2 ngăn), các thông số độ thanh thải ($Cl/F$), thể tích phân bố ($V/F$) và tốc độ hấp thu ($K_a$).
- Mô hình thống kê (Statistical sub-model): Định lượng dao động thông số giữa các cá thể (BSV) theo mô hình cộng, mô hình tỉ lệ hoặc mô hình mũ ($P_i = P_{pop} \times \exp(\eta_i)$); dao động trong cùng một cá thể theo thời gian (BOV); và sai số ngẫu nhiên của phép đo ($RUV$) theo mô hình cộng, mô hình tỉ lệ hoặc mô hình kết hợp ($C_{obs} = C_{pred} + C_{pred} \cdot \varepsilon_1 + \varepsilon_2$).
- Mô hình hiệp biến (Covariate sub-model): Đưa các đặc tính sinh lý, bệnh lý và di truyền vào phương trình để giải thích nguồn gốc biến thiên. Các hiệp biến quan trọng nhất ảnh hưởng đến thanh thải tacrolimus gồm: kiểu gen CYP3A5 (người mang alen $*1$ có biểu hiện enzym nên thanh thải cao hơn người mang kiểu gen đồng hợp tử đột biến $33$), chỉ số hematocrit (HCT - HCT giảm làm tăng tỷ lệ thuốc tự do trong huyết tương, dẫn đến tăng thanh thải tại gan), số ngày hậu phẫu (POD) và cân nặng (WT).
Phương pháp nghiên cứu và thẩm định ngoại
Nghiên cứu sử dụng phương pháp thẩm định ngoại (external validation), đánh giá khả năng dự đoán nồng độ của các mô hình Pop-PK từ y văn trên một bộ dữ liệu lâm sàng độc lập hoàn toàn.
- Thu thập dữ liệu lâm sàng: Nghiên cứu mô tả tiến cứu trên bệnh nhân ngoại trú. Thời gian uống thuốc và thời điểm lấy mẫu máu vào buổi sáng được kiểm soát chính xác. Nồng độ đáy tacrolimus trong máu toàn phần được định lượng bằng phương pháp xét nghiệm miễn dịch vi tiểu phân hóa phát quang (CMIA) trên hệ thống Architect (Abbott Diagnostics, giới hạn phát hiện 0,8 ng/ml).
- Chiến lược sàng lọc y văn: Tìm kiếm trên cơ sở dữ liệu PubMed với chuỗi từ khóa:
("kidney" OR "renal") AND "tacrolimus" AND "Bayesian" AND "population" AND "transplant". Các tiêu chí lựa chọn: nghiên cứu Pop-PK bằng tiếng Anh, trên người lớn ghép thận, dùng phác đồ tacrolimus uống 2 lần/ngày, có kết hợp yếu tố di truyền CYP3A5. Loại trừ các nghiên cứu trên trẻ em, nghiên cứu thẩm định lại, mô hình chứa gen khác hoặc chứa các biến không có trong bộ dữ liệu của bệnh viện, hoặc mô hình thiếu thông số tính toán. - Công cụ và chỉ số đánh giá:
- Phần mềm MONOLIX 2020R1 được sử dụng để chạy mô phỏng dự đoán nồng độ theo hai cách tiếp cận: a priori (dựa trên thông số quần thể và hiệp biến cá thể) và tiếp cận Bayes (MAP-Bayesian, kết hợp thông số mô hình với các điểm nồng độ đo được).
- Phần mềm SPSS 20 dùng để phân tích thống kê mô tả đặc điểm bệnh nhân.
- Phần mềm Microsoft Excel tính toán các chỉ số:
- Sai số dự đoán trung bình (Mean Prediction Error - MPE), phản ánh độ đúng (bias): $$MPE = \frac{\sum_{j=1}^{n} (pred_j - obs_j)}{n}$$
- Trung bình giá trị tuyệt đối sai số dự đoán (Mean Absolute Prediction Error - MAPE), phản ánh độ chính xác (precision): $$MAPE = \frac{\sum_{j=1}^{n} \left| \frac{pred_j - obs_j}{obs_j} \right|}{n} \times 100%$$
- Phân loại mức độ dự đoán: Rất tốt ($MAPE < 10%$), Tốt ($10% \le MAPE < 20%$), Chấp nhận được ($20% \le MAPE < 50%$), Kém ($MAPE \ge 50%$).
- Đánh giá đồ thị khớp mô hình (Goodness of Fit - GOF) giữa nồng độ quan sát và nồng độ dự đoán so với đường chéo chuẩn $y = x$.
Nội dung chính theo từng chương
Tổng quan về tacrolimus trên bệnh nhân ghép thận
Nội dung chương tổng quan hệ thống hóa lịch sử phát triển của tacrolimus từ khi được FDA phê duyệt năm 1997 cho chỉ định ghép thận, khẳng định vai trò thay thế cyclosporin trong phác đồ ức chế miễn dịch phối hợp kinh điển. Chương này làm rõ cơ chế ức chế calcineurin, giải thích chuỗi tương tác phân tử dẫn đến phong tỏa sao chép IL-2 và chọn lọc tế bào lympho T. Mối tương quan giữa phơi nhiễm toàn thân với hiệu quả chống thải ghép và các độc tính được luận giải sâu sắc: nồng độ đáy giảm mỗi 1 ng/ml làm tăng 7,2% nguy cơ thải ghép cấp trong 6 tháng đầu; ngược lại, nồng độ duy trì trên 20 ng/ml làm gia tăng mạnh tỷ lệ độc thận, run tay và đái tháo đường sau ghép.
Chương tổng quan cũng phân tích nguyên lý xây dựng mô hình Pop-PK ảnh hưởng hỗn hợp phi tuyến tính, chỉ rõ vai trò của 4 hiệp biến lâm sàng cốt lõi (POD, CYP3A5, HCT, WT). Tổng kết y văn từ 56 nghiên cứu Pop-PK quốc tế chỉ ra rằng kiểu gen CYP3A5 và HCT ảnh hưởng nhiều nhất đến độ thanh thải $Cl/F$. Đồng thời, nguyên lý hiệu chỉnh liều chính xác dựa trên mô hình (MIPD) và thuật toán ước lượng hậu nghiệm cực đại Bayes (MAPB) được trình bày như một giải pháp khắc phục triệt để hạn chế của phương pháp TDM truyền thống.
Đối tượng và phương pháp nghiên cứu
Chương phương pháp trình bày chi tiết quy trình tiến cứu thu thập dữ liệu lâm sàng trên 45 bệnh nhân ghép thận ngoại trú tại Bệnh viện Quân y 103 và quy trình sàng lọc y văn trên cơ sở dữ liệu PubMed.
Tiêu chuẩn lựa chọn bệnh nhân được kiểm soát nghiêm ngặt: bệnh nhân trưởng thành dùng phác đồ bộ ba với Prograf ổn định liều mỗi 12 giờ trong ít nhất 1 tuần trước khi lấy máu, ghi nhận chính xác giờ uống liều kế trước và giờ trích máu định lượng nồng độ đáy $C_0$. Phương pháp định lượng hóa phát quang CMIA và các thông số cận lâm sàng (creatinin, albumin, men gan, HCT, định nhóm gen CYP3A5) được chuẩn hóa.
Quy trình sàng lọc mô hình Pop-PK y văn được triển khai qua các bước lọc chặt chẽ, loại bỏ các mô hình thiếu thông số hoặc chứa hiệp biến ngoài khả năng thu thập của bệnh viện. Phương pháp phân tích số liệu trên phần mềm MONOLIX 2020R1 và các công thức toán học xác định độ đúng (MPE), độ chính xác (MAPE) cùng biểu đồ phân tán GOF được định nghĩa rõ ràng.
Kết quả nghiên cứu
Kết quả nghiên cứu được trình bày mạch lạc theo hai phần chính: phân tích đặc điểm bộ dữ liệu lâm sàng tại Bệnh viện Quân y 103 và đánh giá hiệu suất dự đoán của các mô hình Pop-PK được chọn lọc từ y văn.
1. Đặc điểm quần thể bệnh nhân tại Bệnh viện Quân y 103
Nghiên cứu thu thập dữ liệu trên 45 bệnh nhân với tổng số 112 điểm đo nồng độ đáy $C_0$ (trung bình 2,5 điểm/bệnh nhân). Quần thể mang đặc trưng rõ nét của nhóm bệnh nhân ghép thận giai đoạn muộn.
| Đặc điểm lâm sàng / cận lâm sàng | Giá trị ghi nhận ($N = 45$) |
|---|---|
| Giới tính: Nam, n (%) | 33 (73,7%) |
| Giới tính: Nữ, n (%) | 12 (26,3%) |
| Tuổi (năm), Trung bình $\pm$ SD | 39,8 $\pm$ 8,3 (khoảng 23 – 60) |
| Cân nặng (kg), Trung vị (Khoảng tứ phân vị) | 58,0 (53,0 – 63,5) |
| Hematocrit (%), Trung bình $\pm$ SD | 42,8 $\pm$ 3,8 |
| Kiểu gen CYP3A5 $11$ hoặc $13$, n (%) | 21 (46,5%) |
| Kiểu gen CYP3A5 $33$, n (%) | 24 (53,5%) |
| Thời gian hậu phẫu (ngày), Trung bình $\pm$ SD | 1429,6 $\pm$ 208,6 |
| Tổng liều tacrolimus/ngày (mg), Trung vị (IQR) | 6,0 (4,0 – 6,9) |
| Nồng độ đáy $C_0$ (ng/ml), Trung bình $\pm$ SD | 5,5 $\pm$ 1,4 (khoảng 2,6 – 12,3) |
2. Kết quả sàng lọc các mô hình Pop-PK từ y văn
Từ 256 bài báo tìm thấy trên PubMed, sau khi loại bỏ bài trùng lặp (còn 160 bài), rà soát tiêu đề/tóm tắt và áp dụng các tiêu chí loại trừ (không phải Pop-PK, nghiên cứu nhi khoa, thiếu yếu tố CYP3A5, chứa hiệp biến không tương thích), tác giả chọn được 6 mô hình Pop-PK phù hợp nhất để đưa vào đánh giá ngoại kiểm.
| Thứ tự | Tác giả và năm công bố | Cấu trúc mô hình | Hiệp biến trong mô hình cuối cùng | Mô hình sai số ($RUV$) | Phương pháp thẩm định gốc |
|---|---|---|---|---|---|
| 1 | Han và cs. (2013) | 1 ngăn, hấp thu và thải trừ bậc 1 | POD, CYP3A5, HCT, WT | Sai số tỉ lệ (40,0%) | Bootstrap, RMSE |
| 2 | Zhang và cs. (2017) | 1 ngăn, hấp thu và thải trừ bậc 1 | POD, CYP3A5, HCT | Sai số cộng (2,33 ng/ml) | GOF, VPC, Bootstrap |
| 3 | Han và cs. (2014) | 1 ngăn, có thời gian trễ $T_{lag} = 0,25$ h | POD, CYP3A5 | Sai số tỉ lệ (5,4%) | VPC, Bootstrap |
| 4 | Woillard và cs. (2011) | 2 ngăn, thời gian trễ, hấp thu Erlang | CYP3A5, HCT | Sai số kết hợp: cộng (0,71 ng/ml) + tỉ lệ (11,3%) | VPC |
| 5 | Zhu và cs. (2018) | 1 ngăn, hấp thu và thải trừ bậc 1 | POD, CYP3A5, HCT, WT | Sai số kết hợp: cộng (2,5 ng/ml) + tỉ lệ (18,8%) | Bootstrap |
| 6 | Ling và cs. (2020) | 1 ngăn, hấp thu và thải trừ bậc 1 | POD, CYP3A5, HCT | Sai số kết hợp: cộng (1,40 ng/ml) + tỉ lệ (24,4%) | GOF, Bootstrap |
3. Đánh giá hiệu suất dự đoán của các mô hình
Hiệu suất dự đoán nồng độ tacrolimus trên bộ dữ liệu 112 điểm nồng độ được kiểm tra qua hai phương thức: tiếp cận ban đầu a priori và tiếp cận hiệu chỉnh Bayes.
| Mô hình nghiên cứu | Tiếp cận a priori: MPE (ng/ml) | Tiếp cận a priori: MAPE (%) | Tiếp cận Bayesian: MPE (ng/ml) | Tiếp cận Bayesian: MAPE (%) | Đánh giá khả năng dự đoán (Bayesian) |
|---|---|---|---|---|---|
| Han và cs. (2013) | 0,97 | 41,14 | -0,48 | 14,46 | Tốt ($10% < MAPE < 20%$) |
| Zhang và cs. (2017) | -3,54 | 67,30 | -3,16 | 61,00 | Kém ($MAPE > 50%$) |
| Han và cs. (2014) | -5,45 | 100,00 | -0,69 | 21,52 | Dự đoán ra nồng độ âm |
| Woillard và cs. (2011) | 1,60 | 52,43 | 0,00 | 13,29 | Tốt nhất ($MAPE$ thấp nhất) |
| Zhu và cs. (2018) | 1,39 | 44,44 | 0,07 | 19,06 | Tốt ($10% < MAPE < 20%$) |
| Ling và cs. (2020) | 4,73 | 96,69 | 1,38 | 32,95 | Chấp nhận được |
Khi sử dụng tiếp cận a priori, tất cả các mô hình đều cho kết quả dự đoán kém hoặc chỉ ở mức chấp nhận được, sai số tuyệt đối rất cao ($MAPE$ từ 41,14% đến 100,00%), đường xu hướng trên biểu đồ GOF lệch xa đường thẳng $y = x$.
Ngược lại, khi áp dụng tiếp cận Bayes, khả năng dự đoán của các mô hình được cải thiện vượt bậc. Ba mô hình đạt mức dự đoán tốt ($10% < MAPE < 20%$) gồm: Woillard và cs. (2011) có $MAPE = 13,29%$ và $MPE = 0,00$; Han và cs. (2013) có $MAPE = 14,46%$ và $MPE = -0,48$; Zhu và cs. (2018) có $MAPE = 19,06%$ và $MPE = 0,07$. Trên đồ thị GOF của 3 mô hình này, các điểm phân bố đồng nhất và đường xu hướng bám sát đường lý tưởng $y = x$. Riêng mô hình của Han và cs. (2014) cho ra một số giá trị nồng độ âm nên bị loại bỏ khỏi ứng dụng; mô hình của Zhang và cs. (2017) tiếp tục dự đoán thấp hơn nhiều so với thực tế ($MPE = -3,16$; $MAPE = 61,00%$).
Bàn luận
Chương bàn luận phân tích chuyên sâu các yếu tố giải thích cho sự sai biệt và tương thích giữa mô hình y văn với dữ liệu lâm sàng thực tế:
- Tác động của yếu tố di truyền CYP3A5 và HCT: Cả 6 mô hình đều chỉ ra vai trò quyết định của enzym CYP3A5 lên độ thanh thải. Tỷ lệ bệnh nhân mang kiểu gen đồng hợp tử đột biến $33$ tại Bệnh viện Quân y 103 chiếm 53,5%, tương đồng với đặc điểm dịch tễ của các nghiên cứu tại châu Á. Những cá thể này có độ thanh thải tacrolimus thấp hơn rõ rệt so với người mang alen $*1$ ($11$ hoặc $13$), do đó cần chế độ liều thấp hơn để đạt cùng một nồng độ đích. Về chỉ số HCT, mức trung bình của quần thể nghiên cứu đạt 42,8%, cao hơn so với các nghiên cứu xây dựng mô hình gốc (30,9% – 39,1%). Vì tacrolimus gắn chủ yếu vào hồng cầu, HCT tăng làm giảm tỷ lệ phân bố thuốc tự do trong huyết tương, dẫn đến giảm thanh thải qua gan.
- Ảnh hưởng của thời gian hậu phẫu (POD): Bộ dữ liệu Bệnh viện Quân y 103 có thời gian sau ghép rất dài (trung bình 1429,6 ngày, tương đương gần 4 năm). Đây là nguyên nhân trực tiếp khiến các mô hình của Zhang và cs. (xây dựng trên bệnh nhân dưới 400 ngày), Han và cs. 2014 (dưới 20 ngày) và Ling và cs. (dưới 31 ngày) thất bại khi ngoại kiểm, dẫn đến hiện tượng nồng độ âm hoặc sai lệch nghiêm trọng. Ngược lại, mô hình của Woillard và cs. (2011) không đưa POD vào mô hình cuối cùng mà chỉ tập trung vào HCT và CYP3A5, giúp mô hình giữ được độ ổn định rất cao khi áp dụng trên bệnh nhân ghép thận lâu năm.
- Giá trị của tiếp cận Bayes và cấu trúc mô hình: Kết quả thẩm định ngoại chứng minh mô hình 2 ngăn kết hợp hấp thu Erlang của Woillard và cs. (2011) mô tả sinh lý phân bố thuốc tốt hơn hẳn các mô hình 1 ngăn giản đơn khi kết hợp công cụ Bayes. Sự kết hợp giữa các điểm nồng độ quan sát thực tế với thông số quần thể giúp thuật toán MAPB hiệu chỉnh chính xác các tham số cá thể, khắc phục triệt để nhược điểm thiếu chính xác của phương pháp dự đoán thuần túy a priori.
- Ưu điểm và hạn chế: Ưu điểm nổi bật của nghiên cứu là thiết kế tiến cứu với việc giám sát chính xác giờ uống thuốc và giờ lấy máu, kiểm soát được sự tuân thủ tốt hơn các nghiên cứu hồi cứu. Hạn chế bao gồm việc cỡ mẫu còn khiêm tốn (45 bệnh nhân), bệnh nhân ngoại trú chỉ đo nồng độ đáy định kỳ 1 tháng/lần, thiếu thông tin về diện tích bề mặt cơ thể (BSA), nồng độ corticosteroid dùng kèm và kiểu gen CYP3A4, cũng như sự khác biệt nhất định giữa phương pháp định lượng CMIA tại viện với LC-MS/MS hoặc MEIA trong các bài báo gốc.
Kết quả và đóng góp
Kết quả định lượng cụ thể
- Đã sàng lọc từ 256 bài báo quốc tế trên PubMed và xác định được 6 mô hình Pop-PK tacrolimus thỏa mãn đầy đủ các tiêu chuẩn hiệp biến trên bệnh nhân ghép thận trưởng thành.
- Thẩm định ngoại trên bộ dữ liệu 45 bệnh nhân (112 điểm nồng độ) tại Bệnh viện Quân y 103, xác định rõ:
- Tiếp cận a priori không đủ độ chính xác để áp dụng lâm sàng ($MAPE$ từ 41,14% đến 100,00%).
- Tiếp cận Bayes cải thiện vượt trội hiệu suất dự đoán, trong đó xác định được 3 mô hình đạt độ chính xác tốt ($MAPE < 20%$): Han và cs. 2013 ($MAPE = 14,46%$), Zhu và cs. 2018 ($MAPE = 19,06%$) và Woillard và cs. 2011 ($MAPE = 13,29%$).
- Mô hình 2 ngăn của Woillard và cs. (2011) là mô hình tối ưu nhất cho quần thể bệnh nhân ghép thận tại Bệnh viện Quân y 103 với độ đúng tuyệt đối ($MPE = 0,00$) và sai số dự đoán thấp nhất ($MAPE = 13,29%$). Mô hình có các thông số cụ thể:
- Hằng số tốc độ chuyển tiếp hấp thu $K_{tr} = 5,11\text{ h}^{-1}$ (IIV 24%), thời gian vận chuyển trung bình $M_{tt} = 0,78\text{ h}$, hằng số hấp thu $K_a = 7,31\text{ h}^{-1}$ (IIV 32,6%).
- Phương trình thanh thải: $$Cl/F\text{ (L/h)} = 21,2 \times \left[\frac{HCT}{35}\right]^{-1,14} \times (2^{CYP3A5})$$ *(trong đó $CYP3A5 = 0$ nếu mang kiểu gen $33$, $CYP3A5 = 1$ nếu mang alen $1$), với biến thiên giữa các cá thể IIV là 28%.
- Thể tích phân bố ngăn trung tâm $V_1/F = 140,94\text{ L}$ (IIV 31%), thanh thải liên ngăn $Q/F = 79\text{ L/h}$ (IIV 54%), thể tích phân bố ngăn ngoại vi $V_2/F = 271\text{ L}$ (IIV 60%).
- Mô hình sai số kết hợp: sai số cộng $0,71\text{ ng/ml}$ và sai số tỉ lệ $11,3%$.
Đóng góp khoa học và thực tiễn
- Đây là công trình đầu tiên tại Việt Nam thực hiện thẩm định ngoại các mô hình Pop-PK tacrolimus có tích hợp yếu tố đa hình di truyền CYP3A5 trên quần thể bệnh nhân ghép thận người Việt.
- Khẳng định tính khả thi của việc ứng dụng thuật toán Bayes kết hợp các mô hình dược động học công bố quốc tế để cá thể hóa chế độ liều tacrolimus mà không nhất thiết phải phát triển lại một mô hình hoàn toàn mới từ đầu.
- Cung cấp bằng chứng thực nghiệm rõ ràng về mô hình tối ưu (Woillard và cs., 2011), tạo tiền đề khoa học để xây dựng công cụ phần mềm hỗ trợ ra quyết định lâm sàng (MIPD) nhằm nâng cao độ an toàn và hiệu quả điều trị ức chế miễn dịch tại các trung tâm ghép tạng.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của nghiên cứu
- Đặc điểm quần thể nghiên cứu thuần nhất ở nhóm bệnh nhân ghép thận lâu năm (thời gian hậu phẫu trung bình gần 4 năm, $POD = 1429,6\text{ ngày}$), chưa bao phủ được nhóm bệnh nhân trong giai đoạn sớm ngay sau ghép (dưới 1 tháng đến 6 tháng đầu) – vốn là giai đoạn có biến động dược động học mạnh nhất và nguy cơ thải ghép cấp cao nhất.
- Cỡ mẫu còn hạn chế (45 bệnh nhân với 112 điểm nồng độ đáy), số điểm đo nồng độ trên từng bệnh nhân còn ít (trung bình 2,5 điểm/người) do thời gian theo dõi tiến cứu chỉ kéo dài 6 tháng.
- Chưa thu thập và đánh giá được một số hiệp biến tiềm năng đã được ghi nhận trong y văn như diện tích bề mặt cơ thể (BSA), nồng độ và liều dùng của corticosteroid, cũng như tình trạng đa hình di truyền của enzym CYP3A4.
- Sự sai biệt về kỹ thuật xét nghiệm định lượng thuốc giữa phương pháp CMIA sử dụng tại Bệnh viện Quân y 103 và phương pháp sắc ký lỏng khối phổ (LC-MS/MS) hoặc xét nghiệm miễn dịch enzym (EMIT, MEIA) ở các nghiên cứu gốc có thể là một nguồn gây sai số ngoại định.
Hướng nghiên cứu tiếp theo
- Mở rộng quy mô nghiên cứu trên tập dữ liệu đa trung tâm với cỡ mẫu lớn hơn, đặc biệt thu thập nhóm bệnh nhân ghép thận mới trong giai đoạn chu phẫu và hậu phẫu sớm.
- Nâng cấp hệ thống số hóa bệnh án điện tử, tích hợp thu thập đầy đủ các dữ liệu nhân trắc học, phác đồ ức chế miễn dịch chi tiết và xét nghiệm mở rộng các kiểu gen chuyển hóa thuốc.
- Phát triển và thử nghiệm lâm sàng công cụ phần mềm tích hợp thuật toán Bayes dựa trên mô hình Woillard 2011 nhằm hướng dẫn bác sĩ và dược sĩ lâm sàng hiệu chỉnh liều chính xác tại giường bệnh.
Giá trị tham khảo
Khóa luận là tài liệu chuyên khảo có giá trị ứng dụng cao cho nhiều nhóm đối tượng:
- Dược sĩ lâm sàng tại các bệnh viện ghép tạng: Cung cấp cơ sở khoa học để chuyển đổi từ việc theo dõi nồng độ đáy $C_0$ thủ công sang ứng dụng công cụ Bayes dự đoán liều và phơi nhiễm toàn thân (AUC), hỗ trợ cá thể hóa liều dùng tacrolimus an toàn.
- Bác sĩ chuyên khoa ghép thận và thận học: Hiểu rõ tác động định lượng của kiểu gen CYP3A5 và chỉ số HCT lên độ thanh thải thuốc, hỗ trợ giải thích các trường hợp nồng độ thuốc dao động bất thường hoặc khó đạt nồng độ mục tiêu.
- Học viên, sinh viên ngành Dược và Y: Là tài liệu tham khảo mẫu mực về phương pháp luận nghiên cứu Dược động học quần thể (Pop-PK), kỹ thuật thẩm định ngoại mô hình (external validation), việc sử dụng phần mềm MONOLIX và các phép kiểm định thống kê y sinh học.
- Nhóm nghiên cứu phát triển công nghệ y tế: Cung cấp bộ thông số toán học chuẩn xác của mô hình Woillard và cs. (2011) để lập trình các phân hệ hiệu chỉnh liều chính xác (MIPD) tích hợp vào phần mềm quản lý bệnh viện (HIS).
Câu hỏi thường gặp
1. Tại sao kiểu gen CYP3A5 lại được chọn làm tiêu chí sàng lọc bắt buộc cho các mô hình Pop-PK trong nghiên cứu này?
Enzym CYP3A5 đóng vai trò chủ đạo trong chuyển hóa tacrolimus in vivo. Người mang ít nhất một alen $*1$ ($11$ hoặc $13$) có biểu hiện enzym hoạt tính đầy đủ, khiến tốc độ thanh thải tacrolimus tăng cao rõ rệt và đòi hỏi mức liều thuốc lớn hơn để đạt cùng nồng độ đáy mục tiêu. Ngược lại, người mang kiểu gen đồng hợp tử biến thể $33$ bị mất chức năng enzym, dẫn đến giảm thanh thải thuốc. Tại quần thể nghiên cứu Bệnh viện Quân y 103, tỷ lệ bệnh nhân mang gen $33$ chiếm tới 53,5%, tương đương với đặc điểm chung của người châu Á, do đó bắt buộc phải chọn mô hình chứa yếu tố CYP3A5 để đảm bảo tính chuẩn xác.
2. Sự khác biệt cốt lõi về hiệu suất dự đoán giữa tiếp cận a priori và tiếp cận Bayesian là gì?
Tiếp cận a priori chỉ dựa vào các thông số đặc trưng cố định của mô hình quần thể kết hợp với các hiệp biến của người bệnh (tuổi, cân nặng, HCT, gen) để suy đoán nồng độ, hoàn toàn chưa có dữ liệu đo thực tế của cá thể đó. Do độ dao động cá thể của tacrolimus rất lớn, tiếp cận này cho kết quả kém với sai số rất cao ($MAPE$ từ 41,14% đến 100,00%). Trong khi đó, tiếp cận Bayes kết hợp thông số mô hình với chính các điểm nồng độ đã đo được trên người bệnh ($C_0$), giúp thuật toán hiệu chỉnh lại tham số riêng của cá thể đó, từ đó giảm sai số ngoạn mục đưa $MAPE$ của 3 mô hình xuống dưới 20%.
3. Vì sao mô hình của Woillard và cs. (2011) lại đạt hiệu suất dự đoán tốt nhất trên bệnh nhân Bệnh viện Quân y 103?
Mô hình Woillard và cs. (2011) đạt kết quả tối ưu nhất với $MPE = 0,00$ và $MAPE = 13,29%$ nhờ hai yếu tố:
- Thứ nhất, cấu trúc 2 ngăn kết hợp hấp thu Erlang và thời gian trễ phản ánh sát thực tế sinh lý phân bố và chuyển tiếp của tacrolimus hơn các mô hình 1 ngăn.
- Thứ hai, mô hình chỉ tích hợp 2 hiệp biến sinh học mang tính ổn định cao là hematocrit và kiểu gen CYP3A5, không đưa số ngày hậu phẫu (POD) vào phương trình. Điều này giúp mô hình không bị ảnh hưởng bởi sự sai lệch thời gian hậu phẫu dài (trung bình gần 4 năm) của bệnh nhân tại Bệnh viện Quân y 103.
4. Tại sao các mô hình Pop-PK của Zhang và cs. (2017) và Han và cs. (2014) lại dự đoán kém trên bộ dữ liệu nghiên cứu?
Cả hai mô hình này đều đưa số ngày hậu phẫu (POD) vào làm hiệp biến tác động đến độ thanh thải và thể tích phân bố. Tuy nhiên, mô hình của Han và cs. (2014) được xây dựng trên bệnh nhân mới ghép dưới 20 ngày, còn mô hình của Zhang và cs. (2017) xây dựng trên bệnh nhân dưới 400 ngày. Khi áp dụng các hàm số mũ hoặc hàm bậc nhất của POD này vào quần thể Bệnh viện Quân y 103 có thời gian hậu phẫu trung bình lên tới 1429,6 ngày, phương trình toán học bị ngoại suy quá mức khỏi khoảng giá trị ban đầu, dẫn đến việc mô hình Han 2014 dự đoán ra nồng độ âm và mô hình Zhang 2017 có sai số rất lớn ($MAPE = 61,00%$).
Kết luận
Khóa luận tốt nghiệp của dược sĩ Nguyễn Minh Thành đã thực hiện thành công việc thẩm định ngoại 6 mô hình Pop-PK tacrolimus từ y văn quốc tế trên bộ dữ liệu lâm sàng tiến cứu gồm 45 bệnh nhân ghép thận tại Bệnh viện Quân y 103. Nghiên cứu chứng minh rằng tiếp cận Bayes vượt trội hoàn toàn so với tiếp cận a priori, giúp cải thiện đáng kể độ chính xác trong dự đoán nồng độ thuốc. Trong các mô hình được đánh giá, mô hình 2 ngăn của Woillard và cs. (2011) kết hợp hiệp biến CYP3A5 và hematocrit cho hiệu suất dự đoán ưu việt nhất ($MAPE = 13,29%$, $MPE = 0,00$), khẳng định khả năng ứng dụng thực tế vào phần mềm hỗ trợ hiệu chỉnh liều cá thể hóa tacrolimus trong lâm sàng ghép tạng tại Việt Nam.