Journal List > Transl Clin Pharmacol > v.27(2) > 1138130

Cho and Yoon: Understanding the pharmacokinetics of reversible metabolism

Abstract

This tutorial introduces the mathematical skills required to obtain exact and approximate solutions for reversible reactions and provides graphical insights to help understand the pharmacokinetics of reversible metabolism. The matrix method provides an easy way to derive the exact solution for the amount of each species as a function of time. The plots of the exact solutions reveal some characteristic features of the pharmacokinetic profiles of the reversible metabolism. We also describe two approximation approaches, steady-state approximation, and equilibrium approximation, to simplify the solutions. The skills and knowledge acquired through this tutorial will provide a basis for understanding more complex reversible reaction systems.

Introduction

Reversible reactions are commonly encountered in pharmacology, chemistry, biology, chemical engineering, etc.[123456] The influence of reversible metabolism on the pharmacokinetics of a drug and its metabolites has been the subject of considerable recent research.[567891011] An interesting example is the administration of chiral drugs. Two enantiomers of a chiral drug, such as ibuprofen[10] and propranolol,[11] often have different potencies and pharmacokinetic profiles.[891011] However, one enantiomer is usually interconverted into the other within the human body.[7891011] When a chiral drug is administered as a racemic mixture, its bioavailability and efficacy can be incorrectly estimated if reversible reactions between two enantiomers are not taken into account.[56] For the determination of an appropriate dosage regimen of a drug, it is therefore important to understand the characteristics of the reversible reaction between a drug and its metabolite.
In this tutorial, we focus on the basics of the pharmacokinetics of reversible metabolism. For a simple reversible reaction scheme, we use the matrix method to get exact solutions and plot them to find characteristic relationships between the pharmacokinetic profiles of a parent drug and its metabolites. In addition, we describe two approximation approaches to deal with complex problems and examine the conditions required for the successful application of each approach.

Theoretical analysis

Exact solution for the kinetics of reversible reaction

Let us consider the following simple reversible reaction scheme
Ak1k2Bk3C,tcp-27-52-e001.jpg
Where the first step is reversible with a forward rate constant k1 and reverse rate constant k2, and the second step is irreversible with the rate constant k3. The rate constant has the units of the amount per unit time, which are omitted for simplicity in this tutorial. The rate equations are:
(1)
dAdt=k1A+k2Btcp-27-52-e002.jpg
(2)
dBdt=k1Ak2+k3Btcp-27-52-e003.jpg
(3)
dCdt=k3Btcp-27-52-e004.jpg
Linear ordinary differential equations can be solved using Laplace transformation or the matrix method.[23712] We used the latter to obtain an exact solution because of its simplicity. The above expressions can be written in matrix form:
(4)
M=k1k20k1k2+k300k30tcp-27-52-e005.jpg
By setting
k1λk20k1k2+k3λ00k3λ=0tcp-27-52-e006.jpg, where λ is the eigenvalue, and solving for λ, we get the following characteristic polynomial:
(5)
λλ2+pλ+q=0tcp-27-52-e007.jpg
Where p = k1 + k2 + k3 and q = k1k3. The solutions of the quadratic equation in Eq. (5) are given by
α=p+p24q2,tcp-27-52-e008.jpg
and
β=pp24q2.tcp-27-52-e009.jpg
The resulting eigenvalues are 0, −α, and −β. These eigenvalues are substituted into the eigenvalue equation to yield the corresponding eigenvectors of
001,βk1k1β,andαk1k1αtcp-27-52-e010.jpg, spectively.
Because the matrix M has three distinct real eigenvalues and their corresponding eigenvectors, we can express the solution as
(6)
ABC=C0001+C1βk1k1βeαt+C2αk1k1αeβttcp-27-52-e011.jpg
At t = 0, A = A0 and B = C = 0, and A + B + C = A0 at all times. Using these boundary conditions, we can solve Eq. (6) to get C0 = A0, C1 = −C2 and C1 = A0/(β – α). Then, the final expressions for the above differential equations (1–3) are
(7)
A=A0βαβk1eαtαk1eβttcp-27-52-e012.jpg
(8)
B=k1A0βαeαteβttcp-27-52-e013.jpg
(9)
C=A0βαβα+αeβtβeαttcp-27-52-e014.jpg
We can compare these results to those for a consecutive two-step reaction, as described in the previous tutorial.[13] We can easily show that the latter has the eigenvalues of 0, −k1, and −k3, by setting k2 = 0 in Eq. (5). The solutions for B and C, equations 8 and 9 in this study, are very similar to those for the consecutive reaction (Equations 7 and 9 in Reference,[13] respectively). Thus, the expressions for the maximum amount of B (Bmax) and the time to reach the peak (tmax) can be obtained in a similar way:
(10)
tmax=1βαlnβαtcp-27-52-e015.jpg
and
(11)
Bmax=A0k1βαβαβα=A0k1ααβββα,tcp-27-52-e016.jpg
respectively.

