Toward Nano-Specific In Silico NAMs: How to Adjust Nano-QSAR to the Recent Advancements of Nanotoxicology?

The rapid development of engineered nanomaterials (ENMs) causes humans to become increasingly exposed to them. Therefore, a better understanding of the health impact of ENMs is highly demanded. Considering the 3Rs (Replacement, Reduction, and Refinement) principle, in vitro and computational methods are excellent alternatives for testing on animals. Among computational methods, nano-quantitative structure-activity relationship (nano-QSAR), which links the physicochemical and structural properties of EMNs with biological activities, is one of the leading method. The nature of toxicological experiments has evolved over the last decades; currently, one experiment can provide thousands of measurements of the organism's functioning at the molecular level. At the same time, the capacity of the in vitro systems to mimic the human organism is also improving significantly. Hence, the authors would like to discuss whether the nano-QSAR approach follows modern toxicological studies and takes full advantage of the opportunities offered by modern toxicological platforms. Challenges and possibilities for improving data integration are underlined narratively, including the need for a consensus built between the in vitro and the QSAR domains.

Hybrid Machine Learning and Experimental Studies of Antiviral Potential of Ionic Liquids against P100, MS2, and Phi6.

Viruses are a group of widespread organisms that are often responsible for very dangerous diseases, as most of them follow a mechanism to multiply and infect their hosts as quickly as possible. Pathogen viruses also mutate regularly, with the result that measures to prevent virus transmission and recover from the disease caused are often limited. The development of new substances is very time-consuming and highly budgeted and requires the sacrifice of many living organisms. Computational chemistry methods allow faster analysis at a much lower cost and, most importantly, reduce the number of living organisms sacrificed experimentally to a minimum. Ionic liquids (ILs) are a group of chemical compounds that could potentially find a wide range of applications due to their potential virucidal activity. In our study, we conducted a complex computational analysis to predict the antiviral activity of ionic liquids against three surrogate viruses: two nonenveloped viruses, Listeria monocytogenes phage P100 and Escherichia coli phage MS2, and one enveloped virus, Pseudomonas syringae phage Phi6. Based on experimental data of toxic activity (logEC90), we assigned activity classes to 154 ILs. Prediction models were created and validated according to the Organization for Economic Co-operation and Development (OECD) recommendations using the Classification Tree method. Further, we performed an external validation of our models through virtual screening on a set of 1277 theoretically generated ionic liquids and then selected 10 active ionic liquids, which were synthesized to verify their activity against the analyzed viruses. Our study proved the effectiveness and efficiency of computational methods to predict the antiviral activity of ionic liquids. Thus, computational models are a cost-effective alternative approach compared with time-consuming experimental studies where live animals are involved.

How Does the Study MD of pH-Dependent Exposure of Nanoparticles Affect Cellular Uptake of Anticancer Drugs?

The lack of knowledge about the uptake of NPs by biological cells poses a significant problem for drug delivery. For this reason, designing an appropriate model is the main challenge for modelers. To address this problem, molecular modeling studies that can describe the mechanism of cellular uptake of drug-loaded nanoparticles have been conducted in recent decades. In this context, we developed three different models for the amphipathic nature of drug-loaded nanoparticles (MTX-SS-γ-PGA), whose cellular uptake mechanism was predicted by molecular dynamics studies. Many factors affect nanoparticle uptake, including nanoparticle physicochemical properties, protein-particle interactions, and subsequent agglomeration, diffusion, and sedimentation. Therefore, the scientific community needs to understand how these factors can be controlled and the NP uptake of nanoparticles. Based on these considerations, in this study, we investigated for the first time the effects of the selected physicochemical properties of the anticancer drug methotrexate (MTX) grafted with hydrophilic-γ-polyglutamic acid (MTX-SS-γ-PGA) on its cellular uptake at different pH values. To answer this question, we developed three theoretical models describing drug-loaded nanoparticles (MTX-SS-γ-PGA) at three different pH values, such as (1) pH 7.0 (the so-called neutral pH model), (2) pH 6.4 (the so-called tumor pH model), and (3) pH 2.0 (the so-called stomach pH model). Exceptionally, the electron density profile shows that the tumor model interacts more strongly with the head groups of the lipid bilayer than the other models due to charge fluctuations. Hydrogen bonding and RDF analyses provide information about the solution of the NPs with water and their interaction with the lipid bilayer. Finally, dipole moment and HOMO-LUMO analysis showed the free energy of the solution in the water phase and chemical reactivity, which are particularly useful for determining the cellular uptake of the NPs. The proposed study provides fundamental insights into molecular dynamics (MD) that will allow researchers to determine the influence of pH, structure, charge, and energetics of NPs on the cellular uptake of anticancer drugs. We believe that our current study will be useful in developing a new model for drug delivery to cancer cells with a much more efficient and less time-consuming model.

