Mathematical and pharmacokinetic modelling of epidermal and dermal transport processes

Skin penetration enhancers: Mechanistic understanding and their selection for formulation and design

The skin functions as a formidable barrier, particularly the stratum corneum, effectively restricting the penetration of most substances, including therapeutic agents. To circumvent this barrier, skin penetration enhancers (SPEs) are frequently employed to transiently increase skin permeability, facilitating drug absorption without causing irritation or damage. Despite advancements in dermal formulation development, a deeper understanding of the fundamental science underpinning drug delivery via SPEs remains essential. This review delivers a critical update on conventional SPEs, exploring their mechanisms in promoting drug permeation across the skin. In addition to offering an overview of percutaneous drug delivery, we examine the prevailing theories on how SPEs enhance drug transport. Furthermore, we address the intricate interplay between SPEs, drugs and the skin, providing valuable insights into how the molecular properties and permeation behaviours of SPEs influence their efficacy. This comprehensive review aims to support the ongoing development of optimised drug delivery systems for dermal applications by elucidating the complexities and challenges involved in using SPEs effectively.

A computational workflow for end-to-end simulation of percutaneous absorption

PURPOSE

Presented is the development of a workflow for end-to-end (e2e) in silico modelling of percutaneous absorption under a range of test conditions, integrating multiple calculation and analysis steps for in-silico simulation of dermal absorption. The aim is to achieve a digital twin that can be used by non-modelling experts to simulate transdermal permeation.

METHODS

A KNIME-based toolbox is used to create the workflow for the E2E in-silico model. The workflow first combines physicochemical property informatics (ChemAxon), molecular dynamics (MD) modelling, and quantitative structure-property relations (QSPRs) to calculate the diffusion and partition properties of permeants in heterogeneous skin layers and complex formulation vehicles. These are then set as input parameters to physiologically based pharmacokinetics (PBPK) model to simulate percutaneous absorption under complex in vitro testing or in vivo exposure conditions set by the end user. Integrated into the PBPK model is the evaporation of volatile permeants and solvents for in vitro unoccluded conditions. The workflow generates a report of the results and records the tested formulation in a database.

RESULTS

The workflow has been tested against several sets of published in vitro permeation test (IVPT) results of percutaneous absorption involving different formulation vehicles. The model predictions of formulation and evaporation effects on percutaneous absorption agreed well with experimental data.

CONCLUSIONS

By automating multiple calculation steps from permeant property, diffusion-partition skin layers and formulation vehicles, to PBPK modelling of dermal absorption, the workflow provides a user-friendly means for non-modelling experts to conduct in-silico simulations of transdermal absorption under various conditions. The workflow is robust to simulate the impact of complex formulation and exposure conditions including evaporation of volatile permeants and solvents on the delivery into the skin.

Mathematical modelling of hollow microneedle-mediated transdermal drug delivery

Hollow microneedles represent a promising approach for overcoming the protective barrier of the stratum corneum, facilitating direct drug infusion into viable skin tissue and thereby enhancing the efficacy of transdermal delivery. However, delivery outcomes across different skin layers and into the systemic circulation can vary substantially due to the diverse properties of drug delivery systems, clinical settings, and environmental factors. The optimal strategies for enhancing the efficiency of hollow microneedle-mediated transdermal drug delivery remain to be elucidated. This study employs mathematical modelling and a reconstructed skin model with realistic anatomical structures to investigate drug transport and accumulation across different skin layers and into the bloodstream under different delivery conditions. The modelling results reveal the crucial role of interstitial fluid flow in determining drug transport in this transdermal delivery. Delivery outcomes of each skin layer and blood exhibit distinct responses to changes in delivery conditions. Specifically, increasing the vascular permeability or nanocarrier diffusivity raises drug concentration in the blood or reticular dermis, respectively, while leading to reductions in other skin layers. The use of microneedles with narrower infusion channels can only enhance drug availability in the viable epidermis. Optimisation requires a tailored approach to several parameters depending on the target skin layer, including drug release rate, infusion rate, infusion duration, and microneedle length. Environmental factors that promote trans-epidermal water loss can increase drug concentration in the viable epidermis but have a limited impact on deeper skin tissues. The findings support the selection or customisation of hollow microneedles and nanocarriers to address specific therapeutic needs, such as targeting specific skin layers or systemic circulation, while minimising the risk of side effects from high drug concentrations in normal tissues. This study provides guidance for optimising transdermal drug delivery systems.

Modelling Hollow Microneedle-Mediated Drug Delivery in Skin Considering Drug Binding

Background/Objectives: Microneedle(MN)-based drug delivery is one of the potential approaches to overcome the limitations of oral and hypodermic needle delivery. An in silico model has been developed for hollow microneedle (HMN)-based drug delivery in the skin and its subsequent absorption in the blood and tissue compartments in the presence of interstitial flow. The drug's reversible specific saturable binding to its receptors and the kinetics of reversible absorption across the blood and tissue compartments have been taken into account. Methods: The governing equations representing the flow of interstitial fluid, the transport of verapamil in the viable skin and the concentrations in the blood and tissue compartments are solved using combined Marker and Cell and Immersed Boundary Methods to gain a quantitative understanding of the model under consideration. Results: The viscoelastic skin is predicted to impede the transport of verapamil in the viable skin and, hence, reduce the concentrations of all forms in the blood and the tissue compartments. The findings reveal that a higher mean concentration in the viable skin is not always associated with a longer MN length. Simulations also predict that the concentrations of verapamil in the blood and bound verapamil in the tissue compartment rise with decreasing tip diameters. In contrast, the concentration of free verapamil in the tissue increases with increasing injection velocities. Conclusions: The novelty of this study includes verapamil metabolism in two-dimensional viscoelastic irregular viable skin and the nonlinear, specific, saturable, and reversible binding of verapamil in the tissue compartment. The tip diameter and the drug's injection velocity are thought to serve as regulatory parameters for the effectiveness and efficacy of MN-mediated therapy if the MN is robust enough to sustain the force needed to penetrate a wider tip into the skin.

Inertial Memory Effects in Molecular Transport Across Nanoporous Membranes

