[Fatigue and sensorimotor instability : Neurologically controlled conversion of post-COVID-19 patients]. / Fatigue und sensomotorische Instabilität : Neurologisch kontrollierte Konversion von Post-COVID-19-Patienten.

BACKGROUND: For the treatment of the symptoms of post-COVID-19 syndrome, no causal therapy is currently widely recommended according to evidence-based criteria. The overarching goal of the intervention study over a 3-year period (Q1-2021-Q4-2023) was to evaluate the changes in the key symptoms of fatigue and sensorimotor instability through individualized stress-controlled training therapy and through intensified cognitive behavioral therapy. MATERIAL AND METHODS: In the 3­year period 407 vaccinated nucleocapsid positive patients were treated at the Post-COVID-19 Center Lausitz (Senftenberg). In 78 (around 19%) fatigue/immunometabolic depression and sensorimotor instability were identified as the leading syndromes. The evaluation of the individualized stress-controlled training therapy was based on the specific post-COVID-19 syndrome and motor fatigue parameters. The secondary psychosomatic syndrome was assessed using cognitive fatigue parameters and cognitive behavioral therapy instruments. The investigation of -parameters influencing behavior took place in Q2-2023-Q4-2023 with a guide-supported qualitative interview among the participants. RESULTS: The post-COVID-19 key symptoms "fatigue," "sensorimotor instability," "neuropsychiatric symptoms," "cardiac/autonomic dysfunction," and "pain" improved significantly in the overall cohort and in the gender-specific analysis. A deterioration occurred in "secondary psychosomatic symptoms". A therapeutic effect was demonstrated for all motor fatigue parameters for the entire cohort using the Cohen's d value. An intensification of cognitive behavioral therapy achieved positive effects through an increasing development of the patients' own activity and their self-control using persuasion and gamification.

A program for real-time surveillance of SARS-CoV-2 genetics.

The COVID-19 pandemic brought forth an urgent need for widespread genomic surveillance for rapid detection and monitoring of emerging SARS-CoV-2 variants. It necessitated design, development, and deployment of a nationwide infrastructure designed for sequestration, consolidation, and characterization of patient samples that disseminates de-identified information to public authorities in tight turnaround times. Here, we describe our development of such an infrastructure, which sequenced 594,832 high coverage SARS-CoV-2 genomes from isolates we collected in the United States (U.S.) from March 13th 2020 to July 3rd 2023. Our sequencing protocol ('Virseq') utilizes wet and dry lab procedures to generate mutation-resistant sequencing of the entire SARS-CoV-2 genome, capturing all major lineages. We also characterize 379 clinically relevant SARS-CoV-2 multi-strain co-infections and ensure robust detection of emerging lineages via simulation. The modular infrastructure, sequencing, and analysis capabilities we describe support the U.S. Centers for Disease Control and Prevention national surveillance program and serve as a model for rapid response to emerging pandemics at a national scale.

A genome-wide screen for variants influencing certolizumab pegol response in a moderate to severe rheumatoid arthritis population.

Certolizumab pegol (CZP) is a PEGylated Fc-free tumor necrosis factor (TNF) inhibitor antibody approved for use in the treatment of rheumatoid arthritis (RA), Crohn's disease, psoriatic arthritis, axial spondyloarthritis and psoriasis. In a clinical trial of patients with severe RA, CZP improved disease symptoms in approximately half of patients. However, variability in CZP efficacy remains a problem for clinicians, thus, the aim of this study was to identify genetic variants predictive of CZP response. We performed a genome-wide association study (GWAS) of 302 RA patients treated with CZP in the REALISTIC trial to identify common single nucleotide polymorphisms (SNPs) associated with treatment response. Whole-exome sequencing was also performed for 74 CZP extreme responders and non-responders within the same population, as well as 1546 population controls. No common SNPs or rare functional variants were significantly associated with CZP response, though a non-significant enrichment in the RA-implicated KCNK5 gene was observed. Two SNPs near spondin-1 and semaphorin-4G approached genome-wide significance. The results of the current study did not provide an unambiguous predictor of CZP response.

Response to Comment on "Evidence of humans in North America during the Last Glacial Maximum".