Graphical insights into the kinetics of reversible reactions

Figure 1 shows representative amount versus time curves of A, B, and C, plotted on a linear and semi-logarithmic scale for three different cases: (a) k1 = 10, k2 = 0.1, and k3 = 1; (b) k1 = 5, k2 = 4, and k3 = 2; (c) k1 = 1, k2 = 5, and k3 = 0.2. When k1 and k3 are relatively larger than k2, that is, the reverse reaction is negligible, linear plots of the amount versus time curves of A and B for the reversible reaction (Fig. 1a) look very similar to those for a consecutive irreversible reaction (see Figure 2 in Reference[13]). The semi-logarithmic plots for the former case (Fig. 1b), however, look different from those for the latter case (see Figure 3a in Reference[13]). In the former case, two curves are always parallel to each other in the elimination phase, whereas they are not in the latter case. The difference can be explained as follows. For an irreversible reaction, the amount of A is described by one exponential term, whereas for the reversible reaction, it is described by two exponential terms. For the latter, the first term in Eq. (7) or (8) decreases faster than the second term over time, because α is always larger than β. This means that the elimination phase is governed by the second term. Therefore, the first term can be neglected when t is sufficiently large. In semi-logarithmic plots of A and B, thus, both slopes for the elimination phase are linear and identical and given by –β/2.303. The y-intercepts for A and B are αk1αβA0tcp-27-52-e017.jpg and k1αβA0tcp-27-52-e018.jpg, respectively. When 2k1 = 2k2 + k3, the intercepts are the same, as shown in Figure 2. If 2k1 < 2k2 + k3, the curves A and B will not cross each other. As a third case, we set k1 = 1, k2 = 5, and k3 = 0.2, which satisfies this condition. It is noteworthy that the semi-logarithmic plots for this case (Fig. 3b) strongly resemble those for the consecutive irreversible reaction (see Figure 3b in Reference[13]). In the latter case, the semi-logarithmic plots of plasma concentrations versus time of A and B become parallel straight lines at large t only when the formation rate of B is much slower than its elimination rate; that is, the metabolism is formation rate-limited (FRL).[131415] In addition, the intercept for [B] becomes greater than that for [A] when the distribution volume of B, VB, is much smaller than that of A, VA (see Figure 3c in Reference[13]). However, we need to be cautious when we apply this feature to the interpretation of the time profiles of the semi-logarithmic plots of a drug and its metabolite because this feature is commonly observed in reversible metabolic reactions. The pharmacokinetic parameters for each case, including tmax and Bmax, are summarized in Table 1.
The area under the curve (AUC) is an important pharmacokinetic parameter, representing total drug exposure over time. The AUC for A and B can be obtained by directly integrating Eqs. (7) and (8) from zero to infinity, respectively, and expressed as
(12)
AUCA=k2+k3A0qtcp-27-52-e019.jpg
and
(13)
AUCB=k1A0q.tcp-27-52-e020.jpg
Thus, the ratio of AUCA to AUCB can be simplified to k2+k3k1tcp-27-52-e021.jpg. The AUC for C can also be calculated by integrating Eq. (9) from zero to t and given by
(14)
AUCC=A0tpq.tcp-27-52-e022.jpg
Using a time integration method,[1416] we can obtain the relationship between AUC and rate constant or clearance, which has the units of volume (V) per unit time (CL = k × V). By separating the variables in Eqs. (7–9) and integrating both sides with respect to their respective variables, we can write the equations as follows:
(15)
VA[A]0dA=k1VA0Adt+k2VB0[B]dt,tcp-27-52-e023.jpg
(16)
VB[B]0dB=k1VA0Adtk2+k3VB0[B]dt,tcp-27-52-e024.jpg
(17)
VC[C]0dC=k3VB0[B]dt.tcp-27-52-e025.jpg
They should satisfy the following boundary conditions: [B] = [C] = 0 at t = 0; [A] = [B] = 0 and VC × [C] = A0 at t = ∞. By substituting k1 = CL1/VA, k2 = CL2/VB, and k3 = CL3/VC, we get
(18)
A0=CL1AUCACL2AUCB,
(19)
CL1AUCA=(CL2+CL3)AUCB
(20)
A0=CL3AUCB
Eq. (19) can be rearranged as
(21)
AUCBAUCA=CL1CL2+CL3.tcp-27-52-e026.jpg
These relationships are helpful in understanding pharmacokinetic features and metabolic pathways.[141718]

