This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Protein 3D Structure Principal Axes

Can anyone suggest some free software for calculating the principal axes for a protein 3d structure? Ideally the input will be the PDB file of 3d coordinates, and the program will return the moment of inertia tensor and the principal axes. I DON'T want the molecule to orient further along the principal axes.

There are a few examples of codes available in the net, but it would be good if someone knows about a well known software to do this job.

Thanks in advance

pdb protein
• 9,661 views • link
updated by Biostar Biostar • written by Woa Woa

3 answers

Here are several options:

  • Paul Grayson wrote a Tcl/Tk-based orient script you can use in VMD. It will calculate the principal axes. For you, just skip the 'move' command. http://www.ks.uiuc.edu/Research/vmd/script_library/scripts/orient/

  • I rewrote the above code in Java. It needs the Colt package for some of the linear algebra functions. http://mem.bio.pitt.edu/?q=node/89

  • The MODELLER software package has a built-in orient command that is very fast: http://salilab.org/modeller/9v4/manual/node180.html For you, if you don't want to actually orient it, just skip the mdl.write() line and use the r.translation() data. Note that since MODELLER is written for homology stuff this will ignore hydrogens by default. To keep your hydrogens you'll want to add this line right after the env = environ() line:

    env.io.hydrogen = True

• 0 views • link
written by Keith Callenberg Keith Callenberg

great!!! many thanks

• 0 views • link
written by Woa Woa

Dear Chris,

Albeit not being a "well known" software, I have a Python script that nicely does the job. In addition, it produces a .pml script for a rendering with the Pymol molecular viewer. I am not (yet) used to github so I posted the script here http://pastebin.com/m8wryusp

I hope it helps.

Cheers

Pierre

• 0 views • link
written by Pierre Poulain • Pierre Poulain

Thanks for the code

• 0 views • link
written by Woa Woa

Can DeepView (Swiss Pdb-Viewer) do this job?

• 0 views • link
written by Dejian Dejian

Log in to answer this question.

More posts like this