Messages

1

General discussions (Q&As) / Small bug in find_pair

« on: April 13, 2008, 07:41:46 pm »

Hi everyone,

I am quite sure that I found the following bug in find_pair.

When running find_pair on 1VSP.pdb ( http://http://www.rcsb.org/pdb/explore/explore.do?structureId=1VSP ), the base pairs are associate to chain “W” instead of “w”. This capitalization is damaging since 1VSP.pdb also have chain “W” which doesn’t contains base pairs.

Did anyone else encounter this ? have anyone created a workaround ?

Thanks.

2

General discussions (Q&As) /

« on: October 09, 2007, 04:43:00 pm »

Hi Janek,

Thanks for reporting back and giving a detailed account on how you get it done. You have provided the community a very useful working example from a user's prospective. This type of communication is what I value the most. Along the same line, I am hoping 'dnaserver' could follow your lead in getting back to 3DNA forum community.

No problem, however few issues need to be clarified before publishing the final version of the methodology:

To make it more accurate, your 4th step: running "frame_mol" with option "-1" is to re-orient the structure with reference to the first base-pair. This step is necessary because 3DNA rebuilt structure (from "rebuild") is always set w.r.t. to the first base-pair, and it has nothing to do with attaching local helical frames.

4. Run frame_mol to attach local helical frames and set the orientation of a structure.
frame_mol -1 ref_frames.dat 2oaaA.pdb 2oaaA_bp1.pdb

Ok, I have copied the description of the frame_mol from the manual without checking what -1 option does. Thanks, I will update the tutorial

Properly aligning the original (crystal) structure and 3DNA rebuilt one is the essence of this method, and it is exactly the point of the following steps in my original post trying to illustrate.

Code: [Select]

rebuild -atomic bp_step_mut.par pd0001_mut.pdb
cat pd0001_mut.pdb pd0001_bp1.pdb > pd0001_cmb.pdb
rasmol pd0001_cmb.pdbIt would be helpful if you could revise your procedures slightly, and re-post   it, with necessary attachments, in the Users' Contributions section. Putting your recipe there will make it standing out so others can easily find it.

Once again, thanks for contributing back to the 3DNA forum!

Xiang-Jun

I cannot uderstand these steps:

Code: [Select]

cat pd0001_mut.pdb pd0001_bp1.pdb > pd0001_cmb.pdb
rasmol pd0001_cmb.pdb
You generate pd0001_cmb.pdb file which you do not use later for the creation of the final structure. What do you do with rasmol apart from visualizing the structure? Why the procedure worked for me without including these two steps?

Thanks for your help, we will make sth suitable for the Users' Contributions section.

Cheers,
Janek

3

General discussions (Q&As) /

« on: October 07, 2007, 10:21:13 am »

Shorting of or in addition to that, please do share with the community on how the following recipe goes with your projects.


Xiang-Jun

Hi,

It worked for me, however two modifications to the procedure had to be added:
1. Preparation of input DNA file for 3DNA is often needed (changing HETATM to ATOM and new bases name to old ones i.e. DG to G, DA to A etc., changing non-standard numbering of residues)
2. Correcting residue and chain identifiers in the final mutated structure and correcting atom numbering. Because I am preparing this procedure for sort of bioinformatic course I wanted not to use any scripts but rather already available tools. I successfuly used VMD for correcting atom numbering.

Thanks a lot
Janek Kosinski

In case anbybody is interested, here is the exact procedure:

Prepare 2oaa.pdb file for mutagenesis with 3DNA program.
1.   Open 2OAA.pdb in Swiss PDB Viewer
2.   Select chain C (the first DNA chain) and renumber residues starting from 1 (Menu Edit -> Rename Current Layer -> Renumber Selected groups from : 1
3.   Select chain D (the second DNA chain) and renumber residues starting from 1 (Menu Edit -> Rename Current Layer -> Renumber Selected groups from : 1
4.   Select chain A, C and D and Save selected residues in new directory DNA_mutagenesis/ under the name 2oaaA.pdb
5.   Close Swiss PDB Viewer
6.   Open 2oaaA.pdb in WordPad
7.   Press Ctrl+H to open Find & Replace dialog. In “Find what” type HETATM and in “Replace with” type ATOM<space><space> (ATOM and two spaces, these spaces are important)
8.   Press Ctrl+H to open Find & Replace dialog. In “Find what” type  <space><space> D and in “Replace with” type <space><space><space>
9.   Save the file.
Make mutagenesis with 3DNA program.
1.   Open Windows Command Prompt. Go to DNA_mutagenesis/ directory.
2.   Execute:
find_pair 2oaaA.pdb 2oaaA.inp
3.   Run analyze
analyze 2oaaA.inp
4.   Run frame_mol to attach local helical frames and set the orientation of a structure.
frame_mol -1 ref_frames.dat 2oaaA.pdb 2oaaA_bp1.pdb
5.   Copy the files:
cp bp_step.par bp_step_mut.par
cp 2oaaA_bp1.pdb 2oaaA_mutok.pdb
6.   Manually mutate bp_step_mut.par
7.   Open bp_step_mut.par in WordPad. Change 6th base pair (T-A) to C-G
8.   Rebuild PDB file
rebuild -atomic bp_step_mut.par 2oaaA_mut.pdb
9.   Open 2oaaA_mut.pdb and 2oaaA_mutok.pdb in WordPad. Manually replace coordinates of the base (including C1`) in residue 6 of chain C of 2oaaA_mutok.pdb with those from residue 6 of chain A of 2oaaA_mut.pdb and in residue 6 of chain D of  2oaaA_mutok.pdb with those from  residue 6 of chain B of 2oaaA_mut.pdb. After this, in 2oaaA_mutok.pdb, adjust the names of residues and chains in mutated residues and residue numbers. Do not care at wrong atom numbering for mutated bases at the moment.
10.   To correct atom numbering open 2oaaA_mutok.pdb in VMD. After opening save the coordinates under new name 2oaaA_mutatedDNA.pdb.

4

General discussions (Q&As) / dna mutagenesis

« on: October 02, 2007, 04:31:16 pm »

I am doing sort of DNA mutagenesis with 3DNA. Basically what I do is
generating base pair parameters file from PDB file with DNA, than
changing DNA sequence there and the rebuilding DNA in PDB format with
rebuild program. However, rebuild program from 3DNA builds its own phosphate backbone which, albeit similar, is different than the original one.

Could someone tell whether it is possible to run rebuilding with 3DNA in a mode which preserves the original backbone conformation? (one
can write a script to copy phosphate backbone from the original PDB file
but maybe it is unnecessary).

Janek