Insights into Dermal Permeation of Skin Oil Oxidation Products from Enhanced Sampling Molecular Dynamics Simulation

TL;DR

This study uses advanced molecular dynamics simulations to analyze how skin oil oxidation products permeate the skin barrier, revealing significant sampling challenges and errors affecting permeability estimates.

Contribution

It demonstrates the necessity of extensive sampling in simulations to accurately assess skin permeation of oxidation products, highlighting limitations of current methods.

Findings

Simulations require microsecond timescales and enhanced sampling.

Structural order and asymmetries cause slow free energy convergence.

Sampling errors significantly impact permeability calculations.

Abstract

The oxidation of human sebum, a lipid mixture covering our skin, generates a range of volatile and semi-volatile carbonyl compounds that contribute largely to indoor air pollution in crowded environments. Kinetic models have been developed to gain a deeper understanding of this complex multiphase chemistry, but they rely partially on rough estimates of kinetic and thermodynamic parameters, especially those describing skin permeation. Here, we employ atomistic molecular dynamics simulations to study the translocation of selected skin oil oxidation products through a model stratum corneum membrane. We find these simulations to be non-trivial, requiring extensive sampling with up to microsecond simulation times, in spite of employing enhanced sampling techniques. We identify the high degree of order and stochastic, long-lived temporal asymmetries in the membrane structure as the leading causes for the slow convergence of the free energy computations. We demonstrate that statistical errors due to insufficient sampling are substantial and propagate to membrane permeabilities. These errors are independent of the enhanced sampling technique employed and very likely independent of the precise membrane model.