Multivariate Statistical Analysis Applied to an IL6 Signal Transduction Model in Hepatocytes

This paper introduces the application of linear multivariate statistical techniques, including partial least squares (PLS), canonical correlation analysis (CCA) and reduced rank regression (RRR), into the area of Systems Biology. This new approach aims to extract the important proteins embedded in complex signal transduction pathway models.The analysis is performed on a model of intracellular signalling along the janus-associated kinases/signal transducers and transcription factors (JAK/STAT) and mitogen activated protein kinases (MAPK) signal transduction pathways in interleukin-6 (IL6) stimulated hepatocytes, which produce signal transducer and activator of transcription factor 3 (STAT3).A region of redundancy within the MAPK pathway that does not affect the STAT3 transcription was identified using CCA. This is the core finding of this analysis and cannot be obtained by inspecting the model by eye. In addition, RRR was found to isolate terms that do not significantly contribute to changes in protein concentrations, while the application of PLS does not provide such a detailed picture by virtue of its construction.This analysis has a similar objective to conventional model reduction techniques with the advantage of maintaining the meaning of the states prior to and after the reduction process. A significant model reduction is performed, with a marginal loss in accuracy, offering a more concise model while maintaining the main influencing factors on the STAT3 transcription.The findings offer a deeper understanding of the reaction terms involved, confirm the relevance of several proteins to the production of Acute Phase Proteins and complement existing findings regarding cross-talk between the two signalling pathways.

ORBITAL MIXING IN CO CHEMISORPTION ON TRANSITION-METAL SURFACES

The chemisorption of CO on metal surfaces is widely accepted as a model for understanding chemical bonding between molecules and solid surfaces, but is nevertheless still a controversial subject. Ab initio total energy calculations using density functional theory with gradient corrections for CO chemisorption on an extended Pd{110} slab yield good agreement with experimental adsorption energies. Examination of the spatial distribution of individual Bloch states demonstrates that the conventional model for CO chemisorption involving charge donation from CO 5 sigma states to metal states and back-donation from metal states into CO 2 pi states is too simplistic, but the computational results provide direct insight into the chemical bonding within the framework of orbital mixing (or hybridisation). The results provide a sound basis for understanding the bonding between molecules and metal surfaces.

Assessment of a zoomed global model for the North Sea by comparison with a conventional nested regional model

Assessment of marine downscaling of global model simulations to the regional scale is a prerequisite for understanding ocean feedback to the atmosphere in regional climate downscaling. Major difficulties arise from the coarse grid resolution of global models, which cannot provide sufficiently accurate boundary values for the regional model. In this study, we first setup a stretched global model (MPIOM) to focus on the North Sea by shifting poles. Second, a regional model (HAMSOM) was performed with higher resolution, while the open boundary values were provided by the stretched global model. In general, the sea surface temperatures (SSTs) in the two experiments are similar. Major SST differences are found in coastal regions (root mean square difference of SST is reaching up to 2°C). The higher sea surface salinity in coastal regions in the global model indicates the general limitation of this global model and its configuration (surface layer thickness is 16 m). By comparison, the advantage of the absence of open lateral boundaries in the global model can be demonstrated, in particular for the transition region between the North Sea and Baltic Sea. On long timescales, the North Atlantic Current (NAC) inflow through the northern boundary correlates well between both model simulations (R~0.9). After downscaling with HAMSOM, the NAC inflow through the northern boundary decreases by ~10%, but the circulation in the Skagerrak is stronger in HAMSOM. The circulation patterns of both models are similar in the northern North Sea. The comparison suggests that the stretched global model system is a suitable tool for long-term free climate model simulations, and the only limitations occur in coastal regions. Regarding the regional studies focusing on the coastal zone, nested regional model can be a helpful alternative.

Controlled release of a model antibacterial drug from a novel self-lubricating silicone biomaterial

