The Virtual Physiological Human: Ten Years After

Biomedical research and clinical practice are struggling to cope with the growing complexity that the progress of health care involves. The most challenging diseases, those with the largest socioeconomic impact (cardiovascular conditions; musculoskeletal conditions; cancer; metabolic, immunity, and neurodegenerative conditions), are all characterized by a complex genotype–phenotype interaction and by a “systemic” nature that poses a challenge to the traditional reductionist approach. In 2005 a small group of researchers discussed how the vision of computational physiology promoted by the Physiome Project could be translated into clinical practice and formally proposed the term Virtual Physiological Human. Our knowledge about these diseases is fragmentary, as it is associated with molecular and cellular processes on the one hand and with tissue and organ phenotype changes (related to clinical symptoms of disease conditions) on the other. The problem could be solved if we could capture all these fragments of knowledge into predictive models and then compose them into hypermodels that help us tame the complexity that such systemic behavior involves. In 2005 this was simply not possible—the necessary methods and technologies were not available. Now, 10 years later, it seems the right time to reflect on the original vision, the results achieved so far, and what remains to be done

Quantitative comparison of myocardial fiber structure between mice, rabbit, and sheep using diffusion tensor cardiovascular magnetic resonance

<p>Abstract</p> <p>Background</p> <p>Accurate interpretations of cardiac functions require precise structural models of the myocardium, but the latter is not available always and for all species. Although scaling or substitution of myocardial fiber information from alternate species has been used in cardiac functional modeling, the validity of such practice has not been tested.</p> <p>Methods</p> <p>Fixed mouse (n = 10), rabbit (n = 6), and sheep (n = 5) hearts underwent diffusion tensor imaging (DTI). The myocardial structures in terms of the left ventricular fiber orientation helix angle index were quantitatively compared between the mouse rabbit and sheep hearts.</p> <p>Results</p> <p>The results show that significant fiber structural differences exist between any two of the three species. Specifically, the subepicardial fiber orientation, and the transmural range and linearity of fiber helix angles are significantly different between the mouse and either rabbit or sheep. Additionally, a significant difference was found between the transmural helix angle range between the rabbit and sheep. Across different circumferential regions of the heart, the fiber orientation was not found to be significantly different.</p> <p>Conclusions</p> <p>The current study indicates that myocardial structural differences exist between different size hearts. An immediate implication of the present findings for myocardial structural or functional modeling studies is that caution must be exercised when extrapolating myocardial structures from one species to another.</p

Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

Background: Cardiovascular magnetic resonance (CMR) can through the two methods 3D FLASH and diffusion tensor imaging (DTI) give complementary information on the local orientations of cardiomyocytes and their laminar arrays. Methods: Eight explanted rat hearts were perfused with Gd-DTPA contrast agent and fixative and imaged in a 9.4T magnet by two types of acquisition: 3D fast low angle shot (FLASH) imaging, voxels 50 × 50 × 50 μm, and 3D spin echo DTI with monopolar diffusion gradients of 3.6 ms duration at 11.5 ms separation, voxels 200 × 200 × 200 μm. The sensitivity of each approach to imaging parameters was explored. Results:The FLASH data showed laminar alignments of voxels with high signal, in keeping with the presumed predominance of contrast in the interstices between sheetlets. It was analysed, using structure-tensor (ST) analysis, to determine the most (v 1 ST ), intermediate (v 2 ST ) and least (v 3 ST ) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e 1 DTI ), intermediate (e 2 DTI ) and least (e 3 DTI ) orthogonal eigenvectors of extent of diffusion. The correspondence between the FLASH and DTI methods was measured and appraised. The most extended direction of FLASH signal (v 1 ST ) agreed well with that of diffusion (e 1 DTI ) throughout the left ventricle (representative discrepancy in the septum of 13.3 ± 6.7°: median ± absolute deviation) and both were in keeping with the expected local orientations of the long-axis of cardiomyocytes. However, the orientation of the least directions of FLASH signal continuity (v 3 ST ) and diffusion (e 3 ST ) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v 3 ST ) and DTI (e 3 DTI ) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v 3 ST and the directly measured FLASH laminar array normal was of 9 ± 7° for the lateral wall and 7 ± 9° for the septum (median ± inter quartile range), and for DTI the discrepancy between the calculated v 3 DTI and the directly measured FLASH laminar array normal was 22 ± 14° and 61 ± 53.4°. DTI was relatively insensitive to the number of diffusion directions and to time up to 72 hours post fixation, but was moderately affected by b-value (which was scaled by modifying diffusion gradient pulse strength with fixed gradient pulse separation). Optimal DTI parameters were b = 1000 mm/s2 and 12 diffusion directions. FLASH acquisitions were relatively insensitive to the image processing parameters explored. Conclusions: We show that ST analysis of FLASH is a useful and accurate tool in the measurement of cardiac microstructure. While both FLASH and the DTI approaches appear promising for mapping of the alignments of myocytes throughout myocardium, marked discrepancies between the cross myocyte anisotropies deduced from each method call for consideration of their respective limitations

