Search

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

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). In inflammatory skin conditions like atopic dermatitis and psoriasis, there are changes in the SC CER composition, such as an increased concentration of a sphingosine-based CER (CER NS) and a reduced concentration of a phytosphingosine-based CER (CER NP). In the present study, a lipid model was created exclusively forming the SPP, to examine whether alterations in the CER NS:CER NP molar ratio would affect the lipid organization. Experimental data were combined with molecular dynamics simulations of lipid models containing CER NS:CER NP at ratios of 1:2 (mimicking a healthy SC ratio) and 2:1 (observed in inflammatory skin diseases), mixed with CHOL and lignoceric acid as the FFA. The experimental findings show that the acyl chains of CER NS and CER NP and the FFA are in close proximity within the SPP unit cell, indicating that CER NS and CER NP adopt a linear conformation, similarly as observed for the LPP. Both the experiments and simulations indicate that the lamellar organization is the same for the two CER NS:CER NP ratios while the SPP NS:NP 1:2 model had a slightly denser hydrogen bonding network than the SPP NS:NP 2:1 model. The simulations show that this might be attributed to intermolecular hydrogen bonding with the additional hydroxide group on the headgroup of CER NP compared with CER NS.

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. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

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., (Copyright © 2017 Elsevier Inc. All rights reserved.)

The tribology of 2-methacryloyloxyethyl phosphorylcholine monolayers in water is studied using molecular dynamics simulations. Our results show two distinct shear regimes where the first is dominated by hydration lubrication, exhibiting near zero friction coefficients, and the second by chain-chain interactions, resembling monomer-supported lubrication. These results provide insight into the hydration lubrication mechanism - a phenomena thought to underlie the extremely efficient lubrication provided by surfaces functionalized with polyzwitterionic polymer brushes and the mammalian synovial joint.

Two different techniques - replica-exchange Wang-Landau (REWL) and statistical temperature molecular dynamics (STMD) - were applied to systematically study the phase transition behavior of self-assembling lipids as a function of temperature using an off-lattice lipid model. Both methods allow the direct calculation of the density of states with improved efficiency compared to the original Wang-Landau method. A 3-segment model of amphiphilic lipids solvated in water has been studied with varied particle interaction energies (ε) and lipid concentrations. The phase behavior of the lipid molecules with respect to bilayer formation has been characterized through the calculation of the heat capacity as a function of temperature, in addition to various order parameters and general visual inspection. The simulations conducted by both methods can go to very low temperatures with the whole system exhibiting well-ordered structures. With optimized parameters, several bilayer phases are observed within the temperature range studied, including gel phase bilayers with frozen water, mixed water (i.e., frozen and liquid water), and liquid water, and a more fluid bilayer with liquid water. The results obtained from both methods, STMD and REWL, are consistently in excellent agreement with each other, thereby validating both the methods and the results.

Because of small surface area to volume ratios nanoscale devices can exhibit dominant surface forces that can quickly degrade unlubricated contacting surfaces. While fluorinated materials have been widely used as lubricants, because of their low critical surface tension and high thermal and mechanical stability, fluorinated monolayer coatings, which are suitable for lubricating nanoscale devices, are less effective as lubricants. Although fluorinated monolayers are more stable than their hydrocarbon counterparts against elevated temperature and humidity, they are known to exhibit higher frictional forces. To overcome this issue, here we study mixed monolayers composed of both hydrocarbon and fluorocarbon chains. Hydrocarbon-based monolayers have been widely studied and shown to improve frictional properties and device life. To investigate the frictional behavior of mixed fluorocarbon/hydrocarbon monolayers, molecular dynamics simulations of pure hydrogenated and fluorinated chains and mixed fluorinated/hydrogenated chains on silica surfaces have been performed. The adhesion and friction between the nanoconfined monolayers as a function of normal load, chain length, and chemical composition of the monolayer coating have been investigated, and mixed fluorocarbon/hydrocarbon monolayers found to outperform both pure fluorocarbon and pure hydrocarbon monolayers. Surface coverage was found to have a significant effect on the performance of all systems studied with higher surface coverage resulting in lower frictional forces. The simulations also show that when the hydrocarbon chains in the monolayer are longer than the fluorocarbon chains, a liquidlike layer is formed by the longer hydrocarbon chains that protrudes above the shorter fluorocarbon chains and aids in friction reduction. A frictional load dependence is also seen in these mixed monolayer systems because of repulsive interactions between the fluorocarbon base layer and the hydrocarbon liquidlike layer. A chain length difference of eight carbons between the base layer and the liquidlike layer was found to provide the lowest friction, while both a larger (because of increased entanglement) and a smaller (insufficient atoms between the contacting base layers to form a liquidlike layer) chain length difference increased friction.

