Molecular Dynamics Simulation Study of Skin Lipids: Effects of the Molar Ratio of Individual Components over a Wide Temperature Range

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.

Mechanisms of the Drug Penetration Enhancer Propylene Glycol Interacting with Skin Lipid Membranes

Very few drugs have the necessary physicochemical properties to cross the skin's main permeability barrier, the stratum corneum (SC), in sufficient amounts. Propylene glycol (PG) is a chemical penetration enhancer that could be included in topical formulations in order to overcome the barrier properties of the skin and facilitate the transport of drugs across it. Experiments have demonstrated that PG increases the mobility and disorder of SC lipids and may extract cholesterol from the SC, but little is known about the molecular mechanisms of drug permeation enhancement by PG. In this work, we have performed molecular dynamics (MD) simulations to investigate the molecular-level effects of PG on the structure and properties of model SC lipid bilayers. The model bilayers were simulated in the presence of PG concentrations over the range of 0-100% w/w PG, using both an all-atom and a united atom force field. PG was found to localize in the hydrophilic headgroup regions at the bilayer interface, to occupy the lipid-water hydrogen-bonding sites, and to slightly increase lipid tail disorder in a concentration-dependent manner. We showed with MD simulation that PG enhances the permeation of small molecules such as water by interacting with the bilayer interface; the results of our study may be used to guide the design of formulations for transdermal drug delivery with enhanced skin permeation, as well as topical formulations and cosmetic products.

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.

Formulation and Skin Permeation of Active-Loaded Lipid Nanoparticles: Evaluation and Screening by Synergizing Molecular Dynamics Simulations and Experiments

Encapsulation into nanoparticles (NPs) is a potential method to deliver pharmaceutical/cosmetic actives deep into the skin. However, understanding the NP formulations and underlying mechanism of active delivery to skin has scarcely been studied. We report a simulation platform that screens, evaluates, formulates, and provides atomic-resolution interpretation of NP-based formulations, and reveals the active permeation mechanism from NPs to skin. First, three actives, namely, ferulic acid (FA), clotrimazole (CZE), and tretinoin (TTN), and five lipid excipients' (Compritol, Precirol, Geleol, Gelot, Gelucire) combinations were screened by MD simulations for the best pairs. For each suggested pair, the actual active and lipid compositions for the synthesis of stable NP formulations were then obtained by experiments. MD simulations demonstrate that in NP formulations, FA and CZE actives are present at the surface of the NPs, whereas TTN actives are present at both the surface and interior of the NP core. The NP shapes obtained by simulation perfectly match with experiments. For each NP, separate MD simulations illustrate that active-loaded NPs approach the skin surface quickly, and then actives translocate from NP surface to skin surface followed by penetration of NPs through skin. The driving force for the translocation which initiates during the penetration process, is the stronger active-skin interaction compared to active-NP interaction. Permeation free energy indicates spontaneous transfer of actives from solution phase to the surface of the skin bilayer. The free energy barriers are increased in the order of FA < TTN < CZE. Significantly lower diffusions of actives are obtained in the main barrier region compared to bulk, and the average diffusion coefficients of actives are in the same order of magnitude (∼10^-6 cm²/s). The estimated permeability coefficients (log P) of actives are mainly governed by free energy barriers. The study would facilitate the development of novel lipid-based NP formulations for personal-care/pharmaceutical applications.

Lycopene, but not zeaxanthin, serves as a skeleton for the formation of an orthorhombic organization of intercellular lipids within the lamellae in the stratum corneum: Molecular dynamics simulations of the hydrated ceramide NS bilayer model

