Biomechanics and Modeling in Mechanobiology最新文献

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier–Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

Architectural parameters of skeletal muscle such as pennation angle provide valuable information on muscle function, since they can be related to the muscle force generating capacity, fiber packing, and contraction velocity. In this paper, we introduce a 3D ultrasound-based workflow for determining 3D fascicle orientations of skeletal muscles. We used a custom-designed automated motor driven 3D ultrasound scanning system for obtaining 3D ultrasound images. From these, we applied a custom-developed multiscale-vessel enhancement filter-based fascicle detection algorithm and determined muscle volume and pennation angle. We conducted trials on a phantom and on the human tibialis anterior (TA) muscle of 10 healthy subjects in plantarflexion (157 ± 7(^circ)), neutral position (109 ± 7(^circ), corresponding to neutral standing), and one resting position in between (145 ± 6(^circ)). The results of the phantom trials showed a high accuracy with a mean absolute error of 0.92 ± 0.59(^circ). TA pennation angles were significantly different between all positions for the deep muscle compartment; for the superficial compartment, angles are significantly increased for neutral position compared to plantarflexion and resting position. Pennation angles were also significantly different between superficial and deep compartment. The results of constant muscle volumes across the 3 ankle joint angles indicate the suitability of the method for capturing 3D muscle geometry. Absolute pennation angles in our study were slightly lower than recent literature. Decreased pennation angles during plantarflexion are consistent with previous studies. The presented method demonstrates the possibility of determining 3D fascicle orientations of the TA muscle in vivo.

We have developed a finite element model to simulate the penetration of nanoneedles into the cellular nucleus. It is found that the nuclear lamina, the primary supporting structure of the nuclear membrane, plays a crucial role in maintaining the integrity of the nuclear envelope and enhancing stress concentration in the nuclear membrane. Notably, nuclear lamina A exhibits a more pronounced effect compared to nuclear lamina B. Subsequently, we further conducted experiments by controlling the time of osteopontin (OPN) treatment to modify the nuclear lamina density, and the results showed that an increase in nuclear lamina density enhances the probability of nanoneedle penetration into the nuclear membrane. Through employing both simulation and experimental techniques, we have gathered compelling evidence indicating that an augmented density of nuclear lamina A can enhance the penetration of nanoneedles into the nuclear membrane.

Vascular tone regulation is a crucial aspect of cardiovascular physiology, with significant implications for overall cardiovascular health. However, the precise physiological mechanisms governing smooth muscle cell contraction and relaxation remain uncertain. The complexity of vascular tone regulation stems from its multiscale and multifactorial nature, involving global hemodynamics, local flow conditions, tissue mechanics, and biochemical pathways. Bridging this knowledge gap and translating it into clinical practice presents a challenge. In this paper, a computational model is presented to integrate chemo-mechano-biological pathways with cardiovascular biomechanics, aiming to unravel the intricacies of vascular tone regulation. The computational framework combines an algebraic description of global hemodynamics with detailed finite element analyses at the scale of vascular segments for describing their passive and active mechanical response, as well as the molecular transport problem linked with chemo-biological pathways triggered by wall shear stresses. Their coupling is accounted for by considering a two-way interaction. Specifically, the focus is on the role of nitric oxide-related molecular pathways, which play a critical role in modulating smooth muscle contraction and relaxation to maintain vascular tone. The computational framework is employed to examine the interplay between localized alterations in the biomechanical response of a specific vessel segment—such as those induced by calcifications or endothelial dysfunction–and the broader global hemodynamic conditions—both under basal and altered states. The proposed approach aims to advance our understanding of vascular tone regulation and its impact on cardiovascular health. By incorporating chemo-mechano-biological mechanisms into in silico models, this study allows us to investigate cardiovascular responses to multifactorial stimuli and incorporate the role of adaptive homeostasis in computational biomechanics frameworks.

We present a class of model-free Data-Driven solvers that effectively enable the utilization of in situ and in vivo imaging data directly in full-scale calculations of the mechanical response of the human brain to sonic and ultrasonic stimulation, entirely bypassing the need for analytical modeling or regression of the data. The well-posedness of the approach and its convergence with respect to data are proven analytically. We demonstrate the approach, including its ability to make detailed spatially resolved patient-specific predictions of wave patterns, using public-domain MRI images, MRE data and commercially available solid-mechanics software.

Cellular traction forces are contractile forces that depend on the material/substrate stiffness and play essential roles in sensing mechanical environments and regulating cell morphology and function. Traction forces are primarily generated by the actin cytoskeleton and transmitted to the substrate through focal adhesions. The cell nucleus is also believed to be involved in the regulation of this type of force; however, the role of the nucleus in cellular traction forces remains unclear. In this study, we explored the effects of nucleus-actin filament coupling on cellular traction forces in human dermal fibroblasts cultured on substrates with varying stiffness (5, 15, and 48 kPa). To investigate these effects, we transfected the cells with a dominant-negative Klarsicht/ANC-1/Syne homology (DN-KASH) protein that was designed to displace endogenous linker proteins and disrupt nucleus-actin cytoskeleton connections. The force that exists between the cytoskeleton and the nucleus (nuclear tension) was also evaluated with a fluorescence resonance energy transfer (FRET)-based tension sensor. We observed a biphasic change in cellular traction forces with a peak at 15 kPa, regardless of DN-KASH expression, that was inversely correlated with the nuclear tension. In addition, the relative magnitude and distribution of traction forces in nontreated wild-type cells were similar across different stiffness conditions, while DN-KASH-transfected cells exhibited a different distribution pattern that was impacted by the substrate stiffness. These results suggest that the nucleus-actin filament coupling play a homeostatic role by maintaining the relative magnitude of cellular traction forces in fibroblasts under different stiffness conditions.

