[gmx-users] Simulated annealing
kay.gottschalk at weizmann.ac.il
Wed Jul 30 21:26:02 CEST 2003
It might very well be that your system is way to flexible; you'd need
experimental (=NMR) restraints, I'd guess. Even a cyclic hexapeptide is
very flexible and basically improbable to predict correctly.
On Thursday, July 31, 2003, at 01:15 AM, Osmany Guirola Cruz wrote:
> my peptide have 11 aminoacids and it is cyclic.
> Marc Ceruso wrote:
> On Wed, 30 Jul 2003, Osmany Guirola Cruz wrote:
> Thanks again Xavier,
> I have this problem:
> I have a peptide whit unknown structure, whit "modeler" i get 10
> diferents structures of my peptide, and then.....
> i use SA whith the ten models to obtain a "superminimized structure"
> my final structures are NOT similar
> How long is the peptide? Why do you think that your structures should
> similar? In any case, if the peptide is small let's say 1-40 amino
> chances are that its structure in solution is not unique, but in any
> the similiraty within your ensemble will depend on the strength of the
> restraints used to generate the structures. If the homology of your
> peptide/protein is low the constraints will be loose reflecting the
> of empirical knowlege regarding your sequence, in others reflecting the
> uncertainty in the structure.
> :-( . i have the idea that my simulations converge to a unique
> (my SA = "1000K" to 300K) you say 10 000 i will probe
> with 10000 to see what happens. ???????
> Xavier Periole wrote:
> I guess you have been reading papers from NMR procedures to
> generate their models corresponding to the data. If that what you
> are doing It would be easier to get a program that NMR people
> use, is is probably implemented
> It seeems difficult to do that automatically in gromacs, at least
> for me.
> I guess you that could modify the strength of the constant on the
> and bonds in the parameter file and do that for each temperature
> need but that's gonna to be a pain in the ...
> I don't know what you are doing but I can assure you that at 10000
> the system is pretty flexible. Try to simulate it for 10 ps at
> 10000 K
> and look at the trajectory !! It should convice you. But I don't
> know if
> it is relevant for your problem.
> _______________________________________________ gmx-users mailing list
> gmx-users at gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the www
> interface or send it to gmx-users-request at gromacs.org.
> gmx-users mailing list
> gmx-users at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
Dr. Kay-E. Gottschalk
Department of Biological Chemistry
Weizmann Institute of Science
Herzl St. 1
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Size: 3473 bytes
Desc: not available
More information about the gromacs.org_gmx-users