Search

We describe recent preliminary studies demonstrating that computational analysis of spoken and written language can provide for highly accurate diagnostics over a wide variety of psychiatric and neurologic conditions, including psychosis, drug abuse, Parkinson's, and Alzheimer's. These results are based on the mathematical formalization of psychiatric qualitative knowledge related to the characterization of the conditions (e.g., "derailment" in psychosis) and drug effects (e.g., increased intimacy/affection with the use of the recreational drug ecstasy), as well as novel linguistic feature extraction approaches. Moreover, we show novel results using publicly available text sources suggesting that 1) it is possible to define an embedding space that allows modeling of the vast language dimension into a smaller space to map and compare different conditions, and 2) computational studies of a public personality (Ronald Reagan) can likely yield novel insights into normal aging and neurodegenerative disorders. Finally, we discuss the implications for mental health (and possibly, computer science) of a systematic application of this methodology, extended to include similar readily available behavioral data such as voice and video. [ABSTRACT FROM AUTHOR]

Background: Electromechanical wave imaging (EWI) is an entirely noninvasive, ultrasound-based imaging method capable of mapping the electromechanical activation sequence of the ventricles in vivo. Given the broad accessibility of ultrasound scanners in the clinic, the application of EWI could constitute a flexible surrogate for the 3-dimensional electrical activation.Objective: The purpose of this report is to reproduce the electromechanical wave (EW) using an anatomically realistic electromechanical model, and establish the capability of EWI to map the electrical activation sequence in vivo when pacing from different locations.Methods: EWI was performed in 1 canine during pacing from 3 different sites. A high-resolution dynamic model of coupled cardiac electromechanics of the canine heart was used to predict the experimentally recorded electromechanical wave. The simulated 3-dimensional electrical activation sequence was then compared with the experimental EW.Results: The electrical activation sequence and the EW were highly correlated for all pacing sites. The relationship between the electrical activation and the EW onset was found to be linear, with a slope of 1.01 to 1.17 for different pacing schemes and imaging angles.Conclusion: The accurate reproduction of the EW in simulations indicates that the model framework is capable of accurately representing the cardiac electromechanics and thus testing new hypotheses. The one-to-one correspondence between the electrical activation and the EW sequences indicates that EWI could be used to map the cardiac electrical activity. This opens the door for further exploration of the technique in assisting in the early detection, diagnosis, and treatment monitoring of rhythm dysfunction. [ABSTRACT FROM AUTHOR]

Reports on the foreign policy of the Georgia Republic and its relations with Russia as of March 2005. Plan of the Georgian government to integrate in the North Atlantic Treaty Organization and the European Union in 2004; Bilateral agreements suggested by Russia to Georgia in May 2004; Development of conflict between Russia and Georgia.

Lawrence Livermore National Laboratory (LLNL) has a long history of working with IBM on Blue GeneA supercomputers. Beginning in November 2001 with the joint announcement of a partnership to expand the Blue Gene research project (including Blue GeneA/L and Blue GeneA/P), the collaboration extends to this day with LLNL planning for the installation of a 96-rack Blue GeneA/Q (called Sequoia) supercomputer. As with previous machines, we envision Blue Gene/Q will be used for a wide array of applications at LLNL, ranging from meeting programmatic requirements for certification to increasing our understanding of basic physical processes. We briefly describe a representative sample of mature codes that span this application space and scale well on Blue Gene hardware. Finally, we describe advances in multi-scale whole-organ modeling of the human heart as an example of breakthrough science that will be enabled with the Blue Gene/Q architecture. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD) is believed to occur when abnormal amounts of the proteins amyloid beta and tau aggregate in the brain, resulting in a progressive loss of neuronal function. Hippocampal neurons in transgenic mice with amyloidopathy or tauopathy exhibit altered intrinsic excitability properties. We used deep hybrid modeling (DeepHM), a recently developed parameter inference technique that combines deep learning with biophysical modeling, to map experimental data recorded from hippocampal CA1 neurons in transgenic AD mice and age-matched wildtype littermate controls to the parameter space of a conductance-based CA1 model. Although mechanistic modeling and machine learning methods are by themselves powerful tools for approximating biological systems and making accurate predictions from data, when used in isolation these approaches suffer from distinct shortcomings: model and parameter uncertainty limit mechanistic modeling, whereas machine learning methods disregard the underlying biophysical mechanisms. DeepHM addresses these shortcomings by using conditional generative adversarial networks to provide an inverse mapping of data to mechanistic models that identifies the distributions of mechanistic modeling parameters coherent to the data. Here, we demonstrated that DeepHM accurately infers parameter distributions of the conductance-based model on several test cases using synthetic data generated with complex underlying parameter structures. We then used DeepHM to estimate parameter distributions corresponding to the experimental data and infer which ion channels are altered in the Alzheimer's mouse models compared to their wildtype controls at 12 and 24 months. We found that the conductances most disrupted by tauopathy, amyloidopathy, and aging are delayed rectifier potassium, transient sodium, and hyperpolarization-activated potassium, respectively., (© 2024. The Author(s).)

