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.

The VERA software: Implementation of the acute fish toxicity endpoint and its application to pharmaceutical compounds.

The Virtual Extensive Read-Across software (VERA) is a new tool for read-across using a global similarity score, molecular groups, and structural alerts to find clusters of similar substances; these clusters are then used to identify suitable similar substances and make an assessment for the target substance. A beta version of VERA GUI is free and available at vegahub.eu; the source code of the VERA algorithm is available on GitHub. In the past we described its use to assess carcinogenicity, a classification endpoint. The aim here is to extend the automated read-across approach to assess continuous endpoints as well. We addressed acute fish toxicity. VERA evaluation on the acute fish toxicity endpoint was done on a dataset containing general substances (pesticides, industrial products, biocides, etc.), obtaining an overall R2 of 0.68. We employed the VERA algorithm also on active pharmaceutical ingredients (APIs). We included a portion of the APIs in the training dataset to predict APIs, successfully achieving an overall R2 of 0.63. VERA evaluates the assessment's reliability, and we reached an R2 of 0.78 and Root Mean Square Error (RMSE) of 0.44 for predictions with high reliability.

Prioritisation of data-poor pharmaceuticals for empirical testing and environmental risk assessment.

There are more than 3,500 active pharmaceutical ingredients (APIs) on the global market for human and veterinary use. Residues of these APIs eventually reach the aquatic environment. Although an environmental risk assessment (ERA) for marketing authorization applications of medicinal products is mandatory in the European Union since 2006, an ERA is lacking for most medicines approved prior to 2006 (legacy APIs). Since it is unfeasible to perform extensive ERA tests for all these legacy APIs, there is a need for prioritization of testing based on the limited data available. Prioritized APIs can then be further investigated to estimate their environmental risk in more detail. In this study, we prioritized more than 1,000 APIs used in Europe based on their predicted risk for aquatic freshwater ecosystems. We determined their risk by combining an exposure estimate (Measured or Predicted Environmental Concentration; MEC or PEC, respectively) with a Predicted No Effect Concentration (PNEC). We developed several procedures to combine the limited empirical data available with in silico data, resulting in multiple API rankings varying in data needs and level of conservativeness. In comparing empirical with in silico data, our analysis confirmed that the PEC estimated with the default parameters used by the European Medicines Agency often - but not always - represents a worst-case scenario. Comparing the ecotoxicological data for the three main taxonomic groups, we found that fish represents the most sensitive species group for most of the APIs in our list. We furthermore show that the use of in silico tools can result in a substantial underestimation of the ecotoxicity of APIs. After combining the different exposure and effect estimates into four risk rankings, the top-ranking APIs were further screened for availability of ecotoxicity data in data repositories. This ultimately resulted in the prioritization of 15 APIs for further ecotoxicological testing and/or exposure assessment.

The VEGA Tool to Check the Applicability Domain Gives Greater Confidence in the Prediction of In Silico Models.

A sound assessment of in silico models and their applicability domain can support the use of new approach methodologies (NAMs) in chemical risk assessment and requires increasing the users' confidence in this approach. Several approaches have been proposed to evaluate the applicability domain of such models, but their prediction power still needs a thorough assessment. In this context, the VEGA tool capable of assessing the applicability domain of in silico models is examined for a range of toxicological endpoints. The VEGA tool evaluates chemical structures and other features related to the predicted endpoints and is efficient in measuring applicability domain, enabling the user to identify less accurate predictions. This is demonstrated with many models addressing different endpoints, towards toxicity of relevance to human health, ecotoxicological endpoints, environmental fate, physicochemical and toxicokinetic properties, for both regression models and classifiers.

Intelligent consensus predictions of bioconcentration factor of pharmaceuticals using 2D and fragment-based descriptors.

