Computational biology approaches in drug repurposing and gene essentiality screening

Indiana University-Purdue University Indianapolis (IUPUI)The rapid innovations in biotechnology have led to an exponential growth of data and electronically accessible scientific literature. In this enormous scientific data, knowledge can be exploited, and novel discoveries can be made. In my dissertation, I have focused on the novel molecular mechanism and therapeutic discoveries from big data for complex diseases. It is very evident today that complex diseases have many factors including genetics and environmental effects. The discovery of these factors is challenging and critical in personalized medicine. The increasing cost and time to develop new drugs poses a new challenge in effectively treating complex diseases. In this dissertation, we want to demonstrate that the use of existing data and literature as a potential resource for discovering novel therapies and in repositioning existing drugs. The key to identifying novel knowledge is in integrating information from decades of research across the different scientific disciplines to uncover interactions that are not explicitly stated. This puts critical information at the fingertips of researchers and clinicians who can take advantage of this newly acquired knowledge to make informed decisions. This dissertation utilizes computational biology methods to identify and integrate existing scientific data and literature resources in the discovery of novel molecular targets and drugs that can be repurposed. In chapters 1 of my dissertation, I extensively sifted through scientific literature and identified a novel interaction between Vitamin A and CYP19A1 that could lead to a potential increase in the production of estrogens. Further in chapter 2 by exploring a microarray dataset from an estradiol gene sensitivity study I was able to identify a potential novel anti-estrogenic indication for the commonly used urinary analgesic, phenazopyridine. Both discoveries were experimentally validated in the laboratory. In chapter 3 of my dissertation, through the use of a manually curated corpus and machine learning algorithms, I identified and extracted genes that are essential for cell survival. These results brighten the reality that novel knowledge with potential clinical applications can be discovered from existing data and literature by integrating information across various scientific disciplines

A translational bioinformatic approach in identifying and validating an interaction between Vitamin A and CYP19A1

INTRODUCTION: One major challenge in personalized medicine research is to identify the environmental factors that can alter drug response, and to investigate their molecular mechanisms. These environmental factors include co-medications, food, and nutrition or dietary supplements. The increasing use of dietary supplements and their potential interactions with cytochrome P450 (CYP450) enzymes is a highly significant personalized medicine research domain, because most of the drugs on the market are metabolized through CYP450 enzymes. METHODS: Initial bioinformatics analysis revealed a number of regulators of CYP450 enzymes from a human liver bank gene expression quantitative loci data set. Then, a compound-gene network was constructed from the curated literature data. This network consisted of compounds that interact with either CYPs and/or their regulators that influence either their gene expression or activity. We further evaluated this finding in three different cell lines: JEG3, HeLa, and LNCaP cells. RESULTS: From a total of 868 interactions we were able to identify an interesting interaction between retinoic acid (i.e. Vitamin A) and the aromatase gene (i.e. CYP19A1). Our experimental results showed that retinoic acid at physiological concentration significantly influenced CYP19A1 gene expressions. CONCLUSIONS: These results suggest that the presence of retinoic acid may alter the efficacy of agents used to suppress aromatase expression

Common genetic polymorphisms of adenosine A2A receptor do not influence response to regadenoson

Aim: Hemodynamic response to regadenoson varies greatly, and underlying mechanisms for variability are poorly understood. We hypothesized that five common variants of adenosine A2A receptor (ADORA2A) are associated with altered response to regadenoson. Methods: Consecutive subjects (n = 357) undergoing resting regadenoson nuclear stress imaging were enrolled. Genotyping was performed using Taqman-based assays for rs5751862, rs2298383, rs3761422, rs2267076 and rs5751876. Results: There was no significant difference in heart rate or blood pressure between different genotypes following regadenoson administration. There was also no significant difference in myocardial ischemia detected by nuclear perfusion imaging as defined by summed difference score, or in self-reported side effects among the genotypes tested. Conclusion: The common A2A variants studied are not associated with variability in hemodynamic response to regadenoson or variability in detection of ischemia with nuclear perfusion stress imaging

AMPD1 polymorphism and response to regadenoson

