Flux balance analysis of biological systems: applications and challenges

Systems level modelling and simulations of biological processes are proving to be invaluable in obtaining a quantitative and dynamic perspective of various aspects of cellular function. In particular, constraint-based analyses of metabolic networks have gained considerable popularity for simulating cellular metabolism, of which flux balance analysis (FBA), is most widely used. Unlike mechanistic simulations that depend on accurate kinetic data, which are scarcely available, FBA is based on the principle of conservation of mass in a network, which utilizes the stoichiometric matrix and a biologically relevant objective function to identify optimal reaction flux distributions. FBA has been used to analyse genome-scale reconstructions of several organisms; it has also been used to analyse the effect of perturbations, such as gene deletions or drug inhibitions in silico. This article reviews the usefulness of FBA as a tool for gaining biological insights, advances in methodology enabling integration of regulatory information and thermodynamic constraints, and finally addresses the challenges that lie ahead. Various use scenarios and biological insights obtained from FBA, and applications in fields such metabolic engineering and drug target identification, are also discussed. Genome-scale constraint-based models have an immense potential for building and testing hypotheses, as well as to guide experimentation.

Convenient construction of tetrahedral finite elements for the quasi-harmonic equation

It is shown that in the finite-element formulation of the general quasi-harmonic equation using tetrahedral elements, for every member of the element family there exists just one numerical universal matrix indpendent of the size, shape and material properties of the element. Thus the element matrix is conveniently constructed by manipulating this single matrix along with a set of reverse sequence codes at the same time accounting for the size, shape and material properties in a simple manner.

Functional dynamics of the nitrogen balance of sorghum. II. Grain filling period.

Maintenance of green leaf area during grain filling can increase grain yield of sorghum grown under terminal water limitation. This 'stay-green' trait has been related to the nitrogen (N) supply-demand balance during grain filling. This study quantifies the N demand of grain and N translocation rates from leaves and stem and explores effects of genotype and N stress on onset and rate of leaf senescence during the grain filling period. Three hybrids differing in potential height were grown at three levels of N supply under well-watered conditions. Vertical profiles of biomass, leaf area, and N% of leaves, stem and grain were measured at regular intervals. Weekly SPAD chlorophyll readings on main shoot leaves were correlated with observed specific leaf nitrogen (SLN) to derive seasonal patterns of leaf N content. For all hybrids, individual grain N demand was sink determined and was initially met through N translocation from the stem and rachis. Only if this was insufficient did leaf N translocation occur. Maximum N translocation rates from leaves and stem were dependent on their N status. However, the supply of N at canopy scale was also related to the amount of leaf area senescing at any one time. This supply-demand framework for N dynamics explained effects of N stress and genotype on the onset and rate of leaf senescence.

Functional dynamics of the nitrogen balance of sorghum: I. N demand of vegetative plant parts.

Stay-green, an important trait for grain yield of sorghum grown under water limitation, has been associated with a high leaf nitrogen content at the start of grain filling. This study quantifies the N demand of leaves and stems and explores effects of N stress on the N balance of vegetative plant parts of three sorghum hybrids differing in potential crop height. The hybrids were grown under well-watered conditions at three levels of N supply. Vertical profiles of biomass and N% of leaves and stems, together with leaf size and number, and specific leaf nitrogen (SLN), were measured at regular intervals. The hybrids had similar minimum but different critical and maximum SLN, associated with differences in leaf size and N partitioning, the latter associated with differences in plant height. N demand of expanding new leaves was represented by critical SLN, and structural stem N demand by minimum stem N%. The fraction of N partitioned to leaf blades increased under N stress. A framework for N dynamics of leaves and stems is developed that captures effects of N stress and genotype on N partitioning and on critical and maximum SLN.

Identity for harmonic oscillator brackets

An identity satisfied by the harmonic oscillator (Talmi-Moshinsky) brackets is derived from two equivalent methods for evaluating an integral often encountered in cluster model studies.

Managing the grass-legume balance in Stylosanthes scabra cv. Seca pastures in central Queensland

The impacts of 4 grazing strategies (year-long grazing, summer grazing, winter grazing and winter grazing plus spring burning) on the grass:legume balance were studied between 2000 and 2006 in a pasture oversown with Stylosanthes scabra cv. Seca (Seca stylo) in central Queensland. Seasonal rainfall throughout the study was generally below average. Total pasture yields in autumn were higher in the 2 winter grazing than the 2 summer grazing treatments, largely reflecting the sampling time relative to when grazing occurred. There were few differences in Seca composition in autumn, although there was a clear trend for Seca composition to be reduced by winter grazing plus burning. Both the frequency of occurrence and plant density of Seca were higher under the 2 summer grazing treatments and there was also a trend for the density of juvenile plants (<5 cm height) to be higher in the 2 summer grazing treatments. Seca soil seed banks were generally low and were reduced in the winter grazing plus burning treatment in spring 2002. The frequency of the palatable perennial grass Pennisetum ciliaris (Biloela buffel grass) was reduced while that of the 'increaser' species Bothriochloa pertusa (Indian couch grass) and Stachytarpheta jamaicensis (snake weed) increased in the 2 summer grazing treatments compared with the 2 winter grazing treatments. Burning in spring increased soil loss in treatments grazed in winter. Differences in Seca frequency and density but not composition were explained by the 2 summer grazing treatments promoting 'gaps' in the pasture which were then colonised by Seca plants and other 'increaser' species. It was reasoned that, with time, mature Seca plants in the 2 winter grazing treatments would die so that Seca composition would eventually become higher under summer grazing regimes than under winter grazing. It was concluded that limiting grazing to particular seasons can alter legume:grass balance and that a time-frame of 5-8 years with average to good rainfall would be necessary to achieve large shifts in composition.

