A patient specific electro-mechanical model of the heart

Effect of pulmonary regurgitation on cardiac functions based on a human bi-ventricle model

BACKGROUND AND OBJECTIVE

Assessing the severity of pulmonary regurgitation (PR) and identifying optimal clinically relevant indicators for its treatment is crucial, yet standards for quantifying PR remain unclear in clinical practice. Computational modelling of the heart is in the process of providing valuable insights and information for cardiovascular physiology research. However, the advancements of finite element computational models have not been widely applied to simulate cardiac outputs in patients with PR. Furthermore, a computational model that incorporates both the left ventricle (LV) and right ventricle (RV) can be valuable in assessing the relationship between left and right ventricular morphometry and septal motion in PR patients. To enhance our understanding of the effect of PR on cardiac functions and mechanical behaviour, we developed a human bi-ventricle model to simulate five cases with varying degrees of PR severity.

METHODS

This bi-ventricle model was built using a patient-specific geometry and a widely used myofibre architecture. The myocardial material properties were described by a hyperelastic passive constitutive law and a modified time-varying elastance active tension model. To simulate realistic cardiac functions and the dysfunction of the pulmonary valve in PR disease cases, open-loop lumped parameter models representing systemic and pulmonary circulatory systems were designed.

RESULTS

In the baseline case, pressures in the aorta and main pulmonary artery and ejection fractions of both the LV and RV were within normal physiological ranges reported in the literature. The end-diastolic volume (EDV) of the RV under varying degrees of PR was comparable to the reported cardiac magnetic resonance imaging data. Moreover, RV dilation and interventricular septum motion from the baseline to the PR cases were clearly observed through the long-axis and short-axis views of the bi-ventricle geometry. The RV EDV in the severe PR case increased by 50.3% compared to the baseline case, while the LV EDV decreased by 18.1%. The motion of the interventricular septum was consistent with the literature. Furthermore, ejection fractions of both the LV and RV decreased as PR became severe, with LV ejection fraction decreasing from 60.5% at baseline to 56.3% in the severe case and RV ejection fraction decreasing from 51.8% to 46.8%. Additionally, the average myofibre stress of the RV wall at end-diastole significantly increased due to PR, from 2.7±12.1 kPa at baseline to 10.9±26.5 kPa in the severe case. The average myofibre stress of the LV wall at end-diastole increased from 3.7±18.1 kPa to 4.3±20.3 kPa.

CONCLUSIONS

This study established a foundation for the computational modelling of PR. The simulated results showed that severe PR leads to reduced cardiac outputs in both the LV and RV, clearly observable septum motion, and a significant increase in the average myofibre stress in the RV wall. These findings demonstrate the potential of the model for further exploration of PR.

Countering reproducibility issues in mathematical models with software engineering techniques: A case study using a one-dimensional mathematical model of the atrioventricular node

One should assume that in silico experiments in systems biology are less susceptible to reproducibility issues than their wet-lab counterparts, because they are free from natural biological variations and their environment can be fully controlled. However, recent studies show that only half of the published mathematical models of biological systems can be reproduced without substantial effort. In this article we examine the potential causes for failed or cumbersome reproductions in a case study of a one-dimensional mathematical model of the atrioventricular node, which took us four months to reproduce. The model demonstrates that even otherwise rigorous studies can be hard to reproduce due to missing information, errors in equations and parameters, a lack in available data files, non-executable code, missing or incomplete experiment protocols, and missing rationales behind equations. Many of these issues seem similar to problems that have been solved in software engineering using techniques such as unit testing, regression tests, continuous integration, version control, archival services, and a thorough modular design with extensive documentation. Applying these techniques, we reimplement the examined model using the modeling language Modelica. The resulting workflow is independent of the model and can be translated to SBML, CellML, and other languages. It guarantees methods reproducibility by executing automated tests in a virtual machine on a server that is physically separated from the development environment. Additionally, it facilitates results reproducibility, because the model is more understandable and because the complete model code, experiment protocols, and simulation data are published and can be accessed in the exact version that was used in this article. We found the additional design and documentation effort well justified, even just considering the immediate benefits during development such as easier and faster debugging, increased understandability of equations, and a reduced requirement for looking up details from the literature.

