The full spectrum of SLC22 OCT1 mutations illuminates the bridge between drug transporter biophysics and pharmacogenomics.

Mutations in transporters can impact an individual's response to drugs and cause many diseases. Few variants in transporters have been evaluated for their functional impact. Here, we combine saturation mutagenesis and multi-phenotypic screening to dissect the impact of 11,213 missense single-amino-acid deletions, and synonymous variants across the 554 residues of OCT1, a key liver xenobiotic transporter. By quantifying in parallel expression and substrate uptake, we find that most variants exert their primary effect on protein abundance, a phenotype not commonly measured alongside function. Using our mutagenesis results combined with structure prediction and molecular dynamic simulations, we develop accurate structure-function models of the entire transport cycle, providing biophysical characterization of all known and possible human OCT1 polymorphisms. This work provides a complete functional map of OCT1 variants along with a framework for integrating functional genomics, biophysical modeling, and human genetics to predict variant effects on disease and drug efficacy.

Functional genomics of OCTN2 variants informs protein-specific variant effect predictor for Carnitine Transporter Deficiency.

Genetic variants in SLC22A5, encoding the membrane carnitine transporter OCTN2, cause the rare metabolic disorder Carnitine Transporter Deficiency (CTD). CTD is potentially lethal but actionable if detected early, with confirmatory diagnosis involving sequencing of SLC22A5. Interpretation of missense variants of uncertain significance (VUSs) is a major challenge. In this study, we sought to characterize the largest set to date (n = 150) of OCTN2 variants identified in diverse ancestral populations, with the goals of furthering our understanding of the mechanisms leading to OCTN2 loss-of-function (LOF) and creating a protein-specific variant effect prediction model for OCTN2 function. Uptake assays with 14C-carnitine revealed that 105 variants (70%) significantly reduced transport of carnitine compared to wild-type OCTN2, and 37 variants (25%) severely reduced function to less than 20%. All ancestral populations harbored LOF variants; 62% of green fluorescent protein (GFP)-tagged variants impaired OCTN2 localization to the plasma membrane of human embryonic kidney (HEK293T) cells, and subcellular localization significantly associated with function, revealing a major LOF mechanism of interest for CTD. With these data, we trained a model to classify variants as functional (>20% function) or LOF (<20% function). Our model outperformed existing state-of-the-art methods as evaluated by multiple performance metrics, with mean area under the receiver operating characteristic curve (AUROC) of 0.895 ± 0.025. In summary, in this study we generated a rich dataset of OCTN2 variant function and localization, revealed important disease-causing mechanisms, and improved upon machine learning-based prediction of OCTN2 variant function to aid in variant interpretation in the diagnosis and treatment of CTD.

Opportunities and challenges for the computational interpretation of rare variation in clinically important genes.

Genome sequencing is enabling precision medicine-tailoring treatment to the unique constellation of variants in an individual's genome. The impact of recurrent pathogenic variants is often understood, however there is a long tail of rare genetic variants that are uncharacterized. The problem of uncharacterized rare variation is especially acute when it occurs in genes of known clinical importance with functionally consequential variants and associated mechanisms. Variants of uncertain significance (VUSs) in these genes are discovered at a rate that outpaces current ability to classify them with databases of previous cases, experimental evaluation, and computational predictors. Clinicians are thus left without guidance about the significance of variants that may have actionable consequences. Computational prediction of the impact of rare genetic variation is increasingly becoming an important capability. In this paper, we review the technical and ethical challenges of interpreting the function of rare variants in two settings: inborn errors of metabolism in newborns and pharmacogenomics. We propose a framework for a genomic learning healthcare system with an initial focus on early-onset treatable disease in newborns and actionable pharmacogenomics. We argue that (1) a genomic learning healthcare system must allow for continuous collection and assessment of rare variants, (2) emerging machine learning methods will enable algorithms to predict the clinical impact of rare variants on protein function, and (3) ethical considerations must inform the construction and deployment of all rare-variation triage strategies, particularly with respect to health disparities arising from unbalanced ancestry representation.

Deorphaning a solute carrier 22 family member, SLC22A15, through functional genomic studies.

