Re: AMBER: drawing repeated units

From: Taufik Al-Sarraj <taufik.alsarraj.utoronto.ca>
Date: Tue, 24 Jun 2008 17:15:39 -0500

Hi,

I tried this few times and it works up to saveamberparm. the error i get
is could not find angle parameter: etc etc
and No torson terms for ....
i built the units properly (using antechamber) but when two units are
connected a new bond is created along with new angles and torsions, how
can describe those? should i set them manually?

I know the mol = sequence { FPR REPU REPU REPU REPU ENDU} is working
because when i type edit mol i see the chain as expected.

I am using amber9 and i usually use xleap -s -f leaprc.ff99SB

Thanks,
Taufik

Ross Walker wrote:
> Hi Taufi,
>
>> Thank you Ross,
>> It worked!
>> Where can i learn to create a simple leaprc file. is it in any of the
>> tutorials?
>
> A leaprc file is just a text file containing the commands exactly as you
> would enter them in leap. E.g. one could be:
>
> source leaprc.ff99SB
> loadamberparams my.frcmod
> loadoff foo.lib
> loadoff bar.lib
> mol = sequence { FOO BAR FOO BAR }
> saveamberparm mol prmtop inpcrd
> quit
>
> Then you could just run this as:
>
> tleap -s -f my.leaprc
>
> and it would just run those commands and quit.
>
> You can also do:
>
> source my.leaprc
>
> to read it from within an active leap session.
>
>> i can also type mymol = sequence { FPRT REPU ... 50 100 150 and 200
>> times.. ENDU } if writing a one -long- line script is too difficult.
>
> You can't do that in leap unfortunately but you could either write a simple
> perl or shell script that writes leaprc files for you - i.e. writes those
> long lines for you. Or you could build some bigger units. I.e. build a unit
> of 10 REPU's - so you would only have to put 10 of these units in to get a
> 100 repeat.
>
> Good luck,
> Ross
>
>
> /\
> \/
> |\oss Walker
>
> | Assistant Research Professor |
> | San Diego Supercomputer Center |
> | Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
> | http://www.rosswalker.co.uk | PGP Key available on request |
>
> Note: Electronic Mail is not secure, has no guarantee of delivery, may not
> be read every day, and should not be used for urgent or sensitive issues.
>
>
>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo.scripps.edu


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu

Received on Wed Jun 25 2008 - 06:07:53 PDT
Custom Search