Bilayer structure and lipid dynamics in a model stratum corneum with oleic acid

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.

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.

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.

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.

Characterisation of the Molecular Mechanism of Permeation of the Prodrug Me-5ALA across the Human Stratum Corneum Using Molecular Dynamics Simulations

The barrier imposed by the outer layer of the skin, the stratum corneum, creates an almost impermeable environment for exogenous substances. Few lipophilic drugs with low molecular mass can passively diffuse through this layer, highlighting the need to develop methods to enable the delivery of more drugs via the transdermal route. The prodrug approach involves modifying the structure of a drug molecule to enhance its permeability across the skin, but it is often difficult to predict how exactly changes in chemical structure affect permeation. This study uses molecular dynamics simulations to predict permeability values and adequately characterise the molecular mechanism of permeation of the prodrugs Me-5ALA and its parent compound 5ALA across a molecular model of the lipid bilayers of the human stratum corneum. The influence of increased hydrophobicity in Me-5ALA on its permeation revealed a reduction in hydrogen bonding capability that enables it to interact more favourably with the hydrophobic region of the bilayer and diffuse at a faster rate with less resistance, thus making it a better permeant compared to its more hydrophilic parent compound. This molecular simulation approach offers a promising route for the rational design of drug molecules that can permeate effectively across the stratum corneum.

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.

Comparison Between Franz Diffusion Cell and a novel Micro-physiological System for In Vitro Penetration Assay Using Different Skin Models

In vitro diffusive models are an important tool to screen the penetration ability of active ingredients in various formulations. A reliable assessment of skin penetration enhancing properties, mechanism of action of carrier systems, and an estimation of a bioavailability are essential for transdermal delivery. Given the importance of testing the penetration kinetics of different compounds across the skin barrier, several in vitro models have been developedThe aim of this study was to compare the Franz Diffusion Cell (FDC) with a novel fluid-dynamic platform (MIVO) by evaluating penetration ability of caffeine, a widely used reference substance, and LIP1, a testing molecule having the same molecular weight but a different lipophilicity in the two diffusion chamber systems. A 0.7% caffeine or LIP1 formulation in either water or propylene glycol (PG) containing oleic acid (OA) was topically applied on the Strat-M® membrane or pig ear skin, according to the infinite-dose experimental condition (780 ul/cm²). The profile of the penetration kinetics was determined by quantify the amount of molecule absorbed at different time-points (1, 2, 4, 6, 8 hours), by means of HPLC analysis. Both diffusive systems show a similar trend for caffeine and LIP1 penetration kinetics. The Strat-M® skin model shows a lower barrier function than the pig skin biopsies, whereby the PGOA vehicle exhibits a higher penetration, enhancing the effect for both diffusive chambers and skin surrogates. Most interestingly, MIVO diffusive system better predicts the lipophilic molecules (i.e. LIP1) permeation through highly physiological fluid flows resembled below the skin 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.

Topical drug delivery: History, percutaneous absorption, and product development

Topical products, widely used to manage skin conditions, have evolved from simple potions to sophisticated delivery systems. Their development has been facilitated by advances in percutaneous absorption and product design based on an increasingly mechanistic understanding of drug-product-skin interactions, associated experiments, and a quality-by-design framework. Topical drug delivery involves drug transport from a product on the skin to a local target site and then clearance by diffusion, metabolism, and the dermal circulation to the rest of the body and deeper tissues. Insights have been provided by Quantitative Structure Permeability Relationships (QSPR), molecular dynamics simulations, and dermal Physiologically Based PharmacoKinetics (PBPK). Currently, generic product equivalents of reference-listed products dominate the topical delivery market. There is an increasing regulatory interest in understanding topical product delivery behavior under 'in use' conditions and predicting in vivo response for population variations in skin barrier function and response using in silico and in vitro findings.

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.

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.

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.

Permeation of Biopolymers Across the Cell Membrane: A Computational Comparative Study on Polylactic Acid and Polyhydroxyalkanoate

Polymeric nanoparticles, which by virtue of their size (1-1000 nm) are able to penetrate even into cells, are attracting increasing interest in the emerging field of nanomedicine, as devices for, e.g., drugs or vaccines delivery. Because of the involved dimensional scale in the nanoparticle/cell membrane interactions, modeling approaches at molecular level are the natural choice in order to understand the impact of nanoparticle formulation on cellular uptake mechanisms. In this work, the passive permeation across cell membrane of oligomers made of two employed polymers in the biomedical field [poly-D,L-lactic acid (PDLA) and poly(3-hydroxydecanoate) (P3HD)] is investigated at fundamental atomic scale through molecular dynamics simulations. The free energy profile related to membrane crossing is computed adopting umbrella sampling. Passive permeation is also investigated using a coarse-grained model with MARTINI force field, adopting well-tempered metadynamics. Simulation results showed that P3HD permeation is favored with respect to PDLA by virtue of its higher hydrophobicity. The free energy profiles obtained at full atomistic and coarse-grained scale are in good agreement each for P3HD, while only a qualitative agreement was obtained for PDLA. Results suggest that a reparameterization of non-bonded interactions of the adopted MARTINI beads for the oligomer is needed in order to obtain a better agreement with more accurate simulations at atomic scale.

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.

l-Ascorbic acid alkyl esters action on stratum corneum model membranes: An insight into the mechanism for enhanced skin permeation

L-ascorbic acid alkyl esters (ASCn) are lipophilic forms of vitamin C, which act as skin permeation enhancers. We investigated the physical changes induced by incorporating ASCn into stratum corneum (SC) lipid membranes and correlated this with the mechanism proposed in the literature for skin permeation enhancement phenomena. We used lipid monolayers to explore the 2D structure and elasticity of the lipid-enhancer systems. As a comparison, the classic permeation enhancer, oleic acid (OA) and the non-enhancer analogue stearic acid (SA) were analysed. The incorporation of ASCn or OA into SC membranes resulted in more liquid-like films, with a dose-dependent lowering of the compressibility modulus. Brewster angle microscopy (BAM) evidenced partial miscibility of the enhancer with SC lipid components, stabilising the liquid-expanded phase. At the nanoscale, AFM showed that SC lipids form heterogeneous membranes, which underwent structural alterations after incorporating ASCn and fatty acids, such as SA and OA. The lower, cholesterol-enriched phase appears to concentrate the enhancers, whilst the higher ceramide-enriched phase concentrated the non-enhancer SA. Our results and previously reported pieces of evidence indicate a strong pattern in which the rheological properties of SC lipid films are determinant for skin permeation phenomena.

Derailed Ceramide Metabolism in Atopic Dermatitis (AD): A Causal Starting Point for a Personalized (Basic) Therapy

Active rebuilding, stabilizing, and maintaining the lipid barrier of the skin is an encouraging disease management and care concept for dry skin, atopic dermatitis (eczema, neurodermatitis), and psoriasis. For decades, corticosteroids have been the mainstay of topical therapy for atopic dermatitis; however, innovations within the scope of basic therapy are rare. In (extremely) dry, irritated, or inflammatory skin, as well as in lesions, an altered (sphingo)lipid profile is present. Recovery of a balanced (sphingo)lipid profile is a promising target for topical and personalized treatment and prophylaxis. New approaches for adults and small children are still lacking. With an ingenious combination of commonly used active ingredients, it is possible to restore and reinforce the dermal lipid barrier and maintain refractivity. Lysosomes and ceramide de novo synthesis play a key role in attenuation of the dermal lipid barrier. Linoleic acid in combination with amitriptyline in topical medication offers the possibility to relieve patients affected by dry and itchy skin, mild to moderate atopic dermatitis lesions, and eczemas without the commonly occurring serious adverse effects of topical corticosteroids or systemic antibody administration.

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.

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.

Clearly Detectable, Kinetically Restricted Solid-Solid Phase Transition in cis-Ceramide Monolayers

