Pith. sign in

REVIEW 4 major objections 4 minor 77 references

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

T0 review · 4 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper shows that free-energy profiles for permeation of skin oil oxidation products through a model stratum corneum converge only on microsecond time scales, and that the resulting sampling errors propagate to membrane permeabilities…

desk verdict Solid, useful warning about slow convergence in SC permeation PMFs, but the lipid-flip mechanism is plausible rather than proven. read the letter →

arxiv 2412.00436 v2 pith:WLEE4FWD submitted 2024-11-30 cond-mat.soft physics.bio-phphysics.chem-ph

Rinto Thomas , Praveen Ranganath Prabhakar , Douglas J. Tobias , Michael von Domaros This is my paper · ORCID
classification cond-mat.softphysics.bio-phphysics.chem-ph
keywords stratumcorneummoleculardynamicsfreeenergyprofileumbrellasamplingwell-temperedmetadynamicslipidflippingskinpermeationindoorairchemistry
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to establish that computing how skin oil oxidation products pass through the outermost skin layer is far more demanding than typical membrane permeation simulations: the free-energy profiles take multiple microseconds to converge even when enhanced sampling is used, and insufficient sampling produces statistical errors that can change predicted permeabilities by up to an order of magnitude. The authors argue that this slow convergence is not a flaw of a particular enhanced-sampling method but a property of the stratum corneum model itself, caused by its rigid, highly ordered lipid matrix and by rare, long-lived flips of cholesterol and free fatty acids that break the membrane's symmetry for microseconds at a time. If this is right, it matters for indoor air chemistry, where kinetic models of skin-oil ozonolysis currently rely on roughly estimated permeation parameters; the simulations can supply those parameters with uncertainties that are acceptable for such models, but only if convergence is checked rigorously. The study therefore carries a practical warning: sampling only half the membrane or force-symmetrizing the profile too early can hide severe errors.

What carries the argument

The argument is carried by two enhanced-sampling techniques applied to a single collective variable, the solute's position $z$ along the membrane normal: umbrella sampling with harmonic restraints and well-tempered metadynamics, both analyzed with force-symmetrized profiles and blocking-based error estimates. The central mechanistic object is the rare, stochastic flipping of cholesterol and free fatty acid lipids, tracked through tilt-angle time series and lipid-solute coordination numbers; these flips create long-lived structural asymmetries that slow free-energy convergence. The connection to permeability is made through the inhomogeneous solubility-diffusion model, which the authors reduce to $P/D = \left[\int e^{\beta \Delta F(z)}\,dz\right]^{-1}$ by assuming a uniform diffusivity, so that errors in the free-energy profile translate exponentially into errors in permeability.

What would settle it

Run the same umbrella-sampling calculation in a stratum corneum model in which cholesterol and free fatty acid flipping is prevented, for example by restraining lipid tilt angles: if the free-energy profile still takes microseconds to converge and still shows the same statistical errors, then lipid flipping is not the leading cause of the slow convergence.

Watch

Extended reading notes

Core claim

The central discovery is empirical and methodological: for a fully hydrated stratum corneum lipid bilayer, potential-of-mean-force profiles for water, acetone, and 6-methyl-5-hepten-2-one converge to their symmetric limits only on microsecond time scales, and both umbrella sampling and well-tempered metadynamics suffer from this bottleneck. The authors show that the slow convergence is tied to the membrane's own slow reorganization: cholesterol and free fatty acid molecules occasionally flip across the bilayer and remain in flipped states for up to microseconds, creating persistent asymmetries that bias the free energy. These sampling errors propagate through the inhomogeneous solubility-diffusion relation into membrane permeabilities; in the worst case examined (acetone with well-tempered metadynamics) the uncertainty spans roughly an order of magnitude. They conclude that the statistical errors are independent of the enhanced-sampling technique and very likely independent of the precise membrane model, and that earlier practice of sampling half the system or symmetrizing prematurely can conceal the problem.

Load-bearing premise

The paper's explanation for why sampling is so slow depends on a small number of visually identified lipid-flipping events; if those flips are an artifact of the computer model rather than genuine behaviour of skin lipids, the mechanism would not generalize even if the slow convergence itself is real.

Editorial extensions

