Determining the Effect of pH on the Partitioning of Neutral, Cationic and Anionic Chemicals to Artificial Sebum: New Physicochemical Insight and QSPR Model

In Silico Study on the Contribution of the Follicular Route to Dermal Permeability of Small Molecules

PURPOSE

This study investigates in silico the contribution of the hair follicle to the overall dermal permeability of small molecules, as published experimental work provides inconclusive information on whether the follicular route favours the permeation of hydrophobic or hydrophilic permeants.

METHOD

A study is conducted varying physico-chemical parameters of permeants such as lipophilicity, molecular weight and protein binding. The simulated data is compared to published experimental data to discuss how those properties can modulate the contribution of the hair follicle to the overall dermal permeation.

RESULTS

The results indicate that the contribution of the follicular route to dermal permeation can range from negligible to notable depending on the combination of lipophilic/hydrophilic properties of the substance filling the follicular route and the permeant.

CONCLUSION

Characterisation of the substance filling the follicular route is required for analysing the experimental data of dermal permeation of small molecules, as changes between in vivo and in vitro due to handling of samples and cessation of vital functions can modify the contribution of the follicular route to overall dermal permeation, hence hindering data interpretation.

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.

Mechanistic Modeling of In Vitro Skin Permeation and Extrapolation to In Vivo for Topically Applied Metronidazole Drug Products Using a Physiologically Based Pharmacokinetic Model

Physiologically based pharmacokinetic (PBPK) modeling has increasingly been employed in dermal drug development and regulatory assessment, providing a framework to integrate relevant information including drug and drug product attributes, skin physiology parameters, and population variability. The current study aimed to develop a stepwise modeling workflow with knowledge gained from modeling in vitro skin permeation testing (IVPT) to describe in vivo exposure of metronidazole locally in the stratum corneum following topical application of complex semisolid drug products. The initial PBPK model of metronidazole in vitro skin permeation was developed using infinite and finite dose aqueous metronidazole solution. Parameters such as stratum corneum lipid-water partition coefficient (K_sclip/water) and stratum corneum lipid diffusion coefficient (D_sclip) of metronidazole were optimized using IVPT data from simple aqueous solutions (infinite) and MetroGel (10 mg/cm² dose application), respectively. The optimized model, when parameterized with physical and structural characteristics of the drug products, was able to accurately predict the mean cumulative amount permeated (cm²/h) and flux (μg/cm²/h) profiles of metronidazole following application of different doses of MetroGel and MetroCream. Thus, the model was able to capture the impact of differences in drug product microstructure and metamorphosis of the dosage form on in vitro metronidazole permeation. The PBPK model informed by IVPT study data was able to predict the metronidazole amount in the stratum corneum as reported in clinical studies. In summary, the proposed model provides an enhanced understanding of the potential impact of drug product attributes in influencing in vitro skin permeation of metronidazole. Key kinetic parameters derived from modeling the metronidazole IVPT data improved the predictions of the developed PBPK model of in vivo local metronidazole concentrations in the stratum corneum. Overall, this work improves our confidence in the proposed workflow that accounts for drug product attributes and utilizes IVPT data toward improving predictions from advanced modeling and simulation tools.

In silico predictions of absorption of MDI substances after dermal or inhalation exposures to support a category based read-across assessment

Methylenediphenyl diisocyanate (MDI) substances used polyurethane production can range from their simplest monomeric forms (e.g., 4,4'-MDI) to mixtures of the monomers with various homologues, homopolymer, and prepolymer derivatives. The relative dermal or inhalation absorption of 39 constituents of these substances in human were predicted using the GastroPlus® program. Predicted dermal uptake and absorption of the three MDI monomers from an acetone vehicle was 84-86% and 1.4-1.5%, respectively, with lower uptake and absorption predicted for the higher MW analogs. Lower absorption was predicted from exposures in a more lipophilic vehicle (1-octanol). Modeled inhalation exposures afforded the highest pulmonary absorption for the MDI monomers (38-54%), with 3-27% for the MW range of 381-751, and <0.1% for the remaining, higher MW derivatives. Predicted oral absorption, representing mucociliary transport, ranged from 5 to 10% for the MDI monomers, 10-25% for constituents of MW 381-751, and ≤3% for constituents with MW > 900. These in silico evaluations should be useful in category-based, worst-case, Read-Across assessments for MDI monomers and modified MDI substances for potential systemic effects. Predictions of appreciable mucociliary transport may also be useful to address data gaps in oral toxicity testing for this category of compounds.

Determination of Solute Diffusion Properties in Artificial Sebum

Despite a number of studies showed that hair follicular pathway contributed significantly to transdermal delivery, there have been limited studies on the diffusion properties of chemicals in sebum. Here, the diffusion property of 17 chemical compounds across artificial sebum has been measured using diffusion cell. The diffusion flux showed 2 types of distinctive behaviors: that reached steady state and that did not. Mathematical models have been developed to fit the experimental data and derive the sebum diffusion and partition coefficients. The models considered the uneven thickness of the sebum film and the additional resistance of the unstirred aqueous boundary layer and the supporting filter. The derived sebum-water partition coefficients agreed well with the experimental data measured previously using equilibrium depletion method. The obtained diffusion coefficients in artificial sebum only depended on the molecular size. Change in pH for ionic chemicals did not affect the diffusion coefficients but influenced their diffusion flux because of the change of sebum-water partition coefficients. Generally, the measured diffusion coefficients of chemicals in artificial sebum are about one order of magnitude higher than those in the stratum corneum lipids, suggesting the hair follicle might have a non-negligible contribution to the overall permeation.