# Merged.hdb file needs improvement by this error appearance?

**URL:** <https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042>\
**Category:** User discussions\
**Created:** [November 4, 2020, 9:25am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042 "2020-11-04T09:25:08Z")\
**Posts on this page:** 17\
**Page:** 1

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 4, 2020, 9:25am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/1 "2020-11-04T09:25:08Z")

</div>

Hey guys,  
I am using CHARMM36 for my catalase peotein.  
I have gt this error

 ![1](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/1X/b3ecaac73e2325d604f532039fa38eb80b352869.jpeg)  
Does it mean that "merged.hdb " file must be edited?  
how this can be solved?  
-Sara

---

<div class="post-metadata">

**Author:** ![cblau](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/cblau/32/5_2.png) [@cblau](https://gromacs.bioexcel.eu/u/cblau)\
**Post date:** [November 4, 2020, 4:33pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/2 "2020-11-04T16:33:34Z")

</div>

Hi Sara, if you did not use it already, `-ignh` might help in forcing GROMACS to add back hyrdrogen atoms that are not found in your input `pdb`.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 4, 2020, 8:42pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/3 "2020-11-04T20:42:55Z")

</div>

Hi,  
I did this but both missing atom warning and hydrogen errors remained. moreover, the simulation stopped in NPT step.  
Thanks  
-SARA

---

<div class="post-metadata">

**Author:** ![alevilla](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/alevilla/32/439_2.png) [@alevilla](https://gromacs.bioexcel.eu/u/alevilla)\
**Post date:** [November 5, 2020, 9:11am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/4 "2020-11-05T09:11:55Z")

</div>

Hi,  
Hydrogen atoms are missing in the HEME group in your pdb file. You can add the hydrogen to the pdb (with the help of visualization program) or build an entry in merge.hdb for HEME group (for syntax see [http://manual.gromacs.org/current/reference-manual/topologies/pdb2gmx-input-files.html#pdb2gmxfiles](http://manual.gromacs.org/current/reference-manual/topologies/pdb2gmx-input-files.html#pdb2gmxfiles))  
Now gmx pdb2gmx does not know where to add the missing hydrogens, thus the option.  
Best regards  
Alessandra

---

<div class="post-metadata">

**Author:** ![yasser.almeida](https://avatars.discourse-cdn.com/v4/letter/y/76d3ee/32.png) [@yasser.almeida](https://gromacs.bioexcel.eu/u/yasser.almeida)\
**Post date:** [November 5, 2020, 9:24am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/5 "2020-11-05T09:24:47Z")

</div>

Hi,

I faced the same problem with a system containing a new residue and did that (`-ignh` and adding the H with another software). However, during the equilibration, the new residue falls apart, so I guess I have to necessarily update the `merged.hdb`, with the new hydrogens, right?

Best

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [November 5, 2020, 2:22pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/6 "2020-11-05T14:22:14Z")

</div>

This is a general problem with any molecule - if the nomenclature of the atoms provided does not match the residue definition, either (1) you have to change the atom names in the input coordinate file or (2) use `-ignh` in conjunction with a suitable `.hdb` entry.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 6, 2020, 1:51pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/7 "2020-11-06T13:51:47Z")

</div>

Dear Justin,  
I changed the atom names based on merged.rtp file. however, the following message appeared.

 ![image](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/1X/f53aecbabe6c3ad47e61c650a63d8b64da05ceea.png)  
Moreover, using -ignh causes the appearance of the missing atom error.

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [November 6, 2020, 4:29pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/8 "2020-11-06T16:29:09Z")

</div>

This is the same problem as before. You’re specifying an H atom in the input coordinates that is named incorrectly.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 9, 2020, 8:48pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/9 "2020-11-09T20:48:57Z")

</div>

I changed the names of heme hydrogens as well as their order bas[Screenshot 2020-11-10 001644|356x499]

 ![Screenshot 2020-11-10 001644](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/1X/8988e20e76140173e208899d7ea1d1cdf383a7be.jpeg)  
. Then, I imbedded them into the original PDB (1dgf). however, an error appeared:  
there was a hydrogen name HMA1 and gromacs couldn’t recognized it. the error said that there is no residue under the name of MA!!!

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [November 9, 2020, 9:15pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/10 "2020-11-09T21:15:38Z")

</div>

Please post the actual error message.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 10, 2020, 9:19am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/11 "2020-11-10T09:19:09Z")

</div>

![Screenshot from 2020-11-10 12-36-41](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/1X/4c7eafcf22fc7e332691c0d22efd7de3d1fb887e.png)

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 10, 2020, 9:20am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/12 "2020-11-10T09:20:04Z")

</div>

Thank you for your response.  
Following the use of -ignh , I got the error of missing atoms in HEME molecules.  
-Sara

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 10, 2020, 9:21am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/13 "2020-11-10T09:21:20Z")

</div>

hi,  
Thank you for your response.  
You mean trying embedding the Heme in my pdb into the merged.hdb?  
-Sara

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [November 10, 2020, 10:56am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/14 "2020-11-10T10:56:14Z")

</div>

The error means you’ve broken the PDB format by changing the required spacing of the columns. Editing the `.hdb` file would be very complicated for a molecule like heme. Renaming the atoms in the PDB is probably the easiest approach but you have to preserve the fixed-format nature of the PDB file.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 10, 2020, 4:05pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/15 "2020-11-10T16:05:05Z")

</div>

you are right dear Justin.  
But the problem is that I have tried different space patterns and I got similar errors. :( for exmple, I followed the spaces in the pdb but I got similar errors. Is it possible there it has something to do with four-character hydrogen names?  
I tried three-character hydrogen names based on what is defined in the merged.hdb file. but again it showed an error. since it seems gromacs can only read the first hydrogen atoms : HA, HB, HC, and HD.

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [November 10, 2020, 6:55pm UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/16 "2020-11-10T18:55:13Z")

</div>

There is no problem with 4-character atom names. But the columns all have to be aligned properly. See, e.g. [https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/tutorials/pdbintro.html](https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/tutorials/pdbintro.html)

Your spacing is simply incorrect for the HEM entries, causing a misinterpretation of the lines.

---

<div class="post-metadata">

**Author:** ![Sara](https://avatars.discourse-cdn.com/v4/letter/s/a5b964/32.png) [@Sara](https://gromacs.bioexcel.eu/u/Sara)\
**Post date:** [November 15, 2020, 11:42am UTC](https://gromacs.bioexcel.eu/t/merged-hdb-file-needs-improvement-by-this-error-appearance/1042/17 "2020-11-15T11:42:03Z")

</div>

Dear Justin,  
I appreciate your kind assistance. Following a precise modification of the PDB, it worked.  
Thank You  
-Sara
