<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Elaine, <div class=""><br class=""></div><div class="">Thank you for getting back with your detailed response. The 6574 is part of the PO4 residue because that is how Maestro ‘<i class="">preserves</i>’ the information that the atom has originally come from a phopsphate? I am using Maestro to do the modelling, essentially breaking a X—H bond and making a Y—H bond with the same hydrogen. </div><div class=""><br class=""></div><div class="">Essentially we are trying to make ‘movies’ of 3D enzyme mechanisms. This means that morph interpolate should be able to know that atom X from residue N goes to position Y of residue M. Am I missing something obvious? </div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Souro </div><div class=""><br class=""></div><div class="">—</div><div class=""><br class=""></div><div class="">Souro Chowdhury</div><div class="">Researcher, Thornton Group </div><div class="">European Bioinformatics Institute, EMBL-EBI </div><div class="">Cambridge, United Kingdom </div><div class=""><a href="http://www.ebi.ac.uk/research/thornton" class="">http://www.ebi.ac.uk/research/thornton</a></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 11 Jul 2017, at 19:09, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Souro,<div class="">Your two input files I2.pdb and I3.pdb have severe errors. Looking at I2.pdb with a text editor, I see that the hydrogen that is supposedly on the histidine is actually part of the PO4 residue according to that file, and the connectivity of the histidine-1082 sidechain is badly messed up. For example, hydrogen HD1 with serial number 3865 is bonded to two other atoms (according to the CONECT section at the bottom of the file) . See image attached with atoms labeled by serial number and “ligand” selected, so that you can see that hydrogen 6574 is part of the ligand residue PO4 even though the CONECT lines have it bonded to the histidine.</div><div class=""><br class=""></div><div class="">When I tried the morphing, it seemed to me that the movement of O3 was reasonable, however (just goes between starting and ending positions).</div><div class=""><br class=""></div><div class="">Initially I was going to answer that Chimera's morphing was developed mainly for the biopolymer parts (protein and nucleic chains) and I wouldn’t necessarily expect good performance on ligand conformations. However, when I saw how badly messed up the input files are, now I think it might work a lot better if those connectivity issues (and putting the H in the histidine residue instead of the PO4 residue) were fixed.</div><div class=""><br class=""></div><div class="">I don’t use Pymol, so I didn’t look at that script. If morphing works in pymol it may be because it is ignoring the connectivity information (CONECT lines at the bottom) in your PDB files.</div><div class=""><br class=""></div><div class="">I hope this clarifies the situation. Best,</div><div class="">Elaine</div><div class=""><div class="">----------<br class="">Elaine C. Meng, Ph.D. <br class="">UCSF Chimera(X) team<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><span id="cid:AD4C8031-B609-4616-96B6-D351D0FEF628@compbio.ucsf.edu"><i2-mess.png></span></div><br class=""><blockquote type="cite" class="">On Jul 11, 2017, at 6:25 AM, Souroprobho Chowdhury <<a href="mailto:souro@ebi.ac.uk" class="">souro@ebi.ac.uk</a>> wrote:<br class=""><br class="">Hello, <br class=""><br class="">I am trying to use morph interpolate on two (for this post) PDBs (I2.pdb and I3.pdb, attached) but something peculiar is happening. The Hydrogen on the HIS 1082 (N1) does a ‘cork-screw’ like movement. Essentially, instead of interpolating directly from position A to B, it is taking a more round-about angle. A similar movement is observed for the O3 in the T1 O ligand. <br class=""><br class="">This problem has come up recurrently in my work and I am yet to find a solution. I have used most of the options that apply on morph interpolate that are mentioned on the chimera website without much improvement. <br class=""><br class="">I am using Maestro by Schrodinger to model my structures. <br class=""><br class="">Please find attached I2.pdb , I3.pdb and mov_all.pdb. mov_all.pdb is the trajectory file, works on Pymol (I am attaching a pymol script called auto-script.pml as well for selection of the correct residues, and a reasonable representation). The chimera interpolation script is also attached. <br class=""><br class="">Thanks! <br class=""><br class="">Souro <br class=""><br class=""><span id="cid:8E8B3B62-4FDF-4DD9-9753-0B0F21659EDF@compbio.ucsf.edu" class=""><I2.pdb></span><br class=""><span id="cid:F963E47E-1F67-41C6-AD20-5E45B34C93D5@compbio.ucsf.edu" class=""><I3.pdb></span><br class=""><span id="cid:26182318-0C30-4CCC-8011-BEB2C5777831@compbio.ucsf.edu" class=""><<a href="http://interpolate.com/" class="">interpolate.com</a>></span><br class=""><span id="cid:1E1DC7CC-C7A8-4A49-B0A6-A631D55BB1DD@compbio.ucsf.edu" class=""><mov_all.pdb></span><br class=""><span id="cid:0413DA44-B96F-4046-9722-F18567171286@compbio.ucsf.edu" class=""><auto-script.pml></span><br class="">_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></blockquote><br class=""></div></div></div></blockquote></div><br class=""></div></body></html>