Correction to: Mechanical factors contributing to the Venus flytrap’s rate-dependent response to stimuli

The sensory hairs of the Venus flytrap (Dionaea muscipula Ellis) detect mechanical stimuli imparted by their prey and fire bursts of electrical signals called action potentials (APs). APs are elicited when the hairs are sufficiently stimulated and two consecutive APs can trigger closure of the trap. Earlier experiments have identified thresholds for the relevant stimulus parameters, namely the angular displacement \(\theta \) and angular velocity \(\omega \). However, these experiments could not trace the deformation of the trigger hair’s sensory cells, which are known to transduce the mechanical stimulus. To understand the kinematics at the cellular level, we investigate the role of two relevant mechanical phenomena: viscoelasticity and intercellular fluid transport using a multi-scale numerical model of the sensory hair. We hypothesize that the combined influence of these two phenomena and \(\omega \) contribute to the flytrap’s rate-dependent response to stimuli. In this study, we firstly perform sustained deflection tests on the hair to estimate the viscoelastic material properties of the tissue. Thereafter, through simulations of hair deflection tests at different loading rates, we were able to establish a multi-scale kinematic link between \(\omega \) and the cell wall stretch \(\delta \). Furthermore, we find that the rate at which \(\delta \) evolves during a stimulus is also proportional to \(\omega \). This suggests that mechanosensitive ion channels, expected to be stretch-activated and localized in the plasma membrane of the sensory cells, could be additionally sensitive to the rate at which stretch is applied.

     

On a three–dimensional model for the description of the passive characteristics of skeletal muscle tissue

In this work, a three–dimensional model was developed to describe the passive mechanical behaviour of anisotropic skeletal muscle tissue. To validate the model, orientation–dependent axial (\(0^\circ\), \(45^\circ\), \(90^\circ\)) and semi–confined compression experiments (mode I, II, III) were performed on soleus muscle tissue from rabbits. In the latter experiments, specimen deformation is prescribed in the loading direction and prevented in an additional spatial direction, fibre compression at \(0^\circ\) (mode I), fibre elongation at \(90^\circ\) (mode II) and a neutral state of the fibres at \(90^\circ\) where their length is kept constant (mode III). Overall, the model can adequately describe the mechanical behaviour with a relatively small number of model parameters. The stiffest tissue response during orientation–dependent axial compression (\(-\,7.7\,\pm \,1.3\) kPa) occurs when the fibres are oriented perpendicular to the loading direction (\(90^\circ\)) and are thus stretched during loading. Semi–confined compression experiments yielded the stiffest tissue (\(-\,36.7\,\pm \,11.2\) kPa) in mode II when the muscle fibres are stretched. The extensive data set collected in this study allows to study the different error measures depending on the deformation state or the combination of deformation states.

     

A thermomechanical framework for reconciling the effects of ultraviolet radiation exposure time and wavelength on connective tissue elasticity

Augmentation of the mechanical properties of connective tissue using ultraviolet (UV) radiation—by targeting collagen cross-linking in the tissue at predetermined UV exposure time \((t)\) and wavelength \((\lambda )\) —has been proposed as a therapeutic method for supporting the treatment for structural-related injuries and pathologies. However, the effects of \(\lambda \) and \(t\) on the tissue elasticity, namely elastic modulus \((E)\) and modulus of resilience \((u_\mathrm{Y})\) , are not entirely clear. We present a thermomechanical framework to reconcile the \(t\) - and \(\lambda \) -related effects on \(E\) and \(u_\mathrm{Y}\) . The framework addresses (1) an energy transfer model to describe the dependence of the absorbed UV photon energy, \(\xi \) , per unit mass of the tissue on \(t\) and \(\lambda \) , (2) an intervening thermodynamic shear-related parameter, \(G\) , to quantify the extent of UV-induced cross-linking in the tissue, (3) a threshold model for the \(G\) versus \(\xi \) relationship, characterized by   \(t_\mathrm{C}\) —the critical \(t\) underpinning the association of \(\xi \) with \(G\) —and (4) the role of \(G\) in the tissue elasticity. We hypothesized that \(G\) regulates \(E\) (UV-stiffening hypothesis) and \(u_\mathrm{Y}\) (UV-resilience hypothesis). The framework was evaluated with the support from data derived from tensile testing on isolated ligament fascicles, treated with two levels of \(\lambda \) (365 and 254 nm) and three levels of \(t\) (15, 30 and 60 min). Predictions from the energy transfer model corroborated the findings from a two-factor analysis of variance of the effects of \(t\) and \(\lambda \) treatments. Student’s t test revealed positive change in \(E\) and \(u_\mathrm{Y}\) with increases in \(G\) —the findings lend support to the hypotheses, implicating the implicit dependence of UV-induced cross-links on \(t\) and \(\lambda \) for directing tissue stiffness and resilience. From a practical perspective, the study is a step in the direction to establish a UV irradiation treatment protocol for effective control of exogenous cross-linking in connective tissues.     

