- ChimeraX-users - RBVI Mailing Lists

Fwd: Unable to load --offscreen and --nogui options
by Sai Janani Ganesan 15 Sep '20

15 Sep '20

Hello All, I am unable to save images from the command line, I get an OpenGL error. I get openGL errors on both Mac and Linux. Any suggestions? ChimeraX --offscreen --nogui test.cxc Traceback (most recent call last): File "/opt/chimerax-1.0-1.fc32/lib/python3.7/runpy.py", line 193, in _run_module_as_main "__main__", mod_spec) File "/opt/chimerax-1.0-1.fc32/lib/python3.7/runpy.py", line 85, in _run_code exec(code, run_globals) File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/ChimeraX_main.py", line 895, in <module> exit_code = init(sys.argv) File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/ChimeraX_main.py", line 453, in init sess = session.Session(app_name, debug=opts.debug, silent=opts.silent) File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/chimerax/core/session.py", line 441, in __init__ self.logger = logger.Logger(self) File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/chimerax/core/session.py", line 501, in __setattr__ if not name.startswith('_') and self.snapshot_methods(value, base_type=StateManager) is not None: File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/chimerax/core/session.py", line 535, in snapshot_methods from chimerax import graphics as gr File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/chimerax/graphics/__init__.py", line 34, in <module> from .opengl import Texture, Lighting, Material File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/chimerax/graphics/opengl.py", line 61, in <module> from OpenGL import GL File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/PyOpenGL-3.1.5-py3.7.egg/OpenGL/GL/__init__.py", line 4, in <module> from OpenGL.GL.VERSION.GL_1_1 import * File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/PyOpenGL-3.1.5-py3.7.egg/OpenGL/GL/VERSION/GL_1_1.py", line 14, in <module> from OpenGL.raw.GL.VERSION.GL_1_1 import * File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/PyOpenGL-3.1.5-py3.7.egg/OpenGL/raw/GL/VERSION/GL_1_1.py", line 7, in <module> from OpenGL.raw.GL import _errors File "/opt/chimerax-1.0-1.fc32/lib/python3.7/site-packages/PyOpenGL-3.1.5-py3.7.egg/OpenGL/raw/GL/_errors.py", line 4, in <module> _error_checker = _ErrorChecker( _p, _p.GL.glGetError ) AttributeError: 'NoneType' object has no attribute 'glGetError'-- Thanks, Sai

3 3

Loading OpenMM platform plugins
by Thomas Mulvaney 15 Sep '20

15 Sep '20

Hi, I've been having fun developing a plugin for ChimeraX and want to firstly thank everyone involved for such a great platform to build upon. I'm starting to introduce a molecular dynamics component and was excited to see that ChimeraX ships with OpenMM support. Unfortunately, the OpenMM Platform plugins (which add CPU, CUDA, OpenCL etc) functionality, do not seem to be loaded by default at least in ChimeraX 1.1 on macOS. I couldn't see any documentation about how to load these plugins, but for anyone running into this issue, the following code got things working for me: ------ import chimerax openmm_plugins_dir = os.path.join(chimerax.app_lib_dir, 'plugins') Platform.loadPluginsFromDirectory(openmm_plugins_dir) ------ Kind regards, Tom

2 2

headless chimerax is built in via --offscreen and --nogui
by Moth, Christopher W 14 Sep '20

14 Sep '20

For other longtime and unimaginative 🙂 legacy chimera users who might be googling/searching for 'headless chimerax', hopefully this note will be picked up by search engines The details you seek: http://www.rbvi.ucsf.edu/pipermail/chimerax-users/2019-April/000472.html No need to reply.

1 0

Choosing solvent exposed residues automatically
by Mutum Yaikhomba 09 Sep '20

09 Sep '20

Hi I am working with a very big pdb model: pdb - 6g2j. I was wondering if there is a quick way to select all residues (surface residues) that are exposed to the solvent (including those facing the membrane)? This is instead of going over each residue one by one and selecting them. Additionally, say you have selected these solvent-exposed residues and now you want to filter them based on Hydrophobicity - is there a way to filter based on the standard Kyte-Doolittle Hydrophobicity index for these residues? If there is, it would be of great help and very much appreciated. Thanks Yaikhomba

2 3

mousemode and ISOLDE
by Alexis Rohou 06 Sep '20

06 Sep '20

Hi all, Not sure whether this is a ChimeraX question, or an ISOLDE question, or a bit of both. I am using ChimeraX 1.0 (2020-06-04) and ISOLDE version ChimeraX_ISOLDE-1.0.1-cp37-cp37m-macosx_10_13_x86_64.whl. When I start ChimeraX, I get the following left mouse bindings, which I believe are default: - ctrl-left click = select - shift-ctrl-left click = select add When I just start ISOLDE via the GUI menu (no further action), this seems to clear the mousemode. Or, at least, ctrl-click does not select anymore, and shift-ctrl-click also doesn't seem to do anything. To remedy this, I issue the following commands: mousemode control left select mousemode control shift left 'select add' This restores ctrl-left to 'select' as expected, but shift-ctrl-left also triggers 'select', not 'select add'. Couple of questions: (1) Is ISOLDE intentionally overriding mousemode when it is launched? I'm 99.9% this didn't use to be the case pre-1.0... (2) How do I restore shift-ctrl-left to 'select add'? I assume I am getting the syntax wrong, but I did read the doc and thought I was following it, so I guess I need some help understanding the syntax. Cheers, Alexis

3 3

Fit Map Question/Articles
by Luis José Castillo Valverde 03 Sep '20

03 Sep '20

Hello, I am currently working on my thesis work (this is why they have received many questions from me these days), this is why I am collecting information on biomolecular structures, and I wanted to know if you know of any article that comments on the pros and cons of its "Fit Map in Map" algorithm? This because I have been searching on this topic and have not been able to find anything. Sorry for the inconvenience and thanks in advance. Best regards, Luis -- Luis José Castillo Valverde Ingeniero en Computación, TEC Costa Rica Estudiante de maestría (MSc) académica en Computación, TEC Costa Rica Telegram: @LCastillo98 LinkedIn: luisjosecastillo98

3 2

Re: [chimerax-users] [Chimera-users] User-defined graphics in chimera: NMR tensor plotting
by Tom Goddard 31 Aug '20

31 Aug '20

Hi Albert, Not sure if you are talking about using a modified version that spherical Python code in Chimera or using it in ChimeraX. I'll assume ChimeraX. http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics <http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics> It is much easier to shift and rotate the vertices than your approach. Python numpy arrays make this very convenient. After the line of code that makes the numpy array of vertices va = array(vertices, float32) you could shift it to a new origin with just va += (32.5, 13, -24.1) Rotation is also easy. from chimerax.geometry import rotation r = rotation(axis = (0,1,0), angle = 45) va = r * va na = r * na Here I rotate the vertices and the normal vectors about 0,0,0. Then I could apply the translation after that. You might want to do a rotation given by a 3 by 3 rotation matrix and the translation in one step, like this from chimerax.geometry import Place t = Place(matrix = ((0,1,0, 32.5), (-1,0,0, 13), (0,0,1, -24.1))) t.transform_points(va, in_place = True) t.transform_vectors(na, in_place = True) Here I used a 3 by 4 matrix where the first 3 columns were the rotation and the 4th column is the translation. I'll update the example code to show this. These geometric transforms are described in the ChimeraX Python programming documentation: https://www.cgl.ucsf.edu/chimerax/docs/devel/bundles/geometry/src/geometry.… As for what went wrong with you np.array([x0,y0,z0]).T approach -- that transpose operation put the vertices in a non-standard memory order -- you created a 3 by N array and then did a transpose to get an N by 3 array but numpy is smart and does not allocate a new array and instead just adjusts the strides and uses the original array. My guess is the calculate_surface_normals() code has a bug in that it does not handle non-standard orders for the vertex array. So the normals were completely wrong and that effects the lighting. In the ChimeraX Python shell (menu Tools / General / Shell) the variable "session" is the ChimeraX session. Tom > On Aug 31, 2020, at 2:43 PM, Albert Smith-Penzel <Albert.Smith-Penzel(a)medizin.uni-leipzig.de> wrote: > > Dear all, > > Thanks for such quick replies and helpful answers. I had two follow-up questions on the approach that Tom suggested in ChimeraX. > > The first is maybe a glitch, or more likely something I don’t fully understand in python/chimera. > > In the code Tom links to (http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics <http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics>), a list of vertices (x,y,z tuples) is created. Suppose I want to shift the center of the spherical harmonic and/or rotate it. One way to do this is to extract x,y,z from the vertices list and modify them, then put them back in the list. My first step is to switch the tuple to a list in the code (switch xyz = (r*sp*ct, r*sp*st, r*cp) to xyz = [r*sp*ct, r*sp*st, r*cp] ). Then I extract the elements, do some editing, and then put them back into a list of vertices. > > x0,y0,z0=np.array(vertices).T > #Some editing of x0,y0,z0 - although editing is not required to produce the strange effects. > vertices=np.array([[x,y,z] for x,y,z in zip(x0,y0,z0)]) > #vertices=np.array([x0,y0,z0]).T > > Finally, I need to reconstruct the vertices list. The first approach (currently uncommented) works just fine. The second approach (commented) doesn’t create any errors, but the lighting is very strange, in that there are alternating dark and light areas on the surface of the harmonic. I would have claimed that the two approaches should produce practically identical numpy arrays, but clearly that’s not the case. > > Anyway, it’s not particularly important for me, since I can work around it, but it does take some time to figure out, because one would think the latter approach should also work. > > The second question is: is it possible to get the “session” object into an interactive python session? Usually, I work interactively with python, and have a number of functions that create various scripts and input files for chimera, and then launch chimera from a shell, calling some script on startup. While this works, it would be really nice if one could manipulate an existing chimera session (maybe a bit off of the original topic now). > > In any case, the new graphics with the tensors come out very nicely- thanks so much for your help. > > Cheers, > Albert > > > >> On 28. Aug 2020, at 20:58, Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> wrote: >> >> WARNUNG: Diese E-Mail kam von außerhalb der Organisation. Klicken Sie nicht auf Links oder öffnen Sie keine Anhänge, es sei denn, Sie kennen den Absender und wissen, dass der Inhalt sicher ist. >> >> Hi Albert, >> >> You need finer phi and theta steps to get a smoother appearance. I found that 10 degree steps gives chunky appearance, while 3 degree steps givs smooth appearance. Chimera does not have code to make an arbitrary surface smoother. >> >> Maybe you are interested in making these surfaces all in Python instead of using BILD. Here is how it is done in our newer ChimeraX program. >> >> http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics <http://plato.cgl.ucsf.edu/trac/chimera/wiki/SphericalHarmonics> >> >> Tom >> >> <harmonics.png> >> >>> On Aug 28, 2020, at 1:54 AM, Albert Smith-Penzel <Albert.Smith-Penzel(a)medizin.uni-leipzig.de <mailto:Albert.Smith-Penzel@medizin.uni-leipzig.de>> wrote: >>> >>> Hello, >>> >>> I have a question about improving user defined graphics in Chimera. >>> >>> For background (maybe not necessary for answering): >>> I’m working with nuclear magnetic resonance and MD simulation, and so I want to plot dynamically averaged spin-spin (H–C) dipole tensors, obtained from MD, on top of the corresponding H–C bonds. I’ve been relatively successful, but while one has a very high-resolution picture of the molecule itself, the tensors themselves are visibly triangulated. >>> >>> How I get this plot: >>> Each tensor is expressed first in spherical coordinate (theta, phi, r). In python, I run Delauney from scipy.spatial on theta and phi to get a list of vertices that can be used for triangulation of the tensor. Then, for each tensor, I convert theta, phi, and r into cartesian coordinates, shift the center so that the tensor sits on top of the correct bond. Finally, I create a .bild file where I write out each triangle using the coordinates and vertices. Tensors have negative and positive values (red and blue in the plot), and so I only take triangles for which r for all vertices is positive, and then in a second step, take only triangles where all vertices are negative. This all gets written to the .bild file, and then opened in chimera on top of the molecule (see picture). >>> >>> The question: >>> The result is certainly acceptable, however, there is clearly a quality difference in the molecule rendering vs. the tensor rendering (obviously, because the tensors are just made of triangles that I’ve created in python). Is there a better way to bring the surface data describing the tensors into chimera that would then allow usage of some sort of smoothing within chimera for a nicer image? (And possibly with a relatively simple data format….) >>> >>> Thanks for your help! >>> >>> Cheers, >>> Albert Smith >>> >>> <example_triangulation.png> >>> _______________________________________________ >>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>> Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users <https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users> >> > > _______________________________________________ > Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu > Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users