Carotenoids play an important role in the protection of biomembranes against oxidative damage. Their function depends on the surroundings and the organization of the lipid membrane they are embedded in. Carotenoids are located parallel or perpendicular to the surface of the lipid bilayer. The influence of carotenoids on the organization of the lipid bilayer in the stratum corneum has not been thoroughly considered. Here, the orientation of the exemplary cutaneous carotenoids lycopene and zeaxanthin in a hydrated ceramide NS24 bilayer model and the influence of carotenoids on the lateral organization of the lipid bilayer model were studied by means of molecular dynamics simulations for 32 °C and 37 °C. The results confirm that lycopene is located parallel and zeaxanthin perpendicular to the surface of the lipid bilayer. The lycopene-loaded lipid bilayer appeared to have a strong orthorhombic organization, while zeaxanthin-loaded and pure lipid bilayers were organized in a disordered hexagonal-like and liquid-like state, respectively. The effect is stronger at 32 °C compared to 37 °C based on p-values. Therefore, it was assumed that carotenoids without hydroxyl polar groups in their structure facilitate the formation of the orthorhombic organization of lipids, which provides the skin barrier function. It was shown that the distance between carotenoid atoms matched the distance between atoms in the lipids, indicating that parallel located carotenoids without hydroxyl groups serve as a skeleton for lipid membranes inside the lamellae. The obtained results provide reasonable prediction of the overall qualitative properties of lipid model systems and show the importance of parallel-oriented carotenoids in the development and maintenance of the skin barrier function.

Molecular dynamics simulations of the human ocular lens with age and cataract

The human ocular lens consists primarily of elongated, static fibers characterized by high stability and low turnover, which differ dramatically in their composition and properties from other biological membranes. Cholesterol (Chol) and sphingolipids (SL) are present at high concentrations, including saturated SLs, such as dihyrosphingomyelin (DHSM). Past molecular dynamics simulations demonstrated that the presence of DHSM and high Chol concentration contributes to higher order in lipid membranes. This current study simulated more complex models of human lens membranes. Models were developed representing physiological compositions in cataractous lenses aged 74 ± 6 years and in healthy lenses aged 22 ± 4, 41 ± 6, and 69 ± 3 years. With older age, Chol and ceramide concentrations increase and glycerophospholipid concentration decreases. With cataract, ceramide concentration increases and Chol and glycerophospholipid concentrations decrease. Surface area per lipid, deuterium order parameters (S_CD), sterol tilt angle, electron density profiles, bilayer thickness, chain interdigitation, two-dimensional radial distribution functions (2D-RDF), lipid clustering, and hydrogen bonding were calculated for all simulations. All systems exhibited low surface area per lipid and high bilayer thickness, indicative of strong vertical packing. S_CD parameters suggest similarly, with saturated tails in the hydrophobic core of the membrane having elevated order. Vertical packing and acyl tail order increased with both age and cataract condition. Lateral diffusion decreased with age and cataracts, with the older and cataractous models demonstrating increased long-range structure by the 2D-RDF analysis. In future work examining the membrane proteins of the lens, these models can serve as a physiologically accurate representation of the lens lipidome.

Investigating the nanostructure of a CER[NP]/CER[AP]-based stratum corneum lipid matrix model: A combined neutron diffraction & molecular dynamics simulations approach

The human skin provides a physiochemical and biological protective barrier due to the unique structure of its outermost layer known as the Stratum corneum. This layer consists of corneocytes and a multi-lamellar lipid matrix forming a composite, which is a major determining factor for the barrier function of the Stratum corneum. A substantiated understanding of this barrier is necessary, as controlled breaching or modulation of the same is also essential for various health and personal care applications such as topical drug delivery and cosmetics to a name few. In this study, we discuss the state-of-the-art of neutron diffraction techniques, using specifically deuterated lipids, combined with the information obtained from molecular models using molecular dynamics simulations, to understand the structure and barrier function of the Stratum corneum lipid matrix. As an example, the effect of ceramide concentration on a lipid lamella system consisting of CER[NP]/CER[AP]/Cholesterol/free fatty acid (deprotonated) is studied. This study demonstrates the usefulness of the combined approach of neutron diffraction and molecular dynamics simulations for effective analysis of the model systems created for the Stratum corneum lipid matrix. The optimization of force fields by comparison with experimental data is furthermore an important step in the direction of providing a predictive quality.

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.

Multiscale Simulation of Ternary Stratum Corneum Lipid Mixtures: Effects of Cholesterol Composition

