Physiologically-Based Pharmacokinetic Modelling of Creatinine-Drug Interactions in the Chronic Kidney Disease Population

Utilising Endogenous Biomarkers in Drug Development to Streamline the Assessment of Drug-Drug Interactions Mediated by Renal Transporters: A Pharmaceutical Industry Perspective

The renal secretion of many drugs is facilitated by membrane transporters, including organic cation transporter 2, multidrug and toxin extrusion protein 1/2-K and organic anion transporters 1 and 3. Inhibition of these transporters can reduce renal excretion of drugs and thereby pose a safety risk. Assessing the risk of inhibition of these membrane transporters by investigational drugs remains a key focus in the evaluation of drug-drug interactions (DDIs). Current methods to predict DDI risk are based on generating in vitro data followed by a clinical assessment using a recommended exogenous probe substrate for the individual drug transporter. More recently, monitoring plasma-based and urine-based endogenous biomarkers to predict transporter-mediated DDIs in early phase I studies represents a promising approach to facilitate, improve and potentially avoid conventional clinical DDI studies. This perspective reviews the evidence for use of these endogenous biomarkers in the assessment of renal transporter-mediated DDI, evaluates how endogenous biomarkers may help to expand the DDI assessment toolkit and offers some potential knowledge gaps. A conceptual framework for assessment that may complement the current paradigm of predicting the potential for renal transporter-mediated DDIs is outlined.

Triglyceride-Glucose Index and Its Correlates: Associations with Serum Creatinine and Estimated Glomerular Filtration Rate in a Cross-Sectional Study from CHARLS 2011-2015

Background: Chronic kidney disease (CKD) has emerged as a significant global public health challenge, and the estimated glomerular filtration rate (eGFR) is widely used due to its convenience, low cost, and broad clinical applicability. Concurrently, insulin resistance (IR) serves as a crucial marker of metabolic disturbance, and alternative indicators have garnered increasing attention in CKD research in recent years. Objective: This study aims to investigate the relationship between IR-related indices (TyG index, TyG-BMI index, and TyG-WC index) and serum creatinine levels, as well as the eGFR, with the intention of uncovering their potential roles in the assessment of renal function. Methods: We analyzed nationally representative cross-sectional data from a cohort of individuals aged 45 and above in China, comprising 11,608 participants. Participants were categorized into different groups based on quartiles of the TyG index, and multiple factors, including gender, age, lifestyle, and co-morbidities, were adjusted for using linear regression models. Results: By linear regression, TyG, TyG-BMI, and TyG-WC indices were significantly positively correlated with serum creatinine and significantly negatively correlated with eGFR. Results showed similar trends when TyG, TyG-BMI, and TyG-WC indices were used as categorical variables. In the fully adjusted model, the highest quartile of serum creatinine was higher than the first quartile for TyG, TyG-BMI, and TyG-WC indices, with β values of 2.673, 3.67, and 1.937 mg/dL, respectively; the highest quartile of eGFR was lower than the first quartile, with β values of -2.4, -2.955, and -1.823 mL/min/1.73 m2. P values were statistically significant. Conclusions: This study indicates a consistent correlation between the TyG index and its related indices with serum creatinine levels and eGFR among the middle aged and elderly population in China. These findings suggest the potential utility of these indices in early screening and management of the risk of chronic kidney disease.

Pharmacological Considerations in the Interpretation of Biochemical Results in Diabetic Patients with Cardiovascular Complications

Diabetes mellitus with cardiovascular diseases is often a multi-systemic disease that requires a multi-therapeutic approach which mostly poses a challenge to laboratory result interpretation. The non-availability of information on many patients due to poor referral, documentation and record keeping has resulted in isolated interpretation of laboratory result of diabetic patients with multisystemic complications. This has led to both analytical and post-analytical errors which has a negative impact on total quality of results. Therefore, this review showed the possible therapeutic treatment of a diabetic patient with cardiovascular disease and how their pharmacological role could affect laboratory result.

Incorporating Uremic Solute-mediated Inhibition of OAT1/3 Improves PBPK Prediction of Tenofovir Renal and Systemic Disposition in Patients with Severe Kidney Disease

BACKGROUND

