GenEpi: gene-based epistasis discovery using machine learning

BACKGROUND

Genome-wide association studies (GWAS) provide a powerful means to identify associations between genetic variants and phenotypes. However, GWAS techniques for detecting epistasis, the interactions between genetic variants associated with phenotypes, are still limited. We believe that developing an efficient and effective GWAS method to detect epistasis will be a key for discovering sophisticated pathogenesis, which is especially important for complex diseases such as Alzheimer's disease (AD).

RESULTS

In this regard, this study presents GenEpi, a computational package to uncover epistasis associated with phenotypes by the proposed machine learning approach. GenEpi identifies both within-gene and cross-gene epistasis through a two-stage modeling workflow. In both stages, GenEpi adopts two-element combinatorial encoding when producing features and constructs the prediction models by L1-regularized regression with stability selection. The simulated data showed that GenEpi outperforms other widely-used methods on detecting the ground-truth epistasis. As real data is concerned, this study uses AD as an example to reveal the capability of GenEpi in finding disease-related variants and variant interactions that show both biological meanings and predictive power.

CONCLUSIONS

The results on simulation data and AD demonstrated that GenEpi has the ability to detect the epistasis associated with phenotypes effectively and efficiently. The released package can be generalized to largely facilitate the studies of many complex diseases in the near future.

Interactions of Oral Molecular Excipients with Breast Cancer Resistance Protein, BCRP

Mechanistic-understanding-based selection of excipients may improve formulation development strategies for generic drug products and potentially accelerate their approval. Our study aimed at investigating the effects of molecular excipients present in orally administered FDA-approved drug products on the intestinal efflux transporter, BCRP (ABCG2), which plays a critical role in drug absorption with potential implications on drug safety and efficacy. We determined the interactions of 136 oral molecular excipients with BCRP in isolated membrane vesicles and identified 26 excipients as BCRP inhibitors with IC₅₀ values less than 5 μM using ³H-cholecystokinin octapeptide (³H-CCK8). These BCRP inhibitors belonged to three functional categories of excipients: dyes, surfactants, and flavoring agents. Compared with noninhibitors, BCRP inhibitors had significantly higher molecular weights and SLogP values. The inhibitory effects of excipients identified in membrane vesicles were also evaluated in BCRP-overexpressing HEK293 cells at similar concentrations. Only 1 of the 26 inhibitors of BCRP identified in vesicles inhibited BCRP-mediated ³H-oxypurinol uptake by more than 50%, consistent with the notion that BCRP inhibition depends on transmembrane or intracellular availability of the inhibitors. Collectively, the results of this study provide new information on excipient selection during the development of drug products with active pharmaceutical ingredients that are BCRP substrates.

Drug-nutrient interactions: discovering prescription drug inhibitors of the thiamine transporter ThTR-2 (SLC19A3)

BACKGROUND

Transporter-mediated drug-nutrient interactions have the potential to cause serious adverse events. However, unlike drug-drug interactions, these drug-nutrient interactions receive little attention during drug development. The clinical importance of drug-nutrient interactions was highlighted when a phase III clinical trial was terminated due to severe adverse events resulting from potent inhibition of thiamine transporter 2 (ThTR-2; SLC19A3).

OBJECTIVE

In this study, we tested the hypothesis that therapeutic drugs inhibit the intestinal thiamine transporter ThTR-2, which may lead to thiamine deficiency.

METHODS

For this exploration, we took a multifaceted approach, starting with a high-throughput in vitro primary screen to identify inhibitors, building in silico models to characterize inhibitors, and leveraging real-world data from electronic health records to begin to understand the clinical relevance of these inhibitors.

RESULTS

Our high-throughput screen of 1360 compounds, including many clinically used drugs, identified 146 potential inhibitors at 200 μM. Inhibition kinetics were determined for 28 drugs with half-maximal inhibitory concentration (IC50) values ranging from 1.03 μM to >1 mM. Several oral drugs, including metformin, were predicted to have intestinal concentrations that may result in ThTR-2-mediated drug-nutrient interactions. Complementary analysis using electronic health records suggested that thiamine laboratory values are reduced in individuals receiving prescription drugs found to significantly inhibit ThTR-2, particularly in vulnerable populations (e.g., individuals with alcoholism).

