Discovery of a NAPE-PLD inhibitor that modulates emotional behavior in mice.

N-acylethanolamines (NAEs), which include the endocannabinoid anandamide, represent an important family of signaling lipids in the brain. The lack of chemical probes that modulate NAE biosynthesis in living systems hamper the understanding of the biological role of these lipids. Using a high-throughput screen, chemical proteomics and targeted lipidomics, we report here the discovery and characterization of LEI-401 as a CNS-active N-acylphosphatidylethanolamine phospholipase D (NAPE-PLD) inhibitor. LEI-401 reduced NAE levels in neuroblastoma cells and in the brain of freely moving mice, but not in NAPE-PLD KO cells and mice, respectively. LEI-401 activated the hypothalamus-pituitary-adrenal axis and impaired fear extinction, thereby emulating the effect of a cannabinoid CB1 receptor antagonist, which could be reversed by a fatty acid amide hydrolase inhibitor. Our findings highlight the distinctive role of NAPE-PLD in NAE biosynthesis in the brain and suggest the presence of an endogenous NAE tone controlling emotional behavior.

Resurrecting the phoenix: When an assay fails.

Understanding protein-small-molecule interactions is a critical component of rational drug-design. Structure-activity relationship (SAR)-guided medicinal chemistry is informed by the biological outcome, as assessed by biochemical activity or cellular effect, of chemical modifications on small molecules. The effectiveness of SAR is reliant on the sturdiness and durability of assay design and the quality of information garnered from assays. Lack of quality data at this step can lead to obstruction of the drug discovery pipeline with profound implications for the timelines of introducing a drug into the market. Hence, it would not be an overstatement to consider biochemical/biological assays as the backbone of drug-discovery. Enzyme assays can fail for many different reasons, with the enzyme and the substrate being the principal players. Lack of clarity can hamper progress and can lead to mounting costs and potentially losing competitive advantage. Although each assay is unique and requires a specific approach to troubleshoot the problem at hand, there are general guidelines that can be followed to maximize the chances of success. This review is a step-by-step attempt at reintroducing fundamental biochemical concepts within the context of an enzyme assay, delineating probable causes for failure and potential approaches to get an assay back up and running.

Mechanistic and Quantitative Understanding of Pharmacokinetics in Zebrafish Larvae through Nanoscale Blood Sampling and Metabolite Modeling of Paracetamol.

Zebrafish larvae are increasingly used for pharmacological research, but internal drug exposure is often not measured. Understanding pharmacokinetics is necessary for reliable translation of pharmacological results to higher vertebrates, including humans. Quantification of drug clearance and distribution requires measurements of blood concentrations. Additionally, measuring drug metabolites is of importance to understand clearance in this model organism mechanistically. We therefore mechanistically studied and quantified pharmacokinetics in zebrafish larvae, and compared this to higher vertebrates, using paracetamol (acetaminophen) as a paradigm compound. A method was developed to sample blood from zebrafish larvae 5 days post fertilization. Blood concentrations of paracetamol and its major metabolites, paracetamol-glucuronide and paracetamol-sulfate, were measured. Blood concentration data were combined with measured amounts in larval homogenates and excreted amounts and simultaneously analyzed through nonlinear mixed-effects modeling, quantifying absolute clearance and distribution volume. Blood sampling from zebrafish larvae was most successful from the posterior cardinal vein, with a median volume (interquartile range) of 1.12 nl (0.676-1.66 nl) per blood sample. Samples were pooled (n = 15-35) to reach measurable levels. Paracetamol blood concentrations at steady state were only 10% of the external paracetamol concentration. Paracetamol-sulfate was the major metabolite, and its formation was quantified using a time-dependent metabolic formation rate. Absolute clearance and distribution volume correlated well with reported values in higher vertebrates, including humans. Based on blood concentrations and advanced data analysis, the mechanistic and quantitative understanding of paracetamol pharmacokinetics in zebrafish larvae has been established. This will improve the translational value of this vertebrate model organism in drug discovery and development. SIGNIFICANCE STATEMENT: In early phases of drug development, new compounds are increasingly screened in zebrafish larvae, but the internal drug exposure is often not taken into consideration. We developed innovative experimental and computational methods, including a blood-sampling technique, to measure the paradigm drug paracetamol (acetaminophen) and its major metabolites and quantify pharmacokinetics (absorption, distribution, elimination) in zebrafish larvae of 5 days post fertilization with a total volume of only 300 nl. These parameter values were scaled to higher vertebrates, including humans.

