S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Knowledge management for systems biology a general and visually driven framework applied to translational medicine

Background: To enhance our understanding of complex biological systems like diseases we need to put all of the available data into context and use this to detect relations, pattern and rules which allow predictive hypotheses to be defined. Life science has become a data rich science with information about the behaviour of millions of entities like genes, chemical compounds, diseases, cell types and organs, which are organised in many different databases and/or spread throughout the literature. Existing knowledge such as genotype - phenotype relations or signal transduction pathways must be semantically integrated and dynamically organised into structured networks that are connected with clinical and experimental data. Different approaches to this challenge exist but so far none has proven entirely satisfactory. Results: To address this challenge we previously developed a generic knowledge management framework, BioXM™, which allows the dynamic, graphic generation of domain specific knowledge representation models based on specific objects and their relations supporting annotations and ontologies. Here we demonstrate the utility of BioXM for knowledge management in systems biology as part of the EU FP6 BioBridge project on translational approaches to chronic diseases. From clinical and experimental data, text-mining results and public databases we generate a chronic obstructive pulmonary disease (COPD) knowledge base and demonstrate its use by mining specific molecular networks together with integrated clinical and experimental data. Conclusions: We generate the first semantically integrated COPD specific public knowledge base and find that for the integration of clinical and experimental data with pre-existing knowledge the configuration based set-up enabled by BioXM reduced implementation time and effort for the knowledge base compared to similar systems implemented as classical software development projects. The knowledgebase enables the retrieval of sub-networks including protein-protein interaction, pathway, gene - disease and gene - compound data which are used for subsequent data analysis, modelling and simulation. Pre-structured queries and reports enhance usability; establishing their use in everyday clinical settings requires further simplification with a browser based interface which is currently under development.

FANCM regulates DNA chain elongation and is stabilized by S-phase checkpoint signalling.

FANCM binds and remodels replication fork structures in vitro. We report that in vivo, FANCM controls DNA chain elongation in an ATPase-dependent manner. In the presence of replication inhibitors that do not damage DNA, FANCM counteracts fork movement, possibly by remodelling fork structures. Conversely, through damaged DNA, FANCM promotes replication and recovers stalled forks. Hence, the impact of FANCM on fork progression depends on the underlying hindrance. We further report that signalling through the checkpoint effector kinase Chk1 prevents FANCM from degradation by the proteasome after exposure to DNA damage. FANCM also acts in a feedback loop to stabilize Chk1. We propose that FANCM is a ringmaster in the response to replication stress by physically altering replication fork structures and by providing a tight link to S-phase checkpoint signalling.

The peroxiredoxin Tpx1 is essential as a H202-scavenger during aerobic growth in fission yeast

Peroxiredoxins are known to interact with hydrogen peroxide (H2O2) and to participate in oxidant scavenging, redox signal transduction, and heat-shock responses. The two-cysteine peroxiredoxin Tpx1 of Schizosaccharomyces pombe has been characterized as the H2O2 sensor that transduces the redox signal to the transcription factor Pap1. Here, we show that Tpx1 is essential for aerobic, but not anaerobic, growth. We demonstrate that Tpx1 has an exquisite sensitivity for its substrate, which explains its participation in maintaining low steady-state levels of H2O2. We also show in vitro and in vivo that inactivation of Tpx1 by oxidation of its catalytic cysteine to a sulfinic acid is always preceded by a sulfinic acid form in a covalently linked dimer, which may be important for understanding the kinetics of Tpx1 inactivation. Furthermore, we provide evidence that a strain expressing Tpx1.C169S, lacking the resolving cysteine, can sustain aerobic growth, and we show that small reductants can modulate the activity of the mutant protein in vitro, probably by supplying a thiol group to substitute for cysteine 169.

MyD88, an adapter protein involved in interleukin-1 signaling.

MyD88 has a modular organization, an N-terminal death domain (DD) related to the cytoplasmic signaling domains found in many members of the tumor necrosis factor receptor (TNF-R) superfamily, and a C-terminal Toll domain similar to that found in the expanding family of Toll/interleukin-1-like receptors (IL-1R). This dual domain structure, together with the following observations, supports a role for MyD88 as an adapter in IL-1 signal transduction; MyD88 forms homodimers in vivo through DD-DD and Toll-Toll interactions. Overexpression of MyD88 induces activation of the c-Jun N-terminal kinase (JNK) and the transcription factor NF-kappaB through its DD. A point mutation in MyD88, MyD88-lpr (F56N), which prevents dimerization of the DD, also blocks induction of these activities. MyD88-induced NF-kappaB activation is inhibited by the dominant negative versions of TRAF6 and IRAK, which also inhibit IL-1-induced NF-kappaB activation. Overexpression of MyD88-lpr or MyD88-Toll (expressing only the Toll domain) acted to inhibit IL-1-induced NF-kappaB and JNK activation in a 293 cell line overexpressing the IL-1RI. MyD88 coimmunoprecipitates with the IL-1R signaling complex in an IL-1-dependent manner.

