Search

Contractions of the uterus play an important role in menstruation and fertility, and contractile dysfunction can lead to chronic diseases such as endometriosis. However, the structure and function of the uterus are difficult to interrogate in humans, and thus animal studies are often employed to understand its function. In rats, anatomical studies of the uterus have typically been based on histological assessment, have been limited to small segments of the uterine structure, and have been time‐consuming to reconstruct at the organ scale. This study used micro‐computed tomography imaging to visualise the muscle structures in the entire non‐pregnant rat uterus and assess its use for 3D virtual histology. An assessment of the rodent uterus is presented to (i) quantify muscle thickness variations along the horns, (ii) identify predominant fibre orientations of the muscles and (iii) demonstrate how the anatomy of the uterus can be mapped to 3D volumetric meshes via virtual histology. Micro‐computed tomography measurements were validated against measurements from histological sections. The average thickness of the myometrium was found to be 0.33 ± 0.11 mm and 0.31 ± 0.09 mm in the left and right horns, respectively. The micro‐computed tomography and histology thickness calculations were found to correlate strongly at different locations in the uterus: at the cervix, r = 0.87, and along the horn from the cervical end to the ovarian end, respectively, r = 0.77, r = 0.89 and r = 0.54, with p < 0.001 in every location. This study shows that micro‐computed tomography can be used to quantify the musculature in the whole non‐pregnant uterus and can be used for 3D virtual histology. [ABSTRACT FROM AUTHOR]

Hypertensive heart disease (HHD) increases risk of ventricular tachycardia (VT) and ventricular fibrillation (VF). The roles of structural vs. electrophysiological remodelling and age vs. disease progression are not fully understood. This cross‐sectional study of cardiac alterations through HHD investigates mechanistic contributions to VT/VF risk. Risk was electrically assessed in Langendorff‐perfused, spontaneously hypertensive rat hearts at 6, 12 and 18 months, and paced optical membrane voltage maps were acquired from the left ventricular (LV) free wall epicardium. Distributions of LV patchy fibrosis and 3D cellular architecture in representative anterior LV mid‐wall regions were quantified from macroscopic and microscopic fluorescence images of optically cleared tissue. Imaging showed increased fibrosis from 6 months, particularly in the inner LV free wall. Myocyte cross‐section increased at 12 months, while inter‐myocyte connections reduced markedly with fibrosis. Conduction velocity decreased from 12 months, especially transverse to the myofibre direction, with rate‐dependent anisotropy at 12 and 18 months, but not earlier. Action potential duration (APD) increased when clustered by age, as did APD dispersion at 12 and 18 months. Among 10 structural, functional and age variables, the most reliably linked were VT/VF risk, general LV fibrosis, a measure quantifying patchy fibrosis, and non‐age clustered APD dispersion. VT/VF risk related to a quantified measure of patchy fibrosis, but age did not factor strongly. The findings are consistent with the notion that VT/VF risk is associated with rate‐dependent repolarization heterogeneity caused by structural remodelling and reduced lateral electrical coupling between LV myocytes, providing a substrate for heterogeneous intramural activation as HHD progresses. Key points: There is heightened arrhythmic risk with progression of hypertensive heart disease.Risk is related to increasing left ventricular fibrosis, but the nature of this relationship has not been quantified.This study is a novel systematic characterization of changes in active electrical properties and fibrotic remodelling during progression of hypertensive heart disease in a well‐established animal disease model.Arrhythmic risk is predicted by several left ventricular measures, in particular fibrosis quantity and structure, and epicardial action potential duration dispersion.Age alone is not a good predictor of risk.An improved understanding of links between arrhythmic risk and fibrotic architectures in progressive hypertensive heart disease aids better interpretation of late gadolinium‐enhanced cardiac magnetic resonance imaging and electrical mapping signals. [ABSTRACT FROM AUTHOR]

The value of digital twins for prototyping controllers or interventions in a sandbox environment are well-established in engineering and physics. However, this is challenging for biophysics trying to seamlessly compose models of multiple spatial and temporal scale behavior into the digital twin. Two challenges stand out as constraining progress: (i) ensuring physical consistency of conservation laws across composite models and (ii) drawing useful and timely clinical and scientific information from conceptually and computationally complex models. Challenge (i) can be robustly addressed with bondgraphs. However, challenge (ii) is exacerbated using this approach. The complexity question can be looked at from multiple angles. First from the perspective of discretizations that reflect underlying biophysics (functional tissue units) and secondly by exploring maximum entropy as the principle guiding multicellular biophysics. Statistical mechanics, long applied to understanding emergent phenomena from atomic physics, coupled with the observation that cellular architecture in tissue is orchestrated by biophysical constraints on metabolism and communication, shows conceptual promise. This architecture along with cell specific properties can be used to define tissue specific network motifs associated with energetic contributions. Complexity can be addressed based on energy considerations and finding mean measures of dependent variables. A probability distribution of the tissue's network motif can be approximated with exponential random graph models. A prototype problem shows how these approaches could be implemented in practice and the type of information that could be extracted. [ABSTRACT FROM AUTHOR]

