Design of greener drugs: aligning parameters in pharmaceutical R&D and drivers for environmental impact

Active pharmaceutical ingredients (APIs) in the environment, primarily resulting from patient excretion, are of concern because of potential risks to wildlife. This has led to more restrictive regulatory policies. Here, we discuss the 'benign-by-design' approach, which encourages the development of environmentally friendly APIs that are also safe and efficacious for patients. We explore the challenges and opportunities associated with identifying chemical properties that influence the environmental impact of APIs. Although a straightforward application of greener properties could hinder the development of new drugs, more nuanced approaches could lead to drugs that benefit both patients and the environment. We advocate for an enhanced dialogue between research and development (R&D) and environmental scientists and development of a toolbox to incorporate environmental sustainability in drug development.

Cross-Species Extrapolation of Biological Data to Guide the Environmental Safety Assessment of Pharmaceuticals-The State of the Art and Future Priorities

The extrapolation of biological data across species is a key aspect of biomedical research and drug development. In this context, comparative biology considerations are applied with the goal of understanding human disease and guiding the development of effective and safe medicines. However, the widespread occurrence of pharmaceuticals in the environment and the need to assess the risk posed to wildlife have prompted a renewed interest in the extrapolation of pharmacological and toxicological data across the entire tree of life. To address this challenge, a biological "read-across" approach, based on the use of mammalian data to inform toxicity predictions in wildlife species, has been proposed as an effective way to streamline the environmental safety assessment of pharmaceuticals. Yet, how effective has this approach been, and are we any closer to being able to accurately predict environmental risk based on known human risk? We discuss the main theoretical and experimental advancements achieved in the last 10 years of research in this field. We propose that a better understanding of the functional conservation of drug targets across species and of the quantitative relationship between target modulation and adverse effects should be considered as future research priorities. This pharmacodynamic focus should be complemented with the application of higher-throughput experimental and computational approaches to accelerate the prediction of internal exposure dynamics. The translation of comparative (eco)toxicology research into real-world applications, however, relies on the (limited) availability of experts with the skill set needed to navigate the complexity of the problem; hence, we also call for synergistic multistakeholder efforts to support and strengthen comparative toxicology research and education at a global level. Environ Toxicol Chem 2024;43:513-525. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Big Question to Developing Solutions: A Decade of Progress in the Development of Aquatic New Approach Methodologies from 2012 to 2022

In 2012, 20 key questions related to hazard and exposure assessment and environmental and health risks of pharmaceuticals and personal care products in the natural environment were identified. A decade later, this article examines the current level of knowledge around one of the lowest-ranking questions at that time, number 19: "Can nonanimal testing methods be developed that will provide equivalent or better hazard data compared with current in vivo methods?" The inclusion of alternative methods that replace, reduce, or refine animal testing within the regulatory context of risk and hazard assessment of chemicals generally faces many hurdles, although this varies both by organism (human-centric vs. other), sector, and geographical region or country. Focusing on the past 10 years, only works that might reasonably be considered to contribute to advancements in the field of aquatic environmental risk assessment are highlighted. Particular attention is paid to methods of contemporary interest and importance, representing progress in (1) the development of methods which provide equivalent or better data compared with current in vivo methods such as bioaccumulation, (2) weight of evidence, or (3) -omic-based applications. Evolution and convergence of these risk assessment areas offer the basis for fundamental frameshifts in how data are collated and used for the protection of taxa across the breadth of the aquatic environment. Looking to the future, we are at a tipping point, with a need for a global and inclusive approach to establish consensus. Bringing together these methods (both new and old) for regulatory assessment and decision-making will require a concerted effort and orchestration. Environ Toxicol Chem 2024;43:559-574. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Towards Precision Ecotoxicology: Leveraging Evolutionary Conservation of Pharmaceutical and Personal Care Product Targets to Understand Adverse Outcomes Across Species and Life Stages

Translation of environmental science to the practice aims to protect biodiversity and ecosystem services, and our future ability to do so relies on the development of a precision ecotoxicology approach wherein we leverage the genetics and informatics of species to better understand and manage the risks of global pollution. A little over a decade ago, a workshop focusing on the risks of pharmaceuticals and personal care products (PPCPs) in the environment identified a priority research question, "What can be learned about the evolutionary conservation of PPCP targets across species and life stages in the context of potential adverse outcomes and effects?" We review the activities in this area over the past decade, consider prospects of more recent developments, and identify future research needs to develop next-generation approaches for PPCPs and other global chemicals and waste challenges. Environ Toxicol Chem 2024;43:526-536. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

The influence of size and surface chemistry on the bioavailability, tissue distribution and toxicity of gold nanoparticles in zebrafish (Danio rerio)

Gold nanoparticles (AuNPs) are widely used in biomedicine and their specific properties including, size, geometrics, and surface coating, will affect their fate and behaviour in biological systems. These properties are well studied for their intended biological targets, but there is a lack of understanding on the mechanisms by which AuNPs interact in non-target organisms when they enter the environment. We investigated the effects of size and surface chemistry of AuNPs on their bioavailability, tissue distribution and potential toxicity using zebrafish (Danio rerio) as an experimental model. Larval zebrafish were exposed to fluorescently tagged AuNPs of different sizes (10-100 nm) and surface modifications (TNFα, NHS/PAMAM and PEG), and uptake, tissue distribution and depuration rates were measured using selective-plane illumination microscopy (SPIM). The gut and pronephric tubules were found to contain detectable levels of AuNPs, and the concentration-dependent accumulation was related to the particle size. Surface addition of PEG and TNFα appeared to enhance particle accumulation in the pronephric tubules compared to uncoated particles. Depuration studies showed a gradual removal of particles from the gut and pronephric tubules, although fluorescence indicating the presence of the AuNPs remained in the pronephros 96 h after exposure. Toxicity assessment using two transgenic zebrafish reporter lines, however, revealed no AuNP-related renal injury or cellular oxidative stress. Collectively, our data show that AuNPs used in medical applications across the size range 40-80 nm, are bioavailable to larval zebrafish and some may persist in renal tissue, although their presence did not result in measurable toxicity with respect to pronephric organ function or cellular oxidative stress for short term exposures.

Active Pharmaceutical Ingredient Uptake by Zebrafish (Danio rerio) Oct2 (slc22a2) Transporter Expressed in Xenopus laevis Oocytes

Uptake of active pharmaceutical ingredients (APIs) across the gill epithelium of fish is via either a passive or facilitated transport process, with the latter being more important at the lower concentrations more readily observed in the environment. The solute carrier (SLC) 22A family, which includes the organic cation transporter OCT2 (SLC22A2), has been shown in mammals to transport several endogenous chemicals and APIs. Zebrafish oct2 was expressed in Xenopus oocytes and the uptake of ranitidine, propranolol, and tetraethylammonium characterized. Uptake of ranitidine and propranolol was time- and concentration-dependent with a k_m and V_max for ranitidine of 246 µM and 45 pmol/(oocyte × min) and for propranolol of 409 µM and 190 pmol/(oocyte × min), respectively. Uptake of tetraethylammonium (TEA) was inhibited by propranolol, amantadine, and cimetidine, known to be human OCT2 substrates, but not quinidine or ranitidine. At external media pH 7 and 8 propranolol uptake was 100-fold greater than at pH 6; pH did not affect ranitidine or TEA uptake. It is likely that cation uptake is driven by the electrochemical gradient across the oocyte. Uptake kinetics parameters, such as those derived in the present study, coupled with knowledge of transporter localization and abundance and API metabolism, can help derive pharmacokinetic models. Environ Toxicol Chem 2022;41:2993-2998. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

GREENER Pharmaceuticals for More Sustainable Healthcare

Medicines are essential to human health but can also impact the aquatic and terrestrial environment after use by patients and release via excreta into wastewater. We highlight the need for a GREENER approach to identify and meet important environmental criteria, which will help reduce the impact of medicinal residues on the environment. These criteria include effect reduction by avoiding nontarget effects or undesirable moieties, exposure reduction via lower emissions or environmental (bio)degradability, no PBT (persistent, bioaccumulative, and toxic) substances, and risk mitigation. With all of these criteria, however, patient health is of primary importance as medicines are required to be safe and efficacious for treating diseases. We discuss the feasibility of including these criteria for green by design active pharmaceutical ingredients in the process of drug discovery and development and which tools or assays are needed to accomplish this. The integrated GREENER approach can be used to accelerate discussions about future innovations in drug discovery and development.

A Generalized Physiologically Based Kinetic Model for Fish for Environmental Risk Assessment of Pharmaceuticals

