ARTINA structure determination
Job ID: fcf110f1d0144d3895e7874426e18e21| Job date/time | |
| Job submitted | April 5, 2022, 3:40 p.m. |
| Job started | April 5, 2022, 4:16 p.m. |
| Job completed | April 5, 2022, 10:44 p.m. |
| Job call details | |
| App identifier | app_csa_struct_calc |
| Queue | public |
| Job status | finished |
| Job hardware details | |
| CPU credits | 420.0 |
| GPU credits | 0.0 |
Download all job results: results.zip
Automated peak picking results
Table 1. Results of automated peak picking. Peak lists with low F1 score (< 0.25) are highlighed in red. Frequently, low F1 scores indicate wrong spectrum annotations in the project summary, such as axes transposition (e.g. C-H vs. H-C) or wrong experiment type specification.
Peak picking results compared to 1000+ spectra benchmark
Figure 1. Evaluation of automated peak picking output in reference to a benchmark of 1000+ NMR spectra. Each box plot presents the distribution of the F1 score calculated for all benchmark spectra of the given type. Red dots present F1 scores calculated for the output of this peak picking application call. Relative comparison of F1 scores serves as a quality verification of the peak picking output. Scores above the median (orange line) indicate that the automated peak picking routine was able to process a spectrum more accurately than 50% of the spectra in the benchmark. Low F1 scores might be related to a low signal-to-noise ratio of the spectrum or to a wrong definition of the spectrum type in the project summary.
Data consistency check: spectra referencing
Proper spectra referencing is a major precondition for automated shift assignment and structure calculation. The table below presents results of the CYANA PeakMatch algorithm (Buchner et al., 2013), which has been executed for each pair of spectra that share at least two axes. Each row in the table contains an estimate of the relative shift between target and reference spectrum. The score presented in the last column (Equation 2, Buchner et al., 2013) represents the confidence of shift identification (the higher the better).
The relative shifts listed in the table below should be verified manually and corrected before submitting an automated shift assignment or structure calculation job. They can be applied to the target spectrum using specific tags (e.g. w1:0.025 to increase the chemical shift coordinates in the 1st dimension of the spectrum and peak list by 0.025 ppm) to define spectrum shifts in the project summary. It is not necessary to repeat automated peak picking after entering such tags in the project summary. The system will apply these corrections automatically.
| Compared peak lists | Relative shift | Score | |
| Reference spectrum | Target spectrum | ||
| N15HSQC | N15NOESY | N=0.0, HN=0.0 | 0.25 |
| HCcoNH_@ALI | N15HSQC | N=0.0, H=0.0 | 0.85 |
| CBCANH | N15HSQC | N=0.0, H=0.0 | 1.32 |
1 Non-zero shifts are significant if greater than 0.2 ppm for 13C or 15N, or 0.015 ppm for 1H.
2 Relative referencing with scores below 0.1 is uncertain.
Table 2. Results of automated spectra referencing. Rows presenting possible referencing errors are highlighted in red. The relative shift is reported for spectra having at least 25 strong cross-peaks (Table 1).
Automated shift assignment results
Figure 1.The diagram illustrates atoms in the protein sequence with color-coded chemical shift information. Dark-blue rectangles represent atoms assigned confidently by the FLYA algorithm, whereas for light-blue ones the assignment is uncertain. The row labeled HN/Hα shows for each residue HN on the left and Hα in the center. The N/Cα/C' row shows for each residue the N, Cα, and C' assignments from left to right. The rows βη show the sidechain assignments for the heavy atoms in the center and for the hydrogen atoms to the left and right. In the case of branched side-chains, the corresponding row is split into an upper part for one branch and a lower part for the other branch.
| Spectrum name | Expected | Assigned expected | Measured | Assigned measured | Ratio | Assigned peak lists |
| C13NOESY_@ALI_@FLYA | 11219 | 4117 (36.7%) | 5605 | 2661 (47.48%) | 1.5 | XEASY Sparky |
| C13NOESY_@ARO_@FLYA | 1538 | 500 (32.51%) | 617 | 216 (35.01%) | 2.3 | XEASY Sparky |
| CBCANH_@FLYA | 363 | 337 (92.84%) | 547 | 309 (56.49%) | 1.1 | XEASY Sparky |
| HCCHTOCSY_@ALI_@FLYA | 2145 | 1354 (63.12%) | 2653 | 1083 (40.82%) | 1.3 | XEASY Sparky |
| HCcoNH_@ALI_@FLYA | 477 | 348 (72.96%) | 432 | 302 (69.91%) | 1.2 | XEASY Sparky |
| N15HSQC_@FLYA | 130 | 119 (91.54%) | 158 | 106 (67.09%) | 1.1 | XEASY Sparky |
| N15NOESY_@FLYA | 3255 | 1745 (53.61%) | 1778 | 1272 (71.54%) | 1.4 | XEASY Sparky |
| ALL | 19127 | 8520 (44.54%) | 11790 | 5949 (50.46%) | 1.4 |
Table 1. Results of the automated shift assignment. Expected: Number of peaks expected in the spectrum. Assigned expected: Number of expected peaks that are assigned that to a measured (picked) peak. The percentage is relative to the number of expected peaks. Measured: Number of measured (picked) peaks. Assigned measured: Number of measured peaks that have one or more expected peaks assigned to it. The percentage is relative to the number of measured peaks. Ratio: The average number of expected peaks that are assigned to the same measured peak. The average is taken over all measured peaks to which at least on expected peak is assigned. Color coding: Rows corresponding to spectra with low percentage of assigned peaks (<15%) are highlighted in red. The most common factors associated with low percentage of assigned peaks are: (a) wrong spectra axes annotation (e.g. transposition of H-HC axes in C13NOESY, more information), (b) systematic shift between spectra (e.g. HCCHTOCSY shifted in C dimension by -2.7 ppm in reference to C13HSQC), (c) low number of signals that are visible in the spectrum (compared to number of signals that are expected to be present in the spectrum given the protein sequence), (d) inverted axis (inverted ppm scale in one of spectra dimensions).
Automated structure calculation results
Fig. Superposition of two structure proposals (proposal 1, proposal 2) determined by NMRtist (click image to show protein visualization at full scale). Structures are shown in three different orientations. The parts of the structure that are not well-defined are drawn with transparency. The two proposals are calculated with a different number of NOESY cross peaks. Download aligned structure proposals [combined.pdb] . If a "Test run" was executed, NMRtist may output only a single structure proposal.
| Property | Proposal 1 | Proposal 2 |
| Structure | structure.pdb | structure.pdb |
| Distance restraints | restraints.upl | restraints.upl |
| Structure calculation details | cyanatable.txt | cyanatable.txt |
| Total number of distance restraints | 2301 | 2506 |
| Number of intraresidual restraints (|i-j| = 0) | 506 | 544 |
| Number of sequential restraints (|i-j| = 1) | 593 | 663 |
| Number of medium-range restraints (1 < |i-j| < 5) | 495 | 534 |
| Number of long-range restraints (|i-j| ≥ 5) | 707 | 765 |
| Number of torsion angle restraints | 144 | 144 |
| CYANA target function value | 2.14 � 0.07 +� | 3.92 � 0.08 +� |
| Distance restraint violations > 0.2 Å | 8 � 1 | 16 � 2 |
| Maximal distance restraint violation | 0.28 � 0.07 + | 0.52 � 0.05 + |
| Angle restraint violations > 5.0° | 0 � 1 | 2 � 0 |
| Maximal angle restraint violation | 4.55 � 0.61 � | 6.74 � 0.20 � |
| Residues in most favored Ramachandran plot regions | 81.4 % | 79.1 % |
| Residues in additionally allowed regions | 18.6 % | 20.8 % |
| Residues in generously allowed regions | 0.0 % | 0.1 % |
| Residues in disallowed regions | 0.0 % | 0.0 % |
Table 1. Structure calculation statistics