Molecular dynamics simulations of mixtures of the ceramide nonhydroxy-sphingosine (NS), cholesterol, and a free fatty acid are performed to gain molecular-level understanding of the structure of the lipids found in the stratum corneum layer of skin. A new coarse-grained force field for cholesterol was developed using the multistate iterative Boltzmann inversion (MS-IBI) method. The coarse-grained cholesterol force field is compatible with previously developed coarse-grained force fields for ceramide NS, free fatty acids, and water and validated against atomistic simulations of these lipids using the CHARMM force field. Self-assembly simulations of multilayer structures using these coarse-grained force fields are performed, revealing that a large fraction of the ceramides adopt extended conformations, which cannot occur in the single bilayer in water structures typically studied using molecular simulation. Cholesterol fluidizes the membrane by promoting packing defects, and an increase in cholesterol content is found to reduce the bilayer thickness due to an increase in interdigitation of the C₂₄ lipid tails, consistent with experimental observations. Using a reverse-mapping procedure, a self-assembled coarse-grained multilayer system is used to construct an equivalent structure with atomistic resolution. Simulations of this atomistic structure are found to closely agree with experimentally derived neutron scattering length density profiles. Significant interlayer hydrogen bonding is observed in the inner layers of the atomistic multilayer structure that are not found in the outer layers in contact with water or in equivalent bilayer structures. This work highlights the importance of simulating multilayer structures, as compared to the more commonly studied bilayer systems, to enable more appropriate comparisons with multilayer experimental membranes. These results also provide validation of the efficacy of the MS-IBI derived coarse-grained force fields and the framework for multiscale simulation.

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.

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.

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.

How Do Urea and Trimethylamine N-Oxide Influence the Dehydration-Induced Phase Transition of a Lipid Membrane?

Living organisms are often exposed to extreme dehydration, which is detrimental to the structure and function of the cell membrane. The lipid membrane undergoes fluid-to-gel phase transition due to dehydration and thus loses fluidity and functionality. To protect the fluid phase of the bilayer these organisms adopt several strategies. Enhanced production of small polar organic solutes (also called osmolytes) is one such strategy. Urea and trimethylamine N-oxide (TMAO) are two osmolytes found in different organisms combating osmotic stress. Previous experiments have found that both these osmolytes have strong effects on lipid membrane under different hydration conditions. Urea prevents the dehydration-induced phase transition of the lipid membrane by directly interacting with the lipids, while TMAO does not inhibit the phase transition. To provide atomistic insights, we have carried out all-atom molecular dynamics (MD) simulation of a lipid membrane under varying hydration levels and studied the effect of these osmolytes on different structural and dynamic properties of the membrane. This study suggests that urea significantly inhibits the dehydration-induced fluid-to-gel phase transition by strongly interacting with the lipid membrane via hydrogen bonds, which balances the reduced lipid hydration due to the decreasing water content. In contrast, TMAO is excluded from the membrane surface due to unfavorable interaction with the lipids. This induces further dehydration of the lipids which reinforces the fluid-to-gel phase transition. We have also studied the counteractive role of TMAO on the effect of urea on lipid membrane when both the osmolytes are present. TMAO draws some urea molecules out of the membrane and thereby reduces the effect of urea on the lipid membrane at lower hydration levels. This is similar to the counteraction of urea's deleterious effects on protein by TMAO. All these observations are consistent with the experimental results and thus provide deep molecular insights into the role of these osmolytes in protecting the fluid phase of the membrane, the key survival strategy against osmotic-stress-induced dehydration.

Two- and Three-Dimensional Physical-Chemical Characterization of CER[AP]: A Study of Stereochemistry and Chain Symmetry

The stratum corneum represents the first skin barrier against chemical and physical damage. These unique properties are based on its peculiar lipid composition with ceramides (CERs) as the main protagonists. In this study, the structural and chemical properties of the α-OH phytosphingosine [AP] CER class have been investigated. α-OH CERs are present in the stratum corneum in their d-forms; however, in most model systems the diastereomer mixture with the synthetically produced l-form is used. The d-form is well-known to form a hydrogen bonding network that helps to reduce the permeability of the lipid matrix, while the l-form does not show any hydrogen bonding network formation. In this paper, 2D (monolayers) and 3D (aqueous dispersions) models have been used to thoroughly study the physical-chemical behaviors of CER[AP] diastereomers taking into account how the symmetry of the chain pattern influences the behavior of the molecules. The chains of both diastereomers arrange in an oblique unit cell, but only the d-CER[AP] forms a supramolecular lattice (subgel phase) in both model systems. Interestingly, the chain pattern does not play any role in structure formation since the hydrogen bonding network dictates the packing properties. The 1:1 mixture of the diastereomers phase separates into two domains: one is composed of practically pure d-form and the other one is composed of a mixture of the l-form with a certain amount of d-form molecules.

