Three-dimensional representation of the pure electric-dipole and the mixed first hyperpolarizabilities: The modified unit sphere representation.

In this work, the theory of the modified unit sphere representation (mUSR) has been proposed as a computational tool suitable for the three-dimensional representation of the pure electric-dipole [ ß λ µ ν ( - 2 ω ; ω , ω ) $$ {\beta}_{\lambda \mu \nu}\left(-2\omega; \omega, \omega \right) $$ ] as well as of the mixed electric-dipole/magnetic-dipole [ α J λ µ ν ( - 2 ω ; ω , ω ) $$ {}^{\alpha }{J}_{\lambda \mu \nu}\left(-2\omega; \omega, \omega \right) $$ and ß J λ µ ν ( - 2 ω ; ω , ω ) $$ {}^{\beta }{J}_{\lambda \mu \nu}\left(-2\omega; \omega, \omega \right) $$ ] or electric-dipole/electric-quadrupole [ α K λ µ ν o ( - 2 ω ; ω , ω ) $$ {}^{\alpha }{K}_{\lambda \mu \nu o}\left(-2\omega; \omega, \omega \right) $$ and ß K λ µ ν o ( - 2 ω ; ω , ω ) $$ {}^{\beta }{K}_{\lambda \mu \nu o}\left(-2\omega; \omega, \omega \right) $$ ] first hyperpolarizabilities. These five quantities are Cartesian tensors and they are responsible for the chiral signal in the chiroptical version of the hyper-Rayleigh scattering (HRS) spectroscopy, namely the HRS optical activity (HRS-OA) spectroscopy. For the first time, for each hyperpolarizability, alongside with the three-dimensional representation of the whole (i.e., reducible) Cartesian tensors, the mUSRs are developed for each of the irreducible Cartesian tensors (ICTs) that constitute them. This scheme has been applied to a series of three (chiral) hexahelicene molecules containing different degrees of electron-withdrawing (quinone) groups and characterized by the same (positive) handedness. For these molecules, the mUSR shows that, upon substitution, the most remarkable qualitative and semi-quantitative (enhancement of the molecular responses) effects are obtained for the pure electric-dipole and for the mixed electric-dipole/magnetic-dipole hyperpolarizabilities.

Free-breathing, fat-corrected T1 mapping of the liver with stack-of-stars MRI, and joint estimation of T1, PDFF, R 2 * , and B 1 + .

PURPOSE: Quantitative T1 mapping has the potential to replace biopsy for noninvasive diagnosis and quantitative staging of chronic liver disease. Conventional T1 mapping methods are confounded by fat and B 1 + $$ {B}_1^{+} $$ inhomogeneities, resulting in unreliable T1 estimations. Furthermore, these methods trade off spatial resolution and volumetric coverage for shorter acquisitions with only a few images obtained within a breath-hold. This work proposes a novel, volumetric (3D), free-breathing T1 mapping method to account for multiple confounding factors in a single acquisition. THEORY AND METHODS: Free-breathing, confounder-corrected T1 mapping was achieved through the combination of non-Cartesian imaging, magnetization preparation, chemical shift encoding, and a variable flip angle acquisition. A subspace-constrained, locally low-rank image reconstruction algorithm was employed for image reconstruction. The accuracy of the proposed method was evaluated through numerical simulations and phantom experiments with a T1/proton density fat fraction phantom at 3.0 T. Further, the feasibility of the proposed method was investigated through contrast-enhanced imaging in healthy volunteers, also at 3.0 T. RESULTS: The method showed excellent agreement with reference measurements in phantoms across a wide range of T1 values (200 to 1000 ms, slope = 0.998 (95% confidence interval (CI) [0.963 to 1.035]), intercept = 27.1 ms (95% CI [0.4 54.6]), r2 = 0.996), and a high level of repeatability. In vivo imaging studies demonstrated moderate agreement (slope = 1.099 (95% CI [1.067 to 1.132]), intercept = -96.3 ms (95% CI [-82.1 to -110.5]), r2 = 0.981) compared to saturation recovery-based T1 maps. CONCLUSION: The proposed method produces whole-liver, confounder-corrected T1 maps through simultaneous estimation of T1, proton density fat fraction, and B 1 + $$ {B}_1^{+} $$ in a single, free-breathing acquisition and has excellent agreement with reference measurements in phantoms.

