Assessing the Performance of In silico Tools and Molecular Dynamics Simulations for Predicting Pharmacogenetic Variant Impact.

The ability of freely available in silico tools to predict the effect of non-synonymous single nucleotide polymorphisms (nsSNPs) in pharmacogenes on protein function is not well defined. We assessed the performance of seven sequence-based (SIFT, PolyPhen2, mutation accessor, FATHMM, PhD-SNP, MutPred2, and SNPs & Go) and five structure-based (mCSM, SDM, DDGun, CupSat, and MAESTROweb) tools in predicting the impact of 118 nsSNPs in the CYP2C19, CYP2C9, CYP2B6, CYP2D6, and DPYD genes with known function (24 normal, one increased, 42 decreased, and 51 no-function). Sequence-based tools had a higher median (IQR) positive predictive value (89% [89-94%] vs. 12% [10-15%], P < 0.001) and lower negative predictive value (30% [24-34%] vs. 90% [80-93%], P < 0.001) than structure-based tools. Accuracy did not significantly differ between sequence-based (59% [37-67%]) and structure-based (34% [23-44%]) tools (P = 0.070). Notably, the no-function CYP2C9*3 allele and decreased function CYP2C9*8 allele were predicted incorrectly as tolerated by 100% of sequenced-based tools and as stabilizing by 60% and 20% of structure-based tools, respectively. As a case study, we performed mutational analysis for the CYP2C9*1, *3 (I359L), and *8 (R150H) proteins through molecular dynamic (MD) simulations using S-warfarin as the substrate. The I359L variant increased the distance of the major metabolic site of S-warfarin to the oxy-ferryl center of CYP2C9, and I359L and R150H caused shifts in the conformation of S-warfarin to a position less favorable for metabolism. These data suggest that MD simulations may better capture the impact of nsSNPs in pharmacogenes than other tools.

Circulating microRNA Biomarkers of Thiazide Response in Hypertension.

BACKGROUND: Thiazide diuretics are the second most frequently prescribed class of antihypertensives, but up to 50% of patients with hypertension have minimal antihypertensive response to thiazides. We explored circulating microRNAs (miRNAs) in search of predictive biomarkers of thiazide response. METHODS AND RESULTS: We profiled 754 miRNAs in baseline plasma samples of 36 hypertensive European American adults treated with hydrochlorothiazide, categorized into responders (n=18) and nonresponders (n=18) on the basis of diastolic blood pressure response to hydrochlorothiazide. miRNAs with ≥2.5-fold differential expression between responders and nonresponders were considered for validation in 3 cohorts (n=50 each): hydrochlorothiazide-treated European Americans, chlorthalidone-treated European Americans, and hydrochlorothiazide-treated Black individuals. Different blood pressure phenotypes including categorical (responder versus nonresponder) and continuous diastolic blood pressure and systolic blood pressure were tested for association with the candidate miRNA expression using multivariate regression analyses adjusting for age, sex, and baseline blood pressure. After quality control, 74 miRNAs were available for screening, 19 of which were considered for validation. In the validation cohort, miR-193b-3p and 30d-5p showed significant associations with continuous SBP or diastolic blood pressure response or both, to hydrochlorothiazide in European Americans at Benjamini-Hochberg adjusted P<0.05. In the combined analysis of validation cohorts, let-7g (odds ratio, 0.6 [95% CI, 0.4-0.8]), miR-142-3p (odds ratio, 1.1 [95% CI, 1.0, 1.2]), and miR-423-5p (odds ratio, 0.7 [95% CI, 0.5-0.9]) associated with categorical diastolic blood pressure response at Benjamini-Hochberg adjusted P<0.05. Predicted target genes of the 5 miRNAs were mapped to key hypertension pathways: lysine degradation, fatty acid biosynthesis, and metabolism. CONCLUSIONS: The above identified circulating miRNAs may have a potential for clinical use as biomarkers for thiazide diuretic selection in hypertension. REGISTRATION: URL: ClinicalTrials.gov. Unique identifiers: NCT00246519, NCT01203852, NCT00005520.

Association of tenofovir diphosphate and lamivudine triphosphate concentrations with HIV and hepatitis B virus viral suppression.

