In Vitro Functional Characterization and in Silico Prediction of Rare Genetic Variation in the Bile Acid and Drug Transporter, Na+-Taurocholate Cotransporting Polypeptide (NTCP, SLC10A1)

Transporter-mediated drug-drug interactions: regulatory guidelines, in vitro and in vivo methodologies and translation, special populations, and the blood-brain barrier

This review, part of a special issue on drug-drug interactions (DDIs) spearheaded by the International Society for the Study of Xenobiotics (ISSX) New Investigators, explores the critical role of drug transporters in absorption, disposition, and clearance in the context of DDIs. Over the past two decades, significant advances have been made in understanding the clinical relevance of these transporters. Current knowledge on key uptake and efflux transporters that affect drug disposition and development is summarized. Regulatory guidelines from the FDA, EMA, and PMDA that inform the evaluation of potential transporter-mediated DDIs are discussed in detail. Methodologies for preclinical and clinical testing to assess potential DDIs are reviewed, with an emphasis on the utility of physiologically based pharmacokinetic (PBPK) modeling. This includes the application of relative abundance and expression factors to predict human pharmacokinetics (PK) using preclinical data, integrating the latest regulatory guidelines. Considerations for assessing transporter-mediated DDIs in special populations, including pediatric, hepatic, and renal impairment groups, are provided. Additionally, the impact of transporters at the blood-brain barrier (BBB) on the disposition of CNS-related drugs is explored. Enhancing the understanding of drug transporters and their role in drug disposition and toxicity can improve efficacy and reduce adverse effects. Continued research is essential to bridge remaining gaps in knowledge, particularly in comparison with cytochrome P450 (CYP) enzymes.

Population genetic admixture and evolutionary history in the Shandong Peninsula inferred from integrative modern and ancient genomic resources

BACKGROUND

Ancient northern East Asians (ANEA) from the Yellow River region, who pioneered millet cultivation, play a crucial role in understanding the origins of ethnolinguistically diverse populations in modern China and the entire landscape of deep genetic structure and variation discovery in modern East Asians. However, the direct links between ANEA and geographically proximate modern populations, as well as the biological adaptive processes involved, remain poorly understood.

RESULTS

Here, we generated genome-wide SNP data for 264 individuals from geographically different Han populations in Shandong. An integrated genomic resource encompassing both modern and ancient East Asians was compiled to examine fine-scale population admixture scenarios and adaptive traits. The reconstruction of demographic history and hierarchical clustering patterns revealed that individuals from the Shandong Peninsula share a close genetic affinity with ANEA, indicating long-term genetic continuity and mobility in the lower Yellow River basin since the early Neolithic period. Biological adaptive signatures, including those related to immune and metabolic pathways, were identified through analyses of haplotype homozygosity and allele frequency spectra. These signatures are linked to complex traits such as height and body mass index, which may be associated with adaptations to cold environments, dietary practices, and pathogen exposure. Additionally, allele frequency trajectories over time and a haplotype network of two highly differentiated genes, ABCC11 and SLC10A1, were delineated. These genes, which are associated with axillary odor and bilirubin metabolism, respectively, illustrate how local adaptations can influence the diversification of traits in East Asians.

CONCLUSIONS

Our findings provide a comprehensive genomic dataset that elucidates the fine-scale genetic history and evolutionary trajectory of natural selection signals and disease susceptibility in Han Chinese populations. This study serves as a paradigm for integrating spatiotemporally diverse ancient genomes in the era of population genomic medicine.

Molecular mechanisms of Na+-driven bile acid transport in human NTCP

Human Na+ taurocholate co-transporting protein (hNTCP) is a key bile salt transporter to maintain enterohepatic circulation and is responsible for the recognition of hepatitis B and D viruses. Despite landmark cryoelectron microscopy studies revealing open-pore and inward-facing states of hNTCP stabilized by antibodies, the transport mechanism remains largely unknown. To address this knowledge gap, we used molecular dynamics and enhanced sampling metadynamics simulations to elucidate the intrinsic mechanism of hNTCP-mediated taurocholate acid (TCA) transport driven by Na+ binding. We uncovered three TCA-binding modes, including one that closely matched the limited cryoelectron microscopy density observed in the open-pore hNTCP. We also captured several key hNTCP conformations in the substrate transport cycle, particularly including an outward-facing, substrate-bound state. Furthermore, we provided thermodynamic evidence supporting that changes in the Na+-binding state drive the TCA transport by exploiting the amphiphilic nature of the substrate and modulating the protein environment, thereby enabling the TCA molecule to flip through. Understanding these mechanistic details of Na+-driven bile acid transport may aid in the development of hNTCP-targeted therapies for liver diseases.

