The drug target genes show higher evolutionary conservation than non-target genes

Comparative proteomic profiling of the ovine and human PBMC inflammatory response

Understanding the cellular and molecular mechanisms of inflammation requires robust animal models. Sheep are commonly used in immune-related studies, yet the validity of sheep as animal models for immune and inflammatory diseases remains to be established. This cross-species comparative study analyzed the in vitro inflammatory response of ovine (oPBMCs) and human PBMCs (hPBMCs) using mass spectrometry, profiling the proteome of the secretome and whole cell lysate. Of the entire cell lysate proteome (oPBMCs: 4217, hPBMCs: 4574 proteins) 47.8% and in the secretome proteome (oPBMCs: 1913, hPBMCs: 1375 proteins) 32.8% were orthologous between species, among them 32 orthologous CD antigens, indicating the presence of six immune cell subsets. Following inflammatory stimulation, 71 proteins in oPBMCs and 176 in hPBMCs showed differential abundance, with only 7 overlapping. Network and Gene Ontology analyses identified 16 shared inflammatory-related terms and 17 canonical pathways with similar activation/inhibition patterns in both species, demonstrating significant conservation in specific immune and inflammatory responses. However, ovine PMBCs also contained a unique WC1+γδ T-cell subset, not detected in hPBMCs. Furthermore, differences in the activation/inhibition trends of seven canonical pathways and the sets of DAPs between sheep and humans, emphasize the need to consider interspecies differences in translational studies and inflammation research.

In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model

Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.

CTRP3 exacerbates tendinopathy by dysregulating tendon stem cell differentiation and altering extracellular matrix composition

Tendinopathy, the most common disorder affecting tendons, is characterized by chronic disorganization of the tendon matrix, which leads to tendon tear and rupture. The goal was to identify a rational molecular target whose blockade can serve as a potential therapeutic intervention for tendinopathy. We identified C1q/TNF-related protein-3 (CTRP3) as a markedly up-regulated cytokine in human and rodent tendinopathy. Overexpression of CTRP3 enhanced the progression of tendinopathy by accumulating cartilaginous proteoglycans and degenerating collagenous fibers in the mouse tendon, whereas CTRP3 knockdown suppressed the tendinopathy pathogenesis. Functional blockade of CTRP3 using a neutralizing antibody ameliorated overuse-induced tendinopathy of the Achilles and rotator cuff tendons. Mechanistically, CTRP3 elicited a transcriptomic pattern that stimulates abnormal differentiation of tendon stem/progenitor cells and ectopic chondrification as an effect linked to activation of Akt signaling. Collectively, we reveal an essential role for CTRP3 in tendinopathy and propose a potential therapeutic strategy for the treatment of tendinopathy.

Evaluating drug targets through human loss-of-function genetic variation

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

Sequence-Derived Markers of Drug Targets and Potentially Druggable Human Proteins

Recent research shows that majority of the druggable human proteome is yet to be annotated and explored. Accurate identification of these unexplored druggable proteins would facilitate development, screening, repurposing, and repositioning of drugs, as well as prediction of new drug–protein interactions. We contrast the current drug targets against the datasets of non-druggable and possibly druggable proteins to formulate markers that could be used to identify druggable proteins. We focus on the markers that can be extracted from protein sequences or names/identifiers to ensure that they can be applied across the entire human proteome. These markers quantify key features covered in the past works (topological features of PPIs, cellular functions, and subcellular locations) and several novel factors (intrinsic disorder, residue-level conservation, alternative splicing isoforms, domains, and sequence-derived solvent accessibility). We find that the possibly druggable proteins have significantly higher abundance of alternative splicing isoforms, relatively large number of domains, higher degree of centrality in the protein-protein interaction networks, and lower numbers of conserved and surface residues, when compared with the non-druggable proteins. We show that the current drug targets and possibly druggable proteins share involvement in the catalytic and signaling functions. However, unlike the drug targets, the possibly druggable proteins participate in the metabolic and biosynthesis processes, are enriched in the intrinsic disorder, interact with proteins and nucleic acids, and are localized across the cell. To sum up, we formulate several markers that can help with finding novel druggable human proteins and provide interesting insights into the cellular functions and subcellular locations of the current drug targets and potentially druggable proteins.

