Increasing the scan-efficiency of pulmonary imaging at 0.55 T using iterative concomitant field and motion-corrected reconstruction

PURPOSE

To develop an iterative concomitant field and motion corrected (iCoMoCo) reconstruction for isotropic high-resolution UTE pulmonary imaging at 0.55 T.

METHODS

A free-breathing golden-angle stack-of-spirals UTE sequence was used to acquire data for 8 min with prototype and commercial 0.55 T MRI scanners. The data was binned into 12 respiratory phases based on superior-inferior navigator readouts. The previously published iterative motion corrected (iMoCo) reconstruction was extended to include concomitant field correction directly in the cost function. The reconstruction was implemented within the Gadgetron framework for inline reconstruction. Data were retrospectively reconstructed to simulate scan times of 2, 4, 6, and 8 min. Image quality was assessed using apparent SNR and image sharpness. The technique was evaluated in healthy volunteers and patients with known lung pathology including coronavirus disease 2019 infection, chronic granulomatous disease, lymphangioleiomyomatosis, and lung nodules.

RESULTS

The technique provided diagnostic-quality images, and image quality was maintained with a slight loss in SNR for simulated scan times down to 4 min. Parenchymal apparent SNR was 4.33 ± 0.57, 5.96 ± 0.65, 7.36 ± 0.64, and 7.87 ± 0.65 using iCoMoCo with scan times of 2, 4, 6, and 8 min, respectively. Image sharpness at the diaphragm was comparable between iCoMoCo and reference images. Concomitant field corrections visibly improved the sharpness of anatomical structures away from the isocenter. Inline image reconstruction and artifact correction were achieved in <5 min.

CONCLUSION

The proposed iCoMoCo pulmonary imaging technique can generate diagnostic quality images with 1.75 mm isotropic resolution in less than 5 min using a 6-min acquisition, on a 0.55 T scanner.

Fully automated contrast selection of joint bright- and black-blood late gadolinium enhancement imaging for robust myocardial scar assessment

PURPOSE

Joint bright- and black-blood MRI techniques provide improved scar localization and contrast. Black-blood contrast is obtained after the visual selection of an optimal inversion time (TI) which often results in uncertainties, inter- and intra-observer variability and increased workload. In this work, we propose an artificial intelligence-based algorithm to enable fully automated TI selection and simplify myocardial scar imaging.

METHODS

The proposed algorithm first localizes the left ventricle using a U-Net architecture. The localized left cavity centroid is extracted and a squared region of interest ("focus box") is created around the resulting pixel. The focus box is then propagated on each image and the sum of the pixel intensity inside is computed. The smallest sum corresponds to the image with the lowest intensity signal within the blood pool and healthy myocardium, which will provide an ideal scar-to-blood contrast. The image's corresponding TI is considered optimal. The U-Net was trained to segment the epicardium in 177 patients with binary cross-entropy loss. The algorithm was validated retrospectively in 152 patients, and the agreement between the algorithm and two magnetic resonance (MR) operators' prediction of TI values was calculated using the Fleiss' kappa coefficient. Thirty focus box sizes, ranging from 2.3mm2 to 20.3cm2, were tested. Processing times were measured.

RESULTS

The U-Net's Dice score was 93.0 ± 0.1%. The proposed algorithm extracted TI values in 2.7 ± 0.1 s per patient (vs. 16.0 ± 8.5 s for the operator). An agreement between the algorithm's prediction and the MR operators' prediction was found in 137/152 patients (κ= 0.89), for an optimal focus box of size 2.3cm2.

CONCLUSION

The proposed fully-automated algorithm has potential of reducing uncertainties, variability, and workload inherent to manual approaches with promise for future clinical implementation for joint bright- and black-blood MRI.

Fully automated planning for anatomical fetal brain MRI on 0.55T

PURPOSE

Widening the availability of fetal MRI with fully automatic real-time planning of radiological brain planes on 0.55T MRI.

METHODS

