Modeling of the MAPK machinery activation in response to various abiotic and biotic stresses in plants by a system biology approach

Optimization of Elicitation Conditions to Enhance the Production of Potent Metabolite Withanolide from Withania somnifera (L.)

This study aimed at optimizing conditions for increased withanolide production in Withania somnifera. The elicitors used for the foliar spray on the aerial parts of the plant were salicylic acid, jasmonic acid, and chitosan for the enhancement of withanolides in Withania somnifera under different environmental regimes. Three different elicitors, i.e., chitosan, jasmonic acid and salicylic acid, were applied on the plants through foliar route every 15th day for 6 months, and later plants were used for sample preparation. Further, the elicitors were used in different concentration, i.e., jasmonic acid (50, 200 and 400 ppm), chitosan (10, 50 and 100 ppm) and salicylic acid (0.5, 1 and 2 ppm). The elicitors were sprayed on the foliar parts of the plant between 10:00-11:00 a.m. on application days. For elicitor spray, a calibrated sprayer was used. The withanolide A/withaferin A was quantified through HPLC. It was found that in an open environment, maximum withaferin A content, i.e., 0.570 mg/g (DW), was recorded with jasmonic acid (50 ppm) treatment in comparison to control (0.067 mg/g DW). Thus, there was an 8.5-fold increase in the withaferin A content. Maximum withanolide A content of 0.352 mg/g (DW) was recorded when chitosan (50 ppm) was sprayed, while in the control, withanolide A content was recorded to be 0.031 mg/g (DW); thus, chitosan application increased the production of withanolide A by 11.3-fold. Under controlled conditions, maximum withaferin A content of 1.659 mg/g (DW) was recorded when plants were sprayed with chitosan (100 ppm), which was 8.1 times greater than the control content of 0.203 mg/g (DW). Maximum withanolide A content of 0.460 mg/g (DW) was recorded when chitosan (100 ppm) was applied, whereas in the control, withanolide A content was found to be 0.061 mg/g (DW). Thus, foliar spraying of elicitors in very low concentrations can serve as a low-cost, eco-friendly, labor-intensive and elegant alternative approach that can be practiced by farmers for the enhancement, consistent production and improved yield of withanolide A/withaferin A. This can be a suitable way to enhance plant productivity, thus increasing the availability of withanolide A and withaferin A for the health and pharma industry.

p38 MAPK Endogenous Inhibition Improves Neurological Deficits in Global Cerebral Ischemia/Reperfusion Mice

Cerebral ischemia/reperfusion (I/R) injury is a complex pathophysiological process that can lead to neurological function damage and the formation of cerebral infarction. The p38 MAPK pathway has attracted considerable attention in cerebral I/R injury (IRI), but little research has been carried out on its direct role in vivo. In this study, to observe the effects of p38 MAPK endogenous inhibition on cerebral IRI, p38 heterozygous knockdown (p38^KI/+) mice were used. We hypothesized that p38 signaling might be involved in I/R injury and neurological damage reduction and that neurological behavioral deficits improve when p38 MAPK is inhibited. First, we examined the neurological damage and neurological behavioral deficit effects of I/R injury in WT mice. Cerebral I/R injury was induced by the bilateral common carotid artery occlusion (BCCAO) method. The cerebral infarction area and volume were assessed and analyzed by 2,3,5-triphenyltetrazolium chloride (TTC) staining. p38 MAPK and caspase-3 were detected by western blotting. Neuronal apoptosis was measured using TUNEL staining. Neurological deficits were detected by behavioral testing. Furthermore, to assess whether these neuroprotective effects occurred when p38 MAPK was inhibited, p38 heterozygous knockdown (p38^KI/+) mice were used. We found that p38 MAPK endogenous inhibition rescued hippocampal cell apoptosis, reduced ischemic penumbra, and improved neurological behavioral deficits. These findings showed that p38 MAPK endogenous inhibition had a neuroprotective effect on IRI and that p38 MAPK may be a potential therapeutic target for cerebral IRI.

