Electrodynamic heart model construction and ECG simulation

Objectives: In this paper, we present a unified electrodynamic heart model that permits simulations of the body surface potentials generated by the heart in motion. The inclusion of motion in the heart model significantly improves the accuracy of the simulated body surface potentials and therefore also the 12-lead ECG. Methods: The key step is to construct an electromechanical heart model. The cardiac excitation propagation is simulated by an electrical heart model, and the resulting cardiac active forces are used to calculate the ventricular wall motion based on a mechanical model. The source-field point relative position changes during heart systole and diastole. These can be obtained, and then used to calculate body surface ECG based on the electrical heart-torso model. Results: An electromechanical biventricular heart model is constructed and a standard 12-lead ECG is simulated. Compared with a simulated ECG based on the static electrical heart model, the simulated ECG based on the dynamic heart model is more accordant with a clinically recorded ECG, especially for the ST segment and T wave of a V1-V6 lead ECG. For slight-degree myocardial ischemia ECG simulation, the ST segment and T wave changes can be observed from the simulated ECG based on a dynamic heart model, while the ST segment and T wave of simulated ECG based on a static heart model is almost unchanged when compared with a normal ECG. Conclusions: This study confirms the importance of the mechanical factor in the ECG simulation. The dynamic heart model could provide more accurate ECG simulation, especially for myocardial ischemia or infarction simulation, since the main ECG changes occur at the ST segment and T wave, which correspond with cardiac systole and diastole phases.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Motion analysis of right ventricular wall based on an electromechanical biventricular model

Based on our previously developed electrical heart model, an electromechanical biventricular model, which couples the electrical property and mechanical property of the heart, was constructed and the right ventricular wall motion and deformation was simulated using this model. The model was developed on the basis of composite material theory and finite element method. The excitation propagation was simulated by electrical heart model, and the resultant active forces were used to study the ventricular wall motion during systole. The simulation results show that: (1) The right ventricular free wall moves towards the septum, and at the same time, the base and middle of free wall move towards the apex, which reduce the volume of right ventricle; (2) The minimum principle strain (E3) is largest at the apex, then at the middle of free wall, and its direction is in the approximate direction of epicardial muscle fibers. These results are in good accordance with solutions obtained from MR tagging images. It suggests that such electromechanical biventricular model can be used to assess the mechanical function of two ventricles.

An electrical heart model incorporating real geometry and motion

Abstract—This paper describes an electrical model of the ventricles incorporating real geometry and motion. Cardiac geometry and motion is obtained from segmentations of multipleslice MRI time sequences. A static heart model developed previously is deformed to match the observed geometry using a novel shape registration algorithm. The resulting electrocardiograms and body surface potential maps are compared to a static simulation in the resting heart. These results demonstrate that introducing motion into the cardiac model modifies the ECG during the T wave at peak contraction of the ventricles.

Hydrodynamic evaluation of kangaroo aortic valve matrices for tissue valve engineering

We evaluated the hydrodynamic performance of kangaroo aortic valve matrices (KMs) (19, 21, and 23 mm), as potential scaffolds in tissue valve engineering using a pulsatile left heart model at low and high cardiac outputs (COs) and heart rates (HRs) of 60 and 90 beats/min. Data were measured in two samples of each type, pooled in two CO levels (2.1 +/- 0.7 and 4.2 +/- 0.6 L/min; mean +/- standard errors on the mean), and analyzed using analysis of variance with CO level, HR, and valve type as fixed factors and compared to similar porcine matrices (PMs). Transvalvular pressure gradient (Delta P) was a function of HR (P < 0.001) and CO (P < 0.001) but not of valve type (P = 0.39). Delta P was consistently lower in KMs but not significantly different from PMs. The effective orifice area and performance index of kangaroo matrices was statistically larger for all sizes at both COs and HRs.

A comparison of different choices for the regularization parameter in inverse electrocardiography models

Calculating the potentials on the heart’s epicardial surface from the body surface potentials constitutes one form of inverse problems in electrocardiography (ECG). Since these problems are ill-posed, one approach is to use zero-order Tikhonov regularization, where the squared norms of both the residual and the solution are minimized, with a relative weight determined by the regularization parameter. In this paper, we used three different methods to choose the regularization parameter in the inverse solutions of ECG. The three methods include the L-curve, the generalized cross validation (GCV) and the discrepancy principle (DP). Among them, the GCV method has received less attention in solutions to ECG inverse problems than the other methods. Since the DP approach needs knowledge of norm of noises, we used a model function to estimate the noise. The performance of various methods was compared using a concentric sphere model and a real geometry heart-torso model with a distribution of current dipoles placed inside the heart model as the source. Gaussian measurement noises were added to the body surface potentials. The results show that the three methods all produce good inverse solutions with little noise; but, as the noise increases, the DP approach produces better results than the L-curve and GCV methods, particularly in the real geometry model. Both the GCV and L-curve methods perform well in low to medium noise situations.

Estimating trends and seasonality in coronary heart disease

We present two methods of estimating the trend, seasonality and noise in time series of coronary heart disease events. In contrast to previous work we use a non-linear trend, allow multiple seasonal components, and carefully examine the residuals from the fitted model. We show the importance of estimating these three aspects of the observed data to aid insight of the underlying process, although our major focus is on the seasonal components. For one method we allow the seasonal effects to vary over time and show how this helps the understanding of the association between coronary heart disease and varying temperature patterns. Copyright (C) 2004 John Wiley Sons, Ltd.

