An analysis of a ‘community-driven’ reconstruction of the human metabolic network

Following a strategy similar to that used in baker’s yeast (Herrgård et al. Nat Biotechnol 26:1155–1160, 2008). A consensus yeast metabolic network obtained from a community approach to systems biology (Herrgård et al. 2008; Dobson et al. BMC Syst Biol 4:145, 2010). Further developments towards a genome-scale metabolic model of yeast (Dobson et al. 2010; Heavner et al. BMC Syst Biol 6:55, 2012). Yeast 5—an expanded reconstruction of the Saccharomyces cerevisiae metabolic network (Heavner et al. 2012) and in Salmonella typhimurium (Thiele et al. BMC Syst Biol 5:8, 2011). A community effort towards a knowledge-base and mathematical model of the human pathogen Salmonella typhimurium LT2 (Thiele et al. 2011), a recent paper (Thiele et al. Nat Biotechnol 31:419–425, 2013). A community-driven global reconstruction of human metabolism (Thiele et al. 2013) described a much improved ‘community consensus’ reconstruction of the human metabolic network, called Recon 2, and the authors (that include the present ones) have made it freely available via a database at http://humanmetabolism.org/ and in SBML format at Biomodels (http://identifiers.org/biomodels.db/MODEL1109130000). This short analysis summarises the main findings, and suggests some approaches that will be able to exploit the availability of this model to advantage.     

Daphnias: from the individual based model to the large population equation

The class of deterministic ‘Daphnia’ models treated by Diekmann et al. (J Math Biol 61:277–318, 2010) has a long history going back to Nisbet and Gurney (Theor Pop Biol 23:114–135, 1983) and Diekmann et al. (Nieuw Archief voor Wiskunde 4:82–109, 1984). In this note, we formulate the individual based models (IBM) supposedly underlying those deterministic models. The models treat the interaction between a general size-structured consumer population (‘Daphnia’) and an unstructured resource (‘algae’). The discrete, size and age-structured Daphnia population changes through births and deaths of its individuals and through their aging and growth. The birth and death rates depend on the sizes of the individuals and on the concentration of the algae. The latter is supposed to be a continuous variable with a deterministic dynamics that depends on the Daphnia population. In this model setting we prove that when the Daphnia population is large, the stochastic differential equation describing the IBM can be approximated by the delay equation featured in (Diekmann et al., loc. cit.).     

A radically non-morphemic approach to bidirectional syncretism

This paper addresses the question of how certain kinds of overlapping syncretisms in inflectional paradigms can be accounted for that Baerman et al. (Language 80:807–824, 2005) refer to as convergent/divergent bidirectional syncretisms (based on earlier work by Stump, Inflectional morphology, 2001). Bidirectional syncretism strongly resists accounts in terms of standard rules of exponence (or similar devices) that correlate inflection markers with (often underspecified) morpho-syntactic specifications (such rules are used in many morphological theories; e.g., Anderson, A-morphous morphology, 1992; Halle and Marantz in The view from building, pp. 111–176, 1993; Aronoff, Morphology by itself, 1994; Wunderlich in Yearbook of morphology 1995, pp. 93–114, 1996; and Stump, Inflectional morphology, 2001). The reason is that it is difficult to capture overlapping distributions by natural classes. In view of this, rules of referral have been proposed to derive bidirectional syncretism (Stump, Inflectional morphology, 2001; Baerman et al. (Language 80:807–824, 2005)). In contrast, I would like to pursue the hypothesis that systematic instances of overlapping syncretism ultimately motivate a new approach to inflectional morphology—one that fully dispenses with the assumption that morphological exponents are paired with morpho-syntactic feature specifications (and that therefore qualifies as radically non-morphemic): First, rules of exponence are replaced with feature co-occurrence restrictions (FCRs; Gazdar et al., Generalized Phrase Structure Grammar, 1985). For phonologically determined natural classes of exponents, FCRs state incompatibilites with morpho-syntactic feature specifications. Second, marker competition is resolved by a principle of Phonology-driven Marker Selection (PMS). PMS takes over the role of the Specificity (Blocking, Elsewhere, Panini) Principle of standard analyses. Empirically, the main focus is on Bonan declension; the analysis is subsequently extended to Gujarati conjugation and Latin o-declension, with further remarks on bidirectional syncretism in other inflectional paradigms.     

