scholarly journals Comparison of seven modelling algorithms for GABA-edited 1H-MRS

2021 ◽  
Author(s):  
Alexander R. Craven ◽  
Pallab K. Bhattacharyya ◽  
William T. Clarke ◽  
Ulrike Dydak ◽  
Richard A. E. Edden ◽  
...  

Edited MRS sequences are widely used for studying GABA in the human brain. Several algorithms are available for modelling these data, deriving metabolite concentration estimates through peak fitting or a linear combination of basis spectra. The present study compares seven such algorithms, using data obtained in a large multi-site study. GABA-edited (GABA+, TE = 68 ms MEGA-PRESS) data from 222 subjects at 20 sites were processed via a standardised pipeline, before modelling with FSL-MRS, Gannet, AMARES, QUEST, LCModel, Osprey and Tarquin, using standardised vendor-specific basis sets (for GE, Philips and Siemens) where appropriate. After referencing metabolite estimates (to water or creatine), systematic differences in scale were observed between datasets acquired on different vendors' hardware, presenting across algorithms. Scale differences across algorithms were also observed. Using the correlation between metabolite estimates and voxel tissue fraction as a benchmark, most algorithms were found to be similarly effective in detecting differences in GABA+. An inter-class correlation across all algorithms showed single-rater consistency for GABA+ estimates of around 0.38, indicating moderate agreement. Upon inclusion of a basis set component explicitly modelling the macromolecule signal underlying the observed 3.0 ppm GABA peaks, single-rater consistency improved to 0.44. Correlation between discrete pairs of algorithms varied, and was concerningly weak in some cases. Our findings highlight the need for consensus on appropriate modelling parameters across different algorithms, and for detailed reporting of the parameters adopted in individual studies to ensure reproducibility and meaningful comparison of outcomes between different studies.

2021 ◽  
Author(s):  
Yulu Song ◽  
Helge J Zollner ◽  
Steve C.N. Hui ◽  
Georg Oeltzschner ◽  
James J Prisciandaro ◽  
...  

Purpose: Two main approaches are used for spectral analysis of edited data: simple peak fitting and linear combination modeling (LCM) with a simulated basis set. Recent consensus recommended LCM as the method of choice for the spectral analysis of edited data. The aim of this study is to compare the performance of simple peak fitting and LCM in a test-retest dataset, hypothesizing that the more sophisticated LCM approach will improve quantification of HERMES data compared with simple peak fitting. Methods: A test-retest dataset was re-analyzed using Gannet (simple peak fitting) and Osprey (LCM). These data were obtained from the dorsal anterior cingulate cortex of twelve healthy volunteers, with TE 80 ms for HERMES and TE 120 ms for MEGA-PRESS of glutathione (GSH). Within-subject coefficients of variance (CVs) were calculated to quantify between-scan reproducibility of each metabolite estimate. Results: The reproducibility of HERMES GSH estimates was substantially improved using LCM compared to simple peak fitting, from a CV of 19.0% to 9.9%. For MEGA-PRESS data, the GSH reproducibility was similar using LCM and simple peak fitting, with CVs of 7.3% and 8.8% respectively. Conclusion: Linear combination modeling with simulated basis functions substantially improves the reproducibility of GSH quantification for HERMES data.


Algorithms ◽  
2020 ◽  
Vol 13 (5) ◽  
pp. 120
Author(s):  
Omid Bazgir ◽  
Eric Walden ◽  
Brian Nutter ◽  
Sunanda Mitra

Developing tools for precise quantification of brain metabolites using magnetic resonance spectroscopy (MRS) is an active area of research with broad application in non-invasive neurodegenerative disease studies. The tools are mainly developed based on black box (data-driven), or basis sets approaches. In this study, we offer a multi-stage framework that integrates data-driven and basis sets methods. We first use truncated Hankel singular value decomposition (HSVD) to decompose free induction decay (FID) signals into single tone FIDs, as the data-driven stage. Subsequently, single tone FIDs are clustered into basis sets while using initialized K-means with prior knowledge of the metabolites, as the basis set stage. The generated basis sets are fitted with the magnetic resonance (MR) spectra while using a linear constrained least square, and then the metabolite concentration is calculated. Prior to using our proposed multi-stage approach, a sequence of preprocessing blocks: water peak removal, phase correction, and baseline correction (developed in house) are used.


2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed


2019 ◽  
Author(s):  
Danilo Carmona ◽  
Pablo Jaque ◽  
Esteban Vöhringer-Martinez

