Master Equation: Biological Applications and Thermodynamic Description

It is well known that many realistic mathematical models of biological systems, such as cell growth, cellular development and differentiation, gene expression, gene regulatory networks, enzyme cascades, synaptic plasticity, aging and population growth need to include stochasticity. These systems are not isolated, but rather subject to intrinsic and extrinsic fluctuations, which leads to a quasi equilibrium state (homeostasis). The natural framework is provided by Markov processes and the Master equation (ME) describes the temporal evolution of the probability of each state, specified by the number of units of each species. The ME is a relevant tool for modeling realistic biological systems and allow also to explore the behavior of open systems. These systems may exhibit not only the classical thermodynamic equilibrium states but also the nonequilibrium steady states (NESS). This thesis deals with biological problems that can be treat with the Master equation and also with its thermodynamic consequences. It is organized into six chapters with four new scientific works, which are grouped in two parts: (1) Biological applications of the Master equation: deals with the stochastic properties of a toggle switch, involving a protein compound and a miRNA cluster, known to control the eukaryotic cell cycle and possibly involved in oncogenesis and with the propose of a one parameter family of master equations for the evolution of a population having the logistic equation as mean field limit. (2) Nonequilibrium thermodynamics in terms of the Master equation: where we study the dynamical role of chemical fluxes that characterize the NESS of a chemical network and we propose a one parameter parametrization of BCM learning, that was originally proposed to describe plasticity processes, to study the differences between systems in DB and NESS.

Thermodynamic description and phase transformation of highly symmetrical monoatomic structures

Based on theoretical considerations an explanation for the temperature dependence of the thermal expansion and the bulk modulus is proposed. A new equation state is also derived. Additionally a physical explanation for the latent heat of fusion is presented. These theoretical predictions are tested against experiments on highly symmetrical monatomic structures. ^ The volume is not an independent variable and must be broken down into its fundamental components when the relationships to the pressure and temperature are defined. Using zero pressure and temperature reference frame, the initial parameters, volume at zero pressure and temperature[V°], bulk modulus at zero temperature [K°] and volume coefficient of thermal expansion at zero pressure[α°] are defined. ^ The new derived EoS is tested against the experiments on perovskite and epsilon iron. The Root-mean-square-deviations (RMSD) of the residuals of the molar volume, pressure, and temperature are in the range of the uncertainty of the experiments. ^ Separating the experiments into 200 K ranges, the new EoS was compared to the most widely used finite strain, interatomic potential, and empirical isothermal EoSs such as the Burch-Murnaghan, the Vinet, and the Roy-Roy respectively. Correlation coefficients, RMSD's of the residuals, and Akaike Information Criteria were used for evaluating the fitting. Based on these fitting parameters, the new p-V-T EoS is superior in every temperature range relative to the investigated conventional isothermal EoS. ^ The new EoS for epsilon iron reproduces the preliminary-reference earth-model (PREM) densities at 6100-7400 K indicating that the presence of light elements might not be necessary to explain the Earth's inner core densities. ^ It is suggested that the latent heat of fusion supplies the energy required for overcoming on the viscous drag resistance of the atoms. The calculated energies for melts formed from highly symmetrical packing arrangements correlate very well with experimentally determined latent heat values. ^ The optical investigation of carhonado-diamond is also part of the dissertation. The collected first complete infrared FTIR absorption spectra for carhonado-diamond confirm the interstellar origin for the most enigmatic diamonds known as carbonado. ^

Vapour composition and activities in Mg-Zn liquid alloy at 923 K

Vapour species effusing from a magnesia Knudsen cell containing Mg-Zn alloy at 923 K were condensed on a water cooled copper plate. The equilibrium composition of the vapour phase over the alloy was determined from chemical analysis of the condensate. The activity coefficients of both components in the alloy have been derived from the data using a modified Gibbs-Duhem relation. The ratio of saturation vapour pressures of pure Zn and Mg obtained from the analysis of alloy data agree well with values from the literature, providing an internal check on the accuracy of data obtained in this study. Both components of the alloy exhibit negative deviations from Raoult's law. The concentration-concentration structure factor of Bhatia and Thomton at zero wave vector, evaluated from the measurements, indicate the presence of MgZn2 type complex in the liquid state. The associated regular solution model has been used for the thermodynamic description of liquid Mg-Zn alloys.