Combined Transcriptome and Metabolome Analysis of Musa nana Laur. Peel Treated With UV-C Reveals the Involvement of Key Metabolic Pathways

An increasing attention is being given to treat fruits with ultraviolet C (UV-C) irradiation to extend shelf-life, senescence, and protection from different diseases during storage. However, the detailed understanding of the pathways and key changes in gene expression and metabolite accumulation related to UV-C treatments are yet to be explored. This study is a first attempt to understand such changes in banana peel irradiated with UV-C. We treated Musa nana Laur. with 0.02 KJ/m2 UV-C irradiation for 0, 4, 8, 12, 15, and 18 days and studied the physiological and quality indicators. We found that UV-C treatment reduces weight loss and decay rate, while increased the accumulation of total phenols and flavonoids. Similarly, our results demonstrated that UV-C treatment increases the activity of defense and antioxidant system related enzymes. We observed that UV-C treatment for 8 days is beneficial for M. nana peels. The peels of M. nana treated with UV-C for 8 days were then subjected to combined transcriptome and metabolome analysis. In total, there were 425 and 38 differentially expressed genes and accumulated metabolites, respectively. We found that UV-C treatment increased the expression of genes in secondary metabolite biosynthesis related pathways. Concomitant changes in the metabolite accumulation were observed. Key pathways that were responsive to UV-C irradiation include flavonoid biosynthesis, phenylpropanoid bios6ynthesis, plant-pathogen interaction, MAPK signaling (plant), and plant hormone signal transduction pathway. We concluded that UV-C treatment imparts beneficial effects on banana peels by triggering defense responses against disease, inducing expression of flavonoid and alkaloid biosynthesis genes, and activating phytohormone and MAPK signaling pathways.

Melatonin increases growth and salt tolerance of Limonium bicolor by improving photosynthetic and antioxidant capacity

BACKGROUND

Soil salinization is becoming an increasingly serious problem worldwide, resulting in cultivated land loss and desertification, as well as having a serious impact on agriculture and the economy. The indoleamine melatonin (N-acetyl-5-methoxytryptamine) has a wide array of biological roles in plants, including acting as an auxin analog and an antioxidant. Previous studies have shown that exogenous melatonin application alleviates the salt-induced growth inhibition in non-halophyte plants; however, to our knowledge, melatonin effects have not been examined on halophytes, and it is unclear whether melatonin provides similar protection to salt-exposed halophytic plants.

RESULTS

We exposed the halophyte Limonium bicolor to salt stress (300 mM) and concomitantly treated the plants with 5 μM melatonin to examine the effect of melatonin on salt tolerance. Exogenous melatonin treatment promoted the growth of L. bicolor under salt stress, as reflected by increasing its fresh weight and leaf area. This increased growth was caused by an increase in net photosynthetic rate and water use efficiency. Treatment of salt-stressed L. bicolor seedlings with 5 μM melatonin also enhanced the activities of antioxidants (superoxide dismutase [SOD], peroxidase [POD], catalase [CAT], and ascorbate peroxidase [APX]), while significantly decreasing the contents of hydrogen peroxide (H₂O₂), superoxide anion (O₂^•-), and malondialdehyde (MDA). To screen for L. bicolor genes involved in the above physiological processes, high-throughput RNA sequencing was conducted. A gene ontology enrichment analysis indicated that genes related to photosynthesis, reactive oxygen species scavenging, the auxin-dependent signaling pathway and mitogen-activated protein kinase (MAPK) were highly expressed under melatonin treatment. These data indicated that melatonin improved photosynthesis, decreased reactive oxygen species (ROS) and activated MAPK-mediated antioxidant responses, triggering a downstream MAPK cascade that upregulated the expression of antioxidant-related genes. Thus, melatonin improves the salt tolerance of L. bicolor by increasing photosynthesis and improving cellular redox homeostasis under salt stress.

CONCLUSIONS

Our results showed that melatonin can upregulate the expression of genes related to photosynthesis, reactive oxygen species scavenging and mitogen-activated protein kinase (MAPK) of L. bicolor under salt stress, which can improve photosynthesis and antioxidant enzyme activities. Thus melatonin can promote the growth of the species and maintain the homeostasis of reactive oxygen species to alleviate salt stress.

