In Silico Study on the Contribution of the Follicular Route to Dermal Permeability of Small Molecules.

PURPOSE: This study investigates in silico the contribution of the hair follicle to the overall dermal permeability of small molecules, as published experimental work provides inconclusive information on whether the follicular route favours the permeation of hydrophobic or hydrophilic permeants. METHOD: A study is conducted varying physico-chemical parameters of permeants such as lipophilicity, molecular weight and protein binding. The simulated data is compared to published experimental data to discuss how those properties can modulate the contribution of the hair follicle to the overall dermal permeation. RESULTS: The results indicate that the contribution of the follicular route to dermal permeation can range from negligible to notable depending on the combination of lipophilic/hydrophilic properties of the substance filling the follicular route and the permeant. CONCLUSION: Characterisation of the substance filling the follicular route is required for analysing the experimental data of dermal permeation of small molecules, as changes between in vivo and in vitro due to handling of samples and cessation of vital functions can modify the contribution of the follicular route to overall dermal permeation, hence hindering data interpretation.

Computational Modelling of the Impact of Evaporation on In-Vitro Dermal Absorption.

PURPOSE: Volatiles are common in personal care products and dermatological drugs. Determining the impact of evaporation of volatiles on skin permeation is crucial to evaluate and understand their delivery, bioavailability, efficacy and safety. We aim to develop an in-silico model to simulate the impact of evaporation on the dermal absorption of volatiles. METHOD: The evaporation of volatile permeants was modelled using vapour pressure as the main factor. This model considers evaporation as a passive diffusion process driven by the concentration gradient between the air-vehicle interface and the ambient environment. The evaporation model was then integrated with a previously published physiologically based pharmacokinetic (PBPK) model of skin permeation and compared with published in vitro permeation test data from the Cosmetics Europe ADME Task Force. RESULTS: The evaporation-PBPK model shows improved predictions when evaporation is considered. In particular, good agreement has been obtained for the distributions in the evaporative loss, and the overall percutaneous absorption. The model is further compared with published in-silico models from the Cosmetics Europe ADME Task Force where favourable results are achieved. CONCLUSION: The evaporation of volatile permeants under finite dose in vitro permeation test conditions has been successfully predicted using a mechanistic model with the intrinsic volatility parameter vapour pressure. Integrating evaporation in PBPK modelling significantly improved the prediction of dermal delivery.

A novel multi-omics data analysis of dose-dependent and temporal changes in regulatory pathways due to chemical perturbation: a case study on caffeine.

Comprehensive analysis of multi-omics data can reveal alterations in regulatory pathways induced by cellular exposure to chemicals by characterizing biological processes at the molecular level. Data-driven omics analysis, conducted in a dose-dependent or dynamic manner, can facilitate comprehending toxicity mechanisms. This study introduces a novel multi-omics data analysis designed to concurrently examine dose-dependent and temporal patterns of cellular responses to chemical perturbations. This analysis, encompassing preliminary exploration, pattern deconstruction, and network reconstruction of multi-omics data, provides a comprehensive perspective on the dynamic behaviors of cells exposed to varying levels of chemical stimuli. Importantly, this analysis is adaptable to any number of omics layers, including site-specific phosphoproteomics. We implemented this analysis on multi-omics data obtained from HepG2 cells exposed to a range of caffeine doses over varying durations and identified six response patterns, along with their associated biomolecules and pathways. Our study demonstrates the effectiveness of the proposed multi-omics data analysis in capturing multidimensional patterns of cellular response to chemical perturbation, enhancing understanding of pathway regulation for chemical risk assessment.

In-Silico Modelling of Transdermal Delivery of Macromolecule Drugs Assisted by a Skin Stretching Hypobaric Device.

