Interaction of menthol with mixed-lipid bilayer of stratum corneum: A coarse-grained simulation study

Curcumin's membrane localization and disruptive effects on cellular processes - insights from neuroblastoma, leukemic cells, and Langmuir monolayers

In therapies, curcumin is now commonly formulated in liposomal form, administered through injections or creams. This enhances its concentration at the cellular level compared to its natural form ingestion. Due to its hydrophobic nature, curcumin is situated in the lipid part of the membrane, thereby modifying its properties and influencing processes The aim of the research was to investigate whether the toxicity of specific concentrations of curcumin, assessed through biochemical tests for the SK-N-SH and H-60 cell lines, is related to structural changes in the membranes of these cells, caused by the localization of curcumin in their hydrophobic regions. Biochemical tests were performed using spectrophotometric methods. Langmuir technique were used to evaluate the interaction of the curcumin with the studied lipids. Direct introduction of curcumin into the membranes alters their physicochemical parameters. The extent of these changes depends on the initial properties of the membrane. In the conducted research, it has been demonstrated that curcumin may exhibit toxicity to human cells. The mechanism of this toxicity is related to its localization in cell membranes, leading to their dysfunction. The sensitivity of cells to curcumin presence depends on the saturation level of their membranes; the more rigid the membrane, the lower the concentration of curcumin causes its disruption.

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.

Unraveling the Molecular Mechanisms of Alcohol-Mediated Skin Permeation Enhancement: Insights from Molecular Dynamics Simulations

The application of alcohols as permeation enhancers in pharmaceutical and cosmetic formulations has attracted considerable attention, owing to their skin permeation-enhancing effect. Nonetheless, the elucidation of the fundamental mechanisms underlying the skin permeation-enhancing effect remains elusive. In this study, molecular dynamics (MD) simulations were employed to investigate the effect of 1,2-propanediol (1,2-PDO), 1,2-butanediol (1,2-BDO), and ethanol (EtOH) on the stratum corneum (SC) model membrane. The results showed that the effect of alcohols on the SC model membrane displayed a concentration-dependent nature. The alcohols can interact with SC lipids and exhibit a remarkable ability to selectively extract free fatty acid (FFA) molecules from the SC model membrane and make the SC looser. Meanwhile, 1,2-BDO and EtOH can penetrate into SC lipid bilayers at higher concentrations, leading to the formation of continuous hydrophilic defects in SC. The FFA extraction and the formation of continuous hydrophilic defects induced ceramide (CER) tail chains to become more disordered and fluid and also weakened the hydrogen bonding (H-bonding) network among SC lipids. Both the FFA extraction and the continuous hydrophilic defect formation endowed alcohols with the permeation-enhancing effect. The constrained simulations revealed that the free energy barriers decreased for the permeation of the hydrophilic model molecule (COL) across the SC model membranes containing alcohols, particularly for 1,2-BDO and EtOH. The possible permeation-enhancing mechanisms of alcohols were proposed correspondingly. This work not only provided a deep understanding of the transdermal permeation-enhancing behavior of alcohols at the molecular level but also provided necessary reference information for designing effective transdermal drug delivery systems in applications.

Evaluation of the in vitro permeation parameters of topical ketoprofen and lidocaine hydrochloride from transdermal Pentravan® products through human skin

In the treatment of pain, especially chronic pain, the rule of multimodal therapy applies, based on various painkillers mechanisms of action. The aim of the conducted study was to evaluate the in vitro penetration of ketoprofen (KET) and lidocaine hydrochloride (LH) through the human skin from a vehicle with transdermal properties. The results obtained with the use of the Franz chamber showed statistically significantly higher penetration of KET from the transdermal vehicle as compared to commercial preparations. It was also shown that the addition of LH to the transdermal vehicle did not change the amount of KET permeated. The study also compared the penetration of KET and LH by adding various excipients to the transdermal vehicle. Comparing the cumulative mass of KET that penetrated after the 24-h study, it was observed that the significantly highest permeation was found for the vehicle containing additionally Tinctura capsici, then for that containing camphor and ethanol, and the vehicle containing menthol and ethanol as compared to that containing Pentravan^® alone. A similar tendency was observed in the case of LH, where the addition of Tinctura capsici, menthol and camphor led to a statistically significant higher penetration. Adding certain drugs such as KET and LH to Pentravan^®, and substances such as menthol, camphor or capsaicin, can be an interesting alternative to administered enteral drugs especially in the group of patients with multiple diseases and polypragmasy.

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.

