Professional Documents
Culture Documents
Coursera BioinfoMethods-II Lab02
Coursera BioinfoMethods-II Lab02
Lab 2
In practice, there are several problems with the yeast two hybrid system. The protein hybrids
must be targetted to the yeast nucleus, so membrane-bound and membrane-associated
1
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
interacting proteins can seldom be identified. Additionally, many proteins are inherently
sticky so the rate of false positives can be quite high this is exacerbated because the hybrid
proteins are typically overexpressed in the yeast cells. Adaptations to the original method have
been devised to obviate some of these problems, and on the plus side, the Y2H system can be
quite sensitive to transient interactions.
To identify protein-protein interactions in their endogenous
context, affinity purification may be used. Here antibodies to
a given protein of interest can be used to capture that protein
and its interactors from cell extracts. Proteins that copurify
with a given target protein are then identified by mass
spectrometry. A variation of this method is the tandem
affinity purification (TAP) tag method, developed in 1999 by
Bertrand Sraphin and colleagues at the EMBL Laboratories
in Heidelberg. A schematic of this method is shown to the
right. Basically, two affinity tags are attached to the protein
of interest, expressed under the control of the native or some
other promoter. These tags are the calmodulin binding
peptide and two IgG binding domains of Protein A from S.
aureus. These are separated by a TEV protease cleavage site. Figure from Rigaut et al. (1999) Nature
Biotechnology 17:1030-1032.
Two rounds of affinity purification are thus possible,
resulting in far fewer false positive interactions. The resultant interacting proteins are then
identified by mass spectrometry. Drawbacks include the fact that the introduced tags may
disrupt potential interactions. In general, the TAP tag method is less sensitive in detecting
transient interactions and is better for identifying proteins in protein complexes.
One final method of note is the identification of potential protein-protein interactions by
orthology. In this method, proteins from a given species whose orthologs have been
identified as interactors can also be assumed to interact. These are sometimes called
interologs, for interacting orthologs. The confidence that interologs are truly interacting can
be increased if the interaction is seen in orthologous proteins across several species.
Other high throughput methods have been developed for identifying interactions between
membrane proteins, and other proteins that are considered difficult to work with. It is
important to recognize the limitations of each of these methods and to be aware of the quality
of PPI data in the databases. Multiple lines of interaction evidence are desirable.
1. DIP
Go to the Database of Interacting Proteins http://dip.doe-mbi.ucla.edu/dip/Main.cgi and click on
SEARCH in the navigation bar on the left. Click on the Node search type and enter BRCA2
in the Name/Description box and select Homo sapiens as the organism under the Node
Annotation search. Then Query DIP using the lower Query button.
2
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Aside: Recall from last weeks lab that BRCA2 is the Breast Cancer Type 2 susceptibility
protein. BRCA2 interacts with RAD51 (among others) in the DNA damage and repair response
pathway where both proteins are critical to its proper function. Mutations in several of the DNA
damage response pathway proteins, including BRCA1, BRCA2 and RAD51 have been linked to
multiple forms of cancer including breast cancer.
When the search is completed, click on the DIP reference number link for the BRCA2 node
(24214N). Click graph in the top right corner of the DIP Node window. Select nodes in the
graph by clicking on them. BRCA2 is shown in red.
Figure 1. Graph output of DIP showing BRCA2 interactions. BRCA2 is coloured red. See Box 2
for an explanation of how to interpret this output.
a. Which proteins have been identified in DIP as interacting with BRCA2 (list some exemplars)?
b. Which interacting protein with BRCA2 has the most identified interactions with other
proteins?
Lab Quiz
Question 1
c. What is this proteins function? (Hint: click on it to see its record in DIP; explore links there.)
3
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Go back to the Node Search Results page and click on the dot under DIP Links. This
provides you with a table of identified interactors with BRCA2. If you click on the Interaction
entries (e.g. DIP:57452E), you can see how these interactions were identified. Check out both the
Binary and Complex tabs.
d. How were the interactions with BRCA2 identified? Do you believe the data?
e. Why are there three edge (interaction) entries for RAD51 in the Binary section?
2. BioGRID
BioGRID (General Respository for Interaction Datasets) is a curated database of over 812,935
non-redundant physical and genetic interactions in dozens species. It was created and is
maintained by Mike Tyers laboratory, formerly in Toronto. Connect to BioGRID at
http://www.thebiogrid.org/index.php. Search for interactions with BRCA2 in Homo sapiens.
Figure 2. Partial output from BioGRID for BRCA2 interactors. Click on the interactor names
(e.g. RAD51) to see the Gene Ontology categories associated with these interactors.
4
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
a. What types of methods have been used to determine protein-protein interactions with
BRCA2? Hint: check out the [details] link.
Lab 2
Lab Quiz
Question 2
5
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Image courtesy of ChaTo, under the GNU Free Documentation License, Version
1.2.
Bioinformatic Methods II
Lab 2
Figure 4: Importing a network into Cytoscape from BioGRID via web services. Search for the
desired term in Step 1, select the appropriate database in Step 2 and finally click on Import.
Configure the Import Network from Web Service dialogue box as shown in Figure 4: type
BRCA2 into the 1. Enter Search Conditions search box, and click Search. In 2. Select
Databases select BioGRID. Click Import and then No on the Manually Merge Networks
dialogue box that will appear after you do this. Close the Import Networks dialogue box.
Explore Cytoscapes interface. You can zoom in on the network that you retrieve by clicking the
magnifying glass icon in the tool bar along the top of the screen. If you click on a node, you will
see that information about that node will appear in the Table Panel at the bottom of the screen
(Figure 5). You will only see a few columns of data, but you can easily add other columns by
clicking on the Show All Columns button:
. You can select multiple nodes by holding the
shift key and clicking the other desired nodes, or by holding the left mouse key and drawing a
box around the nodes of interest. You can also explore different layout options for the network
using the Layout menu option the yFiles Organic layout is shown in Figure 5. You can also
select specific edges (which denote interactions) by clicking on them and then switching to the
Edge Table tab. Unfortunately, we dont get a lot of information about how the interactions were
determined from the BioGRID web service! Its a good thing we explored these in the BioGRID
web interface.
7
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Figure 5: BRCA2 protein-protein interaction network, retrieved from the BioGRID Web Service
Client. BRCA2 has been selected by clicking on it (yellow node). The Layout was set using the
yFiles > Organic layout, and the information for the BRCA2 node retrieved by the web service
call is shown in the Table Panel below the network diagram. The Show All Columns button
was used to display all information about the selected node.
a. How many non-human interactors did we retrieve from BioGRID for BRCA2 and what
organism(s) are they from? (Hint: the default colour scheme for the nodes is by NCBIs
Taxonomy ID; you can find out the corresponding organism by going to
http://www.ncbi.nlm.nih.gov/taxonomy/ and entering the Taxonomy ID).
Delete the non-human proteins by clicking on the corresponding nodes and hitting delete.
2) Lets use a powerful feature of Cytoscape to colour the nodes according to their Gene
Ontology (GO) categories, the VizMapper. First, well need to retrieve the GO terms for the
proteins in the network and well do this by connecting to BioMart, which is a machine-readable
repository for many attributes associated with various bits of data stored at the EBI, the European
counterpart to the NCBI. Go to File>Import>Table>Public Databases. In the first Select
Services Import dialogue choose ENSEMBL GENES 78 (Sanger UK) and click OK. In the
second Import dialogue box, choose ENSEMBL GENES 78 (Sanger UK) Homo sapiens
(GRCh38), change the Key Column in Cytoscapedrop-down to shared name, with the Data
Type as EntrezGene ID(s), see Figure 6 for details. Select GO Term Accession/Definition/
Evidence Code/Name/Domain in the Import Settings list (and just the Definition further up the
list, for good measure). It is important that the identifier that were using to look up the
information at BioMart matches the identifier that BioGRID uses, which is what we make
happen when we change the Key Column in Cytoscape to shared name. Click Import and wait
a couple of minutes for the information to be retrieved. Click OK. Click Show All Columns
(
) to now see the Gene Ontology terms (you will need to scroll to the right to do so).
8
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Figure 6: Retrieving additional node information from BioMart. Here were retrieving Gene
Ontology terms. Click the Show All Columns button again to show the new data.
Now lets colour the nodes according to their Gene Ontology description. Click on the Style tab
in the Control Panel. Open the Fill Color to select options (click on the icon), then select GO
Term Definition as the attribute to use for the colouring. You will see all of the GO Term
Definitions listed. The Mapping Type should be set to Discrete Mapping. Click on Node Color
again, right-click to select Mapping Value Generators then select Rainbow. Voil! See Figure 7.
Figure 7: Using Fill Style to colour nodes according to discretized values (e.g. GO Definition).
9
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
b. Now the node colours correspond to the GO Term Definition categories. What are the general
gene ontology definitions for proteins that interact with BRCA2?
Although you can qualitatively see that there are a lot of
BRCA2-interacting proteins that are in the GO Term
Description categories similar to DNA damage
response, it would be nice to know if there were any
kind of over-representation relative to all gene in the
genome. Fortunately there are some Apps that will tell us
this exact thing! Go to Apps > App Manager and find
BiNGO (either by name or under the Ontology Analysis
tag). Click BiNGO and then Install it. It will take a few
minutes to download and install. When the download has
completed, close the App Manager window (see small
image to the right).
Select all of the nodes in the network by doing Select > Nodes > All Nodes (or just hold the left
mouse button while dragging the box to highlight all the nodes). Next, activate the BiNGO app
by clicking on BiNGO under Apps. Name the Cluster, and be sure to select Homo Sapiens (sic)
as the Organism/Annotation. Start BiNGO!
Figure 8: BiNGO Gene Ontology enrichment analysis with BRCA2-interacting proteins from
BioGRID. All nodes are selected then it is possible to run a Gene Ontology enrichment analysis
for GO Biological Process terms with the BiNGO app using Homo Sapiens as the
organism/annotation.
10
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
After a couple of minutes a table of GO BP terms appears along with their p-values for overrepresentation, along with a network which represents the GO term graph. Well ignore the graph
here and focus instead on the table (Figure 9). The smaller the p-value, the more significantly
enriched is the GO BP category.
Lab Quiz
Question 3
BRCA2 and RAD51 are known to form a critical interaction complex in the DNA damage
response. Disruption of this complex increases susceptibility to various forms of cancer.
d. Is RAD51 in the interaction network?(Hint: use the search bar along the top!)
Lets retrieve all the interactors from BioGRID for RAD51
and add them to the BRCA2 network. Right click on the
RAD51 node (which is labeled HsT16930, unless youve
changed the default labeling) and select Apps > Extend
Network by public interaction database and choose
BioGRID as the data source when prompted.
3) Remove duplicated edges by Edit > Remove Duplicated Edges and then selecting your
network when prompted. Check the Ignore Edge Direction option and click OK.
11
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
4) Organize the interaction graph by doing Layouts > Cytoscape Layouts > Edge-Weighted
Spring Embedded Layout (Biolayout) > All Nodes > (none). You can use the rotate function in
the tool panel to rotate the network to achieve a better fit on the screen.
5) Now lets reduce the graph to those nodes that interact with both BRCA2 and RAD51. First
do Tools > Network Analysis > Analyze Network to generate network statistics that we can filter
on. Treat the network as undirected and combine pair edges when prompted. The kinds of
network statistics that were generated were discussed in the mini-lecture and include node
degree, the number of edges emanating from each node. Then Click on the Select tab in the
Control Panel. Click on the + symbol to add a new filter. Set the values such that only the nodes
with 1 edge are highlighted (see Figure 10) and click on Apply Filter. Then do Edit > Delete
Selected Nodes and Edges.
Figure 10: Filtering a network based on network statistics (here, node degree) or other
parameters. Weve selected nodes having a node degree of 1 and will deleted them to explore the
BRCA2-RAD51 co-interactors further.
e. What do all the remaining nodes have in common?
Note: you may need to save this filtered network using File > Save, shut down and restart
Cytoscape, and then reload the network you generated to see the attributes of selected nodes in
the Data Panel. Even mature software is prone to bugs .
f. If you selected some interacting nodes, how would one use the list of ncbi_gene_id identifiers
in the Table Panel to search for potential interaction domains in these interacting proteins?
12
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
In this lab, weve seen that certain protein-protein interaction (PPI) databases and tools for
viewing PPIs have their strengths and weaknesses. For instance DIP seems to be fairly well
curated and contains links to the papers where the interactions were identified, but has a
somewhat clunky output where the nodes are named by internal identifiers. BioGRID offers a
great summary of the methods used to determine the interactions, but the ability to manipulate
the graphical representation of the network is limited. Finally, Cytoscape allows virtually
unlimited possibilities for the representation of a network in terms of layout options, node
appearance, etc. But, at least with the web service we used, the ability to identify how those
interactions were determined and to be able to access the primary literature concerning them is
limited.
End of Lab!
Where to get it:
Download the Cytoscape executable from http://www.cytoscape.org/download.html. Use the
Platform-Specific Installers to install Cytoscape 3.2.0. You will need to have the appropriate
Java Runtime Environment installed (32-bit or 64-bit for Windows users) first, which you can
get from http://www.java.com. During the set-up/start-up of Cytoscape, permit access with
private networks. Note: this lab has been tested and works with Cytoscape 3.2.0 on Windows,
Mac, and Linux machines.
Lab 2 Objectives
By the end of Lab 2 (comprising the labs including their boxes, and the lectures), you should:
understand why protein-protein interactions are important biologically, and also how they
may be determined experimentally;
be able to assess the advantages and disadvantages of the methods for determining
protein-protein interactions;
know the terminology associated with protein-protein interaction graphs;
be able to use DIP, BioGRID and Cytoscape to identify interacting proteins for your gene
product of interest and to filter and decorate networks based on additional information;
be able to identify the type of support for a given interaction in a given database;
be able to interpret the other types of information (GO categories) provided by the
software tools.
Do not hestitate to use the Coursera discussion forums if you do not understand any of the above
after reading the relevant material.
13
Copyright 2015 by D.S. Guttman and N.J. Provart
Bioinformatic Methods II
Lab 2
Further Reading
Blake JA (2013) Ten Quick Tips for Using the Gene Ontology. PLoS Comput Biol 9(11): e1003343.
doi:10.1371/journal.pcbi.1003343.
Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M,
Gross B, Hanspers K, Isserlin R, Kelley R, Killcoyne S, Lotia S, Maere S, Morris J, Ono K, Pavlovic V, Pico AR,
Vailaya A, Wang PL, Adler A, Conklin BR, Hood L, Kuiper M, Sander C, Schmulevich I, Schwikowski B, Warner
GJ, Ideker T, Bader GD (2007) Nat Protoc. 2(10):2366-82.
Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M (2006). BioGRID: A General Repository for
Interaction Datasets. Nucleic Acids Res. 34:D535-9.
Xenarios I, Rice DW, Salwinski L, Baron MK, Marcotte EM, Eisenberg D (2000). DIP: The Database of Interacting
Proteins. Nucl. Acids Res. 28:289-91.
14
Copyright 2015 by D.S. Guttman and N.J. Provart