<div dir="ltr">It did help!<div><br></div><div>Thanks a lot.</div><div><br></div><div>Best Regards,</div><div>Subha</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 18, 2016 at 5:12 PM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Subha,<br>
Chimera makes its best guess as to the atom types in a nonstandard residue, using the bond lengths and angles. In this case it did not figure out the atom type you expected.  You would need to correct the atom type of that N atom to what you want before adding hydrogens and to do so, you would need to know its atom name (or serial number, or some way of specifying it in the command).  There is no general way to specify the net charge of a nonstandard residue in the addh command.<br>
<br>
The atom types in Chimera are described here:<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/idatm.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/<wbr>chimera/docs/UsersGuide/idatm.<wbr>html</a>><br>
<br>
…you can see the atom-type assignments as labels, for example if the molecule was model 3:<br>
<br>
labelopt info idatmType<br>
label #3<br>
<br>
…and  can use the setattr command to change them, for example something like:<br>
<br>
setattr a idatmType N2  #3@N8<br>
label #3<br>
<br>
I understand you may not know ahead of time about which atoms will have this problem, but this is the only way I can think of to directly control the protonation state of atoms in nonstandard residues with a command.<br>
<br>
I hope this helps,<br>
Elaine<br>
----------<br>
Elaine C. Meng, Ph.D.<br>
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<div><div class="h5"><br>
<br>
> On Aug 18, 2016, at 3:22 PM, Subha Kalyaanamoorthy <<a href="mailto:kalyaana@ualberta.ca">kalyaana@ualberta.ca</a>> wrote:<br>
><br>
> Hi There,<br>
><br>
> Can anybody please suggest me how to use 'setattr' (or any other function) to assign a protonation state for an atom in a non-standard residue.<br>
><br>
> I have a ligand and I am using a python script to add hydrogens and charge to the ligand molecule using chimera no-gui. I see that chimera adds an extra hydrogen to a nitrogen atom in the ligand and assigns a charge of +1 to the molecule. Removing this hydrogen using GUI is simple but not practical for large set of molecules. So, it would be helpful if somebody suggest me a command-line to set specific protonation states for atoms (in non-standard residue).<br>
><br>
> Thanks and Best Regards,<br>
> Subha<br>
<br>
<br>
</div></div>______________________________<wbr>_________________<br>
Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>
Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" rel="noreferrer" target="_blank">http://plato.cgl.ucsf.edu/<wbr>mailman/listinfo/chimera-users</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Dr. Subha Kalyaanamoorthy<div>Post Doctoral Fellow</div><div>Faculty of Pharmacy and Pharmaceutical Sciences</div><div>University of Alberta</div><div>Edmonton, Canada.</div><div><br></div></div></div>
</div>