Bioconcentration factors (BCFs) are markers of chemical substance accumulation in organisms, and they play a significant role in determining the environmental risk of various chemicals. Experiments to obtain BCFs are expensive and time-consuming; therefore, it is better to estimate BCF early in the chemical development process. The current research aims to evaluate the ecotoxicity potential of 122 pharmaceuticals and identify possible important structural attributes using BCF as the determining feature against a group of fish species. We have calculated the theoretical 2D descriptors from the OCHEM platform and SiRMS descriptor calculating software. The regression-based quantitative structure-property relationship (QSPR) modeling was used to identify the chemical features responsible for acute fish bioconcentration. Multiple models with the "intelligent consensus" algorithm were employed for the regression-based approach improving the predictive ability of the models. To ensure the robustness and interpretability of the developed models, rigorous validation was performed employing various statistical internal and external validation metrics. From the developed models, it can be specified that the presence of large lipophilic and electronegative moieties greatly enhances the bioaccumulative potential of pharmaceuticals, whereas the hydrophilic characteristics have shown a negative impact on BCF. Furthermore, the developed models were employed to screen the DrugBank database (https://go.drugbank.com/) for assessing the BCF properties of the entire database. The evidence acquired from the modeled descriptors might be used for aquatic risk assessment in the future, with the added benefit of providing an early caution of their probable negative impact on aquatic ecosystems for regulatory purposes.

Development of new QSAR models for water, sediment, and soil half-life.

Checking the persistence of a chemical in the environment is extremely important. Regulations like REACH, the European one on chemicals, require the measurements or estimates of the half-life of the chemical in water, sediment, and soil. The use of non-testing methods, like quantitative structure-activity relationship (QSAR) models, is encouraged because it reduces costs and time. To our knowledge, there are very few freely available models for these properties and some are for specific chemical classes. Here, we present three new semi-quantitative models, one for each of the required environmental compartments (water, sediment, and soil). Using literature and REACH registration data, we developed three new counter-propagation artificial neural network models using the CPANNatNIC tool. We calculated the VEGA descriptors, and selected the relevant ones using an internal method in R based on the forward selection technique. The best model for each compartment was implemented in two open-source stand-alone tools, the VEGA platform, and the JANUS tool (https://www.vegahub.eu/). These models were also used by ECHA to build their PBT profiler available in the OECD QSAR toolbox (https://qsartoolbox.org/). Screening and prioritization are also our main target. The models perform well, with R2 always above 0.8 in training and validation. The only exception is the validation set of the soil compartment, with R2 0.68, that is above 0.8 only for compounds inside the applicability domain (automatically calculated by the system). The root mean square error (RMSE) is good, 0.34 or less in log units (again, for soil validation it is higher but it reaches 0.21 when considering only the compounds in the applicability domain). Compared with one of the most widely used tools, BIOWIN3, the proposed models give better results in terms of R2 and RMSE. For the classification, the performance is better for water and soil, and comparable or lower for sediment.

The VEGAHUB Platform: The Philosophy and the Tools.

VEGAHUB (www.vegahub.eu) is a repository of freely available, downloadable tools based on computational toxicology methodologies. The main software tool available in VEGAHUB is VEGA QSAR software encoding more than 90 quantitative structure-activity relationship (QSAR) models for tens of endpoints for human toxicology, ecotoxicology, environmental, physico-chemical and toxicokinetic properties. However, beyond VEGA QSAR, VEGAHUB offers several other tools. Here, we present these resources, the possibilities to fully exploit them and the ways in which to integrate results provided by different VEGAHUB tools. Read-across and weight-of-evidence represent a major advantage of VEGAHUB. Integration between hazard and exposure is provided within innovative tools, which are specific for well-defined scenarios, such as those for cosmetic products. Prioritisation can be achieved by integrating results from 48 models. Finally, we highlight how some tools may not only fit predefined endpoints but also could be applied to general problems and research applications in the QSAR field. A couple of examples are provided, in which a critical assessment of the predictions and the documentation associated with the prediction are considered, in order to properly assess the quality of the results. These results may be associated with different levels of uncertainty or even be conflicting.

Paradox of 'ideal correlations': improved model for air half-life of persistent organic pollutants.

The persistence of organic pollutants is an important environmental property due to the extended possibility to have an impact of corresponding substances. In many cases, the experimental values of the thousands of contaminants are missing. The object of the study is novel computational modelling for air pollutions. Quantitative structure-property relationship (QSPR) for air half-life has been built using the Monte Carlo method with applying the index of ideality of correlation (IIC). The basis of the predictive model of air half-life is the representation of the molecular structure by simplifying molecular input-line entry system (SMILES) and numerical data on the above endpoint (expressed by hours) converted to a decimal logarithm. The statistical quality of the model has been checked up with different validation metrics and is quite good. Paradoxically, the improvement of the statistical quality via the IIC for the validation set is done in detriment to the training set. The new model has performed better than those obtained previously on the same set of compounds, for the prediction of new compounds in the validation set. Some semi-quantitative indicators for the mechanistic interpretation of the model are suggested.

New Models to Predict the Acute and Chronic Toxicities of Representative Species of the Main Trophic Levels of Aquatic Environments.

To assess the impact of chemicals on an aquatic environment, toxicological data for three trophic levels are needed to address the chronic and acute toxicities. The use of non-testing methods, such as predictive computational models, was proposed to avoid or reduce the need for animal models and speed up the process when there are many substances to be tested. We developed predictive models for Raphidocelis subcapitata, Daphnia magna, and fish for acute and chronic toxicities. The random forest machine learning approach gave the best results. The models gave good statistical quality for all endpoints. These models are freely available for use as individual models in the VEGA platform and for prioritization in JANUS software.

First report on chemometric modeling of hydrolysis half-lives of organic chemicals.

Hydrolysis is one of the most important processes of transformation of organic chemicals in water. The rates of reactions, final chemical entities of these processes, and half-lives of organic chemicals are of considerable interest to environmental chemists as well as authorities involved in the controlling the processing and disposal of such organic chemicals. In this study, we have proposed QSPR models for the prediction of hydrolysis half-life of organic chemicals as a function of different pH and temperature conditions using only two-dimensional molecular descriptors with definite physicochemical significance. For each model, suitable subsets of variables were elected using a genetic algorithm method; next, the elected subsets of variables were subjected to the best subset selection with a key objective to determine the best combination of descriptors for model generation. Finally, QSPR models were constructed using the best combination of variables employing the partial least squares (PLS) regression technique. Next, every final model was subjected for strict validation employing the internationally accepted internal and external validation parameters. The proposed models could be applicable for data gap filling to determine hydrolysis half-lives of organic chemicals at different environmental conditions. Generally, presence of aliphatic ether and ether functional groups, high percentage of oxygen content in the molecule and presence of O-Si pairs of atoms at topological distance one, results in a shorter hydrolysis half-life of organic chemicals. On the other hand, higher unsaturation content and high percentage of nitrogen content in molecules lead to higher hydrolysis half-life. It is also found that branched and compact molecules will have a lower half-life while straight chain analogues will have a higher half-life. To the best of our knowledge, the presented models are the first reported QSPR models for hydrolysis half-lives of organic chemicals at different pH values.

Chemometric modeling to predict air half-life of persistent organic pollutants (POPs).

We have reported here a quantitative structure-property relationship (QSPR) model for prediction of air half-life of organic chemicals using a dataset of 302 diverse organic chemicals employing only two-dimensional descriptors with definite physicochemical meaning in order to avoid the computational complexity for higher dimensional molecular descriptors. The developed model was rigorously validated using the internationally accepted internal and external validation metrics. The final partial least squares (PLS) regression model obtained at three latent variables comprises six simple and interpretable 2D descriptors. The simple and highly robust model with good quality of predictions explains 66% for the variance of the training set (R2) (64% in terms of LOO variance (Q2)) and 76% for test set variance (R2pred) (prediction quality). This model might be applicable for data gap filling for determination of POPs in the environment, in case of new or untested chemicals falling within the applicability domain of the model. In general, the model indicates that the air half-life of organic chemicals increases with presence of H-bond acceptor atoms, number of halogen atoms and presence of the R-CH-X fragment and lipophilicity, and decreases with presence of a number of halogens on ring C(sp3) (substitution of halogen atoms on a ring).

New in silico models to predict in vitro micronucleus induction as marker of genotoxicity.

The evaluation of genotoxicity is a fundamental part of the safety assessment of chemicals due to the relevance of the potential health effects of genotoxicants. Among the testing methods available, the in vitro micronucleus assay with mammalian cells is one of the most used and required by regulations targeting several industrial sectors such as the cosmetic industry and food-related sectors. As an alternative to the testing methods, in recent years, lots in silico methods were developed to predict the genotoxicity of chemicals, including models for the Ames mutagenicity test, the in vitro chromosomal aberrations and the in vivo micronucleus assay. We developed several in silico models for the prediction of genotoxicity as reflected by the in vitro micronucleus assay. The resulting models include both statistical and knowledge-based models. The most promising model is the one based on fragments extracted with the SARpy platform. More than 100 structural alerts were extracted, including also fragments associated with the non-genotoxic activity. The model is characterized by high accuracy and the lowest false negative rate, making this tool suitable for chemical screening according to the regulators' needs. The SARpy model will be implemented on the VEGA platform (https://www.vegahub.eu) and will be freely available.

Exploring QSAR modeling of toxicity of chemicals on earthworm.

Earthworm provides sustainability towards the agroecosystem which can be degraded day by day by the extensive use of pesticides (e.g., fungicides, insecticides and herbicides). The present study attempts to develop a predictive quantitative structure-activity relationship (QSAR) model for toxicity of pesticides to earthworm in order to give a suitable guidance for designing new analogues with lower toxicity by exploring the important chemical features which are required to develop safer alternatives. The QSAR model was developed by using the negative logarithm of lethal concentration (pLC50) values of pesticides towards earthworm. We have used various 2D descriptors along with extended topochemical atom (ETA) indices as independent variables for the development of the model. The developed partial least squares (PLS) model was subjected to statistical validation tests proving that the model is statistically reliable and robust (R2 = 0.765, Q2 = 0.614, Q2F1 = 0.734, Q2F2 = 0.713). The contributing descriptors in the model suggested that the pesticides may affect the earthworm nucleic acid through various physicochemical interactions including hydrophobicity, hydrogen bonding, electron donor acceptor complex formation, π-π stacking interaction and charge transfer complex formation.

First report on a classification-based QSAR model for chemical toxicity to earthworm.

As the use of the pesticides has increased extensively in the farming fields to have a better agricultural production, the negative impacts of such use have also increased exponentially. Hence, the toxic effects of pesticides along with the targeted organisms affect the non-targeted terrestrial organisms such as earthworm. Therefore, in the present work, we have developed a classification-based quantitative structure-activity relationship (QSAR) model using linear discriminant analysis (LDA) to capture the specific information of pesticides / diverse chemicals in order to determine the structural information responsible for toxicity manifestation towards the non-targeted organism, i.e., earthworm (Eisenia foetida). After variable selection, the model was developed using 2D descriptors only and was subjected to rigorous statistical validation. The best discriminant model obtained with 8 descriptors showed appreciable Wilks' λ value of 0.490, F (Fischer's statistics) value of 14.03, χ2 value of 79.098, canonical regression coefficient (R) value of 0.714 and ρ value of 14.63. The sensitivity, specificity, accuracy, precision and F-measure values of the training set are 90.00, 80.52, 83.76, 70.59 and 79.12 respectively whereas for the test set, these are 58.82, 79.31, 71.74, 62.50 and 60.61 respectively. The insights obtained from the LDA model suggested that lipophilicity, electronrichness, and lower degree of branching of the organic compounds are responsible for earthworm toxicity through various mechanisms. On the other hand, polar and bulky diverse chemicals do not have such toxic effects on earthworm. Hence, this model can be an effective tool to tailor molecular structures of the existing pesticides to develop novel compounds or pesticides which would be less toxic to the non-targeted organisms, specifically earthworm.

QSAR modeling of Daphnia magna and fish toxicities of biocides using 2D descriptors.

In the recent years, ecotoxicological hazard potential of biocidal products has been receiving increasing attention in the industries and regulatory agencies. Biocides/pesticides are currently one of the most studied groups of compounds, and their registration cannot be done without the empirical toxicity information. In view of limited experimental data available for these compounds, we have developed Quantitative Structure-Activity Relationship (QSAR) models for the toxicity of biocides to fish and Daphnia magna following principles of QSAR modeling recommended by the OECD (Organization for Economic Cooperation and Development). The models were developed using simple and interpretable 2D descriptors and validated using stringent tests. Both models showed encouraging statistical quality in terms of determination coefficient R2 (0.800 and 0.648), cross-validated leave-one-out Q2 (0.760 and 0.602) and predictive R2pred or Q2ext (0.875 and 0.817) for fish (nTraining = 66, nTest = 22) and Daphnia magna (nTraining = 100, nTest = 33) toxicity datasets, respectively. These models should be applicable for data gap filling in case of new or untested biocidal compounds falling within the applicability domain of the models. In general, the models indicate that the toxicity increases with lipophilicity and decreases with polarity, branching and unsaturation. We have also developed interspecies toxicity models for biocides using the daphnia and fish toxicity data and used the models for data gap filling.

Ecotoxicological QSAR modeling of organic compounds against fish: Application of fragment based descriptors in feature analysis.

Organic compounds (OCs) constitute an enormously large class of highly persistent and toxic chemicals widely used for various purposes throughout the world. Their increased detection in water bodies, mainly sewage treatment plants via industrial discharge, has rendered them to become a cause for ecological concern. The limited availability of experimental toxicological data has necessitated development of models that can help us identify the most hazardous and potentially toxic compounds thus prioritizing the experiments on the selected chemicals. Computational tools such as quantitative structure-activity relationship (QSAR) can be used to predict the missing data and classify the chemicals based on their acute predicted responses for existing as well as not yet synthesized chemicals. In the current study, novel, externally validated, highly robust local QSAR models for different chemical classes and moderately robust global QSAR models were developed using partial least squares (PLS) regression technique using a large dataset of 1121 OCs for the fish mortality endpoint. For feature selection, genetic algorithm along with stepwise regression was used while model validation was performed using various stringent validation criteria following the strict rules of OECD guidelines of QSAR validation. The variables included in the models were obtained from simplex representation of molecular structures (SiRMS) (Version 4.1.2.270), Dragon (Version 7.0) and PaDEL-descriptor software (Version 2.20). The final developed models were robust, externally predictive and characterized by a large chemical as well as biological domain. The predictive efficiency of the developed models was then compared with the ECOSAR tool in order to justify the applicability of the developed models in ecotoxicological predictions for organic chemicals. Better predictive efficiency of the developed QSAR models compared to the ECOSAR derived predictions signifies their applicability in early risk assessment of known as well as untested chemicals in order to design safer alternatives to the environment.

Correction to: A new semi-automated workflow for chemical data retrieval and quality checking for modeling applications.

It was highlighted that the original article [1] contained an error in the Funding section. This Correction article states the correct and incorrect versions of the Funding section.

A new semi-automated workflow for chemical data retrieval and quality checking for modeling applications.

The quality of data used for QSAR model derivation is extremely important as it strongly affects the final robustness and predictive power of the model. Ambiguous or wrong structures need to be carefully checked, because they lead to errors in calculation of descriptors, hence leading to meaningless results. The increasing amounts of data, however, have often made it hard to check of very large databases manually. In the light of this, we designed and implemented a semi-automated workflow integrating structural data retrieval from several web-based databases, automated comparison of these data, chemical structure cleaning, selection and standardization of data into a consistent, ready-to-use format that can be employed for modeling. The workflow integrates best practices for data curation that have been suggested in the recent literature. The workflow has been implemented with the freely available KNIME software and is freely available to the cheminformatics community for improvement and application to a broad range of chemical datasets.

Criteria and Application on the Use of Nontesting Methods within a Weight of Evidence Strategy.

Nontesting methods (NTM) proved to be a valuable resource for risk assessment of chemical substances. Indeed, they can be particularly useful when the information provided by different sources was integrated to increase the confidence in the final result. This integration can be sometimes difficult because different methods can lead to conflicting results, and because a clear guideline for integrating information from different sources was not available in the recent past. In this chapter, we present and discuss the recently published guideline from EFSA for integrating and weighting evidence for scientific assessment. Moreover, a practical example on the application of these integration principles on evidence from different in silico models was shown for the assessment of bioconcentration factor (BCF). This example represents a demonstration of the suitability and effectiveness of in silico methods for risk assessment, as well as a practical guide to end-users to perform similar analyses on likely hazardous chemicals.