Molecular Structure of the Long Periodicity Phase in the Stratum Corneum

Simulating the Skin Permeation Process of Ionizable Molecules

It is commonly assumed that ionizable molecules, such as drugs, permeate through the skin barrier in their neutral form. By using molecular dynamics simulations of the charged and neutral states separately, we can study the dynamic protonation behavior during the permeation process. We have studied three weak acids and three weak bases and conclude that the acids are ionized to a larger extent than the bases, when passing through the headgroup region of the lipid barrier structure, at pH values close to their pKa. It can also be observed that even if these dynamic protonation simulations are informative, in the cases studied herein they are not necessary for the calculation of permeability coefficients. It is sufficient to base the calculations only on the neutral form, as is commonly done.

Investigating the Free Volumes as Nanospaces in Human Stratum Corneum Lipid Bilayers Using Positron Annihilation Lifetime Spectroscopy (PALS)

This work is the first one that provides not only evidence for the existence of free volumes in the human stratum corneum but also focuses on comparing these experimental data, obtained through the unique positron annihilation lifetime spectroscopy (PALS) method, with theoretical values published in earlier works. The mean free volume of 0.269 nm was slightly lower than the theoretical value of 0.4 nm. The lifetime τ3 (1.83 ns with a coefficient of variation CV of 3.21%) is dependent on the size of open sites in the skin. This information was used to calculate the free volume radius R (0.269 nm with CV 2.14%), free volume size Vf (0.081 nm3 with CV 4.69%), and the intensity I3 (9.01% with CV 10.94%) to estimate the relative fractional free volume fv (1.32 a.u. with CV 13.68%) in human skin ex vivo. The relation between the lifetime of o-Ps (τ3) and the radius of free volume (R) was formulated using the Tao-Eldrup model, which assumes spherical voids and applies to sites with radii smaller than 1 nm. The results indicate that PALS is a powerful tool for confirming the existence of free volumes and determining their size. The studies also focused on describing the probable locations of these nanospaces in SC lipid bilayers. According to the theory, these play an essential role in dynamic processes in biological systems, including the diffusion of low-molecular-weight hydrophobic and moderately hydrophilic molecules. The mechanism of their formation has been determined by the molecular dynamics of the lipid chains.

Molecular Dynamics Investigation into CerENP's Effect on the Lipid Matrix of Stratum Corneum

The extracellular lipid matrix in the stratum corneum (SC) plays a critical role in skin barrier functionality, comprising three primary components: ceramides, cholesterol, and free fatty acids. The diverse ceramides, differentiated by molecular structures such as hydroxylations and varying chain lengths, are essential for the lipid matrix's structural integrity. Recently, a new subclass of ceramide, 1-O-acylceramide NP (CerENP), has been identified; however, its precise role in the lipid matrix of the SC is still elusive. Herein, we investigate the role of CerENP on the structure and permeability of the SC using molecular dynamics simulations. Our findings indicate that CerENP contributes to a compact lipid matrix in the lateral dimension of our SC model with a repeat distance of about 13 nm. Additionally, ethanol permeability assessments show that CerENP effectively reduces molecular penetration through the lipid matrix. This study provides an insight into the role of a new subclass of ceramide in the SC, enhancing our understanding of skin structure and the mechanisms behind barrier dysfunction in skin diseases.

Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study

The construction of the stratum corneum (SC) is crucial to the problems of transdermal drug delivery. SC consists of the keratinocyte layers and the lipid matrix surrounding it. Among them, the lipid matrix is the barrier for many exogenous molecules, mainly composed of ceramides (CERs), free fatty acids (FFA), and cholesterol (CHOL). In this work, we developed single-component (CERs, CER-NS, and CER-EOS) and six three-component models, and each model was simulated by using the GROMOS-54A7 force field. Short-period phase (SPP) and long-period phase (LPP) systems were established separately, and area per lipid (APL), thickness, order of carbon chain (SCD), and density distribution were analyzed. The transition of CER-NS and CER-EOS in LPP was observed. The results of hydrogen bonds in the lipid systems indicated that a strong hydrogen-bond network was formed between the skin-lipid bilayers. Umbrella sampling method simulations were performed to calculate the free energy change of ethanol moving into the skin-lipid bilayer. The results revealed that ethanol molecules pulled some water molecules into the membrane when they passed through SPP-1. Our findings provided some insights and models of the stratum corneum that could be used for the subsequent mechanism of macromolecule permeation through membranes in drugs, cosmetics, and so on.

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.

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.

