Machine Learning Prediction of Tongue Pressure in Elderly Patients with Head and Neck Tumor: A Cross-Sectional Study.

Background: This investigation sought to cross validate the predictors of tongue pressure recovery in elderly patients' post-treatment for head and neck tumors, leveraging advanced machine learning techniques. Methods: By employing logistic regression, support vector regression, random forest, and extreme gradient boosting, the study analyzed an array of variables including patient demographics, surgery types, dental health status, and age, drawn from comprehensive medical records and direct tongue pressure assessments. Results: Among the models, logistic regression emerged as the most effective, demonstrating an accuracy of 0.630 [95% confidence interval (CI): 0.370-0.778], F1 score of 0.688 [95% confidence interval (CI): 0.435-0.853], precision of 0.611 [95% confidence interval (CI): 0.313-0.801], recall of 0.786 [95% confidence interval (CI): 0.413-0.938] and an area under the receiver operating characteristic curve of 0.626 [95% confidence interval (CI): 0.409-0.806]. This model distinctly highlighted the significance of glossectomy (p = 0.039), the presence of functional teeth (p = 0.043), and the patient's age (p = 0.044) as pivotal factors influencing tongue pressure, setting the threshold for statistical significance at p < 0.05. Conclusions: The analysis underscored the critical role of glossectomy, the presence of functional natural teeth, and age as determinants of tongue pressure in logistics regression, with the presence of natural teeth and the tumor site located in the tongue consistently emerging as the key predictors across all computational models employed in this study.

DTox: A deep neural network-based in visio lens for large scale toxicogenomics data.

With the advancement of large-scale omics technologies, particularly transcriptomics data sets on drug and treatment response repositories available in public domain, toxicogenomics has emerged as a key field in safety pharmacology and chemical risk assessment. Traditional statistics-based bioinformatics analysis poses challenges in its application across multidimensional toxicogenomic data, including administration time, dosage, and gene expression levels. Motivated by the visual inspection workflow of field experts to augment their efficiency of screening significant genes to derive meaningful insights, together with the ability of deep neural architectures to learn the image signals, we developed DTox, a deep neural network-based in visio approach. Using the Percellome toxicogenomics database, instead of utilizing the numerical gene expression values of the transcripts (gene probes of the microarray) for dose-time combinations, DTox learned the image representation of 3D surface plots of distinct time and dosage data points to train the classifier on the experts' labels of gene probe significance. DTox outperformed statistical threshold-based bioinformatics and machine learning approaches based on numerical expression values. This result shows the ability of image-driven neural networks to overcome the limitations of classical numeric value-based approaches. Further, by augmenting the model with explainability modules, our study showed the potential to reveal the visual analysis process of human experts in toxicogenomics through the model weights. While the current work demonstrates the application of the DTox model in toxicogenomic studies, it can be further generalized as an in visio approach for multi-dimensional numeric data with applications in various fields in medical data sciences.

Development of a predictive model for nephrotoxicity during tacrolimus treatment using machine learning methods.

AIM: When administering tacrolimus, therapeutic drug monitoring is recommended because nephrotoxicity, an adverse event, occurs at supra-therapeutic whole-blood concentrations of tacrolimus. However, some patients exhibit nephrotoxicity even at the recommended concentrations, therefore establishing a therapeutic range of tacrolimus concentration for the individual patient is necessary to avoid nephrotoxicity. This study aimed to develop a model for individualized prediction of nephrotoxicity in patients administered tacrolimus. METHODS: We collected data, such as laboratory test data at tacrolimus initiation, concomitant drugs and tacrolimus whole-blood concentration, from medical records of patients who received oral tacrolimus. Nephrotoxicity was defined as an increase in serum creatinine levels within 60 days of tacrolimus initiation. We built 13 prediction models based on different machine learning algorithms: logistic regression, support vector machine, gradient-boosting trees, random forest and neural networks. The best performing model was compared with the conventional model, which classifies patients according to the tacrolimus concentration alone. RESULTS: Data from 163 and 41 patients were used to construct models and evaluate the best performing one, respectively. Most of the patients were diagnosed with inflammatory or autoimmune diseases. The best performing model was built using a support vector machine; it showed a high F2 score of 0.750 and outperformed the conventional model (0.500). CONCLUSIONS: A machine learning model to predict nephrotoxicity in patients during tacrolimus treatment was developed using tacrolimus whole-blood concentration and other patient data. This model could potentially assist in identifying high-risk patients who require individualized target therapeutic concentrations of tacrolimus prior to treatment initiation to prevent nephrotoxicity.

