استندت تحليلات الصيدلة الشبكية والمعلوماتية الحيوية إلى قواعد بيانات متاحة للعموم ولم تشمل مشاركين بشريين أو حيوانات. وقد أجريت التجارب in vitro باستخدام خط خلوي معتمد ولم تتطلب موافقة أخلاقية وفقاً للمتطلبات المؤسسية.
فحص المكونات المرشحة واستعادة الهدف
تَمَّ استرداد المكونات المرشحة وأهدافها المقابلة من منصة تحليل وقاعدة بيانات فارماكولوجيا نظم الطب الصيني التقليدي (TCMSP) باستخدام مصطلح البحث “dongchongxiacao” في 13 يوليو 2025. وقد تم الاحتفاظ بالمركبات التي بلغت التوافر الحيوي الفموي المتوقع لها 30% على الأقل، والتشابه الدوائي (drug-likeness) بنسبة 0.18 على الأقل. واستوفت سبعة مكونات مرشحة كلا معياري الفرز. كما سُجل الوزن الجزيئي، والتوافر الحيوي الفموي، والتشابه الدوائي، وغيرها من واصفات الحركية الدوائية المتاحة من TCMSP لإجراء مقارنة وصفية.
تم تنزيل معلومات المركب-الهدف المرتبطة بالمكونات المستبقاة مباشرة من TCMSP. تمت إزالة إدخالات الأهداف المكررة يدويًا عند بناء مجموعة الأهداف غير المتكررة، بينما تم الاحتفاظ بعلاقات المركب-الهدف الأصلية لبناء الشبكة لاحقًا. ونظرًا لأن TCMSP وفرت أسماء الأهداف في صيغتها الكاملة بدلاً من رموز جينية موحدة، فقد تم البحث عن كل هدف في قاعدة بيانات UniProt (تم الدخول إليها في 14 يوليو 2025). وقد تم اختيار الإدخال البشري المقابل، وتحويل اسم الهدف إلى رمزه الجيني الرسمي. واستُخدمت مجموعة بيانات المركب-الهدف الموحدة لتحليل تقاطع المرض-الهدف وتحليل شبكة المركب-الهدف.
تجميع الجينات المرتبطة بالتهاب كبيبات الكلى وتقاطع الأهداف
تم استرجاع الجينات المرتبطة بالمرض من GeneCards وDisGeNET وPharmGKB وTherapeutic Target Database باستخدام مصطلح البحث "glomerulonephritis" في 14 يوليو 2025. وبالنسبة لـ GeneCards، تم الاحتفاظ بالمدخلات التي تزيد درجات صلتها عن الدرجة المتوسطة البالغة 10. كما تم تضمين المدخلات المؤهلة من قواعد البيانات الأخرى وفقاً لتعليقات المرض-الهدف المنسقة الخاصة بها. وتم توحيد أسماء الجينات وفقاً لرموز الجينات البشرية الرسمية، مع إزالة المدخلات المكررة. وأسفر هذا الإجراء عن تحديد 1,134 جيناً فريداً مرتبطاً بـ glomerulonephritis. ونظراً لأن البحث استخدم المصطلح الأعم "glomerulonephritis"، فقد اعتبرت مجموعة البيانات الناتجة ذات صلة بالعمليات الجزيئية المرتبطة بـ chronic glomerulonephritis بدلاً من أن تكون محددة بنمط نسيجي مرضي واحد.
تم تقاطع الأهداف المرتبطة بالمركبات والمستمدة من TCMSP مع مجموعة الجينات المرتبطة بالمرض باستخدام لغة R. وقد تم تحديد ثلاثة وعشرين هدفاً متداخلة. كما تم إنشاء مخططات فين باستخدام حزمة R المسماة venn، وتم بناء شبكة الأهداف المتداخلة للمركبات وتصورها باستخدام برنامج Cytoscape. وتتلخص العلاقات بين المكونات المرشحة والأهداف المتداخلة في الجدول 1.
تحليلات الإثراء الوظيفي والتفاعل بين البروتينات
تم إجراء تحليلات الإثراء لـ Gene Ontology (GO) وKyoto Encyclopedia of Genes and Genomes (KEGG) في بيئة R باستخدام clusterProfiler وorg.Hs.eg.db. وحُللت مصطلحات Gene Ontology ضمن فئات العملية البيولوجية، والمكون الخلوي، والوظيفة الجزيئية. كما تم تعديل قيم P للاختبارات المتعددة باستخدام طريقة Benjamini-Hochberg، وتم الاحتفاظ بالمصطلحات أو المسارات التي كانت قيم P المعدلة لها أقل من 0.05. وعُرِضت المصطلحات والمسارات المثرى باستخدام مخططات شريطية ومخططات نقطية تم إنشاؤها بشكل مستقل.
تم إدخال الأهداف الـ 23 المتداخلة في قاعدة بيانات STRING، مع اختيار Homo sapiens ككائن حي وتحديد الحد الأدنى لدرجة التفاعل بأكثر من 0.4. ثم استوردت شبكة تفاعلات البروتين-بروتين الناتجة إلى برنامج Cytoscape. وتم تقييم طوبولوجيا الشبكة باستخدام CytoNCA بناءً على درجة الاتصال (degree)، والمركزية البينية (betweenness)، والمركزية التقاربية (closeness)، ومركزية المتجه الذاتي (eigenvector centrality)، ومتوسط الاتصال المحلي (local average connectivity)، ومركزية الشبكة (network centrality). عُرِّفت الأهداف المركزية (Hub targets) بأنها العقد التي كانت قيم درجة الاتصال، والمركزية البينية، والمركزية التقاربية، ومركزية المتجه الذاتي، ومتوسط الاتصال المحلي، ومركزية الشبكة فيها جميعاً أكبر من القيم الوسيطة المقابلة لها في الشبكة. وقد أدى هذا التحليل إلى تحديد تسعة أهداف مركزية ذات أولوية للتفسير اللاحق.
الالتحام الجزيئي والتحقق من صحة إعادة الالتحام
أُجري الإرساء الجزيئي باستخدام بنية سيكلوأكسيجيناز-2 المرتبطة بـ celecoxib ذات المعرف PDB ID 3LN1. تم الاحتفاظ بالسلسلة A ومجموعة الهيم المرتبطة بها كمستقبل، بينما أُزيلت جزيئات الماء البلورية ورابط celecoxib الذي تم بلورته مع البروتين أثناء تحضير المستقبل. استُخلص التشكيل الأصلي لـ celecoxib بشكل منفصل من السلسلة A لغرض التحقق من صحة إعادة الإرساء. أُضيفت ذرات الهيدروجين القطبية وشحنات Gasteiger باستخدام AutoDockTools، وحُفظ المستقبل بتنسيق PDBQT. يمثل PDB ID 3LN1 بنية بدقة 2.40 Å لـ COX-2 الفأري المرتبط بـ celecoxib. تم الحصول على بنيات β-sitosterol و linoleyl acetate من PubChem وحُضرت باستخدام Open Babel. أُضيفت ذرات الهيدروجين وشحنات Gasteiger الجزئية، وحُددت الروابط القابلة للدوران، وحُوّل كل رابط إلى تنسيق PDBQT. أُجري إرساء لكل من celecoxib و β-sitosterol و linoleyl acetate باستخدام AutoDock Vina مع استخدام نفس المستقبل ومعايير البحث.
تم تمركز صندوق الإرساء (docking box) على منطقة ارتباط celecoxib الأصلية باستخدام الإحداثيات التالية: center_x = 30.9886 center_y = -22.2836 center_z = -16.5072 size_x = 26 Å size_y = 26 Å size_z = 26 Å. وقد تم ضبط قيمة الشمولية (exhaustiveness) عند 16، وعدد أنماط المخرجات (output modes) عند 10، ونطاق الطاقة عند 3 kcal/mol. وتم الاحتفاظ بالوضعية (pose) الأعلى تصنيفاً بناءً على أدنى طاقة ارتباط متوقعة لغرض التصوير المرئي. وللتحقق من صحة البروتوكول، تمت إعادة إرساء celecoxib الأصلي تحت نفس الظروف، وحُسب متوسط مربع الانحراف (root mean square deviation) للذرات الثقيلة بين وضعيات البلورة والوضعيات المعاد إرساؤها باستخدام برنامج PyMOL. واعتُبر متوسط مربع الانحراف الذي يقل عن 2 Å إعادة إنتاج مقبولة للوضعية المرجعية. كما تم تحديد البقايا (residues) الواقعة ضمن نطاق 4 Å من كل ربيطة (ligand) تم إرساؤها في برنامج PyMOL. وتم تفسير درجات ووضعيات الإرساء كتقديرات للتوافق البنيوي وليس كدليل مباشر على الارتباط بالهدف أو التثبيط الإنزيمي.
زراعة الخلايا وتحضير المركبات
تمت زراعة خلايا ميزانجيوم الكبيبة الكلوية للجرذ من سلالة HBZY-1 في وسط Dulbecco’s modified Eagle المعدل عالي الجلوكوز والمكمل بـ 10% مصل جنين البقري و1% بنسيلين-ستريبتومايسين عند درجة حرارة 37 °C في ظروف رطبة مع 5% CO₂. استُخدمت الخلايا في الممرات من 2–7 لجميع التجارب، وتم التحقق من خلوها من الميكوبلازما قبل الاستخدام. استُخدم كل من β-Sitosterol وcelecoxib كمركبات محددة كيميائياً. أُذيب β-Sitosterol في الإيثانول المطلق لتحضير محلول مخزون بتركيز 10 mM (حوالي 4.15 mg/mL). وعند الحاجة، أجري تسونيك (sonication) قصير والتدفئة إلى 50–60 °C للحصول على محلول رائق بصرياً. أُذيب celecoxib في ثنائي ميثيل سلفوكسيد لتحضير محلول مخزون بتركيز 10 mM (حوالي 3.81 mg/mL). قُسمت المحاليل المخزونة إلى حصص وتم تخزينها عند −20 °C. وقبل المعالجة، خُففت المحاليل المخزونة مباشرة في وسط زراعة مُدفأ مسبقاً ومُزجت على الفور.
تمت مطابقة تركيزات المذيبات (Vehicle) ضمن كل تجربة. في مقايسة حيوية الخلايا، تم ضبط التركيز النهائي للإيثانول ليكون 0.2% v/v في جميع المجموعات لأن أعلى تركيز تم اختباره من β-sitosterol كان 20 µM. وفي تجربة الالتهاب المستحث بواسطة LPS، تم ضبط التركيزات النهائية للإيثانول وdimethyl sulfoxide لتكون 0.1% و0.01% v/v، على التوالي، في جميع المجموعات. أما في تجربة التليف المستحث بواسطة TGF-β1، فقد تم ضبط التركيز النهائي للإيثانول ليكون 0.1% v/v في جميع المجموعات. كما تمت إضافة أحجام مكافئة من المذيبات المقابلة إلى مجموعات Control وLPS وTGF-β1.
مقايسة حيوية الخلايا
وُزعت الخلايا في أطباق ذات 96 بئراً بكثافة 5 × 103 خلية لكل بئر في 100 µL من الوسط الغذائي الكامل، وتُركت لتلتصق لمدة 24 h. ولتقدير نطاق تركيز β-sitosterol، حُضنت الخلايا مع 0.5 أو 1 أو 5 أو 10 أو 20 µM من β-sitosterol لمدة 24 h. أما بالنسبة لاختبار حيوية الخلايا المرتبط بـ LPS، فقد عولجت الخلايا مسبقاً بـ β-sitosterol بتركيز 1 أو 5 أو 10 µM لمدة 2 h، ثم عُرّضت لـ 1 µg/mL من lipopolysaccharide لمدة 24 h إضافية.
كانت مجموعات المعالجة هي المجموعة الضابطة، وβ-sitosterol بتركيز 0.5 µM، وβ-sitosterol بتركيز 1 µM، وβ-sitosterol بتركيز 5 µM، وβ-sitosterol بتركيز 10 µM، وβ-sitosterol بتركيز 20 µM، وLPS، وLPS مع β-sitosterol بتركيز 1 µM، وLPS مع β-sitosterol بتركيز 5 µM، وLPS مع β-sitosterol بتركيز 10 µM. أضيف كاشف Cell Counting Kit-8 بمقدار 10 µL لكل بئر، وحُضنت الأطباق لمدة 2 h. ثم سُجلت الامتصاصية عند 450 nm. تم تحديد حيوية الخلايا بعد طرح قيمة الفراغ الخالي من الخلايا، وعُبّر عنها بالنسبة إلى المجموعة الضابطة. تضمنت كل حالة ستة آبار تقنية في كل تجربة من ست تجارب بيولوجية مستقلة. وتم حساب متوسط قيم الآبار التقنية من كل تجربة لإنتاج قيمة مكرر بيولوجي واحد للتحليل الإحصائي.
التنشيط الالتهابي الناجم عن LPS
بالنسبة لتجارب التنشيط الالتهابي، تم استزراع خلايا HBZY-1 في أطباق ذات 6 آبار بكثافة 2 × 105 خلية لكل بئر، وعولجت بعد وصولها إلى درجة تلاقي تتراوح تقريباً بين 70%–80%. تم الحفاظ على الخلايا في وسط مغذٍ مضاف إليه 1% من مصل جنين البقر لمدة 12 h قبل المعالجة. تم إعطاء β-Sitosterol أو celecoxib قبل 2 h من التحفيز بواسطة LPS. بعد ذلك، عُرِضت الخلايا لـ 1 µg/mL من LPS لمدة 24 h.
تمثلت المجموعات التجريبية في مجموعة ضابطة، ومجموعة β-sitosterol بتركيز 10 µM، ومجموعة LPS، ومجموعة LPS بالإضافة إلى β-sitosterol بتركيز 1 µM، ومجموعة LPS بالإضافة إلى β-sitosterol بتركيز 5 µM، ومجموعة LPS بالإضافة إلى β-sitosterol بتركيز 10 µM، ومجموعة LPS بالإضافة إلى celecoxib بتركيز 1 µM. تم إدراج celecoxib كضابط موجب وظيفي لتثبيط إنتاج البروستاجلاندين المعتمد على COX-2. وفي نهاية فترة المعالجة، جُمعت السوائل الطافية من المزارع لقياس PGE2. كما جُمعت الخلايا بشكل منفصل لإجراء تحليل PCR الكمي في الوقت الحقيقي لـ Tnf و Il6 و Ptgs2، ولإجراء تحليل western blot لـ COX-2.
التنشيط المحفز للتليف بواسطة TGF-β1
تم استخدام نموذج منفصل محفز بواسطة TGF-β1 لتقييم التنشيط المسبب للتليف. زُرعت خلايا HBZY-1 في أطباق ذات 6 آبار بكثافة 2 × 105 خلية لكل بئر، وحُفظت في وسط مُدعم بـ 1% مصل جنين البقر لمدة 12 h قبل المعالجة. عولجت الخلايا مسبقاً بـ β-sitosterol لمدة 2 h، تلا ذلك التحفيز باستخدام 10 ng/mL من TGF-β1 المؤتلف لمدة 48 h.
شملت مجموعات المعالجة مجموعة ضابطة، وβ-sitosterol (10 µM)، وTGF-β1، وTGF-β1 + β-sitosterol (5 µM)، وTGF-β1 + β-sitosterol (10 µM). وفي نهاية فترة المعالجة، تم حصاد الخلايا لإجراء تحليل PCR الكمي في الوقت الحقيقي لكل من Col1a1 وActa2 ولإجراء تحليل لطخة ويسترن لبروتين alpha-smooth muscle actin (α-SMA).
تفاعل البوليميراز المتسلسل الكمي في الوقت الحقيقي
تم عزل الحمض النووي الريبوزي الكلي (Total RNA) من الخلايا المعالجة باستخدام كاشف لاستخلاص الـ RNA يعتمد على الفينول والغوانيدينيوم وفقاً لبروتوكول الشركة المصنعة. وقد تم تحديد تركيز الـ RNA ونقاوته باستخدام مقياس الطيف الضوئي. واستُخدمت فقط العينات التي تراوحت نسبة A260/A280 فيها بين 1.8–2.0 لعملية النسخ العكسي. وبالنسبة لكل عينة، تم تحويل 1 µg من الـ RNA الكلي عن طريق النسخ العكسي إلى حمض نووي ريبوزي منقوص الأكسجين مكمل (complementary DNA) في حجم تفاعل نهائي قدره 20 µL. ثم خُفف الـ complementary DNA الناتج بمقدار خمسة أضعاف باستخدام ماء خالٍ من النيوكلياز قبل إجراء تفاعل البوليميراز المتسلسل الكمي (quantitative PCR). أُجري تفاعل الـ quantitative PCR في حجم كلي قدره 20 µL يحتوي على 10 µL من 2× SYBR Green master mix، و0.4 µL من كل من البادئات الأمامية والعكسية بتركيز 10 µM، و2 µL من الـ complementary DNA المخفف، و7.2 µL من الماء الخالي من النيوكلياز. وتضمنت ظروف التضخيم خطوة مسخ أولية عند 95 °C لمدة 30 s، تلتها 40 دورة من المسخ عند 95 °C لمدة 5 s، والارتباط-التمديد عند 60 °C لمدة 30 s. وأُجري بعد ذلك تحليل لمنحنى الانصهار من 65 °C إلى 95 °C للتحقق من خصوصية التضخيم.
كانت الجينات التي تم قياسها هي Tnf وIl6 وPtgs2 وCol1a1 وActa2 وGapdh. وقد استُخدم Gapdh كجين مرجعي داخلي. تم تحديد التعبير الجيني النسبي باستخدام طريقة 2−ΔΔCt، مع استخدام مجموعة الضبط كمعاير. تتوفر تسلسلات البادئات في الجدول التكميلي 1. تم تحليل كل عينة في ثلاث مكررات تقنية، وأجريت ست تجارب بيولوجية مستقلة. تم حساب متوسط القيم التي تم الحصول عليها من المكررات التقنية قبل إجراء التحليل الإحصائي.
قياس PGE2 بواسطة المقايسة المناعية المرتبطة بالإنزيم
تم جمع طافي المزرعة بعد 24 h من التحفيز بـ LPS وطردها مركزياً عند 1,000 × g لمدة 10 min عند 4 °C لإزالة حطام الخلايا. تم تقدير تركيزات PGE2 باستخدام طقم مقايسة الممتز المناعي المرتبط بالإنزيم متاح تجارياً وفقاً لتعليمات الشركة المصنعة.
تم تحميل المعايير والعينات على لوحة المقايسة في مكررات ثنائية. وبعد الانتهاء من خطوات ربط الأجسام المضادة وتطوير اللون، تم تسجيل الامتصاصية عند 450 nm باستخدام قارئ الأطباق الدقيقة. تم تحديد تركيزات PGE2 من منحنى قياسي لوجستي رباعي المعلمات وتم التعبير عنها بـ pg/mL. كما تم تحليل السوائل الطافية من ست تجارب خلوية مستقلة. وتم حساب متوسط القياسات المكررة لكل عينة بيولوجية قبل إجراء التحليل الإحصائي.
لطخة ويسترن
تم تحلل الخلايا المعالجة في محلول منظم لمقايسة الترسيب المناعي الإشعاعي (radioimmunoprecipitation assay buffer) مضاف إليه خليط من مثبطات البروتياز. حُضنت المحللات على الثلج لمدة 30 min ثم وُضعت في جهاز الطرد المركزي بسرعة 12,000 × g لمدة 15 min عند درجة حرارة 4 °C. تم تحديد تركيزات البروتين في الرواسب الطافية باستخدام مقايسة البروتين بحمض البيسينكونينيك (bicinchoninic acid protein assay). خُلطت كميات متساوية من البروتين (30 µg لكل مسار) مع محلول التحميل ومُسخنت عند 95 °C لمدة 5 min. فُصلت البروتينات المستخلصة من تجارب LPS على هلامات بولي أكريلاميد وبيريدوديسيل كبريتات الصوديوم (sodium dodecyl sulfate-polyacrylamide gels) بتركيز 10%، بينما فُصلت البروتينات من تجارب TGF-β1 على هلامات بتركيز 12%. بعد ذلك، نُقلت البروتينات المفصولة إلى أغشية بولي فينيليدين فلوريد (polyvinylidene fluoride membranes).
تم حجب الأغشية باستخدام 5% من الحليب الخالي من الدسم في محلول ملحي منظم بالتريس يحتوي على 0.1% Tween 20 لمدة 1 h في درجة حرارة الغرفة، ثم حُضنت طوال الليل عند 4 °C مع أجسام مضادة أولية ضد COX-2 (1:1000)، و α-SMA (1:1000)، و GAPDH (1:5000). تم تحليل COX-2 في عينات من تجربة LPS، بينما تم تحليل α-SMA في عينات من تجربة TGF-β1. وبعد الغسل، حُضنت الأغشية مع أجسام مضادة ثانوية مقترنة ببيروكسيديز الفجل (1:5000) لمدة 1 h في درجة حرارة الغرفة. تم تصوير حزم البروتين باستخدام ركيزة التألق الكيميائي المحسنة والتقاطها بواسطة نظام تصوير التألق الكيميائي. كانت أوقات التعريض 10 s تقريبًا لـ COX-2، و 10 s لـ α-SMA، و 2 s لـ GAPDH.
تم تقدير كثافة الحزم باستخدام برنامج ImageJ. كما تم تطبيع إشارات COX-2 و α-SMA بالنسبة لإشارة GAPDH المقابلة. أُجري تحليل Western blotting باستخدام عينات بروتينية تم الحصول عليها من ست تجارب مستقلة لزراعة الخلايا ومعالجتها.
التحليل الإحصائي
أُجريت التحليلات الإحصائية باستخدام برنامج GraphPad Prism. وتُعرض البيانات في صورة المتوسط ± الانحراف المعياري. نُفذت ست تجارب مستقلة لزراعة الخلايا ومعالجتها لكل تحليل كمي، بحيث تمثل كل تجربة مستقلة مكرراً بيولوجياً واحداً. أما قياسات المكررات التقنية التي تم الحصول عليها من نفس التجربة البيولوجية فقد حُسب متوسطها قبل التحليل الإحصائي ولم تُعتبر ملاحظات مستقلة.
بالنسبة لمقايسة حيوية الخلايا، تضمن كل اختبار بيولوجي ست آبار تقنية لكل حالة. وأُجريت تفاعلات PCR الكمي في ثلاث مكررات تقنية لكل عينة بيولوجية، بينما أُجريت قياسات PGE2 في مكررات ثنائية. واستند التكميم في تحليل Western blot إلى ست عينات بروتينية تم تحضيرها بشكل مستقل. وبناءً على ذلك، يحتوي كل شريط تكميمي على قيم ست مكررات بيولوجية فردية.
تم تحليل الفروق بين المجموعات المتعددة باستخدام تحليل التباين أحادي الاتجاه، متبوعاً باختبار توكي للمقارنات المتعددة. وتم استخدام اختبار ثنائي الجانب P قيمة < اعتبرت القيمة 0.05 دالة إحصائياً. وعُرِضت قيم المكررات البيولوجية الفردية جنباً إلى جنب مع المتوسط والانحراف المعياري.