ניתוחי הפארמקולוגיה של הרשת והביואינפורמטיקה התבססו על מאגרי נתונים הזמינים לציבור ולא כללו משתתפים אנושיים או בעלי חיים. ניסויי ה-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. רשומות מתאימות מהמאגרים האחרים נכללו בהתאם להערות המחלה-מטרה (disease-target annotations) המאושרות שלהם. שמות הגנים הותקנו לסמלי גנים רשמיים של בני אדם, ורשומות כפולות הוסרו. הליך זה הניב 1,134 גנים ייחודיים הקשורים לגלומרולוניפריטיס. מכיוון שהחיפוש השתמש במונח הרחב “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), קישוריות ממוצעת מקומית ומרכזיות הרשת. מטרות צומת (Hub targets) הוגדרו כצמתים שערכיהם של הדרגה, מרכזיות הבינתיות, מרכזיות הקרבה, מרכזיות וקטור העצמי, הקישוריות הממוצעת המקומית ומרכזיות הרשת היו כולם גבוהים מערכי החציונים המתאימים של הרשת. ניתוח זה תיעדן תשע מטרות צומת לצורך פרשנות המשך.
עגינה מולקולרית ותיקוף באמצעות עגינה חוזרת (redocking)
ביצוע עגינה מולקולרית (Molecular docking) בוצע באמצעות מבנה של cyclooxygenase-2 קשור ל-celecoxib, מזהה PDB ID 3LN1. שרשרת A וקבוצת ה-heme המשויכת אליה נשמרו כקולטן, בעוד שמולקולות מים קריסטלוגרפיות והליגנד celecoxib שסו-קריסטלו הוסרו במהלך הכנת הקולטן. הקונפורמציה הטבעית של ה-celecoxib חולצה בנפרד משרשרת A לצורך תיקוף של עגינה מחדש (redocking). אטומי מימן פולריים ומטעני Gasteiger נוספו באמצעות AutoDockTools, והקולטן נשמר בפורמט PDBQT. PDB ID 3LN1 הוא מבנה ברזולוציה של 2.40 Å של COX-2 עכברי הקשור ל-celecoxib. המבנים של β-sitosterol ושל linoleyl acetate ה obtained מ-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) המדורגת בראש הרשימה, על סמך אנרגיית הקישור הנמוכה ביותר שנחזתה, נשמרה לצורך ויזואליזציה. לצורך תיקוף הפרוטוקול, בוצע עיגון חוזר (redocking) של celecoxib טבעי תחת אותם תנאים, וחושב באמצעות PyMOL הסטייה הריבועית הממוצעת (RMSD) של האטומים הכבדים בין התנוחה הקריסטלוגרפית לתנוחה שעוגנה מחדש. סטייה ריבועית ממוצעת הפחותה מ-2 Å נחשבה לשחזור מקובל של תנוחת הייחוס. שאריות הממוקמות בטווח של 4 Å מכל ליגנד שעוגן זוהו ב-PyMOL. ציוני העגינה והתנוחות פורשו כהערכות של תאימות מבנית ולא כראיה ישירה לקישור למטרה או לעכבון אנזימטי.
תרבית תאים והכנת תרכובות
תאי מזנגיום גלומרולריים של חולדה מסוג HBZY-1 תורגדו במצע Dulbecco’s modified Eagle בתוספת גלוקוז גבוה, 10% fetal bovine serum ו-1% penicillin-streptomycin בטמפרטורה של 37 °C בתנאים לחים עם 5% CO₂. עבור כל הניסויים נעשה שימוש בתאים במעברים 2–7, ומצבם כנקיים ממיקופלסמה אומת לפני השימוש. β-Sitosterol ו-celecoxib שימשו כתרכובות מוגדרות כימית. β-Sitosterol הומס באתנול מוחלט כדי להכין תמיסת מלאי בריכוז 10 mM (בקירוב 4.15 mg/mL). במידת הצורך, בוצעה סונקציה קצרה וחימום ל-50–60 °C כדי להגיע לתמיסה צלולה באופן ויזואלי. Celecoxib הומס ב-dimethyl sulfoxide להכנת תמיסת מלאי בריכוז 10 mM (בקירוב 3.81 mg/mL). תמיסות המלאי חולקו למנות ואוחסנו בטמפרטורה של −20 °C. לפני הטיפול, תמיסות המלאי דוללו ישירות לתוך מצע תרבית שחומם מראש ועורבבו באופן מיידי.
ריכוזי חומרי הנושא הותאמו בתוך כל ניסוי. במבחן חיוניות התאים, ריכוז האתנול הסופי הותאם ל-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% התמזגות (confluence). התאים הוחזקו במצע המועשר ב-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. תאים נאספו בנפרד לצורך אנליזת quantitative real-time PCR של Tnf, Il6, ו-Ptgs2 ולצורך אנליזת western blot של COX-2.
הפעלה פרו-פיברוטית המושרית על ידי TGF-β1
מודל נפרד המושרה על ידי TGF-β1 שימש להערכת הפעלה פרופיברוטית. תאי HBZY-1 נזרעו בפלטות בעלות 6 באריים בצפיפות של 2 × 105 תאים לבאר ונשמרו במצע המועשר ב-1% fetal bovine serum למשך 12 h לפני הטיפול. התאים טופלו מראש ב-β-sitosterol למשך 2 h, ולאחר מכן עוברוהם בגירוי של 10 ng/mL recombinant TGF-β1 למשך 48 h.
קבוצות הטיפול כללו ביקורת, β-sitosterol (10 µM), TGF-β1, TGF-β1 + β-sitosterol (5 µM), ו-TGF-β1 + β-sitosterol (10 µM). בסיום תקופת הטיפול, התאים נאספו לצורך ניתוח quantitative real-time PCR של Col1a1 ו-Acta2 ולצורך ניתוח western blot של alpha-smooth muscle actin (α-SMA).
PCR כמותי בזמן אמת (Quantitative real-time PCR)
RNA כולל בודד מהתאים המטופלים באמצעות ריאגנט להוצאת RNA מבוסס פנול-גואנידיניום, בהתאם לפרוטוקול היצרן. ריכוז ה-RNA והטוהר שלו נקבעו באמצעות ספקטרופוטומטר. רק דגימות עם יחס A260/A280 של 1.8–2.0 שימשו לשעתוק הפוך. עבור כל דגימה, 1 µg של RNA כולל שעתקה באופן הפוך ל-DNA משלים בנפח תגובה סופי של 20 µL. ה-DNA המשלים שהתקבל דולל פי חמישה במים נקיים מנוקלאז לפני ביצוע qPCR. ה-qPCR בוצע בנפח כולל של 20 µL המכיל 10 µL של 2× SYBR Green master mix, 0.4 µL של פריימרים קדמיים ואחוריים בריכוז 10 µM כל אחד, 2 µL של DNA משלים מדולל, ו-7.2 µL של מים נקיים מנוקלאז. תנאי ההגברה כללו שלב דה-נטורציה ראשוני ב-95 °C למשך 30 s, ולאחריו 40 מחזורים של דה-נטורציה ב-95 °C למשך 5 s והתאמה-הארכה (annealing-extension) ב-60 °C למשך 30 s. לאחר מכן בוצע ניתוח עקומת המסה (melting-curve analysis) מ-65 °C עד 95 °C כדי לאמת את ספציפיות ההגברה.
הגנים שנמדדו היו Tnf, Il6, Ptgs2, Col1a1, Acta2ו- Gapdh. Gapdh שימש כגן הייחוס הפנימי. ביטוי גנים יחסי נקבע באמצעות שיטת ה-2−ΔΔCt שיטה, כאשר קבוצת הביקורת שימשה ככיול. רצפי הפריימרים מופיעים ב- טבלה נלווית 1כל דגימה נותחה בשלוש חזרות טכניות, ובוצעו שישה ניסויים ביולוגיים עצמאיים. הערכים שהתקבלו מהחזרות הטכניות חושבו כממוצע לפני ביצוע הניתוח הסטטיסטי.
מדידת PGE2 באמצעות תגובת ספיחה חיסונית קשורה לאנזים (ELISA)
על־נוזלי תרבית נאספו לאחר 24 h של גירוי ב-LPS וסובבו בcentrifuge במהירות של 1,000 × g למשך 10 min בטמפרטורה של 4 °C להסרת שיירי תאים. ריכוזי PGE2 נקבעו באמצעות ערכה מסחרית של סריקה חיסונית קשורה לאנזים (ELISA) בהתאם להוראות היצרן.
תקנים ודגימות הועמסו על פלטת הבדיקה בשכפולים. לאחר השלמת שלבי קישור הנוגדנים ופיתוח הצבע, נמדדה הבליעה באורך גל של 450 nm באמצעות קורא מיקרו-פלטות. ריכוזי PGE2 נקבעו על סמך עקומת סטנדרט לוגיסטית בעלת ארבעה פרמטרים והובעו כ-pg/mL. נותחו על-נוזלים מששה ניסויי תאים בלתי תלויים. מדידות שכפול עבור כל דגימה ביולוגית חושבו כממוצע לפני ביצוע הניתוח הסטטיסטי.
Western blotting
תאי שהטופלו עברו ליזיס בבופר למסיי radioimmunoprecipitation assay בתוספת תערובת מעכבי פרוטאז. הליזאטים דגרו על קרח למשך 30 min וסנטריפוגו ב-12,000 × g למשך 15 min בטמפרטורה של 4 °C. ריכוזי החלבונים בעיקרי השכבה (supernatants) נקבעו באמצעות בדיקת חלבון bicinchoninic acid. כמויות שוות של חלבון (30 µg לכל נתיב) עורבבו עם בופר טעינה (loading buffer) ועברו דנטורציה ב-95 °C למשך 5 min. חלבונים מניסויי LPS הופרדו על ג'לים של sodium dodecyl sulfate-polyacrylamide בריכוז 10%, בעוד שחלבונים מניסויי TGF-β1 הופרדו על ג'לים בריכוז 12%. החלבונים המופרדים הועברו לאחר מכן לממברנות polyvinylidene fluoride.
את הממברנות חסמו עם 5% חלב דל שומן ב-Tris-buffered saline המכיל 0.1% Tween 20 למשך 1 h בטמפרטורת החדר, ולאחר מכן דגרו אותן למשך לילה ב-4 °C עם נוגדנים ראשוניים נגד COX-2 (1:1000), α-SMA (1:1000) ו-GAPDH (1:5000). את COX-2 ניתוחו בדגימות מהניסוי עם LPS, בעוד שאת α-SMA ניתוחו בדגימות מהניסוי עם TGF-β1. לאחר שטיפה, דגרו את הממברנות עם נוגדנים משניים המשולבים עם horseradish peroxidase (1:5000) למשך 1 h בטמפרטורת החדר. את פסי החלבון הדמיאו באמצעות substrate לכימילומינצנציה משופרת וצולמו במערכת דימות כימילומינצנטית. זמני החשיפה היו בערך 10 s עבור COX-2, 10 s עבור α-SMA ו-2 s עבור GAPDH.
עוצמות הפסים כוונתו באמצעות ImageJ. אותות של COX-2 ו-α-SMA נורמלו לאות GAPDH התואם. Western blotting בוצע באמצעות דגימות חלבון שהתקבלו משישה ניסויים עצמאיים של תרביות תאים וטיפולים.
ניתוח סטטיסטי
ניתוחים סטטיסטיים בוצעו באמצעות GraphPad Prism. הנתונים מוצגים כממוצע ± סטיית תקן. שישה ניסויי תרבית תאים וטיפולים בלתי תלויים בוצעו עבור כל ניתוח כמותי, כאשר כל ניסוי בלתי תלוי ייצג שכפול ביולוגי אחד. מדידות של שכפולים טכניים שהתקבלו מאותו ניסוי ביולוגי מופעו כממוצע לפני הניתוח הסטטיסטי ולא נחשבו לתצפיות בלתי תלויות.
עבור בדיקת חיוניות התאים, כל ניסוי ביולוגי כלל שש בארות טכניות לכל תנאי. תגובות qPCR בוצעו בשלוש חזרות טכניות עבור כל דגימה ביולוגית, ומדידות של PGE2 בוצעו בשתי חזרות. הכימות של ה-Western blot התבססה על שש דגימות חלבון שהוכנו באופן עצמאי. לפיכך, כל עמודה כמותית כוללת שישה ערכים של חזרות ביולוגיות בודדות.
הבדלים בין קבוצות מרובות נותחו באמצעות ניתוח שונות חד-כיווני (one-way analysis of variance), ולאחריו מבחן השוואות מרובות של Tukey. ערך P דו-צדדי < 0.05 נחשב למשמעותי סטטיסטית. ערכים של חזרות ביולוגיות פרטניות הוצגו יחד עם הממוצע וסטיית התקן.