The Role of Corneocytes in Skin Transport Revised—A Combined Computational and Experimental Approach

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.

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.

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

Topical drug delivery: History, percutaneous absorption, and product development

Topical products, widely used to manage skin conditions, have evolved from simple potions to sophisticated delivery systems. Their development has been facilitated by advances in percutaneous absorption and product design based on an increasingly mechanistic understanding of drug-product-skin interactions, associated experiments, and a quality-by-design framework. Topical drug delivery involves drug transport from a product on the skin to a local target site and then clearance by diffusion, metabolism, and the dermal circulation to the rest of the body and deeper tissues. Insights have been provided by Quantitative Structure Permeability Relationships (QSPR), molecular dynamics simulations, and dermal Physiologically Based PharmacoKinetics (PBPK). Currently, generic product equivalents of reference-listed products dominate the topical delivery market. There is an increasing regulatory interest in understanding topical product delivery behavior under 'in use' conditions and predicting in vivo response for population variations in skin barrier function and response using in silico and in vitro findings.

Non-invasive delivery strategies for biologics

Biologics now constitute a significant element of available medical treatments. Owing to their clinical and commercial success, biologics are a rapidly growing class and have become a dominant therapeutic modality. Although most of the successful biologics to date are drugs that bear a peptidic backbone, ranging from small peptides to monoclonal antibodies (~500 residues; 150 kDa), new biologic modalities, such as nucleotide-based therapeutics and viral gene therapies, are rapidly maturing towards widespread clinical use. Given the rise of peptides and proteins in the pharmaceutical landscape, tremendous research and development interest exists in developing less-invasive or non-invasive routes for the systemic delivery of biologics, including subcutaneous, transdermal, oral, inhalation, nasal and buccal routes. This Review summarizes the current status, latest updates and future prospects for such delivery of peptides, proteins and other biologics.

Should we instruct patients to rub topical agents into skin? The evidence

Background: At least 15 factors influence the ability of compounds to penetrate the skin. Massage (rubbing) may be another factor that has gone relatively unrecognized. Method: PubMed, Google Scholar, and EMBASE databases were accessed online in March 2018 in search of studies measuring absorption through skin with and without rubbing or massage. Results: While some studies noted no difference in dermal absorption with regards to rubbing, others have demonstrated the opposite. In general, massage technique does indeed sometimes enhance dermal absorption. In addition to increase skin temperature and blood flow, rubbing likely modifies stratum corneum (SC) structure to enhance diffusion rates and increase retained penetrant amount within the skin. Conclusions: Understanding the mechanism of massage and its role in percutaneous penetration may help elucidate skin barrier function, dermal absorption, skin decontamination, and dermatotoxicology. To achieve such goals, an in vitro model that models in vivo behaviors must first be established. Subsequently, experiments with different penetrants, vehicles, massage time, and other variables may be considered.

How Predictable Are Human Stratum Corneum Lipid/Water Partition Coefficients? Assessment and Useful Correlations for Dermal Absorption

Partition coefficients between human stratum corneum lipids and water (K^sclip/w) are collected or deduced from a variety of sources in a manner that approximately doubles the available data compared to the current state-of-the-art model (Hansen et al., Adv Drug Deliv Rev. 2013;65(2):251-264). An additional datum for water itself in porcine SC that considerably extends the molecular size and lipophilicity range of the data set is considered. The data are analyzed in terms of an extended linear free energy relationship involving octanol/water partition coefficients, Abraham solvation parameters, and a secondary, power law molecular weight dependence. The optimum fit to log K^sclip/w for the full data set reduces the standard error of prediction from 0.50 for a Hansen-like model to 0.39; corresponding multiplicative errors in K^sclip/w are reduced from a factor of 3.1 to one of 2.5. The difference in performance is driven by the water datum, which requires a more complex dependence on molecular size than that afforded by Abraham parameters. In the absence of the water value, the Hansen-like model, which does not include a dependence on molecular size, is essentially optimum. A comparison is presented to fluid-phase phospholipid-water systems, which have a demonstrably different structure-property relationship.

Vehicle effects on human stratum corneum absorption and skin penetration