Identification of molecular apocrine breast tumours by microarray analysis.

Previous microarray studies on breast cancer identified multiple tumour classes, of which the most prominent, named luminal and basal, differ in expression of the oestrogen receptor alpha gene (ER). We report here the identification of a group of breast tumours with increased androgen signalling and a 'molecular apocrine' gene expression profile. Tumour samples from 49 patients with large operable or locally advanced breast cancers were tested on Affymetrix U133A gene expression microarrays. Principal components analysis and hierarchical clustering split the tumours into three groups: basal, luminal and a group we call molecular apocrine. All of the molecular apocrine tumours have strong apocrine features on histological examination (P=0.0002). The molecular apocrine group is androgen receptor (AR) positive and contains all of the ER-negative tumours outside the basal group. Kolmogorov-Smirnov testing indicates that oestrogen signalling is most active in the luminal group, and androgen signalling is most active in the molecular apocrine group. ERBB2 amplification is commoner in the molecular apocrine than the other groups. Genes that best split the three groups were identified by Wilcoxon test. Correlation of the average expression profile of these genes in our data with the expression profile of individual tumours in four published breast cancer studies suggest that molecular apocrine tumours represent 8-14% of tumours in these studies. Our data show that it is possible with microarray data to divide mammary tumour cells into three groups based on steroid receptor activity: luminal (ER+ AR+), basal (ER- AR-) and molecular apocrine (ER- AR+).

Nuclear accumulation of the phytochrome A photoreceptor requires FHY1.

The phytochrome family of red/far-red (R/FR)-responsive photoreceptors plays a key role throughout the life cycle of plants . Arabidopsis has five phytochromes, phyA-phyE, among which phyA and phyB play the most predominant functions . Light-regulated nuclear accumulation of the phytochromes is an important regulatory step of this pathway, but to this date no factor specifically required for this event has been identified . Among all phyA signaling mutants, fhy1 and fhy3 (far-red elongated hypocotyl 1 and 3) have the most severe hyposensitive phenotype, indicating that they play particularly important roles . FHY1 is a small plant-specific protein of unknown function localized both in the nucleus and the cytoplasm . Here we show that FHY1 is specifically required for the light-regulated nuclear accumulation of phyA but not phyB. Moreover, phyA accumulation is only slightly affected in fhy3, indicating that the diminished nuclear accumulation of phyA observed in fhy1 seedlings is not simply a general consequence of reduced phyA signaling. By in vitro pull-down and yeast two-hybrid analyses, we demonstrate that FHY1 physically interacts with phyA, preferentially in its active Pfr form. Furthermore, FHY1 and phyA colocalize in planta. We therefore identify the first component required for light-regulated phytochrome nuclear accumulation.

Fas triggers an alternative, caspase-8-independent cell death pathway using the kinase RIP as effector molecule.

Cell death is achieved by two fundamentally different mechanisms: apoptosis and necrosis. Apoptosis is dependent on caspase activation, whereas the caspase-independent necrotic signaling pathway remains largely uncharacterized. We show here that Fas kills activated primary T cells efficiently in the absence of active caspases, which results in necrotic morphological changes and late mitochondrial damage but no cytochrome c release. This Fas ligand-induced caspase-independent death is absent in T cells that are deficient in either Fas-associated death domain (FADD) or receptor-interacting protein (RIP). RIP is also required for necrotic death induced by tumor necrosis factor (TNF) and TNF-related apoptosis-inducing ligand (TRAIL). In contrast to its role in nuclear factor kappa B activation, RIP requires its own kinase activity for death signaling. Thus, Fas, TRAIL and TNF receptors can initiate cell death by two alternative pathways, one relying on caspase-8 and the other dependent on the kinase RIP.

IRF6 is a mediator of Notch pro-differentiation and tumour suppressive function in keratinocytes.

