VnD: a structure-centric database of disease-related SNPs and drugs

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Role of Structural Bioinformatics in Drug Discovery by Computational SNP Analysis: Analyzing Variation at the Protein Level

With the completion of the human genome project at the beginning of the 21st century, the biological sciences entered an unprecedented age of data generation, and made its first steps toward an era of personalized medicine. This abundance of sequence data has led to the proliferation of numerous sequence-based techniques for associating variation with disease, such as genome-wide association studies and candidate gene association studies. However, these statistical methods do not provide an understanding of the functional effects of variation. Structure-based drug discovery and design is increasingly incorporating structural bioinformatics techniques to model and analyze protein targets, perform large scale virtual screening to identify hit to lead compounds, and simulate molecular interactions. These techniques are fast, cost-effective, and complement existing experimental techniques such as high throughput sequencing. In this paper, we discuss the contributions of structural bioinformatics to drug discovery, focusing particularly on the analysis of nonsynonymous single nucleotide polymorphisms. We conclude by suggesting a protocol for future analyses of the structural effects of nonsynonymous single nucleotide polymorphisms on proteins and protein complexes.

Identifying putative drug targets and potential drug leads: starting points for virtual screening and docking

The availability of 3D models of both drug leads (small molecule ligands) and drug targets (proteins) is essential to molecular docking and computational drug discovery. This chapter describes a simple approach that can be used to identify both drug leads and drug targets using two popular Web-accessible databases: (1) DrugBank and (2) The Human Metabolome Database. First, it is illustrated how putative drug targets and drug leads for exogenous diseases (i.e., infectious diseases) can be readily identified and their 3D structures selected using only the genomic sequences from pathogenic bacteria or viruses as input. The second part illustrates how putative drug targets and drug leads for endogenous diseases (i.e., noninfectious diseases or chronic conditions) can be identified using similar databases and similar sequence input. This chapter is intended to illustrate how bioinformatics and cheminformatics can work synergistically to help provide the necessary inputs for computer-aided drug design.

RNA-Seq approach for genetic improvement of meat quality in pig and evolutionary insight into the substrate specificity of animal carbonyl reductases

Changes in meat quality traits are strongly associated with alterations in postmortem metabolism which depend on genetic variations, especially nonsynonymous single nucleotide variations (nsSNVs) having critical effects on protein structure and function. To selectively identify metabolism-related nsSNVs, next-generation transcriptome sequencing (RNA-Seq) was carried out using RNAs from porcine liver, which contains a diverse range of metabolic enzymes. The multiplex SNV genotyping analysis showed that various metabolism-related genes had different nsSNV alleles. Moreover, many nsSNVs were significantly associated with multiple meat quality traits. Particularly, ch7:g.22112616A>G SNV was identified to create a single amino acid change (Thr/Ala) at the 145th residue of H1.3-like protein, very close to the putative 147th threonine phosphorylation site, suggesting that the nsSNV may affect multiple meat quality traits by affecting the epigenetic regulation of postmortem metabolism-related gene expression. Besides, one nonsynonymous variation, probably generated by gene duplication, led to a stop signal in porcine testicular carbonyl reductase (PTCR), resulting in a C-terminal (E281-A288) deletion. Molecular docking and energy minimization calculations indicated that the binding affinity of wild-type PTCR to 5α-DHT, a C₂₁-steroid, was superior to that of C-terminal-deleted PTCR or human carbonyl reductase, which was very consistent with experimental data, reported previously. Furthermore, P284 was identified as an important residue mediating the specific interaction between PTCR and 5α-DHT, and phylogenetic analysis showed that P284 is an evolutionarily conserved residue among animal carbonyl reductases, which suggests that the C-terminal tails of these reductases may have evolved under evolutionary pressure to increase the substrate specificity for C₂₁-steroids and facilitate metabolic adaptation. Altogether, our RNA-Seq revealed that selective nsSNVs were associated with meat quality traits that could be useful for successful marker-assisted selection in pigs and also represents a useful resource to enhance understanding of protein folding, substrate specificity, and the evolution of enzymes such as carbonyl reductase.

Bioinformatics and variability in drug response: a protein structural perspective

