LipidMS 3.0: an R-package and a web-based tool for LC-MS/MS data processing and lipid annotation.

MOTIVATION: LipidMS was initially envisioned to use fragmentation rules and data-independent acquisition (DIA) for lipid annotation. However, data-dependent acquisition (DDA) remains the most widespread acquisition mode for untargeted LC-MS/MS-based lipidomics. Here, we present LipidMS 3.0, an R package that not only adds DDA and new lipid classes to its pipeline but also the required functionalities to cover the whole data analysis workflow from pre-processing (i.e. peak-peaking, alignment and grouping) to lipid annotation. RESULTS: We applied the new workflow in the data analysis of a commercial human serum pool spiked with 68 representative lipid standards acquired in full scan, DDA and DIA modes. When focusing on the detected lipid standard features and total identified lipids, LipidMS 3.0 data pre-processing performance is similar to XCMS, whereas it complements the annotations returned by MS-DIAL, providing a higher level of structural information and a lower number of incorrect annotations. To extend and facilitate LipidMS 3.0 usage among less experienced R-programming users, the workflow is also implemented as a web-based application. AVAILABILITY AND IMPLEMENTATION: The LipidMS R-package is freely available at https://CRAN.R-project.org/package=LipidMS and as a website at http://www.lipidms.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

LipidMS: An R Package for Lipid Annotation in Untargeted Liquid Chromatography-Data Independent Acquisition-Mass Spectrometry Lipidomics.

High resolution LC-MS untargeted lipidomics using data independent acquisition (DIA) has the potential to increase lipidome coverage, as it enables the continuous and unbiased acquisition of all eluting ions. However, the loss of the link between the precursor and the product ions combined with the high dimensionality of DIA data sets hinder accurate feature annotation. Here, we present LipidMS, an R package aimed to confidently identify lipid species in untargeted LC-DIA-MS. To this end, LipidMS combines a coelution score, which links precursor and fragment ions with fragmentation and intensity rules. Depending on the MS evidence reached by the identification function survey, LipidMS provides three levels of structural annotations: (i) "subclass level", e.g., PG(34:1); (ii) "fatty acyl level", e.g., PG(16:0_18:1); and (iii) "fatty acyl position level", e.g., PG(16:0/18:1). The comparison of LipidMS with freely available data dependent acquisition (DDA) and DIA identification tools showed that LipidMS provides significantly more accurate and structural informative lipid identifications. Finally, to exemplify the utility of LipidMS, we investigated the lipidomic serum profile of patients diagnosed with nonalcoholic steatohepatitis (NASH), which is the progressive form of nonalcoholic fatty liver disease, a disorder underlying a strong lipid dysregulation. As previously published, a significant decrease in lysophosphatidylcholines, phosphatidylcholines and cholesterol esters and an increase in phosphatidylethanolamines were observed in NASH patients. Remarkably, LipidMS allowed the identification of a new set of lipids that may be used for NASH diagnosis. Altogether, LipidMS has been validated as a tool to assist lipid identification in the LC-DIA-MS untargeted analysis of complex biological samples.

RpeakChrom: Novel R package for the automated characterization and optimization of column efficiency in high-performance liquid chromatography analysis.

Characterization of chromatographic columns using the traditional van Deemter method is limited by the necessity of calculating extra-column variance, issue particularly relevant when modeling asymmetrical peaks eluted from monolithic columns. A novel R package that implements Parabolic Variance Modified Gaussian approach for accurate peak modeling, van Deemter equation and two alternatives approaches, based on van Deemter, has been developed to calculate the height equivalent to a theoretical plate (HETP). To assess package capabilities conventional packed reverse-phase and monolithic HPLC columns were characterized. Peaks eluted from the monolithic column showed a high value of factor asymmetry due, in part, to the contribution of extra-column factors. Such deviation can be circumvented by the two alternatives approaches implemented in the R-package. Furthermore, increased values of eddy diffusion and mass transfer kinetics terms in HETP were observed for the packed column, while accuracy was below 9% in all cases. These results showed the usefulness of the R-package for both modeling chromatographic peaks and assessing column efficiency. The RpeakChrom package could become a helpful tool for testing new stationary phases during column development and to evaluate column during its lifetime. This R tool is freely available from CRAN (https://CRAN.R-project.org/package=RpeakChrom).

A lipidomic cell-based assay for studying drug-induced phospholipidosis and steatosis.

Phospholipidosis and steatosis are two toxic effects, which course with overaccumulation of different classes of lipids in the liver. MS-based lipidomics has become a powerful tool for the comprehensive determination of lipids. LC-MS lipid profiling of HepG2 cells is proposed as an in vitro assay to study and anticipate phospholipidosis and steatosis. Cells with and without preincubation with a mixture of free fatty acids (FFA; i.e. oleic and palmitic) were exposed to a set of well-known steatogenic and phospholipidogenic compounds. The use of FFA preloading accelerated the accumulation of phospholipids, thus leading to a better discrimination of phospholipidosis, and magnified the lipidomic alterations induced by steatogenic drugs. Phospholipidosis was characterized by increased levels of phosphatidylcholines, phosphatidylethanolamines, phosphatidylserines, and phosphatidylinositols, while steatosis induced alterations in FA oxidation and triacylglyceride (TG) synthesis pathways (with changes in the levels of FFA, acylcarnitines, monoacylglycerides, diacylglycerides, and TG). Interestingly, palmitic and oleic acids incorporation into lipids differed. A characteristic pattern was observed in the fold of change of particular TG species in the case of steatosis (TG(54:3) > TG(52:2) > TG(50:1) > TG(48:0)). Based on the levels of those lipids containing only palmitic and/or oleic acid moieties a partial least squares-discriminant analysis model was built, which showed good discrimination among nontoxic, phospholipidogenic and steatogenic compounds. In conclusion, it has been shown that the use of FFA preincubation together with intracellular LC-MS based lipid profiling could be a useful approach to identify the potential of drug candidates to induce phospholipidosis and/or steatosis.

Extending metabolome coverage for untargeted metabolite profiling of adherent cultured hepatic cells.

MS-based metabolite profiling of adherent mammalian cells comprises several challenging steps such as metabolism quenching, cell detachment, cell disruption, metabolome extraction, and metabolite measurement. In LC-MS, the final metabolome coverage is strongly determined by the separation technique and the MS conditions used. Human liver-derived cell line HepG2 was chosen as adherent mammalian cell model to evaluate the performance of several commonly used procedures in both sample processing and LC-MS analysis. In a first phase, metabolite extraction and sample analysis were optimized in a combined manner. To this end, the extraction abilities of five different solvents (or combinations) were assessed by comparing the number and the levels of the metabolites comprised in each extract. Three different chromatographic methods were selected for metabolites separation. A HILIC-based method which was set to specifically separate polar metabolites and two RP-based methods focused on lipidome and wide-ranging metabolite detection, respectively. With regard to metabolite measurement, a Q-ToF instrument operating in both ESI (+) and ESI (-) was used for unbiased extract analysis. Once metabolite extraction and analysis conditions were set up, the influence of cell harvesting on metabolome coverage was also evaluated. Therefore, different protocols for cell detachment (trypsinization or scraping) and metabolism quenching were compared. This study confirmed the inconvenience of trypsinization as a harvesting technique, and the importance of using complementary extraction solvents to extend metabolome coverage, minimizing interferences and maximizing detection, thanks to the use of dedicated analytical conditions through the combination of HILIC and RP separations. The proposed workflow allowed the detection of over 300 identified metabolites from highly polar compounds to a wide range of lipids.

LC-MS untargeted metabolomic analysis of drug-induced hepatotoxicity in HepG2 cells.

Hepatotoxicity is the number one cause for agencies not approving and withdrawing drugs for the market. Drug-induced human hepatotoxicity frequently goes undetected in preclinical safety evaluations using animal models. Human-derived in vitro models represent a common alternative to in vivo tests to detect toxic effects during preclinical testing. Most current in vitro toxicity assays rely on the measurement of nonspecific or low sensitive endpoints, which result in poor concordance with human liver toxicity. Therefore, making more accurate predictions of the potential hepatotoxicity of new drugs remains a challenge. Metabolomics, whose aim is to globally assess all the metabolites present in a biological sample, may represent an alternative in the search for sensitive sublethal markers of drug-induced hepatotoxicity. To this end, a comprehensive LC-MS-based untargeted metabolite profiling analysis of HepG2 cells, exposed to a set of well-described model hepatotoxins and innocuous compounds, was performed. It allowed to determine meaningful metabolic changes triggered by a toxic insult and gave a first estimation of the main toxicity-related pathways. Based on these metabolic patterns, a partial least squares-discriminant analysis model, able to discriminate between nontoxic and hepatotoxic compounds, was constructed. The approach described herein may provide an alternative for animal testing in preclinical stages of drug development and a controlled experimental approach to gain a better understanding of the underlying causes of hepatotoxicity.

In vitro/in vivo screening of oxidative homeostasis and damage to DNA, protein, and lipids using UPLC/MS-MS.

Multiple analytical methods are required to comprehensively assess oxidative homeostasis and specific damage to macromolecules. Our aim was to develop a straightforward strategy for the fast assessment of global oxidative status and specific damage to DNA, proteins, and lipids. To this end, an analytical method, based on ultra-performance liquid chromatography coupled to mass spectrometry (UPLC-MS/MS), was developed and validated for the quantification of 16 oxidative stress (OS) biomarkers. Some of these markers were unstable; thus, an easy sample treatment procedure, including fractionation and derivatization, was set up. The method was validated according to Food and Drug Administration (FDA) guidelines, and it provided good results in terms of intra- and inter-day precision (≤17.2 and 16 %, respectively), accuracy (relative error measurement between -16.6 and 19.8 %), and linearity (R (2) > 0.994). The approach was applied to determine the oxidative insult provoked to cultured rat hepatocytes by cumene hydroperoxide and to analyze the liver and serum samples from patients diagnosed with nonalcoholic steatohepatitis. In both studies, significant differences were found if compared to the corresponding control groups; interestingly, ophthalmic acid was shown as an OS biomarker in both models for the first time. A key advantage of the novel approach in comparison with former multi-method approaches is that now a single method is applied to assess the 16 OS biomarkers. Its comprehensive capacity to profile oxidative homeostasis and damage in both in vitro and clinical samples has been illustrated, which indicates that the proposed approach is a good choice to evaluate whether OS is involved in physiological signals, diseases, or toxic events and to what extent.

Immune Checkpoint Inhibitors as First-line Treatment for Brain Metastases in Stage IV NSCLC Patients without Driver Mutations.

Immune checkpoint inhibitors (ICI) therapy with or without chemotherapy has been established as the first-line treatment for patients with non-oncogene addicted advanced Non-Small Cell Lung Cancer (NSCLC). Yet some clinical settings, such as the treatment sequence in patients with brain metastases, have barely been evidenced. Although ICIs cannot directly cross the blood-brain barrier (BBB), evidence suggests that BBB damage could allow ICIs into the central nervous system, or that they can have an indirect effect on the tumor immune microenvironment (TIME) and cause an anti-tumor response. Pivotal phase III trials have included a highly selected population but offer few data on these patients. Here we first review how ICIs can indirectly shape the brain metastases microenvironment through different mechanisms, and some possible causes of ICIs resistance. We also analyze the evidence reported in pivotal phase III trials and phase II trials focused on NSCLC brain metastases for first-line treatment, and the evidence for upfront or delayed local brain therapy. Finally, we discuss the best evidence-based approach to treat NSCLC patients with brain metastases and propose future research.

GLUD1 determines murine muscle stem cell fate by controlling mitochondrial glutamate levels.

Muscle stem cells (MuSCs) enable muscle growth and regeneration after exercise or injury, but how metabolism controls their regenerative potential is poorly understood. We describe that primary metabolic changes can determine murine MuSC fate decisions. We found that glutamine anaplerosis into the tricarboxylic acid (TCA) cycle decreases during MuSC differentiation and coincides with decreased expression of the mitochondrial glutamate deaminase GLUD1. Deletion of Glud1 in proliferating MuSCs resulted in precocious differentiation and fusion, combined with loss of self-renewal in vitro and in vivo. Mechanistically, deleting Glud1 caused mitochondrial glutamate accumulation and inhibited the malate-aspartate shuttle (MAS). The resulting defect in transporting NADH-reducing equivalents into the mitochondria induced compartment-specific NAD+/NADH ratio shifts. MAS activity restoration or directly altering NAD+/NADH ratios normalized myogenesis. In conclusion, GLUD1 prevents deleterious mitochondrial glutamate accumulation and inactivation of the MAS in proliferating MuSCs. It thereby acts as a compartment-specific metabolic brake on MuSC differentiation.

Sex Differences in the Efficacy of Immune Checkpoint Inhibitors in Neoadjuvant Therapy of Non-Small Cell Lung Cancer: A Meta-Analysis.

Non-small cell lung cancer (NSCLC) is one of the world's leading causes of morbidity and mortality. ICIs alone or combined with chemotherapy have become the standard first-line treatment of metastatic NSCLC. The impressive results obtained have stimulated our interest in applying these therapies in early disease stage treatments, as neoadjuvant immunotherapy has shown promising results. Among many of the factors that may influence responses, the role played by sex is attracting increased interest and needs to be addressed. Here, we aim to first review the state of the art regarding neoadjuvant ICIs, whether they are administered in monotherapy or in combination with chemotherapy at stages IB-IIIA, particularly at stage IIIA, before analyzing whether sex may influence responses. To this end, a meta-analysis of publicly available data comparing male and female major pathological responses (MPR) and pathological complete responses (pCR) was performed. In our meta-analysis, MPR was found to be significantly higher in females than in males, with an odds ratio (OR) of 1.82 (95% CI 1.13-2.93; p = 0.01), while pCR showed a trend to be more favorable in females than in males, but the OR of 1.62 was not statistically significant (95% CI 0.97-2.75; p = 0.08). Overall, our results showed that sex should be systematically considered in future clinical trials settings in order to establish the optimal treatment sequence.

Reviewing the metabolome coverage provided by LC-MS: Focus on sample preparation and chromatography-A tutorial.

Metabolomics has become an invaluable tool for both studying metabolism and biomarker discovery. The great technical advances in analytical chemistry and bioinformatics have considerably increased the number of measurable metabolites, yet an important part of the human metabolome remains uncovered. Among the various MS hyphenated techniques available, LC-MS stands out as the most used. Here, we aimed to show the capabilities of LC-MS to uncover part of the metabolome and how to best proceed with sample preparation and LC to maximise metabolite detection. The analyses of various open metabolite databases served us to estimate the size of the already detected human metabolome, the expected metabolite composition of most used human biospecimens and which part of the metabolome can be detected when LC-MS is used. Based on an extensive review and on our experience, we have outlined standard procedures for LC-MS analysis of urine, cells, serum/plasma, tissues and faeces, to guide in the selection of the sample preparation method that best matches with one or more LC techniques in order to get the widest metabolome coverage. These standard procedures may be a useful tool to explore, at a glance, the wide spectrum of possibilities available, which can be a good starting point for most of the LC-MS metabolomic studies.

c-MYC Triggers Lipid Remodelling During Early Somatic Cell Reprogramming to Pluripotency.

Metabolic rewiring and mitochondrial dynamics remodelling are hallmarks of cell reprogramming, but the roles of the reprogramming factors in these changes are not fully understood. Here we show that c-MYC induces biosynthesis of fatty acids and increases the rate of pentose phosphate pathway. Time-course profiling of fatty acids and complex lipids during cell reprogramming using lipidomics revealed a profound remodelling of the lipid content, as well as the saturation and length of their acyl chains, in a c-MYC-dependent manner. Pluripotent cells displayed abundant cardiolipins and scarce phosphatidylcholines, with a prevalence of monounsaturated acyl chains. Cells undergoing cell reprogramming showed an increase in mitochondrial membrane potential that paralleled that of mitochondrial-specific cardiolipins. We conclude that c-MYC controls the rewiring of somatic cell metabolism early in cell reprogramming by orchestrating cell proliferation, synthesis of macromolecular components and lipid remodelling, all necessary processes for a successful phenotypic transition to pluripotency. c-MYC promotes anabolic metabolism, mitochondrial fitness and lipid remodelling early in cell reprogramming. A high rate of aerobic glycolysis is crucial to provide intermediaries for biosynthetic pathways. To ensure the availability of nucleotides, amino acids and lipids for cell proliferation, cells must provide with a constant flux of the elemental building blocks for macromolecule assembly and fulfil the anabolic demands to reach the critical cellular mass levels to satisfactorily undergo cell division. A high rate of aerobic glycolysis is induced by c-MYC, increasing the amounts of intracellular Glucose-6-phosphate (G6P), fructose-6-phosphate (F6P), and glyceraldehyde-3-phosphate (GA3P), which can all enter pentose phosphate pathway (PPP) to produce Ribose-5-Phosphate (R5P) and NADPH, which are necessary for the biosynthesis of biomolecules such as proteins, nucleic acids, or lipids. C-MYC-dependent activation of glucose-6-phosphate dehydrogenase (G6PD) may play a critical role in the shunting of G6P to PPP and generation of NADPH. High glycolytic flux increases the amounts of dihydroxyacetone phosphate (DHAP), which is crucial for biosynthesis of phospholipids and triacylglycerols, and pyruvate (Pyr), which can be converted to citrate (Cit) in the mitochondria and enter the biosynthesis of fatty acids (FA). During cell reprogramming, c-MYC-dependent lipid remodelling leads to Polyunsaturated Fatty Acid (PUFA) downregulation and Monounsaturated Fatty Acid (MUFA) upregulation, which may play critical roles in cytoarchitectural remodelling of cell membrane or non-canonical autophagy, respectively. Cardiolipin (pink dots) rise early in cell reprogramming correlates with an increase in mitochondrial fitness, suggesting that c-MYC may restore proper levels of cardiolipins and antioxidant proteins, such as UCP2, to guarantee an optimal mitochondrial function while upholding ROS levels, reinforcing the idea of cell rejuvenation early in cell reprogramming.

Enhancing Chimeric Antigen Receptor T Cell Anti-tumor Function through Advanced Media Design.

Effective chimeric antigen receptor (CAR)-T cell therapy is dependent on optimal cell culture methods conducive to the activation and expansion of T cells ex vivo, as well as infection with CAR. Media formulations used in CAR-T cell manufacturing have not been optimized for gene delivery, cell expansion, and overall potency. Bioactive components and derivatives that support the generation of functionally-competent T cell progeny with long-lasting persistence are largely undefined. Current media formulations rely on fetal bovine serum (FBS) or human serum (HS), which suffer from a lack of consistency or supply issues. We recognize that components of blood cellular fractions that are absent in serum may have therapeutic value. Here we investigate whether a concentrated growth factor extract, purified from human transfusion grade whole blood fractions, and marketed as PhysiologixTM xeno-free (XF) hGFC (Phx), supports CAR-T cell expansion and function. We show that Phx supports T cell proliferation in clinical and research-grade media. We also show that Phx treatment enhances lentiviral-mediated gene expression across a wide range of multiplicity of infections (MOIs). We compared the ability of anti-GD-2 CAR-T cells expanded ex vivo in medium conditioned with either Phx or HS to clear tumor burden in a human xenograft model of neuroblastoma. We show that T cells expanded in Phx have superior engraftment and potency in vivo, as well as CAR-induced cytolytic activity in vitro. Metabolomic profiling revealed several factors unique to Phx that may have relevance for CAR-T cell preclinical discovery, process development, and manufacturing. In particular, we show that carnosine, a biogenic amine modestly enriched in Phx relative to HS, enhances lentiviral gene delivery in activated T cells. By limiting extracellular acidification, carnosine enhances the metabolic fitness of T cells, shifting their metabolic profile from an acidic, stressed state toward an oxidative, energetic state. These findings are very informative regarding potential derivatives to include in medium customized for gene delivery and overall potency for T cell adoptive immunotherapies.

Endothelin-1-Mediated Drug Resistance in EGFR-Mutant Non-Small Cell Lung Carcinoma.

Progression on therapy in non-small cell lung carcinoma (NSCLC) is often evaluated radiographically, however, image-based evaluation of said therapies may not distinguish disease progression due to intrinsic tumor drug resistance or inefficient tumor penetration of the drugs. Here we report that the inhibition of mutated EGFR promotes the secretion of a potent vasoconstrictor, endothelin-1 (EDN1), which continues to increase as the cells become resistant with a mesenchymal phenotype. As EDN1 and its receptor (EDNR) is linked to cancer progression, EDNR-antagonists have been evaluated in several clinical trials with disappointing results. These trials were based on a hypothesis that the EDN1-EDNR axis activates the MAPK-ERK signaling pathway that is vital to the cancer cell survival; the trials were not designed to evaluate the impact of tumor-derived EDN1 in modifying tumor microenvironment or contributing to drug resistance. Ectopic overexpression of EDN1 in cells with mutated EGFR resulted in poor drug delivery and retarded growth in vivo but not in vitro. Intratumoral injection of recombinant EDN significantly reduced blood flow and subsequent gefitinib accumulation in xenografted EGFR-mutant tumors. Furthermore, depletion of EDN1 or the use of endothelin receptor inhibitors bosentan and ambrisentan improved drug penetration into tumors and restored blood flow in tumor-associated vasculature. Correlatively, these results describe a simplistic endogenous yet previously unrealized resistance mechanism inherent to a subset of EGFR-mutant NSCLC to attenuate tyrosine kinase inhibitor delivery to the tumors by limiting drug-carrying blood flow and the drug concentration in tumors. SIGNIFICANCE: EDNR antagonists can be repurposed to improve drug delivery in VEGFA-secreting tumors, which normally respond to TKI treatment by secreting EDN1, promoting vasoconstriction, and limiting blood and drug delivery.

Serine Catabolism Feeds NADH when Respiration Is Impaired.

NADH provides electrons for aerobic ATP production. In cells deprived of oxygen or with impaired electron transport chain activity, NADH accumulation can be toxic. To minimize such toxicity, elevated NADH inhibits the classical NADH-producing pathways: glucose, glutamine, and fat oxidation. Here, through deuterium-tracing studies in cultured cells and mice, we show that folate-dependent serine catabolism also produces substantial NADH. Strikingly, when respiration is impaired, serine catabolism through methylene tetrahydrofolate dehydrogenase (MTHFD2) becomes a major NADH source. In cells whose respiration is slowed by hypoxia, metformin, or genetic lesions, mitochondrial serine catabolism inhibition partially normalizes NADH levels and facilitates cell growth. In mice with engineered mitochondrial complex I deficiency (NDUSF4-/-), serine's contribution to NADH is elevated, and progression of spasticity is modestly slowed by pharmacological blockade of serine degradation. Thus, when respiration is impaired, serine catabolism contributes to toxic NADH accumulation.

A metabolomics cell-based approach for anticipating and investigating drug-induced liver injury.

In preclinical stages of drug development, anticipating potential adverse drug effects such as toxicity is an important issue for both saving resources and preventing public health risks. Current in vitro cytotoxicity tests are restricted by their predictive potential and their ability to provide mechanistic information. This study aimed to develop a metabolomic mass spectrometry-based approach for the detection and classification of drug-induced hepatotoxicity. To this end, the metabolite profiles of human derived hepatic cells (i.e., HepG2) exposed to different well-known hepatotoxic compounds acting through different mechanisms (i.e., oxidative stress, steatosis, phospholipidosis, and controls) were compared by multivariate data analysis, thus allowing us to decipher both common and mechanism-specific altered biochemical pathways. Briefly, oxidative stress damage markers were found in the three mechanisms, mainly showing altered levels of metabolites associated with glutathione and γ-glutamyl cycle. Phospholipidosis was characterized by a decreased lysophospholipids to phospholipids ratio, suggestive of phospholipid degradation inhibition. Whereas, steatosis led to impaired fatty acids ß-oxidation and a subsequent increase in triacylglycerides synthesis. The characteristic metabolomic profiles were used to develop a predictive model aimed not only to discriminate between non-toxic and hepatotoxic drugs, but also to propose potential drug toxicity mechanism(s).

Ultra-performance liquid chromatography-mass spectrometry targeted profiling of bile acids: application to serum, liver tissue, and cultured cells of different species.

Currently, there is increasing interest in developing accurate methods for the quantitative analysis of bile acids (BAs) in biological samples. We have developed a sensitive, fast, and reproducible UPLC-MRM-MS method for BA profiling in serum, liver tissue, or cultured cells of different species (human, rat, and mouse). This method, validated according to FDA guidelines, allows the quantification of 12 non-conjugated, 8 glycine-conjugated, and 11 taurine-conjugated BAs, using 5 additional deuterated BAs as internal standards in a single analytical run. The main features of this analytical approach are its high sensitivity, low sample requirements, versatility, and comprehensive capacity to profile a considerable number of BAs in samples of different species, which make it a valuable tool with potential applications in many research areas focusing on BAs, particularly in toxicological studies.