<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>On Jul 25, 2010, at 7:15 PM, Thomas Ouldridge wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi, I'm afraid I'm not a crystallographer/protein simulator, and so this question might seem quite naive. I run coarse-grained simulations, and I would like to represent the objects as ellipses in chimera. I have the position (which I am able to use) and the orientation that I would like my ellipses to take, encoded as a matrix relative to (x,y,z) cartesian coordinates. This matrix gives me the principal axes and radii of the ellipse. <br><br>As I understand, Chimera is designed to display a "Thermal ellipsoid" - I would like this ellipsoid to have the axes and radii given by my matrix. However, it seems that the input (via pdb format files) is interpreted as an "anisotropic temperature factor", and I can't find a way to convert my matrix into the anisotropic temperature factor (I tried some papers online but they seemed to be talking in terms of transformations between lab frame and crystallographic frames, but I cannot see how that applies - what is my crystallographic frame?).<br><br>In short, is there a way to convert my matrix into the "anisotropic temperature factor" representation so that Chimera can read it and convert it back? I did try just plugging in the matrix elements to the PDB file but got some weird ellipsoids out.</div></blockquote><br></div><div>Hi Tom,</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>The ANISOU record uses the same coordinate frame as the ATOM records, so you should be able to just plug in your matrix values.  There are a couple of tricks though.  First, the ordering of the six numbers in the ANISOU record is (1,1), (2,2), (3,3), (1,2), (1,3), (2,3), <i>i.e.</i> first the main diagonal:</div><div><br></div><div>(1,1), (2,2), (3,3)</div><div><br></div><div>then filling out the first row:</div><div><br></div><div>(1,2), (1,3)</div><div><br></div><div>then finally the third element of the second row:</div><div><br></div><div>(2,3).</div><div><br></div><div>Also, the values are integers, and are divided by a factor of 10000 to get the final values.</div><div><span class="Apple-tab-span" style="white-space:pre">     </span>Is this what you tried?</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>It should also be possible to directly assign to the atom's anisoU attribute using a Python script if the above doesn't work for some reason.  Let me know if you need that.</div><div><br></div><div>--Eric</div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                       <span class="Apple-converted-space"> </span></span>Eric Pettersen</font></p><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                       <span class="Apple-converted-space"> </span></span>UCSF Computer Graphics Lab</font></p><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                        </span><a href="http://www.cgl.ucsf.edu">http://www.cgl.ucsf.edu</a></font></p><br class="Apple-interchange-newline"></div></span> </div><br></body></html>