An in silico structural approach to characterize human and rainbow trout estrogenicity of mycotoxins: Proof of concept study using zearalenone and alternariol

Human and fish differences in steroid receptors activation: A review

Steroid receptors (SRs) are transcription factors activated by steroid hormones (SHs) that belong to the nuclear receptors (NRs) superfamily. Several studies have shown that SRs are targets of endocrine disrupting chemicals (EDCs), widespread substances in the environment capable of interfering with the endogenous hormonal pathways and causing adverse health effects in living organisms and/or their progeny. Cell lines with SRs reporter gene are currently used for in vitro screening of large quantities of chemicals with suspected endocrine-disrupting activities. However, most of these cell lines express human SRs and therefore the toxicological data obtained are also extrapolated to non-mammalian species. In parallel, in vivo tests have recently been developed on fish species whose data are also extrapolated to mammalian species. As some species-specific differences in SRs activation by natural and synthetic chemicals have been recently reported, the aim of this review is to summarize those between human and fish SRs, as representatives of mammalian and non-mammalian toxicology, respectively. Overall, this literature study aims to improve inter-species extrapolation of toxicological data on EDCs and to understand which reporter gene cell lines expressing human SRs are relevant for the assessment of effects in fish and whether in vivo tests on fish can be properly used in the assessment of adverse effects on human health.

Free fatty acid receptors beyond fatty acids: A computational journey to explore peptides as possible binders of GPR120

Free fatty acids receptors, with members among G protein-coupled receptors (GPCRs), are crucial for biological signaling, including the perception of the so called "fatty taste". In recent years, GPR120, a protein belonging to the GPCR family, drew attention as an interesting pharmacological target to cope with obesity, satiety and diabetes. Apart from long chain fatty acids, which are GPR120 natural agonists, other synthetic molecules were identified as agonists expanding the chemical space of GPR120's ligands. In this scenario, we unveiled peptides as possible GPR120 binders toward a better understanding of this multifaceted and relevant target. This study analyzed a virtual library collecting 531 441 low-polar hexapeptides, providing mechanistic insights on the GPR120 activation and further extending the possible chemical space of GPR120 agonists. The computational pipeline started with a narrow filtering of hexapeptides based on their chemical similarity with known GPR120 agonists. The best hits were tested through docking studies, molecular dynamics and umbrella sampling simulations, which pointed to G[I,L]FGGG as a promising GPR120 agonist sequence. The presence of both peptides in food-related proteins was thoroughly assessed, revealing they may occur in mushrooms, food-grade bacteria and rice. Simulations on the counterparts with D-amino acids were also performed. Umbrella sampling simulations described that GdIFGGG may have a better interaction compared to its all-L counterpart (-13 kCal/mol ΔG and -6 kCal/mol ΔG, respectively). Overall, we obtained a predictive model to better understand the underpinning mechanism of GPR120-hexapeptides interaction, hierarchizing novel potential agonist peptides for further analysis and describing promising food sources worth of further dedicated investigations.

The bitter side of toxicity: A big data analysis spotted the interaction between trichothecenes and bitter receptors

The bitter taste perception evolved in human and animals to rapidly perceive and avoid potential toxic compounds. This is mediated by taste receptors type 2 (TAS2R), expressed in various tissues, which recently proved to be involved in roles beyond the bitter perception itself. With this study, the interaction between food-related toxic compounds and TAS2R46 has been investigated via computational approaches, starting with a virtual screening and moving to molecular docking and dynamics simulations. The virtual screening analysis identified trichothecolone and the trichothecenes class it belongs to, which includes mycotoxins widespread in several commodities raising food safety concerns, as possible TAS2R46 binders. Molecular docking and dynamics simulations were performed to further explore the trichotecenes-TAS2R46 interaction. The results indicated that deoxynivalenol and its 15-acetylated derivative could activate TAS2R46. Eventually, this study provided initial evidence supporting the involvement of TAS2R46 in the underpinning mechanisms of deoxynivalenol action highlighting the need of digging into the involvement of TAS2R46 and TAS2Rs in the adverse effects of deoxynivalenol and congeners.

Using In Silico Approach for Metabolomic and Toxicity Prediction of Alternariol