Modeling Endocrine Control of the Pituitary–Ovarian Axis: Androgenic Influence and Chaotic Dynamics

Mathematical models of the hypothalamus-pituitary-ovarian axis in women were first developed by Schlosser and Selgrade in 1999, with subsequent models of Harris-Clark et al. (Bull. Math. Biol. 65(1):157–173, 2003) and Pasteur and Selgrade (Understanding the dynamics of biological systems: lessons learned from integrative systems biology, Springer, London, pp. 38–58, 2011). These models produce periodic in-silico representation of luteinizing hormone (LH), follicle stimulating hormone (FSH), estradiol (E2), progesterone (P4), inhibin A (InhA), and inhibin B (InhB). Polycystic ovarian syndrome (PCOS), a leading cause of cycle irregularities, is seen as primarily a hyper-androgenic disorder. Therefore, including androgens into the model is necessary to produce simulations relevant to women with PCOS. Because testosterone (T) is the dominant female androgen, we focus our efforts on modeling pituitary feedback and inter-ovarian follicular growth properties as functions of circulating total T levels. Optimized parameters simultaneously simulate LH, FSH, E2, P4, InhA, and InhB levels of Welt et al. (J. Clin. Endocrinol. Metab. 84(1):105–111, 1999) and total T levels of Sinha-Hikim et al. (J. Clin. Endocrinol. Metab. 83(4):1312–1318, 1998). The resulting model is a system of 16 ordinary differential equations, with at least one stable periodic solution. Maciel et al. (J. Clin. Endocrinol. Metab. 89(11):5321–5327, 2004) hypothesized that retarded early follicle growth resulting in “stockpiling” of preantral follicles contributes to PCOS etiology. We present our investigations of this hypothesis and show that varying a follicular growth parameter produces preantral stockpiling and a period-doubling cascade resulting in apparent chaotic menstrual cycle behavior. The new model may allow investigators to study possible interventions returning acyclic patients to regular cycles and guide developments of individualized treatments for PCOS patients.     

Establishment and characterization of METON myoepithelioma cell line derived from human palatal myoepithelioma: apical reference to the diverse differentiation potential

Myoepithelioma is an extremely rare condition that accounts for 1–1.5 % of salivary gland tumors. It was formerly regarded as a subtype of pleomorphic adenoma, in which myoepithelial structural components predominated, but was listed as a separate disease entity in the 1991 World Health Organization classification (Seifert in Histological typing of salivary gland tumours. Springer, Berlin, 1991). Its histology is highly varied and recurrence is frequent (El-Naggar et al. in J Larygol Otol 103:1192–1197, 1989), with cases of malignant transformation having been reported (Seifert in Histological typing of salivary gland tumours. Springer, Berlin, 1991; Barnes et al. in Pathology and Genetics of head and neck tumours. IARC Press, Lyon, 2005), making this a difficult tumor to control in many cases. This is thought to be due to the multiple differentiation potential of myoepithelial cells, but the details are unknown. There have been a number of reports of the establishment of cell lines (Shirasuna et al. Cancer. 45:297–305, 1980; Jaeger et al. Oral Surg Oral Med Oral Pathol Oral Radiol Endod 84:663–667, 1997), but numerous points remain unclear. We established a myoepithelial cell line designated METON, and investigated its characteristics. METON consists of cells with two different morphologies: spindle-shaped cells and epithelial-like cells. Then. we also used single-cell cloning method to establish various subclones (epithelial-like, spindle-like, and mixed epithelial-like/spindle-like cell lines). Among these, pluripotency markers were expressed by the mixed epithelial-like/spindle-like cell lines. The newly established cell line expressing these pluripotency markers will be extremely useful for elucidating the diverse histologies of salivary gland tumors.     

The (Paleo)Geography of Evolution: Making Sense of Changing Biology and Changing Continents

