************************************************************************
********** REPORT OF PROTEIN ANALYSIS by the WHAT IF program **********
************************************************************************
Date : 2006-01-19
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/58g/analyzed_input/anal====_9.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
95 TYR ( 95 )
# 10 # Warning: Phenylalanine convention problem
The phenylalanine residues listed in the table below have their
chi-2 not between -90.0 and 90.0.
50 PHE ( 50 )
53 PHE ( 53 )
132 PHE ( 132 )
143 PHE ( 143 )
# 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 )
28 ASP ( 28 )
30 ASP ( 30 )
37 ASP ( 37 )
88 ASP ( 88 )
106 ASP ( 106 )
123 ASP ( 123 )
125 ASP ( 125 )
163 ASP ( 163 )
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.
9 GLU ( 9 )
41 GLU ( 41 )
70 GLU ( 70 )
91 GLU ( 91 )
105 GLU ( 105 )
129 GLU ( 129 )
133 GLU ( 133 )
134 GLU ( 134 )
# 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 <--> O2
# 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.880
# 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.147
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.416
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.2840
167 HIS ( 167 ) -2.8608
7 HIS ( 7 ) -2.8600
36 LEU ( 36 ) -2.7598
6 HIS ( 6 ) -2.7495
102 ILE ( 102 ) -2.5511
156 LYS ( 156 ) -2.5055
171 GLU ( 171 ) -2.4847
94 GLU ( 94 ) -2.4592
9 GLU ( 9 ) -2.4458
162 ASN ( 162 ) -2.3625
55 LEU ( 55 ) -2.3488
15 LEU ( 15 ) -2.3420
110 LYS ( 110 ) -2.3112
74 LYS ( 74 ) -2.2919
58 ARG ( 58 ) -2.2800
23 TYR ( 23 ) -2.2739
166 SER ( 166 ) -2.1987
97 ARG ( 97 ) -2.1755
107 ASN ( 107 ) -2.1702
172 ASP ( 172 ) -2.1395
121 LYS ( 121 ) -2.0925
70 GLU ( 70 ) -2.0741
113 VAL ( 113 ) -2.0715
104 ARG ( 104 ) -2.0409
4 HIS ( 4 ) -2.0356
169 LEU ( 169 ) -2.0295
13 ASP ( 13 ) -2.0287
# 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.
13 ASP ( 13 ) Poor phi/psi
36 LEU ( 36 ) Poor phi/psi
37 ASP ( 37 ) Poor phi/psi
94 GLU ( 94 ) Poor phi/psi
122 ALA ( 122 ) Poor phi/psi
127 PHE ( 127 ) Poor phi/psi
131 ASN ( 131 ) Poor phi/psi
156 LYS ( 156 ) Poor phi/psi
172 ASP ( 172 ) 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.398
# 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.016
# 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.854
# 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.113
# 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 T SSSS T333THHHHHHHTT SSSSSSSSSSS TTTT
70 80 90 100 110 120
| | | | | |
61 - 120 PCIVPSEKGEVHGKVLMGVTSDELENLDAVEGNEYERVTVGIVREDNSEKMAVKTYMWIN
61 - 120 SSSSSSSSSSSSSSSSS HHHHHHHHHHHTTT SSSSSSSSSTTTT SSSSSSSS TT
130 140 150 160 170
| | | | |
121 - 173 KADPDMFGEWNFEEWKRLHKKKFIETFKKIMECKKKPQGQGNDDISHVLREDQ
121 - 173 TT TTTHHHHHHHHHHHHHHHHHHHHHH 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.
130 TRP ( 130 ) C -- 131 ASN ( 131 ) CG 0.211 2.989 INTRA BF
120 ASN ( 120 ) O -- 121 LYS ( 121 ) C 0.211 2.589 INTRA BF
120 ASN ( 120 ) C -- 121 LYS ( 121 ) CG 0.125 3.075 INTRA
91 GLU ( 91 ) CD -- 116 TYR ( 116 ) CD1 0.115 3.085 INTRA BF
63 ILE ( 63 ) C -- 64 VAL ( 64 ) CG1 0.103 3.097 INTRA
94 GLU ( 94 ) C -- 95 TYR ( 95 ) CD1 0.084 3.116 INTRA BF
31 VAL ( 31 ) CG2 -- 32 ILE ( 32 ) N 0.077 3.023 INTRA BF
162 ASN ( 162 ) C -- 163 ASP ( 163 ) CG 0.076 3.124 INTRA
63 ILE ( 63 ) CG2 -- 64 VAL ( 64 ) N 0.066 3.034 INTRA BF
133 GLU ( 133 ) CG -- 134 GLU ( 134 ) N 0.063 3.037 INTRA
133 GLU ( 133 ) CD -- 134 GLU ( 134 ) N 0.060 3.040 INTRA
130 TRP ( 130 ) O -- 131 ASN ( 131 ) CB 0.060 2.740 INTRA BF
103 VAL ( 103 ) O -- 104 ARG ( 104 ) C 0.059 2.741 INTRA BF
112 ALA ( 112 ) C -- 113 VAL ( 113 ) CG1 0.058 3.142 INTRA BF
140 LYS ( 140 ) O -- 141 LYS ( 141 ) C 0.058 2.742 INTRA
59 LEU ( 59 ) CG -- 60 TYR ( 60 ) N 0.056 3.044 INTRA
104 ARG ( 104 ) CB -- 107 ASN ( 107 ) ND2 0.055 3.045 INTRA
107 ASN ( 107 ) ND2 -- 109 GLU ( 109 ) CG 0.054 3.046 INTRA BF
55 LEU ( 55 ) CD1 -- 56 LYS ( 56 ) N 0.053 3.047 INTRA
91 GLU ( 91 ) C -- 93 ASN ( 93 ) N 0.053 2.847 INTRA
79 VAL ( 79 ) CG1 -- 84 LEU ( 84 ) N 0.052 3.048 INTRA BF
165 ILE ( 165 ) CG1 -- 166 SER ( 166 ) N 0.050 3.050 INTRA
152 GLU ( 152 ) CG -- 153 CYS ( 153 ) N 0.048 3.052 INTRA BF
29 PRO ( 29 ) CB -- 39 THR ( 39 ) CG2 0.045 3.155 INTRA
19 ASN ( 19 ) CG -- 77 MET ( 77 ) CB 0.031 3.169 INTRA
146 THR ( 146 ) CG2 -- 147 PHE ( 147 ) N 0.027 3.073 INTRA BF
35 MET ( 35 ) CE -- 36 LEU ( 36 ) CD2 0.026 3.174 INTRA
144 ILE ( 144 ) O -- 145 GLU ( 145 ) C 0.022 2.778 INTRA
42 ILE ( 42 ) CG1 -- 76 LEU ( 76 ) CD2 0.018 3.182 INTRA BF
131 ASN ( 131 ) O -- 132 PHE ( 132 ) C 0.017 2.783 INTRA
84 LEU ( 84 ) CD1 -- 97 ARG ( 97 ) NH2 0.016 3.084 INTRA
98 VAL ( 98 ) CG1 -- 117 MET ( 117 ) CE 0.014 3.186 INTRA
125 ASP ( 125 ) O -- 126 MET ( 126 ) C 0.012 2.788 INTRA
80 THR ( 80 ) O -- 81 SER ( 81 ) C 0.009 2.791 INTRA
98 VAL ( 98 ) CG2 -- 117 MET ( 117 ) SD 0.007 3.393 INTRA
100 VAL ( 100 ) CG2 -- 115 THR ( 115 ) CG2 0.006 3.194 INTRA BF
94 GLU ( 94 ) CG -- 95 TYR ( 95 ) CE1 0.004 3.196 INTRA
154 LYS ( 154 ) O -- 155 LYS ( 155 ) C 0.003 2.797 INTRA BF
78 GLY ( 78 ) O -- 79 VAL ( 79 ) C 0.003 2.797 INTRA BF
19 ASN ( 19 ) ND2 -- 111 MET ( 111 ) CE 0.003 3.097 INTRA BF
83 GLU ( 83 ) O -- 84 LEU ( 84 ) C 0.002 2.798 INTRA
23 TYR ( 23 ) CG -- 63 ILE ( 63 ) CG1 0.002 3.198 INTRA
124 PRO ( 124 ) O -- 125 ASP ( 125 ) C 0.001 2.799 INTRA
# 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.
160 GLN ( 160 ) -7.34
158 GLN ( 158 ) -7.27
170 ARG ( 170 ) -7.13
8 LEU ( 8 ) -6.71
38 ARG ( 38 ) -6.54
3 HIS ( 3 ) -6.50
2 HIS ( 2 ) -6.46
169 LEU ( 169 ) -6.41
4 HIS ( 4 ) -6.33
6 HIS ( 6 ) -6.31
5 HIS ( 5 ) -6.17
7 HIS ( 7 ) -6.14
60 TYR ( 60 ) -6.11
68 LYS ( 68 ) -6.03
58 ARG ( 58 ) -5.98
9 GLU ( 9 ) -5.93
56 LYS ( 56 ) -5.92
167 HIS ( 167 ) -5.90
67 GLU ( 67 ) -5.84
171 GLU ( 171 ) -5.83
162 ASN ( 162 ) -5.74
165 ILE ( 165 ) -5.56
127 PHE ( 127 ) -5.54
27 GLN ( 27 ) -5.46
23 TYR ( 23 ) -5.25
168 VAL ( 168 ) -5.02
# 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.
1 GLY ( 1 ) --- 10 CYS ( 10 ) -5.98
155 LYS ( 155 ) --- 172 ASP ( 172 ) -5.32
# 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.103
# 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 # Note: Second generation packing environment OK
None of the individual amino acid residues has a bad packing environment.
# 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 # Note: Structural average packing Z-score OK
The structural average for the second generation quality control
value is within normal ranges.
All contacts : Average = -0.542 Z-score = -2.77
BB-BB contacts : Average = -0.299 Z-score = -1.84
BB-SC contacts : Average = -0.371 Z-score = -2.20
SC-BB contacts : Average = -0.579 Z-score = -3.52
SC-SC contacts : Average = -0.516 Z-score = -2.46
# 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!
55 LEU ( 55 ) 3.30 42
102 ILE ( 102 ) 3.09 69
115 THR ( 115 ) 2.68 64
53 PHE ( 53 ) 2.65 80
101 GLY ( 101 ) 2.58 25
26 PHE ( 26 ) 2.57 11
141 LYS ( 141 ) 2.50 35
90 VAL ( 90 ) 2.46 44
18 HIS ( 18 ) 2.46 34
140 LYS ( 140 ) 2.44 59
112 ALA ( 112 ) 2.39 17
150 ILE ( 150 ) 2.35 30
132 PHE ( 132 ) 2.33 17
86 ASN ( 86 ) 2.31 59
32 ILE ( 32 ) 2.26 36
151 MET ( 151 ) 2.23 24
19 ASN ( 19 ) 2.19 24
63 ILE ( 63 ) 2.19 52
146 THR ( 146 ) 2.16 40
51 GLN ( 51 ) 2.14 32
148 LYS ( 148 ) 2.10 58
147 PHE ( 147 ) 2.09 22
96 GLU ( 96 ) 2.02 20
152 GLU ( 152 ) 2.02 48
133 GLU ( 133 ) 1.99 45
116 TYR ( 116 ) 1.98 34
70 GLU ( 70 ) 1.94 21
111 MET ( 111 ) 1.90 31
135 TRP ( 135 ) 1.81 14
91 GLU ( 91 ) 1.80 23
110 LYS ( 110 ) 1.78 67
52 ARG ( 52 ) 1.76 59
134 GLU ( 134 ) 1.58 52
30 ASP ( 30 ) 1.53 39
104 ARG ( 104 ) 1.51 12
# 49 # Warning: Unusual rotamers
The residues listed in the table below have a rotamer that is not
seen very often in the database of solved protein structures. This
option determines for every residue the position specific chi-1
rotamer distribution. Thereafter it verified whether the actual
residue in the molecule has the most preferred rotamer or not. If
the actual rotamer is the preferred one, the score is 1.0. If the
actual rotamer is unique, the score is 0.0. If there are two
preferred rotamers, with a population distribution of 3:2 and your
rotamer sits in the lesser populated rotamer, the score will be
0.667. No value will be given if insufficient hits are found in the
database.
It is not necessarily an error if a few residues have rotamer
values below 0.3, but careful inspection of all residues with these
low values could be worth it.
158 GLN ( 158 ) 0.36
# 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!
23 TYR ( 23 ) 0
36 LEU ( 36 ) 0
37 ASP ( 37 ) 0
56 LYS ( 56 ) 0
58 ARG ( 58 ) 0
77 MET ( 77 ) 0
93 ASN ( 93 ) 0
94 GLU ( 94 ) 0
120 ASN ( 120 ) 0
121 LYS ( 121 ) 0
122 ALA ( 122 ) 0
126 MET ( 126 ) 0
127 PHE ( 127 ) 0
131 ASN ( 131 ) 0
156 LYS ( 156 ) 0
163 ASP ( 163 ) 0
10 CYS ( 10 ) 1
35 MET ( 35 ) 1
59 LEU ( 59 ) 1
68 LYS ( 68 ) 1
130 TRP ( 130 ) 1
155 LYS ( 155 ) 1
164 ASP ( 164 ) 1
171 GLU ( 171 ) 1
12 SER ( 12 ) 2
25 SER ( 25 ) 2
108 SER ( 108 ) 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 : -9.621
# 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 : 38.875
# 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.
7 HIS ( 7 )
33 ASN ( 33 )
162 ASN ( 162 )
167 HIS ( 167 )
# 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.08 HIS-E 0.60
3 HIS ( 3 ) HIS-H 0.05 HIS-E 0.59
4 HIS ( 4 ) HIS-H 0.08 HIS-D 0.55
5 HIS ( 5 ) HIS-H 0.05 HIS-E 0.59
6 HIS ( 6 ) HIS-H 0.07 HIS-D 0.56
7 HIS ( 7 ) HIS-H 0.06 HIS-E 0.60
18 HIS ( 18 ) HIS-H 0.03 HIS-E 0.58
72 HIS ( 72 ) HIS-H 0.06 HIS-E 0.59
139 HIS ( 139 ) HIS-H 0.07 HIS-E 0.58
167 HIS ( 167 ) HIS-H 0.09 HIS-E 0.59
# 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.
26 PHE ( 26 ) N
31 VAL ( 31 ) N
37 ASP ( 37 ) N
50 PHE ( 50 ) N
54 ARG ( 54 ) N
59 LEU ( 59 ) N
60 TYR ( 60 ) N
63 ILE ( 63 ) N
74 LYS ( 74 ) N
96 GLU ( 96 ) N
97 ARG ( 97 ) NE
107 ASN ( 107 ) N
107 ASN ( 107 ) ND2
109 GLU ( 109 ) N
116 TYR ( 116 ) OH
117 MET ( 117 ) N
119 ILE ( 119 ) N
121 LYS ( 121 ) N
125 ASP ( 125 ) N
126 MET ( 126 ) N
129 GLU ( 129 ) N
130 TRP ( 130 ) N
130 TRP ( 130 ) NE1
133 GLU ( 133 ) N
# 57 # Warning: Buried unsatisfied hydrogen bond acceptors
The buried side-chain hydrogen bond acceptors listed in the table
below are not involved in a hydrogen bond in the optimized hydrogen
bond network.
Side-chain hydrogen bond acceptors that are buried inside the
protein normally form hydrogen bonds within the protein. If there
are any not hydrogen bonded in the optimized hydrogen bond network
they will be listed here.
91 GLU ( 91 ) OE1
91 GLU ( 91 ) OE2
# 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 : -4.007 (poor)
2nd generation packing quality : -2.771
Ramachandran plot appearance : -4.398 (bad)
chi-1/chi-2 rotamer normality : -5.854 (bad)
Backbone conformation : -9.621 (bad)
RMS Z-scores, should be close to 1.0:
Bond lengths : 0.147 (tight)
Bond angles : 0.416 (tight)
Omega angle restraints : 0.003 (tight)
Side chain planarity : 0.040 (tight)
Improper dihedral distribution : 0.880
Inside/Outside distribution : 1.113
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).