This study evaluated the effects of three vehicles-ethanol (EtOH), isopropyl alcohol (IPA), and isopropyl myristate (IPM)-on stratum corneum (SC) absorption and diffusion of the [¹⁴C]-model compounds benzoic acid and butenafine hydrochloride to better understand the transport pathways of chemicals passing through and resident in SC. Following application of topical formulations to human dermatomed skin for 30 min, penetration flux was observed for 24 h post dosing, using an in vitro flow-through skin diffusion system. Skin absorption and penetration was compared to the chemical-SC (intact, delipidized, or SC lipid film) binding levels. A significant vehicle effect was observed for chemical skin penetration and SC absorption. IPA resulted in the greatest levels of intact SC/SC lipid absorption, skin penetration, and total skin absorption/penetration of benzoic acid, followed by IPM and EtOH, respectively. For intact SC absorption and total skin absorption/penetration of butenafine, the vehicle that demonstrated the highest level of sorption/penetration was EtOH, followed by IPA and IPM, respectively. The percent doses of butenafine that were absorbed in SC lipid film and penetrated through skin in 24 h were greatest for IPA, followed by EtOH and IPM, respectively. The vehicle effect was consistent between intact SC absorption and total chemical skin absorption and penetration, as well as SC lipid absorption and chemical penetration through skin, suggesting intercellular transport as a main pathway of skin penetration for model chemicals. These results suggest the potential to predict vehicle effects on skin permeability with simple SC absorption assays. As decontamination was applied 30 min after chemical exposure, significant vehicle effects on chemical SC partitioning and percutaneous penetration also suggest that skin decontamination efficiency is vehicle dependent, and an effective decontamination method should act on chemical solutes in the lipid domain.

Proposed human stratum corneum water domain in chemical absorption

Compounds with varying physical and chemical properties may have different affinities to the stratum corneum (SC) and/or its intercellular lipids, keratin protein, and possible water domains. To better understand the mechanism of percutaneous absorption, we utilized 21 carbon-14 labeled chemicals, with wide hydrophilicity (log P = -0.05 to 6.17), and quantified their absorption/adsorption properties for a short incubation time (15 min) with regards to intact SC membrane, delipidized SC membrane and SC lipid. A facile method was developed for SC/lipid absorption, providing a more equivalent procedure and comparable data. SC lipid absorption of chemical solutes positively correlated with the octanol/water partition coefficient (log P). Differences between the percent dose of chemical absorption to intact SC and the total percent dose contributed by the protein and lipid domains suggest the possibility and significance of a water domain. Absorption rate experiments showed a longer lag time for intact SC than for delipidized SC or SC lipid, suggesting that the water domain may delay chemical binding to protein and lipid domains, and may be a factor in the resistance of many chemicals to current decontamination methods. Copyright © 2015 John Wiley & Sons, Ltd.

Mathematical models for dermal drug absorption

INTRODUCTION

Mathematical models of dermal transport offer the advantages of being much faster and less expensive than in vitro or in vivo studies. The number of methods used to create such models has been increasing rapidly, probably due to the steady rise in computational power. Although each of the various approaches has its own virtues and limitations, it may be difficult to decide which approach is best suited to address a given problem.

AREAS COVERED

Here we outline the basic ideas, drawbacks and advantages of compartmental and quantitative structure-activity relationship models, as well as of analytical and numerical approaches for solving the diffusion equation. Examples of special applications of the different approaches are given.

EXPERT OPINION

Although some models are sophisticated and might be used in future to predict transport through damaged or diseased skin, the comparatively low availability of suitable and accurate experimental data limits extensive usage of these models and their predictive accuracy. Due to the lack of experimental data, the possibility of validating mathematical models is limited.

A strategy for in-silico prediction of skin absorption in man