Deep learning-based detection of key brain landmarks on a whole-uterus echo planar imaging scan enables the subsequent fully automatic planning of the radiological single-shot Turbo Spin Echo acquisitions. The landmark detection pipeline was trained on over 120 datasets from varying field strength, echo times, and resolutions and quantitatively evaluated. The entire automatic planning solution was tested prospectively in nine fetal subjects between 20 and 37 weeks. A comprehensive evaluation of all steps, the distance between manual and automatic landmarks, the planning quality, and the resulting image quality was conducted.

RESULTS

Prospective automatic planning was performed in real-time without latency in all subjects. The landmark detection accuracy was 4.2± $$ \pm $$ 2.6 mm for the fetal eyes and 6.5± $$ \pm $$ 3.2 for the cerebellum, planning quality was 2.4/3 (compared to 2.6/3 for manual planning) and diagnostic image quality was 2.2 compared to 2.1 for manual planning.

CONCLUSIONS

Real-time automatic planning of all three key fetal brain planes was successfully achieved and will pave the way toward simplifying the acquisition of fetal MRI thereby widening the availability of this modality in nonspecialist centers.

Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region

BACKGROUND

While the microstructure of the left ventricle (LV) has been largely described, only a few studies investigated the right ventricular insertion point (RVIP). It was accepted that the aggregate cardiomyocytes organization was much more complex due to the intersection of the ventricular cavities but a precise structural characterization in the human heart was lacking even if clinical phenotypes related to right ventricular wall stress or arrhythmia were observed in this region.

METHODS

MRI-derived anatomical imaging (150 µm3) and diffusion tensor imaging (600 µm3) were performed in large mammalian whole hearts (human: N = 5, sheep: N = 5). Fractional anisotropy, aggregate cardiomyocytes orientations and tractography were compared within both species. Aggregate cardiomyocytes orientation on one ex-vivo sheep whole heart was then computed using structure tensor imaging (STI) from 21 µm isotropic acquisition acquired with micro computed tomography (MicroCT) imaging. Macroscopic and histological examination were performed. Lastly, experimental cardiomyocytes orientation distribution was then compared to the usual rule-based model using electrophysiological (EP) modeling. Electrical activity was modeled with the monodomain formulation.

RESULTS

The RVIP at the level of the inferior ventricular septum presented a unique arrangement of aggregate cardiomyocytes. An abrupt, mid-myocardial change in cardiomyocytes orientation was observed, delimiting a triangle-shaped region, present in both sheep and human hearts. FA's histogram distribution (mean ± std: 0.29 ± 0.06) of the identified region as well as the main dimension (22.2 mm ± 5.6 mm) was found homogeneous across samples and species. Averaged volume is 0.34 cm3 ± 0.15 cm3. Both local activation time (LAT) and morphology of pseudo-ECGs were strongly impacted with delayed LAT and change in peak-to-peak amplitude in the simulated wedge model.

CONCLUSION

The study was the first to describe the 3D cardiomyocytes architecture of the basal inferoseptal left ventricle region in human hearts and identify the presence of a well-organized aggregate cardiomyocytes arrangement and cardiac structural discontinuities. The results might offer a better appreciation of clinical phenotypes like RVIP-late gadolinium enhancement or uncommon idiopathic ventricular arrhythmias (VA) originating from this region.

Remotely detuned receiver coil for high-resolution interventional cardiac magnetic resonance imaging

Introduction

Interventional cardiac MRI in the context of the treatment of cardiac arrhythmia requires submillimeter image resolution to precisely characterize the cardiac substrate and guide the catheter-based ablation procedure in real-time. Conventional MRI receiver coils positioned on the thorax provide insufficient signal-to-noise ratio (SNR) and spatial selectivity to satisfy these constraints.

Methods

A small circular MRI receiver coil was developed and evaluated under different experimental conditions, including high-resolution MRI anatomical and thermometric imaging at 1.5 T. From the perspective of developing a therapeutic MR-compatible catheter equipped with a receiver coil, we also propose alternative remote active detuning techniques of the receiver coil using one or two cables. Theoretical details are presented, as well as simulations and experimental validation.

Results

Anatomical images of the left ventricle at 170 µm in-plane resolution are provided on ex vivo beating heart from swine using a 2 cm circular receiver coil. Taking advantage of the increase of SNR at its vicinity (up to 35 fold compared to conventional receiver coils), real-time MR-temperature imaging can reach an uncertainty below 0.1°C at the submillimetric spatial resolution. Remote active detuning using two cables has similar decoupling efficiency to conventional on-site decoupling, at the cost of an acceptable decrease in the resulting SNR.

