I was thinking develop just the protein in AmberTools, use amb2gmx to
convert then solvate and ionize with the gromacs tools pointed at the amber
ports of the relevant water and ion model.
The reason I don't use the amber ports for the proteins is that they have
novel residues.
You don't happen to know another defined tool chain for going from quantum
calculations to gromacs topology do you?
John.
On Mon, Aug 18, 2008 at 5:12 PM, Alan <alanwilter.gmail.com> wrote:
> Good you saw it's hardcoded. Next step you realise it's not difficult to
> tweak your top file to match tip4pew water.
> Alan
>
>
> On Mon, Aug 18, 2008 at 6:04 AM, Grange Hermitage <ghermitage.gmail.com>wrote:
>
>> The definition of the water model is hardcoded in amb2gmx.pl
>>
>>
>>
>> On Mon, Aug 18, 2008 at 11:58 AM, Grange Hermitage <ghermitage.gmail.com>wrote:
>>
>>> Hi all,
>>>
>>> when the following leap script (AmberTools 1.2) is executed and followed
>>> by amb2gmx.pl, the .top file is suitable for TIP3P and the .gro is suitable
>>> for TIP4PEW. The .top file knows nothing about EPW. Gromacs then remarks
>>> that the nunber of atoms in the .top and .gro don't match. Is there
>>> something amiss in the leap file? There are no waters in chp16.pdb .
>>>
>>> Thanks,
>>> John.
>>>
>>>
>>> source leaprc.ff99SB
>>> WAT = T4E
>>> set WAT.1 name "T4E"
>>> loadamberparams frcmod.tip4pew
>>> loadoff ions94.lib
>>> loadamberparams frcmod.ionsjc_tip4pew
>>> loadoff chk.lib
>>> loadoff chd.lib
>>> chp16 = loadpdb chp16.pdb
>>> check chp16
>>> solvatebox chp16 TIP4PEWBOX 5 iso
>>> addions chp16 Na+ 17
>>> addions chp16 Cl- 13
>>> check chp16
>>> saveamberparm chp16 chp16leap.prm chp16leap.crd
>>> quit
>>>
>>>
>>
>
>
> --
> Alan Wilter S. da Silva, D.Sc. - CCPN Research Associate
> Department of Biochemistry, University of Cambridge.
> 80 Tennis Court Road, Cambridge CB2 1GA, UK.
> >>http://www.bio.cam.ac.uk/~awd28<<
>
-----------------------------------------------------------------------
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 Aug 20 2008 - 06:07:18 PDT
