Selection of tracers for 13C-metabolic flux analysis using elementary metabolite units (EMU) basis vector methodology

Metabolic flux analysis (MFA) is a powerful technique for elucidating in vivo fluxes in microbial and mammalian systems. A key step in (13)C-MFA is the selection of an appropriate isotopic tracer to observe fluxes in a proposed network model. Despite the importance of MFA in metabolic engineering and beyond, current approaches for tracer experiment design are still largely based on trial-and-error. The lack of a rational methodology for selecting isotopic tracers prevents MFA from achieving its full potential. Here, we introduce a new technique for tracer experiment design based on the concept of elementary metabolite unit (EMU) basis vectors. We demonstrate that any metabolite in a network model can be expressed as a linear combination of so-called EMU basis vectors, where the corresponding coefficients indicate the fractional contribution of the EMU basis vector to the product metabolite. The strength of this approach is the decoupling of substrate labeling, i.e. the EMU basis vectors, from the dependence on free fluxes, i.e. the coefficients. In this work, we demonstrate that flux observability inherently depends on the number of independent EMU basis vectors and the sensitivities of coefficients with respect to free fluxes. Specifically, the number of independent EMU basis vectors places hard limits on how many free fluxes can be determined in a model. This constraint is used as a guide for selecting feasible substrate labeling. In three example models, we demonstrate that by maximizing the number of independent EMU basis vectors the observability of a system is improved. Inspection of sensitivities of coefficients with respect to free fluxes provides additional constraints for proper selection of tracers. The present contribution provides a fresh perspective on an important topic in metabolic engineering, and gives practical guidelines and design principles for a priori selection of isotopic tracers for (13)C-MFA studies.     

Regional management units for marine turtles: a novel framework for prioritizing conservation and research across multiple scales

Background

Resolving threats to widely distributed marine megafauna requires definition of the geographic distributions of both the threats as well as the population unit(s) of interest. In turn, because individual threats can operate on varying spatial scales, their impacts can affect different segments of a population of the same species. Therefore, integration of multiple tools and techniques — including site-based monitoring, genetic analyses, mark-recapture studies and telemetry — can facilitate robust definitions of population segments at multiple biological and spatial scales to address different management and research challenges.

Methodology/Principal Findings

To address these issues for marine turtles, we collated all available studies on marine turtle biogeography, including nesting sites, population abundances and trends, population genetics, and satellite telemetry. We georeferenced this information to generate separate layers for nesting sites, genetic stocks, and core distributions of population segments of all marine turtle species. We then spatially integrated this information from fine- to coarse-spatial scales to develop nested envelope models, or Regional Management Units (RMUs), for marine turtles globally.

Conclusions/Significance

The RMU framework is a solution to the challenge of how to organize marine turtles into units of protection above the level of nesting populations, but below the level of species, within regional entities that might be on independent evolutionary trajectories. Among many potential applications, RMUs provide a framework for identifying data gaps, assessing high diversity areas for multiple species and genetic stocks, and evaluating conservation status of marine turtles. Furthermore, RMUs allow for identification of geographic barriers to gene flow, and can provide valuable guidance to marine spatial planning initiatives that integrate spatial distributions of protected species and human activities. In addition, the RMU framework — including maps and supporting metadata — will be an iterative, user-driven tool made publicly available in an online application for comments, improvements, download and analysis.     

An elementary metabolite unit (EMU) based method of isotopically nonstationary flux analysis