Topical drug delivery: History, percutaneous absorption, and product development

Topical products, widely used to manage skin conditions, have evolved from simple potions to sophisticated delivery systems. Their development has been facilitated by advances in percutaneous absorption and product design based on an increasingly mechanistic understanding of drug-product-skin interactions, associated experiments, and a quality-by-design framework. Topical drug delivery involves drug transport from a product on the skin to a local target site and then clearance by diffusion, metabolism, and the dermal circulation to the rest of the body and deeper tissues. Insights have been provided by Quantitative Structure Permeability Relationships (QSPR), molecular dynamics simulations, and dermal Physiologically Based PharmacoKinetics (PBPK). Currently, generic product equivalents of reference-listed products dominate the topical delivery market. There is an increasing regulatory interest in understanding topical product delivery behavior under 'in use' conditions and predicting in vivo response for population variations in skin barrier function and response using in silico and in vitro findings.

Translocation of Endo-Functionalized Molecular Tubes across Different Lipid Bilayers: Atomistic Molecular Dynamics Simulation Study

Various artificial receptors, such as calixarenes, cyclodextrins, cucurbit[n]urils, and their acyclic compounds, pliiar[n]arenes, deep cavitands, and molecular tweezers, can permeate the lipid membranes and they are used as drug carriers to improve the drug solubility, stability, and bioavailability. Inspired by these, we have employed atomistic molecular dynamics simulation to examine the effects of endo-functionalized molecular tubes or naphthotubes (host-1a and host-1b) on seven different types of model lipid bilayers and the permeation properties of these receptors through these model lipid bilayers. Lipid types include six model lipid bilayers (POPC, POPE, DOPC, POPG, DPPE, POPE/POPG) and one realistic membrane (Yeast). We observe that these receptors are spontaneously translocated toward these model lipid bilayer head regions and do not proceed further into these lipid bilayer tail regions (reside at the interface between lipid head and lipid tail region), except for the DPPE-containing systems. In the DPPE model lipid bilayer-containing systems (1a-dppe and 1b-dppe), receptor molecules are only adsorbed on the bilayer surface and reside at the interface between lipid head and water. This finding is also supported by the biased free-energy profiles of these translocation processes. Passive transport of these receptors may be possible through these model lipid bilayers (due to low barrier height), except for DPPE bilayer-containing systems (that have a very high energy barrier at the center). The results from these simulations provide insight into the biocompatibility of host-1a or host-1b in microscopic detail. Based on this work, more research is needed to fully comprehend the role of these synthesized receptors as a prospective drug carrier.

Molecular Dynamics Simulations of Micelle Properties and Behaviors of Sodium Lauryl Ether Sulfate Penetrating Ceramide and Phospholipid Bilayers