Molecular dynamics simulations to elucidate translocation and permeation of active from lipid nanoparticle to skin: complemented by experiments

One of the most realistic approaches for delivering actives (pharmaceuticals/cosmetics) deep into skin layers is encapsulation into nanoparticles (NPs). Nonetheless, molecular-level mechanisms related to active delivery from NPs to the skin have scarcely been studied despite the large number of synthesis and characterization studies. We herein report the underlying mechanism of active translocation and permeation through the outermost layer of skin, the stratum corneum (SC), via molecular dynamics (MD) simulations complemented by experimental studies. A SC molecular model is constructed using current state-of-the-art methodology via incorporating the three most abundant skin lipids: ceramides, free fatty acids, and cholesterol. As a potent antioxidant, ferulic acid (FA) is used as the model active, and it is loaded into Gelucire 50/13 NP. MD simulations elucidate that, first, FA-loaded NP approaches the skin surface quickly, followed by slight penetration and adsorption onto the upper skin surface; FA then translocates from the NP surface to the skin surface due to stronger NP-skin interactions compared to the FA-NP interactions; then, once FA is released onto the skin surface, it slowly permeates deep into the skin bilayer. Both the free energy and resistance to permeation not only indicate the spontaneous transfer of FA from the bulk to the skin surface, but they also reveal that the main barrier against permeation exists in the middle of the lipid hydrophobic tails. Significantly lower diffusion of FA is obtained in the main barrier region compared to the bulk. The estimated permeability coefficient (log P) values are found to be higher than the experimental values. Importantly, the permeation process evaluated via MD simulations perfectly matches with experiments. The study suggests a molecular simulation platform that provides various crucial insights relating to active delivery from loaded NP to skin, and it could facilitate the design and development of novel NP-based formulations for transdermal delivery and the topical application of drugs/cosmetics.

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.

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.

A Compliant Ionic Adhesive Electrode with Ultralow Bioelectronic Impedance

Simultaneous implementation of high signal-to-noise ratio (SNR) but low crosstalk is of great importance for weak surface electromyography (sEMG) signals when precisely driving a prosthesis to perform sophisticated activities. However, due to gaps with the curved skin during muscle contraction, many electrodes have poor compliance with skin and suffer from high bioelectrical impedance. This causes serious noise and error in the signals, especially the signals from low-level muscle contractions. Here, the design of a compliant electrode based on an adhesive hydrogel, alginate-polyacrylamide (Alg-PAAm) is reported, which eliminates those large gaps through the strong electrostatic interaction and abundant hydrogen bond with the skin. The obtained compliant electrode, having an ultralow bioelectrical impedance of ≈20 kΩ, can monitor even 2.1% maximal voluntary contraction (MVC) of muscle. Furthermore, benefiting from the high SNR of >5:1 at low-level MVC, the crosstalk from irrelevant muscle is minimized through reducing the electrode size. Finally, a prosthesis is successfully demonstrated to precisely grasp a needle based on a 9 mm² Alg-PAAm compliant electrode. The strategy to design such compliant electrodes provides the potential for improving the quality of dynamically weak sEMG signals to precisely control prosthesis in performing purposefully dexterous activity.

State of the Art in Stratum Corneum Research. Part II: Hypothetical Stratum Corneum Lipid Matrix Models

