Structural and molecular basis of choline uptake into the brain by FLVCR2.

Choline is an essential nutrient that the human body needs in vast quantities for cell membrane synthesis, epigenetic modification and neurotransmission. The brain has a particularly high demand for choline, but how it enters the brain remains unknown1-3. The major facilitator superfamily transporter FLVCR1 (also known as MFSD7B or SLC49A1) was recently determined to be a choline transporter but is not highly expressed at the blood-brain barrier, whereas the related protein FLVCR2 (also known as MFSD7C or SLC49A2) is expressed in endothelial cells at the blood-brain barrier4-7. Previous studies have shown that mutations in human Flvcr2 cause cerebral vascular abnormalities, hydrocephalus and embryonic lethality, but the physiological role of FLVCR2 is unknown4,5. Here we demonstrate both in vivo and in vitro that FLVCR2 is a BBB choline transporter and is responsible for the majority of choline uptake into the brain. We also determine the structures of choline-bound FLVCR2 in both inward-facing and outward-facing states using cryo-electron microscopy. These results reveal how the brain obtains choline and provide molecular-level insights into how FLVCR2 binds choline in an aromatic cage and mediates its uptake. Our work could provide a novel framework for the targeted delivery of therapeutic agents into the brain.

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.

Rosace: a robust deep mutational scanning analysis framework employing position and mean-variance shrinkage.

Deep mutational scanning (DMS) measures the effects of thousands of genetic variants in a protein simultaneously. The small sample size renders classical statistical methods ineffective. For example, p-values cannot be correctly calibrated when treating variants independently. We propose Rosace, a Bayesian framework for analyzing growth-based DMS data. Rosace leverages amino acid position information to increase power and control the false discovery rate by sharing information across parameters via shrinkage. We also developed Rosette for simulating the distributional properties of DMS. We show that Rosace is robust to the violation of model assumptions and is more powerful than existing tools.

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 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.

Membrane transporters in drug development and as determinants of precision medicine.

The effect of membrane transporters on drug disposition, efficacy and safety is now well recognized. Since the initial publication from the International Transporter Consortium, significant progress has been made in understanding the roles and functions of transporters, as well as in the development of tools and models to assess and predict transporter-mediated activity, toxicity and drug-drug interactions (DDIs). Notable advances include an increased understanding of the effects of intrinsic and extrinsic factors on transporter activity, the application of physiologically based pharmacokinetic modelling in predicting transporter-mediated drug disposition, the identification of endogenous biomarkers to assess transporter-mediated DDIs and the determination of the cryogenic electron microscopy structures of SLC and ABC transporters. This article provides an overview of these key developments, highlighting unanswered questions, regulatory considerations and future directions.

Genome-Wide Association Study Identifies Pharmacogenomic Variants Associated With Metformin Glycemic Response in African American Patients With Type 2 Diabetes.

OBJECTIVE: Metformin is the most common treatment for type 2 diabetes (T2D). However, there have been no pharmacogenomic studies for T2D in which a population of color was used in the discovery analysis. This study sought to identify genomic variants associated with metformin response in African American patients with diabetes. RESEARCH DESIGN AND METHODS: Patients in the discovery set were adult, African American participants from the Diabetes Multi-omic Investigation of Drug Response (DIAMOND), a cohort study of patients with T2D from a health system serving southeast Michigan. DIAMOND participants had genome-wide genotype data and longitudinal electronic records of laboratory results and medication fills. The genome-wide discovery analysis identified polymorphisms correlated to changes in glycated hemoglobin (HbA1c) levels among individuals on metformin monotherapy. Lead associations were assessed for replication in an independent cohort of African American participants from Kaiser Permanente Northern California (KPNC) and in European American participants from DIAMOND. RESULTS: The discovery set consisted of 447 African American participants, whereas the replication sets included 353 African American KPNC participants and 466 European American DIAMOND participants. The primary analysis identified a variant, rs143276236, in the gene ARFGEF3, which met the threshold for genome-wide significance, replicated in KPNC African Americans, and was still significant in the meta-analysis (P = 1.17 × 10-9). None of the significant discovery variants replicated in European Americans DIAMOND participants. CONCLUSIONS: We identified a novel and biologically plausible genetic variant associated with a change in HbA1c levels among African American patients on metformin monotherapy. These results highlight the importance of diversity in pharmacogenomic studies.