Tissue-specific genes as an underutilized resource in drug discovery

Tissue-specific genes are believed to be good drug targets due to improved safety. Here we show that this intuitive notion is not reflected in phase 1 and 2 clinical trials, despite the historic success of tissue-specific targets and their 2.3-fold overrepresentation among targets of marketed non-oncology drugs. We compare properties of tissue-specific genes and drug targets. We show that tissue-specificity of the target may also be related to efficacy of the drug. The relationship may be indirect (enrichment in Mendelian disease and PTVesc genes) or direct (elevated betweenness centrality scores for tissue-specifically produced enzymes and secreted proteins). Reduced evolutionary conservation of tissue-specific genes may represent a bottleneck for drug projects, prompting development of novel models with smaller evolutionary gap to humans. We show that the opportunities to identify tissue-specific drug targets are not exhausted and discuss potential use cases for tissue-specific genes in drug research.

Metabolic network-based stratification of hepatocellular carcinoma reveals three distinct tumor subtypes

Hepatocellular carcinoma (HCC) is one of the most frequent forms of liver cancer, and effective treatment methods are limited due to tumor heterogeneity. There is a great need for comprehensive approaches to stratify HCC patients, gain biological insights into subtypes, and ultimately identify effective therapeutic targets. We stratified HCC patients and characterized each subtype using transcriptomics data, genome-scale metabolic networks and network topology/controllability analysis. This comprehensive systems-level analysis identified three distinct subtypes with substantial differences in metabolic and signaling pathways reflecting at genomic, transcriptomic, and proteomic levels. These subtypes showed large differences in clinical survival associated with altered kynurenine metabolism, WNT/β-catenin-associated lipid metabolism, and PI3K/AKT/mTOR signaling. Integrative analyses indicated that the three subtypes rely on alternative enzymes (e.g., ACSS1/ACSS2/ACSS3, PKM/PKLR, ALDOB/ALDOA, MTHFD1L/MTHFD2/MTHFD1) to catalyze the same reactions. Based on systems-level analysis, we identified 8 to 28 subtype-specific genes with pivotal roles in controlling the metabolic network and predicted that these genes may be targeted for development of treatment strategies for HCC subtypes by performing in silico analysis. To validate our predictions, we performed experiments using HepG2 cells under normoxic and hypoxic conditions and observed opposite expression patterns between genes expressed in high/moderate/low-survival tumor groups in response to hypoxia, reflecting activated hypoxic behavior in patients with poor survival. In conclusion, our analyses showed that the heterogeneous HCC tumors can be stratified using a metabolic network-driven approach, which may also be applied to other cancer types, and this stratification may have clinical implications to drive the development of precision medicine.

Evolutionary Perspectives of Genotype-Phenotype Factors in Leishmania Metabolism

The sandfly midgut and the human macrophage phagolysosome provide antagonistic metabolic niches for the endoparasite Leishmania to survive and populate. Although these environments fluctuate across developmental stages, the relative changes in both these environments across parasite generations might remain gradual. Such environmental restrictions might endow parasite metabolism with a choice of specific genotypic and phenotypic factors that can constrain enzyme evolution for successful adaptation to the host. With respect to the available cellular information for Leishmania species, for the first time, we measure the relative contribution of eight inter-correlated predictors related to codon usage, GC content, gene expression, gene length, multi-functionality, and flux-coupling potential of an enzyme on the evolutionary rates of singleton metabolic genes and further compare their effects across three Leishmania species. Our analysis reveals that codon adaptation, multi-functionality, and flux-coupling potential of an enzyme are independent contributors of enzyme evolutionary rates, which can together explain a large variation in enzyme evolutionary rates across species. We also hypothesize that a species-specific occurrence of duplicated genes in novel subcellular locations can create new flux routes through certain singleton flux-coupled enzymes, thereby constraining their evolution. A cross-species comparison revealed both common and species-specific genes whose evolutionary divergence was constrained by multiple independent factors. Out of these, previously known pharmacological targets and virulence factors in Leishmania were identified, suggesting their evolutionary reasons for being important survival factors to the parasite. All these results provide a fundamental understanding of the factors underlying adaptive strategies of the parasite, which can be further targeted.