A decay-modeled compressed sensing reconstruction approach for non-Cartesian hyperpolarized 129Xe MRI.

PURPOSE: Hyperpolarized 129Xe MRI benefits from non-Cartesian acquisitions that sample k-space efficiently and rapidly. However, their reconstructions are complex and burdened by decay processes unique to hyperpolarized gas. Currently used gridded reconstructions are prone to artifacts caused by magnetization decay and are ill-suited for undersampling. We present a compressed sensing (CS) reconstruction approach that incorporates magnetization decay in the forward model, thereby producing images with increased sharpness and contrast, even in undersampled data. METHODS: Radio-frequency, T1, and T 2 * $$ {\mathrm{T}}_2^{\ast } $$ decay processes were incorporated into the forward model and solved using iterative methods including CS. The decay-modeled reconstruction was validated in simulations and then tested in 2D/3D-spiral ventilation and 3D-radial gas-exchange MRI. Quantitative metrics including apparent-SNR and sharpness were compared between gridded, CS, and twofold undersampled CS reconstructions. Observations were validated in gas-exchange data collected from 15 healthy and 25 post-hematopoietic-stem-cell-transplant participants. RESULTS: CS reconstructions in simulations yielded images with threefold increases in accuracy. CS increased sharpness and contrast for ventilation in vivo imaging and showed greater accuracy for undersampled acquisitions. CS improved gas-exchange imaging, particularly in the dissolved-phase where apparent-SNR improved, and structure was made discernable. Finally, CS showed repeatability in important global gas-exchange metrics including median dissolved-gas signal ratio and median angle between real/imaginary components. CONCLUSION: A non-Cartesian CS reconstruction approach that incorporates hyperpolarized 129Xe decay processes is presented. This approach enables improved image sharpness, contrast, and overall image quality in addition to up-to threefold undersampling. This contribution benefits all hyperpolarized gas MRI through improved accuracy and decreased scan durations.

Impact of B 0 $$ {\mathrm{B}}_0 $$ field imperfections correction on BOLD sensitivity in 3D-SPARKLING fMRI data.

PURPOSE: Static and dynamic B 0 $$ {\mathrm{B}}_0 $$ field imperfections are detrimental to functional MRI (fMRI) applications, especially at ultra-high magnetic fields (UHF). In this work, a field camera is used to assess the benefits of retrospectively correcting B 0 $$ {\mathrm{B}}_0 $$ field perturbations on Blood Oxygen Level Dependent (BOLD) sensitivity in non-Cartesian three-dimensional (3D)-SPARKLING fMRI acquisitions. METHODS: fMRI data were acquired at 1 mm 3 $$ {}^3 $$ and for a 2.4s-TR while concurrently monitoring in real-time field perturbations using a Skope Clip-on field camera in a novel experimental setting involving a shorter TR than the required minimal TR of the field probes. Measurements of the dynamic field deviations were used along with a static Δ B 0 $$ \Delta {\mathrm{B}}_0 $$ map to retrospectively correct static and dynamic field imperfections, respectively. In order to evaluate the impact of such a correction on fMRI volumes, a comparative study was conducted on healthy volunteers. RESULTS: Correction of B 0 $$ {\mathrm{B}}_0 $$ deviations improved image quality and yielded between 20% and 30% increase in median temporal signal-to-noise ratio (tSNR).Using fMRI data collected during a retinotopic mapping experiment, we demonstrated a significant increase in sensitivity to the BOLD contrast and improved accuracy of the BOLD phase maps: 44% (resp., 159%) more activated voxels were retrieved when using a significance control level based on a p-value of 0.001 without correcting for multiple comparisons (resp., 0.05 with a false discovery rate correction). CONCLUSION: 3D-SPARKLING fMRI hugely benefits from static and dynamic B 0 $$ {\mathrm{B}}_0 $$ imperfections correction. However, the proposed experimental protocol is flexible enough to be deployed on a large spectrum of encoding schemes, including arbitrary non-Cartesian readouts.

Reducing clustering of readouts in non-Cartesian cine magnetic resonance imaging.

