Search

Hepatitis E, as a zoonotic disease, has been a great challenge to global public health. Therefore, it has important research value and practical significance for the transmission and control of hepatitis E virus (HEV). In the exploration of infectious disease transmission dynamics and optimal control, mathematical models are often applied. Among them, the fractional differential model has become an important and practical tool because of its good memory and genetic characteristics. In this paper, an HEV propagation dynamic model is constructed by the Caputo fractional derivative. First, the properties of the model are analyzed, including the existence, non-negativity, boundedness, and stability of the equilibrium points. Then, from the perspective of fractional optimal control (FOC), control measures were proposed, including improving the awareness and prevention of hepatitis E among susceptible people, strengthening the treatment of infected people, and improving environmental hygiene. Then, an FOC model of HEV was constructed. After analyzing the necessary conditions for optimality, the particle swarm optimization is introduced to optimize the control function. In addition, four control strategies are applied. Finally, the numerical simulation is completed by the fractional Adams–Bashforth–Moulton prediction correction algorithm. The four strategies and no control were compared and analyzed. The numerical simulation results of different fractional orders are also compared and analyzed. The results illustrate that the optimal strategy, compared with no control, reduces the HEV control time by nearly 60 days. Therefore, this method would contribute to the study of HEV transmission dynamics and control mechanisms, thus contributing to the development of global public health. [ABSTRACT FROM AUTHOR]

As a new tumor therapeutic strategy, adaptive therapy involves utilizing the competition between cancer cells to suppress the growth of drug-resistant cells, maintaining a certain tumor burden. However, it is difficult to determine the appropriate time and drug dose. In this paper, we consider the competition model between drug-sensitive cells and drug-resistant cells, propose the problem of drug concentration, and provide two state constraints: the upper limit of the maximum allowable drug concentration and the tumor burden. Using relevant theories, we propose the best treatment strategy. Through a numerical simulation and quantitative analysis, the effects of drug concentrations and different tumor burdens on treatments are studied, and the effects of cell-to-cell competitive advantage on cell changes are taken into account. The clinical dose titration method is further simulated; the results show that our therapeutic regimen can better suppress the growth of drug-resistant cells, control the tumor burden, limit drug toxicity, and extend the effective treatment time. [ABSTRACT FROM AUTHOR]

It is of critical importance to comprehend the biological environment and core tumor populations when trying to design successful therapeutic solutions for fighting cancers. In several diseases, G9a has been recognized as a novel epigenetic therapeutic target, and its blockage can shift tumor cells (TCs) toward tumor propagating cells (TPCs). This study combines mathematical modeling based on ordinary differential equations and dynamical analysis to quantitatively and qualitatively understand the interactions among G9a, TCs, and TPCs, denoted as G9a-TC-TPC. We propose four different dynamical systems with the impact of the strong Allee effect, named the Hill–Hill system, Logistic–Logistic system, Hill–Logistic system and Logistic–Hill system, to simulate different biological processes through the Hill functions and the Logistic functions that are often used in the models of biological systems. Based on theoretical analysis of these models, including the positivity, boundedness and stability of equilibria, we find that the Hill–Logistic system can display bistable states that correspond to the wild-type tumors and the aggressive tumors. Consequently, we use bifurcation analysis and numerical simulations to illustrate the complicated dynamical behavior of this system. It has been shown that under a specific therapy that changes the relative apoptotic rate of TCs (G9a suppresses the apoptosis of TCs), which can affect the bistability and instability of the system, the wide-type state can be obtained. We also discover that the relative handling time of TCs and TPCs can cooperatively enhance bistability, whereas the cooperative coefficient of feedback can contribute to all tumor cells moving from high-level monostability to bistable states in a restricted region, then to low-level stable states. These results offer new insights for more precisely understanding epigenetic therapy treatments with G9a. [ABSTRACT FROM AUTHOR]

