Ecosystem biogeochemistry considered as a distributed metabolic network ordered by maximum entropy production

Author Posting. © The Author(s), 2009. This is the author's version of the work. It is posted here by permission of The Royal Society for personal use, not for redistribution. The definitive version was published in Philosophical Transactions of the Royal Society B 365 (2010): 1417-1427, doi:10.1098/rstb.2009.0272.We examine the application of the maximum entropy production principle for describing ecosystem biogeochemistry. Since ecosystems can be functionally stable despite changes in species composition, we utilize a distributed metabolic network for describing biogeochemistry, which synthesizes generic biological structures that catalyze reaction pathways, but is otherwise organism independent. Allocation of biological structure and regulation of biogeochemical reactions is determined via solution of an optimal control problem in which entropy production is maximized. However, because synthesis of biological structures cannot occur if entropy production is maximized instantaneously, we propose that information stored within the metagenome allows biological systems to maximize entropy production when averaged over time. This differs from abiotic systems that maximize entropy production at a point in space-time, which we refer to as the steepest descent pathway. It is the spatiotemporal averaging that allows biological systems to outperform abiotic processes in entropy production, at least in many situations. A simulation of a methanotrophic system is used to demonstrate the approach. We conclude with a brief discussion on the implications of viewing ecosystems as self organizing molecular machines that function to maximize entropy production at the ecosystem level of organization.The work presented here was funded by the PIE-LTER program (NSF OCE-0423565), as well as from NSF CBET-0756562, NSF EF-0928742 and NASA Exobiology and Evolutionary Biology (NNG05GN61G)

eQuilibrator—the biochemical thermodynamics calculator

