Abstract
The oxidation of human sebum, a lipid mixture covering our skin, generates a range of volatile and semivolatile 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 nontrivial, 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.