AMBER Archive (2007)

Subject: Re: AMBER: perturbed atoms are not bonded on using TIP5P water

From: sugino tatsuro (tsugino_at_comas.frsc.tsukuba.ac.jp)
Date: Mon Jan 15 2007 - 19:40:32 CST


Thank you very much for your advices.

>How did you establish this? Were you looking inside the prmtop file to
>see that the perturbed water had no bonds, but all the others did?
>How did you distinguish the one atom to be perturbed from the others?
>(A different unit name?)

I looked at the inside of the parm.top file by using 'rdparm' to see bond
information. Although perturbed water molecule had no bonds, all the other
water molecules did.
And I distinguished
perturbed one by using VMD and looked it by using movie implemented in VMD.

Here, topology and coordinate files were established by using the following
commands;
> loadamberparams gaff.dat
> source TIP5P.cmd
> loadamberparams parm99EP.dat
> loadamberparams frcmod.tip5p
> model = TP5
> solvatebox model TP5 20
> edit model ( and select a water molecule and set perturbation by GUI )
> saveamberparm model parmpert.top parmpert.crd

I renamed the residue name of the perturbed water / ligand molecules, but
the bonds inside the molecules were not formed; I don't understand what
occurs in the process of the file generation, in which TIP5P water was used
for the water model (in the case of TIP3P, no problems occur).
How can I make bonds among perturbed atoms when TIP5P water model was used?

Any advices are welcome. Thank you very much in advance.

best wishes,
Sugino

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu