BUSTER 2.X: A suite of programs for macromolecular structure refinement by maximum-likelihood and Bayesian methods
Copyright © 1992-2007 Global Phasing Ltd.
All rights reserved.
Authors: G. Bricogne, E. Blanc, M. Brandl, C. Flensburg, P. Keller, W. Paciorek, P. Roversi, O. Smart, C. Vonrhein, T. Womack

Please reference use of the software as:

Version: 2.8.0 for linux-ifort made on Jul 22 2009 at 02:57:53
Licenced to: Global Phasing Ltd (newton-old)


Parsing of the BUSTER cardfile

TITLE : autoBUSTER refinement Cell information (explanation) Symmetry information (explanation) Using CCP4 X-ray form factors. Missing structure chemical composition (explanation) Using ATOM block for chemical composition. Shell information

Parsing of the TNT cardfile


Reading observations file

List of rejected reflexions (explanation) Number of reflexions in input: 21874 Number of accepted reflexions: 21019 Number of rejected reflexions: 855 Reflexions in test set: 1669 (explanation) Reflexions in working set: 19350 (explanation)

1. Generating initial model

2.10 * Shannon sampling rate. NX NY NZ: 128 240 120 Effective Shannon rates are: 2.13 2.55 2.22 Allocated BUFFER of size: 14.282 Mb. Allocated BUFAUX of size: 14.282 Mb. Initial Geometry Information (explanation) No missing atoms declared Summary of model volumes (explanation) No missing atoms declared: no prior computed Computing bulk solvent model scattering NCSMASK PDB radius for solvent mask: 2.15 Angstrom running PDBSET (setting cell and spacegroup P1) ... done running NCSMASK [1] (creating mask, SMOOTH 1.23848) ... done running MAPMASK (converting mask to map) ... done Generation of babslv by theta-filter (explanation) Electron density in frg mask: 0.420 Å-1 (explanation) FOBS put on rough absolute scale by Wilson plot : K = 0.14910E+02 B = -7.219

Cycle 001: MaxLik scaling and refinement

MaxLik Scaling Cycle 001 (explanation) Wilson plot Statistics in resolution bins just after scaling (explanation) : 
     

   <Log-Likelihood> of null-hypothesis (explanation)

Structure Factor Amplitudes (explanation) Average Figure of Merit (explanation) Switching off refinement of overall temperature factor Added B = 2.4326 Å2 to all atoms in fragment Added B = 2.4326 Å2 to B_SOLV