Molecular dynamics (MD) simulations with the umbrella sampling (US) method were used to investigate the properties of micelles formed by sodium lauryl ether sulfate with two ether groups (SLE2S) and behaviors of corresponding surfactants transferring from micelles to ceramide and DMPC bilayer surfaces. Average micelle radii based on the Einstein-Smoluchowski and Stokes-Einstein relations showed excellent agreement with those measured by dynamic light scattering, while those obtained by evaluating the gyration radius or calculating the distance between the micelle sulfur atoms and center of mass overestimate the radii. As an SLE2S micelle was pulled down to the ceramide bilayer surface in a 400 ns constant-force steered MD (cf-SMD) simulation, the micelle was partially deformed on the bilayer surface, and several SLE2S surfactants easily were partitioned from the micelle into the ceramide bilayer. In contrast, a micelle was not deformed on the DMPC bilayer surface, and SLE2S surfactants were not transferred from the micelle to the DMPC bilayer. Potential of mean force (PMF) calculations revealed that the Gibbs free energy required for an SLE2S surfactant monomer to transfer from a micelle to bulk water can be compensated by decreased Gibbs free energy when an SLE2S monomer transfers into the ceramide bilayer from bulk water. In addition, micelle deformation on the ceramide bilayer surface can reduce the Gibbs free energy barrier required for a surfactant to escape the micelle and help the surfactant partition from the micelle into the ceramide bilayer. An SLE2S surfactant partitioning into the ceramide bilayer is attributed to hydrogen bonding and favorable interactions between the hydrophilic surfactant head and ceramide molecules, which are more dominant than the dehydration penalty during bilayer insertion. Such interactions between surfactant and lipid molecule heads are considerably reduced in DMPC bilayers owing to dielectric screening by water molecules deep inside the head/tail boundary between the DMPC bilayer. This computational work demonstrates the distinct behavior of SLE2S surfactant micelles on ceramide and DMPC bilayer surfaces in terms of variation in Gibbs free energy, which offers insight into designing surfactants used in transdermal drug delivery systems and cosmetics.

Menthol Binding to the Human α4β2 Nicotinic Acetylcholine Receptor Facilitated by Its Strong Partitioning in the Membrane

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

Understanding Interactions of Curcumin with Lipid Bilayers: A Coarse-Grained Molecular Dynamics Study

The interactions of curcumin with various lipid bilayers (POPC, DOPC, oxidized POPC, and oxidized DOPC) and model biomembranes (symmetric bacteria and yeast plasma membranes, as well as asymmetric mammalian plasma membrane) are investigated. A nonlinear thinning effect of curcumin with respect to its concentration is demonstrated in PC membranes and in the yeast. Curcumin induces asymmetry to the symmetric yeast membranes but reduces the degree of asymmetry of the mammalian plasma membranes when the molecule is placed facing the outer leaflets. The molecule is found to diffuse through oxidized PC bilayers, POPC bilayers at a curcumin to lipid ratio C/L = 1/5, yeast membranes at C/L = 1/100, and the mammalian plasma membranes at C/L = 1/5 and when the molecule placed facing the outer leaflets. The results of this work demonstrate that the lipid type, the lipid distribution, and curcumin amount play a critical role in defining the interactions of curcumin with the lipids and their transport behavior through the bilayers.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

Permeation-enhancing effects and mechanisms of borneol and menthol on ligustrazine: A multiscale study using in vitro and coarse-grained molecular dynamics simulation methods

Borneol (BO) and menthol (MEN) are two widely used natural permeation enhancers in the transdermal drug delivery system. In previous studies, their permeation enhancement effects and mechanisms of action on the hydrophobic drug osthole (logP = 3.8) and hydrophilic drug 5-fluorouracil (logP = -0.9) have been studied. In this study, ligustrazine (LTZ), whose logP is 1.3, was used as a model drug to provide a comprehensive understanding of the influence of its logP on the permeation-enhancing effects of BO and MEN. Both BO and MEN enhanced the permeation of LTZ through the skin stratum corneum, as determined using the modified Franz diffusion cell experiment. The enhancement mechanisms were illustrated by coarse-grained molecular dynamics simulations as follows: at low concentrations, the enhancing ratio of MEN was higher than that of BO because of the stronger perturbation effects of MEN on the lipid bilayer, making it looser and facilitating LTZ diffusion. However, at high concentrations, in addition to the diffusion mechanism, BO induced the formation of water channels to improve the permeation of LTZ; however, MEN had no significant effects through this mechanism. Their results were different from those found with osthole and 5-fluorouracil and have been discussed in this study.

A Molecular Interpretation on the Different Penetration Enhancement Effect of Borneol and Menthol towards 5-Fluorouracil