Modeling the debonding process of osseointegrated implants due to coupled adhesion and friction

Cementless implants have become widely used for total hip replacement surgery. The long-term stability of these implants is achieved by bone growing around and into the rough surface of the implant, a process called osseointegration. However, debonding of the bone–implant interface can still occur due to aseptic implant loosening and insufficient osseointegration, which may have dramatic consequences. The aim of this work is to describe a new 3D finite element frictional contact formulation for the debonding of partially osseointegrated implants. The contact model is based on a modified Coulomb friction law by Immel et al. (2020), that takes into account the tangential debonding of the bone-implant interface. This model is extended in the direction normal to the bone-implant interface by considering a cohesive zone model, to account for adhesion phenomena in the normal direction and for adhesive friction of partially bonded interfaces. The model is applied to simulate the debonding of an acetabular cup implant. The influence of partial osseointegration and adhesive effects on the long-term stability of the implant is assessed. The influence of different patient- and implant-specific parameters such as the friction coefficient \(\mu _\text {b}\), the trabecular Young’s modulus \(E_\text {b}\), and the interference fit \(I\!F\) is also analyzed, in order to determine the optimal stability for different configurations. Furthermore, this work provides guidelines for future experimental and computational studies that are necessary for further parameter calibration.

     

On real-valued SDE and nonnegative-valued SDE population models with demographic variability

Population dynamics with demographic variability is frequently studied using discrete random variables with continuous-time Markov chain (CTMC) models. An approximation of a CTMC model using continuous random variables can be derived in a straightforward manner by applying standard methods based on the reaction rates in the CTMC model. This leads to a system of Itô stochastic differential equations (SDEs) which generally have the form \( d \mathbf {y} = \varvec{\mu } \, dt + G \, d \mathbf {W},\) where \(\mathbf {y}\) is the population vector of random variables, \(\varvec{\mu }\) is the drift vector, and G is the diffusion matrix. In some problems, the derived SDE model may not have real-valued or nonnegative solutions for all time. For such problems, the SDE model may be declared infeasible. In this investigation, new systems of SDEs are derived with real-valued solutions and with nonnegative solutions. To derive real-valued SDE models, reaction rates are assumed to be nonnegative for all time with negative reaction rates assigned probability zero. This biologically realistic assumption leads to the derivation of real-valued SDE population models. However, small but negative values may still arise for a real-valued SDE model. This is due to the magnitudes of certain problem-dependent diffusion coefficients when population sizes are near zero. A slight modification of the diffusion coefficients when population sizes are near zero ensures that a real-valued SDE model has a nonnegative solution, yet maintains the integrity of the SDE model when sizes are not near zero. Several population dynamic problems are examined to illustrate the methodology.

     

Competition of Multiple Species in Advective Environments

We study the effect of changes in flow speed on competition of an arbitrary number of species living in advective environments, such as streams and rivers. We begin with a spatial Lotka–Volterra model which is described by n reaction–diffusion–advection equations with Danckwerts boundary conditions. Using the dominant eigenvalue \(\lambda \le 0\) of the diffusion–advection operator subject to boundary conditions, we reduce the model to a system of ordinary differential equations. We impose a “transitive arrangement” of the competitors in terms of their interspecific coefficients and growth rates, which means that in the absence of advection, we have the following situation: for all \(1\le i<j\le n\), species i out-competes species j, while species j has higher intrinsic growth rate than species i. Changing advection speed in the original spatial model corresponds to changing the value of \(\lambda \) in the spatially implicit model. Considering the cases of the odd and even n separately, we obtain explicit intervals of the values of \(\lambda \) that allow all n species to be present in the habitat (coexistence interval). Stability of this equilibrium is shown for \(n\le 4\).     

Similarity of decay-associated spectra for tryptophan fluorescence of proteins with different structures