OBJECTIVE: Concentrations of tenofovir diphosphate (TFV-DP) and lamivudine triphosphate (3TC-TP) in cells are correlates of medication adherence and antiviral activity. However, studies have yet to characterize the simultaneous relationship between TFV-DP and 3TC-TP concentrations with HIV and hepatitis B virus (HBV) suppression. METHODS: Individuals with HIV/HBV coinfection on tenofovir disoproxil fumarate (TDF)-containing antiretroviral therapy (ART) were enrolled. Peripheral blood mononuclear cells (PBMCs) and dried blood spots (DBS) samples were collected and steady-state TFV-DP and 3TC-TP concentrations quantified using validated methods. The relationship between patient factors, TFV-DP, and 3TC-TP concentrations in PBMCs and DBS with HBV and HIV viral suppression were examined. RESULTS: Of 138 participants on TDF-containing ART for a median duration (range) of 6 (0.75-15) years, the median age was 43 years and 64% were women. Overall, 128 (92.8%) and 129 (93.5%) had suppressed HIV and HBV viral loads, respectively. Of the 128 participants with suppressed HIV, 122 (95.3%) had suppressed HBV. Self-reported ART adherence, recent change to dolutegravir-based ART, TFV-DP, and 3TC-TP concentrations in PBMCs and DBS were associated with HIV RNA suppression, while HBe antigen positivity, HIV suppression, and TFV-DP concentrations in DBS were associated with HBV DNA suppression (including six persons with HBV nonsuppression and HIV suppression). CONCLUSION: Long-term TDF/3TC-conatining ART was highly efficacious in individuals with HIV/HBV coinfection. Higher TFV-DP concentrations were predictive of suppression for both viruses. Persistent HBV viremia on TDF/3TC-containg ART requires additional research, but may represent poor adherence and the need for adherence interventions or novel antivirals.

A pilot study to identify pharmacogenomic and clinical risk factors associated with opioid related falls and adverse effects in older adults.

Given the high prevalence of pain in older adults and current trends in opioid prescribing, inclusion of genetic information in risk prediction tools may improve opioid risk assessment. Our objectives were to (1) determine the feasibility of recruiting socioeconomically disadvantaged and racially diverse middle aged and older adult populations for a study seeking to identify risk factors for opioid-related falls and other serious adverse effects and (2) explore potential associations between the Risk Index for Overdose or Serious Opioid-induced Respiratory Depression (CIP-RIOSORD) risk class and other patient factors with falls and serious opioid adverse effects. This was an observational study of 44 participants discharged home from the emergency department with an opioid prescription for acute pain and followed for 30 days. We found pain interference may predict opioid-related falls or serious adverse effects within older, opioid-treated patients. If validated, pain interference may prove to be a beneficial marker for risk stratification of older adults initiated on opioids for acute pain.

CYP2D6-guided opioid therapy for adults with cancer pain: A randomized implementation clinical trial.

INTRODUCTION: The CYP2D6 enzyme metabolizes opioids commonly prescribed for cancer-related pain, and CYP2D6 polymorphisms may contribute to variability in opioid response. We evaluated the feasibility of implementing CYP2D6-guided opioid prescribing for patients with cancer and reported pilot outcome data. METHODS: Adult patients from two cancer centers were prospectively enrolled into a hybrid implementation-effectiveness clinical trial and randomized to CYP2D6-genotype-guided opioid selection, with clinical recommendations, or usual care. Implementation metrics, including provider response, medication changes consistent with recommendations, and patient-reported pain and symptom scores at baseline and up to 8 weeks, were assessed. RESULTS: Most (87/114, 76%) patients approached for the study agreed to participate. Of 85 patients randomized, 71% were prescribed oxycodone at baseline. The median (range) time to receive CYP2D6 test results was 10 (3-37) days; 24% of patients had physicians acknowledge genotype results in a clinic note. Among patients with CYP2D6-genotype-guided recommendations to change therapy (n = 11), 18% had a change congruent with recommendations. Among patients who completed baseline and follow-up questionnaires (n = 48), there was no difference in change in mean composite pain score (-1.01 ± 2.1 vs. -0.41 ± 2.5; p = 0.19) or symptom severity at last follow-up (3.96 ± 2.18 vs. 3.47 ± 1.78; p = 0.63) between the usual care arm (n = 26) and genotype-guided arm (n = 22), respectively. CONCLUSION: Our study revealed high acceptance of pharmacogenetic testing as part of a clinical trial among patients with cancer pain. However, provider response to genotype-guided recommendations was low, impacting assessment of pain-related outcomes. Addressing barriers to utility of pharmacogenetics results and clinical recommendations will be critical for implementation success.

Characterizing OPRM1 DNA methylation in prescription opioid users with chronic musculoskeletal pain.