The human solute carrier 22A (SLC22A) family consists of 23 members, representing one of the largest families in the human SLC superfamily. Despite their pharmacological and physiological importance in the absorption and disposition of a range of solutes, eight SLC22A family members remain classified as orphans. In this study, we used a multifaceted approach to identify ligands of orphan SLC22A15. Ligands of SLC22A15 were proposed based on phylogenetic analysis and comparative modeling. The putative ligands were then confirmed by metabolomic screening and uptake assays in SLC22A15 transfected HEK293 cells. Metabolomic studies and transporter assays revealed that SLC22A15 prefers zwitterionic compounds over cations and anions. We identified eight zwitterions, including ergothioneine, carnitine, carnosine, gabapentin, as well as four cations, including MPP+ , thiamine, and cimetidine, as substrates of SLC22A15. Carnosine was a specific substrate of SLC22A15 among the transporters in the SLC22A family. SLC22A15 transport of several substrates was sodium-dependent and exhibited a higher Km for ergothioneine, carnitine, and carnosine compared to previously identified transporters for these ligands. This is the first study to characterize the function of SLC22A15. Our studies demonstrate that SLC22A15 may play an important role in determining the systemic and tissue levels of ergothioneine, carnosine, and other zwitterions.

Understanding drug-drug interaction and pharmacogenomic changes in pharmacokinetics for metabolized drugs.

Here we characterize and summarize the pharmacokinetic changes for metabolized drugs when drug-drug interactions and pharmacogenomic variance are observed. Following multiple dosing to steady-state, oral systemic concentration-time curves appear to follow a one-compartment body model, with a shorter rate limiting half-life, often significantly shorter than the single dose terminal half-life. This simplified disposition model at steady-state allows comparisons of measurable parameters (i.e., area under the curve, half-life, maximum concentration and time to maximum concentration) following drug interaction or pharmacogenomic variant studies to be utilized to characterize whether a drug is low versus high hepatic extraction ratio, even without intravenous dosing. The characteristics of drugs based on the ratios of area under the curve, maximum concentration and half-life are identified with recognition that volume of distribution is essentially unchanged for drug interaction and pharmacogenomic variant studies where only metabolic outcomes are changed and transporters are not significantly involved. Comparison of maximum concentration changes following single dose interaction and pharmacogenomic variance studies may also identify the significance of intestinal first pass changes. The irrelevance of protein binding changes on pharmacodynamic outcomes following oral and intravenous dosing of low hepatic extraction ratio drugs, versus its relevance for high hepatic extraction ratio drugs is re-emphasized.

Building a Predictive PBPK Model for Human OATP Substrates: a Strategic Framework for Early Evaluation of Clinical Pharmacokinetic Variations Using Pitavastatin as an Example.

To select a drug candidate for clinical development, accurately and promptly predicting human pharmacokinetic (PK) profiles, assessing drug-drug interactions (DDIs), and anticipating potential PK variations in disease populations are crucial steps in drug discovery. The complexity of predicting human PK significantly increases when hepatic transporters are involved in drug clearance (CL) and volume of distribution (Vss). A strategic framework is developed here, utilizing pitavastatin as an example. The framework includes the construction of a monkey physiologically-based PK (PBPK) model, model calibration to obtain scaling factors (SF) of in vitro-in vivo extrapolation (IVIVE) for various clearance parameters, human model development and validation, and assessment of DDIs and PK variations in disease populations. Through incorporating in vitro human parameters and calibrated SFs from the monkey model of 3.45, 0.14, and 1.17 for CLint,active, CLint,passive, and CLint,bile, respectively, and together with the relative fraction transported by individual transporters obtained from in vitro studies and the optimized Ki values for OATP inhibition, the model reasonably captured observed pitavastatin PK profiles, DDIs and PK variations in human subjects carrying genetic polymorphisms, i.e., AUC within 20%. Lastly, when applying the functional reduction based on measured OATP1B biomarkers, the model adequately predicted PK changes in the hepatic impairment population. The present study presents a strategic framework for early-stage drug development, enabling the prediction of PK profiles and assessment of PK variations in scenarios like DDIs, genetic polymorphism, and hepatic impairment-related disease states, specifically focusing on OATP substrates.

Effect of Antioxidants in Medicinal Products on Intestinal Drug Transporters.

The presence of mutagenic and carcinogenic N-nitrosamine impurities in medicinal products poses a safety risk. While incorporating antioxidants in formulations is a potential mitigation strategy, concerns arise regarding their interference with drug absorption by inhibiting intestinal drug transporters. Our study screened thirty antioxidants for inhibitory effects on key intestinal transporters-OATP2B1, P-gp, and BCRP in HEK-293 cells (OATP2B1) or membrane vesicles (P-gp, BCRP) using 3H-estrone sulfate, 3H-N-methyl quinidine, and 3H-CCK8 as substrates, respectively. The screen identified that butylated hydroxyanisole (BHA) and carnosic acid inhibited all three transporters (OATP2B1, P-gp, and BCRP), while ascorbyl palmitate (AP) inhibited OATP2B1 by more than 50%. BHA had IC50 values of 71 ± 20 µM, 206 ± 14 µM, and 182 ± 49 µM for OATP2B1, BCRP, and P-gp, respectively. AP exhibited IC50 values of 23 ± 10 µM for OATP2B1. The potency of AP and BHA was tested with valsartan, an OATP2B1 substrate, and revealed IC50 values of 26 ± 17 µM and 19 ± 11 µM, respectively, in HEK-293-OATP2B1 cells. Comparing IC50 values of AP and BHA with estimated intestinal concentrations suggests an unlikely inhibition of intestinal transporters at clinical concentrations of drugs formulated with antioxidants.

Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates.

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17ß-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

