Automated resonance assignment with FLYA (EMBO 2013)

We want to determine the resonance assignments and the structure of a protein.

The protein sequence is stored in the file 'demo.seq'.

The following spectra have been measured:

• [¹H,¹³C]-HSQC, called 'C13H1' in FLYA
• [¹H,¹⁵N]-HSQC, called 'N15H1' in FLYA
• HNCA
• HN(CO)CA, called 'HN_CO_CA' in FLYA
• HNCO
• HN(CA)CO, called 'HN_CA_CO' in FLYA
• CBCANH
• CBCACONH
• HBHACONH
• HCCH-TOCSY, called 'HCCH24' in FLYA
• HCCH-COSY, called 'HCCH7' in FLYA
• C(CO)NH, called 'C_CO_NH' in FLYA
• HC(CO)NH, called 'HC_CO_NH' in FLYA

Peak lists that have been prepared by automatic peak picking with the program NMRView are stored in files XXX.peaks, where XXX denotes the FLYA spectrum type.

Each peak list starts with a header that defines the experiment type and the order of dimensions. For instance, for HNCA.peaks:

# Number of dimensions 3
#FORMAT xeasy3D
#INAME 1 HN
#INAME 2 C
#INAME 3 N
#SPECTRUM HNCA  HN C N
      5   6.475  58.033  98.548 1 U   2.769E+02  0.000E+00 e 0     0     0     0
      6   6.476  62.123  98.126 1 U   2.571E+01  0.000E+00 e 0     0     0     0
      7   6.475  54.017  98.159 1 U   2.547E+01  0.000E+00 e 0     0     0     0

The first line specifies the number of dimensions (3 in this case). The '#SPECTRUM' lines gives the experiment type (HNCA), followed by an identifier for each dimension of the peak list (HN C N) that specifies which chemical shift is stored in the corresponding dimension of the peak list. Below follows one for every peak. For example, the first peak in the HNCA.peaks list has

• Peak number 5
• HN chemical shift 6.475
• C (CA) chemical shift 58.033
• N chemical shift 98.548

The other data are irrelevant for automated chemical shift assignment with FLYA. In particular, the peak volume or intensity (2.769E+02) is not used by the algorithm.