Search

The development of biophysical models for clinical applications is rapidly advancing in the research community, thanks to their predictive nature and their ability to assist the interpretation of clinical data. However, high-resolution and accurate multi-physics computational models are computationally expensive and their personalisation involves fine calibration of a large number of parameters, which may be space-dependent, challenging their clinical translation. In this work, we propose a new approach which relies on the combination of physics-informed neural networks (PINNs) with three-dimensional soft tissue nonlinear biomechanical models, capable of reconstructing displacement fields and estimating heterogeneous patient-specific biophysical properties. The proposed learning algorithm encodes information from a limited amount of displacement and, in some cases, strain data, that can be routinely acquired in the clinical setting, and combines it with the physics of the problem, represented by a mathematical model based on partial differential equations, to regularise the problem and improve its convergence properties. Several benchmarks are presented to show the accuracy and robustness of the proposed method and its great potential to enable the robust and effective identification of patient-specific, heterogeneous physical properties, s.a. tissue stiffness properties. In particular, we demonstrate the capability of the PINN to detect the presence, location and severity of scar tissue, which is beneficial to develop personalised simulation models for disease diagnosis, especially for cardiac applications.

Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D–0D model to a limit cycle under baseline conditions, the model’s ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible., Europäische Kommission 750835, Fonds zur Förderung der Wissenschaftlichen Forschung (DE-588)2054142-9 F3210-N18 I2760-B30

Computational modeling of cardiovascular biomechanics should generally start from a homeostatic state. This is particularly relevant for image-based modeling, where the reference configuration is the loaded in vivo state obtained from imaging. This state includes residual stress of the vascular constituents, as well as anisotropy from the spatially varying orientation of collagen and smooth muscle fibers. Estimation of the residual stress and fiber orientation fields is a formidable challenge in realistic applications. To help address this challenge, we herein develop a growth based Algorithm to recover a residual stress distribution in vascular domains such that the stress state in the loaded configuration is equal to a prescribed homeostatic stress distribution at physiologic pressure. A stress-driven fiber deposition process is included in the framework, which defines the distribution of the fiber alignments in the vascular homeostatic state based on a minimization procedure. Numerical simulations are conducted to test this two-stage homeostasis generation algorithm in both idealized and non-idealized geometries, yielding results that agree favorably with prior numerical and experimental data.

The impact of increased stiffness and pulsatile load on the circulation and their influence on heart performance have been documented not only for cardiovascular events but also for ventricular dysfunctions. For this reason, computer models of cardiac electromechanics (EM) have to integrate effects of the circulatory system on heart function to be relevant for clinical applications. Currently it is not feasible to consider three-dimensional (3D) models of the entire circulation. Instead, simplified representations of the circulation are used, ensuring a satisfactory trade-off between accuracy and computational cost. In this work, we propose a novel and stable strategy to couple a 3D EM model of the heart to a one-dimensional (1D) model of blood flow in the arterial system. A personalised coupled 3D-1D model of LV and arterial system is built and used in a numerical benchmark to demonstrate robustness and accuracy of our scheme over a range of time steps. Validation of the coupled model is performed by investigating the coupled system's physiological response to variations in the arterial system affecting pulse wave propagation, comprising aortic stiffening, aortic stenosis or bifurcations causing wave reflections. Our results show that the coupled 3D-1D model is robust, stable and correctly replicates known physiology. In comparison with standard coupled 3D-0D models, additional computational costs are negligible, thus facilitating the use of our coupled 3D-1D model as a key methodology in studies where wave propagation effects are under investigation.

Fiber-reinforced soft biological tissues are typically modeled as hyperelastic, anisotropic, and nearly incompressible materials. To enforce incompressibility a multiplicative split of the deformation gradient into a volumetric and an isochoric part is a very common approach. However, due to the high stiffness of anisotropic materials in the preferred directions, the finite element analysis of such problems often suffers from severe locking effects and numerical instabilities. In this paper, we present novel methods to overcome locking phenomena for anisotropic materials using stabilized P1-P1 elements. We introduce different stabilization techniques and demonstrate the high robustness and computational efficiency of the chosen methods. In several benchmark problems we compare the approach to standard linear elements and show the accuracy and versatility of the methods to simulate anisotropic, nearly and fully incompressible materials. We are convinced that this numerical framework offers the possibility to accelerate accurate simulations of biological tissues, enabling patient-specfic parameterization studies, which require numerous forward simulations., Comment: This research has received funding from the European Union's Horizon 2020 research and innovation programme under the ERA-NET co-fund action No. 680969 (ERA-CVD SICVALVES) funded by the Austrian Science Fund (FWF), Grant I 4652-B

Computational models of the heart are increasingly being used in the development of devices, patient diagnosis and therapy guidance. While software techniques have been developed for simulating single hearts, there remain significant challenges in simulating cohorts of virtual hearts from multiple patients. To facilitate the development of new simulation and model analysis techniques by groups without direct access to medical data, image analysis techniques and meshing tools, we have created the first publicly available virtual cohort of twenty-four four-chamber hearts. Our cohort was built from heart failure patients, age 67±14 years. We segmented four-chamber heart geometries from end-diastolic (ED) CT images and generated linear tetrahedral meshes with an average edge length of 1.1±0.2mm. Ventricular fibres were added in the ventricles with a rule-based method with an orientation of -60° and 80° at the epicardium and endocardium, respectively. We additionally refined the meshes to an average edge length of 0.39±0.10mm to show that all given meshes can be resampled to achieve an arbitrary desired resolution. We ran simulations for ventricular electrical activation and free mechanical contraction on all 1.1mm-resolution meshes to ensure that our meshes are suitable for electro-mechanical simulations. Simulations for electrical activation resulted in a total activation time of 149±16ms. Free mechanical contractions gave an average left ventricular (LV) and right ventricular (RV) ejection fraction (EF) of 35±1% and 30±2%, respectively, and a LV and RV stroke volume (SV) of 95±28mL and 65±11mL, respectively. By making the cohort publicly available, we hope to facilitate large cohort computational studies and to promote the development of cardiac computational electro-mechanics for clinical applications., Version of record

FEniCS Mechanics is a Python package to facilitate computational mechanics simulations. The Python library dolfin, from the FEniCS Project, is used to formulate and numerically solve the problem in variational form. The general balance laws from continuum mechanics are used to enable rapid prototyping of different material laws. In addition to its generality, FEniCS Mechanics also checks the input provided by users to ensure that problem definitions are physically consistent. In turn, this code enables simulations of custom mechanics problems to be more accessible to those with limited programming or mechanics knowledge., Comment: This work was supported in part from the NSF award 1663747 and the NSF GRFP. Additionally, this project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Action H2020-MSCA-IF-2016 InsiliCardio, GA No. 750835

To observe internalization of the yeast pheromone receptor Ste2 by fluorescence microscopy in live cells in real time, we visualized only those molecules present at the cell surface at the time of agonist engagement (rather than the total cellular pool) by tagging this receptor at its N-terminus with an exocellular fluorogen-activating protein (FAP). A FAP is a single-chain antibody engineered to bind tightly a nonfluorescent, cell-impermeable dye (fluorogen), thereby generating a fluorescent complex. The utility of FAP tagging to study trafficking of integral membrane proteins in yeast, which possesses a cell wall, had not been examined previously. A diverse set of signal peptides and propeptide sequences were explored to maximize expression. Maintenance of the optimal FAP-Ste2 chimera intact required deletion of two, paralogous, glycosylphosphatidylinositol (GPI)-anchored extracellular aspartyl proteases (Yps1 and Mkc7). FAP-Ste2 exhibited a much brighter and distinct plasma membrane signal than Ste2-GFP or Ste2-mCherry yet behaved quite similarly. Using FAP-Ste2, new information was obtained about the mechanism of its internalization, including novel insights about the roles of the cargo-selective endocytic adaptors Ldb19/Art1, Rod1/Art4, and Rog3/Art7.