A microstructure-based model for a full lamellar-interlamellar displacement and shear strain mapping inside human intervertebral disc core

The determinant role of the annulus fibrosus interlamellar zones in the intervertebral disc transversal and volumetric responses and hence on their corresponding three-dimensional conducts have been only revealed and appreciated recently. Their consideration in disc modeling strategies has been proven to be essential for the reproduction of correct local strain and displacement fields inside the disc especially in the unconstrained directions of the disc. In addition, these zones are known to be the starting areas of annulus fibrosus circumferential tears and disc delamination failure mode, which is often judged as one of the most dangerous disc failure modes that could evolve with time leading to disc hernia. For this latter reason, the main goal of the current contribution is to incorporate physically for the first time, the interlamellar zones, at the scale of a complete human lumbar intervertebral disc, in order to allow a correct local vision and replication of the different lamellar-interlamellar interactions and an identification of the interlamellar critical zones. By means of a fully tridimensional chemo-viscoelastic constitutive model, which we implemented into a finite element code, the physical, mechanical and chemical contribution of the interlamellar zones is added to the disc. The chemical-induced volumetric response is accounted by the model for both the interlamellar zones and the lamellae using experimentally-based fluid kinetics. Computational simulations are performed and critically discussed upon different simple and complex physiological movements. The disc core and the interlamellar zones are numerically accessed, allowing the observation of the displacement and shear strain fields that are compared to direct MRI experiments from the literature. Important conclusions about the correct lamellar-interlamellar-nucleus interactions are provided thanks to the developed model. The critical interlamellar spots with the highest delamination potentials are defined, analyzed and related to the local kinetics and microstructure.

A microstructure-based modeling approach to assess aging-sensitive mechanics of human intervertebral disc

BACKGROUND AND OBJECTIVE

The human body soft tissues are hierarchic structures interacting in a complex manner with the surrounding biochemical environment. The loss of soft tissues functionality with age leads to more vulnerability regarding to the external mechanical loadings and increases the risk of injuries. As a main example of the human body soft tissues, the intervertebral disc mechanical response evolution with age is explored. Although the age-dependence of the intervertebral disc microstructure is a well-known feature, no noticeable age effect on the disc stiffness is evidenced in the in-vitro experimental studies of the literature. So, if the disc intrinsic mechanics remains constant, how to explain the correlation of disc degeneration and disc functionality loss with age.

METHODS

A microstructure-based modeling approach was developed to assess in-silico the aging-sensitive mechanics of human intervertebral disc. The model considers the relationship between stress/volumetric macro-response and microstructure along with effective age effects acting at the lamellar and multi-lamellar scales. The stress-stretch and transversal responses of the different disc regions were computed for various age groups (13-18, 36, 58, 69 and 82 years old) and their evolution with age was studied.

RESULTS

While matching with in-vitro experimental data, the predicted stiffness was found to increase while passing from adolescent young discs to mature older discs and then to remain almost constant for the rest of life. Important age-related changes in the disc transversal behavior were also predicted affecting the flexibility of the disc, changing its volumetric behavior, and modifying its dimensions.

CONCLUSION

The developed approach was found able to bring new conclusions about age-dependent mechanical properties including regional dependency. The disc mechanics in terms of rigidity, radial and axial transversal responses were found to alter going from adolescent to middle age where the disc reaches a certain maturity. After reaching maturity, the mechanical properties undergo very slight changes until becoming almost constant with age.

Mathematical modeling and simulation of ventricular activation sequences: implications for cardiac resynchronization therapy

Next to clinical and experimental research, mathematical modeling plays a crucial role in medicine. Biomedical research takes place on many different levels, from molecules to the whole organism. Due to the complexity of biological systems, the interactions between components are often difficult or impossible to understand without the help of mathematical models. Mathematical models of cardiac electrophysiology have made a tremendous progress since the first numerical ECG simulations in the 1960s. This paper briefly reviews the development of this field and discusses some example cases where models have helped us forward, emphasizing applications that are relevant for the study of heart failure and cardiac resynchronization therapy.