Although differential evolution (DE) algorithms perform well on a large variety of complicated optimization problems, only a few theoretical studies are focused on the working principle of DE algorithms. To make the first attempt to reveal the function of binomial crossover, this paper aims to answer whether it can reduce the approximation error of evolutionary algorithms. By investigating the expected approximation error and the probability of not finding the optimum, we conduct a case study comparing two evolutionary algorithms with and without binomial crossover on two classical benchmark problems: OneMax and Deceptive. It is proven that using binomial crossover leads to the dominance of transition matrices. As a result, the algorithm with binomial crossover asymptotically outperforms that without crossover on both OneMax and Deceptive, and outperforms on OneMax, however, not on Deceptive. Furthermore, an adaptive parameter strategy is proposed which can strengthen the superiority of binomial crossover on Deceptive. [ABSTRACT FROM AUTHOR]

With the successful application of single-cell sequencing technology, a large number of single-cell multi-omics sequencing (scMO-seq)data have been generated, which enables researchers to study heterogeneity between individual cells. One prominent problem in single-cell data analysis is the prevalence of dropouts, caused by failures in amplification during the experiments. It is necessary to develop effective approaches for imputing the missing values. Different with general methods imputing single type of single-cell data, we propose an imputation method called scLRTD, using low-rank tensor decomposition based on nuclear norm to impute scMO-seq data and single-cell RNA-sequencing (scRNA-seq)data with different stages, tissues or conditions. Furthermore, four sets of simulated and two sets of real scRNA-seq data from mouse embryonic stem cells and hepatocellular carcinoma, respectively, are used to carry out numerical experiments and compared with other six published methods. Error accuracy and clustering results demonstrate the effectiveness of proposed method. Moreover, we clearly identify two cell subpopulations after imputing the real scMO-seq data from hepatocellular carcinoma. Further, Gene Ontology identifies 7 genes in Bile secretion pathway, which is related to metabolism in hepatocellular carcinoma. The survival analysis using the database TCGA also show that two cell subpopulations after imputing have distinguished survival rates. [ABSTRACT FROM AUTHOR]

In this paper, we focus on the study of evolutionary algorithms for solving multiobjective optimization problems with a large number of objectives. First, a comparative study of a newly developed dynamical muitiobjective evolutionary algorithm (DMOEA) and some modern algorithms, such as the indicator-based evolutionary algorithm, multiple single objective Pareto sampling, and nondominated sorting genetic algorithm II, is presented by employing the convergence metric and relative hypervolume metric. For three scalable test problems (namely, DTLZ1, DTLZ2, and DTLZ6), which represent some of the most difficult problems studied in the literature, the DMOEA shows good performance in both converging to the true Pareto-optimal front and maintaining a widely distributed set of solutions. Second, a new definition of optimality (namely, L-optimality) is proposed in this paper, which not only takes into account the number of improved objective values but also considers the values of improved objective functions if all objectives have the same importance. We prove that L-optimal solutions are subsets of Pareto-optimal solutions. Finally, the new algorithm based on L-optimality (namely, MDMOEA) is developed, and simulation and comparative results indicate that well-distributed L-optimal solutions can be obtained by utilizing the MDMOEA but cannot be achieved by applying L-optimality to make a posteriori selection within the huge Pareto nondominated solutions. We can conclude that our new algorithm is suitable to tackle many-objective problems. Index Terms--Evolutionary algorithm, many-objective optimization, Pareto optimality, performance assessment.

