כל החומרים, כולל מאגרי מידע, תוכנה ומגבים ניסיוניים, ששימשו במחקר זה מופיעים בטבלת החומרים. כל ההליכים הניסיוניים אושרו על ידי ועדת האתיקה לבעלי חיים של בית החולים אגם דושו בסוז'ואו (מספר היתר: 2410008).
סינון רכיבים פעילים ורכישת יעדים של DCQD
הרכיבים הפעילים של צמחי המרפא המרכיבים ריבס, מירביליט, קליפת מגנוליה ותפוז מר לא בוגר במרקם Da-Cheng-Qi (DCQD) נבדקו ממאגר הפרמקולוגיה של מערכות הרפואה הסינית המסורתית (TCMSP) תוך שימוש בקריטריונים הבאים: זמינות ביולוגית דרך הפה (OB) ≥30%) ודמיון תרופה (DL) ≥0.18). המטרות המקבילות של הרכיבים הפעילים שזוהו נאספו לאחר מכן. ספים אלו מאומצים באופן נרחב במחקרי פרמקולוגיה רשתית כדי לסנן תרכובות בעלות תכונות פרמקוקינטיות חיוביות ופוטנציאל גבוה יותר להיות מולקולות דמויות תרופה34,35. כדי לשפר את אמינות בחירת התרכובות ולהפחית הטיה מסינון במאגר יחיד, בוצעה אימות הדדי באמצעות מסד הנתונים BATMAN-TCM. BATMAN-TCM הוא מאגר נתונים אינטגרטיבי המאחסן אינטראקציות ידועות וחזויות בין רכיבי הרפואה הסינית הסינית לחלבוני המטרה, המסייע בחקר מנגנוני הפרמקולוגיה של הרפואה הסינית ובגילוי תרופות. הגרסה המעודכנת 2.0 מציעה מאגר נתונים מורחב משמעותית של TTI ותכונות משופרות, מה שהופך אותה למשאב חשוב להבנת המנגנונים המולקולריים של ה-TCM ולפיתוח טיפולים חדשים למחלות מורכבות. תרכובות ללא מטרות צפויות או עם ציוני מטרות BATMAN-TCM נמוכים מ-20 הוצאו. בנוסף, תרכובות עם ערכי OB נמוכים יחסית אך עם מטרות חזויות ב-BATMAN-TCM ופעילויות פרמקולוגיות מתועדות היטב (כגון אנתראקינונים ופלבונואידים) נשמרו באמצעות אוצרות ספרות. כדי להפחית עוד יותר רכיבים חיוביים שגויים, הורמונים אנדוגניים וחומרים דמויי נוירוטרנסמיטר (כגון פרוגסטרון וסרוטונין) הוצאו ידנית. קבוצת התרכובות הסופית שימשה לחיזוי מטרות ובניית רשת לאחר מכן. תהליך העבודה המפורט מוצג באיור משלים 1. הרשימה הסופית של הרכיבים הפעילים מוצגת בטבלה משלימה 1.
מזהי SMILES (שיטת סימון ליניארית המשתמשת במחרוזות ASCII לתיאור המבנה הכימי של מולקולות. באמצעות שימוש בסמלים וכללים ספציפיים, הוא מקודד אטומים, קשרים כימיים וטופולוגיה מולקולרית לשרשרת רציפה של תווים ללא רווחים, מה שמאפשר אחסון והעברת מידע מולקולרי תמציתי, חד-משמעי וקריא על ידי מחשב) של כל רכיב פעיל נלקח ממאגר PubChem ומטרות פוטנציאליות נחזו באמצעות חיזוי מטרות שווייצריות, כאשר פרמטר האורגניזם הוגדר ל'הומו ספיינס' וחתוך ציון הסתברות של >0. מערכי היעד שהושגו משני מסדי הנתונים אוחדו, ורשומות כפולות הוסרו כדי ליצור את קבוצת היעד המועמדת הסופית ל-DCQD.
חשוב לציין כי Mirabilitum (Mangxiao, Na₂SO₄·10H₂O), תרופה מינרלית ב-DCQD, אינה כלולה ב-TCMSP או במאגרי נתונים דומים למטרות תרכובות. כתוצאה מכך, בשל מגבלת מסד הנתונים הזו, רשת הרכיב-יעד שנבנתה במחקר זה אינה כוללת את התרומה הפוטנציאלית של Mirabilitum. לכן, הרשת מבוססת על המרכיבים הפיטוכימיים של שלושת הרכיבים הצמחיים האחרים: ריאום פאלמטום, מגנוליה אופיצינליס והדר אורנטיום.
ניתוח העשרה של מטרות DCQD
מטרות DCQD יובאו למסד הנתונים המטאסקייפ לניתוח העשרה, כאשר המין מוגבל ל"הומו ספיינס". בוצעו ניתוחי העשרת אונטולוגיית גנים (GO) ואנציקלופדיית קיוטו של גנים וגנומים (KEGG), כאשר הראשון כלל תהליכים ביולוגיים (BP), רכיבים תאיים (CC) ותפקודים מולקולריים (MF). נבחרו ששת התוצאות המובילות בתוצאות העשרה של GO ו-20 הראשונות של תוצאות העשרת KEGG. גרפים של הבר KEGG ו-GO נוצרו באמצעות הכלי המקוון של ביואינפורמטיקה. הנתונים יובאו דרך חבילת DOSE ב-R כדי לבנות את גרף הבר של אונטולוגיית המחלה (DO).
בניית רשת "רכיב-מרפאים-יעד"
הרכיבים הרלוונטיים של DCQD והיעדים המתאימים להם יובאו ל-Cytoscape 3.10.3 כדי לבנות דיאגרמת רשת הממחישה את הקשרים בין "רפואה צמחית-רכיב-מטרה". תוסף CytoNCA שימש לחישוב ערכי דרגות, וחמשת הרכיבים הפעילים המובילים נבדקו לפי דירוג.
שליפת מטרות ARDS ממאגרי מידע ציבוריים
באמצעות "תסמונת מצוקה נשימתית חריפה" כמילת מפתח, מטרות הקשורות ל-ARDS נאספו ממאגרי OMIM ו-GeneCards. בקלפי גנים נבחרו מטרות עם ציון רלוונטיות >5. לאחר מיזוג והסרת כפילויות משני מסדי הנתונים, התקבלה קבוצת יעד מועמדת ל-ARDS.
ניתוח ביואינפורמטיקה
מאגר הנתונים הקשור לספסיס-ARDS GSE32707 (פלטפורמת GPL570, n = 45), הכולל נתוני RNA-seq מתאי דם מונוגרליים היקפיים של 30 חולי ARDS שנגרמו על ידי אלח דם ו-15 ביקורת בריאה, הורד ממאגר ה-GEO של NCBI. הנתונים יובאו באמצעות חבילת GEOquery ב-R, והנרמול בוצע באמצעות חבילת לימה כדי למזער שונות בין דגימות. גנים מבוטאים באופן שונה (DEGs) זוהו באמצעות הספים |log2FC| > 1 ו-p < 0.05. ההדמיה בוצעה באמצעות גרפים של הרי געש ומפות חום. בוצעו ניתוחי העשרת GO ו-KEGG לניתוח המסלול, כאשר התוצאות הוצגו כגרפים בועות וגרפים בר.
חיזוי מטרות פוטנציאליות ל-DCQD בטיפול ב-ARDS
דיאגרמת ון נוצרה באמצעות כלי הביואינפורמטיקה המקוון להשוואת DEGs, מטרות הקשורות ל-ARDS ממאגרי מידע ציבוריים ומטרות תרופות DCQD, ובכך זיהתה מטרות טיפוליות פוטנציאליות ל-DCQD עבור ARDS.
בניית רשת האינטראקציה בין חלבון לחלבון (PPI) וסינון מטרות הליבה
היעדים החופפים יובאו למסד הנתונים של STRING כאשר המין מוגדר ל"הומו ספיינס" וסף ציון אינטראקציה של 0.4. נתוני רשת ה-PPI יוצאו בפורמט TSV והוצגו באמצעות Cytoscape 3.10.3. תוסף CytoHubba יושם לסינון עשרת היעדים המובילים בהתבסס על ציוני MCC, MNC, דרגה, EPC, ביניים, קרבה, רדיאליות ומתח. המפגש בין מטרות אלו הניב 5 מטרות ליבה.
בניית רשת "רפואה צמחית-רכיב-מחלה-מטרה"
בהתבסס על היעדים החופפים, זוהו רכיבי DCQD מתאימים. הרכיבים הרלוונטיים, יחד עם היעדים הקשורים והקשרים המחלות, יובאו ל-Cytoscape 3.10.3 כדי לבנות דיאגרמת רשת הממחישה את הקשרים בין "רפואה צמחית-רכיב-מטרה-מחלה".
ניתוח העשרה של מטרות הליבה ובניית רשת "הרפואה הצמחית-רכיב-מטרה-מחלה-מסלול"
כדי לחקור את הפונקציות הביולוגיות הפוטנציאליות ואת מסלולי האיתות המרכזיים של DCQD בטיפול ב-ARD, היעדים המרכזיים עברו ניתוחי העשרת GO ו-KEGG ב-Metascape (מין: "Homo sapiens"). נבחרו עשרת הראשונים של GO ו-13 הראשונים בתוצאות העשרה של KEGG. ההמחשה בוצעה באמצעות כלי ביואינפורמטיקה ליצירת גרפים בועות KEGG, גרפי GO וגרפי בועה של סאנקי. בנוסף, בהתבסס על תוצאות ההעשרה הפונקציונלית, נבנתה דיאגרמת רשת ב-Cytoscape 3.10.3 כדי להמחיש את הקשרים בין "רפואה צמחית-רכיב-מטרה-מחלה-מסלול".
עגינה מולקולרית והדמיה
המבנים המולקולריים המשניים של התרכובות נאספו ממאגר PubChem, בעוד שהמבנים התלת-ממדיים של מטרות הליבה התקבלו ממסד הנתונים PDB (http://www.rcsb.org/). מולקולות מים ושאריות חומצות אמינו הוסרו באמצעות PyMOL. חלבונים ותרכובות מטרה יובאו ל-AutoDock לצורך עגינה מולקולרית לחישוב אנרגיות קישור. ההמחשה בוצעה באמצעות PyMOL. מפת חום אנרגיית קישור נוצרה באמצעות חבילת pheatmap ב-R, עם מטרות ליבה בציר ה-x, תרכובות בציר ה-y, וגרדיאנט צבע המייצג אנרגיות קישור. דיאגרמות אינטראקציה דו-ממדיות של תוצאות עגינה נוצרו באמצעות Discovery Studio Visualizer 4.5. עגינה מולקולרית בוצעה כדי להעריך את היתכנות הקשר בין תרכובות למטרות, ולא לחזות יעילות מעכבת בפועל.
בניית מודל העכבר ARDS
מחקר זה אושר על ידי ועדת האתיקה הרפואית של בית החולים אגם סוז'ואו דושו (2410008). הושגו עכברים זכרים בריאים C57BL/6J (n = 30; 20 ± 5 גרם). העכברים חולקו לקבוצת ביקורת (n = 10), קבוצת LPS (n = 10, LPS 10 mg/kg), וקבוצת DCQD (n = 10, קבלת DCQD 0.9 גרם לק"ג בשילוב עם LPS 10 mg/kg). עכברים בקבוצת ה-LPS קיבלו הזרקה תוך-פריטונאלית של LPS. עכברים בקבוצת DCQD קיבלו DCQD 0.9 גרם לק"ג באמצעות הערכה אורלית מיד לאחר ההזרקה התוך-פריטוניאלית של LPS. המדידות נלקחו מיד לפני ואחרי הוצאת הריאות מכל העכברים. בהמשך, רקמות הריאה נשמרו בטמפרטורה של -80 מעלות צלזיוס.
ניתוח כתם מערבי (WB)
חלבונים בודדו מליזאטים של רקמת הריאות של העכברים, וריכוזיהם נמדדו באמצעות ערכת כימות חלבון BCA. לאחר מכן, החלבונים עברו דנטור על ידי רתיחה ב-95 מעלות צלזיוס למשך 10 דקות. החלבונים המפוזרים הופרדו על ג'ל SDS-PAGE של 10% שהוכן עם ערכת הכנת ג'ל, ואז הועברו לממברנת PVDF. כדי למנוע קישור לא ספציפי, ממברנות ה-PVDF נחסמו בטמפרטורת החדר למשך 5 דקות באמצעות בופר חסימה מהיר ללא חלבון. לאחר שטיפה עם PBST, הממברנות נדגרו במהלך הלילה בטמפרטורה של 4 מעלות צלזיוס עם נוגדנים ראשוניים נגד MMP9, HSP90, MYC ו-ARG1. לאחר שטיפה נוספת עם PBST, הממברנות דגרו במשך שעתיים בטמפרטורת החדר עם IgG(H+L) אנטי-ארנב מנוגד לעזים בצורת HRP. לבסוף, הוחל סובסטרט ערכת זיהוי כימילומינסנציה רגישה במיוחד, וזוהו רצועות חלבון כימימיומינסנציות באמצעות מערכת הדמיה כימילומינסנציונית BIO-RAD. ניתוח גווני אפור בוצע באמצעות ImageJ. הנוגדנים ששימשו בניסוי זה מופיעים בטבלה משלימה 2.
ניתוח מבחן אימונוסורבנט מקושר לאנזימים (ELISA)
ראשית, נאסף נוזל שטיפה ברונכואלוואולרי (BALF) על ידי ביצוע שטיפה שלוש פעמים עם 2 מ"ל מלח. נוזלי השטיפה שהתקבלו נאספו וצונטריפוגה ב-10,000 גרם במשך 10 דקות להסרת תאים. הסופרנטנט נשאף ואוחסן בהקפאה בחנקן נוזלי. לאחר מכן, נקבעו ריכוזי הציטוקינים המטרה בהתאם להוראות היצרן באמצעות ערכת Mouse IL-6 ELISA, ערכת ELISA של עכבר IL-18 וערכת ELISA TNF-alpha של עכבר. לאחר מכן, לוח התגובה הונח בקורא מיקרופלט, וערך הצפיפות האופטית (OD) של כל באר נמדד באורך גל של 450 ננומטר, כאשר באר הבקרה הריקה שימשה לכיוון אפס. לבסוף, נרשמה עקומה סטנדרטית בהתבסס על הריכוזים וערכי ה-OD המתאימים של התקנים. ריכוז גורם היעד בכל מדגם חושב לפי משוואת העקומה הסטנדרטית.
ניתוח סטטיסטי
הנתונים נותחו באמצעות GraphPad Prism 10.1.2. המשמעות של ההבדלים בין הקבוצות הוערכה סטטיסטית באמצעות ניתוח חד-כיווני של שונות או מבחן t של סטודנט. אם ערך P קטן מ-0.05, ההפרש נחשב משמעותי; אם ערך ה-P קטן מ-0.01, הוא נחשב למשמעותי מאוד.