
#######################################################
#                                                     #
#                        LGA                          #
#                  ---------------                    #
#                                                     #
#               Local-Global Alignment                #
#        A Method for Finding 3-D Similarities        #
#               in Protein Structures                 #
#                                                     #
#                  --------------- 06/2006            #
#                                                     #
#      Adam Zemla (adamz@llnl.gov)                    #
#      Lawrence Livermore National Laboratory, CA     #
#                                                     #
#######################################################

# Molecule1: number of CA atoms   36 (   36),  selected   36 , name T0335TS599_4-D1
# Molecule2: number of CA atoms   36 (  584),  selected   36 , name T0335_D1.pdb
# PARAMETERS: T0335TS599_4-D1.T0335_D1.pdb  -3  -ie  -o1  -sda  -d:4  
# FIXED Atom-Atom correspondence
# GDT and LCS analysis 

LCS - RMSD CUTOFF   5.00      length       segment         l_RMS    g_RMS
  LONGEST_CONTINUOUS_SEGMENT:    36         5 - 40          2.19     2.19
  LCS_AVERAGE:    100.00

LCS - RMSD CUTOFF   2.00      length       segment         l_RMS    g_RMS
  LONGEST_CONTINUOUS_SEGMENT:    34         6 - 39          1.99     2.21
  LONGEST_CONTINUOUS_SEGMENT:    34         7 - 40          2.00     2.25
  LCS_AVERAGE:     94.37

LCS - RMSD CUTOFF   1.00      length       segment         l_RMS    g_RMS
  LONGEST_CONTINUOUS_SEGMENT:    24        16 - 39          0.98     3.64
  LONGEST_CONTINUOUS_SEGMENT:    24        17 - 40          0.94     3.95
  LCS_AVERAGE:     64.89

LCS_GDT    MOLECULE-1    MOLECULE-2    LENGTH_OF_THE
LCS_GDT     RESIDUE       RESIDUE         CONTINUOUS
LCS_GDT   NAME NUMBER   NAME NUMBER          SEGMENT     GLOBAL DISTANCE TEST - GDT_DATA_COLUMNS                             NUMBER OF THE RESIDUES:   36
LCS_GDT     A       5     A       5     21   33   36     9   16   19   24   28   30   32   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K       6     K       6     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     I       7     I       7     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A       8     A       8     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     R       9     R       9     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     I      10     I      10     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     N      11     N      11     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      12     E      12     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     L      13     L      13     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A      14     A      14     22   34   36     9   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A      15     A      15     22   34   36     8   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K      16     K      16     24   34   36     8   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A      17     A      17     24   34   36     6   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K      18     K      18     24   34   36     5   16   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A      19     A      19     24   34   36     4   20   23   26   28   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     G      20     G      20     24   34   36     3   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     V      21     V      21     24   34   36     3   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     I      22     I      22     24   34   36     3   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     T      23     T      23     24   34   36     7   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      24     E      24     24   34   36    10   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      25     E      25     24   34   36    10   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      26     E      26     24   34   36    10   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K      27     K      27     24   34   36    10   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     A      28     A      28     24   34   36    10   20   23   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      29     E      29     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     Q      30     Q      30     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     Q      31     Q      31     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K      32     K      32     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     L      33     L      33     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     R      34     R      34     24   34   36    10   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     Q      35     Q      35     24   34   36     8   20   22   26   29   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     E      36     E      36     24   34   36     8   20   22   26   28   31   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     Y      37     Y      37     24   34   36     8   20   22   26   28   30   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     L      38     L      38     24   34   36     8   20   22   26   28   30   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     K      39     K      39     24   34   36     8   18   22   26   28   30   34   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_GDT     G      40     G      40     24   34   36     8   11   22   25   28   30   32   35   36   36   36   36   36   36   36   36   36   36   36   36 
LCS_AVERAGE  LCS_A:  86.42  (  64.89   94.37  100.00 )

GLOBAL_DISTANCE_TEST (the largest set of residues that can fit under specified DISTANCE_CUTOFF)
GDT DIST_CUTOFF  0.50   1.00   1.50   2.00   2.50   3.00   3.50   4.00   4.50   5.00   5.50   6.00   6.50   7.00   7.50   8.00   8.50   9.00   9.50  10.00
GDT NUMBER_CA     10     20     23     26     29     31     34     35     36     36     36     36     36     36     36     36     36     36     36     36 
GDT PERCENT_CA  27.78  55.56  63.89  72.22  80.56  86.11  94.44  97.22 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00
GDT RMS_LOCAL    0.35   7.10   1.01   1.22   1.56   1.71   1.99   2.08   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19
GDT RMS_ALL_CA   3.34   3.99   3.92   3.38   2.49   2.41   2.21   2.21   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19   2.19

#      Molecule1      Molecule2       DISTANCE
LGA    A       5      A       5          4.924
LGA    K       6      K       6          3.789
LGA    I       7      I       7          3.215
LGA    A       8      A       8          2.789
LGA    R       9      R       9          2.648
LGA    I      10      I      10          1.942
LGA    N      11      N      11          1.301
LGA    E      12      E      12          1.116
LGA    L      13      L      13          1.459
LGA    A      14      A      14          1.071
LGA    A      15      A      15          1.456
LGA    K      16      K      16          1.632
LGA    A      17      A      17          0.833
LGA    K      18      K      18          1.878
LGA    A      19      A      19          3.736
LGA    G      20      G      20          3.211
LGA    V      21      V      21          0.910
LGA    I      22      I      22          1.234
LGA    T      23      T      23          0.964
LGA    E      24      E      24          1.823
LGA    E      25      E      25          2.364
LGA    E      26      E      26          1.534
LGA    K      27      K      27          0.996
LGA    A      28      A      28          1.937
LGA    E      29      E      29          1.861
LGA    Q      30      Q      30          1.026
LGA    Q      31      Q      31          1.173
LGA    K      32      K      32          1.502
LGA    L      33      L      33          1.493
LGA    R      34      R      34          1.026
LGA    Q      35      Q      35          0.875
LGA    E      36      E      36          1.342
LGA    Y      37      Y      37          2.611
LGA    L      38      L      38          2.854
LGA    K      39      K      39          2.452
LGA    G      40      G      40          3.803

#              N1   N2   DIST      N    RMSD    GDT_TS     LGA_S     LGA_Q 
SUMMARY(GDT)   36   36    4.0     35    2.08    81.250    87.846     1.605

LGA_LOCAL      RMSD =  2.080  Number of atoms =   35  DIST =   4.00
LGA_ALL_ATOMS  RMSD =  2.209  Number of atoms =   36 
Std_ALL_ATOMS  RMSD =  2.189  (standard rmsd on all 36 assigned CA atoms) 

Unitary ROTATION matrix and the shift VECTOR superimpose MOLECULES  (1=>2)
  X_new =   0.213298 * X  +  -0.975585 * Y  +  -0.052321 * Z  +   5.182918
  Y_new =   0.749486 * X  +   0.197749 * Y  +  -0.631796 * Z  +  -2.985176
  Z_new =   0.626717 * X  +   0.095547 * Y  +   0.773367 * Z  + -18.693546 

Euler angles from the ROTATION matrix (XYZ convention, two solutions)
  Psi   =   0.122924   -3.018669  [ DEG:     7.0430   -172.9570 ]
  Theta =  -0.677333   -2.464260  [ DEG:   -38.8083   -141.1917 ]
  Phi   =   1.293534   -1.848059  [ DEG:    74.1140   -105.8860 ]
 
