Convex and toric geometry to analyze complex dynamics in chemical reaction systems

Convex cone, toric variety, graph theory, electrochemical catalysis, oxidation of formic acid, feedback-loopsbifurcations, enzymatic catalysis, Peroxidase reaction, Shil'nikov chaos

Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

Studies on nonlinear dynamics of certain reaction systems

The nonlinear dynamics of certain important reaction systems are discussed and analysed in this thesis. The interest in the theoretical and the experimental studies of chemical reactions showing oscillatory dynamics and associated properties is increasing very rapidly. An attempt is made to study some nonlinear phenomena exhibited by the well known chemical oscillator, the BelousovZhabotinskii reaction whose mathematical properties are much in common with the properties of biological oscillators. While extremely complex, this reaction is still much simpler than biological systems at least from the modelling point of view. A suitable model [19] for the system is analysed and the researcher has studied the limit cycle behaviour of the system, for different values of the stoichiometric parameter f, by keeping the value of the reaction rate (k6) fixed at k6 = l. The more complicated three-variable model is stiff in nature.

A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge–Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E. coli, and conclude with a discussion on the significance of this work.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

Acrylamide and pyrazine formation in model systems containing asparagine

The effect of different sugars and glyoxal on the formation of acrylamide in low-moisture starch-based model systems was studied, and kinetic data were obtained. Glucose was more effective than fructose, tagatose, or maltose in acrylamide formation, whereas the importance of glyoxal as a key sugar fragmentation intermediate was confirmed. Glyoxal formation was greater in model systems containing asparagine and glucose rather than fructose. A solid phase microextraction GC-MS method was employed to determine quantitatively the formation of pyrazines in model reaction systems. Substituted pyrazine formation was more evident in model systems containing fructose; however, the unsubstituted homologue, which was the only pyrazine identified in the headspace of glyoxal-asparagine systems, was formed at higher yields when aldoses were used as the reducing sugar. Highly significant correlations were obtained for the relationship between pyrazine and acrylamide formation. The importance of the tautomerization of the asparagine-carbonyl decarboxylated Schiff base in the relative yields of pyrazines and acrylamide is discussed.

Time-resolved gas-phase kinetic study of the germylene addition reaction, CH3C CCH3

Time-resolved kinetic studies of the reaction of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with 2-butyne, CH3C CCH3. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 300-556 K. The second order rate constants obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.46 +/- 10.06) + (5.16 +/- 10.47) kJ mol(-1)/ RT ln 10 Calculations of the energy surface of the GeC4H8 reaction system were carried out employing the additivity principle, by combining previous quantum chemical calculations of related reaction systems. These support formation of 1,2-dimethylvinylgermylene (rather than 2,3-dimethylgermirene) as the end product. RRKM calculations of the pressure dependence of the reaction are in reasonable agreement with this finding. The reactions of GeH2 with C2H2 and with CH3CRCCH3 are compared and contrasted.

Mechanism and model of the oscillatory electro-oxidation of methanol

A mechanism for the kinetic instabilities observed in the galvanostatic electro-oxidation of methanol is suggested and a model developed. The model is investigated using stoichiometric network analysis as well as concepts from algebraic geometry (polynomial rings and ideal theory) revealing the occurrence of a Hopf and a saddle-node bifurcation. These analytical solutions are confirmed by numerical integration of the system of differential equations. (C) 2010 American Institute of Physics

Metabolic analysis of poly(3-hydroxybutyrate) production by recombinant Escherichia coli

Poly(3-hydroxybutyrate) (PHB) production by fermentation was examined under both restricted- and ample-oxygen supply conditions in a single fed-batch fermentation. Recombinant Escherichia coli transformed with the PHB production plasmid pSYL107 was grown to reach high cell density (227 g/l dry cell weight) with a high PHB content (78% of dry cell weight), using a glucose-based minimal medium. A simple flux model containing 12 fluxes was developed and applied to the fermentation data. A superior closure (95%) of the carbon mass balance was achieved. When the data were put into use, the results demonstrated a surprisingly large excretion of formate and lactate. Even though periods of severe oxygen limitation coincided with rapid acetate and lactate excretion, PHB productivity and carbon utilization efficiency were not significantly impaired. These results are very positive in reducing oxygen demand in an industrial PHA fermentation without sacrificing its PHA productivity, thereby reducing overall production costs.

Investigations into the determinations of hydride-forming elements

