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.

The Sphingosine and Phytosphingosine Ceramide Ratio in Lipid Models Forming the Short Periodicity Phase: An Experimental and Molecular Simulation Study.

The lipids located in the outermost layer of the skin, the stratum corneum (SC), play a crucial role in maintaining the skin barrier function. The primary components of the SC lipid matrix are ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs). They form two crystalline lamellar phases: the long periodicity phase (LPP) and the short periodicity phase (SPP).

E() Equivariant Graph Neural Network for Learning Interactional Properties of Molecules.

We have developed a multi-input E() equivariant graph convolution-based model designed for the prediction of chemical properties that result from the interaction of heterogeneous molecular structures. By incorporating spatial features and constraining the functions learned from these features to be equivariant to E() symmetries, the interactional-equivariant graph neural network (IEGNN) can efficiently learn from the 3D structure of multiple molecules. To verify the IEGNN's capability to learn interactional properties, we tested the model's performance on three molecular data sets, two of which are curated in this study and made publicly available for future interactional model benchmarking.

Global Ring Study to Investigate the Comparability of Total Assay Performance of Commercial Claudin 18 Antibodies for Evaluation in Gastric Cancer.

Claudin 18.2 (CLDN18.2), the dominant isoform of CLDN18 in gastric tissues, is a highly specific tight junction protein of the gastric mucosa with variably retained expressions in gastric and gastroesophageal junction cancers.

The skin barrier: An extraordinary interface with an exceptional lipid organization.

The barrier function of the skin is primarily located in the stratum corneum (SC), the outermost layer of the skin. The SC is composed of dead cells with highly organized lipid lamellae in the intercellular space. As the lipid matrix forms the only continuous pathway, the lipids play an important role in the permeation of compounds through the SC.

MoSDeF-GOMC: Python Software for the Creation of Scientific Workflows for the Monte Carlo Simulation Engine GOMC.

MoSDeF-GOMC is a python interface for the Monte Carlo software GOMC to the Molecular Simulation Design Framework (MoSDeF) ecosystem. MoSDeF-GOMC automates the process of generating initial coordinates, assigning force field parameters, and writing coordinate (PDB), connectivity (PSF), force field parameter, and simulation control files. The software lowers entry barriers for novice users while allowing advanced users to create complex workflows that encapsulate simulation setup, execution, and data analysis in a single script.

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.

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.

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning.

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties.

Examining the self-assembly of patchy alkane-grafted silica nanoparticles using molecular simulation.

In this work, molecular dynamics simulations are used to examine the self-assembly of anisotropically coated "patchy" nanoparticles. Specifically, we use a coarse-grained model to examine silica nanoparticles coated with alkane chains, where the poles of the grafted nanoparticle are bare, resulting in strongly attractive patches. Through a systematic screening process, the patchy nanoparticles are found to form dispersed, string-like, and aggregated phases, dependent on the combination of alkane chain length, coating chain density, and the fractional coated surface area.

Examining Tail and Headgroup Effects on Binary and Ternary Gel-Phase Lipid Bilayer Structure.

The structural properties of two- and three-component gel-phase bilayers were studied using molecular dynamics simulations. The bilayers contain distearoylphosphatidylcholine (DSPC) phospholipids mixed with alcohols and/or fatty acids of varying tail lengths, with carbon chain lengths of 12, 16, and 24 studied. Changes in both headgroup chemistry and tail length are found to affect the balance between steric repulsion and van der Waals attraction within the bilayers, manifesting in different bilayer structural properties.

MoSDeF, a Python Framework Enabling Large-Scale Computational Screening of Soft Matter: Application to Chemistry-Property Relationships in Lubricating Monolayer Films.

We demonstrate how the recently developed Python-based Molecular Simulation and Design Framework (MoSDeF) can be used to perform molecular dynamics screening of functionalized monolayer films, focusing on tribological effectiveness. MoSDeF is an open-source package that allows for the programmatic construction and parametrization of soft matter systems and enables TRUE (transferable, reproducible, usable by others, and extensible) simulations. The MoSDeF-enabled screening identifies several film chemistries that simultaneously show low coefficients of friction and adhesion.

Influence of Single-Stranded DNA Coatings on the Interaction between Graphene Nanoflakes and Lipid Bilayers.

Using molecular dynamics simulations, it is demonstrated that a partial coating of single-stranded DNA (ssDNA) reduces the penetration depth of a graphene nanoflake (GNF) into a phospholipid bilayer by attenuating the hydrophobic force that drives the penetration. As the GNF penetrates the bilayer, the ssDNA remains adsorbed to the GNF outside of the bilayer where it shields the graphene from the surrounding water. The penetration depth is found to be controlled by the amount of ssDNA coating the GNF, with a sparser coating resulting in a deeper penetration since the ssDNA shields less of the GNF surface.

Investigation of the Impact of Cross-Polymerization on the Structural and Frictional Properties of Alkylsilane Monolayers Using Molecular Simulation.

