{"title":"Timing of resting zone parathyroid hormone-related protein expression affects maintenance of the growth plate during secondary ossification: a computational study.","authors":"Jorik Stoop, Yuka Yokoyama, Taiji Adachi","doi":"10.1007/s10237-024-01899-3","DOIUrl":"https://doi.org/10.1007/s10237-024-01899-3","url":null,"abstract":"<p><p>Secondary ossification and maintenance of the growth plate are crucial aspects of long bone formation. Parathyroid hormone-related protein (PTHrP) has been implicated as a key factor in maintaining the growth plate, and studies suggest that PTHrP expression in the resting zone is closely related with formation of the secondary ossification center (SOC). However, details of the relationship between resting zone PTHrP expression and preservation of the growth plate remain unclear. In this study, we aim to investigate the role of resting zone PTHrP expression on maintenance of the growth plate using a computational method. We extend an existing continuum-based particle model of tissue morphogenesis to include PTHrP and Indian hedgehog (Ihh) signaling, allowing the model to capture biochemical and mechanical regulation of individual cell activities. Our model indicates that the timing of resting zone PTHrP expression-specifically the rate of increase in production at the onset of SOC formation-is potentially a crucial mechanism for maintenance of the growth plate.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142643457","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 non-intrusive reduced-order model for finite element analysis of implant positioning in total hip replacements.","authors":"Marlis Reiber, Fynn Bensel, Zhibao Zheng, Udo Nackenhorst","doi":"10.1007/s10237-024-01903-w","DOIUrl":"https://doi.org/10.1007/s10237-024-01903-w","url":null,"abstract":"<p><p>Sophisticated high-fidelity simulations can predict bone mass density (BMD) changes around a hip implant after implantation. However, these models currently have high computational demands, rendering them impractical for clinical settings. Model order reduction techniques offer a remedy by enabling fast evaluations. In this work, a non-intrusive reduced-order model, combining proper orthogonal decomposition with radial basis function interpolation (POD-RBF), is established to predict BMD distributions for varying implant positions. A parameterised finite element mesh is morphed using Laplace's equation, which eliminates tedious remeshing and projection of the BMD results on a common mesh in the offline stage. In the online stage, the surrogate model can predict BMD distributions for new implant positions and the results are visualised on the parameterised reference mesh. The computational time for evaluating the final BMD distribution around a new implant position is reduced from minutes to milliseconds by the surrogate model compared to the high-fidelity model. The snapshot data, the surrogate model parameters and the accuracy of the surrogate model are analysed. The presented non-intrusive surrogate model paves the way for on-the-fly evaluations in clinical practice, offering a promising tool for planning and monitoring of total hip replacements.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142611862","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":"Comparison and identification of human coronary plaques with/without erosion using patient-specific optical coherence tomography-based fluid-structure interaction models: a pilot study.","authors":"Yanwen Zhu, Chen Zhao, Zheyang Wu, Akiko Maehara, Dalin Tang, Liang Wang, Zhanqun Gao, Yishuo Xu, Rui Lv, Mengde Huang, Xiaoguo Zhang, Jian Zhu, Haibo Jia, Bo Yu, Minglong Chen, Gary S Mintz","doi":"10.1007/s10237-024-01906-7","DOIUrl":"https://doi.org/10.1007/s10237-024-01906-7","url":null,"abstract":"<p><p>Plaque erosion (PE) with secondary thrombosis is one of the key mechanisms of acute coronary syndrome (ACS) which often leads to drastic cardiovascular events. Identification and prediction of PE are of fundamental significance for disease diagnosis, prevention and treatment. In vivo optical coherence tomography (OCT) data of eight eroded plaques and eight non-eroded plaques were acquired to construct three-dimensional fluid-structure interaction models and obtain plaque biomechanical conditions for investigation. Plaque stenosis severity, plaque burden, plaque wall stress (PWS) and strain (PWSn), flow shear stress (FSS), and ΔFSS (FSS variation in time) were extracted for comparison and prediction. A logistic regression model was used to predict plaque erosion. Our results indicated that the combination of mean PWS and mean ΔFSS gave best prediction (AUC = 0.866, 90% confidence interval (0.717, 1.0)). The best single predictor was max ΔFSS (AUC = 0.819, 90% confidence interval (0.624, 1.0)). The average of maximum FSS values from eroded plaques was 76% higher than that from the non-eroded plaques (127.96 vs. 72.69 dyn/cm<sup>2</sup>) while the average of mean FSS from erosion sites of the eight eroded plaques was 48.6% higher than that from sites without erosion (71.52 vs. 48.11 dyn/cm<sup>2</sup>). The average of mean PWS from plaques with erosion was 22.83% lower than that for plaques without erosion (83.2 kPa vs. 107.8 kPa). This pilot study suggested that combining plaque stress, strain and flow shear stress could help better identify patients with potential plaque erosion, enabling possible early intervention therapy. Further studies are needed to validate our findings.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142611864","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":"Investigating wall shear stress and the static pressure in bone scaffolds: a study of porosity and fluid flow dynamics.","authors":"Vedang Gadgil, Shriram Kumbhojkar, Tushar Sapre, Prathamesh Deshmukh, Pankaj Dhatrak","doi":"10.1007/s10237-024-01904-9","DOIUrl":"https://doi.org/10.1007/s10237-024-01904-9","url":null,"abstract":"<p><p>In bone tissue engineering, scaffolds are crucial as they provide a suitable structure for cell proliferation. Transporting Dulbecco's Modified Eagle Medium (DMEM) to the cells and regulating the scaffold's biocompatibility are both controlled by the dynamics of the fluid passing through the scaffold pores. Scaffold design selection and modeling are thus important in tissue engineering to achieve successful bone regeneration. This study aims to design and analyze three scaffold designs-Face-Centered Cubic (FCC), and two newly developed designs Octagonal Truss and a Square Pyramid with four porosity variations. The research aims to analyze the effect of design and porosity variation on pressure and wall shear stress, essential for analyzing scaffold biocompatibility in tissue engineering. Three scaffold designs with varying porosities with strut diameters ranging from 0.3 to 0.6 mm were modeled to analyze the behavior using BioMed Clear Resin. The fluid dynamics within these scaffolds were then examined using Computational Fluid Dynamics (CFD) to understand how different porosity levels affect fluid flow pressure and wall shear stress. The findings revealed variations in wall shear stress and their influence on cell proliferation. The maximum value of wall shear stress (WSS) is observed in the Square Pyramid model. The analysis shows that WSS at the inlet decreases as strut diameters increase or porosity percentages rise offering valuable insights for the development of effective scaffold designs. It can be concluded from the results that the Square Pyramid design has the highest value of WSS, thus increasing the chances of cell growth. From a biological perspective, the results of this work show promise for creating better scaffolds for tissue engineering.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142542750","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}
Keshav A Kailash, Shamimur R Akanda, Alexandra L Davis, Christie L Crandall, Luis A Castro, Lori A Setton, Jessica E Wagenseil
{"title":"A multiphasic model for determination of mouse ascending thoracic aorta mass transport properties with and without aneurysm.","authors":"Keshav A Kailash, Shamimur R Akanda, Alexandra L Davis, Christie L Crandall, Luis A Castro, Lori A Setton, Jessica E Wagenseil","doi":"10.1007/s10237-024-01897-5","DOIUrl":"https://doi.org/10.1007/s10237-024-01897-5","url":null,"abstract":"<p><p>Thoracic aortic aneurysms (TAAs) are associated with aortic wall remodeling that affects transmural transport or the movement of fluid and solute across the wall. In previous work, we used a Fbln4<sup>E57K/E57K</sup> (MU) mouse model to investigate transmural transport changes as a function of aneurysm severity. We compared wild-type (WT), MU with no aneurysm (MU-NA), MU with aneurysm (MU-A), and MU with an additional genetic mutation that led to increased aneurysm penetrance (MU-XA). We found that all aneurysmal aortas (MU-A and MU-XA) had lower fluid flux compared to WT. Non-aneurysmal aortas (MU-NA) had higher 4 kDa FITC-dextran solute flux than WT, but aneurysmal MU-A and MU-XA aortas had solute fluxes similar to WT. Our experimental results could not isolate competing factors, such as changes in aortic geometry and solid material properties among these mouse models, to determine how intrinsic transport properties change with aneurysm severity. The objective of this study is to use biphasic and multiphasic models to identify changes in transport material properties. Our biphasic model indicates that hydraulic permeability is significantly decreased in the severe aneurysm model (MU-XA) compared to non-aneurysmal aortas (MU-NA). Our multiphasic model shows that effective solute diffusivity is increased in MU-NA aortas compared to all others. Our findings reveal changes in intrinsic transport properties that depend on aneurysm severity and are important for understanding the movement of fluids and solutes that may play a role in the diagnosis, progression, or treatment of TAA.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142542646","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":"Piezoelectricity and flexoelectricity in biological cells: the role of cell structure and organelles.","authors":"Akepogu Venkateshwarlu, Akshayveer, Sundeep Singh, Roderick Melnik","doi":"10.1007/s10237-024-01895-7","DOIUrl":"https://doi.org/10.1007/s10237-024-01895-7","url":null,"abstract":"<p><p>Living tissues experience various external forces on cells, influencing their behaviour, physiology, shape, gene expression, and destiny through interactions with their environment. Despite much research done in this area, challenges remain in our better understanding of the behaviour of the cell in response to external stimuli, including the arrangement, quantity, and shape of organelles within the cell. This study explores the electromechanical behaviour of biological cells, including organelles like microtubules, mitochondria, nuclei, and cell membranes. A two-dimensional bio-electromechanical model for two distinct cell structures has been developed to analyze the behavior of the biological cell to the external electrical and mechanical responses. The piezoelectric and flexoelectric effects have been included via multiphysics coupling for the biological cell. All the governing equations have been discretized and solved by the finite element method. It is found that the longitudinal stress is absent and only the transverse stress plays a crucial role when the mechanical load is imposed on the top side of the cell through compressive displacement. The impact of flexoelectricity is elucidated by introducing a new parameter called the maximum electric potential ratio ( <math><msub><mi>V</mi> <mrow><mi>R</mi> <mo>,</mo> <mtext>max</mtext></mrow> </msub> </math> ). It has been found that <math><msub><mi>V</mi> <mrow><mi>R</mi> <mo>,</mo> <mtext>max</mtext></mrow> </msub> </math> depends upon the orientation angle and shape of the microtubules. The magnitude of <math><msub><mi>V</mi> <mrow><mi>R</mi> <mo>,</mo> <mtext>max</mtext></mrow> </msub> </math> exhibit huge change when we change the shape and orientation of the organelles, which in some cases (boundary condition (BC)-3) can reach to three times of regular shape organelles. Further, the study reveals that the number of microtubules significantly impacts effective elastic and piezoelectric coefficients, affecting cell behavior based on structure, microtubule orientation, and mechanical stress direction. The insight obtained from the current study can assist in advancements in medical therapies such as tissue engineering and regenerative medicine.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-10-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142492502","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}
Di Zuo, Mingji Zhu, Daye Chen, Qiwen Xue, Stéphane Avril, Klaus Hackl, Yiqian He
{"title":"Three-dimensional anisotropic unified continuum model for simulating the healing of damaged soft biological tissues","authors":"Di Zuo, Mingji Zhu, Daye Chen, Qiwen Xue, Stéphane Avril, Klaus Hackl, Yiqian He","doi":"10.1007/s10237-024-01888-6","DOIUrl":"10.1007/s10237-024-01888-6","url":null,"abstract":"<div><p>The soft biological tissues have the ability to heal and self-repair after damage or injury. During the healing process, damaged tissues are replaced by newly produced undamaged tissue to restore homeostasis. Computational modeling serves as an effective tool for simulating the healing process and understanding the underlying mechanisms. In previous work, we developed the first unified continuum damage model for the healing of soft biological tissues. However, the initial theory lacked generalizability to more realistic scenarios and applicability to biomechanical problems due to the simplicity of the isotropic constitutive model and two-dimensional simulations. Therefore, we further improve our approach by developing a three-dimensional anisotropic unified healing model to address more realistic challenges. By using the Holzapfel–Gasser–Ogden model as the hyperelastic term, the influence of the collagen fibers is considered and the reorientation of fibers in healing is simulated. Three numerical examples related to hypertension, aneurysm, and restenosis of the atherosclerotic artery after balloon angioplasty are presented to demonstrate the effectiveness of the proposed model. By comparing numerical solutions and reference solutions, we demonstrate the ability of the proposed model in simulating long-term tissue healing process and analyze the impact of anisotropic terms.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2193 - 2212"},"PeriodicalIF":3.0,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142455155","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}
Daniel Paukner, Jay D. Humphrey, Christian J. Cyron
{"title":"Multiscale homogenized constrained mixture model of the bio-chemo-mechanics of soft tissue growth and remodeling","authors":"Daniel Paukner, Jay D. Humphrey, Christian J. Cyron","doi":"10.1007/s10237-024-01884-w","DOIUrl":"10.1007/s10237-024-01884-w","url":null,"abstract":"<div><p>Constrained mixture models have successfully simulated many cases of growth and remodeling in soft biological tissues. So far, extensions of these models have been proposed to include either intracellular signaling or chemo-mechanical coupling on the organ-scale. However, no version of constrained mixture models currently exists that includes both aspects. Here, we propose such a version that resolves cellular signal processing by a set of logic-gated ordinary differential equations and captures chemo-mechanical interactions between cells by coupling a reaction-diffusion equation with the equations of nonlinear continuum mechanics. To demonstrate the potential of the model, we present 2 case studies within vascular solid mechanics: (i) the influence of angiotensin II on aortic growth and remodeling and (ii) the effect of communication between endothelial and intramural arterial cells via nitric oxide and endothelin-1.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2115 - 2136"},"PeriodicalIF":3.0,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01884-w.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142455154","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Cellular solids and prestressed affine networks as models of the elastic behavior of soft biological structures.","authors":"Dimitrije Stamenović","doi":"10.1007/s10237-024-01894-8","DOIUrl":"https://doi.org/10.1007/s10237-024-01894-8","url":null,"abstract":"<p><p>We reviewed two microstructural models, cellular solid models and prestressed affine network models, that have been used previously in studies of elastic behavior of soft biological materials. These models provide simple and mathematically transparent equations that can be used to interpret experimental data and to obtain quantitative predictions of the elastic properties of biological structures. In both models, volumetric density and elastic properties of the microstructure are key determinants of the macroscopic elastic properties. In the prestressed network model, geometrical rearrangement of the microstructure (kinematic stiffness) is also important. As examples of application of these models, we considered the shear behavior of the cytoskeleton of adherent cells, of the collagen network of articular cartilage, and of the lung parenchymal network since their ability to resist shear is important for their normal biological and physiological functions. All three networks carry a pre-existing stress (prestress). We predicted their shear moduli using the microstructural models and compared those predictions with existing experimental data. Prestressed network models of the cytoskeleton and of the lung parenchyma provided a better correspondence to experimental data than cellular solid models. Both cellular solid and prestressed network models of the cartilage collagen network provided reasonable agreements with experimental values. These findings suggested that the kinematic stiffness and material stiffness of microstructural elements were both important determinants of the shear modulus of the cytoskeleton and of the lung parenchyma, whereas elasticity of collagen fibrils had a predominant role in the cartilage shear behavior.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142455153","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":"Towards a framework for predicting immunotherapy outcome: a hybrid multiscale mathematical model of immune response to vascular tumor growth","authors":"Sayyed Mohammad Ali Mortazavi, Bahar Firoozabadi","doi":"10.1007/s10237-024-01891-x","DOIUrl":"10.1007/s10237-024-01891-x","url":null,"abstract":"<div><p>Studying tumor immune microenvironment (TIME) is pivotal to understand the mechanism and predict the outcome of cancer immunotherapy. Systems biology mathematical models can consider and control various factors of TIME and therefore explore the anti-tumor immune response meticulously. However, the role of tumor vasculature in the recruitment of T cells and the mechanism of T cell migration through TIME have not been studied comprehensively. In this work, we developed a hybrid discrete-continuum multi-scale model to study TIME. The mathematical model includes angiogenesis and T cell recruitment via tumor vasculature. Moreover, solid tumor growth, vascular growth and remodeling, interstitial fluid flow, hemodynamics, and blood rheology are all considered in the model. In addition, different aspects of T cells, including their migration, proliferation, subtype conversion, and interaction with tumor cells are thoroughly included. The model reproduces spatiotemporal distribution of tumor infiltrating T cells that mimics histopathological patterns. Furthermore, TIME model robustly recapitulates different phases of tumor immunoediting. We also examined a number of biomarkers to predict the outcome of immune checkpoint blockade (ICB) treatment. The results demonstrated that although tumor mutational burden (TMB) may predict non-responders to ICB, a combination of different biomarkers is essential to predict the majority of the responders. Based on our results, the ICB response rate varies significantly from 28 to 89% depending on the values of different parameters, even in the cases with high TMB.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2243 - 2264"},"PeriodicalIF":3.0,"publicationDate":"2024-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142379812","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}