While the pro-differentiation and tumour suppressive functions of Notch signalling in keratinocytes are well established, the underlying mechanisms remain poorly understood. We report here that interferon regulatory factor 6 (IRF6), an IRF family member with an essential role in epidermal development, is induced in differentiation through a Notch-dependent mechanism and is a primary Notch target in keratinocytes and keratinocyte-derived SCC cells. Increased IRF6 expression contributes to the impact of Notch activation on growth/differentiation-related genes, while it is not required for induction of 'canonical' Notch targets like p21(WAF1/Cip1), Hes1 and Hey1. Down-modulation of IRF6 counteracts differentiation of primary human keratinocytes in vitro and in vivo, promoting ras-induced tumour formation. The clinical relevance of these findings is illustrated by the strikingly opposite pattern of expression of Notch1 and IRF6 versus epidermal growth factor receptor in a cohort of clinical SCCs, as a function of their grade of differentiation. Thus, IRF6 is a primary Notch target in keratinocytes, which contributes to the role of this pathway in differentiation and tumour suppression.

Dispensable role of myeloid differentiation primary response gene 88 (MyD88) and MyD88-dependent toll-like receptors (TLRs) in a murine model of osteoarthritis.

OBJECTIVES: The aim of our study was to evaluate the role of cell-membrane expressed TLRs and the signaling molecule MyD88 in a murine model of OA induced by knee menisectomy (surgical partial removal of the medial meniscus [MNX]). METHODS: OA was induced in 8-10weeks old C57Bl/6 wild-type (WT) female (n=7) mice and in knockout (KO) TLR-1 (n=7), -2 (n=8), -4 (n=9) -6 (n=5), MyD88 (n=8) mice by medial menisectomy, using the sham-operated contralateral knee as a control. Cartilage destruction and synovial inflammation were evaluated by knee joint histology using the OARSI scoring method. Apoptotic chondrocytes and cartilage metabolism (collagen II synthesis and MMP-mediated aggrecan degradation) were analyzed using immunohistochemistry. RESULTS: Operated knees exhibited OA features at 8weeks post-surgery compared to sham-operated ones. In menisectomized TLR-1, -2, -4, and -6 deficient mice, cartilage lesions, synovial inflammation and cartilage metabolism were similar to that in operated WT mice. Accordingly, using the same approach, we found no significant protection in MyD88-deficient mice in terms of OA progression as compared to WT littermates. CONCLUSIONS: Deficiency of TLRs or their signalling molecule MyD88 did not impact on the severity of experimental OA. Our results demonstrate that MyD88-dependent TLRs are not involved in this murine OA model. Moreover, the dispensable role of MyD88, which is also an adaptor for IL-1 receptor signaling, suggests that IL-1 is not a key mediator in the development of OA. This latter hypothesis is strengthened by the lack of efficiency of IL-1β antagonist in the treatment of OA.

Alum interaction with dendritic cell membrane lipids is essential for its adjuvanticity.

As an approved vaccine adjuvant for use in humans, alum has vast health implications, but, as it is a crystal, questions remain regarding its mechanism. Furthermore, little is known about the target cells, receptors, and signaling pathways engaged by alum. Here we report that, independent of inflammasome and membrane proteins, alum binds dendritic cell (DC) plasma membrane lipids with substantial force. Subsequent lipid sorting activates an abortive phagocytic response that leads to antigen uptake. Such activated DCs, without further association with alum, show high affinity and stable binding with CD4(+) T cells via the adhesion molecules intercellular adhesion molecule-1 (ICAM-1) and lymphocyte function-associated antigen-1 (LFA-1). We propose that alum triggers DC responses by altering membrane lipid structures. This study therefore suggests an unexpected mechanism for how this crystalline structure interacts with the immune system and how the DC plasma membrane may behave as a general sensor for solid structures.

Cutting edge: a Toll-like receptor 2 polymorphism that is associated with lepromatous leprosy is unable to mediate mycobacterial signaling.

Toll-like receptors (TLRs) are key mediators of the innate immune response to microbial pathogens. We investigated the role of TLRs in the recognition of Mycobacterium leprae and the significance of TLR2Arg(677)Trp, a recently discovered human polymorphism that is associated with lepromatous leprosy. In mice, TNF-alpha production in response to M. leprae was essentially absent in TLR2-deficient macrophages. Similarly, human TLR2 mediated M. leprae-dependent activation of NF-kappaB in transfected Chinese hamster ovary and human embryonic kidney 293 cells, with enhancement of this signaling in the presence of CD14. In contrast, activation of NF-kappaB by human TLR2Arg(677)Trp was abolished in response to M. leprae and Mycobacterium tuberculosis. The impaired function of this TLR2 variant provides a molecular mechanism for the poor cellular immune response associated with lepromatous leprosy and may have important implications for understanding the pathogenesis of other mycobacterial infections.