Individualized Pharmacotherapy Utilizing Genetic Biomarkers and Novel In Vitro Systems As Predictive Tools for Optimal Drug Development and Treatment

In the area of drug development and clinical pharmacotherapy, a profound understanding of the pharmacokinetics and potential adverse reactions associated with the drug under investigation is paramount. Essential to this endeavor is a comprehensive understanding about interindividual variations in absorption, distribution, metabolism, and excretion (ADME) genetics and the predictive capabilities of in vitro systems, shedding light on metabolite formation and the risk of adverse drug reactions (ADRs). Both the domains of pharmacogenomics and the advancement of in vitro systems are experiencing rapid expansion. Here we present an update on these burgeoning fields, providing an overview of their current status and illuminating potential future directions. SIGNIFICANCE STATEMENT: There is very rapid development in the area of pharmacogenomics and in vitro systems for predicting drug pharmacokinetics and risk for adverse drug reactions. We provide an update of the current status of pharmacogenomics and developed in vitro systems on these aspects aimed to achieve a better personalized pharmacotherapy.

Translating pharmacogenomic sequencing data into drug response predictions-How to interpret variants of unknown significance

The rapid development of sequencing technologies during the past 20 years has provided a variety of methods and tools to interrogate human genomic variations at the population level. Pharmacogenes are well known to be highly polymorphic and a plethora of pharmacogenomic variants has been identified in population sequencing data. However, so far only a small number of these variants have been functionally characterized regarding their impact on drug efficacy and toxicity and the significance of the vast majority remains unknown. It is therefore of high importance to develop tools and frameworks to accurately infer the effects of pharmacogenomic variants and, eventually, aggregate the effect of individual variations into personalized drug response predictions. To address this challenge, we here first describe the technological advances, including sequencing methods and accompanying bioinformatic processing pipelines that have enabled reliable variant identification. Subsequently, we highlight advances in computational algorithms for pharmacogenomic variant interpretation and discuss the added value of emerging strategies, such as machine learning and the integrative use of omics techniques that have the potential to further contribute to the refinement of personalized pharmacological response predictions. Lastly, we provide an overview of experimental and clinical approaches to validate in silico predictions. We conclude that the iterative feedback between computational predictions and experimental validations is likely to rapidly improve the accuracy of pharmacogenomic prediction models, which might soon allow for an incorporation of the entire pharmacogenetic profile into personalized response predictions.

In vivo versus in silico assessment of potentially pathogenic missense variants in human reproductive genes

Infertility is a heterogeneous condition, with genetic causes thought to underlie a substantial fraction of cases. Genome sequencing is becoming increasingly important for genetic diagnosis of diseases including idiopathic infertility; however, most rare or minor alleles identified in patients are variants of uncertain significance (VUS). Interpreting the functional impacts of VUS is challenging but profoundly important for clinical management and genetic counseling. To determine the consequences of these variants in key fertility genes, we functionally evaluated 11 missense variants in the genes ANKRD31, BRDT, DMC1, EXO1, FKBP6, MCM9, M1AP, MEI1, MSH4 and SEPT12 by generating genome-edited mouse models. Nine variants were classified as deleterious by most functional prediction algorithms, and two disrupted a protein-protein interaction (PPI) in the yeast two hybrid (Y2H) assay. Though these genes are essential for normal meiosis or spermiogenesis in mice, only one variant, observed in the MCM9 gene of a male infertility patient, compromised fertility or gametogenesis in the mouse models. To explore the disconnect between predictions and outcomes, we compared pathogenicity calls of missense variants made by ten widely used algorithms to 1) those annotated in ClinVar and 2) those evaluated in mice. All the algorithms performed poorly in terms of predicting the effects of human missense variants modeled in mice. These studies emphasize caution in the genetic diagnoses of infertile patients based primarily on pathogenicity prediction algorithms and emphasize the need for alternative and efficient in vitro or in vivo functional validation models for more effective and accurate VUS description to either pathogenic or benign categories.

Role of Hepatocyte Transporters in Drug-Induced Liver Injury (DILI)-In Vitro Testing

