The second round of Critical Assessment of Automated Structure Determination of Proteins by NMR: CASD-NMR-2013

The second round of the community-wide initiative Critical Assessment of automated Structure Determination of Proteins by NMR (CASD-NMR-2013) comprised ten blind target datasets, consisting of unprocessed spectral data, assigned chemical shift lists and unassigned NOESY peak and RDC lists, that were made available in both curated (i.e. manually refined) or un-curated (i.e. automatically generated) form. Ten structure calculation programs, using fully automated protocols only, generated a total of 164 three-dimensional structures (entries) for the ten targets, sometimes using both curated and un-curated lists to generate multiple entries for a single target. The accuracy of the entries could be established by comparing them to the corresponding manually solved structure of each target, which was not available at the time the data were provided. Across the entire data set, 71 % of all entries submitted achieved an accuracy relative to the reference NMR structure better than 1.5 Å. Methods based on NOESY peak lists achieved even better results with up to 100 % of the entries within the 1.5 Å threshold for some programs. However, some methods did not converge for some targets using un-curated NOESY peak lists. Over 90 % of the entries achieved an accuracy better than the more relaxed threshold of 2.5 Å that was used in the previous CASD-NMR-2010 round. Comparisons between entries generated with un-curated versus curated peaks show only marginal improvements for the latter in those cases where both calculations converged. Electronic supplementary material The online version of this article (doi:10.1007/s10858-015-9953-4) contains supplementary material, which is available to authorized users.


Introduction
Manual determination of a protein structure from nuclear magnetic resonance (NMR) spectroscopy data is a timeconsuming process, taking weeks up to several months for unfavorable cases. Moreover, as several of the required steps in this process can be complicated, the researcher must possess a high level of skill in order to produce a high quality structure. First the chemical shifts (CSs) observed in multidimensional NMR experiments must be assigned specifically to the originating atom in the molecule, a process called resonance assignment. Next, the thousands of through-space dipolar coupling signals (nuclear Overhauser effects or NOEs) observed in multidimensional NOE spectroscopy (NOESY) experiments must be identified (peak picking). These unassigned NOE peaks are then assigned to interproton interactions in the molecule, and converted into interatomic distance restraints, a process commonly called NOESY assignment. Additional conformational restraints can be obtained also from residual dipolar couplings (RDCs), scalar couplings and/or CS data. Finally, specialized software uses these distance and/or conformational restraints to produce a set of protein conformations, denoted as a structure ensemble of conformers, that is consistent with the experimental data. In the last decade, considerable effort has been expended in automation of these processes (Herrmann et al. 2002a;Linge et al. 2003a;Huang et al. 2006;Donald and Martin 2009;Williamson and Craven 2009;Gossert et al. 2011;Guerry and Herrmann 2011;Güntert and Buchner 2015), with some of these only requiring CSs (Cavalli et al. 2007;Shen et al. 2008) and several strategies are now consistently producing results comparable to a skilled researcher.
In order to evaluate the ability of automated methods to produce protein structures that closely match structures manually determined by experts, the community-wide initiative Critical Assessment of Automated Structure Determination of Proteins by NMR (CASD-NMR) was launched in 2009 (Rosato et al. 2009). In the first round, seven research groups involved in the development of NMR structure determination tools were provided with ten ''blind'' data sets consisting of assigned CS lists and unassigned NOESY peak lists that had been manually curated (i.e. filtered to remove noise peaks) and asked to generate structures using fully automated protocols only, i.e. no manual intervention other than data preparation was allowed. The same data were also used to manually produce reference structures that were only revealed after the submission of the automatically generated entries. The results were encouraging, demonstrating that automated structure-determination methods are feasible and reliable, particularly when they use manually-curated NOESY peak lists ). Subsequently, we will refer to this effort as CASD-NMR-2010.
Here, we present the input data and summarize the results from a second round of CASD-NMR that has recently been completed with the aims of assessing the progress that has been made in the automation of structuredetermination methods and investigating the need for the curation of the NOESY peak lists for accurate structuredetermination. In this round, participants were again given ten ''blind'' data sets. The research groups were also given assigned CS lists and (for some targets) RDC data correlated with HSQC peak positions, and were asked to generate structures from expert curated unassigned NOESY peak lists (as in 2009), un-curated unassigned NOESY peak lists, or unprocessed spectral data, to produce structures for comparison with the reference structures manually determined from the same data. We will refer to this effort as CASD-NMR-2013 to distinguish from the earlier CASD-NMR-2010 effort.
For the discussion in this and the accompanying paper on the structure validation of the results , we first introduce the following definitions: a target comprises the initial data set(s), the restraints generated from these data and the reference ensemble of manually generated conformers describing the three-dimensional structure of the protein of interest; entries denote individual solutions for a target automatically calculated by a specific program using a specific (sub)set of the available data, and comprise an ensemble of conformers and all restraints generated during the calculation of the ensemble.