How the Structure of Per- and Polyfluoroalkyl Substances (PFAS) Influences Their Binding Potency to the Peroxisome Proliferator-Activated and Thyroid Hormone Receptors-An In Silico Screening Study.

In this study, we investigated PFAS (per- and polyfluoroalkyl substances) binding potencies to nuclear hormone receptors (NHRs): peroxisome proliferator-activated receptors (PPARs) α, ß, and γ and thyroid hormone receptors (TRs) α and ß. We have simulated the docking scores of 43 perfluoroalkyl compounds and based on these data developed QSAR (Quantitative Structure-Activity Relationship) models for predicting the binding probability to five receptors. In the next step, we implemented the developed QSAR models for the screening approach of a large group of compounds (4464) from the NORMAN Database. The in silico analyses indicated that the probability of PFAS binding to the receptors depends on the chain length, the number of fluorine atoms, and the number of branches in the molecule. According to the findings, the considered PFAS group bind to the PPARα, ß, and γ only with low or moderate probability, while in the case of TR α and ß it is similar except that those chemicals with longer chains show a moderately high probability of binding.

Transcriptomics-Based and AOP-Informed Structure-Activity Relationships to Predict Pulmonary Pathology Induced by Multiwalled Carbon Nanotubes.

This study presents a novel strategy that employs quantitative structure-activity relationship models for nanomaterials (Nano-QSAR) for predicting transcriptomic pathway level response using lung tissue inflammation, an essential key event (KEs) in the existing adverse outcome pathway (AOP) for lung fibrosis, as a model response. Transcriptomic profiles of mouse lungs exposed to ten different multiwalled carbon nanotubes (MWCNTs) are analyzed using statistical and bioinformatics tools. Three pathways "agranulocyte adhesion and diapedesis," "granulocyte adhesion and diapedesis," and "acute phase signaling," that (1) are commonly perturbed across the MWCNTs panel, (2) show dose response (Benchmark dose, BMDs), and (3) are anchored to the KEs identified in the lung fibrosis AOP, are considered in modelling. The three pathways are associated with tissue inflammation. The results show that the aspect ratio (κ) of MWCNTs is directly correlated with the pathway BMDs. The study establishes a methodology for QSAR construction based on canonical pathways and proposes a MWCNTs grouping strategy based on the κ-values of the specific pathway associated genes. Finally, the study shows how the AOP framework can help guide QSAR modelling efforts; conversely, the outcome of the QSAR modelling can aid in refining certain aspects of the AOP in question (here, lung fibrosis).

Importance of Surface Topography in Both Biological Activity and Catalysis of Nanomaterials: Can Catalysis by Design Guide Safe by Design?

It is acknowledged that the physicochemical properties of nanomaterials (NMs) have an impact on their toxicity and, eventually, their pathogenicity. These properties may include the NMs' surface chemical composition, size, shape, surface charge, surface area, and surface coating with ligands (which can carry different functional groups as well as proteins). Nanotopography, defined as the specific surface features at the nanoscopic scale, is not widely acknowledged as an important physicochemical property. It is known that the size and shape of NMs determine their nanotopography which, in turn, determines their surface area and their active sites. Nanotopography may also influence the extent of dissolution of NMs and their ability to adsorb atoms and molecules such as proteins. Consequently, the surface atoms (due to their nanotopography) can influence the orientation of proteins as well as their denaturation. However, although it is of great importance, the role of surface topography (nanotopography) in nanotoxicity is not much considered. Many of the issues that relate to nanotopography have much in common with the fundamental principles underlying classic catalysis. Although these were developed over many decades, there have been recent important and remarkable improvements in the development and study of catalysts. These have been brought about by new techniques that have allowed for study at the nanoscopic scale. Furthermore, the issue of quantum confinement by nanosized particles is now seen as an important issue in studying nanoparticles (NPs). In catalysis, the manipulation of a surface to create active surface sites that enhance interactions with external molecules and atoms has much in common with the interaction of NP surfaces with proteins, viruses, and bacteria with the same active surface sites of NMs. By reviewing the role that surface nanotopography plays in defining many of the NMs' surface properties, it reveals the need for its consideration as an important physicochemical property in descriptive and predictive toxicology. Through the manipulation of surface topography, and by using principles developed in catalysis, it may also be possible to make safe-by-design NMs with a reduction of the surface properties which contribute to their toxicity.

Risk Governance of Emerging Technologies Demonstrated in Terms of its Applicability to Nanomaterials.

Nanotechnologies have reached maturity and market penetration that require nano-specific changes in legislation and harmonization among legislation domains, such as the amendments to REACH for nanomaterials (NMs) which came into force in 2020. Thus, an assessment of the components and regulatory boundaries of NMs risk governance is timely, alongside related methods and tools, as part of the global efforts to optimise nanosafety and integrate it into product design processes, via Safe(r)-by-Design (SbD) concepts. This paper provides an overview of the state-of-the-art regarding risk governance of NMs and lays out the theoretical basis for the development and implementation of an effective, trustworthy and transparent risk governance framework for NMs. The proposed framework enables continuous integration of the evolving state of the science, leverages best practice from contiguous disciplines and facilitates responsive re-thinking of nanosafety governance to meet future needs. To achieve and operationalise such framework, a science-based Risk Governance Council (RGC) for NMs is being developed. The framework will provide a toolkit for independent NMs' risk governance and integrates needs and views of stakeholders. An extension of this framework to relevant advanced materials and emerging technologies is also envisaged, in view of future foundations of risk research in Europe and globally.

Multi-Objective Genetic Algorithm (MOGA) As a Feature Selecting Strategy in the Development of Ionic Liquids' Quantitative Toxicity-Toxicity Relationship Models.

Quantitative toxicity-toxicity relationship (QTTR) models have a great potential for improving the meaning of toxicological tests conducted on simple organisms. These models allow predicting the toxicological effect of a chemical based on its known toxicological effect in different toxicity tests, even against a different organism. This fact poses a great potential for predicting the toxicity of chemicals against higher organisms based on the results against lower ones. However, the possibility of developing such models is often restricted due to the low availability of data. We present a case study of developing the QTTR model for ionic liquids in different toxicological tests against the same species, in the face of insufficient experimental data (an additional confirmation for a different species is provided in the Supporting Information). In the presented case, we use a series of quantitative structure-activity relationship (QSAR) models developed to deliver the data concerning the toxicity of ionic liquids against human HeLa and MCF-7 cancer cell lines. We use these data to develop a QTTR model with an R2 as high as 0.8. The benefit of applying the multi-objective genetic algorithm (MOGA-a genetic algorithm allowing for selection of the best set of explanatory features for several different dependent variables at the same time) as a QSAR model feature selecting strategy is presented and discussed.

The Essential Elements of a Risk Governance Framework for Current and Future Nanotechnologies.

Societies worldwide are investing considerable resources into the safe development and use of nanomaterials. Although each of these protective efforts is crucial for governing the risks of nanomaterials, they are insufficient in isolation. What is missing is a more integrative governance approach that goes beyond legislation. Development of this approach must be evidence based and involve key stakeholders to ensure acceptance by end users. The challenge is to develop a framework that coordinates the variety of actors involved in nanotechnology and civil society to facilitate consideration of the complex issues that occur in this rapidly evolving research and development area. Here, we propose three sets of essential elements required to generate an effective risk governance framework for nanomaterials. (1) Advanced tools to facilitate risk-based decision making, including an assessment of the needs of users regarding risk assessment, mitigation, and transfer. (2) An integrated model of predicted human behavior and decision making concerning nanomaterial risks. (3) Legal and other (nano-specific and general) regulatory requirements to ensure compliance and to stimulate proactive approaches to safety. The implementation of such an approach should facilitate and motivate good practice for the various stakeholders to allow the safe and sustainable future development of nanotechnology.

Scanning electron microscopy image representativeness: morphological data on nanoparticles.