Patients with coronavirus disease 2019 (COVID-19) often exhibit diverse disease progressions associated with various infectious ability, symptoms, and clinical treatments. To systematically and thoroughly understand the heterogeneous progression of COVID-19, we developed a multi-scale computational model to quantitatively understand the heterogeneous progression of COVID-19 patients infected with severe acute respiratory syndrome (SARS)-like coronavirus (SARS-CoV-2). The model consists of intracellular viral dynamics, multicellular infection process, and immune responses, and was formulated using a combination of differential equations and stochastic modeling. By integrating multi-source clinical data with model analysis, we quantified individual heterogeneity using two indexes, i.e., the ratio of infected cells and incubation period. Specifically, our simulations revealed that increasing the host antiviral state or virus induced type I interferon (IFN) production rate can prolong the incubation period and postpone the transition from asymptomatic to symptomatic outcomes. We further identified the threshold dynamics of T cell exhaustion in the transition between mild-moderate and severe symptoms, and that patients with severe symptoms exhibited a lack of naïve T cells at a late stage. In addition, we quantified the efficacy of treating COVID-19 patients and investigated the effects of various therapeutic strategies. Simulations results suggested that single antiviral therapy is sufficient for moderate patients, while combination therapies and prevention of T cell exhaustion are needed for severe patients. These results highlight the critical roles of IFN and T cell responses in regulating the stage transition during COVID-19 progression. Our study reveals a quantitative relationship underpinning the heterogeneity of transition stage during COVID-19 progression and can provide a potential guidance for personalized therapy in COVID-19 patients. Author summary: Coronavirus disease 2019 (COVID-19) is currently destroying both lives and economies. However, patients infected with severe acute respiratory syndrome (SARS)-like coronavirus (SARS-CoV-2) usually present heterogeneous and complicated progressions, such as different incubation periods (short and long), symptoms (asymptomatic and symptomatic) and severity (mild-moderate and severe). Currently, various clinical data and experimental data are available from different countries, which has great significance for integrating different types of data to comprehensively understand the diverse disease progression in COVID-19 patients and guide individual treatment strategies. Here, we developed a multi-scale computational model to describe the dynamical process of patients infected with SARS-CoV-2, including intracellular viral dynamics, multicellular infection process, and immune responses. By combining data integration, stochastic simulation and quantitative analysis based on the multi-scale mathematical model, we addressed an important question regarding how IFN response and T cell exhaustion quantitatively affect heterogeneous progression in patients with respect to incubation periods, symptoms and severity. Furthermore, the efficacy of various therapeutic strategies for treating COVID-19 patients with different severity degrees was evaluated and validated. The computational framework in this study can also be extended to explore the dynamical process of other coronavirus infections. [ABSTRACT FROM AUTHOR]

Identification of intercellular signaling changes across multiple single-cell RNA-sequencing (scRNA-seq) datasets as well as how intercellular communications affect intracellular transcription factors (TFs) to regulate target genes is crucial in understanding how distinct cell states respond to evolution, perturbations, and diseases. Here, we first generalized our previously developed tool CellChat, enabling flexible comparison analysis of cell–cell communication networks across any number of scRNA-seq datasets from interrelated biological conditions. This greatly facilitates the ready detection of signaling changes of cell–cell communication in response to any biological perturbations. We then investigated how intercellular communications affect intracellular signaling response by inferring a multiscale signaling network which bridges the intercellular communications at the population level and the cell state–specific intracellular signaling network at the molecular level. The latter is constructed by integrating receptor-TF interactions collected from public databases and TF-target gene regulations inferred from a network-regularized regression model. By applying our approaches to three scRNA-seq datasets from skin development, spinal cord injury, and COVID-19, we demonstrated the capability of our approaches in identifying the predominant signaling changes across conditions and the critical signaling mechanisms regulating target gene expression. Together, our work will facilitate the identification of both intercellular and intracellular dysregulated signaling mechanisms responsible for biological perturbations in diverse tissues. [ABSTRACT FROM AUTHOR]

The ongoing coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) has serious influences on human health and economy. The available clinical data suggest that patients infected by SARS-CoV-2 have the possibility of simultaneous infection of bacteria. In this study, we present a data-driven mathematical model for coinfection of SARS-CoV-2 and bacteria to investigate the dynamics of COVID-19 progress. Specifically, based on the statistical analysis of different clinical data from China and some other countries, a system model with ordinary differential equations (ODEs) in four variables, i.e. SARS-CoV-2, bacteria, neutrophils and lymphocytes, is established. We further validate our model through theoretical analysis and fitting of different clinical data. Moreover, through numerical simulations and bifurcation analysis, we find that bacterial infection and immune-related parameters in certain ranges lead to the system transitions among three steady states, i.e. mild, severe and death. We also analyzed the influence of the time it takes for patients to switch from a high-risk area to a low-risk area on the recovery time. These results reveal that the coinfection of viruses and bacteria can explain the changes in neutrophils and lymphocytes, and that initial bacterial infection and immune-related parameters have great influences on the severity degree and mortality in COVID-19 patients. Together, our model and quantitative analysis suggest that preventing bacterial infection and enhancing immune ability during the early phase of infections could be a potential treatment option for high-risk COVID-19 patients. [ABSTRACT FROM AUTHOR]