This review is the second part of a series which presents the state of the art in stratum corneum (SC) lipid matrix (LM) research in depth. In this part, the various hypothetical models which were developed to describe the structure and function of the SC LM as the skin's barrier will be discussed. New as well as a cumulative assortment of older results which change the view on the different models are considered to conclude how well the different models are holding up today. As a final conclusion, a model, factoring in as much of the known data as possible, is concluded, unifying the varying different models into one which can be developed further, as new results are found in the future. So far, the model is described with a single crystalline or gel-like phase with a certain amount of nanocrystallites of concentrated ceramides (CERs) and free fatty acids and more fluid nanodomains caused by a fluidizing effect of the cholesterol. These domains are dynamically resolved and reformed and do not impair the barrier function. The chain conformation is not completely clear yet; however, an equilibrium of fully extended and hairpin-folded CERs with ratios depending on the properties of each individual CER species is proposed as most likely. An overlapping middle layer as described for the tri-layer model in part I of this series would be present for both conformations. The macroscopic broad-narrow-broad layering, observed in electron micrographs, is explained by an external templating by the lipid envelope, and an internal templating by short and long lipid chains each preferentially show a homophilic association, forming thicker and thinner bilayers, respectively. The degree of influence of the very long ω-hydroxy-CERs is discussed as well.

Molecular Dynamics Simulations of Micelle Properties and Behaviors of Sodium Lauryl Ether Sulfate Penetrating Ceramide and Phospholipid Bilayers

Molecular dynamics (MD) simulations with the umbrella sampling (US) method were used to investigate the properties of micelles formed by sodium lauryl ether sulfate with two ether groups (SLE2S) and behaviors of corresponding surfactants transferring from micelles to ceramide and DMPC bilayer surfaces. Average micelle radii based on the Einstein-Smoluchowski and Stokes-Einstein relations showed excellent agreement with those measured by dynamic light scattering, while those obtained by evaluating the gyration radius or calculating the distance between the micelle sulfur atoms and center of mass overestimate the radii. As an SLE2S micelle was pulled down to the ceramide bilayer surface in a 400 ns constant-force steered MD (cf-SMD) simulation, the micelle was partially deformed on the bilayer surface, and several SLE2S surfactants easily were partitioned from the micelle into the ceramide bilayer. In contrast, a micelle was not deformed on the DMPC bilayer surface, and SLE2S surfactants were not transferred from the micelle to the DMPC bilayer. Potential of mean force (PMF) calculations revealed that the Gibbs free energy required for an SLE2S surfactant monomer to transfer from a micelle to bulk water can be compensated by decreased Gibbs free energy when an SLE2S monomer transfers into the ceramide bilayer from bulk water. In addition, micelle deformation on the ceramide bilayer surface can reduce the Gibbs free energy barrier required for a surfactant to escape the micelle and help the surfactant partition from the micelle into the ceramide bilayer. An SLE2S surfactant partitioning into the ceramide bilayer is attributed to hydrogen bonding and favorable interactions between the hydrophilic surfactant head and ceramide molecules, which are more dominant than the dehydration penalty during bilayer insertion. Such interactions between surfactant and lipid molecule heads are considerably reduced in DMPC bilayers owing to dielectric screening by water molecules deep inside the head/tail boundary between the DMPC bilayer. This computational work demonstrates the distinct behavior of SLE2S surfactant micelles on ceramide and DMPC bilayer surfaces in terms of variation in Gibbs free energy, which offers insight into designing surfactants used in transdermal drug delivery systems and cosmetics.

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 Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Size-dependent interaction of hydrophilic/hydrophobic ligand functionalized cationic and anionic nanoparticles with lipid bilayers

We study the nature of nanoparticle (NPs)-membrane interaction as a function of nanoparticle size for different functionalization using molecular dynamics simulation. Zinc sulphide quantum dots of size, 2 nm and 4 nm are used as model NPs, and DLPC and DPPC lipid bilayers are used as model membranes. We use coarse-grained polarizable MARTINI model (MPW) to simulate the NPs and lipid bilayers. Our simulation results show that uncharged bare NPs penetrate the lipid bilayers and embed themselves within the hydrophobic core of the bilayer both in the gel and fluid phases. NPs of size 4 nm are shown to disrupt the bilayer. The bilayer recovers from the damages caused by smaller NPs of size 2 nm. In case of either purely hydrophilic or hybrid (with hydrophilic/hydrophobic ratio of 2:1) ligand-functionalized NPs of smaller size (shell size 2 nm), only cationic NPs bind to the bilayer. However, for larger NPs with a shell size of 4 nm, both anionic and cationic hybrid functionalized NPs bind to the bilayer. The performance of standard Martini (SM) force field for the charged NP/bilayer systems has also been tested and compared with the results obtained using MPW model. Although the overall trend that the cationic NPs interact strongly with the bilayers than their anionic counterparts has been captured correctly using SM, the adsorption behaviour of the functionalized NPs differ significantly in the SM force field. The interaction of anionic NPs with both fluid and gel bilayers has been observed to be least accurately represented in the SM force field.

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.

