DOI: https://doi.org/10.1021/acsestair.5c00412
تاريخ النشر: 2026-01-28
المؤلف: Rinto Thomas وآخرون
الموضوع الرئيسي: التطورات في توصيل الأدوية عبر الجلد
نظرة عامة
تستكشف هذه القسم من ورقة البحث انتشار المنتجات العضوية المتطايرة وشبه المتطايرة من زيوت جلد الإنسان عبر غشاء الدهون في الطبقة القرنية (SC)، وهو أمر حاسم لفهم كيمياء الهواء الداخلي. باستخدام محاكاة الديناميات الجزيئية (MD)، تقوم الدراسة بتحديد معاملات الانتشار المعتمدة على الموقع ($D(z)$) للماء، الأسيتون، و6-ميثيل-5-هيبتين-2-ون (6-MHO). تشير النتائج إلى أن الانتشار يتأثر بشكل أساسي بالحواجز الطاقية بدلاً من حركة الجزيئات، مع اختلافات كبيرة في النفاذية بناءً على مقدر الانتشار المستخدم. تشير النتائج إلى أن ترتيب النفاذيات يتبع مشهد الطاقة الحرة، مما يبرز أهمية القوى الديناميكية الحرارية في نقل الجلد.
يؤكد المؤلفون على الحاجة إلى محاكاة دقيقة على المستوى الجزيئي لإبلاغ النماذج الحركية لجودة الهواء الداخلي، حيث غالبًا ما تفتقر النماذج التجريبية الحالية إلى الدقة لالتقاط تعقيدات التفاعلات الجزيئية في مصفوفة الدهون في SC. يقدمون حدودًا عليا وسفلى لمعاملات الانتشار والنفاذية، والتي يمكن أن تُحسن النماذج الحالية وتعزز الدقة التنبؤية. تشمل اتجاهات البحث المستقبلية استكشاف هياكل SC متعددة الطبقات، وتقييم تأثير مجالات القوة المختلفة، والتحقق من النتائج باستخدام بيانات تجريبية للمواد العضوية المتطايرة ذات الصلة بالبيئات الداخلية. تؤسس هذه العمل إطارًا لربط آليات النقل الذري بنماذج أكبر لعرض الإنسان وجودة الهواء.
مقدمة
في عام 2024، تجاوزت درجة الحرارة العالمية القريبة من السطح المتوسط قبل الثورة الصناعية بأكثر من 1.5 كلفن، مما يبرز التأثيرات الكبيرة لتغير المناخ على صحة الإنسان وسلوكه، خاصة بين الفئات الضعيفة. استجابةً للحرارة الشديدة، قد يسعى الأفراد بشكل متزايد للجوء إلى الداخل، حيث يتم التوسط في تأثيرات تغير المناخ من خلال جودة الهواء الداخلي. تتناول هذه الدراسة الأدلة التجريبية المحدودة حول كيفية تأثير ارتفاع درجات الحرارة الخارجية على البيئات الداخلية من خلال التحقيق في انتشار المركبات العضوية المتطايرة وشبه المتطايرة (VOCs)، وخاصة الأسيتون و6-MHO، عبر جلد الإنسان. يتم توليد هذه المركبات خلال الأوزونوليز لدهون الجلد، مما يجعلها ملوثات ثانوية داخلية.
تستخدم البحث محاكاة الديناميات الجزيئية الذرية (MD) لاستكشاف طاقة الانتشار عبر أغشية الدهون في الطبقة القرنية، التي تعمل كحاجز للاختراق الكيميائي. تكشف الدراسة أن تقارب ملفات الطاقة الحرة لنقل المواد المتخللة بطيء، ويعزى ذلك إلى الهيكل المنظم لأغشية الجلد مقارنةً بغشاء الدهون التقليدي. من خلال تقدير ملفات الانتشار على طول العمود الطبيعي للغشاء باستخدام طريقتين تكميليتين—وظيفة الارتباط الذاتي للسرعة (VACF) ووظيفة الارتباط الذاتي للموقع (PACF)—تقدم المؤلفون رؤى حول التحديات المتعلقة بتوقع معاملات الانتشار بدقة خلال الانتشار البطيء. تسهم النتائج في فهم أفضل لانتشار الجزيئات عبر أغشية الدهون في الجلد وتهدف إلى تحسين النماذج التنبؤية لكيمياء الهواء الداخلي في سياق تغير المناخ.
طرق
في هذه الدراسة، نستخدم محاكاة الديناميات الجزيئية الذرية (MD) للتحقيق في حركة الجزيئات الفردية أثناء عبورها لحاجز الدهون في الجلد النموذجي. تتيح لنا هذه الطريقة، المستندة إلى الميكانيكا الكلاسيكية، تحديد خصائص النقل على مستوى ميكروسكوبي، وهو أمر أساسي لإبلاغ النماذج الأكبر المتعلقة بنقل الجلد وكيمياء الهواء الداخلي. تم توضيح بروتوكول المحاكاة والأنظمة النموذجية المستخدمة في أبحاثنا السابقة.
تشمل الجوانب الرئيسية لمنهجيتنا حساب معاملات الانتشار المعتمدة على الموقع، والتي توفر رؤى حول كيفية اختلاف حركة الجزيئات عبر مناطق مختلفة من الحاجز. بالإضافة إلى ذلك، نقوم بإجراء تحليل للمُروّج لتوضيح ديناميات النقل الجزيئي بشكل أكبر. تسهم هذه النتائج في فهم أعمق للآليات التي تحكم نفاذية الحواجز الشبيهة بالجلد.
نتائج
في هذا القسم، يناقش المؤلفون التحديات المرتبطة باستقراء معاملات الانتشار المعتمدة على الموقع من وظائف الارتباط الذاتي للسرعة (VACFs). بشكل محدد، يبرزون الحاجة إلى حساب معامل انتشار يعتمد على تردد لابلاس، \( D(s) \)، لكل نافذة محاكاة، يليها استقراء إلى تردد صفر. تُعقد هذه العملية بسبب اعتماد \( s^{-1} \) لـ \( D(s) \)، مما يؤدي إلى حد غير محدد يتطلب طريقة استقراء تجريبية. استخدم المؤلفون تقنية الاستقراء الخطي الآلي من Rowley وآخرون، والتي، على الرغم من فعاليتها في منطقة الكتلة المائية، أثبتت أنها أقل موثوقية في مركز الغشاء بسبب الانحناء الملحوظ لـ \( D(s) \). ونتيجة لذلك، اختلفت تقديرات معامل الانتشار بشكل كبير، حيث تراوحت من حوالي \( 0.1 \, \text{Å}^2/\text{ps} \) إلى \( 0.75 \, \text{Å}^2/\text{ps} \).
بالإضافة إلى ذلك، يؤكد المؤلفون على الطبيعة المعتمدة على الحجم للانتشار، خاصة داخل المنطقة المائية المحدودة لنظام غشائهم. أجروا تحليلات لمعاملات الانتشار في خلية محاكاة مستطيلة، مقارنةً بالانتشار أحادي البعد على طول المحور z—الذي يعكس إعداد الغشاء—مع الانتشار ثلاثي الأبعاد، الذي هو أكثر شيوعًا ويوفر عينة إحصائية أفضل. تشير النتائج إلى أن اختيار بعد الانتشار يؤثر بشكل كبير على الانزياحات المربعة المتوسطة (MSD) الملحوظة لمواد مختلفة، بما في ذلك الماء، الأسيتون، و6-MHO، كما هو موضح في مخططات MSD المعتمدة على الزمن.
مناقشة
في هذا القسم، يستكشف المؤلفون انتشار المواد المذابة (الأسيتون، 6-MHO، والماء) عبر غشاء الدهون في الطبقة القرنية (SC) النموذجي، مع التركيز على مرحلة التكرار القصير (SPP) لمصفوفة الدهون. يتميز SC بهيكل معقد “من الطوب والملاط”، يتكون بشكل أساسي من السيراميدات، والأحماض الدهنية الحرة، والكوليسترول، المنظمة في مراحل ميكروفازية لاملامسية متميزة. تستخدم الدراسة محاكاة الديناميات الجزيئية (MD) وتقنيات أخذ العينات المظلية لتقدير معاملات الانتشار المعتمدة على الموقع، \(D(z)\)، التي تعكس الحركة المحلية للمواد المذابة داخل الغشاء. يتم مقارنة طريقتين لحساب \(D(z)\)—وظيفة الارتباط الذاتي للسرعة (VACF) ووظيفة الارتباط الذاتي للموقع (PACF)—مما يكشف عن اختلافات كبيرة في تقديرات معامل الانتشار، خاصة في مركز الغشاء، حيث يتنبأ VACF بمعاملات انتشار أعلى من PACF.
يحلل المؤلفون أيضًا آثار هذه الملفات الخاصة بالانتشار على النفاذية العامة للغشاء، \(P\)، باستخدام نموذج الذوبان-الانتشار غير المتجانس. يجدون أن اختيار طريقة الانتشار يؤدي إلى اختلافات في تقديرات النفاذية، مع اختلافات تصل إلى ترتيب من حيث الحجم عبر المواد المذابة. على الرغم من هذه التباينات، يستنتج المؤلفون أن التغيرات الإحصائية في ملفات الطاقة الحرة لها تأثير ضئيل على توقعات النموذج الحركي. يؤكدون أن معامل الانتشار، والطاقة الحرة، والانتشار هي جوانب مترابطة ولكن متميزة من النقل، محذرين من الخلط بين التأثيرات الديناميكية الحرارية والحركية. بشكل عام، توفر النتائج فهمًا دقيقًا لنقل المواد المذابة عبر أغشية الدهون، مما يبرز أهمية الخيارات المنهجية في تقدير معاملات الانتشار والنفاذية.
DOI: https://doi.org/10.1021/acsestair.5c00412
Publication Date: 2026-01-28
Author(s): Rinto Thomas et al.
Primary Topic: Advancements in Transdermal Drug Delivery
Overview
This section of the research paper investigates the permeation of volatile and semivolatile organic oxidation products from human skin oils through the stratum corneum (SC) lipid membrane, which is crucial for understanding indoor air chemistry. Using molecular dynamics (MD) simulations, the study quantifies position-dependent diffusivities ($D(z)$) for water, acetone, and 6-methyl-5-hepten-2-one (6-MHO). The results indicate that permeation is primarily influenced by energetic barriers rather than molecular mobility, with permeabilities varying significantly based on the diffusivity estimator used. The findings suggest that the order of permeabilities follows the free-energy landscape, highlighting the importance of thermodynamic driving forces in skin transport.
The authors emphasize the need for accurate molecular-scale simulations to inform kinetic models of indoor air quality, as existing empirical models often lack the resolution to capture the complexities of molecular interactions in the SC lipid matrix. They provide upper and lower bounds for diffusivities and permeabilities, which can refine current models and enhance predictive accuracy. Future research directions include exploring multilamellar SC architectures, assessing the impact of different force fields, and validating findings with experimental data for volatile organic compounds relevant to indoor environments. This work establishes a framework for linking atomistic transport mechanisms to larger-scale models of human exposure and air quality.
Introduction
In 2024, the global near-surface temperature surpassed the pre-industrial average by over 1.5 K, highlighting the significant effects of climate change on human health and behavior, particularly among vulnerable populations. As a response to extreme heat, individuals may increasingly seek refuge indoors, where the impacts of climate change are mediated through indoor air quality. This study addresses the limited empirical evidence on how rising outdoor temperatures affect indoor environments by investigating the permeation of volatile and semivolatile organic compounds (VOCs), specifically acetone and 6-MHO, through human skin. These compounds are generated during the ozonolysis of skin lipids, serving as secondary indoor pollutants.
The research employs atomistic molecular dynamics (MD) simulations to explore the permeation energetics across stratum corneum lipid membranes, which act as a barrier to chemical penetration. The study reveals that the convergence of free-energy profiles for permeant translocation is slow, attributed to the ordered structure of skin membranes compared to conventional lipid bilayers. By estimating diffusivity profiles along the membrane normal using two complementary methods—velocity autocorrelation function (VACF) and position autocorrelation function (PACF)—the authors provide insights into the challenges of accurately predicting diffusivities during slow permeation. The findings contribute to a better understanding of molecular permeation across skin lipid membranes and aim to enhance predictive models of indoor air chemistry in the context of climate change.
Methods
In this study, we employ atomistic molecular dynamics (MD) simulations to investigate the motion of individual molecules as they traverse a model skin-lipid barrier. This method, grounded in classical mechanics, enables us to quantify the transport properties at a microscopic level, which are essential for informing larger-scale models related to skin transport and indoor air chemistry. The simulation protocol and model systems utilized have been elaborated upon in our previous research.
Key aspects of our methodology include the computation of position-dependent diffusivities, which provide insights into how molecular movement varies across different regions of the barrier. Additionally, we conduct a propagator analysis to further elucidate the dynamics of molecular transport. These findings contribute to a deeper understanding of the mechanisms governing the permeability of skin-like barriers.
Results
In this section, the authors discuss the challenges associated with extrapolating position-dependent diffusivities from velocity autocorrelation functions (VACFs). Specifically, they highlight the need to compute a Laplace-frequency-dependent diffusivity, \( D(s) \), for each simulation window, followed by an extrapolation to zero frequency. This process is complicated by the \( s^{-1} \) dependence of \( D(s) \), leading to an ill-defined limit that requires an empirical extrapolation method. The authors employed the automated linear extrapolation technique from Rowley et al., which, while effective in the aqueous bulk region, proved less reliable at the membrane center due to the pronounced curvature of \( D(s) \). Consequently, diffusivity estimates varied significantly, ranging from approximately \( 0.1 \, \text{Å}^2/\text{ps} \) to \( 0.75 \, \text{Å}^2/\text{ps} \).
Additionally, the authors emphasize the size-dependent nature of diffusion, particularly within the finite aqueous region of their membrane system. They conducted analyses of bulk diffusivities in a rectangular simulation cell, comparing one-dimensional diffusion along the z-axis—reflecting the membrane setup—with three-dimensional diffusion, which is more common and offers better statistical sampling. The results indicate that the choice of diffusion dimension significantly impacts the observed mean squared displacements (MSD) for different substances, including water, acetone, and 6-MHO, as illustrated in their time-dependent MSD plots.
Discussion
In this section, the authors investigate the permeation of solutes (acetone, 6-MHO, and water) through a model stratum corneum (SC) lipid bilayer, focusing on the short periodicity phase (SPP) of the lipid matrix. The SC is characterized by a complex “brick-and-mortar” architecture, primarily composed of ceramides, free fatty acids, and cholesterol, organized into distinct lamellar microphases. The study employs molecular dynamics (MD) simulations and umbrella sampling techniques to estimate position-dependent diffusivities, \(D(z)\), which reflect the local mobility of solutes within the membrane. Two methods for calculating \(D(z)\)—the velocity autocorrelation function (VACF) and the position autocorrelation function (PACF)—are compared, revealing significant differences in diffusivity estimates, particularly at the membrane center, where the VACF predicts higher diffusivities than the PACF.
The authors further analyze the implications of these diffusivity profiles on the overall membrane permeability, \(P\), using the inhomogeneous solubility-diffusion model. They find that the choice of diffusivity method leads to variations in permeability estimates, with differences reaching an order of magnitude across solutes. Despite these discrepancies, the authors conclude that the statistical variations in free energy profiles have minimal impact on kinetic model predictions. They emphasize that diffusivity, free energy, and permeation are interconnected yet distinct aspects of transport, cautioning against conflating thermodynamic and kinetic influences. Overall, the findings provide a nuanced understanding of solute transport across lipid membranes, highlighting the importance of methodological choices in estimating diffusivities and permeabilities.