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

Delivery of small hydrophilic molecules across the stratum corneum: Identification of model systems and parameters to study topical delivery of free amino acids

Free amino acids (FAAs) constitute the largest component (∼40 %) of the so-called natural moisturizing factors of the skin. Their level declines in dry skin conditions and one strategy to overcome this problem may involve the topical delivery of FAAs through appropriate strategy. The objective of the present study was therefore to identify alternative skin models and study the corneocyte-water partition coefficients (KCOR/W) and permeation coefficient (KP) of 18 FAAs. The KCOR/W was studied using standard protocols and the permeation studies were conducted using Franz diffusion cell. The results indicate that the FAAs have high partitioning behavior to the corneocytes. The KCOR/W values of the human COR and that of pig ear skin were better correlated with each other than that of keratin isolated from chicken feathers. The presence of lipid in the stratum corneum (SC), initial concentration of the FAAs, and permeation enhancers affect the KCOR/W. The FAAs have low permeation into the SC which suggests the need for permeation enhancers in designing dosage form containing these compounds. Even though the investigated mathematical models show good prediction of the Kp values, better prediction could be obtained by considering factors such as the possible entrapment of the FAAs by the CORs.

Understanding Drug Skin Permeation Enhancers Using Molecular Dynamics Simulations

Our skin constitutes an effective permeability barrier that protects the body from exogenous substances but concomitantly severely limits the number of pharmaceutical drugs that can be delivered transdermally. In topical formulation design, chemical permeation enhancers (PEs) are used to increase drug skin permeability. In vitro skin permeability experiments can measure net effects of PEs on transdermal drug transport, but they cannot explain the molecular mechanisms of interactions between drugs, permeation enhancers, and skin structure, which limits the possibility to rationally design better new drug formulations. Here we investigate the effect of the PEs water, lauric acid, geraniol, stearic acid, thymol, ethanol, oleic acid, and eucalyptol on the transdermal transport of metronidazole, caffeine, and naproxen. We use atomistic molecular dynamics (MD) simulations in combination with developed molecular models to calculate the free energy difference between 11 PE-containing formulations and the skin's barrier structure. We then utilize the results to calculate the final concentration of PEs in skin. We obtain an RMSE of 0.58 log units for calculated partition coefficients from water into the barrier structure. We then use the modified PE-containing barrier structure to calculate the PEs' permeability enhancement ratios (ERs) on transdermal metronidazole, caffeine, and naproxen transport and compare with the results obtained from in vitro experiments. We show that MD simulations are able to reproduce rankings based on ERs. However, strict quantitative correlation with experimental data needs further refinement, which is complicated by significant deviations between different measurements. Finally, we propose a model for how to use calculations of the potential of mean force of drugs across the skin's barrier structure in a topical formulation design.

Microneedle arrays integrated with microfluidic systems: Emerging applications and fluid flow modeling

Microneedle arrays are patches of needles at micro- and nano-scale, which are competent and versatile technologies that have been merged with microfluidic systems to construct more capable devices for biomedical applications, such as drug delivery, wound healing, biosensing, and sampling body fluids. In this paper, several designs and applications are reviewed. In addition, modeling approaches used in microneedle designs for fluid flow and mass transfer are discussed, and the challenges are highlighted.

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.

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.

Inverse Mechanistic Modeling of Transdermal Drug Delivery for Fast Identification of Optimal Model Parameters

Transdermal drug delivery systems are a key technology to administer drugs with a high first-pass effect in a non-invasive and controlled way. Physics-based modeling and simulation are on their way to become a cornerstone in the engineering of these healthcare devices since it provides a unique complementarity to experimental data and additional insights. Simulations enable to virtually probe the drug transport inside the skin at each point in time and space. However, the tedious experimental or numerical determination of material properties currently forms a bottleneck in the modeling workflow. We show that multiparameter inverse modeling to determine the drug diffusion and partition coefficients is a fast and reliable alternative. We demonstrate this strategy for transdermal delivery of fentanyl. We found that inverse modeling reduced the normalized root mean square deviation of the measured drug uptake flux from 26 to 9%, when compared to the experimental measurement of all skin properties. We found that this improved agreement with experiments was only possible if the diffusion in the reservoir holding the drug was smaller than the experimentally measured diffusion coefficients suggested. For indirect inverse modeling, which systematically explores the entire parametric space, 30,000 simulations were required. By relying on direct inverse modeling, we reduced the number of simulations to be performed to only 300, so a factor 100 difference. The modeling approach's added value is that it can be calibrated once in-silico for all model parameters simultaneously by solely relying on a single measurement of the drug uptake flux evolution over time. We showed that this calibrated model could accurately be used to simulate transdermal patches with other drug doses. We showed that inverse modeling is a fast way to build up an accurate mechanistic model for drug delivery. This strategy opens the door to clinically ready therapy that is tailored to patients.