Tryptophan fluorescence lifetimes were analyzed for three proteins: human serum albumin, bovine serum albumin, and bacterial luciferase, which contain one, two, and seven tryptophan residues, respectively. For all of the proteins, the fluorescence decays were fitted by three lifetimes: τ₁ = 6–7 ns, τ₂ = 2.0–2.3 ns, and τ₃ ≤ 0.1 ns (the native state), and τ₁ = 4.4–4.6 ns, τ₂ = 1.7–1.8 ns, and τ₃ ≤ 0.1 ns (the denatured state). Corresponding decay-associated spectra had similar peak wavelengths and spectrum half-widths both in the native state (\(\lambda _{\max }^{{\tau _1}} = 324nm\), \(\lambda _{\max }^{{\tau _2}} = 328nm\), and \(\lambda _{\max }^{{\tau _3}} = 315nm\)), and in the denatured state (\(\lambda _{\max }^{{\tau _1}} = 350nm\), \(\lambda _{\max }^{{\tau _2}} = 343nm\), and \(\lambda _{\max }^{{\tau _3}} = 317nm\)). The differences in the steady-state spectra of the studied proteins were accounted for the individual ratio of the lifetime component contributions. The lifetime components were compared with a classification of tryptophan residues in the structure of these proteins within the discrete states model.     

Quantifying interspecific spatial overlap in aquatic macrophyte communities

Levins’s asymmetrical α index quantifies between species overlap over resources more realistically than similar-purpose single-value indices. The associated community-wide \(\bar \alpha\) index expresses the degree of “species packing”. Both indices were formulated upon competing animal (i.e., mobile) organisms and are independent of population densities. However, overlap over resources for nonmobile organisms such as plants may have an impact even below carrying capacity. The proposed \(\hat \alpha\) index, based on Levins’s α index, quantifies spatial overlap for plants integrating information on species spatial distribution and crowding conditions. The \(\hat \alpha\) index is specifically designed for plant distribution data collected in discrete plots with density expressed as percent coverage (%cover) of substratum. We also propose a community-wide \({\hat \alpha_{\text{c}}}\) index, conceptually analogous to \(\bar \alpha\) , but furnished with a measure of dispersion (se \({\hat \alpha_{\text{c}}}\) ). Species importance within the community is inferred from comparisons of pairwise \(\hat \alpha\) ’s with \({\hat \alpha_{\text{c}}}\) . The \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices correlate closely and exponentially with plant density, and correct apparent over- and underestimations of interaction intensity at low and very high crowding by Levins’s α and \(\bar \alpha\) , respectively. Index application to aquatic plant communities gave results consistent with within-community and general ecological patterns, suggesting a high potential of the proposed \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices in basic and applied macrophyte ecological studies and management.     

Glutathione <Emphasis Type="BoldItalic">S</Emphasis>-transferase P1 gene polymorphisms and susceptibility to coronary artery disease in a subgroup of north Indian population

The present study aimed to investigate the association of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 with coronary artery disease (CAD) in a subgroup of north Indian population. In the present case–control study, CAD patients (\(n = 200\)) and age-matched, sex-matched and ethnicity-matched healthy controls (\(n = 200\)) were genotyped for polymorphisms in GSTP1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Genotype distribution of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 gene was significantly different between cases and controls (\(P = 0.005\) and 0.024, respectively). Binary logistic regression analysis showed significant association of A/G (odds ratio (OR): 1.6, 95% CI: 1.08–2.49, \(P = 0.020\)) and G/G (OR: 3.1, 95% CI: 1.41–6.71, P \(=\) 0.005) genotypes of GSTP1 \(\hbox {g}.313\hbox {A}{\!>\!}\hbox {G}\), and C/T (OR: 5.8, 95% CI: 1.26–26.34, \(P = 0.024\)) genotype of GSTP1 \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) with CAD. The A/G and G/G genotypes of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and C/T genotype of \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) conferred 6.5-fold increased risk for CAD (OR: 6.5, 95% CI: 1.37–31.27, \(P = 0.018\)). Moreover, the recessive model of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) is the best fit inheritance model to predict the susceptible gene effect (OR: 2.3, 95% CI: 1.11–4.92, \(P = 0.020\)). In conclusion, statistically significant associations of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) (A/G, G/G) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) (C/T) genotypes with CAD were observed.     

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\(0.36-0.87\)) with a strong correlation with the organ length-scale radial strain (\(R^{2}= 0.95\)); conversely, circumferential cell strains spanned a very narrow range (\(0.75-0.88\)) showing no correlation with the circumferential strain at the organ level (\(R^{2}= 0.02\)). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

A multiphase model of growth factor-regulated atherosclerotic cap formation

