Phase Diagram of a Stratum Corneum Lipid Mixture

Molecular dynamics simulations of lipid composition and its impact on structural and dynamic properties of skin membrane

The stratum corneum (SC) plays the most important role in the absorption of topical and transdermal drugs. In this study, we developed a multi-layered SC model using coarse-grained molecular dynamics (CGMD) simulations of ceramides, cholesterol, and fatty acids in equimolar proportions, starting from two different initial configurations. In the first approach, all ceramide molecules were initially in the hairpin conformation, and the membrane bilayers were pre-formed. In the second approach, ceramide molecules were introduced in either the hairpin or splayed conformation, with the lipid molecules randomly oriented at the start of the simulation. The aim was to evaluate the effects of lipid chain length on the structural and dynamic properties of SC. By incorporating ceramides and fatty acids of different chain lengths, we simulated the SC membrane in healthy and diseased states. We calculated key structural properties including the thickness, normalized lipid area, lipid tail order parameters, and spatial ordering of the lipids from each system. The results showed that systems with higher ordering and structural integrity contained an equimolar ratio of ceramides (chain length of 24 carbon atoms), fatty acids with chain lengths ≥ of 20 carbon atoms, and cholesterol. In these systems, strong apolar interactions between the ceramide and fatty acid long acyl chains restricted the mobility of the lipid molecules, thereby maintaining a compact lipid headgroup region and high order in the lipid tail region. The simulations also revealed distinct flip-flop mechanisms for cholesterol and fatty acid within the multi-layered membrane. Cholesterol is mostly diffused through the tail-tail interface region of the membrane and could flip-flop in the same bilayer. In contrast, fatty acids flip-flopped between adjacent leaflets of two bilayers in which the tails crossed the thinner headgroup region of the membrane. To conclude, our SC model provides mechanistic insights into lipid mobility and is flexible in its design and composition of different lipids, enabling studies of varying skin conditions.

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.

Mechanisms of the Drug Penetration Enhancer Propylene Glycol Interacting with Skin Lipid Membranes

Very few drugs have the necessary physicochemical properties to cross the skin's main permeability barrier, the stratum corneum (SC), in sufficient amounts. Propylene glycol (PG) is a chemical penetration enhancer that could be included in topical formulations in order to overcome the barrier properties of the skin and facilitate the transport of drugs across it. Experiments have demonstrated that PG increases the mobility and disorder of SC lipids and may extract cholesterol from the SC, but little is known about the molecular mechanisms of drug permeation enhancement by PG. In this work, we have performed molecular dynamics (MD) simulations to investigate the molecular-level effects of PG on the structure and properties of model SC lipid bilayers. The model bilayers were simulated in the presence of PG concentrations over the range of 0-100% w/w PG, using both an all-atom and a united atom force field. PG was found to localize in the hydrophilic headgroup regions at the bilayer interface, to occupy the lipid-water hydrogen-bonding sites, and to slightly increase lipid tail disorder in a concentration-dependent manner. We showed with MD simulation that PG enhances the permeation of small molecules such as water by interacting with the bilayer interface; the results of our study may be used to guide the design of formulations for transdermal drug delivery with enhanced skin permeation, as well as topical formulations and cosmetic products.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Phytosphingosine ceramide mainly localizes in the central layer of the unique lamellar phase of skin lipid model systems

Understanding the lipid arrangement within the skin's outermost layer, the stratum corneum (SC), is important for advancing knowledge on the skin barrier function. The SC lipid matrix consists of ceramides (CERs), cholesterol, and free fatty acids, which form unique crystalline lamellar phases, referred to as the long periodicity phase (LPP) and short periodicity phases. As the SC lipid composition is complex, lipid model systems that mimic the properties of native SC are used to study the SC lipid organization and molecular arrangement. In previous studies, such lipid models were used to determine the molecular organization in the trilayer structure of the LPP unit cell. The aim of this study was to examine the location of CER N-(tetracosanoyl)-phytosphingosine (CER NP) in the unit cell of this lamellar phase and compare its position with CER N-(tetracosanoyl)-sphingosine (CER NS). We selected CER NP as it is the most prevalent CER subclass in the human SC, and its location in the LPP is not known. Our neutron diffraction results demonstrate that the acyl chain of CER NP was positioned in the central part of the trilayer structure, with a fraction also present in the outer layers, the same location as determined for the acyl chain of CER NS. In addition, our Fourier transformed infrared spectroscopy results are in agreement with this molecular arrangement, suggesting a linear arrangement for the CER NS and CER NP. These findings provide more detailed insight into the lipid organization in the SC lipid matrix.

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 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.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