Description of the CASD-NMR-2013 targets
As the current effort was aimed at measuring progress since the CASD-NMR-2010 round and aimed to explore the effects of curated versus un-curated NOESY peak lists, the targets were again selected to comprise single domain, relatively small monomeric proteins, for which it is expected that high-resolution NMR spectroscopy can yield high-quality structures when using state-of-the-art technology. All targets were generated by the Northeast Structural Genomics (NESG) consortium of the National Institutes of Health Protein Structure Initiative (Huang et al. 2005a; Fig. 1 and Supplementary Table S1).
Determination of the experimental structures of reference targets targets include six human protein domains: HR6470A, the homeobox domain of homeobox protein Nkx-3.1; HR6430A, the RRM domain of RNA-binding protein FUS; HR5460A, the N-terminal domain of mitotic checkpoint serine/threonine protein kinase BUB1, HR2876B and HR2876C, the N-and C-terminal domains of iron-sulfur cluster scaffold homolog NFU1; and HR8254A, the C terminal domain of DnaJ homolog subfamily C member 2. Besides targets from Homo sapiens, YR313A, the N-terminal NFU1 domain from yeast; StT322, a putative cytoplasmic protein ydiE; and two de novo designed ideal fold proteins (Koga et al. 2012), IF3 like fold design OR135 and Ploop2x3 fold design OR36, were used for the CASD-NMR-2013 study. U-15 N, 5 % 13 C-enriched proteins and U-15 N, U-13 C-enriched proteins were expressed using MJ9 minimal media (Jansson et al. 1996). The U-15 N, 5 % 13 Clabeled proteins were generated for stereo-specific assignments of isopropyl methyl groups of valines and leucines (Neri et al. 1989) and for residual dipolar coupling (RDC) measurements (Prestegard et al. 2004). The final purified protein samples prepared at Rutgers include a short N-terminal tag, with sequence MGH 6 SHM, for targets HR6470A, HR6430A, HR5460A, HR2876B, HR2876C and YR313A, or a short C-terminal LEHHHHHH tag for targets OR135 and OR36. The NMR samples of the two targets StT322 and HR8254 were prepared in Toronto and initially produced with larger purification tags, which were subsequently proteolytically removed to obtain the final NMR samples. The purified proteins were dissolved in 90 % 1 H 2 O/10 % 2 H 2 O NMR buffers optimized for the individual proteins.

