150 likes | 258 Views
The Structured Advanced Query Page. Tomer Altman Mario Latendresse Bioinformatics Research Group SRI International June, 2013. SRI International Bioinformatics. 1. Introduction. BioVelo is a query language Like SQL but simpler and no updates allowed
E N D
The Structured Advanced Query Page Tomer Altman Mario Latendresse Bioinformatics Research Group SRI International June, 2013 SRI International Bioinformatics 1
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 SRI International Bioinformatics 2
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. SRI International Bioinformatics 3
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. SRI International Bioinformatics 4
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). SRI International Bioinformatics 5
The Structure of the SAQP: • Database specification • Class specification • 'Where' constraints on attributes of classes • Output attributes description • Data format (HTML vs TXT) SRI International Bioinformatics 6
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. SRI International Bioinformatics 7
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 tab-delimited text file is generated, with just the table data. SRI International Bioinformatics 8
Example #2: • Find all the proteins of E. coli K-12 for which the DNA-FOOTPRINT-SIZE is smaller than 10. • Display the protein name, and the DNA footprint size. SRI International Bioinformatics 9
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. SRI International Bioinformatics 10
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. SRI International Bioinformatics 11
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 K-12 using the common-name attribute. SRI International Bioinformatics 12
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. SRI International Bioinformatics 13
Introduction to BioVelo • BioVelo is based on set and listcomprehension. • 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 SRI International Bioinformatics 14
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-end-position – g^right-end-position)+1, k < 200 ] • [(r^?name, c^?name) : r<- ecoli^^reactions, c<- r^left, c in r^right] SRI International Bioinformatics 15