User:Wayne Decatur/I-Ppo Morph Methods: Difference between revisions
From Proteopedia
Jump to navigationJump to search
No edit summary |
No edit summary |
||
| Line 5: | Line 5: | ||
I loaded 1a73 structure into swiss-pdb viewer and let it magic fit the apo enzyme in 1evx. I disabled movement on those two structures and loaded the unbound DNA. When I had turned it so it was close, I enable movement on all and turned about 90 degress and then disabled movement again all all but the unbound DNA and moved it into plane with the other structures. I save each layer individually. I edited the aligned 1evx file to match atom numbers with the aligned 1a73 and removed the SO4s. I deleted the 2 MGs in the aligned 1a73 file. | I loaded 1a73 structure into swiss-pdb viewer and let it magic fit the apo enzyme in 1evx. I disabled movement on those two structures and loaded the unbound DNA. When I had turned it so it was close, I enable movement on all and turned about 90 degress and then disabled movement again all all but the unbound DNA and moved it into plane with the other structures. I save each layer individually. I edited the aligned 1evx file to match atom numbers with the aligned 1a73 and removed the SO4s. I deleted the 2 MGs in the aligned 1a73 file. | ||
Next, I needed to renumber the atoms in the unbound DNA to make it easy to compare to file 1a73 for matchin up number of atoms since morphing needs equal numbers. I also needed to change the chain ids from a and b to match c and d in 1a73. (At first, I thought Biopython would do this but after using the [http://biopython.org/DIST/docs/tutorial/Tutorial.html tutorial] to install numpy and Biopython, I found at [http://maple.rsvs.ulaval.ca/mediawiki/index.php/Manipulating_PDB_files_using_BioPython the help site for manipulating pdb files with Biopython] and also the [http://biopython.org/DIST/docs/tutorial/Tutorial.html cookbook]that it doesn't concern itself with atom numbers that I could see and instead just parses all the other information like chain id and residue number to make that easy to work with. In fact I wrote a program with it to easily renumber residues and save file and learned it wasn't what I want - | Next, I needed to renumber the atoms in the unbound DNA to make it easy to compare to file 1a73 for matchin up number of atoms since morphing needs equal numbers. I also needed to change the chain ids from a and b to match c and d in 1a73. (At first, I thought Biopython would do this but after using the [http://biopython.org/DIST/docs/tutorial/Tutorial.html tutorial] to install numpy and Biopython in addition to my already present installation of Python 2.5, I found at [http://maple.rsvs.ulaval.ca/mediawiki/index.php/Manipulating_PDB_files_using_BioPython the help site for manipulating pdb files with Biopython] and also the [http://biopython.org/DIST/docs/tutorial/Tutorial.html cookbook]that it doesn't concern itself with atom numbers that I could see and instead just parses all the other information like chain id and residue number to make that easy to work with. In fact I wrote a program with it to easily renumber residues and save file and learned it wasn't what I want - | ||
the program: | the program: | ||
< | <code> | ||
from Bio.PDB.PDBParser import PDBParser | from Bio.PDB.PDBParser import PDBParser | ||
| Line 29: | Line 29: | ||
w.save('test-r1.pdb') | w.save('test-r1.pdb') | ||
</ | </code> | ||
So I searched for something else thinking alhtough I could break it up and do it with a simple Python program it would be nice to find something to do it.) | So I searched for something else thinking alhtough I could break it up and do it with a simple Python program it would be nice to find something to do it.) | ||
I found [http://dicsoft2.physics.iisc.ernet.in/pdbgoodies/ PDB Goodies] to take care of the chain renumbering. To get started, just a bit down the page I clicked the not overly obvious spot that said 'click here to see available options in PDB Goodies' to start and on the next page uploaded my pdb file. (Note on this page a link to <help> in the bottom right corner that does summarize what PDB Goodies does). And I chose 'Change Identifiers' and changed A and B to C and D easily using the form generated from my uploaded file and getting the results. (Note PDB goodies also renumbers residues numbers, generating a nice form using the file.) | I found [http://dicsoft2.physics.iisc.ernet.in/pdbgoodies/ PDB Goodies] to take care of the chain renumbering. To get started, just a bit down the page I clicked the not overly obvious spot that said 'click here to see available options in PDB Goodies' to start and on the next page uploaded my pdb file. (Note on this page a link to <help> in the bottom right corner that does summarize what PDB Goodies does). And I chose 'Change Identifiers' and changed A and B to C and D easily using the form generated from my uploaded file and getting the results. (Note PDB goodies also renumbers residues numbers, generating a nice form using the file.) | ||
