Structure-activity relationship of pharmacophores and toxicophores: the need for clinical strategy.

OBJECTIVES: Sometimes clinical efficacy and potential risk of therapeutic and toxic agents are difficult to predict over a long period of time. Hence there is need for literature search with a view to assessing cause of toxicity and less efficacy of drugs used in clinical practice. METHOD: Hence literatures were searched for physicochemical properties, chemical formulas, molecular masses, pH values, ionization, receptor type, agonist and antagonist, therapeutic, toxic and structure-activity relationship of chemical compounds with pharmacophore and toxicophore, with a view to identifying high efficacious and relative low toxic agents. Inclusion criteria were manuscripts published on PubMed, Scopus, Web of Science, PubMed Central, Google Scholar among others, between 1960 and 2023. Keywords such as pharmacophore, toxicophore, structure-activity-relationship and disease where also searched. The exclusion criteria were the chemicals that lack pharmacophore, toxicophore and manuscripts published before 1960. RESULTS: Findings have shown that pharmacophore and toxicophore functional groups determine clinical efficacy and safety of therapeutics, but if they overlap therapeutic and toxicity effects go concurrently. Hence the functional groups, dose, co-administration and concentration of drugs at receptor, drug-receptor binding and duration of receptor binding are the determining factors of pharmacophore and toxicophore activity. Molecular mass, chemical configuration, pH value, receptor affinity and binding capacity, multiple pharmacophores, hydrophilic/lipophilic nature of the chemical contribute greatly to functionality of pharmacophore and toxicophore. CONCLUSION: Daily single therapy, avoidance of reversible pharmacology, drugs with covalent adduct, maintenance of therapeutic dose, and the use of multiple pharmacophores for terminal diseases will minimize toxicity and improve efficacy.

Automated detection of toxicophores and prediction of mutagenicity using PMCSFG algorithm.

Maximum common substructures (MCS) have received a lot of attention in the chemoinformatics community. They are typically used as a similarity measure between molecules, showing high predictive performance when used in classification tasks, while being easily explainable substructures. In the present work, we applied the Pairwise Maximum Common Subgraph Feature Generation (PMCSFG) algorithm to automatically detect toxicophores (structural alerts) and to compute fingerprints based on MCS. We present a comparison between our MCS-based fingerprints and 12 well-known chemical fingerprints when used as features in machine learning models. We provide an experimental evaluation and discuss the usefulness of the different methods on mutagenicity data. The features generated by the MCS method have a state-of-the-art performance when predicting mutagenicity, while they are more interpretable than the traditional chemical fingerprints.

Occurrence and Ecological Impact of Chemical Mixtures in a Semiclosed Sea by Suspect Screening Analysis.

Stress from mixtures of synthetic chemicals is among the key issues that have significant adverse impacts on the marine ecosystems. A robust screening workflow integrating toxicological-based ranking schemes is still deficient for comprehensive investigation on the main constituents in chemical mixtures that contribute to the ecological risks. In this study, the presence and compositions of a collection of priority pollutants were monitored by suspect screening analysis of seawater and estuarine water samples from the semiclosed Bohai Sea. In total, 108 organic pollutants in nine use categories were identified. Pesticides, intermediates, plastic additives, and per- and polyfluoroalkyl substances were the extensively detected chemical groups. Varied distribution patterns of the pollutants were illustrated intuitively in distinctive sampling areas by hierarchical cluster analysis, which were mainly influenced by run-off inputs, ocean currents, and chemical use history. Ecological risks of chemicals with quantified residue levels were first assessed by the toxicity-weighted concentration ranking scheme, and pentachlorophenol was found as the main contributor in the investigating areas. By optimization of multiple alternative variables (e.g., instrumental response and detection frequency), extended ranking of all the identified pollutants was plausible under the toxicological priority index framework. Similarity in toxicological endpoints of the prioritized pollutants could further been screened by ToxAlerts. Aromatic amine was highlighted as the most frequently detected structural alert (SA) for genotoxic carcinogenicity and mutagenicity. These findings fully demonstrate rationality of the ranking schemes integrated into the suspect screening analysis for profiling contamination characteristics, assessing ecological risk potentials, and prioritizing SAs.

Use of Lhasa Limited Products for the In Silico Prediction of Drug Toxicity.