Saccharomyces cerevisiae target of rapamycin (TOR) complex 2 (TORC2) is an essential regulator of plasma membrane lipid and protein homeostasis. How TORC2 activity is modulated in response to changes in the status of the cell envelope is unclear. Here we document that TORC2 subunit Avo2 is a direct target of Slt2, the mitogen-activated protein kinase (MAPK) of the cell wall integrity pathway. Activation of Slt2 by overexpression of a constitutively active allele of an upstream Slt2 activator (Pkc1) or by auxin-induced degradation of a negative Slt2 regulator (Sln1) caused hyperphosphorylation of Avo2 at its MAPK phosphoacceptor sites in a Slt2-dependent manner and diminished TORC2-mediated phosphorylation of its major downstream effector, protein kinase Ypk1. Deletion of Avo2 or expression of a phosphomimetic Avo2 allele rendered cells sensitive to two stresses (myriocin treatment and elevated exogenous acetic acid) that the cell requires Ypk1 activation by TORC2 to survive. Thus, Avo2 is necessary for optimal TORC2 activity, and Slt2-mediated phosphorylation of Avo2 down-regulates TORC2 signaling. Compared with wild-type Avo2, phosphomimetic Avo2 shows significant displacement from the plasma membrane, suggesting that Slt2 inhibits TORC2 by promoting Avo2 dissociation. Our findings are the first demonstration that TORC2 function is regulated by MAPK-mediated phosphorylation.

To observe internalization of the yeast pheromone receptor Ste2 by fluorescence microscopy in live cells in real time, we visualized only those molecules present at the cell surface at the time of agonist engagement (rather than the total cellular pool) by tagging this receptor at its N-terminus with an exocellular fluorogen-activating protein (FAP). A FAP is a single-chain antibody engineered to bind tightly a nonfluorescent, cell-impermeable dye (fluorogen), thereby generating a fluorescent complex. The utility of FAP tagging to study trafficking of integral membrane proteins in yeast, which possesses a cell wall, had not been examined previously. A diverse set of signal peptides and propeptide sequences were explored to maximize expression. Maintenance of the optimal FAP-Ste2 chimera intact required deletion of two, paralogous, glycosylphosphatidylinositol (GPI)-anchored extracellular aspartyl proteases (Yps1 and Mkc7). FAP-Ste2 exhibited a much brighter and distinct plasma membrane signal than Ste2-GFP or Ste2-mCherry yet behaved quite similarly. Using FAP-Ste2, new information was obtained about the mechanism of its internalization, including novel insights about the roles of the cargo-selective endocytic adaptors Ldb19/Art1, Rod1/Art4, and Rog3/Art7.

Saccharomyces cerevisiae target of rapamycin (TOR) complex 2 (TORC2) is an essential regulator of plasma membrane lipid and protein homeostasis. How TORC2 activity is modulated in response to changes in the status of the cell envelope is unclear. Here we document that TORC2 subunit Avo2 is a direct target of Slt2, the mitogen-activated protein kinase (MAPK) of the cell wall integrity pathway. Activation of Slt2 by overexpression of a constitutively active allele of an upstream Slt2 activator (Pkc1) or by auxin-induced degradation of a negative Slt2 regulator (Sln1) caused hyperphosphorylation of Avo2 at its MAPK phosphoacceptor sites in a Slt2-dependent manner and diminished TORC2-mediated phosphorylation of its major downstream effector, protein kinase Ypk1. Deletion of Avo2 or expression of a phosphomimetic Avo2 allele rendered cells sensitive to two stresses (myriocin treatment and elevated exogenous acetic acid) that the cell requires Ypk1 activation by TORC2 to survive. Thus, Avo2 is necessary for optimal TORC2 activity, and Slt2-mediated phosphorylation of Avo2 down-regulates TORC2 signaling. Compared with wild-type Avo2, phosphomimetic Avo2 shows significant displacement from the plasma membrane, suggesting that Slt2 inhibits TORC2 by promoting Avo2 dissociation. Our findings are the first demonstration that TORC2 function is regulated by MAPK-mediated phosphorylation.

