Molecular Dynamics Simulation of Skin Lipids: Effect of Ceramide Chain Lengths on Bilayer Properties

Dynamic Simulations of Interaction of the PEG-DPPE Micelle-Encapsulated Short-Chain Ceramides with the Raft-Included Membrane

The lipid raft subdomains in cancer cell membranes play a key role in signal transduction, biomolecule recruitment, and drug transmembrane transport. Augmented membrane rigidity due to the formation of a lipid raft is unfavorable for the entry of drugs, a limiting factor in clinical oncology. The short-chain ceramide (CER) has been reported to promote drug entry into membranes and disrupt lipid raft formation, but the underlying mechanism is not well understood. We recently explored the carrier-membrane fusion dynamics of PEG-DPPE micelles in delivering doxorubicin (DOX). Based on the phase-segregated membrane model composed of DPPC/DIPC/CHOL/GM1/PIP2, we aim to explore the dynamic mechanism of the PEG-DPPE micelle-encapsulating DOXs in association with the raft-included cell membrane modulated by C8 acyl tail CERs. The results show that the lipid raft remains integrated and DOX-resistant subjected to free DOXs and the micelle-encapsulating ones. Addition of CERs disorganizes the lipid raft by pushing CHOL aside from DPPC. It subsequently allows for a good permeability for PEG-DPPE micelle-encapsulated DOXs, which penetrate deeper as CER concentration increases. GM1 is significant in guiding drugs' redistributing between bilayer phases, and the anionic PIP2 further helps DOXs attain the inner bilayer surface. These results elaborate on the perturbing effect of CERs on lipid raft stability, which provides a new comprehensive approach for further design of drug delivery systems.

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.

Self-Assembly of Lipid Molecules under Shear Flows: A Dissipative Particle Dynamics Simulation Study

The self-assembly of lipid molecules in aqueous solution under shear flows was investigated using the dissipative particle dynamics simulation method. Three cases were considered: zero shear flow, weak shear flow and strong shear flow. Various self-assembled structures, such as double layers, perforated double layers, hierarchical discs, micelles, and vesicles, were observed. The self-assembly behavior was investigated in equilibrium by constructing phase diagrams based on chain lengths. Results showed the remarkable influence of chain length, shear flow and solution concentration on the self-assembly process. Furthermore, the self-assembly behavior of lipid molecules was analyzed using the system energy, particle number and shape factor during the dynamic processes, where the self-assembly pathways were observed and analyzed for the typical structures. The results enhance our understanding of biomacromolecule self-assembly in a solution and hold the potential for applications in biomedicine.

Investigation of β-caryophyllene as terpene penetration enhancer: Role of stratum corneum retention

Terpenes are usually used as penetration enhancers (PE) for transdermal drug delivery (TDD) of various molecules. However, TDD of hydrophilic macromolecules is becoming an urgent challenge due to their potent activities. The aim of this study was to investigate the potential application of β-caryophyllene (β-CP), a sequiterpene, as PE for TDD of hydrophilic macromolecules for the first time. Commonly used PEs, namely azone and 1,8-cineole (1,8-CN), were applied as controls. Transepidermal water loss (TEWL) analysis revealed that the reduction of skin barrier function caused by β-CP was reversible. Transdermal experiments showed that when skin was treated with β-CP or azone, there was a significant permeation-enhancing effect on fluorescein isothiocyanate (FITC) and FITC-dextran with different molecular weight (MW) of 4k or 10k. CLSM analysis confirmed that β-CP and azone can facilitate the penetration of FD-4k through epidermis and dermis. However, the cytotoxicity of azone against epidermal keratinocytes was significantly higher than β-CP and 1,8-CN. Additionally, application of β-CP and 1,8-CN didn't increase erythema index (EI) but the EI values of azone group increased significantly and irreversibly, indicating the high biocompatibility of the natural terpenes. β-CP had better permeation-enhancing effect and higher stratum corneum (SC) retention than 1,8-CN due to its increased carbon chain length and lipophilicity, as further demonstrated by molecular dynamics (MD) simulation studies. Skin electrical resistance (SER) and attenuated total reflection fourier transform infrared spectroscopy (ATR-FTIR) studies revealed a significant interfering effect of β-CP on SC lipids. Taken together, β-CP exhibited significant penetration enhancement of hydrophilic macromolecules due to its SC retention and SC lipid fluidization ability.

