Optimal In Silico Target Gene Deletion through Nonlinear Programming for Genetic Engineering

Optimal selection of multiple regulatory genes, known as targets, for deletion to enhance or suppress the activities of downstream genes or metabolites is an important problem in genetic engineering. Such problems become more feasible to address in silico due to the availability of more realistic dynamical system models of gene regulatory and metabolic networks. The goal of the computational problem is to search for a subset of genes to knock out so that the activity of a downstream gene or a metabolite is optimized.Based on discrete dynamical system modeling of gene regulatory networks, an integer programming problem is formulated for the optimal in silico target gene deletion problem. In the first result, the integer programming problem is proved to be NP-hard and equivalent to a nonlinear programming problem. In the second result, a heuristic algorithm, called GKONP, is designed to approximate the optimal solution, involving an approach to prune insignificant terms in the objective function, and the parallel differential evolution algorithm. In the third result, the effectiveness of the GKONP algorithm is demonstrated by applying it to a discrete dynamical system model of the yeast pheromone pathways. The empirical accuracy and time efficiency are assessed in comparison to an optimal, but exhaustive search strategy.Although the in silico target gene deletion problem has enormous potential applications in genetic engineering, one must overcome the computational challenge due to its NP-hardness. The presented solution, which has been demonstrated to approximate the optimal solution in a practical amount of time, is among the few that address the computational challenge. In the experiment on the yeast pheromone pathways, the identified best subset of genes for deletion showed advantage over genes that were selected empirically. Once validated in vivo, the optimal target genes are expected to achieve higher genetic engineering effectiveness than a trial-and-error procedure

MetNetAPI: A flexible method to access and manipulate biological network data from MetNet

<p>Abstract</p> <p>Background</p> <p>Convenient programmatic access to different biological databases allows automated integration of scientific knowledge. Many databases support a function to download files or data snapshots, or a webservice that offers "live" data. However, the functionality that a database offers cannot be represented in a static data download file, and webservices may consume considerable computational resources from the host server.</p> <p>Results</p> <p>MetNetAPI is a versatile Application Programming Interface (API) to the MetNetDB database. It abstracts, captures and retains operations away from a biological network repository and website. A range of database functions, previously only available online, can be immediately (and independently from the website) applied to a dataset of interest. Data is available in four layers: molecular entities, localized entities (linked to a specific organelle), interactions, and pathways. Navigation between these layers is intuitive (e.g. one can request the molecular entities in a pathway, as well as request in what pathways a specific entity participates). Data retrieval can be customized: Network objects allow the construction of new and integration of existing pathways and interactions, which can be uploaded back to our server. In contrast to webservices, the computational demand on the host server is limited to processing data-related queries only.</p> <p>Conclusions</p> <p>An API provides several advantages to a systems biology software platform. MetNetAPI illustrates an interface with a central repository of data that represents the complex interrelationships of a metabolic and regulatory network. As an alternative to data-dumps and webservices, it allows access to a current and "live" database and exposes analytical functions to application developers. Yet it only requires limited resources on the server-side (thin server/fat client setup). The API is available for Java, Microsoft.NET and R programming environments and offers flexible query and broad data- retrieval methods. Data retrieval can be customized to client needs and the API offers a framework to construct and manipulate user-defined networks. The design principles can be used as a template to build programmable interfaces for other biological databases. The API software and tutorials are available at <url>http://www.metnetonline.org/api</url>.</p

A combined model reduction algorithm for controlled biochemical systems

Background: Systems Biology continues to produce increasingly large models of complex biochemical reaction networks. In applications requiring, for example, parameter estimation, the use of agent-based modelling approaches, or real-time simulation, this growing model complexity can present a significant hurdle. Often, however, not all portions of a model are of equal interest in a given setting. In such situations methods of model reduction offer one possible approach for addressing the issue of complexity by seeking to eliminate those portions of a pathway that can be shown to have the least effect upon the properties of interest. Methods: In this paper a model reduction algorithm bringing together the complementary aspects of proper lumping and empirical balanced truncation is presented. Additional contributions include the development of a criterion for the selection of state-variable elimination via conservation analysis and use of an ‘averaged’ lumping inverse. This combined algorithm is highly automatable and of particular applicability in the context of ‘controlled’ biochemical networks. Results: The algorithm is demonstrated here via application to two examples; an 11 dimensional model of bacterial chemotaxis in Escherichia coli and a 99 dimensional model of extracellular regulatory kinase activation (ERK) mediated via the epidermal growth factor (EGF) and nerve growth factor (NGF) receptor pathways. In the case of the chemotaxis model the algorithm was able to reduce the model to 2 state-variables producing a maximal relative error between the dynamics of the original and reduced models of only 2.8% whilst yielding a 26 fold speed up in simulation time. For the ERK activation model the algorithm was able to reduce the system to 7 state-variables, incurring a maximal relative error of 4.8%, and producing an approximately 10 fold speed up in the rate of simulation. Indices of controllability and observability are additionally developed and demonstrated throughout the paper. These provide insight into the relative importance of individual reactants in mediating a biochemical system’s input-output response even for highly complex networks. Conclusions: Through application, this paper demonstrates that combined model reduction methods can produce a significant simplification of complex Systems Biology models whilst retaining a high degree of predictive accuracy. In particular, it is shown that by combining the methods of proper lumping and empirical balanced truncation it is often possible to produce more accurate reductions than can be obtained by the use of either method in isolation

