Advertisement
Guest User

Untitled

a guest
Apr 8th, 2020
348
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Bash 9.21 KB | None | 0 0
  1. maiia@SpectrumAnalayzer:~/Рабочий стол/Protein$ gmx pdb2gmx -f 2OOB.pdb -o peptide.pdb -ff amber99sbmut -water tip3p -ignh
  2.                   :-) GROMACS - gmx pdb2gmx, VERSION 5.1.2 (-:
  3.  
  4.                             GROMACS is written by:
  5.      Emile Apol      Rossen Apostolov  Herman J.C. Berendsen    Par Bjelkmar  
  6.  Aldert van Buuren   Rudi van Drunen     Anton Feenstra   Sebastian Fritsch
  7.   Gerrit Groenhof   Christoph Junghans   Anca Hamuraru    Vincent Hindriksen
  8.  Dimitrios Karkoulis    Peter Kasson        Jiri Kraus      Carsten Kutzner  
  9.     Per Larsson      Justin A. Lemkul   Magnus Lundborg   Pieter Meulenhoff
  10.    Erik Marklund      Teemu Murtola       Szilard Pall       Sander Pronk  
  11.    Roland Schulz     Alexey Shvetsov     Michael Shirts     Alfons Sijbers  
  12.    Peter Tieleman    Teemu Virolainen  Christian Wennberg    Maarten Wolf  
  13.                            and the project leaders:
  14.         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
  15.  
  16. Copyright (c) 1991-2000, University of Groningen, The Netherlands.
  17. Copyright (c) 2001-2015, The GROMACS development team at
  18. Uppsala University, Stockholm University and
  19. the Royal Institute of Technology, Sweden.
  20. check out http://www.gromacs.org for more information.
  21.  
  22. GROMACS is free software; you can redistribute it and/or modify it
  23. under the terms of the GNU Lesser General Public License
  24. as published by the Free Software Foundation; either version 2.1
  25. of the License, or (at your option) any later version.
  26.  
  27. GROMACS:      gmx pdb2gmx, VERSION 5.1.2
  28. Executable:   /usr/bin/gmx
  29. Data prefix:  /usr
  30. Command line:
  31.   gmx pdb2gmx -f 2OOB.pdb -o peptide.pdb -ff amber99sbmut -water tip3p -ignh
  32.  
  33.  
  34. Using the Amber99sbmut force field in directory /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff
  35.  
  36. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.r2b
  37. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.r2b
  38. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.r2b
  39. Reading 2OOB.pdb...
  40. Read 1056 atoms
  41. Analyzing pdb file
  42. Splitting chemical chains based on TER records or chain id changing.
  43. Moved all the water blocks to the end
  44. There are 2 chains and 2 blocks of water and 248 residues with 1056 atoms
  45.  
  46.   chain  #res #atoms
  47.   1 'A'    44    350  
  48.   2 'B'    72    574  
  49.   3 ' '    50     50  (only water)
  50.   4 ' '    82     82  (only water)
  51.  
  52.  
  53. WARNING: there were 0 atoms with zero occupancy and 2 atoms with
  54.          occupancy unequal to one (out of 1056 atoms). Check your pdb file.
  55.  
  56. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/atomtypes.atp
  57. Atomtype 103
  58. Reading residue database... (amber99sbmut)
  59. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.rtp
  60. Residue 98
  61. Sorting it all out...
  62. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.rtp
  63. Residue 114
  64. Sorting it all out...
  65. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/mutres.rtp
  66. Residue 666
  67. Sorting it all out...
  68. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.rtp
  69. Residue 682
  70. Sorting it all out...
  71. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.hdb
  72. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.hdb
  73. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.hdb
  74. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.n.tdb
  75. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.c.tdb
  76. Processing chain 1 'A' (350 atoms, 44 residues)
  77. Identified residue LEU1 as a starting terminus.
  78. Identified residue PHE44 as a ending terminus.
  79. 8 out of 8 lines of specbond.dat converted successfully
  80. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.arn
  81. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.arn
  82. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.arn
  83. Checking for duplicate atoms....
  84. Generating any missing hydrogen atoms and/or adding termini.
  85. Now there are 44 residues with 707 atoms
  86. Chain time...
  87. Making bonds...
  88. Number of bonds was 711, now 710
  89. Generating angles, dihedrals and pairs...
  90. Before cleaning: 1858 pairs
  91. Before cleaning: 1858 dihedrals
  92. Keeping all generated dihedrals
  93. Making cmap torsions...
  94. There are 1858 dihedrals,  139 impropers, 1284 angles
  95.           1846 pairs,      710 bonds and     0 virtual sites
  96. Total mass 4981.738 a.m.u.
  97. Total charge -2.000 e
  98. Writing topology
  99. Processing chain 2 'B' (574 atoms, 72 residues)
  100. Analysing hydrogen-bonding network for automated assignment of histidine
  101.  protonation. 109 donors and 113 acceptors were found.
  102. There are 165 hydrogen bonds
  103. Will use HISE for residue 68
  104. Identified residue MET1 as a starting terminus.
  105. Identified residue ARG72 as a ending terminus.
  106. 8 out of 8 lines of specbond.dat converted successfully
  107. Special Atom Distance matrix:
  108.                     MET1
  109.                      SD7
  110.    HIS68  NE2540   1.673
  111. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.arn
  112. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.arn
  113. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.arn
  114. Checking for duplicate atoms....
  115. Generating any missing hydrogen atoms and/or adding termini.
  116. Now there are 72 residues with 1174 atoms
  117. Chain time...
  118. Making bonds...
  119. Number of bonds was 1181, now 1180
  120. Generating angles, dihedrals and pairs...
  121. Before cleaning: 3128 pairs
  122. Before cleaning: 3148 dihedrals
  123. Keeping all generated dihedrals
  124. Making cmap torsions...
  125. There are 3148 dihedrals,  204 impropers, 2155 angles
  126.           3119 pairs,     1180 bonds and     0 virtual sites
  127. Total mass 8180.452 a.m.u.
  128. Total charge -1.000 e
  129. Writing topology
  130. Processing chain 3 (50 atoms, 50 residues)
  131. Problem with chain definition, or missing terminal residues.
  132. This chain does not appear to contain a recognized chain molecule.
  133. If this is incorrect, you can edit residuetypes.dat to modify the behavior.
  134. 8 out of 8 lines of specbond.dat converted successfully
  135. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.arn
  136. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.arn
  137. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.arn
  138. Checking for duplicate atoms....
  139. Generating any missing hydrogen atoms and/or adding termini.
  140. Now there are 50 residues with 150 atoms
  141. Making bonds...
  142. Number of bonds was 100, now 100
  143. Generating angles, dihedrals and pairs...
  144. Making cmap torsions...
  145. There are    0 dihedrals,    0 impropers,   50 angles
  146.              0 pairs,      100 bonds and     0 virtual sites
  147. Total mass 900.800 a.m.u.
  148. Total charge 0.000 e
  149. Processing chain 4 (82 atoms, 82 residues)
  150. Problem with chain definition, or missing terminal residues.
  151. This chain does not appear to contain a recognized chain molecule.
  152. If this is incorrect, you can edit residuetypes.dat to modify the behavior.
  153. 8 out of 8 lines of specbond.dat converted successfully
  154. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/aminoacids.arn
  155. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/dna.arn
  156. Opening force field file /home/maiia/.local/lib/python2.7/site-packages/pmx/data/mutff45/amber99sbmut.ff/rna.arn
  157. Checking for duplicate atoms....
  158. Generating any missing hydrogen atoms and/or adding termini.
  159. Now there are 82 residues with 246 atoms
  160. Making bonds...
  161. Number of bonds was 164, now 164
  162. Generating angles, dihedrals and pairs...
  163. Making cmap torsions...
  164. There are    0 dihedrals,    0 impropers,   82 angles
  165.              0 pairs,      164 bonds and     0 virtual sites
  166. Total mass 1477.312 a.m.u.
  167. Total charge 0.000 e
  168. Including chain 1 in system: 707 atoms 44 residues
  169. Including chain 2 in system: 1174 atoms 72 residues
  170. Including chain 3 in system: 150 atoms 50 residues
  171. Including chain 4 in system: 246 atoms 82 residues
  172. Now there are 2277 atoms and 248 residues
  173. Total mass in system 15540.302 a.m.u.
  174. Total charge in system -3.000 e
  175.  
  176. Writing coordinate file...
  177.         --------- PLEASE NOTE ------------
  178. You have successfully generated a topology from: 2OOB.pdb.
  179. The Amber99sbmut force field and the tip3p water model are used.
  180.         --------- ETON ESAELP ------------
  181.  
  182. gcq#498: "Problems worthy of attack prove their worth by hitting back." (Piet Hein)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement