<div dir="ltr"><div>Hello again!</div><div>Oh, I got it. So, just two steps are necessary for cleaning heme from all extra H. I try to get an universal protocol for any protein with heme. Indeed I don't know if the Fe ion must be in the heme plane, but one expert of crystallography who works with this protein believes that it must be so. But you are right, even though I delete all extra H, it doesn't help to keep the Fe ion in the heme plane. It's so sad. But thank you very much, Elaine, for your kind help and
<span class="gmail-tlid-translation gmail-translation" lang="en"><span title="" class="gmail-">detailed explanation</span></span>.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">пн, 14 сент. 2020 г. в 23:59, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello,<br>
<br>
HEM[FE] is for specifying the net charge of the Fe ion only.<br>
HEM[non-FE] is for specifying the net charge of the rest of the heme residue (the other atoms except Fe).<br>
<br>
I don't think you needed to do any of the protocol (including defattr) in that previous post, because when I open your structure, delete hydrogens and add them back, commands<br>
<br>
delete H<br>
addh<br>
<br>
....I don't see any extra hydrogens on the ring N. The reason the other person used the defattr file was to fix atom types so that they wouldn't get extra hydrogens.<br>
<br>
I tried minimizing just the heme/Fe from your structure (without the protein to make the test faster) and it still did not go exactly into the plane.<br>
<br>
However:<br>
<br>
(1) are you sure the iron is supposed to be in the plane of the heme? I see some experimental structures where the ion is not in the plane. I'm sure it depends what ligands are binding the metal, etc.<br>
<br>
(2) even if it supposed to be in the plane, you may be expecting too much fine detail from a simple molecular mechanics calculation as done by Chimera. To model small changes from metal oxidation state, you may need to use some kind of quantum-mechanical calculation instead (or in addition to) molecular mechanics, and this cannot be done with Chimera. In Chimera the iron ion is just a ball that does not change size no matter what charge value you put on it, and the charge does not get delocalized no matter what other atoms are around the atom.<br>
<br>
Elaine<br>
<br>
> On Sep 14, 2020, at 11:42 AM, Rio Aquarius <<a href="mailto:rioxel@gmail.com" target="_blank">rioxel@gmail.com</a>> wrote:<br>
> <br>
> On the last step of this protocol I need to specify net charges. Does HEM[FE] mean the charge of the Fe atom or the charge of the heme in state with Fe? But I guess it means the charge of the Fe atom. If it is so can I specify in such a way reduced (Fe+2) and oxidized (Fe+3) states of the heme?<br>
> <br>
> пн, 14 сент. 2020 г. в 21:13, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>>:<br>
> Hello,<br>
> You can put it wherever your computer allows it. In first image, the error is that the file was not found, so you must have not given the same place where you actually put it.<br>
> <br>
> In the second image, it is just a warning message (not an error) that the Render by Attribute dialog will not be shown, but as far as I can tell it did already read the file and assign the attribute correctly. In that case, you can just ignore the warning and continue with your calculations.<br>
> <br>
> I hope this helps,<br>
> Elaine<br>
> -----<br>
> Elaine C. Meng, Ph.D. <br>
> UCSF Chimera(X) team<br>
> Department of Pharmaceutical Chemistry<br>
> University of California, San Francisco<br>
> <br>
> > On Sep 14, 2020, at 11:05 AM, Rio Aquarius <<a href="mailto:rioxel@gmail.com" target="_blank">rioxel@gmail.com</a>> wrote:<br>
> > <br>
> > Thank you for your kind reply. Could you help me a little more? Where must the heme.defattr file be, when I use the 'defattr' command? I have tried to put it in the root directory, but it didn't work and I got such a message:<br>
> > <изображение.png><br>
> > <br>
> > I have also tried to open this file using the tool 'Define Attribute' in Tools/Structure Analysis/, but got such a message:<br>
> > <изображение.png><br>
> > <br>
> > <br>
> > пн, 14 сент. 2020 г. в 19:07, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>>:<br>
> > The heme.defattr file is a link from the URL you already found:<br>
> > <br>
> > > <a href="https://www.cgl.ucsf.edu/pipermail/chimera-users/2019-July/015972.html" rel="noreferrer" target="_blank">https://www.cgl.ucsf.edu/pipermail/chimera-users/2019-July/015972.html</a><br>
> > <br>
> > In the lower part of that page, the link to the defattr file is<br>
> > <br>
> > URL: <<a href="http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20190715/ecb1dce6/attachment.obj" rel="noreferrer" target="_blank">http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20190715/ecb1dce6/attachment.obj</a>><br>
> > <br>
> > ...and you just need to click that and then save it as a plain text file. You can rename that plain text file heme.defattr if you want.<br>
> > <br>
> > I hope this helps,<br>
> > Elaine<br>
> > -----<br>
> > Elaine C. Meng, Ph.D. <br>
> > UCSF Chimera(X) team<br>
> > Department of Pharmaceutical Chemistry<br>
> > University of California, San Francisco<br>
> > <br>
> > > On Sep 12, 2020, at 9:53 AM, Rio Aquarius <<a href="mailto:rioxel@gmail.com" target="_blank">rioxel@gmail.com</a>> wrote:<br>
> > > <br>
> > > Hi,<br>
> > > I have an issue with heme minimization, Fe atom leaves a heme plane. I saw at least one user had the same issue and got help to solve it (<a href="https://www.cgl.ucsf.edu/pipermail/chimera-users/2019-July/015972.html" rel="noreferrer" target="_blank">https://www.cgl.ucsf.edu/pipermail/chimera-users/2019-July/015972.html</a>). I could try the protocol described there, but I don't have the heme.defattr file. I hope you can help me with it too. I have attached a pdb file with this issue.<br>
> > > <CYP51_min300.pdb><br>
<br>
</blockquote></div>