Introduction: Many patients with chronic pain use prescription opioids. Epigenetic modification of the µ-opioid receptor 1 (OPRM1) gene, which codes for the target protein of opioids, may influence vulnerability to opioid abuse and response to opioid pharmacotherapy, potentially affecting pain outcomes. Objective: Our objective was to investigate associations of clinical and sociodemographic factors with OPRM1 DNA methylation in patients with chronic musculoskeletal pain on long-term prescription opioids. Methods: Sociodemographic variables, survey data (Rapid Estimate of Adult Health Literacy in Medicine-Short Form, Functional Comorbidity Index [FCI], PROMIS 43v2.1 Profile, Opioid Risk Tool, and PROMIS Prescription Pain Medication Misuse), and saliva samples were collected. The genomic DNA extracted from saliva samples were bisulfite converted, amplified by polymerase chain reaction, and processed for OPRM1-targeted DNA methylation analysis on a Pyrosequencing instrument (Qiagen Inc, Valencia, CA). General linear models were used to examine the relationships between the predictors and OPRM1 DNA methylation. Results: Data from 112 patients were analyzed. The best-fitted multivariable model indicated, compared with their counterparts, patients with > eighth grade reading level, degenerative disk disease, substance abuse comorbidity, and opioid use < 1 year (compared with >5 years), had average methylation levels that were 7.7% (95% confidence interval [CI] 0.95%, 14.4%), 11.7% (95% CI 2.7%, 21.1%), 21.7% (95% CI 10.7%, 32.5%), and 16.1% (95% CI 3.3%, 28.8%) higher than the reference groups, respectively. Methylation levels were 2.2% (95% CI 0.64%, 3.7%) lower for every 1 unit increase in FCI and greater by 0.45% (95% CI 0.08%, 0.82%) for every fatigue T score unit increase. Conclusions: OPRM1 methylation levels varied by several patient factors. Further studies are warranted to replicate these findings and determine potential clinical utility.

Comparison of Data Normalization Strategies for Array-Based MicroRNA Profiling Experiments and Identification and Validation of Circulating MicroRNAs as Endogenous Controls in Hypertension.

Introduction: MicroRNAs are small noncoding RNAs with potential regulatory roles in hypertension and drug response. The presence of many of these RNAs in biofluids has spurred investigation into their role as possible biomarkers for use in precision approaches to healthcare. One of the major challenges in clinical translation of circulating miRNA biomarkers is the limited replication across studies due to lack of standards for data normalization techniques for array-based approaches and a lack of consensus on an endogenous control normalizer for qPCR-based candidate miRNA profiling studies. Methods: We conducted genome-wide profiling of 754 miRNAs in baseline plasma of 36 European American individuals with uncomplicated hypertension selected from the PEAR clinical trial, who had been untreated for hypertension for at least one month prior to sample collection. After appropriate quality control with amplification score and missingness filters, we tested different normalization strategies such as normalization with global mean of imputed and unimputed data, mean of restricted set of miRNAs, quantile normalization, and endogenous control miRNA normalization to identify the method that best reduces the technical/experimental variability in the data. We identified best endogenous control candidates with expression pattern closest to the mean miRNA expression in the sample, as well as by assessing their stability using a combination of NormFinder, geNorm, Best Keeper and Delta Ct algorithms under the Reffinder software. The suitability of the four best endogenous controls was validated in 50 hypertensive African Americans from the same trial with reverse-transcription-qPCR and by evaluating their stability ranking in that cohort. Results: Among the compared normalization strategies, quantile normalization and global mean normalization performed better than others in terms of reducing the standard deviation of miRNAs across samples in the array-based data. Among the four strongest candidate miRNAs from our selection process (miR-223-3p, 19b, 106a, and 126-5p), miR-223-3p and miR-126-5p were consistently expressed with the best stability ranking in the validation cohort. Furthermore, the combination of miR-223-3p and 126-5p showed better stability ranking when compared to single miRNAs. Conclusion: We identified quantile normalization followed by global mean normalization to be the best methods in reducing the variance in the data. We identified the combination of miR-223-3p and 126-5p as potential endogenous control in studies of hypertension.

A randomized, cross-over trial of metoprolol succinate formulations to evaluate PK and PD end points for therapeutic equivalence.

