- Chimera-users - RBVI Mailing Lists

POV-RAY tracing in Mac OSX (a workaround)
by Miguel Ortiz Lombardia 20 May '09

20 May '09

Hi, I'm almost sure that this problem has been reported here... but couldn't find anything about it. When I try to ray-trace [Chimera 1.3 build 2572 (Darwin/aqua)] I consistently get this error: File Context (5 lines): } object { text { ttf "verdana.ttf" "$" 0.01, 0 texture Parse Error: Cannot open font file. It appears even if I have no labels on the image. I managed to solve the problem but only when I asked for the povray intermediate files not to be deleted and then manually edit them: A) image.ini file: I changed the line: Library_Path="/Applications/Chimera.app/Contents/Resources/share/fonts/" to Library_Path="/Applications/Chimera.app/Contents/Resources/share/fonts/ freefont/" B) image.pov file: All calls to verdana.ttf were changed to FreeSans.ttf It would be great if you point me to an easier solution or if you could make such changes in the code. Thanks! Miguel -- Miguel Ortiz Lombardía Architecture et Fonction des Macromolécules Biologiques UMR6098 ( CNRS, U. de Provence, U. de la Méditerranée ) Case 932 163 Avenue de Luminy 13288 Marseille cedex 9 France Tel : +33(0) 491 82 55 93 Fax: +33(0) 491 26 67 20 e-mail: miguel.ortiz-lombardia(a)afmb.univ-mrs.fr Web: http://www.pangea.org/mol/spip.php?rubrique2

2 1

residue order for match/rmsd
by Klaus Kopec 19 May '09

19 May '09

Hi everybody, at first I would like to thank you for Chimera, it's the best structure viewer I am aware of! Now to my question: I would like to match residues in two structures, but not in the numerical order. My expectation was that if I run "match #0:1,2,3@CA #1:3,2,1@CA" I would end up with least squares fitted coordinate pairs (#0:1, #1:3), (#0:2, #1:2) and (#0:3, #1:1). What I get, however, is a least-squares fit of pairs (#0:1, #1:1), (#0:2, #1:2) and (#0:1, #1:3). Is it possible to use the order of coordinate pairs given in the command line (e.g. pair first specified coordinate of protein 1 with first specified coordinate of protein 2, second with second, etc.) instead of the numerical order of the residue numbers? I would be really glad if someone could help me out on this! Klaus PS: I use Chimera 1.3.2577 under a 64-Bit Kubuntu. I also tried some older (1.2.something) version which gave the same results.

2 2

2d labels and pov-ray option
by Bala subramanian 19 May '09

19 May '09

Friends, I opened a molecule, made some 2d labels and saved the image in jpeg format. But the 2d labels are not apperaing in the final picture. What could be the reason for the same ? When i try to set POV-RAY as True and save the image, the POV-RAY progress till 98% and stays for longtime but dosent save the image. Has anyone come across such problems before ? Thanks, Bala

2 1

residue order for match/rmsd
by Klaus Kopec 19 May '09

19 May '09

Hi everybody, at first I would like to thank you for Chimera, it's the best structure viewer I am aware of! Now to my question: I would like to match residues in two structures, but not in the numerical order. My expectation was that if I run "match #0:1,2,3@CA #1:3,2,1@CA" I would end up with least squares fitted coordinate pairs (#0:1, #1:3), (#0:2, #1:2) and (#0:3, #1:1). What I get, however, is a least-squares fit of pairs (#0:1, #1:1), (#0:2, #1:2) and (#0:1, #1:3). Is it possible to use the order of coordinate pairs given in the command line (e.g. pair first specified coordinate of protein 1 with first specified coordinate of protein 2, second with second, etc.) instead of the numerical order of the residue numbers? I would be really glad if someone could help me out on this! Klaus PS: I use Chimera 1.3.2577 under a 64-Bit Kubuntu. I also tried some older (1.2.something) version which gave the same results.

1 0

Principal axis of inertia
by A. Di Fenza 18 May '09

18 May '09

Hi Tom, I saw your message on the web and I tried to use your script on Chimera 1.3. Unfortunately it does not work, it give me the following message: 0 models opened Do you think it is because of the Chimera version? Do you have suggestion? Thank you very much, Armida Oops. That script I posted for calculating the inertia of a molecule had some problems. The printed matrix of axes was transposed though the correct axes were used for the displayed ellipsoid. Also I was using inertia matrix I_ij = sum_over_k(m_k*r_ki*r_kj) omitting the diagonal term. Should have been I_ij = sum_over_k(m_k*(|r_k| - r_ki*r_kj)). And finally I incorrectly computed the ellipsoid size along each axis. Everything was a little bit wrong. Here's the fixed version. Tested in Chimera version 1.2500 but should work with earlier versions. Tom

5 6

multiple morphs
by Charles Midgett 18 May '09

18 May '09

Dear Chimera developers, I am trying to make a movie of how we think domains of our complex interact. Essentially, there will be three morphs, made from different pdb files, that need to be displayed simultaneously. How can I get the three trajectories to happen at the same time? Sincerely, Charles Midgett

2 1

Saving Render by Attribute Values
by Kevin Wu 18 May '09

18 May '09

Hi all, I'm having trouble saving some data to the session. I've rendered a figure using the "render/select by attribute" by either the occupancy or beta field in a pdb file. However, when I save session and restore it, the colors that I've chosen seem to be erased, and need to be reset again. Is there a way to save this data? Thanks in advance! Kevin Wu

