
Hello all! I'm looking to the scripts ( http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/) and can not see a simple example of how to open a PDB. What I need to my extension is to save a PDB relative to a map and then open it to the model panel. Copying it to the model panel does not make sense to me because I need to fix the positions of the moved PDB relative to the map. Could it be something like this? * # Get opened molecule (the one selected in the menu)* * m = self.object_menu.getvalue()* * # Make copy using the copy_molecule native functionality from Chimera* * from Molecule import copy_molecule* * mc = copy_molecule(m)* * # Set copy name* * mc.name <http://mc.name> = m.name.split('.')[0] + '_' + self.chosen_solution + '.pdb'* * copy_path = self.cwd + mc.name <http://mc.name>* * from Midas import open, write* * write(mc, self.map_menu.volume(), self.cwd)* * open(self.cwd)* Thank you in advance!

Hi Pablo, It sounds like what you want to do is open a PDB file, fit it to a map, and then write out the PDB with the fitted coordinates rather than the original coordinates. Is that right? If so, you should look at the Very Basic Chimera Programming Primer <https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/basicPrimer.html> for an example of looping through a series of data files and performing operations on them. The key thing for your purposes is that you want to use the “write” command with the “relative” option in order to write the PDB file relative to the map — so if the map was model 0 and the PDB was model 1, then “write relative 0 1 output_file_name.pdb”. —Eric Eric Pettersen UCSF Computer Graphics Lab
On Mar 12, 2018, at 1:42 AM, Pablo Solar Rodríguez <pablosolar.r@gmail.com> wrote:
Hello all!
I'm looking to the scripts (http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/ <http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/>) and can not see a simple example of how to open a PDB.
What I need to my extension is to save a PDB relative to a map and then open it to the model panel.
Copying it to the model panel does not make sense to me because I need to fix the positions of the moved PDB relative to the map.
Could it be something like this?
# Get opened molecule (the one selected in the menu) m = self.object_menu.getvalue()
# Make copy using the copy_molecule native functionality from Chimera from Molecule import copy_molecule mc = copy_molecule(m)
# Set copy name mc.name <http://mc.name/> = m.name.split('.')[0] + '_' + self.chosen_solution + '.pdb' copy_path = self.cwd + mc.name <http://mc.name/>
from Midas import open, write write(mc, self.map_menu.volume(), self.cwd) open(self.cwd)
Thank you in advance! _______________________________________________ Chimera-dev mailing list Chimera-dev@cgl.ucsf.edu http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev
participants (2)
-
Eric Pettersen
-
Pablo Solar Rodríguez