# END of job
REMARK  ---------------------------------------------------------- 
REMARK   Citing LGA: 
REMARK   Zemla A., LGA - a Method for Finding 3D Similarities in  
REMARK   Protein Structures, Nucleic Acids Research, 2003, V. 31, 
REMARK   No. 13, pp. 3370-3374. 
REMARK  ---------------------------------------------------------- 
REMARK   Superimposed MOLECULES (1=>2)   Output: 1 molecule      
REMARK     1: T0335TS599_4-D1                               
REMARK     2: T0335_D1.pdb                                  
REMARK   GDT and LCS analysis
REMARK   FIXED Atom-Atom correspondence
REMARK   LGA parameters: T0335TS599_4-D1.T0335_D1.pdb -3 -ie -o1 -sda -d:4 
REMARK   #CA        N1   N2  DIST    N   RMSD   LGA_S  RMSD_global
REMARK   SUMMARY:   36   36   4.0   35   2.08  87.846     2.19
REMARK  ---------------------------------------------------------- 
MOLECULE T0335TS599_4-D1
PFRMAT TS                                                                       
TARGET T0335                                                                    
MODEL  4                                                                        
PARENT N/A                                                                      
ATOM      5  CA  ALA     5      -1.922  -8.176  -9.491  1.00  0.00              
ATOM      6  CA  LYS     6      -0.040  -5.296  -7.892  1.00  0.00              
ATOM      7  CA  ILE     7      -2.980  -2.977  -8.426  1.00  0.00              
ATOM      8  CA  ALA     8      -5.366  -5.489  -6.836  1.00  0.00              
ATOM      9  CA  ARG     9      -3.061  -5.986  -3.865  1.00  0.00              
ATOM     10  CA  ILE    10      -2.697  -2.287  -3.243  1.00  0.00              
ATOM     11  CA  ASN    11      -6.420  -1.860  -3.346  1.00  0.00              
ATOM     12  CA  GLU    12      -6.964  -4.627  -0.866  1.00  0.00              
ATOM     13  CA  LEU    13      -4.404  -3.143   1.545  1.00  0.00              
ATOM     14  CA  ALA    14      -6.002   0.298   1.257  1.00  0.00              
ATOM     15  CA  ALA    15      -9.379  -1.212   2.151  1.00  0.00              
ATOM     16  CA  LYS    16      -7.818  -3.076   5.128  1.00  0.00              
ATOM     17  CA  ALA    17      -6.207   0.151   6.287  1.00  0.00              
ATOM     18  CA  LYS    18      -9.470   2.062   6.005  1.00  0.00              
ATOM     19  CA  ALA    19     -11.271  -0.516   8.108  1.00  0.00              
ATOM     20  CA  GLY    20      -8.610  -0.434  10.783  1.00  0.00              
ATOM     21  CA  VAL    21      -8.076   3.327  10.783  1.00  0.00              
ATOM     22  CA  ILE    22      -4.549   3.340   9.567  1.00  0.00              
ATOM     23  CA  THR    23      -2.044   5.370  11.482  1.00  0.00              
ATOM     24  CA  GLU    24       0.260   7.872   9.902  1.00  0.00              
ATOM     25  CA  GLU    25       3.118   5.303  10.090  1.00  0.00              
ATOM     26  CA  GLU    26       1.011   2.707   8.300  1.00  0.00              
ATOM     27  CA  LYS    27       0.168   5.271   5.567  1.00  0.00              
ATOM     28  CA  ALA    28       3.850   6.153   5.076  1.00  0.00              
ATOM     29  CA  GLU    29       4.780   2.454   4.799  1.00  0.00              
ATOM     30  CA  GLN    30       2.012   1.899   2.229  1.00  0.00              
ATOM     31  CA  GLN    31       3.195   4.811   0.156  1.00  0.00              
ATOM     32  CA  LYS    32       6.808   3.512   0.204  1.00  0.00              
ATOM     33  CA  LEU    33       5.689   0.022  -0.796  1.00  0.00              
ATOM     34  CA  ARG    34       3.625   1.356  -3.665  1.00  0.00              
ATOM     35  CA  GLN    35       6.514   3.424  -4.953  1.00  0.00              
ATOM     36  CA  GLU    36       8.906   0.446  -4.825  1.00  0.00              
ATOM     37  CA  TYR    37       6.471  -1.850  -6.639  1.00  0.00              
ATOM     38  CA  LEU    38       5.853   0.740  -9.292  1.00  0.00              
ATOM     39  CA  LYS    39       9.623   1.256  -9.831  1.00  0.00              
ATOM     40  CA  GLY    40      10.146  -2.509 -10.046  1.00  0.00              
TER                                                                             
END
