In Silico Site-Directed Mutagenesis Informs Species-Specific Predictions of Chemical Susceptibility Derived From the Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) Tool

In vitro assessment of thyroid peroxidase inhibition by chemical exposure: comparison of cell models and detection methods

Disruption of the thyroid hormone (TH) system is connected with diverse adverse health outcomes in wildlife and humans. It is crucial to develop and validate suitable in vitro assays capable of measuring the disruption of the thyroid hormone (TH) system. These assays are also essential to comply with the 3R principles, aiming to replace the ex vivo tests often utilised in the chemical assessment. We compared the two commonly used assays applicable for high throughput screening [Luminol and Amplex UltraRed (AUR)] for the assessment of inhibition of thyroid peroxidase (TPO, a crucial enzyme in TH synthesis) using several cell lines and 21 compounds from different use categories. As the investigated cell lines derived from human and rat thyroid showed low or undetectable TPO expression, we developed a series of novel cell lines overexpressing human TPO protein. The HEK-TPOA7 model was prioritised for further research based on the high and stable TPO gene and protein expression. Notably, the Luminol assay detected significant peroxidase activity and signal inhibition even in Nthy-ori 3-1 and HEK293T cell lines without TPO expression, revealing its lack of specificity. Conversely, the AUR assay was specific to TPO activity. Nevertheless, despite the different specificity, both assays identified similar peroxidation inhibitors. Over half of the tested chemicals with diverse structures and from different use groups caused TPO inhibition, including some widespread environmental contaminants suggesting a potential impact of environmental chemicals on TH synthesis. Furthermore, in silico SeqAPASS analysis confirmed the high similarity of human TPO across mammals and other vertebrate classes, suggesting the applicability of HEK-TPOA7 model findings to other vertebrates.

The NIH Comparative Genomics Resource: addressing the promises and challenges of comparative genomics on human health

Comparative genomics is the comparison of genetic information within and across organisms to understand the evolution, structure, and function of genes, proteins, and non-coding regions (Sivashankari and Shanmughavel, Bioinformation 1:376-8, 2007). Advances in sequencing technology and assembly algorithms have resulted in the ability to sequence large genomes and provided a wealth of data that are being used in comparative genomic analyses. Comparative analysis can be leveraged to systematically explore and evaluate the biological relationships and evolution between species, aid in understanding the structure and function of genes, and gain a better understanding of disease and potential drug targets. As our knowledge of genetics expands, comparative genomics can help identify emerging model organisms among a broader span of the tree of life, positively impacting human health. This impact includes, but is not limited to, zoonotic disease research, therapeutics development, microbiome research, xenotransplantation, oncology, and toxicology. Despite advancements in comparative genomics, new challenges have arisen around the quantity, quality assurance, annotation, and interoperability of genomic data and metadata. New tools and approaches are required to meet these challenges and fulfill the needs of researchers. This paper focuses on how the National Institutes of Health (NIH) Comparative Genomics Resource (CGR) can address both the opportunities for comparative genomics to further impact human health and confront an increasingly complex set of challenges facing researchers.

Weight of evidence for cross-species conservation of androgen receptor-based biological activity

The U.S. Environmental Protection Agency's Endocrine Disruptor Screening Program (EDSP) is tasked with assessing chemicals for their potential to perturb endocrine pathways, including those controlled by androgen receptor (AR). To address challenges associated with traditional testing strategies, EDSP is considering in vitro high-throughput screening assays to screen and prioritize chemicals more efficiently. The ability of these assays to accurately reflect chemical interactions in nonmammalian species remains uncertain. Therefore, a goal of the EDSP is to evaluate how broadly results can be extrapolated across taxa. To assess the cross-species conservation of AR-modulated pathways, computational analyses and systematic literature review approaches were used to conduct a comprehensive analysis of existing in silico, in vitro, and in vivo data. First, molecular target conservation was assessed across 585 diverse species based on the structural similarity of ARs. These results indicate that ARs are conserved across vertebrates and are predicted to share similarly susceptibility to chemicals that interact with the human AR. Systematic analysis of over 5000 published manuscripts was used to compile in vitro and in vivo cross-species toxicity data. Assessment of in vitro data indicates conservation of responses occurs across vertebrate ARs, with potential differences in sensitivity. Similarly, in vivo data indicate strong conservation of the AR signaling pathways across vertebrate species, although sensitivity may vary. Overall, this study demonstrates a framework for utilizing bioinformatics and existing data to build weight of evidence for cross-species extrapolation and provides a technical basis for extrapolating hAR-based data to prioritize hazard in nonmammalian vertebrate species.

Demonstration of the Sequence Alignment to Predict Across Species Susceptibility Tool for Rapid Assessment of Protein Conservation

The US Environmental Protection Agency Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool is a fast, freely available, online screening application that allows researchers and regulators to extrapolate toxicity information across species. For biological targets in model systems such as human cells, mice, rats, and zebrafish, toxicity data are available for a variety of chemicals. Through the evaluation of protein target conservation, this tool can be used to extrapolate data generated from such model systems to thousands of other species lacking toxicity data, yielding predictions of relative intrinsic chemical susceptibility. The latest releases of the tool (versions 2.0-6.1) have incorporated new features that allow for the rapid synthesis, interpretation, and use of the data for publication plus presentation-quality graphics. Among these features are customizable data visualizations and a comprehensive summary report designed to summarize SeqAPASS data for ease of interpretation. This paper describes the protocol to guide users through submitting jobs, navigating the various levels of protein sequence comparisons, and interpreting and displaying the resulting data. New features of SeqAPASS v2.0-6.0 are highlighted. Furthermore, two use-cases focused on transthyretin and opioid receptor protein conservation using this tool are described. Finally, SeqAPASS' strengths and limitations are discussed to define the domain of applicability for the tool and highlight different applications for cross-species extrapolation.

Leveraging ToxCast data and protein sequence conservation to complement aquatic life criteria derivation

The USEPA's 1985 guidelines for the derivation of aquatic life criteria (ALC) are robust but data-intensive. For many chemicals, the extensive in vivo data sets required for ALC derivation are not available. Thus, alternative analyses and processes that can provide provisional values to guide states, tribes, and other stakeholders while data accumulate and more rigorous criteria are derived would be beneficial. The overarching purpose of this study was to assess the feasibility of using data from new approach methodologies (NAMs) like ToxCast to derive first-pass, provisional values to guide chemical prioritization and resource management as a complement to traditional ALC derivation. To address this goal, the study objectives were to (1) estimate chemical potency using data from NAMs for nine compounds with available aquatic benchmarks, (2) evaluate the utility of using NAM data to elucidate potential mechanisms of toxicity to guide problem formulation, and (3) determine the species relevance of toxicity pathways for compounds with clearly defined mechanisms of action as a means to evaluate whether minimum data requirements could potentially be waived when deriving a more formal ALC. Points of departure were derived from ToxCast data based on the fifth percentile of the distribution of activity concentration above cutoff values falling below the cytotoxic burst. Mechanistic inferences were made based on active target hits in ToxCast and, where applicable, assessed for taxonomic conservation using SeqAPASS. ToxCast-based point-of-departure aligned relatively closely (six of nine test chemicals within a factor of 10; eight of nine within a factor of 100) with aquatic benchmarks from the USEPA and US Department of Energy (DOE). Moreover, pathways of toxicity gleaned from NAM data were reflective of in vivo-based findings from the literature. These results, while preliminary, and based on a limited number of substances, support the potential application of NAM data to complement traditional ALC derivation approaches and prioritization. Integr Environ Assess Manag 2023;19:224-238. © 2022 Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Defining the Biologically Plausible Taxonomic Domain of Applicability of an Adverse Outcome Pathway: A Case Study Linking Nicotinic Acetylcholine Receptor Activation to Colony Death

For the majority of developed adverse outcome pathways (AOPs), the taxonomic domain of applicability (tDOA) is typically narrowly defined with a single or a handful of species. Defining the tDOA of an AOP is critical for use in regulatory decision-making, particularly when considering protection of untested species. Structural and functional conservation are two elements that can be considered when defining the tDOA. Publicly accessible bioinformatics approaches, such as the Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool, take advantage of existing and growing databases of protein sequence and structural information to provide lines of evidence toward structural conservation of key events (KEs) and KE relationships (KERs) of an AOP. It is anticipated that SeqAPASS results could readily be combined with data derived from empirical toxicity studies to provide evidence of both structural and functional conservation, to define the tDOA for KEs, KERs, and AOPs. Such data could be incorporated in the AOP-Wiki as lines of evidence toward biological plausibility for the tDOA. We present a case study describing the process of using bioinformatics to define the tDOA of an AOP using an AOP linking the activation of the nicotinic acetylcholine receptor to colony death/failure in Apis mellifera. Although the AOP was developed to gain a particular biological understanding relative to A. mellifera health, applicability to other Apis bees, as well as non-Apis bees, has yet to be defined. The present study demonstrates how bioinformatics can be utilized to rapidly take advantage of existing protein sequence and structural knowledge to enhance and inform the tDOA of KEs, KERs, and AOPs, focusing on providing evidence of structural conservation across species. Environ Toxicol Chem 2023;42:71-87. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Integrative Computational Approaches to Inform Relative Bioaccumulation Potential of Per- and Polyfluoroalkyl Substances Across Species

