Resonance Assignment/CARA: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(10 intermediate revisions by 2 users not shown)
Line 1: Line 1:
== Loading a new template  ==
The topics below provide an overview of resonance assignment of 15N,13C-labeled proteins using CARA. For more information and tutorials please consult [http://www.nmr.ch CARA web site].


To start a new structure determination project in CARA you need to load a template. A template is a CARA repository without project that contains definitions for residue types, spectrum type and LUA scripts. You can either load the default template from the [http://www.cara.ethz.ch/Wiki/TemplatesPage templates page] of the official CARA web-site, or extract a template from an existing CARA repository.
*[[CARA Introduction|Introduction]]
*[[CARA vs Xeasy|Differences from XEASY]] 
*[[CARA Scopes|Scopes - Displaying NMR Spectra]]


In CARA click '''File -> New from template''' and select the appropriate <tt>.cara</tt> file. For detailed instructions see [http://www.cara.ethz.ch/Wiki/ImportingTemplate]


<br>[[Image:CARA New Template Dialogue.png]]  
*[[Resonance Assignment/CARA/Starting a new project|Starting a new project]]
*[[Resonance Assignment/CARA/Backbone assignment|Backbone resonance assignment]]
*Side Chain Assignment
**[[Aliphatic Side Chain Assignment with CARA|Aliphatic side-chain assignment]]
**[[Aromatic Side Chain Assignment with CARA|Aromatic side-chain assignment]]
**[[Amide Side Chain Assignment with CARA|Amide side-chain assignment]]
*[[SSAFromFractional13CLabeledSample|Stereospecific Assignment of Leu and Val Isopropyl Groups]]
*[[Creating NOESY peaklists with CARA|Creating NOESY peaklists]]


CARA Template for GFT spectra [[Media:BoR54_template.cara|BoR54_template.cara]]: Latest GFT template for BoR54 project


== Starting a new project ==
Also see [[Resonance_Assignment/CARA/Backbone_assignment_GFT|backbone assignment with GFT spectra]]
 
Normally you would create a new project based on a protein sequence in XEASY format.
 
If you have the protein sequence in one-letter format you can convert it to XEASY <tt>.seq</tt> file the following way:
#Save the one-letter sequence in a <tt>.aa</tt> file, e.g. <tt>foo.aa</tt>
#Execute '''OneLetterFileToSeqFile''' LUA script from '''Terminal''' tab of the main window. Select your <tt>.aa</tt> file as input and give full path for the output <tt>.seq</tt> file.
 
[[Image:CARA One Letter to Seqfile.png]]
 
To create a new project for you protein in the repository right-lick '''Projects''' -> '''Import Project from Sequence...'''. Select the corresponding <tt>.seq</tt> file in the file explorer and click on the '''Residues''' button in the pop-up menu.
 
[[Image:CARA proj import from seq.png]]
 
This will create a project entry with the same name as the *.seq file, with the sequence field appropriately filled and all other fields empty. The new project will have the same name as the provided <tt>.seq</tt> file, though it can be renamed.
 
[[Image:CARA new project.png]]
 
Even though a repository can contain several projects, it is a good practice to have only one project per repository. For example, many LUA scripts assume that there is a single project in a repository.
 
You can also import a project from a BMRB file. In this case chemical shift data will be imported in addition to the sequence. For more details see [http://www.cara.ethz.ch/Wiki/CreateNewProject]. Also check the page on importing XEASY files: [http://www.cara.ethz.ch/Wiki/WorkingWithOtherProgramsImportExport]
 
== Importing NMR spectra ==
 
The next step is to import NMR spectra into the project. Left-click on the '''Spectra''' node in the project, then right-click in the panel to the right. In the context menu move the mouse cursor to '''Add Spectrum''' and select the appropriate spectrum type, then chose the spectrum file in the file explorer. The following formats can be imported: XEASY (<tt>.3D.param</tt>), Sparky (<tt>.ucsf</tt>), BRUKER (<tt>.rr</tt>), NMRPipe (<tt>.ft</tt>) and Felix (<tt>.mat</tt>).
 
Most common spectra types are defined in the standard template. Please note that in the default template there are distinct spectrum type definitions for aliphatic and aromatic versions 2D [13C,1H] HSQC, 3D 13C-resolved NOESY, and HCCH spectra. You can also define additional spectrum types as described in [http://www.cara.ethz.ch/Wiki/CreateSpectrumType http://www.cara.ethz.ch/Wiki/CreateSpectrumType].
 
The order of spectral dimensions is normally irrelevant as CARA determines the correct mapping to the spectrum type from the type of nucleus associated each dimension. However, if the mapping is not unique (as in the case of 3D HBHA(CO)NH, 3D heteronuclear-resolved NOESY, 3D (H)CCH, 3D H(C)CH spectra, etc.), then a pop-up widow will prompt you to define the correct mapping.
 
[[File:CARA_map_to_spectrum_type.png]]
 
For a given dimension CARA can "guess" the corresponding nucleus type from the ppm range of the dimension. This, however, may still be insufficient to distinguish between 15N and carbonyl 13C dimensions, or between methyl 13C and 1H. For XEASY spectra one solution to this problem is to make sure before importing that <tt>.3D.param</tt> include <tt>Identifier for dimension wX</tt> lines as in this example for a 3D HNCO spectrum:
<pre>
Version ....................... 1
Number of dimensions .......... 3
16 or 8 bit file type ......... 16
Spectrometer frequency in w1 .. 60.7413
Spectrometer frequency in w2 .. 150.7300
Spectrometer frequency in w3 .. 599.4460
Spectral sweep width in w1 .... 27.6576
Spectral sweep width in w2 .... 15.9225
Spectral sweep width in w3 .... 8.3410
Maximum chemical shift in w1 .. 131.7220
Maximum chemical shift in w2 .. 185.6120
Maximum chemical shift in w3 .. 13.1120
Size of spectrum in w1 ........ 128
Size of spectrum in w2 ........ 256
Size of spectrum in w3 ........ 512
Submatrix size in w1 .......... 16
Submatrix size in w2 .......... 32
Submatrix size in w3 .......... 64
Permutation for w1 ............ 3
Permutation for w2 ............ 2
Permutation for w3 ............ 1
Folding in w1 ................. RSH
Folding in w2 ................. RSH
Folding in w3 ................. RSH
Type of spectrum .............. ?
Identifier for dimension w1 ... N
Identifier for dimension w2 ... C
Identifier for dimension w3 ... H
</pre>
 
For more details see [http://www.cara.ethz.ch/Wiki/ImportSpectra]
 
== Referencing and aliasing in NMR spectra ==
 
There could be small systematic shifts even between properly transformed spectra. An easy way to correct these offsets interactively is described in the instructions for using SynchroScope [http://www.cara.ethz.ch/Wiki/SynchroScope]. Another possibility is to right-click on the spectrum of interest in the main window and select '''Calibrate Spectrum...''' from the context menu. You can explicitly enter the required corrections (in ppm) in in the pop-up window.
 
Please note that these are corrections for systematic shifts only - you cannot fix incorrect spectral widths within CARA. These calibration corrections are internal to CARA and the original XEASY <tt>.3D.param</tt> files remain unchanged. To properly use the same spectra and CARA assignments with an external program, such as ATNOSCANDID, you have to export the CARA calibration by right-clicking on the spectrum tab and selecting '''Write Calibration...''' from the context menu, which will update the corresponding <tt>.3D.param</tt> files.
 
In contrast to XEASY, peaks and spins in CARA are always shown at their true chemical shifts and lack folding attributes. To see folded peaks you have to scroll the spectrum beyond its default range in any given display scope (the '''View''' -> '''Show Folded''' checkbox must be selected to display folded regions, otherwise they will be black).
 
Folding attributes of spectral dimensions determine how "unfolded" regions are generated. CARA interprets the dimension folding attributes <tt>Folding in wX</tt> of <tt>.3D.param</tt> files upon importing, which have to be either RSH or TPPI. RSH stands for the most common wrap-around aliasing and TPPI means mirror-reflection folding. You can also change folding attribute after a spectrum is imported. Just expand the spectrum node to see the dimensions, right-click on a dimension tab and select '''Set Folding...''' from the context menu.
 
[[Image:CARA set folding.png]]
 
For additional details and examples see [http://www.cara.ethz.ch/Wiki/FoldingAndAliasing] .
 
 
 
 
#[[CARA Introduction|Introduction]]
#[[Spin System Identification with CARA|Spin System Identification in 2D 15N-HSQC and 3D HNNCO]]
#[[Backbone Assignment with CARA|Backbone Resonance Assignment]]
#[[HA and HB Assignment with CARA|Assignment of HA and HB Resonances with (4,3)D GFT HABCAB(CO)NHN]]
#Side Chain Assignment
##[[Aliphatic Side Chain Assignment with CARA|Aliphatic side-chain assignment]]
##[[Aromatic Side Chain Assignment with CARA|Aromatic side-chain assignment]]
##[[Amide Side Chain Assignment with CARA|Amide side-chain assignment]]

Latest revision as of 20:13, 8 August 2012