Marketed drugs frequently perform worse in clinical practice than in the clinical trials on which their approval is based. Many therapeutic compounds are ineffective for a large subpopulation of patients to whom they are prescribed; worse, a significant fraction of patients experience adverse effects more severe than anticipated. The unacceptable risk-benefit profile for many drugs mandates a paradigm shift towards personalized medicine. However, prior to adoption of patient-specific approaches, it is useful to understand the molecular details underlying variable drug response among diverse patient populations. Over the past decade, progress in structural genomics led to an explosion of available three-dimensional structures of drug target proteins while efforts in pharmacogenetics offered insights into polymorphisms correlated with differential therapeutic outcomes. Together these advances provide the opportunity to examine how altered protein structures arising from genetic differences affect protein-drug interactions and, ultimately, drug response. In this review, we first summarize structural characteristics of protein targets and common mechanisms of drug interactions. Next, we describe the impact of coding mutations on protein structures and drug response. Finally, we highlight tools for analysing protein structures and protein-drug interactions and discuss their application for understanding altered drug responses associated with protein structural variants.

Protein localization as a principal feature of the etiology and comorbidity of genetic diseases

Proteins localized within the same subcellular compartment tend to be functionally associated. This study shows that subcellular localization and network distance between disease-associated proteins provide complementary information explaining patterns of disease comorbidity.

A positive correlation was found between subcellular localization of disease-associated protein pairs and measures of comorbidity.

A higher comorbidity tendency was found for disease-associated protein pairs that are positioned within a shorter distance in the protein interaction network.

The integration of subcellular localization information with protein interaction network sheds light onto the potential molecular connections underlying comorbidity patterns and will help to understand the mechanisms of human disease.

It was shown that the emergence of phenotypically similar diseases are triggered as a result of molecular connections between disease-causing genes (Oti and Brunner, 2007; Zaghloul and Katsanis, 2010). From a genetics, perspective diseases are associated with certain genes (Goh et al, 2007; Feldman et al, 2008), whereas from a proteomics perspective phenotypically similar diseases are connected via biological modules such as protein–protein interactions (PPIs) or molecular pathways (Lage et al, 2007; Jiang et al, 2008; Wu et al, 2008; Linghu et al, 2009; Suthram et al, 2010). These molecular connections between diseases were observed on the population level as well: diseases connected through molecular connections such as shared genes, PPIs, and metabolic pathways tend to show elevated comorbidity (Rzhetsky et al, 2007; Lee et al, 2008; Zhernakova et al, 2009; Park et al, 2009a, 2009b). While these findings constitute a step toward improving our understanding of the mechanism of disease progression, there are still many more molecule-level connections between disease pairs that need to be explored in order to establish a firmer comorbidity association.

Subcellular localization provides spatial information of proteins in the cell; proteins target subcellular localizations to interact with appropriate partners and form functional complexes in signaling pathways and metabolic processes (Au et al, 2007). Abnormal protein localizations are known to lead to the loss of functional effects in diseases (Luheshi et al, 2008; Laurila and Vihinen, 2009). For example, mis-localizations of nuclear/cytoplasmic transport have been detected in many types of carcinoma cells (Kau et al, 2004). A proper identification of protein subcellular localization can hence be useful in discovering disease-associated proteins (Giallourakis et al, 2005; Calvo and Mootha, 2010). With this understanding, we postulate that disease-associated proteins connected by subcellular localizations could also explain the phenotypic similarities between diseases. Furthermore, such connections may also couple to disease progressions that contribute to multiple disease manifestation, that is, comorbidity.

Protein subcellular localization has been extensively studied through various methods to determine a variety of protein functions. To the best of our knowledge, the connection between diseases and subcellular localizations are yet to be studied systematically. To resolve this we constructed, for the first time, a human Disease-associated Protein and subcellular Localization (DPL) matrix (top panel in Box 1). Our DPL matrix provides the ‘cellular localization map of diseases' that represents the spatial index of diseases in the cell. We found that each disease shows unique characteristics of subcellular localization profile in the DPL matrix. We were interested in determining whether subsets of 1284 human diseases exhibit distinct enrichment profiles across subcellular localizations. We calculated pairwise correlations and performed a hierarchical clustering of the enrichments of the 1284 diseases across 10 different subcellular localizations.

