# Error in HADDOCK3 Clustering (RMSD-based \[clustrmsd\]) during Protein-Protein Docking

**URL:** <https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857>\
**Category:** HADDOCK\
**Created:** [February 22, 2024, 3:14pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857 "2024-02-22T15:14:39Z")\
**Posts on this page:** 6\
**Page:** 1

<div class="post-metadata">

**Author:** ![YFeri](https://avatars.discourse-cdn.com/v4/letter/y/b4bc9f/32.png) [@YFeri](https://ask.bioexcel.eu/u/YFeri)\
**Post date:** [February 22, 2024, 3:14pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/1 "2024-02-22T15:14:39Z")

</div>

Hi,  
I’m currently using HADDOCK3 on an HPC for my protein-protein docking, and I’ve opted for the RMSD-based clustering method to cluster my models [clustrmsd]. However, I encountered an error during the docking process. The rigidbody, flexref, emref, and mdref stages run smoothly, but I encountered an issue specifically with clustrmsd.: [2024-02-21 18:02:49,331 **init** INFO] Running [seletop] module  
[2024-02-21 18:02:49,512 **init** INFO] Module [seletop] finished.  
[2024-02-21 18:02:49,512 **init** INFO] [seletop] took 0 seconds  
[2024-02-21 18:02:49,603 **init** INFO] Running [clustrmsd] module  
[2024-02-21 18:02:49,605 libutil ERROR] list index out of range  
I keep the default parameters, and this is how I have them set :  
[clustrmsd]  
criterion = ‘maxclust’  
linkage = ‘average’  
threshold = 4

---

<div class="post-metadata">

**Author:** ![marco.giulini](https://dub1.discourse-cdn.com/flex013/user_avatar/ask.bioexcel.eu/marco.giulini/32/961_2.png) [@marco.giulini](https://ask.bioexcel.eu/u/marco.giulini)\
**Post date:** [February 22, 2024, 3:31pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/2 "2024-02-22T15:31:25Z")

</div>

hi,

you should add a matrix calculation module before clustrmsd, like rmsdmatrix or ilrmsdmatrix. something like this in your workflow:

```auto
[seletop]

[ilrmsdmatrix]

[clustrmsd]

```

My advise is to use rmsdmatrix if you have a good knowledge of the interface, as you can specify it using the resdic parameter (an example here [haddock3/examples/docking-protein-glycan/docking-protein-glycan-test.cfg at main · haddocking/haddock3 · GitHub](https://github.com/haddocking/haddock3/blob/main/examples/docking-protein-glycan/docking-protein-glycan-test.cfg) ). if your knowledge of the interface is not great then ilrmsdmatrix is a better option (see [haddock3/examples/docking-protein-glycan/docking-protein-glycan-ilrmsd-test.cfg at main · haddocking/haddock3 · GitHub](https://github.com/haddocking/haddock3/blob/main/examples/docking-protein-glycan/docking-protein-glycan-ilrmsd-test.cfg)).

Hope it helps!

---

<div class="post-metadata">

**Author:** ![YFeri](https://avatars.discourse-cdn.com/v4/letter/y/b4bc9f/32.png) [@YFeri](https://ask.bioexcel.eu/u/YFeri)\
**Post date:** [February 22, 2024, 3:58pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/3 "2024-02-22T15:58:08Z")

</div>

Hi Marco,  
thanks for getting back to me quickly. I’m currently working on ab-initio protein-protein docking, so I don’t have any information about the interface or interaction surface. I tried checking the ilrmsdmatrix module by running “haddock3-cfg -m ilrmsdmatrix,” but I couldn’t find the module. Is it available on HADDOCK3? Can caprieval provide similar metrics without needing a reference structure? As I mentioned, I’m running in ab-initio mode.

---

<div class="post-metadata">

**Author:** ![marco.giulini](https://dub1.discourse-cdn.com/flex013/user_avatar/ask.bioexcel.eu/marco.giulini/32/961_2.png) [@marco.giulini](https://ask.bioexcel.eu/u/marco.giulini)\
**Post date:** [February 22, 2024, 4:28pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/4 "2024-02-22T16:28:03Z")

</div>

that means that you have to update your HADDOCK codebase (see [haddock3/docs/INSTALL.md at main · haddocking/haddock3 · GitHub](https://github.com/haddocking/haddock3/blob/main/docs/INSTALL.md#5-keep-your-installation-up-to-date)). In the case of ab-initio docking the ilrmsdmatrix is indeed the best choice.

as the code is still evolving fast check also the parameter names once you update the codebase (we changed the names of some of them lately)

---

<div class="post-metadata">

**Author:** ![YFeri](https://avatars.discourse-cdn.com/v4/letter/y/b4bc9f/32.png) [@YFeri](https://ask.bioexcel.eu/u/YFeri)\
**Post date:** [February 26, 2024, 5:25pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/5 "2024-02-26T17:25:53Z")

</div>

Hi @marco.giulini ,  
can I still use haddock3 for my docking ever if the code hasn’t finished yet ? will the code still has a major modification ? does the available code sufficient to use haddock3 ?  
Thanks to you for all the responses that you have given me.

---

<div class="post-metadata">

**Author:** ![amjjbonvin](https://dub1.discourse-cdn.com/flex013/user_avatar/ask.bioexcel.eu/amjjbonvin/32/23_2.png) [@amjjbonvin](https://ask.bioexcel.eu/u/amjjbonvin)\
**Post date:** [February 26, 2024, 6:36pm UTC](https://ask.bioexcel.eu/t/error-in-haddock3-clustering-rmsd-based-clustrmsd-during-protein-protein-docking/4857/6 "2024-02-26T18:36:16Z")

</div>

The code is still evolving, but close to the first production release.

Simply make sure to update your installation in case you run into problems.
