Date: Wed, 30 Mar 94 11:01:42 PST From: case@scripps.edu (David Case) To: amber@cgl.ucsf.EDU Subject: Re: Constraints between moleculesThe restraints in Sander don't care about molecules if you use atom numbers, so try the following:
&rst iat =This will make the constraint. You can use the NSTEP1,NSTEP2 variables to have it active for a while, and then be turned off., , r1=??. r2=??.....etc
From: dap@portal.vpharm.com (David A. Pearlman) Subject: Re: Constraints between moleculesHere's sample input which would do the following:
COM2 is defined by:
C4 in residue 1
O4 in residue 1
N* in residue 1
(These definitions are effected by the igr1/igr2 and grnam1/grnam2
variables; You can use up to 200 atoms to define a center-of-mass
group)
&rst iat=-1,-1,1,1, atnam(1)='DUMM',atnam(2)='DUMM', nstep1=1,nstep2=5000, iresid=1,irstyp=0,ifvari=1,ninc=0,imult=0,ir6=0,ifntyp=0, r1=0.00000E+00,r2=5.0000,r3=5.0000,r4=99.000,rk2=1.0000,rk3=1.0000, r1a=0.00000E+00,r2a=2.5000,r3a=2.5000,r4a=99.000,rk2a=1.0000,rk3a=1.0000, igr1 = 2,3,4,5,0, grnam1(1)='C1',grnam1(2)='C1',grnam1(3)='N2',grnam1(4)='C1', igr2 = 1,1,1,0, grnam2(1)='C4',grnam2(2)='O4',grnam2(3)='N*', &end &rst iat=-1,-1,1,1, atnam(1)='DUMM',atnam(2)='DUMM', nstep1=5001,nstep2=10000, iresid=1,irstyp=0,ifvari=1,ninc=0,imult=0,ir6=0,ifntyp=0, r1=0.00000E+00,r2=2.5000,r3=2.5000,r4=99.000,rk2=1.0000,rk3=1.0000, r1a=0.00000E+00,r2a=2.5000,r3a=2.5000,r4a=99.000,rk2a=1.0000,rk3a=0.0100, igr1 = 2,3,4,5,0, grnam1(1)='C1',grnam1(2)='C1',grnam1(3)='N2',grnam1(4)='C1', igr2 = 1,1,1,0, grnam2(1)='C4',grnam2(2)='O4',grnam2(3)='N*', &end &rst iat=0, &end