A Two-Phase Analysis of Solute Partitioning into the Stratum Corneum

Measurement of phthalates in skin wipes: estimating exposure from dermal absorption

This study has determined the levels of six phthalates (dimethyl phthalate (DMP), diethyl phthalate (DEP), di(isobutyl) phthalate (DiBP), di(n-butyl) phthalate (DnBP), butyl benzyl phthalate (BBzP), and di(2-ethylhexyl) phthalate (DEHP)) in skin wipes; examined factors that might influence the levels, including body location, time of sampling, and hand-washing; and estimated dermal absorption based on the measured levels. Skin wipes were collected from the forehead, forearm, back-of-hand, and palm of 20 participants using gauze pads moistened with isopropanol. DiBP, DnBP, and DEHP were most frequently detected; DEHP levels were substantially higher than DnBP and DiBP levels, and DnBP levels were somewhat lower than DiBP levels. The levels differed at different body locations, with palm > back-of-hand > forearm ≥ forehead. Repeated wipe sampling from six participants over a 1 month period indicated that levels at the same body location did not vary significantly. The estimated median total dermal absorption from skin surface lipids on the palm, back-of-hand, arm, and head are 0.48, 0.68, and 0.66 (μg/kg)/day for DiBP, DnBP, and DEHP, respectively. These estimates are roughly 10-20% of the total uptake reported for Chinese adults and suggest that dermal absorption contributes significantly to the uptake of these phthalates. Washing with soap and water removed more than 50% of the phthalates on the hands and may be a useful tool in decreasing aggregate phthalate exposure.

Predicting dermal absorption of gas-phase chemicals: transient model development, evaluation, and application

A transient model is developed to predict dermal absorption of gas-phase chemicals via direct air-to-skin-to-blood transport under non-steady-state conditions. It differs from published models in that it considers convective mass-transfer resistance in the boundary layer of air adjacent to the skin. Results calculated with this transient model are in good agreement with the limited experimental results that are available for comparison. The sensitivity of the modeled estimates to key parameters is examined. The model is then used to estimate air-to-skin-to-blood absorption of six phthalate esters for scenarios in which (A) a previously unexposed occupant encounters gas-phase phthalates in three different environments over a single 24-h period; (B) the same as 'A', but the pattern is repeated for seven consecutive days. In the 24-h scenario, the transient model predicts more phthalate absorbed into skin and less absorbed into blood than would a steady-state model. In the 7-day scenario, results calculated by the transient and steady-state models converge over a time period that varies between 3 and 4 days for all but the largest phthalate (DEHP). Dermal intake is comparable to or larger than inhalation intake for DEP, DiBP, DnBP, and BBzP in Scenario 'A' and for all six phthalates in Scenario 'B'.

PRACTICAL IMPLICATIONS

Dermal absorption from air has often been overlooked in exposure assessments. However, our transient model suggests that dermal intake of certain gas-phase phthalate esters is comparable to, or larger than, inhalation intake under commonly occurring indoor conditions. This may also be the case for other organic chemicals that have physicochemical properties that favor dermal absorption directly from air. Consequently, this pathway should be included in aggregate exposure and risk assessments. Furthermore, under conditions where the exposure concentrations are changing or there is insufficient time to achieve steady-state, the transient model presented in this study is more appropriate for estimating dermal absorption than is a steady-state model.

Dermic diffusion and stratum corneum: a state of the art review of mathematical models

Transdermal biotechnologies are an ever increasing field of interest, due to the medical and pharmaceutical applications that they underlie. There are several mathematical models at use that permit a more inclusive vision of pure experimental data and even allow practical extrapolation for new dermal diffusion methodologies. However, they grasp a complex variety of theories and assumptions that allocate their use for specific situations. Models based on Fick's First Law found better use in contexts where scaled particle theory Models would be extensive in time-span but the reciprocal is also true, as context of transdermal diffusion of particular active compounds changes. This article reviews extensively the various theoretical methodologies for studying dermic diffusion in the rate limiting dermic barrier, the stratum corneum, and systematizes its characteristics, their proper context of application, advantages and limitations, as well as future perspectives.

