<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"><div dir="ltr"><div dir="ltr"><div>Hi Elaine,<br></div>(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>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">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><span style="font-style:italic">See log for complete Python traceback.</span><br></p>=============================<br></div><div><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></div><div>I was able to see the distance labels in the two structure. 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></div><div>Steve<br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 15, 2019 at 11:31 AM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@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">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>
> On Feb 14, 2019, at 9:49 PM, Steve Chou <<a href="mailto:stevezchou@gmail.com" target="_blank">stevezchou@gmail.com</a>> wrote:<br>
> <br>
> Hi Elaine,<br>
> Thanks for your detailed explanation! <br>
> What's the format of the pseudobond file? Should it look like this?<br>
> ===============================<br>
> #1/A:13@NE2 #1/A:502@MG<br>
> ===============================<br>
> All the best,<br>
> Steve<br>
> <br>
> On Thu, Feb 14, 2019 at 7:43 PM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
> Hi Steve,<br>
> Actually I’m surprised the Mg and water were included (even if present at both ends) because I thought they had to be tethered by a bond or pseudobond (like a distance measurement) to the biopolymer chain, as described here under “Atoms in Common”:<br>
> <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>
> For the case where composition changes, you could add in a fake Mg ion and water to the structure that’s missing them, but you’d have to decide where to put them. Also it was my understanding that they had to be tethered to the biopolymer chain (e.g. add distance measurement, which can then be hidden), but maybe that requirement has changed, or the Mg does already have a metal-coordination pseudobond to the biopolymer.<br>
> <br>
> Or if you don’t care if the Mg or water move during the morph, you could just display those atoms from the original model (that was input to morphing) even though they are not present in the trajectory.<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>
> > On Feb 14, 2019, at 12:51 PM, Steve Chou <<a href="mailto:stevezchou@gmail.com" target="_blank">stevezchou@gmail.com</a>> wrote:<br>
> > <br>
> > Thanks, Elaine! <br>
> > You command "rainbow structures" works beautifully. <br>
> > There is "a Mg ion" and "a water molecule" in the two conformations.<br>
> > (1) If the two conformations have exactly the same compositions (by shifting the structure), the "Mg ion" and the "water molecule" do not follow the same path as the protein. They move but separate from the protein.<br>
> > (2) If the two conformations have slightly different compositions, the "Mg ion" and the "water molecule" do not show up in the morphed pdb ensemble.<br>
> > Is there a way to morph the specific components properly?<br>
> > Steve<br>
> > <br>
> <br>
> <br>
> <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" rel="noreferrer" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br>
<br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Steve Chou</div><div><br></div><br></div></div></div></div>