Bridging experiments, models and simulations : an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatio-temporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and scepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace and refine animal experiments. A fundamental requirement to fulfil these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations between experiments, models and simulations in cardiac electrophysiology. We describe the processes, data and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. Validation must therefore take into account the complex interplay between models, simulations and experiments. Key points for developing strategies for validation are: 1) understanding sources of bio-variability is crucial to the comparison between simulation and experimental results; 2) robustness of techniques and tools is a pre-requisite to conducting physiological investigations using the MSE system; 3) definition and adoption of standards facilitates interoperability of experiments, models and simulations; 4) physiological validation must be understood as an iterative process that defines the specific aspects of electrophysiology the MSE system targets, and is driven by advancements in experimental and computational methods and the combination of both.

Comparison of ocular modulation transfer function determined by a ray-tracing aberrometer and a double-pass system in early cataract patients

BACKGROUND: The evaluation of retinal image quality in cataract eyes has gained importance and the clinical modulation transfer functions (MTF) can obtained by aberrometer and double pass (DP) system. This study aimed to compare MTF derived from a ray tracing aberrometer and a DP system in early cataractous and normal eyes. METHODS: There were 128 subjects with 61 control eyes and 67 eyes with early cataract defined according to the Lens Opacities Classification System III. A laser ray-tracing wavefront aberrometer (iTrace) and a double pass (DP) system (OQAS) assessed ocular MTF for 6.0 mm pupil diameters following dilation. Areas under the MTF (AUMTF) and their correlations were analyzed. Stepwise multiple regression analysis assessed factors affecting the differences between iTrace- and OQAS-derived AUMTF for the early cataract group. RESULTS: For both early cataract and control groups, iTrace-derived MTFs were higher than OQAS-derived MTFs across a range of spatial frequencies (P < 0.01). No significant difference between the two groups occurred for iTrace-derived AUMTF, but the early cataract group had significantly smaller OQAS-derived AUMTF than did the control group (P < 0.01). AUMTF determined from both the techniques demonstrated significant correlations with nuclear opacities, higher-order aberrations (HOAs), visual acuity, and contrast sensitivity functions, while the OQAS-derived AUMTF also demonstrated significant correlations with age and cortical opacity grade. The factors significantly affecting the difference between iTrace and OQAS AUMTF were root-mean-squared HOAs (standardized beta coefficient = -0.63, P < 0.01) and age (standardized beta coefficient = 0.26, P < 0.01). CONCLUSIONS: MTFs determined from a iTrace and a DP system (OQAS) differ significantly in early cataractous and normal subjects. Correlations with visual performance were higher for the DP system. OQAS-derived MTF may be useful as an indicator of visual performance in early cataract eyes.

A method to evaluate the response of the coronary circulation of perfused rat heart to adenosine

Exogenous adenosine causes a monophasic dilation of the coronary vessels in paced, perfused rat heart preparations. Because levels of endogenous adenosine in paced hearts may mask the presence of high potency adenosine receptors, we have developed a method to measure coronary vascular responses in a potassium-arrested heart. Hearts from adult male, Wistar rats were perfused at a constant flow rate of 10 mL/min in the nonrecirculating, Langendorff mode, using Krebs-Henseleit buffer. After 30 min, coronary perfusion pressure was 44 +/- 1 mmHg (mean +/- SEM). Hearts were then perfused with a modified Krebs-Henseleit buffer containing 35 mM potassium. Coronary perfusion pressure increased by 84 +/- 3 mmHg. Adenosine-induced reductions in coronary perfusion pressure were expressed as a percentage of the maximal increase in pressure produced by modified Krebs-Henseleit buffer from the equilibration level. A concentration-response curve for adenosine (n = 6) was biphasic and best described by the presence of two adenosine receptors, with negative log EC50 values of 8.8 +/- 0.3 and 4.3 +/- 0.1, representing 29 +/- 3 and 71 +/- 3%, respectively, of the observed response. Interstitial adenosine sampled by microdialysis during potassium arrest was 25% of the concentration found in paced hearts. Endogenous adenosine in nonarrested hearts may obscure the biphasic response of the coronary vessels to adenosine.

