BioNMR
NMR aggregator & online community since 2003
BioNMR    
Learn or help to learn NMR - get free NMR books!
 

Go Back   BioNMR > NMR community > News from NMR blogs
Advanced Search
Home Forums Wiki NMR feeds Downloads Register Today's Posts



Jobs Groups Conferences Literature Pulse sequences Software forums Programs Sample preps Web resources BioNMR issues


Webservers
NMR processing:
MDD
NMR assignment:
Backbone:
Autoassign
MARS
UNIO Match
PINE
Side-chains:
UNIO ATNOS-Ascan
NOEs:
UNIO ATNOS-Candid
UNIO Candid
ASDP
Structure from NMR restraints:
Ab initio:
GeNMR
Cyana
XPLOR-NIH
ASDP
UNIO ATNOS-Candid
UNIO Candid
Fragment-based:
BMRB CS-Rosetta
Rosetta-NMR (Robetta)
Template-based:
GeNMR
I-TASSER
Refinement:
Amber
Structure from chemical shifts:
Fragment-based:
WeNMR CS-Rosetta
BMRB CS-Rosetta
Homology-based:
CS23D
Simshift
Torsion angles from chemical shifts:
Preditor
TALOS
Promega- Proline
Secondary structure from chemical shifts:
CSI (via RCI server)
TALOS
MICS caps, β-turns
d2D
PECAN
Flexibility from chemical shifts:
RCI
Interactions from chemical shifts:
HADDOCK
Chemical shifts re-referencing:
Shiftcor
UNIO Shiftinspector
LACS
CheckShift
RefDB
NMR model quality:
NOEs, other restraints:
PROSESS
PSVS
RPF scores
iCing
Chemical shifts:
PROSESS
CheShift2
Vasco
iCing
RDCs:
DC
Anisofit
Pseudocontact shifts:
Anisofit
Protein geomtery:
Resolution-by-Proxy
PROSESS
What-If
iCing
PSVS
MolProbity
SAVES2 or SAVES4
Vadar
Prosa
ProQ
MetaMQAPII
PSQS
Eval123D
STAN
Ramachandran Plot
Rampage
ERRAT
Verify_3D
Harmony
Quality Control Check
NMR spectrum prediction:
FANDAS
MestReS
V-NMR
Flexibility from structure:
Backbone S2
Methyl S2
B-factor
Molecular dynamics:
Gromacs
Amber
Antechamber
Chemical shifts prediction:
From structure:
Shiftx2
Sparta+
Camshift
CH3shift- Methyl
ArShift- Aromatic
ShiftS
Proshift
PPM
CheShift-2- Cα
From sequence:
Shifty
Camcoil
Poulsen_rc_CS
Disordered proteins:
MAXOCC
Format conversion & validation:
CCPN
From NMR-STAR 3.1
Validate NMR-STAR 3.1
NMR sample preparation:
Protein disorder:
DisMeta
Protein solubility:
camLILA
ccSOL
Camfold
camGroEL
Zyggregator
Isotope labeling:
UPLABEL
Solid-state NMR:
sedNMR


Reply
 
Thread Tools Search this Thread Rate Thread Display Modes
  #1  
Old 01-05-2014, 04:14 PM
nmrlearner's Avatar
Senior Member
 
Join Date: Jan 2005
Posts: 23,714
Points: 193,617, Level: 100
Points: 193,617, Level: 100 Points: 193,617, Level: 100 Points: 193,617, Level: 100
Level up: 0%, 0 Points needed
Level up: 0% Level up: 0% Level up: 0%
Activity: 50.7%
Activity: 50.7% Activity: 50.7% Activity: 50.7%
Last Achievements
Award-Showcase
NMR Credits: 0
NMR Points: 193,617
Downloads: 0
Uploads: 0
Default Reference Deconvolution

Reference Deconvolution

Introduction

In an idealized situation, according to quantum mechanics theory, NMR transitions in liquid state and excluding dynamic effects such as chemical exchange are of Lorentzian shape [1]. In practice NMR lineshapes are never pure Lorentzians due to a number of reasons [2], ranging from magnetic field inhomogeneity and magnetic field noise to sample temperature gradients, sample spinning or FID weighting, to cite a few.
Another important property that might affect significantly the final observed lineshape is the following: Even in molecules of modest size the number of distinct peaks might be thousands times smaller than that of quantum transitions. As a simple example, the number of transitions of a molecule containing 15 would be 245760 whereas only a few hundreds of peaks would be observed in the spectrum. As a result, an NMR peak is actually an envelope of a distribution of myriads of Lorentzians and its shape is dominated by the coupling pattern of the spin system.
Whilst this kind of line broadening affects signals differently across the spectrum and is very difficult (or even impossible) to resolve by post-processing operations, there are many other distortions that affect all resonances in the spectrum in the same way. These include lineshape distortions caused by poorly shimmed samples and they can be removed by using a post-processing technique known as Reference Deconvolution [3]


Reference Deconvolution


This technique is used to remove the instrumental lineshape distortion by deconvolving the experimental NMR spectrum using a reference signal, usually one within the same spectrum (which should be an isolated singlet) known to be subject to the identical lineshape distortions. Finally, once the lineshape distortion is removed, the spectrum can be reconvoluted with a known lineshape, typically a Lorentzian, so that the result will be a corrected spectrum in which the instrumental distortion has been replaced by the ideal lineshape.
Actually, the concept of deconvolution is very simple: If S(f) is the experimental frequency domain spectrum, it can be decomposed into two main components, the ideal spectrum I(f) and the instrumental distortion D(f). In other words, the observed experimental spectrum is the result of a contaminated ideal spectrum. Mathematically, this contamination is expressed with the concept of convolution which is represented by the symbol *


S(f) = I(f) * D(f)


The goal here is to find the function D(f) so that the ideal spectrum I(f) can be recovered:



I(f) = S(f) [*]-1 D(f)


Where [*]-1 denotes deconvolution. That is to say that the ideal spectrum can be recovered by means of a deconvolution which consists basically in reversing the effects of the convolution.


In practice, this process is more efficiently done in the time domain:


I(t) = S(t) / D(t)


This is possibly by considering the Convolution Theorem which states that point-wise multiplication in one domain (i.e. time domain) is equivalent to convolution in the other FT domain (e.g. frequency domain).


The complete process of Reference Deconvolution will be illustrated with an example using Mnova NMR and one 300MHz 1H-NMR spectrum in deuterated acetone (kindly provided by Gareth Morris) in which the homogeneity of the static field was deliberately perturbed. The spectrum corresponds to ODBC (ortho-dichlorobenzene) and has been folded several times in order to optimize digitization:






First, after issuing command Process/Reference deconvolution in Mnova, the User needs to select a well resolved reference signal in the frequency domain spectrum. In order to avoid numerical instabilities this signal should be a singlet. The reason is that if the reference signal has some multiplicity (i.e. a doublet), inverse FT of this reference signal (remember that Reference Deconvolution takes place in the time domain) might result in an FID with zeroes at regular intervals. As this time domain signal will be used in the denominator of the reference deconvolution function, this would result in severe discontinuities.


In this particular example, as in many others, a convenient reference could be the TMS signal (0 ppm), which in principle should be a singlet (disregarding the 13C and 29Si satellites, more about this in a moment), but as it can be noticed, it shows lineshape errors and spinning sidebands due to a combination of poor shimming and sample spinning.
In the figure below, the result of selecting the reference signal with Mnova is depicted.







In the Reference Deconvolution dialog box, there are two check boxes, 29Si and 13C satellites and the explanation is this: The TMS reference signal comes with the presence of small 29Si and 13C satellites flanking the central peak at 3.3 and 59 Hz, respectively. Since the reference line is supposed to be representative of all signals in the spectrum, any fine structure which is unique to the reference should generally be removed before deconvolution is performed. The 13C satellites are not usually a concern as they are quite distant from the main TMS peak, but the 29Si satellites are more problematic owing to their close proximity to the central signal.
So when, for example, the 29Si satellites option is ticket, the software will automatically synthesize the peaks corresponding to the 29Si satellites which in turn will be part of the reference FID model.







Once the reference region is selected, the software calculates a reference FID, Sr(t), by zeroing all the spectrum except the selected region followed by an inverse FT (there are some additional processing steps required to avoid the negative effects of the long tails of the imaginary components – The interested reader is referred to the [2-3] for further details).


Having calculated the reference FID, Sr(t), an ideal reference FID Si(t) can be computed by simply simulating an FID using the set of frequencies and amplitudes for the parent signal (e.g. TMS) and any attendant satellites, and decay rate which depends on the target lineshape, a value that can be specified in the reference deconvolution dialog box (in this example a value of 0.35 Hz has been used). With these two reference FIDs, it is possible to calculate the correction function c(t) by simply taking the (complex) ratio of both:


c(t) = Si(t) / Sr(t)


This correction function c(t) simply represents the (inverse of the) instrumental function responsible of the lineshape distortion. Multiplying the original experimental FID, Se(t), by this function yields a corrected FID, Sc(t) which may then Fourier Transformed to yield a corrected spectrum Fc having the lineshape of the specified ideal lineshape:


Sc(t) = c(t) x Se(t)
Fc = FT[Sc(t)]


