A geometrical model for diffusion of hydrophilic compounds in human stratum corneum

Permeable Cornified Envelope Layer Regulates the Solute Transport in Human Stratum Corneum

To unravel the diffusion mechanisms of percutaneous drug delivery, suitable numerical analysis of stratum corneum structure is essential. In this research paper, we accounted for the permeable envelope layer in the brick-and-mortar finite element models of human stratum corneum. Both penetration and desorption experiments for tritiated water were simulated by transient finite element analysis. Rivet-shaped corneodesmosomes were included in the brick and mortar model. Results showed that cornified lipid permeability (Penv) is a determinant in desorption of the solute, while lipid transverse diffusion coefficient (Dlip-trans) is prominent during penetration. These two major unknowns (Penv and Dlip-trans) were obtained by extensive fitting of the finite element model to the experimental water data. Penv and Dlip-trans were determined to be 1×10-2 cm/s and 5.7×10-10 cm2/s, respectively.

Membrane permeability based on mesh analysis

The main function of a membrane is to control the exchange of matter between the surrounding regions. As such, accurate modeling of membranes is important to properly describe their properties. In many cases in both biological systems and technical applications, the membranes are composite structures where transport properties may vary between the different sub-regions of the membrane. In this work we develop a method based on Mesh analysis that is asymptotically exact and can describe diffusion in composite membrane structures. We do this by first reformulating a generalized Fick's law to include the effects from activity coefficient, diffusion coefficient, and solubility using a single condensed parameter. We then use the derived theory and Mesh analysis to, in essence, retrieve a finite element method approach. The calculated examples are based on a membrane structure that reassembles that of the brick and mortar structure of stratum corneum, the upper layer of our skin. Resulting concentration profiles from this procedure are then compared to experimental results for the distribution of different probes within intact stratum corneum, showing good agreement. Based on the derived approach we further investigate the impact from a gradient in the fluidity of the stratum corneum mortar lipids across the membrane, and find that it is substantial. We also show that anisotropic organisation of the lipid mortar can have large impact on the effective permeability compared to isotropic mortar lipids. Finally, we examine the effects of corneocyte swelling, and their lateral arrangement in the membrane on the overall membrane permeability.

Transdermal and lateral effective diffusivities for drug transport in stratum corneum from a microscopic anisotropic diffusion model

This paper presents a computational model of molecular diffusion through the interfollicular stratum corneum. Specifically, it extends an earlier two-dimensional microscopic model for the permeability in two ways: (1) a microporous leakage pathway through the intercellular lipid lamellae allows slow permeation of highly hydrophilic permeants through the tissue; and (2) the model yields explicit predictions of both lateral (D‾_‖^sc) and transdermal (D‾_⊥^sc) effective (average, homogenized) diffusivities of solutes within the tissue. We present here the mathematical framework for the analysis and a comparison of the predictions with experimental data on desorption of both hydrophilic and lipophilic solutes from human stratum corneum in vitro. Diffusion in the lipid lamellae is found to make the effective diffusivity highly anisotropic, with the predicted ratio D‾_‖^sc/D‾_⊥^sc ranging from 34-39 for fully hydrated skin and 150 to more than 1000 for partially hydrated skin. The diffusivities and their ratio are in accord with both experimental data and the results of mathematical analyses performed by others.

Compartmental modeling of skin absorption and desorption kinetics: Donor solvent evaporation, variable diffusion/partition coefficients, and slow equilibration process within stratum corneum

This work expands the recently developed compartmental model for skin transport to model variable diffusion and/or partition coefficients, and the presence of slow equilibration/slow binding kinetics within stratum corneum. The model was validated by comparing it with the diffusion model which was solved numerically using the finite element method. It was found that the new compartmental model predictions agreed well with that of the diffusion model, providing a sufficient number of compartments was used. The compartmental model was applied to two previously published experimental data sets: water penetration and desorption data and the finite dose dermal penetration of testosterone. Significant improvement of the fitting quality for all these data sets was achieved using the compartmental model.

Relating transdermal delivery plasma pharmacokinetics with in vitro permeation test (IVPT) findings using diffusion and compartment-in-series models

Increasing emphasis is being placed on using in vitro permeation test (IVPT) results for topical products as a surrogate for their in vivo behaviour. This study sought to relate in vivo plasma concentration - time pharmacokinetic (PK) profiles after topical application of drug products to IVPT findings with mechanistic diffusion and compartment models that are now widely used to describe permeation of solutes across the main skin transport barrier, the stratum corneum. Novel in vivo forms of the diffusion and compartment-in-series models were developed by combining their IVPT model forms with appropriate in vivo disposition functions. Available in vivo and IVPT data were then used with the models in data analyses, including the estimation of prediction intervals for in vivo plasma concentrations derived from IVPT data. The resulting predicted in vivo plasma concentration - time profiles for the full models corresponded closely with the observed results for both nitroglycerin and rivastigmine at all times. In contrast, reduced forms of these in vivo models led to discrepancies between model predictions and observed results at early times. A two-stage deconvolution procedure was also used to estimate the in vivo cumulative amount absorbed and shown to be linearly related to that from IVPT, with an acceptable prediction error. External predictability was also shown using a separate set of in vitro and in vivo data for different nitroglycerin patches. This work suggests that mechanistic and physiologically based pharmacokinetic models can be used to predict in vivo behaviour from IVPT data for topical products.