Phase Diagram for the System RuO2-TiO2 in Air

There are conflicting reports in the literature regarding solid solubility in the system RuO2-TiO2. To resolve this issue a few experiments were conducted in air at 1673, 1723, and 1773 K. The results show limited terminal solid solubility. There is an extended solid-state miscibility gap that intersects the decomposition curve for the RuO2-rich solid solution generating a peritectoid reaction at 1698 K. The measured equilibrium compositions of the solid solutions are used to develop a thermodynamic description of the oxide solid solution with rutile structure. Using the subregular solution model, the enthalpy of mixing can be represented by the expression, Delta H-M/J center dot mol(-1) = XTiO2XRuO2 ( 34,100X(TiO2) + 30,750X(RuO2)). The binodal and spinodal curves and T-X phase diagram in air are computed using this datum and Gibbs energy of formation of RuO2 available in the literature. The computed results suggest that equilibrium was not attained during solubility measurements at lower temperatures reported in the literature.

Water in silicate glasses

The speciation of water in a variety of hydrous silicate glasses, including simple and rhyolitic compositions, synthesized over a range of experimental conditions with up to 11 weight percent water has been determined using infrared spectroscopy. This technique has been calibrated with a series of standard glasses and provides a precise and accurate method for determining the concentrations of molecular water and hydroxyl groups in these glasses.

For all the compositions studied, most of the water is dissolved as hydroxyl groups at total water contents less than 3-4 weight percent; at higher total water contents, molecular water becomes the dominant species. For total water contents above 3-4 weight percent, the amount of water dissolved as hydroxyl groups is approximately constant at about 2 weight percent and additional water is incorporated as molecular water. Although there are small but measurable differences in the ratio of molecular water to hydroxyl groups at a given total water content among these silicate glasses, the speciation of water is similar over this range of composition. The trends in the concentrations of the H-bearing species in the hydrous glasses included in this study are similar to those observed in other silicate glasses using either infrared or NMR spectroscopy.

The effects of pressure and temperature on the speciation of water in albitic glasses have been investigated. The ratio of molecular water to hydroxyl groups at a given total water content is independent of the pressure and temperature of equilibration for albitic glasses synthesized in rapidly quenching piston cylinder apparatus at temperatures greater than 1000°C and pressures greater than 8 kbar. For hydrous glasses quenched from melts cooled at slower rates (i.e., in internally heated or in air-quench cold seal pressure vessels), there is an increase in the ratio of molecular water to hydroxyl group content that probably reflects reequilibration of the melt to lower temperatures during slow cooling.

Molecular water and hydroxyl group concentrations in glasses provide information on the dissolution mechanisms of water in silicate liquids. Several mixing models involving homogeneous equilibria of the form H_2O + O = 20H among melt species have been explored for albitic melts. These models can account for the measured species concentrations if the effects of non-ideal behavior or mixing of polymerized units are included, or by allowing for the presence of several different types of anhydrous species.

A thermodynamic model for hydrous albitic melts has been developed based on the assumption that the activity of water in the melt is equal to the mole fraction of molecular water determined by infrared spectroscopy. This model can account for the position of the watersaturated solidus of crystalline albite, the pressure and temperature dependence of the solubility of water in albitic melt, and the volumes of hydrous albitic melts. To the extent that it is successful, this approach provides a direct link between measured species concentrations in hydrous albitic glasses and the macroscopic thermodynamic properties of the albite-water system.

The approach taken in modelling the thermodynamics of hydrous albitic melts has been generalized to other silicate compositions. Spectroscopic measurements of species concentrations in rhyolitic and simple silicate glasses quenched from melts equilibrated with water vapor provide important constraints on the thermodynamic properties of these melt-water systems. In particular, the assumption that the activity of water is equal to the mole fraction of molecular water has been tested in detail and shown to be a valid approximation for a range of hydrous silicate melts and the partial molar volume of water in these systems has been constrained. Thus, the results of this study provide a useful thermodynamic description of hydrous melts that can be readily applied to other melt-water systems for which spectroscopic measurements of the H-bearing species are available.

The Mode of Action of Recombinant Mycobacterium tuberculosis Shikimate Kinase: Kinetics and Thermodynamics Analyses

Tuberculosis remains as one of the main cause of mortality worldwide due to a single infectious agent, Mycobacterium tuberculosis. The aroK-encoded M. tuberculosis Shikimate Kinase (MtSK), shown to be essential for survival of bacilli, catalyzes the phosphoryl transfer from ATP to the carbon-3 hydroxyl group of shikimate (SKH), yielding shikimate-3-phosphate and ADP. Here we present purification to homogeneity, and oligomeric state determination of recombinant MtSK. Biochemical and biophysical data suggest that the chemical reaction catalyzed by monomeric MtSK follows a rapid-equilibrium random order of substrate binding, and ordered product release. Isothermal titration calorimetry (ITC) for binding of ligands to MtSK provided thermodynamic signatures of non-covalent interactions to each process. A comparison of steady-state kinetics parameters and equilibrium dissociation constant value determined by ITC showed that ATP binding does not increase the affinity of MtSK for SKH. We suggest that MtSK would more appropriately be described as an aroL-encoded type II shikimate kinase. Our manuscript also gives thermodynamic description of SKH binding to MtSK and data for the number of protons exchanged during this bimolecular interaction. The negative value for the change in constant pressure heat capacity (ΔCp) and molecular homology model building suggest a pronounced contribution of desolvation of non-polar groups upon binary complex formation. Thermodynamic parameters were deconvoluted into hydrophobic and vibrational contributions upon MtSK:SKH binary complex formation. Data for the number of protons exchanged during this bimolecular interaction are interpreted in light of a structural model to try to propose the likely amino acid side chains that are the proton donors to bulk solvent following MtSK:SKH complex formation. © 2013 Rosado et al.

Effects of photobioreactor configuration, nitrogen source and light intensity on the fed-batch cultivation of Arthrospira (Spirulina) platensis. Bioenergetic aspects

Bioenergetic analysis may be applied in order to predict microbial growth yields, based on the Gibbs energy dissipation and mass conservation principles of the overall growth reaction. The bioenergetics of the photoautotrophic growth of the cyanobacterium Arthrospira (Spirulina) platensis was investigated in different bioreactor configurations (tubular photobioreactor and open ponds) using different nitrogen sources (nitrate and urea) and under different light intensity conditions to determine the best growing conditions in terms of Gibbs energy dissipation, number of photons to sustain cell growth and phototrophic energy yields distribution in relation to the ATP and NADPH formation, and release of heat. Although an increase in the light intensity increased the Gibbs energy dissipated for cell growth and maintenance with both nitrogen sources, it did not exert any appreciable influence on the moles of photons absorbed by the system to produce one C-mol biomass. On the other hand, both bioenergetic parameters were higher in cultures with nitrate than with urea, likely because of the higher energy requirements needed to reduce the former nitrogen source to ammonia. They appreciably increased also when open ponds were substituted by the tubular photobioreactor, where a more efficient light distribution ensured a remarkably higher cell mass concentration. The estimated percentages of the energy absorbed by the cell showed that, compared with nitrate, the use of urea as nitrogen source allowed the system to address higher energy fractions to ATP production and light fixation by the photosynthetic apparatus, as well as a lower fraction released as heat. The best energy yields values on Gibbs energy necessary for cell growth and maintenance were achieved in up to 4-5 days of cultivation, indicating that it would be the optimum range to maintain cell growth. Thanks to this better bioenergetic situation, urea appears to be a quite promising low-cost, alternative nitrogen source for Arthrospira platensis cultures in photobioreactors. (C) 2011 Elsevier Ltd. All rights reserved.

Proceedings of the Fifteenth Annual biochemical Engineering Symposium

