Raimond L. Winslow

Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging

Diffusion tensor magnetic resonance imaging (MRI) is a possible new means of elucidating the anatomic structure of the myocardium. It enjoys several advantages over traditional histological approaches, including the ability to rapidly measure fiber organization in isolated, perfused, arrested hearts, thereby avoiding fixation and sectioning of artifacts. However, quantitative validation of this MRI method has been lacking. Here, fiber orientations estimated in the same locations in the same heart using both diffusion tensor MRI and histology are compared in a total of two perfused rabbit hearts. Fiber orientations were statistically similar for both methods and differed on average by 12 degrees at any single location. This is similar to the 10 degrees uncertainty in fiber orientation achieved with histology. In addition, imaging studies performed in a total of seven hearts support a level of organization beyond the myofiber, the recently described laminar organization of the ventricular myocardium.Diffusion tensor magnetic resonance imaging (MRI) is a possible new means of elucidating the anatomic structure of the myocardium. It enjoys several advantages over traditional histological approaches, including the ability to rapidly measure fiber organization in isolated, perfused, arrested hearts, thereby avoiding fixation and sectioning of artifacts. However, quantitative validation of this MRI method has been lacking. Here, fiber orientations estimated in the same locations in the same heart using both diffusion tensor MRI and histology are compared in a total of two perfused rabbit hearts. Fiber orientations were statistically similar for both methods and differed on average by 12° at any single location. This is similar to the 10° uncertainty in fiber orientation achieved with histology. In addition, imaging studies performed in a total of seven hearts support a level of organization beyond the myofiber, the recently described laminar organization of the ventricular myocardium.

An Integrated Model of Cardiac Mitochondrial Energy Metabolism and Calcium Dynamics

We present an integrated thermokinetic model describing control of cardiac mitochondrial bioenergetics. The model describes the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and mitochondrial Ca(2+) handling. The kinetic component of the model includes effectors of the TCA cycle enzymes regulating production of NADH and FADH(2), which in turn are used by the electron transport chain to establish a proton motive force (Delta mu(H)), driving the F(1)F(0)-ATPase. In addition, mitochondrial matrix Ca(2+), determined by Ca(2+) uniporter and Na(+)/Ca(2+) exchanger activities, regulates activity of the TCA cycle enzymes isocitrate dehydrogenase and alpha-ketoglutarate dehydrogenase. The model is described by twelve ordinary differential equations for the time rate of change of mitochondrial membrane potential (Delta Psi(m)), and matrix concentrations of Ca(2+), NADH, ADP, and TCA cycle intermediates. The model is used to predict the response of mitochondria to changes in substrate delivery, metabolic inhibition, the rate of adenine nucleotide exchange, and Ca(2+). The model is able to reproduce, qualitatively and semiquantitatively, experimental data concerning mitochondrial bioenergetics, Ca(2+) dynamics, and respiratory control. Significant increases in oxygen consumption (V(O(2))), proton efflux, NADH, and ATP synthesis, in response to an increase in cytoplasmic Ca(2+), are obtained when the Ca(2+)-sensitive dehydrogenases are the main rate-controlling steps of respiratory flux. These responses diminished when control is shifted downstream (e.g., the respiratory chain or adenine nucleotide translocator). The time-dependent behavior of the model, under conditions simulating an increase in workload, closely reproduces experimentally observed mitochondrial NADH dynamics in heart trabeculae subjected to changes in pacing frequency. The steady-state and time-dependent behavior of the model support the hypothesis that mitochondrial matrix Ca(2+) plays an important role in matching energy supply with demand in cardiac myocytes.

Cardiac Ca2+ dynamics: the roles of ryanodine receptor adaptation and sarcoplasmic reticulum load.

We construct a detailed mathematical model for Ca2+ regulation in the ventricular myocyte that includes novel descriptions of subcellular mechanisms based on recent experimental findings: 1) the Keizer-Levine model for the ryanodine receptor (RyR), which displays adaptation at elevated Ca2+; 2) a model for the L-type Ca2+ channel that inactivates by mode switching; and 3) a restricted subspace into which the RyRs and L-type Ca2+ channels empty and interact via Ca2+. We add membrane currents from the Luo-Rudy Phase II ventricular cell model to our description of Ca2+ handling to formulate a new model for ventricular action potentials and Ca2+ regulation. The model can simulate Ca2+ transients during an action potential similar to those seen experimentally. The subspace [Ca2+] rises more rapidly and reaches a higher level (10-30 microM) than the bulk myoplasmic Ca2+ (peak [Ca2+]i approximately 1 microM). Termination of sarcoplasmic reticulum (SR) Ca2+ release is predominately due to emptying of the SR, but is influenced by RyR adaptation. Because force generation is roughly proportional to peak myoplasmic Ca2+, we use [Ca2+]i in the model to explore the effects of pacing rate on force generation. The model reproduces transitions seen in force generation due to changes in pacing that cannot be simulated by previous models. Simulation of such complex phenomena requires an interplay of both RyR adaptation and the degree of SR Ca2+ loading. This model, therefore, shows improved behavior over existing models that lack detailed descriptions of subcellular Ca2+ regulatory mechanisms.

