Summary
פרוטוקול זה מדגים את זרימת העבודה של פלטפורמת SILICOFCM ליצירה אוטומטית של מודל פרמטרי של החדר השמאלי מתמונות אולטרסאונד ספציפיות למטופל על ידי יישום מודל אלקטרומכני רב-ממדי של הלב. פלטפורמה זו מאפשרת בניסויים קליניים בסיליקו שנועדו להפחית ניסויים קליניים אמיתיים ולמקסם את התוצאות הטיפוליות החיוביות.
Abstract
פרויקט SILICOFCM נועד בעיקר לפתח פלטפורמה חישובית לניסויים קליניים בסיליקו של קרדיומיופתיות משפחתיות (FCMs). המאפיין הייחודי של הפלטפורמה הוא שילוב של נתוני הדמיה ביולוגיים, גנטיים וקליניים ספציפיים לחולה. הפלטפורמה מאפשרת בדיקה ואופטימיזציה של טיפול רפואי כדי למקסם את התוצאות הטיפוליות החיוביות. לפיכך, ניתן להימנע מתופעות לוואי ואינטראקציות בין תרופתיות, למנוע מוות לבבי פתאומי, ולקצר את הזמן שבין תחילת הטיפול התרופתי לבין התוצאה הרצויה. מאמר זה מציג מודל פרמטרי של החדר השמאלי שנוצר באופן אוטומטי מתמונות אולטרסאונד ספציפיות למטופל על ידי יישום מודל אלקטרומכני של הלב. השפעות תרופתיות נקבעו דרך תנאי גבול ספציפיים לזרימת כניסה ויציאה, מדידות א.ק.ג. ותפקוד סידן לתכונות שריר הלב. נתונים גנטיים מחולים שולבו באמצעות התכונה החומרית של דופן החדר. ניתוח ראייה אפיקלית כולל פילוח החדר השמאלי באמצעות מסגרת U-net שהוכשרה בעבר וחישוב המלבן הגובל בהתבסס על אורך החדר השמאלי במחזור הדיאסטולי והסיסטולי. ניתוח תצוגה במצב M כולל גבול של האזורים האופייניים של החדר השמאלי בתצוגת מצב M. לאחר חילוץ הממדים של החדר השמאלי, נוצרה רשת אלמנטים סופיים המבוססת על אפשרויות רשת, והופעלה סימולציה של ניתוח אלמנטים סופיים עם מהירויות כניסה ויציאה שסופקו על ידי המשתמש. משתמשים יכולים לדמיין ישירות על הפלטפורמה תוצאות סימולציה שונות כגון לחץ-נפח, מתח-לחץ, ודיאגרמות זמן עבודה שריר הלב, כמו גם אנימציות של שדות שונים כגון תזוזות, לחצים, מהירות ומתחי גזירה.
Introduction
ההתפתחות המהירה של טכנולוגיות מידע, חבילות תוכנה לסימולציה ומכשור רפואי בשנים האחרונות מספקת הזדמנות לאיסוף כמות גדולה של מידע קליני. יצירת כלים חישוביים מקיפים ומפורטים הפכה, אם כן, לחיונית לעיבוד מידע ספציפי משפע הנתונים הזמינים.
מנקודת מבטם של הרופאים, יש חשיבות עליונה להבחין בין פנוטיפים "נורמליים" לעומת "לא נורמליים" בחולה ספציפי כדי להעריך את התקדמות המחלה, תגובות טיפוליות וסיכונים עתידיים. מודלים חישוביים עדכניים שיפרו באופן משמעותי את ההבנה האינטגרטיבית של התנהגות שרירי הלב בקרדיומיופתיות היפרטרופיות (HCM) ומורחבות (DCM)1. חיוני להשתמש במודל ברזולוציה גבוהה, מפורט ומדויק אנטומית של פעילות חשמלית בלב שלם, המחייב זמני חישוב מסיביים, תוכנות ייעודיות ומחשבי-על 1,2,3. מתודולוגיה למודל לב תלת-ממדי אמיתי פותחה לאחרונה באמצעות מודל חומר אלסטי ואורתוטרופי ליניארי המבוסס על ניסויי הולצפפל, שיכולים לחזות במדויק את שדה הובלת האותות החשמליים והתזוזה בתוך לב4. פיתוח גישות חדשניות של מודלים אינטגרטיביים יכול להיות כלי יעיל להבחנה בין סוג וחומרת הסימפטומים בחולים עם הפרעות מולטיגניות ולהערכת מידת הפגיעה בפעילות גופנית רגילה.
עם זאת, ישנם אתגרים חדשים רבים עבור מודלים ספציפיים למטופל. התכונות הפיזיקליות והביולוגיות של הלב האנושי אינן ניתנות לקביעה מלאה. מדידות לא פולשניות כוללות בדרך כלל נתונים רועשים שמהם קשה להעריך פרמטרים ספציפיים עבור המטופל הבודד. חישוב בקנה מידה גדול דורש זמן רב כדי לרוץ, בעוד מסגרת הזמן הקלינית מוגבלת. יש לנהל את הנתונים האישיים של המטופלים באופן כזה שניתן יהיה לעשות שימוש חוזר במטה-נתונים שנוצרו מבלי להתפשר על סודיות המטופל. למרות אתגרים אלה, מודלים רב-ממדיים של לב יכולים לכלול רמת פירוט מספקת כדי להשיג תחזיות העוקבות מקרוב אחר תגובות חולפות שנצפו, ובכך לספק הבטחה ליישומים קליניים פוטנציאליים.
עם זאת, ללא קשר למאמץ המדעי המשמעותי של מספר מעבדות מחקר ולכמות המשמעותית של תמיכה במענקים, נכון לעכשיו, יש רק חבילת תוכנה אחת זמינה מסחרית לסימולציות לב רב-ממדיות ושלמות, הנקראת SIMULIA Living Heart Model5. הוא כולל סימולציה אלקטרו-מכנית דינמית, גיאומטריה של הלב המעודן, מודל זרימת דם ואפיון מלא של רקמת הלב, כולל מאפיינים פסיביים ואקטיביים, טבע סיבי ומסלולים חשמליים. מודל זה מיועד לשימוש ברפואה מותאמת אישית, אך אפיון החומר הפעיל מבוסס על מודל פנומנולוגי שהוצג על ידי Guccione et al.6,7. לכן, SIMULIA אינו יכול לתרגם באופן ישיר ומדויק את השינויים במאפיינים התפקודיים של חלבון התכווצות שנצפו במחלות לב רבות. שינויים אלה נגרמים על ידי מוטציות וחריגות אחרות ברמה המולקולרית והתת-תאית6. השימוש המוגבל בתוכנת SIMULIA למספר קטן של יישומים בקליניקה הוא דוגמה מצוינת למאבקים של ימינו בפיתוח מודלים רב-ממדיים של לב אנושי ברמה גבוהה יותר. מצד שני, היא מניעה את הפיתוח של דור חדש של חבילות תוכנית רב-ממדיות שיכולות לעקוב אחר ההשפעות של מוטציות מקנה המידה המולקולרי לסולם האיברים.
המטרה העיקרית של אלקטרופיזיולוגיה של הלב היא לקבוע את התפשטות האות בתוך פלג הגוף העליון ואת התכונות של כל התאים 4,5,6. פרויקט SILICOFCM8 מנבא התפתחות מחלות קרדיומיופתיה באמצעות נתוני הדמיה ביולוגיים, גנטיים וקליניים ספציפיים לחולה. היא מושגת באמצעות מודלים רב-ממדיים של המערכת הסרקומרית הריאליסטית, הפרופיל הגנטי של המטופל, כיוון סיבי השריר, אינטראקציה בין נוזל למבנה וצימוד אלקטרופיזיולוגי. ההשפעות של עיוות החדר השמאלי, תנועת המסתם המיטרלי והמודינמיקה מורכבת נותנות התנהגות תפקודית מפורטת של מצבי הלב בחולה ספציפי.
מאמר זה מדגים את השימוש בפלטפורמת SILICOFCM עבור מודל פרמטרי של החדר השמאלי (LV) שנוצר באופן אוטומטי מתמונות אולטרסאונד ספציפיות למטופל באמצעות מודל לב בעל מבנה נוזל עם צימוד אלקטרומכני. ניתוחי תצוגה אפיקלית ותצוגת מצב M של LV נוצרו באמצעות אלגוריתם למידה עמוקה. לאחר מכן, באמצעות מחולל הרשת, מודל האלמנטים הסופיים נבנה באופן אוטומטי כדי לדמות תנאי גבול שונים של המחזור המלא עבור התכווצות LV9. בפלטפורמה זו, משתמשים יכולים לדמיין ישירות את תוצאות הסימולציה כגון נפח לחץ, מתח לחץ ודיאגרמות זמן עבודה שריר הלב, כמו גם אנימציות של תחומים שונים כגון תזוזות, לחצים, מהירות ומתחי גזירה. פרמטרי קלט ממטופלים ספציפיים הם גיאומטריה מתמונות אולטרסאונד, פרופיל מהירות בתנאי זרימת גבול הקלט והפלט עבור LV, וטיפול תרופתי ספציפי (למשל, אנטרסטו, דיגוקסין, mavacamten וכו ').
Subscription Required. Please recommend JoVE to your librarian.
Protocol
הפרוטוקול במחקר זה אושר על ידי רשות המחקר לבריאות של שירותי הבריאות הלאומיים של בריטניה צפון מזרח-טיין ו- Wear South ועדת אתיקה למחקר עם מספר הפניה 18/NE/0318 ב-6 בפברואר 2019 ואומץ על ידי מועצת הביקורת המוסדית של כל מרכז משתתף. המחקר נערך במסגרת עקרונות הפרקטיקה הקלינית הטובה ובעקבות הצהרת הלסינקי. הסכמה מדעת התקבלה מכל הנבדקים המעורבים במחקר. פרטי המטופל נשמרים אנונימיים.
1. זרימת עבודה עבור אולטרסאונד M-mode או תצוגה Apical DICOM ניתוח תמונה וחילוץ פרמטרים LV
הערה: כדי להתחיל פרוטוקול זה, על המשתמש להיכנס לפלטפורמת SILICOFCM8 ולבחור את זרימת העבודה המתאימה (כלומר, ניתוח תמונת אולטרסאונד באמצעות מצב M או תצוגה Apical). זרימת העבודה עבור אולטרסאונד M-mode DICOM ניתוח תמונה וחילוץ פרמטר LV, למשל, כוללת מספר שלבים. השלב הראשון הוא התאמת תבניות. תבניות להתאמה כוללות את כל הגבולות הדרושים הרלוונטיים. יש להדגיש כי עבור כל מערך נתונים חדש, יש לחלץ את התבנית באופן ידני פעם אחת בלבד עבור כל מערך נתונים המתאים למכשיר אולטרסאונד ספציפי. דיאגרמת זרימת העבודה מוצגת באיור 1. האזור התואם לתבנית יופק מהתמונה המנותחת לאחר שאלגוריתם פילוח התמונה היעיל מבוסס הגרף של Felsenszwalb ימצא גבולות "חזקים" התואמים את גבולות המחיצה והקיר של ה- LV. בהתבסס על גבולות אלה והמקומות שבהם קוטרי הלב הם הגדולים ביותר (מתאים לקוטר בדיאסטולה) והקטנים ביותר (מתאים לקוטר בסיסטולה), יחושבו מידות LV שונות. על המשתמש להגדיר אם התצוגה היא במצב M או תצוגה Apical.
- היכנס לפלטפורמה עם שם משתמש וסיסמה. תחת מודול אוכלוסייה וירטואלית, בחר מצב אולטרסאונד M-mode או זרימת עבודה של תצוגה Apical.
- ברשימת זרימות העבודה הזמינות, בחר את זרימת העבודה הממוזגת באולטרסאונד .
- במקטע העלאת קובץ , העלה תמונות וקבצי DICOM המאוחסנים באופן מקומי במחשב המשתמש (איור 2).
הערה: קלטים לתצוגת DICOM של מודול M-mode DICOM ניתוח תמונה ממטופל ספציפי וחילוץ פרמטר LV הם תמונת הבדיקה (המסופקת על-ידי המשתמש בתבנית DICOM המייצגת מצב M או תצוגה Apical) וקבצים סטטיים (איור 1). הקבצים הסטטיים מהווים א) קובץ הקלט.txt , שהוא קובץ תבנית המשמש לדריסת מידות הפלט המחושבות מהמודול (קובצי פלט חדשים נוצרים, בתורם, מקובץ זה), וב) קובץ תמונת התבנית, המשמש להתאמת תבניות (כלומר, template_dicom_GEMS.jpg). - בחר תיקיות פרטיות או תיקיות ציבוריות כתיקיית היעד עבור הקבצים. הקלד את ההערה או ההערה הרצויה במקטע הערות לפני הפעלת זרימת העבודה. בחר את התצוגה Apical ואת תמונות האולטרסאונד LV במצב M וקבצי DICOM שהועלו בעבר (שלב 1.3.) (איור 3).
הערה: כלי סגמנטציה של לב אולטרסאונד מורכב משני תת-מודולים: התצוגה apical ותצוגת מצב M. הגישה לשני המודולים מתבצעת באמצעות הפלטפורמה. הניתוח כולל את חישוב הפרמטרים האופייניים הנראים בחתכים המתאימים. בשני המקרים, הערכים הדרושים נכתבים לקובץ שנלקח כקלט למודל הלב הפרמטרי מצד המשתמש (ראה את הדיון לפרטים). - לחץ על הלחצן בצע . הפלטפורמה האינטראקטיבית מודיעה למשתמש כאשר זרימת העבודה הפועלת מסתיימת.
- לדמיין את הגיאומטריה שנוצרה של LV ישירות על הפלטפורמה (מודל 3D יופיע באופן אוטומטי על המסך, אשר ניתן לסובב באמצעות העכבר). האפשרויות הזמינות כוללות מודלים של 'מוצלל ' ו'מסגרת תיל' לתצוגה חזותית.
הערה: לאחר חילוץ ממדי LV, בהתבסס על אפשרויות רשת שינוי, נוצרת רשת השינוי של אלמנטים סופיים, והדמיית ניתוח אלמנטים סופיים מופעלת עם מהירויות כניסה ויציאה שנקבעו על-ידי המשתמש (ראה להלן).
2. פותר אלמנטים סופיים PAK לסימולציות של מבנה זורם (FS)
הערה: כלי זה יכול לשמש לניתוח אלמנטים סופיים של בעיות נוזל מוצק מצומד. הוא תומך בצימוד חזק ומשוחרר בין מוצק לנוזל. אלמנטים יכולים להיות משושה או טטרהדרונים, עם או בלי צומת נוסף במרכז האלמנט. לפותר זה יש מודלים חומריים מובנים כגון מודל הולצפפל, מודל שריר צייד וכו '. דיאגרמת זרימת המידע PAK-FS מתוארת באיור 4. זה מתחיל מקובץ הקלט ומעבד PAK. כלי הקדם-מעבד PAK מפיק פלט של קובץ DAT, שישמש כקובץ קלט עבור פותר האלמנטים הסופיים. הפלט הסופי מהפותר הוא קבצי VTK הכוללים את התוצאות של סימולציית האלמנט הסופי: מהירויות, לחצים, מתחים ולחצים בחדר השמאלי.
- הורד את קבצי התבנית עבור תנאי הגבול של מהירויות הכניסה והיציאה באמצעות הלחצנים בחלק התחתון (איור 3); אם תנאי גבול זרימה ספציפיים למטופל זמינים, הורד קבצים אלה והשתמש בהם. הורד את אפשרויות הרשת על ידי לחיצה על הכפתור המתאים בחלק התחתון. שמור קבצים אלה בתיקיות פרטיות או ציבוריות.
- העלה קבצים אלה באופן דומה להעלאת תמונות (איור 3). מהירויות הכניסה והיציאה שנקבעו מדמות את מצב התרופה, בעוד שאפשרויות הרשת שולטות בצפיפות רשת היסוד הסופית. להדמיית תנאים ספציפיים למטופל, שנה את ערכי ברירת המחדל של לחץ, זרימה, תכונות החומר ותפקוד הסידן.
- לחץ על הלחצן בצע . זרימת עבודה פועלת חדשה תופיע ברשימה. אם אחד המקטעים של זרימת העבודה אינו ברור, לחץ על לחצן קובץ עזרה (איור 3) כדי להציג הוראות מפורטות לגבי אופן השימוש בזרימת עבודה זו, וכן כדי לפרש את התוצאות.
הערה: אם הכל מתבצע ללא שגיאות, מצב זרימת העבודה ישתנה מ"פועל" ל"סיים בסדר". - הצג את התוצאות עם מספר אפשרויות. לחלופין, הורד את התוצאות; תיקיית התוצאות מכילה קבצי VTK, קבצי CSV והנפשות.
- לחץ על כפתור העין כדי להציג את חלק הפליטה ואת ערכי יעילות העבודה הגלובליים, כמו גם דיאגרמות של לחץ לעומת נפח, לחץ לעומת מאמץ ועבודה שריר הלב לעומת זמן.
- לחץ על כפתור המצלמה כדי להציג בתצוגה מקדימה ולהפעיל אנימציות של שדות התזוזה, הלחץ, הגזירה והמהירות.
- לחץ על כפתור ההדמיה התלת-ממדית כדי לדמיין את הפלטים באופן מקוון ב - ParaView Glance.
- טען מספר קבצי VTK שהורדו בעבר כתוצאה מכך. ראה כמה פרמטרים מעניינים ושנה את השדה, למשל, למהירות להדמיה.
- סובב את הדגם או שנה את ערכת הצבעים. בחרו ' משטח עם קצוות ' או ' מסגרת תיל' לייצוג המשטח. החל את אותה מתודולוגיה על כל קובץ VTK שנטען.
3. קביעת רצף הפעלה חדרית ממדידת א.ק.ג.
הערה: מודל FitzHugh-Nagumo שונה של תא הלב יושם. המוליכים המקדימים עוצבו עם שש האלקטרודות הסטנדרטיות. הפוטנציאל של הלב היה מותאם עם ECG הפוך. החל מההפעלה בצומת הסינואטרלי (שהיא פונקציה של זמן), עם פוטנציאל פעולה הטרוגני דרך הלב ופלג הגוף העליון, המשתמש יכול לקבל פעילות חשמלית במודל פלג הגוף העליון. דיאגרמת זרימת המידע PAK-TORSO ניתנת באיור 5. המשתמש מספק שינוי קנה מידה של מודל פלג הגוף העליון לכל הכיוונים (x,y,z) ופונקציית אות אק"ג. לאחר מכן נוצר המודל המוגדל, והתנהגותו מדומה באמצעות פותר PAK-FS. המשתמש מספק ערכי קלט אלה בקובץ טקסט. הפלט של הסימולציה הוא קובץ VTK עם הפעילות החשמלית של הלב בסביבה משובצת פלג הגוף העליון.
- בדף הבית, עבור אל הפעל זרימת עבודה ולאחר מכן בחר torse-cwl ברשימת זרימות העבודה הזמינות. הוסף הערה או הערה במקטע הערות לפני ביצוע זרימת העבודה.
- לחץ על כפתור קובץ תבנית קלט ושמור את התוכן המוצג בדף האינטרנט כקובץ קלט.txt , שישמש לדגם פלג הגוף העליון. בשדה קובץ קלט , בחר את קובץ הקלט.txt שהורדת. לאחר ייבוא הקובץ, לחץ על לבצע כפתור כדי להתחיל את החישוב.
הערה: זרימת עבודה חדשה תופיע בפינה הימנית התחתונה, כאשר מצבה מוצג כ"פועל"; החישוב נמשך כ 0.5 שעות, ולאחר מכן המצב ישתנה "סיים אישור". - לחץ על כפתורי העין או המצלמה בפינה השמאלית התחתונה כדי לדמיין את דוחות הסימולציה או ההנפשות הזמינים ישירות בפלטפורמה.
הערה: התוצאות כוללות התפלגויות של השדה החשמלי, מהירות, לחץ, לחץ גזירה ועיוות בכל שלב זמן. ניתן לצפות באנימציות עבור כל אחת מההתפלגויות הללו עבור מחזור הלב הכולל. - לחלופין, לחץ על כפתור ההדמיה התלת-ממדית כדי לדמיין את הפלטים באופן מקוון ב - ParaView Glance.
- בחר בלחצן פתח קובץ , עבור אל הכרטיסייה GIRDER , הזן אישורי משתמש אם תתבקש ופתח את התיקיה פרטי.
- בעמוד הבא, בחר את תיקיית פלטי זרימת העבודה ופתח את תיקיית פלג הגוף העליון . פתח את התיקיה הראשונה ברשימה.
- עיין ברשימת קבצי VTK המייצגים את תוצאות הסימולציה. בחר קובץ אחד או יותר ולחץ על בחר כפתור לטעינת הקובץ ב - ParaView Glance.
- טפל בגיאומטריית המודל (כלומר, הזז, סובב, התקרב או התרחק וכו ') באמצעות העכבר.
- בחר אפשרויות שונות לייצוג המודל באופן הבא.
- בחר באפשרות מסגרת תיל כדי לראות את פנים פלג הגוף העליון עם לב משולב בתוך פלג הגוף העליון. בחרו באפשרות ' נקודות' כדי להציג ייצוג מנוקד של דגם פלג הגוף העליון עם רשת לב מלאה.
- התאם את הערך Point Size כדי לשנות את תוצאות התצוגה. התאם את ערך האטימות כדי לראות את פנים פלג הגוף העליון ולהציג תוצאות בתוך רשת הלב.
- לחץ על התפריט הנפתח צבע לפי ובחר את האפשרות הרצויה, למשל, פוטנציאל חשמלי. שנה את סרגל הצבעים של ברירת המחדל לכל אחת מהאפשרויות המפורטות.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
לדוגמה, זרימת העבודה עבור ניתוח תמונה DICOM במצב M אולטרסאונד וחילוץ פרמטרים LV מוצגת באיור 1. מצב M ותצוגה apical ניתן לבדוק בנפרד או אחד אחרי השני, בהתאם לפרמטרים של עניין. אם נבדקים בזה אחר זה, התוצאות מצורפות לקובץ משותף אחד (בנפרד עבור שלבי סיסטולה ודיאסטולה). אם רק תצוגה אחת נבדקת, הערכים עבור פרמטרים לא ידועים נלקחים מקלט הקובץ המוגדר כברירת מחדל .txt (איור 1). הכלי SILICOFCM מנתח את התמונות בהתאם למצב התצוגה, באופן הבא.
ניתוח תצוגה Apical כולל פילוח ה- LV באמצעות מסגרת U-net שהוכשרה בעבר וחישוב המלבן הגובל, שעל בסיסו מחושבים הפרמטרים אורך LV, דיאסטולי, דו-ממדי (LVLd[cm]) ואורך LV, סיסטולי, דו-ממדי (LVLs[cm]). על המשתמש להגדיר אם התצוגה מייצגת את השלב הסיסטולי או הדיאסטולי. קובצי הפלט הם כדלקמן: א) input_parametric_diastole.stl או input_parametric_systole.stl, המייצג את קובץ הקלט עבור פותר PAK עם גיאומטריה שנוצרה, כולל הערך המחושב עבור LVLd[cm] או LVLs[cm] בהתאם לקלט המשתמש של אם התמונה מייצגת דיאסטולה או סיסטולה; ב) output_Apical_view.txt, המכיל את הנתונים עבור LVLd[cm] או LVLs[cm], בהתאם לקלט המשתמש של אם התמונה מייצגת דיאסטולה או סיסטולה.
ניתוח תצוגה במצב M כולל את הגבולות של האזורים האופייניים של LV בתצוגה זו. בהתבסס על אזורים אלה מחושבים פרמטרים שונים. קבצי הפלט הם כדלקמן: א) input_parametric_diastole.stl ו- input_parametric_systole.stl, המייצגים את גיאומטריית הקלט עבור פותר PAK, כולל ערכים מחושבים עבור עובי מחיצה בין-חדרית, דיאסטולית, מצב M (IVSd[cm]), עובי מחיצה בין-חדרית, סיסטולי, מצב M (IVSs[cm]), ממד פנימי LV, דיאסטולי, מצב M (LVIDd[cm]), ממד פנימי LV, סיסטולי, מצב M (LVIDs [cm]), עובי דופן אחורית LV, דיאסטולי, מצב M (LVPWd[cm]), ועובי דופן אחורית LV, סיסטולי, מצב M (LVPWs[cm]); ב) output_M_mode.txt, המכיל את הנתונים עבור IVSd [cm], IVSs [cm], LVIDd [cm], LVIDs [cm], LVPWd [cm] ו- LVPWs [cm].
כדי ליצור דוגמה של LV, מציאותי ככל האפשר, סוג פרמטרי של מודל עם חלקים ספציפיים, כלומר חלק הבסיס, שסתומים (אבי העורקים ומיטרלי), וחלק חיבור (כלומר, החיבור בין בסיס ושסתומים), נוצר באמצעות סימולציות אלמנט סופי PAK (שלב 2) באמצעות הפרמטרים המחושבים של גיאומטריית LV. לכל חלק יש את האורך והרדיוס שלו, כמו גם את מספר השכבות (חלוקה). מודל האלמנטים הסופיים מורכב מתחום זורם המוקף בקיר מוצק המחובר לבסיס ולחלק משכבת החיבור. כדי לדמות את ההתנהגות המציאותית של מודל LV, פונקציות שנקבעו למהירויות במפרצון (מיטרלי; איור 6A) ושקע (אבי העורקים; איור 6B) נעשה שימוש בשסתומים. אלגוריתם לחישוב אוטומטי של כיוון הסיבים הוחל על מודל אלמנט סופי זה. התוצאות עבור שכבה אחת וייצוג מוצק תלת-שכבתי מוצגות באיור 7.
איור 8 מציג את התפלגות הלחץ בתוך המודל הפרמטרי במהלך מחזור זמן של 1 שניות, המחולק ל-10 שלבי זמן. במהלך חמשת השלבים הראשונים, אין שינויי לחץ משמעותיים עד להתכווצות, כאשר הוא מגיע לערך הגבוה ביותר על הסקאלה ואת נפח המודל מתחיל להגדיל באותו זמן. איור 9 מראה כיצד המהירות מתפזרת בתוך החלק הנוזלי של מודל LV. כפי שמוצג, ישנם שיאי ערך בולטים בתוך הענפים, הנגרמים על ידי זרימת נוזלים במהלך מחזור הטעינה/פריקה. איור 10 מראה כיצד התזוזות מתפלגות לאורך המודל. בדומה לשינוי הלחץ, במהלך שני השלבים הראשונים, התזוזות זניחות עד להתכווצות, כאשר הן הופכות למקסימום בחלק התחתון של המודל. לאורך כל הזמן הנותר, המודל חוזר לאט למצבו הלא מעוות. איור 11 מציג את דיאגרמת נפח הלחץ (PV) המתקבלת. איור 12 מראה את הייצוג הווקטורי של מהירויות בתוך החלק הנוזלי של מודל החדרים. וקטורים של מהירות פאזה נמצאים בתוך המסתם המיטרלי במהלך זרימת הנוזל, ולאחר שהחדר מתמלא בנוזל מ-0.7 שניות, הנוזל יוצא דרך החלק האאורטלי.
משוואות מונודומין עבור מערכת פלג הגוף העליון של הלב המצומדות במלואן יושמו על מנת לדמות אלקטרופיזיולוגיה. מוליכות שריר הלב ופלג הגוף העליון מוגדרות בטבלה 1 10,11,12,13,14. שטף אפס שימש עבור Vm כתנאי גבול לכל הגבולות הפנימיים במגע עם הריאות, פלג הגוף העליון והלב. לכן, −n · Γ = 0, כאשר n הוא וקטור נורמלי חיצוני על הגבול ו- Γ הוא וקטור השטף דרך גבול זה עבור המתח התוך-תאי, שווה ל- Γ = − σ· ∂Vm/∂n. לאחר מכן, יושמו גישות קלאסיות לפתרון הבעיה ההופכית של א.ק.ג. באמצעות ניסוח הפוטנציאל האפיקרדיאלי. השיטות בהן נעשה שימוש היו משפחת שיטות טיחונוב ושיטות מבוססות רגולריזציה L 10,11,12,13,14.
פעילות חשמלית של לב שלם בסביבה המשובצת בפלג הגוף העליון מוצגת באיור 13. מערכת הולכה מיוחדת עם מורפולוגיות פוטנציאליות פעולה הטרוגניות ברחבי הלב שולבה בייזום ההפעלה בצומת הסינואטריאלי. איור 14 מראה את המפות של פוטנציאל פני השטח של הגוף אצל נבדק בריא במהלך התקדמות ההפעלה החדרית בתשעה רצפים התואמים את אות האק"ג שנמדד.
איור 1: דיאגרמת זרימת מידע לניתוח תמונה DICOM במצב אולטרסאונד M וחילוץ פרמטרים של LV. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 2: העלאת קבצים חדשים לניתוח. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 3: ממשק משתמש עבור אולטרסאונד אפיקלי תצוגה DICOM ניתוח תמונה וחילוץ פרמטרים LV. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 4: דיאגרמת זרימת מידע PAK-FS. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 5: דיאגרמת זרימת מידע PAK-TORSO. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 6: פרופיל מהירות . (A) פונקציית כניסה של מהירות בחתך מסתם מיטרלי, ו-(B) פונקציית מהירות מוצא בחתך מסתם אבי העורקים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 7: מודל LV פרמטרי עם כיוון סיבים . (A) דופן מוצקה חד-שכבתית, ו-(B) ייצוג קיר מוצק תלת-שכבתי. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 8: שדה לחץ בתוך מודל LV הפרמטרי. חמישה שלבי זמן שונים מוצגים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 9: שדה מהירות בתוך מודל LV פרמטרי. חמישה שלבי זמן שונים מוצגים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 10: תזוזות במודל LV הפרמטרי. חמישה שלבי זמן שונים מוצגים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 11: דיאגרמת לחץ לעומת נפח (PV) עבור מודל אינטראקציה בין נוזל למבנה LV. משך המחזור שנקבע של 1 שניות מחולק ל -10 צעדי זמן. חמישה צעדי זמן מייצגים מוצגים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 12: ייצוג וקטורי של מהירויות במודל הפרמטרי של החדר השמאלי. ארבעה צעדי זמן מייצגים מוצגים. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 13: סימולציה של הפעלת לב שלם בנקודות זמן שונות על אות האק"ג המוביל II. רצפי הפעלה 1-9 ב-(B) תואמים לאות האק"ג ב-(A). פוטנציאל הטרנסממברנה ב- mV מסומן על ידי סרגל הצבע. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 14: מפות של פוטנציאל פני השטח של הגוף בנבדק בריא. התקדמות הפעלת החדר בתשעה רצפים (פאנל תחתון) המתאים לאות א.ק.ג. (פאנל עליון). טווח פעילות הלב ב- mV מסומן על ידי סרגל הצבע. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
| פרמטר | סן | אטריה | א.בן | שלו | BNL | Purkinje | חדרים | |
| A | -0.6 | 0.13 | 0.13 | 0.13 | 0.13 | 0.13 | 0.13 | |
| B | -0.3 | 0 | 0 | 0 | 0 | 0 | 0 | |
| c1(AsV−1 m−3) | 1000 | 2.6 | 2.6 | 2.6 | 2.6 | 2.6 | 2.6 | |
| c2 (AsV−1 m−3) | 1 | 1 | 1 | 1 | 1 | 1 | 1 | |
| D | 0 | 1 | 1 | 1 | 1 | 1 | 1 | |
| e | 0.066 | 0.0132 | 0.0132 | 0.005 | 0.0022 | 0.0047 | 0.006 | |
| א (mV) | 33 | 140 | 140 | 140 | 140 | 140 | 140 | |
| B (mV) | -22 | -85 | -85 | -85 | -85 | -85 | -85 | |
| k | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | |
| σ (mS·m-1) | 0.5 | 8 | 0.5 | 10 | 15 | 35 | 8 | |
טבלה 1: פרמטרים למודל מונודומין עם משוואות FitzHugh-Nahumo שהשתנו.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
פרויקט SILICOFCM הוא פלטפורמה לניסויים קליניים בסיליקו לתכנון אוכלוסיות מטופלים וירטואליות לחיזוי סיכונים, בדיקת ההשפעות של טיפול תרופתי והפחתת ניסויים בבעלי חיים וניסויים קליניים בבני אדם. בדיקת ההשפעות של טיפול תרופתי עוצבה על פי תנאי זרימת גבול כניסה/יציאה שנקבעו, תפקוד סידן ותכונות דופן החומר. פלטפורמה זו משלבת שיטות רב-ממדיות ברמה הסרקומית עם ביצועי לב שלם ורמת אופטימיזציה תפקודית לחיזוי הסיכון של מצבים ספציפיים לחולה במהלך התקדמות מחלת קרדיומיופתיה.
מודלים לב עבור קרדיומיופתיה וצימוד אלקטרומכני של LV בפרויקט SILICOFCM8 הוצג. הגיאומטריה של מודל הלב כללה שבעה אזורים שונים: 1) צומת סינואטרלי; 2) אטריה; 3) צומת אטריובנטריקולרי; 4) צרורו; 5) צרור סיבים; 6) סיבי פורקיניה; ו-7) שריר הלב החדרי. הוצגו מפות פוטנציאליות של פני השטח של הגוף בנושא בריא במהלך התקדמות ההפעלה החדרית בתשעה רצפים המתאימים לאות האק"ג.
לאחר חילוץ הממדים של החדר השמאלי, בהתבסס על אפשרויות רשת, נוצרת הרשת של אלמנטים סופיים והדמיית ניתוח אלמנטים סופיים מופעלת עם מהירויות כניסה ויציאה שנקבעו על ידי המשתמש. משתמשים יכולים לדמיין פתרונות ישירות על הפלטפורמה על ידי עיון באנימציות ובדיאגרמות הזמינות. משתמשים יכולים לדמיין דיאגרמות של נפח לחץ, מתח לחץ וזמן עבודה שריר הלב על-ידי לחיצה על סמל העין. אם משתמשים לוחצים על כפתור המצלמה, מופיעה רשימת האנימציות הזמינות של שדות שונים (תזוזות, לחצים, מהירות, מתחי גזירה).
יש כמה מגבלות למחקר. גיאומטריית LV פרמטרית חולצה מתמונות אולטרסאונד. הגרסה העתידית תיכנס לשחזור גיאומטרי מפורט יותר. קשר ישיר עם תרופות ספציפיות לקרדיומיופתיה אינו מוצג בכתב יד זה. במחקר זה, הוא נשלט באמצעות תנאי גבול לזרימה ולחץ. נתונים גנטיים משולבים כיום באמצעות תפקוד שריר הלב סידן ותכונת החומר הלא ליניארי של הקיר. בגרסה העתידית יילקחו בחשבון פרטים נוספים משחזור תמונה ונתונים גנטיים למחלת קרדיומיופתיה ספציפית לחולה (וריאנט גנטי היפרטרופי ומורחב).
פלטפורמת החישובSILICOFCM 8 תפתח אפיק חדש לניסויים קליניים בסיליקו , במיוחד עבור מחלות לב וחיזוי סיכונים למצב הספציפי למטופל. תקן הזהב לפרקטיקה הקלינית של ימינו לחיזוי סיכונים הוא מחשבון ההישרדות הסטנדרטי לחולי קרדיומיופתיה. הפלטפורמה המשמשת כאן יכולה לתת מידע נוסף בהשוואה לתקן הרפואי הנוכחי, שכן המידול כולל לא רק את הסמנים הביולוגיים אלא גם את הגיאומטריה הספציפית למטופל, תנאי המודינמיקה של זרימה ולחץ, תכונות חומר הקיר (מעיוות בתמונות), ותגובת תרופות עם שילובים שונים של תנאי גבול, תפקוד סידן ותכונות חומר.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
למחברים אין ניגודי עניינים.
Acknowledgments
מחקר זה נתמך על ידי תוכנית המחקר והחדשנות Horizon 2020 של האיחוד האירופי במסגרת הסכם המענקים SILICOFCM 777204 ומשרד החינוך, המדע והפיתוח הטכנולוגי של הרפובליקה של סרביה באמצעות חוזים מס '451-03-68/2022-14/200107. מאמר זה משקף רק את דעות המחברים. הנציבות האירופית אינה אחראית לכל שימוש שעשוי להיעשות במידע שהמאמר מכיל.
Materials
| Name | Company | Catalog Number | Comments |
| SILICOFCM project | www.silicofcm.eu | open access for registered users |
References
- Gibbons Kroeker, C. A., Adeeb, S., Tyberg, J. V., Shrive, N. G. A 2D FE model of the heart demonstrates the role of the pericardium in ventricular deformation. American Journal of Physiology. 291 (5), 2229-2236 (2006).
- Pullan, A. J., Buist, M. L., Cheng, L. K. Mathematically Modelling the Electrical Activity of the Heart - From Cell To Body Surface and Back Again. , World Scientific. (2005).
- Trudel, M. -C., Dub´e, B., Potse, M., Gulrajani, R. M., Leon, L. J. Simulation of QRST integral maps with a membrane based computer heart model employing parallel processing. IEEE Transactions on Biomedical Engineering. 51 (8), 1319-1329 (2004).
- Kojic, M., et al. Smeared multiscale finite element models for mass transport and electrophysiology coupled to muscle mechanics. Frontiers in Bioengineering and Biotechnology. 7, 381 (2019).
- Baillargeon, B., Rebelo, N., Fox, D. D., Taylor, R. L., Kuhl, E. The Living Heart Project: A robust and integrative simulator for human heart function. European Journal of Mechanics - A/Solids. 48, 38-47 (2014).
- Guccione, J. M., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part I--Constitutive relations for fiber stress that describe deactivation. TheJournal of Biomechanical Engineering. 115, 72-81 (1993).
- Guccione, J. M., Waldman, L. K., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part II--Cylindrical models of the systolic left ventricle. The Journal of Biomechanical Engineering. 115, 82-90 (1993).
- H2020 project SILICOFCM: In Silico. trials for drug tracing the effects of sarcomeric protein mutations leading to familial cardiomyopathy. , Available from: www.silicofcm.eu 2018-2022 (2022).
- Kojic, M., et al. Smeared multiscale finite element model for electrophysiology and ionic transport in biological tissue. Computers in Biology and Medicine. 108, 288-304 (2019).
- Wang, Y., Rudy, Y. Application of the method of fundamental solutions to potential-based inverse electrocardiography. Annals of Biomedical Engineering. 34 (8), 1272-1288 (2006).
- Van Oosterom, A. The use of the spatial covariance in computing pericardial potentials. IEEE Transactions on Biomedical Engineering. 46 (7), 778-787 (1999).
- Van Oosterom, A. The spatial covariance used in computing the pericardial potential distribution. Computational Inverse Problems in Electrocardiography. , 1-50 (2001).
- Van Oosterom, A.
- Van Oosterom, A. The equivalent double layer: source models for repolarization. Comprehensive Electrocardiology. , Springer. 227-246 (2010).