Nanoporous membranes are heterogeneous structures, with heterogeneity manifesting at the microscale. In examining particle transport through such media, it has been observed that this transport deviates from classical diffusion, as described by Fick's second law. Moreover, the classical model is physically unsustainable, as it is non-causal and predicts an infinite speed of concentration perturbation propagation through a substantial medium. In this work, we have derived two causal models as extensions of Fick's second law, where causality is linked to the effects of inertial memory in the nanoporous membrane. The results of the derived models have been compared with each other and with those obtained from the classical model. It has been demonstrated that both causal models, one with exponentially fading inertial memory and the other with power-law fading memory, predict that the concentration perturbation propagates as a damped wave, leading to an increased time required for the cumulative amount of molecules passing through the membrane to reach a steady state compared to the classical model. The power-law fading memory model predicts a longer time required to achieve a stationary state. These findings have significant implications for understanding cell physiology, developing drug delivery systems, and designing nanoporous membranes for various applications.

An advanced in silico model of the oral mucosa reveals the impact of extracellular spaces on chemical permeation

Accurately predicting the permeation of chemicals through human epithelial tissues is crucial for pharmaceutical therapeutic design and toxicology. Current mathematical models of multi-layered stratified squamous epithelium such as those in the oral cavity use simplistic 'bricks and mortar' geometries that do not fully account for the complex cellular architecture that may affect chemical permeation in these tissues. Here we aimed to develop a new, advanced mechanistic mathematical model of the human epithelium that more accurately represents chemical tissue permeation. Using measurements of cell size and tortuosity from micrograph images of both human oral (buccal) and tissue-engineered buccal mucosa along with mechanistic mathematical modelling, we show that the convoluted geometry of the extracellular spaces within the epithelium significantly impacts chemical permeation. We next developed an advanced histologically and physiologically-relevant in silico model of buccal mucosal chemical permeation using partial differential equations, fitted to chemical permeation from in vitro assay data derived from tissue-engineered buccal mucosal models and chemicals with known physiochemical properties. Our novel in silico model can predict epithelial permeation kinetics for chemicals with different physicochemical properties in the absence or presence of permeability enhancers. This in vitro - in silico approach constitutes a step-change in the modelling of chemical tissue permeation and has the potential to expedite pharmaceutical innovation by improved and more rapid screening of chemical entities whilst reducing the need for in vivo animal experiments.

Is personal physiology-based rapid prediction digital twin for minimal effective fentanyl dose better than standard practice: a pilot study protocol

INTRODUCTION

Patients with advanced cancer frequently suffer from chronic, severe disabling pain. Opioids such as morphine and fentanyl are commonly used to manage this pain. Transdermal drug delivery systems are important technologies for administering drugs in a non-invasive, continuous and controlled manner. Due to the narrow therapeutic range of fentanyl, individualised dosing is essential to avoid underdosing or overdosing. Standard clinical calculation tools for opioid rotation however do not include important patient characteristics that account for interindividual variability of opioid pharmacology.

METHODS AND ANALYSIS

We developed a clinical protocol to optimise individual fentanyl dosing in patients with advanced cancer switching from oral or intravenous opioids to transdermal fentanyl by using a physics-based digital twin (DT) that is fed by important clinical and physiological parameters. Individual tailoring of transdermal fentanyl therapy is an approach with the potential for personalised and effective care with an improved benefit-risk ratio. However, clinical validation of physics-based digital twins (PBDT) dosing is crucial to proving clinical benefit.Therapeutic drug monitoring will allow to validate the accuracy of PBDT predictions. Additional monitoring for breathing dynamics, sequential pain levels and fentanyl-related adverse events will contribute to evaluating the performance of PBDT-based dosing of transdermal fentanyl. The primary objective of the study is to develop an experimental protocol to validate DT-guided fentanyl dosing in patients with advanced cancer. This clinical study will bring individualised opioid dosing closer to clinical practice.

ETHICS AND DISSEMINATION

Study documents have been approved by the responsible Ethics Committee and study initiation is planned for late summer 2024. Data will be shared with the scientific community no more than 1 year following completion of the study and data assembly.

Exploring the thermally-controlled fentanyl transdermal therapy to provide constant drug delivery by physics-based digital twins

Transdermal drug delivery is suitable for low-molecular-weight drugs with specific lipophilicity, like fentanyl, which is widely used for cancer-induced pain management. However, fentanyl's transdermal therapy displays high intra-individual variability. Factors like skin characteristics at application sites and ambient temperature contribute to this variation. In this study, we developed a physics-based digital twin of the human body to cope with this variability and propose better adapted setups. This twin includes an in-silico skin model for drug penetration, a pharmacokinetic model, and a pharmacodynamic model. Based on the results of our simulations, applying the patch on the flank (side abdominal area) showed a 15.3 % higher maximum fentanyl concentration in the plasma than on the chest. Additionally, the time to reach this maximum concentration when delivered through the flank was 19.8 h, which was 10.3 h earlier than via the upper arm. Finally, this variation led to an 18 % lower minimum pain intensity for delivery via the flank than the chest. Moreover, the impact of seasonal changes on ambient temperature and skin temperature by considering the activity level was investigated. Based on our result, the fentanyl uptake flux by capillaries increased by up to 11.8 % from an inactive state in winter to an active state in summer. We also evaluated the effect of controlling fentanyl delivery by adjusting the temperature of the patch to alleviate the pain to reach a mild pain intensity (rated three on the VAS scale). By implementing this strategy, the average pain intensity decreased by 1.1 points, and the standard deviation for fentanyl concentration in plasma and average pain intensity reduced by 37.5 % and 33.3 %, respectively. Therefore, our digital twin demonstrated the efficacy of controlled drug release through temperature regulation, ensuring the therapy toward the intended target outcome and reducing therapy outcome variability. This holds promise as a potentially useful tool for physicians.

The Influence of the Variable Wettability Characteristics of Layers on the Transport of Nanoparticles in the Context of Drug Delivery in Skin Structures

The rapid development of nanotechnology has offered the possibility of creating nanosystems that can be used as drug carriers. The use of such carriers offers real opportunities for the development of non-invasive drug delivery through skin structures. However, in addition to the ability to create suitable nanocarriers, it is also necessary to know how they move through dermal layers. The human skin consists of layers with different wettability characteristics, which greatly complicates how introduced substances move through it. In this work, an experimental study of the diffusion process of nanoparticles through partitions with different wettability properties was carried out. Conventional diffusion tests using Franz chambers were used for this purpose. We quantified how the wettability of the barrier, the number of layers, and their mutual configuration affect the transport of nanoparticles. Based on the results, an analysis of the phenomena taking place, depending on the wettability of the partition, was carried out. A model relationship was also proposed to determine the effective diffusion coefficient, taking into account the influence of the wettability and porosity of the barrier.