Discussion

This study shows the potential of small dimension surface coils for minimally invasive therapy of cardiac arrhythmia intraoperatively guided by MRI. The proposed remote decoupling approaches may simplify the construction process and reduce the cost of such single-use devices.

Neural fatigue by passive induction: repeated stimulus exposure results in cognitive fatigue and altered representations in task-relevant networks

Cognitive fatigue is defined by a reduced capacity to perform mental tasks. Despite its pervasiveness, the underlying neural mechanisms remain elusive. Specifically, it is unclear whether prolonged effort affects performance through alterations in over-worked task-relevant neuronal assemblies. Our paradigm based on repeated passive visual stimulation discerns fatigue effects from the influence of motivation, skill and boredom. We induced performance loss and observed parallel alterations in the neural blueprint of the task, by mirroring behavioral performance with multivariate neuroimaging techniques (MVPA) that afford a subject-specific approach. Crucially, functional areas that responded the most to repeated stimulation were also the most affected. Finally, univariate analysis revealed clusters displaying significant disruption within the extrastriate visual cortex. In sum, here we show that repeated stimulation impacts the implicated brain areas' activity and causes tangible behavioral repercussions, providing evidence that cognitive fatigue can result from local, functional, disruptions in the neural signal induced by protracted recruitment.

3D motion strategy for online volumetric thermometry using simultaneous multi-slice EPI at 1.5T: an evaluation study

PURPOSE

In presence of respiratory motion, temperature mapping is altered by in-plane and through-plane displacements between successive acquisitions together with periodic phase variations. Fast 2D Echo Planar Imaging (EPI) sequence can accommodate intra-scan motion, but limited volume coverage and inter-scan motion remain a challenge during free-breathing acquisition since position offsets can arise between the different slices.

METHOD

To address this limitation, we evaluated a 2D simultaneous multi-slice EPI sequence with multiband (MB) acceleration during radiofrequency ablation on a mobile gel and in the liver of a volunteer (no heating). The sequence was evaluated in terms of resulting inter-scan motion, temperature uncertainty and elevation, potential false-positive heating and repeatability. Lastly, to account for potential through-plane motion, a 3D motion compensation pipeline was implemented and evaluated.

RESULTS

In-plane motion was compensated whatever the MB factor and temperature distribution was found in agreement during both the heating and cooling periods. No obvious false-positive temperature was observed under the conditions being investigated. Repeatability of measurements results in a 95% uncertainty below 2 °C for MB1 and MB2. Uncertainty up to 4.5 °C was reported with MB3 together with the presence of aliasing artifacts. Lastly, fast simultaneous multi-slice EPI combined with 3D motion compensation reduce residual out-of-plane motion.

CONCLUSION

Volumetric temperature imaging (12 slices/700 ms) could be performed with 2 °C accuracy or less, and offer tradeoffs in acquisition time or volume coverage. Such a strategy is expected to increase procedure safety by monitoring large volumes more rapidly for MR-guided thermotherapy on mobile organs.

The interplay of sex-specific electrophysiology and repolarization heterogeneity governs ventricular arrhythmia sustainability in women

Background

The role of sex-specific electrophysiology in the development of lethal ventricular arrhythmia is poorly understood. Since female ventricles have a longer action potential duration (APD) than male ventricles, we hypothesized that female hearts are more vulnerable to sustained ventricular arrhythmia in the presence of known arrhythmogenic substrates.

Purpose

To determine if prolonged APD in females leads to sustained ventricular arrhythmia in the presence of left ventricular (LV) hypertrophy, diffuse fibrosis, and repolarization heterogeneity.

Methods

