Merge branch 'prerelease-3.2.1' of mmka:unres into prerelease-3.2.1
[unres.git] / examples / unres / MINIM / 1L2Y_min-rand.out_GB
1 --------------------------------------------------------------------------------
2                               FILE ASSIGNMENT
3 --------------------------------------------------------------------------------
4  Input file                      : 1L2Y_min-rand.inp
5  Output file                     : 1L2Y_min-rand.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 ********************************************************************************
46                     Options in energy minimization:
47 ********************************************************************************
48 MaxMin: 2000 MaxFun: 5000MinMin: 2000 MinFun: 2000 TolF:   1.00000E-02 RTolF:   1.00000E-04
49
50 Energy-term weights (unscaled):
51
52 WSCC=     1.352790 (SC-SC)
53 WSCP=     1.593040 (SC-p)
54 WELEC=    0.715340 (p-p electr)
55 WVDWPP=   0.113710 (p-p VDW)
56 WBOND=    1.000000 (stretching)
57 WANG=     1.138730 (bending)
58 WSCLOC=   0.162580 (SC local)
59 WTOR=     1.985990 (torsional)
60 WTORD=    1.570690 (double torsional)
61 WSTRAIN=  1.000000 (SS bridges & dist. cnstr.)
62 WEL_LOC=  0.160360 (multi-body 3-rd order)
63 WCORR4=   0.428870 (multi-body 4th order)
64 WCORR5=   0.000000 (multi-body 5th order)
65 WCORR6=   0.000000 (multi-body 6th order)
66 WSCCOR=   0.000000 (back-scloc correlation)
67 WTURN3=   1.687220 (turns, 3rd order)
68 WTURN4=   0.662300 (turns, 4th order)
69 WTURN6=   0.000000 (turns, 6th order)
70
71 Hydrogen-bonding correlation between contact pairs of peptide groups
72
73 Scaling factor of 1,4 SC-p interactions:   0.400
74 General scaling factor of SC-p interactions:   1.000
75
76 Energy-term weights (scaled):
77
78 WSCC=     1.352790 (SC-SC)
79 WSCP=     1.593040 (SC-p)
80 WELEC=    0.715340 (p-p electr)
81 WVDWPP=   0.113710 (p-p VDW)
82 WBOND=    1.000000 (stretching)
83 WANG=     1.138730 (bending)
84 WSCLOC=   0.162580 (SC local)
85 WTOR=     1.985990 (torsional)
86 WTORD=    1.570690 (double torsional)
87 WSTRAIN=  1.000000 (SS bridges & dist. cnstr.)
88 WEL_LOC=  0.160360 (multi-body 3-rd order)
89 WCORR4=   0.428870 (multi-body 4th order)
90 WCORR5=   0.000000 (multi-body 5th order)
91 WCORR6=   0.000000 (multi-body 6th order)
92 WSCCOR=   0.000000 (back-scloc correlatkion)
93 WTURN3=   1.687220 (turns, 3rd order)
94 WTURN4=   0.662300 (turns, 4th order)
95 WTURN6=   0.000000 (turns, 6th order)
96  Reference temperature for weights calculation:   300.000000000000     
97  Parameters of the SS-bond potential:
98  D0CM   3.78000000000000       AKCM   15.1000000000000       AKTH
99    11.0000000000000       AKCT   12.0000000000000     
100  V1SS  -1.08000000000000       V2SS   7.61000000000000       V3SS
101    13.7000000000000     
102  EBR  -5.50000000000000     
103 PDB data will be read from file 1L2Y.pdb
104  Nres:    21
105 Backbone and SC coordinates as read from the PDB
106   1 21  D    -9.841   4.399  -5.051       -9.841   4.399  -5.051
107   2 14  ASN  -8.608   3.135  -1.618      -10.407   3.153  -2.437
108   3  5  LEU  -4.923   4.002  -2.452       -4.618   6.091  -1.850
109   4  8  TYR  -3.690   2.738   0.981       -1.959   3.143   3.797
110   5  4  ILE  -5.857  -0.449   0.613       -7.484  -0.369   1.074
111   6 13  GLN  -4.122  -1.167  -2.743       -5.089  -1.450  -4.853
112   7  7  TRP  -0.716  -0.631  -0.993        1.727   0.440   1.450
113   8  5  LEU  -1.641  -2.932   1.963       -2.244  -2.097   3.799
114   9 19  LYS  -3.024  -5.791  -0.269       -3.820  -5.527  -3.146
115  10 16  ASP   0.466  -6.016  -1.905        0.653  -5.125  -3.676
116  11 10  GLY   2.060  -6.618   1.593        2.060  -6.618   1.593
117  12 10  GLY   2.626  -2.967   2.723        2.626  -2.967   2.723
118  13 20  PRO   6.333  -2.533   3.806        5.724  -2.372   5.058
119  14 12  SER   7.049  -6.179   2.704        6.757  -6.938   3.675
120  15 12  SER   6.389  -5.315  -1.015        5.245  -5.350  -1.546
121  16 10  GLY   9.451  -3.116  -1.870        9.451  -3.116  -1.870
122  17 18  ARG   7.289   0.084  -2.054        5.225  -1.826  -3.986
123  18 20  PRO   6.782   3.088   0.345        7.458   3.741  -0.688
124  19 20  PRO   3.287   4.031   1.686        4.025   4.206   2.856
125  20 20  PRO   1.185   6.543  -0.353        0.358   5.421  -0.430
126  21 12  SER   0.852  10.027   1.285        1.151  10.636   1.488
127  22 21  D    -1.250  12.539  -0.754       -1.250  12.539  -0.754
128 nsup= 20 nstart_sup=  2
129  ITEL
130            1          21           0
131            2          14           1
132            3           5           1
133            4           8           1
134            5           4           1
135            6          13           1
136            7           7           1
137            8           5           1
138            9          19           1
139           10          16           1
140           11          10           1
141           12          10           2
142           13          20           1
143           14          12           1
144           15          12           1
145           16          10           1
146           17          18           2
147           18          20           2
148           19          20           2
149           20          20           1
150           21          12           0
151  ns=           0  iss:
152 Boundaries in phi angle sampling:
153 D      1    -180.0     180.0
154 ASN    2    -180.0     180.0
155 LEU    3    -180.0     180.0
156 TYR    4    -180.0     180.0
157 ILE    5    -180.0     180.0
158 GLN    6    -180.0     180.0
159 TRP    7    -180.0     180.0
160 LEU    8    -180.0     180.0
161 LYS    9    -180.0     180.0
162 ASP   10    -180.0     180.0
163 GLY   11    -180.0     180.0
164 GLY   12    -180.0     180.0
165 PRO   13    -180.0     180.0
166 SER   14    -180.0     180.0
167 SER   15    -180.0     180.0
168 GLY   16    -180.0     180.0
169 ARG   17    -180.0     180.0
170 PRO   18    -180.0     180.0
171 PRO   19    -180.0     180.0
172 PRO   20    -180.0     180.0
173 SER   21    -180.0     180.0
174 D     22    -180.0     180.0
175 nsup= 20
176  nsup=          20  nstart_sup=           2  nstart_seq=           2
177  NZ_START=           2  NZ_END=          21
178  IZ_SC=           0
179  Contact order:  0.000000000000000E+000
180  Shifting contacts:           2           2
181 Random-generated initial geometry.
182
183 Geometry of the virtual chain.
184   Res           d     Theta     Gamma       Dsc     Alpha      Beta 
185 D     1     0.000     0.000     0.000     0.000     0.000     0.000
186 ASN   2     3.800     0.000     0.000     1.684    94.624  -178.093
187 LEU   3     3.800    92.239     0.000     1.939   102.878   -79.236
188 TYR   4     3.800    92.239  -180.000     2.484   163.195    -7.440
189 ILE   5     3.800    90.357    45.849     1.776   144.011  -104.516
190 GLN   6     3.800    89.090    55.194     2.240   170.986  -139.318
191 TRP   7     3.800    88.657    49.396     2.605   123.650   -21.913
192 LEU   8     3.800    93.032    48.298     1.939   151.723  -105.899
193 LYS   9     3.800    94.826    46.843     2.541   106.974   -55.960
194 ASP  10     3.800    87.966    59.623     1.709   153.447  -128.646
195 GLY  11     3.800    89.908    56.679     0.000     0.000     0.000
196 GLY  12     3.800   100.181   -87.304     0.000     0.000     0.000
197 PRO  13     3.800   109.073  -127.499     1.345   101.771  -115.580
198 SER  14     3.800    89.537     4.261     1.150   144.515  -129.218
199 SER  15     3.800    91.815    66.108     1.150   161.047  -100.177
200 GLY  16     3.800   101.784    70.140     0.000     0.000     0.000
201 ARG  17     3.800    92.201  -108.949     3.020   139.846  -132.716
202 PRO  18     3.800   133.225   103.824     1.345   115.610  -118.024
203 PRO  19     3.800   121.502  -122.527     1.345   118.575  -122.417
204 PRO  20     3.800   117.950   -90.285     1.345   118.959  -126.207
205 SER  21     3.800   114.201  -108.328     1.150   128.925   -37.341
206 D    22     3.800   114.201   180.000     0.000     0.000     0.000
207 Energy evaluation or minimization calculation.
208
209 Conformations will be energy-minimized.
210 ********************************************************************************
211
212  Time for energy evaluation  0.000000000000000E+000
213
214 Virtual-chain energies:
215
216 EVDW=     -2.897067E+01 WEIGHT=    1.352790D+00 (SC-SC)
217 EVDW2=     5.054652E+01 WEIGHT=    1.593040D+00 (SC-p)
218 EES=      -1.018887E+02 WEIGHT=    7.153400D-01 (p-p)
219 EVDWPP=   -2.665276E+01 WEIGHT=    1.137100D-01 (p-p VDW)
220 ESTR=      3.636156E-27 WEIGHT=    1.000000D+00 (stretching)
221 EBE=      -3.261471E+01 WEIGHT=    1.138730D+00 (bending)
222 ESC=       5.789477E+01 WEIGHT=    1.625800D-01 (SC local)
223 ETORS=     1.479510E+01 WEIGHT=    1.985990D+00 (torsional)
224 ETORSD=    1.006444E+00 WEIGHT=    1.570690D+00 (double torsional)
225 EHBP=      0.000000E+00 WEIGHT=    1.000000D+00 (SS bridges & dist. cnstr.)
226 ECORR4=   -7.263447E+01 WEIGHT=    4.288700D-01 (multi-body)
227 ECORR5=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
228 ECORR6=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
229 EELLO=    -4.196810E+00 WEIGHT=    1.603600D-01 (electrostatic-local)
230 ETURN3=    1.924843E+01 WEIGHT=    1.687220D+00 (turns, 3rd order)
231 ETURN4=    7.171768E-01 WEIGHT=    6.623000D-01 (turns, 4th order)
232 ETURN6=    0.000000E+00 WEIGHT=    0.000000D+00 (turns, 6th order)
233 ESCCOR=    0.000000E+00 WEIGHT=    0.000000D+00 (backbone-rotamer corr)
234 EDIHC=     0.000000E+00 (dihedral angle constraints)
235 ESS=       0.000000E+00 (disulfide-bridge intrinsic energy)
236 UCONST=     0.000000E+00 (Constraint energy)
237 ETOT=     -3.021986E+01 (total)
238
239 SUMSL return code:   4 energy      -79.51931
240
241
242 Virtual-chain energies:
243
244 EVDW=     -5.120726E+01 WEIGHT=    1.352790D+00 (SC-SC)
245 EVDW2=     4.630730E+01 WEIGHT=    1.593040D+00 (SC-p)
246 EES=      -9.334462E+01 WEIGHT=    7.153400D-01 (p-p)
247 EVDWPP=   -2.041364E+01 WEIGHT=    1.137100D-01 (p-p VDW)
248 ESTR=      3.636156E-27 WEIGHT=    1.000000D+00 (stretching)
249 EBE=      -3.685480E+01 WEIGHT=    1.138730D+00 (bending)
250 ESC=       4.338842E+01 WEIGHT=    1.625800D-01 (SC local)
251 ETORS=     1.030614E+01 WEIGHT=    1.985990D+00 (torsional)
252 ETORSD=    3.510843E-01 WEIGHT=    1.570690D+00 (double torsional)
253 EHBP=      0.000000E+00 WEIGHT=    1.000000D+00 (SS bridges & dist. cnstr.)
254 ECORR4=   -6.845424E+01 WEIGHT=    4.288700D-01 (multi-body)
255 ECORR5=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
256 ECORR6=    0.000000E+00 WEIGHT=    0.000000D+00 (multi-body)
257 EELLO=    -5.185692E+00 WEIGHT=    1.603600D-01 (electrostatic-local)
258 ETURN3=    1.856072E+01 WEIGHT=    1.687220D+00 (turns, 3rd order)
259 ETURN4=   -3.252103E+00 WEIGHT=    6.623000D-01 (turns, 4th order)
260 ETURN6=    0.000000E+00 WEIGHT=    0.000000D+00 (turns, 6th order)
261 ESCCOR=    0.000000E+00 WEIGHT=    0.000000D+00 (backbone-rotamer corr)
262 EDIHC=     0.000000E+00 (dihedral angle constraints)
263 ESS=       0.000000E+00 (disulfide-bridge intrinsic energy)
264 UCONST=     0.000000E+00 (Constraint energy)
265 ETOT=     -7.951931E+01 (total)
266
267 Geometry of the virtual chain.
268   Res           d     Theta     Gamma       Dsc     Alpha      Beta 
269 D     1     0.000     0.000     0.000     0.000     0.000     0.000
270 ASN   2     3.800     0.000     0.000     1.684   145.636   -16.341
271 LEU   3     3.800   115.891     0.000     1.939   137.106   -20.479
272 TYR   4     3.800    91.470   164.571     2.484   172.221    20.304
273 ILE   5     3.800    91.922    85.116     1.776   143.240   -90.802
274 GLN   6     3.800    90.239    44.027     2.240   114.569  -126.866
275 TRP   7     3.800    89.768    44.854     2.605   125.459   -26.243
276 LEU   8     3.800    89.846    45.610     1.939   153.751  -173.458
277 LYS   9     3.800    90.114    64.235     2.541    83.102   -75.669
278 ASP  10     3.800    91.773    38.260     1.709   135.700  -110.240
279 GLY  11     3.800    92.683    65.542     0.000   180.000   180.000
280 GLY  12     3.800   109.008   -65.773     0.000   180.000   180.000
281 PRO  13     3.800   149.210  -162.848     1.345   128.883  -155.795
282 SER  14     3.800    91.775   -38.027     1.150   133.697  -101.096
283 SER  15     3.800    91.305    50.569     1.150   117.134   -81.693
284 GLY  16     3.800   123.414    71.938     0.000   180.000   180.000
285 ARG  17     3.800    97.656   -79.129     3.020   126.102   -91.006
286 PRO  18     3.800   134.525    87.337     1.345   121.053  -126.854
287 PRO  19     3.800   122.674  -109.034     1.345   109.296  -110.132
288 PRO  20     3.800   117.844  -106.682     1.345   118.782  -122.215
289 SER  21     3.800   118.382  -134.825     1.150   146.519  -119.631
290 D    22     3.800   120.888  -118.739     0.000   180.000   180.000
291 SUMSL return code:  4
292 # of energy evaluations:                 520
293 # of energy evaluations/sec:        3460.000
294
295
296 ***** Computation time:    0 hours  0 minutes  0 seconds *****