Objective: To develop, train and test neural networks for predicting heart surface potentials (HSPs) from body surface potentials (BSPs). The method re-frames traditional inverse problems of electrocardiography into regression problems, constraining the solution space by decomposing signals with multidimensional Gaussian impulse basis functions. Methods: Impulse HSPs were generated with single Gaussian basis functions at discrete heart surface locations and projected to corresponding BSPs using a volume conductor torso model. Both BSP (inputs) and HSP (outputs) were mapped to regular 2D surface meshes and used to train a neural network. Predictive capabilities of the network were tested with unseen synthetic and experimental data. Results: A dense full connected single hidden layer neural network was trained to map body surface impulses to heart surface Gaussian basis functions for reconstructing HSP. Synthetic pulses moving across the heart surface were predicted from the neural network with root mean squared error of $9.1\pm 1.4$ %. Predicted signals were robust to noise up to 20 dB and errors due to displacement and rotation of the heart within the torso were bounded and predictable. A shift of the heart 40 mm toward the spine resulted in a 4% increase in signal feature localization error. The set of training impulse function data could be reduced, and prediction error remained bounded. Recorded HSPs from in-vitro pig hearts were reliably decomposed using space-time Gaussian basis functions. Activation times calculated from predicted HSPs for left-ventricular pacing had a mean absolute error of $10.4\pm 11.4$ ms. Other pacing scenarios were analyzed with similar success. Conclusion: Impulses from Gaussian basis functions are potentially an effective and robust way to train simple neural network data models for reconstructing HSPs from decomposed BSPs. Significance: The HSPs predicted by the neural network can be used to generate activation maps that non-invasively identify features of cardiac electrical dysfunction and can guide subsequent treatment options. [ABSTRACT FROM AUTHOR]

Objective: Cardiovascular Implantable Electronic Devices (CIEDs) are used extensively for treating life-threatening conditions such as bradycardia, atrioventricular block and heart failure. The complicated heterogeneous physical dynamics of patients provide distinct challenges to device development and validation. We address this problem by proposing a device testing framework within the in-silico closed-loop context of patient physiology. Methods: We develop an automated framework to validate CIEDs in closed-loop with a high-level physiologically based computational heart model. The framework includes test generation, execution and evaluation, which automatically guides an integrated stochastic optimization algorithm for exploration of physiological conditions. Conclusion: The results show that using a closed loop device-heart model framework can achieve high system test coverage, while the heart model provides clinically relevant responses. The simulated findings of pacemaker mediated tachycardia risk evaluation agree well with the clinical observations. Furthermore, we illustrate how device programming parameter selection affects the treatment efficacy for specific physiological conditions. Significance: This work demonstrates that incorporating model based closed-loop testing of CIEDs into their design provides important indications of safety and efficacy under constrained physiological conditions. [ABSTRACT FROM AUTHOR]

Objective: Evaluating and testing cardiac electrical devices in a closed-physiologic-loop can help design safety, but this is rarely practical or comprehensive. Furthermore, in silico closed-loop testing with biophysical computer models cannot meet the requirements of time-critical cardiac device systems, while simplified models meeting time-critical requirements may not have the necessary dynamic features. We propose a new high-level (abstracted) physiologically-based computational heart model that is time-critical and dynamic. Methods: The model comprises cardiac regional cellular-electrophysiology types connected by a path model along a conduction network. The regional electrophysiology and paths are modeled with hybrid automata that capture non-linear dynamics, such as action potential and conduction velocity restitution and overdrive suppression. The hierarchy of pacemaker functions is incorporated to generate sinus rhythms, while abnormal automaticity can be introduced to form a variety of arrhythmias such as escape ectopic rhythms. Model parameters are calibrated using experimental data and prior model simulations. Conclusion: Regional electrophysiology and paths in the model match human action potentials, dynamic behavior, and cardiac activation sequences. Connected in closed loop with a pacing device in DDD mode, the model generates complex arrhythmia such as atrioventricular nodal reentry tachycardia. Such device-induced outcomes have been observed clinically and we can establish the key physiological features of the heart model that influence the device operation. Significance: These findings demonstrate how an abstract heart model can be used for device validation and to design personalized treatment. [ABSTRACT FROM AUTHOR]

Altered electrical behavior alongside healed myocardial infarcts (MIs) is associated with increased risk of sudden cardiac death. However, the multidimensional mechanisms are poorly understood and described. This study characterizes, for the first time, the intramural spread of electrical activation in the peri-infarct region of chronic reperfusion MIs. Four sheep were studied 13 wk after antero-apical reperfusion infarction. Extracellular potentials (ECPs) were recorded in a ~20 × 20-mm2 region adjacent to the infarct boundary (25 plunge needles <0.5-mm diameter with 15 electrodes at 1-mm centers) during multisite stimulation. Infarct geometry and electrode locations were reconstructed from magnetic resonance images. Three-dimensional activation spread was characterized by local activation times and interpolated ECP fields (n = 191 records). Control data were acquired in 4 non-infarcted sheep (n = 96 records). Electrodes were distributed uniformly around 15 ± 5% of the intramural infarct boundary. There were marked changes in pacing success and ECP morphology across a functional border zone (BZ) ±2 mm from the boundary. Stimulation adjacent to the infarct boundary was associated with low-amplitude electrical activity within the BZ and delayed activation of surrounding myocardium. Bulk tissue depolarization occurred 3.5–14.6 mm from the pacing site for 39% of stimuli with delays of 4–37 ms, both significantly greater than control (P < 0.0001). Conduction velocity (CV) adjacent to the infarct was not reduced compared with control, consistent with structure-only computer model results. Insignificant CV slowing, irregular stimulus-site specific activation delays, and obvious indirect activation pathways strongly suggest that the substrate for conduction abnormalities in chronic MI is predominantly structural in nature. NEW & NOTEWORTHY Intramural in vivo measurements of peri-infarct electrical activity were not available before this study. We use pace-mapping in a three-dimensional electrode array to show that a subset of stimuli in the peri-infarct region initiates coordinated myocardial activation some distance from the stimulus site with substantial associated time delays. This is site dependent and heterogeneous and occurs for <50% of ectopic stimuli in the border zone. Furthermore, once coordinated activation is initiated, conduction velocity adjacent to the infarct boundary is not significantly different from control. These results give new insights to peri-infarct electrical activity and do not support the widespread view of uniform electrical remodeling in the border zone of chronic myocardial infarcts, with depressed conduction velocity throughout. [ABSTRACT FROM AUTHOR]

Key points: Vagal reflexes slow heart rate and can change where the heartbeat originates within the sinoatrial node (SAN). The mechanisms responsible for this process – termed leading pacemaker (LP) shift – have not been investigated fully.We used optical mapping to measure the effects of baroreflex, chemoreflex and carbachol on pacemaker entrainment and electrical conduction across the SAN.All methods of stimulation triggered shifts in LP site from the central SAN to one or two caudal pacemaker regions. These shifts were associated with reduced current generation capacity centrally and increased electrical load caudally.Previous studies suggest LP shift is a rate‐dependent phenomenon whereby acetylcholine slows central pacemaker rate disproportionately, enabling caudal cells that are less acetylcholine sensitive to assume control. However, our findings indicate the LP region is defined by both pacemaker rate and capacity to drive activation.Shifts in LP site provide an important homeostatic mechanism for rapid switches in heart rate. Reflex vagal activity causes abrupt heart rate slowing with concomitant caudal shifts of the leading pacemaker (LP) site within the sinoatrial node (SAN). However, neither the mechanisms responsible nor their dynamics have been investigated fully. Therefore, the objective of this study was to elucidate the mechanisms driving cholinergic LP shift. Optical maps of right atrial activation were acquired in a rat working heart–brainstem preparation during baroreflex and chemoreflex stimulation or with carbachol. All methods of stimulation triggered shifts in LP site from the central SAN to caudal pacemaker regions, which were positive for HCN4 and received uniform cholinergic innervation. During baroreflex onset, the capacity of the central region to drive activation declined with a decrease in amplitude and gradient of optical action potentials (OAPs) in the surrounding myocardium. Accompanying this decline, there was altered entrainment in the caudal SAN as shown by decreased conduction velocity, OAP amplitude, gradient and activation time. Atropine abolished these responses. Chemoreflex stimulation produced similar effects but central capacity to drive activation was preserved before the LP shift. In contrast, carbachol produced a prolonged period of reduced capacity to drive and altered entrainment. Previous studies suggest LP shift is a rate‐dependent phenomenon whereby acetylcholine slows central pacemaker rate disproportionately, enabling caudal cells that are less acetylcholine sensitive to assume control. Our findings indicate that cholinergic LP shifts are also determined by altered electrical source‐to‐sink balance in the SAN. We conclude that the LP region is defined by both rate and capacity to drive atrial activation. Key points: Vagal reflexes slow heart rate and can change where the heartbeat originates within the sinoatrial node (SAN). The mechanisms responsible for this process – termed leading pacemaker (LP) shift – have not been investigated fully.We used optical mapping to measure the effects of baroreflex, chemoreflex and carbachol on pacemaker entrainment and electrical conduction across the SAN.All methods of stimulation triggered shifts in LP site from the central SAN to one or two caudal pacemaker regions. These shifts were associated with reduced current generation capacity centrally and increased electrical load caudally.Previous studies suggest LP shift is a rate‐dependent phenomenon whereby acetylcholine slows central pacemaker rate disproportionately, enabling caudal cells that are less acetylcholine sensitive to assume control. However, our findings indicate the LP region is defined by both pacemaker rate and capacity to drive activation.Shifts in LP site provide an important homeostatic mechanism for rapid switches in heart rate. [ABSTRACT FROM AUTHOR]

