[gmx-users] Capping (neuralising) of peptide

David van der Spoel spoel at xray.bmc.uu.se
Thu Sep 30 15:11:16 CEST 2004


On Thu, 2004-09-30 at 14:54, Alok wrote:
> Hello
>       Thanks for your suggestion.
> i used that options but now i am getting following Fatal error.
> i used  pdb2gmx -ignh -ter -f peptide.pdb -o peptide.gro -p peptide.top
> 
> 
your pdb file is severely crippled. you have multiple peptides in it all
marked with labels in the wrong place. get rid of the B's and C's
and put only one peptide in the file.

> 
> 
> Opening library file /usr/local/share/gromacs/top/FF.dat
> 
> Select the Force Field:
>  0: GROMOS96 43a1 Forcefield (official distribution)
>  1: GROMOS96 43b1 Vacuum Forcefield (official distribution)
>  2: GROMOS96 43a2 Forcefield (development) (improved alkane dihedrals)
>  3: OPLS-AA/L all-atom force field (2001 aminoacid dihedrals)
>  4: Gromacs Forcefield (see manual)
>  5: gmx Forcefield with hydrogens for NMR stuff (Do NOT use for new runs)
> 0
> Looking whether force field file ffG43a1.rtp exists
> Opening library file /usr/local/share/gromacs/top/ffG43a1.rtp
> Opening library file /usr/local/share/gromacs/top/aminoacids.dat
> Reading peptide.pdb...
> Read 64 atoms
> Opening library file /usr/local/share/gromacs/top/xlateat.dat
> 23 out of 23 lines of xlateat.dat converted succesfully
> Analyzing pdb file
> There are 1 chains and 0 blocks of water and 16 residues with 64 atoms
> 
>   chain  #res #atoms
>   1 ' '    16     64
> 
> WARNING: there were 1 atoms with zero occupancy and 0 atoms with
>          occupancy unequal to one (out of 64 atoms). Check your pdb file.
> Opening library file /usr/local/share/gromacs/top/ffG43a1.atp
> Atomtype 50
> Reading residue database... (ffG43a1)
> Opening library file /usr/local/share/gromacs/top/ffG43a1.rtp
> Residue 96
> Sorting it all out...
> Opening library file /usr/local/share/gromacs/top/ffG43a1.hdb
> Opening library file /usr/local/share/gromacs/top/ffG43a1-n.tdb
> Opening library file /usr/local/share/gromacs/top/ffG43a1-c.tdb
> 
> Back Off! I just backed up peptide.top to ./#peptide.top.16#
> Processing chain 1 (64 atoms, 16 residues)
> Opening library file /usr/local/share/gromacs/top/specbond.dat
> 5 out of 5 lines of specbond.dat converted succesfully
> There are 12 donors and 16 acceptors
> There are 19 hydrogen bonds
> Checking for duplicate atoms....
> Select N-terminus type (start)
>  0: GLY-NH3+
>  1: NH2
>  2: None
> 1
> N-terminus: NH2
> Select C-terminus type (end)
>  0: COO-
>  1: COOH
>  2: None
> 1
> C-terminus: COOH
> Fatal error: Atom C not found in residue GLY16 while adding hydrogens
> 
> 
> 
> 
> My PDB file :-
> 
> 
> ATOM      1  N   GLY     1      15.055  26.899  27.745  1.00  0.00        
>   N
> ATOM      2  CA  GLY     1      14.262  27.999  27.204  1.00  0.00        
>   C
> ATOM      3  C   GLY     1      14.753  29.328  27.725  1.00  0.00        
>   C
> ATOM      4  O   GLY     1      15.669  29.411  28.537  1.00  0.00        
>   O
> ATOM      5 1H   GLY     1      15.783  27.276  28.371  1.00  0.00        
>   H
> ATOM      6 2H   GLY     1      14.444  26.262  28.276  1.00  0.00        
>   H
> ATOM      7 1HA  GLY     1      14.313  28.000  26.100  1.00  0.00        
>   H
> ATOM      8 2HA  GLY     1      13.198  27.873  27.476  1.00  0.00        
>   H
> ATOM      9  N   PRO     1B     14.138  30.394  27.257  1.00  0.00        
>   N
> ATOM     10  CA  PRO     1B     14.502  31.745  27.655  1.00  0.00        
>   C
> ATOM     11  C   PRO     1B     14.016  32.068  29.078  1.00  0.00        
>   C
> ATOM     12  O   PRO     1B     12.967  31.584  29.504  1.00  0.00        
>   O
> ATOM     13  CB  PRO     1B     13.737  32.594  26.628  1.00  0.00        
>   C
> ATOM     14  CG  PRO     1B     12.506  31.737  26.294  1.00  0.00        
>   C
> ATOM     15  CD  PRO     1B     13.059  30.311  26.295  1.00  0.00        
>   C
> ATOM     16  HA  PRO     1B     15.585  31.876  27.616  1.00  0.00        
>   H
> ATOM     17 1HB  PRO     1B     13.474  33.584  27.007  1.00  0.00        
>   H
> ATOM     18 2HB  PRO     1B     14.352  32.706  25.731  1.00  0.00        
>   H
> ATOM     19 1HG  PRO     1B     11.763  31.839  27.088  1.00  0.00        
>   H
> ATOM     20 2HG  PRO     1B     12.057  32.005  25.335  1.00  0.00        
>   H
> ATOM     21 1HD  PRO     1B     12.312  29.578  26.607  1.00  0.00        
>   H
> ATOM     22 2HD  PRO     1B     13.464  30.046  25.315  1.00  0.00        
>   H
> ATOM     23  N   GLY     1C     14.726  32.860  29.852  1.00  0.00        
>   N
> ATOM     24  CA  GLY     1C     14.238  33.314  31.151  1.00  0.00        
>   C
> ATOM     25  C   GLY     1C     13.769  34.748  31.082  1.00  0.00        
>   C
> ATOM     26  O   GLY     1C     13.445  35.378  32.087  1.00  0.00        
>   O
> ATOM     27  H   GLY     1C     15.654  33.172  29.532  1.00  0.00        
>   H
> ATOM     28 1HA  GLY     1C     13.403  32.675  31.489  1.00  0.00        
>   H
> ATOM     29 2HA  GLY     1C     15.036  33.230  31.910  1.00  0.00        
>   H
> ATOM     30  O   GLY     2      13.704  35.373  29.853  1.00  0.00        
>   O
> ATOM     31  HO  GLY     2      13.391  36.273  29.966  1.00  0.00        
>   H
> TER
> ATOM     32  N   GLY     1      56.800  33.657  40.211  1.00  0.00        
>   N
> ATOM     33  CA  GLY     1      55.946  34.770  39.811  1.00  0.00        
>   C
> ATOM     34  C   GLY     1      56.520  36.088  40.276  1.00  0.00        
>   C
> ATOM     35  O   GLY     1      57.546  36.150  40.944  1.00  0.00        
>   O
> ATOM     36 1H   GLY     1      57.615  34.017  40.728  1.00  0.00        
>   H
> ATOM     37 2H   GLY     1      56.269  33.014  40.815  1.00  0.00        
>   H
> ATOM     38 1HA  GLY     1      55.833  34.791  38.711  1.00  0.00        
>   H
> ATOM     39 2HA  GLY     1      54.932  34.647  40.234  1.00  0.00        
>   H
> ATOM     40  N   PRO     1B     55.851  37.165  39.921  1.00  0.00        
>   N
> ATOM     41  CA  PRO     1B     56.282  38.506  40.283  1.00  0.00        
>   C
> ATOM     42  C   PRO     1B     56.014  38.808  41.767  1.00  0.00        
>   C
> ATOM     43  O   PRO     1B     55.036  38.322  42.335  1.00  0.00        
>   O
> ATOM     44  CB  PRO     1B     55.382  39.378  39.395  1.00  0.00        
>   C
> ATOM     45  CG  PRO     1B     54.108  38.535  39.232  1.00  0.00        
>   C
> ATOM     46  CD  PRO     1B     54.642  37.105  39.127  1.00  0.00        
>   C
> ATOM     47  HA  PRO     1B     57.349  38.632  40.087  1.00  0.00        
>   H
> ATOM     48 1HB  PRO     1B     55.187  40.363  39.824  1.00  0.00        
>   H
> ATOM     49 2HB  PRO     1B     55.859  39.502  38.419  1.00  0.00        
>   H
> ATOM     50 1HG  PRO     1B     53.490  38.628  40.128  1.00  0.00        
>   H
> ATOM     51 2HG  PRO     1B     53.524  38.823  38.355  1.00  0.00        
>   H
> ATOM     52 1HD  PRO     1B     53.942  36.372  39.534  1.00  0.00        
>   H
> ATOM     53 2HD  PRO     1B     54.896  36.856  38.094  1.00  0.00        
>   H
> ATOM     54  N   GLY     1C     56.838  39.581  42.441  1.00  0.00        
>   N
> ATOM     55  CA  GLY     1C     56.550  40.015  43.806  1.00  0.00        
>   C
> ATOM     56  C   GLY     1C     56.089  41.453  43.830  1.00  0.00        
>   C
> ATOM     57  O   GLY     1C     55.923  42.068  44.883  1.00  0.00        
>   O
> ATOM     58  H   GLY     1C     57.711  39.892  41.993  1.00  0.00        
>   H
> ATOM     59 1HA  GLY     1C     55.769  39.375  44.253  1.00  0.00        
>   H
> ATOM     60 2HA  GLY     1C     57.451  39.913  44.437  1.00  0.00        
>   H
> TER
> ATOM     61  O   GLY     2      55.849  42.100  42.635  1.00  0.00        
>   O
> ATOM     62  HO  GLY     2      55.564  43.000  42.808  1.00  0.00        
>   H
> ATOM     63  N   GLY     1      51.433  80.760  35.677  1.00  0.00        
>   N
> ATOM     64  CA  GLY     1      50.640  81.859  35.136  1.00  0.00        
>   C
> ATOM     65  C   GLY     1      51.130  83.189  35.657  1.00  0.00        
>   C
> ATOM     66  O   GLY     1      52.047  83.271  36.469  1.00  0.00        
>   O
> ATOM     67 1H   GLY     1      52.160  81.136  36.303  1.00  0.00        
>   H
> ATOM     68 2H   GLY     1      50.822  80.123  36.208  1.00  0.00        
>   H
> ATOM     69 H    GLY     1      50.690  81.860  34.032  1.00  0.00        
>   H
> ATOM     70 2HA  GLY     1      49.575  81.733  35.408  1.00  0.00        
>   H
> ATOM     71  N   PRO     1B     50.515  84.254  35.189  1.00  0.00        
>   N
> ATOM     72  CA  PRO     1B     50.880  85.605  35.587  1.00  0.00        
>   C
> ATOM     73  C   PRO     1B     50.393  85.929  37.009  1.00  0.00        
>   C
> ATOM     74  O   PRO     1B     49.345  85.444  37.436  1.00  0.00        
>   O
> ATOM     75  CB  PRO     1B     50.115  86.454  34.560  1.00  0.00        
>   C
> ATOM     76  CG  PRO     1B     48.884  85.597  34.226  1.00  0.00        
>   C
> ATOM     77  CD  PRO     1B     49.437  84.171  34.227  1.00  0.00        
>   C
> ATOM     78  HA  PRO     1B     51.963  85.737  35.548  1.00  0.00        
>   H
> ATOM     79 1HB  PRO     1B     49.852  87.444  34.939  1.00  0.00        
>   H
> ATOM     80 2HB  PRO     1B     50.730  86.566  33.663  1.00  0.00        
>   H
> ATOM     81 1HG  PRO     1B     48.140  85.700  35.020  1.00  0.00        
>   H
> ATOM     82 2HG  PRO     1B     48.434  85.866  33.267  1.00  0.00        
>   H
> ATOM     83 1HD  PRO     1B     48.689  83.438  34.538  1.00  0.00        
>   H
> ATOM     84 2HD  PRO     1B     49.841  83.907  33.247  1.00  0.00        
>   H
> ATOM     85  N   GLY     1C     51.104  86.721  37.784  1.00  0.00        
>   N
> ATOM     86  CA  GLY     1C     50.615  87.174  39.083  1.00  0.00        
>   C
> ATOM     87  C   GLY     1C     50.147  88.609  39.014  1.00  0.00        
>   C
> ATOM     88  O   GLY     1C     49.823  89.239  40.019  1.00  0.00        
>   O
> ATOM     89  H   GLY     1C     52.032  87.032  37.464  1.00  0.00        
>   H
> ATOM     90 1HA  GLY     1C     49.780  86.535  39.421  1.00  0.00        
>   H
> ATOM     91 2HA  GLY     1C     51.414  87.090  39.842  1.00  0.00        
>   H
> ATOM     92  O   GLY     2      50.081  89.233  37.785  1.00  0.00        
>   O
> ATOM     93  HO  GLY     2      49.769  90.133  37.898  1.00  0.00        
>   H
> TER
> ATOM     94  N   GLY     1      11.093  76.967  23.302  1.00  0.00        
>   N
> ATOM     95  CA  GLY     1      10.300  78.066  22.761  1.00  0.00        
>   C
> ATOM     96  C   GLY     1      10.791  79.396  23.282  1.00  0.00        
>   C
> ATOM     97  O   GLY     1      11.708  79.478  24.094  1.00  0.00        
>   O
> ATOM     98 1H   GLY     1      11.821  77.343  23.928  1.00  0.00        
>   H
> ATOM     99 2H   GLY     1      10.482  76.330  23.833  1.00  0.00        
>   H
> ATOM    100 1HA  GLY     1      10.351  78.068  21.657  1.00  0.00        
>   H
> ATOM    101 2HA  GLY     1       9.236  77.941  23.032  1.00  0.00        
>   H
> ATOM    102  N   PRO     1B     10.176  80.462  22.814  1.00  0.00        
>   N
> ATOM    103  CA  PRO     1B     10.540  81.812  23.212  1.00  0.00        
>   C
> ATOM    104  C   PRO     1B     10.054  82.136  24.634  1.00  0.00        
>   C
> ATOM    105  O   PRO     1B      9.005  81.651  25.061  1.00  0.00        
>   O
> ATOM    106  CB  PRO     1B      9.775  82.661  22.185  1.00  0.00        
>   C
> ATOM    107  CG  PRO     1B      8.544  81.804  21.851  1.00  0.00        
>   C
> ATOM    108  CD  PRO     1B      9.097  80.378  21.852  1.00  0.00        
>   C
> ATOM    109  HA  PRO     1B     11.623  81.944  23.173  1.00  0.00        
>   H
> ATOM    110 1HB  PRO     1B      9.512  83.652  22.564  1.00  0.00        
>   H
> ATOM    111 2HB  PRO     1B     10.390  82.773  21.288  1.00  0.00        
>   H
> ATOM    112 1HG  PRO     1B      7.801  81.907  22.645  1.00  0.00        
>   H
> ATOM    113 2HG  PRO     1B      8.095  82.073  20.892  1.00  0.00        
>   H
> ATOM    114 1HD  PRO     1B      8.350  79.645  22.163  1.00  0.00        
>   H
> ATOM    115 2HD  PRO     1B      9.502  80.114  20.872  1.00  0.00        
>   H
> ATOM    116  N   GLY     1C     10.764  82.928  25.408  1.00  0.00        
>   N
> ATOM    117  CA  GLY     1C     10.276  83.382  26.708  1.00  0.00        
>   C
> ATOM    118  C   GLY     1C      9.807  84.816  26.639  1.00  0.00        
>   C
> ATOM    119  O   GLY     1C      9.483  85.446  27.644  1.00  0.00        
>   O
> ATOM    120  H   GLY     1C     11.693  83.239  25.089  1.00  0.00        
>   H
> ATOM    121 1HA  GLY     1C      9.441  82.743  27.046  1.00  0.00        
>   H
> ATOM    122 2HA  GLY     1C     11.074  83.297  27.467  1.00  0.00        
>   H
> ATOM    123  O   GLY      2       9.742  85.440  25.410  1.00  0.00       
>    O
> ATOM    124  HO  GLY      2       9.429  86.341  25.523  1.00  0.00       
>    H
> END
> 
> 
> 
> I also tryied to find out GLY16 residue but i didnt find anywhere in pdb
> file.
> 
> Thanks!
> 
>                                        Alok
> 
> 
> 
> _______________________________________________
> gmx-users mailing list
> gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  	75124 Uppsala, Sweden
phone:	46 18 471 4205		fax: 46 18 511 755
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++





More information about the gromacs.org_gmx-users mailing list