Using in silico tools to predict flame retardant metabolites for more informative exposomics-based approaches.

Introduction: The positive identification of xenobiotics and their metabolites in human biosamples is an integral aspect of exposomics research, yet challenges in compound annotation and identification continue to limit the feasibility of comprehensive identification of total chemical exposure. Nonetheless, the adoption of in silico tools such as metabolite prediction software, QSAR-ready structural conversion workflows, and molecular standards databases can aid in identifying novel compounds in untargeted mass spectral investigations, permitting the assessment of a more expansive pool of compounds for human health hazard. This strategy is particularly applicable when it comes to flame retardant chemicals. The population is ubiquitously exposed to flame retardants, and evidence implicates some of these compounds as developmental neurotoxicants, endocrine disruptors, reproductive toxicants, immunotoxicants, and carcinogens. However, many flame retardants are poorly characterized, have not been linked to a definitive mode of toxic action, and are known to share metabolic breakdown products which may themselves harbor toxicity. As U.S. regulatory bodies begin to pursue a subclass- based risk assessment of organohalogen flame retardants, little consideration has been paid to the role of potentially toxic metabolites, or to expanding the identification of parent flame retardants and their metabolic breakdown products in human biosamples to better inform the human health hazards imposed by these compounds. Methods: The purpose of this study is to utilize publicly available in silico tools to 1) characterize the structural and metabolic fates of proposed flame retardant classes, 2) predict first pass metabolites, 3) ascertain whether metabolic products segregate among parent flame retardant classification patterns, and 4) assess the existing coverage in of these compounds in mass spectral database. Results: We found that flame retardant classes as currently defined by the National Academies of Science, Engineering and Medicine (NASEM) are structurally diverse, with highly variable predicted pharmacokinetic properties and metabolic fates among member compounds. The vast majority of flame retardants (96%) and their predicted metabolites (99%) are not present in spectral databases, posing a challenge for identifying these compounds in human biosamples. However, we also demonstrate the utility of publicly available in silico methods in generating a fit for purpose synthetic spectral library for flame retardants and their metabolites that have yet to be identified in human biosamples. Discussion: In conclusion, exposomics studies making use of fit-for-purpose synthetic spectral databases will better resolve internal exposure and windows of vulnerability associated with complex exposures to flame retardant chemicals and perturbed neurodevelopmental, reproductive, and other associated apical human health impacts.

Similarities and Differences in Gene Expression Networks Between the Breast Cancer Cell Line Michigan Cancer Foundation-7 and Invasive Human Breast Cancer Tissues.

Failure to adequately characterize cell lines, and understand the differences between in vitro and in vivo biology, can have serious consequences on the translatability of in vitro scientific studies to human clinical trials. This project focuses on the Michigan Cancer Foundation-7 (MCF-7) cells, a human breast adenocarcinoma cell line that is commonly used for in vitro cancer research, with over 42,000 publications in PubMed. In this study, we explore the key similarities and differences in gene expression networks of MCF-7 cell lines compared to human breast cancer tissues. We used two MCF-7 data sets, one data set collected by ARCHS4 including 1032 samples and one data set from Gene Expression Omnibus GSE50705 with 88 estradiol-treated MCF-7 samples. The human breast invasive ductal carcinoma (BRCA) data set came from The Cancer Genome Atlas, including 1212 breast tissue samples. Weighted Gene Correlation Network Analysis (WGCNA) and functional annotations of the data showed that MCF-7 cells and human breast tissues have only minimal similarity in biological processes, although some fundamental functions, such as cell cycle, are conserved. Scaled connectivity-a network topology metric-also showed drastic differences in the behavior of genes between MCF-7 and BRCA data sets. Finally, we used canSAR to compute ligand-based druggability scores of genes in the data sets, and our results suggested that using MCF-7 to study breast cancer may lead to missing important gene targets. Our comparison of the networks of MCF-7 and human breast cancer highlights the nuances of using MCF-7 to study human breast cancer and can contribute to better experimental design and result interpretation of study involving this cell line.

Mapping Chemical Respiratory Sensitization: How Useful Are Our Current Computational Tools?

Chemical respiratory sensitization is an immunological process that manifests clinically mostly as occupational asthma and is responsible for 1 in 6 cases of adult asthma, although this may be an underestimate of the prevalence, as it is under-diagnosed. Occupational asthma results in unemployment for roughly one-third of those affected due to severe health issues. Despite its high prevalence, chemical respiratory sensitization is difficult to predict, as there are currently no validated models and the mechanisms are not entirely understood, creating a significant challenge for regulatory bodies and industry alike. The Adverse Outcome Pathway (AOP) for respiratory sensitization is currently incomplete. However, some key events have been identified, and there is overlap with the comparatively well-characterized AOP for dermal sensitization. Because of this, and the fact that dermal sensitization is often assessed by in vivo, in chemico, or in silico methods, regulatory bodies are defaulting to the dermal sensitization status of chemicals as a proxy for respiratory sensitization status when evaluating chemical safety. We identified a data set of known human respiratory sensitizers, which we used to investigate the accuracy of a structural alert model, Toxtree, designed for skin sensitization and the Centre for Occupational and Environmental Health (COEH)'s model, a model developed specifically for occupational asthma. While both models had a reasonable level of accuracy, the COEH model achieved the highest balanced accuracy at 76%; when the models agreed, the overall accuracy was 87%. There were important differences between the models: Toxtree had superior performance for some structural alerts and some categories of well-characterized skin sensitizers, while the COEH model had high accuracy in identifying sensitizers that lacked identified skin sensitization reactivity domains. Overall, both models achieved respectable accuracy. However, neither model addresses potency, which, along with data quality, remains a hurdle, and the field must prioritize these issues to move forward.

Author Correction: Functionally Enigmatic Genes in Cancer: Using TCGA Data to Map the Limitations of Annotations.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Functionally Enigmatic Genes in Cancer: Using TCGA Data to Map the Limitations of Annotations.

Cancer is a comparatively well-studied disease, yet despite decades of intense focus, we demonstrate here using data from The Cancer Genome Atlas that a substantial number of genes implicated in cancer are relatively poorly studied. Those genes will likely be missed by any data analysis pipeline, such as enrichment analysis, that depends exclusively on annotations for understanding biological function. There is no indication that the amount of research - indicated by number of publications - is correlated with any objective metric of gene significance. Moreover, these genes are not missing at random but reflect that our information about genes is gathered in a biased manner: poorly studied genes are more likely to be primate-specific and less likely to have a Mendelian inheritance pattern, and they tend to cluster in some biological processes and not others. While this likely reflects both technological limitations as well as the fact that well-known genes tend to gather more interest from the research community, in the absence of a concerted effort to study genes in an unbiased way, many genes (and biological processes) will remain opaque.