Magnetic resonance imaging (MRI) data were used to develop computational biventricular models of two aortic stenosis patients (male, 54 years; female, 64 years) before (LV hypertrophy) and 3 months after (regressed LV hypertrophy) aortic valve replacement (AVR). These patients were chosen since they exhibit reduced LV mass, less fibrosis, and no ventricular arrhythmias after AVR. The models were assigned sex-specific cellular electrophysiology and an apicobasal repolarization gradient derived from human data. Simulations of apical pacing with a cycle length (CL) of 750 ms were performed while APD at 90% repolarization (APD90) and repolarization time (activation time + APD90) were calculated. Reentrant arrhythmia was induced by rapidly pacing the LV with a 2 mm wide apicobasal line electrode at CL≤265 ms. The simulations were repeated in the presence of diffuse myocardial fibrosis derived from the MRI patient data, and repeated again with an apicobasal repolarization gradient reduced by ≈35%. Ventricular mass was calculated as the product of model volume and myocardial density.

Results

In hypertrophic ventricles with default apicobasal repolarization gradient, rapid pacing led to reentry that lasted longer in the female (&gt;25 s) than in the male (1.0 s) ventricles. The duration of reentry increased in the presence of fibrosis in the male ventricles (male: 6.6 s; female &gt;25 s) and decreased in the female ventricles when the apicobasal repolarization gradient was reduced (Figure). Sustained reentry (&gt;25 s) was not observed in the ventricles with regressed hypertrophy, regardless of the presence of fibrosis or repolarization gradient. Average APD90 was longer in the female ventricles both before (male: 267±35 ms; female: 318±34 ms) and after (male: 272±35 ms; female: 328±38 ms) AVR. The female heart had lower myocardial mass than the male heart both before (280 vs. 406 g) and after AVR (227 vs. 310 g).

Conclusion

Longer APD and steeper apicobasal repolarization gradients in hypertrophic female ventricles resulted in sustained arrhythmia, despite their lower mass when compared to males. The absence of sustained reentry in the post-AVR models indicates a potential electrophysiological benefit of AVR and encourages future clinical studies to quantify sex-specific differences in spatial repolarization heterogeneity.

Funding Acknowledgement

Type of funding sources: Public grant(s) – EU funding. Main funding source(s): European Union's Horizon 2020 research and innovation program under the ERA-NET co-fund action No. 680969 with ANR (ERA-CVD SICVALVES) and Agence Nationale de la Recherche

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

A fast MR-thermometry method for quantitative assessment of temperature increase near an implanted wire

PURPOSE

To propose a MR-thermometry method and associated data processing technique to predict the maximal RF-induced temperature increase near an implanted wire for any other MRI sequence.

METHODS

A dynamic single shot echo planar imaging sequence was implemented that interleaves acquisition of several slices every second and an energy deposition module with adjustable parameters. Temperature images were processed in real time and compared to invasive fiber-optic measurements to assess accuracy of the method. The standard deviation of temperature was measured in gel and in vivo in the human brain of a volunteer. Temperature increases were measured for different RF exposure levels in a phantom containing an inserted wire and then a MR-conditional pacemaker lead. These calibration data set were fitted to a semi-empirical model allowing estimation of temperature increase of other acquisition sequences.

RESULTS

The precision of the measurement obtained after filtering with a 1.6x1.6 mm2 in plane resolution was 0.2°C in gel, as well as in the human brain. A high correspondence was observed with invasive temperature measurements during RF-induced heating (0.5°C RMSE for a 11.5°C temperature increase). Temperature rises of 32.4°C and 6.5°C were reached at the tip of a wire and of a pacemaker lead, respectively. After successful fitting of temperature curves of the calibration data set, temperature rise predicted by the model was in good agreement (around 5% difference) with measured temperature by a fiber optic probe, for three other MRI sequences.

CONCLUSION

This method proposes a rapid and reliable quantification of the temperature rise near an implanted wire. Calibration data set and resulting fitting coefficients can be used to estimate temperature increase for any MRI sequence as function of its power and duration.

Assessment of left ventricle magnetic resonance temperature stability in patients in the presence of arrhythmias

BACKGROUND

Magnetic resonance (MR) thermometry allows visualization of lesion formation in real-time during cardiac radiofrequency (RF) ablation. The present study was performed to evaluate the precision of MR thermometry without RF heating in patients exhibiting cardiac arrhythmia in a clinical setting. The evaluation relied on quantification of changes in temperature measurements caused by noise and physiological motion.