For some time, in-silico models to address substance transport into and through the skin are gaining more and more importance in different fields of science and industry. In particular, the mathematical prediction of in-vivo skin absorption is of great interest to overcome ethical and economical issues. The presented work outlines a strategy to address this problem and in particular, investigates in-vitro and in-vivo skin penetration experiments of the model compound flufenamic acid solved in an ointment by means of a mathematical model. Experimental stratum corneum concentration-depth profiles (SC-CDP) for various time intervals using two different in-vitro systems (Franz diffusion cell, Saarbruecken penetration model) were examined and simulated with the help of a highly optimized three compartment numerical diffusion model and compared to the findings of SC-CDPs of the in-vivo scenario. Fitted model input parameters (diffusion coefficient and partition coefficient with respect to the stratum corneum) for the in-vitro infinite dose case could be used to predict the in-use conditions in-vitro. Despite apparent differences in calculated partition coefficients between in-vivo and in-vitro studies, prediction of in-vivo scenarios from input parameters calculated from the in-vitro case yielded reasonable results.

Monounsaturated fatty acids reduce the barrier of stratum corneum lipid membranes by enhancing the formation of a hexagonal lateral packing

The effectiveness of the skin barrier underlies the outer layer of the skin: the stratum corneum (SC). However, in several skin diseases this barrier is impaired. In two inflammatory skin diseases, atopic eczema and Netherton syndrome, an increased level of monounsaturated fatty acids (MUFAs) has been observed as opposed to healthy skin. In the present study, we aimed to investigate the effect of MUFAs on the lipid organization and skin lipid barrier using an in vitro model membrane system, the stratum corneum substitute (SCS), mimicking the SC lipid composition and organization. To achieve our goal, the SCS has been prepared with increasing levels of MUFAs using various chain length. Permeation studies and trans-epidermal water loss measurements show that an increment of MUFAs reduces the lipid barrier in the SCS. The increased level of unsaturation exerts its effect by reducing the packing density in the lipid organization, while the lamellar phases are not affected. Our findings indicate that increased levels of MUFAs may contribute to the impaired skin barrier in diseased skin.

Active ingredients against human epidermal aging

The decisive role of the epidermis in maintaining body homeostasis prompted studies to evaluate the changes in epidermal structure and functionality over the lifetime. This development, along with the identification of molecular mechanisms of epidermal signaling, maintenance, and differentiation, points to a need for new therapeutic alternatives to treat and prevent skin aging. In addition to recovering age- and sun-compromised functions, proper treatment of the epidermis has important esthetic implications. This study reviews active ingredients capable of counteracting symptoms of epidermal aging, organized according to the regulation of specific age-affected epidermal functions: (1) several compounds, other than retinoids and derivatives, act on the proliferation and differentiation of keratinocytes, supporting the protective barrier against mechanical and chemical insults; (2) natural lipidic compounds, as well as glycerol and urea, are described as agents for maintaining water-ion balance; (3) regulation of immunological pathogen defense can be reinforced by natural extracts and compounds, such as resveratrol; and (4) antioxidant exogenous sources enriched with flavonoids and vitamin C, for example, improve solar radiation protection and epidermal antioxidant activity. The main objective is to provide a functional classification of active ingredients as regulatory elements of epidermal homeostasis, with potential cosmetic and/or dermatological applications.

A microscopic multiphase diffusion model of viable epidermis permeability

A microscopic model of passive transverse mass transport of small solutes in the viable epidermal layer of human skin is formulated on the basis of a hexagonal array of cells (i.e., keratinocytes) bounded by 4-nm-thick, anisotropic lipid bilayers and separated by 1-μm layers of extracellular fluid. Gap junctions and tight junctions with adjustable permeabilities are included to modulate the transport of solutes with low membrane permeabilities. Two keratinocyte aspect ratios are considered to represent basal and spinous cells (longer) and granular cells (more flattened). The diffusion problem is solved in a unit cell using a coordinate system conforming to the hexagonal cross section, and an efficient two-dimensional treatment is applied to describe transport in both the cell membranes and intercellular spaces, given their thinness. Results are presented in terms of an effective diffusion coefficient, D¯(epi), and partition coefficient, K¯(epi/w), for a homogenized representation of the microtransport problem. Representative calculations are carried out for three small solutes-water, L-glucose, and hydrocortisone-covering a wide range of membrane permeability. The effective transport parameters and their microscopic interpretation can be employed within the context of existing three-layer models of skin transport to provide more realistic estimates of the epidermal concentrations of topically applied solutes.

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.

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.

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.