The Human Cell Atlas (HCA) is a large project that aims to identify all cell types in the human body. The dimension reduction and clustering for identification of cell types from single-cell RNA-sequencing (scRNA-seq) data have become foundational approaches to HCA. The major challenges of current computational analyses are of poor performance on large scale data and sensitive to initial data. We present a new ensemble framework called Adaptive Slice KNNs (scASK) to address the challenges for analyzing scRNA-seq data with high dimensionality. scASK consists of three innovational modules, called DAS (Data Adaptive Slicing), MCS (Meta Classifiers Selecting) and EMS (Ensemble Mode Switching), respectively, which facilitate scASK to approximate a bias-variance tradeoff beyond classification. Thirteen real scRNA-seq datasets are used to evaluate the performance of scASK. Compared with five popular classification algorithms, our experimental results indicate that scASK achieves the best accuracy and robustness among all competing methods. In conclusion, adaptive slicing is an effective structural reduction procedure, and meanwhile scASK provides novel and robust ensemble framework especially for classifying cell types based on scRNA-seq data. scASK is now publically available at https://github.com/liubo2358/scASKcmd. [ABSTRACT FROM AUTHOR]

Alternative splicing makes significant contributions to functional diversity of transcripts and proteins. Many alternatively spliced gene isoforms have been shown to perform specific biological functions under different contexts. In addition to gene-level expression, the advances of high-throughput sequencing offer a chance to estimate isoform-specific exon expression with a high resolution, which is informative for studying splice variants with network analysis. In this study, we propose a novel network-based analysis framework to predict isoform-specific functions from exon-level RNA-Seq data. In particular, based on exon-level expression data, we firstly propose a unified framework, referred to as Iso-Net , to integrate two new mathematical methods (named MINet and RVNet) that infer co-expression networks at different data scenarios. We demonstrate the superior prediction accuracy of Iso-Net over the existing methods for most simulation data, especially in two extreme cases: sample size is very small and exon numbers of two isoforms are quite different. Furthermore, by defining relevant quantitative measures (e.g., Jaccard correlation coefficient) and combining differential co-expression network analysis and GO functional enrichment analysis, a co-expression network analysis framework is developed to predict functions of isoforms and further, to discover their distinct functions within the same gene. We apply Iso-Net to study gene isoforms for several important transcription factors in human myeloid differentiation with the exon-level RNA-Seq data from three different cell lines. Iso-Net is open source and freely available from https://github.com/Dingjie-Wang/Iso-Net. [ABSTRACT FROM AUTHOR]

