|
|||||||||||||||||||||||||||||||||
AMBER Archive (2007)Subject: Re: AMBER: heme forcefield atom types
From: L. W. Yang (lwy1_at_pitt.edu)
Dear John, Vlad, Professor Case and Amber experts,
To study the dynamics of deoxy-Myoglobin, I wanted to take advantage of Amber 9_patch34's forcefield, dat/contrib/heme that is dedicated specifically for heme contained molecules. I'm new to Amber and also to heme-contained molecules. However, thanks for John and Vlad's conversation earlier and the suggestions from 0README.txt in the heme folder, I chose the united-atom version, frcmod.hemuni and heme_uni.in for my molecule.
heme_uni.in not only defines the force and connectivity for heme but also the histidine that coordinates the FE(II) as a whole, so I wonder how I should tell the tleap that my histidine 93 actually belongs to part of the HEM (heme) group.
If I simply run tleap on my PDB file as it is (except I delete the 'CONECT' cards in the PDB file), such as the following three-line commands, tleap was easy on me.
>>source leaprc.ff03ua
>>loadamberprep /home/software/amber/amber9_patch34/dat/contrib/heme/heme_uni.in
However, when I did the following 2 lines, the errors and warning messages appeared.
>>mbdo=loadpdb ./1MBD.pdb (I have deleted data following the 'CONECT' cards)
-------- The Error message was too long. I put it at http://www.pitt.edu/~lwy1/1stcheck.txt -------
Can people suggest me how I should modify my PDB file for tleap?
Best regards,
Lee
-----------------------------------------------------------------------
| |||||||||||||||||||||||||||||||||
|