NMR data collection
All NMR spectra were recorded at 25°C using cryogenic NMR probes at Rutgers University or University of Toronto. Triple resonance NMR data were collected on a Varian INOVA 600 MHz spectrometer or Bruker AVANCE 800 MHz spectrometers, while simultaneous 3D 15 N/ 13 C aliphatic / 13 C aromatic -edited NOESY (mixing time: 100 ms) and 3D 13 C-edited aromatic NOESY (mixing time: 100 ms) spectra were acquired on a Bruker AVANCE 800 MHz spectrometer. 2D constant-time 1 H-13 C HSQC spectra, with 28 and 42 ms constant-time delays, were recorded for a U-15 N, 5 % 13 C-enriched samples on the Varian INOVA 600 MHz spectrometer in order to obtain stereo-specific assignments for isopropyl groups of valines and leucines (Neri et al. 1989). Data recorded on targets StT322 and HR8254 used non-uniform sampling methods in the indirect dimensions (Orekhov et al. 2003;Gutmanas Fig. 1 Side by side superimposed backbone ribbon traces and cartoon representations for the ten manually-determined CASD-NMR-2013 reference structures, labeled with PDB codes and coloured blue to red from N-to C-terminus. Ill-defined regions are shown in light grey J Biomol NMR (2015) 62: 413-424 415 et al. 2002). All NMR data were processed using the program NMRPipe (Delaglio et al. 1995) and analysed using the program XEASY (Bartels et al. 1995). Spectra were referenced to external DSS. Sequence-specific resonance assignments were determined as described previously (Baran et al. 2004;Moseley et al. 2001). Targets StT322 and HR8254 employed the ABACUS semi-automated assignment strategy (Lemak et al. 2010). Chemical shift data were deposited in the Biological Magnetic Resonance Bank (Ulrich et al. 2008; cf. Table 1). Residual dipolar coupling measurements were made at University of Georgia on a 600 MHz Varian INOVA spectrometer. Samples were aligned with polyacrylamide gel or polyethylene-glycolalkyl bicelles and RDCs were collected using either interleaved HSQC-TROSY or J-modulation sequences as described elsewhere Lee et al. 2010).

Manual NMR structure calculations
For NMR structure calculations, initial NOESY peak lists containing expected intra-residue, sequential, and a-helical medium-range NOE peaks were first generated from the resonance assignments and then manually edited by visual inspection of the NOESY spectra. Subsequent manual peak picking was then used to identify remaining peaks that were not easily identified from these simulated spectra, which arise primarily from long-range NOEs. Backbone dihedral angle constraints were derived from chemical shift data using the program TALOS? (Shen et al. 2009) for residues located in well-defined secondary structure elements. The program CYANA (Herrmann et al. 2002b) was used to automatically assign NOEs and to calculate the NMR structure ensemble. Subsequently, RPF-DP analysis ) was used to guide manual peak list editing, including iterative cycles of noise/artefact peak removal, peak picking and NOE assignments. The output of the RPF-DP program was used to iteratively refine the NOESY peaks list, followed by cycles of automated analysis and structure generation, together with RDC data, with CYANA. In the final cycle of structure generation calculations, the 20 conformers with the lowest target function value were refined with RDC data in the presence of explicit water solvent (Linge et al. 2003b) using the program CNS (Brunger et al. 1998). The resulting structure coordinates and restraints were deposited in the Protein Data Bank (Berman et al. 2003), and relevant metrics are summarized in Supplementary Table S2.