If this is right

  • Future permeation simulations of stratum corneum models should not assume rapid convergence or sample only half the bilayer; force-symmetrized free-energy profiles must be checked for time convergence.
  • For the solutes studied, the converged profiles span orders of magnitude in reduced permeability, with water and 6-MHO agreeing between umbrella sampling and well-tempered metadynamics, while acetone differs substantially, indicating it needs even longer sampling.
  • When computational cost is not a concern, umbrella sampling appears to converge more reliably for permeability in these systems than well-tempered metadynamics, at an 8- to 15-fold increase in total simulation time.
  • The observed order-of-magnitude permeability uncertainty is judged acceptable for indoor-air kinetic models, so the molecular dynamics approach remains a viable source of the transport parameters those models require.
  • The paper recommends that any future permeation study verify convergence of the full free-energy profile, since eliminating asymmetries through premature symmetrization can lead to severe errors.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Inference: If the observed lipid flips are an artifact of the CHARMM36 force field rather than a genuine property of stratum corneum lipids, the mechanistic explanation would not generalize, even though the empirical slow-convergence finding might survive in other atomistic models; a direct test would compare convergence in simulations with flipping artificially suppressed.
  • Inference: Because flipping is rare, stochastic, and long-lived, biasing only the solute's position may be insufficient; efficient convergence may require biasing membrane degrees of freedom such as lipid tilt or leaflet asymmetry, an approach the paper notes is nontrivial.
  • Inference: The same slow-convergence problem may affect other highly ordered multicomponent lipid models, including long-periodicity-phase stratum corneum models, though the paper only tests a short-periodicity-phase bilayer and a POPC control membrane.
  • Inference: A practical consequence for kinetic modeling is that computed permeabilities for volatile products like acetone should be reported with error bars of roughly an order of magnitude, and single-number estimates from short simulations should be treated with caution.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. This manuscript reports atomistic MD simulations of water, acetone, and 6-MHO permeating a fully hydrated model stratum corneum (SC) short-periodicity-phase bilayer (the Wang–Klauda model) and, for comparison, a POPC bilayer. Free-energy profiles are computed using umbrella sampling (US) and well-tempered metadynamics (WTM), and reduced permeabilities P/D are obtained from the inhomogeneous solubility–diffusion model. The central empirical finding is that convergence of the FE profiles to the symmetry required of a symmetric bilayer is extremely slow for the SC model—requiring tens of microseconds of US sampling and multi-microsecond WTM runs—whereas POPC converges much faster. The paper attributes this slow convergence to lipid-flip-induced long-lived leaflet asymmetries and argues that insufficient sampling produces substantial errors that propagate to permeabilities.

Significance. The empirical convergence analysis is valuable and methodologically important. The use of two enhanced-sampling methods, two membrane models, and three solutes, with consistent trends across most systems, provides credible evidence that ordered SC-like bilayers are a hard case for permeation PMF calculations. The explicit propagation of FE uncertainty into P/D, the check against literature water permeabilities, and the warning against half-system symmetrization are useful and actionable. However, the paper's stronger claims—that lipid flipping is the leading cause and that the sampling bottleneck is independent of the membrane model—go beyond what the presented evidence supports; these should be reframed as hypotheses or supported with targeted analysis.

major comments (4)
  1. [Long-lived asymmetries arise due to lipid flipping / Table S1] The causal claim that stochastic lipid flipping is a leading cause of the slow FE convergence is not quantitatively supported. The authors state that 'despite µs simulation times, flipping events were found to be rare and stochastic, and not enough data could be gathered to determine rates,' and the evidence base in Table S1 is about 18–19 events, with no correlation established between the timing, duration, or leaflet location of flips and the magnitude of leaflet asymmetry in the evolving FE profiles. Many listed flips last only a few nanoseconds to a few hundred nanoseconds, while the convergence timescale is microseconds. The paper should either soften this to a plausible hypothesis or provide direct evidence (e.g., conditional asymmetry metrics or replicate runs). This also affects the abstract's generalization to other membrane models.
  2. [Abstract / Conclusions] The abstract's statement that the statistical errors 'are independent of the enhanced sampling technique employed and very likely independent of the precise membrane model' is not supported by the data. The WTM SC/Acetone force-symmetrized profile has not converged by 6.7 µs (Figure 4, bottom row), and the resulting reduced permeability differs substantially from the US result (Figures 12 and 13). The membrane-model independence rests on a single comparison between one SC model and one POPC model. The data support the narrower conclusion that slow convergence was observed with both methods in this SC model and was less severe in POPC; the stronger generalizations should be removed or explicitly labeled as speculation.
  3. [Table 1 and Computational Details] Table 1 is internally inconsistent. The header says all simulation times are in µs, and the US per-window column lists 150 for POPC. With 65 windows at 150 ps per window, the total is 9.75 ns, not 9.8 µs, and the Computational Details text says '150ps per window' while Figure S3 axis labels are in ns. This inconsistency makes the quantitative SC-vs-POPC comparison—which is central to the claim that slow convergence is a feature of the SC model—ambiguous. Please correct the units and ensure the table, text, and figures agree.
  4. [Results (Permeabilities) / Figures 12 and 13] After stating that the WTM force-symmetrized FE profiles have not converged for water and acetone (Figure 4, bottom row), the paper nonetheless uses the final WTM acetone profile to compute a reduced permeability and then interprets the US/WTM acetone discrepancy as a possible methodological difference. Because the WTM acetone simulation is explicitly unconverged, the discrepancy cannot be attributed to the sampling technique itself; it may be residual sampling error. The comparison should be restricted to converged quantities, or the acetone WTM result should be labeled as an unconverged estimate.
