Scoring functions for de novo protein structure prediction revisited

Exploring Polypharmacology in Drug Discovery and Repurposing Using the CANDO Platform

BACKGROUND

Traditional drug discovery approaches focus on a limited set of target molecules for treatment against specific indications/diseases. However, drug absorption, dispersion, metabolism, and excretion (ADME) involve interactions with multiple protein systems. Drugs approved for particular indication(s) may be repurposed as novel therapeutics for others. The severely declining rate of discovery and increasing costs of new drugs illustrate the limitations of the traditional reductionist paradigm in drug discovery.

METHODS

We developed the Computational Analysis of Novel Drug Opportunities (CANDO) platform based on a hypothesis that drugs function by interacting with multiple protein targets to create a molecular interaction signature that can be exploited for therapeutic repurposing and discovery. We compiled a library of compounds that are human ingestible with minimal side effects, followed by an 'all-compounds' vs 'all-proteins' fragment-based multitarget docking with dynamics screen to construct compound-proteome interaction matrices that were then analyzed to determine similarity of drug behavior. The proteomic signature similarity of drugs is then ranked to make putative drug predictions for all indications in a shotgun manner.

RESULTS

We have previously applied this platform with success in both retrospective benchmarking and prospective validation, and to understand the effect of druggable protein classes on repurposing accuracy. Here we use the CANDO platform to analyze and determine the contribution of multitargeting (polypharmacology) to drug repurposing benchmarking accuracy. Taken together with the previous work, our results indicate that a large number of protein structures with diverse fold space and a specific polypharmacological interactome is necessary for accurate drug predictions using our proteomic and evolutionary drug discovery and repurposing platform.

CONCLUSION

These results have implications for future drug development and repurposing in the context of polypharmacology.

Identification of Peptide Inhibitors of Enveloped Viruses Using Support Vector Machine

The peptides derived from envelope proteins have been shown to inhibit the protein-protein interactions in the virus membrane fusion process and thus have a great potential to be developed into effective antiviral therapies. There are three types of envelope proteins each exhibiting distinct structure folds. Although the exact fusion mechanism remains elusive, it was suggested that the three classes of viral fusion proteins share a similar mechanism of membrane fusion. The common mechanism of action makes it possible to correlate the properties of self-derived peptide inhibitors with their activities. Here we developed a support vector machine model using sequence-based statistical scores of self-derived peptide inhibitors as input features to correlate with their activities. The model displayed 92% prediction accuracy with the Matthew's correlation coefficient of 0.84, obviously superior to those using physicochemical properties and amino acid decomposition as input. The predictive support vector machine model for self- derived peptides of envelope proteins would be useful in development of antiviral peptide inhibitors targeting the virus fusion process.

3dRNAscore: a distance and torsion angle dependent evaluation function of 3D RNA structures

Model evaluation is a necessary step for better prediction and design of 3D RNA structures. For proteins, this has been widely studied and the knowledge-based statistical potential has been proved to be one of effective ways to solve this problem. Currently, a few knowledge-based statistical potentials have also been proposed to evaluate predicted models of RNA tertiary structures. The benchmark tests showed that they can identify the native structures effectively but further improvements are needed to identify near-native structures and those with non-canonical base pairs. Here, we present a novel knowledge-based potential, 3dRNAscore, which combines distance-dependent and dihedral-dependent energies. The benchmarks on different testing datasets all show that 3dRNAscore are more efficient than existing evaluation methods in recognizing native state from a pool of near-native states of RNAs as well as in ranking near-native states of RNA models.

A human homeotic transformation resulting from mutations in PLCB4 and GNAI3 causes auriculocondylar syndrome

Auriculocondylar syndrome (ACS) is a rare, autosomal-dominant craniofacial malformation syndrome characterized by variable micrognathia, temporomandibular joint ankylosis, cleft palate, and a characteristic "question-mark" ear malformation. Careful phenotypic characterization of severely affected probands in our cohort suggested the presence of a mandibular patterning defect resulting in a maxillary phenotype (i.e., homeotic transformation). We used exome sequencing of five probands and identified two novel (exclusive to the patient and/or family studied) missense mutations in PLCB4 and a shared mutation in GNAI3 in two unrelated probands. In confirmatory studies, three additional novel PLCB4 mutations were found in multigenerational ACS pedigrees. All mutations were confirmed by Sanger sequencing, were not present in more than 10,000 control chromosomes, and resulted in amino-acid substitutions located in highly conserved protein domains. Additionally, protein-structure modeling demonstrated that all ACS substitutions disrupt the catalytic sites of PLCB4 and GNAI3. We suggest that PLCB4 and GNAI3 are core signaling molecules of the endothelin-1-distal-less homeobox 5 and 6 (EDN1-DLX5/DLX6) pathway. Functional studies demonstrated a significant reduction in downstream DLX5 and DLX6 expression in ACS cases in assays using cultured osteoblasts from probands and controls. These results support the role of the previously implicated EDN1-DLX5/6 pathway in regulating mandibular specification in other species, which, when disrupted, results in a maxillary phenotype. This work defines the molecular basis of ACS as a homeotic transformation (mandible to maxilla) in humans.

GPU-Q-J, a fast method for calculating root mean square deviation (RMSD) after optimal superposition

BACKGROUND

Calculation of the root mean square deviation (RMSD) between the atomic coordinates of two optimally superposed structures is a basic component of structural comparison techniques. We describe a quaternion based method, GPU-Q-J, that is stable with single precision calculations and suitable for graphics processor units (GPUs). The application was implemented on an ATI 4770 graphics card in C/C++ and Brook+ in Linux where it was 260 to 760 times faster than existing unoptimized CPU methods. Source code is available from the Compbio website http://software.compbio.washington.edu/misc/downloads/st_gpu_fit/ or from the author LHH.

FINDINGS

The Nutritious Rice for the World Project (NRW) on World Community Grid predicted de novo, the structures of over 62,000 small proteins and protein domains returning a total of 10 billion candidate structures. Clustering ensembles of structures on this scale requires calculation of large similarity matrices consisting of RMSDs between each pair of structures in the set. As a real-world test, we calculated the matrices for 6 different ensembles from NRW. The GPU method was 260 times faster that the fastest existing CPU based method and over 500 times faster than the method that had been previously used.

CONCLUSIONS

GPU-Q-J is a significant advance over previous CPU methods. It relieves a major bottleneck in the clustering of large numbers of structures for NRW. It also has applications in structure comparison methods that involve multiple superposition and RMSD determination steps, particularly when such methods are applied on a proteome and genome wide scale.