Approximation of the kinetics of reversible reactions

To get a good approximate solution for a reversible reaction scheme, approximation approaches have been developed and are widely used.[121920] Depending upon the reactivity of the intermediate species, we can use either of two approximation approaches. One is a steady-state approximation, in which the intermediate species B is very reactive and hence short lived. This approximation is useful for explaining reactions involving radical species. The other is an equilibrium approximation, in which B is quite stable and long-lived because the equilibrium reaction proceeds faster than the product formation reaction; that is, k1, k2 >> k3. For both approaches, the following condition should be satisfied:
(22)
dBdt0tcp-27-52-e027.jpg
In the steady-state approximation, B is extremely unstable and very short-lived and thus is present in most of the time in extremely small amounts. By assuming that B ≈ 0, we can get a new boundary condition: A + C = A0. From Eq. (22), Eq. (2) can be rearranged to yield:
(23)
BA=k1k2+k31tcp-27-52-e028.jpg
For the steady-state approximation, this ratio should be very small because the amount of B is much smaller than that of A. Substituting Eq. (23) into Eq. (1) and re-arranging it for A, we have
(24)
dAdt=k1k3k2+k3Atcp-27-52-e029.jpg
By separating the variables and integrating the resulting equation, we get
(25)
A=A0ek1k3k2+k3ttcp-27-52-e030.jpg
Using the above expression and the new boundary condition, A + C = A0, we obtain
(26)
C=A01ek1k3k2+k3ttcp-27-52-e031.jpg
Now we can solve Eq. (2) for B using Eq. (26), to get
(27)
B=A0k1k2+k3ek1k3k2+k3ttcp-27-52-e032.jpg
We plot the exact and approximate solutions for the case where k1 = 0.1, k2 = 1, and k3 = 5 (Fig. 4a). Because k1k2+k3=0.11+5~0.0171tcp-27-52-e033.jpg, the requirement for the steady-state approximation is satisfied. As the reaction begins, the amount of B quickly increases but is still extremely small (Bmax = 0.016 at tmax = 0.72) for most of the time, as expected. As shown in Figure 4a, the approximate solutions (dashed lines) agree well with the exact solutions (solid lines).
For the equilibrium approximation, the reaction needs to be in pseudo-equilibrium. In such a case, we can express
(28)
B=KeqA=k1k2A,tcp-27-52-e034.jpg
Where Keq is an equilibrium constant. Because
(29)
dAdt+dBdt+dCdt=0,tcp-27-52-e035.jpg
we also have
dAdt0ordCdt0.tcp-27-52-e036.jpg
Because of the above condition, it is not feasible to directly solve Eqs. (1) and (2). Instead, we need to combine the two equations to get
(30)
d(A+B)dt=k3B.tcp-27-52-e037.jpg
If we define X = A + B = A0C, from Eq. (28), we can express A and B as
(31)
A=k2k1+k2XandB=k1k1+k2X.tcp-27-52-e038.jpg
Now Eq. (30) can be solved for X to get
(32)
X=A0ek1k3k1+k2ttcp-27-52-e039.jpg
Then, we have
(33)
A=A0k2k1+k2ek1k3k1+k2t,tcp-27-52-e040.jpg
(34)
B=A0k1k1+k2ek1k3k1+k2t,tcp-27-52-e041.jpg
(35)
C=A0(1ek1k3k1+k2t).tcp-27-52-e042.jpg
In Figure 4b, we show the plots of the exact and approximate solutions for the case where k1 = 10, k2 = 1, and k3 = 0.1. Unlike in Figure 4a, the amount of B is not negligible (Bmax = 0.87 at tmax = 0.44). The reaction rapidly reaches quasi-equilibrium, after which the amount of B slowly decreases. As shown in Figure 4b, the approximate solutions (dashed lines) agree well with the exact solutions (solid lines).
The difference between steady-state and equilibrium approximate approaches lies in the range of their applicability. [20] When we introduce an additional constraint, k2 >> k1 into the equilibrium approach and k2 >> k3 into the steady state approach, both methods will yield the same solution.[20] This relationship can be further explained by considering an approximate expression for β. Using the following binomial theorem
1+xn=1+nx1!+nn1x22!+,tcp-27-52-e043.jpg
we can obtain the binomial expansion of β, which yields
(36)
β=k1+k2+k3(k1+k2+k3)24k1k32=k1k3k1+k2+k3+k1k32k1+k2+k33+2k1k33k1+k2+k35+.tcp-27-52-e044.jpg
If k1 + k2 >> k3 or k2 + k3 >> k1, the second and higher order terms in the expansion can be negligible, and the first term can be further simplified to be the same as the exponent in the equilibrium approximation [see Eqs. (33–35)] or in the steady-state approximation [see Eqs. (25–27)], respectively. For a refined approximation, we can differentiate Eq. (23) for A and B with respect to time t to get
(37)
dBdt=k1k2+k3dAdttcp-27-52-e045.jpg
By substituting Eqs. (3), (23), and (37) into Eq. (29) and rearranging the resulting equation, we get
(38)
dAdt=k1k3k1+k2+k3Atcp-27-52-e046.jpg
We expect a better approximate solution from the above equation compared with that produced by Eq. (24).