Aberrant morphometric networks in Alzheimer's disease have hemispheric asymmetry and age dependence.

Alzheimer's disease (AD) is associated with abnormal accumulations of hyperphosphorylated tau and amyloid-ß proteins, resulting in unique patterns of atrophy in the brain. We aimed to elucidate some characteristics of the AD's morphometric networks constructed by associating different morphometric features among brain areas and evaluating their relationship to Mini-Mental State Examination total score and age. Three-dimensional T1-weighted (3DT1) image data scanned by the same 1.5T magnetic resonance imaging (MRI) were obtained from 62 AD patients and 41 healthy controls (HCs) and were analysed by using FreeSurfer. The associations of the extracted six morphometric features between regions were estimated by correlation coefficients. The global and local graph theoretical measures for this network were evaluated. Associations between graph theoretical measures and age, sex and cognition were evaluated by multiple regression analysis in each group. Global measures of integration: global efficiency and mean information centrality were significantly higher in AD patients. Local measures of integration: node global efficiency and information centrality were significantly higher in the entorhinal cortex, fusiform gyrus and posterior cingulate cortex of AD patients but only in the left hemisphere. All global measures were correlated with age in AD patients but not in HCs. The information centrality was associated with age in AD's broad brain regions. Our results showed that altered morphometric networks due to AD are left-hemisphere dominant, suggesting that AD pathogenesis has a left-right asymmetry. Ageing has a unique impact on the morphometric networks in AD patients. The information centrality is a sensitive graph theoretical measure to detect this association.

Polysulfide Serves as a Hallmark of Desmoplastic Reaction to Differentially Diagnose Ductal Carcinoma In Situ and Invasive Breast Cancer by SERS Imaging.

Pathological examination of formalin-fixed paraffin-embedded (FFPE) needle-biopsied samples by certified pathologists represents the gold standard for differential diagnosis between ductal carcinoma in situ (DCIS) and invasive breast cancers (IBC), while information of marker metabolites in the samples is lost in the samples. Infrared laser-scanning large-area surface-enhanced Raman spectroscopy (SERS) equipped with gold-nanoparticle-based SERS substrate enables us to visualize metabolites in fresh-frozen needle-biopsied samples with spatial matching between SERS and HE staining images with pathological annotations. DCIS (n = 14) and IBC (n = 32) samples generated many different SERS peaks in finger-print regions of SERS spectra among pathologically annotated lesions including cancer cell nests and the surrounding stroma. The results showed that SERS peaks in IBC stroma exhibit significantly increased polysulfide that coincides with decreased hypotaurine as compared with DCIS, suggesting that alterations of these redox metabolites account for fingerprints of desmoplastic reactions to distinguish IBC from DCIS. Furthermore, the application of supervised machine learning to the stroma-specific multiple SERS signals enables us to support automated differential diagnosis with high accuracy. The results suggest that SERS-derived biochemical fingerprints derived from redox metabolites account for a hallmark of desmoplastic reaction of IBC that is absent in DCIS, and thus, they serve as a useful method for precision diagnosis in breast cancer.

Sub-GOFA: A tool for Sub-Gene Ontology function analysis in clonal mosaicism using semantic (logical) similarity.