Quantitative profiling of endocannabinoids and related N-acylethanolamines in human CSF using nano LC-MS/MS.

Endocannabinoids, a class of lipid messengers, have emerged as crucial regulators of synaptic communication in the CNS. Dysregulation of these compounds has been implicated in many brain disorders. Although some studies have identified and quantified a limited number of target compounds, a method that provides comprehensive quantitative information on endocannabinoids and related N-acylethanolamines (NAEs) in cerebrospinal fluid (CSF) is currently lacking, as measurements are challenging due to low concentrations under normal physiological conditions. Here we developed and validated a high-throughput nano LC-ESI-MS/MS platform for the simultaneous quantification of endocannabinoids (anandamide and 2-arachidonoylglycerol), ten related NAEs, and eight additional putatively annotated NAEs in human CSF. Requiring only 200 µl of CSF, our method has limits of detection from 0.28 to 61.2 pM with precisions of relative SD <15% for most compounds. We applied our method to CSF from 45 healthy humans and demonstrated potential age and gender effects on concentrations of endocannabinoids and NAEs. Notably, our results show that docosahexaenoylethanolamide concentrations increase with age in males. Our method may offer new opportunities to gain insight into regulatory functions of endocannabinoids in the context of (ab)normal brain function.

Highly Selective, Reversible Inhibitor Identified by Comparative Chemoproteomics Modulates Diacylglycerol Lipase Activity in Neurons.

Diacylglycerol lipase (DAGL)-α and -ß are enzymes responsible for the biosynthesis of the endocannabinoid 2-arachidonoylglycerol (2-AG). Selective and reversible inhibitors are required to study the function of DAGLs in neuronal cells in an acute and temporal fashion, but they are currently lacking. Here, we describe the identification of a highly selective DAGL inhibitor using structure-guided and a chemoproteomics strategy to characterize the selectivity of the inhibitor in complex proteomes. Key to the success of this approach is the use of comparative and competitive activity-based proteome profiling (ABPP), in which broad-spectrum and tailor-made activity-based probes are combined to report on the inhibition of a protein family in its native environment. Competitive ABPP with broad-spectrum fluorophosphonate-based probes and specific ß-lactone-based probes led to the discovery of α-ketoheterocycle LEI105 as a potent, highly selective, and reversible dual DAGL-α/DAGL-ß inhibitor. LEI105 did not affect other enzymes involved in endocannabinoid metabolism including abhydrolase domain-containing protein 6, abhydrolase domain-containing protein 12, monoacylglycerol lipase, and fatty acid amide hydrolase and did not display affinity for the cannabinoid CB1 receptor. Targeted lipidomics revealed that LEI105 concentration-dependently reduced 2-AG levels, but not anandamide levels, in Neuro2A cells. We show that cannabinoid CB1-receptor-mediated short-term synaptic plasticity in a mouse hippocampal slice model can be reduced by LEI105. Thus, we have developed a highly selective DAGL inhibitor and provide new pharmacological evidence to support the hypothesis that "on demand biosynthesis" of 2-AG is responsible for retrograde signaling.

Development of a Series of Pyrrolopyridone MAT2A Inhibitors.

The optimization of an allosteric fragment, discovered by differential scanning fluorimetry, to an in vivo MAT2a tool inhibitor is discussed. The structure-based drug discovery approach, aided by relative binding free energy calculations, resulted in AZ'9567 (21), a potent inhibitor in vitro with excellent preclinical pharmacokinetic properties. This tool showed a selective antiproliferative effect on methylthioadenosine phosphorylase (MTAP) KO cells, both in vitro and in vivo, providing further evidence to support the utility of MAT2a inhibitors as potential anticancer therapies for MTAP-deficient tumors.

Discovery, Optimization, and Biological Evaluation of Arylpyridones as Cbl-b Inhibitors.

Casitas B-lymphoma proto-oncogene-b (Cbl-b), a member of the Cbl family of RING finger E3 ubiquitin ligases, has been demonstrated to play a central role in regulating effector T-cell function. Multiple studies using gene-targeting approaches have provided direct evidence that Cbl-b negatively regulates T, B, and NK cell activation via a ubiquitin-mediated protein modulation. Thus, inhibition of Cbl-b ligase activity can lead to immune activation and has therapeutic potential in immuno-oncology. Herein, we describe the discovery and optimization of an arylpyridone series as Cbl-b inhibitors by structure-based drug discovery to afford compound 31. This compound binds to Cbl-b with an IC50 value of 30 nM and induces IL-2 production in T-cells with an EC50 value of 230 nM. Compound 31 also shows robust intracellular target engagement demonstrated through inhibition of Cbl-b autoubiquitination, inhibition of ubiquitin transfer to ZAP70, and the cellular modulation of phosphorylation of a downstream signal within the TCR axis.