Lhasa Limited have had a role in the in silico prediction of drug and other chemical toxicity for over 30 years. This role has always been multifaceted, both as a provider of predictive software such as Derek Nexus, and as an honest broker for the sharing of proprietary chemical and toxicity data. A changing regulatory environment and the drive for the Replacement, Reduction and Refinement (the 3Rs) of animal testing have led both to increased acceptance of in silico predictions and a desire for the sharing of data to reduce duplicate testing. The combination of these factors has led to Lhasa Limited providing a suite of products and coordinating numerous data-sharing consortia that do indeed facilitate a significant reduction in the testing burden that companies would otherwise be laboring under. Many of these products and consortia can be organized into workflows for specific regulatory use cases, and it is these that will be used to frame the narrative in this chapter.

The mitochondrial paradox.

A structural motif that is found in two cancer drugs may be responsible for their ability to tackle cancers and for the side-effects caused by the drugs.

Identification of a novel toxicophore in anti-cancer chemotherapeutics that targets mitochondrial respiratory complex I.

Disruption of mitochondrial function selectively targets tumour cells that are dependent on oxidative phosphorylation. However, due to their high energy demands, cardiac cells are disproportionately targeted by mitochondrial toxins resulting in a loss of cardiac function. An analysis of the effects of mubritinib on cardiac cells showed that this drug did not inhibit HER2 as reported, but directly inhibits mitochondrial respiratory complex I, reducing cardiac-cell beat rate, with prolonged exposure resulting in cell death. We used a library of chemical variants of mubritinib and showed that modifying the 1H-1,2,3-triazole altered complex I inhibition, identifying the heterocyclic 1,3-nitrogen motif as the toxicophore. The same toxicophore is present in a second anti-cancer therapeutic carboxyamidotriazole (CAI) and we demonstrate that CAI also functions through complex I inhibition, mediated by the toxicophore. Complex I inhibition is directly linked to anti-cancer cell activity, with toxicophore modification ablating the desired effects of these compounds on cancer cell proliferation and apoptosis.

The pharmaceutical industry needs to make safe and effective drugs. At the same time this industry is under pressure to keep the costs of developing these drugs at an acceptable level. Drugs work by interacting with and typically blocking a specific target, such as a protein in a particular type of cell. Sometimes, however, drugs also bind other unexpected targets. These "off-target" effects can be the reason for a drug's toxicity, and it is important ­ both for the benefit of patients and the money that can be saved when developing drugs ­ to identify how drugs cause toxic side effects. The earlier researchers detect off-target effects, the better. Recent data has suggested that an anti-cancer drug called mubritinib has off-target effects on the compartments within cells that provide the cell with most of their energy, the mitochondria. This drug's intended target is a protein called HER2, which is found in large amounts on the surfaces of some breast cancer cells. Yet if mubritinib has this off-target effect on mitochondria, it may be harmful to other cells including heart cells because the heart is an organ that needs a large amount of energy from its mitochondria. Stephenson et al. have now performed experiments to show that mubritinib does not actually interact with HER2 at all, but only targets mitochondria. The effect of mubritinib as an anti-cancer drug is therefore only due to its activity against mitochondria. Digging deeper into the chemistry revealed the small parts of its chemical structure that was responsible for mubritinib's toxicity against heart cells, the so-called toxic substructure. Another anti-cancer drug called carboxyamidotriazole also has the same toxic substructure. Carboxyamidotriazole is supposed to stop cells from taking up calcium ions, but a final set of experiments demonstrated that this drug also only acts by inhibiting mitochondria. Often there is not enough information about many drugs' substructures, meaning off-target effects and toxicities cannot be predicted. The pharmaceutical industry will now be able to benefit from this new knowledge about the toxic substructures within some drugs. This research may also help patients who take mubritinib or carboxyamidotriazole, because their doctors will have to check for side effects on the heart more carefully.

Elucidating the aryl hydrocarbon receptor antagonism from a chemical-structural perspective.

