Dear amber users,
I have read tutorial concerning use of gibbs to perform free energy of 2
states; Valine to Alanine, at
http://web.njit.edu/all_topics/Prog_Lang_Docs/html/amber7/gibbs_box/val-to-ala.html#xleap_done
This is free energy between 2 different molecules. However, I would like to
calculate the free energy of dimer with 2 conformations in explicit solvent;
(I). state 0 := monomer1 is on top of monomer2 as parallel conformation
(II). state 1 := monomer1 is on top of monomer2, but twist angle between 2
monomers is 180 degree (antiparallel conformation)
As I understand, I should do as following
1). Load pdb file of state 0 to xleap
2). Edit state 0 and perturb. Since both monomers are identical, thus
PERT.name, PERT.type and PERT.charge should be the same as unperturb state.
Am I right? If yes, topology files of unperturb and perturb states should be
the same.
3). Minimization and dynamics the unperturb system to equilibrate solvent.
4). Run GIBBS
To run gibbs, I have to use restrt file from step (3) and topology file of
perturb state which is the same as unperturb state. Because pdb file of
antiparallel conformation has never used, how can I specify the different
between these two states in GIBBS?
Many thanks in advance.
Sincerely Yours,
Xioling
_________________________________________________________________
Add photos to your messages with MSN 8. Get 2 months FREE*.
http://join.msn.com/?page=features/featuredemail
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Jun 18 2004 - 14:53:00 PDT