BACKGROUND: Non-Cartesian magnetic resonance imaging trajectories at golden angle increments have the advantage of allowing motion correction and gating using intermediate real-time reconstructions. However, when the acquired data are cardiac binned for cine imaging, trajectories can cluster together at certain heart rates (HR) causing image artifacts. Here, we demonstrate an approach to reduce clustering by inserting additional angular increments within the trajectory, and optimizing them while still allowing for intermediate reconstructions. METHODS: Three acquisition models were simulated under constant and variable HR: golden angle (Mtrd), random additional angles (Mrnd), and optimized additional angles (Mopt). The standard deviations of trajectory angular differences (STAD) were compared through their interquartile ranges (IQR) and the Kolmogorov-Smirnov test (significance level: p = 0.05). Agreement between an image reconstructed with uniform sampling and images from Mtrd, Mrnd, and Mopt was analyzed using the structural similarity index measure (SSIM). Mtrd and Mopt were compared in three adults at high, low, and no HR variability. RESULTS: STADs from Mtrd were significantly different (p < 0.05) from Mopt and Mrnd. STAD (IQR × 10-2 rad) showed that Mopt (0.5) and Mrnd (0.5) reduced clustering relative to Mtrd (1.9) at constant HR. For variable HR, Mopt (0.5) and Mrnd (0.5) outperformed Mtrd (0.9). The SSIM (IQR) showed that Mopt (0.011) produced the best image quality, followed by Mrnd (0.014), and Mtrd (0.030). Mopt outperformed Mtrd at reduced HR variability in in-vivo studies. At high HR variability, both models performed well. CONCLUSION: This approach reduces clustering in k-space and improves image quality.

Large-scale 3D non-Cartesian coronary MRI reconstruction using distributed memory-efficient physics-guided deep learning with limited training data.

OBJECT: To enable high-quality physics-guided deep learning (PG-DL) reconstruction of large-scale 3D non-Cartesian coronary MRI by overcoming challenges of hardware limitations and limited training data availability. MATERIALS AND METHODS: While PG-DL has emerged as a powerful image reconstruction method, its application to large-scale 3D non-Cartesian MRI is hindered by hardware limitations and limited availability of training data. We combine several recent advances in deep learning and MRI reconstruction to tackle the former challenge, and we further propose a 2.5D reconstruction using 2D convolutional neural networks, which treat 3D volumes as batches of 2D images to train the network with a limited amount of training data. Both 3D and 2.5D variants of the PG-DL networks were compared to conventional methods for high-resolution 3D kooshball coronary MRI. RESULTS: Proposed PG-DL reconstructions of 3D non-Cartesian coronary MRI with 3D and 2.5D processing outperformed all conventional methods both quantitatively and qualitatively in terms of image assessment by an experienced cardiologist. The 2.5D variant further improved vessel sharpness compared to 3D processing, and scored higher in terms of qualitative image quality. DISCUSSION: PG-DL reconstruction of large-scale 3D non-Cartesian MRI without compromising image size or network complexity is achieved, and the proposed 2.5D processing enables high-quality reconstruction with limited training data.

Model-based determination of the synchronization delay between MRI and trajectory data.

PURPOSE: Real-time monitoring of dynamic magnetic fields has recently become a commercially available option for measuring MRI k-space trajectories and magnetic fields induced by eddy currents in real time. However, for accurate image reconstructions, sub-microsecond synchronization between the MRI data and field dynamics (ie, k-space trajectory plus other spatially varying fields) is required. In this work, we introduce a new model-based algorithm to automatically perform this synchronization using only the MRI data and field dynamics. METHODS: The algorithm works by enforcing consistency among the MRI data, field dynamics, and receiver sensitivity profiles by iteratively alternating between convex optimizations for (a) the image and (b) the synchronization delay. A healthy human subject was scanned at 7 T using a transmit-receive coil with integrated field probes using both single-shot spiral and EPI, and reconstructions with various synchronization delays were compared with the result of the proposed algorithm. The accuracy of the algorithm was also investigated using simulations, in which the acquisition delays for simulated acquisitions were determined using the proposed algorithm and compared with the known ground truth. RESULTS: In the in vivo scans, the proposed algorithm minimized artifacts related to synchronization delay for both spiral and EPI acquisitions, and the computation time required was less than 30 s. The simulations demonstrated accuracy to within tens of nanoseconds. CONCLUSIONS: The proposed algorithm can automatically determine synchronization delays between MRI data and field dynamics measured using a field probe system.

Ultra-short T2 components imaging of the whole brain using 3D dual-echo UTE MRI with rosette k-space pattern.