Atherosclerosis is characterised by the growth of fatty plaques in the inner artery wall. In mature plaques, vascular smooth muscle cells (SMCs) are recruited from adjacent tissue to deposit a collagenous cap over the fatty plaque core. This cap isolates the thrombogenic plaque content from the bloodstream and prevents the clotting cascade that leads to myocardial infarction or stroke. Despite the protective role of the cap, the mechanisms that regulate cap formation and maintenance are not well understood. It remains unclear why some caps become stable, while others become vulnerable to rupture. We develop a multiphase PDE model with non-standard boundary conditions to investigate collagen cap formation by SMCs in response to diffusible growth factor signals from the endothelium. Platelet-derived growth factor stimulates SMC migration, proliferation and collagen degradation, while transforming growth factor (TGF)-\(\beta \) stimulates SMC collagen synthesis and inhibits collagen degradation. The model SMCs respond haptotactically to gradients in the collagen phase and have reduced rates of migration and proliferation in dense collagenous tissue. The model, which is parameterised using in vivo and in vitro experimental data, reproduces several observations from plaque growth in mice. Numerical and analytical results demonstrate that a stable cap can be formed by a relatively small SMC population and emphasise the critical role of TGF-\(\beta \) in effective cap formation. These findings provide unique insight into the mechanisms that may lead to plaque destabilisation and rupture. This work represents an important step towards the development of a comprehensive in silico plaque model.

     

Modulation of the cAMP Response by G\alpha _i and G\beta \gamma : A Computational Study of G Protein Signaling in Immune Cells

Cyclic AMP is important for the resolution of inflammation, as it promotes anti-inflammatory signaling in several immune cell lines. In this paper, we present an immune cell specific model of the cAMP signaling cascade, paying close attention to the specific isoforms of adenylyl cyclase (AC) and phosphodiesterase that control cAMP production and degradation, respectively, in these cells. The model describes the role that G protein subunits, including G \(\alpha _s\) , G \(\alpha _i\) , and G \(\beta \gamma \) , have in regulating cAMP production. Previously, G \(\alpha _i\) activation has been shown to increase the level of cAMP in certain immune cell types. This increase in cAMP is thought to be mediated by \(\beta \gamma \) subunits which are released upon G \(\alpha \) activation and can directly stimulate specific isoforms of AC. We conduct numerical experiments in order to explore the mechanisms through which G \(\alpha _i\) activation can increase cAMP production. An important conclusion of our analysis is that the relative abundance of different G protein subunits is an essential determinant of the cAMP profile in immune cells. In particular, our model predicts that limited availability of \(\beta \gamma \) subunits may both \((i)\) enable immune cells to link inflammatory G \(\alpha _i\) signaling to anti-inflammatory cAMP production thereby creating a balanced immune response to stimulation with low concentrations of PGE2, and \((ii)\) prohibit robust anti-inflammatory cAMP signaling in response to stimulation with high concentrations of PGE2.     

Stability of a Stochastic Model of an SIR Epidemic with Vaccination

We prove almost sure exponential stability for the disease-free equilibrium of a stochastic differential equations model of an SIR epidemic with vaccination. The model allows for vertical transmission. The stochastic perturbation is associated with the force of infection and is such that the total population size remains constant in time. We prove almost sure positivity of solutions. The main result concerns especially the smaller values of the diffusion parameter, and describes the stability in terms of an analogue \(\mathcal{R}_\sigma\) of the basic reproduction number \(\mathcal{R}_0\) of the underlying deterministic model, with \(\mathcal{R}_\sigma \le \mathcal{R}_0\). We prove that the disease-free equilibrium is almost sure exponentially stable if \(\mathcal{R}_\sigma <1\).     

Cellular mechanosensitivity to substrate stiffness decreases with increasing dissimilarity to cell stiffness

Computational modelling has received increasing attention to investigate multi-scale coupled problems in micro-heterogeneous biological structures such as cells. In the current study, we investigated for a single cell the effects of (1) different cell-substrate attachment (2) and different substrate modulus \(\textit{E}_\mathrm{s}\) on intracellular deformations. A fibroblast was geometrically reconstructed from confocal micrographs. Finite element models of the cell on a planar substrate were developed. Intracellular deformations due to substrate stretch of \(\lambda =1.1\), were assessed for: (1) cell-substrate attachment implemented as full basal contact (FC) and 124 focal adhesions (FA), respectively, and \(\textit{E}_\mathrm{s}\,=\,\)140 KPa and (2) \(\textit{E}_\mathrm{s}\,=\,10\), 140, 1000, and 10,000 KPa, respectively, and FA attachment. The largest strains in cytosol, nucleus and cell membrane were higher for FC (1.35\(\text {e}^{-2}\), 0.235\(\text {e}^{-2}\) and 0.6\(\text {e}^{-2}\)) than for FA attachment (0.0952\(\text {e}^{-2}\), 0.0472\(\text {e}^{-2}\) and 0.05\(\text {e}^{-2}\)). For increasing \(\textit{E}_\mathrm{s}\), the largest maximum principal strain was 4.4\(\text {e}^{-4}\), 5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\) and 5.3\(\text {e}^{-4}\) in the membrane, 9.5\(\text {e}^{-4}\), 1.1\(\text {e}^{-4}\), 1.2\(\text {e}^{-3}\) and 1.2\(\text {e}^{-3}\) in the cytosol, and 4.5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\), 5.7\(\text {e}^{-4}\) and 5.7\(\text {e}^{-4}\) in the nucleus. The results show (1) the importance of representing FA in cell models and (2) higher cellular mechanical sensitivity for substrate stiffness changes in the range of cell stiffness. The latter indicates that matching substrate stiffness to cell stiffness, and moderate variation of the former is very effective for controlled variation of cell deformation. The developed methodology is useful for parametric studies on cellular mechanics to obtain quantitative data of subcellular strains and stresses that cannot easily be measured experimentally.     