Predictive toxicology is increasingly reliant on innovative computational methods to address pressing questions in chemicals assessment. Of importance is the evaluation of contaminant impact differences across species to inform ecosystem protection and identify appropriate model species for human toxicity studies. Here we evaluated 2 complementary tools to predict cross-species differences in binding affinity between per- and polyfluoroalkyl substances (PFAS) and the liver fatty acid-binding protein (LFABP): the Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool and molecular dynamics (MD). SeqAPASS determined that the structure of human LFABP, a key determinant of PFAS bioaccumulation, was conserved in the majority of vertebrate species, indicating these species would have similar PFAS bioaccumulation potentials. Level 3 SeqAPASS evaluation identified several potentially destabilizing amino acid differences across species, which were generally supported by DUET stability change predictions. Nine single-residue mutations and 7 whole species sequences were selected for MD evaluation. One mutation (F50V for PFNA) showed a statistically significant difference with stronger affinity than wild-type human LFABP. Predicted binding affinities for 9 different PFAS across 7 species showed human, rat, chicken, and rainbow trout had similar binding affinities to one another for each PFAS, whereas Japanese medaka and fathead minnow had significantly weaker LFABP-binding affinity for some PFAS. Based on these analyses, the combined use of SeqAPASS and MD provides rapid screening for potential species differences with deeper structural insight. This approach can be easily extended to other important biological receptors and potential ligands.

Quantitative Chemical Proteomics Reveals Interspecies Variations on Binding Schemes of L-FABP with Perfluorooctanesulfonate

Evaluating interspecies toxicity variation is a long-standing challenge for chemical hazard assessment. This study developed a quantitative interspecies thermal shift assay (QITSA) for in situ, quantitative, and modest-throughput investigation of chemical-protein interactions in cell and tissue samples across species. By using liver fatty acid binding protein (L-FABP) as a case study, the QITSA method was benchmarked with six per- and polyfluoroalkyl substances, and thermal shifts (ΔT_m) were inversely related to their dissociation constants (R² = 0.98). The QITSA can also distinguish binding modes of chemicals exemplified by palmitic acid. The QITSA was applied to determine the interactions between perfluorooctanesulfonate (PFOS) and L-FABP in liver cells or tissues from humans, mice, rats, and zebrafish. The largest thermal stability enhancement by PFOS was observed for human L-FABP followed by the mouse, rat, and zebrafish. While endogenous ligands were revealed to partially contribute to the large interspecies variation, recombinant proteins were employed to confirm the high binding affinity of PFOS to human L-FABP, compared to the rat and mouse. This study implemented an experimental strategy to characterize chemical-protein interactions across species, and future application of QITSA to other chemical contaminants is of great interest.

Off-Target Stoichiometric Binding Identified from Toxicogenomics Explains Why Some Species Are More Sensitive than Others to a Widely Used Neonicotinoid

Neonicotinoids are currently licensed for use in 120 countries, making accurate nontarget species sensitivity predictions critical. Unfortunately, such predictions are fraught with uncertainty, as sensitivity is extrapolated from only a few test species and neonicotinoid sensitivities can differ greatly between closely related taxa. Combining classical toxicology with de novo toxicogenomics could greatly improve sensitivity predictions and identify unexpectedly susceptible species. We show that there is a >30-fold differential species sensitivity (DSS) for the neonicotinoid imidacloprid between five earthworm species, a critical nontarget taxon. This variation could not be explained by differential toxicokinetics. Furthermore, comparing key motif expression in subunit genes of the classical nicotinic acetylcholine receptor (nAChR) target predicts only minor differences in the ligand binding domains (LBDs). In contrast, predicted dissimilarities in LBDs do occur in the highly expressed but nonclassical targets, acetylcholine binding proteins (AChBPs). Critically, the predicted AChBP divergence is capable of explaining DSS. We propose that high expression levels of putative nonsynaptic AChBPs with high imidacloprid affinities reduce imidacloprid binding to critical nAChRs involved in vital synaptic neurotransmission. This study provides a clear example of how pragmatic interrogation of key motif expression in complex multisubunit receptors can predict observed DSS, thereby informing sensitivity predictions for essential nontarget species.