PURPOSE: This study aimed to develop a new 3D dual-echo rosette k-space trajectory, specifically designed for UTE MRI applications. The imaging of the ultra-short transverse relaxation time (uT2 ) of brain was acquired to test the performance of the proposed UTE sequence. THEORY AND METHODS: The rosette trajectory was developed based on rotations of a "petal-like" pattern in the kx -ky plane, with oscillated extensions in the kz -direction for 3D coverage. Five healthy volunteers underwent 10 dual-echo 3D rosette UTE scans with various TEs. Dual-exponential complex model fitting was performed on the magnitude data to separate uT2 signals, with the output of uT2 fraction, uT2 value, and long-T2 value. RESULTS: The 3D rosette dual-echo UTE sequence showed better performance than a 3D radial UTE acquisition. More significant signal intensity decay in white matter than gray matter was observed along with the TEs. The white matter regions had higher uT2 fraction values than gray matter (10.9% ± 1.9% vs. 5.7% ± 2.4%). The uT2 value was approximately 0.10 ms in white matter . CONCLUSION: The higher uT2 fraction value in white matter compared to gray matter demonstrated the ability of the proposed sequence to capture rapidly decaying signals.

Free-breathing 2D radial cine MRI with respiratory auto-calibrated motion correction (RAMCO).

PURPOSE: To develop a free-breathing (FB) 2D radial balanced steady-state free precession cine cardiac MRI method with 100% respiratory gating efficiency using respiratory auto-calibrated motion correction (RAMCO) based on a motion-sensing camera. METHODS: The signal from a respiratory motion-sensing camera was recorded during a FB retrospectively electrocardiogram triggered 2D radial balanced steady-state free precession acquisition using pseudo-tiny-golden-angle ordering. With RAMCO, for each acquisition the respiratory signal was retrospectively auto-calibrated by applying different linear translations, using the resulting in-plane image sharpness as a criterium. The auto-calibration determines the optimal magnitude of the linear translations for each of the in-plane directions to minimize motion blurring caused by bulk respiratory motion. Additionally, motion-weighted density compensation was applied during radial gridding to minimize through-plane and non-bulk motion blurring. Left ventricular functional parameters and sharpness scores of FB radial cine were compared with and without RAMCO, and additionally with conventional breath-hold Cartesian cine on 9 volunteers. RESULTS: FB radial cine with RAMCO had similar sharpness scores as conventional breath-hold Cartesian cine and the left ventricular functional parameters agreed. For FB radial cine, RAMCO reduced respiratory motion artifacts with a statistically significant difference in sharpness scores (P < 0.05) compared to reconstructions without motion correction. CONCLUSION: 2D radial cine imaging with RAMCO allows evaluation of left ventricular functional parameters in FB with 100% respiratory efficiency. It eliminates the need for breath-holds, which is especially valuable for patients with no or impaired breath-holding capacity. Validation of the proposed method on patients is warranted.

Self-calibrated through-time spiral GRAPPA for real-time, free-breathing evaluation of left ventricular function.

PURPOSE: Through-time spiral GRAPPA is a real-time imaging technique that enables ungated, free-breathing evaluation of the left ventricle. However, it requires a separate fully-sampled calibration scan to calculate GRAPPA weights. A self-calibrated through-time spiral GRAPPA method is proposed that uses a specially designed spiral trajectory with interleaved arm ordering such that consecutive undersampled frames can be merged to form calibration data, eliminating the separate fully-sampled acquisition. THEORY AND METHODS: The proposed method considers the time needed to acquire data at all points in a GRAPPA calibration kernel when using interleaved arm ordering. Using this metric, simulations were performed to design a spiral trajectory for self-calibrated GRAPPA. Data were acquired in healthy volunteers using the proposed method and a comparison electrocardiogram-gated and breath-held cine scan. Left ventricular functional values and image quality are compared. RESULTS: A 12-arm spiral trajectory was designed with a temporal resolution of 32.72 ms/cardiac phase with an acceleration factor of 3. Functional values calculated using the proposed method and the gold-standard method were not statistically significantly different (paired t-test, p < 0.05). Image quality ratings were lower for the proposed method, with statistically significantly different ratings (Wilcoxon signed rank test, p < 0.05) for two of five image quality aspects rated (level of artifact, blood-myocardium contrast). CONCLUSIONS: A self-calibrated through-time spiral GRAPPA reconstruction can enable ungated, free-breathing evaluation of the left ventricle in 71 s. Functional values are equivalent to a gold-standard cine technique, although some aspects of image quality may be inferior due to the real-time nature of the data collection.