Bile acids and bile salts (BA/BS) are substrates of both influx and efflux transporters on hepatocytes. Canalicular efflux transporters, such as BSEP and MRP2, are crucial for the removal of BA/BS to the bile. Basolateral influx transporters, such as NTCP, OATP1B1/1B3, and OSTα/β, cooperate with canalicular transporters in the transcellular vectorial flux of BA/BS from the sinusoids to the bile. The blockage of canalicular transporters not only impairs the bile flow but also causes the intracellular accumulation of BA/BS in hepatocytes that contributes to, or even triggers, liver injury. In the case of BA/BS overload, the efflux of these toxic substances back to the blood via MRP3, MRP4, and OST α/β is considered a relief function. FXR, a key regulator of defense against BA/BS toxicity suppresses de novo bile acid synthesis and bile acid uptake, and promotes bile acid removal via increased efflux. In drug development, the early testing of the inhibition of these transporters, BSEP in particular, is important to flag compounds that could potentially inflict drug-induced liver injury (DILI). In vitro test systems for efflux transporters employ membrane vesicles, whereas those for influx transporters employ whole cells. Additional in vitro pharmaceutical testing panels usually include cellular toxicity tests using hepatocytes, as well as assessments of the mitochondrial toxicity and accumulation of reactive oxygen species (ROS). Primary hepatocytes are the cells of choice for toxicity testing, with HepaRG cells emerging as an alternative. Inhibition of the FXR function is also included in some testing panels. The molecular weight and hydrophobicity of the drug, as well as the steady-state total plasma levels, may positively correlate with the DILI potential. Depending on the phase of drug development, the physicochemical properties, dosing, and cut-off values of BSEP IC₅₀ ≤ 25-50 µM or total C_ss,plasma/BSEP IC₅₀ ≥ 0.1 may be an indication for further testing to minimize the risk of DILI liability.

A paradigm shift in pharmacogenomics: From candidate polymorphisms to comprehensive sequencing

Genetic factors have long been recognized as important determinants of interindividual variability in drug efficacy and toxicity. However, despite the increasing number of established gene-drug associations, candidate polymorphisms can only explain a fraction of the genetically encoded functional variability in drug disposition. Advancements in genetic profiling methods now allow to analyse the landscape of human pharmacogenetic variations comprehensively, which opens new opportunities to identify novel factors that could explain the "missing heritability." Here, we provide an updated overview of the landscape of pharmacogenomic variability based on recent analyses of population-scale sequencing projects. We then summarize the current state-of-the-art how the functional consequences of variants with unknown effects can be quantitatively estimated while discussing challenges and peculiarities that are specific to pharmacogenes. In the last sections, we discuss the importance of considering ethnogeographic diversity to provide equitable benefits of pharmacogenomics and summarize current roadblocks for the implementation of sequencing-based guidance of clinical decision-making. Based on the current state of the field, we conclude that testing is likely to gradually shift from the interrogation of selected candidate polymorphisms to comprehensive sequencing, which allows to consider the full spectrum of pharmacogenomic variations for a true personalization of genomic prescribing.

Hepatic Expression of the Na+-Taurocholate Cotransporting Polypeptide Is Independent from Genetic Variation

The hepatic Na+-taurocholate cotransporting polypeptide NTCP/SLC10A1 is important for the uptake of bile salts and selected drugs. Its inhibition results in increased systemic bile salt concentrations. NTCP is also the entry receptor for the hepatitis B/D virus. We investigated interindividual hepatic SLC10A1/NTCP expression using various omics technologies. SLC10A1/NTCP mRNA expression/protein abundance was quantified in well-characterized 143 human livers by real-time PCR and LC-MS/MS-based targeted proteomics. Genome-wide SNP arrays and SLC10A1 next-generation sequencing were used for genomic analyses. SLC10A1 DNA methylation was assessed through MALDI-TOF MS. Transcriptomics and untargeted metabolomics (UHPLC-Q-TOF-MS) were correlated to identify NTCP-related metabolic pathways. SLC10A1 mRNA and NTCP protein levels varied 44-fold and 10.4-fold, respectively. Non-genetic factors (e.g., smoking, alcohol consumption) influenced significantly NTCP expression. Genetic variants in SLC10A1 or other genes do not explain expression variability which was validated in livers (n = 50) from The Cancer Genome Atlas. The identified two missense SLC10A1 variants did not impair transport function in transfectants. Specific CpG sites in SLC10A1 as well as single metabolic alterations and pathways (e.g., peroxisomal and bile acid synthesis) were significantly associated with expression. Inter-individual variability of NTCP expression is multifactorial with the contribution of clinical factors, DNA methylation, transcriptional regulation as well as hepatic metabolism, but not genetic variation.