An increasing number of pharmaceuticals found in the environment potentially impose adverse effects on organisms such as fish. Physiologically based kinetic (PBK) models are essential risk assessment tools, allowing a mechanistic approach to understanding chemical effects within organisms. However, fish PBK models have been restricted to a few species, limiting the overall applicability given the countless species. Moreover, many pharmaceuticals are ionizable, and fish PBK models accounting for ionization are rare. Here, we developed a generalized PBK model, estimating required parameters as functions of fish and chemical properties. We assessed the model performance for five pharmaceuticals (covering neutral and ionic structures). With biotransformation half-lives (HLs) from EPI Suite, 73 and 41% of the time-course estimations were within a 10-fold and a 3-fold difference from measurements, respectively. The performance improved using experimental biotransformation HLs (87 and 59%, respectively). Estimations for ionizable substances were more accurate than any of the existing species-specific PBK models. The present study is the first to develop a generalized fish PBK model focusing on mechanism-based parameterization and explicitly accounting for ionization. Our generalized model facilitates its application across chemicals and species, improving efficiency for environmental risk assessment and supporting an animal-free toxicity testing paradigm.

Are synthetic glucocorticoids in the aquatic environment a risk to fish?

The glucocorticosteroid, or glucocorticoid (GC), system is largely conserved across vertebrates and plays a central role in numerous vital physiological processes including bone development, immunomodulation, and modification of glucose metabolism and the induction of stress-related behaviours. As a result of their wide-ranging actions, synthetic GCs are widely prescribed for numerous human and veterinary therapeutic purposes and consequently have been detected extensively within the aquatic environment. Synthetic GCs designed for humans are pharmacologically active in non-mammalian vertebrates, including fish, however they are generally detected in surface waters at low (ng/L) concentrations. In this review, we assess the potential environmental risk of synthetic GCs to fish by comparing available experimental data and effect levels in fish with those in mammals. We found the majority of compounds were predicted to have insignificant risk to fish, however some compounds were predicted to be of moderate and high risk to fish, although the dataset of compounds used for this analysis was small. Given the common mode of action and high level of inter-species target conservation exhibited amongst the GCs, we also give due consideration to the potential for mixture effects, which may be particularly significant when considering the potential for environmental impact from this class of pharmaceuticals. Finally, we also provide recommendations for further research to more fully understand the potential environmental impact of this relatively understudied group of commonly prescribed human and veterinary drugs.

A comprehensive aquatic risk assessment of the beta-blocker propranolol, based on the results of over 600 research papers

A comprehensive aquatic environmental risk assessment (ERA) of the human pharmaceutical propranolol was conducted, based on all available scientific literature. Over 200 papers provided information on environmental concentrations (77 of which provided river concentrations) and 98 dealt with potential environmental effects. The median concentration of propranolol in rivers was 7.1 ng/L (range of median values of individual studies 0.07 to 89 ng/L), and the highest individual value was 590 ng/L. Sixty-eight EC50 values for 35 species were available. The lowest EC50 value was 0.084 mg/L. A species sensitivity distribution (SSD) provided an HC50 value of 6.64 mg/L and an HC5 value of 0.22 mg/L. Thus, there was a difference of nearly 6 orders of magnitude between the median river concentration and the HC50 value, and over 4 orders of magnitude between the median river concentration and the HC5 value. Even if an assessment factor of 100 was applied to the HC5 value, to provide considerable protection to all species, the safety margin is over 100-fold. However, nearly half of all papers reporting effects of propranolol did not provide an EC50 value. Some reported that very low concentrations of propranolol caused effects. The lowest concentration reported to cause an effect - in fact, a range of biochemical and physiological effects on mussels - was 0.3 ng/L. In none of these 'low concentration' papers was a sigmoidal concentration-response relationship obtained. Although inclusion of data from these papers in the ERA cause a change in the conclusion reached, we are sceptical of the repeatability of these 'low concentration' results. We conclude that concentrations of propranolol present currently in rivers throughout the world do not constitute a risk to aquatic organisms. We discuss the need to improve the quality of ecotoxicology research so that more robust ERAs acceptable to all stakeholders can be completed.

Effect of Water pH on the Uptake of Acidic (Ibuprofen) and Basic (Propranolol) Drugs in a Fish Gill Cell Culture Model

Water pH is predicted to affect the uptake of ionizable pharmaceuticals in fish. The current study used an in vitro primary fish gill cell culture system to assess the effect of pH values in the range of 4.5-8.75 on the uptake rates of the base propranolol (pKa 9.42) and the acid ibuprofen (pKa 4.59). The rate-limiting step in the uptake was the diffusive supply flux of the unionized form from the water to the apical membrane, with subsequent rapid transfer across the epithelium. Computed uptake rate based on the unionized fraction best described the uptake of propranolol and ibuprofen over the range of pH values 5-8 and 6-8.75, respectively. For ibuprofen, the computed uptake rate overestimated the uptake below pH 6 where the unionized fraction increased from 4% at pH 6 to 55% at pH 4.5. As the unionized fraction increased, the uptake rate plateaued suggesting a saturation of the transport process. For both drugs, large variations in the uptake occur with only small fluctuations in pH values. This occurs between pH values 6 and 8, which is the pH range acceptable in regulatory test guidelines and seen in most of our freshwaters.

Developing in vitro models to assess fish gill excretion of emerging contaminants

Advances in analytical methods have enabled the detection of emerging contaminants at ever lower concentrations in freshwaters. However, such measurements need to be linked to effect-based assays to identify risks. The bioconcentration factor (BCF) forms part of a chemical's environmental risk assessment (ERA), and current regulatory testing guidelines to calculate fish BCFs use hundreds of fish per chemical. Due to ethical concerns a reduction in the numbers of animals used is desired, and there is a need to identify in vitro or in silico alternatives which meet regulatory acceptance. This study describes the successful demonstration of a FIsh Gill Cell culture System (FIGCS) to assess an often overlooked parameter in pharmacokinetics: the excretion of drugs across the gill. The FIGCS tolerates the application of natural waters on its apical surface, mimicking the situation of the live fish, and thus in combination with advanced analytical methods, offers an opportunity to take lab-based testing used for ERA, such as compound uptake, biotransformation or excretion directly into field for validation with natural waters. Here we used the basic drug propranolol and the acidic ibuprofen as a demonstration of the FIGCS utility in three separate experiments. Excretion across the apical membrane showed saturation kinetics, suggesting the involvement of carrier-mediated processes. Both propranolol and ibuprofen were excreted across the epithelium from the media (internal blood equivalent) to the water, with ibuprofen excretion being considerably slower than propranolol excretion. Further studies indicate that ibuprofen may be complexing with fetal bovine serum (FBS) reducing bioavailability; in contrast propranolol efflux rate was unaffected, indicating that drugs behave differently in the presence of FBS and other plasma proteins. A key issue in future ERA is to better understand the effects of mixtures of different pollutant classes found in environmental samples, and this model offers an ethical path to do this.

Pharmacology-informed prediction of the risk posed to fish by mixtures of non-steroidal anti-inflammatory drugs (NSAIDs) in the environment

The presence of non-steroidal anti-inflammatory drugs (NSAIDs) in the aquatic environment has raised concern that chronic exposure to these compounds may cause adverse effects in wild fish populations. This potential scenario has led some stakeholders to advocate a stricter regulation of NSAIDs, especially diclofenac. Considering their global clinical importance for the management of pain and inflammation, any regulation that may affect patient access to NSAIDs will have considerable implications for public health. The current environmental risk assessment of NSAIDs is driven by the results of a limited number of standard toxicity tests and does not take into account mechanistic and pharmacological considerations. Here we present a pharmacology-informed framework that enables the prediction of the risk posed to fish by 25 different NSAIDs and their dynamic mixtures. Using network pharmacology approaches, we demonstrated that these 25 NSAIDs display a significant mechanistic promiscuity that could enhance the risk of target-mediated mixture effects near environmentally relevant concentrations. Integrating NSAIDs pharmacokinetic and pharmacodynamic features, we provide highly specific predictions of the adverse phenotypes associated with exposure to NSAIDs, and we developed a visual multi-scale model to guide the interpretation of the toxicological relevance of any given set of NSAIDs exposure data. Our analysis demonstrated a non-negligible risk posed to fish by NSAID mixtures in situations of high drug use and low dilution of waste-water treatment plant effluents. We anticipate that this predictive framework will support the future regulatory environmental risk assessment of NSAIDs and increase the effectiveness of ecopharmacovigilance strategies. Moreover, it can facilitate the prediction of the toxicological risk posed by mixtures via the implementation of mechanistic considerations and could be readily extended to other classes of chemicals.

