Exploring the parameter space of a rabbit ventricular action potential model to investigate the effect of variation on action potential and calcium transients

Computational models for cardiomyocyte action potentials (AP) often make use of a large parameter set. This parameter set can contain some elements that are fitted to experimental data independently of any other element, some elements that are derived concurrently with other elements to match experimental data, and some elements that are derived purely from phenomenological fitting to produce the desired AP output. Furthermore, models can make use of several different data sets, not always derived for the same conditions or even the same species. It is consequently uncertain whether the parameter set for a given model is physiologically accurate. Furthermore, it is only recently that the possibility of degeneracy in parameter values in producing a given simulation output has started to be addressed. In this study, we examine the effects of varying two parameters (the L-type calcium current (I(CaL)) and the delayed rectifier potassium current (I(Ks))) in a computational model of a rabbit ventricular cardiomyocyte AP on both the membrane potential (V(m)) and calcium (Ca(2+)) transient. It will subsequently be determined if there is degeneracy in this model to these parameter values, which will have important implications on the stability of these models to cell-to-cell parameter variation, and also whether the current methodology for generating parameter values is flawed. The accuracy of AP duration (APD) as an indicator of AP shape will also be assessed.

Phenomenological modeling of cell-to-cell and beat-to-beat variability in isolated Guinea Pig ventricular myocytes

Experimental action potential (AP) recordings in isolated ventricular myoctes display significant temporal beat-to-beat variability in morphology and duration. Furthermore, significant cell-to-cell differences in AP also exist even for isolated cells originating from the same region of the same heart. However, current mathematical models of ventricular AP fail to replicate the temporal and cell-to-cell variability in AP observed experimentally. In this study, we propose a novel mathematical framework for the development of phenomenological AP models capable of capturing cell-to-cell and temporal variabilty in cardiac APs. A novel stochastic phenomenological model of the AP is developed, based on the deterministic Bueno-Orovio/Fentonmodel. Experimental recordings of AP are fit to the model to produce AP models of individual cells from the apex and the base of the guinea-pig ventricles. Our results show that the phenomenological model is able to capture the considerable differences in AP recorded from isolated cells originating from the location. We demonstrate the closeness of fit to the available experimental data which may be achieved using a phenomenological model, and also demonstrate the ability of the stochastic form of the model to capture the observed beat-to-beat variablity in action potential duration.

Frank-starling control of a left ventricular assist device

A physiological control system was developed for a rotary left ventricular assist device (LVAD) in which the target pump flow rate (LVADQ) was set as a function of left atrial pressure (LAP), mimicking the Frank-Starling mechanism. The control strategy was implemented using linear PID control and was evaluated in a pulsatile mock circulation loop using a prototyped centrifugal pump by varying pulmonary vascular resistance to alter venous return. The control strategy automatically varied pump speed (2460 to 1740 to 2700 RPM) in response to a decrease and subsequent increase in venous return. In contrast, a fixed-speed pump caused a simulated ventricular suction event during low venous return and higher ventricular volumes during high venous return. The preload sensitivity was increased from 0.011 L/min/mmHg in fixed speed mode to 0.47L/min/mmHg, a value similar to that of the native healthy heart. The sensitivity varied automatically to maintain the LAP and LVADQ within a predefined zone. This control strategy requires the implantation of a pressure sensor in the left atrium and a flow sensor around the outflow cannula of the LVAD. However, appropriate pressure sensor technology is not yet commercially available and so an alternative measure of preload such as pulsatility of pump signals should be investigated.

Evaluation of inflow cannulation site for implantation of right-sided rotary ventricular assist device

Right heart dysfunction is one of the most serious complications following implantation of a left ventricular assist device (LVAD), often leading to the requirement for short or long term right ventricular support (RVAD). The inflow cannulation site induces major haemodynamic changes and so there is a need to optimize the site used depending on the patient's condition. Therefore, this study evaluated and compared the haemodynamic influence of right atrial (RAC) and right ventricular (RVC) inflow cannulation sites. An in-vitro, variable heart failure, mock circulation loop was used to compare RAC and RVC in mild and severe biventricular heart failure (BHF) conditions. In the severe BHF condition, higher ventricular ejection fraction (RAC: 13.6%, RVC: 32.7%) and thus improved heart chamber and RVAD washout was observed with RVC, which suggested this strategy might be preferable for long term support (ie. bridge to transplant or destination therapy) to reduce the risk of thrombus formation. In the mild BHF condition, higher pulmonary valve flow (RAC: 3.33 L/min, RVC: 1.97 L/min) and lower right ventricular stroke work (RAC: 0.10 W, RVC: 0.13 W) and volumes were recorded with RAC. These results indicate an improved potential for myocardial recovery, thus RAC should be chosen in this condition. This in-vitro study suggests that RVAD inflow cannulation site should be chosen on a patient-specific basis with a view to the support strategy to promote myocardial recovery or reduce the risk of long-term complications.