OBJECTIVES: To develop a simulation model to explore the interplay between mechanical stretch and diffusion of large molecules into the skin under locally applied hypobaric pressure, a novel penetration enhancement method. METHODS: Finite element method was used to model the skin mechanical deformation and molecular diffusion processes, with validation against in-vitro transdermal permeation experiments. Simulations and experimental data were used together to investigate the transdermal permeation of large molecules under local hypobaric pressure. RESULTS: Mechanical simulations resulted in skin stretching and thinning (20%-26% hair follicle diameter increase, and 21%-27% skin thickness reduction). Concentration of dextrans in the stratum corneum was below detection limit with and without hypobaric pressure. Concentrations in viable epidermis and dermis were not affected by hypobaric pressure (approximately 2 µg [Formula: see text] cm-2). Permeation into the receptor fluid was substantially enhanced from below the detection limit at atmospheric pressure to up to 6 µg [Formula: see text] cm-2 under hypobaric pressure. The in-silico simulations compared satisfactorily with the experimental results at atmospheric conditions. Under hypobaric pressure, satisfactory comparison was attained when the diffusion coefficients of dextrans in the skin layers were increased from [Formula: see text] 10 µm2 [Formula: see text] s-1 to between 200-500 µm2 [Formula: see text] s-1. CONCLUSIONS: Application of hypobaric pressure induces skin mechanical stretching and enlarges the hair follicle. This enlargement alone cannot satisfactorily explain the increased transdermal permeation into the receptor fluid under hypobaric pressure. The results from the in-silico simulations suggest that the application of hypobaric pressure increases diffusion in the skin, which leads to improved overall transdermal permeation.

A regression analysis using simple descriptors for multiple dermal datasets: Going from individual membranes to the full skin.

In silico methods to estimate and/or quantify skin absorption of chemicals as a function of chemistry are needed to realistically predict pharmacological, occupational, and environmental exposures. The Potts-Guy equation is a well-established approach, using multi-linear regression analysis describing skin permeability (Kp) in terms of the octanol/water partition coefficient (logP) and molecular weight (MW). In this work, we obtained regression equations for different human datasets relevant to environmental and cosmetic chemicals. Since the Potts-Guy equation was published in 1992, we explored recent datasets that include different skin layers, such as dermatomed (including dermis to a defined thickness) and full skin. Our work was consistent with others who have observed that fits to the Potts-Guy equation are stronger for experiments focused on the epidermis. Permeability estimates for dermatomed skin and full skin resulted in low regression coefficients when compared to epidermis datasets. An updated regression equation uses a combination of fitted permeability values obtained with a published 2D compartmental model previously evaluated. The resulting regression equation was: logKp = -2.55 + 0.65logP - 0.0085MW, R2 = 0.91 (applicability domain for all datasets: MW ranges from 18 to >584 g/mol and -4 to >5 for logP). This approach demonstrates the advantage of combining mechanistic with structural activity relationships in a single modeling approach. This combination approach results in an improved regression fit when compared to permeability estimates obtained using the Potts-Guy approach alone. The analysis presented in this work assumes a one-compartment skin absorption route; future modeling work will consider adding multiple compartments.

Prediction of the Liquid-Liquid Extraction Properties of Imidazolium-Based Ionic Liquids for the Extraction of Aromatics from Aliphatics.

Liquid-liquid extraction (LLE) is an important technique to separate aromatics from aliphatics since these compounds have very similar boiling points and cannot be separated by distillation. Ionic liquids (ILs) are considered as potential extractants to extract aromatics from aliphatics. In this paper, molecular dynamics (MD) simulations were used to predict the extraction property (i.e., capacity and selectivity) of ILs for the LLE of aromatics from aliphatics. The extraction properties of seven different ILs including [C2mim][Tf2N], [C2mim][TFO], [C2mim][SCN], [C2mim][DCA], [C2mim][TCM], [C4mim][Tf2N], and [C8mim][Tf2N] were investigated. Results show that ILs with shorter alkyl chain cations and [Tf2N]- anion exhibit better extraction efficiency than other ILs, which is in agreement with previously reported experimental data on the extraction of toluene from aliphatics and further validated the reliability of the proposed model. The binding energies between ILs and organic molecules were calculated by the density functional theory, which help explain the different extraction behaviors of different ILs. The symmetry-adapted perturbation theory analysis was performed to further understand the interaction mechanisms between ILs and organics. Our study shows that the [Tf2N]- anion also has the best extraction capability for heavier aromatics (o-xylene, m-xylene, and p-xylene) from common aliphatics (heptane and octane). The MD modeling approach can be a low-cost in silico tool for the high-throughput fast screening of ILs for the LLE of aromatics from aliphatics.

Numerical analysis of the strain distribution in skin domes formed upon the application of hypobaric pressure.

