Go to the About Us section
Go to the Courses section
Go to the Research section
Go to the Specialist Services & Consultancy section
Go to the Outreach section
Go to the New Students section
You can find details about all the staff in the department below.
return to full list
Department of Biological and Medical Sciences
Faculty of Health and Life Sciences
After studying biochemistry at Oxford University, followed by a DPhil on the physical biochemistry of yeast pyruvate kinase, David started lecturing at Oxford Polytechnic. His research gradually moved from experimental biochemistry into computer simulation and theoretical analysis of metabolic control, and he has written the only textbook on metabolic control analysis, Understanding the Control of Metabolism. In 2001, he helped to found and became part-time Chief Scientific Officer of the Oxford company, Physiomics plc, which is using computer simulation of cellular systems for the development and analysis of therapeutic strategies for the pharmaceutical industry.
Physiomics has firmly established itself as a leading light in systems biology approaches to drug discovery and latterly in therapy design, demonstrable through contracts with three major international pharmaceutical companies. Through its strong advocacy of this approach the sector has invested in and adopted new computational biology processes. As Physiomics has continued to grow, it has expanded its own specialist research team, in many cases recruiting scientists trained within Fell’s Brookes-based research group at Brookes.
David is chairman of the Policy Committee of the Biochemical Society, and has been a member of several panels and committees of the Biotechnology and Biological Sciences Research Council.
David’s group formed nearly thirty years ago with initial interests in computer simulation of metabolism and the theory of metabolic control. To these it has since added interests in modelling signal transduction, in various different approaches to network analysis of metabolism, and in reconstructing metabolic networks from genomic data. In the course of this research, he has addressed problems in microbial, plant and mammalian metabolism, often in conjunction with collaborators who have contributed experimental results.
His work forms part of the emerging field of Systems Biology, in that we are concerned with understanding how biological function arises from the interactions between many components, and with building predictive models. Potential applications of our work include the design of changes in cellular metabolism to improve the output of product such as antibiotics, detecting vulnerable sites in cellular networks that could be targets for drugs to control disease-causing organisms, and improved understanding of how organisms manage to adjust their metabolism in response to environmental changes and other signals.
Analysis of the impact of photorespiration on plant metabolism is usually based on manual inspection of small network diagrams. Here we create a structural metabolic model that contains the reactions that participate in photorespiration in the plastid, peroxisome, mitochondrion and cytosol and the metabolite exchanges between them. This model was subjected to elementary flux modes analysis, a technique that enumerates all the component, minimal pathways of a network. Any feasible photorespiratory metabolism in the plant will be some combination of the elementary flux modes (EFMs) that contain the Rubisco oxygenase reaction. Amongst the EFMs we obtained was the classic photorespiratory cycle, but there were also modes that involve photorespiration coupled with mitochondrial metabolism and ATP production, the glutathione‐ascorbate (GSH‐ASC) cycle and nitrate reduction to ammonia. The modes analysis demonstrated the underlying basis of the metabolic linkages with photorespiration that have been inferred experimentally. The set of reactions common to all the elementary modes showed good agreement with the gene products of mutants that have been reported to have a defective phenotype in photorespiratory conditions. Finally, the set of modes provided a formal demonstration that photorespiration itself does not impact on the CO2:O2 ratio (assimilation quotient, AQ), except in those modes associated with concomitant nitrate reduction.
Previously we have used a genome scale model of rice metabolism to describe how metabolism reconfigures at different light intensities in an expanding leaf of rice. Although this established that the metabolism of the leaf was adequately represented, in the model, the scenario was not that of the typical function of the leaf—to provide material for the rest of the plant. Here we extend our analysis to explore the transition to a source leaf as export of photosynthate increases at the expense of making leaf biomass precursors, again as a function of light intensity. In particular we investigate whether, when the leaf is making a smaller range of compounds for export to the phloem, the same changes occur in the interactions between mitochondrial and chloroplast metabolism as seen in biomass synthesis for growth when light intensity increases. Our results show that the same changes occur qualitatively, though there are slight quantitative differences reflecting differences in the energy and redox requirements for the different metabolic outputs.
Produce rich in phytochemicals may alter postprandial glucose and insulin responses by interacting with the pathways that regulate glucose uptake and insulin secretion in humans. The aims of the present study were to assess the phytochemical constituents of red beetroot juice and to measure the postprandial glucose and insulin responses elicited by either 225 ml beetroot juice (BEET), a control beverage matched for macronutrient content (MCON) or a glucose beverage in healthy adults. Beetroot juice was a particularly rich source of betalain degradation compounds. The orange/yellow pigment neobetanin was measured in particularly high quantities (providing 1·3 g in the 225 ml). A total of sixteen healthy individuals were recruited, and consumed the test meals in a controlled single-blind cross-over design. Results revealed a significant lowering of the postprandial insulin response in the early phase (0–60 min) (P < 0·05) and a significantly lower glucose response in the 0–30 min phase (P < 0·05) in the BEET treatment compared with MCON. Betalains, polyphenols and dietary nitrate found in the beetroot juice may each contribute to the observed differences in the postprandial insulin concentration.
We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.
Flux balance models of metabolism generally utilize synthesis of biomass as the main determinant of intracellular fluxes. However, the biomass constraint alone is not sufficient to predict realistic fluxes in central heterotrophic metabolism of plant cells because of the major demand on the energy budget due to transport costs and cell maintenance. This major limitation can be addressed by incorporating transport steps into the metabolic model and by implementing a procedure that uses Pareto optimality analysis to explore the trade-off between ATP and NADPH production for maintenance. This leads to a method for predicting cell maintenance costs on the basis of the measured flux ratio between the oxidative steps of the oxidative pentose phosphate pathway and glycolysis. We show that accounting for transport and maintenance costs substantially improves the accuracy of fluxes predicted from a flux balance model of heterotrophic Arabidopsis cells in culture, irrespective of the objective function used in the analysis. Moreover, when the new method was applied to cells under control, elevated temperature and hyper-osmotic conditions, only elevated temperature led to a substantial increase in cell maintenance costs. It is concluded that the hyper-osmotic conditions tested did not impose a metabolic stress, in as much as the metabolic network is not forced to devote more resources to cell maintenance.
The active, yet energetically inefficient electron transport chain of the ethanologenic bacterium Zymomonas mobilis could be used in metabolic engineering for redox-balancing purposes during synthesis of certain products. Although several reconstructions of Z. mobilis metabolism have been published, important aspects of redox balance and aerobic catabolism have not previously been considered. Here, annotated genome sequences and metabolic reconstructions have been combined with existing biochemical evidence to yield a medium-scale model of Z. mobilis central metabolism in the form of COBRA Toolbox model files for flux balance analysis (FBA). The stoichiometric analysis presented here suggests the feasibility of several metabolic engineering strategies for obtaining high-value products, such as glycerate, succinate, and glutamate that would use the electron transport chain to oxidize the excess NAD(P)H, generated during synthesis of these metabolites. Oxidation of the excess NAD(P)H would also be needed for synthesis of ethanol from glycerol. Maximum product yields and the byproduct spectra have been estimated for each product, with glucose, xylose, or glycerol as the carbon substrates. These novel pathways represent targets for future metabolic engineering, as they would exploit both the rapid Entner–Doudoroff glycolysis, and the energetically uncoupled electron transport of Z. mobilis.
Organisms share a common core to their metabolic networks. But what determined this: chance, chemical necessity, or evolutionary optimization? In this issue of Molecular Cell, Noor et al. (2010) provide new evidence for selection of a network with optimal features from a broader set of possibilities.
Reconstructing a model of the metabolic network of an organism from its annotated genome sequence would seem, at first sight, to be one of the most straightforward tasks in functional genomics, even if the various data sources required were never designed with this application in mind. The number of genome-scale metabolic models is, however, lagging far behind the number of sequenced genomes and is likely to continue to do so unless the model-building process can be accelerated. Two aspects that could usefully be improved are the ability to find the sources of error in a nascent model rapidly, and the generation of tenable hypotheses concerning solutions that would improve a model. We will illustrate these issues with approaches we have developed in the course of building metabolic models of Streptococcus agalactiae and Arabidopsis thaliana.
Flux is a key measure of the metabolic phenotype. Recently, complete (genome-scale) metabolic network models have been established for Arabidopsis (Arabidopsis thaliana), and flux distributions have been predicted using constraints-based modeling and optimization algorithms such as linear programming. While these models are useful for investigating possible flux states under different metabolic scenarios, it is not clear how close the predicted flux distributions are to those occurring in vivo. To address this, fluxes were predicted for heterotrophic Arabidopsis cells and compared with fluxes estimated in parallel by 13C-metabolic flux analysis (MFA). Reactions of the central carbon metabolic network (glycolysis, the oxidative pentose phosphate pathway, and the tricarboxylic acid [TCA] cycle) were independently analyzed by the two approaches. Net fluxes in glycolysis and the TCA cycle were predicted accurately from the genome-scale model, whereas the oxidative pentose phosphate pathway was poorly predicted. MFA showed that increased temperature and hyperosmotic stress, which altered cell growth, also affected the intracellular flux distribution. Under both conditions, the genome-scale model was able to predict both the direction and magnitude of the changes in flux: namely, increased TCA cycle and decreased phosphoenolpyruvate carboxylase flux at high temperature and a general decrease in fluxes under hyperosmotic stress. MFA also revealed a 3-fold reduction in carbon-use efficiency at the higher temperature. It is concluded that constraints-based genome-scale modeling can be used to predict flux changes in central carbon metabolism under stress conditions.
Motivation: In recent years, several methods have been proposed for determining metabolic pathways in an automated way based on network topology. The aim of this work is to analyse these methods by tackling a concrete example relevant in biochemistry. It concerns the question whether even-chain fatty acids, being the most important constituents of lipids, can be converted into sugars at steady state. It was proved five decades ago that this conversion using the Krebs cycle is impossible unless the enzymes of the glyoxylate shunt (or alternative bypasses) are present in the system. Using this example, we can compare the various methods in pathway analysis. Results: Elementary modes analysis (EMA) of a set of enzymes corresponding to the Krebs cycle, glycolysis and gluconeogenesis supports the scientific evidence showing that there is no pathway capable of converting acetyl-CoA to glucose at steady state. This conversion is possible after the addition of isocitrate lyase and malate synthase (forming the glyoxylate shunt) to the system. Dealing with the same example, we compare EMA with two tools based on graph theory available online, PathFinding and Pathway Hunter Tool. These automated network generating tools do not succeed in predicting the conversions known from experiment. They sometimes generate unbalanced paths and reveal problems identifying side metabolites that are not responsible for the carbon net flux. This shows that, for metabolic pathway analysis, it is important to consider the topology (including bimolecular reactions) and stoichiometry of metabolic systems, as is done in EMA.
A decrease in retinoic acid levels due to alcohol consumption has been proposed as a contributor to such conditions as fetal alcohol spectrum diseases and ethanol-induced cancers. One molecular mechanism, competitive inhibition by ethanol of the catalytic activity of human alcohol dehydrogenase (EC 18.104.22.168) (ADH) on all-trans-retinol oxidation has been shown for the ADH7 isoform. Ethanol metabolism also causes an increase in the free reduced nicotinamide adenine dinucleotide (NADH) in cells, which might reasonably be expected to decrease the retinol oxidation rate by product inhibition of ADH isoforms. To understand the relative importance o