During the voyage of the H.M.S. Beagle, Charles Darwin quickly realized that geographic isolation led to significant changes in the adaptation of local flora and fauna (Darwin 1859). Genetic isolation is one of the well-known mechanisms by which adaptation (allopatric speciation) can occur (Palumbi, Annu Rev Ecol Syst 25:547?C72, 1994; Ricklefs, J Avian Biol 33:207?C11, 2002; Burns et al., Evolution 56:1240?C52, 2002; Hendry et al., Science 290:516?C8, 2009). Evolutionary changes can also occur when landmasses converge or are ??bridged.?? An important and relatively recent (Pliocene Epoch) example known as the ??Great American Biotic Interchange?? allowed for the migration of previously isolated species into new ecological niches between North and South America (Webb 1985, Ann Mo Bot Gard 93:245?C57, 2006; Kirby and MacFadden, Palaeogeogr Palaeoclimatol Palaeoecol 228:193?C202, 2005). Geographic isolation (vicariance) or geographic merging (geodispersal) can occur for a variety of reasons (sea level rise, splitting of continents, mountain building). In addition, the growth of a large supercontinent (or breakup) may change the climatic zonation on the globe and form a different type of barrier for species migration. This short review paper focuses on changing paleogeography throughout the Phanerozoic and the close ties between paleogeography and the evolutionary history of life on Earth.     

Simple signaling games of sexual selection (Grafen’s revisited)

We investigate several versions of a simple game of sexual selection, to explore the role of secondary sexual characters (the “handicap paradox”) with the tools of signaling theory. Our models admit closed form solutions. They are very much inspired by Grafen’s (J Theor Biol 144:517–546, 1990a; J Theor Biol 144:473–516, 1990b) seminal companion papers. By merging and simplifying his two approaches, we identify a not so minor artifact in the seminal study. We propose an alternative model to start with Grafen’s sexual selection theory, with several similarities with Getty (Anim Behav 56:127–130, 1998).     

Phenomenology of Psychoanalytic Data. A Biosemiotic Framework

In my continuing efforts to build a bridge between psychoanalytic findings and biosemiotics here, as in previous works, ‘biosemiotic’ refers to the hierarchy of meaning-forms (from biological to semiotic-organizations) underlying an updated psychoanalytic model of mind. Within this framework I present a broad range of bio-semiotic phenomena, processes, dynamics, defenses, and universal and unique internalized interpersonal patterns, that in psychoanalysis all commonly fall under the broad heading of the “Unconscious.” Reconceptualized as interpretive data within the purview of a psychoanalytic discourse-semantic this biosemiotic framework posits an epigenetic continuum of human meaning-organizations originating at basic organic levels, moving upward through biological, psycho-somatic and affective expression, proto-semiotic transmissions, represented forms, and finally to explicit linguistic signs and complex symbol systems. In addition to assuming an uninterrupted epigenetic continuum crystallizing in hierarchic organization, this framework accentuates the multilayered and increasingly condensed quality of higher more elaborate organizations of meaning in human communication, drawing attention to persisting biological undercurrents in implied sense, intent, and motivation, all of which impact on repressive/defensive mechanisms. Drawing from previous works (Aragno 1997, 2005, 2008a, b, Psychoanalytic Inquiry 29(1):30–47, 2009, Biosemiotics 3:57–77, 2010, 2011a, Signs 5:71–74, 2011b, Journal of the American Psychoanalytic Association, Centennial Paper, Special Centennial Issue 59(2):239–288, 2011b, Signs 5:29–70, 2011c) in which I labored to update and revise Freud’s first topographical theory of mind, this paper presents the phenomenology of unconscious ‘data’ for the purpose of introducing a diverse range of non-linguistic signifying forms from which psychoanalysts infer mental processes and ‘interpret’ meanings. An important underlying premise regarding psychoanalytic data and its relation to the basic biosemiotic ‘agenda’ is that until grounded in an updated developmental theory of mind inclusive of pre- and proto-semiotic-forms, that is evolutionarily plausible, epistemologically based, and correlates with contemporary neuroscience, the term “sign” is merely an abstract linguistic ‘label’ rather than a mental act with antecedent developmental stages manifesting meanings through different forms and modes of expression. Drawn from the yields of the psychoanalytic method and semantic this revised metatheoretical approach provides insights into the sensory-emotive, bodily origins of unconscious layers of non-linguistic signification thereby expanding our understanding of the formative stages of the ‘semiotic function’ in human evolution. This being the third in a series of papers integrating the yields of psychoanalytic methodology with the underlying premises of ‘Biosemiotics,’ some familiarity with the background knowledge provided in the previous two is strongly recommended.     