Formulation Development and Molecular Mechanism Characterization of Long-Acting Patches of Asenapine for Efficient Delivery by Combining API-ILs Strategy and Controlled-Release Polymers

The objective of our study, which combined API-ILs strategy and controlled-release polymers, was to prepare a 72 h long-acting drug-in-adhesive patch for optimum delivery of asenapine (ASE). Special attention was paid to the permeation promotion mechanism and the controlled release behavior of ASE-ILs in pressure sensitive adhesives (PSA). Formulation factors were investigated by ex vivo transdermal experiments. The optimized patch was evaluated by pharmacokinetics study and skin irritation test. The obtained formulation was as follows, 15% w/w ASE-MA (about 1136 μg/cm2 ASE, 413 μg/cm2 MA), AACONH2 (Amide adhesive) as the matrix, 80 μm thickness, backing film of CoTran™ 9733. The optimized patch displayed satisfactory ex vivo and in vivo performance with Q 72 h of 620 ± 44 µg/cm2 and Fabs of 62.4%, which utilization rate (54.6%) was significantly higher than the control group (38.3%). By using the classical shake flask method, 13C NMR, DSC, and FTIR, the physicochemical properties and structure of ILs were characterized. log Do/w, ATR-FTIR, Raman, and molecular dynamics simulation results confirmed that ASE-MA (MA: 3-Methoxypropionic acid) had appropriate lipophilicity, and affected lipid fluidity as well as the conformation of keratin to improve the skin permeation. The FTIR, MDSC, rheology, and molecular docking results revealed that hydrogen bond (H-bond), were formed between ASE-MA and PSA, and the drug increased the molecular mobility of polymer chains. In summary, the 72 h long-acting patch of ASE was successfully prepared and it supplied a reference for the design of long-acting patches with ASE.

The effects of size and surface functionalization of polystyrene nanoplastics on stratum corneum model membranes: An experimental and computational study

Nanoplastics are mainly generated from the decomposition of plastic waste and artificial production and have attracted much attention due to their wide distribution in the environment and the potential risk for humans. As the largest organ of the human body, the skin is inevitably in contact with nanoplastics. Stratum corneum is the first barrier when the skin is exposed to nanoplastics. However, little is known about the interactions between nanoplastics and stratum corneum. Here, the effects of particle size and surface functionalization (amino-modified and carboxy-modified) of polystyrene nanoplastics on the stratum corneum models were studied by Langmuir monolayer and molecular dynamics simulations. An equimolar mixture of ceramide/cholesterol/free fatty acid was used to mimic stratum corneum intercellular lipids. The Langmuir monolayer studies demonstrated that the larger size and surface functionalization of polystyrene nanoplastics significantly reduced the stability of stratum corneum lipid monolayer in a concentration-dependent fashion. Simulation results elucidated that functionalized polystyrene oligomers had a stronger interaction with lipid components of the stratum corneum model membrane. The cell experiments also indicated that functionalized polystyrene nanoplastics, especially for amino-modified polystyrene nanoplastics, had significant cytotoxicity on normal human dermal fibroblast cells. Our results provide fundamental information and the basis for a deeper understanding of the health risks of nanoplastics to humans.

The Intriguing Molecular Dynamics of Cer[EOS] in Rigid Skin Barrier Lipid Layers Requires Improvement of the Model

Omega-O-acyl ceramides such as 32-linoleoyloxydotriacontanoyl sphingosine (Cer[EOS]) are essential components of the lipid skin barrier, which protects our body from excessive water loss and the penetration of unwanted substances. These ceramides drive the lipid assembly to epidermal-specific long periodicity phase (LPP), structurally much different than conventional lipid bilayers. Here, we synthesized Cer[EOS] with selectively deuterated segments of the ultralong N-acyl chain or deuterated or ¹³C-labeled linoleic acid and studied their molecular behavior in a skin lipid model. Solid-state ²H NMR data revealed surprising molecular dynamics for the ultralong N-acyl chain of Cer[EOS] with increased isotropic motion towards the isotropic ester-bound linoleate. The sphingosine moiety of Cer[EOS] is also highly mobile at skin temperature, in stark contrast to the other LPP components, N-lignoceroyl sphingosine acyl, lignoceric acid and cholesterol, which are predominantly rigid. The dynamics of the linoleic chain is quantitatively described by distributions of correlation times and using dynamic detector analysis. These NMR results along with neutron diffraction data suggest an LPP structure with alternating fluid (sphingosine chain-rich), rigid (acyl chain-rich), isotropic (linoleate-rich), rigid (acyl-chain rich), and fluid layers (sphingosine chain-rich). Such an arrangement of the skin barrier lipids with rigid layers separated with two different dynamic "fillings" i) agrees well with ultrastructural data, ii) satisfies the need for simultaneous rigidity (to ensure low permeability) and fluidity (to ensure elasticity, accommodate enzymes or antimicrobial peptides), and iii) offers a straightforward way to remodel the lamellar body lipids into the final lipid barrier.