Dose modification of renally secreted drugs in patients with chronic kidney disease (CKD) has relied on serum creatinine concentration as a biomarker to estimate glomerular filtration (GFR) under the assumption that filtration and secretion decline in parallel. A discrepancy between actual renal clearance and predicted renal clearance based on GFR alone is observed in severe CKD patients with tenofovir, a compound secreted by renal OAT1/3. Uremic solutes that inhibit OAT1/3 may play a role in this divergence.

METHODS

To examine the impact of transporter inhibition by uremic solutes on tenofovir renal clearance, we determined the inhibitory potential of uremic solutes hippuric acid, indoxyl sulfate, and p-cresol sulfate. The inhibition parameters (IC50) were incorporated into a previously validated mechanistic kidney model; simulated renal clearance and plasma PK profile were compared to data from clinical studies.

RESULTS

Without the incorporation of uremic solute inhibition, the PBPK model failed to capture the observed data with an absolute average fold error (AAFE) > 2. However, when the inhibition of renal uptake transporters and uptake transporters in the slow distribution tissues were included, the AAFE value was within the pre-defined twofold model acceptance criterion, demonstrating successful model extrapolation to CKD patients.

CONCLUSION

A PBPK model that incorporates inhibition by uremic solutes has potential to better predict renal clearance and systemic disposition of secreted drugs in patients with CKD. Ongoing research is warranted to determine if the model can be expanded to include other OAT1/3 substrate drugs and to evaluate how these findings can be translated to clinical guidance for drug selection and dose optimization in patients with CKD.

Quantitative Consideration of Clinical Increases in Serum Creatinine Caused by Renal Transporter Inhibition

Creatinine is a common biomarker of renal function and is secreted in the renal tubular cells via drug transporters, such as organic cation transporter 2 and multidrug and toxin extrusion (MATE) 1/2-K. To differentiate between drug-induced acute kidney injury (AKI) and drug interactions through the renal transporter, it has been examined whether these transporter inhibitions quantitatively explained increases in serum creatinine (SCr) at their clinically relevant concentrations using drugs without any changes in renal function. For such renal transporter inhibitors and recently approved tyrosine kinase inhibitors (TKIs), this mini-review describes clinical increases in SCr and inhibitory potentials against the renal transporters. Most cases of SCr elevations can be explained by considering the renal transporter inhibitions based on unbound maximum plasma concentrations, except for drugs associated with obvious changes in renal function. SCr increases for cobicistat, dolutegravir, and dronedarone, and some TKIs were significantly underestimated, and these underestimations were suggested to be associated with low plasma unbound fractions. Sensitivity analysis of SCr elevations regarding inhibitory potentials of MATE1/2-K demonstrated that typical inhibitors such as cimetidine, DX-619, pyrimethamine, and trimethoprim could give false interpretations of AKI according to the criteria based on relative or absolute levels of SCr elevations. Recent progress and current challenges of physiologically-based pharmacokinetics modeling for creatinine disposition were also summarized. Although it should be noted for the potential impact of in vitro assay designs on clinical translatability of transporter inhibitions data, mechanistic approaches could support decision-making in clinical development to differentiate between AKI and creatinine-drug interactions. SIGNIFICANCE STATEMENT: Serum creatinine (SCr) is widely used as an indicator of kidney function, but it increases due to inhibitions of renal transporters, such as multidrug and toxin extrusion protein 1/2-K despite no functional changes in the kidney. Such SCr elevations were quantitatively explained by renal transporter inhibitions except for some drugs with high protein binding. The present analysis demonstrated that clinically relevant inhibitors of the renal transporters could cause SCr elevations above levels corresponding to acute kidney injury criteria.

Emerging Roles of Uremic Toxins and Inflammatory Cytokines in the Alteration of Hepatic Drug Disposition in Patients with Kidney Dysfunction