Cardiac contraction and relaxation are regulated by distinct subcellular cAMP pools.

Cells interpret a variety of signals through G-protein-coupled receptors (GPCRs) and stimulate the generation of second messengers such as cyclic adenosine monophosphate (cAMP). A long-standing puzzle is deciphering how GPCRs elicit different physiological responses despite generating similar levels of cAMP. We previously showed that some GPCRs generate cAMP from both the plasma membrane and the Golgi apparatus. Here we demonstrate that cardiomyocytes distinguish between subcellular cAMP inputs to elicit different physiological outputs. We show that generating cAMP from the Golgi leads to the regulation of a specific protein kinase A (PKA) target that increases the rate of cardiomyocyte relaxation. In contrast, cAMP generation from the plasma membrane activates a different PKA target that increases contractile force. We further validated the physiological consequences of these observations in intact zebrafish and mice. Thus, we demonstrate that the same GPCR acting through the same second messenger regulates cardiac contraction and relaxation dependent on its subcellular location.

SLCO1B1 functional variants and statin-induced myopathy in people with recent genealogical ancestors from Africa: a population-based real-world study.

Background: Clinical pharmacogenetic implementation guidelines for statin therapy are derived from evidence of primarily Eurocentric study populations. Functional SLCO1B1 variants that are rare in these study populations have not been investigated as a determinant of statin myotoxicity and are thus missing from guideline inclusion. Objective: Determine the relationship between candidate functional SLCO1B1 variants and statin-induced myopathy in people with recent genealogical ancestors from Africa. Design: Population-based pharmacogenetic study using real-world evidence from electronic health record-linked biobanks. Setting: Various health care settings. Participants: Self-identified white and Black statin users with genome-wide genotyping data available. Measurements: Primarily, the odds of statin-induced myopathy + rhabdomyolysis. Secondarily, total bilirubin levels. Thirdly, cell-based functional assay results. Results: Meta-analyses results demonstrated an increased risk of statin-induced myopathy + rhabdomyolysis with c.481+1G>T (odds ratio [OR] = 3.27, 95% confidence interval [CI] 1.43-7.46, P =.005) and c.1463G>C (OR = 2.45, 95% CI 1.04-5.78, P =.04) for Black participants. For White participants, c.521T>C was also significantly associated with increased risk of statin-induced myopathy + rhabdomyolysis (OR = 1.41, 95% CI 1.20-1.67, P =5.4x10 -5 ). This effect size for c.521T>C was similar in the Black participants, but did not meet the level of statistical significance (OR = 1.47, 95% CI 0.58-3.73, P =0.41). Supporting evidence using total bilirubin as an endogenous biomarker of SLCO1B1 function as well as from cell-based functional studies corroborated these findings. Limitations: Data limited to severe statin myotoxicity events. Conclusion: Our findings implicate Afrocentric SLCO1B1 variants on preemptive pharmacogenetic testing panels, which could have an instant impact on reducing the risk of statin-associated myotoxicity in historically excluded groups. Primary Funding Source: National Institutes of Health, Office of the Director - All of Us (OD-AoURP).

Structural and molecular basis of choline uptake into the brain by FLVCR2.

Choline is an essential nutrient that the human body needs in vast quantities for cell membrane synthesis, epigenetic modification, and neurotransmission. The brain has a particularly high demand for choline, but how it enters the brain has eluded the field for over fifty years. The MFS transporter FLVCR1 was recently determined to be a choline transporter, and while this protein is not highly expressed at the blood-brain barrier (BBB), its relative FLVCR2 is. Previous studies have shown that mutations in human Flvcr2 cause cerebral vascular abnormalities, hydrocephalus, and embryonic lethality, but the physiological role of FLVCR2 is unknown. Here, we demonstrate both in vivo and in vitro that FLVCR2 is a BBB choline transporter and is responsible for the majority of choline uptake into the brain. We also determine the structures of choline-bound FLVCR2 in the inward- and outward-facing states using cryo-electron microscopy to 2.49 and 2.77 Å resolution, respectively. These results reveal how the brain obtains choline and provide molecular-level insights into how FLVCR2 binds choline in an aromatic cage and mediates its uptake. Our work could provide a novel framework for the targeted delivery of neurotherapeutics into the brain.

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.