Nutrient mass balance for Lockyer vegetable production

The potential for fertiliser use in Lockyer Valleys intensive vegetable production to impact on the Moreton Bay Waterways (MBW) is not well defined. Notwithstanding nutrient runoff through soil erosion of agricultural lands has been identified as a process that significantly contributes artificial fertiliser to the MBW (SEQ Healthy Waterways Draft Strategy 2006). In order to better understand this issue the present study undertakes a nutrient mass balance to evaluate nitrogen use efficiency in the intensive horticultural industry of the Lockyer Valley.

Heat balance analysis of single stage Gifford-McMahon cycle cryorefrigerator

A heat balance analysis of single stage Gifford-McMahon cycle cryorefrigerator is presented. Ideal refrigeration, actual refrigeration, net refrigeration and the various losses are tabulated. It is observed that pressure-volume losses account for a major fraction of the total losses.

Analyzing hysteresis behavior of capacitance-voltage characteristics of IZO/C₆₀/pentacene/Au diodes with a hole-transport electron-blocking polyterpenol layer by electric-field-induced optical second-harmonic generation measurement

By using electric-field-induced optical second-harmonic generation (EFISHG) measurement, we analyzed hysteresis behavior of capacitance-voltage (C-V) characteristics of IZO/polyterpenol (PT)/C₆₀/pentacene/Au diodes, where PT layer is actively working as a hole-transport electron-blocking layer. The EFISHG measurement verified the presence of interface accumulated charges in the diodes, and showed that a space charge electric field from accumulated excess electrons (holes) that remain at the PT/C₆₀ (C₆₀/pentacene) interface is responsible for the hysteresis loop observed in the C-V characteristics.

Investigation of interfacial charging and discharging in double-layer pentacene-based metal-insulator-metal device with polyterpenol blocking layer using electric field induced second harmonic generation

Time-resolved electric field induced second harmonic generation technique was used to probe the carrier transients within double-layer pentacene-based MIM devices. Polyterpenol thin films fabricated from non-synthetic environmentally sustainable source were used as a blocking layer to assist in visualisation of single-species carrier transportation during charging and discharging under different bias conditions. Results demonstrated that carrier transients were comprised of charging on electrodes, followed by carrier injection and charging of the interface. Polyterpenol was demonstrated to be a sound blocking material and can therefore be effectively used for probing of double-layer devices using EFISHG.

Soil carbon modelling as a tool for carbon balance studies in forestry

Soils represent a remarkable stock of carbon, and forest soils are estimated to hold half of the global stock of soil carbon. Topical concern about the effects of climate change and forest management on soil carbon as well as practical reporting requirements set by climate conventions have created a need to assess soil carbon stock changes reliably and transparently. The large spatial variability of soil carbon commensurate with relatively slow changes in stocks hinders the assessment of soil carbon stocks and their changes by direct measurements. Models therefore widely serve to estimate carbon stocks and stock changes in soils. This dissertation aimed to develop the soil carbon model YASSO for upland forest soils. The model was aimed to take into account the most important processes controlling the decomposition in soils, yet remain simple enough to ensure its practical applicability in different applications. The model structure and assumptions were presented and the model parameters were defined with empirical measurements. The model was evaluated by studying the sensitivities of the model results to parameter values, by estimating the precision of the results with an uncertainty analysis, and by assessing the accuracy of the model by comparing the predictions against measured data and to the results of an alternative model. The model was applied to study the effects of intensified biomass extraction on the forest carbon balance and to estimate the effects of soil carbon deficit on net greenhouse gas emissions of energy use of forest residues. The model was also applied in an inventory based method to assess the national scale forest carbon balance for Finland’s forests from 1922 to 2004. YASSO managed to describe sufficiently the effects of both the variable litter and climatic conditions on decomposition. When combined with the stand models or other systems providing litter information, the dynamic approach of the model proved to be powerful for estimating changes in soil carbon stocks on different scales. The climate dependency of the model, the effects of nitrogen on decomposition and forest growth as well as the effects of soil texture on soil carbon stock dynamics are areas for development when considering the applicability of the model to different research questions, different land use types and wider geographic regions. Intensified biomass extraction affects soil carbon stocks, and these changes in stocks should be taken into account when considering the net effects of forest residue utilisation as energy. On a national scale, soil carbon stocks play an important role in forest carbon balances.

The harmonic torsional oscillations of a thin disk submerged in a fluid with a surfactant surface layer

A formulation in terms of a Fredholm integral equation of the first kind is given for the axisymmetric problem of a disk oscillating harmonically in a viscous fluid whose surface is contaminated with a surfactant film. The equation of the first kind is converted to a pair of coupled integral equations of the second kind, which are solved numerically. The resistive torque on the disk is evaluated and surface velocity profiles are computed for varying values of the ratio of the coefficient of surface shear viscosity to the coefficient of viscosity of the substrate fluid, and the depth of the disk below the surface.

Rail power conditioner technology in Australian Heavy Haul Railway: A case study

Aurizon, Australia's largest freight railway operator, is investigating the use of Rail Power Conditioner (RPC) technology for load balancing, reactive power compensation and harmonic filtering. The new technology has the capability of replacing Static VAr Compensators (SVC) and Harmonic Filters, and is expected to have a significant impact on the overall costs of railway electrification. This paper presents the theoretical analysis of the real and reactive power flows in an RPC used to balance active powers in an existing V/V feeder station. This informed an RPC feasibility study undertaken at four existing Aurizon's feeder stations with V/V connected transformers.