Madsen et al. question the reliability of calibrated radiocarbon ages associated with human footprints discovered recently in White Sands National Park, New Mexico, USA. On the basis of the geologic, hydrologic, stratigraphic, and chronologic evidence, we maintain that the ages are robust and conclude that the footprints date to between ~23,000 and 21,000 years ago.Madsen et al. (1) question the veracity of calibrated radiocarbon ages used to constrain the antiquity of human trackways discovered recently at White Sands National Park (WHSA) Locality 2, New Mexico, USA (2). The ages were derived from seeds of the aquatic plant Ruppia cirrhosa, which they suggest may suffer from hard-water (or reservoir) effects, making them too old, potentially by thousands of years. We were well aware of this possibility, investigated it, and presented several lines of evidence that argued against such a problem. Here we respond to each of their four primary points.

Evidence of humans in North America during the Last Glacial Maximum.

Archaeologists and researchers in allied fields have long sought to understand human colonization of North America. Questions remain about when and how people migrated, where they originated, and how their arrival affected the established fauna and landscape. Here, we present evidence from excavated surfaces in White Sands National Park (New Mexico, United States), where multiple in situ human footprints are stratigraphically constrained and bracketed by seed layers that yield calibrated radiocarbon ages between ~23 and 21 thousand years ago. These findings confirm the presence of humans in North America during the Last Glacial Maximum, adding evidence to the antiquity of human colonization of the Americas and providing a temporal range extension for the coexistence of early inhabitants and Pleistocene megafauna.

3-D radar imaging unlocks the untapped behavioral and biomechanical archive of Pleistocene ghost tracks.

Footprint evidence of human-megafauna interactions remains extremely rare in the archaeological and palaeontological records. Recent work suggests ancient playa environments may hold such evidence, though the prints may not be visible. These so-called "ghost tracks" comprise a rich archive of biomechanical and behavioral data that remains mostly unexplored. Here we present evidence for the successful detection and 3-D imaging of such footprints via ground-penetrating radar (GPR), including co-associated mammoth and human prints. Using GPR we have found that track density and faunal diversity may be much greater than realized by the unaided human eye. Our data further suggests that detectable subsurface consolidation below mammoth tracks correlates with typical plantar pressure patterns from extant elephants. This opens future potential for more sophisticated biomechanical studies on the footprints of other extinct land vertebrates. Our approach allows rapid detection and documentation of footprints while enhancing the data available from these fossil archives.

An exploratory genome-wide analysis of genetic risk for alcoholic hepatitis.

OBJECTIVES: To elucidate the genetic variability between heavy drinkers with and without alcoholic hepatitis (AH). MATERIALS AND METHODS: An exploratory genome-wide association study (GWAS; NCT02172898) was conducted comparing 90 AH cases with 93 heavy drinking matched controls without liver disease in order to identify variants or genes associated with risk for AH. Individuals were genotyped using the multi-ethnic genotyping array, after which the data underwent conventional quality control. Using bioinformatics tools, pathways associated with AH were explored on the basis of individual variants, and based on genes with a higher 'burden' of functional variation. RESULTS: Although no single variant reached genome-wide significance, an association signal was observed for PNPLA3 rs738409 (p = .01, OR 1.9, 95% CI 1.1-3.1), a common single nucleotide polymorphism that has been associated with a variety of liver-related pathologies including alcoholic cirrhosis. Using the improved gene set enrichment analysis for GWAS tool, it was shown that, based on the single variants' trait-association p-values, multiple pathways were associated with risk for AH with high confidence (false discovery rate [FDR] < 0.05), including several pathways involved in lymphocyte activation and chemokine signaling, which coincides with findings from other research groups. Several Tox Functions and Canonical Pathways were highlighted using Ingenuity Pathway Analysis, with an especially conspicuous role for pathways related to ethanol degradation, which is not surprising considering the phenotype of the genotyped individuals. CONCLUSION: This preliminary analysis suggests a role for PNPLA3 variation and several gene sets/pathways that may influence risk for AH among heavy drinkers.

Minocycline hepatotoxicity: Clinical characterization and identification of HLA-B∗35:02 as a risk factor.