Overall, the result of applying reference deconvolution to the ODBC spectrum used in this example is shown in the figure below.







Conclusion

Reference deconvolution is a powerful processing method to remove some distortions that affect all the peaks in a spectrum in the same way. In practice, this is done by extracting the distorted component from a reference signal and deconvolving the whole imperfect spectrum.
Present implementation of the algorithm in Mnova 9.0 supports 1D spectra only thus far, but extensions to 2D spectra are planned.


Biblography

[1] Why are spectral lines Lorentzian? http://www.ebyte.it/stan/blog10b.html#10May16


[2] Metz, K. R., Lam, M. M., & Webb, A. G. (2000). Reference deconvolution: A simple and effective method for resolution enhancement in nuclear magnetic resonance spectroscopy. Concepts in Magnetic Resonance, 12(1), 21–42. (link)


[3] Morris, G. A., Barjat, H., & Home, T. J. (1997). Reference deconvolution methods. Progress in Nuclear Magnetic Resonance Spectroscopy, 31(2), 197–257.




More...

Source: NMR-analysis blog
Reply With Quote


Did you find this post helpful? Yes | No

Reply
Similar Threads
Thread Thread Starter Forum Replies Last Post
Scientist III-NMR, Reference Standards Laboratory | US Pharmacopeia
Scientist III-NMR, Reference Standards Laboratory | US Pharmacopeia US - MD - Rockville, Minimum Requirements: Ph.D. in chemistry or related field OR M.S. in chemistry or related field AND 5 years of relevant laboratory experience. BRUKER experience a PLUS. Knowledge, Skills, Abilitie More...
nmrlearner Job marketplace 0 08-09-2012 08:10 AM
[NMR paper] Elimination of 13Calpha splitting in protein NMR spectra by deconvolution with maximu
Elimination of 13Calpha splitting in protein NMR spectra by deconvolution with maximum entropy reconstruction. Related Articles Elimination of 13Calpha splitting in protein NMR spectra by deconvolution with maximum entropy reconstruction. J Am Chem Soc. 2003 Mar 5;125(9):2382-3 Authors: Shimba N, Stern AS, Craik CS, Hoch JC, Dötsch V Homonuclear 13C-13C couplings can significantly reduce the sensitivity and resolution of multidimensional NMR experiments. The most important of these couplings is the 13Calpha-13Cbeta coupling, and several...
nmrlearner Journal club 0 11-24-2010 09:01 PM
Conformations of stevastelin C3 analogs: computational deconvolution of NMR data reve
Conformations of stevastelin C3 analogs: computational deconvolution of NMR data reveals conformational heterogeneity and novel motifs. Related Articles Conformations of stevastelin C3 analogs: computational deconvolution of NMR data reveals conformational heterogeneity and novel motifs. Biopolymers. 2010 Nov;93(11):968-76 Authors: Jogalekar AS The stevastelins are depsipeptide natural products that show valuable immunomodulatory and phosphatase inhibitory activity. A previous report described the synthesis, conformational analysis, and biological...
nmrlearner Journal club 0 11-06-2010 11:02 AM
[Question from NMRWiki Q&A forum] How to reference indirect dimensions in Topspin?
How to reference indirect dimensions in Topspin? Hello, When I process 3d data on topspin with the ft3d command, I get the three-dimensional spectra, but the chemical shifts of the the indirect third dimension are not right (even though the sr parameter is correct) We do 3d ncc experiments on proteins (solid-state NMR). Anyone does 3d data processing on tospin and it actually works? Check if somebody has answered this question on NMRWiki QA forum
nmrlearner News from other NMR forums 0 08-22-2010 02:30 AM
[NMR analysis blog] Fighting against peak overlap – Introducing Global Spectral Deconvolution (GSD)
Fighting against peak overlap – Introducing Global Spectral Deconvolution (GSD) 1H NMR is for sure the most powerful technique for structure elucidation, especially for small organic molecules. Typically, an organic chemist uses the chemical shift, coupling constants and integration information contained in an 1H-NMR spectrum to either verify or elucidate an unknown compound. Of course, it’s quite common that a simple 1H-NMR spectrum is not enough to unambiguously confirm a structure and thus other NMR experiments (e.g. 13C-NMR, HSQC, COSY, etc) are used to get more structural information....
nmrlearner News from NMR blogs 0 08-21-2010 09:12 PM



Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are Off
Refbacks are Off



BioNMR advertisements to pay for website hosting and domain registration. Nobody does it for us.



Powered by vBulletin® Version 3.7.3
Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.
Copyright, BioNMR.com, 2003-2013
Search Engine Friendly URLs by vBSEO 3.6.0

All times are GMT. The time now is 05:21 PM.


Map