AMBER: details of new aminoacid

From: Gustavo Pierdominici Sottile <gsottile.unq.edu.ar>
Date: Tue, 18 Jan 2005 10:01:15 -0300

Hi, I send the details of the problems with the developing of a new
aminoacid. Suppose I have to designate a new amino acid. I create the
archive .in that looks like this (it is not complete)
  0 0 2
POCOIMAGIN
trp1.db4
 AVO INT 1
 CORR OMIT DU BEG
   0.00000
   1 DUMM DU M 0 -1 -2 0.000 0.000 0.000 0.00000
   2 DUMM DU M 1 0 -1 1.449 0.000 0.000 0.00000
   3 DUMM DU M 2 1 0 1.522 111.100 0.000 0.00000
   4 N N M 3 2 1 1.335 116.600 180.000 -0.41570
   5 H H E 4 3 2 1.020 119.800 0.000 0.27190
   6 CA CT M 4 3 2 1.470 121.900 180.000 -0.05270
   7 HA H1 E 6 4 3 1.090 107.300 300.000 0.08660
   8 CB CT 3 6 4 3 1.550 110.500 60.000 0.03550
   9 HB2 HC E 8 6 4 1.090 108.400 173.900 0.05730
  10 HB3 HC E 8 6 4 1.090 109.800 57.500 0.05730
  11 CG C* S 8 6 4 1.510 112.600 292.400 -0.09950
  12 CD1 CW B 11 8 6 1.370 129.100 126.700 -0.05060
  13 NE1 NA B 12 11 8 1.390 107.100 190.000 -0.31550
  14 HE1 H E 13 12 11 1.020 121.700 167.400 0.31710
  15 CE2 CN S 13 12 11 1.380 113.200 350.600 0.07850
  16 CZ2 CA B 15 13 12 1.400 133.900 183.700 -0.15890


Then I execute the command tleap:
tleap -f leaprc.ff94
Once it is executed I write
loadamberprep "filename" (the name of the file which part is shown above)
mad= loadpdb "name of the pdb file"
when I write saveamberparm mad top.l crd.l, the following lines appears

Could not find bond parameter for: CW - C
Could not find angle parameter: C* - CW - C

The problem is that there are no atoms of this type joined in the new
aminoacid. I do not have two atoms of type CW and C joined. The same happens
with angles and dihedral.
It also appears one warning for each water molecule. All of then states the
same:

WARNING: There is a bond of 3.464102 angstroms between:
------- .R<WAT 7341>.A<H2 2> and .R<WAT 7341>.A<O 3>



This is not true because all the bonds in water molecule are smaller then 1
angstroms.


Thanks in advance
Gustavo


   

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Jan 18 2005 - 13:53:00 PST
Custom Search