Vision of a near future: Bridging the human health-environment divide. Toward an integrated strategy to understand mechanisms across species for chemical safety assessment

There is a growing recognition that application of mechanistic approaches to understand cross-species shared molecular targets and pathway conservation in the context of hazard characterization, provide significant opportunities in risk assessment (RA) for both human health and environmental safety. Specifically, it has been recognized that a more comprehensive and reliable understanding of similarities and differences in biological pathways across a variety of species will better enable cross-species extrapolation of potential adverse toxicological effects. Ultimately, this would also advance the generation and use of mechanistic data for both human health and environmental RA. A workshop brought together representatives from industry, academia and government to discuss how to improve the use of existing data, and to generate new NAMs data to derive better mechanistic understanding between humans and environmentally-relevant species, ultimately resulting in holistic chemical safety decisions. Thanks to a thorough dialogue among all participants, key challenges, current gaps and research needs were identified, and potential solutions proposed. This discussion highlighted the common objective to progress toward more predictive, mechanistically based, data-driven and animal-free chemical safety assessments. Overall, the participants recognized that there is no single approach which would provide all the answers for bridging the gap between mechanism-based human health and environmental RA, but acknowledged we now have the incentive, tools and data availability to address this concept, maximizing the potential for improvements in both human health and environmental RA.

Testing the Translational Power of the Zebrafish: An Interspecies Analysis of Responses to Cardiovascular Drugs

The zebrafish is rapidly emerging as a promising alternative in vivo model for the detection of drug-induced cardiovascular effects. Despite its increasing popularity, the ability of this model to inform the drug development process is often limited by the uncertainties around the quantitative relevance of zebrafish responses compared with nonclinical mammalian species and ultimately humans. In this test of concept study, we provide a comparative quantitative analysis of the in vivo cardiovascular responses of zebrafish, rat, dog, and human to three model compounds (propranolol, losartan, and captopril), which act as modulators of two key systems (beta-adrenergic and renin–angiotensin systems) involved in the regulation of cardiovascular functions. We used in vivo imaging techniques to generate novel experimental data of drug-mediated cardiovascular effects in zebrafish larvae. These data were combined with a database of interspecies mammalian responses (i.e., heart rate, blood flow, vessel diameter, and stroke volume) extracted from the literature to perform a meta-analysis of effect size and direction across multiple species. In spite of the high heterogeneity of study design parameters, our analysis highlighted that zebrafish and human responses were largely comparable in >80% of drug/endpoint combinations. However, it also revealed a high intraspecies variability, which, in some cases, prevented a conclusive interpretation of the drug-induced effect. Despite the shortcomings of our study, the meta-analysis approach, combined with a suitable data visualization strategy, enabled us to observe patterns of response that would likely remain undetected with more traditional methods of qualitative comparative analysis. We propose that expanding this approach to larger datasets encompassing multiple drugs and modes of action would enable a rigorous and systematic assessment of the applicability domain of the zebrafish from both a mechanistic and phenotypic standpoint. This will increase the confidence in its application for the early detection of adverse drug reactions in any major organ system.

ECOdrug: a database connecting drugs and conservation of their targets across species

Pharmaceuticals are designed to interact with specific molecular targets in humans and these targets generally have orthologs in other species. This provides opportunities for the drug discovery community to use alternative model species for drug development. It also means, however, there is potential for mode of action related effects in non-target wildlife species as many pharmaceuticals reach the environment through patient use and manufacturing wastes. Acquiring insight in drug target ortholog predictions across species and taxonomic groups has proven difficult because of the lack of an optimal strategy and because necessary information is spread across multiple and diverse sources and platforms. We introduce a new research platform tool, ECOdrug, that reliably connects drugs to their protein targets across divergent species. It harmonizes ortholog predictions from multiple sources via a simple user interface underpinning critical applications for a wide range of studies in pharmacology, ecotoxicology and comparative evolutionary biology. ECOdrug can be used to identify species with drug targets and identify drugs that interact with those targets. As such, it can be applied to support intelligent targeted drug safety testing by ensuring appropriate and relevant species are selected in ecological risk assessments. ECOdrug is freely accessible and available at: http://www.ecodrug.org.

Pharmacology beyond the patient - The environmental risks of human drugs

BACKGROUND

The presence of pharmaceuticals in the environment is a growing global concern and although environmental risk assessment is required for approval of new drugs in Europe and the USA, the adequacy of the current triggers and the effects-based assessments has been questioned.

OBJECTIVE

To provide a comprehensive analysis of all regulatory compliant aquatic ecotoxicity data and evaluate the current triggers and effects-based environmental assessments to facilitate the development of more efficient approaches for pharmaceuticals toxicity testing.

METHODS

Publicly-available regulatory compliant ecotoxicity data for drugs targeting human proteins was compiled together with pharmacological information including drug targets, Cmax and lipophilicity. Possible links between these factors and the ecotoxicity data for effects on, growth, mortality and/or reproduction, were evaluated. The environmental risks were then assessed based on a combined analysis of drug toxicity and predicted environmental concentrations based on European patient consumption data.

RESULTS

For most (88%) of the of 975 approved small molecule drugs targeting human proteins a complete set of regulatory compliant ecotoxicity data in the public domain was lacking, highlighting the need for both intelligent approaches to prioritize legacy human drugs for a tailored environmental risk assessment and a transparent database that captures environmental data. We show that presence/absence of drug-target orthologues are predictive of susceptible species for the more potent drugs. Drugs that target the endocrine system represent the highest potency and greatest risk. However, for most drugs (>80%) with a full set of ecotoxicity data, risk quotients assuming worst-case exposure assessments were below one in all European countries indicating low environmental risks for the endpoints assessed.

CONCLUSION

We believe that the presented analysis can guide improvements to current testing procedures, and provide valuable approaches for prioritising legacy drugs (i.e. those registered before 2006) for further ecotoxicity testing. For drugs where effects of possible concern (e.g. behaviour) are not captured in regulatory tests, additional mechanistic testing may be required to provide the highest confidence for avoiding environmental impacts.

The Use of Molecular Descriptors To Model Pharmaceutical Uptake by a Fish Primary Gill Cell Culture Epithelium

Modeling approaches such as quantitative structure-activity relationships (QSARs) use molecular descriptors to predict the bioavailable properties of a compound in biota. However, these models have mainly been derived based on empirical data for lipophilic neutral compounds and may not predict the uptake of ionizable compounds. The majority of pharmaceuticals are ionizable, and freshwaters can have a range of pH values that affect speciation. In this study, we assessed the uptake of 10 pharmaceuticals (acetazolamide, beclomethasone, carbamazepine, diclofenac, gemfibrozil, ibuprofen, ketoprofen, norethindrone, propranolol, and warfarin) with differing modes of action and physicochemical properties (p K_a, log S, log D, log K_ow, molecular weight (MW), and polar surface area (PSA)) by an in vitro primary fish gill cell culture system (FIGCS) for 24 h in artificial freshwater. Principal component analysis (PCA) and partial least-squares (PLS) regression was used to determine the molecular descriptors that influence the uptake rates. Ionizable drugs were taken up by FIGCS; a strong positive correlation was observed between log S and the uptake rate, and a negative correlation was observed between p K_a, log D, and MW and the uptake rate. This approach shows that models can be derived on the basis of the physicochemical properties of pharmaceuticals and the use of an in vitro gill system to predict the uptake of other compounds. There is a need for a robust and validated model for gill uptake that could be used in a tiered risk assessment to prioritize compounds for experimental testing.

Environmental risk assessment of metformin and its transformation product guanylurea: II. Occurrence in surface waters of Europe and the United States and derivation of predicted no-effect concentrations

Metformin (MET), CAS 1115-70-4 (Metformin hydrochloride), is an antidiabetic drug with high usage in North America and Europe and has become the subject of regulatory interest. A pharmaceutical industry working group investigated environmental risks of MET. Environmental fate and chronic effects data were collated across the industry for the present risk assessment. Predicted environmental concentrations (PECs) for MET were modeled for the USA and Europe using the PhATE and GREAT-ER models, respectively. PECs were compared with measured environmental concentrations (MECs) for the USA and Europe. A predicted no effect concentration (PNEC) of 1 mg/L for MET was derived by deterministic procedures, applying an assessment factor of 10 to the lowest no observed effect concentration (i.e., 10 mg/L) from multiple chronic studies with algae, daphnids and fish. The PEC/PNEC and MEC/PNEC risk characterization ratios were <1, indicating no significant risk for MET with high Margins of Safety (MOS) of >868. MET is known to degrade during wastewater treatment to guanylurea (GUU, CAS 141-83-3), which we have shown to further degrade. There are no GUU toxicity data in the literature; hence, chronic studies for GUU were conducted to derive a PNEC of 0.16 mg/L. PECs were derived for GUU as for MET, plus MECs were retrieved from the literature. The PEC/PNEC and MEC/PNEC risk characterization ratios for GUU were also <1, with an MOS of >6.5. Based on standard risk assessment procedures for both MET and its transformation product GUU, there is no significant risk to aquatic life.