The aryl hydrocarbon receptor (AhR) plays an important role in several biological processes such as reproduction, immunity and homoeostasis. However, little is known on the chemical-structural and physicochemical features that influence the activity of AhR antagonistic modulators. In the present report, in vitro AhR antagonistic activity evaluations, based on a chemical-activated luciferase gene expression (AhR-CALUX) bioassay, and an extensive literature review were performed with the aim of constructing a structurally diverse database of contaminants and potentially toxic chemicals. Subsequently, QSAR models based on Linear Discriminant Analysis and Logistic Regression, as well as two toxicophoric hypotheses were proposed to model the AhR antagonistic activity of the built dataset. The QSAR models were rigorously validated yielding satisfactory performance for all classification parameters. Likewise, the toxicophoric hypotheses were validated using a diverse set of 350 decoys, demonstrating adequate robustness and predictive power. Chemical interpretations of both the QSAR and toxicophoric models suggested that hydrophobic constraints, the presence of aromatic rings and electron-acceptor moieties are critical for the AhR antagonism. Therefore, it is hoped that the deductions obtained in the present study will contribute to elucidate further on the structural and physicochemical factors influencing the AhR antagonistic activity of chemical compounds.

Fluoromethylcyclopropylamine derivatives as potential in vivo toxicophores - A cautionary disclosure.

Fluorination of metabolic hotspots in a molecule is a common medicinal chemistry strategy to improve in vivo half-life and exposure and, generally, this strategy offers significant benefits. Here, we report the application of this strategy to a series of poly-ADP ribose glycohydrolase (PARG) inhibitors, resulting in unexpected in vivo toxicity which was attributed to this single-atom modification.

The use of high-throughput screening techniques to evaluate mitochondrial toxicity.

Toxicologists and chemical regulators depend on accurate and effective methods to evaluate and predict the toxicity of thousands of current and future compounds. Robust high-throughput screening (HTS) experiments have the potential to efficiently test large numbers of chemical compounds for effects on biological pathways. HTS assays can be utilized to examine chemical toxicity across multiple mechanisms of action, experimental models, concentrations, and lengths of exposure. Many agricultural, industrial, and pharmaceutical chemicals classified as harmful to human and environmental health exert their effects through the mechanism of mitochondrial toxicity. Mitochondrial toxicants are compounds that cause a decrease in the number of mitochondria within a cell, and/or decrease the ability of mitochondria to perform normal functions including producing adenosine triphosphate (ATP) and maintaining cellular homeostasis. Mitochondrial dysfunction can lead to apoptosis, necrosis, altered metabolism, muscle weakness, neurodegeneration, decreased organ function, and eventually disease or death of the whole organism. The development of HTS techniques to identify mitochondrial toxicants will provide extensive databases with essential connections between mechanistic mitochondrial toxicity and chemical structure. Computational and bioinformatics approaches can be used to evaluate compound databases for specific chemical structures associated with toxicity, with the goal of developing quantitative structure-activity relationship (QSAR) models and mitochondrial toxicophores. Ultimately these predictive models will facilitate the identification of mitochondrial liabilities in consumer products, industrial compounds, pharmaceuticals and environmental hazards.

Why are most phospholipidosis inducers also hERG blockers?

Recent reports have noted that a number of compounds that block the human Ether-à-go-go related gene (hERG) ion channel also induce phospholipidosis (PLD). To explore a hypothesis explaining why most PLD inducers are also hERG inhibitors, a modeling approach was undertaken with data sets comprised of 4096 compounds assayed for hERG inhibition and 5490 compounds assayed for PLD induction. To eliminate the chemical domain effect, a filtered data set of 567 compounds tested in quantitative high-throughput screening (qHTS) format for both hERG inhibition and PLD induction was constructed. Partial least squares (PLS) modeling followed by 3D-SDAR mapping of the most frequently occurring bins and projection on to the chemical structure suggested that both adverse effects are driven by similar structural features, namely two aromatic rings and an amino group forming a three-center toxicophore. Non-parametric U-tests performed on the original 3D-SDAR bins indicated that the distance between the two aromatic rings is the main factor determining the differences in activity; at distances of up to about 5.5 Å, a phospholipidotic compound would also inhibit hERG, while at longer distances, a sharp reduction of the PLD-inducing potential leaves only a well-pronounced hERG blocking effect. The hERG activity itself diminishes after the distance between the centroids of the two aromatic rings exceeds 12.5 Å. Further comparison of the two toxicophores revealed that the almost identical aromatic rings to amino group distances play no significant role in distinguishing between PLD and hERG activity. The hypothesis that the PLD toxicophore appears to be a subset of the hERG toxicophore explains why about 80% of all phospholipidotic chemicals (the remaining 20% are thought to act via a different mechanism) also inhibit the hERG ion channel. These models were further validated in large-scale qHTS assays testing 1085 chemicals for their PLD-inducing potential and 1570 compounds for hERG inhibition. After removal of the modeling and experimental inconclusive compounds, the area under the receiver-operating characteristic (ROC) curve was 0.92 for the PLD model and 0.87 for the hERG model. Due to the exceptional ability of these models to recognize safe compounds (negative predictive values of 0.99 for PLD and 0.94 for hERG were achieved), their use in regulatory settings might be particularly useful.

