Mechanisms of lipid extraction from skin lipid bilayers by sebum triglycerides.

The skin surface, our first barrier against the external environment, is covered by the sebum oil, a lipid film composed of sebaceous and epidermal lipids, which is important in the regulation of the hydration level of our skin. Here, we investigate the pathways leading to the transfer of epidermal lipids from the skin lipid bilayer to the sebum. We show that the sebum triglycerides, a major component of sebum, interact strongly with the epidermal lipids and extract them from the bilayer. Using microsecond time scale molecular dynamics simulations, we identify and quantify the free energy associated with the skin lipid extraction process.

Understanding the interactions between sebum triglycerides and water: a molecular dynamics simulation study.

In recent years, sebum oil has been found to play a key role in the regulation of the hydration of the outermost layer of the skin, the stratum corneum. Understanding how a major component of the sebum oil, the triglyceride tri-cis-6-hexadecenoin (TG), interacts with water is an important step in gaining insight into the water regulation function of the sebum oil. Here we use molecular dynamics simulations to investigate the structural and interfacial properties of TG in bulk and at the air and water interface. Our model performs very well in reproducing experimental results, such as density, surface tensions and surface pressure area isotherms. We show that triglyceride molecules in the liquid phase assemble together, through the glycerol group, forming a single percolating network. TG-air interfaces orient the lipids with the interface enriched with the hydrophobic tails and the glycerol groups buried inside. When in contact with water, the TG molecules at the interface orient the glycerol group towards the water phase and adopt a characteristic trident conformation. Water is shown to penetrate the TG layer thanks to the interaction with the oxygen atoms of the TG molecules, which acts as a pathway for water diffusion. The activation energy for the passage of water is found to be ≈9.5kBT at 310 K, showing that the layer is permeable to water diffusion.

Molecular modeling and structural characterization of a high glycine-tyrosine hair keratin associated protein.

High glycine-tyrosine (HGT) proteins are an important constituent of the keratin associated proteins (KAPs) present in human hair. The glassy state physics of hair fibres are thought to be largely regulated by KAPs, which exist in an amorphous state and are readily affected by environmental conditions. However, there are no studies characterizing the individual KAPs. In this paper, we present the first step to fill this gap by computational modeling and experimental studies on a HGT protein, KAP8.1. In particular, we have modeled the three-dimensional structure of this 63-residue protein using homology information from an anti-freeze protein in snow flea. The model for KAP8.1 is characterized by four strands of poly-proline II (or PPII) type helical secondary structures, held together by two cysteine disulphide bridges. Computer simulations confirm the stability of the modelled structure and show that the protein largely samples the PPII and ß-sheet conformations during the molecular dynamics simulations. Spectroscopic studies including Raman, IR and vibrational circular dichroism have also been performed on synthesized KAP8.1. The experimental studies suggest that KAP8.1 is characterised by ß-sheet and PPII structures, largely consistent with the simulation studies. The model built in this work is a good starting point for further simulations to study in greater depth the glassy state physics of hair, including its water sorption isotherms, glass transition, and the effect of HGT proteins on KAP matrix plasticization. These results are a significant step towards our goal of understanding how the properties of hair can be affected and manipulated under different environmental conditions of temperature, humidity, ageing and small molecule additives.

Dissipative particle dynamics: Systematic parametrization using water-octanol partition coefficients.

We present a systematic, top-down, thermodynamic parametrization scheme for dissipative particle dynamics (DPD) using water-octanol partition coefficients, supplemented by water-octanol phase equilibria and pure liquid phase density data. We demonstrate the feasibility of computing the required partition coefficients in DPD using brute-force simulation, within an adaptive semi-automatic staged optimization scheme. We test the methodology by fitting to experimental partition coefficient data for twenty one small molecules in five classes comprising alcohols and poly-alcohols, amines, ethers and simple aromatics, and alkanes (i.e., hexane). Finally, we illustrate the transferability of a subset of the determined parameters by calculating the critical micelle concentrations and mean aggregation numbers of selected alkyl ethoxylate surfactants, in good agreement with reported experimental values.

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.

Fast cholesterol flip-flop and lack of swelling in skin lipid multilayers.