Mathematical modeling of transdermal delivery of topical drug formulations in a dynamic microfluidic diffusion chamber in health and disease

Mathematical models of epidermal and dermal transport are essential for optimization and development of products for percutaneous delivery both for local and systemic indication and for evaluation of dermal exposure to chemicals for assessing their toxicity. These models often help directly by providing information on the rate of drug penetration through the skin and thus on the dermal or systemic concentration of drugs which is the base of their pharmacological effect. The simulations are also helpful in analyzing experimental data, reducing the number of experiments and translating the in vitro investigations to an in-vivo setting. In this study skin penetration of topically administered caffeine cream was investigated in a skin-on-a-chip microfluidic diffusion chamber at room temperature and at 32°C. Also the transdermal penetration of caffeine in healthy and diseased conditions was compared in mouse skins from intact, psoriatic and allergic animals. In the last experimental setup dexamethasone, indomethacin, piroxicam and diclofenac were examined as a cream formulation for absorption across the dermal barrier. All the measured data were used for making mathematical simulation in a three-compartmental model. The calculated and measured results showed a good match, which findings indicate that our mathematical model might be applied for prediction of drug delivery through the skin under different circumstances and for various drugs in the novel, miniaturized diffusion chamber.

A compartment model to predict in vitro finite dose absorption of chemicals by human skin

Dermal uptake is an important and complex exposure route for a wide range of chemicals. Dermal exposure can occur due to occupational settings, pharmaceutical applications, environmental contamination, or consumer product use. The large range of both chemicals and scenarios of interest makes it difficult to perform generalizable experiments, creating a need for a generic model to simulate various scenarios. In this study, a model consisting of a series of four well-mixed compartments, representing the source solution (vehicle), stratum corneum, viable tissue, and receptor fluid, was developed for predicting dermal absorption. The model considers experimental conditions including small applied doses as well as evaporation of the vehicle and chemical. To evaluate the model assumptions, we compare model predictions for a set of 26 chemicals to finite dose in-vitro experiments from a single laboratory using steady-state permeability coefficient and equilibrium partition coefficient data derived from in-vitro experiments of infinite dose exposures to these same chemicals from a different laboratory. We find that the model accurately predicts, to within an order of magnitude, total absorption after 24 h for 19 of these chemicals. In combination with key information on experimental conditions, the model is generalizable and can advance efficient assessment of dermal exposure for chemical risk assessment.

Model-Informed Drug Development: In Silico Assessment of Drug Bioperformance following Oral and Percutaneous Administration

The pharmaceutical industry has faced significant changes in recent years, primarily influenced by regulatory standards, market competition, and the need to accelerate drug development. Model-informed drug development (MIDD) leverages quantitative computational models to facilitate decision-making processes. This approach sheds light on the complex interplay between the influence of a drug's performance and the resulting clinical outcomes. This comprehensive review aims to explain the mechanisms that control the dissolution and/or release of drugs and their subsequent permeation through biological membranes. Furthermore, the importance of simulating these processes through a variety of in silico models is emphasized. Advanced compartmental absorption models provide an analytical framework to understand the kinetics of transit, dissolution, and absorption associated with orally administered drugs. In contrast, for topical and transdermal drug delivery systems, the prediction of drug permeation is predominantly based on quantitative structure-permeation relationships and molecular dynamics simulations. This review describes a variety of modeling strategies, ranging from mechanistic to empirical equations, and highlights the growing importance of state-of-the-art tools such as artificial intelligence, as well as advanced imaging and spectroscopic techniques.

Strat-M® positioning for skin permeation studies: A comparative study including EpiSkin® RHE, and human skin

In the development and optimization of dermatological products, In Vitro Permeation Testing (IVPT) is pivotal for controlled study of skin penetration. To enhance standardization and replicate human skin properties reconstructed human skin and synthetic membranes are explored as alternatives. Strat-M® is a membrane designed to mimic the multi-layered structure of human skin for IVPT. For instance, in Strat-M®, the steady-state fluxes (JSS) of resorcinol in formulations free of permeation enhancers were found to be 41 ± 5 µg/cm2·h for the aqueous solution, 42 ± 6 µg/cm2·h for the hydrogel, and 40 ± 6 µg/cm2·h for the oil-in-water emulsion. These results were closer to excised human skin (5 ± 3, 9 ± 2, 13 ± 6 µg/cm2·h) and surpassed the performance of EpiSkin® RHE (138 ± 5, 142 ± 6, and 162 ± 11 µg/cm2·h). While mass spectrometry and Raman microscopy demonstrated the qualitative molecular similarity of EpiSkin® RHE to human skin, it was the porous and hydrophobic polymer nature of Strat-M® that more faithfully reproduced the skin's diffusion-limiting barrier. Further validation through similarity factor analysis (∼80-85%) underscored Strat-M®'s significance as a reliable substitute for human skin, offering a promising approach to enhance realism and reproducibility in dermatological product development.

Poroelastic behavior and water permeability of human skin at the nanoscale

Topical skin care products and hydrating compositions (moisturizers or injectable fillers) have been used for years to improve the appearance of, for example facial wrinkles, or to increase "plumpness". Most of the studies have addressed these changes based on the overall mechanical changes associated with an increase in hydration state. However, little is known about the water mobility contribution to these changes as well as the consequences to the specific skin layers. This is important as the biophysical properties and the biochemical composition of normal stratum corneum, epithelium, and dermis vary tremendously from one another. Our current studies and results reported here have focused on a novel approach (dynamic atomic force microscopy-based nanoindentation) to quantify biophysical characteristics of individual layers of ex vivo human skin. We have discovered that our new methods are highly sensitive to the mechanical properties of individual skin layers, as well as their hydration properties. Furthermore, our methods can assess the ability of these individual layers to respond to both compressive and shear deformations. In addition, since human skin is mechanically loaded over a wide range of deformation rates (frequencies), we studied the biophysical properties of skin over a wide frequency range. The poroelasticity model used helps to quantify the hydraulic permeability of the skin layers, providing an innovative method to evaluate and interpret the impact of hydrating compositions on water mobility of these different skin layers.

Mathematical modelling of microneedle-mediated transdermal delivery of drug nanocarriers into skin tissue and circulatory system