A personal view of skin permeation (1960-2013)

I want to thank Mike Roberts for arranging this issue of the journal and Annette Bunge and Richard Guy for reviewing this paper. This is, first, a history of my introduction to the field of skin permeability, and then an attempt to recall (with the inaccuracies that implies) the highlights of my research (1960-1975) on skin permeation, and finally a reflection by an old-timer on more recent concepts.

Dermal uptake of 18 dilute aqueous chemicals: in vivo disappearance-method measures greatly exceed in vitro-based predictions

Average rates of total dermal uptake (Kup ) from short-term (e.g., bathing) contact with dilute aqueous organic chemicals (DAOCs) are typically estimated from steady-state in vitro diffusion-cell measures of chemical permeability (Kp ) through skin into receptor solution. Widely used ("PCR-vitro") methods estimate Kup by applying diffusion theory to increase Kp predictions made by a physico-chemical regression (PCR) model that was fit to a large set of Kp measures. Here, Kup predictions for 18 DAOCs made by three PCR-vitro models (EPA, NIOSH, and MH) were compared to previous in vivo measures obtained by methods unlikely to underestimate Kup . A new PCR model fit to all 18 measures is accurate to within approximately threefold (r = 0.91, p < 10(-5) ), but the PCR-vitro predictions (r > 0.63) all tend to underestimate the Kup measures by mean factors (UF, and p value for testing UF = 1) of 10 (EPA, p < 10(-6) ), 11 (NIOSH, p < 10(-8) ), and 6.2 (MH, p = 0.018). For all three PCR-vitro models, log(UF) correlates negatively with molecular weight (r(2) = 0.31 to 0.84, p = 0.017 to < 10(-6) ) but not with log(vapor pressure) as an additional predictor (p > 0.05), so vapor pressure appears not to explain the significant in vivo/PCR-vitro discrepancy. Until this discrepancy is explained, careful in vivo measures of Kup should be obtained for more chemicals, the expanded in vivo database should be compared to in vitro-based predictions, and in vivo data should be considered in assessing aqueous dermal exposure and its uncertainty.

Modeling the human skin barrier--towards a better understanding of dermal absorption

Many drugs are presently delivered through the skin from products developed for topical and transdermal applications. Underpinning these technologies are the interactions between the drug, product and skin that define drug penetration, distribution, and elimination in and through the skin. Most work has been focused on modeling transport of drugs through the stratum corneum, the outermost skin layer widely recognized as presenting the rate-determining step for the penetration of most compounds. However, a growing body of literature is dedicated to considering the influence of the rest of the skin on drug penetration and distribution. In this article we review how our understanding of skin physiology and the experimentally observed mechanisms of transdermal drug transport inform the current models of drug penetration and distribution in the skin. Our focus is on models that have been developed to describe particular phenomena observed at particular sites of the skin, reflecting the most recent directions of investigation.

Recent advances in predicting skin permeability of hydrophilic solutes

Understanding the permeation of hydrophilic molecules is of relevance to many applications including transdermal drug delivery, skin care as well as risk assessment of occupational, environmental, or consumer exposure. This paper reviews recent advances in modeling skin permeability of hydrophilic solutes, including quantitative structure-permeability relationships (QSPR) and mechanistic models. A dataset of measured human skin permeability of hydrophilic and low hydrophobic solutes has been compiled. Generally statistically derived QSPR models under-estimate skin permeability of hydrophilic solutes. On the other hand, including additional aqueous pathway is necessary for mechanistic models to improve the prediction of skin permeability of hydrophilic solutes, especially for highly hydrophilic solutes. A consensus yet has to be reached as to how the aqueous pathway should be modeled. Nevertheless it is shown that the contribution of aqueous pathway can constitute to more than 95% of the overall skin permeability. Finally, future prospects and needs in improving the prediction of skin permeability of hydrophilic solutes are discussed.