Stochastic optimization of three-dimensional non-Cartesian sampling trajectory.

PURPOSE: Optimizing three-dimensional (3D) k-space sampling trajectories is important for efficient MRI yet presents a challenging computational problem. This work proposes a generalized framework for optimizing 3D non-Cartesian sampling patterns via data-driven optimization. METHODS: We built a differentiable simulation model to enable gradient-based methods for sampling trajectory optimization. The algorithm can simultaneously optimize multiple properties of sampling patterns, including image quality, hardware constraints (maximum slew rate and gradient strength), reduced peripheral nerve stimulation (PNS), and parameter-weighted contrast. The proposed method can either optimize the gradient waveform (spline-based freeform optimization) or optimize properties of given sampling trajectories (such as the rotation angle of radial trajectories). Notably, the method can optimize sampling trajectories synergistically with either model-based or learning-based reconstruction methods. We proposed several strategies to alleviate the severe nonconvexity and huge computation demand posed by the large scale. The corresponding code is available as an open-source toolbox. RESULTS: We applied the optimized trajectory to multiple applications including structural and functional imaging. In the simulation studies, the image quality of a 3D kooshball trajectory was improved from 0.29 to 0.22 (NRMSE) with Stochastic optimization framework for 3D NOn-Cartesian samPling trajectorY (SNOPY) optimization. In the prospective studies, by optimizing the rotation angles of a stack-of-stars (SOS) trajectory, SNOPY reduced the NRMSE of reconstructed images from 1.19 to 0.97 compared to the best empirical method (RSOS-GR). Optimizing the gradient waveform of a rotational EPI trajectory improved participants' rating of the PNS from "strong" to "mild." CONCLUSION: SNOPY provides an efficient data-driven and optimization-based method to tailor non-Cartesian sampling trajectories.

Single-shot spiral diffusion-weighted imaging at 7T using expanded encoding with compressed sensing.

PURPOSE: The expanded encoding model incorporates spatially- and time-varying field perturbations for correction during reconstruction. To date, these reconstructions have used the conjugate gradient method with early stopping used as implicit regularization. However, this approach is likely suboptimal for low-SNR cases like diffusion or high-resolution MRI. Here, we investigate the extent that ℓ 1 $$ {\ell}_1 $$ -wavelet regularization, or equivalently compressed sensing (CS), combined with expanded encoding improves trade-offs between spatial resolution, readout time and SNR for single-shot spiral DWI at 7T. The reconstructions were performed using our open-source graphics processing unit-enabled reconstruction toolbox, "MatMRI," that allows inclusion of the different components of the expanded encoding model, with or without CS. METHODS: In vivo accelerated single-shot spirals were acquired with five acceleration factors (R) (2×-6×) and three in-plane spatial resolutions (1.5, 1.3, and 1.1 mm). From the in vivo reconstructions, we estimated diffusion tensors and computed fractional anisotropy maps. Then, simulations were used to quantitatively investigate and validate the impact of CS-based regularization on image quality when compared to a known ground truth. RESULTS: In vivo reconstructions revealed improved image quality with retainment of small features when CS was used. Simulations showed that the joint use of the expanded encoding model and CS improves accuracy of image reconstructions (reduced mean-squared error) over the range of R investigated. CONCLUSION: The expanded encoding model and CS regularization are complementary tools for single-shot spiral diffusion MRI, which enables both higher spatial resolutions and higher R.

Deep learning-assisted model-based off-resonance correction for non-Cartesian SWI.

PURPOSE: Patient-induced inhomogeneities in the static magnetic field cause distortions and blurring (off-resonance artifacts) during acquisitions with long readouts such as in SWI. Conventional versatile correction methods based on extended Fourier models are too slow for clinical practice in computationally demanding cases such as 3D high-resolution non-Cartesian multi-coil acquisitions. THEORY: Most reconstruction methods can be accelerated when performing off-resonance correction by reducing the number of iterations, compressed coils, and correction components. Recent state-of-the-art unrolled deep learning architectures could help but are generally not adapted to corrupted measurements as they rely on the standard Fourier operator in the data consistency term. The combination of correction models and neural networks is therefore necessary to reduce reconstruction times. METHODS: Hybrid pipelines using UNets were trained stack-by-stack over 99 SWI 3D SPARKLING 20-fold accelerated acquisitions at 0.6 mm isotropic resolution using different off-resonance correction methods. Target images were obtained using slow model-based corrections based on self-estimated Δ B 0 $$ \Delta {B}_0 $$ field maps. The proposed strategies, tested over 11 volumes, are compared to model-only and network-only pipelines. RESULTS: The proposed hybrid pipelines achieved scores competing with two to three times slower baseline methods, and neural networks were observed to contribute both as pre-conditioner and through inter-iteration memory by allowing more degrees of freedom over the model design. CONCLUSION: A combination of model-based and network-based off-resonance correction was proposed to significantly accelerate conventional methods. Different promising synergies were observed between acceleration factors (iterations, coils, correction) and model/network that could be expanded in the future.