Automated 3D mapping & shape analysis of the lateral ventricles via fluid registration of multiple surface-based atlases

We developed and validated a new method to create automated 3D parametric surface models of the lateral ventricles, designed for monitoring degenerative disease effects in clinical neuroscience studies and drug trials. First we used a set of parameterized surfaces to represent the ventricles in a manually labeled set of 9 subjects' MRIs (atlases). We fluidly registered each of these atlases and mesh models to a set of MRIs from 12 Alzheimer's disease (AD) patients and 14 matched healthy elderly subjects, and we averaged the resulting meshes for each of these images. Validation experiments on expert segmentations showed that (1) the Hausdorff labeling error rapidly decreased, and (2) the power to detect disease-related alterations monotonically improved as the number of atlases, N, was increased from 1 to 9. We then combined the segmentations with a radial mapping approach to localize ventricular shape differences in patients. In surface-based statistical maps, we detected more widespread and intense anatomical deficits as we increased the number of atlases, and we formulated a statistical stopping criterion to determine the optimal value of N. Anterior horn anomalies in Alzheimer's patients were only detected with the multi-atlas segmentation, which clearly outperformed the standard single-atlas approach.

4D deformation modeling of cortical disease progression in Alzheimer's dementia

This work describes the development of a model of cerebral atrophic changes associated with the progression of Alzheimer's disease (AD). Linear registration, region-of-interest analysis, and voxel-based morphometry methods have all been employed to elucidate the changes observed at discrete intervals during a disease process. In addition to describing the nature of the changes, modeling disease-related changes via deformations can also provide information on temporal characteristics. In order to continuously model changes associated with AD, deformation maps from 21 patients were averaged across a novel z-score disease progression dimension based on Mini Mental State Examination (MMSE) scores. The resulting deformation maps are presented via three metrics: local volume loss (atrophy), volume (CSF) increase, and translation (interpreted as representing collapse of cortical structures). Inspection of the maps revealed significant perturbations in the deformation fields corresponding to the entorhinal cortex (EC) and hippocampus, orbitofrontal and parietal cortex, and regions surrounding the sulci and ventricular spaces, with earlier changes predominantly lateralized to the left hemisphere. These changes are consistent with results from post-mortem studies of AD.

MR image-based measurement of rates of change in volumes of brain structures. Part II: Application to a study of Alzheimer's disease and normal aging

We present global and regional rates of brain atrophy measured on serially acquired Tl-weighted brain MR images for a group of Alzheimer's disease (AD) patients and age-matched normal control (NC) subjects using the analysis procedure described in Part I. Three rates of brain atrophy: the rate of atrophy in the cerebrum, the rate of lateral ventricular enlargement and the rate of atrophy in the region of temporal lobes, were evaluated for 14 AD patients and 14 age-matched NC subjects. All three rates showed significant differences between the two groups. However, the greatest separation of the two groups was obtained when the regional rates were combined. This application has demonstrated that rates of brain atrophy, especially in specific regions of the brain, based on MR images can provide sensitive measures for evaluating the progression of AD. These measures will be useful for the evaluation of therapeutic effects of novel therapies for AD.

Letter: Brain natriuretic peptide is a potentially useful screening tool for the detection of cardiovascular disease in patients with rheumatoid arthritis

Patients with rheumatoid arthritis (RA) have a significantly higher risk of coronary heart disease, despite being less likely to report symptoms of angina, and are more likely to experience unrecognised myocardial infarction and sudden cardiac death than non-RA controls.1 Furthermore, left ventricular diastolic dysfunction has been described in up to 40% of patients with RA.2...

Flow evaluation and hemolysis analysis of BVAD centrifugal blood pump by computational fluids dynamics

