<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 Catherine,<div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Remind me how exactly you are creating the axis in Chimera, and send me whatever script you are using to compute it yourself…</div><div class=""><br class=""></div><div class="">—Eric</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Jul 18, 2017, at 8:24 AM, Catherine Jenifer Rajam Rajendran <<a href="mailto:catrajen@umail.iu.edu" class="">catrajen@umail.iu.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hi Eric,<div class=""><br class=""></div><div class="">I am able to generate the axis. Like you said, we used the x,y,z axis subtracting it from centroid position and did the SVD.<br class=""><div class="">For most cases, the axis is same as Chimera but for few it is little different. I checked the number of atoms used between both and it is same. <br class=""></div><div class="">Are you using all the atoms or only specific atoms or is there any other thing I am missing?<br class=""></div></div><div class=""><br class=""><div class="">Thanks,</div><div class="">Catherine</div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Fri, Jun 30, 2017 at 7:13 PM, Eric Pettersen <span dir="ltr" class=""><<a target="_blank" href="mailto:pett@cgl.ucsf.edu" class="">pett@cgl.ucsf.edu</a>></span> wrote:<br class=""><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><div style="word-wrap:break-word" class="">Chimera simply calls numpy.linalg.svd (singular value decomposition) on an array of the xyz coordinate (after the centroid position has been subtracted from them). The documentation for numpy.linalg.svd is here:<div class=""><br class=""></div><div class=""><a target="_blank" href="https://docs.scipy.org/doc/numpy/reference/generated/numpy.linalg.svd.html" class="">https://docs.scipy.org/doc/num<wbr class="">py/reference/generated/numpy.<wbr class="">linalg.svd.html</a></div><div class=""><br class=""></div><div class="">It looks like it in turn calls LAPACK routine_gesdd, which is documented, in general, here:</div><div class=""><br class=""></div><div class=""><a target="_blank" href="https://software.intel.com/en-us/mkl-developer-reference-c-singular-value-decomposition-lapack-computational-routines?language=es" class="">https://software.intel.com/en-<wbr class="">us/mkl-developer-reference-c-s<wbr class="">ingular-value-decomposition-la<wbr class="">pack-computational-routines?la<wbr class="">nguage=es</a></div><div class=""><br class=""></div><div class="">It’s an implementation of singular value decomposition, described in general terms here:</div><div class=""><br class=""></div><div class=""><a target="_blank" href="https://en.wikipedia.org/wiki/Singular_value_decomposition#Relation_to_eigenvalue_decomposition" class="">https://en.wikipedia.org/wiki/<wbr class="">Singular_value_decomposition#R<wbr class="">elation_to_eigenvalue_decompos<wbr class="">ition</a></div><div class=""><br class=""></div><div class="">The math is pretty complicated.</div><div class=""><br class=""></div><div class="">—Eric</div><div class=""><br class=""><div class="">
<div style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; word-wrap: break-word;" class=""><div class=""><span style="white-space:pre-wrap" class="gmail-m_4101912461970981994m_2254719073479732994Apple-tab-span"> </span>Eric Pettersen</div><div class=""><span style="white-space:pre-wrap" class="gmail-m_4101912461970981994m_2254719073479732994Apple-tab-span"> </span>UCSF Computer Graphics Lab</div><div class=""><br class=""></div></div><br class="gmail-m_4101912461970981994m_2254719073479732994Apple-interchange-newline"><br class="gmail-m_4101912461970981994m_2254719073479732994Apple-interchange-newline">
</div>
<br class=""><div class=""><blockquote type="cite" class=""><div class=""><div class="gmail-m_4101912461970981994h5"><div class="">On Jun 30, 2017, at 10:09 AM, Catherine Jenifer Rajam Rajendran <<a target="_blank" href="mailto:catrajen@umail.iu.edu" class="">catrajen@umail.iu.edu</a>> wrote:</div><br class="gmail-m_4101912461970981994m_2254719073479732994Apple-interchange-newline"></div></div><div class=""><div class=""><div class="gmail-m_4101912461970981994h5"><div dir="ltr" class="">Hi Elaine,<div class=""><br class=""></div><div class="">Thank you for helping me understand it. Can you give me more details about how PDB file is decomposed to find the plane and normal. Like, how the eigenvectors are calculated from the atomic coordinates.</div><div class=""><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Catherine</div></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Thu, Jun 29, 2017 at 7:23 PM, Elaine Meng <span dir="ltr" class=""><<a target="_blank" href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>></span> wrote:<br class=""><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><br class="">
<br class="">
> On Jun 29, 2017, at 2:14 PM, Catherine Jenifer Rajam Rajendran <<a target="_blank" href="mailto:catrajen@umail.iu.edu" class="">catrajen@umail.iu.edu</a>> wrote:<br class="">
><br class="">
> Hi,<br class="">
><br class="">
Hi Catherine,<br class="">
<span class=""><br class="">
> I am trying to understand how the plane is defined in chimera. Can you please help me with that?<br class="">
><br class="">
> "Define a plane for the specified atoms. Eigenvectors/values are calculated from the atomic coordinates after subtracting the position of their non-mass-weighted centroid. The plane is anchored at the centroid and aligned with the first two eigenvectors (the third eigenvector is normal to the plane)."<br class="">
><br class="">
> Define a plane for the specified atoms.<br class="">
</span>> - specified atoms, i use 'sel protein’<br class="">
<br class="">
OK, so if using the graphical tool (Axes/Planes/Centroids) it uses the selected atoms, meaning all the protein atoms if you had used “sel protein.” If using the “define” command, you can just specify “protein” directly in the command to use all the protein atoms, without making a selection.<br class="">
<span class="">><br class="">
> Eigenvectors/values are calculated from the atomic coordinates<br class="">
> - atomic coordinates, so x,y,z coordinates of all atoms. (Can you please explain a bit more like how and what you read from PDB file to define plane).<br class="">
<br class="">
</span>For the atoms you specified, there are ATOM lines in the PDB file with X,Y,Z coordinates. Each atom is a point in 3D (X,Y,Z).<br class="">
<span class=""><br class="">
><br class="">
> after subtracting the position of their non-mass-weighted centroid<br class="">
> - how is non-mass-weighted centroid is calculated<br class="">
<br class="">
</span>It is the centroid of the points (atomic positions), not sure how else to say it. The average position of all the atomic points is the centroid. Non-mass-weighted just means all the atomic points are treated equally, instead of heavier atoms being treated as more important than lighter atoms.<br class="">
<span class="">><br class="">
> The plane is anchored at the centroid<br class="">
> - which centroid is this? Is the protein moved to the origin(0,0,0) position?<br class="">
><br class="">
</span>The centroid is the same centroid as in your previous question. Nothing is moved.<br class="">
><br class="">
> Thanks,<br class="">
> Catherine<br class="">
<br class="">
You’re welcome. I hope this helps,<br class="">
Elaine<br 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="">
<br class="">
<br class="">
</blockquote></div><br class=""></div></div></div>
______________________________<wbr class="">_________________<br class="">Chimera-users mailing list: <a target="_blank" href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a target="_blank" href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mail<wbr class="">man/listinfo/chimera-users</a><br class=""></div></blockquote></div><br class=""></div></div></blockquote></div><br class=""></div></div>
</div></blockquote></div><br class=""></div></body></html>