You are here

The question of ligand docking with constraints

3 posts / 0 new
Last post
The question of ligand docking with constraints


I read the paper Rosetta and the Design of Ligand Binding Sites

and it enlightens me.

According to the method of this paper,  I design my enzyme.

However,  I find these are a few unreasonable results which show reaction atom of ligand  do not close to  catalytic triad,

but  the opposite side of ligand close to  catalytic triad.

So I want to  make docking with  constraints

I think it can achieved with ConstraintSetMover, but  do not understand how write the Constraint File.


For example,  I want to set the distance betweent  atom CA of  resid 516 (one of catalytic triad)  and C1 of ligand, and the distance is 3 ± 0.5 A .

AtomPair  516 CA   585  C1 HARMONIC  Func_Def 

I  read  the  equation f(x)=[(x-x0)/sd] ^2 , but I do not konw how to set  3 ± 0.5 A in HARMONIC function.

I think " AtomPair  516 CA   585  C1 HARMONIC  3  0.5 " is wrong

Could  give a right sample that how  to set distance  (3 ± 0.5 A)?




Post Situation: 
Mon, 2018-06-11 09:39

Since you found the harmonic equation I assume you found the constraints doc page:

You can't set a "range" of 3 +- 0.5 with harmonic.  It's a parabola, it doesn't have a range, it is defined over all possible X.  You can play with the SD parameter such that the parabola is narrower or wider, which will have the effect of disfavoring  / being permissive of distances further from 3 angstrom, but you can't forbid values outside the 2.5-3.5 range.

Square well functions don't work terribly well with Rosetta because they are unminimizable, but reading the doc page I do see FADE which may be of some use.  Find FADE on the page and see if meets your needs better?

I always encourage people to plot out constraint energy value vs. x deviation in the plotting software of their choice to get a feel for how strong a constraint they are creating.  Remember that "good residues" are in the range of -3 to -2 energy units.

Mon, 2018-06-11 19:27

Thanks for your suggestion.

Function of constraint energy is complicated.

I will have a try  ~~

Thu, 2018-06-14 04:01