There are limited comparison data throughout the dosing interval for generic versus brand metoprolol extended-release (ER) tablets. We compared the pharmacokinetics (PKs) and pharmacodynamics of brand name versus two generic formulations (drugs 1 and 2) of metoprolol ER tablets with different time to maximum concentration (Tmax ) in adults with hypertension. Participants were randomized to equal drug doses (50-150 mg/day) administered in one of two sequences (brand-drug1-brand-drug2 or brand-drug2-brand-drug1) and completed 24-h PK, digital heart rate (HR), ambulatory blood pressure (BP), and HR studies after taking each formulation for greater than or equal to 7 days. Metoprolol concentrations were determined by liquid chromatography tandem mass spectrometry, with noncompartmental analysis performed to obtain PK parameters in Phoenix WinNonlin. Heart rate variability (HRV) low-to-high frequency ratio was determined per quartile over the 24-h period. Thirty-six participants completed studies with the brand name and at least one generic product. Among 30 participants on the 50 mg dose, the primary PK end points of area under the concentration-time curve and Cmax were similar between products; Tmax was 6.1 ± 3.6 for the brand versus 3.5 ± 4.9 for drug 1 (p = 0.019) and 9.6 ± 3.2 for drug 2 (p < 0.001). Among all 36 participants, 24-h BPs and HRs were similar between products. Mean 24-h HRV low-to-high ratio was also similar for drug 1 (2.04 ± 1.35), drug 2 (1.86 ± 1.35), and brand (2.04 ± 1.77), but was more sustained over time for the brand versus drug 1 (drug × quartile interaction p = 0.017). Differences in Tmax between metoprolol ER products following repeated doses may have implications for drug effects on autonomic balance over the dosing interval.

Molecular docking and mouse modeling suggest CMKLR1 and INSR as targets for improving PCOS phenotypes by minocycline.

Polycystic ovary syndrome (PCOS) is the most common cause of women's infertility. Some inflammatory pathways play a pivotal role in the pathogenesis of PCOS. This study aimed to investigate the possible beneficial effects of minocycline on chemokine-like receptor 1 (CMKLR1) and Insulin Receptor (INSR) in a PCOS model. A molecular docking study was implemented using Molecular Operating Environment (MOE) software. The PCOS was induced in NMRI mice (mean body weight 14.47±0.23) by 28 days estradiol valerate injection (2 mg/kg/day). The mice were then divided into six groups (n=8 per group, mean body weight 17.77± 0.26): control (received normal saline), PCOS model, control for minocycline, minocycline treated PCOS (50 mg/kg), letrozole treated PCOS (0.5 mg/kg), and metformin-treated PCOS (300 mg/kg). Serum FSH, LH, estradiol (E2), and testosterone were detected by ELISA. The ovarian tissues were stained by hematoxylin and eosin. The CMKLR1 and INSR expression levels were determined by Real-time-PCR. The molecular docking studies showed scores of -10.92 and -9.30 kcal/mol, respectively, for minocycline with CMKLR1 and INSR. Estradiol valerate treatment led to a significant increase in E2, graffian follicle, and decrease in corpus luteum (CL) numbers (P<0.05), while minocycline treatment improved these PCOS features. The minocycline treatment significantly decreased the CMKLR1 expression and increased the INSR expression (P<0.05) while the CMKLR1 expression was increased in PCOS model. Minocycline may improve ovulation in PCOS model by returning E2 to a normal level and increasing CL number (ovulation signs). These beneficial outcomes may be related to the changes in CMKLR1 and INSR gene expression involved in glucose metabolism and inflammation.

Effect of CYP2C19 genetic variants on bleeding and major adverse cardiovascular events in a cohort of Arab patients undergoing percutaneous coronary intervention and stent implantation.

INTRODUCTION: One-third of patients have clopidogrel resistance that may lead to major adverse cardiac events (MACEs). By contrast, it was found that some clopidogrel-treated patients have hyperresponsive platelets that are associated with higher bleeding risk. Several studies have shown that polymorphisms in the gene encoding the CYP2C19 contribute to the variability in response to clopidogrel. Data on genetic and nongenetic factors affecting clopidogrel response in the Arab population are scarce. In this prospective cohort study, we sought to assess the association between the increased function allele (CYP2C19*17) and bleeding events, and validate the effect of the CYP2C19 genetic variants and nongenetic factors on the incidence of MACEs. METHODS: Blood samples were collected from patients that were undergoing percutaneous coronary intervention and receiving clopidogrel at the Heart Hospital, a specialist tertiary hospital in Doha, Qatar. Patients were followed for 12 months. Genotyping was performed for CYP2C19*2, *3, and *17 using TaqMan assays. RESULTS: In 254 patients, the minor allele frequencies were 0.13, 0.004, and 0.21 for *2, *3, and *17, respectively. Over a 12-month follow-up period, there were 21 bleeding events (8.5 events/100 patient-year). CYP2C19*17 carriers were found to be associated with increased risk of bleeding (OR, 21.6; 95% CI, 4.8-96.8; P < 0.0001). CYP2C19*2 or *3 carriers were found to be associated with increased risk of baseline and incident MACE combined (OR, 8.4; 95% CI, 3.2-23.9; P < 0.0001). CONCLUSION: This study showed a significant association between CYP2C19*17 allele and the increased risk of bleeding, and CYP2C19*2 or *3 with MACE outcomes.

