forked from oxpig/ANARCI
-
Notifications
You must be signed in to change notification settings - Fork 0
Antibody Numbering and Antigen Receptor ClassIfication
License
dstcr/ANARCI
Folders and files
Name | Name | Last commit message | Last commit date | |
---|---|---|---|---|
Repository files navigation
ANARCI \\ // Antibody Numbering and Antigen Receptor ClassIfication \\ // || (c) Oxford Protein Informatics Group (OPIG). 2015-20 || Usage: ANARCI -i <inputsequence or fasta file> Requirements: - HMMER3 version 3.1b1 or higher - http://hmmer.janelia.org/ e.g. ANARCI -i Example_sequence_files/12e8.fasta This will number the files in 12e8.fasta with imgt numbering scheme and print to stdout. ANARCI -i Example_sequence_files/sequences.fasta -o Numbered_sequences.anarci -ht hit_tables.txt -s chothia -r ig This will number the files in sequences.fasta with chothia numbering scheme only if they are an antibody chain (ignore TCRs). It will put the numbered sequences in Numbered_sequences.anarci and the alignment statistics in hit_tables.txt ANARCI -i Example_sequence_files/lysozyme.fasta No antigen receptors should be found. The program will just list the names of the sequences. ANARCI -i EVQLQQSGAEVVRSGASVKLSCTASGFNIKDYYIHWVKQRPEKGLEWIGWIDPEIGDTEYVPKFQGKATMTADTSSNTAYLQLSSLTSEDTAVYYCNAGHDYDRGRFPYWGQGTLVTVSA Or just give a single sequence to be numbered. optional arguments: -h, --help show this help message and exit --sequence INPUTSEQUENCE, -i INPUTSEQUENCE A sequence or an input fasta file --outfile OUTFILE, -o OUTFILE The output file to use. Default is stdout --scheme {kabat,aho,wolfguy,imgt,a,c,chothia,i,k,m,w,martin}, -s {kabat,aho,wolfguy,imgt,a,c,chothia,i,k,m,w,martin} Which numbering scheme should be used. i, k, c, m, w and a are shorthand for IMGT, Kabat, Chothia, Martin (Extended Chothia), Wolfguy and Aho respectively. Default IMGT --restrict {ig,tr,heavy,light,H,K,L,A,B} [{ig,tr,heavy,light,H,K,L,A,B} ...], -r {ig,tr,heavy,light,H,K,L,A,B} [{ig,tr,heavy,light,H,K,L,A,B} ...] Restrict ANARCI to only recognise certain types of receptor chains. --csv Write the output in csv format. Outfile must be specified. A csv file is written for each chain type <outfile>_<chain_type>.csv. Kappa and lambda are considered together. --outfile_hits HITFILE, -ht HITFILE Output file for domain hit tables for each sequence. Otherwise not output. --hmmerpath HMMERPATH, -hp HMMERPATH The path to the directory containing hmmer programs. (including hmmscan) --ncpu NCPU, -p NCPU Number of parallel processes to use. Default is 1. --assign_germline Assign the v and j germlines to the sequence. The most sequence identical germline is assigned. --use_species {alpaca,rabbit,rhesus,pig,rat,human,mouse} Use a specific species in the germline assignment. --bit_score_threshold BIT_SCORE_THRESHOLD Change the bit score threshold used to confirm an alignment should be used. Author: James Dunbar ([email protected]) Charlotte Deane ([email protected]) Contact: [email protected] -------------------------------------------------------------------------------------- Output files: The numbering file. The numbering file (--outfile or stdout) reports the numbering for all sequences given in the sequence file. Each record is separated by "//". Those chains for which no significant alignment was found report the name as in the fasta file. e.g: # 1A14:N|PDBID|CHAIN|SEQUENCE // Those sequences where a significant alignment has been found report as below: # 1A14:H|PDBID|CHAIN|SEQUENCE # ANARCI numbered # Domain 1 of 1 # Most significant HMM hit #|species|chain_type|e-value|score|seqstart_index|seqend_index| #|mouse|H|8.6e-58|184.9|0|119| # Scheme = imgt H 1 Q H 2 V H 3 Q H 4 L H 5 Q . . . // where: species = The species of the most significant aligned HMM chain_type = The chain type of the most significant aligned HMM e-value = The e-value of the alignment to the most significant aligned HMM score = The bit-score of the alignment to the most significant aligned HMM seqstart_index = The index in the sequence of the first numbered residue seqend_index = The index in the sequence of the last numbered residue Scheme = The numbering scheme used to number the sequence Then follows the numbering. Chain type (H, L (for both kappa(K) and lambda(L) chain types) , A (alpha), B (beta)) If the "assign_germline" option has been specified the further following lines are added to the header. e.g. # Most sequence-identical germlines #|species|v_gene|v_identity|j_gene|j_identity| #|mouse|IGHV1-12*01|0.86|IGHJ2*01|0.79| where: species = The species of the most sequence identical germline v_gene = The identifier of the most sequence identical germline over the v-region v_identity = The sequence identity over the v-region to the most sequence identical germline j_gene = The identifier of the most sequence identical germline over the j-region j_identity = The sequence identity over the j-region to the most sequence identical germline The csv format output file. When the --csv option is specified, numbered sequences are output into separate comma separated value files depending on their chain type. This provides a horizontal output format and contains all the properties detailed above. In addition, sequences are aligned according to the numbering scheme. The hit file. The hit file reports the statistics for all alignments to each HMM in the database even if the sequence was not numbered. Each record is separated by "//". The corresponding hit table for the numbered entry above looks like: """ NAME 1a14_H mol:protein length:120 NC10 FV (HEAVY CHAIN) SEQUENCE QVQLQQSGAELVKPGASVRMSCKASGYTFTNYNMYWVKQSPGQGLEWIGIFYPGNGDTSYNQKFKDKATLT SEQUENCE ADKSSNTAYMQLSSLTSEDSAVYYCARSGGSYRYDGGFDYWGQGTTVTV id description evalue bitscore bias best_dom_evalue best_dom_bitscore best_dom_bias domain_exp_num domain_obs_num mouse_H 1.1e-57 184.5 1.5 1.3e-57 184.4 1.5 1.0 1 human_H 7.8e-53 169.0 1.9 8.7e-53 168.8 1.9 1.0 1 rat_H 4.7e-47 150.2 2.2 5.2e-47 150.0 2.2 1.0 1 rabbit_H 3.7e-37 118.2 0.7 4e-37 118.1 0.7 1.0 1 pig_H 1.5e-35 113.3 2.7 1.6e-35 113.1 2.7 1.0 1 rhesus_H 4.4e-32 101.5 1.8 4.9e-32 101.4 1.8 1.0 1 mouse_B 2.4e-19 60.6 0.7 2.6e-19 60.5 0.7 1.0 1 human_B 4.2e-19 59.7 0.9 4.6e-19 59.5 0.9 1.0 1 mouse_A 8.7e-19 58.5 1.1 9.6e-19 58.4 1.1 1.0 1 human_A 1.7e-18 57.6 0.9 1.9e-18 57.5 0.9 1.0 1 mouse_D 5.1e-17 53.3 0.7 5.9e-17 53.1 0.7 1.1 1 rhesus_L 1.6e-16 51.7 2.7 1.9e-16 51.4 2.7 1.1 1 human_L 1.7e-15 48.3 3.5 2e-15 48.0 3.5 1.1 1 human_D 6.7e-15 46.1 0.2 7.4e-15 45.9 0.2 1.0 1 rhesus_K 3.9e-13 40.6 1.7 5.1e-13 40.2 1.7 1.2 1 mouse_G 4.1e-13 40.3 0.0 4.3e-13 40.2 0.0 1.0 1 rabbit_L 6.1e-13 40.0 2.8 8.1e-13 39.6 2.8 1.2 1 rat_K 3.9e-12 37.4 1.4 4.4e-12 37.2 1.4 1.1 1 pig_L 4.2e-12 37.5 1.0 4.7e-12 37.3 1.0 1.1 1 mouse_K 1.2e-11 35.7 2.6 1.3e-11 35.6 2.6 1.1 1 human_K 2.2e-11 34.8 2.9 3.5e-11 34.2 2.9 1.4 1 mouse_L 1.9e-10 31.8 2.2 3.4e-10 30.9 2.2 1.4 1 rat_L 2.5e-10 31.7 1.2 2.9e-10 31.5 1.2 1.1 1 pig_K 3.2e-10 31.1 1.9 4.5e-10 30.6 1.9 1.3 1 human_G 2.9e-09 27.8 0.8 4.9e-09 27.1 0.8 1.4 1 rabbit_K 2.5e-06 18.4 5.8 4.2e-06 17.7 5.8 1.4 1 // """ We therefore get a ranking of the alignments to each chain type. -------------------------------------------------------------------------------------- Schemes: Currently implemented schemes: IMGT Chothia (IGs only) Kabat (IGs only) Martin / Enhanced Chothia (IGs only) AHo Wolfguy (IGs only) Currently recognisable species (chains): Human (heavy, kappa, lambda, alpha, beta) Mouse (heavy, kappa, lambda, alpha, beta) Rat (heavy, kappa, lambda) Rabbit (heavy, kappa, lambda) Pig (heavy, kappa, lambda) Rhesus Monkey (heavy, kappa) Other species may still be numbered correctly and the chain type recognised but the species be incorrect. e.g. llama VHH. IMGT - has 128 possible positions for *all* antigen receptor types. These are supposed to be structurally equivalent. In theory these are supposed to account for all possible positions. However, insertions are possible especially at CDR3. ANARCI gives the insertion codes as letters. Insertions at CDR3 are placed symmetrically about imgt positions 111 and 112. e.g. 111-ABCD DCBA-112. Kabat - is defined for heavy and light chain antibody chains only (in ANARCI). Positions in the two chain types are not equivalent. Insertions occur at specific positions and can occur in both the framework and the CDRs. They are annotated from A->Z. e.g 100ABCDEFGH 101. Chothia - is defined for heavy and light chain antibody chains only (in ANARCI). Numbering in the two chain types are not equivalent. Insertions occur at specific positions and can occur in both the framework and the CDRs. They are annotated from A->Z. e.g 100ABCDEFGH 101. It differs to the Kabat scheme by the position it places the insertions at CDRH1. Martin - is defined for heavy and light chain antibody chains only. Numbering in the two chain types are not equivalent. Insertions occur at specific positions and can occur in both the framework and the CDRs. They are annotated from A->Z. e.g 100ABCDEFGH 101. It differs to the Chothia scheme by the position it places the certain insertions in the framework. It is also referred to as the enhanced Chothia scheme. AHo - has 149 possible for *all* antigen receptor types. These are supposed to be structurally equivalent. The AHo scheme's large number of positions is supposed to account for all possible positions without the need for specifying insertion positions. In ANARCI, extra residues in the framework may be represented by insertions although these are unlikely to occur in natural sequences. Wolfguy - is defined for heavy and light antibody chains. Numbering in the two chain types are not equivalent. Different regions of the domain are denoted by a range of numbers. Many possible positions in the CDRs mean that insertion codes are not required. In ANARCI, extra residues in the framework may be represented by insertions although these are unlikely to occur in natural sequences. The CDRs are numbered in an 'up' and 'down' direction. The annotations of CDRL1 is defined according to the canonical structure. In ANARCI this is recognised by taking a sequence similarity to hard coded sequence motifs for different lengths. -------------------------------------------------------------------------------------
About
Antibody Numbering and Antigen Receptor ClassIfication
Resources
License
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published
Languages
- Python 99.5%
- Shell 0.5%