Metabolic Network-Based Identification and Prioritization of Anticancer Targets Based on Expression Data in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a deadly form of liver cancer with high mortality worldwide. Unfortunately, the large heterogeneity of this disease makes it difficult to develop effective treatment strategies. Cellular network analyses have been employed to study heterogeneity in cancer, and to identify potential therapeutic targets. However, the existing approaches do not consider metabolic growth requirements, i.e., biological network functionality, to rank candidate targets while preventing toxicity to non-cancerous tissues. Here, we developed an algorithm to overcome these issues based on integration of gene expression data, genome-scale metabolic models, network controllability, and dispensability, as well as toxicity analysis. This method thus predicts and ranks potential anticancer non-toxic controlling metabolite and gene targets. Our algorithm encompasses both objective-driven and-independent tasks, and uses network topology to finally rank the predicted therapeutic targets. We employed this algorithm to the analysis of transcriptomic data for 50 HCC patients with both cancerous and non-cancerous samples. We identified several potential targets that would prevent cell growth, including 74 anticancer metabolites, and 3 gene targets (PRKACA, PGS1, and CRLS1). The predicted anticancer metabolites showed good agreement with existing FDA-approved cancer drugs, and the 3 genes were experimentally validated by performing experiments in HepG2 and Hep3B liver cancer cell lines. Our observations indicate that our novel approach successfully identifies therapeutic targets for effective treatment of cancer. This approach may also be applied to any cancer type that has tumor and non-tumor gene or protein expression data.

Construction and analysis of gene-gene dynamics influence networks based on a Boolean model

BACKGROUND

Identification of novel gene-gene relations is a crucial issue to understand system-level biological phenomena. To this end, many methods based on a correlation analysis of gene expressions or structural analysis of molecular interaction networks have been proposed. They have a limitation in identifying more complicated gene-gene dynamical relations, though.

RESULTS

To overcome this limitation, we proposed a measure to quantify a gene-gene dynamical influence (GDI) using a Boolean network model and constructed a GDI network to indicate existence of a dynamical influence for every ordered pair of genes. It represents how much a state trajectory of a target gene is changed by a knockout mutation subject to a source gene in a gene-gene molecular interaction (GMI) network. Through a topological comparison between GDI and GMI networks, we observed that the former network is denser than the latter network, which implies that there exist many gene pairs of dynamically influencing but molecularly non-interacting relations. In addition, a larger number of hub genes were generated in the GDI network. On the other hand, there was a correlation between these networks such that the degree value of a node was positively correlated to each other. We further investigated the relationships of the GDI value with structural properties and found that there are negative and positive correlations with the length of a shortest path and the number of paths, respectively. In addition, a GDI network could predict a set of genes whose steady-state expression is affected in E. coli gene-knockout experiments. More interestingly, we found that the drug-targets with side-effects have a larger number of outgoing links than the other genes in the GDI network, which implies that they are more likely to influence the dynamics of other genes. Finally, we found biological evidences showing that the gene pairs which are not molecularly interacting but dynamically influential can be considered for novel gene-gene relationships.

CONCLUSION

Taken together, construction and analysis of the GDI network can be a useful approach to identify novel gene-gene relationships in terms of the dynamical influence.