CONCLUSIONS

Our comprehensive analysis of prescription drugs suggests that several marketed drugs inhibit ThTR-2, which may contribute to thiamine deficiency, especially in at-risk populations.

Unraveling the functional role of the orphan solute carrier, SLC22A24 in the transport of steroid conjugates through metabolomic and genome-wide association studies

Variation in steroid hormone levels has wide implications for health and disease. The genes encoding the proteins involved in steroid disposition represent key determinants of interindividual variation in steroid levels and ultimately, their effects. Beginning with metabolomic data from genome-wide association studies (GWAS), we observed that genetic variants in the orphan transporter, SLC22A24 were significantly associated with levels of androsterone glucuronide and etiocholanolone glucuronide (sentinel SNPs p-value <1x10-30). In cells over-expressing human or various mammalian orthologs of SLC22A24, we showed that steroid conjugates and bile acids were substrates of the transporter. Phylogenetic, genomic, and transcriptomic analyses suggested that SLC22A24 has a specialized role in the kidney and appears to function in the reabsorption of organic anions, and in particular, anionic steroids. Phenome-wide analysis showed that functional variants of SLC22A24 are associated with human disease such as cardiovascular diseases and acne, which have been linked to dysregulated steroid metabolism. Collectively, these functional genomic studies reveal a previously uncharacterized protein involved in steroid homeostasis, opening up new possibilities for SLC22A24 as a pharmacological target for regulating steroid levels.

A Comprehensive Analysis of Ontogeny of Renal Drug Transporters: mRNA Analyses, Quantitative Proteomics, and Localization

Human renal membrane transporters play key roles in the disposition of renally cleared drugs and endogenous substrates, but their ontogeny is largely unknown. Using 184 human postmortem frozen renal cortical tissues (preterm newborns to adults) and a subset of 62 tissue samples, we measured the mRNA levels of 11 renal transporters and the transcription factor pregnane X receptor (PXR) with quantitative real‐time polymerase chain reaction, and protein abundance of nine transporters using liquid chromatography tandem mass spectrometry selective reaction monitoring, respectively. Expression levels of p‐glycoprotein, urate transporter 1, organic anion transporter 1, organic anion transporter 3, and organic cation transporter 2 increased with age. Protein levels of multidrug and toxin extrusion transporter 2‐K and breast cancer resistance protein showed no difference from newborns to adults, despite age‐related changes in mRNA expression. Multidrug and toxin extrusion transporter 1, glucose transporter 2, multidrug resistance‐associated protein 2, multidrug resistance‐associated protein 4 (MRP4), and PXR expression levels were stable. Using immunohistochemistry, we found that MRP4 localization in pediatric samples was similar to that in adult samples. Collectively, our study revealed that renal drug transporters exhibited different rates and patterns of maturation, suggesting that renal handling of substrates may change with age.

l-Type amino acid transporter 1 activity of 1,2,3-triazolyl analogs of l-histidine and l-tryptophan

A series of 1,2,3-triazole analogs of the amino acids l-histidine and l-tryptophan were modeled, synthesized and tested for l-type amino acid transporter 1 (LAT1; SLC7A5) activity to guide the design of amino acid-drug conjugates (prodrugs). These triazoles were conveniently prepared by the highly convergent Huisgen 1,3-dipolar cycloaddition (Click Chemistry). Despite comparable predicted binding modes, triazoles generally demonstrated reduced cell uptake and LAT1 binding potency relative to their natural amino acid counterparts. The structure-activity relationship (SAR) data for these triazoles has important ramifications for treating cancer and brain disorders using amino acid prodrugs or LAT1 inhibitors.

Genome-Wide Association and Functional Studies Reveal Novel Pharmacological Mechanisms for Allopurinol

