Molecular dynamics simulations of stratum corneum lipid mixtures: A multiscale perspective

Molecular Dynamics Investigation into CerENP's Effect on the Lipid Matrix of Stratum Corneum

The extracellular lipid matrix in the stratum corneum (SC) plays a critical role in skin barrier functionality, comprising three primary components: ceramides, cholesterol, and free fatty acids. The diverse ceramides, differentiated by molecular structures such as hydroxylations and varying chain lengths, are essential for the lipid matrix's structural integrity. Recently, a new subclass of ceramide, 1-O-acylceramide NP (CerENP), has been identified; however, its precise role in the lipid matrix of the SC is still elusive. Herein, we investigate the role of CerENP on the structure and permeability of the SC using molecular dynamics simulations. Our findings indicate that CerENP contributes to a compact lipid matrix in the lateral dimension of our SC model with a repeat distance of about 13 nm. Additionally, ethanol permeability assessments show that CerENP effectively reduces molecular penetration through the lipid matrix. This study provides an insight into the role of a new subclass of ceramide in the SC, enhancing our understanding of skin structure and the mechanisms behind barrier dysfunction in skin diseases.

The molecular arrangement of ceramides in the unit cell of the long periodicity phase of stratum corneum models shows a high adaptability to different ceramide head group structures

The stratum corneum (SC) lipid matrix, composed primarily of ceramides (CERs), cholesterol and free fatty acids (FFA), has an important role for the skin barrier function. The presence of the long periodicity phase (LPP), a unique lamellar phase, is characteristic for the SC. Insight into the lipid molecular arrangement within the LPP unit cell is imperative for understanding the relationship between the lipid subclasses and the skin barrier function. In this study, the impact of the CER head group structure on the lipid arrangement and barrier functionality was investigated using lipid models forming the LPP. The results demonstrate that the positions of CER N-(tetracosanoyl)-sphingosine (CER NS) and CER N-(tetracosanoyl)-phytosphingosine (CER NP), two essentials CER subclasses, are not influenced by the addition of another CER subclass (N-(tetracosanoyl)-dihydrosphingosine (CER NdS), N-(2R-hydroxy-tetracosanoyl)-sphingosine (CER AS) or D-(2R-hydroxy-tetracosanoyl)-phytosphingosine (CER AP)). However, differences are observed in the lipid organization and the hydrogen bonding network of the three different models. A similar localization of CER NP and CER NS is also observed in a more complex lipid model, with the CER subclass composition mimicking that of human SC. These studies show the adaptability and insensitivity of the LPP unit cell structure to changes in the lipid head group structures of the CER subclasses.

Multifaceted Activity of Fabimycin: Insights from Molecular Dynamics Studies on Bacterial Membrane Models

Membranes─cells' essential scaffolds─are valid molecular targets for substances with an antimicrobial effect. While certain substances, such as octenidine, have been developed to target membranes for antimicrobial purposes, the recently reported molecule, fabimycin (F2B)─a novel agent targeting drug-resistant Gram-negative bacteria─has not received adequate attention regarding its activity on membranes in the literature. The following study aims to investigate the effects of F2B on different bacterial membrane models, including simple planar bilayers and more complex bilayer systems that mimic the Escherichia coli shell equipped with double inner and outer bilayers. Our results show that F2B exhibited more pronounced interactions with bacterial membrane systems compared to the control PC system. Furthermore, we observed significant changes in local membrane property homeostasis in both the inner and outer membrane models, specifically in the case of lateral diffusion, membrane thickness, and membrane resilience (compressibility, tilt). Finally, our results showed that the effect of F2B differed in a complex system and a single membrane system. Our study provides new insights into the multifaceted activity of F2B, demonstrating its potential to disrupt bacterial membrane homeostasis, indicating that its activity extends the currently known mechanism of FabI enzyme inhibition. This disruption, coupled with the ability of F2B to penetrate the outer membrane layers, sheds new light on the behavior of this antimicrobial molecule. This highlights the importance of the interaction with the membrane, crucial in combating bacterial infections, particularly those caused by drug-resistant strains.

Permeability of Dermatological Solutes through the Short Periodicity Phase of Human Stratum Corneum Lipid Bilayers