SV40-induced expression of calretinin protects mesothelial cells from asbestos cytotoxicity and may be a key factor contributing to mesothelioma pathogenesis.

The calcium-binding protein calretinin has emerged as a useful marker for the identification of mesotheliomas of the epithelioid and mixed types, but its putative role in tumor development has not been addressed previously. Although exposure to asbestos fibers is considered the main cause of mesothelioma, undoubtedly, not all mesothelioma patients have a history of asbestos exposure. The question as to whether the SV40 virus is involved as a possible co-factor is still highly debated. Here we show that increased expression of SV40 early gene products in the mesothelial cell line MeT-5A induces the expression of calretinin and that elevated calretinin levels strongly correlate with increased resistance to asbestos cytotoxicity. Calretinin alone mediates a significant part of this protective effect because cells stably transfected with calretinin cDNA were clearly more resistant to the toxic effects of crocidolite than mock-transfected control cells. Down-regulation of calretinin by antisense methods restored the sensitivity to asbestos toxicity to a large degree. The protective effect observed in clones with higher calretinin expression levels could be eliminated by phosphatidylinositol 3-kinase (PI3K) inhibitors, implying an important role for the PI3K/AKT signaling (survival) pathway in mediating the protective effect. Up-regulation of calretinin, resulting from either asbestos exposure or SV40 oncoproteins, may be a common denominator that leads to increased resistance to asbestos cytotoxicity and thereby contributes to mesothelioma carcinogenesis.

Calcineurin signaling as a negative determinant of keratinocyte cancer stem cell potential and carcinogenesis.

Calcineurin is the only known serine-threonine phosphatase under calcium-calmodulin control and key regulator of the immune system. Treatment of patients with calcineurin-inhibitory drugs like cyclosporin A and FK506 to prevent graft rejection dramatically increases the risk of cutaneous squamous cell carcinoma, which is a major cause of death after organ transplants. Recent evidence indicates that suppression of calcineurin signaling, together with its impact on the immune system, exerts direct tumor-promoting effects in keratinocytes, enhancing cancer stem cell potential. The underlying mechanism involves interruption of a double negative regulatory axis, whereby calcineurin and nuclear factors of activated T-cell signaling inhibits expression of ATF3, a negative regulator of p53. The resulting suppression of keratinocyte cancer cell senescence is of likely clinical significance for the many patients under treatment with calcineurin inhibitors and may be of relevance for other cancer types in which altered calcium-calcineurin signaling plays a role.

Syk kinase signalling couples to the Nlrp3 inflammasome for anti-fungal host defence.

Fungal infections represent a serious threat, particularly in immunocompromised patients. Interleukin-1beta (IL-1beta) is a key pro-inflammatory factor in innate antifungal immunity. The mechanism by which the mammalian immune system regulates IL-1beta production after fungal recognition is unclear. Two signals are generally required for IL-1beta production: an NF-kappaB-dependent signal that induces the synthesis of pro-IL-1beta (p35), and a second signal that triggers proteolytic pro-IL-1beta processing to produce bioactive IL-1beta (p17) via Caspase-1-containing multiprotein complexes called inflammasomes. Here we demonstrate that the tyrosine kinase Syk, operating downstream of several immunoreceptor tyrosine-based activation motif (ITAM)-coupled fungal pattern recognition receptors, controls both pro-IL-1beta synthesis and inflammasome activation after cell stimulation with Candida albicans. Whereas Syk signalling for pro-IL-1beta synthesis selectively uses the Card9 pathway, inflammasome activation by the fungus involves reactive oxygen species production and potassium efflux. Genetic deletion or pharmalogical inhibition of Syk selectively abrogated inflammasome activation by C. albicans but not by inflammasome activators such as Salmonella typhimurium or the bacterial toxin nigericin. Nlrp3 (also known as NALP3) was identified as the critical NOD-like receptor family member that transduces the fungal recognition signal to the inflammasome adaptor Asc (Pycard) for Caspase-1 (Casp1) activation and pro-IL-1beta processing. Consistent with an essential role for Nlrp3 inflammasomes in antifungal immunity, we show that Nlrp3-deficient mice are hypersusceptible to Candida albicans infection. Thus, our results demonstrate the molecular basis for IL-1beta production after fungal infection and identify a crucial function for the Nlrp3 inflammasome in mammalian host defence in vivo.