Atomistic simulations were performed on hydrated model lipid multilayers that are representative of the lipid matrix in the outer skin (stratum corneum). We find that cholesterol transfers easily between adjacent leaflets belonging to the same bilayer via fast orientational diffusion (tumbling) in the inter-leaflet disordered region, while at the same time there is a large free energy cost against swelling. This fast flip-flop may play an important role in accommodating the variety of curvatures that would be required in the three dimensional arrangement of the lipid multilayers in skin, and for enabling mechanical or hydration induced strains without large curvature elastic costs.

Free energy predictions of ligand binding to an α-helix using steered molecular dynamics and umbrella sampling simulations.

Free energy prediction of ligand binding to macromolecules using explicit solvent molecular dynamics (MD) simulations is computationally very expensive. Recently, we reported a linear correlation between the binding free energy obtained via umbrella sampling (US) versus the rupture force from steered molecular dynamics (SMD) simulations for epigallocatechin-3-gallate (EGCG) binding to α-helical-rich keratin. This linear correlation suggests a potential route for fast free energy predictions using SMD alone. In this work, the generality of the linear correlation is further tested for several ligands interacting with the α-helical motif of keratin. These molecules have significantly varying properties, i.e., octanol/water partition coefficient (log P), and/or overall charges (oleic acid, catechin, Fe(2+), citric acid, hydrogen citrate, dihydrogen citrate, and citrate). Using the constant loading rate of our previous study of the keratin-EGCG system, we observe that the linear correlation for keratin-EGCG can be extended to other uncharged molecules where interactions are governed by hydrogen bonds and/or a combination of hydrogen bonds and hydrophobic forces. For molecules where interactions with the keratin helix are governed primarily by electrostatics between charged molecules, a second, alternative linear correlation model is derived. While further investigations are needed to expand the molecular space and build a fully predictive model, the current approach represents a promising methodology for fast free energy predictions based on short SMD simulations (requiring picoseconds to nanoseconds of sampling) for defined biomolecular systems.

Lamellar and inverse micellar structures of skin lipids: effect of templating.

The outermost layer of skin comprises rigid nonviable cells (corneocytes) in a layered lipid matrix. Using atomistic simulations we find that the equilibrium phase of the skin lipids is inverse micellar. A model of the corneocyte is used to demonstrate that lamellar layering is induced by the patterned corneocyte wall. The inverse micellar phase is consistent with in vivo observations in regions where corneocyte walls are well separated (lacunar spaces) and in the inner layers of skin, and suggests a functional role in the lipid synthesis pathway in vivo.

Quantitative Prediction of the Structure and Viscosity of Aqueous Micellar Solutions of Ionic Surfactants: A Combined Approach Based on Coarse-Grained MARTINI Simulations Followed by Reverse-Mapped All-Atom Molecular Dynamics Simulations.

We address the problem of the quantitative prediction of micelle formation in dilute aqueous solutions of ionic surfactants using sodium dodecyl sulfate (SDS) as a model system through a computational approach that involves three steps: (a) execution of coarse-grained simulations based on the MARTINI force field (with slightly modified parameters to afford the formation of large micelles); (b) reverse mapping of the final self-assembled coarse-grained configuration into an all-atom configuration; and (c) final relaxation of this all-atom configuration through short-time (on the order of a few tens of nanoseconds), detailed isothermal-isobaric molecular dynamics simulations using the CHARMM36 force field. For a given concentration of the solution in SDS molecules, the modified MARTINI-based coarse-grained simulations lead to the formation of large micelles characterized by mean aggregation numbers above the experimentally observed ones. However, by reintroducing the detailed chemical structure through a strategy that solves a well-defined geometric problem and re-equilibrating, these large micellar aggregates quickly dissolve to smaller ones and equilibrate to sizes that perfectly match the average micelle size measured experimentally at the given surfactant concentration. From the all-atom molecular dynamics simulations, we also deduce the surfactant diffusivity DSDS and the zero-shear rate viscosity, η0, of the solution, which are observed to compare very favorably with the few experimental values that we were able to find in the literature.

Simulation studies of stratum corneum lipid mixtures.

