Nano-QSAR in cell biology: Model of cell viability as a mathematical function of available eclectic data

In Silico Simulation of Impacts of Metal Nano-Oxides on Cell Viability in THP-1 Cells Based on the Correlation Weights of the Fragments of Molecular Structures and Codes of Experimental Conditions Represented by Means of Quasi-SMILES

A simulation of the effect of metal nano-oxides at various concentrations (25, 50, 100, and 200 milligrams per millilitre) on cell viability in THP-1 cells (%) based on data on the molecular structure of the oxide and its concentration is proposed. We used a simplified molecular input-line entry system (SMILES) to represent the molecular structure. So-called quasi-SMILES extends usual SMILES with special codes for experimental conditions (concentration). The approach based on building up models using quasi-SMILES is self-consistent, i.e., the predictive potential of the model group obtained by random splits into training and validation sets is stable. The Monte Carlo method was used as a basis for building up the above groups of models. The CORAL software was applied to building the Monte Carlo calculations. The average determination coefficient for the five different validation sets was R² = 0.806 ± 0.061.

Nano-QSAR modeling for predicting the cytotoxicity of metallic and metal oxide nanoparticles: A review

Given the rapid development of nanotechnology, it is crucial to understand the effects of nanoparticles on living organisms. However, it is laborious to perform toxicological tests on a case-by-case basis. Quantitative structure-activity relationship (QSAR) is an effective computational technique because it saves time, costs, and animal sacrifice. Therefore, this review presents general procedures for the construction and application of nano-QSAR models of metal-based and metal-oxide nanoparticles (MBNPs and MONPs). We also provide an overview of available databases and common algorithms. The molecular descriptors and their roles in the toxicological interpretation of MBNPs and MONPs are systematically reviewed and the future of nano-QSAR is discussed. Finally, we address the growing demand for novel nano-specific descriptors, new computational strategies to address the data shortage, in situ data for regulatory concerns, a better understanding of the physicochemical properties of NPs with bioactivity, and, most importantly, the design of nano-QSAR for real-life environmental predictions rather than laboratory simulations.

Use of quasi-SMILES to build models based on quantitative results from experiments with nanomaterials

Quasi-SMILES deviate from traditional SMILES (simplified molecular input-line entry system) by the extension of additional symbols that encode for conditions of an experiment. Descriptors calculated with SMILES are useful for the development of quantitative structure-property/activity relationships (QSPRs/QSARs), while descriptors calculated with quasi-SMILES can be useful for the development of quantitative models of experimental results obtained under different conditions. Here, this approach has been applied for the development of generalized models using aquatic nanotoxicity data (i.e., related to fish and daphnia). The statistical quality of the above models (pLC50) is quite good with a determination coefficient for the external validation set ranging from 0.62 to 0.71 and RMSE ranging from 0.58 to 0.60. The principle of the approach includes splitting the experimental data into three random distributions defining training, calibration, and validation sets.

Nanomaterials: Quasi-SMILES as a flexible basis for regulation and environmental risk assessment

Basic principles and problems of the systematization of data on nanomaterials are discussed. The eclectic character of nanomaterials is defined as the key difference between nanomaterials and traditional substances. The quasi-SMILES technique is described and discussed. The possible role of the approach is bridging between experimentalists and developers of models for endpoints related to nanomaterials. The use of models on the possible impact of nanomaterials on the environment and human health has been collected and compared. The new criteria of the predictive potential for the above models are discussed. The advantage of the statistical criteria sensitive simultaneously to both the correlation coefficient and the root mean square error noted. The rejection of the border between the effect of the biochemical reality of substances at a molecular level and the effect of experiment conditions at the macro level gives the possibility to develop models that are epistemologically more reliable in the comparison with traditional models based exclusively on the molecular structure-biological activity interdependence (without taking into account experimental conditions). Models of the physicochemical and biochemical behaviour of nanomaterials are necessary in order to develop and apply new industrial achievements, everyday comfort species, medicine, cosmetics, and foods without negative effects on ecology and human health. The CORAL (abbreviation CORrelation And Logic) software provides the user with the possibility to build up nano-QSAR models as a mathematical function of so-called correlation weights of fragments of quasi-SMILES. These models are built up via the Monte Carlo method. Apparently, the quasi-SMILES is a universal representation of nano-reality since there is no limitation to choose the list of eclectic data able to have an impact on nano-phenomena. This paradigm is a convenient language to the conversation of experimentalists and developers of models for nano-phenomena.