Microneedle-mediated transdermal delivery using nanocarriers can successfully overcome the barrier of the stratum corneum and protect drugs from elimination in skin tissues. However, the effectiveness of drug delivery to different layers of skin tissues and the circulatory system varies considerably, subject to the properties of the drug delivery system and delivery regime. How to maximise delivery outcomes remains unclear. In this study, mathematical modelling is employed to investigate this transdermal delivery under various conditions, using the skin model that is reconstructed based on the realistic skin anatomical structure. Treatment efficacy is evaluated in terms of drug exposure over time. The modelling results demonstrate the complex dependence of drug accumulation and distribution on the nanocarrier properties, microneedle properties and environment in different skin layers and blood. Specifically, delivery outcomes in the entire skin and blood can be improved by increasing the loading dose and reducing microneedle spacing. However, several parameters need to be optimised with respect to the specific location of the target site in the tissue for better treatment; these include the drug release rate, nanocarrier diffusivity in microneedle and skin tissue, nanocarrier transvascular permeability, nanocarrier partition coefficient between tissue and microneedle, microneedle length, wind speed and relative humidity. The delivery is less sensitive to the diffusivity and physical degradation rate of free drugs in microneedle, and their partition coefficient between tissue and microneedle. Results obtained from this study can be used to improve the design of the microneedle-nanocarrier combined drug delivery system and delivery regime.

Novel Insights into the Dermal Bioaccessibility and Human Exposure to Brominated Flame Retardant Additives in Microplastics

In this study, we optimized and applied an in vitro physiologically based extraction test to investigate the dermal bioaccessibility of polybrominated diphenyl ethers (PBDEs) and hexabromocyclododecane (HBCDD), incorporated as additives in different types of microplastics (MPs), and assess human dermal exposure to these chemicals. The dermal bioaccessibility of PBDEs in polyethylene (PE) MPs was significantly higher (P < 0.05) than in polypropylene (PP) MPs. Both log K_ow and water solubility influenced the dermal bioaccessibility of PBDEs. For HBCDDs in polystyrene MPs, the dermally bioaccessible fractions were 1.8, 2.0, and 1.6% of the applied dose for α-, β-, and γ-HBCDDs, respectively. MP particle size and the presence of cosmetic formulations (antiperspirant, foundation, moisturizer and sunscreen) influenced the bioaccessibility of PBDEs and HBCDDs in MP matrices at varying degrees of significance. Human exposure to ∑PBDEs and ∑HBCDDs via dermal contact with MPs ranged from 0.02 to 22.2 and 0.01 to 231 ng (kg bw)^-1 d^-1 and from 0.02 to 6.27 and 0.2 to 65 ng (kg bw)^-1 d^-1 for adults and toddlers, respectively. Dermal exposure to PBDEs and HBCDDs in MPs is substantial, highlighting for the first time the significance of the dermal pathway as a major route of human exposure to additive chemicals in microplastics.

Dynamic Modeling of Carbon Dioxide Transport through the Skin Using a Capnometry Wristband

The development of a capnometry wristband is of great interest for monitoring patients at home. We consider a new architecture in which a non-dispersive infrared (NDIR) optical measurement is located close to the skin surface and is combined with an open chamber principle with a continuous circulation of air flow in the collection cell. We propose a model for the temporal dynamics of the carbon dioxide exchange between the blood and the gas channel inside the device. The transport of carbon dioxide is modeled by convection-diffusion equations. We consider four compartments: blood, skin, the measurement cell and the collection cell. We introduce the state-space equations and the associated transition matrix associated with a Markovian model. We define an augmented system by combining a first-order autoregressive model describing the supply of carbon dioxide concentration in the blood compartment and its inertial resistance to change. We propose to use a Kalman filter to estimate the carbon dioxide concentration in the blood vessels recursively over time and thus monitor arterial carbon dioxide blood pressure in real time. Four performance factors with respect to the dynamic quantification of the CO₂ blood concentration are considered, and a simulation is carried out based on data from a previous clinical study. These demonstrate the feasibility of such a technological concept.

Migration of polycyclic aromatic hydrocarbons from a polymer surrogate through the stratum corneum layer of the skin

In this study, we determined partition (K_sc/m) and diffusion (D_sc) coefficients of five different polycyclic aromatic hydrocarbons (PAH) migrating from squalane into and through the stratum corneum (s.c.) layer of the skin. Carcinogenic PAH have previously been detected in numerous polymer-based consumer products, especially those dyed with carbon black. Upon dermal contact with these products, PAH may penetrate into and through the viable layers of the skin by passing the s.c. and thus may become bioavailable. Squalane, a frequent ingredient in cosmetics, has also been used as a polymer surrogate matrix in previous studies. K_sc/m and D_sc are relevant parameters for risk assessment because they allow estimating the potential of a substance to become bioavailable upon dermal exposure. We developed an analytical method involving incubation of pigskin with naphthalene, anthracene, pyrene, benzo[a]pyrene and dibenzo[a,h]pyrene in Franz diffusion cell assays under quasi-infinite dose conditions. PAH were subsequently quantified within individual s.c. layers by gas chromatography coupled to tandem mass spectrometry. The resulting PAH depth profiles in the s.c. were fitted to a solution of Fick's second law of diffusion, yielding K_sc/m and D_sc. The decadic logarithm logK_sc/m ranged from -0.43 to +0.69 and showed a trend to higher values for PAH with higher molecular masses. D_sc, on the other hand, was similar for the four higher molecular mass PAH but about 4.6-fold lower than for naphthalene. Moreover, our data suggests that the s.c./viable epidermis boundary layer represents the most relevant barrier for the skin penetration of higher molecular mass PAH. Finally, we empirically derived a mathematical description of the concentration depth profiles that better fits our data. We correlated the resulting parameters to substance specific constants such as the logarithmic octanol-water partition coefficient logP, K_sc/m and the removal rate at the s.c./viable epidermis boundary layer.

In-Silico Modelling of Transdermal Delivery of Macromolecule Drugs Assisted by a Skin Stretching Hypobaric Device

OBJECTIVES

To develop a simulation model to explore the interplay between mechanical stretch and diffusion of large molecules into the skin under locally applied hypobaric pressure, a novel penetration enhancement method.

METHODS

Finite element method was used to model the skin mechanical deformation and molecular diffusion processes, with validation against in-vitro transdermal permeation experiments. Simulations and experimental data were used together to investigate the transdermal permeation of large molecules under local hypobaric pressure.

RESULTS