Total temporomandibular joint replacement (TMJR) surgery is the established treatment for severe temporomandibular joint disorders. While TMJR surgery is known to increase mouth-opening capacity, reduce pain and improve quality of life, little is known about post-surgical jaw function during activities of daily living such as biting and chewing. The aim of this study was to use subject-specific 3D bite force measurements to evaluate the magnitude and direction of joint loading in unilateral total TMJR patients and compare these data to those in healthy control subjects. An optoelectronic tracking system was used to measure jaw kinematics while biting a rubber sample for 5 unilateral total TMJR patients and 8 controls. Finite element simulations driven by the measured kinematics were employed to calculate the resultant bite force generated when compressing the rubber between teeth during biting tasks. Subject-specific musculoskeletal models were subsequently used to calculate muscle and TMJ loading. Unilateral total TMJR patients generated a bite force of 249.6 ± 24.4 N and 164.2 ± 62.3 N when biting on the contralateral and ipsilateral molars, respectively. In contrast, controls generated a bite force of 317.1 ± 206.6 N. Unilateral total TMJR patients biting on the contralateral molars had a significantly higher lateral TMJ force direction (median difference: 63.6°, p = 0.028) and a significantly lower ratio of working TMJ force to bite force (median difference: 0.17, p = 0.049) than controls. Results of this study may guide TMJ prosthesis design and evaluation of dental implants.

The present study investigates the multiphasic nature of the mechanical behavior of human dermis. Motivated by experimental observations and by consideration of its composition, a quadriphasic model of the dermis is proposed, distinguishing solid matrix components, interstitial fluid and charged constituents moving within the fluid, i.e., anions and cations. Compression and tensile experiments with and without change of osmolarity of the bath are performed to characterize the chemo-mechanical coupling in the dermis. Model parameters are determined through inverse analysis. The computations predict a dominant role of the permeability in the determination of the temporal evolution of the mechanical response of the tissue. In line with the previous studies on other tissues, the analysis shows that an ideal model based on Donnan’s equilibrium overestimates the osmotic pressure in skin for the case of very dilute solutions. The quadriphasic model is applied to predict changes in dermal cell environment and therefore alterations in what is called the “mechanome,” associated with skin stretch. The simulations indicate that skin deformation causes a variation in several local variables, including in particular the electric field associated with a deformation-induced non-homogeneous distribution of fixed charges.

Aneurysmal rupture is associated with wall thinning, but the mechanism is poorly understood. This study aimed to characterize the three-dimensional wall-thickness distributions of unruptured intracranial aneurysms. Five aneurysmal tissues were investigated using micro-computed tomography. First, the wall thickness was related to the aneurysmal wall appearances during surgery. The median wall thicknesses of the translucent and non-translucent walls were 50.56 and 155.93 µm, respectively (p < 0.05) with significant variation in the non-translucent wall thicknesses (p < 0.05). The three-dimensional observations characterized the spatial variation of wall thicknesses. Thin walls showed a uniform thickness profile ranging from 10 to 40 µm, whereas thick walls presented a peaked thickness profile ranging from 300 to 500 µm. In transition walls, the profile undulated due to the formation of focal thin/thick spots. Overall, the aneurysmal wall thicknesses were strongly site-dependent and spatially varied by 10 to 40 times within individual cases. Aneurysmal walls are exposed to wall stress driven by blood pressure. In theory, the magnitude of wall stress is inversely proportional to wall thickness. Thus, the observed spatial variation of wall thickness may increase the spatial variation of wall stress to a similar extent. The irregular wall thickness may yield stress concentration. The observed thin walls and focal thin spots may be caused by excessive wall stresses at the range of mechanical failure inducing wall injuries, such as microscopic tears, during aneurysmal enlargement. The present results suggested that blood pressure (wall stress) may have a potential of acting as a trigger of aneurysmal wall injury.

Understanding the viscoelastic behavior of pediatric brain tissue is critical to interpret how external mechanical forces affect head injury in children. However, knowledge of the viscoelastic properties of pediatric brain tissue is limited, and this reduces the biofidelity of developed numeric simulations of the pediatric head in analysis of brain injury. Thus, it is essential to characterize the viscoelastic behavior of pediatric brain tissue in various loading conditions and to identify constitutive models. In this study, the pediatric porcine brain tissue was investigated in compression with determine the viscoelasticity under small and large strain, respectively. A range of frequencies between 0.1 and 40 Hz was applied to determine frequency-dependent viscoelastic behavior via dynamic mechanical analysis, while brain samples were divided into three strain rate groups of 0.01/s, 1/s and 10/s for compression up to 0.3 strain level and stress relaxation to obtain time-dependent viscoelastic properties. At frequencies above 20 Hz, the storage modulus did not increase, while the loss modulus increased continuously. With strain rate increasing from 0.01/s to 10/s, the mean stress at 0.1, 0.2 and 0.3 strain increased to approximate 6.8, 5.6 and 4.4 times, respectively. The brain compressive response was sensitive to strain rate and frequency. The characterization of brain tissue will be valuable for development of head protection systems and prediction of brain injury.