BACKGROUND & AIMS: Minocycline hepatotoxicity can present with prominent autoimmune features in previously healthy individuals. The aim of this study was to identify genetic determinants of minocycline drug-induced liver injury (DILI) in a well-phenotyped cohort of patients. METHODS: Caucasian patients with minocycline DILI underwent genome-wide genotyping and were compared to unexposed population controls. Human leukocyte antigen (HLA) binding of minocycline was assessed using AutoDock Vina. RESULTS: Among the 25 cases, 80% were female, median age was 19years and median latency from drug start to DILI onset was 318days. At presentation, 76% had acute hepatocellular liver injury, median ALT 1,077U/L (range: 63 to 2,333), median bilirubin 4.5mg/dl (range: 0.2 to 16.7), and 90% had a +ANA. During follow-up, 50% were treated with corticosteroids and no participants died or required a liver transplant. A significant association was noted between HLA-B∗35:02 and risk for minocycline DILI; a 16% carrier frequency in DILI cases compared to 0.6% in population controls (odds ratio: 29.6, 95% CI: 7.8-89.8, p=2.5×10-8). Verification of HLA-B∗35:02 imputation was confirmed by sequence-based HLA typing. HLA-B∗35:02 carriers had similar presenting features and outcomes compared to non-carriers. In silico modeling studies support the hypothesis that direct binding of minocycline to this novel HLA risk allele might be an important initiating event in minocycline DILI. CONCLUSION: HLA-B∗35:02 is a rare HLA allele that was more frequently identified in the 25 minocycline DILI cases compared to population controls. If confirmed in other cohorts, this HLA allele may prove to be a useful diagnostic marker of minocycline DILI. LAY SUMMARY: Development of liver injury following prolonged use of minocycline for acne is a rare but potentially severe form of drug-induced liver injury. Our study demonstrates that individuals who are HLA-B∗35:02 carriers are at increased risk of developing minocycline related liver injury. These results may help doctors more rapidly and confidently diagnose affected patients and possibly reduce the risk of liver injury in individuals receiving minocycline going forward.

Association of Liver Injury From Specific Drugs, or Groups of Drugs, With Polymorphisms in HLA and Other Genes in a Genome-Wide Association Study.

BACKGROUND & AIMS: We performed a genome-wide association study (GWAS) to identify genetic risk factors for drug-induced liver injury (DILI) from licensed drugs without previously reported genetic risk factors. METHODS: We performed a GWAS of 862 persons with DILI and 10,588 population-matched controls. The first set of cases was recruited before May 2009 in Europe (n = 137) and the United States (n = 274). The second set of cases were identified from May 2009 through May 2013 from international collaborative studies performed in Europe, the United States, and South America. For the GWAS, we included only cases with patients of European ancestry associated with a particular drug (but not flucloxacillin or amoxicillin-clavulanate). We used DNA samples from all subjects to analyze HLA genes and single nucleotide polymorphisms. After the discovery analysis was concluded, we validated our findings using data from 283 European patients with diagnosis of DILI associated with various drugs. RESULTS: We associated DILI with rs114577328 (a proxy for A*33:01 a HLA class I allele; odds ratio [OR], 2.7; 95% confidence interval [CI], 1.9-3.8; P = 2.4 × 10-8) and with rs72631567 on chromosome 2 (OR, 2.0; 95% CI, 1.6-2.5; P = 9.7 × 10-9). The association with A*33:01 was mediated by large effects for terbinafine-, fenofibrate-, and ticlopidine-related DILI. The variant on chromosome 2 was associated with DILI from a variety of drugs. Further phenotypic analysis indicated that the association between DILI and A*33:01 was significant genome wide for cholestatic and mixed DILI, but not for hepatocellular DILI; the polymorphism on chromosome 2 was associated with cholestatic and mixed DILI as well as hepatocellular DILI. We identified an association between rs28521457 (within the lipopolysaccharide-responsive vesicle trafficking, beach and anchor containing gene) and only hepatocellular DILI (OR, 2.1; 95% CI, 1.6-2.7; P = 4.8 × 10-9). We did not associate any specific drug classes with genetic polymorphisms, except for statin-associated DILI, which was associated with rs116561224 on chromosome 18 (OR, 5.4; 95% CI, 3.0-9.5; P = 7.1 × 10-9). We validated the association between A*33:01 terbinafine- and sertraline-induced DILI. We could not validate the association between DILI and rs72631567, rs28521457, or rs116561224. CONCLUSIONS: In a GWAS of persons of European descent with DILI, we associated HLA-A*33:01 with DILI due to terbinafine and possibly fenofibrate and ticlopidine. We identified polymorphisms that appear to be associated with DILI from statins, as well as 2 non-drug-specific risk factors.