Prediction of bioconcentration factors in fish and invertebrates using machine learning

The application of machine learning has recently gained interest from ecotoxicological fields for its ability to model and predict chemical and/or biological processes, such as the prediction of bioconcentration. However, comparison of different models and the prediction of bioconcentration in invertebrates has not been previously evaluated. A comparison of 24 linear and machine learning models is presented herein for the prediction of bioconcentration in fish and important factors that influenced accumulation identified. R² and root mean square error (RMSE) for the test data (n = 110 cases) ranged from 0.23-0.73 and 0.34-1.20, respectively. Model performance was critically assessed with neural networks and tree-based learners showing the best performance. An optimised 4-layer multi-layer perceptron (14 descriptors) was selected for further testing. The model was applied for cross-species prediction of bioconcentration in a freshwater invertebrate, Gammarus pulex. The model for G. pulex showed good performance with R² of 0.99 and 0.93 for the verification and test data, respectively. Important molecular descriptors determined to influence bioconcentration were molecular mass (MW), octanol-water distribution coefficient (logD), topological polar surface area (TPSA) and number of nitrogen atoms (nN) among others. Modelling of hazard criteria such as PBT, showed potential to replace the need for animal testing. However, the use of machine learning models in the regulatory context has been minimal to date and is critically discussed herein. The movement away from experimental estimations of accumulation to in silico modelling would enable rapid prioritisation of contaminants that may pose a risk to environmental health and the food chain.

Spheroid Size Does not Impact Metabolism of the β-blocker Propranolol in 3D Intestinal Fish Model

Compared to two-dimensional (2D) cell culture, cellular aggregates or spheroids (3D) offer a more appropriate alternative in vitro system where individual cell-cell communication and micro-environment more closely represent the in vivo organ; yet we understand little of the physiological conditions at this scale. The relationship between spheroid size and oxygen microenvironment, an important factor influencing the metabolic capacity of cells, was first established using the fish intestine derived RTgutGC cell line. Subsequently, pharmaceutical metabolism (Propranolol), as determined by high performance liquid chromatography, in this intestinal model was examined as a function of spheroid size. Co-efficient of variation between spheroid size was below 12% using the gyratory platform method, with the least variation observed in the highest cell seeding density. The viable, high oxygen micro-environment of the outer rim of the spheroid, as determined by electron paramagnetic resonance (EPR) oximetry, decreased over time, and the hypoxic zone increased as a function of spheroid size. Despite a trend of higher metabolism in smaller spheroids, the formation of micro-environments (quiescent, hypoxic or anoxic) did not significantly affect metabolism or function of an environmentally relevant pharmaceutical in this spheroid model.

A review of the pharmaceutical exposome in aquatic fauna

Pharmaceuticals have been considered 'contaminants of emerging concern' for more than 20 years. In that time, many laboratory studies have sought to identify hazard and assess risk in the aquatic environment, whilst field studies have searched for targeted candidates and occurrence trends using advanced analytical techniques. However, a lack of a systematic approach to the detection and quantification of pharmaceuticals has provided a fragmented literature of serendipitous approaches. Evaluation of the extent of the risk for the plethora of human and veterinary pharmaceuticals available requires the reliable measurement of trace levels of contaminants across different environmental compartments (water, sediment, biota - of which biota has been largely neglected). The focus on pharmaceutical concentrations in surface waters and other exposure media have therefore limited both the characterisation of the exposome in aquatic wildlife and the understanding of cause and effect relationships. Here, we compile the current analytical approaches and available occurrence and accumulation data in biota to review the current state of research in the field. Our analysis provides evidence in support of the 'Matthew Effect' and raises critical questions about the use of targeted analyte lists for biomonitoring. We provide six recommendations to stimulate and improve future research avenues.

Establishment and long-term maintenance of primary intestinal epithelial cells cultured from the rainbow trout, Oncorhynchus mykiss

A novel method for the establishment and long-term maintenance of ex vivo cultures from intestinal regions of the rainbow trout, Oncorhynchus mykiss (Walbaum), is reported. Adherence of cells was observed within hours, epithelial island formation recorded at 48 h and rapid proliferation with confluence achieved between 9-14 days. In addition to metabolic characterisation, basic morphology of growing cells was characterised using histology, immunofluorescence, transmission electron microscopy (TEM) and transepithelial electrical resistance (TEER). Regional differences in intestinal ethoxyresorufin-O-deethylase (EROD) and 7-ethoxycoumarin-O-deethylation (ECOD) activities in these primary grown enterocytes were compared following exposure to model inducers [i.e. α-NF, β-NF, B(a)P] which demonstrated significant differences. Regional differences in dietary uptake and metabolism of contaminants can therefore be studied in this in vitro system to increase our understanding of fundamental processes, while concurrently providing a means to reduce the number of fish required for biological studies in line with the principles of the 3Rs (Reduce, Refine and Replace).

This article has an associated First Person interview with the first author of the paper.

Summary: Understanding chemical uptake from the diet is difficult in live fish: we developed long-term intestinal cell cultures that enables the science and provides an alternative method.

Assessing the impact of benzo[a]pyrene with the in vitro fish gut model: An integrated approach for eco-genotoxicological studies

In vitro models are emerging tools for reducing reliance on traditional toxicity tests, especially in areas where information is sparse. For studies of fish, this is especially important for extrahepatic organs, such as the intestine, which, until recently, have been largely overlooked in favour of the liver or gill. Considering the importance of dietary uptake of contaminants, the rainbow trout (Oncorhynchus mykiss) intestine-derived cell line RTgutGC was cultured, to test its suitability as a high-throughput in vitro model. Benzo[a]pyrene (B[a]P) is an important contaminant and a model polycyclic aromatic hydrocarbon (PAH). Over 48 h exposure, a range of endpoints and xenobiotic metabolism rates were examined at three different pH levels indicative of the in vitro (pH 7.5) and in vivo mid-gut (pH 7.7) and hind-gut (pH 7.4) regions as a function of time. These endpoints included (i) cell viability: acid phosphatase (APH) and lactate dehydrogenase (LDH) assays; (ii) glucose uptake; (iii) cytochrome P450 enzyme activity: 7-ethoxyresoorufin-O-deethylase (EROD) assay; (iv) glutathione transferase (GST) activity; (v) genotoxic damage determined using the comet assay. Absence of cell viability loss, in parallel with decrease in the parent compound (B[a]P) in the medium and its subsequent increase in the cells suggested active sequestration, biotransformation, and removal of this representative PAH. With respect to genotoxic response, significant differences were observed at both the sampling times and the two highest concentrations of B[a]P. No significant differences were observed for the different pH conditions. Overall, this in vitro xenobiotic metabolism system appears to be a robust model, providing a basis for further development to evaluate metabolic and toxicological potential of contaminants without use of animals.

Investigations to extend viability of a rainbow trout primary gill cell culture

The primary culture of fish gill cells can provide functional, cell diverse, model in vitro platforms able to tolerate an aqueous exposure analogous to in vivo tissues. The utility of such models could be extended to a variety of longer term exposure scenarios if a method could be established to extend culture viability when exposed to water for longer periods. Here we report findings of a series of experiments to establish increased longevity, as monitored by culture transepithelial electrical resistance (TEER) and concurrent histological developments. Experimental cultures improved TEER during apical freshwater exposure for a mean of twelve days, compared to previous viabilities of up to 3 days. Cultures with larger surface areas and the use of trout serum rather than foetal bovine serum (FBS) contributed to the improvement, while perfusion of the intact gill prior to cell harvest resulted in a significantly faster preparation. Detailed scanning electron microscopy analysis of cultures revealed diverse surface structures that changed with culture age. Cultures grown on membranes with an increased porosity, collagen coating or 3D structure were of no benefit compared to standard membranes. Increased culture longevity, achieved in this study and reported for the first time, is a significant breakthrough and opens up a variety of future experimentation that has previously not been possible. The extended viability facilitates exploration of in vitro chronic or pulse-exposure test paradigms, longer term physiological and environmental monitoring studies and the potential for interactive co-culture with other organoid micro-tissues.

Application of the rainbow trout derived intestinal cell line (RTgutGC) for ecotoxicological studies: molecular and cellular responses following exposure to copper