Computational analysis of microarray data of Arabidopsis thaliana challenged with Alternaria brassicicola for identification of key genes in Brassica

Background

Alternaria blight, a recalcitrant disease caused by Alternaria brassicae and Alternaria brassicicola, has been recognized for significant losses of oilseed crops especially rapeseed-mustard throughout the world. Till date, no resistance source is available against the disease; hence, plant breeding methods cannot be used to develop disease-resistant varieties. Therefore, in the present study, efforts have been made to identify resistance and defense-related genes as well as key components of JA-SA-ET-mediated pathway involved in resistance against Alternaria brasscicola through computational analysis of microarray data and network biology approach. Microarray profiling data from wild type and mutant Arabidopsis plants challenged with Alternaria brassicicola along with control plant were obtained from the Gene Expression Omnibus (GEO) database. The data analysis, including DEGs extraction, functional enrichment, annotation, and network analysis, was used to identify genes associated with disease resistance and defense response.

Results

A total of 2854 genes were differentially expressed in WT9C9; among them, 1327 genes were upregulated and 1527 genes were downregulated. A total of 1159 genes were differentially expressed in JAM9C9; among them, 809 were upregulated and 350 were downregulated. A total of 2516 genes were differentially expressed in SAM9C9; among them, 1355 were upregulated and 1161 were downregulated. A total of 1567 genes were differentially expressed in ETM9C9; among them, 917 were upregulated and 650 were downregulated. Besides, a total of 2965 genes were differentially expressed in contrast WT24C24; among them, 1510 genes were upregulated and 1455 genes were downregulated. A total of 4598 genes were differentially expressed in JAM24C24; among them, 2201 were upregulated and 2397 were downregulated. A total of 3803 genes were differentially expressed in SAM24C24; among them, 1819 were upregulated and 1984 were downregulated. A total of 4164 genes were differentially expressed in ETM24C24; among them, 1895 were upregulated and 2269 were downregulated. The upregulated genes of Arabidopsis thaliana were mapped and annotated with CDS sequences of Brassica rapa obtained from PlantGDB database. Additionally, PPI network of these genes were constructed to investigate the key components of hormone-mediated pathway involved in resistance during pathogenesis.

Conclusion

The obtained information from present study can be used to engineer resistance to Alternaria blight caused by Alternaria brasscicola through molecular breeding or genetic manipulation-based approaches for improving Brassica oilseed productivity.

Transcriptome analysis of Actinidia chinensis in response to Botryosphaeria dothidea infection

Ripe rot caused by Botryosphaeria dothidea causes extensive production losses in kiwifruit (Actinidia chinensis Planch.). Our previous study showed that kiwifruit variety “Jinyan” is resistant to B. dothidea while “Hongyang” is susceptible. For a comparative analysis of the response of these varieties to B. dothidea infection, we performed a transcriptome analysis by RNA sequencing. A total of 305.24 Gb of clean bases were generated from 36 libraries of which 175.76 Gb was from the resistant variety and 129.48 Gb from the susceptible variety. From the libraries generated, we identified 44,656 genes including 39,041 reference genes, 5,615 novel transcripts, and 13,898 differentially expressed genes (DEGs). Among these, 2,373 potentially defense-related genes linked to calcium signaling, mitogen-activated protein kinase (MAPK), cell wall modification, phytoalexin synthesis, transcription factors, pattern-recognition receptors, and pathogenesis-related proteins may regulate kiwifruit resistance to B. dothidea. DEGs involved in calcium signaling, MAPK, and cell wall modification in the resistant variety were induced at an earlier stage and at higher levels compared with the susceptible variety. Thirty DEGs involved in plant defense response were strongly induced in the resistant variety at all three time points. This study allowed the first comprehensive understanding of kiwifruit transcriptome in response to B. dothidea and may help identify key genes required for ripe rot resistance in kiwifruit.

In silico identification of mimicking molecule(s) triggering von Willebrand factor in human: a molecular drug target for regulating coagulation pathway