Feature extraction from spike trains with Bayesian binning: ‘Latency is where the signal starts’

The peristimulus time histogram (PSTH) and its more continuous cousin, the spike density function (SDF) are staples in the analytic toolkit of neurophysiologists. The former is usually obtained by binning spike trains, whereas the standard method for the latter is smoothing with a Gaussian kernel. Selection of a bin width or a kernel size is often done in an relatively arbitrary fashion, even though there have been recent attempts to remedy this situation (DiMatteo et al., Biometrika 88(4):1055–1071, 2001; Shimazaki and Shinomoto 2007a, Neural Comput 19(6):1503–1527, 2007b, c; Cunningham et al. 2008). We develop an exact Bayesian, generative model approach to estimating PSTHs. Advantages of our scheme include automatic complexity control and error bars on its predictions. We show how to perform feature extraction on spike trains in a principled way, exemplified through latency and firing rate posterior distribution evaluations on repeated and single trial data. We also demonstrate using both simulated and real neuronal data that our approach provides a more accurate estimates of the PSTH and the latency than current competing methods. We employ the posterior distributions for an information theoretic analysis of the neural code comprised of latency and firing rate of neurons in high-level visual area STSa. A software implementation of our method is available at the machine learning open source software repository (www.mloss.org, project ‘binsdfc’).     

On the ‘transmission sense of information’

In order to illuminate the role of information in biology, Bergstrom and Rosvall (Biol Philos 26:159–176, 2011a; Biol Philos 26:195–200, 2011b) propose a ‘transmission sense of information’ which builds on Shannon’s theory. At the core of the transmission sense is an appeal to the reduction in uncertainty in receivers and to etiological function. I explore several ways of cashing out uncertainty reduction as well as the consequences of appealing to function.     

Mathematical Model of Growth Factor Driven Haptotaxis and Proliferation in a Tissue Engineering Scaffold

Motivated by experimental work (Miller et al. in Biomaterials 27(10):2213–2221, 2006, 32(11):2775–2785, 2011) we investigate the effect of growth factor driven haptotaxis and proliferation in a perfusion tissue engineering bioreactor, in which nutrient-rich culture medium is perfused through a 2D porous scaffold impregnated with growth factor and seeded with cells. We model these processes on the timescale of cell proliferation, which typically is of the order of days. While a quantitative representation of these phenomena requires more experimental data than is yet available, qualitative agreement with preliminary experimental studies (Miller et al. in Biomaterials 27(10):2213–2221, 2006) is obtained, and appears promising. The ultimate goal of such modeling is to ascertain initial conditions (growth factor distribution, initial cell seeding, etc.) that will lead to a final desired outcome.     

Allocation in life cycle inventory: partial set of solutions to an ill-posed problem

Background, aim, and scope

The primary aim of this paper is to indicate that partitioning allocation methods yields only a small subset of solutions to an ill-posed problem that has potentially infinitely many exact solutions. It will be shown that each of the existing partitioning methods arrives at just one particular solution from among infinitely many solutions of an underdetermined system of linear equations.

Materials and methods

Some life cycle inventories fall into a class of functions called estimable functions in linear model framework, in which case they are invariant to allocation assumptions. This class of functions unites results described by Heijungs and Frischknecht (Int J Life Cycle Assess 3:321–332, 1998) and Heijungs and Suh (2002, Conjecture 1, p. 91). The inventories for non-estimable functions obtained through allocation are, in fact, derived under a set of additional implicit equality constraints called side conditions, often resulting in inventory results which differ greatly from one allocation to the next.

