النمذجة متعددة المقاييس لإصابة الكبد الناتجة عن الأدوية من العضو إلى الفصيص
Multiscale modeling of drug-induced liver injury from organ to lobule

المجلة: npj Digital Medicine، المجلد: 8، العدد: 1
DOI: https://doi.org/10.1038/s41746-025-01736-6
PMID: https://pubmed.ncbi.nlm.nih.gov/40550856
تاريخ النشر: 2025-06-23
المؤلف: Alon Malka-Markovitz وآخرون
الموضوع الرئيسي: التسمم الكبدي الناتج عن المخدرات والحماية

نظرة عامة

تحدد هذه الفقرة التحديات التي تطرحها إصابة الكبد الناتجة عن الأدوية في كل من تطوير الأدوية والإعدادات السريرية. تقدم الدراسة بحثًا سابقًا من خلال تقديم إطار عمل ديناميكي سوائل حاسوبي متعدد المقاييس يستخدم بيانات تشريحية محددة للمرضى للتنبؤ بإصابة الكبد الناتجة عن الأسيتامينوفين. يدمج هذا النموذج بفعالية المقاييس الوعائية واللوبيولية والخلوية لمحاكاة تدفق الدم الديناميكي، ونقل الأدوية، وآليات الإصابة، مما يعكس التباين المكاني الملحوظ في السيناريوهات السريرية.

تشير النتائج إلى أن النموذج يعكس بدقة ديناميات تدفق الدم ويتنبأ بتلف الخلايا الكبدية، مما يوضح إمكانيته كإطار عمل قابل للتوسع للتحقيق في التباين المكاني في مختلف أمراض الكبد. تمهد هذه الأبحاث الطريق لمنهجية شاملة لمحاكاة الكبد الافتراضي، والتي يمكن أن تكون أداة قيمة في تعزيز السلامة أثناء تطوير العلاجات والعلاج السريري، مما يقلل من الاعتماد على نماذج الحيوانات في اختبار الأدوية قبل السريرية.

طرق

تحدد فقرة “الطرق” المواد والمنهجيات المستخدمة في البحث. توضح تصميم التجربة، بما في ذلك اختيار المواد، وتحضير العينات، والتقنيات المحددة المستخدمة لجمع البيانات وتحليلها. تؤكد الفقرة على صرامة المنهجية لضمان إمكانية تكرار النتائج وموثوقيتها.

تُبرز النتائج الرئيسية من تطبيق هذه الطرق، مما يظهر فعاليتها في معالجة الأسئلة البحثية المطروحة. قد تناقش الفقرة أيضًا أي تحليلات إحصائية تم إجراؤها، بما في ذلك النماذج المستخدمة لتفسير البيانات، مما يوفر فهمًا شاملاً للإطار التجريبي وآثاره على استنتاجات الدراسة.

نتائج

كشفت نتائج محاكاة الوريد البابي، استنادًا إلى أشكال هندسية محددة للمرضى، عن ضغوط دخول قدرها $1150 \pm 50$ با، ومعدلات تدفق تتراوح بين $9.48$ و $11.43$ سم³/ث، مما يتماشى جيدًا مع الدراسات السابقة. أشارت أنماط السرعة إلى تدفق غير منتظم بسبب التركيب المعقد للوريد، مع سرعات أعلى في الوريد المركزي وسرعات أقل بالقرب من الجدران، متوافقة مع قياسات التصوير بالرنين المغناطيسي 4D في الجسم الحي. تراوحت معدلات تدفق المخرج من $0.6$ إلى $2.1$ سم³/ث، مع سرعات تتراوح بين $0.035 \pm 0.007$ م/ث و $0.10 \pm 0.02$ م/ث. أظهرت المحاكاة نقلًا فعالًا لـ APAP عبر جميع المخرجات، مع كسر كتلة قدره $0.3 \times 10^{-3}$، وأكدت تحليلات الحساسية أن التغيرات في تركيز الدواء كان لها تأثير ضئيل على ديناميات التدفق.