Blood coagulation is a complex and dynamic process wherein the body activates its emergency mechanism to stop bleeding and wound healing via the interactions of prothrombotic and antithrombotic agents. von Willebrand factor (VWF) is a complex glycoprotein and initial component of the hemostasis pathway which serves a multipurpose role in blood coagulation process. There are reports of various plants that contain several bioactive compounds possessing properties of inducing blood coagulation directly or indirectly. In the present study, efforts have been made to identify bioactive compounds that may play a significant role in regulation of the coagulation cascade by accelerating VWF and thus enhance the hemostasis process. An antidiuretic peptide drug, Desmopressin, works on VWF and releases them in circulation. Forty-seven compounds from different plant sources were screened through molecular docking, out of which two compounds, Emodin and Peruvianoside II, showed more binding affinity than the reference drug Desmopressin. Emodin and Peruvianoside II showed binding energies -7.2 and -7.0 kcal/mol, respectively, when docked with VWF, whereas Desmopressin displayed less binding energy (-6.9 kcal/mol). Emodin belongs to anthraquinone from Rumex hastasus and Peruvianoside II belongs to flavanone glycosides from Thevetia peruviana. The mimicking potential of top identified molecules with respect to the drug was confirmed through simulation analysis. Besides, the molecular dynamics simulation (MDS) study (for 20 ns) showed that the Peruvianoside II protein complex was energetically more stable than Emodin protein complex. Based on the results, Peruvianoside II possesses great potential and thus may be considered for development of drugs for hemostasis.

Systems biology of seeds: decoding the secret of biochemical seed factories for nutritional security

Seeds serve as biochemical factories of nutrition, processing, bio-energy and storage related important bio-molecules and act as a delivery system to transmit the genetic information to the next generation. The research pertaining towards delineating the complex system of regulation of genes and pathways related to seed biology and nutrient partitioning is still under infancy. To understand these, it is important to know the genes and pathway(s) involved in the homeostasis of bio-molecules. In recent past with the advent and advancement of modern tools of genomics and genetic engineering, multi-layered 'omics' approaches and high-throughput platforms are being used to discern the genes and proteins involved in various metabolic, and signaling pathways and their regulations for understanding the molecular genetics of biosynthesis and homeostasis of bio-molecules. This can be possible by exploring systems biology approaches via the integration of omics data for understanding the intricacy of seed development and nutrient partitioning. These information can be exploited for the improvement of biologically important chemicals for large-scale production of nutrients and nutraceuticals through pathway engineering and biotechnology. This review article thus describes different omics tools and other branches that are merged to build the most attractive area of research towards establishing the seeds as biochemical factories for human health and nutrition.

Augmentation of crop productivity through interventions of omics technologies in India: challenges and opportunities

With the continuous increase in the population of developing countries and decline of natural resources, there is an urgent need to qualitatively and quantitatively augment crop productivity by using new tools and technologies for improvement of agriculturally important traits. The new scientific and technological omics-based approaches have enabled us to deal with several issues and challenges faced by modern agricultural system and provided us novel opportunities for ensuring food and nutritional security. Recent developments in sequencing techniques have made available huge amount of genomic and transcriptomic data on model and cultivated crop plants including Arabidopsis thaliana, Oryza sativa, Triticum aestivum etc. The sequencing data along with other data generated through several omics platforms have significantly influenced the disciplines of crop sciences. Gene discovery and expression profiling-based technologies are offering enormous opportunities to the scientific community which can now apply marker-assisted selection technology to assess and enhance diversity in their collected germplasm, introgress essential traits from new sources and investigate genes that control key traits of crop plants. Utilization of omics science and technologies for crop productivity, protection and management has recently been receiving a lot of attention; the majority of the efforts have been put into signifying the possible applications of various omics technologies in crop plant sciences. This article highlights the background of challenges and opportunities for augmentation of crop productivity through interventions of omics technologies in India.

