AMBER Archive (2005)

Subject: Re: AMBER: read pdb back to xleap after solvation of chloroform

From: David A. Case (case_at_scripps.edu)
Date: Thu May 19 2005 - 20:38:49 CDT


On Fri, May 20, 2005, chaiann ng wrote:
>
> I am trying to read the pdb back to leap and output prmtop and inpcrd after
> the dynamic of solvation of chloroform box.
>
> I read in leaprc.gaff; leaprc.ff03; loadOff by lib, load frcmod.chcl3, when
                                      ^^^^^^^^^^^^^^

I' not sure what you mean by the above...

> i loadpdb in and try to output inpcrd and prmtop, this is what i get.
>
> Creating new UNIT for residue: CL3 sequence: 403
> Created a new atom named: H1 within residue: .R<CL3 403>
> Created a new atom named: C1 within residue: .R<CL3 403>
> Created a new atom named: CL1 within residue: .R<CL3 403>
> Created a new atom named: CL2 within residue: .R<CL3 403>
> Created a new atom named: CL3 within residue: .R<CL3 403>
>
> If I also load chcl3.in then i get this.

You definitely have to do this, so that leap knows about the CL3 unit.

> +Currently only Sp3-Sp3/Sp3-Sp2/Sp2-Sp2 are supported
> +---Tried to superimpose torsions for: *-H1-C1-*
> +--- With Sp0 - Sp0
> +--- Sp0 probably means a new atom type is involved
> +--- which needs to be added via addAtomTypes

Check your output to make sure leap is not trying to add some x-H1-C1-x
torsion; be sure to put TER cards in your pdb file between each separate
molecule. It could be that leap is trying to make a covalent link between the
chloroform and some previous residue.

Why are you trying to do this, anyway?

...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu