On systems thinking, systems biology and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

Modeling hybridoma cell metabolism using a generic genome-scale metabolic model of Mus musculus

The reconstructed cellular metabolic network of Mus musculus, based on annotated genomic data, pathway databases, and currently available biochemical and physiological information, is presented. Although incomplete, it represents the first attempt to collect and characterize the metabolic network of a mammalian cell on the basis of genomic data. The reaction network is generic in nature and attempts to capture the carbon, energy, and nitrogen metabolism of the cell. The metabolic reactions were compartmentalized between the cytosol and the mitochondria, including transport reactions between the compartments and the extracellular medium. The reaction list consists of 872 internal metabolites involved in a total of 1220 reactions, whereof 473 relate to known open reading frames. Initial in silico analysis of the reconstructed model is presented.

Studies on the membrane interactions of the cyclotides kalata B1 and kalata B6 on model membrane systems by surface plasmon resonance

In this study we have demonstrated the interactions of kalata B1 and its naturally occurring analogue kalata B6 with five model lipid membranes and have analyzed the binding kinetics using surface plasmon resonance. Two kalata peptides showed a higher affinity for the phosphatidylethanolamine-containing membranes, indicating that the peptides would bind selectively to bacterial membranes. Also we have optimized the procedure for the immobilization of five liposome mixtures and have shown that the procedure provides reproducible levels of immobilized liposomes and could be used to screen the selective binding of putative antimicrobial peptides to model mammalian or microbial phospholipid membranes. (C) 2004 Elsevier Inc. All rights reserved.

Anaerobic metabolism of propionate by polyphosphate-accumulating organisms in enhanced biological phosphorus removal systems

Propionate, a carbon substrate abundant in many prefermenters, has been shown in several previous studies to be a more favorable substrate than acetate for enhanced biological phosphorus removal (EBPR). The anaerobic metabolism of propionate by polyphosphate accumulating organisms (PAOs) is studied in this paper. A metabolic model is proposed to characterize the anaerobic biochemical transformations of propionate uptake by PAOs. The model is demonstrated to predict very well the experimental data from a PAO culture enriched in a laboratory-scale reactor with propionate as the sole carbon source. Quantitative fluorescence in-situ hybridization (FISH) analysis shows that Candidatus Accumulibacter phosphatis, the only identified PAO to date, constitute 63% of the bacterial population in this culture. Unlike the anaerobic metabolism of acetate by PAOs, which induces mainly poly-beta-hydroxybutyrate (PHB) production, the major fractions of poly-beta-hydroxyalkanoate (PHA) produced with propionate as the carbon source are poly-beta-hydroxyvalerate (PHV) and poly-beta-hydroxy-2-methylvalerate (PH2MV). PHA formation correlates very well with a selective (or nonrandom) condensation of acetyl-CoA and propionyl-CoA molecules. The maximum specific propionate uptake rate by PAOs found in this study is 0.18 C-mol/C-mol-biomass h, which is very similar to the maximum specific acetate uptake rate reported in literature. The energy required for transporting 1 carbon-mole of propionate across the PAO cell membrane is also determined to be similar to the transportation of 1 carbon-mole of acetate. Furthermore, the experimental results suggest that PAOs possess a similar preference toward acetate and propionate uptake on a carbon-mole basis. (c) 2005 Wiley Periodicals, Inc.

Optimisation of poly-b-hydroxyalkanoate analysis using gas chromatography for enhanced biological phosphorus removal systems

Poly-beta-hydroxyalkanoate (PHA) is a polymer commonly used in carbon and energy storage for many different bacterial cells. Polyphosphate accumulating organisms (PAOs) and glycogen accumulating organisms (GAOs), store PHA anaerobically through metabolism of carbon substrates such as acetate and propionate. Although poly-beta-hydroxybutyrate (PHB)and poly-beta-hydroxyvalerate (PHV) are commonly quantified using a previously developed gas chromatography (GC) method, poly-beta-hydroxy-2-methyl valerate (PH2MV) is seldom quantified despite the fact that it has been shown to be a key PHA fraction produced when PAOs or GAOs metabolise propionate. This paper presents two GC-based methods modified for extraction and quantification of PHB, PHV and PH2MV from enhanced biological phosphorus removal (EBPR) systems. For the extraction Of PHB and PHV from acetate fed PAO and GAO cultures, a 3% sulfuric acid concentration and a 2-20 h digestion time is recommended, while a 10% sulfuric acid solution digested for 20 h is recommended for PHV and PH2MV analysis from propionate fed EBPR systems. (c) 2005 Elsevier B.V. All rights reserved.

A windows program for the derivation of steady-state equations in enzyme systems

Despite the number of computer-assisted methods described for the derivation of steady-state equations of enzyme systems, most of them are focused on strict steady-state conditions or are not able to solve complex reaction mechanisms. Moreover, many of them are based on computer programs that are either not readily available or have limitations. We present here a computer program called WinStes, which derives equations for both strict steady-state systems and those with the assumption of rapid equilibrium, for branched or unbranched mechanisms, containing both reversible and irreversible conversion steps. It solves reaction mechanisms involving up to 255 enzyme species, connected by up to 255 conversion steps. The program provides all the advantages of the Windows programs, such as a user-friendly graphical interface, and has a short computation time. WinStes is available free of charge on request from the authors. (c) 2006 Elsevier Inc. All rights reserved.

Stochastic delay differential equations for genetic regulatory networks

Time delay is an important aspect in the modelling of genetic regulation due to slow biochemical reactions such as gene transcription and translation, and protein diffusion between the cytosol and nucleus. In this paper we introduce a general mathematical formalism via stochastic delay differential equations for describing time delays in genetic regulatory networks. Based on recent developments with the delay stochastic simulation algorithm, the delay chemical masterequation and the delay reaction rate equation are developed for describing biological reactions with time delay, which leads to stochastic delay differential equations derived from the Langevin approach. Two simple genetic regulatory networks are used to study the impact of' intrinsic noise on the system dynamics where there are delays. (c) 2006 Elsevier B.V. All rights reserved.

Advanced computing for systems biology

Systems biology is based on computational modelling and simulation of large networks of interacting components. Models may be intended to capture processes, mechanisms, components and interactions at different levels of fidelity. Input data are often large and geographically disperse, and may require the computation to be moved to the data, not vice versa. In addition, complex system-level problems require collaboration across institutions and disciplines. Grid computing can offer robust, scaleable solutions for distributed data, compute and expertise. We illustrate some of the range of computational and data requirements in systems biology with three case studies: one requiring large computation but small data (orthologue mapping in comparative genomics), a second involving complex terabyte data (the Visible Cell project) and a third that is both computationally and data-intensive (simulations at multiple temporal and spatial scales). Authentication, authorisation and audit systems are currently not well scalable and may present bottlenecks for distributed collaboration particularly where outcomes may be commercialised. Challenges remain in providing lightweight standards to facilitate the penetration of robust, scalable grid-type computing into diverse user communities to meet the evolving demands of systems biology.

Incorporating level set methods in Geographical Information Systems (GIS) for land-surface process modelling

Land-surface processes include a broad class of models that operate at a landscape scale. Current modelling approaches tend to be specialised towards one type of process, yet it is the interaction of processes that is increasing seen as important to obtain a more integrated approach to land management. This paper presents a technique and a tool that may be applied generically to landscape processes. The technique tracks moving interfaces across landscapes for processes such as water flow, biochemical diffusion, and plant dispersal. Its theoretical development applies a Lagrangian approach to motion over a Eulerian grid space by tracking quantities across a landscape as an evolving front. An algorithm for this technique, called level set method, is implemented in a geographical information system (GIS). It fits with a field data model in GIS and is implemented as operators in map algebra. The paper describes an implementation of the level set methods in a map algebra programming language, called MapScript, and gives example program scripts for applications in ecology and hydrology.

Simultaneous chemical reaction and distillation of formaldehyde

The thesis is concerned with the development and testing of a mathematical model of a distillation process in which the components react chemically. The formaldehyde-methanol-water system was selected and only the reversible reactions between formaldehyde and water giving methylene glycol and between formaldehyde and methanol producing hemiformal were assumed to occur under the distillation conditions. Accordingly the system has been treated as a five component system. The vapour-liquid equilibrium calculations were performed by solving iteratively the thermodynamic relationships expressing the phase equilibria with the stoichiometric equations expressing the chemical equilibria. Using optimisation techniques, the Wilson single parameters and Henry's constants were calculated for binary systems containing formaldehyde which was assumed to be a supercritical component whilst Wilson binary parameters were calculated for the remaining binary systems. Thus the phase equilibria for the formaldehyde system could be calculated using these parameters and good accuracy was obtained when calculated values were compared with experimental values. The distillation process was modelled using the mass and energy balance equations together with the phase equilibria calculations. The plate efficiencies were obtained from a modified A.I.Ch.E. Bubble Tray method. The resulting equations were solved by an iterative plate to plate calculation based on the Newton Raphson method. Experiments were carried out in a 76mm I.D., eight sieve plate distillation column and the results were compared with the mathematical model calculations. Overall, good agreement was obtained but some discrepancies were observed in the concentration profiles and these may have been caused by the effect of limited physical property data and a limited understanding of the reactions mechanism. The model equations were solved in the form of modular computer programs. Although they were written to describe the steady state distillation with simultaneous chemical reaction of the formaldehyde system, the approach used may be of wider application.

Liposomes:Formulation and characterisation as contrast agents and as vaccine delivery systems

