# Restraints violation

**URL:** <https://ask.bioexcel.eu/t/restraints-violation/4266>\
**Category:** HADDOCK\
**Created:** [April 16, 2023, 11:59am UTC](https://ask.bioexcel.eu/t/restraints-violation/4266 "2023-04-16T11:59:04Z")\
**Posts on this page:** 10\
**Page:** 1

<div class="post-metadata">

**Author:** ![Miriam](https://avatars.discourse-cdn.com/v4/letter/m/c89c15/32.png) [@Miriam](https://ask.bioexcel.eu/u/Miriam)\
**Post date:** [April 16, 2023, 11:59am UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/1 "2023-04-16T11:59:04Z")

</div>

Hi,  
I want to dock a peptide and protein.  
I have some ambiguous restraints that I see from the NOESY NMR spectrum.  
But for some reason the HADDOCK results do not match the restraints and the restraints violation in the result page equals 0.  
How can I fix this?  
I didn’t add any unambiguous restraints and made sure it calculates all hydrogens.  
Also, how many iterations is recommended?

my restraints are:  
assign(segid A and resid 24 and name HD)(segid B and resid 3 and name HG\*) 3.5 2.0 2.5  
assign(segid A and resid 24 and name HD)(segid B and resid 3 and name HD1) 3.5 2.0 2.5  
assign(segid A and resid 24 and name HD)(segid B and resid 3 and name HD2) 3.5 2.0 2.5  
assign(segid A and resid 68 and name HD)(segid B and resid 3 and name HD1) 3.5 2.0 2.5  
assign(segid A and resid 68 and name HD)(segid B and resid 3 and name HD2) 3.5 2.0 2.5  
assign(segid A and resid 52 and name HH)(segid B and resid 5 and name HG\*) 3.5 2.0 2.5  
assign(segid A and resid 52 and name HH)(segid B and resid 5 and name HB) 3.5 2.0 2.5  
assign(segid A and resid 24 and name HE)(segid B and resid 3 and name HD1) 3.5 2.0 2.5  
assign(segid A and resid 24 and name HE)(segid B and resid 3 and name HD2) 3.5 2.0 2.5  
assign(segid A and resid 52 and name HZ3)(segid B and resid 5 and name HG\*) 3.5 2.0 2.5  
assign(segid A and resid 52 and name HZ2)(segid B and resid 5 and name HG\*) 3.5 2.0 2.5  
assign(segid A and resid 52 and name HZ2)(segid B and resid 5 and name HB) 3.5 2.0 2.5  
assign(segid A and resid 48 and name HD\*)(segid B and resid 11 and name HZ2) 33.5 2.0 2.5  
assign(segid A and resid 48 and name HD\*)(segid B and resid 11 and name HZ3) 3.5 2.0 2.5  
assign(segid A and resid 48 and name HD\*)(segid B and resid 11 and name HH2) 3.5 2.0 2.5

Is this not enough?

Thx,  
Miriam

---

<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:** [April 16, 2023, 1:18pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/2 "2023-04-16T13:18:23Z")

</div>

Two things:

1. Turn on the option to keep all hydrogens as by default all non-polar hydrogens are removed

2. Make sure that the naming matches the naming in the PDB file

Or as solution to both points: define the restraints to the heavy atoms (with a small correction on the upper limit.

---

<div class="post-metadata">

**Author:** ![Miriam](https://avatars.discourse-cdn.com/v4/letter/m/c89c15/32.png) [@Miriam](https://ask.bioexcel.eu/u/Miriam)\
**Post date:** [April 16, 2023, 2:12pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/3 "2023-04-16T14:12:03Z")

</div>

Thx!  
how can I define the restraints to the heavy atoms with a small correction on the upper limit?

---

<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:** [April 16, 2023, 4:17pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/4 "2023-04-16T16:17:23Z")

</div>

Simply change the atom name to the heavy atom to which the hydrogen it attached to.  
And add 1A to the upper limit

---

<div class="post-metadata">

**Author:** ![Miriam](https://avatars.discourse-cdn.com/v4/letter/m/c89c15/32.png) [@Miriam](https://ask.bioexcel.eu/u/Miriam)\
**Post date:** [April 28, 2023, 2:43pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/5 "2023-04-28T14:43:37Z")

</div>

Hi Thank you  
I’ve tried this many times and it still does nit work and some of the restraints are not kept.  
is there a problem with my format? maybe the \* is confusing HADDOCK? spaces? or any other reason?  
I cant think of anything else.  
Here are the list when the are between carbon atoms:  
assign(segid A and resid 24 and name CD\*)(segid B and resid 3 and name CG\*) 4.5 2.0 2.5  
assign(segid A and resid 24 and name CD\*)(segid B and resid 3 and name CD1) 4.5 2.0 2.5  
assign(segid A and resid 24 and name CD\*)(segid B and resid 3 and name CD2) 4.5 2.0 2.5  
assign(segid A and resid 68 and name CD\*)(segid B and resid 3 and name CD1) 4.5 2.0 2.5  
assign(segid A and resid 68 and name CD\*)(segid B and resid 3 and name CD2) 4.5 2.0 2.5  
assign(segid A and resid 52 and name CH\*)(segid B and resid 5 and name CG\*) 4.5 2.0 2.5  
assign(segid A and resid 52 and name CH\*)(segid B and resid 5 and name CB) 4.5 2.0 2.5  
assign(segid A and resid 24 and name CE\*)(segid B and resid 3 and name CD1) 4.5 2.0 2.5  
assign(segid A and resid 24 and name CE\*)(segid B and resid 3 and name CD2) 4.5 2.0 2.5  
assign(segid A and resid 52 and name CZ3)(segid B and resid 5 and name CG\*) 4.5 2.0 2.5  
assign(segid A and resid 52 and name CZ2)(segid B and resid 5 and name CG\*) 4.5 2.0 2.5  
assign(segid A and resid 52 and name CZ2)(segid B and resid 5 and name CB) 4.5 2.0 2.5  
assign(segid A and resid 48 and name CD\*)(segid B and resid 11 and name CZ2) 4.5 2.0 2.5  
assign(segid A and resid 48 and name CD\*)(segid B and resid 11 and name CZ3) 4.5 2.0 2.5  
assign(segid A and resid 48 and name CD\*)(segid B and resid 11 and name CH2) 4.5 2.0 2.5

---

<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:** [April 29, 2023, 5:31pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/6 "2023-04-29T17:31:24Z")

</div>

Can you share the link to the run, so that we can check if all restraints are properly read in.

If this is the case it could be you have some inconsistent restraints or more conformational changes are required

---

<div class="post-metadata">

**Author:** ![Miriam](https://avatars.discourse-cdn.com/v4/letter/m/c89c15/32.png) [@Miriam](https://ask.bioexcel.eu/u/Miriam)\
**Post date:** [April 29, 2023, 6:06pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/7 "2023-04-29T18:06:02Z")

</div>

yes, thank you!  
I ran a lot of them (with both restraints on hydrogens and carbons), attached are two different structures I tried to dock (they have different restraints):  
[https://wenmr.science.uu.nl](https://wenmr.science.uu.nl)

[https://wenmr.science.uu.nl](https://wenmr.science.uu.nl)

Thank you!

---

<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:** [April 29, 2023, 6:56pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/8 "2023-04-29T18:56:42Z")

</div>

The links are to the home page of our server and not the specific runs…

---

<div class="post-metadata">

**Author:** ![Miriam](https://avatars.discourse-cdn.com/v4/letter/m/c89c15/32.png) [@Miriam](https://ask.bioexcel.eu/u/Miriam)\
**Post date:** [April 29, 2023, 7:08pm UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/9 "2023-04-29T19:08:47Z")

</div>

Sorry about that  
[https://wenmr.science.uu.nl/haddock2.4/result/7747164191/234918-QKC\_1H\_2804](https://wenmr.science.uu.nl/haddock2.4/result/7747164191/234918-QKC_1H_2804)

[https://wenmr.science.uu.nl/haddock2.4/result/7747164191/234916-QKC\_2H\_gui](https://wenmr.science.uu.nl/haddock2.4/result/7747164191/234916-QKC_2H_gui)

---

<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:** [May 1, 2023, 11:01am UTC](https://ask.bioexcel.eu/t/restraints-violation/4266/10 "2023-05-01T11:01:34Z")

</div>

To check that your restraints are properly read, edit e.g. the `complex_run1_it0_refine_1.out` file which you can find after downloading and unpacking the full archive of the run.

Then search for: `NOE>assign`

and check that atoms are selected for all your restraints. A quick check of the QKC\_1H\_2804 runs gives the following problems:

> assign(segid A and resid 24 and name HE\*)(segid B and resid 11 and name HG\*) 3.5 2.0 2.5  
> SELRPN: 2 atoms have been selected out of 1100  
> SELRPN: 0 atoms have been selected out of 1100

> assign(segid A and resid 48 and name HD\*)(segid B and resid 1 and name HH2) 3.5 2.0 2.5  
> SELRPN: 3 atoms have been selected out of 1100  
> SELRPN: 0 atoms have been selected out of 1100

Check the atom naming for those two restraints.

Further the restraint energy of your final models is rather low.