BACKGROUND: Suction cups are widely used in applications such as in measurement of mechanical properties of skin in vivo, in drug delivery devices or in acupuncture treatment. Understanding mechanical response of skin under hypobaric pressure is of great importance for users of suction cups. The aim of this work is to predict the hypobaric pressure induced 3D stretching of the skin. METHODS: Experimental skin tensile tests were carried out for mechanical property characterization. Both linear elasticity and hyperelasticity parameters were determined and implemented in Finite Element modelling. Skin suction tests were performed in both experiments and FEM simulations for model validation. 3D skin stretching is then visualized in detail in FEM simulations. RESULTS: The simulations showed that the skin was compressed consistently along the thickness direction, leading to reduced thickness. At the center of the dome, the radial and angular strain decreases from the top surface to the bottom surface, although always in tension. Hyperelasticity modelling showed superiority over linear elasticity modelling while predicting the strain distribution because the stretch ratio reaches values exceeding the initial linear elastic stage of the stress-strain curve for skin. CONCLUSION: Hyperelasticity modelling is an effective approach to predict the 3D strain distribution, which paves a way to accurately design safe commercial products that interface with the skin.

In Silico Simulation of Simultaneous Percutaneous Absorption and Xenobiotic Metabolism: Model Development and a Case Study on Aromatic Amines.

PURPOSE: To advance physiologically-based pharmacokinetic modelling of xenobiotic metabolism by integrating metabolic kinetics with percutaneous absorption. METHOD: Kinetic rate equations were proposed to describe the metabolism of a network of reaction pathways following topical exposure and incorporated into the diffusion-partition equations of both xenobiotics and metabolites. The published ex vivo case study of aromatic amines was simulated. Diffusion and partition properties of xenobiotics and subsequent metabolites were determined using physiologically-based quantitative structure property relationships. Kinetic parameters of metabolic reactions were best fitted from published experimental data. RESULTS: For aromatic amines, the integrated transdermal permeation and metabolism model produced data closely matched by experimental results following limited parameter fitting of metabolism rate constants and vehicle:water partition coefficients. The simulation was able to produce dynamic concentration data for all the dermal layers, as well as the vehicle and receptor fluid. CONCLUSION: This mechanistic model advances the dermal in silico functionality. It provides improved quantitative spatial and temporal insight into exposure of xenobiotics, enabling the isolation of governing features of skin. It contributes to accurate modelling of concentrations of xenobiotics reaching systemic circulation and additional metabolite concentrations. This is vital for development of both pharmaceuticals and cosmetics.

In Silico Prediction of the Thermodynamic Equilibrium of Solute Partition in Multiphase Complex Fluids: A Case Study of Oil-Water Microemulsion.

Multiphase complex fluids such as micelles, microemulsions, and dispersions are ubiquitous in product formulations of foods, pharmaceuticals, cosmetics, and fine chemicals. Quantifying how active solutes partition in the microstructure of such multiphase fluids is necessary for designing formulations that can optimally deliver the benefits of functional actives. In this paper, we at first predict the structure of a heptane/butanol/sodium dodecyl sulfate droplet in water that self-assembled to form a microemulsion through the molecular dynamics (MD) simulation and subsequently investigate the thermodynamic equilibrium of solute partitioning using COSMOmic. To our knowledge, this is the first time that the MD/COSMOmic approach is used for predicting solute partitioning in a microemulsion. The predicted partition coefficients are compared to experimental values derived from retention measurements of the same microemulsion. We show that the experimental data of droplet-water partition coefficients (Kdroplet/w) can be reliably predicted by the method that combines MD simulations with COSMOmic.

Determining the Effect of pH on the Partitioning of Neutral, Cationic and Anionic Chemicals to Artificial Sebum: New Physicochemical Insight and QSPR Model.

PURPOSE: Sebum is an important shunt pathway for transdermal permeation and targeted delivery, but there have been limited studies on its permeation properties. Here we report a measurement and modelling study of solute partition to artificial sebum. METHODS: Equilibrium experiments were carried out for the sebum-water partition coefficients of 23 neutral, cationic and anionic compounds at different pH. RESULTS: Sebum-water partition coefficients not only depend on the hydrophobicity of the chemical but also on pH. As pH increases from 4.2 to 7.4, the partition of cationic chemicals to sebum increased rapidly. This appears to be due to increased electrostatic attraction between the cationic chemical and the fatty acids in sebum. Whereas for anionic chemicals, their sebum partition coefficients are negligibly small, which might result from their electrostatic repulsion to fatty acids. Increase in pH also resulted in a slight decrease of sebum partition of neutral chemicals. CONCLUSIONS: Based on the observed pH impact on the sebum-water partition of neutral, cationic and anionic compounds, a new quantitative structure-property relationship (QSPR) model has been proposed. This mathematical model considers the hydrophobic interaction and electrostatic interaction as the main mechanisms for the partition of neutral, cationic and anionic chemicals to sebum.