There is an acknowledged need for in vitro fish intestinal model to help understand dietary exposure to chemicals in the aquatic environment. The presence and use of such models is however largely restrictive due to technical difficulties in the culturing of enterocytes in general and the availability of appropriate established cell lines in particular. In this study, the rainbow trout (Oncorhynchus mykiss) intestinal derived cell line (RTgutGC) was used as a surrogate for the "gut sac" method. To facilitate comparison, RTgutGC cells were grown as monolayers (double-seeded) on permeable Transwell supports leading to a two-compartment intestinal model consisting of polarised epithelium. This two-compartment model divides the system into an upper apical (lumen) and a lower basolateral (portal blood) compartment. In our studies, these cells stained weakly for mucosubstances, expressed the tight junction protein ZO-1 in addition to E-cadherin and revealed the presence of polarised epithelium in addition to microvilli protrusions. The cells also revealed a comparable transepithelial electrical resistance (TEER) to the in vivo situation. Importantly, the cell line tolerated apical saline (1:1 ratio) thus mimicking the intact organ to allow assessment of uptake of compounds across the intestine. Following an exposure over 72 h, our study demonstrated that the RTgutGC cell line under sub-lethal concentrations of copper sulphate (Cu) and modified saline solutions demonstrated uptake of the metal with saturation levels comparable to short term ex situ gut sac preparations. Gene expression analysis revealed no significant influence of pH or time on mRNA expression levels of key stress related genes (i.e. CYP3A, GST, mtA, Pgp and SOD) in the Transwell model. However, significant positive correlations were found between all genes investigated suggesting a co-operative relationship amongst the genes studied. When the outlined characteristics of the cell line are combined with the division of compartments, the RTgutGC double seeded model represents a potential animal replacement model for ecotoxicological studies. Overall, this model could be used to study the effects and predict aquatic gastrointestinal permeability of metals and other environmentally relevant contaminants in a cost effective and high throughput manner.

Uptake, biotransformation and elimination of selected pharmaceuticals in a freshwater invertebrate measured using liquid chromatography tandem mass spectrometry

Methods were developed to assess uptake and elimination kinetics in Gammarus pulex of nine pharmaceuticals (sulfamethazine, carbamazepine, diazepam, temazepam, trimethoprim, warfarin, metoprolol, nifedipine and propranolol) using targeted LC-MS/MS to determine bioconcentration factors (BCFs) using a 96 h toxicokinetic exposure and depuration period. The derived BCFs for these pharmaceuticals did not trigger any regulatory thresholds and ranged from 0 to 73 L kg^-1 (sulfamethazine showed no bioconcentration). Metabolism of chemicals can affect accurate BCF determination through parameterisation of the kinetic models. The added selectivity of LC-MS/MS allowed us to develop confirmatory methods to monitor the biotransformation of propranolol, carbamazepine and diazepam in G. pulex. Varying concentrations of the biotransformed products; 4-hydroxypropranolol sulphate, carbamazepine-10,11-epoxide, nordiazepam, oxazepam and temazepam were measured following exposure of the precursor compounds. For diazepam, the biotransformation product nordiazepam was present at higher concentrations than the parent compound at 94 ng g^-1 dw. Overall, the results indicate that pharmaceutical accumulation is low in these freshwater amphipods, which can potentially be explained by the rapid biotransformation and excretion.

Species specific anaesthetics for fish anaesthesia and euthanasia

There is a need to ensure that the care and welfare for fish maintained in the laboratory are to the highest standards. This extends to the use of anaesthetics for both scientific study, humane killing and euthanasia at end of life. An anaesthetic should not induce negative behaviours and fish should not seek to avoid the anaesthetic. Surprisingly little information is available to facilitate a humane choice of anaesthetic agent for fish despite over 100 years of use and the millions of fish currently held in thousands of laboratories worldwide. Using a chemotaxic choice chamber we found different species specific behavioural responses among four closely related fish species commonly held in the laboratory, exposed to three widely used anaesthetic agents. As previously found for zebrafish (Danio rerio), the use of MS-222 and benzocaine also appears to induce avoidance behaviours in medaka (Oryzias latipes); but etomidate could provide an alternative choice. Carp (Cyprinus carpio), although closely related to zebrafish showed avoidance behaviours to etomidate, but not benzocaine or MS-222; and rainbow trout (Oncorhynchus mykiss) showed no avoidance to the three agents tested. We were unable to ascertain avoidance responses in fathead minnows (Pimephales promelas) and suggest different test paradigms are required for that species.

The Role of Omics in the Application of Adverse Outcome Pathways for Chemical Risk Assessment

In conjunction with the second International Environmental Omics Symposium (iEOS) conference, held at the University of Liverpool (United Kingdom) in September 2014, a workshop was held to bring together experts in toxicology and regulatory science from academia, government and industry. The purpose of the workshop was to review the specific roles that high-content omics datasets (eg, transcriptomics, metabolomics, lipidomics, and proteomics) can hold within the adverse outcome pathway (AOP) framework for supporting ecological and human health risk assessments. In light of the growing number of examples of the application of omics data in the context of ecological risk assessment, we considered how omics datasets might continue to support the AOP framework. In particular, the role of omics in identifying potential AOP molecular initiating events and providing supportive evidence of key events at different levels of biological organization and across taxonomic groups was discussed. Areas with potential for short and medium-term breakthroughs were also discussed, such as providing mechanistic evidence to support chemical read-across, providing weight of evidence information for mode of action assignment, understanding biological networks, and developing robust extrapolations of species-sensitivity. Key challenges that need to be addressed were considered, including the need for a cohesive approach towards experimental design, the lack of a mutually agreed framework to quantitatively link genes and pathways to key events, and the need for better interpretation of chemically induced changes at the molecular level. This article was developed to provide an overview of ecological risk assessment process and a perspective on how high content molecular-level datasets can support the future of assessment procedures through the AOP framework.

4-dimensional functional profiling in the convulsant-treated larval zebrafish brain

Functional neuroimaging, using genetically-encoded Ca²⁺ sensors in larval zebrafish, offers a powerful combination of high spatiotemporal resolution and higher vertebrate relevance for quantitative neuropharmacological profiling. Here we use zebrafish larvae with pan-neuronal expression of GCaMP6s, combined with light sheet microscopy and a novel image processing pipeline, for the 4D profiling of chemoconvulsant action in multiple brain regions. In untreated larvae, regions associated with autonomic functionality, sensory processing and stress-responsiveness, consistently exhibited elevated spontaneous activity. The application of drugs targeting different convulsant mechanisms (4-Aminopyridine, Pentylenetetrazole, Pilocarpine and Strychnine) resulted in distinct spatiotemporal patterns of activity. These activity patterns showed some interesting parallels with what is known of the distribution of their respective molecular targets, but crucially also revealed system-wide neural circuit responses to stimulation or suppression. Drug concentration-response curves of neural activity were identified in a number of anatomically-defined zebrafish brain regions, and in vivo larval electrophysiology, also conducted in 4dpf larvae, provided additional measures of neural activity. Our quantification of network-wide chemoconvulsant drug activity in the whole zebrafish brain illustrates the power of this approach for neuropharmacological profiling in applications ranging from accelerating studies of drug safety and efficacy, to identifying pharmacologically-altered networks in zebrafish models of human neurological disorders.

International survey on the use and welfare of zebrafish Danio rerio in research

A survey was conducted regarding zebrafish Danio rerio use for scientific research with a focus on: anaesthesia and euthanasia; housing and husbandry; breeding and production; refinement opportunities. A total of 98 survey responses were received from laboratories in 22 countries in Europe, North America, South America, Asia and Australia. There appears a clear and urgent need to identify the most humane methods of anaesthesia and euthanasia. Aversive responses to MS-222 were widely observed raising concerns about the use of this anaesthetic for D. rerio. The use of anaesthesia in fin clipping for genetic identification is widely practised and there appears to be an opportunity to further develop less invasive methods and refine this process. Optimization (and potentially standardization) of feeding is an area for further investigation. Given that diet and body condition can have such profound effects on results of experiments, differences in practice could have significant scientific implications. Further research into transition between dark and light phases in the laboratory appears to represent an opportunity to establish best practice. Plants and gravel were not considered practical by many laboratories. The true value and benefits need to be established and communicated. Overproduction is a concern both from ethical and financial viewpoints. There is an opportunity to further reduce wastage of D. rerio. There are clear concerns and opportunities for the scientific community to work together to further improve the welfare of these important laboratory models.

Pharmaceutical Metabolism in Fish: Using a 3-D Hepatic In Vitro Model to Assess Clearance