Organ level simulation of bioelectric behavior in the body benefits from flexible and efficient models of cellular membrane potential. These computational organ and cell models can be used to study the impact of pharmaceutical drugs, test hypotheses, assess risk and for closed-loop validation of medical devices. To move closer to the real-time requirements of this modeling a new flexible Fourier based general membrane potential model, called as a Resonant model, is developed that is computationally inexpensive. The new model accurately reproduces non-linear potential morphologies for a variety of cell types. Specifically, the method is used to model human and rabbit sinoatrial node, human ventricular myocyte and squid giant axon electrophysiology. The Resonant models are validated with experimental data and with other published models. Dynamic changes in biological conditions are modeled with changing model coefficients and this approach enables ionic channel alterations to be captured. The Resonant model is used to simulate entrainment between competing sinoatrial node cells. These models can be easily implemented in low-cost digital hardware and an alternative, resource-efficient implementations of sine and cosine functions are presented and it is shown that a Fourier term is produced with two additions and a binary shift. [ABSTRACT FROM AUTHOR]

Background: Gastric slow waves regulate peristalsis, and gastric dysrhythmias have been implicated in functional motility disorders. To accurately define slow wave patterns, it is currently necessary to collect high-resolution serosal recordings during open surgery. We therefore developed a novel gastric slow wave mapping device for use during laparoscopic procedures.Methods: The device consists of a retractable catheter constructed of a flexible nitinol core coated with Pebax. Once deployed through a 5-mm laparoscopic port, the spiral head is revealed with 32 electrodes at 5 mm intervals. Recordings were validated against a reference electrode array in pigs and tested in a human patient.Results: Recordings from the device and a reference array in pigs were identical in frequency (2.6 cycles per minute; p = 0.91), and activation patterns and velocities were consistent (8.9 ± 0.2 vs 8.7 ± 0.1 mm s-1; p = 0.2). Device and reference amplitudes were comparable (1.3 ± 0.1 vs 1.4 ± 0.1 mV; p = 0.4), though the device signal-to-noise ratio was higher (17.5 ± 0.6 vs 12.8 ± 0.6 dB; P < 0.0001). In the human patient, corpus slow waves were recorded and mapped (frequency 2.7 ± 0.03 cycles per minute, amplitude 0.8 ± 0.4 mV, velocity 2.3 ± 0.9 mm s-1).Conclusion: In conclusion, the novel laparoscopic device achieves high-quality serosal slow wave recordings. It can be used for laparoscopic diagnostic studies to document slow wave patterns in patients with gastric motility disorders. [ABSTRACT FROM AUTHOR]