Development of a Two-Dimensional Model for Predicting Transdermal Permeation with the Follicular Pathway: Demonstration with a Caffeine Study.

PURPOSE: The development of a new two-dimensional (2D) model to predict follicular permeation, with integration into a recently reported multi-scale model of transdermal permeation is presented. METHODS: The follicular pathway is modelled by diffusion in sebum. The mass transfer and partition properties of solutes in lipid, corneocytes, viable dermis, dermis and systemic circulation are calculated as reported previously [Pharm Res 33 (2016) 1602]. The mass transfer and partition properties in sebum are collected from existing literature. None of the model input parameters was fit to the clinical data with which the model prediction is compared. RESULTS: The integrated model has been applied to predict the published clinical data of transdermal permeation of caffeine. The relative importance of the follicular pathway is analysed. Good agreement of the model prediction with the clinical data has been obtained. The simulation confirms that for caffeine the follicular route is important; the maximum bioavailable concentration of caffeine in systemic circulation with open hair follicles is predicted to be 20% higher than that when hair follicles are blocked. CONCLUSIONS: The follicular pathway contributes to not only short time fast penetration, but also the overall systemic bioavailability. With such in silico model, useful information can be obtained for caffeine disposition and localised delivery in lipid, corneocytes, viable dermis, dermis and the hair follicle. Such detailed information is difficult to obtain experimentally.

In Silico Modelling of Transdermal and Systemic Kinetics of Topically Applied Solutes: Model Development and Initial Validation for Transdermal Nicotine.

PURPOSE: The purpose was to develop a mechanistic mathematical model for predicting the pharmacokinetics of topically applied solutes penetrating through the skin and into the blood circulation. The model could be used to support the design of transdermal drug delivery systems and skin care products, and risk assessment of occupational or consumer exposure. METHODS: A recently reported skin penetration model [Pharm Res 32 (2015) 1779] was integrated with the kinetic equations for dermis-to-capillary transport and systemic circulation. All model parameters were determined separately from the molecular, microscopic and physiological bases, without fitting to the in vivo data to be predicted. Published clinical studies of nicotine were used for model demonstration. RESULTS: The predicted plasma kinetics is in good agreement with observed clinical data. The simulated two-dimensional concentration profile in the stratum corneum vividly illustrates the local sub-cellular disposition kinetics, including tortuous lipid pathway for diffusion and the "reservoir" effect of the corneocytes. CONCLUSIONS: A mechanistic model for predicting transdermal and systemic kinetics was developed and demonstrated with published clinical data. The integrated mechanistic approach has significantly extended the applicability of a recently reported microscopic skin penetration model by providing prediction of solute concentration in the blood.

In silico prediction of percutaneous absorption and disposition kinetics of chemicals.

PURPOSE: To develop in-silico model for predicting percutaneous absorption and disposition kinetics of chemicals in skin layers so as to facilitate the design of transdermal drug delivery systems and skin care products, and risk assessment of occupational or consumer exposure. METHODS: A general-purpose computer model for simulating skin permeation, absorption and disposition kinetics in the stratum corneum, viable dermis and dermis has been developed. Equations have been proposed for determining the partition and diffusion properties of chemicals by considering molecular partition, binding and mobility in skin layers. In vitro skin penetration data of 12 chemicals was used to validate the model. RESULTS: The observed and simulated permeation and disposition in skin layers were compared for 12 tested chemicals. For most tested chemicals, the experimental and model results are in good agreement with the coefficient of determination >0.80 and relative root mean squared error <1.20. The disposition kinetic parameters of the maximum concentration and the area under the curve in the viable epidermis and dermis initially increased with hydrophobicity, but reached maxima and then decreased with further increase of hydrophobicity. CONCLUSIONS: By considering skin physiological structure and composition, the partition and diffusion properties of chemicals in skin layers are determined. This allows in-silico simulation of percutaneous permeation, absorption and disposition kinetics of wide chemical space. The model produced results in good agreement with experimental data of 12 chemicals, suggesting a much improved framework to support transdermal delivery of drug and cosmetic actives as well as integrated risk assessment.

Free energy predictions of ligand binding to an α-helix using steered molecular dynamics and umbrella sampling simulations.