In Silico Prediction of Stratum Corneum Partition Coefficients via COSMOmic and Molecular Dynamics Simulations

Stratum corneum (SC) is the main barrier of human skin where the inter-corneocytes lipids provide the main pathway for transdermal permeation of functional actives of skin care and health. Molecular dynamics (MD) has been increasingly used to simulate the SC lipid bilayer structure so that the barrier property and its affecting factors can be elucidated. Among reported MD simulation studies, solute partition in the SC lipids, an important parameter affecting SC permeability, has received limited attention. In this work, we combine MD simulation with COSMOmic to predict the partition coefficients of dermatologically relevant solutes in SC lipid bilayer. Firstly, we run MD simulations to obtain equilibrated SC lipid bilayers with different lipid types, compositions, and structures. Then, the simulated SC lipid bilayer structures are fed to COSMOmic to calculate the partition coefficients of the solutes. The results show that lipid types and bilayer geometries play a minor role in the predicted partition coefficients. For the more lipophilic solutes, the predicted results of solute partition in SC lipid bilayers agree well with reported experimental values of solute partition in extracted SC lipids. For the more hydrophilic molecules, there is a systematical underprediction. Nevertheless, the MD/COSMOmic approach correctly reproduces the phenomenological correlation between the SC lipid/water partition coefficients and the octanol/water partition coefficients. Overall, the results show that the MD/COSMOmic approach is a fast and valid method for predicting solute partitioning into SC lipids and hence supporting the assessment of percutaneous absorption of skin care ingredients, dermatological drugs as well as environmental pollutants.

Biomimetic Stratum Corneum Liposome Models: Lamellar Organization and Permeability Studies

The stratum corneum (SC), the outer layer of the skin, plays a crucial role as a barrier protecting the underlying cells from external stress. The SC comprises three key components: ceramide (CER), free fatty acid (FFA), and cholesterol, along with small fractions of cholesterol sulfate and cholesterol ester. In order to gain a deeper understanding about the interdependence of the two major components, CER and FFA, on the organizational, structural, and functional properties of the SC layer, a library of SC lipid liposome (SCLL) models was developed by mixing CER (phytosphingosine or sphingosine), FFA (oleic acid, palmitic acid, or stearic acid), cholesterol, and cholesterol sulfate. Self-assembly of the SC lipids into lamellar phases was first confirmed by small-angle X-ray scattering. Short periodicity and long periodicity phases were identified for SCLLs containing phytosphingosines and sphingosine CERs, respectively. Furthermore, unsaturation in the CER acyl and FFA chains reduced the lipid conformational ordering and packing density of the liposomal bilayer, which were measured by differential scanning calorimetry and Fourier transform infrared spectroscopy. The introduction of unsaturation in the CER and/or FFA chains also impacted the lamellar integrity and permeability. This extensive library of SCLL models exhibiting physiologically relevant lamellar phases with defined structural and functional properties may potentially be used as a model system for screening pharmaceuticals or cosmetic agents.

Skin permeability prediction with MD simulation sampling spatial and alchemical reaction coordinates

A molecular-level understanding of skin permeation may rationalize and streamline product development, and improve quality and control, of transdermal and topical drug delivery systems. It may also facilitate toxicity and safety assessment of cosmetics and skin care products. Here, we present new molecular dynamics simulation approaches that make it possible to efficiently sample the free energy and local diffusion coefficient across the skin's barrier structure to predict skin permeability and the effects of chemical penetration enhancers. In particular, we introduce a new approach to use two-dimensional reaction coordinates in the accelerated weight histogram method, where we combine sampling along spatial coordinates with an alchemical perturbation virtual coordinate. We present predicted properties for 20 permeants, and demonstrate how our approach improves correlation with ex vivo/in vitro skin permeation data. For the compounds included in this study, the obtained log KPexp-calc mean square difference was 0.9 cm2 h-2.

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.

