Green Chemistry in the Synthesis of Pharmaceuticals.

The principles of green chemistry (GC) can be comprehensively implemented in green synthesis of pharmaceuticals by choosing no solvents or green solvents (preferably water), alternative reaction media, and consideration of one-pot synthesis, multicomponent reactions (MCRs), continuous processing, and process intensification approaches for atom economy and final waste reduction. The GC's execution in green synthesis can be performed using a holistic design of the active pharmaceutical ingredient's (API) life cycle, minimizing hazards and pollution, and capitalizing the resource efficiency in the synthesis technique. Thus, the presented review accounts for the comprehensive exploration of GC's principles and metrics, an appropriate implication of those ideas in each step of the reaction schemes, from raw material to an intermediate to the final product's synthesis, and the final execution of the synthesis into scalable industry-based production. For real-life examples, we have discussed the synthesis of a series of established generic pharmaceuticals, starting with the raw materials, and the intermediates of the corresponding pharmaceuticals. Researchers and industries have thoughtfully instigated a green synthesis process to control the atom economy and waste reduction to protect the environment. We have extensively discussed significant reactions relevant for green synthesis, one-pot cascade synthesis, MCRs, continuous processing, and process intensification, which may contribute to the future of green and sustainable synthesis of APIs.

Report of the First ONTOX Stakeholder Network Meeting: Digging Under the Surface of ONTOX Together With the Stakeholders.

The first Stakeholder Network Meeting of the EU Horizon 2020-funded ONTOX project was held on 13-14 March 2023, in Brussels, Belgium. The discussion centred around identifying specific challenges, barriers and drivers in relation to the implementation of non-animal new approach methodologies (NAMs) and probabilistic risk assessment (PRA), in order to help address the issues and rank them according to their associated level of difficulty. ONTOX aims to advance the assessment of chemical risk to humans, without the use of animal testing, by developing non-animal NAMs and PRA in line with 21st century toxicity testing principles. Stakeholder groups (regulatory authorities, companies, academia, non-governmental organisations) were identified and invited to participate in a meeting and a survey, by which their current position in relation to the implementation of NAMs and PRA was ascertained, as well as specific challenges and drivers highlighted. The survey analysis revealed areas of agreement and disagreement among stakeholders on topics such as capacity building, sustainability, regulatory acceptance, validation of adverse outcome pathways, acceptance of artificial intelligence (AI) in risk assessment, and guaranteeing consumer safety. The stakeholder network meeting resulted in the identification of barriers, drivers and specific challenges that need to be addressed. Breakout groups discussed topics such as hazard versus risk assessment, future reliance on AI and machine learning, regulatory requirements for industry and sustainability of the ONTOX Hub platform. The outputs from these discussions provided insights for overcoming barriers and leveraging drivers for implementing NAMs and PRA. It was concluded that there is a continued need for stakeholder engagement, including the organisation of a 'hackathon' to tackle challenges, to ensure the successful implementation of NAMs and PRA in chemical risk assessment.

Does the accounting of the local symmetry fragments in quasi-SMILES improve the predictive potential of the QSAR models of toxicity toward tadpoles?

Models of toxicity to tadpoles have been developed as single parameters based on special descriptors which are sums of correlation weights, molecular features, and experimental conditions. This information is presented by quasi-SMILES. Fragments of local symmetry (FLS) are involved in the development of the model and the use of FLS correlation weights improves their predictive potential. In addition, the index of ideality correlation (IIC) and correlation intensity index (CII) are compared. These two potential predictive criteria were tested in models built through Monte Carlo optimization. The CII was more effective than IIC for the models considered here.

CORAL: Model of Ecological Impact of Heavy Metals on Soils via the Study of Modification of Concentration of Biomolecules in Earthworms (Eisenia fetida).

The traditional application for quantitative structure-property/activity relationships (QSPRs/QSARs) in the fields of thermodynamics, toxicology or drug design is predicting the impact of molecular features using data on the measurable characteristics of substances. However, it is often necessary to evaluate the influence of various exposure conditions and environmental factors, besides the molecular structure. Different enzyme-driven processes lead to the accumulation of metal ions by the worms. Heavy metals are sequestered in these organisms without being released back into the soil. In this study, we propose a novel approach for modeling the absorption of heavy metals, such as mercury and cobalt by worms. The models are based on optimal descriptors calculated for the so-called quasi-SMILES, which incorporate strings of codes reflecting experimental conditions. We modeled the impact on the levels of proteins, hydrocarbons, and lipids in an earthworm's body caused by different combinations of concentrations of heavy metals and exposure time observed over two months of exposure with a measurement interval of 15 days.

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.

A KNIME Workflow to Assist the Analogue Identification for Read-Across, Applied to Aromatase Activity.