Submillimeter lung MRI at 0.55 T using balanced steady-state free precession with half-radial dual-echo readout (bSTAR).

PURPOSE: To demonstrate the feasibility of high-resolution morphologic lung MRI at 0.55 T using a free-breathing balanced steady-state free precession half-radial dual-echo imaging technique (bSTAR). METHODS: Self-gated free-breathing bSTAR (TE1 /TE2 /TR of 0.13/1.93/2.14 ms) lung imaging in five healthy volunteers and a patient with granulomatous lung disease was performed using a 0.55 T MR-scanner. A wobbling Archimedean spiral pole (WASP) trajectory was used to ensure a homogenous coverage of k-space over multiple breathing cycles. WASP uses short-duration interleaves randomly tilted by a small polar angle and rotated by a golden angle about the polar axis. Data were acquired continuously over 12:50 min. Respiratory-resolved images were reconstructed off-line using compressed sensing and retrospective self-gating. Reconstructions were performed with a nominal resolution of 0.9 mm and a reduced isotropic resolution of 1.75 mm corresponding to shorter simulated scan times of 8:34 and 4:17 min, respectively. Analysis of apparent SNR was performed in all volunteers and reconstruction settings. RESULTS: The technique provided artifact-free morphologic lung images in all subjects. The short TR of bSTAR in conjunction with a field strength of 0.55 T resulted in a complete mitigation of off-resonance artifacts in the chest. Mean SNR values in healthy lung parenchyma for the 12:50 min scan were 3.6 ± 0.8 and 24.9 ± 6.2 for 0.9 mm and 1.75 mm reconstructions, respectively. CONCLUSION: This study demonstrates the feasibility of morphologic lung MRI with a submillimeter isotropic spatial resolution in human subjects with bSTAR at 0.55 T.

Design and development of a novel flexible ultra-short echo time (FUSE) sequence.

PURPOSE: To present the validation of a new Flexible Ultra-Short Echo time (FUSE) pulse sequence using a short-T2 phantom. METHODS: FUSE was developed to include a range of RF excitation pulses, trajectories, dimensionalities, and long-T2 suppression techniques, enabling real-time interchangeability of acquisition parameters. Additionally, we developed an improved 3D deblurring algorithm to correct for off-resonance artifacts. Several experiments were conducted to validate the efficacy of FUSE, by comparing different approaches for off-resonance artifact correction, variations in RF pulse and trajectory combinations, and long-T2 suppression techniques. All scans were performed on a 3 T system using an in-house short-T2 phantom. The evaluation of results included qualitative comparisons and quantitative assessments of the SNR and contrast-to-noise ratio. RESULTS: Using the capabilities of FUSE, we demonstrated that we could combine a shorter readout duration with our improved deblurring algorithm to effectively reduce off-resonance artifacts. Among the different RF and trajectory combinations, the spiral trajectory with the regular half-inc pulse achieves the highest SNRs. The dual-echo subtraction technique delivers better short-T2 contrast and superior suppression of water and agar signals, whereas the off-resonance saturation method successfully suppresses water and lipid signals simultaneously. CONCLUSION: In this work, we have validated the use of our new FUSE sequence using a short T2 phantom, demonstrating that multiple UTE acquisitions can be achieved within a single sequence. This new sequence may be useful for acquiring improved UTE images and the development of UTE imaging protocols.

Improving spreading projection algorithm for rapid k-space sampling trajectories through minimized off-resonance effects and gridding of low frequencies.