Accelerating the Validation of Endogenous On-Target Engagement and In Cellulo Kinetic Assessment for Covalent Inhibitors of KRASG12C in Early Drug Discovery.

Covalent inhibition is a valuable modality in drug discovery because of its potential ability in decoupling pharmacokinetics from pharmacodynamics by prolonging the residence time of the drug on the target of interest. This increase in target occupancy is limited only by the rate of target turnover. However, a limitation in such studies is to translate the in vitro inhibition assessment to the appropriate in cellulo target engagement parameter by covalent probes. Estimation of such parameters is often impeded by the low-throughput nature of current probe-free approaches. In this study, an ultra-performance liquid chromatography-multiple reaction monitoring mass spectrometry platform was utilized to develop a targeted proteomics workflow that can evaluate cellular on-target engagement of covalent molecules in an increased throughput manner. This workflow enabled a throughput increase of 5-10 fold when compared to traditional nanoLC-based proteomics studies. To demonstrate the applicability of the method, KRASG12C was used as a model system to investigate the interaction of an irreversible covalent small molecule, compound 25, both in vitro and in cellulo. Initial biochemical studies confirmed that the small molecule forms an adduct with the targeted cysteine on the protein, as assessed at the level of both intact protein and on the target peptide. In cellulo studies were carried out to quantify target engagement and allele selectivity assessment for the small molecule in the heterozygous NCI-H358 cell line for KRASG12C with respect to the WT type protein. The workflow enabled evaluation of in vitro and in cellulo target engagement kinetics, providing mechanistic insights into the irreversible mode of inhibition. In summary, the method has the potential for target agnostic application in the assessment of on-target engagement of covalent probes compatible with the high-throughput requirements of early drug discovery.

Structure-Activity Relationship Studies of α-Ketoamides as Inhibitors of the Phospholipase A and Acyltransferase Enzyme Family.

The phospholipase A and acyltransferase (PLAAT) family of cysteine hydrolases consists of five members, which are involved in the Ca2+-independent production of N-acylphosphatidylethanolamines (NAPEs). NAPEs are lipid precursors for bioactive N-acylethanolamines (NAEs) that are involved in various physiological processes such as food intake, pain, inflammation, stress, and anxiety. Recently, we identified α-ketoamides as the first pan-active PLAAT inhibitor scaffold that reduced arachidonic acid levels in PLAAT3-overexpressing U2OS cells and in HepG2 cells. Here, we report the structure-activity relationships of the α-ketoamide series using activity-based protein profiling. This led to the identification of LEI-301, a nanomolar potent inhibitor for the PLAAT family members. LEI-301 reduced the NAE levels, including anandamide, in cells overexpressing PLAAT2 or PLAAT5. Collectively, LEI-301 may help to dissect the physiological role of the PLAATs.

Impact of post-hatching maturation on the pharmacokinetics of paracetamol in zebrafish larvae.

Zebrafish larvae are increasingly used in pharmacological and toxicological studies, but it is often overlooked that internal exposure to exogenous compounds, rather than the incubation medium concentration, is driving observed effects. Moreover, as the zebrafish larva is a developing organism, continuous physiological changes impact pharmacokinetic or toxicokinetic processes like the absorption and elimination of exogenous compounds, influencing the interpretation of observations and conclusions drawn from experiments at different larval ages. Here, using paracetamol as paradigm compound, mathematical modelling is used to quantify absorption and elimination rates from internal exposure over time profiles after waterborne treatment, as well as changes in these parameters in post-hatching larvae of 3, 4, and 5 days post fertilisation (dpf). An increase of 106% in absorption rate was observed between 3 and 4 dpf, but no further increase at 5 dpf, and an increase of 17.5% in elimination rate for each dpf. Paracetamol clearance, determined from elimination rate constants and reported total larval volumes of 253, 263, and 300 nL at 3, 4, and 5 dpf respectively, correlates best with higher vertebrates at 5 dpf. This suggests that when studying direct effects of exogenous compounds, experiments with zebrafish larvae are best performed at 5 dpf.

Identification of α,ß-Hydrolase Domain Containing Protein 6 as a Diacylglycerol Lipase in Neuro-2a Cells.