We present atomistic molecular dynamics results for fully hydrated bilayers composed of ceramide NS-24:0, free fatty acid 24:0 and cholesterol, to address the effect of the different components in the stratum corneum (the outermost layer of skin) lipid matrix on its structural properties. Bilayers containing ceramide molecules show higher in-plane density and hence lower rate of passive transport compared to phospholipid bilayers. At physiological temperatures, for all composition ratios explored, the lipids are in a gel phase with ordered lipid tails. However, the large asymmetry in the lengths of the two tails of the ceramide molecule leads to a fluidlike environment at the bilayer midplane. The lateral pressure profiles show large local variations across the bilayer for pure ceramide or any of the two-component mixtures. Close to the skin composition ratio, the lateral pressure fluctuations are greatly suppressed, the ceramide tails from the two leaflets interdigitate significantly, the depression in local density at the interleaflet region is lowered, and the bilayers have lowered elastic moduli. This indicates that the observed composition ratio in the stratum corneum lipid layer is responsible for both the good barrier properties and the stability of the lipid structure against mechanical stresses.

Simulations of skin barrier function: free energies of hydrophobic and hydrophilic transmembrane pores in ceramide bilayers.

Transmembrane pore formation is central to many biological processes such as ion transport, cell fusion, and viral infection. Furthermore, pore formation in the ceramide bilayers of the stratum corneum may be an important mechanism by which penetration enhancers such as dimethylsulfoxide (DMSO) weaken the barrier function of the skin. We have used the potential of mean constraint force (PMCF) method to calculate the free energy of pore formation in ceramide bilayers in both the innate gel phase and in the DMSO-induced fluidized state. Our simulations show that the fluid phase bilayers form archetypal water-filled hydrophilic pores similar to those observed in phospholipid bilayers. In contrast, the rigid gel-phase bilayers develop hydrophobic pores. At the relatively small pore diameters studied here, the hydrophobic pores are empty rather than filled with bulk water, suggesting that they do not compromise the barrier function of ceramide membranes. A phenomenological analysis suggests that these vapor pores are stable, below a critical radius, because the penalty of creating water-vapor and tail-vapor interfaces is lower than that of directly exposing the strongly hydrophobic tails to water. The PMCF free energy profile of the vapor pore supports this analysis. The simulations indicate that high DMSO concentrations drastically impair the barrier function of the skin by strongly reducing the free energy required for pore opening.

Interaction of oleic acid with dipalmitoylphosphatidylcholine (DPPC) bilayers simulated by molecular dynamics.

The fatty acid oleic acid (OA) is known to modulate the structure of membranes, which forms the basis for a number of its important applications including its use as a therapeutic supplement to reduce the risk of cardiovascular disease, in molecule delivery systems such as liposomes, and as a skin permeability enhancer. While a number of studies have investigated the effect of OA on lipid membranes, our understanding of its mechanisms of action at the molecular level remains rudimentary. We have carried out molecular dynamics simulations using coarse-grained models to investigate the interactions of OA at a range of concentrations with a dipalmitoylphosphatidylcholine (DPPC) bilayer in the liquid-crystalline phase. We have also investigated the relative permeability of the bilayers to model hydrophilic and hydrophobic penetrants by means of chemical potential calculations. The results indicate that OA is able to disperse homogeneously into the bilayer at all concentrations without much perturbation. OA appears to slightly weaken the lateral forces between lipid headgroups, and as the concentration of OA increases this manifests itself as a slight decrease in the area compressibility modulus and a minor increase in the diffusion rate of the OA molecules. While the chemical potential profiles showed little or no variation as a function of OA concentration, the frequency of water permeation events was found to double, indicating some OA-induced permeability enhancement. The study suggests that physiological effects of OA are probably more subtle rather than via gross perturbation of the structure, or that its significant effects are restricted to more condensed membrane structures such as the gel phase.

Toward a Standard Protocol for Micelle Simulation.

In this paper, we present protocols for simulating micelles using dissipative particle dynamics (and in principle molecular dynamics) that we expect to be appropriate for computing micelle properties for a wide range of surfactant molecules. The protocols address challenges in equilibrating and sampling, specifically when kinetics can be very different with changes in surfactant concentration, and with minor changes in molecular size and structure, even using the same force field parameters. We demonstrate that detection of equilibrium can be automated and is robust, for the molecules in this study and others we have considered. In order to quantify the degree of sampling obtained during simulations, metrics to assess the degree of molecular exchange among micellar material are presented, and the use of correlation times are prescribed to assess sampling and for statistical uncertainty estimates on the relevant simulation observables. We show that the computational challenges facing the measurement of the critical micelle concentration (CMC) are somewhat different for high and low CMC materials. While a specific choice is not recommended here, we demonstrate that various methods give values that are consistent in terms of trends, even if not numerically equivalent.