Elastic behavior of a red blood cell with the membrane’s nonuniform natural state: equilibrium shape,motion transition under shear flow,and elongation during tank-treading motion

Direct numerical simulations of the mechanics of a single red blood cell (RBC) were performed by considering the nonuniform natural state of the elastic membrane. A RBC was modeled as an incompressible viscous fluid encapsulated by an elastic membrane. The in-plane shear and area dilatation deformations of the membrane were modeled by Skalak constitutive equation, while out-of-plane bending deformation was formulated by the spring model. The natural state of the membrane with respect to in-plane shear deformation was modeled as a sphere ( \(\alpha =0\) ), biconcave disk shape ( \(\alpha =1\) ) and their intermediate shapes ( \(0<\alpha <1\) ) with the nonuniformity parameter \(\alpha \) , while the natural state with respect to out-of-plane bending deformation was modeled as a flat plane. According to the numerical simulations, at an experimentally measured in-plane shear modulus of \(2.5\times 10^{-6}\,\hbox {N}/\hbox {m}\) and an out-of-plane bending rigidity of \(2.0\times 10^{-19}\,\hbox {N}\cdot \hbox {m}\) of the cell membrane, the following results were obtained. (i) The RBC shape at equilibrium was biconcave discoid for \(\alpha >0.22\) and cupped otherwise; (ii) the experimentally measured fluid shear stress at the transition between tumbling and tank-treading motions under shear flow was reproduced for \(0.05<\alpha <0.34\) ; (iii) the elongation deformation of the RBC during tank-treading motion from the simulation was consistent with that from in vitro experiments, irrespective of the \(\alpha \) value. Based on our RBC modeling, the three phenomena (i), (ii), and (iii) were mechanically consistent for \(0.22<\alpha <0.34\) . The condition \(0.05<\alpha <0.22\) precludes a biconcave discoid shape at equilibrium (i); however, it gives appropriate fluid shear stress at the motion transition under shear flow (ii), suggesting that a combined effect of \(\alpha \) and the natural state with respect to out-of-plane bending deformation is necessary for understanding details of the RBC mechanics at equilibrium. Our numerical results demonstrate that moderate nonuniformity in a membrane’s natural state with respect to in-plane shear deformation plays a key role in RBC mechanics.     

Generalized quasispecies model on finite metric spaces: isometry groups and spectral properties of evolutionary matrices

The quasispecies model introduced by Eigen in 1971 has close connections with the isometry group of the space of binary sequences relative to the Hamming distance metric. Generalizing this observation we introduce an abstract quasispecies model on a finite metric space X together with a group of isometries \(\Gamma \) acting transitively on X. We show that if the domain of the fitness function has a natural decomposition into the union of tG-orbits, G being a subgroup of \(\Gamma \), then the dominant eigenvalue of the evolutionary matrix satisfies an algebraic equation of degree at most \(t\cdot \mathrm{rk}_{\mathbf {Z}} R\), where R is the orbital ring that is defined in the text. The general theory is illustrated by three detailed examples. In the first two of them the space X is taken to be the metric space of vertices of a regular polytope with the natural “edge” metric, these are the cases of a regular m-gon and of a hyperoctahedron; the final example takes as X the quotient rings \(\mathbf {Z}/p^n\mathbf {Z}\) with p-adic metric.

     

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were \(C_1=3000\pm 1700\,\hbox {Pa}\) and \(\alpha =45\pm 25\) (\(b_f=25\pm 14\), \(b_{ft}=11\pm 6\), \(b_{t}=9\pm 5\)) for the HF cases and \(C_1=1700\,\hbox {Pa}\) and \(\alpha =15\) (\(b_f=8\), \(b_{ft}=4\), \(b_{t}=3\)) for the healthy case.