capra-eval


Summary: capra-eval is a handy utility for comparing a model built by TEXTALTM with a true structure (if known). Essentially, it pairs up nearest atoms between the two structures, reports the matches it makes (and those that remain unmatched), and then computes an RMSD (root-mean-sequare distance) between the matched atoms, enabling the goodness of a model to be assessed.

Comments: If an atom in one model does not have a unique neighbor in the other model within a cutoff of 3.0Å, then it is left unassigned and excluded from the RMSD calculation. This makes the RMSD value less sensitive to being skewed by "over-hanging" or un-matched regions. Obviously, capra-eval only works if you have a model of the true structure, for example, when running tests on TEXTALTM with a protein whose structure is already known.

Usage: capra-eval -a|-c <textal-model> <true-structure>

Command-line arguments:

  1. (Mode) -a | -c (-a means compute RMSD over all atom types; -c means just match up C-alpha's)
  2. (Input) textal model (PDB file) - generated by TEXTALTM
  3. (Input) true structure (PDB file) - supplied by user, if known

Explanation of Output:

The output of capra-eval contains the following information:

Example:
57 tropical> trim-chains ../capra/czra/czra-chains.pdb ../capra/czra.pdb 
# command line:
# capra-eval czra1-imp5.pdb czra1.pdb 
#
# Date:Fri Aug 17 13:04:49 2001
#
# Compilation options:
# COMPRESSED_MAPS
#
# !textal global parameters
# 0 ! mode (see textal.h for settings)
# 1 ! Seed for RNG.  If non-zero, seed with this, otherwise use time.
# 4 ! number of radii
# 3.000000 ! radius 1
# 4.000000 ! radius 2
# 5.000000 ! radius 3
# 6.000000 ! radius 4
# 0 ! should we peakmatch (0/1)?
# 9 ! integer, number of peaks to use in peak-matching
# 5.000000 ! float, radius to use for density correlation calculations
# 0.700000 ! float, CC, good correlation cutoff (0.7)
# 400 ! int, K: number of top features to calc correlations on
# 0 ! int, connected should we use only connected density(0/1)?
# 1 !use sidechain axis method for CC calculations?
# 0 !use error_correcting CC?
# !end textal global parameters
#
CA 14 in CAPRA model matches CA 49 in true struct, dist=0.143221
LEU LEU $
CA 69 in CAPRA model matches CA 17 in true struct, dist=0.148825
THR THR $
CA 70 in CAPRA model matches CA 18 in true struct, dist=0.155348
SER GLU $
CA 51 in CAPRA model matches CA 90 in true struct, dist=0.166043
ARG MET $
CA 71 in CAPRA model matches CA 19 in true struct, dist=0.173582
ILE ILE $
CA 67 in CAPRA model matches CA 15 in true struct, dist=0.180827
ILE ARG $
CA 18 in CAPRA model matches CA 53 in true struct, dist=0.1826
GLU GLN $
CA 68 in CAPRA model matches CA 16 in true struct, dist=0.191828
THR VAL $
CA 25 in CAPRA model matches CA 60 in true struct, dist=0.191938
ASN LEU $
CA 82 in CAPRA model matches CA 29 in true struct, dist=0.195699
THR ILE $
CA 34 in CAPRA model matches CA 69 in true struct, dist=0.20645
THR VAL $
CA 52 in CAPRA model matches CA 91 in true struct, dist=0.233409
VAL LEU $
CA 73 in CAPRA model matches CA 21 in true struct, dist=0.247099
LYS LYS $
CA 20 in CAPRA model matches CA 55 in true struct, dist=0.248034
LEU ASN $
CA 72 in CAPRA model matches CA 20 in true struct, dist=0.267033
ALA PHE $
CA 64 in CAPRA model matches CA 12 in true struct, dist=0.271284
SER THR $
CA 23 in CAPRA model matches CA 58 in true struct, dist=0.276682
PHE HIS $
CA 35 in CAPRA model matches CA 70 in true struct, dist=0.28461
VAL LYS $
CA 22 in CAPRA model matches CA 57 in true struct, dist=0.285651
THR SER $
CA 10 in CAPRA model matches CA 45 in true struct, dist=0.294118
VAL ILE $
CA 21 in CAPRA model matches CA 56 in true struct, dist=0.296282
VAL VAL $
CA 74 in CAPRA model matches CA 22 in true struct, dist=0.302792
ALA ALA $
CA 86 in CAPRA model matches CA 33 in true struct, dist=0.305709
VAL GLU $
CA 84 in CAPRA model matches CA 31 in true struct, dist=0.310475
VAL ILE $
CA 53 in CAPRA model matches CA 92 in true struct, dist=0.310539
ALA LYS $
CA 8 in CAPRA model matches CA 43 in true struct, dist=0.31382
ALA GLY $
CA 32 in CAPRA model matches CA 67 in true struct, dist=0.316397
LEU HIS $
CA 83 in CAPRA model matches CA 30 in true struct, dist=0.320017
LYS ARG $
CA 24 in CAPRA model matches CA 59 in true struct, dist=0.323953
VAL GLN $
CA 19 in CAPRA model matches CA 54 in true struct, dist=0.325659
SER SER $
CA 81 in CAPRA model matches CA 28 in true struct, dist=0.327022
SER ARG $
CA 33 in CAPRA model matches CA 68 in true struct, dist=0.344496
ILE LEU $
CA 75 in CAPRA model matches CA 23 in true struct, dist=0.353048
LEU LEU $
CA 56 in CAPRA model matches CA 95 in true struct, dist=0.359246
VAL ILE $
CA 85 in CAPRA model matches CA 32 in true struct, dist=0.363176
GLN MET $
CA 13 in CAPRA model matches CA 48 in true struct, dist=0.363553
VAL GLN $
CA 2 in CAPRA model matches CA 37 in true struct, dist=0.374253
ILE VAL $
CA 31 in CAPRA model matches CA 66 in true struct, dist=0.383704
ILE VAL $
CA 27 in CAPRA model matches CA 62 in true struct, dist=0.384827
SER LEU $
CA 9 in CAPRA model matches CA 44 in true struct, dist=0.395598
SER HIS $
CA 26 in CAPRA model matches CA 61 in true struct, dist=0.420359
ALA LYS $
CA 28 in CAPRA model matches CA 63 in true struct, dist=0.43851
VAL LEU $
CA 37 in CAPRA model matches CA 72 in true struct, dist=0.451946
ILE LYS $
CA 3 in CAPRA model matches CA 38 in true struct, dist=0.461032
VAL SER $
CA 54 in CAPRA model matches CA 93 in true struct, dist=0.463501
MET GLN $
CA 63 in CAPRA model matches CA 11 in true struct, dist=0.473032
LEU ASP $
CA 50 in CAPRA model matches CA 89 in true struct, dist=0.476538
VAL THR $
CA 1 in CAPRA model matches CA 36 in true struct, dist=0.483934
SER SER $
CA 87 in CAPRA model matches CA 34 in true struct, dist=0.489684
VAL LEU $
CA 7 in CAPRA model matches CA 42 in true struct, dist=0.494246
ALA VAL $
CA 36 in CAPRA model matches CA 71 in true struct, dist=0.53035
ALA ALA $
CA 66 in CAPRA model matches CA 14 in true struct, dist=0.530894
ARG GLU $
CA 62 in CAPRA model matches CA 10 in true struct, dist=0.537468
ALA THR $
CA 80 in CAPRA model matches CA 27 in true struct, dist=0.550415
ALA ASN $
CA 16 in CAPRA model matches CA 51 in true struct, dist=0.556051
ILE LEU $
CA 44 in CAPRA model matches CA 83 in true struct, dist=0.558243
LEU ASP $
CA 61 in CAPRA model matches CA 9 in true struct, dist=0.564629
ALA ASN $
CA 65 in CAPRA model matches CA 13 in true struct, dist=0.606639
VAL LEU $
CA 17 in CAPRA model matches CA 52 in true struct, dist=0.627004
SER SER $
CA 55 in CAPRA model matches CA 94 in true struct, dist=0.627106
LEU ALA $
CA 42 in CAPRA model matches CA 81 in true struct, dist=0.653717
SER SER $
CA 60 in CAPRA model matches CA 99 in true struct, dist=0.663103
ALA ASN $
CA 12 in CAPRA model matches CA 47 in true struct, dist=0.704399
ASP HIS $
CA 43 in CAPRA model matches CA 82 in true struct, dist=0.730131
GLU LEU $
CA 79 in CAPRA model matches CA 26 in true struct, dist=0.736839
LEU TYR $
CA 29 in CAPRA model matches CA 64 in true struct, dist=0.779634
LEU LYS $
CA 11 in CAPRA model matches CA 46 in true struct, dist=0.784573
VAL SER $
CA 38 in CAPRA model matches CA 73 in true struct, dist=0.799915
ALA ARG $
CA 30 in CAPRA model matches CA 65 in true struct, dist=0.845913
SER SER $
CA 45 in CAPRA model matches CA 84 in true struct, dist=0.846603
GLY ASP $
CA 41 in CAPRA model matches CA 80 in true struct, dist=0.863744
LYS TYR $
CA 15 in CAPRA model matches CA 50 in true struct, dist=0.894648
ALA ASN $
CA 88 in CAPRA model matches CA 35 in true struct, dist=0.915986
VAL LEU $
CA 57 in CAPRA model matches CA 96 in true struct, dist=0.946286
SER HIS $
CA 46 in CAPRA model matches CA 85 in true struct, dist=0.960305
VAL ILE $
CA 47 in CAPRA model matches CA 86 in true struct, dist=1.04083
ALA HIS $
CA 78 in CAPRA model matches CA 25 in true struct, dist=1.13521
LEU ASP $
CA 49 in CAPRA model matches CA 88 in true struct, dist=1.22177
ALA ALA $
CA 6 in CAPRA model matches CA 41 in true struct, dist=1.31865
VAL SER $
CA 4 in CAPRA model matches CA 39 in true struct, dist=1.32312
VAL GLU $
CA 48 in CAPRA model matches CA 87 in true struct, dist=1.53654
ALA VAL $
CA 5 in CAPRA model matches CA 40 in true struct, dist=1.65693
MET ALA $
CA 76 in CAPRA model matches CA 24 in true struct, dist=1.78613
SER GLY $
CA 59 in CAPRA model matches CA 98 in true struct, dist=1.90765
PRO ALA $
CA 58 in CAPRA model matches CA 97 in true struct, dist=2.97927
GLN HIS $
CA 40 in CAPRA model matches CA 79 in true struct, dist=5.78428
GLY ILE $
CA 39 in CAPRA model matches CA 78 in true struct, dist=5.87647
GLY MET $
rms=1.14881 (n=87)
CAPRA atoms not matched...
0 77 
atoms in true struct. not matched...
74 75 76 77 100 101 102 103