Free energy prediction of ligand binding to macromolecules using explicit solvent molecular dynamics (MD) simulations is computationally very expensive. Recently, we reported a linear correlation between the binding free energy obtained via umbrella sampling (US) versus the rupture force from steered molecular dynamics (SMD) simulations for epigallocatechin-3-gallate (EGCG) binding to α-helical-rich keratin. This linear correlation suggests a potential route for fast free energy predictions using SMD alone. In this work, the generality of the linear correlation is further tested for several ligands interacting with the α-helical motif of keratin. These molecules have significantly varying properties, i.e., octanol/water partition coefficient (log P), and/or overall charges (oleic acid, catechin, Fe(2+), citric acid, hydrogen citrate, dihydrogen citrate, and citrate). Using the constant loading rate of our previous study of the keratin-EGCG system, we observe that the linear correlation for keratin-EGCG can be extended to other uncharged molecules where interactions are governed by hydrogen bonds and/or a combination of hydrogen bonds and hydrophobic forces. For molecules where interactions with the keratin helix are governed primarily by electrostatics between charged molecules, a second, alternative linear correlation model is derived. While further investigations are needed to expand the molecular space and build a fully predictive model, the current approach represents a promising methodology for fast free energy predictions based on short SMD simulations (requiring picoseconds to nanoseconds of sampling) for defined biomolecular systems.

Progress in Topical and Transdermal Drug Delivery Research-Focus on Nanoformulations.

Skin is the largest organ and a multifunctional interface between the body and its environment. It acts as a barrier against cold, heat, injuries, infections, chemicals, radiations or other exogeneous factors, and it is also known as the mirror of the soul. The skin is involved in body temperature regulation by the storage of fat and water. It is an interesting tissue in regard to the local and transdermal application of active ingredients for prevention or treatment of pathological conditions. Topical and transdermal delivery is an emerging route of drug and cosmetic administration. It is beneficial for avoiding side effects and rapid metabolism. Many pharmaceutical, technological and cosmetic innovations have been described and patented recently in the field. In this review, the main features of skin morphology and physiology are presented and are being followed by the description of classical and novel nanoparticulate dermal and transdermal drug formulations. The biophysical aspects of the penetration of drugs and cosmetics into or across the dermal barrier and their investigation in diffusion chambers, skin-on-a-chip devices, high-throughput measuring systems or with advanced analytical techniques are also shown. The current knowledge about mathematical modeling of skin penetration and the future perspectives are briefly discussed in the end, all also involving nanoparticulated systems.

In Silico Prediction of Stratum Corneum Partition Coefficients via COSMOmic and Molecular Dynamics Simulations.

Stratum corneum (SC) is the main barrier of human skin where the inter-corneocytes lipids provide the main pathway for transdermal permeation of functional actives of skin care and health. Molecular dynamics (MD) has been increasingly used to simulate the SC lipid bilayer structure so that the barrier property and its affecting factors can be elucidated. Among reported MD simulation studies, solute partition in the SC lipids, an important parameter affecting SC permeability, has received limited attention. In this work, we combine MD simulation with COSMOmic to predict the partition coefficients of dermatologically relevant solutes in SC lipid bilayer. Firstly, we run MD simulations to obtain equilibrated SC lipid bilayers with different lipid types, compositions, and structures. Then, the simulated SC lipid bilayer structures are fed to COSMOmic to calculate the partition coefficients of the solutes. The results show that lipid types and bilayer geometries play a minor role in the predicted partition coefficients. For the more lipophilic solutes, the predicted results of solute partition in SC lipid bilayers agree well with reported experimental values of solute partition in extracted SC lipids. For the more hydrophilic molecules, there is a systematical underprediction. Nevertheless, the MD/COSMOmic approach correctly reproduces the phenomenological correlation between the SC lipid/water partition coefficients and the octanol/water partition coefficients. Overall, the results show that the MD/COSMOmic approach is a fast and valid method for predicting solute partitioning into SC lipids and hence supporting the assessment of percutaneous absorption of skin care ingredients, dermatological drugs as well as environmental pollutants.

An examination of published datasets of skin permeability and partition coefficients.