Application of the Reversible Reaction Scheme

The reversible reaction scheme is also used to describe drug distribution after rapid intravenous injection. Drug distribution within the body is often described by a two-compartment model with central and peripheral compartments.[21] In this model, an injected drug will be either eliminated or distributed within the central compartment, representing intravascular space, and the peripheral compartment, representing non-metabolizing body tissues. The distribution phase can be expressed by two linear differential equations similar to Eqs. 1 and 2.[21] The solution for drug concentration in each compartment can be obtained using a similar approach and described by two exponential terms.
To explain the kinetics of an enzyme-catalyzed reaction, a Michaelis-Menten model has been widely used.[222324] When an enzyme and a substrate form a complex, the binding of the substrate to the active site of the enzyme is non-covalent and reversible. If the conversion of the enzyme-substrate complex to a product is the rate-limiting process, the concentration of the complex can be considered to be constant. Then, the equilibrium approximation will give a hyperbolic relationship between the initial reaction rate, v_0, and substrate concentration [S], known as the Michaelis-Menten equation,
v0=Vmax[S]KS(orKm)+[S]tcp-27-52-e047.jpg
Where Vmax is the maximum reaction rate achieved when all of the active sites are occupied by substrates. The Michaelis constant, in the original article,[2223] is given by KS=k2k1tcp-27-52-e048.jpg. When a steady-state approximation is used, the constant is given by Km=k2+k3k1tcp-27-52-e049.jpg. [24]

Concluding remarks

We have described the mathematical approaches needed to determine the reaction kinetics of a reversible reaction system and provided graphical insights into the results. The knowledge acquired through this tutorial may help deepen the understanding of the pharmacokinetics of reversible metabolism and ultimately help with the prediction of the appropriate dose of a drug that is subject to equilibrium in a human body.

Acknowledgments

This study was supported by a grant of the Korean Health Technology R&D Project, Ministry of Health & Welfare, Republic of Korea (HI14C2750), as well as by the Industrial Core Technology Development Program (10051129, Development of the system for ADME assessment using radiolabeled compounds) funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea).

Notes

Reviewer: This article was invited and reviewed by the peer experts who are not TCP editors.

Conflict of interest: - Authors: The authors have no conflicts of interest to declare.

- Reviewers: Nothing to declare

- Editors: Nothing to declare

References

