Dear all,<br>              i am doing a sampling method and in a stage of completion But unfortunatelymy topol.top file is missing and am really helpless .i cant continue the simulation.i have the tpr file adn pdb file at last step.i tried to create the top file from <br>
<br><b>g_x2top -f beta_md17.tpr -o topol.top -ff select</b><br>when it prompts for a force field i selected charmm27 which is what i used for the simulation it give me the following error<br>Reading file beta_md17.tpr, VERSION 4.5.3 (single precision)<br>
Reading file beta_md17.tpr, VERSION 4.5.3 (single precision)<br><br>-------------------------------------------------------<br>Program g_x2top, VERSION 4.5.3<br>Source code file: g_x2top.c, line: 505<br><br>Fatal error:<br>
No or incorrect atomname2type.n2t file found (looking for charmm27.ff)<br>For more information and tips for troubleshooting, please check the GROMACS<br>website at <a href="http://www.gromacs.org/Documentation/Errors">http://www.gromacs.org/Documentation/Errors</a><br>
<br>i<b> tried to use pdb2gmx as well<br> pdb2gmx -f  2.pdb -o 2.gro -p topool.top </b><br>End terminus GLU-56: COO-<br>Opening force field file /sw/gromacs-4.5.3/share/gromacs/top/charmm27.ff/aminoacids.arn<br>Opening force field file /sw/gromacs-4.5.3/share/gromacs/top/charmm27.ff/dna.arn<br>
Opening force field file /sw/gromacs-4.5.3/share/gromacs/top/charmm27.ff/rna.arn<br>Checking for duplicate atoms....<br>Now there are 16 residues with 247 atoms<br>Chain time...<br>Making bonds...<br>Number of bonds was 250, now 250<br>
Generating angles, dihedrals and pairs...<br>Before cleaning: 649 pairs<br>Before cleaning: 654 dihedrals<br>Keeping all generated dihedrals<br>Making cmap torsions...There are   14 cmap torsion pairs<br>There are  654 dihedrals,   35 impropers,  444 angles<br>
           640 pairs,      250 bonds and     0 virtual sites<br>Total mass 1860.928 a.m.u.<br>Total charge -3.000 e<br>Writing topology<br>Processing chain 2 (3 atoms, 3 residues)<br>There are 0 donors and 0 acceptors<br>
There are 0 hydrogen bonds<br>Warning: Starting residue NA2139 in chain not identified as Protein/RNA/DNA.<br>Warning: Starting residue NA2140 in chain not identified as Protein/RNA/DNA.<br>Warning: Starting residue NA2141 in chain not identified as Protein/RNA/DNA.<br>
Problem with chain definition, or missing terminal residues.<br>This chain does not appear to contain a recognized chain molecule.<br>If this is incorrect, you can edit residuetypes.dat to modify the behavior.<br>8 out of 8 lines of specbond.dat converted successfully<br>
<br>-------------------------------------------------------<br>Program pdb2gmx, VERSION 4.5.3<br>Source code file: resall.c, line: 581<br><br>Fatal error:<br>Residue &#39;NA&#39; not found in residue topology database<br>
For more information and tips for troubleshooting, please check the GROMACS<br>website at <a href="http://www.gromacs.org/Documentation/Errors">http://www.gromacs.org/Documentation/Errors</a><br><br>