This work represents the proceedings of the fifteenth symposium which convened at Colorado State University on May 24, 1985. The two day meeting was scheduled one month later than usual, i.e., after the spring semester, so that travelers from the Midwest (Iowa State University, Kansas State University and University of Missouri) could enjoy the unique mountain setting provided at Pingree Park. The background of the photograph on the cover depicts the beauty of the area. ContentsGreg Sinton and S.M. Leo, KSU. Models for the Biodegration of 2.4-D and Related Xenobiotic Compounds. V. Bringi, CSU. Intrinsic Kinetics from a Novel Immobilized Cell CSTR. Steve Birdsell, CU. Novel Microbial Separation Techniques. Mark Smith, MU. Kinetic Characterization of Growth of E. coli on Glucose. Michael M. Meagher, ISU. Kinetic Parameters of Di- and Trisaccharaide Hydrolysis by Glucoamylase II. G.T. Jones and A.K. Ghosh Hajra, KSU. Modeling and Simulation of Legume Modules with Reactive Cores and Inert Shells. S.A. Patel and C.H. Lee, KSU. Energetic Analysis and Liquid Circulation in an Airlift Fermenter. Rod R. Fisher, ISU. The Effects of Mixing during Acid Addition of Fractionally Precipitated Protein. Mark M. Paige, CSU. Fed-batch Fermentations of Clostridium acetobutylicum. Michael K. Dowd, ISU. A Nonequilibirium Thermodynamic Description of the Variation of Contractile Velocity and Energy Use in Muscle. David D. Drury, CSU. Analysis of Hollow Fiber Bioreactor Performance for MAmmalian Cells by On-Line MMR. H.Y. Lee, KSU. Process Analysis of Photosynthetic Continuous Culture Systems. C.J. Wang, MU. Kinetic Consideration in Fermentation of Cheese Whey to Ethanol.

Description of the behavior of dichloroalkanes-containing solutions with three [bXmpy][4] isomer, using the experimental information of thermodynamic properties, 1H-NMR spectral and the COSMO-RS methodology

[EN]This work studies the binaries of 1-butyl-X-methylpyridinium tetrafluoroborate [bXmpy][BF4] (X = 2, 3, and 4) with four 1,ω-dichloroalkanes, ω = 1−4, using the results obtained for the mixing properties hE and v E at two temperatures. The three isomers of the ionic liquid (IL) are weakly miscible with the 1,ω-dichloroalkanes when ω ≥ 5 and moderately soluble for ω = 4. The vE s of all the binaries present contractive effects, v E < 0, which are more pronounced with increasing temperature; the variation in vE with ω is positive, although this changes after ω = 4 due to problems of immiscibility

Betwixt kinetic and thermodynamic control:'Pseudo-Thermodynamic control'

A third mechanistic category of organic reactions, conceptually in-between kinetic and thermodynamic control, needs to be recognized, in order to obtain a proper mechanistic description of certain highly important organic reactions, but also to avoid controversy. In fact, an instance of debatable conclusions having been reached is already known. The term �pseudo-thermodynamic control� is proposed for this newly-defined category, which would include any reaction in which intermediates along different routes mutually interconvert. The new category would thus encompass a very large number of reactions, many of them of great utility.

Representation of thermodynamic properties of ternary systems and its application to the system silver-gold-copper at 1350 K

The article presents a generalized analytical expression for description of the integral excess Gibbs free energy of mixing of a ternary system. Twelve constants of the equation are assessed by the least mean squares regressional analysis of the experimental integral excess data of the constituent binaries; three ternary parameters are evaluated by a regressional analysis based on the partial experimental data of a component of the ternary system. The assessed values of the ternary parameters describe the nature of the ternary interaction in the system. Activities and isoactivities of the components in the Ag-Au-Cu system at 1350 K are calculated and found to be in good agreement with the experimental data. This analytical treatment is particularly useful to ternary systems where the thermodynamic data are available from different sources.

Activities in the spinel solid solution, phase equilibria and thermodynamic properties of ternary phases in the system Cu Fe O