A sample of a nanomaterial contains a distribution of nanoparticles of various shapes and/or sizes. A scanning electron microscopy image of such a sample often captures only a fragment of the morphological variety present in the sample. In order to quantitatively analyse the sample using scanning electron microscope digital images, and, in particular, to derive numerical representations of the sample morphology, image content has to be assessed. In this work, we present a framework for extracting morphological information contained in scanning electron microscopy images using computer vision algorithms, and for converting them into numerical particle descriptors. We explore the concept of image representativeness and provide a set of protocols for selecting optimal scanning electron microscopy images as well as determining the smallest representative image set for each of the morphological features. We demonstrate the practical aspects of our methodology by investigating tricalcium phosphate, Ca3 (PO4 )2 , and calcium hydroxyphosphate, Ca5 (PO4 )3 (OH), both naturally occurring minerals with a wide range of biomedical applications.

Compilation of Data and Modelling of Nanoparticle Interactions and Toxicity in the NanoPUZZLES Project.

The particular properties of nanomaterials have led to their rapidly increasing use in diverse fields of application. However, safety assessment is not keeping pace and there are still gaps in the understanding of their hazards. Computational models predicting nanotoxicity, such as (quantitative) structure-activity relationships ((Q)SARs), can contribute to safety evaluation, in line with general efforts to apply alternative methods in chemical risk assessment. Their development is highly dependent on the availability of reliable and high quality experimental data, both regarding the compounds' properties as well as the measured toxic effects. In particular, "nano-QSARs" should take the nano-specific characteristics into account. The information compiled needs to be well organized, quality controlled and standardized. Integrating the data in an overarching, structured data collection aims to (a) organize the data in a way to support modelling, (b) make (meta)data necessary for modelling available, and (c) add value by making a comparison between data from different sources possible.Based on the available data, specific descriptors can be derived to parameterize the nanomaterial-specific structure and physico-chemical properties appropriately. Furthermore, the interactions between nanoparticles and biological systems as well as small molecules, which can lead to modifications of the structure of the active nanoparticles, need to be described and taken into account in the development of models to predict the biological activity and toxicity of nanoparticles. The EU NanoPUZZLES project was part of a global cooperative effort to advance data availability and modelling approaches supporting the characterization and evaluation of nanomaterials.

Metabolomic Biomarkers in Urine of Cushing's Syndrome Patients.

Cushing's syndrome (CS) is a disease which results from excessive levels of cortisol in the human body. The disorder is associated with various signs and symptoms which are also common for the general population not suffering from compound hypersecretion. Thus, more sensitive and selective methods are required for the diagnosis of CS. This follow-up study was conducted to determine which steroid metabolites could serve as potential indicators of CS and possible subclinical hypercortisolism in patients diagnosed with so called non-functioning adrenal incidentalomas (AIs). Urine samples from negative controls (n = 37), patients with CS characterized by hypercortisolism and excluding iatrogenic CS (n = 16), and patients with non-functioning AIs with possible subclinical Cushing's syndrome (n = 25) were analyzed using gas chromatography-mass spectrometry (GC/MS) and gas chromatograph equipped with flame ionization detector (GC/FID). Statistical and multivariate methods were applied to investigate the profile differences between examined individuals. The analyses revealed hormonal differences between patients with CS and the rest of examined individuals. The concentrations of selected metabolites of cortisol, androgens, and pregnenetriol were elevated whereas the levels of tetrahydrocortisone were decreased for CS when opposed to the rest of the study population. Moreover, after analysis of potential confounding factors, it was also possible to distinguish six steroid hormones which discriminated CS patients from other study subjects. The obtained discriminant functions enabled classification of CS patients and AI group characterized by mild hypersecretion of cortisol metabolites. It can be concluded that steroid hormones selected by applying urinary profiling may serve the role of potential biomarkers of CS and can aid in its early diagnosis.

Towards the Application of Structure-Property Relationship Modeling in Materials Science: Predicting the Seebeck Coefficient for Ionic Liquid/Redox Couple Systems.

