SAQP - SRI International
Download
Report
Transcript SAQP - SRI International
The Structured Advanced Query
Page
Tomer Altman
Mario Latendresse
Bioinformatics Research Group
SRI, International
May 17, 2010
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 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
1
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.
1
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.
1
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).
1
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)
1
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.
1
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.
1
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.
1
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.
1
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.
1
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.
1
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.
1
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
1
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]
1
SRI International Bioinformatics