In-vitro characterization of coding variants with predicted functional implications in the efflux transporter multidrug resistance protein 4 (MRP4, ABCC4)

MRP4 (gene ABCC4) is a polymorphic efflux transporter that has been implicated in drug-induced toxicity. We selected ten commonly observed MRP4 coding variants among Europeans for experimental characterization including nine variants predicted to be deleterious or functional (combined annotation-dependent depletion score >15). We assessed protein localization and activity by quantifying intracellular accumulation of two prototypic substrates, taurocholic acid (TCA) and estradiol 17-β-glucuronide (E217βG), in HEK293T over-expressing MRP4 wildtype or variant where cellular substrate loading was optimized through co-transfection with an uptake transporter. V458M, a novel variant not previously studied, and T1142M, showed reduced activity compared to MRP4 wildtype for E217βG and TCA (P < 0.01), while L18I, G187W, K293E, and R531Q moderately increased activity in a substrate-dependent manner. Protein expression analysis indicated reduced cell surface expression for V458M (P < 0.01) but not T1142M compared to wildtype. Reduced activity may result from altered surface expression (V458M) or intrinsic activity as both variants map within the nucleotide-binding domains of MRP4. G187W showed a trend for reduced surface expression (P = 0.054) despite transport comparable or increased to wildtype suggesting enhanced intrinsic activity. Our findings suggest moderately altered MRP4 activity in six out of nine predicted functional variants with likely different mechanisms and substrate-specific effects. Cell-based studies using multiple known substrates are warranted to more accurately predict functional variants in this clinically important transporter.

Structural Plasticity Is a Feature of Rheostat Positions in the Human Na+/Taurocholate Cotransporting Polypeptide (NTCP)

In the Na⁺/taurocholate cotransporting polypeptide (NTCP), the clinically relevant S267F polymorphism occurs at a "rheostat position". That is, amino acid substitutions at this position ("S267X") lead to a wide range of functional outcomes. This result was particularly striking because molecular models predicted the S267X side chains are buried, and thus, usually expected to be less tolerant of substitutions. To assess whether structural tolerance to buried substitutions is widespread in NTCP, here we used Rosetta to model all 19 potential substitutions at another 13 buried positions. Again, only subtle changes in the calculated stabilities and structures were predicted. Calculations were experimentally validated for 19 variants at codon 271 ("N271X"). Results showed near wildtype expression and rheostatic modulation of substrate transport, implicating N271 as a rheostat position. Notably, each N271X substitution showed a similar effect on the transport of three different substrates and thus did not alter substrate specificity. This differs from S267X, which altered both transport kinetics and specificity. As both transport and specificity may change during protein evolution, the recognition of such rheostat positions may be important for evolutionary studies. We further propose that the presence of rheostat positions is facilitated by local plasticity within the protein structure. Finally, we note that identifying rheostat positions may advance efforts to predict new biomedically relevant missense variants in NTCP and other membrane transport proteins.

Pharmacogenomics of statins: lipid response and other outcomes in Brazilian cohorts

Statins are inhibitors of 3-hydroxy-3-methylglutaryl-CoA reductase, a key enzyme in cholesterol biosynthesis, that are highly effective in reducing plasma low-density lipoprotein (LDL) cholesterol and decreasing the risk of cardiovascular events. In recent years, a multitude of variants in genes involved in pharmacokinetics (PK) and pharmacodynamics (PD) have been suggested to influence the cholesterol-lowering response. However, the vast majority of studies have analyzed the pharmacogenetic associations in populations in Europe and the USA, whereas data in other populations, including Brazil, are mostly lacking. This narrative review provides an update of clinical studies on statin pharmacogenomics in Brazilian cohorts exploring lipid-lowering response, adverse events and pleiotropic effects. We find that variants in drug transporter genes (SLCO1B1 and ABCB1) positively impacted atorvastatin and simvastatin response, whereas variants in genes of drug metabolizing enzymes (CYP3A5) decreased response. Furthermore, multiple associations of variants in PD genes (HMGCR, LDLR and APOB) with statin response were identified. Few studies have explored statin-related adverse events, and only ABCB1 but not SLCO1B1 variants were robustly associated with increased risk in Brazil. Statin-related pleiotropic effects were shown to be influenced by variants in PD (LDLR, NR1H2) and antioxidant enzyme (NOS3, SOD2, MTHFR, SELENOP) genes. The findings of these studies indicate that statin pharmacogenomic associations are distinctly different in Brazil compared to other populations. This review also discusses the clinical implications of pharmacogenetic studies and the rising importance of investigating rare variants to explore their association with statin response.