This work focuses on determining the influence of both ionic-liquid (IL) type and redox couple concentration on Seebeck coefficient values of such a system. The quantitative structure-property relationship (QSPR) and read-across techniques are proposed as methods to identify structural features of ILs (mixed with LiI/I2 redox couple), which have the most influence on the Seebeck coefficient (Se ) values of the system. ILs consisting of small, symmetric cations and anions with high values of vertical electron binding energy are recognized as those with the highest values of Se . In addition, the QSPR model enables the values of Se to be predicted for each IL that belongs to the applicability domain of the model. The influence of the redox-couple concentration on values of Se is also quantitatively described. Thus, it is possible to calculate how the value of Se will change with changing redox-couple concentration. The presence of the LiI/I2 redox couple in lower concentrations increases the values of Se , as expected.

Geometry optimization method versus predictive ability in QSPR modeling for ionic liquids.

Computational techniques, such as Quantitative Structure-Property Relationship (QSPR) modeling, are very useful in predicting physicochemical properties of various chemicals. Building QSPR models requires calculating molecular descriptors and the proper choice of the geometry optimization method, which will be dedicated to specific structure of tested compounds. Herein, we examine the influence of the ionic liquids' (ILs) geometry optimization methods on the predictive ability of QSPR models by comparing three models. The models were developed based on the same experimental data on density collected for 66 ionic liquids, but with employing molecular descriptors calculated from molecular geometries optimized at three different levels of the theory, namely: (1) semi-empirical (PM7), (2) ab initio (HF/6-311+G*) and (3) density functional theory (B3LYP/6-311+G*). The model in which the descriptors were calculated by using ab initio HF/6-311+G* method indicated the best predictivity capabilities ([Formula: see text] = 0.87). However, PM7-based model has comparable values of quality parameters ([Formula: see text] = 0.84). Obtained results indicate that semi-empirical methods (faster and less expensive regarding CPU time) can be successfully employed to geometry optimization in QSPR studies for ionic liquids.

Development of a novel in silico model of zeta potential for metal oxide nanoparticles: a nano-QSPR approach.

Once released into the aquatic environment, nanoparticles (NPs) are expected to interact (e.g. dissolve, agglomerate/aggregate, settle), with important consequences for NP fate and toxicity. A clear understanding of how internal and environmental factors influence the NP toxicity and fate in the environment is still in its infancy. In this study, a quantitative structure-property relationship (QSPR) approach was employed to systematically explore factors that affect surface charge (zeta potential) under environmentally realistic conditions. The nano-QSPR model developed with multiple linear regression (MLR) was characterized by high robustness [Formula: see text] and external predictivity [Formula: see text] The results clearly showed that zeta potential values varied markedly as functions of the ionic radius of the metal atom in the metal oxides, confirming that agglomeration and the extent of release of free MexOy largely depend on their intrinsic properties. A developed nano-QSPR model was successfully applied to predict zeta potential in an ionized solution of NPs for which experimentally determined values of response have been unavailable. Hence, the application of our model is possible when the values of zeta potential in the ionized solution for metal oxide nanoparticles are undetermined, without the necessity of performing more time consuming and expensive experiments. We believe that our studies will be helpful in predicting the conditions under which MexOy is likely to become problematic for the environment and human health.

Extrapolating between toxicity endpoints of metal oxide nanoparticles: Predicting toxicity to Escherichia coli and human keratinocyte cell line (HaCaT) with Nano-QTTR.

Synthesis of novel nanoparticles should always be accompanied by a comprehensive assessment of risk to human health and to ecosystem. Application of in silico models is encouraged by regulatory authorities to fill the data gaps related to the properties of nanoparticles affecting the environment and human health. Interspecies toxicity correlations provide a tool for estimation of contaminant's sensitivity with known levels of uncertainty for a diverse pool of species. We propose here first interspecies cytotoxicity correlation models between Escherichia coli (prokaryotic system) and human keratinocyte cell line (HaCaT) (eukaryotic system) to assess the discriminatory features for cytotoxicity of metal oxide nanoparticles. The nano-QTTR models can be employed for extrapolating cytotoxicity to E. coli and human keratinocyte cell line (HaCaT) for metal nanoparticles when the data for the other species are available. Informative illustrations of the contributing mechanisms of toxic action of the metal oxide nanoparticles to the HaCaT cell line as well as to the E. coli are identified from the developed nano quantitative toxicity-toxicity relationship (nano-QTTR) models.

Comparing the CORAL and Random Forest approaches for modelling the in vitro cytotoxicity of silica nanomaterials.