Classifying Gene Expression Profiles from Pairwise mRNA Comparisons

We present a new approach to molecular classification based on mRNA comparisons. Our method, referred to as the top-scoring pair(s) (TSP) classifier, is motivated by current technical and practical limitations in using gene expression microarray data for class prediction, for example to detect disease, identify tumors or predict treatment response. Accurate statistical inference from such data is difficult due to the small number of observations, typically tens, relative to the large number of genes, typically thousands. Moreover, conventional methods from machine learning lead to decisions which are usually very difficult to interpret in simple or biologically meaningful terms. In contrast, the TSP classifier provides decision rules which i) involve very few genes and only relative expression values (e.g., comparing the mRNA counts within a single pair of genes); ii) are both accurate and transparent; and iii) provide specific hypotheses for follow-up studies. In particular, the TSP classifier achieves prediction rates with standard cancer data that are as high as those of previous studies which use considerably more genes and complex procedures. Finally, the TSP classifier is parameter-free, thus avoiding the type of over-fitting and inflated estimates of performance that result when all aspects of learning a predictor are not properly cross-validated.

Measuring and mapping cardiac fiber and laminar architecture using diffusion tensor MR imaging

Abstract: The ventricular myocardium is known to exhibit a complex spatial organization, with fiber orientation varying as a function of transmural location. It is now well established that diffusion tensor magnetic resonance imaging (DTMRI) may be used to measure this fiber orientation at high spatial resolution. Cardiac fibers are also known to be organized in sheets with surface orientation varying throughout the ventricles. This article reviews results on use of DTMRI for measuring ventricular fiber orientation, as well as presents new results providing strong evidence that the tertiary eigenvector of the diffusion tensor is aligned locally with the cardiac sheet surface normal. Considered together, these data indicate that DTMRI may be used to reconstruct both ventricular fiber and sheet organization. This article also presents the large deformation diffeomorphic metric mapping (LDDMM) algorithm and shows that this algorithm may be used to bring ensembles of imaged and reconstructed hearts into correspondence (e.g., registration) so that variability of ventricular geometry, fiber, and sheet orientation may be quantified. Ventricular geometry and fiber structure is known to be remodeled in a range of disease processes; however, descriptions of this remodeling have remained subjective and qualitative. We anticipate that use of DTMRI for reconstruction of ventricular anatomy coupled with application of the LDDMM method for image volume registration will enable the detection and quantification of changes in cardiac anatomy that are characteristic of specific disease processes in the heart. Finally, we show that epicardial electrical mapping and DTMRI imaging may be performed in the same hearts. The anatomic data may then be used to simulate electrical conduction in a computational model of the very same heart that was mapped electrically. This facilitates direct comparison and testing of model versus experimental results and opens the door to quantitative measurement, modeling, and analysis of the ways in which remodeling of ventricular microanatomy influences electrical conduction in the heart.

Role of the calcium-independent transient outward current I(to1) in shaping action potential morphology and duration

The Kv4.3-encoded current (IKv4.3) has been identified as the major component of the voltage-dependent Ca2+-independent transient outward current (Ito1) in human and canine ventricular cells. Experimental evidence supports a correlation between Ito1 density and prominence of the phase 1 notch; however, the role of Ito1 in modulating action potential duration (APD) remains unclear. To help resolve this role, Markov state models of the human and canine Kv4.3- and Kv1.4-encoded currents at 35°C are developed on the basis of experimental measurements. A model of canine Ito1 is formulated as the combination of these Kv4.3 and Kv1.4 currents and is incorporated into an existing canine ventricular myocyte model. Simulations demonstrate strong coupling between L-type Ca2+ current and IKv4.3 and predict a bimodal relationship between IKv4.3 density and APD whereby perturbations in IKv4.3 density may produce either prolongation or shortening of APD, depending on baseline Ito1 current level.

An integrative model of the cardiac ventricular myocyte incorporating local control of Ca2+ release.

