Virtual clinical QT exposure-response studies - A translational computational approach.

BACKGROUND AND PURPOSE: A recent paradigm shift in proarrhythmic risk assessment suggests that the integration of clinical, non-clinical, and computational evidence can be used to reach a comprehensive understanding of the proarrhythmic potential of drug candidates. While current computational methodologies focus on predicting the incidence of proarrhythmic events after drug administration, the objective of this study is to predict concentration-response relationships of QTc as a clinical endpoint. EXPERIMENTAL APPROACH: Full heart computational models reproducing human cardiac populations were created to predict the concentration-response relationship of changes in the QT interval as recommended for clinical trials. The concentration-response relationship of the QT-interval prolongation obtained from the computational cardiac population was compared against the relationship from clinical trial data for a set of well-characterized compounds: moxifloxacin, dofetilide, verapamil, and ondansetron. KEY RESULTS: Computationally derived concentration-response relationships of QT interval changes for three of the four drugs had slopes within the confidence interval of clinical trials (dofetilide, moxifloxacin and verapamil) when compared to placebo-corrected concentration-ΔQT and concentration-ΔQT regressions. Moxifloxacin showed a higher intercept, outside the confidence interval of the clinical data, demonstrating that in this example, the standard linear regression does not appropriately capture the concentration-response results at very low concentrations. The concentrations corresponding to a mean QTc prolongation of 10 ms were consistently lower in the computational model than in clinical data. The critical concentration varied within an approximate ratio of 0.5 (moxifloxacin and ondansetron) and 1 times (dofetilide, verapamil) the critical concentration observed in human clinical trials. Notably, no other in silico methodology can approximate the human critical concentration values for a QT interval prolongation of 10 ms. CONCLUSION AND IMPLICATIONS: Computational concentration-response modelling of a virtual population of high-resolution, 3-dimensional cardiac models can provide comparable information to clinical data and could be used to complement pre-clinical and clinical safety packages. It provides access to an unlimited exposure range to support trial design and can improve the understanding of pre-clinical-clinical translation.

Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models.

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O'Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studies.

Study protocol: MyoFit46-the cardiac sub-study of the MRC National Survey of Health and Development.

BACKGROUND: The life course accumulation of overt and subclinical myocardial dysfunction contributes to older age mortality, frailty, disability and loss of independence. The Medical Research Council National Survey of Health and Development (NSHD) is the world's longest running continued surveillance birth cohort providing a unique opportunity to understand life course determinants of myocardial dysfunction as part of MyoFit46-the cardiac sub-study of the NSHD. METHODS: We aim to recruit 550 NSHD participants of approximately 75 years+ to undertake high-density surface electrocardiographic imaging (ECGI) and stress perfusion cardiovascular magnetic resonance (CMR). Through comprehensive myocardial tissue characterization and 4-dimensional flow we hope to better understand the burden of clinical and subclinical cardiovascular disease. Supercomputers will be used to combine the multi-scale ECGI and CMR datasets per participant. Rarely available, prospectively collected whole-of-life data on exposures, traditional risk factors and multimorbidity will be studied to identify risk trajectories, critical change periods, mediators and cumulative impacts on the myocardium. DISCUSSION: By combining well curated, prospectively acquired longitudinal data of the NSHD with novel CMR-ECGI data and sharing these results and associated pipelines with the CMR community, MyoFit46 seeks to transform our understanding of how early, mid and later-life risk factor trajectories interact to determine the state of cardiovascular health in older age. TRIAL REGISTRATION: Prospectively registered on ClinicalTrials.gov with trial ID: 19/LO/1774 Multimorbidity Life-Course Approach to Myocardial Health- A Cardiac Sub-Study of the MCRC National Survey of Health and Development (NSHD).

Concomitant Respiratory Failure Can Impair Myocardial Oxygenation in Patients with Acute Cardiogenic Shock Supported by VA-ECMO.