Determining the permeability of drug-like solutes through the densely packed and heterogeneous stratum corneum lipid layer presents a significant challenge. In this study, we employed umbrella sampling with a periodic weighing function applied to the center of mass of the lipid bilayers. Precise umbrella sampling was conducted with an interframe distance of 0.5 Å, spanning from the bilayer center to the water phase, and each frame was simulated for at least 20 ns. Autocorrelation functions, potential of mean force (PMF), and diffusivity profiles were analyzed for six solutes (testosterone, benzene, caffeine, ethanol, mannitol, and histidine). The results revealed that autocorrelations were dependent solely on the medium, whether water or lipid phase. Diffusivity and PMF profiles along the reaction coordinate were influenced by the hydrophilicity of the solute rather than its size. For hydrophobic solutes, the PMF curves exhibited a minimum at the bilayer center, while for hydrophilic solutes, the PMFs peaked at the bilayer center and lipid tails (where the lipid tails are not interacting with the cholesterol). Diffusivity curves were low at the bilayer center and water phase, with peaks observed at the headgroup or the boundary between fatty acid and cholesterol (1 nm from the bilayer center). The quantitative findings presented in this work hold significance for pharmacists and drug designers.

Lycopene, but not zeaxanthin, serves as a skeleton for the formation of an orthorhombic organization of intercellular lipids within the lamellae in the stratum corneum: Molecular dynamics simulations of the hydrated ceramide NS bilayer model

Carotenoids play an important role in the protection of biomembranes against oxidative damage. Their function depends on the surroundings and the organization of the lipid membrane they are embedded in. Carotenoids are located parallel or perpendicular to the surface of the lipid bilayer. The influence of carotenoids on the organization of the lipid bilayer in the stratum corneum has not been thoroughly considered. Here, the orientation of the exemplary cutaneous carotenoids lycopene and zeaxanthin in a hydrated ceramide NS24 bilayer model and the influence of carotenoids on the lateral organization of the lipid bilayer model were studied by means of molecular dynamics simulations for 32 °C and 37 °C. The results confirm that lycopene is located parallel and zeaxanthin perpendicular to the surface of the lipid bilayer. The lycopene-loaded lipid bilayer appeared to have a strong orthorhombic organization, while zeaxanthin-loaded and pure lipid bilayers were organized in a disordered hexagonal-like and liquid-like state, respectively. The effect is stronger at 32 °C compared to 37 °C based on p-values. Therefore, it was assumed that carotenoids without hydroxyl polar groups in their structure facilitate the formation of the orthorhombic organization of lipids, which provides the skin barrier function. It was shown that the distance between carotenoid atoms matched the distance between atoms in the lipids, indicating that parallel located carotenoids without hydroxyl groups serve as a skeleton for lipid membranes inside the lamellae. The obtained results provide reasonable prediction of the overall qualitative properties of lipid model systems and show the importance of parallel-oriented carotenoids in the development and maintenance of the skin barrier function.

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.

Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants

A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration c_T equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.

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 C₂₄ 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.

Complex nanoemulsion for vitamin delivery: droplet organization and interaction with skin membranes

Lipid nanoemulsions are promising nanomaterials for drug delivery applications in food, pharmaceutical and cosmetic industries. Despite the noteworthy commercial interest, little is known about their supramolecular organization, especially about how such multicomponent formulations interact with cell membranes. In the present work, coarse-grained molecular dynamics simulations have been employed to study the self-assembly of a 15-component lipid nanoemulsion droplet containing vitamins A and E for skin delivery. Our results display aspects of the unique "onion-like" agglomeration between the chemical constituents in the different layers of the lipid nanodroplet. Vitamin E molecules are more concentrated in the center of the droplet together with other hydrophobic constituents such as the triglycerides with long tails. On the other hand, vitamin A occupies an intermediate layer between the core and the co-emulsifier surface of the nanodroplet, together with lecithin phospholipids. Coarse-grained molecular dynamics simulations were also performed to provide insight into the first steps involved in absorption and penetration of the nanodroplet through skin membrane models, representing an intracellular (hair follicle infundibulum) and intercellular pathway (stratum corneum) through the skin. Our data provide a first view on the complex organization of commercial nanoemulsion and its interaction with skin membranes. We expect our results to open the way towards the rational design of such nanomaterials.

Modeling the Percutaneous Absorption of Solvent-deposited Solids Over a Wide Dose Range