Sphingosine [(2 S,3 R,4 E)-2-amino-4-octadecene-1,3-diol] is the most common sphingoid base in mammals. Ceramides are N-acyl sphingosines. Numerous small variations on this canonical structure are known, including the 1-deoxy, the 4,5-dihydro, and many others. However, whenever there is a Δ4 double bond, it adopts the trans (or E) configuration. We synthesized a ceramide containing 4 Z-sphingosine and palmitic acid ( cis-pCer) and studied its behavior in the form of monolayers extended on an air-water interface. cis-pCer acted very differently from the trans isomer in that, upon lateral compression of the monolayer, a solid-solid transition was clearly observed at a mean molecular area ≤44 Ų·molecule^-1, whose characteristics depended on the rate of compression. The solid-solid transition, as well as states of domain coexistence, could be imaged by atomic force microscopy and by Brewster-angle microscopy. Atomistic molecular dynamics simulations provided results compatible with the experimentally observed differences between the cis and trans isomers. The data can help in the exploration of other solid-solid transitions in lipids, both in vitro and in vivo, that have gone up to now undetected because of their less obvious change in surface properties along the transition, as compared to cis-pCer.

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.

Caco-2 Cell Conditions Enabling Studies of Drug Absorption from Digestible Lipid-Based Formulations

Purpose

To identify conditions allowing the use of cell-based models for studies of drug absorption during in vitro lipolysis of lipid-based formulations (LBFs).

Methods

Caco-2 was selected as the cell-based model system. Monolayer integrity was evaluated by measuring mannitol permeability after incubating Caco-2 cells in the presence of components available during lipolysis. Pure excipients and formulations representing the lipid formulation classification system (LFCS) were evaluated before and after digestion. Porcine mucin was evaluated for its capacity to protect the cell monolayer.

Results

Most undigested formulations were compatible with the cells (II-LC, IIIB-LC, and IV) although some needed mucin to protect against damaging effects (II-MC, IIIB-MC, I-LC, and IIIA-LC). The pancreatic extract commonly used in digestion studies was incompatible with the cells but the Caco-2 monolayers could withstand immobilized recombinant lipase. Upon digestion, long chain formulations caused more damage to Caco-2 cells than their undigested counterparts whereas medium chain formulations showed better tolerability after digestion.

Conclusions

Most LBFs and components thereof (undigested and digested) are compatible with Caco-2 cells. Pancreatic enzyme is not tolerated by the cells but immobilized lipase can be used in combination with the cell monolayer. Mucin is beneficial for critical formulations and digestion products.

Electronic supplementary material

The online version of this article (10.1007/s11095-017-2327-8) contains supplementary material, which is available to authorized users.

Effects of carotenoids on lipid bilayers

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. However, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity.

Permeation pathways through lateral domains in model membranes of skin lipids

Lateral organisation of skin lipids in membranes produces regions with different permeability; water permeation is favoured through cholesterol-rich regions.

Getting Drugs Across Biological Barriers

The delivery of drugs to a target site frequently involves crossing biological barriers. The degree and nature of the impediment to flux, as well as the potential approaches to overcoming it, depend on the tissue, the drug, and numerous other factors. Here an overview of approaches that have been taken to crossing biological barriers is presented, with special attention to transdermal drug delivery. Technology and knowledge pertaining to addressing these issues in a variety of organs could have a significant clinical impact.

In vivo confocal Raman spectroscopy and molecular dynamics analysis of penetration of retinyl acetate into stratum corneum

OBJECTIVE

The purpose of this study is to elucidate the behavior of retinyl acetate in penetrating human skin without the presence of enhancers by using confocal Raman spectroscopy and molecular dynamics simulation.

METHODS

In this study, in vivo confocal Raman spectroscopy was combined with molecular dynamics simulation to investigate the transdermal permeation of the aqueous suspension of retinyl acetate.

RESULTS

Permeation was measured after 30min, and retinyl acetate was found up to 20μm deep inside the stratum corneum. The delivery of retinyl acetate inside a skin membrane model was studied by molecular dynamics. The membrane model that was used represented normal young skin containing a lipid bilayer with 25% ceramide, 36% fatty acid, 30% cholesterol, and 6% cholesterol sulfate.

CONCLUSION

Spectroscopy data indicate that retinyl acetate permeates into the stratum corneum. Molecular dynamics data showed that retinyl acetate permeates in the membrane model and that their final location is deep inside the lipid bilayer. We showed, for the first time, a correlation between Raman permeation data and computational data.

Protein-dependent Membrane Interaction of A Partially Disordered Protein Complex with Oleic Acid: Implications for Cancer Lipidomics