Effect of Ceramide Tail Length on the Structure of Model Stratum Corneum Lipid Bilayers

Lipid bilayers composed of non-hydroxy sphingosine ceramide (CER NS), cholesterol (CHOL), and free fatty acids (FFAs), which are components of the human skin barrier, are studied via molecular dynamics simulations. Since mixtures of these lipids exist in dense gel phases with little molecular mobility at physiological conditions, care must be taken to ensure that the simulations become decorrelated from the initial conditions. Thus, we propose and validate an equilibration protocol based on simulated tempering, in which the simulation takes a random walk through temperature space, allowing the system to break out of metastable configurations and hence become decorrelated from its initial configuration. After validating the equilibration protocol, which we refer to as random-walk molecular dynamics, the effects of the lipid composition and ceramide tail length on bilayer properties are studied. Systems containing pure CER NS, CER NS + CHOL, and CER NS + CHOL + FFA, with the CER NS fatty acid tail length varied within each CER NS-CHOL-FFA composition, are simulated. The bilayer thickness is found to depend on the structure of the center of the bilayer, which arises as a result of the tail-length asymmetry between the lipids studied. The hydrogen bonding between the lipid headgroups and with water is found to change with the overall lipid composition, but is mostly independent of the CER fatty acid tail length. Subtle differences in the lateral packing of the lipid tails are also found as a function of CER tail length. Overall, these results provide insight into the experimentally observed trend of altered barrier properties in skin systems where there are more CERs with shorter tails present.

Structural and barrier properties of the skin ceramide lipid bilayer: a molecular dynamics simulation study

Skin provides excellent protection against the harsh external environment and foreign substances. The lipid matrix of the stratum corneum, which contains various kinds of ceramides, plays a major role in the barrier function of the skin. Here we report a study of the effects of ceramide type on the structural and transport properties of ceramide bilayers using molecular dynamics (MD) simulations. Specifically, the effects of headgroup chemistry (number and positions of hydroxyl groups) and tail structure (unsaturation of the sphingoid moiety) on the structural and transport properties of various ceramide bilayers at 310 K were analyzed. Theoretical results for structural properties such as area per lipid, bilayer thickness, lateral arrangement, order parameter, and hydrogen bonding are reported here and compared with corresponding experimental data. Our study revealed that the presence of a double bond disrupts the bilayer packing, which leads to a low area compressibility modulus, a large area per lipid, and low bilayer thickness. Furthermore, the effect of structural changes on water permeation was studied using steered MD simulations. Water permeation was found to be influenced by headgroup polarity, chain packing, and the ability of the water to hydrogen bond with the ceramides. The molecular-level information obtained from the current study should aid the design of mixed bilayer systems with desired properties and provide the basis for the development of higher order coarse-grained models.

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.

Mechanisms of lipid extraction from skin lipid bilayers by sebum triglycerides

Microsecond computations identify the pathways leading to the extraction of skin lipids by sebum triglycerides and the associated energetic costs.

State of the art in Stratum Corneum research: The biophysical properties of ceramides

This review is summarizing an important part of the state of the art in stratum corneum research. A complete overview on discoveries about the general biophysical and physicochemical properties of the known ceramide species' is provided. The ceramides are one of the three major components of the lipid matrix and mainly govern its properties and structure. They are shown to exhibit very little redundancy, despite the minor differences in their chemical structure. The results are discussed, compared to each other as well as the current base of knowledge. New interesting aspects and concepts are concluded or suggested. A novel interpretation of the 3-dimensional structure of the lipid matrix and its influence on the barrier function will be discussed. The most important conclusion is the presentation of a new and up to date theoretical model of the nanostructure of the short periodicity phase. The model suggests three perpendicular layers: The rigid head group region, the rigid chain region and, a liquid-like overlapping middle layer. The general principle of the skin barrier function is highlighted in regard to this structure and the ceramides biophysical and physicochemical properties. As a result of these considerations, the entropy vs. enthalpy principle is introduced, shedding light on the function as well as the effectiveness of the skin barrier. Additionally, general ideas to effectively overcome this barrier principle for dermal and transdermal delivery of actives or how to use it for specific targeting of the stratum corneum are proposed.

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.