Allopurinol, which lowers uric acid (UA) concentration, is increasingly being recognized for its benefits in cardiovascular and renal disease. However, response to allopurinol is variable. We gathered samples from 4,446 multiethnic subjects for a genome-wide association study of allopurinol response. Consistent with previous studies, we observed that the Q141K variant in ABCG2 (rs2231142), which encodes the efflux pump breast cancer resistance protein (BCRP), associated with worse response to allopurinol. However, for the first time this association reached genome-wide level significance (P = 8.06 × 10^-11 ). Additionally, we identified a novel association with a variant in GREM2 (rs1934341, P = 3.22 × 10^-6 ). In vitro studies identified oxypurinol, the active metabolite of allopurinol, as an inhibitor of the UA transporter GLUT9, suggesting that oxypurinol may modulate UA reabsorption. These results provide strong evidence for a role of BCRP Q141K in allopurinol response, and suggest that allopurinol may have additional hypouricemic effects beyond xanthine oxidase inhibition.

Functional and structural analysis of rare SLC2A2 variants associated with Fanconi-Bickel syndrome and metabolic traits

Deleterious variants in SLC2A2 cause Fanconi-Bickel Syndrome (FBS), a glycogen storage disorder, whereas less common variants in SLC2A2 associate with numerous metabolic diseases. Phenotypic heterogeneity in FBS has been observed, but its causes remain unknown. Our goal was to functionally characterize rare SLC2A2 variants found in FBS and metabolic disease-associated variants to understand the impact of these variants on GLUT2 activity and expression and establish genotype-phenotype correlations. Complementary RNA-injected Xenopus laevis oocytes were used to study mutant transporter activity and membrane expression. GLUT2 homology models were constructed for mutation analysis using GLUT1, GLUT3, and XylE as templates. Seventeen FBS variants were characterized. Only c.457_462delCTTATA (p.Leu153_Ile154del) exhibited residual glucose uptake. Functional characterization revealed that only half of the variants were expressed on the plasma membrane. Most less common variants (except c.593 C>A (p.Thr198Lys) and c.1087 G>T (p.Ala363Ser)) exhibited similar GLUT2 transport activity as the wild type. Structural analysis of GLUT2 revealed that variants affect substrate-binding, steric hindrance, or overall transporter structure. The mutant transporter that is associated with a milder FBS phenotype, p.Leu153_Ile154del, retained transport activity. These results improve our overall understanding of the underlying causes of FBS and impact of GLUT2 function on various clinical phenotypes ranging from rare to common disease.

Organic Anion Transporter Polypeptide 1B1 Polymorphism Modulates the Extent of Drug-Drug Interaction and Associated Biomarker Levels in Healthy Volunteers

Understanding transporter‐mediated drug–drug interactions is an integral part of risk assessment in drug development. Recent studies support the use of hexadecanedioate (HDA), tetradecanedioate (TDA), coproporphyrin (CP)‐I, and CP‐III as clinical biomarkers for evaluating organic anion‐transporting polypeptide (OATP)1B1 (SLCO1B1) inhibition. The current study investigated the effect of OATP1B1 genotype c.521T>C (OATP1B1‐Val174Ala) on the extent of interaction between cyclosporin A (CsA) and pravastatin, and associated endogenous biomarkers of the transporter (HDA, TDA, CP‐I, and CP‐III), in 20 healthy volunteers. The results show that the levels of each clinical biomarker and pravastatin were significantly increased in plasma samples of the volunteers following administration of pravastatin plus CsA compared with pravastatin plus placebo. The overall fold change in the area under the concentration–time curve (AUC) and maximum plasma concentration (C_max) was similar among the four biomarkers (1.8–2.5‐fold, paired t‐test P value < 0.05) in individuals who were homozygotes or heterozygotes of the major allele, c.521T. However, the fold change in AUC and C_max for HDA and TDA was significantly abolished in the subjects who were c.521‐CC, whereas the respective fold change in AUC and C_max for pravastatin and CP‐I and CP‐III were slightly weaker in individuals who were c.521‐CC compared with c.521‐TT/TC genotypes. In addition, this study provides the first evidence that SLCO1B1 c.521T>C genotype is significantly associated with CP‐I but not CP‐III levels. Overall, these results suggest that OATP1B1 genotype can modulate the effects of CsA on biomarker levels; the extent of modulation differs among the biomarkers.

