Dear Dr. Case,
Thank you very much for your prompt reply.
I am trying to put constraints such that restraints are put such that symmetrical groups are free rotate but the main group cannot move away from its position.
My restrained mask in my minimization script reads as follows :
Initial Structural Minimization with explicit water
&cntrl
imin=1, maxcyc=5000, ncyc=500, ntmin=1,
ntb=1, cut=15, ntpr=10, ntr=1,
RESTRAINT_WT=1.0,
RESTRAINTMASK="(:2-150) & ! (:PHE.CE=,CD= | :TYR.CD=,CE= | :ASP.OD= | :GLU.OE= | :ARG.NH= | :ASN.ND2,OD1 :GLN.NE2,OE1 | :HIE.ND1,CD2,CE1,NE2 | :HID.ND1,CD2,CE1,NE2 | :HIP.ND1,CD2,CE1,NE2 | .H= | :150.O)",
/
In my mini.out file, the restrained mask (in bold) I see that is being processed is only the first 80 characters or so :
-------------------------------------------------------
Amber 9 SANDER 2006
-------------------------------------------------------
| Run on 03/09/2009 at 14:52:07
[-O]verwriting output
File Assignments:
| MDIN: mini.in
| MDOUT: output/mini.out
|INPCRD: bak.inpcrd
| PARM: bak.prmtop
|RESTRT: output/mini.rst
| REFC: bak.inpcrd
| MDVEL: mdvel
| MDEN: mden
| MDCRD: output/mini.mdcrd
|MDINFO: output/mdinfo
|INPDIP: inpdip
|RSTDIP: rstdip
|INPTRA: inptraj
|
Here is the input file:
Initial Structural Minimization with explicit water
&cntrl
imin=1, maxcyc=5000, ncyc=500, ntmin=1,
ntb=1, cut=15, ntpr=10, ntr=1,
RESTRAINT_WT=1.0,
RESTRAINTMASK="(:2-23,33-150) & ! (:PHE.CE=,CD= | :TYR.CD=,CE= | :ASP.OD= | :
/
As you see, I am using the Sander module available in Amber 9.
Thanks once again,
Cheers,
Sampath.
----- Original Message ----
> From: David A. Case <case.biomaps.rutgers.edu>
> To: AMBER Mailing List <amber.ambermd.org>
> Sent: Monday, March 9, 2009 2:43:34 PM
> Subject: Re: [AMBER] Restraint Mask length restriction !!
>
> On Mon, Mar 09, 2009, Sampath Koppole wrote:
> >
> > Is there a way to work around the "restraint mask length being limited
> > to 80 characters" limitation ? I have a complex a complex restrained
> > mask which is longer than 80 characters (~120 characters !!).
>
> The restraintmask limit is currently 256 characters in the code. But
> there might be a namelist limitation(?). Can you give details (exact
> input you are using, which program and which version of Amber)?
>
> ...thx...dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 11 2009 - 01:10:33 PDT
