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

An Investigation into the Transdermal Behavior of Active Ingredients by Combination of Experiments and Multiscale Simulations

Transdermal behavior is a critical aspect of studying delivery systems and evaluating the efficacy of cosmetics. However, existing methods face challenges such as lengthy experiments, high cost, and limited model accuracy. Therefore, developing accurate transdermal models is essential for formulation development and effectiveness assessment. In this study, we developed a multiscale model to describe the transdermal behavior of active ingredients in the stratum corneum. Molecular dynamics simulations were used to construct lipid bilayers and determine the diffusion coefficients of active ingredients in different regions of these bilayers. These diffusion coefficients were integrated into a multilayer lipid pathway model using finite element simulations. The simulation results were in close agreement with our experimental results for three active ingredients (mandelic acid (MAN), nicotinamide (NIC), and pyruvic acid (PYR)), demonstrating the effectiveness of our multiscale model. This research provides valuable insights for advancing transdermal delivery methods.

An Antidehydration Hydrogel Based on Zwitterionic Oligomers for Bioelectronic Interfacing

Hydrogels are ideal interfacing materials for on-skin healthcare devices, yet their susceptibility to dehydration hinders their practical use. While incorporating hygroscopic metal salts can prevent dehydration and maintain ionic conductivity, concerns arise regarding metal toxicity due to the passage of small ions through the skin barrier. Herein, an antidehydration hydrogel enabled by the incorporation of zwitterionic oligomers into its network is reported. This hydrogel exhibits exceptional water retention properties, maintaining ≈88% of its weight at 40% relative humidity, 25 °C for 50 days and about 84% after being heated at 50 °C for 3 h. Crucially, the molecular weight design of the embedded oligomers prevents their penetration into the epidermis, as evidenced by experimental and molecular simulation results. The hydrogel allows stable signal acquisition in electrophysiological monitoring of humans and plants under low-humidity conditions. This research provides a promising strategy for the development of epidermis-safe and biocompatible antidehydration hydrogel interfaces for on-skin devices.

Unraveling the Molecular Mechanisms of Alcohol-Mediated Skin Permeation Enhancement: Insights from Molecular Dynamics Simulations

The application of alcohols as permeation enhancers in pharmaceutical and cosmetic formulations has attracted considerable attention, owing to their skin permeation-enhancing effect. Nonetheless, the elucidation of the fundamental mechanisms underlying the skin permeation-enhancing effect remains elusive. In this study, molecular dynamics (MD) simulations were employed to investigate the effect of 1,2-propanediol (1,2-PDO), 1,2-butanediol (1,2-BDO), and ethanol (EtOH) on the stratum corneum (SC) model membrane. The results showed that the effect of alcohols on the SC model membrane displayed a concentration-dependent nature. The alcohols can interact with SC lipids and exhibit a remarkable ability to selectively extract free fatty acid (FFA) molecules from the SC model membrane and make the SC looser. Meanwhile, 1,2-BDO and EtOH can penetrate into SC lipid bilayers at higher concentrations, leading to the formation of continuous hydrophilic defects in SC. The FFA extraction and the formation of continuous hydrophilic defects induced ceramide (CER) tail chains to become more disordered and fluid and also weakened the hydrogen bonding (H-bonding) network among SC lipids. Both the FFA extraction and the continuous hydrophilic defect formation endowed alcohols with the permeation-enhancing effect. The constrained simulations revealed that the free energy barriers decreased for the permeation of the hydrophilic model molecule (COL) across the SC model membranes containing alcohols, particularly for 1,2-BDO and EtOH. The possible permeation-enhancing mechanisms of alcohols were proposed correspondingly. This work not only provided a deep understanding of the transdermal permeation-enhancing behavior of alcohols at the molecular level but also provided necessary reference information for designing effective transdermal drug delivery systems in applications.

Multiscale modeling of skin mechanical Behavior: Effect of dehydrating agent on Collagen's mechanical properties

The dermis, second layer of human skin, is mainly responsible for mechanical response of the skin. The unique viscoelastic nature of this layer arises from the characteristic hierarchical structure of collagen at various length scales. The effect of topical formulation on skin's mechanical properties of great importance for several personal-care applications. Understanding the transport of an active ingredient across skin layer and its effects on the structure of collagen assembly is crucial for successful design of these applications. In this study, we report a multiscale modelling framework mimicking the skin's mechanical behavior. The framework captures the details from the nanoscale (tropocollagen) to microscale (fibers). At first, atomistic molecular dynamics simulations (MDS) of tropocollagen (TC) molecules of various lengths (∼100 nm) were performed to obtain the molecular modulus of TC. The stress-strain response data obtained from these simulations, were utilized in macroscopic models of fibrils and fibers. The modulus obtained from the mentioned framework was in good agreement with earlier reported experimental data. Further, we have utilized this framework to show the effect of dehydrating agent on skin's mechanical response. The hydration effect is utilized in many anti-ageing strategies to improve the overall mechanical property of skin. We showed that on incorporation of hydrating agent, the collagen structure changes significantly at molecular scale which effects the overall response of the skin at macroscopic scale. The reported multiscale framework can further be explored to gain insights into interlinked properties of collagen at much larger scales without extensive molecular simulations and detailed experiments.

Using molecular simulation to understand the skin barrier

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

In silico predictions of absorption of MDI substances after dermal or inhalation exposures to support a category based read-across assessment