Organic cation transporter 3 (Oct3) is a distinct catecholamines clearance route in adipocytes mediating the beiging of white adipose tissue

Beiging of white adipose tissue (WAT) is a particularly appealing target for therapeutics in the treatment of metabolic diseases through norepinephrine (NE)-mediated signaling pathways. Although previous studies report NE clearance mechanisms via SLC6A2 on sympathetic neurons or proinflammatory macrophages in adipose tissues (ATs), the low catecholamine clearance capacity of SLC6A2 may limit the cleaning efficiency. Here, we report that mouse organic cation transporter 3 (Oct3; Slc22a3) is highly expressed in WAT and displays the greatest uptake rate of NE as a selective non-neural route of NE clearance in white adipocytes, which differs from other known routes such as adjacent neurons or macrophages. We further show that adipocytes express high levels of NE degradation enzymes Maoa, Maob, and Comt, providing the molecular basis on NE clearance by adipocytes together with its reuptake transporter Oct3. Under NE administration, ablation of Oct3 induces higher body temperature, thermogenesis, and lipolysis compared with littermate controls. After prolonged cold challenge, inguinal WAT (ingWAT) in adipose-specific Oct3-deficient mice shows much stronger browning characteristics and significantly elevated expression of thermogenic and mitochondrial biogenesis genes than in littermate controls, and this response involves enhanced β-adrenergic receptor (β-AR)/protein kinase A (PKA)/cyclic adenosine monophosphate (cAMP)-responsive element binding protein (Creb) pathway activation. Glycolytic genes are reprogrammed to significantly higher levels to compensate for the loss of ATP production in adipose-specific Oct3 knockout (KO) mice, indicating the fundamental role of glucose metabolism during beiging. Inhibition of β-AR largely abolishes the higher lipolytic and thermogenic activities in Oct3-deficient ingWAT, indicating the NE overload in the vicinity of adipocytes in Oct3 KO adipocytes. Of note, reduced functional alleles in human OCT3 are also identified to be associated with increased basal metabolic rate (BMR). Collectively, our results demonstrate that Oct3 governs β-AR activity as a NE recycling transporter in white adipocytes, offering potential therapeutic applications for metabolic disorders.

Reevaluating the Substrate Specificity of the L-Type Amino Acid Transporter (LAT1)

The L-type amino acid transporter 1 (LAT1, SLC7A5) transports essential amino acids across the blood-brain barrier (BBB) and into cancer cells. To utilize LAT1 for drug delivery, potent amino acid promoieties are desired, as prodrugs must compete with millimolar concentrations of endogenous amino acids. To better understand ligand-transporter interactions that could improve potency, we developed structural LAT1 models to guide the design of substituted analogues of phenylalanine and histidine. Furthermore, we evaluated the structure-activity relationship (SAR) for both enantiomers of naturally occurring LAT1 substrates. Analogues were tested in cis-inhibition and trans-stimulation cell assays to determine potency and uptake rate. Surprisingly, LAT1 can transport amino acid-like substrates with wide-ranging polarities including those containing ionizable substituents. Additionally, the rate of LAT1 transport was generally nonstereoselective even though enantiomers likely exhibit different binding modes. Our findings have broad implications to the development of new treatments for brain disorders and cancer.

Transporters in Drug Development: 2018 ITC Recommendations for Transporters of Emerging Clinical Importance

This white paper provides updated International Transporter Consortium (ITC) recommendations on transporters that are important in drug development following the 3^rd ITC workshop. New additions include prospective evaluation of organic cation transporter 1 (OCT1) and retrospective evaluation of organic anion transporting polypeptide (OATP)2B1 because of their important roles in drug absorption, disposition, and effects. For the first time, the ITC underscores the importance of transporters involved in drug-induced vitamin deficiency (THTR2) and those involved in the disposition of biomarkers of organ function (OAT2 and bile acid transporters).