The endocannabinoid 2-arachidonoylglycerol (2-AG) is involved in neuronal differentiation. This study aimed to identify the biosynthetic enzymes responsible for 2-AG production during retinoic acid (RA)-induced neurite outgrowth of Neuro-2a cells. First, we confirmed that RA stimulation of Neuro-2a cells increases 2-AG production and neurite outgrowth. The diacylglycerol lipase (DAGL) inhibitor DH376 blocked 2-AG production and reduced neuronal differentiation. Surprisingly, CRISPR/Cas9-mediated knockdown of DAGLα and DAGLß in Neuro-2a cells did not reduce 2-AG levels, suggesting another enzyme capable of producing 2-AG in this cell line. Chemical proteomics revealed DAGLß and α,ß-hydrolase domain containing protein (ABHD6) as the only targets of DH376 in Neuro-2a cells. Biochemical, genetic and lipidomic studies demonstrated that ABHD6 possesses DAGL activity in conjunction with its previously reported monoacylglycerol lipase activity. RA treatment of Neuro-2a cells increased by three-fold the amount of active ABHD6. Our study shows that ABHD6 exhibits significant DAG lipase activity in Neuro-2a cells in addition to its known MAG lipase activity and suggest it is involved in neuronal differentiation.

Activity-Based Protein Profiling Identifies α-Ketoamides as Inhibitors for Phospholipase A2 Group XVI.

Phospholipase A2, group XVI (PLA2G16) is a thiol hydrolase from the HRASLS family that regulates lipolysis in adipose tissue and has been identified as a host factor enabling the cellular entry of picornaviruses. Chemical tools are essential to visualize and control PLA2G16 activity, but they have not been reported to date. Here, we show that MB064, which is a fluorescent lipase probe, also labels recombinant and endogenously expressed PLA2G16. Competitive activity-based protein profiling (ABPP) using MB064 enabled the discovery of α-ketoamides as the first selective PLA2G16 inhibitors. LEI110 was identified as a potent PLA2G16 inhibitor ( Ki = 20 nM) that reduces cellular arachidonic acid levels and oleic acid-induced lipolysis in human HepG2 cells. Gel-based ABPP and chemical proteomics showed that LEI110 is a selective pan-inhibitor of the HRASLS family of thiol hydrolases (i.e., PLA2G16, HRASLS2, RARRES3 and iNAT). Molecular dynamic simulations of LEI110 in the reported crystal structure of PLA2G16 provided insight in the potential ligand-protein interactions to explain its binding mode. In conclusion, we have developed the first selective inhibitor that can be used to study the cellular role of PLA2G16.

Gene Expression of Endocannabinoid System Components in Skeletal Muscle and Adipose Tissue of South Asians and White Caucasians with Overweight.

OBJECTIVE: The study aimed to investigate whether markers of endocannabinoid signaling differed between men with overweight of South Asian and white Caucasian descent. METHODS: We included South Asian (n = 10) and white Caucasian (n = 10) men with overweight and prediabetes aged 35 to 50 years. Plasma samples were analyzed for endocannabinoids, their congeners, and lipids. In white adipose tissue (WAT) and skeletal muscle biopsies, mRNA expression of genes involved in the endocannabinoid system (ECS) was assessed using quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Fasting lipid oxidation and glucose oxidation were determined with indirect calorimetry. RESULTS: Compared to white Caucasians, South Asians had higher levels of plasma 2-linoleoyl glycerol (P < 0.01) and N-linoleoylethanolamine (P < 0.05). Interestingly, in skeletal muscle of South Asians, expression of cannabinoid receptors CB1 and CB2 was 10-fold lower (P < 0.001) and that of the endocannabinoid degradation enzyme fatty acid amide hydrolase 2 (FAAH2) was 5-fold lower (P < 0.001) compared to white Caucasians. Expression of genes involved in the ECS in WAT were not different between the two ethnicities. After pooling of both ethnicities, plasma 2-arachidonoylglycerol (2-AG) positively correlated with plasma triglycerides (R = 0.77, P < 0.001) and lipid oxidation (R = 0.55, P < 0.05). CONCLUSIONS: South Asian men with overweight have higher plasma 2-linoleoyl glycerol and N-linoleoylethanolamine levels and lower expression of CB receptors and the endocannabinoid degradation enzyme FAAH2 in skeletal muscle compared to white Caucasians.

Caloric restriction lowers endocannabinoid tonus and improves cardiac function in type 2 diabetes.