An Exome Sequencing Study to Assess the Role of Rare Genetic Variation in Pulmonary Fibrosis.

RATIONALE: Idiopathic pulmonary fibrosis (IPF) is an increasingly recognized, often fatal lung disease of unknown etiology. OBJECTIVES: The aim of this study was to use whole-exome sequencing to improve understanding of the genetic architecture of pulmonary fibrosis. METHODS: We performed a case-control exome-wide collapsing analysis including 262 unrelated individuals with pulmonary fibrosis clinically classified as IPF according to American Thoracic Society/European Respiratory Society/Japanese Respiratory Society/Latin American Thoracic Association guidelines (81.3%), usual interstitial pneumonia secondary to autoimmune conditions (11.5%), or fibrosing nonspecific interstitial pneumonia (7.2%). The majority (87%) of case subjects reported no family history of pulmonary fibrosis. MEASUREMENTS AND MAIN RESULTS: We searched 18,668 protein-coding genes for an excess of rare deleterious genetic variation using whole-exome sequence data from 262 case subjects with pulmonary fibrosis and 4,141 control subjects drawn from among a set of individuals of European ancestry. Comparing genetic variation across 18,668 protein-coding genes, we found a study-wide significant (P < 4.5 × 10-7) case enrichment of qualifying variants in TERT, RTEL1, and PARN. A model qualifying ultrarare, deleterious, nonsynonymous variants implicated TERT and RTEL1, and a model specifically qualifying loss-of-function variants implicated RTEL1 and PARN. A subanalysis of 186 case subjects with sporadic IPF confirmed TERT, RTEL1, and PARN as study-wide significant contributors to sporadic IPF. Collectively, 11.3% of case subjects with sporadic IPF carried a qualifying variant in one of these three genes compared with the 0.3% carrier rate observed among control subjects (odds ratio, 47.7; 95% confidence interval, 21.5-111.6; P = 5.5 × 10-22). CONCLUSIONS: We identified TERT, RTEL1, and PARN-three telomere-related genes previously implicated in familial pulmonary fibrosis-as significant contributors to sporadic IPF. These results support the idea that telomere dysfunction is involved in IPF pathogenesis.

Variant Inosine Triphosphatase Phenotypes Are Associated With Increased Ribavirin Triphosphate Levels.

Individuals with lower inosine triphosphatase (ITPA) enzyme activity have a reduced likelihood of experiencing hemolytic anemia during hepatitis C virus (HCV) treatment containing ribavirin (RBV). Because ITPA degrades purines and RBV is a purine analogue, it is conceivable that ITPA activity may affect intracellular RBV concentrations. Here we assessed the association between ITPA activity phenotype and concentrations of RBV triphosphate (RBV-TP) in red blood cells (RBCs) during HCV treatment. RBV-TP was quantified in the RBCs of 177 HCV-infected individuals at a median (range) of 84 (19 to 336) days into HCV treatment that included RBV. Mean (SD) RBV-TP concentrations were 92.8 (51.6), 101.3 (53.5), 184.8 (84.5), and 197.7 (64.6) pmol/106 cells for 100%, 60%, 30%, and ≤10% ITPA activity groups, respectively. Overall, RBV-TP was approximately 2-fold higher in patients with ≤30% ITPA activity compared to 100% activity (P < .0001). Despite higher RBV-TP levels, individuals with variant ITPA phenotypes had less anemia. The 100% activity group had, on average, a -2.20 g/dL drop in hemoglobin vs -1.43 g/dL (P = .04) for 60% activity, -1.14 g/dL (P = .008) for 30% activity, and -0.70 g/dL (P = .06) for ≤10% activity. This finding of higher RBV-TP concentrations in RBCs in ITPA variants was unexpected given that ITPA activity-deficient individuals have a reduced likelihood of RBV-induced anemia. It also refutes the hypothesis that the mechanism by which ITPA variants are protected against anemia is due to lower RBV-TP levels in RBCs.