Formulation and evaluation of transdermal nanogel for delivery of artemether

Artemether (ART) is second to artesunate in being the most widely used derivatives of artemisinin in combination therapy of malaria. Nanostructured lipid carrier (NLC) formulations were prepared following our previous report using optimized ART concentration of 0.25 g dissolved in 5% w/v mixture of solid (Gelucire 43/01 and Phospholipon 85G) and liquid (Transcutol) lipids at 90 °C. An aqueous surfactant phase at 90 °C was added (dropwise) under magnetic stirring (1000 rpm) for 5 min. The pre-emulsion was speedily homogenized at 28,000 rpm for 15 min and further probe sonicated at 60% amplitude (15 min). Resultant sample was cooled at room temperature and frozen at - 80 °C prior to lyophilization. The freeze-dried sample was used for solid-state characterization as well as in the formulation of transdermal nanogels using three polymers (Carbopol 971P, Poloxamer 407, and Prosopis africana peel powder) to embed the ART-NLC, using ethanol as a penetration enhancer. Transdermal ART-nanogels were characterized accordingly (physical examination, pH, drug content, rheology, spreadability, stability, particle size and morphology, skin irritation, in vitro and ex vivo skin permeation, and analysis of permeation data), P < 0.05. Results indicated that ART nanogels showed good encapsulation, drug release, pH-dependent swelling, stability, and tolerability. Overall, ART nanogels prepared from Poloxamer 407 showed the most desirable drug permeation, pH, swellability, spreadability, viscosity, and transdermal antiplasmodial properties superior to PAPP-ANG > C971P-ANG. A two-patch/week concurrent application of the studied nanogels could offer 100% cure of malaria as a lower-dose (50 mg ART) patient-friendly regimen devoid of the drug's many side effects.

Influence of keratolytics on cutaneous pharmacokinetics of glucocorticoids

BACKGROUND

Keratolytics are often used to accelerate and improve the therapeutic response of hyperkeratotic dermatoses. Keratolytics are a chemically inhomogeneous group of substances and substance mixtures that clinically lead to a decrease in symptoms of a cornification disorder, but mediate different effects. Thus, keratolytic, keratoplastic, keratoemulsifying and keratodiluting effects are distinguished. The physicochemical effects or pharmacological efficacy of the respective keratolytics result in different mechanism as well as risks with regard to local or systemic compatibility. Until now, only little attention has been paid upon selection of keratolytics to the immediate consequences regarding diffusion conditions and pharmacokinetics of sequentially applied topicals, in particular of glucocorticoids.

PATIENTS AND METHODS

This paper deals with the influence of keratolytics on the penetration-time profile of betamethasone dipropionate in sequential application. For this purpose, cutaneous bioavailability was investigated with the Franz chamber test using a tritium-labeled drug depending on the previous application of a keratolytic agent. Comparative data analyses were performed.

RESULTS

It was shown that keratoplastic substances significantly promote diffusion of the glucocorticoid. Keratoemulsifying substance mixtures exert no relevant effects in this regard, while keratodiluting substance mixtures inhibit penetration.

CONCLUSIONS

The targeted selection of a keratolytic can optimize the therapeutic effect and influence the bioavailability of sequentially applied topicals.

Predicting Transdermal Fentanyl Delivery Using Mechanistic Simulations for Tailored Therapy

