Low-energy defibrillation using a base-apex epicardial electrode.

BACKGROUND: Current implantable cardioverter defibrillators (ICDs) require electric conduction with high voltage and high energy, which can impair cardiac function and induce another malignant arrhythmia. As a result, there has been a demand for an ICD that can effectively operate with lower energy to mitigate the risks of a strong electric shock. METHODS: A pair of sheet-shaped electrodes covering the heart were analyzed in three configurations (top-bottom, left-right, and front-back) using a heart simulator. We also varied the distance between the two electrodes (clearance) to identify the electrode shape with the lowest defibrillation threshold (DFT). We also investigated the ICD shock waveform, shock direction, and the effect of the backside insulator of the electrode. RESULTS: The DFT was high when the clearance was too small and the DFT was high even when the clearance was too large, suggesting that an optimal value clearance. The top-bottom electrodes with optimal clearance showed the lowest DFT when the biphasic shocks set the top electrode to a high potential first and then the bottom electrode was set to a high potential. An interval between a first shock waveform and a second shock waveform should be provided for low-energy defibrillation. Because the insulator prevents unnecessary current flow to the backside, the DFT of the electrodes with insulators is less than those without insulators. CONCLUSION: Painless defibrillation using sheet-shaped electrodes on the epicardium is predicated on the basis of results using a heart simulator.

Transition mechanisms from atrial flutter to atrial fibrillation during anti-tachycardia pacing therapy.

BACKGROUND: Atrial anti-tachycardia pacing (aATP) has been shown to be effective for the termination of atrial tachyarrhythmias, but its success rate is still not high enough. OBJECTIVE: The main objective of this study was to investigate the mechanisms of atrial flutter (AFL) termination by aATP and the transition from AFL to atrial fibrillation (AF) during aATP. METHODS: We developed a multi-scale model of the human atrium based on magnetic resonance images and examined the atrial electrophysiology of AFL during aATP with a ramp protocol. RESULTS: In successful cases of aATP, paced excitation entered the excitable gap and collided with the leading edge of the reentrant wave front. Furthermore, the excitation propagating in the opposite direction collided with the trailing edge of the reentrant wave to terminate AFL. The second collision was made possible by the distribution of the wave propagation velocity in the atria. The detailed analysis revealed that the slowing of propagation velocity occurred at the exit of the sub-Eustachian isthmus, probably due to source-sink mismatch. During the transition from AFL to AF, the excitation collided with the refractory zone of the preceding wave and broke into multiple wave fronts to induce AF. A similar observation was made for the transition from AF to sinus rhythm. In both cases, the complex anatomy of the atria played an essential role. CONCLUSION: The complex anatomy of atria plays an essential role in the maintenance of stable AFL and its termination by aATP, which were revealed by the realistic three-dimensional simulation model.

An application of a patient-specific cardiac simulator for the prediction of outcomes after mitral valve replacement: a pilot study.

Despite advancements in preoperative prediction of patient outcomes, determination of the most appropriate surgical treatments for patients with severely impaired cardiac function remains a challenge. "UT-Heart" is a multi-scale, multi-physics heart simulator, which can be used to assess the effects of treatment without imposing any burden on the patients. This retrospective study aimed to assess whether UT-Heart can function as a tool that aids decision making for performing mitral valve replacements (MVR) in patients with severe mitral regurgitation (MR) and impaired left ventricular (LV) function. We used preoperative clinical data to create a patient-specific heart model using UT-Heart for a patient who had dilated cardiomyopathy with severe MR. After confirming that this heart model reproduced the preoperative state of the patient, we performed an in silico MVR operation without changing any parameters, such as the end-diastolic volume of the left ventricle, systemic vascular resistance, and the number of myocardiocytes. Among the functional changes introduced by in silico surgery, we found two indices, forward flow and the mechanical efficiency of the work done to the systemic circulation, which may relate positively to the favorable outcome observed in the real world. Thus, multi-scale, multi-physics heart simulators can reproduce the pathophysiology of MR with impaired LV function. By performing in silico MVR and examining the resultant functional changes, we identified two indices, whose usefulness should be tested in future studies.

Patient-specific heart simulation can identify non-responders to cardiac resynchronization therapy.