Patients with kidney dysfunction exhibit distinct pharmacokinetic profiles compared to those with normal kidney function. Hence, it is desirable to monitor the drug efficacy and toxicity caused by fluctuations in plasma drug concentrations associated with kidney dysfunction. Recently, pharmacokinetic information of drugs excreted mainly through the urine of patients with kidney dysfunction has been reported via drug-labeling information. Pharmacokinetic changes in drugs mainly eliminated by the liver cannot be overlooked as drug metabolism and/or transport activity in the liver may also be altered in patients with kidney dysfunction; however, the underlying mechanisms remain unclear. To plan an appropriate dosage regimen, it is necessary to clarify the underlying processes of functional changes in pharmacokinetic proteins. In recent years, uremic toxins have been shown to reduce the activity and/or expression of renal and hepatic transporters. This inhibitory effect has been reported to be time-dependent. In addition, inflammatory cytokines, such as interleukin-6, released from immune cells activated by uremic toxins and/or kidney injury can reduce the expression levels of drug-metabolizing enzymes and transporters in human hepatocytes. In this mini-review, we have summarized the renal and hepatic pharmacokinetic changes as well as the potential underlying mechanisms in kidney dysfunction, such as the chronic kidney disease and acute kidney injury. SIGNIFICANCE STATEMENT: Patients with kidney dysfunction exhibit distinct pharmacokinetic profiles compared to those with normal kidney function. Increased plasma concentrations of uremic toxins and inflammatory cytokines during kidney disease may potentially affect the activities and/or expression levels of drug-metabolizing enzymes and transporters in the liver and kidneys.

Coproporphyrin I as an Endogenous Biomarker to Detect Reduced OATP1B Activity and Shift in Elimination Route in Chronic Kidney Disease

Coproporphyrin I (CPI) is an endogenous biomarker of organic anion transporting polypeptide 1B transporter (OATP1B). CPI plasma baseline was reported to increase with severity of chronic kidney disease (CKD). Further, ratio of CPI area under the plasma concentration-time curve (AUCR) in the presence/absence of OATP1B inhibitor rifampin was higher in patients with CKD compared with healthy participants, in contrast to pitavastatin (a clinical OATP1B probe). This study investigated mechanism(s) contributing to altered CPI baseline in patients with CKD by extending a previously developed physiologically-based pharmacokinetic (PBPK) model to this patient population. CKD-related covariates were evaluated in a stepwise manner on CPI fraction unbound in plasma (f_u,p ), OATP1B-mediated hepatic uptake clearance (CL_active ), renal clearance (CL_R ), and endogenous synthesis (k_syn ). The CPI model successfully recovered increased baseline and rifampin-mediated AUCR in patients with CKD by accounting for the following disease-related changes: 13% increase in f_u,p , 29% and 39% decrease in CL_active in mild and moderate to severe CKD, respectively, decrease in CL_R proportional to decline in glomerular filtration rate, and 27% decrease in k_syn in severe CKD. Almost complete decline in CPI renal elimination in severe CKD increased its fraction transported by OATP1B, rationalizing differences in the CPI-rifampin interaction observed between healthy participants and patients with CKD. In conclusion, mechanistic modeling performed here supports CKD-related decrease in OATP1B function to inform prospective PBPK modeling of OATP1B-mediated drug-drug interaction in these patients. Monitoring of CPI allows detection of CKD-drug interaction risk for OATP1B drugs with combined hepatic and renal elimination which may be underestimated by extrapolating the interaction risk based on pitavastatin data in healthy participants.

Effect of Chronic Kidney Disease on the Renal Secretion via Organic Anion Transporters 1/3: Implications for Physiologically-Based Pharmacokinetic Modeling and Dose Adjustment

There is growing evidence that active tubular secretory clearance (CL_s) may not decline proportionally with the glomerular filtration rate (GFR) in chronic kidney disease (CKD), leading to the overestimation of renal clearance (CL_r) when using solely GFR to approximate disease effect on renal elimination. The clinical pharmacokinetic data of 33 renally secreted OAT1/3 substrates were collated to investigate the impact of mild, moderate, and severe CKD on CL_r, tubular secretion and protein binding (f_u,p). The f_u,p of the collated substrates ranged from 0.0026 to 1.0 in healthy populations; observed CKD‐related increase in the f_u,p (up to 2.7‐fold) of 8 highly bound substrates (f_u,p ≤ 0.2) was accounted for in the analysis. Use of prediction equation based on disease‐related changes in albumin resulted in underprediction of the CKD‐related increase in f_u,p of highly bound substrates, highlighting the necessity to measure protein binding in severe CKD. The critical analysis of clinical data for 33 OAT1/3 probes established that decrease in OAT1/3 activity proportional to the changes in GFR was insufficient to recapitulate effects of severe CKD on unbound tubular secretion clearance. OAT1/3‐mediated CL_s was estimated to decline by an additional 50% relative to the GFR decline in severe CKD, whereas change in active secretion in mild and moderate CKD was proportional to GFR. Consideration of this additional 50% decline in OAT1/3‐mediated CL_s is recommended for physiologically‐based pharmacokinetic models and dose adjustment of OAT1/3 substrates in severe CKD, especially for substrates with high contribution of the active secretion to CL_r.

