HELANAL — analysis of protein helices

Author:Lily Wang
Copyright:GNU Public License v3

New in version 1.0.0.

This module contains code to analyse protein helices using the HELANAL algorithm ([Bansal2000] , [Sugeta1967] ).

HELANAL quantifies the geometry of helices in proteins on the basis of their Cα atoms. It can determine local structural features such as the local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes.

[Sugeta1967]Sugeta, H. and Miyazawa, T. 1967. General method for calculating helical parameters of polymer chains from bond lengths, bond angles and internal rotation angles. Biopolymers 5 673 - 679
[Bansal2000]Bansal M, Kumar S, Velavan R. 2000. HELANAL - A program to characterise helix geometry in proteins. J Biomol Struct Dyn. 17(5):811-819. Example use

You can pass in a single selection:

import MDAnalysis as mda
from MDAnalysis.tests.datafiles import PSF, DCD
from MDAnalysis.analysis import helix_analysis as hel
u = mda.Universe(PSF, DCD)
helanal = hel.HELANAL(u, select='name CA and resnum 161-187')


Alternatively, you can analyse several helices at once by passing in multiple selection strings:

helanal2 = hel.HELANAL(u, select=('name CA and resnum 100-160',
                                  'name CA and resnum 200-230'))

The helix_analysis() function will carry out helix analysis on atom positions, treating each row of coordinates as an alpha-carbon equivalent:

hel_xyz = hel.helix_analysis(u.atoms.positions, ref_axis=[0, 0, 1])