# A Transversal Approach Combining In Silico, In Vitro and In Vivo Models to Describe the Metabolism of the Receptor Interacting Protein 1 Kinase Inhibitor Sibiriline

Romain Pelletier, Thomas Gicquel, Mélanie Simoes Eugenio, Pierre-Jean Ferron, Isabelle Morel, Claire Delehouzé, Marie-Thérèse Dimanche-Boitrel, Morgane Rousselot, Brendan Le Daré

► To cite this version:

Romain Pelletier, Thomas Gicquel, Mélanie Simoes Eugenio, Pierre-Jean Ferron, Isabelle Morel, et al. A Transversal Approach Combining In Silico, In Vitro and In Vivo Models to Describe the Metabolism of the Receptor Interacting Protein 1 Kinase Inhibitor Sibiriline. Pharmaceutics, 2022, 14 (12), pp.2665. 10.3390/pharmaceutics14122665. [hal-03930266](https://hal.science/hal-03930266)

HAL Id: hal-03930266

https://hal.science/hal-03930266

Submitted on 9 Jan 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# A Transversal Approach Combining In Silico, In Vitro and In Vivo Models to Describe the Metabolism of the Receptor Interacting Protein 1 Kinase Inhibitor Sibiriline

Romain Pelletier <fnref n="1" />,<fnref n="2" /> ID, Thomas Gicquel <fnref n="1" />,<fnref n="2" /> ID, Mélanie Simoes Eugenio <fnref n="3" />,<fnref n="4" />, Pierre-Jean Ferron <fnref n="1" />, Isabelle Morel <fnref n="1" />,<fnref n="2" />, Claire Delehouzé <fnref n="4" />, Marie-Thérèse Dimanche-Boitrel <fnref n="3" />,<fnref n="4" />, Morgane Rousselot <fnref n="4" /> ID and Brendan Le Daré <fnref n="1" />,<fnref n="5" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Pelletier, R.; Gicquel, T.; Simoes Eugenio, M.; Ferron, P.-J.; Morel, I.; Delehouzé, C.; Dimanche-Boitrel, M.-T.; Rousselot, M.; Le Daré, B. A Transversal Approach Combining In Silico, In Vitro and In Vivo Models to Describe the Metabolism of the Receptor Interacting Protein 1 Kinase Inhibitor Sibiriline. *Pharmaceutics* **2022**, *14*, 2665. https://doi.org/10.3390/pharmaceutics14122665

Academic Editor: Im-Sook Song

Received: 17 October 2022

Accepted: 26 November 2022

Published: 30 November 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Sibiriline is a novel drug inhibiting receptor-interacting protein 1 kinase (RIPK1) and necroptosis, a regulated form of cell death involved in several disease models. In this study, we aimed to investigate the metabolic fate of sibiriline in a cross-sectional manner using an in silico prediction, coupled with in vitro and in vivo experiments. In silico predictions were performed using GLORYx and Biotransformer 3.0 freeware; in vitro incubation was performed on differentiated human HepaRG cells, and in vivo experiments including a pharmacokinetic study were performed on mice treated with sibiriline. HepaRG culture supernatants and mice plasma samples were analyzed with ultra-high-performance liquid chromatography, coupled with tandem mass spectrometry (LC-HRMS/MS). The molecular networking bioinformatics tool applied to LC-HRMS/MS data allowed us to visualize the sibiriline metabolism kinetics. Overall, 14 metabolites, mostly produced by Phase II transformations (glucuronidation and sulfation) were identified. These data provide initial reassurance regarding the toxicology of this new RIPK1 inhibitor, although further studies are required.

**Keywords:** molecular networking; Sibiriline; metabolism; LC-HRMS/MS

## 1. Introduction

Receptor-interacting protein 1 kinase (RIPK1) inhibitors are emerging as a new pharmaceutical class proposed for inflammatory, autoimmune, and neurodegenerative diseases such as rheumatoid arthritis, ulcerative colitis, psoriasis, Alzeihmer disease, amyotrophic lateral sclerosis, or multiple sclerosis [1-4]. Among them, Sibiriline appeared as a novel drug inhibiting necroptosis (a regulated form of cell death involved in several disease models). Sibiriline was shown to protect liver functions on a RIPK1-dependent concanavalin A-induced hepatitis, positioning it as a potential drug candidate for the treatment of immune-dependent hepatitis [5]. However, the therapeutic potential of Sibiriline will depend on its risk/benefit ratio and pharmacokinetic profile, including knowledge of its metabolism and potential drug-drug interactions. To date, nothing is known about the metabolism of Sibiriline, and its comprehensive exploration would require the use of cross-sectional analytical approaches.

In vitro, numerous models exist to study metabolism, including acellular models such as human microsomes and cell-based models such as the HepaRG line. Differentiated HepaRG cells (a bipotent cell line capable of differentiating into either cholangiocyte- or hepatocyte-like cells under specific culture conditions) have been described as a robust model to explore both toxicity and metabolism of xenobiotics [6-9]. In particular, this

<footnote n="1">¹ Institut NuMeCan (Nutrition, Metabolism and Cancer), INSERM, INRAE, CHU Rennes, PREVITOX Network, 35033 Rennes, France</footnote>

<footnote n="2">² Clinical and Forensic Toxicology Laboratory, Rennes University Hospital, Univ Rennes, 35033 Rennes, France</footnote>

<footnote n="3">³ IRSET (Institut de Recherche en Santé, Environnement et Travail)—UMR_S 1085, Univ Rennes, Inserm, EHESP, 35000 Rennes, France</footnote>

<footnote n="4">⁴ SeaBeLife SAS, Station Biologique, Place Georges Teissier, 29680 Roscoff, France</footnote>

<footnote n="5">⁵ Pharmacy Department, Rennes University Hospital, 35033 Rennes, France</footnote>

<footnote>* Correspondence: brendan.le.dare@chu-rennes.fr</footnote>
model expresses the majority of the drug processing genes including major CYPs [10], SULTs [11], and UGTs [11,12] and has long-term functional stability whereas the human hepatocytes in primary culture lose their differentiated phenotype and drug metabolism over time. Moreover, HepaRG cells have shown greater metabolite production compared with human liver microsomes motivating its growing use [13].

To complete the in vitro data, several in vivo models are also described and allow the addition of some pharmacokinetic notions to the metabolism study [14,15]. In particular, the rodent model is currently widely used to study the metabolism of xenobiotics for which we do not have human in vivo samples, as in the case of new psychoactive products [16,17]. While the quantitative aspect of rodent metabolism cannot be extrapolated to humans, the qualitative identification of metabolites is particularly relevant.

In order to reprocess complex analytical data collected in vitro and in vivo, innovative bioinformatic approaches have recently emerged. This is notably the case of molecular networking (MN) which allows the organization and representation of untargeted tandem mass spectrometry (MS/MS) data in a graphical form [18]. Each node represents an ion and its associated fragmentation spectrum, the links between the nodes indicating similarities between spectra. By allowing structural information to be propagated through the network and facilitating sample-to-sample comparison, the MN approach provides valuable insights into drug metabolism [19]. Consequently, a semi-quantitative visualization of molecule repartition in different matrix samples can be provided by a multi-matrix approach. MN has already proven its interest in drug metabolism analysis for in vitro, in vivo clinical, or forensic purposes [9,19–24].

Going further, the knowledge evolution of metabolism mechanisms has allowed the development of in silico metabolism prediction algorithms, which are a convenient, open access, time-saving, and inexpensive tool to broaden the metabolite search and confirm data found in vivo or in vitro [25]. Several approaches can be used in metabolism studies to produce in silico systems and include models based on (i) the quantitative structure-activity relationship (QSAR), which assumes that molecules with similar structures potentially have similar metabolic properties, or (ii) quantum mechanical calculations used to predict reactivity and/or (iii) docking of potential substrates into the active site of the enzyme [26–30].

In the present study, we aimed to investigate the metabolic fate of Sibiriline in a cross-sectional manner using an in silico prediction, coupled with in vitro and in vivo experiments.

## 2. Materials and Methods

### 2.1. *In Silico* Metabolite Prediction

Sibiriline putative metabolites were predicted using Sibiriline SMILES string (OC(C=C1)=CC=C1C2NC3=NC=CC=C3C2) through the free online software GLORYx [31] and Biotransformer 3.0 [32]. Briefly, GLORYx integrates machine learning-based site of metabolism prediction with reaction rule sets to predict and classify the putative structures of metabolites that could be formed by phase I and/or phase II metabolism. Biotransformer 3.0 is an open-access software tool that supports the rapid, accurate, and comprehensive prediction of the metabolism of small molecules in both mammals and environmental microorganisms. All metabolite structures were generated using ChemDraw software (PerkinElmer, Inc., Waltham, MA, USA).

### 2.2. Pharmacokinetic Study in Mice

Mice treatment. The in vivo pharmacokinetic study of Sibiriline was performed with male Swiss mice. A dose of 5 mg/kg Sibiriline was administered to 18 mice intraperitoneally (IP). After the times of 15 min, 30 min, 1 h, 2 h, 3 h, and 4 h, the mice were sacrificed. Blood was collected and plasma was separated by centrifugation. The samples were frozen and stored at -80 °C before processing.

Samples extraction. Each mouse plasma sample (400 µL) was mixed with 1 mL acetonitrile to precipitate the proteins and extract the compound. The samples were then
vortexed for 5 min and placed in an ultrasonic bath for 1 min. The precipitated proteins were sedimented by centrifugation (15,000× g, 5 min at 16 °C) and the supernatants were transferred to a microplate for analysis by LC-MS/MS.

Analytical method for the quantitation of Sibiriline in mice. All samples are analysed with a validated Ultra-High Performance Liquid Chromatography method coupled to a Shimadzu LC-MS 8030 triple quadrupole. Solutions of a standard range, prepared with plasma from untreated mice, were analysed in the same series of injections.

LC separation was performed using a gradient on a C18 Phenomenex Kinetex column (50 × 2.1 mm, 2.6 µm). The mobile phase consisted of water 0.5% formic acid (solvent A) and acetonitrile (solvent B) delivered at a flow rate of 0.5 mL/min using the following stepwise gradient elution program: initial conditions of 95:5 (A:B) run maintained from 0 to 1.2 min, run from 5:95 (A:B) at 1.2 min maintained from 1.2 to 1.4 min, run from 95:5 (A:B) at 1.4 min, maintained from 1.4 to 2.8 min, Injections of 2 µL into the LC-HR-MS system were performed in a thermostated column at 40 °C. The tuning parameters of mass spectrometry (MS) were optimized by direct infusion of Sibiriline at the concentration of 10 µg/mL in a 1:1 acetonitrile/water mixture into the ionization probe in positive ion mode. The detection parameters were optimized to increase the sensitivity and signal stability by infusion of Sibiriline. Multiple reaction monitoring transitions (i.e., 211.0 → 118.0 m/z and 211.0 → 183.1 m/z) were used respectively for Sibiriline quantitation and confirmation. These analytical parameters provided satisfactory separation and peak shapes of Sibiriline in less than 3 min, and retention times (RTs) for Sibiriline was 0.86 min.

## 2.3. *In Vitro and In Vivo Metabolism Studies*

Material. William's E medium (ref: 12551032) was purchased from Gibco (ThermoFischer Scientific, San Jose, CA, USA). Penicillin-streptomycin was obtained from Life Technologies (Grand Island, NY, USA). Fetal Bovine Serum (FBS) was purchased from Eurobio (Courtaboeuf, France) and from Hyclone GE Healthcare Life Sciences (Logan, UT, USA). Hydrocortisone hemisuccinate was purchased from Serb (Paris, France). Dimethyl sulfoxide (DMSO), formic acid, and insulin were obtained from Sigma-Aldrich (Saint Louis, MO, USA). Sibiriline was generously donated by SeaBeLife Biotech (Roscoff, France).

Cell culture and treatment. Progenitor HepaRG cells were cultured as already described [10]. Briefly, HepaRG cells were seeded at a density of 105 cells/well in 96-well plates and cultured for two weeks in culture medium (William's E medium (1X) supplemented with 10% FBS, 50 U/mL penicillin, 50 µg/mL streptomycin, 5 µg/mL insulin, 2 mM glutamine, 50 µM sodium hydrocortisone hemisuccinate and 2% DMSO). Cells were then cultured for two more weeks in the same medium supplemented with 2% DMSO to induce cell differentiation into cholangiocyte- and hepatocyte-like cells [6]. The detection of Sibiriline and its metabolites was performed using this model as already described [9,33]. Differentiated HepaRG cells were incubated with 100 µL of Sibiriline (10 µM) during H0, H8, H24, or H48.

Mice treatment. The metabolism study was performed in male C57BL/6J mice (Janvier Labs, Le Genest Saint Isle, France). A dose of 5 mg/kg Sibiriline was administered to 10 mice intraperitoneally (IP). After the times 1 h and 3 h, the mice were sacrificed. Blood was collected and plasma was separated by centrifugation. The samples were frozen and stored at -80 °C before processing. Mice were housed in cages in a conventional animal facility with a 12 h dark-light cycle with controlled temperature between 19 °C and 20 °C. Experiments were done in compliance with French laws and the institution's guidelines for animal welfare. Authors were authorized to conduct animal experimentation by "la Direction des Services Vétérinaires" (License M Samson #A3523840). The project was authorized by the "Comité Rennais d'Ethique en matière d'Experimentation Animal" [CREAA] and the license was given by the "Ministère de l'Enseignement Supérieur, de la Recherche et de l'Innovation", #32246-2021061616397414).

Samples extraction. In vitro samples (200 µL) and in vivo plasma samples (100 µL) obtained from mice at the Rennes University Hospital were extracted as previously de-
scribed [20,33]. Briefly, samples were supplemented with 500 µL of methanol containing internal standard (risperidone-D4) and then extracted with 300 µL of 0.1 M zinc sulfate solution. After supernatant evaporation, the residue was dissolved in 200 µL of LC-MS grade water and transferred into chromatographic vials for LC-HR-MS analysis.

LC-MS settings. Non-targeted screening LC-HRMS/MS method used for MN building was performed as already described [9,33]. Briefly, liquid chromatography was performed on an Accucore Phenyl Hexyl (100 × 2.1 mm, 2.6 µm) (Thermo Scientific, San Jose, CA, USA) at 40 °C using an elution gradient at a flow rate of 500 µL/min during 15 min with 10 µL as injection volume. An orbitrap mass spectrometer was operated in positive ESI mode and the acquisition range was 100–700 *m/z*. Ionic precursor selection was performed in a “data-dependent” mode of operation, where the 5 most intense ions from the previous scan were selected for fragmentation (Top N of 5).

MN generation. Spectral data allowed us to generate MN using a semi-quantitative approach. Data acquisition, processing (i.e., MS data conversion, preprocessing, MS1 annotation, and generation of molecular networks), visualization, and network analysis have already been described [19]. Briefly, raw data were converted to an open MS format (.mzXML) with ProteoWizard's MSConvert module [34]. The mzXML files were then preprocessed (deconvolution, de-isotoping, alignment, gap-filling) with MZmine 2 software [35]. The single mgf output file was then loaded on the Global Natural Products Social networking (GNPS) web-based platform in order to generate the multi-matrix molecular network [18]. For the use of high-resolution data, the basic parameters were modified to *m/z* 0.02 for the mass tolerance of precursor and fragment ions used for MS/MS spectral library searching, and *m/z* 0.02 for the mass tolerance of fragment ions used for MN. The minimum cluster size was set to 1. In addition, links between nodes were created when the cosine score was greater than 0.70, and the minimum number of common fragment ions shared by two MS/MS spectra was 6. Links between two nodes were only kept in the network if each node was in the top 10 most similar nodes.

The molecular network was visualized using Cytoscape 3.8.0 software [36]. The nodes were annotated by comparison with reference standards by spectral matching with the curated GNPS, mzCloud online mass spectral libraries, and information propagation [37].

Figure 1 shows the overall methodology of the cross-sectional metabolic study approach used in this study.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Overview methodology for cross-sectional study of sibiriline metabolism used in this study. LC-HRMS/MS = Liquid Chromatography coupled to a High-Resolution Tandem Mass Spectrometer.</figcaption></figure>
## 3. Results

### 3.1. *In Silico* Prediction of Sibiriline Metabolism

In silico prediction of Sibiriline metabolism was performed on GloryX and Biotransformer 3.0 freeware and results are shown in Table 1. Predicted Sibiriline derivatives with an assigned GloryX score were considered probable (P1 to P14), while derivatives without an assigned score were considered minority and less likely (P15 to P21). Overall, Phase I and Phase II metabolites were predicted. Major Phase I metabolites included hydroxylated (P2, P4, P5, P7, P13, P14, P16, P17, P18, P19, P20) and oxidated (P10, P12, P15) derivatives. Major Phase II metabolites included glucuronoconjugated (P1, P6, P9, P16, P17, P18, P19), sulfoconjugated (P3, P20) and methylated (P8, P21) derivatives.

<table><caption>**Table 1.** In silico predicted sibiriline metabolites and associated *m/z*, molecular structure, elemental composition, metabolic reaction, exact mass shift, predictive score (GLORYx software) and enzyme predicted to be involved in biotransformation (BT3: Biotransformer 3.0).</caption><thead><tr><th>Predicted Metabolite (m/z)<br>Elemental<br>Composition<br>Metabolite Reaction</th><th>Putative<br>Biotransformation<br>Reaction<br>Exact Mass Shift</th><th>Molecular Structure</th><th>Prediction<br>Software</th><th>Score and/or Enzyme<br>Involved</th></tr></thead><tbody><tr><td>Sibiriline (m/z 211.0871)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<br>None</td><td>Not applicable</td><td><img src="" alt="Sibiriline molecular structure"></td><td>Not applicable</td><td>Not applicable</td></tr><tr><td>P1 (m/z 387.1192)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>7</sub><br>O-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>6</sub><br>+176.032</td><td><img src="" alt="P1 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.948<br>Predicted (UDP-<br>glucuronosyltransferase)</td></tr><tr><td>P2 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation</td><td>+O<br>+15.994</td><td><img src="" alt="P2 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.656<br>Predicted (CYP1A2)</td></tr><tr><td>P3 (m/z 291.0439)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>4</sub>S<br>Sulfation</td><td>+SO<sub>3</sub><br>+79.956</td><td><img src="" alt="P3 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.628<br>Predicted (Aryl<br>sulfotransferase)</td></tr><tr><td>P4 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation</td><td>+O<br>+15.994</td><td><img src="" alt="P4 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.348<br>Not predicted</td></tr><tr><td>P5 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation (para<br>to carbon)</td><td>+O<br>+15.994</td><td><img src="" alt="P5 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.348<br>Not predicted</td></tr><tr><td>P6 (m/z 387.1192)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>7</sub><br>N-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>6</sub><br>+176.032</td><td><img src="" alt="P6 molecular structure"></td><td>GLORYx<br>BT3</td><td>0.308<br>Not predicted</td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Predicted Metabolite (m/z)<br>Elemental<br>Composition<br>Metabolite Reaction</th><th>Putative<br>Biotransformation<br>Reaction<br>Exact Mass Shift</th><th>Molecular Structure</th><th>Prediction<br>Software</th><th>Score and/or Enzyme<br>Involved</th></tr></thead><tbody><tr><td>P7 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation (para<br>to nitrogen)</td><td>+O<br>+15.994</td><td><img src="" alt="Chemical structure of p7 (m/z 227.0820)"></td><td>GLORYx<br>BT3</td><td>0.295<br>Predicted (CYP2C9)</td></tr><tr><td>P8 (m/z 225.1027)<br>C<sub>14</sub>H<sub>12</sub>N<sub>2</sub>O<br>Methylation (aromatic OH)</td><td>+CH<sub>2</sub><br>+14.015</td><td><img src="" alt="Chemical structure of p8 (m/z 225.1027)"></td><td>GLORYx<br>BT3</td><td>0.28<br>Not predicted</td></tr><tr><td>P9 (m/z 387.1192)<br>C<sub>19</sub>H<sub>19</sub>N<sub>2</sub>O<sub>7</sub>+<br>N-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>6</sub><br>+176.032</td><td><img src="" alt="Chemical structure of p9 (m/z 387.1192)"></td><td>GLORYx<br>BT3</td><td>0.224<br>Not predicted</td></tr><tr><td>P10 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>N-oxidation</td><td>+O<br>+15.994</td><td><img src="" alt="Chemical structure of p10 (m/z 227.0820)"></td><td>GLORYx<br>BT3</td><td>0.212<br>Predicted (CYP2C8)</td></tr><tr><td>P11 (m/z 119.0609)<br>C<sub>7</sub>H<sub>6</sub>N<sub>2</sub></td><td>Unknown</td><td><img src="" alt="Chemical structure of p11 (m/z 119.0609)"></td><td>GLORYx<br>BT3</td><td>0.1312<br>Not predicted</td></tr><tr><td>P12 (m/z 109.0289)<br>C<sub>6</sub>H<sub>4</sub>O<sub>2</sub><br>Oxidation of 4-substituted<br>phenol to quinone</td><td>Unknown</td><td><img src="" alt="Chemical structure of p12 (m/z 109.0289)"></td><td>GLORYx<br>BT3</td><td>0.1312<br>Not predicted</td></tr><tr><td>P13 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation</td><td>+O<br>+15.994</td><td><img src="" alt="Chemical structure of p13 (m/z 227.0820)"></td><td>GLORYx<br>BT3</td><td>0.072<br>Predicted (CYP2C9)</td></tr><tr><td>P14 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Aromatic hydroxylation</td><td>+O<br>+15.994</td><td><img src="" alt="Chemical structure of p14 (m/z 227.0820)"></td><td>GLORYx<br>BT3</td><td>0.028<br>Predicted (CYP1A2)</td></tr><tr><td>P15 (m/z 227.0820)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub><br>Oxidation</td><td>+O<br>+15.994</td><td><img src="" alt="Chemical structure of p15 (m/z 227.0820)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (CYP2E1)</td></tr><tr><td>P16 (m/z 403.1141)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>8</sub><br>Hydroxylation +<br>O-glucuronidation</td><td>C<sub>6</sub>H<sub>8</sub>O<sub>7</sub><br>+192.027</td><td><img src="" alt="Chemical structure of p16 (m/z 403.1141)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (UDP-glucuronosyltransferase)</td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Predicted Metabolite (m/z)<br>Elemental<br>Composition<br>Metabolite Reaction</th><th>Putative<br>Biotransformation<br>Reaction<br>Exact Mass Shift</th><th>Molecular Structure</th><th>Prediction<br>Software</th><th>Score and/or Enzyme<br>Involved</th></tr></thead><tbody><tr><td>P17 (m/z 403.1141)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>8</sub><br>Hydroxylation +<br>O-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>7</sub><br>+192.027</td><td><img src="" alt="Chemical structure of P17 (m/z 403.1141)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (UDP-<br>glucuronosyltransferase)</td></tr><tr><td>P18 (m/z 403.1141)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>8</sub><br>Hydroxylation +<br>O-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>7</sub><br>+192.027</td><td><img src="" alt="Chemical structure of P18 (m/z 403.1141)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (UDP-<br>glucuronosyltransferase)</td></tr><tr><td>P19 (m/z 403.1141)<br>C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>8</sub><br>Hydroxylation +<br>O-glucuronidation</td><td>+C<sub>6</sub>H<sub>8</sub>O<sub>7</sub><br>+192.027</td><td><img src="" alt="Chemical structure of P19 (m/z 403.1141)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (UDP-<br>glucuronosyltransferase)</td></tr><tr><td>P20 (m/z 307.0388)<br>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>5</sub>S<br>Hydroxylation + O-sulfation</td><td>+SO<sub>4</sub><br>+95.951</td><td><img src="" alt="Chemical structure of P20 (m/z 307.0388)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (Aryl<br>sulfotransferase)</td></tr><tr><td>P21 (m/z 241.0977)<br>C<sub>14</sub>H<sub>12</sub>N<sub>2</sub>O<sub>2</sub><br>O-methylation</td><td>+CH<sub>2</sub><br>+14.015</td><td><img src="" alt="Chemical structure of P21 (m/z 241.0977)"></td><td>GLORYx<br>BT3</td><td>None<br>Predicted (Catechol<br>O-methyltransferase)</td></tr></tbody></table>

### 3.2. *In Vitro Sibiriline Metabolism*

To assess in vitro Sibiriline metabolite production, differentiated human HepaRG cells were incubated with Sibiriline for 24 h. Analysis of culture media by LC-HRMS/MS allowed us to generate a molecular network. Nodes are labelled with the exact protonated mass (m/z) and the links are labelled with the exact mass shift. Nodes were linked together in cluster according to their MS² spectral similarities (Figure 2).

In the sibiriline-containing cluster (Figure 2B), a specific color was assigned to each type of biotransformation reaction (glucuronidation in blue, sulfation in green, oxidation + sulfation in black, oxidation + glucuronidation in yellow, 2 × oxidation + 2 × glucuronidation in orange, methylation + sulfation reaction in pink and methylation + glucuroconjugation in salmon). Sibiriline was associated with 18 structurally related molecules. Using information propagation, a node identified using a spectral library can be used as a starting point to identify another node in the same cluster. The aim is to determine the
structure of an unknown molecule using structural information from known neighboring nodes. Especially, spectrally related molecules may have mass shifts corresponding to well-known biotransformation reactions as shown in Table 1 [20].

<figure><img src="image_6.png" /><figcaption>**Figure 2.** In vitro Sibiriline metabolism. Differentiated HepaRG were incubated with Sibiriline 10 µM during 24 h and culture supernatants were analyzed using LC-HR/MSMS. (A) The overall molecular network. (B) In this sibiriline-containing cluster, each type of biotransformation reaction was depicted in different colors (glucuronidation in bleu, sulfation in green, oxidation + sulfation in black, oxidation + glucuronidation in yellow, 2 × oxidation + 2 × glucuronidation in orange, methylation + sulfation reaction in pink and methylation + glucuroconjugation in salmon). Nodes are labelled with the exact protonated mass (m/z), retention time (RT) in minute and the links are labelled with the exact mass shift.</figcaption></figure>

Figure 2 shows two compounds with the same exact mass (m/z 387.116) but with two different retention times (RT: 4.8 and 3.7 min). These nodes are linked to Sibiriline (m/z 211.087) with a mass shift of +176.032, corresponding to glucuronidation reaction. From these nodes, we found three compounds with m/z 403.113 that could correspond to a hydroxylated derivative of Sibiriline glucuronide (RT: 4.1, 4.3, and 5 min) due to the mass delta of +15.994. We also found mass shifts that could correspond to sulfation (+79.956), giving rise to two m/z 291.043 compounds (RT: 5.2 and 4.4 min) (Figure 2) as well as mass shifts that could correspond to oxidation + sulfation (+95.951), giving rise to two m/z 307.038 compounds (RT: 5.2 and 4.4 min). Due to the mass shift of +28.031 that might correspond to a
double methylation, the *m/z* 415.150 and *m/z* 319.075 nodes could be methylated derivatives of Sibiriline glucuronide and Sibiriline sulfate, respectively. Lastly, we identified a mass shift of 192.060 that could correspond to a hydroxylation + glucuronidation of *m/z* 403.113 compound (sibiriline + 2 × hydroxylation + 2 × glucuronidation). Structure and annotation proposals were made based on the *m/z* found during in silico prediction. Taken together, these results suggest that these Sibiriline spectrally related molecules are putative Sibiriline metabolites.

### 3.3. Sibiriline Pharmacokinetic Study in Mice

To determine how Sibiriline is metabolized in vivo, a pharmacokinetic study was performed by treating Swiss mice with an IP injection of sibiriline (5 mg/kg) (Figure 3). Sibiriline rapidly reached a peak plasma concentration within 15 min. The half-life was estimated at 21 min. The volume of distribution at 3247 mL/kg is indicative of a hydrophilic molecule (Figure 3).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** In vivo pharmacokinetics study of sibiriline. Evolution of the plasma concentration of sibiriline after intraperitoneal administration of sibiriline (5 mg/kg) in Swiss mice.</figcaption></figure>

### 3.4. *In Vivo Sibiriline Metabolism*

In order to strengthen the confidence level of the potential metabolites that we found in vitro, we performed an MN approach following LC-HRMS/MS on plasma samples from mice treated with Sibiriline (5 mg/kg IP) (Figure 4). As Sibiriline is rapidly metabolised in mice (half-life: 21 min), we tested for the presence of metabolites at 1 h and 3 h after the Sibiriline IP injection. Analysis of plasma samples at different times after Sibiriline treatment allowed us to generate a corresponding multi-matrix molecular network. A specific color was assigned to each time (H1 in salmon, H3 in red) so that we can visualize the fate of molecules structurally close to Sibiriline. The different colored areas in each node represent concentrations of the corresponding compound in each condition in a semi-quantitative manner.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** In vivo Sibiriline metabolism kinetic. Mice plasma sample collected 1 h or 3 h after intraperitoneal injection of Sibiriline (5 mg/kg) were analyzed using LC-HRMS/MS. In this sibiriline-containing cluster, each type of biotransformation reaction was depicted in different colors (glucuronidation in bleu, sulfation in green, oxidation + glucuronidation in yellow, and oxidation/hydroxylation in grey). Nodes are labelled with the exact protonated mass (m/z), putative chemical structures and the links are labelled with the exact mass shift. RT = retention time in minute.</figcaption></figure>

Sibiriline-containing cluster visual analysis revealed that all putative metabolites found in the mice plasma samples were found in our previous in vitro experiment (Sibiriline glucuronide; Sibiriline sulfate, hydroxy-sibiriline glucuronide and hydroxy-sibiriline sulfate), except *m/z* 227.081 compounds which were not found in vitro. Due to the mass shift of +15.994, these compounds could be hydroxylated derivatives of Sibiriline, as predicted in silico. Interestingly, we found that almost all Sibiriline structurally related molecules are in the highest concentrations at 1 h of incubation, compared with 3 h of incubation. Only *m/z* 403.114 (RT: 4.3 min), which was found outside the sibiriline-containing cluster, showed higher concentrations at 3 h, compared with 1 h. These results show that Sibiriline and its structurally related compounds disappear rapidly over time.

Taken together, these results suggest that the Sibiriline structurally related compounds found in vitro and in vivo using molecular networking are putative metabolites. Also, we show that differentiated HepaRG is a powerful model in metabolism studies. Table 2 reports putative Sibiriline metabolites found in silico, in vitro and in vivo and provides annotations from M1 to M8, with the letters a, b, and c denoting isomers.
<table><caption>**Table 2.** Putative identified compounds or metabolites contained in the Sibiriline molecular network. RT: retention time; BT: Biotransformer 3; X = presence; - = absence.</caption><thead><tr><th rowspan="2">Metabolite Proposal</th><th rowspan="2">Formula Proposal</th><th rowspan="2">Biotransformation</th><th rowspan="2">m/z</th><th rowspan="2">Mass Shift</th><th rowspan="2">RT (min)</th><th rowspan="2">In Vivo</th><th rowspan="2">In Vitro</th><th colspan="2">In Silico</th></tr><tr><th>GLORYx</th><th>BT3</th></tr></thead><tbody><tr><td>-</td><td>C<sub>6</sub>H<sub>4</sub>O<sub>2</sub></td><td>Unknown</td><td>109.029</td><td>-102.058</td><td>-</td><td>-</td><td>-</td><td>X</td><td>-</td></tr><tr><td>-</td><td>C<sub>7</sub>H<sub>6</sub>N<sub>2</sub></td><td>Unknown</td><td>119.061</td><td>-92.026</td><td>-</td><td>-</td><td>-</td><td>X</td><td>-</td></tr><tr><td>-</td><td>C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O</td><td>Parent molecule (sibiriline)</td><td>211.087</td><td>0</td><td>5.8</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>M1a</td><td rowspan="2">C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>2</sub></td><td rowspan="2">Hydroxylation</td><td rowspan="2">227.082</td><td rowspan="2">15.995</td><td>5.0</td><td>X</td><td>-</td><td rowspan="2">X</td><td rowspan="2">X</td></tr><tr><td>M1b</td><td>5.2</td><td>X</td><td>-</td></tr><tr><td>M2a</td><td rowspan="2">C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>4</sub>S</td><td rowspan="2">Sulfation</td><td rowspan="2">291.044</td><td rowspan="2">79.957</td><td>4.4</td><td>-</td><td>X</td><td rowspan="2">X</td><td rowspan="2">X</td></tr><tr><td>M2b</td><td>5.2</td><td>X</td><td>X</td></tr><tr><td>M3a</td><td rowspan="2">C<sub>13</sub>H<sub>10</sub>N<sub>2</sub>O<sub>5</sub>S</td><td rowspan="2">Hydroxylation + sulfation</td><td rowspan="2">307.038</td><td rowspan="2">95.952</td><td>4.4</td><td>-</td><td>X</td><td rowspan="2">-</td><td rowspan="2">X</td></tr><tr><td>M3b</td><td>5.2</td><td>X</td><td>X</td></tr><tr><td>M4</td><td>C<sub>15</sub>H<sub>14</sub>N<sub>2</sub>O<sub>4</sub>S</td><td>2 × methylation + sulfation</td><td>319.075</td><td>107.988</td><td>5.1</td><td>-</td><td>X</td><td>-</td><td>-</td></tr><tr><td>M5a</td><td rowspan="2">C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>7</sub></td><td rowspan="2">Glucuronidation</td><td rowspan="2">387.119</td><td rowspan="2">176.033</td><td>3.7</td><td>X</td><td>X</td><td rowspan="2">X</td><td rowspan="2">X</td></tr><tr><td>M5b</td><td>4.8</td><td>X</td><td>X</td></tr><tr><td>M6a</td><td rowspan="3">C<sub>19</sub>H<sub>18</sub>N<sub>2</sub>O<sub>8</sub></td><td rowspan="3">Hydroxylation + glucuronidation</td><td rowspan="3">403.113</td><td rowspan="3">192.027</td><td>4.1</td><td>X</td><td>X</td><td rowspan="3">-</td><td rowspan="3">X</td></tr><tr><td>M6b</td><td>4.3</td><td>X</td><td>X</td></tr><tr><td>M6c</td><td>5.0</td><td>X</td><td>X</td></tr><tr><td>M7</td><td>C<sub>21</sub>H<sub>22</sub>N<sub>2</sub>O<sub>7</sub></td><td>2 × methylation + glucuronidation</td><td>415.150</td><td>204.063</td><td>4.3</td><td>-</td><td>X</td><td>-</td><td>-</td></tr><tr><td>M8</td><td>C<sub>25</sub>H<sub>26</sub>N<sub>2</sub>O<sub>15</sub></td><td>2 × hydroxylation + 2 × glucuronidation</td><td>595.183</td><td>384.096</td><td>5.9</td><td>-</td><td>X</td><td>-</td><td>-</td></tr></tbody></table>

## 4. Discussion

In this study, we aimed to investigate the metabolic fate of Sibiriline in a cross-sectional manner using an in silico prediction, coupled with in vitro and in vivo experiments.

Overall, this cross-sectional approach allowed us to identify 8 Sibiriline putative metabolites in mice plasma (in vivo) and 12 putative metabolites in a human in vitro model (M1 to M8, Table 2). The higher number of metabolites found in vitro compared with in vivo could probably be due to the accumulation of metabolites in the culture medium (unlike the mouse model which includes all pharmacokinetic steps). As we have already shown that the metabolites identified in the HepaRG supernatant were quite similar to those found in urine, a comparison with mouse urine would be interesting in this context [33]. Interestingly, the in silico predictions proposed the main compounds found in vitro/in vivo showing that this approach is suitable for Sibiriline and strengthening the confidence of our in vitro/in vivo results. Biotransformer 3 outperformed GLORYx in predicting Sibiriline metabolites based on in vivo/in vitro data. However, two inconsistent metabolites (m/z 109.029 and m/z 119.061) were predicted by the in silico approach, leading to caution in the interpretation of the results. Most significantly, these complementary approaches allow (i) to come closer to the comprehensiveness of a metabolic study and (ii) to increase the degree of confidence in the separately generated results.

In terms of biotransformation reactions, we showed in silico, in vivo, or in vitro that the metabolism of Sibiriline is mainly performed by phase II enzymes (m/z 291.044, m/z 307.038, m/z 319.075, m/z 387.119, m/z 403.113, m/z 415.149, m/z 595.183), despite phase I metabolites were also identified (m/z 227.082). These results constitute a first reassuring basis from a toxicology point of view, as phase II metabolism is widely regarded as detoxifying. These reactive metabolites can exert initial cellular stress, through many mechanisms, such as glutathione depletion, or binding to proteins, nucleic acid lipids, and other cellular structures. Toxicological studies on sibiriline will be interesting in this context.
Phase II metabolism produces very hydrophilic compounds and may explain the rapid elimination of conjugated compounds shown in vivo in this study. While commonly used laboratory animals such as mice have many anatomical, physiological, and biochemical similarities to humans, there are marked differences between the species in their ability to process a drug [38–41]. Therefore, the interpretation of the pharmacokinetic data in mice presented in this study must be weighed against these considerations and a pharmacokinetic study in humans would be interesting here to evaluate Sibiriline elimination. Also, future pharmacokinetic studies of repeated-dose Sibiriline should take into account the selected animal species because the pharmacokinetics of the molecules may differ between species [42]. In terms of metabolism and toxicity, further studies are also needed, due to the fact that a molecule may be toxic due to a particular metabolism in one species, but not toxic to other species [43].

In addition, this would be of great importance in apprehending the administration modalities if used for therapeutic purposes. A short half-life should necessitate the use of strategies to prolong the effect of the molecule and could include the use of continuous intravenous infusion [44], the use of liposomal forms [45], intramuscular [46] or subcutaneous [47] administration, or the use of sustained-release oral forms [48]. Also, the pharmacokinetic and pharmacodynamic parameters of each of these strategies should be able to be based on trough concentration or area under the curve to further optimize the administration modalities [49].

This study presents several limitations. First, pharmacokinetic data in mice should be interpreted with caution due to interspecies differences in drug processing as drug metabolism is generally faster in rodents than in humans. Also, the pharmacokinetic study was performed by IP which is not the route of administration in acute liver failure. Second, the analysis of metabolites in mouse plasma does not necessarily allow the visualisation of metabolites excreted in urine, and an analysis in this matrix could be of interest. Third, since we analysed the culture supernatant of differentiated HepaRGs in vitro, it cannot be excluded that reactive metabolites have conjugated to intracellular proteins, not allowing their detection. Lastly, since the data presented in this study are mostly qualitative, it is difficult to determine which are the main metabolites.

**5. Conclusions**

This cross-sectional approach to studying metabolism using in silico, in vitro, and in vivo models brings us closer to a comprehensive exploration of xenobiotic metabolism and is a valuable tool in the study of new drug candidates. In addition, we applied an original MN-based approach to make the best use of the complex analytical data acquired by LC-HRMS/MS. In this study, we propose 14 putative metabolites that have mainly undergone phase II biotransformations. Pharmacokinetic data in mice showed rapid elimination of Sibiriline and its metabolites. Overall, these data provide initial reassurance regarding the toxicology of this new RIPK-1 inhibitor, although further studies are required.

**Author Contributions:** Conceptualization, R.P., T.G., M.S.E., P.-J.F., M.-T.D.-B., M.R. and B.L.D.; methodology, R.P., T.G., M.S.E., P.-J.F., C.D., M.-T.D.-B., M.R. and B.L.D.; data curation, R.P., T.G., M.S.E., P.-J.F., M.-T.D.-B., M.R. and B.L.D.; writing—original draft preparation, R.P., T.G., B.L.D.; writing—review and editing, All authors; visualization, R.P., T.G., B.L.D.; supervision, B.L.D., T.G.; project administration, M.R., C.D.; funding acquisition, M.R., C.D. All authors have read and agreed to the published version of the manuscript.

**Funding:** In vivo experiments in mice and article processing charges were funded by SeaBeLife SAS (Roscoff, France).

**Institutional Review Board Statement:** All animal work was completed in accordance with institutional guidelines and approved by the French Ministère de l'Enseignement Supérieur, de la recherche et de l'innovation and local ethical committee (CREMEAS) according the APAFiS national procedure (agreement number #19987-2019032715282421, #31122-2021041916575097, #31303-202104271424385, #32246-2021061616397414 and #34292-202112082339929).
**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** We thank the ADME plateform from Institut Pasteur (Lille) & PCBIS—TechMedILL CNRS—UMS 3286 for the pharmacokinetic data and assays of sibiriline in mice.

**Conflicts of Interest:** The authors declare the following financial interests/personal relationships that can be considered as potential competing interests: Morgane Rousselot, Claire Delehouzé and Marie-Thérèse Dimanche-Boitrel are founders and hold shares in SeaBeLife, the company that produces sibiriline. Mélanie Simoes Eugenio is financed by a collaboration contract between Rennes 1 and SeaBeLife. The other authors declare that they have no conflicts of interest.

## References

1. Harris, P.A.; Berger, S.B.; Jeong, J.U.; Nagilla, R.; Bandyopadhyay, D.; Campobasso, N.; Capriotti, C.A.; Cox, J.A.; Dare, L.; Dong, X.; et al. Discovery of a First-in-Class Receptor Interacting Protein 1 (RIP1) Kinase Specific Clinical Candidate (GSK2982772) for the Treatment of Inflammatory Diseases. *J. Med. Chem.* **2017**, *60*, 1247–1261. [CrossRef]

2. Weisel, K.; Scott, N.E.; Tompson, D.J.; Votta, B.J.; Madhavan, S.; Povey, K.; Wolstenholme, A.; Simeoni, M.; Rudo, T.; Richards-Peterson, L.; et al. Randomized Clinical Study of Safety, Pharmacokinetics, and Pharmacodynamics of RIPK1 Inhibitor GSK2982772 in Healthy Volunteers. *Pharmacol. Res. Perspect.* **2017**, *5*, e00365. [CrossRef]

3. Sheridan, C. Death by Inflammation: Drug Makers Chase the Master Controller. *Nat. Biotechnol.* **2019**, *37*, 111–113. [CrossRef]

4. Martens, S.; Hofmans, S.; Declercq, W.; Augustyns, K.; Vandenabeele, P. Inhibitors Targeting RIPK1/RIPK3: Old and New Drugs. *Trends Pharmacol. Sci.* **2020**, *41*, 209–224. [CrossRef]

5. Le Cann, F.; Delehouzé, C.; Leverrier-Penna, S.; Filliol, A.; Comte, A.; Delalande, O.; Desban, N.; Baratte, B.; Gallais, I.; Piquet-Pellorce, C.; et al. Sibiriline, a New Small Chemical Inhibitor of Receptor-Interacting Protein Kinase 1, Prevents Immune-Dependent Hepatitis. *FEBS J.* **2017**, *284*, 3050–3068. [CrossRef]

6. Cerec, V.; Glaise, D.; Garnier, D.; Morosan, S.; Turlin, B.; Drenou, B.; Gripon, P.; Kremsdorf, D.; Guguen-Guillouzo, C.; Corlu, A. Transdifferentiation of Hepatocyte-like Cells from the Human Hepatoma HepaRG Cell Line through Bipotent Progenitor. *Hepatology* **2007**, *45*, 957–967. [CrossRef]

7. Guguen-Guillouzo, C.; Guillouzo, A. Setup and Use of HepaRG Cells in Cholestasis Research. In *Experimental Cholestasis Research*; Vinken, M., Ed.; Methods in Molecular Biology; Springer: New York, NY, USA, 2019; Volume 1981, pp. 291–312, ISBN 978-1-4939-9419-9.

8. Guillouzo, A.; Corlu, A.; Aninat, C.; Glaise, D.; Morel, F.; Guguen-Guillouzo, C. The Human Hepatoma HepaRG Cells: A Highly Differentiated Model for Studies of Liver Metabolism and Toxicity of Xenobiotics. *Chem. Biol. Interact.* **2007**, *168*, 66–73. [CrossRef]

9. Le Daré, B.; Ferron, P.-J.; Allard, P.-M.; Clément, B.; Morel, I.; Gicquel, T. New Insights into Quetiapine Metabolism Using Molecular Networking. *Sci. Rep.* **2020**, *10*, 19921. [CrossRef]

10. Aninat, C.; Piton, A.; Glaise, D.; Le Charpentier, T.; Langouët, S.; Morel, F.; Guguen-Guillouzo, C.; Guillouzo, A. Expression of Cytochromes P450, Conjugating Enzymes and Nuclear Receptors in Human Hepatoma HepaRG Cells. *Drug Metab. Dispos.* **2006**, *34*, 75–83. [CrossRef]

11. Hugbart, C.; Verres, Y.; Le Daré, B.; Bucher, S.; Vène, E.; Bodin, A.; Lagente, V.; Fromenty, B.; Bouvet, R.; Morel, I.; et al. Non-Oxidative Ethanol Metabolism in Human Hepatic Cells in Vitro: Involvement of Uridine Diphospho-Glucuronosyltransferase 1A9 in Ethylglucuronide Production. *Toxicol. Vitr.* **2020**, *66*, 104842. [CrossRef]

12. Quesnot, N.; Bucher, S.; Gade, C.; Vlach, M.; Vene, E.; Valença, S.; Gicquel, T.; Holst, H.; Robin, M.-A.; Loyer, P. Production of Chlorzoxazone Glucuronides via Cytochrome P4502E1 Dependent and Independent Pathways in Human Hepatocytes. *Arch. Toxicol.* **2018**, *92*, 3077–3091. [CrossRef]

13. Richeval, C.; Gicquel, T.; Hugbart, C.; Le Dare, B.; Allorge, D.; Morel, I.; Gaulier, J. In Vitro Characterization of NPS Metabolites Produced by Human Liver Microsomes and the HepaRG Cell Line Using Liquid Chromatography high Resolution Mass Spectrometry (LC-HRMS) Analysis: Application to Furanyl Fentanyl. *Curr. Pharm. Biotechnol.* **2018**, *18*, 806–814. [CrossRef]

14. Katoh, M.; Sawada, T.; Soeno, Y.; Nakajima, M.; Tateno, C.; Yoshizato, K.; Yokoi, T. In Vivo Drug Metabolism Model for Human Cytochrome P450 Enzyme Using Chimeric Mice with Humanized Liver. *J. Pharm. Sci.* **2007**, *96*, 428–437. [CrossRef]

15. Park, Y.M.; Meyer, M.R.; Müller, R.; Herrmann, J. Drug Administration Routes Impact the Metabolism of a Synthetic Cannabinoid in the Zebrafish Larvae Model. *Molecules* **2020**, *25*, 4474. [CrossRef]

16. Caspar, A.T.; Gaab, J.B.; Michely, J.A.; Brandt, S.D.; Meyer, M.R.; Maurer, H.H. Metabolism of the Tryptamine-Derived New Psychoactive Substances 5-MeO-2-Me-DALT, 5-MeO-2-Me-ALCHT, and 5-MeO-2-Me-DIPT and Their Detectability in Urine Studied by GC-MS, LC-MS<sup>n</sup>, and LC-HR-MS/MS: Metabolism of 5-MeO-2-Me-DALT, 5-MeO-2-Me-ALCHT, and 5-MeO-2-Me-DIPT. *Drug Test. Anal.* **2018**, *10*, 184–195. [CrossRef]

17. Wink, C.S.D.; Michely, J.A.; Jacobsen-Bauer, A.; Zapp, J.; Maurer, H.H. Diphenidine, a New Psychoactive Substance: Metabolic Fate Elucidated with Rat Urine and Human Liver Preparations and Detectability in Urine Using GC-MS, LC-MS<sup>n</sup>, and LC-HR-MS<sup>n</sup>: Metabolism of Diphenidine. *Drug Test. Anal.* **2016**, *8*, 1005–1014. [CrossRef]
18. Wang, M.; Carver, J.J.; Phelan, V.V.; Sanchez, L.M.; Garg, N.; Peng, Y.; Nguyen, D.D.; Watrous, J.; Kapono, C.A.; Luzzatto-Knaan, T.; et al. Sharing and Community Curation of Mass Spectrometry Data with Global Natural Products Social Molecular Networking. *Nat. Biotechnol.* **2016**, *34*, 828–837. [CrossRef]

19. Allard, S.; Allard, P.-M.; Morel, I.; Gicquel, T. Application of a Molecular Networking Approach for Clinical and Forensic Toxicology Exemplified in Three Cases Involving 3-MeO-PCP, Doxylamine, and Chlormequat. *Drug Test. Anal.* **2018**, *11*, 669–677. [CrossRef]

20. Le Daré, B.; Allard, S.; Bouvet, R.; Baert, A.; Allard, P.-M.; Morel, I.; Gicquel, T. A Case of Fatal Acebutolol Poisoning: An Illustration of the Potential of Molecular Networking. *Int. J. Leg. Med.* **2019**, *13*, 251–256. [CrossRef]

21. Le Daré, B.; Ferron, P.-J.; Couette, A.; Ribault, C.; Morel, I.; Gicquel, T. In Vivo and in Vitro α-Amanitin Metabolism Studies Using Molecular Networking. *Toxicol. Lett.* **2021**, *346*, 1–6. [CrossRef]

22. Ferron, P.-J.; Le Daré, B.; Bronsard, J.; Steichen, C.; Babina, E.; Pelletier, R.; Hauet, T.; Morel, I.; Tarte, K.; Reizine, F.; et al. Molecular Networking for Drug Toxicities Studies: The Case of Hydroxychloroquine in COVID-19 Patients. *IJMS* **2021**, *23*, 82. [CrossRef]

23. Magny, R.; Auzeil, N.; Lefrère, B.; Mégarbane, B.; Houzé, P.; Labat, L. Molecular Network-Based Identification of Tramadol Metabolites in a Fatal Tramadol Poisoning. *Metabolites* **2022**, *12*, 665. [CrossRef]

24. Allard, S.; Le Daré, B.; Allard, P.-M.; Morel, I.; Gicquel, T. Comparative Molecular Networking Analysis of a Rauwolfia Plant Powder and Biological Matrices in a Fatal Ingestion Case. *Forensic Toxicol.* **2020**, *38*, 447–454. [CrossRef]

25. Kirchmair, J.; Göller, A.H.; Lang, D.; Kunze, J.; Testa, B.; Wilson, I.D.; Glen, R.C.; Schneider, G. Predicting Drug Metabolism: Experiment and/or Computation? *Nat. Rev. Drug Discov.* **2015**, *14*, 387–404. [CrossRef]

26. Gertrudes, J.C.; Maltarollo, V.G.; Silva, R.A.; Oliveira, P.R.; Honorio, K.M.; da Silva, A.B.F. Machine Learning Techniques and Drug Design. *CMC* **2012**, *19*, 4289–4297. [CrossRef]

27. Du, Q.-S.; Huang, R.-B.; Chou, K.-C. Recent Advances in QSAR and Their Applications in Predicting the Activities of Chemical Molecules, Peptides and Proteins for Drug Design. *CPPS* **2008**, *9*, 248–259. [CrossRef]

28. Kazmi, S.R.; Jun, R.; Yu, M.-S.; Jung, C.; Na, D. In Silico Approaches and Tools for the Prediction of Drug Metabolism and Fate: A Review. *Comput. Biol. Med.* **2019**, *106*, 54–64. [CrossRef]

29. Tyzack, J.D.; Kirchmair, J. Computational Methods and Tools to Predict Cytochrome P450 Metabolism for Drug Discovery. *Chem. Biol. Drug Des.* **2019**, *93*, 377–386. [CrossRef]

30. Di Trana, A.; Brunetti, P.; Giorgetti, R.; Marinelli, E.; Zaami, S.; Busardò, F.P.; Carlier, J. In Silico Prediction, LC-HRMS/MS Analysis, and Targeted/Untargeted Data-Mining Workflow for the Profiling of Phenylfentanyl in Vitro Metabolites. *Talanta* **2021**, *235*, 122740. [CrossRef]

31. de Bruyn Kops, C.; Šícho, M.; Mazzolari, A.; Kirchmair, J. GLORYx: Prediction of the Metabolites Resulting from Phase 1 and Phase 2 Biotransformations of Xenobiotics. *Chem. Res. Toxicol.* **2021**, *34*, 286–299. [CrossRef]

32. Djoumbou-Feunang, Y.; Fiamoncini, J.; Gil-de-la-Fuente, A.; Greiner, R.; Manach, C.; Wishart, D.S. BioTransformer: A Comprehensive Computational Tool for Small Molecule Metabolism Prediction and Metabolite Identification. *J. Cheminform.* **2019**, *11*, 2. [CrossRef]

33. Gicquel, T.; Pelletier, R.; Richeval, C.; Gish, A.; Hakim, F.; Ferron, P.; Mesli, V.; Allorge, D.; Morel, I.; Gaulier, J. Metabolite Elucidation of 2-fluoro-deschloroketamine (2F-DCK) Using Molecular Networking across Three Complementary in Vitro and in Vivo Models. *Drug Test. Anal.* **2021**, *14*, 144–153. [CrossRef]

34. Kessner, D.; Chambers, M.; Burke, R.; Agus, D.; Mallick, P. ProteoWizard: Open Source Software for Rapid Proteomics Tools Development. *Bioinformatics* **2008**, *24*, 2534–2536. [CrossRef]

35. Pluskal, T.; Castillo, S.; Villar-Briones, A.; Orešič, M. MZmine 2: Modular Framework for Processing, Visualizing, and Analyzing Mass Spectrometry-Based Molecular Profile Data. *BMC Bioinform.* **2010**, *11*, 395. [CrossRef]

36. Shannon, P. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. *Genome Res.* **2003**, *13*, 2498–2504. [CrossRef]

37. HighChem LLC MzCloud-Advanced Mass Spectral Database. Available online: https://www.mzcloud.org/ (accessed on 12 November 2018).

38. Bogaards, J.J.P.; Bertrand, M.; Jackson, P.; Oudshoorn, M.J.; Weaver, R.J.; Van Bladeren, P.J.; Walther, B. Determining the Best Animal Model for Human Cytochrome P450 Activities: A Comparison of Mouse, Rat, Rabbit, Dog, Micropig, Monkey and Man. *Xenobiotica* **2000**, *30*, 1131–1152. [CrossRef]

39. Guengerich, F.P. Comparisons of Catalytic Selectivity of Cytochrome P450 Subfamily Enzymes from Different Species. *Chem. Biol. Interact.* **1997**, *106*, 161–182. [CrossRef]

40. Lin, J.H. Applications and Limitations of Interspecies Scaling and in Vitro Extrapolation in Pharmacokinetics. *Drug Metab. Dispos.* **1998**, *26*, 1202–1212.

41. Lin, J.H. Species Similarities and Differences in Pharmacokinetics. *Drug Metab. Dispos.* **1995**, *23*, 1008–1021.

42. Jeon, J.-H.; Lee, J.; Choi, M.-K.; Song, I.-S. Pharmacokinetics of Ginsenosides Following Repeated Oral Administration of Red Ginseng Extract Significantly Differ between Species of Experimental Animals. *Arch. Pharm. Res.* **2020**, *43*, 1335–1346. [CrossRef]

43. Natsch, A.; Nordone, A.; Adamson, G.M.; Laue, H. A Species Specific Metabolism Leading to Male Rat Reprotoxicity of Cyclamen Aldehyde: In Vivo and in Vitro Evaluation. *Food Chem. Toxicol.* **2021**, *153*, 112243. [CrossRef] [PubMed]
44. Arensdorff, L.; Boillat-Blanco, N.; Decosterd, L.; Buclin, T.; de Vallière, S. Adequate Plasma Drug Concentrations Suggest That Amoxicillin Can Be Administered by Continuous Infusion Using Elastomeric Pumps. *J. Antimicrob. Chemother.* **2017**, *72*, 2613–2615.
[CrossRef] [PubMed]

45. Powell, J.S. Liposomal Approach towards the Development of a Longer-Acting Factor VIII. *Haemophilia* **2007**, *13*, 23–28. [CrossRef]

46. Pacchiarotti, I.; Tiihonen, J.; Kotzalidis, G.D.; Verdolini, N.; Murru, A.; Goikolea, J.M.; Valentí, M.; Aedo, A.; Vieta, E. Long-Acting Injectable Antipsychotics (LAIs) for Maintenance Treatment of Bipolar and Schizoaffective Disorders: A Systematic Review. *Eur. Neuropsychopharmacol.* **2019**, *29*, 457–470. [CrossRef] [PubMed]

47. Hall, S.; Isaacs, D.; Clements, J.N. Pharmacokinetics and Clinical Implications of Semaglutide: A New Glucagon-Like Peptide (GLP)-1 Receptor Agonist. *Clin. Pharmacokinet.* **2018**, *57*, 1529–1538. [CrossRef] [PubMed]

48. Bialer, M.; Friedman, M.; Dubrovsky, J. Relation between Absorption Half-Life Values of Four Novel Sustained-Release Dosage Forms of Valproic Acid in Dogs and Human. *Biopharm. Drug Dispos.* **1986**, *7*, 495–500. [CrossRef]

49. Gunaydin, H.; Altman, M.D.; Ellis, J.M.; Fuller, P.; Johnson, S.A.; Lahue, B.; Lapointe, B. Strategy for Extending Half-Life in Drug Design and Its Significance. *ACS Med. Chem. Lett.* **2018**, *9*, 528–533. [CrossRef]