4 4

fit map --nogui
by sean connell 18 May '09

18 May '09

Hello, I am having some problems using/adapting a script I found by searching the mailing list. The script is pasted below but to summarize it should print out the transformation matrix used to fit one map to another. In principle I would like to do the same but in addition to writing out the matrix I would also like to apply it and save the resulting map. I am currently using chimera 1.4 build 27496 (Aqua on OSX 10.5.7) the first problem I have is just invoking the python interpreter..... if I use: Chimera.app/Contents/MacOS/chimera --nogui fittomap.py Executing fittomap.py... Traceback (most recent call last): File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ __main__.py", line 65, in <module> value = chimeraInit.init(sys.argv) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimeraInit.py", line 407, in init chimera.openModels.open(a, prefixableType=1) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimera/__init__.py", line 1357, in open models = func(filename, *args, **kw) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimera/__init__.py", line 794, in _openPython loadFunc(sandboxName, fileName, f) File "fittomap.py", line 70, in <module> fit_map_in_map(p1, p2) File "fittomap.py", line 24, in fit_map_in_map from VolumeViewer import Data_Region, full_region, Rendering_Options ImportError: cannot import name Data_Region Error while processing fittomap.py: ImportError: cannot import name Data_Region (see reply log for Python traceback info) # ----------------------------------------------------------------------------- # Example script for fitting one map in another without the graphical user # interface. # # chimera --nogui fitnogui.py # # It can also be run using the graphical Chimera interface using File/ Open. # def fit_map_in_map(map1_path, map2_path, initial_map1_transform = None, map1_threshold = None, ijk_step_size_min = 0.01, # Grid index units ijk_step_size_max = 0.5, # Grid index units max_steps = 100, optimize_translation = True, optimize_rotation = True): # Files have to have file suffix indicating volume format. from VolumeData import open_file data1 = open_file(map1_path) data2 = open_file(map2_path) from VolumeViewer import Data_Region, full_region, Rendering_Options ro = Rendering_Options() map1 = Data_Region(data1, full_region(data1.size), '', ro) map2 = Data_Region(data2, full_region(data2.size), '', ro) if initial_map1_transform: from PDBmatrices.matrices import chimera_xform xf = chimera_xform(initial_map1_transform) map1.surface_model().openState.globalXform(xf) use_threshold = (map1_threshold != None) from FitMap import map_points_and_weights, motion_to_maximum points, point_weights = map_points_and_weights(map1, use_threshold) if len(points) == 0: if use_threshold: print 'No grid points above map threshold.' else: print 'Map has no non-zero values.' return move_tf, stats = motion_to_maximum(points, point_weights, map2, max_steps, ijk_step_size_min, ijk_step_size_max, optimize_translation, optimize_rotation) if initial_map1_transform: from PDBmatrices.matrices import multiply_matrices move_tf = multiply_matrices(move_tf, initial_map1_transform) header = ('\nFit map %s in map %s using %d points\n' % (map1.full_name(), map2.full_name(), stats['points']) + ' correlation = %.4g, overlap = %.4g\n' % (stats['correlation'], stats['overlap']) + ' steps = %d, shift = %.3g, angle = %.3g degrees\n' % (stats['steps'], stats['shift'], stats['angle'])) print header from FitMap.gui import transformation_description tfs = transformation_description(move_tf) print tfs # ----------------------------------------------------------------------------- # Example call fitting map into itself. # p1 = p2 = '/usr/local/src/chimera-demos/volume/examples/groel.mrc' fit_map_in_map(p1, p2)

2 1

Solvent Occupancy Calculation
by Saikumar Ramadagu 18 May '09

18 May '09

Hi Elaine, I know people in the past have asked a question about the solvent Occupancy Calculation, but please help me. I have an oligosaccharide simulaiton trajectory from AMBER9. I used a 100ps trajectory for solvent occupancy calculation using Chimera. Everything looks absolutely fine. The questions I have are as follows: I have 0-220 in my levels section of volume viewer. If I fix the level of 80, what does that mean? When I'm sliding the bar of levels from left-right what information about solvent density am I obtaining? For example does a value of 80 mean low density when compared with a value of 70? or is it the other way? When I have two different oligosaccharide simulations and I want to say that oligosaccharide 1 has more water density around it, what value of levels should I choose? What I mean by this statement is do I select same value of levels? or do I calculate the percentage i.e., 80/220 for first model gives me 36% and for second sugar if the levels are 0-300 the i should choose a level of 109 which gives 36%? Also is there any reference paper by which Chimera calculates the solvent occupancy? Thanks very much for your time Sai

2 1

Formation of coordination bond
by Enrico Morelli 18 May '09

18 May '09

Dear all, We would like to know if is it possible, when we have two molecules docked together, to form a coordination bond between a metal ion present in one of the two molecules and a donor atom of the other molecule, in such way forming a covalent bond between the two proteins. Thanks -- ------------------------------------------------------------------- (o_ (o_ //\ Coltivate Linux che tanto Windows si pianta da solo. (/)_ V_/_ +------------------------------------------------------------------+ | ENRICO MORELLI | email: morelli(a)CERM.UNIFI.IT | | * * * * | phone: +39 055 4574269 | | University of Florence | fax : +39 055 4574253 | | CERM - via Sacconi, 6 - 50019 Sesto Fiorentino (FI) - ITALY | +------------------------------------------------------------------+

3 2