Abstract There is considerable interest in developing medical devices that provide controlled delivery of biologically active agents, for example, to reduce the incidence of device-related infection. Silicone elastomers are one of the commonest biomaterials used in medical device production. However, they have a relatively high coefficient of friction and the resulting lack of lubricity can cause pain and tissue damage on device insertion and removal. Novel silicone cross-linking agents have recently been reported that produce inherently ‘self-lubricating’ silicone elastomers with very low coefficients of friction. In this study, the model antibacterial drug metronidazole has been incorporated into these self-lubricating silicone elastomers to produce a novel bioactive biomaterial. The in vitro release characteristics of the bioactive component were evaluated as a function of cross-linker composition and drug loading. Although conventional matrix-type release kinetics were observed for metronidazole from the silicone systems, it was also observed that increasing the concentration of the cross-linking agent responsible for the lubricious character (tetra(oleyloxy)silane) relative to that of the standard non-lubricious cross-linking agent (tetrapropoxysilane) produced an increase in the metronidazole flux rate by up to 65% for a specified drug loading. The results highlight the potential for developing lubricious silicone medical devices with enhanced drug release characteristics.

Towards An Acoustic Model-Based Poroelastic Imaging Method: I. Theoretical Foundation

The ultrasonic measurement and imaging of tissue elasticity is currently under wide investigation and development as a clinical tool for the assessment of a broad range of diseases, but little account in this field has yet been taken of the fact that soft tissue is porous and contains mobile fluid. The ability to squeeze fluid out of tissue may have implications for conventional elasticity imaging, and may present opportunities for new investigative tools. When a homogeneous, isotropic, fluid-saturated poroelastic material with a linearly elastic solid phase and incompressible solid and fluid constituents is subjected to stress, the behaviour of the induced internal strain field is influenced by three material constants: the Young's modulus (E(s)) and Poisson's ratio (nu(s)) of the solid matrix and the permeability (k) of the solid matrix to the pore fluid. New analytical expressions were derived and used to model the time-dependent behaviour of the strain field inside simulated homogeneous cylindrical samples of such a poroelastic material undergoing sustained unconfined compression. A model-based reconstruction technique was developed to produce images of parameters related to the poroelastic material constants (E(s), nu(s), k) from a comparison of the measured and predicted time-dependent spatially varying radial strain. Tests of the method using simulated noisy strain data showed that it is capable of producing three unique parametric images: an image of the Poisson's ratio of the solid matrix, an image of the axial strain (which was not time-dependent subsequent to the application of the compression) and an image representing the product of the aggregate modulus E(s)(1-nu(s))/(1+nu(s))(1-2nu(s)) of the solid matrix and the permeability of the solid matrix to the pore fluid. The analytical expressions were further used to numerically validate a finite element model and to clarify previous work on poroelastography.

Speech Recognition with unknown partial feature corruption - a review of the union model

This paper provides a summary of our studies on robust speech recognition based on a new statistical approach – the probabilistic union model. We consider speech recognition given that part of the acoustic features may be corrupted by noise. The union model is a method for basing the recognition on the clean part of the features, thereby reducing the effect of the noise on recognition. To this end, the union model is similar to the missing feature method. However, the two methods achieve this end through different routes. The missing feature method usually requires the identity of the noisy data for noise removal, while the union model combines the local features based on the union of random events, to reduce the dependence of the model on information about the noise. We previously investigated the applications of the union model to speech recognition involving unknown partial corruption in frequency band, in time duration, and in feature streams. Additionally, a combination of the union model with conventional noise-reduction techniques was studied, as a means of dealing with a mixture of known or trainable noise and unknown unexpected noise. In this paper, a unified review, in the context of dealing with unknown partial feature corruption, is provided into each of these applications, giving the appropriate theory and implementation algorithms, along with an experimental evaluation.

On velocity-based local model networks for nonlinear identification

This paper exposes the strengths and weaknesses of the recently proposed velocity-based local model (LM) network. The global dynamics of the velocity-based blended representation are directly related to the dynamics of the underlying local models, an important property in the design of local controller networks. Furthermore, the sub-models are continuous-time and linear providing continuity with established linear theory and methods. This is not true for the conventional LM framework, where the global dynamics are only weakly related to the affine sub-models. In this paper, a velocity-based multiple model network is identified for a highly nonlinear dynamical system. The results show excellent dynamical modelling performances, highlighting the value of the velocity-based approach for the design and analysis of LM based control. Three important practical issues are also addressed. These relate to the blending of the velocity-based local models, the use of normalised Gaussian basis functions and the requirement of an input derivative.

A hierarchical physical model-based approach to predictive control of a thermal power plant for efficient plant-wide disturbance rejection

Extending the work presented in Prasad et al. (IEEE Proceedings on Control Theory and Applications, 147, 523-37, 2000), this paper reports a hierarchical nonlinear physical model-based control strategy to account for the problems arising due to complex dynamics of drum level and governor valve, and demonstrates its effectiveness in plant-wide disturbance handling. The strategy incorporates a two-level control structure consisting of lower-level conventional PI regulators and a higher-level nonlinear physical model predictive controller (NPMPC) for mainly set-point manoeuvring. The lower-level PI loops help stabilise the unstable drum-boiler dynamics and allow faster governor valve action for power and grid-frequency regulation. The higher-level NPMPC provides an optimal load demand (or set-point) transition by effective handling of plant-wide interactions and system disturbances. The strategy has been tested in a simulation of a 200-MW oil-fired power plant at Ballylumford in Northern Ireland. A novel approach is devized to test the disturbance rejection capability in severe operating conditions. Low frequency disturbances were created by making random changes in radiation heat flow on the boiler-side, while condenser vacuum was fluctuating in a random fashion on the turbine side. In order to simulate high-frequency disturbances, pulse-type load disturbances were made to strike at instants which are not an integral multiple of the NPMPC sampling period. Impressive results have been obtained during both types of system disturbances and extremely high rates of load changes, right across the operating range, These results compared favourably with those from a conventional state-space generalized predictive control (GPC) method designed under similar conditions.

A heuristic model for the simulation of the deformation of elastic and spongy material for virtual reality applications

This paper presents a practical algorithm for the simulation of interactive deformation in a 3D polygonal mesh model. The algorithm combines the conventional simulation of deformation using a spring-mass-damping model, solved by explicit numerical integration, with a set of heuristics to describe certain features of the transient behaviour, to increase the speed and stability of solution. In particular, this algorithm was designed to be used in the simulation of synthetic environments where it is necessary to model realistically, in real time, the effect on non-rigid surfaces being touched, pushed, pulled or squashed. Such objects can be solid or hollow, and have plastic, elastic or fabric-like properties. The algorithm is presented in an integrated form including collision detection and adaptive refinement so that it may be used in a self-contained way as part of a simulation loop to include human interface devices that capture data and render a realistic stereoscopic image in real time. The algorithm is designed to be used with polygonal mesh models representing complex topology, such as the human anatomy in a virtual-surgery training simulator. The paper evaluates the model behaviour qualitatively and then concludes with some examples of the use of the algorithm.

A new model based approach for the prediction and optimisation of thermal homogeneity in single screw extrusion

In polymer extrusion, delivery of a melt which is homogenous in composition and temperature is important for good product quality. However, the process is inherently prone to temperature fluctuations which are difficult to monitor and control via single point based conventional thermo- couples. In this work, the die melt temperature profile was monitored by a thermocouple mesh and the data obtained was used to generate a model to predict the die melt temperature profile. A novel nonlinear model was then proposed which was demonstrated to be in good agreement with training and unseen data. Furthermore, the proposed model was used to select optimum process settings to achieve the desired average melt temperature across the die while improving the temperature homogeneity. The simulation results indicate a reduction in melt temperature variations of up to 60%.

Order Flow and the Monetary Model of Exchange Rates: Evidence from a Novel Data Set

We propose an exchange rate model that is a hybrid of the conventional specification with monetary fundamentals and the Evans–Lyons microstructure approach. We estimate a model augmented with order flow variables, using a unique data set: almost 100 monthly observations on interdealer order flow on dollar/euro and dollar/yen. The augmented macroeconomic, or “hybrid,” model exhibits greater in-sample stability and out of sample forecasting improvement vis-à-vis the basic macroeconomic and random walk specifications.

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments. (C) 2004 by Radiation Research Society.

Optimal diet choice for large herbivores: an extended contingency model

1. A more general contingency model of optimal diet choice is developed, allowing for simultaneous searching and handling, which extends the theory to include grazing and browsing by large herbivores.</p><p>2. Foraging resolves into three modes: purely encounter-limited, purely handling-limited and mixed-process, in which either a handling-limited prey type is added to an encounter-limited diet, or the diet becomes handling-limited as it expands.</p><p>3. The purely encounter-limited diet is, in general, broader than that predicted by the conventional contingency model,</p><p>4. As the degree of simultaneity of searching and handling increases, the optimal diet expands to the point where it is handling-limited, at which point all inferior prey types are rejected,</p><p>5. Inclusion of a less profitable prey species is not necessarily independent of its encounter rate and the zero-one rule does not necessarily hold: some of the less profitable prey may be included in the optimal diet. This gives an optimal foraging explanation for herbivores' mixed diets.</p><p>6. Rules are shown for calculating the boundary between encounter-limited and handling-limited diets and for predicting the proportion of inferior prey to be included in a two-species diet,</p><p>7. The digestive rate model is modified to include simultaneous searching and handling, showing that the more they overlap, the more the predicted diet-breadth is likely to be reduced.</p>

A material model for multiaxial stretching and stress relaxation of polypropylene under process conditions

Polypropylene sheets have been stretched at 160 °C to a state of large biaxial strain of extension ratio 3, and the stresses then allowed to relax at constant strain. The state of strain is reached via a path consisting of two sequential planar extensions, the second perpendicular to the first, under plane stress conditions with zero stress acting normal to the sheet. This strain path is highly relevant to solid phase deformation processes such as stretch blow moulding and thermoforming, and also reveals fundamental aspects of the flow rule required in the constitutive behaviour of the material. The rate of decay of stress is rapid, and such as to be highly significant in the modelling of processes that include stages of constant strain. A constitutive equation is developed that includes Eyring processes to model both the stress relaxation and strain rate dependence of the stress. The axial and transverse stresses observed during loading show that the use of a conventional Levy-Mises flow rule is ineffective, and instead a flow rule is used that takes account of the anisotropic state of the material via a power law function of the principal extension ratios. Finally the constitutive model is demonstrated to give quantitatively useful representation of the stresses both in loading and in stress relaxation.

Conventional in vivo irradiation procedures are insufficient to accurately determine tumor responses to non-uniform radiation fields

Purpose: To determine differences in overall tumor responses measured by volumetric assessment and bioluminescence imaging (BLI) following exposure to uniform and non-uniform radiation fields in an ectopic prostate tumor model.

Materials and methods: Bioluminescent human prostate tumor xenografts were established by subcutaneous implantation into male mice. Tumors were irradiated with uniform or non-uniform field configurations using conventional in vivo irradiation procedures performed using a 225 kVp generator with custom lead shielding. Tumor responses were measured using Vernier calipers and by BLI using an in vivo imaging system. Survival was defined as the time to quadroupling of pre-treatment tumor volume. 

Results: The correlation between BLI and tumor volume measurements was found to be different for un-irradiated (R = 0.61), uniformly irradiated (R = 0.34) and partially irradiated (R = 0.30) tumors. Uniformly irradiated tumors resulted in an average tumor growth delay of 60 days with median survival of 75 days, compared to partially irradiated tumors which showed an average growth delay of 24 days and median survival of 38 days. 

Conclusions: Correlation between BLI and tumor volume measurements is lower for partially irradiated tumors than those exposed to uniform dose distributions. The response of partially irradiated tumors suggests non-uniformity in response beyond physical dose distribution within the target volume. Dosimetric uncertainty associated with conventional in vivo irradiation procedures prohibits their ability to accurately determine tumor response to non-uniform radiation fields and stresses the need for image guided small animal radiation research platforms.