The reduction and replacement of in vivo tests have become crucial in terms of resources and animal benefits. The read-across approach reduces the number of substances to be tested, exploiting existing experimental data to predict the properties of untested substances. Currently, several tools have been developed to perform read-across, but other approaches, such as computational workflows, can offer a more flexible and less prescriptive approach. In this paper, we are introducing a workflow to support analogue identification for read-across. The implementation of the workflow was performed using a database of azole chemicals with in vitro toxicity data for human aromatase enzymes. The workflow identified analogues based on three similarities: structural similarity (StrS), metabolic similarity (MtS), and mechanistic similarity (McS). Our results showed how multiple similarity metrics can be combined within a read-across assessment. The use of the similarity based on metabolism and toxicological mechanism improved the predictions in particular for sensitivity. Beyond the results predicting a large population of substances, practical examples illustrate the advantages of the proposed approach.

Using the Correlation Intensity Index to Build a Model of Cardiotoxicity of Piperidine Derivatives.

The assessment of cardiotoxicity is a persistent problem in medicinal chemistry. Quantitative structure-activity relationships (QSAR) are one possible way to build up models for cardiotoxicity. Here, we describe the results obtained with the Monte Carlo technique to develop hybrid optimal descriptors correlated with cardiotoxicity. The predictive potential of the cardiotoxicity models (pIC50, Ki in nM) of piperidine derivatives obtained using this approach provided quite good determination coefficients for the external validation set, in the range of 0.90-0.94. The results were best when applying the so-called correlation intensity index, which improves the predictive potential of a model.

The System of Self-Consistent Models: The Case of Henry's Law Constants.

Data on Henry's law constants make it possible to systematize geochemical conditions affecting atmosphere status and consequently triggering climate changes. The constants of Henry's law are desired for assessing the processes related to atmospheric contaminations caused by pollutants. The most important are those that are capable of long-term movements over long distances. This ability is closely related to the values of Henry's law constants. Chemical changes in gaseous mixtures affect the fate of atmospheric pollutants and ecology, climate, and human health. Since the number of organic compounds present in the atmosphere is extremely large, it is desirable to develop models suitable for predictions for the large pool of organic molecules that may be present in the atmosphere. Here, we report the development of such a model for Henry's law constants predictions of 29,439 compounds using the CORAL software (2023). The statistical quality of the model is characterized by the value of the coefficient of determination for the training and validation sets of about 0.81 (on average).

The system of self-consistent models for pesticide toxicity to Daphnia magna.

Quantitative structure-property/activity relationships (QSPRs/QSARs) are a tool of modern theoretical and computational chemistry. The self-consistent model system is both a method to build up a group of QSPR/QSAR models and an approach to checking the reliability of these models. Here, a group of models of pesticide toxicity toward Daphnia magna for different distributions into training and test sub-sets is compared. This comparison is the basis for formulating the system of self-consistent models. The so-called index of the ideality of correlation (IIC) has been used to improve the above models' predictive potential of pesticide toxicity. The predictive potential of the suggested models should be classified as high since the average value of the determination coefficient for the validation sets is 0.841, and the dispersion is 0.033 (on all five models). The best model (number 4) has an average determination coefficient of 0.89 for the external validation sets (related to all five splits).

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.

Prediction of the Neurotoxic Potential of Chemicals Based on Modelling of Molecular Initiating Events Upstream of the Adverse Outcome Pathways of (Developmental) Neurotoxicity.

Developmental and adult/ageing neurotoxicity is an area needing alternative methods for chemical risk assessment. The formulation of a strategy to screen large numbers of chemicals is highly relevant due to potential exposure to compounds that may have long-term adverse health consequences on the nervous system, leading to neurodegeneration. Adverse Outcome Pathways (AOPs) provide information on relevant molecular initiating events (MIEs) and key events (KEs) that could inform the development of computational alternatives for these complex effects. We propose a screening method integrating multiple Quantitative Structure-Activity Relationship (QSAR) models. The MIEs of existing AOP networks of developmental and adult/ageing neurotoxicity were modelled to predict neurotoxicity. Random Forests were used to model each MIE. Predictions returned by single models were integrated and evaluated for their capability to predict neurotoxicity. Specifically, MIE predictions were used within various types of classifiers and compared with other reference standards (chemical descriptors and structural fingerprints) to benchmark their predictive capability. Overall, classifiers based on MIE predictions returned predictive performances comparable to those based on chemical descriptors and structural fingerprints. The integrated computational approach described here will be beneficial for large-scale screening and prioritisation of chemicals as a function of their potential to cause long-term neurotoxic effects.

Monte Carlo Models for Sub-Chronic Repeated-Dose Toxicity: Systemic and Organ-Specific Toxicity.