Clinical Implications of Altered Drug Transporter Abundance/Function and PBPK Modeling in Specific Populations: An ITC Perspective

The role of membrane transporters on pharmacokinetics (PKs), drug-drug interactions (DDIs), pharmacodynamics (PDs), and toxicity of drugs has been broadly recognized. However, our knowledge of modulation of transporter expression and/or function in the diseased patient population or specific populations, such as pediatrics or pregnancy, is still emerging. This white paper highlights recent advances in studying the changes in transporter expression and activity in various diseases (i.e., renal and hepatic impairment and cancer) and some specific populations (i.e., pediatrics and pregnancy) with the focus on clinical implications. Proposed alterations in transporter abundance and/or activity in diseased and specific populations are based on (i) quantitative transporter proteomic data and relative abundance in specific populations vs. healthy adults, (ii) clinical PKs, and emerging transporter biomarker and/or pharmacogenomic data, and (iii) physiologically-based pharmacokinetic modeling and simulation. The potential for altered PK, PD, and toxicity in these populations needs to be considered for drugs and their active metabolites in which transporter-mediated uptake/efflux is a major contributor to their absorption, distribution, and elimination pathways and/or associated DDI risk. In addition to best practices, this white paper discusses current challenges and knowledge gaps to study and quantitatively predict the effects of modulation in transporter activity in these populations, together with the perspectives from the International Transporter Consortium (ITC) on future directions.

Renal Transporter-Mediated Drug-Biomarker Interactions of the Endogenous Substrates Creatinine and N1 -Methylnicotinamide: A PBPK Modeling Approach

Endogenous biomarkers for transporter-mediated drug-drug interaction (DDI) predictions represent a promising approach to facilitate and improve conventional DDI investigations in clinical studies. This approach requires high sensitivity and specificity of biomarkers for the targets of interest (e.g., transport proteins), as well as rigorous characterization of their kinetics, which can be accomplished utilizing physiologically-based pharmacokinetic (PBPK) modeling. Therefore, the objective of this study was to develop PBPK models of the endogenous organic cation transporter (OCT)2 and multidrug and toxin extrusion protein (MATE)1 substrates creatinine and N¹ -methylnicotinamide (NMN). Additionally, this study aimed to predict kinetic changes of the biomarkers during administration of the OCT2 and MATE1 perpetrator drugs trimethoprim, pyrimethamine, and cimetidine. Whole-body PBPK models of creatinine and NMN were developed utilizing studies investigating creatinine or NMN exogenous administration and endogenous synthesis. The newly developed models accurately describe and predict observed plasma concentration-time profiles and urinary excretion of both biomarkers. Subsequently, models were coupled to the previously built and evaluated perpetrator models of trimethoprim, pyrimethamine, and cimetidine for interaction predictions. Increased creatinine plasma concentrations and decreased urinary excretion during the drug-biomarker interactions with trimethoprim, pyrimethamine, and cimetidine were well-described. An additional inhibition of NMN synthesis by trimethoprim and pyrimethamine was hypothesized, improving NMN plasma and urine interaction predictions. To summarize, whole-body PBPK models of creatinine and NMN were built and evaluated to better assess creatinine and NMN kinetics while uncovering knowledge gaps for future research. The models can support investigations of renal transporter-mediated DDIs during drug development.

Imeglimin population pharmacokinetics and dose adjustment predictions for renal impairment in Japanese and Western patients with type 2 diabetes

