This is a MIME-encapsulated message
--==_10300106002765.ams1.ulsop.ac.uk==_
Content-Type: text/plain; charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
Dear all,
I need to use leap to remove a large number of atoms from a structure,=
before reading in another PDB file and generating prmtop and prmcrd file.=
(I'm creating a ligand covalently bound to DNA).=20
I would like to do this in a script (leap.in file), but I have found that
leap crashes with a segmentation fault if I try to remove an atom (using=
'remove' text command) when it is bonded to any other atom. If I remove=
all the relevant bonds ('deletebond') first, so that the atom is isolated,=
then the 'remove' command works. But identifying all the bonds correctly=
is a bit tedious, and not easily transferrable to other molecules, which=
was the point of using a script.
I guess that this is a bug otherwise there would be no seg fault.=
Presumably it is the desired behaviour that removing an atom also removes=
all bonds to it. If anyone can suggest an alternative approach, or is=
interested in trying to fix the bug, then I'd be very grateful.
I am using AMBER7 on RedHat 9, compiled with ifort8+gcc3.2.2
Dave Evans
London School of Pharmacy
--==_10300106002765.ams1.ulsop.ac.uk==_
Content-Type: text/directory; charset="iso-8859-1"; profile="vcard"; name="david_evans.vcf"
Content-Transfer-Encoding: quoted-printable
Content-Disposition: attachment; filename="david_evans.vcf"
BEGIN:VCARD
VERSION:3.0
PRODID:-//Gordano//NONSGML GMS 9.01.3158//EN
REV:2004-01-22T09:31:19Z
FN:
N:;;;;
EMAIL;TYPE=3DINTERNET:david.evans.ulsop.ac.uk
END:VCARD
--==_10300106002765.ams1.ulsop.ac.uk==_--
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Mar 04 2004 - 10:53:00 PST
