In silico biotechnology. Era of reconstruction and interrogation

Functional characterization of alternate optimal solutions of Escherichia coli's transcriptional and translational machinery

The constraint-based reconstruction and analysis approach has recently been extended to describe Escherichia coli's transcriptional and translational machinery. Here, we introduce the concept of reaction coupling to represent the dependency between protein synthesis and utilization. These coupling constraints lead to a significant contraction of the feasible set of steady-state fluxes. The subset of alternate optimal solutions (AOS) consistent with maximal ribosome production was calculated. The majority of transcriptional and translational reactions were active for all of these AOS, showing that the network has a low degree of redundancy. Furthermore, all calculated AOS contained the qualitative expression of at least 92% of the known essential genes. Principal component analysis of AOS demonstrated that energy currencies (ATP, GTP, and phosphate) dominate the network's capability to produce ribosomes. Additionally, we identified regulatory control points of the network, which include the transcription reactions of sigma70 (RpoD) as well as that of a degradosome component (Rne) and of tRNA charging (ValS). These reactions contribute significant variance among AOS. These results show that constraint-based modeling can be applied to gain insight into the systemic properties of E. coli's transcriptional and translational machinery.

Filling kinetic gaps: dynamic modeling of metabolism where detailed kinetic information is lacking

BACKGROUND

Integrative analysis between dynamical modeling of metabolic networks and data obtained from high throughput technology represents a worthy effort toward a holistic understanding of the link among phenotype and dynamical response. Even though the theoretical foundation for modeling metabolic network has been extensively treated elsewhere, the lack of kinetic information has limited the analysis in most of the cases. To overcome this constraint, we present and illustrate a new statistical approach that has two purposes: integrate high throughput data and survey the general dynamical mechanisms emerging for a slightly perturbed metabolic network.

METHODOLOGY/PRINCIPAL FINDINGS

This paper presents a statistic framework capable to study how and how fast the metabolites participating in a perturbed metabolic network reach a steady-state. Instead of requiring accurate kinetic information, this approach uses high throughput metabolome technology to define a feasible kinetic library, which constitutes the base for identifying, statistical and dynamical properties during the relaxation. For the sake of illustration we have applied this approach to the human Red blood cell metabolism (hRBC) and its capacity to predict temporal phenomena was evaluated. Remarkable, the main dynamical properties obtained from a detailed kinetic model in hRBC were recovered by our statistical approach. Furthermore, robust properties in time scale and metabolite organization were identify and one concluded that they are a consequence of the combined performance of redundancies and variability in metabolite participation.

CONCLUSIONS/SIGNIFICANCE

In this work we present an approach that integrates high throughput metabolome data to define the dynamic behavior of a slightly perturbed metabolic network where kinetic information is lacking. Having information of metabolite concentrations at steady-state, this method has significant relevance due its potential scope to analyze others genome scale metabolic reconstructions. Thus, I expect this approach will significantly contribute to explore the relationship between dynamic and physiology in other metabolic reconstructions, particularly those whose kinetic information is practically nulls. For instances, I envisage that this approach can be useful in genomic medicine or pharmacogenomics, where the estimation of time scales and the identification of metabolite organization may be crucial to characterize and identify (dis)functional stages.

A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory

Background

Pseudomonas putida is the best studied pollutant degradative bacteria and is harnessed by industrial biotechnology to synthesize fine chemicals. Since the publication of P. putida KT2440's genome, some in silico analyses of its metabolic and biotechnology capacities have been published. However, global understanding of the capabilities of P. putida KT2440 requires the construction of a metabolic model that enables the integration of classical experimental data along with genomic and high-throughput data. The constraint-based reconstruction and analysis (COBRA) approach has been successfully used to build and analyze in silico genome-scale metabolic reconstructions.

Results