Nonstationary metabolic flux analysis (NMFA) is at present a very computationally intensive exercise, especially for large reaction networks. We applied elementary metabolite unit (EMU) theory to NMFA, dramatically reducing computational difficulty. We also introduced block decoupling, a new method that systematically and comprehensively divides EMU systems of equations into smaller subproblems to further reduce computational difficulty. These improvements led to a 5000-fold reduction in simulation times, enabling an entirely new and more complicated set of problems to be analyzed with NMFA. We simulated a series of nonstationary and stationary GC/MS measurements for a large E. coli network that was then used to estimate parameters and their associated confidence intervals. We found that fluxes could be successfully estimated using only nonstationary labeling data and external flux measurements. Addition of near-stationary and stationary time points increased the precision of most parameters. Contrary to prior reports, the precision of nonstationary estimates proved to be comparable to the precision of estimates based solely on stationary data. Finally, we applied EMU-based NMFA to experimental nonstationary measurements taken from brown adipocytes and successfully estimated fluxes and some metabolite concentrations. By using NFMA instead of traditional MFA, the experiment required only 6 h instead of 50 (the time necessary for most metabolite labeling to reach 99% of isotopic steady state).     

Characterization of a novel (±)-abscisic acid metabolite

By application of (±)-abscisic acid and (±)-[1-¹⁴C]-abscisic acid to pea seedlings ( Pisum sativum L. cv. Kleine Rheinländerin), a new abscisic acid metabolite (Pisumic acid, PISA) could be isolated and structurally characterized by combined gas chromatography-mass spectrometry. From data of exact mass determination, it is suggested that the metabolite has the tentative structure of 4'-dihydroabscisic acid with a hydroxylated methyl group at C-6'. That this could be evidence for an alternative pathway of abscisic acid metabolism which was suggested earlier by Walton et al. [Planta 112: 87–90 (1973)] is discussed.     

肠道菌群及其代谢产物氧化三甲胺——冠心病治疗的新靶点

冠心病 (Coronary artery disease,CAD) 是全球发病率和死亡率最高的一种心血管疾病,冠心病和肠道菌群失调密切相关,肠道菌群可能是未来冠心病的重要诊断标志物,改善肠道菌群微环境有望成为治疗冠心病的新途径。作为肠道菌群参与合成的活性代谢产物,氧化三甲胺 (Trimethylamine-N-oxide,TMAO) 水平的升高与心血管疾病患病风险、全因死亡率的增加有关;基础研究表明TMAO可能具有促动脉粥样硬化特性;这些研究提示TMAO可作为预防和治疗冠心病的潜在靶点。文中分析了当前调控肠道菌群及其代谢产物TMAO治疗冠心病的临床及基础性研究,以期为冠心病的治疗提供帮助。     

Following the genes: a framework for animal modeling of psychiatric disorders

Background

The caspase family, which plays a central role in apoptosis in metazoans, has undergone an expansion in amphioxus, increasing to 45 members through domain recombination and shuffling.

Results

In order to shed light on the conservation and uniqueness of this family in amphioxus, we cloned three representative caspase genes, designated as bbtCaspase-8, bbtCaspase-1/2 and bbtCaspase3-like, from the amphioxus Branchiostoma belcheri tsingtauense. We found that bbtCaspase-8 with conserved protein architecture is involved in the Fas-associated death domain-Caspase-8 mediated pro-apoptotic extrinsic pathway, while bbtCaspase3-like may mediate a nuclear apoptotic pathway in amphioxus. Also, bbtCaspase-1/2 can co-localize with bbtFADD2 in the nucleus, and be recruited to the cytoplasm by amphioxus apoptosis associated speck-like proteins containing a caspase recruitment domain, indicating that bbtCaspase-1/2 may serve as a switch between apoptosis and caspase-dependent innate immune response in invertebrates. Finally, amphioxus extrinsic apoptotic pathway related caspases played important roles in early embryogenesis.

Conclusions

Our study not only demonstrates the conservation of bbtCaspase-8 in apoptosis, but also reveals the unique features of several amphioxus caspases with novel domain architectures arose some 500 million years ago.     

An activated sludge modeling framework for xenobiotic trace chemicals (ASM-X): Assessment of diclofenac and carbamazepine