ß1-receptor polymorphisms and junctional ectopic tachycardia in children after cardiac surgery.

Junctional ectopic tachycardia (JET) is a potentially life-threatening postoperative arrhythmia in children with specific congenital heart defects and can contribute significantly to postoperative morbidity for at-risk populations. In adults, ß1-adrenergic receptor (ADRB1) and ß2-adrenergic receptor (ADRB2) genotypes have been associated with increased risk for arrhythmias. However, their association with arrhythmia risk in children is unknown. We aimed to test associations between ADRB1 and ADRB2 genotypes and postoperative JET in patients with congenital heart defects. Children who underwent cardiac surgery were genotyped for the ADRB1 p.Ser49Gly (rs1801252; c.145A>G), p.Arg389Gly (rs1801253; c.1165C>G), ADRB2 p.Arg16Gly (rs1042713; c.46A>G), and p.Glu27Gln (rs1042714; c.79G>C) polymorphisms. The occurrence of postoperative JET was assessed via cardiologist-interpreted electrocardiograms. Genotype associations with JET were analyzed via logistic regression, adjusted for clinical variables associated with JET, with separate analysis in patients not on a ß-blocker. Of the 343 children included (median age 8 months, 53% boys, 69% European ancestry), 45 (13%) developed JET. The Arg389Arg genotype was not significantly associated with JET in the overall population (odds ratio [OR] = 1.96, 95% confidence interval [CI] = 0.96-4.03, p = 0.064), but was nominally associated in patients not taking a ß-blocker (n = 324, OR = 2.25, 95% CI = 1.05-4.80. p = 0.034). None of the other variants were associated with JET. These data suggest that the ADRB1 Arg389Arg genotype may predict risk for JET following cardiac surgery in pediatric patients in the absence of ß-blockade. Whether treatment with a ß-blocker ameliorates this association requires further research.

Exploring VEGF-Linked Pathways: Investigating Multiple miRNAs for Their Therapeutic Potential in Angiogenesis Targets and as Biomarkers in Recurrent Glioblastoma Multiforme.

Alternative pathways frequently operate as the origins of resistance to drugs that block the vascular endothelial growth factor (VEGF) pathway. To find possible therapeutic targets and indicators, this study explored the VEGF pathway and how miRNAs control it in recurrent glioblastoma multiforme (rGBM). Differentially expressed miRNAs (DEmiRNAs) were identified by using GBM GSE profiles (GSE32466). To find pathways containing DEmiRNAs, VEGF pathway genes, and their related genes, DIANA-miRPath v3.0 and the ToppGene database were utilized. miRNAs linked to VEGF signaling pathway genes, interactional genes, and DEmiRNAs were discovered by extracting common pathways. The ability of these miRNAs to distinguish rGBM patients from those with primary GBM was assessed using ROC analysis. The study revealed that in rGBM, 30 miRNAs were significantly up-regulated and 49 miRNAs were considerably down-regulated. Among them, the VEGF pathway was connected to 22 up-regulated miRNAs and 29 down-regulated miRNAs. The MAPK pathway shared the most genes with the VEGF pathway, accounting for 1,014 of the interacting genes, which were discovered to have interactions with VEGF signaling pathway genes. Furthermore, 14 miRNAs were identified as having a great deal of potential as molecular biomarkers and therapeutic targets for rGBM. The results indicate that the VEGF pathway in rGBM is regulated by a number of interrelated pathways. The discovered miRNAs hold promise as rGBM biomarkers and therapeutic targets, offering possibilities for novel therapy strategies and aiding rGBM diagnosis and prognosis.

Genome-wide Association Study Identified Chromosome 8 Locus Associated with Medication-Related Osteonecrosis of the Jaw.