METHODS

Fourteen patients referred for cardiovascular magnetic resonance imaging underwent an extra sequence to test the temperature mapping stability during free-breathing acquisition. Phase images were acquired using a multi-slice, cardiac-triggered, single-shot echo planar imaging sequence. Temperature maps were calculated and displayed in real-time while the electrocardiogram (ECG) was recorded. The precision of temperature measurement was assessed by measuring the temporal standard deviation and temporal mean of consecutive temperature maps over a period of three minutes. The cardiac cycle was analyzed from ECG recordings to quantify the impact of arrhythmia events on the precision of temperature measurement. Finally, two retrospective strategies were tested to remove acquisition dynamics related either to arrhythmia events or sudden breathing motion.

RESULTS

ECG synchronization allowed categorization of inter-beat intervals (RR) into distinct beat morphologies. Five patients were in stable sinus rhythm, while nine patients showed irregular RR intervals due to ectopic beats. An average temporal standard deviation of temperature of 1.6°C was observed in patients under sinus rhythm with a frame rate corresponding to the heart rate of the patient. The temporal standard deviation rose to 2.5°C in patients with arrhythmia. The retrospective rejection strategies increased the temperature precision measurement while maintaining a sufficient frame rate.

CONCLUSIONS

Our results indicated that real-time cardiac MR thermometry shows good precision in patients under clinical conditions, even in the presence of arrhythmia. By providing real-time visualization of temperature distribution within the myocardium during RF delivery, MR thermometry could prevent insufficient or excessive heating and thus improve safety and efficacy.

High-resolution 3D diffusion tensor MRI of anesthetized rhesus macaque brain at 3T

Diffusion Magnetic Resonance Imaging (dMRI) has been widely used to investigate human brain microstructure and connectivity and its abnormalities in a variety of brain deficits, whether acute, neurodevelopmental or neurodegenerative. However, the biological interpretation and validation of dMRI data modelling is still a crucial challenge in the field. In this respect, achieving high spatial resolution in-vivo dMRI in the non-human primate to compare these observations both with human dMRI on the one hand and 'ground truth' microstructural and histological data on the other hand is of outmost importance. Here, we developed a dMRI pulse sequence based on 3D-multishot Echo Planar Imaging (3D-msEPI) on a 3T human clinical scanner. We demonstrate the feasibility of cerebral dMRI at an isotropic resolution of 0.5 mm in 4 anesthetized macaque monkeys. The added value of the high-resolution dMRI is illustrated by focusing on two aspects. First, we show an enhanced descriptive power of the fine substructure of the hippocampus. Second, we show a more physiological description of the interface between cortex grey matter, superficial and deep white matter. Overall, the high spatial resolution dMRI acquisition method proposed in this study is a significant achievement with respect to the state of the art of dMRI on anesthetized monkeys. This study highlights also the potential of very high-resolution dMRI to precisely capture the microstructure of thin cerebral structures such as the hippocampus and superficial white matter.

Real-time 3D ultrasound based motion tracking for the treatment of mobile organs with MR-guided high-intensity focused ultrasound

INTRODUCTION

Magnetic resonance-guided high-intensity focused ultrasound (MRgHIFU) treatments of mobile organs require locking the HIFU beam on the targeted tissue to maximise heating efficiency. We propose to use a standalone 3 D ultrasound (US)-based motion correction technique using the HIFU transducer in pulse-echo mode. Validation of the method was performed in vitro and in vivo in the liver of pig under MR-thermometry.

METHODS

3 D-motion estimation was implemented using ultrasonic speckle-tracking between consecutive acquisitions. Displacement was estimated along four sub-apertures of the HIFU transducer by computing the normalised cross-correlation of backscattered signals followed by a triangulation algorithm. The HIFU beam was steered accordingly and energy was delivered under real-time MR-thermometry (using the proton resonance frequency shift method with online motion compensation and correction of associated susceptibility artefacts). An MR-navigator echo was used to assess the quality of the US-based motion correction.

RESULTS

Displacement estimations from US measurements were in good agreement with 1 D MR-navigator echo readings. In vitro, the maximum temperature increase was improved by 37% as compared to experiments performed without motion correction and temperature distribution remained much more focussed. Similar results were reported in vivo, with an increase of 35% on the maximum temperature using this US-based HIFU target locking.