The transient absorption of two skin care agents, niacinamide (nicotinamide, NA) and methyl nicotinate (MN), solvent-deposited on ex vivo human skin mounted in Franz diffusion cells has been analyzed according to a new variation on a recently published mechanistic skin permeability model (Yu et al. 2020. J Pharm Sci 110:2149-56). The model follows the absorption and evaporation of two components, solute and solvent, and it includes both a follicular transport component and a dissolution rate limitation for high melting, hydrophilic solids deposited on the skin. Explicit algorithms for improving the simulation of transient diffusion of solvent-deposited solids are introduced. The simulations can account for the ex vivo skin permeation time course of both NA and MN over a dose range exceeding 4.5 orders of magnitude. The model allows one to describe on a mechanistic basis why the percutaneous absorption rate of NA is approximately 60-fold lower than that of its lower melting, more lipophilic analog, MN. It furthermore suggests that MN perturbs stratum corneum barrier lipids and increases their permeability while NA does not, presenting a challenge to molecular modelers engaged in simulating biological lipid barriers.

Surface science of cosmetic substrates, cleansing actives and formulations

The development of shampoo and cleansing formulations in cosmetics is at a crossroads due to consumer demands for better performing, more natural products and also the strong commitment of cosmetic companies to improve the sustainability of cosmetic products. In order to go beyond traditional formulations, it is of great importance to clearly establish the science behind cleansing technologies and appreciate the specificity of cleansing biological surfaces such as hair and skin. In this review, we present recent advances in our knowledge of the physicochemical properties of the hair surface from both an experimental and a theoretical point of view. We discuss the opportunities and challenges that newer, sustainable formulations bring compared to petroleum-based ingredients. The inevitable evolution towards more bio-based, eco-friendly ingredients and sustainable formulations requires a complete rethink of many well-known physicochemical principles. The pivotal role of digital sciences and modelling in the understanding and conception of new ingredients and formulations is discussed. We describe recent numerical approaches that take into account the specificities of the hair surface in terms of structuration, different methods that study the adsorption of formulation ingredients and finally the success of new data-driven approaches. We conclude with practical examples on current formulation efforts incorporating bio-surfactants, controlling foaming and searching for new rheological properties.

Mechanisms of skin moisturization with hyperharmonized hydroxyl modified fullerene substance

BACKGROUND

Hyper Harmonized Hydroxyl Modified Fullerene Substance (3HFWC⁺ ) establishes hydrogen bonds with the surrounding water molecules and organizes them in clusters with the liquid crystalline state, similar to the properties of water surrounding the biomolecules.

AIMS

To investigate the moisturizing properties of hyperharmonized fullerenol-3HFWC⁺ as an emulsion O/W ingredient on the skin.

PATIENTS/METHODS

We have analyzed the reflexion of the blue light from the skin, in vivo, with different levels of moisturization and compared the influence of three groups of cosmetic products (with various active ingredients, 3HFWC⁺ or water in same percentage in the "vehiculum") on skin moisturization by measuring paramagnetic/diamagnetic properties by Optomagnetic Imaging Spectroscopy.

RESULTS

Regenerating Cream and Body Lotion have shown statistically significant increase of diamagnetic features predominantly in all 3 types of creams, while the Hand Cream formulations have not shown statistically significant results after 4 weeks of treatment. Anti-aging cream with 3HFWC⁺ demonstrated statistically significant increase of paramagnetic properties of the skin (p⁺ /p^- from 0.90 to 0.62) and showed beneficial effects on information transfer and water content in stratum corneum.

CONCLUSION

Although it does not act as a classical moisturizing agent on the skin, 3HFWC⁺ can increase moisturization as a result of hydrogen bonds established among 3HFWC⁺ and water or biomolecules, liquid crystalline state of water in 3HFWC⁺ resembling the properties of water surrounding the biomolecules, and finally, liquid crystalline state of water in the products.

Coexistence of Lipid Phases Stabilizes Interstitial Water in the Outer Layer of Mammalian Skin

The lipid matrix in the outer layer of mammalian skin, the stratum corneum, has been previously investigated by multiple biophysical techniques aimed at identifying hydrophilic and lipophilic pathways of permeation. Although consensus is developing over the microscopic structure of the lipid matrix, no molecular-resolution model describes the permeability of all chemical species simultaneously. Using molecular dynamics simulations of a model mixture of skin lipids, the self-assembly of the lipid matrix lamellae has been studied. At higher humidity, the resulting lamellar phase is maintained by partitioning excess water into isolated droplets of controlled size and spatial distribution. The droplets may fuse together to form intralamellar water channels, thereby providing a pathway for the permeation of hydrophilic species. These results reconcile competing data on the outer skin's structure and broaden the scope of molecular-based methods to improve the safety of topical products and to advance transdermal drug delivery.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.