Design and performance of a spreadsheet-based model for estimating bioavailability of chemicals from dermal exposure

A comprehensive transient model of chemical penetration through the stratum corneum, viable epidermis and dermis formulated in terms of an Excel™ spreadsheet and associated add-in is presented. The model is a one-dimensional homogenization of underlying microscopic transport models for stratum corneum and dermis; viable epidermis is treated as unperfused dermis. The model's salient features are a detailed structural description of the skin layers, a combination of first-principles based transport equations and empirical partition and diffusion coefficients, and the capability of simulating a variety of exposure scenarios. Model predictions are compared with representative in vitro skin permeation data obtained from the literature using as summary parameters total absorption (Q(abs)), maximum flux (J(max)) and skin permeability coefficient (k(p)). The results of this evaluation demonstrate the current state-of-the-art in prediction of transient skin absorption and highlight areas in which further elaborations are needed to obtain satisfactory predictions.

Improved input parameters for diffusion models of skin absorption

To use a diffusion model for predicting skin absorption requires accurate estimates of input parameters on model geometry, affinity and transport characteristics. This review summarizes methods to obtain input parameters for diffusion models of skin absorption focusing on partition and diffusion coefficients. These include experimental methods, extrapolation approaches, and correlations that relate partition and diffusion coefficients to tabulated physico-chemical solute properties. Exhaustive databases on lipid-water and corneocyte protein-water partition coefficients are presented and analyzed to provide improved approximations to estimate lipid-water and corneocyte protein-water partition coefficients. The most commonly used estimates of lipid and corneocyte diffusion coefficients are also reviewed. In order to improve modeling of skin absorption in the future diffusion models should include the vertical stratum corneum heterogeneity, slow equilibration processes, the absorption from complex non-aqueous formulations, and an improved representation of dermal absorption processes. This will require input parameters for which no suitable estimates are yet available.

Detailed modeling of skin penetration--an overview

In recent years, the combination of computational modeling and experiments has become a useful tool that is proving increasingly powerful for explaining biological complexity. As computational power is increasing, scientists are able to explore ever more complex models in finer detail and to explain very complex real world data. This work provides an overview of one-, two- and three-dimensional diffusion models for penetration into mammalian skin. Besides diffusive transport this includes also binding of substances to skin proteins and metabolism. These models are based on partial differential equations that describe the spatial evolution of the transport process through the biological barrier skin. Furthermore, the work focuses on analytical and numerical techniques for this type of equations such as discretization schemes or homogenization (upscaling) techniques. Finally, the work compares different geometry models with respect to the permeability.

Mathematical and pharmacokinetic modelling of epidermal and dermal transport processes

Topical delivery to the various regions of the skin and underlying tissues, transdermal drug delivery and dermal exposure to environmental chemicals are important areas of research. Mathematical models of epidermal and dermal transport, involving penetration of a solute through various layers of the skin, metabolism in the skin and its subsequent distribution and clearance into systemic circulation from underlying tissues, play an essential role in this research area and are reviewed in this work.

SVOC exposure indoors: fresh look at dermal pathways

This paper critically examines indoor exposure to semivolatile organic compounds (SVOCs) via dermal pathways. First, it demonstrates that--in central tendency--an SVOC's abundance on indoor surfaces and in handwipes can be predicted reasonably well from gas-phase concentrations, assuming that thermodynamic equilibrium prevails. Then, equations are developed, based upon idealized mass-transport considerations, to estimate transdermal penetration of an SVOC either from its concentration in skin-surface lipids or its concentration in air. Kinetic constraints limit air-to-skin transport in the case of SVOCs that strongly sorb to skin-surface lipids. Air-to-skin transdermal uptake is estimated to be comparable to or larger than inhalation intake for many SVOCs of current or potential interest indoors, including butylated hydroxytoluene, chlordane, chlorpyrifos, diethyl phthalate, Galaxolide, geranyl acetone, nicotine (in free-base form), PCB28, PCB52, Phantolide, Texanol and Tonalide. Although air-to-skin transdermal uptake is anticipated to be slow for bisphenol A, we find that transdermal permeation may nevertheless be substantial following its transfer to skin via contact with contaminated surfaces. The paper concludes with explorations of the influence of particles and dust on dermal exposure, the role of clothing and bedding as transport vectors, and the potential significance of hair follicles as transport shunts through the epidermis.

