
Dear Chimera people, I would like to save a lot of pdb. I have for example 500 pdb. Each saved pdb file should contains 2 pdbs of these 500 pdbs. I can write a script for this but the command write does not work for more than two pdb. Is there any other way? thank you Julien

Dear Julien, It is difficult to answer without knowing more about your situation, but I will describe what can be done in Chimera. The "write" command creates one PDB file containing atoms from one model ID number. If you opened your 500 structures from 500 separate PDB files, they will have 500 different model ID numbers (#0 #1 #2 ...), so "write" cannot be used to combine them. However, if you opened the 500 structures from a single PDB file, where MODEL/ENDMDL records are used to indicate the start and stop of each structure, they will be opened as submodels but with the same model ID number (#0.1 #0.2 #0.3 ...). For examples of this single- file input, see any NMR structure ensemble from the Protein DataBank, such as 1plx. In this single-file case, you can select any set of the submodels and write their coordinates to a single PDB file with "write" and the "selected" keyword: open 1plx select #0.1 #0.3 #0.5 write selected 0 1-3-5.pdb select #0.2 #0.25 write selected 0 2-25.pdb ... Disadvantages: (1) you still need to do a separate selection and write command for each pair; depending on how many pairs you planned to write, you might want to write a script to create the Chimera script. (2) the written files will still contain the MODEL and ENDMDL lines for submodels for which no atomic coordinates are included, which you might want to remove depending on what you are going to do with the files. The effect of not removing them is that when you open the file later in Chimera, it will still say there are 80 submodels (or whatever the total number you started with was) in the Model Panel even though some contain 0 atoms. The display will only show the coordinates you wanted, however. Given these issues, you might prefer to use your own scripts. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng@cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Dec 4, 2007, at 1:24 PM, Julien Orts wrote:
Dear Chimera people, I would like to save a lot of pdb. I have for example 500 pdb. Each saved pdb file should contains 2 pdbs of these 500 pdbs. I can write a script for this but the command write does not work for more than two pdb. Is there any other way? thank you Julien
participants (2)
-
Elaine Meng
-
Julien Orts