November 2015 - Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] RMSD script
by Elaine Meng 18 Nov '15

18 Nov '15

Yes, as described in more detail in the “RMSD Analysis” description you can select specific atoms (for example command “select protein & @n,ca,c,o” to select the peptide backbone atoms) and use the “Restrict map to current selection, if any” option. That option is on by default. If it was on and you didn’t select any atoms, it used all atoms. > See MD Movie manpage, especially sections “Startup and Input” and “RMSD Analysis" > <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.…> More about selection: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/selection.html> Best, Elaine On Nov 18, 2015, at 2:47 PM, Ruth Nelson <Ruth.Nelson.180(a)nd.edu> wrote: > Hi Elaine, thank you for your response. > > Using MD Movie worked to calculated the RMSDs, but I am unsure if the output values are full atom or for the backbone (Ca, C, N, O) only. Is there some way I can specify this? > > Thank you very much, > Ruthie > > On Wed, Nov 18, 2015 at 5:06 PM, Elaine Meng <meng(a)cgl.ucsf.edu> wrote: > Hi Ruthie, > If all of the PDB files have exactly the same atoms with the same names, and no END lines, you could try opening all of them as one trajectory in MD Movie. That tool has an all-by-all RMSD calculation option and you can save the results to a text file. You would still have to process that file into CSV, however, as there is no CSV output option, but this would greatly decress the number of steps you have to script yourself, and it would be done outside of Chimera in whatever type of scripting you prefer. > > I.e. you would start MD Movie ( in menu under Tools… MD/Ensemble Analysis) and then choose to open the trajectory as PDB, multiple files. They have be named similarly with numbers but the numbers don’t have to be consecutive. If loading the trajectory is successful, the MD Movie playback dialog will appear, and from that dialog you could choose Analysis…RMSD Map. The RMSD can be calculated from all atoms or only specify ones that you have selected. Then there is a button to “Save RMSDs” on the resulting map window. > > See MD Movie manpage, especially sections “Startup and Input” and “RMSD Analysis" > <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.…> > > I hope this helps, > Elaine > ----- > Elaine C. Meng, Ph.D. > UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab > Department of Pharmaceutical Chemistry > University of California, San Francisco > > On Nov 18, 2015, at 9:38 AM, Ruth Nelson <Ruth.Nelson.180(a)nd.edu> wrote: > > > Hello, > > > > I have a matrix of 10x10 .pdb files, and I need to calculate the RMSD from each PDB file to every other PDB file in that matrix. When I do this in the command line of Chimera, I type something like "match #0 #1" or "match #0 #2" for each of my PDB files, but I would like to automate this process so that it does not take as long. I have hundreds of files and they are all are named similarly (file_1.pdb, for example) but they are not numbered consecutively (file_1.pdb, file_17.pdb, file_43.pdb, for example). If possible, I would like the RMSD values to be output into a CSV file so that I can open it in excel. > > > > I have very little experience writing scripts and no idea where to begin to write this one, so I would greatly appreciate any help. I have been reading some of the responses given to other people regarding similar scripts, but I am having a hard time translating that to what I would like to do. > > > > Thank you very much, > > Ruthie > > > > > > > -- > Ruth Nelson, Graduate Student > Department of Chemistry and Biochemistry > University of Notre Dame > rnelson9(a)nd.edu

1 0

RMSD script
by Ruth Nelson 18 Nov '15

18 Nov '15

Hello, I have a matrix of 10x10 .pdb files, and I need to calculate the RMSD from each PDB file to every other PDB file in that matrix. When I do this in the command line of Chimera, I type something like "match #0 #1" or "match #0 #2" for each of my PDB files, but I would like to automate this process so that it does not take as long. I have hundreds of files and they are all are named similarly (file_1.pdb, for example) but they are not numbered consecutively (file_1.pdb, file_17.pdb, file_43.pdb, for example). If possible, I would like the RMSD values to be output into a CSV file so that I can open it in excel. I have very little experience writing scripts and no idea where to begin to write this one, so I would greatly appreciate any help. I have been reading some of the responses given to other people regarding similar scripts, but I am having a hard time translating that to what I would like to do. Thank you very much, Ruthie

2 1

Changing default framerate of Animation tool
by Eiros Zamora, Juan 17 Nov '15

17 Nov '15

Dear Chimera experts, I would like to know if there’s a way to increase the default 25 fps framerate when using the ‘Animation' tool to make movies. I tried increasing it from the ‘Movie Recorder’ menu, but it’s still producing the video at 25 fps. Also, I am not sure about what does the ‘Nth frame for discrete transitions’ property of items in the timeline mean. Should I increase this to get more FPS? I would like to do this because I am trying to do a video of an MD trajectory showing at the beginning the water box (and then hiding it), but at the moment the videos that I am getting feel very ‘clunky’ Cheers Juan

2 1

how to prepare ligand by command line?
by Albert 17 Nov '15

17 Nov '15

Hello: I've got a ligand in 2D .sdf format without any hydrogens. I am just wondering how can we prepare this compound for docking? Is it possible to run it by command from Chimera? Thanks a lot Albert

3 2

Problem using Chimera
by BRUXELLES Tatiana (ETU VIE) 17 Nov '15

17 Nov '15

Dear Sir, dear Madam I am a Master student at the Université de Strasbourg in France and I am currently doing an internship in a bioinformatics lab : the LBGI in Strasbourg. One of my tasks during this internship is to make Chimera read about 80 files with the extensions .scf. The problem is that software doesn't to be able to read my command file. I can't even make it read a print instruction. To execute those files I use the command : read File_name in command line on chimera. I am gonna join the files so you can see. Could you suggest a solution to my problem ? Kind regards, Bruxelles Tatiana

2 1

helices
by A G Szabo 17 Nov '15

17 Nov '15

I sent a message yesterday regarding how I can change helices from ribbons to smooth wires. I am pleased that I can indicate that I figured it out.. Arthur G Szabo Professor Emeritus Wilfrid Laurier University

1 0

smooth backbones
by A G Szabo 17 Nov '15

17 Nov '15

I have no problems with ribbons. However I would like to show some helices as backbones only. How can I achieve the result shown on the attachment for this helix. Don't know where to start looking. Any " Actions" only gave me this. Thank you Arthur G. Szabo Professor Emeritus Wilfrid Laurier University

1 0

isosurface and blobs
by Raffaele Raucci 13 Nov '15

13 Nov '15

Hi all, I imported in chimera a map file with .cube extension. I was able to show isosurface after setting specific isovalue. Can you tell me if it is possible, by command line, count the number of different blobs displayed as isosurface? Is it also possible determine the area and volume of each blob? Thanks in advance, Raphael

3 4

Volume transformation based on Phenix symmetry operators
by Morris, Kyle 13 Nov '15

13 Nov '15

Dear Chimera developers and community, I have a volume and an associated symmetry operator identified by another software package for this map. I would like to transform the map according to the symmetry operator to inspect what the symmetry related map looks like, does anyone know how to do this based on the appended matrix? Thanks in advance for your help! Best wishes, Kyle ================================== >From Phenix: new_operator rota_matrix -0.0877 -0.5114 0.8549 rota_matrix 0.9251 -0.3601 -0.1206 rota_matrix 0.3695 0.7803 0.5047 tran_orth 578.6845 434.1862 -508.9168 center_orth 918.6495 790.3525 1034.0200

2 1

ESP
by Putkey, John 13 Nov '15

13 Nov '15

I have a question about electrostatic surface potential coloring. It works fine when only one .pqr and .dx file is open. I can uncolor, recolor and adjust the ESP color scale, but odd things happen when muliple .pqr and .dx files are open. Strange color patterns are obtained even though the correct MSMS surface and .dx file are selected in the Surface Color dialog box. Its as if the program is reading two .dx files. Is it me or a bug with the system? Cheers John John A. Putkey Professor Biochemistry and Molecular Biology University of Texas Medical School 6431 Fannin Houston Texas, 77030 John.Putkey(a)uth.tmc.edu<mailto:John.Putkey@uth.tmc.edu> (713) 500-6061 Office (713) 500-6062 Laboratory (713) 500-5073 NMR Laboratory

2 1