Single-cell RNA transcriptome data present a tremendous opportunity for studying the cellular heterogeneity. Identifying subpopulations based on scRNA-seq data is a hot topic in recent years, although many researchers have been focused on designing elegant computational methods for identifying new cell types; however, the performance of these methods is still unsatisfactory due to the high dimensionality, sparsity and noise of scRNA-seq data. In this study, we propose a new cell type detection method by learning a robust and accurate similarity matrix, named SCCLRR. The method simultaneously captures both global and local intrinsic properties of data based on a low rank representation (LRR) framework mathematical model. The integrated normalized Euclidean distance and cosine similarity are used to balance the intrinsic linear and nonlinear manifold of data in the local regularization term. To solve the non-convex optimization model, we present an iterative optimization procedure using the alternating direction method of multipliers (ADMM) algorithm. We evaluate the performance of the SCCLRR method on nine real scRNA-seq datasets and compare it with seven state-of-the-art methods. The simulation results show that the SCCLRR outperforms other methods and is robust and effective for clustering scRNA-seq data. (The code of SCCLRR is free available for academic https://github.com/wzhangwhu/SCCLRR). [ABSTRACT FROM AUTHOR]

Essential proteins are indispensable for maintaining normal cellular functions. Identification of essential proteins from Protein-protein interaction (PPI) networks has become a hot topic in recent years. Traditionally biological experimental based approaches are time-consuming and expensive, although lots of computational based methods have been developed in the past years; however, the prediction accuracy is still unsatisfied. In this research, by introducing the protein sub-cellular localization information, we define a new measurement for characterizing the protein's subcellular localization essentiality, and a new data fusion based method is developed for identifying essential proteins, named TEGS, based on integrating network topology, gene expression profile, GO annotation information, and protein subcellular localization information. To demonstrate the efficiency of the proposed method TEGS, we evaluate its performance on two $Saccharomyces cerevisiae$ S a c c h a r o m y c e s c e r e v i s i a e datasets and compare with other seven state-of-the-art methods (DC, BC, NC, PeC, WDC, SON, and TEO) in terms of true predicted number, jackknife curve, and precision-recall curve. Simulation results show that the TEGS outperforms the other compared methods in identifying essential proteins. The source code of TEGS is freely available at https://github.com/wzhangwhu/TEGS. [ABSTRACT FROM AUTHOR]

Disease development and cell differentiation both involve dynamic changes; therefore, the reconstruction of dynamic gene regulatory networks (DGRNs) is an important but difficult problem in systems biology. With recent technical advances in single-cell RNA sequencing (scRNA-seq), large volumes of scRNA-seq data are being obtained for various processes. However, most current methods of inferring DGRNs from bulk samples may not be suitable for scRNA-seq data. In this work, we present scPADGRN, a novel DGRN inference method using "time-series" scRNA-seq data. scPADGRN combines the preconditioned alternating direction method of multipliers with cell clustering for DGRN reconstruction. It exhibits advantages in accuracy, robustness and fast convergence. Moreover, a quantitative index called Differentiation Genes' Interaction Enrichment (DGIE) is presented to quantify the interaction enrichment of genes related to differentiation. From the DGIE scores of relevant subnetworks, we infer that the functions of embryonic stem (ES) cells are most active initially and may gradually fade over time. The communication strength of known contributing genes that facilitate cell differentiation increases from ES cells to terminally differentiated cells. We also identify several genes responsible for the changes in the DGIE scores occurring during cell differentiation based on three real single-cell datasets. Our results demonstrate that single-cell analyses based on network inference coupled with quantitative computations can reveal key transcriptional regulators involved in cell differentiation and disease development. Author summary: Single-cell RNA sequencing (scRNA-seq) data are gaining popularity for providing access to cell-level measurements. Currently, time-series scRNA-seq data allow researchers to study dynamic changes during biological processes. This work proposes a novel method, scPADGRN, for application to time-series scRNA-seq data to construct dynamic gene regulatory networks, which are informative for investigating dynamic changes during disease development and cell differentiation. The proposed method shows satisfactory performance on both simulated data and three real datasets concerning cell differentiation. To quantify network dynamics, we present a quantitative index, DGIE, to measure the degree of activity of a certain set of genes in a regulatory network. Quantitative computations based on dynamic networks identify key regulators in cell differentiation and reveal the activity states of the identified regulators. Specifically, Bhlhe40, Msx2, Foxa2 and Dnmt3l might be important regulatory genes involved in differentiation from mouse ES cells to primitive endoderm (PrE) cells. For differentiation from mouse embryonic fibroblast cells to myocytes, Scx, Fos and Tcf12 are suggested to be key regulators. Sox5, Meis2, Hoxb3, Tcf7l1 and Plagl1 critically contribute during differentiation from human ES cells to definitive endoderm cells. These results may guide further theoretical and experimental efforts to understand cell differentiation processes and explore cell heterogeneity. [ABSTRACT FROM AUTHOR]

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the “fate decisions” of transition cells within the transition. The present study provides new insights into transition events during differentiation and development. [ABSTRACT FROM AUTHOR]

The prediction of protein complexes based on the protein interaction network is a fundamental task for the understanding of cellular life as well as the mechanisms underlying complex disease. A great number of methods have been developed to predict protein complexes based on protein–protein interaction (PPI) networks in recent years. However, because the high throughput data obtained from experimental biotechnology are incomplete, and usually contain a large number of spurious interactions, most of the network-based protein complex identification methods are sensitive to the reliability of the PPI network. In this paper, we propose a new method, Identification of Protein Complex based on Refined Protein Interaction Network (IPC-RPIN), which integrates the topology, gene expression profiles and GO functional annotation information to predict protein complexes from the reconstructed networks. To demonstrate the performance of the IPC-RPIN method, we evaluated the IPC-RPIN on three PPI networks of Saccharomycescerevisiae and compared it with four state-of-the-art methods. The simulation results show that the IPC-RPIN achieved a better result than the other methods on most of the measurements and is able to discover small protein complexes which have traditionally been neglected. [ABSTRACT FROM AUTHOR]

The identification of essential proteins in protein-protein interaction (PPI) networks is of great significance for understanding cellular processes. With the increasing availability of large-scale PPI data, numerous centrality measures based on network topology have been proposed to detect essential proteins from PPI networks. However, most of the current approaches focus mainly on the topological structure of PPI networks, and largely ignore the gene ontology annotation information. In this paper, we propose a novel centrality measure, called TEO, for identifying essential proteins by combining network topology, gene expression profiles, and GO information. To evaluate the performance of the TEO method, we compare it with five other methods (degree, betweenness, NC, Pec, and CowEWC) in detecting essential proteins from two different yeast PPI datasets. The simulation results show that adding GO information can effectively improve the predicted precision and that our method outperforms the others in predicting essential proteins. [ABSTRACT FROM PUBLISHER]

Gag, as the major structural protein of HIV-1, is necessary for the assembly of the HIV-1 sphere shell. An in-depth understanding of its trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. We developed a mathematical model to simulate two biophysical processes, specifically Gag monomer and dimer transport in the cytoplasm and the polymerization of monomers to form a hexamer underneath the plasma membrane. Using experimental data, an optimization approach was utilized to identify the model parameters, and the identifiability and sensitivity of these parameters were then analyzed. Using our model, we analyzed the weight of the pathways involved in the polymerization reactions and concluded that the predominant pathways for the formation of a hexamer might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. We then deduced that the dimer and trimer intermediates might be crucial in hexamer formation. We also explored four theoretical combined methods for Gag suppression, and hypothesized that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. This work serves as a guide for future theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization, and might help accelerate the efficiency of anti-AIDS drug design. [ABSTRACT FROM AUTHOR]

The controllability of multilayer networks has become increasingly important in many areas of science and engineering. In this paper, we identify the general rules that determine the controllability and control energy cost of multilayer networks. First, we quantitatively estimate the control energy cost of multilayer networks and investigate the impacts of different coupling strength and coupling patterns on the control energy cost for multilayer networks. The results indicate that the average energy and the coupling strength have an approximately linear relationship in multilayer networks with two layers. Second, we study how the coupling strength and the connection patterns between different layers affect the controllability of multilayer networks from both theoretical and numerical aspects. The obtained piecewise functional relations between the controllability's measure and coupling strength reveal the existence of an optimal coupling strength for the different interconnection strategies in multilayer networks. In particular, the numerical experiments demonstrate that there exists a tradeoff between the optimal controllability and optimal control energy for selecting interlayer connection patterns in multilayer networks. These results provide a comprehensive understanding of the impact of interlayer couplings on the controllability and control energy cost for multilayer networks and provide a methodology for selecting the control nodes and coupling strength to maximize the controllability and minimize the control energy cost. [ABSTRACT FROM AUTHOR]

With the increasing availability of high-throughput data, various computational methods have recently been developed for understanding the cell through protein-protein interaction (PPI) networks at a systems level. However, due to the incompleteness of the original PPI networks those efforts have been significantly hindered. In this paper, we propose a two stage method to predict underlying links between two originally unlinked protein pairs. First, we measure gene expression and gene functional similarly between unlinked protein pairs on Saccharomyces cerevisiae benchmark network and obtain new constructed networks. Then, we select the significant part of the new predicted links by analyzing the difference between essential proteins that have been identified based on the new constructed networks and the original network. Furthermore, we validate the performance of the new method by using the reliable and comprehensive PPI dataset obtained from the STRING database and compare the new proposed method with four other random walk-based methods. Comparing the results indicates that the new proposed strategy performs well in predicting underlying links. This study provides a general paradigm for predicting new interactions between protein pairs and offers new insights into identifying essential proteins. [ABSTRACT FROM AUTHOR]