Proteins employ aromatic residues for carbohydrate binding in a wide range of biological functions. Glycoside hydrolases, which are ubiquitous in nature, typically exhibit tunnels, clefts, or pockets lined with aromatic residues for processing carbohydrates. Mutation of these aromatic residues often results in significant activity differences on insoluble and soluble substrates. However, the thermodynamic basis and molecular level role of these aromatic residues remain unknown. Here, we calculate the relative ligand binding free energy by mutating tryptophans in the Trichoderma reesei family 6 cellulase (Cel6A) to alanine. Removal of aromatic residues near the catalytic site has little impact on the ligand binding free energy, suggesting that aromatic residues immediately upstream of the active site are not directly involved in binding, but play a role in the glucopyranose ring distortion necessary for catalysis. Removal of aromatic residues at the entrance and exit of the Cel6A tunnel, however, dramatically impacts the binding affinity, suggesting that these residues play a role in chain acquisition and product stabilization, respectively. The roles suggested from differences in binding affinity are confirmed by molecular dynamics and normal mode analysis. Surprisingly, our results illustrate that aromatic-carbohydrate interactions vary dramatically depending on the position in the enzyme tunnel. As aromatic-carbohydrate interactions are present in all carbohydrate-active enzymes, these results have implications for understanding protein structure-function relationships in carbohydrate metabolism and recognition, carbon turnover in nature, and protein engineering strategies for biomass utilization. Generally, these results suggest that nature employs aromatic-carbohydrate interactions with a wide range of binding affinities for diverse functions.

The cellular membrane functions as a regulating barrier between the intracellular and extracellular regions. For a molecule to reach the interior of the cell from the extracellular fluid, it must diffuse across the membrane, via either active or passive transport. The rigid structure of lipid bilayers, which are a key component of cellular membranes, prohibit simple diffusion of most particles, while vital nutrients are transported to the interior by specific mechanisms, such as ion channels and transport proteins. Although the cellular membrane provides the cell with protection against unwanted toxins that may be in the extracellular medium, some foreign particles can reach the interior of the cell, resulting in irregularities in cellular function. This behavior is particularly noted for permeants with compact molecular structure, suggesting that common nanoscale building blocks, such as fullerenes, may enter into the interior of a cell. To gauge the propensity for such particles to cross the membrane, we have computed the Gibbs free energy of transfer along the axis normal to the bilayer surface for two nanoscale building blocks, C(60) and a hydrogen-terminated polyhedral oligomeric silsequioxane (H-POSS) monomer, in a hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer using molecular dynamics simulations and potential of mean force calculations. The studies show that C(60) has a substantial energetic preference for the soft polymer region of the lipid bilayer system, below the water/bilayer interface, with a transition energy from bulk water of -19.8 kcal/mol. The transition of C(60) from the bulk water to the center of the bilayer, while also energetically favorable, has to overcome a +5.9 kcal/mol energetic barrier in the hydrophobic lipid tail region. The H-POSS simulations indicate an energy minimum at the water-bilayer interface, with an energy of -10.9 kcal/mol; however, a local minimum of -2.7 kcal/mol is also observed in the hydrophobic dense aliphatic region. The energy barrier seen in the hydrophobic core region of the C(60) study is likely due to the significant penalty associated with inserting the relatively large particle into such a dense region. In contrast, whereas H-POSS is found to be subject to an energetic penalty upon insertion into the bilayer, the relatively small size of the H-POSS solute renders this penalty less significant. The energy barrier seen in the soft polymer region for the H-POSS monomer is primarily attributed to the lack of favorable solute-bilayer electrostatic interactions, which are present in the interfacial region, and fewer van der Waals interactions in the soft polymer region than the dense aliphatic region. These results indicate that C(60) may partition into the organic phase of the DPPC/water system, given the favorable free energies in the soft polymer and dense aliphatic regions of the bilayer, and H-POSS is likely to partition near the water-bilayer interface, where the particle has low-energy electrostatic interactions with the polar head groups of the bilayer.

Abstract: Molecular dynamics (MD) simulations play an important role in materials design. However, the effective use of the most widely used MD simulators require significant expertise of the scientific domain and deep knowledge of the given software tool itself. In this paper, we present a tool that offers an intuitive, component-oriented approach to design complex molecular systems and set up initial conditions of the simulations. We integrate this tool into a web- and cloud-based software infrastructure, called MetaMDS, that lowers the barrier of entry into MD simulations for practitioners. The web interface makes it possible for experts to build a rich library of simulation components and for ordinary users to create full simulations by parameterizing and composing the components. A visual programming interface makes it possible to create optimization workflows where the simulators are invoked multiple times with various parameter configurations based on results of earlier simulation runs. Simulation configurations including the various parameters, the version of tools utilized and the results are stored in a database to support searching and browsing of existing simulation outputs and facilitating the reproducibility of scientific results. [Copyright &y& Elsevier]