Cross-linked chemisorbed -alkylsilane (CH(CH)Si(OH)) monolayers on amorphous silica surfaces have been studied and their structural properties and frictional performance were compared to those of equivalent monolayers without cross-linkages. The simulations isolated for the first time the effects of both siloxane cross-linkages and the fraction of chains chemisorbed to the surface, providing insight into a longstanding fundamental question in the literature regarding molecular-level structure. The results demonstrate that both cross-linkages and the fraction of chemisorbed chains affect monolayer structure in small but measurable ways, particularly for monolayers constructed from short chains; however, these changes do not appear to have a significant impact on frictional performance.

A Transferable, Multi-Resolution Coarse-Grained Model for Amorphous Silica Nanoparticles.

Despite the ubiquity of nanoparticles in modern materials research, computational scientists are often forced to choose between simulations featuring detailed models of only a few nanoparticles or simplified models with many nanoparticles. Herein, we present a coarse-grained model for amorphous silica nanoparticles with parameters derived via potential matching to atomistic nanoparticle data, thus enabling large-scale simulations of realistic models of silica nanoparticles. Interaction parameters are optimized to match a range of nanoparticle diameters in order to increase transferability with nanoparticle size.

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.

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.

Investigating Alkylsilane Monolayer Tribology at a Single-Asperity Contact with Molecular Dynamics Simulation.

Chemisorbed monolayer films are known to possess favorable characteristics for nanoscale lubrication of micro- and nanoelectromechanical systems (MEMS/NEMS). Prior studies have shown that the friction observed for monolayer-coated surfaces features a strong dependence on the geometry of contact. Specifically, tip-like geometries have been shown to penetrate into monolayer films, inducing defects in the monolayer chains and leading to plowing mechanisms during shear, which result in higher coefficients of friction (COF) than those observed for planar geometries.

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.

Perfluoropolyethers: Development of an All-Atom Force Field for Molecular Simulations and Validation with New Experimental Vapor Pressures and Liquid Densities.

A force field for perfluoropolyethers (PFPEs) based on the general optimized potentials for liquid simulations all-atom (OPLS-AA) force field has been derived in conjunction with experiments and ab initio quantum mechanical calculations. Vapor pressures and densities of two liquid PFPEs, perfluorodiglyme (CF-O-(CF-CF-O)-CF) and perfluorotriglyme (CF-O-(CF-CF-O)-CF), have been measured experimentally to validate the force field and increase our understanding of the physical properties of PFPEs. Force field parameters build upon those for related molecules (e.

Structural Properties of Phospholipid-based Bilayers with Long-Chain Alcohol Molecules in the Gel Phase.

The structural properties of two-component gel-phase bilayers of distearylphosphatidylcholine (DSPC) and alcohol molecules with different compositions and chain lengths (12-24 carbons long) are studied via molecular dynamics simulations. Several bilayer properties, including area per lipid, tilt angle, chain interdigitation, and headgroup offset, are studied for each system and compared, revealing important structural implications depending upon headgroup size and chain length. While tail tilt is the primary mechanism for single-component bilayers to balance tail attraction and headgroup repulsion, our results demonstrate that the lipid mixtures studied adjust this balance via an offset between the depths of the different molecular species in the bilayer; this behavior is found to depend both on composition and on the length of alcohol molecules relative to the length of DSPC tails.

Investigating the Structure of Multicomponent Gel-Phase Lipid Bilayers.

Single- and multicomponent lipid bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphatidylcholine (DSPC), isostearyl isostearate, and heptadecanoyl heptadecanoate in the gel phase are studied via molecular dynamics simulations. It is shown that the structural properties of multicomponent bilayers can deviate strongly from the structures of their single-component counterparts. Specifically, the lipid mixtures are shown to adopt a compact packing by offsetting the positioning depths at which different lipid species are located in the bilayer.

Influence of Surface Morphology on the Shear-Induced Wear of Alkylsilane Monolayers: Molecular Dynamics Study.

Chemisorbed alkylsilane monolayer coatings have been shown to possess favorable lubrication properties; however, film degradation prevents the widespread use of these materials as lubricants in micro- and nanoelectromechanical systems (MEMS/NEMS). In this work, molecular dynamics (MD) simulations are used to provide insight into the conditions that promote the degradation and wear of these materials. This is achieved through removal of interfacial chain-substrate bonds during shear and the examination of the mobility of the resulting free, unbound chains.

Development of a coarse-grained water forcefield via multistate iterative Boltzmann inversion.

A coarse-grained water model is developed using multistate iterative Boltzmann inversion. Following previous work, the -means algorithm is used to dynamically map multiple water molecules to a single coarse-grained bead, allowing the use of structure-based coarse-graining methods. The model is derived to match the bulk and interfacial properties of liquid water and improves upon previous work that used single state iterative Boltzmann inversion.