Left ventricular volume and valve gradient analysis using an Apple Macintosh computer

A description of a computer program to analyse cine angiograms of the heart and pressure waveforms to calculate valve gradients.

Quantitative genetic modeling of lateral ventricular shape and volume using multi-atlas fluid image alignment in twins

Despite substantial progress in measuring the 3D profile of anatomical variations in the human brain, their genetic and environmental causes remain enigmatic. We developed an automated system to identify and map genetic and environmental effects on brain structure in large brain MRI databases . We applied our multi-template segmentation approach ("Multi-Atlas Fluid Image Alignment") to fluidly propagate hand-labeled parameterized surface meshes into 116 scans of twins (60 identical, 56 fraternal), labeling the lateral ventricles. Mesh surfaces were averaged within subjects to minimize segmentation error. We fitted quantitative genetic models at each of 30,000 surface points to measure the proportion of shape variance attributable to (1) genetic differences among subjects, (2) environmental influences unique to each individual, and (3) shared environmental effects. Surface-based statistical maps revealed 3D heritability patterns, and their significance, with and without adjustments for global brain scale. These maps visualized detailed profiles of environmental versus genetic influences on the brain, extending genetic models to spatially detailed, automatically computed, 3D maps.

Mapping genetic influences on ventricular structure in twins

Despite substantial progress in measuring the anatomical and functional variability of the human brain, little is known about the genetic and environmental causes of these variations. Here we developed an automated system to visualize genetic and environmental effects on brain structure in large brain MRI databases. We applied our multi-template segmentation approach termed "Multi-Atlas Fluid Image Alignment" to fluidly propagate hand-labeled parameterized surface meshes, labeling the lateral ventricles, in 3D volumetric MRI scans of 76 identical (monozygotic, MZ) twins (38 pairs; mean age = 24.6 (SD = 1.7)); and 56 same-sex fraternal (dizygotic, DZ) twins (28 pairs; mean age = 23.0 (SD = 1.8)), scanned as part of a 5-year research study that will eventually study over 1000 subjects. Mesh surfaces were averaged within subjects to minimize segmentation error. We fitted quantitative genetic models at each of 30,000 surface points to measure the proportion of shape variance attributable to (1) genetic differences among subjects, (2) environmental influences unique to each individual, and (3) shared environmental effects. Surface-based statistical maps, derived from path analysis, revealed patterns of heritability, and their significance, in 3D. Path coefficients for the 'ACE' model that best fitted the data indicated significant contributions from genetic factors (A = 7.3%), common environment (C = 38.9%) and unique environment (E = 53.8%) to lateral ventricular volume. Earlier-maturing occipital horn regions may also be more genetically influenced than later-maturing frontal regions. Maps visualized spatially-varying profiles of environmental versus genetic influences. The approach shows promise for automatically measuring gene-environment effects in large image databases.

Automated ventricular mapping with multi-atlas fluid image alignment reveals genetic effects in Alzheimer's disease

We developed and validated a new method to create automated 3D parametric surface models of the lateral ventricles in brain MRI scans, providing an efficient approach to monitor degenerative disease in clinical studies and drug trials. First, we used a set of parameterized surfaces to represent the ventricles in four subjects' manually labeled brain MRI scans (atlases). We fluidly registered each atlas and mesh model to MRIs from 17 Alzheimer's disease (AD) patients and 13 age- and gender-matched healthy elderly control subjects, and 18 asymptomatic ApoE4-carriers and 18 age- and gender-matched non-carriers. We examined genotyped healthy subjects with the goal of detecting subtle effects of a gene that confers heightened risk for Alzheimer's disease. We averaged the meshes extracted for each 3D MR data set, and combined the automated segmentations with a radial mapping approach to localize ventricular shape differences in patients. Validation experiments comparing automated and expert manual segmentations showed that (1) the Hausdorff labeling error rapidly decreased, and (2) the power to detect disease- and gene-related alterations improved, as the number of atlases, N, was increased from 1 to 9. In surface-based statistical maps, we detected more widespread and intense anatomical deficits as we increased the number of atlases. We formulated a statistical stopping criterion to determine the optimal number of atlases to use. Healthy ApoE4-carriers and those with AD showed local ventricular abnormalities. This high-throughput method for morphometric studies further motivates the combination of genetic and neuroimaging strategies in predicting AD progression and treatment response. © 2007 Elsevier Inc. All rights reserved.