PURPOSE: Non-Cartesian MRI with long arbitrary readout directions are susceptible to off-resonance artifacts due to patient induced B 0 $$ {B}_0 $$ inhomogeneities. This results in degraded image quality with strong signal losses and blurring. Current solutions to address this issue involve correcting the off-resonance artifacts during image reconstruction or reducing inhomogeneities through improved shimming. THEORY: The recently developed SPARKLING algorithm is extended to drastically diminish off-resonance artifacts by generating temporally smooth k-space sampling patterns. For doing so, the cost function which is optimized in SPARKLING is modified using a temporal weighting factor. Additionally, oversampling of the center of k-space beyond the Nyquist criteria is prevented through the use of gridded sampling in the region, enforced with affine constraints. METHODS: Prospective k-space data was acquired at 3 T on new trajectories, and we show robustness to B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities through in silico experiments by adding Δ B 0 $$ \Delta {\mathrm{B}}_0 $$  through artificial degradation of system B 0 $$ {\mathrm{B}}_0 $$ shimming. Later on, in vivo experiments were carried out to optimize parameters of the new improvements and benchmark the gain in performance. RESULTS: The improved trajectories allowed for the recovery of signal dropouts observed on original SPARKLING acquisitions at larger B 0 $$ {\mathrm{B}}_0 $$ field inhomogeneities. Furthermore, imposing gridded sampling at the center of k-space provided improved reconstructed image quality with limited artifacts. CONCLUSION: These advancements allowed us for nearly 4 . 62 × $$ 4.62\times $$ shorter scan time compared to GRAPPA-p4x1, allowing us to reach 600 µm isotropic resolution in 3D T 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ -w imaging in just 3.3 min at 3 T with negligible degradation in image quality.

Alteration of Chain-Length Selectivity and Thermostability of Rhizopus oryzae Lipase via Virtual Saturation Mutagenesis Coupled with Disulfide Bond Design.

Rhizopus oryzae lipase (ROL) is one of the most important enzymes used in the food, biofuel, and pharmaceutical industries. However, the highly demanding conditions of industrial processes can reduce its stability and activity. To seek a feasible method to improve both the catalytic activity and the thermostability of this lipase, first, the structure of ROL was divided into catalytic and noncatalytic regions by identifying critical amino acids in the crevice-like binding pocket. Second, a mutant screening library aimed at improvement of ROL catalytic performance by virtual saturation mutagenesis of residues in the catalytic region was constructed based on Rosetta's Cartesian_ddg protocol. A double mutant, E265V/S267W (with an E-to-V change at residue 265 and an S-to-W change at residue 267), with markedly improved catalytic activity toward diverse chain-length fatty acid esters was identified. Then, computational design of disulfide bonds was conducted for the noncatalytic amino acids of E265V/S267W, and two potential disulfide bonds, S61C-S115C and E190C-E238C, were identified as candidates. Experimental data validated that the variant E265V/S267W/S61C-S115C/E190C-E238C had superior stability, with an increase of 8.5°C in the melting temperature and a half-life of 31.7 min at 60°C, 4.2-fold longer than that of the wild-type enzyme. Moreover, the variant improved the lipase activity toward five 4-nitrophenyl esters by 1.5 to 3.8 times, exhibiting a potential to modify the catalytic efficiency. IMPORTANCE Rhizopus oryzae lipase (ROL) is very attractive in biotechnology and industry as a safe and environmentally friendly biocatalyst. Functional expression of ROL in Escherichia coli facilitates effective high-throughput screening for positive variants. This work highlights a method to improve both selectivity and thermostability based on a combination of virtual saturation mutagenesis in the substrate pocket and disulfide bond prediction in the noncatalytic region. Using the method, ROL thermostability and activity to diverse 4-nitrophenyl esters could be substantially improved. The strategy of rational introduction of multiple mutations in different functional domains of the enzyme is a great prospect in the modification of biocatalysts.

Customizing Anterolateral Thigh Flap With Magnetic Resonance Angiography Differential Subsampling With Cartesian Ordering Imaging for Individualized Reconstruction of Extremity Defects.

