************************************************************************ ********** REPORT OF PROTEIN ANALYSIS by the WHAT IF program ********** ************************************************************************ Date : 2006-02-04 This report was created by WHAT IF version 20030529-0952 INTRODUCTION ------------ This document contains a report of findings by the WHAT IF program during the analysis of one or more proteins. It contains a separate section for each of the proteins that have been analysed. Each reported fact has an assigned severity, one of: * error: severe errors encountered during the analyses. Items marked as errors are considered severe problems requiring immediate attention. * warning: Either less severe problems or uncommon structural features. These still need special attention. * note: Statistical values, plots, or other verbose results of tests and analyses that have been performed. If alternate conformations are present, only the first is evaluated. Hydrogen atoms are only included if explicitly requested, and even then they are not used by all checks. Legend ------ Some notations need a little explanation: RESIDUE: Residues in tables are normally given in 3-5 parts: - A number. This is the internal sequence number of the residue used by WHAT IF. - The residue name. Normally this is a three letter amino acid name. - The sequence number, between brackets. This is the residue number as it was given in the input file. It can be followed by the insertion code. - The chain identifier. A single character. If no chain identifier was given in the input file, this will be invisible. - A model number (only for NMR structures). Z-VALUE: To indicate the normality of a score, the score may be expressed as a Z-value or Z-score. This is just the number of standard deviations that the score deviates from the expected value. A property of Z-values is that the root-mean-square of a group of Z-values (the RMS Z-value) is expected to be 1.0. Z-values above 4.0 and below -4.0 are very uncommon. If a Z-score is used in WHAT IF, the accompanying text will explain how the expected value and standard deviation were obtained. ======================================================================== ==== Compound code /u/lytle/at5g39720/9valid/77/analyzed_input/analy====8.pdb ======================================================================== # 1 # Error: Missing unit cell information No SCALE matrix is given in the PDB file. # 2 # Error: Missing symmetry information Problem: No CRYST1 card is given in the PDB file. # 3 # Note: No rounded coordinates detected No significant rounding of atom coordinates has been detected. # 4 # Note: Valine nomenclature OK No errors were detected in valine nomenclature. # 5 # Note: Threonine nomenclature OK No errors were detected in threonine nomenclature. # 6 # Note: Isoleucine nomenclature OK No errors were detected in isoleucine nomenclature. # 7 # Note: Leucine nomenclature OK No errors were detected in leucine nomenclature. # 8 # Note: Arginine nomenclature OK No errors were detected in arginine nomenclature. # 9 # Warning: Tyrosine convention problem The tyrosine residues listed in the table below have their chi-2 not between -90.0 and 90.0 60 TYR ( 60 ) # 10 # Warning: Phenylalanine convention problem The phenylalanine residues listed in the table below have their chi-2 not between -90.0 and 90.0. 26 PHE ( 26 ) 147 PHE ( 147 ) # 11 # Warning: Aspartic acid convention problem The aspartic acid residues listed in the table below have their chi-2 not between -90.0 and 90.0, or their proton on OD1 instead of OD2. 13 ASP ( 13 ) 30 ASP ( 30 ) 88 ASP ( 88 ) 106 ASP ( 106 ) 164 ASP ( 164 ) # 12 # Warning: Glutamic acid convention problem The glutamic acid residues listed in the table below have their chi-3 outside the -90.0 to 90.0 range, or their proton on OE1 instead of OE2. 41 GLU ( 41 ) 70 GLU ( 70 ) 83 GLU ( 83 ) 85 GLU ( 85 ) 94 GLU ( 94 ) 96 GLU ( 96 ) 109 GLU ( 109 ) 129 GLU ( 129 ) 133 GLU ( 133 ) 134 GLU ( 134 ) 152 GLU ( 152 ) # 13 # Warning: Heavy atom naming problem The atoms listed in the table below have nonstandard names in the input file. (Be aware that we sometomes consider an asterix and an apostrophe identical, and thus do not warn for the use of asterixes. Swapped OP1 and OP2 on nucleic acid phosphors also are reported elsewhere). 173 GLN ( 173 ) O <--> OT1 # 14 # Note: Chirality OK All protein atoms have proper chirality. # 15 # Note: Improper dihedral angle distribution OK The RMS Z-score for all improper dihedrals in the structure is within normal ranges. Improper dihedral RMS Z-score : 0.862 # 16 # Note: Chain names are OK All chain names assigned to polymer molecules are unique, and all residue numbers are strictly increasing within each chain. # 17 # Note: Weights checked OK All atomic occupancy factors ('weights') fall in the 0.0--1.0 range. # 18 # Note: No missing atoms detected All expected atoms are present. # 19 # Note: OXT check OK All required C-terminal oxygen atoms are present. # 20 # Note: No extra C-terminal groups found No C-terminal groups are present for non C-terminal residues # 21 # Note: All bond lengths OK All bond lengths are in agreement with standard bond lengths using a tolerance of 4 sigma (both standard values and sigma for amino acid residues have been taken from Engh and Huber [REF], for DNA/RNA from Parkinson et al [REF]) # 22 # Warning: Low bond length variability Bond lengths were found to deviate less than normal from the mean Engh and Huber [REF] and/or Parkinson et al [REF] standard bond lengths. The RMS Z-score given below is expected to be around 1.0 for a normally restrained data set. The fact that it is lower than 0.667 in this structure might indicate that too-strong constraints have been used in the refinement. This can only be a problem for high resolution X-ray structures. RMS Z-score for bond lengths: 0.148 RMS-deviation in bond distances: 0.003 # 23 # Note: No bond length directionality Comparison of bond distances with Engh and Huber [REF] standard values for protein residues and Parkinson et al [REF] values for DNA/RNA does not show significant systematic deviations. # 24 # Note: All bond angles OK All bond angles are in agreement with standard bond angles using a tolerance of 4 sigma (both standard values and sigma for protein residues have been taken from Engh and Huber [REF], for DNA/RNA from Parkinson et al. [REF]). Please note that disulphide bridges are neglected. # 25 # Warning: Low bond angle variability Bond angles were found to deviate less than normal from the standard bond angles (normal values for protein residues were taken from Engh and Huber [REF], for DNA/RNA from Parkinson et al [REF]). The RMS Z-score given below is expected to be around 1.0 for a normally restrained data set. More common values are around 1.55. The fact that it is lower than 0.667 in this structure might indicate that too-strong constraints have been used in the refinement. This can only be a problem for high resolution X-ray structures. RMS Z-score for bond angles: 0.415 RMS-deviation in bond angles: 0.552 # 26 # Note: Side chain planarity OK All of the side chains of residues that have a planar group are planar within expected RMS deviations. # 27 # Note: Atoms connected to aromatic rings OK All of the atoms that are connected to planar aromatic rings in side chains of amino-acid residues are in the plane within expected RMS deviations. # 28 # Note: PRO puckering amplitude OK Puckering amplitudes for all PRO residues are within normal ranges. # 29 # Note: PRO puckering phases OK Puckering phases for all PRO residues are normal # 30 # Warning: Torsion angle evaluation shows unusual residues The residues listed in the table below contain bad or abnormal torsion angles. These scores give an impression of how ``normal'' the torsion angles in protein residues are. All torsion angles except omega are used for calculating a `normality' score. Average values and standard deviations were obtained from the residues in the WHAT IF database. These are used to calculate Z-scores. A residue with a Z-score of below -2.0 is poor, and a score of less than -3.0 is worrying. For such residues more than one torsion angle is in a highly unlikely position. 127 PHE ( 127 ) -3.3094 156 LYS ( 156 ) -2.8058 170 ARG ( 170 ) -2.7236 36 LEU ( 36 ) -2.6486 165 ILE ( 165 ) -2.5390 108 SER ( 108 ) -2.5330 56 LYS ( 56 ) -2.4525 4 HIS ( 4 ) -2.4169 58 ARG ( 58 ) -2.2725 166 SER ( 166 ) -2.2655 94 GLU ( 94 ) -2.2641 14 SER ( 14 ) -2.2560 57 GLY ( 57 ) -2.2316 171 GLU ( 171 ) -2.2141 3 HIS ( 3 ) -2.2090 23 TYR ( 23 ) -2.2018 115 THR ( 115 ) -2.1970 113 VAL ( 113 ) -2.1803 6 HIS ( 6 ) -2.1783 97 ARG ( 97 ) -2.1669 15 LEU ( 15 ) -2.1042 102 ILE ( 102 ) -2.0647 52 ARG ( 52 ) -2.0267 111 MET ( 111 ) -2.0214 # 31 # Warning: Backbone torsion angle evaluation shows unusual conformations The residues listed in the table below have abnormal backbone torsion angles. Residues with ``forbidden'' phi-psi combinations are listed, as well as residues with unusual omega angles (deviating by more than 3 sigma from the normal value). Please note that it is normal if about 5 percent of the residues is listed here as having unusual phi-psi combinations. 5 HIS ( 5 ) Poor phi/psi 13 ASP ( 13 ) Poor phi/psi 37 ASP ( 37 ) Poor phi/psi 56 LYS ( 56 ) Poor phi/psi 108 SER ( 108 ) Poor phi/psi 122 ALA ( 122 ) Poor phi/psi 127 PHE ( 127 ) Poor phi/psi 155 LYS ( 155 ) Poor phi/psi 156 LYS ( 156 ) Poor phi/psi 163 ASP ( 163 ) Poor phi/psi 165 ILE ( 165 ) Poor phi/psi 166 SER ( 166 ) Poor phi/psi 169 LEU ( 169 ) Poor phi/psi 170 ARG ( 170 ) Poor phi/psi # 32 # Error: Ramachandran Z-score very low The score expressing how well the backbone conformations of all residues are corresponding to the known allowed areas in the Ramachandran plot is very low. Ramachandran Z-score : -4.372 # 33 # Warning: Omega angles too tightly restrained The omega angles for trans-peptide bonds in a structure are expected to give a gaussian distribution with the average around +178 degrees and a standard deviation around 5.5 degrees. These expected values were obtained from very accurately determined structures. Many protein structures are too tightly constrained. This seems to be the case with the current structure, as the observed standard deviation is below 4.0 degrees. Standard deviation of omega values : 0.162 # 34 # Error: chi-1/chi-2 angle correlation Z-score very low The score expressing how well the chi-1/chi-2 angles of all residues are corresponding to the populated areas in the database is very low. chi-1/chi-2 correlation Z-score : -5.902 # 35 # Note: Ramachandran plot In this Ramachandran plot X-signs represent glycines, squares represent prolines and small plus-signs represent the other residues. If too many plus-signs fall outside the contoured areas then the molecule is poorly refined (or worse). In a colour picture, the residues that are part of a helix are shown in blue, strand residues in red. "Allowed" regions for helical residues are drawn in blue, for strand residues in red, and for all other residues in green. In the TeX file, a plot has been inserted here Chain without chain identifier # 36 # Note: Inside/Outside residue distribution normal The distribution of residue types over the inside and the outside of the protein is normal. inside/outside RMS Z-score : 1.139 # 37 # Note: Inside/Outside RMS Z-score plot The Inside/Outside distribution normality RMS Z-score over a 15 residue window is plotted as function of the residue number. High areas in the plot (above 1.5) indicate unusual inside/outside patterns. In the TeX file, a plot has been inserted here Chain without chain identifier # 38 # Note: Secondary structure This is the secondary structure according to DSSP. Only helix (H), strand (S), turn (T) and coil (blank) are shown. [REF] Secondary structure assignment 10 20 30 40 50 60 | | | | | | 1 - 60 GHHHHHHLECSSDSLQLHNVFVYGSFQDPDVINVMLDRTPEIVSATLPGFQRFRLKGRLY 1 - 60 T SSSS TTTT HHHHHHHTT SSSSSSSSSSS TT 70 80 90 100 110 120 | | | | | | 61 - 120 PCIVPSEKGEVHGKVLMGVTSDELENLDAVEGNEYERVTVGIVREDNSEKMAVKTYMWIN 61 - 120 SSSSSSSSSSSSSSSSS HHHHHHHHHHHTTT SSSSSSSSSTTT SSSSSSSS TT 130 140 150 160 170 | | | | | 121 - 173 KADPDMFGEWNFEEWKRLHKKKFIETFKKIMECKKKPQGQGNDDISHVLREDQ 121 - 173 TT TTT THHHHHHHHTHHHHHHHHHHHHHHT T # 39 # Error: Abnormally short interatomic distances The pairs of atoms listed in the table below have an unusually short distance. The contact distances of all atom pairs have been checked. Two atoms are said to `bump' if they are closer than the sum of their Van der Waals radii minus 0.40 Angstrom. For hydrogen bonded pairs a tolerance of 0.55 Angstrom is used. The first number in the table tells you how much shorter that specific contact is than the acceptable limit. The second distance is the distance between the centers of the two atoms. The last text-item on each line represents the status of the atom pair. The text `INTRA' means that the bump is between atoms that are explicitly listed in the PDB file. `INTER' means it is an inter-symmetry bump. If the final column contains the text 'HB', the bump criterium was relaxed because there could be a hydrogen bond. Similarly relaxed criteria are used for 1--3 and 1--4 interactions (listed as 'B2' and 'B3', respectively). If the last column is 'BF', the sum of the B-factors of the atoms is higher than 80, which makes the appearance of the bump somewhat less severe because the atoms probably aren't there anyway. Bumps between atoms for which the sum of their occupancies is lower than one are not reported. In any case, each bump is listed in only one direction. 90 VAL ( 90 ) CG1 -- 91 GLU ( 91 ) N 0.155 2.945 INTRA 163 ASP ( 163 ) C -- 164 ASP ( 164 ) CG 0.148 3.052 INTRA 163 ASP ( 163 ) CG -- 164 ASP ( 164 ) N 0.120 2.980 INTRA 49 GLY ( 49 ) C -- 50 PHE ( 50 ) CD1 0.113 3.087 INTRA BF 94 GLU ( 94 ) C -- 95 TYR ( 95 ) CD1 0.101 3.099 INTRA BF 50 PHE ( 50 ) C -- 71 VAL ( 71 ) CG2 0.090 3.110 INTRA 67 GLU ( 67 ) C -- 68 LYS ( 68 ) CD 0.081 3.119 INTRA BF 112 ALA ( 112 ) C -- 113 VAL ( 113 ) CG1 0.081 3.119 INTRA 140 LYS ( 140 ) O -- 141 LYS ( 141 ) C 0.055 2.745 INTRA BF 59 LEU ( 59 ) CG -- 60 TYR ( 60 ) N 0.055 3.045 INTRA 149 LYS ( 149 ) CG -- 162 ASN ( 162 ) N 0.046 3.054 INTRA 124 PRO ( 124 ) O -- 125 ASP ( 125 ) C 0.042 2.758 INTRA 73 GLY ( 73 ) C -- 74 LYS ( 74 ) CG 0.042 3.158 INTRA 63 ILE ( 63 ) C -- 64 VAL ( 64 ) CG1 0.039 3.161 INTRA BF 29 PRO ( 29 ) CB -- 39 THR ( 39 ) CG2 0.034 3.166 INTRA 91 GLU ( 91 ) C -- 93 ASN ( 93 ) N 0.028 2.872 INTRA BF 85 GLU ( 85 ) CG -- 86 ASN ( 86 ) N 0.027 3.073 INTRA 63 ILE ( 63 ) CG2 -- 64 VAL ( 64 ) N 0.027 3.073 INTRA 103 VAL ( 103 ) O -- 104 ARG ( 104 ) C 0.027 2.773 INTRA 130 TRP ( 130 ) C -- 132 PHE ( 132 ) N 0.024 2.876 INTRA 146 THR ( 146 ) CG2 -- 147 PHE ( 147 ) N 0.021 3.079 INTRA 130 TRP ( 130 ) O -- 131 ASN ( 131 ) C 0.021 2.779 INTRA 51 GLN ( 51 ) N -- 71 VAL ( 71 ) CG2 0.012 3.088 INTRA 54 ARG ( 54 ) CD -- 130 TRP ( 130 ) CD1 0.008 3.192 INTRA 18 HIS ( 18 ) ND1 -- 84 LEU ( 84 ) CD2 0.006 3.094 INTRA BF 118 TRP ( 118 ) CD1 -- 126 MET ( 126 ) SD 0.006 3.394 INTRA BF 144 ILE ( 144 ) O -- 145 GLU ( 145 ) C 0.005 2.795 INTRA BF 33 ASN ( 33 ) CG -- 39 THR ( 39 ) OG1 0.005 2.795 INTRA BF 155 LYS ( 155 ) CG -- 156 LYS ( 156 ) N 0.004 3.096 INTRA BF 60 TYR ( 60 ) CB -- 61 PRO ( 61 ) CD 0.004 3.196 INTRA 139 HIS ( 139 ) O -- 140 LYS ( 140 ) C 0.004 2.796 INTRA 97 ARG ( 97 ) CZ -- 114 LYS ( 114 ) CD 0.002 3.198 INTRA BF 100 VAL ( 100 ) CG1 -- 115 THR ( 115 ) CG2 0.002 3.198 INTRA 100 VAL ( 100 ) CG2 -- 115 THR ( 115 ) CG2 0.002 3.198 INTRA 23 TYR ( 23 ) CG -- 63 ILE ( 63 ) CG1 0.001 3.199 INTRA 23 TYR ( 23 ) CB -- 63 ILE ( 63 ) CG1 0.001 3.199 INTRA 118 TRP ( 118 ) CD1 -- 126 MET ( 126 ) CE 0.001 3.199 INTRA 51 GLN ( 51 ) CA -- 71 VAL ( 71 ) CG2 0.000 3.200 INTRA BF 163 ASP ( 163 ) O -- 164 ASP ( 164 ) CG 0.000 2.800 INTRA BF # 40 # Warning: Abnormal packing environment for some residues The residues listed in the table below have an unusual packing environment. The packing environment of the residues is compared with the average packing environment for all residues of the same type in good PDB files. A low packing score can indicate one of several things: Poor packing, misthreading of the sequence through the density, crystal contacts, contacts with a co-factor, or the residue is part of the active site. It is not uncommon to see a few of these, but in any case this requires further inspection of the residue. 170 ARG ( 170 ) -7.84 158 GLN ( 158 ) -7.83 160 GLN ( 160 ) -7.09 58 ARG ( 58 ) -6.92 60 TYR ( 60 ) -6.76 2 HIS ( 2 ) -6.45 5 HIS ( 5 ) -6.29 3 HIS ( 3 ) -6.19 169 LEU ( 169 ) -6.16 38 ARG ( 38 ) -6.13 7 HIS ( 7 ) -6.10 4 HIS ( 4 ) -5.98 10 CYS ( 10 ) -5.87 8 LEU ( 8 ) -5.85 67 GLU ( 67 ) -5.73 162 ASN ( 162 ) -5.67 56 LYS ( 56 ) -5.47 171 GLU ( 171 ) -5.34 6 HIS ( 6 ) -5.26 27 GLN ( 27 ) -5.24 155 LYS ( 155 ) -5.23 55 LEU ( 55 ) -5.19 156 LYS ( 156 ) -5.18 165 ILE ( 165 ) -5.14 121 LYS ( 121 ) -5.05 # 41 # Warning: Abnormal packing environment for sequential residues A stretch of at least three sequential residues with a questionable packing environment was found. This could indicate that these residues are part of a strange loop, but might also be an indication of misthreading. The table below lists the first and last residue in each stretch found, as well as the average residue score of the series. 2 HIS ( 2 ) --- 8 LEU ( 8 ) -6.02 54 ARG ( 54 ) --- 56 LYS ( 56 ) -4.99 154 LYS ( 154 ) --- 156 LYS ( 156 ) -4.88 158 GLN ( 158 ) --- 160 GLN ( 160 ) -6.32 164 ASP ( 164 ) --- 172 ASP ( 172 ) -5.23 # 42 # Error: Abnormal average packing environment The average quality control value for the structure is very low. A molecule is certain to be incorrect if the average quality score is below -3.0. Poorly refined molecules, very well energy minimized misthreaded molecules and low homology models give values between -2.0 and -3.0. The average quality of 200 highly refined Xray structures was -0.5+/-0.4 [REF]. Average for range 1 - 173 : -2.073 # 43 # Note: Quality value plot The quality value smoothed over a 10 residue window is plotted as function of the residue number. Low areas in the plot (below -2.0) indicate "unusual" packing. In the TeX file, a plot has been inserted here Chain without chain identifier # 44 # Warning: Low packing Z-score for some residues The residues listed in the table below have an unusual packing environment according to the 2nd generation quality check. The score listed in the table is a packing normality Z-score: positive means better than average, negative means worse than average. Only residues scoring less than -2.50 are listed here. These are the "unusual" residues in the structure, so it will be interesting to take a special look at them. 59 LEU ( 59 ) -2.70 # 45 # Note: No series of residues with abnormal new packing environment There are no stretches of four or more residues each having a quality control Z-score worse than -1.75. # 46 # Warning: Structural average packing Z-score a bit worrysome The structural 2nd generation average quality control value is a bit low. The protein is probably threaded correctly, but either poorly refined, or it is just a protein with an unusual (but correct) structure. The average quality of properly refined Xray structures is 0.0+/-1.0. All contacts : Average = -0.599 Z-score = -3.05 BB-BB contacts : Average = -0.343 Z-score = -2.10 BB-SC contacts : Average = -0.445 Z-score = -2.61 SC-BB contacts : Average = -0.629 Z-score = -3.81 SC-SC contacts : Average = -0.482 Z-score = -2.29 # 47 # Note: Second generation quality Z-score plot The second generation quality Z-score smoothed over a 10 residue window is plotted as function of the residue number. Low areas in the plot (below -1.3) indicate "unusual" packing. In the TeX file, a plot has been inserted here Chain without chain identifier # 48 # Warning: Backbone oxygen evaluation The residues listed in the table below have an unusual backbone oxygen position. For each of the residues in the structure, a search was performed to find 5-residue stretches in the WHAT IF database with superposable C-alpha coordinates, and some constraints on the neighboring backbone oxygens. In the following table the RMS distance between the backbone oxygen positions of these matching structures in the database and the position of the backbone oxygen atom in the current residue is given. If this number is larger than 1.5 a significant number of structures in the database show an alternative position for the backbone oxygen. If the number is larger than 2.0 most matching backbone fragments in the database have the peptide plane flipped. A manual check needs to be performed to assess whether the experimental data can support that alternative as well. The number in the last column is the number of database hits (maximum 80) used in the calculation. It is "normal" that some glycine residues show up in this list, but they are still worth checking! 106 ASP ( 106 ) 3.01 16 53 PHE ( 53 ) 2.97 80 22 VAL ( 22 ) 2.84 27 70 GLU ( 70 ) 2.74 49 151 MET ( 151 ) 2.58 22 141 LYS ( 141 ) 2.56 41 19 ASN ( 19 ) 2.55 27 18 HIS ( 18 ) 2.48 50 115 THR ( 115 ) 2.43 46 90 VAL ( 90 ) 2.33 31 101 GLY ( 101 ) 2.20 28 63 ILE ( 63 ) 2.19 52 147 PHE ( 147 ) 2.17 20 138 LEU ( 138 ) 2.12 80 146 THR ( 146 ) 2.06 22 51 GLN ( 51 ) 2.03 27 152 GLU ( 152 ) 2.02 64 32 ILE ( 32 ) 1.99 32 96 GLU ( 96 ) 1.98 17 136 LYS ( 136 ) 1.92 80 133 GLU ( 133 ) 1.92 62 132 PHE ( 132 ) 1.90 16 52 ARG ( 52 ) 1.89 41 148 LYS ( 148 ) 1.84 80 17 LEU ( 17 ) 1.81 80 45 ALA ( 45 ) 1.65 16 150 ILE ( 150 ) 1.64 28 102 ILE ( 102 ) 1.62 80 159 GLY ( 159 ) 1.57 14 104 ARG ( 104 ) 1.54 10 9 GLU ( 9 ) 1.51 10 # 49 # Note: Rotamers checked OK None of the residues that have a normal backbone environment have abnormal rotamers. # 50 # Warning: Unusual backbone conformations For the residues listed in the table below, the backbone formed by itself and two neighboring residues on either side is in a conformation that is not seen very often in the database of solved protein structures. The number given in the table is the number of similar backbone conformations in the database with the same amino acid in the center. For this check, backbone conformations are compared with database structures using C-alpha superpositions with some restraints on the backbone oxygen positions. A residue mentioned in the table can be part of a strange loop, or there might be something wrong with it or its directly surrounding residues. There are a few of these in every protein, but in any case it is worth looking at! 12 SER ( 12 ) 0 14 SER ( 14 ) 0 23 TYR ( 23 ) 0 24 GLY ( 24 ) 0 26 PHE ( 26 ) 0 36 LEU ( 36 ) 0 37 ASP ( 37 ) 0 56 LYS ( 56 ) 0 58 ARG ( 58 ) 0 68 LYS ( 68 ) 0 94 GLU ( 94 ) 0 121 LYS ( 121 ) 0 122 ALA ( 122 ) 0 126 MET ( 126 ) 0 127 PHE ( 127 ) 0 155 LYS ( 155 ) 0 156 LYS ( 156 ) 0 158 GLN ( 158 ) 0 162 ASN ( 162 ) 0 163 ASP ( 163 ) 0 165 ILE ( 165 ) 0 169 LEU ( 169 ) 0 170 ARG ( 170 ) 0 77 MET ( 77 ) 1 154 LYS ( 154 ) 1 164 ASP ( 164 ) 1 168 VAL ( 168 ) 1 38 ARG ( 38 ) 2 # 51 # Error: Backbone conformation Z-score very low A comparison of the backbone conformation with database proteins shows that the backbone fold in this structure is very unusual. Backbone conformation Z-score : -8.859 # 52 # Warning: Average B-factor problem The average B-factor for all buried protein atoms normally lies between 10--20. Values around 3--5 are expected for X-ray studies performed at liquid nitrogen temperature. Because of the extreme value for the average B-factor, no further analysis of the B-factors is performed. Average B-factor for buried atoms : 37.079 # 53 # Note: B-factor plot The average atomic B-factor per residue is plotted as function of the residue number. In the TeX file, a plot has been inserted here Chain without chain identifier # 54 # Error: HIS, ASN, GLN side chain flips Listed here are Histidine, Asparagine or Glutamine residues for which the orientation determined from hydrogen bonding analysis are different from the assignment given in the input. Either they could form energetically more favorable hydrogen bonds if the terminal group was rotated by 180 degrees, or there is no assignment in the input file (atom type 'A') but an assignment could be made. If a residue is marked ``flexible'' the flipped conformation is only slightly better than the non-flipped conformation. 4 HIS ( 4 ) 7 HIS ( 7 ) 19 ASN ( 19 ) 33 ASN ( 33 ) 86 ASN ( 86 ) 93 ASN ( 93 ) # 55 # Note: Histidine type assignments For all complete HIS residues in the structure a tentative assignment to HIS-D (protonated on ND1), HIS-E (protonated on NE2), or HIS-H (protonated on both ND1 and NE2, positively charged) is made based on the hydrogen bond network. A second assignment is made based on which of the Engh and Huber [REF] histidine geometries fits best to the structure. In the table below all normal histidine residues are listed. The assignment based on the geometry of the residue is listed first, together with the RMS Z-score for the fit to the Engh and Huber parameters. For all residues where the H-bond assignment is different, the assignment is listed in the last columns, together with its RMS Z-score to the Engh and Huber parameters. As always, the RMS Z-scores should be close to 1.0 if the residues were restrained to the Engh and Huber parameters during refinement. Please note that because the differences between the geometries of the different types are small it is possible that the geometric assignment given here does not correspond to the type used in refinement. This is especially true if the RMS Z-scores are much higher than 1.0. If the two assignments differ, or the ``geometry'' RMS Z-score is high, it is advisable to verify the hydrogen bond assignment, check the HIS type used during the refinement and possibly adjust it. 2 HIS ( 2 ) HIS-H 0.07 HIS-E 0.59 3 HIS ( 3 ) HIS-H 0.05 HIS-E 0.58 4 HIS ( 4 ) HIS-H 0.07 HIS-D 0.56 5 HIS ( 5 ) HIS-H 0.06 HIS-E 0.60 6 HIS ( 6 ) HIS-H 0.06 HIS-E 0.60 7 HIS ( 7 ) HIS-H 0.07 HIS-D 0.55 18 HIS ( 18 ) HIS-H 0.06 HIS-E 0.58 72 HIS ( 72 ) HIS-H 0.07 HIS-D 0.56 139 HIS ( 139 ) HIS-H 0.04 HIS-E 0.58 167 HIS ( 167 ) HIS-H 0.08 HIS-E 0.60 # 56 # Warning: Buried unsatisfied hydrogen bond donors The buried hydrogen bond donors listed in the table below have a hydrogen atom that is not involved in a hydrogen bond in the optimized hydrogen bond network. Hydrogen bond donors that are buried inside the protein normally use all of their hydrogens to form hydrogen bonds within the protein. If there are any non hydrogen bonded buried hydrogen bond donors in the structure they will be listed here. In very good structures the number of listed atoms will tend to zero. 19 ASN ( 19 ) N 23 TYR ( 23 ) N 24 GLY ( 24 ) N 27 GLN ( 27 ) N 31 VAL ( 31 ) N 37 ASP ( 37 ) N 46 THR ( 46 ) OG1 52 ARG ( 52 ) NE 55 LEU ( 55 ) N 56 LYS ( 56 ) N 58 ARG ( 58 ) N 59 LEU ( 59 ) N 60 TYR ( 60 ) N 63 ILE ( 63 ) N 70 GLU ( 70 ) N 74 LYS ( 74 ) N 86 ASN ( 86 ) ND2 93 ASN ( 93 ) N 97 ARG ( 97 ) NE 104 ARG ( 104 ) NE 106 ASP ( 106 ) N 107 ASN ( 107 ) N 109 GLU ( 109 ) N 118 TRP ( 118 ) NE1 119 ILE ( 119 ) N 125 ASP ( 125 ) N 126 MET ( 126 ) N 128 GLY ( 128 ) N # 57 # Note: Buried hydrogen bond acceptors OK All buried polar side-chain hydrogen bond acceptors are involved in a hydrogen bond in the optimized hydrogen bond network. # 58 # Note: Summary report for users of a structure This is an overall summary of the quality of the structure as compared with current reliable structures. This summary is most useful for biologists seeking a good structure to use for modelling calculations. The second part of the table mostly gives an impression of how well the model conforms to common refinement constraint values. The first part of the table shows a number of constraint-independent quality indicators. Structure Z-scores, positive is better than average: 1st generation packing quality : -3.933 (poor) 2nd generation packing quality : -3.052 (poor) Ramachandran plot appearance : -4.372 (bad) chi-1/chi-2 rotamer normality : -5.902 (bad) Backbone conformation : -8.859 (bad) RMS Z-scores, should be close to 1.0: Bond lengths : 0.148 (tight) Bond angles : 0.415 (tight) Omega angle restraints : 0.029 (tight) Side chain planarity : 0.041 (tight) Improper dihedral distribution : 0.862 Inside/Outside distribution : 1.139 REFERENCES ========== WHAT IF G.Vriend, WHAT IF: a molecular modelling and drug design program, J. Mol. Graph. 8, 52--56 (1990). WHAT_CHECK (verification routines from WHAT IF) R.W.W.Hooft, G.Vriend, C.Sander and E.E.Abola, Errors in protein structures Nature 381, 272 (1996). Bond lengths and angles, protein residues R.Engh and R.Huber, Accurate bond and angle parameters for X-ray protein structure refinement, Acta Crystallogr. A47, 392--400 (1991). Bond lengths and angles, DNA/RNA G.Parkinson, J.Voitechovsky, L.Clowney, A.T.Bruenger and H.Berman, New parameters for the refinement of nucleic acid-containing structures Acta Crystallogr. D52, 57--64 (1996). DSSP W.Kabsch and C.Sander, Dictionary of protein secondary structure: pattern recognition of hydrogen bond and geometrical features Biopolymers 22, 2577--2637 (1983). Hydrogen bond networks R.W.W.Hooft, C.Sander and G.Vriend, Positioning hydrogen atoms by optimizing hydrogen bond networks in protein structures PROTEINS, 26, 363--376 (1996). Matthews' Coefficient B.W.Matthews Solvent content of Protein Crystals J. Mol. Biol. 33, 491--497 (1968). Protein side chain planarity R.W.W. Hooft, C. Sander and G. Vriend, Verification of protein structures: side-chain planarity J. Appl. Cryst. 29, 714--716 (1996). Puckering parameters D.Cremer and J.A.Pople, A general definition of ring puckering coordinates J. Am. Chem. Soc. 97, 1354--1358 (1975). Quality Control G.Vriend and C.Sander, Quality control of protein models: directional atomic contact analysis, J. Appl. Cryst. 26, 47--60 (1993). Ramachandran plot G.N.Ramachandran, C.Ramakrishnan and V.Sasisekharan, Stereochemistry of Polypeptide Chain Conformations J. Mol. Biol. 7, 95--99 (1963). Symmetry Checks R.W.W.Hooft, C.Sander and G.Vriend, Reconstruction of symmetry related molecules from protein data bank (PDB) files J. Appl. Cryst. 27, 1006--1009 (1994).