3D-SDAR modeling of hERG potassium channel affinity: A case study in model design and toxicophore identification.

A dataset of 237 human Ether-à-go-go Related Gene (hERG) potassium channel inhibitors (180 of which were used for model building and validation, whereas 57 constituted the "true" external prediction set) collected from 22 literature sources was modeled by 3D-SDAR. To produce reliable and reproducible classification models for hERG blocking, the initial set of 180 chemicals was split into two subsets: a balanced modeling set consisting of 118 compounds and an unbalanced validation set comprised of 62 compounds. A PLS bagging-like algorithm written in Matlab was used to process the data and assign each compound to one of the two (hERG+ or hERG-) activity classes. The best predictive model evaluated on the basis of a fully randomized hold-out test set (comprising 20% of the modeling set) used 4 latent variables and a grid of 6ppm×6ppm×1Å in the C-C region, 6ppm×30ppm×1Å in the C-N region, and 30ppm×30ppm×1Å in the N-N region. An overall accuracy of 0.84 was obtained for both the hold-out test set and the validation set. Further, an external prediction set consisting of 57 drugs and drug derivatives was used to estimate the true predictive power of the reported 3D-SDAR model - a slight reduction of the overall accuracy down to 0.77 was observed. 3D-SDAR map of the most frequently occurring bins and their projection on the standard coordinate space of the chemical structures allowed identification of a three-center toxicophore composed of two aromatic rings and an amino group. A U test along the distance axis of the most frequently occurring 3D-SDAR bins was used to set the distance limits of the toxicophore. This toxicophore was found to be similar to an earlier reported phospholipidosis (PLD) toxicophore.

Toxicophore exploration as a screening technology for drug design and discovery: techniques, scope and limitations.

Toxicity is a common drawback of newly designed chemotherapeutic agents. With the exception of pharmacophore-induced toxicity (lack of selectivity at higher concentrations of a drug), the toxicity due to chemotherapeutic agents is based on the toxicophore moiety present in the drug. To date, methodologies implemented to determine toxicophores may be broadly classified into biological, bioanalytical and computational approaches. The biological approach involves analysis of bioactivated metabolites, whereas the computational approach involves a QSAR-based method, mapping techniques, an inverse docking technique and a few toxicophore identification/estimation tools. Being one of the major steps in drug discovery process, toxicophore identification has proven to be an essential screening step in drug design and development. The paper is first of its kind, attempting to cover and compare different methodologies employed in predicting and determining toxicophores with an emphasis on their scope and limitations. Such information may prove vital in the appropriate selection of methodology and can be used as screening technology by researchers to discover the toxicophoric potentials of their designed and synthesized moieties. Additionally, it can be utilized in the manipulation of molecules containing toxicophores in such a manner that their toxicities might be eliminated or removed.

Evaluating the risk of phospholipidosis using a new multidisciplinary pipeline approach.

Phospholipidosis (PLD) is an undesirable potential side-effect of drugs, and cationic amphiphilic drugs (CADs) represent the main class of PLD inducers. A CADs toxicophore has been recently proposed, although the CADs definition is far from being trivial. In this work we derive a three-dimensional CADs toxicophore (here named PLD-phore) using a molecular interaction field approach, and test its suitability to discriminate between PLD inducers and non-inducers in a virtual screening approach. Ten commercially available compounds predicted to be PLD inducers and non-inducers based on their similarity to the PLD-phore were experimentally tested for PLD induction using two cell-based in vitro assays (fluorescent lipid uptake, activity of secreted lysosomal ß-hexosaminidase). When a positive effect was observed, the PLD induction was also confirmed by transmission electron microscopy. Two exceptions to the general statement about CADs and PLD induction were detected and discussed, and for one compound the cell-based in-vitro assays lead to different outcomes.