Implications of surfactant hydrophobic chain architecture on the Surfactant-Skin lipid model interaction

Although surfactants have been widely used in skin care and other related applications, our knowledge about how surfactants interact with stratum corneum (SC) lipids remains limited. This work reports how surfactants interact with a lipid SC model by neutron diffraction and molecular dynamics (MD) simulations, focusing on examining the impact of surfactant molecular architecture. The surfactant-SC mixed membrane was constructed by an equimolar mixture of ceramide/cholesterol/fatty acids and surfactant at 1% molar ratio of total lipids. The arrangements of water and surfactant molecules in the membrane were obtained through neutron scattering length density (NSLD) profiles via contrast variation method, meanwhile, MD simulation clearly demonstrated the mechanism of hydration change in the surfactant-model SC mixed membrane. No drastic difference was detected in the repeating distance of the short periodicity phase (SPP) upon adding surfactants, however, it significantly enhanced the membrane hydration and reduced the amount of phase separated crystalline cholesterol, showing a strong dependence on surfactant chain length, branching and double bond. This work clearly demonstrates how surfactant architecture affects its interaction with the SC membrane, providing useful guidance for either choosing an existing surfactant or designing a new one for surfactant-based transdermal application.

Cholesterol Sulfate Fluidizes the Sterol Fraction of the Stratum Corneum Lipid Phase and Increases its Permeability

Desulfation of cholesterol sulfate (CholS) to cholesterol (Chol) is an important event in epidermal homeostasis and necessary for stratum corneum (SC) barrier function. The CholS/Chol ratio decreases during SC maturation but remains high in pathological conditions, such as X-linked ichthyosis, characterized by dry and scaly skin. The aim of this study was to characterize the influence of the CholS/Chol molar ratio on the structure, dynamics, and permeability of SC lipid model mixtures. We synthesized deuterated CholS and investigated lipid models with specifically deuterated components using ²H solid-state NMR spectroscopy at temperatures from 25°C to 80°C. Although the rigid acyl chains in ceramides and fatty acids remained essentially rigid upon variation of the CholS/Chol ratio, both sterols were increasingly fluidized in lipid models containing higher CholS concentrations. We also show the X-ray repeat distance of the lipid lamellar phase (105 Å) and the orthorhombic chain packing of the ceramide’s acyl chains and long free fatty acids did not change upon the variation of the CholS content. However, the Chol phase separation visible in models with high Chol concentration disappeared at the 50:50 CholS/Chol ratio. This increased fluidity resulted in higher permeabilities to model markers of these SC models. These results reveal that a high CholS/Chol ratio fluidizes the sterol fraction and increases the permeability of the SC lipid phase while maintaining the lamellar lipid arrangement with an asymmetric sterol distribution.

All-Atom Modeling of Complex Cellular Membranes

Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.

A Systematic Quantitative Evaluation of Permeation Enhancement Window: Transdermal Permeation Enhancing Dynamics Establishment and Molecular Mechanisms Characterization of Permeation Enhancer

At present, transdermal permeation enhancing dynamics studies on permeation enhancers are still limited. In this study, these dynamics were established based on the content of enhancer Plurol Oleique CC in skin (C_POCC) and the increment of drug permeation amount (ΔQ). A new concept deemed "permeation enhancement window" (ΔC_POCC), comprised of a threshold dose (C^thr), maximal dose (C^max) and permeation enhancement efficiency (Eff) was used to evaluate the enhancement effect of POCC for different drugs. According to results of FT-IR, ATR-FTIR and DSC analyses, the higher CPOCC of patches containing acidic drugs vs. basic drugs resulted from their stronger interaction with pressure-sensitive adhesives, leading to more free POCC and a greater disturbing effect on stratum corneum (SC) lipids. Below C^thr, a longer lag phase for acidic drugs resulted from more POCC required to compete with ceramide. When CPOCC exceeded C^max by about 400 μg/g, plateau phases for all drugs were reached due to the upper limit of SC lipid fluidity, as confirmed by SAXS and Raman imaging. In summary, the differences in the permeation enhancement window for the test drugs resulted from the varied interaction strengths among POCC, drugs and adhesives, as well as changeable SC lipid fluidity.