gelly refinement module
 Developed by O. S. Smart, C. Flensburg,  W. Paciorek, C. Vonrhein, 
              M. Brandl, T. O. Womack and G. Bricogne
 Copyright © 2006, 2007 Global Phasing Ltd. 
 All rights reserved.
 gelly module compilation date Jul 22 2009 02:55:39
 
 gelly uses routines from TNT geometry. References: 
     Tronrud, D. E., Ten Eyck, L. F., Matthews, B. W. (1987)
     "An Efficient General-Purpose Least-Squares Refinement Program 
     for Macromolecular Structures". Acta Crystallogr A, 43:489-501.
     Tronrud, D. E. (1997). "The TNT Refinement Package", 
     in Macromolecular Crystallography, Part B, Eds Charlie Carter 
     and Robert Sweet, Methods in Enzymology, 277:306-319.
 wildcard matching function by Julian Robichaux http://www.nsftools.com

 Control flags set to:  -autoncs -screen \
                     100 -screen_sigma \
                     3.0 
 
 Setting up tnt routines including license check:

                   TNT Stereochemistry Program
 Version: 5.20.0.0
 Licence for Global Phasing Ltd (newton-old) is valid
 
 Using TNT geometry routines to load the tnt control file: gelly.tnt
 Using TNT routines to load initial coordinates from file fragment.cor
 Have copied information for    4217 atoms from TNT to gelly's store 
 
 The TNT WEIGHT cards used this run:
 -----WEIGHT       BOND     2.0
 -----WEIGHT       ANGLE    2.0
 -----WEIGHT       IMPROPER 0.0
 -----WEIGHT       TORSION  2.0
 -----WEIGHT       PSEUDO   0.0
 -----WEIGHT       TRIGONAL 2.0
 -----WEIGHT       PLANE    5.0
 -----WEIGHT       CONTACT  5.0
 -----WEIGHT       BCORREL  20.0
 -----WEIGHT       CHIRAL   5.0
 -----WEIGHT       NCS      50.0
 -----WEIGHT       IDEAL    4.0
 -----WEIGHT       RFACTOR  3.67
 Dictionary versions picked up from NOTE BUSTER_RESTRAINT_DICT cards as: 
 ------ protgeo_eh99.dat (V1.6) 20090612  STANDARD AMINO ACID DICTIONARY. BONDS AND ANGLES FROM ENGH AND HUBER EH99. OTHER VALUES BASED ON PREVIOUS TNT OR TAKEN FROM CCP4.
 ------ exoticaa.dat (V1.3.2.1) 20090706  COLLECTION OF NON-STANDARD AMINO ACIDS, MAINLY EH91 WITHOUT IDEAL DISTANCE INFO
 ------ nuclgeo.dat (V1.13) 20090617
 ------ bcorrel.dat (V1.15) 20080423
 ------ contact.dat (V1.15) 20070207
 ------ idealdist_contact.dat (V1.3.2.1) 20090706  IDEAL-DISTANCE CONTACT TERM DATA AS USED IN PROLSQ. VALUES USED HERE ARE BASED ON THE REFMAC 5.5 IMPLEMENTATION.
 ------ assume.dat (V1.9) 20090602
 Have picked up spacegroup name for pdb output as 'P 21 21 2 '
 Setting up Chain/Residue arrays
 
 Have command line specification -autoncs*
     Have found       1 NCS relations
     Pairwise matching residue number and type between chains: 
                    A        B    
          A         .      all   
          B                  .   
       Table key--                                                       
          all   : two chains have 100% match and same number of residues 
         multi  : ditto and they contribute to a multi chain ncs cluster 
          ***%  : two chains have >80% matching residue number/type      
           -    : chains appear to be unrelated                          
      chain A     have identified ncs to     1 other chains. So will weight each autoncs by (       2.000/#)=      2.0000
      chain B     have identified ncs to     1 other chains. So will weight each autoncs by (       2.000/#)=      2.0000
     For .Gelly cards produced see below - after pruning
 
 Setting up BUSTER_SET groups
 Have found    34 BUSTER_SET cards
 Now silently applying automatic/dictionary set definitions. For a verbose listing use '-verbose_set' flag.
 Applying user defined set definitions:
 =====NOTE BUSTER_SET ncsautoXcld = Water --> set 'ncsautoXcld' now has     0 atoms 
 =====NOTE BUSTER_SET ncsautoABset = Chain_A \ ncsautoXcld --> set 'ncsautoABset' now has  2109 atoms 
 Have defined    14 sets: 
      set: 'All' has   4217 atoms in    562 residues - automatically defined
      set: 'Empty' has      0 atoms in      0 residues - automatically defined
      set: 'ExoticAA' has      0 atoms in      0 residues - defined by dictionary file
      set: 'NotHET' has   4215 atoms in    560 residues - automatically defined
      set: 'Chain_A' has   2109 atoms in    281 residues - automatically defined
      set: 'Chain_B' has   2108 atoms in    281 residues - automatically defined
      set: 'Water' has      0 atoms in      0 residues - defined by dictionary file
      set: 'StdProtein' has   4217 atoms in    562 residues - defined by dictionary file
      set: 'Protein' has   4217 atoms in    562 residues - defined by dictionary file
      set: 'Back' has   2241 atoms in    562 residues - defined by dictionary file
      set: 'Side' has   1976 atoms in    532 residues - defined by dictionary file
      set: 'Other' has      0 atoms in      0 residues - defined by dictionary file
      set: 'ncsautoXcld' has      0 atoms in      0 residues - user specified
      set: 'ncsautoABset' has   2109 atoms in    281 residues - user specified
 
 Interpreting CONSTANT cards: 
 =====CONSTANT     OCC
 Have found and interpreted      1 CONSTANT cards in the TNT input. 
 Before any COMBINE - Number of atoms with fixed xyz =     0
 Before any COMBINE - Number of atoms with fixed B   =     0
 Before any COMBINE - Number of atoms with fixed OCC =  4217
 Before any COMBINE - Number of atoms with fixed BIJ =  4217
                            cf total number of atoms =  4217
 
 There are no COMBINE cards to be found in the TNT card input 
 
 There are no NOTE BUSTER_COMBINE_RES cards to be found in the TNT card input 
 RANGE allowed for B variables in refinement:   from      3.00000 to    300.00000
 No 'NOTE BUSTER_TARGET' cards or -target?? args found (used to read coordinates files for soft ncs to external)
 No 'NOTE BUSTER_OCCSUM' cards found - so not using sum of occupancy restraints
 No 'NOTE BUSTER_DISTANCE' cards found - so not setting up any additional atom-atom distance restraints
 No 'NOTE BUSTER_UTILANGLE' cards found - so not setting up any additional restraints on any angle
 No 'NOTE BUSTER_UTILTOR' cards found - so not setting up any additional sinusoidal restraints on torsions
 
 --- Have some similarity definitions 
     Similarity restraints on RMSD use superposition routines from: 
     Coutsias, E.A., Seok, C., Dill, K.A.(2004)
     "Using quaternions to calculate RMSD", J. Comput. Chem., 25:1849-1857 

 Have found    1 'NOTE BUSTER_SIM_DEFINE' cards 
     Maximum number of chains involved in any SIM group (excluding the template)=     1
 Interpreting card: NOTE BUSTER_SIM_DEFINE ncsautoAB ncsautoABset B
     Number of atoms in template set=  2109
     Number of related chains (other than template) in SIM set=    1
     List of chains to be restrained (excluding template)=   B   
     Pairwise RMSD, XYZ in Angstroms, for SIM group: ncsautoAB
                 A         B   
     A         -----     ----- 
     B         0.059     ----- 
     Overall RMSD,XYZ for the SIM group is     0.059 Angstroms
     Pairwise rms deviation of [individual B's minus their set/chain average] in Angs**2
                 A         B         Bavg
     A         -----     -----     32.71
     B          1.32     -----     33.26
     Overall RMSD,B   for the SIM group is     1.320 Angs**2
     The number of common atoms for the above fits:
                 A         B   
     A         -----     ----- 
     B          2108     ----- 
     (compared to the template set that has    2109 atoms in it)
     List of atoms that are absent in one or more chains for this group: 
                                    A    B   
                    N       7       +    -
     (key  +  atom present for this chain 
           -  atom absent  for this chain)
 transforms for each of the superpositions orthogonal matrix r11 r12 r13 r21 r22 r23 r31 r32 r33 tx ty tz (PDBSET syntax): 
TRANS_AB -0.999595 -0.028070 -0.004602  0.028101 -0.999582 -0.006834 -0.004408 -0.006960  0.999966   121.428   256.003   -35.416
 reverse transforms 
TRANS_BA -0.999595  0.028101 -0.004408 -0.028070 -0.999582 -0.006960 -0.004602 -0.006834  0.999966   114.029   259.058    37.723
 
 Setting up LSSR interatomic distance restraints
 Interpreting card: NOTE BUSTER_SIM_RESTRAIN_LSSR ncsautoAB 2.000
        setting up interatomic distance pair restraints with:
        Dmin=         0.00
        Dmax=         5.50
        sigma=        0.20
        Vmax=        10.00
        weight=       2.00
        number of distance pairs defined for this set is      24104
 Have not supplied either '-sim_swap_equiv' or '-sim_swap_equiv_plus' argument 
   so will not try swaping equivalent atoms to improve NCS and/or target agreement.
 
 NCS by LSSR outlier analysis
 -autoncs option automatic NCS prune routine -------- 
 Residues pruned out by gradient condition:
           A|23      A|24      A|25      A|27      A|83      A|84      A|85
           A|86      A|87      A|88      A|89      A|90      A|91      A|92
           A|93     A|181     A|184     A|246     A|247     A|288      B|23
           B|24      B|25      B|27      B|83      B|84      B|85      B|86
           B|87      B|88      B|89      B|90      B|91      B|92      B|93
          B|181     B|184     B|246     B|247     B|288  
 Residues pruned out by function/max condition:
       
 Total number of residues pruned out is       40 out of      561 participating
 Number of LSSR restraints pruned out is      1985 this is    8.2% of the total
  
 -autoncs .Gelly cards suitable for adaptation for manual definition/pruning
           (note auto pruning is done in a slightly different way so exact results may differ slightly)
 -autoncs option .Gelly cards also written to file  auto_gelly_cards.txt
NOTE BUSTER_SET ncsautoXcld = Water    # card added - produced by -autoncs
NOTE BUSTER_SET ncsautoXcld = ncsautoXcld + { A|23 A|24 A|25 A|27 A|83 A|84 }    # card added - produced by -autoncs
NOTE BUSTER_SET ncsautoXcld = ncsautoXcld + { A|85 A|86 A|87 A|88 A|89 A|90 }    # card added - produced by -autoncs
NOTE BUSTER_SET ncsautoXcld = ncsautoXcld + { A|91 A|92 A|93 A|181 A|184 }    # card added - produced by -autoncs
NOTE BUSTER_SET ncsautoXcld = ncsautoXcld + { A|246 A|247 A|288 }    # card added - produced by -autoncs
NOTE BUSTER_SET ncsautoABset = Chain_A \ ncsautoXcld    # card added - produced by -autoncs
NOTE BUSTER_SIM_DEFINE ncsautoAB ncsautoABset B    # card added - produced by -autoncs
NOTE BUSTER_SIM_RESTRAIN_LSSR ncsautoAB 2.000    # card added - produced by -autoncs

 LSSR function contribution before pruning       827.6
                             after pruning        62.2
 Graph of which residues pruned from LSSR: ncs_pruned.mtv
 Graph of LSSR function outliers ncs_function_outliers.mtv
 Graph of LSSR gradient outliers ncs_gradient_outliers.mtv
 LSSR - analysis of distance difference distribution for NCS
         DeltaD less than 0.5 Angs:  proportion of population=100.00% average DeltaD of subset     0.002 Angs
 Histogram of LSSR deltaD distribution ncs_histogram.mtv
 
 No 'NOTE BUSTER_TLS_SET' cards found - so no TLS constraints will be used
 Number of parameters to describe molecular conformation  =   16868
 Number of free parameters to be adjusted in refinement   =   16868
 No command line argument '-tnt_fn' so will use gelly to do geometry function calculations
      Getting symmetry operators from TNT. 
          gelly will classify symmetry using pdb-like convention: 
               SYMOP   SYMMETRY
              NNNMMM   OPERATOR
                1555   X,Y,Z
                2555   -X,-Y,Z
                3555   1/2-X,1/2+Y,-Z
                4555   1/2+X,1/2-Y,-Z
          where NNN -> operator number and MMM -> translation vector
 Ideal-distance contact term setup:
    Restraint function constants (normally values got from NOTE BUSTER_IDEAL_CONTACT_* cards)
    Sigma for general contacts      0.200 Angs
    Sigma for 1-4 contacts          0.200 Angs (contacts across a torsion)
    Sigma for HBOND pairs           0.200 Angs (contacts between HB_types A and D or B's)
    Contact distance increment for each atom in HBOND          -0.150 Angs 
    Contact distance increment for 1-4 with an O atom          -0.100 Angs 
    Contact distance increment for 1-4 with an N atom          -0.100 Angs 
    Contact distance increment for 1-4 with an C atom          -0.150 Angs 
    Contact distance increment for 1-4 with any other atom     -0.150 Angs 
    Have matched       21 residue-type records and will use      20 ccp4 energy atomtypes
    All atoms matched for ccp4 atom type - so ideal contact will be used for all contacts (rather than TNT short).
   Using gelly special position handling. (For TNT instead rerun with argument '-special_tnt')
       TERM      weight      number      TNT_fn    Gelly_fn        diff
       BOND       2.000        4311    2617.220    2617.220       0.000
       ANGL       2.000        5778    6304.804    6304.804       0.000
       IMPR       0.000           0       0.000       0.000       0.000
       TORS       2.000        1400    2877.756    2877.757       0.000
       TRIG       2.000          64      20.590      20.590       0.000
       PLAN       5.000         610     846.290     846.290       0.000
       CONT       5.000           0    5665.763       0.000    5665.763
       BCOR      20.000        4311     574.450     574.450       0.000

   Switching torsion functionial form to sinusoidal rather than split harmonic functional.
      if you want this not to happen then rerun with '-torsharm'
      Weight for sinusoidal torsion=     2.000
      Torsion angles split into peptide omega and other
         Number of peptide omega angles is       556
         Number of other torsion is              844
      Peptide omega angles are monitored rather than directly restrained (PLANE used instead).
 Now apply the fix to get correct restraints for cis PRO residues
   Have got       4 bonds records and       8 angles
   from GEOMETRY cards for the special PROC residue
   Apply the fix to cis PRO residue A|102 
   because omega torsion A|101:CA(PRO)=A|101:C(PRO)=A|102:N(PRO)=A|102:CA(PRO) is       -0.90 degrees
     Fix up bond A|102:CA=CB (PRO)
                    ideal from 1.531 to 1.533 actual= 1.533,  sigma from 0.020 to 0.018 Angs
     Fix up bond A|102:CB=CG (PRO)
                    ideal from 1.495 to 1.506 actual= 1.510,  sigma from 0.050 to 0.039 Angs
     Fix up bond A|102:CG=CD (PRO)
                    ideal from 1.502 to 1.512 actual= 1.513,  sigma from 0.033 to 0.027 Angs
     Fix up bond A|102:CD=N (PRO)
                    ideal from 1.474 to 1.474 actual= 1.478,  sigma from 0.014 to 0.014 Angs
     Fix up angle A|101:C(PRO)=A|102:N(PRO)=A|102:CA(PRO)
                    ideal from 119.3 to 127.0 actual= 128.9,  sigma from   1.5 to   2.4 Degs
     Fix up angle A|101:C(PRO)=A|102:N(PRO)=A|102:CD(PRO)
                    ideal from 128.4 to 120.6 actual= 119.7,  sigma from   2.1 to   2.2 Degs
     Fix up angle A|102:N=CA=CB (PRO)
                    ideal from 103.3 to 102.6 actual= 103.1,  sigma from   1.2 to   1.1 Degs
     Fix up angle A|102:CB=CA=C (PRO)
                    ideal from 111.7 to 112.0 actual= 113.3,  sigma from   2.1 to   2.5 Degs
     Fix up angle A|102:CA=CB=CG (PRO)
                    ideal from 104.8 to 104.0 actual= 105.2,  sigma from   1.9 to   1.9 Degs
     Fix up angle A|102:CB=CG=CD (PRO)
                    ideal from 106.5 to 105.4 actual= 107.6,  sigma from   3.9 to   2.3 Degs
     Fix up angle A|102:CG=CD=N (PRO)
                    ideal from 103.2 to 103.8 actual= 104.7,  sigma from   1.5 to   1.2 Degs
     Fix up angle A|102:CA=N=CD (PRO)
                    ideal from 111.7 to 111.5 actual= 111.4,  sigma from   1.4 to   1.4 Degs
   Apply the fix to cis PRO residue B|102 
   because omega torsion B|101:CA(PRO)=B|101:C(PRO)=B|102:N(PRO)=B|102:CA(PRO) is       -0.74 degrees
     Fix up bond B|102:CA=CB (PRO)
                    ideal from 1.531 to 1.533 actual= 1.534,  sigma from 0.020 to 0.018 Angs
     Fix up bond B|102:CB=CG (PRO)
                    ideal from 1.495 to 1.506 actual= 1.511,  sigma from 0.050 to 0.039 Angs
     Fix up bond B|102:CG=CD (PRO)
                    ideal from 1.502 to 1.512 actual= 1.514,  sigma from 0.033 to 0.027 Angs
     Fix up bond B|102:CD=N (PRO)
                    ideal from 1.474 to 1.474 actual= 1.478,  sigma from 0.014 to 0.014 Angs
     Fix up angle B|101:C(PRO)=B|102:N(PRO)=B|102:CA(PRO)
                    ideal from 119.3 to 127.0 actual= 128.8,  sigma from   1.5 to   2.4 Degs
     Fix up angle B|101:C(PRO)=B|102:N(PRO)=B|102:CD(PRO)
                    ideal from 128.4 to 120.6 actual= 119.7,  sigma from   2.1 to   2.2 Degs
     Fix up angle B|102:N=CA=CB (PRO)
                    ideal from 103.3 to 102.6 actual= 103.1,  sigma from   1.2 to   1.1 Degs
     Fix up angle B|102:CB=CA=C (PRO)
                    ideal from 111.7 to 112.0 actual= 113.3,  sigma from   2.1 to   2.5 Degs
     Fix up angle B|102:CA=CB=CG (PRO)
                    ideal from 104.8 to 104.0 actual= 105.2,  sigma from   1.9 to   1.9 Degs
     Fix up angle B|102:CB=CG=CD (PRO)
                    ideal from 106.5 to 105.4 actual= 107.6,  sigma from   3.9 to   2.3 Degs
     Fix up angle B|102:CG=CD=N (PRO)
                    ideal from 103.2 to 103.8 actual= 104.7,  sigma from   1.5 to   1.2 Degs
     Fix up angle B|102:CA=N=CD (PRO)
                    ideal from 111.7 to 111.5 actual= 111.5,  sigma from   1.4 to   1.4 Degs
   Weight for ideal-distance contact term=     4.000

   MaxLik Scaling Cycle 0002 (explanation)
 
 Initial functional value = 2.901124E+05 rms gradient=      36.3488
 X-ray contribution to function value = 2.493038E+05
 Geometry contribution (inc ncs)      =      40808.6
 Geometry function summary
     term                           number     weight       rms        GooF      function     fn/numb
     ================================================================================================
     BOND bond lengths (angs)       4311       2.000      0.0109       0.551      2616.36       0.607
     ANGL bond angles (degs)        5778       2.000        1.18       0.725      6079.28       1.052
     TORS split-harmonic tors        844       0.000       20.32       1.306         0.00       0.000
     SINTOR sinusoidal tors          844       2.000       20.32                  2015.37       2.388
     TRIG.onal planes (angs)          64       2.000      0.0080       0.401        20.59       0.322
     PLAN general plane (angs)       610       5.000      0.0105       0.527       846.29       1.387
     BCORrelations (angs*2)         4311      20.000       0.474       0.082       574.45       0.133
     CONT (bad contacts)               0       5.000                                 0.00
     IDEAL(ideal-dist contact)      4331       4.000                             28594.07
     CHIRAL (gelly semiharmon)       557       5.000                                 0.00
     SIM similarity restraint for NCS and/or TARGET                                 62.20
     OCCSUM restraint on sum of occupanies                                           0.00       as none defined
     DISTAN utility restraints interatomic distance                                  0.00       as none defined
     UTILANGLE utility angle restraints                                              0.00       as none defined
     UTILTOR utility restraints on torsion angles                                    0.00       as none defined
     Peptide omega torsion angles                           3.70 degs - monitor only
 
 WARNING: Have |delta/sigma| deviations > 5.0 sigma. Number of outliers for each term: 
 WARNING:         4 bond lengths. Worst is    10.7 sigs     1.80 Angs A|89:CA=C (GLU)
 WARNING:         1 bond angles.  Worst is     5.4 sigs   125.37 degs B|89:CA=CB=CG (GLU)
 WARNING:        16 idealD contacts. Worst    13.7 sigs     1.16 Angs A|89:CA=CA (GLU) symm: 1555=2675

 Weighted rms      0.0109         1.181  (temporary rms bond and angle for autobuster to pick up!)

 Initial Geometry Sanity Check

 Check for duplicate bonds: no duplicates found
 
 Check for duplicate or missing bond angles: all are OK
   There are    10 bad contacts where Delta/sigma is greater than     6.5
   Looking at each one in turn to see whether it is likely to indicate 
   a setup problem with geometry restraints:

        1 dist=   1.16 idealD=   3.90 sigma=   0.20 abs(Delta/sigma)=   13.69
          Between atoms A|89:CA (GLU) A|89:CA (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        2 dist=   1.10 idealD=   3.70 sigma=   0.20 abs(Delta/sigma)=   12.99
          Between atoms A|89:CA (GLU) A|89:C (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        3 dist=   1.48 idealD=   3.35 sigma=   0.20 abs(Delta/sigma)=    9.33
          Between atoms A|89:N (GLU) A|89:C (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        4 dist=   1.80 idealD=   3.55 sigma=   0.20 abs(Delta/sigma)=    8.77
          Between atoms A|89:CA (GLU) A|90:N (PRO)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        5 dist=   1.80 idealD=   3.55 sigma=   0.20 abs(Delta/sigma)=    8.73
          Between atoms A|89:N (GLU) A|89:CA (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        6 dist=   1.61 idealD=   3.20 sigma=   0.20 abs(Delta/sigma)=    7.97
          Between atoms A|89:N (GLU) A|90:N (PRO)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        7 dist=   1.99 idealD=   3.52 sigma=   0.20 abs(Delta/sigma)=    7.64
          Between atoms A|89:N (GLU) A|89:CB (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        8 dist=   2.13 idealD=   3.50 sigma=   0.20 abs(Delta/sigma)=    6.86
          Between atoms A|88:C (GLY) A|89:C (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

        9 dist=   2.51 idealD=   3.87 sigma=   0.20 abs(Delta/sigma)=    6.81
          Between atoms A|89:CA (GLU) A|89:CB (GLU)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.

       10 dist=   2.20 idealD=   3.52 sigma=   0.20 abs(Delta/sigma)=    6.62
          Between atoms A|89:CB (GLU) A|90:N (PRO)
          The contact is between two atoms via a symmetry operator,
          hopefully refinement will fix this.
 GEOMETRY CHECK OVERALL:  WARN

 TAO: Toolkit for Advanced Optimization. Reference: 
     Benson, S.J., McInnes, L.C., More, J., Sarich, J. (2005)
     "TAO User Manual (Revision 1.8)", Mathematics and 
     Computer Science Division, Argonne National Laboratory
     ANL/MCS-TM-242, http://www.mcs.anl.gov/tao
 
 Will write writing pdb file every 10 its during opt (overide by -keeppdb value)
 Will apply TAO optimization 
 Maximum number of iterations       40
 Convergence test on Grms       4.0000 (-glim)
 Limit on rms displacement to original position (rmsD2init): NONE test inactive (-dlim)

   MaxLik Scaling Cycle 0003 (explanation)
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init       Opttag_H
       0  2.901124E+05  2.493038E+05       40808.6       36.3488  0.000000E+00  0.000000E+00       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0004 (explanation)

   MaxLik Scaling Cycle 0005 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 001 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       1  2.899431E+05  2.493013E+05       40641.8       44.4385  8.767856E-04     0.0116109       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0006 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 002 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       2  2.897336E+05  2.492975E+05       40436.1       38.0133    0.00189257     0.0131602       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0007 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 003 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       3  2.894051E+05  2.492784E+05       40126.7       49.4711    0.00543372     0.0433422       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0008 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 004 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       4  2.892589E+05  2.492627E+05       39996.1       52.0392    0.00773707     0.0643794       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0009 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 005 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       5  2.891127E+05  2.492617E+05       39851.0       26.6888    0.00741581     0.0624999       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0010 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 006 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       6  2.889685E+05  2.492517E+05       39716.8       21.1519    0.00809443     0.0693602       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0011 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 007 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       7  2.888048E+05  2.492293E+05       39575.5       32.2533     0.0102141     0.0871238       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0012 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 008 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       8  2.886510E+05  2.491847E+05       39466.3       56.8004     0.0142409      0.124506       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0013 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 009 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
       9  2.884546E+05  2.491645E+05       39290.1       27.2903     0.0162354      0.144506       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0014 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 010 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      10  2.883047E+05  2.491418E+05       39162.9       23.4983     0.0180927      0.164955       Opttag_D
   Writing pdb file to gelly_output_iter000010.pdb
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0015 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 011 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      11  2.881555E+05  2.491106E+05       39044.9       30.5112     0.0202837      0.190596       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0016 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 012 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      12  2.880555E+05  2.490390E+05       39016.5       69.8180     0.0252390      0.247186       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0017 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 013 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      13  2.878420E+05  2.490268E+05       38815.3       27.4684     0.0257841      0.254659       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0018 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 014 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      14  2.877260E+05  2.490110E+05       38715.1       23.6547     0.0264845      0.263935       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0019 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 015 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      15  2.875925E+05  2.489767E+05       38615.8       28.8305     0.0284071      0.286561       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0020 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 016 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      16  2.873615E+05  2.489088E+05       38452.8       30.8788     0.0323398      0.331653       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0021 (explanation)

   MaxLik Scaling Cycle 0022 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 017 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      17  2.872360E+05  2.488453E+05       38390.6       57.0554     0.0359462      0.372657       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0023 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 018 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      18  2.869872E+05  2.487552E+05       38232.0       29.2559     0.0413505      0.432361       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0024 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 019 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      19  2.868161E+05  2.486972E+05       38119.0       20.6198     0.0444573      0.466914       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0025 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 020 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      20  2.866043E+05  2.486122E+05       37992.1       25.9156     0.0487073      0.513110       Opttag_D
   Writing pdb file to gelly_output_iter000020.pdb
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0026 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 021 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      21  2.864968E+05  2.485115E+05       37985.3       61.4226     0.0533587      0.560924       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0027 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 022 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      22  2.863049E+05  2.484764E+05       37828.5       28.9451     0.0548841      0.575436       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0028 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 023 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      23  2.861522E+05  2.484417E+05       37710.5       25.0305     0.0562476      0.586900       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0029 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 024 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      24  2.859741E+05  2.483640E+05       37610.0       30.1095     0.0597476      0.616772       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0030 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 025 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      25  2.857396E+05  2.482377E+05       37501.9       38.1979     0.0655430      0.661550       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0031 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 026 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      26  2.857359E+05  2.479715E+05       37764.5       87.4720     0.0782711      0.751989       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0032 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 027 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      27  2.854004E+05  2.480241E+05       37376.3       22.9263     0.0755637      0.729729       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0033 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 028 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      28  2.853196E+05  2.480212E+05       37298.4       22.1368     0.0755302      0.725658       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0034 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 029 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      29  2.851397E+05  2.479521E+05       37187.7       29.0154     0.0783051      0.747479       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0035 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 030 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      30  2.848753E+05  2.478068E+05       37068.4       34.0534     0.0841594      0.797116       Opttag_D
   Writing pdb file to gelly_output_iter000030.pdb
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0036 (explanation)

   MaxLik Scaling Cycle 0037 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 031 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      31  2.847201E+05  2.476508E+05       37069.4       46.7978     0.0907776      0.852962       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0038 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 032 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      32  2.845102E+05  2.475005E+05       37009.7       25.3381     0.0972346      0.906769       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0039 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 033 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      33  2.843584E+05  2.474049E+05       36953.5       28.0982      0.101113      0.937480       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0040 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 034 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      34  2.841786E+05  2.472767E+05       36901.9       36.7958      0.106558      0.980840       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0041 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 035 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      35  2.839318E+05  2.471264E+05       36805.4       35.0296      0.112811       1.02778       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0042 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 036 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      36  2.836683E+05  2.469402E+05       36728.1       39.7705      0.120210       1.09201       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0043 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 037 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      37  2.835105E+05  2.467801E+05       36730.4       50.4147      0.126735       1.15219       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0044 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 038 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      38  2.833520E+05  2.467817E+05       36570.3       29.3506      0.126428       1.14888       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0045 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 039 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      39  2.831803E+05  2.467396E+05       36440.7       29.2154      0.127741       1.16154       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 

   MaxLik Scaling Cycle 0046 (explanation)

   Wilson plot

   Statistics in resolution bins after refinement cycle 040 (explanation) : 
  
     
   Iters     Tot_Funct    Xray_Funct    Geom_Funct          Grms     rmsD2init     maxD2init
      40  2.829445E+05  2.466115E+05       36333.0       34.7915      0.132455       1.20584       Opttag_D
   Geometry Restraint Deviations SCREEN Information

 


 TLS information. Total number of groups:     0


 Group   Natoms    FixRB   FixTLS   Group Tag    Group Spec

 For final position, Geometry function summary
     term                           number     weight       rms        GooF      function     fn/numb
     ================================================================================================
     BOND bond lengths (angs)       4311       2.000      0.0100       0.507      2217.63       0.514
     ANGL bond angles (degs)        5778       2.000        1.00       0.613      4342.20       0.752
     TORS split-harmonic tors        844       0.000       19.63       1.262         0.00       0.000
     SINTOR sinusoidal tors          844       2.000       19.63                  1945.22       2.305
     TRIG.onal planes (angs)          64       2.000      0.0050       0.251         8.03       0.126
     PLAN general plane (angs)       610       5.000      0.0088       0.438       584.53       0.958
     BCORrelations (angs*2)         4311      20.000       0.468       0.081       559.48       0.130
     CONT (bad contacts)               0       5.000                                 0.00
     IDEAL(ideal-dist contact)      4291       4.000                             26635.78
     CHIRAL (gelly semiharmon)       557       5.000                                 4.30
     SIM similarity restraint for NCS and/or TARGET                                 35.79
     OCCSUM restraint on sum of occupanies                                           0.00       as none defined
     DISTAN utility restraints interatomic distance                                  0.00       as none defined
     UTILANGLE utility angle restraints                                              0.00       as none defined
     UTILTOR utility restraints on torsion angles                                    0.00       as none defined
     Peptide omega torsion angles                           3.84 degs - monitor only
 
 WARNING: Have |delta/sigma| deviations > 5.0 sigma. Number of outliers for each term: 
 WARNING:         3 bond lengths. Worst is    16.8 sigs     1.96 Angs A|89:CA=C (GLU)
 WARNING:        12 idealD contacts. Worst    15.2 sigs     0.86 Angs A|89:CA=CA (GLU) symm: 1555=2675

 Weighted rms      0.0100         0.998  (temporary rms bond and angle for autobuster to pick up!)
 Outputing geometry progress graph in plotmtv format to file:Geom_Res.mtv
 s/r gelly normal termination
Statistics and scale factors vs. cycle number  (explanation) : 
     
   Final Geometry Information (explanation)


2. Final results

Statistics and scale factors vs. cycle number  (explanation) : 
     

   Average Figure of Merit     (explanation)
   Structure Factor Amplitudes (explanation)



Final refined coordinates (explanation)

Final Fourier amplitudes and phases (explanation)

# normal termination Program stopped. closesok message: time spent in tree = 88.48 (user) 3.78 (sys)


Normal termination


 Resources usage summary:
     Time spent in client ...................   150.46 (user)        4.10 (system)
     Time spent in scripts (   104 calls) ...     1.18 (user)        0.34 (system)
     Time spent in server ...................     0.01 (user)        0.01 (system)

     Max resident memory used (kb) .......... 0
     Number of page faults .................. 0 (major)    240314 (minor)
     Number of context switches ............. 23978 (niced)    1403 (other)