ACoM: A classification method for elementary flux modes based on motif finding

Elementary flux mode analysis is a powerful tool for the theoretical study of metabolic networks. However, when the networks are complex, the determination of elementary flux modes leads to combinatorial explosion of their number which prevents from drawing simple conclusions from their analysis. To deal with this problem we have developed a method based on the Agglomeration of Common Motifs (ACoM) for classifying elementary flux modes. We applied this algorithm to describe the decomposition into elementary flux modes of the central carbon metabolism in Bacillus subtilis and of the yeast mitochondrial energy metabolism. ACoM helps to give biological meaning to the different elementary flux modes and to the relatedness between reactions. ACoM, which can be viewed as a bi-clustering method, can be of general use for sets of vectors with values 0, +1 or −1.     

Combinatorial complexity of pathway analysis in metabolic networks

Elementary flux mode analysis is a promising approach for a pathway-oriented perspective of metabolic networks. However, in larger networks it is hampered by the combinatorial explosion of possible routes. In this work we give some estimations on the combinatorial complexity including theoretical upper bounds for the number of elementary flux modes in a network of a given size. In a case study, we computed the elementary modes in the central metabolism of Escherichia coli while utilizing four different substrates. Interestingly, although the number of modes occurring in this complex network can exceed half a million, it is still far below the upper bound. Hence, to a certain extent, pathway analysis of central catabolism is feasible to assess network properties such as flexibility and functionality.     

Headspace analyses of acetone: A rapid method for measuring the H-labeling of water

Measurements of a ²H-labeling of water in biological fluids are required for determining the rates of biochemical flux and for estimating body composition. We have been using the method which relies on the base-catalyzed exchange of hydrogen (deuterium) between water and acetone. ²H-labeling of acetone is then determined using GCMS. Although not noted in the original paper, when chloroform is used to extract the acetone there is slow but substantial backexchange between [²H]acetone and solvent (unpublished observations). We report herein on a refinement of the assay that utilizes headspace analysis, which minimizes the number of transfers and decreases sample preparation time and instrument run time.     

Robustness of circadian rhythms in the presence of molecular fluctuations: an investigation based on a mechanistic, statistical theory and a simulation algorithm

After a very brief introduction to a mechanistic and statistical theory of molecular fluctuations in chemical reactions developed by Joel Keizer, we explore the robustness of a circadian rhythm model by using the theory and the exact stochastic simulation (ESS). The comparative study shows that the theory reflects the effects of the dynamics of the model on the robustness more than ESS does. Even though the theory is a macroscopic one, the robustness of the model compares well with that computed from the ESS when the system size is larger than 50. The robustness increases nonlinearly with the system size and it reaches an asymptotic value at higher system sizes. As we can expect from the dynamics of the system, the robustness is minimum near the bifurcation point and as the most sensitive parameter increases away from the bifurcation point the robustness according to the theory as well as the ESS increases and then reaches to a steady value.     

Isotopic non-stationary 13C gluconate tracer method for accurate determination of the pentose phosphate pathway split-ratio in Penicillium chrysogenum

Current (13)C labeling experiments for metabolic flux analysis (MFA) are mostly limited by either the requirement of isotopic steady state or the extremely high computational effort due to the size and complexity of large metabolic networks. The presented novel approach circumvents these limitations by applying the isotopic non-stationary approach to a local metabolic network. The procedure is demonstrated in a study of the pentose phosphate pathway (PPP) split-ratio of Penicillium chrysogenum in a penicillin-G producing chemostat-culture grown aerobically at a dilution rate of 0.06h(-1) on glucose, using a tracer amount of uniformly labeled [U-(13)C(6)] gluconate. The rate of labeling inflow can be controlled by using different cell densities and/or different fractions of the labeled tracer in the feed. Due to the simplicity of the local metabolic network structure around the 6-phosphogluconate (6pg) node, only three metabolites need to be measured for the pool size and isotopomer distribution. Furthermore, the mathematical modeling of isotopomer distributions for the flux estimation has been reduced from large scale differential equations to algebraic equations. Under the studied cultivation condition, the estimated split-ratio (41.2+/-0.6%) using the novel approach, shows statistically no difference with the split-ratio obtained from the originally proposed isotopic stationary gluconate tracing method.     

A method for the selection of production media for actinomycete strains based on their metabolite HPLC profiles

The manipulation of growth conditions of microorganisms is a common strategy used by pharmaceutical companies to improve the quantities and spectra of secondary metabolites with potential therapeutic interest. In this work, the effects of fermentation media on secondary metabolite production from a set of Actinomycetes was statistically compared. For this purpose, we created an automated method for comparing the ability of microorganisms to produce different secondary metabolites. HPLC analyses guided the selection of those media in which a wider chemical diversity was obtained from microorganisms inoculated in a wide spectrum of production media. Fermented media yielding a better secondary metabolite profile were included in subsequent drug discovery screening.     