The laws of thermodynamics constrain the action of biochemical systems. However, thermodynamic data on biochemical compounds can be difficult to find and is cumbersome to perform calculations with manually. Even simple thermodynamic questions like ‘how much Gibbs energy is released by ATP hydrolysis at pH 5?’ are complicated excessively by the search for accurate data. To address this problem, eQuilibrator couples a comprehensive and accurate database of thermodynamic properties of biochemical compounds and reactions with a simple and powerful online search and calculation interface. The web interface to eQuilibrator (http://equilibrator.weizmann.ac.il) enables easy calculation of Gibbs energies of compounds and reactions given arbitrary pH, ionic strength and metabolite concentrations. The eQuilibrator code is open-source and all thermodynamic source data are freely downloadable in standard formats. Here we describe the database characteristics and implementation and demonstrate its use

anNET: a tool for network-embedded thermodynamic analysis of quantitative metabolome data

Background: Compared to other omics techniques, quantitative metabolomics is still at its infancy. Complex sample preparation and analytical procedures render exact quantification extremely difficult. Furthermore, not only the actual measurement but also the subsequent interpretation of quantitative metabolome data to obtain mechanistic insights is still lacking behind the current expectations. Recently, the method of network-embedded thermodynamic (NET) analysis was introduced to address some of these open issues. Building upon principles of thermodynamics, this method allows for a quality check of measured metabolite concentrations and enables to spot metabolic reactions where active regulation potentially controls metabolic flux. So far, however, widespread application of NET analysis in metabolomics labs was hindered by the absence of suitable software. Results: We have developed in Matlab a generalized software called 'anNET' that affords a user-friendly implementation of the NET analysis algorithm. anNET supports the analysis of any metabolic network for which a stoichiometric model can be compiled. The model size can span from a single reaction to a complete genome-wide network reconstruction including compartments. anNET can (i) test quantitative data sets for thermodynamic consistency, (ii) predict metabolite concentrations beyond the actually measured data, (iii) identify putative sites of active regulation in the metabolic reaction network, and (iv) help in localizing errors in data sets that were found to be thermodynamically infeasible. We demonstrate the application of anNET with three published Escherichia coli metabolome data sets. Conclusion: Our user-friendly and generalized implementation of the NET analysis method in the software anNET allows users to rapidly integrate quantitative metabolome data obtained from virtually any organism. We envision that use of anNET in labs working on quantitative metabolomics will provide the systems biology and metabolic engineering communities with a mean to proof the quality of metabolome data sets and with all further benefits of the NET analysis approach.

Multiplicity Distributions and Rapidity Gaps

I examine the phenomenology of particle multiplicity distributions, with special emphasis on the low multiplicities that are a background in the study of rapidity gaps. In particular, I analyze the multiplicity distribution in a rapidity interval between two jets, using the HERWIG QCD simulation with some necessary modifications. The distribution is not of the negative binomial form, and displays an anomalous enhancement at zero multiplicity. Some useful mathematical tools for working with multiplicity distributions are presented. It is demonstrated that ignoring particles with pt<0.2 has theoretical advantages, in addition to being convenient experimentally.Comment: 24 pages, LaTeX, MSUHEP/94071

Oncocytic carcinoma of parotid gland: a case report with clinical, immunohistochemical and ultrastructural features

BACKGROUND: Oncocytic carcinoma is an extremely rare neoplasm of the salivary glands. We report a case of oncocytic carcinoma arising in a parotid gland in a 66-year-old female. METHOD: An excisional biopsy of the parotid tumor was performed. The specimen was submitted for histology and after fixation in formalin solution and inclusion in paraffin, 3–5 μm sections were stained with hematoxylin and eosin for conventional evaluation and Periodic acid Schiff stain. Immunohistochemical studies were performed using antibodies against mitochondrial antigen, keratin, S-100, alpha-actin, vimentin, alpha-1-antichymotrypsin as well as an ultrastructural analysis was performed. RESULTS: Frozen sections revealed an infiltrative growth pattern and the diagnosis of a malignant epithelial lesion was made. Permanent sections stained with haematoxylin and eosin revealed a neoplasm that had replaced a wide area of the parotid gland and had invaded subcutaneous adipose tissue. Perineural invasion was evident, but vascular invasion was not found. Neoplastic elements were large, round or polyhedral cells and were arranged in solid sheets, islands and cords. The cytoplasm was abundant, eosinophilic and finely granular. The nuclei were large and located centrally or peripherally. The nucleoli were distinct and large. Periodic acid Schiff stain demonstrated a granular cytoplasm. Immunohistochemistry demonstrated mithochondrial antigen, keratin, and chymotrypsin immunoreactivity in the neoplastic cells. Ultrastructural analysis revealed numerous mitochondria packed into the cytoplasm of the neoplastic cells. Thus, the final diagnosis was that of oncocytic carcinoma of parotid gland. CONCLUSION: This neoplasm shows clinical, microscopical, histological and ultrastructural features of oncocytic carcinoma and this must be considered in the differential diagnosis of other proliferations in the parotid gland with abundant granular cytoplasm and metastatic oncocytic carcinomas

Conservation laws and work fluctuation relations in chemical reaction networks

We formulate a nonequilibrium thermodynamic description for open chemical reaction networks (CRN) described by a chemical master equation. The topological properties of the CRN and its conservation laws are shown to play a crucial role. They are used to decompose the entropy production into a potential change and two work contributions, the first due to time dependent changes in the externally controlled chemostats concentrations and the second due to flows maintained across the system by nonconservative forces. These two works jointly satisfy a Jarzynski and Crooks fluctuation theorem. In absence of work, the potential is minimized by the dynamics as the system relaxes to equilibrium and its equilibrium value coincides with the maximum entropy principle. A generalized Landauer's principle also holds: the minimal work needed to create a nonequilibrium state is the relative entropy of that state to its equilibrium value reached in absence of any work.Comment: revtex format: 30 pages (25 + 5 for appendices), 9 figures, 3 tables. v2: published versio

Numerical approach to a model for quasistatic damage with spatial BV-regularization

We address a model for rate-independent, partial, isotropic damage in quasistatic small strain linear elasticity, featuring a damage variable with spatial BV-regularization. Discrete solutions are obtained using an alternate time-discrete scheme and the Variable-ADMM algorithm to solve the constrained nonsmooth optimization problem that determines the damage variable at each time step. We prove convergence of the method and show that discrete solutions approximate a semistable energetic solution of the rate-independent system. Moreover, we present our numerical results for two benchmark problems

Including metabolite concentrations into flux balance analysis: thermodynamic realizability as a constraint on flux distributions in metabolic networks

<p>Abstract</p> <p>Background</p> <p>In recent years, constrained optimization – usually referred to as flux balance analysis (FBA) – has become a widely applied method for the computation of stationary fluxes in large-scale metabolic networks. The striking advantage of FBA as compared to kinetic modeling is that it basically requires only knowledge of the stoichiometry of the network. On the other hand, results of FBA are to a large degree hypothetical because the method relies on plausible but hardly provable optimality principles that are thought to govern metabolic flux distributions.</p> <p>Results</p> <p>To augment the reliability of FBA-based flux calculations we propose an additional side constraint which assures thermodynamic realizability, i.e. that the flux directions are consistent with the corresponding changes of Gibb's free energies. The latter depend on metabolite levels for which plausible ranges can be inferred from experimental data. Computationally, our method results in the solution of a mixed integer linear optimization problem with quadratic scoring function. An optimal flux distribution together with a metabolite profile is determined which assures thermodynamic realizability with minimal deviations of metabolite levels from their expected values. We applied our novel approach to two exemplary metabolic networks of different complexity, the metabolic core network of erythrocytes (30 reactions) and the metabolic network iJR904 of <it>Escherichia coli </it>(931 reactions). Our calculations show that increasing network complexity entails increasing sensitivity of predicted flux distributions to variations of standard Gibb's free energy changes and metabolite concentration ranges. We demonstrate the usefulness of our method for assessing critical concentrations of external metabolites preventing attainment of a metabolic steady state.</p> <p>Conclusion</p> <p>Our method incorporates the thermodynamic link between flux directions and metabolite concentrations into a practical computational algorithm. The weakness of conventional FBA to rely on intuitive assumptions about the reversibility of biochemical reactions is overcome. This enables the computation of reliable flux distributions even under extreme conditions of the network (e.g. enzyme inhibition, depletion of substrates or accumulation of end products) where metabolite concentrations may be drastically altered.</p

Abiotic ammonium formation in the presence of Ni-Fe metals and alloys and its implications for the Hadean nitrogen cycle

Experiments with dinitrogen-, nitrite-, nitrate-containing solutions were conducted without headspace in Ti reactors (200°C), borosilicate septum bottles (70°C) and HDPE tubes (22°C) in the presence of Fe and Ni metal, awaruite (Ni80Fe20) and tetrataenite (Ni50Fe50). In general, metals used in this investigation were more reactive than alloys toward all investigated nitrogen species. Nitrite and nitrate were converted to ammonium more rapidly than dinitrogen, and the reduction process had a strong temperature dependence. We concluded from our experimental observations that Hadean submarine hydrothermal systems could have supplied significant quantities of ammonium for reactions that are generally associated with prebiotic synthesis, especially in localized environments. Several natural meteorites (octahedrites) were found to contain up to 22 ppm Ntot. While the oxidation state of N in the octahedrites was not determined, XPS analysis of metals and alloys used in the study shows that N is likely present as nitride (N3-). This observation may have implications toward the Hadean environment, since, terrestrial (e.g., oceanic) ammonium production may have been supplemented by reduced nitrogen delivered by metal-rich meteorites. This notion is based on the fact that nitrogen dissolves into metallic melts