Computational identification of a phospholipidosis toxicophore using (13)C and (15)N NMR-distance based fingerprints.

Modified 3D-SDAR fingerprints combining (13)C and (15)N NMR chemical shifts augmented with inter-atomic distances were used to model the potential of chemicals to induce phospholipidosis (PLD). A curated dataset of 328 compounds (some of which were cationic amphiphilic drugs) was used to generate 3D-QSDAR models based on tessellations of the 3D-SDAR space with grids of different density. Composite PLS models averaging the aggregated predictions from 100 fully randomized individual models were generated. On each of the 100 runs, the activities of an external blind test set comprised of 294 proprietary chemicals were predicted and averaged to provide composite estimates of their PLD-inducing potentials (PLD+ if PLD is observed, otherwise PLD-). The best performing 3D-QSDAR model utilized a grid with a density of 8ppm×8ppm in the C-C region, 8ppm×20ppm in the C-N region and 20ppm×20ppm in the N-N region. The classification predictive performance parameters of this model evaluated on the basis of the external test set were as follows: accuracy=0.70, sensitivity=0.73 and specificity=0.66. A projection of the most frequently occurring bins on the standard coordinate space suggested a toxicophore composed of an aromatic ring with a centroid 3.5-7.5Å distant from an amino-group. The presence of a second aromatic ring separated by a 4-5Å spacer from the first ring and at a distance of between 5.5Å and 7Å from the amino-group was also associated with a PLD+ effect. These models provide comparable predictive performance to previously reported models for PLD with the added benefit of being based entirely on non-confidential, publicly available training data and with good predictive performance when tested in a rigorous, external validation exercise.

Exploring QSTR modeling and toxicophore mapping for identification of important molecular features contributing to the chemical toxicity in Escherichia coli.

Biodiversity deprivation can affect functions and services of the ecosystem. Changes in biodiversity alter ecosystem processes and change the resilience of ecosystems to ecological changes. Bacterial communities are the main form of biomass in the ecosystem and one of largest populations on the planet. Bacterial communities provide important services to biodiversity. They break down pollutants, municipal waste and ingested food, and they are the primary route for recycling of organic matter to plants and other autotrophs, conversion of inorganic matter into new biological tissue using sunlight, management of energy crisis through use of biofuel. In the present study, computational chemistry and statistical modeling have been used to develop mathematical equations which can be applied to calculate toxicity of new/unknown chemicals/biofuels/metabolites in Escherichia coli. 2D and 3D descriptors were generated from molecular structure of compounds and mathematical models have been developed using genetic function approximation followed by multiple linear regression (GFA-MLR) method. Model validity was checked through defined internal (R(2)=0.751 and Q(2)=0.711), and external (Rpred(2)=0.773) statistical parameters. Molecular features responsible for toxicity were also assessed through 3D toxicophore study. The toxicophore-based model was validated (R=0.785) using qualitative statistical metrics and randomization test (Fischer validation).

High-throughput respirometric assay identifies predictive toxicophore of mitochondrial injury.

Many environmental chemicals and drugs negatively affect human health through deleterious effects on mitochondrial function. Currently there is no chemical library of mitochondrial toxicants, and no reliable methods for predicting mitochondrial toxicity. We hypothesized that discrete toxicophores defined by distinct chemical entities can identify previously unidentified mitochondrial toxicants. We used a respirometric assay to screen 1760 compounds (5 µM) from the LOPAC and ChemBridge DIVERSet libraries. Thirty-one of the assayed compounds decreased uncoupled respiration, a stress test for mitochondrial dysfunction, prior to a decrease in cell viability and reduced the oxygen consumption rate in isolated mitochondria. The mitochondrial toxicants were grouped by chemical similarity and two clusters containing four compounds each were identified. Cheminformatic analysis of one of the clusters identified previously uncharacterized mitochondrial toxicants from the ChemBridge DIVERSet. This approach will enable the identification of mitochondrial toxicants and advance the prediction of mitochondrial toxicity for both drug discovery and risk assessment.