Mechanical simulations resulted in skin stretching and thinning (20%-26% hair follicle diameter increase, and 21%-27% skin thickness reduction). Concentration of dextrans in the stratum corneum was below detection limit with and without hypobaric pressure. Concentrations in viable epidermis and dermis were not affected by hypobaric pressure (approximately 2 μg [Formula: see text] cm-2). Permeation into the receptor fluid was substantially enhanced from below the detection limit at atmospheric pressure to up to 6 μg [Formula: see text] cm-2 under hypobaric pressure. The in-silico simulations compared satisfactorily with the experimental results at atmospheric conditions. Under hypobaric pressure, satisfactory comparison was attained when the diffusion coefficients of dextrans in the skin layers were increased from [Formula: see text] 10 μm2 [Formula: see text] s-1 to between 200-500 μm2 [Formula: see text] s-1.

CONCLUSIONS

Application of hypobaric pressure induces skin mechanical stretching and enlarges the hair follicle. This enlargement alone cannot satisfactorily explain the increased transdermal permeation into the receptor fluid under hypobaric pressure. The results from the in-silico simulations suggest that the application of hypobaric pressure increases diffusion in the skin, which leads to improved overall transdermal permeation.

Multi-phase multi-layer mechanistic dermal absorption (MPML MechDermA) model to predict local and systemic exposure of drug products applied on skin

Physiologically-based pharmacokinetic models combine knowledge about physiology, drug product properties, such as physicochemical parameters, absorption, distribution, metabolism, excretion characteristics, formulation attributes, and trial design or dosing regimen to mechanistically simulate drug pharmacokinetics (PK). The current work describes the development of a multiphase, multilayer mechanistic dermal absorption (MPML MechDermA) model within the Simcyp Simulator capable of simulating uptake and permeation of drugs through human skin following application of drug products to the skin. The model was designed to account for formulation characteristics as well as body site- and sex- population variability to predict local and systemic bioavailability. The present report outlines the structure and assumptions of the MPML MechDermA model and includes results from simulations comparing absorption at multiple body sites for two compounds, caffeine and benzoic acid, formulated as solutions. Finally, a model of the Feldene (piroxicam) topical gel, 0.5% was developed and assessed for its ability to predict both plasma and local skin concentrations when compared to in vivo PK data.

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.

The Combined Diffusion and Adsorption Concept for Prediction of Nanoparticles Transport through Dermal Layers Based on Experiments in Membranes

The non-invasive introduction of active substances into the human body is a top challenge for researchers in medicine, pharmacology, and cosmetology. Development of nanotechnology and possibilities of creating more and more complex drug carriers on a nanoscale give a more realistic prospect of meeting this challenge. However, in the absence of sufficient knowledge of the mechanisms of such systems' transport through the human skin structure, it is necessary to look deeper into these issues. There are several models describing nanoparticles transport through the skin, but they are mainly based on diffusion process analysis. In this work, a model was proposed to predict nanoparticles transport through the skin, based on the combined diffusion and adsorption concept. This approach was based on experimental studies of silver and copper nanoparticles' diffusion process through different filtration membrane layers. Dependence of the degree of adsorption on the surface parameter was described using modified Langmuir equation. Then, these considerations were related to the structure of the stratum corneum, which made it possible to predict the changes in the mass of penetrating nanoparticles as a function of transport path length. A discussion of the presented model, depending on such parameters as nanoparticle size, skin cell thickness, or viscosity of the "intercellular cement", was also performed.

Microemulsions and Nanoemulsions in Skin Drug Delivery

Microemulsions and nanoemulsions are lipid-based pharmaceutical systems with a high potential to increase the permeation of drugs through the skin. Although being isotropic dispersions of two nonmiscible liquids (oil and water), significant differences are encountered between microemulsions and nanoemulsions. Microemulsions are thermodynamically stable o/w emulsions of mean droplet size approximately 100–400 nm, whereas nanoemulsions are thermodynamically unstable o/w emulsions of mean droplet size approximately 1 to 100 nm. Their inner oil phase allows the solubilization of lipophilic drugs, achieving high encapsulation rates, which are instrumental for drug delivery. In this review, the importance of these systems, the key differences regarding their composition and production processes are discussed. While most of the micro/nanoemulsions on the market are held by the cosmetic industry to enhance the activity of drugs used in skincare products, the development of novel pharmaceutical formulations designed for the topical, dermal and transdermal administration of therapeutic drugs is being considered. The delivery of poorly water-soluble molecules through the skin has shown some advantages over the oral route, since drugs escape from first-pass metabolism; particularly for the treatment of cutaneous diseases, topical delivery should be the preferential route in order to reduce the number of drugs used and potential side-effects, while directing the drugs to the site of action. Thus, nanoemulsions and microemulsions represent versatile options for the delivery of drugs through lipophilic barriers, and many synthetic and natural compounds have been formulated using these delivery systems, aiming to improve stability, delivery and bioactivity. Detailed information is provided concerning the most relevant recent scientific publications reporting the potential of these delivery systems to increase the skin permeability of drugs with anti-inflammatory, sun-protection, anticarcinogenic and/or wound-healing activities. The main marketed skincare products using emulsion-based systems are also presented and discussed.

Predicting transdermal fentanyl delivery using physics-based simulations for tailored therapy based on the age

Transdermal fentanyl patches are an effective alternative to the sustained release of oral morphine for chronic pain management. Due to the narrow therapeutic range of fentanyl, the concentration of fentanyl in the blood needs to be carefully monitored. Only then can effective pain relief be achieved while avoiding adverse effects such as respiratory depression. This study developed a physics-based digital twin of a patient by implementing drug uptake, pharmacokinetics, and pharmacodynamics models. The twin was employed to predict the in-silico effect of conventional fentanyl transdermal in a 20–80-year-old virtual patient. The results show that, with increasing age, the maximum transdermal fentanyl flux and maximum concentration of fentanyl in the blood decreased by 11.4% and 7.0%, respectively. However, the results also show that as the patient's age increases, the pain relief increases by 45.2%. Furthermore, the digital twin was used to propose a tailored therapy based on the patient's age. This predesigned therapy customized the duration of applying the commercialized fentanyl patches. According to this therapy, a 20-year-old patient needs to change the patch 2.1 times more frequently than conventional therapy, which leads to 30% more pain relief and 315% more time without pain. In addition, the digital twin was updated by the patient's pain intensity feedback. Such therapy increased the patient's breathing rate while providing effective pain relief, so a safer treatment. We quantified the added value of a patient's physics-based digital twin and sketched the future roadmap for implementing such twin-assisted treatment into the clinics.

Development of a diffusion-weighed mathematical model for intradermal drainage quantification