Conversion of Estrone to 17β-Estradiol: A Potential Confounding Factor in Assessing Risks of Environmental Estrogens to Fish

Feminization of male fish and the role of endocrine-active chemicals in this phenomenon has been an area of intense study for many years. Estrone (E1), a natural steroid, is found in aquatic environments sometimes at high concentrations relative to the estrogenic steroids 17β-estradiol (E2) and 17α-ethynylestradiol. However, E1 has been less thoroughly studied than E2 or 17α-ethynylestradiol due in part to a relatively lower potency in metabolically limited estrogen receptor (ER) binding/activation assays. Recent evidence suggests that in vivo biotransformation of E1 to E2 may occur in fathead minnows (Pimephales promelas) residing in environments with high concentrations of E1, such as near wastewater treatment plants. The enzymes likely responsible for this biotransformation, 17β-hydroxysteroid dehydrogenases (17βHSDs), have been well characterized in mammals but to a lesser extent in fish species. In the present study, a novel systematic analysis of amino acid sequence data from the National Center for Biotechnology Information database demonstrated that multiple 17βHSD isoforms are conserved across different fish species. Experimentally, we showed that metabolically active hepatic cytosolic preparations from 2 commercially important salmonid species, rainbow trout and lake trout, biotransformed E1 to E2 to a degree sufficient to alter results of competitive ER binding assays. These results from in silico and in vitro analyses indicate that E1 and biotransformation may play a significant role in adverse effects on development and reproduction of a variety of fish species in contaminated aquatic environments. Environ Toxicol Chem 2020;39:2028-2040. Published 2020. This article is a US Government work and is in the public domain in the USA.

Neuropathy target esterase (NTE/PNPLA6) and organophosphorus compound-induced delayed neurotoxicity (OPIDN)

Systemic inhibition of neuropathy target esterase (NTE) with certain organophosphorus (OP) compounds produces OP compound-induced delayed neurotoxicity (OPIDN), a distal degeneration of axons in the central nervous system (CNS) and peripheral nervous system (PNS), thereby providing a powerful model for studying a spectrum of neurodegenerative diseases. Axonopathies are important medical entities in their own right, but in addition, illnesses once considered primary neuronopathies are now thought to begin with axonal degeneration. These disorders include Alzheimer's disease, Parkinson's disease, and motor neuron diseases such as amyotrophic lateral sclerosis (ALS). Moreover, conditional knockout of NTE in the mouse CNS produces vacuolation and other degenerative changes in large neurons in the hippocampus, thalamus, and cerebellum, along with degeneration and swelling of axons in ascending and descending spinal cord tracts. In humans, NTE mutations cause a variety of neurodegenerative conditions resulting in a range of deficits including spastic paraplegia and blindness. Mutations in the Drosophila NTE orthologue SwissCheese (SWS) produce neurodegeneration characterized by vacuolization that can be partially rescued by expression of wild-type human NTE, suggesting a potential therapeutic approach for certain human neurological disorders. This chapter defines NTE and OPIDN, presents an overview of OP compounds, provides a rationale for NTE research, and traces the history of discovery of NTE and its relationship to OPIDN. It then briefly describes subsequent studies of NTE, including practical applications of the assay; aspects of its domain structure, subcellular localization, and tissue expression; abnormalities associated with NTE mutations, knockdown, and conventional or conditional knockout; and hypothetical models to help guide future research on elucidating the role of NTE in OPIDN.

Homology models of mouse and rat estrogen receptor-α ligand-binding domain created by in silico mutagenesis of a human template: molecular docking with 17ß-estradiol, diethylstilbestrol, and paraben analogs

