SPECS is a side-chain-orientation-included protein model-native similarity metric for improved evaluation of protein structural models. SPECS stands for Superposition-based Protein Embedded CA SC score. It combines side-chain orientation and global distance based measures in an integrated framework using the united-residue model of polypeptide conformation for computing model-native similarity. SPECS captures both global and local quality aspects when evaluating structural similarity and is sensitive to minute variations in side-chain, thereby being a robust evaluation metric covering a wide range of modeling scenarios and various aspects of structural similarity. SPECS has the value in (0,1], where 1 indicates a perfect match between two structures.
- Linux system: SPECS is tested on on x86_64 redhat linux system
- GCC compiler
You can run SPECS from SPECS web-server. However, for large-scale benchmarkingm we strongly recommend that you install and run SPECS locally. SPECS is a stand-alone application. To download and install SPECS, please type,
$ git clone https://github.com/Bhattacharya-Lab/SPECS.git
$ cd SPECS
$ g++ -o SPECS SPECS.cpp
To run SPECS, type
$ ./SPECS
You will see the following output
**********************************************************************
* SPECS *
* Superposition-based Protein Embedded CA SC score *
* Range of SPECS: *
* 0.0 <= SPECS <= 1.0, higher scores indicate better similarity *
* For comments, please email to bhattacharyad@auburn.edu *
**********************************************************************
Error! Model file must be provided
Usage: ./SPECS -m model -n native
-m model : model pdb file
-n native : native pdb file
-h help : this message
Example commands to run SPECS
$ cd SPECS
$ ./SPECS -m example/sample_model.pdb -n example/sample_native.pdb
SPECS offer dynamic input validation and offers instantaneous output. Upon running the above command, you will see following output
**********************************************************************
* SPECS *
* Superposition-based Protein Embedded CA SC score *
* Range of SPECS: *
* 0.0 <= SPECS <= 1.0, higher scores indicate better similarity *
* For comments, please email to bhattacharyad@auburn.edu *
**********************************************************************
Structure1: example/sample_model.pdb Length = 72
Structure2: example/sample_native.pdb Length = 72
Number of residues in common = 72
RMSD of common residues = 3.323
SPECS = 0.3823 (dCA = 0.4722 rSC = 0.3823 angl = 0.1597 ang2 = 0.3903 tors = 0.1417)
TM-score = 0.6443 (d0 = 2.97)
MaxSub-score = 0.6197 (d0 = 3.50)
GDT-TS-score = 0.6840 %(d<1) = 0.3611 %(d<2) = 0.5556 %(d<4) = 0.8333 %(d<8) = 0.9861
GDT-HA-score = 0.4722 %(d<0.5) = 0.1389 %(d<1) = 0.3611 %(d<2) = 0.5556 %(d<4) = 0.8333
-------- rotation matrix to rotate Chain-1 to Chain-2 ------
i t(i) u(i,1) u(i,2) u(i,3)
1 -23.6888008118 0.7470999956 0.6341999769 0.1993000060
2 12.4891996384 -0.0847999975 -0.2064999938 0.9747999907
3 45.8860015869 0.6593000293 -0.7451000214 -0.1005000025
Superposition in the TM-score: Length(d<5.0) = 62 RMSD = 2.26
(':' denotes the residue pairs of distance < 5.0 Amstrong)
RLALSDAHFRRICQLIYQRAGIVLADHKRDMVYNRLVRRLRALGLDDFGRYLSMLEANQNSAEWQAFINALT
:::::::::::::::::::::: ::::::::::::::::::::::::::::::::: :::: :::
RLALSDAHFRRICQLIYQRAGIVLADHKRDMVYNRLVRRLRALGLDDFGRYLSMLEANQNSAEWQAFINALT
123456789012345678901234567890123456789012345678901234567890123456789012
For question and comment, please contact bhattacharyad@auburn.edu