PRACTICAL IMPLICATIONS

Human exposure to indoor pollutants can occur through dietary and nondietary ingestion, inhalation, and dermal absorption. Many factors influence the relative importance of these pathways, including physical and chemical properties of the pollutants. This paper argues that exposure to indoor semivolatile organic compounds (SVOCs) through the dermal pathway has often been underestimated. Transdermal permeation of SVOCs can be substantially greater than is commonly assumed. Transport of SVOCs from the air to and through the skin is typically not taken into account in exposure assessments. Yet, for certain SVOCs, intake through skin is estimated to be substantially larger than intake through inhalation. Exposure scientists, risk assessors, and public health officials should be mindful of the dermal pathway when estimating exposures to indoor SVOCs. Also, they should recognize that health consequences vary with exposure pathway. For example, an SVOC that enters the blood through the skin does not encounter the same detoxifying enzymes that an ingested SVOC would experience in the stomach, intestines, and liver before it enters the blood.

Kinetics and equilibrium of solute diffusion into human hair

The uptake kinetics of five molecules by hair has been measured and the effects of pH and physical chemical properties of molecules were investigated. A theoretical model is proposed to analyze the experimental data. The results indicate that the binding affinity of solute to hair, as characterized by hair-water partition coefficient, scales to the hydrophobicity of the solute and decreases dramatically as the pH increases to the dissociation constant. The effective diffusion coefficient of solute depended not only on the molecular size as most previous studies suggested, but also on the binding affinity as well as solute dissociation. It appears that the uptake of molecules by hair is due to both hydrophobic interaction and ionic charge interaction. Based on theoretical considerations of the cellular structure, composition and physical chemical properties of hair, quantitative-structure-property-relationships (QSPR) have been proposed to predict the hair-water partition coefficient (PC) and the effective diffusion coefficient (D (e)) of solute. The proposed QSPR models fit well with the experimental data. This paper could be taken as a reference for investigating the adsorption properties for polymeric materials, fibres, and biomaterials.

Determination of partition and binding properties of solutes to stratum corneum

The binding property of a number of relatively hydrophilic solutes to native and delipidized stratum corneum (SC) and their partition coefficients to extracted lipid have been measured by equilibration experiments to expand the current database which consisted of mostly hydrophobic solutes. Using the extended database, quantitative structure property relationships (QSPR) have been proposed for predicting the partition and binding coefficients of both hydrophobic and hydrophilic solutes to the SC protein, and lipid. Solute partition to the SC lipid is best fitted by PC(lip/w)=K(ow)(0.69) and solute binding to the SC protein is best described by PC(pro/w)=4.2K(ow)(0.31). The two QSPR models of solute partition to the SC lipid and binding to the SC protein have been further combined into a two-phase model to predict the overall partition coefficient of solutes to the stratum corneum (K(sc/w)). Our study not only extends the database of solute partition and binding properties of the SC to include hydrophilic solutes, but also demonstrates that the thermodynamic equilibrium properties of the SC partition and binding can be fitted with good accuracy by combining QSPR models with the multiphase and heterogeneous structures of the SC.

Toward a better understanding of pesticide dermal absorption: diffusion model analysis of parathion absorption in vitro and in vivo

