Complex Biological Networks Research Articles

Pharmacologic Treatment of Type 2-High Severe Asthma

Asthma is a common and heterogeneous disease whose treatment has considerably changed over the last decades. While inhaled therapies based on inhaled corticosteroids and long acting β2 agonists are effective in controlling asthma in the majority of patients, about 5% of asthmatics poorly respond to inhaled steroids or inhaled steroids/long acting β2 agonists combinations. These patients are affected by “severe asthma”, which is associated with need of oral corticosteroids, progression of the disease, increased use of healthcare services, deterioration of quality of life, and a significant economic burden on society. Asthma is no longer considered a single disease, but a respiratory syndrome with complex biological network of distinct and interrelating inflammatory and remodeling pathways (endotypes) that are associated with different clinical manifestations (phenotypes) both in the lungs (asthma) and other organs (e.g. nose and skin). Severe asthma endotypes may be broadly regarded as Type 2-high and Type 2-low, a model that has become central to asthma management with the development of novel treatments for the Type 2-high endotypes. The hallmark feature of Type 2-high asthma is eosinophilic inflammation, often associated with increased serum IgE, increased exhaled nitric oxide (FeNO) and blood eosinophilia. The discovery of the main key drivers of Type 2-high inflammation (IgE, cytokines such as interleukin IL-5, -4 and -13) enabled the development of new biological agents directed towards specific molecular targets. These advances have shifted the existing paradigm “one drug fits all” to “patient-tailored” novel therapies. The monoclonal antibodies direct to IgE (omalizumab), IL-5 and IL-5 receptor (mepolizumab, benralizumab, reslizumab), and to the α chain of the IL-4 and IL-13 combined receptor (dupilumab), and more recently, to the thymic stromal lymphopoietin (tezepelumab) have been shown in both clinical trials and real-life studies to control symptoms, reduce asthma exacerbations and improve lung function in severe asthmatics not controlled by full inhalation therapies. More recently, the Single Inhaler Triple Therapy (SITT) containing inhaled steroids, long acting β2 agonists and muscarinic antagonists has been developed, slightly improving the effectiveness/safety of the inhalation therapy. This report aims to review available therapeutic opportunities for patients with severe asthma focusing on patients with Type 2-high severe asthma and how to position these new therapeutic alternatives in clinical practice.

Attractor-Transition Control of Complex Biological Networks: A Constant Control Approach.

This article presents attractor-transition control of complex biological networks represented by Boolean networks (BNs) wherein the BN is steered from a prescribed initial attractor toward a desired one. The proposed approach leverages the similarity between attractors and Boolean algebraic properties embedded in the underlying state transition equations. To enhance the clarity of expression regarding stabilization toward the desired attractor, a simple coordinate transformation is performed on the considered BN. Based on the characteristics of transformed state equations, self-stabilizing state variables requiring no control efforts are derived in the first. Next, by applying the feedback vertex set (FVS) control scheme, control inputs stabilizing the remaining state variables are determined. The proposed control scheme exhibits versatility by accommodating both fixed-point and cyclic attractors. We validate the effectiveness of the proposed strategy through extensive numerical experiments conducted on random BNs as well as complex biological systems. In adherence to the reproducible research initiative, detailed results of numerical experiments and all the implementation codes are provided on the authors' website: https://github.com/choonlog/AttractorTransition.

An adaptable in silico ensemble model of the arachidonic acid cascade.

Eicosanoids are a family of bioactive lipids, including derivatives of the ubiquitous fatty acid arachidonic acid (AA). The intimate involvement of eicosanoids in inflammation motivates the development of predictive in silico models for a systems-level exploration of disease mechanisms, drug development and replacement of animal models. Using an ensemble modelling strategy, we developed a computational model of the AA cascade. This approach allows the visualisation of plausible and thermodynamically feasible predictions, overcoming the limitations of fixed-parameter modelling. A quality scoring method was developed to quantify the accuracy of ensemble predictions relative to experimental data, measuring the overall uncertainty of the process. Monte Carlo ensemble modelling was used to quantify the prediction confidence levels. Model applicability was demonstrated using mass spectrometry mediator lipidomics to measure eicosanoids produced by HaCaT epidermal keratinocytes and 46BR.1N dermal fibroblasts, treated with stimuli (calcium ionophore A23187), (ultraviolet radiation, adenosine triphosphate) and a cyclooxygenase inhibitor (indomethacin). Experimentation and predictions were in good qualitative agreement, demonstrating the ability of the model to be adapted to cell types exhibiting differences in AA release and enzyme concentration profiles. The quantitative agreement between experimental and predicted outputs could be improved by expanding network topology to include additional reactions. Overall, our approach generated an adaptable, tuneable ensemble model of the AA cascade that can be tailored to represent different cell types and demonstrated that the integration of in silico and in vitro methods can facilitate a greater understanding of complex biological networks such as the AA cascade.