Our DPL matrix revealed that 778 diseases (∼62%, P=1.40 × 10⁻³) are enriched in a single localization and 273 diseases (∼21%, P=3.45 × 10⁻³) are enriched in dual localizations. In the DPL matrix, certain disease-associated proteins are likely to be found in membrane-bounded organelles such as mitochondria, lysosome, and peroxisome, indicating that the mutations of proteins localized to these compartments are connected to the pathophysiological conditions of those organelles. Meanwhile, certain disease-associated proteins in the DPL matrix are enriched in dual localizations, such as extracellular/plasma membrane or endoplasmic reticulum/Golgi. Although these two pairs of subcellular localizations appear to be distinct compartments at first, they are functionally related compartments in close proximity during protein translocation process in the cell, and thus are likely to share interacting protein partners (Gandhi et al, 2006).

Comorbidity represents the co-occurrence of multiple diseases in the same individual (Lee et al, 2008; Hidalgo et al, 2009; Park et al, 2009a). Many comorbid disease pairs have been shown to share common genes in the human disease network. For example, Diabetes and Alzheimer's disease share a risk factor in angiotensin I converting enzyme, and frequently occur together in an individual. In such instances, comorbidity can be partially attributed to the disease connections on the molecular level. To explore the impact of protein subcellular localization on comorbidity, we hypothesized that certain disease pairs could also be connected via subcellular localization by the molecular connections between the disease-associated proteins (bottom panel in Box 1).

We found a positive correlation between subcellular localization similarity and relative risk (Figure 3B, Pearson's correlation coefficient between relative risk and subcellular localization similarity=0.81, P=2.96 × 10⁻⁵). The subcellular localization similarity represents the correlation of subcellular localization profiles between disease pairs. To our surprise, when we compared the relative risk of disease pairs linked via various molecular connections, we found that disease pairs connected by subcellular localization showed a near three-fold higher comorbidity tendency (with link distances equal to 2 or 3) when compared with random pairs (Figure 3E).

We then assessed quantitatively the impact of network distances and subcellular localizations on the comorbidity tendency of disease pairs. We expected the proteins associated with comorbid disease pairs to be located closely in the protein interaction network via fewer links compared with random disease pairs. Indeed, a higher comorbidity tendency was found when two disease-associated proteins were positioned within a shorter distance (gray plots in Figure 3F). Moreover, when subcellular localization information was combined with small network distances, the comorbidity tendency increased dramatically (orange plots in Figure 3F). It suggests that subcellular localization and close network distances, two conceptually distinct molecular connections, contributed synergistically to the comorbidity tendency.

Disease progression is not restricted to the mutation of disease-causing genes, but also affected by molecular connections in ‘disease modules,' resulting in comorbidity (Fraser, 2006; Lee et al, 2008). In this study, for the first time we applied subcellular localization information to elucidate the molecular connections between comorbid diseases. We believe that, based on our finding, our approach helps to define the boundaries of ‘disease modules.' Taken together, integration of diverse molecular connections should improve the molecular level understanding of hitherto unexplained comorbid disease pairs and help us in expanding the scope of our knowledge of the mechanism of human disease progression.

Proteins targeting the same subcellular localization tend to participate in mutual protein–protein interactions (PPIs) and are often functionally associated. Here, we investigated the relationship between disease-associated proteins and their subcellular localizations, based on the assumption that protein pairs associated with phenotypically similar diseases are more likely to be connected via subcellular localization. The spatial constraints from subcellular localization significantly strengthened the disease associations of the proteins connected by subcellular localizations. In particular, certain disease types were more prevalent in specific subcellular localizations. We analyzed the enrichment of disease phenotypes within subcellular localizations, and found that there exists a significant correlation between disease classes and subcellular localizations. Furthermore, we found that two diseases displayed high comorbidity when disease-associated proteins were connected via subcellular localization. We newly explained 7584 disease pairs by using the context of protein subcellular localization, which had not been identified using shared genes or PPIs only. Our result establishes a direct correlation between protein subcellular localization and disease association, and helps to understand the mechanism of human disease progression.