1. Lowry TM, John WT. CCLXIX. Studies of dynamic isomerism. Part XII. The equations for two consecutive unimolecular changes. J Chem Soc Trans. 1910; 97:2634–2645. DOI: 10.1039/CT9109702634.
2. Chrastil J. Determination of the first-order consecutive reversible reaction kinetics. Comput Chem. 1993; 17:103–106. DOI: 10.1016/0097-8485(93)80035-C.
crossref
3. Binti Hasan NAS, Balasubramanian P. Exact Solution for the Kinetic Equations of First Order Reversible Reaction Systems through Flow Graph Theory Approach. Ind Eng Chem Res. 2013; 52:10594–10600. DOI: 10.1021/ie303501t.
crossref
4. DiStefano JJ 3rd. Concepts, properties, measurements, and computation of clearance rates of hormones and other substances in biological systems. Ann Biomed Eng. 1976; 4:302–319. PMID: 984535.
5. Hwang S, Kwan KC, Albert KS. A linear mode of reversible metabolism and its application to bioavailability assessment. J Pharmacokinet Biopharm. 1981; 9:693–709. PMID: 7341755.
6. Cheng H, Jusko WJ. Pharmacokinetics of reversible metabolic systems. Biopharm Drug Dispos. 1993; 14:721–766. PMID: 8298069.
crossref
7. Hallak HO, Wedlund PJ. Reversible metabolism of vitamin K-vitamin K epoxide: modeling considerations and limitations. J Pharmacokinet Biopharm. 1992; 20:1–18. PMID: 1588501.
crossref
8. Eriksson T, Björkman S, Roth B, Höglund P. Intravenous formulations of the enantiomers of thalidomide: pharmacokinetic and initial pharmacodynamic characterization in man. J Pharm Pharmacol. 2000; 52:807–817. PMID: 10933131.
crossref
9. Ebling WF, Jusko WJ. The determination of essential clearance, volume, and residence time parameters of recirculating metabolic systems: the reversible metabolism of methylprednisolone and methylprednisone in rabbits. J Pharmacokinet Biopharm. 1986; 14:557–599. PMID: 3820091.
crossref
10. Baillie TA, Adams WJ, Kaiser DG, Olanoff LS, Halstead GW, Harpootlian H, et al. Mechanistic studies of the metabolic chiral inversion of (R)-ibuprofen in humans. J Pharmacol Exp Ther. 1989; 249:517–523. PMID: 2724138.
11. Stoschitzky K, Lindner W, Egginger G, Brunner F, Obermayer-Pietsch B, Passath A, et al. Racemic (R,S)-propranolol versus half-dosed optically pure (S)-propranolol in humans at steady state: Hemodynamic effects, plasma concentrations, and influence on thyroid hormone levels. Clin Pharmacol Ther. 1992; 51:445–453. PMID: 1563214.
crossref
12. Volk L, Richardson W, Lau KH, Hall M, Lin SH. Steady state and equilibrium approximations in reaction kinetics. J Chem Educ. 1977; 54:95–97. DOI: 10.1021/ed054p95.
crossref
13. Cho S, Yoon YR. Understanding the pharmacokinetics of prodrug and metabolite. Transl Clin Pharmacol. 2018; 26:1–5. DOI: 10.12793/tcp.2018.26.1.
crossref
14. Rowland M, Tozer TN. Clinical pharmacokinetics and pharmacodynamics: concepts and applications. 4th ed. Philadelphia: Wolters Kluwer Health/Lippincott William & Wilkins;2011. p. 603–632.
15. Smith PC. Pharmacokinetics of drug metabolites. In : Pearson PG, Wienkers LC, editors. Handbook of drug metabolism. 2nd ed. New York: Informa Healthcare;2009. p. 17–59.
16. Mucientes AE, de la Peña MA. Kinetic Analysis of Parallel-Consecutive First-Order Reactions with a Reversible Step: Concentration–Time Integrals Method. J Chem Educ. 2009; 86:390–392. DOI: 10.1021/ed086p390.
crossref
17. Pang KS. A review of metabolite kinetics. J Pharmacokinet Biopharm. 1985; 13:633–662. PMID: 3914545.
crossref
18. de Campos ML, Padilha EC, Peccinini RG. A review of pharmacokinetic parameters of metabolites and prodrugs. Drug Metab Lett. 2014; 7:105–116. PMID: 24628402.
crossref
19. McDaniel DH, Smoot CR. Approximations in the kinetics of consecutive reactions. J Phys Chem. 1956; 60:966–969. DOI: 10.1021/j150541a035.
crossref
20. Pyun CW. Steady-state and equilibrium approximations in chemical kinetics. J Chem Educ. 1971; 48:194–196. DOI: 10.1021/ed048p194.
crossref
21. Atkinson Jr AJ. Compartment analysis of drug distribution. In : Atkinson AJ, Huang S-M, Lertora JJL, Markey SP, editors. Principles of Clinical Pharmacology. 3rd ed. Amsterdam: Academic Press;2012. p. 27–39.
22. Michaelis L, Menten ML. Die Kinetik der Invertinwirkung. Biochem Z. 1913; 49:333–369.
23. Michaelis L, Menten ML, Johnson KA, Goody RS. The Original Michaelis Constant: Translation of the 1913 Michaelis–Menten Paper. Biochemistry. 2011; 50:8264–8269. DOI: 10.1021/bi201284u. PMID: 21888353.
crossref
24. Briggs GE, Haldane JB. A note on the kinetics of enzyme action. Biochem J. 1925; 19:338–339. PMID: 16743508.
crossref
Figure 1