Regulatory network of genes associated with stimuli sensing, signal transduction and physiological transformation of appressorium in Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is the most destructive disease affecting the rice production (Oryza sativa), with an average global loss of 10-30% per annum. Recent reports have indicated that the fungus also inflicts blast disease on wheat (Triticum aestivum) posing a serious threat to the wheat production. Due to its easily detected infectious process and manoeuvrable genetic manipulation, M. oryzae is considered a model organism for exploring the molecular mechanism underlying fungal pathogenicity during the pathogen-host interaction. M. oryzae utilises an infectious structure called appressorium to breach the host surface by generating high turgor pressure. The appressorium development is induced by physical and chemical cues which are coordinated by the highly conserved cAMP/PKA, MAPK and calcium signalling cascades. Genes involved in the appressorium development have been identified and well studied in M. oryzae, a summary of the working gene network linking stimuli sensing and physiological transformation of appressorium is needed. This review provides a comprehensive discussion regarding the regulatory networks underlying appressorium development with particular emphasis on sensing of appressorium inducing stimuli, signal transduction, transcriptional regulation and the corresponding developmental and physiological responses. We also discussed the crosstalk and interaction of various pathways during the appressorium development.

Modeling of the jasmonate signaling pathway in Arabidopsis thaliana with respect to pathophysiology of Alternaria blight in Brassica

The productivity of Oilseed Brassica, one of the economically important crops of India, is seriously affected by the disease, Alternaria blight. The disease is mainly caused by two major necrotrophic fungi, Alternaria brassicae and Alternaria brassicicola which are responsible for significant yield losses. Till date, no resistant source is available against Alternaria blight, hence plant breeding methods can not be used to develop disease resistant varieties. Jasmonate mediated signalling pathway, which is known to play crucial role during defense response against necrotrophs, could be strengthened in Brassica plants to combat the disease. Since scanty information is available in Brassica-Alternaria pathosystems at molecular level therefore, in the present study efforts have been made to model jasmonic acid pathway in Arabidopsis thaliana to simulate the dynamic behaviour of molecular species in the model. Besides, the developed model was also analyzed topologically for investigation of the hubs node. COI1 is identified as one of the promising candidate genes in response to Alternaria and other linked components of plant defense mechanisms against the pathogens. The findings from present study are therefore informative for understanding the molecular basis of pathophysiology and rational management of Alternaria blight for securing food and nutritional security.

Protective effects of bellidifolin in hypoxia-induced in pheochromocytoma cells (PC12) and underlying mechanisms

Bellidifolin, a xanthone compound derived from plants of Gentiana species, is known to exert a variety of pharmacological activities including anti-oxidation, anti-inflammatory and antitumor actions as well as a protective effect on cerebral ischemic nerve injury. The aim of this study was to examine the protective effects of bellidifolin on nerve injury produced by hypoxia and possible underlying mechanisms using pheochromocytoma cells (PC12). Data showed that the viability of PC12 cells subjected to hypoxia resulted in a significant decrease; however; pretreatment with certain concentrations of bellidifolin (20 or 40 μmol/L) prior to hypoxia significantly increased the survival rate. The results of immunohistochemical staining analysis revealed that there were no marked alterations in the expression of pERK protein between all bellidifolin groups while the expression of p-p38MAPK protein was significantly enhanced by hypoxia. Pretreatment with different concentrations of bellidifolin followed by hypoxia significantly decreased the expression of p-p38MAPK protein. The results of western blot analysis showed that hypoxia induced the expression of the MAPK signaling pathway downstream of the key apoptosis factor caspase-3. Compared to hypoxia, the expression of caspase-3 in the presence of belliidifolin was significantly lower. Data suggest that bellidifolin may contribute to the protective effects associated with nerve injury initiated by hypoxia by mechanisms related to inhibition of cell apoptosis independent of the ERK pathway, but may involve blockade of p38MAPK signaling pathway activation and downstream caspase-3 expression.

In Silico Identification of Mimicking Molecules as Defense Inducers Triggering Jasmonic Acid Mediated Immunity against Alternaria Blight Disease in Brassica Species