AIMS: AMPD1 c.34C > T (rs17602729) polymorphism results in AMPD1 deficiency. We examined the association of AMPD1 deficiency and variability of hemodynamic response to regadenoson. SUBJECTS & METHODS: Genotyping for c.34C>T was performed in 267 patients undergoing regadenoson cardiac stress testing. RESULTS: Carriers of c.34C >T variant exhibited higher relative changes in systolic blood pressure (SBP) compared with wild-type subjects ([%] SBP change to peak: 12 ± 25 vs 5 ± 13%; p = 0.01) ([%] SBP change to nadir: -3 ± 15 vs -7 ± 11%; p = 0.04). Change in heart rate was similar between groups, but side effects were more common in carriers of the variant (+LR = 4.2; p = 0.04). CONCLUSION: AMPD1 deficiency may be involved in the modulation of regadenoson's systemic effects

Age-Related Changes in MicroRNA Expression and Pharmacogenes in Human Liver

Developmental changes in the liver can significantly impact drug disposition. Due to the emergence of microRNAs (miRNAs) as important regulators of drug disposition gene expression, we studied age-dependent changes in miRNA expression. Expression of 533 miRNAs was measured in 90 human liver tissues (fetal, pediatric [1-17 years], and adult [28-80 years]; n = 30 each). In all, 114 miRNAs were upregulated and 72 were downregulated from fetal to pediatric, and 2 and 3, respectively, from pediatric to adult. Among the developmentally changing miRNAs, 99 miRNA-mRNA interactions were predicted or experimentally validated (e.g., hsa-miR-125b-5p-CYP1A1; hsa-miR-34a-5p-HNF4A). In human liver samples (n = 10 each), analyzed by RNA-sequencing, significant negative correlations were observed between the expression of >1,000 miRNAs and mRNAs of drug disposition and regulatory genes. Our data suggest a mechanism for the marked changes in hepatic gene expression between the fetal and pediatric developmental periods, and support a role for these age-dependent miRNAs in regulating drug disposition

Whole Genome Amplification of DNA for Genotyping Pharmacogenetics Candidate Genes

Whole genome amplification (WGA) technologies can be used to amplify genomic DNA when only small amounts of DNA are available. The Multiple Displacement Amplification Phi polymerase based amplification has been shown to accurately amplify DNA for a variety of genotyping assays; however, it has not been tested for genotyping many of the clinically relevant genes important for pharmacogenetic studies, such as the cytochrome P450 genes, that are typically difficult to genotype due to multiple pseudogenes, copy number variations, and high similarity to other related genes. We evaluated whole genome amplified samples for Taqman™ genotyping of SNPs in a variety of pharmacogenetic genes. In 24 DNA samples from the Coriell human diversity panel, the call rates, and concordance between amplified (∼200-fold amplification) and unamplified samples was 100% for two SNPs in CYP2D6 and one in ESR1. In samples from a breast cancer clinical trial (Trial 1), we compared the genotyping results in samples before and after WGA for three SNPs in CYP2D6, one SNP in CYP2C19, one SNP in CYP19A1, two SNPs in ESR1, and two SNPs in ESR2. The concordance rates were all >97%. Finally, we compared the allele frequencies of 143 SNPs determined in Trial 1 (whole genome amplified DNA) to the allele frequencies determined in unamplified DNA samples from a separate trial (Trial 2) that enrolled a similar population. The call rates and allele frequencies between the two trials were 98 and 99.7%, respectively. We conclude that the whole genome amplified DNA is suitable for Taqman™ genotyping for a wide variety of pharmacogenetically relevant SNPs

Identification of germline cancer predisposition variants during clinical ctDNA testing

Next-generation sequencing of circulating tumor DNA (ctDNA) is a non-invasive method to guide therapy selection for cancer patients. ctDNA variant allele frequency (VAF) is commonly reported and may aid in discerning whether a variant is germline or somatic. We report on the fidelity of VAF in ctDNA as a predictor for germline variant carriage. Two patient cohorts were studied. Cohort 1 included patients with known germline variants. Cohort 2 included patients with any variant detected by the ctDNA assay with VAF of 40–60%. In cohort 1, 36 of 91 (40%) known germline variants were identified through ctDNA analysis with a VAF of 39–87.6%. In cohort 2, 111 of 160 (69%) variants identified by ctDNA analysis with a VAF between 40 and 60% were found to be germline. Therefore, variants with a VAF between 40 and 60% should induce suspicion for germline status but should not be used as a replacement for germline testing

Functional characterization of a genetic polymorphism in the promoter of the ESR2 gene