Heterogeneity of fractional anisotropy and mean diffusivity measurements by in vivo diffusion tensor imaging in normal human hearts

Background: Cardiac diffusion tensor imaging (cDTI) by cardiovascular magnetic resonance has the potential to assess microstructural changes through measures of fractional anisotropy (FA) and mean diffusivity (MD). However, normal variation in regional and transmural FA and MD is not well described. Methods: Twenty normal subjects were scanned using an optimised cDTI sequence at 3T in systole. FA and MD were quantified in 3 transmural layers and 4 regional myocardial walls. Results: FA was higher in the mesocardium (0.46 ±0.04) than the endocardium (0.40 ±0.04, p≤0.001) and epicardium (0.39 ±0.04, p≤0.001). On regional analysis, the FA in the septum was greater than the lateral wall (0.44 ±0.03 vs 0.40 ±0.05 p = 0.04). There was a transmural gradient in MD increasing towards the endocardium (epicardium 0.87 ±0.07 vs endocardium 0.91 ±0.08×10-3 mm2/s, p = 0.04). With the lateral wall (0.87 ± 0.08×10-3 mm2/s) as the reference, the MD was higher in the anterior wall (0.92 ±0.08×10-3 mm2/s, p = 0.016) and septum (0.92 ±0.07×10-3 mm2/s, p = 0.028). Transmurally the signal to noise ratio (SNR) was greatest in the mesocardium (14.5 ±2.5 vs endocardium 13.1 ±2.2, p<0.001; vs epicardium 12.0 ± 2.4, p<0.001) and regionally in the septum (16.0 ±3.4 vs lateral wall 11.5 ± 1.5, p<0.001). Transmural analysis suggested a relative reduction in the rate of change in helical angle (HA) within the mesocardium. Conclusions: In vivo FA and MD measurements in normal human heart are heterogeneous, varying significantly transmurally and regionally. Contributors to this heterogeneity are many, complex and interactive, but include SNR, variations in cardiac microstructure, partial volume effects and strain. These data indicate that the potential clinical use of FA and MD would require measurement standardisation by myocardial region and layer, unless pathological changes substantially exceed the normal variation identified

Non-uniform recovery of left ventricular transmural mechanics in ST-segment elevation myocardial infarction

<p>Abstract</p> <p>Background</p> <p>After a transient ischemic episode, the subendocardial region is more severely injured than outer subepicardial layers and may regain a proportionately greater degree of mechanical function in the longitudinal direction. We sought to explore left ventricular (LV) transmural mechanics in patients with ST-segment elevation myocardial infarction (STEMI) for determining the mechanism underlying recovery of global LV function after primary percutaneous coronary intervention (PCI).</p> <p>Methods</p> <p>A total of 42 patients (62 ± 11 years old, 71% male) with a first STEMI underwent serial assessments of LV longitudinal, circumferential and radial strains (LS, CS and RS) by selective tracking of subendocardial and subepicardial regions within 48 hours and a median of 5 months after PCI. LV mechanical parameters were compared with sixteen age and gender matched normal controls.</p> <p>Results</p> <p>In comparison with controls, endocardial and epicardial LS were markedly attenuated at 48 hours following PCI (P < 0.001). An improvement in LV ejection fraction (EF > 5%) following PCI was seen in 24 (57%) patients and was associated with improvement in endocardial and epicardial LS (P < 0.001 and P = 0.003, respectively) and endocardial CS (P = 0.01). Radial strain and wall motion score index, however, remained persistently abnormal. The change in endocardial LS (OR 1.2, 95% CI 1.03 to 1.42, P = 0.01) and the change in epicardial LS (OR 1.2, 95% 1.03 to 1.46, P = 0.02) were significantly associated with the improvement in LVEF, independent of the location of STEMI and the presence of underlying multivessel disease.</p> <p>Conclusions</p> <p>In patients with STEMI treated by PCI, the recovery of LV subendocardial shortening strain seen in the longitudinal direction underlies the improvement in LV global function despite persistent abnormalities in radial mechanics and wall motion score index.</p

The Role of the Frank–Starling Law in the Transduction of Cellular Work to Whole Organ Pump Function: A Computational Modeling Analysis