Human skin absorption of radiolabeled parathion was studied in vitro at specific doses (mass loadings) of 0.4, 4.0, 41, or 117 microg/cm(2), with and without occlusion. The compound was applied in small volumes of acetone solution to split-thickness skin. Permeation of radiolabel into the receptor solutions was monitored for 76 h, after which the tissue was dissected and analyzed for residual radioactivity. For the 3 lower doses, cumulative permeation after 76 h was approximately dose-proportional, ranging from 28.5-30.5% of applied dose (unoccluded) to 45.5-55.7% (occluded). Total absorption, calculated as receptor fluid plus dermis content, followed a similar pattern. Both permeation rate and total absorption continued to increase up to the highest dose tested, consistent with results from other laboratories. These results are compared with predictions from a previously developed skin diffusion model (Kasting et al., 2008a). The model predicted total absorption to within a factor of 1.4 at 0.4 microg/cm(2) and 1.6 at 4 microg/cm(2), but substantially underpredicted absorption at the 2 higher doses. The analysis showed that parathion partitioned more favorably into the stratum corneum than the diffusion model prediction. Nevertheless, comparison of the model predictions to a previously reported human study showed that the skin absorption model, when corrected for surface losses occurring in vivo, satisfactorily described in vivo dermal absorption of parathion applied at 4 microg/cm(2) to various body sites.

A theoretical model for transdermal drug delivery from emulsions and its dependence upon formulation

This article presents a theoretical model of transdermal drug delivery from an emulsion-type vehicle that addresses the vehicle heterogeneity and incorporates the prediction of drug transport parameters as function of the vehicle composition. The basic mass transfer model considers interfacial and diffusion resistances within the emulsion and partition/diffusion phenomena across two skin compartments in series. Drug transport parameters are predicted as follows: partition coefficients are derived from regular solutions theory, drug diffusivity in the continuous phase is computed from a free volume theory with segmental motion, and permeability of the surfactant layer around droplets is estimated based on a free surface area model. These relationships are incorporated within the basic mass transfer model, so that the overall model is able to predict temporal profiles of drug release from the vehicle and of drug concentration in plasma, as a function of vehicle composition. In this way, the proposed model provides a sound physicochemical basis to support the development of new formulations and the planning of experiments. A simulated case study regarding a nitroglycerin ointment is presented in detail, illustrating how thermodynamic and kinetic factors inherent to the emulsion vehicle can modulate drug release and subsequent systemic absorption.

The Role of Non-Covalent Protein Binding in Skin Sensitisation Potency of Chemicals

Skin sensitisation is a delayed hypersensitivity reaction caused by repeated exposure to common natural and synthetic chemical allergens. It is thought that small chemical sensitisers (haptens) are required to form a strong irreversible bond with a self protein/peptide and generate an immunogenic hapten-protein complex in order to be recognised by the immune system and stimulate T cell proliferation. The sensitisers are usually electrophilic chemicals that are directly reactive with proteins or reactive intermediates (metabolites) of chemically inert compounds (prohaptens). Sensitising chemicals are also capable of weak, non-covalent association with proteins and there is an ongoing debate about the role of weak interactions of chemicals and proteins in the chemistry of allergy. The non-covalent interactions are reversible and thus have a major impact on skin/epidermal bioavailability of chemical/reactive metabolites. We investigated the relationship between the relative level of non-covalent association to a model protein and their relative potencies as determined by the EC3 values in the murine local lymph node assay (LLNA) for a number of chemicals. Using human serum albumin as a model protein, we determined that no observable relationship exists between the two parameters for the chemicals tested. Therefore, at least for this model protein, non-covalent interactions appear not to be a key determinant of allergen potency.

A spreadsheet-based method for estimating the skin disposition of volatile compounds: application to N,N-diethyl-m-toluamide (DEET)