CO2 uptake prediction of metal–organic frameworks using quasi-SMILES and Monte Carlo optimization

The first report of quasi-SMILES-based QSPR models for CO2 capture of MOFs based on experimental data.

The QSAR-search of effective agents towards coronaviruses applying the Monte Carlo method

Perhaps there is some similarity between the coronavirus of 2017 and the COVID-19. Consequently, a predictive model for the antiviral activity for the Middle East respiratory syndrome coronavirus (MERS-CoV, 2017) could be useful for designing the strategy and tactics in the struggle with coronaviruses in general and with COVID 19 in particular. Quantitative structure-activity relationships (QSARs) of inhibitory activity to MERS-CoV were developed. The index of ideality of correlation was applied to build up these models for the antiviral activity. The statistical quality of the best model is quite good (r² = 0.84). A mechanistic interpretation of these models based on the molecular features with strong positive (i.e. promoters for endpoint increase) and strong negative (i.e. promoters for endpoint decrease) influence on the inhibitory activity is suggested. A collection of possible biologically active compounds, constructed using data on the above molecular features which are statistically reliable promoters of increase or decrease of the activity, is presented.

Application of quasi-SMILES to the model of gold-nanoparticles uptake in A549 cells

Cell death is critical to human health and is associated with a variety of medical conditions. Therefore, new controllers of cell death are needed for the treatment of diverse diseases. In particular, nanoparticles (NP) are now regularly used in various applications, including a variety of products and medicines. Gold nanoparticles (GNPs) are widely used in the medical field against A549 lung carcinoma cells. The present study is devoted to developing computational models of the cellular uptake potentials by A549 cells of gold nanoparticles (GNPs) under various conditions. Simplified molecular input-line entry system (SMILES) is an efficient tool to represent the molecular structure by a sequence of symbols. Quasi-SMILES represents an extended version of SMILES where symbols to denote physicochemical and/or biochemical conditions are added. In other words, the quasi-SMILES represents a biochemical (medical) phenomenon related to the whole matter (not only molecular structure). We developed models for the cellular utpake potential of gold nanoparticles (GNPs) in A549 [10-11 g Au/Cell] under various conditions based on quasi-SMILES using the Monte Carlo method. The statistical quality of these models is quite good.

Using quasi-SMILES for the predictive modeling of the safety of 574 metal oxide nanoparticles measured in different experimental conditions

The production of nanomaterials continues its rapid growth; however, newly manufactured nanomaterials' environmental and health safety are among the most significant concerns. A safety assessment is usually a lengthy and costly process, so computational studies are often used to complement experimental testing. One of the most time-efficient techniques is structure-activity relationships (SAR) modeling. In this project, we analyzed the Sustainable Nanotechnology (S2NANO) dataset that contains 574 experimental cell viability and toxicity datapoints for Al₂O₃, CuO, Fe₂O₃, Fe₃O₄, SiO₂, TiO₂, and ZnO measured in different conditions. We aimed to develop classification- and regression-based structure-activity relationship models using quasi-SMILES molecular representation. Introduced quasi-SMILES took into consideration all available information, including structural features of nanoparticles (molecular structure, core size, etc.) and related experimental parameters (cell line, dose, exposure time, assay, hydrodynamic size, surface charge, etc.). Resultant regression models demonstrated sufficient predictive power, while classification models demonstrated higher accuracy.