Alternariol is a metabolite produced by Alternaria fungus that can contaminate a variety of food and feed materials. The objective of the present paper was to provide a prediction of Phase I and II metabolites of alternariol and a detailed ADME/Tox profile for alternariol and its metabolites using an in silico working model based on the MetaTox, SwissADME, pKCMS, and PASS online computational programs. A number of 12 metabolites were identified as corresponding to the metabolomic profile of alternariol. ADME profile for AOH and predicted metabolites indicated a moderate or high intestinal absorption probability but a low probability to penetrate the blood-brain barrier. In addition to cytotoxic, mutagenic, carcinogenic, and endocrine disruptor effects, the computational model has predicted other toxicological endpoints for the analyzed compounds, such as vascular toxicity, haemato-toxicity, diarrhea, and nephrotoxicity. AOH and its metabolites have been predicted to act as a substrate for different isoforms of phase I and II drug-metabolizing enzymes and to interact with the response to oxidative stress. In conclusion, in silico methods can represent a viable alternative to in vitro and in vivo tests for the prediction of mycotoxins metabolism and toxicity.

A Computational Inter-Species Study on Safrole Phase I Metabolism-Dependent Bioactivation: A Mechanistic Insight into the Study of Possible Differences among Species

Safrole, a 162.2 Da natural compound belonging to the alkenylbenzenes class, is classified as a possible carcinogen to humans by IARC (group IIB) and has proven to be genotoxic and carcinogenic to rodents. Despite its use as a food or feed additive, it is forbidden in many countries due to its documented toxicity; yet, it is still broadly present within food and feed and is particularly abundant in spices, herbs and essential oils. Specifically, safrole may exert its toxicity upon bioactivation to its proximate carcinogen 1'-hydroxy-safrole via specific members of the cytochrome P450 protein family with a certain inter/intra-species variability. To investigate this variability, an in-silico workflow based on molecular modelling, docking and molecular dynamics has been successfully applied. This work highlighted the mechanistic basis underpinning differences among humans, cats, chickens, goats, sheep, dogs, mice, pigs, rats and rabbits. The chosen metric to estimate the likeliness of formation of 1'-hydroxy-safrole by the species-specific cytochrome P450 under investigation allowed for the provision of a knowledge-based ground to rationally design and prioritise further experiments and deepen the current understanding of alkenylbenzenes bioactivation and CYPs mechanics. Both are crucial for a more informed framework of analysis for safrole toxicity.

Deciphering the origin of total estrogenic activity of complex mixtures

Introduction

Identifying compounds with endocrine properties in food is getting increasingly important. Current chemical analysis methodology is mainly focused on the identification of known substances without bringing insight for biological activity. Recently, the application of bioassays has been promoted for their potential to detect unknown bioactive substances and to provide information on possible interactions between molecules. From the toxicological perspective, measuring endocrine activity cannot inform on endocrine disruption and/or health risks without sufficient knowledge on the nature of the responsible factors.

Methods

The present study addresses a promising approach using High Performance Thin-Layer Chromatography (HPTLC) coupled to bioassays were analyzed using the Liquid Chromatography Mass-Spectrometry (LC-MS). The estrogen receptor activation was assessed using the transcription activation Estrogen Receptor Alpha Chemical Activated LUciferase gene eXpression assay (ERα- CALUX) and the HPTLC coupled to the Estrogen Screen Yeast assay (p-YES).

Results

Seven isoflavones were identified in the soy isolates. Estrogen receptor activation was assessed for both, the identified isoflavones and the soy isolates with ERα-CALUX test. Correlation between the soy isolates extracts and the identified isoflavones was shown. Moreover, p-YES revealed the presence of an estrogenic bioactive zone. Analysis of the bioactive zone through LCHRMS highlighted signals corresponding to several isoflavones already detected in the isolates as well as two additional ones. For all detected isoflavones, an estrogenic activity dose-response was established in both bioassays.

Conclusion

Finally, genistein, daidzein, and naringenin were found as the most active substances. A concordance analysis integrating the analytical and bioassay data indicated that genistein and daidzein were the drivers of the estrogenic activity of these soy protein isolates. Altogether, these data suggest that the integration of HPTLC-bioassay together with chemical analysis is a powerful approach to characterize the endocrine activity of complex mixtures.

