Semi-algebraic methods for symbolic analysis of complex reaction networks

The identification of chemical mechanism that can exhibit oscillatory phenomena in reaction networks are currently of intense interest. In particular, the parametric question of the existence of Hopf bifurcations has gained increasing popularity due to its relation to the oscillatory behavior around the fixed points. However, the detection of oscillations in high-dimensional systems and systems with constraints by the available symbolic methods has proven to be difficult. The development of new efficient methods are therefore required to tackle the complexity caused by the high-dimensionality and non-linearity of these systems. In this thesis, we mainly present efficient algorithmic methods to detect Hopf bifurcation fixed points in (bio)-chemical reaction networks with symbolic rate constants, thereby yielding information about their oscillatory behavior of the networks. The methods use the representations of the systems on convex coordinates that arise from stoichiometric network analysis. One of the methods called HoCoQ reduces the problem of determining the existence of Hopf bifurcation fixed points to a first-order formula over the ordered field of the reals that can then be solved using computational-logic packages. The second method called HoCaT uses ideas from tropical geometry to formulate a more efficient method that is incomplete in theory but worked very well for the attempted high-dimensional models involving more than 20 chemical species. The instability of reaction networks may lead to the oscillatory behaviour. Therefore, we investigate some criterions for their stability using convex coordinates and quantifier elimination techniques. We also study Muldowney's extension of the classical Bendixson-Dulac criterion for excluding periodic orbits to higher dimensions for polynomial vector fields and we discuss the use of simple conservation constraints and the use of parametric constraints for describing simple convex polytopes on which periodic orbits can be excluded by Muldowney's criteria. All developed algorithms have been integrated into a common software framework called PoCaB (platform to explore bio- chemical reaction networks by algebraic methods) allowing for automated computation workflows from the problem descriptions. PoCaB also contains a database for the algebraic entities computed from the models of chemical reaction networks.

The periplasmic domain of the barrel assembly machinery protein A (BamA) from Escherichia coli assists folding of outer membrane protein A

The assembly of outer membranes of the cell wall of Gram-negative bacteria and of various organelles of eukaryotic cells requires the evolutionarily conserved β-barrel-assembly machinery (BAM) complex. This thesis describes the biochemical and biophysical properties of the periplasmic domain of the β-barrel assembly machinery protein A (PD-BamA) of the E. coli BAM complex, its effect on insertion and folding of the Outer membrane protein A (OmpA) into lipid bilayers and the identification of regions of PD-BamA that may be involved in protein-protein interactions. The secondary structure of PD-BamA in mixed lipid bilayers, analyzed by Circular dichroism (CD) spectroscopy, contained less β-sheet at an increased content of phosphatidylglycerol (PG) in the lipid membrane. This result showed membrane binding, albeit only in the presence of negatively charged lipids. Fluorescence spectroscopy demonstrated that PD-BamA only binds to lipid bilayers containing the negatively charged DOPG, confirming the results of CD spectroscopy. PD-BamA did not bind to zwitterionic but overall neutral lipid bilayers. PD-BamA bound to OmpA at a stoichiometry of 1:1. PD-BamA strongly facilitated insertion and folding of OmpA into lipid membranes. Kinetics of PD-BamA mediated folding of OmpA was well described by two parallel folding processes, a fast folding process and a slow folding process, differing by 2-3 orders of magnitude in their rate constants. The folding yields of OmpA depended on the concentration of lipid membranes and also on the lipid head groups. The presence of PD-BamA resulted in increased folding yields of OmpA in negatively charged DOPG, but PD-BamA did not affect the folding kinetics of OmpA into bilayers of zwitterionic but overall neutral lipids. The efficiency of folding and insertion of OmpA into lipid bilayers strongly depended on the ratio PD-BamA/OmpA and was optimal at equimolar concentrations of PD-BamA and OmpA. To examine complexes of unfolded OmpA with PD-BamA in more detail, site-directed spectroscopy was used to explore contact regions in both, PD-BamA and OmpA. Similarly, contact regions were also investigated for another protein complex formed by PD-BamA and the lipoprotein BamD. The obtained data suggest, that the site of interaction on PD-BamA for OmpA might be oriented towards the exterior environment away from the preceding POTRA domains, but that PD-BamA is oriented with its short α-helix α1 of POTRA domain 5 towards the C-terminal end of BamD.