The quantitative assessment of lymphatic dermal clearance using NIR fluorescent tracers is particularly important for the early diagnosis of several potential disabling diseases. Currently, half-life values are computed using a mono-exponential mathematical model, neglecting diffusion of the tracer within the dermis after injection. The size and position of the region of interest are subjectively manually selected around the point of injection on the skin surface where the fluorescence signal intensity is averaged, neglecting any spatial information contained in the image. In this study we present and test a novel mathematical model allowing the objective quantification of dermal clearance, taking into consideration potential dermal diffusion. With only two parameters, this "clearance-diffusion" model is simple enough to be applied in a variety of settings and requires almost no prior information about the system. We demonstrate that if dermal diffusion is low, the mono-exponential approach is suitable but still lacking objectivity. However, if dermal diffusion is substantial, the clearance-diffusion model is superior and allows the accurate calculation of half-life values.

Development and Evaluation of an In Silico Dermal Absorption Model Relevant for Children

The higher skin surface area to body weight ratio in children and the prematurity of skin in neonates may lead to higher chemical exposure as compared to adults. The objectives of this study were: (i) to provide a comprehensive review of the age-dependent anatomical and physiological changes in pediatric skin, and (ii) to construct and evaluate an age-dependent pediatric dermal absorption model. A comprehensive review was conducted to gather data quantifying the differences in the anatomy and physiology of child and adult skin. Maturation functions were developed for model parameters that were found to be age-dependent. A pediatric dermal absorption model was constructed by updating a MoBi implementation of the Dancik et al. 2013 skin permeation model with these maturation functions. Using a workflow for adult-to-child model extrapolation, the predictive performance of the model was evaluated by comparing its predicted rates of flux of diamorphine, phenobarbital and buprenorphine against experimental observations using neonatal skin. For diamorphine and phenobarbital, the model provided reasonable predictions. The ratios of predicted:observed flux in neonates for diamorphine ranged from 0.55 to 1.40. For phenobarbital, the ratios ranged from 0.93 to 1.26. For buprenorphine, the model showed acceptable predictive performance. Overall, the physiologically based pediatric dermal absorption model demonstrated satisfactory prediction accuracy. The prediction of dermal absorption in neonates using a model-based approach will be useful for both drug development and human health risk assessment.

Mathematical modelling of nanoparticle-mediated topical drug delivery to skin tissue

Nanoparticles have been extensively studied to improve drug delivery outcomes, however, their use in topical delivery remains controversial. Although the feasibility to cross the human skin barrier has been demonstrated in experiments, the risk of low drug concentration in deep tissue still limits the application. In this study, mathematical modelling is employed to examine the performance of nanoparticle-mediated topical delivery for sending drugs into the deep skin tissue. The pharmacokinetic effect is evaluated based on the drug exposure over time. As compared to the delivery using plain drugs, nanoparticle-mediated topical delivery has the potential to significantly improve the drug exposure in deep skin tissue. Modelling predictions denote that the importance of sufficient long-term drug-skin contact in achieving effective drug deposition in the deep skin tissue. The delivery outcomes are highly sensitive to the release rate. Accelerating the release from nanoparticles in stratum corneum is able to improve the drug exposure in stratum corneum and viable epidermis while resulting in the reductions in dermis and blood. The release rate in stratum corneum and viable epidermis should be well-designed below a threshold for generating effective drug accumulation in dermis and blood. A more localised drug accumulation can be achieved in the capillary-rich region of dermis by increasing the local release rate. The release rate in dermis needs to be optimised to increase the drug exposure in the dermis region where there are fewer blood and lymphatics capillaries. Results from this study can be used to improve the regimen of topical delivery for localised treatment.

High-Payload Nanosuspension of Centella asiatica Extract for Improved Skin Delivery with No Irritation

Background

The titrated extract of Centella asiatica (CA) has received much attention as a cosmeceutical ingredient owing to its anti-wrinkle effect. However, due to the low solubility and high molecular weight of pharmacologically active constituents, including asiatic acid (AA), madecassic acid (MA), and asiaticoside (AS), it is challenging to fabricate high-payload topical preparations of CA with satisfactory skin absorption profiles.

Purpose

This study aimed to design a high-payload topical preparation of CA using nanocrystallization technique and to evaluate its skin absorption profile and local tolerability.

Methods

High-payload nanocrystal suspensions (NSs) were prepared using lab-scale bead-milling technology, by adjusting the type and amount of suspending agent, CA content, type of vehicle, and milling speed. CA-loaded NSs were characterized in terms of morphology, particle size, crystallinity, and in vitro dissolution pattern. Skin absorption of CA nanocrystals was evaluated using a vertical Franz diffusion cell mounted with porcine skin. In vivo skin irritation following topical application of high-payload NS was assessed in normal rats.

Results

The optimized NS system, composed of 10% (w/v) CA, 0.5% polyvinylpyrrolidone (PVP) K30 as steric stabilizer, and 89.5% of distilled water, was characterized as follows: spherical or elliptical in shape, 200 nm in size, with low crystallinity. The in vitro dissolution of AA or MA from NSs was markedly faster compared to raw material, under sink condition. Penetration of AA, MA, and AS in the porcine skin was markedly elevated using the high-payload NS formula, providing 5-, 4-, and 4.5-fold higher accumulation in skin layer, compared to that of the marketed cream formula (CA 1%, Madeca cream). Moreover, topical application of high-payload NS was tolerable, showing neither erythema nor oedema in normal rats.

Conclusion

The novel NS system is expected to be a virtuous approach for offering a better skin absorption of CA, without using an excess quantity of solubilizers.

Ability of mathematical models to predict human in vivo percutaneous penetration of steroids

Human skin is a common route for topical steroids to enter the body. To aid with risk management of therapeutic steroid usage, the US Environmental Protection Agency estimates percutaneous penetration using mathematical models. However, it is unclear how accurate are mathematical models in estimating percutaneous penetration/absorption of steroids. In this study, accuracy of predicted flux (penetration/absorption) by the main mathematical model used by the EPA, the Potts and Guy model based on in vitro data is compared to actual human in vivo data from our laboratory of percutaneous absorption of topical steroids. We focused on steroids due to the availability of steroid in vivo human data in our laboratory. For most steroids the flux was underestimated by a factor 10-60. However, within the group itself, there was an association between the Potts and Guy model and experimental human in vivo data (Pearson Correlation = 0.8925, p = 0.000041). Additionally, some physiochemical parameters used in the Potts and Guy equation, namely log Kp (Pearson Correlation = 0.7307, p = 0.0046) and molecular weight (Pearson correlation = -0.6807, p = 0.0105) correlated significantly with in vivo flux. Current mathematical models used in estimating percutaneous penetration/absorption did not accurately predict in vivo flux of steroids. Why? Proposed limitations to mathematical models currently used include: not accounting for volatility, lipid solubility, hydrogen bond effects, drug metabolism, as well as protein binding. Further research is needed in order to increase the predictive nature of such models for in vivo flux.