Arrangement of Ceramides in the Skin: Sphingosine Chains Localize at a Single Position in Stratum Corneum Lipid Matrix Models

Understanding the structure of the stratum corneum (SC) is essential to understand the skin barrier process. The long periodicity phase (LPP) is a unique trilayer lamellar structure located in the SC. Adjustments in the composition of the lipid matrix, as in many skin abnormalities, can have severe effects on the lipid organization and barrier function. Although the location of individual lipid subclasses has been identified, the lipid conformation at these locations remains uncertain. Contrast variation experiments via small-angle neutron diffraction were used to investigate the conformation of ceramide (CER) N-(tetracosanoyl)-sphingosine (NS) within both simplistic and porcine mimicking LPP models. To identify the lipid conformation of the twin chain CER NS, the chains were individually deuterated, and their scattering length profiles were calculated to identify their locations in the LPP unit cell. In the repeating trilayer unit of the LPP, the acyl chain of CER NS was located in the central and outer layers, while the sphingosine chain was located exclusively in the middle of the outer layers. Thus, for the CER NS with the acyl chain in the central layer, this demonstrates an extended conformation. Electron density distribution profiles identified that the lipid structure remains consistent regardless of the lipid's lateral packing phase, this may be partially due to the anchoring of the extended CER NS. The presented results provide a more detailed insight on the internal arrangement of the LPP lipids and how they are expected to be arranged in healthy skin.

Predicting asymmetric phospholipid microstructures in solutions

Asymmetric phospholipid microstructures, such as asymmetric phospholipid membranes, have potential applications in biological and medicinal processes. Here, we used the dissipative particle dynamics simulation method to predict the asymmetric phospholipid microstructures in aqueous solutions. The asymmetric phospholipid membranes, tubes and vesicles are determined and characterized by the chain density distributions and order parameters. The phase diagrams are constructed to evaluate the effects of the chain length on the asymmetric structure formations at equilibrium states, while the average radius of gyration and shape factors are calculated to analyze the asymmetric structure formations in the non-equilibrium processes. Meanwhile, we predicted the mechanical properties of the asymmetric membranes by analyzing the spatial distributions of the interface tensions and osmotic pressures in solutions.

Twisting and tilting of a mechanosensitive molecular probe detects order in membranes

Lateral forces in biological membranes affect a variety of dynamic cellular processes. Recent synthetic efforts have introduced fluorescent "flippers" as environment-sensitive planarizable push-pull probes that can detect lipid packing and membrane tension, and respond to lipid-induced mechanical forces by a shift in their spectroscopic properties. Herein, we investigate the molecular origin of the mechanosensitivity of the best known flipper, Flipper-TR, by an extended set of molecular dynamics (MD) simulations in membranes of increasing complexity and under different physicochemical conditions, revealing unprecedented details of the sensing process. Simulations enabled by accurate refinement of Flipper-TR force field using quantum mechanical calculations allowed us to unambiguously correlate the planarization of the two fluorescent flippers to spectroscopic response. In particular, Flipper-TR conformation exhibits bimodal distribution in disordered membranes and a unimodal distribution in highly ordered membranes. Such dramatic change was associated with a shift in Flipper-TR excitation spectra, as supported both by our simulated and experimentally-measured spectra. Flipper-TR sensitivity to phase-transition is confirmed by a temperature-jump protocol that alters the lipid phase of an ordered membrane, triggering an instantaneous mechanical twisting of the probe. Simulations show that the probe is also sensitive to surface tension, since even in a naturally disordered membrane, the unimodal distribution of coplanar flippers can be achieved if a sufficiently negative surface tension is applied to the membrane. MD simulations in ternary mixtures containing raft-like nanodomains show that the probe can discriminate lipid domains in phase-separated complex bilayers. A histogram-based approach, called DOB-phase classification, is introduced that can differentiate regions of disordered and ordered lipid phases by comparing dihedral distributions of Flipper-TR. Moreover, a new sensing mechanism involving the orientation of Flipper-TR is elucidated, corroborating experimental evidence that the probe tilt angle is strongly dependent on lipid ordering. The obtained atomic-resolution description of Flipper-TR mechanosensitivity is key to the interpretation of experimental data and to the design of novel mechanosensors with improved spectroscopic properties.

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.