Mapping hippocampal and ventricular change in Alzheimer disease

We developed an anatomical mapping technique to detect hippocampal and ventricular changes in Alzheimer disease (AD). The resulting maps are sensitive to longitudinal changes in brain structure as the disease progresses. An anatomical surface modeling approach was combined with surface-based statistics to visualize the region and rate of atrophy in serial MRI scans and isolate where these changes link with cognitive decline. Fifty-two high-resolution MRI scans were acquired from 12 AD patients (age: 68.4 ± 1.9 years) and 14 matched controls (age: 71.4 ± 0.9 years), each scanned twice (2.1 ± 0.4 years apart). 3D parametric mesh models of the hippocampus and temporal horns were created in sequential scans and averaged across subjects to identify systematic patterns of atrophy. As an index of radial atrophy, 3D distance fields were generated relating each anatomical surface point to a medial curve threading down the medial axis of each structure. Hippocampal atrophic rates and ventricular expansion were assessed statistically using surface-based permutation testing and were faster in AD than in controls. Using color-coded maps and video sequences, these changes were visualized as they progressed anatomically over time. Additional maps localized regions where atrophic changes linked with cognitive decline. Temporal horn expansion maps were more sensitive to AD progression than maps of hippocampal atrophy, but both maps correlated with clinical deterioration. These quantitative, dynamic visualizations of hippocampal atrophy and ventricular expansion rates in aging and AD may provide a promising measure to track AD progression in drug trials.

Corrigendum to "Mapping hippocampal and ventricular change in Alzheimer disease" [NeuroImage 22 (2004) 1754-1766] (DOI:10.1016/j.neuroimage.2004.03.040)

We recently noticed an error in the demographic data in this article. The validity of the findings and the conclusions of the paper is not affected. However, there is an error in the reported sample size and in the means and standard deviations of the subjects’ ages and MMSE scores. We would like to correct this error, which came to light when we were re-analyzing the data for a meta-analysis. The error occurred because an older version of a spreadsheet was incorrectly used when reporting the sample composition. Instead of examining 12 Alzheimer's disease patients and 14 healthy elderly controls, we in fact examined 17 Alzheimer’s disease patients and 14 healthy elderly controls. All maps and morphometric data reported in the paper are correct, except that the sample size was in fact slightly higher than that originally reported, and the maps computed in the paper were based on the larger sample (which included five more subjects in the Alzheimer’s disease group). All of the maps and figures in the paper are correct, and the conclusions of the paper are unchanged. We apologize for this error, which falls under the sole responsibility of the first author. The corrected demographic information appears below.

Feasibility of perflutren microsphere contrast transthoracic echocardiography in the visualization of ventricular endocardium during venovenous extracorporeal membrane oxygenation in a validated ovine model

Background: Transthoracic echocardiography (TTE) during extra corporeal membrane oxygenation (ECMO) is important but can be technically challenging. Contrast-specific TTE can improve imaging in suboptimal studies. These contrast microspheres are hydrodynamically labile structures. This study assessed the feasibility of contrast echocardiography (CE) during venovenous (VV) ECMO in a validated ovine model. Method: Twenty-four sheep were commenced on VV ECMO. Parasternal long-axis (Plax) and short-axis (Psax) views were obtained pre- and postcontrast while on VV ECMO. Endocardial definition scores (EDS) per segment were graded: 1 = good, 2 = suboptimal 3 = not seen. Endocardial border definition score index (EBDSI) was calculated for each view. Endocardial length (EL) in the Plax view for the left ventricle (LV) and right ventricle (RV) was measured. Results: Summation EDS data for the LV and RV for unenhanced TTE (UE) versus CE TTE imaging: EDS 1 = 289 versus 346, EDS 2 = 38 versus 10, EDS 3 = 33 versus 4, respectively. Wilcoxon matched-pairs rank-sign tests showed a significant ranking difference (improvement) pre- and postcontrast for the LV (P < 0.0001), RV (P < 0.0001) and combined ventricular data (P < 0.0001). EBDSI for CE TTE was significantly lower than UE TTE for the LV (1.05 ± 0.17 vs. 1.22 ± 0.38, P = 0.0004) and RV (1.06 ± 0.22 vs. 1.42 ± 0.47, P = 0.0.0006) respectively. Visualized EL was significantly longer in CE versus UE for both the LV (58.6 ± 11.0 mm vs. 47.4 ± 11.7 mm, P < 0.0001) and the RV (52.3 ± 8.6 mm vs. 36.0 ± 13.1 mm, P < 0.0001), respectively. Conclusions: Despite exposure to destructive hydrodynamic forces, CE is a feasible technique in an ovine ECMO model. CE results in significantly improved EDS and increased EL.