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.

Monooxygenase Diversity and Oxygen Activation within Polyketide Antibiotic Biosynthesis

Molecular oxygen (O2) is a key component in cellular respiration and aerobic life. Through the redox potential of O2, the amount of free energy available to organisms that utilize it is greatly increased. Yet, due to the nature of the O2 electron configuration, it is non-reactive to most organic molecules in the ground state. For O2 to react with most organic compounds it must be activated. By activating O2, oxygenases can catalyze reactions involving oxygen incorporation into organic compounds. The oxygen activation mechanisms employed by many oxygenases to have been studied, and they often include transition metals and selected organic compounds. Despite the diversity of mechanisms for O2 activation explored in this thesis, all of the monooxygenases studied in the experimental part activate O2 through a transient carbanion intermediate. One of these enzymes is the small cofactorless monooxygenase SnoaB. Cofactorless monooxygenases are unusual oxygenases that require neither transition metals nor cofactors to activate oxygen. Based on our biochemical characterization and the crystal structure of this enzyme, the mechanism most likely employed by SnoaB relies on a carbanion intermediate to activate oxygen, which is consistent with the proposed substrate-assisted mechanism for this family of enzymes. From the studies conducted on the two-component system AlnT and AlnH, both the functions of the NADH-dependent flavin reductase, AlnH, and the reduced flavin dependent monooxygenase, AlnT, were confirmed. The unusual regiochemistry proposed for AlnT was also confirmed on the basis of the structure of a reaction product. The mechanism of AlnT, as with other flavin-dependent monooxygenases, is likely to involve a caged radical pair consisting of a superoxide anion and a neutral flavin radical formed from an initial carbanion intermediate. In the studies concerning the engineering of the S-adenosyl-L-methionine (SAM) dependent 4-O-methylase DnrK and the homologous atypical 10-hydroxylase RdmB, our data suggest that an initial decarboxylation of the substrate is catalyzed by both of these enzymes, which results in the generation of a carbanion intermediate. This intermediate is not essential for the 4-O-methylation reaction, but it is important for the 10-hydroxylation reaction, since it enables substrate-assisted activation of molecular oxygen involving a single electron transfer to O2 from a carbanion intermediate. The only role for SAM in the hydroxylation reaction is likely to be stabilization of the carbanion through the positive charge of the cofactor. Based on the DnrK variant crystal structure and the characterizations of several DnrK variants, the insertion of a single amino acid in DnrK (S297) is sufficient for gaining a hydroxylation function, which is likely caused by carbanion stabilization through active site solvent restriction. Despite large differences in the three-dimensional structures of the oxygenases and the potential for multiple oxygen activation mechanisms, all the enzymes in my studies rely on carbanion intermediates to activate oxygen from either flavins or their substrates. This thesis provides interesting examples of divergent evolution and the prevalence of carbanion intermediates within polyketide biosynthesis. This mechanism appears to be recurrent in aromatic polyketide biosynthesis and may reflect the acidic nature of these compounds, propensity towards hydrogen bonding and their ability to delocalize π-electrons.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.