No abstract

In non-alcoholic fatty liver disease (NAFLD), lipid build-up and the resulting damage is known to occur more severely in pericentral cells. Due to the complexity of studying individual regions of the sinusoid, the causes of this zone specificity and its implications on treatment are largely ignored. In this study, a computational model of liver glucose and lipid metabolism is presented which treats the sinusoid as the repeating unit of the liver rather than the single hepatocyte. This allows for inclusion of zonated enzyme expression by splitting the sinusoid into periportal to pericentral compartments. By simulating insulin resistance (IR) and high intake diets leading to the development of steatosis in the model, we identify key differences between periportal and pericentral cells accounting for higher susceptibility to pericentral steatosis. Secondly, variation between individuals is seen in both susceptibility to steatosis and in its development across the sinusoid. Around 25% of obese individuals do not show excess liver fat, whilst 16% of lean individuals develop NAFLD. Furthermore, whilst pericentral cells tend to show higher lipid levels, variation is seen in the predominant location of steatosis from pericentral to pan-sinusoidal or azonal. Sensitivity analysis was used to identify the processes which have the largest effect on both total hepatic triglyceride levels and on the sinusoidal location of steatosis. As is seen in vivo, steatosis occurs when simulating IR in the model, predominantly due to increased uptake, along with an increase in de novo lipogenesis. Additionally, concentrations of glucose intermediates including glycerol-3-phosphate increased when simulating IR due to inhibited glycogen synthesis. Several differences between zones contributed to a higher susceptibility to steatosis in pericentral cells in the model simulations. Firstly, the periportal zonation of both glycogen synthase and the oxidative phosphorylation enzymes meant that the build-up of glucose intermediates was less severe in the periportal hepatocyte compartments. Secondly, the periportal zonation of the enzymes mediating β-oxidation and oxidative phosphorylation resulted in excess fats being metabolised more rapidly in the periportal hepatocyte compartments. Finally, the pericentral expression of de novo lipogenesis contributed to pericentral steatosis when additionally simulating the increase in sterol-regulatory element binding protein 1c (SREBP-1c) seen in NAFLD patients in vivo. The hepatic triglyceride concentration was predicted to be most sensitive to inter-individual variation in the activity of enzymes which, either directly or indirectly, determine the rate of free fatty acid (FFA) oxidation. The concentration was most strongly dependent on the rate constants for β-oxidation and oxidative phosphorylation. It also showed moderate sensitivity to the rate constants for processes which alter the allosteric inhibition of β-oxidation by acetyl-CoA. The predominant sinusoidal location of steatosis meanwhile was most s...

One of the main challenges in the development of mathematical and computational models of biological systems is the precise estimation of parameter values. Understanding the effects of uncertainties in parameter values on model behaviour is crucial to the successful use of these models. Global sensitivity analysis (SA) can be used to quantify the variability in model predictions resulting from the uncertainty in multiple parameters and to shed light on the biological mechanisms driving system behaviour. We present a new methodology for global SA in systems biology which is computationally efficient and can be used to identify the key parameters and their interactions which drive the dynamic behaviour of a complex biological model. The approach combines functional principal component analysis with established global SA techniques. The methodology is applied to a model of the insulin signalling pathway, defects of which are a major cause of type 2 diabetes and a number of key features of the system are identified.

No abstract

A computational model of the glucagon/insulin-driven liver glucohomeostasis function, focusing on the buffering of glucose into glycogen, has been developed. The model exemplifies an 'engineering' approach to modelling in systems biology, and was produced by linking together seven component models of separate aspects of the physiology. The component models use a variety of modelling paradigms and degrees of simplification. Model parameters were determined by an iterative hybrid of fitting to high-scale physiological data, and determination from small-scale in vitro experiments or molecular biological techniques. The component models were not originally designed for inclusion within such a composite model, but were integrated, with modification, using our published modelling software and computational frameworks. This approach facilitates the development of large and complex composite models, although, inevitably, some compromises must be made when composing the individual models. Composite models of this form have not previously been demonstrated.

No abstract

This paper presents the results of the life cycle assessment (LCA) of 10 integrated waste management systems (IWMSs) for 3 potential post-event site design scenarios of the London Olympic Park. The aim of the LCA study is to evaluate direct and indirect emissions resulting from various treatment options of municipal solid waste (MSW) annually generated on site together with avoided emissions resulting from energy, materials and nutrients recovery. IWMSs are modelled using GaBi v6.0 Product Sustainability software and results are presented based on the CML (v.Nov-10) characterisation method. The results show that IWMSs with advanced thermal treatment (ATT) and incineration with energy recovery have the lowest Global Warming Potential (GWP) than IWMSs where landfill is the primary waste treatment process. This is due to higher direct emissions and lower avoided emissions from the landfill process compared to the emissions from the thermal treatment processes. LCA results demonstrate that significant environmental savings are achieved through substitution of virgin materials with recycled ones. The results of the sensitivity analysis carried out for IWMS 1 shows that increasing recycling rate by 5%, 10% and 15% compared to the baseline scenario can reduce GWP by 8%, 17% and 25% respectively. Sensitivity analysis also shows how changes in waste composition affect the overall result of the system. The outcomes of such assessments provide decision-makers with fundamental information regarding the environmental impacts of different waste treatment options necessary for sustainable waste management planning.

The high-pressure homogenization of Escherichia coli, strain JM101, containing inclusion bodies of recombinant porcine somatotropin was investigated. A novel technique employing an analytical disc centrifuge was used to monitor the disruption. This a direct technique which measures cell disintegration rather than soluble protein release. The technique is particularly suited to measurements where the disruption approaches 100%. The disk centrifuge provides a size distribution of the homogenate, and furnishes evidence for the preferential disruption of larger cells. For E. coli containing inclusion bodies, and increase in the cell feed concentration from 145 g/L (wet weight) to 330 g/L resulted is poorer homogenization. Poorer disruption was also obtained by lowering the feed temperature from 20 degrees C to 5 degrees C. Only slight variations in performance were obtained by increasing the feed pH from 7.5 to 9.0 or by storing the feed at 4 degrees C for 24 h prior to disruption. Comparison with uninduced E. coli strain JM101, showed that the disruption obtained is higher for bacteria containing a recombinant inclusion body.