Borneol and menthol are terpenes that are widely used as penetration enhancers in transdermal drug delivery. To explore their penetration-enhancement effects on hydrophilic drugs, 5-fluorouracil (5-FU) was selected as a model drug. An approach that combined in vitro permeation studies and coarse-grained molecular dynamics was used to investigate their penetration-enhancement effect on 5-FU. The results showed that although both borneol and menthol imparted penetration-enhancement effects on 5-FU, these differed in terms of their mechanism, which may account for the observed variations in penetration-enhancement effects. The main mechanism of action of menthol involves the disruption of the stratum corneum (SC) bilayer, whereas borneol involves multiple mechanisms, including the disruption of the SC bilayer, increasing the diffusion coefficient of 5-FU, and inducing the formation of transient pores. The findings of the present study improve our understanding of the molecular mechanism that is underlying 5-FU penetration-enhancement by borneol and menthol, which may be utilized in future investigations and applications.

Chirality-Dependent Interaction of d- and l-Menthol with Biomembrane Models

Chirality plays a vital role in biological membranes and has a significant effect depending on the type and arrangement of the isomer. Menthol has two typical chiral forms, d- and l-, which exhibit different behaviours. l-Menthol is known for its physiological effect on sensitivity (i.e. a cooling effect), whereas d-menthol causes skin irritation. Menthol molecules may affect not only the thermoreceptors on biomembranes, but also the membrane itself. Membrane heterogeneity (lipid rafts, phase separation) depends on lipid packing and acyl chain ordering. Our interest is to elaborate the chirality dependence of d- and l-menthol on membrane heterogeneity. We revealed physical differences between the two optical isomers of menthol on membrane heterogeneity by studying model membranes using nuclear magnetic resonance and microscopic observation.

Cutaneous Penetration-Enhancing Effect of Menthol: Calcium Involvement

Menthol is a naturally occurring terpene used as a penetration enhancer in topical and transdermal formulations. Literature shows a growing interest in menthol's interactions with the transient receptor potential melastatin 8. A decrease in extracellular Ca²⁺ due to the activation of the transient receptor potential melastatin 8 receptor produces inhibition of E-cadherin expression that is responsible for cell-cell adhesion. Because calcium is present in the entire epidermis, the purpose of this study is to evaluate whether the aforementioned properties of menthol are also related to its penetration-enhancing effects. We formulated 16 gels: (i) drug-alone (diphenhydramine or lidocaine), (ii) drug with menthol, (iii) drug, menthol, and calcium channel blocker (CCB; verapamil or diltiazem), and (iv) drug and CCB. In vitro studies showed no effect of the CCB on the release of the drugs either with or without menthol. In vivo experiments were performed for each drug/menthol/CCB combination gel by applying 4 formulations on a shaved rabbit's dorsum on the same day. Dermis concentration profiles were assessed with microdialysis. The gels containing menthol showed higher penetration of drugs than those without whereas the addition of the CCB consistently inhibited the penetration-enhancing effects of menthol. In summary, these findings strongly support the involvement of calcium in the penetration-enhancing effect of menthol.

Natural Terpenes as Penetration Enhancers for Transdermal Drug Delivery

The greatest hindrance for transdermal drug delivery (TDD) is the barrier property of skin, especially the stratum corneum (SC). Various methodologies have been investigated and developed to enhance the penetration of drugs through the skin. Among them, the most popular approach is the application of penetration enhancers (PEs), including natural terpenes, a very safe and effective class of PEs. In the present paper, we focused on terpenes as skin PEs for TDD. The mechanism of their action, the factors affecting their penetration enhancement effect, as well as their possible skin toxicity were discussed. Terpenes abundant in nature have great potential in the development of PEs. Compared to synthetic PEs, natural terpenes have been proved to possess higher enhancement activity. Interaction with SC intercellular lipids is the main mechanism of action for terpenes. The key factor affecting the enhancement effect is the lipophilicity of both terpenes and drug molecules. In addition, a lot of terpenes have also been proved to be much less toxic compared to azone, the classic synthetic PE. In summary, terpenes may be preferred over the chemically synthesized compounds as safe and effective PEs to promote the percutaneous absorption of drugs.

