AAindex: Difference between revisions
Jump to navigation
Jump to search
m (py3) |
(psico) |
||
Line 4: | Line 4: | ||
|author = [[User:Speleo3|Thomas Holder]] | |author = [[User:Speleo3|Thomas Holder]] | ||
|license = BSD | |license = BSD | ||
}} | |||
{{Infobox psico | |||
|module = psico.aaindex | |||
}} | }} | ||
Line 38: | Line 42: | ||
== PyMOL Example == | == PyMOL Example == | ||
Solvent-accessible surface coloring by [https://www.genome.jp/dbget-bin/www_bget?aaindex:KYTJ820101 Hydropathy index (Kyte-Doolittle, 1982)] | |||
<source lang="python"> | <source lang="python"> | ||
run aaindex.py | |||
aaindex2b KYTJ820101 | aaindex2b KYTJ820101 | ||
spectrum b, | spectrum b, white yellow forest | ||
set surface_solvent | |||
show surface | show surface | ||
</source> | </source> |
Latest revision as of 12:23, 22 March 2021
Type | Python Script |
---|---|
Download | aaindex.py |
Author(s) | Thomas Holder |
License | BSD |
This code has been put under version control in the project Pymol-script-repo |
Included in psico | |
Module | psico.aaindex |
---|
Introduction
AAindex is a database of numerical indices representing various physicochemical and biochemical properties of amino acids and pairs of amino acids. See http://www.genome.jp/aaindex/
This script is a python parser for the AAindex flat files which will be downloaded from ftp://ftp.genome.jp/pub/db/community/aaindex/
The script provides two PyMOL commands (but can also be used without PyMOL).
- aaindex2b: Loads numerical indices from aaindex1 as b-factors into your structure
- pmf: Potential of Mean Force (aaindex3)
Python Example
Consider the script is called aaindex.py, it is placed somewhere in your PYTHONPATH and the aaindex flatfiles are found in the current directory.
import aaindex
aaindex.init(path='.')
aaindex.grep('volume')
x = aaindex.get('KRIW790103')
print(x)
print(x.get('A'))
aaindex.init(index='2')
aaindex.grep('blosum')
x = aaindex.get('HENS920102')
print(x.get('A', 'K'))
PyMOL Example
Solvent-accessible surface coloring by Hydropathy index (Kyte-Doolittle, 1982)
run aaindex.py
aaindex2b KYTJ820101
spectrum b, white yellow forest
set surface_solvent
show surface
See Also
- Protscale from rTools does a similar job in coloring by amino acid properties