Toxicity of zearalenone and its nutritional intervention by natural products

Zearalenone (ZEN) is a toxic secondary metabolite mainly produced by fungi of the genus Fusarium, and is often present in various food and feed ingredients such as corn and wheat. The structure of ZEN is similar to that of natural estrogen, and it can bind to estrogen receptors and has estrogenic activity. Therefore, it can cause endocrine-disrupting effects and promote the proliferation of estrogen receptor-positive cell lines. In addition, ZEN can cause oxidative damage, endoplasmic reticulum stress, apoptosis, and other hazards, resulting in systemic toxic effects, including reproductive toxicity, hepatotoxicity, and immunotoxicity. In the past few decades, researchers have tried many ways to remove ZEN from food and feed, but it is still a challenge to eliminate it. In recent years, natural compounds have become of interest for their excellent protective effects on human health from food contaminants. Researchers have discovered that natural compounds often used as dietary supplements can effectively alleviate ZEN-induced systemic toxic effects. Most of the compounds mitigate ZEN-induced toxicity through antioxidant effects. In this article, the contamination of food and feed by ZEN and the various toxic effects and mechanisms of ZEN are reviewed, as well as the mitigation effects of natural compounds on ZEN-induced toxicity.

In silico analysis and in planta production of recombinant ccl21/IL1β protein and characterization of its in vitro anti-tumor and immunogenic activity

CCL21 has an essential role in anti-tumor immune activity. Epitopes of IL1β have adjuvant activity without causing inflammatory responses. CCR7 and its ligands play a vital role in the immune balance; specifically, in transport of T lymphocytes and antigen-presenting cells such as dendritic cells to the lymph nodes. This study aimed to produce epitopes of CCL21 and IL1β as a recombinant protein and characterize its in vitro anti-tumor and immunogenic activity. A codon-optimized ccl21/IL1β gene was designed and synthesized from human genes. Stability and binding affinity of CCL21/IL1β protein and CCR7 receptor were examined through in silico analyses. The construct was introduced into N. tabacum to produce this recombinant protein and the structure and function of CCL21/IL1β were examined. Purified protein from transgenic leaves generated a strong signal in SDS PAGE and western blotting assays. FTIR measurement and MALDI-TOF/TOF mass spectrography showed that ccl21/IL-1β was correctly expressed in tobacco plants. Potential activity of purified CCL21/IL1β in stimulating the proliferation and migration of MCF7 cancer cell line was investigated using the wound healing method. The results demonstrated a decrease in survival rate and metastasization of cancer cells in the presence of CCL21/IL1β, and IC50 of CCL21 on MCF7 cells was less than that of non-recombinant protein. Agarose assay on PBMCsCCR7+ showed that CCL21/IL1β has biological activity and there is a distinguishable difference between chemokinetic (CCL21) and chemotactic (FBS) movements. Overall, the results suggest that CCL21/IL1β could be considered an effective adjuvant in future in vivo and clinical tests.

A synergism of in silico and statistical approaches to discover new potential endocrine disruptor mycotoxins

Mycotoxins are secondary metabolites produced by pathogenic fungi. They are found in a variety of different products, such as spices, cocoa, and cereals, and they can contaminate fields before and/or after harvest and during storage. Mycotoxins negatively impact human and animal health, causing a variety of adverse effects, ranging from acute poisoning to long-term effects. Given a large number of mycotoxins (currently more than 300 are known), it is impossible to use in vitro/in vivo methods to detect the potentially harmful effects to human health of all of these. To overcome this problem, this work aims to present a new robust computational approach, based on a combination of in silico and statistical methods, in order to screen a large number of molecules against the nuclear receptor family in a cost and time-effective manner and to discover the potential endocrine disruptor activity of mycotoxins. The results show that a high number of mycotoxins is predicted as a potential binder of nuclear receptors. In particular, ochratoxin A, zearalenone, α- and β-zearalenol, aflatoxin B1, and alternariol have been shown to be putative endocrine disruptors chemicals for nuclear receptors.

