how to identify CDR region in antibody sequence
2
3
Entering edit mode
12 days ago
reany ▴ 50

I want to extract CDR region form an antibody sequence or numbered antibody sequence. Because SCALOP will miss H-CDR3, is there annother tool could identify CDR region? After numbering the antibody sequence by ANARCI, whether it's the right way to extract CDR region according to the position of CDR definitions?

I am looking for some stand-alone tools or correct rules to extract CDR after numbering antibody sequence. I am not sure the definition is enough to do this.

Any other suggestions would be appreciated.

CDR antobidy • 433 views
ADD COMMENT
0
Entering edit mode

What species are you analyzing? IMGT is a great resource.

ADD REPLY
1
Entering edit mode

Sorry about that the problem was not clearly described. Yes IMGT is useful, but i am looking for some stand-alone tools or correct rules to extract CDR after numbering antibody sequence. I am not sure the definition is enough to do this.

ADD REPLY
1
Entering edit mode
11 days ago
Jeremy ▴ 340

I would try MIXCR. You can look for other tools on b-t.cr under the Software category. If your main focus is antibodies, you could also consider joining the AIRR Slack channel and posting a question there.

CDR H3 starts one amino acid after the conserved cysteine at the end of the VH region and ends 1 amino acid before the conserved tryptophan at the beginning of the JH region. (See Figure 5 in the paper below.)

Antibody Diversity Paper

ADD COMMENT
1
Entering edit mode

MIXCR seems like situed for massively parallel sequencing data while my input is relatively simple, for example, just a protein sequence of light chain. I'm still going to look into MIXCR and hopefully learn something from it. Other suggestions are also ueful and thanks very much.

ADD REPLY
0
Entering edit mode

If you can provide an example sequence, I can probably write some R code to detect CDR H3.

ADD REPLY
1
Entering edit mode

The AIRR folks are on slack? I had no idea. Thanks for mentioning that.

ADD REPLY
0
Entering edit mode
11 days ago
Jesse ▴ 90

I generally just use IgBLAST when I need antibody sequence annotations like CDR3. There's a command-line version and you can have arbitrary species and gene references (though in that case you need to jump through some hoops and create an "auxiliary data file" to get it to report the CDR3 info). It can also give AIRR-compatible TSV output so you can extract sequences directly from the specific columns you want, like cdr3, without worrying about a particular numbering scheme (e.g. kabat) and extracting subsequences yourself, but there are a ton of columns with position info as well. The web interface uses IMGT references by default, too. IMGT's own V-QUEST tool can give similar info, but is web-only and I don't believe supports custom references. There are a bunch of tools out there for bulk data but IgBLAST scales up and down nicely to even just one or a few sequences.

Also I'd reiterate what Jeremy said about the junction (where the sequence for the conserved amino acids at each end are included) versus the CDR3 (which leaves out those bits). Some texts jumble those definitions and that's tripped me up before.

ADD COMMENT
0
Entering edit mode

For igblastp, auxiliary_data is not support to get accurate CDR3 aligment although it will report H-CDR3 sometimes. Is it a challenge for existing tools to delimit H-CDR3 or have other reasons?

ADD REPLY
0
Entering edit mode

Oh, sorry, I haven't tried it for amino acid sequences so I don't know about igblastp's behavior. I'm a bit surprised any of these tools would have much trouble labeling CDR3 though I suppose there's less to go by with just the amino acid sequence. The IgBLAST docs say igblastp doesn't search D and J which matches what you're saying. Do you have nucleotide sequence you could use too, or just amino acid?

ADD REPLY
0
Entering edit mode

Only protein sequence aviliable. By changing some settings, H-CDR3 can be report with SCALOP finally although that may be inaccurate. Thanks anyway.

ADD REPLY

Login before adding your answer.

Traffic: 817 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6