Identifying quantitative operation principles in metabolic pathways: a systematic method for searching feasible enzyme activity patterns leading to cellular adaptive responses

Background: Optimization methods allow designing changes in a system so that specific goals are attained. These techniques are fundamental for metabolic engineering. However, they are not directly applicable for investigating the evolution of metabolic adaptation to environmental changes. Although biological systems have evolved by natural selection and result in well-adapted systems, we can hardly expect that actual metabolic processes are at the theoretical optimum that could result from an optimization analysis. More likely, natural systems are to be found in a feasible region compatible with global physiological requirements. Results: We first present a new method for globally optimizing nonlinear models of metabolic pathways that are based on the Generalized Mass Action (GMA) representation. The optimization task is posed as a nonconvex nonlinear programming (NLP) problem that is solved by an outer- approximation algorithm. This method relies on solving iteratively reduced NLP slave subproblems and mixed-integer linear programming (MILP) master problems that provide valid upper and lower bounds, respectively, on the global solution to the original NLP. The capabilities of this method are illustrated through its application to the anaerobic fermentation pathway in Saccharomyces cerevisiae. We next introduce a method to identify the feasibility parametric regions that allow a system to meet a set of physiological constraints that can be represented in mathematical terms through algebraic equations. This technique is based on applying the outer-approximation based algorithm iteratively over a reduced search space in order to identify regions that contain feasible solutions to the problem and discard others in which no feasible solution exists. As an example, we characterize the feasible enzyme activity changes that are compatible with an appropriate adaptive response of yeast Saccharomyces cerevisiae to heat shock Conclusion: Our results show the utility of the suggested approach for investigating the evolution of adaptive responses to environmental changes. The proposed method can be used in other important applications such as the evaluation of parameter changes that are compatible with health and disease states.

Sink limitations to yield in wheat : how could it be reduced?

Further genetic gains in wheat yield are required to match expected increases in demand. This may require the identification of physiological attributes able to produce such improvement, as well as the genetic bases controlling those traits in order to facilitate their manipulation. In the present paper, a theoretical framework of source and sink limitation to wheat yield is presented and the fine-tuning of crop development as an alternative for increasing yield potential is discussed. Following a top-down approach, most crop physiologists have agreed that the main attribute explaining past genetic gains in yield was harvest index (HI). By virtue of previous success, no further gains may be expected in HI and an alternative must be found. Using a bottom-up approach, the present paper firstly provides evidence on the generalized sink-limited condition of grain growth, determining that for further increases in yield potential, sink strength during grain filling has to be increased. The focus should be on further increasing grain number per m2, through fine-tuning pre-anthesis developmental patterns. The phase of rapid spike growth period (RSGP) is critical for grain number determination and increasing spike growth during pre-anthesis would result in an increased number of grains. This might be achieved by lengthening the duration of the phase (though without altering flowering time), as there is genotypic variation in the proportion of pre-anthesis time elapsed either before or after the onset of the stem elongation phase. Photoperiod sensitivity during RSGP could be then used as a genetic tool to further increase grain number, since slower development results in smoother floret development and more floret primordia achieve the fertile floret stage, able to produce a grain. Far less progress has been achieved on the genetic control of this attribute. None of the well-known major Ppd alleles seems to be consistently responsible for RSGP sensitivity. Alternatives for identifying the genetic factors responsible for this sensitivity (e.g. quantitative trait locus (QTL) identification in mapping populations) are being considered.

A web ontologies framework for digital rights management

In order to improve the management of copyright in the Internet, known as Digital Rights Management, there is the need for a shared language for copyright representation. Current approaches are based on purely syntactic solutions, i.e. a grammar that defines a rights expression language. These languages are difficult to put into practise due to the lack of explicit semantics that facilitate its implementation. Moreover, they are simple from the legal point of view because they are intended just to model the usage licenses granted by content providers to end-users. Thus, they ignore the copyright framework that lies behind and the whole value chain from creators to end-users. Our proposal is to use a semantic approach based on semantic web ontologies. We detail the development of a copyright ontology in order to put this approach into practice. It models the copyright core concepts for creation, rights and the basic kinds of actions that operate on content. Altogether, it allows building a copyright framework for the complete value chain. The set of actions operating on content are our smaller building blocks in order to cope with the complexity of copyright value chains and statements and, at the same time, guarantee a high level of interoperability and evolvability. The resulting copyright modelling framework is flexible and complete enough to model many copyright scenarios, not just those related to the economic exploitation of content. The ontology also includes moral rights, so it is possible to model this kind of situations as it is shown in the included example model for a withdrawal scenario. Finally, the ontology design and the selection of tools result in a straightforward implementation. Description Logic reasoners are used for license checking and retrieval. Rights are modelled as classes of actions, action patterns are modelled also as classes and the same is done for concrete actions. Then, to check if some right or license grants an action is reduced to check for class subsumption, which is a direct functionality of these reasoners.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.

Effects of a post-weaning cafeteria diet in young rats: metabolic syndrome, reduced activity and low anxiety-like behaviour.

