************************************************************************
********** 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====6.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.
21 PHE ( 21 )
26 PHE ( 26 )
53 PHE ( 53 )
127 PHE ( 127 )
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 )
37 ASP ( 37 )
82 ASP ( 82 )
88 ASP ( 88 )
106 ASP ( 106 )
125 ASP ( 125 )
163 ASP ( 163 )
172 ASP ( 172 )
# 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 )
83 GLU ( 83 )
85 GLU ( 85 )
91 GLU ( 91 )
94 GLU ( 94 )
105 GLU ( 105 )
145 GLU ( 145 )
171 GLU ( 171 )
# 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.871
# 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.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.
2 HIS ( 2 ) -3.2389
121 LYS ( 121 ) -2.7713
36 LEU ( 36 ) -2.7114
56 LYS ( 56 ) -2.5975
108 SER ( 108 ) -2.5776
115 THR ( 115 ) -2.4531
5 HIS ( 5 ) -2.4295
167 HIS ( 167 ) -2.3512
54 ARG ( 54 ) -2.3072
102 ILE ( 102 ) -2.1826
110 LYS ( 110 ) -2.1729
126 MET ( 126 ) -2.1528
156 LYS ( 156 ) -2.0982
129 GLU ( 129 ) -2.0906
113 VAL ( 113 ) -2.0797
171 GLU ( 171 ) -2.0396
100 VAL ( 100 ) -2.0203
57 GLY ( 57 ) -2.0102
70 GLU ( 70 ) -2.0002
# 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.
2 HIS ( 2 ) Poor phi/psi
24 GLY ( 24 ) Poor phi/psi
37 ASP ( 37 ) Poor phi/psi
56 LYS ( 56 ) Poor phi/psi
168 VAL ( 168 ) 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.104
# 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.115
# 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.856
# 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.135
# 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 SSSSSTTTT THHHHHHHTT SS SSSSSSSS TTT
70 80 90 100 110 120
| | | | | |
61 - 120 PCIVPSEKGEVHGKVLMGVTSDELENLDAVEGNEYERVTVGIVREDNSEKMAVKTYMWIN
61 - 120 SSSSSSSSSS SSSSSS HHHHHHHHHHHTTT SSSSSSSSSTTTTSSSSSSSSS TT
130 140 150 160 170
| | | | |
121 - 173 KADPDMFGEWNFEEWKRLHKKKFIETFKKIMECKKKPQGQGNDDISHVLREDQ
121 - 173 333 HHHHHHHHHHHHHHHHHHHHHHHHT
# 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.
120 ASN ( 120 ) O -- 121 LYS ( 121 ) C 0.224 2.576 INTRA BF
120 ASN ( 120 ) C -- 121 LYS ( 121 ) CD 0.125 3.075 INTRA BF
107 ASN ( 107 ) CG -- 109 GLU ( 109 ) CG 0.107 3.093 INTRA BF
50 PHE ( 50 ) C -- 71 VAL ( 71 ) CG2 0.105 3.095 INTRA BF
25 SER ( 25 ) C -- 26 PHE ( 26 ) CG 0.096 3.104 INTRA
112 ALA ( 112 ) C -- 113 VAL ( 113 ) CG1 0.089 3.111 INTRA
49 GLY ( 49 ) C -- 50 PHE ( 50 ) CD1 0.078 3.122 INTRA BF
31 VAL ( 31 ) CG1 -- 32 ILE ( 32 ) N 0.067 3.033 INTRA
168 VAL ( 168 ) CG1 -- 169 LEU ( 169 ) N 0.061 3.039 INTRA BF
140 LYS ( 140 ) O -- 141 LYS ( 141 ) C 0.055 2.745 INTRA BF
23 TYR ( 23 ) CD2 -- 24 GLY ( 24 ) N 0.054 3.046 INTRA BF
19 ASN ( 19 ) ND2 -- 111 MET ( 111 ) CE 0.053 3.047 INTRA
146 THR ( 146 ) CG2 -- 147 PHE ( 147 ) N 0.053 3.047 INTRA BF
75 VAL ( 75 ) CG1 -- 76 LEU ( 76 ) N 0.050 3.050 INTRA
91 GLU ( 91 ) C -- 93 ASN ( 93 ) N 0.049 2.851 INTRA BF
73 GLY ( 73 ) C -- 74 LYS ( 74 ) CG 0.044 3.156 INTRA BF
96 GLU ( 96 ) CG -- 97 ARG ( 97 ) N 0.038 3.062 INTRA
67 GLU ( 67 ) C -- 68 LYS ( 68 ) CG 0.038 3.162 INTRA
120 ASN ( 120 ) C -- 121 LYS ( 121 ) CG 0.032 3.168 INTRA BF
59 LEU ( 59 ) CG -- 60 TYR ( 60 ) N 0.026 3.074 INTRA BF
25 SER ( 25 ) C -- 26 PHE ( 26 ) CD1 0.019 3.181 INTRA
54 ARG ( 54 ) CB -- 130 TRP ( 130 ) CD1 0.014 3.186 INTRA BF
35 MET ( 35 ) CE -- 87 LEU ( 87 ) CD2 0.010 3.190 INTRA BF
98 VAL ( 98 ) CG1 -- 117 MET ( 117 ) CE 0.006 3.194 INTRA BF
155 LYS ( 155 ) O -- 156 LYS ( 156 ) C 0.005 2.795 INTRA BF
62 CYS ( 62 ) SG -- 126 MET ( 126 ) SD 0.003 3.597 INTRA
29 PRO ( 29 ) CB -- 39 THR ( 39 ) CG2 0.003 3.197 INTRA BF
144 ILE ( 144 ) O -- 145 GLU ( 145 ) C 0.002 2.798 INTRA
50 PHE ( 50 ) O -- 69 GLY ( 69 ) C 0.002 2.798 INTRA BF
26 PHE ( 26 ) CE2 -- 35 MET ( 35 ) SD 0.001 3.399 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.
60 TYR ( 60 ) -7.48
55 LEU ( 55 ) -7.02
15 LEU ( 15 ) -6.53
58 ARG ( 58 ) -6.52
2 HIS ( 2 ) -6.50
170 ARG ( 170 ) -6.40
127 PHE ( 127 ) -6.25
4 HIS ( 4 ) -6.17
56 LYS ( 56 ) -6.15
23 TYR ( 23 ) -6.06
67 GLU ( 67 ) -5.91
7 HIS ( 7 ) -5.84
160 GLN ( 160 ) -5.84
5 HIS ( 5 ) -5.83
169 LEU ( 169 ) -5.82
38 ARG ( 38 ) -5.71
167 HIS ( 167 ) -5.67
68 LYS ( 68 ) -5.47
171 GLU ( 171 ) -5.21
8 LEU ( 8 ) -5.05
110 LYS ( 110 ) -5.01
# 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 ) --- 5 HIS ( 5 ) -5.40
7 HIS ( 7 ) --- 10 CYS ( 10 ) -5.03
54 ARG ( 54 ) --- 56 LYS ( 56 ) -5.99
58 ARG ( 58 ) --- 60 TYR ( 60 ) -6.12
159 GLY ( 159 ) --- 161 GLY ( 161 ) -4.80
169 LEU ( 169 ) --- 172 ASP ( 172 ) -5.52
# 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.019
# 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 # 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.600 Z-score = -3.06
BB-BB contacts : Average = -0.338 Z-score = -2.07
BB-SC contacts : Average = -0.478 Z-score = -2.80
SC-BB contacts : Average = -0.621 Z-score = -3.77
SC-SC contacts : Average = -0.531 Z-score = -2.53
# 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!
164 ASP ( 164 ) 3.27 10
170 ARG ( 170 ) 3.11 27
165 ILE ( 165 ) 2.89 24
10 CYS ( 10 ) 2.81 57
3 HIS ( 3 ) 2.81 18
102 ILE ( 102 ) 2.76 78
18 HIS ( 18 ) 2.66 39
19 ASN ( 19 ) 2.64 25
125 ASP ( 125 ) 2.57 25
141 LYS ( 141 ) 2.49 43
32 ILE ( 32 ) 2.46 36
140 LYS ( 140 ) 2.44 58
90 VAL ( 90 ) 2.40 38
8 LEU ( 8 ) 2.40 80
101 GLY ( 101 ) 2.40 24
70 GLU ( 70 ) 2.33 36
51 GLN ( 51 ) 2.26 33
147 PHE ( 147 ) 2.22 19
52 ARG ( 52 ) 2.21 65
133 GLU ( 133 ) 2.21 69
22 VAL ( 22 ) 2.20 40
9 GLU ( 9 ) 2.16 38
153 CYS ( 153 ) 2.08 80
148 LYS ( 148 ) 2.06 72
139 HIS ( 139 ) 2.05 80
135 TRP ( 135 ) 1.96 41
21 PHE ( 21 ) 1.94 69
146 THR ( 146 ) 1.90 21
39 THR ( 39 ) 1.88 54
104 ARG ( 104 ) 1.82 13
4 HIS ( 4 ) 1.69 11
53 PHE ( 53 ) 1.67 80
114 LYS ( 114 ) 1.66 10
144 ILE ( 144 ) 1.62 30
30 ASP ( 30 ) 1.58 16
91 GLU ( 91 ) 1.58 21
171 GLU ( 171 ) 1.58 20
17 LEU ( 17 ) 1.57 28
134 GLU ( 134 ) 1.55 80
# 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.
171 GLU ( 171 ) 0.35
111 MET ( 111 ) 0.38
# 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
25 SER ( 25 ) 0
36 LEU ( 36 ) 0
37 ASP ( 37 ) 0
56 LYS ( 56 ) 0
68 LYS ( 68 ) 0
93 ASN ( 93 ) 0
94 GLU ( 94 ) 0
120 ASN ( 120 ) 0
121 LYS ( 121 ) 0
127 PHE ( 127 ) 0
167 HIS ( 167 ) 0
168 VAL ( 168 ) 0
26 PHE ( 26 ) 1
27 GLN ( 27 ) 1
38 ARG ( 38 ) 1
55 LEU ( 55 ) 1
58 ARG ( 58 ) 1
77 MET ( 77 ) 1
122 ALA ( 122 ) 1
156 LYS ( 156 ) 1
24 GLY ( 24 ) 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 : -5.999
# 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 : 39.702
# 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.
16 GLN ( 16 )
19 ASN ( 19 )
86 ASN ( 86 )
93 ASN ( 93 )
158 GLN ( 158 )
# 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.06 HIS-E 0.59
3 HIS ( 3 ) HIS-H 0.07 HIS-E 0.60
4 HIS ( 4 ) HIS-H 0.08 HIS-E 0.58
5 HIS ( 5 ) HIS-H 0.04 HIS-E 0.58
6 HIS ( 6 ) HIS-H 0.06 HIS-E 0.58
7 HIS ( 7 ) HIS-H 0.05 HIS-E 0.60
18 HIS ( 18 ) HIS-H 0.05 HIS-E 0.59
72 HIS ( 72 ) HIS-H 0.05 HIS-D 0.56
139 HIS ( 139 ) HIS-H 0.05 HIS-E 0.59
167 HIS ( 167 ) HIS-H 0.04 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.
19 ASN ( 19 ) N
23 TYR ( 23 ) N
24 GLY ( 24 ) N
26 PHE ( 26 ) N
37 ASP ( 37 ) N
45 ALA ( 45 ) N
54 ARG ( 54 ) N
55 LEU ( 55 ) N
59 LEU ( 59 ) N
60 TYR ( 60 ) N
63 ILE ( 63 ) N
78 GLY ( 78 ) N
86 ASN ( 86 ) ND2
96 GLU ( 96 ) N
104 ARG ( 104 ) NE
119 ILE ( 119 ) N
121 LYS ( 121 ) N
122 ALA ( 122 ) N
125 ASP ( 125 ) N
130 TRP ( 130 ) N
132 PHE ( 132 ) N
154 LYS ( 154 ) N
156 LYS ( 156 ) 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 ) 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 : -3.797 (poor)
2nd generation packing quality : -3.056 (poor)
Ramachandran plot appearance : -4.104 (bad)
chi-1/chi-2 rotamer normality : -5.856 (bad)
Backbone conformation : -5.999 (bad)
RMS Z-scores, should be close to 1.0:
Bond lengths : 0.147 (tight)
Bond angles : 0.415 (tight)
Omega angle restraints : 0.021 (tight)
Side chain planarity : 0.038 (tight)
Improper dihedral distribution : 0.871
Inside/Outside distribution : 1.135
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).