Results and discussions

This paper explicates (1) identification of all estimable functions from any given technology matrix and (2) recovery of side conditions imposed on non-estimable functions through partitioning. These methods are illustrated in a simple example, and their relation to least squares techniques for allocation explored by Marvuglia et al. (Int J Life Cycle Assess 15:1020–1040, 2010) ;(Int J Agr Environ Inf Syst 3:51–71, 2012) are discussed.

Conclusions and outlook

Recommendations are made that may lead to more meaningful ways to obtain additional data or include additional information in life cycle inventories in the future.     

Mitochondrial complex I inhibitors, acetogenins, induce HepG2 cell death through the induction of the complete apoptotic mitochondrial pathway

The development of new anti-neoplastic drugs is a key issue for cancer chemotherapy due to the reality that, most likely, certain cancer cells are resistant to current chemotherapy. The past two decades have witnessed tremendous advances in our understanding of the pathogenesis of cancer. These advances have allowed identification new targets as oncogenes, tumor supressor genes and the possible implication of the mitochondria (Fulda et al. Nat Rev Drug Discov 9:447–464, 2010). Annonaceous Acetogenins (ACGs) have been described as the most potent inhibitors of the respiratory chain because of their interaction with mitochondrial Complex I (Degli Esposti and Ghelli Biochim Biophys Acta 1187:116–120, 1994; Zafra-Polo et al. Phytochemistry 42:253–271, 1996; Miyoshi et al. Biochim Biophys Acta 1365:443–452, 1998; Tormo et al. Arch Biochem Biophys 369:119–126, 1999; Motoyama et al. Bioorg Med Chem Lett 12:2089–2092, 2002). To explore a possible application of natural products from Annonaceous plants to cancer treatment, we have selected four bis-tetrahydrofuranic ACGs, three from Annona cherimolia (cherimolin-1, motrilin and laherradurin) and one from Rollinia mucosa (rollinianstatin-1) in order to fully describe their mechanisms responsible within the cell (Fig. 1). In this study, using a hepato-carcinoma cell line (HepG2) as a model, we showed that the bis-THF ACGs caused cell death through the induction of the apoptotic mitochondrial pathway. Their potency and behavior were compared with the classical mitochondrial respiratory chain Complex I inhibitor rotenone in every apoptotic pathway step.

Metalloproteins and the Pyrite-based Origin of Life: A Critical Assessment

We critically examine the proposal by Wächtershäuser (Prokaryotes 1:275?C283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787?C1808, 2006b) that putative transition metal binding sites in protein components of the translation machinery of hyperthermophiles provide evidence of a direct relationship with the FeS clusters of pyrite and thus indicate an autotrophic origin of life in volcanic environments. Analysis of completely sequenced cellular genomes of Bacteria, Archaea and Eucarya does not support the suggestion by Wächtershäuser (Prokaryotes 1:275?C283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787?C1808, 2006b) that aminoacyl-tRNA synthetases and ribosomal proteins bear sequence signatures typical of strong covalent metal bonding whose absence in mesophilic species reveals a process of adaptation towards less extreme environments.     

A Mathematical and Numerical Investigation of the Hemodynamical Origins of Oscillations in Microvascular Networks