Physiologically-Based Pharmacokinetic Modeling to Support Determination of Bioequivalence for Dermatological Drug Products: Scientific and Regulatory Considerations

Physiologically-based pharmacokinetic (PBPK) modeling and simulation provides mechanism-based predictions of the pharmacokinetics of an active ingredient following its administration in humans. Dermal PBPK models describe the skin permeation and disposition of the active ingredient following the application of a dermatological product on the skin of virtual healthy and diseased human subjects. These models take into account information on product quality attributes, physicochemical properties of the active ingredient and skin (patho)physiology, and their interplay with each other. Regulatory and product development decision makers can leverage these quantitative tools to identify factors impacting local and systemic exposure. In the realm of generic drug products, the number of US Food and Drug Administratioin (FDA) interactions that use dermal PBPK modeling to support alternative bioequivalence (BE) approaches is increasing. In this report, we share scientific considerations on the development, verification and validation (V&V), and application of PBPK models within the context of a virtual BE assessment for dermatological drug products. We discuss the challenges associated with model V&V for these drug products stemming from the fact that target-site active ingredient concentrations are typically not measurable. Additionally, there are no established relationships between local and systemic PK profiles, when the latter are quantifiable. To that end, we detail a multilevel model V&V approach involving validation for the model of the drug product of interest coupled with the overall assessment of the modeling platform in use while leveraging in vitro and in vivo data related to local and systemic bioavailability.

Assessment of Vehicle Volatility and Deposition Layer Thickness in Skin Penetration Models

Systemic disposition of dermally applied chemicals is often formulation-dependent. Rapid evaporation of the vehicle can result in crystallization of active compounds, limiting their degree of skin penetration. In addition, the choice of vehicle can affect the permeant’s degree of penetration into the stratum corneum. The aim of this study is to build a predictive, mechanistic, dermal absorption model that accounts for vehicle-specific effects on the kinetics of permeant transport into skin. An existing skin penetration model is extended to explicitly include the effect of vehicle volatility over time. Using in vitro measurements of skin penetration by chemicals applied in both a saline and an ethanol solvent, the model is optimized to learn two vehicle-specific quantities: the solvent evaporation rate and the extent of permeant deposition into the upper stratum corneum immediately following application. The dermal disposition estimates of the trained model are subsequently compared against those of the original model using further in vitro measurements. The trained model showed a 1.5-fold improvement and a 19-fold improvement in overall goodness of fit among compounds tested in saline and ethanol solvents, respectively. The proposed model structure can thus form a basis for in vitro to in vivo extrapolations of dermal disposition for skin formulations containing volatile components.

Predicting viable skin concentration: Diffusional and convective drug transport

Viable skin drug transport is an important concept to consider as it can have a significant impact on the local concentration of a drug. The concentration becomes even more critical for toxicological issues when implementing different permeability enhancement techniques. For this reason, it is important to develop models that can predict drug transport in the viable skin. This paper expands upon previous capillary modeling by representing the convective transport of a solute that has permeated into the capillary loops. As a result, convective transport caused the concentration profile to plateau within the deeper dermal layers, effectively matching the trend of previous experimental data. Furthermore, the new model also has a significantly quicker transient profile as the time required to reach steady-state is five-fold faster than predicted in previous homogenous models.

Pharmacokinetics and pharmacodynamics in the treatment of cutaneous leishmaniasis - challenges and opportunities

Pharmacological efficacy is obtained when adequate concentrations of a potent drug reach the target site. In cutaneous leishmaniasis, a heterogeneous disease characterised by a variety of skin manifestations from simple nodules, skin discoloration, plaques to extensive disseminated forms, the parasites are found in the dermal layers of the skin. Treatment thus involves the release of the active compound from the formulation (administered either topically or systemically), it's permeation into the skin, accumulation by the local macrophages and further transport into the phagolysosome of the macrophage. The pharmacodynamic activity of a drug against the parasite is relatively straight forward to evaluate both in vivo and in vitro. The pharmacokinetic processes taking place inside the skin are more complex to elucidate due to the multi-lamellar structure of the skin, heterogeneous distribution of drugs within the tissue, the difficulty of accessing the site of infection complicating sampling and the lack of surrogate markers reflecting the activity of a drug in the skin. This review will discuss the difficulties encountered when investigating drug distribution, PK PD relationships and efficacy in the skin with a focus on cutaneous leishmaniasis treatment.

The strategies for the modelling of the passive mass transport through porous membranes: Applicability to transdermal delivery systems

The paper concerns the modelling of the passive solute transport through porous membranes. A general scheme for the mass transport has been developed upon the mixed diffusion-advection-reaction model. The passive advection has been introduced as a certain simplification of the Navier-Stokes problem, involving a pressure gradient-induced creeping flow of an incompressible Newtonian fluid. Nine scenarios for the drug transport process have been tested versus two experimental datasets acquired earlier (photoacoustic depth-profiling and contact angle surface wettability techniques) for the characterization of bulk and interfacial processes in a model pharmaceutical system: the synthetic dodecanol-collodion porous membrane in contact with a photodegradable pigment dithranol. The scenarios considered include three mass transport models (the diffusion-advection-reaction, diffusion-advection and diffusion-reaction models) under three distinct types of the lower (the donor/acceptor interface) boundary conditions: the Dirichlet-type instantaneous source, the Dirichlet-type interface relaxation, and the Neumann-type concentration gradient. The results obtained indicate a considerable agreement between the experimental data and predictions of the diffusion-reaction and the general models for long times, however, some deviations were exhibited at the initial stages of the permeation process. It is considered, that the discrepancies originate from a specific penetrant behaviour at the interfaces, which violates boundary transfer schemes classically employed for the mass transport phenomena quantification. Moreover, an additional mixing process taking place close to the interface related to the liquid flow driven by the surface tension gradients (so-called classic and thermal Marangoni effect) could play a still underestimated role in the trans-interfacial mass transport.