Aldo Keto Reductase 1B7 and Prostaglandin F2α Are Regulators of Adrenal Endocrine Functions

Prostaglandin F2α (PGF2α), represses ovarian steroidogenesis and initiates parturition in mammals but its impact on adrenal gland is unknown. Prostaglandins biosynthesis depends on the sequential action of upstream cyclooxygenases (COX) and terminal synthases but no PGF2α synthases (PGFS) were functionally identified in mammalian cells. In vitro, the most efficient mammalian PGFS belong to aldo-keto reductase 1B (AKR1B) family. The adrenal gland is a major site of AKR1B expression in both human (AKR1B1) and mouse (AKR1B3, AKR1B7). Thus, we examined the PGF2α biosynthetic pathway and its functional impact on both cortical and medullary zones. Both compartments produced PGF2α but expressed different biosynthetic isozymes. In chromaffin cells, PGF2α secretion appeared constitutive and correlated to continuous expression of COX1 and AKR1B3. In steroidogenic cells, PGF2α secretion was stimulated by adrenocorticotropic hormone (ACTH) and correlated to ACTH-responsiveness of both COX2 and AKR1B7/B1. The pivotal role of AKR1B7 in ACTH-induced PGF2α release and functional coupling with COX2 was demonstrated using over- and down-expression in cell lines. PGF2α receptor was only detected in chromaffin cells, making medulla the primary target of PGF2α action. By comparing PGF2α-responsiveness of isolated cells and whole adrenal cultures, we demonstrated that PGF2α repressed glucocorticoid secretion by an indirect mechanism involving a decrease in catecholamine release which in turn decreased adrenal steroidogenesis. PGF2α may be regarded as a negative autocrine/paracrine regulator within a novel intra-adrenal feedback loop. The coordinated cell-specific regulation of COX2 and AKR1B7 ensures the generation of this stress-induced corticostatic signal

Long Chain Polyunsaturated Fatty Acid Supplementation in Infancy Reduces Heart Rate and Positively Affects Distribution of Attention

A double-blind, randomized, controlled, parallel-group prospective trial was conducted to determine whether a dose-response existed for four different levels of docosahexaenoic acid (DHA) supplementation on the cognitive performance of infants. A total of 122 term infants were fed one of four different formulas varying in their DHA composition (0.00%, 0.32%, 0.64% and 0.96% of total fatty acids as DHA) from birth to 12 months. The three DHA-supplemented formulas also contained 0.64% of total fatty acids as arachidonic acid (ARA, 20:4n-6). Infants were tested at 4, 6, and 9 months of age on a visual habituation protocol that yielded both behavioral and psychophysiological indices of attention. Infants in all DHA+ARA-supplemented conditions had lower heart rates than those in the unsupplemented condition; there was no dose-response for this effect. The distribution of time that infants spent in different phases of attention (a cognitive index derived from the convergence of behavioral and cardiac responses) varied as a function of dosage. Infants supplemented at the two lower DHA doses spent proportionately more time engaged in active stimulus processing than infants fed the unsupplemented formula, while infants fed the highest dose were intermediate and did not differ from any other group

Loss of Col3a1, the Gene for Ehlers-Danlos Syndrome Type IV, Results in Neocortical Dyslamination

It has recently been discovered that Collagen III, the encoded protein of the type IV Ehlers-Danlos Syndrome (EDS) gene, is one of the major constituents of the pial basement membrane (BM) and serves as the ligand for GPR56. Mutations in GPR56 cause a severe human brain malformation called bilateral frontoparietal polymicrogyria, in which neurons transmigrate through the BM causing severe mental retardation and frequent seizures. To further characterize the brain phenotype of Col3a1 knockout mice, we performed a detailed histological analysis. We observed a cobblestone-like cortical malformation, with BM breakdown and marginal zone heterotopias in Col3a1−/− mouse brains. Surprisingly, the pial BM appeared intact at early stages of development but starting as early as embryonic day (E) 11.5, prominent BM defects were observed and accompanied by neuronal overmigration. Although collagen III is expressed in meningeal fibroblasts (MFs), Col3a1−/− MFs present no obvious defects. Furthermore, the expression and posttranslational modification of α-dystroglycan was undisturbed in Col3a1−/− mice. Based on the previous finding that mutations in COL3A1 cause type IV EDS, our study indicates a possible common pathological pathway linking connective tissue diseases and brain malformations

Different Mechanisms Underlie Post-menarchial Increase in Depression and Weight