High-resolution and anatomically realistic computer models of biological soft tissues play a significant role in the understanding of the function of cardiovascular components in health and disease. However, the computational effort to handle fine grids to resolve the geometries as well as sophisticated tissue models is very challenging. One possibility to derive a strongly scalable parallel solution algorithm is to consider finite element tearing and interconnecting (FETI) methods. In this study, we propose and investigate the application of FETI methods to simulate the elastic behavior of biological soft tissues. As one particular example, we choose the artery which is—as most other biological tissues—characterized by anisotropic and nonlinear material properties. We compare two specific approaches of FETI methods, classical and all-floating, and investigate the numerical behavior of different preconditioning techniques. In comparison with classical FETI, the all-floating approach not only has advantages concerning the implementation but also has advantages concerning the convergence of the global iterative solution method. This behavior is illustrated with numerical examples. We present results of linear elastic simulations to show convergence rates, as expected from the theory, and results from the more sophisticated nonlinear case where we apply a well-known anisotropic model to the realistic geometry of an artery. Although the FETI methods have a great applicability on artery simulations, we will also discuss some limitations concerning the dependence on material parameters.

High-resolution and anatomically realistic computer models of biological soft tissues play a significant role in the understanding of the function of cardiovascular components in health and disease. However, the computational effort to handle fine grids to resolve the geometries as well as sophisticated tissue models is very challenging. One possibility to derive a strongly scalable parallel solution algorithm is to consider finite element tearing and interconnecting (FETI) methods. In this study, we propose and investigate the application of FETI methods to simulate the elastic behavior of biological soft tissues. As one particular example, we choose the artery which is—as most other biological tissues—characterized by anisotropic and nonlinear material properties. We compare two specific approaches of FETI methods, classical and all-floating, and investigate the numerical behavior of different preconditioning techniques. In comparison with classical FETI, the all-floating approach not only has advantages concerning the implementation but also has advantages concerning the convergence of the global iterative solution method. This behavior is illustrated with numerical examples. We present results of linear elastic simulations to show convergence rates, as expected from the theory, and results from the more sophisticated nonlinear case where we apply a well-known anisotropic model to the realistic geometry of an artery. Although the FETI methods have a great applicability on artery simulations, we will also discuss some limitations concerning the dependence on material parameters.

Computational modeling of ventricular electromechanics is considered to be among the most promising approaches to gain novel insight into cardiac function in health and disease. Such models allow to integrate the wealth of available experimental data into a mechanistic framework which allows to study complex cause-effect relationships across several scales of biological organization – ranging from subcellular processes such as cellular force genera-tion up to the organ scale. However, due to the enormous challenges involved in constructing such models, the number of available multiscale models is very limited. Their implementation resorts to very crude approximations which ignores many aspects of known biophysical details and/or uses extremely coarse spatio-temporal discretizations to reduce computational complexity. In this study we focus on a subset of the parametrization and validation issues in models of cardiac electromechanics. An experimental setup for measuring passive stress-strain relationships as well as isometric force transients in thin preparations of ventricular trabeculae or papillary muscles was developed along with a matching computer model to parametrize and validate passive and active mechanical model responses at the tissue scale., G. Plank, A.J. Prassl, R. Arnold, Y. Rezk, T.E. Fastl, E. Hofer and C.M. Augustin

Computational modeling of ventricular electromechanics is considered to be among the most promising approaches to gain novel insight into cardiac function in health and disease. Such models allow to integrate the wealth of available experimental data into a mechanistic framework which allows to study complex cause-effect relationships across several scales of biological organization – ranging from subcellular processes such as cellular force genera-tion up to the organ scale. However, due to the enormous challenges involved in constructing such models, the number of available multiscale models is very limited. Their implementation resorts to very crude approximations which ignores many aspects of known biophysical details and/or uses extremely coarse spatio-temporal discretizations to reduce computational complexity. In this study we focus on a subset of the parametrization and validation issues in models of cardiac electromechanics. An experimental setup for measuring passive stress-strain relationships as well as isometric force transients in thin preparations of ventricular trabeculae or papillary muscles was developed along with a matching computer model to parametrize and validate passive and active mechanical model responses at the tissue scale., G. Plank, A.J. Prassl, R. Arnold, Y. Rezk, T.E. Fastl, E. Hofer and C.M. Augustin

Durch Christoph Augustin/ ..., Schlüsselseiten aus Exemplar AB 56497 (3) (ULB Halle)

Durch Christoph Augustin/ ..., Schlüsselseiten aus Exemplar AB 56497 (3) (ULB Halle)