Transdermal drug delivery is a key technology for administering drugs. However, most devices are “one-size-fits-all”, even though drug diffusion through the skin varies significantly from person-to-person. For next-generation devices, personalization for optimal drug release would benefit from an augmented insight into the drug release and percutaneous uptake kinetics. Our objective was to quantify the changes in transdermal fentanyl uptake with regards to the patient’s age and the anatomical location where the patch was placed. We also explored to which extent the drug flux from the patch could be altered by miniaturizing the contact surface area of the patch reservoir with the skin. To this end, we used validated mechanistic modeling of fentanyl diffusion, storage, and partitioning in the epidermis to quantify drug release from the patch and the uptake within the skin. A superior spatiotemporal resolution compared to experimental methods enabled in-silico identification of peak concentrations and fluxes, and the amount of stored drug and bioavailability. The patients’ drug uptake showed a 36% difference between different anatomical locations after 72 h, but there was a strong interpatient variability. With aging, the drug uptake from the transdermal patch became slower and less potent. A 70-year-old patient received 26% less drug over the 72-h application period, compared to an 18-year-old patient. Additionally, a novel concept of using micron-sized drug reservoirs was explored in silico. These reservoirs induced a much higher local flux (µg cm^-2 h^-1) than conventional patches. Up to a 200-fold increase in the drug flux was obtained from these small reservoirs. This effect was mainly caused by transverse diffusion in the stratum corneum, which is not relevant for much larger conventional patches. These micron-sized drug reservoirs open new ways to individualize reservoir design and thus transdermal therapy. Such computer-aided engineering tools also have great potential for in-silico design and precise control of drug delivery systems. Here, the validated mechanistic models can serve as a key building block for developing digital twins for transdermal drug delivery systems.

Mathematical Modelling, Simulation and Optimisation of Microneedles for Transdermal Drug Delivery: Trends and Progress

In the last two decades, microneedles (MNs) have received significant interest due to their potential for painless transdermal drug delivery (TDD) and minimal skin damage. MNs have found applications in a range of research and development areas in drug delivery. They have been prepared using a variety of materials and fabrication techniques resulting in MN arrays with different dimensions, shapes, and geometries for delivery of a variety of drug molecules. These parameters play crucial roles in determining the drug release profiles from the MNs. Developing mathematical modelling, simulation, and optimisation techniques is vital to achieving the desired MN performances. These will then be helpful for pharmaceutical and biotechnological industries as well as professionals working in the field of regulatory affairs focusing on MN based TDD systems. This is because modelling has a great potential to reduce the financial and time cost of both the MNs' studies and manufacturing. For example, a number of robust mathematical models for predicting the performance of the MNs in vivo have emerged recently which incorporate the roles of the structural and mechanical properties of the skin. In addressing these points, this review paper aims to highlight the current status of the MN modelling research, in particular, the modelling, simulation and optimisation of the systems for drug delivery. The theoretical basis for the simulation of MN enhanced diffusion is discussed within this paper. Thus, this review paper provides a better understanding of the modelling of the MN mediated drug delivery process.

Salt water and skin interactions: new lines of evidence

In Health Resort Medicine, both balneotherapy and thalassotherapy, salt waters and their peloids, or mud products are mainly used to treat rheumatic and skin disorders. These therapeutic agents act jointly via numerous mechanical, thermal, and chemical mechanisms. In this review, we examine a new mechanism of action specific to saline waters. When topically administered, this water rich in sodium and chloride penetrates the skin where it is able to modify cellular osmotic pressure and stimulate nerve receptors in the skin via cell membrane ion channels known as "Piezo" proteins. We describe several models of cutaneous adsorption/desorption and penetration of dissolved ions in mineral waters through the skin (osmosis and cell volume mechanisms in keratinocytes) and examine the role of these resources in stimulating cutaneous nerve receptors. The actions of salt mineral waters are mediated by a mechanism conditioned by the concentration and quality of their salts involving cellular osmosis-mediated activation/inhibition of cell apoptotic or necrotic processes. In turn, this osmotic mechanism modulates the recently described mechanosensitive piezoelectric channels.

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

PURPOSE

Sebum is an important shunt pathway for transdermal permeation and targeted delivery, but there have been limited studies on its permeation properties. Here we report a measurement and modelling study of solute partition to artificial sebum.

METHODS

Equilibrium experiments were carried out for the sebum-water partition coefficients of 23 neutral, cationic and anionic compounds at different pH.

RESULTS

Sebum-water partition coefficients not only depend on the hydrophobicity of the chemical but also on pH. As pH increases from 4.2 to 7.4, the partition of cationic chemicals to sebum increased rapidly. This appears to be due to increased electrostatic attraction between the cationic chemical and the fatty acids in sebum. Whereas for anionic chemicals, their sebum partition coefficients are negligibly small, which might result from their electrostatic repulsion to fatty acids. Increase in pH also resulted in a slight decrease of sebum partition of neutral chemicals.

CONCLUSIONS

Based on the observed pH impact on the sebum-water partition of neutral, cationic and anionic compounds, a new quantitative structure-property relationship (QSPR) model has been proposed. This mathematical model considers the hydrophobic interaction and electrostatic interaction as the main mechanisms for the partition of neutral, cationic and anionic chemicals to sebum.

How Sensitive Are Transdermal Transport Predictions by Microscopic Stratum Corneum Models to Geometric and Transport Parameter Input?

