Stratum Corneum (SC), the outermost layer of skin, is mainly responsible for skin's barrier function. The complex lipid matrix of SC determines these barrier properties. In this study, the lipid matrix is modeled as an equimolar mixture of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA). The permeation of water, oxygen, ethanol, acetic acid, urea, butanol, benzene, dimethyl sulfoxide (DMSO), toluene, phenol, styrene, and ethylbenzene across this layer is studied using a constrained MD simulations technique. Several long constrained simulations are performed at a skin temperature of 310 K under NPT conditions. The free energy profiles and diffusion coefficients along the bilayer normal have been calculated for each molecule. Permeability coefficients are also calculated and compared with experimental data. The main resistance for the permeation of hydrophilic and hydrophobic permeants has been found to be in the interior of the lipid bilayer and near the lipid-water interface, respectively. The obtained permeability is found to be a few orders of magnitude higher than experimental values for hydrophilic molecules while for hydrophobic molecules more discrepancy was observed. Overall, the qualitative ranking is consistent with the experiments.

Accurate in-silico models are required to predict the release of drug molecules through skin in order to supplement the in-vivo experiments for faster development/testing of drugs. The upper most layer of the skin, stratum corneum (SC), offers the main resistance for permeation of actives. Most of the SC's molecular level models comprise cholesterol and phospholipids only, which is far from reality. In this study we have implemented a multiscale modeling framework to obtain the release profile of three drugs, namely, caffeine, fentanyl, and naphthol, through skin SC. We report for the first time diffusion of drugs through a realistic skin molecular model comprised of ceramides, cholesterol, and free fatty acid. The diffusion coefficients of drugs in the SC lipid matrix were determined from multiple constrained molecular dynamics simulations. The calculated diffusion coefficients were then used in the macroscopic models to predict the release profiles of drugs through the SC. The obtained release profiles were in good agreement with available experimental data. The partition coefficient exhibits a greater effect on the release profiles. The reported multiscale modeling framework would provide insight into the delivery mechanisms of the drugs through the skin and shall act as a guiding tool in performing targeted experiments to come up with a suitable delivery system.

The PMF between two nanoparticles (of any chemistry and coverage) can be used as a measure to predict their stability.