[gmx-users] grompp -r again

hugo verli hugo at acd.ufrj.br
Wed Oct 23 16:05:46 CEST 2002


David,

I follow your suggestion to use free energy perturbation to transform one
conformation in another conformation. I fix the boat conformation using improper
dihedrals.
To test if these improper works I simulate a chair conformation with the boat
topology and the transformation occur quickly, in a 50ps simulation. So I go to
the FEP calculations, using Dlambda step of 0.05 for 50ps. However in the end of
the simulation I did not obtain the the boat conformation. The ring adopt a
distorted conformation very different from the intented boat conformation.

Did you understand what is happening? Am I doing something wrong? I tested
different combinations of improper dihedrals and I see the same profile. How can
I fix it?

Many thanks in advance,

Hugo Verli.

> On Thu, 2002-10-03 at 17:40, hugo verli wrote:
> > Thanks, David.
> > 
> > But I am using for now a system of just one sugar ring to test the
> conversion.
> > However, I need to apply this transformation to a polyssacharide
> system
> > complexed to one protein, so the system syze will became prohibitive
> to free
> > enegy perturbation, dont'it?
> Not necessarily. If you only want to modify a single ring in the whole
> sugar chain it should still work, at least if the binding mode remains
> the same.
> > Besides the two conformations coexist in solution, so the energy
> barrier between
> > them should not be soo high. 
> Well, that both exist means that the free energy difference is not very
> large; the barrier may be in principle, although it doesn't seem too
> likely in this case.
> 
> 
> -- 
> Groeten, David.
> ________________________________________________________________________
> Dr. David van der Spoel, 	Biomedical center, Dept. of Biochemistry
> Husargatan 3, Box 576,  	75123 Uppsala, Sweden
> phone:	46 18 471 4205		fax: 46 18 511 755
> spoel at xray.bmc.uu.se	spoel at gromacs.org   http://zorn.bmc.uu.se/~spoel
> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
> _______________________________________________
> 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.
> 



More information about the gromacs.org_gmx-users mailing list