Complete Structure of an Epithelial Keratin Dimer: Implications for Intermediate Filament Assembly.

Keratins are cytoskeletal proteins that hierarchically arrange into filaments, starting with the dimer sub-unit. They are integral to the structural support of cells, in skin, hair and nails. In skin, keratin is thought to play a critical role in conferring the barrier properties and elasticity of skin. In general, the keratin dimer is broadly described by a tri-domain structure: a head, a central rod and a tail. As yet, no atomistic-scale picture of the entire dimer structure exists; this information is pivotal for establishing molecular-level connections between structure and function in intermediate filament proteins. The roles of the head and tail domains in facilitating keratin filament assembly and function remain as open questions. To address these, we report results of molecular dynamics simulations of the entire epithelial human K1/K10 keratin dimer. Our findings comprise: (1) the first three-dimensional structural models of the complete dimer unit, comprising of the head, rod and tail domains; (2) new insights into the chirality of the rod-domain twist gained from analysis of the full domain structure; (3) evidence for tri-subdomain partitioning in the head and tail domains; and, (4) identification of the residue characteristics that mediate non-covalent contact between the chains in the dimer. Our findings are immediately applicable to other epithelial keratins, such as K8/K18 and K5/K14, and to intermediate filament proteins in general.

A study on Fe(2+) - α-helical-rich keratin complex formation using isothermal titration calorimetry and molecular dynamics simulation.

Iron binding to protein is common in biological processes of dioxygen transport, electron transfer as well as in stabilizing drug-protein complexes. α-Helix is the most prevalent secondary structure of proteins. In this study, Fe(2+) binding to α-helix has been studied by isothermal titration calorimetry (ITC) and explicitly solvated molecular dynamics (MD) simulation. Ferrous gluconate and α-helix-rich keratin are used for the ITC study and the results revealed followed one set of identical sites binding model. The MD simulations further revealed that only the acidic side-chain functional groups and η(2) (O,O) coordination modes are involved in the binding of Fe(2+) to α-helix. The ITC results also showed that the binding of ferrous gluconate to keratin was entropy driven and the higher the temperature, the stronger the binding free energy. The favorable entropy of Fe(2+) binding to keratin was attributed to the displacement of water molecules on the α-helix surface, and was confirmed via MD simulations. The most stable coordination states of Fe(2+) and α-helix were identified via simulation: Fe(2+) stacks between two glutamic acid side chain carboxylate groups, displacing water molecules. The binding free energies calculated using MD simulation and the theoretical values were in excellent agreement with the ITC results.

Self-consistent field theory for the interactions between keratin intermediate filaments.

BACKGROUND: Keratins are important structural proteins found in skin, hair and nails. Keratin Intermediate Filaments are major components of corneocytes, nonviable horny cells of the Stratum Corneum, the outermost layer of skin. It is considered that interactions between unstructured domains of Keratin Intermediate Filaments are the key factor in maintaining the elasticity of the skin. RESULTS: We have developed a model for the interactions between keratin intermediate filaments based on self-consistent field theory. The intermediate filaments are represented by charged surfaces, and the disordered terminal domains of the keratins are represented by charged heteropolymers grafted to these surfaces. We estimate the system is close to a charge compensation point where the heteropolymer grafting density is matched to the surface charge density. Using a protein model with amino acid resolution for the terminal domains, we find that the terminal chains can mediate a weak attraction between the keratin surfaces. The origin of the attraction is a combination of bridging and electrostatics. The attraction disappears when the system moves away from the charge compensation point, or when excess small ions and/or NMF-representing free amino acids are added. CONCLUSIONS: These results are in concordance with experimental observations, and support the idea that the interaction between keratin filaments, and ultimately in part the elastic properties of the keratin-containing tissue, is controlled by a combination of the physico-chemical properties of the disordered terminal domains and the composition of the medium in the inter-filament region.

Dual-resolution molecular dynamics simulation of antimicrobials in biomembranes.

