<div dir="ltr">Dear Sir,<div><br></div><div>Thank you for your clarification.</div><div><br></div><div>Just a follow up question: What makes the 'add hydrogen' option protonate an aliphatic amine and skip when it comes to an aromatic proton? Since in both the cases, N is valency wise saturated (trivalent), so what prompts/ guides this differential protonation? Is it pH guided?</div><div><br></div><div>For your kind reference, I am sending image (as attached tiff file) of such a substrate. </div><div><br></div><div>With reference to the image, N1 (in red) is a part of aromatic framework, N2 (in blue) is an aliphatic amine moiety. Both are trivalent (and hence valency wise saturated). Upon activating the 'add hydrogen' option, it is N2 (the aliphatic one) that gets protonated and hence goes to the hypervalent state (as it now becomes tetravalent) but N1 remains unprotonated.</div><div><br></div><div>My question is that upon what basis is this differential protonation done by chimera? Is is guided by their (N1 and N2's) relative pKb (basic strengths) values under a particular pH?</div><div><br></div><div>I will be grateful if you can explain me this issue.</div><div><br></div><div>Thanking you in advance,</div><div><br></div><div>Best Regards,</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="garamond, serif" size="2"><b>Dr. Prabuddha Bhattacharya</b></font></div><div><font face="garamond, serif" size="2"><b>Dept. Of Chemistry</b></font></div><div><b style="font-family:garamond,serif">Adamas University</b></div><div><b style="font-family:garamond,serif">Kolkata-700126</b><br></div><div><b style="font-family:garamond,serif">India</b></div><div><br><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Apr 14, 2020 at 1:34 AM Eric Pettersen <<a href="mailto:pett@cgl.ucsf.edu">pett@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">If you want Chimera to use the protonation state of your structures “as is” and not add any hydrogens of its own, then you need to uncheck the “add hydrogens” check box of Dock Prep (which is also a substep of minimization).  Adding hydrogens will not remove any hydrogens, but will usually add protons to <i>e.g.</i> aliphatic amines, as per your experience.  This is because some PDB structures are only partially protonated, and therefore Chimera does not assume that the presence of pre-existing hydrogens means that adding hydrogens should be skipped.<div><br></div><div><div>--Eric</div><div><br></div><div><span style="white-space:pre-wrap">       </span>Eric Pettersen</div><div><span style="white-space:pre-wrap">   </span>UCSF Computer Graphics Lab</div></div><div><br><div><br><blockquote type="cite"><div>On Apr 13, 2020, at 12:21 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:</div><br><div><div>Dear Dr. Prabuddha Bhattacharya, <br>You didn't say if you used Chimera to add hydrogens or whether the structures already had hydrogens when they were opened.<br><br>If you used Chimera to add hydrogens, it tries to generate states reasonable at neutral pH.  For example, it will not add hydrogens to aspartic acid and glutamic acid side chains of a protein. (If Chimera is adding the -COOH hydrogen, it may be that the atom-type guessing is wrong because the -CO bond lengths are not consistent with a carboxylate group.)<br><br>How AddH in Chimera works is explained here:<br><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/addh/addh.html" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/addh/addh.html</a>><br><br>However, if your molecules already had the hydrogens on them when they were opened (or fetched from PubChem or by SMILES string), Chimera doesn't delete the hydrogens.... it assumes that is what you want.  You should manually delete the hydrogens you don't want.  Or if it's easier, you could try deleting all hydrogens (e.g. command: delete H) and then adding them back with Chimera, but it may not work the way you want if the bond lengths are not consistent with the negatively charged -COO group.<br><br>The force field does not do anything to change protonation states, it just tries to work with what it is given: whatever hydrogens were read in with the structure, or were added in Chimera before charge calculation.<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><blockquote type="cite">On Apr 13, 2020, at 11:51 AM, Prabuddha Bhattacharya <<a href="mailto:b.prabuddha3@gmail.com" target="_blank">b.prabuddha3@gmail.com</a>> wrote:<br><br>Dear Sir,<br><br>I had a query regarding charge assignment to any ligand in chimera. <br><br>While minimizing the energy of a particular ligand (Tools --> Structure Editing --> Minimize Structure) and also during Dock Prep, I found that for most of the aliphatic and alicyclic amines, they are getting protonated with a residual positive charge; while the aromatic amines remain unprotonated (and hence no residual charge). I understand that this is because of their differential pKb values (basic strengths). <br><br>On the other side, deprotonation never occurs with the carboxylic acid (-COOH) group.<br><br>Does the Chimera (force fields) control this issue of protonation and deprotonation considering the physiological pH ?<br><br>I will be extremely grateful if some one may kindly clarify my doubts.<br><br>Thanking you in advance,<br><br>With Best Regards,<br><br>Dr. Prabuddha Bhattacharya <br>Dept. Of Chemistry<br>Adamas University<br>Kolkata 700126<br>India<br><br>_______________________________________________<br>Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br>Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote><br><br>_______________________________________________<br>Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br>Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></div></div></blockquote></div><br></div></div></blockquote></div>