Evolutionary Diagnosis of non-synonymous variants involved in differential drug response

Global spectrum of population-specific common missense variation in cytochrome P450 pharmacogenes

Next-generation sequencing technology has afforded the discovery of many novel variants that are of significance to inheritable pharmacogenomics (PGx) traits but a large proportion of them have unknown consequences. These include missense variants resulting in single amino acid substitutions in cytochrome P450 (CYP) proteins that can impair enzyme function, leading to altered drug efficacy and toxicity. While most unknown variants are rare, an overlooked minority are variants that are collectively rare but enriched in specific populations. Here, we analyzed sequence variation data in 141,456 individuals from across eight study populations in gnomAD for 38 CYP genes to identify such variants in addition to common variants. By further comparison with data from two PGx-specific databases (PharmVar and PharmGKB) and ClinVar, we identified 234 missense variants in 35 CYP genes, of which 107 were unknown to these databases. Most unknown variants (n = 83) were population-specific common variants and several (n = 7) were found in important CYP pharmacogenes (CYP2D6, CYP4F2, and CYP2C19). Overall, 29% (n = 31) of 107 unknown variants were predicted to affect CYP enzyme function although further biochemical characterization is necessary. These variants may elucidate part of the unexplained interpopulation differences observed in drug response.

Common Treatment, Common Variant: Evolutionary Prediction of Functional Pharmacogenomic Variants

Pharmacogenomics holds the promise of personalized drug efficacy optimization and drug toxicity minimization. Much of the research conducted to date, however, suffers from an ascertainment bias towards European participants. Here, we leverage publicly available, whole genome sequencing data collected from global populations, evolutionary characteristics, and annotated protein features to construct a new in silico machine learning pharmacogenetic identification method called XGB-PGX. When applied to pharmacogenetic data, XGB-PGX outperformed all existing prediction methods and identified over 2000 new pharmacogenetic variants. While there are modest pharmacogenetic allele frequency distribution differences across global population samples, the most striking distinction is between the relatively rare putatively neutral pharmacogene variants and the relatively common established and newly predicted functional pharamacogenetic variants. Our findings therefore support a focus on individual patient pharmacogenetic testing rather than on clinical presumptions about patient race, ethnicity, or ancestral geographic residence. We further encourage more attention be given to the impact of common variation on drug response and propose a new ‘common treatment, common variant’ perspective for pharmacogenetic prediction that is distinct from the types of variation that underlie complex and Mendelian disease. XGB-PGX has identified many new pharmacovariants that are present across all global communities; however, communities that have been underrepresented in genomic research are likely to benefit the most from XGB-PGX’s in silico predictions.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.

Biological relevance of computationally predicted pathogenicity of noncoding variants

Computational prediction of the phenotypic propensities of noncoding single nucleotide variants typically combines annotation of genomic, functional and evolutionary attributes into a single score. Here, we evaluate if the claimed excellent accuracies of these predictions translate into high rates of success in addressing questions important in biological research, such as fine mapping causal variants, distinguishing pathogenic allele(s) at a given position, and prioritizing variants for genetic risk assessment. A significant disconnect is found to exist between the statistical modelling and biological performance of predictive approaches. We discuss fundamental reasons underlying these deficiencies and suggest that future improvements of computational predictions need to address confounding of allelic, positional and regional effects as well as imbalance of the proportion of true positive variants in candidate lists.

Researchers can make use of a variety of computational tools to prioritize genetic variants and predict their pathogenicity. Here, the authors evaluate the performance of six of these tools in three typical biological tasks and find generally low concordance of predictions and experimental confirmation.

Use of Germline Genetic Variability for Prediction of Chemoresistance and Prognosis of Breast Cancer Patients