Nanotechnology is one of the most important technological developments of the 21st century. In silico methods to predict toxicity, such as quantitative structure-activity relationships (QSARs), promote the safe-by-design approach for the development of new materials, including nanomaterials. In this study, a set of cytotoxicity experimental data corresponding to 19 data points for silica nanomaterials were investigated, to compare the widely employed CORAL and Random Forest approaches in terms of their usefulness for developing so-called 'nano-QSAR' models. 'External' leave-one-out cross-validation (LOO) analysis was performed, to validate the two different approaches. An analysis of variable importance measures and signed feature contributions for both algorithms was undertaken, in order to interpret the models developed. CORAL showed a more pronounced difference between the average coefficient of determination (R²) for training and for LOO (0.83 and 0.65 for training and LOO, respectively), compared to Random Forest (0.87 and 0.78 without bootstrap sampling, 0.90 and 0.78 with bootstrap sampling), which may be due to overfitting. With regard to the physicochemical properties of the nanomaterials, the aspect ratio and zeta potential were found to be the two most important variables for Random Forest, and the average feature contributions calculated for the corresponding descriptors were consistent with the clear trends observed in the data set: less negative zeta potential values and lower aspect ratio values were associated with higher cytotoxicity. In contrast, CORAL failed to capture these trends.

The performance of selected semi-empirical and DFT methods in studying C60 fullerene derivatives.

The capability of reproducing the open circuit voltages (V(oc)) of 15 representative C60 fullerene derivatives was tested using the selected quantum mechanical methods (B3LYP, PM6, and PM7) together with the two one-electron basis sets. Certain theoretical treatments (e.g. PM6) were found to be satisfactory for preliminary estimates of the open circuit voltages (V(oc)), whereas the use of the B3LYP/6-31G(d) approach has been proven to assure highly accurate results. We also examined the structural similarity of 19 fullerene derivatives by employing principle component analysis (PCA). In order to express the structural features of the studied compounds we used molecular descriptors calculated with semi-empirical (PM6 and PM7) and density functional (B3LYP/6-31G(d)) methods separately. In performing PCA, we noticed that semi-empirical methods (i.e. PM6 and PM7) seem satisfactory for molecules, in which one can distinguish the aromatic and the aliphatic parts in the cyclopropane ring of PCBM (phenyl-C61-buteric acid methyl ester) and they significantly overestimate the energy of the highest occupied molecular orbital (E(HOMO)). The use of the B3LYP functional, however, is recommended for studying methanofullerenes, which closely resemble the structure of PCBM, and for their modifications.

Novel approach for efficient predictions properties of large pool of nanomaterials based on limited set of species: nano-read-across.

Creating suitable chemical categories and developing read-across methods, supported by quantum mechanical calculations, can be an effective solution to solving key problems related to current scarcity of data on the toxicity of various nanoparticles. This study has demonstrated that by applying a nano-read-across, the cytotoxicity of nano-sized metal oxides could be estimated with a similar level of accuracy as provided by quantitative structure-activity relationship for nanomaterials (nano-QSAR model(s)). The method presented is a suitable computational tool for the preliminary hazard assessment of nanomaterials. It also could be used for the identification of nanomaterials that may pose potential negative impact to human health and the environment. Such approaches are especially necessary when there is paucity of relevant and reliable data points to develop and validate nano-QSAR models.

A new metric for long-range transport potential of chemicals.

We propose a new metric for long-range transport potential (LRTP), GIF, based on source-receptor analyses and evaluate the LRTP and persistence of a wide variety of chlorinated and brominated organic compounds using GIF and overall persistence (POV), respectively. We calculated GIF and POV using our global 3D dynamic multimedia model (FATE). Physicochemical properties were obtained from quantitative structure-property relationship (QSPR) models. The FATE-QSPR combined model enabled us to systematically investigate the LRTP and persistence of a wide variety of chemical substances. On average, the estimated GIF and POV for chlorinated compounds were larger than those for their brominated counterparts, with the largest and smallest values found for polychlorinated biphenyls and polybrominated dibenzodioxins, respectively. We also compared GIF with four differently defined LRTP metrics and two LRTP metrics obtained from a simple model. The results of our analyses indicate that the LRTP ranks can differ considerably among LRTP metrics, the differences being dependent on the governing environmental processes, relevant physicochemical properties, and multimedia model.