שיפור הנדסת הביו-אינפורמטיקה באמצעות ניצול התכונות התרפויטיות הגרפיות של תרופות אנטיטוקסין המאושרות קלינית במחלות זואונוטיות

מדוע זה חשוב עבור תרופות בעתיד

אנתראקס נתפס לעיתים כמחלת חוות ישנה, אך הוא עדיין מאיים על משקים, על אנשים באזורים מתפתחים ואף על ביטחון לאומי מכיוון שנבגיו ניתנים להשתמש ככלי נשק. הטיפול באנתראקס מסתמך על קבוצת קטנה של אנטיביוטיקה חזקה ותרופות אנטיטוקסין. המחקר הזה שואל שאלה שנראית פשוטה אבל בעלת השלכות רחבות: האם ניתן להשתמש ברעיונות מתוך מתמטיקה של רשתות — תורת הגרפים — כדי לקרוא את הצורות של תרופות אלה ולחזות במהירות תכונות פיזיקליות מרכזיות המשפיעות על התנהגותן בגוף? אם כן, אותה גישה עשויה להאיץ את החיפוש אחרי תרופות אנטיטוקסין מדור הבא.

מנבגים קטלניים לתכניות תרופות

אנתראקס נגרם על ידי Bacillus anthracis, חיידק היוצר נבגים עמידים ומדביק בעיקר בעלי חיים מרעיים כמו בקר וכבשים, אך יכול גם לעבור לבני אדם. ברגע שנכנס לגוף יכול להתקיף את העור, הריאות או המעי, ולהוביל מתסמינים דמויי שפעת למחלה פתאומית ומסכנת חיים. כיום רופאים מסתמכים על סט מוגבל של תרופות שאושרו קלינית — כגון ציפרופלוקסצין, תרכובות מבוססות פניצילין וכן מספר אנטיביוטיקה דמוית טטרציקלין — כדי לעצור את הזיהום ואת הרעלנים שהוא מייצר. כל אחת מהתרופות הללו היא סבך תלת־ממדי של אטומים. האופן שבו האטומים מחוברים קובע תכונות בסיסיות כמו המסה המולקולרית, כמה בקלות גלולה מתמוססת במים, וכיצד היא נעה דרך הגוף. הבנת הקישורים האלה בין מבנה להתנהגות חיונית לשינוי בטוח של תרופות ישנות או לפיתוח תרופות חדשות.

הפיכת מולקולות לרשתות של נקודות וקווים

המחברים משתמשים בתורת הגרפים הכימית, תחום המתייחס למולקולה כרשת: אטומים הופכים לנקודות וקשרים הופכים לקווים. מהרשת הזו הם מחשבים סכומים מספריים הנקראים אינדקסים טופולוגיים, התופסים דפוסים של הסתעפות וחיבוריות. בעבודה זו הם מתמקדים במשפחה של אינדקסים "הפוכים" — וריאנטים של מדדים קלאסיים כמו אינדקסי זגרב והאינטגרציה אטום–קשר — ובאובייקט אלגברי קומפקטי הנקרא פולינום M, שמקודד עד כמה דפוסי חיבור שונים מיוצגים במבנה של התרופה. עבור ארבע־עשרה תרופות הקשורות לאנתראקס שאושרו קלינית הם השתמשו בכלי פייתון (RDKit ו‑NetworkX) לבניית הגרפים, חישוב פולינומי M, וגזירת פאנל של אינדקסים הפוכים. לאחר מכן אספו תכונות שנמדדו בניסויים מתוך מאגר PubChem: נקודות התכה ורתיחה, מסיסות במים, מסה מולקולרית, צפיפות ותיאורים נוספים שמשפיעים על איך תרופה נספגת, מופצת, מתמטבת ומופרשת.

התאמת דפוסי מבנה לתכונות מדידות

