Developing a logical model of yeast metabolism

With the completion of the sequencing of genomes of increasing numbers of organisms, the focus of biology is moving to determining the role of these genes (functional genomics). To this end it is useful to view the cell as a biochemical machine: it consumes simple molecules to manufacture more complex ones by chaining together biochemical reactions into long sequences referred to as em metabolic pathways. Such metabolic pathways are not linear but often interesect to form complex networks. Genes play a fundamental role in these networks by providing the information to synthesise the enzymes that catalyse biochemical reactions. Although developing a complete model of metabolism is of fundamental importance to biology and medicine, the size and complexity of the network has proven beyond the capacity of human reasoning. This paper presents the first results of the Robot Scientist research programme that aims to automatically discover the function of genes in the metabolism of the yeast em Saccharomyces cerevisiae. Results include: (1) the first logical model of metabolism;(2) a method to predict phenotype by deductive inference; and (3) a method to infer reactions and gene function by aductive inference. We describe the em in vivo experimental set-up which will allow these em in silico predictions to be automatically tested by a laboratory robot

Combining inductive logic programming, active learning and robotics to discover the function of genes

The paper is addressed to AI workers with an interest in biomolecular genetics and also to biomolecular geneticists interested in what AI tools may do for them. The authors are engaged in a collaborative enterprise aimed at partially automating some aspects of scientific work. These aspects include the processes of forming hypotheses, devising trials to discriminate between these competing hypotheses, physically performing these trials and then using the results of these trials to converge upon an accurate hypothesis. As a potential component of the reasoning carried out by an "artificial scientist" this paper describes ASE-Progol, an Active Learning system which uses Inductive Logic Programming to construct hypothesised first-order theories and uses a CART-like algorithm to select trials for eliminating ILP derived hypotheses. In simulated yeast growth tests ASE-Progol was used to rediscover how genes participate in the aromatic amino acid pathway of Saccharomyces cerevisiae. The cost of the chemicals consumed in converging upon a hypothesis with an accuracy of around 88% was reduced by five orders of magnitude when trials were selected by ASE-Progol rather than being sampled at random. While the naive strategy of always choosing the cheapest trial from the set of candidate trials led to lower cumulative costs than ASE-Progol, both the naive strategy and the random strategy took significantly longer to converge upon a final hypothesis than ASE-Progol. For example to reach an accuracy of 80%, ASE-Progol required 4 days while random sampling required 6 days and the naive strategy required 10 days

Local and global modes of drug action in biochemical networks

It becomes increasingly accepted that a shift is needed from the traditional target-based approach of drug development to an integrated perspective of drug action in biochemical systems. We here present an integrative analysis of the interactions between drugs and metabolism based on the concept of drug scope. The drug scope represents the set of metabolic compounds and reactions that are potentially affected by a drug. We constructed and analyzed the scopes of all US approved drugs having metabolic targets. Our analysis shows that the distribution of drug scopes is highly uneven, and that drugs can be classified into several categories based on their scopes. Some of them have small scopes corresponding to localized action, while others have large scopes corresponding to potential large-scale systemic action. These groups are well conserved throughout different topologies of the underlying metabolic network. They can furthermore be associated to specific drug therapeutic properties

The edges of understanding

A culture's icons are a window onto its soul. Few would disagree that, in the culture of molecular biology that dominated much of the life sciences for the last third of the 20th century, the dominant icon was the double helix. In the present, post-modern, 'systems biology' era, however, it is, arguably, the hairball

Scientific discovery as a combinatorial optimisation problem: How best to navigate the landscape of possible experiments?