Medication-related osteonecrosis of the jaw (MRONJ) is a rare but serious drug-related adverse event. To identify pharmacogenomic markers of MRONJ associated with bisphosphonate therapy, we conducted a genomewide association study (GWAS) meta-analysis followed by functional analysis of 5,008 individuals of European ancestry treated with bisphosphonates, which includes the largest number of MRONJ cases to date (444 cases and 4,564 controls). Discovery GWAS was performed in randomly selected 70% of the patients with cancer and replication GWAS was performed in the remaining 30% of the patients with cancer treated with intravenous bisphosphonates followed by meta-analysis of all 3,639 patients with cancer. GWAS was also performed in 1,369 patients with osteoporosis treated with oral bisphosphonates. The lead single-nucleotide polymorphism (SNP), rs2736308 on chromosome 8, was associated with an increased risk of MRONJ with an odds ratio (OR) of 2.71 and 95% confidence interval (CI) of 1.90-3.86 (P = 3.57*10-8 ) in the meta-analysis of patients with cancer. This SNP was validated in the MRONJ GWAS in patients with osteoporosis (OR: 2.82, 95% CI: 1.55-4.09, P = 6.84*10-4 ). The meta-analysis combining patients with cancer and patients with osteoporosis yielded the same lead SNP rs2736308 on chromosome 8 as the top SNP (OR: 2.74, 95% CI: 2.09-3.39, P = 9.65*10-11 ). This locus is associated with regulation of the BLK, CTSB, and FDFT1 genes, which had been associated with bone mineral density. FDFT1 encodes a membrane-associated enzyme, which is implicated in the bisphosphonate pathway. This study provides insights into the potential mechanism of MRONJ.

Lactate Dehydrogenase B and Pyruvate Oxidation Pathway Associated With Carfilzomib-Related Cardiotoxicity in Multiple Myeloma Patients: Result of a Multi-Omics Integrative Analysis.

Multiple myeloma (MM) is the second most frequent hematologic cancer in the United States. Carfilzomib (CFZ), an irreversible proteasome inhibitor being used to treat relapsed and refractory MM, has been associated with cardiotoxicity, including heart failure. We hypothesized that a multi-omics approach integrating data from different omics would provide insights into the mechanisms of CFZ-related cardiovascular adverse events (CVAEs). Plasma samples were collected from 13 MM patients treated with CFZ (including 7 with CVAEs and 6 with no CVAEs) at the University of Florida Health Cancer Center. These samples were evaluated in global metabolomic profiling, global proteomic profiling, and microRNA (miRNA) profiling. Integrative pathway analysis was performed to identify genes and pathways differentially expressed between patients with and without CVAEs. The proteomics analysis identified the up-regulation of lactate dehydrogenase B (LDHB) [fold change (FC) = 8.2, p = 0.01] in patients who experienced CVAEs. The metabolomics analysis identified lower plasma abundance of pyruvate (FC = 0.16, p = 0.0004) and higher abundance of lactate (FC = 2.4, p = 0.0001) in patients with CVAEs. Differential expression analysis of miRNAs profiling identified mir-146b to be up-regulatein (FC = 14, p = 0.046) in patients with CVAE. Pathway analysis suggested that the pyruvate fermentation to lactate pathway is associated with CFZ-CVAEs. In this pilot multi-omics integrative analysis, we observed the down-regulation of pyruvate and up-regulation of LDHB among patients who experienced CVAEs, suggesting the importance of the pyruvate oxidation pathway associated with mitochondrial dysfunction. Validation and further investigation in a larger independent cohort are warranted to better understand the mechanisms of CFZ-CVAEs.

Pharmacogenetic predictors of nevirapine pharmacokinetics in Ghanaian children living with HIV with or without TB coinfection.

Nevirapine (NVP) is a non-nucleoside reverse transcriptase inhibitor that is used in the treatment of human immunodeficiency virus (HIV) infection in children younger than 3 years old. Identifying genetic predictors of NVP pharmacokinetics (PK) in young children is important because inter-individual variability in NVP concentrations contributes to variable treatment response and the information may be used to individualize dosing decisions. We examined the relationship between genetic variations in relevant drug disposition genes and NVP PK parameters in Ghanaian children living with HIV eligible to initiate NVP-based antiretroviral therapy. Participants received NVP plus zidovudine and lamivudine or abacavir and lamivudine twice daily, and those with tuberculosis (TB) coinfection received concurrent anti-TB therapy with NVP. Pharmacokinetic sampling was performed after at least 4 weeks of antiretroviral therapy. Nevirapine minimum concentration (Cmin), area under the concentration-time curve from time 0 to 12 h (AUC0-12h), and apparent clearance (CL/F) were calculated using non-compartmental analysis using Phoenix v8.0 software. Genotyping for CYP2B6, CYP2A6, CYP3A5, ABCB1, NR1I2, and NR1I3 single nucleotide polymorphism (SNPs) was performed by TaqMan® allelic discrimination method. The median (range) NVP dose received was 10 (7-14) mg/kg. Of the 53 participants, the median (range) Cmin was 3.3 (0.0-14.0) mg/L and AUC0-12h was 56.0 (16.7-202.6) mg.hr/L. Using step-wise regression, CYP2B6 rs3745274 and NR1I2 rs6785049 SNPs were independent as well as joint predictors of NVP AUC0-12h, Cmin, and CL/F. We concluded that genotyping for CYP2B6 rs3745274, and the NR1I2 rs6785049 G > A SNP (which encodes the transcriptional factor, pregnane X receptor), could improve prediction of NVP PK for individualized therapy.