Measuring the stratum corneum reservoir: desorption kinetics from keratin

High keratin binding and slow desorption kinetics are assumed to be responsible for the formation of the stratum corneum (SC) reservoir. We measured equilibrium binding coefficients (K(b)) and desorption rate constants (k(off)) with bovine hoof/horn keratin and six solutes with similar molecular weight (180-288 Da) and varying lipophilicities [expressed as octanol-water distribution coefficient, i.e., a partition coefficient corrected for pH (log K(pH))-0.13 to 3.8]. Two ionizable solutes within this set were tested at different pH values as degree of ionization and lipophilicity were expected to influence equilibrium binding and desorption kinetics. The unbound fraction at equilibrium varied between 18% and 93%. All solutes exhibited linear binding isotherms within the investigated concentration range. Equilibrium binding and the rate of desorption are both functions of solute lipophilicity [log K(b) = 1.23 + 0.32 log K(pH); log k(off) = 1/(25.75 + 8.35 K(pH) (0.34))]. Our results prove that slow desorption from keratin may be a major contributor to the SC reservoir. Also, they prove that reservoir formation is relevant for lipophilic solutes independent of drug class, thus allowing new options for topical pharmacotherapy.

Nanoparticles for transcutaneous vaccination

The living epidermis and dermis are rich in antigen presenting cells (APCs). Their activation can elicit a strong humoral and cellular immune response as well as mucosal immunity. Therefore, the skin is a very attractive site for vaccination, and an intradermal application of antigen may be much more effective than a subcutaneous or intramuscular injection. However, the stratum corneum (SC) is a most effective barrier against the invasion of topically applied vaccines. Products which have reached the stage of clinical testing, avoid this problem by injecting the nano‐vaccine intradermally or by employing a barrier disrupting method and applying the vaccine to a relatively large skin area. Needle‐free vaccination is desirable from a number of aspects: ease of application, improved patient acceptance and less risk of infection among them. Nanocarriers can be designed in a way that they can overcome the SC. Also incorporation into nanocarriers protects instable antigen from degradation, improves uptake and processing by APCs, and facilitates endosomal escape and nuclear delivery of DNA vaccines. In addition, sustained release systems may build a depot in the tissue gradually releasing antigen which may avoid booster doses. Therefore, nanoformulations of vaccines for transcutaneous immunization are currently a very dynamic field of research. Among the huge variety of nanocarrier systems that are investigated hopes lie on ultra‐flexible liposomes, superfine rigid nanoparticles and nanocarriers, which are taken up by hair follicles. The potential and pitfalls associated with these three classes of carriers will be discussed.

Drug delivery through the skin: molecular simulations of barrier lipids to design more effective noninvasive dermal and transdermal delivery systems for small molecules, biologics, and cosmetics

The delivery of drugs through the skin provides a convenient route of administration that is often preferable to injection because it is noninvasive and can typically be self-administered. These two factors alone result in a significant reduction of medical complications and improvement in patient compliance. Unfortunately, a significant obstacle to dermal and transdermal drug delivery alike is the resilient barrier that the epidermal layers of the skin, primarily the stratum corneum, presents for the diffusion of exogenous chemical agents. Further advancement of transdermal drug delivery requires the development of novel delivery systems that are suitable for modern, macromolecular protein and nucleotide therapeutic agents. Significant effort has already been devoted to obtain a functional understanding of the physical barrier properties imparted by the epidermis, specifically the membrane structures of the stratum corneum. However, structural observations of membrane systems are often hindered by low resolutions, making it difficult to resolve the molecular mechanisms related to interactions between lipids found within the stratum corneum. Several models describing the molecular diffusion of drug molecules through the stratum corneum have now been postulated, where chemical permeation enhancers are thought to disrupt the underlying lipid structure, resulting in enhanced permeability. Recent investigations using biphasic vesicles also suggested a possibility for novel mechanisms involving the formation of complex polymorphic lipid phases. In this review, we discuss the advantages and limitations of permeation-enhancing strategies and how computational simulations, at the atomic scale, coupled with physical observations can provide insight into the mechanisms of diffusion through the stratum corneum.

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.