A considerable number of areas of bioscience, including gene and drug discovery, metabolic engineering for the biotechnological improvement of organisms, and the processes of natural and directed evolution, are best viewed in terms of a ‘landscape’ representing a large search space of possible solutions or experiments populated by a considerably smaller number of actual solutions that then emerge. This is what makes these problems ‘hard’, but as such these are to be seen as combinatorial optimisation problems that are best attacked by heuristic methods known from that field. Such landscapes, which may also represent or include multiple objectives, are effectively modelled in silico, with modern active learning algorithms such as those based on Darwinian evolution providing guidance, using existing knowledge, as to what is the ‘best’ experiment to do next. An awareness, and the application, of these methods can thereby enhance the scientific discovery process considerably. This analysis fits comfortably with an emerging epistemology that sees scientific reasoning, the search for solutions, and scientific discovery as Bayesian processes

Edible crabs “Go West”: migrations and incubation cycle of Cancer pagurus revealed by electronic tags

Crustaceans are key components of marine ecosystems which, like other exploited marine taxa, show seasonable patterns of distribution and activity, with consequences for their availability to capture by targeted fisheries. Despite concerns over the sustainability of crab fisheries worldwide, difficulties in observing crabs’ behaviour over their annual cycles, and the timings and durations of reproduction, remain poorly understood. From the release of 128 mature female edible crabs tagged with electronic data storage tags (DSTs), we demonstrate predominantly westward migration in the English Channel. Eastern Channel crabs migrated further than western Channel crabs, while crabs released outside the Channel showed little or no migration. Individual migrations were punctuated by a 7-month hiatus, when crabs remained stationary, coincident with the main period of crab spawning and egg incubation. Incubation commenced earlier in the west, from late October onwards, and brooding locations, determined using tidal geolocation, occurred throughout the species range. With an overall return rate of 34%, our results demonstrate that previous reluctance to tag crabs with relatively high-cost DSTs for fear of loss following moulting is unfounded, and that DSTs can generate precise information with regards life-history metrics that would be unachievable using other conventional means

Efficient discovery of anti-inflammatory small-molecule combinations using evolutionary computing

The control of biochemical fluxes is distributed, and to perturb complex intracellular networks effectively it is often necessary to modulate several steps simultaneously. However, the number of possible permutations leads to a combinatorial explosion in the number of experiments that would have to be performed in a complete analysis. We used a multiobjective evolutionary algorithm to optimize reagent combinations from a dynamic chemical library of 33 compounds with established or predicted targets in the regulatory network controlling IL-1β expression. The evolutionary algorithm converged on excellent solutions within 11 generations, during which we studied just 550 combinations out of the potential search space of ~9 billion. The top five reagents with the greatest contribution to combinatorial effects throughout the evolutionary algorithm were then optimized pairwise. A p38 MAPK inhibitor together with either an inhibitor of IκB kinase or a chelator of poorly liganded iron yielded synergistic inhibition of macrophage IL-1β expression. Evolutionary searches provide a powerful and general approach to the discovery of new combinations of pharmacological agents with therapeutic indices potentially greater than those of single drugs

Mining metabolites: extracting the yeast metabolome from the literature

Text mining methods have added considerably to our capacity to extract biological knowledge from the literature. Recently the field of systems biology has begun to model and simulate metabolic networks, requiring knowledge of the set of molecules involved. While genomics and proteomics technologies are able to supply the macromolecular parts list, the metabolites are less easily assembled. Most metabolites are known and reported through the scientific literature, rather than through large-scale experimental surveys. Thus it is important to recover them from the literature. Here we present a novel tool to automatically identify metabolite names in the literature, and associate structures where possible, to define the reported yeast metabolome. With ten-fold cross validation on a manually annotated corpus, our recognition tool generates an f-score of 78.49 (precision of 83.02) and demonstrates greater suitability in identifying metabolite names than other existing recognition tools for general chemical molecules. The metabolite recognition tool has been applied to the literature covering an important model organism, the yeast Saccharomyces cerevisiae, to define its reported metabolome. By coupling to ChemSpider, a major chemical database, we have identified structures for much of the reported metabolome and, where structure identification fails, been able to suggest extensions to ChemSpider. Our manually annotated gold-standard data on 296 abstracts are available as supplementary materials. Metabolite names and, where appropriate, structures are also available as supplementary materials