<div>Dear All<br>I used the following command sequentially to prepare file for energy<br>minimization and subsequent MD run.<br>1. pdb2gmx  -f *.pdb -o seq.gro -p seq.top<br>2. editconf -f seq.gro -o seq_box.gro -d 1.0 -bt cubic<br>
3. genbox   -cp seq_box.gro -cs spc216.gro -o seq_b4ion.gro -p seq.top<br>4. grompp   -c seq_b4ion.gro -p seq.top -o seq_b4ion.tpr -f em.mdp<br>grompp gives following error.processing topology...<br>Opening library file /usr/local/gromacs/share/gromacs/top/ffG43a1.itp<br>
Opening library file /usr/local/gromacs/share/gromacs/top/ffG43a1nb.itp<br>Opening library file /usr/local/gromacs/share/gromacs/top/ffG43a1bon.itp<br>Opening library file /usr/local/gromacs/share/gromacs/top/ff_dum.itp<br>
Generated 279 of the 1225 non-bonded parameter combinations<br><br>ERROR 1 [file seq.top, line 1965]:<br> No default G96Bond types<br><br><br>ERROR 2 [file seq.top, line 5271]:<br> No default G96Angle types<br><br><br>ERROR 3 [file seq.top, line 5272]:<br>
 No default G96Angle types<br><br><br>ERROR 4 [file seq.top, line 5648]:<br> No default G96Angle types<br><br><br>ERROR 5 [file seq.top, line 5653]:<br> No default G96Angle types<br><br><br>ERROR 6 [file seq.top, line 5654]:<br>
 No default G96Angle types<br><br><br>ERROR 7 [file seq.top, line 5655]:<br> No default G96Angle types<br><br><br>ERROR 8 [file seq.top, line 5656]:<br> No default G96Angle types<br><br><br>ERROR 9 [file seq.top, line 6201]:<br>
 No default Proper Dih. types<br><br>Opening library file /usr/local/gromacs/share/gromacs/top/spc.itp<br>Opening library file /usr/local/gromacs/share/gromacs/top/ions.itp<br>Excluding 3 bonded neighbours molecule type &#39;Protein_A&#39;<br>
Excluding 2 bonded neighbours molecule type &#39;SOL&#39;<br>Excluding 2 bonded neighbours molecule type &#39;SOL&#39;<br><br>NOTE 1 [file seq.top, line 6932]:<br> System has non-zero total charge: 7.000001e+00<br><br><br>
<br>processing coordinates...<br>double-checking input for internal consistency...<br>renumbering atomtypes...<br>converting bonded parameters...<br><br>There was 1 note<br><br>-------------------------------------------------------<br>
Program grompp, VERSION 4.0.7<br>Source code file: grompp.c, line: 986<br><br>Fatal error:<br>There were 9 errors in input file(s)<br>-------------------------------------------------------<br> pdb2gmx works properly using ff43a1 forcefield. My protein contains<br>
Heme.  This heme is bonded to MET SD and NE2 HIS as well as NA NB NC ND of heme porphyrins. It is actually these six bond which is created in topology file with pdb2gmx using specbond.dat, giving error in grompp. I tried to find these bonds in in built files of gromacs i.e. /share/top/ but I couldnt. Can any one help me and tell me where can I find these bonds/angle parameters and what should I add in ffG43a1 .itp .rtp files</div>

<div>please help.</div>
<div>Shahid Nayeem<br><br><br> </div>