Hi Frido,
Indeed we are also in need of similar stuff. Davide Bau (a postdoc in my lab and also in the IMP-Dev list) has already implemented a sort of a MORSE potential in IMP restraints, which is not what you ask but may be similar. I implemented many many many years ago something similar in CHARMM (http://sgu.bioinfo.cipf.es/martirenom/Thesis/pdf/ Thesis.pdf, Chapter 4 page 86).
Unfortunately, we are NOT yet ready to commit the changes in the code to the SVN since we are in the middle of production phase (ie, making models) thus would prefer not to update IMP. But I am sure we could do soon.
Davide, any comments? :-)
m
------- Marc A. Marti-Renom, Head of the Structural Genomics Unit Bioinformatics & Genomics Department, Prince Felipe Research Center Avda. Autopista del Saler, 16 , 46012 Valencia, Spain Tel +34 96 328 96 80 (ext. 1015) Fax +34 96 328 97 01 email mmarti@cipf.es web http://sgu.bioinfo.cipf.es
*** Consider the environment. Print responsibly. ***
On Jan 8, 2009, at 10:00 AM, Friedrich Foerster wrote:
> hi all, > > when dealing with experimental restraints i often feel that harmonic > restraints are somewhat unsatisfactory: > experimental data can be wrong and thus there is not necessarily a > solution that fulfills all restraints (score=0). in those cases, a > model with a lower score is not necessarily any more correct than a > higher scoring one. thus, i think it'd be good to have an additional > scoring function that levels off if the feature exceeds a certain > value. for example, i'd suggest the following two functions > alternatively to Harmonic: > 1. a function that is a Harmonic for x<std and constant for x>std > 2. a function that is a Harmonic for x<std and asymptotically > approaches a constant value for x-> inf: > e.g.: f=x^2 for x<1 and f=2-x^-2 for x>1 (would even be > continuous in 1st derivative) > > is anyone else also interested in something like that? > any better suggestions? > is any IMP core expert interested and willing to code that > (including analogous functions for Upper and Lower Harmonic)? > probably i'd be one of the slowest and clumsiest persons in doing > that... > > thanks > > frido > > > -- > > Friedrich Foerster > Max-Planck Institut fuer Biochemie > Am Klopferspitz 18 > D-82152 Martinsried > > Tel: +49 89 8578 2651 > Fax: +49 89 8578 2641 > > foerster@biochem.mpg.de > > www.tomotronic.org > > > > > > > _______________________________________________ > IMP-dev mailing list > IMP-dev@salilab.org > https://salilab.org/mailman/listinfo/imp-dev >