Clonal mosaicism (a detectable post-zygotic mutational event in cellular subpopulations) is common in cancer patients. Detected segments of clonal mosaicism are usually bundled into large-locus regions for statistical analysis. However, low-frequency genes are overlooked and are not sufficient to elucidate qualitative differences between cancer patients and non-patients. Therefore, it is of interest to develop and describe a tool named Sub-GOFA for Sub-Gene Ontology function analysis in clonal mosaicism using semantic similarity. Sub-GOFA measures the semantic (logical) similarity among patients using the sub-GO network structures of various sizes segmented from the gene ontology (GO) for clustering analysis. The sub-GO's root-terms with significant differences are extracted as disease-associated genetic functions. Sub-GOFA selected a high ratio of cancer-associated genes under validation with acceptable threshold.

A Prehospital Triage System to Detect Traumatic Intracranial Hemorrhage Using Machine Learning Algorithms.

Importance: An adequate system for triaging patients with head trauma in prehospital settings and choosing optimal medical institutions is essential for improving the prognosis of these patients. To our knowledge, there has been no established way to stratify these patients based on their head trauma severity that can be used by ambulance crews at an injury site. Objectives: To develop a prehospital triage system to stratify patients with head trauma according to trauma severity by using several machine learning techniques and to evaluate the predictive accuracy of these techniques. Design, Setting, and Participants: This single-center retrospective cohort study was conducted by reviewing the electronic medical records of consecutive patients who were transported to Tokyo Medical and Dental University Hospital in Japan from April 1, 2018, to March 31, 2021. Patients younger than 16 years with cardiopulmonary arrest on arrival or with a significant amount of missing data were excluded. Main Outcomes and Measures: Machine learning-based prediction models to detect the presence of traumatic intracranial hemorrhage were constructed. The predictive accuracy of the models was evaluated with the area under the receiver operating curve (ROC-AUC), area under the precision recall curve (PR-AUC), sensitivity, specificity, and other representative statistics. Results: A total of 2123 patients (1527 male patients [71.9%]; mean [SD] age, 57.6 [19.8] years) with head trauma were enrolled in this study. Traumatic intracranial hemorrhage was detected in 258 patients (12.2%). Among several machine learning algorithms, extreme gradient boosting (XGBoost) achieved the mean (SD) highest ROC-AUC (0.78 [0.02]) and PR-AUC (0.46 [0.01]) in cross-validation studies. In the testing set, the ROC-AUC was 0.80, the sensitivity was 74.0% (95% CI, 59.7%-85.4%), and the specificity was 74.9% (95% CI, 70.2%-79.3%). The prediction model using the National Institute for Health and Care Excellence (NICE) guidelines, which was calculated after consultation with physicians, had a sensitivity of 72.0% (95% CI, 57.5%-83.8%) and a specificity of 73.3% (95% CI, 68.7%-77.7%). The McNemar test revealed no statistically significant differences between the XGBoost algorithm and the NICE guidelines for sensitivity or specificity (P = .80 and P = .55, respectively). Conclusions and Relevance: In this cohort study, the prediction model achieved a comparatively accurate performance in detecting traumatic intracranial hemorrhage using only the simple pretransportation information from the patient. Further validation with a prospective multicenter data set is needed.

Toxicity Analysis of Pentachlorophenol Data with a Bioinformatics Tool Set.

Rapid progress in technologies opened the new era of computer-leaded analytics, leaving humans more space for experimental design and decision making. Here we demonstrate the machine learning analysis workflow represented by spectral clustering, elucidation of evolutionary conserved transcription regulation, and network analysis using reverse engineering. Analysis of genes induced by the Pentachlorophenol toxic chemical revealed two subnetworks, one orchestrated by Interferon and another by Nuclear receptor factor 2 (NRF2) gene. Furthermore, network-inference based analysis identified a gene network module composed of genes associated with interferon signaling and their regulatory interaction with downstream genes, especially TRIM family proteins involved in responses of innate immune systems.

Artificial intelligence-based computational framework for drug-target prioritization and inference of novel repositionable drugs for Alzheimer's disease.

