KIMMDY: אמולטור תגובות ביומולקולרי

צפייה במולקולות כשהן יוצרות ושוברות קשרים

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

דרך חדשה לחקות תגובות כימיות

סימולציות מולקולריות מסורתיות מצטיינות בהצגת תנועתן וקיפולן של מולקולות ביומולקולריות גדולות כמו חלבונים ו‑DNA, אך הן מתייחסות לקשרים כימיים כבלתי שבירים. סימולציות תגובתיות קיימות, אך דרישות החישוב שלהן כה גבוהות שלרוב אינן מאפשרות הגעה לזמני אמת הרלוונטיים לביולוגיה, במיוחד כשמעורבים סוגי ואתרי תגובה רבים. KIMMDY (קיצור של Kinetic Monte Carlo Molecular Dynamics) נוקט גישה שונה. הוא מאפשר לסימולציות המסורתיות לטפל בתנועה המולקולרית, ואז משתמש במודול נוסף כדי לחפש מקומות שבהם תגובות עלולות להתרחש, לאמוד כמה מהר ייתכנו, ולבחור באקראי איזו תגובה מתבצעת לאחר מכן. באמצעות מחזוריות של תנועה, חיפוש תגובה ובחירת תגובה, KIMMDY עוקב אחר קסדות תגובה מורכבות על פני זמנים הנעים מחלקי שנייה ועד זמן רב בהרבה, בלי לחשב מסלול תגובה מלא ברמת הקוונטים.

איך האמולטור של תגובות עובד

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

מעקב אחר נזק רדיקלי במסגרת תבניות חלבון

כדי לבדוק ש‑KIMMDY פועל כפי שמצופה, המחברים החילו אותו תחילה על קטגוריה פשוטה של תגובות ברדיקלים אורגניים קטנים שבהם קיימים נתונים ניסיוניים. האמולטור שחזר אילו העברות מימן הן הסבירות ביותר, גם אם הקטין את מהירותן המוחלטת. לאחר מכן הם עברו למקרה מסובך הרבה יותר: קולגן, החלבון הסיבי המעניק לרקמות החיבור את חוזקן. תחת מתח מכני, קשרים בקולגן יכולים להיפרד בצורה הומוליטית, ויוצרים רדיקלים מאוד תגובתיים שמדדוורים מאתר לאתר על ידי העברת אטומי מימן לאורך שרשרת. בסיב קולגן המכיל 2.6 מיליון אטומים, KIMMDY עקב אחר מאות העברות כאלה והראה כי רדיקלים יכולים להירתק מספר ננו־מטרים מהמקום שבו נוצרו. הסימולציות חשפו כי צילוב יוצא דופן (המכונה PYD) וחומצת אמינו שעברה שינוי (DOPA) פועלים שניהם כמלכודות רדיקליות מצוינות, מסקנה שנתמכה על ידי חישובים תרמודינמיים ועל ידי פרשנות מחודשת של ספקטרות תהודה ספין אלקטרוני מוקדמות.

כאשר תגובות שונות מתחרות

מערכות ביולוגיות רבות יכולות לשבור את אותו הקשר ביותר מדרך אחת, מה שמוביל או להיווצרות רדיקלים או למוצרים "סגורים‑מעטפת" קונבנציונליים יותר. בקולגן, למשל, קשר פפטידי יכול להישבר באופן סימטרי לשני רדיקלים (הומוליזה) או להגיב עם מים ולהישבר באופן לא סימטרי (הידרוליזה). סימולציות קוונטיות ישירות של שתי האפשרויות בסביבה ריאליסטית יהיו בלתי אפשריות מבחינת עלות חישובית. עם KIMMDY, המחברים שילבו מודל פיזיקלי של הומוליזה מסייעת‑כוח עם מודל אירגוני המבוסס על ניסויים להידרוליזה, ותירגמו את התחרות בין שתי התגובות בשני מצבים: פפטיד בודד ו‑פיבריל קולגן צפוף. הם מצאו שבעוד שבשרשרת מבודדת תחת מתיחה ההידרוליזה דומיננטית בבירור, בפיבריל צפוף ומחובר צומת המתח המכני מרוכז באזורים מסוימים כל כך עד שההומוליזה הופכת לתחרותית ואף מהירה יותר, ומסבירה מדוע תצפיות של יצירת רדיקלים נצפות ברקמות אמיתיות.

נזק מונע‑אור וצילובים ב‑DNA

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

מדוע זה חשוב לביולוגיה ולעיצוב

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

ציטוט: Hartmann, E., Buhr, J., Riedmiller, K. et al. KIMMDY: a biomolecular reaction emulator. Nat Commun 17, 3500 (2026). https://doi.org/10.1038/s41467-026-71955-2

מילות מפתח: סימולציות ביומולקולריות, קינטיקה של תגובות, מודלינג מולקולרי, נזק ל-DNA, מכנוכימיה