At high internal doses, pharmaceuticals have the potential for inducing biological/pharmacological effects in fish. One particular concern for the environment is their potential to bioaccumulate and reach pharmacological levels; the study of these implications for environmental risk assessment has therefore gained increasing attention. To avoid unnecessary testing on animals, in vitro methods for assessment of xenobiotic metabolism could aid in the ecotoxicological evaluation. Here we report the use of a 3-D in vitro liver organoid culture system (spheroids) derived from rainbow trout to measure the metabolism of seven pharmaceuticals using a substrate depletion assay. Of the pharmaceuticals tested, propranolol, diclofenac and phenylbutazone were metabolised by trout liver spheroids; atenolol, metoprolol, diazepam and carbamazepine were not. Substrate depletion kinetics data was used to estimate intrinsic hepatic clearance by this spheroid model, which was similar for diclofenac and approximately 5 fold higher for propranolol when compared to trout liver microsomal fraction (S9) data. These results suggest that liver spheroids could be used as a relevant and metabolically competent in vitro model with which to measure the biotransformation of pharmaceuticals in fish; and propranolol acts as a reproducible positive control.

Testing the "read-across hypothesis" by investigating the effects of ibuprofen on fish

Human pharmaceuticals present in the environment have the potential to cause adverse effects on non-target organisms. The "read-across hypothesis" stipulates that pharmaceuticals will exhibit similar biological effects across species (e.g. human and fish) if the molecular target has been conserved and the effective drug concentrations are reached (Cmax). We tested this hypothesis by evaluating if ibuprofen, a non-selective inhibitor of prostaglandins and the cyclooxygenase (COX) enzyme, can mimic its primary effect in humans, on fish, at comparable plasma concentrations. The endpoints, prostaglandin E metabolite (PGEM) levels and the mRNA expression of COX (ptgs) gene, were measured in the gills of control and exposed fathead minnows (Pimephales promelas), using enzyme-immunoassay and quantitative real-time PCR (qPCR). Fish were exposed, for 24-72 h, to measured water concentrations of 9 (n = 12), 370 (n = 40) and 470 μg ibuprofen/L (n = 12). Water and blood plasma concentrations were determined using LC-MS/MS. Results showed that PGEM levels in fish exposed to 370 and 470 μg ibuprofen/L were significantly decreased compared to control fish, when mean plasma ibuprofen concentrations were 1.8-5.6-fold below the Cmax. The plasma ibuprofen concentrations and PGEM levels varied greatly between individuals. In fish exposed to 9 μg ibuprofen/L, when the mean plasma ibuprofen concentration was 224-fold below Cmax, no change in PGEM levels was observed. These data provide evidence for the read-across hypothesis, but suggest establishing a direct dose-response between internal plasma and PGEM is difficult, and would require significantly larger numbers of fish to overcome the inter-individual variation.

The First Attempt at Non-Linear in Silico Prediction of Sampling Rates for Polar Organic Chemical Integrative Samplers (POCIS)

Modeling and prediction of polar organic chemical integrative sampler (POCIS) sampling rates (Rs) for 73 compounds using artificial neural networks (ANNs) is presented for the first time. Two models were constructed: the first was developed ab initio using a genetic algorithm (GSD-model) to shortlist 24 descriptors covering constitutional, topological, geometrical and physicochemical properties and the second model was adapted for Rs prediction from a previous chromatographic retention model (RTD-model). Mechanistic evaluation of descriptors showed that models did not require comprehensive a priori information to predict Rs. Average predicted errors for the verification and blind test sets were 0.03 ± 0.02 L d(-1) (RTD-model) and 0.03 ± 0.03 L d(-1) (GSD-model) relative to experimentally determined Rs. Prediction variability in replicated models was the same or less than for measured Rs. Networks were externally validated using a measured Rs data set of six benzodiazepines. The RTD-model performed best in comparison to the GSD-model for these compounds (average absolute errors of 0.0145 ± 0.008 L d(-1) and 0.0437 ± 0.02 L d(-1), respectively). Improvements to generalizability of modeling approaches will be reliant on the need for standardized guidelines for Rs measurement. The use of in silico tools for Rs determination represents a more economical approach than laboratory calibrations.

Correction: Direct Measurements of Oxygen Gradients in Spheroid Culture System Using Electron Parametric Resonance Oximetry

[This corrects the article DOI: 10.1371/journal.pone.0149492.].

High-Content and Semi-Automated Quantification of Responses to Estrogenic Chemicals Using a Novel Translucent Transgenic Zebrafish

Rapid embryogenesis, together with genetic similarities with mammals, and the desire to reduce mammalian testing, are major incentives for using the zebrafish model in chemical screening and testing. Transgenic zebrafish, engineered for identifying target gene expression through expression of fluorophores, have considerable potential for both high-content and high-throughput testing of chemicals for endocrine activity. Here we generated an estrogen responsive transgenic zebrafish model in a pigment-free "Casper" phenotype, facilitating identification of target tissues and quantification of these responses in whole intact fish. Using the ERE-GFP-Casper model we show chemical type and concentration dependence for green fluorescent protein (GFP) induction and both spatial and temporal responses for different environmental estrogens tested. We also developed a semiautomated (ArrayScan) imaging and image analysis system that we applied to quantify whole body fluorescence responses for a range of different estrogenic chemicals in the new transgenic zebrafish model. The zebrafish model developed provides a sensitive and highly integrative system for identifying estrogenic chemicals, their target tissues and effect concentrations for exposures in real time and across different life stages. It thus has application for chemical screening to better direct health effects analysis of environmental estrogens and for investigating the functional roles of estrogens in vertebrates.

Individuals Maintain Similar Rates of Protein Synthesis over Time on the Same Plane of Nutrition under Controlled Environmental Conditions

Consistent individual differences in animal performance drive individual fitness under variable environmental conditions and provide the framework through which natural selection can operate. Underlying this concept is the assumption that individuals will display consistent levels of performance in fitness-related traits and interest has focused on individual variation and broad sense repeatability in a range of behavioural and physiological traits. Despite playing a central role in maintenance and growth, and with considerable inter-individual variation documented, broad sense repeatability in rates of protein synthesis has not been assessed. In this study we show for the first time that juvenile flounder Platichthys flesus reared under controlled environmental conditions on the same plane of nutrition for 46 days maintain consistent whole-animal absolute rates of protein synthesis (A_s). By feeding meals containing ¹⁵N-labelled protein and using a stochastic end-point model, two non-terminal measures of protein synthesis were made 32 days apart (d₁₄ and d₄₆). A_s values (mass-corrected to a standard mass of 12 g) showed 2- to 3-fold variation between individuals on d₁₄ and d₄₆ but individuals showed similar A_s values on both days with a broad sense repeatability estimate of 0.684 indicating significant consistency in physiological performance under controlled experimental conditions. The use of non-terminal methodologies in studies of animal ecophysiology to make repeat measures of physiological performance enables known individuals to be tracked across changing conditions. Adopting this approach, repeat measures of protein synthesis under controlled conditions will allow individual ontogenetic changes in protein metabolism to be assessed to better understand the ageing process and to determine individual physiological adaptive capacity, and associated energetic costs of adaptation, to global environmental change.

Assessing the reliability of uptake and elimination kinetics modelling approaches for estimating bioconcentration factors in the freshwater invertebrate, Gammarus pulex

This study considers whether the current standard toxicokinetic methods are an accurate and applicable assessment of xenobiotic exposure in an aquatic freshwater invertebrate. An in vivo exposure examined the uptake and elimination kinetics for eight pharmaceutical compounds in the amphipod crustacean, Gammarus pulex by measuring their concentrations in both biological material and in the exposure medium over a 96 h period. Selected pharmaceuticals included two anti-inflammatories (diclofenac and ibuprofen), two beta-blockers (propranolol and metoprolol), an anti-depressant (imipramine), an anti-histamine (ranitidine) and two beta-agonists (formoterol and terbutaline). Kinetic bioconcentration factors (BCFs) for the selected pharmaceuticals were derived from a first-order one-compartment model using either the simultaneous or sequential modelling methods. Using the simultaneous method for parameter estimation, BCF values ranged from 12 to 212. In contrast, the sequential method for parameter estimation resulted in bioconcentration factors ranging from 19 to 4533. Observed toxicokinetic plots showed statistically significant lack-of-fits and further interrogation of the models revealed a decreasing trend in the uptake rate constant over time for ranitidine, diclofenac, imipramine, metoprolol, formoterol and terbutaline. Previous published toxicokinetic data for 14 organic micro-pollutants were also assessed and similar trends were identified to those observed in this study. The decreasing trend of the uptake rate constant over time highlights the need to interpret modelled data more comprehensively to ensure uncertainties associated with uptake and elimination parameters for determining bioconcentration factors are minimised.