Bovine α-lactalbumin (BLA) forms cytotoxic complexes with oleic acid (OA) that perturbs tumor cell membranes, but molecular determinants of these membrane-interactions remain poorly understood. Here, we aim to obtain molecular insights into the interaction of BLA/BLA-OA complex with model membranes. We characterized the folding state of BLA-OA complex using tryptophan fluorescence and resolved residue-specific interactions of BLA with OA using molecular dynamics simulation. We integrated membrane-binding data using a voltage-sensitive probe and molecular dynamics (MD) to demonstrate the preferential interaction of the BLA-OA complex with negatively charged membranes. We identified amino acid residues of BLA and BLA-OA complex as determinants of these membrane interactions using MD, functionally corroborated by uptake of the corresponding α-LA peptides across tumor cell membranes. The results suggest that the α-LA component of these cytotoxic complexes confers specificity for tumor cell membranes through protein interactions that are maintained even in the lipid complex, in the presence of OA.

A Coarse-Grained Model of Stratum Corneum Lipids: Free Fatty Acids and Ceramide NS

Ceramide (CER)-based biological membranes are used both experimentally and in simulations as simplified model systems of the skin barrier. Molecular dynamics studies have generally focused on simulating preassembled structures using atomistically detailed models of CERs, which limit the system sizes and time scales that can practically be probed, rendering them ineffective for studying particular phenomena, including self-assembly into bilayer and lamellar superstructures. Here, we report on the development of a coarse-grained (CG) model for CER NS, the most abundant CER in human stratum corneum. Multistate iterative Boltzmann inversion is used to derive the intermolecular pair potentials, resulting in a force field that is applicable over a range of state points and suitable for studying ceramide self-assembly. The chosen CG mapping, which includes explicit interaction sites for hydroxyl groups, captures the directional nature of hydrogen bonding and allows for accurate predictions of several key structural properties of CER NS bilayers. Simulated wetting experiments allow the hydrophobicity of CG beads to be accurately tuned to match atomistic wetting behavior, which affects the whole system, since inaccurate hydrophobic character is found to unphysically alter the lipid packing in hydrated lamellar states. We find that CER NS can self-assemble into multilamellar structures, enabling the study of lipid systems more representative of the multilamellar lipid structures present in the skin barrier. The coarse-grained force field derived herein represents an important step in using molecular dynamics to study the human skin barrier, which gives a resolution not available through experiment alone.

Molecular Dynamics Simulation Study of Permeation of Molecules through Skin Lipid Bilayer

Stratum Corneum (SC), the outermost layer of skin, is mainly responsible for skin's barrier function. The complex lipid matrix of SC determines these barrier properties. In this study, the lipid matrix is modeled as an equimolar mixture of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA). The permeation of water, oxygen, ethanol, acetic acid, urea, butanol, benzene, dimethyl sulfoxide (DMSO), toluene, phenol, styrene, and ethylbenzene across this layer is studied using a constrained MD simulations technique. Several long constrained simulations are performed at a skin temperature of 310 K under NPT conditions. The free energy profiles and diffusion coefficients along the bilayer normal have been calculated for each molecule. Permeability coefficients are also calculated and compared with experimental data. The main resistance for the permeation of hydrophilic and hydrophobic permeants has been found to be in the interior of the lipid bilayer and near the lipid-water interface, respectively. The obtained permeability is found to be a few orders of magnitude higher than experimental values for hydrophilic molecules while for hydrophobic molecules more discrepancy was observed. Overall, the qualitative ranking is consistent with the experiments.

The physics of stratum corneum lipid membranes

The stratum corneum (SC), the outermost layer of skin, comprises rigid corneocytes (keratin-filled dead cells) in a specialized lipid matrix. The continuous lipid matrix provides the main barrier against uncontrolled water loss and invasion of external pathogens. Unlike all other biological lipid membranes (such as intracellular organelles and plasma membranes), molecules in the SC lipid matrix show small hydrophilic groups and large variability in the length of the alkyl tails and in the numbers and positions of groups that are capable of forming hydrogen bonds. Molecular simulations provide a route for systematically probing the effects of each of these differences separately. In this article, we present the results from atomistic molecular dynamics of selected lipid bilayers and multi-layers to probe the effect of these polydispersities. We address the nature of the tail packing in the gel-like phase, the hydrogen bond network among head groups, the bending moduli expected for leaflets comprising SC lipids and the conformation of very long ceramide lipids in multi-bilayer lipid assemblies.This article is part of the themed issue 'Soft interfacial materials: from fundamentals to formulation'.

