Predicting Human Dermal Drug Concentrations Using PBPK Modeling and Simulation: Clobetasol Propionate Case Study

Quantitative in silico tools may be leveraged to mechanistically predict the dermato-pharmacokinetics of compounds delivered from topical and transdermal formulations by integrating systems of rate equations that describe permeation through the formulation and layers of skin and pilo-sebaceous unit, and exchange with systemic circulation via local blood flow. Delivery of clobetasol-17 propionate (CP) from DermovateTM cream was simulated using the Transdermal Compartmental Absorption & Transit (TCATTM) Model in GastroPlus®. The cream was treated as an oil-in-water emulsion, with model input parameters estimated from publicly available information and quantitative structure-permeation relationships. From the ranges of values available for model input parameters, a set of parameters was selected by comparing model outputs to CP dermis concentration-time profiles measured by dermal open-flow microperfusion (Bodenlenz et al. Pharm Res. 33(9):2229-38, 2016). Predictions of unbound dermis CP concentrations were reasonably accurate with respect to time and skin depth. Parameter sensitivity analyses revealed considerable dependence of dermis CP concentration profiles on drug solubility in the emulsion, relatively less dependence on dispersed phase volume fraction and CP effective diffusivity in the continuous phase of the emulsion, and negligible dependence on dispersed phase droplet size. Effects of evaporative water loss from the cream and corticosteroid-induced vasoconstriction were also assessed. This work illustrates the applicability of computational modeling to predict sensitivity of dermato-pharmacokinetics to changes in thermodynamic and transport properties of a compound in a topical formulation, particularly in relation to rate-limiting steps in skin permeation. Where these properties can be related to formulation composition and processing, such a computational approach may support the design of topically applied formulations.

Modeling the percutaneous absorption of solvent-deposited solids over a wide dose range: II. Weak electrolytes

Dermal absorption of weak electrolytes applied to skin from pharmaceutical and cosmetic compositions is an important consideration for both their efficacy and skin safety. We developed a mechanistic, physics-based framework that simulates this process for leave on applications following solvent deposition. We incorporated this framework into our finite dose computational skin permeation model previously tested with nonelectrolytes to generate quantitative predictions for weak electrolytes. To test the model, we analyzed experimental data from an in vitro human skin permeation study of a weak acid (benzoic acid) and a weak base (propranolol) and their sodium and hydrochloride salts from simple, ethanol/water vehicles as a function of dose and ionization state. Key factors controlling absorption are the pH and buffer capacity of the dose solution, the dissolution rate of precipitated solids into a lipid boundary layer and the rate of conversion of the deposited solid to its conjugate form as the nonionized component permeates and (sometimes) evaporates from the skin surface. The resulting framework not only describes the current test data but has the potential to predict the absorption of other weak electrolytes following topical application.

Impact of solvent dry down, vehicle pH and slowly reversible keratin binding on skin penetration of cosmetic relevant compounds: I. Liquids

To measure progress and evaluate performance of the newest UB/UC/P&G skin penetration model we simulated an 18-compound subset of finite dose in vitro human skin permeation data taken from a solvent-deposition study of cosmetic-relevant compounds (Hewitt et al., J. Appl. Toxicol. 2019, 1-13). The recent model extension involved slowly reversible binding of solutes to stratum corneum keratins. The selected subset was compounds that are liquid at skin temperature. This set was chosen to distinguish between slow binding and slow dissolution effects that impact solid phase compounds. To adequately simulate the physical experiments there was a need to adjust the evaporation mass transfer coefficient to better represent the diffusion cell system employed in the study. After this adjustment the model successfully predicted both dermal delivery and skin surface distribution of 12 of the 18 compounds. Exceptions involved compounds that were cysteine-reactive, highly water-soluble or highly ionized in the dose solution. Slow binding to keratin, as presently parameterized, was shown to significantly modify the stratum corneum kinetics and diffusion lag times, but not the ultimate disposition, of the more lipophilic compounds in the dataset. Recommendations for further improvement of both modeling methods and experimental design are offered.

Modeling the Percutaneous Absorption of Solvent-deposited Solids Over a Wide Dose Range

The transient absorption of two skin care agents, niacinamide (nicotinamide, NA) and methyl nicotinate (MN), solvent-deposited on ex vivo human skin mounted in Franz diffusion cells has been analyzed according to a new variation on a recently published mechanistic skin permeability model (Yu et al. 2020. J Pharm Sci 110:2149-56). The model follows the absorption and evaporation of two components, solute and solvent, and it includes both a follicular transport component and a dissolution rate limitation for high melting, hydrophilic solids deposited on the skin. Explicit algorithms for improving the simulation of transient diffusion of solvent-deposited solids are introduced. The simulations can account for the ex vivo skin permeation time course of both NA and MN over a dose range exceeding 4.5 orders of magnitude. The model allows one to describe on a mechanistic basis why the percutaneous absorption rate of NA is approximately 60-fold lower than that of its lower melting, more lipophilic analog, MN. It furthermore suggests that MN perturbs stratum corneum barrier lipids and increases their permeability while NA does not, presenting a challenge to molecular modelers engaged in simulating biological lipid barriers.

Computer Simulation of Skin Permeability of Hydrophobic and Hydrophilic Chemicals - Influence of Follicular Pathway

A recently published mechanistic skin permeability model (Kasting et al., 2019. J Pharm Sci 108:337-349) that included a follicular diffusion pathway has been extended to describe transient diffusion and finite dose applications. The model follows the disposition of two components, solute and solvent, so that solvent deposition processes can be explicitly represented. Experimentally-calibrated permeability characteristics of the follicular pathway leading to the permeation of highly hydrophilic permeants are further refined. Details of the refinements and a comparison with the earlier model using two large experimental datasets are presented. An example calculation shows the marked difference between the time scales for achievement of near steady-state diffusion for large hydrophilic and lipophilic compounds, with the former being more than 100-fold faster than the latter. However, the true steady state for the hydrophilic compound is not reached until much later due to the very slow filling of the corneocyte phase.