RESUMO
In patients with dyssynchronous heart failure (DHF), cardiac conduction abnormalities cause the regional distribution of myocardial work to be non-homogeneous. Cardiac resynchronization therapy (CRT) using an implantable, programmed biventricular pacemaker/defibrillator, can improve the synchrony of contraction between the right and left ventricles in DHF, resulting in reduced morbidity and mortality and increased quality of life. Since regional work depends on wall stress, which cannot be measured in patients, we used computational methods to investigate regional work distributions and their changes after CRT. We used three-dimensional multi-scale patient-specific computational models parameterized by anatomic, functional, hemodynamic, and electrophysiological measurements in eight patients with heart failure and left bundle branch block (LBBB) who received CRT. To increase clinical translatability, we also explored whether streamlined computational methods provide accurate estimates of regional myocardial work. We found that CRT increased global myocardial work efficiency with significant improvements in non-responders. Reverse ventricular remodeling after CRT was greatest in patients with the highest heterogeneity of regional work at baseline, however the efficacy of CRT was not related to the decrease in overall work heterogeneity or to the reduction in late-activated regions of high myocardial work. Rather, decreases in early-activated regions of myocardium performing negative myocardial work following CRT best explained patient variations in reverse remodeling. These findings were also observed when regional myocardial work was estimated using ventricular pressure as a surrogate for myocardial stress and changes in endocardial surface area as a surrogate for strain. These new findings suggest that CRT promotes reverse ventricular remodeling in human dyssynchronous heart failure by increasing regional myocardial work in early-activated regions of the ventricles, where dyssynchrony is specifically associated with hypoperfusion, late systolic stretch, and altered metabolic activity and that measurement of these changes can be performed using streamlined approaches.
RESUMO
Electrical dyssynchrony leads to prestretch in late-activated regions and alters the sequence of mechanical contraction, although prestretch and its mechanisms are not well defined in the failing heart. We hypothesized that in heart failure, fiber prestretch magnitude increases with the amount of early-activated tissue and results in increased end-systolic strains, possibly due to length-dependent muscle properties. In five failing dog hearts with scars, three-dimensional strains were measured at the anterolateral left ventricle (LV). Prestretch magnitude was varied via ventricular pacing at increasing distances from the measurement site and was found to increase with activation time at various wall depths. At the subepicardium, prestretch magnitude positively correlated with the amount of early-activated tissue. At the subendocardium, local end-systolic strains (fiber shortening, radial wall thickening) increased proportionally to prestretch magnitude, resulting in greater mean strain values in late-activated compared with early-activated tissue. Increased fiber strains at end systole were accompanied by increases in preejection fiber strain, shortening duration, and the onset of fiber relengthening, which were all positively correlated with local activation time. In a dog-specific computational failing heart model, removal of length and velocity dependence on active fiber stress generation, both separately and together, alter the correlations between local electrical activation time and timing of fiber strains but do not primarily account for these relationships.
Assuntos
Insuficiência Cardíaca/fisiopatologia , Contração Miocárdica , Miocárdio/patologia , Taquicardia Ventricular/fisiopatologia , Função Ventricular Esquerda , Animais , Fenômenos Biomecânicos , Estimulação Cardíaca Artificial , Modelos Animais de Doenças , Cães , Eletrocardiografia , Técnicas Eletrofisiológicas Cardíacas , Análise de Elementos Finitos , Insuficiência Cardíaca/complicações , Insuficiência Cardíaca/patologia , Hemodinâmica , Imageamento por Ressonância Magnética , Modelos Cardiovasculares , Volume Sistólico , Sístole , Taquicardia Ventricular/complicações , Taquicardia Ventricular/patologia , Fatores de Tempo , Pressão VentricularRESUMO
AIMS: To test whether a functional growth law leads to asymmetric hypertrophy and associated changes in global and regional cardiac function when integrated with a computational model of left bundle branch block (LBBB). METHODS AND RESULTS: In recent studies, we proposed that cardiac myocytes grow longer when a threshold of maximum fibre strain is exceeded and grow thicker when the smallest maximum principal strain in the cellular cross-sectional plane exceeds a threshold. A non-linear cardiovascular model of the beating canine ventricles was combined with the cellular growth law. After inducing LBBB, the ventricles were allowed to adapt in shape over time in response to mechanical stimuli. When subjected to electrical dyssynchrony, the combined model of ventricular electromechanics, haemodynamics, and growth led to asymmetric hypertrophy with a faster increase of wall mass in the left ventricular (LV) free wall (FW) than the septum, increased LV end-diastolic and end-systolic volumes, and decreased LV ejection fraction. Systolic LV pressure decreased during the acute phase of LBBB and increased at later stages. The relative changes of these parameters were similar to those obtained experimentally. Most of the dilation was due to radial and axial fibre growth, and hence altered shape of the LVFW. CONCLUSION: Our previously proposed growth law reproduced measured dyssynchronously induced asymmetric hypertrophy and the associated functional changes, when combined with a computational model of the LBBB heart. The onset of LBBB leads to a step increase in LV mechanical discoordination that continues to increase as the heart remodels despite the constant electrical dyssynchrony.
Assuntos
Fibrilação Atrial/fisiopatologia , Acoplamento Excitação-Contração , Sistema de Condução Cardíaco/fisiopatologia , Ventrículos do Coração/fisiopatologia , Modelos Cardiovasculares , Contração Miocárdica , Disfunção Ventricular Esquerda/fisiopatologia , Animais , Fibrilação Atrial/complicações , Simulação por Computador , Cães , Humanos , Disfunção Ventricular Esquerda/etiologiaRESUMO
Adult cardiac muscle adapts to mechanical changes in the environment by growth and remodeling (G&R) via a variety of mechanisms. Hypertrophy develops when the heart is subjected to chronic mechanical overload. In ventricular pressure overload (e.g. due to aortic stenosis) the heart typically reacts by concentric hypertrophic growth, characterized by wall thickening due to myocyte radial growth when sarcomeres are added in parallel. In ventricular volume overload, an increase in filling pressure (e.g. due to mitral regurgitation) leads to eccentric hypertrophy as myocytes grow axially by adding sarcomeres in series leading to ventricular cavity enlargement that is typically accompanied by some wall thickening. The specific biomechanical stimuli that stimulate different modes of ventricular hypertrophy are still poorly understood. In a recent study, based on in-vitro studies in micropatterned myocyte cell cultures subjected to stretch, we proposed that cardiac myocytes grow longer to maintain a preferred sarcomere length in response to increased fiber strain and grow thicker to maintain interfilament lattice spacing in response to increased cross-fiber strain. Here, we test whether this growth law is able to predict concentric and eccentric hypertrophy in response to aortic stenosis and mitral valve regurgitation, respectively, in a computational model of the adult canine heart coupled to a closed loop model of circulatory hemodynamics. A non-linear finite element model of the beating canine ventricles coupled to the circulation was used. After inducing valve alterations, the ventricles were allowed to adapt in shape in response to mechanical stimuli over time. The proposed growth law was able to reproduce major acute and chronic physiological responses (structural and functional) when integrated with comprehensive models of the pressure-overloaded and volume-overloaded canine heart, coupled to a closed-loop circulation. We conclude that strain-based biomechanical stimuli can drive cardiac growth, including wall thickening during pressure overload.
RESUMO
Recently, attention has been focused on comparing left ventricular (LV) endocardial (ENDO) with epicardial (EPI) pacing for cardiac resynchronization therapy. However, the effects of ENDO and EPI lead placement at multiple sites have not been studied in failing hearts. We hypothesized that differences in the improvement of ventricular function due to ENDO vs. EPI pacing in dyssynchronous (DYSS) heart failure may depend on the position of the LV lead in relation to the original activation pattern. In six nonfailing and six failing dogs, electrical DYSS was created by atrioventricular sequential pacing of the right ventricular apex. ENDO was compared with EPI biventricular pacing at five LV sites. In failing hearts, increases in the maximum rate of LV pressure change (dP/dt; r = 0.64), ejection fraction (r = 0.49), and minimum dP/dt (r = 0.51), relative to DYSS, were positively correlated (P < 0.01) with activation time at the LV pacing site during ENDO but not EPI pacing. ENDO pacing at sites with longer activation delays led to greater improvements in hemodynamic parameters and was associated with an overall reduction in electrical DYSS compared with EPI pacing (P < 0.05). These findings were qualitatively similar for nonfailing hearts. Improvement in hemodynamic function increased with activation time at the LV pacing site during ENDO but not EPI pacing. At the anterolateral wall, end-systolic transmural function was greater with local ENDO compared with EPI pacing. ENDO pacing and intrinsic activation delay may have important implications for management of DYSS heart failure.
Assuntos
Estimulação Cardíaca Artificial/métodos , Endocárdio/fisiologia , Insuficiência Cardíaca/fisiopatologia , Insuficiência Cardíaca/terapia , Coração/fisiologia , Animais , Valva Aórtica/fisiologia , Fenômenos Biomecânicos , Cinerradiografia , Cães , Eletrocardiografia , Eletrodos Implantados , Ventrículos do Coração , Hemodinâmica/fisiologia , Técnicas In Vitro , Valva Mitral/fisiologia , Inclusão em Parafina , Taquicardia/fisiopatologia , Função Ventricular Esquerda/fisiologia , Função Ventricular Direita/fisiologiaRESUMO
Changes in muscle fiber orientation across the wall of the left ventricle (LV) cause the apex of the heart to turn 10-15 deg in opposition to its base during systole and are believed to increase stroke volume and lower wall stress in healthy hearts. Studies show that cardiac torsion is sensitive to various disease states, which suggests that it may be an important aspect of cardiac function. Modern imaging techniques have sparked renewed interest in cardiac torsion dynamics, but no work has been done to determine whether mechanically augmented apical torsion can be used to restore function to failing hearts. In this report, we discuss the potential advantages of this approach and present evidence that turning the cardiac apex by mechanical means can displace a clinically significant volume of blood from failing hearts. Computational models of normal and reduced-function LVs were created to predict the effects of applied apical torsion on ventricular stroke work and wall stress. These same conditions were reproduced in anesthetized pigs with drug-induced heart failure using a custom apical torsion device programmed to rotate over various angles during cardiac systole. Simulations of applied 90 deg torsion in a prolate spheroidal computational model of a reduced-function pig heart produced significant increases in stroke work (25%) and stroke volume with reduced fiber stress in the epicardial region. These calculations were in substantial agreement with corresponding in vivo measurements. Specifically, the computer model predicted torsion-induced stroke volume increases from 13.1 to 14.4 mL (9.9%) while actual stroke volume in a pig heart of similar size and degree of dysfunction increased from 11.1 to 13.0 mL (17.1%). Likewise, peak LV pressures in the computer model rose from 85 to 95 mm Hg (11.7%) with torsion while maximum ventricular pressures in vivo increased in similar proportion, from 55 to 61 mm Hg (10.9%). These data suggest that: (a) the computer model of apical torsion developed for this work is a fair and accurate predictor of experimental outcomes, and (b) supra-physiologic apical torsion may be a viable means to boost cardiac output while avoiding blood contact that occurs with other assist methods.
Assuntos
Débito Cardíaco/fisiologia , Modelos Cardiovasculares , Contração Miocárdica/fisiologia , Volume Sistólico/fisiologia , Torção Mecânica , Animais , Simulação por Computador , Coração/fisiopatologia , Insuficiência Cardíaca/terapia , Ventrículos do Coração/patologia , Rotação , Suínos , Sístole/fisiologia , Disfunção Ventricular Esquerda/terapiaRESUMO
The excitation-contraction coupling properties of cardiac myocytes isolated from different regions of the mammalian left ventricular wall have been shown to vary considerably, with uncertain effects on ventricular function. We embedded a cell-level excitation-contraction coupling model with region-dependent parameters within a simple finite element model of left ventricular geometry to study effects of electromechanical heterogeneity on local myocardial mechanics and global haemodynamics. This model was compared with one in which heterogeneous myocyte parameters were assigned randomly throughout the mesh while preserving the total amount of each cell subtype. The two models displayed nearly identical transmural patterns of fibre and cross-fibre strains at end-systole, but showed clear differences in fibre strains at earlier points during systole. Haemodynamic function, including peak left ventricular pressure, maximal rate of left ventricular pressure development and stroke volume, were essentially identical in the two models. These results suggest that in the intact ventricle heterogeneously distributed myocyte subtypes primarily impact local deformation of the myocardium, and that these effects are greatest during early systole.
Assuntos
Modelos Cardiovasculares , Miócitos Cardíacos/fisiologia , Função Ventricular Esquerda/fisiologia , Potenciais de Ação , Animais , Fenômenos Biomecânicos , Cães , Fenômenos Eletrofisiológicos , Análise de Elementos Finitos , Contração Miocárdica/fisiologiaRESUMO
Advances in computer power, novel diagnostic and therapeutic medical technologies, and an increasing knowledge of pathophysiology from gene to organ systems make it increasingly feasible to apply multiscale patient-specific modeling based on proven disease mechanisms. Such models may guide and predict the response to therapy in many areas of medicine. This is an exciting and relatively new approach, for which efficient methods and computational tools are of the utmost importance. Investigators have designed patient-specific models in almost all areas of human physiology. Not only will these models be useful in clinical settings to predict and optimize the outcome from surgery and non-interventional therapy, but they will also provide pathophysiologic insights from the cellular level to the organ system level. Models, therefore, will provide insight as to why specific interventions succeed or fail.
Assuntos
Simulação por Computador , Cardiopatias , Modelos Teóricos , Inteligência Artificial , Protocolos Clínicos , Cardiopatias/diagnóstico , Cardiopatias/terapia , Humanos , SoftwareRESUMO
INTRODUCTION: The postnatal heart grows mostly in response to increased hemodynamic load. However, the specific biomechanical stimuli that stimulate cardiac growth as a reaction to increased hemodynamic load are still poorly understood. It has been shown that isolated neonatal rat cardiac myocytes normalize resting sarcomere length by adding sarcomeres in series when subjected to uniaxial static strain. Because there is experimental evidence that myocytes can distinguish the direction of stretch, it was postulated that myocytes also may normalize interfilament lattice spacing as a response to cross-fiber stretch. METHODS: A growth law was proposed in which fiber axial growth was stimulated by fiber strain deviating from zero and fiber radial growth by cross-fiber strain (parallel to the wall surface) deviating from zero. Fiber radial growth rate constant was 1/3 of the fiber axial growth rate constant. The growth law was implemented in a finite element model of the newborn Sprague-Dawley rat residually stressed left ventricle (LV). The LV was subjected to an end-diastolic pressure of 1 kPa and about 25 weeks of normal growth was simulated. RESULTS: Most cellular and chamber dimension changes in the model matched experimentally measured ones: LV cavity and wall volume increased from 2.3 and 54 µl, respectively, in the newborn to 276 µl and 1.1 ml, respectively, in the adult rat; LV shape became more spherical; internal LV radius increased faster than wall thickness; and unloaded sarcomere lengths exhibited a transmural gradient. The major discrepancy with experiments included a reversed transmural gradient of cell length in the older rat. CONCLUSION: A novel strain-based growth law has been presented that reproduced physiological postnatal growth in the rat LV.
Assuntos
Coração/crescimento & desenvolvimento , Animais , Fenômenos Biomecânicos/fisiologia , Pressão Sanguínea/fisiologia , Simulação por Computador , Hemodinâmica/fisiologia , Ratos , Ratos Sprague-Dawley , Sarcômeros/fisiologia , Estresse Mecânico , Função Ventricular Esquerda/fisiologiaRESUMO
The development and clinical use of patient-specific models of the heart is now a feasible goal. Models have the potential to aid in diagnosis and support decision-making in clinical cardiology. Several groups are now working on developing multi-scale models of the heart for understanding therapeutic mechanisms and better predicting clinical outcomes of interventions such as cardiac resynchronization therapy. Here we describe the methodology for generating a patient-specific model of the failing heart with a myocardial infarct and left ventricular bundle branch block. We discuss some of the remaining challenges in developing reliable patient-specific models of cardiac electromechanical activity, and identify some of the main areas for focusing future research efforts. Key challenges include: efficiently generating accurate patient-specific geometric meshes and mapping regional myofiber architecture to them; modeling electrical activation patterns based on cellular alterations in human heart failure, and estimating regional tissue conductivities based on clinically available electrocardiographic recordings; estimating unloaded ventricular reference geometry and material properties for biomechanical simulations; and parameterizing systemic models of circulatory dynamics from available hemodynamic measurements.
Assuntos
Insuficiência Cardíaca/patologia , Insuficiência Cardíaca/fisiopatologia , Modelos Biológicos , Idoso , Fenômenos Biomecânicos , Fenômenos Eletrofisiológicos , Insuficiência Cardíaca/complicações , Ventrículos do Coração/patologia , Ventrículos do Coração/fisiopatologia , Hemodinâmica , Humanos , Masculino , Modelos Anatômicos , Contração Muscular , Infarto do Miocárdio/complicações , Medicina de PrecisãoRESUMO
BACKGROUND: Heart failure (HF) in combination with mechanical dyssynchrony is associated with a high mortality rate. To quantify contractile dysfunction in patients with HF, investigators have proposed several indices of mechanical dyssynchrony, including percentile range of time to peak shortening (WTpeak), circumferential uniformity ratio estimate (CURE), and internal stretch fraction (ISF). The goal of this study was to compare the sensitivity of these indices to 4 major abnormalities responsible for cardiac dysfunction in dyssynchronous HF: dilation, negative inotropy, negative lusitropy, and dyssynchronous activation. METHODS AND RESULTS: All combinations of these 4 major abnormalities were included in 3D computational models of ventricular electromechanics. Compared with a nonfailing heart model, ventricles were dilated, inotropy was reduced, twitch duration was prolonged, and activation sequence was changed from normal to left bundle branch block. In the nonfailing heart, CURE, ISF, and WTpeak were 0.97+/-0.004, 0.010+/-0.002, and 78+/-1 milliseconds, respectively. With dilation alone, CURE decreased 2.0+/-0.07%, ISF increased 58+/-47%, and WTpeak increased 31+/-3%. With dyssynchronous activation alone, CURE decreased 15+/-0.6%, ISF increased 14-fold (+/-3), and WTpeak increased 121+/-4%. With the combination of dilation and dyssynchronous activation, CURE decreased 23+/-0.8%, ISF increased 20-fold (+/-5), and WTpeak increased 147+/-5%. CONCLUSIONS: Dilation and left bundle branch block combined synergistically decreased regional cardiac function. CURE and ISF were sensitive to this combination, but WTpeak was not. CURE and ISF also reflected the relative nonuniform distribution of regional work better than WTpeak. These findings might explain why CURE and ISF are better predictors of reverse remodeling in cardiac resynchronization therapy.
Assuntos
Simulação por Computador , Insuficiência Cardíaca Diastólica/etiologia , Insuficiência Cardíaca Sistólica/etiologia , Imageamento Tridimensional , Modelos Cardiovasculares , Disfunção Ventricular Esquerda/complicações , Disfunção Ventricular Esquerda/fisiopatologia , Bloqueio de Ramo/complicações , Bloqueio de Ramo/diagnóstico , Dilatação Patológica/complicações , Dilatação Patológica/diagnóstico , Eletrofisiologia , Insuficiência Cardíaca Diastólica/diagnóstico , Insuficiência Cardíaca Sistólica/diagnóstico , Testes de Função Cardíaca , Humanos , Interpretação de Imagem Assistida por Computador , Contração Miocárdica , Valores de Referência , Sensibilidade e Especificidade , Processamento de Sinais Assistido por Computador , Volume Sistólico/fisiologia , Disfunção Ventricular Esquerda/diagnósticoRESUMO
OBJECTIVES: To study the impact of biventricular pacing (BiV) and scar size on left ventricular (LV) regional and global function using a detailed finite element model of three-dimensional electromechanics in the failing canine heart. BACKGROUND: Cardiac resynchronization therapy (CRT) clinical trials have demonstrated that up to 30% of patients may be classified as non-responders. The presence of a scar appears to contribute to those that do not respond to CRT. A recent study in patients with myocardial scar showed that LV dyssynchrony was the sole independent predictor of reverse remodeling, and not scar location or size. METHODS: Two activation sequences were simulated: left bundle branch block (LBBB) and acute simultaneous BiV (with leads in the left and right ventricle) in hearts with chronic scars of various sizes. The dependence of regional function (mean fiber ejection strain, variance of fiber isovolumic strain and fraction of tissue stretched during ejection) and global function (left ventricular dP/dt(max), ejection fraction, stroke work) on scar size and pacing protocol was tested. RESULTS: Global function and regional function averaged over the whole LV during LBBB and BiV decreased as a function of scar size. In the non-scarred regions, however, regional function was largely independent of scar size for a fixed pacing site. CONCLUSIONS: The model results suggest that uniformity of mechanical contraction in non-scarred regions in the failing heart during biventricular pacing is independent of scar size for a fixed pacing site.
Assuntos
Bloqueio de Ramo/prevenção & controle , Bloqueio de Ramo/fisiopatologia , Estimulação Cardíaca Artificial/métodos , Cicatriz/fisiopatologia , Insuficiência Cardíaca/prevenção & controle , Insuficiência Cardíaca/fisiopatologia , Modelos Cardiovasculares , Animais , Bloqueio de Ramo/complicações , Simulação por Computador , Cães , Módulo de Elasticidade , Técnicas de Imagem por Elasticidade/métodos , Insuficiência Cardíaca/etiologia , Humanos , Resultado do Tratamento , Disfunção Ventricular Esquerda/etiologia , Disfunção Ventricular Esquerda/fisiopatologia , Disfunção Ventricular Esquerda/prevenção & controleRESUMO
Here we describe new computational models of cardiac electromechanics starting from the cellular scale and building to the tissue, organ and system scales. We summarize application to human genetic diseases (LQT1 and LQT3) and to modeling of congestive heart failure.
Assuntos
Insuficiência Cardíaca/fisiopatologia , Coração/fisiopatologia , Mutação , Algoritmos , Simulação por Computador , Eletrofisiologia/métodos , Desenho de Equipamento , Sistema de Condução Cardíaco , Humanos , Síndrome do QT Longo , Modelos Cardiovasculares , Modelos Genéticos , Células Musculares/patologia , Contração Miocárdica , SoftwareRESUMO
In this study we present a novel, robust method to couple finite element (FE) models of cardiac mechanics to systems models of the circulation (CIRC), independent of cardiac phase. For each time step through a cardiac cycle, left and right ventricular pressures were calculated using ventricular compliances from the FE and CIRC models. These pressures served as boundary conditions in the FE and CIRC models. In succeeding steps, pressures were updated to minimize cavity volume error (FE minus CIRC volume) using Newton iterations. Coupling was achieved when a predefined criterion for the volume error was satisfied. Initial conditions for the multi-scale model were obtained by replacing the FE model with a varying elastance model, which takes into account direct ventricular interactions. Applying the coupling, a novel multi-scale model of the canine cardiovascular system was developed. Global hemodynamics and regional mechanics were calculated for multiple beats in two separate simulations with a left ventricular ischemic region and pulmonary artery constriction, respectively. After the interventions, global hemodynamics changed due to direct and indirect ventricular interactions, in agreement with previously published experimental results. The coupling method allows for simulations of multiple cardiac cycles for normal and pathophysiology, encompassing levels from cell to system.
Assuntos
Velocidade do Fluxo Sanguíneo/fisiologia , Pressão Sanguínea/fisiologia , Circulação Colateral/fisiologia , Modelos Cardiovasculares , Circulação Pulmonar/fisiologia , Fluxo Pulsátil/fisiologia , Função Ventricular , Animais , Circulação Sanguínea , Simulação por Computador , CãesRESUMO
INTRODUCTION: For efficient pump function, contraction of the heart should be as synchronous as possible. Ventricular pacing induces asynchrony of depolarization and contraction. The degree of asynchrony depends on the position of the pacing electrode. The aim of this study was to extend an existing numerical model of electromechanics in the left ventricle (LV) to the application of ventricular pacing. With the model, the relation between pacing site and patterns of depolarization and contraction was investigated. METHODS AND RESULTS: The LV was approximated by a thick-walled ellipsoid with a realistic myofiber orientation. Propagation of the depolarization wave was described by the eikonal-diffusion equation, in which five parameters play a role: myocardial and subendocardial velocity of wave propagation along the myofiber cm and ce; myocardial and subendocardial anisotropy am and ae; and parameter k, describing the influence of wave curvature on wave velocity. Parameters cm, ae, and k were taken from literature. Parameters am and ce were estimated by fitting the model to experimental data, obtained by pacing the canine left ventricular free wall (LVFW). The best fit was found with cm = 0.75 m/s, ce = 1.3 m/s, am = 2.5, ae = 1.5, and k = 2.1 x 10(-4) m2/s. With these parameter settings, for right ventricular apex (RVA) pacing, the depolarization times were realistically simulated as also shown by the wavefronts and the time needed to activate the LVFW. The moment of depolarization was used to initiate myofiber contraction in a model of LV mechanics. For both pacing situations, mid-wall circumferential strains and onset of myofiber shortening were obtained. CONCLUSION: With a relatively simple model setup, simulated depolarization timing patterns agreed with measurements for pacing at the LVFW and RVA in an LV. Myocardial cross-fiber wave velocity is estimated to be 0.40 times the velocity along the myofiber direction (0.75 m/s). Subendocardial wave velocity is about 1.7 times faster than in the rest of the myocardium, but about 3 times slower than as found in Purkinje fibers. Furthermore, model and experiment agreed in the following respects. (1) Ventricular pacing decreased both systolic pressure and ejection fraction relative to natural sinus rhythm. (2) In early depolarized regions, early shortening was observed in the isovolumic contraction phase; in late depolarized regions, myofibers were stretched in this phase. Maps showing timing of onset of shortening were similar to previously measured maps in which wave velocity of contraction appeared similar to that of depolarization.