Permeability and partition coefficients of the skin barrier are important for assessing dermal absorption, bioavailability, and safety of cosmetics and medicine. We use the Potts and Guy equation to analyse the dependence of skin permeability on the hydrophobicity of permeants and highlight the significant differences in published datasets. Correlations of solute partition to skin are examined to understand the likely causes of the differences in the skin permeability datasets. Recently published permeability datasets show weak correlation and low dependence on hydrophobicity. As expected, early datasets show good correlation with hydrophobicity due to the related derivation. The weaker correlation of later datasets cannot be explained by the partition to skin lipids. All the datasets of solute partition to skin lipid showed a similar correlation to hydrophobicity where the log-linear correlation coefficient of partition is almost the same of the log-linear coefficient of Potts and Guy equation. Weak correlation and dependence of the late permeability datasets with SC lipid/water partition and that they are significantly under predicted by the Potts and Guy equation suggests either additional non-lipid pathway at play or a weaker skin barrier property.

Evaluation of Constrained and Restrained Molecular Dynamics Simulation Methods for Predicting Skin Lipid Permeability.

Recently, molecular dynamics (MD) simulations have been utilized to investigate the barrier properties of human skin stratum corneum (SC) lipid bilayers. Different MD methods and force fields have been utilized, with predicted permeabilities varying by few orders of magnitude. In this work, we compare constrained MD simulations with restrained MD simulations to obtain the potential of the mean force and the diffusion coefficient profile for the case of a water molecule permeating across an SC lipid bilayer. Corresponding permeabilities of the simulated lipid bilayer are calculated via the inhomogeneous solubility diffusion model. Results show that both methods perform similarly, but restrained MD simulations have proven to be the more robust approach for predicting the potential of the mean force profile. Critical to both methods are the sampling of the whole trans-bilayer axis and the following symmetrization process. Re-analysis of the previously reported free energy profiles showed that some of the discrepancies in the reported permeability values is due to misquotation of units, while some are due to the inaccurately obtained potential of the mean force. By using the existing microscopic geometrical models via the intercellular lipid pathway, the permeation through the whole SC is predicted from the MD simulation results, and the predicted barrier properties have been compared to experimental data from the literature with good agreement.

Inhibition of the intestinal postprandial glucose transport by gallic acid and gallic acid derivatives.

Inhibition of glucose uptake in the intestine through sodium-dependent glucose transporter 1 (SGLT1) or glucose transporter 2 (GLUT2) may be beneficial in controlling postprandial blood glucose levels. Gallic acid and ten of its derivatives were identified in the active fractions of Terminalia chebula Retz. fructus immaturus, a popular edible plant fruit which has previously been associated with the inhibition of glucose uptake. Gallic acid derivatives (methyl gallate, ethyl gallate, pentyl gallate, 3,4,6-tri-O-galloyl-ß-d-glucose, and corilagin) showed good glucose transport inhibition with inhibitory rates of 72.1 ± 1.6%, 71.5 ± 1.4%, 79.9 ± 1.2%, 44.7 ± 1.2%, and 75.0 ± 0.7% at 5 mM d-glucose and/or 56.3 ± 2.3, 52.1 ± 3.2%, 70.2 ± 1.7%, 15.6 ± 1.6%, and 37.1 ± 0.8% at 25 mM d-glucose. However, only 3,4,6-tri-O-galloyl-ß-d-glucose and corilagin were confirmed GLUT2-specific inhibitors. Whilst some tea flavonoids demonstrated minimal glucose transport inhibition, their gallic acid derivatives strongly inhibited transport effect with GLUT2 specificity. This suggests that gallic acid structures are crucial for glucose transport inhibition. Plants, such as T. chebula, which contain high levels of gallic acid and its derivatives, show promise as natural functional ingredients for inclusion in foods and drinks designed to control postprandial glucose levels.

Predicting Partition Coefficients of Neutral and Charged Solutes in the Mixed SLES-Fatty Acid Micellar System.

Sodium laureth sulfate (SLES) and fatty acids are common ingredients in many cosmetic products. Understanding how neutral and charged fatty acid compounds partition between micellar and water phases is crucial to achieve the optimal design of the product formulation. In this paper, we first study the formation of mixed SLES and fatty acid micelles using molecular dynamics (MD) simulations. Micelle/water partition coefficients of neutral and charged fatty acids are then calculated using COSMOmic as well as a MD approach based on the potential of mean force (PMF) calculations performed using umbrella sampling (US). The combined US/PMF approach was performed with both the additive, non-polarizable CHARMM general force field (CGenFF) and the classical Drude polarizable force field. The partition coefficients for the neutral solutes are shown to be accurately calculated with the COSMOmic and additive CGenFF US/PMF approaches, while only the US/PMF approach with the Drude polarizable force field accurately calculated the experimental partition coefficient of the charged solute. These results indicate the utility of the Drude polarizable force field as a tool for the rational development of mixed micelles.