CONCLUSION

This standalone 3D US-based motion correction technique is robust and allows maintaining the HIFU focal spot in the presence of motion without adding any burden or complexity to MR thermal imaging. In vitro and in vivo results showed about 35% improvement in heating efficiency when focus position was locked on the target using the proposed technique.

Feasibility of real-time MR thermal dose mapping for predicting radiofrequency ablation outcome in the myocardium in vivo

BACKGROUND

Clinical treatment of cardiac arrhythmia by radiofrequency ablation (RFA) currently lacks quantitative and precise visualization of lesion formation in the myocardium during the procedure. This study aims at evaluating thermal dose (TD) imaging obtained from real-time magnetic resonance (MR) thermometry on the heart as a relevant indicator of the thermal lesion extent.

METHODS

MR temperature mapping based on the Proton Resonance Frequency Shift (PRFS) method was performed at 1.5 T on the heart, with 4 to 5 slices acquired per heartbeat. Respiratory motion was compensated using navigator-based slice tracking. Residual in-plane motion and related magnetic susceptibility artifacts were corrected online. The standard deviation of temperature was measured on healthy volunteers (N = 5) in both ventricles. On animals, the MR-compatible catheter was positioned and visualized in the left ventricle (LV) using a bSSFP pulse sequence with active catheter tracking. Twelve MR-guided RFA were performed on three sheep in vivo at various locations in left ventricle (LV). The dimensions of the thermal lesions measured on thermal dose images, on 3D T1-weighted (T1-w) images acquired immediately after the ablation and at gross pathology were correlated.

RESULTS

MR thermometry uncertainty was 1.5 °C on average over more than 96% of the pixels covering the left and right ventricles, on each volunteer. On animals, catheter repositioning in the LV with active slice tracking was successfully performed and each ablation could be monitored in real-time by MR thermometry and thermal dosimetry. Thermal lesion dimensions on TD maps were found to be highly correlated with those observed on post-ablation T1-w images (R = 0.87) that also correlated (R = 0.89) with measurements at gross pathology.

CONCLUSIONS

Quantitative TD mapping from real-time rapid CMR thermometry during catheter-based RFA is feasible. It provides a direct assessment of the lesion extent in the myocardium with precision in the range of one millimeter. Real-time MR thermometry and thermal dosimetry may improve safety and efficacy of the RFA procedure by offering a reliable indicator of therapy outcome during the procedure.

Combination of principal component analysis and optical-flow motion compensation for improved cardiac MR thermometry

The use of magnetic resonance (MR) thermometry for the monitoring of thermal ablation is rapidly expanding. However, this technique remains challenging for the monitoring of the treatment of cardiac arrhythmia by radiofrequency ablation due to the heart displacement with respiration and contraction. Recent studies have addressed this problem by compensating in-plane motion in real-time with optical-flow based tracking technique. However, these algorithms are sensitive to local variation of signal intensity on magnitude images associated with tissue heating. In this study, an optical-flow algorithm was combined with a principal component analysis method to reduce the impact of such effects. The proposed method was integrated to a fully automatic cardiac MR thermometry pipeline, compatible with a future clinical workflow. It was evaluated on nine healthy volunteers under free breathing conditions, on a phantom and in vivo on the left ventricle of a sheep. The results showed that local intensity changes in magnitude images had lower impact on motion estimation with the proposed method. Using this strategy, the temperature mapping accuracy was significantly improved.

Improved cardiac magnetic resonance thermometry and dosimetry for monitoring lesion formation during catheter ablation

PURPOSE

A new real-time MR-thermometry pipeline was developed to measure multiple temperature images per heartbeat with 1.6×1.6×3 mm³ spatial resolution. The method was evaluated on 10 healthy volunteers and during radiofrequency ablation (RFA) in sheep.

METHODS