Computational fluid dynamics (CFD) and particle image velocimetry (PIV) are commonly used techniques to evaluate the flow characteristics in the development stage of blood pumps. CFD technique allows rapid change to pump parameters to optimize the pump performance without having to construct a costly prototype model. These techniques are used in the construction of a bi-ventricular assist device (BVAD) which combines the functions of LVAD and RVAD in a compact unit. The BVAD construction consists of two separate chambers with similar impellers, volutes, inlet and output sections. To achieve the required flow characteristics of an average flow rate of 5 l/min and different pressure heads (left – 100mmHg and right – 20mmHg), the impellers were set at different rotating speeds. From the CFD results, a six-blade impeller design was adopted for the development of the BVAD. It was also observed that the fluid can flow smoothly through the pump with minimum shear stress and area of stagnation which are related to haemolysis and thrombosis. Based on the compatible Reynolds number the flow through the model was calculated for the left and the right pumps. As it was not possible to have both the left and right chambers in the experimental model, the left and right pumps were tested separately.

Carvedilol induces greater control of β2- than β1-adrenoceptor-mediated inotropic and lusitropic effects by PDE3, while PDE4 has no effect in human failing myocardium

The beta-blockers carvedilol and metoprolol provide important therapeutic strategies for heart failure treatment. Therapy with metoprolol facilitates the control by phosphodiesterase PDE3, but not PDE4, of inotropic effects of catecholamines in human failing ventricle. However, it is not known whether carvedilol has the same effect. We investigated whether the PDE3-selective inhibitor cilostamide (0.3 mu M) or PDE4-selective inhibitor rolipram (1 mu M) modified the positive inotropic and lusitropic effects of catecholamines in ventricular myocardium of heart failure patients treated with carvedilol. Right ventricular trabeculae from explanted hearts of nine carvedilol-treated patients with terminal heart failure were paced to contract at 1 Hz. The effects of (-)-noradrenaline, mediated through beta(1)-adrenoceptors (beta(2)-adrenoceptors blocked with ICI118551), and (-)-adrenaline, mediated through beta(2)-adrenoceptors (beta(1)-adrenoceptors blocked with CGP20712A), were assessed in the absence and presence of the PDE inhibitors. The inotropic potency, estimated from -logEC(50)s, was unchanged for (-)-noradrenaline but decreased 16-fold for (-)-adrenaline in carvedilol-treated compared to non-beta-blocker-treated patients, consistent with the previously reported beta(2)-adrenoceptor-selectivity of carvedilol. Cilostamide caused 2- to 3-fold and 10- to 35-fold potentiations of the inotropic and lusitropic effects of (-)-noradrenaline and (-)-adrenaline, respectively, in trabeculae from carvedilol-treated patients. Rolipram did not affect the inotropic and lusitropic potencies of (-)-noradrenaline or (-)-adrenaline. Treatment of heart failure patients with carvedilol induces PDE3 to selectively control the positive inotropic and lusitropic effects mediated through ventricular beta(2)-adrenoceptors compared to beta(1)-adrenoceptors. The beta(2)-adrenoceptor-selectivity of carvedilol may provide protection against beta(2)-adrenoceptor-mediated ventricular overstimulation in PDE3 inhibitor-treated patients. PDE4 does not control beta(1)- and beta(2)-adrenoceptor-mediated inotropic and lusitropic effects in carvedilol-treated patients.

Robust iris recognition using decision fusion and degradation modelling

This thesis investigates the use of fusion techniques and mathematical modelling to increase the robustness of iris recognition systems against iris image quality degradation, pupil size changes and partial occlusion. The proposed techniques improve recognition accuracy and enhance security. They can be further developed for better iris recognition in less constrained environments that do not require user cooperation. A framework to analyse the consistency of different regions of the iris is also developed. This can be applied to improve recognition systems using partial iris images, and cancelable biometric signatures or biometric based cryptography for privacy protection.

A biomechanical approach to iris normalization