עם שתי סדרות המספרים ביד — האינדקסים המבוססים על גרפים והתכונות מהעולם האמיתי — הצוות ערך ניתוח רגרסיה, שיטה להתאמת עקומות מתמטיות לנתונים. בהנחיית צורות פני השטח של פולינום M הם בחנו שתי משפחות של מודלים: עקומות לוגריתמיות, שגדלות במהירות ואז מתמזגות, ועקומות קוביות, שיכולות להתעקל בצורה דרמטית יותר. עבור כל אינדקס וכל תכונה הם העריכו עד כמה עקומה נתונה מסבירה את הנתונים (באמצעות סטטיסטיקת R² המוכרת) ועד כמה החיזוי יציב כאשר ויתרו שיטתית על נקודות נתונים (באמצעות מדד חוצה‑אימות מחמיר בשם Q²).

מה הצליח, מה נכשל, ולמה זה חשוב

הדפוס הבולט ביותר היה שלא כל התכונות ניתנות ללמידה באופן שווה מתוך המבנה. עבור תכונות תרמודינמיות כמו נקודת התכה או מדד אהדת־שומן LogP, האינדקסים ההפוכים הופיעו כחלשים: ערכי ה‑R² של המודלים נותרו נמוכים, מה שמעיד על חוסר כוח חיזוי ממשי. לעומת זאת, מספר אינדקסים — במיוחד מדד שנקרא אינדקס זגרב שני מותאם (mM2) ואינדקס הפוך של חיבור אטום–קשר — הראו קשרים חזקים מאוד עם המסה המולקולרית, מדד יסוד לגודל המולקולרי. מודל לוגריתמי פשוט שקישר בין mM2 למסה המולקולרית השיג גם התאמה גבוהה וגם חיזוי חוצה‑אימות חזק (R² בערך 0.97 ו‑Q² סביב 0.99), אפילו אחרי שהמחברים בדקו בקפדנות בעיות של התאמה יתרה ומקריות באמצעות מבחני השמטת‑אחד, ניתוח תחום‑הישימות ו‑Y‑randomization. מודלים קוביים מורכבים יותר התאימו כמעט טוב מדי לנתונים הקיימים אך נכשלו בבדיקות היציבות האלה, מה שמדגים כמה ערכות נתונים קטנות יכולות להטעות כאשר משתמשים בנוסחאות מסובכות מדי.

כיצד זה מסייע בעיצוב תרופות אנטיטוקסין טובות יותר

עבור לא־מומחים, המסקנה המרכזית היא שהמחברים בנו דרך מהירה ומבוססת־מתמטית לאמוד את המסה המולקולרית של תרופת אנטראקס ישירות מתוך דפוס החיבורים האטומיים שלה — בלי צורך בהת characterized ניסויית מלאה. המסה המולקולרית אינה מדד ליעילות ההרג של החיידק, אבל היא מסנן מרכזי בעיצוב תרופות, קשורה ליכולתו של תרכובת להיספג, להתפזר ולהיות מופרשת מהגוף. בזיהוי אילו אינדקסים מבוססי‑גרף עוקבים בעקביות אחר גודל המולקולה, ואילו לא מנבאים משמעותית תכונות עדינות יותר כמו נקודת התכה או מסיסות שומן, עבודה זו מחדדת את כלי העבודה של עיצוב תרופות בעזרת מחשב. בעתיד, מודלים תורת‑גרפים דומים יכולים לסייע לחוקרים לבצע סריקה מהירה של מאגרי מועמדים גדולים של מולקולות אנטיטוקסין, ולסנן החוצה אלה שהגודל או המורכבות שלהם חורגים מהטווח הרצוי הרבה לפני שמתחילים בבדיקות מעבדה יקרות.

ציטוט: Imran, M., Aqib, M., Malik, M.A. et al. Enhancing bioinformatics engineering by utilizing graph therapeutic properties for clinically approved antitoxin drugs in zoonotic diseases. Sci Rep 16, 8590 (2026). https://doi.org/10.1038/s41598-026-36036-w

מילות מפתח: תרופות לאנתראקס, תאוריית הגרפים הכימית, אינדקסים טופולוגיים, דוגמנות QSPR, חיזוי מסה מולקולרית