Objective: A flexible, efficient, and verifiable pacemaker cell model is essential to the design of real-time virtual hearts that can be used for closed-loop validation of cardiac devices. A new parametric model of pacemaker action potential is developed to address this need. Methods: The action potential phases are modeled using hybrid automaton with one piecewise-linear continuous variable. The model can capture rate-dependent dynamics, such as action potential duration restitution, conduction velocity restitution, and overdrive suppression by incorporating nonlinear update functions. Simulated dynamics of the model compared well with previous models and clinical data. Conclusion: The results show that the parametric model can reproduce the electrophysiological dynamics of a variety of pacemaker cells, such as sinoatrial node, atrioventricular node, and the His-Purkinje system, under varying cardiac conditions. Significance: This is an important contribution toward closed-loop validation of cardiac devices using real-time heart models. [ABSTRACT FROM PUBLISHER]

Models of the cardiac conduction system are usually at two extremes: (1) high fidelity models with excellent precision but lacking a real-time response for emulation (hardware in the loop simulation); or (2) models amenable for emulation, but that do not exhibit appropriate dynamic response, which is necessary for arrhythmia susceptibility.We introduce two abstractions to remedy the situation. The first abstraction is a new cell model, which is a semi-linear hybrid automata. The proposed model is as computationally efficient as current state-of-the-art cell models amenable for emulation. Yet, unlike these models, it is also able to capture the dynamic response of the cardiac cell like the higher-fidelity models. The second abstraction is the use of smooth-tokens to develop a new path model, connecting cells, which is efficient in terms of memory consumption. Moreover, the memory requirements of the path model can be statically bounded and are invariant to the emulation step size. Results show that the proposed semi-linear abstraction for the cell reduces the execution time by up to 44%. Furthermore, the smooth-tokens based path model reduces the memory consumption by 40 times when compared to existing path models. This paves the way for the emulation of complex cardiac conduction systems, using hardware code-generators. [ABSTRACT FROM AUTHOR]

Hybrid systems combine discrete controllers with adjoining physical processes. While many approaches exist for simulating hybrid systems, there are few approaches for their emulation, especially when the actual physical plant is not available. This paper develops the first formal framework for emulation along with a new compiler that enables large-scale (1000+ components) simulation. We propose a formal model called Synchronous Emulation Automaton (SEA) specifically for modular compilation and parallel execution. SEA combines Linear Time Invariant (LTI) systems with discrete mode switches and has the following semantic differences with Hybrid Automata: 1) the Ordinary Differential Equations are solved analytically and the solutions are sampled at the Worst-Case Reaction Time of the model and 2) we develop a new composition semantics, which allows individual SEAs to execute in parallel with each other. The proposed semantics eliminates: a) the need for dynamic numerical solvers, and b) the Zeno-phenomenon by construction. Experimental results show that process models designed using our tool (Piha) give a 3.6 times execution speedup over Simulink®, and upto 26 times speedup on manycore architectures. [ABSTRACT FROM AUTHOR]

Orthotropic Intramural Spread of LV Electrical Potentials Introduction Extracellular potentials measured on the heart surfaces are used to infer events that originate deep within the heart wall. We have reconstructed intramural potentials in three dimensions for the first time, and compare these with epicardial and endocardial surface potentials and cardiac microstructure. Methods and Results Extracellular potentials from intramural point stimulation were measured from a high density 3-D electrode array in the in vivo pig LV. MR and extended volume imaging were used to register electrode locations and characterize fiber and laminar orientations throughout the recording volume. Measured potentials were compared with predictions of tissue-specific bidomain computer activation models. Positive potentials recorded in the LV wall preceded the depolarization wavefront as it spread in the fiber direction. Transverse to this, passive and active potentials spread preferentially in the laminar direction (anisotropy ratio ∼1.6:1). Epicardial surface potentials reflect initial intramural propagation at the stimulus location, but endocardial potentials do not, particularly adjacent to papillary muscles. Measured 3-D potentials were consistently better captured by computer models that incorporate three distinct conductivities aligned with local microstructural axes, but the preferential spread of potentials in the laminar direction was not fully predicted. Conclusions This study provides evidence for preferential transmural conduction and raises questions about the extent to which intramural electrical events can be inferred from endocardial potentials. [ABSTRACT FROM AUTHOR]