Alternaria brassicae and Alternaria brassicicola are two major phytopathogenic fungi which cause Alternaria blight, a recalcitrant disease on Brassica crops throughout the world, which is highly destructive and responsible for significant yield losses. Since no resistant source is available against Alternaria blight, therefore, efforts have been made in the present study to identify defense inducer molecules which can induce jasmonic acid (JA) mediated defense against the disease. It is believed that JA triggered defense response will prevent necrotrophic mode of colonization of Alternaria brassicae fungus. The JA receptor, COI1 is one of the potential targets for triggering JA mediated immunity through interaction with JA signal. In the present study, few mimicking compounds more efficient than naturally occurring JA in terms of interaction with COI1 were identified through virtual screening and molecular dynamics simulation studies. A high quality structural model of COI1 was developed using the protein sequence of Brassica rapa. This was followed by virtual screening of 767 analogs of JA from ZINC database for interaction with COI1. Two analogs viz. ZINC27640214 and ZINC43772052 showed more binding affinity with COI1 as compared to naturally occurring JA. Molecular dynamics simulation of COI1 and COI1-JA complex, as well as best screened interacting structural analogs of JA with COI1 was done for 50 ns to validate the stability of system. It was found that ZINC27640214 possesses efficient, stable, and good cell permeability properties. Based on the obtained results and its physicochemical properties, it is capable of mimicking JA signaling and may be used as defense inducers for triggering JA mediated resistance against Alternaria blight, only after further validation through field trials.

Systems Biology for Smart Crops and Agricultural Innovation: Filling the Gaps between Genotype and Phenotype for Complex Traits Linked with Robust Agricultural Productivity and Sustainability

In recent years, rapid developments in several omics platforms and next generation sequencing technology have generated a huge amount of biological data about plants. Systems biology aims to develop and use well-organized and efficient algorithms, data structure, visualization, and communication tools for the integration of these biological data with the goal of computational modeling and simulation. It studies crop plant systems by systematically perturbing them, checking the gene, protein, and informational pathway responses; integrating these data; and finally, formulating mathematical models that describe the structure of system and its response to individual perturbations. Consequently, systems biology approaches, such as integrative and predictive ones, hold immense potential in understanding of molecular mechanism of agriculturally important complex traits linked to agricultural productivity. This has led to identification of some key genes and proteins involved in networks of pathways involved in input use efficiency, biotic and abiotic stress resistance, photosynthesis efficiency, root, stem and leaf architecture, and nutrient mobilization. The developments in the above fields have made it possible to design smart crops with superior agronomic traits through genetic manipulation of key candidate genes.

The wheat MAP kinase phosphatase 1 alleviates salt stress and increases antioxidant activities in Arabidopsis

Mitogen-activated protein kinase phosphatases (MKPs) are important negative regulators in the MAPK signaling pathways, which play crucial roles in plant growth, development and stress responses. We have previously shown that the heterologous expression of a durum wheat MKP, TMKP1, results in increased tolerance to salt stress in yeast but its particular contribution in salt stress tolerance in plants was not investigated. Here, TMKP1 was overexpressed in Arabidopsis thaliana and physiological changes were assessed in transgenic plants exposed to stress conditions. Under salt stress and especially LiCl, the TMKP1 overexpressors displayed higher germination rates in comparison to wild type plants. The enhancement of salt stress tolerance was accompanied by increased antioxidant enzyme activities, namely superoxide dismutase, catalase and peroxydases. Such increases in antioxidant activities were concomitant with lower malondialdehyde, superoxide anion O2(-) and hydrogen peroxide levels in the TMKP1 transgenic seedlings. Moreover, we provide evidence that, in contrast to the Arabidopsis ortholog AtMKP1, TMKP1 acts as a positive regulator of salt stress tolerance via its ectopic expression in the Arabidopsis mkp1 mutant.

ZBIT Bioinformatics Toolbox: A Web-Platform for Systems Biology and Expression Data Analysis