Investigation on electronic structure, vibrational spectra, NBO analysis, and molecular docking studies of aflatoxins and selected emerging mycotoxins against wild-type androgen receptor

The geometry, frontier molecular orbitals (FMOs), vibrational, NBO analysis, and molecular docking simulations of aflatoxins (B1, B2, M1, M2, G1, G2), zearalenone (ZEA) emodin (EMO), alternariol (AOH), alternariol monoethyl ether (AMME), and tenuazonic acid (TeA) mycotoxins have been extensively theoretically studied and discussed based on quantum density functional theory calculations using Gaussian 16 software package. The theoretical computation for the geometry optimization, NBOs, and the molecular docking interaction was conducted using Density Functional Theory with B3LYP/6-31+G(d,p), NBO program, and AutoDock Vina tools respectively. Charge delocalization patterns and second-order perturbation energies of the most interacting natural bond orbitals (NBOs) of these mycotoxins have also been computed and predicted. Interestingly, among the mycotoxins investigated, aflatoxin G1 is seen to give the strongest stabilization energy while Zearalenone shows the highest tendency to accept electron(s) and emodin, an emerging mycotoxin gave the best binding pose within the androgen receptor pocket with a mean binding affinity of -7.40 kcal/mol.

Computational Study on Temperature Driven Structure-Function Relationship of Polysaccharide Producing Bacterial Glycosyl Transferase Enzyme

Glycosyltransferase (GTs) is a wide class of enzymes that transfer sugar moiety, playing a key role in the synthesis of bacterial exopolysaccharide (EPS) biopolymer. In recent years, increased demand for bacterial EPSs has been observed in pharmaceutical, food, and other industries. The application of the EPSs largely depends upon their thermal stability, as any industrial application is mainly reliant on slow thermal degradation. Keeping this in context, EPS producing GT enzymes from three different bacterial sources based on growth temperature (mesophile, thermophile, and hyperthermophile) are considered for in silico analysis of the structural–functional relationship. From the present study, it was observed that the structural integrity of GT increases significantly from mesophile to thermophile to hyperthermophile. In contrast, the structural plasticity runs in an opposite direction towards mesophile. This interesting temperature-dependent structural property has directed the GT–UDP-glucose interactions in a way that thermophile has finally demonstrated better binding affinity (−5.57 to −10.70) with an increased number of hydrogen bonds (355) and stabilizing amino acids (Phe, Ala, Glu, Tyr, and Ser). The results from this study may direct utilization of thermophile-origin GT as best for industrial-level bacterial polysaccharide production.

Thinking Out of the Box: On the Ability of Zea mays L. to Biotrasform Aflatoxin B1 Into Its Modified Forms

While aflatoxin metabolism in animals has been clarified, very limited information is so far available on the possible biotransformation occurring in plants. Therefore, this work aimed at investigating whether AFB1 metabolites could occur in field-grown infected maize and the putative role of Zea mays L. metabolism in their production. For such scope, asymptomatic in vitro-grown plantlets and in silico evaluations of plant transforming enzymes were used to pinpoint how plants may handle these compounds. Our data demonstrated the role of maize plants in the production of Phase I hydroxylated aflatoxins, including, among others, AFM1, AFM2, and aflatoxicol, and suggest that plant cytochromes may be involved in this biotransformation of AFB1.

Molecular Aspects of Mycotoxins-A Serious Problem for Human Health

Mycotoxins are toxic fungal secondary metabolities formed by a variety of fungi (moulds) species. Hundreds of potentially toxic mycotoxins have been already identified and are considered a serious problem in agriculture, animal husbandry, and public health. A large number of food-related products and beverages are yearly contaminated by mycotoxins, resulting in economic welfare losses. Mycotoxin indoor environment contamination is a global problem especially in less technologically developed countries. There is an ongoing effort in prevention of mould growth in the field and decontamination of contaminated food and feed in order to protect human and animal health. It should be emphasized that the mycotoxins production by fungi (moulds) species is unavoidable and that they are more toxic than pesticides. Human and animals are exposed to mycotoxin via food, inhalation, or contact which can result in many building-related illnesses including kidney and neurological diseases and cancer. In this review, we described in detail the molecular aspects of main representatives of mycotoxins, which are serious problems for global health, such as aflatoxins, ochratoxin A, T-2 toxin, deoxynivalenol, patulin, and zearalenone.