Construction of a Label-Detection Integrated Visual Probe to Reveal the Double-Edged Sword Principle of Ferroptosis in Liver Injury.

Ferroptosis has been confirmed as a potential mediator and an indicator of the severity of liver injury. Despite the fruitful results, there are still two deficiencies in the research on the association between ferroptosis and liver injury. First, iron ions are usually selected as the target bioanalyte, but its detection based on a fluorescent probe is interfered with by specific chemical reaction mechanisms, leading to low sensitivity and poor physiological stability. Second, more efforts were focused on the harmful effects of ferroptosis on liver injury and less involved in the therapeutic value of ferroptosis for liver injury. Hence, in this work, we proposed a new nonreactive analyte (mitochondrial viscosity) as an analysis marker, which can circumvent the challenges caused by specific reaction mechanisms of iron ions. Meanwhile, we constructed a novel label-detection integrated visual probe (VPF) to explore the feasibility of ferroptosis in the treatment of liver injury. As expected, we not only successfully traced the dynamic changes in mitochondrial viscosity but also visualized the changes in cell morphology during induced and inhibited ferroptosis. Conspicuously, this work revealed that liver injury can be alleviated by regulating ferroptosis, confirming the therapeutic value of ferroptosis in liver injury. In addition, a complex biological communication network between ferroptosis and liver injury was constructed by western blotting, providing an important theoretical mechanism for revealing their double-edged sword relationship. This study not only provides a new strategy for studying the complex relationship between ferroptosis and liver injury but also facilitates the future treatment of liver injury.

Integrated Network Analysis of microRNAs, mRNAs, and Proteins Reveals the Regulatory Interaction between hsa-mir-200b and CFL2 Associated with Advanced Stage and Poor Prognosis in Patients with Intestinal Gastric Cancer.

Intestinal gastric cancer (IGC) carcinogenesis results from a complex interplay between environmental and molecular factors, ultimately contributing to disease development. We used integrative bioinformatic analysis to investigate IGC high-throughput molecular data to uncover interactions among differentially expressed genes, microRNAs, and proteins and their roles in IGC. An integrated network was generated based on experimentally validated microRNA-gene/protein interaction data, with three regulatory circuits involved in a complex network contributing to IGC progression. Key regulators were determined, including 23 microRNA and 15 gene/protein hubs. The regulatory circuit networks were associated with hallmarks of cancer, e.g., cell death, apoptosis and the cell cycle, the immune response, and epithelial-to-mesenchymal transition, indicating that different mechanisms of gene regulation impact similar biological functions. Altered expression of hubs was related to the clinicopathological characteristics of IGC patients and showed good performance in discriminating tumors from adjacent nontumor tissues and in relation to T stage and overall survival (OS). Interestingly, expression of upregulated hub hsa-mir-200b and its downregulated target hub gene/protein CFL2 were related not only to pathological T staging and OS but also to changes during IGC carcinogenesis. Our study suggests that regulation of CFL2 by hsa-miR-200b is a dynamic process during tumor progression and that this control plays essential roles in IGC development. Overall, the results indicate that this regulatory interaction is an important component in IGC pathogenesis. Also, we identified a novel molecular interplay between microRNAs, proteins, and genes associated with IGC in a complex biological network and the hubs closely related to IGC carcinogenesis as potential biomarkers.

Ggkegg: analysis and visualization of KEGG data utilizing grammar of graphics.

The Kyoto Encyclopedia of Genes and Genomes (KEGG) database serves as a valuable systems biology resource and is widely utilized in diverse research fields. However, existing software does not allow flexible visualization and network analyses of the vast and complex KEGG data. We developed ggkegg, an R package that integrates KEGG information with ggplot2. ggkegg enables enhanced visualization and network analyses of KEGG data. We demonstrate the utility of the package by providing examples of its application in single-cell, bulk transcriptome, and microbiome analyses. ggkegg may empower researchers to analyze complex biological networks and present their results effectively. The package and user documentation are available at: https://github.com/noriakis/ggkegg. Supplementary data are available at Bioinformatics online.

Kinetic Modeling and Parameter Estimation of a Prebiotic Peptide Reaction Network.

