A roadmap for interpreting 13C metabolite labeling patterns from cells

Measuring intracellular metabolism has increasingly led to important insights in biomedical research. [superscript 13]C tracer analysis, although less information-rich than quantitative [superscript 13]C flux analysis that requires computational data integration, has been established as a time-efficient method to unravel relative pathway activities, qualitative changes in pathway contributions, and nutrientcontributions. Here, we review selected key issues in interpreting [superscript 13]C metabolite labeling patterns, with the goal of drawing accurate conclusions from steady state and dynamic stable isotopic tracer experiment

Bioinformatics and Handwriting/Speech Reconition: Uncoventional Applications of Similarity Search Tools

This work introduces two unconventional applications for sequence alignment algorithms outside the domain of bioinformatics: handwriting recognition and speech recognition. In each application we treated data samples, such as the path of a and written pen stroke, as a protein sequence and use the FastA sequence alignment tool to classify unknown data samples, such as a written character. That is, we handle the handwriting and speech recognition problems like the protein annotation problem: given a sequence of unknown function, we annotate the sequence via sequence alignment. This approach achieves classification rates of 99.65% and 93.84% for the handwriting and speech recognition respectively. In addition, we provide a framework for applying sequence alignment to a variety of other non–traditional problems.Singapore-MIT Alliance (SMA

Activation Ratios For Reconstruction Of Signal Transduction Networks

We have developed a novel framework that can be applied for the analysis of signal transduction networks, both to facilitate reconstruction of the network structure and quantitatively characterize the interaction between network components. This approach, termed activation ratio analysis, involves the ratio between active and inactive forms of signaling intermediates at steady state. The activation ratio of an intermediate is shown to depend linearly upon the concentration of the activating enzyme. The slope of the line is defined as the activation factor, and is determined by the kinetic parameters of activation and inactivation. When activation ratios for simple signaling systems are considered, a set of rules develop that can be used to transform a set of experimental data to a proposed model network structure, with activation factors yielding a measure of activation potential between intermediates.Singapore-MIT Alliance (SMA

A Functional Protein Chip for Combinatorial Pathway Optimization and In Vitro Metabolic Engineering

Pathway optimization is, in general, a very demanding task due to the complex, nonlinear and largely unknown interactions of enzymes, regulators and metabolites. While in vitro reconstruction and pathway analysis is a viable alternative, a major limitation of this approach is the availability of the pathway enzymes for reliable pathway reconstruction. Here, we report the application of RNA display methods for the construction of fusion (chimeric) molecules, comprising mRNA and the protein they express, that can be used for the above purpose. The chimeric molecule is immobilized via hybridization of its mRNA end with homologous capture DNA spotted on a substrate surface. We show that the protein (enzyme) end of the fusion molecule retains its function under immobilized conditions and that the enzymatic activity is proportional to the amount of capture DNA spotted on the surface of a microarray or 96-well microplate. The relative amounts of all pathway enzymes can thus be changed at will by changing the amount of the corresponding capture DNA. Hence, entire pathways can be reconstructed and optimized in vitro from genomic information alone by generating chimeric molecules for all pathway enzymes in a single in vitro translation step and hybridizing on 96-well microplates where each well contains a different combination of capture DNA. We provide validation of this concept with the sequential reactions catalyzed by luciferase and nucleoside diphosphate kinase and further illustrate this method with the optimization of the five-step pathway for trehalose synthesis. Multi-enzyme pathways leading to the synthesis of specialty molecules can thus be optimized from genomic information about the pathway enzymes, provided the latter retain their activity under the in vitro immobilized conditions.Singapore-MIT Alliance (SMA

Machine Learning Approaches to Modeling the Physiochemical Properties of Small Peptides

Peptide and protein sequences are most commonly represented as a strings: a series of letters selected from the twenty character alphabet of abbreviations for the naturally occurring amino acids. Here, we experiment with representations of small peptide sequences that incorporate more physiochemical information. Specifically, we develop three different physiochemical representations for a set of roughly 700 HIV–I protease substrates. These different representations are used as input to an array of six different machine learning models which are used to predict whether or not a given peptide is likely to be an acceptable substrate for the protease. Our results show that, in general, higher–dimensional physiochemical representations tend to have better performance than representations incorporating fewer dimensions selected on the basis of high information content. We contend that such representations are more biologically relevant than simple string–based representations and are likely to more accurately capture peptide characteristics that are functionally important.Singapore-MIT Alliance (SMA

Multiresolution analysis in statistical mechanics. II. The wavelet transform as a basis for Monte Carlo simulations on lattices

In this paper, we extend our analysis of lattice systems using the wavelet transform to systems for which exact enumeration is impractical. For such systems, we illustrate a wavelet-accelerated Monte Carlo (WAMC) algorithm, which hierarchically coarse-grains a lattice model by computing the probability distribution for successively larger block spins. We demonstrate that although the method perturbs the system by changing its Hamiltonian and by allowing block spins to take on values not permitted for individual spins, the results obtained agree with the analytical results in the preceding paper, and ``converge'' to exact results obtained in the absence of coarse-graining. Additionally, we show that the decorrelation time for the WAMC is no worse than that of Metropolis Monte Carlo (MMC), and that scaling laws can be constructed from data performed in several short simulations to estimate the results that would be obtained from the original simulation. Although the algorithm is not asymptotically faster than traditional MMC, because of its hierarchical design, the new algorithm executes several orders of magnitude faster than a full simulation of the original problem. Consequently, the new method allows for rapid analysis of a phase diagram, allowing computational time to be focused on regions near phase transitions.Comment: 11 pages plus 7 figures in PNG format (downloadable separately

Inverse Metabolic Engineering of Synechocystis PCC 6803 for Improved Growth Rate and Poly-3-hydroxybutyrate Production

Synechocystis PCC 6803 is a photosynthetic bacterium that has the potential to make bioproducts from carbon dioxide and light. Biochemical production from photosynthetic organisms is attractive because it replaces the typical bioprocessing steps of crop growth, milling, and fermentation, with a one-step photosynthetic process. However, low yields and slow growth rates limit the economic potential of such endeavors. Rational metabolic engineering methods are hindered by limited cellular knowledge and inadequate models of Synechocystis. Instead, inverse metabolic engineering, a scheme based on combinatorial gene searches which does not require detailed cellular models, but can exploit sequence data and existing molecular biological techniques, was used to find genes that (1) improve the production of the biopolymer poly-3-hydroxybutyrate (PHB) and (2) increase the growth rate. A fluorescence activated cell sorting assay was developed to screen for high PHB producing clones. Separately, serial sub-culturing was used to select clones that improve growth rate. Novel gene knock-outs were identified that increase PHB production and others that increase the specific growth rate. These improvements make this system more attractive for industrial use and demonstrate the power of inverse metabolic engineering to identify novel phenotype-associated genes in poorly understood systems.Singapore-MIT Alliance (SMA

Parameter estimation of kinetic models from metabolic profiles: two-phase dynamic decoupling method

Motivation: Time-series measurements of metabolite concentration have become increasingly more common, providing data for building kinetic models of metabolic networks using ordinary differential equations (ODEs). In practice, however, such time-course data are usually incomplete and noisy, and the estimation of kinetic parameters from these data is challenging. Practical limitations due to data and computational aspects, such as solving stiff ODEs and finding global optimal solution to the estimation problem, give motivations to develop a new estimation procedure that can circumvent some of these constraints. Results: In this work, an incremental and iterative parameter estimation method is proposed that combines and iterates between two estimation phases. One phase involves a decoupling method, in which a subset of model parameters that are associated with measured metabolites, are estimated using the minimization of slope errors. Another phase follows, in which the ODE model is solved one equation at a time and the remaining model parameters are obtained by minimizing concentration errors. The performance of this two-phase method was tested on a generic branched metabolic pathway and the glycolytic pathway of Lactococcus lactis. The results showed that the method is efficient in getting accurate parameter estimates, even when some information is missing. Contact: [email protected] Supplementary information: Supplementary data are available at Bioinformatics onlin

Downstream reactions and engineering in the microbially reconstituted pathway for Taxol

Taxol (a trademarked product of Bristol-Myers Squibb) is a complex isoprenoid natural product which has displayed potent anticancer activity. Originally isolated from the Pacific yew tree (Taxus brevifolia), Taxol has been mass-produced through processes reliant on plant-derived biosynthesis. Recently, there have been alternative efforts to reconstitute the biosynthetic process through technically convenient microbial hosts, which offer unmatched growth kinetics and engineering potential. Such an approach is made challenging by the need to successfully introduce the significantly foreign enzymatic steps responsible for eventual biosynthesis. Doing so, however, offers the potential to engineer more efficient and economical production processes and the opportunity to design and produce tailored analog compounds with enhanced properties. This mini review will specifically focus on heterologous biosynthesis as it applies to Taxol with an emphasis on the challenges associated with introducing and reconstituting the downstream reaction steps needed for final bioactivity.National Institutes of Health (U.S.) (GM085323)Milheim Foundation (2006-2017

Metabolomic and 13C-metabolic flux analysis of a xylose-consuming Saccharomyces cerevisiae strain expressing xylose isomerase

Over the past two decades, significant progress has been made in the engineering of xylose-consuming Saccharomyces cerevisiae strains for production of lignocellulosic biofuels. However, the ethanol productivities achieved on xylose are still significantly lower than those observed on glucose for reasons that are not well understood. We have undertaken an analysis of central carbon metabolite pool sizes and metabolic fluxes on glucose and on xylose under aerobic and anaerobic conditions in a strain capable of rapid xylose assimilation via xylose isomerase in order to investigate factors that may limit the rate of xylose fermentation. We find that during xylose utilization the flux through the non-oxidative Pentose Phosphate Pathway (PPP) is high but the flux through the oxidative PPP is low, highlighting an advantage of the strain employed in this study. Furthermore, xylose fails to elicit the full carbon catabolite repression response that is characteristic of glucose fermentation in S. cerevisiae. We present indirect evidence that the incomplete activation of the fermentation program on xylose results in a bottleneck in lower glycolysis, leading to inefficient re-oxidation of NADH produced in glycolysis.Shell Oil CompanyNational Institute of General Medical Sciences (U.S.) Biotechnology Training Progra