The richness of the iris texture and its variability across individuals make it a useful biometric trait for personal authentication. One of the key stages in classical iris recognition is the normalization process, where the annular iris region is mapped to a dimensionless pseudo-polar coordinate system. This process results in a rectangular structure that can be used to compensate for differences in scale and variations in pupil size. Most iris recognition methods in the literature adopt linear sampling in the radial and angular directions when performing iris normalization. In this paper, a biomechanical model of the iris is used to define a novel nonlinear normalization scheme that improves iris recognition accuracy under different degrees of pupil dilation. The proposed biomechanical model is used to predict the radial displacement of any point in the iris at a given dilation level, and this information is incorporated in the normalization process. Experimental results on the WVU pupil light reflex database (WVU-PLR) indicate the efficacy of the proposed technique, especially when matching iris images with large differences in pupil size.

Comparison between computational models of abdominal aortic aneurysm

Objective: To compare the differences in the hemodynamic parameters of abdominal aortic aneurysm (AAA) between fluid-structure interaction model (FSIM) and fluid-only model (FM), so as to discuss their application in the research of AAA. Methods: An idealized AAA model was created based on patient-specific AAA data. In FM, the flow, pressure and wall shear stress (WSS) were computed using finite volume method. In FSIM, an Arbitrary Lagrangian-Eulerian algorithm was used to solve the flow in a continuously deforming geometry. The hemodynamic parameters of both models were obtained for discussion. Results: Under the same inlet velocity, there were only two symmetrical vortexes in the AAA dilation area for FSIM. In contrast, four recirculation areas existed in FM; two were main vortexes and the other two were secondary flow, which were located between the main recirculation area and the arterial wall. Six local pressure concentrations occurred in the distal end of AAA and the recirculation area for FM. However, there were only two local pressure concentrations in FSIM. The vortex center of the recirculation area in FSIM was much more close to the distal end of AAA and the area was much larger because of AAA expansion. Four extreme values of WSS existed at the proximal of AAA, the point of boundary layer separation, the point of flow reattachment and the distal end of AAA, respectively, in both FM and FSIM. The maximum wall stress and the largest wall deformation were both located at the proximal and distal end of AAA. Conclusions: The number and center of the recirculation area for both models are different, while the change of vortex is closely associated with the AAA growth. The largest WSS of FSIM is 36% smaller than that of FM. Both the maximum wall stress and largest wall displacement shall increase with the outlet pressure increasing. FSIM needs to be considered for studying the relationship between AAA growth and shear stress.

Exercise and vascular function in child obesity: A meta-analysis

CONTEXT: Conduit artery flow-mediated dilation (FMD) is a noninvasive index of preclinical atherosclerosis in humans. Exercise interventions can improve FMD in both healthy and clinical populations. OBJECTIVE: This systematic review and meta-analysis aimed to summarize the effect of exercise training on FMD in overweight and obese children and adolescents as well as investigate the role of cardiorespiratory fitness (peak oxygen consumption [Vo2peak]) on effects observed. DATA SOURCES: PubMed, Medline, Embase, and Cinahl databases were searched from the earliest available date to February 2015. STUDY SELECTION: Studies of children and/or adolescents who were overweight or obese were included. DATA EXTRACTION: Standardized data extraction forms were used for patient and intervention characteristics, control/comparator groups, and key outcomes. Procedural quality of the studies was assessed using a modified version of the Physiotherapy Evidence Base Database scale. RESULTS: A meta-analysis involving 219 participants compared the mean difference of pre- versus postintervention vascular function (FMD) and Vo2peak between an exercise training intervention and a control condition. There was a significantly greater improvement in FMD (mean difference 1.54%, P < .05) and Vo2peak (mean difference 3.64 mL/kg/min, P < .05) after exercise training compared with controls. LIMITATIONS: Given the diversity of exercise prescriptions, participant characteristics, and FMD measurement protocols, varying FMD effect size was noted between trials. CONCLUSIONS: Exercise training improves vascular function in overweight and obese children, as indicated by enhanced FMD. Further research is required to establish the optimum exercise program for maintenance of healthy vascular function in this at-risk pediatric population.