examples/unres/MINIM/outputs/1L2Y_ene.out_GB

   1 --------------------------------------------------------------------------------
   2                               FILE ASSIGNMENT
   3 --------------------------------------------------------------------------------
   4  Input file                      : 1L2Y_ene.inp
   5  Output file                     : 1L2Y_ene.out_GB
   6
   7  Sidechain potential file        :
   8  /users/adam/unres/PARAM/sc_GB_opt.1gab_3S_qclass5no310-shan2-sc-16-10-8k
   9  SCp potential file              : /users/adam/unres/PARAM/scp.parm
  10  Electrostatic potential file    : /users/adam/unres/PARAM/electr_631Gdp.parm
  11  Cumulant coefficient file       :
  12  /users/adam/unres/PARAM/fourier_opt.parm.1igd_hc_iter3_3
  13  Torsional parameter file        : /users/adam/unres/PARAM/torsion_631Gdp.parm
  14  Double torsional parameter file :
  15  /users/adam/unres/PARAM/torsion_double_631Gdp.parm
  16  SCCOR parameter file : /users/adam/unres/PARAM/rotcorr_AM1.parm
  17  Bond & inertia constant file    : /users/adam/unres/PARAM/bond.parm
  18  Bending parameter file          : /users/adam/unres/PARAM/thetaml.5parm
  19  Rotamer parameter file          : /users/adam/unres/PARAM/scgauss.parm
  20  Threading database              : /users/adam/unres/PARAM/patterns.cart
  21 --------------------------------------------------------------------------------
  22 ********************************************************************************
  23 United-residue force field calculation - serial job.
  24 ********************************************************************************
  25  ### LAST MODIFIED  11/03/09 1:19PM by czarek
  26  ++++ Compile info ++++
  27  Version MINI energy and minimization only
  28
  29 Potential is GB , exponents are   6 12
  30
  31 Disulfide bridge parameters:
  32 S-S bridge energy:      -5.50
  33 d0cm:      3.78 akcm:     15.10
  34 akth:     11.00 akct:     12.00
  35 v1ss:     -1.08 v2ss:      7.61 v3ss:     13.70
  36 RMSDBC =        3.0
  37 RMSDBC1 =        0.5
  38 RMSDBC1MAX =        1.5
  39 DRMS    =        0.1
  40 RMSDBCM =        3.0
  41 Time limit (min):     960.0
  42  RESCALE_MODE           2
  43 Library  routine used to diagonalize matrices.
  44
  45 Energy-term weights (unscaled):
  46
  47 WSCC=     1.352790 (SC-SC)
  48 WSCP=     1.593040 (SC-p)
  49 WELEC=    0.715340 (p-p electr)
  50 WVDWPP=   0.113710 (p-p VDW)
  51 WBOND=    1.000000 (stretching)
  52 WANG=     1.138730 (bending)
  53 WSCLOC=   0.162580 (SC local)
  54 WTOR=     1.985990 (torsional)
  55 WTORD=    1.570690 (double torsional)
  56 WSTRAIN=  1.000000 (SS bridges & dist. cnstr.)
  57 WEL_LOC=  0.160360 (multi-body 3-rd order)
  58 WCORR4=   0.428870 (multi-body 4th order)
  59 WCORR5=   0.000000 (multi-body 5th order)
  60 WCORR6=   0.000000 (multi-body 6th order)
  61 WSCCOR=   0.000000 (back-scloc correlation)
  62 WTURN3=   1.687220 (turns, 3rd order)
  63 WTURN4=   0.662300 (turns, 4th order)
  64 WTURN6=   0.000000 (turns, 6th order)
  65
  66 Hydrogen-bonding correlation between contact pairs of peptide groups
  67
  68 Scaling factor of 1,4 SC-p interactions:   0.400
  69 General scaling factor of SC-p interactions:   1.000
  70
  71 Energy-term weights (scaled):
  72
  73 WSCC=     1.352790 (SC-SC)
  74 WSCP=     1.593040 (SC-p)
  75 WELEC=    0.715340 (p-p electr)
  76 WVDWPP=   0.113710 (p-p VDW)
  77 WBOND=    1.000000 (stretching)
  78 WANG=     1.138730 (bending)
  79 WSCLOC=   0.162580 (SC local)
  80 WTOR=     1.985990 (torsional)
  81 WTORD=    1.570690 (double torsional)
  82 WSTRAIN=  1.000000 (SS bridges & dist. cnstr.)
  83 WEL_LOC=  0.160360 (multi-body 3-rd order)
  84 WCORR4=   0.428870 (multi-body 4th order)
  85 WCORR5=   0.000000 (multi-body 5th order)
  86 WCORR6=   0.000000 (multi-body 6th order)
  87 WSCCOR=   0.000000 (back-scloc correlatkion)
  88 WTURN3=   1.687220 (turns, 3rd order)
  89 WTURN4=   0.662300 (turns, 4th order)
  90 WTURN6=   0.000000 (turns, 6th order)
  91  Reference temperature for weights calculation:   300.000000000000
  92  Parameters of the SS-bond potential:
  93  D0CM   3.78000000000000       AKCM   15.1000000000000       AKTH
  94    11.0000000000000       AKCT   12.0000000000000
  95  V1SS  -1.08000000000000       V2SS   7.61000000000000       V3SS
  96    13.7000000000000
  97  EBR  -5.50000000000000
  98 PDB data will be read from file 1L2Y.pdb
  99  Nres:    21
 100 Backbone and SC coordinates as read from the PDB
 101   1 21  D    -9.841   4.399  -5.051       -9.841   4.399  -5.051
 102   2 14  ASN  -8.608   3.135  -1.618      -10.407   3.153  -2.437
 103   3  5  LEU  -4.923   4.002  -2.452       -4.618   6.091  -1.850
 104   4  8  TYR  -3.690   2.738   0.981       -1.959   3.143   3.797
 105   5  4  ILE  -5.857  -0.449   0.613       -7.484  -0.369   1.074
 106   6 13  GLN  -4.122  -1.167  -2.743       -5.089  -1.450  -4.853
 107   7  7  TRP  -0.716  -0.631  -0.993        1.727   0.440   1.450
 108   8  5  LEU  -1.641  -2.932   1.963       -2.244  -2.097   3.799
 109   9 19  LYS  -3.024  -5.791  -0.269       -3.820  -5.527  -3.146
 110  10 16  ASP   0.466  -6.016  -1.905        0.653  -5.125  -3.676
 111  11 10  GLY   2.060  -6.618   1.593        2.060  -6.618   1.593
 112  12 10  GLY   2.626  -2.967   2.723        2.626  -2.967   2.723
 113  13 20  PRO   6.333  -2.533   3.806        5.724  -2.372   5.058
 114  14 12  SER   7.049  -6.179   2.704        6.757  -6.938   3.675
 115  15 12  SER   6.389  -5.315  -1.015        5.245  -5.350  -1.546
 116  16 10  GLY   9.451  -3.116  -1.870        9.451  -3.116  -1.870
 117  17 18  ARG   7.289   0.084  -2.054        5.225  -1.826  -3.986
 118  18 20  PRO   6.782   3.088   0.345        7.458   3.741  -0.688
 119  19 20  PRO   3.287   4.031   1.686        4.025   4.206   2.856
 120  20 20  PRO   1.185   6.543  -0.353        0.358   5.421  -0.430
 121  21 12  SER   0.852  10.027   1.285        1.151  10.636   1.488
 122  22 21  D    -1.250  12.539  -0.754       -1.250  12.539  -0.754
 123 nsup= 20 nstart_sup=  2
 124  ITEL
 125            1          21           0
 126            2          14           1
 127            3           5           1
 128            4           8           1
 129            5           4           1
 130            6          13           1
 131            7           7           1
 132            8           5           1
 133            9          19           1
 134           10          16           1
 135           11          10           1
 136           12          10           2
 137           13          20           1
 138           14          12           1
 139           15          12           1
 140           16          10           1
 141           17          18           2
 142           18          20           2
 143           19          20           2
 144           20          20           1
 145           21          12           0
 146  ns=           0  iss:
 147 Boundaries in phi angle sampling:
 148 D      1    -180.0     180.0
 149 ASN    2    -180.0     180.0
 150 LEU    3    -180.0     180.0
 151 TYR    4    -180.0     180.0
 152 ILE    5    -180.0     180.0
 153 GLN    6    -180.0     180.0
 154 TRP    7    -180.0     180.0
 155 LEU    8    -180.0     180.0
 156 LYS    9    -180.0     180.0
 157 ASP   10    -180.0     180.0
 158 GLY   11    -180.0     180.0
 159 GLY   12    -180.0     180.0
 160 PRO   13    -180.0     180.0
 161 SER   14    -180.0     180.0
 162 SER   15    -180.0     180.0
 163 GLY   16    -180.0     180.0
 164 ARG   17    -180.0     180.0
 165 PRO   18    -180.0     180.0
 166 PRO   19    -180.0     180.0
 167 PRO   20    -180.0     180.0
 168 SER   21    -180.0     180.0
 169 D     22    -180.0     180.0
 170 nsup= 20
 171  nsup=          20  nstart_sup=           2  nstart_seq=           2
 172  NZ_START=           2  NZ_END=          21
 173  IZ_SC=           0
 174  Contact order:  0.000000000000000E+000
 175  Shifting contacts:           2           2
 176 Initial geometry will be read in.
 177
 178 Geometry of the virtual chain.
 179   Res           d     Theta     Gamma       Dsc     Alpha      Beta
 180 D     1     0.000     0.000     0.000     0.000     0.000     0.000
 181 ASN   2     3.800     0.000     0.000     1.684   102.356   -82.317
 182 LEU   3     3.800   104.929     0.000     1.939   120.092   -56.685
 183 TYR   4     3.800    87.695  -100.430     2.484   152.364    85.090
 184 ILE   5     3.800    81.904    72.426     1.776   134.976   -88.666
 185 GLN   6     3.800    82.466    66.026     2.240   122.343  -140.945
 186 TRP   7     3.800    84.466    51.434     2.605   152.178    38.024
 187 LEU   8     3.800    83.945    53.508     1.939   159.052   179.471
 188 LYS   9     3.800    85.396    60.569     2.541   100.558   -73.090
 189 ASP  10     3.800    91.449    44.177     1.709   139.961  -144.797
 190 GLY  11     3.800    94.602    67.604     0.000     0.000     0.000
 191 GLY  12     3.800   101.862   -72.814     0.000     0.000     0.000
 192 PRO  13     3.800   119.363   -61.916     1.345   117.453  -133.163
 193 SER  14     3.800    94.363   -75.894     1.150   137.025  -106.659
 194 SER  15     3.800    96.264    67.358     1.150   146.290  -130.305
 195 GLY  16     3.800   138.119   129.701     0.000     0.000     0.000
 196 ARG  17     3.800    96.299   -95.571     3.020    93.901  -102.747
 197 PRO  18     3.800   129.702    63.972     1.345   101.025  -111.641
 198 PRO  19     3.800   109.445   -74.504     1.345   113.043  -122.044
 199 PRO  20     3.800   106.349  -122.306     1.345    93.778  -102.374
 200 SER  21     3.800   106.042  -134.605     1.150   153.835  -143.303
 201 D    22     3.800   108.718    92.113     0.000     0.000     0.000
 202 Energy evaluation or minimization calculation.
 203
 204 ********************************************************************************
 205
 206  Time for energy evaluation  0.000000000000000E+000
 207
 208 Virtual-chain energies:
 209
 210 EVDW=     -4.954316E+01 WEIGHT=    1.352790D+00 (SC-SC)
 211 EVDW2=     5.228224E+01 WEIGHT=    1.593040D+00 (SC-p)
 212 EES=      -1.132072E+02 WEIGHT=    7.153400D-01 (p-p)
 213 EVDWPP=    2.375178E+01 WEIGHT=    1.137100D-01 (p-p VDW)
 214 ESTR=      6.779273E-27 WEIGHT=    1.000000D+00 (stretching)
 215 EBE=      -1.651995E+01 WEIGHT=    1.138730D+00 (bending)
 216 ESC=       7.185904E+01 WEIGHT=    1.625800D-01 (SC local)
 217 ETORS=     8.243370E+00 WEIGHT=    1.985990D+00 (torsional)
 218 ETORSD=    2.390466E+00 WEIGHT=    1.570690D+00 (double torsional)
 219 EHBP=      0.000000E+00 WEIGHT=    1.000000D+00 (SS bridges & dist. cnstr.)
 220 ECORR4=   -8.338619E+01 WEIGHT=    4.288700D-01 (multi-body)
 221 ECORR5=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
 222 ECORR6=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
 223 EELLO=    -5.753258E+00 WEIGHT=    1.603600D-01 (electrostatic-local)
 224 ETURN3=    1.901845E+01 WEIGHT=    1.687220D+00 (turns, 3rd order)
 225 ETURN4=   -7.304754E+00 WEIGHT=    6.623000D-01 (turns, 4th order)
 226 ETURN6=    0.000000E+00 WEIGHT=    0.000000D+00 (turns, 6th order)
 227 ESCCOR=    0.000000E+00 WEIGHT=    0.000000D+00 (backbone-rotamer corr)
 228 EDIHC=     0.000000E+00 (dihedral angle constraints)
 229 ESS=       0.000000E+00 (disulfide-bridge intrinsic energy)
 230 UCONST=     0.000000E+00 (Constraint energy)
 231 ETOT=     -5.845169E+01 (total)
 232
 233
 234 ***** Computation time:    0 hours  0 minutes  0 seconds *****