The ESR2 gene encodes the estrogen receptor beta protein. Several studies have shown that genetic variants in the ESR2 gene are associated with a variety of clinical phenotypes. However, very little is known about the functional significance of ESR2 genetic variants. We used a bioinformatics approach to identify regions of the ESR2 promoter that is evolutionarily conserved across the genomes of several species. We resequenced 1.6 kb of the ESR2 gene which included 0.8 kb of the promoter, 0.3 kb of exon ON, and 0.5 kb of the following intron. We identified five single-nucleotide polymorphisms (SNPs) in the ESR2 promoter and one SNP in the intron. Phase analysis indicated that the SNPs likely exist in 11 different haplotypes. Three of the SNPs (rs8008187, rs3829768, rs35036378) were predicted to alter transcription factor binding sites in the ESR2 promoter. All three were detected only in African American subjects. The rs35036378 SNP was in the TATA box and was highly conserved across species. ESR2 promoter reporter assays in LNCaP and SKBR3 cell lines showed that the variant construct containing the rs35036378 SNP allele had approximately 50% less activity relative to the wild-type construct. We conclude that the rs35036378 SNP appears to cause a reduced promoter activity of the ESR2 gene

Implications of Incidental Germline Findings Identified In the Context of Clinical Whole Exome Sequencing for Guiding Cancer Therapy

PURPOSE Identification of incidental germline mutations in the context of next-generation sequencing is an unintended consequence of advancing technologies. These data are critical for family members to understand disease risks and take action. PATIENTS AND METHODS A retrospective cohort analysis was conducted of 1,028 adult patients with metastatic cancer who were sequenced with tumor and germline whole exome sequencing (WES). Germline variant call files were mined for pathogenic/likely pathogenic (P/LP) variants using the ClinVar database and narrowed to high-quality submitters. RESULTS Median age was 59 years, with 16% of patients ≤ 45 years old. The most common tumor types were breast cancer (12.5%), colorectal cancer (11.5%), sarcoma (9.3%), prostate cancer (8.4%), and lung cancer (6.6%). We identified 3,427 P/LP variants in 471 genes, and 84% of patients harbored one or more variant. One hundred thirty-two patients (12.8%) carried a P/LP variant in a cancer predisposition gene, with BRCA2 being the most common (1.6%). Patients with breast cancer were most likely to carry a P/LP variant (19.2%). One hundred ten patients (10.7%) carried a P/LP variant in a gene that would be recommended by the American College of Medical Genetics and Genomics to be reported as a result of clinical actionability, with the most common being ATP7B (2.7%), BRCA2 (1.6%), MUTYH (1.4%), and BRCA1 (1%). Of patients who carried a P/LP variant in a cancer predisposition gene, only 53% would have been offered correct testing based on current clinical practice guidelines. Of 471 mutated genes, 231 genes had a P/LP variant identified in one patient, demonstrating significant genetic heterogeneity. CONCLUSION The majority of patients undergoing clinical cancer WES harbor a pathogenic germline variation. Identification of clinically actionable germline findings will create additional burden on oncology clinics as broader WES becomes common

Clinical Opportunities for Germline Pharmacogenetics and Management of Drug-Drug Interactions in Patients With Advanced Solid Cancers.

PURPOSE: Precision medicine approaches, including germline pharmacogenetics (PGx) and management of drug-drug interactions (DDIs), are likely to benefit patients with advanced cancer who are frequently prescribed multiple concomitant medications to treat cancer and associated conditions. Our objective was to assess the potential opportunities for PGx and DDI management within a cohort of adults with advanced cancer. METHODS: Medication data were collected from the electronic health records for 481 subjects since their first cancer diagnosis. All subjects were genotyped for variants with clinically actionable recommendations in Clinical Pharmacogenetics Implementation Consortium guidelines for 13 pharmacogenes. DDIs were defined as concomitant prescription of strong inhibitors or inducers with sensitive substrates of the same drug-metabolizing enzyme and were assessed for six major cytochrome P450 (CYP) enzymes. RESULTS: Approximately 60% of subjects were prescribed at least one medication with Clinical Pharmacogenetics Implementation Consortium recommendations, and approximately 14% of subjects had an instance for actionable PGx, defined as a prescription for a drug in a subject with an actionable genotype. The overall subject-level prevalence of DDIs and serious DDIs were 50.3% and 34.8%, respectively. Serious DDIs were most common for CYP3A, CYP2D6, and CYP2C19, occurring in 24.9%, 16.8%, and 11.7% of subjects, respectively. When assessing PGx and DDIs together, approximately 40% of subjects had at least one opportunity for a precision medicine-based intervention and approximately 98% of subjects had an actionable phenotype for at least one CYP enzyme. CONCLUSION: Our findings demonstrate numerous clinical opportunities for germline PGx and DDI management in adults with advanced cancer