
Dear Elaine, Thank you very much for your detailed reply. I worked on the PDB id 1A6M as suggested by you and followed the same instructions. I used +2 for FE and -4 for HEME and found following message at the end. The details about the calculations are also attached in this mail along with mol2 file. Before calculation, I removed the water, oxy and so4 in het entry and did the calculations with protein and heme alone. Here is what I am getting at the end using BC-AM1 Charges for residue HEM[non-FE] determined Total charge for #0: 2.0001 Unable to find GAFF type for #0:154.het@FE what is this message? log file along with the charge snapshot is attached with this mail. So, Can I use this mol2 file in DOCK 6 calculations? I am lil worried since the pre-calculated charges for amber are different from the the one I am getting using chimera and antechamber combination. I will highly appreciate your reply and suggestions. Thanks in advance, s On Nov 6, 2007 12:16 PM, Elaine Meng <meng@cgl.ucsf.edu> wrote:
Hi "S" We are glad you like Chimera.
(1) HEM and FE charges. AMBER charges are assigned for the protein, but not for the nonstandard residue HEM/FE. Since FE is only one atom, it will put whatever integer charge you specify on that atom (without trying to "calculate" the charge). You can see the assigned charge by labeling the atom with its charge value (see below) or looking in the charge column of the output Mol2 file, if you wrote it out. The Add Charge dialog will always show +2 for your FE because it always uses the same method to guess what the charge should be, but that does not mean the currently assigned charge is +2. I just tested this by specifying FE charge +7 and then labeling the atom; Add Charge did assign +7. For the example structure PDB 1a6m, however, +2 is the correct charge for FE and -4 the correct charge for the rest of HEM, protoporphyrin IX.
To label atoms with their charges, select the atoms of interest, then either use the commands:
labelopt info charge label sel
or use the menu Actions... Label... other... and choose the attribute "charge" from the pulldown menu in the labeling dialog.
(2) Charge set. Only all-atom charges can be assigned by Chimera, not united-atom (where the nonpolar CH CH2 and CH3 are treated as single points).
I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng@cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html
On Nov 6, 2007, at 7:58 AM, snoze pa wrote:
Hi, I am new in this mailing list but found chimera more useful than any other software. I was building a protein that contain HEME. When I was calculating the charge using chimera interface, then it is asking for HEM[FE] and HEM[nonFE] formal charge. I know chimera is using MMTK and AMBER forcefield. If I will give the default value then it is assigning +2 charge to Fe, which is not a match based on AMBER forcefield. Meanwhile I try to change the Formal charge manually, but still it gives the same +2 for Fe. Anybody used chimera to prepare a protein having heme for dock or geometry build.
My another question is: Is it possible in chimera to merge nonpolar hydrogen charges on the heavy atoms, the technique used in autodock.
I will highly appreciate your help and suggestions in this matter. thanks in advance. s