Internal exposure dynamics drive the Adverse Outcome Pathways of synthetic glucocorticoids in fish

The Adverse Outcome Pathway (AOP) framework represents a valuable conceptual tool to systematically integrate existing toxicological knowledge from a mechanistic perspective to facilitate predictions of chemical-induced effects across species. However, its application for decision-making requires the transition from qualitative to quantitative AOP (qAOP). Here we used a fish model and the synthetic glucocorticoid beclomethasone dipropionate (BDP) to investigate the role of chemical-specific properties, pharmacokinetics, and internal exposure dynamics in the development of qAOPs. We generated a qAOP network based on drug plasma concentrations and focused on immunodepression, skin androgenisation, disruption of gluconeogenesis and reproductive performance. We showed that internal exposure dynamics and chemical-specific properties influence the development of qAOPs and their predictive power. Comparing the effects of two different glucocorticoids, we highlight how relatively similar in vitro hazard-based indicators can lead to different in vivo risk. This discrepancy can be predicted by their different uptake potential, pharmacokinetic (PK) and pharmacodynamic (PD) profiles. We recommend that the development phase of qAOPs should include the application of species-specific uptake and physiologically-based PK/PD models. This integration will significantly enhance the predictive power, enabling a more accurate assessment of the risk and the reliable transferability of qAOPs across chemicals.

Direct Measurements of Oxygen Gradients in Spheroid Culture System Using Electron Parametric Resonance Oximetry

Advanced in vitro culture from tissues of different origin includes three-dimensional (3D) organoid micro structures that may mimic conditions in vivo. One example of simple 3D culture is spheroids; ball shaped structures typically used as liver and tumour models. Oxygen is critically important in physiological processes, but is difficult to quantify in 3D culture: and the question arises, how small does a spheroid have to be to have minimal micro-environment formation? This question is of particular importance in the growing field of 3D based models for toxicological assessment. Here, we describe a simple non-invasive approach modified for the quantitative measurement and subsequent evaluation of oxygen gradients in spheroids developed from a non-malignant fish cell line (i.e. RTG-2 cells) using Electron Paramagnetic Resonance (EPR) oximetry. Sonication of the paramagnetic probe Lithium phthalocyanine (LiPc) allows for incorporation of probe particulates into spheroid during its formation. Spectra signal strength after incorporation of probe into spheroid indicated that a volume of 20 μl of probe (stock solution: 0.10 mg/mL) is sufficient to provide a strong spectra across a range of spheroid sizes. The addition of non-toxic probes (that do not produce or consume oxygen) report on oxygen diffusion throughout the spheroid as a function of size. We provide evidence supporting the use of this model over a range of initial cell seeding densities and spheroid sizes with the production of oxygen distribution as a function of these parameters. In our spheroid model, lower cell seeding densities (∼2,500 cells/spheroid) and absolute size (118±32 μm) allow control of factors such as pre-existing stresses (e.g. ∼ 2% normoxic/hypoxic interface) for more accurate measurement of treatment response. The applied methodology provides an elegant, widely applicable approach to directly characterize spheroid (and other organoid) cultures in biomedical and toxicological research.

Pharmaceuticals in the freshwater invertebrate, Gammarus pulex, determined using pulverised liquid extraction, solid phase extraction and liquid chromatography-tandem mass spectrometry

The development, characterisation and application of a new analytical method for multi-residue PPCP determination in the freshwater amphipod, Gammarus pulex are presented. Analysis was performed using pulverised liquid extraction (PuLE), solid phase extraction (SPE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS). Qualitative method performance offered excellent limits of detection at <20 ng g(-1) for 18 out of 29 compounds. For quantitative application, linearity and precision were considered acceptable for 10 compounds across the ng-μg g(-1) range (R2≥0.99; ≤20% relative standard deviation respectively). The method was applied to the analysis of G. pulex and river water sourced from six tributaries of the River Thames. Carbamazepine, diazepam, nimesulide, trimethoprim and warfarin were determined in G. pulex samples at low ng g(-1) (dry weight) concentrations across these sites. Temazepam and diclofenac were also detected, but were not quantifiable. Six pharmaceuticals were quantified in surface waters across the eight sites at concentrations ranging from 3 to 344 ng L(-1). The possibility for confirmatory detection and subsequent quantification of pharmaceutical residues in benthic organisms such as G. pulex will enable further understanding on the susceptibility and ecological effects of PPCPs in the aquatic environment.

Effects of the lipid regulating drug clofibric acid on PPARα-regulated gene transcript levels in common carp (Cyprinus carpio) at pharmacological and environmental exposure levels

In mammals, the peroxisome proliferator-activated receptor α (PPARα) plays a key role in regulating various genes involved in lipid metabolism, bile acid synthesis and cholesterol homeostasis, and is activated by a diverse group of compounds collectively termed peroxisome proliferators (PPs). Specific PPs have been detected in the aquatic environment; however little is known on their pharmacological activity in fish. We investigated the bioavailability and persistence of the human PPARα ligand clofibric acid (CFA) in carp, together with various relevant endpoints, at a concentration similar to therapeutic levels in humans (20mg/L) and for an environmentally relevant concentration (4μg/L). Exposure to pharmacologically-relevant concentrations of CFA resulted in increased transcript levels of a number of known PPARα target genes together with increased acyl-coA oxidase (Acox1) activity, supporting stimulation of lipid metabolism pathways in carp which are known to be similarly activated in mammals. Although Cu,Zn-superoxide dismutase (Sod1) activity was not affected, mRNA levels of several biotransformation genes were also increased, paralleling previous reports in mammals and indicating a potential role in hepatic detoxification for PPARα in carp. Importantly, transcription of some of these genes (and Acox1 activity) were affected at exposure concentrations comparable with those reported in effluent discharges. Collectively, these data suggest that CFA is pharmacologically active in carp and has the potential to invoke PPARα-related responses in fish exposed in the environment, particularly considering that CFA may represent just one of a number of PPAR-active compounds present to which wild fish may be exposed.

Functional xenobiotic metabolism and efflux transporters in trout hepatocyte spheroid cultures

Prediction of xenobiotic fate in fish is important for the regulatory assessment of chemicals under current legislation.

Prediction of xenobiotic fate in fish is important for the regulatory assessment of chemicals under current legislation. Trout hepatocyte spheroids are a promising in vitro model for this assessment. In this investigation, the gene expression and function for xenobiotic metabolism and cellular efflux were characterised. Using fluorescence, transport and real time PCR analysis, the expression and functionality of a variety of genes related to xenobiotic metabolism and drug efflux were assessed in a range of trout hepatocyte culture preparations. Significantly greater levels of expression of genes involved in xenobiotic metabolism and efflux were measured in spheroids (which have been shown to remain viable in excess of 30 days), compared to hepatocytes cultured using conventional suspension and monolayer culture techniques. A transient decline in the expression of genes related to both xenobiotic metabolism and transport was determined during spheroid development, with a subsequent recovery in older spheroids. The most mature spheroids also exhibited an expression profile most comparable to that reported in vivo. Functionality of efflux transporters in spheroids was also demonstrated using fluorescent markers and specific inhibitors. In conclusion, the more physiologically relevant architecture in spheroid cultures provides a high functional integrity in relation to xenobiotic metabolism and efflux. Together with the enhanced gene expression and longevity of the model, hepatocytes in spheroid culture may prove to be an accurate alternative model to study the mechanisms of these processes in fish liver and provide an assay to determine the bioaccumulation potential of environmental contaminants.

A primary fish gill cell culture model to assess pharmaceutical uptake and efflux: evidence for passive and facilitated transport

The gill is the principle site of xenobiotic transfer to and from the aqueous environment. To replace, refine or reduce (3Rs) the large numbers of fish used in in vivo uptake studies an effective in vitro screen is required that mimics the function of the teleost gill. This study uses a rainbow trout (Oncorhynchus mykiss) primary gill cell culture system grown on permeable inserts, which tolerates apical freshwater thus mimicking the intact organ, to assess the uptake and efflux of pharmaceuticals across the gill. Bidirectional transport studies in media of seven pharmaceuticals (propranolol, metoprolol, atenolol, formoterol, terbutaline, ranitidine and imipramine) showed they were transported transcellularly across the epithelium. However, studies conducted in water showed enhanced uptake of propranolol, ranitidine and imipramine. Concentration-equilibrated conditions without a concentration gradient suggested that a proportion of the uptake of propranolol and imipramine is via a carrier-mediated process. Further study using propranolol showed that its transport is pH-dependent and at very low environmentally relevant concentrations (ng L(-1)), transport deviated from linearity. At higher concentrations, passive uptake dominated. Known inhibitors of drug transport proteins; cimetidine, MK571, cyclosporine A and quinidine inhibited propranolol uptake, whilst amantadine and verapamil were without effect. Together this suggests the involvement of specific members of SLC and ABC drug transporter families in pharmaceutical transport.

Quantitative cross-species extrapolation between humans and fish: the case of the anti-depressant fluoxetine

Fish are an important model for the pharmacological and toxicological characterization of human pharmaceuticals in drug discovery, drug safety assessment and environmental toxicology. However, do fish respond to pharmaceuticals as humans do? To address this question, we provide a novel quantitative cross-species extrapolation approach (qCSE) based on the hypothesis that similar plasma concentrations of pharmaceuticals cause comparable target-mediated effects in both humans and fish at similar level of biological organization (Read-Across Hypothesis). To validate this hypothesis, the behavioural effects of the anti-depressant drug fluoxetine on the fish model fathead minnow (Pimephales promelas) were used as test case. Fish were exposed for 28 days to a range of measured water concentrations of fluoxetine (0.1, 1.0, 8.0, 16, 32, 64 µg/L) to produce plasma concentrations below, equal and above the range of Human Therapeutic Plasma Concentrations (H_TPCs). Fluoxetine and its metabolite, norfluoxetine, were quantified in the plasma of individual fish and linked to behavioural anxiety-related endpoints. The minimum drug plasma concentrations that elicited anxiolytic responses in fish were above the upper value of the H_TPC range, whereas no effects were observed at plasma concentrations below the H_TPCs. In vivo metabolism of fluoxetine in humans and fish was similar, and displayed bi-phasic concentration-dependent kinetics driven by the auto-inhibitory dynamics and saturation of the enzymes that convert fluoxetine into norfluoxetine. The sensitivity of fish to fluoxetine was not so dissimilar from that of patients affected by general anxiety disorders. These results represent the first direct evidence of measured internal dose response effect of a pharmaceutical in fish, hence validating the Read-Across hypothesis applied to fluoxetine. Overall, this study demonstrates that the qCSE approach, anchored to internal drug concentrations, is a powerful tool to guide the assessment of the sensitivity of fish to pharmaceuticals, and strengthens the translational power of the cross-species extrapolation.

Bioavailability of the imidazole antifungal agent clotrimazole and its effects on key biotransformation genes in the common carp (Cyprinus carpio)

Clotrimazole (CTZ) is a persistent imidazole antifungal agent which is frequently detected in the aquatic environment and predicted to bio-concentrate in fish. Common carp (Cyprinus carpio) were exposed to mean measured concentrations of either 1.02 or 14.63μgl(-1) CTZ for 4 and 10 days, followed by a depuration period of 4 days in a further group of animals. Following each exposure regimen, plasma and liver CTZ concentrations were measured. Mean measured plasma concentrations of CTZ in animals exposed to the lower concentration of CTZ were 30 and 44μgl(-1) on days 4 and 10, respectively, and in the higher concentration were 318 and 336μgl(-1). Mean measured liver levels in the same animals were 514, 1725, 2111 and 7017μgl(-1) suggesting progressive hepatic accumulation. Measurement of CTZ in plasma after depuration suggested efficient elimination within 4 days, but appreciable levels of CTZ remained in the liver after depuration suggesting a degree of persistence in this tissue. In addition we measured responses of a number of key hepatic detoxification gene targets in the liver associated with the transcription factor pregnane X receptor (PXR); namely cyp450s 2k and 3a, glutathione-S-transferases a and p (gsta and p), and drug transporters multidrug resistance protein1 (mdr1), and MDR-related protein2 (mrp2). CTZ is a potent ligand of the PXR in humans and there is some evidence of PXR activation following exposure to CTZ in fish. The highest concentration of CTZ was adopted to explore the potential for alterations to detoxification gene expression in fish at a pharmacologically relevant dose level, and the lower concentration is within the range reported in effluents from waste water treatment works (WWTW). The genes for all biotransformation enzymes were up-regulated after exposure to the higher concentration of CTZ for 10 days, and alterations in expression occurred for the drug transporter genes mdr1 and mrp2 following exposure to the lower concentration of 1.02μgl(-1) CTZ (mean measured concentration). These data support the potential for CTZ to induce alterations in biotransformation and drug transporter genes associated with PXR in fish at concentrations measured in some WWTW effluents.

The read-across hypothesis and environmental risk assessment of pharmaceuticals

Pharmaceuticals in the environment have received increased attention over the past decade, as they are ubiquitous in rivers and waterways. Concentrations are in sub-ng to low μg/L, well below acute toxic levels, but there are uncertainties regarding the effects of chronic exposures and there is a need to prioritise which pharmaceuticals may be of concern. The read-across hypothesis stipulates that a drug will have an effect in non-target organisms only if the molecular targets such as receptors and enzymes have been conserved, resulting in a (specific) pharmacological effect only if plasma concentrations are similar to human therapeutic concentrations. If this holds true for different classes of pharmaceuticals, it should be possible to predict the potential environmental impact from information obtained during the drug development process. This paper critically reviews the evidence for read-across, and finds that few studies include plasma concentrations and mode of action based effects. Thus, despite a large number of apparently relevant papers and a general acceptance of the hypothesis, there is an absence of documented evidence. There is a need for large-scale studies to generate robust data for testing the read-across hypothesis and developing predictive models, the only feasible approach to protecting the environment.

Interactive effects of inbreeding and endocrine disruption on reproduction in a model laboratory fish

Inbreeding depression is expected to be more severe in stressful environments. However, the extent to which inbreeding affects the vulnerability of populations to environmental stressors, such as chemical exposure, remains unresolved. Here we report on the combined impacts of inbreeding and exposure to an endocrine disrupting chemical (the fungicide clotrimazole) on zebrafish (Danio rerio). We show that whilst inbreeding can negatively affect reproductive traits, not all traits are affected equally. Inbreeding depression frequently only became apparent when fish were additionally stressed by chemical exposure. Embryo viability was significantly reduced in inbred exposed fish and there was a tendency for inbred males to sire fewer offspring when in direct competition with outbred individuals. Levels of plasma 11-ketotestosterone, a key male sex hormone, showed substantial inbreeding depression that was unaffected by addition of the fungicide. In contrast, there was no effect of inbreeding or clotrimazole exposure on egg production. Overall, our data provide evidence that stress may amplify the effects of inbreeding on key reproductive traits, particularly those associated with male fitness. This may have important implications when considering the consequences of exposure to chemical pollutants on the fitness of wild populations.

Towards a more representative in vitro method for fish ecotoxicology: morphological and biochemical characterisation of three-dimensional spheroidal hepatocytes

The use of fish primary cells and cell lines offer an in vitro alternative for assessment of chemical toxicity and the evaluation of environmental samples in ecotoxicology. However, their uses are not without limitations such as short culture periods and loss of functionality, particularly with primary tissue. While three-dimensional (spheroid) technology is now established for in vitro mammalian toxicity studies, to date it has not been considered for environmental applications in a model aquatic species. In this study we report development of a reproducible six-well plate, gyratory-mediated method for rainbow trout (Oncorhynchus mykiss) hepatocyte spheroid culture and compare their functional and biochemical status with two-dimensional (2D) monolayer hepatocytes. Primary liver spheroid formation was divided into two stages, immature (1-5 days) and mature (≥6 days) according to size, shape and changes in functional and biochemical parameters (protein, glucose, albumin and lactate dehydrogenase). Mature spheroids retained the morphological characteristics (smooth outer surface, tight cell-cell contacts) previously described for mammalian spheroids as demonstrated by light and scanning electron microscopy. Glucose production and albumin synthesis were significantly higher in mature spheroids when compared to conventional 2D monolayer cultures (P < 0.01) and increased as spheroids matured (P < 0.01). Basal lactate dehydrogenase (LDH) leakage significantly decreased during spheroid formation and was significantly lower than 2D cultures (P < 0.01). It is therefore suggested that mature spheroids can maintain a high degree of functional, biochemical and morphological status over-time in culture that is superior to conventional 2D models and can provide realistic organotypic responses in vitro. Trout spheroids that take ~6-8 days to reach maturity would be suitable for use in acute toxicological tests and since it is possible to culture individual spheroids for over a month, there is potential for this work to lead towards in vitro bioaccumulation alternatives and to conduct high throughput screens of chronic exposure. This is an important step forward for developing alternative in vitro tools in future fish ecotoxicological studies.