We have developed a multi-scale biophysical electromechanics model of the rat left ventricle at room temperature. This model has been applied to investigate the relative roles of cellular scale length dependent regulators of tension generation on the transduction of work from the cell to whole organ pump function. Specifically, the role of the length dependent Ca2+ sensitivity of tension (Ca50), filament overlap tension dependence, velocity dependence of tension, and tension dependent binding of Ca2+ to Troponin C on metrics of efficient transduction of work and stress and strain homogeneity were predicted by performing simulations in the absence of each of these feedback mechanisms. The length dependent Ca50 and the filament overlap, which make up the Frank-Starling Law, were found to be the two dominant regulators of the efficient transduction of work. Analyzing the fiber velocity field in the absence of the Frank-Starling mechanisms showed that the decreased efficiency in the transduction of work in the absence of filament overlap effects was caused by increased post systolic shortening, whereas the decreased efficiency in the absence of length dependent Ca50 was caused by an inversion in the regional distribution of strain

Validation of diffusion tensor MRI measurements of cardiac microstructure with structure tensor synchrotron radiation imaging.

Background Diffusion tensor imaging (DTI) is widely used to assess tissue microstructure non-invasively. Cardiac DTI enables inference of cell and sheetlet orientations, which are altered under pathological conditions. However, DTI is affected by many factors, therefore robust validation is critical. Existing histological validation is intrinsically flawed, since it requires further tissue processing leading to sample distortion, is routinely limited in field-of-view and requires reconstruction of three-dimensional volumes from two-dimensional images. In contrast, synchrotron radiation imaging (SRI) data enables imaging of the heart in 3D without further preparation following DTI. The objective of the study was to validate DTI measurements based on structure tensor analysis of SRI data. Methods One isolated, fixed rat heart was imaged ex vivo with DTI and X-ray phase contrast SRI, and reconstructed at 100 μm and 3.6 μm isotropic resolution respectively. Structure tensors were determined from the SRI data and registered to the DTI data. Results Excellent agreement in helix angles (HA) and transverse angles (TA) was observed between the DTI and structure tensor synchrotron radiation imaging (STSRI) data, where HADTI-STSRI = −1.4° ± 23.2° and TADTI-STSRI = −1.4° ± 35.0° (mean ± 1.96 standard deviation across all voxels in the left ventricle). STSRI confirmed that the primary eigenvector of the diffusion tensor corresponds with the cardiomyocyte long-axis across the whole myocardium. Conclusions We have used STSRI as a novel and high-resolution gold standard for the validation of DTI, allowing like-with-like comparison of three-dimensional tissue structures in the same intact heart free of distortion. This represents a critical step forward in independently verifying the structural basis and informing the interpretation of cardiac DTI data, thereby supporting the further development and adoption of DTI in structure-based electro-mechanical modelling and routine clinical applications

Quantifying the effect of uncertainty in input parameters in a simplified bidomain model of partial thickness ischaemia

Reduced blood flow in the coronary arteries can lead to damaged heart tissue (myocardial ischaemia). Although one method for detecting myocardial ischaemia involves changes in the ST segment of the electrocardiogram, the relationship between these changes and subendocardial ischaemia is not fully understood. In this study, we modelled ST-segment epicardial potentials in a slab model of cardiac ventricular tissue, with a central ischaemic region, using the bidomain model, which considers conduction longitudinal, transverse and normal to the cardiac fibres. We systematically quantified the effect of uncertainty on the input parameters, fibre rotation angle, ischaemic depth, blood conductivity and six bidomain conductivities, on outputs that characterise the epicardial potential distribution. We found that three typical types of epicardial potential distributions (one minimum over the central ischaemic region, a tripole of minima, and two minima flanking a central maximum) could all occur for a wide range of ischaemic depths. In addition, the positions of the minima were affected by both the fibre rotation angle and the ischaemic depth, but not by changes in the conductivity values. We also showed that the magnitude of ST depression is affected only by changes in the longitudinal and normal conductivities, but not by the transverse conductivities

Parameter estimation in a Holzapfel–Ogden law for healthy myocardium

A central problem in biomechanical studies of personalized human left ventricular (LV) modelling is to estimate material properties from in vivo clinical measurements. In this work we evaluate the passive myocardial mechanical properties inversely from the in vivo LV chamber pressure–volume and strain data. The LV myocardium is described using a structure-based orthotropic Holzapfel–Ogden constitutive law with eight parameters. In the first part of the paper we demonstrate how to use a multi-step non-linear least-squares optimization procedure to inversely estimate the parameters from the pressure–volume and strain data obtained from a synthetic LV model in diastole. In the second part, we show that to apply this procedure to clinical situations with limited in vivo data, additional constraints are required in the optimization procedure. Our study, based on three different healthy volunteers, demonstrates that the parameters of the Holzapfel–Ogden law could be extracted from pressure–volume and strain data with a suitable multi-step optimization procedure. Although the uniqueness of the solution cannot be addressed using our approaches, the material response is shown to be robustly determined