BACKGROUND: Identifying novel therapeutic targets is crucial for the successful development of drugs. However, the cost to experimentally identify therapeutic targets is huge and only approximately 400 genes are targets for FDA-approved drugs. As a result, it is inevitable to develop powerful computational tools that can identify potential novel therapeutic targets. Fortunately, the human protein-protein interaction network (PIN) could be a useful resource to achieve this objective. METHODS: In this study, we developed a deep learning-based computational framework that extracts low-dimensional representations of high-dimensional PIN data. Our computational framework uses latent features and state-of-the-art machine learning techniques to infer potential drug target genes. RESULTS: We applied our computational framework to prioritize novel putative target genes for Alzheimer's disease and successfully identified key genes that may serve as novel therapeutic targets (e.g., DLG4, EGFR, RAC1, SYK, PTK2B, SOCS1). Furthermore, based on these putative targets, we could infer repositionable candidate-compounds for the disease (e.g., tamoxifen, bosutinib, and dasatinib). CONCLUSIONS: Our deep learning-based computational framework could be a powerful tool to efficiently prioritize new therapeutic targets and enhance the drug repositioning strategy.

Disease-related cellular protein networks differentially affected under different EGFR mutations in lung adenocarcinoma.

It is unclear how epidermal growth factor receptor EGFR major driver mutations (L858R or Ex19del) affect downstream molecular networks and pathways. This study aimed to provide information on the influences of these mutations. The study assessed 36 protein expression profiles of lung adenocarcinoma (Ex19del, nine; L858R, nine; no Ex19del/L858R, 18). Weighted gene co-expression network analysis together with analysis of variance-based screening identified 13 co-expressed modules and their eigen proteins. Pathway enrichment analysis for the Ex19del mutation demonstrated involvement of SUMOylation, epithelial and mesenchymal transition, ERK/mitogen-activated protein kinase signalling via phosphorylation and Hippo signalling. Additionally, analysis for the L858R mutation identified various pathways related to cancer cell survival and death. With regard to the Ex19del mutation, ROCK, RPS6KA1, ARF1, IL2RA and several ErbB pathways were upregulated, whereas AURK and GSKIP were downregulated. With regard to the L858R mutation, RB1, TSC22D3 and DOCK1 were downregulated, whereas various networks, including VEGFA, were moderately upregulated. In all mutation types, CD80/CD86 (B7), MHC, CIITA and IFGN were activated, whereas CD37 and SAFB were inhibited. Costimulatory immune-checkpoint pathways by B7/CD28 were mainly activated, whereas those by PD-1/PD-L1 were inhibited. Our findings may help identify potential therapeutic targets and develop therapeutic strategies to improve patient outcomes.

A dual controllability analysis of influenza virus-host protein-protein interaction networks for antiviral drug target discovery.

BACKGROUND: Host factors of influenza virus replication are often found in key topological positions within protein-protein interaction networks. This work explores how protein states can be manipulated through controllability analysis: the determination of the minimum manipulation needed to drive the cell system to any desired state. Here, we complete a two-part controllability analysis of two protein networks: a host network representing the healthy cell state and an influenza A virus-host network representing the infected cell state. In this context, controllability analyses aim to identify key regulating host factors of the infected cell's progression. This knowledge can be utilized in further biological analysis to understand disease dynamics and isolate proteins for study as drug target candidates. RESULTS: Both topological and controllability analyses provide evidence of wide-reaching network effects stemming from the addition of viral-host protein interactions. Virus interacting and driver host proteins are significant both topologically and in controllability, therefore playing important roles in cell behavior during infection. Functional analysis finds overlap of results with previous siRNA studies of host factors involved in influenza replication, NF-kB pathway and infection relevance, and roles as interferon regulating genes. 24 proteins are identified as holding regulatory roles specific to the infected cell by measures of topology, controllability, and functional role. These proteins are recommended for further study as potential antiviral drug targets. CONCLUSIONS: Seasonal outbreaks of influenza A virus are a major cause of illness and death around the world each year with a constant threat of pandemic infection. This research aims to increase the efficiency of antiviral drug target discovery using existing protein-protein interaction data and network analysis methods. These results are beneficial to future studies of influenza virus, both experimental and computational, and provide evidence that the combination of topology and controllability analyses may be valuable for future efforts in drug target discovery.