BACKGROUND/OBJECTIVES: Endocannabinoids (ECs) are associated with obesity and ectopic fat accumulation, both of which play a role in the development of cardiovascular disease (CVD) in type 2 diabetes (T2D). The effect of prolonged caloric restriction on ECs in relation to fat distribution and cardiac function is still unknown. Therefore, our aim was to investigate this relationship in obese T2D patients with coronary artery disease (CAD). SUBJECTS/METHODS: In a prospective intervention study, obese T2D patients with CAD (n = 27) followed a 16 week very low calorie diet (VLCD; 450-1000 kcal/day). Cardiac function and fat accumulation were assessed with MRI and spectroscopy. Plasma levels of lipid species, including ECs, were measured using liquid chromatography-mass spectrometry. RESULTS: VLCD decreased plasma levels of virtually all measured lipid species of the class of N-acylethanolamines including the EC anandamide (AEA; -15%, p = 0.016), without decreasing monoacylglycerols including the EC 2-arachidonoylglycerol (2-AG). Baseline plasma AEA levels strongly correlated with the volume of subcutaneous white adipose tissue (SAT; R2 = 0.44, p < 0.001). VLCD decreased the volume of SAT (-53%, p < 0.001), visceral white adipose tissue (VAT) (-52%, p < 0.001), epicardial white adipose tissue (-15%, p < 0.001) and paracardial white adipose tissue (-28%, p < 0.001). VLCD also decreased hepatic (-86%, p < 0.001) and myocardial (-33%, p < 0.001) fat content. These effects were accompanied by an increased left ventricular ejection fraction (54.8 ± 8.7-56.2 ± 7.9%, p = 0.016). CONCLUSIONS: Caloric restriction in T2D patients with CAD decreases AEA levels, but not 2-AG levels, which is paralleled by decreased lipid accumulation in adipose tissue, liver and heart, and improved cardiovascular function. Interestingly, baseline AEA levels strongly correlated with SAT volume. We anticipate that dietary interventions are worthwhile strategies in advanced T2D, and that reduction in AEA may contribute to the improved cardiometabolic phenotype induced by weight loss.

High Fat Diet Increases Circulating Endocannabinoids Accompanied by Increased Synthesis Enzymes in Adipose Tissue.

The endocannabinoid system (ECS) controls energy balance by regulating both energy intake and energy expenditure. Endocannabinoid levels are elevated in obesity suggesting a potential causal relationship. This study aimed to elucidate the rate of dysregulation of the ECS, and the metabolic organs involved, in diet-induced obesity. Eight groups of age-matched male C57Bl/6J mice were randomized to receive a chow diet (control) or receive a high fat diet (HFD, 45% of calories derived from fat) ranging from 1 day up to 18 weeks before euthanasia. Plasma levels of the endocannabinoids 2-arachidonoylglycerol (2-AG) and anandamide (N-arachidonoylethanolamine, AEA), and related N-acylethanolamines, were quantified by UPLC-MS/MS and gene expression of components of the ECS was determined in liver, muscle, white adipose tissue (WAT) and brown adipose tissue (BAT) during the course of diet-induced obesity development. HFD feeding gradually increased 2-AG (+132% within 4 weeks, P < 0.05), accompanied by upregulated expression of its synthesizing enzymes Daglα and ß in WAT and BAT. HFD also rapidly increased AEA (+81% within 1 week, P < 0.01), accompanied by increased expression of its synthesizing enzyme Nape-pld, specifically in BAT. Interestingly, Nape-pld expression in BAT correlated with plasma AEA levels (R 2 = 0.171, ß = 0.276, P < 0.001). We conclude that a HFD rapidly activates adipose tissue depots to increase the synthesis pathways of endocannabinoids that may aggravate the development of HFD-induced obesity.

Integration of pharmacometabolomics with pharmacokinetics and pharmacodynamics: towards personalized drug therapy.