The long periodicity phase (LPP) controversy part I: The influence of a natural-like ratio of the CER[EOS] analogue [EOS]-br in a CER[NP]/[AP] based stratum corneum modelling system: A neutron diffraction study

This study used neutron diffraction to investigate a ceramide-[NP] C24/[AP] C24 /[EOS]-br C30/cholesterol/lignoceric acid (0.6: 0.3: 0.1: 0.7: 1) based stratum corneum modelling system. By adding specifically deuterated ceramides-[NP]-D₃, [AP]-D₃, and [EOS]-br-D₃, detailed information on the lamellar and the nanostructure of the system was obtained. For the short periodicity phase a natural-like lamellar repeat distance of 5.47 ± 0.02 nm was observed, similar to the [NP]/[AP] base system without the [EOS]-br. Unlike in this system the ceramides here were slightly tilted, hinting towards a slightly less natural arrangement. Due to the deuteration it was possible to observe that the long ceramide chains were overlapping in the lamellar mid-plane. This is considered to be an important feature for the natural stratum corneum. Despite the presence of a ceramide [EOS] analogue - able to form a long phase arrangement - no distinct long periodicity phase was formed, despite a slightly higher than natural ω-acyl ceramide ratio of 10 mol%. The deuterated variant of this ceramide determined that the very long ceramide was integrated into the short periodicity phase, spanning multiple layers instead. The - compared to the base system - unchanged repeat distance highlights the stability of this structure. Furthermore, the localisation of the very long ceramide in the short periodicity phase indicates the possibility of a crosslinking effect and thus a multilayer stabilizing role for the ceramide [EOS]. It can be concluded, that additionally to the mere presence of ceramide-[EOS] more complex conditions have to be met in order to form this long phase. This has to be further investigated in the future.

Impact of the ceramide subspecies on the nanostructure of stratum corneum lipids using neutron scattering and molecular dynamics simulations. Part I: impact of CER[NS]

For this study mixtures based on the ceramides [NS] (NS = non-hydroxy-sphingosine) and [AP] (AP = α-hydroxy-phytosphingosine) in a 2:1 and 1:2 ratio, together with cholesterol and lignoceric acid, were investigated. These mixtures are modelling the uppermost skin layer, the stratum corneum. Neutron diffraction, utilizing specifically deuterated ceramide molecules, was used to obtain a maximum amount of experimental detail. Highly detailed molecular dynamics simulations were used to generate even more information from the experimental data. It was possible to observe a single lamellar phase for both systems. They had a lamellar repeat distance of 5.43 ± 0.05 nm for the [NS]/[AP] 2:1 and a slightly shorter one of 5.34 ± 0.05 nm for the 1:2 system. The structure and water content was uninfluenced by excess humidity. Both the experimental and simulation data indicated slightly tilted ceramides, with their C24 chains overlapping in the lamellar mid-plane. This arrangement is well comparable to systems investigated before. The structure of both systems, except for the differing repeat distance, looks similar at first. However, on a smaller scale there were various distinct differences, demonstrating only low redundancy between the different ceramide species, despite only minor chemical differences. The mainly ceramide [AP] determined 1:2 system has a slightly smaller repeat distance. This is a result of a tighter arrangement of the lipids chain along the bilayer normal and increased overlapping of the long chains in the lamellar middle. For the CER[NS] some novel features could be shown, despite it being the overall most investigated ceramide. These include the low adaptability to changed lateral interactions, leading to an increased chain opening. This effect could explain its low miscibility with other lipids. The investigated model systems allows it to directly compare results from the literature which have used ceramide [NS] to the most recent studies using the phytosphingosine ceramides such as ceramide [AP].

Electroporation of Skin Stratum Corneum Lipid Bilayer and Molecular Mechanism of Drug Transport: A Molecular Dynamics Study