Quasi-SMILES as a basis for the development of models for the toxicity of ZnO nanoparticles

The application of nanomaterials is expanding. Therefore, it is necessary to investigate the relationship between the structure and toxicity of different nanomaterials. Quasi-SMILES is a line of symbols which are codes of corresponding conditions of experiments aimed to estimate the toxicity of ZnO nanoparticles towards the rat via intraperitoneal injections. By means of the Monte Carlo method, the so-called correlation weights for fragments of quasi-SMILES can be calculated. Having the numerical data on the correlation weights one can build up a one-variable model for the toxicity. The checking up of the approach with five random splits of all available data on results of thirty-six experiments into a sub-system of training and sub-system of validation has confirmed the significance of the statistical quality of models obtained with the above approach. The average determination coefficient equal to 0.957 (dispersion 0.010) and average root mean square error equal to 7.25 [mg/kg] (dispersion 0.59 [mg/kg]).

Ecosystem ecology: Models for acute toxicity of pesticides towards Daphnia magna

Quantitative structure - activity relationships (QSARs) which are obtained with a representation of the molecular architecture via simplified molecular input-line entry system (SMILES) are applied to build up predictive models of acute toxicity of pesticides towards Daphnia magna. The acute toxicity towards Daphnia magna is an adequate measure of the ecological impact of various substances. The Monte Carlo technique is the basis to build up the above QSAR models. The statistical quality of suggested models is good: the best model is characterized by n = 103, R² = 0.76, RMSE = 0.91 (training set); n = 53, R² = 0.82, RMSE = 0.87 (validation set). The approach provides the mechanistic interpretation (e.g. aromaticity and branching of carbon skeleton are promoters of increase for toxicity towards Daphnia magna in the case of the examined set of pesticides). The approach is attractive to build up predictive models since instead of a large number of different molecular descriptors the corresponding model is based on solely one optimal descriptor calculated with SMILES and all necessary calculations can be done using the CORAL software available on the Internet (http://ww.insilico.eu/coral).

The cytotoxicity of nanomaterials: Modeling multiple human cells uptake of functionalized magneto-fluorescent nanoparticles via nano-QSAR

The vast majority of nanomaterials have attracted an upsurge of interest since their discovery and considerable researches are being carried out about their adverse outcomes for human health and the environment. In this study, two regression-based quantitative structure-activity relationship models for nanoparticles (nano-QSAR) were established to predict the cellular uptakes of 109 functionalized magneto-fluorescent nanoparticles to pancreatic cancer cells (PaCa2) and human umbilical vein endothelial cells (HUVEC) lines, respectively. The improved SMILES-based optimal descriptors encoded with certain easily available physicochemical properties were proposed to describe the molecular structure characteristics of the involved nanoparticles, and the Monte Carlo method was used for calculating the improved SMILES-based optimal descriptors. Both developed nano-QSAR models for cellular uptake prediction provided satisfactory statistical results, with the squared correlation coefficient (R²) being 0.852 and 0.905 for training sets, and 0.822 and 0.885 for test sets, respectively. Both models were rigorously validated and further extensively compared to literature models. Predominant physicochemical features responsible for cellular uptake were identified by model interpretation. The proposed models could be reasonably expected to provide guidance for synthesizing or choosing safer, more suitable surface modifiers of desired properties prior to their biomedical applications.

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.

QSPR/QSAR: State-of-Art, Weirdness, the Future

Ability of quantitative structure-property/activity relationships (QSPRs/QSARs) to serve for epistemological processes in natural sciences is discussed. Some weirdness of QSPR/QSAR state-of-art is listed. There are some contradictions in the research results in this area. Sometimes, these should be classified as paradoxes or weirdness. These points are often ignored. Here, these are listed and briefly commented. In addition, hypotheses on the future evolution of the QSPR/QSAR theory and practice are suggested. In particular, the possibility of extending of the QSPR/QSAR problematic by searching for the "statistical similarity" of different endpoints is suggested and illustrated by an example for relatively "distanced each from other" endpoints, namely (i) mutagenicity, (ii) anticancer activity, and (iii) blood-brain barrier.

Mathematical modeling of cytotoxicity of metal oxide nanoparticles using the index of ideality correlation criteria

Several types of metal oxide nanoparticles (MO-NPs) are often utilized as one of the novel class of materials in the pharmaceutical industry and human health. The wide use of MO-NPs forces an enhanced understanding of their potential impact on human health and the environment. The research aims to investigate and develop a nano-QFAR (nano-quantitative feature activity relationship) model applying the quasi-SMILES such as cell line, assay, time exposition, concentration, nanoparticles size and metal oxide type for prediction of cell viability (%) of MO-NPs. The total set of 83 quasi-SMILES of MO-NPs divided into training, validation and test sets randomly three times. The statistical model results based on the balance of correlation target function (TF₁) and index of ideality correlation target function (TF₂) and the Monte Carlo optimization were compared. The comparison of two target function results indicated that TF₂ improves the predictability of models. The significance of various eclectic features of both increase and decrease of cell viability (%) is provided. Mechanistic interpretation of significant factors for the model are proposed as well. The sufficient statistical quality of three nano-QFAR models based on TF₂ reveals that the developed models can be efficiency for predictions of the cell viability (%) of MO-NPs.

Toward Rigorous Materials Production: New Approach Methodologies Have Extensive Potential to Improve Current Safety Assessment Practices

Advanced material development, including at the nanoscale, comprises costly and complex challenges coupled to ensuring human and environmental safety. Governmental agencies regulating safety have announced interest toward acceptance of safety data generated under the collective term New Approach Methodologies (NAMs), as such technologies/approaches offer marked potential to progress the integration of safety testing measures during innovation from idea to product launch of nanomaterials. Divided in overall eight main categories, searchable databases for grouping and read across purposes, exposure assessment and modeling, in silico modeling of physicochemical structure and hazard data, in vitro high-throughput and high-content screening assays, dose-response assessments and modeling, analyses of biological processes and toxicity pathways, kinetics and dose extrapolation, consideration of relevant exposure levels and biomarker endpoints typify such useful NAMs. Their application generally agrees with articulated stakeholder needs for improvement of safety testing procedures. They further fit for inclusion and add value in nanomaterials risk assessment tools. Overall 37 of 50 evaluated NAMs and tiered workflows applying NAMs are recommended for considering safer-by-design innovation, including guidance to the selection of specific NAMs in the eight categories. An innovation funnel enriched with safety methods is ultimately proposed under the central aim of promoting rigorous nanomaterials innovation.

In silico design of diacylglycerol acyltransferase-1 (DGAT1) inhibitors based on SMILES descriptors using Monte-Carlo method

Diabetes, obesity and other diseases related to metabolism are worldwide health problems. These syndromes can be well treated when a particular enzyme-based therapy is developed. Diacylglycerol acyltransferase (DGAT; EC 2.3.1.20) is a microsomal enzyme which is responsible for the synthesis of triglycerides from 1,2-diacylglycerol by catalyzing the acyl-CoA-dependent acylation. The obesity and type-II diabetes can be checked by the inhibition of DGAT1 enzyme. Quantitative structure-activity relationship (QSAR) modelling is an essential technique in drug design and development. To study the aspect of DGAT1 inhibitors, Monte-Carlo method-based QSAR was developed for 197 DGAT1 inhibitors. QSAR models were derived by using the optimal descriptor based on SMILES notation. Different statistical parameters including the novel index of ideality of correlation were applied to validate the generated QSAR models. Four random splits were prepared from the data set. The statistical criteria r² = 0.8129, CCC = 0.8979 and Q² = 0.7962 of the validation set of split 1 were the best; therefore, the developed QSAR model of split 1 was decided to be the leading model. The molecular fragments, which were promoter of endpoint increase or decrease were also determined. Thirteen new DGAT1 inhibitors were designed from the lead compound DGAT011.

Importance of Choosing Relevant Biological End Points To Predict Nanoparticle Toxicity with Computational Approaches for Human Health Risk Assessment

Because it is impossible to assess in vitro or in vivo the toxicity of all nanoparticles available on the market on a case-by-case basis, computational approaches have been proposed as useful alternatives to predict in silico the hazard potential of engineered nanoparticles. Despite promising results, a major issue associated with these mathematical models lies in the a priori choice of the physicochemical descriptors and the biological end points. We performed a thorough bibliographic survey on the biological end points used for nanotoxicology purposes and compared them between experimental and computational approaches. They were found to be disparate: while conventional in vitro nanotoxicology assays usually investigate a large array of biological effects using eukaryotic cells (cytotoxicity, pro-inflammatory response, oxidative stress, genotoxicity), computational studies mostly focus on cell viability and also include studies on prokaryotic cells. We may thus wonder the relevance of building complex mathematical models able to predict accurately a biological end point if this latter is not the most relevant to support human health risk assessment. The choice of biological end points clearly deserves to be more carefully discussed. This could bridge the gap between experimental and computational nanotoxicology studies and allow in silico predictive models to reach their full potential.

Quasi-QSAR for predicting the cell viability of human lung and skin cells exposed to different metal oxide nanomaterials

A quasi-QSAR model was developed to predict the cell viability of human lung (BEAS-2B) and skin (HaCaT) cells exposed to 21 types of metal oxide nanomaterials. A wide range of toxicity datasets obtained from the S2NANO (www.s2nano.org) database was used. The data of descriptors representing the physicochemical properties and experimental conditions were coded to quasi-SMILES. In particular, hierarchical cluster analysis (HCA) and min-max normalization method were respectively used in assigning alphanumeric codes for numerical descriptors (e.g., core size, hydrodynamic size, surface charge, and dose) and then quasi-QSAR model performances for both methods were compared. The quasi-QSAR models were developed using CORAL software (www.insilico.eu/coral). Quasi-QSAR model built using quasi-SMILES generated by means of HCA showed better performance than the min-max normalization method. The model showed satisfactory statistical results (R_adj² for the training dataset: 0.71-0.73; R_adj² for the calibration dataset: 0.74-0.82; and R_adj² for the validation dataset: 0.70-0.76).

Towards a generalized toxicity prediction model for oxide nanomaterials using integrated data from different sources

A generalized toxicity classification model for 7 different oxide nanomaterials is presented in this study. A data set extracted from multiple literature sources and screened by physicochemical property based quality scores were used for model development. Moreover, a few more preprocessing techniques, such as synthetic minority over-sampling technique, were applied to address the imbalanced class problem in the data set. Then, classification models using four different algorithms, such as generalized linear model, support vector machine, random forest, and neural network, were developed and their performances were compared to find the best performing preprocessing methods as well as algorithms. The neural network model built using the balanced data set was identified as the model with best predictive performance, while applicability domain was defined using k-nearest neighbours algorithm. The analysis of relative attribute importance for the built neural network model identified dose, formation enthalpy, exposure time, and hydrodynamic size as the four most important attributes. As the presented model can predict the toxicity of the nanomaterials in consideration of various experimental conditions, it has the advantage of having a broader and more general applicability domain than the existing quantitative structure-activity relationship model.

Quasi-SMILES-Based Nano-Quantitative Structure-Activity Relationship Model to Predict the Cytotoxicity of Multiwalled Carbon Nanotubes to Human Lung Cells

Quantitative structure-activity relationship (QSAR) models for nanomaterials (nano-QSAR) were developed to predict the cytotoxicity of 20 different types of multiwalled carbon nanotubes (MWCNTs) to human lung cells by using quasi-SMILES. The optimal descriptors, recorded as quasi-SMILES, were encoded to represent the physicochemical properties and experimental conditions for the MWCNTs from 276 data records collected from previously published studies. The quasi-SMILES used to build the optimal descriptors were (i) diameter, (ii) length, (iii) surface area, (iv) in vitro toxicity assay, (v) cell line, (vi) exposure time, and (vii) dose. The model calculations were performed by using the Monte Carlo method and computed with CORAL software ( www.insilico.eu/coral ). The quasi-SMILES-based nano-QSAR model provided satisfactory statistical results ( R² for internal validation data sets: 0.60-0.80; R²_pred for external validation data sets: 0.81-0.88). The model showed potential for use in the estimation of human lung cell viability after exposure to MWCNTs with the following properties: diameter, 12-74 nm; length, 0.19-20.25 μm; surface area, 11.3-380.0 m²/g; and dose, 0-200 ppm.

Decrypting Strong and Weak Single-Walled Carbon Nanotubes Interactions with Mitochondrial Voltage-Dependent Anion Channels Using Molecular Docking and Perturbation Theory

The current molecular docking study provided the Free Energy of Binding (FEB) for the interaction (nanotoxicity) between VDAC mitochondrial channels of three species (VDAC1-Mus musculus, VDAC1-Homo sapiens, VDAC2-Danio rerio) with SWCNT-H, SWCNT-OH, SWCNT-COOH carbon nanotubes. The general results showed that the FEB values were statistically more negative (p < 0.05) in the following order: (SWCNT-VDAC2-Danio rerio) > (SWCNT-VDAC1-Mus musculus) > (SWCNT-VDAC1-Homo sapiens) > (ATP-VDAC). More negative FEB values for SWCNT-COOH and OH were found in VDAC2-Danio rerio when compared with VDAC1-Mus musculus and VDAC1-Homo sapiens (p < 0.05). In addition, a significant correlation (0.66 > r² > 0.97) was observed between n-Hamada index and VDAC nanotoxicity (or FEB) for the zigzag topologies of SWCNT-COOH and SWCNT-OH. Predictive Nanoparticles-Quantitative-Structure Binding-Relationship models (nano-QSBR) for strong and weak SWCNT-VDAC docking interactions were performed using Perturbation Theory, regression and classification models. Thus, 405 SWCNT-VDAC interactions were predicted using a nano-PT-QSBR classifications model with high accuracy, specificity, and sensitivity (73–98%) in training and validation series, and a maximum AUROC value of 0.978. In addition, the best regression model was obtained with Random Forest (R² of 0.833, RMSE of 0.0844), suggesting an excellent potential to predict SWCNT-VDAC channel nanotoxicity. All study data are available at https://doi.org/10.6084/m9.figshare.4802320.v2.

Nano-QSAR models for predicting cytotoxicity of metal oxide nanoparticles (MONPs) to E. coli

Nanotechnology has been applied to many aspects of human life. Meanwhile, concerns regarding the toxicity of engineered nanomaterials to the environment have also been growing. Herein, an economic and convenient approach based on quantitative structure–activity relationship for nanomaterials (nano-QSAR) was proposed to evaluate the cytotoxicity of metal oxide nanoparticles (MONPs) to E. coli. Six molecular descriptors of 17 MONPs were selected and calculated using Gaussian98 software and DFT-B3LYP method on the LANL2DZ basis set. Two multivariable models, linear and nonlinear, were built based on the calculated molecular descriptors using multiple linear regression (MLR) and support vector machine (SVM) methods, respectively. Results demonstrated that both models presented high reliability, good predictive performance, and fine generalization ability, with all R2 values greater than 0.84. It was also revealed that the lowest unoccupied molecular orbital (LUMO) and molar heat capacity (Cp) were the two key descriptors influencing the cytotoxicity of MONPs.