Validation of the experimental NMR reference structures
NMR data statistics, structural statistics, and global structure quality factors including Verify3D (Lüthy et al. 1992), ProsaII (Sippl 1993), PROCHECK (Laskowski et al. 1993), and MolProbity (Chen et al. 2009) raw and statistical Z-scores were computed using the PSVS version 1.4 software package (Bhattacharya et al. 2006). The global goodness-of-fit of the final structure ensembles with the NOESY peak list data was determined using the RPF analysis program (Huang et al. 2005b. The NESG minimum standard scores for each global structure-quality score were used as an initial guide to assess the quality of each structure. In addition, a closer examination of the local structure quality was performed to identify potential problem areas. If either potential global or local structural problems were identified, then this information was reported back to the researcher performing the structure determination, and the researcher was asked to carefully reexamine the data and to resolve any problematic issues. The structural ensembles of the ten targets have excellent validation statistics, both with respect to structural criteria and model versus data. Superpositions of the NMR ensembles of the reference structures are shown in Fig. 1 and the structural statistics are presented in Supplementary Information (Table S2). A detailed analysis of the structural quality of the reference protein structures is also given in the accompanying paper .

Production of NOESY peak lists
Two sets of NOESY peak list data were provided for each protein target in CASD-NMR-2013: (1) initial automatically-picked un-curated (or raw) NOESY peak lists and (2) the final manually-curated NOESY peak lists used for generating the reference structures that are deposited in the PDB. These NOESY peak lists and the final chemical shifts were provided in both XEASY (CYANA) and Sparky formats. The following protocol was used to prepare these ''un-curated raw'' NOESY peaks lists. First, peak lists for 2D 15 N-HSQC, aliphatic 13 C-HSQC and aromatic 13 C-HSQC spectra were simulated, without spectral folding, from the complete NMR resonance assignments. These resonance frequencies were then manually curated by comparison with 3D 13 C, 15 N-edited NOESY spectra using the interactive spectral visualization program XEASY. This process ensures good matching between the chemical shifts in these simulated HSQC spectra and the experimental NOESY peak list. These refined simulated HSQC peak lists were then used to guide automatic peak picking of the entire experimental NOESY spectrum. The ''restricted peak picking'' module of the program Sparky (Goddard and Keller, Sparky 3, University of California, San Francisco) was used to automatically peak pick the 3D 13 C, 15 N-edited NOESY spectra. The 15 N, 13 C and direct 1 H tolerances were chosen based on data resolution; typically 0.3-0.4 ppm for 15 N and 13 C dimensions, and Well-defined ranges as determined by CyRange (Kirchner and Guentert 2011) 5 Closest homologue in the PDB dated prior to the release date of the target. Where several entries had the same homology we have given priority to structure ensembles for systems and conditions as similar to the target as possible 0.02-0.03 ppm for the direct 1 H dimension. For the indirect 1 H dimension, the match tolerance was set to the spectral width, as all frequencies are possible. Additional details on this simple peak picking protocol are available at http://wiki.nesg.org.

CASD-NMR-2013 data sets
The data sets for CASD-NMR-2013 comprised CS assignments in both version 2.1 and 3.1 of the NMRSTAR format (Markley et al. 2003) and unassigned NOESY peak lists in SPARKY and/or XEASY/CARA format. The data and their metadata were made available on the WeNMR (Wassenaar et al. 2012) CASD-NMR website (http://www. wenmr.eu/wenmr/casd-nmr). For all targets, raw NOESY spectral data were also made available. CS assignments, unrefined peak lists and raw data were released simultaneously to the participants, 6-8 weeks ahead of the release of the target structure from the PDB. Refined lists were released 4 weeks after the release of the initial unrefined data, so that the participants had a minimum of 2 weeks to use them for structure calculation. With this design, a period of 4 weeks was available to generate an entry with unrefined peak lists or raw data, and two to four additional weeks were available to generate an entry with refined lists, while the manually solved structure was on hold in the PDB. The entries were deposited directly by the participants into a password-protected database, again via the CASD-NMR website. An overview of all targets and the data made available to the CASD-NMR-2013 round is given in Table 1.

Description of the CASD-NMR-2013 entries
Current structure generation protocols are based either on NOESY-based distance restraints, on CS data, or a combination of the two. In CASD-NMR-2013, the data were used by a total of twelve different automated protocols, subsequently referred to as 'Programs' for NMR structure generation. The CASD-NMR-2013 effort is a near complete representation of the array of computational methods currently available to the NMR spectroscopist for generating three-dimensional atomic-resolution structures from NMR data. Detailed descriptions of the different participating programs and of the protocols used to prepare the entries are available as published or in subsequent papers (Table 2). A total of 164 entries were submitted covering all ten targets (Fig. 2). UNIO and Ponderosa were the only programs that submitted entries calculated directly from raw, unprocessed spectral data. However, not all programs submitted entries for all targets and not all the programs that used NOESY data submitted entries for both un-curated and curated peak lists. This could potentially distort the apparent success rates of the programs. Queries with the submitting groups revealed that 36 calculations that might otherwise have been expected had never been attempted (or, in a couple of cases, could not be carried out in time because of temporary file I/O problems). Eleven calculations were attempted but did not converge: Five ARIA and one Cheshire-YAPP entries for un-curated peaks were deposited but qualified as 'incorrect' by the submitters, and these served as controls, but were excluded from all analyses. A further three CYANA calculations started from unrefined peaks lists and two Rosetta web server calculations did not converge and were never deposited. Finally, nine entries from three different programs were missing without it being possible to determine the reason.

Results and discussion
The results of CASD-NMR-2013 provide a comprehensive comparison of the performance of most currently available programs for automated protein structure generation from NMR data. We have evaluated both methods that primarily rely on NOESY peak lists as the input data and methods that use CS data as the primary input data. The latter are sometimes augmented by the use of the NOESY peak lists or RDCs, usually for the purpose of selection from the initially generated ensemble of structures. We also assessed the impact of using un-curated rather than curated NOESY peak lists. Furthermore, we have evaluated structures generated by methods that use CS data only (Cheshire and Rosetta) or raw spectral data (Ponderosa and UNIO) as the input data.  The wwPDB NMR Validation Task Force has recommended that the data for ill-defined regions should be part of the deposited dataset, but should not be included in the validation analysis (Montelione et al. 2013). It also recommended the program CyRange (Kirchner and Guentert 2011) as a robust tool for determining the well-defined and ill-defined regions. Following this recommendation, we established the well-defined regions for the ten targets (Table 1) and used these to superpose targets and entries. Assuming that the manually determined reference ensemble constitutes the correct representation of the three-dimensional structure of the target protein, the average pairwise backbone RMSD between the target ensemble and the entry ensemble were used as a measure of accuracy (RMSD bias). The use of CYANA to calculate the reference structures could in theory bias the results, but no such effect was evident from the data. A closer investigation would require recalculating structures with different programs using identical, manually curated input data. An entry is considered to be indistinguishable from the target when the average RMSD between the two ensembles is less than the sum of their ensemble convergence. We established a threshold of 1.5 Å to be a reasonable criterion, above which any ensembles can be regarded as describing structures with differences beyond experimental uncertainty (see Ragan et al. 2015). This threshold is somewhat more restrictive than that used to identify accurate structures in CASD-NMR-2010. We also used a slightly different method to compute the RMSD bias, i.e. as the average RMSD between the conformers in the two ensembles instead of the RMSD between the two average conformers, as was done for CASD-NMR-2010. Thus, a level of accuracy comparable to CASD-NMR-2010 is in between 2 and 2.5 Å . Figure 3 shows the fraction of entries provided by each method using un-curated peak lists or raw data in relation to accuracy (cf. Table 3). Because of the similar outcomes, we grouped entries calculated with or without RDC restraints. The median accuracy over the entire dataset is 1.14 Å , with 71 % of the entries below the 1.5 Å threshold. Approaches using NOESY peak lists (curated or un-curated) achieved the best results, with a median accuracy of 1.07 Å and 80 % of the entries below the threshold (increasing up to 100 % for some programs). In contrast, calculations based on either raw spectral data or CS-only data performed less well, both approaches yielding a median accuracy of 1.5 Å with 50 % of ensembles below the threshold. These statistics are discussed in more detail in the accompanying paper .
CASD-NMR-2013 had a specific focus on the structure calculation performance achievable starting from automatically generated un-curated peak lists or raw spectral data, in addition to evaluating the progress of CS-only approaches. For the former set of calculations, in most cases the fraction of accurate entries (RMSD bias B1.5 Å ) was 50 % or more of all submitted entries, up to 100 % for the NOESY-based methods (ARIA, CYANA UNIO; Table 3). Note however that some of these NOESY-based methods did not submit entries for all targets, precluding head-to-head comparisons of the methods for more than a few targets (Table 3). In absolute terms, NOESY-based methods provided accurate entries for up to seven of the ten targets using un-curated lists (ASDP; UNIO). This value increases to ten (ASDP) and a success rate of 100 % Fig. 2 Targets submitted per program. A target is counted as submitted if there is at least one entry; programs often submitted multiple entries for a single target. Calculations that did not converge are ignored. Colour (see legend) encodes the targets calculated using different input data sets; e.g. using curated NOESY peak lists only (light blue) or two entries for one target with one using curated and one using curated peaklists (dark blue). Results for ASDP were provided using two different refinement methods (CNS and Rosetta), but at least one submission was provided for each of the 10 targets J Biomol NMR (2015) 62:413-424 419 assuming a more relaxed accuracy threshold of *2.5 Å , comparable to the previous CASD-NMR-2010 round . For approaches based on raw data, the performance is close to that observed for traditional methods, with UNIO performing better than Ponderosa. Chemical-shift based tools can be run using CS-only data, or in conjunction with additional un-curated NOESY or RDC peak lists to filter or bias the results obtained with CS data. Cheshire-YAPP and Autonoe-Rosetta, which integrate NOESY data, provided four and six correct structures, respectively, corresponding to 50 and 60 % of the submissions (Fig. 3a). The number of correct structures was lower for CS-only based calculations, ranging between three and four for the various tools. CS-HM-Rosetta, which exploits RDC's and restraints derived from protein homology analysis, provided four accurate structures; all of its eight submissions were within 2.5 Å from the target structure. I-Tasser joined CASD-NMR at an advanced stage and thus could provide blind entries only for the two latest targets; both entries were accurate. BE-Metadynamics similarly participated only for the last target, yielding a structure with an accuracy of 1.68 Å . While it is difficult to compare these methods with other methods assessed because only a few targets were submitted, their results were included in our analysis to encourage the participation of a wide range of methods in the CASD-NMR experiment. Figure 3b shows a comparison of accuracies of the entries obtained from either the use of un-curated NOESY peak lists or curated peak lists. It reveals that, in cases where results were submitted for both curated and un-curated peak lists, there is no advantage to using curated peak lists as the two distributions are overlapping with very similar median accuracies of 1.08 Å (79 % of entries below the threshold) and 1.05 Å (80 % of the entries below the threshold), respectively. It appears that the iterative procedures implemented in the programs are efficient at filtering the peaks for consistent information. It should be noted however, that five out of ten ARIA calculations and three out of eight CYANA calculations with un-curated peak lists failed to converge, where the corresponding calculations with curated peak lists converged to good Fig. 3 Accuracy of CASD-NMR-2013 entries based on un-curated information. a Number of CASD-NMR2013 entries based on uncurated NOESY, raw spectral data, or chemical shift (CS)-only data, colour-coded based on their accuracy. For each program, we include all targets for which there was at least one entry submitted. Each column is colour-coded based on the average accuracy of the submitted entry(ies) (compare to Table 2). Green high accuracy (RMSD bias to the reference \1.5 Å ); orange intermediate accuracy (RMSD bias to the reference\2.5 Å ); red low accuracy. b Histograms of the accuracy of entries using curated (blue, n = 55) or un-curated (orange, n = 63) peak lists. The RMSD values were calculated as the average of the pairwise root mean square deviation of the backbone atoms between the conformers in the target and entry ensembles using the well-defined regions as determined by CyRange (cf. Table 2) quality structures. The convergence rate seems to depend on the quality of the un-curated peak lists. We estimated the quality of the un-curated peak list as the ratio r of the number of curated to un-curated peaks (Table 1). Three targets with r [ 0.8 converged for both the CYANA and ARIA programs, and two targets with r \ 0.2 failed to converge for both programs. The remaining five targets had 0.5 \ r \ 0.7; of these three failed to converge for ARIA, one failed to converge for CYANA, and two were not attempted by CYANA for unrelated reasons. All ten targets were submitted by ASDP, using either CNS or Rosetta refinement. We also analysed the relation between the number of un-curated peaks and structure accuracy. For the entries that do converge, the protein length and the number of assigned chemical shifts are only correlated with the final structure accuracy for the methods starting from raw data ( Table S3), suggesting that these approaches are still sensitive to protein size. The number of unrefined peaks Structure precision is also shown in parentheses. When multiple structures were submitted for a given target based on the same method, the average accuracy and precision are given. Structures marked as not converged or incorrect at the time of submission were excluded. At the bottom of the table, we report the median accuracy of each tool. In addition, the number of structures with an accuracy better than a fixed threshold are reported, together with the percentage of entries that were better than (i.e. correct within) the threshold. Tools were grouped according to the main input data used J Biomol NMR (2015) 62:413-424 421 from the 13 C and 15 N NOESY, however, display a weak correlation with the final structure accuracy (Supplementary Figure S1 and Table S3). This is an additional indication that the combination of the quality and complexity of the NOESY spectra, which determines the number of un-curated NOESY peaks, may be relevant for the final structure accuracy. Two of the targets were particularly challenging, StT322 and HR8254A. StT322 includes a large ill-defined region (Fig. 1), and HR8254A has a long third-helix extending outside the core, which makes RMSD metrics very sensitive to the precise determination of helical tilt angles (Fig. 1, see also Fig. 3 of Huang et al. 2015). These two data sets were the only ones to include non-uniform sampled NMR data, had the lowest ratio of curated to uncurated peaks (r \ 0.2) and did not have RDC data (Table 1). Targets StT322 and HR8254A failed to converge in both CYANA and ARIA for un-curated peaks. Indeed, only one of the methods based on un-curated NOESY peaks lists, ASDP-Rosetta, reported results for both StT322 and HR8254A (Table 3); high RMSDs to these reference structures for the ASDP entries (2.56 and 1.73 Å ) were the highest across all the methods based on un-curated NOESY peak lists (Table 3), significantly increasing the median RMSD score for ASDP results. Pondorosa, using raw spectra, also submitted results for StT322 and HR8254A with significantly higher RMSDs (3.69 and 3.01 Å , respectively) compared with most of the other Ponderosa results (Table 3). Entries submitted for the StT322 and HR8254A also had the worst accuracy scores among the targets submitted for Autonoe-Rosetta, while Cheshire did reasonably well with un-curated peaks provided for target StT322 (1.30 Å ). As only a subset of the methods submitted results for these particularly challenging targets StT322 and HR8254A care should be taken in interpreting median RMSDs and percentages of structures in various accuracy ranges, summarized in Table 3, in comparing the performance of the several methods of CASD-NMR-2013.
In conclusion, the results from CASD-NMR-2013 demonstrate that for small, single domain proteins automated structure determination protocols are capable of reliably producing structures of comparable accuracy to those generated by a skilled researcher. Compared to CASD-NMR-2010, a significantly larger number of entries were successful and the overall accuracy was significantly better. As in the previous round, the performance of NOESY-based methods was superior to that of CS-only methods. Augmenting the latter with NOESY or other information resulted in an appreciable improvement. The present data show that the use of curated peak lists only yields marginal improvements over un-curated, automatically picked lists, provided that the calculations converge.
High quality peak lists, with a ratio of curated to un-curated peaks above 0.5, greatly aids convergence, indicating that by using a conservative peak-picking approach the manual filtering step is no longer required for the successful application of the automated protocols, potentially broadening their utility. A more detailed comparison of the performance of individual programs can be found in the accompanying paper on validation of the structures . The data collected for CASD-NMR-2013 will continue to be available as they provide a valuable resource for methods development in NMR Lange 2014;Buchner and Güntert 2015). Subsequent rounds of CASD-NMR are expected to focus on automated structure determination from more challenging datasets e.g. for larger and/or multi-domain proteins, datasets containing various assignment errors or datasets lacking CS assignments. Another topic that merits further investigation is the incorporation of robust, fully automatic peak picking, as done in the present set of calculations by UNIO and PONDEROSA. We will also upgrade the data submission process for CASD-NMR and center it on the new NMR Exchange Format (NEF) (Gutmanas et al. 2015), which should facilitate the participation of a broader range of computational groups in future CASD-NMR experiments.