Imeglimin is an orally administered first-in-class drug to treat type 2 diabetes mellitus (T2DM) and is mainly excreted unchanged by the kidneys. The present study aimed to define the pharmacokinetic (PK) characteristics of imeglimin using population PK analysis and to determine the optimal dosing regimen for Japanese patients with T2DM and chronic kidney disease (CKD). Imeglimin plasma concentrations in Japanese and Western healthy volunteers, and patients with T2DM, including patients with mild to severe CKD with an estimated glomerular filtration rate (eGFR) > 14 mL/min/1.73 m² were included in a population PK analysis. Pharmacokinetic simulations were conducted using a population PK model, and the area under concentration-time curve (AUC) was extrapolated with power regression analysis to lower eGFR. The influence of eGFR, weight, and age on apparent clearance and of dose on relative bioavailability were quantified by population PK analysis. Simulations and extrapolation revealed that the recommended dosing regimen based on the AUC was 500 mg twice daily (b.i.d.) for patients with eGFR 15-45 mL/min/1.73 m² , and 500 mg with a longer dosing interval was suggested for those with eGFR < 15. Simulations revealed that differences in plasma AUCs between Japanese and Western patients at the same dose were mainly driven by a difference in the eGFR and that the plasma AUC after 1,000 and 1,500 mg b.i.d. in Japanese and Western patients, respectively, was comparable in the phase IIb studies. These results indicate suitable dosages of imeglimin in the clinical setting of T2DM with renal impairment.

PBPK Simulation-Based Evaluation of Ganciclovir Crystalluria Risk Factors: Effect of Renal Impairment, Old Age, and Low Fluid Intake

Dosing guidance is often lacking for chronic kidney disease (CKD) due to exclusion of such patients from pivotal clinical trials. Physiologically based pharmacokinetic (PBPK) modelling supports model-informed dosing when clinical data are lacking, but application of these approaches to patients with impaired renal function is not yet at full maturity. In the current study, a ganciclovir PBPK model was developed for patients with normal renal function and extended to CKD population. CKD-related changes in tubular secretion were explored in the mechanistic kidney model and implemented either as proportional or non-proportional decline relative to GFR. Crystalluria risk was evaluated in different clinical settings (old age, severe CKD and low fluid intake) by simulating ganciclovir medullary collecting duct (MCD) concentrations. The ganciclovir PBPK model captured observed changes in systemic pharmacokinetic endpoints in mild-to-severe CKD; these trends were evident irrespective of assumed pathophysiological mechanism of altered active tubular secretion in the model. Minimal difference in simulated ganciclovir MCD concentrations was noted between young adult and geriatric populations with normal renal function and urine flow (1 mL/min), with lower concentrations predicted for severe CKD patients. High crystalluria risk was identified at reduced urine flow (0.1 mL/min) as simulated ganciclovir MCD concentrations exceeded its solubility (2.6–6 mg/mL), irrespective of underlying renal function. The analysis highlighted the importance of appropriate distribution of virtual subjects’ systems data in CKD populations. The ganciclovir PBPK model illustrates the ability of this translational tool to explore individual and combined effects of age, urine flow, and renal impairment on local drug renal exposure.

Quantitative Proteomic Map of Enzymes and Transporters in the Human Kidney: Stepping Closer to Mechanistic Kidney Models to Define Local Kinetics

The applications of translational modeling of local drug concentrations in various organs had a sharp increase over the last decade. These are part of the model-informed drug development initiative, adopted by the pharmaceutical industry and promoted by drug regulatory agencies. With respect to the kidney, the models serve as a bridge for understanding animal vs. human observations related to renal drug disposition and any consequential adverse effects. However, quantitative data on key drug-metabolizing enzymes and transporters relevant for predicting renal drug disposition are limited. Using targeted and global quantitative proteomics, we determined the abundance of multiple enzymes and transporters in 20 human kidney cortex samples. Nine enzymes and 22 transporters were quantified (8 for the first time in the kidneys). In addition, > 4,000 proteins were identified and used to form an open database. CYP2B6, CYP3A5, and CYP4F2 showed comparable, but generally low expression, whereas UGT1A9 and UGT2B7 levels were the highest. Significant correlation between abundance and activity (measured by mycophenolic acid clearance) was observed for UGT1A9 (Rs = 0.65, P = 0.004) and UGT2B7 (Rs = 0.70, P = 0.023). Expression of P-gp ≈ MATE-1 and OATP4C1 transporters were high. Strong intercorrelations were observed between several transporters (P-gp/MRP4, MRP2/OAT3, and OAT3/OAT4); no correlation in expression was apparent for functionally related transporters (OCT2/MATEs). This study extends our knowledge of pharmacologically relevant proteins in the kidney cortex, with implications on more prudent use of mechanistic kidney models under the general framework of quantitative systems pharmacology and toxicology.

