220.127.116.11. HELANAL — analysis of protein helices¶
|Copyright:||GNU Public License v3|
New in version 1.0.0.
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.|
18.104.22.168.1. 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') helanal.run() print(helanal.summary)
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'))
helix_analysis() function will carry out helix analysis on
atom positions, treating each row of coordinates as an alpha-carbon
hel_xyz = hel.helix_analysis(u.atoms.positions, ref_axis=[0, 0, 1])