Evidence is presented to show that self-sustained oscillations of purely hemodynamical origin are possible in some arcade-type microvascular networks supplied with steady boundary conditions, but that in others the oscillations disappear with sufficient reduction of the time step Δt, showing them to be numerical artefacts. In an attempt to elucidate the mechanisms involved in the onset of fluctuations, we proceed to perform a linear stability analysis for the convective model of Kiani et al. (Microvasc. Res. 45:219–232, 1993; Am. J. Physiol. 266(35):H1822–H1828, 1994), and show that this leads via a system of delay differential equations to a nonlinear eigenvalue problem. This result generalises the characteristic equation obtained by Carr et al. (Ann. Biomed. Eng. 33:764–771, 2005) and Geddes et al. (SIAM J. Appl. Dyn. Syst. 6(4):694–727, 2007) who solved a special case in a two node network. An implicit numerical method is proposed for the computation of blood flows in networks using the convective model. In a moderate size subnetwork of one of the networks chosen by Kiani et al. (Am. J. Physiol. 266(35):H1822–H1828, 1994), the topology, vessel lengths, and diameters of which were based on microvascular networks in the rat mesentery, we compare results generated using the original explicit numerical method of Kiani et al. (Am. J. Physiol. 266(35):H1822–H1828, 1994) with those from our implicit scheme. From the linear stability theory, a critical value D _RBC,crit of a red blood cell diameter parameter D _RBC in the plasma skimming model of Fenton et al. (Pflügers Arch. 403:396–401, 1985b) is identified for the onset of oscillations about steady state and both the explicit and implicit methods are used to calculate the inflow hematocrit solutions in all vessels of the subnetwork at the critical parameter value, subject to perturbed initial conditions. The results of the implicit method are demonstrated to be in excellent and superior agreement with the predictions of the linear analysis in this case. For values of D _RBC slightly larger than D _RBC,crit the bifurcating periodic solutions calculated using either the explicit or implicit schemes are characteristic of those of a supercritical Hopf bifurcation and the graphs of D _RBC vs. oscillation amplitude would seem to converge as Δt→0.     

Leges sine moribus vanae: does language make moral thinking possible?

Does language make moral cognition possible? Some authors like Andy Clark have argued for a positive answer whereby language and the ways people use it mark a fundamental divide between humans and all other animals with respect to moral thinking (Clark, Mind and morals: essays on cognitive science and ethics. MIT Press, Cambridge, MA, 1996; Moral Epistemol Nat Can J Philos Suppl XXVI, 2000a; Moral Epistemol Nat Can J Philos Suppl XXVI, 2000b; Philosophy of mental representation. Oxford University Press, Oxford, pp 37–43 and discussion, 44–61, 2002). I take issue with Clark’s view and argue that language is probably unnecessary for the emergence of moral cognition. I acknowledge, however, that humans unlike other animals seem to posses what Haugeland in Philosophy of mental representation. Oxford University Press, Oxford (2002) terms ‘norm-hungriness’: an idiosyncratic need or desire to create and abide by a multitude of norms. Our peculiar norm-hungriness, I suggest, depends on what can be called florid control rather than on language.     

Taxonomic review of Kyphosus (Pisces: Kyphosidae) in the Atlantic and Eastern Pacific Oceans

The taxonomy of the Atlantic and Eastern Pacific species of Kyphosus is reviewed with K. bosquii (Lacepède 1802), K. incisor (Cuvier 1831), K. analogus (Gill 1862) and K. elegans (Peters 1869) considered valid, and K. atlanticus sp. nov. newly described. Kyphosus bosquii and K. atlanticus are both characterized by 12 dorsal- and 11 anal-fin soft rays, but differ in the number of longitudinal scale rows along the midbody (61–66, mode 63 vs. 50–56, mode 54). Kyphosus incisor and K. analogus, characterized by 14 dorsal- and 13 anal-fin soft rays, similarly differ from each other in midbody longitudinal scale row counts (57–64, mode 60 vs. 68–74, mode 70 or 72). Kyphosus elegans is characterized by 13 dorsal- and 12 anal-fin soft rays, and 51–57 midbody longitudinal scale rows. Kyphosus bosquii, K. atlanticus and K. incisor are distributed in the Atlantic Ocean, K. analogus and K. elegans occurring in the Eastern Pacific. The holotype of Pimelepterus flavolineatus Poey 1866, here regarded as a junior synonym of K. incisor, was located within a collection of Cuban fishes donated to the Smithsonian Institution by Poey in 1873. A neotype is designated here for K. analogus. Pimelepterus gallveii Cunningham 1910, Kyphosus palpebrosus Miranda-Ribeiro 1919 and K. metzelaari Jordan and Evermann 1927 are recognized as junior synonyms of K. bosquii. Pimelepterus sandwicensis Sauvage 1880 is a junior synonym of K. elegans. Perca saltatrix Linnaeus 1758, together with the replacement name Perca sectatrix Linnaeus 1766, is regarded as nomina dubia.