Although our understanding of how life emerged on Earth from simple organic precursors is speculative, early precursors likely included amino acids. The polymerization of amino acids into peptides and interactions between peptides are of interest because peptides and proteins participate in complex interaction networks in extant biology. However, peptide reaction networks can be challenging to study because of the potential for multiple species and systems-level interactions between species. We developed and employed a computational network model to describe reactions between amino acids to form di-, tri-, and tetra-peptides. Our experiments were initiated with two of the simplest amino acids, glycine and alanine, mediated by trimetaphosphate-activation and drying to promote peptide bond formation. The parameter estimates for bond formation and hydrolysis reactions in the system were found to be poorly constrained due to a network property known as sloppiness. In a sloppy model, the behavior mostly depends on only a subset of parameter combinations, but there is no straightforward way to determine which parameters should be included or excluded. Despite our inability to determine the exact values of specific kinetic parameters, we could make reasonably accurate predictions of model behavior. In short, our modeling has highlighted challenges and opportunities toward understanding the behaviors of complex prebiotic chemical experiments.

Robust discovery of gene regulatory networks from single-cell gene expression data by Causal Inference Using Composition of Transactions.

Gene regulatory networks (GRNs) drive organism structure and functions, so the discovery and characterization of GRNs is a major goal in biological research. However, accurate identification of causal regulatory connections and inference of GRNs using gene expression datasets, more recently from single-cell RNA-seq (scRNA-seq), has been challenging. Here we employ the innovative method of Causal Inference Using Composition of Transactions (CICT) to uncover GRNs from scRNA-seq data. The basis of CICT is that if all gene expressions were random, a non-random regulatory gene should induce its targets at levels different from the background random process, resulting in distinct patterns in the whole relevance network of gene-gene associations. CICT proposes novel network features derived from a relevance network, which enable any machine learning algorithm to predict causal regulatory edges and infer GRNs. We evaluated CICT using simulated and experimental scRNA-seq data in a well-established benchmarking pipeline and showed that CICT outperformed existing network inference methods representing diverse approaches with many-fold higher accuracy. Furthermore, we demonstrated that GRN inference with CICT was robust to different levels of sparsity in scRNA-seq data, the characteristics of data and ground truth, the choice of association measure and the complexity of the supervised machine learning algorithm. Our results suggest aiming at directly predicting causality to recover regulatory relationships in complex biological networks substantially improves accuracy in GRN inference.

Implication of DNA methylation during lifestyle mediated weight loss.

Over the past 50 years, the number of overweight/obese people increased significantly, making obesity a global public health challenge. Apart from rare monogenic forms, obesity is a multifactorial disease, most likely resulting from a concerted interaction of genetic, epigenetic and environmental factors. Although recent studies opened new avenues in elucidating the complex genetics behind obesity, the biological mechanisms contributing to individual's risk to become obese are not yet fully understood. Non-genetic factors such as eating behaviour or physical activity are strong contributing factors for the onset of obesity. These factors may interact with genetic predispositions most likely via epigenetic mechanisms. Epigenome-wide association studies or methylome-wide association studies are measuring DNA methylation at single CpGs across thousands of genes and capture associations to obesity phenotypes such as BMI. However, they only represent a snapshot in the complex biological network and cannot distinguish between causes and consequences. Intervention studies are therefore a suitable method to control for confounding factors and to avoid possible sources of bias. In particular, intervention studies documenting changes in obesity-associated epigenetic markers during lifestyle driven weight loss, make an important contribution to a better understanding of epigenetic reprogramming in obesity. To investigate the impact of lifestyle in obesity state specific DNA methylation, especially concerning the development of new strategies for prevention and individual therapy, we reviewed 19 most recent human intervention studies. In summary, this review highlights the huge potential of targeted interventions to alter disease-associated epigenetic patterns. However, there is an urgent need for further robust and larger studies to identify the specific DNA methylation biomarkers which influence obesity.

An investigation into various visualization tools for complex biological networks

Network biology has become crucial to understanding the complex structural characteristics of biological systems. Consequently, advanced visualization approaches are needed to support the investigation of such structures, and several network visualization tools have subsequently been developed to help researchers analyze intricate biological networks. While these tools support a range of analytical and interactive features, it is sometimes unclear to a data analyst or visualization designer which features are of most relevance to biologists. Thus, this study investigates and identifies essential factors for the visualization of complex biological networks using a mixed methodology approach. Based on the findings, essential factors were categorized as either generic and heuristic, where the former concern different analytical and interactive functionalities, such as an efficient layout, advanced search capabilities, plugin availability, graph analysis and user-friendliness, while the latter concern usability, such as information coding, flexibility, orientation and help.1Furthermore, the findings indicate that 12 of the 15 generic factors identified were moderately important, while all 10 heuristic factors identified herein were moderately important.