minor comments (4)
  1. [Supporting Information] The phrase 'A vail able' appears where 'available' is intended; please fix this typo in the Supporting Information section.
  2. [Computational Details] There are several missing spaces and inconsistent formatting, e.g., 'for10ns', 'of2Å', '50000', and '100fs'. A light copyedit would improve readability.
  3. [Figure 5 and Figure S5] The statement that lines and error bars have been 'slightly smoothed' is vague. Since the barrier heights in Table 2 are read from these profiles, please describe the smoothing procedure quantitatively or state that the unsmoothed values do not differ beyond the reported confidence intervals.
  4. [Conclusions] The concluding paragraph says the largest error interval (acetone, WTM) is 'likely good enough' for indoor-air kinetic modeling. This is reasonable, but it would be helpful to state explicitly that this judgment applies to the reduced permeability P/D rather than to a full permeability coefficient, since D(z) is set constant in this work.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the free-energy profiles and permeabilities are simulation outputs, not fitted inputs.

full rationale

The paper's central deliverables, the free-energy profiles and reduced permeabilities, are simulation outputs processed through the inhomogeneous solubility-diffusion model (Eqs. 3 and 4). No parameter is fitted to a target permeability, and no literature permeability is used as an input; the agreement of the water permeabilities in SC and POPC with prior studies is an external consistency check rather than a fitted constraint. The bilayer-symmetry criterion is used as a physical consistency condition: force-symmetrization averages the two independently sampled leaflets, and convergence is additionally assessed from the time evolution of the raw and symmetrized profiles and of the permeabilities (Figs. 3, 4, and 13). The symmetry constraint therefore does not by construction produce the reported values. The mechanistic claim that lipid flipping drives slow convergence is inferred from trajectory inspection, tilt-angle time series, and contact analysis; although the evidence base is limited, with the authors noting that flipping rates could not be determined from the scarce events catalogued in Table S1, this is an evidentiary weakness rather than a circular reduction. No load-bearing self-citation or imported uniqueness theorem appears; the self-citations provide context but do not carry the derivation. No specific circular step can be exhibited in which a quoted equation or fitted parameter reduces to the claimed result.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claims rest on no fitted constants. The membrane model, force field, and the constant-diffusivity reduction are inputs from prior literature or simplifying choices, and the error estimation uses established methods. The only new interpretive entity, lipid flipping, is an observed trajectory feature, not a postulated degree of freedom.

assumptions (4)
  • domain assumption The Wang-Klauda SPP bilayer is a representative model of the stratum corneum lipid barrier.
    The central claims about SC permeation are made for this specific three-component model; the paper notes it is 'rather simplistic' but argues the convergence problems very likely generalize. Section: Computational Details, Model system and force field.
  • domain assumption CHARMM36/CGenFF force field parameters, including auto-assigned parameters for 6-MHO with a charge penalty of 6.377, accurately represent solute-lipid interactions.
    Permeation free energy barriers and the observed lipid flipping behavior are entirely determined by the force field; the 6-MHO parameters were auto-generated and not validated against experiment or high-level quantum chemistry. Section: Computational Details, force field paragraph.
  • domain assumption Reduced permeabilities P/D computed with a constant diffusivity profile D(z)=D capture the thermodynamic contribution to permeation.
    The authors explicitly assume D(z)=D to isolate thermodynamic effects (Equation 4), so the reported P/D values are not full permeabilities and ignore position-dependent diffusivity. Section: Results, Permeabilities are sensitive to free energy differences.
  • standard math WHAM and Jonsson automated blocking provide valid error estimates for the free energy profiles.
    Error bars from Zhu-Hummer propagation and blocking are standard tools; they are used to define the 95% confidence intervals that drive the convergence conclusions. Section: Computational Details and SI section 1.1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Insights into Dermal Permeation of Skin Oil Oxidation Products from Enhanced Sampling Molecular Dynamics Simulation." pith.science (2026). https://pith.science/paper/WLEE4FWD

@misc{pith2026241200436,
  author       = {Pith},
  title        = {Pith review of: Insights into Dermal Permeation of Skin Oil Oxidation Products from Enhanced Sampling Molecular Dynamics Simulation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WLEE4FWD}},
  note         = {Machine review of arXiv:2412.00436}
}
read the original 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.

Figures

Figures reproduced from arXiv: 2412.00436 by the authors.