Effect of monoglycerides and fatty acids on a ceramide bilayer

Monoglycerides and unsaturated fatty acids, naturally present in trace amounts in the stratum corneum (top layer of skin) lipid matrix, are commonly used in pharmaceutical, cosmetic and health care formulations. However, a detailed molecular understanding of how the oil additives get incorporated into the skin lipids from topical application and, once incorporated, how they affect the properties and integrity of the lipid matrix remains unexplored. Using ceramide 2 bilayers as skin lipid surrogates, we use a series of molecular dynamics simulations with six different natural oil ingredients at multiple concentrations to investigate the effect of the oils on the properties and stability of the bilayers. The six oils: monoolein, monostearin, monoelaidin, oleic acid, stearic acid and linoleic acid - all having the same length of the alkyl chain, C18, but a varying degree of saturation, allow us to systematically address the effect of unsaturation in the additives. Our results show that at low oil concentration (∼5%) the mixed bilayers containing any of the oils and ceramide 2 (CER2) become more rigid than pure CER2 bilayers due to more efficient lipid packing. Better packing also results in the formation of larger numbers of hydrogen bonds between the lipids, which occurs at the expense of the hydrogen bonds between lipids and water. The mixed bilayers with saturated or trans-unsaturated oils remain stable over the whole range of oil concentration. In contrast, the presence of the oils with at least one cis-double bond leads to bilayer instability and complete loss of bilayer structure at the oil content of about 50-65%. Two cis-double bonds in the lipid tail induce bilayer disruption at even lower concentration (∼30%). The mixed bilayers remain in the gel phase (without melting to a fluid phase) until the phase transition to a non-bilayer phase occurs. We also demonstrate that the stability of the bilayer strongly correlates with the order parameter of the lipid tails.

Chemical penetration enhancers in stratum corneum - Relation between molecular effects and barrier function

Skin is attractive for drug therapy because it offers an easily accessible route without first-pass metabolism. Transdermal drug delivery is also associated with high patient compliance and through the site of application, the drug delivery can be locally directed. However, to succeed with transdermal drug delivery it is often required to overcome the low permeability of the upper layer of the skin, the stratum corneum (SC). One common strategy is to employ so-called penetration enhancers that supposedly act to increase the drug passage across SC. Still, there is a lack of understanding of the molecular effects of so-called penetration enhancers on the skin barrier membrane, the SC. In this study, we provide a molecular characterization of how different classes of compounds, suggested as penetration enhancers, influence lipid and protein components in SC. The compounds investigated include monoterpenes, fatty acids, osmolytes, surfactant, and Azone. We employ natural abundance (13)C polarization transfer solid-state nuclear magnetic resonance (NMR) on intact porcine SC. With this method it is possible to detect small changes in the mobility of the minor fluid lipid and protein SC components, and simultaneously obtain information on the major fraction of solid SC components. The balance between fluid and solid components in the SC is essential to determine macroscopic material properties of the SC, including barrier and mechanical properties. We study SC at different hydration levels corresponding to SC in ambient air and under occlusion. The NMR studies are complemented with diffusion cell experiments that provide quantitative data on skin permeability when treated with different compounds. By correlating the effects on SC molecular components and SC barrier function, we aim at deepened understanding of diffusional transport in SC, and how this can be controlled, which can be utilized for optimal design of transdermal drug delivery formulations.

Functionalized lipids and surfactants for specific applications

Synthetic lipids and surfactants that do not exist in biological systems have been used for the last few decades in both basic and applied science. The most notable applications for synthetic lipids and surfactants are drug delivery, gene transfection, as reporting molecules, and as support for structural lipid biology. In this review, we describe the potential of the synergistic combination of computational and experimental methodologies to study the behavior of synthetic lipids and surfactants embedded in lipid membranes and liposomes. We focused on select cases in which molecular dynamics simulations were used to complement experimental studies aiming to understand the structure and properties of new compounds at the atomistic level. We also describe cases in which molecular dynamics simulations were used to design new synthetic lipids and surfactants, as well as emerging fields for the application of these compounds. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Coarse-grain simulations of skin ceramide NS with newly derived parameters clarify structure of melted phase