Conventional models for predicting the fate of xenobiotic organic trace chemicals, identified, and calibrated using data obtained in batch experiments spiked with reference substances, can be limited in predicting xenobiotic removal in wastewater treatment plants (WWTPs). At stake is the level of model complexity required to adequately describe a general theory of xenobiotic removal in WWTPs. In this article, we assess the factors that influence the removal of diclofenac and carbamazepine in activated sludge, and evaluate the complexity required for the model to effectively predict their removal. The results are generalized to previously published cases. Batch experimental results, obtained under anoxic and aerobic conditions, were used to identify extensions to, and to estimate parameter values of the activated sludge modeling framework for Xenobiotic trace chemicals (ASM‐X). Measurement and simulation results obtained in the batch experiments, spiked with the diclofenac and carbamazepine content of preclarified municipal wastewater shows comparably high biotransformation rates in the presence of growth substrates. Forward dynamic simulations were performed using full‐scale data obtained from Bekkelaget WWTP (Oslo, Norway) to evaluate the model and to estimate the level of re‐transformable xenobiotics present in the influent. The results obtained in this study demonstrate that xenobiotic loading conditions can significantly influence the removal capacity of WWTPs. We show that the trace chemical retransformation in upstream sewer pipes can introduce considerable error in assessing the removal efficiency of a WWTP, based only on parent compound concentration measurements. The combination of our data with those from the literature shows that solids retention time (SRT) can enhance the biotransformation of diclofenac, which was not the case for carbamazepine. Model approximation of the xenobiotic concentration, detected in the solid phase, suggest that between approximately 1% and 16% of the total solid carbamazepine and diclofenac concentrations, respectively, is due to sorption—the remainder being non‐bioavailable and sequestered. We demonstrate the effectiveness of the model's predictive power over conventional tools in a statistical analysis, performed at four levels of structural complexity. To assess WWTP retrofitting needs to remove xenobiotic trace chemicals, we suggest using mechanistic models, e.g., ASM‐X, in regional risk assessments. For preliminary evaluations, we present operating charts that can be used to estimate average xenobiotic removal rates in WWTPs as a function of SRT and the xenobiotics mass loads normalised to design treatment capacity. Biotechnol. Bioeng. 2012; 109: 2757–2769. © 2012 Wiley Periodicals, Inc.     

A novel conotoxin framework with a helix-loop-helix (Cs alpha/alpha) fold

Venomous predatory animals, such as snakes, spiders, scorpions, sea anemones, and cone snails, produce a variety of highly stable cystine-constrained peptide scaffolds as part of their neurochemical strategy for capturing prey. Here we report a new family of four-cystine, three-loop conotoxins (designated framework 14). Three peptides of this family (flf14a-c) were isolated from the venom of Conus floridanus floridensis, and one (vil14a) was isolated from the venom of Conus villepinii, two worm-hunting Western Atlantic cone snail species. The primary structure for these peptides was determined using Edman degradation sequencing, and their cystine pairing was assessed by limited hydrolysis with a combination of CNBr and chymotrypsin under nonreducing, nonalkylating conditions in combination with MALDI-TOF MS analysis of the resulting peptidic fragments. CD spectra and nanoNMR spectroscopy of these conotoxins directly isolated from the cone snails revealed a highly helical secondary structure for the four conotoxins. Sequence-specific nanoNMR analysis at room temperature revealed a well-defined helix-loop-helix tertiary structure that resembles that of the Cs alpha/alpha scorpion toxins kappa-hefutoxin, kappa-KTx1.3, and Om-toxins, which adopt a stable three-dimensional fold where the two alpha-helices are linked by the two disulfide bridges. One of these conotoxins (vil14a) has a Lys/Tyr dyad, separated by approximately 6A, which is a conserved structural feature in K(+) channel blockers. The presence of this framework in scorpions and in cone snails indicates a common molecular imprint in the venom of apparently unrelated predatory animals and suggests a common ancestral genetic origin.     

Metabolism of 2-(2-thienyl)allylamine hydrochloride in the rat: identification of a novel metabolite