Contains fulltext : 99254.pdf (publisher's version ) (Open Access)Background Depression and being overweight are correlated health problems in adulthood. Adolescence is a significant period for the onset and increase of depression and obesity, especially among girls. Pubertal development also occurs with concomitant increases in weight. Thus, it is not yet clear whether the association between depression and being overweight can be explained by pubertal development. Purpose We examined the association between depressive mood, body weight, and pubertal status in adolescent girls. Method The design was cross-sectional. In 962 young adolescent Dutch girls (age range, 11.9-15.9) weight and height measurements were used to calculate height, age, and gender-standardized body weight (zBMI). Questionnaires assessed depressive mood (the Center for Epidemiological Studies-Depression, CES-D, inventory) and menarcheal status (pre or post). Results The correlation between menarcheal status and body weight (r = 0.34, p < 0.001) was not affected by depressive mood, and the correlation between menarcheal status and depressive mood (r = 0.20, p < 0.001) was not affected by body weight. A small correlation between depressive mood and body weight (r = 0.12, p < 0.01) largely disappeared after controlling for menarche. Conclusion Menarcheal status largely explains the association between weight and depression. It is independently associated with both BMI and depression, suggesting that different mechanisms underlie the post-menarcheal increased prevalence of depression and overweight.6 p

Changes in the geographical distribution and abundance of the tick Ixodes ricinus during the past 30 years in Sweden

<p>Abstract</p> <p>Background</p> <p><it>Ixodes ricinus </it>is the main vector in Europe of human-pathogenic Lyme borreliosis (LB) spirochaetes, the tick-borne encephalitis virus (TBEV) and other pathogens of humans and domesticated mammals. The results of a previous 1994 questionnaire, directed at people living in Central and North Sweden (Svealand and Norrland) and aiming to gather information about tick exposure for humans and domestic animals, suggested that <it>Ixodes ricinus </it>ticks had become more widespread in Central Sweden and the southern part of North Sweden from the early 1980s to the early 1990s. To investigate whether the expansion of the tick's northern geographical range and the increasing abundance of ticks in Sweden were still occurring, in 2009 we performed a follow-up survey 16 years after the initial study.</p> <p>Methods</p> <p>A questionnaire similar to the one used in the 1994 study was published in Swedish magazines aimed at dog owners, home owners, and hunters. The questionnaire was published together with a popular science article about the tick's biology and role as a pathogen vector in Sweden. The magazines were selected to get information from people familiar with ticks and who spend time in areas where ticks might be present.</p> <p>Results</p> <p>Analyses of data from both surveys revealed that during the near 30-year period from the early 1980s to 2008, <it>I. ricinus </it>has expanded its distribution range northwards. In the early 1990s ticks were found in new areas along the northern coastline of the Baltic Sea, while in the 2009 study, ticks were reported for the first time from many locations in North Sweden. This included locations as far north as 66°N and places in the interior part of North Sweden. During this 16-year period the tick's range in Sweden was estimated to have increased by 9.9%. Most of the range expansion occurred in North Sweden (north of 60°N) where the tick's coverage area doubled from 12.5% in the early 1990s to 26.8% in 2008. Moreover, according to the respondents, the abundance of ticks had increased markedly in LB- and TBE-endemic areas in South (Götaland) and Central Sweden.</p> <p>Conclusions</p> <p>The results suggest that <it>I. ricinus </it>has expanded its range in North Sweden and has become distinctly more abundant in Central and South Sweden during the last three decades. However, in the northern mountain region <it>I. ricinus </it>is still absent. The increased abundance of the tick can be explained by two main factors: First, the high availability of large numbers of important tick maintenance hosts, i.e., cervids, particularly roe deer (<it>Capreolus capreolus</it>) during the last three decades. Second, a warmer climate with milder winters and a prolonged growing season that permits greater survival and proliferation over a larger geographical area of both the tick itself and deer. High reproductive potential of roe deer, high tick infestation rate and the tendency of roe deer to disperse great distances may explain the range expansion of <it>I. ricinus </it>and particularly the appearance of new TBEV foci far away from old TBEV-endemic localities. The geographical presence of LB in Sweden corresponds to the distribution of <it>I. ricinus</it>. Thus, LB is now an emerging disease risk in many parts of North Sweden. Unless countermeasures are undertaken to keep the deer populations, particularly <it>C. capreolus </it>and <it>Dama dama</it>, at the relatively low levels that prevailed before the late 1970s - especially in and around urban areas where human population density is high - by e.g. reduced hunting of red fox (<it>Vulpes vulpes</it>) and lynx (<it>Lynx lynx</it>), the incidences of human LB and TBE are expected to continue to be high or even to increase in Sweden in coming decades.</p

Targeted anti-vascular therapies for ovarian cancer: current evidence

Ovarian cancer presents at advanced stage in around 75% of women, and despite improvements in treatments such as chemotherapy, the 5-year survival from the disease in women diagnosed between 1996 and 1999 in England and Wales was only 36%. Over 80% of patients with advanced ovarian cancer will relapse and despite a good chance of remission from further chemotherapy, they will usually die from their disease. Sequential treatment strategies are employed to maximise quality and length of life but patients eventually become resistant to cytotoxic agents. The expansion in understanding of the molecular biology that characterises cancer cells has led to the rapid development of new agents to target important pathways but the heterogeneity of ovarian cancer biology means that there is no predominant defect. This review attempts to discuss progress to date in tackling a more general target applicable to ovary cancer-angiogenesis