Multi-Omics Techniques in Genetic Studies and Breeding of Forest Plants

In recent years, the ecological and economic values of forest plants have been gradually recognized worldwide. However, the growing global demand for new forest plant varieties with higher wood production capacity and better stress tolerance cannot be satisfied by conventional phenotype-based breeding, marker-assisted selection, and genomic selection. In the recent past, diverse omics technologies, including genomics, transcriptomics, epigenomics, proteomics, and metabolomics, have been developed rapidly, providing powerful tools for the precision genetic breeding of forest plants. Genomics lays a solid foundation for understanding complex biological regulatory networks, while other omics technologies provide different perspectives at different levels. Multi-omics integration combines the different omics technologies, becoming a powerful tool for genome-wide functional element identification in forest plant breeding. This review summarizes the recent progress of omics technologies and their applications in the genetic studies on forest plants. It will provide forest plant breeders with an elementary knowledge of multi-omics techniques for future breeding programs.

Topological feature generation for link prediction in biological networks.

Graph or network embedding is a powerful method for extracting missing or potential information from interactions between nodes in biological networks. Graph embedding methods learn representations of nodes and interactions in a graph with low-dimensional vectors, which facilitates research to predict potential interactions in networks. However, most graph embedding methods suffer from high computational costs in the form of high computational complexity of the embedding methods and learning times of the classifier, as well as the high dimensionality of complex biological networks. To address these challenges, in this study, we use the Chopper algorithm as an alternative approach to graph embedding, which accelerates the iterative processes and thus reduces the running time of the iterative algorithms for three different (nervous system, blood, heart) undirected protein-protein interaction (PPI) networks. Due to the high dimensionality of the matrix obtained after the embedding process, the data are transformed into a smaller representation by applying feature regularization techniques. We evaluated the performance of the proposed method by comparing it with state-of-the-art methods. Extensive experiments demonstrate that the proposed approach reduces the learning time of the classifier and performs better in link prediction. We have also shown that the proposed embedding method is faster than state-of-the-art methods on three different PPI datasets.

Random search immune algorithm for community detection

Community detection is a prominent research topic in Complex Network Analysis, and it constitutes an important research field on all those areas where complex networks represent a powerful interpretation tool for describing and understanding systems involved in neuroscience, biology, social science, economy, and many others. A challenging approach to uncover the community structure in complex network, and then revealing the internal organization of nodes, is Modularity optimization. In this research paper, we present an immune optimization algorithm (opt-IA) developed to detect community structures, with the main aim to maximize the modularity produced by the discovered communities. In order to assess the performance of opt-IA, we compared it with an overall of 20 heuristics and metaheuristics, among which one Hyper-Heuristic method, using social and biological complex networks as data set. Unlike these algorithms, opt-IA is entirely based on a fully random search process, which in turn is combined with purely stochastic operators. According to the obtained outcomes, opt-IA shows strictly better performances than almost all heuristics and metaheuristics to which it was compared; whilst it turns out to be comparable with the Hyper-Heuristic method. Overall, it can be claimed that opt-IA, even if driven by a purely random process, proves to be reliable and with efficient performance. Furthermore, to prove the latter claim, a sensitivity analysis of the functionality was conducted, using the classic metrics NMI, ARI and NVI.

On biological networks capable of robust adaptation in the presence of uncertainties: A linear systems-theoretic approach.