The electroporation technique has been used significantly to increase drug permeation through the skin. This technique relies on the application of short-timed (microseconds to millisecond) electric fields (generally, order of 50--300 V) on the skin to create microscopic pores. However, the molecular mechanism of pore formation, resulting in an enhanced flux of active molecules through the skin, remains poorly understood. In this study, extensive atomistic molecular dynamics simulation of skin lipids [made up of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA)] has been performed at various external electric fields. We show for the first time the pore formation in the skin lipid bilayer during electroporation. We show the effect of the applied external electrical field (0.6-1.0 V/nm) on the pore formation dynamics in the lipid bilayer of different sizes (154, 616, and 2464 lipids) and compositions (CER/CHOL/FFA, 1:0:0, 1:0:1, 1:1:0, 1:1:1). The pore formation and resealing kinetics were different and were found to be highly dependent on the composition of the skin lipid bilayer. The pore formation time decreased with increase in the bilayer size. The pore sustaining electric field was found to be in the range of 0.20-0.25 V/nm for equimolar CER, CHOL, and FFA lipid bilayers. The skin lipid bilayer (1:1:1) sealed itself within 20 ns after the removal of the external electric field. We also present the molecular mechanism of enhancement of drug permeation in the presence of external field as compared to the passive diffusion. The molecular-level understanding obtained here could help in optimizing/designing the electroporation experiments for effective drug delivery. For a given skin composition and size of the drug molecule, the combination of pore formation time and pore growth model can be used to know a priori the desired electric field and time for the application of the electric field.

Molecular dynamics simulations of stratum corneum lipid mixtures: A multiscale perspective

The lipid matrix of the stratum corneum (SC) layer of skin is essential for human survival; it acts as a barrier to prevent rapid dehydration while keeping potentially hazardous material outside the body. While the composition of the SC lipid matrix is known, the molecular-level details of its organization are difficult to infer experimentally, hindering the discovery of structure-property relationships. To this end, molecular dynamics simulations, which give molecular-level resolution, have begun to play an increasingly important role in understanding these relationships. However, most simulation studies of SC lipids have focused on preassembled bilayer configurations, which, owing to the slow dynamics of the lipids, may influence the final structure and hence the calculated properties. Self-assembled structures would avoid this dependence on the initial configuration, however, the size and length scales involved make self-assembly impractical to study with atomistic models. Here, we report on the development of coarse-grained models of SC lipids designed to study self-assembly. Building on previous work, we present the interactions between the headgroups of ceramide and free fatty acid developed using the multistate iterative Boltzmann inversion method. Validation of the new interactions is performed with simulations of preassembled bilayers and good agreement between the atomistic and coarse-grained models is found for structural properties. The self-assembly of mixtures of ceramide and free fatty acid is investigated and both bilayer and multilayer structures are found to form. This work therefore represents a necessary step in studying SC lipid systems on multiple time and length scales.

In-silico design of nanoparticles for transdermal drug delivery application

Nanoparticles are used in the medical field for various applications like cell imaging, drug delivery, gene and si-RNA delivery, to name a few. Designing nanoparticles for a given application, purely based on the trial and error experimentation, requires a lot of time and effort. In this study we show that computer simulations could help in designing nanoparticles for drug delivery thus reducing the time and cost associated with their design, development and deployment. The permeation of nanoparticles, having various surface chemistries and patterns, through the skin lipid bilayer was studied using constrained and unconstrained molecular dynamics simulations. Interestingly, the permeation mechanism of nanoparticles having the same surface chemistry but different patterns was found to be completely different. Nanoparticles (NPs) were screened based on the free energy of permeation through the skin lipid bilayer. The behavior of the screened NPs was further validated with unconstrained simulations using the skin lipid bilayer. Nanoparticles thus screened through both of the techniques were further used for the co-delivery of a model protein into the skin lipid bilayer. It was observed that the nanoparticles having a 2 : 1 homogeneous ratio of hydrophobic to hydrophilic regions were the most promising in transdermal delivery of proteins. The obtained results are in line with the results of recent permeation experiments on cell and plasma membrane. Our study could help in in-silico design of nanoparticles for delivery of actives through skin. These in-silico experiments thus could help speed up the development process by guiding formulation chemists.