Reduced physiologically-based pharmacokinetic model of dabigatran etexilate-dabigatran and its application for prediction of intestinal P-gp-mediated drug-drug interactions

BACKGROUND

Dabigatran etexilate (DABE) has been suggested as a clinical probe for intestinal P-glycoprotein (P-gp)-mediated drug-drug interaction (DDI) studies and, as an alternative to digoxin. Clinical DDI data with various P-gp inhibitors demonstrated a dose-dependent inhibition of P-gp with DABE. The aims of this study were to develop a joint DABE (prodrug)-dabigatran reduced physiologically-based-pharmacokinetic (PBPK) model and to evaluate its ability to predict differences in P-gp DDI magnitude between a microdose and a therapeutic dose of DABE.

METHODS

A joint DABE-dabigatran PBPK model was developed with a mechanistic intestinal model accounting for the regional P-gp distribution in the gastrointestinal tract. Model input parameters were estimated using DABE and dabigatran pharmacokinetic (PK) clinical data obtained after administration of DABE alone or with a strong P-gp inhibitor, itraconazole, and over a wide range of DABE doses (from 375 µg to 400 mg). Subsequently, the model was used to predict extent of DDI with additional P-gp inhibitors and with different DABE doses.

RESULTS

The reduced DABE-dabigatran PBPK model successfully described plasma concentrations of both prodrug and metabolite following administration of DABE at different dose levels and when co-administered with itraconazole. The model was able to capture the dose dependency in P-gp mediated DDI. Predicted magnitude of itraconazole P-gp DDI was higher at the microdose (predicted vs. observed median fold-increase in AUC_+inh/AUC_control (min-max) = 5.88 (4.29-7.93) vs. 6.92 (4.96-9.66) ) compared to the therapeutic dose (predicted median fold-increase in AUC_+inh/AUC_control = 3.48 (2.37-4.84) ). In addition, the reduced DABE-dabigatran PBPK model predicted successfully the extent of DDI with verapamil and clarithromycin as P-gp inhibitors. Model-based simulations of dose staggering predicted the maximum inhibition of P-gp when DABE microdose was concomitantly administered with itraconazole solution; simulations also highlighted dosing intervals required to minimise the DDI risk depending on the DABE dose administered (microdose vs. therapeutic).

CONCLUSIONS

This study provides a modelling framework for the evaluation of P-gp inhibitory potential of new molecular entities using DABE as a clinical probe. Simulations of dose staggering and regional differences in the extent of intestinal P-gp inhibition for DABE microdose and therapeutic dose provide model-based guidance for design of prospective clinical P-gp DDI studies.

Population pharmacokinetic modeling and simulation to support qualification of pyridoxic acid as endogenous biomarker of OAT1/3 renal transporters

Renal clearance of many drugs is mediated by renal organic anion transporters OAT1/3 and inhibition of these transporters may lead to drug‐drug interactions (DDIs). Pyridoxic acid (PDA) and homovanillic acid (HVA) were indicated as potential biomarkers of OAT1/3. The objective of this study was to develop a population pharmacokinetic model for PDA and HVA to support biomarker qualification. Simultaneous fitting of biomarker plasma and urine data in the presence and absence of potent OAT1/3 inhibitor (probenecid, 500 mg every 6 h) was performed. The impact of study design (multiple vs. single dose of OAT1/3 inhibitor) and ability to detect interactions in the presence of weak/moderate OAT1/3 inhibitors was investigated, together with corresponding power calculations. The population models developed successfully described biomarker baseline and PDA/HVA OAT1/3‐mediated interaction data. No prominent effect of circadian rhythm on PDA and HVA individual baseline levels was evident. Renal elimination contributed greater than 80% to total clearance of both endogenous biomarkers investigated. Estimated probenecid unbound in vivo OAT inhibitory constant was up to 6.4‐fold lower than in vitro values obtained with PDA as a probe. The PDA model was successfully verified against independent literature reported datasets. No significant difference in power of DDI detection was found between multiple and single dose study design when using the same total daily dose of 2000 mg probenecid. Model‐based simulations and power calculations confirmed sensitivity and robustness of plasma PDA data to identify weak, moderate, and strong OAT1/3 inhibitors in an adequately powered clinical study to support optimal design of prospective clinical OAT1/3 interaction studies.