This script will calculate the solvent-accessible surface (SASA) for each amino acid within a protein and point out which amino acids are in the wrong layer.
This script does the following:
- Prints out PYMOL commands to select spesific amino acids in a protein depending on their layer calculated by their SASA (solvent-accessible surface area). These commands allow the user to colour different parts of a protein and better visualise it in order to take better decisions on refining it.
- Prints out a sequence alignment which shows which amino acid belongs to which layers in the protein accodring to its SASA value, in which secondary structure does each amino acid belong to, the protein's amino acid sequence, and points out with a * which amino acids belong to the wrong layer and should be mutated.
This script was written to run on GNU/Linux, it was not tested in Windows or MacOS. This script will mostly be useful to refine proteins after a Rosetta FFL (Fold From Loop) computation, but can still be used to refine any protein. Contact the author at sari.sabban@gmail.com for any questions regarding this script.
To use in python 2 follow these steps:
- Install biopython by running the following command in terminal (pip install biopython) you need pip to be installed, if you do not have it you can install it in linux by running the following command in terminal (sudo apt-get install python-pip).
- Install DSSP in linux by running the following command in terminal (sudo apt-get install dssp).
- Run by typeing in command line:
python2 RefineAssist.py FILENAME.pdb
- File must be in the same directory as this script.
To use in python 3 follow these steps:
- Install biopython by running the following command in terminal (python3 -m pip install biopython) you need pip to be installed, if you do not have it you can install it in linux by running the following command in terminal (sudo apt-get install python3-pip).
- Install DSSP in linux by running the following command in terminal (sudo apt-get install dssp).
- Install numpy (python3 -m pip install numpy).
- Run by typeing in command line:
python3 RefineAssist3.py FILENAME.pdb
- File must be in the same directory as this script.
- Refere to the paper by (Koga et.al., 2012 - PMID: 23135467) Methods section's Sequence Design Protocol for more explanation on protein refinement and each layer's SASA parameters.
- Refere to the paper by (Correia et.al., 2014 - PMID: 24499818) for details about the Rosetta Fold From Loop (FFL) protocol.
- Refere to this webpage for details about Rosetta LayerDesign Protocol and which amino acids should be in which layer.
- Refere to this webpage for how to use BioPython.
- Refere to the references (Cock et.al., 2009 - PMID: 19304878) and (Hamelryck and Manderick, 2003 - PMID: 14630660) regarding the applications used here from biopython.
- Refere to the references (Kabsch W. and Sander C., 1983 - PMID: 6667333) regading DSSP.
- Refere to the reference (Tien et.al., 2013 - PMID: 24278298) regarding Wilke SASA parameters.