ניתוח פרמקולוגיה רשתית
"אוטיזם" שימש כמילה המפתח לחיפוש במאגר GeneCards (https://www.genecards.org/) אחר מטרות הקשורות לאוטיזם 12,13,14. מטרות שנאספו ממאגרי הנתונים אוחדו לאחר שהוסרו רשומות כפולות. מטרות ידועות של רכיבים פעילים שלא תפסו על ידי תחזיות מסד הנתונים הושלמו בהתבסס על דוחות ספרות15. מטרות המחלה והיעדים הפוטנציאליים של רכיבי התרופות הוסדרו באופן אחיד לסמלי גנים באמצעות מאגר החלבונים UniProt (https://www.uniprot.org/), ושתי קבוצות היעדים הללו מופו כדי לזהות את היעדים הטיפוליים הפוטנציאליים של KSZZD לאוטיזם16.
בהתבסס על בסיס הנתונים והפלטפורמת ניתוח של מערכות פרמקולוגיה מסורתית לרפואה סינית (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php), התרכובות הביואקטיביות הפוטנציאליות ב-KSZZD והיעדים המתאימים שלהן נבדקו בהתאם לקריטריונים של זמינות ביולוגית דרך הפה (OB ≥ 0.30) ומדד דמיון תרופות (DL ≥ 0.18)17. מזהי SMILES של התרכובות נאספו ממסד הנתונים PubChem (https://pubchem.ncbi.nlm.nih.gov/), ו-Swiss Target Prediction (http://swisstargetprediction.ch/) שימש עוד כדי לחקור את היעדים הפוטנציאליים שאינם כלולים בפלטפורמת TCMSP18. בהמשך, מאגר הנתונים HERB (http://herb.ac.cn) שימש לסינון משלים בהתאם לכלל החמשת של ליפינסקי, עם הקריטריונים הבאים: משקל מולקולרי (MW ≤ 500 Da), מקדם חלוקת אוקטנול-מים (AlogP ≤ 5), מספר תורמי קשרי מימן (Hdon ≤ 5), מספר קבלי קשרי מימן (Hacc ≤ 10), ומספר הקשרים הניתנים לסיבוב (RBN ≤ 10). לאחר מכן, בוצעה סקר מטרות באמצעות מאגר הנתונים BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/) עם הקריטריונים הבאים: סף ציון (≥0.84), ציון יכולת התרופות (≥0.10), וערך P (≤0.05). לבסוף, מסד הנתונים UniProt (https://www.uniprot.org) שימש להמרת שמות היעד שנבחרו לסמלי גנים סטנדרטיים.
היעדים החופפים בין יעדי התרופות לאוטיזם זוהו באמצעות כלי Draw Venn Diagram המקוון (http://bioinformatics.psb.ugent.be). מטרות חופפות אלו יובאו למסד הנתונים STRING (https://string-db.org/) לבניית רשת אינטראקציה בין חלבון לחלבון (PPI)19. המין הוגדר כ-Homo sapiens, ואינטראקציות חלבון-חלבון עם ציוני ביטחון מתחת ל-0.40 סיננו החוצה. הרשת שהתקבלה יובאה לתוכנת Cytoscape 3.10.0 לניתוח חזותי16, כאשר תרופות, תרכובות ומטרות יוצגו על ידי יהלומים אדומים, עיגולים כחולים ומשולשים ירוקים, בהתאמה, ומשקלי הקצוות שיקפו את דרגת המרכזיות של הצמתים. תוסף CytoNCA שימש לאחר מכן לחישוב ערכי מרכזיות דרגות, לדרג את תרכובות הליבה ולזיהוי תרכובות הליבה העיקריות. ניתוח העשרת ביואינפורמטיקה של הגנים המטרה בוצע באמצעות פלטפורמת Metascape20, כולל ניתוח GO (תהליך ביולוגי, BP; תפקוד מולקולרי, MF; רכיב תאי, CC) וניתוח מסלולי KEGG.
עגינה מולקולרית
חלבוני הקולטנים הקשורים לאוטיזם שזוהו לעיל עובדו מראש לפני העגינה על ידי תיקון שאריות חסרות, אופטימיזציה של מצב הפרוטונציה והסרת מולקולות מים גבישיות לשמירה על שלמות מבנית של הקולטן. שיטת העגינה החצי-גמישה במודולCDOCKER 21 שימשה לאחר מכן. שאריות בטווח של 10 Å מהליגנד הקו-קריסטלי הוגדרו ככיס הפעיל, וקשירת הקולטנים ליגנד דומה באזור זה.
תוצאות העגינה הוערכו תוך שימוש באנרגיית האינטראקציה של CDOCKER כמדד המפתח, כאשר ערכי אנרגיה נמוכים יותר מצביעים על קשירה יציבה יותר בין הליגנד לקולטן. אמינות מצב העגינה הוערכה על ידי השוואת קונפורמציית הקישור המרחבית של רכיב הליבה לזו של הליגנד הקו-קריסטלי בכיס הפעיל. הקומפלקס עם אנרגיית הקשירה האופטימלית והאינטראקציה החזקה ביותר נבחר כקונפורמציה הראשונית לסימולציה הדינמיקה המולקולרית הבאה.
סימולציית דינמיקה מולקולרית
בהתבסס על ממצאי העגינה המולקולרית, בוצעה סימולציה של דינמיקה מולקולרית כדי לחקור את מנגנון הקשר בין קומפלקס הקוורצטין ל-TNF-α. גישה זו מדמה תנועה מולקולרית ואינטראקציות ברמה האטומית ומנתחת שינויים דינמיים בחלבונים, ליגנדים והסביבה הסובבת, ובכך מספקת מידע על שינוי קונפורמציה מולקולרי, יציבות קישור ודינמיקת חלבון-ליגנאד. קומפלקס הקולטן-ליגנד הומס באמצעות מודל המים TIP3P, עם מרחק בופר של לפחות 12 Å בין הקומפלקס לגבול המערכת, כדי להבטיח פתרון מלא ולהפחית את השפעות הגבול. הריכוז היוני הוגדר ל-0.154 מ', ויוני Na⁺ ו-Cl⁻ נוספו כדי לנטרל את מטען המערכת. לשיפור דיוק הסימולציה, אומץ שדה כוח חלבוני Amber14SB22 , המתאר ביעילות אינטראקציות בין-מולקולריות לא קשורות ומצבי קישור ומתאים במיוחד למחקרים מורכבים של חלבון-ליגנד23.
למזעור אנרגיה ראשוני, שיטת הירידה התלולה ביותר הופעלה במשך 5000 צעדים להסרת מגעים לא סבירים וקונפורמציות באנרגיה גבוהה, תוך שימוש בסף התכנסות של 10 kJ/mol/nm להבטחת הרפיית פרמטרי שדה הכוח. שיטת הגרדיאנט הצמוד הופעלה לאחר מכן במשך 2,000 שלבים כדי לאופטימיזציה נוספת של המצב התרמודינמי ולהבטיח יציבות מערכת24.
במהלך שיווי משקל, בוצעה לראשונה סימולציית אנסמבל NVT ל-100 ps, עם צעד זמן של 2 fs. המערכת חוממה בהדרגה ל-300 K כדי להפחית את השפעת המבנה ההתחלתי ולהגיע לאיזון תרמודינמי. האנסמבל הועבר לאחר מכן ל-NPT, ובוצעה איזון נוסף של 100 פס"ס בלחץ קבוע של בר אחד לייצוב הצפיפות והלחץ. שלב זה שימש להבאת המערכת למצב תרמודינמי יציב לפני סימולציית ייצור25. הסימולציה הפורמלית של דינמיקה מולקולרית בוצעה למשך 20 ננו-שניות, כאשר הטמפרטורה נשמרה על 300 קלווין, לחץ של בר אחד, וצעד זמן של 2 fs. המסלול נשמר כל 10 ps. כדי להבטיח יציבות ודיוק סימולציה, פרמטרים פיזיקליים, כולל טמפרטורה, לחץ ונפח, נבדקו באופן קבוע כדי לוודא שהם נשארים בטווחים הצפויים.
סריקה גמישה של אלנין
בהתבסס על הקונפורמציה היציבה שהושגה מסימולציית דינמיקה מולקולרית, בוצעה סריקת אלנין על כל שאריות חומצות האמינו ברדיוס של 3 Å מממשק קישור הליגנד. בהליך זה, שאריות המטרה מוחלפות באופן שיטתי באלנין, ובכך חותכות את שרשרת הצד תוך שמירה על מבנה השרשרת הראשית והסרת אינטראקציות ספציפיות המתווכות בשרשרת הצדדית. משקל התרומה של כל שארית לאפיניות הקישור נמדד על ידי חישוב השינוי באנרגיית הקישור החופשית בין הקומפלקסים מסוג פרא למוטנטים. בניגוד למודלים סטטיים מסורתיים, מחקר זה הציג מנגנון הרפיה גמיש בשרשרת צדדית, שאפשר לסביבה סביב אתר המוטציה לעבור הרפיה מבנית ולדמות בצורה ריאליסטית יותר את התגובה הדינמית של ממשק הקישור. ניתוח זה נועד לזהות שאריות מועמדות לנקודות חמות השומרות על יציבות מורכבת, ולספק טביעת אצבע אנרגטית לאופטימיזציה של תרכובות עופרת המיועדות לחלבונים הקשורים לאוטיזם.
ניסויים בבעלי חיים
בעלי חיים ניסיוניים
נבחרו חולדות Sprague-Dawley (SD) בריאות בדרגת SPF (שלושה זכרים ושלוש נקבות, בני 3 חודשים), שנולדו וגודלו בתנאים זהים. טמפרטורת החדר נשלטה ב-18–22 מעלות צלזיוס, הלחות היחסית נשמרה ב-60%–70%, ומחזור האור היה 12 שעות: 12 שעות (אור: חושך). כל הניתוחים הניסיוניים בבעלי חיים אושרו על ידי ועדת האתיקה הניסיונית של בית החולים הראשון של העם של זוני (מספר אישור: LunShen (2025)-2-362).
הזדווגות בעלי חיים, זיהוי הריון וקיבוץ
כל העכברים הוזנו באופן אדפטיבי בסביבה עם SPF במשך שבוע לאחר הרכישה. עכברוש נקבה וזכר אחד היו כלואים יחד בשעה 18:00 בכל אחר הצהריים. בדיקת פקק וגינלי בוצעה בשעה 8:00 בבוקר למחרת (12 שעות לאחר הכלוב). נוכחות פקק וגינלי נחשבה להזדווגות מוצלחת, ואותו יום סומן כיום הריון 0.5 (GD0.5). חולדות בהריון אוחסנו בכלובים נפרדים. משקל הגוף של חולדות בהריון נמדד ונועד מדי יום. משקל הגוף של חולדות בהריון עלה ברצף, עם עלייה יומית ממוצעת של 2–5 גרם, ומשקל הגוף הכולל יכול היה לעלות בכ-30% לפני הלידה. לאחר כ-10 ימי הריון, הייתה בולטת אסימטרית טיפוסית בצורת אגס של הבטן, וניתן היה להרגיש גושים קשים של עוברים במישוש, שהובחנה מהיקף הבטן האחיד והרך שנגרם מהשמנה.
הקמת מודלים והקמת קבוצת בקרה:
בקבוצת המודל, שתי חולדות בהריון נבחרו באקראי והוזרקו לתוך-פריטונית בתמיסת VPA (600 מ"ג/ק"ג) פעם אחת ביום ההריון 12.5 (GD12.5). VPA ניתנה באמצעות הזרקה תוך-פריטוניאלית אחת במינון של 600 מ"ג/ק"ג ביום ההריון ה-12.5. משטר זה נבחר בהתבסס על מחקר פורץ דרך של שניידר ופרז'בלוצקי26, שקבעו כי חשיפה ל-VPA בנקודת זמן הריון ספציפית זו משחזרת הן את התכונות הנוירואנטומיות והן ההתנהגותיות של ASD בבני אדם. פרוטוקול זה הפך מאז למודל הסטנדרטי ואומת בעקביות במחקרים פרמקולוגיים אחרונים באמצעות פרמטרים זהים27. בקבוצת הביקורת הריקה, נבחרה חולדה בהריון אחת והוזרק תוך-פריטונית בנפח שווה של 0.9% תמיסת מלח תקינה באותו זמן. חולדות בהריון ילדו באופן טבעי, ויום הלידה של הצאצא נרשם כיום 0 לאחר הלידה (PND0). כל הצאצאים של העכברים גמלו והוחזקו בנפרד לפי מגדר ב-PND21.
קיבוץ הצאצאים והתערבות
ב-PND28, 12 גורים זכרים נבחרו באקראי מצאצאי חולדות הרות שנחשפו ל-VPA (6 בקבוצת המודל ו-6 בקבוצת ההתערבות של הקוורצטין), ו-6 גורים זכרים נבחרו באקראי מצאצאי חולדות הרות שנחשפו למלחים (6 בקבוצת הריק). הובטח שלא יהיה הבדל משמעותי במשקל הגוף בין הגורים בכל קבוצה. צאצאים שנולדו לעכברים בהריון שטופלו במי מלח תקינים הוקצו לקבוצת הריק, אלו שנולדו לעכברים בהריון שטופלו ב-VPA לקבוצת המודל, ואלו שנולדו לחולדות הרות שטופלו ב-VPA לקבוצת ההתערבות של קוורצטין.
התערבות רציפה בוצעה במשך 4 שבועות החל מ-PND28. בקבוצת ההתערבות בקוורצטין, תוספי קוורצטין ניתן תוך-קיבה בשעה קבועה כל יום במינון של 100 מ"ג/ק"ג/יום. מינון זה נבחר בהתבסס על הראיות המשולבות הבאות: (i) מחקר קודם שערך מינון זיהה 100 mg/kg כמנה האופטימלית להקלה על התנהגויות חרדה ולהפחתת ציטוקינים פרו-דלקתיים במודל עכברוש נוירודלקתי המושרה על ידי LPS28; (ii) קוורצטין במינון 50 מ"ג/ק"ג הוכח כמונע ליקויים באינטראקציה חברתית ונזק מוחי חמצוני בעכברוש אוטיזם שנגרם על ידי VPA במודל29; ו-(iii) קוורצטין דרך הפה הוכח לאחרונה כמפחית רמות TNF-α במוח ומשפר התנהגויות דמויות אוטיזם בדגם30 בעכבר אוטיסטי שנגרם מחומצה פרופיונית. ביחד, אימות עצמאי זה תומך בבחירת 100 מ"ג/ק"ג כדי להבטיח מעורבות איתנה של מסלול הדלקתי המתווך על ידי TNF-α במודל ה-ASD המושרה על ידי VPA לאחר הלידה הנוכחי. בקבוצות הבלנק והמודל ניתנה תמיסת תמי מלח רגילה של 0.5% CMC-Na באופן תוך-קיבה בכל יום. לכל בעלי החיים הייתה גישה חופשית למזון ומים במהלך תקופת ההתערבות, ומשקל גופם נמדד מדי שבוע כדי להתאים את נפח הטיפול בהתאם למשקל גופם.
מבחן שדה פתוח
לאחר 4 שבועות של התערבות (בערך PND56), בוצעה בדיקת שטח פתוח להערכת פעילות ספונטנית ורמות חרדה. מבחן השדה הפתוח (OFT) הוא ניסוי התנהגותי קלאסי להערכת התנהגות הקשורה לחרדה בבעלי חיים ניסיוניים. המדדים שנמדדו היו יכולת התנועה הספונטנית של חולדות בסביבה פתוחה והזמן שבילו במרכז השדה הפתוח. מכשיר ה-OFT לחולדות היה בגובה 30 ס"מ, אורך 50 ס"מ ורוחב 50 ס"מ בתחתית, עם קירות פנימיים לבנים, וחולק באופן מלאכותי ל-16 רשתות קטנות, כולל 4 רשתות באזור הפנימי ו-12 רשתות באזור החיצוני. אתר הניסוי נשמר בשקט כדי למנוע גירוי קולי שעלול להשפיע על דיוק תוצאות הניסוי. כל עכברוש הונח במרכז תחתית הקופסה, והקלטת וידאו ותזמון בוצעו בו זמנית. שדה הראייה של המצלמה כיסה את כל השדה הפתוח ותיעד את התנועה הספונטנית של העכברים ואת מספר המעברים בין הרשתות. כל בדיקה נמשכה 5 דקות, ולאחר מכן הופסקה ההקלטה. הקיר הפנימי והתחתית של הקופסה הפתוחה נוקו ב-75% אלכוהול כדי למנוע מהצואה וריח הגוף שנותר על ידי בעל חיים אחד להשפיע על תוצאות הבדיקות של החיה הבאה. הניתוח חזר על עצמו לאחר החלפת העכברים עד שכל העכברים השלימו את הבדיקה.
שיטת זיהוי פתולוגית
לאחר קיבוע עם 4% פאראפורמלדהיד, רקמת המוח עברה התייבשות גרדיאנטית באמצעות מייבש אוטומטי מלא: 75% אתנול לשעתיים, 85% אתנול לשעה, 95% אתנול לשעה, ואתנול מוחלט I-IV למשך 20 דקות כל אחד. לאחר מכן ניקו את הרקמות עם חומר ניקוי I למשך 25 דקות וסוכן ניקוי II למשך 30 דקות, ואחריהם הטמעת פרפין. חתכים בעובי 5 מיקרון פורקו עם תמיסת הסרת שעווה I ו-II למשך 30 דקות כל אחד והודרו מחדש באתנול גרדיאנטי. החלקים צובקו בהמטוקסילין במשך 5–10 דקות, הובדלו באלכוהול חומצה כלורית למשך 3 שניות, וצבעו במים אלקליים כדי להפוך לכחול. לאחר מכן צובעו בצבע נגד באאוזין מסיס באלכוהול למשך 3 דקות, יובשו באתנול גרדיאנטי, נוקו, הוסרו עם מדיום התקנה ניטרלי, ונצפו תחת מיקרוסקופ.
זיהוי רמות TNF-α בסרום וברקמת המוח
לאחר המתן האחרון, נאסף דם אאורטה בטני מחולדות וצנטריפוגה ב-3,000 r/min בטמפרטורה של 4°C למשך 15 דקות, ונאסף הסופרנטנט בסרום. בינתיים, ננתקו רקמות היפוקמפוס וקורטיקלית, והוספה PBS מקוררת מראש ביחס של 1:9 להומוגניזציה מכנית. ההומוגנט עבר צנטריפוגה ב-12,000 r/min בטמפרטורה של 4°C למשך 20 דקות כדי לאסוף את הסופרנטנט. שיטת BCA שימשה לכימות חלבוני רקמה. בהתאם להוראות ערכת ELISA, הדגימות לבדיקה והנוגדנים הביוטיניליים נוספו למיקרופלטה ברצף. לאחר דגירה ושטיפה, נוסף ודוגר קומפלקס אבידין-פרוקסידאז. לאחר שטיפה נוספת, נוספה תמיסת המצע להתפתחות צבע, התגובה הסתיימה בתמיסת עצירה, והספיגה נמדדה ב-450 ננומטר באמצעות קורא מיקרופלט. ריכוזי TNF-α בסרום (pg/mL) וברקמת המוח (pg/mg prot) חושבו באמצעות עקומת הסטנדרטית.
ניתוח סטטיסטי
פונקציית RANDBETWEEN ב-Microsoft Excel שימשה לאקראיות את העכברים וליצירת מספרים אקראיים להקצאת חולדות. כל הנתונים הניסיוניים במחקר זה הובאו לידי ביטוי כסטיית תקן ממוצעת ± (SD). לפני ניתוח סטטיסטי, בוצעו תחילה בדיקות נורמליות (מבחן שפירו-וילק) ובדיקת הומוגניות שונות (מבחן לוון) על נתוני כל קבוצה. עבור נתונים שתאמו את התפלגות הנורמלית וההומוגניות של השנות, נעשה שימוש בניתוח חד-כיווני של שונות (ANOVA חד-כיווני). אם ההבדל היה מובהק סטטיסטית, שיטת LSD שימשה בהשוואות מרובות נוספות. ערך של *p < 0.05, **p < 0.01, ו-***p < 0.001 נחשב מובהק סטטיסטית.