|
|
DOMAINNR documentation |
ID D1CS4A_
XX
EN 1CS4
XX
TY SCOP
XX
SI 53931 CL; 54861 FO; 55073 SF; 55074 FA; 55077 DO; 55078 SO; 39418 DD;
XX
CL Alpha and beta proteins (a+b)
XX
FO Ferredoxin-like
XX
SF Adenylyl and guanylyl cyclase catalytic domain
XX
FA Adenylyl and guanylyl cyclase catalytic domain
XX
DO Adenylyl cyclase VC1, domain C1a
XX
OS Dog (Canis familiaris)
XX
DS SEQUENCE 52 AA; 5817 MW; D8CCAE0E1FC0849A CRC64;
ADIEGFTSLA SQCTAQELVM TLNELFARFD KLAAENHCLR IKILGDCYYC VS
XX
NC 1
XX
CN [1]
XX
CH A CHAIN; . START; . END;
//
ID D1II7A_
XX
EN 1II7
XX
TY SCOP
XX
SI 53931 CL; 56299 FO; 56300 SF; 64427 FA; 64428 DO; 64429 SO; 62415 DD;
XX
CL Alpha and beta proteins (a+b)
XX
FO Metallo-dependent phosphatases
XX
SF Metallo-dependent phosphatases
XX
FA DNA double-strand break repair nuclease
XX
DO Mre11
XX
OS Archaeon Pyrococcus furiosus
XX
DS SEQUENCE 65 AA; 7396 MW; 75FBE75B22FD3678 CRC64;
MKFAHLADIH LGYEQFHKPQ REEEFAEAFK NALEIAVQEN VDFILIAGDL FHSSRPSPGT
LKKAI
XX
NC 1
XX
CN [1]
XX
CH A CHAIN; . START; . END;
//
|
ID D1CS4A_
XX
EN 1CS4
XX
TY SCOP
XX
SI 53931 CL; 54861 FO; 55073 SF; 55074 FA; 55077 DO; 55078 SO; 39418 DD;
XX
CL Alpha and beta proteins (a+b)
XX
FO Ferredoxin-like
XX
SF Adenylyl and guanylyl cyclase catalytic domain
XX
FA Adenylyl and guanylyl cyclase catalytic domain
XX
DO Adenylyl cyclase VC1, domain C1a
XX
OS Dog (Canis familiaris)
XX
DS SEQUENCE 52 AA; 5817 MW; D8CCAE0E1FC0849A CRC64;
ADIEGFTSLA SQCTAQELVM TLNELFARFD KLAAENHCLR IKILGDCYYC VS
XX
NC 1
XX
CN [1]
XX
CH A CHAIN; . START; . END;
//
ID D1II7A_
XX
EN 1II7
XX
TY SCOP
XX
SI 53931 CL; 56299 FO; 56300 SF; 64427 FA; 64428 DO; 64429 SO; 62415 DD;
XX
CL Alpha and beta proteins (a+b)
XX
FO Metallo-dependent phosphatases
XX
SF Metallo-dependent phosphatases
XX
FA DNA double-strand break repair nuclease
XX
DO Mre11
XX
OS Archaeon Pyrococcus furiosus
XX
DS SEQUENCE 65 AA; 7396 MW; 75FBE75B22FD3678 CRC64;
MKFAHLADIH LGYEQFHKPQ REEEFAEAFK NALEIAVQEN VDFILIAGDL FHSSRPSPGT
LKKAI
XX
NC 1
XX
CN [1]
XX
CH A CHAIN; . START; . END;
//
|
ID D1II7A_
XX
EN 1II7
XX
TY SCOP
XX
SI 53931 CL; 56299 FO; 56300 SF; 64427 FA; 64428 DO; 64429 SO; 62415 DD;
XX
CL Alpha and beta proteins (a+b)
XX
FO Metallo-dependent phosphatases
XX
SF Metallo-dependent phosphatases
XX
FA DNA double-strand break repair nuclease
XX
DO Mre11
XX
OS Archaeon Pyrococcus furiosus
XX
DS SEQUENCE 65 AA; 7396 MW; 75FBE75B22FD3678 CRC64;
MKFAHLADIH LGYEQFHKPQ REEEFAEAFK NALEIAVQEN VDFILIAGDL FHSSRPSPGT
LKKAI
XX
NC 1
XX
CN [1]
XX
CH A CHAIN; . START; . END;
//
|
Classes are non-redundant 5% redundancy threshold // Alpha and beta proteins (a+b) Retained D1II7A_ Rejected D1CS4A_ |
Standard (Mandatory) qualifiers (* if not always prompted):
[-dcfinfile] infile This option specifies name of DCF file
(domain classification file) (input). A
'domain classification file' contains
classification and other data for domains
from SCOP or CATH, in DCF format
(EMBL-like). The files are generated by
using SCOPPARSE and CATHPARSE. Domain
sequence information can be added to the
file by using DOMAINSEQS.
-retain toggle [N] This option specifies whether to write
redundant domains to a separate file. If
this option is selected, redundant domains
are written to a separate output file.
-node menu [1] This option specifies the node for
redundancy removal. Redundancy can be
removed at any specified node in the SCOP or
CATH hierarchies. For example by selecting
'Class' entries belonging to the same Class
will be non-redundant. (Values: 1 (Class
(SCOP)); 2 (Fold (SCOP)); 3 (Superfamily
(SCOP)); 4 (Family (SCOP)); 5 (Class
(CATH)); 6 (Architecture (CATH)); 7
(Topology (CATH)); 8 (Homologous Superfamily
(CATH)); 9 (Family (CATH)))
-mode menu [1] This option specifies whether to remove
redundancy at a single threshold % sequence
similarity or remove redundancy outside a
range of acceptable threshold % similarity.
All permutations of pair-wise sequence
alignments are calculated for each domain
family in turn using the EMBOSS
implementation of the Needleman and Wunsch
global alignment algorithm. Redundant
sequences are removed in one of two modes as
follows: (i) If a pair of proteins achieve
greater than a threshold percentage sequence
similarity (specified by the user) the
shortest sequence is discarded. (ii) If a
pair of proteins have a percentage sequence
similarity that lies outside an acceptable
range (specified by the user) the shortest
sequence is discarded. (Values: 1 (Remove
redundancy at a single threshold % sequence
similarity); 2 (Remove redundancy outside a
range of acceptable threshold % similarity))
* -threshold float [95.0] This option specifies the % sequence
identity redundancy threshold, which
determines the redundancy calculation. If a
pair of proteins achieve greater than this
threshold the shortest sequence is
discarded. (Any numeric value)
* -threshlow float [30.0] This option specifies the % sequence
identity redundancy threshold, which
determines the redundancy calculation. If a
pair of proteins have a percentage sequence
similarity that lies outside an acceptable
range the shortest sequence is discarded.
(Any numeric value)
* -threshup float [90.0] This option specifies the % sequence
identity redundancy threshold, which
determines the redundancy calculation. If a
pair of proteins have a percentage sequence
similarity that lies outside an acceptable
range the shortest sequence is discarded.
(Any numeric value)
[-dcfoutfile] outfile [test.scop] This option specifies the name
of non-redundant DCF file (domain
classification file) (output). A 'domain
classification file' contains classification
and other data for domains from SCOP or
CATH, in DCF format (EMBL-like). The files
are generated by using SCOPPARSE and
CATHPARSE. Domain sequence information can
be added to the file by using DOMAINSEQS.
* -redoutfile outfile [*.domainnr] This option specifies the name
of DCF file (domain classification file) for
redundant sequences (output). A 'domain
classification file' contains classification
and other data for domains from SCOP or
CATH, in DCF format (EMBL-like). The files
are generated by using SCOPPARSE and
CATHPARSE. Domain sequence information can
be added to the file by using DOMAINSEQS.
-logfile outfile [domainnr.log] This option specifies the
name of log file for the build. The log file
contains messages about any errors arising
while domainnr ran.
Additional (Optional) qualifiers:
-datafile matrixf [EBLOSUM62] This option specifies the
residue substitution matrix. This is used
for sequence comparison.
-gapopen float [10.0 for any sequence] This option
specifies the gap insertion penalty. This is
the score taken away when a gap is created.
The best value depends on the choice of
comparison matrix. The default value assumes
you are using the EBLOSUM62 matrix for
protein sequences, and the EDNAFULL matrix
for nucleotide sequences. (Floating point
number from 1.0 to 100.0)
-gapextend float [0.5 for any sequence] This option specifies
the gap extension penalty. This is added to
the standard gap penalty for each base or
residue in the gap. This is how long gaps
are penalized. Usually you will expect a few
long gaps rather than many short gaps, so
the gap extension penalty should be lower
than the gap penalty. (Floating point number
from 0.0 to 10.0)
Advanced (Unprompted) qualifiers: (none)
Associated qualifiers:
"-dcfoutfile" associated qualifiers
-odirectory2 string Output directory
"-redoutfile" associated qualifiers
-odirectory string Output directory
"-logfile" associated qualifiers
-odirectory string Output directory
General qualifiers:
-auto boolean Turn off prompts
-stdout boolean Write standard output
-filter boolean Read standard input, write standard output
-options boolean Prompt for standard and additional values
-debug boolean Write debug output to program.dbg
-verbose boolean Report some/full command line options
-help boolean Report command line options. More
information on associated and general
qualifiers can be found with -help -verbose
-warning boolean Report warnings
-error boolean Report errors
-fatal boolean Report fatal errors
-die boolean Report dying program messages
| Standard (Mandatory) qualifiers | Allowed values | Default | |||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| [-dcfinfile] (Parameter 1) |
This option specifies name of DCF file (domain classification file) (input). A 'domain classification file' contains classification and other data for domains from SCOP or CATH, in DCF format (EMBL-like). The files are generated by using SCOPPARSE and CATHPARSE. Domain sequence information can be added to the file by using DOMAINSEQS. | Input file | Required | ||||||||||||||||||
| -retain | This option specifies whether to write redundant domains to a separate file. If this option is selected, redundant domains are written to a separate output file. | Toggle value Yes/No | No | ||||||||||||||||||
| -node | This option specifies the node for redundancy removal. Redundancy can be removed at any specified node in the SCOP or CATH hierarchies. For example by selecting 'Class' entries belonging to the same Class will be non-redundant. |
|
1 | ||||||||||||||||||
| -mode | This option specifies whether to remove redundancy at a single threshold % sequence similarity or remove redundancy outside a range of acceptable threshold % similarity. All permutations of pair-wise sequence alignments are calculated for each domain family in turn using the EMBOSS implementation of the Needleman and Wunsch global alignment algorithm. Redundant sequences are removed in one of two modes as follows: (i) If a pair of proteins achieve greater than a threshold percentage sequence similarity (specified by the user) the shortest sequence is discarded. (ii) If a pair of proteins have a percentage sequence similarity that lies outside an acceptable range (specified by the user) the shortest sequence is discarded. |
|
1 | ||||||||||||||||||
| -threshold | This option specifies the % sequence identity redundancy threshold, which determines the redundancy calculation. If a pair of proteins achieve greater than this threshold the shortest sequence is discarded. | Any numeric value | 95.0 | ||||||||||||||||||
| -threshlow | This option specifies the % sequence identity redundancy threshold, which determines the redundancy calculation. If a pair of proteins have a percentage sequence similarity that lies outside an acceptable range the shortest sequence is discarded. | Any numeric value | 30.0 | ||||||||||||||||||
| -threshup | This option specifies the % sequence identity redundancy threshold, which determines the redundancy calculation. If a pair of proteins have a percentage sequence similarity that lies outside an acceptable range the shortest sequence is discarded. | Any numeric value | 90.0 | ||||||||||||||||||
| [-dcfoutfile] (Parameter 2) |
This option specifies the name of non-redundant DCF file (domain classification file) (output). A 'domain classification file' contains classification and other data for domains from SCOP or CATH, in DCF format (EMBL-like). The files are generated by using SCOPPARSE and CATHPARSE. Domain sequence information can be added to the file by using DOMAINSEQS. | Output file | test.scop | ||||||||||||||||||
| -redoutfile | This option specifies the name of DCF file (domain classification file) for redundant sequences (output). A 'domain classification file' contains classification and other data for domains from SCOP or CATH, in DCF format (EMBL-like). The files are generated by using SCOPPARSE and CATHPARSE. Domain sequence information can be added to the file by using DOMAINSEQS. | Output file | <*>.domainnr | ||||||||||||||||||
| -logfile | This option specifies the name of log file for the build. The log file contains messages about any errors arising while domainnr ran. | Output file | domainnr.log | ||||||||||||||||||
| Additional (Optional) qualifiers | Allowed values | Default | |||||||||||||||||||
| -datafile | This option specifies the residue substitution matrix. This is used for sequence comparison. | Comparison matrix file in EMBOSS data path | EBLOSUM62 | ||||||||||||||||||
| -gapopen | This option specifies the gap insertion penalty. This is the score taken away when a gap is created. The best value depends on the choice of comparison matrix. The default value assumes you are using the EBLOSUM62 matrix for protein sequences, and the EDNAFULL matrix for nucleotide sequences. | Floating point number from 1.0 to 100.0 | 10.0 for any sequence | ||||||||||||||||||
| -gapextend | This option specifies the gap extension penalty. This is added to the standard gap penalty for each base or residue in the gap. This is how long gaps are penalized. Usually you will expect a few long gaps rather than many short gaps, so the gap extension penalty should be lower than the gap penalty. | Floating point number from 0.0 to 10.0 | 0.5 for any sequence | ||||||||||||||||||
| Advanced (Unprompted) qualifiers | Allowed values | Default | |||||||||||||||||||
| (none) | |||||||||||||||||||||
% domainnr
Removes redundant domains from a DCF file.
Domain classification file: ../domainseqs-keep/all_s.scop
Write redundant domains to separate file. [N]:
Node at which to remove redundancy
1 : Class (SCOP)
2 : Fold (SCOP)
3 : Superfamily (SCOP)
4 : Family (SCOP)
5 : Class (CATH)
6 : Architecture (CATH)
7 : Topology (CATH)
8 : Homologous Superfamily (CATH)
9 : Family (CATH)
Select number. [1]: 1
Redundancy removal options
1 : Remove redundancy at a single threshold % sequence similarity
2 : Remove redundancy outside a range of acceptable threshold % similarity
Select number. [1]: 1
The % sequence identity redundancy threshold. [95.0]: 5
Domain classification output file [test.scop]: all_nr.scop
Domainatrix log output file [domainnr.log]:
Warning: Bad args passed to ajDomainWrite
// Alpha and beta proteins (a+b)
D1CS4A_
D1II7A_
|
Go to the input files for this example
Go to the output files for this example
| FILE TYPE | FORMAT | DESCRIPTION | CREATED BY | SEE ALSO |
| Domain classification file (for SCOP) | DCF format (EMBL-like format for domain classification data). | Classification and other data for domains from SCOP. | SCOPPARSE | Domain sequence information can be added to the file by using DOMAINSEQS. |
| Domain classification file (for CATH) | DCF format (EMBL-like format for domain classification data). | Classification and other data for domains from CATH. | CATHPARSE | Domain sequence information can be added to the file by using DOMAINSEQS. |
| Program name | Description |
|---|---|
| aaindexextract | Extract data from AAINDEX |
| allversusall | Sequence similarity data from all-versus-all comparison |
| cathparse | Generates DCF file from raw CATH files |
| cutgextract | Extract data from CUTG |
| domainer | Generates domain CCF files from protein CCF files |
| domainseqs | Adds sequence records to a DCF file |
| domainsse | Add secondary structure records to a DCF file |
| hetparse | Converts heterogen group dictionary to EMBL-like format |
| pdbparse | Parses PDB files and writes protein CCF files |
| pdbplus | Add accessibility & secondary structure to a CCF file |
| pdbtosp | Convert swissprot:PDB codes file to EMBL-like format |
| printsextract | Extract data from PRINTS |
| prosextract | Build the PROSITE motif database for use by patmatmotifs |
| rebaseextract | Extract data from REBASE |
| scopparse | Generate DCF file from raw SCOP files |
| seqnr | Removes redundancy from DHF files |
| sites | Generate residue-ligand CON files from CCF files |
| ssematch | Search a DCF file for secondary structure matches |
| tfextract | Extract data from TRANSFAC |
Excerpt from DOMAINNR log file Families are non-redundant 95% redundancy threshold // Homeodomain Retained D2HDDA_ D1AKHA_ D1MNMC_ Rejected D2HDDB_ D1ENH__ D3HDDA_ WARN d3hdda_.pxyz not found // Di-haem cytohrome c peroxidase WARN ds005__.pxyz not found WARN Empty family // Nuclear receptor coactivator Src-1 Retained D2PRGC_ Rejected |
See also http://emboss.sourceforge.net/