Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

Minimization of Incompletely Specified Sequential Machines

A simple yet efficient method for the minimization of incompletely specified sequential machines (ISSMs) is proposed. Precise theorems are developed, as a consequence of which several compatibles can be deleted from consideration at the very first stage in the search for a minimal closed cover. Thus, the computational work is significantly reduced. Initial cardinality of the minimal closed cover is further reduced by a consideration of the maximal compatibles (MC's) only; as a result the method converges to the solution faster than the existing procedures. "Rank" of a compatible is defined. It is shown that ordering the compatibles, in accordance with their rank, reduces the number of comparisons to be made in the search for exclusion of compatibles. The new method is simple, systematic, and programmable. It does not involve any heuristics or intuitive procedures. For small- and medium-sized machines, it canle used for hand computation as well. For one of the illustrative examples used in this paper, 30 out of 40 compatibles can be ignored in accordance with the proposed rules and the remaining 10 compatibles only need be considered for obtaining a minimal solution.

Purification and properties of a DNA polymerase from Mycobacterium tuberculosis H37Rv

DNA polymerase has been purified approximately 2000-fold from Mycobacterium tuberculosis H37Rv. The purified preparation was homogeneous by electrophoretic criteria and has a molecular weight of 135 000. The purified enzyme resembles Escherichia coli polymerase I in its properties, being insensitive to sulfhydryl drugs and possessing 5′,3′-exonuclease activity in addition to polymerase and 3′,5′-exonuclease activities. However, it differs from the latter in its sensitivity to higher salt concentration and DNA intercalating agents such as 8-aminoquinoline. The polymerase exhibited maximal activity between 37–42°C and pH 8.8–9.5. The polymerase was stable for several months below 0°C. However, the 5′,3′-exonuclease activity was more labile. The effects of different metal ions, polyamines and drugs on the polymerase activity are presented.

Dynamics of Circulating Concentrations of Gonadotropins and Ovarian Hormones Throughout the Menstrual Cycle in the Bonnet Monkey: Role of Inhibin A in the Regulation of Follicle-Stimulating Hormone Secretion

In higher primates, increased circulating follicle-stimulating hormone (FSH) levels seen during late menstrual cycle and during menstruation has been suggested to be necessary for initiation of follicular growth, recruitment of follicles and eventually culminating in ovulation of a single follicle. With a view to establish the dynamics of circulating FSH secretion with that of inhibin A (INH A) and progesterone (P-4)secretions during the menstrual cycle, blood was collected daily from bonnet monkeys beginning day 1 of the menstrual cycle up to 35 days. Serum INH A levels were low during early follicular phase, increased significantly coinciding with the mid cycle luteinizing hormone (LH) surge to reach maximal levels during the mid luteal phase before declining at the late luteal phase, essentially paralleling the pattern Of P-4 secretion seen throughout the luteal phase. Circulating FSH levels were low during early and mid luteal phases, but progressively increased during the late luteal phase and remained high for few days after the onset of menses. In another experiment, lutectomy performed during the mid luteal phase resulted in significant decrease in INH A concentration within 2 hr (58.3 +/- 2 vs. 27.3 +/- 3 pg/mL), and a 2- to 3-fold rise in circulating FSH levels by 24 hr (0.20 +/- 0.02 vs. 0.53 +/- 0.14 ng/mL) that remained high until 48 hr postlutectomy. Systemic administration of Cetrorelix (150 mu g/kg body weight), a gonadotropin releasing hormone receptor antagonist, at mid luteal phase in monkeys led to suppression of serum INH A and P-4 concentrations 24 hr post treatment, but circulating FSH levels did not change. Administration of exogenous LH, but not FSH, significantly increased INH A concentration. The results taken together suggest a tight coupling between LH and INH A secretion and that INH A is largely responsible for maintenance of low FSH concentration seen during the luteal phase. Am. J. Primatol. 71:817-824, 2009.

Mechanism of aromatic hydroxylation *1, , *2, , *3: Properties of a model for pyridine nucleotide-dependent flavoprotein hydroxylases

A model (NADH-phenazine methosulfate-O2) formally similar to pyridine nucleotide-dependent flavoprotein hydroxylases catalyzed the hydroxylation of several aromatic compounds. The hydroxylation was maximal at acid pH and was inhibited by ovine Superoxide dismutase, suggesting that perhydroxyl radicals might be intermediates in this process. The stoichiometry of the reaction indicated that a univalent reduction of oxygen was occurring. The correlation between the concentration of semiquinone and hydroxylation, and the inhibition of hydroxylation by ethanol which inhibited semiquinone oxidation, suggested the involvement of phenazine methosulfate-semiquinone. Activation of hydroxylation by Fe3+ and Cu2+ supported the contention that univalently reduced species of oxygen was involved in hydroxylation. Catalase was without effect on the hydroxylation by the model, ruling out H2O2 as an intermediate. A reaction sequence, involving a two-electron reduction of phenazine methosulfate to reduced phenazine methosulfate followed by disproportionation with phenazine methosulfate to generate the semiquinone, was proposed. The semiquinone could donate an electron to O2 to generate O2 which could be subsequently protonated to form the perhydroxyl radical.

Specific requirement of lysophosphatidylcholine for palmitoyl-CoA synthetase of chicken intestine

The presence of palmitoyl-CoA synthetase (EC 6.2.1.3) in the brush borderfree particulate fraction of chicken intestinal mucosa is demonstrated. The enzyme was dependent on the simultaneous presence of lysophosphatidylcholine and Triton X-100 as well as ATP, CoA and Mg2+ for maximal activity. Lysophosphatidylcholine could not be replaced by other lipids. Enzyme preparations solubilized by Triton X-100 or lysophosphatidylcholine were still dependent on the presence of detergents for maximal activity.

State Minimization of Incompletely Specified Sequential Machines

A simple procedure for the state minimization of an incompletely specified sequential machine whose number of internal states is not very large is presented. It introduces the concept of a compatibility graph from which the set of maximal compatibles of the machine can be very conveniently derived. Primary and secondary implication trees associated with each maximal compatible are then constructed. The minimal state machine covering the incompletely specified machine is then obtained from these implication trees.

Succinate and malate improve development of hamster eight-cell embryos in vitro: Confirmation of viability by embryo transfer

The in vitro development of hamster preimplantation embryos is supported by non-glucose energy substrates. To investigate the importance of embryonic metabolism, influence of succinate and malate on the development of hamster 8-cell embryos to blastocysts was examined using a chemically defined protein-free modified hamster embryo culture medium-2 (HECM-2m). There was a dose-dependent influence of succinate on blastocyst development; 0.5 mM succinate was optimal (85.1% ± 3.9 vs. 54.5% ± 3.5). In succinate-supplemented HECM-2m, blastocyst development was reduced by omission of lactate (68.5% ± 7.2), but not pyruvate (85.8% ± 6.2) or glutamine (84.1% ± 2.1). Succinate along with either glutamine or lactate or pyruvate poorly supported blastocyst development (28%-58%). Malate also stimulated blastocyst development; 0.01 mM malate was optimal (86.3% ± 2.8). Supplementation of both succinate and malate to HECM-2m supported maximal (100%) blastocyst development, which was inhibited 4-fold by the addition of glucose/phosphate. The mean cell numbers (MCN) of blastocysts cultured in succinate-supplemented HECM-2m was higher (28.3 ± 1.1) than it was for those cultured in the absence of glutamine or pyruvate (range 20-24). The MCN was the highest (33.4 ± 1.6) for blastocysts cultured in succinate-malate-supplemented HECM-2m followed by those in succinate (28.3 ± 1.1) or malate (24.7 ± 0.5) supplemented HECM-2m. Embryo transfer experiments showed that 29.8% (±4.5) of transferred blastocysts cultured in succinate-malate-supplemented HECM-2m produced live births, similar (P > 0.1) to the control transfers of freshly recovered 8-cells (33.5% ± 2.0) or blastocysts (28.9% ± 3.0). These data show that supplementation of succinate and malate to HECM-2m supports 100% development of hamster 8-cell embryos to high quality viable blastocysts and that non-glucose oxidizable energy substrates are the most preferred components in hamster embryo culture medium. Mol. Reprod. Dev. 47:440-447, 1997. © 1997 Wiley-Liss, Inc.

Calculating the bound spectrum by path summation in action-angle variables

The density of states n(E) is calculated for a bound system whose classical motion is integrable, starting from an expression in terms of the trace of the time-dependent Green function. The novel feature is the use of action-angle variables. This has the advantages that the trace operation reduces to a trivial multiplication and the dependence of n(E) on all classical closed orbits with different topologies appears naturally. The method is contrasted with another, not applicable to integrable systems except in special cases, in which quantization arises from a single closed orbit which is assumed isolated and the trace taken by the method of stationary phase.

Thermodynamic characterization of the reversible, two-state unfolding of maltose binding protein, a large two-domain protein

The folding and stability of maltose binding protein (MBP) have been investigated as a function of pH and temperature by intrinsic tryptophan fluorescence, far- and near-UV circular dichroism, and high-sensitivity differential scanning calorimetric measurements. MBP is a monomeric, two-domain protein containing 370 amino acids. The protein is stable in the pH range of 4-10.5 at 25 degrees C. The protein exhibits reversible, two-state, thermal and guanidine hydrochloride-mediated denaturation at neutral pH. The thermostability of MBP is maximal at pH 6, with a Tm of 64.9 degrees C and a deltaHm of 259.7 kcal mol(-1). The linear dependence of deltaHm on Tm was used to estimate a value of deltaCp of 7.9 kcal mol(-1) K(-1) or 21.3 cal (mol of residue)(-1) K(-1). These values are higher than the corresponding deltaCp's for most globular proteins studied to date. However, the extrapolated values of deltaH and deltaS (per mole of residue) at 110 degrees C are similar to those of other globular proteins. These data have been used to show that the temperature at which a protein undergoes cold denaturation depends primarily on the deltaCp (per mol of residue) and that this temperature increases with an increase in deltaCp. The predicted decrease in stability of MBP at low temperatures was experimentally confirmed by carrying out denaturant-mediated unfolding studies at neutral pH at 2 and 28 degrees C.

On the ecology of cold-water phytoplankton in the Baltic Sea

Increased anthropogenic loading of nitrogen (N) and phosphorus (P) has led to an eutrophication problem in the Baltic Sea, and the spring bloom is a key component in the biological uptake of increased nutrient concentrations. The spring bloom in the Baltic Sea is dominated by both diatoms and dinoflagellates. However, the sedimentation of these groups is different: diatoms tend to sink to the sea floor at the end of the bloom, while dinoflagellates to a large degree are been remineralized in the euphotic zone. Understanding phytoplankton competition and species specific ecological strategies is thus of importance for assessing indirect effects of phytoplankton community composition on eutrophication problems. The main objective of this thesis was to describe some basic physiological and ecological characteristics of the main cold-water diatoms and dinoflagellates in the Baltic Sea. This was achieved by specific studies of: (1) seasonal vertical positioning, (2) dinoflagellate life cycle, (3) mixotrophy, (4) primary production, respiration and growth and (5) diatom silicate uptake, using cultures of common cold-water diatoms: Chaetoceros wighamii, C. gracilis, Pauliella taeniata, Thalassiosira baltica, T. levanderi, Melosira arctica, Diatoma tenuis, Nitzschia frigida, and dinoflagellates: Peridiniella catenata, Woloszynskia halophila and Scrippsiella hangoei. The diatoms had higher primary production capacity and lower respiration rate compared with the dinoflagellates. This difference was reflected in the maximum growth rate, which for the examined diatoms range from 0.6 to 1.2 divisions d-1, compared with 0.2 to 0.3 divisions d-1 for the dinoflagellates. Among diatoms there were species specific differences in light utilization and uptake of silicate, and C. wighamii had the highest carbon assimilation capacity and maximum silicate uptake. The physiological properties of diatoms and dinoflagellates were used in a model of the onset of the spring bloom: for the diatoms the model could predict the initiation of the spring bloom; S. hangoei, on the other hand, could not compete successfully and did not obtain positive growth in the model. The other dinoflagellates did not have higher growth rates or carbon assimilation rates and would thus probably not perform better than S. hangoei in the model. The dinoflagellates do, however, have competitive advantages that were not included in the model: motility and mixotrophy. Previous investigations has revealed that the chain-forming P. catenata performs diurnal vertical migration (DVM), and the results presented here suggest that active positioning in the water column, in addition to DVM, is a key element in this species' life strategy. There was indication of mixotrophy in S. hangoei, as it produced and excreted the enzyme leucine aminopeptidase (LAP). Moreover, there was indirect evidence that W. halophila obtains carbon from other sources than photosynthesis when comparing increase in cell numbers with in situ carbon assimilation rates. The results indicate that mixotrophy is a part of the strategy of vernal dinoflagellates in the Baltic Sea. There were also indications that the seeding of the spring bloom is very important for the dinoflagellates to succeed. In mesocosm experiments dinoflagellates could not compete with diatoms when their initial numbers were low. In conclusion, this thesis has provided new information about the basic physiological and ecological properties of the main cold-water phytoplankton in the Baltic Sea. The main phytoplankton groups, diatoms and dinoflagellates, have different physiological properties, which clearly separate their life strategies. The information presented here could serve as further steps towards better prognostic models of the effects of eutrophication in the Baltic Sea.

Steepest ascent tariff reform

The policy reform literature is primarily concerned with the construction of reforms that yield welfare gains. By contrast, this paper’s contribution is to develop a theoretical concept for which the focus is upon the sizes of welfare gains accruing from policy reforms rather than upon their signs. In undertaking this task, and by focusing on tariff reforms, we introduce the concept of a steepest ascent policy reform, which is a locally optimal reform in the sense that it achieves the highest marginal gain in utility of any feasible local reform. We argue that this reform presents itself as a natural benchmark for the evaluation of the welfare effectiveness of other popular tariff reforms such as the proportional tariff reduction and the concertina rules, since it provides the maximal welfare gain of all possible local reforms. We derive properties of the steepest ascent tariff reform, construct an index to measure the relative welfare effectiveness of any given tariff reform, determine conditions under which proportional and concertina reforms are locally optimal and provide illustrative examples.

Actin Dynamics in Muscle Cells

In every cell, actin is a key component involved in migration, cytokinesis, endocytosis and generation of contraction. In non-muscle cells, actin filaments are very dynamic and regulated by an array of proteins that interact with actin filaments and/or monomeric actin. Interestingly, in non-muscle cells the barbed ends of the filaments are the predominant assembly place, whereas in muscle cells actin dynamics was reported to predominate at the pointed ends of thin filaments. The actin-based thin filament pointed (slow growing) ends extend towards the middle of the sarcomere's M-line where they interact with the thick filaments to generate contraction. The actin filaments in muscle cells are organized into a nearly crystalline array and are believed to be significantly less dynamic than the ones in other cell types. However, the exact mechanisms of the sarcomere assembly and turnover are largely unknown. Interestingly, although sarcomeric actin structures are believed to be relatively non-dynamic, many proteins promoting actin dynamics are expressed also in muscle cells (e.g ADF/cofilin, cyclase-associated protein and twinfilin). Thus, it is possible that the muscle-specific isoforms of these proteins promote actin dynamics differently from their non-muscle counterparts, or that actin filaments in muscle cells are more dynamic than previously thought. To study protein dynamics in live muscle cells, I used primary cell cultures of rat cardiomyocytes. My studies revealed that a subset of actin filaments in cardiomyocyte sarcomeres displays rapid turnover. Importantly, I discovered that the turnover of actin filaments depends on contractility of the cardiomyocytes and that the contractility-induced actin dynamics plays an important role in sarcomere maturation. Together with previous studies those findings suggest that sarcomeres undergo two types of actin dynamics: (1) contractility-dependent turnover of whole filaments and (2) regulatory pointed end monomer exchange to maintain correct thin filament length. Studies involving an actin polymerization inhibitor suggest that the dynamic actin filament pool identified here is composed of filaments that do not contribute to contractility. Additionally, I provided evidence that ADF/cofilins, together with myosin-induced contractility, are required to disassemble non-productive filaments in developing cardiomyocytes. In addition, during these studies we learned that isoforms of actin monomer binding protein twinfilin, Twf-1 and Twf-2a localise to myofibrils in cardiomyocytes and may thus contribute to actin dynamics in myofibrils. Finally, in collaboration with Roberto Dominguez s laboratory we characterized a new actin nucleator in muscle cells - leiomodin (Lmod). Lmod localises towards actin filament pointed ends and its depletion by siRNA leads to severe sarcomere abnormalities in cardiomyocytes. The actin filament nucleation activity of Lmod is enhanced by interactions with tropomyosin. We also revealed that Lmod expression correlates with the maturation of myofibrils, and that it associates with sarcomeres only at relatively late stages of myofibrillogenesis. Thus, Lmod is unlikely to play an important role in myofibril formation, but rather might be involved in the second step of the filament arrangement and/or maintenance through its ability to promote tropomyosin-induced actin filament nucleation occurring at the filament pointed ends. The results of these studies provide valuable new information about the molecular mechanisms underlying muscle sarcomere assembly and turnover. These data offer important clues to understanding certain physiological and pathological behaviours of muscle cells. Better understanding of the processes occurring in muscles might help to find strategies for determining, diagnosis, prognosis and therapy in heart and skeletal muscles diseases.

Acute fluid shifts influence the assessment of serum vitamin D status in critically ill patients

Introduction Recent reports have highlighted the prevalence of vitamin D deficiency and suggested an association with excess mortality in critically ill patients. Serum vitamin D concentrations in these studies were measured following resuscitation. It is unclear whether aggressive fluid resuscitation independently influences serum vitamin D. Methods Nineteen patients undergoing cardiopulmonary bypass were studied. Serum 25(OH)D3, 1α,25(OH)2D3, parathyroid hormone, C-reactive protein (CRP), and ionised calcium were measured at five defined timepoints: T1 - baseline, T2 - 5 minutes after onset of cardiopulmonary bypass (CPB) (time of maximal fluid effect), T3 - on return to the intensive care unit, T4 - 24 hrs after surgery and T5 - 5 days after surgery. Linear mixed models were used to compare measures at T2-T5 with baseline measures. Results Acute fluid loading resulted in a 35% reduction in 25(OH)D3 (59 ± 16 to 38 ± 14 nmol/L, P < 0.0001) and a 45% reduction in 1α,25(OH)2D3 (99 ± 40 to 54 ± 22 pmol/L P < 0.0001) and i(Ca) (P < 0.01), with elevation in parathyroid hormone (P < 0.0001). Serum 25(OH)D3 returned to baseline only at T5 while 1α,25(OH)2D3 demonstrated an overshoot above baseline at T5 (P < 0.0001). There was a delayed rise in CRP at T4 and T5; this was not associated with a reduction in vitamin D levels at these time points. Conclusions Hemodilution significantly lowers serum 25(OH)D3 and 1α,25(OH)2D3, which may take up to 24 hours to resolve. Moreover, delayed overshoot of 1α,25(OH)2D3 needs consideration. We urge caution in interpreting serum vitamin D in critically ill patients in the context of major resuscitation, and would advocate repeating the measurement once the effects of the resuscitation have abated.

Strategies for efficient disruption of metabolism in Mycobacterium tuberculosis from network analysis

Tuberculosis continues to be a major health challenge, warranting the need for newer strategies for therapeutic intervention and newer approaches to discover them. Here, we report the identification of efficient metabolism disruption strategies by analysis of a reactome network. Protein-protein dependencies at a genome scale are derived from the curated metabolic network, from which insights into the nature and extent of inter-protein and inter-pathway dependencies have been obtained. A functional distance matrix and a subsequent nearness index derived from this information, helps in understanding how the influence of a given protein can pervade to the metabolic network. Thus, the nearness index can be viewed as a metabolic disruptability index, which suggests possible strategies for achieving maximal metabolic disruption by inhibition of the least number of proteins. A greedy approach has been used to identify the most influential singleton, and its combination with the other most pervasive proteins to obtain highly influential pairs, triplets and quadruplets. The effect of deletion of these combinations on cellular metabolism has been studied by flux balance analysis. An obvious outcome of this study is a rational identification of drug targets, to efficiently bring down mycobacterial metabolism.