A novel metabolite, 2-(2-thienyl)propionic acid, is formed in vivo from 2-(2-thienyl)allylamine hydrochloride. Mass spectral analysis suggested 2-(2-thienyl)propionic acid formation involves loss of the amine moiety followed by reduction of the olefinic group.     

A novel framework for spatial normalization of dose distributions in voxel-based analyses of brain irradiation outcomes

PurposeTo devise a novel Spatial Normalization framework for Voxel-based analysis (VBA) in brain radiotherapy. VBAs rely on accurate spatial normalization of different patients’ planning CTs on a common coordinate system (CCS). The cerebral anatomy, well characterized by MRI, shows instead poor contrast in CT, resulting in potential inaccuracies in VBAs based on CT alone.MethodsWe analyzed 50 meningioma patients treated with proton-therapy, undergoing planning CT and T1-weighted (T1w) MRI. The spatial normalization pipeline based on MR and CT images consisted in: intra-patient registration of CT to T1w, inter-patient registration of T1w to MNI space chosen as CCS, doses propagation to MNI.The registration quality was compared with that obtained by Statistical Parametric Mapping software (SPM), used as benchmark. To evaluate the accuracy of dose normalization, the dose organ overlap (DOO) score was computed on gray matter, white matter and cerebrospinal fluid before and after normalization. In addition, the trends in the DOOs distribution were investigated by means of cluster analysis.ResultsThe registration quality was higher for the proposed method compared to SPM (p < 0.001). The DOO scores showed a significant improvement after normalization (p < 0.001). The cluster analysis highlighted 2 clusters, with one of them including the majority of data and exhibiting acceptable DOOs.ConclusionsOur study presents a robust tool for spatial normalization, specifically tailored for brain dose VBAs. Furthermore, the cluster analysis provides a formal criterion for patient exclusion in case of non-acceptable normalization results. The implemented framework lays the groundwork for future reliable VBAs in brain irradiation studies.     

Detection of a novel cyclooxygenase metabolite produced by human promyelocytic leukemia (HL-60) cells

Arachidonic acid metabolism via the lipoxygenase pathway was examined in HL-60 cells before and after N,N-dimethylformamide induced differentiation along granulocytic lines. Untreated HL-60 cells produced small amounts of the 5-lipoxygenase products, 5-hydroxy-eicosatetraenoic acid and leukotriene B4 upon stimulation with calcium ionophore A23187. N,N-dimethylformamide treatment, caused a 10 to 20 fold increase in the amount of ionophore A23187-induced 5-lipoxygenase metabolites. An additional, and as yet unidentified arachidonic acid metabolite was routinely observed during reverse-phase high pressure liquid chromatography analyses of lipoxygenase products. Sensitivity to inhibition by less than 10(-7)M indomethacin coupled with other characteristics of its production, strongly suggest the compound is a cyclooxygenase product. The unusual UV absorbance and chromatographic elution pattern, however, suggest that it is not a typical prostaglandin, thromboxane or prostacyclin product.     

MollDE: a homology modeling framework you can click with

SUMMARY: Molecular Integrated Development Environment (MolIDE) is an integrated application designed to provide homology modeling tools and protocols under a uniform, user-friendly graphical interface. Its main purpose is to combine the most frequent modeling steps in a semi-automatic, interactive way, guiding the user from the target protein sequence to the final three-dimensional protein structure. The typical basic homology modeling process is composed of building sequence profiles of the target sequence family, secondary structure prediction, sequence alignment with PDB structures, assisted alignment editing, side-chain prediction and loop building. All of these steps are available through a graphical user interface. MolIDE's user-friendly and streamlined interactive modeling protocol allows the user to focus on the important modeling questions, hiding from the user the raw data generation and conversion steps. MolIDE was designed from the ground up as an open-source, cross-platform, extensible framework. This allows developers to integrate additional third-party programs to MolIDE. AVAILABILITY: http://dunbrack.fccc.edu/molide/molide.php CONTACT: rl_dunbrack@fccc.edu.     

