Nonlinear model structure detection using optimum experimental design and orthogonal least squares

A very efficient learning algorithm for model subset selection is introduced based on a new composite cost function that simultaneously optimizes the model approximation ability and model adequacy. The derived model parameters are estimated via forward orthogonal least squares, but the subset selection cost function includes an A-optimality design criterion to minimize the variance of the parameter estimates that ensures the adequacy and parsimony of the final model. An illustrative example is included to demonstrate the effectiveness of the new approach.

Neurofuzzy network model construction using Bézier-Bernstein polynomial functions

Neurofuzzy modelling systems combine fuzzy logic with quantitative artificial neural networks via a concept of fuzzification by using a fuzzy membership function usually based on B-splines and algebraic operators for inference, etc. The paper introduces a neurofuzzy model construction algorithm using Bezier-Bernstein polynomial functions as basis functions. The new network maintains most of the properties of the B-spline expansion based neurofuzzy system, such as the non-negativity of the basis functions, and unity of support but with the additional advantages of structural parsimony and Delaunay input space partitioning, avoiding the inherent computational problems of lattice networks. This new modelling network is based on the idea that an input vector can be mapped into barycentric co-ordinates with respect to a set of predetermined knots as vertices of a polygon (a set of tiled Delaunay triangles) over the input space. The network is expressed as the Bezier-Bernstein polynomial function of barycentric co-ordinates of the input vector. An inverse de Casteljau procedure using backpropagation is developed to obtain the input vector's barycentric co-ordinates that form the basis functions. Extension of the Bezier-Bernstein neurofuzzy algorithm to n-dimensional inputs is discussed followed by numerical examples to demonstrate the effectiveness of this new data based modelling approach.

c-Jun N-terminal kinase (JNK)-mediated modulation of brain mitochondria function: new target proteins for JNK signalling in mitochondrion-dependent apoptosis

The molecular mechanisms underlying the initiation and control of the release of cytochrome c during mitochondrion-dependent apoptosis are thought to involve the phosphorylation of mitochondrial Bcl-2 and Bcl-x(L). Although the c-Jun N-terminal kinase (JNK) has been proposed to mediate the phosphorylation of Bcl-2/Bcl-x(L) the mechanisms linking the modification of these proteins and the release of cytochrome c remain to be elucidated. This study was aimed at establishing interdependency between JNK signalling and mitochondrial apoptosis. Using an experimental model consisting of isolated, bioenergetically competent rat brain mitochondria, these studies show that (i) JNK catalysed the phosphorylation of Bcl-2 and Bcl-x(L) as well as other mitochondrial proteins, as shown by two-dimensional isoelectric focusing/SDS/PAGE; (ii) JNK-induced cytochrome c release, in a process independent of the permeability transition of the inner mitochondrial membrane (imPT) and insensitive to cyclosporin A; (iii) JNK mediated a partial collapse of the mitochondrial inner-membrane potential (Deltapsim) in an imPT- and cyclosporin A-independent manner; and (iv) JNK was unable to induce imPT/swelling and did not act as a co-inducer, but as an inhibitor of Ca-induced imPT. The results are discussed with regard to the functional link between the Deltapsim and factors influencing the permeability transition of the inner and outer mitochondrial membranes. Taken together, JNK-dependent phosphorylation of mitochondrial proteins including, but not limited to, Bcl-2/Bcl-x(L) may represent a potential of the modulation of mitochondrial function during apoptosis.

Assessment of a climate model to reproduce rainfall variability and extremes over Southern Africa

It is increasingly accepted that any possible climate change will not only have an influence on mean climate but may also significantly alter climatic variability. A change in the distribution and magnitude of extreme rainfall events (associated with changing variability), such as droughts or flooding, may have a far greater impact on human and natural systems than a changing mean. This issue is of particular importance for environmentally vulnerable regions such as southern Africa. The sub-continent is considered especially vulnerable to and ill-equipped (in terms of adaptation) for extreme events, due to a number of factors including extensive poverty, famine, disease and political instability. Rainfall variability and the identification of rainfall extremes is a function of scale, so high spatial and temporal resolution data are preferred to identify extreme events and accurately predict future variability. The majority of previous climate model verification studies have compared model output with observational data at monthly timescales. In this research, the assessment of ability of a state of the art climate model to simulate climate at daily timescales is carried out using satellite-derived rainfall data from the Microwave Infrared Rainfall Algorithm (MIRA). This dataset covers the period from 1993 to 2002 and the whole of southern Africa at a spatial resolution of 0.1° longitude/latitude. This paper concentrates primarily on the ability of the model to simulate the spatial and temporal patterns of present-day rainfall variability over southern Africa and is not intended to discuss possible future changes in climate as these have been documented elsewhere. Simulations of current climate from the UKMeteorological Office Hadley Centre’s climate model, in both regional and global mode, are firstly compared to the MIRA dataset at daily timescales. Secondly, the ability of the model to reproduce daily rainfall extremes is assessed, again by a comparison with extremes from the MIRA dataset. The results suggest that the model reproduces the number and spatial distribution of rainfall extremes with some accuracy, but that mean rainfall and rainfall variability is underestimated (over-estimated) over wet (dry) regions of southern Africa.