أثبتت محاكاة إضافية باستخدام نموذج الكابيلاري البديل ثنائي الأبعاد (CSM) قابلية التوسع ودقة النموذج في التقاط ديناميات الأوعية الدقيقة. أظهرت محاكاة CSM الكبيرة انخفاضًا في سرعات الدخول بنسبة $16-9\%$ عبر المخرجات بسبب التفرع، بينما أظهرت CSM الصغيرة انخفاضًا بنسبة $12-8\%$. أكدت كلا النموذجين نقل APAP عبر جميع المخرجات. أشارت محاكاة تدفق الدواء اللوبيولي إلى سرعات دخول تتراوح من $140$ إلى $1260$ ميكرومتر/ث، مع قيم نفاذية تتماشى مع الأدبيات. أظهرت ملفات السرعة علاقات خطية مع سرعات الدخول، مما يؤكد البيانات التجريبية. بالإضافة إلى ذلك، أشار نموذج تلف الخلايا إلى أن جرعة زائدة قدرها $27.1$ غرام أدت إلى تلف يتماشى مع النتائج السابقة. تم استكشاف قدرات النموذج التنبؤية من خلال محاكاة إصابة الكبد الناتجة عن الأسيتامينوفين، مما يكشف أن زيادة قدرها $2.5\%$ في ضغط مخرج الأوعية يمكن أن تؤثر بشكل كبير على نقل الدواء ومقاومة اللوبيول، مما يبرز التفاعل المعقد بين مقاومة الأوعية وتوزيع الدواء في سياق إصابة الكبد.

مناقشة

تقدم الدراسة نموذجًا شاملًا متعدد المقاييس يدمج محاكاة ديناميكا السوائل الحاسوبية (CFD) عبر تسلسل الأوعية الدموية في الكبد، مما يسهل التنبؤات المحددة للمرضى حول توزيع الدواء، والتمثيل الغذائي، والسمية الكبدية. باستخدام مخطط نموذج الكابيلاري البديل المزدوج (CSM)، يربط الإطار بفعالية بين الديناميات على المستويين الكلي والجزئي، مما يوفر رؤى جديدة حول نقل الدواء وتلف الأنسجة أثناء امتصاص الأسيتامينوفين (APAP). تؤكد مصادقة النموذج ضد القياسات في الجسم الحي والنماذج السابقة في الحاسوب قدرته على تكرار الظروف الديناميكية المعروفة، على الرغم من أن القيود في دقة التصوير تعيق المصادقة في المقاطع الوعائية الأصغر. من الجدير بالذكر أن المحاكاة تكشف عن تباين مكاني كبير في تدفق الدم والضغط، خاصة في الفص الأيمن من الكبد، مما يبرز ضرورة أخذ التباين الإقليمي في الاعتبار عند التنبؤ بديناميات الدواء.

ينجح نموذج CSM في توليد بيانات تدفق مفصلة لحوالي 100,000 مخرج لوبيولي، مما يربط ديناميات الوريد البابي على المستوى الكلي بتدفقات الجيوب على المستوى الجزئي. تسمح هذه القابلية للتوسع بتقييم التلف اللوبيولي المحلي على صحة الكبد العامة، مما يلتقط حلقات التغذية الراجعة الحرجة التي تؤثر على إعادة توزيع تدفق الدم استجابةً للإصابة. يوضح النموذج العلاقة بين سرعة الدخول والتدفق داخل الأنسجة اللوبيولية، متماشيًا مع البيانات التجريبية، ويظهر أن حلقة واحدة من إصابة الكبد الحادة لا تغير بشكل كبير ديناميات توزيع الدواء. بينما يلتقط النموذج الحالي بشكل فعال تباين تدفق الأوعية، فإن التحسينات المستقبلية ضرورية لدمج شروط الحدود الديناميكية وإعادة البناء ثلاثي الأبعاد، مما يعزز الدقة التنبؤية ويوسع قابليته للتطبيق على مختلف أمراض الكبد واستراتيجيات العلاج. بشكل عام، يمثل هذا النهج في النمذجة متعددة المقاييس تقدمًا كبيرًا في فهم وظيفة الكبد وإصابة الكبد الناتجة عن الأدوية، مما يمهد الطريق للطب الشخصي.

Journal: npj Digital Medicine, Volume: 8, Issue: 1
DOI: https://doi.org/10.1038/s41746-025-01736-6
PMID: https://pubmed.ncbi.nlm.nih.gov/40550856
Publication Date: 2025-06-23
Author(s): Alon Malka-Markovitz et al.
Primary Topic: Drug-Induced Hepatotoxicity and Protection

Overview

The section outlines the challenges posed by drug-induced liver injury in both drug development and clinical settings. The study advances previous research by introducing a Multiscale Computational Fluid Dynamics framework that utilizes patient-specific anatomical data to predict acetaminophen-induced liver injury. This model effectively integrates vascular, lobular, and cellular scales to simulate dynamic blood flow, drug transport, and injury mechanisms, capturing the spatial heterogeneity observed in clinical scenarios.