Methylenediphenyl diisocyanate (MDI) substances used polyurethane production can range from their simplest monomeric forms (e.g., 4,4'-MDI) to mixtures of the monomers with various homologues, homopolymer, and prepolymer derivatives. The relative dermal or inhalation absorption of 39 constituents of these substances in human were predicted using the GastroPlus® program. Predicted dermal uptake and absorption of the three MDI monomers from an acetone vehicle was 84-86% and 1.4-1.5%, respectively, with lower uptake and absorption predicted for the higher MW analogs. Lower absorption was predicted from exposures in a more lipophilic vehicle (1-octanol). Modeled inhalation exposures afforded the highest pulmonary absorption for the MDI monomers (38-54%), with 3-27% for the MW range of 381-751, and <0.1% for the remaining, higher MW derivatives. Predicted oral absorption, representing mucociliary transport, ranged from 5 to 10% for the MDI monomers, 10-25% for constituents of MW 381-751, and ≤3% for constituents with MW > 900. These in silico evaluations should be useful in category-based, worst-case, Read-Across assessments for MDI monomers and modified MDI substances for potential systemic effects. Predictions of appreciable mucociliary transport may also be useful to address data gaps in oral toxicity testing for this category of compounds.

Does skin permeation kinetics influence efficacy of topical dermal drug delivery system?: Assessment, prediction, utilization, and integration of chitosan biomacromolecule for augmenting topical dermal drug delivery in skin

Skin permeation is an integral part of penetration of topical therapeutics. Zero order in addition to Higuchi permeation kinetic is usually preferred in topical drug delivery cargo. Penetration of therapeutic entities through epidermal barrier is a major challenge for scientific fraternity. Furthermore, penetration of therapeutic entities determines the transportation and ultimately therapeutic efficacy of topical dermal dosage forms. Apart from experimentation models, mathematical equations, in silico docking, molecular dynamics (MDs), and artificial neural network (Neural) techniques are being used to assess free energies and prediction of electrostatic attractions in order to predict the permeation phenomena of therapeutic entities. Therefore, in the present review, we have summarized the significance of kinetic equations, in silico docking, MDs, and ANN in assessing and predicting the penetration behavior of topical therapeutics through dermal dosage form. In addition, the role of chitosan biomacromolecule in modulating permeation of topical therapeutics in skin has also been illustrated using computational techniques.

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.

Integrated “Generate, Make, and Test” for Formulated Products using Knowledge Graphs

In the multi-billion dollar formulated product industry, state of the art continues to rely heavily on experts during the “generate, make and test” steps of formulation design. We propose automation aids to each step with a knowledge graph of relevant information as the central artifact. The generate step usually focuses on coming up with new recipes for intended formulation. We propose to aid the experts who generally carry out this step manually by providing a recommendation system and a templating system on top of the knowledge graph. Using the former, the expert can create a recipe from scratch using historical formulations and related data. With the latter, the expert starts with a recipe template created by our system and substitutes the requisite constituents to form a recipe. In the current state of practice, the three steps mentioned above operate in a fragmented manner wherein observations from one step do not aid other steps in a streamlined manner. Instead of manually operated labs for the make and test steps, we assume automated or robotic labs and in-silico testing, respectively. Using two formulations, namely face cream and an exterior coating, we show how the knowledge graph may help integrate and streamline the communication between the generate, the make, and the test steps. Our initial exploration shows considerable promise.

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.

Not only in silico drug discovery: Molecular modeling towards in silico drug delivery formulations

The use of methods at molecular scale for the discovery of new potential active ligands, as well as previously unknown binding sites for target proteins, is now an established reality. Literature offers many successful stories of active compounds developed starting from insights obtained in silico and approved by Food and Drug Administration (FDA). One of the most famous examples is raltegravir, a HIV integrase inhibitor, which was developed after the discovery of a previously unknown transient binding area thanks to molecular dynamics simulations. Molecular simulations have the potential to also improve the design and engineering of drug delivery devices, which are still largely based on fundamental conservation equations. Although they can highlight the dominant release mechanism and quantitatively link the release rate to design parameters (size, drug loading, et cetera), their spatial resolution does not allow to fully capture how phenomena at molecular scale influence system behavior. In this scenario, the "computational microscope" offered by simulations at atomic scale can shed light on the impact of molecular interactions on crucial parameters such as release rate and the response of the drug delivery device to external stimuli, providing insights that are difficult or impossible to obtain experimentally. Moreover, the new paradigm brought by nanomedicine further underlined the importance of such computational microscope to study the interactions between nanoparticles and biological components with an unprecedented level of detail. Such knowledge is a fundamental pillar to perform device engineering and to achieve efficient and safe formulations. After a brief theoretical background, this review aims at discussing the potential of molecular simulations for the rational design of drug delivery systems.

The Implementation of the Three Rs in Regulatory Toxicity and Biosafety Assessment: The Indian Perspective

Animal models have long served as a basis for scientific experimentation, biomedical research, drug development and testing, disease modelling and toxicity studies, as they are widely thought to provide meaningful, human-relevant predictions. However, many of these systems are resource intensive and time-consuming, have low predictive value and are associated with great social and ethical dilemmas. Often drugs appear to be effective and safe in these classical animal models, but later prove to be ineffective and/or unsafe in clinical trials. These issues have paved the way for a paradigm shift from the use of in vivo approaches, toward the 'science of alternatives'. This has fuelled several research and regulatory initiatives, including the ban on the testing of cosmetics on animals. The new paradigm has been shifted toward increasing the relevance of the models for human predictivity and translational efficacy, and this has resulted in the recent development of many new methodologies, from 3-D bio-organoids to bioengineered 'human-on-a-chip' models. These improvements have the potential to significantly advance medical research globally. This paper offers a stance on the existing strategies and practices that utilise alternatives to animals, and outlines progress on the incorporation of these models into basic and applied research and education, specifically in India. It also seeks to provide a strategic roadmap to streamline the future directions for the country's policy changes and investments. This strategic roadmap could be a useful resource to guide research institutions, industries, regulatory agencies, contract research organisations and other stakeholders in transitioning toward modern approaches to safety and risk assessment that could replace or reduce the use of animals without compromising the safety of humans or the environment.

Challenges and opportunities for small volumes delivery into the skin

Each individual's skin has its own features, such as strength, elasticity, or permeability to drugs, which limits the effectiveness of one-size-fits-all approaches typically found in medical treatments. Therefore, understanding the transport mechanisms of substances across the skin is instrumental for the development of novel minimal invasive transdermal therapies. However, the large difference between transport timescales and length scales of disparate molecules needed for medical therapies makes it difficult to address fundamental questions. Thus, this lack of fundamental knowledge has limited the efficacy of bioengineering equipment and medical treatments. In this article, we provide an overview of the most important microfluidics-related transport phenomena through the skin and versatile tools to study them. Moreover, we provide a summary of challenges and opportunities faced by advanced transdermal delivery methods, such as needle-free jet injectors, microneedles, and tattooing, which could pave the way to the implementation of better therapies and new methods.

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.

Multiscale Modeling of Skin Electroporation

Human skin, the largest external organ of the body, provides a selective barrier to therapeutics applied topically. The molecules having specific chemical and physical properties can only penetrate the deeper layer of the skin. However, the lag time for reaching a steady state in the deeper layer is generally of the order of hours. In order to deliver higher-molecular-weight, charged, and hydrophilic therapeutics in the deeper layer, the skin barrier must be breached. Electroporation is one of the methods used to breach the skin barrier for enhancement of drug permeation and reduction of lag time. However, the underlying mechanism responsible for the enhancement of drug permeation is not well understood. In this study, a multiscale model of skin electroporation is developed by connecting molecular phenomena to a macroscopic model. At the atomic scale, molecular dynamics simulations of the lipid matrix of the human stratum corneum (SC) were performed under the influence of an external electric field. The pores get formed during the electroporation process and the transport properties (diffusivity) of drug molecules are computed. The diffusion coefficient obtained during electroporation was found to be higher than passive diffusion. However, this alone could not explain the multifold increase in the drug flux on application of an electric field as observed in the experiments. Hence, a finite element method (FEM) model of the skin SC is also developed. The release of fentanyl through this model is compared with the available experimental results. Both experimental and simulated results of pore formation on application of an electric field and many folds' increase in drug flux are comparable. Once validated, the framework was used for the design of skin electroporation experiments (in silico) by changing the electric pulse parameters such as voltage, pulse duration, and number of pulses. This multiscale modeling framework provides valuable insight at the molecular and macroscopic levels to design the electroporation experiments. The framework can be utilized as a design tool for skin electroporation applications.

Molecular mechanism of the skin permeation enhancing effect of ethanol: a molecular dynamics study

Ethanol is widely used in various pharmaceutical and cosmetic formulations in order to enhance skin penetration of active ingredients. While it is well known that ethanol partitions into the skin and enhances the permeation of both polar and nonpolar molecules, the exact mechanisms by which it enhances skin permeability are not fully understood. Several mechanisms have been proposed including lipid extraction from the stratum corneum (SC), fluidisation of SC lipid bilayer, alteration of SC protein conformation and enhancement of the drug solubility in the SC lipids. In this study, we performed molecular dynamics (MD) simulations of SC lipid bilayers comprised of an equimolar mixture of ceramides, cholesterol and free fatty acid in the presence of aqueous mixtures of ethanol. Various unrestrained MD simulations were performed in the presence of aqueous ethanol solution at molar ratios (x) ranging from x = 0 to x = 1. It was found that ethanol enhances bilayer permeability by dual actions (a) extraction of the skin lipids and (b) enhancing the mobility of lipid chains. Ethanol's permeation enhancing effect arises from its superior ability to form hydrogen bonds with headgroup atoms of skin lipids. Further, the free energy of extraction of ceramides (CER) and fatty acids (FFA) from the lipid bilayer was studied using umbrella sampling simulations. The free energy of extraction of CER was found to be much higher compared to FFA for all ethanol concentrations which shows that CER are difficult to extract as compared to FFA. Finally, the permeation of benzoic acid drug molecules through the skin lipid bilayer is shown in presence of ethanol molecules. It was found that ethanol selectively targets the FFA of the skin lipid bilayer and extracts it out of the lipid bilayer within few microseconds. Further, ethanol penetrates inside the lipid layer and creates the channels from which drug molecules can easily cross the lipid layer. Our observations (both in unrestrained and umbrella sampling simulations) are consistent with the experimental findings reported in the literature. The simulation methodology could be used for design and testing of permeation enhancers (acting on skin SC lipid lamella) for topical and transdermal drug delivery applications.

Concentration dependent action of mechanism of ethanol on skin SC lipid barrier.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

Molecular dynamic simulations of oxidized skin lipid bilayer and permeability of reactive oxygen species

Lipid peroxidation by reactive oxygen species (ROS) during oxidative stress is non-enzymatic damage that affects the integrity of biological membrane, and alters the fluidity and permeability. We conducted molecular dynamic simulation studies to evaluate the structural properties of the bilayer after lipid peroxidation and to measure the permeability of distinct ROS. The oxidized membrane contains free fatty acid, ceramide, cholesterol, and 5α-hydroperoxycholesterol (5α-CH). The result of unconstrained molecular dynamic simulations revealed that lipid peroxidation causes area-per-lipid of the bilayer to increase and bilayer thickness to decrease. The simulations also revealed that the oxidized group of 5α-CH (-OOH) moves towards the aqueous layer and its backbone tilts causing lateral expansion of the bilayer membrane. These changes are detrimental to structural and functional properties of the membrane. The measured free energy profile for different ROS (H₂O₂, HO₂, HO, and O₂) across the peroxidized lipid bilayer showed that the increase in lipid peroxidation resulted in breaching barrier decrease for all species, allowing easy traversal of the membrane. Thus, lipid peroxidation perturbs the membrane barrier and imposes oxidative stress resulting into apoptosis. The collective insights increase the understanding of oxidation stress at the atomic level.

Effect of Chemical Permeation Enhancers on Skin Permeability: In silico screening using Molecular Dynamics simulations

Breaching of the skin barrier is essential for delivering active pharmaceutical ingredients (APIs) for pharmaceutical, dermatological and aesthetic applications. Chemical permeation enhancers (CPEs) are molecules that interact with the constituents of skin’s outermost and rate limiting layer stratum corneum (SC), and increase its permeability. Designing and testing of new CPEs is a resource intensive task, thus limiting the rate of discovery of new CPEs. In-silico screening of CPEs in a rigorous skin model could speed up the design of CPEs. In this study, we performed coarse grained (CG) molecule dynamics (MD) simulations of a multilayer skin lipid matrix in the presence of CPEs. The CPEs are chosen from different chemical functionalities including fatty acids, esters, and alcohols. A multi-layer in-silico skin model was developed. The CG parameters of permeation enhancers were also developed. Interactions of CPEs with SC lipids was studied in silico at three different CPE concentrations namely, 1% w/v, 3% w/v and 5% w/v. The partitioning and diffusion coefficients of CPEs in the SC lipids were found to be highly size- and structure-dependent and these dependencies are explained in terms of structural properties such as radial distribution function, area per lipid and order parameter. Finally, experimentally reported effects of CPEs on skin from the literature are compared with the simulation results. The trends obtained using simulations are in good agreement with the experimental measurements. The studies presented here validate the utility of in-silico models for designing, screening and testing of novel and effective CPEs.

Predicting Skin Permeability by Means of Computational Approaches: Reliability and Caveats in Pharmaceutical Studies

The skin is the main barrier between the internal body environment and the external one. The characteristics of this barrier and its properties are able to modify and affect drug delivery and chemical toxicity parameters. Therefore, it is not surprising that permeability of many different compounds has been measured through several in vitro and in vivo techniques. Moreover, many different in silico approaches have been used to identify the correlation between the structure of the permeants and their permeability, to reproduce the skin behavior, and to predict the ability of specific chemicals to permeate this barrier. A significant number of issues, like interlaboratory variability, experimental conditions, data set building rationales, and skin site of origin and hydration, still prevent us from obtaining a definitive predictive skin permeability model. This review wants to show the main advances and the principal approaches in computational methods used to predict this property, to enlighten the main issues that have arisen, and to address the challenges to develop in future research.

Insight from Molecular dynamic simulation of reactive oxygen species in oxidized skin membrane

Non-enzymatic lipid peroxidation of the skin-lipid bilayer causes perturbations that affect the biomembrane structure, function, and permeability of reactive oxygen species (ROS). In the present study, we employed molecular dynamics simulations to study the effect of lipid peroxidation on the bilayer structural properties and permeability of various ROS. The oxidized skin-lipid bilayer was composed of ceramide, cholesterol, free fatty acid, and 5α-hydroperoxycholesterol (5α-CH). The simulation showed that, upon oxidation, the oxidized group (-OOH) of 5α-CH migrates towards the aqueous phase and the backbone of 5α-CH tilts, which causes the membrane to expand laterally. Measurements of the permeability of all ROS along the oxidized skin-lipid bilayer revealed a decreased breaching barrier for all the species as the degree of peroxidation increased, with a resulting easy passage across the membrane. The insights from the simulations indicate that lipid peroxidation might perturb the membrane barrier, thereby inflicting oxidative stress that leads to apoptosis. This study helps to understand oxidative stress at the atomic level. To our knowledge, this is the first reported molecular dynamics simulation study on oxidized skin-lipid bilayer and permeability of ROS.

Insight Into the Molecular Dynamic Simulation Studies of Reactive Oxygen Species in Native Skin Membrane

In recent years, the role of reactive oxygen species (ROS) in regulating cancer cell apoptosis, inflammation, cell ischemia, and cell signaling pathways has been well established. The most common sources of intracellular ROS are the mitochondrial electron transport system, NADH oxidase, and cytochrome P450. In this study, we investigated the dynamics and permeability of ROS using molecular dynamics (MD) simulations on native skin-lipid bilayer membranes. Native skin-lipid bilayers are composed of ceramide, cholesterol, and free fatty acid in an almost equal molar ratio (1:1:1). Dynamic distribution studies on ROS, i.e., hydrogen peroxide (H₂O₂) and O₂ (¹O₂ by analogy), revealed that these species interact with cholesterol as a primary target in lipid peroxidation of the skin-lipid bilayer. Moreover, the permeability of ROS, i.e., H₂O₂, hydroxyl radicals (HO), hydroperoxy radical (HOO), and O₂, along the skin-lipid bilayer was measured using free energy profiles (FEPs). The FEPs showed that in spite of high-energy barriers, ROS traveled through the membrane easily. Breaching the free energy barriers, these ROS permeated into the membrane, inflicting oxidative stress, and causing apoptosis. Collectively, the insight acquired from simulations may result in a better understanding of oxidative stress at the atomic level.