TICBase: Integrated Resource for Data on Drug and Environmental Chemical Interactions with Mammalian Drug Transporters.

Environmental health science seeks to predict how environmental toxins, chemical toxicants, and prescription drugs accumulate and interact within the body. Xenobiotic transporters of the ATP-binding cassette (ABC) and solute carrier (SLC) superfamilies are major determinants of the uptake and disposition of xenobiotics across the kingdoms of life. The goal of this study was to integrate drug and environmental chemical interactions of mammalian ABC and SLC proteins in a centralized, integrative database. We built upon an existing publicly accessible platform-the "TransPortal"-which was updated with novel data and searchable features on transporter-interfering chemicals from manually curated literature data. The integrated resource TransPortal-TICBase (https://transportal.compbio.ucsf.edu) now contains information on 46 different mammalian xenobiotic transporters of the ABC- and SLC-type superfamilies, including 13 newly added rodent and 2 additional human drug transporters, 126 clinical drug-drug interactions, and a more than quadrupled expansion of the initial in vitro chemical interaction data from 1,402 to 6,296 total interactions. Based on our updated database, environmental interference with major human and rodent drug transporters occurs across the ABC- and SLC-type superfamilies, with kinetics indicating that some chemicals, such as the ionic liquid 1-hexylpyridinium chloride and the antiseptic chlorhexidine, can act as strong inhibitors with potencies similar or even higher than pharmacological model inhibitors. The new integrated web portal serves as a central repository of current and emerging data for interactions of prescription drugs and environmental chemicals with human drug transporters. This archive has important implications for predicting adverse drug-drug and drug-environmental chemical interactions and can serve as a reference website for the broader scientific community of clinicians and researchers.

The Impacts of Slc19a3 Deletion and Intestinal SLC19A3 Insertion on Thiamine Distribution and Brain Metabolism in the Mouse.

The Thiamine Transporter 2 (THTR2) encoded by SLC19A3 plays an ill-defined role in the maintenance of tissue thiamine, thiamine monophosphate, and thiamine diphosphate (TDP) levels. To evaluate the impact of THTR2 on tissue thiamine status and metabolism, we expressed the human SLC19A3 transgene in the intestine of total body Slc19a3 knockout (KO) mice. Male and female wildtype (WT) and transgenic (TG) mice were fed either 17 mg/kg (1×) or 85 mg/kg (5×) thiamine hydrochloride diet, while KOs were only fed the 5× diet. Thiamine vitamers in plasma, red blood cells, duodenum, brain, liver, kidney, heart, and adipose tissue were measured. Untargeted metabolomics were performed on the brain tissues of groups with equivalent plasma thiamine. KO mice had ~two- and ~three-fold lower plasma and brain thiamine levels than WT on the 5× diet. Circulating vitamers were sensitive to diet and equivalent in TG and WT mice. However, TG had 60% lower thiamine but normal brain TDP levels regardless of diet, with subtle differences in the heart and liver. The loss of THTR2 reduced levels of nucleic acid and amino acid derivatives in the brain. Therefore, mutation or inhibition of THTR2 may alter the brain metabolome and reduce the thiamine reservoir for TDP biosynthesis.

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.

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.

The Effect of Trimethoprim on Thiamine Absorption: A Transporter-Mediated Drug-Nutrient Interaction.

Trimethoprim is predicted to inhibit several thiamine transporters, including the primary thiamine intestinal absorptive transporter, ThTR-2, and the hepatic and renal organic cation transporters, OCT1, OCT2, and MATEs. To investigate the effect of trimethoprim on thiamine absorption, studies were conducted in cells, mice, and healthy volunteers and supported by use of real-world data. In a randomized, crossover clinical study, seven healthy volunteers were given a single oral dose of thiamine or thiamine plus trimethoprim, followed by blood sampling. The thiamine area under the curve (AUC) increased with trimethoprim co-administration (P value = 0.031). Similar results were seen in mice. Trimethoprim appeared to act on thiamine absorption through inhibition of hepatic OCT1 as evidenced from its ability to modulate levels of isobutyrylcarnitine and propionylcarnitine, OCT1 biomarkers identified from metabolomic analyses. Real-world data further supported this finding, showing an association between trimethoprim use and higher levels of triglycerides, LDL cholesterol, and total cholesterol, consistent with OCT1 inhibition (P values: 2.2 × 10-16 , 5.75 × 10-7 , and 5.82 × 10-7 , respectively). These findings suggest that trimethoprim increases plasma levels of thiamine by inhibiting hepatic OCT1. Trimethoprim reduced urinary excretion and clearance of biomarkers for OCT2 and MATEs, consistent with inhibition of renal organic cation transporters. This inhibition did not appear to play a role in the observed increases in thiamine levels. This study highlights the potential for drug-nutrient interactions involving transporters, in addition to transporters' established role in drug-drug interactions.

Characterization of P-glycoprotein orthologs from human, sheep, pig, dog, and cat.

The ATP-binding cassette transporter P-glycoprotein (P-gp) limits the oral bioavailability of many drugs. Although P-gp has been well studied in humans and mice, little is known about the substrate specificities of many of its species orthologs. To address this, we performed in vitro analysis of P-gp transporter function using HEK293 cells stably expressing human, ovine, porcine, canine, and feline P-gp. We also employed a human physiologically based pharmacokinetic (PBPK) model to assess variations in digoxin exposure resulting from altered P-gp function. Compared to human P-gp, sheep P-gp had significantly less digoxin efflux (2.3-fold ±0.04 vs. 1.8-fold ±0.03, p < .0001) and all species orthologs had significantly less quinidine efflux compared with human P-gp (p < .05). Human P-gp also had significantly greater efflux of talinolol compared to sheep and dog P-gp (1.9-fold ±0.04 vs. 1.6-fold ±0.06, p = .003 and 1.6-fold ±0.05, p = .0002, respectively). P-gp expression protected all lines against paclitaxel-induced toxicity, with sheep P-gp being significantly less protective. The inhibitor verapamil demonstrated dose-dependent inhibition of all P-gp orthologs. Finally, a PBPK model showed digoxin exposure was sensitive to altered P-gp activity. Overall, our study found that species differences in this major drug transporter exist and that the appropriate species ortholog of P-gp should be evaluated during veterinary drug development.

Identification of Genetic Variation Influencing Metformin Response in a Multiancestry Genome-Wide Association Study in the Diabetes Prevention Program (DPP).

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, ß = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, ß = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Initial Insights into the Genetic Variation Associated with Metformin Treatment Failure in Youth with Type 2 Diabetes.

Metformin is the first-line treatment for type 2 diabetes (T2D) in youth but with limited sustained glycemic response. To identify common variants associated with metformin response, we used a genome-wide approach in 506 youth from the Treatment Options for Type 2 Diabetes in Adolescents and Youth (TODAY) study and examined the relationship between T2D partitioned polygenic scores (pPS), glycemic traits, and metformin response in these youth. Several variants met a suggestive threshold (P < 1 × 10-6), though none including published adult variants reached genome-wide significance. We pursued replication of top nine variants in three cohorts, and rs76195229 in ATRNL1 was associated with worse metformin response in the Metformin Genetics Consortium (n = 7,812), though statistically not being significant after Bonferroni correction (P = 0.06). A higher ß-cell pPS was associated with a lower insulinogenic index (P = 0.02) and C-peptide (P = 0.047) at baseline and higher pPS related to two insulin resistance processes were associated with increased C-peptide at baseline (P = 0.04,0.02). Although pPS were not associated with changes in glycemic traits or metformin response, our results indicate a trend in the association of the ß-cell pPS with reduced ß-cell function over time. Our data show initial evidence for genetic variation associated with metformin response in youth with T2D.