Microphase transitions of Langmuir-Blodgett lipid-assembled monolayers with new types of ceramides, ultra-long-chain ceramide and 1-O-acylceramide

HYPOTHESIS

Intercellular lipid lamellae, consisting of ceramide, cholesterol, and free fatty acids, are the primary pathways for substances in the stratum corneum (SC). The microphase transition of lipid-assembled monolayers (LAMs), mimicking an initial layer of the SC, would be affected by new types of ceramides: ceramide with ultra-long chain (CULC) and 1-O-acylceramide (CENP) with three chains in different direction.

EXPERIMENTS

The LAMs were fabricated with varying the mixing ratio of CULC (or CENP) against base ceramide via Langmuir-Blodgett assembly. Surface pressure-area isotherms and elastic modulus-surface pressure plots were obtained to characterize π-dependent microphase transitions. The surface morphology of LAMs was observed by atomic force microscopy.

FINDINGS

The CULCs favored lateral lipid packing, and the CENPs hindered the lateral lipid packing by lying alignment, which was due to their different molecular structures and conformations. The sporadic clusters and empty spaces in the LAMs with CULC were presumably due to the short-range interactions and self-entanglements of ultra-long alkyl chains following the freely jointed chain model, respectively, which was not noticeably observed in the neat LAM films and the LAM films with CENP. The addition of surfactants disrupted the lateral packing of lipids, thus weakening the LAM elasticity. These findings allowed us to understand the role of CULC and CENP in the lipid assemblies and microphase transition behaviors in an initial layer of SC.

Conformations of Three Types of Ultra-Long-Chain Fatty Acids in Multicomponent Lipid Bilayers

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in certain types of tissues, but their biological roles remain unknown. Here, we report how the conformation of ULCFAs depends on the length and unsaturated-bond ratio of the ultra-long chains and the composition of the host bilayer membrane using molecular dynamics simulations. The ultra-long chain of ULCFAs flips between the two leaflets and fluctuates among three conformations: elongated, L-shaped, and turned. Furthermore, we found that the saturated ultra-long chain exhibited an elongated conformation more frequently than the unsaturated chain. In addition, the truncation of the ultra-long chain at C26 had little effect on the remaining ULCFAs. ULCFAs respond to lipid-density differences in the two leaflets, and the ratio of the elongated and turned conformations changed to reduce this difference. However, in cholesterol-containing membranes, ULCFAs exhibit no density difference after the flip-flop of cholesterol removes the difference.

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.

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.

Asymmetric Lipid Membranes under Shear Flows: A Dissipative Particle Dynamics Study

We investigate the phase behavior of the asymmetric lipid membranes under shear flows, using the dissipative particle dynamics simulation. Two cases, the weak and strong shear flows, are considered for the asymmetric lipid microstructures. Three typical asymmetric structures, the membranes, tubes, and vesicle, are included in the phase diagrams, where the effect of two different types of lipid chain length on the formation of asymmetric membranes is evaluated. The dynamic processes are demonstrated for the asymmetric membranes by calculating the average radius of gyration and shape factor. The result indicates that different shear flows will affect the shape of the second type of lipid molecules; the shape of the first type of lipid molecules is more stable than that of the second type of lipid molecules. The mechanical properties are investigated for the asymmetric membranes by analyzing the interface tension. The results reveal an absolute pressure at the junctions of different types of particles under the weak shear flow; the other positions are almost in a state of no pressure; there is almost no pressure inside the asymmetric lipid membrane structure under the strong shear flow. The findings will help us to understand the potential applications of asymmetric lipid microstructures in the biological and medical fields.

A novel mineralocorticoid receptor antagonist, 7,3',4'-trihydroxyisoflavone improves skin barrier function impaired by endogenous or exogenous glucocorticoids