Newly diagnosed cardiovascular disease in patients treated with immune checkpoint inhibitors: a retrospective analysis of patients at an academic tertiary care center.

BACKGROUND: Immune checkpoint inhibitors (ICIs) are a novel class of anticancer agents that have demonstrated clinical response for both solid and hematological malignancies. ICIs are associated with development of immune-related adverse events including cardiotoxicity. We estimated the incidence of newly diagnosed cardiovascular disease in patients treated with ICIs at a large, tertiary care center. METHODS: All patients with a cancer diagnosis who received any ICI treatment in the University of Florida's Integrated Data Repository from 2011 to 2017 were included. Cardiovascular disease was defined as a new ICD diagnosis code for cardiomyopathy, heart failure, arrhythmia, heart block, pericardial disease, or myocarditis after initiation of ICI treatment. RESULTS: Of 102,701 patients with a diagnosis of malignancy, 424 patients received at least one ICI. Sixty-two (14.6%) patients were diagnosed with at least one new cardiovascular disease after initiation of ICI therapy. Of the 374 patients receiving one ICI, 21 (5.6%) developed heart failure. Of the 49 patients who received two ICIs sequentially, three (6.1%) developed heart failure and/or cardiomyopathy. Incident cardiovascular disease was diagnosed at a median of 63 days after initial ICI exposure. One patient developed myocarditis 28 days after receiving nivolumab. Mortality in ICI treated patients with a concomitant diagnosis of incident cardiovascular disease was higher compared to those who did not (66.1% vs. 41.4%, odds ratio = 2.77, 1.55-4.95, p = 0.0006). CONCLUSIONS: This study suggests a high incidence of newly diagnosed cardiovascular disease after the initiation of ICI therapy in a real-world clinical setting.

Implications of Polymorphisms in the BCKDK and GATA-4 Gene Regions on Stable Warfarin Dose in African Americans.

VKORC1 and CYP2C9 genotypes explain less variability in warfarin dose requirements in African Americans compared with Europeans. Variants in BCKDK and GATA-4 gene regions, purported to regulate VKORC1 and CYP2C9 expression, have been shown to play an important role in warfarin dose requirements in Europeans and Asians, respectively. We sought to determine whether rs56314408 near BCKDK or GATA-4 rs2645400 influence warfarin dose requirements in 200 African Americans. Unlike the strong linkage disequilibrium (LD) between rs56314408 and VKORC1 rs9923231 in Europeans, they were not in LD in African Americans. No associations were found on univariate analysis. On multivariable analysis, rs56314408 was associated (P = 0.027) with dose in a regression model excluding VKORC1 rs9923231, and GATA-4 rs2645400 was associated (P = 0.032) with dose in a model excluding CYP2C (CYP2C9*2, *3, *5, *6, *8, and *11, CYP2C rs12777823) variants. Neither variant contributed to dose in the model that included both VKORC1 rs9923231 and CYP2C variants. Our results do not support contributions of the studied variants to warfarin dose requirements in African Americans. However, they illustrate the value of studies in African descent populations, who have low LD in their genome, in teasing out genetic variation underlying drug response associations. They also emphasize the importance of confirming associations in persons of African ancestry.

SIRT1 Gene SNP rs932658 Is Associated With Medication-Related Osteonecrosis of the Jaw.