Using zero-norm constraint for sparse probability density function estimation

A new sparse kernel probability density function (pdf) estimator based on zero-norm constraint is constructed using the classical Parzen window (PW) estimate as the target function. The so-called zero-norm of the parameters is used in order to achieve enhanced model sparsity, and it is suggested to minimize an approximate function of the zero-norm. It is shown that under certain condition, the kernel weights of the proposed pdf estimator based on the zero-norm approximation can be updated using the multiplicative nonnegative quadratic programming algorithm. Numerical examples are employed to demonstrate the efficacy of the proposed approach.

Limitations of the MRL mouse as a model for cardiac regeneration

Objective Myocardial repair following injury in mammals is restricted such that damaged areas are replaced by scar tissue, impairing cardiac function. MRL mice exhibit exceptional regenerative healing in an ear punch wound model. Some myocardial repair with restoration of heart function has also been reported following cryoinjury. Increased cardiomyocyte proliferation and a foetal liver stem cell population were implicated. We investigated molecular mechanisms facilitating myocardial repair in MRL mice to identify potential therapeutic targets in non-regenerative species. Methods Expressions of specific cell-cycle regulators that might account for regeneration (CDKs 1, 2, 4 and 6; cyclins A, E, D1 and B1; p21, p27 and E2F5) were compared by immunoblotting in MRL and control C57BL/6 ventricles during development. Flow cytometry was used to investigate stem cell populations in livers from foetal mice, and infarct sizes were compared in coronary artery-ligated and sham-treated MRL and C57BL/6 adult mice. Key findings No differences in the expressions of cell cycle regulators were observed between the two strains. Expressions of CD34+Sca1+ckit-, CD34+Sca1+ckit+ and CD34+Sca1-ckit+ increased in livers from C57BL/6 vs MRL mice. No differences were observed in infarct sizes, levels of fibrosis, Ki67 staining or cardiac function between MRL and C57BL/6 mice. Conclusions No intrinsic differences were observed in cell cycle control molecules or stem cell populations between MRL and control C57BL mouse hearts. Pathophysiologically relevant ischaemic injury is not repaired more efficiently in MRL myocardium, questioning the use of the MRL mouse as a reliable model for cardiac regeneration in response to pathophysiologically relevant forms of injury.

Automated tertiary structure prediction with accurate local model quality assessment using the IntFOLD-TS method

The IntFOLD-TS method was developed according to the guiding principle that the model quality assessment would be the most critical stage for our template based modelling pipeline. Thus, the IntFOLD-TS method firstly generates numerous alternative models, using in-house versions of several different sequence-structure alignment methods, which are then ranked in terms of global quality using our top performing quality assessment method – ModFOLDclust2. In addition to the predicted global quality scores, the predictions of local errors are also provided in the resulting coordinate files, using scores that represent the predicted deviation of each residue in the model from the equivalent residue in the native structure. The IntFOLD-TS method was found to generate high quality 3D models for many of the CASP9 targets, whilst also providing highly accurate predictions of their per-residue errors. This important information may help to make the 3D models that are produced by the IntFOLD-TS method more useful for guiding future experimental work

A model for exchange rates with crawling bands: an application to the Colombian peso

This paper builds upon previous research on currency bands, and provides a model for the Colombian peso. Stochastic differential equations are combined with information related to the Colombian currency band to estimate competing models of the behaviour of the Colombian peso within the limits of the currency band. The resulting moments of the density function for the simulated returns describe adequately most of the characteristics of the sample returns data. The factor included to account for the intra-marginal intervention performed to drive the rate towards the Central Parity accounts only for 6.5% of the daily change, which supports the argument that intervention, if performed by the Central Bank, it is not directed to push the currency towards the limits. Moreover, the credibility of the Colombian Central Bank, Banco de la República’s ability to defend the band seems low.