While predictive models of transdermal transport have the potential to reduce human and animal testing, microscopic stratum corneum (SC) model output is highly dependent on idealized SC geometry, transport pathway (transcellular vs. intercellular), and penetrant transport parameters (e.g., compound diffusivity in lipids). Most microscopic models are limited to a simple rectangular brick-and-mortar SC geometry and do not account for variability across delivery sites, hydration levels, and populations. In addition, these models rely on transport parameters obtained from pure theory, parameter fitting to match in vivo experiments, and time-intensive diffusion experiments for each compound. In this work, we develop a microscopic finite element model that allows us to probe model sensitivity to variations in geometry, transport pathway, and hydration level. Given the dearth of experimentally-validated transport data and the wide range in theoretically-predicted transport parameters, we examine the model's response to a variety of transport parameters reported in the literature. Results show that model predictions are strongly dependent on all aforementioned variations, resulting in order-of-magnitude differences in lag times and permeabilities for distinct structure, hydration, and parameter combinations. This work demonstrates that universally predictive models cannot fully succeed without employing experimentally verified transport parameters and individualized SC structures.

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.

In-Silico Skin Model: A Multiscale Simulation Study of Drug Transport

Accurate in-silico models are required to predict the release of drug molecules through skin in order to supplement the in-vivo experiments for faster development/testing of drugs. The upper most layer of the skin, stratum corneum (SC), offers the main resistance for permeation of actives. Most of the SC's molecular level models comprise cholesterol and phospholipids only, which is far from reality. In this study we have implemented a multiscale modeling framework to obtain the release profile of three drugs, namely, caffeine, fentanyl, and naphthol, through skin SC. We report for the first time diffusion of drugs through a realistic skin molecular model comprised of ceramides, cholesterol, and free fatty acid. The diffusion coefficients of drugs in the SC lipid matrix were determined from multiple constrained molecular dynamics simulations. The calculated diffusion coefficients were then used in the macroscopic models to predict the release profiles of drugs through the SC. The obtained release profiles were in good agreement with available experimental data. The partition coefficient exhibits a greater effect on the release profiles. The reported multiscale modeling framework would provide insight into the delivery mechanisms of the drugs through the skin and shall act as a guiding tool in performing targeted experiments to come up with a suitable delivery system.

Surging footprints of mathematical modeling for prediction of transdermal permeability

In vivo skin permeation studies are considered gold standard but are difficult to perform and evaluate due to ethical issues and complexity of process involved. In recent past, a useful tool has been developed by combining the computational modeling and experimental data for expounding biological complexity. Modeling of percutaneous permeation studies provides an ethical and viable alternative to laboratory experimentation. Scientists are exploring complex models in magnificent details with advancement in computational power and technology. Mathematical models of skin permeability are highly relevant with respect to transdermal drug delivery, assessment of dermal exposure to industrial and environmental hazards as well as in developing fundamental understanding of biotransport processes. Present review focuses on various mathematical models developed till now for the transdermal drug delivery along with their applications.

Comparison of protocols measuring diffusion and partition coefficients in the stratum corneum

