$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
השלד הציטולוגי מורכב מחוטים בתוך התא ומולקולות נלוות כגון מנועים מולקולריים, המהווים לרוב רשת דינמית בעלת תכונות מכניות יוצאות דופן. השלד הציטולוגי קיים בתצורות שונות בסוגי תאים שונים כמעט בכל צורות החיים. תפקודו הנכון חיוני לתהליכים תאיים בסיסיים כמו חלוקה, תנועתיות וקיטוב. הוא גם שולט באינטראקציות מכניות בין תאים, ובכך משפיע על המורפוגנזה של רקמות ואורגניזמים. השלד הציטולוגי עומד בבסיס מספר פונקציות ומתבטא בתהליכים ביולוגיים רבים. לדוגמא, התכווצות השרירים קשורה לשבץ הכוח של מנועים מולקולריים של מיוזין על חוטי אקטין. דוגמה נוספת היא תחזוקה של נוירונים, המסתמכת על תנועות של מנועי קינזין לאורך מיקרו-צינורות הממוקמים בתוך האקסונים של נוירונים אלה. אקטין ומיקרו-צינורות הם שני סוגים בולטים של חוטי שלד תא, ובלעדיהם, החיים כפי שאנו מכירים אותם יהיו בלתי אפשריים.
השלד הציטולוגי הוא בעצם מערכת ביומכנית, שלא ניתן לצמצם אותה רק לכימיה שלה. מיקרו-צינורות או חוטי אקטין בנויים מאלפי מונומרים ומשתרעים על פני מספר מיקרומטרים. הקונפורמציות של חוטים אלה בחלל והכוחות שהם יכולים להעביר לקרום הפלזמה, לגרעין או לאברונים אחרים הם היבטים מרכזיים בתפקידם בתא. לדוגמה, רשת של חוטי אקטין ומנועי מיוזין, הנקראת אקטומיוסיןקליפת המוח 1, מייצרת כוחות לשמירה על תנועתיות התאים ושינויים מורפולוגיים בתאי בעלי חיים. סידור שונה מאוד נראה בתאי צמחים, כאשר מיקרו-צינורות קליפת המוח מכוונים את שקיעת סיבי התאית, ובכך שולטים בארכיטקטורת דופן התא, מה שבסופו של דבר קובע כיצד תאים אלה יגדלו בעתיד2.
בעוד שמכניקה ממלאת תפקיד משמעותי בפעולות השלד התא, כימיה חשובה לא פחות. חוטים גדלים בתהליך הרכבה עצמית לפיו מונומרים מוצאים את אתר העגינה שלהם בקצה החוט לאחר פיזור דרך הציטופלזמה3. בקנה מידה מולקולרי, הרכבה ופירוק בקצה החוטים נקבעים, אם כן, על ידי זיקות מולקולריות4. באופן דומה, חלבונים של השלד הציטולוגי מפוזרים, ושיעורי קשירה וביטול קשירה קובעים את זיקתם לחוטים שהם נתקלים בהם. במקרה של מנועים מולקולריים, מחזורים של תגובות כימיות הכוללות הידרוליזה של ATP קשורים לתנועות לאורך החוטים ואולי לכוחות המלווים אותם5. למרבה הפלא, השלד הציטולוגי מציע אתגרים יוצאי דופן רבים ומגוון גדול של תהליכים הכוללים רכיבים דומים. זהו מגרש משחקים עשיר בממשק שבין ביולוגיה, כימיה ופיזיקה.
מערכות שלד תא ניתנות למודלים מתמטיים. למעשה, הודות למחקר מצוין שנעשה בעשורים האחרונים, סביר להניח שהמרכיבים המולקולריים העיקריים כבר זוהו, כפי שמודגם באנדוציטוזיס6. באורגניזמים מודלים, כמו שמרים, התכונות של היסודות האלה ידועות, כמו גם הרכב המערכת של חלק מהתהליכים שלהם. לדוגמה, המבנה ותכונות החומר של מיקרו-צינורות7, כמו גם מספרם ואורכם הממוצע בשלבים שונים של הציר המיטוטי, תוארו8. מספר הקינזינים המחברים מיקרו-צינורות למבנה מכני קוהרנטי ידוע לעתים קרובות9. המהירויות של מנועים רבים נמדדו במבחנה10. בנוסף, ניסויים יכולים לצפות ולכמת את המערכות הללו in vivo בתנאים פראיים או מוטציות. שילוב תיאוריה לצד ניסויים in vivo ו-in vitro מאפשר לחוקרים לבחון אם הידע הנוכחי על מערכת שלד תא מספיק כדי להסביר את התנהגותה הנצפית. השימוש בכלים מתמטיים וחישוביים מאפשר לנו גם להסיק מסקנות כיצד רכיבים פועלים באופן קולקטיבי על בסיס הנחות הנגזרות מתצפיות בקנה מידה מולקולרי, בדרך כלל במצבים פשוטים (למשל, ניסויים של מולקולה בודדת).
ניתן להמחיש את תפקידה של התיאוריה באמצעות דוגמה מעשית: הכאת ריסים. מכה זו נובעת מתנועת מנועי דינאין לאורך מיקרו-צינורות בריסים. ניתן לשאול מה קובע את מהירות מנוע הדיניין במערכת זו. תשובה אפשרית אחת היא שהמהירות המרבית מוגבלת על ידי הדרישה לשמור על דפוס פעימה מסוים. זה היה מובן אם המכות היו תחת ברירה טבעית. במקרה כזה, אם המנועים ינועו מהר יותר, התהליך יאבד את התכונות הרצויות שלו - הריסים לא יפעמו באותה יעילות או אפילו ייכשלו לחלוטין. למרות שזה אפשרי, חלופה שנייה היא שגורם מהותי כלשהו יכול להגביל את מהירות הדינין.
לדוגמה, ייתכן שלתא אין מספיק ATP כדי להפוך את הדיניין למהיר יותר, או שפשוט לא ניתן להאיץ את תנועות החלבון הנדרשות לפעילות הדיניין. במקרה כזה, אם ניתן היה להפוך את המנועים למהירים יותר למרות המגבלות הפיזיות, המכות היו משתפרות. אפשרות שלישית, כמובן, היא ששינוי המהירות אינו משפיע על התהליך באופן משמעותי, מה שעשוי להועיל לאורגניזם בכך שהוא מספק "חוסן" מסוים כנגד גורמים בלתי נשלטים. מבין שלוש האפשרויות הללו, ניתן לזהות את הנכונה על ידי חישוב דפוס הפעימות מתכונות הדינין. אכן, מודל מתמטי מתאים אמור לחזות כיצד דפוס הפעימות מושפע ממהירות משתנה של דינין ואינו כפוף למגבלות הקיימות בעולם הפיזיקלי. מטבע הדברים, יש לאמת את תוקף המודל, אך גם מודלים "שגויים" יכולים לייצר רעיונות מעניינים.
המודל יכול ללבוש צורה של מסגרת אנליטית או להיות סימולציה מספרית של המערכת. כך או כך, הפער בין הסקאלה המולקולרית לסולם הפונקציונלי נותר מכשול, ופיתוח מודלים אלה אינו משימה פשוטה, שכן יש לשלב מספר תהליכים מכניים וכימיים במשוואות המתארות את המערכת הביולוגית. התיאוריה מגיעה בצורות שונות, ומציעה פשרות שונות בין פשטות לריאליזם. הגדלת מידת הפרטים במודל אינה תמיד מועילה מכיוון שהיא עלולה להגביל את יכולתנו לפתור את המשוואות או, במילים אחרות, לגזור את התחזיות של התיאוריה. אותה פשרה קיימת עבור סימולציות. המודלים יצטרכו לבחור את מרכיבי המערכת שיש לקחת בחשבון תוך התעלמות מהיבטים מסוימים. החלטות מפתח אלה יהיו תלויות מאוד במטרת המחקר. כיום, השיפורים יוצאי הדופן בחומרת המחשב מאפשרים לדמות מערכות שלד תא רבות עם מספיק פרטים לאורך זמן מספיק כדי לנתח את ההתאמה שלהן. זה לעתים קרובות ייצר רעיונות בלתי צפויים וכיוונים חדשים במחקר. לדוגמה, סימולציות דומות לאלו שישמשו בפרוטוקול זה הובילו לחישוב מאחורי המעטפה שיכול לחזות את ההתכווצות של רשת על סמך הרכבה11.
שיטות נומריות נמצאות בכל מקום בהנדסה ובמדעי הפיזיקה, והשימוש בהן בביולוגיה הולך וגדל. כיום, כמעט כל האמצעים הטכנולוגיים שלנו (שעונים, טלפונים, מכוניות ומחשבים) נוצרו לראשונה במחשב, וקיימות תוכנות חזקות כדי לעשות זאת. בהינתן מערכת שלד תא מאופיינת היטב ובהנחה שנקבעה רמת תיאור מתאימה, עדיין יש לפתור מספר בעיות לפני שניתן יהיה לדמות אותה. עבור הבעיות הפשוטות ביותר, דרך הפעולה המתאימה ביותר עשויה להיות כתיבת סימולציה "על ידי קידוד מאפס", במילים אחרות, להתחיל בשפת תכנות גנרית או פלטפורמה מתמטית כמו MATLAB. יש לכך יתרון שלמחבר הקוד יהיה ידע אינטימי על מה שיושם ויודע בדיוק כיצד התוכנה עובדת. עם זאת, מסלול זה אינו נטול סיכונים, ואין זה נדיר לראות דוקטורנטים מבלים את רוב זמן עבודתם בכתיבת קוד במקום לעסוק בשאלות מדעיות.
האלטרנטיבה היא להשתמש בתוכנה שנוצרה על ידי אחרים, אבל גם זה לא נטול סיכונים; כל קוד מקור גדול נוטה לרכוש באופן ספונטני תכונות של קופסה שחורה בלתי חדירה, למרות המאמצים הראויים להערצה של מחבריהם למנוע זאת. שימוש בקופסאות שחורות הוא בוודאי לא חלומו של מדען. קוד מקור גדול יכול גם להפוך לנטל, וייתכן שיהיה מהיר יותר להתחיל מאפס מאשר לשנות בסיס קוד קיים כדי לגרום לו לעשות משהו שונה. כדי להקל על בעיה זו, תמיד אפשר להזמין את מחברי התוכנה לעזור, אך ייתכן שזה לא מספיק. לעתים קרובות, יש הבדל בתרבות המדעית בין מחברי התוכנה לבין האנשים שרוצים להשתמש בה, מה שאומר שצריך להבהיר הנחות מרומזות רבות. על ידי הפיכת הקוד לקוד פתוח, צפוי שיותר אנשים יהיו מעורבים בפיתוח התוכנה ולשמור על התיעוד שלה, ובכך לשפר את איכותה. כל אלה הם נושאים חשובים שיש לתת להם את הדעת הראויה לפני ביצוע השקעה כלשהי. עם זאת, הדרך היחידה להתקדם בטווח הארוך היא לקדם פתרונות תוכנה מוצקים, המשמשים ומתוחזקים על ידי קהילה רחבה עם אינטרסים מדעיים משותפים.
למרות שפרוטוקול זה משתמש ב-Cytosim, ישנם כלי קוד פתוח אחרים שעשויים להיות מסוגלים לדמות את אותה מערכת, למשל, AFINES12, MEDYAN13, CyLaKS14, aLENS15 ו-AKYT16, אם להזכיר כמה. למרבה הצער, השוואה בין מיזמים אלה היא מעבר לתחום המאמר. כאן ניתנות הוראות שלב אחר שלב להדמיית רשת אקטומיוזין דו-ממדית מתכווצת. מערכת זו פשוטה ועושה שימוש ביכולות המבוססות יותר של Cytosim. Cytosim בנוי סביב מנוע ליבה חוצה פלטפורמות שיכול להריץ סימולציות בדו-ממד או בתלת-ממד. יש לו בסיס קוד מודולרי, מה שהופך אותו להתאמה אישית בקלות לביצוע משימות מסוימות. ציטוסים יציב ויעיל באותה מידה בתלת מימד ושימש בהצלחה בעבר לחקירת בעיות מגוונות הכוללות מיקרו-צינורות וחוטי אקטין: שיוך של שני אסטרים של מיקרו-צינורות17, תנועת הגרעינים בתאים18,19, אנדוציטוזיס6, ציטוקינזיס20, היווצרות הציר המיטוטי21, תנועות הציר המיטוטי22, לכידת כרומוזומים23, התכווצות רשתות האקטומיוזין 11,24, והמכניקה של טבעת המיקרו-צינורות בטסיות הדם25, והיכולות שפותחו עבור פרויקטים אלה נשמרו בקוד. ניתן להתאים את זרימת העבודה המתוארת כאן לבעיות רבות אחרות. הוא עושה שימוש בשורת הפקודה של יוניקס, שאולי אינה מוכרת לחלק מהקוראים. עם זאת, השימוש בשורת הפקודה הוא הדרך הניידת והנוחה ביותר להפוך את תהליך הפעלת הסימולציות לאוטומטי. ממשקי משתמש גרפיים משולבים שואפים להציע גישה קלה ואינטואיטיבית לתוכנה, אך לעתים קרובות זה בא על חשבון הכלליות. מטרת מאמר זה היא להמחיש גישה שניתן לשנות או להתאים בקלות לבעיות אחרות. הערות ניתנות כדי להסביר את משמעות הפקודות.
כדי לדמות רשת אקטומיוזין, חוטים מודלים כקווים מכוונים ומיוצגים על ידי קודקודים המפוזרים לאורכם (איור 1). זוהי רמת ביניים של תיאור, הנפוצה בפיזיקה של פולימרים, המתעלמת מהאופי התלת-ממדי האמיתי של החוטים אך מאפשרת לחשב את הכיפוף. חוטים עשויים לגדול ולהתכווץ בקצותיהם, בעקבות מודלים שונים המכסים גם פנומנולוגיה של אקטין וגם של מיקרו-צינורות. בתאים, חוטים מאורגנים בעיקר באמצעות אינטראקציות המגבילות את תנועתם, למשל, היצמדות לחוטים אחרים או פשוט כליאה בתוך התא. ב-Cytosim, כל האינטראקציות הללו הן ליניאריות ומשולבות במטריצה גדולה26. המשוואות המתארות את התנועה של כל קודקודי הנימה נגזרות ממטריצה זו, בהנחה של תווך צמיג ומונחים משתנים אקראיים המייצגים תנועה בראונית. משוואות אלו נפתרות מספרית כדי להשיג את תנועת החוטים יחד עם כל הכוחות הפועלים עליהם בצורה עקבית ויעילה26. על גבי מנוע מכני זה, ישנו מנוע סטוכסטי המדמה אירועים בדידים, כגון חיבורים וניתוק של מנועים מולקולריים או דינמיקת הרכבה של חוטים. לסיכום, Cytosim משתמש תחילה בדינמיקה מדומה כדי לחשב את המכניקה של רשת חוטים, המחוברים באופן שרירותי כלשהו, ושנית, שיטות סטוכסטיות כדי לדמות את הקשירה, הניתוק והדיפוזיה של חלבונים המחברים או משפיעים על החוטים.
זרימת העבודה המומחשת כאן בוצעה לעתים קרובות כדי לחקור תחילה מערכת המשתמשת ב-Cytosim. השלב הקריטי עבור משתמשים פוטנציאליים רבים עשוי להיות התקנת רכיבי התוכנה. הפצת התוכנה כקוד מקור ממלאת את הציוויים של המדע הפתוח, אך היא מועדת לטעויות מכיוון שלמפתחי התוכנה יש גישה רק למאגר מוגבל של ארכיטקטורה לבדיקת התוכנית. ההידור עלול להיכשל מכיוון שמערכות ההפעלה שונות. ההוראות המובאות כאן עשויות להתיישן עם התפתחות מערכות המחשב וקודי המקור. לפיכך, בדיקה תקופתית של ההוראות העדכניות ביותר באינטרנט היא חיונית. במקרה של בעיה, מומלץ מאוד שמשתמשים ידווחו בחזרה על ידי פרסום בערוץ המשוב הרלוונטי (כרגע דף הבית של Cytosim ב-Gitlab) כדי לעזור לפתור את הבעיה.