The aim of our study was to set up a panel for targeted sequencing of chemoresistance genes and the main transcription factors driving their expression and to evaluate their predictive and prognostic value in breast cancer patients. Coding and regulatory regions of 509 genes, selected from PharmGKB and Phenopedia, were sequenced using massive parallel sequencing in blood DNA from 105 breast cancer patients in the testing phase. In total, 18,245 variants were identified of which 2565 were novel variants (without rs number in dbSNP build 150) in the testing phase. Variants with major allele frequency over 0.05 were further prioritized for validation phase based on a newly developed decision tree. Using emerging in silico tools and pharmacogenomic databases for functional predictions and associations with response to cytotoxic therapy or disease-free survival of patients, 55 putative variants were identified and used for validation in 805 patients with clinical follow up using KASP^TM technology. In conclusion, associations of rs2227291, rs2293194, and rs4376673 (located in ATP7A, KCNAB1, and DFFB genes, respectively) with response to neoadjuvant cytotoxic therapy and rs1801160 in DPYD with disease-free survival of patients treated with cytotoxic drugs were validated and should be further functionally characterized.

DROIDS 1.20: A GUI-Based Pipeline for GPU-Accelerated Comparative Protein Dynamics

Traditional informatics in comparative genomics work only with static representations of biomolecules (i.e., sequence and structure), thereby ignoring the molecular dynamics (MD) of proteins that define function in the cell. A comparative approach applied to MD would connect this very short timescale process, defined in femtoseconds, to one of the longest in the universe: molecular evolution measured in millions of years. Here, we leverage advances in graphics-processing-unit-accelerated MD simulation software to develop a comparative method of MD analysis and visualization that can be applied to any two homologous Protein Data Bank structures. Our open-source pipeline, DROIDS (Detecting Relative Outlier Impacts in Dynamic Simulations), works in conjunction with existing molecular modeling software to convert any Linux gaming personal computer into a "comparative computational microscope" for observing the biophysical effects of mutations and other chemical changes in proteins. DROIDS implements structural alignment and Benjamini-Hochberg-corrected Kolmogorov-Smirnov statistics to compare nanosecond-scale atom bond fluctuations on the protein backbone, color mapping the significant differences identified in protein MD with single-amino-acid resolution. DROIDS is simple to use, incorporating graphical user interface control for Amber16 MD simulations, cpptraj analysis, and the final statistical and visual representations in R graphics and UCSF Chimera. We demonstrate that DROIDS can be utilized to visually investigate molecular evolution and disease-related functional changes in MD due to genetic mutation and epigenetic modification. DROIDS can also be used to potentially investigate binding interactions of pharmaceuticals, toxins, or other biomolecules in a functional evolutionary context as well.

A Computational Approach to Identify a Potential Alternative Drug With Its Positive Impact Toward PMP22

Mutations in the Peripheral Myelin Protein 22 (PMP22) leads to Charcot Marie Tooth type 1A (CMT1A, a subtype of CMT1) disease which is the most common inherited neuropathy of peripheral nervous system. In the present study, we used series of in silico prediction methods to screen and identify the most deleterious non-synonymous SNPs (nsSNPs) in PMP22 gene. Out of 48 nsSNPs, five nsSNPs (L16P, L19P, T23R, W28R, and L147R) associated with PMP22 were predicted to be highly deleterious and destabilizing the protein. To explore the possible structure-function relationship, we employed abinitio modeling strategy using the CABS-fold server to predict the three-dimensional structure models in the absence of crystallized structures in PMP22 protein. We used Cytoscape 3.4.0 plugin Integrated Complex Traits Networks interface (iCTNet) to identify the probable drug-gene interactions in PMP22 gene. A total of 22 chemical compounds yielded from the aforementioned tool was subjected to Molinspiration and OSIRIS program to screen and identify the potent drug molecules for further analysis. Five chemical compounds with excellent bioavailability and drug relevant property were selected for molecular docking simulation study. We modeled five mutant structures at their corresponding positions and performed molecular docking simulation analysis using AutoDock Tools (ADT) version 1.5.6 and ArgusLab 4.0.1 tools to analyze their interaction patterns and binding efficacy. Based on the results obtained from the computational study, we predict that estradiol could be a potential drug of choice for treating patients with CMT1A which needs larger attention from biologists in the near future. J. Cell. Biochem. 118: 3730-3743, 2017. © 2017 Wiley Periodicals, Inc.