Partition (K) and diffusion (D) coefficients are important to measure for the modelling of skin penetration of chemicals through the stratum corneum (SC). We compared the feasibility of three protocols for the testing of 50 chemicals in our main studies, using three cosmetics-relevant model chemicals with a wide range of logP values. Protocol 1: SC concentration-depth profile using tape-stripping (measures KSC/v and DSC /HSC2 , where HSC is the SC thickness); Protocol 2A: incubation of isolated SC with chemical (direct measurement of KSC/v only) and Protocol 2B: diffusion through isolated SC mounted on a Franz cell (measures KSC/v and DSC /HSC2 , and is based on Fick's laws). KSC/v values for caffeine and resorcinol using Protocol 1 and 2B were within 30% of each other, values using Protocol 2A were ~two-fold higher, and all values were within 10-fold of each other. Only indirect determination of KSC/v by Protocol 2B was different from the direct measurement of KSC/v by Protocol 2A and Protocol 1 for 7-EC. The variability of KSC/v for all three chemicals using Protocol 2B was higher compared to Protocol 1 and 2A. DSC /HSC2 values for the three chemicals were of the same order of magnitude using all three protocols. Additionally, using Protocol 1, there was very little difference between parameters measured in pig and human SC. In conclusion, KSC/v, and DSC values were comparable using different methods. Pig skin might be a good surrogate for human skin for the three chemicals tested. Copyright © 2017 The Authors Journal of Applied Toxicology published by John Wiley & Sons Ltd.

Effect of stratum corneum heterogeneity, anisotropy, asymmetry and follicular pathway on transdermal penetration

The impact of the complex structure of the stratum corneum on transdermal penetration is not yet fully described by existing models. A quantitative and thorough study of skin permeation is essential for chemical exposure assessment and transdermal delivery of drugs. The objective of this study is to analyze the effects of heterogeneity, anisotropy, asymmetry, follicular diffusion, and location of the main barrier of diffusion on percutaneous permeation. In the current study, the solution of the transient diffusion through a two-dimensional-anisotropic brick-and-mortar geometry of the stratum corneum is obtained using the commercial finite element program COMSOL Multiphysics. First, analytical solutions of an equivalent multilayer geometry are used to determine whether the lipids or corneocytes constitute the main permeation barrier. Also these analytical solutions are applied for validations of the finite element solutions. Three illustrative compounds are analyzed in these sections: diethyl phthalate, caffeine and nicotine. Then, asymmetry with depth and follicular diffusion are studied using caffeine as an illustrative compound. The following findings are drawn from this study: the main permeation barrier is located in the lipid layers; the flux and lag time of diffusion through a brick-and-mortar geometry are almost identical to the values corresponding to a multilayer geometry; the flux and lag time are affected when the lipid transbilayer diffusivity or the partition coefficients vary with depth, but are not affected by depth-dependent corneocyte diffusivity; and the follicular contribution has significance for low transbilayer lipid diffusivity, especially when flux between the follicle and the surrounding stratum corneum is involved. This study demonstrates that the diffusion is primarily transcellular and the main barrier is located in the lipid layers.

Development of a Two-Dimensional Model for Predicting Transdermal Permeation with the Follicular Pathway: Demonstration with a Caffeine Study

Purpose

The development of a new two-dimensional (2D) model to predict follicular permeation, with integration into a recently reported multi-scale model of transdermal permeation is presented.

Methods

The follicular pathway is modelled by diffusion in sebum. The mass transfer and partition properties of solutes in lipid, corneocytes, viable dermis, dermis and systemic circulation are calculated as reported previously [Pharm Res 33 (2016) 1602]. The mass transfer and partition properties in sebum are collected from existing literature. None of the model input parameters was fit to the clinical data with which the model prediction is compared.

Results

The integrated model has been applied to predict the published clinical data of transdermal permeation of caffeine. The relative importance of the follicular pathway is analysed. Good agreement of the model prediction with the clinical data has been obtained. The simulation confirms that for caffeine the follicular route is important; the maximum bioavailable concentration of caffeine in systemic circulation with open hair follicles is predicted to be 20% higher than that when hair follicles are blocked.

Conclusions

The follicular pathway contributes to not only short time fast penetration, but also the overall systemic bioavailability. With such in silico model, useful information can be obtained for caffeine disposition and localised delivery in lipid, corneocytes, viable dermis, dermis and the hair follicle. Such detailed information is difficult to obtain experimentally.

Direct estimation of the permeation of topical excipients through artificial membranes and human skin with non-invasive Terahertz time-domain techniques

BACKGROUND

Drug permeation through skin, or a synthetic membrane, from locally acting pharmaceutical products can be influenced by the permeation behaviour of pharmaceutical excipients.

OBJECTIVE

Terahertz time-domain technology is investigated as a non-invasive method for a direct and accurate measurement of excipients permeation through synthetic membranes or human skin.

METHODS

A series of in-vitro release and skin permeation experiments of liquid excipients (e.g. propylene glycol and polyethylene glycol 400) has been conducted with vertical diffusion cells. The permeation profiles of excipients through different synthetic membranes or skin were obtained using Terahertz pulses providing a direct measurement. Corresponding permeation flux and permeability coefficient values were calculated based on temporal changes of the terahertz pulses.

RESULTS

The influence of different experimental conditions, such as the polarity of the membrane and the viscosity of the permeant, was assessed in release experiments. Specific transmembrane flux values of those excipients were directly calculated with statistical differences between cases. Finally, an attempt to estimate the skin permeation of propylene glycol with this technique was also achieved. All these permeation results were likely comparable to those obtained by other authors with usual analytical techniques.

CONCLUSION

Terahertz time-domain technology is shown to be a suitable technique for an accurate and non-destructive measurement of the permeation of liquid substances through different synthetic membranes or even human skin.

Anatomically Detailed and Large-Scale Simulations Studying Synapse Loss and Synchrony Using NeuroBox

The morphology of neurons and networks plays an important role in processing electrical and biochemical signals. Based on neuronal reconstructions, which are becoming abundantly available through databases such as NeuroMorpho.org, numerical simulations of Hodgkin-Huxley-type equations, coupled to biochemical models, can be performed in order to systematically investigate the influence of cellular morphology and the connectivity pattern in networks on the underlying function. Development in the area of synthetic neural network generation and morphology reconstruction from microscopy data has brought forth the software tool NeuGen. Coupling this morphology data (either from databases, synthetic, or reconstruction) to the simulation platform UG 4 (which harbors a neuroscientific portfolio) and VRL-Studio, has brought forth the extendible toolbox NeuroBox. NeuroBox allows users to perform numerical simulations on hybrid-dimensional morphology representations. The code basis is designed in a modular way, such that e.g., new channel or synapse types can be added to the library. Workflows can be specified through scripts or through the VRL-Studio graphical workflow representation. Third-party tools, such as ImageJ, can be added to NeuroBox workflows. In this paper, NeuroBox is used to study the electrical and biochemical effects of synapse loss vs. synchrony in neurons, to investigate large morphology data sets within detailed biophysical simulations, and used to demonstrate the capability of utilizing high-performance computing infrastructure for large scale network simulations. Using new synapse distribution methods and Finite Volume based numerical solvers for compartment-type models, our results demonstrate how an increase in synaptic synchronization can compensate synapse loss at the electrical and calcium level, and how detailed neuronal morphology can be integrated in large-scale network simulations.

A microneedle-based method for the characterization of diffusion in skin tissue using doxorubicin as a model drug

Hollow microneedles can overcome the stratum corneum (SC) barrier and deposit a compound directly into the viable epidermis or the dermis, unlike adhesive patches that rely on drug diffusion across the SC. The traditional one-dimensional methods used to study the diffusivity of drugs across the skin layers are not very accurate for hollow microneedles, since the ejection of compounds out of microneedle lumens resembles a point-source spreading in all directions and is highly dependent on injection depth. This paper presents a technique that is useful for studying drug injection using hollow microneedles at various depths below the SC. This technique uses confocal microscopy to image the distribution of a fluorescent compound in the skin after injection. The fluorescence distribution in the skin is observed over time and applied to a spherical Gaussian diffusion model for limited source diffusion to determine the diffusion coefficient of the compound in the skin. Applied to freshly excised pig skin, the diffusion coefficient for the anti-cancer drug doxorubicin was measured as 4.61 × 10(-9) cm(2)/s, while the diffusion coefficient in previously refrigerated or frozen pig skin was 1.31 × 10(-8) cm(2)/s and 4.21 × 10(-8) cm(2)/s, respectively. Our data suggests that skin storage conditions can substantially alter the diffusion of drugs. The use of refrigerated and, even more so, previously frozen skin should be avoided for quantitative transdermal drug delivery studies.

Implementation of a single quad MS detector in high-throughput transdermal research of plant extracts

In this study, a new type of single quadrupole mass spectrometric detector was implemented in transdermal research. The local skin pharmacokinetic properties of the plant N-alkylamides (NAAs) pellitorine and anacycline, present in an Anacyclus pyrethrum extract, and spilanthol, present in a Spilanthes acmella extract were investigated. This single quad MS detection method showed great advantages compared to the traditional UV detector. The NAAs could be identified and quantified in the samples with an ultra performance liquid chromatography (UPLC)-single quad MS detection system, even if they were not separated, which is a requirement when using an UV-detector. Another advantage of the UPLC-MS system is that lower limit of detection values could be obtained allowing a more accurate and precise determination of the experimental lag time in the in vitro skin permeation experiments. To conclude, this single quad MS detector coupled to UPLC is a useful analytical tool with improved performance compared to high performance liquid chromatography (HPLC)-UV for biomedical-pharmaceutical purposes in transdermal research.

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.

Quantitative detection of caffeine in human skin by confocal Raman spectroscopy--A systematic in vitro validation study

For rational development and evaluation of dermal drug delivery, the knowledge of rate and extent of substance penetration into the human skin is essential. However, current analytical procedures are destructive, labor intense and lack a defined spatial resolution. In this context, confocal Raman microscopy bares the potential to overcome current limitations in drug depth profiling. Confocal Raman microscopy already proved its suitability for the acquisition of qualitative penetration profiles, but a comprehensive investigation regarding its suitability for quantitative measurements inside the human skin is still missing. In this work, we present a systematic validation study to deploy confocal Raman microscopy for quantitative drug depth profiling in human skin. After we validated our Raman microscopic setup, we successfully established an experimental procedure that allows correlating the Raman signal of a model drug with its controlled concentration in human skin. To overcome current drawbacks in drug depth profiling, we evaluated different modes of peak correlation for quantitative Raman measurements and offer a suitable operating procedure for quantitative drug depth profiling in human skin. In conclusion, we successfully demonstrate the potential of confocal Raman microscopy for quantitative drug depth profiling in human skin as valuable alternative to destructive state-of-the-art techniques.

Development of artemether-loaded nanostructured lipid carrier (NLC) formulation for topical application

NLC topical formulation as an alternative to oral and parenteral (IM) delivery of artemether (ART), a poorly water-soluble drug was designed. A Phospholipon 85G-modified Gelucire 43/01 based NLC formulation containing 75% Transcutol was chosen from DSC studies and loaded with gradient concentration of ART (100-750 mg). ART-loaded NLCs were stable (-22 to -40 mV), polydispersed (0.4-0.7) with d90 size distribution range of 247-530 nm without microparticles up to one month of storage. The encapsulation efficiency (EE%) for ART in the NLC was concentration independent as 250 mg of ART loading achieved ∼61%. DSC confirmed molecular dispersion of ART due to low matrix crystallinity (0.028J/g). Ex vivo study showed detectable ART amounts after 20h which gradually increased over 48h achieving ∼26% cumulative amount permeated irrespective of the applied dose. This proves that ART permeates excised human epidermis, where the current formulation served as a reservoir to gradually control drug release over an extended period of time. Full thickness skin study therefore may confirm if this is a positive signal to hope for a topical delivery system of ART.

1D-3D hybrid modeling-from multi-compartment models to full resolution models in space and time

Investigation of cellular and network dynamics in the brain by means of modeling and simulation has evolved into a highly interdisciplinary field, that uses sophisticated modeling and simulation approaches to understand distinct areas of brain function. Depending on the underlying complexity, these models vary in their level of detail, in order to cope with the attached computational cost. Hence for large network simulations, single neurons are typically reduced to time-dependent signal processors, dismissing the spatial aspect of each cell. For single cell or networks with relatively small numbers of neurons, general purpose simulators allow for space and time-dependent simulations of electrical signal processing, based on the cable equation theory. An emerging field in Computational Neuroscience encompasses a new level of detail by incorporating the full three-dimensional morphology of cells and organelles into three-dimensional, space and time-dependent, simulations. While every approach has its advantages and limitations, such as computational cost, integrated and methods-spanning simulation approaches, depending on the network size could establish new ways to investigate the brain. In this paper we present a hybrid simulation approach, that makes use of reduced 1D-models using e.g., the NEURON simulator—which couples to fully resolved models for simulating cellular and sub-cellular dynamics, including the detailed three-dimensional morphology of neurons and organelles. In order to couple 1D- and 3D-simulations, we present a geometry-, membrane potential- and intracellular concentration mapping framework, with which graph- based morphologies, e.g., in the swc- or hoc-format, are mapped to full surface and volume representations of the neuron and computational data from 1D-simulations can be used as boundary conditions for full 3D simulations and vice versa. Thus, established models and data, based on general purpose 1D-simulators, can be directly coupled to the emerging field of fully resolved, highly detailed 3D-modeling approaches. We present the developed general framework for 1D/3D hybrid modeling and apply it to investigate electrically active neurons and their intracellular spatio-temporal calcium dynamics.

Numerical modelling of transdermal delivery from matrix systems: parametric study and experimental validation with silicone matrices

A model is presented for transdermal drug delivery from single-layered silicone matrix systems. The work is based on our previous results that, in particular, extend the well-known Higuchi model. Recently, we have introduced a numerical transient model describing matrix systems where the drug dissolution can be non-instantaneous. Furthermore, our model can describe complex interactions within a multi-layered matrix and the matrix to skin boundary. The power of the modelling approach presented here is further illustrated by allowing the possibility of a donor solution. The model is validated by a comparison with experimental data, as well as validating the parameter values against each other, using various configurations with donor solution, silicone matrix and skin. Our results show that the model is a good approximation to real multi-layered delivery systems. The model offers the ability of comparing drug release for ibuprofen and diclofenac, which cannot be analysed by the Higuchi model because the dissolution in the latter case turns out to be limited. The experiments and numerical model outlined in this study could also be adjusted to more general formulations, which enhances the utility of the numerical model as a design tool for the development of drug-loaded matrices for trans-membrane and transdermal delivery.

Physical signatures of drug transport through an artificial skin barrier--a proposed model and its validation

The objective of the present study is to better characterize the system acting as a model for the penetration of a pharmaceutical drug into the skin. With a new mathematical formalism, the transport of the drug (dithranol) from a semisolid vaseline suspension into an artificial membrane was described. In our novel approach, we have taken into account not only diffusion but also other effects dependent on chemical reactivity of drug, medium structure, and drug-matrix interactions. The transport equation was solved with two methods: the Laplace transform and the reflection-and-superposition. Despite the applied method, the three-dimensional calculations were found to be, in major parts, in good agreement with the experimental data resulting from the photoacoustic depth profiling. For the first time, from the depth profiling we were capable of estimating not only the diffusion coefficient but also other parameters of the permeation phenomenon as the pore velocity, the first and the zero order reaction rate coefficients.

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.

Application of numerical methods for diffusion-based modeling of skin permeation

The application of numerical methods for mechanistic, diffusion-based modeling of skin permeation is reviewed. Methods considered here are finite difference, method of lines, finite element, finite volume, random walk, cellular automata, and smoothed particle hydrodynamics. First the methods are briefly explained with rudimentary mathematical underpinnings. Current state of the art numerical models are described, and then a chronological overview of published models is provided. Key findings and insights of reviewed models are highlighted. Model results support a primarily transcellular pathway with anisotropic lipid transport. Future endeavors would benefit from a fundamental analysis of drug/vehicle/skin interactions.

Finite and infinite dosing: difficulties in measurements, evaluations and predictions

Due to the increased demand for reliable data regarding penetration into and permeation across human skin, assessment of the absorption of xenobiotics has been gaining in importance steadily. In vitro experiments allow for determining these data faster and more easily than in vivo experiments. However, the experiments described in literature and the subsequent evaluation procedures differ considerably. Here we will give an overview on typical finite and infinite dose experiments performed in fundamental research and on the evaluation of the data. We will point out possible difficulties that may arise and give a short overview on attempts at predicting skin absorption in vitro and in vivo.

Improved input parameters for diffusion models of skin absorption

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

Detailed modeling of skin penetration--an overview

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

Towards drug quantification in human skin with confocal Raman microscopy

Understanding the penetration behaviour of drugs into human skin is a prerequisite for the rational development and evaluation of effective dermal drug delivery. The general procedure for the acquisition of quantitative drug penetration profiles in human skin is performed by sequential segmentation and extraction. Unfortunately, this technique is destructive, laborious and lacks spatial resolution. Confocal Raman microscopy bares the potential of a chemically selective, label free and nondestructive analysis. However, the acquisition of quantitative drug depth profiles within skin by Raman microscopy is impeded by imponderable signal attenuation inside the tissue. In this study, we present a chemical semi-solid matrix system simulating the optical properties of human skin. This system serves as a skin surrogate for investigation of Raman signal attenuation under controlled conditions. Caffeine was homogeneously incorporated within the skin surrogate, and Raman intensity depth profiles were acquired. A mathematical algorithm describing the Raman signal attenuation within the surrogate was derived from these profiles. Human skin samples were incubated with caffeine, and Raman intensity depth profiles were similarly acquired. The surrogate algorithm was successfully applied to correct the drug profiles in human skin for signal attenuation. For the first time, a mathematical algorithm was established, which allows correction of Raman signal attenuation in human skin, thus facilitating reliable drug quantification in human skin by confocal Raman spectroscopy.

Finite dose skin mass balance including the lateral part: comparison between experiment, pharmacokinetic modeling and diffusion models

This work investigates in vitro finite dose skin absorption of the model compounds flufenamic acid and caffeine experimentally and mathematically. The mass balance in different skin compartments (donor, stratum corneum (SC), deeper skin layers (DSL), lateral skin parts and acceptor) is analyzed as a function of time. For both substances high amounts were found in the lateral skin compartment after 6h of incubation, which emphasizes not to elide these parts in the modeling. Here, three different mathematical models were investigated and tested with the experimental data: a pharmacokinetic model (PK), a detailed microscopic two-dimensional diffusion model (MICRO) and a macroscopic homogenized diffusion model (MACRO). While the PK model was fitted to the experimental data, the MICRO and the MACRO models employed input parameters derived from infinite dose studies to predict the underlying diffusion process. All models could satisfyingly predict or describe the experimental data. The PK model and MACRO model also feature the lateral parts.