The full spectrum of OCT1 (SLC22A1) mutations bridges transporter biophysics to drug pharmacogenomics.

Membrane transporters play a fundamental role in the tissue distribution of endogenous compounds and xenobiotics and are major determinants of efficacy and side effects profiles. Polymorphisms within these drug transporters result in inter-individual variation in drug response, with some patients not responding to the recommended dosage of drug whereas others experience catastrophic side effects. For example, variants within the major hepatic Human organic cation transporter OCT1 (SLC22A1) can change endogenous organic cations and many prescription drug levels. To understand how variants mechanistically impact drug uptake, we systematically study how all known and possible single missense and single amino acid deletion variants impact expression and substrate uptake of OCT1. We find that human variants primarily disrupt function via folding rather than substrate uptake. Our study revealed that the major determinants of folding reside in the first 300 amino acids, including the first 6 transmembrane domains and the extracellular domain (ECD) with a stabilizing and highly conserved stabilizing helical motif making key interactions between the ECD and transmembrane domains. Using the functional data combined with computational approaches, we determine and validate a structure-function model of OCT1s conformational ensemble without experimental structures. Using this model and molecular dynamic simulations of key mutants, we determine biophysical mechanisms for how specific human variants alter transport phenotypes. We identify differences in frequencies of reduced function alleles across populations with East Asians vs European populations having the lowest and highest frequency of reduced function variants, respectively. Mining human population databases reveals that reduced function alleles of OCT1 identified in this study associate significantly with high LDL cholesterol levels. Our general approach broadly applied could transform the landscape of precision medicine by producing a mechanistic basis for understanding the effects of human mutations on disease and drug response.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

New and Emerging Research on Solute Carrier and ATP Binding Cassette Transporters in Drug Discovery and Development: Outlook From the International Transporter Consortium.

Enabled by a plethora of new technologies, research in membrane transporters has exploded in the past decade. The goal of this state-of-the-art article is to describe recent advances in research on membrane transporters that are particularly relevant to drug discovery and development. This review covers advances in basic, translational, and clinical research that has led to an increased understanding of membrane transporters at all levels. At the basic level, we describe the available crystal structures of membrane transporters in both the solute carrier (SLC) and ATP binding cassette superfamilies, which has been enabled by the development of cryogenic electron microscopy methods. Next, we describe new research on lysosomal and mitochondrial transporters as well as recently deorphaned transporters in the SLC superfamily. The translational section includes a summary of proteomic research, which has led to a quantitative understanding of transporter levels in various cell types and tissues and new methods to modulate transporter function, such as allosteric modulators and targeted protein degraders of transporters. The section ends with a review of the effect of the gut microbiome on modulation of transporter function followed by a presentation of 3D cell cultures, which may enable in vivo predictions of transporter function. In the clinical section, we describe new genomic and pharmacogenomic research, highlighting important polymorphisms in transporters that are clinically relevant to many drugs. Finally, we describe new clinical tools, which are becoming increasingly available to enable precision medicine, with the application of tissue-derived small extracellular vesicles and real-world biomarkers.

Differential toxicity to murine small and large intestinal epithelium induced by oncology drugs.

Gastrointestinal toxicity is a major concern in the development of drugs. Here, we establish the ability to use murine small and large intestine-derived monolayers to screen drugs for toxicity. As a proof-of-concept, we applied this system to assess gastrointestinal toxicity of ~50 clinically used oncology drugs, encompassing diverse mechanisms of action. Nearly all tested drugs had a deleterious effect on the gut, with increased sensitivity in the small intestine. The identification of differential toxicity between the small and large intestine enabled us to pinpoint differences in drug uptake (antifolates), drug metabolism (cyclophosphamide) and cell signaling (EGFR inhibitors) across the gut. These results highlight an under-appreciated distinction between small and large intestine toxicity and suggest distinct tissue properties important for modulating drug-induced gastrointestinal toxicity. The ability to accurately predict where and how drugs affect the murine gut will accelerate preclinical drug development.