A review of the structural and thermodynamic information and phase equilibria in the Cu-Fe-O system suggested that a consistent, quantitative description of the system is hampered by lack of data on activities in the spinel solid solution CuFe2O4-Fe3O4. Therefore the activity of Fe3O4 in this solid solution is derived from measurements of the oxygen potentials established at 1000°C by mixtures containing Fe2O3 and spinel solid solutions of known composition. The oxygen pressures were measured manometrically for solid solutions rich in CuFe2O4, while for Fe3O4-rich compositions the oxygen potentials were obtained by an emf technique. The activities show significant negative deviations from Raoult’s law. The compositions of the spinel solid solutions in equilibrium with CuO + CuFeO2 and Cu + CuFeO2 were obtained from chemical analysis of the solid solution after magnetic separation. The oxygen potential of the three-phase mixture Cu + CuFeO2 + Fe3O4(spinel s.s.) was determined by a solid oxide galvanic cell. From these measurements a complete phase diagram and consistent thermodynamic data on the ternary condensed phases, CuFeO2 and CuFeO2O4, were obtained. An analysis of the free energy of mixing of the spinel solid solution furnished information on the distribution of cations and their valencies between the tetrahedral and octahedral sites of the spinel lattice, which is consistent with X-ray diffraction, magnetic and Seebeck coefficient measurements.

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Thermodynamical description of the interaction between holographic dark energy and dark matter

We present a thermodynamical description of the interaction between holographic dark energy and dark matter. If holographic dark energy and dark matter evolve separately, each of them remains in thermodynamic equilibrium. A small interaction between them may be viewed as a stable thermal fluctuation that brings a logarithmic correction to the equilibrium entropy. From this correction we obtain a physical expression for the interaction which is consistent with phenomenological descriptions and passes reasonably well the observational tests: (c) 2008 Elsevier B.V. All rights reserved.

Scale bridging concepts in molecular simulation : coarse-graining and thermodynamic coupling

In dieser Arbeit wurden Simulation von Flüssigkeiten auf molekularer Ebene durchgeführt, wobei unterschiedliche Multi-Skalen Techniken verwendet wurden. Diese erlauben eine effektive Beschreibung der Flüssigkeit, die weniger Rechenzeit im Computer benötigt und somit Phänomene auf längeren Zeit- und Längenskalen beschreiben kann.rnrnEin wesentlicher Aspekt ist dabei ein vereinfachtes (“coarse-grained”) Modell, welches in einem systematischen Verfahren aus Simulationen des detaillierten Modells gewonnen wird. Dabei werden ausgewählte Eigenschaften des detaillierten Modells (z.B. Paar-Korrelationsfunktion, Druck, etc) reproduziert.rnrnEs wurden Algorithmen untersucht, die eine gleichzeitige Kopplung von detaillierten und vereinfachten Modell erlauben (“Adaptive Resolution Scheme”, AdResS). Dabei wird das detaillierte Modell in einem vordefinierten Teilvolumen der Flüssigkeit (z.B. nahe einer Oberfläche) verwendet, während der Rest mithilfe des vereinfachten Modells beschrieben wird.rnrnHierzu wurde eine Methode (“Thermodynamische Kraft”) entwickelt um die Kopplung auch dann zu ermöglichen, wenn die Modelle in verschiedenen thermodynamischen Zuständen befinden. Zudem wurde ein neuartiger Algorithmus der Kopplung beschrieben (H-AdResS) der die Kopplung mittels einer Hamilton-Funktion beschreibt. In diesem Algorithmus ist eine zur Thermodynamischen Kraft analoge Korrektur mit weniger Rechenaufwand möglich.rnrnAls Anwendung dieser grundlegenden Techniken wurden Pfadintegral Molekulardynamik (MD) Simulationen von Wasser untersucht. Mithilfe dieser Methode ist es möglich, quantenmechanische Effekte der Kerne (Delokalisation, Nullpunktsenergie) in die Simulation einzubeziehen. Hierbei wurde zuerst eine Multi-Skalen Technik (“Force-matching”) verwendet um eine effektive Wechselwirkung aus einer detaillierten Simulation auf Basis der Dichtefunktionaltheorie zu extrahieren. Die Pfadintegral MD Simulation verbessert die Beschreibung der intra-molekularen Struktur im Vergleich mit experimentellen Daten. Das Modell eignet sich auch zur gleichzeitigen Kopplung in einer Simulation, wobei ein Wassermolekül (beschrieben durch 48 Punktteilchen im Pfadintegral-MD Modell) mit einem vereinfachten Modell (ein Punktteilchen) gekoppelt wird. Auf diese Weise konnte eine Wasser-Vakuum Grenzfläche simuliert werden, wobei nur die Oberfläche im Pfadintegral Modell und der Rest im vereinfachten Modell beschrieben wird.