We present a genome-scale reconstruction of P. putida KT2440's metabolism, iJN746, which was constructed based on genomic, biochemical, and physiological information. This manually-curated reconstruction accounts for 746 genes, 950 reactions, and 911 metabolites. iJN746 captures biotechnologically relevant pathways, including polyhydroxyalkanoate synthesis and catabolic pathways of aromatic compounds (e.g., toluene, benzoate, phenylacetate, nicotinate), not described in other metabolic reconstructions or biochemical databases. The predictive potential of iJN746 was validated using experimental data including growth performance and gene deletion studies. Furthermore, in silico growth on toluene was found to be oxygen-limited, suggesting the existence of oxygen-efficient pathways not yet annotated in P. putida's genome. Moreover, we evaluated the production efficiency of polyhydroxyalkanoates from various carbon sources and found fatty acids as the most prominent candidates, as expected.

Conclusion

Here we presented the first genome-scale reconstruction of P. putida, a biotechnologically interesting all-surrounder. Taken together, this work illustrates the utility of iJN746 as i) a knowledge-base, ii) a discovery tool, and iii) an engineering platform to explore P. putida's potential in bioremediation and bioplastic production.

Metabolic reconstruction and modeling of nitrogen fixation in Rhizobium etli

Rhizobiaceas are bacteria that fix nitrogen during symbiosis with plants. This symbiotic relationship is crucial for the nitrogen cycle, and understanding symbiotic mechanisms is a scientific challenge with direct applications in agronomy and plant development. Rhizobium etli is a bacteria which provides legumes with ammonia (among other chemical compounds), thereby stimulating plant growth. A genome-scale approach, integrating the biochemical information available for R. etli, constitutes an important step toward understanding the symbiotic relationship and its possible improvement. In this work we present a genome-scale metabolic reconstruction (iOR363) for R. etli CFN42, which includes 387 metabolic and transport reactions across 26 metabolic pathways. This model was used to analyze the physiological capabilities of R. etli during stages of nitrogen fixation. To study the physiological capacities in silico, an objective function was formulated to simulate symbiotic nitrogen fixation. Flux balance analysis (FBA) was performed, and the predicted active metabolic pathways agreed qualitatively with experimental observations. In addition, predictions for the effects of gene deletions during nitrogen fixation in Rhizobia in silico also agreed with reported experimental data. Overall, we present some evidence supporting that FBA of the reconstructed metabolic network for R. etli provides results that are in agreement with physiological observations. Thus, as for other organisms, the reconstructed genome-scale metabolic network provides an important framework which allows us to compare model predictions with experimental measurements and eventually generate hypotheses on ways to improve nitrogen fixation.

Understanding the roadmap of metabolism by pathway analysis

The theoretical investigation of the structure of metabolic systems has recently attracted increasing interest. In this chapter, the basic concepts of metabolic pathway analysis are described and various applications are outlined. In particular, the concepts of nullspace and elementary flux modes are explained. The presentation is illustrated by a simple example from tyrosine metabolism and a system describing lysine production in Corynebacterium glutamicum. The latter system gives rise to 37 elementary modes, 36 of which produce lysine with different molar yields. The examples illustrate that metabolic pathway analysis is a useful tool for better understanding the complex architecture of intracellular metabolism, for determining the pathways on which the molar conversion yield of a substrate-product pair under study is maximal, and for assigning functions to orphan genes (functional genomics). Moreover, problems emerging in the modeling of large networks are discussed. An outlook on current trends in the field concludes the chapter.

Application of global sensitivity analysis to determine goals for design of experiments: an example study on antibody-producing cell cultures

Global sensitivity analysis (GSA) can be used to quantify the importance of model parameters and their interactions with respect to model output. In this study, the Sobol' method for GSA is applied to a dynamic model of monoclonal antibody-producing mammalian cell cultures in order to identify the parameters that need to be accurately determined experimentally. Our results show that most parameters have low sensitivity indices and exhibit strong interactions with one another. These parameters can be set at their nominal values and unnecessary experimentation can therefore be avoided. In contrast, certain parameters are identified as sensitive, necessitating their estimation given sufficiently rich experimental data. Moreover, parameter sensitivity varies during culture time in a biologically meaningful manner. In conclusion, GSA can serve as an excellent precursor to optimal experiment design.

In silico reconstruction of the metabolic pathways of Lactobacillus plantarum: comparing predictions of nutrient requirements with those from growth experiments

On the basis of the annotated genome we reconstructed the metabolic pathways of the lactic acid bacterium Lactobacillus plantarum WCFS1. After automatic reconstruction by the Pathologic tool of Pathway Tools (http://bioinformatics.ai.sri.com/ptools/), the resulting pathway-genome database, LacplantCyc, was manually curated extensively. The current database contains refinements to existing routes and new gram-positive bacterium-specific reactions that were not present in the MetaCyc database. These reactions include, for example, reactions related to cell wall biosynthesis, molybdopterin biosynthesis, and transport. At present, LacplantCyc includes 129 pathways and 704 predicted reactions involving some 670 chemical species and 710 enzymes. We tested vitamin and amino acid requirements of L. plantarum experimentally and compared the results with the pathways present in LacplantCyc. In the majority of cases (32 of 37 cases) the experimental results agreed with the final reconstruction. LacplantCyc is the most extensively curated pathway-genome database for gram-positive bacteria and is open to the microbiology community via the World Wide Web (www.lacplantcyc.nl). It can be used as a reference pathway-genome database for gram-positive microbes in general and lactic acid bacteria in particular.

Metabolic pathways reconstruction by frequency and amplitude response to forced glycolytic oscillations in yeast

The hypothesis that frequency and amplitude response can be used in a complicated metabolic pathway kinetics model for optimal parameter estimation, as speculated by its successful prior usage for a mechanical oscillator and a heterogeneous chemical system, is tested here. Given the complexity of the glycolysis model of yeast chosen, this question is limited to three kinetics parameters of the 87 in the in vitro model developed in the literature. The direct application of the approach, used with the uninformed selection of operating conditions for the oscillation of external glucose concentration, led to miring the data assimilation process in local minima. Application of linear systems theory, however, identified two natural resonant frequencies that, when excited by external forced oscillations of the same frequency, result in the expression of many harmonics in the Fourier spectra, that is, information-rich experiments. A single such information-rich experiment at one of the resonant frequencies was sufficient to break away from the local minima to find the optimum kinetics parameter estimates. The resonant frequencies themselves represent oscillation modes in glycolysis akin to those previously observed. Furthermore, operation of the bioreactor with large amplitude oscillations of glucose feed (25%) leads to enhanced ethanol average yield by 1.6% at the resonant frequency.

Reconstructing the metabolic network of a bacterium from its genome

The prospect of understanding the relationship between the genome and the physiology of an organism is an important incentive to reconstruct metabolic networks. The first steps in the process can be automated and it does not take much effort to obtain an initial metabolic reconstruction from a genome sequence. However, such a reconstruction is certainly not flawless and correction of the many imperfections is laborious. It requires the combined analysis of the available information on protein sequence, phylogeny, gene-context and co-occurrence but is also aided by high-throughput experimental data. Simultaneously, the reconstructed network provides the opportunity to visualize the "omics" data within a relevant biological functional context and thus aids the interpretation of those data.

Microarray reality checks in the context of a complex disease

A problem in analyzing microarray-based gene expression data is the separation of genes causally involved in a disease from innocent bystander genes, whose expression levels have been secondarily altered by primary changes elsewhere. To investigate this issue systematically in the context of a class of complex human diseases, we have compared microarray-based gene expression data with non-microarray-based clinical and biological data about the schizophrenias to ask whether these two approaches prioritize the same genes. We find that genes whose expression changes are deemed to be of importance from microarrays are rarely those classified as of importance from clinical, in situ, molecular, single-nucleotide polymorphism (SNP) association, knockout and drug perturbation data. This disparity is not limited to the schizophrenias but characterizes other human disease data sets. It also extends to biological validation of microarray data in model organisms, in which genome-wide phenotypic data have been systematically compared with microarray data. In addition, different bioinformatic protocols applied to the same microarray data yield quite different gene sets and thus make clinical decisions less straightforward. We discuss how progress may be improved in the clinical area by the assignment of high-quality phenotypic values to each member of a microarray-assigned gene set.