Figure 1
Figure 1. Squalene and its volatile and semivolatile ozonolysis products. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Schematic representation of a perme￾ation path through the SC model employed in this study. Lipid carbon atoms are colored or￾ange (CHL), green (CER) and yellow (FFA). Oxygen atoms are shown in red and nitrogen in blue. The water phase is represented as a trans￾parent surface. The collective variable (CV) is the projection of the path onto the membrane normal, with z = 0 corresponding to the center of the membrane. … view at source ↗
Figure 3
Figure 3. Top row: Convergence of FE profiles toward symmetry for US. Bottom row: Convergence [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (10 more)
Figure 4
Figure 4. Figure 4: Top row: Convergence of FE profiles toward symmetry for WTM. Bottom row: Convergence [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: Final, force-symmetrized FE profiles for all investigated solutes in the SC model. Lines and error bars (representing 95 % con￾fidence intervals) have been slightly smoothed for visual clarity. The raw data is shown in Figure S5 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: Carbon–hydrogen order parameters for the acyl and free fatty acid chains found in SC [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: Snapshots of a tagged cholesterol molecule flipping to the opposite layer of the SC membrane in one of the US simulations. Visual inspection of the trajectories revealed the occasional flipping of lipids, that is, the de￾parture of the lipids from their average positio…
Figure 8
Figure 8. Figure 8: Unit vectors describing molecular ori￾entation are placed along the dashed lines for CHL (a), FFA (b), and CER (c). Only heavy atoms are shown in the structural representa￾tions. Despite µs simulation times, flipping events were found to be rare and stochastic, and not…
Figure 10
Figure 10. Figure 10: Coordination number between SC lipids and acetone, stratified by the position along the membrane normal. A similar trend was observed for the other solutes (Figure S9). cholesterol and free fatty acid molecules, on the other hand, we observe a significant degree of co…
Figure 9
Figure 9. Figure 9: Exemplary tilt angle time series for FFA and CHL in one of the WTM simulations of the SC membrane. Contact analysis reveals a significant de￾gree of lipid flipping To further characterize the flipping of the lipids, we calculated the number of lipids coordinating to a …
Figure 11
Figure 11. Figure 11: Schematic representation of an acetone-driven flipping of a free fatty acid lipid (a), the [PITH_FULL_IMAGE:figures/full_fig_p012_11.png]
Figure 12
Figure 12. Figure 12: Reduced permeability P/D for all investigated solutes in the SC membrane. Error bars represent 95 % confidence intervals. Force￾symmetrized FE profiles ( [PITH_FULL_IMAGE:figures/full_fig_p012_12.png]
Figure 13
Figure 13. Figure 13: Reduced permeability P/D as a function of simulation time. Top row: US, bottom row: WTM. Force-symmetrized FE profiles were used in the computation. US simulation times are per window. Total simulation times can be obtained by multiplying with the number of windows (N…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

77 extracted references · 74 canonical work pages

  1. [1]

    Weschler, C. J. Chemistry in Indoor Environments: 20 Years of Research. Indoor Air 2011, 21, 205--218

    work page 2011

  2. [2]

    J.; Carslaw, N

    Weschler, C. J.; Carslaw, N. Indoor Chemistry. Environ. Sci. Technol. 2018, 52, 2419--2428

    work page 2018

  3. [3]

    The Past, Present, and Future of Indoor Air Chemistry

    Bek \"o , G.; Carslaw, N.; Fauser, P.; Kauneliene, V.; Nehr, S.; Phillips, G.; Saraga, D.; Schoemaecker, C.; Wierzbicka, A.; Querol, X. The Past, Present, and Future of Indoor Air Chemistry. Indoor Air 2020, 30, 373--376

    work page 2020

  4. [4]

    Abbatt, J. P. D.; Wang, C. The Atmospheric Chemistry of Indoor Environments. Environ. Sci. Processes Impacts 2020, 22, 25--48

    work page 2020

  5. [5]

    G.; Dudzinska, M.; Wiesen, P.; Nehr, S.; Huttunen, K.; Querol, X

    Carslaw, N.; Bek \"o , G.; Langer, S.; Schoemaecker, C.; Mihucz, V. G.; Dudzinska, M.; Wiesen, P.; Nehr, S.; Huttunen, K.; Querol, X. et al. A New Framework for Indoor Air Chemistry Measurements: Towards a Better Understanding of Indoor Air Pollution. Indoor Environ. 2024, 1, 100001

    work page 2024

  6. [6]

    P.; Grassian, V

    Ault, A. P.; Grassian, V. H.; Carslaw, N.; Collins, D. B.; Destaillats, H.; Donaldson, D. J.; Farmer, D. K.; Jimenez, J. L.; McNeill, V. F.; Morrison, G. C. et al. Indoor Surface Chemistry: Developing a Molecular Picture of Reactions on Indoor Interfaces. Chem 2020, 6, 3203--3218

    work page 2020

  7. [7]

    von Domaros , M.; Tobias, D. J. Molecular Dynamics Simulations of the Interactions of Organic Compounds at Indoor Relevant Surfaces. Annu. Rev. Phys. Chem. 2025, in press

    work page 2025

  8. [8]

    Farmer, D.; E

    K. Farmer, D.; E. Vance, M. Indoor Air: Sources, Chemistry and Health Effects. Environ. Sci. Processes Impacts 2019, 21, 1227--1228

    work page 2019

Show all 77 references
  1. [9]

    How Do Breath and Skin Emissions Impact Indoor Air Chemistry? Indoor Air 2019, 29, 369--379

    Kruza, M.; Carslaw, N. How Do Breath and Skin Emissions Impact Indoor Air Chemistry? Indoor Air 2019, 29, 369--379

    2019

  2. [10]

    J.; Wisthaler, A.; Cowlin, S.; Tamás, G.; Strøm-Tejsen, P.; Hodgson, A

    Weschler, C. J.; Wisthaler, A.; Cowlin, S.; Tamás, G.; Strøm-Tejsen, P.; Hodgson, A. T.; Destaillats, H.; Herrington, J.; Zhang, J. J.; Nazaroff, W. W. Ozone-Initiated Chemistry in an Occupied Simulated Aircraft Cabin. Environ. Sci. Technol. 2007, 41, 6177--6184

    2007

  3. [11]

    Wisthaler, A.; Weschler, C. J. Reactions of Ozone with Human Skin Lipids: Sources of Carbonyls, Dicarbonyls, and Hydroxycarbonyls in Indoor Air. Proc. Natl. Acad. Sci. 2010, 107, 6568--6575

    2010

  4. [12]

    J.; Nazaroff, W

    Weschler, C. J.; Nazaroff, W. W. Human Skin Oil: A Major Ozone Reactant Indoors. Environ. Sci. Atmos. 2023, 3, 640--661

    2023

  5. [13]

    E.; Franko, J.; Jackson, L

    Anderson, S. E.; Franko, J.; Jackson, L. G.; Wells, J. R.; Ham, J. E.; Meade, B. J. Irritancy and Allergic Responses Induced by Exposure to the Indoor Air Chemical 4-Oxopentanal. Toxicol. Sci. 2012, 127, 371--381

    2012

  6. [14]

    T.; Hammer, M.; Kofoed-S rensen , V.; Clausen, P

    Wolkoff, P.; Larsen, S. T.; Hammer, M.; Kofoed-S rensen , V.; Clausen, P. A.; Nielsen, G. D. Human Reference Values for Acute Airway Effects of Five Common Ozone-Initiated Terpene Reaction Products in Indoor Air. Toxicol. Lett. 2013, 216, 54--64

    2013

  7. [15]

    G.; Weschler, C

    Bek \"o , G.; Allen, J. G.; Weschler, C. J.; Vallarino, J.; Spengler, J. D. Impact of Cabin Ozone Concentrations on Passenger Reported Symptoms in Commercial Aircraft. PLOS ONE 2015, 10, e0128454

    2015

  8. [16]

    Indoor Air Humidity, Air Quality, and Health -- an Overview

    Wolkoff, P. Indoor Air Humidity, Air Quality, and Health -- an Overview. Int. J. Hyg. Envir. Heal. 2018, 221, 376--390

    2018

  9. [17]

    Lakey, P. S. J.; Wisthaler, A.; Berkemeier, T.; Mikoviny, T.; P \"o schl, U.; Shiraiwa, M. Chemical Kinetics of Multiphase Reactions between Ozone and Human Skin Lipids: Implications for Indoor Air Quality and Health Effects. Indoor Air 2017, 27, 816--828

    2017

  10. [18]

    Lakey, P. S. J.; Morrison, G. C.; Won, Y.; Parry, K. M.; von Domaros , M.; Tobias, D. J.; Rim, D.; Shiraiwa, M. The Impact of Clothing on Ozone and Squalene Ozonolysis Products in Indoor Environments. Comm. Chem. 2019, 2, 1--8

    2019

  11. [19]

    B.; Nitsche, J

    Wang, T.-F.; Kasting, G. B.; Nitsche, J. M. A Multiphase Microscopic Diffusion Model for Stratum Corneum Permeability. II . Estimation of Physicochemical Parameters, and Application to a Large Permeability Database. J. Pharm. Sci. 2007, 96, 3024--3051

    2007

  12. [20]

    von Domaros , M.; Lakey, P. S. J.; Shiraiwa, M.; Tobias, D. J. Multiscale Modeling of Human Skin Oil-Induced Indoor Air Chemistry: Combining Kinetic Models and Molecular Dynamics. J. Phys. Chem. B 2020, 124, 3836--3843

    2020

  13. [21]

    J.; Waring, M

    Shiraiwa, M.; Carslaw, N.; Tobias, D. J.; Waring, M. S.; Rim, D.; Morrison, G.; Lakey, P. S. J.; Kruza, M.; von Domaros, M.; Cummings, B. E. et al. Modelling Consortium for Chemistry of Indoor Environments ( MOCCIE ): Integrating Chemical Processes from Molecular to Room Scale...

    2019

  14. [22]

    M.; Traynor, M

    Pecoraro, B.; Tutone, M.; Hoffman, E.; Hutter, V.; Almerico, A. M.; Traynor, M. Predicting Skin Permeability by Means of Computational Approaches: Reliability and Caveats in Pharmaceutical Studies. J. Chem. Inf. Model. 2019, 59, 1759--1771

    2019

  15. [23]

    Dissecting the Formation, Structure and Barrier Function of the Stratum Corneum

    Matsui, T.; Amagai, M. Dissecting the Formation, Structure and Barrier Function of the Stratum Corneum. Int. Immunol. 2015, 27, 269--280

    2015

  16. [24]

    Determination of Stratum Corneum Lipid Profile by Tape Stripping in Combination with High-Performance Thin-Layer Chromatography

    Weerheim, A.; Ponec, M. Determination of Stratum Corneum Lipid Profile by Tape Stripping in Combination with High-Performance Thin-Layer Chromatography. Arch. Dermatol. Res. 2001, 293, 191--199

    2001

  17. [25]

    A.; Gooris, G

    Bouwstra, J. A.; Gooris, G. S.; van der Spek , J. A.; Bras, W. Structural Investigations of Human Stratum Corneum by Small-Angle X-Ray Scattering. J. Invest. Dermatol. 1991, 97, 1005--1012

    1991

  18. [26]

    Schmitt, T.; Neubert, R. H. State of the Art in Stratum Corneum Research. Part II : Hypothetical Stratum Corneum Lipid Matrix Models. Skin Pharmacol. Physi. 2020, 33, 213--230

    2020

  19. [27]

    A.; N a d a ban, A.; Bras, W.; MCabe, C.; Bunge, A.; Gooris, G

    Bouwstra, J. A.; N a d a ban, A.; Bras, W.; MCabe, C.; Bunge, A.; Gooris, G. S. The Skin Barrier: An Extraordinary Interface with an Exceptional Lipid Organization. Prog. Lipid Res. 2023, 92, 101252

    2023

  20. [28]

    R.; Gooris, G

    Shamaprasad, P.; N a d a ban, A.; Iacovella, C. R.; Gooris, G. S.; Bunge, A. L.; Bouwstra, J. A.; McCabe, C. The Phase Behavior of Skin-Barrier Lipids: A Combined Approach of Experiments and Simulations. Biophys. J. 2024, 123, 3188--3204

    2024

  21. [29]

    C.; Xia, D.; Iacovella, C

    Shamaprasad, P.; Moore, T. C.; Xia, D.; Iacovella, C. R.; Bunge, A. L.; McCabe, C. Multiscale Simulation of Ternary Stratum Corneum Lipid Mixtures: Effects of Cholesterol Composition. Langmuir 2022, 38, 7496--7511

    2022

  22. [30]

    O.; Moore, T

    Shamaprasad, P.; Frame, C. O.; Moore, T. C.; Yang, A.; Iacovella, C. R.; Bouwstra, J. A.; Bunge, A. L.; MCabe, C. Using Molecular Simulation to Understand the Skin Barrier. Prog. Lipid Res. 2022, 88, 101184

    2022

  23. [31]

    G.; Bunge, A

    Mitragotri, S.; Anissimov, Y. G.; Bunge, A. L.; Frasch, H. F.; Guy, R. H.; Hadgraft, J.; Kasting, G. B.; Lane, M. E.; Roberts, M. S. Mathematical Models of Skin Permeability: An Overview. Int. J. Pharmaceut. 2011, 418, 115--129

    2011

  24. [32]

    F.; Barbero, A

    Frasch, H. F.; Barbero, A. M. Application of Numerical Methods for Diffusion-Based Modeling of Skin Permeation. Adv. Drug Deliver. Rev. 2013, 65, 208--220

    2013

  25. [33]

    M.; Kasting, G

    Wang, J.; Nitsche, J. M.; Kasting, G. B.; Wittum, G.; Nägel, A. Transdermal and Lateral Effective Diffusivities for Drug Transport in Stratum Corneum from a Microscopic Anisotropic Diffusion Model. Eur. J. Pharm. Biopharm. 2023, 188, 271--286

    2023

  26. [34]

    Wang, E.; Klauda, J. B. Models for the Stratum Corneum Lipid Matrix: Effects of Ceramide Concentration, Ceramide Hydroxylation, and Free Fatty Acid Protonation. J. Phys. Chem. B 2018, 122, 11996--12008

    2018

  27. [35]

    G.; Im, W

    Jo, S.; Kim, T.; Iyer, V. G.; Im, W. CHARMM - GUI : A Web-Based Graphical User Interface for CHARMM . J. Comput. Chem. 2008, 29, 1859--1865

    2008

  28. [36]

    B.; Venable, R

    Klauda, J. B.; Venable, R. M.; Freites, J. A.; O’Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A. D. J.; Pastor, R. W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114, 7830--7843

    2010

  29. [37]

    Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I. et al. CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields. J. Comp...

    2010

  30. [38]

    L.; Chandrasekhar, J.; Madura, J

    Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926--935

    1983

  31. [39]

    MacKerell, A. D. Jr .; Bashford, D.; Bellott, M.; Dunbrack, R. L. Jr .; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S. et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586--3616

    1998

  32. [40]

    D.; Friesner, R

    Humphreys, D. D.; Friesner, R. A.; Berne, B. J. A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules. J. Phys. Chem. 1994, 98, 6885--6892

    1994

  33. [41]

    Miyamoto, S.; Kollman, P. A. Settle: An Analytical Version of the SHAKE and RATTLE Algorithm for Rigid Water Models. J. Comput. Chem. 1992, 13, 952--962

    1992

  34. [42]

    Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C. Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes. J. Comput. Phys. 1977, 23, 327--341

    1977

  35. [43]

    Canonical Sampling through Velocity Rescaling

    Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101

    2007

  36. [44]

    L.; Darden, T.; Lee, H.; Pedersen, L

    Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577--8593

    1995

  37. [45]

    C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R

    Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kalé, L.; Schulten, K. Scalable Molecular Dynamics with NAMD . J. Comput. Chem. 2005, 26, 1781--1802

    2005

  38. [46]

    C.; Hardy, D

    Phillips, J. C.; Hardy, D. J.; Maia, J. D. C.; Stone, J. E.; Ribeiro, J. V.; Bernardi, R. C.; Buch, R.; Fiorin, G.; Hénin, J.; Jiang, W. et al. Scalable Molecular Dynamics on CPU and GPU Architectures with NAMD . J. Chem. Phys. 2020, 153, 044130

    2020

  39. [47]

    VMD : Visual Molecular Dynamics

    Humphrey, W.; Dalke, A.; Schulten, K. VMD : Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33--38

    1996

  40. [48]

    J.; Woolf, T

    Michaud-Agrawal, N.; Denning, E. J.; Woolf, T. B.; Beckstein, O. MDAnalysis : A Toolkit for the Analysis of Molecular Dynamics Simulations. J. Comput. Chem. 2011, 32, 2319--2327

    2011

  41. [49]

    G owers; M ax L inke; J onathan B arnoud; T yler J

    R ichard J . G owers; M ax L inke; J onathan B arnoud; T yler J . E . R eddy; M anuel N . M elo; S ean L . S eyler; J an D omański; D avid L . D otson; S ébastien B uchoux; I an M . K enney et al. MDAnalysis : A Python Package for the Rapid Analysis of Molecular Dynamics Simul...

    2016

  42. [50]

    J.; Tobias, D

    Martyna, G. J.; Tobias, D. J.; Klein, M. L. Constant Pressure Molecular Dynamics Algorithms. J. Chem. Phys. 1994, 101, 4177--4189

    1994

  43. [51]

    E.; Zhang, Y.; Pastor, R

    Feller, S. E.; Zhang, Y.; Pastor, R. W.; Brooks, B. R. Constant Pressure Molecular Dynamics Simulation: The Langevin Piston Method. J. Chem. Phys. 1995, 103, 4613--4621

    1995

  44. [52]

    M.; Valleau, J

    Torrie, G. M.; Valleau, J. P. Nonphysical Sampling Distributions in Monte Carlo Free-Energy Estimation: Umbrella Sampling. J. Comput. Phys. 1977, 23, 187--199

    1977

  45. [53]

    M.; Bouzida, D.; Swendsen, R

    Kumar, S.; Rosenberg, J. M.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A. The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules. I . The Method. J. Comput. Chem. 1992, 13, 1011--1021

    1992

  46. [54]

    WHAM : The Weighted Histogram Analysis Method, Version 2.0.11

    Grossfield, A. WHAM : The Weighted Histogram Analysis Method, Version 2.0.11. http://membrane.urmc.rochester.edu/wordpress/?page_id=126 (accessed 2025-01-06)

    2025

  47. [55]

    Convergence and Error Estimation in Free Energy Calculations Using the Weighted Histogram Analysis Method

    Zhu, F.; Hummer, G. Convergence and Error Estimation in Free Energy Calculations Using the Weighted Histogram Analysis Method. J. Comput. Chem. 2011, 33, 453

    2011

  48. [56]

    Standard Error Estimation by an Automated Blocking Method

    Jonsson, M. Standard Error Estimation by an Automated Blocking Method. Phys. Rev. E 2018, 98, 043304

    2018

  49. [57]

    Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method

    Barducci, A.; Bussi, G.; Parrinello, M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100, 020603

    2008

  50. [58]

    L.; Hénin, J

    Fiorin, G.; Klein, M. L.; Hénin, J. Using Collective Variables to Drive Molecular Dynamics Simulations. Mol. Phys. 2013, 111, 3345--3362

    2013

  51. [59]

    Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup–Solute Interactions

    Pokhrel, N.; Maibaum, L. Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup–Solute Interactions. J. Chem. Theory Comput. 2018, 14, 1762--1771

    2018

  52. [60]

    Awoonor-Williams , E.; Rowley, C. N. Molecular Simulation of Nonfacilitated Membrane Permeation. BBA-Biomembranes 2016, 1858, 1672--1687

    2016

  53. [61]

    Permeability across Lipid Membranes

    Shinoda, W. Permeability across Lipid Membranes. BBA-Biomembranes 2016, 1858, 2254--2265

    2016

  54. [62]

    J.; Poolman, B

    Frallicciardi, J.; Melcr, J.; Siginou, P.; Marrink, S. J.; Poolman, B. Membrane Thickness, Lipid Phase and Sterol Type Are Determining Factors in the Permeability of Membranes to Small Solutes. Nat. Commun. 2022, 13, 1605

    2022

  55. [63]

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

    Piasentin, N.; Lian, G.; Cai, Q. Evaluation of Constrained and Restrained Molecular Dynamics Simulation Methods for Predicting Skin Lipid Permeability. ACS Omega 2021, 6, 35363--35374

    2021

  56. [64]

    Sajadi, F.; Rowley, C. N. Simulations of Lipid Bilayers Using the CHARMM36 Force Field with the TIP3P - FB and TIP4P - FB Water Models. PeerJ 2018, 6, e5472

    2018

  57. [65]

    M.; Klauda, J

    Krämer, A.; Ghysels, A.; Wang, E.; Venable, R. M.; Klauda, J. B.; Brooks, B. R.; Pastor, R. W. Membrane Permeability of Small Molecules from Unbiased Molecular Dynamics Simulations. J. Chem. Phys. 2020, 153, 124107

    2020

  58. [66]

    D.; Noro, M

    Das, C.; Olmsted, P. D.; Noro, M. G. Water Permeation through Stratum Corneum Lipid Bilayers from Atomistic Simulations. Soft Matter 2009, 5, 4549--4555

    2009

  59. [67]

    Das, C.; Olmsted, P. D. The Physics of Stratum Corneum Lipid Membranes. Philos. Trans. R. Soc. Math. Phys. Eng. Sci. 2016, 374, 20150126

    2016

  60. [68]

    J.; Allison, J

    Piggot, T. J.; Allison, J. R.; Sessions, R. B.; Essex, J. W. On the Calculation of Acyl Chain Order Parameters from Lipid Simulations. J. Chem. Theory Comput. 2017, 13, 5683--5696

    2017

  61. [69]

    Effect of a Single Cis Double Bond on the Structure of a Phospholipid Bilayer

    Seelig, A.; Seelig, J. Effect of a Single Cis Double Bond on the Structure of a Phospholipid Bilayer. Biochemistry 1977, 16, 45--50

    1977

  62. [70]

    Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study

    Jiang, Z.; Liu, S.; Yuan, S.; Zhang, H.; Yuan, S. Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study. J. Phys. Chem. B 2024, 128, 3833--3843

    2024

  63. [71]

    M.; Katz, Y

    Diamond, J. M.; Katz, Y. Interpretation of Nonelectrolyte Partition Coefficients between Dimyristoyl Lecithin and Water. J. Membrane Biol. 1974, 17, 121--154 mcitethebibliography manuscript.tex0000664000000000000000000013146114744167366012515 0ustar rootroot [journal=jpcbfk,ma...

    1974

  64. [72]

    Convergence and Error Estimation in Free Energy Calculations Using the Weighted Histogram Analysis Method

    Zhu, F.; Hummer, G. Convergence and Error Estimation in Free Energy Calculations Using the Weighted Histogram Analysis Method. J. Comput. Chem. 2012, 33, 453--465

    2012

  65. [73]

    O.; Kyrychenko, A

    Posokhov, Y. O.; Kyrychenko, A. Effect of Acetone Accumulation on Structure and Dynamics of Lipid Membranes Studied by Molecular Dynamics Simulations. Comput. Biol. Chem. 2013, 46, 23--31

    2013

  66. [74]

    Octanol--Water Partition Coefficient Measurement by a Simple ^1 H NMR Method

    Cumming, H.; Rücker, C. Octanol--Water Partition Coefficient Measurement by a Simple ^1 H NMR Method. ACS Omega 2017, 2, 6244--6249

    2017

  67. [75]

    Molecular Order in Cis and Trans Unsaturated Phospholipid Bilayers

    Seelig, J.; Waespe-Sarcevic , N. Molecular Order in Cis and Trans Unsaturated Phospholipid Bilayers. Biochem. 1978, 17, 3310--3315

    1978

  68. [76]

    T.; Vattulainen, I

    Niemelä, P.; Hyvönen, M. T.; Vattulainen, I. Structure and Dynamics of Sphingomyelin Bilayer: Insight Gained through Systematic Comparison to Phosphatidylcholine. Biophys. J. 2004, 87, 2976--2989

    2004

  69. [77]

    )( Q ѾK s

    Merz, K. M. Molecular Dynamics Simulations of Lipid Bilayers. Curr. Opin. Struct. Biol. 1997, 7, 511--517 mcitethebibliography si.tex0000664000000000000000000002776114744167367010753 0ustar rootroot article [sc] mathpazo 1.05 [utf8] inputenc [T1] fontenc [final] microtype [let...

    1997

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.