Multislice, electrocardiogram-triggered, echo-planar imaging was combined with parallel imaging, under free breathing conditions. In-plane respiratory motion was corrected on magnitude images by an optical flow algorithm. Motion-related susceptibility artifacts were compensated on phase images by an algorithm based on Principal Component Analysis. Correction of phase drift and temporal filter were included in the pipeline implemented in the Gadgetron framework. Contact electrograms were recorded simultaneously with MR thermometry by an MR-compatible ablation catheter.

RESULTS

The temporal standard deviation of temperature in the left ventricle remained below 2 °C on each volunteer. In sheep, focal heated regions near the catheter tip were observed on temperature images (maximal temperature increase of 38 °C) during RFA, with contact electrograms of acceptable quality. Thermal lesion dimensions at gross pathology were in agreement with those observed on thermal dose images.

CONCLUSION

This fully automated MR thermometry pipeline (five images/heartbeat) provides direct assessment of lesion formation in the heart during catheter-based RFA, which may improve treatment of cardiac arrhythmia by ablation. Magn Reson Med 77:673-683, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Predictive atomic resolution descriptions of intrinsically disordered hTau40 and α-synuclein in solution from NMR and small angle scattering

The development of molecular descriptions of intrinsically disordered proteins (IDPs) is essential for elucidating conformational transitions that characterize common neurodegenerative disorders. We use nuclear magnetic resonance, small angle scattering, and molecular ensemble approaches to characterize the IDPs Tau and α-synuclein. Ensemble descriptions of IDPs are highly underdetermined due to the inherently large number of degrees of conformational freedom compared with available experimental measurements. Using extensive cross-validation we show that five different types of independent experimental parameters are predicted more accurately by selected ensembles than by statistical coil descriptions. The improvement increases in regions whose local sampling deviates from statistical coil, validating the derived conformational description. Using these approaches we identify enhanced polyproline II sampling in aggregation-nucleation sites, supporting suggestions that this region of conformational space is important for aggregation.

Phosphorylation of human Tau protein by microtubule affinity-regulating kinase 2

Tau protein plays an important role in neuronal physiology and Alzheimer's neurodegeneration. Its abilities to aggregate abnormally, to bind to microtubules (MTs), and to promote MT assembly are all influenced by phosphorylation. Phosphorylation of serine residues in the KXGS motifs of Tau's repeat domain, crucial for MT interactions and aggregation, is facilitated most efficiently by microtubule-associated protein/microtubule affinity-regulating kinases (MARKs). Here we applied high-resolution nuclear magnetic resonance analysis to study the kinetics of phosphorylation of Tau by MARK2 and its impact on the structure and microtubule binding of Tau. We demonstrate that MARK2 binds to the N-terminal tail of Tau and selectively phosphorylates three major and five minor serine residues in the repeat domain and C-terminal tail. Structural changes induced by phosphorylation of Tau by MARK2 are highly localized in the proximity of the phosphorylation site and do not affect the global conformation, in contrast to phosphorylation in the proline-rich region. Furthermore, single-residue analysis of binding of Tau to MTs provides support for a model in which Tau's hot spots of MT interaction bind independently of each other and are differentially affected by phosphorylation.

Flexible-meccano: a tool for the generation of explicit ensemble descriptions of intrinsically disordered proteins and their associated experimental observables

MOTIVATION

Intrinsically disordered proteins (IDPs) represent a significant fraction of the human proteome. The classical structure function paradigm that has successfully underpinned our understanding of molecular biology breaks down when considering proteins that have no stable tertiary structure in their functional form. One convenient approach is to describe the protein in terms of an equilibrium of rapidly inter-converting conformers. Currently, tools to generate such ensemble descriptions are extremely rare, and poorly adapted to the prediction of experimental data.

RESULTS

We present flexible-meccano-a highly efficient algorithm that generates ensembles of molecules, on the basis of amino acid-specific conformational potentials and volume exclusion. Conformational sampling depends uniquely on the primary sequence, with the possibility of introducing additional local or long-range conformational propensities at an amino acid-specific resolution. The algorithm can also be used to calculate expected values of experimental parameters measured at atomic or molecular resolution, such as nuclear magnetic resonance (NMR) and small angle scattering, respectively. We envisage that flexible-meccano will be useful for researchers who wish to compare experimental data with those expected from a fully disordered protein, researchers who see experimental evidence of deviation from 'random coil' behaviour in their protein, or researchers who are interested in working with a broad ensemble of conformers representing the flexibility of the IDP of interest.