Ceramides are lipids that are involved in numerous biologically important structures (e.g., the stratum corneum and ceramide-rich platforms) and processes (e.g., signal transduction and membrane fusion), but their behavior is not fully understood. We report coarse-grain force field parameters for N-lignocerylsphingosine (ceramide NS, also known as ceramide 2) that are consistent with the Martini force field. These parameters were optimized for simulations in the gel phase and validated against atomistic simulations. Coarse-grained simulations with our parameters provide areas per lipid, membrane thicknesses, and electron density profiles that are in good agreement with atomistic simulations. Properties of the simulated membranes are compared with available experimental data. The obtained parameters were used to model the phase behavior of ceramide NS as a function of temperature and hydration. At low water content and above the main phase transition temperature, the bilayer melts into an irregular phase, which may correspond to the unstructured melted-chain phase observed in X-ray diffraction experiments. The developed parameters also reproduce the extended conformation of ceramide, which may occur in the stratum corneum. The parameters presented herein will facilitate studies on important complex functional structures such as the uppermost layer of the skin and ceramide-rich platforms in phospholipid membranes.

The effects of solvent composition on the affinity of a peptide towards hair keratin: experimental and molecular dynamics data

Molecular dynamics simulations with a developed hair protofibril model demonstrated the ability to improve peptide uptake by hair shafts.

Tristearin bilayers: structure of the aqueous interface and stability in the presence of surfactants

Molecular dynamics simulations predict that sodium dodecylbenzene sulphonate surfactant molecules embed themselves in a tristearin bilayer, packing commensurate with the hexagonally packed lattice formed by the acyl tails of tristearin.

Rationalization of reduced penetration of drugs through ceramide gel phase membrane

Since computing resources have advanced enough to allow routine molecular simulation studies of drug molecules interacting with biologically relevant membranes, a considerable amount of work has been carried out with fluid phospholipid systems. However, there is very little work in the literature on drug interactions with gel phase lipids. This poses a significant limitation for understanding permeation through the stratum corneum where the primary pathway is expected to be through a highly ordered lipid matrix. To address this point, we analyzed the interactions of p-aminobenzoic acid (PABA) and its ethyl (benzocaine) and butyl (butamben) esters with two membrane bilayers, which differ in their fluidity at ambient conditions. We considered a dioleoylphosphatidylcholine (DOPC) bilayer in a fluid state and a ceramide 2 (CER2, ceramide NS) bilayer in a gel phase. We carried out unbiased (100 ns long) and biased z-constraint molecular dynamics simulations and calculated the free energy profiles of all molecules along the bilayer normal. The free energy profiles converged significantly slower for the gel phase. While the compounds have comparable affinities for both membranes, they exhibit penetration barriers almost 3 times higher in the gel phase CER2 bilayer. This elevated barrier and slower diffusion in the CER2 bilayer, which are caused by the high ordering of CER2 lipid chains, explain the low permeability of the gel phase membranes. We also compared the free energy profiles from MD simulations with those obtained from COSMOmic. This method provided the same trends in behavior for the guest molecules in both bilayers; however, the penetration barriers calculated by COSMOmic did not differ between membranes. In conclusion, we show how membrane fluid properties affect the interaction of drug-like molecules with membranes.

Molecular interactions of plant oil components with stratum corneum lipids correlate with clinical measures of skin barrier function