Crystal structures exist for human, but not rodent, estrogen receptor-α ligand-binding domain (ERα-LBD). Consequently, rodent studies involving binding of compounds to ERα-LBD are limited in their molecular-level interpretation and extrapolation to humans. Because the sequences of rodent and human ERα-LBDs are > 95% identical, we expected their 3D structures and ligand binding to be highly similar. To test this hypothesis, we used the human ERα-LBD structure (PDB 3UUD) as a template to produce rat and mouse homology models. Employing the rodent models and human structure, we generated docking poses of 23 Group A ligands (17ß-estradiol, diethylstilbestrol, and 21 paraben analogs) in AutoDock Vina for interspecies comparisons. Ligand RMSDs (Å) (median, 95% CI) were 0.49 (0.21-1.82) (human-mouse) and 1.19 (0.22-1.82) (human-rat), well below the 2.0-2.5 Å range for equivalent docking poses. Numbers of interspecies ligand-receptor residue contacts were highly similar, with Sorensen Sc (%) = 96.8 (90.0-100) (human-mouse) and 97.7 (89.5-100) (human-rat). Likewise, numbers of interspecies ligand-receptor residue contacts were highly correlated: Pearson r = 0.913 (human-mouse) and 0.925 (human-rat). Numbers of interspecies ligand-receptor atom contacts were even more tightly correlated: r = 0.979 (human-mouse) and 0.986 (human-rat). Pyramid plots of numbers of ligand-receptor atom contacts by residue exhibited high interspecies symmetry and had Spearman r s = 0.977 (human-mouse) and 0.966 (human-rat). Group B ligands included 15 ring-substituted parabens recently shown experimentally to exhibit decreased binding to human ERα and to exert increased antimicrobial activity. Ligand efficiencies calculated from docking ligands into human ERα-LBD were well correlated with those derived from published experimental data (Pearson partial r p = 0.894 and 0.918; Groups A and B, respectively). Overall, the results indicate that our constructed rodent ERα-LBDs interact with ligands in like manner to the human receptor, thus providing a high level of confidence in extrapolations of rodent to human ligand-receptor interactions.

A unified approach for protecting listed species and ecosystem services in isolated wetlands using community-level protection goals

The protection of listed species through the Ecological Risk Assessment (ERA) process is encumbered by the number and diversity of species that need protection and the limited data available to inform assessments. Ecological communities within isolated ecosystems often contain a number of biologically diverse endemic, endangered, and threatened species, as well as providing numerous ecosystem services (ES). We propose an approach that develops community-level protection goals using isolated wetlands that includes both listed species and Service Providing Units (SPUs) that drive ES for ecological risk assessments (ERAs). Community-level protection goals are achieved by developing a protection community and weighing lines of evidence to determine a set of focal species within that community upon which to base the assessment. Lines of evidence include chemical mechanism of action, likely routes of exposure, and taxa susceptibility, as well as relationships among species, and other ecological factors. We demonstrate the process using case studies of chlorpyrifos in California vernal pools and coal ash effluent in Carolina bays. In the California vernal pool case study, listed species were the primary SPUs for the ES provided by the critical habitat. The weight of evidence demonstrated the honey bee as the focal species for the terrestrial environment and the vernal pool fairy shrimp as the focal species for the aquatic environment. The protection community within the Carolina bay case study was more taxonomically diverse than vernal pools for both listed species and SPUs, with amphibians identified as the focal species for which to target mitigation goals and hazard levels. The approach presented here will reduce the time and resource investment required for assessment of risk to listed species and adds an ES perspective to demonstrate value of assessments beyond listed species concerns.

Evidence for Cross Species Extrapolation of Mammalian-Based High-Throughput Screening Assay Results

High-throughput screening (HTS) and computational technologies have emerged as important tools for chemical hazard identification. The US Environmental Protection Agency (EPA) launched the Toxicity ForeCaster (ToxCast) Program, which has screened thousands of chemicals in hundreds of mammalian-based HTS assays for biological activity. The data are being used to prioritize toxicity testing on those chemicals likely to lead to adverse effects. To use HTS assays in predicting hazard to both humans and wildlife, it is necessary to understand how broadly these data may be extrapolated across species. The US EPA Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS; https://seqapass.epa.gov/seqapass/ ) tool was used to assess conservation of the 484 protein targets represented in the suite of ToxCast assays and other HTS assays. To demonstrate the utility of the SeqAPASS data for guiding extrapolation, case studies were developed which focused on targets of interest to the US Endocrine Disruptor Screening Program and the Organisation for Economic Cooperation and Development. These case studies provide a line of evidence for conservation of endocrine targets across vertebrate species, with few exceptions, and demonstrate the utility of SeqAPASS for defining the taxonomic domain of applicability for HTS results and identifying organisms for suitable follow-up toxicity tests.