Bioinformatics Programming
dssp_parser: A new Python package to extract helices from DSSP files.
A new Python package named ‘dssp_parser‘ is developed to parse DSSP files. This package fetches all helices including long and short ones from DSSP files.
dssp_parser consists of Python scripts that extract long and short helices in separate files. You have to assign secondary structures using DSSP [1] to the PDB [2] structures from which you want to get helices.
Requirements
This package requires Python3.
Prepare your files
Download the package and save all DSSP files in the ‘dssp_files‘ directory or provide a full path in the scripts if saved elsewhere. Similarly, you can provide the path to your output directories in the scripts.
Usage
$ python3 parse_longH.py
or
$ python3 parse_shortH.py
The program will fetch long and short helices in separate directories. It will name the helix files according to their occurrence in the .dssp file.
Availability
The package is available to download from the GitHub account. The test data is also provided in the given directories.
References
- Frishman, D., & Argos, P. (1995). Knowledge‐based protein secondary structure assignment. Proteins: Structure, Function, and Bioinformatics, 23(4), 566-579.
- https://www.rcsb.org
With the new year, we are going to start with a very simple yet complicated topic (for beginners) in bioinformatics. In this tutorial, we provide a simple code to plot swarm boxplot using matplotlib and seaborn. (more…)
Previously, we provided a similar script to download ligand SMILES from PDB ligand IDs. In this article, we are downloading PDB ligand structures from their corresponding IDs. (more…)
Fetching SMILE strings for a given number of SDF files of chemical compounds is not such a trivial task. We can quickly obtain them using RDKit or OpenBabel. But what if you don’t have SDF files of ligands in the first place? All you have is Ligand IDs from PDB. If they are a few then you can think of downloading SDF files manually but still, it seems time-consuming, especially when you have multiple compounds to work with. Therefore, we provide a Python script that will read all Ligand IDs and fetch their SDF files, and will finally convert them into SMILE strings. (more…)
You must be logged in to post a comment Login