Download The Structured Advanced Query Page

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts
no text concepts found
Transcript
The Structured Advanced Query
Page
Tomer Altman
Mario Latendresse
Bioinformatics Research Group
SRI International
June, 2013
1
SRI International Bioinformatics
Introduction
BioVelo
is a query language
Like SQL but simpler and no updates allowed
Documentation: http://biocyc.org/bioveloLanguage.shtml
Free-Form Advanced Query Page (FFAQP) allows Web
submission of BioVelo queries
Structured
Advanced Query Page (SAQP)
Web page for interactively constructing advanced and precise
queries to PGDBs
Queries are translated to BioVelo and sent to the server for
processing
SAQP: http://biocyc.org/query.shtml
Documentation: http://biocyc.org/webQueryDoc.shtml
2
SRI International Bioinformatics
Why a query interface?
 Allow
a structured way to access the rich data
representation stored in a PGDB.
 Most advanced databases have a high-level,
declarative method of access (i.e., SQL).
 Provides an intermediate level of access between
graphically browsing the PGDB and programmatically
processing the data using Lisp.
3
SRI International Bioinformatics
The Structured Advanced Query
Page
 'Advanced',
in that it allows you to ask more advanced
and complicated queries than the basic search
interface.
 In other words, the SAQP allows you to search for a
precise set of answers given simple or complex
conditions
 'Structured', in that it is a dynamic HTML form, that
provides greater ease in crafting queries, but trades
flexibility and power for simplicity (FFAQP).
 'Page', in that it is accessed via the Web interface for
BioCyc (www.biocyc.org/query.shtml), or from your
own Pathway Tools Web server.
4
SRI International Bioinformatics
SAQP Architecture
 The
SAQP is built on top of a high-level functional
declarative language called BioVelo (Mario
Latendresse, SRI), which is built on top of Pathway
Tools.
 On every result page, you will see the equivalent
BioVelo code that was generated from the SAQP,
which, in turn, generated the results.
 You don't need to know anything about BioVelo to use
the SAQP, but it might be helpful later if you need the
ability to write even more complicated queries using
the Free Form Advanced Query Page (FFAQP).
5
SRI International Bioinformatics
The Structure of the SAQP:
 Database
specification
 Class specification
 'Where' constraints on attributes of classes
 Output attributes description
 Data format (HTML vs TXT)
6
SRI International Bioinformatics
Example #1:
A
simple query usually consists of querying a
particular database about a particular class.
 Find all the proteins in E. coli K-12.
 Display the protein names.
7
SRI International Bioinformatics
Structure of the Results
A
line that shows the equivalent BioVelo expression
that the SAQP generated to answer the query.
 A HTML table of the results, with the corresponding
entries hyperlinked to the matching Pathway Tools
Web pages.
 If a text data format was requested, then a tabdelimited text file is generated, with just the table data.
8
SRI International Bioinformatics
Example #2:
 Find
all the proteins of E. coli K-12 for which the DNAFOOTPRINT-SIZE is smaller than 10.
 Display the protein name, and the DNA footprint size.
9
SRI International Bioinformatics
Example #3:
 In
EcoCyc, display polypeptides constrained by
experimentally determined molecular weight and
isoelectric point.
 The experimental molecular weight should be between
50 and 100 kD.
 The pI should be less than 7.
 Display the polypeptide name, the experimental
molecular weight, and the pI.
10
SRI International Bioinformatics
Example #4:
 The
SAQP allows for specifying quantifiers on
relations between PGDB classes.
 Extending example #3, now we want only proteins
where at least one of the genes that encodes the
protein to be within the first 500 kilobases of the E. coli
chromosome.
11
SRI International Bioinformatics
Example #5: Queries with
Several Components
A
second search component will search potentially
another database and another class of objects for each
element found in the first search component.
 It is called a 'cross-product' search.
 Any number of search components can be added. In
general, the new search component is done for each
set of objects found in the previous components.
 Some restraints is needed not to build a query that
takes too long to answer. (The server gives a limit of a
few minutes for a query.)
 Example: Search for MetaCyc pathways in the
taxonomic range of Bacteria that also exist in E. coli K12 using the common-name attribute.
12
SRI International Bioinformatics
Exercises
1) Find all genes of E. coli that contain “trp” in their
name.
2) Find all genes in MetaCyc that have more than one
product. Output the gene names and product
names.
3) Find all reactions in E. coli which have the reactant
(i.e., the left side) “acetaldehyde”.
4) Find all monomers in E. coli. A monomer has no
components.
5) Find all reactions in MetaCyc that have more than 4
reactants.
6) Find all metabolic pathways, in MetaCyc, that have
more than 5 reactions. Output the reaction lists as
well as the pathway names.
13
SRI International Bioinformatics
Introduction to BioVelo
 BioVelo
is based on set and list comprehension.
 In Mathematics, a set comprehension describes a set of
values as in: {x | x in Prime, x > 100}
 The output is 'x', the body has a generator 'x in Prime' and
a condition 'x > 100'. Several conditions and several
generators could be used.
 BioVelo used a concise syntax:
1) [ output-expression : generator, condition, ... ]
2) a generator has the form v ← database^^class
3) a condition uses logical and relational operators
14
SRI International Bioinformatics
Examples of BioVelo Queries
 [r : r <- ecoli^^reactions]
 [p^name : p <- ecoli^^proteins]
 [p^?name : p<- ecoli^^proteins]
 [p^?name : p <- ecoli^^proteins, p^dna-footprint-size < 10]
 [(g^?name, g^left-end-position): g <- ecoli^^genes,
g^left-end-position < 153000]
 [(g^?name, k): g<- ecoli^^genes, k := abs(g^left-endposition – g^right-end-position)+1, k < 200 ]
 [(r^?name, c^?name) : r<- ecoli^^reactions, c<- r^left, c in
r^right]
15
SRI International Bioinformatics