Evaluation of Constrained and Restrained Molecular Dynamics Simulation Methods for Predicting Skin Lipid Permeability

Recently, molecular dynamics (MD) simulations have been utilized to investigate the barrier properties of human skin stratum corneum (SC) lipid bilayers. Different MD methods and force fields have been utilized, with predicted permeabilities varying by few orders of magnitude. In this work, we compare constrained MD simulations with restrained MD simulations to obtain the potential of the mean force and the diffusion coefficient profile for the case of a water molecule permeating across an SC lipid bilayer. Corresponding permeabilities of the simulated lipid bilayer are calculated via the inhomogeneous solubility diffusion model. Results show that both methods perform similarly, but restrained MD simulations have proven to be the more robust approach for predicting the potential of the mean force profile. Critical to both methods are the sampling of the whole trans-bilayer axis and the following symmetrization process. Re-analysis of the previously reported free energy profiles showed that some of the discrepancies in the reported permeability values is due to misquotation of units, while some are due to the inaccurately obtained potential of the mean force. By using the existing microscopic geometrical models via the intercellular lipid pathway, the permeation through the whole SC is predicted from the MD simulation results, and the predicted barrier properties have been compared to experimental data from the literature with good agreement.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

An investigation on the effect of drug physicochemical properties on the enhancement strength of enhancer: The role of drug-skin-enhancer interactions

The aim of present work was to investigate the influence of drug physicochemical properties on the enhancement effect of enhancers, which guided the application of enhancers in different drug transdermal prescriptions. Firstly, Polyglyceryl-3 dioleate (POCC) was selected as a model enhancer and its enhancement effect on ten drugs was assessed by in vitro skin permeation experiment. Secondly, the correlation analysis of physicochemical properties of drugs was carried out from the aspects of partition and permeation. The interactions of drug-skin-POCC were elucidated by FT-IR, molecular docking, solubility parameters calculation, ATR-FTIR, Raman study, molecular dynamics simulation and confocal laser scanning microscopy (CLSM). The results showed that the enhancement ratio (ER) of drugs was ranging from 2.23 to 7.45. On one hand, the miscibility between drugs with low polar surface area (P.S.A) and donor solution was decreased more pronounced by the addition of POCC because of the drug was difficult to form hydrogen-bond with POCC, facilitating the vehicle/SC partition of drugs. On the other hand, the permeation of drugs with low P.S.A and polarizability was enhanced more significantly by POCC because the drug was less likely to interact with skin lipids compared to others, causing that POCC had more chance to interact with skin lipids to improve permeation drugs across the SC more easily. In conclusion, the different strength of drug-skin-POCC interactions was the main reason for the discrepancy in the enhancement effect of the POCC on ten drugs, which laid a basis for the research on the drug-specific molecular mechanisms of enhancers.

How to control interactions of cellulose-based biomaterials with skin: the role of acidity in the contact area

Being able to control the interactions of biomaterials with living tissues and skin is highly desirable for many biomedical applications. This is particularly the case for cellulose-based materials which provide one of the most versatile platforms for tissue engineering due to their strength, biocompatibility and abundance. Achieving such control, however, requires detailed molecular-level knowledge of the dominant interaction mechanisms. Here, we employed both biased and unbiased atomic-scale molecular dynamics simulations to explore how cellulose crystals interact with model stratum corneum bilayers, ternary mixtures of ceramides, cholesterol, and free fatty acids. Our findings show that acidity in the contact area directly affects binding between cellulose and the stratum corneum bilayer: Protonation of free fatty acids in the bilayer promotes attractive cellulose-bilayer interactions. We identified two major factors that control the cellulose-skin interactions: (i) the electrostatic repulsion between a cellulose crystal and the charged (anionic due to deprotonated fatty acids) surface of a stratum corneum bilayer and (ii) the cellulose-stratum corneum hydrogen bonding. When less than half of the fatty acids in the bilayer are protonated, the first factor dominates and there is no binding to skin. At a larger degree of fatty acid protonation the cellulose-stratum corneum hydrogen bonding prevails yielding a tight binding. Remarkably, we found that ceramide molecules are the key component in hydrogen bonding with cellulose. Overall, our findings highlight the critical role of fatty acid protonation in biomaterial-stratum corneum interactions and can be used for optimizing the surface properties of cellulose-based materials aimed at biomedical applications such as wound dressings.