Predictions for physical systems often rely upon knowledge acquired from ensembles of entities, e.g. ensembles of cells in biological sciences. For qualitative and quantitative analysis, these ensembles are simulated with parametric families of mechanistic models (MMs). Two classes of methodologies, based on Bayesian inference and population of models, currently prevail in parameter estimation for physical systems. However, in Bayesian analysis, uninformative priors for MM parameters introduce undesirable bias. Here, we propose how to infer parameters within the framework of stochastic inverse problems (SIPs), also termed data-consistent inversion, wherein the prior targets only uncertainties that arise due to MM non-invertibility. To demonstrate, we introduce new methods to solve SIPs based on rejection sampling, Markov chain Monte Carlo, and generative adversarial networks (GANs). In addition, to overcome limitations of SIPs, we reformulate SIPs based on constrained optimization and present a novel GAN to solve the constrained optimization problem., Competing Interests: We declare we have no competing interests., (© 2023 The Authors.)

Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM's mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force-calcium (F-Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system., (© 2021. The Author(s).)

Background: Up to 30-50% of chronic heart failure patients who underwent cardiac resynchronization therapy (CRT) do not respond to the treatment. Therefore, patient stratification for CRT and optimization of CRT device settings remain a challenge. Objective: The main goal of our study is to develop a predictive model of CRT outcome using a combination of clinical data recorded in patients before CRT and simulations of the response to biventricular (BiV) pacing in personalized computational models of the cardiac electrophysiology. Materials and Methods: Retrospective data from 57 patients who underwent CRT device implantation was utilized. Positive response to CRT was defined by a 10% increase in the left ventricular ejection fraction in a year after implantation. For each patient, an anatomical model of the heart and torso was reconstructed from MRI and CT images and tailored to ECG recorded in the participant. The models were used to compute ventricular activation time, ECG duration and electrical dyssynchrony indices during intrinsic rhythm and BiV pacing from the sites of implanted leads. For building a predictive model of CRT response, we used clinical data recorded before CRT device implantation together with model-derived biomarkers of ventricular excitation in the left bundle branch block mode of activation and under BiV stimulation. Several Machine Learning (ML) classifiers and feature selection algorithms were tested on the hybrid dataset, and the quality of predictors was assessed using the area under receiver operating curve (ROC AUC). The classifiers on the hybrid data were compared with ML models built on clinical data only. Results: The best ML classifier utilizing a hybrid set of clinical and model-driven data demonstrated ROC AUC of 0.82, an accuracy of 0.82, sensitivity of 0.85, and specificity of 0.78, improving quality over that of ML predictors built on clinical data from much larger datasets by more than 0.1. Distance from the LV pacing site to the post-infarction zone and ventricular activation characteristics under BiV pacing were shown as the most relevant model-driven features for CRT response classification. Conclusion: Our results suggest that combination of clinical and model-driven data increases the accuracy of classification models for CRT outcomes., Competing Interests: VG was employed by IBM Research. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Khamzin, Dokuchaev, Bazhutina, Chumarnaya, Zubarev, Lyubimtseva, Lebedeva, Lebedev, Gurev and Solovyova.)