Venous-arterial extracorporeal membrane oxygenation (VA-ECMO) treatment for acute cardiogenic shock in patients who also have acute lung injury predisposes development of a serious complication called "north-south syndrome" (NSS) which causes cerebral hypoxia. NSS is poorly characterized and hemodynamic studies have focused on cerebral perfusion ignoring the heart. We hypothesized in NSS the heart would be more likely to receive hypoxemic blood than the brain due to the proximity of the coronary arteries to the aortic annulus. To test this, we conducted a computational fluid dynamics simulation of blood flow in a human supported by VA-ECMO. Simulations quantified the fraction of blood at each aortic branching vessel originating from residual native cardiac output versus VA-ECMO. As residual cardiac function was increased, simulations demonstrated myocardial hypoxia would develop prior to cerebral hypoxia. These results illustrate the conditions where NSS will develop and the relative cardiac function that will lead to organ-specific hypoxia. Illustration of the impact of north-south syndrome on organ-specific oxygen delivery. Patients on VA-ECMO have two sources of blood flow, one from the VA-ECMO circuit and one from the residual cardiac function. When there is no residual cardiac function, all organs are perfused with oxygenated blood. As myocardial recovery progresses, blood supply from the two sources will begin to mix resulting in non-homogeneous mixing and differential oxygenation based upon the anatomical site of branching vessels.

Time-efficient three-dimensional transmural scar assessment provides relevant substrate characterization for ventricular tachycardia features and long-term recurrences in ischemic cardiomyopathy.

Delayed gadolinium-enhanced cardiac magnetic resonance (LGE-CMR) imaging requires novel and time-efficient approaches to characterize the myocardial substrate associated with ventricular arrhythmia in patients with ischemic cardiomyopathy. Using a translational approach in pigs and patients with established myocardial infarction, we tested and validated a novel 3D methodology to assess ventricular scar using custom transmural criteria and a semiautomatic approach to obtain transmural scar maps in ventricular models reconstructed from both 3D-acquired and 3D-upsampled-2D-acquired LGE-CMR images. The results showed that 3D-upsampled models from 2D LGE-CMR images provided a time-efficient alternative to 3D-acquired sequences to assess the myocardial substrate associated with ischemic cardiomyopathy. Scar assessment from 2D-LGE-CMR sequences using 3D-upsampled models was superior to conventional 2D assessment to identify scar sizes associated with the cycle length of spontaneous ventricular tachycardia episodes and long-term ventricular tachycardia recurrences after catheter ablation. This novel methodology may represent an efficient approach in clinical practice after manual or automatic segmentation of myocardial borders in a small number of conventional 2D LGE-CMR slices and automatic scar detection.

Three-dimensional cardiac fibre disorganization as a novel parameter for ventricular arrhythmia stratification after myocardial infarction.

AIMS: Myocardial infarction (MI) alters cardiac fibre organization with unknown consequences on ventricular arrhythmia. We used diffusion tensor imaging (DTI) of three-dimensional (3D) cardiac fibres and scar reconstructions to identify the main parameters associated with ventricular arrhythmia inducibility and ventricular tachycardia (VT) features after MI. METHODS AND RESULTS: Twelve pigs with established MI and three controls underwent invasive electrophysiological characterization of ventricular arrhythmia inducibility and VT features. Animal-specific 3D scar and myocardial fibre distribution were obtained from ex vivo high-resolution contrast-enhanced T1 mapping and DTI sequences. Diffusion tensor imaging-derived parameters significantly different between healthy and scarring myocardium, scar volumes, and left ventricular ejection fraction (LVEF) were included for arrhythmia risk stratification and correlation analyses with VT features. Ventricular fibrillation (VF) was the only inducible arrhythmia in 4 out of 12 infarcted pigs and all controls. Ventricular tachycardia was also inducible in the remaining eight pigs during programmed ventricular stimulation. A DTI-based 3D fibre disorganization index (FDI) showed higher disorganization within dense scar regions of VF-only inducible pigs compared with VT inducible animals (FDI: 0.36; 0.36-0.37 vs. 0.32; 0.26-0.33, respectively, P = 0.0485). Ventricular fibrillation induction required lower programmed stimulation aggressiveness in VF-only inducible pigs than VT inducible and control animals. Neither LVEF nor scar volumes differentiated between VF and VT inducible animals. Re-entrant VT circuits were localized within areas of highly disorganized fibres. Moreover, the FDI within heterogeneous scar regions was associated with the median VT cycle length per animal (R2 = 0.5320). CONCLUSION: The amount of scar-related cardiac fibre disorganization in DTI sequences is a promising approach for ventricular arrhythmia stratification after MI.

Implications of bipolar voltage mapping and magnetic resonance imaging resolution in biventricular scar characterization after myocardial infarction.

AIMS: We aimed to study the differences in biventricular scar characterization using bipolar voltage mapping compared with state-of-the-art in vivo delayed gadolinium-enhanced cardiac magnetic resonance (LGE-CMR) imaging and ex vivo T1 mapping. METHODS AND RESULTS: Ten pigs with established myocardial infarction (MI) underwent in vivo scar characterization using LGE-CMR imaging and high-density voltage mapping of both ventricles using a 3.5-mm tip catheter. Ex vivo post-contrast T1 mapping provided a high-resolution reference. Voltage maps were registered onto the left and right ventricular (LV and RV) endocardium, and epicardium of CMR-based geometries to compare voltage-derived scars with surface-projected 3D scars. Voltage-derived scar tissue of the LV endocardium and the epicardium resembled surface projections of 3D in vivo and ex vivo CMR-derived scars using 1-mm of surface projection distance. The thinner wall of the RV was especially sensitive to lower resolution in vivo LGE-CMR images, in which differences between normalized low bipolar voltage areas and CMR-derived scar areas did not decrease below a median of 8.84% [interquartile range (IQR) (3.58, 12.70%)]. Overall, voltage-derived scars and surface scar projections from in vivo LGE-CMR sequences showed larger normalized scar areas than high-resolution ex vivo images [12.87% (4.59, 27.15%), 18.51% (11.25, 24.61%), and 9.30% (3.84, 19.59%), respectively], despite having used optimized surface projection distances. Importantly, 43.02% (36.54, 48.72%) of voltage-derived scar areas from the LV endocardium were classified as non-enhanced healthy myocardium using ex vivo CMR imaging. CONCLUSION: In vivo LGE-CMR sequences and high-density voltage mapping using a conventional linear catheter fail to provide accurate characterization of post-MI scar, limiting the specificity of voltage-based strategies and imaging-guided procedures.

Complex Congenital Heart Disease Associated With Disordered Myocardial Architecture in a Midtrimester Human Fetus.

BACKGROUND: In the era of increasingly successful corrective interventions in patients with congenital heart disease (CHD), global and regional myocardial remodeling are emerging as important sources of long-term morbidity/mortality. Changes in organization of the myocardium in CHD, and in its mechanical properties, conduction, and blood supply, result in altered myocardial function both before and after surgery. To gain a better understanding and develop appropriate and individualized treatment strategies, the microscopic organization of cardiomyocytes, and their integration at a macroscopic level, needs to be completely understood. The aim of this study is to describe, for the first time, in 3 dimensions and nondestructively the detailed remodeling of cardiac microstructure present in a human fetal heart with complex CHD. METHODS AND RESULTS: Synchrotron X-ray phase-contrast imaging was used to image an archival midgestation formalin-fixed fetal heart with right isomerism and complex CHD and compare with a control fetal heart. Analysis of myocyte aggregates, at detail not accessible with other techniques, was performed. Macroanatomic and conduction system changes specific to the disease were clearly observable, together with disordered myocyte organization in the morphologically right ventricle myocardium. Electrical activation simulations suggested altered synchronicity of the morphologically right ventricle. CONCLUSIONS: We have shown the potential of X-ray phase-contrast imaging for studying cardiac microstructure in the developing human fetal heart at high resolution providing novel insight while preserving valuable archival material for future study. This is the first study to show myocardial alterations occur in complex CHD as early as midgestation.

Fully coupled fluid-electro-mechanical model of the human heart for supercomputers.

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

Evaluating the roles of detailed endocardial structures on right ventricular haemodynamics by means of CFD simulations.

Computational modelling plays an important role in right ventricular (RV) haemodynamic analysis. However, current approaches use smoothed ventricular anatomies. The aim of this study is to characterise RV haemodynamics including detailed endocardial structures like trabeculae, moderator band, and papillary muscles. Four paired detailed and smoothed RV endocardium models (2 male and 2 female) were reconstructed from ex vivo human hearts high-resolution magnetic resonance images. Detailed models include structures with ≥1 mm2 cross-sectional area. Haemodynamic characterisation was done by computational fluid dynamics simulations with steady and transient inflows, using high-performance computing. The differences between the flows in smoothed and detailed models were assessed using Q-criterion for vorticity quantification, the pressure drop between inlet and outlet, and the wall shear stress. Results demonstrated that detailed endocardial structures increase the degree of intra-ventricular pressure drop, decrease the wall shear stress, and disrupt the dominant vortex creating secondary small vortices. Increasingly turbulent blood flow was observed in the detailed RVs. Female RVs were less trabeculated and presented lower pressure drops than the males. In conclusion, neglecting endocardial structures in RV haemodynamic models may lead to inaccurate conclusions about the pressures, stresses, and blood flow behaviour in the cavity.