Biological adaptation, the tendency of every living organism to regulate its essential activities in environmental fluctuations, is a well-studied functionality in systems and synthetic biology. In this work, we present a generic methodology inspired by systems theory to discover the design principles for robust adaptation, perfect and imperfect, in two different contexts: (1)in the presence of deterministic external and parametric disturbances and (2)in a stochastic setting. In all the cases, firstly, we translate the necessary qualitative conditions for adaptation to mathematical constraints using the language of systems theory, which we then map back as design requirements for the underlying networks. Thus, contrary to the existing approaches, the proposed methodologies provide an exhaustive set of admissible network structures without resorting to computationally burdensome brute-force techniques. Further, the proposed frameworks do not assume prior knowledge about the particular rate kinetics, thereby validating the conclusions for a large class of biological networks. In the deterministic setting, we show that unlike the incoherent feed-forward network structures (IFFLP), the modules containing negative feedback with buffer action (NFBLB) are robust to parametric fluctuations when a specific part of the network is assumed to remain unaffected. To this end, we propose a sufficient condition for imperfect adaptation and show that adding negative feedback in an IFFLP topology improves the robustness concerning parametric fluctuations. Further, we propose a stricter set of necessary conditions for imperfect adaptation. Turning to the stochastic scenario, we adopt a Wiener-Kolmogorov filter strategy to tune the parameters of a given network structure towards minimum output variance. We show that both NFBLB and IFFLP can be used as a reduced-order W-K filter. Further, we define the notion of nearest neighboring motifs to compare the output variances across different network structures. We argue that the NFBLB achieves adaptation at the cost of a variance higher than its nearest neighboring motifs whereas the IFFLP topology produces locally minimum variance while compared with its nearest neighboring motifs. We present numerical simulations to support the theoretical results. Overall, our results present a generic, systematic, and robust framework for advancing the understanding of complex biological networks.

Machine Learning Methodologies for Clustering Gene Expression Data in Cancer

Gene expression data hide vital information required to understand the biological process that takes place in a particular organism. Extracting the hidden patterns in gene expression data helps to strengthen the understanding of functional genomics. The complexity of biological networks and the volume of genes present increase the challenges of comprehending and interpretation of the resulting mass of data, which consists of millions of measurements; these data also inhibit vagueness, imprecision, and noise. Therefore, thousands of gens can be analyzed at a time using clustering techniques is a first step toward addressing these challenges, which is essential in the data mining process to understand natural structures and identify interesting patterns in the underlying gene expression data [2]. The clustering of gene expression data has been proven to be useful in making known the natural structure inherent in gene expression data, understanding gene functions, cellular processes, and subtypes of cells, finding useful information from noisy data, and understanding gene regulation. The other benefit of clustering gene expression data is the identification of homology, which is very important in drug design. Clustering is a useful method that groups items based on certain similarity measures for understanding the structures, functions, regulation of genes, and cellular processes obtained from gene expression data and providing more insight on a given data set [13].
 This review examines the various clustering algorithms applicable to the gene expression data in order to discover and provide useful knowledge of the appropriate clustering technique that will guarantee stability and high degree of accuracy in its analysis procedure.

A linear reciprocal relationship between robustness and plasticity in homeostatic biological networks.

In physics of living systems, a search for relationships of a few macroscopic variables that emerge from many microscopic elements is a central issue. We evolved gene regulatory networks so that the expression of core genes (partial system) is insensitive to environmental changes. Then, we found the expression levels of the remaining genes autonomously increase to provide a plastic (sensitive) response. A feedforward structure from the non-core to core genes evolved autonomously. Negative proportionality was observed between the average changes in core and non-core genes, reflecting reciprocity between the macroscopic robustness of homeostatic genes and plasticity of regulator genes. The proportion coefficient between those genes is represented by their number ratio, as in the "lever principle", whereas the decrease in the ratio results in a transition from perfect to partial adaptation, in which only a portion of the core genes exhibits robustness against environmental changes. This reciprocity between robustness and plasticity was satisfied throughout the evolutionary course, imposing an evolutionary constraint. This result suggests a simple macroscopic law for the adaptation characteristic in evolved complex biological networks.

Dauer fate in a Caenorhabditis elegans Boolean network model.

Cellular fates are determined by genes interacting across large, complex biological networks. A critical question is how to identify causal relationships spanning distinct signaling pathways and underlying organismal phenotypes. Here, we address this question by constructing a Boolean model of a well-studied developmental network and analyzing information flows through the system. Depending on environmental signals Caenorhabditis elegans develop normally to sexual maturity or enter a reproductively delayed, developmentally quiescent 'dauer' state, progressing to maturity when the environment changes. The developmental network that starts with environmental signal and ends in the dauer/no dauer fate involves genes across 4 signaling pathways including cyclic GMP, Insulin/IGF-1, TGF-β and steroid hormone synthesis. We identified three stable motifs leading to normal development, each composed of genes interacting across the Insulin/IGF-1, TGF-β and steroid hormone synthesis pathways. Three genes known to influence dauer fate, daf-2, daf-7 and hsf-1, acted as driver nodes in the system. Using causal logic analysis, we identified a five gene cyclic subgraph integrating the information flow from environmental signal to dauer fate. Perturbation analysis showed that a multifactorial insulin profile determined the stable motifs the system entered and interacted with daf-12 as the switchpoint driving the dauer/no dauer fate. Our results show that complex organismal systems can be distilled into abstract representations that permit full characterization of the causal relationships driving developmental fates. Analyzing organismal systems from this perspective of logic and function has important implications for studies examining the evolution and conservation of signaling pathways.