To identify non-responders to cardiac resynchronization therapy (CRT), various biomarkers have been proposed, but these attempts have not been successful to date. We tested the clinical applicability of computer simulation of CRT for the identification of non-responders. We used the multi-scale heart simulator "UT-Heart," which can reproduce the electrophysiology and mechanics of the heart based on a molecular model of the excitation-contraction mechanism. Patient-specific heart models were created for eight heart failure patients who were treated with CRT, based on the clinical data recorded before treatment. Using these heart models, bi-ventricular pacing simulations were performed at multiple pacing sites adopted in clinical practice. Improvement in pumping function measured by the relative change of maximum positive derivative of left ventricular pressure (%ΔdP/dtmax) was compared with the clinical outcome. The operators of the simulation were blinded to the clinical outcome. In six patients, the relative reduction in end-systolic volume exceeded 15% in the follow-up echocardiogram at 3 months (responders) and the remaining two patients were judged as non-responders. The simulated %ΔdP/dtmax at the best lead position could identify responders and non-responders successfully. With further refinement of the model, patient-specific simulation could be a useful tool for identifying non-responders to CRT.

Longitudinal dissociation and transition in thickness of the His-Purkinje system cause various QRS waveforms of surface ECG under His bundle pacing: A simulation study based on clinical observations.

AIMS: His bundle pacing (HBP) is a feasible and reliable alternative to conventional right ventricular pacing (RVP), but associated ECG (electrocardiogram) changes have not been well-studied. This study aimed to determine the mechanisms underlying ECG changes associated with HBP using patient-specific multiscale heart simulations. METHODS: ECGs were recorded in two patients who were treated by HBP under a native rhythm and HBP at high and low voltages. We created patient-specific multiscale simulation heart models of these patients and performed ECG simulation under these conditions. Using these results and detailed information on the electrical field around the pacing lead, we investigated mechanisms underlying the observed ECG changes. RESULTS: Heart simulations successfully reproduced ECGs under a native rhythm for both cases. In case 1, nonselective HBP produced a left bundle branch (LBB) block pattern, which was reproduced as a selective right bundle branch (RBB) pacing. However, in case 2, ECG under nonselective HBP showed an RBB block pattern, which could not be reproduced by the commonly used framework. Findings on the electrical field and anatomy of the His bundle and its branches suggested that longitudinal dissociation of the His bundle and transition of thickness in the stem of the LBB caused a conduction delay in the RBB to produce these ECG changes in this patient. CONCLUSION: Variations in the anatomy of the His bundle and its branches may underlie the diverse ECG responses to HBP. These variations should be taken into account when performing this therapy.

Clinical and pharmacological application of multiscale multiphysics heart simulator, UT-Heart.