1 0

Fit in Map - Overlap
by Luis José Castillo Valverde 28 Aug '20

28 Aug '20

Hi, As I understand it, the overlay resulting from the "Fit in Map" command corresponds to the sum over the fit map grid points of the product of the fit map value and the reference map value at that point, determined by trilinear interpolation. I want to calculate the percentage of overlap of two maps, it would be as simple as calculating the formula "min ((sum of the grid points of map 1), (sum of the grid points of map 2)) / (overlap obtained from the fit ) "the problem I have that I do not know how to obtain those values (" sum of the grid points of map 1 "," sum of the grid points of map 2 "); I tried to do "Fit in Map" of map 1 with map 1 and use the overlap that returns as the sum of the points of map 1 but it didn't work. So I wanted to know if you can tell me how to get the sum of the grid points of a given map. Regards, Luis -- Luis José Castillo Valverde Ingeniero en Computación, TEC Costa Rica Estudiante de maestría (MSc) académica en Computación, TEC Costa Rica Telegram: @LCastillo98 LinkedIn: luisjosecastillo98

2 1

MD trajectories in ChimeraX
by David Slochower 24 Aug '20

24 Aug '20

Hi, are there any plans to implement a "metafile" interface (e.g., http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#…) for loading trajectories with ChimeraX? In particular, it would be great to be able to load `prmtop` and `nc` trajectories, even if it were available only through the command line and without a GUI. (This can be considered a follow up to this message: http://www.rbvi.ucsf.edu/pipermail/chimerax-users/2019-December/000773.html) I'm aware of the Molecular Dynamics Viewer in the Toolshed, but that appears to be geared for multi-PDB trajectories. Thanks.

2 3

Stack sidebar tools as tabs
by Sean Connell 22 Aug '20

22 Aug '20

Hello, I was wondering if there is a command to automatically stack the tools (i.e. log, volume viewer) in the sidebar as TABs? For example to put in the preferences so it starts this way or to make an alias so it can be done quickly. Thanks in advance for any tips you have. Sean

2 1