Providing therapies tailored to each patient is the vision of precision medicine, enabled by the increasing ability to capture extensive data about individual patients. In this position paper, we argue that the second enabling pillar towards this vision is the increasing power of computers and algorithms to learn, reason, and build the 'digital twin' of a patient. Computational models are boosting the capacity to draw diagnosis and prognosis, and future treatments will be tailored not only to current health status and data, but also to an accurate projection of the pathways to restore health by model predictions. The early steps of the digital twin in the area of cardiovascular medicine are reviewed in this article, together with a discussion of the challenges and opportunities ahead. We emphasize the synergies between mechanistic and statistical models in accelerating cardiovascular research and enabling the vision of precision medicine., (© The Author(s) 2020. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Computational models of the cardiovascular system and specifically heart function are currently being investigated as analytic tools to assist medical practice and clinical trials. To achieve clinical utility, models should be able to assimilate the diagnostic multi-modality data available for each patient and generate consistent representations of the underlying cardiovascular physiology. While finite element models of the heart can naturally account for patient-specific anatomies reconstructed from medical images, optimizing the many other parameters driving simulated cardiac functions is challenging due to computational complexity. With the goal of streamlining parameter adaptation, in this paper we present a novel, multifidelity strategy for model order reduction of 3-D finite element models of ventricular mechanics. Our approach is centered around well established findings on the similarity between contraction of an isolated muscle and the whole ventricle. Specifically, we demonstrate that simple linear transformations between sarcomere strain (tension) and ventricular volume (pressure) are sufficient to reproduce global pressure-volume outputs of 3-D finite element models even by a reduced model with just a single myocyte unit. We further develop a procedure for congruency training of a surrogate low-order model from multi-scale finite elements, and we construct an example of parameter optimization based on medical images. We discuss how the presented approach might be employed to process large datasets of medical images as well as databases of echocardiographic reports, paving the way towards application of heart mechanics models in the clinical practice., Competing Interests: PA, JP, JK and VG are employees of IBM. This does not alter our adherence to PLOS ONE policies on sharing data and materials. The cardiac mechanics models were implemented in the IBM Cardiac Modeling Platform. The IBM Cardiac Modeling Platform software is experimental. Readers are therefore encouraged to contact the authors if interested in using the software.

Multiscale computational models of the heart are being extensively investigated for improved assessment of drug-induced torsades de pointes (TdP) risk, a fatal side effect of many drugs. Model-derived metrics such as action potential duration and net charge carried by ionic currents ( qNet ) have been proposed as potential candidates for TdP risk stratification after being tested on small datasets. Unlike purely statistical approaches, model-derived metrics are thought to provide mechanism-based classification. In particular, qNet has been recently proposed as a surrogate metric for early afterdepolarizations (EADs), which are known to be cellular triggers of TdP. Analysis of critical model components and of the ion channels that have major impact on model-derived metrics can lead to improvements in the confidence of the prediction. In this paper, we analyze large populations of virtual drugs to systematically examine the influence of different ion channels on model-derived metrics that have been proposed for proarrhythmic risk assessment. We demonstrate via global sensitivity analysis (GSA) that model-derived metrics are most sensitive to different sets of input parameters. Similarly, important differences in sensitivity to specific channel blocks are highlighted when classifying drugs into different risk categories by either qNet or a metric directly based on simulated EADs. In particular, the higher sensitivity of qNet to the block of the late sodium channel might explain why its classification accuracy is better than that of the EAD-based metric, as shown for a small set of known drugs. Our results highlight the need for a better mechanistic interpretation of promising metrics like qNet based on a formal analysis of models. GSA should, therefore, constitute an essential component of the in silico workflow for proarrhythmic risk assessment, as an improved understanding of the structure of model-derived metrics could increase confidence in model-predicted risk., (Copyright © 2019 Parikh, Di Achille, Kozloski and Gurev.)

Patient specific models of ventricular mechanics require the optimization of their many parameters under the uncertainties associated with imaging of cardiac function. We present a strategy to reduce the complexity of parametric searches for 3-D FE models of left ventricular contraction. The study employs automatic image segmentation and analysis of an image database to gain geometric features for several classes of patients. Statistical distributions of geometric parameters are then used to design parametric studies investigating the effects of: (1) passive material properties during ventricular filling, and (2) infarct geometry on ventricular contraction in patients after a heart attack. Gaussian Process regression is used in both cases to build statistical models trained on the results of biophysical FEM simulations. The first statistical model estimates unloaded configurations based on either the intraventricular pressure or the end-diastolic fiber strain. The technique provides an alternative to the standard fixed-point iteration algorithm, which is more computationally expensive when used to unload more than 10 ventricles. The second statistical model captures the effects of varying infarct geometries on cardiac output. For training, we designed high resolution models of non-transmural infarcts including refinements of the border zone around the lesion. This study is a first effort in developing a platform combining HPC models and machine learning to investigate cardiac function in heart failure patients with the goal of assisting clinical diagnostics.

Ectopic heartbeats can trigger reentrant arrhythmias, leading to ventricular fibrillation and sudden cardiac death. Such events have been attributed to perturbed Ca2+ handling in cardiac myocytes leading to spontaneous Ca2+ release and delayed afterdepolarizations (DADs). However, the ways in which perturbation of specific molecular mechanisms alters the probability of ectopic beats is not understood. We present a multiscale model of cardiac tissue incorporating a biophysically detailed three-dimensional model of the ventricular myocyte. This model reproduces realistic Ca2+ waves and DADs driven by stochastic Ca2+ release channel (RyR) gating and is used to study mechanisms of DAD variability. In agreement with previous experimental and modeling studies, key factors influencing the distribution of DAD amplitude and timing include cytosolic and sarcoplasmic reticulum Ca2+ concentrations, inwardly rectifying potassium current (IK1) density, and gap junction conductance. The cardiac tissue model is used to investigate how random RyR gating gives rise to probabilistic triggered activity in a one-dimensional myocyte tissue model. A novel spatial-average filtering method for estimating the probability of extreme (i.e. rare, high-amplitude) stochastic events from a limited set of spontaneous Ca2+ release profiles is presented. These events occur when randomly organized clusters of cells exhibit synchronized, high amplitude Ca2+ release flux. It is shown how reduced IK1 density and gap junction coupling, as observed in heart failure, increase the probability of extreme DADs by multiple orders of magnitude. This method enables prediction of arrhythmia likelihood and its modulation by alterations of other cellular mechanisms.

While pre-clinical Torsades de Pointes (TdP) risk classifiers had initially been based on drug-induced block of hERG potassium channels, it is now well established that improved risk prediction can be achieved by considering block of non-hERG ion channels. The current multi-channel TdP classifiers can be categorized into two classes. First, the classifiers that take as input the values of drug-induced block of ion channels (direct features). Second, the classifiers that are built on features extracted from output of the drug-induced multi-channel blockage simulations in the in-silico models (derived features). The classifiers built on derived features have thus far not consistently provided increased prediction accuracies, and hence casts doubt on the value of such approaches given the cost of including biophysical detail. Here, we propose a new two-step method for TdP risk classification, referred to as Multi-Channel Blockage at Early After Depolarization (MCB@EAD). In the first step, we classified the compound that produced insufficient hERG block as non-torsadogenic. In the second step, the role of non-hERG channels to modulate TdP risk are considered by constructing classifiers based on direct or derived features at critical hERG block concentrations that generates EADs in the computational cardiac cell models. MCB@EAD provides comparable or superior TdP risk classification of the drugs from the direct features in tests against published methods. TdP risk for the drugs highly correlated to the propensity to generate EADs in the model. However, the derived features of the biophysical models did not improve the predictive capability for TdP risk assessment.

Models of cardiac mechanics are increasingly used to investigate cardiac physiology. These models are characterized by a high level of complexity, including the particular anisotropic material properties of biological tissue and the actively contracting material. A large number of independent simulation codes have been developed, but a consistent way of verifying the accuracy and replicability of simulations is lacking. To aid in the verification of current and future cardiac mechanics solvers, this study provides three benchmark problems for cardiac mechanics. These benchmark problems test the ability to accurately simulate pressure-type forces that depend on the deformed objects geometry, anisotropic and spatially varying material properties similar to those seen in the left ventricle and active contractile forces. The benchmark was solved by 11 different groups to generate consensus solutions, with typical differences in higher-resolution solutions at approximately 0.5%, and consistent results between linear, quadratic and cubic finite elements as well as different approaches to simulating incompressible materials. Online tools and solutions are made available to allow these tests to be effectively used in verification of future cardiac mechanics software.

Modeling of the heart ventricles is one of the most challenging tasks in soft tissue mechanics because cardiac tissue is a strongly anisotropic incompressible material with an active component of stress. In most current approaches with active force, the number of degrees of freedom (DOF) is limited by the direct method of solution of linear systems of equations. We develop a new approach for high-resolution heart models with large numbers of DOF by: (1) developing a hex-dominant finite element mixed formulation and (2) developing a Krylov subspace iterative method that is able to solve the system of linearized equations for saddle-point problems with active stress. In our approach, passive cardiac tissue is modeled as a hyperelastic, incompressible material with orthotropic properties, and mixed pressure-displacement finite elements are used to enforce incompressibility. Active stress is generated by a model with force dependence on length and velocity of muscle shortening. The ventricles are coupled to a lumped circulatory model. For efficient solution of linear systems, we use Flexible GMRES with a nonlinear preconditioner based on block matrix decomposition involving the Schur complement. Three methods for approximating the inverse of the Schur complement are evaluated: inverse of the pressure mass matrix; least squares commutators; and sparse approximate inverse. The sub-matrix corresponding to the displacement variables is preconditioned by a V-cycle of hybrid geometric-algebraic multigrid followed by correction with several iterations of GMRES preconditioned by sparse approximate inverse. The overall solver is demonstrated on a high-resolution two ventricle mesh based on a human anatomy with roughly 130 K elements and 1.7 M displacement DOF. Effectiveness of the numerical method for active contraction is shown. To the best of our knowledge, this solver is the first to efficiently model ventricular contraction using an iterative linear solver for the mesh size demonstrated and therefore opens the possibility for future very high-resolution models. In addition, several relatively simple benchmark problems are designed for a verification exercise to show that the solver is functioning properly and correctly solves the underlying mathematical model. Here, the output of the newly designed solver is compared to that of the mechanics component of Chaste ('Cancer, Heart and Soft Tissue Environment'). These benchmark tests may be used by other researchers to verify their newly developed methods and codes.

Background: Cardiac resynchronization therapy (CRT) has been demonstrated to lead to restoration of oxygen consumption homogeneity throughout the left ventricle (LV), which is important for long-term reverse remodeling of the ventricles. However, research has focused exclusively on identifying the LV pacing sites that led to acute hemodynamic improvements. It remains unclear whether there exist LV pacing sites that could both improve the hemodynamics and result in ATP consumption homogeneity throughout the LV, thus maximizing both CRT short-term and long-term benefits., Objective: The purpose of this study was to demonstrate the feasibility of optimizing CRT pacing locations to achieve maximal improvement in both ATPCTHI (an ATP consumption heterogeneity index) and stroke work., Methods: We used an magnetic resonance image-based electromechanical model of the dyssynchronous heart failure (DHF) canine ventricles. ATPCTHI and stroke work improvement were determined for each of 34 CRT pacing sites evenly spaced over the LV epicardium., Results: Results demonstrated the feasibility of determining the optimal LV pacing site that achieves simultaneous maximum improvements in ATPCTHI and stroke work. The optimal LV CRT pacing sites in the DHF canine ventricles were located midway between apex and base. The improvement in ATPCTHI decreased more rapidly with the distance from the optimal sites compared to stroke work improvement. CRT from the optimal sites homogenized ATP consumption by increasing septal ATP consumption and decreasing that of the lateral wall., Conclusion: Simulation results using a canine heart failure model demonstrated that CRT can be optimized to achieve improvements in both ATPCTHI and stroke work., (Copyright © 2014 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Background: The acute response to cardiac resynchronization therapy (CRT) has been shown to be due to 3 mechanisms: resynchronization of ventricular contraction, efficient preloading of the ventricles by a properly timed atrial contraction, and mitral regurgitation reduction. However, the contribution of each of the 3 mechanisms to the acute response to CRT, specifically stroke work improvement, has not been quantified., Objective: To use a magnetic resonance image-based anatomically accurate 3-dimensional model of failing canine ventricular electromechanics to quantify the contribution of each of the 3 mechanisms to stroke work improvement and identify the predominant mechanisms., Methods: An MRI-based electromechanical model of the failing canine ventricles assembled previously by our group was further developed and modified. Three different protocols were used to dissect the contribution of each of the 3 mechanisms to stroke work improvement., Results: Resynchronization of ventricular contraction did not lead to a significant stroke work improvement. Efficient preloading of the ventricles by a properly timed atrial contraction was the predominant mechanism underlying stroke work improvement. Stroke work improvement peaked at an intermediate atrioventricular delay, as it allowed ventricular filling by atrial contraction to occur at a low diastolic left ventricular pressure but also provided adequate time for ventricular filling before ventricular contraction. Reduction of mitral regurgitation by CRT led to stroke work worsening instead of improvement., Conclusion: Efficient preloading of the ventricles by a properly timed atrial contraction is responsible for a significant stroke work improvement in the acute CRT response., (Copyright © 2013 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, -60 mV and -10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of -60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of -10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of -10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.

We have developed the capability to rapidly simulate cardiac electrophysiological phenomena in a human heart discretised at a resolution comparable with the length of a cardiac myocyte. Previous scientific investigation has generally invoked simplified geometries or coarse-resolution hearts, with simulation duration limited to 10s of heartbeats. Using state-of-the-art high-performance computing techniques coupled with one of the most powerful computers available (the 20 PFlop/s IBM BlueGene/Q at Lawrence Livermore National Laboratory), high-resolution simulation of the human heart can now be carried out over 1200 times faster compared with published results in the field. We demonstrate the utility of this capability by simulating, for the first time, the formation of transmural re-entrant waves in a 3D human heart. Such wave patterns are thought to underlie Torsades de Pointes, an arrhythmia that indicates a high risk of sudden cardiac death. Our new simulation capability has the potential to impact a multitude of applications in medicine, pharmaceuticals and implantable devices.

Left ventricular-assist devices (LVADs) are used to supply blood to the body of patients with heart failure. Pressure unloading is greater for counter-pulsating LVADs than for continuous LVADs. However, several clinical trials have demonstrated that myocardial recovery is similar for both types of LVAD. This study examined the contractile energy consumption of the myocardium with continuous and counter-pulsating LVAD support to ascertain the effect of the different LVADs on myocardial recovery. We used a three-dimensional electromechanical model of canine ventricles, with models of the circulatory system and an LVAD. We compared the left ventricular peak pressure (LVPP) and contractile ATP consumption between pulsatile and continuous LVADs. With the continuous and counter-pulsating LVAD, the LVPP decreased to 46 and 10%, respectively, and contractile ATP consumption decreased to 60 and 50%. The small difference between the contractile ATP consumption of these two types of LVAD may explain the comparable effects of the two types on myocardial recovery.

A three-dimensional (3D) finite element electromechanical model of the heart is employed in simulations of seismocardiograms (SCGs). To simulate SCGs, a previously developed 3D model of ventricular contraction is extended by adding the mechanical interaction of the heart with the chest and internal organs. The proposed model reproduces the major peaks of seismocardiographic signals during the phases of the cardiac cycle. Results indicate that SCGs record the pressure of the heart acting on the ribs. In addition, the model reveals that the rotation of the rib with respect to the heart has a minor effect on seismocardiographic signal morphology during the respiratory cycle. SCGs are obtained from 24 human volunteers and their morphology is analyzed. Experimental results demonstrate that the peak of the maximum acceleration of blood in the aorta occurs at the same time as the global minimum of the SCG. It is confirmed that the first SCG peak after the electrocardiogram R-wave corresponds to aortic valve opening, as determined from the impedance cardiogram (p = 0.92). The simulation results reveal that the SCG peaks corresponding to aortic valve opening and the maximum acceleration of blood in the aorta result from ventricular contraction in the longitudinal direction of the ventricles and a decrease in the dimensions of the ventricles due to the ejection of blood, respectively.

Computational modeling has traditionally played an important role in dissecting the mechanisms for cardiac dysfunction. Ventricular electromechanical models, likely the most sophisticated virtual organs to date, integrate detailed information across the spatial scales of cardiac electrophysiology and mechanics and are capable of capturing the emergent behavior and the interaction between electrical activation and mechanical contraction of the heart. The goal of this review is to provide an overview of the latest advancements in multiscale electromechanical modeling of the ventricles. We first detail the general framework of multiscale ventricular electromechanical modeling and describe the state of the art in computational techniques and experimental validation approaches. The powerful utility of ventricular electromechanical models in providing a better understanding of cardiac function is then demonstrated by reviewing the latest insights obtained by these models, focusing primarily on the mechanisms by which mechanoelectric coupling contributes to ventricular arrythmogenesis, the relationship between electrical activation and mechanical contraction in the normal heart, and the mechanisms of mechanical dyssynchrony and resynchronization in the failing heart. Computational modeling of cardiac electromechanics will continue to complement basic science research and clinical cardiology and holds promise to become an important clinical tool aiding the diagnosis and treatment of cardiac disease.

Current multi-scale computational models of ventricular electromechanics describe the full process of cardiac contraction on both the micro- and macro- scales including: the depolarization of cardiac cells, the release of calcium from intracellular stores, tension generation by cardiac myofilaments, and mechanical contraction of the whole heart. Such models are used to reveal basic mechanisms of cardiac contraction as well as the mechanisms of cardiac dysfunction in disease conditions. In this paper, we present a methodology to construct finite element electromechanical models of ventricular contraction with anatomically accurate ventricular geometry based on magnetic resonance and diffusion tensor magnetic resonance imaging of the heart. The electromechanical model couples detailed representations of the cardiac cell membrane, cardiac myofilament dynamics, electrical impulse propagation, ventricular contraction, and circulation to simulate the electrical and mechanical activity of the ventricles. The utility of the model is demonstrated in an example simulation of contraction during sinus rhythm using a model of the normal canine ventricles.

One of the most important components of mechanoelectric coupling is stretch-activated channels, sarcolemmal channels that open upon mechanical stimuli. Uncovering the mechanisms by which stretch-activated channels contribute to ventricular arrhythmogenesis under a variety of pathologic conditions is hampered by the lack of experimental methodologies that can record the 3-dimensional electromechanical activity simultaneously at high spatiotemporal resolution. Computer modeling provides such an opportunity. The goal of this review is to illustrate the utility of sophisticated, physiologically realistic, whole heart computer simulations in determining the role of mechanoelectric coupling in ventricular arrhythmogenesis. We first present the various ways by which stretch-activated channels have been modeled and demonstrate how these channels affect cardiac electrophysiologic properties. Next, we use an electrophysiologic model of the rabbit ventricles to understand how so-called commotio cordis, the mechanical impact to the precordial region of the heart, can initiate ventricular tachycardia via the recruitment of stretch-activated channels. Using the same model, we also provide mechanistic insight into the termination of arrhythmias by precordial thump under normal and globally ischemic conditions. Lastly, we use a novel anatomically realistic dynamic 3-dimensional coupled electromechanical model of the rabbit ventricles to gain insight into the role of electromechanical dysfunction in arrhythmogenesis during acute regional ischemia., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Rationale: Although ventricular premature beats (VPBs) during acute regional ischemia have been linked to mechanical stretch of ischemic tissue, whether and how ischemia-induced mechanical dysfunction can induce VPBs and facilitate their degradation into reentrant arrhythmias has not been yet addressed., Objective: This study used a novel multiscale electromechanical model of the rabbit ventricles to investigate the origin of and the substrate for spontaneous arrhythmias arising from ischemia-induced electrophysiological and mechanical changes., Methods and Results: Two stages of ischemia were simulated. Dynamic mechanoelectrical feedback was modeled as spatially and temporally nonuniform membrane currents through mechanosensitive channels, the conductances of which depended on local strain rate. Our results reveal that both strains and strain rates were significantly larger in the central ischemic zone than in the border zone. However, in both ischemia stages, a VPB originated from the ischemic border in the left ventricular apical endocardium because of mechanically induced suprathreshold depolarizations. It then traveled fully intramurally until emerging from the ischemic border on the anterior epicardium. Reentry was formed only in the advanced ischemia stage as the result of a widened temporal excitable gap. Mechanically induced delayed afterdepolarization-like events contributed to the formation of reentry by further decreasing the already reduced-by-hyperkalemia local excitability, causing extended conduction block lines and slowed conduction in the ischemic region., Conclusions: Mechanically induced membrane depolarizations in the ischemic region are the mechanism by which mechanical activity contributes to both the origin of and substrate for spontaneous arrhythmias under the conditions of acute regional ischemia.

The simulation of cardiac electrical function is an example of a successful integrative multiscale modelling approach that is directly relevant to human disease. Today we stand at the threshold of a new era, in which anatomically detailed, tomographically reconstructed models are being developed that integrate from the ion channel to the electromechanical interactions in the intact heart. Such models hold high promise for interpretation of clinical and physiological measurements, for improving the basic understanding of the mechanisms of dysfunction in disease, such as arrhythmias, myocardial ischaemia and heart failure, and for the development and performance optimization of medical devices. The goal of this article is to present an overview of current state-of-art advances towards predictive computational modelling of the heart as developed recently by the authors of this article. We first outline the methodology for constructing electrophysiological models of the heart. We then provide three examples that demonstrate the use of these models, focusing specifically on the mechanisms for arrhythmogenesis and defibrillation in the heart. These include: (1) uncovering the role of ventricular structure in defibrillation; (2) examining the contribution of Purkinje fibres to the failure of the shock; and (3) using magnetic resonance imaging reconstructed heart models to investigate the re-entrant circuits formed in the presence of an infarct scar.

Cardiovascular diseases are leading causes of death, reduce life quality and consume almost half a trillion dollars in healthcare expenses in the USA alone. Cardiac modeling and simulation technologies hold promise as important tools to improve cardiac care and are already in use to elucidate the fundamental mechanisms of cardiac physiology and pathophysiology. However, the emphasis has been on simulating average or exemplar cases. This report describes two classes of cardiac modeling efforts. First, electrophysiological models of channelopathies simulate the altered electrical activity that is thought to promote arrhythmias. Second, electromechanical models attempt to capture both the electrophysiological and mechanical aspects of heart function. One goal of the community is to develop models with sufficient patient customization to assist in personalized treatment planning. Some model aspects can be customized with existing data collection techniques to more closely represent individual patients while other model aspects will likely remain based on generic data. Despite important challenges, cardiac models hold promise to be important enablers of personalized medicine.

We introduce and compare three different modalities to study seismocardiogram (SCG) and its correlation with cardiac events. We used an accelerometer attached to the subject sternum to get a reference measure. Cardiac events were then approximately identified using echocardiography. As an alternative approximation, we used consecutive Cine-MRI images of the heart to capture cardiac movements and compared them with the experimental SCG. We also employed an anatomically accurate, finite element base electromechanical model with geometry built completely from DT-MRI to simulate a portion of the cardiac cycle as observed in the SCG signal. The preliminary results demonstrate the usability of these newly proposed methods to investigate the mechanism of SCG waves and also demonstrate the usability of echocardiograph in interpretation of these results in terms of correlating them to underlying cardiac cycle events.

Experimental and clinical studies have shown that ventricular dilatation is associated with increased arrhythmogenesis and elevated defibrillation threshold; however, the underlying mechanisms remain poorly understood. The goal of the present study was to test the hypothesis that (1) stretch-activated channel (SAC) recruitment and (2) geometrical deformations in organ shape and fiber architecture lead to increased arrhythmogenesis by electric shocks following acute ventricular dilatation. To elucidate the contribution of these two factors, the study employed, for the first time, a combined electro-mechanical simulation approach. Acute dilatation was simulated in a model of rabbit ventricular mechanics by raising the LV end-diastolic pressure from 0.6 (control) to 4.2 kPa (dilated). The output of the mechanics model was used in the electrophysiological model. Vulnerability to shocks was examined in the control, the dilated ventricles, and in the dilated ventricles that also incorporated currents through SAC as a function of local strain, by constructing vulnerability grids. Results showed that dilatation-induced deformation alone decreased upper limit of vulnerability (ULV) slightly and did not result in increased vulnerability. With SAC recruitment in the dilated ventricles, the number of shock-induced arrhythmia episodes increased by 37% (from 41 to 56) and the lower limit of vulnerability (LLV) decreased from 9 to 7 V/cm, while ULV did not change. The heterogeneous activation of SAC caused by the heterogeneous fiber strain in the ventricular walls was the main reason for increased vulnerability to electric shocks since it caused dispersion of electrophysiological properties in the tissue, resulting in postshock unidirectional block and establishment of reentry.

Ventricular dilatation increases the defibrillation threshold (DFT). In order to elucidate the mechanisms responsible for this increase, the present article investigates changes in the postshock behavior of the myocardium upon stretch. A two-dimensional electro-mechanical model of cardiac tissue incorporating heterogeneous fiber orientation was used to explore the effect of sustained stretch on postshock behavior via (a) recruitment of mechanosensitive channels (MSC) and (b) tissue deformation and concomitant changes in tissue conductivities. Recruitment of MSC had no influence on vulnerability to electric shocks as compared to control, but increased the complexity of postshock VF patterns. Stretch-induced deformation and changes in tissue conductivities resulted in a decrease in vulnerability to electric shocks.