Global analysis of metabolites in rat and human urine based on gas chromatography/time-of-flight mass spectrometry

Sediment in urine may contain low-molecular-weight compounds that should be included in the analysis. To date, no systematic investigation has addressed this issue. We investigated three primary factors that influence the extraction efficiency of metabolites during preparation of urine samples for metabolomic research: centrifugation, pH, and extraction solvents. Obtained with the use of gas chromatography/time-of-flight mass spectrometry (GC/TOFMS) technique and principal component analysis (PCA), our results indicate that (1) conventional centrifugation causes an apparent loss of some metabolites, indicating that urine samples for metabolomic research should not be centrifuged before procedures are undertaken to recover the metabolites; (2) pH adjustment has a large impact on the recovery of metabolites and is therefore not encouraged; (3) with design of experiment analysis, methanol and water yield the optimal extraction efficiency. Differences between rat and human urine were observed and are discussed. Ninety-nine metabolites identified in rat and human urine are presented. An efficient protocol is proposed for the pretreatment of urine samples.     

基于分组重量编码的蛋白质功能预测

从蛋白质序列出发,采用分组重量编码(Encoding Based on Grouped Weight,简记EBGW),并结合最近邻居算法对蛋白质功能进行预测。对酵母(Saccharomyces cerevisiae)蛋白质的1826条序列进行预测,整体预测准确率与其他基于序列信息的蛋白质功能预测方法相当。实验结果表明基于EBGW编码方案的新方法可有效地应用于蛋白质功能预测。     

13C‐metabolic flux analysis for batch culture of Escherichia coli and its pyk and pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites

Since most bio‐production processes are conducted in a batch or fed‐batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied ¹³C‐metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild‐type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild‐type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose‐6‐phosphate dehydrogenase reaction was the rate‐limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild‐type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

一种基于PCR分析多样性的小鼠肠道微生物宏基因组提取方法

目的建立一种基于PCR分析分子多样性的小鼠肠道菌群宏基因组提取方法。方法比较、综合国内外小鼠肠道菌群宏基因组的提取方法后建立一种新方法,小鼠肠道内容物经丙酮洗涤,差速离心,溶菌酶、SDS裂解,CTAB处理,酚／氯仿抽提后可得到高质量的DNA,通过紫外分光光度计、琼脂糖凝胶电泳、细菌通用引物PCR和扩增核糖体限制性酶切片段分析（ARDRA）等检测该方法的实用性。结果该方法获得的小鼠肠道菌群宏基因组DNA大小在23kb左右,A260／A280在1．8—2．0,经细菌通用引物PCR后能得到适用于ARDRA的目的产物。结论该方法经济适用性较强,具备一定的应用价值。     

A new method for non-destructive measurement of biomass, growth rates, vertical biomass distribution and dry matter content based on digital image analysis

BACKGROUND AND AIMS: Biomass is an important trait in functional ecology and growth analysis. The typical methods for measuring biomass are destructive. Thus, they do not allow the development of individual plants to be followed and they require many individuals to be cultivated for repeated measurements. Non-destructive methods do not have these limitations. Here, a non-destructive method based on digital image analysis is presented, addressing not only above-ground fresh biomass (FBM) and oven-dried biomass (DBM), but also vertical biomass distribution as well as dry matter content (DMC) and growth rates. METHODS: Scaled digital images of the plants silhouettes were taken for 582 individuals of 27 grass species (Poaceae). Above-ground biomass and DMC were measured using destructive methods. With image analysis software Zeiss KS 300, the projected area and the proportion of greenish pixels were calculated, and generalized linear models (GLMs) were developed with destructively measured parameters as dependent variables and parameters derived from image analysis as independent variables. A bootstrap analysis was performed to assess the number of individuals required for re-calibration of the models. KEY RESULTS: The results of the developed models showed no systematic errors compared with traditionally measured values and explained most of their variance (R(2) > or = 0.85 for all models). The presented models can be directly applied to herbaceous grasses without further calibration. Applying the models to other growth forms might require a re-calibration which can be based on only 10-20 individuals for FBM or DMC and on 40-50 individuals for DBM. CONCLUSIONS: The methods presented are time and cost effective compared with traditional methods, especially if development or growth rates are to be measured repeatedly. Hence, they offer an alternative way of determining biomass, especially as they are non-destructive and address not only FBM and DBM, but also vertical biomass distribution and DMC.     

A pretreatment method for the determination of nitrate in brackish water and seawater based on the hydrazinium reduction technique