The findings indicate that the model accurately reflects blood flow dynamics and predicts hepatocellular damage, demonstrating its potential as a scalable framework for investigating spatial heterogeneity in various hepatic pathologies. This research lays the groundwork for a comprehensive virtual liver simulation methodology, which could serve as a valuable tool in enhancing safety during therapeutic development and clinical treatment, thereby reducing the dependence on animal models in preclinical drug testing.

Methods

The “Methods” section outlines the materials and methodologies employed in the research. It details the experimental design, including the selection of materials, sample preparation, and the specific techniques used for data collection and analysis. The section emphasizes the rigor of the methodology to ensure reproducibility and reliability of the results.

Key findings from the application of these methods are highlighted, demonstrating their effectiveness in addressing the research questions posed. The section may also discuss any statistical analyses performed, including the models used to interpret the data, thereby providing a comprehensive understanding of the experimental framework and its implications for the study’s conclusions.

Results

The results of the portal vein simulations, based on patient-specific geometries, revealed inlet pressures of $1150 \pm 50$ Pa and flow rates between $9.48$ and $11.43$ cm³/s, aligning well with previous studies. The velocity patterns indicated non-uniform flow due to the complex vein structure, with higher velocities in the central vein and lower velocities near the walls, consistent with in vivo 4D MRI measurements. Outlet flow rates ranged from $0.6$ to $2.1$ cm³/s, with velocities between $0.035 \pm 0.007$ m/s and $0.10 \pm 0.02$ m/s. The simulations demonstrated effective APAP transport across all outlets, with a mass fraction of $0.3 \times 10^{-3}$, and sensitivity analyses confirmed that variations in drug concentration had minimal impact on flow dynamics.

Further simulations using the 2D Capillary Surrogate Model (CSM) validated the scalability and accuracy of the model in capturing microvascular dynamics. The Large CSM simulations showed inlet velocities decreasing by $16-9\%$ across outlets due to bifurcation, while the Small CSM exhibited a $12-8\%$ decrease. Both models confirmed APAP transport across all outlets. Lobular drug flow simulations indicated inlet velocities from $140$ to $1260$ μm/s, with permeability values consistent with literature. The velocity profiles demonstrated linear relationships with inlet velocities, corroborating experimental data. Additionally, the cellular damage model indicated that a $27.1$ g overdose resulted in damage consistent with previous findings. The model’s predictive capabilities were further explored through simulations of acetaminophen-induced liver injury, revealing that a $2.5\%$ increase in vascular outlet pressure could significantly affect drug transport and lobular resistance, highlighting the complex interplay between vascular resistance and drug distribution in the context of liver injury.

Discussion

The study introduces a comprehensive multiscale model that integrates computational fluid dynamics (CFD) simulations across the liver’s vascular hierarchy, facilitating patient-specific predictions of drug distribution, metabolism, and hepatotoxicity. Utilizing a dual Capillary Surrogate Model (CSM) scheme, the framework effectively bridges macro- and microscale dynamics, yielding novel insights into drug transport and tissue damage during acetaminophen (APAP) absorption. The model’s validation against in vivo measurements and previous in silico models demonstrates its capability to replicate known hemodynamic conditions, although limitations in imaging resolution hinder validation in smaller vascular segments. Notably, the simulations reveal significant spatial heterogeneity in blood flow and pressure, particularly in the right liver lobe, underscoring the necessity of accounting for regional variability in predicting drug dynamics.

The CSM successfully generates detailed flow data for approximately 100,000 lobule outlets, linking macro-level portal vein dynamics to micro-level sinusoidal flows. This scalability allows for the assessment of localized lobular damage on overall liver health, capturing critical feedback loops that influence blood flow redistribution in response to injury. The model elucidates the relationship between inlet velocity and flow within lobular tissue, aligning with experimental data, and demonstrates that a single episode of acute liver injury does not significantly alter drug distribution dynamics. While the current model effectively captures vascular inflow heterogeneity, future refinements are necessary to incorporate dynamic boundary conditions and 3D reconstructions, enhancing predictive accuracy and expanding its applicability to various liver diseases and therapeutic strategies. Overall, this multiscale modeling approach represents a significant advancement in understanding liver function and drug-induced liver injury, paving the way for personalized medicine.