The variability of tropical ice cloud properties as a function of the large-scale context from ground-based radar-lidar observations over Darwin, Australia

The high complexity of cloud parameterizations now held in models puts more pressure on observational studies to provide useful means to evaluate them. One approach to the problem put forth in the modelling community is to evaluate under what atmospheric conditions the parameterizations fail to simulate the cloud properties and under what conditions they do a good job. It is the ambition of this paper to characterize the variability of the statistical properties of tropical ice clouds in different tropical "regimes" recently identified in the literature to aid the development of better process-oriented parameterizations in models. For this purpose, the statistical properties of non-precipitating tropical ice clouds over Darwin, Australia are characterized using ground-based radar-lidar observations from the Atmospheric Radiation Measurement (ARM) Program. The ice cloud properties analysed are the frequency of ice cloud occurrence, the morphological properties (cloud top height and thickness), and the microphysical and radiative properties (ice water content, visible extinction, effective radius, and total concentration). The variability of these tropical ice cloud properties is then studied as a function of the large-scale cloud regimes derived from the International Satellite Cloud Climatology Project (ISCCP), the amplitude and phase of the Madden-Julian Oscillation (MJO), and the large-scale atmospheric regime as derived from a long-term record of radiosonde observations over Darwin. The vertical variability of ice cloud occurrence and microphysical properties is largest in all regimes (1.5 order of magnitude for ice water content and extinction, a factor 3 in effective radius, and three orders of magnitude in concentration, typically). 98 % of ice clouds in our dataset are characterized by either a small cloud fraction (smaller than 0.3) or a very large cloud fraction (larger than 0.9). In the ice part of the troposphere three distinct layers characterized by different statistically-dominant microphysical processes are identified. The variability of the ice cloud properties as a function of the large-scale atmospheric regime, cloud regime, and MJO phase is large, producing mean differences of up to a factor 8 in the frequency of ice cloud occurrence between large-scale atmospheric regimes and mean differences of a factor 2 typically in all microphysical properties. Finally, the diurnal cycle of the frequency of occurrence of ice clouds is also very different between regimes and MJO phases, with diurnal amplitudes of the vertically-integrated frequency of ice cloud occurrence ranging from as low as 0.2 (weak diurnal amplitude) to values in excess of 2.0 (very large diurnal amplitude). Modellers should now use these results to check if their model cloud parameterizations are capable of translating a given atmospheric forcing into the correct statistical ice cloud properties.

A novel method for assessing risks to pollinators from plant protection products using honeybees as a model species

Pollination is one of the most important ecosystem services in agroecosystems and supports food production. Pollinators are potentially at risk being exposed to pesticides and the main route of exposure is direct contact, in some cases ingestion, of contaminated materials such as pollen, nectar, flowers and foliage. To date there are no suitable methods for predicting pesticide exposure for pollinators, therefore official procedures to assess pesticide risk are based on a Hazard Quotient. Here we develop a procedure to assess exposure and risk for pollinators based on the foraging behaviour of honeybees (Apis mellifera) and using this species as indicator representative of pollinating insects. The method was applied in 13 European field sites with different climatic, landscape and land use characteristics. The level of risk during the crop growing season was evaluated as a function of the active ingredients used and application regime. Risk levels were primarily determined by the agronomic practices employed (i.e. crop type, pest control method, pesticide use), and there was a clear temporal partitioning of risks through time. Generally the risk was higher in sites cultivated with permanent crops, such as vineyard and olive, than in annual crops, such as cereals and oil seed rape. The greatest level of risk is generally found at the beginning of the growing season for annual crops and later in June–July for permanent crops.

Multiscale evolving complex network model of functional connectivity in neuronal cultures

Cultures of cortical neurons grown on multielectrode arrays exhibit spontaneous, robust and recurrent patterns of highly synchronous activity called bursts. These bursts play a crucial role in the development and topological selforganization of neuronal networks. Thus, understanding the evolution of synchrony within these bursts could give insight into network growth and the functional processes involved in learning and memory. Functional connectivity networks can be constructed by observing patterns of synchrony that evolve during bursts. To capture this evolution, a modelling approach is adopted using a framework of emergent evolving complex networks and, through taking advantage of the multiple time scales of the system, aims to show the importance of sequential and ordered synchronization in network function.

Comparative effects of six probiotic strains on immune function in vitro