Analytical methods for the determination of trace amounts of germanium, tin and arsenic were established using hydride generation coupled with direct current plasma atomic emission spectrometry. A continuous gas flowing batch system for the hydride generation was investigated and was applied to the determination of germanium(Ge), tin(Sn), antimony(Sb) and lead(Pb) (Preliminary results suggest that it is also applicable to arsenic)As) ). With this system, the reproducibility of signals was improved and the determination was speeded up, compared with the conventional batch type hydride generation system. Each determination was complete within one minute. Interferences from a number of transition metal ions, especially from Pd(II), Pt(IV), Ni(II), Cu(II), Co(II), and Fe(II, III), have proven to be very serious under normal conditions, in the determination of germanium, tin, and arsenic. These interference effects were eliminated or significantly reduced in the presence of L-cystine or L-cysteine. Thus, a 10-1000 fold excess of Ni(II), Cu(II), Co(II), Fe(II), Pt(IV), Pd(II), etc. can be tolerated without interference, In the presence of L-cystine or L-cysteine, compared with absence of interference reducing agent. The methods for the determination of Ge, Sn, and As were examined by the analyses of standard reference materials. Interference effects from the sample matrix, for example, in transition metal-rich samples, copper, iron and steel, were eliminated by L-cystine (for As and Sn) and by LI cysteine (for Ge). The analysis of a number of standard reference materials gave excellent results of As and Sn contents in agreement with the certified values, showing there was no systematic interference. The detection limits for both germanium and tin were 20 pg ml- I . Preliminary studies were carried out for the determination of antimony and lead. Antimony was found to react with NaBH4, remaInIng from the previous determinations, giving an analytical signal. A reversed injection manner, i.e., injection of the NaBH4 solution prior to the analyte solution was used to avoid uncertainty caused by residual NaBH4 present and to ensure that an excess of NaB H4 was available. A solution of 0.4% L-cysteine was found to reduce the interference from selected transition metal ions, Co(II), Cu(II), Ni(II) and Pt(IV). Hydrochloric acid - hydrogen peroxide, nitric acid - ammonium persulphate, and potassium dichromate malic acid reaction systems for lead hydride generation were compared. The potassium dichromate - malic acid reaction medium proved to be the best with respect to reproducibility and minimal interference. Cu(II), Ni(II), and Fe(II) caused strong interference In lead determinations, which was not reduced by L-cysteine or Lcystine. Sodium citrate, ascorbic acid, dithizone, thiosemicarbazide and penicillamine reduced interferences to some extent. Further interference reduction studies were carried out uSIng a number of amino acids, glycine, alanine, valine, leucine and histidine, as possible interference reducing agents in the determination of germanium. From glycine, alanine, valine to leucine, the interference reduction effect in germanium determinations decreased. Histidine II was found to be very promising In the reduction of interference. In fact, histidine proved more efficient than L-cystine and L-cysteine In the reduction of interference from Ni(II) in the determination of germanium. Signal enhancement by easily ionized elements (EIEs), usually regarded as an interference effect in analysis by DCP-AES, was studied and successfully applied to advantage in improving the sensitivity and detection limit in the determination of As, Ge, Sn, Sb, and Pb. The effect of alkali and alkaline-earth elements on the determination of the above five hydride forming elements was studied. With the appropriate EIE, a signal enhancement of 40-115% was achieved. Linear calibration and good reproducibility were also obtained in the presence of EIEs. III

Strong hydrogen bonding in organic synthesis

The preparation of phenacyl and para-phenylphenacyl esters, the reactions of carboxylic acids, phenols, 2-nitropropane and alcohols with alkyl halides in the presence of fluoride anion are described. The reactions are thought to be accelerated by the formation of hydrogen bonds between the fluoride anion and the organic electron acceptor. The fluoride ,carboxylic acids, fluoride-phenols and fluoride-2-nitropropane are better reaction systems than the fluoride-alcohol. The source of the fluoride anion and the choice of solvents are also discussed.

Desarrollo de catalizadores en sistemas de reacción avanzados para la producción de hidrógeno a partir de biogás

[ES]En el presente documento se expone el estudio experimental que consiste en la comparación de diferentes catalizadores en sistemas de reacción avanzados, microrreactores, para la producción de hidrógeno a partir de biogás. El hidrógeno es un vector energético que puede emplearse como combustible, y por tanto, un candidato perfecto para sustituir aquellos combustibles provenientes de fuentes fósiles. Este Proyecto de Fin de Grado se ha basado en desarrollar el proceso conocido como tri-reformado de biogás. Esta técnica, al contrario que las técnicas convencionales de producción, presenta numerosas ventajas ya que la materia prima empleada, el biogás, es una fuente de origen renovable. Para llevar a cabo este estudio se han preparado distintos catalizadores, basados todos ellos en platino. La elección de este metal noble es debido a su alto grado de reactividad, especialmente en sistemas de reacción micro-estructurados. Con el objetivo de reducir los costes asociados al propio catalizador, se añadieron, junto con el platino, otros metales más baratos a fin de analizar su idoneidad en el proceso del tri-reformado. Por tanto, partiendo de un catalizador de referencia que contiene un 5% de platino, se prepararon otros catalizadores mantenido constante, en un 5%, la composición de metal total. Estos catalizadores se han denominado de la siguiente manera: 2.5(Pt-Me), siendo Me cada especie metálica diferente. Tras preparar los catalizadores, se llevaron a cabo diferentes experimentos con un reactor convencional de lecho fijo. El objetivo fue establecer unas condiciones de operación tal que asegurasen una adecuada comparación entre catalizadores. De este modo, se realizaron ensayos de actividad catalítica a diferentes temperaturas y velocidades espaciales para someter a los catalizadores a condiciones de operación extremas, bajo las cuales, las diferencias entre catalizadores fueran más notorias. Una vez detalladas las condiciones de operación adecuadas para la comparación de catalizadores, se llevaron a cabo los correspondientes ensayos con los catalizadores impregnados en los sistemas de reacción avanzados, denominados micorreactores. Estos ensayos se basaron en operar a una velocidad espacial constante para estudiar la influencia de la temperatura. Además, se llevaron a cabo ensayos de estabilidad de hasta 110 horas en las condiciones de operación más desfavorables. Por último se compararon las conversiones de CH4 y CO2 y el rendimiento de H2 de todas las formulaciones catalíticas preparadas, y se concluyó que el mejor candidato para el proceso del tri-reformado es el catalizador 2.5(Pt-Pd).

Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.