Proteins as fuzzy controllers: Auto tuning a biological fuzzy inference system to predict protein dynamics in complex biological networks

Biological systems such as mammalian cell cycle are complex systems consisting of a large number of molecular species interacting in ways that produce complex nonlinear systems dynamics. Discrete models such as Boolean models and continuous models such as Ordinary Differential Equations (ODEs) have been widely used to study these systems. Boolean models are simple and can capture qualitative systems behaviour, but they cannot capture the continuous trends of protein concentrations, while ODE models capture continuous trends but require kinetics parameters that are limited. Further, as systems get larger, complexity of these models becomes an issue for parameterization, analysis and interpretation. Also, molecular systems operate under the conditions of uncertainty and noise and our understanding of molecular processes in general is more at a qualitative level characterised by vagueness, imprecision and ambiguity. Hence, as more data are generated, there is a greater need for simpler data driven methods that can approximate continuous system behaviour while representing vagueness and ambiguity without requiring kinetic parameters. Fuzzy inferencing is one such promising method with the ability to work with qualitative vague/imprecise biological knowledge. In this study, we propose a fuzzy inference system for representing continuous behaviour of proteins and apply to some key proteins in the mammalian cell cycle system. The methods we introduced here is novel to protein interaction systems and cell cycle proteins. Our study proposes a three-stage approach to develop fuzzy protein controllers. In stage one, protein system is studied for interactions. We studied some significant core controllers of mammalian cell cycle and their producers and degraders as presented in a published ODE model. Based on the observations from a dataset generated from it, we developed Fuzzy inference systems (FIS) in the second stage, that involved deriving fuzzy IF-THEN rules and their processing, and manually tuned the FIS to predict the dynamics of individual proteins. In stage three, we employed Particle Swarm Optimisation (PSO) for optimising the FIS to further enhance prediction accuracy. Systems dynamics simulation results of the optimised FIS models were in close agreement with the benchmark ODE model results. The results show that the FIS models provide a close approximation to the comprehensive benchmark model in robustly representing continuous protein dynamics while representing the control of protein behavior in an intuitive and transparent format without requiring kinetic parameters. Therefore, FIS models can be an alternative to ODEs in network modelling. Further, FIS models can be assembled to develop large complex systems without losing information or accuracy.

Drug screening and biomarker gene investigation in cancer therapy through the human transcriptional regulatory network

A complex and vast biological network regulates all biological functions in the human body in a sophisticated manner, and abnormalities in this network can lead to disease and even cancer. The construction of a high-quality human molecular interaction network is possible with the development of experimental techniques that facilitate the interpretation of the mechanisms of drug treatment for cancer. We collected 11 molecular interaction databases based on experimental sources and constructed a human protein-protein interaction (PPI) network and a human transcriptional regulatory network (HTRN). A random walk-based graph embedding method was used to calculate the diffusion profiles of drugs and cancers, and a pipeline was constructed by using five similarity comparison metrics combined with a rank aggregation algorithm, which can be implemented for drug screening and biomarker gene prediction. Taking NSCLC as an example, curcumin was identified as a potentially promising anticancer drug from 5450 natural small molecules, and combined with differentially expressed genes, survival analysis, and topological ranking, we obtained BIRC5 (survivin), which is both a biomarker for NSCLC and a key target for curcumin. Finally, the binding mode of curcumin and survivin was explored using molecular docking. This work has a guiding significance for antitumor drug screening and the identification of tumor markers.

EPURE_JSBML: A Tool for Constructing a Deterministic Model of a Reconstituted Escherichia coli Protein Translation System with a User-Specified Nucleic Acid Sequence.

A protein synthesis system is one of the most important and complex biological networks, which translates DNA-encoded information into specific functions. Here, ePURE_JSBML, a tool for constructing biologically relevant large-scale and detailed computational models based on a reconstituted cell-free protein synthesis system, is presented; the user can specify the mRNA sequence, initial component concentration, and decoding rule. Model construction is based on Systems Biology Markup Language (SBML) using JSBML, a pure Java programming library. The tool generates simulation files, executable with Matlab, that enable a variety of simulation experiments including the synthesis of proteins of a few hundred residues.