High-pressure hydrocephalus: a novel analytical modeling approach

Hydrocephalus is an abnormal accumulation of cerebrospinal fluid (CSF) within ventricles and subarachnoid space (SAS) as a result of disturbances in secretion or absorption procedures. It is believed that arachnoid villi cells, which are microscopic projections of pia-arachnoid mater that extend into venous channels in sagittal sinus, are the main sites for CSF absorption, but it is tempting to speculate that a significant portion of CSF is removed from the SAS by nasal lymphatic vessels around olfactory nerve. Thus, in this paper, we propose an analytical model of CSF-lymphatic-blood circulation, in which these two output pathways for CSF absorption have been considered. Mathematical relations governing the pressures in different interacting compartments of the brain are considered. In addition, for increasing the similarity of our model to the physiological conditions, the bulk flow mechanism, which is supposed to occur during CSF absorption, has been considered in our model. We used our model to simulate hydrocephalus. The results indicate that the lymphatic disorders have more considerable effect in decreasing CSF absorption, compared to the disturbances in arachnoid villi cells. Based on our modeling, we believe that disorders in lymphatic pathway may be a cause of high-pressure hydrocephalus. Surely experimental studies are required to validate our hypothesis.     

Envelope: interactive software for modeling and fitting complex isotope distributions

Background  

An important aspect of proteomic mass spectrometry involves quantifying and interpreting the isotope distributions arising from mixtures of macromolecules with different isotope labeling patterns. These patterns can be quite complex, in particular with in vivo metabolic labeling experiments producing fractional atomic labeling or fractional residue labeling of peptides or other macromolecules. In general, it can be difficult to distinguish the contributions of species with different labeling patterns to an experimental spectrum and difficult to calculate a theoretical isotope distribution to fit such data. There is a need for interactive and user-friendly software that can calculate and fit the entire isotope distribution of a complex mixture while comparing these calculations with experimental data and extracting the contributions from the differently labeled species.     

Evolutionary models for insertions and deletions in a probabilistic modeling framework

Background  

Probabilistic models for sequence comparison (such as hidden Markov models and pair hidden Markov models for proteins and mRNAs, or their context-free grammar counterparts for structural RNAs) often assume a fixed degree of divergence. Ideally we would like these models to be conditional on evolutionary divergence time.     

Dragon Promoter Mapper (DPM): a Bayesian framework for modelling promoter structures

SUMMARY: Dragon Promoter Mapper (DPM) is a tool to model promoter structure of co-regulated genes using methodology of Bayesian networks. DPM exploits an exhaustive set of motif features (such as motif, its strand, the order of motif occurrence and mutual distance between the adjacent motifs) and generates models from the target promoter sequences, which may be used to (1) detect regions in a genomic sequence which are similar to the target promoters or (2) to classify other promoters as similar or not to the target promoter group. DPM can also be used for modelling of enhancers and silencers. AVAILABILITY: http://defiant.i2r.a-star.edu.sg/projects/BayesPromoter/ CONTACT: vlad@sanbi.ac.za SUPPLEMENTARY INFORMATION: Manual for using DPM web server is provided at http://defiant.i2r.a-star.edu.sg/projects/BayesPromoter/html/manual/manual.htm.     

GC/MS based metabolomics: development of a data mining system for metabolite identification by using soft independent modeling of class analogy (SIMCA)

Background  

The goal of metabolomics analyses is a comprehensive and systematic understanding of all metabolites in biological samples. Many useful platforms have been developed to achieve this goal. Gas chromatography coupled to mass spectrometry (GC/MS) is a well-established analytical method in metabolomics study, and 200 to 500 peaks are routinely observed with one biological sample. However, only ~100 metabolites can be identified, and the remaining peaks are left as "unknowns".     

Transducers: an emerging probabilistic framework for modeling indels on trees

Contact: ihh{at}berkeley.edu Associate Editor: Alex Bateman