The phase behavior of skin-barrier lipids: A combined approach of experiments and simulations.

Skin barrier function is localized in its outermost layer, the stratum corneum (SC), which is comprised of corneocyte cells embedded in an extracellular lipid matrix containing ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs). The unique structure and composition of this lipid matrix are important for skin barrier function. In this study, experiments and molecular dynamics simulation were combined to investigate the structural properties and phase behavior of mixtures containing nonhydroxy sphingosine CER (CER NS), CHOL, and FFA. X-ray scattering for mixtures with varying CHOL levels revealed the presence of the 5.4 nm short periodicity phase in the presence of CHOL. Bilayers in coarse-grained multilayer simulations of the same compositions contained domains with thicknesses of approximately 5.3 and 5.8 nm that are associated with elevated levels, respectively, of CER sphingosine chains with CHOL, and CER acyl chains with FFA chains. The prevalence of the thicker domain increased with decreasing CHOL content. This might correspond to a phase with ∼5.8 nm spacing observed by x-rays (other details unknown) in mixtures with lower CHOL content. Scissoring and stretching frequencies from Fourier transform infrared spectroscopy (FTIR) also indicate interaction between FFA and CER acyl chains and little interaction between CER acyl and CER sphingosine chains, which requires CER molecules to adopt a predominantly extended conformation. In the simulated systems, neighbor preferences of extended CER chains align more closely with the FTIR observations than those of CERs with hairpin ceramide chains. Both FTIR and atomistic simulations of reverse mapped multilayer membranes detect a hexagonal to fluid phase transition between 65 and 80°C. These results demonstrate the utility of a collaborative experimental and simulation effort in gaining a more comprehensive understanding of SC lipid membranes.

Using molecular simulation to understand the skin barrier.

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

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

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

Towards Molecular Simulations that are Transparent, Reproducible, Usable By Others, and Extensible (TRUE).

Systems composed of soft matter (e.g., liquids, polymers, foams, gels, colloids, and most biological materials) are ubiquitous in science and engineering, but molecular simulations of such systems pose particular computational challenges, requiring time and/or ensemble-averaged data to be collected over long simulation trajectories for property evaluation. Performing a molecular simulation of a soft matter system involves multiple steps, which have traditionally been performed by researchers in a "bespoke" fashion, resulting in many published soft matter simulations not being reproducible based on the information provided in the publications. To address the issue of reproducibility and to provide tools for computational screening, we have been developing the open-source Molecular Simulation and Design Framework (MoSDeF) software suite. In this paper, we propose a set of principles to create Transparent, Reproducible, Usable by others, and Extensible (TRUE) molecular simulations. MoSDeF facilitates the publication and dissemination of TRUE simulations by automating many of the critical steps in molecular simulation, thus enhancing their reproducibility. We provide several examples of TRUE molecular simulations: All of the steps involved in creating, running and extracting properties from the simulations are distributed on open-source platforms (within MoSDeF and on GitHub), thus meeting the definition of TRUE simulations.