In silico assessment of biomedical products: the conundrum of rare but not so rare events in two case studies

In silico clinical trials, defined as “The use of individualized computer simulation in the development or regulatory evaluation of a medicinal product, medical device, or medical intervention,” have been proposed as a possible strategy to reduce the regulatory costs of innovation and the time to market for biomedical products. We review some of the the literature on this topic, focusing in particular on those applications where the current practice is recognized as inadequate, as for example, the detection of unexpected severe adverse events too rare to be detected in a clinical trial, but still likely enough to be of concern. We then describe with more details two case studies, two successful applications of in silico clinical trial approaches, one relative to the University of Virginia/Padova simulator that the Food and Drug Administration has accepted as possible replacement for animal testing in the preclinical assessment of artificial pancreas technologies, and the second, an investigation of the probability of cardiac lead fracture, where a Bayesian network was used to combine in vivo and in silico observations, suggesting a whole new strategy of in silico-augmented clinical trials, to be used to increase the numerosity where recruitment is impossible, or to explore patients’ phenotypes that are unlikely to appear in the trial cohort, but are still frequent enough to be of concern

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

Exploring the Fundamental Dynamics of Error-Based Motor Learning Using a Stationary Predictive-Saccade Task

The maintenance of movement accuracy uses prior performance errors to correct future motor plans; this motor-learning process ensures that movements remain quick and accurate. The control of predictive saccades, in which anticipatory movements are made to future targets before visual stimulus information becomes available, serves as an ideal paradigm to analyze how the motor system utilizes prior errors to drive movements to a desired goal. Predictive saccades constitute a stationary process (the mean and to a rough approximation the variability of the data do not vary over time, unlike a typical motor adaptation paradigm). This enables us to study inter-trial correlations, both on a trial-by-trial basis and across long blocks of trials. Saccade errors are found to be corrected on a trial-by-trial basis in a direction-specific manner (the next saccade made in the same direction will reflect a correction for errors made on the current saccade). Additionally, there is evidence for a second, modulating process that exhibits long memory. That is, performance information, as measured via inter-trial correlations, is strongly retained across a large number of saccades (about 100 trials). Together, this evidence indicates that the dynamics of motor learning exhibit complexities that must be carefully considered, as they cannot be fully described with current state-space (ARMA) modeling efforts

Regulation of Coronary Blood Flow

The heart is uniquely responsible for providing its own blood supply through the coronary circulation. Regulation of coronary blood flow is quite complex and, after over 100 years of dedicated research, is understood to be dictated through multiple mechanisms that include extravascular compressive forces (tissue pressure), coronary perfusion pressure, myogenic, local metabolic, endothelial as well as neural and hormonal influences. While each of these determinants can have profound influence over myocardial perfusion, largely through effects on end-effector ion channels, these mechanisms collectively modulate coronary vascular resistance and act to ensure that the myocardial requirements for oxygen and substrates are adequately provided by the coronary circulation. The purpose of this series of Comprehensive Physiology is to highlight current knowledge regarding the physiologic regulation of coronary blood flow, with emphasis on functional anatomy and the interplay between the physical and biological determinants of myocardial oxygen delivery. © 2017 American Physiological Society. Compr Physiol 7:321-382, 2017

Computational evidence for protein-mediated fatty acid transport across the sarcolemma

Long-chain fatty acids (FAs) are important substrates used by the heart to fulfil its energy requirements. Prior to mitochondrial oxidation, blood-borne FAs must pass through the cell membrane of the cardiac myocyte (sarcolemma). The mechanism underlying the sarcolemmal transport of FAs is incompletely understood. The aim of the present study was to estimate the trans-sarcolemmal FA uptake rate using a comprehensive computer model, in which the most relevant mechanisms proposed for cardiac FA uptake were incorporated. Our in silico findings show that diffusion of FA, present in its unbound form (uFA) in close proximity to the outer leaflet of the sarcolemma and serving as sole FA source, is insufficient to account for the physiological FA uptake rate. The inclusion of a hypothetical membrane-associated FA-TFPC (FA-transport-facilitating protein complex) in the model calculations substantially increased the FA uptake rate across the sarcolemma. The model requires that the biological properties of the FA-TFPC allow for increasing the rate of absorption of FA into the outer leaflet and the ‘flip-flop’ rate of FA from the outer to the inner leaflet of the sarcolemma. Experimental studies have identified various sarcolemma-associated proteins promoting cardiac FA uptake. It remains to be established whether these proteins possess the properties predicted by our model. Our findings also indicate that albumin receptors located on the outer leaflet of the sarcolemma facilitate the transfer of FA across the membrane to a significant extent. The outcomes of the computer simulations were verified with physiologically relevant FA uptake rates as assessed in the intact, beating heart in experimental studies

The shunt problem: control of functional shunting in normal and tumour vasculature

Networks of blood vessels in normal and tumour tissues have heterogeneous structures, with widely varying blood flow pathway lengths. To achieve efficient blood flow distribution, mechanisms for the structural adaptation of vessel diameters must be able to inhibit the formation of functional shunts (whereby short pathways become enlarged and flow bypasses long pathways). Such adaptation requires information about tissue metabolic status to be communicated upstream to feeding vessels, through conducted responses. We propose that impaired vascular communication in tumour microvascular networks, leading to functional shunting, is a primary cause of dysfunctional microcirculation and local hypoxia in cancer. We suggest that anti-angiogenic treatment of tumours may restore vascular communication and thereby improve or normalize flow distribution in tumour vasculature