The disposition of N,N-diethyl-3-methylbenzamide (DEET) applied to split-thickness human cadaver skin was measured in modified Franz cells maintained at 32 degrees C and fitted with a vapor trap. Ethanolic solutions of DEET (1% w/w) spiked with (14)C radiolabel were applied to skin at a dose of 10 microL per cell, corresponding to a DEET dose of 127 microg/cm(2). Room air was drawn over the skin at velocities ranging from 10-100 mL/min. Evaporation of radiolabel from the skin surface and absorption into the receptor solution were monitored for 24 hr post-dose. The percentage of radioactivity collected in the vapor trap after 24 hr increased with airflow, ranging from 16 +/- 4% at 10 mL/min to 59 +/- 7% at 70 mL/min. The percentage of radioactivity absorbed through the skin after 24 hours decreased with increasing airflow, ranging from 69 +/- 7% at 10 mL/min to 20 +/- 1% at 80 mL/min. Tissue retention after 24 hr was 6-14% of the radioactive dose with no clear correlation to airflow. This data as well as DEET absorption data from two previous in vitro studies in which dose and location (fume hood or bench top) was varied were analyzed in terms of a recently developed diffusion/evaporation model for skin implemented on an Excel spreadsheet. A priori model calculations based on independently estimated transport parameters (Model 1) were compared with calculations based on fitted parameters (Models 2 and 3). The analysis of the combined dataset (n = 272 observations) showed that the Model 1 estimates matched the cumulative disposition profiles to within a root mean square error of 12.4% of the applied dose (r(2) = 0.65), whereas the Model 2 and Model 3 fits matched to within 9.4% (r(2) = 0.80) and 6.5% (r(2) = 0.91), respectively. The Model 3 fits were obtained using a concentration-dependent diffusivity of DEET in the stratum corneum, the value of which increased 3.4-fold between low concentrations and saturation. This result was consistent with the mild skin penetration enhancement effect for DEET reported elsewhere. [Supplementary materials are available for this article. Go to the publisher's online edition of Journal of Occupational and Environmental Hygiene for the following free supplemental resource: a word document containing tables and figures including more information on the spreadsheet skin absorption model.]

Permeation of tecnazene through human skin in vitro as assessed by HS-SPME and GC-MS

Permeation of tecnazene into and through human cadaver skin in vitro was assessed using a CC-MS method employing HS-SPME for receptor solution analyses. Two doses of tecnazene dissolved in acetone, corresponding to 103 and 864 microg/cm2 of tecnazene, were applied to skin mounted on Franz diffusion cells and placed in a fume hood. Cells were either occluded with aluminum foil or left unoccluded. Total absorption of tecnazene (dermis + receptor fluid) after 48 h was 2.2-6.1% of the applied dose for the unoccluded treatments and 22-33% for the occluded treatments. Potentially absorbed dose including all tecnazene that may have eventually permeated the skin ranged from 10% unoccluded to 42-53% occluded. Accumulation in the receptor solutions was satisfactorily described by a working diffusion model after upward adjustment of the partition coefficient for tecnazene in all skin layers by a factor of 5-16 versus a priori values. However, residual amounts of tecnazene in both the epidermis and dermis were higher than those estimated from the model, suggesting the existence of tissue binding not accounted for in the calculation. The results indicate that the diffusion model as presently calibrated may significantly underestimate both systemic absorption and skin concentrations of highly lipophilic compounds, as predicted from data generated from in vitro skin permeation assays. Model predictions could be improved by better accounting for partitioning into the epidermis and dermis.

The transient dermal exposure: theory and experimental examples using skin and silicone membranes

A diffusion model is presented to account for the disposition of chemicals applied to skin as transient exposures. Two conditions are considered that apply to the skin surface following the exposure period, which are applicable to chemicals exhibiting two extremes of chemical volatility. For one case, representing highly volatile compounds, the solution is generalized to apply to multiple transient exposures. For both cases, algebraic expressions are derived to calculate the total amount of chemical that penetrates the skin. The theory is applied to experimental measurements of the in vitro penetration of diethyl phthalate applied to hairless guinea pig (HGP) skin and silicone rubber membranes (SRMs) as transient exposures. The transient exposure theory ably models the experimental data, with coefficients of determination greater than 0.97 (HGP) and greater than 0.99 (SRM). The ability of parameters derived from concurrent infinite dose experiments to predict the time course of absorption from transient exposures is explored. Discrepancies were found between measured cumulative penetration of chemical from transient exposure experiments and penetration predicted from parameters derived from infinite dose experiments, particularly for HGP. Possible reasons are explored. The current model may provide a realistic framework for estimating absorption from occupational, environmental and pharmaceutical dermal exposures.

First-principles, structure-based transdermal transport model to evaluate lipid partition and diffusion coefficients of hydrophobic permeants solely from stratum corneum permeation experiments

To account for the effect of branched, parallel transport pathways in the intercellular domain of the stratum corneum (SC) on the passive transdermal transport of hydrophobic permeants, we have developed, from first-principles, a new theoretical model-the Two-Tortuosity Model. This new model requires two tortuosity factors to account for: (1) the effective diffusion path length, and (2) the total volume of the branched, parallel transport pathways present in the SC intercellular domain, both of which may be evaluated from known values of the SC structure. After validating the Two-Tortuosity model with simulated SC diffusion experiments in FEMLAB (a finite element software package), the vehicle-bilayer partition coefficient, K(b), and the lipid bilayer diffusion coefficient, D(b), in untreated human SC were evaluated using this new model for two hydrophobic permeants, naphthol (K(b) = 225 +/- 42, D(b) = 1.7 x 10(-7) +/- 0.3 x 10(-7) cm(2)/s) and testosterone (K(b) = 92 +/- 29, D(b) = 1.9 x 10(-8) +/- 0.5 x 10(-8) cm(2)/s). The results presented in this paper demonstrate that this new method to evaluate K(b) and D(b) is comparable to, and simpler than, previous methods, in which SC permeation experiments were combined with octanol-water partition experiments, or with SC solute release experiments, to evaluate K(b) and D(b).

A multiphase microscopic diffusion model for stratum corneum permeability. II. estimation of physicochemical parameters, and application to a large permeability database

The full parameterization for the stratum corneum biphasic microtransport model presented previously in this Journal [95:620-648 (2006)] is developed through a combination of fundamental transport theory and calibration with existing data. Of the five microscopic transport properties, four (D(cor), K(cor/w), D(lip), K(lip/w)) are developed from sources independent of the existing steady-state permeability database. The fifth parameter, k(trans) (the mass transfer coefficient for transbilayer hopping), is derived from a fit of the model to the permeability data according to a modified free surface area function of the form log(10) k(trans) = A-B x (MW)(1/3). Examination of the experimental data in terms of the two dimensionless groups, R and sigma, arising from the analysis leads to the conclusion that SC permeation for most compounds is dominated by the transcellular pathway regardless of their lipophilicity, a striking departure from recent skin permeability models. Overall fit of the developed model(s) to the permeability data is somewhat better than for the Potts-Guy equation and variants thereof; however, marked improvement is seen in the estimation of lag times and the related potential for predicting skin hydration effects and transient skin permeation profiles. Simple approximations to the full numerical solution are presented that allow the developed model(s) to be implemented on a spreadsheet.

Effect of hydration on the permeation of ketoconazole through human nail plate in vitro

The impact of hydration on the permeation of the antifungal drug, ketoconazole, through excised human nails in vitro was evaluated in diffusion cell studies. Nails treated with [(3)H]ketoconazole solvent-deposited onto the dorsal surface were maintained in incubators at 32 degrees C and exposed sequentially to relative humidities (dorsal side) of 15, 40, 80 and 100% over a period of 40 days. The ventral side was bathed in a pH 7.4 phosphate buffer. Ascending and descending humidity regimens were tested. Increasing the ambient RH from 15 to 100% enhanced permeation of radiolabel associated with [(3)H]ketoconazole by a factor of three. Diffusivities estimated from these data and the associated nail water contents (estimated in a separate study) can be described by a free volume theory. Therefore, formulations or treatments, which increase nail hydration, have potential to improve topical therapy for onychomycosis, if a favorable balance between drug delivery and growth conditions for the dermatophytes can be achieved.

In-silico model of skin penetration based on experimentally determined input parameters. Part I: experimental determination of partition and diffusion coefficients

Mathematical modeling of skin transport is considered a valuable alternative of in-vitro and in-vivo investigations especially considering ethical and economical questions. Mechanistic diffusion models describe skin transport by solving Fick's 2nd law of diffusion in time and space; however models relying entirely on a consistent experimental data set are missing. For a two-dimensional model membrane consisting of a biphasic stratum corneum (SC) and a homogeneous epidermal/dermal compartment (DSL) methods are presented to determine all relevant input parameters. The data were generated for flufenamic acid (M(W) 281.24g/mol; logK(Oct/H2O) 4.8; pK(a) 3.9) and caffeine (M(W) 194.2g/mol; logK(Oct/H2O) -0.083; pK(a) 1.39) using female abdominal skin. K(lip/don) (lipid-donor partition coefficient) was determined in equilibration experiments with human SC lipids. K(cor/lip) (corneocyte-lipid) and K(DSL/lip) (DSL-lipid) were derived from easily available experimental data, i.e. K(SC/don) (SC-donor), K(lip/don) and K(SC/DSL) (SC-DSL) considering realistic volume fractions of the lipid and corneocyte phases. Lipid and DSL diffusion coefficients D(lip) and D(DSL) were calculated based on steady state flux. The corneocyte diffusion coefficient D(cor) is not accessible experimentally and needs to be estimated by simulation. Based on these results time-dependent stratum corneum concentration-depth profiles were simulated and compared to experimental profiles in an accompanying study.

Transcellular route of diffusion through stratum corneum: results from finite element models

Insight into the stratum corneum (SC) permeation pathway for hydrophilic compounds is gained by comparing experimental measurements of permeability and lag time (tlag) with the predictions of a finite element (FE) model. A database of permeability and lag time measurements (n=27) of hydrophilic compounds was compiled from the literature. Transcellular and lateral lipid diffusion pathways were modeled within a brick-and-mortar geometry representing fully hydrated human SC. Modeled tlag's for the lipid pathway are too brief to account for the experimental quantities, whereas the transcellular pathway with preferential corneocyte partitioning does account for them. Measured tlag's are highly correlated (p<0.0001) with the compound's octanol-water partition coefficient, supporting the hypothesis of an aqueous-lipid partition mechanism in the permeation of hydrophilic compounds. The importance of the lag time for identifying the diffusion pathway is demonstrated.

A multiphase microscopic diffusion model for stratum corneum permeability. I. Formulation, solution, and illustrative results for representative compounds

A two-dimensional microscopic transport model of the stratum corneum (SC) incorporating corneocytes of varying hydration and permeability embedded in an anisotropic lipid matrix is presented. Results are expressed in terms of a dimensionless permeability (P(SC/w)(comp), which is a function of two dimensionless parameters, R and sigma. R is a ratio of transbilayer to lateral molecular flows within a lipid bilayer and sigma is the ratio of (lateral) permeability in the lipid phase, D(lip)K(lip/w), to that in the corneocyte phase, D(cor)K(cor/w.) The shape of the dimensionless permeability surface is also governed by the arrangement of the SC lipids, where Model 1 represents the extreme in which lipid-phase transport can occur with no transbilayer transport, whereas Model 2 entails maximum transbilayer transport. Model calculations are exemplified by characterizing the skin permeability of four representative permeants: water, ethanol, nicotinamide, and testosterone. A comparison with experimental steady state permeability and partition data supports that the transport properties of the SC lipids are highly anisotropic, with lateral diffusivities several orders of magnitude higher than the equivalent diffusivity calculated from transbilayer hopping. Nevertheless, the calculations suggest that corneocyte-phase transport plays a major role for all four permeants. These results confirm our previous calculations on water permeability and present a marked contrast to the commonly stated doctrine that the SC transport pathway is primarily intercellular.