Personalized medicine, in modern drug therapy, aims at a tailored drug treatment accounting for inter-individual variations in drug pharmacology to treat individuals effectively and safely. The inter-individual variability in drug response upon drug administration is caused by the interplay between drug pharmacology and the patients' (patho)physiological status. Individual variations in (patho)physiological status may result from genetic polymorphisms, environmental factors (including current/past treatments), demographic characteristics, and disease related factors. Identification and quantification of predictors of inter-individual variability in drug pharmacology is necessary to achieve personalized medicine. Here, we highlight the potential of pharmacometabolomics in prospectively informing on the inter-individual differences in drug pharmacology, including both pharmacokinetic (PK) and pharmacodynamic (PD) processes, and thereby guiding drug selection and drug dosing. This review focusses on the pharmacometabolomics studies that have additional value on top of the conventional covariates in predicting drug PK. Additionally, employing pharmacometabolomics to predict drug PD is highlighted, and we suggest not only considering the endogenous metabolites as static variables but to include also drug dose and temporal changes in drug concentration in these studies. Although there are many endogenous metabolite biomarkers identified to predict PK and more often to predict PD, validation of these biomarkers in terms of specificity, sensitivity, reproducibility and clinical relevance is highly important. Furthermore, the application of these identified biomarkers in routine clinical practice deserves notable attention to truly personalize drug treatment in the near future.

Endocannabinoid tone is higher in healthy lean South Asian than white Caucasian men.

South Asians have a higher risk to develop obesity and related disorders compared to white Caucasians. This is likely in part due to their lower resting energy expenditure (REE) as related with less energy-combusting brown adipose tissue (BAT). Since overactivation of the endocannabinoid system is associated with obesity and low BAT activity, we hypothesized that South Asians have a higher endocannabinoid tone. Healthy lean white Caucasian (n = 10) and South Asian (n = 10) men were cold-exposed to activate BAT. Before and after cooling, REE was assessed and plasma was collected for analysis of endocannabinoids and lipids. At thermoneutrality, South Asians had higher plasma levels of 2-arachidonoylglycerol (2-AG; 11.36 vs 8.19 pmol/mL, p < 0.05), N-arachidonylethanolamine (AEA; 1.04 vs 0.89 pmol/mL, p = 0.05) and arachidonic acid (AA; 23.24 vs 18.22 nmol/mL, p < 0.001). After pooling of both ethnicities, plasma 2-AG but not AEA positively correlated with triglycerides (R2 = 0.32, p < 0.05) and body fat percentage (R2 = 0.18, p < 0.05). Interestingly, AA negative correlated with REE (R2 = 0.46, p < 0.001) and positively with body fat percentage (R2 = 0.33, p < 0.01). Cooling increased endocannabinoids. In conclusion, South Asian compared to white Caucasian men have higher endocannabinoid tone. This suggests that endocannabinoids may, at least in part, underlie the disadvantageous metabolic phenotype of South Asians later in life.

Activity-based protein profiling reveals off-target proteins of the FAAH inhibitor BIA 10-2474.

A recent phase 1 trial of the fatty acid amide hydrolase (FAAH) inhibitor BIA 10-2474 led to the death of one volunteer and produced mild-to-severe neurological symptoms in four others. Although the cause of the clinical neurotoxicity is unknown, it has been postulated, given the clinical safety profile of other tested FAAH inhibitors, that off-target activities of BIA 10-2474 may have played a role. Here we use activity-based proteomic methods to determine the protein interaction landscape of BIA 10-2474 in human cells and tissues. This analysis revealed that the drug inhibits several lipases that are not targeted by PF04457845, a highly selective and clinically tested FAAH inhibitor. BIA 10-2474, but not PF04457845, produced substantial alterations in lipid networks in human cortical neurons, suggesting that promiscuous lipase inhibitors have the potential to cause metabolic dysregulation in the nervous system.

Pharmacokinetic Modeling of Paracetamol Uptake and Clearance in Zebrafish Larvae: Expanding the Allometric Scale in Vertebrates with Five Orders of Magnitude.

Zebrafish larvae (Danio rerio) are increasingly used to translate findings regarding drug efficacy and safety from in vitro-based assays to vertebrate species, including humans. However, the limited understanding of drug exposure in this species hampers its implementation in translational research. Using paracetamol as a paradigm compound, we present a novel method to characterize pharmacokinetic processes in zebrafish larvae, by combining sensitive bioanalytical methods and nonlinear mixed effects modeling. The developed method allowed quantification of paracetamol and its two major metabolites, paracetamol-sulfate and paracetamol-glucuronide in pooled samples of five lysed zebrafish larvae of 3 days post-fertilization. Paracetamol drug uptake was quantified to be 0.289 pmole/min and paracetamol clearance was quantified to be 1.7% of the total value of the larvae. With an average volume determined to be 0.290 µL, this yields an absolute clearance of 2.96 × 107 L/h, which scales reasonably well with clearance rates in higher vertebrates. The developed methodology will improve the success rate of drug screens in zebrafish larvae and the translation potential of findings, by allowing the establishment of accurate exposure profiles and thereby also the establishment of concentration-effect relationships.