Features of Autoimmune Hepatitis in Patients With Drug-induced Liver Injury.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) has features similar to those of other liver diseases including autoimmune hepatitis (AIH). We aimed to characterize the clinical and autoimmune features of liver injury caused by nitrofurantoin, minocycline, methyldopa, or hydralazine. METHODS: We analyzed data from 88 cases of DILI attributed to nitrofurantoin, minocycline, methyldopa, or hydralazine included in the Drug-Induced Liver Injury Network prospective study from 2004 through 2014. Sera were collected from patients at baseline and follow-up examination and tested for levels of immunoglobulin G (IgG), antibodies to nuclear antigen (ANA), smooth muscle (SMA), and soluble liver antigen (SLA). An autoimmune score was derived on the basis of increases in levels of IgG, ANA, SMA, and SLA (assigned values of 0, 1+, or 2+). AIH-associated HLA-DRB1*03:01 and HLA-DRB1*04:01 allele frequencies were compared with those of the general population (controls). RESULTS: Of the 88 cases, 80 were women (91%), 74% had hepatocellular injury, and 25% had severe injury. At the onset of DILI, 39% of cases had increased levels of IgG, 72% had increased levels of ANA, 60% had increased levels of SMA, and none had increases in SLA. A phenotype of autoimmunity (autoimmune score ≥2) was observed in 82% of cases attributed to nitrofurantoin and 73% of cases attributed to minocycline (73%) but only 55% of cases attributed to methyldopa and 43% of cases attributed to hydralazine (P = .16 for nitrofurantoin and minocycline vs methyldopa and hydralazine). We observed a decrease in numbers of serum samples positive for ANA (P = .01) or SMA (P < .001) and in autoimmune scores (P < .001) between DILI onset and follow-up. Similar percentages of patients with DILI had HLA-DRB1*03:01 (15%) and HLA-DRB1*04:01 (9%) as controls (12% and 9%, respectively). CONCLUSIONS: In analysis of data from the DILIN prospective study, we found that most cases of DILI attributed to nitrofurantoin or minocycline and about half of cases that were due to methyldopa and hydralazine have a phenotype of autoimmunity similar to AIH. These features decrease with recovery of the injury and are not associated with the typical HLA alleles found in patients with idiopathic AIH.

Effects of Age, Sex, Body Weight, and Quantity of Alcohol Consumption on Occurrence and Severity of Alcoholic Hepatitis.

BACKGROUND & AIMS: Only a minority of heavy drinking individuals develop alcoholic hepatitis (AH), for unclear reasons. We analyzed data from the Translational Research and Evolving Alcoholic Hepatitis Treatment cohort, consisting of subjects who drink heavily with normal results from liver tests (controls) and patients with AH. We examined risk factors for the development of AH including body mass index (BMI), drinking pattern and quantity, and sex. METHODS: We compared data from 145 patients with AH and 124 controls based on BMI when they joined the cohort; groups were matched for sex and race. Drinking patterns were assessed using the timeline followback method, the Alcohol Use Disorders Identification Test, and the National Institute of Alcohol Abuse and Alcoholism 6-question survey. We performed univariable and multivariable analyses to assess the effects of these factors and their interaction in increasing the risk for AH. We also explored the association between PNPLA3 variants and AH. RESULTS: Cases with AH were older (47 vs 44 y; P = .03). For nearly all measures of quantity of alcohol consumed or frequency of binge drinking, controls drank more heavily than cases with AH. We did not find an association between BMI, sex, drinking patterns, and the presence of AH. Age and BMI were independent predictors for the severity of AH. When we analyzed cases and controls of European ancestry, the PNPLA3 single-nucleotide polymorphism rs738409 was associated with risk for AH (odds ratio, 1.89; P = .007). CONCLUSIONS: Compared with heavy drinkers without liver disease, subjects with AH consumed lower levels of alcohol and had less binge drinking, suggesting an increased sensitivity to the toxic effects of alcohol. The risk for AH may be associated with the PNPLA3 rs738409 polymorphism.