AVAILABILITY

A fully documented multi-platform executable is provided, with examples, at http://www.ibs.fr/science-213/scientific-output/software/flexible-meccano/

CONTACT

martin.blackledge@ibs.fr.

Towards a robust description of intrinsic protein disorder using nuclear magnetic resonance spectroscopy

In order to understand the conformational behaviour of Intrinsically Disordered Proteins (IDPs), it is essential to develop a molecular representation of the partially folded state. Due to the very large number of degrees of conformational freedom available to such a disordered system, this problem is highly underdetermined. Characterisation therefore requires extensive experimental data, and novel analytical tools are required to exploit the specific conformational sensitivity of different experimental parameters. In this review we concentrate on the use of nuclear magnetic resonance (NMR) spectroscopy for the study of conformational behaviour of IDPs at atomic resolution. Each experimental NMR parameter is sensitive to different aspects of the structural and dynamic behaviour of the disordered state and requires specific consideration of the relevant averaging properties of the physical interaction. In this review we present recent advances in the description of disordered proteins and the selection of representative ensembles on the basis of experimental data using statistical coil sampling from flexible-meccano and ensemble selection using ASTEROIDS. Using these tools we aim to develop a unified molecular representation of the disordered state, combining complementary data sets to extract a meaningful description of the conformational behaviour of the protein.

NMR characterization of long-range order in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are predicted to represent a significant fraction of the human genome, and the development of meaningful molecular descriptions of these proteins remains a key challenge for contemporary structural biology. In order to describe the conformational behavior of IDPs, a molecular representation of the disordered state based on diverse sources of structural data that often exhibit complex and very different averaging behavior is required. In this study, we propose a combination of paramagnetic relaxation enhancements (PREs) and residual dipolar couplings (RDCs) to define both long-range and local structural features of IDPs in solution. We demonstrate that ASTEROIDS, an ensemble selection algorithm, faithfully reproduces intramolecular contacts, even in the presence of highly diffuse, ill-defined target interactions. We also show that explicit modeling of spin-label mobility significantly improves the reproduction of experimental PRE data, even in the case of highly disordered proteins. Prediction of the effects of transient long-range contacts on RDC profiles reveals that weak intramolecular interactions can induce a severe distortion of the profiles that compromises the description of local conformational sampling if it is not correctly taken into account. We have developed a solution to this problem that involves efficiently combining RDC and PRE data to simultaneously determine long-range and local structure in highly flexible proteins. This combined analysis is shown to be essential for the accurate interpretation of experimental data from alpha-synuclein, an important IDP involved in human neurodegenerative disease, confirming the presence of long-range order between distant regions in the protein.

Quantitative description of backbone conformational sampling of unfolded proteins at amino acid resolution from NMR residual dipolar couplings

An atomic resolution characterization of the structural properties of unfolded proteins that explicitly invokes the highly dynamic nature of the unfolded state will be extremely important for the development of a quantitative understanding of the thermodynamic basis of protein folding and stability. Here we develop a novel approach using residual dipolar couplings (RDCs) from unfolded proteins to determine conformational behavior on an amino acid specific basis. Conformational sampling is described in terms of ensembles of structures selected from a large pool of conformers. We test this approach, using extensive simulation, to determine how well the fitting of RDCs to reduced conformational ensembles containing few copies of the molecule can correctly reproduce the backbone conformational behavior of the protein. Having established approaches that allow accurate mapping of backbone dihedral angle conformational space from RDCs, we apply these methods to obtain an amino acid specific description of ubiquitin denatured in 8 M urea at pH 2.5. Cross-validation of data not employed in the fit verifies that an ensemble size of 200 structures is appropriate to characterize the highly fluctuating backbone. This approach allows us to identify local conformational sampling properties of urea-unfolded ubiquitin, which shows that the backbone sampling of certain types of charged or polar amino acids, in particular threonine, glutamic acid, and arginine, is affected more strongly by urea binding than amino acids with hydrophobic side chains. In general, the approach presented here establishes robust procedures for the study of all denatured and intrinsically disordered states.