I'm not familiar with the restraint definition that you used, are you
following the manual
or one of the tutorials? I don't see any use of "iatms" in the manual. Check
page
178 in the Amber9 manual for correct variables.
I also do not think Amber supports dihedral restraints using center
of mass, only distances between centers of mass.
Have you already used Amber to perform normal (non-PMF) simulations?
You'll want to make sure everything is set up properly before you start
trying
restraints.
good luck
On 5/8/07, Catein Catherine <askamber23.hotmail.com> wrote:
>
> Dear All,
>
> I am trying to calculate the PMF for a system, with the four domains:
> Domain 1 to domain 4.
>
> Domain1---Domain2
> |
> Domain4---Domain3
>
>
> I tried to define the in file as follows (with the center of mass defined
> for each domain):
>
> &wt type='DUMPFREQ', istep1=10 /
> &wt type='END' /
> DISANG=xxx.RST
> DUMPAVE=xxx_vs_t.10
>
> =====xxx.RST to study the Gibbs free energy vs different dihedral angle of
> different domains======
>
> # Dihedral restrainst
> &rst iatms(1)=1/33, iatms(2)=34/64, iatms(3)=65/94, iatms(4)=95/126,
> iatms(5)=127/177, iatms=1,2,3,4, r1=0., r2=90., r3=180., r4=360., rk2 =
> 30.,
> rk3= 30.0, /
>
> =====xxx.RST to study the Gibbs free energy vs different distance of
> domain
> 1 and 2========
>
> # Distance restrainst
> &rst iatms(1)=1/33, iatms(2)=34/64, iatms=1,2, r1=0., r2=2., r3=5.,
> r4=10.,
> rk2 = 30., rk3= 30.0, /
>
>
> =============================================================================
>
> It is the right way to define the RST and in files for the PMF
> calculations?
> Did I miss a lot of things here, as I cannot get the xxx_vs_t.10 files
> after submit the jobs?
>
> Best regards and many thanks!
> Cat.
>
> _________________________________________________________________
> Get 10Mb extra storage for MSN Hotmail. Subscribe Now!
> http://join.msn.com/?pgmarket=en-hk
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed May 09 2007 - 06:07:27 PDT