Left Ventricular Trabeculations Decrease the Wall Shear Stress and Increase the Intra-Ventricular Pressure Drop in CFD Simulations.

The aim of the present study is to characterize the hemodynamics of left ventricular (LV) geometries to examine the impact of trabeculae and papillary muscles (PMs) on blood flow using high performance computing (HPC). Five pairs of detailed and smoothed LV endocardium models were reconstructed from high-resolution magnetic resonance images (MRI) of ex-vivo human hearts. The detailed model of one LV pair is characterized only by the PMs and few big trabeculae, to represent state of art level of endocardial detail. The other four detailed models obtained include instead endocardial structures measuring ≥1 mm2 in cross-sectional area. The geometrical characterizations were done using computational fluid dynamics (CFD) simulations with rigid walls and both constant and transient flow inputs on the detailed and smoothed models for comparison. These simulations do not represent a clinical or physiological scenario, but a characterization of the interaction of endocardial structures with blood flow. Steady flow simulations were employed to quantify the pressure drop between the inlet and the outlet of the LVs and the wall shear stress (WSS). Coherent structures were analyzed using the Q-criterion for both constant and transient flow inputs. Our results show that trabeculae and PMs increase the intra-ventricular pressure drop, reduce the WSS and disrupt the dominant single vortex, usually present in the smoothed-endocardium models, generating secondary small vortices. Given that obtaining high resolution anatomical detail is challenging in-vivo, we propose that the effect of trabeculations can be incorporated into smoothed ventricular geometries by adding a porous layer along the LV endocardial wall. Results show that a porous layer of a thickness of 1.2·10-2 m with a porosity of 20 kg/m2 on the smoothed-endocardium ventricle models approximates the pressure drops, vorticities and WSS observed in the detailed models.

Attenuation of wave reflection by wave entrapment creates a "horizon effect" in the human aorta.

Wave reflection is thought to be important in the augmentation of blood pressure. However, identification of distal reflections sites remains unclear. One possible explanation for this is that wave reflection is predominately determined by an amalgamation of multiple proximal small reflections rather than large discrete reflections originating from the distal peripheries. In 19 subjects (age, 35-73 years), sensor-tipped intra-arterial wires were used to measure pressure and Doppler velocity at 10-cm intervals along the aorta, starting at the aortic root. Incident and reflected waves were identified and timings and magnitudes quantified using wave intensity analysis. Mean wave speed increased along the length of the aorta (proximal, 6.8±0.9 m/s; distal, 10.7±1.5 m/s). The incident wave was tracked moving along the aorta, taking 55±4 ms to travel from the aortic root to the distal aorta. However, the timing to the refection site distance did not differ between proximal and distal aortic measurement sites (proximal aorta, 48±5 ms versus distal aorta, 42±4 ms; P=0.3). We performed a second analysis using aortic waveforms in a nonlinear model of pulse-wave propagation. This demonstrated very similar results to those observed in vivo and also an exponential attenuation in reflection magnitude. There is no single dominant refection site in or near the distal aorta. Rather, there are multiple reflection sites along the aorta, for which the contributions are attenuated with distance. We hypothesize that rereflection of reflected waves leads to wave entrapment, preventing distal waves being seen in the proximal aorta.

Patient-specific modeling of dyssynchronous heart failure: a case study.

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.

The arterial reservoir pressure increases with aging and is the major determinant of the aortic augmentation index.

The augmentation index predicts cardiovascular mortality and is usually explained as a distally reflected wave adding to the forward wave generated by systole. We propose that the capacitative properties of the aorta (the arterial reservoir) also contribute significantly to the augmentation index and have calculated the contribution of the arterial reservoir, independently of wave reflection, and assessed how these contributions change with aging. In 15 subjects (aged 53 +/- 10 yr), we measured pressure and Doppler velocity simultaneously in the proximal aorta using intra-arterial wires. We calculated the components of augmentation pressure in two ways: 1) into forward and backward (reflected) components by established separation methods, and 2) using an approach that accounts for an additional reservoir component. When the reservoir was ignored, augmentation pressure (22.7 +/- 13.9 mmHg) comprised a small forward wave (peak pressure = 6.5 +/- 9.4 mmHg) and a larger backward wave (peak pressure = 16.2 +/- 7.6 mmHg). After we took account of the reservoir, the contribution to augmentation pressure of the backward wave was reduced by 64% to 5.8 +/- 4.4 mmHg (P < 0.001), forward pressure was negligible, and reservoir pressure was the largest component (peak pressure = 19.8 +/- 9.3 mmHg). With age, reservoir pressure increased progressively (9.9 mmHg/decade, r = 0.69, P < 0.001). In conclusion, the augmentation index is principally determined by aortic reservoir function and other elastic arteries and only to a minor extent by reflected waves. Reservoir function rather than wave reflection changes markedly with aging, which accounts for the age-related changes in the aortic pressure waveform.

Effect of transmurally heterogeneous myocyte excitation-contraction coupling on canine left ventricular electromechanics.

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.

Differences in cardiac microcirculatory wave patterns between the proximal left mainstem and proximal right coronary artery.

Despite having almost identical origins and similar perfusion pressures, the flow-velocity waveforms in the left and right coronary arteries are strikingly different. We hypothesized that pressure differences originating from the distal (microcirculatory) bed would account for the differences in the flow-velocity waveform. We used wave intensity analysis to separate and quantify proximal- and distal-originating pressures to study the differences in velocity waveforms. In 20 subjects with unobstructed coronary arteries, sensor-tipped intra-arterial wires were used to measure simultaneous pressure and Doppler velocity in the proximal left main stem (LMS) and proximal right coronary artery (RCA). Proximal- and distal-originating waves were separated using wave intensity analysis, and differences in waves were examined in relation to structural and anatomic differences between the two arteries. Diastolic flow velocity was lower in the RCA than in the LMS (35.1 +/- 21.4 vs. 56.4 +/- 32.5 cm/s, P < 0.002), and, consequently, the diastolic-to-systolic ratio of peak flow velocity in the RCA was significantly less than in the LMS (1.00 +/- 0.32 vs. 1.79 +/- 0.48, P < 0.001). This was due to a lower distal-originating suction wave (8.2 +/- 6.6 x 10(3) vs. 16.0 +/- 12.2 x 10(3) W.m(-2).s(-1), P < 0.01). The suction wave in the LMS correlated positively with left ventricular pressure (r = 0.6, P < 0.01) and in the RCA with estimated right ventricular systolic pressure (r = 0.7, P = 0.05) but not with the respective diameter in these arteries. In contrast to the LMS, where coronary flow velocity was predominantly diastolic, in the proximal RCA coronary flow velocity was similar in systole and diastole. This difference was due to a smaller distal-originating suction wave in the RCA, which can be explained by differences in elastance and pressure generated between right and left ventricles.

Reservoir-wave separation and wave intensity analysis applied to carotid arteries: a hybrid 1D model to understand haemodynamics.

Pressure waveforms measured at different locations in the cardiovascular system present a very similar diastolic decay. Previous work has shown the cardiovascular system can be modelled as a Windkessel and wave system. This concept has been extended to any arbitrary location in the cardiovascular system. We suggest that it is possible to calculate a time-varying reservoir pressure P(t) and a distance- and time-varying wave pressure p(x, t) by fitting an exponential function to the diastolic decay of the measured pressure P; defining that the measured pressure P(x, t) = P(t)+p(x, t). Velocity waveforms U can also be separated into its reservoir, U , and wave, u,components as U(x, t) = U (x, t) + u(x, t).In this study we explore the implications of applying are servoir-wave separation and wave intensity analysis techniques to understand the haemodynamics of in-vivo, noninvasive measurements of P and U in the carotid arteries of normal human subjects. Wave intensity analysis reveals a particular wave pattern where reflections can be estimated easily, but foremost, it shows that reflections are a lot smaller than previously thought.We suggest through the use of this model that the heart is the main wave generator of the cardiovascular system. The arterial system instead of impeding the flow, it stores it and distributes it throughout the arteries towards the tissue during diastole. There are some wave reflections, mainly during systole,that contribute to the changes in the pressure and velocity waveforms, however, they are small and are more evident as the measurements get further away from the ascending aorta.The application of wave intensity analysis to non-invasively measured data can provide a good insight on the physiology and the local and global properties of the cardiovascular system in health and disease in the clinical setting. This study shows preliminary results and the potential of the technique for analysing non-invasive measures, and could be particularly useful to understand and quantify the effects of therapeutic drugs in the cardiovascular system.