Co-Occurrence and Levels of Mycotoxins in Fish Feeds in Kenya

This study determined the presence, levels and co-occurrence of mycotoxins in fish feeds in Kenya. Seventy-eight fish feeds and ingredients were sampled from fish farms and fish feed manufacturing plants and analysed for 40 mycotoxins using high-performance liquid chromatography-high resolution mass spectrometry. Twenty-nine (73%) mycotoxins were identified with 76 (97%) samples testing positive for mycotoxins presence. Mycotoxins with the highest prevalences were enniatin B (91%), deoxynivalenol (76%) and fumonisin B1 (54%) while those with the highest maximum levels were sterigmatocystin (<30.5-3517.1 µg/kg); moniliformin (<218.9-2583.4 µg/kg) and ergotamine (<29.3-1895.6 µg/kg). Mycotoxin co-occurrence was observed in 68 (87%) samples. Correlations were observed between the fumonisins; enniatins B and zearalenone and its metabolites. Fish dietary exposure estimates ranged between <0.16 and 43.38 µg/kg body weight per day. This study shows evidence of mycotoxin presence and co-occurrence in fish feeds and feed ingredients in Kenya. Fish exposure to these levels of mycotoxins over a long period of time may lead to adverse health effects due to their possible additive, synergistic or antagonist toxic effects. Measures to reduce fish feed mycotoxin contamination should be taken to avoid mycotoxicosis in fish and subsequently in humans and animals through residues.

Preventing the Interaction between Coronaviruses Spike Protein and Angiotensin I Converting Enzyme 2: An In Silico Mechanistic Case Study on Emodin as a Potential Model Compound

Emodin, a widespread natural anthraquinone, has many biological activities including health-protective and adverse effects. Amongst beneficial effects, potential antiviral activity against coronavirus responsible for the severe acute respiratory syndrome outbreak in 2002–2003 has been described associated with the inhibition of the host cells target receptors recognition by the viral Spike protein. However, the inhibition mechanisms have not been fully characterized, hindering the rational use of emodin as a model compound to develop more effective analogues. This work investigates emodin interaction with the Spike protein to provide a mechanistic explanation of such inhibition. A 3D molecular modeling approach consisting of docking simulations, pharmacophoric analysis and molecular dynamics was used. The plausible mechanism is described as an interaction of emodin at the protein–protein interface which destabilizes the viral protein-target receptor complex. This analysis has been extended to the Spike protein of the coronavirus responsible for the current pandemic hypothesizing emodin’s functional conservation. This solid knowledge-based foothold provides a possible mechanistic rationale of the antiviral activity of emodin as a future basis for the potential development of efficient antiviral cognate compounds. Data gaps and future work on emodin-related adverse effects in parallel to its antiviral pharmacology are explored.

Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens

In this review, applications of various molecular modelling methods in the study of estrogens and xenoestrogens are summarized. Selected biomolecules that are the most commonly chosen as molecular modelling objects in this field are presented. In most of the reviewed works, ligand docking using solely force field methods was performed, employing various molecular targets involved in metabolism and action of estrogens. Other molecular modelling methods such as molecular dynamics and combined quantum mechanics with molecular mechanics have also been successfully used to predict the properties of estrogens and xenoestrogens. Among published works, a great number also focused on the application of different types of quantitative structure-activity relationship (QSAR) analyses to examine estrogen's structures and activities. Although the interactions between estrogens and xenoestrogens with various proteins are the most commonly studied, other aspects such as penetration of estrogens through lipid bilayers or their ability to adsorb on different materials are also explored using theoretical calculations. Apart from molecular mechanics and statistical methods, quantum mechanics calculations are also employed in the studies of estrogens and xenoestrogens. Their applications include computation of spectroscopic properties, both vibrational and Nuclear Magnetic Resonance (NMR), and also in quantum molecular dynamics simulations and crystal structure prediction. The main aim of this review is to present the great potential and versatility of various molecular modelling methods in the studies on estrogens and xenoestrogens.