Excess glucocorticoids (GCs) with either endogenous or exogenous origins deteriorate skin barrier function. GCs bind to mineralocorticoid and GC receptors (MRs and GRs) in normal human epidermal keratinocytes (NHEKs). Inappropriate MR activation by GCs mediates various GC-induced cutaneous adverse events. We examined whether MR antagonists can ameliorate GC-mediated skin barrier dysfunction in NHEKs, reconstructed human epidermis (RHE), and subjects under psychological stress (PS). In a preliminary clinical investigation, topical MR antagonists improved skin barrier function in topical GC-treated subjects. In NHEKs, cortisol induced nuclear translocation of GR and MR, and GR and MR antagonists inhibited cortisol-induced reductions of keratinocyte differentiation. We identified 7,3',4'-trihydroxyisoflavone (7,3',4'-THIF) as a novel compound that inhibits MR transcriptional activity by screening 30 cosmetic compounds. 7,3',4'-THIF ameliorated the cortisol effect which decreases keratinocyte differentiation in NHEKs and RHE. In a clinical study on PS subjects, 7,3',4'-THIF (0.1%)-containing cream improved skin barrier function, including skin surface pH, barrier recovery rate, and stratum corneum lipids. In conclusion, skin barrier dysfunction owing to excess GC is mediated by MR and GR; thus, it could be prevented by treatment with MR antagonists. Therefore, topical MR antagonists are a promising therapeutic option for skin barrier dysfunction after topical GC treatment or PS.

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

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

Conformation of ultra-long-chain fatty acid in lipid bilayer: Molecular dynamics study

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in the restricted tissues such as retina, testis, and skin. The conformation of a single ULCFA, in which the sn-1 unsaturated chain has 32 carbons, in three types of phospholipid bilayers is studied by molecular dynamics simulations. It is found that the ultra-long tail of the ULCFA flips between two leaflets and fluctuates among an elongation into the opposite leaflet, lies between two leaflets, and turns back. As the number ratio of lipids in the opposite leaflet increases, the ratio of the elongated shape linearly decreases in all three cases. Thus, ULCFAs can sense the density differences between the two leaflets and respond to these changes.

Development and application of coarse-grained MARTINI model of skin lipid ceramide [AP]

Stratum corneum (SC), the outermost layer of the skin, contains large variety of lipids, endowing them with the amphiphilic properties, needed to fulfil their key role in skin's barrier function. The individual role of lipid types in the barrier function is difficult to understand due to the immense heterogeneity and complexity of the lipid's organization within the SC. The lipid organization is being explored using both computational (molecular dynamics simulations) and experimental (neutron diffraction) techniques. Even though atomistic simulations provide unprecedented atomic level details, the major limitation is time and length scale that can be achieved with decent computational facility. Alternatively, coarse-grain (CG) models are currently being used to capture physics at bigger time and length scale without losing essential underlined structural information. In this study, a CG model of α-hydroxy phytosphingosines (CER[AP]) is developed based on philosophy of MARTINI force field. At first, the model is validated with various atomistic simulations and available experimental data. Later on, the model's compatibility with other major skin lipids, cholesterol, and free fatty acid (palmitic acid) is checked by simulating a mixture of lipid multilayer in presence and absence of water. The developed model of CER[AP] is able to predict key structural properties within the acceptable error limits. The phenomena of ceramide conformation transformation, cholesterol flip-flop, and specificity of lipid arrangement within the multilayered systems is observed during the simulation. This signifies the importance of model in capturing higher order structural transformations.

Predicting asymmetric phospholipid microstructures in solutions

Asymmetric phospholipid microstructures, such as asymmetric phospholipid membranes, have potential applications in biological and medicinal processes. Here, we used the dissipative particle dynamics simulation method to predict the asymmetric phospholipid microstructures in aqueous solutions. The asymmetric phospholipid membranes, tubes and vesicles are determined and characterized by the chain density distributions and order parameters. The phase diagrams are constructed to evaluate the effects of the chain length on the asymmetric structure formations at equilibrium states, while the average radius of gyration and shape factors are calculated to analyze the asymmetric structure formations in the non-equilibrium processes. Meanwhile, we predicted the mechanical properties of the asymmetric membranes by analyzing the spatial distributions of the interface tensions and osmotic pressures in solutions.

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.

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.

