A tracked approach for automated NMR assignments in proteins (TATAPRO)

A novel automated approach for the sequence specific NMR assignments of 1HN, 13Calpha, 13Cbeta, 13C'/1Halpha and 15N spins in proteins, using triple resonance experimental data, is presented. The algorithm, TATAPRO (Tracked AuTomated Assignments in Proteins) utilizes the protein primary sequence and peak lists from a set of triple resonance spectra which correlate 1HN and 15N chemical shifts with those of 13Calpha, 13Cbeta and 13C'/1Halpha. The information derived from such correlations is used to create a 'master-_list' consisting of all possible sets of 1HN(i), 15N(i)13Calpha(i),13Cbeta(i) 13C'beta(i)/1Halpha(i), 13Calpha(i-1), 13Cbeta(i-1) and 13C'(i-1)/1Halpha(i-1) chemical shifts. On the basis of an extensive statistical analysis of 13Calpha and 13Cbeta chemical shift data of proteins derived from the BioMagResBank (BMRB), it is shown that the 20 amino acid residues can be grouped into eight distinct categories, each of which is assigned a unique two-digit code. Such a code is used to tag individual sets of chemical shifts in the master_list and also to translate the protein primary sequence into an array called pps_array. The program then uses the master_list to search for neighbouring partners of a given amino acid residue along the polypeptide chain and sequentially assigns a maximum possible stretch of residues on either side. While doing so. each assigned residue is tracked in an array called assig_array, with the two-digit code assigned earlier. The assig_array is then mapped onto the pps_array for sequence specific resonance assignment. The program has been tested using experimental data on a calcium binding protein from Entamoeba histolytica (Eh-CaBP, 15 kDa) having substantial internal sequence homology and using published data on four other proteins in the molecular weight range of 18-42 kDa. In all the cases, nearly complete sequence specific resonance assignments (> 95%) are obtained. Furthermore, the reliability of the program has been tested by deleting sets of chemical shifts randomly from the master_list created for the test proteins.