There is considerable interest in the strain specificity of immune modulation by probiotics. The present study compared the immunomodulatory properties of six probiotic strains of different species and two genera in a human peripheral blood mononuclear cell (PBMC) model in vitro. Live cells of lactobacilli (Lactobacillus casei Shirota, L. rhamnosus GG, L. plantarum NCIMB 8826 and L. reuteri NCIMB 11951) and bifidobacteria (Bifidobacterium longum SP 07/3 and B. bifidum MF 20/5) were individually incubated with PBMC from seven healthy subjects for 24 h. Probiotic strains increased the proportion of CD69+ on lymphocytes, T cells, T cell subsets and natural killer (NK) cells, and increased the proportion of CD25+, mainly on lymphocytes and NK cells. The effects on activation marker expression did not appear to be strain specific. NK cell activity was significantly increased by all six strains, without any significant difference between strains. Probiotic strains increased production of IL-1β, IL-6, IL-10, TNF-α, granulocyte-macrophage colony-stimulating factor and macrophage inflammatory protein 1α to different extents, but had no effect on the production of IL-2, IL-4, IL-5 or TNF-β. The cytokines that showed strain-specific modulation included IL-10, interferon-γ, TNF-α, IL-12p70, IL-6 and monocyte chemotactic protein-1. The Lactobacillus strains tended to promote T helper 1 cytokines, whereas bifidobacterial strains tended to produce a more anti-inflammatory profile. The results suggest that there was limited evidence of strain-specific effects of probiotics with respect to T cell and NK cell activation or NK cell activity, whereas production of some cytokines was differentially influenced by probiotic strains.

Investigation of the factors influencing the survival of Bifidobacterium longum in model acidic solutions and fruit juices

The survival of Bifidobacterium longum NCIMB 8809 was studied during refrigerated storage for 6 weeks in model solutions, based on which a mathematical model was constructed describing cell survival as a function of pH, citric acid, protein and dietary fibre. A Central Composite Design (CCD) was developed studying the influence of four factors at three levels, i.e., pH (3.2–4), citric acid (2–15 g/l), protein (0–10 g/l), and dietary fibre (0–8 g/l). In total, 31 experimental runs were carried out. Analysis of variance (ANOVA) of the regression model demonstrated that the model fitted well the data. From the regression coefficients it was deduced that all four factors had a statistically significant (P < 0.05) negative effect on the log decrease [log10N0 week−log10N6 week], with the pH and citric acid being the most influential ones. Cell survival during storage was also investigated in various types of juices, including orange, grapefruit, blackcurrant, pineapple, pomegranate and strawberry. The highest cell survival (less than 0.4 log decrease) after 6 weeks of storage was observed in orange and pineapple, both of which had a pH of about 3.8. Although the pH of grapefruit and blackcurrant was similar (pH ∼3.2), the log decrease of the former was ∼0.5 log, whereas of the latter was ∼0.7 log. One reason for this could be the fact that grapefruit contained a high amount of citric acid (15.3 g/l). The log decrease in pomegranate and strawberry juices was extremely high (∼8 logs). The mathematical model was able to predict adequately the cell survival in orange, grapefruit, blackcurrant, and pineapple juices. However, the model failed to predict the cell survival in pomegranate and strawberry, most likely due to the very high levels of phenolic compounds in these two juices.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

The platelet-surface thiol isomerase enzyme ERp57 modulates platelet function

Background: Thiol isomerases are a family of endoplasmic reticulum enzymes which orchestrate redox-based modifications of protein disulphide bonds. Previous studies have identified important roles for the thiol isomerases PDI and ERp5 in the regulation of normal platelet function. Objectives: Recently, we demonstrated the presence of a further five thiol isomerases at the platelet surface. In this report we aim to report the role of one of these enzymes - ERp57 in the regulation of platelet function. Methods/Results: Using enzyme activity function blocking antibodies, we demonstrate a role for ERp57 in platelet aggregation, dense granule secretion, fibrinogen binding, calcium mobilisation and thrombus formation under arterial conditions. In addition to the effects of ERp57 on isolated platelets, we observe the presence of ERp57 in the developing thrombus in vivo. Furthermore the inhibition of ERp57 function was found to reduce laser-injury induced arterial thrombus formation in a murine model of thrombosis. Conclusions: These data suggest that ERp57 is important for normal platelet function and opens up the possibility that the regulation of platelet function by a range of cell surface thiol isomerases may represent a broad paradigm for the regulation of haemostasis and thrombosis.