Pressure reservoir-wave separation applied to the coronary arterial data.

Pressure waveforms measured at different locations in the coronary arteries are similar. Previous work has shown the cardiovascular system can be modelled as a Windkessel and wave system. [1] We now suggest that simultaneously measured coronary pressure, P, and velocity, U, can be used to calculate the separated reservoir, P, and wave, p, pressure such that P = P + p. This separation enables us to deduce the resistance of the coronary microcirculation as it varies with time; resistance during systole being almost double the resistance during diastole. Wave intensity analysis is used to separate forward and backward waves. This requires knowledge of the local wave speed, c, which is determined from the data using the sum of squares technique. [2] Wave timings differ at different sites in the coronary tree, revealing the complexity of wave propagation. [3] Their mapping give a fairly complete picture of wave travel in the coronary network and insight into the physiology of coronary artery flow in health and disease.

A computational study on the influence of catheter-delivered intravascular probes on blood flow in a coronary artery model.

Catheter-delivered intravascular probes are widely used in clinical practice to measure coronary arterial velocity and pressure, but the artefactual effect of the probe on the variables being measured is not well characterised. A coronary artery was simulated with a 180 degrees curved tube 3mm in diameter and the effect of catheters of different diameters was modelled numerically under pulsatile flow conditions. The presence of a catheter increased pressure by 1.3-4.3 mmHg depending on its diameter, and reduced velocity-pressure phase-lag. For an ultrasound sample volume 5mm downstream from the probe tip, the underestimation in velocity measurement attributed to catheter blockage is approximately 15-21% for an average inlet velocity of 0.1m/s. The velocity measurement error is lower at higher mean flow velocity. Accuracy of clinical velocity measurements could be improved by moving the sample volume farther downstream from the probe tip, because the centrifugal pressure gradient intrinsic to the curvature promotes re-development of flow.

Evidence of a dominant backward-propagating "suction" wave responsible for diastolic coronary filling in humans, attenuated in left ventricular hypertrophy.

BACKGROUND: Coronary blood flow peaks in diastole when aortic blood pressure has fallen. Current models fail to completely explain this phenomenon. We present a new approach-using wave intensity analysis-to explain this phenomenon in normal subjects and to evaluate the effects of left ventricular hypertrophy (LVH). METHOD AND RESULTS: We measured simultaneous pressure and Doppler velocity with intracoronary wires in the left main stem, left anterior descending, and circumflex arteries of 20 subjects after a normal coronary arteriogram. Wave intensity analysis was used to identify and quantify individual pressure and velocity waves within the coronary artery circulation. A consistent pattern of 6 predominating waves was identified. Ninety-four percent of wave energy, accelerating blood forward along the coronary artery, came from 2 waves: first a pushing wave caused by left ventricular ejection-the dominant forward-traveling pushing wave; and later a suction wave caused by relief of myocardial microcirculatory compression-the dominant backward-traveling suction wave. The dominant backward-traveling suction wave (18.2+/-13.7 x 10(3) W m(-2)s(-1), 30%) was larger than the dominant forward-traveling pushing wave (14.3+/-17.6 x 10(3) W m(-2) s(-1), 22.3%, P =0.001) and was associated with a substantially larger increment in coronary blood flow velocity (0.51 versus 0.14 m/s, P <0.001). In LVH, the dominant backward-traveling suction wave percentage was significantly decreased (33.1% versus 26.9%, P =0.01) and inversely correlated with left ventricular septal wall thickness (r =-0.52, P <0.02). CONCLUSIONS: Six waves predominantly drive human coronary blood flow. Coronary flow peaks in diastole because of the dominance of a "suction" wave generated by myocardial microcirculatory decompression. This is significantly reduced in LVH.