Triclocarban and triclosan, two potent antibacterial molecules present in many consumer products, have been subject to growing debate on a number of issues, particularly in relation to their possible role in causing microbial resistance. In this computational study, we present molecular-level insights into the interaction between these antimicrobial agents and hydrated phospholipid bilayers (taken as a simple model for the cell membrane). Simulations are conducted by a novel 'dual-resolution' molecular dynamics approach which combines accuracy with efficiency: the antimicrobials, modelled atomistically, are mixed with simplified (coarse-grain) models of lipids and water. A first set of calculations is run to study the antimicrobials' transfer free energies and orientations as a function of depth inside the membrane. Both molecules are predicted to preferentially accumulate in the lipid headgroup-glycerol region; this finding, which reproduces corresponding experimental data, is also discussed in terms of a general relation between solute partitioning and the intramembrane distribution of pressure. A second set of runs involves membranes incorporated with different molar concentrations of antimicrobial molecules (up to one antimicrobial per two lipids). We study the effects induced on fundamental membrane properties, such as the electron density, lateral pressure and electrical potential profiles. In particular, the analysis of the spontaneous curvature indicates that increasing antimicrobial concentrations promote a 'destabilizing' tendency towards non-bilayer phases, as observed experimentally. The antimicrobials' influence on the self-assembly process is also investigated. The significance of our results in the context of current theories of antimicrobial action is discussed.

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

The stratum corneum is the uppermost layer of the skin and acts as a barrier to keep out contaminants and retain moisture. Understanding the molecular structure and behavior of this layer will provide guidance for optimizing its biological function. In this study we use a model mixture comprised of equimolar portions of ceramide NS (24:0), lignoceric acid, and cholesterol to model the effect of the addition of small amounts of oleic acid to the bilayer at 300 and 340 K. Five systems at each temperature have been simulated with concentrations between 0 and 0.1 mol % oleic acid. Our major finding is that subdiffusive behavior over the 200 ns time scale is evident in systems at 340 K, with cholesterol diffusion being enhanced with increased oleic acid. Importantly, cholesterol and other species diffuse faster when radial densities indicate nearest neighbors include more cholesterol. We also find that, with the addition of oleic acid, the bilayer midplane and interfacial densities are reduced and there is a 3% decrease in total thickness occurring mostly near the hydrophilic interface at 300 K with reduced overall density at 340 K. Increased interdigitation occurs independent of oleic acid with a temperature increase. Slight ordering of the long non-hydroxy fatty acid of the ceramide occurs near the hydrophilic interface as a function of the oleic acid concentration, but no significant impact on hydrogen bonding is seen in the chosen oleic acid concentrations.

The permeability enhancing mechanism of DMSO in ceramide bilayers simulated by molecular dynamics.

The lipids of the topmost layer of the skin, the stratum corneum, represent the primary barrier to molecules penetrating the skin. One approach to overcoming this barrier for the purpose of delivery of active molecules into or via the skin is to employ chemical permeability enhancers, such as dimethylsulfoxide (DMSO). How these molecules exert their effect at the molecular level is not understood. We have investigated the interaction of DMSO with gel-phase bilayers of ceramide 2, the predominant lipid in the stratum corneum, by means of molecular dynamics simulations. The simulations satisfactorily reproduce the phase behavior and the known structural parameters of ceramide 2 bilayers in water. The effect of DMSO on the gel-phase bilayers was investigated at various concentrations over the range 0.0-0.6 mol fraction DMSO. The DMSO molecules accumulate in the headgroup region and weaken the lateral forces between the ceramides. At high concentrations of DMSO (> or =0.4 mol fraction), the ceramide bilayers undergo a phase transition from the gel phase to the liquid crystalline phase. The liquid-crystalline phase of ceramides is expected to be markedly more permeable to solutes than the gel phase. The results are consistent with the experimental evidence that high concentrations of DMSO fluidize the stratum corneum lipids and enhance permeability.

Molecular basis for dimethylsulfoxide (DMSO) action on lipid membranes.

Dimethylsulfoxide (DMSO) is an aprotic solvent that has the ability to induce cell fusion and cell differentiation and enhance the permeability of lipid membranes. It is also an effective cryoprotectant. Insights into how this molecule modulates membrane structure and function would be invaluable toward regulating the above processes and for developing chemical means for enhancing or hindering the absorption of biologically active molecules, in particular into or via the skin. We show here by means of molecular simulations that DMSO can induce water pores in dipalmitoyl-phosphatidylcholine bilayers and propose this to be a possible pathway for the enhancement of penetration of actives through lipid membranes. DMSO also causes the membrane to become floppier, which would enhance permeability, facilitate membrane fusion, and enable the cell membrane to accommodate osmotic and mechanical stresses during cryopreservation.