Precise Temporal Profiling of Signaling Complexes in Primary Cells Using SWATH Mass Spectrometry.

Spatiotemporal organization of protein interactions in cell signaling is a fundamental process that drives cellular functions. Given differential protein expression across tissues and developmental stages, the architecture and dynamics of signaling interaction proteomes is, likely, highly context dependent. However, current interaction information has been almost exclusively obtained from transformed cells. In this study, we applied an advanced and robust workflow combining mouse genetics and affinity purification (AP)-SWATH mass spectrometry to profile the dynamics of 53 high-confidence protein interactions in primary T cells, using the scaffold protein GRB2 as a model. The workflow also provided a sufficient level of robustness to pinpoint differential interaction dynamics between two similar, but functionally distinct, primary T cell populations. Altogether, we demonstrated that precise and reproducible quantitative measurements of protein interaction dynamics can be achieved in primary cells isolated from mammalian tissues, allowing resolution of the tissue-specific context of cell-signaling events.

A machine learning approach for the identification of key markers involved in brain development from single-cell transcriptomic data.

BACKGROUND: The ability to sequence the transcriptomes of single cells using single-cell RNA-seq sequencing technologies presents a shift in the scientific paradigm where scientists, now, are able to concurrently investigate the complex biology of a heterogeneous population of cells, one at a time. However, till date, there has not been a suitable computational methodology for the analysis of such intricate deluge of data, in particular techniques which will aid the identification of the unique transcriptomic profiles difference between the different cellular subtypes. In this paper, we describe the novel methodology for the analysis of single-cell RNA-seq data, obtained from neocortical cells and neural progenitor cells, using machine learning algorithms (Support Vector machine (SVM) and Random Forest (RF)). RESULTS: Thirty-eight key transcripts were identified, using the SVM-based recursive feature elimination (SVM-RFE) method of feature selection, to best differentiate developing neocortical cells from neural progenitor cells in the SVM and RF classifiers built. Also, these genes possessed a higher discriminative power (enhanced prediction accuracy) as compared commonly used statistical techniques or geneset-based approaches. Further downstream network reconstruction analysis was carried out to unravel hidden general regulatory networks where novel interactions could be further validated in web-lab experimentation and be useful candidates to be targeted for the treatment of neuronal developmental diseases. CONCLUSION: This novel approach reported for is able to identify transcripts, with reported neuronal involvement, which optimally differentiate neocortical cells and neural progenitor cells. It is believed to be extensible and applicable to other single-cell RNA-seq expression profiles like that of the study of the cancer progression and treatment within a highly heterogeneous tumour.

Network analyses based on comprehensive molecular interaction maps reveal robust control structures in yeast stress response pathways.