A Multiscale Study on the Penetration Enhancement Mechanism of Menthol to Osthole

Menthol is a widely used penetration enhancer in clinical medicine due to its high efficiency and relative safety. However, details of the penetration enhancement mechanism of menthol on the molecular level is rarely involved in the discussion. In this work, the penetration enhancement (PE) mechanism of menthol is explored by a multiscale method containing molecular dynamics simulations, in vitro penetration experiments, and transmission electron microscopy. Osthole is chosen to be the tested drug due to its common use in external preparations and because it often accompanies menthol as a PE in the preparations. The results show that menthol in each testing concentration can impair the lipid packing of stratum corneum (SC) and promote osthole permeating into SC, and the penetration promoting effect has an optimal concentration. At a low concentration, menthol causes the bilayer to relax with a reduction in thickness and increment in the lipid headgroup area. At a high concentration, menthol destroys the bilayer structure of SC and causes lipids to form a reversed micelle structure. The penetration enhancement mechanism of menthol is characterized mainly by the disruption of the highly ordered SC lipid in low concentrations and an improvement in the partitioning of drugs into the SC in high concentrations. The results can provide some assistance for additional studies and applications of menthol as a penetration enhancer.

Effects of Concentrations on the Transdermal Permeation Enhancing Mechanisms of Borneol: A Coarse-Grained Molecular Dynamics Simulation on Mixed-Bilayer Membranes

Borneol is a natural permeation enhancer that is effective in drugs used in traditional clinical practices as well as in modern scientific research. However, its molecular mechanism is not fully understood. In this study, a mixed coarse-grained model of stratum corneum (SC) lipid bilayer comprised of Ceramide-N-sphingosine (CER NS) 24:0, cholesterol (CHOL) and free fatty acids (FFA) 24:0 (2:2:1) was used to examine the permeation enhancing mechanism of borneol on the model drug osthole. We found two different mechanisms that were dependent on concentrations levels of borneol. At low concentrations, the lipid system maintained a bilayer structure. The addition of borneol made the lipid bilayer loosen and improved drug permeation. The "pull" effect of borneol also improved drug permeation. However, for a strongly hydrophobic drug like osthole, the permeation enhancement of borneol was limited. When most borneol molecules permeated into bilayers and were located at the hydrophobic tail region, the spatial competition effect inhibited drug molecules from permeating deeper into the bilayer. At high concentrations, borneol led to the formation of water pores and long-lived reversed micelles. This improved the permeation of osthole and possibly other hydrophobic or hydrophilic drugs through the SC. Our simulation results were supported by Franz diffusion tests and transmission electron microscope (TEM) experiments.

The physics of stratum corneum lipid membranes

The stratum corneum (SC), the outermost layer of skin, comprises rigid corneocytes (keratin-filled dead cells) in a specialized lipid matrix. The continuous lipid matrix provides the main barrier against uncontrolled water loss and invasion of external pathogens. Unlike all other biological lipid membranes (such as intracellular organelles and plasma membranes), molecules in the SC lipid matrix show small hydrophilic groups and large variability in the length of the alkyl tails and in the numbers and positions of groups that are capable of forming hydrogen bonds. Molecular simulations provide a route for systematically probing the effects of each of these differences separately. In this article, we present the results from atomistic molecular dynamics of selected lipid bilayers and multi-layers to probe the effect of these polydispersities. We address the nature of the tail packing in the gel-like phase, the hydrogen bond network among head groups, the bending moduli expected for leaflets comprising SC lipids and the conformation of very long ceramide lipids in multi-bilayer lipid assemblies.This article is part of the themed issue 'Soft interfacial materials: from fundamentals to formulation'.

A Coarse Grained Model for a Lipid Membrane with Physiological Composition and Leaflet Asymmetry

The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.