Hello Santhos--

> 
> My concern is this :
> 
> The force on the conformation whose weight is 0.1 is too low to drive it to
> the minimum. The possible solution is that increase the force constant by
> using rdc.setScale or something similar but that would mean the conformation
> with weight 0.9 would have a large force and am not sure if that would not
> cause any problems. The multiplication of the weight factor might cause
> problems in setting the right force constant. 

The prefactors for each ensemble contribution is what it is, and
changing that in the gradient calculation will only cause a worse
minimum to be obtained (because the gradient would be incorrect).

> 
> While what is being done is mathematically correct, in my
> experiments the protocol did not work well. We took 2 distinct
> conformations of a protein and assigned them weights of 0.8, 0.2
> while generating the artificial RDC's and NOEs from those. When
> starting with the global minimum [the original 2 conformations  with
>  weights of 0.8, 0.2], the second member moved away from the global
> minimum. [I hope I was clear in my explanation, please ask if
> something is not clear].
> 

ok. That means that something is wrong or not accounted for. I'm happy
to help if I know more details. You might upload a complete example
which shows the poor behavior here (after logging in):

  http://nmr.cit.nih.gov/xplor-nih/download.cgi

best regards--
Charles

_______________________________________________
Xplor-nih mailing list
Xplor-nih@cake.cit.nih.gov
http://cake.cit.nih.gov/mailman/listinfo/xplor-nih

Reply via email to