Synergistic Effect of Chemical Penetration Enhancers on Lidocaine Permeability Revealed by Coarse-Grained Molecular Dynamics Simulations

The search for new formulations for transdermal drug delivery (TDD) is an important field in medicine and cosmetology. Molecules with specific physicochemical properties which can increase the permeability of active ingredients across the stratum corneum (SC) are called chemical penetration enhancers (CPEs), and it was shown that some CPEs can act synergistically. In this study, we performed coarse-grained (CG) molecular dynamics (MD) simulations of the lidocaine delivery facilitated by two CPEs-linoleic acid (LA) and ethanol-through the SC model membrane containing cholesterol, N-Stearoylsphingosine (DCPE), and behenic acid. In our simulations, we probed the effects of individual CPEs as well as their combination on various properties of the SC membrane and the lidocaine penetration across it. We demonstrated that the addition of both CPEs decreases the membrane thickness and the order parameters of the DPCE hydrocarbon chains. Moreover, LA also enhances diffusion of the SC membrane components, especially cholesterol. The estimated potential of mean force (PMF) profiles for the lidocaine translocation across SC in the presence/absence of two individual CPEs and their combination demonstrated that while ethanol lowers the free energy barrier for lidocaine to enter SC, LA decreases the depth of the free energy minima for lidocaine inside SC. These two effects supposedly result in synergistic penetration enhancement of drugs. Altogether, the present simulations provide a detailed molecular picture of CPEs' action and their synergistic effect on the penetration of small molecular weight therapeutics that can be beneficial for the design of novel drug and cosmetics formulations.

Acidic pH Is Required for the Multilamellar Assembly of Skin Barrier Lipids In Vitro

Lipid membrane remodeling belongs to the most fundamental processes in the body. The skin barrier lipids, which are ceramide dominant and highly rigid, must attain an unusual multilamellar nanostructure with long periodicity to restrict water loss and prevent the entry of potentially harmful environmental factors. Our data suggest that the skin acid mantle, apart from regulating enzyme activities and keeping away pathogens, may also be a prerequisite for the multilamellar assembly of the skin barrier lipids. Atomic force microscopy on monolayers composed of synthetic or human stratum corneum lipids showed multilayer formation (approximately 10-nm step height) in an acidic but not in a neutral environment. X-ray diffraction, Fourier transform infrared spectroscopy, and permeability studies showed markedly altered lipid nanostructure and increased water loss at neutral pH compared with that at acidic pH. These findings are consistent with the data on the altered organization of skin lipids and increased transepidermal water loss under conditions such as inadequate skin acidification, for example, in neonates, the elderly, and patients with atopic dermatitis.

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.

Very Long Chain Lipids Favor the Formation of a Homogeneous Phase in Stratum Corneum Model Membranes

The stratum corneum (SC), the outermost layer of mammal epidermis, acts as a barrier dictating the rate of absorption of exogenous molecules through the skin, as well as to prevent excessive water loss from the body. The SC consists of protein-rich corneocytes embedded into a complex lipid mixture. The lipid fraction is mainly constituted of an equimolar mixture of ceramides (Cer), free fatty acids (FFA), and cholesterol (Chol), forming a solid phase in the intracellular space; this lipid phase is supposed to play a fundamental role in the SC barrier function. An unusual characteristic of this biological membrane is that its lipids generally bear very long acyl chains, with the 24-carbon long ones being the most abundant. In this work, we used Raman microspectroscopy and infrared spectroscopy to study the influence of the acyl chain length on the lipid mixing properties in SC model membranes. Our results revealed that the combination of ceramides and FFA bearing a very long chain is required for the formation of homogeneous lipid mixtures, while lipids with shorter chains (16-carbon and 20-carbon atom long) lead to domains with micrometer dimensions. It is proposed that the biological machinery necessary for acyl chain elongation occurring at the mammalian skin level is required to inhibit lipid phase separation, a critical feature in the proper barrier functioning.

The Sphingosine and Acyl Chains of Ceramide [NS] Show Very Different Structure and Dynamics That Challenge Our Understanding of the Skin Barrier