Modeling drug transport within the viable skin - a review

INTRODUCTION

In the past, mathematical modeling of the transport of transdermal drugs has been primarily focused on the stratum corneum. However, the development of pharmaceutical technologies, such as chemical enhancers, iontophoresis, and microneedles, has led to two outcomes; an increase in permeability in the stratum corneum or the ability to negate the layer entirely. As a result, these outcomes have made the transport of a solute in the viable skin far more critical when studying transdermal drug delivery.

AREAS COVERED

The review will explicitly show the various attempts to model drug transport within the viable skin. Furthermore, a brief review will be conducted on the different models that explain stratum corneum transport, microneedle dynamics and estimation of the diffusion coefficient.

EXPERT OPINION

Future development of mathematical models requires the focus to be changed from traditional diffusion-based tissue models to more sophisticated three-dimensional models that incorporate the physiology of the skin.

Diffusion through skin in the light of a fractional derivative approach: progress and challenges

This review is focussed on modelling the transport processes of different drugs across the intact human skin by introducing a memory formalism based on the fractional derivative approach. The fundamental assumption of the classic transport equation in the light of the Fick's law is that the skin barrier behaves as a pseudo-homogeneous membrane and that its properties, summarized by the diffusion coefficient D,  do not vary with time and position. This assumption does not hold in the case of a highly heterogeneous system as the skin is, whose outermost layer (the stratum corneum) is comprised of a multi-layered structure of keratinocytes embedded in a lamellar matrix of hydrophobic lipids, followed by the dermis that contains a network of capillaries that connect to the systemic circulation. A possible way to overcome these difficulties resides in the introduction of mathematical models which involve fractional derivatives to describe complex systems with interactions in space and time, following the model originally developed by Caputo in order to consider the memory effects in materials. Although the introduction of fractional derivatives to model memory effects is completely phenomenological, i.e., characterized by a single parameter, i.e., the fractional derivative order [Formula: see text] a number of authors have found that this approach can provide a better comparison to experimental data and that this technique may be alternative to integer-order derivative models. In this review, we aim to summarize some our recent results, concerning the transport of different diffusing compounds of different structural complexity across the intact skin.

The application of label-free imaging technologies in transdermal research for deeper mechanism revealing

The penetration behavior of topical substances in the skin not only relates to the transdermal delivery efficiency but also involves the safety and therapeutic effect of topical products, such as sunscreen and hair growth products. Researchers have tried to illustrate the transdermal process with diversified theories and technologies. Directly observing the distribution of topical substances on skin by characteristic imaging is the most convincing approach. Unfortunately, fluorescence labeling imaging, which is commonly used in biochemical research, is limited for transdermal research for most topical substances with a molecular mass less than 500 Da. Label-free imaging technologies possess the advantages of not requiring any macromolecular dyes, no tissue destruction and an extensive substance detection capability, which has enabled rapid development of such technologies in recent years and their introduction to biological tissue analysis, such as skin samples. Through the specific identification of topical substances and endogenous tissue components, label-free imaging technologies can provide abundant tissue distribution information, enrich theoretical and practical guidance for transdermal drug delivery systems. In this review, we expound the mechanisms and applications of the most popular label-free imaging technologies in transdermal research at present, compare their advantages and disadvantages, and forecast development prospects.

Methods to Evaluate Skin Penetration In Vitro

Dermal and transdermal drug therapy is increasing in importance nowadays in drug development. To completely utilize the potential of this administration route, it is necessary to optimize the drug release and skin penetration measurements. This review covers the most well-known and up-to-date methods for evaluating the cutaneous penetration of drugs in vitro as a supporting tool for pharmaceutical research scientists in the early stage of drug development. The aim of this article is to present various experimental models used in dermal/transdermal research and summarize the novel knowledge about the main in vitro methods available to study skin penetration. These techniques are: Diffusion cell, skin-PAMPA, tape stripping, two-photon microscopy, confocal laser scanning microscopy, and confocal Raman microscopic method.

Percutaneous penetration of drugs applied in transdermal delivery systems: an in vivo based approach for evaluating computer generated penetration models

Human skin is a viable pathway for administration of therapeutics. Transdermal delivery systems (TDS) have been approved by the US-FDA since 1981. To enable the risk assessment of dermal exposure, predictive mathematical models are used. In this work the accuracy of predicted flux of the models is compared to experimental human in vivo data of drugs applied in US-FDA approved TDS. A database of pharmacokinetic data of drugs applied in TDS was used and updated. Three mathematical models (QSAR) were used to calculate predicted fluxes, and compared to the human in vivo data. For more than half of the drugs applied in TDS, the predicted flux by the mathematical models was an underestimation compared to the flux calculated with the experimental in vivo data. The flux was over- or underestimated by a factor 10-100. All mathematical models were significantly correlated with the experimental in vivo data. The process of percutaneous penetration has several influencing factors, TDS minimize some of these factors. Limitations are discussed. Further research is needed, with a focus on validation and standardization of the permeability coefficient. This work offers observations that should give a stimulus for refinement on the appropriate usage and limitations of mathematical models.

Space- and time-resolved investigation on diffusion kinetics of human skin following macromolecule delivery by microneedle arrays

Microscale medical devices are being developed for targeted skin delivery of vaccines and the extraction of biomarkers, with the potential to revolutionise healthcare in both developing and developed countries. The effective clinical development of these devices is dependent on understanding the macro-molecular diffusion properties of skin. We hypothesised that diffusion varied according to specific skin layers. Using three different molecular weights of rhodamine dextran (RD) (MW of 70, 500 and 2000 kDa) relevant to the vaccine and therapeutic scales, we deposited molecules to a range of depths (0-300 µm) in ex vivo human skin using the Nanopatch device. We observed significant dissipation of RD as diffusion with 70 and 500 kDa within the 30 min timeframe, which varied with MW and skin layer. Using multiphoton microscopy, image analysis and a Fick's law analysis with 2D cartesian and axisymmetric cylindrical coordinates, we reported experimental trends of epidermal and dermal diffusivity values ranging from 1-8 µm2 s-1 to 1-20 µm2 s-1 respectively, with a significant decrease in the dermal-epidermal junction of 0.7-3 µm2 s-1. In breaching the stratum corneum (SC) and dermal-epidermal junction barriers, we have demonstrated practical application, delivery and targeting of macromolecules to both epidermal and dermal antigen presenting cells, providing a sound knowledge base for future development of skin-targeting clinical technologies in humans.