<div><div><div><p>Peroxides play a central role in many chemical and biological pro- cesses such as the Fenton reaction. The relevance of these compounds lies in the low stability of the O–O bond which upon dissociation results in radical species able to initiate various chemical or biological processes. In this work, a set of 64 DFT functional-basis set combinations has been validated in terms of their capability to describe bond dissociation energies (BDE) for the O–O bond in a database of 14 ROOH peroxides for which experimental values ofBDE are available. Moreover, the electronic contributions to the BDE were obtained for four of the peroxides and the anion H2O2− at the CBS limit at CCSD(T) level with Dunning’s basis sets up to triple–ζ quality provid- ing a reference value for the hydrogen peroxide anion as a model. Almost all the functionals considered here yielded mean absolute deviations around 5.0 kcal mol−1. The smallest values were observed for the ωB97 family and the Minnesota M11 functional with a marked basis set dependence. Despite the mean deviation, order relations among BDE experimental values of peroxides were also considered. The ωB97 family was able to reproduce the relations correctly whereas other functionals presented a marked dependence on the chemical nature of the R group. Interestingly, M11 functional did not show a very good agreement with the established order despite its good performance in the mean error. The obtained results support the use of similar validation strategies for proper prediction of BDE or other molecular properties by DF Tmethods in subsequent related studies.</p></div></div></div>


2019 ◽  
Author(s):  
Kamal Batra ◽  
Stefan Zahn ◽  
Thomas Heine

<p>We thoroughly benchmark time-dependent density- functional theory for the predictive calculation of UV/Vis spectra of porphyrin derivatives. With the aim to provide an approach that is computationally feasible for large-scale applications such as biological systems or molecular framework materials, albeit performing with high accuracy for the Q-bands, we compare the results given by various computational protocols, including basis sets, density-functionals (including gradient corrected local functionals, hybrids, double hybrids and range-separated functionals), and various variants of time-dependent density-functional theory, including the simplified Tamm-Dancoff approximation. An excellent choice for these calculations is the range-separated functional CAM-B3LYP in combination with the simplified Tamm-Dancoff approximation and a basis set of double-ζ quality def2-SVP (mean absolute error [MAE] of ~0.05 eV). This is not surpassed by more expensive approaches, not even by double hybrid functionals, and solely systematic excitation energy scaling slightly improves the results (MAE ~0.04 eV). </p>


2007 ◽  
Vol 62 (12) ◽  
pp. 711-715 ◽  
Author(s):  
Ahmad Seif ◽  
Mahmoud Mirzaei ◽  
Mehran Aghaie ◽  
Asadollah Boshra

Density functional theory (DFT) calculations were performed to calculate the electric field gradient (EFG) tensors at the sites of aliminium (27Al) and nitrogen (14N) nuclei in an 1 nm of length (6,0) single-walled aliminium nitride nanotube (AlNNT) in three forms of the tubes, i. e. hydrogencapped, aliminium-terminated and nitrogen-terminated as representatives of zigzag AlNNTs. At first, each form was optimized at the level of the Becke3,Lee-Yang-Parr (B3LYP) method, 6-311G∗∗ basis set. After, the EFG tensors were calculated at the level of the B3LYP method, 6-311++G∗∗ and individual gauge for localized orbitals (IGLO-II and IGLO-III) types of basis sets in each of the three optimized forms and were converted to experimentally measurable nuclear quadrupole resonance (NQR) parameters, i. e. quadrupole coupling constant (qcc) and asymmetry parameter (ηQ). The evaluated NQR parameters revealed that the considered model of AlNNT can be divided into four equivalent layers with similar electrostatic properties.With the exception of Al-1, all of the three other Al layers have almost the same properties, however, N layers show significant differences in the magnitudes of the NQR parameters in the length of the nanotube. Furthermore, the evaluated NQR parameters of Al-1 in the Al-terminated form and N-1 in the N-terminated form revealed the different roles of Al (base agent) and of N (acid agent) in AlNNT. All the calculations were carried out using the GAUSSIAN 98 package program.


2018 ◽  
Vol 71 (4) ◽  
pp. 238 ◽  
Author(s):  
Manoj K. Kesharwani ◽  
Amir Karton ◽  
Nitai Sylvetsky ◽  
Jan M. L. Martin

The S66 benchmark for non-covalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 ‘high-level corrections’ are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66×8 problem set in particular, and for accurate calculations on non-covalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, whereas half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.


1985 ◽  
Vol 63 (7) ◽  
pp. 1487-1491 ◽  
Author(s):  
Giuseppe Del Re ◽  
Sándor Fliszár ◽  
Michel Comeau ◽  
Claude Mijoule

Net charges and valence AO's for ammonia, methylamine, dimethylamine, and trimethylamine were calculated using extended basis sets. Superposition effects, evaluated by replacing Pople's standard 6-31G* basis by an extended form in which the basis of the ammonia H atoms and of the methyl groups of trimethylamine are retained in the treatment of each molecule, indicate that the quality of the treatment of amine nitrogen atoms is strongly dependent on the number of methyl groups. A new, augmented basis is proposed for the hydrogens, which appears to be reasonably well balanced: comparison with familiar (e.g., 6-31G*) calculations illustrates in what manner the treatment of nitrogen is worsened when even just one methyl group is replaced by hydrogen unless the impoverishment of the basis is suitably taken care of.


Sign in / Sign up

Export Citation Format

Share Document