Medication-related osteonecrosis of the jaw (MRONJ) is a rare but serious adverse drug reaction. Our previous whole-exome sequencing study found SIRT1 intronic region single-nucleotide polymorphism (SNP) rs7896005 to be associated with MRONJ in cancer patients treated with intravenous (iv) bisphosphonates (BPs). This study aimed to identify causal variants for this association. In silico analyses identified three SNPs (rs3758391, rs932658, and rs2394443) in the SIRT1 promoter region that are in high linkage disequilibrium (r2 > 0.8) with rs7896005. To validate the association between these SNPs and MRONJ, we genotyped these three SNPs on the germline DNA from 104 cancer patients of European ancestry treated with iv BPs (46 cases and 58 controls). Multivariable logistic regression analysis showed the minor alleles of these three SNPs were associated with lower odds for MRONJ. The odds ratios (95% confidence interval) and p values were 0.351 (0.164-0.751; p = 0.007) for rs3758391, 0.351 (0.164-0.751; p = 0.007) for rs932658, and 0.331 (0.157-0.697; p = 0.0036) for rs2394443, respectively. In the reporter gene assays, constructs containing rs932658 with variant allele A had higher luciferase activity than the reference allele, whereas constructs containing SNP rs3758391 and/or rs2394443 did not significantly affect activity. These results indicate that the promoter SNP rs932658 regulates the expression of SIRT1 and presumably lowers the risk of MRONJ by increasing SIRT1 expression. © 2020 American Society for Bone and Mineral Research (ASBMR).

Examination of Metoprolol Pharmacokinetics and Pharmacodynamics Across CYP2D6 Genotype-Derived Activity Scores.

Recent CYP2D6 phenotype standardization efforts by CYP2D6 activity score (AS) are based on limited pharmacokinetic (PK) and pharmacodynamic (PD) data. Using data from two independent clinical trials of metoprolol, we compared metoprolol PK and PD across CYP2D6 AS with the goal of determining whether the PK and PD data support the new phenotype classification. S-metoprolol apparent oral clearance (CLo), adjusted for clinical factors, was correlated with CYP2D6 AS (P < 0.001). The natural log of CLo was lower with an AS of 1 (7.6 ± 0.4 mL/minute) vs. 2-2.25 (8.3 ± 0.6 mL/minute; P = 0.012), similar between an AS of 1 and 1.25-1.5 (7.8 ± 0.5 mL/minute; P = 0.702), and lower with an AS of 1.25-1.5 vs. 2-2.25 (P = 0.03). There was also a greater reduction in heart rate with metoprolol among study participants with AS of 1 (-10.8 ± 5.5) vs. 2-2.25 (-7.1 ± 5.6; P < 0.001) and no significant difference between those with an AS of 1 and 1.25-1.5 (-9.2 ± 4.7; P = 0.095). These data highlight linear trends among CYP2D6 AS and metoprolol PK and PD, but inconsistencies with the phenotypes assigned by AS based on the current standards. Overall, this case study with metoprolol suggests that utilizing CYP2D6 AS, instead of collapsing AS into phenotype categories, may be the most precise approach for utilizing CYP2D6 pharmacogenomics in clinical practice.

Age-Related Changes in miRNA Expression Influence GSTZ1 and Other Drug Metabolizing Enzymes.

Previous work has shown that hepatic levels of human glutathione transferase zeta 1 (GSTZ1) protein, involved in tyrosine catabolism and responsible for metabolism of the investigational drug dichloroacetate, increase in cytosol after birth before reaching a plateau around age 7. However, the mechanism regulating this change of expression is still unknown, and previous studies showed that GSTZ1 mRNA levels did not correlate with GSTZ1 protein expression. In this study, we addressed the hypothesis that microRNAs (miRNAs) could regulate expression of GSTZ1. We obtained liver samples from donors aged less than 1 year or older than 13 years and isolated total RNA for use in a microarray to identify miRNAs that were downregulated in the livers of adults compared with children. From a total of 2578 human miRNAs tested, 63 miRNAs were more than 2-fold down-regulated in adults, of which miR-376c-3p was predicted to bind to the 3' untranslated region of GSTZ1 mRNA. There was an inverse correlation of miR-376c-3p and GSTZ1 protein expression in the liver samples. Using cell culture, we confirmed that miR-376c-3p could downregulate GSTZ1 protein expression. Our findings suggest that miR-376c-3p prevents production of GSTZ1 through inhibition of translation. These experiments further our understanding of GSTZ1 regulation. Furthermore, our array results provide a database resource for future studies on mechanisms regulating human hepatic developmental expression. SIGNIFICANCE STATEMENT: Hepatic glutathione transferase zeta 1 (GSTZ1) is responsible for metabolism of the tyrosine catabolite maleylacetoacetate as well as the investigational drug dichloroacetate. Through examination of microRNA (miRNA) expression in liver from infants and adults and studies in cells, we showed that expression of GSTZ1 is controlled by miRNA. This finding has application to the dosing regimen of the drug dichloroacetate. The miRNA expression profiles are provided and will prove useful for future studies of drug-metabolizing enzymes in infants and adults.