Biomechanics and Modeling in Mechanobiology最新文献

{"title":"Mathematical and numerical tumour development modelling for personalised treatment planning.","authors":"J M Bajkowski, H Piotrzkowska-Wróblewska, B Dyniewicz, C I Bajer","doi":"10.1007/s10237-025-01946-7","DOIUrl":"https://doi.org/10.1007/s10237-025-01946-7","url":null,"abstract":"<p><p>This paper presents a mathematical and numerical framework for modelling and parametrising tumour evolution dynamics to enhance computer-aided diagnosis and personalised treatment. The model comprises six differential equations describing cancer cell and blood vessel concentrations, tissue stiffness, Ki- 67 marker distribution, and the apparent velocity of marker propagation. These equations are coupled through S-functions with adjustable coefficients. An inverse problem approach calibrates the model by fitting adjustable coefficients to patient-specific clinical data, thereby enabling disease progression and treatment response simulations. By integrating historical and prospective patient data supported by machine learning algorithms, this framework holds promise as a robust decision-support tool for optimising therapeutic strategies.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143809986","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stress analysis method for ascending aortic aneurysm based on unloaded geometry with non-uniform thickness distribution.","authors":"Xiaoyu Liu, Zhihong Lin, Shihua Zhao, Fei Li, Qi Gao","doi":"10.1007/s10237-025-01949-4","DOIUrl":"https://doi.org/10.1007/s10237-025-01949-4","url":null,"abstract":"<p><p>Using finite element method (FEM) to compute wall stress is now a common way to assess ascending thoracic aortic aneurysms (ATAA) severity. Medical images can provide aortic geometry for FEM, but thickness information is lacked and the geometry is at loaded state. Therefore, in this study, an unloaded geometry with a non-uniform thickness distribution is reconstructed. Measurements of wall thickness are taken from resected tissue to accurately replicate the thickness distribution. Subsequently, a novel method, derived from the existing fixed-point iterative (FPI) approach, is developed and applied to estimate the unloaded aortic geometry. This new method involves updating the relaxation factor at each iteration to improve robustness by constraining it within a threshold and normalizing it. Compared to the traditional FPI method, this novel approach is better tailored to the aortic geometries examined in this study. The study compares stress results obtained from models with uniform and non-uniform aortic wall thickness, both with and without assuming unloaded conditions. Findings indicate that stress distribution of non-uniform geometry matches better to the measured damage extent. Stress distribution of unloaded geometry is similar to that of loaded geometry, while the use of unloaded geometry enhances the stress gradient. The stress analysis reveals variations across different directions and regions, with the second principal stress (SPS) magnitude that is more sensitive to the circumferential region than the first principal stress (FPS) and von Mises stress (VMS). There is an overlap area between the high SPS region and the most expanded region. The most dilated area usually matched with high SPS region for loaded and unloaded geometry or uniform and non-uniform geometry. Thus, although magnitude of SPS is smaller than that of FPS and of VMS, it is suggested to pay more attention to SPS in severity assessment of ATAA aneurysm.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143802173","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multistep model reduction of coagulation schemes.","authors":"Junyi Chen, Quentin Cazères, Eleonore Riber, Franck Nicoud","doi":"10.1007/s10237-025-01944-9","DOIUrl":"https://doi.org/10.1007/s10237-025-01944-9","url":null,"abstract":"<p><p>This study introduces a comprehensive multistep reduction technique for coagulation models, specifically targeting the dynamics of thrombin generation. By employing a synergistic approach that combines direct relation graph with error propagation, chemical lumping, quasi-steady-state assumption, and conservation analysis, the method efficiently reduces the complexity of original coagulation models without compromising accuracy. Applied to both extrinsic and intrinsic coagulation pathway schemes, this approach significantly diminishes the number of species and reactions, and the resulting reduced schemes appear to be robust to changes in initial conditions relevant to hemophilia A. The findings underscore the potential of this reduction method to facilitate more efficient computational simulations that retain the essential characteristics of different coagulation models.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143802132","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Smooth leaflets with curved belly and attachment edge profiles promote adaptive remodeling in tissue-engineered heart valves: an in silico study.","authors":"Valery L Visser, Sarah E Motta, Simon P Hoerstrup, Frank P T Baaijens, Sandra Loerakker, Maximilian Y Emmert","doi":"10.1007/s10237-025-01937-8","DOIUrl":"https://doi.org/10.1007/s10237-025-01937-8","url":null,"abstract":"<p><p>Tissue-engineered heart valves (TEHVs) are promising valve replacements due to their potential to regenerate into living heart valves, capable of growth and adaptation. Previous TEHVs showed promising results, but often developed progressive leaflet retraction in the long term. In a prior proof-of-concept study, we demonstrated that a novel geometry with more native-like mechanical behavior could give rise to more adaptive remodeling, thereby minimizing leaflet retraction in vivo. In the current study, we aimed to systematically analyze the impact of TEHV geometry on in vivo remodeling under both pulmonary and aortic conditions. Using a bio-inspired in silico framework, we predicted cell-driven, mechano-mediated remodeling in TEHVs post-implantation. Two parameterized valve designs were evaluated under both pulmonary and aortic pressure conditions. The results indicate that a valve design with smooth leaflets, a curved belly profile, and medium to wide attachment edge effectively minimizes stress concentrations and reduces the risk of valve insufficiency in both conditions. Additionally, this design should be tailored to specific hemodynamic conditions to prevent retraction in pulmonary applications and excessive stress concentrations in aortic applications. These insights provide essential guidelines for optimizing TEHV designs, aiming to promote functional remodeling and maintain valve functionality over time, thereby advancing the development of next-generation TEHVs with enhanced long-term outcomes.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143778652","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Role of player-specific white matter parcellation and scaling in impact-induced strain responses.","authors":"Véronique Bouvette, Samuel Guay, Louis De Beaumont, Yvan Petit, Sophie-Andrée Vinet, Eric Wagnac","doi":"10.1007/s10237-025-01945-8","DOIUrl":"https://doi.org/10.1007/s10237-025-01945-8","url":null,"abstract":"<p><p>Head finite element models (hFEMs) are valuable in understanding injury mechanisms in head impacts. Personalizing hFEMs is important for capturing individualized brain responses, with brain volume scaling proving effective. However, the role of refined white matter (WM) parcellation in hFEMs for evaluating brain strain responses, particularly important in the context of subconcussive head impacts (SHIs) often assessed through changes in WM integrity, remains relatively underexplored. This study evaluated the effect of refined subject-specific WM parcellation in 34 WM segments on responses variability due to brain volume variations, using peak maximum principal strain (95MPS) and strain rate (95MPSr) as injury predictive metrics. Data from diffusion-weighted imaging of 21 Canadian varsity football players were utilized to personalize 21 hFEMs. Simulating four different head impacts, representing 50th and 99th percentile resultant accelerations in frontal and angled-top-right directions, refined player-specific WM parcellation better captured variability of strain responses compared to baseline parcellation. Up to 75.71% of 95MPS and 77.14% of 95MPSr responses were deemed different across refined WM segments for players, compared to a maximum of 16.19% of responses with baseline parcellation. These results suggest that player-specific refined WM parcellation improves the ability to capture player-specific responses. Both impact direction and intensity influenced variations in strain response, with angled-top head impacts combined with high intensity showing greater player-specificity compared to lower intensity and frontal head impacts. These findings highlight the potential benefit of model scaling along with player-specific refined WM parcellation in hFEMs for comprehensively evaluating strain responses. Detailed WM parcellation in hFEMs is important for comprehensive injury assessment, enhancing the alignment of hFEMs with imaging studies evaluating changes in WM integrity across segments. The simple and straightforward method presented herein to achieve player-specific strain response is promising for future SHI studies.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143770830","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Computational model of coarctation of the aorta in rabbits suggests persistent ascending aortic remodeling post-correction","authors":"Ashley A. Hiebing,&nbsp;Matthew A. Culver,&nbsp;John F. LaDisa Jr.,&nbsp;Colleen M. Witzenburg","doi":"10.1007/s10237-025-01933-y","DOIUrl":"10.1007/s10237-025-01933-y","url":null,"abstract":"<div><p>Coarctation of the aorta (CoA) is a common congenital cardiovascular lesion that presents as a localized narrowing of the proximal descending aorta. While improvements in surgical and catheter-based techniques have increased short-term survival, there is a high long-term risk of hypertension and a reduced average lifespan despite correction. Computational models can be used to estimate aortic remodeling and peripheral vascular compensation, potentially serving as key tools in developing a mechanistic understanding of the interplay between pre-treatment dynamics, post-treatment recovery, and long-term hypertension risk. In this study, we developed a lumped-parameter model of the heart and circulation to simulate CoA. After fitting model parameters using imaging and catheterization data from healthy rabbits, we then used the model to estimate differences in ascending aortic compliance and peripheral resistance between the healthy group and rabbits with both untreated and corrected CoA using their imaging and catheterization data. CoA was defined by the current putative clinical treatment threshold (a pressure gradient &gt; 20 mm Hg). Model inputs were fitted such that outputs matched reported stroke volume, ejection fraction, systolic and diastolic aortic pressure, peak aortic flow, mean and peak blood pressure gradients, and upper-to-lower body flow split, with all results falling within one standard deviation of the data for all groups. In the untreated CoA and corrected simulations, a decrease in ascending aortic compliance was necessary to match reported hemodynamics. This suggests exposure to a pressure gradient &gt; 20 mm Hg results in vascular remodeling that persists after repair, a process strongly correlated with hypertension.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"24 2","pages":"683 - 700"},"PeriodicalIF":3.0,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143668568","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical simulation of voluntary respiration in a model of the whole human lower airway","authors":"Xinying Ou,&nbsp;Jiahuan Meng,&nbsp;Chen Ma,&nbsp;Huajing Wan,&nbsp;Yu Chen,&nbsp;Fengming Luo","doi":"10.1007/s10237-025-01932-z","DOIUrl":"10.1007/s10237-025-01932-z","url":null,"abstract":"<div><p>The lung model construction is limited to the local scale, and the numerical simulation of autonomous breathing is mostly computed from top to bottom in recent research. In this study, models of the entire lower airway from G0–G23 were constructed, and computational simulations were performed for the alveolar model using coupled fluid–solid analysis with pressure changes on the wall and for the rigid bronchial model using computational fluid dynamics by transmitting the boundary conditions step from bottom to top. This paper provides the results under spontaneous respiration, including the ventilation volume of the tracheobronchial tree, the situation of the internal flow field, and the mechanical characteristics of the lung tissues. The mechanical characteristics and the lung functions computed by the models were consistent with clinical or experimental data. This model could provide quantitative analysis results of respiratory mechanics in the lower respiratory tract of the human, which offers a reference for mechanical studies, such as the morphological changes and differentiation of cell types induced by force stimulation and tumor induction. Furthermore, various pathological models can be developed based on this model.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"24 2","pages":"663 - 681"},"PeriodicalIF":3.0,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143646768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Advances in computational modeling of cytokine and growth factor dynamics in bone healing: a scoping review.","authors":"Ahmad Hedayatzadeh Razavi, Nazanin Nafisi, Maria Velasquez-Hammerle, Mohammad Javad Shariyate, Mohammad Khak, Alireza Mirahmadi, Megan McNichol, Edward K Rodrogiuez, Ara Nazarian","doi":"10.1007/s10237-025-01938-7","DOIUrl":"https://doi.org/10.1007/s10237-025-01938-7","url":null,"abstract":"<p><p>Bone healing is a complex process regulated by intricate biological and mechanical factors and spatially varied regions over time. This scoping review synthesizes current computational models that incorporate cytokines and growth factors, examining their role in bone healing. Through a systematic analysis of 71 studies, this review identifies and categorizes the modeling methodologies used, including mathematical, finite element, agent-based, mechanobiological, pharmacobiological, and hybrid approaches. The findings highlight the predominant use of mathematical models while noting a recent shift toward more sophisticated techniques like finite element and agent-based models. Key cytokines and growth factors, such as TGF-β, RANK-RANKL-OPG, and PTH, are repeatedly used, underscoring their essential roles in regulating cellular processes. This review also analyzes parameter selection and validation strategies, identifying gaps in current practices and emphasizing the need for high-quality experimental validation to improve model reliability. Some bibliometric analyses provide insights into citation networks and keyword co-occurrence, illustrating influential studies in the field and central themes. The findings offer a foundation for future research to enhance model accuracy, aiming toward more predictive and clinically relevant models accounting for biology and mechanics in bone healing.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143630137","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The order of precedence in treatment of multiple intracranial aneurysms: insights from a fluid–structure interaction study","authors":"Kornelia M. Kliś,&nbsp;Jerzy Gąsowski,&nbsp;Antoni Cierniak,&nbsp;Borys M. Kwinta,&nbsp;Krzysztof Stachura,&nbsp;Tadeusz J. Popiela,&nbsp;Igor Szydłowski,&nbsp;Bartłomiej Łasocha,&nbsp;Karolina Piotrowicz,&nbsp;Tomasz Grodzicki,&nbsp;Roger M. Krzyżewski","doi":"10.1007/s10237-025-01928-9","DOIUrl":"10.1007/s10237-025-01928-9","url":null,"abstract":"<div><p>Treatment strategy for multiple intracranial aneurysms is challenging, as in many cases the choice of the order in which to treat aneurysms is not based on high-quality evidence. We aimed to digitally simulate clinical scenarios of two different orders in which multiple intracranial aneurysms were treated and analyze changes in hemodynamics after first stage of treatment. We prospectively included patients with two intracranial aneurysms, with order of treatment difficult to determine based on clinical data alone. For each patient we prepared three models of arteries harboring aneurysms: with both aneurysms present and with one of them removed. Computational modeling of blood flow using fluid–structure interaction methodology was performed for each model. Hemodynamic parameters of aneurysm domes were compared between models with both aneurysm present, and models in with aneurysms were removed in changing order. In 25 included patients, the calculated hemodynamic parameters such as Time-Averaged Wall Shear Stress (0.46 ± 0.40 vs. 0.54 ± 0.44 Pa; <i>p</i> &lt; 0.01) and surface vortex fraction (12.73% ± 7.92% vs. 14.26% ± 7.46%; <i>p</i> = 0.02) decreased after first stage of treatment, while Time-Averaged Wall Shear Stress Gradient (1.44 ± 0.41 vs. 1.34 ± 0.46 Pa; <i>p</i> = 0.04) and percentage of wall shear stress &lt; 0.5 Pa (50.13% ± 33.01% vs. 44.08% ± 34.16%; <i>p</i> &lt; 0.01) increased. Changes of wall shear stress in remaining aneurysm dome were independently correlated with dome-to-neck ratio of both removed and remaining aneurysms. Hemodynamics of untreated aneurysm worsens after first stage of treatment. Dome-to-neck ratio of both treated and untreated aneurysm was the strongest and independent predictor of that change.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"24 2","pages":"589 - 598"},"PeriodicalIF":3.0,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143630140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A combined 4D flow MR imaging and fluid-structure interaction analysis of ascending thoracic aortic aneurysms.","authors":"Yu Zhu, Chlöe Armour, Binghuan Li, Selene Pirola, Yousuf Salmasi, Thanos Athanasiou, Declan P O'Regan, Xiao Yun Xu","doi":"10.1007/s10237-025-01939-6","DOIUrl":"https://doi.org/10.1007/s10237-025-01939-6","url":null,"abstract":"<p><p>This study aimed to characterize the altered hemodynamics and wall mechanics in ascending thoracic aortic aneurysms (ATAA) by employing fully coupled two-way fluid-structure interaction (FSI) analyses. Our FSI models incorporated hyperelastic wall mechanical properties, prestress, and patient-specific inlet velocity profiles (IVP) extracted from 4D flow magnetic resonance imaging (MRI). By performing FSI analyses on 7 patient-specific ATAA models and 6 healthy aortas, the primary objective of the study was to compare hemodynamic and biomechanical features in ATAA versus healthy controls. A secondary objective was to examine the need for 4D flow MRI-derived IVP in FSI simulations by comparing results with those using two commonly adopted idealized IVPs: Flat-IVP and Para-IVP for selected cases. Our results show that, compared to the healthy aortas, the ATAA models exhibited highly disturbed blood flow in the ascending aorta. Consequently, maximum turbulent kinetic energy (TKE) at peak systole (155.0 ± 188.4 Pa) and maximum time-averaged wall shear stress (TAWSS) (8.6 ± 6.5 Pa) were significantly higher in the ATAA cohort, compared to 0.6 ± 0.5 Pa and 2.8 ± 0.7 Pa in the healthy aortas. Peak wall stress was also nearly doubled in the ATAA group (414 ± 108 kPa vs. 215 ± 31 kPa). Additionally, comparisons of simulation results across models with different IVPs underscore the importance of prescribing 3D-IVP at the inlet, especially for ATAA cases. Using idealized IVPs in two selected ATAA models (P1 and P7) substantially reduced the maximum TKE from 571 Pa to 0.01 Pa (Flat-IVP) and 0.02 Pa (Para-IVP) in P1 and from 73 Pa to 0.01 Pa (Flat-IVP) and 0.08 Pa (Para-IVP) in P7, while the maximum TAWSS in the ascending aorta decreased from 9.6 Pa to 0.7 Pa (Flat-IVP) and 0.9 Pa (Para-IVP) in P1, and from 3.6 Pa to 1.2 Pa and 0.9 Pa, respectively, in P7. Moreover, idealized IVPs also caused the peak wall stress to reduce by up to 11.5% in P1 with severe aortic valve stenosis, and by up to 2% in P7 with mild aortic regurgitation. These results highlight the importance of FSI simulations combined with 4D flow MRI in capturing realistic hemodynamic and biomechanical changes in aneurysmal aortas.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143603215","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}