<div>Hi Eric</div><div>Thanks for your suggestions. I tried the mentioned steps but still facing the same issue (see attached file). I forgot to mention that I am minimizing the modelled structure, maybe that's why it's creating this problem. In modelling I modelled the structure with heme and also with out heme (and later on pasted the heme coordinates from template to model structure). I minimized both structures and ended up with the same problem. I am not sure where am I making mistake. I'd be grateful if you could guide me more.</div><div>Regards</div><div>Sara</div><br><br><blockquote class="ReferenceQuote" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0" name="replyContent">-----Original Messages-----<br>
<b>From:</b><span id="rc_from">"Eric Pettersen" <pett@cgl.ucsf.edu></span><br>
<b>Sent Time:</b><span id="rc_senttime">2019-07-16 06:55:59 (Tuesday)</span><br>
<b>To:</b> sara20@mail.ustc.edu.cn<br>
<b>Cc:</b> "chimera-users@cgl.ucsf.edu BB" <chimera-users@cgl.ucsf.edu><br>
<b>Subject:</b> Re: [Chimera-users] energy minimization<br><br>Hi Sara,<div><span class="Apple-tab-span" style="white-space:pre">       </span>There are a couple of problems here.  One is that the heme coordinates / bond lengths aren’t very good, so Chimera’s estimates for the atom types is wrong for some of the atoms.  The other problem is that because the iron is displaced out of the porphyrin ring, addh is adding hydrogens to the inside of the ring that would normally be prevented by the presence of the iron.  Those hydrogens repel the iron of course.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>So, to correct things, you need to:</div><div><br></div><div>1) Delete all hydrogens</div><div>2) Correct the atom types</div><div>3) Add the hydrogens back</div><div>4) Delete inward-facing porphyrin-ring hydrogens</div><div>5) Add charges</div><div>6) Minimize</div><div><br></div><div>The commands (Favorites→Command Line) corresponding to those steps are:</div><div><br></div><div>1) “del H”</div><div>2) “defattr /path/to/heme.defattr raise false”  (I have attached the heme.defattr file) </div><div>3) “addh”</div><div>4) “del :hem@hd”</div><div>5) “addch”  (this will take awhile with AM1-BCC charges)</div><div>6) Minimize however you were doing it before.</div><div><br></div><div>I have also attached a session with all the above steps performed.  The iron is not yet planar with the ring, but additional minimization will undoubtedly get it there.</div><div><br></div><div>—Eric</div><div><br></div><div><div><span class="Apple-tab-span" style="white-space: pre;">    </span>Eric Pettersen</div><div><span class="Apple-tab-span" style="white-space: pre;">     </span>UCSF Computer Graphics Lab</div><div><br></div></div><div></div></blockquote>