The hydrazinium reduction technique has so far been inapplicable to the determination of nitrate in brackish water and seawater due to interference from magnesium ions. We developed a pretreatment method for brackish water and seawater samples for the determination of nitrate based on the hydrazinium reduction technique. Magnesium ions in water samples were converted to a precipitate of a complex with oxine (8-quinolinol) at pH 9.5, and then the precipitate was centrifuged at 3000 rpm for 20 min. The extra oxine in the resulting sample (the supernatant liquid), which inhibits the reduction of nitrate to nitrite, was removed using a Sep-Pak C18 cartridge. Thus we achieved the preparation of a magnesium-free water sample. Using the hydrazinium reduction technique with the proposed pretreatment method, nitrate in brackish water and seawater as well as fresh water was quantitatively determined with high accuracy. Received: July 21, 1999 / Accepted: September 26, 1999     

A method for reverse interactome analysis: High‐resolution mapping of interdomain interaction network in Dam1 complex and its specific disorganization based on the interaction domain expression

An experimental methodology that facilitates functional analysis of numerous protein–protein interactions, which have been found in genome‐wide interactome researches, has long been awaited. We propose herein an antagonistic inhibition‐based approach. The antagonizing polypeptide is generated in the course of interaction domain mapping based on yeast 2‐hybrid (Y2H) screening coupled with in vitro convergence of the Y2H‐selected fragments, which is performed in a formatted procedure. Using the coupled methodology, we first performed a high‐resolution mapping of an interdomain interaction network within budding yeast's Dam1 complex. Dam1 complex is a kinetochore protein complex composed of 10 essential subunits including Spc34p and Spc19p. The high‐resolution mapping revealed the overall network structure within the complex for the first time: Dam1 components form into two separated subnetworks on N‐terminal scaffolding domains of Spc34p and Spc19p, and the coiled‐coil interaction in their C‐terminal domains connects the subnetworks. Secondly, we show that the domain fragments converged in the high‐resolution mapping acted as potent inhibitors for the endogenous interactions when episomally overexpressed. The in vivo Dam1 interaction targeting with the fragments conferred a similar phenotype on the host cells; a critical and irreversible damage, which was accompanied with disturbed budding and chromosome mis‐segregation as a result of disorganized spindle. These phenotypes were strongly related to the cellular function of the Dam1 complex. The results and approach we demonstrated herein not only shed light on the Dam1 molecular architecture but also pave the road to reverse‐interactome analysis and discoveries of novel drugs that target disease‐related protein–protein interactions. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

A simple liquid chromatographic method based on intramolecular excimer-forming derivatization and fluorescence detection for the determination of tyrosine and tyramine in urine

A liquid chromatographic (LC) method for sensitive and selective fluorometric determination of p-hydroxyphenylethylamino group containing compounds is described. This method is based on an intramolecular excimer-forming fluorescence derivatization with a pyrene reagent, 4-(1-pyrene)butanoyl chloride, followed by reversed-phase LC. The analytes, containing an amino moiety and a phenolic hydroxyl moiety in a molecule, were converted to the corresponding dipyrene-labeled derivatives by one-step derivatization. The dipyrene-labeled derivatives afforded intramolecular excimer fluorescence (440-540 nm), which can clearly be discriminated from the normal fluorescence (360-420 nm) emitted from reagent blanks. The derivatives of tyrosine and tyramine could be separated by reversed-phase LC on ODS column under conditions of isocratic elution. The detection limits (signal-to-noise ratio = 3) for tyrosine and tyramine were 4.5 and 2.6 fmol per 20 microL injection, which corresponded to analyte concentrations of 0.9 and 0.5 nM, respectively.     

A sensitive method for quantification of aceticlastic methanogens and estimation of total methanogenic cells in natural environments based on an analysis of ether-linked glycerolipids

Abstract A highly sensitive method for the quantification of methanogens in anaerobic digestor sludges was developed, based on an analysis of ether-linked glycerolipids. Core lipids were prepared from total lipids by HF treatment and mild methanolysis, and these core lipids were quantified as the corresponding 9-anthroyl derivatives by high-performance liquid chromatography with fluorescence detection. The amounts, in terms of cell carbon content, of Methanosaeta and Methanosarcina were proportional to the amounts of α-hydroxyarchaeol and β-hydroxyarchaeol, respectively. Moreover, the total amount of core lipids was well correlated with the cell mass of aceticlastic and H₂/CO₂-consuming methanogens. The limit of detection for Methanosaeta concilii was 17 ng of cell carbon when the signal/noise ratio was 3. This method allowed us to quantitate aceticlastic methanogens with high accuracy and to make a rough estimate of total methanogenic cells without any interference by the multifarious impurities that are present in anaerobic sludges. These results suggest that the present method will be a useful tool for investigations of methanogenic ecosystems.