Garlic Tutorial

Lesson 9 - Venn Diagram


In this simple lesson, you will learn how to prepare the Venn diagram with some residue counts. There is nothing special in this routine and experts will not find it very useful. It is primarily intended for students and other newbies, which are interested in the residue composition of the given protein.

The only information required to prepare the Venn diagram is the protein sequence. In this example, the sequence will be extracted from PDB data. Another approach will be to load the sequence from some FASTA file. The PDB file used here is 1HUC.pdb, but you can replace it with any PDB file which contains the protein structure.

Step 1 - start garlic (if not started already):

garlic

Step 2 - load the structure.

load 1HUC.pdb

Step 3 - discard all hetero atoms:

sel het
sel com

Step 4 - restrict your selection to chains A and B:

res a,b/*/*/*

Step 5 - extract sequence:

seq from 1

Step 6 - draw Venn diagram with residue counts:

venn

Step 7 - discard the structure:

dis all

Step 8 - clear the sequence buffer. You can omit this command.

seq reset