Bioinformatics analysis has become an integral part of research in biology. However, installation and use of scientific software can be difficult and often requires technical expert knowledge. Reasons are dependencies on certain operating systems or required third-party libraries, missing graphical user interfaces and documentation, or nonstandard input and output formats. In order to make bioinformatics software easily accessible to researchers, we here present a web-based platform. The Center for Bioinformatics Tuebingen (ZBIT) Bioinformatics Toolbox provides web-based access to a collection of bioinformatics tools developed for systems biology, protein sequence annotation, and expression data analysis. Currently, the collection encompasses software for conversion and processing of community standards SBML and BioPAX, transcription factor analysis, and analysis of microarray data from transcriptomics and proteomics studies. All tools are hosted on a customized Galaxy instance and run on a dedicated computation cluster. Users only need a web browser and an active internet connection in order to benefit from this service. The web platform is designed to facilitate the usage of the bioinformatics tools for researchers without advanced technical background. Users can combine tools for complex analyses or use predefined, customizable workflows. All results are stored persistently and reproducible. For each tool, we provide documentation, tutorials, and example data to maximize usability. The ZBIT Bioinformatics Toolbox is freely available at https://webservices.cs.uni-tuebingen.de/.

SBMLsqueezer 2: context-sensitive creation of kinetic equations in biochemical networks

BACKGROUND

The size and complexity of published biochemical network reconstructions are steadily increasing, expanding the potential scale of derived computational models. However, the construction of large biochemical network models is a laborious and error-prone task. Automated methods have simplified the network reconstruction process, but building kinetic models for these systems is still a manually intensive task. Appropriate kinetic equations, based upon reaction rate laws, must be constructed and parameterized for each reaction. The complex test-and-evaluation cycles that can be involved during kinetic model construction would thus benefit from automated methods for rate law assignment.

RESULTS

We present a high-throughput algorithm to automatically suggest and create suitable rate laws based upon reaction type according to several criteria. The criteria for choices made by the algorithm can be influenced in order to assign the desired type of rate law to each reaction. This algorithm is implemented in the software package SBMLsqueezer 2. In addition, this program contains an integrated connection to the kinetics database SABIO-RK to obtain experimentally-derived rate laws when desired.

CONCLUSIONS

The described approach fills a heretofore absent niche in workflows for large-scale biochemical kinetic model construction. In several applications the algorithm has already been demonstrated to be useful and scalable. SBMLsqueezer is platform independent and can be used as a stand-alone package, as an integrated plugin, or through a web interface, enabling flexible solutions and use-case scenarios.

Atlas of Cancer Signalling Network: a systems biology resource for integrative analysis of cancer data with Google Maps

Cancerogenesis is driven by mutations leading to aberrant functioning of a complex network of molecular interactions and simultaneously affecting multiple cellular functions. Therefore, the successful application of bioinformatics and systems biology methods for analysis of high-throughput data in cancer research heavily depends on availability of global and detailed reconstructions of signalling networks amenable for computational analysis. We present here the Atlas of Cancer Signalling Network (ACSN), an interactive and comprehensive map of molecular mechanisms implicated in cancer. The resource includes tools for map navigation, visualization and analysis of molecular data in the context of signalling network maps. Constructing and updating ACSN involves careful manual curation of molecular biology literature and participation of experts in the corresponding fields. The cancer-oriented content of ACSN is completely original and covers major mechanisms involved in cancer progression, including DNA repair, cell survival, apoptosis, cell cycle, EMT and cell motility. Cell signalling mechanisms are depicted in detail, together creating a seamless ‘geographic-like' map of molecular interactions frequently deregulated in cancer. The map is browsable using NaviCell web interface using the Google Maps engine and semantic zooming principle. The associated web-blog provides a forum for commenting and curating the ACSN content. ACSN allows uploading heterogeneous omics data from users on top of the maps for visualization and performing functional analyses. We suggest several scenarios for ACSN application in cancer research, particularly for visualizing high-throughput data, starting from small interfering RNA-based screening results or mutation frequencies to innovative ways of exploring transcriptomes and phosphoproteomes. Integration and analysis of these data in the context of ACSN may help interpret their biological significance and formulate mechanistic hypotheses. ACSN may also support patient stratification, prediction of treatment response and resistance to cancer drugs, as well as design of novel treatment strategies.