Plant-derived oils consisting of triglycerides and small amounts of free fatty acids (FFAs) are commonly used in skincare regimens. FFAs are known to disrupt skin barrier function. The objective of this study was to mechanistically study the effects of FFAs, triglycerides and their mixtures on skin barrier function. The effects of oleic acid (OA), glyceryl trioleate (GT) and OA/GT mixtures on skin barrier were assessed in vivo through measurement of transepidermal water loss (TEWL) and fluorescein dye penetration before and after a single application. OA's effects on stratum corneum (SC) lipid order in vivo were measured with infrared spectroscopy through application of perdeuterated OA (OA-d34 ). Studies of the interaction of OA and GT with skin lipids included imaging the distribution of OA-d34 and GT ex vivo with IR microspectroscopy and thermodynamic analysis of mixtures in aqueous monolayers. The oil mixtures increased both TEWL and fluorescein penetration 24 h after a single application in an OA dose-dependent manner, with the highest increase from treatment with pure OA. OA-d34 penetrated into skin and disordered SC lipids. Furthermore, the ex vivo IR imaging studies showed that OA-d34 permeated to the dermal/epidermal junction while GT remained in the SC. The monolayer experiments showed preferential interspecies interactions between OA and SC lipids, while the mixing between GT and SC lipids was not thermodynamically preferred. The FFA component of plant oils may disrupt skin barrier function. The affinity between plant oil components and SC lipids likely determines the extent of their penetration and clinically measurable effects on skin barrier functions.

Use of Awamori-pressed Lees and Tofu Lees as Feed Ingredients for Growing Male Goats

Awamori is produced by fermenting steamed indica rice. Awamori-pressed lees is a by-product of the Awamori production process. Tofu lees is a by-product of the Tofu production process. Research was conducted to test if dried Awamori-pressed lees and Tofu lees can be used as a mixed feed ingredient for raising male goats. Eighteen male kids were divided into three groups of six animals (control feed group (CFG), Awamori-pressed lees mixed feed group (AMFG), Tofu lees mixed feed group (TMFG)). The CFG used feed containing 20% soybean meal as the main protein source, while the AMFG and TMFG used feed mixed with 20% dried Awamori-pressed lees or dried Tofu lees. The groups were fed mixed feed (volume to provide 100 g/d increase in body weight) and alfalfa hay cubes (2.0 kg/d) twice a day (10:00, 16:00). Klein grass hay and water was given ad libitum. Hay intake was measured at 10:00 and 16:00. Body weight and size measurements were taken once a month. At the end of the experiment, a blood sample was drawn from the jugular vein of each animal and the carcass characteristics, the physical and chemical characteristics of loin were analyzed. DCP and TDN intakes in AMFG and TMFG showed no significant difference to the CFG. Cumulative measurements of growth in body weight and size over the 10 mo period in the AMFG and TMFG were similar to the CFG. Blood parameter values were similar to those in normal goats. Dressing carcass weight and percentages, and total weight of meat in the AMFG were similar to that in the CFG, but smaller in the TMFG. The compressed meat juice ratio was higher in both the TMFG and AMFG than the CFG. While the fat in corn, Awamori-pressed lees, and Tofu lees contains more than 50% linoleic acid, the loin fat in both the AMFG and TMFG was very low in linoleic acid due to the increase in the content of oleic acid, stearic acid, and palmitic acid. This indicates that feeding on AMF and TMF does not inhibit hydrogenation by ruminal microorganisms. As in the CFG, the total essential and non-essential amino acids in the loin of the AMFG and TMFG were well balanced. Compared to the CFG, the AMFG and TMFG were high in taurine and carnosine. The results indicate dried Awamori-pressed lees and Tofu lees can be used as a feed ingredient for raising male goats.

Localization of cholesterol and fatty acid in a model lipid membrane: a neutron diffraction approach

The intercellular lipid matrix of the skin's stratum corneum serves to protect the body against desiccation and simultaneously limits the passage of drugs and other xenobiotics into the body. The matrix is made up of ceramides, free fatty acids, and cholesterol, which are organized as two coexisting crystalline lamellar phases. In studies reported here, we sought to use the technique of neutron diffraction, together with the device of isotopic (H/D) substitution, to determine the molecular architecture of the lamellar phase having a repeat distance of 53.9 ± 0.3 Å. Using hydrogenous samples as well as samples incorporating perdeuterated (C24:0) fatty acids and selectively deuterated cholesterol, the diffraction data obtained were used to construct neutron scattering length density profiles. By this means, the locations within the unit cell were determined for the cholesterol and fatty acids. The cholesterol headgroup was found to lie slightly inward from the unit cell boundary and the tail of the molecule located 6.2 ± 0.2 Å from the unit cell center. The fatty acid headgroups were located at the unit cell boundary with their acyl chains straddling the unit cell center. Based on these results, a molecular model is proposed for the arrangement of the lipids within the unit cell.