[gmx-users] A very strange problem about version of pdb2gmx

英雄不再寂寞 xiaowu759 at qq.com
Sun Oct 17 09:07:35 CEST 2010


Dear gmxers,
  Yesterday, I posted a mail titled "pdb2gmx stop at "AtomType 1" " where I mentioned that I used gmx-4.5.1. The resulting problem is again posted below. Today, I turn to the gmx-4.0.7, and use the same pdb and rtp files with the name modified to compatible to the version. Strangely, it is found that the pdb2gmx of this version can generate the top file without any errors. Could you give me clues to why the latest version can not whereas the older version can? To get the clear help, my pdb and rtp files are also posted below. Thanks a lot for any reply.
 Chaofu Wu, Dr.
  
 My command line and the screen output:
xiaowu759 at linux-s38y:~/workshop> pdb2gmx -f deta.pdb -o deta.gro -p deta.top -i deta.itp -ter
                         :-)  G  R  O  M  A  C  S  (-:
               Gromacs Runs One Microsecond At Cannonball Speeds
                            :-)  VERSION 4.5.1  (-:
......
Opening force field file oplsaa.ff/aminoacids.r2b
Reading deta.pdb...
Read 20 atoms
Analyzing pdb file
Splitting PDB chains based on TER records or changing chain id.
There are 1 chains and 0 blocks of water and 1 residues with 20 atoms
  chain  #res #atoms
  1 ' '     1     20  
All occupancies are one
Opening force field file oplsaa.ff/atomtypes.atp
Atomtype 1
 
My pdb file:
ATOM      1  N1  MOL     1      -5.199   1.343   0.334  1.00  0.00           N    
ATOM      2  C2  MOL     1      -4.582   0.639   1.470  1.00  0.00           C    
ATOM      3  C3  MOL     1      -3.897   1.462   2.593  1.00  0.00           C    
ATOM      4  N4  MOL     1      -3.194   0.615   3.574  1.00  0.00           N    
ATOM      5  C5  MOL     1      -2.373   1.331   4.566  1.00  0.00           C    
ATOM      6  C6  MOL     1      -1.520   0.397   5.464  1.00  0.00           C    
ATOM      7  N7  MOL     1      -0.684   1.133   6.432  1.00  0.00           N    
ATOM      8  H12 MOL     1      -4.614   2.067  -0.035  1.00  0.00           H    
ATOM      9  H11 MOL     1      -6.128   1.694   0.483  1.00  0.00           H    
ATOM     10  H22 MOL     1      -3.826  -0.053   1.056  1.00  0.00           H    
ATOM     11  H21 MOL     1      -5.345  -0.016   1.926  1.00  0.00           H    
ATOM     12  H31 MOL     1      -4.628   2.144   3.065  1.00  0.00           H    
ATOM     13  H32 MOL     1      -3.141   2.116   2.124  1.00  0.00           H    
ATOM     14  H41 MOL     1      -3.854   0.018   4.043  1.00  0.00           H    
ATOM     15  H51 MOL     1      -2.980   2.003   5.204  1.00  0.00           H    
ATOM     16  H52 MOL     1      -1.679   1.992   4.016  1.00  0.00           H    
ATOM     17  H62 MOL     1      -0.857  -0.180   4.795  1.00  0.00           H    
ATOM     18  H61 MOL     1      -2.188  -0.346   5.939  1.00  0.00           H    
ATOM     19  H72 MOL     1      -0.023   0.541   6.916  1.00  0.00           H    
ATOM     20  H71 MOL     1      -1.253   1.554   7.143  1.00  0.00           H    
TER
 My rtp file:
[ DET ]
 [ atoms ]
    N1    opls_900    -0.900      1
   H11    opls_909     0.360      1
   H12    opls_909     0.360      1
    C2    opls_906     0.060      1
   H21    opls_911     0.060      1
   H22    opls_911     0.060      1
    C3    opls_907     0.080      2
   H31    opls_911     0.060      2
   H32    opls_911     0.060      2
    N4    opls_901    -0.780      2
   H41    opls_910     0.380      2
    C5    opls_907     0.080      2
   H51    opls_911     0.060      2
   H52    opls_911     0.060      2
    C6    opls_906     0.060      3
   H61    opls_911     0.060      3
   H62    opls_911     0.060      3
    N7    opls_900    -0.900      3
   H71    opls_909     0.360      3
   H72    opls_909     0.360      3
 [ bonds ]
    N1   H11
    N1   H12
    N1    C2
    C2   H21
    C2   H22
    C2    C3
    C3   H31
    C3   H32
    C3    N4
    N4   H41
    N4    C5
    C5   H51
    C5   H52
    C5    C6
    C6   H61
    C6   H62
    C6    N7
    N7   H71
    N7   H72
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20101017/b0535a6a/attachment.html>


More information about the gromacs.org_gmx-users mailing list