Liposome systems are well reported for their activity as vaccine adjuvants; however novel lipid-based microbubbles have also been reported to enhance the targeting of antigens into dendritic cells (DCs) in cancer immunotherapy (Suzuki et al 2009). This research initially focused on the formulation of gas-filled lipid coated microbubbles and their potential activation of macrophages using in vitro models. Further studies in the thesis concentrated on aqueous-filled liposomes as vaccine delivery systems. Initial work involved formulating and characterising four different methods of producing lipid-coated microbubbles (sometimes referred to as gas-filled liposomes), by homogenisation, sonication, a gas-releasing chemical reaction and agitation/pressurisation in terms of stability and physico-chemical characteristics. Two of the preparations were tested as pressure probes in MRI studies. The first preparation composed of a standard phospholipid (DSPC) filled with air or nitrogen (N2), whilst in the second method the microbubbles were composed of a fluorinated phospholipid (F-GPC) filled with a fluorocarbon saturated gas. The studies showed that whilst maintaining high sensitivity, a novel contrast agent which allows stable MRI measurements of fluid pressure over time, could be produced using lipid-coated microbubbles. The F-GPC microbubbles were found to withstand pressures up to 2.6 bar with minimal damage as opposed to the DSPC microbubbles, which were damaged at above 1.3 bar. However, it was also found that DSPC-filled with N2 microbubbles were also extremely robust to pressure and their performance was similar to that of F-GPC based microbubbles. Following on from the MRI studies, the DSPC-air and N2 filled lipid-based microbubbles were assessed for their potential activation of macrophages using in vitro models and compared to equivalent aqueous-filled liposomes. The microbubble formulations did not stimulate macrophage uptake, so studies thereafter focused on aqueous-filled liposomes. Further studies concentrated on formulating and characterising, both physico-chemically and immunologically, cationic liposomes based on the potent adjuvant dimethyldioctadecylammonium (DDA) and immunomodulatory trehalose dibehenate (TDB) with the addition of polyethylene glycol (PEG). One of the proposed hypotheses for the mechanism behind the immunostimulatory effect obtained with DDA:TDB is the ‘depot effect’ in which the liposomal carrier helps to retain the antigen at the injection site thereby increasing the time of vaccine exposure to the immune cells. The depot effect has been suggested to be primarily due to their cationic nature. Results reported within this thesis demonstrate that higher levels of PEG i.e. 25 % were able to significantly inhibit the formation of a liposome depot at the injection site and also severely limit the retention of antigen at the site. This therefore resulted in a faster drainage of the liposomes from the site of injection. The versatility of cationic liposomes based on DDA:TDB in combination with different immunostimulatory ligands including, polyinosinic-polycytidylic acid (poly (I:C), TLR 3 ligand), and CpG (TLR 9 ligand) either entrapped within the vesicles or adsorbed onto the liposome surface was investigated for immunogenic capacity as vaccine adjuvants. Small unilamellar (SUV) DDA:TDB vesicles (20-100 nm native size) with protein antigen adsorbed to the vesicle surface were the most potent in inducing both T cell (7-fold increase) and antibody (up to 2 log increase) antigen specific responses. The addition of TLR agonists poly(I:C) and CpG to SUV liposomes had small or no effect on their adjuvanticity. Finally, threitol ceramide (ThrCer), a new mmunostimulatory agent, was incorporated into the bilayers of liposomes composed of DDA or DSPC to investigate the uptake of ThrCer, by dendritic cells (DCs), and presentation on CD1d molecules to invariant natural killer T cells. These systems were prepared both as multilamellar vesicles (MLV) and Small unilamellar (SUV). It was demonstrated that the IFN-g secretion was higher for DDA SUV liposome formulation (p<0.05), suggesting that ThrCer encapsulation in this liposome formulation resulted in a higher uptake by DCs.

Evaluation of the importance of astrocytes when screening for acute toxicity in neuronal cell systems

Reliable, high throughput, in vitro preliminary screening batteries have the potential to greatly accelerate the rate at which regulatory neurotoxicity data is generated. This study evaluated the importance of astrocytes when predicting acute toxic potential using a neuronal screening battery of pure neuronal (NT2.N) and astrocytic (NT2.A) and integrated neuronal/astrocytic (NT2.N/A) cell systems derived from the human NT2.D1 cell line, using biochemical endpoints (mitochondrial membrane potential (MMP) depolarisation and ATP and GSH depletion). Following exposure for 72 h, the known acute human neurotoxicants trimethyltin-chloride, chloroquine and 6-hydroxydopamine were frequently capable of disrupting biochemical processes in all of the cell systems at non-cytotoxic concentrations. Astrocytes provide key metabolic and protective support to neurons during toxic challenge in vivo and generally the astrocyte containing cell systems showed increased tolerance to toxicant insult compared with the NT2.N mono-culture in vitro. Whilst there was no consistent relationship between MMP, ATP and GSH log IC(50) values for the NT2.N/A and NT2.A cell systems, these data did provide preliminary evidence of modulation of the acute neuronal toxic response by astrocytes. In conclusion, the suitability of NT2 neurons and astrocytes as cell systems for acute toxicity screening deserves further investigation.

The mechanochemical modification of polyolefins using reactive modifier systems

The principal objective of this work was to improve the mechanical properties of glass fibre reinforced polypropylene (PP) composites by the mechanochemical modification of the PP. The modification of the PP was carried out by reactive processing of the PP with a modifier in a Buss Ko-Kneader. Two main types of modifier were evaluated one type based on N-substituted maleimides the others based on 2-allylamino-4,6-dichloro-1,3,5-triazine (ACCT). The modification of the PP was carried out in two stages. Firstly the PP was reactively processed with the modifier and a free radical initiator. The objective of this stage was to bind the modifier to the PP. In the second stage the modified PP was reactively processed with the glass fibre. The objective in this stage was to form a chemical bond between the bound modifier and the silane coupling agent on the surface of the glass. Two silane coupling agents were evaluated these had a aliphatic chloro group and an aliphatic amino group respectively available for reaction with the modifier. The modifiers synthesised for this work had two main functional groups. The first was a double bond for free radical addition to the PP. The second was an organic group chosen for its potential reactivity to the silane coupling agent. A preliminary investigation was carried out using maleic anhydride (MA) as the modifier, this is reactive to the amino silane coupled glass. Studies of a commercially available system were also carried out for comparison purposes. During the work it was found that the amino silane coupled glass fibres produced, without any modification being made to the PP, mechanical properties comparable to the commercial system. Further any modification added to the amino silane system failed to improve the mechanical performance and in some cases acted in the opposite fashion. This failure was evident even when a chemical bond between glass fibre and PP could be shown. In the case of the chloro silane coupled glass fibres the mechanical properties of the composite without modification were poorer than those of the commercial system. It was found that the mechanical properties of these systems could be enhanced by the modifiers, however, no system tested significantly out performed the commercial system. Of the two modifier systems tested those based on the n-substituted maleimides were more successful at enhancing mechanical properties than those based on ACCT. This was attributed to the Poor chemical binding of the ACCT based modifiers to the PP. During the work it was found that several of the modifiers improved the properties of the PP when no glass fibres were present, particularly the % elongation and impact strength. It is possible that these modifiers could be used to improve the impact performance of PP, this may be of particular interest in recycling. These modifiers have only been tested for improving the properties of glass fibre composites. The N-substituted maleimide based modifiers could be used as compatibleisers for alloys of PP and other polymers. These could function by the formation of the bond with PP via the double bond whilst the group attached to the nitrogen atom could react with the alloying polymer.

Studies in three phase gas-liquid-solid fluidised systems

The work is a logical continuation of research started at Aston some years ago when studies were conducted on fermentations in bubble columns. The present work highlights typical design and operating problems that could arise in such systems as waste water, chemical, biochemical and petroleum operations involving three-phase, gas-liquid-solid fluidisation; such systems are in increasing use. It is believed that this is one of few studies concerned with `true' three-phase, gas-liquid-solid fluidised systems, and that this work will contribute significantly to closing some of the gaps in knowledge in this area. The research work was mainly experimentally based and involved studies of the hydrodynamic parameters, phase holdups (gas and solid), particle mixing and segregation, and phase flow dynamics (flow regime and circulation patterns). The studies have focused particularly on the solid behaviour and the influence of properties of solids present on the above parameters in three-phase, gas-liquid-solid fluidised systems containing single particle components and those containing binary and ternary mixtures of particles. All particles were near spherical in shape and two particle sizes and total concentration levels were used. Experiments were carried out in two- and three-dimensional bubble columns. Quantitative results are presented in graphical form and are supported by qualitative results from visual studies which are also shown as schematic diagrams and in photographic form. Gas and solid holdup results are compared for air-water containing single, binary and ternary component particle mixtures. It should be noted that the criteria for selection of the materials used are very important if true three-phase fluidisation is to be achieved: this is very evident when comparing the results with those in the literature. The fluid flow and circulation patterns observed were assessed for validation of the generally accepted patterns, and the author believes that the present work provides more accurate insight into the modelling of liquid circulation in bubble columns. The characteristic bubbly flow at low gas velocity in a two-phase system is suppressed in the three-phase system. The degree of mixing within the system is found to be dependent on flow regime, liquid circulation and the ratio of solid phase physical properties. Evidence of strong `trade-off' of properties is shown; the overall solid holdup is believed to be a major parameter influencing the gas holdup structure.