INTRODUCTION: Magnetic resonance angiography (MRA) with the differential subsampling with cartesian ordering (DISCO) imaging technique is rarely used in anterolateral thigh (ALT) flap. In our series, MRA DISCO imaging technique is used as a tool to customize ALT flaps. The aim of this study was to report the accuracy of cutaneous perforators identified by the MRA DISCO imaging. METHODS: Nineteen patients underwent the MRA DISCO imaging for perforator mapping before the ALT flap transfer. A total of 38 ALT regions were studied on the MRA DISCO images. Flap thinning was performed under the guidance of MRA DISCO imaging. RESULTS: The lateral circumflex femoral artery (LCFA) most commonly stems from the deep femoral artery (84.2%), followed by the common femoral artery (15.8%). The average number of perforator vessels per LCFA was 10.2 ± 1.7. The distinct oblique branch was observed in 16 out of the 38 ALT regions (42.1%). Among the 19 ALT flaps harvested, 5 were septocutaneous perforator flaps and 14 musculocutaneous perforator flaps. Ten were harvested based on the descending branch, and 3 used the oblique branch as the flap vascular pedicle. In addition, the displayed course and types of perforator vessels on the DISCO images of the 18 skin flaps were consistent with the intraoperative findings, with an accuracy of 94.7%. CONCLUSIONS: The state of the cutaneous perforators of LCFA can be identified on the MRA DISCO images. The 3D-CE-MRA DISCO imaging is a practical method, which can ameliorate the design and customization of ALT flap for an individualized reconstruction.

Concomitant field compensation of spiral turbo spin-echo at 0.55 T.

OBJECTIVE: Diagnostic-quality neuroimaging methods are vital for widespread clinical adoption of low field MRI. Spiral imaging is an efficient acquisition method that can mitigate the reduced signal-to-noise ratio at lower field strengths. As concomitant field artifacts are worse at lower field, we propose a generalizable quadratic gradient-field nulling as an echo-to-echo compensation and apply it to spiral TSE at 0.55 T. MATERIALS AND METHODS: A spiral in-out TSE acquisition was developed with a compensation for concomitant field variation between spiral interleaves, by adding bipolar gradients around each readout to minimize phase differences at each refocusing pulse. Simulations were performed to characterize concomitant field compensation approaches. We demonstrate our proposed compensation method in phantoms and (n = 8) healthy volunteers at 0.55 T. RESULTS: Spiral read-outs with integrated spoiling demonstrated strong concomitant field artifacts but were mitigated using the echo-to-echo compensation. Simulations predicted a decrease of concomitant field phase RMSE between echoes of 42% using the proposed compensation. Spiral TSE improved SNR by 17.2 ± 2.3% compared to reference Cartesian acquisition. DISCUSSION: We demonstrated a generalizable approach to mitigate concomitant field artifacts for spiral TSE acquisitions via the addition of quadratic-nulling gradients, which can potentially improve neuroimaging at low-field through increased acquisition efficiency.

Utility and Usability of Two Forms of Supplemental Vibrotactile Kinesthetic Feedback for Enhancing Movement Accuracy and Efficiency in Goal-Directed Reaching.

Recent advances in wearable sensors and computing have made possible the development of novel sensory augmentation technologies that promise to enhance human motor performance and quality of life in a wide range of applications. We compared the objective utility and subjective user experience for two biologically inspired ways to encode movement-related information into supplemental feedback for the real-time control of goal-directed reaching in healthy, neurologically intact adults. One encoding scheme mimicked visual feedback encoding by converting real-time hand position in a Cartesian frame of reference into supplemental kinesthetic feedback provided by a vibrotactile display attached to the non-moving arm and hand. The other approach mimicked proprioceptive encoding by providing real-time arm joint angle information via the vibrotactile display. We found that both encoding schemes had objective utility in that after a brief training period, both forms of supplemental feedback promoted improved reach accuracy in the absence of concurrent visual feedback over performance levels achieved using proprioception alone. Cartesian encoding promoted greater reductions in target capture errors in the absence of visual feedback (Cartesian: 59% improvement; Joint Angle: 21% improvement). Accuracy gains promoted by both encoding schemes came at a cost in terms of temporal efficiency; target capture times were considerably longer (1.5 s longer) when reaching with supplemental kinesthetic feedback than without. Furthermore, neither encoding scheme yielded movements that were particularly smooth, although movements made with joint angle encoding were smoother than movements with Cartesian encoding. Participant responses on user experience surveys indicate that both encoding schemes were motivating and that both yielded passable user satisfaction scores. However, only Cartesian endpoint encoding was found to have passable usability; participants felt more competent using Cartesian encoding than joint angle encoding. These results are expected to inform future efforts to develop wearable technology to enhance the accuracy and efficiency of goal-directed actions using continuous supplemental kinesthetic feedback.