The lipid phase of the uppermost human skin layer is thought to comprise highly rigid lipids in an orthorhombic phase state to protect the body against the environment. By synthesizing sphingosine-d28 deuterated N-lignoceroyl-d-erythro-sphingosine (ceramide [NS]), we compare the structure and dynamics of both chains of that lipid in biologically relevant mixtures using X-ray diffraction, 2 H NMR analysis, and infrared spectroscopy. Our results reveal a substantial fraction of sphingosine chains in a fluid and dynamic phase state at physiological temperature. These findings prompt revision of our current understanding of the skin lipid barrier, where an extended ceramide [NS] conformation is preferred and a possible domain structure is proposed. Mobile lipid chains may be crucial for skin elasticity and the translocation of physiologically important molecules.

Raman and AFM-IR chemical imaging of stratum corneum model membranes

Stratum corneum (SC), the outermost layer of the epidermis, is the primary barrier to percutaneous absorption. The diffusion of substances through the skin occurs through the SC lipid fraction, which is essentially constituted of an equimolar mixture of ceramides, free fatty acids, and cholesterol. The lipid constituents of SC are mainly forming continuous multilamellar membranes in the solid/crystalline state. However, recent findings suggest the presence of a highly disordered (liquid) phase formed by the unsaturated C18 chain of ceramide EOS, surrounded by a highly ordered lipid environment. The aim of the present work was to study the lipid spatial distribution of model SC membranes composed of ceramide EOS, ceramide NS, a mixture of free fatty acids, and cholesterol, using Raman microspectroscopy and AFM-IR spectroscopy techniques. The enhanced spatial resolution at the tens of nanometers scale of the AFM-IR technique revealed that the lipid matrix is overall homogeneous, with the presence of small, slightly enriched, and depleted regions in a lipid component. No liquid domains of ceramide EOS were observed at this scale, a result that is consistent with the model proposing that the oleate nanodrops are concentrated in the central layer of the three-layer organization of the SC membranes forming the long periodicity phase. In addition, both Raman microspectroscopy and AFM-IR techniques confirmed the fluid nature of the unsaturated chain of ceramide EOS while the rest of the lipid matrix was found highly ordered.

Arrangement of Ceramides in the Skin: Sphingosine Chains Localize at a Single Position in Stratum Corneum Lipid Matrix Models

Understanding the structure of the stratum corneum (SC) is essential to understand the skin barrier process. The long periodicity phase (LPP) is a unique trilayer lamellar structure located in the SC. Adjustments in the composition of the lipid matrix, as in many skin abnormalities, can have severe effects on the lipid organization and barrier function. Although the location of individual lipid subclasses has been identified, the lipid conformation at these locations remains uncertain. Contrast variation experiments via small-angle neutron diffraction were used to investigate the conformation of ceramide (CER) N-(tetracosanoyl)-sphingosine (NS) within both simplistic and porcine mimicking LPP models. To identify the lipid conformation of the twin chain CER NS, the chains were individually deuterated, and their scattering length profiles were calculated to identify their locations in the LPP unit cell. In the repeating trilayer unit of the LPP, the acyl chain of CER NS was located in the central and outer layers, while the sphingosine chain was located exclusively in the middle of the outer layers. Thus, for the CER NS with the acyl chain in the central layer, this demonstrates an extended conformation. Electron density distribution profiles identified that the lipid structure remains consistent regardless of the lipid's lateral packing phase, this may be partially due to the anchoring of the extended CER NS. The presented results provide a more detailed insight on the internal arrangement of the LPP lipids and how they are expected to be arranged in healthy skin.

Coexistence of Lipid Phases Stabilizes Interstitial Water in the Outer Layer of Mammalian Skin

The lipid matrix in the outer layer of mammalian skin, the stratum corneum, has been previously investigated by multiple biophysical techniques aimed at identifying hydrophilic and lipophilic pathways of permeation. Although consensus is developing over the microscopic structure of the lipid matrix, no molecular-resolution model describes the permeability of all chemical species simultaneously. Using molecular dynamics simulations of a model mixture of skin lipids, the self-assembly of the lipid matrix lamellae has been studied. At higher humidity, the resulting lamellar phase is maintained by partitioning excess water into isolated droplets of controlled size and spatial distribution. The droplets may fuse together to form intralamellar water channels, thereby providing a pathway for the permeation of hydrophilic species. These results reconcile competing data on the outer skin's structure and broaden the scope of molecular-based methods to improve the safety of topical products and to advance transdermal drug delivery.