HLA-DRB1*16: 01-DQB1*05: 02 is a novel genetic risk factor for flupirtine-induced liver injury.

OBJECTIVE: Flupirtine is a nonopioid analgesic with regulatory approval in a number of European countries. Because of the risk of serious liver injury, its use is now limited to short-term pain management. We aimed to identify genetic risk factors for flupirtine-related drug-induced liver injury (DILI) as these are unknown. MATERIALS AND METHODS: Six flupirtine-related DILI patients from Germany were included in a genome-wide association study (GWAS) involving a further 614 European cases of DILI because of other drugs and 10,588 population controls. DILI was diagnosed by causality assessment and expert review. Human leucocyte antigen (HLA) and single nucleotide polymorphism genotypes were imputed from the GWAS data, with direct HLA typing performed on selected cases to validate HLA predictions. Four replication cases that were unavailable for the GWAS were genotyped by direct HLA typing, yielding an overall total of 10 flupirtine DILI cases. RESULTS: In the six flupirtine DILI cases included in the GWAS, we found a significant enrichment of the DRB1*16:01-DQB1*05:02 haplotype compared with the controls (minor allele frequency cases 0.25 and minor allele frequency controls 0.013; P=1.4 × 10(-5)). We estimated an odds ratio for haplotype carriers of 18.7 (95% confidence interval 2.5-140.5, P=0.002) using population-specific HLA control data. The result was replicated in four additional cases, also with a haplotype frequency of 0.25. In the combined cohort (six GWAS plus four replication cases), the haplotype was also significant (odds ratio 18.7, 95% confidence interval 4.31-81.42, P=6.7 × 10(-5)). CONCLUSION: We identified a novel HLA class II association for DILI, confirming the important contribution of HLA genotype towards the risk of DILI generally.

IFNL3 mRNA structure is remodeled by a functional non-coding polymorphism associated with hepatitis C virus clearance.

Polymorphisms near the interferon lambda 3 (IFNL3) gene strongly predict clearance of hepatitis C virus (HCV) infection. We analyzed a variant (rs4803217 G/T) located within the IFNL3 mRNA 3' untranslated region (UTR); the G allele (protective allele) is associated with elevated therapeutic HCV clearance. We show that the IFNL3 3' UTR represses mRNA translation and the rs4803217 allele modulates the extent of translational regulation. We analyzed the structures of IFNL3 variant mRNAs at nucleotide resolution by SHAPE-MaP. The rs4803217 G allele mRNA forms well-defined 3' UTR structure while the T allele mRNA is more dynamic. The observed differences between alleles are among the largest possible RNA structural alterations that can be induced by a single nucleotide change and transform the UTR from a single well-defined conformation to one with multiple dynamic interconverting structures. These data illustrate that non-coding genetic variants can have significant functional effects by impacting RNA structure.

Prevalence of genetic variants of keratins 8 and 18 in patients with drug-induced liver injury.

BACKGROUND: Keratin 8 and 18 (K8/K18) cytoskeletal proteins protect hepatocytes from undergoing apoptosis and their mutations predispose to adverse outcomes in acute liver failure (ALF). All known K8/K18 variants occur at relatively non-conserved residues and do not cause keratin cytoskeleton reorganization, whereas epidermal keratin-conserved residue mutations disrupt the keratin cytoskeleton and cause severe skin disease. The aim of our study was to identify keratin variants in idiosyncratic drug-induced liver injury (DILI). METHODS: Genomic DNA was isolated from 800 patients enrolled in an ongoing US multicenter study, with DILI attributed to a wide range of drugs. Specific K8/K18 exonic regions were PCR-amplified and screened by denaturing HPLC followed by DNA sequencing. The functional impact of keratin variants was assessed using cell transfection and immune staining. RESULTS: Heterozygous and compound amino acid-altering K8/K18 variants were identified in 86 DILI patients and non-coding variants in 15 subjects. Five novel amino acid-altering (K8 Lys393Arg, K8 Ala351Val, K8 Ala358Val, K8 Ile346Val, K18 Asp89His) and two non-coding variants were observed. Several variants segregated with specific ethnic backgrounds but were found at similar frequencies in DILI subjects and ethnically matched population controls. Notably, variants in highly conserved residues of K8 Lys393Arg (ezetimibe/simvastatin-related) and K18 Asp89His (isoniazid-related) were found in patients with fatal DILI. These novel variants also led to keratin network disruption in transfected cells. CONCLUSIONS: Novel K8/K18 cytoskeleton-disrupting variants were identified in two patients and segregated with fatal DILI. Other non-cytoskeleton-disrupting keratin variants did not preferentially associate with DILI.