High-resolution (HR) mapping has been used to study gastric slowwave activation; however, the specific characteristics of antral electrophysiology remain poorly defined. This study applied HR mapping and computational modeling to define functional human antral physiology. HR mapping was performed in 10 subjects using flexible electrode arrays (128-192 electrodes; 16-24 cm2) arranged from the pylorus to mid-corpus. Anatomical registration was by photographs and anatomical landmarks. Slow-wave parameters were computed, and resultant data were incorporated into a computational fluid dynamics (CFD) model of gastric flow to calculate impact on gastric mixing. In all subjects, extracellular mapping demonstrated normal aboral slow-wave propagation and a region of increased amplitude and velocity in the prepyloric antrum. On average, the high-velocity region commenced 28 mm proximal to the pylorus, and activation ceased 6 mm from the pylorus. Within this region, velocity increased 0.2 mm/s per mm of tissue, from the mean 3.3 - 0.1 mm/s to 7.5 - 0.6 mm/s (P < 0.001), and extracellular amplitude increased from 1.5 ± 0.1 mV to 2.5 ± 0.1 mV (P < 0.001). CFD modeling using representative parameters quantified a marked increase in antral recirculation, resulting in an enhanced gastric mixing, due to the accelerating terminal antral contraction. The extent of gastric mixing increased almost linearly with the maximal velocity of the contraction. In conclusion, the human terminal antral contraction is controlled by a short region of rapid high-amplitude slow-wave activity. Distal antral wave acceleration plays a major role in antral flow and mixing, increasing particle strain and trituration. [ABSTRACT FROM AUTHOR]

Goal: Gastric motility disorders have been associated with abnormal slow wave electrical activity (gastric dysrhythmias). Gastric pacing is a potential therapy for gastric dysrhythmias; however, new pacing protocols are required that can effectively modulate motility patterns, while being power efficient. This study presents a novel comprehensive 3-D multiscale modeling framework of the human stomach, including anisotropic conduction, capable of evaluating pacing strategies. Methods: A high-resolution anatomically realistic mesh was generated from CT images taken from a human stomach. Principal conduction axes were calculated and embedded within this model based on a modified Laplace–Dirichlet rule-based algorithm. A continuum-based tridomain formulation was implemented and evaluated for performance and used to model the slow-wave propagation, which takes into account the two main cell types present in gastric musculature. Model parameters were found by matching predicted normal slow-wave activity to experimental observation and data. These simulation parameters were applied while modeling an external pacing event to entrain slow-wave patterns. Results: The proposed formulation was found to be two times more efficient than a previous formulation for a normal slow-wave simulation. Convergence analysis showed that a mesh resolution of \approx0.4\-0.5 \text mm is required for an accurate solution process. Conclusion: The effect of different pacing frequencies on entrainment demonstrated that the pacing protocols are limited by the frequency of the native propagation and the refractory period of the cellular activity. Significance: The model is expected to become an important tool in studying pacing protocols for both efficiency and effectiveness. [ABSTRACT FROM AUTHOR]

Ectopic electrical activity that originates in the peri-infarct region can give rise to potentially lethal re-entrant arrhythmias. The spatial variation in electrotonic loading that results from structural remodelling in the infarct border zone may increase the probability that focal activity will trigger electrical capture, but this has not previously been investigated systematically. This study uses in-silico experiments to examine the structural modulation of effective refractory period on ectopic beat capture. Informed by 3-D reconstructions of myocyte organization in the infarct border zone, a region of rapid tissue expansion is abstracted to an idealized representation. A novel metric is introduced that defines the local electrotonic loading as a function of passive tissue properties and boundary conditions. The effective refractory period correlates closely with local electrotonic loading, while the action potential duration, conduction, and upstroke velocity reduce in regions of increasing electrotonic load. In the presence of focal ectopic stimuli, spatial variation in effective refractory period can cause unidirectional conduction block providing a substrate for reentrant arrhythmias. Consequently, based on the observed results, a possible novel mechanism for arrhythmogenesis in the infarct border zone is proposed. [ABSTRACT FROM AUTHOR]

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 (v1 ST), intermediate (v2 ST) and least (v3 ST) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e1 DTI), intermediate (e2 DTI) and least (e3 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 (v1 ST) agreed well with that of diffusion (e1 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 (v3 ST) and diffusion (e3 ST) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v3 ST) and DTI (e3 DTI) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v3 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 v3 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. [ABSTRACT FROM AUTHOR]

A theoretical model of the cross-linking topology of venoicular muscle tissue is developed. Using parameter estimation the terms of the theoretical model are estimated for normal and pathological conditions. The model represents the anisotropic structure of the tissue, reproduces published experimental data and characterizes the role of different tissue components in the observed macroscopic behavior. Changes in the material parameters are consistent with expected structural changes and the model is extended to reproduce force-Calcium relationships. Model results are invoked to argue that semisoft behavior and the material axis anisotropy arise from the constraints on the extracellular matrix cross-linking topology. [ABSTRACT FROM AUTHOR]