The risk-characterization of chemicals requires the determination of repeated-dose toxicity (RDT). This depends on two main outcomes: the no-observed-adverse-effect level (NOAEL) and the lowest-observed-adverse-effect level (LOAEL). These endpoints are fundamental requirements in several regulatory frameworks, such as the Registration, Evaluation, Authorization and Restriction of Chemicals (REACH) and the European Regulation of 1223/2009 on cosmetics. The RDT results for the safety evaluation of chemicals are undeniably important; however, the in vivo tests are time-consuming and very expensive. The in silico models can provide useful input to investigate sub-chronic RDT. Considering the complexity of these endpoints, involving variable experimental designs, this non-testing approach is challenging and attractive. Here, we built eight in silico models for the NOAEL and LOAEL predictions, focusing on systemic and organ-specific toxicity, looking into the effects on the liver, kidney and brain. Starting with the NOAEL and LOAEL data for oral sub-chronic toxicity in rats, retrieved from public databases, we developed and validated eight quantitative structure-activity relationship (QSAR) models based on the optimal descriptors calculated by the Monte Carlo method, using the CORAL software. The results obtained with these models represent a good achievement, to exploit them in a safety assessment, considering the importance of organ-related toxicity.

Virtual Extensive Read-Across: A New Open-Access Software for Chemical Read-Across and Its Application to the Carcinogenicity Assessment of Botanicals.

Read-across applies the principle of similarity to identify the most similar substances to represent a given target substance in data-poor situations. However, differences between the target and the source substances exist. The present study aims to screen and assess the effect of the key components in a molecule which may escape the evaluation for read-across based only on the most similar substance(s) using a new open-access software: Virtual Extensive Read-Across (VERA). VERA provides a means to assess similarity between chemicals using structural alerts specific to the property, pre-defined molecular groups and structural similarity. The software finds the most similar compounds with a certain feature, e.g., structural alerts and molecular groups, and provides clusters of similar substances while comparing these similar substances within different clusters. Carcinogenicity is a complex endpoint with several mechanisms, requiring resource intensive experimental bioassays and a large number of animals; as such, the use of read-across as part of new approach methodologies would support carcinogenicity assessment. To test the VERA software, carcinogenicity was selected as the endpoint of interest for a range of botanicals. VERA correctly labelled 70% of the botanicals, indicating the most similar substances and the main features associated with carcinogenicity.

Integrated In Silico Models for the Prediction of No-Observed-(Adverse)-Effect Levels and Lowest-Observed-(Adverse)-Effect Levels in Rats for Sub-chronic Repeated-Dose Toxicity.

Repeated-dose toxicity (RDT) is a critical endpoint for hazard characterization of chemicals and is assessed to derive safe levels of exposure for human health. Here we present the first attempt to model simultaneously no-observed-(adverse)-effect level (NO(A)EL) and lowest-observed-(adverse)-effect level (LO(A)EL). Classification and regression models were derived based on rat sub-chronic repeated dose toxicity data for 327 compounds from the Fraunhofer RepDose database. Multi-category classification models were built for both NO(A)EL and LO(A)EL though a consensus of statistics- and fragment-based algorithms, while regression models were based on quantitative relationships between the endpoints and SMILES-based attributes. NO(A)EL and LO(A)EL models were integrated, and predictions were compared to exclude inconsistent values. This strategy improved the performance of single models, leading to R2 greater than 0.70, root-mean-square error (RMSE) lower than 0.60 (for regression models), and accuracy of 0.61-0.73 (for classification models) on the validation set, based on the endpoint and the threshold applied for selecting predictions. This study confirms the effectiveness of the modeling strategy presented here for assessing RDT of chemicals using in silico models.

Structures of Endocrine-Disrupting Chemicals Correlate with the Activation of 12 Classic Nuclear Receptors.