Cellular stress responses require exquisite coordination between intracellular signaling molecules to integrate multiple stimuli and actuate specific cellular behaviors. Deciphering the web of complex interactions underlying stress responses is a key challenge in understanding robust biological systems and has the potential to lead to the discovery of targeted therapeutics for diseases triggered by dysregulation of stress response pathways. We constructed large-scale molecular interaction maps of six major stress response pathways in Saccharomyces cerevisiae (baker's or budding yeast). Biological findings from over 900 publications were converted into standardized graphical formats and integrated into a common framework. The maps are posted at http://www.yeast-maps.org/yeast-stress-response/ for browse and curation by the research community. On the basis of these maps, we undertook systematic analyses to unravel the underlying architecture of the networks. A series of network analyses revealed that yeast stress response pathways are organized in bow-tie structures, which have been proposed as universal sub-systems for robust biological regulation. Furthermore, we demonstrated a potential role for complexes in stabilizing the conserved core molecules of bow-tie structures. Specifically, complex-mediated reversible reactions, identified by network motif analyses, appeared to have an important role in buffering the concentration and activity of these core molecules. We propose complex-mediated reactions as a key mechanism mediating robust regulation of the yeast stress response. Thus, our comprehensive molecular interaction maps provide not only an integrated knowledge base, but also a platform for systematic network analyses to elucidate the underlying architecture in complex biological systems.

Harnessing diversity towards the reconstructing of large scale gene regulatory networks.

Elucidating gene regulatory network (GRN) from large scale experimental data remains a central challenge in systems biology. Recently, numerous techniques, particularly consensus driven approaches combining different algorithms, have become a potentially promising strategy to infer accurate GRNs. Here, we develop a novel consensus inference algorithm, TopkNet that can integrate multiple algorithms to infer GRNs. Comprehensive performance benchmarking on a cloud computing framework demonstrated that (i) a simple strategy to combine many algorithms does not always lead to performance improvement compared to the cost of consensus and (ii) TopkNet integrating only high-performance algorithms provide significant performance improvement compared to the best individual algorithms and community prediction. These results suggest that a priori determination of high-performance algorithms is a key to reconstruct an unknown regulatory network. Similarity among gene-expression datasets can be useful to determine potential optimal algorithms for reconstruction of unknown regulatory networks, i.e., if expression-data associated with known regulatory network is similar to that with unknown regulatory network, optimal algorithms determined for the known regulatory network can be repurposed to infer the unknown regulatory network. Based on this observation, we developed a quantitative measure of similarity among gene-expression datasets and demonstrated that, if similarity between the two expression datasets is high, TopkNet integrating algorithms that are optimal for known dataset perform well on the unknown dataset. The consensus framework, TopkNet, together with the similarity measure proposed in this study provides a powerful strategy towards harnessing the wisdom of the crowds in reconstruction of unknown regulatory networks.

A comprehensive map of the influenza A virus replication cycle.

BACKGROUND: Influenza is a common infectious disease caused by influenza viruses. Annual epidemics cause severe illnesses, deaths, and economic loss around the world. To better defend against influenza viral infection, it is essential to understand its mechanisms and associated host responses. Many studies have been conducted to elucidate these mechanisms, however, the overall picture remains incompletely understood. A systematic understanding of influenza viral infection in host cells is needed to facilitate the identification of influential host response mechanisms and potential drug targets. DESCRIPTION: We constructed a comprehensive map of the influenza A virus ('IAV') life cycle ('FluMap') by undertaking a literature-based, manual curation approach. Based on information obtained from publicly available pathway databases, updated with literature-based information and input from expert virologists and immunologists, FluMap is currently composed of 960 factors (i.e., proteins, mRNAs etc.) and 456 reactions, and is annotated with ~500 papers and curation comments. In addition to detailing the type of molecular interactions, isolate/strain specific data are also available. The FluMap was built with the pathway editor CellDesigner in standard SBML (Systems Biology Markup Language) format and visualized as an SBGN (Systems Biology Graphical Notation) diagram. It is also available as a web service (online map) based on the iPathways+ system to enable community discussion by influenza researchers. We also demonstrate computational network analyses to identify targets using the FluMap. CONCLUSION: The FluMap is a comprehensive pathway map that can serve as a graphically presented knowledge-base and as a platform to analyze functional interactions between IAV and host factors. Publicly available webtools will allow continuous updating to ensure the most reliable representation of the host-virus interaction network. The FluMap is available at http://www.influenza-x.org/flumap/.

Origin and evolution of genes related to ABA metabolism and its signaling pathways.

Since plants cannot move to avoid stress, they have sophisticated acclimation mechanisms against a variety of abiotic stresses. The phytohormone abscisic acid (ABA) plays essential roles in abiotic stress tolerances in land plants. Therefore, it is interesting to address the evolutionary origins of ABA metabolism and its signaling pathways in land plants. Here, we focused on 48 ABA-related Arabidopsis thaliana genes with 11 protein functions, and generated 11 orthologous clusters of ABA-related genes from A. thaliana, Arabidopsis lyrata, Populus trichocarpa, Oryza sativa, Selaginella moellendorffii, and Physcomitrella patens. Phylogenetic analyses suggested that the common ancestor of these six species possessed most of the key protein functions of ABA-related genes. In two species (A. thaliana and O. sativa), duplicate genes related to ABA signaling pathways contribute to the expression variation in different organs or stress responses. In particular, there is significant expansion of gene families related to ABA in evolutionary periods associated with morphological divergence. Taken together, these results suggest that expansion of the gene families related to ABA signaling pathways may have contributed to the sophisticated stress tolerance mechanisms of higher land plants.

Difference in gene duplicability may explain the difference in overall structure of protein-protein interaction networks among eukaryotes.

BACKGROUND: A protein-protein interaction network (PIN) was suggested to be a disassortative network, in which interactions between high- and low-degree nodes are favored while hub-hub interactions are suppressed. It was postulated that a disassortative structure minimizes unfavorable cross-talks between different hub-centric functional modules and was positively selected in evolution. However, by re-examining yeast PIN data, several researchers reported that the disassortative structure observed in a PIN might be an experimental artifact. Therefore, the existence of a disassortative structure and its possible evolutionary mechanism remains unclear. RESULTS: In this study, we investigated PINs from the yeast, worm, fly, human, and malaria parasite including four different yeast PIN datasets. The analyses showed that the yeast, worm, fly, and human PINs are disassortative while the malaria parasite PIN is not. By conducting simulation studies on the basis of a duplication-divergence model, we demonstrated that a preferential duplication of low- and high-degree nodes can generate disassortative and non-disassortative networks, respectively. From this observation, we hypothesized that the difference in degree dependence on gene duplications accounts for the difference in assortativity of PINs among species. Comparison of 55 proteomes in eukaryotes revealed that genes with lower degrees showed higher gene duplicabilities in the yeast, worm, and fly, while high-degree genes tend to have high duplicabilities in the malaria parasite, supporting the above hypothesis. CONCLUSIONS: These results suggest that disassortative structures observed in PINs are merely a byproduct of preferential duplications of low-degree genes, which might be caused by an organism's living environment.

Structure of protein interaction networks and their implications on drug design.

Protein-protein interaction networks (PINs) are rich sources of information that enable the network properties of biological systems to be understood. A study of the topological and statistical properties of budding yeast and human PINs revealed that they are scale-rich and configured as highly optimized tolerance (HOT) networks that are similar to the router-level topology of the Internet. This is different from claims that such networks are scale-free and configured through simple preferential-attachment processes. Further analysis revealed that there are extensive interconnections among middle-degree nodes that form the backbone of the networks. Degree distributions of essential genes, synthetic lethal genes, synthetic sick genes, and human drug-target genes indicate that there are advantageous drug targets among nodes with middle- to low-degree nodes. Such network properties provide the rationale for combinatorial drugs that target less prominent nodes to increase synergetic efficacy and create fewer side effects.

Non-uniform survival rate of heterodimerization links in the evolution of the yeast protein-protein interaction network.

Protein-protein interaction networks (PINs) are scale-free networks with a small-world property. In a small-world network, the average cluster coefficient () is much higher than in a random network, but the average shortest path length () is similar between the two networks. To understand the evolutionary mechanisms shaping the structure of PINs, simulation studies using various network growth models have been performed. It has been reported that the heterodimerization (HD) model, in which a new link is added between duplicated nodes with a uniform probability, could reproduce scale-freeness and a high . In this paper, however, we show that the HD model is unsatisfactory, because (i) to reproduce the high in the yeast PIN, a much larger number (n(HI)) of HD links (links between duplicated nodes) are required than the estimated number of n(HI) in the yeast PIN and (ii) the spatial distribution of triangles in the yeast PIN is highly skewed but the HD model cannot reproduce the skewed distribution. To resolve these discrepancies, we here propose a new model named the non-uniform heterodimerization (NHD) model. In this model, an HD link is preferentially attached between duplicated nodes when they share many common neighbors. Simulation studies demonstrated that the NHD model can successfully reproduce the high , the low n(HI), and the skewed distribution of triangles in the yeast PIN. These results suggest that the survival rate of HD links is not uniform in the evolution of PINs, and that an HD link between high-degree nodes tends to be evolutionarily conservative. The non-uniform survival rate of HD links can be explained by assuming a low mutation rate for a high-degree node, and thus this model appears to be biologically plausible.