User:Wayne Decatur/I-Ppo Morph Methods: Difference between revisions

From Proteopedia
Jump to navigationJump to search
Wayne Decatur (talk | contribs)
No edit summary
Wayne Decatur (talk | contribs)
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 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 -  
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:


Line 18: Line 19:


for model in structure:<br>
for model in structure:<br>
     for chain in model:<br>
<nowiki>
     for chain in model:</nowiki><br>
         i = 2493<br>
         i = 2493<br>
         for residue in chain:<br>
         for residue in chain:<br>