(a) Linear and (b) semi-logarithmic plots of the amount-time curves of A, B, or C in the reversible reaction, Ak1k2Bk3Ctcp-27-52-e050.jpg, where k1 = 10, k2 = 0.1, and k3 = 1. The abscissa and ordinate denote the time and the amount of each species, respectively, in arbitrary unit. See text for details.

tcp-27-52-g001
Figure 2

(a) Linear and (b) semi-logarithmic plots of the amount-time curves of A, B, or C in the reversible reaction, Ak1k2Bk3Ctcp-27-52-e051.jpg, where k1 = 5, k2 = 4, and k3 = 2. The abscissa and ordinate denote time and the amount of each species, respectively, in arbitrary unit. See text for details.

tcp-27-52-g002
Figure 3

(a) Linear and (b) semi-logarithmic plots of the amount-time curves of A, B, or C in the reversible reaction, Ak1k2Bk3Ctcp-27-52-e052.jpg, where k1 = 1, k2 = 5, and k3 = 0.2. The abscissa and ordinate denote time and the amount of each species, respectively, in arbitrary unit. See text for details.

tcp-27-52-g003
Figure 4

Linear plots of the amount-time curves of A, B, or C for the cases of (a) k1 = 0.1, k2 = 1, and k3 = 5, corresponding to the steady-state approximation, and (b) k1 = 10, k2 = 1, and k3 = 0.1, corresponding to the equilibrium approximation. The solid and dashed lines represent the exact and approximate solutions, respectively. The abscissa and ordinate denote time and the amount of each species, respectively, in arbitrary unit. See text for details.

tcp-27-52-g004
Table 1

Pharmacokinetic parameters for the following cases: (a) k1 = 10, k2 = 0.1, and k3 = 1 (as in Fig. 1); (b) k1 = 5, k2 = 4, and k3 = 2 (as in Figure 2); (c) k1 = 1, k2 = 5, and k3 = 0.2 (as in Figure 3); (d) k1 = 0.1, k2 = 1, and k3 = 5, (as in Figure 4a); (e) k1 = 10, k2 = 1, and k3 = 0.1 (as in Figure 4b)

tcp-27-52-i001
(a) (b) (c) (d) (e)
tmax (B) 0.245 0.256 0.855 0.722 0.439
Bmax 0.768 0.387 0.158 0.0157 0.873
AUC0–∞ (A) 0.11 0.6 26 12.0 1.1
AUC0–∞ (B) 1.00 0.5 5 0.2 10.0
AUC0–100 (C) 98.89 98.9 69 87.8 88.9
Intercept (A)* 0.0122 0.556 0.842 0.998 0.0924
Intercept (B)* 1.10 0.556 0.163 0.0169 0.915
Slope (A & B)* 0.429 0.434 0.014 0.0361 0.0394

AUC0–∞ and AUC0–100 represent the areas under the curve from t = 0 to t = ∞ and 100, respectively. All the intercepts and slopes (*) were obtained from semi-logarithmic plots.

TOOLS
Similar articles