The Effects of Genetic Mutations and Drugs on the Activity of the Thiamine Transporter, SLC19A2.

A rare cause of megaloblastic anemia (MA) is thiamine-responsive megaloblastic anemia (TRMA), a genetic disorder caused by mutations in SLC19A2 (encoding THTR1), a thiamine transporter. The study objectives were to (1) functionally characterize selected TRMA-associated SLC19A2 variants and (2) determine whether current prescription drugs associated with drug-induced MA (DIMA) may act via inhibition of SLC19A2. Functional characterization of selected SLC19A2 variants was performed by confocal microscopy and isotopic uptake studies of [3H]-thiamine in HEK293 cells. Sixty-three drugs associated with DIMA were screened for SLC19A2 inhibition in isotopic uptake studies. Three previously uncharacterized SLC19A2 variants identified in TRMA patients exhibited disrupted localization to the plasma membrane along with near-complete loss-of-function. Ten of 63 drugs inhibited SLC19A2-mediated thiamine transport ≥ 50% at screening concentrations; however, with the exception of erythromycin, none was predicted to inhibit SLC19A2 at clinically relevant unbound plasma concentrations. Data from electronic health records revealed reduced levels of thiamine pyrophosphate (TPP) in patients prescribed erythromycin, consistent with inhibition of SLC19A2-mediated thiamine transport. Here, we confirmed the role of three SLC19A2 variants in TRMA pathology. Additionally, we report that inhibition of SLC19A2 is a potential, but uncommon mechanism for DIMA.

Drugs in COVID-19 Clinical Trials: Predicting Transporter-Mediated Drug-Drug Interactions Using In Vitro Assays and Real-World Data.

Numerous drugs are currently under accelerated clinical investigation for the treatment of coronavirus disease 2019 (COVID-19); however, well-established safety and efficacy data for these drugs are limited. The goal of this study was to predict the potential of 25 small molecule drugs in clinical trials for COVID-19 to cause clinically relevant drug-drug interactions (DDIs), which could lead to potential adverse drug reactions (ADRs) with the use of concomitant medications. We focused on 11 transporters, which are targets for DDIs. In vitro potency studies in membrane vesicles or HEK293 cells expressing the transporters coupled with DDI risk assessment methods revealed that 20 of the 25 drugs met the criteria from regulatory authorities to trigger consideration of a DDI clinical trial. Analyses of real-world data from electronic health records, including a database representing nearly 120,000 patients with COVID-19, were consistent with several of the drugs causing transporter-mediated DDIs (e.g., sildenafil, chloroquine, and hydroxychloroquine). This study suggests that patients with COVID-19, who are often older and on various concomitant medications, should be carefully monitored for ADRs. Future clinical studies are needed to determine whether the drugs that are predicted to inhibit transporters at clinically relevant concentrations, actually result in DDIs.

Tetrodotoxin-Sensitive Neuronal-Type Na+ Channels: A Novel and Druggable Target for Prevention of Atrial Fibrillation.

Background Atrial fibrillation (AF) is a comorbidity associated with heart failure and catecholaminergic polymorphic ventricular tachycardia. Despite the Ca2+-dependent nature of both of these pathologies, AF often responds to Na+ channel blockers. We investigated how targeting interdependent Na+/Ca2+ dysregulation might prevent focal activity and control AF. Methods and Results We studied AF in 2 models of Ca2+-dependent disorders, a murine model of catecholaminergic polymorphic ventricular tachycardia and a canine model of chronic tachypacing-induced heart failure. Imaging studies revealed close association of neuronal-type Na+ channels (nNav) with ryanodine receptors and Na+/Ca2+ exchanger. Catecholamine stimulation induced cellular and in vivo atrial arrhythmias in wild-type mice only during pharmacological augmentation of nNav activity. In contrast, catecholamine stimulation alone was sufficient to elicit atrial arrhythmias in catecholaminergic polymorphic ventricular tachycardia mice and failing canine atria. Importantly, these were abolished by acute nNav inhibition (tetrodotoxin or riluzole) implicating Na+/Ca2+ dysregulation in AF. These findings were then tested in 2 nonrandomized retrospective cohorts: an amyotrophic lateral sclerosis clinic and an academic medical center. Riluzole-treated patients adjusted for baseline characteristics evidenced significantly lower incidence of arrhythmias including new-onset AF, supporting the preclinical results. Conclusions These data suggest that nNaVs mediate Na+-Ca2+ crosstalk within nanodomains containing Ca2+ release machinery and, thereby, contribute to AF triggers. Disruption of this mechanism by nNav inhibition can effectively prevent AF arising from diverse causes.