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

!!!   License for LGA program has expired  !!!
#   For updated version of the LGA program   #
# please contact Adam Zemla:  adamz@llnl.gov #

# Molecule1: number of CA atoms    5 (   35),  selected    5 , name one
# Molecule2: number of CA atoms   63 (  471),  selected    5 , name two
# PARAMETERS: input.pdb  -ie  -3  -d:4  -sda  -swap  -ch1:A  -ch2:A  
# FIXED Atom-Atom correspondence
# GDT and LCS analysis 

LCS - RMSD CUTOFF   5.00      length       segment         l_RMS    g_RMS
  LONGEST_CONTINUOUS_SEGMENT:     5      39_A - 43_A        2.03     2.03
  LCS_AVERAGE:      7.94

LCS - RMSD CUTOFF   2.00      length       segment         l_RMS    g_RMS
  LONGEST_CONTINUOUS_SEGMENT:     4      39_A - 42_A        1.58     2.80
  LONGEST_CONTINUOUS_SEGMENT:     4      40_A - 43_A        1.77     2.65
  LCS_AVERAGE:      6.35

LCS - RMSD CUTOFF   1.00      length       segment         l_RMS    g_RMS
# WARNING! The change of the parameter DIST cutoff may give you better result.

  LONGEST_CONTINUOUS_SEGMENT:     3      39_A - 41_A        0.55     3.92
  LONGEST_CONTINUOUS_SEGMENT:     3      41_A - 43_A        0.04     5.24
  LCS_AVERAGE:      4.76

LCS_GDT    MOLECULE-1    MOLECULE-2     LCS_DETAILS     GDT_DETAILS                                                    TOTAL NUMBER OF RESIDUE PAIRS:    5
LCS_GDT     RESIDUE       RESIDUE       SEGMENT_SIZE    GLOBAL DISTANCE TEST COLUMNS: number of residues under the threshold assigned to each residue pair
LCS_GDT   NAME NUMBER   NAME NUMBER    1.0  2.0  5.0    0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5  6.0  6.5  7.0  7.5  8.0  8.5  9.0  9.5 10.0
LCS_GDT     G    39_A     G    39_A      3    4    5      1    3    3    4    4    5    5    5    5    5    5    5    5    5    5    5    5    5    5    5 
LCS_GDT     T    40_A     T    40_A      3    4    5      1    3    3    4    4    5    5    5    5    5    5    5    5    5    5    5    5    5    5    5 
LCS_GDT     Y    41_A     Y    41_A      3    4    5      3    3    3    4    4    5    5    5    5    5    5    5    5    5    5    5    5    5    5    5 
LCS_GDT     G    42_A     G    42_A      3    4    5      3    3    3    4    4    5    5    5    5    5    5    5    5    5    5    5    5    5    5    5 
LCS_GDT     N    43_A     N    43_A      3    4    5      3    3    3    4    4    5    5    5    5    5    5    5    5    5    5    5    5    5    5    5 
LCS_AVERAGE  LCS_A:   6.35  (   4.76    6.35    7.94 )

GLOBAL_DISTANCE_TEST (summary information about detected largest sets of residues (represented by selected AToms) that can fit under specified thresholds)
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_AT      3      3      3      4      4      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5 
GDT PERCENT_AT   4.76   4.76   4.76   6.35   6.35   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94   7.94
GDT RMS_LOCAL    0.04   0.04   0.04   1.47   1.47   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03
GDT RMS_ALL_AT   5.24   5.24   5.24   2.23   2.23   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03   2.03

# Checking swapping
#   possible swapping detected:  Y    41_A      Y    41_A

#      Molecule1      Molecule2  DISTANCE    Mis    MC     All    Dist_max   GDC_mc  GDC_all
LGA    G    39_A      G    39_A     1.701     0    0.356   0.356     3.805   61.548   61.548
LGA    T    40_A      T    40_A     0.698     0    0.176   0.210     2.455   85.952   79.184
LGA    Y    41_A      Y    41_A     2.906     0    0.505   1.279    12.641   61.190   27.659
LGA    G    42_A      G    42_A     1.865     0    0.234   0.234     2.354   70.952   70.952
LGA    N    43_A      N    43_A     2.298     0    0.594   1.350     5.576   63.810   50.119

# RMSD_GDC results:       CA      MC common percent     ALL common percent   GDC_mc  GDC_all
NUMBER_OF_ATOMS_AA:        5      20     20  100.00      35     35  100.00                63
SUMMARY(RMSD_GDC):     2.029          2.100                  2.776            5.452    4.595

#CA            N1   N2   DIST      N    RMSD    GDT_TS    LGA_S3     LGA_Q 
SUMMARY(GDT)    5   63    4.0      5    2.03     6.746     6.837     0.235

LGA_LOCAL      RMSD:   2.029  Number of atoms:    5  under DIST:   4.00
LGA_ASGN_ATOMS RMSD:   2.029  Number of assigned atoms:    5 
Std_ASGN_ATOMS RMSD:   2.029  Standard rmsd on all 5 assigned CA atoms 

Unitary ROTATION matrix and the SHIFT vector superimpose molecules  (1=>2)
  X_new =  -0.636461 * X  +  -0.453314 * Y  +   0.624039 * Z  + -39.814526
  Y_new =  -0.131038 * X  +   0.860853 * Y  +   0.491693 * Z  +  -9.593656
  Z_new =  -0.760097 * X  +   0.231170 * Y  +  -0.607300 * Z  +  41.140179 

Euler angles from the ROTATION matrix. Conventions XYZ and ZXZ:
           Phi     Theta       Psi   [DEG:       Phi     Theta       Psi ]
XYZ: -2.938544  0.863462  2.777876   [DEG: -168.3662   49.4727  159.1606 ]
ZXZ:  2.238128  2.223454 -1.275552   [DEG:  128.2353  127.3945  -73.0838 ]
 
# END of job
