<div dir="auto">By covalent bonds, do you mean that I should add the CONNECT keyword for an atom pair between a HETATM and an regular ATOM in the PDB?</div><div dir="auto">Steve</div><div><br><div class="gmail_quote"><div>On Fri, Feb 15, 2019 at 7:32 PM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">… and now I’m told that currently in ChimeraX, adding pseudobonds wouldn’t help anyway, only covalent bonds.  Sorry, the functionality and documentation were ported from Chimera (where this pseudobond bond stuff did apply), and I was not informed of the difference!<br>
<br>
Apologies for the red herring(s),<br>
Elaine<br>
<br>
<br>
> On Feb 15, 2019, at 4:16 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
> <br>
> I’m told that the pseudobonds may not be needed if:<br>
> <br>
> (A) the two endpoint structures have exactly and only the same atoms (with same names, res numbers, chain IDs)<br>
> -or-<br>
> (B) the “same true” option is used, when some atoms are different but the atoms in question have the same names, res numbers, and chain IDs.  <br>
> <br>
> That may explain why your Mg was included in the first place if it wasn’t pseudobonded to the protein chain.<br>
> <br>
> However, that is orthogonal to the problem of poor movement path, which as far as I know is not solved.  Sorry about that.<br>
> Elaine<br>
> <br>
>> On Feb 15, 2019, at 3:40 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
>> <br>
>> Hi Steve,<br>
>> The problem with the pseudobond file is the blank line.  The reader is not very smart.<br>
>> <br>
>> The pseudobonds (or distances) were just to have the HET atoms simply included in the morph at all, as I understood it and as written here:<br>
>> <<a href="http://rbvi.ucsf.edu/chimerax/docs/user/commands/morph.html#pairing" rel="noreferrer" target="_blank">http://rbvi.ucsf.edu/chimerax/docs/user/commands/morph.html#pairing</a>><br>
>> <br>
>> How the HET atoms move along the trajectory is a separate matter, and it sounds like they are being handled poorly.  Sorry about the limitations.  You could use Help… Report a Bug and attach your session as an example for us if these structures are shareable (not private data).<br>
>> <br>
>> Elaine<br>
>> <br>
>>> On Feb 15, 2019, at 3:20 PM, Steve Chou <<a href="mailto:stevezchou@gmail.com" target="_blank">stevezchou@gmail.com</a>> wrote:<br>
>>> <br>
>>> <br>
>>> Hi Elaine,<br>
>>> (1) I tried to read the pseudobonds from a file (.pb) , <br>
>>> =========================<br>
>>> #1/A:13@NE2  #1/A:704@O<br>
>>> #1/A:13@NE2  #1/A:502@MG<br>
>>> #1/A:13@NE2  #1/A:501@PB<br>
>>> <br>
>>> #2/A:13@NE2  #2/A:704@O<br>
>>> #2/A:13@NE2  #2/A:502@MG<br>
>>> #2/A:13@NE2  #2/A:501@PB<br>
>>> =========================<br>
>>> the program reported an error.<br>
>>> =============================<br>
>>> ValueError: not enough values to unpack (expected 2, got 0)<br>
>>> File "/opt/UCSF/ChimeraX-daily/lib/python3.7/site-packages/chimerax/read_pbonds/readpbonds.py", line 27, in read_pseudobond_file<br>
>>> aspec1, aspec2 = line.split()[:2]<br>
>>> See log for complete Python traceback.<br>
>>> =============================<br>
>>> <br>
>>> (2) Then I tried to add distances on the ChimeraX terminal<br>
>>> =============================<br>
>>> distance #1/A:13@NE2  #1/A:704@O<br>
>>> distance #1/A:13@NE2  #1/A:502@MG<br>
>>> distance #1/A:13@NE2  #1/A:501@PB<br>
>>> <br>
>>> distance #2/A:13@NE2  #2/A:704@O<br>
>>> distance #2/A:13@NE2  #2/A:502@MG<br>
>>> distance #2/A:13@NE2  #2/A:501@PB<br>
>>> =============================<br>
>>> I was able to see the distance labels in the two structures. And I saw an new model (name: distances; ID:3) show up in the "Models" panels. Then I ran morphing again. The Mg and water molecules moved, but still didn't move together with the protein.<br>
>>> Steve<br>
>>> <br>
>>> On Fri, Feb 15, 2019 at 11:31 AM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
>>> Hi Steve,<br>
>>> Yes, if those are the correct specifiers for atoms in your structure.  You'd make a plain text file with name *.pb with lines like that and open it.<br>
>>> <br>
>>> You don’t have to read in pseudobonds from a file, however.  If not too many, probably easier just to add distance monitors.<br>
>>> <<a href="http://rbvi.ucsf.edu/chimerax/docs/user/commands/distance.html" rel="noreferrer" target="_blank">http://rbvi.ucsf.edu/chimerax/docs/user/commands/distance.html</a>><br>
>>> <<a href="http://rbvi.ucsf.edu/chimerax/docs/user/selection.html#context" rel="noreferrer" target="_blank">http://rbvi.ucsf.edu/chimerax/docs/user/selection.html#context</a>><br>
>>> <br>
>>> Elaine<br>
>>> <br>
>> <br>
>> _______________________________________________<br>
>> ChimeraX-users mailing list<br>
>> <a href="mailto:ChimeraX-users@cgl.ucsf.edu" target="_blank">ChimeraX-users@cgl.ucsf.edu</a><br>
>> Manage subscription:<br>
>> <a href="http://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a><br>
>> <br>
> <br>
> <br>
> _______________________________________________<br>
> ChimeraX-users mailing list<br>
> <a href="mailto:ChimeraX-users@cgl.ucsf.edu" target="_blank">ChimeraX-users@cgl.ucsf.edu</a><br>
> Manage subscription:<br>
> <a href="http://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a><br>
> <br>
<br>
</blockquote></div></div><div dir="ltr">-- <br></div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Steve Chou</div><div><br></div><br></div></div></div></div>