Mechanistic skin penetration model by the COSMOperm method: Routes of permeation, vehicle effects and skin variations in the healthy and compromised skin

MemCross: Accelerated Weight Histogram method to assess membrane permeability

Passive permeation events across biological membranes are determining steps in the pharmacokinetics of xenobiotics. To reach an accurate and rapid prediction of membrane permeation coefficients of drugs is a complex challenge, which can efficiently support drug discovery. Such predictions are indeed highly valuable as they may guide the selection of potential leads with optimum bioavailabilities prior to synthesis. Theoretical models exist to predict these coefficients. Many of them are based on molecular dynamics (MD) simulations, which allow calculation of permeation coefficients through the evaluation of both the potential of mean force (PMF) and the diffusivity profiles. However, these simulations still require intensive computational efforts, and novel methodologies should be developed and benchmarked. Free energy perturbation (FEP) method was recently shown to estimate PMF with a significantly reduced computational cost compared to the adaptive biasing force method. This benchmarking was achieved with small molecules, namely short-chain alcohols. Here, we show that to estimate the PMF of bulkier, drug-like xenobiotics, conformational sampling is a critical issue. To reach a sufficient sampling with FEP calculations requires a relatively long time-scale, which can lower the benefits related to the computational gain. In the present work, the Accelerated Weight Histogram (AWH) method was employed for the first time in all-atom membrane models. The AWH-based protocol, named MemCross, appears affordable to estimate PMF profiles of a series of drug-like xenobiotics, compared to other enhanced sampling methods. The continuous exploration of the crossing pathway by MemCross also allows modeling subdiffusion by computing fractional diffusivity profiles. The method is also versatile as its input parameters are largely insensitive to the molecule properties. It also ensures a detailed description of the molecule orientations along the permeation pathway, picturing all intermolecular interactions at an atomic resolution. Here, MemCross was applied on a series of 12 xenobiotics, including four weak acids, and a coherent structure-activity relationship was established.

In Silico Prediction of Skin Permeability Using a Two-QSAR Approach

Topical and transdermal drug delivery is an effective, safe, and preferred route of drug administration. As such, skin permeability is one of the critical parameters that should be taken into consideration in the process of drug discovery and development. The ex vivo human skin model is considered as the best surrogate to evaluate in vivo skin permeability. This investigation adopted a novel two-QSAR scheme by collectively incorporating machine learning-based hierarchical support vector regression (HSVR) and classical partial least square (PLS) to predict the skin permeability coefficient and to uncover the intrinsic permeation mechanism, respectively, based on ex vivo excised human skin permeability data compiled from the literature. The derived HSVR model functioned better than PLS as represented by the predictive performance in the training set, test set, and outlier set in addition to various statistical estimations. HSVR also delivered consistent performance upon the application of a mock test, which purposely mimicked the real challenges. PLS, contrarily, uncovered the interpretable relevance between selected descriptors and skin permeability. Thus, the synergy between interpretable PLS and predictive HSVR models can be of great use for facilitating drug discovery and development by predicting skin permeability.

Concentration of Antioxidant Compounds from Calendula officinalis through Sustainable Supercritical Technologies, and Computational Study of Their Permeability in Skin for Cosmetic Use

The growing interest in the cosmetic industry in using compounds of natural and sustainable origin that are safe for humans is encouraging the development of processes that can satisfy these needs. Chlorogenic acid (CHA), caffeic acid (CAF) and ferulic acid (FA) are three compounds widely used within the cosmetic industry due to their functionalities as antioxidants, collagen modifiers or even as radiation protectors. In this work, two advanced separation techniques with supercritical CO2 are used to obtain these three compounds from Calendula officinalis, and these are then evaluated using a computational skin permeability model. This model is encompassed by the COSMO-RS model, the calculations of which make it possible to study the behaviour of the compounds in the epidermis. The results show that both CAF and FA are retained in the stratum corneum, while CHA manages to penetrate to the stratum spinosum. These compounds were concentrated by antisolvent fractionation with super-critical CO2 using a Response Surface Methodology to study the effect of pressure and CO2 flow rate. CHA, CAF and FA were completely retained in the precipitation vessel, with concentrations between 40% and 70% greater than in the original extract. The conditions predicted that the optimal overall yield and enrichment achieved would be 153 bar and 42 g/min.

Modeling the Percutaneous Absorption of Solvent-deposited Solids Over a Wide Dose Range

The transient absorption of two skin care agents, niacinamide (nicotinamide, NA) and methyl nicotinate (MN), solvent-deposited on ex vivo human skin mounted in Franz diffusion cells has been analyzed according to a new variation on a recently published mechanistic skin permeability model (Yu et al. 2020. J Pharm Sci 110:2149-56). The model follows the absorption and evaporation of two components, solute and solvent, and it includes both a follicular transport component and a dissolution rate limitation for high melting, hydrophilic solids deposited on the skin. Explicit algorithms for improving the simulation of transient diffusion of solvent-deposited solids are introduced. The simulations can account for the ex vivo skin permeation time course of both NA and MN over a dose range exceeding 4.5 orders of magnitude. The model allows one to describe on a mechanistic basis why the percutaneous absorption rate of NA is approximately 60-fold lower than that of its lower melting, more lipophilic analog, MN. It furthermore suggests that MN perturbs stratum corneum barrier lipids and increases their permeability while NA does not, presenting a challenge to molecular modelers engaged in simulating biological lipid barriers.

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.

Role of pharmacokinetic consideration for the development of drug delivery systems: A historical overview

Drug delivery system is defined as a system or technology to achieve optimum therapeutic effects of drugs through precise control of their movements in the body. In order to optimize function of drug delivery systems aiming at targeting, their whole-body distribution profiles should be systematically evaluated and analyzed, where pharmacokinetic analysis based on the clearance concepts plays important role. Organ perfusion experiments combined with statistical moment analysis further supply detailed information on drug disposition at organ and cellular levels. Based on general relationship between physicochemical properties and distribution profile, macromolecular prodrugs or polymer conjugates of proteins are rationally designed and further introduction of ligand structure brings cell-specific delivery for them. These approaches are also applicable for particulate carriers such as liposomes and offer various opportunities for biological drugs such as nucleic acid drugs for their delivery. Mechanistic approach for dermal absorption analysis based on physiological skin model offers another opportunity in rational design of drug delivery. Potential of drug delivery technology in future medicines such as cell therapy and nanomaterial platform application is further discussed in relation to pharmacokinetic consideration.

Bio-Functional Textiles: Combining Pharmaceutical Nanocarriers with Fibrous Materials for Innovative Dermatological Therapies

In the field of pharmaceutical technology, significant attention has been paid on exploiting skin as a drug administration route. Considering the structural and chemical complexity of the skin barrier, many research works focused on developing an innovative way to enhance skin drug permeation. In this context, a new class of materials called bio-functional textiles has been developed. Such materials consist of the combination of advanced pharmaceutical carriers with textile materials. Therefore, they own the possibility of providing a wearable platform for continuous and controlled drug release. Notwithstanding the great potential of these materials, their large-scale application still faces some challenges. The present review provides a state-of-the-art perspective on the bio-functional textile technology analyzing the several issues involved. Firstly, the skin physiology, together with the dermatological delivery strategy, is keenly described in order to provide an overview of the problems tackled by bio-functional textiles technology. Secondly, an overview of the main dermatological nanocarriers is provided; thereafter the application of these nanomaterial to textiles is presented. Finally, the bio-functional textile technology is framed in the context of the different dermatological administration strategies; a comparative analysis that also considers how pharmaceutical regulation is conducted.

COSMOplex: self-consistent simulation of self-organizing inhomogeneous systems based on COSMO-RS

During the past 20 years, the efficient combination of quantum chemical calculations with statistical thermodynamics by the COSMO-RS method has become an important alternative to force-field based simulations for the accurate prediction of free energies of molecules in liquid systems. While it was originally restricted to homogeneous liquids, it later has been extended to the prediction of the free energy of molecules in inhomogeneous systems such as micelles, biomembranes, or liquid interfaces, but these calculations were based on external input about the structure of the inhomogeneous system. Here we report the rigorous extension of COSMO-RS to a self-consistent prediction of the structure and the free energies of molecules in self-organizing inhomogeneous systems. This extends the application range to many new areas, such as the prediction of micellar structures and critical micelle concentrations, finite loading effects in micelles and biomembranes, the free energies and structure of liquid interfaces, microemulsions, and many more related topics, which often are of great practical importance.