An integrated cognitive model of panic disorder: The role of positive and negative cognitions

One reason for the neglect of the role of positive factors in cognitive-behavioural therapy (CBT) may relate to a failure to develop cognitive models that integrate positive and negative cognitions. Bandura [Psychol. Rev. 84 (1977) 191; Anxiety Res. 1 (1988) 77] proposed that self-efficacy beliefs mediate a range of emotional and behavioural outcomes. However, in panic disorder, cognitively based research to date has largely focused on catastrophic misinterpretation of bodily sensations. Although a number of studies support each of the predictions associated with the account of panic disorder that is based on the role of negative cognitions, a review of the literature indicated that a cognitively based explanation of the disorder may be considerably strengthened by inclusion of positive cognitions that emphasize control or coping. Evidence to support an Integrated Cognitive Model (ICM) of panic disorder was examined and the theoretical implications of this model were discussed in terms of both schema change and compensatory skills accounts of change processes in CBT. (C) 2004 Elsevier Ltd. All rights reserved.

Determinants of functional capacity in patients with chronic heart failure: Role of filling pressure and systolic and diastolic function

Background Previous work suggesting a better correlation of diastolic than systolic function with exercise capacity in heart failure may reflect the -relative insensitivity and load-dependence of ejection fraction (EF). We sought the correlation of new and more sensitive methods of quantifying systolic and diastolic function and filling pressure with functional capacity. Methods We studied 155 consecutive exercise tests on 95 patients with congestive heart failure (81 male, aged 62 +/- 10 years), who underwent resting 2-climensional echocardiography and tissue Doppler imaging before and after measurement of maximum oxygen uptake (peak VO2)Results The resting EF was 3 1 % 10% and a peak VO(2)was 13 +/- 5 mL/kg/min; the majority of these patients (80%) had an ischemic cardiornyopathy. Resting EF (r 0.14, P =.09) correlated poorly with peak VO2 and mean systolic (r = 0.23, P =.004) and diastolic tissue velocities (r 0.18, P =.02). Peak EF was weakly correlated with the mean systolic (r = 0.18, P =.02) and diastolic velocities (r = 0.16, P <.04). The mean sum of systolic and diastolic velocities in both annuli (r = 0.30, P <.001) and E/Ea ratio (r 0.31, P <.001) were better correlated with peak VO2 Prediction of peak VO2 was similar with models based on models of filling pressure (R = 0.61), systolic factors (R = 0.63), and diastolic factors (R 0.59), although a composite model of filling pressure, systolic and diastolic function was a superior predictor of peak VO2 (R 0.69; all P<.001). Conclusions The reported association of diastolic rather than systolic function with functional capacity may have reflected the limitations of EF. Functional capacity appears related not only to diastolic function, but also to systolic function and filling pressure, and is most closely associated with a combination of these factors.

Physiological and haematological findings and clinical observations in a model of acute systemic anaphylaxis in Dirofilaria immitis-sensitised cats

Objective This study aims to understand the pathophysiology of anaphylaxis in Dirofflaria immitis-sensitised cats by analysing objective physiological and haematological measurements after challenge. Design Nineteen healthy D immitis-naive cats were sensitised using weekly injections of aluminium hydroxide-adjuvanted D immitis antigen, administered subcutaneously over 6 weeks. After sensitisation, cats (n = 16) were anaesthetised and challenged with intravenous D immitis antigen. A control group (n = 3) was sham-challenged using intravenous sterile 0.9% saline. Systolic blood pressure (measured non-invasively/indirectly), respiratory rate, degree of dyspnoea, blood 0, saturation, expired CO2, and heart rate and were measured immediately before and at 10 to 15 min intervals after challenge until terminal apnoea occurred or euthanasia at 140 mins post-challenge. Blood was collected for complete blood count immediately before and at 10, 20 and 35 mins after challenge. Clinical observations were recorded as they occurred. Results Antigen-challenged cats were divided into two groups: acute (apnoea occurred within 25 mins of challenge) and subacute (breathing at 25 mins after challenge). In both groups, the degree of dyspnoea increased and blood O-2 saturation and blood pressure decreased. Respiratory rate increased in the subacute group. Expired CO2 decreased in both Ag-challenged and control groups. Haematocrit increased in the subacute group. Neutrophil count decreased in the acute group and platelet count decreased in the subacute group. Eosinophil count decreased in the subacute and control groups. Sustained dyspnoea and gastrointestinal signs were the most common clinical manifestations of anaphylaxis in the antigen-challenged cats. Conclusions Intravenous challenge with D immitis antigen in sensitised cats results in dyspnoea, hypoxaemia and systemic hypotension accompanied by haemoconcentration.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Structural ambiguity of the Chinese version of the Hospital Anxiety and Depression Scale in patients with coronary heart disease

Background The Hospital Anxiety and Depression Scale (HADS) is a widely used screening tool designed as a case detector for clinically relevant anxiety and depression. Recent studies of the HADS in coronary heart disease (CHD) patients in European countries suggest it comprises three, rather than two, underlying sub-scale dimensions. The factor structure of the Chinese version of the HADS was evaluated in patients with CHD in mainland China. Methods Confirmatory factor analysis (CFA) was conducted on self-report HADS forms from 154 Chinese CHD patients. Results Little difference was observed in model fit between best performing three-factor and two-factor models. Conclusion The current observations are inconsistent with recent studies highlighting a dominant underlying tri-dimensional structure to the HADS in CHD patients. The Chinese version of the HADS may perform differently to European language versions of the instrument in patients with CHD.