Toward predicting drug-induced liver injury: parallel computational approaches to identify multidrug resistance protein 4 and bile salt export pump inhibitors.

Drug-induced liver injury (DILI) is an important cause of drug toxicity. Inhibition of multidrug resistance protein 4 (MRP4), in addition to bile salt export pump (BSEP), might be a risk factor for the development of cholestatic DILI. Recently, we demonstrated that inhibition of MRP4, in addition to BSEP, may be a risk factor for the development of cholestatic DILI. Here, we aimed to develop computational models to delineate molecular features underlying MRP4 and BSEP inhibition. Models were developed using 257 BSEP and 86 MRP4 inhibitors and noninhibitors in the training set. Models were externally validated and used to predict the affinity of compounds toward BSEP and MRP4 in the DrugBank database. Compounds with a score above the median fingerprint threshold were considered to have significant inhibitory effects on MRP4 and BSEP. Common feature pharmacophore models were developed for MRP4 and BSEP with LigandScout software using a training set of nine well characterized MRP4 inhibitors and nine potent BSEP inhibitors. Bayesian models for BSEP and MRP4 inhibition/noninhibition were developed with cross-validated receiver operator curve values greater than 0.8 for the test sets, indicating robust models with acceptable false positive and false negative prediction rates. Both MRP4 and BSEP inhibitor pharmacophore models were characterized by hydrophobic and hydrogen-bond acceptor features, albeit in distinct spatial arrangements. Similar molecular features between MRP4 and BSEP inhibitors may partially explain why various drugs have affinity for both transporters. The Bayesian (BSEP, MRP4) and pharmacophore (MRP4, BSEP) models demonstrated significant classification accuracy and predictability.

Population pharmacokinetic modeling of plasma and intracellular ribavirin concentrations in patients with chronic hepatitis C virus infection.

Ribavirin, a guanosine analog, is a broad-spectrum antiviral agent. Ribavirin has been a fundamental component of the treatment of hepatitis C virus (HCV) infection for decades, but there is a very limited understanding of the clinical pharmacology of this drug. Furthermore, it is associated with a major dose-limiting toxicity, hemolytic anemia. Ribavirin undergoes intracellular phosphorylation by host enzymes to ribavirin monophosphate (RMP), ribavirin diphosphate (RDP), and ribavirin triphosphate (RTP). The intracellular forms have been associated with antiviral and toxic effects in vitro, but the kinetics of these phosphorylated moieties have not been fully elucidated in vivo. We developed a model to characterize the plasma pharmacokinetics of ribavirin and the difference between intracellular phosphorylation kinetics in red cells (nonnucleated) and in peripheral blood mononuclear cells (nucleated). A time-independent two-compartment model with first-order absorption described the plasma data well. The cellular phosphorylation kinetics was described by a one-compartment model for RMP, with the formation rate driven by plasma concentrations and the first-order degradation rate. RDP and RTP rapidly reached equilibrium with RMP. Concomitant telaprevir use, inosine triphosphatase genetics, creatinine clearance, weight, and sex were significant covariates. The terminal ribavirin half-life in plasma and phosphorylated anabolites in cells was approximately 224 h. We found no evidence of time-dependent kinetics. These data provide a foundation for uncovering concentration-effect associations for ribavirin and determining the optimal dose and duration of this drug for use in combination with newer direct-acting HCV agents. (This study has been registered at ClinicalTrials.gov under registration no. NCT01097395.).