Genetic Variants in CPA6 and PRPF31 Are Associated With Variation in Response to Metformin in Individuals With Type 2 Diabetes

Metformin is the first-line treatment for type 2 diabetes (T2D). Although widely prescribed, the glucose-lowering mechanism for metformin is incompletely understood. Here, we used a genome-wide association approach in a diverse group of individuals with T2D from the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial to identify common and rare variants associated with HbA_1c response to metformin treatment and followed up these findings in four replication cohorts. Common variants in PRPF31 and CPA6 were associated with worse and better metformin response, respectively (P < 5 × 10^-6), and meta-analysis in independent cohorts displayed similar associations with metformin response (P = 1.2 × 10^-8 and P = 0.005, respectively). Previous studies have shown that PRPF31(+/-) knockout mice have increased total body fat (P = 1.78 × 10^-6) and increased fasted circulating glucose (P = 5.73 × 10^-6). Furthermore, rare variants in STAT3 associated with worse metformin response (q <0.1). STAT3 is a ubiquitously expressed pleiotropic transcriptional activator that participates in the regulation of metabolism and feeding behavior. Here, we provide novel evidence for associations of common and rare variants in PRPF31, CPA6, and STAT3 with metformin response that may provide insight into mechanisms important for metformin efficacy in T2D.

Influence of Transporter Polymorphisms on Drug Disposition and Response: A Perspective From the International Transporter Consortium

Advances in genomic technologies have led to a wealth of information identifying genetic polymorphisms in membrane transporters, specifically how these polymorphisms affect drug disposition and response. This review describes the current perspective of the International Transporter Consortium (ITC) on clinically important polymorphisms in membrane transporters. ITC suggests that, in addition to previously recommended polymorphisms in ABCG2 (BCRP) and SLCO1B1 (OATP1B1), polymorphisms in the emerging transporter, SLC22A1 (OCT1), be considered during drug development. Collectively, polymorphisms in these transporters are important determinants of interindividual differences in the levels, toxicities, and response to many drugs.

ITC Commentary on Metformin Clinical Drug-Drug Interaction Study Design That Enables an Efficacy- and Safety-Based Dose Adjustment Decision

Metformin drug-drug interaction (DDI) studies are conducted during development of drugs that inhibit organic cation transporters and/or multidrug and toxin extrusion proteins (OCTs/MATEs). Monitoring solely changes in systemic exposure, the typical DDI study endpoint appears inadequate for metformin, which is metabolically stable, has poor passive membrane permeability, and undergoes transporter-mediated tissue distribution and clearance. Evaluation of renal clearance, antihyperglycemic effects, and potentially lactate as an exploratory safety marker, can support rational metformin dose adjustment. The proposed DDI study design aims to adequately inform metformin dosing during comedication.

Pharmacogenetics of Antidiabetic Drugs

Pharmacogenetic studies of antidiabetic drugs have so far focused largely on response to metformin, which is the first-line therapy for treatment of type 2 diabetes (T2D). The first studies of metformin pharmacogenetics were focused on candidate genes that were implicated in metformin pharmacokinetics and transport. Since 2011, genome-wide association studies have been conducted in large cohorts of individuals with T2D identifying genes that are associated with glycemic response to metformin. There have been fewer pharmacogenetic studies of other antidiabetic drugs, and those have been largely limited to candidate gene studies with small sample sizes. Understanding the pharmacogenetics of antidiabetes medications is important for the integration of genetic screening into therapeutic decision making, and to achieve the goal of "precision medicine" for patients with T2D. In this chapter, we provide a review of the pharmacogenetics investigations of metformin and other antidiabetes medications. In addition, we highlight the importance of collaborative efforts with large sample size and representation from multiple ethnic groups in pharmacogenetics studies.