Endocrine-disrupting chemicals (EDCs) can inadvertently interact with 12 classic nuclear receptors (NRs) that disrupt the endocrine system and cause adverse effects. There is no widely accepted understanding about what structural features make thousands of EDCs able to activate different NRs as well as how these structural features exert their functions and induce different outcomes at the cellular level. This paper applies the hierarchical characteristic fragment methodology and high-throughput screening molecular docking to comprehensively explore the structural and functional features of EDCs for the 12 NRs based on more than 7000 chemicals from curated datasets. EDCs share three levels of key fragments. The primary and secondary fragments are associated with the binding of EDCs to four groups of receptors: steroidal nuclear receptors (SNRs, including androgen, estrogen, glucocorticoid, mineralocorticoid, and progesterone), retinoic acid receptors, thyroid hormone receptors, and vitamin D receptors. The tertiary fragments determine the activity type by interacting with two key locations in the ligand-binding domains of NRs (N-H5-H3-C and N-H7-H11-C for SNRs and N-H5-H5'-H2'-H3-C and N-H6'-H11-C for non-SNRs). The resulting compiled structural fragments of EDCs together with elucidated compound NR binding modes provide a framework for understanding the interactions between EDCs and NRs, facilitating faster and more accurate screening of EDCs for multiple NRs in the future.

Pesticides, cosmetics, drugs: identical and opposite influences of various molecular features as measures of endpoints similarity and dissimilarity.

The similarity is an important category in natural sciences. A measure of similarity for a group of various biochemical endpoints is suggested. The list of examined endpoints contains (1) toxicity of pesticides towards rainbow trout; (2) human skin sensitization; (3) mutagenicity; (4) toxicity of psychotropic drugs; and (5) anti HIV activity. Further applying and evolution of the suggested approach is discussed. In particular, the conception of the similarity (dissimilarity) of endpoints can play the role of a "useful bridge" between quantitative structure property/activity relationships (QSPRs/QSARs) and read-across technique.

Defining the Human-Biota Thresholds of Toxicological Concern for Organic Chemicals in Freshwater: The Proposed Strategy of the LIFE VERMEER Project Using VEGA Tools.

Several tons of chemicals are released every year into the environment and it is essential to assess the risk of adverse effects on human health and ecosystems. Risk assessment is expensive and time-consuming and only partial information is available for many compounds. A consolidated approach to overcome this limitation is the Threshold of Toxicological Concern (TTC) for assessment of the potential health impact and, more recently, eco-TTCs for the ecological aspect. The aim is to allow a safe assessment of substances with poor toxicological characterization. Only limited attempts have been made to integrate the human and ecological risk assessment procedures in a "One Health" perspective. We are proposing a strategy to define the Human-Biota TTCs (HB-TTCs) as concentrations of organic chemicals in freshwater preserving both humans and ecological receptors at the same time. Two sets of thresholds were derived: general HB-TTCs as preliminary screening levels for compounds with no eco- and toxicological information, and compound-specific HB-TTCs for chemicals with known hazard assessment, in terms of Predicted No effect Concentration (PNEC) values for freshwater ecosystems and acceptable doses for human health. The proposed strategy is based on freely available public data and tools to characterize and group chemicals according to their toxicological profiles. Five generic HB-TTCs were defined, based on the ecotoxicological profiles reflected by the Verhaar classes, and compound-specific thresholds for more than 400 organic chemicals with complete eco- and toxicological profiles. To complete the strategy, the use of in silico models is proposed to predict the required toxicological properties and suitable models already available on the VEGAHUB platform are listed.

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.

Structures of Endocrine-Disrupting Chemicals Determine Binding to and Activation of the Estrogen Receptor α and Androgen Receptor.

Endocrine-disrupting chemicals (EDCs) can interact with nuclear receptors, including estrogen receptor α (ERα) and androgen receptor (AR), to affect the normal endocrine system function, causing severe symptoms. Limited studies queried the EDC mechanisms, focusing on limited chemicals or a set of structurally similar compounds. It remained uncertain how hundreds of diverse EDCs could bind to ERα and AR and cause distinct functional consequences. Here, we employed a series of computational methodologies to investigate the structural features of EDCs that bind to and activate ERα and AR based on more than 4000 compounds. We used molecular docking and molecular dynamics simulations to elucidate the functional consequences and validated structure-function correlations experimentally using a time-resolved fluorescence resonance energy-transfer assay. We found that EDCs share three levels of key fragments. Primary (20 for ERα and 18 for AR) and secondary fragments (38 for ERα and 29 for AR) are responsible for the binding to receptors, and tertiary fragments determine the activity type (agonist, antagonist, or mixed). In summary, our study provides a general mechanism for the EDC function. Discovering the three levels of key fragments may drive fast screening and evaluation of potential EDCs from large sets of commercially used synthetic compounds.

Value and limitation of structure-based profilers to characterize developmental and reproductive toxicity potential.

The uncertainty regarding the safety of chemicals leaching from food packaging triggers attention. In silico models provide solutions for screening of these chemicals, since many are toxicologically uncharacterized. For hazard assessment, information on developmental and reproductive toxicity (DART) is needed. The possibility to apply in silico toxicology to identify and quantify DART alerts was investigated. Open-source models and profilers were applied to 195 packaging chemicals and analogues. An approach based on DART and estrogen receptor (ER) binding profilers and molecular docking was able to identify all except for one chemical with documented DART properties. Twenty percent of the chemicals in the database known to be negative in experimental studies were classified as positive. The scheme was then applied to 121 untested chemicals. Alerts were identified for sixteen of them, five being packaging substances, the others structural analogues. Read-across was then developed to translate alerts into quantitative toxicological values. They can be used to calculate margins of exposure (MoE), the size of which reflects safety concern. The application of this approach appears valuable for hazard characterization of toxicologically untested packaging migrants. It is an alternative to the use of default uncertainty factor (UF) applied to animal chronic toxicity value to handle absence of DART data in hazard characterization.