Among adolescents, overweight, obesity and metabolic syndrome are rapidly increasing in recent years as a consequence of unhealthy palatable diets. Animal models of diet-induced obesity have been developed, but little is known about the behavioural patterns produced by the consumption of such diets. The aim of the present study was to determine the behavioural and biochemical effects of a cafeteria diet fed to juvenile male and female rats, as well as to evaluate the possible recovery from these effects by administering standard feeding during the last week of the study. Two groups of male and female rats were fed with either a standard chow diet (ST) or a cafeteria (CAF) diet from weaning and for 8 weeks. A third group of males (CAF withdrawal) was fed with the CAF diet for 7 weeks and the ST in the 8th week. Both males and females developed metabolic syndrome as a consequence of the CAF feeding, showing overweight, higher adiposity and liver weight, increased plasma levels of glucose, insulin and triglycerides, as well as insulin resistance, in comparison with their respective controls. The CAF diet reduced motor activity in all behavioural tests, enhanced exploration, reduced anxiety-like behaviour and increased social interaction; this last effect was more pronounced in females than in males. When compared to animals only fed with a CAF diet, CAF withdrawal increased anxiety in the open field, slightly decreased body weight, and completely recovered the liver weight, insulin sensitivity and the standard levels of glucose, insulin and triglycerides in plasma. In conclusion, a CAF diet fed to young animals for 8 weeks induced obesity and metabolic syndrome, and produced robust behavioural changes in young adult rats, whereas CAF withdrawal in the last week modestly increased anxiety, reversed the metabolic alterations and partially reduced overweight.

Phenomenological models of holographic superconductors and hall currents

We study general models of holographic superconductivity parametrized by four arbitrary functions of a neutral scalar field of the bulk theory. The models can accommodate several features of real superconductors, like arbitrary critical temperatures and critical exponents in a certain range, and perhaps impurities or boundary or thickness effects. We find analytical expressions for the critical exponents of the general model and show that they satisfy the Rushbrooke identity. An important subclass of models exhibit second order phase transitions. A study of the specific heat shows that general models can also describe holographic superconductors undergoing first, second and third (or higher) order phase transitions. We discuss how small deformations of the HHH model can lead to the appearance of resonance peaks in the conductivity, which increase in number and become narrower as the temperature is gradually decreased, without the need for tuning mass of the scalar to be close to the Breitenlohner-Freedman bound. Finally, we investigate the inclusion of a generalized ¿theta term¿ producing Hall effect without magnetic field.

Beyond pure public and pure private management models: Partial privatization in the European airport industry

The use of private funding and management is enjoying an increasing trend in airports. The literature has not paid enough attention to the mixed management models in this industry, although many European airports take the form of mixed public-private companies, where ownership is shared between public and private sectors. We examine the determinants of the degree of private participation in the European airport sector. Drawing on a sample of the 100 largest European airports, we estimate a multivariate equation in order to determine the role of airport characteristics, fiscal variables, and political factors on the extent of private involvement. Our results confirm the alignment between public and private interests in partially privatized airports. Fiscal constraints and market attractiveness promote private participation. Integrated governance models and the share of network carriers prevent the presence of private ownership, while the degree of private participation appears to be pragmatic rather than ideological.

Beyond pure public and pure private management models: Mixed firms in the European Airport Industry

The use of private funding and management enjoys an increasing trend in airports. The literature has not paid enough attention to the mixed management models in this industry, although many European airports take the form of mixed firms or Institutional PPP, where ownership is shared between public and private sectors. We examine the determinants of the degree of private participation in the European airport sector. Drawing on a sample of the 100 largest European airports we estimate a multivariate equation in order to determine the role of airport characteristics, fiscal variables and political factors on the extent of private involvement. Our results confirm the alignment between public and private interests in PPPs. Fiscal constraints and market attractiveness promote private participation. Integrated governance models and the share of network carriers prevent the presence of private ownership, while the degree of private participation appears to be pragmatic rather than ideological.

New approximation to the third-order density : application to the calculation of correlated multicenter indices

In this work we present the formulas for the calculation of exact three-center electron sharing indices (3c-ESI) and introduce two new approximate expressions for correlated wave functions. The 3c-ESI uses the third-order density, the diagonal of the third-order reduced density matrix, but the approximations suggested in this work only involve natural orbitals and occupancies. In addition, the first calculations of 3c-ESI using Valdemoro's, Nakatsuji's and Mazziotti's approximation for the third-order reduced density matrix are also presented for comparison. Our results on a test set of molecules, including 32 3c-ESI values, prove that the new approximation based on the cubic root of natural occupancies performs the best, yielding absolute errors below 0.07 and an average absolute error of 0.015. Furthemore, this approximation seems to be rather insensitive to the amount of electron correlation present in the system. This newly developed methodology provides a computational inexpensive method to calculate 3c-ESI from correlated wave functions and opens new avenues to approximate high-order reduced density matrices in other contexts, such as the contracted Schrödinger equation and the anti-Hermitian contracted Schrödinger equation