Physiologically based pharmacokinetic modelling

Physiologically based pharmacokinetic (PBPK) modeling is a mathematical modeling technique for predicting the absorption, distribution, metabolism and excretion (ADME) of synthetic or natural chemical substances in humans and other animal species. PBPK modeling is used in pharmaceutical research and drug development, and in health risk assessment for cosmetics or general chemicals. d Q i d t = F i ( C a r t − Q i P i V i ) {displaystyle {dQ_{i} over dt}=F_{i}(C_{art}-{{Q_{i}} over {P_{i}V_{i}}})} d Q g d t = F g ( C a r t − Q g P g V g ) + K a Q i n g {displaystyle {dQ_{g} over dt}=F_{g}(C_{art}-{{Q_{g}} over {P_{g}V_{g}}})+K_{a}Q_{ing}} d Q i n g d t = − K a Q i n g {displaystyle {dQ_{ing} over dt}=-K_{a}Q_{ing}} d Q g d t = F g ( C a r t − Q g P g V g ) + R i n g {displaystyle {dQ_{g} over dt}=F_{g}(C_{art}-{{Q_{g}} over {P_{g}V_{g}}})+R_{ing}} d Q s c u d t = K p × S s × ( 1 − f S e ) × ( Q s u P s c V s c u − C s c u ) {displaystyle {dQ_{{sc}_{u}} over dt}=K_{p} imes S_{s} imes (1-f_{S_{e}}) imes ({Q_{s_{u}} over {P_{sc}V_{{sc}_{u}}}}-C_{{sc}_{u}})} d Q s u d t = F s ( 1 − f S e ) ( C a r t − Q s u P s V s u ) − d Q s c u d t {displaystyle {dQ_{s_{u}} over dt}=F_{s}(1-f_{S_{e}})(C_{art}-{{Q_{s_{u}}} over {P_{s}V_{s_{u}}}})-{dQ_{{sc}_{u}} over dt}} d Q s c e d t = K p × S s × f S e × ( Q s e P s c V s c e − C s c e ) {displaystyle {dQ_{{sc}_{e}} over dt}=K_{p} imes S_{s} imes f_{S_{e}} imes ({Q_{s_{e}} over {P_{sc}V_{{sc}_{e}}}}-C_{{sc}_{e}})} d Q s e d t = F s f S e ( C a r t − Q s e P s V s e ) − d Q s c e d t {displaystyle {dQ_{s_{e}} over dt}=F_{s}f_{S_{e}}(C_{art}-{{Q_{s_{e}}} over {P_{s}V_{s_{e}}}})-{dQ_{{sc}_{e}} over dt}} Physiologically based pharmacokinetic (PBPK) modeling is a mathematical modeling technique for predicting the absorption, distribution, metabolism and excretion (ADME) of synthetic or natural chemical substances in humans and other animal species. PBPK modeling is used in pharmaceutical research and drug development, and in health risk assessment for cosmetics or general chemicals. PBPK models strive to be mechanistic by mathematically transcribing anatomical, physiological, physical, and chemical descriptions of the phenomena involved in the complex ADME processes. A large degree of residual simplification and empiricism is still present in those models, but they have an extended domain of applicability compared to that of classical, empirical function based, pharmacokinetic models. PBPK models may have purely predictive uses, but other uses, such as statistical inference, have been made possible by the development of Bayesian statistical tools able to deal with complex models. That is true for both toxicity risk assessment and therapeutic drug development. PBPK models try to rely a priori on the anatomical and physiological structure of the body, and to a certain extent, on biochemistry. They are usually multi-compartment models, with compartments corresponding to predefined organs or tissues, with interconnections corresponding to blood or lymph flows (more rarely to diffusions). A system of differential equations for concentration or quantity of substance on each compartment can be written, and its parameters represent blood flows, pulmonary ventilation rate, organ volumes etc., for which information is available in scientific publications. Indeed, the description they make of the body is simplified and a balance needs to be struck between complexity and simplicity. Besides the advantage of allowing the recruitment of a priori information about parameter values, these models also facilitate inter-species transpositions or extrapolation from one mode of administration to another (e.g., inhalation to oral). An example of a 7-compartment PBPK model, suitable to describe the fate of many solvents in the mammalian body, is given in the Figure on the right. The first pharmacokinetic model described in the scientific literaturewas in fact a PBPK model. It led, however, to computations intractable at that time. The focus shifted then to simpler models,for which analytical solutions could be obtained (such solutions were sums of exponential terms, which led to further simplifications.) The availability of computers and numerical integration algorithms marked a renewed interest in physiological models in the early 1970s.For substances with complex kinetics, or when inter-species extrapolations were required, simple models were insufficient and research continued on physiological models.By 2010, hundreds of scientific publications have described and used PBPK models, and at least two private companies are basing their business on their expertise in this area. The model equations follow the principles of mass transport, fluid dynamics, and biochemistry in order to simulate the fate of a substance in the body.Compartments are usually defined by grouping organs or tissues with similar blood perfusion rate and lipid content (i.e. organs for which chemicals' concentration vs. time profiles will be similar). Ports of entry (lung, skin, intestinal tract...), ports of exit (kidney, liver...) and target organs for therapeutic effect or toxicity are often left separate. Bone can be excluded from the model if the substance of interest does not distribute to it. Connections between compartment follow physiology (e.g., blood flow in exit of the gut goes to liver, etc.) Drug distribution into a tissue can be rate-limited by either perfusion or permeability. Perfusion-rate-limited kinetics apply when the tissue membranes present no barrier to diffusion. Blood flow, assuming that the drug is transported mainly by blood, as is often the case, is then the limiting factor to distribution in the various cells of the body. That is usually true for small lipophilic drugs. Under perfusion limitation, the instantaneous rate of entry for the quantity of drug in a compartment is simply equal to (blood) volumetric flow rate through the organ times the incoming blood concentration. In that case; for a generic compartment i, the differential equation for the quantity Qi of substance, which defines the rate of change in this quantity, is: where Fi is blood flow (noted Q in the Figure above), Cart incoming arterial blood concentration, Pi the tissue over blood partition coefficient and Vi the volume of compartment i.