Difference between revisions of "Annotate v"
(A simple-to-use script for annotation of VH or VL sequence of an antibody.)
|Line 57:||Line 57:|
Revision as of 00:36, 8 February 2020
|This code has been put under version control in the project Pymol-script-repo|
A simple-to-use script for annotation of VH or VL sequence of an antibody. It creates a Pymol object for each FR or CDR region. It utilizes the REST API interface of Abnum (http://www.bioinf.org.uk/abs/abnum/) from Dr Andrew Martin's group at UCL
Currently supports Kabat, Chothia, Contact, IMGT1
1Definitions for Kabat, Chothia, Contact, and IIMGT are same as listed in the table (http://www.bioinf.org.uk/abs/info.html#kabatnum), except that for IMGT, H-CDR2 is defined as H51-H57 in this script, as opposed to of H51-H56 in the table. This slight revision generates result that matches that from IMGT website (http://www.imgt.org/)
Dependencies and Limitations
1. Import request module
2. Relies on internet connection to Abnum
3. Incomplete VH or VL sequence might not be annotated
How to use
1. Create a selection object for the V region, (such as VH and VL shown here). Only 1 V-region (VH OR VL) can be annotated each time. Alternatively, simply select the residues into the default <sele> group.
2. Copy the entire script and create a `.py` file. Run script.
3. The general syntax is:
selection_group_name: name of the selection group where you have the input sequence. Must be single-letter amino acid coded. Either VH or VL but not both
scheme: currently supports Kabat, Chothia, Contact, IMGT. Must be lowercase
annotate_v("VH", "kabat") #input sequence from VH selection group, annotate using "kabat" annotate_v("sele", "imgt") #input sequence from the default sele group, annotate using "imgt". You can also try "contact", "chothia"
4. In the output window, the script will print the FR and CDR regions (if found). It also automatically creates a selection group for each of the FR and CDR region.
Bugs? Questions? Email email@example.com