A heart simulator, UT-Heart, is a finite element model of the human heart that can reproduce all the fundamental activities of the working heart, including propagation of excitation, contraction, and relaxation and generation of blood pressure and blood flow, based on the molecular aspects of the cardiac electrophysiology and excitation-contraction coupling. In this paper, we present a brief review of the practical use of UT-Heart. As an example, we focus on its application for predicting the effect of cardiac resynchronization therapy (CRT) and evaluating the proarrhythmic risk of drugs. Patient-specific, multiscale heart simulation successfully predicted the response to CRT by reproducing the complex pathophysiology of the heart. A proarrhythmic risk assessment system combining in vitro channel assays and in silico simulation of cardiac electrophysiology using UT-Heart successfully predicted druginduced arrhythmogenic risk. The assessment system was found to be reliable and efficient. We also developed a comprehensive hazard map on the various combinations of ion channel inhibitors. This in silico electrocardiogram database (now freely available at http://ut-heart.com/) can facilitate proarrhythmic risk assessment without the need to perform computationally expensive heart simulation. Based on these results, we conclude that the heart simulator, UT-Heart, could be a useful tool in clinical medicine and drug discovery.

Multi-scale, tailor-made heart simulation can predict the effect of cardiac resynchronization therapy.

BACKGROUND: The currently proposed criteria for identifying patients who would benefit from cardiac resynchronization therapy (CRT) still need to be optimized. A multi-scale heart simulation capable of reproducing the electrophysiology and mechanics of a beating heart may help resolve this problem. The objective of this retrospective study was to test the capability of patient-specific simulation models to reproduce the response to CRT by applying the latest multi-scale heart simulation technology. METHODS AND RESULTS: We created patient-specific heart models with realistic three-dimensional morphology based on the clinical data recorded before treatment in nine patients with heart failure and conduction block treated by biventricular pacing. Each model was tailored to reproduce the surface electrocardiogram and hemodynamics of each patient in formats similar to those used in clinical practice, including electrocardiography (ECG), echocardiography, and hemodynamic measurements. We then performed CRT simulation on each heart model according to the actual pacing protocol and compared the results with the clinical data. CRT simulation improved the ECG index and diminished wall motion dyssynchrony in each patient. These results, however, did not correlate with the actual response. The best correlation was obtained between the maximum value of the time derivative of ventricular pressure (dP/dtmax) and the clinically observed improvement in the ejection fraction (EF) (r=0.94, p<0.01). CONCLUSIONS: By integrating the complex pathophysiology of the heart, patient-specific, multi-scale heart simulation could successfully reproduce the response to CRT. With further verification, this technique could be a useful tool in clinical decision making.

Distinct functional roles of cardiac mitochondrial subpopulations revealed by a 3D simulation model.

Experimental characterization of two cardiac mitochondrial subpopulations, namely, subsarcolemmal mitochondria (SSM) and interfibrillar mitochondria (IFM), has been hampered by technical difficulties, and an alternative approach is eagerly awaited. We previously developed a three-dimensional computational cardiomyocyte model that integrates electrophysiology, metabolism, and mechanics with subcellular structure. In this study, we further developed our model to include intracellular oxygen diffusion, and determined whether mitochondrial localization or intrinsic properties cause functional variations. For this purpose, we created two models: one with equal SSM and IFM properties and one with IFM having higher activity levels. Using these two models to compare the SSM and IFM responses of [Ca(2+)], tricarboxylic acid cycle activity, [NADH], and mitochondrial inner membrane potential to abrupt changes in pacing frequency (0.25-2 Hz), we found that the reported functional differences between these subpopulations appear to be mostly related to local [Ca(2+)] heterogeneity, and variations in intrinsic properties only serve to augment these differences. We also examined the effect of hypoxia on mitochondrial function. Under normoxic conditions, intracellular oxygen is much higher throughout the cell than the half-saturation concentration for oxidative phosphorylation. However, under limited oxygen supply, oxygen is mostly exhausted in SSM, leaving the core region in an anoxic condition. Reflecting this heterogeneous oxygen environment, the inner membrane potential continues to decrease in IFM, whereas it is maintained to nearly normal levels in SSM, thereby ensuring ATP supply to this region. Our simulation results provide clues to understanding the origin of functional variations in two cardiac mitochondrial subpopulations and their differential roles in maintaining cardiomyocyte function as a whole.

Mitochondrial colocalization with Ca2+ release sites is crucial to cardiac metabolism.

In cardiomyocyte subcellular structures, colocalization of mitochondria with Ca2+ release sites is implicated in regulation of cardiac energetics by facilitating Ca2+ influx into mitochondria to modulate the tricarboxylic acid (TCA) cycle. However, current experimental techniques limit detailed examination of this regulatory mechanism. Earlier, we developed a three-dimensional (3D) finite-element cardiomyocyte model featuring a subcellular structure that integrates excitation-contraction coupling and energy metabolism. Here, using this model, we examined the influence of distance between mitochondria and Ca2+ release sites by comparing a normal (50-nm) distance model and a large (200-nm) distance model (LD). The influence of distance was minimal under a low pacing rate (0.25 Hz), but under a higher pacing rate (2 Hz), lower levels of mitochondrial Ca2+ and NADH, elevated phosphate, and suppressed force generation became apparent in the LD model. Such differences became greater when functional impairments (reduced TCA cycle activity, uncoupling effect, and failing excitation-contraction coupling) were additionally imposed. We concluded that juxtaposition of the mitochondria and the Ca2+ release sites is crucial for rapid signal transmission to maintain cardiac-energy balance. The idealized 3D model of cardiac excitation-contraction and metabolism is a powerful tool to study cardiac energetics.

Effects of an hERG activator, ICA-105574, on electrophysiological properties of canine hearts.

In short QT syndrome, inherited gain-of-function mutations in the human ether a-gogo-related gene (hERG) K(+) channel have been associated with development of fatal arrhythmias. This implies that drugs that activate hERG as a side effect may likewise pose significant arrhythmia risk. hERG activators have been found to have diverse mechanisms of activation, which may reflect their distinct binding sites. Recently, the new hERG activator ICA-105574 was introduced, which disables inactivation of the hERG channel with very high potency. We explored characteristics of this new drug in several experimental models. Patch clamp experiments were used to verify activation of hERG channels by ICA-105574 in human embryonic kidney cells stably-expressing hERG channels. ICA-105574 significantly shortened QT and QTc intervals and monophasic action potential duration (MAP(90)) in Langendorff-perfused guinea-pig hearts. We also administered ICA-105574 to anesthetized dogs while recording ECG and drug plasma concentrations. ICA-105574 (10 mg/kg) significantly shortened QT and QTc intervals, with a free plasma concentration of approximately 1.7 µM at the point of maximal effect. Our data showed that unbound ICA-105574 caused QT shortening in dogs at concentrations comparable to the half maximal effective concentration (EC(50), 0.42 µM) of hERG activation in the patch clamp studies.

Patient specific simulation of body surface ECG using the finite element method.

BACKGROUND: Recent studies, supported by advances in computer science, have successfully simulated the excitation and repolarization processes of the heart, based on detailed cell models of electrophysiology and implemented with realistic morphology. METHODS: In this study, we extend these approaches to simulate the body surface electrocardiogram (ECG) of specific individuals. Patient-specific finite element models of the heart and torso are created for four patients with various heart diseases, based on clinical data including computer tomography, while the parallel multi-grid method is used to solve the dynamic bi-domain problem. Personalization procedures include demarcation of nonexcitable tissue, allocation of the failing myocyte model of electrophysiology, and modification of the excitation sequence. In particular, the adjustment of QRS morphology requires iterative computations, facilitated by the simultaneous visualization of the propagation of excitation in the heart, average QRS vector in the torso, and 12-lead ECG. RESULTS: In all four cases we obtained reasonable agreement between the simulated and actual ECGs. Furthermore, we also simulated the ECGs of three of the patients under bi-ventricular pacing, and once again successfully reproduced the actual ECG morphologies. Since no further adjustments were made to the heart models in the pacing simulations, the good agreement provides strong support for the validity of the models. CONCLUSIONS: These results not only help us understand the cellular basis of the body surface ECG, but also open the possibility of heart simulation for clinical applications.

UT-Heart: A Finite Element Model Designed for the Multiscale and Multiphysics Integration of our Knowledge on the Human Heart.

To fully understand the health and pathology of the heart, it is necessary to integrate knowledge accumulated at molecular, cellular, tissue, and organ levels. However, it is difficult to comprehend the complex interactions occurring among the building blocks of biological systems across these scales. Recent advances in computational science supported by innovative high-performance computer hardware make it possible to develop a multiscale multiphysics model simulating the heart, in which the behavior of each cell model is controlled by molecular mechanisms and the cell models themselves are arranged to reproduce elaborate tissue structures. Such a simulator could be used as a tool not only in basic science but also in clinical settings. Here, we describe a multiscale multiphysics heart simulator, UT-Heart, which uses unique technologies to realize the abovementioned features. As examples of its applications, models for cardiac resynchronization therapy and surgery for congenital heart disease will be also shown.

A thermodynamically consistent monte carlo cross-bridge model with a trapping mechanism reveals the role of stretch activation in heart pumping.

Changes in intracellular calcium concentrations regulate heart beats. However, the decline in the left ventricular pressure during early diastole is much sharper than that of the Ca2+ transient, resulting in a rapid supply of blood to the left ventricle during the diastole. At the tissue level, cardiac muscles have a distinct characteristic, known as stretch activation, similar to the function of insect flight muscles. Stretch activation, which is a delayed increase in force following a rapid muscle length increase, has been thought to be related to autonomous control in these muscles. In this numerical simulation study, we introduced a molecular mechanism of stretch activation and investigated the role of this mechanism in the pumping function of the heart, using the previously developed coupling multiple-step active stiffness integration scheme for a Monte Carlo (MC) cross-bridge model and a bi-ventricular finite element model. In the MC cross-bridge model, we introduced a mechanism for trapping the myosin molecule in its post-power stroke state. We then determined the rate constants of transitions for trapping and escaping in a thermodynamically consistent manner. Based on our numerical analysis, we draw the following conclusions regarding the stretch activation mechanism: (i) the delayed force becomes larger than the original isometric force because the population of trapped myosin molecules and their average force increase after stretching; (ii) the delayed force has a duration of more than a few seconds owing to a fairly small rate constant of escape from the trapped state. For the role of stretch activation in heart pumping, we draw the following conclusions: (iii) for the regions in which the contraction force decreases earlier than the neighboring region in the end-systole phase, the trapped myosin molecules prevent further lengthening of the myocytes, which then prevents further shortening of neighboring myocytes; (iv) as a result, the contraction forces are sustained longer, resulting in a larger blood ejection, and their degeneration is synchronized.

A three-dimensional simulation model of cardiomyocyte integrating excitation-contraction coupling and metabolism.

Recent studies have revealed that Ca(2+) not only regulates the contraction of cardiomyocytes, but can also function as a signaling agent to stimulate ATP production by the mitochondria. However, the spatiotemporal resolution of current experimental techniques limits our investigative capacity to understand this phenomenon. Here, we created a detailed three-dimensional (3D) cardiomyocyte model to study the subcellular regulatory mechanisms of myocardial energetics. The 3D cardiomyocyte model was based on the finite-element method, with detailed subcellular structures reproduced, and it included all elementary processes involved in cardiomyocyte electrophysiology, contraction, and ATP metabolism localized to specific loci. The simulation results were found to be reproducible and consistent with experimental data regarding the spatiotemporal pattern of cytosolic, intrasarcoplasmic-reticulum, and mitochondrial changes in Ca(2+); as well as changes in metabolite levels. Detailed analysis suggested that although the observed large cytosolic Ca(2+) gradient facilitated uptake by the mitochondrial Ca(2+) uniporter to produce cyclic changes in mitochondrial Ca(2+) near the Z-line region, the average mitochondrial Ca(2+) changes slowly. We also confirmed the importance of the creatine phosphate shuttle in cardiac energy regulation. In summary, our 3D model provides a powerful tool for the study of cardiac function by overcoming some of the spatiotemporal limitations of current experimental approaches.

Transmural and apicobasal gradients in repolarization contribute to T-wave genesis in human surface ECG.

The cellular basis of the T-wave morphology of surface ECG remains controversial in clinical cardiology. We examined the effect of action potential duration (APD) distribution on T-wave morphology using a realistic model of the human ventricle and torso. We developed a finite-element model of the ventricle consisting of ∼26 million elements, including the conduction system, each implemented with the ion current model of cardiomyocytes. This model was embedded in a torso model with distinct organ structures to obtain the standard ECG leads. The APD distribution was changed in the transmural direction by locating the M cells in either the endocardial or epicardial region. We also introduced apicobasal gradients by modifying the ion channel parameters. Both the transmural gradient (with M cells on the endocardial side) and the apicobasal gradient produced positive T waves, although a very large gradient was required for the apicobasal gradient. By contrast, T waves obtained with the transmural gradient were highly symmetric and, therefore, did not represent the true physiological state. Only combination of the transmural and the moderate apicobasal gradients produced physiological T waves in surface ECG. Positive T waves in surface ECG mainly originated from the transmural distribution of APD with M cells on the endocardial side, although the apicobasal gradient was also required to attain the physiological waveform.

Semi-Implicit Time Integration with Hessian Eigenvalue Corrections for a Larger Time Step in Molecular Dynamics Simulations.

In molecular dynamics simulations, the limited time step size has been a barrier to simulating long-time behaviors. Implicit time integration methods allow markedly larger time steps than the standard explicit time method, although they have major drawbacks such as overheads solving linear systems and instability of Newton iterations. To overcome these issues, we propose a semi-implicit time integration scheme, the semi-implicit Hessian correction (SimHec) scheme, for overdamped Langevin dynamics. The method focuses on the Hessian matrices of bonded and nonbonded interactions, where components with large negative Hessian eigenvalues are cut off in the linear approximation of momentum equations to avoid instability. The narrow band Hessian matrix enables an efficient parallelized linear solution with an overlapping approximation. We tested SimHec for the interdomain fluctuations in adenylate kinase and the powerstroke transition of myosin II using a coarse-grained protein model. SimHec reproduced the same dynamics as the explicit method, although the transition dynamics tended to be accelerated and fluctuations in bonded potentials were slightly reduced. These deviations were corrected using a hybrid method, SimHec-H, which adds explicit time steps after the semi-implicit time step. The proposed scheme allowed us to use time steps 50-200 times larger than those in explicit time integration, which resulted in a speedup factor of 7-30 taking the overhead into account.

Chloroquine and hydroxychloroquine provoke arrhythmias at concentrations higher than those clinically used to treat COVID-19: A simulation study.

The risk of fatal arrhythmias is the major concern for using chloroquine (CQ) or hydroxychloroquine (HCQ) to treat coronavirus disease 2019 (COVID-19), but the reported number of life-threatening arrhythmic events or deaths is relatively small. The objective of this study was to assess the arrhythmogenic risk of these two drugs using a multiscale heart simulation, which allows testing even at high concentrations, including those that cause fatal arrhythmias. We measured the inhibitory action of CQ, HCQ, and HCQ with 30 µM azithromycin (AZ) on six ion currents (fast [INa] and late [INa,L] components of the sodium current, L-type calcium current [ICa,L], rapid [IKr/hERG], and slow [IKs] components of delayed rectifier potassium, and inward rectifier potassium [IK1]) over a wide range of concentrations using the automated patch-clamp system. Using the concentration-inhibition relationship that was thus obtained, we simulated the drug effects while increasing the concentration until the life-threatening arrhythmia, torsade de pointes (TdP), was observed. The obtained threshold concentrations for TdP were 12.5, 35, and 22.5 µM for CQ, HCQ, and HCQ with AZ, respectively. Adding therapeutic concentrations of mexiletine or verapamil successfully prevented the occurrence of TdP, and verapamil was more effective. CQ, HCQ, and HCQ with AZ thresholds for TdP were larger than both antiviral concentrations that were reported by in vitro experiments and free plasma concentrations that were attained by the clinically used dosage. The current simulation data provided a safety margin to the currently used clinical dose for CQ and HCQ/AZ. Study Highlights WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC? Despite the potent in vitro antiviral effect, clinical trials have failed to show the therapeutic effects of chloroquine (CQ) and hydroxychloroquine (HCQ)/azithromycin (AZ) to treat coronavirus disease 2019. Torsadogenic potentials may limit the dosage of these drugs, but the reported incidence of fatal arrhythmias is rare. WHAT QUESTION DID THIS STUDY ADDRESS? Our objective was to assess the arrhythmogenicity of CQ and HCQ/AZ over a wide range of drug concentrations using a multiscale heart simulation. WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE? Our study showed that CQ and HCQ/AZ do not induce fatal arrhythmias even at concentrations much higher than in vitro antiviral half-maximal effective concentration (EC50 ) values at which QT prolongation exceeds 150 ms. We also found that estimated free plasma concentrations of CQ and HCQ/AZ achieved by currently used dosing protocols are lower than the antiviral EC50 for these drugs. HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE? Our simulation data provided a safety margin to the currently used clinical dose for CQ and HCQ/AZ.

A Multiple Step Active Stiffness Integration Scheme to Couple a Stochastic Cross-Bridge Model and Continuum Mechanics for Uses in Both Basic Research and Clinical Applications of Heart Simulation.

In a multiscale simulation of a beating heart, the very large difference in the time scales between rapid stochastic conformational changes of contractile proteins and deterministic macroscopic outcomes, such as the ventricular pressure and volume, have hampered the implementation of an efficient coupling algorithm for the two scales. Furthermore, the consideration of dynamic changes of muscle stiffness caused by the cross-bridge activity of motor proteins have not been well established in continuum mechanics. To overcome these issues, we propose a multiple time step scheme called the multiple step active stiffness integration scheme (MusAsi) for the coupling of Monte Carlo (MC) multiple steps and an implicit finite element (FE) time integration step. The method focuses on the active tension stiffness matrix, where the active tension derivatives concerning the current displacements in the FE model are correctly integrated into the total stiffness matrix to avoid instability. A sensitivity analysis of the number of samples used in the MC model and the combination of time step sizes confirmed the accuracy and robustness of MusAsi, and we concluded that the combination of a 1.25 ms FE time step and 0.005 ms MC multiple steps using a few hundred motor proteins in each finite element was appropriate in the tradeoff between accuracy and computational time. Furthermore, for a biventricular FE model consisting of 45,000 tetrahedral elements, one heartbeat could be computed within 1.5 h using 320 cores of a conventional parallel computer system. These results support the practicality of MusAsi for uses in both the basic research of the relationship between molecular mechanisms and cardiac outputs, and clinical applications of perioperative prediction.

Using Systolic Local Mechanical Load to Predict Fiber Orientation in Ventricles.

A simple rule adopted for myofiber reorientation in the ventricles is pursued by taking the microscopic branching network of myocytes into account. The macroscopic active tension generated on the microscopic branching structure is modeled by a multidirectional active stress tensor, which is defined as a function of the strains in the branching directions. In our reorientation algorithm, the principal direction of the branching network is updated so that it turns in the direction of greater active tension in the isovolumetric systole. Updates are performed step-by-step after the mechanical equilibrium has been attained with the current fiber structure. Starting from a nearly flat distribution of the principal fiber orientation along the circumferential direction, the reoriented fiber helix angles range from 70 to 40° at epicardium and from 60 to 80° at endocardium, in agreement with experimental observations. The helical ventricular myocardial band of Torrent-Guasp's model and the apical spiral structure of Rushmer's model are also reconstructed by our algorithm. Applying our algorithm to the infarcted ventricle model, the fiber structure near the infarcted site is remodeled so that the helix angle becomes steeper with respect to the circumferential direction near the epicardial surface. Based on our numerical analysis, we draw the following conclusions. (i) The multidirectional active tension based on the microscopic branching network is potentially used to seek tighter connection with neighboring aggregates. (ii) The thickening and thinning transitions in response to active tension in each myocyte allow the macroscopic principal fiber orientation of the microscopic branching network to move toward the direction of greater active tension. (iii) The force-velocity relationship is the key factor in transferring the fiber shortening strain to the magnitude of active tensions used in the myofiber reorientation. (iv) The algorithm naturally leads to homogeneity in the macroscopic active tension and the fiber shortening strain, and results in near-optimal pumping performance. (v) However, the reorientation mechanism may degrade the pumping performance if there is severely inhomogeneous contractility resulting from infarction. Our goal is to provide a tool to predict the fiber architecture of various heart disease patients for numerical simulations of their treatment plans.

Personalized Perioperative Multi-scale, Multi-physics Heart Simulation of Double Outlet Right Ventricle.

For treatment of complex congenital heart disease, computer simulation using a three-dimensional heart model may help to improve outcomes by enabling detailed preoperative evaluations. However, no highly integrated model that accurately reproduces a patient's pathophysiology, which is required for this simulation has been reported. We modelled a case of complex congenital heart disease, double outlet right ventricle with ventricular septal defect and atrial septal defect. From preoperative computed tomography images, finite element meshes of the heart and torso were created, and cell model of cardiac electrophysiology and sarcomere dynamics was implemented. The parameter values of the heart model were adjusted to reproduce the patient's electrocardiogram and haemodynamics recorded preoperatively. Two options of in silico surgery were performed using this heart model, and the resulting changes in performance were examined. Preoperative and postoperative simulations showed good agreement with clinical records including haemodynamics and measured oxyhaemoglobin saturations. The use of a detailed sarcomere model also enabled comparison of energetic efficiency between the two surgical options. A novel in silico model of congenital heart disease that integrates molecular models of cardiac function successfully reproduces the observed pathophysiology. The simulation of postoperative state by in silico surgeries can help guide clinical decision-making.