Computer-Aided Design of Nanoparticles for Transdermal Drug Delivery

Human skin provides an excellent opportunity for drug delivery application. However, the delivery of hydrophilic drug and big protein molecules is challenging due to barrier provided by the top layer of skin known as stratum corneum (SC). The chemical permeation enhancers or specialized carriers such as nanoparticles (NPs) are needed which can deliver drug molecules into the deeper layer.Here, we describe the in silico design of nanoparticle carriers using molecular dynamics (MD) simulations for the transdermal drug delivery application. At first, setup of a skin lipid bilayer model is demonstrated. Further, nanoparticles are designed based on the Monte Carlo simulation technique. These nanoparticles are then tested on skin model using various MD simulation techniques.

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.

Molecular Dynamics Simulations of Membrane Permeability

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

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.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Extracellular vesicles derived from Malassezia furfur stimulate IL-6 production in keratinocytes as demonstrated in in vitro and in vivo models

BACKGROUND

Malassezia is one of the commensal microorganisms colonized on human skin and has been shown to be related to several inflammatory cutaneous disorders. Previous studies indicated that Malassezia. sympodialis (M. sympodialis) can produce extracellular vesicles, however, the immunoregulatory function of Malassezia extracellular vesicles on keratinocytes has not been studied.

OBJECTIVE

To investigate the extracellular vesicular production capability of Malassezia. furfur (M. furfur) and examine their immunoregulatory effects both in vitro and in vivo.

METHODS

Extracellular vesicles derived from M. furfur were isolated by sequential ultracentrifugation procedure. Their structure and diameter were determined by negative stain TEM and NTA, respectively. Confocal microscopy was used to visualize the internalization of these nanoparticles into HaCaT cells and mice epidermal keratinocytes. The expressions of inflammatory cytokines were screened using PCR Array assay and validated in vitro by qPCR and ELISA assays. In vivo cytokine production was measured by the IHC method. The role of NF-κB in such process was evaluated in HaCaT cells by western blot assay.

RESULTS

Our results showed that M. furfur produced ovoid-shaped nanoparticles, which could be then internalized into HaCaT cells, as well as mice epidermal keratinocytes. IL-6 expression was significantly enhanced in response to extracellular vesicular stimulation both in vitro and in vivo, in which process the activation of NF-κB was involved.

CONCLUSION

M. furfur has the ability to release extracellular vesicles, which can be internalized into keratinocytes and promote the production of IL-6 with the involvement of NF-κB dependent pathway. Such findings reveal some important new insights into Malassezia pathogenesis and therapy.

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.

Composition Dependence of Water Permeation Across Multicomponent Gel-Phase Bilayers

The permeability of multicomponent phospholipid bilayers in the gel phase is investigated via molecular dynamics simulation. The physical role of the different molecules is probed by comparing multiple mixed-component bilayers containing distearylphosphatidylcholine (DSPC) with varying amounts of either the emollient isostearyl isostearate or long-chain alcohol (dodecanol, octadecanol, or tetracosanol) molecules. Permeability is found to depend on both the tail packing density and hydrogen bonding between lipid headgroups and water. Whereas the addition of emollient or alcohol molecules to a gel-phase DSPC bilayer can increase the tail packing density, it also disturbed the hydrogen-bonding network, which in turn can increase interfacial water dynamics. These phenomena have opposing effects on bilayer permeability, which is found to depend on the balance between enhanced tail packing and decreased hydrogen bonding.

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.

Transdermal cellular membrane penetration of proteins with gold nanoparticles: a molecular dynamics study

Transdermal delivery, where the skin acts as the route for local or systemic distribution, presents a lot of advantages over conventional routes such as oral and intravenous and intramuscular injections. However, the delivery of large biomolecules like proteins through the skin is challenging due to their size and structural properties. A molecular level understanding of their transport across the skin barrier is desirable to design successful formulations. We have employed constrained and unconstrained coarse grained molecular dynamics simulation techniques to obtain the molecular mechanism of penetration of the horseradish peroxidase (HRP) protein into the skin, in the presence and absence of gold nanoparticles (AuNPs). Unconstrained simulations show that HRP, when considered individually, was not able to breach the skin barrier, while in the presence of AuNPs, it first binds to the AuNPs and then breaches the barrier. The constrained simulations revealed that there was a free energy barrier for HRP to permeate inside the skin lipid layer when taken alone, while in the presence of gold nanoparticles, no barrier was found. Our study opens up the field of computational modeling based design of nanoparticle carriers for a given protein's transdermal delivery.

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