Bile formation and secretion: An update

Bile formation is a fundamental physiological process that is vital to the survival of all vertebrates. However, little was known about the mechanisms of this secretion until after World War II. Initial studies involved classic physiologic studies in animal models and humans, which progressed to include studies in isolated cells and membrane vesicles. The advent of molecular biology then led to the identification of specific transport systems that are the determinants of this secretion. Progress in this field was reviewed in the American Physiologic Society's series on "Comprehensive Physiology" in 2013. Herein, we provide an in-depth update of progress since that time.

Pharmacogenomics in the era of next generation sequencing - from byte to bedside

Pharmacogenetic research has resulted in the identification of a multitude of genetic variants that impact drug response or toxicity. These polymorphisms are mostly common and have been included as actionable information in the labels of numerous drugs. In addition to common variants, recent advances in Next Generation Sequencing (NGS) technologies have resulted in the identification of a plethora of rare and population-specific pharmacogenetic variations with unclear functional consequences that are not accessible by conventional forward genetics strategies. In this review, we discuss how comprehensive sequencing information can be translated into personalized pharmacogenomic advice in the age of NGS. Specifically, we provide an update of the functional impacts of rare pharmacogenetic variability and how this information can be leveraged to improve pharmacogenetic guidance. Furthermore, we critically discuss the current status of implementation of pharmacogenetic testing across drug development and layers of care. We identify major gaps and provide perspectives on how these can be minimized to optimize the utilization of NGS data for personalized clinical decision-support.

Hyperosmolarity stimulates transporter-mediated insertion of estrone sulfate into the plasma membrane, but inhibits the uptake by SLC10A1 (NTCP)

Many drugs are largely hydrophobic molecules; a transporter might conceivably insert these into the plasma membrane. At least 18 transporters from diverse families have been reported to transport the model compound estrone sulfate alias estrone-3-sulfate (E3S). Out of these, we recently examined SLC22A11 (OAT4). We concluded from a comparison of E3S and uric acid transport that SLC22A11 does not translocate E3S into the cytosol, but into the plasma membrane. Here we present a hyperosmolarity alias hypertonicity assay to differentiate transport mechanisms. Human transporters were expressed heterologously in 293 cells. Solute uptake into intact cells was measured by LC-MS. Addition of mannitol or sucrose led to rapid cell shrinkage, but cell viability after 60 min in hyperosmolar buffer was not impaired. A decrease in substrate accumulation with increasing osmolarity as observed here for several substrates and the transporters SLC22A11, ETT (SLC22A4), OCT2 (SLC22A2), OAT3 (SLC22A8), and MATE1 (SLC47A1) suggests regular substrate translocation into the cytosol. An increase as observed for E3S transport by SLC22A11, OAT3, MATE1, SLC22A9, and SLC10A6 implies insertion into the membrane. In marked contrast to the other E3S transporters, the bile acid transporter SLC10A1 (NTCP, Na⁺ taurocholate co-transporting polypeptide) showed a decrease in the accumulation of E3S in hyperosmolar buffer; the same was observed with taurocholic acid. Indeed, our data from several functional assays strongly suggest that the transport mechanism is identical for both substrates. Apparently, a unique transport mechanism has been established for SLC10A1 by evolution that ensures the transport of amphipathic, detergent-like molecules into the cytosol.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Variant discovery using next-generation sequencing and its future role in pharmacogenetics

Next-generation sequencing (NGS) has enabled the discovery of a multitude of novel and mostly rare variants in pharmacogenes that may alter a patient’s therapeutic response to drugs. In addition to single nucleotide variants, structural variation affecting the number of copies of whole genes or parts of genes can be detected. While current guidelines concerning clinical implementation mostly act upon well-documented, common single nucleotide variants to guide dosing or drug selection, in silico and large-scale functional assessment of rare variant effects on protein function are at the forefront of pharmacogenetic research to facilitate their clinical integration. Here, we discuss the role of NGS in variant discovery, paving the way for more comprehensive genotype-guided pharmacotherapy that can translate to improved clinical care.