The local control theory of excitation-contraction (EC) coupling in cardiac muscle asserts that L-type Ca(2+) current tightly controls Ca(2+) release from the sarcoplasmic reticulum (SR) via local interaction of closely apposed L-type Ca(2+) channels (LCCs) and ryanodine receptors (RyRs). These local interactions give rise to smoothly graded Ca(2+)-induced Ca(2+) release (CICR), which exhibits high gain. In this study we present a biophysically detailed model of the normal canine ventricular myocyte that conforms to local control theory. The model formulation incorporates details of microscopic EC coupling properties in the form of Ca(2+) release units (CaRUs) in which individual sarcolemmal LCCs interact in a stochastic manner with nearby RyRs in localized regions where junctional SR membrane and transverse-tubular membrane are in close proximity. The CaRUs are embedded within and interact with the global systems of the myocyte describing ionic and membrane pump/exchanger currents, SR Ca(2+) uptake, and time-varying cytosolic ion concentrations to form a model of the cardiac action potential (AP). The model can reproduce both the detailed properties of EC coupling, such as variable gain and graded SR Ca(2+) release, and whole-cell phenomena, such as modulation of AP duration by SR Ca(2+) release. Simulations indicate that the local control paradigm predicts stable APs when the L-type Ca(2+) current is adjusted in accord with the balance between voltage- and Ca(2+)-dependent inactivation processes as measured experimentally, a scenario where common pool models become unstable. The local control myocyte model provides a means for studying the interrelationship between microscopic and macroscopic behaviors in a manner that would not be possible in experiments.

Ex vivo 3D diffusion tensor imaging and quantification of cardiac laminar structure

A three‐dimensional (3D) diffusion‐weighted imaging (DWI) method for measuring cardiac fiber structure at high spatial resolution is presented. The method was applied to the ex vivo reconstruction of the fiber architecture of seven canine hearts. A novel hypothesis‐testing method was developed and used to show that distinct populations of secondary and tertiary eigenvalues may be distinguished at reasonable confidence levels (P ≤ 0.01) within the canine ventricle. Fiber inclination and sheet angles are reported as a function of transmural depth through the anterior, lateral, and posterior left ventricle (LV) free wall. Within anisotropic regions, two consistent and dominant orientations were identified, supporting published results from histological studies and providing strong evidence that the tertiary eigenvector of the diffusion tensor (DT) defines the sheet normal. Magn Reson Med, 2005. Published 2005 Wiley‐Liss, Inc.

Single-tone intensity discrimination based on auditory-nerve rate responses in backgrounds of quiet, noise, and with stimulation of the crossed olivocochlear bundle

We use simple statistical models of the firing patterns of high, medium, and low spontaneous rate auditory-nerve fibers to study mechanisms which determine the overall dynamic range of the auditory periphery. The models relate experimentally measured rate response properties of fibers with best frequency (BF) near 8.0 kHz to their ability to encode changes in BF tone level by changes in discharge rate in backgrounds of quiet and noise, with and without electrical stimulation of the crossed olivocochlear bundle (COCB). Application of the models to the BF tone rate responses of auditory-nerve fibers in backgrounds of quiet shows that optimum processing of the rate responses of fibers with BF near 8.0 kHz yields performance in the intensity discrimination task meeting or exceeding that of human subjects over an 80 dB range of levels. By defining a statistical measure of dynamic range, we confirm the results of Costalupes et al. (1984) demonstrating that masking noise shifts the dynamic range of auditory-nerve fibers to higher stimulus levels, thus preventing rate saturation. However, model analysis shows that masking noise also produces large reductions of dynamic range as well as large increases in the minimum intensity difference that can be encoded by the rate responses of single and ensembles of fibers. Electrical stimulation of the COCB can restore auditory-nerve fiber dynamic range and sensitivity to changes in BF tone level in noise backgrounds, in some cases to roughly that observed in backgrounds of quiet.

Direct histological validation of diffusion tensor MRI in formaldehyde- fixed myocardium

Diffusion tensor MRI is emerging as a rapid, nondestructive method to map myocardial fiber organization. It accurately measures myofiber orientation in hearts bathed in or perfused with cardioplegic solution. This study shows it also accurately maps the fibrous architecture of formalin‐fixed hearts. Fiber orientations obtained by MRI and histology at the same locations in an excised portion of rabbit ventricle differed on average by 3.7° (SD = 6.4°, N = 70), a closer correspondence than achieved with previous preparations. The longer acquisition times afforded by fixed‐heart imaging provides better accuracy, and should enable high‐resolution reconstruction of the entire ventricular architecture. Magn Reson Med 44:157–161, 2000.