Accurate in-silico models are required to predict the release of drug molecules through skin in order to supplement the in-vivo experiments for faster development/testing of drugs. The upper most layer of the skin, stratum corneum (SC), offers the main resistance for permeation of actives. Most of the SC's molecular level models comprise cholesterol and phospholipids only, which is far from reality. In this study we have implemented a multiscale modeling framework to obtain the release profile of three drugs, namely, caffeine, fentanyl, and naphthol, through skin SC. We report for the first time diffusion of drugs through a realistic skin molecular model comprised of ceramides, cholesterol, and free fatty acid. The diffusion coefficients of drugs in the SC lipid matrix were determined from multiple constrained molecular dynamics simulations. The calculated diffusion coefficients were then used in the macroscopic models to predict the release profiles of drugs through the SC. The obtained release profiles were in good agreement with available experimental data. The partition coefficient exhibits a greater effect on the release profiles. The reported multiscale modeling framework would provide insight into the delivery mechanisms of the drugs through the skin and shall act as a guiding tool in performing targeted experiments to come up with a suitable delivery system.

Effect of Size and Surface Charge of Gold Nanoparticles on their Skin Permeability: A Molecular Dynamics Study

Molecular level understanding of permeation of nanoparticles through human skin establishes the basis for development of novel transdermal drug delivery systems and design and formulation of cosmetics. Recent experiments suggest that surface coated nano-sized gold nanoparticles (AuNPs) can penetrate the rat and human skin. However, the mechanisms by which these AuNPs penetrate are not well understood. In this study, we have carried out coarse grained molecular dynamics simulations to explore the permeation of dodecanethiol coated neutral hydrophobic AuNPs of different sizes (2–5 nm) and surface charges (cationic and anionic) through the model skin lipid membrane. The results indicate that the neutral hydrophobic AuNPs disrupted the bilayer and entered in it with in ~200 ns, while charged AuNPs were adsorbed on the bilayer headgroup. The permeation free energy calculation revealed that at the head group of the bilayer, a very small barrier existed for neutral hydrophobic AuNP while a free energy minimum was observed for charged AuNPs. The permeability was maximum for neutral 2 nm gold nanoparticle (AuNP) and minimum for 3 nm cationic AuNP. The obtained results are aligned with recent experimental findings. This study would be helpful in designing customized nanoparticles for cosmetic and transdermal drug delivery application.

Molecular dynamics simulation study of translocation of fullerene C60 through skin bilayer: effect of concentration on barrier properties

The molecular level permeation mechanism of fullerenes and its derivatives through human skin could open a vast area for designing novel nanoparticles for cosmetics and drug delivery applications. In this study, we report the permeation mechanism of pristine fullerene C₆₀ for the first time through the skin lipid layer, as determined via prolonged unconstrained and constrained coarse-grained molecular dynamics simulations. The skin layer was modelled as an equimolar ratio of ceramides, cholesterol and free fatty acids. It was observed that at lower concentrations fullerenes formed small clusters (3 or 5 molecules) in the aqueous phase, which further spontaneously permeated inside the bilayer and remained dispersed inside the bilayer interior. On the other hand, at higher concentrations fullerenes aggregated in the aqueous layer, penetrated in that form and remained aggregated in the bilayer interior. Lower concentrations of fullerenes did not induce significant structural changes in the bilayer, whereas at higher concentrations undulations were observed. The permeability of fullerene molecules was found to be concentration-dependent and was explained in terms of their free energy of permeation (thermodynamics) and diffusivity (dynamics). On the basis of the aggregation and dispersion of fullerenes, an optimum fullerene concentration was determined, which could be used for drug delivery and cosmetic applications.